CassiniSoldner.hpp

Go to the documentation of this file.

 
#if !defined(GEOGRAPHICLIB_CASSINISOLDNER_HPP)
#define GEOGRAPHICLIB_CASSINISOLDNER_HPP 1
 
#include <GeographicLib/Geodesic.hpp>
#include <GeographicLib/GeodesicLine.hpp>
#include <GeographicLib/Constants.hpp>
 
namespace GeographicLib {
 
  class GEOGRAPHICLIB_EXPORT CassiniSoldner {
  private:
    typedef Math::real real;
    Geodesic _earth;
    GeodesicLine _meridian;
    real _sbet0, _cbet0;
    static const unsigned maxit_ = 10;
 
  public:
 
    explicit CassiniSoldner(const Geodesic& earth = Geodesic::WGS84());
 
    CassiniSoldner(real lat0, real lon0,
                   const Geodesic& earth = Geodesic::WGS84());
 
    void Reset(real lat0, real lon0);
 
    void Forward(real lat, real lon,
                 real& x, real& y, real& azi, real& rk) const;
 
    void Reverse(real x, real y,
                 real& lat, real& lon, real& azi, real& rk) const;
 
    void Forward(real lat, real lon,
                 real& x, real& y) const {
      real azi, rk;
      Forward(lat, lon, x, y, azi, rk);
    }
 
    void Reverse(real x, real y,
                 real& lat, real& lon) const {
      real azi, rk;
      Reverse(x, y, lat, lon, azi, rk);
    }
 
 
    bool Init() const { return _meridian.Init(); }
 
    Math::real LatitudeOrigin() const
    { return _meridian.Latitude(); }
 
    Math::real LongitudeOrigin() const
    { return _meridian.Longitude(); }
 
    Math::real MajorRadius() const { return _earth.MajorRadius(); }
 
    Math::real Flattening() const { return _earth.Flattening(); }
 
  };
 
} // namespace GeographicLib
 
#endif  // GEOGRAPHICLIB_CASSINISOLDNER_HPP