TransverseMercatorExact.hpp

Go to the documentation of this file.

 
#if !defined(GEOGRAPHICLIB_TRANSVERSEMERCATOREXACT_HPP)
#define GEOGRAPHICLIB_TRANSVERSEMERCATOREXACT_HPP 1
 
#include <GeographicLib/Constants.hpp>
#include <GeographicLib/EllipticFunction.hpp>
 
namespace GeographicLib {
 
  class GEOGRAPHICLIB_EXPORT TransverseMercatorExact {
  private:
    typedef Math::real real;
    static const int numit_ = 10;
    real tol_, tol2_, taytol_;
    real _a, _f, _k0, _mu, _mv, _e;
    bool _extendp;
    EllipticFunction _Eu, _Ev;
 
    void zeta(real u, real snu, real cnu, real dnu,
              real v, real snv, real cnv, real dnv,
              real& taup, real& lam) const;
 
    void dwdzeta(real u, real snu, real cnu, real dnu,
                 real v, real snv, real cnv, real dnv,
                 real& du, real& dv) const;
 
    bool zetainv0(real psi, real lam, real& u, real& v) const;
    void zetainv(real taup, real lam, real& u, real& v) const;
 
    void sigma(real u, real snu, real cnu, real dnu,
               real v, real snv, real cnv, real dnv,
               real& xi, real& eta) const;
 
    void dwdsigma(real u, real snu, real cnu, real dnu,
                  real v, real snv, real cnv, real dnv,
                  real& du, real& dv) const;
 
    bool sigmainv0(real xi, real eta, real& u, real& v) const;
    void sigmainv(real xi, real eta, real& u, real& v) const;
 
    void Scale(real tau, real lam,
               real snu, real cnu, real dnu,
               real snv, real cnv, real dnv,
               real& gamma, real& k) const;
 
  public:
 
    TransverseMercatorExact(real a, real f, real k0, bool extendp = false);
 
    void Forward(real lon0, real lat, real lon,
                 real& x, real& y, real& gamma, real& k) const;
 
    void Reverse(real lon0, real x, real y,
                 real& lat, real& lon, real& gamma, real& k) const;
 
    void Forward(real lon0, real lat, real lon,
                 real& x, real& y) const {
      real gamma, k;
      Forward(lon0, lat, lon, x, y, gamma, k);
    }
 
    void Reverse(real lon0, real x, real y,
                 real& lat, real& lon) const {
      real gamma, k;
      Reverse(lon0, x, y, lat, lon, gamma, k);
    }
 
 
    Math::real MajorRadius() const { return _a; }
 
    Math::real Flattening() const { return _f; }
 
    Math::real CentralScale() const { return _k0; }
 
    static const TransverseMercatorExact& UTM();
  };
 
} // namespace GeographicLib
 
#endif  // GEOGRAPHICLIB_TRANSVERSEMERCATOREXACT_HPP