[gdal-dev] Create a vrt accounting for multiple raster bands

[fargalaxy1]

Hi,

if i do

/gdal_translate -of VRT -gcp 447.0 419.0 764294.728 905306.256 -gcp 447.0
931.0 764294.728 904999.056 -gcp 447.0 1442.0 764294.728 904691.856 -gcp
448.0 1952.0 764294.728 904384.656 -gcp 448.0 2463.0 764294.728 904077.456 
input_rgb.tif input_rgb.vrt
/
the output  /input_rgb.vrt/ will have just one VRTRasterBand block, even if
the input tif is in rgb.

/ <VRTRasterBand dataType="Byte" band="1">
    <Metadata />
    <SimpleSource>
      <SourceFilename relativeToVRT="0">input_rgb.tif</SourceFilename>
      <SourceBand>*1*</SourceBand>
      [....]
    </SimpleSource>
  </VRTRasterBand>
/
but I need to have this block repeated for each band:

/ <VRTRasterBand dataType="Byte" band="1">
    <Metadata />
    <SimpleSource>
      <SourceFilename relativeToVRT="0">input_rgb.tif</SourceFilename>
      <SourceBand>*1*</SourceBand>
      [....]
    </SimpleSource>
  </VRTRasterBand>

 <VRTRasterBand dataType="Byte" band="1">
    <Metadata />
    <SimpleSource>
      <SourceFilename relativeToVRT="0">input_rgb.tif</SourceFilename>
      <SourceBand>*2*</SourceBand>
      [....]
    </SimpleSource>
  </VRTRasterBand>
 <VRTRasterBand dataType="Byte" band="1">
    <Metadata />
    <SimpleSource>
      <SourceFilename relativeToVRT="0">input_rgb.tif</SourceFilename>
      <SourceBand>*3*</SourceBand>
      [....]
    </SimpleSource>
  </VRTRasterBand>
/
How can I modify the /gdal_translate -of VRT/ [...]  command line to add in
the vrt a VRTRasterBand block for each band?

thank you for your help,

M.



--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Create-a-vrt-accounting-for-multiple-raster-bands-tp5082573.html
Sent from the GDAL - Dev mailing list archive at Nabble.com.