Chris@16: // Boost.Geometry (aka GGL, Generic Geometry Library)
Chris@16: 
Chris@16: // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
Chris@16: // Copyright (c) 2008-2012 Bruno Lalande, Paris, France.
Chris@16: // Copyright (c) 2009-2012 Mateusz Loskot, London, UK.
Chris@16: 
Chris@16: // Parts of Boost.Geometry are redesigned from Geodan's Geographic Library
Chris@16: // (geolib/GGL), copyright (c) 1995-2010 Geodan, Amsterdam, the Netherlands.
Chris@16: 
Chris@16: // Use, modification and distribution is subject to the Boost Software License,
Chris@16: // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
Chris@16: // http://www.boost.org/LICENSE_1_0.txt)
Chris@16: 
Chris@16: #ifndef BOOST_GEOMETRY_ALGORITHMS_ENVELOPE_HPP
Chris@16: #define BOOST_GEOMETRY_ALGORITHMS_ENVELOPE_HPP
Chris@16: 
Chris@101: #include <vector>
Chris@101: 
Chris@101: #include <boost/numeric/conversion/cast.hpp>
Chris@16: #include <boost/range.hpp>
Chris@16: 
Chris@101: #include <boost/variant/apply_visitor.hpp>
Chris@101: #include <boost/variant/static_visitor.hpp>
Chris@101: #include <boost/variant/variant_fwd.hpp>
Chris@16: 
Chris@16: #include <boost/geometry/algorithms/assign.hpp>
Chris@16: #include <boost/geometry/algorithms/expand.hpp>
Chris@16: #include <boost/geometry/algorithms/not_implemented.hpp>
Chris@16: #include <boost/geometry/core/cs.hpp>
Chris@16: #include <boost/geometry/core/exterior_ring.hpp>
Chris@101: #include <boost/geometry/core/point_type.hpp>
Chris@101: #include <boost/geometry/core/tags.hpp>
Chris@101: 
Chris@16: #include <boost/geometry/geometries/concepts/check.hpp>
Chris@16: 
Chris@16: 
Chris@16: namespace boost { namespace geometry
Chris@16: {
Chris@16: 
Chris@16: #ifndef DOXYGEN_NO_DETAIL
Chris@16: namespace detail { namespace envelope
Chris@16: {
Chris@16: 
Chris@16: 
Chris@16: /// Calculate envelope of an 2D or 3D segment
Chris@16: struct envelope_expand_one
Chris@16: {
Chris@16:     template<typename Geometry, typename Box>
Chris@16:     static inline void apply(Geometry const& geometry, Box& mbr)
Chris@16:     {
Chris@16:         assign_inverse(mbr);
Chris@16:         geometry::expand(mbr, geometry);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: /// Iterate through range (also used in multi*)
Chris@16: template<typename Range, typename Box>
Chris@16: inline void envelope_range_additional(Range const& range, Box& mbr)
Chris@16: {
Chris@16:     typedef typename boost::range_iterator<Range const>::type iterator_type;
Chris@16: 
Chris@16:     for (iterator_type it = boost::begin(range);
Chris@16:         it != boost::end(range);
Chris@16:         ++it)
Chris@16:     {
Chris@16:         geometry::expand(mbr, *it);
Chris@16:     }
Chris@16: }
Chris@16: 
Chris@16: 
Chris@16: 
Chris@16: /// Generic range dispatching struct
Chris@16: struct envelope_range
Chris@16: {
Chris@16:     /// Calculate envelope of range using a strategy
Chris@16:     template <typename Range, typename Box>
Chris@16:     static inline void apply(Range const& range, Box& mbr)
Chris@16:     {
Chris@16:         assign_inverse(mbr);
Chris@16:         envelope_range_additional(range, mbr);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@101: 
Chris@101: struct envelope_multi_linestring
Chris@101: {
Chris@101:     template<typename MultiLinestring, typename Box>
Chris@101:     static inline void apply(MultiLinestring const& mp, Box& mbr)
Chris@101:     {
Chris@101:         assign_inverse(mbr);
Chris@101:         for (typename boost::range_iterator<MultiLinestring const>::type
Chris@101:                 it = mp.begin();
Chris@101:             it != mp.end();
Chris@101:             ++it)
Chris@101:         {
Chris@101:             envelope_range_additional(*it, mbr);
Chris@101:         }
Chris@101:     }
Chris@101: };
Chris@101: 
Chris@101: 
Chris@101: // version for multi_polygon: outer ring's of all polygons
Chris@101: struct envelope_multi_polygon
Chris@101: {
Chris@101:     template<typename MultiPolygon, typename Box>
Chris@101:     static inline void apply(MultiPolygon const& mp, Box& mbr)
Chris@101:     {
Chris@101:         assign_inverse(mbr);
Chris@101:         for (typename boost::range_const_iterator<MultiPolygon>::type
Chris@101:                 it = mp.begin();
Chris@101:             it != mp.end();
Chris@101:             ++it)
Chris@101:         {
Chris@101:             envelope_range_additional(exterior_ring(*it), mbr);
Chris@101:         }
Chris@101:     }
Chris@101: };
Chris@101: 
Chris@101: 
Chris@16: }} // namespace detail::envelope
Chris@16: #endif // DOXYGEN_NO_DETAIL
Chris@16: 
Chris@16: #ifndef DOXYGEN_NO_DISPATCH
Chris@16: namespace dispatch
Chris@16: {
Chris@16: 
Chris@16: 
Chris@16: template
Chris@16: <
Chris@16:     typename Geometry,
Chris@16:     typename Tag = typename tag<Geometry>::type
Chris@16: >
Chris@16: struct envelope: not_implemented<Tag>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Point>
Chris@16: struct envelope<Point, point_tag>
Chris@16:     : detail::envelope::envelope_expand_one
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Box>
Chris@16: struct envelope<Box, box_tag>
Chris@16:     : detail::envelope::envelope_expand_one
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Segment>
Chris@16: struct envelope<Segment, segment_tag>
Chris@16:     : detail::envelope::envelope_expand_one
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Linestring>
Chris@16: struct envelope<Linestring, linestring_tag>
Chris@16:     : detail::envelope::envelope_range
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Ring>
Chris@16: struct envelope<Ring, ring_tag>
Chris@16:     : detail::envelope::envelope_range
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Polygon>
Chris@16: struct envelope<Polygon, polygon_tag>
Chris@16:     : detail::envelope::envelope_range
Chris@16: {
Chris@16:     template <typename Box>
Chris@16:     static inline void apply(Polygon const& poly, Box& mbr)
Chris@16:     {
Chris@16:         // For polygon, inspecting outer ring is sufficient
Chris@16:         detail::envelope::envelope_range::apply(exterior_ring(poly), mbr);
Chris@16:     }
Chris@16: 
Chris@16: };
Chris@16: 
Chris@16: 
Chris@101: template <typename Multi>
Chris@101: struct envelope<Multi, multi_point_tag>
Chris@101:     : detail::envelope::envelope_range
Chris@101: {};
Chris@101: 
Chris@101: 
Chris@101: template <typename Multi>
Chris@101: struct envelope<Multi, multi_linestring_tag>
Chris@101:     : detail::envelope::envelope_multi_linestring
Chris@101: {};
Chris@101: 
Chris@101: 
Chris@101: template <typename Multi>
Chris@101: struct envelope<Multi, multi_polygon_tag>
Chris@101:     : detail::envelope::envelope_multi_polygon
Chris@101: {};
Chris@101: 
Chris@101: 
Chris@16: } // namespace dispatch
Chris@16: #endif
Chris@16: 
Chris@16: 
Chris@101: namespace resolve_variant {
Chris@101: 
Chris@101: template <typename Geometry>
Chris@101: struct envelope
Chris@101: {
Chris@101:     template <typename Box>
Chris@101:     static inline void apply(Geometry const& geometry, Box& box)
Chris@101:     {
Chris@101:         concept::check<Geometry const>();
Chris@101:         concept::check<Box>();
Chris@101: 
Chris@101:         dispatch::envelope<Geometry>::apply(geometry, box);
Chris@101:     }
Chris@101: };
Chris@101: 
Chris@101: template <BOOST_VARIANT_ENUM_PARAMS(typename T)>
Chris@101: struct envelope<boost::variant<BOOST_VARIANT_ENUM_PARAMS(T)> >
Chris@101: {
Chris@101:     template <typename Box>
Chris@101:     struct visitor: boost::static_visitor<void>
Chris@101:     {
Chris@101:         Box& m_box;
Chris@101: 
Chris@101:         visitor(Box& box): m_box(box) {}
Chris@101: 
Chris@101:         template <typename Geometry>
Chris@101:         void operator()(Geometry const& geometry) const
Chris@101:         {
Chris@101:             envelope<Geometry>::apply(geometry, m_box);
Chris@101:         }
Chris@101:     };
Chris@101: 
Chris@101:     template <typename Box>
Chris@101:     static inline void
Chris@101:     apply(boost::variant<BOOST_VARIANT_ENUM_PARAMS(T)> const& geometry,
Chris@101:           Box& box)
Chris@101:     {
Chris@101:         boost::apply_visitor(visitor<Box>(box), geometry);
Chris@101:     }
Chris@101: };
Chris@101: 
Chris@101: } // namespace resolve_variant
Chris@101: 
Chris@101: 
Chris@16: /*!
Chris@16: \brief \brief_calc{envelope}
Chris@16: \ingroup envelope
Chris@16: \details \details_calc{envelope,\det_envelope}.
Chris@16: \tparam Geometry \tparam_geometry
Chris@16: \tparam Box \tparam_box
Chris@16: \param geometry \param_geometry
Chris@16: \param mbr \param_box \param_set{envelope}
Chris@16: 
Chris@16: \qbk{[include reference/algorithms/envelope.qbk]}
Chris@16: \qbk{
Chris@16: [heading Example]
Chris@16: [envelope] [envelope_output]
Chris@16: }
Chris@16: */
Chris@16: template<typename Geometry, typename Box>
Chris@16: inline void envelope(Geometry const& geometry, Box& mbr)
Chris@16: {
Chris@101:     resolve_variant::envelope<Geometry>::apply(geometry, mbr);
Chris@16: }
Chris@16: 
Chris@16: 
Chris@16: /*!
Chris@16: \brief \brief_calc{envelope}
Chris@16: \ingroup envelope
Chris@16: \details \details_calc{return_envelope,\det_envelope}. \details_return{envelope}
Chris@16: \tparam Box \tparam_box
Chris@16: \tparam Geometry \tparam_geometry
Chris@16: \param geometry \param_geometry
Chris@16: \return \return_calc{envelope}
Chris@16: 
Chris@16: \qbk{[include reference/algorithms/envelope.qbk]}
Chris@16: \qbk{
Chris@16: [heading Example]
Chris@16: [return_envelope] [return_envelope_output]
Chris@16: }
Chris@16: */
Chris@16: template<typename Box, typename Geometry>
Chris@16: inline Box return_envelope(Geometry const& geometry)
Chris@16: {
Chris@16:     Box mbr;
Chris@101:     resolve_variant::envelope<Geometry>::apply(geometry, mbr);
Chris@16:     return mbr;
Chris@16: }
Chris@16: 
Chris@16: }} // namespace boost::geometry
Chris@16: 
Chris@16: #endif // BOOST_GEOMETRY_ALGORITHMS_ENVELOPE_HPP