[pdal] ground filtering

Matthew Hanson mhanson at appliedgeosolutions.com
Fri Feb 6 13:00:06 PST 2015


Hello,

While pdal ground seems to work fine, I'm trying to get it to work with
pdal pcl and PCL Block JSON, so we can test out different parameters.

However, when I run it using the simple example from the online docs:

$ pdal pcl -i input.las -o output.las -p PMFsimple.json -v4

with this JSON*:

{
    "pipeline":
    {
        "name": "ProgressiveMorphologicalFilterExample",
        "filters":
        [
            {
                "name": "ProgressiveMorphologicalFilter",
                "setMaxWindowSize": 200
            }
        ]
    }}


I get an error:
####################################################
NAME:   ProgressiveMorphologicalFilterExample ()
HELP:
AUTHOR:
--------------------------------------------------------------------------------
process tile 0 through the pipeline

   Step 1) ProgressiveMorphologicalFilter

      max window size: 200
      slope: 1.000000
      max distance: 2.500000
      initial distance: 0.150000
      cell size: 1.000000
      base: 2.000000
      exponential: true
      negative: false
      Iteration 0 (height threshold = 0.150000, window size =
3.000000)...ground now has 108129 points
      Iteration 1 (height threshold = 2.150000, window size =
5.000000)...ground now has 79939 points
      Iteration 2 (height threshold = 2.500000, window size =
9.000000)...ground now has 78532 points
      Iteration 3 (height threshold = 2.500000, window size =
17.000000)...ground now has 78243 points
      Iteration 4 (height threshold = 2.500000, window size =
33.000000)...ground now has 77934 points
      Iteration 5 (height threshold = 2.500000, window size =
65.000000)...ground now has 77537 points
      Iteration 6 (height threshold = 2.500000, window size =
129.000000)...ground now has 77127 points
      Iteration 7 (height threshold = 2.500000, window size =
257.000000)...ground now has 76664 points
      1069983 points filtered to 76664 following progressive morphological
filter
tile 0 filtered to 76664 points, adding to output, which now has 76664
points

.simplepfm-test.las: Found invalid value of '0' for point's return number.
simplepfm-test.las: Found invalid value of '0' for point's number of
returns.
.............................................simplepfm-test.las: Found
invalid value of '28' for point's return number.
simplepfm-test.las: Found invalid value of '224' for point's number of
returns.
Caught PDAL exception: Unable to fetch data and convert as requested:
ScanAngleRank:float(-4.72698e+19) -> signed char
###############################################

Which is similar to the error I was getting with pdal ground a few weeks
ago before it was kindly fixed by the devs.   So I'm wondering if perhaps
this JSON is out of date and there is some additional information required?

Thanks in advance!



* there's a syntax error in the JSON provided on the website with the
commas, this one is fixed



Matthew Hanson
Applied GeoSolutions
(603) 659-3363 x91
http://appliedgeosolutions.com
mhanson at appliedgeosolutions.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20150206/5d85e0f6/attachment.html>


More information about the pdal mailing list