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 <[email protected]>wrote: > Hi, > > On Thu, Aug 16, 2012 at 11:02 PM, Eli Adam <[email protected]> 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
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
