transform_from_viewpoint.cpp
Go to the documentation of this file.
00001 /*
00002  * Software License Agreement (BSD License)
00003  *
00004  *  Point Cloud Library (PCL) - www.pointclouds.org
00005  *  Copyright (c) 2010-2012, Willow Garage, Inc.
00006  *  
00007  *  All rights reserved.
00008  *
00009  *  Redistribution and use in source and binary forms, with or without
00010  *  modification, are permitted provided that the following conditions
00011  *  are met:
00012  *
00013  *   * Redistributions of source code must retain the above copyright
00014  *     notice, this list of conditions and the following disclaimer.
00015  *   * Redistributions in binary form must reproduce the above
00016  *     copyright notice, this list of conditions and the following
00017  *     disclaimer in the documentation and/or other materials provided
00018  *     with the distribution.
00019  *   * Neither the name of the copyright holder(s) nor the names of its
00020  *     contributors may be used to endorse or promote products derived
00021  *     from this software without specific prior written permission.
00022  *
00023  *  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
00024  *  "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
00025  *  LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
00026  *  FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
00027  *  COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
00028  *  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
00029  *  BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
00030  *  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
00031  *  CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
00032  *  LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
00033  *  ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
00034  *  POSSIBILITY OF SUCH DAMAGE.
00035  *
00036  * $Id$
00037  *
00038  */
00039 
00040 #include <pcl/PCLPointCloud2.h>
00041 #include <pcl/io/pcd_io.h>
00042 #include <pcl/features/fpfh.h>
00043 #include <pcl/console/print.h>
00044 #include <pcl/console/parse.h>
00045 #include <pcl/console/time.h>
00046 #include <pcl/common/transforms.h>
00047 
00048 using namespace pcl;
00049 using namespace pcl::io;
00050 using namespace pcl::console;
00051 
00052 Eigen::Vector4f    translation;
00053 Eigen::Quaternionf orientation;
00054 
00055 void
00056 printHelp (int, char **argv)
00057 {
00058   print_error ("Syntax is: %s input.pcd output.pcd\n", argv[0]);
00059 }
00060 
00061 bool
00062 loadCloud (const std::string &filename, pcl::PCLPointCloud2 &cloud)
00063 {
00064   TicToc tt;
00065   print_highlight ("Loading "); print_value ("%s ", filename.c_str ());
00066 
00067   tt.tic ();
00068   if (loadPCDFile (filename, cloud, translation, orientation) < 0)
00069     return (false);
00070   print_info ("[done, "); print_value ("%g", tt.toc ()); print_info (" ms : "); print_value ("%d", cloud.width * cloud.height); print_info (" points]\n");
00071   print_info ("Available dimensions: "); print_value ("%s\n", getFieldsList (cloud).c_str ());
00072 
00073   return (true);
00074 }
00075 
00076 void
00077 transform (const pcl::PCLPointCloud2::ConstPtr &input, pcl::PCLPointCloud2 &output)
00078 {
00079   // Check for 'normals'
00080   bool has_normals = false;
00081   for (size_t i = 0; i < input->fields.size (); ++i)
00082     if (input->fields[i].name == "normals")
00083       has_normals = true;
00084 
00085   // Estimate
00086   TicToc tt;
00087   tt.tic ();
00088   print_highlight (stderr, "Transforming ");
00089 
00090   // Convert data to PointCloud<T>
00091   if (has_normals)
00092   {
00093     PointCloud<PointNormal> xyznormals;
00094     fromPCLPointCloud2 (*input, xyznormals);
00095     pcl::transformPointCloud<PointNormal> (xyznormals, xyznormals, translation.head<3> (), orientation);
00096     // Copy back the xyz and normals
00097     pcl::PCLPointCloud2 output_xyznormals;
00098     toPCLPointCloud2 (xyznormals, output_xyznormals);
00099     concatenateFields (*input, output_xyznormals, output);
00100   }
00101   else
00102   {
00103     PointCloud<PointXYZ> xyz;
00104     fromPCLPointCloud2 (*input, xyz);
00105     pcl::transformPointCloud<PointXYZ> (xyz, xyz, translation.head<3> (), orientation);
00106     // Copy back the xyz and normals
00107     pcl::PCLPointCloud2 output_xyz;
00108     toPCLPointCloud2 (xyz, output_xyz);
00109     concatenateFields (*input, output_xyz, output);
00110   }
00111 
00112   translation = Eigen::Vector4f::Zero ();
00113   orientation = Eigen::Quaternionf::Identity ();
00114 
00115   print_info ("[done, "); print_value ("%g", tt.toc ()); print_info (" ms : "); print_value ("%d", output.width * output.height); print_info (" points]\n");
00116 }
00117 
00118 void
00119 saveCloud (const std::string &filename, const pcl::PCLPointCloud2 &output)
00120 {
00121   TicToc tt;
00122   tt.tic ();
00123 
00124   print_highlight ("Saving "); print_value ("%s ", filename.c_str ());
00125   
00126   PCDWriter writer;
00127   writer.writeBinaryCompressed (filename, output, translation, orientation);
00128   
00129   print_info ("[done, "); print_value ("%g", tt.toc ()); print_info (" ms : "); print_value ("%d", output.width * output.height); print_info (" points]\n");
00130 }
00131 
00132 /* ---[ */
00133 int
00134 main (int argc, char** argv)
00135 {
00136   print_info ("Take the input point cloud and transform it according to its stored VIEWPOINT information. For more information, use %s -h\n", argv[0]);
00137   bool help = false;
00138   parse_argument (argc, argv, "-h", help);
00139   if (argc < 3 || help)
00140   {
00141     printHelp (argc, argv);
00142     return (-1);
00143   }
00144 
00145   // Parse the command line arguments for .pcd files
00146   std::vector<int> p_file_indices;
00147   p_file_indices = parse_file_extension_argument (argc, argv, ".pcd");
00148   if (p_file_indices.size () != 2)
00149   {
00150     print_error ("Need one input PCD file and one output PCD file to continue.\n");
00151     return (-1);
00152   }
00153 
00154   // Load the first file
00155   pcl::PCLPointCloud2::Ptr cloud (new pcl::PCLPointCloud2);
00156   if (!loadCloud (argv[p_file_indices[0]], *cloud)) 
00157     return (-1);
00158 
00159   // Perform the feature estimation
00160   pcl::PCLPointCloud2 output;
00161   transform (cloud, output);
00162 
00163   // Save into the second file
00164   saveCloud (argv[p_file_indices[1]], output);
00165 }
00166 


pcl
Author(s): Open Perception
autogenerated on Wed Aug 26 2015 15:36:53