[pdal] pdal pipeline error

[Giuseppe Falcone]

Hi to all,

I have a problem with a simple pdal pipeline.
I have to filter point from las file that are into a 13m radius buffer from
given point.

json = u'''{  "pipeline":[    "data.las",

    {
      "distance": 13,
      "type": "filters.crop",
      "point": "POINT(626715,4481770,155)"
    }}''


p = pdal.Pipeline(json)p.validate() p.loglevel = 4count =
p.execute()data = p.arrays[0]


data is an array of 5321 element.


If I elaborate data with arcmap (view las data a clip with circle),

I have 6203 element into point buffer.

I rewrite same code with R and I have also 6203 element as intersection result.


Where is the error with pdal pipeline?

I have compiled pdal 1.5.0 with pcl 1.8.1 on ubuntu 16.04 LTS and
using it with python 2.7.12


Thanks.


Giuseppe
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20171003/66a9a2ff/attachment.html>