PCL点云(平面点云)分割:Plane Model Segmentation

背景:

pcl官方教程:

http://www.pointclouds.org/documentation/tutorials/planar_segmentation.php#planar-segmentation


注:

这个算法并不是真正意义上的分割算法,因为它只能分割出平面的点云。这种情况和PCL地面点云分割算法:Progressive Morphological Filter很像,后者只能分割出地面点云,详见我的另一文章点击打开链接。但是呢,这个算法能够把地面墙壁这些平面给找出来并过滤掉,方便后面的物体的点云分割。


1.使用感受

很简单的一个算法,但是效果很好。只有一个参数(distance threshold)要设置,调参容易。


2.算法细节

它是基于Ransac的来做平面拟合的。首先Ransac是什么东东?我觉得zinnc的这篇文章点击打开链接讲的就很清楚。在这个Plane Model Segmentation算法里,Ransac为了找到点云的平面,不停的改变平面模型(ax+by+cz+d=0)的参数:a, b, c和d。经过多次调参后,找出哪一组的参数能使得这个模型一定程度上拟合最多的点。这个程度就是由distance threshold这个参数来设置。那找到这组参数后,这些能够被拟合的点就是平面的点。


3.调参

把distance threshold调大,离平面更远的点也被算进平面来。distance threshold 可以等同于平面厚度阈值。


4.跑自己的数据集

跑了两次,把墙壁和地板给过滤掉了。


5. 总结

这是一个很好用很容易调参的算法,但它只能找平面的点云。如果拿来过滤墙壁和地板就很好用。

  • 3
    点赞
  • 38
    收藏
    觉得还不错? 一键收藏
  • 6
    评论
以下是一个基于PCL库的平面分割示例代码: ```c++ #include <iostream> #include <pcl/io/pcd_io.h> #include <pcl/point_types.h> #include <pcl/sample_consensus/method_types.h> #include <pcl/sample_consensus/model_types.h> #include <pcl/segmentation/sac_segmentation.h> int main(int argc, char** argv) { // 加载点云数据 pcl::PointCloud<pcl::PointXYZ>::Ptr cloud(new pcl::PointCloud<pcl::PointXYZ>); pcl::io::loadPCDFile<pcl::PointXYZ>("plane.pcd", *cloud); // 创建分割pcl::ModelCoefficients::Ptr coefficients(new pcl::ModelCoefficients); pcl::PointIndices::Ptr inliers(new pcl::PointIndices); pcl::SACSegmentation<pcl::PointXYZ> seg; seg.setOptimizeCoefficients(true); seg.setModelType(pcl::SACMODEL_PLANE); seg.setMethodType(pcl::SAC_RANSAC); seg.setDistanceThreshold(0.01); // 执行分割 seg.setInputCloud(cloud); seg.segment(*inliers, *coefficients); // 输出分割结果 std::cerr << "Model coefficients: " << coefficients->values[0] << " " << coefficients->values[1] << " " << coefficients->values[2] << " " << coefficients->values[3] << std::endl; std::cerr << "Inliers: " << inliers->indices.size() << std::endl; for (std::size_t i = 0; i < inliers->indices.size(); ++i) std::cerr << inliers->indices[i] << " " << cloud->points[inliers->indices[i]].x << " " << cloud->points[inliers->indices[i]].y << " " << cloud->points[inliers->indices[i]].z << std::endl; return 0; } ``` 这个示例代码首先加载了一个点云文件,并创建了一个`SACSegmentation`对象作为分割器。然后通过设置模型类型以及分割方法类型等参数,最后执行了分割操作。最终输出了分割得到的平面模型系数以及属于该平面点云索引。
评论 6
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值