[pdal] Crop SRS not working as expected

[Andrew Bell]

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.

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.

On Tue, Aug 17, 2021 at 1:45 PM Andreas Yankopolus <andreas at yank.to> wrote:

> Andrew,
>
> > 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?).
>
> 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?
>
> Thanks,
>
> Andreas
>
>

-- 
Andrew Bell
andrew.bell.ia at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20210817/f94427f6/attachment.html>