|
本帖最后由 just_rookie 于 2013-10-29 10:32 编辑
我写了一段代码,关于marchingcubesRBF,却提示下标越界,请问怎么解决?谢谢!我在网上搜索了一下,有人说是代码有问题,但是我不知道代码错在哪?求助!
代码如下:
#include <iostream>
#include <pcl/common/common.h>
#include <pcl/io/pcd_io.h>
#include <pcl/point_types.h>
#include <pcl/point_cloud.h>
#include <pcl/surface/marching_cubes.h>
#include <pcl/surface/marching_cubes_rbf.h>//RBF;
#include <pcl/features/normal_3d.h>
#include <pcl/search/search.h>//searching method;
#include <pcl/search/kdtree.h>
using namespace std;
int main()
{
pcl::PointCloud<pcl::PointXYZ>::Ptr cloud(new pcl::PointCloud<pcl::PointXYZ>);
if (pcl::io::loadPCDFile<pcl::PointXYZ> ("cloud_6.pcd", *cloud) == -1)
{
PCL_ERROR ("Couldn't read file test_pcd.pcd \n");
return (-1);
}
pcl::search::KdTree<pcl::PointXYZ>::Ptr tree(new pcl::search::KdTree<pcl::PointXYZ>);
tree->setInputCloud(cloud);
pcl::PointCloud<pcl::Normal>::Ptr normal(new pcl::PointCloud<pcl::Normal>); pcl::NormalEstimation<pcl::PointXYZ,pcl::Normal> ne;
ne.setInputCloud(cloud);
ne.setSearchMethod(tree);
ne.setKSearch(20);
ne.compute(*normal);
pcl::PointCloud<pcl::PointNormal>::Ptr cloudWithNormal(new pcl::PointCloud<pcl::PointNormal>); pcl::concatenateFields(*cloud,*normal,*cloudWithNormal);
pcl::search::KdTree<pcl::PointNormal>::Ptr tree1(new pcl::search::KdTree<pcl::PointNormal>);
tree1->setInputCloud(cloudWithNormal);
pcl::MarchingCubesRBF<pcl::PointNormal> rbf;
pcl::PolygonMesh mesh;
rbf.setInputCloud(cloudWithNormal);
rbf.setSearchMethod(tree1);
rbf.reconstruct(mesh);
cout<<"OK!"<<endl;
return 0;
}
|
|