目录

上一主题

10.6. 利用图形建模自动执行复杂的工作流程

下一主题

11. 专题图制作

本教程由OSGeo中国中心维护。

10.7. QGIS创建和利用高程面模型

10.7.1. 数字高程模型

数字高程模型DEM（Digital Elevation Model）,以数字形式按一定结构组织在一起的，表示实际地形特征空间分布的一种实体地面模型， 是地形形状、大小和起伏的数字描述。Z=f（x,y),Z为（x,y）处的高程。用数字的形式模拟二维表面连续变化的属性特征，如温度、 降水量、土地利用、土壤类型等。

10.7.2. 读取和显示高程点数据

转换为shapefile文件的数字高程模型是高程点数据，间隔为0.4 s（约10m），因此可用于地形分析等。启动QGIS后， 从主菜单中选择 "Project-Open" ，然后打开我们创建的“空间分析”文件夹的“空间分析3.qgs”。 然后，从主菜单中选择 Layer ‣ Add Layer ‣ Add Vector Layer ，然后调出 "Data Source Manager/Vector" 窗口。

在此窗口中，首先在 Source 中选择 "File" ，然后将 Encoding 设置为 "UTF-8" 。接下来， 单击 "Vector Dataset(s)" 右侧的 ... 按钮以显示 "Open OGR Supported Vector Dataset(s)" 窗口， 并在窗口右下角的菜单中选择文件类型。"All files"在文件夹 "DEM" 中指定 "demlOm.shp" 作为“文件名”。 然后，单击"Open OGR Supported Vector Dataset(s)" 窗口中的 Open 按钮，然后单击 "Data Source Manager/Vector" 窗口中的 Add 按钮，在地图视图中显示高程点数据，如 图 10.106 。

_images/img_1_x5m.png

图 10.106 读取数据高程数据

10.7.3. 转换为高程表面模型

当选择主菜单 "Raster-Conversion-Rasterize (vector to raster)" 时，将显示 "Rasterize (vector to raster)" 窗口。 在此窗口中，在 "Input layer" 中选择 "dem10m" ，在“属性字段”中选择“高程”。

然后，单击 "Rasterized" 上的 ... 按钮，将出现 "Save file" 窗口，打开 "DEM" 。在这里，选择 "Geo Tiff（.tif，.TIFF）" 作为“文件类型”，单击“保存”按钮作为“文件名”作为 "dem1OmR.tif" 。输出尺寸或者显示消息“需要指定分辨率”， 单击 OK 按钮。

当您返回 "Rasterize (vector to raster)" 窗口时，选中“以地图单位像素为单位的栅格分辨率”， 并为 "Width/Horizontal resolution" 和 "Height/Vertical resolution" 输入 "000101111" 。 如上所述，这里处理的数值高程模型数据是0.4秒（约10米）的间隔。 由于需要在 "Rasterize (vector to raster)" 窗口中输入要以度为单位输入的数值， 将0.4秒转换为度（将0.4秒除以3,600秒（=1度））并计算约0.000111111。并将此数字放在输入字段中。

此外，要进行其他设置，请单击窗口右侧的 Edit 按钮。 为了将栅格转换后不存在高程点的像素的高度值（节点值）设置为 "-999999" ， 将空格字符添加到文本字符串 "-a_nodata-999999"（ "-999999" ）的末尾，不要忘记按回车确认。 “地图单位像素中的栅格分辨率”的 "Width/Horizontal resolution" 和 "Height/Vertical resolution" 输入字段中的位数有一个限制，因此可以直接编辑该文本字符串并提高精度。 在这种情况下，分辨率从 "000111111" 更正为 "0.0001111111" ，如 图 10.107 。

_images/img_2_x51.png

图 10.107 Rasterize (vector to raster)设置窗口

单击 Help 按钮可以显示文本字符串的详细信息。另外，如上所述，请注意， 如果此设置中的文件夹名称和文件名中包含全角字符，则会发生错误。

设置完成后，单击 OK 按钮，当转换完成消息显示一段时间后，单击 OK 按钮。 此外，单击 "Rasterize (vector to raster)" 窗口中的 Close 按钮时， 地图视图将显示图层 "demlOmR" 作为转换结果。

10.7.4. 颜色设置更改

双击 "demlOmR" 图层，调出 "Layer Properties" 窗口，选择选项卡菜单 "Symbology" ， 然后在 "Render type" 中选择 "Singleband pseudocolor" 。 在 "Min/Max Vaiue Settings" 中对 "Interpolation" 与 "Color ramp" 进行设置， 然后单击 OK 按钮，如 图 10.108 。 在地图视图中显示按海拔颜色编码的地图，如 图 10.109 。

_images/img_3_x05.png

图 10.108 颜色设置更改

_images/img_4_xaw.png

图 10.109 颜色设置结果

如果在工件中设置了另一个坐标系，请单击状态栏上的 CRS 按钮。 打开窗口，在左侧的选项卡菜单中选择 "CRS" ，然后首先检查“启用在线”CRS转换。 接下来，在“世界空间参考系统”屏幕上选择WGS84，然后单击 OK 按钮。