目录
1. 案例背景
2. 出图展示
3. 数据说明
topo_map.jpg:雷尼尔山国家公园的扫描地形图。雷尼尔山国家公园位于美国华盛顿州,地理坐标系是 WGS1984,投影是 UTM 分度带投影。雷尼尔山国家公园的经度范围在西经 121 度-西经 122 度之间。扫描图上的高程值单位是英尺,1 英尺=0.3048 米。
4. 要求
4.1 使扫描图有正确的坐标值和坐标系(20 分)
①扫描图上的横线和竖线是方里格网,间隔是 2km,左上角的交叉点坐标值是(596000,5194000 )。
②请创建一份正确的扫描地形图数据,命名为“TopoMap”。
③请在您提交的步骤文档中,一定要附上 RMS 误差表。
4.2 绘制冰川区域、等高线、高程点和河流数据(30 分)
请对地图进行数字化工作,绘制冰川区域、等高线、高程点和河流数据,具体要求如下:
①冰川区域和河流不需要额外添加属性字段;
②等高线只需要绘制计曲线(即加粗的等高线),请添加高程值属性字段;
③高程点需添加高程值属性字段。
4.3 创建 DEM 数据并制作地形图(30 分)
①创建 DEM 栅格数据,分辨率为 20 米,像元值代表高程(米)。
②使用创建好的矢量和栅格数据,重新制作地形图。
a.要求和扫描图尽可能一致,请仔细观察地图。
b.没有矢量化的露营点、首曲线、步道等矢量信息不需要体现。
c.山峰、冰川和河流的名字需要在地图上显示出来。
d.不需要排版出图,请您直接将地图截图放入提交的步骤文档当中。
4.4创建该地区的三维模型(20 分)
①请创建如下图所示的三维模型。请仔细观察,完全还原。
②示例图中我们看到的表面是遥感影像,本试题请使用扫描地形图代替。
5. 步骤
5.1 数据预处理
5.1.1 定义投影
首先加载显示原始数据,然后查看原始数据的空间参考信息,如下图所示。
从原始数据的属性中可以得知该扫描图不包含空间参考信息,既没有基准面,也没有坐标信息。波段数为3,像素类型为无符号整型、位深为8位。为了使扫描图具有正确的坐标值和坐标系,接下来需要对其定义投影和地理配准。
已知雷尼尔山国家公园的地理坐标系是 WGS_1984坐标系,投影是 UTM 分度带投影。雷尼尔山国家公园的经度范围在西经121度~西经122度之间。因此,可推算出雷尼尔山国家公园位于WGS_1984_UTM_Zone_10N坐标系。
定义投影使用到的是【定义投影】工具,输入数据集或要素类设置为扫描图数据,坐标系设置为WGS_1984_UTM_Zone_10N坐标系,如下图所示。
5.1.2 地理配准
定义完投只是让扫描图数据拥有了正确的坐标系,还没有拥有正确的坐标值,现在需要对其进行地理配准来使其拥有正确的坐标值。
由于扫描图上的方里格网间隔已知为2km,左上角的交叉点坐标值已知为(596000,5194000),且当前数据框所使用显示单位为米,故可由此推算用于地理配准的其它控制点的坐标值。根据方里格网间隔值和左上角交叉点坐标值计算另外3个角的交叉点点的坐标值分别为:右上角交叉点坐标(606000,5194000)、左下角交叉点坐标(596000,5192000)、右下角交叉点坐标(606000,5192000),将这4个点都作为控制点。
使用【地理配准】工具依次添加这四个控制点坐标,方法是把添加控制点的光标精确放置于控制点所在的像元位置,然后右键鼠标就可以输入控制点坐标了,如下图所示。
然后点击应用就可以进行地理配准了,最后还需要打开RMS误差表查看地理配准精度是否足够高,如果地理配准误差过大就需要重新进行地理配准。地理配准误差过大的原因可能是输入的坐标值与添加控制点时的光标所处的位置不对应,也可能是没有把图片放到足够大再去添加控制点导致位置偏差过大。最后,如果需要的话可以把控制点坐标导出做备份,如下图所示。
根据题目要求需要创建一份正确的扫描地形图数据,命名为“TopoMap”,因此还需要将地理配准后的扫描地形图数据导出并重命名,如下图所示。
5.2 矢量化地形要素
5.2.1 矢量化操作
该部分任务是绘制冰川区域、等高线、高程点和河流数据。结合题目要求,通过目视判读地形图可知冰川区域的几何类型为面,河流、等高线的几何类型为线,高程点的几何类型为点。通过目视判断可知等高线计曲线之间的等高距为400米,该信息可用于推算未被标注高程值的等高线对应的高程值。由于原始地形图中的高度单位为英尺,而地图文档中显示单位为米,因此需要进行高度单位的换算。(1英尺=0.3048米)
首先在数据库下新建要素数据集用来存储矢量化的冰川区域、等高线、高程点和河流数据,把要素数据集命名为DataSet,坐标系选择WGS_1984_UTM_Zone_10N坐标系,然后再在要素数据集下分别新建要素类,分别命名为冰川区域、等高线、高程点和河流,要素类类型分别为面、线、点、线,各要素类的坐标系也选择WGS_1984_UTM_Zone_10N坐标系,如下图所示。
然后就可以开始矢量化冰川区域、等高线、高程点和河流数据了,如下图所示。(注意,在矢量化等高线的时候每一条等高线需要分开进行矢量化,也就是说在矢量化完一条等高线以后右键选择的是完成草图而不是完成部件,因为后续还需要在属性表中为每一条等高线添加高程数据,如果不同的等高线矢量化为同一条线数据,那么就没有办法为等高线添加高程数据。)
5.2.2 换算高程单位
按照题目要求需要为等高线和高程点数据添加高程数据,那么就需要分别为等高线和高程点数据各新建两个字段,分别为整型字段“高程(英尺)”和双精度字段“高程(米)”。然后通过仔细判读扫描地形图的内容来获取和推算每一个高程点和每一条等高线的高程数据,并依次把以英尺为单位的高程数据录入“高程(英尺)”字段中。接下来需要换算高程单位,使用的是【计算字段】功能,右键“高程(米)”字段选择【计算字段】,按照1英尺=0.3048米来计算出“高程(米)”字段的数值,如下图所示。
5.3 创建DEM与制作地形图
5.3.1 创建DEM
首先,根据已有的等高线和高程点矢量数据创建不规则三角网TIN。其次,由TIN转出栅格,得到分辨率为20米的DEM栅格数据。创建不规则三角网TIN使用的是【创建TIN】工具,坐标系选择WGS_1984_UTM_Zone_10N坐标系,输入要素类分别输入高程点和等高线的矢量数据,高度字段都选择“高程(米)”字段,高程点的类型选择Mass_Points,等高线的类型选择Hard_Line,标签字段选择<None>,如下图所示。
由不规则三角网TIN得到DEM栅格数据使用的是【TIN转栅格】工具,输入TIN设置为上一步得到的不规则三角网TIN数据,输出数据类型选择浮点型,方法选择线性,采样距离选择像元大小,采样值输入20,如下图所示。
5.3.2 制作地形图
首先,通过插入文本在原始栅格地形图上分别找到山峰的位置、冰川区域与河流的文本标注位置,以文本形式标注并设置其标注字体、放置属性,使其和原始地形图基本保持一致。然后,标注图层中的要素,为高程点和等高线设置标注字段为“高程(英尺)”。接下来调整文本标注的字体、颜色和放置的方向与位置。最后,修改DEM栅格数据的色带,如下图所示。
5.4 创建三维模型
该操作使用的是ArcScene软件。
5.4.1 设置DEM基本高度属性
首先打开ArcScene软件,加载之前创建好的DEM数据,在DEM的图层属性窗口中选择“基本高度”选项卡,将“从表面获取的高程”设置为在自定义表面上浮动,并选择自定义表面为该DEM数据,点击应用,如下图所示。
5.4.2 提取栅格范围并拉伸显示
首先,使用【栅格范围】工具将DEM数据作为输入栅格,输出要素类型选择LINE,然后在其图层属性窗口中,勾选“拉伸图层中的要素”,并设置拉伸值为1000,设置拉伸方式为“将其用作要素的拉伸数值”,如下图所示。
5.4.3 获取底部面要素
上面得到的数据只有侧面而没有底面,现在需要设法获取底面的数据,如下图所示。
首先,使用【栅格范围】工具,将DEM数据作为输入栅格,输出要素类型选择POLYGON,得到和DEM数据的高度保持一致的面要素。然后,在其图层属性窗口中选择“基本高度”选项卡,并修改“从要素获取的高程”为常量值1000。最后调整侧面和底面的颜色属性,如下图所示。
5.4.4 制作三维地形图表面效果
将TopoMap.jpg地形图数据添加到ArcScene软件中,在其图层属性窗口中选择“基本高度”选项卡,设置自定义浮动表面为DEM数据,这样即可得到的雷尼尔山国家公园的三维模型,如下图所示。