空间数据获取与处理.ppt

《空间数据获取与处理.ppt》由会员分享，可在线阅读，更多相关《空间数据获取与处理.ppt（46页珍藏版）》请在三一办公上搜索。

1、1,空间数据获取在GIS中的地位,GIS,硬件软件数据=127,2,第六章 空间数据获取与处理,1、地图数字化1.1 概述1.2 地图数据处理1.3 数字化仪数字化1.4 扫描矢量化及常用算法2、空间数据录入后的处理2.1 图形坐标变换2.2 图形拼接2.3 拓扑生成,3,1、地图数字化,1.1 概述地图是重要的信息形式纸质地图的弊端电子地图（GIS的基本功能）,当纸地图经过计算机图形图像系统光电转换量化为点阵数字图像，经图像处理和曲线矢量化，或者直接进行手扶跟踪数字化后，生成可以为地理信息系统显示、修改、标注、漫游、计算、管理和打印的矢量地图数据文件，这种与纸地图相对应的地图称为矢量化电子地

2、图。,4,矢量电子地图的优点：可视化：放大、缩小、漫游，分层显示地图；智能化：多尺度表达，自动配置注记，最短路径查询等空间分析功能；交互性：编辑，计算距离和标注地名等；时态性：监控，线路回放等。,5,电子地图多比例尺表达。比例尺越大，表达越详细。,6,航标监控中心,电子航道图与航标遥测监控系统,GPS定位技术,电子地图技术,通讯技术,传感技术,7,1.3 数字化仪数字化输入点、线及多边形地物的坐标工作量非常繁重,先确定需要数字化哪些信息，确定图层及图层包含的信息。确定三个定位点,8,两种方式：点方式流方式：数据量大，可以实时采样（化简）来减少数据量,9,曲线化简算法：Douglas-Peuck

3、er算法：,10,1.4 扫描矢量化及常用算法扫描矢量化及其处理过程,扫描仪,栅格转矢量的运算，一般称为扫描矢量化过程。在实际应用中，常常采用交互跟踪矢量化，或者称为半自动矢量化。,11,图像二值化平滑（去噪）细化链式编码矢量线提取,矢量化过程：,next,12,1.4 扫描矢量化及常用算法扫描矢量化及其处理过程,扫描仪,栅格转矢量的运算，一般称为扫描矢量化过程。在实际应用中，常常采用交互跟踪矢量化，或者称为半自动矢量化。,13,图像拼接：,14,图像裁剪：,15,图像细化预处理二值图像平滑：,去毛刺模板,去孔洞模板,16,图像细化：要求：1、保证细化后曲线的连通性2、细化结果是原曲线的中心线

4、3、保留细线端点,N(p)=2T(p)=1pW=0,pE=0,pS=1,pN=0,17,“剥皮”图像细化算法：算法思想：层层剥皮，从线条边缘开始一层一层向里剥夺,直到线条剩下一个象素的为止过程：1：对于栅格图像中的每个点p：如果2=N(p)=6并且T(p)=1并且pNpEpS=0并且pWpEpS=0则标志p点；2：将所有被标志的栅格点赋值为0，如果没有被标志的点，则算法结束；3：pNpSpW=0并且pWpEpN=04：5：转到第一步。,18,19,1 pNpEpS=0且 pWpEpS=023 pEpWpN=0且 pSpWpN=0450-第2步被被赋值0-第4步被被赋值,20,链码（Freema

5、n），矢量线生成（Douglas-Peucher）用曲线出发点坐标和线的斜率来描述二值线图形的一种方法。C(SPx,SPy)a1a2.an，0ai7,21,流程回顾,22,2、数据录入后的处理,2.1 图形坐标变换（几何纠正，投影变换）概念:利用一套控制点和变换方程,将数字地图或图像从一种坐标系转换成另一种坐标系的过程.2.1.1 几何纠正1)高次变换需要有6对以上控制点的坐标和理论值2)二次变换需要有5对以上控制点的坐标和理论值3)仿射变换 需要有3对以上控制点的坐标和理论值,23,3)仿射变换A 基本坐标变换：平移，缩放，旋转,平移：X=X+TxY=Y+Ty,24,缩放变换是指对p点相对于

6、坐标原点沿x方向放缩Sx倍，沿y方向放缩Sy倍。其中Sx和Sy称为比例系数。,Y,X,缩放变换(Sx=2,Sy=3),P(4,3),P(2,1),缩放：X=XSxY=YSy,25,二维旋转是指将p点绕坐标原点转动某个角度（顺时针为正，逆时针为负）得到新的点p的重定位过程。,旋转：X=Xcos+YsinY=-Xsin+Ycos,26,综合考虑图形的平移、旋转和缩放，则其坐标变换式如下：,仿射变换公式：,特点：1、直线变换后仍为直线；2、平行线变换后仍为平行线；3、不同方向上的长度比发生变化。,27,仿射变换的作用:将数据的空间坐标从数字化的坐标系统转换为真实世界中的坐标系统,28,均方根误差RM

7、S控制点的实际位置（由控制点经纬度数值投影到输出地图上的）与估算位置（输入地图上数字化的）之间偏差的估量一个控制点的RMS=sqrt(actX-estX)2+(actY-estY)2平均RMS=sqrt(在X和Y中偏差的平方和)/(控制点的数量)为确保几何变换的精确度，控制点的RMS误差必须控制在一定的容差值内。容差值依据源地图的精度和比例尺而异。,29,2.1.2 投影变换（1）解析变换法 找出两投影间坐标变换的解析计算公式 A 正解 B 反解,30,（2）数值变换法 如果原投影点的坐标解析式不知道，或不易求出两投影之间坐标的直接关系，可以采用多项式逼近的方法，即用数值变换法来建立两投影间的

8、变换关系式。例如，可采用二元三次多项式进行变换。二元三次多项式为：,31,（3）数值解析变换法当已知新投影的公式，但不知原投影的公式时，可先通过数值变换求出原投影点的地理坐标，然后代入新投影公式中，求出新投影点的坐标。即:,32,2.2 图幅拼接图幅尺寸标准分幅地形图一般先进行投影变换，从地形图使用的高斯克吕格投影转换到经纬度坐标系中，然后再进行拼接。,33,图幅自动拼接(1),拉框选择多个图幅,34,图幅自动拼接(2),35,2.3 拓扑生成(topology)2.3.1 图形修改2.3.2 建立拓扑关系,36,2.3.1 图形修改造成数字化错误的具体原因：1）遗漏某些实体；2）某些实体重复

9、录入；3）定位的不准确。,37,数字化后的地图上，错误的具体表现形式：1）伪节点（Pseudo Node）2）悬挂节点（Dangling Node）,38,数字化后的地图上，错误的具体表现形式：3）“碎屑”多边形或“条带”多边形（Sliver Polygon）4）不正规的多边形,39,DIME结构，双重独立地图编码（美国人口普查局于1970年提出）,a、多边形-组成弧段（多边形表）b、弧段-左右多边形，两端节点（弧段表）c、节点-相邻的弧段（节点表）,2.3.2 建立拓扑关系,40,几个概念：1、弧段有方向，弧段A的起始节点称为首节点NS(A)，而终止节点为尾节点NE(A)2、沿弧段前进方向，

10、将其相邻的多边形分别定义为左多边形和右多边形PL(A)和PR(A)3、N0 搜寻的起始节点，NC 当前节点,41,多边形拓扑建立过程,PL(A)=null？,得到一条弧段A，所有弧段处理完毕，算法结束,N,PR(A)=null？,N,创建多边形PPL(A)=PN0=NS(A)NC=NE(A)A为P的一条弧段,Y,创建多边形PPR(A)=PN0=NE(A)NC=NS(A)A为P的一条弧段,Y,N0=NC？,Y,检查当前节点连接的，当前弧段的下一条弧段 AA为P的一条弧段,NC=NS(A)？,PL(A)=P NC=NE(A),Y,PR(A)=P NC=NS(A),N,N,当前节点是否为当前弧段的起

11、点?,42,拓扑关系的重要性：1）维护空间数据一致性2）空间分析,道路交通应用中，一些道路虽然在平面上相交，但是实际上并不连通，如立交桥，这是需要手工修改，将连通的节点删除。,43,网络拓扑关系：在输入道路、水系、管网、通信线路等信息时，为了进行流量以及连通性分析，需要确定线实体之间的连接关系。,44,要求掌握,掌握：常用的两种地图数字化方法掌握：空间数据录入后要经过哪些处理?理解各处理过程的意义。掌握:曲线离散化的道格拉斯算法理解：仿射变换的原理和作用,45,思 考 题,1.栅格图像转换为矢量地图一般需要经过哪些步骤?2.请结合图形描述一下道格拉斯曲线离散化算法.,46,1)高次变换其中A、B代表二次以上高次项之和。上式是高次变换方程，符合上式的变换称为高次变换。在进行高次变换时，需要有6对以上控制点的坐标和理论值，才能求出待定系数。2)二次变换当不考虑高次变换方程中的A和B时，则变成二次变换方程，称为二次变换。二次变换适用于原图有非线性变形的情况，至少需要5对控制点的坐标及其理论值，才能求出待定系数。,了解,