00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037 #ifndef MODEL_REGISTRATION_H_
00038 #define MODEL_REGISTRATION_H_
00039
00040 #include <pcl/registration/registration.h>
00041 #include <pcl/registration/icp.h>
00042 #include <pcl/registration/icp_nl.h>
00043
00044 #include "tabletop_object_detector/model_fitter.h"
00045
00046 namespace tabletop_object_detector
00047 {
00049 template <typename T>
00050 class ICPRegistrationFitter : public ModelToCloudFitter
00051 {
00052 public:
00053
00054 void initializeFromMesh(const geometric_shapes_msgs::Shape &mesh)
00055 {
00056 std::vector<btVector3> btVectors;
00057 sampleMesh(mesh, btVectors, 0.002 );
00058 cloudmodel_.points.resize (btVectors.size ());
00059 for (int i = 0; i < (int)btVectors.size (); ++i)
00060 {
00061 cloudmodel_.points[i].x = btVectors[i].x ();
00062 cloudmodel_.points[i].y = btVectors[i].y ();
00063 cloudmodel_.points[i].z = btVectors[i].z ();
00064 }
00065 }
00066
00067 ModelFitInfo fitPointCloud (const T &cloud)
00068 {
00069 pcl::IterativeClosestPoint<T, T> reg;
00070 reg.setMaximumIterations (10000);
00071 reg.setTransformationEpsilon (1e-5);
00072
00073
00074
00075
00076 T cloud_out;
00077 reg.compute (cloud_out);
00078
00079 double score = reg.getFitnessScore ();
00080
00081 Eigen::Matrix4f tma = reg.getFinalTransformationMatrix ();
00082 geometry_msgs::Pose p;
00083 p.position.x = tma (0, 3);
00084 p.position.y = tma (1, 3);
00085 p.position.z = tma (2, 3);
00086 Eigen::Quaternion<double> q (tma.block<3,3>(0,0));
00087 p.orientation.x = q.x ();
00088 p.orientation.y = q.y ();
00089 p.orientation.z = q.z ();
00090 p.orientation.w = q.w ();
00091
00092 return ModelFitInfo(model_id_, p, score);
00093 }
00094
00095 private:
00096 T cloudmodel_;
00097 };
00098 }
00099
00100 #endif