第2 空间邻近度分析

1、第第2 2部分部分 空间邻近度分析空间邻近度分析空间缓冲区类型矢量缓冲区的建立栅格缓冲区确定空间缓冲区半径的模型Voronoi多边形分析要求掌握的：o 定义o 缓冲区分类o 缓冲区建立方法、步骤o 生成缓冲区注意的问题o 上机怎样生成缓冲区o 泰森多边形一、矢量数据的缓冲区分析1、定义o 邻近度邻近度（ProximityProximity）描述了地理空间中两个地物距离相近的）描述了地理空间中两个地物距离相近的程度度，其确定是空间分析的一个重要手段。程度度，其确定是空间分析的一个重要手段。o 交通沿线或河流沿线的地物有其独特的重要性，公共设施的服交通沿线或河流沿线的地物有其独特的重要性，公共设施

2、的服务半径，大型水库建设引起的搬迁，铁路、公路以及航运河道务半径，大型水库建设引起的搬迁，铁路、公路以及航运河道对其所穿过区域经济的发展的重要性等，都是一个邻近度问题。对其所穿过区域经济的发展的重要性等，都是一个邻近度问题。o 所谓缓冲区就是地理空间目标的一种影响范围或服务范围。所谓缓冲区就是地理空间目标的一种影响范围或服务范围。o 缓冲区分析是解决邻近度问题的空间分析工具之一。缓冲区分析是解决邻近度问题的空间分析工具之一。o 从数学的角度看，缓冲区分析的基本思想是给定一个空间对象从数学的角度看，缓冲区分析的基本思想是给定一个空间对象或集合，确定它们的邻域，邻域的大小有邻域半径或集合，确定它们

3、的邻域，邻域的大小有邻域半径R R确定确定o 空间缓冲区分析的方法 方案一 方案二计算模型参数f0和d0设定di计算ri和Fi缓冲区图形表达计算模型参数f0和d0设定Fi计算di缓冲区图形表达1 缓冲区定义 对象 O i 缓冲区可定义为： B i = x:d ( x, O i )=R既对象的半径为R的缓冲区为距离对象距离d小于或等于R的全部点的集合。d一般是欧氏距离，也可以是其他定义的距离。对于对象集合O = O i : i=1,2,n其半径为R的缓冲区是各个对象缓冲区的并。i=1n点的缓冲区线的缓冲区面的缓冲区二、缓冲区分类三三 矢量缓冲区的建立矢量缓冲区的建立 矢量缓冲区常见的有角平分法和

4、叠置角平分法和叠置算法。角平分法由三步组成，即逐个线段计算简单平行线，尖角光滑矫正和自相交处理。尖角光滑矫正除角平分线法之外，还可采取圆弧法，但矫正过程都很复杂，难以完备地实现。叠置方法分两步完成。首先求出点、线段等基本元素的缓冲区，然后通过对基本元素缓冲区的叠置运算，求解折线、面边界等复杂目标的缓冲区。 1 角分线法难以保证双线的等宽性2 凸角圆弧法BdRd=r/sin(B/2)保证了平行曲线的等宽性折点凸凹性的判断把相邻两个线段看成两个向量，其方向取坐标点序方向。若前一个向量以最小角度转向第二个向量时呈逆时钟方向，则为凸顶点。反之为凹顶点。3 凸角圆弧法四四 栅格空间邻域运算栅格空间邻域运

5、算（一）（一） 邻域大小与类型（分析窗口的类型）邻域大小与类型（分析窗口的类型）（二） 邻域统计值o 统计运算（平均值、最大值、最小值、中值、求和、标准差等）o 测度运算（范围统计、多数统计、少数统计、种类统计）o 函数运算（滤波、地形参数）o 追踪分析miminjnjijijijcfC)(统计： 最大值、最小值 平均值、标准差 值域、总和、模、测度： 多数、少数、种类、 （三）（三） 邻域运算函数邻域运算函数四四 邻域运算的应用：数据平滑邻域运算的应用：数据平滑低通滤波 平均值高通滤波 原图-低通滤波四四 邻域运算应用：滤波邻域运算应用：滤波（四）（四） 分区运算分区运算（五）（五） 自然距

6、离量测运算自然距离量测运算距离量算：计算规则距离量算：计算规则1.414 X Grid CellGrid Cell1 1 自然距离量算自然距离量算n单位：格网单元 n类型： l 连续距离对源单元建立缓冲，形成距离带 l 对格网中每个单元确定于最近源单元的自然距离 nArcgis实现 l 连续：Find Distance l 最近：Assign proximity自然距离量算例自然距离量算例o 在arcgis里生成距离图：（1）在空间分析下拉菜单选择“选项”，确定分析淹没即分析范围；（2）空间分析下拉菜单选择“距离直线”（3）在直线对话框选择像元大小以及输出图像名，确定；（4）对生成的距离图重分

7、类（可选可不选）移经每个每个单元的成本或阻抗作为距离单位的距离量移经每个每个单元的成本或阻抗作为距离单位的距离量测分析方法测分析方法 最小成本表面、最小成本路径最小成本表面、最小成本路径2 2 成本距离量算成本距离量算二二 成本距离量算成本距离量算12121415236713442124153673365.74.54.9 3.54.2成本距离计算：连接成本矩阵成本距离计算：连接成本矩阵成本距离计算成本距离计算BA1212141523671344源点矩阵成本矩阵最小累计成本计算示例最小累计成本计算示例31212141526713442.834.92.553.5 5.72.8链接矩阵最小累计成本计

8、算示例最小累计成本计算示例1.5B4.211.52.8A211.51.52.02.84.21212141523671344最小累计成本计算示例最小累计成本计算示例1.5B4.211.52.86.74.5A2411.51.52.02.84.21.51.52.02.84.06.74.512121415236713443.51.5B35.74.011.52.86.74.5A21.51.52.02.84.04.56.72.02.84.05.74.56.73.03.51212141523671344最小累计成本计算示例最小累计成本计算示例最小累计成本计算示例最小累计成本计算示例3.51.5B35.74.

9、011.52.86.74.5A25.52.84.05.74.56.73.03.512121415236713442.02.84.05.74.56.73.03.55.53.51.5B35.74.011.52.86.74.5A25.52.84.05.74.56.73.03.512121415236713445.54.05.74.56.73.03.55.5最小累计成本计算示例最小累计成本计算示例43.51.5B35.74.011.52.86.74.5A25.512121415236713444.05.74.56.73.03.55.54.05.54.56.73.55.55.5最小累计成本计算示例最小累

10、计成本计算示例43.51.5B35.74.011.52.86.74.5A25.59.5311B35173487A55712345678最小累计成本矩阵单元指派矩阵五五 应注意的问题应注意的问题 （1 1）缓冲区发生重叠时的处理）缓冲区发生重叠时的处理 缓冲区的重叠包括多个特缓冲区的重叠包括多个特征缓冲区之间的重叠以及同一特征缓冲区图形的重叠。对征缓冲区之间的重叠以及同一特征缓冲区图形的重叠。对于前者，首先通过拓扑分析的方法，自动地识别出落在某于前者，首先通过拓扑分析的方法，自动地识别出落在某个缓冲区内部的那些线段或弧段，然后删除这些线段或弧个缓冲区内部的那些线段或弧段，然后删除这些线段或弧段，

11、得到经处理后的连通缓冲区；对于后者，可通过缓冲段，得到经处理后的连通缓冲区；对于后者，可通过缓冲区边界曲线逐条线段求交。区边界曲线逐条线段求交。 （2 2）缓冲区宽度不同时的处理）缓冲区宽度不同时的处理 在进行缓冲区分析时，经在进行缓冲区分析时，经常发生不同级别的同一类要求具有不同的缓冲区大小。例常发生不同级别的同一类要求具有不同的缓冲区大小。例如，在城市土地地价评估时，沿主要街道两侧的通达度、如，在城市土地地价评估时，沿主要街道两侧的通达度、繁华度的辐射范围大，而小街道较小，这与要素的类型和繁华度的辐射范围大，而小街道较小，这与要素的类型和特点有关。在建立这种缓冲区时，首先应建立要素属性表，

12、特点有关。在建立这种缓冲区时，首先应建立要素属性表，根据不同属性确定不同的缓冲区宽度，然后产生缓冲区根据不同属性确定不同的缓冲区宽度，然后产生缓冲区 。 （3 3）复杂图形缓冲区的内外标）复杂图形缓冲区的内外标识处理识处理 复杂图形经缓冲复杂图形经缓冲区分析后会产生许多多边形。区分析后会产生许多多边形。为了标识哪些区域是缓冲区，为了标识哪些区域是缓冲区，哪些是缓冲区外，应在这些哪些是缓冲区外，应在这些多边形中加入特征属性。如多边形中加入特征属性。如ARC/INFOARC/INFO在缓冲区分析后的在缓冲区分析后的多边形属性表中加入了一项多边形属性表中加入了一项INSIDEINSIDE，INSID

13、EINSIDE值为值为1 1表示该表示该多边形在缓冲区内，值为多边形在缓冲区内，值为0 0则则在缓冲区外。在缓冲区外。 五上机生成缓冲区o 设置地图单位和距离单位o 在arctoolbox里生成缓冲区（只能生成单一缓冲区）o 添加缓冲区工具（可生成指定距离的缓冲区、多重缓冲区、根据某个字段记录里指定距离的缓冲区）六 确定空间缓冲区半径的模型o 线性函数模型o 幂函数模型o 指数函数模型七、泰森多边形分析七、泰森多边形分析 泰森多边形泰森多边形：首先由荷兰气象学家：首先由荷兰气象学家A.H.ThiessenA.H.Thiessen提出的，最提出的，最初用于根据离散分布气象站的降雨量来计算区域的平

14、均降雨量，初用于根据离散分布气象站的降雨量来计算区域的平均降雨量，现在已广泛应用于区域的地学分析和现在已广泛应用于区域的地学分析和GISGIS的其它应用中。的其它应用中。 基本方法基本方法：将所有相邻的气象站连成三角形，作这些三角：将所有相邻的气象站连成三角形，作这些三角形各边的垂直平分线，于是，每个气象站周围若干条垂直平分线形各边的垂直平分线，于是，每个气象站周围若干条垂直平分线便围成一个多边形，用这个多边形内所包含的一个唯一气象站的便围成一个多边形，用这个多边形内所包含的一个唯一气象站的降雨强度就可以用来表示这个多边形区域内的降雨强度，并称这降雨强度就可以用来表示这个多边形区域内的降雨强度

15、，并称这个多边形为泰森多边形。个多边形为泰森多边形。 如图所示，区域如图所示，区域D D内有几个离散点内有几个离散点Pi(Xi,Yi),iPi(Xi,Yi),i=1,2,=1,2,n n。将该区域用一组直线段分成几个互相。将该区域用一组直线段分成几个互相邻接的多边形，并且：邻接的多边形，并且： 1 1、每个多边形内含有且仅有一个离散数据点；、每个多边形内含有且仅有一个离散数据点； 2 2、若区域内的任意一点、若区域内的任意一点K(Xk,YkK(Xk,Yk) )位于离散点位于离散点（Xi,YiXi,Yi）的多边的多边形内，则不等式形内，则不等式 (Xk-Xi)(Xk-Xi)2 2+(Yk-Yi)

16、+(Yk-Yi)2 2 1/21/2(Xk-Xj)(Xk-Xj)2 2+(Yk-Yj)+(Yk-Yj)2 2 1/21/2 当当i i j j时，恒成立；时，恒成立； 3 3、若、若点点M(XmM(Xm、YmYm) )位于含离散点位于含离散点（Xi,YiXi,Yi）与含与含（Xj,YjXj,Yj）的两的两个多边形的公共边上，则有个多边形的公共边上，则有 (Xm-Xi)(Xm-Xi)2 2+(Ym-Yi)+(Ym-Yi)2 2 1/21/2=(Xm-Xj)=(Xm-Xj)2 2+(Ym-Yj)+(Ym-Yj)2 2 1/21/2 这样的多边形就叫泰森多边形，连接相邻多边形离散点，这样的多边形就叫泰森多边形，连接相邻多边形离散点，构成的三角形叫做泰森三角形构成的三角形叫做泰森三角形。 o 应当指出，泰森三角形和泰森多边形都不是唯一的，若应当指出，泰森三角形和泰森多边形都不是唯一的，若生成的泰森多边形面积差异过大，说明观测站点分布不生成的泰森多边形面积差异过大，说明观测站点分布不合理，需要调整。合理，需要调整。o 泰森多边形分析方法及构成的多边形和三角网可