<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);"><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 11pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">Meghan,</span></div>
<div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);"><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 11pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);"> Looks like you are on Mars and in your warp
command your +a and +b appear to be in kilometers -- when they should be meters.</span></div>
<div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);"><span style="display: inline !important; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);"><br>
</span></div>
<div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);"><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 11pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); display: inline !important;">Now I'm not sure
where those radii values are coming from. I guess I would recommend using the same as the input image though.</span></div>
<div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);"><span style="display: inline !important; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
<div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);"><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 11pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); display: inline !important;">+a=</span><span style="display: inline !important; background-color: rgb(255, 255, 255);"><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 11pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); display: inline !important;">3396000 </span></span><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 11pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); display: inline !important;">+b=</span><span style="display: inline !important; background-color: rgb(255, 255, 255);"><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 11pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); display: inline !important;">3396000</span></span><br>
</div>
<div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);"><span style="display: inline !important; background-color: rgb(255, 255, 255);"><span style="display: inline !important; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);"><br>
</span></span></div>
<div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);"><span style="display: inline !important; background-color: rgb(255, 255, 255);"><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 11pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); display: inline !important;">-Trent</span></span></div>
<div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);"></div>
<br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 11pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
<span style="font-family: "Segoe UI", "Segoe UI Web (West European)", "Segoe UI", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); display: inline !important;"><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 11pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
<span style="font-family: "Segoe UI", "Segoe UI Web (West European)", "Segoe UI", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); display: inline !important;"><br>
</span></div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> gdal-dev <gdal-dev-bounces@lists.osgeo.org> on behalf of Meghan Burris <emmburris@gmail.com><br>
<b>Sent:</b> Friday, May 8, 2020 11:36 AM<br>
<b>To:</b> gdal-dev@lists.osgeo.org <gdal-dev@lists.osgeo.org><br>
<b>Subject:</b> [EXTERNAL] [gdal-dev] gdalwarp results differ in different gdal versions?</font>
<div> </div>
</div>
<div>
<div dir="ltr">Hello!
<div><br>
</div>
<div>I'm hoping someone can shine some light in my direction because I can't find any answers. </div>
<div><br>
</div>
<div>I'm trying to warp a Polar Stereographic image into an Equal Cylindrical projection. When I run the gdalwarp command on older versions (specifically 1.8.1), my image looks as I would expect. When I run the same command in gdal versions 2.2.4 or 2.4.1,
the pole is clipped out and missing. I've tried adding different parameters such as -te, and multiple -wo parameters and nothing works for me. Am I missing something? Any help would be greatly appreciated because I can't use the older versions anymore in my
script. My original warp command is below as well as the original polar header. Although, I have tried multiple images and I still don't get the results that I was getting in version 1.8.</div>
<div><br>
</div>
<div>Thank you in advance!!</div>
<div><br>
</div>
<div>Meghan</div>
<div><br>
</div>
<div>gdalwarp -overwrite -of GTiff -t_srs "+proj=eqc +a=3375.0 +b=3375.0 +towgs=0,0,0,0,0,0" input.tif output.tif <br>
</div>
<div><br>
</div>
<div>Driver: GTiff/GeoTIFF<br>
Files: SP_MY34_summer_18mpp_gd2.tif<br>
Size is 35928, 35815<br>
Coordinate System is:<br>
PROJCS["unnamed",<br>
GEOGCS["unnamed ellipse",<br>
DATUM["unknown",<br>
SPHEROID["unnamed",3396000,0]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433]],<br>
PROJECTION["Polar_Stereographic"],<br>
PARAMETER["latitude_of_origin",-90],<br>
PARAMETER["central_meridian",0],<br>
PARAMETER["scale_factor",1],<br>
PARAMETER["false_easting",0],<br>
PARAMETER["false_northing",0],<br>
UNIT["metre",1,<br>
AUTHORITY["EPSG","9001"]]]<br>
Origin = (-383491.061348880291916,410399.872797137824818)<br>
Pixel Size = (17.999204982112094,-17.999204982112094)<br>
Metadata:<br>
AREA_OR_POINT=Area<br>
Image Structure Metadata:<br>
COMPRESSION=DEFLATE<br>
INTERLEAVE=BAND<br>
Corner Coordinates:<br>
Upper Left ( -383491.061, 410399.873) ( 43d 3'31.36"W, 80d32'41.89"S)<br>
Lower Left ( -383491.061, -234241.654) (121d25' 1.79"W, 82d25'46.00"S)<br>
Upper Right ( 263184.375, 410399.873) ( 32d40'17.59"E, 81d47'18.78"S)<br>
Lower Right ( 263184.375, -234241.654) (131d40'12.00"E, 84d 3'39.60"S)<br>
Center ( -60153.343, 88079.110) ( 34d19'51.59"W, 88d12' 2.26"S)<br>
Band 1 Block=256x256 Type=Byte, ColorInterp=Gray<br>
Min=2146.000 Max=10112.000<br>
Minimum=2146.000, Maximum=10112.000, Mean=5382.365, StdDev=1638.523<br>
NoData Value=0<br>
Overviews: 17964x17908, 8982x8954, 4491x4477, 2246x2239, 1123x1120, 562x560, 281x280, 141x140, 71x70<br>
Metadata:<br>
STATISTICS_MAXIMUM=10112<br>
STATISTICS_MEAN=5382.3651107267<br>
STATISTICS_MINIMUM=2146<br>
STATISTICS_STDDEV=1638.5226548658<br>
</div>
<div><br>
</div>
</div>
</div>
</body>
</html>