<div dir="ltr">Here's a quick update on point cloud processing in VTK. We've got an initial base of decent capability going now, including (and I am absolutely delighted about this) a surface extraction capability based on TSDF (truncated signed distance functions).<div><br></div><div><b>Classes now running (in VTK/Remote/vtkPointCloud):</b></div><div><div>vtkEuclideanClusterExtraction.h<br></div><div>vtkExtractHierarchicalBins.h</div><div>vtkExtractPointCloudPiece.h</div><div>vtkExtractPoints.h</div><div>vtkExtractSurface.h</div><div>vtkFitImplicitFunction.h</div><div>vtkHierarchicalBinningFilter.h</div><div>vtkPCACurvatureEstimation.h</div><div>vtkPCANormalEstimation.h</div><div>vtkPointCloudFilter.h</div><div>vtkRadiusOutlierRemoval.h</div><div>vtkSignedDistance.h</div><div>vtkStatisticalOutlierRemoval.h</div><div>vtkVoxelGrid.h</div></div><div><br></div><div><b>Surface extraction (in VTK/Remote/vtkPointCloud):</b></div><div>We've build a pipeline of objects that perform rudimentary surface extraction from a point cloud. There are three major parts:</div><div>1. Generating normals - currently use vtkPCANormalEstimation, including graph traversal to create consistently oriented normals.</div><div>2. Compute the TSDF using vtkSignedDistance. </div><div>3. Extract zero-crossing isocontour to generate surface (vtkExtractSurface).</div><div><br></div><div>Note that each of these parts will eventually be expanded by writing new filters, for example I want to use in #1 Brad King's tensor voting (from his <a href="http://www.kitware.com/publications/item/view/1253" target="_blank">PhD thesis</a>) to compute better normals. For #2, based loosely on the <a href="https://graphics.stanford.edu/papers/volrange/volrange.pdf" target="_blank">Curless & Levoy</a> paper, we can do better by using different interpolation kernels (see PointInterpolation framework below) rather than the ad hoc probability weights they use (to be added soon). Also for #2 we need to support incremental updating. Finally, #3 is based on a modified <a href="https://www.researchgate.net/publication/282975362_Flying_Edges_A_High-Performance_Scalable_Isocontouring_Algorithm" target="_blank">Flying Edges</a> algorithm so it is very fast :-)</div><div><br></div><div><b>Point Interpolation (in VTK/Filters/Points)</b></div><div>We just pushed this in, it looks to be a very powerful, quite general framework for performing data interpolation across point samples. It's threaded, uses fast threaded locators, and has an abstraction for a family of interpolation kernels (e.g., linear, gaussian, shepard, voronoi, hopefully SPH) at some point soon, with more on the way.</div><div><br></div><div>As usual any feedback is appreciated. Make sure you update both VTK and the vtkPointCloud remote module. Please don't be shy, this is happening really fast and I want to catch the inevitable stupid stuff as soon as possible.</div><div><br></div><div>Best,</div><div>W</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 28, 2016 at 9:12 AM, Will Schroeder <span dir="ltr"><<a href="mailto:will.schroeder@kitware.com" target="_blank">will.schroeder@kitware.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>FYI- I have committed an initial set of filters for performing point cloud processing. Any feedback or suggestions are welcome as this is an initial prototype. The work is currently available as a remote module to VTK (vtkPointCloud) via this repository:</div><div><a href="https://gitlab.kitware.com/vtk/point-cloud.git" target="_blank">https://gitlab.kitware.com/vtk/point-cloud.git</a><br></div><div><br></div><div>A couple of notes:</div><div>+ Right now I am using vtkPolyData to represent the point cloud via a vtkPoints instance. There are no vtkVertex, vtkPolyVertex cells created to save on memory.</div><div>+ The classes will process as input any vtkPointSet dataset</div><div>+ There is a general framework for filtering point clouds via the class vtkPointCloudFilter. Besides their filtered cloud output, these filters also have an optional, second output which contains any points removed from the input.</div><div>+ Current filters include vtkRadiusOutlierRemoval, vtkStatisticalOutlierRemoval, vtkExtractPoints (extract points using an implicit function). Some of  these names are inspired by <a href="http://pointclouds.org/" target="_blank">PCL</a> names.</div><div>+ All filters are threaded using vtkSMPTools using a threaded locator (vtkStaticPointLocator) so I believe that this is relatively fast, although I have not done much testing.</div><div>+ I'm using vtkPointGaussianMapper in the tests, a class that Ken wrote that is very fast. </div><div><br></div><div>As usual comments and suggestions are requested. In particular any suggestions for other filters to write are welcome (to round out some of the core functionality). The repository is in flux as I try crazy ideas and try to educate myself, so be forewarned.</div><div><br></div><div>Best,<br>W</div><div><br></div><div><br></div>
</div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div>William J. Schroeder, PhD<br>Kitware, Inc. - Building the World's Technical Computing Software<br>28 Corporate Drive<br>Clifton Park, NY 12065<br><a href="mailto:will.schroeder@kitware.com" target="_blank">will.schroeder@kitware.com</a><br><a href="http://www.kitware.com" target="_blank">http://www.kitware.com</a><br>(518) 881-4902</div></div></div>
</div>