[pdal] LAS 1.4 & WKT vertical coordinate system and

Andrew Bell andrew.bell.ia at gmail.com
Mon Dec 7 16:36:11 PST 2015


On Fri, Dec 4, 2015 at 3:26 PM, Robert Dzur <rdzur at bhinc.com> wrote:

> Hello,
>
>
>
> 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.
>
>
>
> 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.
>
>
>
> 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"
>
>
>
> 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.
>
>
>
> Thanks in advance for any advice.
>

Thanks for locating my bugs.  Hopefully these fixes address your problems.

https://github.com/PDAL/PDAL/issues/1063
https://github.com/PDAL/PDAL/issues/1064
https://github.com/PDAL/PDAL/issues/1065

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


More information about the pdal mailing list