### 使用 Python 读取 DEM 文件中的高程数据 为了使用 Python 读取 DEM 文件中的高程数据，可以借助 GDAL 库来处理地理空间栅格数据。下面展示了具体的操作流程以及一段简洁的代码示例。 #### 安装 GDAL 包 GDAL 是一个用于读写栅格和矢量地理数据格式的强大库，在开始之前需要先安装此库。有两种方式可完成安装操作[^1]： - 方法一：通过 pip 工具直接在线安装； - 方法二：下载对应版本的 wheel 文件离线安装； #### 示例代码展示 这段代码实现了基于给定经纬度查询特定位置处的海拔高度的功能。这里以湖北省某一点为例进行了说明[^2]: ```python from osgeo import gdal def get_elevation_from_dem(dem_file_path, longitude, latitude): """ 获取指定经度纬度下的DEM高程值 参数: dem_file_path (str): DEM文件路径 longitude (float): 经度 latitude (float): 纬度 返回: float: 对应位置的高度值 """ try: # 打开DEM数据集 dataset = gdal.Open(dem_file_path) if not dataset: raise Exception("无法打开DEM文件") band = dataset.GetRasterBand(1) elevation_array = band.ReadAsArray() transform = dataset.GetGeoTransform() xOrigin = transform[0] yOrigin = transform[3] pixelWidth = transform[1] pixelHeight = -transform[5] col = int((longitude - xOrigin) / pixelWidth) row = int((yOrigin - latitude) / pixelHeight) return elevation_array[row][col] except Exception as e: print(f"发生错误:{e}") return None if __name__ == "__main__": file_path = './data/hubei_wgs84.tif' lat = 31.15 lon = 111.24 height = get_elevation_from_dem(file_path, lon, lat) if height is not None: print(f'({lon}, {lat}) 的高程为 {height} 米') ``` 上述函数 `get_elevation_from_dem` 接受三个参数——DEM 数据文件的位置、目标地点的经度与纬度，并返回该点对应的海拔高度。如果遇到任何异常情况，则会捕获并报告相应的错误信息。