网上相关文章不多，参考了这一篇https://blog.csdn.net/liminlu0314/article/details/8522725 可能旧了点，需要更新一下
1.关于gdaldem的帮助文档要看，才知道参数如何填写https://gdal.org/programs/gdaldem.html
示例
gdaldem hillshade 8933.dem hillshade.tif -z 5 -s 10 -az 315 -alt 45
gdaldem color-relief 8933.dem rain.txt color.tif
rain.txt是色表，格式跟globalmapper导出的clr一样，支持高程和百分比两种形式
0 153 102 255
200 0 0 255
400 0 255 0
600 255 255 0
800 255 102 0
1000 255 0 0
2.合并图像脚本，下载地址 https://github.com/postmates/gdal/blob/master/samples/hsv_merge.py
需要注意代码是基于python2版本的，如果是py3则会报错，原因是print语法变了，py3需要加一对括号
python hsv_merge.py color.tif hillshade.tif merge.tif
3.去除黑边
nearblack merge.tif -setalpha -of GTiff -o noblack.tif
