PCL点云分割：Conditional Euclidean Clustering

背景：

pcl官方教程:http://www.pointclouds.org/documentation/tutorials/conditional_euclidean_clustering.php#conditional-euclidean-clustering

1.使用感受

这个算法实际就是region growing segmentation的扩展版。你可以先看看我的另一个文章，关于region growing segmentation的：点击打开链接。Conditional Euclidean Clustering实际就是放开那个蔓延的condition, 你可以写一个自己的condition函数，而官方教程中提供了以下这几种condition函数：enforceIntensitySimilarity, enforceCurvatureOrIntensitySimilarity 还有customRegionGrowing。而且一些阈值你需要在你的condition函数里设置好，你的condition函数变成一个回调函数setConditionFunction的输入。至于调参的感受和我的那篇region growing segmentation的文章里说的一样，调参起来还是相对容易的。主要需要调节蔓延时候用到的阈值。

如果你对你的场景里的墙壁地板这些物体的点的数目有个先验的信息，你可以在用户setMaxClusterSize()来过滤掉他们。再用setMinClusterSize()来过滤一些噪音。还可以调用getRemoveCluster()来获取这些被过滤掉的点云。

2. 算法细节

为了了解设置这些参数，让我们来了解一下如何内部细节。从官方的介绍里，算法的流程是写的很仔细的。为了让你更直观的明白，我概括的来说，就是从曲率小的面播种，从种子的位置出发，开始往四周搜索点，然后先对比点与点之间的距离是否小于距离阈值，如果是然后再比对点于点之间的法线/Intensity/Color/Curvature等的差距，如果差距小于阈值就视为同一个cluster。如果一个cluster无法再蔓延，在剩下的点云里再播种，然后继续重复以上步骤找新的cluster, 直到最后遍历完毕。

3.调参

略。可以参考我的关于region growing segmentation的文章。

你需要在你自己的condition函数里设置好阈值。例如：

bool
customRegionGrowing (const PointTypeFull& point_a, const PointTypeFull& point_b, float squared_distance)
{
  Eigen::Map<const Eigen::Vector3f> point_a_normal = point_a.getNormalVector3fMap (), point_b_normal = point_b.getNormalVector3fMap ();
  if (squared_distance < 10000)
  {
    if (fabs (point_a.intensity - point_b.intensity) < 8.0f)
      return (true);
    if (fabs (point_a_normal.dot (point_b_normal)) < 0.02)
      return (true);
  }
  else
  {
    if (fabs (point_a.intensity - point_b.intensity) < 3.0f)
      return (true);
  }
  return (false);
}

这里的法线和intensity阈值8.0f, 0.02和3.0f就需要你根据实际调节。那么距离阈值在哪调节？是在main函数里用setClusterTolerance()来设置。那么官方提供的customRegionGrowing里面为什么还加一个

if (squared_distance < 10000)

经过我的实验，我发现condition函数是首先经过setClusterTolerance()设置的距离阈值过滤后才运行的。像customRegionGrowing里那样再加一个距离阈值，其实是为了根据距离来自动调节法线和intensity阈值。

可以参考点击打开链接里面的参数解释。

4.跑数据集

本来我想跑自己的数据集的，但我发现我的数据集少了一个维度的信息：intensity。Intensity数据在我看看来就是，激光雷达发出激光打到物体后反射回来的激光强度。

那就跑跑他们的数据吧，同时在代码最底部加了可视化：

#include <pcl/visualization/cloud_viewer.h> 
  pcl::visualization::CloudViewer viewer ("Cluster viewer");
  viewer.showCloud(cloud_out);
  while (!viewer.wasStopped ())
  {
  }

效果：