GDAL：DEM数据分析

（1）地形分析

GDAL自带了gdaldem用来处理DEM数据，主要是包括山体阴影（Hillshade）、坡度（slope）、坡向（aspect）、色彩浮雕（Color-relief）、地形耐用指数Terrain Ruggedness Index (TRI)、地形位指数Topographic Position Index (TPI)、地表粗糙度Roughness。但是需要注意，DEM的经纬度和高程应该都是一样的单位，如果DEM的坐标系是地理坐标，那么就需要使用-s指定比例系数，将度转为米，这个系数约为111120，一度约111千米。

以创建hillshade为例，使用hillshade参数，利用之前合并的DEM数据（WGS 84坐标），创建一个具有良好阴影的8位栅格浮雕效果。

 # gdaldem，这里也可以不使用cmd /r
 cmd /r gdaldem hillshade merged.tif hillshade.tif -s 111120

并且gdaldem支持多种阴影算法，例如：组合阴影(-combined) 斜坡和斜面阴影的组合；多向阴影(-multidirectional) 从225度，270度，315度和360度方位角照明的山丘阴影的组合，以多向阴影为例：

cmd /r gdaldem hillshade merged.tif hillshade_multidir.tif -s 111120 -multidirectional

（2）DEM渲染

DEM高程还可以使用不同的颜色进行渲染，直接将渲染信息固定到tif中，color-relief参数可以创建包含颜色和高程范围的彩色地貌地图，如下所示的色彩信息保存到colormap.txt文件中，放在DEM同一个路径下。

1000,101,146,82
1500,190,202,130
2000,241,225,145
2500,244,200,126
3000,197,147,117
4000,204,169,170
5000,251,238,253
6000,255,255,255

使用代码进行渲染，渲染后看到一个三波段的影像。

cmd /r gdaldem color-relief merged.tif colormap.txt colorized.tif