Hello all,<div>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. <br>
<br></div><div>I have the SRTM3 DEM file for my area of interest. Notice it only has 1 color band. Gray</div>
<div><br></div><div><div>Driver: GTiff/GeoTIFF</div><div>Files: hillshade_cut_33N.tif</div><div>Size is 13, 17</div><div>Coordinate System is:</div><div>PROJCS["WGS 84 / UTM zone 33N",</div><div> GEOGCS["WGS 84",</div>
<div> DATUM["WGS_1984",</div><div> SPHEROID["WGS 84",6378137,298.257223563,</div><div> AUTHORITY["EPSG","7030"]],</div><div> AUTHORITY["EPSG","6326"]],</div>
<div> PRIMEM["Greenwich",0],</div><div> UNIT["degree",0.0174532925199433],</div><div> AUTHORITY["EPSG","4326"]],</div><div> PROJECTION["Transverse_Mercator"],</div>
<div> PARAMETER["latitude_of_origin",0],</div><div> PARAMETER["central_meridian",15],</div><div> PARAMETER["scale_factor",0.9996],</div><div> PARAMETER["false_easting",500000],</div>
<div> PARAMETER["false_northing",0],</div><div> UNIT["metre",1,</div><div> AUTHORITY["EPSG","9001"]],</div><div> AUTHORITY["EPSG","32633"]]</div>
<div>Pixel Size = (75.631606558927956,-75.631606558927956)</div>
<div>Band 1 Block=13x17 Type=Byte, ColorInterp=Gray</div><div> NoData Value=255</div></div><div><br></div><div>My geotiff mosiac:</div><div><br></div><div><div>Driver: GTiff/GeoTIFF</div><div>Files: mosaic_33N.tif</div>
<div>
Size is 7946, 6763</div><div>Coordinate System is:</div><div>PROJCS["WGS 84 / UTM zone 33N",</div><div> GEOGCS["WGS 84",</div><div> DATUM["WGS_1984",</div><div> SPHEROID["WGS 84",6378137,298.257223563,</div>
<div> AUTHORITY["EPSG","7030"]],</div><div> AUTHORITY["EPSG","6326"]],</div><div> PRIMEM["Greenwich",0],</div><div> UNIT["degree",0.0174532925199433],</div>
<div> AUTHORITY["EPSG","4326"]],</div><div> PROJECTION["Transverse_Mercator"],</div><div> PARAMETER["latitude_of_origin",0],</div><div> PARAMETER["central_meridian",15],</div>
<div> PARAMETER["scale_factor",0.9996],</div><div> PARAMETER["false_easting",500000],</div><div> PARAMETER["false_northing",0],</div><div> UNIT["metre",1,</div><div> AUTHORITY["EPSG","9001"]],</div>
<div> AUTHORITY["EPSG","32633"]]</div><div>Pixel Size = (0.026845850153315,-0.026845850153315)</div>
<div>Band 1 Block=7946x1 Type=Byte, ColorInterp=Red</div><div> NoData Value=255</div><div>Band 2 Block=7946x1 Type=Byte, ColorInterp=Green</div><div> NoData Value=255</div><div>Band 3 Block=7946x1 Type=Byte, ColorInterp=Blue</div>
<div> NoData Value=255</div></div><div><br></div><div>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.</div>
<div><br></div><div>Many thanks,</div>
<div>Jon-Pierre</div>