### 遥感图像裁剪 在 Python 中，可以利用 `rasterio` 和 `GDAL` 库来实现遥感图像的裁剪功能。以下是基于 `rasterio` 的裁剪方法： 通过读取影像并定义一个边界框或多边形区域，可完成裁剪操作。具体来说，`rasterio.mask.mask()` 方法支持按照矢量几何对象裁剪栅格数据。 #### 裁剪代码示例 ```python import rasterio from rasterio.mask import mask import json # 定义裁剪范围（GeoJSON格式） geojson_polygon = { "type": "Polygon", "coordinates": [ [(658185.0, -203985.0), (658185.0, -436815.0), (886215.0, -436815.0), (886215.0, -203985.0), (658185.0, -203985.0)] ] } with rasterio.open("input_image.tif") as src: out_image, out_transform = mask(src, [geojson_polygon], crop=True) out_meta = src.meta.copy() out_meta.update({ "driver": "GTiff", "height": out_image.shape[1], "width": out_image.shape[2], "transform": out_transform, }) # 输出裁剪后的文件 with rasterio.open("output_clipped_image.tif", "w", **out_meta) as dest: dest.write(out_image) ``` 此代码片段展示了如何使用 `mask()` 函数根据指定的 GeoJSON 多边形裁剪输入影像[^2]。 --- ### 遥感图像拼接 对于遥感图像的拼接，通常涉及多个场景的数据融合成一幅完整的地图。这可以通过以下两种方式实现： 1. 使用 GDAL 工具中的虚拟镶嵌 (`vrt`) 文件； 2. 利用 NumPy 数组手动堆叠影像像素值。 #### 基于 VRT 的拼接方法 GDAL 提供了一个命令行工具 `gdalbuildvrt` 来创建虚拟镶嵌文件 (.vrt)，随后将其导出为 TIFF 格式。 ##### 示例代码 ```bash gdalbuildvrt mosaic.vrt input_image_1.tif input_image_2.tif ... gdal_translate -of GTiff mosaic.vrt output_mosaic.tif ``` 如果需要纯 Python 实现，则可通过 `rasterio.merge.merge()` 方法完成相同的功能。 #### 合并代码示例 ```python import rasterio from rasterio.merge import merge files_to_merge = ["image1.tif", "image2.tif"] src_files_to_mosaic = [] for fp in files_to_merge: src = rasterio.open(fp) src_files_to_mosaic.append(src) mosaic, out_trans = merge(src_files_to_mosaic) # 更新元数据 out_meta = src.meta.copy() out_meta.update( {"driver": "GTiff", "height": mosaic.shape[1], "width": mosaic.shape[2], "transform": out_trans} ) # 将结果保存到新文件中 with rasterio.open("merged_output.tif", "w", **out_meta) as dest: dest.write(mosaic) ``` 这段脚本实现了两个或更多 TIF 图像的自动拼接，并生成一个新的地理配准输出文件[^1]。 --- ### 总结 - 对于裁剪任务，推荐使用 `rasterio.mask.mask()` 或者直接调用 GDAL 的 API。 - 拼接方面，既可以依赖 GDAL 创建虚拟镶嵌文件后再转存为实际影像，也可以借助 `rasterio.merge.merge()` 动态组合多张图片。 ---