第六章空间数据的分析

1、第六章第六章 空间数据的分析空间数据的分析6.1 6.1 空间数据分析的内容空间数据分析的内容 6.2 6.2 空间数据的基本运算空间数据的基本运算6. 3 空间数据的查询空间数据的查询6.4 6.4 空间数据的叠置分析空间数据的叠置分析6.5 6.5 空间数据的邻域分析空间数据的邻域分析6.6 6.6 空间数据的网络分析空间数据的网络分析6.7 6.7 空间数据的统计分析空间数据的统计分析6.16.1空间分析的内容空间分析的内容一、空间分析概述一、空间分析概述空间分析通过空间数据的分析算法，获取地理对象的空间位置、空间分析通过空间数据的分析算法，获取地理对象的空间位置、空间分布、空间形态、空

2、间演变等信息。空间分布、空间形态、空间演变等信息。GISGIS中空间分析的二个层次中空间分析的二个层次1 1、咨询式空间分析、咨询式空间分析实现空间数据的查询检索及集合分析。实现空间数据的查询检索及集合分析。2 2、产生式空间分析、产生式空间分析基于基于GISGIS中拓朴关系和空间操作运算、空间统计分析，及将中拓朴关系和空间操作运算、空间统计分析，及将GISGIS作作为通用工具同其它专业模型结合，实现空间数据的模拟和分析。为通用工具同其它专业模型结合，实现空间数据的模拟和分析。二、空间分析的主要内容二、空间分析的主要内容查询检索查询检索 拓扑查询拓扑查询 位置查询位置查询 属性查询属性查询 区

3、域查询；区域查询；形态分析形态分析 面积量算面积量算 距离量算距离量算 质心计算质心计算 周长量算；周长量算；地形分析地形分析 等值线分析等值线分析 坡度坡向分析坡度坡向分析 分水岭分析分水岭分析 视域分视域分 析析 剖面分析剖面分析 ； 叠置分析叠置分析 视觉信息复合视觉信息复合 条件叠置条件叠置 无条件叠置；无条件叠置；邻域分析邻域分析 缓冲器分析缓冲器分析 泰森多边形分析泰森多边形分析 插值拟合分析；插值拟合分析；网络分析网络分析 最短或最佳路经分析最短或最佳路经分析 空间规划；空间规划；图象分析图象分析 图象增强图象增强 图象分割图象分割 图象细化；图象细化；应用模型分析应用模型分析

4、数学模型数学模型 统计模型统计模型 逻辑模型。逻辑模型。、 6. 2 空间数据的基本运算空间数据的基本运算 1 、几何计算、几何计算1）长度的计算）长度的计算（1） 矢量数据的长度的计算矢量数据的长度的计算 两点（两点（x1,y1) (x2,y2)之间距离的计算之间距离的计算 d= (x2- x1)2 +( y2 -y1) 2 点（点（x0,y0)到直线到直线 Ax+By+C=0 之间欧氏距离的计算之间欧氏距离的计算 d= |Ax0+By0+C| A2 + B2 线目标的长度线目标的长度 L=（2）栅格数据的长度的计算）栅格数据的长度的计算用用8邻域方向累加地物骨架线通过的栅格数目邻域方向累加

5、地物骨架线通过的栅格数目,对角线方向对角线方向* 2 niid1 2)面积的计算面积的计算空间数据的自动量算是空间数据的自动量算是GIS 的重要功能，也是进行空间分析的的重要功能，也是进行空间分析的定量化基础。定量化基础。 面积的量算（辛普森（面积的量算（辛普森（Simposion)公式）公式）（1） 矢量格式下面积的量算通常用辛普森公式：矢量格式下面积的量算通常用辛普森公式： Y Si =1/2 ( yi+1 + yi ) ( xi - xi+1 ) S = Si + ( yn + y 1) ( x1 - xn ) /2 S =1/2 yi (xi+1 - xi -1 ) （2） 栅格数据面

6、积的计算栅格数据面积的计算 统计相同属性的格网数目统计相同属性的格网数目i=1n-1n123n-1i=1nxn +1 = x1x0 = xnX3）辛普森公式的用处）辛普森公式的用处（1）编程时自动计算面积的公式编程时自动计算面积的公式（2）判断多边形闭合坐标链走向）判断多边形闭合坐标链走向计算出面积为正，坐标链走向顺时针；计算出面积为正，坐标链走向顺时针；计算出面积为负，坐标链走向逆时针。计算出面积为负，坐标链走向逆时针。（3）判断线段（矢量）和空间的关系）判断线段（矢量）和空间的关系计算计算ABPA的面积，如为正，坐标链走向顺时针，的面积，如为正，坐标链走向顺时针，P在右面在右面计算计算AB

7、PA的面积，如为负，坐标链走向逆时针，的面积，如为负，坐标链走向逆时针，P在左面在左面辛普森公式的通用性辛普森公式的通用性 。 P4) 重重(质质)心的计算心的计算 区域的质心是目标保持均匀分布的平衡点，它通过对目标坐标值加区域的质心是目标保持均匀分布的平衡点，它通过对目标坐标值加权平均求得权平均求得 。质心可能在多边形内，也可能在多边形外。质心可能在多边形内，也可能在多边形外。 质心通常是指一个多边形或面积的几何中心。在有些情况下质心质心通常是指一个多边形或面积的几何中心。在有些情况下质心描述的不是几何中心，而是分布中心，或加权的平均中心。描述的不是几何中心，而是分布中心，或加权的平均中心。

8、 xicx= yii=1nni=1nncy=x , y分别是多边形顶点坐标分别是多边形顶点坐标; n 为顶点数为顶点数;cx,cy为多边形的重心为多边形的重心,也称质心也称质心,注意重心不一定在多边型内。注意重心不一定在多边型内。 2、线状目标的基本操作运算、线状目标的基本操作运算（1）直线）直线AB和和CD之间求交的计算之间求交的计算（2）曲线光滑的处理）曲线光滑的处理（3）曲线从简的处理删除多余点）曲线从简的处理删除多余点（4）曲线平行线的计算）曲线平行线的计算（5）直角平差的计算方法）直角平差的计算方法 如将直角平差。如将直角平差。3、多边形目标的基本操作运算、多边形目标的基本操作运算（

9、1）点在多边形内的判断）点在多边形内的判断（2）线与多边形交点的计算）线与多边形交点的计算（3）多边形之间交点的计算）多边形之间交点的计算（4）多边形内实现区域填充的算法）多边形内实现区域填充的算法 （2）线与多边形交点的计算）线与多边形交点的计算线与多边形相交否线与多边形相交否两点在多边形内否两点在多边形内否线与多边形相交否线与多边形相交否求线段与多边形的交点求线段与多边形的交点端点在多边形内又端点在多边形内又与多边形不相交与多边形不相交线在多边形内线在多边形内NY相交相交不相交不相交端点在多边形内又端点在多边形内又与多边形不相交与多边形不相交线与多边形相离线与多边形相离用来判断线与多边形是

10、否分离、包含、相交用来判断线与多边形是否分离、包含、相交（3）多边形之间交点的计算）多边形之间交点的计算 用来判断两个多边形是分离、包含、相交、相邻用来判断两个多边形是分离、包含、相交、相邻判断两个多边形的所有线段间是否有相交判断两个多边形的所有线段间是否有相交有相交无相交有公共边界否有公共边界否NY相邻相邻分离分离包含包含相交相交多边形分离还是包含多边形分离还是包含（4）多边形内实现区域填充的算法）多边形内实现区域填充的算法矢量数据结构的多边形区域内填充包括符号、颜色、晕线填充。矢量数据结构的多边形区域内填充包括符号、颜色、晕线填充。填充的算法主要有两大类填充的算法主要有两大类1、种子填充、

11、种子填充以种子点向外填充以种子点向外填充2、扫描线算法、扫描线算法主要分求交、排序、填充主要分求交、排序、填充6.3 空间数据的查询分析空间数据的查询分析一、空间数据查询的类型一、空间数据查询的类型空间数据的查询属于咨询式分析，是空间数据的查询属于咨询式分析，是GIS的最基本最重要的功能的最基本最重要的功能，也是，也是GIS面向用户的直接窗口，通常要求交互式进行，其结果面向用户的直接窗口，通常要求交互式进行，其结果要求动态地通过两个视窗进行显示。要求动态地通过两个视窗进行显示。1、基于属性（非空间）特征的查询、基于属性（非空间）特征的查询属性特征的查询主要在属性数据库中完成，这种查询通常基于标

12、属性特征的查询主要在属性数据库中完成，这种查询通常基于标准的准的SQL查询语言查询语言实现，之后按照属性数据和空间数据的对应关实现，之后按照属性数据和空间数据的对应关系显示图形。系显示图形。101102103104标号标号 植被植被 土壤土壤101 林地102 农地103 林地104 工业地2、基于空间特性的查询、基于空间特性的查询该查询通常指以图形图象或符号为语言元素的可视化查询。目前该查询通常指以图形图象或符号为语言元素的可视化查询。目前GIS都能实现对空间实体的简单查询，并显示该实体的属性，进行都能实现对空间实体的简单查询，并显示该实体的属性，进行相关的统计分析。相关的统计分析。如如 空

13、间几何数据查询，查询空间目标的坐标点、线长、面积、周空间几何数据查询，查询空间目标的坐标点、线长、面积、周长等基本参数。长等基本参数。 空间位置查询空间位置查询 简单的空间点击查询，区域（圆、多边形、矩形）简单的空间点击查询，区域（圆、多边形、矩形）内的查询。内的查询。空间关系查询空间关系查询 基于拓朴关系的查询如基于拓朴关系的查询如 邻接性查询、邻接性查询、 包含性查询、包含性查询、 穿越性查询穿越性查询 、 落入性查询、落入性查询、 方向性查询方向性查询 等。等。3、结合空间特性和非空间（属性）特征的查询、结合空间特性和非空间（属性）特征的查询 这种查询涉及空间特征和属性特征。如查长江以南

14、（空间）、人这种查询涉及空间特征和属性特征。如查长江以南（空间）、人口数大于口数大于100万的万的 （查属性）同北京的距离（空间距离）（查属性）同北京的距离（空间距离） 小于小于1000公里的城市。公里的城市。二、空间数据查询的方法二、空间数据查询的方法1）SQL查询查询如如 Select 所需数据项所需数据项 From 属性表属性表 Where 条件表达式条件表达式3）基于空间查询语言（）基于空间查询语言（Spatial Query Language)的查询的查询用标准查询语言用标准查询语言SQL，无法表达空间关系及空间运算操作，为此对，无法表达空间关系及空间运算操作，为此对原有的原有的SQ

15、L查询进行扩充或改造，从而提出了查询进行扩充或改造，从而提出了空间查询语言空间查询语言称扩充称扩充SQL和空间和空间SQL。扩充。扩充SQL或或空间查询语言均未空间查询语言均未标准化。标准化。如如Mapinfo提供的提供的空间查询例空间查询例从世界地图从世界地图world中查中国的邻国有哪些中查中国的邻国有哪些?空间查询语空间查询语 句句world.obj intersect(select obj from world where country=“China”) 基于空间操作基于空间操作6.4 空间数据的叠置分析空间数据的叠置分析一、叠置的分类一、叠置的分类1、按叠置的实质分、按叠置的实质分

16、 1）视觉叠置；）视觉叠置；2）条件叠置）条件叠置; 3）无条件叠置）无条件叠置2、按叠置的对象类型分、按叠置的对象类型分 1）点和面的叠置；）点和面的叠置；2）线和面的叠置；）线和面的叠置；3）面和面的叠置。）面和面的叠置。3、按叠置对象的数据结构分、按叠置对象的数据结构分 1）矢量叠置；）矢量叠置；2）栅格叠置）栅格叠置二、叠置分析的数学基础二、叠置分析的数学基础 空间逻辑运算空间逻辑运算为讨论方便将空间图层为讨论方便将空间图层A，B，C定义为二值图象定义为二值图象1、空间逻辑并（或）运算；、空间逻辑并（或）运算；AB =X XA XA 或或 X B X B 2、空间逻辑交（与）运算；、空

17、间逻辑交（与）运算；AB = X XA XA 且且 XB XB 3、空间逻辑差运算；、空间逻辑差运算；A - B =X XA XA 且且 XB XB 4、空间包含；、空间包含；A B三、视觉叠置分析三、视觉叠置分析视觉叠置分析的实质是将不同层面的图信息进行叠加显示，以视觉叠置分析的实质是将不同层面的图信息进行叠加显示，以判断其间的空间关系。判断其间的空间关系。视觉叠置分析的对象种类：视觉叠置分析的对象种类：1、不同图层（如点、线、面）自身及之间的叠加显示；、不同图层（如点、线、面）自身及之间的叠加显示；2、RS和专题图之间的叠加显示；和专题图之间的叠加显示；3、数字高程图同各种图之间的叠加显示

18、。、数字高程图同各种图之间的叠加显示。 四、基于矢量数据结构的叠置分析四、基于矢量数据结构的叠置分析1、点与多边形的叠置、点与多边形的叠置*1 *3 *4 *2A DB C点点 属性属性1234多边形多边形 属性属性1 属性属性2ABCD点点 多边形多边形 点属性点属性 面属性面属性1 面属性面属性21 A2 D3 C4 B2、线与多边形的叠置、线与多边形的叠置 新线ID 老线ID 多边形 线属性 多边形属性 12 22” 213 3AABCD1233、多边形之间的叠置、多边形之间的叠置 原始多边形的属性原始多边形的属性新多边形的属性新多边形的属性AXYZ新多边形新多边形 层层1多边形属性多边

19、形属性 层层2多边形属性多边形属性1 A 02 0 X3 A X4 0 X5 A 06 A 07 A Y8 A Z9 0 Y10 0 Z11 0 YID 属性属性101 AID 属性属性 Y Z 1 X37910145611 2层层1层层2新层新层层层1层层2多边形叠置流程多边形叠置流程 取本底多边形取本底多边形取上叠多边形取上叠多边形比较两个多边形相交否比较两个多边形相交否两个多边形是包含关系否两个多边形是包含关系否多边形求交多边形求交,重建拓朴关系重建拓朴关系建立包含关系建立包含关系属性抽取建立关联连接属性抽取建立关联连接上叠多边形结束否上叠多边形结束否 END本底多边形结束否本底多边形结

20、束否NYNY属性抽取属性抽取五、基于栅格数据结构的叠置分析五、基于栅格数据结构的叠置分析 AB123A1A2A3B1B2标号 地貌A 阳坡B 阴坡标号 植被1 林地2 农地3 牧地标号 综合属性A1 阳坡 林地A2 阳坡 农地A3 阳坡 牧地B1 阴坡 林地B2 阴坡 农地1、叠置条件的确定叠置条件的确定根据求解问题确定关系表达式，如根据水稻种植条件，求某区域适根据求解问题确定关系表达式，如根据水稻种植条件，求某区域适合种植水稻的地区。合种植水稻的地区。积温积温 3200度度;降雨量降雨量 800mm;坡度坡度200 天；天；关系表达式：关系表达式：（积温（积温 3200） （降雨量（降雨量

21、800） （坡度（坡度200） 2、叠置算法的确定叠置算法的确定（二值非权重二值非权重模型模型）单个栅格为单元的叠置；单个栅格为单元的叠置；游程编码叠置。游程编码叠置。设求解：某地区降雨量设求解：某地区降雨量=1000mm,土厚土厚=50cm的区域的区域 降雨量图降雨量图 土厚图土厚图游程号 降雨游程 游程最右列（终止码）1 600mm 2002 800mm 4003 1000mm 6804 1200mm 800游程号 土厚游程 游程最右列（终止码）1 30cm 170 2 20cm 3603 50cm 8000 400 800 0 400 80030cm1000mm1200mm800mm50

22、cm30cm40cm20cm600mm第K行属性的获取属性的获取游程号游程号 降雨土厚游程降雨土厚游程 游程最右列（终止码）游程最右列（终止码）1 600mm 30cm 170 2 600mm 20cm 2003 800mm 20cm 3604 800mm 50cm 4005 1000mm 50cm 6806 1200mm 50cm 800第第K行条件叠置后的游程编码表行条件叠置后的游程编码表游程号游程号 降雨土厚游程降雨土厚游程 游程最右列（终止码）游程最右列（终止码）1 0 4002 1000mm 50cm 6803 0 800六、二值权重模型六、二值权重模型 Wi* Ci S = Wi这

23、样一来，可使区域内任意一点的输出值在这样一来，可使区域内任意一点的输出值在0和和1之间，从而进行分之间，从而进行分级。级。nnii如如C1：积温：积温 3200度度;C2：降雨量：降雨量 800mm;C3：坡度：坡度200 天；天；6.5 空间数据的空间数据的邻域分析邻域分析一、一、缓冲区分析缓冲区分析1、缓冲区的类型、缓冲区的类型 点缓冲区；线缓冲区；面缓冲区；复杂缓冲区；点缓冲区；线缓冲区；面缓冲区；复杂缓冲区；可变距离缓冲区。可变距离缓冲区。2、缓冲区的建立、缓冲区的建立 点缓冲区的建立点缓冲区的建立 等距离的点缓冲区是一个圆；等距离的点缓冲区是一个圆；线缓冲区的建立；线缓冲区的建立；面

24、缓冲区的建立。面缓冲区的建立。*3、线缓冲区和面缓冲区建立的基本算法、线缓冲区和面缓冲区建立的基本算法1）角平分线法）角平分线法角平分线法的实质是通过求平行线建立线缓冲区和面缓冲区角平分线法的实质是通过求平行线建立线缓冲区和面缓冲区在求算过程中对出现的凸角作特殊处理。在求算过程中对出现的凸角作特殊处理。如下如下凸角出现过长的尖角凸角出现过长的尖角2）凸角圆弧法）凸角圆弧法 ABC3、缓冲区分析的过程、缓冲区分析的过程1）建立缓冲区）建立缓冲区;建立缓冲区实际上要复杂得多，一方面常规方法建立的缓冲区会建立缓冲区实际上要复杂得多，一方面常规方法建立的缓冲区会形成很多重叠区，必须解决；另一方面实际上

25、常要求建立不同距形成很多重叠区，必须解决；另一方面实际上常要求建立不同距离的缓冲区。离的缓冲区。2）进行进一步的分析）进行进一步的分析,如叠置分析或图形的运算。如叠置分析或图形的运算。4、缓冲区分析例、缓冲区分析例如已知一湖泊，要求在它周围如已知一湖泊，要求在它周围5000m 内必需禁止任何污染性工业企内必需禁止任何污染性工业企业存在，在它周围业存在，在它周围500m 内必需禁止建筑任何永久性建筑物。内必需禁止建筑任何永久性建筑物。（1）先建立缓冲区；）先建立缓冲区；（2）同现有污染性工业企业图叠置，显示在范围内应禁止的污染）同现有污染性工业企业图叠置，显示在范围内应禁止的污染性工业企业；性工

26、业企业；（3）同现有永久性建筑物图叠置，显示在范围内应禁止的永久性）同现有永久性建筑物图叠置，显示在范围内应禁止的永久性建筑物。建筑物。二、泰森多边形分析二、泰森多边形分析泰森多边形分析是菏兰气象学家泰森多边形分析是菏兰气象学家A.H.Thiessen提出的一种分析方提出的一种分析方法，最初用于计算平均降雨量。法，最初用于计算平均降雨量。泰森多边形由一批具有一定分布的样本点数据生成。其实质是将泰森多边形由一批具有一定分布的样本点数据生成。其实质是将每个样本点周围划出一个地块，认为这地块可以用此样本点数据每个样本点周围划出一个地块，认为这地块可以用此样本点数据作为平均数据，代表该地块的属性。如用

27、分布在各地的气象站测作为平均数据，代表该地块的属性。如用分布在各地的气象站测得的降雨量，得到该地区的降雨量分布图。得的降雨量，得到该地区的降雨量分布图。 三角形内任一点到顶三角形内任一点到顶点的距离点的距离 Ai* Ri R = Aii=1i=1nn平均降雨量平均降雨量6.6 空间数据的网络分析空间数据的网络分析(Network Analysis）GIS最早的应用领域是土地管理和城市规划，之后，在城市交通规最早的应用领域是土地管理和城市规划，之后，在城市交通规划管理划管理,城市地下管线管理城市地下管线管理,城市电力供应，城市通信等部门，开始城市电力供应，城市通信等部门，开始大量使用大量使用GI

28、S。这些地方使用的数据是以点和线组成的网状数据，为全面的描述网这些地方使用的数据是以点和线组成的网状数据，为全面的描述网状事物及其相互之间的关系，形成一类空间分析称网络分析。状事物及其相互之间的关系，形成一类空间分析称网络分析。从数学的角度看，网络分析的基础是图论和运筹学；从数学的角度看，网络分析的基础是图论和运筹学；从数据结构的角度，网络分析的基础是非线性图数据结构。从数据结构的角度，网络分析的基础是非线性图数据结构。网络分析的实质是通过研究网络的状态，模拟和分析资源在网络上网络分析的实质是通过研究网络的状态，模拟和分析资源在网络上的流动和分配，对网络结构及资源等的优化问题进行研究。的流动和

29、分配，对网络结构及资源等的优化问题进行研究。一、空间网络分析基础一、空间网络分析基础1、网络图论概念网络图论概念有向图和无向图有向图和无向图 v1v4v3v2v50 1 1 1 11 0 1 0 01 1 0 1 01 0 1 0 11 0 0 1 0e1e2e3e4e5e6e7vv1 0 1 0 1 0 11 1 0 0 0 0 0 0 1 1 1 0 0 00 0 0 1 1 1 00 0 0 0 0 1 1ev邻接矩阵邻接矩阵 关联矩阵关联矩阵 2、空间网络分析的类型、空间网络分析的类型空间网络的拓朴分类空间网络的拓朴分类二维平面网络二维平面网络非平面网络非平面网络线型流系统线型流系统线

30、型栅格系统线型栅格系统线型立体系统线型立体系统道道路路型型树树型型交交错错型型环环型型细细胞胞型型3、空间网络分析的方法、空间网络分析的方法1）路经分析（）路经分析（path analysis)从网络模型的角度看最佳路经问题是求解网络中任意两点之间从网络模型的角度看最佳路经问题是求解网络中任意两点之间阻抗强度最小的路经。阻抗强度最小的路经。大量用于交通、救灾、抢险、消防及信息传输中。大量用于交通、救灾、抢险、消防及信息传输中。2）地理编码）地理编码二、最佳路径问题二、最佳路径问题 路径分析来源于图论中的网络优化。路径分析来源于图论中的网络优化。Dijkstra算法（比较经典的算法）算法（比较经

31、典的算法）1959年由年由E.W.Dijkstra提出提出.Dijkstra算法执行过程算法执行过程: t1 t2 t3 t4 v1 0 v2 (v1,v2), 10 v3 ( v1, v2, v3), 60 ( v1, v4, v3), 50 v4 (v1,v4), 30 (v1,v4), 30 v5 (v1,v5), 100 (v1,v5), 100 (v1,v4,v5), 90 (v1,v4,v3,v5), 60 v2, 10 v4, 30 v3, 60 v5 ,60V1V4V3V5V21030605020100 0 10 30 100 30 100 0 50 0 50 0 10 0 1

32、0 20 0 60 20 0 60 0 0有向图的有向图的邻接矩阵邻接矩阵三、地理编码（三、地理编码（Geocoding)1、地理编码的概念、地理编码的概念2、 地理编码的种类地理编码的种类1）按地址地理编码）按地址地理编码2）按边界地理编码）按边界地理编码6.7 空间统计分析空间统计分析 空间统计分析空间统计分析变量筛选分析变量筛选分析变量聚类分析变量聚类分析主主因因子子分分析析 主主成成分分分分析析 关关键键变变量量分分析析判判别别分分析析 聚聚类类分分析析 空间分析综合应用例空间分析综合应用例例一、利用例一、利用GIS进行土地适宜性评价进行土地适宜性评价目标目标对玉米地进行土地适宜性评价

33、，通过评价将研究区分成四种适宜对玉米地进行土地适宜性评价，通过评价将研究区分成四种适宜性等级，性等级，S1、S2、S3、S4方法方法1）对评价对象的生态条件进行调查，得出结论：）对评价对象的生态条件进行调查，得出结论：玉米属于禾本科，为一年生草本植物；玉米属于禾本科，为一年生草本植物；主导生态条件是：喜高温、水量大、土壤肥、土层松、根系要防主导生态条件是：喜高温、水量大、土壤肥、土层松、根系要防水土流失；水土流失；2）确定玉米地土地适宜性评价的主要因素和因子及等级标准；）确定玉米地土地适宜性评价的主要因素和因子及等级标准；3）确定各因素的权重和贡献函数值；）确定各因素的权重和贡献函数值；4）计

34、算玉米地土地适宜性评价等级，并输出分级图。）计算玉米地土地适宜性评价等级，并输出分级图。主要因素和因子主要因素和因子 因素因素F1供水供水因素因素F2供肥供肥因素因素F3供氧供氧因素因素F4土壤侵蚀土壤侵蚀因素因素5温度温度水系图水系图 土壤图土壤图数字化数字化 数字化数字化缓冲区缓冲区 属性分级属性分级 叠置分析叠置分析 供水分级供水分级 因子因子1 因子因子1一级一级 离水源离水源 50cm二级二级 离水源离水源 100-200m 50-20cm三级三级 离水源离水源 200-500m 20-10cm四级四级 离水源离水源 500m 10cm 最后求得供水等级、供氧等级、土壤侵蚀等级的分级标准，最后求得供水等级、供氧等级、土壤侵蚀等级的分级标准，（温度反映在季节上）。（温度反映在季节上）。3）确定各因素的权重和贡献函数值，建立评价指标表）确定各因素的权重和贡献函数值，建立评价指标表（1）建立