### 如何使用 Python 创建栅格地图 #### 推荐使用的库 为了创建栅格地图，`numpy` 和 `rasterio` 是两个非常重要的库。`numpy` 可用于处理数组数据，而 `rasterio` 则提供了强大的功能来读取和写入地理空间光栅数据集[^2]。 #### 示例代码 下面是一个简单的例子，展示了如何利用这两个库创建一个基本的栅格文件： ```python import numpy as np import rasterio from rasterio.transform import from_origin def create_raster(file_path, width, height, transform, crs='EPSG:4326'): """ 创建一个新的栅格文件 参数: file_path (str): 输出栅格文件路径. width (int): 栅格宽度. height (int): 栅格高度. transform (affine.Affine): 像素坐标转换矩阵. crs (str or CRS object): 投影坐标系，默认为WGS84经纬度投影. 返回值: None """ # 构建随机数值作为示例栅格数据 array = np.random.rand(height, width) with rasterio.open( file_path, 'w', driver='GTiff', # GeoTIFF 驱动程序名称 height=height, width=width, count=1, # 波段数量 dtype=array.dtype, crs=crs, transform=transform) as dst: dst.write(array, indexes=1) if __name__ == "__main__": output_file = "example.tif" # 定义图像尺寸 img_width = 100 img_height = 100 # 设置左上角起点位置及像元大小 pixel_size_x = 0.01 pixel_size_y = -0.01 upper_left_x = -117.0 upper_left_y = 34.0 # 计算仿射变换参数 affine_transform = from_origin(upper_left_x, upper_left_y, pixel_size_x, abs(pixel_size_y)) create_raster(output_file, img_width, img_height, affine_transform) ``` 这段代码定义了一个名为 `create_raster()` 的函数，它接受几个参数以指定要创建的新栅格文件的位置、维度和其他属性。这里还提供了一组默认设置以便快速测试此功能的工作情况。注意，在实际应用中可能需要调整这些参数以适应具体需求。