<div dir="ltr">Hi,<div>I've been using pdal info --boundary to generate bounds for datasets. I take the boundary_json element and feed it to ogr2ogr to use for clipping. That has worked well. However, if the data is in geographic, the coordinates get flipped. Instead of xy, the entries are lat lon. For example (in Florida):</div><div><br></div><div>{"type": "MultiPolygon", "coordinates": [[[[26.0972695, -80.41166909], ...</div><div><br></div><div>That seems like it would match with the GDAL expectation since GDAL 3  honors by default authority (EPSG) axis order, so that means lat/lon for geographic CRS. But I think it conflicts with GeoJSON in the spec <a href="https://datatracker.ietf.org/doc/html/rfc7946#section-3.1.1">section 3.1.1</a>. In any case, when I convert that JSON of the bounds with something like:</div><div><br></div><div>ogr2ogr -a_srs EPSG:4759 bounds.gpkg bounds.json</div><div><br></div><div>I'll get the order flipped and it will be south of Africa instead of in Florida. If I flip all the coordinates in the json and then run it through ogr2ogr, it lands in the right place. Is this a PDAL issue or am I just doing something wrong?</div><div><br></div><div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Kirk</div></div></div></div></div>