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)
参考
- https://gdal.org/api/python/utilities.html
- python GDAL 简单操作
