<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Thanks for the data.<div class=""><br class=""></div><div class="">PDAL filters.crop doesn't first buffer the point and then do point in polygon on that buffered polygon. It instead computes the circular (2d) or spherical (3d) distance. See <a href="https://github.com/PDAL/PDAL/blob/master/filters/CropFilter.cpp#L255" class="">https://github.com/PDAL/PDAL/blob/master/filters/CropFilter.cpp#L255</a> for more information on what's done. </div><div class=""><br class=""></div><div class="">If you want to use a buffer, you'll need to pass one explicitly. </div><div class=""><br class=""></div><div class="">Presumably, this is the reason for the discrepancy. I will update our docs to be more explicit about this behavior.</div><div class=""><br class=""></div><div class="">Howard</div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Oct 3, 2017, at 8:40 AM, Howard Butler <<a href="mailto:howard@hobu.co" class="">howard@hobu.co</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html charset=us-ascii" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Can you post example data? Presumably there is an issue with how the buffer is generated from the point.<div class=""><br class=""></div><div class="">Howard</div><div class=""><br class=""><div class=""><blockquote type="cite" class=""><div class="">On Oct 3, 2017, at 8:38 AM, Giuseppe Falcone <<a href="mailto:falcone.giuseppe@gmail.com" class="">falcone.giuseppe@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class="">Hi to all,<div class=""><br class=""></div><div class="">I have a problem with a simple pdal pipeline.</div><div class="">I have to filter point from las file that are into a 13m radius buffer from given point.</div><div class=""><br class=""></div><div class=""><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;color:rgb(51,51,51);padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><span class="gmail-n" style="box-sizing:border-box">json</span> <span class="gmail-o" style="box-sizing:border-box;color:rgb(102,102,102)">=</span> <span class="gmail-sa" style="box-sizing:border-box">u</span><span class="gmail-s1" style="box-sizing:border-box;color:rgb(186,33,33)">'''</span>
<span class="gmail-s1" style="box-sizing:border-box;color:rgb(186,33,33)">{</span>
<span class="gmail-s1" style="box-sizing:border-box;color:rgb(186,33,33)"> "pipeline":[</span>
<span class="gmail-s1" style="box-sizing:border-box;color:rgb(186,33,33)"> "data.las",</span></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""> {
"distance": 13,
"type": "filters.crop",
"point": "POINT(626715,4481770,155)"
}
<span class="gmail-s1" style="color:rgb(186,33,33);box-sizing:border-box">}''</span></pre></div><div class=""><br class=""></div><div class=""><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><span style="color:rgb(51,51,51)" class=""><span class="gmail-n" style="box-sizing:border-box">p</span> <span class="gmail-o" style="box-sizing:border-box;color:rgb(102,102,102)">=</span> <span class="gmail-n" style="box-sizing:border-box">pdal</span><span class="gmail-o" style="box-sizing:border-box;color:rgb(102,102,102)">.</span><span class="gmail-n" style="box-sizing:border-box">Pipeline</span><span class="gmail-p" style="box-sizing:border-box">(</span><span class="gmail-n" style="box-sizing:border-box">json</span><span class="gmail-p" style="box-sizing:border-box">)</span>
<span class="gmail-n" style="box-sizing:border-box">p</span><span class="gmail-o" style="box-sizing:border-box;color:rgb(102,102,102)">.</span><span class="gmail-n" style="box-sizing:border-box">validate</span><span class="gmail-p" style="box-sizing:border-box">()</span>
<span class="gmail-n" style="box-sizing:border-box">p</span><span class="gmail-o" style="box-sizing:border-box;color:rgb(102,102,102)">.</span><span class="gmail-n" style="box-sizing:border-box">loglevel</span> <span class="gmail-o" style="box-sizing:border-box;color:rgb(102,102,102)">=</span> </span><font color="#666666" class="">4</font><font color="#333333" class="">
</font><span class="gmail-n" style="color:rgb(51,51,51);box-sizing:border-box">count</span><font color="#333333" class=""> </font><span class="gmail-o" style="color:rgb(102,102,102);box-sizing:border-box">=</span><font color="#333333" class=""> </font><span class="gmail-n" style="color:rgb(51,51,51);box-sizing:border-box">p</span><span class="gmail-o" style="color:rgb(102,102,102);box-sizing:border-box">.</span><span class="gmail-n" style="color:rgb(51,51,51);box-sizing:border-box">execute</span><span class="gmail-p" style="color:rgb(51,51,51);box-sizing:border-box">()</span><font color="#333333" class="">
</font><span class="gmail-n" style="color:rgb(51,51,51);box-sizing:border-box">data</span><font color="#333333" class=""> </font><span class="gmail-o" style="color:rgb(102,102,102);box-sizing:border-box">=</span><font color="#333333" class=""> </font><span class="gmail-n" style="color:rgb(51,51,51);box-sizing:border-box">p</span><span class="gmail-o" style="color:rgb(102,102,102);box-sizing:border-box">.</span><span class="gmail-n" style="color:rgb(51,51,51);box-sizing:border-box">arrays</span><span class="gmail-p" style="color:rgb(51,51,51);box-sizing:border-box">[</span><span class="gmail-mi" style="color:rgb(102,102,102);box-sizing:border-box">0</span><span class="gmail-p" style="color:rgb(51,51,51);box-sizing:border-box">]</span></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><span class="gmail-p" style="color:rgb(51,51,51);box-sizing:border-box"><br class=""></span></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><span class="gmail-p" style="color:rgb(51,51,51);box-sizing:border-box">data is an array of 5321 element.</span></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><span class="gmail-p" style="color:rgb(51,51,51);box-sizing:border-box"><br class=""></span></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><span class="gmail-p" style="color:rgb(51,51,51);box-sizing:border-box">If I elaborate data with arcmap (view las data a clip with circle),</span></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><span class="gmail-p" style="color:rgb(51,51,51);box-sizing:border-box">I have 6203 element into point buffer.</span></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><span class="gmail-p" style="color:rgb(51,51,51);box-sizing:border-box">I rewrite same code with R and I have also 6203 element as intersection result.</span></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><br class=""></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><font color="#333333" class="">Where is the error with pdal pipeline?</font></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><font color="#333333" class="">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</font></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><font color="#333333" class=""><br class=""></font></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><font color="#333333" class="">Thanks.</font></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><font color="#333333" class=""><br class=""></font></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><font color="#333333" class="">Giuseppe</font></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><span class="gmail-p" style="color:rgb(51,51,51);box-sizing:border-box"><br class=""></span></pre><pre style="box-sizing:border-box;overflow:auto;font-size:14px;line-height:inherit;padding:0px;margin-top:0px;margin-bottom:0px;word-break:break-all;word-wrap:break-word;background-color:rgb(247,247,247);border:none;border-radius:2px" class=""><span class="gmail-p" style="color:rgb(51,51,51);box-sizing:border-box"><br class=""></span></pre></div></div>
_______________________________________________<br class="">pdal mailing list<br class=""><a href="mailto:pdal@lists.osgeo.org" class="">pdal@lists.osgeo.org</a><br class=""><a href="https://lists.osgeo.org/mailman/listinfo/pdal" class="">https://lists.osgeo.org/mailman/listinfo/pdal</a></div></blockquote></div><br class=""></div></div></div></blockquote></div><br class=""></div></body></html>