GIS课程(空间数据处理).ppt

《GIS课程(空间数据处理).ppt》由会员分享，可在线阅读，更多相关《GIS课程(空间数据处理).ppt（136页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、地理信息系统原理GIS第三章 空间数据的处理第一节第一节空间数据的变换空间数据的变换第二节第二节空间数据结构的转换空间数据结构的转换第三节第三节多元空间数据的融合多元空间数据的融合第四节第四节空间数据的压缩与重分类空间数据的压缩与重分类第五节第五节空间数据的内插方法空间数据的内插方法第六节拓扑关系的编辑第六节拓扑关系的编辑地理信息系统原理GIS数据处理的内容数据变换：数据变换：指数据从一种数学状态到另一种数学状态的指数据从一种数学状态到另一种数学状态的变换，包括几何纠正、投影转换和辐射纠正等，以解变换，包括几何纠正、投影转换和辐射纠正等，以解决空间数据的几何配准。决空间数据的几何配准。数据重构

2、：数据重构：指数据从一种格式到另一种格式的转换，包指数据从一种格式到另一种格式的转换，包括结构转换、格式转换、类型替换等，以解决空间数括结构转换、格式转换、类型替换等，以解决空间数据在结构、格式和类型上的统一，实现多源数据和异据在结构、格式和类型上的统一，实现多源数据和异构数据的联接与融合。构数据的联接与融合。数据提取：数据提取：指对数据进行某种条件取舍，包括类型提取、指对数据进行某种条件取舍，包括类型提取、窗口提取、空间内插等，以适应不同用户对数据的特窗口提取、空间内插等，以适应不同用户对数据的特定需求。定需求。地理信息系统原理GIS第一节 空间数据的变换空间数据的变换即空间坐标系的变换。空

3、间数据的变换即空间坐标系的变换。实质：实质：是建立是建立两个平面点之间的一一对应关系，包括几何纠正和投两个平面点之间的一一对应关系，包括几何纠正和投影转换，它们是空间数据处理的基本内容之一。包括：影转换，它们是空间数据处理的基本内容之一。包括：几何纠正几何纠正和和投影转换投影转换。地理信息系统原理GIS地理信息系统原理GIS一、几何纠正一、几何纠正几何纠正：几何纠正：是为了实现对数字化数据的坐标系转换和是为了实现对数字化数据的坐标系转换和图纸变形误差的改正。图纸变形误差的改正。几何纠正一般包括：几何纠正一般包括：仿射变换、相似变换、二次变仿射变换、相似变换、二次变换换等功能。本教材主要介绍了等

4、功能。本教材主要介绍了仿射变换仿射变换：地理信息系统原理GIS特点：特点：1、直线变换后仍为直线。2、平行线变换后认为平行线3、不同方向上的长度比发生变化。仿射变换仿射变换仿射变换可以对坐标数据在仿射变换可以对坐标数据在x x和和y y方向进行不同比例的方向进行不同比例的缩放、旋转、平移。缩放、旋转、平移。地理信息系统原理GIS两层的数据不匹配两层的数据不匹配地理信息系统原理GISX方向Y 方向平移平移 旋转旋转 缩小缩小 地理信息系统原理GIS0yxP（x，y）P(x，y)xyx=x+xy=y+y1 1 1 1、平移变换、平移变换、平移变换、平移变换地理信息系统原理GISyP（x，y）0 x

5、P(X，Y)、旋转变换、旋转变换、旋转变换、旋转变换地理信息系统原理GISX=X=xcos-y sin Y=Y=xsin+y cosX=xX=x0 0+(x-x0)cos-(y-y y0 0)sinY=y0+(x-x0)sin+(y-y0)cos地理信息系统原理GIS点可以通过对其P P（x x，y y）坐标分别乘以各自的比例因子S Sx x和S Sy y来改变它们到坐标原点的距离。x=xSxy=ySyx=x0+(x-x0)Sx y=y0+(y-y0)Sy、比例变换（图形缩放）、比例变换（图形缩放）、比例变换（图形缩放）、比例变换（图形缩放）地理信息系统原理GIS仿射变换仿射变换转换坐标与理转

6、换坐标与理论坐标之差论坐标之差法方程组法方程组中位点误差中位点误差地理信息系统原理GIS二、投影投影及其转换二、投影投影及其转换（一）地图投影的基本原理（一）地图投影的基本原理就是建立地球椭球面上的点的地理坐标（，）与平面上对应点的平面坐标（x，y）之间的函数关系：地理信息系统原理GIS地图投影：投影变形 将不可展的地球椭球面展开成平面，并且不能有断裂，则图形必将在某些地方被拉伸，某些地方被压缩，故投影变形是不可避免的。长度变形长度变形 面积变形面积变形 角度变形角度变形（二）投影类型（二）投影类型地理信息系统原理GIS2.变形椭圆 取地面上一个微分圆（小到可忽略地球曲面的取地面上一个微分圆（

7、小到可忽略地球曲面的影响，把它当作平面看待），它投影到平面上通常影响，把它当作平面看待），它投影到平面上通常会变为椭圆，通过对这个椭圆的研究，分析地图投会变为椭圆，通过对这个椭圆的研究，分析地图投影的变形状况。这种图解方法就叫影的变形状况。这种图解方法就叫变形椭圆变形椭圆。为经线长度比；为经线长度比；为纬线长度比为纬线长度比地理信息系统原理GIS微小圆变形椭圆该方程证明该方程证明:地球面上的微小圆，地球面上的微小圆，投影后通常会变为椭圆，即：投影后通常会变为椭圆，即：以以O为原点，以相交成为原点，以相交成q q角的两共角的两共轭直径为坐标轴的椭圆方程式。轭直径为坐标轴的椭圆方程式。代入：代入：

8、X2+Y2=1，得得地理信息系统原理GIS特特别方向：方向：变形椭圆上相互垂直的两个方向及经向和纬向 长轴方向（极大值）长轴方向（极大值）a短轴方向（极小值）短轴方向（极小值）b经线方向方向m；纬线方向纬线方向n统称统称主方向主方向据据阿波隆尼定理阿波隆尼定理阿波隆尼定理阿波隆尼定理，有，有m2+n2=a2+b2mnsinq q=ab地理信息系统原理GIS地理信息系统原理GIS 在分析地图投影时，可借助对变形椭圆和在分析地图投影时，可借助对变形椭圆和微小圆的比较，说明变形的性质和大小。椭圆微小圆的比较，说明变形的性质和大小。椭圆半径与小圆半径之比，可说明长度变形。很显半径与小圆半径之比，可说明

9、长度变形。很显然，长度变形随方向的变化而变化，其中有一然，长度变形随方向的变化而变化，其中有一个极大值，即椭圆长轴方向，一个极小值，即个极大值，即椭圆长轴方向，一个极小值，即椭圆短轴方向。这两个方向是相互垂直的，称椭圆短轴方向。这两个方向是相互垂直的，称为主方向。椭圆面积与小圆面积之比，可说明为主方向。椭圆面积与小圆面积之比，可说明面积变形。椭圆上两方向线的夹角和小圆上相面积变形。椭圆上两方向线的夹角和小圆上相应两方向线的夹角的比较，可说明角度变形。应两方向线的夹角的比较，可说明角度变形。地理信息系统原理GIS变形椭圆与投影变形的关系 a=b：等角投影等角投影（投影面上某点的任意两方向线夹角与

10、椭球面上相应两线段夹角相等，即角度变形为零=0。制作航海图、洋流图、风向图）。ab=r*r:等积投影等积投影（投影面与椭球面上相应区域的面积相等，即面积变形为零。制作地质图、土壤图、行政区划图等）。其他：任意投影任意投影（投影图上，长度、面积和角度都有变形，它既不等角又不等积。其中，等距投影是在特定方向上没有长度变形的任意投影）。地理信息系统原理GIS地理信息系统原理GIS按投影面的形状分为按投影面的形状分为圆锥投影圆锥投影、圆柱投影圆柱投影和和方位投影；方位投影；按投影面与地球的相对位置关系分为按投影面与地球的相对位置关系分为正轴投影正轴投影、斜轴投斜轴投影影、和、和横轴投影横轴投影；按投影

11、面和地球的空间逻辑关系分为按投影面和地球的空间逻辑关系分为相切相切和和相割相割两类投两类投影。影。投影的分类投影的分类地理信息系统原理GIS根据地图投影变形情况，地图比例尺分为：主比例尺主比例尺：在投影面上没有变形的点或线上的比例尺。局部比例尺：局部比例尺：在投影面上有变形处的比例尺。图解直线比例尺图解直线比例尺图解复式比例尺图解复式比例尺地理信息系统原理GIS（三）地理信息系统常用的地图投影（三）地理信息系统常用的地图投影1 1、高斯、高斯克吕格投影克吕格投影地理信息系统原理GIS高斯克吕格投影高斯投影是一种横轴等角切椭圆柱投影，其条件为：1）中央经线和地球赤道投影成为直线且为投影的对 称轴

12、；2）等角投影；3）中央经线上没有长度变形。由公式可分析出高斯投影变形具有以下特点：1）中央经线上无变形 2)中央经线上的长度比为1，其他任何点上的长度比大于1.3）同一条纬线上，离中央经线越远，变形越大；4）同一条经线上，纬度越低，变形越大；5）投影属于等角性质，面积比为长度比的平方。6）等变形线为平行于中央经线的直线。地理信息系统原理GIS2 2 墨卡托墨卡托(Mercator)(Mercator)投影投影特点特点:1 1、无角度变形，但面积变、无角度变形，但面积变形较大。形较大。2 2、经线和纬线是两组相互、经线和纬线是两组相互垂直的平行直线，经线间隔垂直的平行直线，经线间隔相等，纬线间

13、隔由赤道向两相等，纬线间隔由赤道向两极逐渐扩大。极逐渐扩大。3 3、保持方向和相对位置的、保持方向和相对位置的正确。正确。属于等角正切圆柱投影。属于等角正切圆柱投影。地理信息系统原理GIS墨卡托投影（正轴等角圆柱投影墨卡托投影（正轴等角圆柱投影）地理信息系统原理GIS墨卡托投影（正轴等角圆柱投影墨卡托投影（正轴等角圆柱投影）墨卡托投影常用来制作航海图和航空图墨卡托投影常用来制作航海图和航空图地理信息系统原理GIS3 UTM3 UTM投影（横轴墨卡托投影）投影（横轴墨卡托投影）UTM投影是“等角横轴割圆柱投影”，圆柱割地球于南纬80度、北纬84度两条等高圈，投影后两条割线上没有变形，中央经线上长

14、度比 0.9996。（高斯高斯-克吕格投影是克吕格投影是“等角横切椭圆柱投影等角横切椭圆柱投影”）我国卫星地图、我国卫星地图、Google MapGoogle Map采用采用UTMUTM投影。投影。地理信息系统原理GIS4.4.兰勃特兰勃特(Lambert)(Lambert)投影投影兰勃特等角投影，在双标准纬线下是一兰勃特等角投影，在双标准纬线下是一“等角正轴等角正轴割圆锥投影割圆锥投影”.我国我国1 1：100100万地形图采用了兰勃特投影。万地形图采用了兰勃特投影。采用双标准纬线，投影变形小而均匀特点：采用双标准纬线，投影变形小而均匀特点：1 1、角度没有变形。、角度没有变形。2 2、两条

15、标准纬线上没有任何变形。、两条标准纬线上没有任何变形。3 3、等变形线和纬度一致，同一条纬线上的变形、等变形线和纬度一致，同一条纬线上的变形处处相等。处处相等。4 4、两标准纬线外侧为正变形，而两标准纬线之、两标准纬线外侧为正变形，而两标准纬线之间为负变形。间为负变形。5 5、同一纬线上等经差的线段长度相等，两条纬、同一纬线上等经差的线段长度相等，两条纬线间的经纬度长度处处相等。线间的经纬度长度处处相等。地理信息系统原理GIS5.5.阿尔伯斯阿尔伯斯(Albers)(Albers)投影投影阿尔伯斯阿尔伯斯(Albers)(Albers)投影是等面积正割圆锥投影，我投影是等面积正割圆锥投影，我国

16、部分省图采用这种投影国部分省图采用这种投影.地理信息系统原理GIS四、投影转换及其转换四、投影转换及其转换主要研究从一种地图投影变为另一种地图投影的理论和主要研究从一种地图投影变为另一种地图投影的理论和方法。其实质是建立两平面之间点的一一对应关系。方法。其实质是建立两平面之间点的一一对应关系。x=f1(x=f1(，)y=f2(y=f2(，)X X=f3(=f3(，)Y Y=f4(=f4(，)X=F1(x,y)X=F1(x,y)Y=F2(x,y)Y=F2(x,y)地理信息系统原理GIS投影变换：投影A（x，y）投影B（X，Y）正解变换：解析函数关系X=f(x,y)，Y=g(x,y)反解变换：经纬

17、度B=f(x,y),L=g(x,y)X=F(B,L),Y=G(B,L)从投影转换的方式来分：地理信息系统原理GIS、正解变换、正解变换、正解变换、正解变换通过建立资料地图的投影坐标通过建立资料地图的投影坐标数据到目标地图投影坐标数据的严数据到目标地图投影坐标数据的严密或近似的解析关系式，直接由资密或近似的解析关系式，直接由资料地图投影坐标数据料地图投影坐标数据x x,y y转换为目标转换为目标投影的直角坐标投影的直角坐标X X，Y Y。两个不同投影平面场上的点可对应两个不同投影平面场上的点可对应写成：写成：X X=f1(=f1(x x,y y)Y Y=f2=f2（x x,y y）式中式中f1f

18、1，f2f2为定域内单值、连续的为定域内单值、连续的函数。函数。X=f1(x,y)Y=f2（x,y）投影转换的方式有两种：投影转换的方式有两种：正解变换正解变换、反解变换反解变换地理信息系统原理GIS、反解变换、反解变换、反解变换、反解变换将资料地图的投影坐标数据将资料地图的投影坐标数据x x,y y反解出地理坐标反解出地理坐标，然后再将地理坐标代入到目标地图的投影坐标公，然后再将地理坐标代入到目标地图的投影坐标公式中，从而实现投影坐标的转换。式中，从而实现投影坐标的转换。对前后两种地图投影，可分别有如下表达形式：对前后两种地图投影，可分别有如下表达形式：x=f1(x=f1(，)y=f2()y

19、=f2(，)X X=f3(=f3(，)Y Y=f4(=f4(，)根据资料地图的投影公式求反解，对前一投影则有：根据资料地图的投影公式求反解，对前一投影则有：=f1f1 (x x,y y)=f2()=f2(x x,y y)代入目标地图的投影方程即有：代入目标地图的投影方程即有：X X=f3f1(=f3f1(x x,y y),f2(),f2(x x,y y)Y Y=f4f1(=f4f1(x x,y y),),f2(f2(x x,y y)这就是地图投影反解变换的数学模型。这就是地图投影反解变换的数学模型。地理信息系统原理GIS(，)（x x,y y）（X,X,Y Y）=f1f1 (x x,y y)=

20、f2(=f2(x x,y y)X X=f3(=f3(，)Y Y=f4(=f4(，)X X=f3f1(=f3f1(x x,y y),f2(),f2(x x,y y)Y Y=f4f1(=f4f1(x x,y y),f2(),f2(x x,y y)、反解变换、反解变换、反解变换、反解变换地理信息系统原理GIS投影转换的方法：投影转换的方法：（1 1）解析变换：）解析变换：投影坐标（投影坐标（x x、y y）变换为地理坐标（）变换为地理坐标（B B、L L），再由），再由地理坐标变换为另一种投影坐标（地理坐标变换为另一种投影坐标（X X、Y Y）。）。（2 2）数字变换：一种投影的数字化坐标）数字变换

21、：一种投影的数字化坐标x x、y y变换到另一种投影的直变换到另一种投影的直角坐标角坐标X X、Y Y。（3 3）解析）解析-数值变换：当原数据投影公式不知道的情况下，反解出数值变换：当原数据投影公式不知道的情况下，反解出地理坐标，再由代入已知的新数据投影公式中进行计算。地理坐标，再由代入已知的新数据投影公式中进行计算。地理信息系统原理GIS第二节 空间数据结构的转换优优点点 缺点缺点 矢矢量量数数据据结结构构 1.便于面向便于面向实实体的数据表达体的数据表达2.数据数据结结构构紧紧凑，冗余度低凑，冗余度低3.拓扑拓扑结结构有利于于网构有利于于网络络分析、分析、空空间间分析分析1.数据数据结结

22、构构处处理理较较复复杂杂2.软软件件实现实现的技的技术术要求比要求比较较高。高。3.多多边边形叠置等分析相形叠置等分析相对对困困难难；栅栅格格数数据据结结构构 1.数据数据结结构相构相对简单对简单，易于算法，易于算法实现实现；2.空空间间分析分析较较容易容易实现实现；3.有利于与遥感数据的匹配有利于与遥感数据的匹配应应用用和分析；和分析；1.图图形数据量大，冗余度高，需要形数据量大，冗余度高，需要压缩处压缩处理；理；2.定位精度比矢量低；定位精度比矢量低；3.拓扑关系拓扑关系难难以表达；以表达；地理信息系统原理GIS第二节 空间数据结构的转换 应用矢量数据与栅格数据的一般原则：应用矢量数据与栅

23、格数据的一般原则：数据采集采用矢量数据结构，有利于保证空间实体的几何精度和数据采集采用矢量数据结构，有利于保证空间实体的几何精度和拓朴特性的描述拓朴特性的描述空间分析则主要采用栅格数据结构，有利于加快系统的运行速度空间分析则主要采用栅格数据结构，有利于加快系统的运行速度和分析应用的进程和分析应用的进程因此，在地理信息系统中，需要进行数据结构的转换。因此，在地理信息系统中，需要进行数据结构的转换。地理信息系统原理GIS一、由矢量向栅格的转换一、由矢量向栅格的转换 矢量向栅格转换处理的根本任务：矢量向栅格转换处理的根本任务：通过就是把点、线或面的矢量数据，转换成对应的栅格数据。这一过程称为栅格化。

24、栅格化首先要建立矢量数据的平面坐标系和栅格行列坐标系之间的对应关系。0 0 1 1 2 2 3 3 4 4 5 5 6 6 7 71 12 23 34 45 56 67 78 89 98 8 9 9X XY Y地理信息系统原理GIS设矢量数据的一坐标点值为(x，y)，转成栅格数据其行列值为(I，J)。表示矢量数据的X最小值和Y最小值（一）点的转换地理信息系统原理GISMinX,MinY=1000,1000MinX,MinY=1000,1000X,Y=4600,2300X,Y=4600,2300XXYY0 1 2 3 4 5 6 7 8 90 1 2 3 4 5 6 7 8 9 9 9 8 8

25、7 7 6 6 5 5 4 4 3 3 2 2 1 1 0 0dx,dy=1000,1000dx,dy=1000,1000 9 9 8 8 7 7 6 6 5 5 4 4 3 3 2 2 1 1 0 00 1 2 3 4 5 6 7 8 90 1 2 3 4 5 6 7 8 90 1 2 3 4 5 6 7 8 90 1 2 3 4 5 6 7 8 9 9 9 8 8 7 7 6 6 5 5 4 4 3 3 2 2 1 1 0 0点要素栅格化的过程点要素栅格化的过程0 1 2 3 4 5 6 7 8 90 1 2 3 4 5 6 7 8 9地理信息系统原理GIS1）用点栅格化方法，实现直线的起

26、点和终点坐标点栅格化 用以上点栅格计算公式分别求出矢量数据中直线端点a、b的栅格行列值(I1、J1)和(I2、J2)。2）求出直线段所对应的栅格单元的行列值范围 这里直线段ab所对应的栅格单元的行范围为(I2-I1)；列范围为(J2-J1)。（二）线的转换地理信息系统原理GIS3）如果行列数差I、J|,分两种情况考虑：（1）列数大于行数 J I 地理信息系统原理GIS3）如果行列数差I、J|,分两种情况考虑：(2)行数大于列数I J地理信息系统原理GIS（三）面的栅格化1 1、基于弧段数据的栅格化方法、基于弧段数据的栅格化方法1 12 23 3两种不同的方法：基于弧段数据的栅格化和基于多边形数

27、据的栅格化。两种不同的方法：基于弧段数据的栅格化和基于多边形数据的栅格化。地理信息系统原理GIS交点交点 列列 左多边形左多边形 右多边形右多边形 1 0 0 11 0 0 1 1 12 23 3 2 3 1 22 3 1 2 3 6 2 3 3 6 2 3 4 9 3 0 4 9 3 0 地理信息系统原理GIS 主要针对实体结构的多边形矢量数据栅格化。就是在矢量表示的多边形边界内部的所有栅格上赋予相应的多边形属性值，从而形成栅格数据阵列。转换方法转换方法边界代数算法内点填充法包含检验法多边形数据的栅格化方法2 2、基于多边形数据的栅格化方法、基于多边形数据的栅格化方法地理信息系统原理GIS

28、边界代数法基于积分求多边形的思想，通过简单的代数运算，实现多边形的矢栅转换。该算法简单可靠，被大量使用。假定沿边界前进方向Y值下降为下行，Y值上升为上行。沿着多边形实体的边界环绕多边形一圈当向上环绕时，把边界左边一行中所有的栅格单元的数值都减去属性值。向下环绕时，把边界左边一行中所有的栅格单元的数值加属性值。(1)边界代数算法这样，多边形外部的栅格正负抵消，而多边形内部的栅这样，多边形外部的栅格正负抵消，而多边形内部的栅格被赋予属性值。格被赋予属性值。地理信息系统原理GIS地理信息系统原理GIS 首先按线的栅格化方法把多边形的边界栅格化，然后在多边形的内部找一个内点，从该点出发，向外填从多边形

29、区域，直到边界为止。(2)内点填充法内点填充法地理信息系统原理GIS判别方法：判别方法：由待判点对每个多边形的夹角和。如果夹角和为2，则该待判点属于此多边形，赋予多边形编号(纪录属性)；如果积分值为0，则该待判点在此多边形外部。1 1检验夹角法检验夹角法(3)包含检验法包含检验法地理信息系统原理GIS角度角度方向方向V0,V0,角度为逆时针方向；角度为逆时针方向；V0,V0,角度为顺时针方向角度为顺时针方向地理信息系统原理GIS (a)(a)夹角和为零夹角和为零 (b)(b)夹角和夹角和360360度度地理信息系统原理GIS由待判点向图外某点引射线，判断该射线与多边形所有边由待判点向图外某点引

30、射线，判断该射线与多边形所有边界相交的总次数。界相交的总次数。判别方法：判别方法：如相交偶数次，则待判点在该多边形的外部，如相交偶数次，则待判点在该多边形的外部，如相交奇数次，则待判点在该多边形内部。如相交奇数次，则待判点在该多边形内部。2 2 铅垂线法铅垂线法地理信息系统原理GIS二、由栅格向矢量的转换二、由栅格向矢量的转换目的目的一是将扫描仪获取的图像栅格数据存入矢量形式的空间数据库；一是将扫描仪获取的图像栅格数据存入矢量形式的空间数据库；二是将栅格数据进行数据压缩，将面状栅格数据转换为由矢量数据表二是将栅格数据进行数据压缩，将面状栅格数据转换为由矢量数据表示的多边形边界；示的多边形边界；

31、栅格数据结构向矢量数据结构的转换又称为矢量化。栅格数据结构向矢量数据结构的转换又称为矢量化。地理信息系统原理GIS边界栅格数据边界栅格数据填充栅格数据填充栅格数据方法方法：根据数据文件的不同，分别采用不同的算法：根据数据文件的不同，分别采用不同的算法：(一）图像数据的转换方法一）图像数据的转换方法 （二）栅格数据的转换方法（二）栅格数据的转换方法地理信息系统原理GIS 对对任任意意栅栅格格点点数数据据P，假假设设其其坐坐标标数数据据为为(I，J)，按按下下图图所所示示坐坐标标，将将其其转转换换为为矢矢量量数数据据，其其中中心心点点坐坐标标(x，y)计算公式为：计算公式为：点的矢量化点的矢量化

32、9 9 8 8 7 7 6 6 5 5 4 4 3 3 2 2 1 1 0 00 1 2 3 4 5 6 7 8 90 1 2 3 4 5 6 7 8 9地理信息系统原理GIS 图像数据是不同灰阶的影像，一般通过扫描仪按一定分图像数据是不同灰阶的影像，一般通过扫描仪按一定分辨率扫描采样，得到不同灰度值（辨率扫描采样，得到不同灰度值（0255）表示的数据。目）表示的数据。目前扫描仪的分辨率可达到前扫描仪的分辨率可达到0.0125mm，因此对一般粗度（例，因此对一般粗度（例如如0.1mm）的线条，其横断面扫描后平均也有）的线条，其横断面扫描后平均也有8个像元，而个像元，而矢量化的要求只能允许横断面

33、保持一个栅格的宽度，因此矢量化的要求只能允许横断面保持一个栅格的宽度，因此需要进行二值化、细化和跟踪等矢量化步骤。需要进行二值化、细化和跟踪等矢量化步骤。（一）基于图像数据的矢量化方法（一）基于图像数据的矢量化方法具体步骤如下图，分3步：二值化二值化 、细化、更踪、细化、更踪地理信息系统原理GIS（1）二值化）二值化 将用256级或128级灰度值G(i,j)量度的图像扫描数据的灰阶压缩到2个灰阶，即0和1两级。具体方法：在最大与最小灰阶之间定义一个阈值T，则根据下式得到二值图：地理信息系统原理GIS具体过程举例具体过程举例地理信息系统原理GIS（2）细化）细化细化细化是消除线划横断面栅格数的差

34、异，使得第一是消除线划横断面栅格数的差异，使得第一条线只保留代表其轴线的单个栅格的宽度。条线只保留代表其轴线的单个栅格的宽度。细化的方法可分为细化的方法可分为“剥皮法剥皮法”和和“骨架化骨架化”两大类。两大类。剥皮法剥皮法剥皮法的实质：剥皮法的实质：从曲线的边缘开始，每次剥掉等于一个栅格宽度的从曲线的边缘开始，每次剥掉等于一个栅格宽度的一层，直到最后留下彼此连通的、由单个栅格点组一层，直到最后留下彼此连通的、由单个栅格点组成的图形。成的图形。解决办法：根据待剥栅格为中心的解决办法：根据待剥栅格为中心的33栅格组合图栅格组合图来决定。来决定。地理信息系统原理GIS3X33X3栅格组合图栅格组合图

35、1 12 23 34 45 56 67 78 89 9101011111212131314141515161617171818191920202121222223232424252526262727282829293030313132323333343435353636373738383939404041414242434344444545464647474848494950505151地理信息系统原理GIS地理信息系统原理GIS地理信息系统原理GIS（3）跟踪）跟踪目的：将栅格数据转换为线段和闭合的线条，并以矢目的：将栅格数据转换为线段和闭合的线条，并以矢量形式存储线段的坐标。量形式存储线段

36、的坐标。方法：从起始点开始，根据八邻域搜索下一个相邻点，记方法：从起始点开始，根据八邻域搜索下一个相邻点，记录其坐标，直到完成全部栅格数据的矢量化。录其坐标，直到完成全部栅格数据的矢量化。地理信息系统原理GIS（二）栅格数据的矢量化方法（二）栅格数据的矢量化方法栅格数据的的矢量化常针对栅格数据中的多边栅格数据的的矢量化常针对栅格数据中的多边形进行。形进行。步骤：步骤：1 1、首先搜索多边形弧段相交的节点位置。、首先搜索多边形弧段相交的节点位置。2 2、任选一个起始点，顺着栅格单元属性值不同、任选一个起始点，顺着栅格单元属性值不同的两个栅格单元之间进行多边形边界弧段的跟的两个栅格单元之间进行多边

37、形边界弧段的跟踪矢量化。踪矢量化。3 3、将跟踪得到的弧段数据连接组织成多边形。、将跟踪得到的弧段数据连接组织成多边形。地理信息系统原理GIS第一步第一步搜索多边形弧段相交的节点位置搜索多边形弧段相交的节点位置第二步第二步跟踪弧段跟踪弧段第三步第三步构建多边形构建多边形地理信息系统原理GIS第三节 多源空间数据的融合 地理数据的多语义性、多时空性、多尺度性、获取手段的多样性、存储格式的不同以及数据模型与数据结构的差异等，导致了多源数据的产生，给数据的集成和信息共享带来了困难。因此，需要进行多源数据的融合。遥感影像DEM数据行政界线地理信息系统原理GIS一、遥感与GIS数据的融合遥感与遥感与GI

38、SGIS数据融合的具体方法表现为：数据融合的具体方法表现为：（1 1）遥感影像与数字线画图）遥感影像与数字线画图(DLG)(DLG)的融合。的融合。（2 2）遥感影像与数字地形模型（）遥感影像与数字地形模型（DEMDEM）的融合）的融合（3 3）遥感图像与数字栅格图（）遥感图像与数字栅格图（DRG)DRG)的融合的融合借助遥感技术获得的信息具有借助遥感技术获得的信息具有周期动态性、信息丰富、周期动态性、信息丰富、获取效率高等获取效率高等优点。优点。GISGIS具有具有空间数据管理和灵活的空空间数据管理和灵活的空间数据综合分析能力。间数据综合分析能力。两者的结合：有利于增强数据两者的结合：有利于

39、增强数据的复合能力，改善的复合能力，改善GISGIS数据的及时性和可靠性，便于利数据的及时性和可靠性，便于利用遥感影像辅助用遥感影像辅助GISGIS空间数据的获取与更新，提高数据空间数据的获取与更新，提高数据的使用率。的使用率。地理信息系统原理GIS遥遥感感影影像像与与数数字字线线画画图图(D DL LG G)的的融融合合地理信息系统原理GIS遥感影像与专题地图的复合遥感影像与专题地图的复合地理信息系统原理GIS遥感影像与DEM复合地理信息系统原理GISIRS-1CIRS-1C全色卫星影像（分辨率全色卫星影像（分辨率5.8m5.8m）同一地区的原始地形图（同一地区的原始地形图（DRGDRG）修

40、测后的地形图（修测后的地形图（1 1：5000050000）遥感影像与数字栅格图（遥感影像与数字栅格图（DRG)DRG)的融合的融合地理信息系统原理GIS二、不同格式数据的融合二、不同格式数据的融合目前目前GIS软件系统使用的空间数据格式主要有：软件系统使用的空间数据格式主要有：ESRI公司的公司的ARCINFOCoverage、ArcShapeFiles、E00格式；格式；Autodesk公司有公司有DXF和和DWG格式；格式；MapInfo公司的公司的TAB和和MIF格式格式Intergraph公司的公司的DGN格式格式不同格式数据融合的几种主要方法：不同格式数据融合的几种主要方法：（1

41、1）基于转换器的融合）基于转换器的融合（2 2）基于数据标准的数据融合）基于数据标准的数据融合（3 3）基于公共接口的数据融合）基于公共接口的数据融合（4 4）基于直接访问的数据融合）基于直接访问的数据融合地理信息系统原理GIS（1）基于转换器的数据融合）基于转换器的数据融合 在这种模式下，数据转换一般通过交换格式进行。在这种模式下，数据转换一般通过交换格式进行。地理信息系统原理GIS（2 2）基于数据标准的数据融合）基于数据标准的数据融合采用同一个标准来转换数据。采用一种空间数据的转换标准来实现多源GIS数据的融合。例如：美国国家空间数据协会（NSDI)制定了统一的空间数据格式规范SDTS(

42、Spatial Data Transformation Standard),包括几何坐标、投影、拓扑关系、属性数据、数据字典等。国土资源部在国家标准地球空间数据交换格式（简称CNSDTF）基础上制定了矢量数据交换格式VCT 2.0数据交换格式。地理信息系统原理GIS（3 3）基于公共接口的数据融合）基于公共接口的数据融合又称为数据互操作模式，接口就是一种标准。各个系统通过接口相互联系.OGC（Open GIS Consortium）为数据互操作制定了统一的规范。根据OGC的规范，可以将提供数据源的软件成为数据服务器（Data Servers)，使用数据的软件称为数据客户(Data Client

43、s)。地理信息系统原理GIS（4 4）基于直接访问的数据融合）基于直接访问的数据融合地理信息系统原理GIS（一）空间数据压缩与综合的意义（一）空间数据压缩与综合的意义（1）数据采集系统获得的坐标数据量极其巨大数据采集系统获得的坐标数据量极其巨大（2）建立无级比例尺数据库）建立无级比例尺数据库，随着比例尺的缩小，需要对要素进，随着比例尺的缩小，需要对要素进行简化。行简化。数据压缩的定义数据压缩的定义 所谓数据压缩，即从所取得的数据集合所谓数据压缩，即从所取得的数据集合S中抽出一个子集中抽出一个子集A，使这，使这个在规定的精度范围内最好地逼近原数据集合，而又取得尽可能大个在规定的精度范围内最好地逼

44、近原数据集合，而又取得尽可能大的压缩比。的压缩比。一、空间数据的压缩第四节 空间数据的压缩与重分类地理信息系统原理GIS曲线压缩的示意图（比例尺不断变小）曲线压缩的示意图（比例尺不断变小）地理信息系统原理GIS压缩比压缩比压缩比：压缩比：表示信息载体减少的程度。表示信息载体减少的程度。以一条线实体为例，设数据集以一条线实体为例，设数据集S中曲线的原来点序列为：中曲线的原来点序列为：A：A1，A2，An压缩处理后，获得新的子序列为：压缩处理后，获得新的子序列为：A：As1，As2，Asma值的大小，既与曲线的复杂程度、缩小倍数、精值的大小，既与曲线的复杂程度、缩小倍数、精度要求、数字化取点的密度

45、等因素有关，又与数据度要求、数字化取点的密度等因素有关，又与数据压缩技术有关。压缩技术有关。m1a=n压缩比为：压缩比为：地理信息系统原理GIS基于矢量的压缩通常是对线状实体上点的数基于矢量的压缩通常是对线状实体上点的数量的压缩，其中最常用的是：量的压缩，其中最常用的是：（二）基于矢量的压缩（二）基于矢量的压缩1 1、道格拉斯、道格拉斯-佩克算法佩克算法2 2、垂距法、垂距法 3 3、光栏法、光栏法地理信息系统原理GIS1 1、道格拉斯、道格拉斯佩克法佩克法(Douglas(DouglasPeucker)Peucker)基本思路与方法：基本思路与方法：基本思路与方法：基本思路与方法：1 1 1

46、 1）对对对对每每每每一一一一条条条条曲曲曲曲线线线线的的的的首首首首末末末末点点点点虚虚虚虚连连连连一一一一条条条条直直直直线线线线，求求求求所所所所有有有有点点点点与与与与直直直直线的距离，并找出最大距离值线的距离，并找出最大距离值线的距离，并找出最大距离值线的距离，并找出最大距离值dmaxdmaxdmaxdmax，用，用，用，用dmaxdmaxdmaxdmax与限差与限差与限差与限差D D D D相比；相比；相比；相比；dmax地理信息系统原理GIS2)2)若若dmaxdmaxD D，这条曲线上的中间点全部舍去；，这条曲线上的中间点全部舍去；地理信息系统原理GIS3)3)若若dmaxDd

47、maxD，保保留留dmaxdmax对对应应的的坐坐标标点点，并并以以该该点点为为界界，把把曲曲线线分分为为两两部部分分，对对这这两两部分重复使用该方法。部分重复使用该方法。道格拉斯道格拉斯佩克法佩克法(Douglas(DouglasPeucker)Peucker)地理信息系统原理GIS地理信息系统原理GIS2、垂距法垂距法 垂距法的基本思路是垂距法的基本思路是:1)1)每每次次顺顺序序取取曲曲线线上上的的三三个个点点，计计算算中中间间点点与与其其它它两两点点连连线线的的垂垂线线距距离离d d，并并与与限限差差D D比较。比较。d地理信息系统原理GIS2)2)2)2)若若若若d d d dD D

48、 D D，则中间点去掉，则中间点去掉，则中间点去掉，则中间点去掉,然后再依次取下三个点继续处理；然后再依次取下三个点继续处理；然后再依次取下三个点继续处理；然后再依次取下三个点继续处理；地理信息系统原理GIS3)3)3)3)若若若若dDdDdDdD，则中间点保留，则中间点保留，则中间点保留，则中间点保留,然后顺序取下三个点然后顺序取下三个点然后顺序取下三个点然后顺序取下三个点继续处理，直到这条线结束。继续处理，直到这条线结束。继续处理，直到这条线结束。继续处理，直到这条线结束。地理信息系统原理GIS3、光栏法光栏法 光光光光栏栏栏栏法法法法的的的的基基基基本本本本思思思思想想想想是是是是：定定

49、定定义义义义一一一一个个个个扇扇扇扇形形形形区区区区域域域域，通通通通过过过过判判判判断断断断曲曲曲曲线线线线上上上上的的的的点点点点在在在在扇扇扇扇形形形形外外外外还还还还是是是是在在在在扇扇扇扇形形形形内内内内，确确确确定定定定保保保保留留留留还还还还是是是是舍舍舍舍去去去去。设设设设曲曲曲曲线线线线上上上上的的的的点点点点列列列列为为为为p p p pi i i i，i i i i1 1 1 1，2 2 2 2，n n n n，光光光光栏栏栏栏口口口口径径径径为为为为d d d d，可可可可根根根根据据据据压压压压缩缩缩缩量量量量的的的的大小自己定义，则光栏法的实施步骤可描述为：大小自己

50、定义，则光栏法的实施步骤可描述为：大小自己定义，则光栏法的实施步骤可描述为：大小自己定义，则光栏法的实施步骤可描述为：地理信息系统原理GIS1 1 1 1）连连连连接接接接p p p p1 1 1 1和和和和p p p p2 2 2 2点点点点，过过过过p p p p2 2 2 2点点点点作作作作一一一一条条条条垂垂垂垂直直直直于于于于p p p p1 1 1 1p p p p2 2 2 2的的的的直直直直线线线线，在在在在该该该该垂垂垂垂线线线线上上上上取取取取两两两两点点点点a a a a1 1 1 1和和和和a a a a2 2 2 2，使使使使a a a a1 1 1 1p p p p