<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none"><!--P{margin-top:0;margin-bottom:0;} p, li
{white-space:pre-wrap}--></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Thanks for your reply.</p>
<p><br>
</p>
<p>I have these differents only with images using driver SRP/Standard Raster Product (ASRP/USRP)<br>
</p>
<p><br>
</p>
<p>i did it with gdalwrap with 2.1.0 and 2.2.2 :</p>
<p>>>gdalwrap input.IMG output.IMG -overwrite -t_srs EPSG:4326 </p>
<p><br>
</p>
<p>and here are the 2 gdalinfo i got in outputs :</p>
<p><br>
</p>
<p><br>
Driver: GTiff/GeoTIFF<br>
Files: OUT_2.2.2.IMG<br>
Size is 2165, 188<br>
Coordinate System is:<br>
GEOGCS["WGS 84",<br>
DATUM["WGS_1984",<br>
SPHEROID["WGS 84",6378137,298.257223563,<br>
AUTHORITY["EPSG","7030"]],<br>
AUTHORITY["EPSG","6326"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4326"]]<br>
Origin = (1.000000000000000,-84.998318287963713)<br>
Pixel Size = (0.000376928984429,-0.000376928984429)<br>
Metadata:<br>
AREA_OR_POINT=Area<br>
SRP_CREATIONDATE=20070417<br>
SRP_EDN=0<br>
SRP_NAM=APOLES01<br>
SRP_PRODUCT=ASRP<br>
SRP_REVISIONDATE=20070417<br>
SRP_SCA=50000<br>
SRP_ZNA=18<br>
Image Structure Metadata:<br>
INTERLEAVE=BAND<br>
Corner Coordinates:<br>
Upper Left ( 1.0000000, -84.9983183) ( 1d 0' 0.00"E, 84d59'53.95"S)<br>
Lower Left ( 1.0000000, -85.0691809) ( 1d 0' 0.00"E, 85d 4' 9.05"S)<br>
Upper Right ( 1.8160513, -84.9983183) ( 1d48'57.78"E, 84d59'53.95"S)<br>
Lower Right ( 1.8160513, -85.0691809) ( 1d48'57.78"E, 85d 4' 9.05"S)<br>
Center ( 1.4080256, -85.0337496) ( 1d24'28.89"E, 85d 2' 1.50"S)<br>
Band 1 Block=2165x3 Type=Byte, ColorInterp=Palette<br>
NoData Value=0<br>
Color Table (RGB with 256 entries)<br>
0: 0,0,0,0<br>
1: 0,0,17,255<br>
</p>
<p>-------------------------------------------------------------------------------------------------<br>
</p>
<p><br>
</p>
<p>gdalinfo OUT_2.1.0.IMG <br>
Driver: GTiff/GeoTIFF<br>
Size is 11865, 2172<br>
Coordinate System is:<br>
GEOGCS["WGS 84",<br>
DATUM["WGS_1984",<br>
SPHEROID["WGS 84",6378137,298.257223563,<br>
AUTHORITY["EPSG","7030"]],<br>
AUTHORITY["EPSG","6326"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4326"]]<br>
Origin = (-71.565051177077976,-75.410804142370040)<br>
Pixel Size = (0.004477910774178,-0.004477910774178)<br>
Metadata:<br>
AREA_OR_POINT=Area<br>
SRP_CREATIONDATE=20070417<br>
SRP_EDN=0<br>
SRP_NAM=APOLES01<br>
SRP_PRODUCT=ASRP<br>
SRP_REVISIONDATE=20070417<br>
SRP_SCA=50000<br>
SRP_ZNA=18<br>
Image Structure Metadata:<br>
INTERLEAVE=BAND<br>
Corner Coordinates:<br>
Upper Left ( -71.5650512, -75.4108041) ( 71d33'54.18"W, 75d24'38.89"S)<br>
Lower Left ( -71.5650512, -85.1368263) ( 71d33'54.18"W, 85d 8'12.57"S)<br>
Upper Right ( -18.4346398, -75.4108041) ( 18d26' 4.70"W, 75d24'38.89"S)<br>
Lower Right ( -18.4346398, -85.1368263) ( 18d26' 4.70"W, 85d 8'12.57"S)<br>
Center ( -44.9998455, -80.2738152) ( 44d59'59.44"W, 80d16'25.73"S)<br>
Band 1 Block=11865x1 Type=Byte, ColorInterp=Palette<br>
NoData Value=0<br>
Color Table (RGB with 256 entries)<br>
0: 0,0,0,255<br>
</p>
<p>-------------------------------------------------------------------------------------------------<br>
</p>
<p><br>
</p>
<p><br>
</p>
<div style="font-family:'Sans Serif'; font-size:9pt; font-weight:400; font-style:normal">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="Calibri, sans-serif"><b>De :</b> Even Rouault <even.rouault@spatialys.com><br>
<b>Envoyé :</b> mardi 31 octobre 2017 12:49<br>
<b>À :</b> gdal-dev@lists.osgeo.org<br>
<b>Cc :</b> Vincent BERNAT<br>
<b>Objet :</b> Re: [gdal-dev] [gdal 2.2.2] different result with gdal 2.1.0</font>
<div> </div>
</div>
<div>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
On mardi 31 octobre 2017 07:47:50 CET Vincent BERNAT wrote:</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> Hello all,</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> </p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> </p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> I was working with gdal 2.1.0 and i decide to update to 2.2.2</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> </p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> I work on Centos.</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> </p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> </p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> Everything is ok but on one image, i don't get the same result (and the</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> result with 2.2.2 seems incorrect).</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> </p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> </p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> One of the most important difference is the result of the call of this</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> function :</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> </p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
> GDALSuggestedWarpOutput</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
You don't give enough information for other people to be able to reproduce. Which arguments do you pass to GDALSuggestedWarpOutput ?</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Can you reproduce with gdalwarp ?</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Here's what I've tried:</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
* generate an image with equivalent size and georeferencing as yours</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
gdal_translate byte.tif test.tif -a_srs 'PROJCS["ARC_System_Zone_18",</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
GEOGCS["GCS_Sphere",</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
DATUM["D_Sphere",</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
SPHEROID["Sphere",6378137.0,0.0]],</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
PRIMEM["Greenwich",0],</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
UNIT["degree",0.0174532925199433]],</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
PROJECTION["Azimuthal_Equidistant"],</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
PARAMETER["latitude_of_center",-90],</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
PARAMETER["longitude_of_center",0],</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
PARAMETER["false_easting",0],</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
PARAMETER["false_northing",0]]' \</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
-outsize 1536 1536 -a_ullr 9713.965 556512.681 17400.989 548825.658</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
* wrap it to WGS84</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
gdalwarp test.tif out.tif -overwrite -t_srs EPSG:4326</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
With GDAL trunk or 2.1, I get the exact same result:</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
$ gdalinfo -checksum out.tif</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Driver: GTiff/GeoTIFF</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Files: out.tif</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Size is 2165, 188</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Coordinate System is:</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
GEOGCS["WGS 84",</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
DATUM["WGS_1984",</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
SPHEROID["WGS 84",6378137,298.257223563,</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
AUTHORITY["EPSG","7030"]],</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
AUTHORITY["EPSG","6326"]],</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
PRIMEM["Greenwich",0],</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
UNIT["degree",0.0174532925199433],</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
AUTHORITY["EPSG","4326"]]</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Origin = (1.000000002034652,-84.998318292160022)</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Pixel Size = (0.000376928999761,-0.000376928999761)</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Metadata:</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
AREA_OR_POINT=Area</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Image Structure Metadata:</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
INTERLEAVE=BAND</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Corner Coordinates:</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Upper Left ( 1.0000000, -84.9983183) ( 1d 0' 0.00"E, 84d59'53.95"S)</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Lower Left ( 1.0000000, -85.0691809) ( 1d 0' 0.00"E, 85d 4' 9.05"S)</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Upper Right ( 1.8160513, -84.9983183) ( 1d48'57.78"E, 84d59'53.95"S)</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Lower Right ( 1.8160513, -85.0691809) ( 1d48'57.78"E, 85d 4' 9.05"S)</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Center ( 1.4080256, -85.0337496) ( 1d24'28.89"E, 85d 2' 1.50"S)</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Band 1 Block=2165x3 Type=Byte, ColorInterp=Gray</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Checksum=51779</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Even</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
-- </p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Spatialys - Geospatial professional services</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
http://www.spatialys.com</p>
</div>
</div>
</body>
</html>