/
githubmirror
/
omim
Обзор
Документация
Войти
/
githubmirror
/
omim
Код
Запросы
0
Пакеты
0
Релизы
0
Аналитика
Безопасность
master
geometry/distance_on_sphere.cpp
34 строки
985 B
Mikhail Gorbushin
[geometry] Fix area geometry tests
29 окт 2019, 12:00
29 окт 2019, 12:00
16566ae
Код
Авторство
О чём код?
#include "geometry/distance_on_sphere.hpp" #include "geometry/point3d.hpp" #include "base/assert.hpp" #include "base/math.hpp" #include <algorithm> #include <cmath> using namespace std; namespace ms { double DistanceOnSphere(double lat1Deg, double lon1Deg, double lat2Deg, double lon2Deg) { double const lat1 = base::DegToRad(lat1Deg); double const lat2 = base::DegToRad(lat2Deg); double const dlat = sin((lat2 - lat1) * 0.5); double const dlon = sin((base::DegToRad(lon2Deg) - base::DegToRad(lon1Deg)) * 0.5); double const y = dlat * dlat + dlon * dlon * cos(lat1) * cos(lat2); return 2.0 * atan2(sqrt(y), sqrt(max(0.0, 1.0 - y))); } double DistanceOnEarth(double lat1Deg, double lon1Deg, double lat2Deg, double lon2Deg) { return kEarthRadiusMeters * DistanceOnSphere(lat1Deg, lon1Deg, lat2Deg, lon2Deg); } double DistanceOnEarth(LatLon const & ll1, LatLon const & ll2) { return DistanceOnEarth(ll1.m_lat, ll1.m_lon, ll2.m_lat, ll2.m_lon); } } // namespace ms