相当于gdalbuildvrt的Python

有没有一种方法可以使用GDAL Python绑定执行与gdalbuildvrt实用程序相同的任务？到目前为止，除了创建单个数据集的vrt和手动编辑xml之外，我没有找到任何其他方法。我想从多个栅格创建vrt（本质上是执行镶嵌）。使用纯Python是否可能？我的另一个选择是使用子进程简单地调用gdalbuildvrt。

老实说，通过在或中使用gdalbuildvrt可以更容易地做到这一点。subprocessos.system

如果您希望通过Python做到这一点，那就可以做到。使用GDAL Python中的标准数据集创建方法，我们可以轻松创建基本数据集VRT。

from osgeo import gdal

drv = gdal.GetDriverByName("VRT")
vrt = drv.Create("test.vrt", x_size, y_size, 0)

请注意，我们最初创建的数据集没有波段。从有关VRT 的文档中可以看出，VRT数据集是少数可以接受AddBand参数的数据集类型之一。

vrt.AddBand(gdal.GDT_Float32)
band = vrt.GetRasterBand(1)

现在，对于每个波段，我们必须手动设置元数据项：

simple_source = '<SourceFilename relativeToVRT="1">%s</SourceFilename>' % source_path + \
    '<SourceBand>%i</SourceBand>' % source_band + \
    '<SourceProperties RasterXSize="%i" RasterYSize="%i" DataType="Real" BlockXSize="%i" BlockYSize="%i"/>' % (x_size, y_size, x_block, y_block) + \
    '<SrcRect xOff="%i" yOff="%i" xSize="%i" ySize="%i"/>' % (x_offset, y_offset, x_source_size, y_source_size) + \
    '<DstRect xOff="%i" yOff="%i" xSize="%i" ySize="%i"/>' % (dest_x_offset, dest_y_offset, x_dest_size, y_dest_size)
band.SetMetadataItem("SimpleSource", simple_source)
band.SetMetadataItem("NoDataValue", -9999)

SetMetadatItem有两个参数，第一个是元数据项的字符串，第二个是项本身。这意味着您不能子集元数据项，因此对于数据源，必须将整个内容设置为字符串。

请注意，我们可以使用此方法来创建复杂的源（ComplexSource），其中包含值的查找表KernelFilteredSource，任意大小和形状的内核过滤器源（）和蒙版带（MaskBand）。

从GDAL 2.1开始，CLI工具可作为库函数使用，实际上，这就是CLI工具现在在内部调用的功能。

例如：

gdalbuildvrt -r cubic -addalpha my.vrt one.tif two.tif

等效于：

from osgeo import gdal

vrt_options = gdal.BuildVRTOptions(resampleAlg='cubic', addAlpha=True)
gdal.BuildVRT('my.vrt', ['one.tif', 'two.tif'], options=vrt_options)

在可用的CLI选项直接映射到的参数BuildVRTOptions，再加上有像进度回调一些额外可用。

@rcoup的答案仅对我有用，如果将其修改如下：

from osgeo import gdal 

vrt_options = gdal.BuildVRTOptions(resampleAlg='cubic', addAlpha=True)
my_vrt = gdal.BuildVRT('my.vrt', ['one.tif', 'two.tif'], options=vrt_options)
my_vrt = None

否则，文件不会写入磁盘。