[gdal-dev] [gdal 2.2.2] different result with gdal 2.1.0

[Even Rouault]

On mardi 31 octobre 2017 07:47:50 CET Vincent BERNAT wrote:
> Hello all,
> 
> 
> I was working with gdal 2.1.0 and i decide to update to 2.2.2
> 
> I work on Centos.
> 
> 
> Everything is ok but on one image, i don't get the same result (and the
> result with 2.2.2 seems incorrect).
> 
> 
> One of the most important difference is the result of the call of this
> function :
> 
> GDALSuggestedWarpOutput

You don't give enough information for other people to be able to reproduce. Which 
arguments do you pass to GDALSuggestedWarpOutput ?

Can you reproduce with gdalwarp ?

Here's what I've tried:

* generate an image with equivalent size and georeferencing as yours

gdal_translate byte.tif test.tif -a_srs 'PROJCS["ARC_System_Zone_18",
    GEOGCS["GCS_Sphere",
        DATUM["D_Sphere",
            SPHEROID["Sphere",6378137.0,0.0]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Azimuthal_Equidistant"],
    PARAMETER["latitude_of_center",-90],
    PARAMETER["longitude_of_center",0],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0]]' \
   -outsize 1536 1536 -a_ullr 9713.965 556512.681 17400.989 548825.658

* wrap it to WGS84

gdalwarp test.tif out.tif -overwrite -t_srs EPSG:4326


With GDAL trunk or 2.1, I get the exact same result:

$ gdalinfo -checksum out.tif
Driver: GTiff/GeoTIFF
Files: out.tif
Size is 2165, 188
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 = (1.000000002034652,-84.998318292160022)
Pixel Size = (0.000376928999761,-0.000376928999761)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (   1.0000000, -84.9983183) (  1d 0' 0.00"E, 84d59'53.95"S)
Lower Left  (   1.0000000, -85.0691809) (  1d 0' 0.00"E, 85d 4' 9.05"S)
Upper Right (   1.8160513, -84.9983183) (  1d48'57.78"E, 84d59'53.95"S)
Lower Right (   1.8160513, -85.0691809) (  1d48'57.78"E, 85d 4' 9.05"S)
Center      (   1.4080256, -85.0337496) (  1d24'28.89"E, 85d 2' 1.50"S)
Band 1 Block=2165x3 Type=Byte, ColorInterp=Gray
  Checksum=51779

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20171031/18fd2b79/attachment.html>