<div dir="ltr"><div><div>I'm not sure whether this is a real issue or not but I thought I'd bring it up, at the very least I'll learn something <br><br></div>When using "gdalwarp -cutline <shapefile> -crop_to_cutline" on an input raster that is in EPSG:3857 with square-pixels the output raster, while still EPSG:3857, now has non-square pixels.  <br><br>They're aren't terribly non-square but aren't they supposed to be completely square for EPSG:3857?<br><br></div>(FWIW, I see that there are tickets opened that reference similar issues but they reference the output being shifted or the origin changing, not the pixel shape changing)<br><br><div>For example:<br><div style="margin-left:40px">#Warp our original .tif to EPSG:3857<br></div><div style="margin-left:40px">$gdalwarp \<br>        -t_srs EPSG:3857 \<br>        -dstalpha \<br>        -co TILED=YES \<br>        "ENR_L33.tif" \<br>        "./2.tif"<br><br></div><div style="margin-left:40px">#See that the output pixels are square<br></div><div style="margin-left:40px">$gdalinfo 2.tif<br></div><div style="margin-left:80px">Origin = (-8577554.996301921084523,5421778.172851986251771)<br>Pixel Size = (43.677179501975118,-43.677179501975118)<br><br></div><div style="margin-left:40px"><br></div><div style="margin-left:40px">#Now crop the image to a cutline<br></div><div style="margin-left:40px">$gdalwarp \<br>        -crop_to_cutline \<br>        -dstalpha \<br>        -cutline "./ENR_L33.shp" \<br>        -cblend 10 \<br>        -co TILED=YES \<br>        "./2.tif" \<br>        "./3.tif"<br><br></div><div style="margin-left:40px">#See that output pixels are not square<br></div><div style="margin-left:40px">$ gdalinfo 3.tif <br></div><div style="margin-left:80px">Origin = (-8480047.445924906060100,5366376.137789577245712)<br>Pixel Size = (43.678439570399853,-43.675457269061347)<br><br></div><br>I ran some more tests and without the -crop_to_cutline option the output pixels remain square<br><br></div><div>Thanks,<br></div><div>  Jesse<br></div></div>