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