CH4 空间数据处理.ppt

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

1、第四章第四章第四章第四章 空间数据处理空间数据处理空间数据处理空间数据处理靖娟利土木工程系主要内容4.1 4.1 坐标变换坐标变换vv坐标变换坐标变换坐标变换坐标变换的实质的实质的实质的实质是建立两个平面是建立两个平面是建立两个平面是建立两个平面点之间的一一对点之间的一一对点之间的一一对点之间的一一对应关系，包括应关系，包括应关系，包括应关系，包括几几几几何纠正和投影转何纠正和投影转何纠正和投影转何纠正和投影转换。换。换。换。vv如右图所示，旋如右图所示，旋如右图所示，旋如右图所示，旋转、平移、比例转、平移、比例转、平移、比例转、平移、比例缩放都属于坐标缩放都属于坐标缩放都属于坐标缩放都属于坐

2、标变换。变换。变换。变换。4.1 4.1 坐标变换坐标变换一、几何纠正一、几何纠正v为什么进行几何纠正？为什么进行几何纠正？地图变形地图变形(均匀变形、非均匀变形）均匀变形、非均匀变形）分幅扫描分幅扫描数字化中的位置移动数字化中的位置移动遥感影像本身存在几何变形遥感影像本身存在几何变形投影方式不同投影方式不同v目的目的为了实现对数字化数据的坐标系转换和图纸变形误差的改正为了实现对数字化数据的坐标系转换和图纸变形误差的改正v实质实质建立纠正图象与标准地图的一一对应关系建立纠正图象与标准地图的一一对应关系v变换方法变换方法精确方法：仿射变换精确方法：仿射变换(相似变换相似变换)、双线性变换、平方变

3、换、立方变换、双线性变换、平方变换、立方变换等等近似方法：近似方法：橡皮板变换橡皮板变换v纠正步骤纠正步骤纠正点纠正点数据采集数据采集函数建立函数建立逐点或网格纠正逐点或网格纠正橡皮板变换函数变换4.1 4.1 坐标变换坐标变换v1、相似变换：当图纸在图纸在x,y两个两个方向上的变形比例尺相同。方向上的变形比例尺相同。v2、仿射变换：设图纸变形引起设图纸变形引起x,y两个方向比例尺不同，此时两个方向比例尺不同，此时图纸存在仿射变形。图纸存在仿射变形。根据图形变换原理，得出坐根据图形变换原理，得出坐根据图形变换原理，得出坐根据图形变换原理，得出坐标公式标公式标公式标公式:p px xy yX X

4、Y YA A0 0B B0 0X X yyOO4.1 4.1 坐标变换坐标变换vv化简为化简为化简为化简为:vv误差方程为误差方程为误差方程为误差方程为:其中：其中：其中：其中：X X、Y Y为已知的理论坐标。为已知的理论坐标。为已知的理论坐标。为已知的理论坐标。vv由由由由QxQx2最小和最小和最小和最小和QyQy2最小的条件可得到两组法方程：最小的条件可得到两组法方程：最小的条件可得到两组法方程：最小的条件可得到两组法方程：其中其中其中其中n n为控制点个数，为控制点个数，为控制点个数，为控制点个数，x,yx,y为控制点的数字化坐标，为控制点的数字化坐标，为控制点的数字化坐标，为控制点的数

5、字化坐标，X X，Y Y为控制点的理论为控制点的理论为控制点的理论为控制点的理论值，值，值，值，a1a1，a2a2，a3a3，b1b1，b2b2，b3b3为待定系数。为待定系数。为待定系数。为待定系数。4.1 4.1 坐标变换坐标变换二、投影转换投影A（x，y）投影B（X，Y）正解变换：解析函数关系X=f(x,y)，Y=g(x,y)反解变换：经纬度B=f(x,y),L=g(x,y)X=F(B,L),Y=G(B,L)数值变换：数学方法4.1 4.1 坐标变换坐标变换v正解变换正解变换：通过建立一种投影变换与另一种投影变换的严密或近似的解析关系式，直接由一种投影的数字化坐标x，y变换为另一种投影的

6、直角坐标X，Y。v反解变换反解变换：由一种投影的坐标反解出地理坐标（x、y B、L），然后将地理坐标代入另一种投影的坐标公式中（B、Lx、y），实现由一种投影的坐标到另一种投影坐标的变换（x、yX、Y）。v数值变换数值变换：根据两种投影在变换区内的若干同名数字化点，采用插值法，或有限差分法，或有限元法，或待定系数法等，实现由一种投影的坐标到另一种投影坐标的变换。vv以上三种变换方法中，前两种是在已知原投影和新投以上三种变换方法中，前两种是在已知原投影和新投影解析式条件下，来求得两者之间的解析关系式；后影解析式条件下，来求得两者之间的解析关系式；后一种是不知原投影的解析式或不知原投影方程的情况一

7、种是不知原投影的解析式或不知原投影方程的情况下，采用多项式来确定两者之间的关系。下，采用多项式来确定两者之间的关系。4.2 4.2 图形编辑图形编辑v图形编辑图形编辑又叫数据编辑、数字化编辑，是又叫数据编辑、数字化编辑，是指对地图资料数字化后的数据进指对地图资料数字化后的数据进行编辑加工；行编辑加工；主要的目的是在改正数据差错的主要的目的是在改正数据差错的同时，相应地改正数字化资料的同时，相应地改正数字化资料的图形。图形。4.2 4.2 图形编辑图形编辑一、图形数字化常见错误一、图形数字化常见错误v伪节点（伪节点（Pseudo Node）伪节点使一条完整的线变成两段伪节点使一条完整的线变成两段

8、 v 悬挂节点（悬挂节点（Dangling Node）如果一个节点只与一条线相连接，那么该节点称为悬如果一个节点只与一条线相连接，那么该节点称为悬挂节点挂节点伪节点伪节点悬挂节点悬挂节点4.2 4.2 图形编辑图形编辑v“碎屑碎屑”多边形或多边形或“条带条带”多边形（多边形（Sliver Polygon）一般由于重复录入引起，由于前后两次录入同一条线的一般由于重复录入引起，由于前后两次录入同一条线的位置不可能完全一致，造成了位置不可能完全一致，造成了“碎屑碎屑”多边形。多边形。v不正规的多边形（不正规的多边形（Weird Polygon）是由于输入线时，点的次序倒置或者位置不准确引起的是由于输

9、入线时，点的次序倒置或者位置不准确引起的4.2 4.2 图形编辑图形编辑二、结点的编辑结点的编辑v1 1、结点吻合、结点吻合(Snap)(Snap)，或称结点匹配、和结点附合。，或称结点匹配、和结点附合。结点结点移动移动鼠标鼠标拉框拉框求交点求交点编辑编辑方法方法自动匹自动匹配配4.2 4.2 图形编辑图形编辑v2、结点与线的吻合、结点与线的吻合在数字化过程中，常遇到一个结点与一个线状目标的中间相交。在数字化过程中，常遇到一个结点与一个线状目标的中间相交。由于测量或数字化误差，它不可能完全交于线目标上，需要进行由于测量或数字化误差，它不可能完全交于线目标上，需要进行编辑，称为结点与线的吻合。编

10、辑，称为结点与线的吻合。编辑编辑方法方法结点结点移动移动线段线段求交求交自动自动编辑编辑A AB BD DC CE E无结点无结点4.2 4.2 图形编辑图形编辑v3、清除假结点由仅有两个线目标相关联的结点成为假结点。由仅有两个线目标相关联的结点成为假结点。有些系统要将这种假结点清除掉（如有些系统要将这种假结点清除掉（如ARC/INFO），即），即将目标将目标A 和和B合并成一条，使它们之间不存在结点合并成一条，使它们之间不存在结点;但有些系统并不要求清除假结点，如但有些系统并不要求清除假结点，如Geostar,因为它因为它们并不影响空间查询、分析和制图。们并不影响空间查询、分析和制图。A A

11、B B4.2 4.2 图形编辑图形编辑v三、图形编辑包括用鼠标增加或删除一个点、线、面实体，移动、旋转一个点、包括用鼠标增加或删除一个点、线、面实体，移动、旋转一个点、线、面实体。线、面实体。边界弧段顶边界弧段顶点的编辑点的编辑删除与增加删除与增加一个顶点一个顶点删除一条删除一条弧段弧段移动一个移动一个顶点顶点j jk kj jk ka ab bL3L3L1L1L2L2j jk k4.2 4.2 图形编辑图形编辑v四、数据检查与清理四、数据检查与清理数据检查指拓扑关系的检查，结点是否匹配，是否存数据检查指拓扑关系的检查，结点是否匹配，是否存在悬挂弧段，多边形是否封闭，是否有假结点。在悬挂弧段，

12、多边形是否封闭，是否有假结点。要求系统能将有错误或不正确的拓扑关系的点、线和要求系统能将有错误或不正确的拓扑关系的点、线和面用不同的颜色和符号表示出来，以便于人工检查和面用不同的颜色和符号表示出来，以便于人工检查和修改。修改。数据清理则是用自动的方法清除空间数据的错误数据清理则是用自动的方法清除空间数据的错误.例如给定一个结点吻合的容差使该容差范围内的结点自例如给定一个结点吻合的容差使该容差范围内的结点自动吻合在一起，并建立拓扑关系。给定悬挂弧段容差，动吻合在一起，并建立拓扑关系。给定悬挂弧段容差，将小于该容差的短弧自动删除。将小于该容差的短弧自动删除。4.2 4.2 图形编辑图形编辑v五、关

13、键算法点的捕捉点的捕捉123线的捕捉线的捕捉多边形的捕捉多边形的捕捉4.2 4.2 图形编辑图形编辑v1、点的捕捉点的捕捉设光标点为设光标点为S(x,yS(x,y)，某一点状，某一点状要素的坐标为要素的坐标为A(XA(X，Y)Y)；可设一捕捉半径可设一捕捉半径D D。若。若S S和和A A的的距离距离d d小于小于D D则认为捕捉成功，则认为捕捉成功，即认为找到的点是即认为找到的点是A A，否则失，否则失败，继续搜索其它点；败，继续搜索其它点；乘方运算影响了搜索的速度，因乘方运算影响了搜索的速度，因此，把距离此，把距离d d的计算改为：的计算改为：捕捉范围由圆改为矩形，这可捕捉范围由圆改为矩形

14、，这可大大加快搜索速度。大大加快搜索速度。4.2 4.2 图形编辑图形编辑v2 2、线的捕捉、线的捕捉设光标点坐标为设光标点坐标为S(x,yS(x,y)，D D为捕捉半径，为捕捉半径，线的坐标为线的坐标为(x1,y1),(x2,y2),(x1,y1),(x2,y2),(xn,ynxn,yn)。通过计。通过计算算S S到该线的每个直线段的距离到该线的每个直线段的距离d d。若若min(d1,d2,min(d1,d2,dn-1)dn-1)D D，则认为光标，则认为光标S S捕捉到了该条线，否则为未捕捉到。捕捉到了该条线，否则为未捕捉到。加快线捕捉的速度的方法：加快线捕捉的速度的方法：每计算一个距离

15、每计算一个距离didi就进行一次比较，若就进行一次比较，若didiD D，则捕捉成功，不需再进行下面，则捕捉成功，不需再进行下面直线段到点直线段到点S S的距离计算了；的距离计算了；把不可能被光标捕捉到的线，用简单算把不可能被光标捕捉到的线，用简单算法去除。法去除。对于线段也采用类似的方法处理。对于线段也采用类似的方法处理。简化距离公式：点简化距离公式：点S(x,yS(x,y)到直线段到直线段(x1,y1),(x2,y2)(x1,y1),(x2,y2)的距离的距离d d的计算公式为：的计算公式为：点到线的距离点到线的距离简化简化方法方法简化距离公式简化距离公式4.2 4.2 图形编辑图形编辑v

16、3 3、面的捕捉、面的捕捉实际上就是判断光标点实际上就是判断光标点S(x,yS(x,y)是否在是否在多边形内，若在多边形内则说明捕捉多边形内，若在多边形内则说明捕捉到。到。判断点是否在多边形内的算法主要有判断点是否在多边形内的算法主要有垂线法或转角法垂线法或转角法。垂线法的基本思想是从光标点引垂线垂线法的基本思想是从光标点引垂线(实际上可以是任意方向的射线实际上可以是任意方向的射线)，计，计算与多边形的交点个数。算与多边形的交点个数。若交点个数为奇数则说明该点在多边若交点个数为奇数则说明该点在多边形内；若交点个数为偶数，则该点在形内；若交点个数为偶数，则该点在多边形外。多边形外。加快速度的方法

17、：加快速度的方法：找出该多边形的外接矩形，若光标点找出该多边形的外接矩形，若光标点落在该矩形中，才有可能捕捉到该面，落在该矩形中，才有可能捕捉到该面，否则放弃对该多边形的进一步计算和否则放弃对该多边形的进一步计算和判断。判断。对不可能有交点的线段应通过简单的对不可能有交点的线段应通过简单的坐标比较迅速去除。坐标比较迅速去除。运用计算交点的技巧。运用计算交点的技巧。4.3 4.3 拓扑关系的建立拓扑关系的建立点线点线(弧弧)拓扑关系的建立拓扑关系的建立12在图形采集和编辑中实时建立在图形采集和编辑中实时建立多边形拓扑关系自动建立多边形拓扑关系自动建立在图形采集和编辑之后自动建立在图形采集和编辑之

18、后自动建立4.3 4.3 拓扑关系的建立拓扑关系的建立v一、点线点线(弧弧)拓扑关系的建立拓扑关系的建立在图形采集和编辑中实时建立在图形采集和编辑中实时建立结点结点-弧段表弧段表Oid起结点终结点a1a2N1N2N2N3Oid弧段号N1N2N3a1a1,a2a2弧段弧段-结点表结点表Oid起结点终结点a1a2a3N1N2N2N2N3N4Oid弧段号N1N2N3N4a1a1,a2,a3a2a3Oid起结点终结点a1a2a3a4N1N2N2N4N2N3N4N3Oid弧段号N1N2N3N4a1a1,a2,a3a2,a4a3,a4a1a1a2a2 N1 N1 N2N2N3N3a1a1a2a2 N1 N

19、1 N2N2N4N4a3a3a1a1a2a2 N1 N1 N2N2N3N3N4N4a3a3a4a44.3 4.3 拓扑关系的建立拓扑关系的建立v二、多边形拓扑关系自动建立多边形拓扑关系自动建立链的组织链的组织结点匹配结点匹配拓扑检查拓扑检查 找出在链的中间找出在链的中间相交的情况，自相交的情况，自动切成新链；动切成新链；把链按一定顺序把链按一定顺序存储，并把链按存储，并把链按顺序编号。顺序编号。把一定限差内把一定限差内的链的端点作为的链的端点作为一个结点，其坐一个结点，其坐标值取多个端点标值取多个端点的平均值。的平均值。对结点顺序编对结点顺序编号。号。通过判断一条通过判断一条链的端点是否链的端

20、点是否有与之匹配的有与之匹配的端点来进行端点来进行.创建多边形创建多边形在所有的错误在所有的错误都改正的情况都改正的情况下，建立多边下，建立多边形。形。4.4 4.4 图形的裁剪、合并图形的裁剪、合并v一、图形裁剪（开窗处理开窗处理 ）方式：方式：正窗，提取窗口内的数据。开负窗，提取窗口外正窗，提取窗口内的数据。开负窗，提取窗口外的数据子集。的数据子集。矩形窗和多边形窗。矩形窗和多边形窗。算法：算法：点、线、面的窗口裁剪点、线、面的窗口裁剪-计算机图形学。计算机图形学。不规则多边形开窗不规则多边形开窗-相当于多边形叠置处理。相当于多边形叠置处理。4.4 4.4 图形的裁剪、合并图形的裁剪、合并

21、v二、图形合并（图形合并（数据文件合并数据文件合并）一幅图内的多层数据合并在一起一幅图内的多层数据合并在一起;或将相邻的多幅图的同一层数据或将相邻的多幅图的同一层数据合并合并.涉及到空间拓扑关系的重建。对于多边形，由于同一个目标在两涉及到空间拓扑关系的重建。对于多边形，由于同一个目标在两幅图内已形成独立的多边形，合并时，需幅图内已形成独立的多边形，合并时，需去除公共边界去除公共边界，属性合，属性合并，具体算法，删去共同线段。并，具体算法，删去共同线段。实际处理过程是先删除两个多边实际处理过程是先删除两个多边形，解除空间关系后，删除公共边，再重建拓扑。形，解除空间关系后，删除公共边，再重建拓扑。

22、p pL1L1p p A AA Ap pL1L1p p A AA Ap pL1L1p p A AA A去除去除公共边界公共边界属性属性合并合并4.5 4.5 图幅接边图幅接边v一、为什么进行接边处理？几何裂缝几何裂缝：指由数据文件边界分开的一个地物的两部：指由数据文件边界分开的一个地物的两部分不能精确地衔接。分不能精确地衔接。逻辑裂缝逻辑裂缝：同一地物地物编码不同或具有不同的属性：同一地物地物编码不同或具有不同的属性信息，如公路的宽度，等高线高程等。信息，如公路的宽度，等高线高程等。v二、图幅接边二、图幅接边识别或提取相邻图幅识别或提取相邻图幅（当图幅编号合理的情况下）（当图幅编号合理的情况下

23、）几何接边几何接边逻辑接边逻辑接边4.5 4.5 图幅接边图幅接边v几何接边人人工工接接边边接边接边直接移动，突变直接移动，突变回缩回缩2-32-3个点个点减少突变减少突变 4.5 4.5 图幅接边图幅接边v逻辑接边目的目的：检查同一地物在相邻图幅检查同一地物在相邻图幅的地物编码和属性值是否的地物编码和属性值是否一致，不一致，进行人工一致，不一致，进行人工编辑。编辑。将同一地物在相邻图幅的将同一地物在相邻图幅的空间数据在逻辑上连在一空间数据在逻辑上连在一起。起。方法方法：索引文件，建立双向指针。索引文件，建立双向指针。关键字，空间操作的方法。关键字，空间操作的方法。Oid指针A1AOid指针A

24、2AOid指针A3AOid指针AA1A2A3Oid KeyA1AOid KeyA2AOid KeyA3A索引文件索引文件关键字关键字A3A3A AA1A1A2A24.6 4.6 数据压缩与光滑数据压缩与光滑v一、数据压缩一、数据压缩一、数据压缩一、数据压缩（矢量数据压缩）（矢量数据压缩）（矢量数据压缩）（矢量数据压缩）目的目的目的目的：减少数据存储量：减少数据存储量：减少数据存储量：减少数据存储量 前提前提前提前提：位置的正确性，力求压缩后的曲线与压缩前的曲线偏差位置的正确性，力求压缩后的曲线与压缩前的曲线偏差位置的正确性，力求压缩后的曲线与压缩前的曲线偏差位置的正确性，力求压缩后的曲线与压缩

25、前的曲线偏差最小；最小；最小；最小；弯曲特征的正确性，力求保留反映曲线形态特征的特征弯曲特征的正确性，力求保留反映曲线形态特征的特征弯曲特征的正确性，力求保留反映曲线形态特征的特征弯曲特征的正确性，力求保留反映曲线形态特征的特征点。点。点。点。方法方法方法方法：距离限差（垂距法、距离限差（垂距法、距离限差（垂距法、距离限差（垂距法、DouglasDouglasPeuckerPeucker）角度限差角度限差角度限差角度限差4.6 4.6 数据压缩与光滑数据压缩与光滑v垂距法基本思路：基本思路：每次顺序取曲线上的每次顺序取曲线上的三个三个点点，计算中间点与其它两，计算中间点与其它两点连线的垂线距离

26、点连线的垂线距离d，并，并与限差与限差D比较。若比较。若dD，则中间点去掉；若则中间点去掉；若dD，则中间点保留。然后顺序则中间点保留。然后顺序取下三个点继续处理，直取下三个点继续处理，直到这条线结束。到这条线结束。特点特点：压缩算法好，可在数压缩算法好，可在数字化时实时处理，每次判断字化时实时处理，每次判断下一个数字化的点，且计算下一个数字化的点，且计算量较小。量较小。4.6 4.6 数据压缩与光滑数据压缩与光滑vDouglasPeucker 基本思路基本思路基本思路基本思路：对每一条曲线的首末点虚连对每一条曲线的首末点虚连一条直线，求所有点与直线一条直线，求所有点与直线的距离，并找出最大距

27、离值的距离，并找出最大距离值dmax，用，用dmax与限差与限差D相相比比;若若dmaxD，这条曲线上的，这条曲线上的中间点全部舍去；中间点全部舍去；若若dmaxD，保留，保留dmax对应对应的坐标点，并以该点为界，的坐标点，并以该点为界，把曲线分为两部分，对这两把曲线分为两部分，对这两部分重复使用该方法。部分重复使用该方法。特点特点特点特点：压缩效果好，但必须在：压缩效果好，但必须在对整条曲线数字化完成后才能进对整条曲线数字化完成后才能进行，且计算量较大。行，且计算量较大。4.6 4.6 数据压缩与光滑数据压缩与光滑v二、曲线光滑二、曲线光滑二、曲线光滑二、曲线光滑（拟合）（拟合）（拟合）（

28、拟合）实质实质实质实质：是假想曲线为一是假想曲线为一组离散点，寻找形式较组离散点，寻找形式较简单、性能良好的曲线简单、性能良好的曲线解析式。解析式。曲线光滑处理，曲线光滑处理，也称采样点曲线拟合。也称采样点曲线拟合。方法方法：张力样条函数、：张力样条函数、分段三次多项式、分段分段三次多项式、分段圆弧和切线方向内插等。圆弧和切线方向内插等。4.6 4.6 空间数据格式转换空间数据格式转换v矢量向栅格的转换矢量向栅格的转换 栅格化栅格化点的栅格化点的栅格化点的栅格化点的栅格化123线的栅格化线的栅格化线的栅格化线的栅格化面的栅格化面的栅格化面的栅格化面的栅格化4.6 4.6 空间数据格式转换空间数

29、据格式转换yx010109 98 87 76 65 54 43 32 21 12 23 34 45 56 67 78 89 91010dxdxdydyx xminmin,y,yminminI IJ J坐标系统坐标系统4.6 4.6 空间数据格式转换空间数据格式转换v一、点的栅格化一、点的栅格化 设：矢量坐标点（设：矢量坐标点（设：矢量坐标点（设：矢量坐标点（x x，y y）转换后的单元行列值为转换后的单元行列值为转换后的单元行列值为转换后的单元行列值为I I，J J 栅格单元大小：栅格单元大小：栅格单元大小：栅格单元大小：dxdx，dydy 图形最小坐标值：图形最小坐标值：图形最小坐标值：图形

30、最小坐标值：x xminmin，y yminmin 则有：则有：则有：则有：4.6 4.6 空间数据格式转换空间数据格式转换v二、线的栅格化二、线的栅格化 设：线段两端点坐标分别为：设：线段两端点坐标分别为：设：线段两端点坐标分别为：设：线段两端点坐标分别为：（x1x1，y1y1）和（和（和（和（x2x2，y2y2），转换后的单元行列值分别为：），转换后的单元行列值分别为：），转换后的单元行列值分别为：），转换后的单元行列值分别为：(I1I1,J1J1),(),(I2I2,J2 J2)行数差：行数差：行数差：行数差：I=|I2 I1|I=|I2 I1|，列数差：列数差：列数差：列数差：J=|J

31、2 J1|J=|J2 J1|分两种情况：分两种情况：分两种情况：分两种情况：列数差大于行数差列数差大于行数差列数差大于行数差列数差大于行数差 J J II 行数差大于列数差行数差大于列数差行数差大于列数差行数差大于列数差 I I JJ4.6 4.6 空间数据格式转换空间数据格式转换v三、面的栅格化三、面的栅格化 又称为边界填充，就是在矢量表示的多边形界内部的又称为边界填充，就是在矢量表示的多边形界内部的所有栅格上赋予相应的多边形编号。所有栅格上赋予相应的多边形编号。基于弧段数据的栅格化基于弧段数据的栅格化基于弧段数据的栅格化基于弧段数据的栅格化基于多边形数据的栅格化基于多边形数据的栅格化基于多

32、边形数据的栅格化基于多边形数据的栅格化方法：方法：方法：方法：4.6 4.6 空间数据格式转换空间数据格式转换vv基于弧段数据的栅格化基于弧段数据的栅格化基于弧段数据的栅格化基于弧段数据的栅格化（带扫描算法）（带扫描算法）（带扫描算法）（带扫描算法）对每条弧段求与各扫描线交点对每条弧段求与各扫描线交点记录每个交点的左右区码记录每个交点的左右区码逐段生成栅格逐段生成栅格111111122211111222221111322222113333222233333332223333333322弧段弧段弧段弧段1 1弧段弧段弧段弧段2 2弧段弧段弧段弧段3 311322(xi,yi)(xi+1,yi+1

33、)(x2,y2)(x3,y3)(x1,y1)y ymm4.6 4.6 空间数据格式转换空间数据格式转换vv基于多边形数据的栅格化基于多边形数据的栅格化基于多边形数据的栅格化基于多边形数据的栅格化-边界代数法（边界代数法（边界代数法（边界代数法（BAFBAF）基本思想基本思想基本思想基本思想:对每幅地图的全部具有左右多边形编号的边界弧段，沿对每幅地图的全部具有左右多边形编号的边界弧段，沿其前进的方向逐条搜索；其前进的方向逐条搜索；当弧段上行时，将该弧段与左图框之间的栅格加上一个当弧段上行时，将该弧段与左图框之间的栅格加上一个值值（该弧段的左多边形编号）（该弧段的右多边形（该弧段的左多边形编号）（

34、该弧段的右多边形编号）；编号）；当弧段下行时，将该弧段与左图框之间的栅格加上一个当弧段下行时，将该弧段与左图框之间的栅格加上一个值值（该弧段的右多边形编号）（该弧段的左多边形（该弧段的右多边形编号）（该弧段的左多边形编号）；编号）；当弧段平行与栅格行行走时，不作运算。当弧段平行与栅格行行走时，不作运算。4.6 4.6 空间数据格式转换空间数据格式转换4.6 4.6 空间数据格式转换空间数据格式转换v栅格向矢量的转换栅格向矢量的转换-矢量化矢量化矢量化矢量化点的矢量化点的矢量化点的矢量化点的矢量化123线的矢量化线的矢量化线的矢量化线的矢量化面的矢量化面的矢量化面的矢量化面的矢量化4.6 4.6

35、 空间数据格式转换空间数据格式转换v一、点的矢量化点的矢量化 设：单元行列值为设：单元行列值为设：单元行列值为设：单元行列值为I I，J J 转换后的中心坐标点（转换后的中心坐标点（转换后的中心坐标点（转换后的中心坐标点（x x，y y）栅格单元大小：栅格单元大小：栅格单元大小：栅格单元大小：dxdx，dydy 图形最小坐标值：图形最小坐标值：图形最小坐标值：图形最小坐标值：XminXmin，Y Yminmin 则有：则有：则有：则有：4.6 4.6 空间数据格式转换空间数据格式转换v二、线的矢量化线的矢量化 基于图像数据的矢量化方法基于图像数据的矢量化方法基于图像数据的矢量化方法基于图像数据

36、的矢量化方法 对扫描仪获得的地图图像中线段中心线进行矢量跟踪对扫描仪获得的地图图像中线段中心线进行矢量跟踪对扫描仪获得的地图图像中线段中心线进行矢量跟踪对扫描仪获得的地图图像中线段中心线进行矢量跟踪 手工跟踪屏幕数字化手工跟踪屏幕数字化手工跟踪屏幕数字化手工跟踪屏幕数字化 自动跟踪矢量化自动跟踪矢量化自动跟踪矢量化自动跟踪矢量化4.6 4.6 空间数据格式转换空间数据格式转换v三、面的矢量化面的矢量化 提取以相同编号的栅格集合表示的多边形区域的边界和边界的拓提取以相同编号的栅格集合表示的多边形区域的边界和边界的拓扑关系，并表示成多个小直线断的矢量格式边界线的过程。扑关系，并表示成多个小直线断的

37、矢量格式边界线的过程。步骤：步骤：多边形边界提取多边形边界提取多边形边界提取多边形边界提取：图像二值化；：图像二值化；边界线追踪边界线追踪边界线追踪边界线追踪：对每个边界弧段由一个结点向另一个结点搜索，：对每个边界弧段由一个结点向另一个结点搜索，通常对每个已知边界点需沿除了进入方向的其他通常对每个已知边界点需沿除了进入方向的其他7 7个方向搜索个方向搜索下一个边界点，直到连成边界弧段；下一个边界点，直到连成边界弧段；拓扑关系生成拓扑关系生成拓扑关系生成拓扑关系生成：对于矢量表示的边界弧段数据，判断其与原图：对于矢量表示的边界弧段数据，判断其与原图上各多边形的空间关系，以形成完整的拓扑结构并建立

38、与属性上各多边形的空间关系，以形成完整的拓扑结构并建立与属性数据的联系；数据的联系；去除多余点及曲线光滑去除多余点及曲线光滑去除多余点及曲线光滑去除多余点及曲线光滑：去除多余点，进行曲线光滑。：去除多余点，进行曲线光滑。4.6 4.6 空间数据格式转换空间数据格式转换v基于图像处理的矢量化方法基于图像处理的矢量化方法步骤：步骤：二值化二值化细化细化跟踪跟踪创建拓扑创建拓扑4.6 4.6 空间数据格式转换空间数据格式转换二值化：是在一个设定的灰度阈值的基础上，对扫描获得的二值化：是在一个设定的灰度阈值的基础上，对扫描获得的灰度图像进行灰度图像进行0或或1的简化处理；的简化处理；（a）扫描前的矢量

39、数据 （b）扫描得到的灰度值 （c）二值化后得到的栅格数据4.6 4.6 空间数据格式转换空间数据格式转换细化：细化：将二值图像象元阵列逐步剥除轮廓边缘的点，使之成为将二值图像象元阵列逐步剥除轮廓边缘的点，使之成为线划宽度只有一个象元的骨架图形。线划宽度只有一个象元的骨架图形。细化后的图形骨架既保留了原图形的绝大部分特征，又细化后的图形骨架既保留了原图形的绝大部分特征，又便于下一步的跟踪处理。便于下一步的跟踪处理。细化的代表性方法：剥皮法、骨架法细化的代表性方法：剥皮法、骨架法扫描地图细化处理应符合下列基本要求：扫描地图细化处理应符合下列基本要求：保持原线划的连续性；线宽只为一个象元；细划后的

40、骨架应是原线划的中心线；保持图形的原有特征。4.6 4.6 空间数据格式转换空间数据格式转换跟综：把骨架转换为矢量图形的坐标序列。跟综：把骨架转换为矢量图形的坐标序列。基本步骤：从左向右，从上向下搜索线划起始点，并记下坐标；朝该点的8个方向追踪点，若没有，则本条线的追踪结束，转1进行下条线的追踪；否则记下坐标。把搜索点移到新取的点上，转2。创建拓扑创建拓扑 4.6 4.6 空间数据格式转换空间数据格式转换v基于窗口匹配的矢量化算法基于窗口匹配的矢量化算法（双边界直接搜索算法，（双边界直接搜索算法，DBDFDBDF）基本思想：基本思想：通过边界提取，将左右多边形信息保存在边界点上，每条边界弧段由

41、两通过边界提取，将左右多边形信息保存在边界点上，每条边界弧段由两个并行的边界链组成，分别记录该边界弧段的左右多边形编号个并行的边界链组成，分别记录该边界弧段的左右多边形编号。转化步骤：转化步骤：边界点和节点提取：边界点和节点提取：采用采用2*22*2栅格阵列作为窗口顺序沿行、列方向对栅栅格阵列作为窗口顺序沿行、列方向对栅格图像全图扫描，如果窗口内四个栅格有且仅有两个不同的编号，则该格图像全图扫描，如果窗口内四个栅格有且仅有两个不同的编号，则该四个栅格表示为边界点；如果窗口内四个栅格有三个以上不同编号，则四个栅格表示为边界点；如果窗口内四个栅格有三个以上不同编号，则标识为结点，保持各栅格原多边形

42、编号信息。对于对角线上栅格两两相标识为结点，保持各栅格原多边形编号信息。对于对角线上栅格两两相同的情，由于造成了多边形的不连通，也当作结点处理。同的情，由于造成了多边形的不连通，也当作结点处理。节点的节点的8 8种结构种结构 边界点的边界点的6 6种结构种结构 4.6 4.6 空间数据格式转换空间数据格式转换边界线搜索与左右多边形信息记录：边界线搜索与左右多边形信息记录：首先记录开始边界点的两个首先记录开始边界点的两个多边形编号，作为该弧段的左右多边形，下一点组的搜索方向则多边形编号，作为该弧段的左右多边形，下一点组的搜索方向则由进入当前点的搜索方向和该点组的可能走向决定，由进入当前点的搜索方

43、向和该点组的可能走向决定，每个边界点每个边界点组只能有两个走向，一个是前点组进入的方向，另一个则可确定组只能有两个走向，一个是前点组进入的方向，另一个则可确定为将要搜索后续点组的方向。为将要搜索后续点组的方向。小结小结v本章主要讨论了空间数据处理的基本方法。包括坐标变换、图形编辑、创建拓扑、数据压缩、图幅接边、数据格式转换等。v重点掌握坐标变换的原理、图形编辑的基本方法、矢量数据压缩方法、图幅接边的基本原理、数据转换的主要方法。习题习题v1如何发现进入GIS中的数据有错误？v2一般从扫描仪上直接得到的地图有什么问题？如何改正？v3如果两个作业小组各自从数字化仪上得到两张相邻图幅的地图数据在GIS中不能准确对接该怎么办？v4空间数据几何纠正的常用方法有哪些？v5假设一条矢量等高线上的点太过于密集了，如何减少占用系统的存储空间？你能给出多少方法？各有什么适用范围？v6.为什么要进行图幅拼接？v7.空间数据格式转换的主要方法？