[gdal-dev] gdalwarp giving two different results with the same input

Fri Oct 24 01:28:17 PDT 2014

Le vendredi 24 octobre 2014 09:35:30, vous avez écrit :
> Hi,
> 
> I am using GDAL (GDAL 1.9.2, released 2012/10/08) to change of EPSG
> referential on a GeoTiff file. I do the exact same command on two
> different computers (with the same version of GDAL and tiff libraries).
> And I get two different results... I am out of options, any idea ? Any
> place I should have a look ?

Pierre,

In case you did several attempts on computer B, did you erase output.tif 
before re-running gdalwarp ? gdalwarp will not do it (this is the mosaicing 
mode onto an existing target), unless you specify -overwrite.

Is the O.S., architecture and binary distribution of GDAL the same on both 
computers ?

Otherwise the other potential suspects might be the GDAL_DATA environment 
variable not defined or pointing to a wrong directory. Or different proj.4 
versions.

Even

> 
> Command :
> 
> gdalwarp -s_srs EPSG:3975 -t_srs EPSG:4326 input.tif output.tif
> 
> Input file :
> 
> Driver: GTiff/GeoTIFF
> Files: input.tif
> Size is 1388, 584
> Coordinate System is `'
> Origin = (-17321655.543400000780821,7332251.365589999593794)
> Pixel Size = (24959.157843515851710,-25110.449882157532556)
> Metadata:
> lat#_FillValue=9.96921e+36
> lat#long_name=latitude
> lat#units=degrees_north
> lon#_FillValue=9.96921e+36
> lon#long_name=longitude
> lon#units=degrees_east
> NC_GLOBAL#conventions=CF-1.4
> NC_GLOBAL#creation_date=2014-10-20T16:12:01
> NC_GLOBAL#datum=+ellps = WGS84
> NC_GLOBAL#ease_global=yes
> NC_GLOBAL#ease_origin_lat=-83.517136
> NC_GLOBAL#ease_origin_lon=-179.87032
> NC_GLOBAL#ease_projection=cylindrical
> NC_GLOBAL#ease_resolution=25
> NC_GLOBAL#grip_mapping=projection
> NC_GLOBAL#history=none
> NC_GLOBAL#institution=SMOS CATDS Processing Chain
> NC_GLOBAL#netcdf_version_id=3.6.2
> NC_GLOBAL#product_version=1.0
> NC_GLOBAL#proj4text=+proj=cea +lon_0=0 +lat_ts=30 +x_0=0 +y_0=0
> +ellps=WGS84 +units=m +no_defs NC_GLOBAL#srid=EPSG:3975
> Image Structure Metadata:
> INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  (-17321655.543, 7332251.366)
> Lower Left  (-17321655.543,-7332251.366)
> Upper Right (17321655.543, 7332251.366)
> Lower Right (17321655.543,-7332251.366)
> Center      (   0.0000000,   0.0000000)
> Band 1 Block=1388x1 Type=Float64, ColorInterp=Gray
> 
> On computer A (the result I am expecting) :
> 
> Driver: GTiff/GeoTIFF
> Files: output.tif
> Size is 1355, 657
> Coordinate System is:
> GEOGCS["WGS 84",
>     DATUM["WGS_1984",
>         SPHEROID["WGS 84",6378137,298.257223563,
>             AUTHORITY["EPSG","7030"]],
>         AUTHORITY["EPSG","6326"]],
>     PRIMEM["Greenwich",0],
>     UNIT["degree",0.0174532925199433],
>     AUTHORITY["EPSG","4326"]]
> Origin = (-179.524544819822239,87.025797325816669)
> Pixel Size = (0.264973369536817,-0.264973369536817)
> Metadata:
>   AREA_OR_POINT=Area
> Image Structure Metadata:
>   INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  (-179.5245448,  87.0257973) (179d31'28.36"W, 87d 1'32.87"N)
> Lower Left  (-179.5245448, -87.0617065) (179d31'28.36"W, 87d 3'42.14"S)
> Upper Right ( 179.5143709,  87.0257973) (179d30'51.74"E, 87d 1'32.87"N)
> Lower Right ( 179.5143709, -87.0617065) (179d30'51.74"E, 87d 3'42.14"S)
> Center      (  -0.0050870,  -0.0179546) (  0d 0'18.31"W,  0d 1' 4.64"S)
> Band 1 Block=1355x1 Type=Float64, ColorInterp=Gray
> 
> On computer B :
> 
> Driver: GTiff/GeoTIFF
> Files: output.tif
> Size is 1358, 651
> Coordinate System is:
> GEOGCS["WGS 84",
>     DATUM["WGS_1984",
>         SPHEROID["WGS 84",6378137,298.257223563,
>             AUTHORITY["EPSG","7030"]],
>         AUTHORITY["EPSG","6326"]],
>     PRIMEM["Greenwich",0],
>     UNIT["degree",0.0174532925199433],
>     AUTHORITY["EPSG","4326"]]
> Origin = (-179.674959335984454,86.219877359659435)
> Pixel Size = (0.264688310439600,-0.264688310439600)
> Metadata:
>   AREA_OR_POINT=Area
> Image Structure Metadata:
>   INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  (-179.6749593,  86.2198774) (179d40'29.85"W, 86d13'11.56"N)
> Lower Left  (-179.6749593, -86.0922127) (179d40'29.85"W, 86d 5'31.97"S)
> Upper Right ( 179.7717662,  86.2198774) (179d46'18.36"E, 86d13'11.56"N)
> Lower Right ( 179.7717662, -86.0922127) (179d46'18.36"E, 86d 5'31.97"S)
> Center      (   0.0484035,   0.0638323) (  0d 2'54.25"E,  0d 3'49.80"N)
> Band 1 Block=1358x1 Type=Float64, ColorInterp=Gray
> 
> Thanks for your help !
> Regards,
> Pierre.
> This message contains information that may be privileged or confidential
> and is the property of the Capgemini Group. It is intended only for the
> person to whom it is addressed. If you are not the intended recipient, you
> are not authorized to read, print, retain, copy, disseminate, distribute,
> or use this message or any part thereof. If you receive this message in
> error, please notify the sender immediately and delete all copies of this
> message.
> 
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com