第十二讲：空间分析.pptVIP

十二讲: 空间缓冲分析(Buffer)(P138) 中山大学 遥感与地理信息工程系 主讲人：张新长 缓冲区分析 邻近度（Proximity）描述了地理空间中两个地物距离相近的程度，其确定是空间分析的一个重要手段。 交通沿线或河流沿线的地物有其独特的重要性，公共设施的服务半径，大型水库建设引起的搬迁，铁路、公路以及航运河道对其所穿过区域经济的发展的重要性等，均是一个邻近度问题。 缓冲区分析是解决邻近度问题的空间分析工具之一。 缓冲区定义 所谓缓冲区就是地理空间目标的一种影响范围或服务范围。 从数学的角度看，缓冲区分析的基本思想是给定一个空间对象或集合，确定它们的邻域，邻域的大小根据邻域半径R值确定。 缓冲区定义 对象 O i 缓冲区可定义为： B i = {x:d ( x, O i )=R} 既对象的半径为R的缓冲区为对象距离d小于或等于R的全部点的集合。d一般是欧氏距离，也可以是其他定义的距离。对于对象集合 O = {O i : i=1,2,…,n} 其半径为R的缓冲区是各个对象缓冲区的并。 i=1 n B= U B i 实例 另外，还有特殊形态的缓冲区，如点对象有三角形、矩形、圆形；线对象有双侧对称，双侧不对称或是单侧缓冲区，对于面对象有内侧和外侧缓冲区。 栅格缓冲区的建立 将栅格数据表示为一个二值(0,1)矩阵(M×N)，其中“0”像元为空白位置，“1”元素为空间物体所占据的位置。经过距离变换，计算出每个“0”元素与最近的“1”元素的距离，即背景像元与空间物体的最小距离。假设缓冲区的宽度为d，则缓冲区边界就是距离为d的各个背景像元的集合。 距离的计算 某像元Pij与“1”像元的欧氏距离的计算可通过其行号差aij与列号差bij得到，dij = 。欧氏距离变换的方法是，首先设“1”元素Pij的aij=bij=0，设“0”元素的aij=bij=MAX(M, N)；然后计算各个像元及其周围8个像元的欧氏距离值并刷新aij和bij值，这时aij和bij表示了该像元与最邻近的“1”像元的行号差及列号差；最后通过公式dij= 计算它与“1”像元（空间物体）的最小距离。 距离的计算 (1)对Pij=1的所有像元，置aij=bij=0，否则置aij=bij=MAX(M,N)。 (2)按照从上到下，从左到右的次序计算dij并刷新aij、bij的值： 首先计算d’ij： 1 0 距离的计算 d’ij=dij= d’i-1,j-1= d’i-1,j = d’i-1,j+1= d’i,j-1= d’min=min(d’ij,d’i-1,j-1,d’i-1,j,d’i-1,j+1,d’i,j-1) 距离的计算 然后刷新aij,bij： 距离的计算 (3)类似于(2)，按从下到上，从右到左的次序计算dij并刷新aij，bij的值。 首先计算d’ij： d’ij=dij= d’i+1,j+1= d’i+1,j = d’i+1,j-1= d’i,j+1= 距离的计算 然后刷新aij,bij： 距离的计算 (4)对任一像元Pij，计算其距离值，即 dij = 距离的计算 欧氏距离变换的精度受栅格尺寸的影响，可以通过减小栅格的尺寸而获得较高的精度。其计算速度也较快。 栅格方法原理简单，但精度较低，而且内存开销较大，难以实现大数据量的缓冲区分析。 矢量缓冲区的建立 矢量缓冲区常见的有角平分法和叠置算法。角平分法由三步组成，即逐个线段计算简单平行线，尖角光滑矫正和自相交处理。尖角光滑矫正除角平分线法之外，还可采取圆弧法，但矫正过程都很复杂，难以完全地实现。叠置方法分两步完成。首先求出点、线段等基本元素的缓冲区，然后通过对基本元素缓冲区的叠置运算，求解折线、面边界等复杂目标的缓冲区。 角分线法 难以保证双线的等宽性 叠置算法-基元的生成 假设圆半径是r,A的坐标为(Ax,Ay), B的坐标为(Bx,By)。 AB的倾角α=arctg((By-Ay)/(Bx-Ax))。Δx=|BD|=rsinα, Δy=|Db|=rcosα。基元矩形框顶点a,b,c,d的坐标为： a.x=Ax+rsinα; a.y=Ay- rcosα b.x=Bx+rsinα; b.y=By- rcosα c.x=Bx- rsinα; c.y=By+rcosα d.x=Ax- rsinα; d.y=Ax+rcosα 叠置算法-基元的合并 方法是在交点处将基元边界元素分裂打断，再判断其是否落入其它基元内部，并删除落入基元内部的边界元