[gdal-dev] Potential bug or misunderstanding. Please help

[Even Rouault]

Billy,

bandList is currently ignored in the separate=True case. This would be a 
reasonable enhancement request that you could file on 
https://github.com/OSGeo/gdal/issues/new

One of the workaround for now is to postedit the generate VRT file and 
substitute the <SourceBand>1</SourceBand> element with 
<SourceBand>2</SourceBand>

You could also use use gdal.Translate() to produce a VRT with the second 
band of each source file, and run gdal.BuildVRT() on those files, but 
that would produce temporary files. If you don't need to serialize the 
output_vrt_fname, but use the returned dataset object with other APIs in 
the same process, you could avoid actual file creation by using an empty 
name for those intermediate VRT files, get the returned dataset object, 
and pass an array with those dataset objects instead of 
lst_of_input_geotiff_files.

Even

Le 19/01/2022 à 16:51, Billy Babis a écrit :

> Hi there.
>
> Given 3 input geotiffs (each with 2 bands), I want to create 1 geotiff 
> with 3 bands (the 2nd band from each aforementioned geotiff). This 
> should be quite simple using gdal.BuildVRT and gdal.translate. But, 
> gdal.BuildVRT will not allow me to select the 2nd band from each input 
> image. The following code produces the following warning:
>
> *Code:*
> gdal.BuildVRT(output_vrt_fname, lst_of_input_geotiff_files, 
> separate=True, bandList=[2])
>
> *Warning:*
>  "fname.tif has 2 bands. Only the first one will be taken into account 
> for the -separate case"
>
> Of course, I do not want to only use the first band, I want to use the 
> 2nd band to create the VRT. Please let me know what I'm doing wrong 
> and how I can fix it.
>
> Best,
> Billy B
>
> -- 
> Billy Babis
> PhD Candidate: Environmental Economics & Policy (ABD)
> Computer Science B.S.
> 1-(908)-403-2006
>
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev

-- 
http://www.spatialys.com
My software is free, but my time generally not.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20220119/ec51a54e/attachment.html>