Chris@16: // Boost.Geometry (aka GGL, Generic Geometry Library)
Chris@16: 
Chris@101: // Copyright (c) 2007-2014 Barend Gehrels, Amsterdam, the Netherlands.
Chris@101: // Copyright (c) 2008-2014 Bruno Lalande, Paris, France.
Chris@101: // Copyright (c) 2009-2014 Mateusz Loskot, London, UK.
Chris@101: 
Chris@101: // This file was modified by Oracle on 2014.
Chris@101: // Modifications copyright (c) 2014, Oracle and/or its affiliates.
Chris@101: 
Chris@101: // Contributed and/or modified by Menelaos Karavelas, on behalf of Oracle
Chris@101: // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
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_APPEND_HPP
Chris@16: #define BOOST_GEOMETRY_ALGORITHMS_APPEND_HPP
Chris@16: 
Chris@16: 
Chris@101: #include <boost/range.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@101: 
Chris@16: #include <boost/geometry/algorithms/num_interior_rings.hpp>
Chris@16: #include <boost/geometry/algorithms/detail/convert_point_to_point.hpp>
Chris@16: #include <boost/geometry/core/access.hpp>
Chris@16: #include <boost/geometry/core/mutable_range.hpp>
Chris@16: #include <boost/geometry/core/point_type.hpp>
Chris@16: #include <boost/geometry/core/tags.hpp>
Chris@16: #include <boost/geometry/geometries/concepts/check.hpp>
Chris@16: #include <boost/geometry/geometries/variant.hpp>
Chris@101: #include <boost/geometry/util/range.hpp>
Chris@16: 
Chris@16: 
Chris@16: namespace boost { namespace geometry
Chris@16: {
Chris@16: 
Chris@16: 
Chris@16: #ifndef DOXYGEN_NO_DETAIL
Chris@16: namespace detail { namespace append
Chris@16: {
Chris@16: 
Chris@16: template <typename Geometry, typename Point>
Chris@16: struct append_no_action
Chris@16: {
Chris@16:     static inline void apply(Geometry& , Point const& ,
Chris@16:                 int = 0, int = 0)
Chris@16:     {
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: template <typename Geometry, typename Point>
Chris@16: struct append_point
Chris@16: {
Chris@16:     static inline void apply(Geometry& geometry, Point const& point,
Chris@16:                 int = 0, int = 0)
Chris@16:     {
Chris@16:         typename geometry::point_type<Geometry>::type copy;
Chris@16:         geometry::detail::conversion::convert_point_to_point(point, copy);
Chris@16:         traits::push_back<Geometry>::apply(geometry, copy);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: template <typename Geometry, typename Range>
Chris@16: struct append_range
Chris@16: {
Chris@16:     typedef typename boost::range_value<Range>::type point_type;
Chris@16: 
Chris@16:     static inline void apply(Geometry& geometry, Range const& range,
Chris@16:                 int = 0, int = 0)
Chris@16:     {
Chris@16:         for (typename boost::range_iterator<Range const>::type
Chris@16:             it = boost::begin(range);
Chris@16:              it != boost::end(range);
Chris@16:              ++it)
Chris@16:         {
Chris@16:             append_point<Geometry, point_type>::apply(geometry, *it);
Chris@16:         }
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: template <typename Polygon, typename Point>
Chris@16: struct point_to_polygon
Chris@16: {
Chris@16:     typedef typename ring_type<Polygon>::type ring_type;
Chris@16: 
Chris@16:     static inline void apply(Polygon& polygon, Point const& point,
Chris@16:                 int ring_index, int = 0)
Chris@16:     {
Chris@16:         if (ring_index == -1)
Chris@16:         {
Chris@16:             append_point<ring_type, Point>::apply(
Chris@16:                         exterior_ring(polygon), point);
Chris@16:         }
Chris@16:         else if (ring_index < int(num_interior_rings(polygon)))
Chris@16:         {
Chris@16:             append_point<ring_type, Point>::apply(
Chris@101:                         range::at(interior_rings(polygon), ring_index), point);
Chris@16:         }
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: template <typename Polygon, typename Range>
Chris@16: struct range_to_polygon
Chris@16: {
Chris@16:     typedef typename ring_type<Polygon>::type ring_type;
Chris@16: 
Chris@16:     static inline void apply(Polygon& polygon, Range const& range,
Chris@101:                 int ring_index, int = 0)
Chris@16:     {
Chris@16:         if (ring_index == -1)
Chris@16:         {
Chris@16:             append_range<ring_type, Range>::apply(
Chris@16:                         exterior_ring(polygon), range);
Chris@16:         }
Chris@16:         else if (ring_index < int(num_interior_rings(polygon)))
Chris@16:         {
Chris@16:             append_range<ring_type, Range>::apply(
Chris@101:                         range::at(interior_rings(polygon), ring_index), range);
Chris@16:         }
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: }} // namespace detail::append
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: namespace splitted_dispatch
Chris@16: {
Chris@16: 
Chris@16: template <typename Tag, typename Geometry, typename Point>
Chris@16: struct append_point
Chris@16:     : detail::append::append_no_action<Geometry, Point>
Chris@16: {};
Chris@16: 
Chris@16: template <typename Geometry, typename Point>
Chris@16: struct append_point<linestring_tag, Geometry, Point>
Chris@16:     : detail::append::append_point<Geometry, Point>
Chris@16: {};
Chris@16: 
Chris@16: template <typename Geometry, typename Point>
Chris@16: struct append_point<ring_tag, Geometry, Point>
Chris@16:     : detail::append::append_point<Geometry, Point>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Polygon, typename Point>
Chris@16: struct append_point<polygon_tag, Polygon, Point>
Chris@16:         : detail::append::point_to_polygon<Polygon, Point>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Tag, typename Geometry, typename Range>
Chris@16: struct append_range
Chris@16:     : detail::append::append_no_action<Geometry, Range>
Chris@16: {};
Chris@16: 
Chris@16: template <typename Geometry, typename Range>
Chris@16: struct append_range<linestring_tag, Geometry, Range>
Chris@16:     : detail::append::append_range<Geometry, Range>
Chris@16: {};
Chris@16: 
Chris@16: template <typename Geometry, typename Range>
Chris@16: struct append_range<ring_tag, Geometry, Range>
Chris@16:     : detail::append::append_range<Geometry, Range>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Polygon, typename Range>
Chris@16: struct append_range<polygon_tag, Polygon, Range>
Chris@16:         : detail::append::range_to_polygon<Polygon, Range>
Chris@16: {};
Chris@16: 
Chris@16: } // namespace splitted_dispatch
Chris@16: 
Chris@16: 
Chris@16: // Default: append a range (or linestring or ring or whatever) to any geometry
Chris@16: template
Chris@16: <
Chris@16:     typename Geometry, typename RangeOrPoint,
Chris@16:     typename TagRangeOrPoint = typename tag<RangeOrPoint>::type
Chris@16: >
Chris@16: struct append
Chris@16:     : splitted_dispatch::append_range<typename tag<Geometry>::type, Geometry, RangeOrPoint>
Chris@16: {};
Chris@16: 
Chris@16: // Specialization for point to append a point to any geometry
Chris@16: template <typename Geometry, typename RangeOrPoint>
Chris@16: struct append<Geometry, RangeOrPoint, point_tag>
Chris@16:     : splitted_dispatch::append_point<typename tag<Geometry>::type, Geometry, RangeOrPoint>
Chris@16: {};
Chris@16: 
Chris@101: } // namespace dispatch
Chris@101: #endif // DOXYGEN_NO_DISPATCH
Chris@101: 
Chris@101: #ifndef DOXYGEN_NO_DETAIL
Chris@101: namespace detail { namespace append
Chris@101: {
Chris@101: 
Chris@101: template <typename MultiGeometry, typename RangeOrPoint>
Chris@101: struct append_to_multigeometry
Chris@101: {
Chris@101:     static inline void apply(MultiGeometry& multigeometry,
Chris@101:                              RangeOrPoint const& range_or_point,
Chris@101:                              int ring_index, int multi_index)
Chris@101:     {
Chris@101: 
Chris@101:         dispatch::append
Chris@101:             <
Chris@101:                 typename boost::range_value<MultiGeometry>::type,
Chris@101:                 RangeOrPoint
Chris@101:             >::apply(range::at(multigeometry, multi_index), range_or_point, ring_index);
Chris@101:     }
Chris@101: };
Chris@101: 
Chris@101: }} // namespace detail::append
Chris@101: #endif // DOXYGEN_NO_DETAIL
Chris@101: 
Chris@101: #ifndef DOXYGEN_NO_DISPATCH
Chris@101: namespace dispatch
Chris@101: {
Chris@101: 
Chris@101: namespace splitted_dispatch
Chris@101: {
Chris@101: 
Chris@101: template <typename Geometry, typename Point>
Chris@101: struct append_point<multi_point_tag, Geometry, Point>
Chris@101:     : detail::append::append_point<Geometry, Point>
Chris@101: {};
Chris@101: 
Chris@101: template <typename Geometry, typename Range>
Chris@101: struct append_range<multi_point_tag, Geometry, Range>
Chris@101:     : detail::append::append_range<Geometry, Range>
Chris@101: {};
Chris@101: 
Chris@101: template <typename MultiGeometry, typename RangeOrPoint>
Chris@101: struct append_point<multi_linestring_tag, MultiGeometry, RangeOrPoint>
Chris@101:     : detail::append::append_to_multigeometry<MultiGeometry, RangeOrPoint>
Chris@101: {};
Chris@101: 
Chris@101: template <typename MultiGeometry, typename RangeOrPoint>
Chris@101: struct append_range<multi_linestring_tag, MultiGeometry, RangeOrPoint>
Chris@101:     : detail::append::append_to_multigeometry<MultiGeometry, RangeOrPoint>
Chris@101: {};
Chris@101: 
Chris@101: template <typename MultiGeometry, typename RangeOrPoint>
Chris@101: struct append_point<multi_polygon_tag, MultiGeometry, RangeOrPoint>
Chris@101:     : detail::append::append_to_multigeometry<MultiGeometry, RangeOrPoint>
Chris@101: {};
Chris@101: 
Chris@101: template <typename MultiGeometry, typename RangeOrPoint>
Chris@101: struct append_range<multi_polygon_tag, MultiGeometry, RangeOrPoint>
Chris@101:     : detail::append::append_to_multigeometry<MultiGeometry, RangeOrPoint>
Chris@101: {};
Chris@101: 
Chris@101: } // namespace splitted_dispatch
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 Geometry>
Chris@101: struct append
Chris@16: {
Chris@16:     template <typename RangeOrPoint>
Chris@16:     static inline void apply(Geometry& geometry,
Chris@16:                              RangeOrPoint const& range_or_point,
Chris@16:                              int ring_index,
Chris@16:                              int multi_index)
Chris@16:     {
Chris@16:         concept::check<Geometry>();
Chris@101:         dispatch::append<Geometry, RangeOrPoint>::apply(geometry,
Chris@101:                                                         range_or_point,
Chris@101:                                                         ring_index,
Chris@101:                                                         multi_index);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: template <BOOST_VARIANT_ENUM_PARAMS(typename T)>
Chris@101: struct append<boost::variant<BOOST_VARIANT_ENUM_PARAMS(T)> >
Chris@16: {
Chris@16:     template <typename RangeOrPoint>
Chris@16:     struct visitor: boost::static_visitor<void>
Chris@16:     {
Chris@16:         RangeOrPoint const& m_range_or_point;
Chris@16:         int m_ring_index;
Chris@16:         int m_multi_index;
Chris@101: 
Chris@16:         visitor(RangeOrPoint const& range_or_point,
Chris@16:                 int ring_index,
Chris@16:                 int multi_index):
Chris@16:             m_range_or_point(range_or_point),
Chris@16:             m_ring_index(ring_index),
Chris@16:             m_multi_index(multi_index)
Chris@16:         {}
Chris@101: 
Chris@16:         template <typename Geometry>
Chris@16:         void operator()(Geometry& geometry) const
Chris@16:         {
Chris@101:             append<Geometry>::apply(geometry,
Chris@101:                                     m_range_or_point,
Chris@101:                                     m_ring_index,
Chris@101:                                     m_multi_index);
Chris@16:         }
Chris@16:     };
Chris@16: 
Chris@16:     template <typename RangeOrPoint>
Chris@16:     static inline void apply(boost::variant<BOOST_VARIANT_ENUM_PARAMS(T)>& variant_geometry,
Chris@16:                              RangeOrPoint const& range_or_point,
Chris@16:                              int ring_index,
Chris@16:                              int multi_index)
Chris@16:     {
Chris@16:         apply_visitor(
Chris@16:             visitor<RangeOrPoint>(
Chris@16:                 range_or_point,
Chris@16:                 ring_index,
Chris@16:                 multi_index
Chris@16:             ),
Chris@16:             variant_geometry
Chris@16:         );
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@101: } // namespace resolve_variant;
Chris@16: 
Chris@16: 
Chris@16: /*!
Chris@16: \brief Appends one or more points to a linestring, ring, polygon, multi-geometry
Chris@16: \ingroup append
Chris@16: \tparam Geometry \tparam_geometry
Chris@16: \tparam RangeOrPoint Either a range or a point, fullfilling Boost.Range concept or Boost.Geometry Point Concept
Chris@16: \param geometry \param_geometry
Chris@16: \param range_or_point The point or range to add
Chris@16: \param ring_index The index of the ring in case of a polygon:
Chris@16:     exterior ring (-1, the default) or  interior ring index
Chris@101: \param multi_index The index of the geometry to which the points are appended
Chris@16: 
Chris@16: \qbk{[include reference/algorithms/append.qbk]}
Chris@16: }
Chris@16:  */
Chris@16: template <typename Geometry, typename RangeOrPoint>
Chris@16: inline void append(Geometry& geometry, RangeOrPoint const& range_or_point,
Chris@16:                    int ring_index = -1, int multi_index = 0)
Chris@16: {
Chris@101:     resolve_variant::append<Geometry>
Chris@101:                    ::apply(geometry, range_or_point, ring_index, multi_index);
Chris@16: }
Chris@16: 
Chris@16: 
Chris@16: }} // namespace boost::geometry
Chris@16: 
Chris@16: 
Chris@16: #endif // BOOST_GEOMETRY_ALGORITHMS_APPEND_HPP