<br><br><div class="gmail_quote">On Thu, Aug 16, 2012 at 4:32 PM, Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@mines-paris.org" target="_blank">even.rouault@mines-paris.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Selon Margherita Di Leo <<a href="mailto:diregola@gmail.com">diregola@gmail.com</a>>:<br>
<div class="im"><br>
> Hi All,<br>
><br>
> On Mon, Jul 30, 2012 at 7:31 PM, Casper Børgesen (CABO) <<a href="mailto:CABO@niras.dk">CABO@niras.dk</a>>wrote:<br>
><br>
</div>> > Hi! ****<br>
> ><br>
> > ** **<br>
<div class="im">> ><br>
> > I’m trying to create a VRT file from 4 GeoTIFFs using gdalbuildvrt. Three<br>
> > of the GeoTIFFs are created using OrthoVista and the last one created with<br>
> > an earlier version of the same program. When I run gdalbuildvrt I get the<br>
> > following error message when I try to add the file created by the early<br>
</div>> > version: ****<br>
> ><br>
> > ** **<br>
<div class="im">> ><br>
> > Warning 6: gdalbuildvrt does not support heterogenous band<br>
</div>> > characteristics. Skipping...****<br>
> ><br>
> > ** **<br>
<div class="im">> ><br>
> > I get the same error message, so may be related... I've tried to make a<br>
> virtual mosaic, using SRTM DEMs,<br>
><br>
> gdalbuildvrt -input_file_list list.txt STM_CP-DEMS_EU_mosaic.vrt<br>
><br>
> and got:<br>
><br>
> 0...Warning 6: gdalbuildvrt does not support heterogenous band<br>
> characteristics. Skipping STM_CP-DEMS_2500015000.tif<br>
> 10.Warning 6: gdalbuildvrt does not support heterogenous band<br>
> characteristics. Skipping STM_CP-DEMS_2500025000.tif<br>
> ..Warning 6: gdalbuildvrt does not support heterogenous band<br>
> characteristics. Skipping STM_CP-DEMS_2500035000.tif<br>
> [....]<br>
><br>
> Hence, I made my mosaic using:<br>
> gdalwarp --optfile list.txt STM_CP-DEMS_EU_mosaic.tif<br>
> It took a bit more time obviously, but it didn't complain.<br>
><br>
> So.. Could it be a bug of gdalbuildvrt, or.. ?<br>
<br>
</div>I don't think it is a bug. gdalbuildvrt is a "lightweight" but restrictive<br>
mosaicing solution, whereas gdalwarp is "heavier" but more capable.<br>
<br>
gdalbuildvrt makes various checks to ensure that the tiles to be mosaiced are<br>
very similar. Perhaps some of the tests are too restrictive. You should try to<br>
compare the output of gdalinfo between a file that is integrated into the VRT<br>
and one that is skipped to see which differences there are (band data type,<br>
presence/absence of color table, color interpretation, etc...)<br></blockquote><div><br>Even, thank you for replying. Actually the files are all skipped, except the first one of course. About this file, gdalinfo says:<br>
<br>gdalinfo STM_CP-DEMS_1790010475.tif<br>Driver: GTiff/GeoTIFF<br>Files: STM_CP-DEMS_1790010475.tif<br>Size is 5000, 2250<br>Coordinate System is:<br>PROJCS["ETRS89 / LAEA Europe",<br> GEOGCS["ETRS89",<br>
DATUM["European_Terrestrial_Reference_System_1989",<br> SPHEROID["GRS 1980",6378137,298.2572221000027,<br> AUTHORITY["EPSG","7019"]],<br> AUTHORITY["EPSG","6258"]],<br>
PRIMEM["Greenwich",0],<br> UNIT["degree",0.0174532925199433],<br> AUTHORITY["EPSG","4258"]],<br> PROJECTION["Lambert_Azimuthal_Equal_Area"],<br> PARAMETER["latitude_of_center",52],<br>
PARAMETER["longitude_of_center",10],<br> PARAMETER["false_easting",4321000],<br> PARAMETER["false_northing",3210000],<br> UNIT["metre",1,<br> AUTHORITY["EPSG","9001"]],<br>
AUTHORITY["EPSG","3035"]]<br>Origin = (1540000.000000000000000,1160000.000000000000000)<br>Pixel Size = (100.000000000000000,-100.000000000000000)<br>Metadata:<br> AREA_OR_POINT=Area<br> TIFFTAG_DATETIME=2010:08:09 09:37:53<br>
TIFFTAG_DOCUMENTNAME=STM_CP-DEMS_1790010475.tif<br> TIFFTAG_IMAGEDESCRIPTION=File written by egcs_wrgtif 2.1<br> TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)<br> TIFFTAG_SOFTWARE=IDL 7.1.1, ITT Visual Information Solutions<br>
TIFFTAG_XRESOLUTION=100<br> TIFFTAG_YRESOLUTION=100<br>Image Structure Metadata:<br> COMPRESSION=LZW<br> INTERLEAVE=BAND<br>Corner Coordinates:<br>Upper Left ( 1540000.000, 1160000.000) ( 18d35'52.39"W, 28d46'34.65"N)<br>
Lower Left ( 1540000.000, 935000.000) ( 17d57'14.30"W, 26d50'46.78"N)<br>Upper Right ( 2040000.000, 1160000.000) ( 13d41'27.23"W, 30d17'11.28"N)<br>Lower Right ( 2040000.000, 935000.000) ( 13d 8'11.03"W, 28d19' 3.39"N)<br>
Center ( 1790000.000, 1047500.000) ( 15d51'19.76"W, 28d35'33.72"N)<br>Band 1 Block=5000x1 Type=Int16, ColorInterp=Gray<br><br>So I take a skipped file for comparison:<br><br>gdalinfo STM_CP-DEMS_2500015000.tif<br>
Driver: GTiff/GeoTIFF<br>Files: STM_CP-DEMS_2500015000.tif<br>Size is 10000, 10000<br>Coordinate System is:<br>PROJCS["ETRS89 / LAEA Europe",<br> GEOGCS["ETRS89",<br> DATUM["European_Terrestrial_Reference_System_1989",<br>
SPHEROID["GRS 1980",6378137,298.2572221000027,<br> AUTHORITY["EPSG","7019"]],<br> AUTHORITY["EPSG","6258"]],<br> PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br> AUTHORITY["EPSG","4258"]],<br> PROJECTION["Lambert_Azimuthal_Equal_Area"],<br> PARAMETER["latitude_of_center",52],<br>
PARAMETER["longitude_of_center",10],<br> PARAMETER["false_easting",4321000],<br> PARAMETER["false_northing",3210000],<br> UNIT["metre",1,<br> AUTHORITY["EPSG","9001"]],<br>
AUTHORITY["EPSG","3035"]]<br>Origin = (2000000.000000000000000,2000000.000000000000000)<br>Pixel Size = (100.000000000000000,-100.000000000000000)<br>Metadata:<br> AREA_OR_POINT=Area<br> TIFFTAG_DATETIME=2010:08:06 10:58:28<br>
TIFFTAG_DOCUMENTNAME=STM_CP-DEMS_2500015000.tif<br> TIFFTAG_IMAGEDESCRIPTION=File written by egcs_wrgtif 2.1<br> TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)<br> TIFFTAG_SOFTWARE=IDL 7.1.1, ITT Visual Information Solutions<br>
TIFFTAG_XRESOLUTION=100<br> TIFFTAG_YRESOLUTION=100<br>Image Structure Metadata:<br> COMPRESSION=LZW<br> INTERLEAVE=BAND<br>Corner Coordinates:<br>Upper Left ( 2000000.000, 2000000.000) ( 16d36'13.25"W, 37d23'21.20"N)<br>
Lower Left ( 2000000.000, 1000000.000) ( 13d41' 3.88"W, 28d46'45.49"N)<br>Upper Right ( 3000000.000, 2000000.000) ( 5d28'46.07"W, 39d52'33.70"N)<br>Lower Right ( 3000000.000, 1000000.000) ( 3d40'26.43"W, 30d56'55.67"N)<br>
Center ( 2500000.000, 1500000.000) ( 9d49'34.83"W, 34d25'40.17"N)<br>Band 1 Block=10000x1 Type=Float32, ColorInterp=Gray<br><br>I can't see anything obvious why these two files can't be mosaiced.. <br>
<br>Thanks,<br></div><br><br><br></div>-- <br>Dr. Margherita Di Leo<br>