编写Python如何生成数字地表模型的代码

生成数字地表模型的代码可以使用Python中的地理信息系统库，如GDAL或Rasterio。以下是一个示例代码：

import numpy as np
import rasterio

# 定义地图范围和分辨率
xmin, ymin, xmax, ymax = (0, 0, 10, 10)
res = 1

# 创建数组
cols = int(np.ceil((xmax - xmin) / res))
rows = int(np.ceil((ymax - ymin) / res))
arr = np.random.randint(0, 100, size=(rows, cols))

# 写入为GeoTIFF文件
transform = rasterio.transform.from_bounds(xmin, ymin, xmax, ymax, cols, rows)
crs = {'init': 'epsg:4326'}
with rasterio.open('output.tif', 'w', driver='GTiff', height=rows, width=cols, count=1, dtype=arr.dtype, crs=crs, transform=transform) as dst:
    dst.write(arr, 1)

这段代码会生成一个随机数字的地表模型，并将其保存为GeoTIFF格式的文件。