TransverseMercator.hpp

Go to the documentation of this file.

 
#if !defined(GEOGRAPHICLIB_TRANSVERSEMERCATOR_HPP)
#define GEOGRAPHICLIB_TRANSVERSEMERCATOR_HPP 1
 
#include <GeographicLib/Constants.hpp>
 
#if !defined(GEOGRAPHICLIB_TRANSVERSEMERCATOR_ORDER)
 
#  define GEOGRAPHICLIB_TRANSVERSEMERCATOR_ORDER \
  (GEOGRAPHICLIB_PRECISION == 2 ? 6 : \
   (GEOGRAPHICLIB_PRECISION == 1 ? 4 : 8))
#endif
 
namespace GeographicLib {
 
  class GEOGRAPHICLIB_EXPORT TransverseMercator {
  private:
    typedef Math::real real;
    static const int maxpow_ = GEOGRAPHICLIB_TRANSVERSEMERCATOR_ORDER;
    static const int numit_ = 5;
    real _a, _f, _k0, _e2, _es, _e2m,  _c, _n;
    // _alp[0] and _bet[0] unused
    real _a1, _b1, _alp[maxpow_ + 1], _bet[maxpow_ + 1];
    friend class Ellipsoid;           // For access to taupf, tauf.
  public:
 
    TransverseMercator(real a, real f, real k0);
 
    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 TransverseMercator& UTM();
  };
 
} // namespace GeographicLib
 
#endif  // GEOGRAPHICLIB_TRANSVERSEMERCATOR_HPP