1 // Copyright (C) 2006 Mathias Froehlich - Mathias.Froehlich@web.de
3 // This library is free software; you can redistribute it and/or
4 // modify it under the terms of the GNU Library General Public
5 // License as published by the Free Software Foundation; either
6 // version 2 of the License, or (at your option) any later version.
8 // This library is distributed in the hope that it will be useful,
9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
11 // Library General Public License for more details.
13 // You should have received a copy of the GNU General Public License
14 // along with this program; if not, write to the Free Software
15 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
23 // Hard numbers from the WGS84 standard.
24 static const double EQURAD;
25 static const double iFLATTENING;
26 static const double SQUASH;
27 static const double STRETCH;
28 static const double POLRAD;
30 /// Takes a cartesian coordinate data and returns the geodetic
32 static void SGCartToGeod(const SGVec3<double>& cart, SGGeod& geod);
34 /// Takes a geodetic coordinate data and returns the cartesian
36 static void SGGeodToCart(const SGGeod& geod, SGVec3<double>& cart);
38 /// Takes a geodetic coordinate data and returns the sea level radius.
39 static double SGGeodToSeaLevelRadius(const SGGeod& geod);
41 /// Takes a cartesian coordinate data and returns the geocentric
43 static void SGCartToGeoc(const SGVec3<double>& cart, SGGeoc& geoc);
45 /// Takes a geocentric coordinate data and returns the cartesian
47 static void SGGeocToCart(const SGGeoc& geoc, SGVec3<double>& cart);
49 // Geodetic course/distance computation
50 static bool direct(const SGGeod& p1, double course1,
51 double distance, SGGeod& p2, double& course2);
53 static bool inverse(const SGGeod& p1, const SGGeod& p2, double& course1,
54 double& course2, double& distance);
56 static double courseDeg(const SGGeod& from, const SGGeod& to);
57 static double distanceM(const SGGeod& from, const SGGeod& to);
58 static double distanceNm(const SGGeod& from, const SGGeod& to);
60 // Geocentric course/distance computation
61 static void advanceRadM(const SGGeoc& geoc, double course, double distance,
63 static double courseRad(const SGGeoc& from, const SGGeoc& to);
64 static double distanceRad(const SGGeoc& from, const SGGeoc& to);
65 static double distanceM(const SGGeoc& from, const SGGeoc& to);
68 * compute the intersection of two (true) radials (in degrees), or return false
69 * if no intersection culd be computed.
71 static bool radialIntersection(const SGGeoc& a, double aRadial,
72 const SGGeoc& b, double bRadial, SGGeoc& result);
74 static bool radialIntersection(const SGGeod& a, double aRadial,
75 const SGGeod& b, double bRadial, SGGeod& result);