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

classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

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

Rose Phillips

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 [hidden email]) 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
[hidden email]
https://lists.osgeo.org/mailman/listinfo/pdal
Reply | Threaded
Open this post in threaded view
|

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

andrew.bell.ia@gmail.com
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 <[hidden email]> 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 [hidden email]) 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
[hidden email]
https://lists.osgeo.org/mailman/listinfo/pdal


--
Andrew Bell
[hidden email]

_______________________________________________
pdal mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/pdal