// Boost.Geometry (aka GGL, Generic Geometry Library) // Copyright (c) 2007-2015 Barend Gehrels, Amsterdam, the Netherlands. // This file was modified by Oracle on 2015. // Modifications copyright (c) 2015 Oracle and/or its affiliates. // Contributed and/or modified by Menelaos Karavelas, 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_SYM_DIFFERENCE_HPP #define BOOST_GEOMETRY_ALGORITHMS_SYM_DIFFERENCE_HPP #include #include #include #include #include #include namespace boost { namespace geometry { #ifndef DOXYGEN_NO_DETAIL namespace detail { namespace sym_difference { template struct compute_difference { template < typename Geometry1, typename Geometry2, typename RobustPolicy, typename OutputIterator, typename Strategy > static inline OutputIterator apply(Geometry1 const& geometry1, Geometry2 const& geometry2, RobustPolicy const& robust_policy, OutputIterator out, Strategy const& strategy) { return geometry::dispatch::intersection_insert < Geometry1, Geometry2, GeometryOut, overlay_difference, geometry::detail::overlay::do_reverse < geometry::point_order::value >::value, geometry::detail::overlay::do_reverse < geometry::point_order::value, true >::value >::apply(geometry1, geometry2, robust_policy, out, strategy); } }; template struct sym_difference_generic { template < typename RobustPolicy, typename OutputIterator, typename Strategy > static inline OutputIterator apply(Geometry1 const& geometry1, Geometry2 const& geometry2, RobustPolicy const& robust_policy, OutputIterator out, Strategy const& strategy) { out = compute_difference < GeometryOut >::apply(geometry1, geometry2, robust_policy, out, strategy); return compute_difference < GeometryOut >::apply(geometry2, geometry1, robust_policy, out, strategy); } }; template struct sym_difference_areal_areal { template < typename RobustPolicy, typename OutputIterator, typename Strategy > static inline OutputIterator apply(Areal1 const& areal1, Areal2 const& areal2, RobustPolicy const& robust_policy, OutputIterator out, Strategy const& strategy) { typedef geometry::model::multi_polygon < GeometryOut > helper_geometry_type; helper_geometry_type diff12, diff21; std::back_insert_iterator oit12(diff12); std::back_insert_iterator oit21(diff21); compute_difference < GeometryOut >::apply(areal1, areal2, robust_policy, oit12, strategy); compute_difference < GeometryOut >::apply(areal2, areal1, robust_policy, oit21, strategy); return geometry::dispatch::union_insert < helper_geometry_type, helper_geometry_type, GeometryOut >::apply(diff12, diff21, robust_policy, out, strategy); } }; }} // namespace detail::sym_difference #endif // DOXYGEN_NO_DETAIL #ifndef DOXYGEN_NO_DISPATCH namespace dispatch { template < typename Geometry1, typename Geometry2, typename GeometryOut, typename TagIn1 = typename geometry::tag_cast < typename tag::type, areal_tag >::type, typename TagIn2 = typename geometry::tag_cast < typename tag::type, areal_tag >::type, typename TagOut = typename geometry::tag::type > struct sym_difference_insert : detail::sym_difference::sym_difference_generic < GeometryOut, Geometry1, Geometry2 > {}; template < typename Areal1, typename Areal2, typename GeometryOut, typename TagOut > struct sym_difference_insert < Areal1, Areal2, GeometryOut, areal_tag, areal_tag, TagOut > : detail::sym_difference::sym_difference_areal_areal < GeometryOut, Areal1, Areal2 > {}; } // namespace dispatch #endif // DOXYGEN_NO_DISPATCH #ifndef DOXYGEN_NO_DETAIL namespace detail { namespace sym_difference { /*! \brief \brief_calc2{symmetric difference} \brief_strategy \ingroup sym_difference \details \details_calc2{symmetric difference, spatial set theoretic symmetric difference (XOR)} \brief_strategy. \details_insert{sym_difference} \tparam GeometryOut output geometry type, must be specified \tparam Geometry1 \tparam_geometry \tparam Geometry2 \tparam_geometry \tparam Strategy \tparam_strategy_overlay \param geometry1 \param_geometry \param geometry2 \param_geometry \param out \param_out{difference} \param strategy \param_strategy{difference} \return \return_out \qbk{distinguish,with strategy} */ template < typename GeometryOut, typename Geometry1, typename Geometry2, typename RobustPolicy, typename OutputIterator, typename Strategy > inline OutputIterator sym_difference_insert(Geometry1 const& geometry1, Geometry2 const& geometry2, RobustPolicy const& robust_policy, OutputIterator out, Strategy const& strategy) { concepts::check(); concepts::check(); concepts::check(); return dispatch::sym_difference_insert < Geometry1, Geometry2, GeometryOut >::apply(geometry1, geometry2, robust_policy, out, strategy); } /*! \brief \brief_calc2{symmetric difference} \ingroup sym_difference \details \details_calc2{symmetric difference, spatial set theoretic symmetric difference (XOR)} \details_insert{sym_difference} \tparam GeometryOut output geometry type, must be specified \tparam Geometry1 \tparam_geometry \tparam Geometry2 \tparam_geometry \param geometry1 \param_geometry \param geometry2 \param_geometry \param out \param_out{difference} \return \return_out */ template < typename GeometryOut, typename Geometry1, typename Geometry2, typename RobustPolicy, typename OutputIterator > inline OutputIterator sym_difference_insert(Geometry1 const& geometry1, Geometry2 const& geometry2, RobustPolicy const& robust_policy, OutputIterator out) { concepts::check(); concepts::check(); concepts::check(); typedef intersection_strategies < typename cs_tag::type, Geometry1, Geometry2, typename geometry::point_type::type, RobustPolicy > strategy_type; return sym_difference_insert(geometry1, geometry2, robust_policy, out, strategy_type()); } }} // namespace detail::sym_difference #endif // DOXYGEN_NO_DETAIL /*! \brief \brief_calc2{symmetric difference} \ingroup sym_difference \details \details_calc2{symmetric difference, spatial set theoretic symmetric difference (XOR)}. \tparam Geometry1 \tparam_geometry \tparam Geometry2 \tparam_geometry \tparam Collection output collection, either a multi-geometry, or a std::vector / std::deque etc \param geometry1 \param_geometry \param geometry2 \param_geometry \param output_collection the output collection \qbk{[include reference/algorithms/sym_difference.qbk]} */ template < typename Geometry1, typename Geometry2, typename Collection > inline void sym_difference(Geometry1 const& geometry1, Geometry2 const& geometry2, Collection& output_collection) { concepts::check(); concepts::check(); typedef typename boost::range_value::type geometry_out; concepts::check(); typedef typename geometry::rescale_overlay_policy_type < Geometry1, Geometry2 >::type rescale_policy_type; rescale_policy_type robust_policy = geometry::get_rescale_policy(geometry1, geometry2); detail::sym_difference::sym_difference_insert( geometry1, geometry2, robust_policy, range::back_inserter(output_collection)); } }} // namespace boost::geometry #endif // BOOST_GEOMETRY_ALGORITHMS_SYM_DIFFERENCE_HPP