第五讲空间数据处理

（优选）第五讲空间数据处理目前一页\总数五十二页\编于十八点基本算法设在平面笛卡尔坐标系中的两点a（x1，y1）和b（x2，y2）两点的距离两点的方向线不考虑线的方向考虑线的方向，怎么办？目前二页\总数五十二页\编于十八点基本算法点到线目标的距离设有一直线段l，两端点的坐标为（xA，yA）和（xB，yB）,另一给定点P的坐标为（xP，yP）。根据解析几何可知直线方程为：

ax+by+c=0其中a=yB-yAb=xA-xBc=yAxB-xAyB

此种方法求得的为点到直线的垂直距离，垂足可能在线的延长线上，这不是我们所需要的？如何解决？目前三页\总数五十二页\编于十八点基本算法线的长度或者面的周长线目标的长度和面目标的周长，实际上是所有线段两点之间距离的总和

面的面积是所有边与X或者Y轴构成的直角梯形的代数和

目前四页\总数五十二页\编于十八点基本算法多边形的重心重心是分别求多边形顶点x坐标和y坐标的平均值重心不一定是在多边形内，有时重心在多边形外

目前五页\总数五十二页\编于十八点基本算法线相交（交点）设有两条线段AB和CD，它们的端点坐标分别为（xA，yA），（xB，yB），（xC，yC），（xD，yD）

求出的交点可能不是在两条线之间，而y在它们的延长线上也可能两条线平行，而没有交点如何解决？目前六页\总数五十二页\编于十八点基本算法曲线光滑处理张力样条分段三次多项式分段圆弧切线方向内插法五点法目前七页\总数五十二页\编于十八点基本算法曲线化简（压缩）

Douglas—Peucker压缩效果好，但必须在对整条曲线数字化完成后才能进行，且计算量较大目前八页\总数五十二页\编于十八点基本算法曲线化简（压缩）垂距法每次顺序取曲线上的三个点，计算中间点与其它两点连线的垂线距离d，并与限差D比较。若d＜D，则中间点去掉；若d≥D，则中间点保留。然后顺序取下三个点继续处理，直到这条线结束。目前九页\总数五十二页\编于十八点基本算法曲线化简（压缩）光栏法定义一个扇形区域，通过判断曲线上的点在扇形外还是在扇形内，确定保留还是舍去。目前十页\总数五十二页\编于十八点基本算法平行线处理直线段的平行线若线段AB的斜率为k=tgα其中α是AB的方向角，则距AB为d的平行线与过A、B两点法线的交点为A′、B′，则：其中（xa，ya），（xb，yb），（xa′，ya′），（xb′，yb′）分别是A、B、A′和B′的坐标。

目前十一页\总数五十二页\编于十八点基本算法点在多边形中的判断点在多边形内的判别最直接的方法是铅重线法或者说平行线法或者说射线法，即从需判别的点开始划一任一方向的直线，（该直线可以是铅直线或平行线），然后计算它所通过多边形的交点，当交点的个数是奇数时，该点在多边形内，若是偶数，表明它在多边形外

目前十二页\总数五十二页\编于十八点基本算法点在多边形中的判断但是使用射线法有时候可能失效，产生判断错误。当射线通过多边形的拐点或某一条边时，这时按统计通过多边形边界交点的奇偶数，产生错误的判断结果目前十三页\总数五十二页\编于十八点基本算法线与多边形求交线与多边形是否相交，需要判断每条线段与多边形的边界线段是否有交点如果没有任何交点，再判断端点是在多边形内还是多边形外，如果两端点在多边形外，线段又与多边形不相交，则该线段相离多边形，如果两点都在多边形内，并且与多边形边界没有交点，则该线段在多边形内如果有一个或多个交点，该线段与多边形相交，部分在多边形内，部分在多边形外即使两个点都在多边形内如GH或都在多边形外如IJ，它们都可能与多边形相交所以判断线与多边形是否相交，仅判别端点是不够的，必须判断线状目标的每一段与多边形边界的每一段是否有交点

目前十四页\总数五十二页\编于十八点基本算法多边形与多边形相交判断两个多边形是否相交需要判断两个多边形边界的所有线段相互之间是否有交点。如果没有任何交点，它们可能相分离，也可以一个多边形在另一个多边形之内两个多边形边界线段只要存在一个交点则表明两个多边形相交如果它们公共一条边界，则它们相邻目前十五页\总数五十二页\编于十八点基本算法区域填充种子法目前十六页\总数五十二页\编于十八点基本算法区域填充扫描线法8目前十七页\总数五十二页\编于十八点图形编辑图形编辑又叫数据编辑、数字化编辑，是指对地图资料数字化后的数据进行编辑加工其主要的目的是在改正数据差错的同时，相应地改正数字化资料的图形图形编辑是一交互处理过程，GIS具备的图形编辑功能的要求是：1）具有友好的人机界面，即操作灵活、易于理解、响应迅速等；2）具有对几何数据和属性编码的修改功能，如点、线、面的增加、删除、修改等；3）具有分层显示和窗口操作功能，便于用户的使用。目前十八页\总数五十二页\编于十八点图形编辑点的选择设光标点为S(x,y),某一点状要素的坐标为A(X，Y)可设一捕捉半径D(通常为3～5个象素，这主要由屏幕的分辩率和屏幕的尺寸决定)若S和A的距离d小于D则认为捕捉成功，即认为找到的点是A，否则失败，继续搜索其它点。乘方运算影响了搜索的速度，因此，把距离d的计算改为：捕捉范围由圆改为矩形，这可大大加快搜索速度。

目前十九页\总数五十二页\编于十八点图形编辑线的选择设光标点坐标为S(x,y)，D为捕捉半径，线的坐标为(x1,y1),(x2,y2),…(xn,yn)。通过计算S到该线的每个直线段的距离d。.

若min(d1,d2,…dn-1)＜D，则认为光标S捕捉到了该条线，否则为未捕捉到。目前二十页\总数五十二页\编于十八点图形编辑面的选择实际上就是判断光标点S(x,y)是否在多边形内，若在多边形内则说明捕捉到。垂线法垂线法的基本思想是从光标点引垂线(实际上可以是任意方向的射线)，计算与多边形的交点个数。若交点个数为奇数则说明该点在多边形内；若交点个数为偶数，则该点在多边形外。加快速度的方法：找出该多边形的外接矩形，若光标点落在该矩形中，才有可能捕捉到该面，否则放弃对该多边形的进一步计算和判断。目前二十一页\总数五十二页\编于十八点图形编辑结点咬合－结点匹配结点移动用鼠标拉一个矩形，落入在这矩形内的结点坐标符合成一致，即求它们的中点坐标通过求交点的方法自动匹配目前二十二页\总数五十二页\编于十八点图形编辑结点与线的咬合在数字化过程中，经常遇到一个结点与一个线状目标的中间相交，这时由于测量误差，它也可能不完全交于线目标上，而需要进行编辑，称为结点与线的咬合结点移动线段求交自动编辑的方法目前二十三页\总数五十二页\编于十八点图形编辑伪结点的删除由仅有两个线目标相关联的结点称为假结点有些系统要将这种假结点清除掉，有些系统并不要求清除假结点这些所谓的假结点并不影响空间查询、空间分析和制图方法：线段或者弧断合并需要考虑两种情况A、要求坐标一致，而不建立拓扑关系；如高架桥（不需打断，直接移动）B、不仅坐标一致，且要建立之间的空间关联关系；如道路交叉口（需要打断）目前二十四页\总数五十二页\编于十八点图形编辑伪结点的删除由仅有两个线目标相关联的结点称为假结点有些系统要将这种假结点清除掉，有些系统并不要求清除假结点这些所谓的假结点并不影响空间查询、空间分析和制图方法：线段或者弧断合并目前二十五页\总数五十二页\编于十八点图形编辑删除与增加结点目前二十六页\总数五十二页\编于十八点图形编辑移动一个结点目前二十七页\总数五十二页\编于十八点图形编辑删除一条弧断目前二十八页\总数五十二页\编于十八点图形编辑数据检查与清理目前二十九页\总数五十二页\编于十八点图形编辑RedoandUndo对于一个实用的地理信息系统软件，Undo与Redo功能是必要的然而这一功能的实现对于空间数据的编辑来说是相当困难的因为它不象字处理系统那样，在进行文本编辑时只要记录编辑操作和编辑的字符即可进行Undo与Redo这里除了要记录编辑的操作和编辑目标的空间坐标以外，还要记录它的目标标识，以及空间拓扑关系目前三十页\总数五十二页\编于十八点属性编辑类似于关系数据库的编辑目前三十一页\总数五十二页\编于十八点图形的裁剪矩形裁剪排除没有交点的方法目前三十二页\总数五十二页\编于十八点图形的裁剪矩形裁剪求交目前三十三页\总数五十二页\编于十八点图形的裁剪多边形裁剪目前三十四页\总数五十二页\编于十八点图形的裁剪多边形裁剪逐边裁剪法目前三十五页\总数五十二页\编于十八点图形的合并线线合并方法面面合并方法P1P2

P目前三十六页\总数五十二页\编于十八点图幅接边由于空间数据采集的误差和人工操作的误差，两个相邻图幅的地图的空间数据在结合处可能出现逻辑裂隙与几何裂隙。逻辑裂隙指的是当一个地物在一幅图的数据文件中具有地物编码A，而在另一幅图的数据文件中却具有地物编码B，或者同一个物体在这两个数据文件中具有不同的属性信息，如公路的宽度，等高线的高程等。几何裂隙指的是由数据文件边界分开的一个地物的两部分不能精确地衔接。在地理信息系统和机助制图中，需要把单独数字化的相邻图幅的空间数据在逻辑上和几何上融成一个连续一致的数据体图幅接边包括几何接边和逻辑接边目前三十七页\总数五十二页\编于十八点图幅接边几何接边调出需要接边的两幅或多幅图数据，以其中的一个作为活动图幅（或称活动工作区），其它图幅的作为参考，沿图幅的边缘选取一定范围例如5厘米的空间目标。这些目标（主要是弧段）一般都终结于图廓边附近，以活动工作区的目标为基准，根据图廓边上弧段的结点坐标查找相邻图幅对应弧段，如果它们的地物编码相同，结点坐标在一定的容差范围内，则将两边的结点坐标取中数自动吻合，空间关系不变。如果地物编码不同，或超过接边的匹配容差，则需要进行人工编辑与接边

超限需人工接边目前三十八页\总数五十二页\编于十八点图幅接边逻辑接边逻辑接边包括两方面的含义，一是检查同一目标在相邻图幅的地物编码和属性赋值是否一致，如果不一致，则进行人工编辑修改。这种逻辑接边容易处理。另一种逻辑接边的含义是将同一目标在相邻图幅的空间实体数据在逻辑上连在一起。例如长江可能跨越多个图幅，当要进行查询时，点取到某幅图的一段目标时要能够同时将多幅图内的长江一起显示出来，这就要在逻辑上建立某种联系。否则，由于每幅图的数据是单独存贮（见第六章），一般来说只能查询到该图幅内的空间数据（全关系型无缝数据库除外）。目前三十九页\总数五十二页\编于十八点图幅接边逻辑接边为了进行空间目标的逻辑接边，可以有两种方案，一种是在图幅数据文件的上一层，将有逻辑联系的空间目标，建立一个新的文件，即索引到它在每幅图的子目标，并建立双向指针（目标标识）。当在某一幅图点取子目标时，通过指针，指向上一层总目标文件的记录，这一条记录记录了所有该目标的子目标的目标标识，通过它即可显示整个目标

目前四十页\总数五十二页\编于十八点图幅接边逻辑接边逻辑接边的另外一种方法是不建立总目标文件，也不在每幅图的空间目标的数据文件中为逻辑接边的子目标建立索引，而是通过空间操作的方法，根据每个关键字如“长江”，让系统自动在周边图幅的文件中搜索到同一目标，从而在效果上，等同于建立了跨图幅空间目标的逻辑关系。

目前四十一页\总数五十二页\编于十八点坐标变换产生原因：坐标系的不统一、数据变形、地图投影的差异等。实质是建立两个平面点之间的一一对应关系。包括几何纠正和投影转换。目前四十二页\总数五十二页\编于十八点坐标变换相似变换设XOY为新的平面直角坐标系如地面大地坐标系，xoy为旧的平面直角坐标系如数字化仪坐标系，两坐标系之间的坐标轴夹角为，0′相对于XOY坐标系原点的平移距离为A0，B0，两坐标系之间坐标的比例因子为m，则根据坐标变换原理，可写出变换公式为：X=m（xcos-ysin）+A0Y=m（xsin+ycos）+B0

目前四十三页\总数五十二页\编于十八点坐标变换仿射变换

如果坐标在X、Y方向的比例因子不一致，或者说图纸存在仿射变形，此时需要采用仿射变换公式。令m1和m2分别表示X和Y方向的比例尺，则变换公式为：X=（m1cosα）x-（m1sinα）y+A0Y=（m2sinα）x+（m2cosα）y+B0目前四十四页\总数五十二页\编于十八点其他高次变换：X=A0+A1x+A2y+A11x2+A12xy+A22y2+AY=B0+B1x+B2x+B11x2+B12xy+B22y2+B二次变换X=A0+A1x+A2y+A11x2+A12xy+A22y2Y=B0+B1x+B2x+B11x2+B12xy+B22y2坐标变换目前四十五页\总数五十二页\编于十八点当系统所使用的数据是来自不同地图投影的图幅时，必须按需求统一地图投影，需要进行地图投影变换。基本原理：X=f1(x,y)Y=f2(x,y)投影变换目前四十六页\总数五十二页\编于十八点基本方法：解析变换法：找出两投影间的解析计算公式正解变化法:又称直接变换法,是通过两种投影之间的严密数学解析关系，直接由一种投影的数字化坐标精确变化到令一种投影的坐标反解变换法：又称间接变化法，以地理坐标（φ，λ）为中间媒介，先将一种投影坐标反解出地理坐标，然后再将其地理坐标代入另一种投影的坐标计算公式中，计算出新的投影坐标x,yφ,λX,Y解析