<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">All—<div class=""><br class=""></div><div class="">I’m cropping and processing USGS LPC data into WGS84 raster files, but I can’t get the cropping to work out correctly: My pdal crop stage specifies a square region in WGS8, but the resulting files display as rhombuses. Here’s a minimal example:</div><div class=""><br class=""></div><div class=""><font face="Menlo" class=""><span style="font-style: normal;" class="">{<br class=""> "pipeline": [<br class=""> "/home/ayank/Documents/geodata/lpc/*.laz",<br class=""> {<br class=""> "type": "filters.crop",<br class=""> "bounds": "([-117.18,-117.13],[32.68,32.73])",<br class=""> "a_srs": "EPSG:4326"<br class=""> },<br class=""> {<br class=""> "type": "filters.range",<br class=""> "limits": "Classification[2:2]"<br class=""> },<br class=""> {<br class=""> "filename": "/home/ayank/Documents/geodata/export/sd_coord_test.laz",<br class=""> "type": "writers.las",<br class=""> "compression": "laszip",<br class=""> "extra_dims": "all"<br class=""> },<br class=""> {<br class=""> "filename": "/home/ayank/Documents/geodata/export/sd_coord_test.tiff",<br class=""> "type": "writers.gdal",<br class=""> "gdaldriver": "GTiff",<br class=""> "dimension": "Z",<br class=""> "output_type": "idw",<br class=""> "resolution": "10"<br class=""> }</span></font></div><div class=""><font face="Menlo" class=""><span style="font-style: normal;" class=""> ]</span></font></div><div class=""><font face="Menlo" class=""><span style="font-style: normal;" class="">}</span></font></div><div class=""><br class=""></div><div class="">My expectation is that the output files, when read into QGIS with WGS84 (EPSG:4326) as the project CRS, would render as squares. But they’re both a rhombus.</div><div class=""><br class=""></div><div class="">I’ve also tried adding these two stages to the pipeline:</div><div class=""><br class=""></div><div class=""><font face="Menlo" class=""><span style="font-style: normal;" class=""> {<br class=""> "type": "filters.reprojection",<br class=""> "out_srs": "EPSG:4326"<br class=""> },<br class=""><br class=""> {<br class=""> "filename": "/home/ayank/Documents/geodata/export/sd_coord_test_4326.tiff",<br class=""> "type": "writers.gdal",<br class=""> "gdaldriver": "GTiff",<br class=""> "dimension": "Z",<br class=""> "output_type": "idw",<br class=""> "resolution": "9.25925925925925925925925e-5"<br class=""> }</span></font></div><div class=""><br class=""></div><div class="">Running gdalinfo on the resulting TIFF shows it as not a square:</div><div class=""><br class=""></div><div class=""><font face="Menlo" class=""><span style="font-style: normal;" class="">Size is 552, 541<br class=""><br class=""></span></font></div><div class=""><font face="Menlo" class=""><span style="font-style: normal;" class="">Corner Coordinates:<br class="">Upper Left (-117.1805780, 32.7300307) (117d10'50.08"W, 32d43'48.11"N)<br class="">Lower Left (-117.1805780, 32.6799381) (117d10'50.08"W, 32d40'47.78"N)<br class="">Upper Right (-117.1294669, 32.7300307) (117d 7'46.08"W, 32d43'48.11"N)<br class="">Lower Right (-117.1294669, 32.6799381) (117d 7'46.08"W, 32d40'47.78"N)<br class="">Center (-117.1550224, 32.7049844) (117d 9'18.08"W, 32d42'17.94”N)<br class=""></span></font><br class=""></div><div class="">It overlaps with the previous two output files perfectly.</div><div class=""><br class=""></div><div class="">What am I doing wrong?</div><div class=""><br class=""></div><div class="">My thinking has been to keep the LPC data in its original CRS to avoid dealing with scaling factors and offsets. It appears to combine NAD83 and survey feet.</div><div class=""><br class=""></div><div class="">Thanks,</div><div class=""><br class=""></div><div class="">Andreas</div></body></html>