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@101: // Copyright (c) 2014 Adam Wulkiewicz, Lodz, Poland.
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_CONVERT_HPP
Chris@16: #define BOOST_GEOMETRY_ALGORITHMS_CONVERT_HPP
Chris@16: 
Chris@16: 
Chris@16: #include <cstddef>
Chris@16: 
Chris@16: #include <boost/numeric/conversion/cast.hpp>
Chris@16: #include <boost/range.hpp>
Chris@16: #include <boost/type_traits/is_array.hpp>
Chris@101: #include <boost/type_traits/remove_reference.hpp>
Chris@101: 
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/arithmetic/arithmetic.hpp>
Chris@16: #include <boost/geometry/algorithms/not_implemented.hpp>
Chris@16: #include <boost/geometry/algorithms/append.hpp>
Chris@16: #include <boost/geometry/algorithms/clear.hpp>
Chris@16: #include <boost/geometry/algorithms/for_each.hpp>
Chris@16: #include <boost/geometry/algorithms/detail/assign_values.hpp>
Chris@16: #include <boost/geometry/algorithms/detail/assign_box_corners.hpp>
Chris@16: #include <boost/geometry/algorithms/detail/assign_indexed_point.hpp>
Chris@16: #include <boost/geometry/algorithms/detail/convert_point_to_point.hpp>
Chris@16: #include <boost/geometry/algorithms/detail/convert_indexed_to_indexed.hpp>
Chris@101: #include <boost/geometry/algorithms/detail/interior_iterator.hpp>
Chris@16: 
Chris@16: #include <boost/geometry/views/closeable_view.hpp>
Chris@16: #include <boost/geometry/views/reversible_view.hpp>
Chris@16: 
Chris@101: #include <boost/geometry/util/range.hpp>
Chris@101: 
Chris@16: #include <boost/geometry/core/cs.hpp>
Chris@16: #include <boost/geometry/core/closure.hpp>
Chris@16: #include <boost/geometry/core/point_order.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: // Silence warning C4127: conditional expression is constant
Chris@16: // Silence warning C4512: assignment operator could not be generated
Chris@16: #if defined(_MSC_VER)
Chris@101: #pragma warning(push)
Chris@16: #pragma warning(disable : 4127 4512)
Chris@16: #endif
Chris@16: 
Chris@16: 
Chris@16: #ifndef DOXYGEN_NO_DETAIL
Chris@16: namespace detail { namespace conversion
Chris@16: {
Chris@16: 
Chris@16: template
Chris@16: <
Chris@16:     typename Point,
Chris@16:     typename Box,
Chris@16:     std::size_t Index,
Chris@16:     std::size_t Dimension,
Chris@16:     std::size_t DimensionCount
Chris@16: >
Chris@16: struct point_to_box
Chris@16: {
Chris@16:     static inline void apply(Point const& point, Box& box)
Chris@16:     {
Chris@16:         typedef typename coordinate_type<Box>::type coordinate_type;
Chris@16: 
Chris@16:         set<Index, Dimension>(box,
Chris@16:                 boost::numeric_cast<coordinate_type>(get<Dimension>(point)));
Chris@16:         point_to_box
Chris@16:             <
Chris@16:                 Point, Box,
Chris@16:                 Index, Dimension + 1, DimensionCount
Chris@16:             >::apply(point, box);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: template
Chris@16: <
Chris@16:     typename Point,
Chris@16:     typename Box,
Chris@16:     std::size_t Index,
Chris@16:     std::size_t DimensionCount
Chris@16: >
Chris@16: struct point_to_box<Point, Box, Index, DimensionCount, DimensionCount>
Chris@16: {
Chris@16:     static inline void apply(Point const& , Box& )
Chris@16:     {}
Chris@16: };
Chris@16: 
Chris@16: template <typename Box, typename Range, bool Close, bool Reverse>
Chris@16: struct box_to_range
Chris@16: {
Chris@16:     static inline void apply(Box const& box, Range& range)
Chris@16:     {
Chris@16:         traits::resize<Range>::apply(range, Close ? 5 : 4);
Chris@16:         assign_box_corners_oriented<Reverse>(box, range);
Chris@16:         if (Close)
Chris@16:         {
Chris@101:             range::at(range, 4) = range::at(range, 0);
Chris@16:         }
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: template <typename Segment, typename Range>
Chris@16: struct segment_to_range
Chris@16: {
Chris@16:     static inline void apply(Segment const& segment, Range& range)
Chris@16:     {
Chris@16:         traits::resize<Range>::apply(range, 2);
Chris@16: 
Chris@16:         typename boost::range_iterator<Range>::type it = boost::begin(range);
Chris@16: 
Chris@16:         assign_point_from_index<0>(segment, *it);
Chris@16:         ++it;
Chris@16:         assign_point_from_index<1>(segment, *it);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@101: template
Chris@16: <
Chris@101:     typename Range1,
Chris@101:     typename Range2,
Chris@16:     bool Reverse = false
Chris@16: >
Chris@16: struct range_to_range
Chris@16: {
Chris@16:     typedef typename reversible_view
Chris@16:         <
Chris@101:             Range1 const,
Chris@16:             Reverse ? iterate_reverse : iterate_forward
Chris@16:         >::type rview_type;
Chris@16:     typedef typename closeable_view
Chris@16:         <
Chris@101:             rview_type const,
Chris@16:             geometry::closure<Range1>::value
Chris@16:         >::type view_type;
Chris@16: 
Chris@16:     static inline void apply(Range1 const& source, Range2& destination)
Chris@16:     {
Chris@16:         geometry::clear(destination);
Chris@16: 
Chris@16:         rview_type rview(source);
Chris@16: 
Chris@16:         // We consider input always as closed, and skip last
Chris@16:         // point for open output.
Chris@16:         view_type view(rview);
Chris@16: 
Chris@101:         typedef typename boost::range_size<Range1>::type size_type;
Chris@101:         size_type n = boost::size(view);
Chris@16:         if (geometry::closure<Range2>::value == geometry::open)
Chris@16:         {
Chris@16:             n--;
Chris@16:         }
Chris@16: 
Chris@101:         // If size == 0 && geometry::open <=> n = numeric_limits<size_type>::max()
Chris@101:         // but ok, sice below it == end()
Chris@101: 
Chris@101:         size_type i = 0;
Chris@16:         for (typename boost::range_iterator<view_type const>::type it
Chris@16:             = boost::begin(view);
Chris@16:             it != boost::end(view) && i < n;
Chris@16:             ++it, ++i)
Chris@16:         {
Chris@16:             geometry::append(destination, *it);
Chris@16:         }
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: template <typename Polygon1, typename Polygon2>
Chris@16: struct polygon_to_polygon
Chris@16: {
Chris@16:     typedef range_to_range
Chris@16:         <
Chris@101:             typename geometry::ring_type<Polygon1>::type,
Chris@16:             typename geometry::ring_type<Polygon2>::type,
Chris@16:             geometry::point_order<Polygon1>::value
Chris@16:                 != geometry::point_order<Polygon2>::value
Chris@16:         > per_ring;
Chris@16: 
Chris@16:     static inline void apply(Polygon1 const& source, Polygon2& destination)
Chris@16:     {
Chris@16:         // Clearing managed per ring, and in the resizing of interior rings
Chris@16: 
Chris@101:         per_ring::apply(geometry::exterior_ring(source),
Chris@16:             geometry::exterior_ring(destination));
Chris@16: 
Chris@16:         // Container should be resizeable
Chris@16:         traits::resize
Chris@16:             <
Chris@16:                 typename boost::remove_reference
Chris@16:                 <
Chris@16:                     typename traits::interior_mutable_type<Polygon2>::type
Chris@16:                 >::type
Chris@16:             >::apply(interior_rings(destination), num_interior_rings(source));
Chris@16: 
Chris@101:         typename interior_return_type<Polygon1 const>::type
Chris@101:             rings_source = interior_rings(source);
Chris@101:         typename interior_return_type<Polygon2>::type
Chris@101:             rings_dest = interior_rings(destination);
Chris@16: 
Chris@101:         typename detail::interior_iterator<Polygon1 const>::type
Chris@101:             it_source = boost::begin(rings_source);
Chris@101:         typename detail::interior_iterator<Polygon2>::type
Chris@101:             it_dest = boost::begin(rings_dest);
Chris@16: 
Chris@16:         for ( ; it_source != boost::end(rings_source); ++it_source, ++it_dest)
Chris@16:         {
Chris@16:             per_ring::apply(*it_source, *it_dest);
Chris@16:         }
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@101: template <typename Single, typename Multi, typename Policy>
Chris@101: struct single_to_multi: private Policy
Chris@101: {
Chris@101:     static inline void apply(Single const& single, Multi& multi)
Chris@101:     {
Chris@101:         traits::resize<Multi>::apply(multi, 1);
Chris@101:         Policy::apply(single, *boost::begin(multi));
Chris@101:     }
Chris@101: };
Chris@101: 
Chris@101: 
Chris@101: 
Chris@101: template <typename Multi1, typename Multi2, typename Policy>
Chris@101: struct multi_to_multi: private Policy
Chris@101: {
Chris@101:     static inline void apply(Multi1 const& multi1, Multi2& multi2)
Chris@101:     {
Chris@101:         traits::resize<Multi2>::apply(multi2, boost::size(multi1));
Chris@101: 
Chris@101:         typename boost::range_iterator<Multi1 const>::type it1
Chris@101:                 = boost::begin(multi1);
Chris@101:         typename boost::range_iterator<Multi2>::type it2
Chris@101:                 = boost::begin(multi2);
Chris@101: 
Chris@101:         for (; it1 != boost::end(multi1); ++it1, ++it2)
Chris@101:         {
Chris@101:             Policy::apply(*it1, *it2);
Chris@101:         }
Chris@101:     }
Chris@101: };
Chris@101: 
Chris@16: 
Chris@16: }} // namespace detail::conversion
Chris@16: #endif // DOXYGEN_NO_DETAIL
Chris@16: 
Chris@16: 
Chris@16: #ifndef DOXYGEN_NO_DISPATCH
Chris@16: namespace dispatch
Chris@16: {
Chris@16: 
Chris@16: template
Chris@16: <
Chris@16:     typename Geometry1, typename Geometry2,
Chris@16:     typename Tag1 = typename tag_cast<typename tag<Geometry1>::type, multi_tag>::type,
Chris@16:     typename Tag2 = typename tag_cast<typename tag<Geometry2>::type, multi_tag>::type,
Chris@16:     std::size_t DimensionCount = dimension<Geometry1>::type::value,
Chris@16:     bool UseAssignment = boost::is_same<Geometry1, Geometry2>::value
Chris@16:                          && !boost::is_array<Geometry1>::value
Chris@16: >
Chris@101: struct convert: not_implemented<Tag1, Tag2, boost::mpl::int_<DimensionCount> >
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template
Chris@16: <
Chris@16:     typename Geometry1, typename Geometry2,
Chris@16:     typename Tag,
Chris@16:     std::size_t DimensionCount
Chris@16: >
Chris@16: struct convert<Geometry1, Geometry2, Tag, Tag, DimensionCount, true>
Chris@16: {
Chris@16:     // Same geometry type -> copy whole geometry
Chris@16:     static inline void apply(Geometry1 const& source, Geometry2& destination)
Chris@16:     {
Chris@16:         destination = source;
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: template
Chris@16: <
Chris@16:     typename Geometry1, typename Geometry2,
Chris@16:     std::size_t DimensionCount
Chris@16: >
Chris@16: struct convert<Geometry1, Geometry2, point_tag, point_tag, DimensionCount, false>
Chris@16:     : detail::conversion::point_to_point<Geometry1, Geometry2, 0, DimensionCount>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template
Chris@16: <
Chris@16:     typename Box1, typename Box2,
Chris@16:     std::size_t DimensionCount
Chris@16: >
Chris@16: struct convert<Box1, Box2, box_tag, box_tag, DimensionCount, false>
Chris@16:     : detail::conversion::indexed_to_indexed<Box1, Box2, 0, DimensionCount>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template
Chris@16: <
Chris@16:     typename Segment1, typename Segment2,
Chris@16:     std::size_t DimensionCount
Chris@16: >
Chris@16: struct convert<Segment1, Segment2, segment_tag, segment_tag, DimensionCount, false>
Chris@16:     : detail::conversion::indexed_to_indexed<Segment1, Segment2, 0, DimensionCount>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Segment, typename LineString, std::size_t DimensionCount>
Chris@16: struct convert<Segment, LineString, segment_tag, linestring_tag, DimensionCount, false>
Chris@16:     : detail::conversion::segment_to_range<Segment, LineString>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Ring1, typename Ring2, std::size_t DimensionCount>
Chris@16: struct convert<Ring1, Ring2, ring_tag, ring_tag, DimensionCount, false>
Chris@16:     : detail::conversion::range_to_range
Chris@101:         <
Chris@101:             Ring1,
Chris@16:             Ring2,
Chris@16:             geometry::point_order<Ring1>::value
Chris@16:                 != geometry::point_order<Ring2>::value
Chris@16:         >
Chris@16: {};
Chris@16: 
Chris@16: template <typename LineString1, typename LineString2, std::size_t DimensionCount>
Chris@16: struct convert<LineString1, LineString2, linestring_tag, linestring_tag, DimensionCount, false>
Chris@16:     : detail::conversion::range_to_range<LineString1, LineString2>
Chris@16: {};
Chris@16: 
Chris@16: template <typename Polygon1, typename Polygon2, std::size_t DimensionCount>
Chris@16: struct convert<Polygon1, Polygon2, polygon_tag, polygon_tag, DimensionCount, false>
Chris@16:     : detail::conversion::polygon_to_polygon<Polygon1, Polygon2>
Chris@16: {};
Chris@16: 
Chris@16: template <typename Box, typename Ring>
Chris@16: struct convert<Box, Ring, box_tag, ring_tag, 2, false>
Chris@16:     : detail::conversion::box_to_range
Chris@16:         <
Chris@101:             Box,
Chris@101:             Ring,
Chris@16:             geometry::closure<Ring>::value == closed,
Chris@16:             geometry::point_order<Ring>::value == counterclockwise
Chris@16:         >
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Box, typename Polygon>
Chris@16: struct convert<Box, Polygon, box_tag, polygon_tag, 2, false>
Chris@16: {
Chris@16:     static inline void apply(Box const& box, Polygon& polygon)
Chris@16:     {
Chris@16:         typedef typename ring_type<Polygon>::type ring_type;
Chris@16: 
Chris@16:         convert
Chris@16:             <
Chris@16:                 Box, ring_type,
Chris@16:                 box_tag, ring_tag,
Chris@16:                 2, false
Chris@16:             >::apply(box, exterior_ring(polygon));
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: template <typename Point, typename Box, std::size_t DimensionCount>
Chris@16: struct convert<Point, Box, point_tag, box_tag, DimensionCount, false>
Chris@16: {
Chris@16:     static inline void apply(Point const& point, Box& box)
Chris@16:     {
Chris@16:         detail::conversion::point_to_box
Chris@16:             <
Chris@16:                 Point, Box, min_corner, 0, DimensionCount
Chris@16:             >::apply(point, box);
Chris@16:         detail::conversion::point_to_box
Chris@16:             <
Chris@16:                 Point, Box, max_corner, 0, DimensionCount
Chris@16:             >::apply(point, box);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: template <typename Ring, typename Polygon, std::size_t DimensionCount>
Chris@16: struct convert<Ring, Polygon, ring_tag, polygon_tag, DimensionCount, false>
Chris@16: {
Chris@16:     static inline void apply(Ring const& ring, Polygon& polygon)
Chris@16:     {
Chris@16:         typedef typename ring_type<Polygon>::type ring_type;
Chris@16:         convert
Chris@16:             <
Chris@16:                 Ring, ring_type,
Chris@16:                 ring_tag, ring_tag,
Chris@16:                 DimensionCount, false
Chris@16:             >::apply(ring, exterior_ring(polygon));
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: template <typename Polygon, typename Ring, std::size_t DimensionCount>
Chris@16: struct convert<Polygon, Ring, polygon_tag, ring_tag, DimensionCount, false>
Chris@16: {
Chris@16:     static inline void apply(Polygon const& polygon, Ring& ring)
Chris@16:     {
Chris@16:         typedef typename ring_type<Polygon>::type ring_type;
Chris@16: 
Chris@16:         convert
Chris@16:             <
Chris@16:                 ring_type, Ring,
Chris@16:                 ring_tag, ring_tag,
Chris@16:                 DimensionCount, false
Chris@16:             >::apply(exterior_ring(polygon), ring);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@101: // Dispatch for multi <-> multi, specifying their single-version as policy.
Chris@101: // Note that, even if the multi-types are mutually different, their single
Chris@101: // version types might be the same and therefore we call boost::is_same again
Chris@101: 
Chris@101: template <typename Multi1, typename Multi2, std::size_t DimensionCount>
Chris@101: struct convert<Multi1, Multi2, multi_tag, multi_tag, DimensionCount, false>
Chris@101:     : detail::conversion::multi_to_multi
Chris@101:         <
Chris@101:             Multi1,
Chris@101:             Multi2,
Chris@101:             convert
Chris@101:                 <
Chris@101:                     typename boost::range_value<Multi1>::type,
Chris@101:                     typename boost::range_value<Multi2>::type,
Chris@101:                     typename single_tag_of
Chris@101:                                 <
Chris@101:                                     typename tag<Multi1>::type
Chris@101:                                 >::type,
Chris@101:                     typename single_tag_of
Chris@101:                                 <
Chris@101:                                     typename tag<Multi2>::type
Chris@101:                                 >::type,
Chris@101:                     DimensionCount
Chris@101:                 >
Chris@101:         >
Chris@101: {};
Chris@101: 
Chris@101: 
Chris@101: template <typename Single, typename Multi, typename SingleTag, std::size_t DimensionCount>
Chris@101: struct convert<Single, Multi, SingleTag, multi_tag, DimensionCount, false>
Chris@101:     : detail::conversion::single_to_multi
Chris@101:         <
Chris@101:             Single,
Chris@101:             Multi,
Chris@101:             convert
Chris@101:                 <
Chris@101:                     Single,
Chris@101:                     typename boost::range_value<Multi>::type,
Chris@101:                     typename tag<Single>::type,
Chris@101:                     typename single_tag_of
Chris@101:                                 <
Chris@101:                                     typename tag<Multi>::type
Chris@101:                                 >::type,
Chris@101:                     DimensionCount,
Chris@101:                     false
Chris@101:                 >
Chris@101:         >
Chris@101: {};
Chris@101: 
Chris@101: 
Chris@101: } // namespace dispatch
Chris@101: #endif // DOXYGEN_NO_DISPATCH
Chris@101: 
Chris@101: 
Chris@101: namespace resolve_variant {
Chris@101: 
Chris@16: template <typename Geometry1, typename Geometry2>
Chris@101: struct convert
Chris@16: {
Chris@16:     static inline void apply(Geometry1 const& geometry1, Geometry2& geometry2)
Chris@16:     {
Chris@16:         concept::check_concepts_and_equal_dimensions<Geometry1 const, Geometry2>();
Chris@101:         dispatch::convert<Geometry1, Geometry2>::apply(geometry1, geometry2);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: template <BOOST_VARIANT_ENUM_PARAMS(typename T), typename Geometry2>
Chris@101: struct convert<boost::variant<BOOST_VARIANT_ENUM_PARAMS(T)>, Geometry2>
Chris@16: {
Chris@16:     struct visitor: static_visitor<void>
Chris@16:     {
Chris@16:         Geometry2& m_geometry2;
Chris@16: 
Chris@16:         visitor(Geometry2& geometry2)
Chris@16:             : m_geometry2(geometry2)
Chris@16:         {}
Chris@16: 
Chris@16:         template <typename Geometry1>
Chris@16:         inline void operator()(Geometry1 const& geometry1) const
Chris@16:         {
Chris@101:             convert<Geometry1, Geometry2>::apply(geometry1, m_geometry2);
Chris@16:         }
Chris@16:     };
Chris@16: 
Chris@16:     static inline void apply(
Chris@16:         boost::variant<BOOST_VARIANT_ENUM_PARAMS(T)> const& geometry1,
Chris@16:         Geometry2& geometry2
Chris@16:     )
Chris@16:     {
Chris@16:         apply_visitor(visitor(geometry2), geometry1);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@101: }
Chris@16: 
Chris@16: 
Chris@16: /*!
Chris@16: \brief Converts one geometry to another geometry
Chris@16: \details The convert algorithm converts one geometry, e.g. a BOX, to another
Chris@16: geometry, e.g. a RING. This only works if it is possible and applicable.
Chris@101: If the point-order is different, or the closure is different between two
Chris@101: geometry types, it will be converted correctly by explicitly reversing the
Chris@16: points or closing or opening the polygon rings.
Chris@16: \ingroup convert
Chris@16: \tparam Geometry1 \tparam_geometry
Chris@16: \tparam Geometry2 \tparam_geometry
Chris@16: \param geometry1 \param_geometry (source)
Chris@16: \param geometry2 \param_geometry (target)
Chris@16: 
Chris@16: \qbk{[include reference/algorithms/convert.qbk]}
Chris@16:  */
Chris@16: template <typename Geometry1, typename Geometry2>
Chris@16: inline void convert(Geometry1 const& geometry1, Geometry2& geometry2)
Chris@16: {
Chris@101:     resolve_variant::convert<Geometry1, Geometry2>::apply(geometry1, geometry2);
Chris@16: }
Chris@16: 
Chris@16: #if defined(_MSC_VER)
Chris@16: #pragma warning(pop)
Chris@16: #endif
Chris@16: 
Chris@16: }} // namespace boost::geometry
Chris@16: 
Chris@16: #endif // BOOST_GEOMETRY_ALGORITHMS_CONVERT_HPP