PCL實現泊松表面重建

pcl實現的泊松重建算法，輸出爲重建後的曲面模型數據。具體內容可以參照文獻Possion surface reconstruction.#include <pcl/point_types.h>

#include <pcl/io/pcd_io.h>

#include <pcl/kdtree/kdtree_flann.h>

#include <pcl/features/normal_3d_omp.h>

#include <pcl/features/normal_3d.h>

#include <pcl/surface/gp3.h>

#include <pcl/surface/poisson.h>

#include <pcl/visualization/pcl_visualizer.h>

#include <boost/thread/thread.hpp>

#include <fstream>

void main()

{  pcl::PointCloud<pcl::PointXYZ>::Ptr cloud(new pcl::PointCloud<pcl::PointXYZ>) ;  

if (pcl::io::loadPCDFile("bun000.pcd" , *cloud) == -1)  

{ PCL_ERROR("Could not read pcd file!\n") ;  

return -1;  

}  

pcl::PointCloud<pcl::PointNormal>::Ptr cloud_with_normals(new pcl::PointCloud<pcl::PointNormal>) ;  

pcl::NormalEstimation<pcl::PointXYZ , pcl::Normal> n ;//法線估計對象  

pcl::PointCloud<pcl::Normal>::Ptr normals(new pcl::PointCloud<pcl::Normal>) ;//存儲估計的法線  

pcl::search::KdTree<pcl::PointXYZ>::Ptr tree(new pcl::search::KdTree<pcl::PointXYZ>) ;  

tree->setInputCloud(cloud) ;  

n.setInputCloud(cloud) ; 

n.setSearchMethod(tree) ;  

n.setKSearch(20) ;  

n.compute(*normals) ;  

pcl::concatenateFields(*cloud , *normals , *cloud_with_normals) ;  

pcl::search::KdTree<pcl::PointNormal>::Ptr tree2(new pcl::search::KdTree<pcl::PointNormal>) ;  

tree2->setInputCloud(cloud_with_normals) ;  

pcl::Poisson<pcl::PointNormal> pn ;

pn.setSearchMethod(tree2) ; 

pn.setInputCloud(cloud_with_normals) ; 

pn.setConfidence(false) ;//設置置信標誌，爲true時，使用法線向量長度作爲置信度信息，false則需要對法線進行歸一化處理  

pn.setManifold(false) ;//設置流行標誌，如果設置爲true，則對多邊形進行細分三角話時添加重心，設置false則不添加  

pn.setOutputPolygons(false) ;//設置是否輸出爲多邊形  

pn.setIsoDivide(8) ;  

pn.setSamplesPerNode(3) ;//設置每個八叉樹節點上最少採樣點數目  

pcl::PolygonMesh mesh ; 

pn.performReconstruction(mesh) ; 

boost::shared_ptr<pcl::visualization::PCLVisualizer> viewer(new pcl::visualization::PCLVisualizer("3D viewer")) ; 

viewer->setBackgroundColor(0 , 0 , 0) ; 

viewer->addPolygonMesh(mesh , "my") ; 

viewer->initCameraParameters() ;  

while (!viewer->wasStopped())  

{ viewer->spinOnce(100) ;  

boost::this_thread::sleep(boost::posix_time::microseconds(100000)) ;  

}  

return 0;

}