3 // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
5 // This file was modified by Oracle on 2014-2020.
6 // Modifications copyright (c) 2014-2020 Oracle and/or its affiliates.
8 // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
10 // Use, modification and distribution is subject to the Boost Software License,
11 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
12 // http://www.boost.org/LICENSE_1_0.txt)
14 #ifndef BOOST_GEOMETRY_STRATEGIES_GEOGRAPHIC_SIDE_HPP
15 #define BOOST_GEOMETRY_STRATEGIES_GEOGRAPHIC_SIDE_HPP
17 #include <boost/geometry/core/cs.hpp>
18 #include <boost/geometry/core/access.hpp>
19 #include <boost/geometry/core/radian_access.hpp>
20 #include <boost/geometry/core/radius.hpp>
22 #include <boost/geometry/formulas/spherical.hpp>
24 #include <boost/geometry/srs/spheroid.hpp>
26 #include <boost/geometry/util/math.hpp>
27 #include <boost/geometry/util/promote_floating_point.hpp>
28 #include <boost/geometry/util/select_calculation_type.hpp>
30 #include <boost/geometry/strategy/geographic/envelope.hpp>
32 #include <boost/geometry/strategies/geographic/disjoint_segment_box.hpp>
33 #include <boost/geometry/strategies/geographic/parameters.hpp>
34 #include <boost/geometry/strategies/side.hpp>
35 #include <boost/geometry/strategies/spherical/point_in_point.hpp>
36 //#include <boost/geometry/strategies/concepts/side_concept.hpp>
39 namespace boost { namespace geometry
43 namespace strategy { namespace side
48 \brief Check at which side of a segment a point lies
49 left of segment (> 0), right of segment (< 0), on segment (0)
51 \tparam FormulaPolicy Geodesic solution formula policy.
52 \tparam Spheroid Reference model of coordinate system.
53 \tparam CalculationType \tparam_calculation
57 [link geometry.reference.srs.srs_spheroid srs::spheroid]
62 typename FormulaPolicy = strategy::andoyer,
63 typename Spheroid = srs::spheroid<double>,
64 typename CalculationType = void
69 typedef geographic_tag cs_tag;
71 typedef strategy::envelope::geographic
76 > envelope_strategy_type;
78 inline envelope_strategy_type get_envelope_strategy() const
80 return envelope_strategy_type(m_model);
83 typedef strategy::disjoint::segment_box_geographic
88 > disjoint_strategy_type;
90 inline disjoint_strategy_type get_disjoint_strategy() const
92 return disjoint_strategy_type(m_model);
95 typedef strategy::within::spherical_point_point equals_point_point_strategy_type;
96 static inline equals_point_point_strategy_type get_equals_point_point_strategy()
98 return equals_point_point_strategy_type();
104 explicit geographic(Spheroid const& model)
108 template <typename P1, typename P2, typename P>
109 inline int apply(P1 const& p1, P2 const& p2, P const& p) const
111 typedef typename promote_floating_point
113 typename select_calculation_type_alt
120 typedef typename FormulaPolicy::template inverse
121 <calc_t, false, true, false, false, false> inverse_formula;
123 calc_t a1p = azimuth<calc_t, inverse_formula>(p1, p, m_model);
124 calc_t a12 = azimuth<calc_t, inverse_formula>(p1, p2, m_model);
126 return formula::azimuth_side_value(a1p, a12);
129 Spheroid const& model() const
135 template <typename ResultType,
136 typename InverseFormulaType,
140 static inline ResultType azimuth(Point1 const& point1, Point2 const& point2,
143 return InverseFormulaType::apply(get_as_radian<0>(point1),
144 get_as_radian<1>(point1),
145 get_as_radian<0>(point2),
146 get_as_radian<1>(point2),
154 }} // namespace strategy::side
157 }} // namespace boost::geometry
160 #endif // BOOST_GEOMETRY_STRATEGIES_GEOGRAPHIC_SIDE_HPP