实验3、影像配准及矢量化.doc

1实验 3、影像配准及矢量化一、实验目的1. 利用影像配准(Georeferencing) 工具进行影像数据的地理配准2. 编辑器的使用（点要素、线要素、多边形要素的数字化） 注意：在基于 ArcMap 的操作过程中请注意保存地图文档二、实验准备基本概念：1) 地图—Map (ArcMap document)在 ArcGIS 中，一个地图存储了数据源的表达方式( 地图，图表, 表格) 以及空间参考在ArcMap 中保存一个地图时， ArcMap 将创建与数据的链接，并把这些链接与具体的表达方式保存起来当打开一个地图时，它会检查数据链接，并且用存储的表达方式显示数据一个保存的地图并不真正存储显示的空间数据！2) 数据框架—Data Frame在“New Map File(新建地图)”操作中，系统自动创建了一个名称为“Layers ”的数据框架在 ArcMap 中，一个数据框架显示统一地理区域的多层信息一个地图中可以包含多个数据框架，同时一个数据框架中可以包含多个图层例如，一个数据框架包含中国的行政区域等信息，另一个数据框架表示中国在世界的位置但在数据操作时，只能有一个数据框架处于活动状态。

在 Data View 只能显示当前活动的数据框架，而在 Layout View 可以同时显示多个数据框架，而且它们在版面布局也是可以任意调整的3）组图层-- New Group Layer有时需要把一组数据源组织到一个图层中，把它们看作 Contents 窗口中的一个实体例如，有时需要把一个地图中的所有图层放在一起或者把与交通相关的图层（如道路、铁路和站点等）放在一起，以方便管理4）数据层ArcMap 可以将多种数据类型作为数据层进行加载，诸如 AutoCAD 矢量数据DWG，ArcGIS 的矢量数据 Coverage、GeoDatabase、TIN 和栅格数据 GRID，ArcView 的矢量数据 ShapeFile，ERDAS 的栅格数据 ImageFile，USDS 的栅格数据 DEM 等注意：Coverage 不能直接编辑，要编辑需要将 Coverage 转换成 ShapeFile软件准备：ArcGIS Desktop ---ArcMap实验数据：海燕乡土地利用现状图.jpg，在发送的压缩包“Exec03”中可以找到2三、实验内容及步骤 第 1 步 地形图的配准——加载数据和影像配准工具所有图件扫描后都必须经过扫描配准，对扫描后的栅格图进行检查，以确保矢量化工作顺利进行。

 打开 ArcMap，添加“Georeferencing（影像配准） ”工具栏 把需要进行配准的影像“海燕乡土地利用现状图.jpg”增加到 ArcMap 中，会发现“Georeferencing”工具栏中的工具被激活 首先将 Georeferncing 工具条 Georeferncing 菜单下 Auto Adjust 不选择第 2 步 输入控制点在配准中我们需要知道一些特殊点的坐标通过读图，我们可以得到一些控件点――公里网格的交点，我们可以从图中均匀的取几个点一般在实际中，这些点应该能够均匀分布 在“Georeferencing”工具栏上，点击“Add Control Points(添加控制点)”按钮 使用该工具在扫描图上精确到找一个控制点点击，然后鼠标右击输入该点实际的坐标位置，如下图所示：3 用相同的方法，在影像上增加多个控制点（一般是 4 个） ，输入它们的实际坐标点击“Georeferencing”工具栏上的“View Link Table(查看链接表)”按钮注意：在链接表对话框中点击“Save…”按钮，可以将当前的控制点保存为磁盘上的文件，以备使用检查控制点的 Residual(残差)和 RMS(均方差)。

如果某点残差特别大，需删除并重新选取控制点转换方式设定为“1st Order Polynomial（或仿射）” 因为目前我们只用了四个控制点，ArcGIS 只能执行一个 1st Order Polynomial (affine)（二次或仿射变换）4第 3 步 设定数据框的属性 增加所有控制点，并检查均方差（RMS）后，在“Georeferencing”菜单下，点击“Update Display(更新显示)” 执行菜单命令“View”－“Data Frame Properties(数据框属性)” ，设定数据框属性在“General(常规)”选项页中，将地图显示单位设置为“Meters(米)”5在“Coordinate System”选项页中，设定数据框的坐标系统为“Xian 1980 GK CM 111E”（西安 80 投影坐标系，6 度分带，中央经线在东经 111 度，横坐标前不加带号） ，与待矢量化的影像图的坐标系一致由于我们的影像图是二手资料，确切的地理坐标参数没办法获取，只能根据坐标换算得到的地理坐标范围（大约北纬 20°13′－20°17′，东经 110°34′－110°38′）进行推断其中央经线为 111°，坐标系统因为是 2003 年左右制作地图转换而来的影像图，一般用的是西安 80 坐标系统。

 更新后，就变成真实的坐标预定义_Projected Coordinate System_Gauss Kruger_Xian 19806第 4 步 矫正并重采样栅格生成新的栅格文件 在“Georeferencing”菜单下，点击“Rectify(矫正)” ，打开 Save as 对话框，将更新后的影像另存其中参数 Cell Size 选择与更新前相同，Resample Type(重采样方法)有 3 种：Nearest Neighbor(最近邻法)、Bilinear Interpolation(双线性内插法)和 Cubic Convolution(立方卷积法)，一般对栅格地图来说采用最近邻法为好在Output Location 中指明输出图像存储路径，Name 中输入名称“海燕乡土地利用现状图5”，Format 采用 JPG 格式输出图像中保存了配准信息 后面我们的数字化工作是对这个配准和重新采样后的影像“海燕乡土地利用现状图 5”进行操作在数字化之前，我们需要建立一些数据层来存储从影像上提取的要素信息，下面介绍在 ArcCatalog 中如何新建数据层？第 5 步 在 ArcCatalog 新建数据层该数据采用的是西安 80 坐标系统、6 度分带。

1) 打开 ArcCatalog在指定目录鼠标右击，在“New”中，选择“Shapefile…”，如下图：7(2) 输入数据层 Name，如“行政界线” ，选择要素类型，如 Polyline，点按钮“Edit…”选择空间参照系在“Spatial Reference Properties(空间参考属性)”对话框中选择合适的坐标系统，点击“Select”按钮在(Projected Coordinate Systems 目录下，选择 Gauss Kruger--Xian 1980-- Xian 1980 GK CM 111E.prj)点击“Add” ，现在这些坐标系统信息应该如下图所示：(3) 点击“确定” ，再按“OK” 这样，我们就创建了一个线状要素的数据层注意：上面我们只是介绍了土地利用现状要素中线要素“地类界线”的建立方法，点、面要素的建立方法类似第 6 步 以海燕乡土地利用现状图为背景，进行分层数字化关于数据分层：由于一幅土地利用现状图包含的要素有许多（如地类图斑、行政界线、铁路、公路、单线河流等等） ，因此在通过数字化提取相应信息之前我们要对地图数据进行预处理——分层一般存储以上要素的图层有三个：地类图斑（面） 、行政界线（线）、线状地物（线） 。

不过，为了使地图上所有的线划要素均只数字化一次，数字化过程不直接构面（当然 Arcmap 中可以直接构面） ，而是先将现状图上所有的线要素数字化，这个过程也要考虑数据分层，根据上面的分析我们可分为如下三个图层：行政界线（XZJX） 、线状地物（XZDW） 、地类界线（DLJX） 以上三层线要素数字化完后，再利用 XZJX、XZDW 及 DLJX 层构面，就得到面要8素“地类图斑” 关于数字化的顺序：刚刚前面提到：地图上所有的线划要素均只数字化一次，如何实现呢？实际上在数字化之前，我们先确定好采集要素的顺序，即可做到地图上所有的线划要素均只数字化一次数字化的顺序是先采集行政界线、然后线状地物（道路、单线河） 、最后再地类界线注意：跟踪时要保证线的交点处尽量闭合或相交，目的是保证线构面时候的成功率关于数字化的步骤：(1) 打开 ArcMap 中，将已配准的影像“海燕乡土地利用现状图 5”和新建的线要素“行政界线”加入数据框中在主菜单中选择“File-Add Data…”，或者在standard 工具条中选择 ，按住 shift 键，选择“海燕乡土地利用现状图 5.JPG”、“行政界线.shp” ，得到结果如图：注意：在 ArcMap 中，创建新地图并向其中加载数据层时，第 1 个被加载的数据层的坐标系统就作为该数据组默认的坐标系统，随后加载的数据层，无论其坐标系统如何，只要含有坐标信息，满足坐标转换的需要，都将被自动的转换成该数据组的坐标系统。

当然，这种转换不影响数据层所对应的数据文件本身2) 在开始矢量化之前，保存一个地图文档 Exec03.mxd，这样每次开始矢量化之前只需双击“Exec03.mxd” ，即可打开 ArcMap 及工作数据3) 选取“Editor” ，打开编辑工具栏，如图：9(4) 选取数据层“行政界线” ，按 Editor 的下拉键，点取“Start Editing”，利用编辑工具中的“Sketch　tool” ，就可数字化行政要素的边界线前面讲过，线的交点处要注意尽量闭合或相交，目的是保证线构面时候的成功率下面介绍线相交的技巧：线的起始点应尽量跨越已经画好的线要素，如下图比如地类界线数字化应该这样操作下面介绍线闭合的技巧（进行捕捉环境设置即可）：线闭合可通过直接捕捉已经存在要素的结点、终点、边线来实现具体操作方法：在 editor 下拉菜单选择“Start Editing”（如果已经有图层处于编辑状态，则可不进行上述操作） ，然后在 editor 下拉菜单选择“snapping” ， 打开 Snapping Environment 对话框，如下：在 Snapping Environment 对话框的上部列表中设置数据层抓点特性：可以设置抓取结点、终点、边线；下部列表中编辑草图的三个钩都选上，如上图，在对图层起始点10XZJX 进行编辑时，会自动捕捉顶点（vertex） 、边（edge）和终点（end） 。

5) 添加要素的属性项点取要加要素的数据层，如“行政界线”，点右键，“Open Attribute Table”,出现属性表，再按“Option”中的“Add Field”,可增加所需的属性项，在 Name 中输入“界线类型” ，在 Type 中选“Text” ，Length 设为20（10 个中文字符长） ，用来存储各类行政界线的名称如下图：删除属性项，鼠标放在属性项上，点右键，出现下拉菜单，点“Delete Field”注意：当数据层处于图形编辑(Start Editing)状况下，“Add Field”变灰,不可用6) 增加属性值，设置数据层处于编辑状况，按 Editor 的下拉键，点取“Start Editing”点 Edit Tool ，选取某要素，点右键，出现下拉菜单，点“Attribute…”，如下左图进入属性编辑窗口，如下图，即可输入或修改属性值11注意：并不是所有数字化采集的线要素都要增加属性字段、赋属性值，只有行政界线和线状地物需要。