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

2022-01-19 Thread Billy Babis
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

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

2022-01-19 Thread 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 1 element with 2 You could a

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

2022-01-19 Thread Matt.Wilkie
@lists.osgeo.org Subject: Re: [gdal-dev] Potential bug or misunderstanding. Please help 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<https://imsva91-ctp.trendmicro.

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

2022-01-19 Thread Matt.Wilkie
Would an in-memory /vsi also work as an intermediate? Expressed as windows CMD batch file I have handy: gdalwarp -of vrt –t_srs epsg:3578 %1 /vsistdout | gdal_translate /vsistdin/ %2 I should add that I use this like: project-to-local.bat input.tif outfile.tif -Matt _