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
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
