summaryrefslogtreecommitdiff log msg author committer range
 ```1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 ``` ``````// Boost.Geometry (aka GGL, Generic Geometry Library) // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands. // This file was modified by Oracle on 2014, 2016, 2017. // Modifications copyright (c) 2014-2017, Oracle and/or its affiliates. // Contributed and/or modified by Vissarion Fysikopoulos, on behalf of Oracle // Contributed and/or modified by Adam Wulkiewicz, 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_DETAIL_AZIMUTH_HPP #define BOOST_GEOMETRY_ALGORITHMS_DETAIL_AZIMUTH_HPP #include #include #include #include #include #include #include #include #include namespace boost { namespace geometry { // An azimuth is an angle between a vector/segment from origin to a point of // interest and a reference vector. Typically north-based azimuth is used. // North direction is used as a reference, angle is measured clockwise // (North - 0deg, East - 90deg). For consistency in 2d cartesian CS // the reference vector is Y axis, angle is measured clockwise. // http://en.wikipedia.org/wiki/Azimuth #ifndef DOXYGEN_NO_DISPATCH namespace detail_dispatch { template struct azimuth : not_implemented {}; template struct azimuth { template static inline ReturnType apply(P1 const& p1, P2 const& p2, Spheroid const& spheroid) { return geometry::formula::vincenty_inverse().apply ( get_as_radian<0>(p1), get_as_radian<1>(p1), get_as_radian<0>(p2), get_as_radian<1>(p2), spheroid ).azimuth; } template static inline ReturnType apply(P1 const& p1, P2 const& p2) { return apply(p1, p2, srs::spheroid()); } }; template struct azimuth { template static inline ReturnType apply(P1 const& p1, P2 const& p2, Sphere const& /*unused*/) { return geometry::formula::spherical_azimuth ( get_as_radian<0>(p1), get_as_radian<1>(p1), get_as_radian<0>(p2), get_as_radian<1>(p2)).azimuth; } template static inline ReturnType apply(P1 const& p1, P2 const& p2) { return apply(p1, p2, 0); // dummy model } }; template struct azimuth : azimuth {}; template struct azimuth { template static inline ReturnType apply(P1 const& p1, P2 const& p2, Plane const& /*unused*/) { ReturnType x = get<0>(p2) - get<0>(p1); ReturnType y = get<1>(p2) - get<1>(p1); // NOTE: azimuth 0 is at Y axis, increasing right // as in spherical/geographic where 0 is at North axis return atan2(x, y); } template static inline ReturnType apply(P1 const& p1, P2 const& p2) { return apply(p1, p2, 0); // dummy model } }; } // detail_dispatch #endif // DOXYGEN_NO_DISPATCH #ifndef DOXYGEN_NO_DETAIL namespace detail { /// Calculate azimuth between two points. /// The result is in radians. template inline ReturnType azimuth(Point1 const& p1, Point2 const& p2) { return detail_dispatch::azimuth < ReturnType, typename geometry::cs_tag::type >::apply(p1, p2); } /// Calculate azimuth between two points. /// The result is in radians. template inline ReturnType azimuth(Point1 const& p1, Point2 const& p2, Model const& model) { return detail_dispatch::azimuth < ReturnType, typename geometry::cs_tag::type >::apply(p1, p2, model); } } // namespace detail #endif // DOXYGEN_NO_DETAIL }} // namespace boost::geometry #endif // BOOST_GEOMETRY_ALGORITHMS_DETAIL_AZIMUTH_HPP ``````