本小节描述了怎样使用贪婪投影三角化算法对有向点云进行三角化，具体方法是先将有向点云投影到某一局部二维坐标平面内，再在坐标平面内进行平面内的三角化，再根据平面内三位点的拓扑连接关系获得一个三角网格曲面模型。
原理介绍
贪婪投影三角化算法原理是处理一系列可以使网格“生长扩大”的点（边缘点），延伸这些点直到所有符合几何正确性和拓扑正确性的点都被连上。该算法的优点是可以处理来自一个或者多个扫描仪扫描得到并且有多个连接处的散乱点云。但该算法也有一定的局限性，它更适用于采样点云来自于表面连续光滑的曲面并且点云密度变化比较均匀的情况。
参数介绍：
该算法的三角化过程是局部进行的，首先沿着一点的法线将该点投影到局部二维坐标平面内并连接其他悬空点，然后在进行下一点。所以这里我们设置如下参数：
1)函数
SetMaximumNearestNeighbors(unsigned)和SetMu(double)，这两个函数的作用是控制搜索邻域大小。前者定义了可搜索的邻域个数，后者规定了被样本点搜索其邻近点的最远距离，（是为了适应点云密度的变化），特征值一般是50-100和2.5-3（或者1.5每栅格）。
2）函数SetSearchRadius(double)，该函数设置了三角化后得到的每个三角形的最大可能边长。
3）函数SetMinimumAngle(double)和SetMaximumAngle(double)，这两个函数是三角化后每个三角形的最大角和最小角。两者至少要符合一个。典型值分别是10和120度（弧度）。
4）函数SetMaximumSurfaceAgle(double)和SetNormalConsistency(bool)，这两个函数是为了处理边缘或者角很尖锐以及一个表面的两边非常靠近的情况。为了处理这些特殊情况，函数SetMaximumSurfaceAgle(double)规定如果某点法线方向的偏离超过指定角度（注：大多数表面法线估计方法可以估计出连续变化的表面法线方向，即使在尖锐的边缘条件下），该点就不连接到样本点上。该角度是通过计算法向线段（忽略法线方向）之间的角度。函数SetNormalConsistency(bool)保证法线朝向，如果法线方向一致性标识没有设定，就不能保证估计出的法线都可以始终朝向一致。第一个函数特征值为45度（弧度）、第二个函数缺省值为false。
示例代码
将一个xyz点类型的pcd文件打开并存储到对象PointCloud中,由于例子中用到的pcd文件只有XYZ坐标，所以我们把它加载到对象PointCloud< PointXYZ>中，上面代码打开pcd文件，并将点云存储到cloud指针对象中。
pcl::PointCloud<pcl::PointXYZ>::Ptr cloud(new pcl::PointCloud<pcl::PointXYZ>); //定义点云对象指针
#pcl#::PointCloud2 cloud_blob;
pcl::io::loadPCDFile("table_scene_lms400_downsampled.pcd", cloud_blob); //加载bun0.pcd文件
pcl::#fromPCLPointCloud2 #(cloud_blob, *cloud);//数据最终存储在cloud中`
由于本例中使用的三角化算法输入必须为有向点云，所以需要使用PCL中的法线估计方法预先估计出数据中每个点的法线。
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>);//定义kd树指针tree->setInputCloud (cloud); //用cloud构建tree对象n.setInputCloud (cloud); //为法线估计对象设置输入点云
n.setSearchMethod (tree); //设置搜索方法
n.setKSearch (20); //设置k搜索的k值为20
n.compute (*normals); //估计法线存储结果到normals中
由于XYZ坐标字段和法线字段需要在相同PointCloud对象中，所以创建一个新的PointNormal类型的点云来存储坐标字段和法线连接后的点云。
pcl::PointCloud<pcl::PointNormal>::Ptr cloud_with_normals (new pcl::PointCloud<pcl::PointNormal>);
pcl::concatenateFields (*cloud, *normals, *cloud_with_normals);//连接字段，cloud_with_normals存储有向点云
对三角化对象相关变量进行定义。
pcl::search::KdTree<pcl::PointNormal>::Ptr tree2 (new pcl::search::KdTree<pcl::PointNormal>);//定义搜索树对象
tree2->setInputCloud (cloud_with_normals); //利用点云构建搜索树
pcl::GreedyProjectionTriangulation<pcl::PointNormal> gp3;//定义三角化对象
pcl::PolygonMesh triangles; //存储最终三角化的网格模型
设置参数特征值和实际三角化的过程。
gp3.setSearchRadius (0.025);//设置连接点之间的最大距离（即为三角形最大边长）为0.025gp3.setMu (2.5);//设置被样本点搜索其邻近点的最远距离为2.5，为了适应点云密度的变化gp3.setMaximumNearestNeighbors (100); //设置样本点可搜索的邻域个数为100gp3.setMaximumSurfaceAngle(M_PI/4); //设置某点法线方向偏离样本点法线方向的最大角度为45度
gp3.setMinimumAngle(M_PI/18); //设置三角化后得到三角形内角最小角度为10度gp3.setMaximumAngle(2*M_PI/3); //设置三角化后得到三角形内角最大角度为120度gp3.setNormalConsistency(false); //设置该参数保证法线朝向一致gp3.setInputCloud (cloud_with_normals);//设置输入点云为有向点云cloud_with_normals gp3.setSearchMethod (tree2); //设置搜索方式为tree2gp3.reconstruct (triangles); //重建提取三角化
对每个点来说，ID字段中中含有连接组件和该点本身的“状态”(例如 gp3.FREE, gp3.BOUNDARY或 gp3.COMPLETED)可以被检索到。
std::vector<int> parts = gp3.getPartIDs();
std::vector<int> states = gp3.getPointStates();