[gdal-dev] Question regarding tif mosaic (RGB and Gray)

[Stoermer, Jon-Pierre]

Hello all,
Thanks in advance for your time. This may be a stupid question but I can't
figure out how to mosaic two geotiff files with different color bands? Is
this possible? I get an error using gdalbuildvrt saying it ignores the
image with only color band.

I have the SRTM3 DEM file for my area of interest. Notice it only has 1
color band. Gray

Driver: GTiff/GeoTIFF
Files: hillshade_cut_33N.tif
Size is 13, 17
Coordinate System is:
PROJCS["WGS 84 / UTM zone 33N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",15],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32633"]]
Pixel Size = (75.631606558927956,-75.631606558927956)
Band 1 Block=13x17 Type=Byte, ColorInterp=Gray
  NoData Value=255

My geotiff mosiac:

Driver: GTiff/GeoTIFF
Files: mosaic_33N.tif
 Size is 7946, 6763
Coordinate System is:
PROJCS["WGS 84 / UTM zone 33N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",15],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32633"]]
Pixel Size = (0.026845850153315,-0.026845850153315)
Band 1 Block=7946x1 Type=Byte, ColorInterp=Red
  NoData Value=255
Band 2 Block=7946x1 Type=Byte, ColorInterp=Green
  NoData Value=255
Band 3 Block=7946x1 Type=Byte, ColorInterp=Blue
  NoData Value=255

I would like to drape the mosaic over top of the DEM background tiff, which
is slightly larger. Is it possible to mosaic these images together using
gdalbuildvrt? I realize the pixel size of the DEM image is much larger, my
goal is to generate the same sort of image as you would see viewing both of
these tifs together in ArcGIS or QGIS. Hope this makes sense.

Many thanks,
Jon-Pierre
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20120327/6f833507/attachment.html