### 使用Python批量合并GeoTIFF文件 对于批量处理和合并GeoTIFF文件的任务，`rasterio` 和 `gdal` 是两个非常强大的库。下面介绍一种利用这些工具来实现批处理的方式。 #### 利用Rasterio读取并写入GeoTIFF数据 当涉及到操作栅格图像时，`rasterio` 提供了一个简洁而直观的接口用于打开、编辑以及保存地理空间光栅数据集。为了加载多个GeoTIFF文件并将它们组合成单个输出文件，可以按照如下方式编写脚本： ```python import os import rasterio from rasterio.merge import merge from glob import glob # 定义输入路径模式匹配所有待拼接的GeoTIFF文件 input_pattern = "path/to/tiffs/*.tif" output_file = 'merged.tif' # 获取所有符合条件的文件列表 files_to_mosaic = glob(input_pattern) # 创建一个空列表存储已打开的数据源对象 mosaic_files = [] for file_path in files_to_mosaic: src = rasterio.open(file_path) mosaic_files.append(src) # 执行融合过程 mosaic_array, out_transform = merge(mosaic_files) # 关闭所有的数据源连接 for f in mosaic_files: f.close() # 设置元数据参数准备写出最终结果 out_meta = src.meta.copy() out_meta.update({"driver": "GTiff", "height": mosaic_array.shape[1], "width": mosaic_array.shape[2], "transform": out_transform}) with rasterio.open(output_file, "w", **out_meta) as dest: dest.write(mosaic_array) ``` 这段代码展示了如何通过遍历指定目录下的所有`.tif` 文件，并调用 `merge()` 函数完成实际的空间镶嵌工作[^1]。 #### GDAL命令行工具与Python API结合使用 除了纯编程方法外，GDAL还提供了丰富的命令行实用程序可以直接应用于shell环境中执行复杂的GIS任务。如果倾向于更灵活的工作流设计，则可以在Python中嵌套运行外部命令来进行大批量的操作。例如，可以通过subprocess模块调用 gdal_merge.py 工具: ```python import subprocess command = ['gdal_merge.py', '-o', 'merged_output.tif'] tifs = ["file1.tif", "file2.tif"] # 替换为真实的文件名数组 command.extend(tifs) result = subprocess.run(command, check=True) if result.returncode != 0: raise Exception('Failed to execute gdal_merge') ``` 这种方法允许用户快速测试不同的选项而不必修改核心逻辑，同时也支持多线程加速大范围区域内的影像合成效率。