码迷,mamicode.com
首页 > 其他好文 > 详细

ArcGIS影像配准与空间配准

时间:2014-12-18 18:13:35      阅读:415      评论:0      收藏:0      [点我收藏+]

标签:style   blog   http   ar   io   os   使用   sp   on   

ArcGIS影像配准与空间配准

   

   

ArcGIS影像配准与空间配准 

   

   

地图配准可分为影像配准和空间配准。影像配准的对象是raster图,譬如TIFF图。配准后的图可以保存为ESRI GRID, TIFF,或ERDAS IMAGINE格式。空间配准(Spatial Adjustment)是对矢量数据配准。

   

一、影像配准

   

在ArcGIS中配准:

   

1.打开ArcMap,增加Georeferencing工具条。

2. 把需要进行纠正的影像增加到ArcMap中,会发现Georeferencing工具条中的工具被激活。在view/data  frame  properties的coordinate  properties中选择坐标系。如果是大地(投影)坐标系选择predefined中的Projected coordinate system,坐标单位一般为米。如果是地理坐标系(坐标用经纬度表示)表示则选择Geographic coordinate system。

3.纠正前可以去掉"auto adjust"前的勾。在校正中我们需要知道一些特殊点的坐标。如公里网格的交点,我们从图中均匀的取几个点,不少于7个。在实际中,这些点要能够均匀分布在图中。

4.首先将Georeferencing工具条的Georeferencing菜单下Auto Adjust不选择。

5.在Georeferencing工具条上,点击Add Control Point按钮。

6.使用该工具在扫描图上精确到找一个控制点点击,然后鼠标右击,Input X and Y输入该点实际的坐标位置。采用地理坐标系时应输入经纬度,经纬度用小数表示,如110°30‘30‘应写成 110.508(=110+30.5/60)。

7.用相同的方法,在影像上增加多个控制点,输入它们的实际坐标。

8.增加所有控制点后,在Georeferencing菜单下,点击Update Display。

9.更新后,就变成真实的坐标。

10.在Georeferencing菜单下,点击Rectify,将校准后的影像另存。

   

在arcview 中进行配准:

   

软件准备:arcview,必须安装了image wape扩展模块

数据准备:栅格图像,必须有坐标(大地坐标或经纬座标)

步骤:

1、找到栅格图像上的坐标点,至少要四个,坐标点分布应均匀。

2、打开arcview,单击chart,单击new,创建一张新表

3、输入X,Y坐标(栅格图上的点坐标),

4、点击保存。打开view窗口,单击theme,add event theme,设定字段为刚刚新建的表格里的X,Y字段。

5、这时导入到视图中,你可以看到几个点。

6、再convert to shp,把它转换为点.shp文件,加入到视图中。

7、在view中,propertise里设置你想要的投影方式。

8、这时控制点文件就做好了。

9、接下来,关闭所有的窗口,单击file-extension,image wape,点击此模块,这时菜单栏多了一行image wape

10、点击image wape——image wape session,选择你想要配准的栅格图像到第一个框。再选择刚刚生成的点文件到第二个框。

二、空间配准(ArcGIS)

1、首先在Arc catalog下面对需要配准的矢量图设置投影方式(右键点击该图,在shapefile属性表的shape字段属性中,定义相应的投影方式;

2、通过excel或记事本建立一个文本文件,输入几个控制点的坐标值,x表示经向,y为纬向,可以是公里网格坐标,和上面的定义的投影方式对应,保存;

3、在ArcMap里面打开矢量图,在"编辑器"(Editor)里设为"开始编辑";

4、在"空间配准"(Spatial Adjustment)里把"设置配准数据"(set adjust data)定为"选择这些层中的所有要素"(All features in these),然后在"链接"(link)里打开上面建立的控制点文件;

5、先双击控制点文件中的坐标值,再在矢量图中找到对应的点双击,即建立了第一个链接,依照此法,建立数个链接;

6、点击"空间配准"下面的"配准"(已由灰变黑可用),则配准完成。

   

   

   

   

   

zz地理坐标的精确配准

   

   

矢量配准用Spatial Adjustment

   

栅格配准用Georeferncing

   

地理坐标的精确配准

   

ArcMap里面配准

   

http://hi.baidu.com/sinogis/blog ... 7186264e4aea02.html

   

可能和很多朋友一样,我们只能拿到纸质地形图,还得自己扫描,之后再在ArcMap里面配准。对于刚开始接触的人来说,真是一头雾水。

   

我以前的误区:"刚接触ArcGIS,感觉挺简单的,不就先扫描,然后在ArcMap中加载数据,再Georeferencing配准。而实地采集的GPS数据就是通过excel转换到Access数据库后加载到arcmap里面,再display data。这样,点就投到图上了。"

   

呵呵,刚开始,还真以为是这样做呢。不过到Layout view界面,Insert Scale bar后,非常困惑,才发现怎么插入的比例尺和原图中的对不上啊。经过多方学习,才知道其中还有许多深奥的知识,任何一点没注意都会出差错。下面就与大家分享一下我的学习总结。

   

首先要知道什么是地理坐标和大地坐标,简单的说地理坐标就是球面坐标,而大地坐标是平面坐标。要深入了解这两者的区别,大家可以到网上搜相关资料。

   

地理坐标到大地坐标的过程叫做投影。

   

用1:50000地形图来举个例子,是北京1954坐标系,6度分带的高斯克吕格投影。

   

上面有两种坐标,一种是常说的方里网,就是间隔为1公里的网格。还有就是四角上的经纬度。

   

我们常用的GPS采集的数据都是WGS84坐标系的,是常见的度分秒格式。这里要特别注意,地形图上的经纬度不是WGS-84的经纬度,而是北京1954基准面下的经纬度。而且,这个经纬度也不是球面坐标,而是经过高斯克吕格投影的平面坐标。论坛上常提及的wgs-84到beijing 1954的转换就是这个问题,两者的换算涉及到图幅所在分带保密点的问题,转换参数并未公开。在我国,这两个坐标系同一点经纬度有几十米的差别。

   

了解了上面的知识后,就开始在arcmap中进行图像配准了,这里假设大家都知道基本操作步骤,因此我只注重于工作内容的阐述。

   

1, 首先把图复制到工作目录下

   

2, 在ArcCatalog中找到需要配准的图,给他定义一个投影系统(注意和投影变换的区别)。这里我们用的是1:50000的地形图,它是基于北京1954坐标系,6度分带的高斯克吕格投影。在地形图方里网上可以看出本图幅位于哪个分度带,这里假设是19。因此我们要选择的是Beijing 1954 GK Zone 19.prj。同时目录里面还有一个Beijing 1954 GK Zone 19N.prj,这个是用于没有分度带号的。而我们的图幅是包括分度带号。这个可以在网上搜索 "arcgis 坐标文件" 来进一步了解。

   

3, 给配准的图定义了一个投影系统后,把图加载到ArcMap里面,打开GeoReferencing工具,直接利用方里网交点进行配准。这里要注意的问题是,地形图上的方里网坐标为公里,而需要输入的应该是米。所以要在方里网对应坐标后面加000。如地形图上读出一个交点为(19387, 3420),19387的19为分带号,也要一并输入,那么这个点应该输入(19387000, 3420000)。其余配准过程和网上教程一致。

   

4, 配准后可以在Lay Properity里面把显示单位改成度分秒,于是地图就以经纬度格式显示了。此时的经纬度是基于北京1954基准面的,可以把光标指向四角的经纬度标记,以核对配准的精确度。

   

到此为止,地形图已经精确配准了。接下来往里面加载GPS数据。

   

1, 用你能想到的方式把GPS接收机上的数据导出来,除了一个一个手工输入,呵呵。

   

2, 尽量建立起数据库吧,这样方便些

   

3, ArcGIS中加载所得到的数据库,用Defined Query来选择自己所需要显示的数据

   

4, 如果需要精确地在地形图上显示GPS数据,那么还需要一个WGS-84到北京1954坐标的转换。有转换参数最好,如果没有,但的确需要,就去网上查一下手持GPS接收机转换参数估算这篇文章。我们这里是假设不需要精确显示,直接把WGS-84数据当成北京1954来用。前面也说了,有一定误差,但是不大。因此直接跳到5。

   

5, Display X,Y Datas, 定义为WGS-84坐标系统。

   

(其实此处我有疑问,因为现在的地形图是投影后的图,按我的理解,也应该把是球面坐标的WGS-84数据进行投影变换后才能用。但是查了很多资料,都是在这里直接定义数据来源是WGS-84就可以了,也许是ArcGIS在投点的时候,可以自动进行变换。)

   

6, OK,现在点也投上去了,ArcGIS工作中最基础的工作也就完成了。

   

再转到Layout view 界面,插入一个比例尺,发现和地形图上的比例尺非常吻合:-),把显示单位改成kilometer,用measure tool 就可以量地形图上任意点的距离了,试了一下公里网格,非常精准。

   

   

   

   

   

   

地理坐标,大地坐标

   

   

地理坐标:为球面坐标。 参考平面地是 椭球面。坐标单位:经纬度

大地坐标:为平面坐标。参考平面地是水平面  坐标单位:米、千米等。

   

地理坐标转换到大地坐标的过程可理解为投影。  (投影:将不规则的地球曲面转换为平面)

在ArcGIS中预定义了两套坐标系:地理坐标系(Geographic coordinate system)投影坐标系(Projected coordinate system),

   

1、首先理解地理坐标系(Geographic coordinate system),Geographic coordinate system直译为

地理坐标系统,是以经纬度为地图的存储单位的。很明显,Geographic coordinate syst

em是球面坐标系统。我们要将地球上的数字化信息存放到球面坐标系统上,如何进行操作

呢?地球是一个不规则的椭球,如何将数据信息以科学的方法存放到椭球上?这必然要求

我们找到这样的一个椭球体。这样的椭球体具有特点:可以量化计算的。具有长半轴,短

半轴,偏心率。以下几行便是Krasovsky_1940椭球及其相应参数。

Spheroid: Krasovsky_1940

Semimajor Axis: 6378245.000000000000000000

Semiminor Axis: 6356863.018773047300000000

Inverse Flattening(扁率): 298.300000000000010000

然而有了这个椭球体以后还不够,还需要一个大地基准面将这个椭球定位。在坐标系统描

述中,可以看到有这么一行:

Datum: D_Beijing_1954

表示,大地基准面是D_Beijing_1954。

--------------------------------------------------------------------------------

有了Spheroid和Datum两个基本条件,地理坐标系统便可以使用。

完整参数:

Alias:

Abbreviation:

Remarks:

Angular Unit: Degree (0.017453292519943299)

Prime Meridian(起始经度): Greenwich (0.000000000000000000)

Datum(大地基准面): D_Beijing_1954

Spheroid(参考椭球体): Krasovsky_1940

Semimajor Axis: 6378245.000000000000000000

Semiminor Axis: 6356863.018773047300000000

Inverse Flattening: 298.300000000000010000

   

2、接下来便是Projection coordinate system(投影坐标系统),首先看看投影坐

标系统中的一些参数。

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

Linear Unit: Meter (1.000000)

Geographic Coordinate System:

Name: GCS_Beijing_1954

Alias:

Abbreviation:

Remarks:

Angular Unit: Degree (0.017453292519943299)

Prime Meridian: Greenwich (0.000000000000000000)

Datum: D_Beijing_1954

Spheroid: Krasovsky_1940

Semimajor Axis: 6378245.000000000000000000

Semiminor Axis: 6356863.018773047300000000

Inverse Flattening: 298.300000000000010000

从参数中可以看出,每一个投影坐标系统都必定会有Geographic Coordinate System。

投影坐标系统,实质上便是平面坐标系统,其地图单位通常为米。

那么为什么投影坐标系统中要存在坐标系统的参数呢?

这时候,又要说明一下投影的意义:将球面坐标转化为平面坐标的过程便称为投影。

好了,投影的条件就出来了:

a、球面坐标

b、转化过程(也就是算法)

也就是说,要得到投影坐标就必须得有一个"拿来"投影的球面坐标,然后才能使用算法

去投影!

即每一个投影坐标系统都必须要求有Geographic Coordinate System参数。

   

关于北京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=20 745 921.8m。

在Coordinate Systems\Projected Coordinate Systems\Gauss Kruger\Beijing 1954目录中,我们可以看到四种不同的命名方式:

Beijing 1954 3 Degree GK CM 75E.prj

Beijing 1954 3 Degree GK Zone 25.prj

Beijing 1954 GK Zone 13.prj

Beijing 1954 GK Zone 13N.prj

对它们的说明分别如下:

三度分带法的北京54坐标系,中央经线在东75度的分带坐标,横坐标前不加带号

三度分带法的北京54坐标系,中央经线在东75度的分带坐标,横坐标前加带号

六度分带法的北京54坐标系,分带号为13,横坐标前加带号

六度分带法的北京54坐标系,分带号为13,横坐标前不加带号

在Coordinate Systems\Projected Coordinate Systems\Gauss Kruger\Xian 1980目录中,文件命名方式又有所变化:

Xian 1980 3 Degree GK CM 75E.prj

Xian 1980 3 Degree GK Zone 25.prj

Xian 1980 GK CM 75E.prj

Xian 1980 GK Zone 13.prj

西安80坐标文件的命名方式、含义和北京54前两个坐标相同,但没有出现"带号+N"这种形式,为什么没有采用统一的命名方式?让人看了有些费解。

   

=======================================

大地坐标(Geodetic Coordinate):大地测量中以参考椭球面为基准面的坐标。地面点P的位置用大地经度L、大地纬度B和大地高H表示。当点在参考椭球面上时,仅用大地经度和大地纬度表示。大地经度是通过该点的大地子午面与起始大地子午面之间的夹角,大地纬度是通过该点的法线与赤道面的夹角,大地高是地面点沿法线到参考椭球面的距离。

方里网:是由平行于投影坐标轴的两组平行线所构成的方格网。因为是每隔整公里绘出坐标纵线和坐标横线,所以称之为方里网,由于方里线同时又是平行于直角坐标轴的坐标网线,故又称直角坐标网。

在1:1万——1:20万比例尺的地形图上,经纬线只以图廓线的形式直接表现出来,并在图角处注出相应度数。为了在用图时加密成网,在内外图廓间还绘有加密经纬网的加密分划短线(图式中称"分度带"),必要时对应短线相连就可以构成加密的经纬线网。1:2 5万地形图上,除内图廓上绘有经纬网的加密分划外,图内还有加密用的十字线。

我国的1:50万——1:100万地形图,在图面上直接绘出经纬线网,内图廓上也有供加密经纬线网的加密分划短线。

直角坐标网的坐标系以中央经线投影后的直线为X轴,以赤道投影后的直线为Y轴,它们的交点为坐标原点。这样,坐标系中就出现了四个象限。纵坐标从赤道算起向北为正、向南为负;横坐标从中央经线算起,向东为正、向西为负。

虽然我们可以认为方里网是直角坐标,大地坐标就是球面坐标。但是我们在一副地形图上经常见到方里网和经纬度网,我们很习惯的称经纬度网为大地坐标,这个时候的大地坐标不是球面坐标,她与方里网的投影是一样的(一般为高斯投影),也是平面坐标

   

   

   

   

   

   

ArcGIS中图象配准经验总结

   

   

可能和很多朋友一样,我们只能拿到纸质地形图,还得自己扫描,之后再在ArcMap里面配准。对于刚开始接触的人来说,真是一头雾水。

我以前的误区:"刚接触ArcGIS,感觉挺简单的,不就先扫描,然后在ArcMap中加载数据,再Georeferencing配准。而实地采集的GPS数据就是通过excel转换到Access数据库后加载到arcmap里面,再display data。这样,点就投到图上了。"

   

呵呵,刚开始,还真以为是这样做呢。不过到Layout view界面,Insert Scale bar后,非常困惑,才发现怎么插入的比例尺和原图中的对不上啊。经过多方学习,才知道其中还有许多深奥的知识,任何一点没注意都会出差错。下面就与大家分享一下我的学习总结。

   

首先要知道什么是地理坐标和大地坐标,简单的说地理坐标就是球面坐标,而大地坐标是平面坐标。要深入了解这两者的区别,大家可以到网上搜相关资料。

地理坐标到大地坐标的过程叫做投影。

用1:50000地形图来举个例子,是北京1954坐标系,6度分带的高斯克吕格投影。

上面有两种坐标,一种是常说的方里网,就是间隔为1公里的网格。还有就是四角上的经纬度。

我们常用的GPS采集的数据都是WGS84坐标系的,是常见的度分秒格式。这里要特别注意,地形图上的经纬度不是WGS-84的经纬度,而是北京1954基准面下的经纬度。而且,这个经纬度也不是球面坐标,而是经过高斯克吕格投影的平面坐标。论坛上常提及的wgs-84到beijing 1954的转换就是这个问题,两者的换算涉及到图幅所在分带保密点的问题,转换参数并未公开。在我国,这两个坐标系同一点经纬度有几十米的差别。中国3S吧 3s8.cn

   

了解了上面的知识后,就开始在arcmap中进行图像配准了,这里假设大家都知道基本操作步骤,因此我只注重于工作内容的阐述。

1, 首先把图复制到工作目录下

2, 在ArcCatalog中找到需要配准的图,给他定义一个投影系统(注意和投影变换的区别)。这里我们用的是1:50000的地形图,它是基于北京1954坐标系,6度分带的高斯克吕格投影。在地形图方里网上可以看出本图幅位于哪个分度带,这里假设是19。因此我们要选择的是Beijing 1954 GK Zone 19.prj。同时目录里面还有一个Beijing 1954 GK Zone 19N.prj,这个是用于没有分度带号的。而我们的图幅是包括分度带号。这个可以在网上搜索 "arcgis 坐标文件" 来进一步了解。

3, 给配准的图定义了一个投影系统后,把图加载到ArcMap里面,打开GeoReferencing工具,直接利用方里网交点进行配准。这里要注意的问题是,地形图上的方里网坐标为公里,而需要输入的应该是米。所以要在方里网对应坐标后面加000。如地形图上读出一个交点为(19387, 3420),19387的19为分带号,也要一并输入,那么这个点应该输入(19387000, 3420000)。其余配准过程和网上教程一致。

4, 配准后可以在Lay Properity里面把显示单位改成度分秒,于是地图就以经纬度格式显示了。此时的经纬度是基于北京1954基准面的,可以把光标指向四角的经纬度标记,以核对配准的精确度。中国3S吧 3s8.cn

   

到此为止,地形图已经精确配准了。接下来往里面加载GPS数据。

1, 用你能想到的方式把GPS接收机上的数据导出来,除了一个一个手工输入,呵呵。

2, 尽量建立起数据库吧,这样方便些

3, ArcGIS中加载所得到的数据库,用Defined Query来选择自己所需要显示的数据

4, 如果需要精确地在地形图上显示GPS数据,那么还需要一个WGS-84到北京1954坐标的转换。有转换参数最好,如果没有,但的确需要,就去网上查一下手持GPS接收机转换参数估算这篇文章。我们这里是假设不需要精确显示,直接把WGS-84数据当成北京1954来用。前面也说了,有一定误差,但是不大。因此直接跳到5。

5, Display X,Y Datas, 定义为WGS-84坐标系统。

(其实此处我有疑问,因为现在的地形图是投影后的图,按我的理解,也应该把是球面坐标的WGS-84数据进行投影变换后才能用。但是查了很多资料,都是在这里直接定义数据来源是WGS-84就可以了,也许是ArcGIS在投点的时候,可以自动进行变换。)

6, OK,现在点也投上去了,ArcGIS工作中最基础的工作也就完成了。:)

   

再转到Layout view 界面,插入一个比例尺,发现和地形图上的比例尺非常吻合:-),把显示单位改成kilometer,用measure tool 就可以量地形图上任意点的距离了,试了一下公里网格,非常精准。

   

ArcGIS影像配准与空间配准

标签:style   blog   http   ar   io   os   使用   sp   on   

原文地址:http://www.cnblogs.com/tuncaysanli/p/4172018.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!