如何基于DEM数据来计算河道方向上的坡度

      工具：

      使用ArcGIS中【Polyline To Raster】以及【Raster Calculator】实现。注意：河道方向上的坡度是指仅考虑在河道方向的纵剖面上进行坡度计算，与常规的区域内坡度值的计算方法相区别。因此，必须提取出河道的高程值才能计算出河道方向上的坡度值。本例讲述内容包括矢量河道数据的栅格化处理、河道高程值的提取以及河道方向上坡度的计算。

      实验步骤：

      1.栅格化河道数据

      载入数据后，打开【Conversion Tools】|【To Raster】|【Polyline To Raster】，输入要素设置为“rivers”，值字段设置为“GRID_CODE”，输出栅格数据集设置为“str”，像元大小设置为与“dem”相同；确定，得到栅格化后的河道数据。注意：①河道栅格值由GRID_CODE字段决定，即等于1；在进行栅格分析时，必须保证所有栅格图层的分辨率一致。②建议将文件（str）存放在默认数据库，或者将文件夹设置为数据库，不然有些情况下转化会显示失败。操作如下图。

      2.提取河道高程值

      打开【Spatial Analyst Tools】|【Map Algebra】|【Raster Calculator】工具，输入公式【"str" * "dem"】，输出栅格设置为“str_dem”；确定，得到河道的高程数据。注意：第一步中生成的栅格河道图层中，河道栅格值为1，非河道栅格值为nodata，故执行公式计算后，河道栅格值变为高程值，而非河道栅格值仍为nodata。操作如下图。

      3.计算坡度数据

      打开【Spatial Analyst Tools】|【Surface】|【slope】工具，输入栅格设置为“str_dem”，输出栅格设置为“str_slope”；确定，得到河道方向上的坡度数据。注意：由于仅河道栅格处有值，非河道栅格值为nodata，在坡度计算时，nodata不参与计算，故结果坡度为河道方向上的坡度值。操作如下图。