第3章--空间数据的处理.pptVIP

其转换公式为：⊿X=(xmax-xmin)/J⊿Y=(ymax-ymin)/I(5-4)式中：⊿X，⊿Y分别表示每个栅格单元的边长。xmax,xmin分别表示矢量坐标中x的最大值和最小值。ymax,ymin分别表示矢量坐标中y的最大值和最小值。I，J分别表示栅格的行数和列数。例如：已知某一地区x方向为15km，y方向为30km，现要把该地区的地块图转换成栅格数据。要求栅格的分辨率为30mx30m，则由式5-4可知：行数I=30km/30m=1000列数J=15km/30m=5001.2线的转换由于曲线可用折线来表示，也就是当折线上取点足够多时，所画的折线在视觉上成为曲线。因此，线的变换实质上是完成相邻两点之间直线的转换。若已知一直线AB其两端点坐标分别为A(x1,y1)和B(x2,y2)，则其转换过程不仅包括标点A，B分别从点矢量数据转换成栅格数据，还包括求出直线AB所经过的中间栅格数据。其过程如下：(1)利用上述点转换法，将点A(x1,y1)，B(x2,y2)分别转换成栅格数据，求出相应的栅格的行列值。(2)由上述行列值求出直线所在行列值的范围。(3)确定直线经过的中间栅格点。若从直线两端点转换中，求出该直线经过的起始行号为i1,终止行号im,其中间点行号必定i2,i3,....im-1。接着求出相应行号相交于直线的列号。第七节空间数据的插值1.点插值点插值针对连续空间，用连续的平滑的数学面描述，分为：整体拟合：模型由研究区域内所有采样点上的全部特征观测值建立。如趋势面拟合技术。局部拟合：用邻近的数据点来估计未知点的值。在GIS中大量使用。1.1趋势面拟合技术假定数据点的空间坐标(x,y)为独立变量，而表示特征值的z坐标为因变量。当数据为一维时，可用回归线近似表示：z=a0+a1xN4N2弧上行，左-右=1-2=-1，弧左边栅格值加-1N4N3弧下行，右-左=3-2=1，弧右边栅格值加1。第四节栅格数据向矢量数据的转换转换的实质是将相同属性代码的栅格集合转变成由少量数据组成的边界弧段以及区域边界的拓扑关系。1.栅格数据矢量化的过程1.1前处理目的是把栅格图预处理成近似线划图的二值图形(0和1)，使每条线只有一个象元宽度。1.2矢量化将线划图从栅格数据转变成直角坐标数据，并通过边界线搜索生成拓扑关系，建立与属性数据的关系。1.3后处理目的一方面除去栅格矢量化产生的多余点，减少数据的冗余；另一方面采用插值算法对曲线进行光滑处理。下面以双边界直接搜索法为例，说明栅格数据向矢量数据矢量化的一种方法。2.双边界直接搜索法该方法的基本思想是通过边界提取，将左右多边形信息保存在边界点上。每条边界弧段由两个并行的边界链组成，分别记录该边界弧段的左右多边形号。2.1双边界直接搜索法的原则以2x2栅格窗口为单位，在单位窗口内四个栅格数据的模式可确定下一个窗口的搜索方向及弧段的拓扑关系。搜索的原则如图所示。其要点为：若4个栅格仅有2个不同编号且对角上编号不相同，为边界点。[2]有3个或4个不同的编号为结点。[3]有2个不同编号，但对角线编号两两相同，也是结点。2.2双边界直接搜索法的步骤如下:①提取结点和边界线。用2x2栅格窗口沿行列方向扫描全图，遇到边界点及结点栅格窗口时将该窗口内栅格元素作出标记，结点标识符前加负号;②边界线跟踪及左右多边形信息的获取。在提取结点和边界线基础上，通过边界的搜索获取结点文件、弧段文件及多边形文件。记录边界点的两个多边形编号作为被搜索边界的左右多边形号。第五节多元空间数据的融合GIS技术经过近40年的发展和应用，已经积累了大量的数据资源。但是，由于地理数据的多语义性、多时空型、多尺度性、获取手段的多样性、存储格式的不同以及数据模型与数据结构的差异等，，导致多元数据的产生，给数据的继承和信息共享困难。为了实现空间数据的共享，特别是随因特网的发展、数字地球的兴起和GIS应用的日益深入，多元数据的融合已成为GIS设计者和用户的共同要求。1.遥感与GIS数据的融合1.1最常见的几种融合方法1.1.1遥感图像与图形的融合。经过正射纠正后的遥感影像，与数字地图信息融合，可产生影像地图。这种影像地图具有一定的数字基础，有丰富的光谱信息与几何信息，又有行政界限和属性信息，直接提高了用户的可视化效果。1.1.2遥感数据与DEM的融合。DEM代表精确的地形信息，它与遥感数据