我有一个hdf5文件，其中包含两个数据集s_i和s_q，我需要在一个栅格中合并两个波段，一个波段用于s_i，一个波段用于s_q
现在我找到了一种方法，可以在vrt（'MEM'）中复制数据集，但这需要大量的内存和时间
`ds_tmpI = gdal.Open(fnameI)
ds_tmpQ = gdal.Open(fnameQ)`
tmp_ds = gdal.GetDriverByName('MEM').CreateCopy('', ds_tmpI, 0)
tmpQ = ds_tmpQ.ReadAsArray()
tmp_ds.AddBand()
tmp_ds.GetRasterBand(2).WriteArray(tmpQ)
vrt_fname = os.path.join(xxxxxxx)
gds = drv.CreateCopy(vrt_fname, tmp_ds)`
最终栅格有两个波段，每个波段对应一个数据集