/usr/include/osmium/geometry/haversine.hpp is in libosmium-dev 0.0~20140910-9a069af-1+b1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 | #ifndef OSMIUM_GEOMETRY_HAVERSINE_HPP
#define OSMIUM_GEOMETRY_HAVERSINE_HPP
/*
Copyright 2012 Jochen Topf <jochen@topf.org> and others (see README).
This file is part of Osmium (https://github.com/joto/osmium).
Osmium is free software: you can redistribute it and/or modify it under the
terms of the GNU Lesser General Public License or (at your option) the GNU
General Public License as published by the Free Software Foundation, either
version 3 of the Licenses, or (at your option) any later version.
Osmium is distributed in the hope that it will be useful, but WITHOUT ANY
WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A
PARTICULAR PURPOSE. See the GNU Lesser General Public License and the GNU
General Public License for more details.
You should have received a copy of the Licenses along with Osmium. If not, see
<http://www.gnu.org/licenses/>.
*/
#include <osmium/osm/position.hpp>
#include <osmium/osm/way_node_list.hpp>
namespace Osmium {
namespace Geometry {
/**
* @brief Functions to calculate arc distance on Earth using the haversine formula.
*
* See http://en.wikipedia.org/wiki/Haversine_formula
*
* Implementation derived from
* http://blog.julien.cayzac.name/2008/10/arc-and-distance-between-two-points-on.html
*/
namespace Haversine {
/// @brief The usual PI/180 constant
static const double DEG_TO_RAD = 0.017453292519943295769236907684886;
/// @brief Earth's quadratic mean radius for WGS84
static const double EARTH_RADIUS_IN_METERS = 6372797.560856;
inline double distance(const double x1, const double y1, const double x2, const double y2) {
const double lon_arc = (x1 - x2) * DEG_TO_RAD;
const double lat_arc = (y1 - y2) * DEG_TO_RAD;
double lonh = sin(lon_arc * 0.5);
lonh *= lonh;
double lath = sin(lat_arc * 0.5);
lath *= lath;
const double tmp = cos(y1 * DEG_TO_RAD) * cos(y2 * DEG_TO_RAD);
return 2.0 * EARTH_RADIUS_IN_METERS * asin(sqrt(lath + tmp*lonh));
}
inline double distance(const Osmium::OSM::Position& from, const Osmium::OSM::Position& to) {
return distance(from.lon(), from.lat(), to.lon(), to.lat());
}
inline double distance(const Osmium::OSM::WayNodeList& wnl) {
double sum_length=0;
if (wnl.size() < 2) {
return 0.0;
}
for (Osmium::OSM::WayNodeList::const_iterator it = wnl.begin(); it != wnl.end()-1; ++it) {
sum_length += Osmium::Geometry::Haversine::distance(it->position(), (it+1)->position());
}
return sum_length;
}
} // namespace Haversine
} // namespace Geometry
} // namespace Osmium
#endif // OSMIUM_GEOMETRY_HAVERSINE_HPP
|