GDAL栅格数据操作拼接

本文主要是介绍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 globtemp_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 简单操作

这篇关于GDAL栅格数据操作拼接的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！

---