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. _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev