使用pcl::ApproximateVoxelGrid时出现异常现象求助

环境与操作流程 #

使用PCL 1.11.1 + VS2019，通过CloudCompare生成规整网格状的合成点云，采用ApproximateVoxelGrid进行下采样，代码如下：

pcl::PointCloud<pcl::PointXYZ>::Ptr cloud_in(new pcl::PointCloud<pcl::PointXYZ>);
pcl::io::loadPLYFile("data_path", *cloud_in);

pcl::ApproximateVoxelGrid<pcl::PointXYZ> AVG_filter;
AVG_filter.setInputCloud(cloud_in);
AVG_filter.setLeafSize(2.0f, 2.0f, 2.0f);
AVG_filter.filter(*cloud_in);

异常现象 #

原始点云为规则排列的网格结构，滤波后可视化输出出现大量离散、无规则的异常点，点云原有结构完全破坏。

问题原因与解决方案 #

核心原因：内存读写冲突 #

代码中直接将滤波结果写入输入点云对象cloud_in，而ApproximateVoxelGrid在滤波过程中会同时读取输入云数据并修改其内存空间，这种读写同一块内存的操作会导致数据错乱，进而出现异常点。

修复方案 #

新建独立的输出点云对象，避免覆盖输入云：

pcl::PointCloud<pcl::PointXYZ>::Ptr cloud_in(new pcl::PointCloud<pcl::PointXYZ>);
pcl::io::loadPLYFile("data_path", *cloud_in);

// 创建独立的输出点云指针
pcl::PointCloud<pcl::PointXYZ>::Ptr cloud_filtered(new pcl::PointCloud<pcl::PointXYZ>);
pcl::ApproximateVoxelGrid<pcl::PointXYZ> AVG_filter;
AVG_filter.setInputCloud(cloud_in);
AVG_filter.setLeafSize(2.0f, 2.0f, 2.0f);
// 将滤波结果写入新对象
AVG_filter.filter(*cloud_filtered);

// 可视化时使用cloud_filtered

额外排查项 #

1. 过滤无效点：检查CloudCompare导出的PLY文件是否包含NaN值点，这类点会干扰滤波计算，可提前过滤：

pcl::PointCloud<pcl::PointXYZ>::Ptr cloud_clean(new pcl::PointCloud<pcl::PointXYZ>);
std::vector<int> valid_indices;
pcl::removeNaNFromPointCloud(*cloud_in, *cloud_clean, valid_indices);
// 使用cloud_clean作为滤波输入

2. 体素尺寸适配：若2.0f的体素尺寸远大于原始点云的点间距，可适当调小LeafSize值，避免体素内采样逻辑异常。

内容的提问来源于stack exchange，提问作者felix Zhao