当前位置：首页 > news >正文

基于地形数据计算山体阴影

news 来源：原创 2025/6/13 21:24:11

1 基础知识

3 种最常见的激光雷达派生数据产品：

Digital Terrain Model (or DTM)：地面高程。
Digital Surface Model (or DSM)：表面顶部(想象一下在森林的树冠上覆盖一张床单)。
Canopy Height Model (CHM)：地面与物体顶部离地之间的高度或残差距离。这包括树木、建筑物和地球表面任何其他物体的实际高度。CHM 通过DSM 中减去 DTM 来创建。

2 基于地形数据计算山体阴影

1 读取地形数据

dtm_dataset = rio.open(os.path.join('.\data','NEON_D17_TEAK_DP3_320000_4092000_DTM.tif'))
dtm_data = dtm_dataset.read(1)

2 计算山体阴影

山体阴影用于可视化给定指定光源的表面上每个像素的假设照明值（从 0 到 255）。

要计算山体阴影，需要照明源的天顶（高度）和方位角，以及地形的坡度和坡向。

山体阴影的公式为：

其中所有角度均以弧度为单位。

# function to caluclate hillshade
def hillshade(array,azimuth,angle_altitude):azimuth = 360.0 - azimuth x, y = np.gradient(array)slope = np.pi/2. - np.arctan(np.sqrt(x*x + y*y))aspect = np.arctan2(-x, y)azm_rad = azimuth*np.pi/180. #azimuth in radiansalt_rad = angle_altitude*np.pi/180. #altitude in radiansshaded = np.sin(alt_rad)*np.sin(slope) + np.cos(alt_rad)*np.cos(slope)*np.cos((azm_rad - np.pi/2.) - aspect)return 255*(shaded + 1)/2hs_data = hillshade(dtm_data,225,45)

DTM+山体阴影：

3 计算冠层高度

dsm_dataset = rio.open(os.path.join('.\data','NEON_D17_TEAK_DP3_320000_4092000_DSM.tif'))
dsm_data = dsm_dataset.read(1)
chm_data = dsm_data - dtm_data;