第六章：空间数据获取.PPT

地理信息系统导论 多边形拓扑的建立，要注意多边形带“岛”的情况，按照上述算法，对于带“岛”的多边形，或者称为环，其包含的弧段构成了多个闭合曲线，并且“岛”的弧段排序是顺时针的（图6-20）（实际上，从环状多边形内部看，它仍然是逆时针的）。 图6-20：带“岛”的多边形建立拓扑的结果。 2．空间数据录入后的处理 2）网络拓扑关系的建立 在输入道路、水系、管网、通信线路等信息时，为了进行流量以及连通性分析，需要确定线实体之间的连接关系。网络拓扑关系的建立包括确定节点与连接线之间的关系。 图6-21：将实际不连通的线路形成的节点删除 2．空间数据录入后的处理 * 算法步骤如下： 1：对于栅格图像中的每个点p，进行如下操作： 如果2?N(p)?6并且T(p)=1并且pNpSpE=0并且pWpEpS=0 则标志p点； 2：将所有被标志的栅格点赋值为0，如果没有被标志的点，则算法结束； 3：对于栅格图像中的每个点p，进行如下操作： 如果2?N(p)?6并且T(p)=1并且pNpSpW=0并且pWpEpN=0 则标志p点； 4：将所有被标志的栅格点赋值为0，如果没有被标志的点，则算法结束； 5：转到第一步。 1．地图数字化 图6-9显示了采用该算法细化的过程和结果。 图6-9：线状地物的细化 1．地图数字化 1．4．5链码（弗里曼码） 链码是由弗里曼（Freeman）提出的用曲线出发点坐标和线的斜率来描述二值线图形的一种方法。图6-10（a）所示是链码的八个方向及它们的序号。图6-10（b）的细线的链码为(3,0)21100066567，其中（3,0）为起始点坐标，之后的数值序列描述了方向。 任意一条细线都可用链码序列表示为下式： C＝a1a2...an，0≤ai≤7 如果始点a1和终点an重合，则说明曲线是闭合的。 1．地图数字化 图6-10：链码及其对细线的表示 (a)：链码的8个方向；(b)细线 1．地图数字化 1．4．5矢量线生成 扫描矢量化的最后一步是生成矢量线，可以很方便地将链式编码的每一条链转换成为一条矢量线。自然地，弯曲的矢量线比直线需要更多的点，这还取决于要求的精度。在矢量线生成过程中，可以使用Douglas-Peucher算法。 1．地图数字化 2．空间数据录入后的处理 2．1图形坐标变换 在地图录入完毕后，经常需要进行投影变换，得到经纬度参照系下的地图。对各种投影进行坐标变换的原因主要是输入时地图是一种投影，而输出的地图产物是另外一种投影。进行投影变换有两种方式，一种是利用多项式拟合，类似于图像几何纠正；另一种是直接应用投影变换公式进行变换。 2．1．2 基本坐标变换 在投影变换过程中，有以下三种基本的操作：平移、旋转和缩放。 1）平移 平移是将图形的一部分或者整体移动到笛卡尔坐标系中另外的位置(图6-11-a)，其变换公式如下： X’=X+Tx Y’=Y+Ty (a)平移 2．空间数据录入后的处理 2）缩放 缩放操作可以用于输出大小不同的图形（图6-11-b），其公式为： X’=XSx Y’=YSy (b)缩放 2．空间数据录入后的处理 3）旋转 在地图投影变换中，经常要应用旋转操作（图6-11-c），实现旋转操作要用到三角函数，假定顺时针旋转角度为θ，其公式为*： X’=Xcosθ+Ysinθ Y’=-Xsinθ+Ycosθ (c)图形旋转 2．空间数据录入后的处理 2．1．2 仿射变换 综合考虑图形的平移、旋转和缩放，坐标变换式如下： 上式是一个正交变换，其更为一般的形式是： 2．空间数据录入后的处理 后者被称为二维的仿射变换，仿射变换在不同的方向可以有不同的压缩和扩张，可以将球变为椭球，将正方形变为平行四边形（图6-12）。 图6-12：仿射变换 2．空间数据录入后的处理 2．2图形拼接 在对底图进行数字化以后，由于图幅比较大或者使用小型数字化仪时，难以将研究区域的底图以整幅的形式来完成，这是需要将整个图幅划分成几部分分别输入。在所有部分都输入完毕并进行