[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