123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
- // Use, modification and distribution is subject to the Boost Software License,
- // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
- // http://www.boost.org/LICENSE_1_0.txt)
- #ifndef BOOST_GEOMETRY_ALGORITHMS_DETAIL_OVERLAY_COPY_SEGMENT_POINT_HPP
- #define BOOST_GEOMETRY_ALGORITHMS_DETAIL_OVERLAY_COPY_SEGMENT_POINT_HPP
- #include <boost/array.hpp>
- #include <boost/mpl/assert.hpp>
- #include <boost/range.hpp>
- #include <boost/geometry/core/assert.hpp>
- #include <boost/geometry/core/ring_type.hpp>
- #include <boost/geometry/core/exterior_ring.hpp>
- #include <boost/geometry/core/interior_rings.hpp>
- #include <boost/geometry/core/tags.hpp>
- #include <boost/geometry/algorithms/convert.hpp>
- #include <boost/geometry/algorithms/detail/signed_size_type.hpp>
- #include <boost/geometry/geometries/concepts/check.hpp>
- #include <boost/geometry/util/range.hpp>
- #include <boost/geometry/iterators/ever_circling_iterator.hpp>
- #include <boost/geometry/views/closeable_view.hpp>
- #include <boost/geometry/views/reversible_view.hpp>
- namespace boost { namespace geometry
- {
- #ifndef DOXYGEN_NO_DETAIL
- namespace detail { namespace copy_segments
- {
- template <typename Range, bool Reverse, typename SegmentIdentifier, typename PointOut>
- struct copy_segment_point_range
- {
- static inline bool apply(Range const& range,
- SegmentIdentifier const& seg_id, int offset,
- PointOut& point)
- {
- typedef typename closeable_view
- <
- Range const,
- closure<Range>::value
- >::type cview_type;
- typedef typename reversible_view
- <
- cview_type const,
- Reverse ? iterate_reverse : iterate_forward
- >::type rview_type;
- cview_type cview(range);
- rview_type view(cview);
- typedef typename boost::range_iterator<rview_type>::type iterator;
- geometry::ever_circling_iterator<iterator> it(boost::begin(view), boost::end(view),
- boost::begin(view) + seg_id.segment_index, true);
- for (signed_size_type i = 0; i < offset; ++i, ++it)
- {
- }
- geometry::convert(*it, point);
- return true;
- }
- };
- template <typename Polygon, bool Reverse, typename SegmentIdentifier, typename PointOut>
- struct copy_segment_point_polygon
- {
- static inline bool apply(Polygon const& polygon,
- SegmentIdentifier const& seg_id, int offset,
- PointOut& point)
- {
- // Call ring-version with the right ring
- return copy_segment_point_range
- <
- typename geometry::ring_type<Polygon>::type,
- Reverse,
- SegmentIdentifier,
- PointOut
- >::apply
- (
- seg_id.ring_index < 0
- ? geometry::exterior_ring(polygon)
- : range::at(geometry::interior_rings(polygon), seg_id.ring_index),
- seg_id, offset,
- point
- );
- }
- };
- template <typename Box, bool Reverse, typename SegmentIdentifier, typename PointOut>
- struct copy_segment_point_box
- {
- static inline bool apply(Box const& box,
- SegmentIdentifier const& seg_id, int offset,
- PointOut& point)
- {
- signed_size_type index = seg_id.segment_index;
- for (int i = 0; i < offset; i++)
- {
- index++;
- }
- boost::array<typename point_type<Box>::type, 4> bp;
- assign_box_corners_oriented<Reverse>(box, bp);
- point = bp[index % 4];
- return true;
- }
- };
- template
- <
- typename MultiGeometry,
- typename SegmentIdentifier,
- typename PointOut,
- typename Policy
- >
- struct copy_segment_point_multi
- {
- static inline bool apply(MultiGeometry const& multi,
- SegmentIdentifier const& seg_id, int offset,
- PointOut& point)
- {
- BOOST_GEOMETRY_ASSERT
- (
- seg_id.multi_index >= 0
- && seg_id.multi_index < int(boost::size(multi))
- );
- // Call the single-version
- return Policy::apply(range::at(multi, seg_id.multi_index), seg_id, offset, point);
- }
- };
- }} // namespace detail::copy_segments
- #endif // DOXYGEN_NO_DETAIL
- #ifndef DOXYGEN_NO_DISPATCH
- namespace dispatch
- {
- template
- <
- typename Tag,
- typename GeometryIn,
- bool Reverse,
- typename SegmentIdentifier,
- typename PointOut
- >
- struct copy_segment_point
- {
- BOOST_MPL_ASSERT_MSG
- (
- false, NOT_OR_NOT_YET_IMPLEMENTED_FOR_THIS_GEOMETRY_TYPE
- , (types<GeometryIn>)
- );
- };
- template <typename LineString, bool Reverse, typename SegmentIdentifier, typename PointOut>
- struct copy_segment_point<linestring_tag, LineString, Reverse, SegmentIdentifier, PointOut>
- : detail::copy_segments::copy_segment_point_range
- <
- LineString, Reverse, SegmentIdentifier, PointOut
- >
- {};
- template <typename Ring, bool Reverse, typename SegmentIdentifier, typename PointOut>
- struct copy_segment_point<ring_tag, Ring, Reverse, SegmentIdentifier, PointOut>
- : detail::copy_segments::copy_segment_point_range
- <
- Ring, Reverse, SegmentIdentifier, PointOut
- >
- {};
- template <typename Polygon, bool Reverse, typename SegmentIdentifier, typename PointOut>
- struct copy_segment_point<polygon_tag, Polygon, Reverse, SegmentIdentifier, PointOut>
- : detail::copy_segments::copy_segment_point_polygon
- <
- Polygon, Reverse, SegmentIdentifier, PointOut
- >
- {};
- template <typename Box, bool Reverse, typename SegmentIdentifier, typename PointOut>
- struct copy_segment_point<box_tag, Box, Reverse, SegmentIdentifier, PointOut>
- : detail::copy_segments::copy_segment_point_box
- <
- Box, Reverse, SegmentIdentifier, PointOut
- >
- {};
- template
- <
- typename MultiGeometry,
- bool Reverse,
- typename SegmentIdentifier,
- typename PointOut
- >
- struct copy_segment_point
- <
- multi_polygon_tag,
- MultiGeometry,
- Reverse,
- SegmentIdentifier,
- PointOut
- >
- : detail::copy_segments::copy_segment_point_multi
- <
- MultiGeometry,
- SegmentIdentifier,
- PointOut,
- detail::copy_segments::copy_segment_point_polygon
- <
- typename boost::range_value<MultiGeometry>::type,
- Reverse,
- SegmentIdentifier,
- PointOut
- >
- >
- {};
- template
- <
- typename MultiGeometry,
- bool Reverse,
- typename SegmentIdentifier,
- typename PointOut
- >
- struct copy_segment_point
- <
- multi_linestring_tag,
- MultiGeometry,
- Reverse,
- SegmentIdentifier,
- PointOut
- >
- : detail::copy_segments::copy_segment_point_multi
- <
- MultiGeometry,
- SegmentIdentifier,
- PointOut,
- detail::copy_segments::copy_segment_point_range
- <
- typename boost::range_value<MultiGeometry>::type,
- Reverse,
- SegmentIdentifier,
- PointOut
- >
- >
- {};
- } // namespace dispatch
- #endif // DOXYGEN_NO_DISPATCH
- /*!
- \brief Helper function, copies a point from a segment
- \ingroup overlay
- */
- template<bool Reverse, typename Geometry, typename SegmentIdentifier, typename PointOut>
- inline bool copy_segment_point(Geometry const& geometry,
- SegmentIdentifier const& seg_id, int offset,
- PointOut& point_out)
- {
- concepts::check<Geometry const>();
- return dispatch::copy_segment_point
- <
- typename tag<Geometry>::type,
- Geometry,
- Reverse,
- SegmentIdentifier,
- PointOut
- >::apply(geometry, seg_id, offset, point_out);
- }
- /*!
- \brief Helper function, to avoid the same construct several times,
- copies a point, based on a source-index and two geometries
- \ingroup overlay
- */
- template
- <
- bool Reverse1, bool Reverse2,
- typename Geometry1, typename Geometry2,
- typename SegmentIdentifier,
- typename PointOut
- >
- inline bool copy_segment_point(Geometry1 const& geometry1, Geometry2 const& geometry2,
- SegmentIdentifier const& seg_id, int offset,
- PointOut& point_out)
- {
- concepts::check<Geometry1 const>();
- concepts::check<Geometry2 const>();
- BOOST_GEOMETRY_ASSERT(seg_id.source_index == 0 || seg_id.source_index == 1);
- if (seg_id.source_index == 0)
- {
- return dispatch::copy_segment_point
- <
- typename tag<Geometry1>::type,
- Geometry1,
- Reverse1,
- SegmentIdentifier,
- PointOut
- >::apply(geometry1, seg_id, offset, point_out);
- }
- else if (seg_id.source_index == 1)
- {
- return dispatch::copy_segment_point
- <
- typename tag<Geometry2>::type,
- Geometry2,
- Reverse2,
- SegmentIdentifier,
- PointOut
- >::apply(geometry2, seg_id, offset, point_out);
- }
- // Exception?
- return false;
- }
- /*!
- \brief Helper function, to avoid the same construct several times,
- copies a point, based on a source-index and two geometries
- \ingroup overlay
- */
- template
- <
- bool Reverse1, bool Reverse2,
- typename Geometry1, typename Geometry2,
- typename SegmentIdentifier,
- typename PointOut
- >
- inline bool copy_segment_points(Geometry1 const& geometry1, Geometry2 const& geometry2,
- SegmentIdentifier const& seg_id,
- PointOut& point1, PointOut& point2)
- {
- concepts::check<Geometry1 const>();
- concepts::check<Geometry2 const>();
- return copy_segment_point<Reverse1, Reverse2>(geometry1, geometry2, seg_id, 0, point1)
- && copy_segment_point<Reverse1, Reverse2>(geometry1, geometry2, seg_id, 1, point2);
- }
- /*!
- \brief Helper function, copies three points: two from the specified segment
- (from, to) and the next one
- \ingroup overlay
- */
- template
- <
- bool Reverse1, bool Reverse2,
- typename Geometry1, typename Geometry2,
- typename SegmentIdentifier,
- typename PointOut
- >
- inline bool copy_segment_points(Geometry1 const& geometry1, Geometry2 const& geometry2,
- SegmentIdentifier const& seg_id,
- PointOut& point1, PointOut& point2, PointOut& point3)
- {
- concepts::check<Geometry1 const>();
- concepts::check<Geometry2 const>();
- return copy_segment_point<Reverse1, Reverse2>(geometry1, geometry2, seg_id, 0, point1)
- && copy_segment_point<Reverse1, Reverse2>(geometry1, geometry2, seg_id, 1, point2)
- && copy_segment_point<Reverse1, Reverse2>(geometry1, geometry2, seg_id, 2, point3);
- }
- }} // namespace boost::geometry
- #endif // BOOST_GEOMETRY_ALGORITHMS_DETAIL_OVERLAY_COPY_SEGMENT_POINT_HPP
|