[pdal] Adding WKT string to LAS1.4 Point Format ID 6 File

[Andrew Bell]

The problem is that your quotes aren't escaped, I think.

If all you're trying to do is replace an existing SRS (or fill one in where
it doesn't exist), you might just try:

pdal translate in.las out.las
--readers.las.spatialreference="EPSG:2193+7839" --writers.las.forward=all


On Mon, Jun 17, 2019 at 6:20 PM Rose Phillips <RPhillips at linz.govt.nz>
wrote:

> Hi guys,
>
>
>
> I have a question about manually adding a WKT string - EPSG:2193 + 7839
> projection information to a LAS 1.4 Point Format ID 6 file.
>
>
>
> I’ve tried these different methods and had the corresponding errors:-
>
>
>
> Method One
>
>
>
> Input (Command Line)
>
> pdal translate --writers.las.a_srs="COMPD_CS["NZGD2000 / New Zealand
> Transverse Mercator 2000 + NZVD2016 height",PROJCS["NZGD2000 / New Zealand
> Transverse Mercator
> 2000",GEOGCS["NZGD2000",DATUM["New_Zealand_Geodetic_Datum_2000",SPHEROID["GRS
> 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6167"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4167"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",173],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",1600000],PARAMETER["false_northing",10000000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","2193"]],VERT_CS["NZVD2016
> height",VERT_DATUM["New Zealand Vertical Datum
> 2016",2005,AUTHORITY["EPSG","1169"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Up",UP],AUTHORITY["EPSG","7839"]]]"
> H:/RoseBackup/other_datasets/Gisbourne/test_wkt/original/CL2_BF44_2018_1000_4619.las
> H:/RoseBackup/other_datasets/Gisbourne/test_wkt/original/test_pdal_WKT.las
> --writers.las.forward="all"
>
>
>
> Output
> PDAL: Cannot determine reader for input file: /
>
>
>
> Method Two (JSON)
>
> pdal pipeline H:/RoseBackup/other_datasets/Gisbourne/wkt_2193_7839.json
> --readers.las.filename=H:/RoseBackup/other_datasets/Gisbourne/test_wkt/original/CL2_BF44_2018_1000_4619.las
> --writers.las.filename=H:/RoseBackup/other_datasets/Gisbourne/test_wkt/original/test_pdal_WKT.las
>
>
>
> {
>
>   "pipeline" : [
>
>     {
>
>                            "type" : "readers.las",
>
>                            "filename" : "input.las"
>
>     },
>
>     {
>
>                            "type" : "writers.las",
>
>                            "minor_version" : 4,
>
>                            "forward" : "all",
>
>                            "a_srs" : "COMPD_CS["NZGD2000 / New Zealand
> Transverse Mercator 2000 + NZVD2016 height",PROJCS["NZGD2000 / New Zealand
> Transverse Mercator
> 2000",GEOGCS["NZGD2000",DATUM["New_Zealand_Geodetic_Datum_2000",SPHEROID["GRS
> 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6167"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4167"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",173],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",1600000],PARAMETER["false_northing",10000000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","2193"]],VERT_CS["NZVD2016
> height",VERT_DATUM["New Zealand Vertical Datum
> 2016",2005,AUTHORITY["EPSG","1169"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Up",UP],AUTHORITY["EPSG","7839"]]]",
>
>                            "filename" : "output.las"
>
>
>
>     }
>
>  ]
>
> }
>
>
>
> Output
>
> PDAL: JSON pipeline: Unable to parse pipeline:
>
> * Line 10, Column 24
>
>   Missing ',' or '}' in object declaration
>
> ----------------------------------------------------------------
>
> Reference:-
>
>
> http://osgeo-org.1560.x6.nabble.com/pdal-LAS-1-4-amp-WKT-vertical-coordinate-system-and-td5240123.html
>
> ----------------------------------------------------------------
>
> I was wondering if you guys could point me in the direction where I’ve
> gone wrong?
>
> Many thanks,
>
> Rose Phillips
>
> ------------------------------
>
> This message contains information, which may be in confidence and may be
> subject to legal privilege. If you are not the intended recipient, you must
> not peruse, use, disseminate, distribute or copy this message. If you have
> received this message in error, please notify us immediately (Phone 0800
> 665 463 or info at linz.govt.nz) and destroy the original message. LINZ
> accepts no responsibility for changes to this email, or for any
> attachments, after its transmission from LINZ. Thank You.
> _______________________________________________
> pdal mailing list
> pdal at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/pdal



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