ArcGIS中配准栅格图像并且矢量化(转载)
栅格图像的校正和坐标系确定
启动ArcMap,新建一个新工程,右键Layers选择Add Data…添加TIF图像,将出现如下提示(如果提示无法加载rester data时请安装ArcGIS Desktop SP3补丁),单击Yes确定,加载图像后提示图像没有进行配准,确定然后配准图像。
图像加载后即可看到图像内容,右键工具栏打开Georeferencing工具条,进行图像的配准工作,在配准之前最好先保存工程。
在File菜单下打开Map Properties编辑地图属性,Data Source Options可设置保存地图文件的相对路径和绝对路径。(这里选择相对路径以确保将工程复制到其他机器可用)。
配准前要先读懂地图,望都县土地利用现状图采用1954北京坐标系,比例尺1:40000,查阅河北省地图发现望都县位于东经115度附近,那么按6度分带属于20带中央经线117度,按3度分带属于38带。从图框看到的公里数发现没有带号,应该是公里数。
这里只找了4个点进行配置(可以找更多的点),从左到右从下到上,逆时针编号为1、2、3、4;在ArcMap中单击Georefercning工具条上的Add Control Ponit工具(先掉Auto Adjuest选项),添加4个点控制点。
然后编辑Link Table中的4个控制点的代表的公里数,然后单击“Georeferecning下拉菜单的Auto Adjuest”图像即进行校正这时可看到参差值这里是0.00175(Total RMS)非常小说明配准较为精确。单击Save按钮可将控制点信息保存到文件,单击Load按钮可从文件加载控制点坐标。
给校准后的地图选择适合的坐标系,右键Layers打开Properties对话框属性对话框选择投影坐标系,(Prokected Coordinate Systems)展开Predefined/ Prokected Coordinate Systems/Gauss Kruger/Beijing 1954下找Beijing 1954 GK Zone 20坐标系(高斯克里克投影20带无带号),单击确定保存工程;这时配准工作即完成,在状态栏就可以看到正确的坐标单位了。
最后保存校正重新生成采样数据,单击“Georeferencing”工具条的“Rectify”菜单矫正并重采样栅格生成新的栅格文件;
打开“Save As”对话框设置保存参数和保存的校正栅格图像格式、保存位置等信息;其中“Resample Type”指定重新采样的类型,有三种选项Nearest Neighbor(for discreate data)自然邻近内插(不连续数据)、Bilinear Interpolation(for continuous data)双线性内插(连续数据)、Cubic Convolution(for continuous data)立方卷积内插(连续数据);这里选择了第二项,确定生成新的栅格数据,然后加载到ArcMap中进行下一步矢量化工作。
栅格图像的矢量化工作
单击ArcMap工具条上的ArcCatalog按钮打开ArcCatalog程序(ArcGIS的地理信息资源都这里完成创建、删除、复制等管理工作),在Catalog树下展开工程所在位置,鼠标右键菜单中选择New子菜单的Shapefile…新建一个Shape格式的地理要素文件(地理要素可存储为其他格式);
在“Create New Shapefile”对话框中给新的要素命名,在“Feature Type”要素类型下拉列表框中选择创建要素的类型(一个Shape文件只能表示一种要素),如:Ponit点、Polyline多边型线、Polygon多边形面和MultiPoint、MultiPatch,这里我们创建一个“Ponit点”文件;在“Spatial Reference”框中没有指定坐标系,因此单击“Edit…”按钮给新建的要素类指定坐标系;
单击“Edit…”按钮后打开“Spatial Reference Properties”对话框,然后单击“Select…”按钮选择一个合适的投影坐标系,这里我们选择“Beijing 1954 GK Zone 20”坐标系,最后确定返回完成Shapefile的创建;
回到ArcMap中添加新建的Shapefile点文件到Layers下;
将栅格图对准到矢量图上,先将两个图层打开,并且将georeferencing工具条调出来(在标题栏空白处右键即可找到)。然后在矢量图上分别在上下左右中间等各找一个参照点(以易于分辨为准)并且纪录她们的坐标值。如果栅格和矢量的坐标系统不一样,或者两个的单位不一样,要将栅格的坐标转一下。使二者一致。一般要把参考点的经纬度的度分秒形式坐标改成带小数的度的坐标。这几个点坐标形式的转化可通过自己编写一个小程序来实现。
找好参照点(至少要五个)之后,在georeferencing中的layer选择栅格图层,点击add control point图标(一个一头红一头绿的双箭头图标)在栅格图层中找到你的参照点,单击之后再击右键选择inputX andY,在弹出的对话框中输入刚刚采集到的点的坐标即可,依次将各个参照点全部输入。即完成了配准的工作。
配准前的准备工作:
要配准的栅格和矢量图层必须具有相同的投应坐标系。如果不相同,要修改其中一个的投影坐标系。修改(添加)投影要在arccatlog里面进行。在目录中找到文件点击右键打开属性表,在Spatial Reference中点击编辑,打开对话框,点import导入参考图层即可加载参考图层的投影坐标系统。应用并确定即可。
1 图幅拼接针对ARCGIS 9.0。
a) 保证待拼接图层的属性项名称和类型都相同,不相同则修改(修改方法见属性表操作);
b) 本操作假定图层已做过接边,如没有接边请先在“arcedit”中使用“edgematch”工具接边;
c) 在“ArcToolbox”>>“Data Management Tools”的“General”中找到“Append”工具,以待拼接图层的任意一个图层作为目标图层,把其它图层加入这个图层,“Schema Type (Optional)”选择“TEST”,然后点击“OK”按钮;
d) 进入“ArcInfo Workstation”,在“Arc”提示符(Arc:)后依次输入“&wo …”,“ae”,“disp 9999”,“de all”,“ec …”,“draw”,“ef arc”,“sel all”,“unsplit none”命令消除伪结点;
e) 这样目标图层就是拼接好的图层(但是仍然可能存在未消除的伪结点,有打断的线实体,可多使用“unsplit none”命令几次);
2 图幅裁剪针对ARCGIS 9.0。
a) 在“ArcTools”中选择“Command Tools”点击“OK”,出现“Command Tools”菜单;
b) 在菜单中选择“Analysis”>>“Map Sheets”的“Clip”工具;
c) 依次填写“Clip feature”、“Input cover”、“Clip cover”、“Output cover”等参数,点击“OK”按钮即完成
注意:
如使用“ArcToolbox”中的“Clip”工具则ARCGIS有可能会自动改变属性某些属性项的类型。
3 属性表操作针对ARCGIS 9.0。
在“ArcCatalog”、“ArcMap”、“ArcInfo Workstation”中都可以对图层属性表进行操作,推荐使用“ArcInfo Workstation”。
a) 进入“ArcInfo Workstation”,在“Arc”提示符(Arc:)后输入“tables”进入属性表操作工具(提示符变成Tables:);
b) 使用“select …”(如“select terlk.aat”)命令选择要编辑的属性表;
c) “items …”命令显示所有的属性字段定义,“additem …”增加属性项,“dropitem …”删除属性项,“alter”改变属性项;
d) “cal …”命令可以录入或修改属性项值;