使用GDAL对DEM进行彩色渲染
在之前的博客中，我们已经看到了gdal对dem数据的强大的处理功能，其中除了坡度坡向等，还有一个比较厉害的，就是使用DEM生成一个彩色的图像。之前关于这方面也翻译了几篇博客，详见《使用GDAL对DEM渲染处理流程》、《使用阿富汗和巴基斯坦地区的SRTM数据生成山体阴影和彩色地形图》和《使用gdaldem创建彩色地形图和坡度阴影——thematicmapping.org译文（三）》，不过这些都是使用的说明。
把gdaldem的源代码进行了整理，整理了一个使用dem数据生成彩色图像的函数，函数的定义如下。关于第四个参数，sColorMode是个枚举值，一共有三个值，默认的是使用线性插值方式，还有两个是最邻近和严格插值。线性插值就是根据你指定的颜色映射关系，把没有指定的高程值，使用颜色映射关系中最近的两个颜色按照线性插值进行处理，得到最后的颜色；而最邻近就是找到最近的一个颜色，而严格插值就是所有的高程数据都必须对应一个颜色值，要不然就是黑色的。
图1 原始DEM数据
