#ifndef SGGeodesy_H
#define SGGeodesy_H
-class SGGeoc;
-class SGGeod;
-
-template<typename T>
-class SGVec3;
-
class SGGeodesy {
public:
// Hard numbers from the WGS84 standard.
/// Takes a geocentric coordinate data and returns the cartesian
/// coordinates.
static void SGGeocToCart(const SGGeoc& geoc, SGVec3<double>& cart);
+
+ // Geodetic course/distance computation
+ static bool direct(const SGGeod& p1, double course1,
+ double distance, SGGeod& p2, double& course2);
+
+ static bool inverse(const SGGeod& p1, const SGGeod& p2, double& course1,
+ double& course2, double& distance);
+
+ static double courseDeg(const SGGeod& from, const SGGeod& to);
+ static double distanceM(const SGGeod& from, const SGGeod& to);
+ static double distanceNm(const SGGeod& from, const SGGeod& to);
+
+ // Geocentric course/distance computation
+ static void advanceRadM(const SGGeoc& geoc, double course, double distance,
+ SGGeoc& result);
+ static double courseRad(const SGGeoc& from, const SGGeoc& to);
+ static double distanceRad(const SGGeoc& from, const SGGeoc& to);
+ static double distanceM(const SGGeoc& from, const SGGeoc& to);
+
+ /**
+ * compute the intersection of two (true) radials (in degrees), or return false
+ * if no intersection culd be computed.
+ */
+ static bool radialIntersection(const SGGeoc& a, double aRadial,
+ const SGGeoc& b, double bRadial, SGGeoc& result);
+
+ static bool radialIntersection(const SGGeod& a, double aRadial,
+ const SGGeod& b, double bRadial, SGGeod& result);
};
#endif