Chris@16: // Boost.Geometry (aka GGL, Generic Geometry Library)
Chris@16: 
Chris@101: // Copyright (c) 2007-2015 Barend Gehrels, Amsterdam, the Netherlands.
Chris@101: // Copyright (c) 2008-2015 Bruno Lalande, Paris, France.
Chris@101: // Copyright (c) 2009-2015 Mateusz Loskot, London, UK.
Chris@101: // Copyright (c) 2014-2015 Adam Wulkiewicz, Lodz, Poland.
Chris@101: 
Chris@101: // This file was modified by Oracle on 2014, 2015.
Chris@101: // Modifications copyright (c) 2014-2015 Oracle and/or its affiliates.
Chris@101: 
Chris@101: // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
Chris@101: // Contributed and/or modified by Menelaos Karavelas, 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_EQUALS_HPP
Chris@16: #define BOOST_GEOMETRY_ALGORITHMS_EQUALS_HPP
Chris@16: 
Chris@16: 
Chris@16: #include <cstddef>
Chris@16: #include <vector>
Chris@16: 
Chris@16: #include <boost/range.hpp>
Chris@16: 
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/core/access.hpp>
Chris@16: #include <boost/geometry/core/coordinate_dimension.hpp>
Chris@101: #include <boost/geometry/core/geometry_id.hpp>
Chris@16: #include <boost/geometry/core/reverse_dispatch.hpp>
Chris@101: #include <boost/geometry/core/tags.hpp>
Chris@16: 
Chris@16: #include <boost/geometry/geometries/concepts/check.hpp>
Chris@16: 
Chris@101: #include <boost/geometry/algorithms/detail/equals/point_point.hpp>
Chris@16: #include <boost/geometry/algorithms/detail/not.hpp>
Chris@16: #include <boost/geometry/algorithms/not_implemented.hpp>
Chris@16: 
Chris@16: // For trivial checks
Chris@16: #include <boost/geometry/algorithms/area.hpp>
Chris@16: #include <boost/geometry/algorithms/length.hpp>
Chris@16: #include <boost/geometry/util/math.hpp>
Chris@16: #include <boost/geometry/util/select_coordinate_type.hpp>
Chris@16: #include <boost/geometry/util/select_most_precise.hpp>
Chris@16: 
Chris@16: #include <boost/geometry/algorithms/detail/equals/collect_vectors.hpp>
Chris@101: #include <boost/geometry/algorithms/detail/relate/relate.hpp>
Chris@16: 
Chris@101: #include <boost/geometry/views/detail/indexed_point_view.hpp>
Chris@16: 
Chris@16: 
Chris@16: namespace boost { namespace geometry
Chris@16: {
Chris@16: 
Chris@16: #ifndef DOXYGEN_NO_DETAIL
Chris@16: namespace detail { namespace equals
Chris@16: {
Chris@16: 
Chris@16: 
Chris@16: template
Chris@16: <
Chris@16:     std::size_t Dimension,
Chris@16:     std::size_t DimensionCount
Chris@16: >
Chris@16: struct box_box
Chris@16: {
Chris@16:     template <typename Box1, typename Box2>
Chris@16:     static inline bool apply(Box1 const& box1, Box2 const& box2)
Chris@16:     {
Chris@16:         if (!geometry::math::equals(get<min_corner, Dimension>(box1), get<min_corner, Dimension>(box2))
Chris@16:             || !geometry::math::equals(get<max_corner, Dimension>(box1), get<max_corner, Dimension>(box2)))
Chris@16:         {
Chris@16:             return false;
Chris@16:         }
Chris@16:         return box_box<Dimension + 1, DimensionCount>::apply(box1, box2);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: template <std::size_t DimensionCount>
Chris@16: struct box_box<DimensionCount, DimensionCount>
Chris@16: {
Chris@16:     template <typename Box1, typename Box2>
Chris@16:     static inline bool apply(Box1 const& , Box2 const& )
Chris@16:     {
Chris@16:         return true;
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@101: struct segment_segment
Chris@101: {
Chris@101:     template <typename Segment1, typename Segment2>
Chris@101:     static inline bool apply(Segment1 const& segment1, Segment2 const& segment2)
Chris@101:     {
Chris@101:         return equals::equals_point_point(
Chris@101:                     indexed_point_view<Segment1 const, 0>(segment1),
Chris@101:                     indexed_point_view<Segment2 const, 0>(segment2) )
Chris@101:                 ? equals::equals_point_point(
Chris@101:                     indexed_point_view<Segment1 const, 1>(segment1),
Chris@101:                     indexed_point_view<Segment2 const, 1>(segment2) )
Chris@101:                 : ( equals::equals_point_point(
Chris@101:                         indexed_point_view<Segment1 const, 0>(segment1),
Chris@101:                         indexed_point_view<Segment2 const, 1>(segment2) )
Chris@101:                  && equals::equals_point_point(
Chris@101:                         indexed_point_view<Segment1 const, 1>(segment1),
Chris@101:                         indexed_point_view<Segment2 const, 0>(segment2) )
Chris@101:                   );
Chris@101:     }
Chris@101: };
Chris@101: 
Chris@101: 
Chris@16: struct area_check
Chris@16: {
Chris@16:     template <typename Geometry1, typename Geometry2>
Chris@16:     static inline bool apply(Geometry1 const& geometry1, Geometry2 const& geometry2)
Chris@16:     {
Chris@16:         return geometry::math::equals(
Chris@16:                 geometry::area(geometry1),
Chris@16:                 geometry::area(geometry2));
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: struct length_check
Chris@16: {
Chris@16:     template <typename Geometry1, typename Geometry2>
Chris@16:     static inline bool apply(Geometry1 const& geometry1, Geometry2 const& geometry2)
Chris@16:     {
Chris@16:         return geometry::math::equals(
Chris@16:                 geometry::length(geometry1),
Chris@16:                 geometry::length(geometry2));
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: template <typename TrivialCheck>
Chris@16: struct equals_by_collection
Chris@16: {
Chris@16:     template <typename Geometry1, typename Geometry2>
Chris@16:     static inline bool apply(Geometry1 const& geometry1, Geometry2 const& geometry2)
Chris@16:     {
Chris@16:         if (! TrivialCheck::apply(geometry1, geometry2))
Chris@16:         {
Chris@16:             return false;
Chris@16:         }
Chris@16: 
Chris@16:         typedef typename geometry::select_most_precise
Chris@16:             <
Chris@16:                 typename select_coordinate_type
Chris@16:                     <
Chris@16:                         Geometry1, Geometry2
Chris@16:                     >::type,
Chris@16:                 double
Chris@16:             >::type calculation_type;
Chris@16: 
Chris@16:         typedef std::vector<collected_vector<calculation_type> > v;
Chris@16:         v c1, c2;
Chris@16: 
Chris@16:         geometry::collect_vectors(c1, geometry1);
Chris@16:         geometry::collect_vectors(c2, geometry2);
Chris@16: 
Chris@16:         if (boost::size(c1) != boost::size(c2))
Chris@16:         {
Chris@16:             return false;
Chris@16:         }
Chris@16: 
Chris@16:         std::sort(c1.begin(), c1.end());
Chris@16:         std::sort(c2.begin(), c2.end());
Chris@16: 
Chris@16:         // Just check if these vectors are equal.
Chris@16:         return std::equal(c1.begin(), c1.end(), c2.begin());
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@101: template<typename Geometry1, typename Geometry2>
Chris@101: struct equals_by_relate
Chris@101:     : detail::relate::relate_base
Chris@101:         <
Chris@101:             detail::relate::static_mask_equals_type,
Chris@101:             Geometry1,
Chris@101:             Geometry2
Chris@101:         >
Chris@101: {};
Chris@16: 
Chris@16: }} // namespace detail::equals
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: template
Chris@16: <
Chris@16:     typename Geometry1,
Chris@16:     typename Geometry2,
Chris@16:     typename Tag1 = typename tag<Geometry1>::type,
Chris@16:     typename Tag2 = typename tag<Geometry2>::type,
Chris@16:     std::size_t DimensionCount = dimension<Geometry1>::type::value,
Chris@16:     bool Reverse = reverse_dispatch<Geometry1, Geometry2>::type::value
Chris@16: >
Chris@16: struct equals: not_implemented<Tag1, Tag2>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: // If reversal is needed, perform it
Chris@16: template
Chris@16: <
Chris@16:     typename Geometry1, typename Geometry2,
Chris@16:     typename Tag1, typename Tag2,
Chris@16:     std::size_t DimensionCount
Chris@16: >
Chris@16: struct equals<Geometry1, Geometry2, Tag1, Tag2, DimensionCount, true>
Chris@16:     : equals<Geometry2, Geometry1, Tag2, Tag1, DimensionCount, false>
Chris@16: {
Chris@16:     static inline bool apply(Geometry1 const& g1, Geometry2 const& g2)
Chris@16:     {
Chris@16:         return equals
Chris@16:             <
Chris@16:                 Geometry2, Geometry1,
Chris@16:                 Tag2, Tag1,
Chris@16:                 DimensionCount,
Chris@16:                 false
Chris@16:             >::apply(g2, g1);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: 
Chris@16: template <typename P1, typename P2, std::size_t DimensionCount, bool Reverse>
Chris@16: struct equals<P1, P2, point_tag, point_tag, DimensionCount, Reverse>
Chris@16:     : geometry::detail::not_
Chris@16:         <
Chris@16:             detail::disjoint::point_point<P1, P2, 0, DimensionCount>
Chris@16:         >
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Box1, typename Box2, std::size_t DimensionCount, bool Reverse>
Chris@16: struct equals<Box1, Box2, box_tag, box_tag, DimensionCount, Reverse>
Chris@16:     : detail::equals::box_box<0, DimensionCount>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Ring1, typename Ring2, bool Reverse>
Chris@16: struct equals<Ring1, Ring2, ring_tag, ring_tag, 2, Reverse>
Chris@16:     : detail::equals::equals_by_collection<detail::equals::area_check>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Polygon1, typename Polygon2, bool Reverse>
Chris@16: struct equals<Polygon1, Polygon2, polygon_tag, polygon_tag, 2, Reverse>
Chris@16:     : detail::equals::equals_by_collection<detail::equals::area_check>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Polygon, typename Ring, bool Reverse>
Chris@16: struct equals<Polygon, Ring, polygon_tag, ring_tag, 2, Reverse>
Chris@16:     : detail::equals::equals_by_collection<detail::equals::area_check>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Ring, typename Box, bool Reverse>
Chris@16: struct equals<Ring, Box, ring_tag, box_tag, 2, Reverse>
Chris@16:     : detail::equals::equals_by_collection<detail::equals::area_check>
Chris@16: {};
Chris@16: 
Chris@16: 
Chris@16: template <typename Polygon, typename Box, bool Reverse>
Chris@16: struct equals<Polygon, Box, polygon_tag, box_tag, 2, Reverse>
Chris@16:     : detail::equals::equals_by_collection<detail::equals::area_check>
Chris@16: {};
Chris@16: 
Chris@101: template <typename Segment1, typename Segment2, std::size_t DimensionCount, bool Reverse>
Chris@101: struct equals<Segment1, Segment2, segment_tag, segment_tag, DimensionCount, Reverse>
Chris@101:     : detail::equals::segment_segment
Chris@101: {};
Chris@101: 
Chris@101: template <typename LineString1, typename LineString2, bool Reverse>
Chris@101: struct equals<LineString1, LineString2, linestring_tag, linestring_tag, 2, Reverse>
Chris@101:     //: detail::equals::equals_by_collection<detail::equals::length_check>
Chris@101:     : detail::equals::equals_by_relate<LineString1, LineString2>
Chris@101: {};
Chris@101: 
Chris@101: template <typename LineString, typename MultiLineString, bool Reverse>
Chris@101: struct equals<LineString, MultiLineString, linestring_tag, multi_linestring_tag, 2, Reverse>
Chris@101:     : detail::equals::equals_by_relate<LineString, MultiLineString>
Chris@101: {};
Chris@101: 
Chris@101: template <typename MultiLineString1, typename MultiLineString2, bool Reverse>
Chris@101: struct equals<MultiLineString1, MultiLineString2, multi_linestring_tag, multi_linestring_tag, 2, Reverse>
Chris@101:     : detail::equals::equals_by_relate<MultiLineString1, MultiLineString2>
Chris@101: {};
Chris@101: 
Chris@101: 
Chris@101: template <typename MultiPolygon1, typename MultiPolygon2, bool Reverse>
Chris@101: struct equals
Chris@101:     <
Chris@101:         MultiPolygon1, MultiPolygon2,
Chris@101:         multi_polygon_tag, multi_polygon_tag,
Chris@101:         2,
Chris@101:         Reverse
Chris@101:     >
Chris@101:     : detail::equals::equals_by_collection<detail::equals::area_check>
Chris@101: {};
Chris@101: 
Chris@101: 
Chris@101: template <typename Polygon, typename MultiPolygon, bool Reverse>
Chris@101: struct equals
Chris@101:     <
Chris@101:         Polygon, MultiPolygon,
Chris@101:         polygon_tag, multi_polygon_tag,
Chris@101:         2,
Chris@101:         Reverse
Chris@101:     >
Chris@101:     : detail::equals::equals_by_collection<detail::equals::area_check>
Chris@101: {};
Chris@101: 
Chris@101: 
Chris@101: } // namespace dispatch
Chris@101: #endif // DOXYGEN_NO_DISPATCH
Chris@101: 
Chris@101: 
Chris@101: namespace resolve_variant {
Chris@16: 
Chris@16: template <typename Geometry1, typename Geometry2>
Chris@101: struct equals
Chris@16: {
Chris@16:     static inline bool apply(Geometry1 const& geometry1,
Chris@16:                              Geometry2 const& geometry2)
Chris@16:     {
Chris@16:         concept::check_concepts_and_equal_dimensions
Chris@16:         <
Chris@16:             Geometry1 const,
Chris@16:             Geometry2 const
Chris@16:         >();
Chris@101: 
Chris@101:         return dispatch::equals<Geometry1, Geometry2>
Chris@101:                        ::apply(geometry1, geometry2);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: template <BOOST_VARIANT_ENUM_PARAMS(typename T), typename Geometry2>
Chris@101: struct equals<boost::variant<BOOST_VARIANT_ENUM_PARAMS(T)>, Geometry2>
Chris@16: {
Chris@16:     struct visitor: static_visitor<bool>
Chris@16:     {
Chris@16:         Geometry2 const& m_geometry2;
Chris@16: 
Chris@16:         visitor(Geometry2 const& geometry2)
Chris@16:             : m_geometry2(geometry2)
Chris@16:         {}
Chris@16: 
Chris@16:         template <typename Geometry1>
Chris@16:         inline bool operator()(Geometry1 const& geometry1) const
Chris@16:         {
Chris@101:             return equals<Geometry1, Geometry2>
Chris@101:                    ::apply(geometry1, m_geometry2);
Chris@16:         }
Chris@16: 
Chris@16:     };
Chris@16: 
Chris@16:     static inline bool apply(
Chris@16:         boost::variant<BOOST_VARIANT_ENUM_PARAMS(T)> const& geometry1,
Chris@16:         Geometry2 const& geometry2
Chris@16:     )
Chris@16:     {
Chris@16:         return apply_visitor(visitor(geometry2), geometry1);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: template <typename Geometry1, BOOST_VARIANT_ENUM_PARAMS(typename T)>
Chris@101: struct equals<Geometry1, boost::variant<BOOST_VARIANT_ENUM_PARAMS(T)> >
Chris@16: {
Chris@16:     struct visitor: static_visitor<bool>
Chris@16:     {
Chris@16:         Geometry1 const& m_geometry1;
Chris@16: 
Chris@16:         visitor(Geometry1 const& geometry1)
Chris@16:             : m_geometry1(geometry1)
Chris@16:         {}
Chris@16: 
Chris@16:         template <typename Geometry2>
Chris@16:         inline bool operator()(Geometry2 const& geometry2) const
Chris@16:         {
Chris@101:             return equals<Geometry1, Geometry2>
Chris@101:                    ::apply(m_geometry1, geometry2);
Chris@16:         }
Chris@16: 
Chris@16:     };
Chris@16: 
Chris@16:     static inline bool apply(
Chris@16:         Geometry1 const& geometry1,
Chris@16:         boost::variant<BOOST_VARIANT_ENUM_PARAMS(T)> const& geometry2
Chris@16:     )
Chris@16:     {
Chris@16:         return apply_visitor(visitor(geometry1), geometry2);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@16: template <
Chris@16:     BOOST_VARIANT_ENUM_PARAMS(typename T1),
Chris@16:     BOOST_VARIANT_ENUM_PARAMS(typename T2)
Chris@16: >
Chris@101: struct equals<
Chris@16:     boost::variant<BOOST_VARIANT_ENUM_PARAMS(T1)>,
Chris@16:     boost::variant<BOOST_VARIANT_ENUM_PARAMS(T2)>
Chris@16: >
Chris@16: {
Chris@16:     struct visitor: static_visitor<bool>
Chris@16:     {
Chris@16:         template <typename Geometry1, typename Geometry2>
Chris@16:         inline bool operator()(Geometry1 const& geometry1,
Chris@16:                                Geometry2 const& geometry2) const
Chris@16:         {
Chris@101:             return equals<Geometry1, Geometry2>
Chris@101:                    ::apply(geometry1, geometry2);
Chris@16:         }
Chris@16: 
Chris@16:     };
Chris@16: 
Chris@16:     static inline bool apply(
Chris@16:         boost::variant<BOOST_VARIANT_ENUM_PARAMS(T1)> const& geometry1,
Chris@16:         boost::variant<BOOST_VARIANT_ENUM_PARAMS(T2)> const& geometry2
Chris@16:     )
Chris@16:     {
Chris@16:         return apply_visitor(visitor(), geometry1, geometry2);
Chris@16:     }
Chris@16: };
Chris@16: 
Chris@101: } // namespace resolve_variant
Chris@16: 
Chris@16: 
Chris@16: /*!
Chris@16: \brief \brief_check{are spatially equal}
Chris@101: \details \details_check12{equals, is spatially equal}. Spatially equal means
Chris@16:     that the same point set is included. A box can therefore be spatially equal
Chris@101:     to a ring or a polygon, or a linestring can be spatially equal to a
Chris@16:     multi-linestring or a segment. This only works theoretically, not all
Chris@16:     combinations are implemented yet.
Chris@16: \ingroup equals
Chris@16: \tparam Geometry1 \tparam_geometry
Chris@16: \tparam Geometry2 \tparam_geometry
Chris@16: \param geometry1 \param_geometry
Chris@16: \param geometry2 \param_geometry
Chris@16: \return \return_check2{are spatially equal}
Chris@16: 
Chris@16: \qbk{[include reference/algorithms/equals.qbk]}
Chris@16: 
Chris@16:  */
Chris@16: template <typename Geometry1, typename Geometry2>
Chris@16: inline bool equals(Geometry1 const& geometry1, Geometry2 const& geometry2)
Chris@16: {
Chris@101:     return resolve_variant::equals<Geometry1, Geometry2>
Chris@101:                           ::apply(geometry1, geometry2);
Chris@16: }
Chris@16: 
Chris@16: 
Chris@16: }} // namespace boost::geometry
Chris@16: 
Chris@16: 
Chris@16: #endif // BOOST_GEOMETRY_ALGORITHMS_EQUALS_HPP
Chris@16: