Python中如何创建数字高程模型

在Python中，可以使用GDAL库来创建数字高程模型。具体步骤如下：

1. 导入GDAL库：import gdal

2. 打开DEM数据集：dataset = gdal.Open("dem.tif")

3. 获取DEM数据集的投影信息和地理变换信息：proj = dataset.GetProjection()，geotransform = dataset.GetGeoTransform()

4. 获取DEM数据集的行数和列数：rows = dataset.RasterYSize，cols = dataset.RasterXSize

5. 读取DEM数据集中的高程值：band = dataset.GetRasterBand(1)，elevation = band.ReadAsArray(0, 0, cols, rows)

6. 创建一个新的高程模型数据集：driver = gdal.GetDriverByName("GTiff")，out_dataset = driver.Create("dem_out.tif", cols, rows, 1, gdal.GDT_Float32)

7. 将投影信息和地理变换信息写入新的高程模型数据集：out_dataset.SetProjection(proj)，out_dataset.SetGeoTransform(geotransform)

8. 将高程值写入新的高程模型数据集：out_band = out_dataset.GetRasterBand(1)，out_band.WriteArray(elevation)

9. 关闭数据集：dataset = None，out_dataset = None

这样就可以创建一个数字高程模型了。