123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // Copyright (c) 2007-2014 Barend Gehrels, Amsterdam, the Netherlands.
- // This file was modified by Oracle on 2014, 2017, 2018, 2019.
- // Modifications copyright (c) 2014-2019 Oracle and/or its affiliates.
- // Contributed and/or modified by Menelaos Karavelas, on behalf of Oracle
- // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
- // 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_UNION_HPP
- #define BOOST_GEOMETRY_ALGORITHMS_UNION_HPP
- #include <boost/range/metafunctions.hpp>
- #include <boost/geometry/core/is_areal.hpp>
- #include <boost/geometry/core/point_order.hpp>
- #include <boost/geometry/core/reverse_dispatch.hpp>
- #include <boost/geometry/geometries/concepts/check.hpp>
- #include <boost/geometry/algorithms/not_implemented.hpp>
- #include <boost/geometry/algorithms/detail/overlay/overlay.hpp>
- #include <boost/geometry/policies/robustness/get_rescale_policy.hpp>
- #include <boost/geometry/strategies/default_strategy.hpp>
- #include <boost/geometry/util/range.hpp>
- #include <boost/geometry/algorithms/detail/overlay/linear_linear.hpp>
- #include <boost/geometry/algorithms/detail/overlay/pointlike_pointlike.hpp>
- namespace boost { namespace geometry
- {
- #ifndef DOXYGEN_NO_DISPATCH
- namespace dispatch
- {
- template
- <
- typename Geometry1, typename Geometry2, typename GeometryOut,
- typename TagIn1 = typename tag<Geometry1>::type,
- typename TagIn2 = typename tag<Geometry2>::type,
- typename TagOut = typename tag<GeometryOut>::type,
- bool Areal1 = geometry::is_areal<Geometry1>::value,
- bool Areal2 = geometry::is_areal<Geometry2>::value,
- bool ArealOut = geometry::is_areal<GeometryOut>::value,
- bool Reverse1 = detail::overlay::do_reverse<geometry::point_order<Geometry1>::value>::value,
- bool Reverse2 = detail::overlay::do_reverse<geometry::point_order<Geometry2>::value>::value,
- bool ReverseOut = detail::overlay::do_reverse<geometry::point_order<GeometryOut>::value>::value,
- bool Reverse = geometry::reverse_dispatch<Geometry1, Geometry2>::type::value
- >
- struct union_insert: not_implemented<TagIn1, TagIn2, TagOut>
- {};
- // If reversal is needed, perform it first
- template
- <
- typename Geometry1, typename Geometry2, typename GeometryOut,
- typename TagIn1, typename TagIn2, typename TagOut,
- bool Areal1, bool Areal2, bool ArealOut,
- bool Reverse1, bool Reverse2, bool ReverseOut
- >
- struct union_insert
- <
- Geometry1, Geometry2, GeometryOut,
- TagIn1, TagIn2, TagOut,
- Areal1, Areal2, ArealOut,
- Reverse1, Reverse2, ReverseOut,
- true
- >: union_insert<Geometry2, Geometry1, GeometryOut>
- {
- template <typename RobustPolicy, typename OutputIterator, typename Strategy>
- static inline OutputIterator apply(Geometry1 const& g1,
- Geometry2 const& g2,
- RobustPolicy const& robust_policy,
- OutputIterator out,
- Strategy const& strategy)
- {
- return union_insert
- <
- Geometry2, Geometry1, GeometryOut
- >::apply(g2, g1, robust_policy, out, strategy);
- }
- };
- template
- <
- typename Geometry1, typename Geometry2, typename GeometryOut,
- typename TagIn1, typename TagIn2, typename TagOut,
- bool Reverse1, bool Reverse2, bool ReverseOut
- >
- struct union_insert
- <
- Geometry1, Geometry2, GeometryOut,
- TagIn1, TagIn2, TagOut,
- true, true, true,
- Reverse1, Reverse2, ReverseOut,
- false
- > : detail::overlay::overlay
- <Geometry1, Geometry2, Reverse1, Reverse2, ReverseOut, GeometryOut, overlay_union>
- {};
- // dispatch for union of non-areal geometries
- template
- <
- typename Geometry1, typename Geometry2, typename GeometryOut,
- typename TagIn1, typename TagIn2, typename TagOut,
- bool Reverse1, bool Reverse2, bool ReverseOut
- >
- struct union_insert
- <
- Geometry1, Geometry2, GeometryOut,
- TagIn1, TagIn2, TagOut,
- false, false, false,
- Reverse1, Reverse2, ReverseOut,
- false
- > : union_insert
- <
- Geometry1, Geometry2, GeometryOut,
- typename tag_cast<TagIn1, pointlike_tag, linear_tag>::type,
- typename tag_cast<TagIn2, pointlike_tag, linear_tag>::type,
- TagOut,
- false, false, false,
- Reverse1, Reverse2, ReverseOut,
- false
- >
- {};
- // dispatch for union of linear geometries
- template
- <
- typename Linear1, typename Linear2, typename LineStringOut,
- bool Reverse1, bool Reverse2, bool ReverseOut
- >
- struct union_insert
- <
- Linear1, Linear2, LineStringOut,
- linear_tag, linear_tag, linestring_tag,
- false, false, false,
- Reverse1, Reverse2, ReverseOut,
- false
- > : detail::overlay::linear_linear_linestring
- <
- Linear1, Linear2, LineStringOut, overlay_union
- >
- {};
- // dispatch for point-like geometries
- template
- <
- typename PointLike1, typename PointLike2, typename PointOut,
- bool Reverse1, bool Reverse2, bool ReverseOut
- >
- struct union_insert
- <
- PointLike1, PointLike2, PointOut,
- pointlike_tag, pointlike_tag, point_tag,
- false, false, false,
- Reverse1, Reverse2, ReverseOut,
- false
- > : detail::overlay::union_pointlike_pointlike_point
- <
- PointLike1, PointLike2, PointOut
- >
- {};
- } // namespace dispatch
- #endif // DOXYGEN_NO_DISPATCH
- #ifndef DOXYGEN_NO_DETAIL
- namespace detail { namespace union_
- {
- /*!
- \brief_calc2{union}
- \ingroup union
- \details \details_calc2{union_insert, spatial set theoretic union}.
- \details_insert{union}
- \tparam GeometryOut output geometry type, must be specified
- \tparam Geometry1 \tparam_geometry
- \tparam Geometry2 \tparam_geometry
- \tparam OutputIterator output iterator
- \param geometry1 \param_geometry
- \param geometry2 \param_geometry
- \param out \param_out{union}
- \return \return_out
- */
- template
- <
- typename GeometryOut,
- typename Geometry1,
- typename Geometry2,
- typename OutputIterator
- >
- inline OutputIterator union_insert(Geometry1 const& geometry1,
- Geometry2 const& geometry2,
- OutputIterator out)
- {
- concepts::check<Geometry1 const>();
- concepts::check<Geometry2 const>();
- concepts::check<GeometryOut>();
- typename strategy::intersection::services::default_strategy
- <
- typename cs_tag<GeometryOut>::type
- >::type strategy;
- typedef typename geometry::rescale_overlay_policy_type
- <
- Geometry1,
- Geometry2
- >::type rescale_policy_type;
- rescale_policy_type robust_policy
- = geometry::get_rescale_policy<rescale_policy_type>(
- geometry1, geometry2, strategy);
- return dispatch::union_insert
- <
- Geometry1, Geometry2, GeometryOut
- >::apply(geometry1, geometry2, robust_policy, out, strategy);
- }
- }} // namespace detail::union_
- #endif // DOXYGEN_NO_DETAIL
- namespace resolve_strategy {
- struct union_
- {
- template
- <
- typename Geometry1,
- typename Geometry2,
- typename Collection,
- typename Strategy
- >
- static inline void apply(Geometry1 const& geometry1,
- Geometry2 const& geometry2,
- Collection & output_collection,
- Strategy const& strategy)
- {
- typedef typename boost::range_value<Collection>::type geometry_out;
- typedef typename geometry::rescale_overlay_policy_type
- <
- Geometry1,
- Geometry2,
- typename Strategy::cs_tag
- >::type rescale_policy_type;
- rescale_policy_type robust_policy
- = geometry::get_rescale_policy<rescale_policy_type>(
- geometry1, geometry2, strategy);
- dispatch::union_insert
- <
- Geometry1, Geometry2, geometry_out
- >::apply(geometry1, geometry2, robust_policy,
- range::back_inserter(output_collection),
- strategy);
- }
- template
- <
- typename Geometry1,
- typename Geometry2,
- typename Collection
- >
- static inline void apply(Geometry1 const& geometry1,
- Geometry2 const& geometry2,
- Collection & output_collection,
- default_strategy)
- {
- typedef typename strategy::relate::services::default_strategy
- <
- Geometry1,
- Geometry2
- >::type strategy_type;
- apply(geometry1, geometry2, output_collection, strategy_type());
- }
- };
- } // resolve_strategy
- namespace resolve_variant
- {
-
- template <typename Geometry1, typename Geometry2>
- struct union_
- {
- template <typename Collection, typename Strategy>
- static inline void apply(Geometry1 const& geometry1,
- Geometry2 const& geometry2,
- Collection& output_collection,
- Strategy const& strategy)
- {
- concepts::check<Geometry1 const>();
- concepts::check<Geometry2 const>();
- concepts::check<typename boost::range_value<Collection>::type>();
- resolve_strategy::union_::apply(geometry1, geometry2,
- output_collection,
- strategy);
- }
- };
- template <BOOST_VARIANT_ENUM_PARAMS(typename T), typename Geometry2>
- struct union_<variant<BOOST_VARIANT_ENUM_PARAMS(T)>, Geometry2>
- {
- template <typename Collection, typename Strategy>
- struct visitor: static_visitor<>
- {
- Geometry2 const& m_geometry2;
- Collection& m_output_collection;
- Strategy const& m_strategy;
-
- visitor(Geometry2 const& geometry2,
- Collection& output_collection,
- Strategy const& strategy)
- : m_geometry2(geometry2)
- , m_output_collection(output_collection)
- , m_strategy(strategy)
- {}
-
- template <typename Geometry1>
- void operator()(Geometry1 const& geometry1) const
- {
- union_
- <
- Geometry1,
- Geometry2
- >::apply(geometry1, m_geometry2, m_output_collection, m_strategy);
- }
- };
-
- template <typename Collection, typename Strategy>
- static inline void
- apply(variant<BOOST_VARIANT_ENUM_PARAMS(T)> const& geometry1,
- Geometry2 const& geometry2,
- Collection& output_collection,
- Strategy const& strategy)
- {
- boost::apply_visitor(visitor<Collection, Strategy>(geometry2,
- output_collection,
- strategy),
- geometry1);
- }
- };
- template <typename Geometry1, BOOST_VARIANT_ENUM_PARAMS(typename T)>
- struct union_<Geometry1, variant<BOOST_VARIANT_ENUM_PARAMS(T)> >
- {
- template <typename Collection, typename Strategy>
- struct visitor: static_visitor<>
- {
- Geometry1 const& m_geometry1;
- Collection& m_output_collection;
- Strategy const& m_strategy;
-
- visitor(Geometry1 const& geometry1,
- Collection& output_collection,
- Strategy const& strategy)
- : m_geometry1(geometry1)
- , m_output_collection(output_collection)
- , m_strategy(strategy)
- {}
-
- template <typename Geometry2>
- void operator()(Geometry2 const& geometry2) const
- {
- union_
- <
- Geometry1,
- Geometry2
- >::apply(m_geometry1, geometry2, m_output_collection, m_strategy);
- }
- };
-
- template <typename Collection, typename Strategy>
- static inline void
- apply(Geometry1 const& geometry1,
- variant<BOOST_VARIANT_ENUM_PARAMS(T)> const& geometry2,
- Collection& output_collection,
- Strategy const& strategy)
- {
- boost::apply_visitor(visitor<Collection, Strategy>(geometry1,
- output_collection,
- strategy),
- geometry2);
- }
- };
- template <BOOST_VARIANT_ENUM_PARAMS(typename T1), BOOST_VARIANT_ENUM_PARAMS(typename T2)>
- struct union_<variant<BOOST_VARIANT_ENUM_PARAMS(T1)>, variant<BOOST_VARIANT_ENUM_PARAMS(T2)> >
- {
- template <typename Collection, typename Strategy>
- struct visitor: static_visitor<>
- {
- Collection& m_output_collection;
- Strategy const& m_strategy;
-
- visitor(Collection& output_collection, Strategy const& strategy)
- : m_output_collection(output_collection)
- , m_strategy(strategy)
- {}
-
- template <typename Geometry1, typename Geometry2>
- void operator()(Geometry1 const& geometry1,
- Geometry2 const& geometry2) const
- {
- union_
- <
- Geometry1,
- Geometry2
- >::apply(geometry1, geometry2, m_output_collection, m_strategy);
- }
- };
-
- template <typename Collection, typename Strategy>
- static inline void
- apply(variant<BOOST_VARIANT_ENUM_PARAMS(T1)> const& geometry1,
- variant<BOOST_VARIANT_ENUM_PARAMS(T2)> const& geometry2,
- Collection& output_collection,
- Strategy const& strategy)
- {
- boost::apply_visitor(visitor<Collection, Strategy>(output_collection,
- strategy),
- geometry1, geometry2);
- }
- };
-
- } // namespace resolve_variant
- /*!
- \brief Combines two geometries which each other
- \ingroup union
- \details \details_calc2{union, spatial set theoretic union}.
- \tparam Geometry1 \tparam_geometry
- \tparam Geometry2 \tparam_geometry
- \tparam Collection output collection, either a multi-geometry,
- or a std::vector<Geometry> / std::deque<Geometry> etc
- \tparam Strategy \tparam_strategy{Union_}
- \param geometry1 \param_geometry
- \param geometry2 \param_geometry
- \param output_collection the output collection
- \param strategy \param_strategy{union_}
- \note Called union_ because union is a reserved word.
- \qbk{distinguish,with strategy}
- \qbk{[include reference/algorithms/union.qbk]}
- */
- template
- <
- typename Geometry1,
- typename Geometry2,
- typename Collection,
- typename Strategy
- >
- inline void union_(Geometry1 const& geometry1,
- Geometry2 const& geometry2,
- Collection& output_collection,
- Strategy const& strategy)
- {
- resolve_variant::union_
- <
- Geometry1,
- Geometry2
- >::apply(geometry1, geometry2, output_collection, strategy);
- }
- /*!
- \brief Combines two geometries which each other
- \ingroup union
- \details \details_calc2{union, spatial set theoretic union}.
- \tparam Geometry1 \tparam_geometry
- \tparam Geometry2 \tparam_geometry
- \tparam Collection output collection, either a multi-geometry,
- or a std::vector<Geometry> / std::deque<Geometry> etc
- \param geometry1 \param_geometry
- \param geometry2 \param_geometry
- \param output_collection the output collection
- \note Called union_ because union is a reserved word.
- \qbk{[include reference/algorithms/union.qbk]}
- */
- template
- <
- typename Geometry1,
- typename Geometry2,
- typename Collection
- >
- inline void union_(Geometry1 const& geometry1,
- Geometry2 const& geometry2,
- Collection& output_collection)
- {
- resolve_variant::union_
- <
- Geometry1,
- Geometry2
- >::apply(geometry1, geometry2, output_collection, default_strategy());
- }
- }} // namespace boost::geometry
- #endif // BOOST_GEOMETRY_ALGORITHMS_UNION_HPP
|