affordances.h

Go to the documentation of this file.

/*
 * Software License Agreement (BSD License)
 *
 *  Copyright (c) 2014, Andreas ten Pas
 *  All rights reserved.
 *
 *  Redistribution and use in source and binary forms, with or without
 *  modification, are permitted provided that the following conditions
 *  are met:
 *
 *   * Redistributions of source code must retain the above copyright
 *     notice, this list of conditions and the following disclaimer.
 *   * Redistributions in binary form must reproduce the above
 *     copyright notice, this list of conditions and the following
 *     disclaimer in the documentation and/or other materials provided
 *     with the distribution.
 *
 *  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
 *  "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
 *  LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
 *  FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
 *  COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
 *  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
 *  BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
 *  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
 *  CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
 *  LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
 *  ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 *  POSSIBILITY OF SUCH DAMAGE.
 */

#ifndef AFFORDANCES_H
#define AFFORDANCES_H

#include <fstream>
#include <iostream>
#include <omp.h>
#include <pcl/features/feature.h>
#include <pcl/features/normal_3d.h>
#include <pcl/features/normal_3d_omp.h>
#include <pcl/kdtree/kdtree_flann.h>
#include <pcl/point_cloud.h>
#include <pcl/point_types.h>
#include <ros/ros.h>
#include <stdlib.h>
#include <string>
#include <tf/transform_datatypes.h>
#include "curvature_estimation_taubin.h"
#include "curvature_estimation_taubin.hpp"
#include "cylindrical_shell.h"


typedef pcl::PointCloud<pcl::PointXYZRGBA> PointCloud;


// workspace limits of the robot
struct WorkspaceLimits
{
        double min_x;
        double max_x;
        double min_y;
        double max_y;
        double min_z;
        double max_z;
};


class Affordances
{
        public:

                void 
    initParams(ros::NodeHandle node);
    
                PointCloud::Ptr 
    maxRangeFilter(const PointCloud::Ptr &cloud_in);
    
                PointCloud::Ptr 
    workspaceFilter(const PointCloud::Ptr &cloud_in, tf::StampedTransform *transform = NULL);

    std::vector<CylindricalShell> 
    searchAffordances(const PointCloud::Ptr &cloud, tf::StampedTransform *transform = NULL);
    
    std::vector<CylindricalShell> 
    searchAffordances(const PointCloud::Ptr &cloud, const std::vector<int> &indices);
    
    std::vector<CylindricalShell> 
    searchAffordancesTaubin(const PointCloud::Ptr &cloud, const Eigen::MatrixXd &samples, 
      bool is_logging = true);
        
    std::vector< std::vector<CylindricalShell> > 
    searchHandles(const PointCloud::Ptr &cloud, std::vector<CylindricalShell> shells);
    
    std::vector<int> 
    createRandomIndices(const PointCloud::Ptr &cloud, int size);
    
                bool 
    isPointInWorkspace(double x, double y, double z, tf::StampedTransform *transform = NULL);
    
    std::string getPCDFile() { return this->file; }
    
    int getNumSamples() { return this->num_samples; }
    
    double getTargetRadius() { return this->target_radius; }
      
  
        private:    
  
    void 
    estimateNormals(const PointCloud::Ptr &cloud, 
                    const pcl::PointCloud<pcl::Normal>::Ptr &cloud_normals);
    
                void 
    estimateCurvatureAxisPCA(const PointCloud::Ptr &cloud, int nn_center_idx, 
                            std::vector<int> nn_indices, Eigen::Vector3d &axis, 
                            Eigen::Vector3d &normal);

    void 
    estimateCurvatureAxisNormals(const pcl::PointCloud<pcl::Normal>::Ptr &cloud_normals, 
                                const std::vector<int> &nn_indices, 
                                Eigen::Vector3d &axis, Eigen::Vector3d &normal);
  
    std::vector<CylindricalShell> 
    searchAffordancesNormalsOrPCA(const PointCloud::Ptr &cloud, 
      tf::StampedTransform *transform = NULL);
                                
    std::vector<CylindricalShell> 
    searchAffordancesTaubin(const PointCloud::Ptr &cloud, tf::StampedTransform *transform = NULL);
    
    void findBestColinearSet(const std::vector<CylindricalShell> &list, 
                            std::vector<int> &inliersMaxSet, 
                            std::vector<int> &outliersMaxSet);
                            
                int numInFront(const PointCloud::Ptr &cloud, int center_index, double radius);

    // parameters (read-in from ROS launch file)
                double target_radius;
                double radius_error;
                double handle_gap;
                int num_samples;
                double max_range;
                bool use_clearance_filter;
                bool use_occlusion_filter;
                int curvature_estimator;
                int alignment_runs;
                int alignment_min_inliers;
                double alignment_dist_radius;
                double alignment_orient_radius;
                double alignment_radius_radius;
                WorkspaceLimits workspace_limits;
                int num_threads;
    std::string file;
                
                // standard parameters
                static const int CURVATURE_ESTIMATOR; // curvature axis estimation method
                static const int NUM_SAMPLES; // number of neighborhoods
                static const int NUM_NEAREST_NEIGHBORS; // number of nearest neighbors to be found
                static const double NEIGHBOR_RADIUS;
                static const int MAX_NUM_IN_FRONT; // max. threshold of allowed points in front of a neighborhood center point (occlusion filtering)
                static const double TARGET_RADIUS; // approx. radius of the target handle
                static const double RADIUS_ERROR; // allowed deviation from target radius
                static const double HANDLE_GAP; // min. gap around affordances
                static const double MAX_RANGE; // max. range of robot arms
                static const bool USE_CLEARANCE_FILTER; // whether the clearance filter is used
                static const bool USE_OCCLUSION_FILTER; // whether the occlusion filter is used
                static const int ALIGNMENT_RUNS; // number of RANSAC runs
                static const int ALIGNMENT_MIN_INLIERS; // min. number of inliers for colinear cylinder set
                static const double ALIGNMENT_DIST_RADIUS; // distance threshold
                static const double ALIGNMENT_ORIENT_RADIUS; // orientation threshold
                static const double ALIGNMENT_RADIUS_RADIUS; // radius threshold
                static const double WORKSPACE_MIN;
                static const double WORKSPACE_MAX;
};

#endif