mesh_segmentation.h

Go to the documentation of this file.

#ifndef CGAL_SURFACE_MESH_SEGMENTATION_MESH_SEGMENTATION_H
#define CGAL_SURFACE_MESH_SEGMENTATION_MESH_SEGMENTATION_H

#include <CGAL/internal/Surface_mesh_segmentation/Surface_mesh_segmentation.h>
#include <boost/config.hpp>

namespace CGAL {

  
template <class Polyhedron, class SDFPropertyMap, class GeomTraits
#ifndef BOOST_NO_FUNCTION_TEMPLATE_DEFAULT_ARGS
  = typename Polyhedron::Traits
#endif
>
std::pair<double, double> 
compute_sdf_values(const Polyhedron& polyhedron, 
                        SDFPropertyMap sdf_values, 
                        double cone_angle = 2.0 / 3.0 * CGAL_PI,
                        int number_of_rays = 25,
                        GeomTraits traits = GeomTraits())
{
    internal::Surface_mesh_segmentation<Polyhedron, GeomTraits> algorithm(polyhedron, traits);
    return algorithm.calculate_sdf_values(cone_angle, number_of_rays, sdf_values);
}


template <class Polyhedron, class SDFPropertyMap, class SegmentPropertyMap, class GeomTraits
#ifndef BOOST_NO_FUNCTION_TEMPLATE_DEFAULT_ARGS
  = typename Polyhedron::Traits
#endif
>
int 
segment_from_sdf_values(const Polyhedron& polyhedron, 
                                           SDFPropertyMap sdf_values, 
                                           SegmentPropertyMap segment_ids,
                                           int number_of_levels = 5, 
                                           double smoothing_lambda = 0.26,
                                           GeomTraits traits = GeomTraits())
{
    smoothing_lambda = (std::max)(0.0, (std::min)(1.0, smoothing_lambda)); // clip into [0-1]
    
    internal::Surface_mesh_segmentation<Polyhedron, GeomTraits> algorithm(polyhedron, traits);
    return algorithm.partition(number_of_levels, smoothing_lambda, sdf_values, segment_ids);
}


template < class Polyhedron, class SegmentPropertyMap, class GeomTraits
#ifndef BOOST_NO_FUNCTION_TEMPLATE_DEFAULT_ARGS
  = typename Polyhedron::Traits
#endif
>
int 
compute_sdf_values_and_segment(const Polyhedron& polyhedron,                                
                           SegmentPropertyMap segment_ids,
                           double cone_angle = 2.0 / 3.0 * CGAL_PI,
                           int number_of_rays = 25,
                           int number_of_levels = 5, 
                           double smoothing_lambda = 0.26,
                           GeomTraits traits = GeomTraits())
{
    smoothing_lambda = (std::max)(0.0, (std::min)(1.0, smoothing_lambda)); // clip into [0-1]
    
    typedef std::map< typename Polyhedron::Facet_const_handle, double> Facet_double_map;
    Facet_double_map internal_sdf_map;
    boost::associative_property_map<Facet_double_map> sdf_property_map(internal_sdf_map);
    
    compute_sdf_values<Polyhedron, boost::associative_property_map<Facet_double_map>, GeomTraits>
        (polyhedron, sdf_property_map, cone_angle, number_of_rays, traits);
    return segment_from_sdf_values<Polyhedron, boost::associative_property_map<Facet_double_map>, SegmentPropertyMap, GeomTraits>
        (polyhedron, sdf_property_map, segment_ids, number_of_levels, smoothing_lambda, traits);
}


#ifdef BOOST_NO_FUNCTION_TEMPLATE_DEFAULT_ARGS
template <class Polyhedron, class SDFPropertyMap>
std::pair<double, double> 
compute_sdf_values(const Polyhedron& polyhedron, 
                        SDFPropertyMap sdf_values, 
                        double cone_angle = 2.0 / 3.0 * CGAL_PI,
                        int number_of_rays = 25,
                        typename Polyhedron::Traits traits = typename Polyhedron::Traits())
{
    return compute_sdf_values<Polyhedron, SDFPropertyMap, typename Polyhedron::Traits>
        (polyhedron, sdf_values, cone_angle, number_of_rays, traits);
}

template <class Polyhedron, class SDFPropertyMap, class SegmentPropertyMap>
int 
segment_from_sdf_values(const Polyhedron& polyhedron, 
                                           SDFPropertyMap sdf_values, 
                                           SegmentPropertyMap segment_ids,
                                           int number_of_levels = 5, 
                                           double smoothing_lambda = 0.26,
                                           typename Polyhedron::Traits traits = typename Polyhedron::Traits())
{
    return segment_from_sdf_values<Polyhedron, SDFPropertyMap, SegmentPropertyMap, typename Polyhedron::Traits>
        (polyhedron, sdf_values, segment_ids, number_of_levels, smoothing_lambda, traits);
}

template <class Polyhedron, class SegmentPropertyMap>
int 
compute_sdf_values_and_segment(const Polyhedron& polyhedron,                                
                           SegmentPropertyMap segment_ids,
                           double cone_angle = 2.0 / 3.0 * CGAL_PI,
                           int number_of_rays = 25,
                           int number_of_levels = 5, 
                           double smoothing_lambda = 0.26,
                           typename Polyhedron::Traits traits = typename Polyhedron::Traits())
{
    return compute_sdf_values_and_segment<Polyhedron, SegmentPropertyMap, typename Polyhedron::Traits>
        (polyhedron, segment_ids, cone_angle, number_of_rays, number_of_levels, smoothing_lambda, traits);
}
#endif

}//namespace CGAL

#endif // CGAL_SURFACE_MESH_SEGMENTATION_MESH_SEGMENTATION_H //