Merge pull request #1633 from sgiraudot/Point_set_processing_3-Point_set_structuring-GF

Point Set Structuring

Point_set_processing_3/examples/Point_set_processing_3/CMakeLists.txt

@@ -55,6 +55,7 @@ if ( CGAL_FOUND )
   create_single_source_cgal_program( "remove_outliers_example.cpp" )
   create_single_source_cgal_program( "wlop_simplify_and_regularize_point_set_example.cpp" )
   create_single_source_cgal_program( "edge_aware_upsample_point_set_example.cpp" )
+  create_single_source_cgal_program( "structuring_example.cpp" )
   
   # Use Eigen or BLAS and LAPACK (optional)
   find_package(Eigen3 3.1.0) #(requires 3.1.0 or greater)

Point_set_processing_3/examples/Point_set_processing_3/structuring_example.cpp

@@ -0,0 +1,70 @@
+#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
+#include <CGAL/IO/read_xyz_points.h>
+#include <CGAL/IO/write_xyz_points.h>
+#include <CGAL/Point_with_normal_3.h>
+#include <CGAL/property_map.h>
+
+#include <CGAL/Shape_detection_3.h>
+
+#include <CGAL/structure_point_set.h>
+
+#include <iostream>
+#include <fstream>
+
+// Type declarations
+typedef CGAL::Exact_predicates_inexact_constructions_kernel  Kernel;
+typedef Kernel::Point_3                                      Point;
+typedef std::pair<Kernel::Point_3, Kernel::Vector_3>         Point_with_normal;
+typedef std::vector<Point_with_normal>                       Pwn_vector;
+typedef CGAL::First_of_pair_property_map<Point_with_normal>  Point_map;
+typedef CGAL::Second_of_pair_property_map<Point_with_normal> Normal_map;
+
+// Efficient RANSAC types
+typedef CGAL::Shape_detection_3::Efficient_RANSAC_traits
+  <Kernel, Pwn_vector, Point_map, Normal_map>                Traits;
+typedef CGAL::Shape_detection_3::Efficient_RANSAC<Traits>    Efficient_ransac;
+typedef CGAL::Shape_detection_3::Plane<Traits>               Plane;
+
+int main (int argc, char** argv)
+{
+  // Points with normals.
+  Pwn_vector points;
+
+  // Loading point set from a file. 
+  std::ifstream stream(argc>1 ? argv[1] : "data/cube.pwn");
+
+  if (!stream || 
+    !CGAL::read_xyz_points_and_normals(stream,
+      std::back_inserter(points),
+      Point_map(),
+      Normal_map()))
+  {
+      std::cerr << "Error: cannot read file cube.pwn" << std::endl;
+      return EXIT_FAILURE;
+  }
+
+  std::cerr << points.size() << " point(s) read." << std::endl;
+  
+  // Shape detection
+  Efficient_ransac ransac;
+  ransac.set_input(points);
+  ransac.add_shape_factory<Plane>();
+  ransac.detect();
+
+  Pwn_vector structured_pts;
+
+  CGAL::structure_point_set (points.begin (), points.end (), // input points
+                             std::back_inserter (structured_pts),
+                             ransac, // shape detection engine
+                             0.015); // epsilon for structuring points
+
+  std::cerr << structured_pts.size ()
+            << " structured point(s) generated." << std::endl;
+
+  std::ofstream out ("out.pwn");
+  CGAL::write_xyz_points_and_normals (out, structured_pts.begin(), structured_pts.end(),
+                                      Point_map(), Normal_map());
+  out.close();
+  
+  return EXIT_SUCCESS;
+}