*点群からマーチングキューブス(Hoppe)で面を作成(PCLライブラリ)  [#k0eb16de]

[[PCL:http://pointclouds.org/]]ライブラリを使用して、マーチングキューブス法(Hoppe)で点群から面を作成します。

ここでは、ねこの点群を入力にします。

#ref(neko_cloud_ponits.png,center,nowrap,60%,ネコ点群);

では早速、下のプログラムでマーチングキューブス法で面を作成してみます。

#ref(neko_surface1.png,center,nowrap,60%,ネコ面);

うまくいきません。。ねこの周りに、いらない面が残ってしまいます。
うまくいきません。。ねこの周りに、いらない面が残ってしまいます。0の等値面はどのように定義されているのでしょうか?

**pcl_marching_cubes_hoppe.cxx [#y88ae58d]
#code(c){{
#include <pcl/point_types.h>
#include <pcl/io/pcd_io.h>
#include <pcl/io/vtk_io.h>
#include <pcl/features/normal_3d.h>
#include <pcl/surface/marching_cubes_hoppe.h>

int
main (int argc, char** argv)
{
  std::string input_file, output_file;

  if (argc < 3)
  {
    printf ("\nUsage: marching_cubes_hoppe pcd<PointXYZ>-in-file vtk-out-file\n\n");
    exit(0);
  }
  input_file = argv[1];
  output_file = argv[2];

  // Load input file into a PointCloud<T> with an appropriate type
  pcl::PointCloud<pcl::PointXYZ>::Ptr cloud (new pcl::PointCloud<pcl::PointXYZ>);
  pcl::PCLPointCloud2 cloud_blob;
  pcl::io::loadPCDFile (input_file, cloud_blob);
  pcl::fromPCLPointCloud2 (cloud_blob, *cloud);
  //* the data should be available in cloud

  // Normal estimation*
  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);
  //* normals should not contains the point normals + surface curvatures

  // Concatenate the XYZ and normal fields*
  pcl::PointCloud<pcl::PointNormal>::Ptr cloud_with_normals (new pcl::PointCloud<pcl::PointNormal>);
  pcl::concatenateFields (*cloud, *normals, *cloud_with_normals);
  //* cloud_with_normals = cloud + normals

  // Create search tree
  pcl::search::KdTree<pcl::PointNormal>::Ptr tree2 (new pcl::search::KdTree<pcl::PointNormal>);
  tree2->setInputCloud (cloud_with_normals);

  // Initialize objects
  pcl::MarchingCubesHoppe<pcl::PointNormal> mc;
  pcl::PolygonMesh triangles;

  // Set parameters
  mc.setIsoLevel (0.0);
  mc.setGridResolution(100, 100, 100);
  mc.setPercentageExtendGrid(0.3f);

  // Get result
  mc.setInputCloud (cloud_with_normals);
  mc.setSearchMethod (tree2);
  mc.reconstruct (triangles);

  pcl::io::saveVTKFile (output_file, triangles);

  // Finish
  return (0);
}

}}

**ダウンロードとビルド [#k551ac14]
ソースコードと'''CMakeLists.txt'''ファイル:
#ref(pcl_marching_cubes_hoppe.zip)

トップ   編集 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS