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 就可以量地形图上任意点的距离了，试了一下公里网格，非常精准。