<div dir="ltr">The problem is that your quotes aren't escaped, I think.<div><br></div><div>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:</div><div><br></div><div>pdal translate in.las out.las --readers.las.spatialreference="EPSG:2193+7839" --writers.las.forward=all</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jun 17, 2019 at 6:20 PM Rose Phillips <<a href="mailto:RPhillips@linz.govt.nz">RPhillips@linz.govt.nz</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div lang="EN-NZ">
<div class="gmail-m_-4434244311925071290WordSection1">
<p class="MsoNormal">Hi guys,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">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.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I’ve tried these different methods and had the corresponding errors:-
<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Method One<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Input (Command Line) <u></u><u></u></p>
<p class="MsoNormal">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"<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Output <br>
PDAL: Cannot determine reader for input file: /<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Method Two (JSON) <u></u><u></u></p>
<p class="MsoNormal">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<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">{<u></u><u></u></p>
<p class="MsoNormal"> "pipeline" : [<u></u><u></u></p>
<p class="MsoNormal"> { <u></u><u></u></p>
<p class="MsoNormal"> "type" : "readers.las",<u></u><u></u></p>
<p class="MsoNormal"> "filename" : "input.las"<u></u><u></u></p>
<p class="MsoNormal"> },<u></u><u></u></p>
<p class="MsoNormal"> {<u></u><u></u></p>
<p class="MsoNormal"> "type" : "writers.las",<u></u><u></u></p>
<p class="MsoNormal"> "minor_version" : 4, <u></u><u></u></p>
<p class="MsoNormal"> "forward" : "all",<u></u><u></u></p>
<p class="MsoNormal"> "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"]]]",<u></u><u></u></p>
<p class="MsoNormal"> "filename" : "output.las"<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> }<u></u><u></u></p>
<p class="MsoNormal"> ]<u></u><u></u></p>
<p class="MsoNormal">}<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Output <u></u><u></u></p>
<p class="MsoNormal">PDAL: JSON pipeline: Unable to parse pipeline:<u></u><u></u></p>
<p class="MsoNormal">* Line 10, Column 24<u></u><u></u></p>
<p class="MsoNormal"> Missing ',' or '}' in object declaration<u></u><u></u></p>
<p class="MsoNormal">----------------------------------------------------------------<u></u><u></u></p>
<p class="MsoNormal">Reference:-<u></u><u></u></p>
<p class="MsoNormal"><a href="http://osgeo-org.1560.x6.nabble.com/pdal-LAS-1-4-amp-WKT-vertical-coordinate-system-and-td5240123.html" target="_blank">http://osgeo-org.1560.x6.nabble.com/pdal-LAS-1-4-amp-WKT-vertical-coordinate-system-and-td5240123.html</a><u></u><u></u></p>
<p class="MsoNormal">----------------------------------------------------------------<u></u><u></u></p>
<p class="MsoNormal">I was wondering if you guys could point me in the direction where I’ve gone wrong?<u></u><u></u></p>
<p class="MsoNormal">Many thanks,<u></u><u></u></p>
<p class="MsoNormal">Rose Phillips<u></u><u></u></p>
</div>
<br>
<hr>
<font face="Verdana" color="Black" size="2"><br>
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 <a href="mailto:info@linz.govt.nz" target="_blank">info@linz.govt.nz</a>) 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.<br>
</font>
</div>
_______________________________________________<br>
pdal mailing list<br>
<a href="mailto:pdal@lists.osgeo.org" target="_blank">pdal@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/pdal" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/pdal</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature">Andrew Bell<br><a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a></div>