坐标知识和ArcGIS配准.docx
- 文档编号:6824030
- 上传时间:2023-01-10
- 格式:DOCX
- 页数:11
- 大小:335.14KB
坐标知识和ArcGIS配准.docx
《坐标知识和ArcGIS配准.docx》由会员分享,可在线阅读,更多相关《坐标知识和ArcGIS配准.docx(11页珍藏版)》请在冰豆网上搜索。
坐标知识和ArcGIS配准
坐标知识和ArcGIS配准
ArcGIS配准坐标转换
正文
1、坐标系统
地理坐标,大地坐标
地理坐标:
为球面坐标。
参考平面地是椭球面。
坐标单位:
经纬度
大地坐标:
为平面坐标。
参考平面地是水平面坐标单位:
米、千米等。
地理坐标转换到大地坐标的过程可理解为投影。
(投影:
将不规则的地球曲面转换为平面)
在ArcGIS中预定义了两套坐标系:
地理坐标系(Geographiccoordinatesystem)投影坐标系(Projectedcoordinatesystem),
1、首先理解地理坐标系(Geographiccoordinatesystem),Geographiccoordinatesystem直译为地理坐标系统,是以经纬度为地图的存储单位的。
很明显,Geographiccoordinatesystem是球面坐标系统。
我们要将地球上的数字化信息存放到球面坐标系统上,如何进行操作呢地球是一个不规则的椭球,如何将数据信息以科学的方法存放到椭球上这必然要求我们找到这样的一个椭球体。
这样的椭球体具有特点:
可以量化计算的。
具有长半轴,短半轴,偏心率。
以下几行便是Krasovsky_1940椭球及其相应参数。
Spheroid:
Krasovsky_1940
SemimajorAxis:
6378245.000000000000000000
SemiminorAxis:
6356863.018773047300000000
InverseFlattening(扁率):
298.300000000000010000
然而有了这个椭球体以后还不够,还需要一个大地基准面将这个椭球定位。
在坐标系统描述中,可以看到有这么一行:
Datum:
D_Beijing_1954
表示,大地基准面是D_Beijing_1954。
有了Spheroid和Datum两个基本条件,地理坐标系统便可以使用。
完整参数:
Alias:
Abbreviation:
Remarks:
AngularUnit:
Degree(0.017453292519943299)
PrimeMeridian(起始经度):
Greenwich(0.000000000000000000)
Datum(大地基准面):
D_Beijing_1954
Spheroid(参考椭球体):
Krasovsky_1940
SemimajorAxis:
6378245.000000000000000000
SemiminorAxis:
6356863.018773047300000000
InverseFlattening:
298.300000000000010000
2、接下来便是Projectioncoordinatesystem(投影坐标系统),首先看看投影坐标系统中的一些参数。
Projection:
Gauss_Kruger
Parameters:
False_Easting:
500000.000000
False_Northing:
0.000000
Central_Meridian:
117.000000
Scale_Factor:
1.000000
Latitude_Of_Origin:
0.000000
LinearUnit:
Meter(1.000000)
GeographicCoordinateSystem:
Name:
GCS_Beijing_1954
Alias:
Abbreviation:
Remarks:
AngularUnit:
Degree(0.017453292519943299)
PrimeMeridian:
Greenwich(0.000000000000000000)
Datum:
D_Beijing_1954
Spheroid:
Krasovsky_1940
SemimajorAxis:
6378245.000000000000000000
SemiminorAxis:
6356863.018773047300000000
InverseFlattening:
298.300000000000010000
从参数中可以看出,每一个投影坐标系统都必定会有GeographicCoordinateSystem。
投影坐标系统,实质上便是平面坐标系统,其地图单位通常为米。
那么为什么投影坐标系统中要存在坐标系统的参数呢?
这时候,又要说明一下投影的意义:
将球面坐标转化为平面坐标的过程便称为投影。
好了,投影的条件就出来了:
a、球面坐标
b、转化过程(也就是算法)
也就是说,要得到投影坐标就必须得有一个“拿来”投影的球面坐标,然后才能使用算法去投影!
即每一个投影坐标系统都必须要求有GeographicCoordinateSystem参数。
关于北京54和西安80是我们使用最多的坐标系
先简单介绍高斯-克吕格投影的基本知识,了解就直接跳过,我国大中比例尺地图均采用高斯-克吕格投影,其通常是按6度和3度分带投影,1:
2.5万-1:
50万比例尺地形图采用经差6度分带,1:
1万比例尺的地形图采用经差3度分带。
具体分带法是:
6度分带从本初子午线开始,按经差6度为一个投影带自西向东划分,全球共分60个投影带,带号分别为1-60;3度投影带是从东经1度30秒经线开始,按经差3度为一个投影带自西向东划分,全球共分120个投影带。
为了便于地形图的测量作业,在高斯-克吕格投影带内布置了平面直角坐标系统,具体方法是,规定中央经线为X轴,赤道为Y轴,中央经线与赤道交点为坐标原点,x值在北半球为正,南半球为负,y值在中央经线以东为正,中央经线以西为负。
由于我国疆域均在北半球,x值均为正值,为了避免y值出现负值,规定各投影带的坐标纵轴均西移500km,中央经线上原横坐标值由0变为500km。
为了方便带间点位的区分,可以在每个点位横坐标y值的百千米位数前加上所在带号,如20带内A点的坐标可以表示为YA=20745921.8m。
在CoordinateSystems\ProjectedCoordinateSystems\GaussKruger\Beijing1954目录中,我们可以看到四种不同的命名方式:
Beijing19543DegreeGKCM75E.prj
Beijing19543DegreeGKZone25.prj
Beijing1954GKZone13.prj
Beijing1954GKZone13N.prj
对它们的说明分别如下:
三度分带法的北京54坐标系,中央经线在东75度的分带坐标,横坐标前不加带号
三度分带法的北京54坐标系,中央经线在东75度的分带坐标,横坐标前加带号
六度分带法的北京54坐标系,分带号为13,横坐标前加带号
六度分带法的北京54坐标系,分带号为13,横坐标前不加带号
在CoordinateSystems\ProjectedCoordinateSystems\GaussKruger\Xian1980目录中,文件命名方式又有所变化:
Xian19803DegreeGKCM75E.prj
Xian19803DegreeGKZone25.prj
Xian1980GKCM75E.prj
Xian1980GKZone13.prj
西安80坐标文件的命名方式、含义和北京54前两个坐标相同,但没有出现“带号+N”这种形式,为什么没有采用统一的命名方式?
让人看了有些费解。
大地坐标(GeodeticCoordinate):
大地测量中以参考椭球面为基准面的坐标。
地面点P的位置用大地经度L、大地纬度B和大地高H表示。
当点在参考椭球面上时,仅用大地经度和大地纬度表示。
大地经度是通过该点的大地子午面与起始大地子午面之间的夹角,大地纬度是通过该点的法线与赤道面的夹角,大地高是地面点沿法线到参考椭球面的距离。
方里网:
是由平行于投影坐标轴的两组平行线所构成的方格网。
因为是每隔整公里绘出坐标纵线和坐标横线,所以称之为方里网,由于方里线同时又是平行于直角坐标轴的坐标网线,故又称直角坐标网。
在1:
1万——1:
20万比例尺的地形图上,经纬线只以图廓线的形式直接表现出来,并在图角处注出相应度数。
为了在用图时加密成网,在内外图廓间还绘有加密经纬网的加密分划短线(图式中称“分度带”),必要时对应短线相连就可以构成加密的经纬线网。
1:
25万地形图上,除内图廓上绘有经纬网的加密分划外,图内还有加密用的十字线。
我国的1:
50万——1:
100万地形图,在图面上直接绘出经纬线网,内图廓上也有供加密经纬线网的加密分划短线。
直角坐标网的坐标系以中央经线投影后的直线为X轴,以赤道投影后的直线为Y轴,它们的交点为坐标原点。
这样,坐标系中就出现了四个象限。
纵坐标从赤道算起向北为正、向南为负;横坐标从中央经线算起,向东为正、向西为负。
虽然我们可以认为方里网是直角坐标,大地坐标就是球面坐标。
但是我们在一副地形图上经常见到方里网和经纬度网,我们很习惯的称经纬度网为大地坐标,这个时候的大地坐标不是球面坐标,她与方里网的投影是一样的(一般为高斯投影),也是平面坐标。
2、影像配准(栅格配准)
在ArcGIS中配准:
1.打开ArcMap,增加Georeferencing工具条。
2.把需要进行纠正的影像增加到ArcMap中,会发现Georeferencing工具条中的工具被激活。
在view/dataframeproperties的coordinateproperties中选择坐标系。
如果是大地(投影)坐标系选择predefined中的Projectedcoordinatesystem,坐标单位一般为米。
如果是地理坐标系(坐标用经纬度表示)表示则选择Geographiccoordinatesystem。
3.纠正前可以去掉“autoadjust”前的勾。
在校正中我们需要知道一些特殊点的坐标。
如公里网格的交点,我们从图中均匀的取几个点,不少于7个。
在实际中,这些点要能够均匀分布在图中。
4.首先将Georeferencing工具条的Georeferencing菜单下AutoAdjust不选择。
5.在Georeferencing工具条上,点击AddControlPoint按钮。
6.使用该工具在扫描图上精确到找一个控制点点击,然后鼠标右击,InputXandY输入该点实际的坐标位置。
采用地理坐标系时应输入经纬度,经纬度用小数表示,如110°30'30'应写成110.508(=110+30.5/60)。
7.用相同的方法,在影像上增加多个控制点,输入它们的实际坐标。
8.增加所有控制点后,在Georeferencing菜单下,点击UpdateDisplay。
9.更新后,就变成真实的坐标。
10.在Georeferencing菜单下,点击Rectify,将校准后的影像另存。
3、矢量配准(矢量空间校正)
有时,我们拿到的数据是没有坐标信息的,为了将这些数据定位到一定的坐标下,需要对它们进行配准或者校正操作。
对栅格数据进行配准和对矢量数据进行校正的原理是一样的。
下面笔者利用ArcGIS对一幅shape格式的数据进行空间校正为例,将没有坐标信息的海图数据,校正到地理经纬度坐标下。
步骤1:
查看图幅范围
目录“未校正数据”中的矢量数据是没有坐标信息的。
将数据导入ArcMap中,查看右下角鼠标位置栏里,可看到显示的不是正确的地理经纬度坐标。
坐标位置信息
在图层ocldnta上点击右键,选择“Data-ViewMetadata”,在spatial页面中,可查看ocldnta图层的图幅范围。
查看图层元数据
图层四角范围
步骤2:
建立控制点文件
事先得知该图幅的对应经纬度范围为99.023293~148.031029°E,-9.939420~53.060567°N,由此可推算出如下四个用于校正的控制点:
10.0000000.00000099.023293-9.939420
267.5415340.000000148.031029-9.939420
367.54153497.370674148.03102953.060567
40.00000097.37067499.02329353.060567
将以上控制点复制到文本文件中,保存为“控制点文件.txt”,待用于空间校正之中。
步骤3:
启动ArcMap的编辑器和空间校正模块
启动ArcMap中的编辑器Editor,点击StartEditing。
启动空间校正模块SpatialAdjustment,点击SetAdjustData,选择Allfeaturesintheselayers,以便对所有图层数据一并进行校正。
步骤4:
利用控制点文件进行空间校正
在SpatialAdjustment中,点击Links-OpenLinksFile…,在打开文件对话框中,选择“控制点文件.txt”。
空间校正方法采取仿射校正方式Affine。
在空间校正模块中打开控制点文件
选择控制点文件
采取仿射校正方式
在SpatialAdjustment中,点击Adjust,实行空间校正,得到下图所示结果:
校正完毕,此时能够正确显示经纬度坐标了。
在编辑器Editor中,点击SaveEditors,保存校正结果。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 坐标 知识 ArcGIS 配准