123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
- // Copyright (c) 2008-2012 Bruno Lalande, Paris, France.
- // Copyright (c) 2009-2012 Mateusz Loskot, London, UK.
- // Copyright (c) 2014-2017 Adam Wulkiewicz, Lodz, Poland.
- // This file was modified by Oracle on 2017.
- // Modifications copyright (c) 2017 Oracle and/or its affiliates.
- // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
- // Parts of Boost.Geometry are redesigned from Geodan's Geographic Library
- // (geolib/GGL), copyright (c) 1995-2010 Geodan, 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_CORRECT_HPP
- #define BOOST_GEOMETRY_ALGORITHMS_CORRECT_HPP
- #include <algorithm>
- #include <cstddef>
- #include <functional>
- #include <boost/mpl/assert.hpp>
- #include <boost/range.hpp>
- #include <boost/type_traits/remove_reference.hpp>
- #include <boost/variant/apply_visitor.hpp>
- #include <boost/variant/static_visitor.hpp>
- #include <boost/variant/variant_fwd.hpp>
- #include <boost/geometry/algorithms/correct_closure.hpp>
- #include <boost/geometry/algorithms/detail/interior_iterator.hpp>
- #include <boost/geometry/core/closure.hpp>
- #include <boost/geometry/core/cs.hpp>
- #include <boost/geometry/core/exterior_ring.hpp>
- #include <boost/geometry/core/interior_rings.hpp>
- #include <boost/geometry/core/mutable_range.hpp>
- #include <boost/geometry/core/ring_type.hpp>
- #include <boost/geometry/core/tags.hpp>
- #include <boost/geometry/geometries/concepts/check.hpp>
- #include <boost/geometry/algorithms/area.hpp>
- #include <boost/geometry/algorithms/detail/multi_modify.hpp>
- #include <boost/geometry/util/order_as_direction.hpp>
- namespace boost { namespace geometry
- {
- // Silence warning C4127: conditional expression is constant
- #if defined(_MSC_VER)
- #pragma warning(push)
- #pragma warning(disable : 4127)
- #endif
- #ifndef DOXYGEN_NO_DETAIL
- namespace detail { namespace correct
- {
- template <typename Geometry>
- struct correct_nop
- {
- template <typename Strategy>
- static inline void apply(Geometry& , Strategy const& )
- {}
- };
- template <typename Box, std::size_t Dimension, std::size_t DimensionCount>
- struct correct_box_loop
- {
- typedef typename coordinate_type<Box>::type coordinate_type;
- static inline void apply(Box& box)
- {
- if (get<min_corner, Dimension>(box) > get<max_corner, Dimension>(box))
- {
- // Swap the coordinates
- coordinate_type max_value = get<min_corner, Dimension>(box);
- coordinate_type min_value = get<max_corner, Dimension>(box);
- set<min_corner, Dimension>(box, min_value);
- set<max_corner, Dimension>(box, max_value);
- }
- correct_box_loop
- <
- Box, Dimension + 1, DimensionCount
- >::apply(box);
- }
- };
- template <typename Box, std::size_t DimensionCount>
- struct correct_box_loop<Box, DimensionCount, DimensionCount>
- {
- static inline void apply(Box& )
- {}
- };
- // Correct a box: make min/max correct
- template <typename Box>
- struct correct_box
- {
- template <typename Strategy>
- static inline void apply(Box& box, Strategy const& )
- {
- // Currently only for Cartesian coordinates
- // (or spherical without crossing dateline)
- // Future version: adapt using strategies
- correct_box_loop
- <
- Box, 0, dimension<Box>::type::value
- >::apply(box);
- }
- };
- // Close a ring, if not closed
- template <typename Ring, template <typename> class Predicate>
- struct correct_ring
- {
- typedef typename point_type<Ring>::type point_type;
- typedef typename coordinate_type<Ring>::type coordinate_type;
- typedef detail::area::ring_area
- <
- order_as_direction<geometry::point_order<Ring>::value>::value,
- geometry::closure<Ring>::value
- > ring_area_type;
- template <typename Strategy>
- static inline void apply(Ring& r, Strategy const& strategy)
- {
- // Correct closure if necessary
- detail::correct_closure::close_or_open_ring<Ring>::apply(r);
- // Check area
- typedef typename area_result<Ring, Strategy>::type area_result_type;
- Predicate<area_result_type> predicate;
- area_result_type const zero = 0;
- if (predicate(ring_area_type::apply(r, strategy), zero))
- {
- std::reverse(boost::begin(r), boost::end(r));
- }
- }
- };
- // Correct a polygon: normalizes all rings, sets outer ring clockwise, sets all
- // inner rings counter clockwise (or vice versa depending on orientation)
- template <typename Polygon>
- struct correct_polygon
- {
- typedef typename ring_type<Polygon>::type ring_type;
-
- template <typename Strategy>
- static inline void apply(Polygon& poly, Strategy const& strategy)
- {
- correct_ring
- <
- ring_type,
- std::less
- >::apply(exterior_ring(poly), strategy);
- typename interior_return_type<Polygon>::type
- rings = interior_rings(poly);
- for (typename detail::interior_iterator<Polygon>::type
- it = boost::begin(rings); it != boost::end(rings); ++it)
- {
- correct_ring
- <
- ring_type,
- std::greater
- >::apply(*it, strategy);
- }
- }
- };
- }} // namespace detail::correct
- #endif // DOXYGEN_NO_DETAIL
- #ifndef DOXYGEN_NO_DISPATCH
- namespace dispatch
- {
- template <typename Geometry, typename Tag = typename tag<Geometry>::type>
- struct correct: not_implemented<Tag>
- {};
- template <typename Point>
- struct correct<Point, point_tag>
- : detail::correct::correct_nop<Point>
- {};
- template <typename LineString>
- struct correct<LineString, linestring_tag>
- : detail::correct::correct_nop<LineString>
- {};
- template <typename Segment>
- struct correct<Segment, segment_tag>
- : detail::correct::correct_nop<Segment>
- {};
- template <typename Box>
- struct correct<Box, box_tag>
- : detail::correct::correct_box<Box>
- {};
- template <typename Ring>
- struct correct<Ring, ring_tag>
- : detail::correct::correct_ring
- <
- Ring,
- std::less
- >
- {};
- template <typename Polygon>
- struct correct<Polygon, polygon_tag>
- : detail::correct::correct_polygon<Polygon>
- {};
- template <typename MultiPoint>
- struct correct<MultiPoint, multi_point_tag>
- : detail::correct::correct_nop<MultiPoint>
- {};
- template <typename MultiLineString>
- struct correct<MultiLineString, multi_linestring_tag>
- : detail::correct::correct_nop<MultiLineString>
- {};
- template <typename Geometry>
- struct correct<Geometry, multi_polygon_tag>
- : detail::multi_modify
- <
- Geometry,
- detail::correct::correct_polygon
- <
- typename boost::range_value<Geometry>::type
- >
- >
- {};
- } // namespace dispatch
- #endif // DOXYGEN_NO_DISPATCH
- namespace resolve_variant {
- template <typename Geometry>
- struct correct
- {
- template <typename Strategy>
- static inline void apply(Geometry& geometry, Strategy const& strategy)
- {
- concepts::check<Geometry const>();
- dispatch::correct<Geometry>::apply(geometry, strategy);
- }
- };
- template <BOOST_VARIANT_ENUM_PARAMS(typename T)>
- struct correct<boost::variant<BOOST_VARIANT_ENUM_PARAMS(T)> >
- {
- template <typename Strategy>
- struct visitor: boost::static_visitor<void>
- {
- Strategy const& m_strategy;
- visitor(Strategy const& strategy): m_strategy(strategy) {}
- template <typename Geometry>
- void operator()(Geometry& geometry) const
- {
- correct<Geometry>::apply(geometry, m_strategy);
- }
- };
- template <typename Strategy>
- static inline void
- apply(boost::variant<BOOST_VARIANT_ENUM_PARAMS(T)>& geometry, Strategy const& strategy)
- {
- boost::apply_visitor(visitor<Strategy>(strategy), geometry);
- }
- };
- } // namespace resolve_variant
- /*!
- \brief Corrects a geometry
- \details Corrects a geometry: all rings which are wrongly oriented with respect
- to their expected orientation are reversed. To all rings which do not have a
- closing point and are typed as they should have one, the first point is
- appended. Also boxes can be corrected.
- \ingroup correct
- \tparam Geometry \tparam_geometry
- \param geometry \param_geometry which will be corrected if necessary
- \qbk{[include reference/algorithms/correct.qbk]}
- */
- template <typename Geometry>
- inline void correct(Geometry& geometry)
- {
- typedef typename point_type<Geometry>::type point_type;
- typedef typename strategy::area::services::default_strategy
- <
- typename cs_tag<point_type>::type
- >::type strategy_type;
- resolve_variant::correct<Geometry>::apply(geometry, strategy_type());
- }
- /*!
- \brief Corrects a geometry
- \details Corrects a geometry: all rings which are wrongly oriented with respect
- to their expected orientation are reversed. To all rings which do not have a
- closing point and are typed as they should have one, the first point is
- appended. Also boxes can be corrected.
- \ingroup correct
- \tparam Geometry \tparam_geometry
- \tparam Strategy \tparam_strategy{Area}
- \param geometry \param_geometry which will be corrected if necessary
- \param strategy \param_strategy{area}
- \qbk{distinguish,with strategy}
- \qbk{[include reference/algorithms/correct.qbk]}
- */
- template <typename Geometry, typename Strategy>
- inline void correct(Geometry& geometry, Strategy const& strategy)
- {
- resolve_variant::correct<Geometry>::apply(geometry, strategy);
- }
- #if defined(_MSC_VER)
- #pragma warning(pop)
- #endif
- }} // namespace boost::geometry
- #endif // BOOST_GEOMETRY_ALGORITHMS_CORRECT_HPP
|