<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Dec 4, 2015 at 3:26 PM, Robert Dzur <span dir="ltr"><<a href="mailto:rdzur@bhinc.com" target="_blank">rdzur@bhinc.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">





<div lang="EN-US" link="#0563C1" vlink="#954F72">
<div>
<p class="MsoNormal">Hello,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">With pdal translate we are trying to add WKT vertical coordinate system information to a 1.4 LAS point cloud without modifying any other aspects of the point cloud. 
<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">We are utilizing the ‘forward’ option on the writers.las in the following pdal translate command to maintain point cloud properties, however, in the output file it appears that we are getting some modifications that we would like to avoid.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">pdal translate --writers.las.a_srs="PROJCS[\"NAD83(HARN) / New Mexico Central (ftUS)\",GEOGCS[\"NAD83(HARN)\",DATUM[\"NAD83_High_Accuracy_Reference_Network\",SPHEROID[\"GRS 1980\",6378137,298.257222101,AUTHORITY[\"EPSG\",\"7019\"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY[\"EPSG\",\"6152\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4152\"]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"latitude_of_origin\",31],PARAMETER[\"central_meridian\",-106.25],PARAMETER[\"scale_factor\",0.9999],PARAMETER[\"false_easting\",1640416.667],PARAMETER[\"false_northing\",0],UNIT[\"US
 survey foot\",0.3048006096012192,AUTHORITY[\"EPSG\",\"9003\"]],AXIS[\"X\",EAST],AXIS[\"Y\",NORTH],AUTHORITY[\"EPSG\",\"2903\"]]VERTCS[\"North American Vertical Datum of 1988\",VERT_DATUM[\"North American Vertical Datum 1988\",2005,AUTHORITY[\"EPSG\",\"5103\"]],UNIT[\"US
 survey foot\",1.0],AXIS[\"Gravity-related height\",UP],AUTHORITY[\"EPSG\",\"5703\"]]" pc16901815_flagged_fme_combined.las test_pdal_WKT.las --writers.las.forward="all"<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">The input las “pc16901815_flagged_fme_combined.las” has GPS Type Bit Field set to 1 which appears to get reset to 0 on output (as noted in the altered output Global Encoding value).  Also added are a number of VLR (GeoTIFF included)—and
 we’d like to maintain only WKT definitions.  Also noteworthy is that the “overlap” bit flag on all the points appears to get reset on output as well so that the point cloud drops the definitions for those points that were defined with an overlap classification
 flag in the input las file.  <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Thanks in advance for any advice.</p></div></div></blockquote><div><br></div><div>Thanks for locating my bugs.  Hopefully these fixes address your problems.</div><div><br></div><div><a href="https://github.com/PDAL/PDAL/issues/1063">https://github.com/PDAL/PDAL/issues/1063</a><br></div><div><a href="https://github.com/PDAL/PDAL/issues/1064">https://github.com/PDAL/PDAL/issues/1064</a><br></div><div><a href="https://github.com/PDAL/PDAL/issues/1065">https://github.com/PDAL/PDAL/issues/1065</a><br></div><div><br></div></div>-- <br><div class="gmail_signature">Andrew Bell<br><a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a></div>
</div></div>