GDAL栅格数据操作拼接

GDAL栅格数据操作

使用python GDAL 对下载好的瓦片（一定处理后）进行拼接。

栅格数据操作

osgeo.gdal.BuildVRT(destName, srcDSOrSrcDSTab, **kwargs)
Build a VRT from a list of datasets.
从数据集列表生成 VRT。
VRT即虚拟栅格文件。

osgeo.gdal.Translate(destName, srcDS, **kwargs)
Convert a dataset.
转换数据集合。

示例

对文件夹下的所有栅格数据构建VRT。

import os
from osgeo import gdal
import glob

temp_dir = os.path.join(os.getcwd(), 'temp')
vrt_path = temp_dir + "/tiles.vrt"
input_pattern=temp_dir + '/*.tif'
gdal.BuildVRT(vrt_path, glob.glob(input_pattern))

可以使用其他工具对VRT进行操作，如下，对栅格数据进行拼接。

output_path=output_dir + '/merged.tif'
gdal.Translate(output_path, vrt_path)

参考

1. https://gdal.org/api/python/utilities.html
2. python GDAL 简单操作