123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260 |
- #ifndef BOOST_GEOMETRY_STRATEGIES_AGNOSTIC_POINT_IN_BOX_BY_SIDE_HPP
- #define BOOST_GEOMETRY_STRATEGIES_AGNOSTIC_POINT_IN_BOX_BY_SIDE_HPP
- #include <array>
- #include <boost/core/ignore_unused.hpp>
- #include <boost/geometry/algorithms/assign.hpp>
- #include <boost/geometry/core/access.hpp>
- #include <boost/geometry/core/coordinate_dimension.hpp>
- #include <boost/geometry/strategies/covered_by.hpp>
- #include <boost/geometry/strategies/geographic/side.hpp>
- #include <boost/geometry/strategies/side.hpp>
- #include <boost/geometry/strategies/spherical/ssf.hpp>
- #include <boost/geometry/strategies/within.hpp>
- namespace boost { namespace geometry { namespace strategy
- {
- namespace within
- {
- #ifndef DOXYGEN_NO_DETAIL
- namespace detail
- {
- struct decide_within
- {
- static inline bool apply(int side, bool& result)
- {
- if (side != 1)
- {
- result = false;
- return false;
- }
- return true;
- }
- };
- struct decide_covered_by
- {
- static inline bool apply(int side, bool& result)
- {
- if (side != 1)
- {
- result = side >= 0;
- return false;
- }
- return true;
- }
- };
- template <typename Decide, typename Point, typename Box, typename Strategy>
- inline bool point_in_box_by_side(Point const& point, Box const& box,
- Strategy const& strategy)
- {
- boost::ignore_unused(strategy);
-
-
-
- std::array<typename point_type<Box>::type, 5> bp;
- geometry::detail::assign_box_corners_oriented<true>(box, bp);
- bp[4] = bp[0];
- bool result = true;
- for (int i = 1; i < 5; i++)
- {
- int const side = strategy.apply(point, bp[i - 1], bp[i]);
- if (! Decide::apply(side, result))
- {
- return result;
- }
- }
- return result;
- }
- }
- #endif
- template <typename CalculationType = void>
- struct cartesian_point_box_by_side
- {
- template <typename Point, typename Box>
- static inline bool apply(Point const& point, Box const& box)
- {
- using side_strategy_type
- = typename strategy::side::services::default_strategy
- <cartesian_tag, CalculationType>::type;
- return within::detail::point_in_box_by_side
- <
- within::detail::decide_within
- >(point, box, side_strategy_type());
- }
- };
- template <typename CalculationType = void>
- struct spherical_point_box_by_side
- {
- template <typename Point, typename Box>
- static inline bool apply(Point const& point, Box const& box)
- {
- return within::detail::point_in_box_by_side
- <
- within::detail::decide_within
- >(point, box,
- strategy::side::spherical_side_formula<CalculationType>());
- }
- };
- template
- <
- typename FormulaPolicy = strategy::andoyer,
- typename Spheroid = srs::spheroid<double>,
- typename CalculationType = void
- >
- struct geographic_point_box_by_side
- {
- geographic_point_box_by_side() = default;
- explicit geographic_point_box_by_side(Spheroid const& spheroid)
- : m_side(spheroid)
- {}
- template <typename Point, typename Box>
- bool apply(Point const& point, Box const& box) const
- {
- return within::detail::point_in_box_by_side
- <
- within::detail::decide_within
- >(point, box, m_side);
- }
- Spheroid const& model() const
- {
- return m_side.model();
- }
- private:
- strategy::side::geographic
- <
- FormulaPolicy, Spheroid, CalculationType
- > m_side;
- };
- }
- namespace covered_by
- {
- template <typename CalculationType = void>
- struct cartesian_point_box_by_side
- {
- template <typename Point, typename Box>
- static bool apply(Point const& point, Box const& box)
- {
- using side_strategy_type
- = typename strategy::side::services::default_strategy
- <cartesian_tag, CalculationType>::type;
- return within::detail::point_in_box_by_side
- <
- within::detail::decide_covered_by
- >(point, box, side_strategy_type());
- }
- };
- template <typename CalculationType = void>
- struct spherical_point_box_by_side
- {
- template <typename Point, typename Box>
- static bool apply(Point const& point, Box const& box)
- {
- return within::detail::point_in_box_by_side
- <
- within::detail::decide_covered_by
- >(point, box,
- strategy::side::spherical_side_formula<CalculationType>());
- }
- };
- template
- <
- typename FormulaPolicy = strategy::andoyer,
- typename Spheroid = srs::spheroid<double>,
- typename CalculationType = void
- >
- struct geographic_point_box_by_side
- {
- geographic_point_box_by_side() = default;
- explicit geographic_point_box_by_side(Spheroid const& spheroid)
- : m_side(spheroid)
- {}
- template <typename Point, typename Box>
- bool apply(Point const& point, Box const& box) const
- {
- return within::detail::point_in_box_by_side
- <
- within::detail::decide_covered_by
- >(point, box, m_side);
- }
- Spheroid const& model() const
- {
- return m_side.model();
- }
- private:
- strategy::side::geographic
- <
- FormulaPolicy, Spheroid, CalculationType
- > m_side;
- };
- }
- }}}
- #endif
|