PCL泊松重建结果异常，请求排查原因

我尝试加载类似bunny的模型文件并执行泊松重建，使用以下代码实现，但得到了异常的网格结果。请问该异常是什么原因导致的？

#include <pcl/point_types.h>
#include <pcl/point_cloud.h>
#include <pcl/point_representation.h>
#include <pcl/io/pcd_io.h>
#include <iostream>
#include <pcl/io/obj_io.h>
#include <pcl/io/ply_io.h>
#include <pcl/surface/poisson.h>
#include <pcl/features/normal_3d_omp.h>
#include <pcl/features/normal_3d.h>
#include <pcl/visualization/pcl_visualizer.h>

using namespace pcl;
using namespace std;

int testRecon()
{
    
    string plyfile = R"(E:\bunny2.obj)";
    PointCloud<PointNormal>::Ptr cloud_in(new PointCloud<PointNormal>());
    int ret = pcl::io::loadOBJFile(plyfile, *cloud_in);
    std::cout << "Ret: " << ret;

    Poisson<PointNormal> poisson;
    poisson.setDepth(8);
    poisson.setInputCloud(cloud_in);
    PolygonMesh mesh;
    poisson.reconstruct(mesh);

    string file_all = R"(E:\mesh.ply)";
    pcl::io::savePLYFile(file_all, mesh, true);

    std::cout << "done\n";
}

异常原因 #

• 法向量缺失/无效：泊松重建算法严重依赖点云的法向量信息判断表面朝向，但OBJ文件本身只存储顶点坐标，不包含法向量数据。你用loadOBJFile直接加载到PointCloud<PointNormal>中时，法向量字段是未初始化的垃圾值，完全破坏了重建逻辑，导致生成畸形网格。
• 文件加载逻辑不匹配：代码中变量名标注为plyfile却指向OBJ文件，且错误地用OBJ加载函数直接适配带法向量的点云类型，进一步加剧了数据无效性。

修正方案 #

1. 先加载纯坐标点云，再计算法向量：OBJ文件无自带法向量，需先加载为PointCloud<PointXYZ>，再用法向量估计算法生成法向量，最后合并为PointCloud<PointNormal>作为泊松重建的输入。
2. 匹配文件格式与加载函数：加载OBJ就用loadOBJFile对应PointXYZ，加载PLY则用loadPLYFile（若PLY带法向量可直接加载到PointNormal）。

修正后的示例代码：

#include <pcl/point_types.h>
#include <pcl/point_cloud.h>
#include <pcl/io/obj_io.h>
#include <pcl/io/ply_io.h>
#include <pcl/surface/poisson.h>
#include <pcl/features/normal_3d_omp.h>
#include <iostream>

using namespace pcl;
using namespace std;

int testRecon()
{
    // 加载OBJ文件到纯坐标点云
    string objfile = R"(E:\bunny2.obj)";
    PointCloud<PointXYZ>::Ptr cloud_xyz(new PointCloud<PointXYZ>());
    int ret = pcl::io::loadOBJFile(objfile, *cloud_xyz);
    if (ret == -1) {
        cerr << "加载OBJ文件失败！" << endl;
        return -1;
    }
    cout << "加载了 " << cloud_xyz->size() << " 个点" << endl;

    // 并行计算法向量（根据点云密度调整搜索半径）
    NormalEstimationOMP<PointXYZ, Normal> ne;
    ne.setInputCloud(cloud_xyz);
    ne.setRadiusSearch(0.01);
    PointCloud<Normal>::Ptr normals(new PointCloud<Normal>());
    ne.compute(*normals);

    // 合并坐标与法向量
    PointCloud<PointNormal>::Ptr cloud_in(new PointCloud<PointNormal>());
    concatenateFields(*cloud_xyz, *normals, *cloud_in);

    // 泊松重建
    Poisson<PointNormal> poisson;
    poisson.setDepth(8);
    poisson.setInputCloud(cloud_in);
    PolygonMesh mesh;
    poisson.reconstruct(mesh);

    // 保存重建结果
    string file_all = R"(E:\mesh.ply)";
    pcl::io::savePLYFile(file_all, mesh, true);

    cout << "完成\n";
    return 0;
}

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