[gdal-dev] gdalbuildvrt problem with heterogenous band characteristics

Margherita Di Leo diregola at gmail.com
Tue Aug 21 02:35:45 PDT 2012


Dear all,

in regard of this different behavior performed by gdalwarp and gdalbuildvrt
respect to different data type entries, I have some considerations. The
error raised by gdalbuildvrt made me wonder about my data and then with
your help I discovered what was wrong with them. On the other hand,
gdalwarp didn't raise any complaint, and i happily created my mosaic,
causing a loss of data (from float to int) which i was not aware of. In my
opinion it would be good to have a warning message in both gdalbuildvrt and
gdalwarp if input data are not homogeneous in type.

Thanks,


On Fri, Aug 17, 2012 at 8:56 AM, Margherita Di Leo <diregola at gmail.com>wrote:

> Hi,
>
> On Thu, Aug 16, 2012 at 11:02 PM, Eli Adam <eadam at co.lincoln.or.us> wrote:
>
>> > Band 1 Block=5000x1 Type=Int16, ColorInterp=Gray
>> ...
>> > Band 1 Block=10000x1 Type=Float32, ColorInterp=Gray
>> >
>> > I can't see anything obvious why these two files can't be mosaiced..
>> >
>>
>> It looks like one is 16 bit and the other 32.  gdalbuildvrt probably
>> can't determine the type of output to create.
>>
>> What is the output type of gdalwarp on the same files? Did gdalwarp do
>> anything to make the input files match in the output?
>>
>
> This is the output of the mosaic made by gdalwarp:
>
> gdalinfo STM_CP-DEMS_EU_mosaic.tif
> Driver: GTiff/GeoTIFF
> Files: /home/leomarg/DATA/SRTM_DEMs/_ALL/STM_CP-DEMS_EU_mosaic.tif
> Size is 64600, 50650
>
> Coordinate System is:
> PROJCS["ETRS89 / LAEA Europe",
>     GEOGCS["ETRS89",
>         DATUM["European_Terrestrial_Reference_System_1989",
>             SPHEROID["GRS 1980",6378137,298.2572221010002,
>
>                 AUTHORITY["EPSG","7019"]],
>             AUTHORITY["EPSG","6258"]],
>         PRIMEM["Greenwich",0],
>         UNIT["degree",0.0174532925199433],
>         AUTHORITY["EPSG","4258"]],
>     PROJECTION["Lambert_Azimuthal_Equal_Area"],
>     PARAMETER["latitude_of_center",52],
>     PARAMETER["longitude_of_center",10],
>     PARAMETER["false_easting",4321000],
>     PARAMETER["false_northing",3210000],
>     UNIT["metre",1,
>         AUTHORITY["EPSG","9001"]],
>     AUTHORITY["EPSG","3035"]]
> Origin = (1540000.000000000000000,6000000.000000000000000)
>
> Pixel Size = (100.000000000000000,-100.000000000000000)
> Metadata:
>   AREA_OR_POINT=Area
> Image Structure Metadata:
>   INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  ( 1540000.000, 6000000.000) ( 57d31'59.79"W, 63d25'17.15"N)
>
> Lower Left  ( 1540000.000,  935000.000) ( 17d57'14.30"W, 26d50'46.78"N)
> Upper Right ( 8000000.000, 6000000.000) ( 86d 0'31.83"E, 56d26'45.47"N)
> Lower Right ( 8000000.000,  935000.000) ( 46d17'14.56"E, 23d32'41.83"N)
> Center      ( 4770000.000, 3467500.000) ( 16d52'39.75"E, 54d 7'19.36"N)
> Band 1 Block=64600x1 Type=Int16, ColorInterp=Gray
>
> Thanks,
>
>
> --
> Dr. Margherita Di Leo
>



-- 
Dr. Margherita Di Leo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20120821/19039979/attachment-0001.html>


More information about the gdal-dev mailing list