1、一、初始数据分析 本案例的要求为根据项目需求，进行数据搜集和空间分析，最终给该物种的潜在分布位置，最终制作一张珍稀物种待查栖息地位置分布图以及相应图表以供调查人员参考。 提供的数据基础 有一张热振森林公园界限图（ Coverage 格式，数据名称为“ fanwei” ） 、由国家 1： 100 万公路数据得出的研究范围内等级公路数据（ shapefile 格式，数据名称为“ road”）、数字高程数据一份（数据名称为“ dem”）以及一份研究范围栅格地图（其中包含居民点、景点、河流和植被类型数据，数据名称为“热振森林公园图”） 。 打开 Arcgis，新建项目，在未定义坐标系的空图层中打开四份
2、原始数据 ，可得到 数据的以下 信息：首先四份数据虽然为同一研究范围但是在同一图层中打开后却不在同一位置，其次分别查看四份数据的属性信息 ，如图 1.1， 可得 road数据有基于 1940 克拉索夫斯基椭球（ GCS_Krasovsky_1940）的地理坐标系，其范围显示也符合 GIS 案例介绍中关于热振公园范围的介绍（大致为东经 91度 23 分 91 度 42分，北纬 30 度 12 分 30度 30 分），因此 判断 该数据有正确的空间信息； fanwei 数据没有坐标系信息且其范围数 值很大，目前暂时无法判别数据的范围是否正确；热振森林公园图 .jpg 数据也没有坐标系参考，其范围数
3、值很大，但大体与 fanwei 数据相吻合； dem 数据有坐标系信息，其空间参考信息为：基于 1940 克拉索夫斯基 椭球中央子午圈经度 93的 横轴墨卡托投影（ Krasovsky_1940_Transverse Mercator） ，其范围数值与 fanwei 图层一致，且在数据框中两个数据 相互重叠，因此认为 fanwei 图层虽没有坐标系信息但其坐标数值正确。 （注：将没有空间参考信息的图层添加入数据框时会提示未知空间参考，如图 1.2。向不具有金字塔的 ArcGIS 应用程序添加栅格数据集时，系统将提示您构建金字塔。金字塔十分有用，因为金字塔提升了分辨率低于其全分辨率的栅格数据集的
4、绘制速度。创建金字塔有三种重采样技术：最 邻近 法 - 用于标称数据或具有色彩映射表（如土地利用或伪彩色图像）的栅格数据集。双线性插值法 - 用于卫星影像或航空摄影等连续数据。 三次卷积插值法 - 用于卫星影像或航空摄影等连续数据。它与双线性插值法类似；不过，它使用较大的矩阵对数据进行重采样。） 图 1.1 四份原始数据的空间信息属性 图 1.2 添加数据无空间参考时弹出对话框提示 图 1.3 首次添加栅格数据时提示是否为数据创建金字塔 二、实验数据预处理 本次案例提供数据中的 fanwei 数据为 Coverage 格式，在当前 Arcgis10.2的环境下无法对其进行投影等操作，因此需要将
5、其转换为 Shapefile 格式 后才能对其进行投影操作。其次，原始数据中的热振森林公园图栅格数据中有经纬度格网，但经纬度格网交点处经纬度数值与数据框右下角数值稍有差异，因此本着学习知识的态度而且尽量避免前人使用过数据对本次实验造成影响，预处理时删除栅格数据中除了栅格图以外的其他信息， 之后 再 对数据重新进行地理配准。 2.1 将 fanwei 数据转为 Shapefile 格式 方法一 （在 ArcCatalog中导出为 Shapefile） ： 直接在 ArcCatalog 中 右键点击 fanwei 图层中的 polygon 要素 ，选择“导出 转为 Shapefile” ，如图 2
6、.1， 弹出“要素类至要素”对话框，选择输入要素和输出位置 并 为输出的 Shapefile 文件命名 ， 点击确定 后即可完成转换 ，如图 2.2。 图 2.1 转为 Shapefile 方法一 图 2.2 要素至要素类对 话框 方法二 （在 ArcMap 中导出为 Shapefile） ： 先将数据添加至数据框中，右键单击数据，选择“数据 导出数据”选项打开“ 导出 数据”对话框， 如图 2.3， 选择导出所有要素，坐标系信息目前先不设置，直接选源数据（源数据无空间参考信息），点击“浏览”选择输出要素类位置时 打开 “ 保存数据 ” 对话框，选择保存位置，输入新文件名称，并将保存类型选择为
7、“ Shapefile”， 如图 2.4， 点击确定后即可将数据转换为 Shapefile 格式。 图 2.3 转为 Shapefile 方法二 图 2.4 保存数据对话框 将 fanwei 数据中的 Polygon 要素转换为 Shapefile 格式后，将 arc 要素也转换为了 Shapefile 格式， 实验中 很可能用得到。 注： Coverage 格式数据的要素类型包括基本要素类型、复合要素类型、辅助要素类型三大类，基本要素类型包括标注点 Label Point、线 Arc、多边形 Polygon三类，其中多边形 Polygon 表示面状区域，边界由 Arc 构成。 2.2 删除
8、“热振森林公园图”栅格数据 的多余信息 将栅格数据复制到该步骤对应文件夹中，仅保留“ JPG 文件”和“ XML 文档”格式的两个文件作为后 续试验原始数据，其余全部删除 ，此时若将栅格数据添加入数据框中，则栅格数据左上角第一个像元中心点的坐标将变为（ 0,0） ，如图2.5。 图 2.5 栅格数据左上角第一个像素中心坐标值为（ 0,0） 三、地理配准 一般通过以下方法获取栅格数据：扫描地图、收集航空像片和卫星影像。扫描的地图数据集通常不包含空间参考信息（嵌入于文件中或作为单独的文件）。航空摄影和卫星影像提供的位置信息通常不够充分，无法与其他现有数据完全对齐。因此，要将这些栅格数据集与其他空间
9、数据结合使用，通常需要将这些数据对齐或 配准 到某个 地图坐标系 ，这个过程称为地理配准。 本实验将对“热振森林公园图”进行地理配准。 本例中的栅格数据“热振森林公园图”便没有空间参考信息（不具有坐标系信息且预处理过后栅格图的坐标不正确，如图 3.1，因此需要进行地理配准操作 。 ） 图 3.1 未配准之前坐标信息不正确 3.1 方法一 （创建地理参考） ： 3.1.1 新建保存点位坐标文档 首先从栅格数据中可读出经纬格网中四个经纬度交点的经纬度信息，并将之写入 TXT 文本或者 Excel 文档中，如图 3.2所示（单位选用“十进制度”，本例中 ABCD 四点分别为栅格图像由左 上点、右上点
10、、右下点、左下点）。 图 3.2 在 Excel 文档中输入控制点坐标 3.1.2 显示地理参考点位并保存为 Shapefile格式 点击 ArcMap 菜单中的“文件 添加数据 添加 XY数据”选项，如图 3.3，弹出“添加 XY 数据 ”对话框，设置 X 字段为经度 L， Y字段为纬度 B，如图 3.4。或者可以右键单击内容列表中的图层，选择“添加数据”将 Excel 表添加入图层中，再右键单击 Excel 表格并选择“显示 XY 数据” ，如图 3.5，弹出“显示 XY数据”对话框中的设置内容与“添加 XY 数据”中相同，两种方法均可将 准确坐标值的点位导入 ArcMap 数据框中 ，导
11、入 ArcMap 数据框中的结果及其空间参考属性信息如图 3.7。 图 3.3 添加 XY 数据 图 3.4 添加 XY 数据对话框 3.5显示 XY 数据 3.5 显示 XY 数据 3.6 提示没有 Object-ID字段 此时创建的空间参考点文件具有栅格数据格网交点的正确 坐标数值信息但是不具备 坐标系信息，而且由于属性表没有 Object-ID 字段因此无法选择或编辑图中的点数据，但可以用将其导出为 ShapeFile 格式的方法为其自动穿件Object-ID 字 段和属性，之后便可以对数据进行查询、编辑和使用。 图 3.7 添加至数据框的 XY 数据 四个控制点在数据框中显示出之后按照
12、上文中“将 fanwei 数据转为Shapefile 格式”中的方法二步骤将四个控制点保存为 Shapefile 文件。 3.1.3 为地理参考图层文件建立空间参考信息 为控制点文件定义空间参考。使用 ArcToolBox 中 “数据管理工具 投影和变换 定义投影” 工具 将控制点 Shapefile 文件（ “ Pointshp” ）定义投影为道路矢量数据“ road”的地理坐标系（可以在选择坐标系时选择 导入坐标系 ，如图3.8，或者在同一数据框中同时打开 “ Pointshp” 和 “ road” 两份数据，选择坐标系是选择“图层”分类中 road 的空间参考信息 如图 3.9，或定义数
13、据框的空间参考为 “ road”数据 的地理坐标系，再打开 “ Pointshp” 数据，右键点击“ Pointshp” 数据，选择“数据 导出数据”选项，选择“使用数据框相同的坐标系” ，如图 3.10）。 图 3.8 导入坐标系 图 3.9选图层中数据坐标系 图 3.10 选择数据框坐标系 对“ Pointshp”文件进行投影。使用 ArcToolBox 中 “数据管理工具 投影和变换 要素 投影”工具将“ Pointshp”文件的地理坐标系投影为“ dem”数据所使用的投影坐标系，以便将栅格数据直接配准到“ dem”数据所用的投影坐标系，选择输出坐标系的方法与上一步“为控制点文件定义空间参考”中的方法相同 （此处选择导入 dem 空间参考方法有时会发生中央子午圈经度发生变化的现象，应多加留意），投影后文件命名为“ Pointprj” 。
