123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // Copyright (c) 2007-2015 Barend Gehrels, Amsterdam, the Netherlands.
- // Copyright (c) 2008-2015 Bruno Lalande, Paris, France.
- // Copyright (c) 2009-2015 Mateusz Loskot, London, UK.
- // This file was modified by Oracle on 2015-2018.
- // Modifications copyright (c) 2015-2018, 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_STRATEGIES_CARTESIAN_POINT_IN_BOX_HPP
- #define BOOST_GEOMETRY_STRATEGIES_CARTESIAN_POINT_IN_BOX_HPP
- #include <boost/geometry/core/access.hpp>
- #include <boost/geometry/core/coordinate_dimension.hpp>
- #include <boost/geometry/core/cs.hpp>
- #include <boost/geometry/strategies/covered_by.hpp>
- #include <boost/geometry/strategies/within.hpp>
- #include <boost/geometry/util/normalize_spheroidal_coordinates.hpp>
- namespace boost { namespace geometry { namespace strategy
- {
- namespace within
- {
- #ifndef DOXYGEN_NO_DETAIL
- namespace detail
- {
- struct within_coord
- {
- template <typename Value1, typename Value2>
- static inline bool apply(Value1 const& value, Value2 const& min_value, Value2 const& max_value)
- {
- return value > min_value && value < max_value;
- }
- };
- struct covered_by_coord
- {
- template <typename Value1, typename Value2>
- static inline bool apply(Value1 const& value, Value2 const& min_value, Value2 const& max_value)
- {
- return value >= min_value && value <= max_value;
- }
- };
- template <typename Geometry, std::size_t Dimension, typename CSTag>
- struct within_range
- : within_coord
- {};
- template <typename Geometry, std::size_t Dimension, typename CSTag>
- struct covered_by_range
- : covered_by_coord
- {};
- // NOTE: the result would be the same if instead of structs defined below
- // the above xxx_range were used with the following arguments:
- // (min_value + diff_min, min_value, max_value)
- struct within_longitude_diff
- {
- template <typename CalcT>
- static inline bool apply(CalcT const& diff_min, CalcT const& min_value, CalcT const& max_value)
- {
- CalcT const c0 = 0;
- return diff_min > c0
- && (min_value + diff_min < max_value
- /*|| max_value - diff_min > min_value*/);
- }
- };
- struct covered_by_longitude_diff
- {
- template <typename CalcT>
- static inline bool apply(CalcT const& diff_min, CalcT const& min_value, CalcT const& max_value)
- {
- return min_value + diff_min <= max_value
- /*|| max_value - diff_min >= min_value*/;
- }
- };
- template <typename Geometry,
- typename CoordCheck,
- typename DiffCheck>
- struct longitude_range
- {
- template <typename Value1, typename Value2>
- static inline bool apply(Value1 const& value, Value2 const& min_value, Value2 const& max_value)
- {
- typedef typename select_most_precise
- <
- Value1, Value2
- >::type calc_t;
- typedef typename geometry::detail::cs_angular_units<Geometry>::type units_t;
- typedef math::detail::constants_on_spheroid<calc_t, units_t> constants;
- if (CoordCheck::apply(value, min_value, max_value))
- {
- return true;
- }
- // min <= max <=> diff >= 0
- calc_t const diff_ing = max_value - min_value;
- // if containing covers the whole globe it contains all
- if (diff_ing >= constants::period())
- {
- return true;
- }
- // calculate positive longitude translation with min_value as origin
- calc_t const diff_min = math::longitude_distance_unsigned<units_t, calc_t>(min_value, value);
- return DiffCheck::template apply<calc_t>(diff_min, min_value, max_value);
- }
- };
- // spherical_equatorial_tag, spherical_polar_tag and geographic_cat are casted to spherical_tag
- template <typename Geometry>
- struct within_range<Geometry, 0, spherical_tag>
- : longitude_range<Geometry, within_coord, within_longitude_diff>
- {};
- template <typename Geometry>
- struct covered_by_range<Geometry, 0, spherical_tag>
- : longitude_range<Geometry, covered_by_coord, covered_by_longitude_diff>
- {};
- template
- <
- template <typename, std::size_t, typename> class SubStrategy,
- typename CSTag, // cartesian_tag or spherical_tag
- std::size_t Dimension,
- std::size_t DimensionCount
- >
- struct relate_point_box_loop
- {
- template <typename Point, typename Box>
- static inline bool apply(Point const& point, Box const& box)
- {
- if (! SubStrategy<Point, Dimension, CSTag>::apply(get<Dimension>(point),
- get<min_corner, Dimension>(box),
- get<max_corner, Dimension>(box))
- )
- {
- return false;
- }
- return relate_point_box_loop
- <
- SubStrategy,
- CSTag,
- Dimension + 1, DimensionCount
- >::apply(point, box);
- }
- };
- template
- <
- template <typename, std::size_t, typename> class SubStrategy,
- typename CSTag,
- std::size_t DimensionCount
- >
- struct relate_point_box_loop<SubStrategy, CSTag, DimensionCount, DimensionCount>
- {
- template <typename Point, typename Box>
- static inline bool apply(Point const& , Box const& )
- {
- return true;
- }
- };
- } // namespace detail
- #endif // DOXYGEN_NO_DETAIL
- struct cartesian_point_box
- {
- template <typename Point, typename Box>
- static inline bool apply(Point const& point, Box const& box)
- {
- return detail::relate_point_box_loop
- <
- detail::within_range,
- cartesian_tag,
- 0, dimension<Point>::value
- >::apply(point, box);
- }
- };
- struct spherical_point_box
- {
- template <typename Point, typename Box>
- static inline bool apply(Point const& point, Box const& box)
- {
- return detail::relate_point_box_loop
- <
- detail::within_range,
- spherical_tag,
- 0, dimension<Point>::value
- >::apply(point, box);
- }
- };
- #ifndef DOXYGEN_NO_STRATEGY_SPECIALIZATIONS
- namespace services
- {
- template <typename Point, typename Box>
- struct default_strategy
- <
- Point, Box,
- point_tag, box_tag,
- pointlike_tag, areal_tag,
- cartesian_tag, cartesian_tag
- >
- {
- typedef within::cartesian_point_box type;
- };
- // spherical_equatorial_tag, spherical_polar_tag and geographic_cat are casted to spherical_tag
- template <typename Point, typename Box>
- struct default_strategy
- <
- Point, Box,
- point_tag, box_tag,
- pointlike_tag, areal_tag,
- spherical_tag, spherical_tag
- >
- {
- typedef within::spherical_point_box type;
- };
- } // namespace services
- #endif // DOXYGEN_NO_STRATEGY_SPECIALIZATIONS
- } // namespace within
- namespace covered_by
- {
- struct cartesian_point_box
- {
- template <typename Point, typename Box>
- static inline bool apply(Point const& point, Box const& box)
- {
- return within::detail::relate_point_box_loop
- <
- within::detail::covered_by_range,
- cartesian_tag,
- 0, dimension<Point>::value
- >::apply(point, box);
- }
- };
- struct spherical_point_box
- {
- template <typename Point, typename Box>
- static inline bool apply(Point const& point, Box const& box)
- {
- return within::detail::relate_point_box_loop
- <
- within::detail::covered_by_range,
- spherical_tag,
- 0, dimension<Point>::value
- >::apply(point, box);
- }
- };
- #ifndef DOXYGEN_NO_STRATEGY_SPECIALIZATIONS
- namespace services
- {
- template <typename Point, typename Box>
- struct default_strategy
- <
- Point, Box,
- point_tag, box_tag,
- pointlike_tag, areal_tag,
- cartesian_tag, cartesian_tag
- >
- {
- typedef covered_by::cartesian_point_box type;
- };
- // spherical_equatorial_tag, spherical_polar_tag and geographic_cat are casted to spherical_tag
- template <typename Point, typename Box>
- struct default_strategy
- <
- Point, Box,
- point_tag, box_tag,
- pointlike_tag, areal_tag,
- spherical_tag, spherical_tag
- >
- {
- typedef covered_by::spherical_point_box type;
- };
- } // namespace services
- #endif // DOXYGEN_NO_STRATEGY_SPECIALIZATIONS
- } // namespace covered_by
- }}} // namespace boost::geometry::strategy
- #endif // BOOST_GEOMETRY_STRATEGIES_CARTESIAN_POINT_IN_BOX_HPP
|