<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=Windows-1252">
<META name=GENERATOR content="MSHTML 8.00.6001.19258"></HEAD>
<BODY>
<DIV><FONT color=#0000ff size=2 face=Arial><SPAN
class=934212315-16082012>Hello,</SPAN></FONT></DIV>
<DIV><FONT color=#0000ff size=2 face=Arial><SPAN
class=934212315-16082012> The problem comes from the data
type: Int16 for the first file, Float32 for the second file.</SPAN></FONT></DIV>
<DIV><FONT color=#0000ff size=2 face=Arial><SPAN
class=934212315-16082012></SPAN></FONT> </DIV>
<DIV><FONT color=#0000ff size=2 face=Arial><SPAN
class=934212315-16082012>Nicolas</SPAN></FONT></DIV>
<DIV dir=ltr class=OutlookMessageHeader align=left><FONT size=2
face=Tahoma>-----Message d'origine-----<BR><B>De :</B>
gdal-dev-bounces@lists.osgeo.org [mailto:gdal-dev-bounces@lists.osgeo.org]<B>De
la part de</B> Margherita Di Leo<BR><B>Envoyé :</B> jeudi 16 août 2012
17:15<BR><B>À :</B> Even Rouault<BR><B>Cc :</B>
gdal-dev@lists.osgeo.org<BR><B>Objet :</B> Re: [gdal-dev] gdalbuildvrt
problem with heterogenous bandcharacteristics<BR><BR></FONT></DIV><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
style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex"
class=gmail_quote>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></BODY></HTML>