地形图的拼接
地形图的拼接
2011-01-11 11:58
在Arcgis9.3 支持下进行的地形图的拼接。
地形图的数据是jpg 格式的图片，比例尺为国家基本地形图系列。
思路：配准------》裁切-----》拼接
1. 配准
（1）控制点
采用高斯平面直角坐标系下的坐标进行配准。选好控制点，一般选择公里网的交点作为控制点，其坐标也能很好的在地形图上读出。
X 坐标；Y 坐标。在Arcgis9.3里该格网点的坐标为（34364000，
2776000），单位为米。34为带号。在实际的高斯-克吕格投影下坐标为（2776000，34364000）单位为米。这是为什么呢？
控制点的要求：大于7个。在输入前面3个控制点时要非常严格精确才行，同时不能在一条直线上，后面的点会自动根据前面输入的点进行设置，只有进行微小调整就行。
检查残差看是否符合要求，删除删除残差特别大的控制点并重新选取控制点。待残差符合要求后保存控制点文件。
（2）数据框属性
添加采集完所有控制点并检查均方差（RMS ）后，在“Georeferencing （影像配准）”菜单下，点击“Update Display（更新显示）”。
执行菜单命令View （视图）/Data Frame Properties（数据框架属性）设定数据框架属性。在对话框“General （常规）”标签中，在“Map （地图单位）”下拉列表中选择Decimal Degrees ，将“Display （地图显示单位）”设置为“Meters （米）”。
在“Coordinate System （坐标系统）”标签中，设定数据框架的坐标系统为“Xian 1980 3 Degree GK CM 102E（西安1980投影坐标系，3度分带，中央经线东经102度）”，与原始地图的坐标系一致。
更新后，就变成真实的坐标，大家可从状态栏中看到鼠标所在处的坐标值。
（3）矫正Rectify
在矫正方式上选择2nd Order Polynomial，设置好输出位置和输出格式。
2. 裁切
地形图右边会有图例，当两幅图进行直接拼接时会挡住右半边的地形图。先对原始地形图进行裁切，去掉其图例部分。
在裁切之前，先要建一个面的图层。在Arccatalog 里建立。以地形图想要留下来的部分作为面图层的范围。
在Arctoolbox 里Spatial Analyst Tools/Extraction /Extract by Mask 这个工具来实现。假设原始地图的文件名为source ，建立的面图层的文件名为pol.
双击“Extract by Mask”，在对话框中进行如下设置：在“Input Rasters ”下拉列表中选择添加source 到文本框中；在“Input raster or feature mask data”下拉列表中选择pol ；在“Output raster ”文本框中输入D:\PAGIS\GIS02\Exec3\Temp\Extract_dem_1。完成后点按“OK ”完成对地形图的裁切。
3. 拼接
Data Management Tools/Raster/RasterDataset / Mosaic
双击“Mosaic ”，在对话框中进行如下设置：在“Input Rasters ”下拉列表中选择添加dem1和dem2到文本框中；在“Output Location”文本框中输入D:\PAGIS\；在“Raster dataset name with extension ”文本框中输入dem_he；在“Pixel type ”下拉列表中选择16_BIT_UNSIGNED。其它设置不改变。完成后点按“OK ”完成DEM 数据的拼接。