<div dir="ltr"><div><br></div>This is complicated. I could be wrong here, but what I think would be happening is that after you crop, you have a rhombus in your original CRS. Then when you run the GDAL writer, it makes a *rectangular* tiff (because that's what it does, even if you're missing data in some places. Then when you transform that rectangle in your CRS to 4326, you're back to something non-rectangular.<br><div><br></div><div>Maybe someone else can chime in, but I would think the easiest thing to do would be to transform your points to 4326 right after you read them since that's your ultimate goal, I think.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Aug 17, 2021 at 1:45 PM Andreas Yankopolus <<a href="mailto:andreas@yank.to">andreas@yank.to</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">Andrew,<br>
<br>
> If the SRS of the bounding region (you specified 4326) doesn't match that of the input, it will be transformed to match that of the points before cropping occurs, which probably means that what you specified as a box becomes something else (rhombus?).<br>
<br>
I’m tracking here—the square in WGS84 is not a square in the LPC SRS—but shouldn’t the resulting cropped region become a square when read in by GDAL and transformed into WGS84?<br>
<br>
Thanks,<br>
<br>
Andreas<br>
<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature">Andrew Bell<br><a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a></div>