DEM数据用python输出

DEM（Digital Elevation Model）数据是一种数字地形模型，常用于地理信息系统中表示地面高程信息。在Python中处理DEM数据并输出，通常会涉及到读取、操作（如计算坡度、面积等）和可视化。

1. **读取DEM数据**：可以使用GDAL库（Geospatial Data Abstraction Library）或者像Rasterio这样的专门处理遥感数据的库来读取各种格式的DEM文件，例如.TIF或.HGT。

from rasterio import open

dem_file = "path_to_your_dem.tif"
with open(dem_file, 'r') as src:
    dem_data = src.read(1)

2. **处理和分析**：对读取的数据进行处理，比如计算栅格单元的平均值、统计高程分布等。numpy库常常在这个阶段发挥作用。

import numpy as np
average_elevation = np.mean(dem_data)

3. **输出结果**：你可以将处理后的数据保存回文件，或者直接绘制地图。matplotlib和geopandas都是常用的绘图工具。

import matplotlib.pyplot as plt
plt.imshow(dem_data, cmap='terrain')
plt.colorbar()
plt.show()

# 或者保存到新的TIF文件
with rasterio.open("output_dem.tif", "w", driver="GTiff", ... , height=dem_data.shape[0], width=dem_data.shape[1]) as dst:
    dst.write(dem_data, 1)