[pdal] Help by reprojection of laz file

[Connor Manning]

Your transformation is working properly, but you are viewing the data in a
viewer which interprets your data as cartesian.  Your XYZ range of
(degrees, degrees, meters) is probably something like (0.0000001,
0.0000001, 100), and since your Z-values are several orders of magnitude
larger than the XY, you get a vertical line when you interpret this as
cartesian.  You can use any non-geographic coordinate system for a better
viewing experience (or find a viewer that understands geographic
coordinates), but there's nothing wrong with your data transformation
itself.

El mié, 31 ene 2024 a la(s) 10:12 a.m., Heiko Steinfeld via pdal (
pdal at lists.osgeo.org) escribió:

> Hello,
>
> i try to reproject from EPSG:28992 to EPSG:4326. This doesn’t work. The
> result is a straight line.
> Original:
> [image: image.png]
>
> Projected:
> [image: image.png]
>
> This is the command i use: pdal translate ./Test/original.laz
> ./Test/projected.laz reprojection
> --filters.reprojection.in_srs="EPSG:28992"
> --filters.reprojection.out_srs="EPSG:4326" --writers.las.scale_x=0.0000001
> --writers.las.scale_y=0.0000001 --writers.las.offset_x="Auto"
> --writers.las.offset_y="Auto"
>
>
>
> The spatial reference of the original file is: PROJCS[\"Amersfoort / RD
> New\",GEOGCS[\"Amersfoort\",DATUM[\"Amersfoort\",SPHEROID[\"Bessel
> 1841\",6377397.155,299.1528128,AUTHORITY[\"EPSG\",\"7004\"]],AUTHORITY[\"EPSG\",\"6289\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],
>
>
> UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4289\"]],PROJECTION[\"Oblique_Stereographic\"],PARAMETER[\"latitude_of_origin\",52.1561605555556],
>
>
> PARAMETER[\"central_meridian\",5.38763888888889],PARAMETER[\"scale_factor\",0.9999079],PARAMETER[\"false_easting\",155000],PARAMETER[\"false_northing\",463000],
>
>
> UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH],AUTHORITY[\"EPSG\",\"28992\"]]
>
>
>
> I use the pdal version: *2.6.0*
>
> But i also tried 2.6.2 with same result.
>
>
>
> When i look at epsg.io the transformation from EPSG:28992 to EPSG:4326
> should be possible.
>
>
>
> I don’t know what i’m doing wrong. I tried with several laz-files but the
> result is the same for all of them.
>
> I hope anybody can give me a hint.
>
>
>
> Kind regards
>
> Heiko
> _______________________________________________
> pdal mailing list
> pdal at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/pdal
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20240131/bb0ef930/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 328695 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20240131/bb0ef930/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 8852 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20240131/bb0ef930/attachment-0003.png>