matlab求坡向,GDAL使用DEM数据计算山体阴影（Hillshade）

求坡向可以使用Matlab的DEM处理工具箱中的`demflow`函数。该函数可以计算坡向、坡度和累积高程等DEM相关参数。以下是一个简单的示例：

% 读取DEM数据
[Z,R] = arcgridread('dem.asc');

% 计算坡向和坡度
[aspect,slope] = demflow(Z,R,'aspect','slope');

% 可视化坡向
figure
mapshow(aspect,R,'DisplayType','surface')
colormap(jet)
colorbar

计算山体阴影可以使用GDAL的`gdaldem`工具。以下是一个简单的示例：

gdaldem hillshade dem.tif hillshade.tif -z 1.0 -s 1.0 -az 315.0 -alt 45.0

其中，`dem.tif`是输入的DEM数据，`hillshade.tif`是输出的山体阴影数据，`-z`指定DEM数据的垂直缩放因子，`-s`指定DEM数据的水平缩放因子，`-az`和`-alt`指定太阳的方位和高度角。