数字高程模型生成热力图python代码

以下是使用Python和GDAL库生成数字高程模型的热力图的示例代码：

# 导入需要的库
import gdal
import numpy as np
import matplotlib.pyplot as plt

# 打开数字高程模型数据
dem = gdal.Open("dem.tif")

# 读取数字高程模型数据的高程值
dem_band = dem.GetRasterBand(1)
dem_array = dem_band.ReadAsArray()

# 设置热力图参数
heatmap, xedges, yedges = np.histogram2d(dem_array.flatten(), np.zeros_like(dem_array.flatten()), bins=100)

# 绘制热力图
extent = [xedges[0], xedges[-1], yedges[0], yedges[-1]]
plt.imshow(heatmap.T, extent=extent, origin='lower', cmap=plt.cm.jet)
plt.colorbar()
plt.show()

需要替换示例代码中的`dem.tif`为实际数字高程模型数据的文件路径。另外，还可以根据需要调整热力图的分析参数，如bins等。生成的热力图可以通过调整cmap参数来修改颜色方案。