Port the used bits from CheapRuler, drop dependencies

Author: DennisOSRM

CMakeLists.txt

@@ -432,11 +432,6 @@ include_directories(SYSTEM ${RAPIDJSON_INCLUDE_DIR})
 set(MICROTAR_INCLUDE_DIR "${CMAKE_CURRENT_SOURCE_DIR}/third_party/microtar/src")
 include_directories(SYSTEM ${MICROTAR_INCLUDE_DIR})
 
-set(MBXGEOM_INCLUDE_DIR "${CMAKE_CURRENT_SOURCE_DIR}/third_party/geometry.hpp-0.9.2/include")
-include_directories(SYSTEM ${MBXGEOM_INCLUDE_DIR})
-set(CHEAPRULER_INCLUDE_DIR "${CMAKE_CURRENT_SOURCE_DIR}/third_party/cheap-ruler-cpp-2778eb8/include")
-include_directories(SYSTEM ${CHEAPRULER_INCLUDE_DIR})
-
 add_library(MICROTAR OBJECT "${CMAKE_CURRENT_SOURCE_DIR}/third_party/microtar/src/microtar.c")
 set_property(TARGET MICROTAR PROPERTY POSITION_INDEPENDENT_CODE ON)
 target_no_warning(MICROTAR unused-variable)

include/util/cheap_ruler.hpp

@@ -0,0 +1,88 @@
+#pragma once
+
+#include <cassert>
+#include <cmath>
+#include <cstdint>
+#include <limits>
+#include <tuple>
+#include <utility>
+
+namespace mapbox {
+
+namespace geometry {
+    template <typename T>
+    struct point
+    {
+        using coordinate_type = T;
+
+        constexpr point()
+            : x(), y()
+        {}
+        constexpr point(T x_, T y_)
+            : x(x_), y(y_)
+        {}
+
+        T x;
+        T y;
+    };
+}
+
+namespace cheap_ruler {
+
+using point             = geometry::point<double>;
+
+class CheapRuler {
+
+    // Values that define WGS84 ellipsoid model of the Earth
+    static constexpr double RE = 6378.137; // equatorial radius
+    static constexpr double FE = 1.0 / 298.257223563; // flattening
+
+    static constexpr double E2 = FE * (2 - FE);
+    static constexpr double RAD = M_PI / 180.0;
+
+public:
+    explicit CheapRuler(double latitude) {
+        // Curvature formulas from https://en.wikipedia.org/wiki/Earth_radius#Meridional
+        double mul = RAD * RE * 1000;
+        double coslat = std::cos(latitude * RAD);
+        double w2 = 1 / (1 - E2 * (1 - coslat * coslat));
+        double w = std::sqrt(w2);
+
+        // multipliers for converting longitude and latitude degrees into distance
+        kx = mul * w * coslat;        // based on normal radius of curvature
+        ky = mul * w * w2 * (1 - E2); // based on meridonal radius of curvature
+    }
+
+    double squareDistance(point a, point b) const {
+        auto dx = longDiff(a.x, b.x) * kx;
+        auto dy = (a.y - b.y) * ky;
+        return dx * dx + dy * dy;
+    }
+
+    //
+    // Given two points of the form [x = longitude, y = latitude], returns the distance.
+    //
+    double distance(point a, point b) const {
+        return std::sqrt(squareDistance(a, b));
+    }
+
+    //
+    // Returns the bearing between two points in angles.
+    //
+    double bearing(point a, point b) const {
+        auto dx = longDiff(b.x, a.x) * kx;
+        auto dy = (b.y - a.y) * ky;
+
+        return std::atan2(dx, dy) / RAD;
+    }
+
+private:
+    double ky;
+    double kx;
+    static double longDiff(double a, double b) {
+        return std::remainder(a - b, 360);
+    }
+};
+
+} // namespace cheap_ruler
+} // namespace mapbox

src/util/coordinate_calculation.cpp

@@ -1,12 +1,11 @@
 #include "util/coordinate_calculation.hpp"
+#include "util/cheap_ruler.hpp"
 #include "util/coordinate.hpp"
 #include "util/trigonometry_table.hpp"
 #include "util/web_mercator.hpp"
 
 #include <boost/assert.hpp>
 
-#include <mapbox/cheap_ruler.hpp>
-
 #include <algorithm>
 #include <iterator>
 #include <limits>
@@ -32,7 +31,7 @@ class CheapRulerContainer
         for (int n = 0; n < number_of_rulers; n++)
         {
             cheap_ruler_cache[n] = mapbox::cheap_ruler::CheapRuler(
-                step * (n + 0.5) / COORDINATE_PRECISION, mapbox::cheap_ruler::CheapRuler::Meters);
+                step * (n + 0.5) / COORDINATE_PRECISION);
         }
     };