数据准备:山东省30m分辨率的ASTER GDEM数据
所用软件:ArcGIS 10.2
一、Python实现
按照“箜——Kong的博客走了一遍流程https://liminlu.blog.youkuaiyun.com/article/details/8522725”,不过最后一点儿没实现,什么时候实现了再补上
1.生成山体阴影图-彩色地形图-合并
from osgeo import gdal, osr
from osgeo.gdalconst import *
from ReadWriteTIF import IMAGE
import numpy as np
# 获取影像信息
infoDEM = gdal.Info(r"D:\ProfessionalProfile\DEMdata\2_OutputArcMap\demsd0330UTM.tif")
# 生成山体阴影图像
hillshade = gdal.DEMProcessing(r'D:\ProfessionalProfile\DEMdata\6_DEMXuanRan\hillshadePy.tif',
r"D:\ProfessionalProfile\DEMdata\2_OutputArcMap\demsd0330UTM.tif", "hillshade",
format='GTiff', # 指定图像的格式
zFactor=1, azimuth=315,
a