实验二 地图配准文档格式.docx
- 文档编号:20263620
- 上传时间:2023-01-21
- 格式:DOCX
- 页数:31
- 大小:7.11MB
实验二 地图配准文档格式.docx
《实验二 地图配准文档格式.docx》由会员分享,可在线阅读,更多相关《实验二 地图配准文档格式.docx(31页珍藏版)》请在冰豆网上搜索。
10万
5万
2.5万
1万
5000
图幅
范围
经
差
6°
3°
1°
30′
15′
7′30″
3′45″
1′52.5″
纬
4°
2°
20′
10′
5′
2′30″
1′15″
行列
关系
行
数
1
2
4
12
24
48
96
192
列
图
幅
间
量
关
系
16
144
576
2304
9216
36864
36
9
64
256
比例尺代码
A
B
C
D
E
F
G
H
2.2地图投影
地图配准是为了使得影像数据可以和GIS矢量数据集成在一起,而为影像数据指定一个参考坐标系的过程,因此在学习如何进行地图配准之前,本节有必要对我国常用的地图投影及SuperMapGIS6R中坐标系类型、投影设置等内容进行介绍。
2.2.1我国常用地图投影
1、高斯—克吕格投影
(1)基本概念
如图3.1所示,假想有一个椭圆柱面横套在地球椭球体外面,并与某一条子午线(此子午线称为中央子午线或轴子午线)相切,椭圆柱的中心轴通过椭球体中心,然后用一定投影方法,将中央子午线两侧各一定经差范围内的地区投影到椭圆柱面上,再将此柱面展开即成为投影面,如图3.2所示,此投影为高斯投影。
高斯投影是正形投影的一种。
图3.1高斯投影示意
图3.2中央子午线附近投影示意
(2)分带投影
由于高斯—克吕格投影中变形随着离中央经线的距离增大而递增,为了使投影的变形不致过大,所以采用分带投影。
我国的1:
2.5万~1:
50万地形图采用6°
分带;
1:
1万及更大比例尺图采用3°
分带。
分带是从零子午线起,每隔6°
为一带,全球共分为60个6°
带。
分带是从1°
30′起,每隔3°
为一带,全球共分为120个3°
高斯投影6°
带:
自0°
子午线起每隔经差6°
自西向东分带,依次编号1,2,3….。
我国6°
带中央子午线的经度,由75°
起每隔6°
至135°
,共计11带(13~23带),带号用n表示,中央子午线的经度用L0表示,它们的关系式L0=6n-3。
高斯投影3°
它的中央子午线一部分同6°
带中央子午线重合,一部分同6°
带的分界子线重合,如用nˊ表示3°
带的带号,L表示3°
带中央子午线经度,它们的关系L=3nˊ,我国3°
带共计22带(24~45带)。
如图3.3所示为6°
带与3°
带投影示意图。
图3.36°
带投影示意图
2、兰勃特(Lambert)投影
我国1:
100万地形图采用兰勃特(Lambert)投影,其分幅原则与国际地理学会规定的全球统一使用的国际百万分之一地图投影保持一致。
纬度按纬差4°
分带,从南到北共分15个投影带,每隔投影带单独计算坐标,每带两条标准纬线,第一标准纬线为图幅南端纬度加30′的纬线,第二标准纬线为图幅北端纬度减30′的纬线,这样处于同一投影带中的各图幅的坐标成果完全相同,不同带的图幅变形值接近相等,因此每投影带只需计算其中一幅图(纬差4°
,经差6°
)的投影成果即可。
由于是纬差4°
分带投影的,所以当沿着纬线方向拼接地图时,不论多少图幅,均不会产生裂隙;
但是,当沿着经线方向拼接时,因拼接线分别处于上下不同的投影带,投影后的曲率不同,致使拼接时产生裂隙。
2.2.2SuperMap坐标系类型
在SuperMapGIS6R中,数据的坐标系分为三类,平面坐标系、地理坐标系、投影坐标系。
1、平面坐标系
一般用来作为与地理位置无关的数据的坐标参考,也是默认新建数据的坐标参考,如CAD设计图、纸质地图扫描后的图片、与地理位置无关的示意图等。
平面坐标系是一个二维坐标系,原点坐标(0,0),数据中每一个点的坐标是由其距水平和垂直的X轴和Y轴的距离确定。
一般城市规划中用到的大比例尺地图,如1:
500,1:
1000等的道路施工图、建筑设计图等多采用平面坐标系。
2、地理坐标系
地理坐标系,也可称真实世界的坐标系,是用于确定地物在地球上位置的坐标系,其通过经纬度坐标来表示椭球上任意一点的坐标。
我国三大常用地理坐标系有WGS-84、北京54、西安80等。
3、投影坐标系
投影坐标系使用二维平面坐标(X,Y)来表示地物位置。
投影坐标系由地理坐标系和投影方法来确定。
一般经过投影的地理数据,可进行地图量算、各种空间分析、制图表达等。
2.2.3SuperMap投影设置
图3.4投影设置下拉按钮
1、“开始”选项卡→“数据”组→“投影设置”按钮用于设置与管理当前工作空间中数据源或数据集的坐标投影信息。
注:
该按钮在工作空间管理其中有被选中的数据源时可以使用。
2、单击“投影设置”下拉按钮,下拉菜单中包含“常用投影”和“投影设置”两项,常用投影显示收藏夹内的投影文件。
若收藏夹内无任何投影文件,“常用投影”项无内容显示。
图3.5常用投影及投影设置
3、单击“投影设置”项,打开“投影设置”窗口,可在该窗口中设置当前选中数据源或数据集的投影信息。
图3.6投影设置窗口
功能区域:
主要用于进行投影配置文件的浏览、设置与管理等功能。
投影信息管理目录树:
主要用于管理投影信息文件的分组以及浏览投影信息文件,其作用和呈现形式类似于windows的资源管理器。
文件列表区域:
主要用于浏览、设置于管理投影配置文件以及收藏夹中包含的投影信息文件。
在投影信息管理目录树中选中文件夹结点不同,文件列表区域显示的内容不同。
描述信息区域:
主要用于显示投影信息文件的描述信息。
4、平面坐标系设置
(1)在投影设置窗口的投影信息管理目录树中,选择“平面坐标系”。
(2)在投影设置窗口的文件列表区域,通过双击修改平面坐标系的单位(默认为米)。
(3)在弹出的对话框中选择“是”。
(4)点击“应用”按钮。
至此完成平面坐标系单位的修改。
图3.7平面坐标系设置
5、设置地理坐标系
(1)系统预定义地理坐标系
在投影设置窗口的投影信息管理目录树中,在“地理坐标系”节点下选择“Default”文件夹,接下来在右侧的文件列表中选择需要的坐标系,如“Beijing54”,接下来点击“应用”按钮完成预定义地理坐标系的设置。
图3.8设置预定义地理坐标系
(2)自定义地理坐标系
在投影设置窗口的投影信息管理目录树中,在“地理坐标系”节点下选择“Default”文件夹,接下来在右侧的文件列表中单击鼠标右键,在弹出的下拉框中选择“自定义坐标系”。
在弹出的自定义地理坐标系对话框中完成自定义地理坐标系。
1)设置相关参数
若选择系统提供的地理坐标系,大地参照系和中央经线等参数是固定不可以编辑的;
若用户选择自定义(UserDefined),大地参照系的类型下拉列表框则被激活,为可编辑状态。
大地参照系的类型下拉列表框、椭球参数的类型下拉列表框和中央经线下拉列表框都是类似的,选择了UserDefined后,可以编辑其相关参数。
●大地参考系:
基于地球椭球体建立的,确定了地球椭球体相对于地球球心的位置,为地表地物的测量提供了一个参照框架,确定了地表经纬网线的原点和方向,简言之,就是椭球体的定位和定向。
每个国家或地区均有各自的基准面,我国常用的两个大地参考系包括:
北京54坐标系(Beijing1954)、西安80坐标系(Xian1980)。
用户也可根据实际需要,设置大地类型为自定义(UserDefined),即可自定义椭球参数和中央经线的相关参数。
●椭球参数:
地球椭球体表面是一个用数据公式表达的规则的数学表面,在测量和制图中常用地球椭球体表面用来代替地球的自然表面。
椭球参数"
"
在应用时应该根据各个国家或地区的具体情况选择合适的地球椭球体。
我国常用的椭球体包括:
Krasovsky1940(赤道半径为6378245.0米,扁率为1/298.3)、International1975(赤道半径为6378140,扁率为1/298.257)等。
用户也可根据实际需要,设置椭球类型为自定义(UserDefined),即可自定义赤道半径和扁率的具体数值。
●中央经线:
定义地理坐标系横坐标的起算位置。
在大部分坐标系统中,中央经线是指经过英国伦敦格林威治的经线(Greenwich,经度为0)。
用户也可根据实际需要,设置中央经线类型为自定义(UserDefined),即可自定义其经度的具体数值。
图3.9自定义地理坐标系
6、投影坐标系设置
(1)系统预定义投影坐标系
在投影设置窗口的投影信息管理目录树中,在“投影坐标系”节点下选择所需的投影方法如Gauss-kruger(BEIJING1954),接下来在右侧的文件列表中选择分带投影,如“GKZone19N”,接下来点击“应用”按钮完成预定义投影坐标系的设置。
图3.10预定义投影坐标系设置
(2)系统自定义投影坐标系
在右侧的文件列表中单击鼠标右键,在弹出的下拉框中选择“自定义坐标系”。
在弹出的自定义地理坐标系对话框中完成自定义投影坐标系。
图3.11自定义投影坐标系设置
2.3地图要素
地图要素是构成地图的基本内容,分为数学要素、地理要素和辅助要素。
数学要素包括地图投影、制图网、比例尺、大地控制基础以及方位标等。
地理要素是指地图内容,包括自然地理要素与社会经济要素,据其性质可分为:
自然要素,包括海洋、陆地水系以及地质、地球物理、地貌、气象、水文、土质与植被、动物等;
社会经济要素,包括居民地、交通网、境界以及政治、行政、人口、城市、历史、文化和经济等方面的现象或物体;
环境要素,是指人类生活的环境状况,包括自然灾害、自然保护、环境污染及其保护与治理、疾病与医疗等。
辅助要素是有利于读图、用图方面的内容,包括:
图名、图号、接图表、接合图号、图廓、图廓间注记、图例、数字比例尺与图解比例尺、坡度尺、三北方向图以及其他附图、资料和成图方法说明等。
本节用本书示例数据中的SuperMap影像配准及数据采集中的“地图数据采集.udb”中数据源中J48G024008栅格影像来完成,通过本节使读者能够准确把握扫描地形图配准时所需要的地理信息。
2.3.1数学要素
获取地形图投影信息是能够正确进行地图配准的基础。
如图3.12所示为新荣县1:
1万标准图幅地形图。
根据地形图正下方的比例尺可以知道该地形图为1:
1万地形图;
根据地形图左下角的文字可以看出,该地形图在制作时采用的是西安80地理坐标系;
根据左上角的经纬度及格网坐标值可以看出,该地形图采用的投影方法为高斯-克吕格。
根据3.1节中的内容我国1:
1万地形图均采用高斯-克吕格3度带投影,因此可以判读该地形图投影信息为:
地理坐标采用的是西安80,投影方法采用高斯-克吕格3度带,带号为34带,中央经线为102度。
图3.12数学要素
2.3.2地理要素
正确判读地形图地理要素对于地图数据采集具有十分重要的意义。
如图3.13所示为该地形图上部分地理要素示意图。
从地形图上可以看出,该地形图包括河流、居民地、等高线、交通等地理要素。
图3.13地理要素
2.3.3辅助要素
如图3.14所示,该地形图辅助要素包括图例、接图表、图名、图幅号及其他附图、资料和成图方法说明等。
图3.14辅助要素
2.4地图配准
本节用本书示例数据中的SuperMap影像配准及数据采集中的“地图数据采集.udb”中数据源中J48G024008栅格影像、交通线两个图层来完成,通过本节使读者能够掌握地图配准的基本原理及SuperMapGIS6R中地图配准的操作方法。
2.4.1配准基本原理
1、基本概念
(1)数据配准:
数据配准是通过参考数据集(图层)对配准数据集(图层)进行空间位置纠正和变换的过程。
通过确定的配准算法和控制点信息,对配准数据集进行配准,可以得到与参考数据集(图层)空间位置一致的配准结果数据集。
(2)配准图层:
是指没有空间信息的图层,待配准图层。
(3)参考图层:
是指配准过程中,通过在配准图层中选择合适的配准点,同时在参考图层的相应位置上也需要选择控制点。
配准过程会将配准图层中配准点位置通过一定的配准算法,转换到与参考图层一致的空间位置。
2、配准算法简介
(1)线性配准
也称仿射变换。
这种配准方法假设地图因变形而引起的实际比例尺在X和Y方向上不相同,因此,具有纠正地图变形的功能。
假设x,y为输入图像的原始坐标,X,Y为纠正后的输出坐标,则坐标变换的公式如下:
上式中有6个参数:
A、B、C、D、E、F,可以通过选择配准图层和参考图层的同名控制点来确定。
理论上要求上述6个待定参数,至少需要知道不在同一直线上的3对控制点坐标,才能实现线性变换。
如果取常数项参数C、F为0,则至少需要4个点,所以实际应用中,通常利用4个或者4个以上的控制点来进行线性配准。
待定系数是通过最小二乘法的原理进行求解的。
线性变换是最常用的一种配准方法,由于同时考虑了X和Y方向上的变形,所以纠正后的坐标在不同的方向上的长度比会不同,表现为原始坐标会发生如缩放、旋转、平移等变化后得到输出坐标。
(2)二次多项式配准
是常用的精度较高的配准方法。
多项式纠正把原始图像变形看成是某种曲面,输出图像为规则平面。
从理论上讲,任何曲面都能以适当高次的多项式来拟合。
一般用二次二项式就可以对变形比较严重的图像进行纠正,并得到较高的精度。
二次二项式的纠正公式如下:
上式中共有12个参数,理论上需要6对控制点就可实现转换,但为了得到比较高的精度,一般要求二次多项式纠正的控制点为至少7对,适当增加控制点的个数,可以明显提高影像配准的精度。
多项式系数是用所选定的控制点坐标,按照最小二乘法求得的。
对比多项式纠正公式与线性纠正的公式,可以看出线性配准其实是一次的多项式配准,二次多项式纠正比线性配准需要更多的控制点对,也会纠正更多的变形。
(3)矩形配准
实质上是一种特殊的,有限定条件的线性配准。
如果原图像为规则矩形,纠正后的图像坐标仍是规则矩形,则选择两个相对的角点就可以确定矩形4个角点的坐标,在线性配准的纠正公式中的常数项取作0,则已知2个相对角点坐标就可以求解得到纠正公式,从而实现几何纠正。
这种方法既方便省时,也避免了由于选择多个控制点时造成的误差累积。
矩形配准是一种简单方便的配准纠正方法,但是因为输出结果不会计算误差,所以其配准的精度不可知,是一种精度不高的粗纠正方法。
(4)偏移配准
偏移配准仅需要一组控制点和参考点,分别对X坐标和Y坐标求差值,再利用差值对原数据集所有组坐标点进行偏移。
2.4.2无参考图层地图配准
1、打开“地图数据采集与编辑.udb”数据源,该数据源下存在新荣县1:
1万标准地形图数据“J48G024008”,该地形图数据没有投影信息,双击“J48G024008”,在地图窗口中打开该地形图。
图3.15地图窗口
2、新建配准。
单击功能区→“数据”选项卡→“配准”组→“新建配准”按钮,弹出“配准数据设置”对话框。
图3.16配准数据设置对话框
3、在配准图层处选择数据集为“J48G024008”;
参考图层数据集为空;
勾选“另存数据集”,在数据源处选择“地图数据采集与编辑”,校正后数据集名称为“J48G0024008_adjust”。
图3.17配准数据设置
4、单击“确定”按钮,弹出配准窗口。
图3.18配准窗口
5、选择配准算法。
单击功能区→“配准”选项卡→“算法”组→在配准算法下拉框中选择常用的二次多项式配准算法,如图3.19所示。
图3.19选择配准算法
6、设置标注色和高亮色。
单击功能区→“配准”选项卡→“控制点设置”组→单击“标注色”下拉按钮,选择合适的颜色。
这里使用默认红色。
单击“高亮色”下拉按钮,在弹出的颜色面板中,选择合适的颜色。
这里使用默认蓝色。
7、刺点
(1)单击功能区→“配准”选项卡→“控制点设置”组→单击“刺点”按钮。
图3.20刺点工具
(2)当配准窗口中的鼠标状态变为
时,找到公里格网的交点,单击鼠标左键完成刺点,此时控制点用蓝色十字丝标注。
图3.21刺点
(3)双击目标点X下方的蓝色区域,弹出如右图的输入控制点对话框。
图3.22双击目标点X下方处
图3.23输入控制点对话框
(4)在输入控制点对话框目标点X坐标处,根据公里格网坐标输入控制点实际X坐标538000。
(5)在输入控制点对话框目标点Y坐标处,根据公里格网坐标输入控制点实际Y坐标4323000。
(6)单击确定按钮,此时控制点为红色十字丝标注。
图3.24刺点过程
(7)重复以上步骤完成均匀地选取至少7个控制点,如图3.25所示。
图3.25均匀分布控制点
(8)控制点删除:
在控制点列表中选中需要删除的控制点,单击鼠标右键,在弹出的右键菜单中选择删除。
图3.26删除控制点
(9)控制点保存:
在控制点列表中用Shift键选中需要保存的控制点,单击鼠标右键在弹出的右键菜单中选择导出配准信息,在弹出的另存为对话框中将配准信息存储为.drfu格式文件。
图3.27导出配准信息命令
图3.28另存为配准信息
(10)控制点导入:
在控制点列表中空白处单击鼠标右键在弹出的右键菜单中选择导入配准信息,在弹出的打开对话框将配准文件加载到控制点列表中来。
图3.29导入配准信息
图3.30配准信息导入后
8、设置栅格数据重采样
图3.31栅格数据重采样
栅格/影像数据集在完成配准以后,原始图像由于拉伸、旋转等操作,像元的中心位置发生了移动,此时需要对配准后的栅格/影像数据集进行重采样。
“栅格/影像”组中的模式下拉按钮,提供了三种重采样的方法:
最邻近法、双线性内插法、三次卷积法。
9、计算误差
在“配准”选项卡的“运算”组中,单击“计算误差”按钮,则应用程序会自动计算所有控制点的误差,包括X残差、Y残差以及均方根误差。
计算结果会显示在控制点列表中“X残差”、“Y残差”以及“均方根误差”列中,如下图所示。
配准窗口中的状态栏会同时会输出总误差值,即各个控制点的均方根误差之和。
图3.32控制点列表
10、执行配准
当控制点的误差在允许范围内时,在“配准”选项卡的“运算”组中,单击“配准”按钮,执行配准操作。
当配准方式为线性配准或者二次多项式配准,在配准结束后,应用程序会在输出窗口中显示配准转换的公式及各个参数值,以便用户查阅。
如图3.24所示。
图3.33输出窗口
如果用户选择保存配准结果,则此时在相应数据源结点下会生成配准结果数据集。
11、为校正后影像赋投影信息
(1)右击“J48G024008_校正后”在弹出的右键菜单中点击属性,弹出属性对话框。
图3.34右键选择属性
图3.35属性对话框
(2)选择投影信息,可以看到“J48G024008_校正后”的坐标系信息为“平面坐标系无投影参数”,再点击重新设定坐标系。
图3.36重新设定坐标按钮
(3)在弹出的投影设置对话框中,点击“重置”按钮,将SuperMap的坐标系统恢复到默认设置。
图3.37重置坐标系统
(4)在默认投影坐标系中没有基于西安1980地理坐标系的投影信息。
因此这里需要采用自定义坐标系来实现本数据投影信息的选取,点击任投影坐标系下任意一个文件夹节点,会发现重置附近有一个自定义坐标系按钮。
图3.38自定义坐标系按钮
(5)点击自定义坐标系按钮,弹出自定义投影坐标系对话框。
图3.39自定义投影坐标系对话框
(6)在投影方式选择“GaussKruger”,在地理坐标系类型处选择“Xian80”,在名称处选择Xian1980DegreeGk34N,单击确定按钮完成投影方式的自定义。
图3.40自定义投影坐标系
(7)下拉投影设置对话框右侧的滚动条至底,选中自定义的Xian19803DegreeGk34N,单击应用按钮完成配准图层的投影信息设置。
图3.41选择自定义投影坐标系
(8)至此完成了配准图层的投影信息设置。
图3.42投影信息设置后
2.4.2基于参考图层地图配准
图3.43地图窗口
单击功能区→“数据”选项卡→“配准”组→“新建配准”按钮,弹出“配准数据设置”对话框中。
图3.44配准数据设置对话框
参考图层数据集为交通线;
图3.45配准数据设置
图3.46配准窗口
单击功能区→“配准”选项卡→“算法”组→在配准算法下拉框中选择常用的二次多项式配准算法。
图3.47选择配准算法
这里使用
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验二 地图配准 实验 地图