simgear/math/SGGeodesy.hxx

   1 // Copyright (C) 2006  Mathias Froehlich - Mathias.Froehlich@web.de
   2 //
   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.
   7 //
   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.
  12 //
  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.
  16 //
  17
  18 #ifndef SGGeodesy_H
  19 #define SGGeodesy_H
  20
  21 class SGGeodesy {
  22 public:
  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;
  29
  30   /// Takes a cartesian coordinate data and returns the geodetic
  31   /// coordinates.
  32   static void SGCartToGeod(const SGVec3<double>& cart, SGGeod& geod);
  33
  34   /// Takes a geodetic coordinate data and returns the cartesian
  35   /// coordinates.
  36   static void SGGeodToCart(const SGGeod& geod, SGVec3<double>& cart);
  37
  38   /// Takes a geodetic coordinate data and returns the sea level radius.
  39   static double SGGeodToSeaLevelRadius(const SGGeod& geod);
  40
  41   /// Takes a cartesian coordinate data and returns the geocentric
  42   /// coordinates.
  43   static void SGCartToGeoc(const SGVec3<double>& cart, SGGeoc& geoc);
  44
  45   /// Takes a geocentric coordinate data and returns the cartesian
  46   /// coordinates.
  47   static void SGGeocToCart(const SGGeoc& geoc, SGVec3<double>& cart);
  48
  49   // Geodetic course/distance computation
  50   static bool direct(const SGGeod& p1, double course1,
  51                      double distance, SGGeod& p2, double& course2);
  52
  53   static bool inverse(const SGGeod& p1, const SGGeod& p2, double& course1,
  54                       double& course2, double& distance);
  55
  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);
  59
  60   // Geocentric course/distance computation
  61   static void advanceRadM(const SGGeoc& geoc, double course, double distance,
  62                           SGGeoc& result);
  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);
  66
  67   /**
  68    * compute the intersection of two (true) radials (in degrees), or return false
  69    * if no intersection culd be computed.
  70    */
  71   static bool radialIntersection(const SGGeoc& a, double aRadial,
  72     const SGGeoc& b, double bRadial, SGGeoc& result);
  73
  74   static bool radialIntersection(const SGGeod& a, double aRadial,
  75     const SGGeod& b, double bRadial, SGGeod& result);
  76 };
  77
  78 #endif