影像图的校正

hk163 · 发表于 2016-8-22 16:00

地形图的校正

——地形图的粗校正

本文是专供初学ArcGIS 10参考用的。

纸质地形图在扫描中会受热，产生膨胀，造成公里网格不是正常的10cm长，这种膨胀也不一定是均匀的，因而各公里网格大小不一定相等，公里线也不是正常的一条直线；一般而言，1：10000纸质地形图内框下边是水平的，这就意味着其公里网格横线同水平线有一个约5°的夹角，因而在GIS系统中应用扫描地形图是必须先校正的。

ArcGIS 10可以较好地校正扫描地形图。因大多数时候我们并不知道公里网与内框交点的坐标，ArcGIS 10也就无法直接校正地形图的边框，未校正边框的地形图，其边缘部分误差较大。

本文试图分地形图粗校正与地形图精校正两部分来介绍个人校正地形图的一点浅见。

file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image002.jpg下图是要举例校正的图。

1、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image004.jpg打开ArcMAP；

2、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image006.jpg点击；

3、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image008.jpg找到要校正的文件；

4、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image010.jpgfile:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image012.jpg点击“添加”；

5、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image014.jpg点击“是”；

6、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image016.jpg点击“确定”；

7、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image018.jpgfile:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image020.jpg点击“自定义”菜单，调出“地理配准”菜单；

8、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image021.giffile:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image023.jpg点击“地理配准”右边小三角符“ ”；

9、点击，取消“自动校正”前“√”；

10、点击“file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image025.jpg”；

从地形图左下角开始，逆时针方向校输入公里网交点坐标。file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image027.jpg先记住左下角公里网交点坐标，如“39416，3289”，放大左下角公里网交点在十字交点正中，

file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image029.jpg左键点击一下，出现“十”字点，然后点右键，选“输入x,y坐标”，出现如下图：

11、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image031.jpg输入“416000，3289000”，如下图：

这里要解释一下，在地形图上该点坐标是 “39416，3289”，为什么却输入“416000，3289000”呢？这是因为“39416”前面39是带号，对一个县域来讲通常只有一个带号，日常应用中可忽略，后面补三个“0”，是精确到米（原是千米），所有的公里网交点坐标均作上述处理。

12、点击“确定”，返回输入坐标状态，一个公里网交点坐标就输入完成了。

13、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image033.jpg如上所述，输入所有公里网交点坐标；

14、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image035.jpg点击“file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image037.jpg”，可以看见输入的坐标，检查一下，如有错误，可直接在这里改正；

15、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image039.jpg点击“更新显示”

file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image041.jpg这时，你看到的是一个空白，不要紧，点一下“file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image043.jpg”就出现了全图。如下：

这里还要说明一下，有的人建议在这里，要注意残差，残差过大的应重新输入坐标（点击要删除的链接，再点“file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image045.jpg”即可删除，然后重新输入坐标），个人认为：这对一个已校好的地形图来说，理当如此做，但对于未校正过的扫描图来说，本就要求校正地形图，如无残差，还需校正吗？

16、点击“自动校正”，将“√”打上；

17、点击菜单“视图”，选“数据框属性”；

file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image047.jpg

18、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image049.jpg在打开和窗口中选“坐标系”

19、打开“投影坐标系”，选Gauss Kruger；

file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image051.jpg

file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image053.jpgfile:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image055.jpg选BeiJing1954

20、选“Beijing 1954 3 Degree GK CM117E”；

这里需说明一下：目前林业部门用的“蓝图”大多是北京54坐标系，高斯-克吕格投影。1：10000的地形图是3度带（2.5万以上是6度带），3度带的中央经度算法为：带号×3（带号是地形图上经度坐标前两位数），如这幅图中央经度为：39×3=117。

Beijing 1954 3 Degree GK CM117E，是指北京54坐标系，3度带，中央经度117度，所以选它。

21、点“确定”；

22、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image057.jpg点“校正”；

file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image059.jpg

23、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image061.jpg将“重采样类型”改为“双线性”；如下图：

选好“输出位置”，命好“名称”，格式最好选“TIFF”，点击“保存”，一张地形图就校正好了。

24、 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtml1/01/clip_image063.jpg重启ArcMAP,将校正后的地形图载入：