<div dir="ltr">That is correct, my image is a Mars southern pole image. Those radii are pulled from another file but I have run the command with the +a/b being in meters, specifically <span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:14.6667px">3396000 and I still get the same result (missing polar information).</span><div><font color="#000000" face="Calibri, Arial, Helvetica, sans-serif"><span style="font-size:14.6667px"><br></span></font></div><div><font color="#000000" face="Calibri, Arial, Helvetica, sans-serif"><span style="font-size:14.6667px">I have also run the command using -wo SAMPLE_GRID=YES and nothing changed. I just tried the command with the -wo SAMPLE_GRID and the +a/b being 337500 and I still get the pole clipped.</span></font></div><div><font color="#000000" face="Calibri, Arial, Helvetica, sans-serif"><span style="font-size:14.6667px"><br></span></font></div><div><font color="#000000" face="Calibri, Arial, Helvetica, sans-serif"><span style="font-size:14.6667px">I would be happy to submit a ticket. Is this the location that I would do that : </span></font>
<a href="https://github.com/OSGeo/gdal">https://github.com/OSGeo/gdal</a> </div><div><br></div><div>Thank you Trent and Even! <font color="#000000" face="Calibri, Arial, Helvetica, sans-serif"><span style="font-size:14.6667px"><br></span></font><div><span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:14.6667px"><br></span></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, May 8, 2020 at 11:55 AM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><u></u>
<div style="font-family:monospace;font-size:9pt;font-weight:400;font-style:normal">
<p style="margin:0px;text-indent:0px">Meghan,</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">Handling poles tend to be fragile, so perhaps something that worked previously no longer works.</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">Did you try -wo SAMPLE_GRID=YES ?</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">If that still doesn't work, the best is you open a ticket with an input image attached (dimensions reduced as much as possible while still reproducing the issue), and the exact command line to reproduce.</p>
<p style="margin:0px;text-indent:0px">I'd note that in your below attempt there's a likely mismatch in the units of the target ellipsoid w.r.t the initial one. This should be a/b=3375000.</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">Even</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">> Hello!</p>
<p style="margin:0px;text-indent:0px">> </p>
<p style="margin:0px;text-indent:0px">> I'm hoping someone can shine some light in my direction because I can't</p>
<p style="margin:0px;text-indent:0px">> find any answers.</p>
<p style="margin:0px;text-indent:0px">> </p>
<p style="margin:0px;text-indent:0px">> I'm trying to warp a Polar Stereographic image into an Equal Cylindrical</p>
<p style="margin:0px;text-indent:0px">> projection. When I run the gdalwarp command on older versions (specifically</p>
<p style="margin:0px;text-indent:0px">> 1.8.1), my image looks as I would expect. When I run the same command in</p>
<p style="margin:0px;text-indent:0px">> gdal versions 2.2.4 or 2.4.1, the pole is clipped out and missing. I've</p>
<p style="margin:0px;text-indent:0px">> tried adding different parameters such as -te, and multiple -wo parameters</p>
<p style="margin:0px;text-indent:0px">> and nothing works for me. Am I missing something? Any help would be greatly</p>
<p style="margin:0px;text-indent:0px">> appreciated because I can't use the older versions anymore in my script. My</p>
<p style="margin:0px;text-indent:0px">> original warp command is below as well as the original polar header.</p>
<p style="margin:0px;text-indent:0px">> Although, I have tried multiple images and I still don't get the results</p>
<p style="margin:0px;text-indent:0px">> that I was getting in version 1.8.</p>
<p style="margin:0px;text-indent:0px">> </p>
<p style="margin:0px;text-indent:0px">> Thank you in advance!!</p>
<p style="margin:0px;text-indent:0px">> </p>
<p style="margin:0px;text-indent:0px">> Meghan</p>
<p style="margin:0px;text-indent:0px">> </p>
<p style="margin:0px;text-indent:0px">> gdalwarp -overwrite -of GTiff -t_srs "+proj=eqc +a=3375.0 +b=3375.0</p>
<p style="margin:0px;text-indent:0px">> +towgs=0,0,0,0,0,0" input.tif output.tif</p>
<p style="margin:0px;text-indent:0px">> </p>
<p style="margin:0px;text-indent:0px">> Driver: GTiff/GeoTIFF</p>
<p style="margin:0px;text-indent:0px">> Files: SP_MY34_summer_18mpp_gd2.tif</p>
<p style="margin:0px;text-indent:0px">> Size is 35928, 35815</p>
<p style="margin:0px;text-indent:0px">> Coordinate System is:</p>
<p style="margin:0px;text-indent:0px">> PROJCS["unnamed",</p>
<p style="margin:0px;text-indent:0px">> GEOGCS["unnamed ellipse",</p>
<p style="margin:0px;text-indent:0px">> DATUM["unknown",</p>
<p style="margin:0px;text-indent:0px">> SPHEROID["unnamed",3396000,0]],</p>
<p style="margin:0px;text-indent:0px">> PRIMEM["Greenwich",0],</p>
<p style="margin:0px;text-indent:0px">> UNIT["degree",0.0174532925199433]],</p>
<p style="margin:0px;text-indent:0px">> PROJECTION["Polar_Stereographic"],</p>
<p style="margin:0px;text-indent:0px">> PARAMETER["latitude_of_origin",-90],</p>
<p style="margin:0px;text-indent:0px">> PARAMETER["central_meridian",0],</p>
<p style="margin:0px;text-indent:0px">> PARAMETER["scale_factor",1],</p>
<p style="margin:0px;text-indent:0px">> PARAMETER["false_easting",0],</p>
<p style="margin:0px;text-indent:0px">> PARAMETER["false_northing",0],</p>
<p style="margin:0px;text-indent:0px">> UNIT["metre",1,</p>
<p style="margin:0px;text-indent:0px">> AUTHORITY["EPSG","9001"]]]</p>
<p style="margin:0px;text-indent:0px">> Origin = (-383491.061348880291916,410399.872797137824818)</p>
<p style="margin:0px;text-indent:0px">> Pixel Size = (17.999204982112094,-17.999204982112094)</p>
<p style="margin:0px;text-indent:0px">> Metadata:</p>
<p style="margin:0px;text-indent:0px">> AREA_OR_POINT=Area</p>
<p style="margin:0px;text-indent:0px">> Image Structure Metadata:</p>
<p style="margin:0px;text-indent:0px">> COMPRESSION=DEFLATE</p>
<p style="margin:0px;text-indent:0px">> INTERLEAVE=BAND</p>
<p style="margin:0px;text-indent:0px">> Corner Coordinates:</p>
<p style="margin:0px;text-indent:0px">> Upper Left ( -383491.061, 410399.873) ( 43d 3'31.36"W, 80d32'41.89"S)</p>
<p style="margin:0px;text-indent:0px">> Lower Left ( -383491.061, -234241.654) (121d25' 1.79"W, 82d25'46.00"S)</p>
<p style="margin:0px;text-indent:0px">> Upper Right ( 263184.375, 410399.873) ( 32d40'17.59"E, 81d47'18.78"S)</p>
<p style="margin:0px;text-indent:0px">> Lower Right ( 263184.375, -234241.654) (131d40'12.00"E, 84d 3'39.60"S)</p>
<p style="margin:0px;text-indent:0px">> Center ( -60153.343, 88079.110) ( 34d19'51.59"W, 88d12' 2.26"S)</p>
<p style="margin:0px;text-indent:0px">> Band 1 Block=256x256 Type=Byte, ColorInterp=Gray</p>
<p style="margin:0px;text-indent:0px">> Min=2146.000 Max=10112.000</p>
<p style="margin:0px;text-indent:0px">> Minimum=2146.000, Maximum=10112.000, Mean=5382.365, StdDev=1638.523</p>
<p style="margin:0px;text-indent:0px">> NoData Value=0</p>
<p style="margin:0px;text-indent:0px">> Overviews: 17964x17908, 8982x8954, 4491x4477, 2246x2239, 1123x1120,</p>
<p style="margin:0px;text-indent:0px">> 562x560, 281x280, 141x140, 71x70</p>
<p style="margin:0px;text-indent:0px">> Metadata:</p>
<p style="margin:0px;text-indent:0px">> STATISTICS_MAXIMUM=10112</p>
<p style="margin:0px;text-indent:0px">> STATISTICS_MEAN=5382.3651107267</p>
<p style="margin:0px;text-indent:0px">> STATISTICS_MINIMUM=2146</p>
<p style="margin:0px;text-indent:0px">> STATISTICS_STDDEV=1638.5226548658</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">-- </p>
<p style="margin:0px;text-indent:0px">Spatialys - Geospatial professional services</p>
<p style="margin:0px;text-indent:0px"><a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a></p></div></blockquote></div>