0

I have two images with 3 and 5 bands respectively. How do I retain each image's individual bands after merging?

from osgeo import gdal

img_list = ['img1.tif', 'img2.tif']

vrt = gdal.BuildVRT("merged.vrt", img_list, separate=True)
gdal.Translate('merge_img.tif', vrt)

The code above will result in merge image with 2 bands. I want it to be 8 bands, that is: image1 (3 bands) + image2 (5 bands) = 8 bands.

Mark Rotteveel
  • 100,966
  • 191
  • 140
  • 197
Umar Yusuf
  • 926
  • 2
  • 13
  • 30

1 Answers1

0

gdal.BuildVRT can't do this, instead I used gdal_merge as follow;-

# pip install gdal-utils
import osgeo_utils.gdal_merge

parameters = ['', '-o', output_file_path] + img_list + ['-separate', '-ot', 'Float32']
osgeo_utils.gdal_merge.main(parameters)
Umar Yusuf
  • 926
  • 2
  • 13
  • 30