Transform4326ToUTM Class Reference

#include <Utm.h>

Public Member Functions
DEVICE ALWAYS_INLINE	Transform4326ToUTM (unsigned const srid, double const x, double const y)
DEVICE ALWAYS_INLINE double	calculateX () const
DEVICE ALWAYS_INLINE double	calculateY () const

Private Attributes
unsigned const	srid_
double	eta_
double	xi_
bool	small_dlambda_

Detailed Description

Definition at line 120 of file Utm.h.

Constructor & Destructor Documentation

DEVICE ALWAYS_INLINE Transform4326ToUTM::Transform4326ToUTM ( unsigned const  srid, double const  x, double const  y  )

inline

Definition at line 134 of file Utm.h.

References eta_, shared::fastAtanh(), shared::fastCos(), shared::fastSin(), anonymous_namespace{Utm.h}::n, math_consts::rad_div_deg, small_dlambda_, srid_, and xi_.

      : srid_(srid) {
    unsigned const zone = srid_ % 100u;
    double const x0 = zone * 6.0 - 183;
    double const dlambda = (x - x0) * math_consts::rad_div_deg;
    small_dlambda_ = -1.0 / 12 <= dlambda && dlambda <= 1.0 / 12;
    double const phi = y * math_consts::rad_div_deg;
    double const c = 2 * sqrt(n) / (1 + n);  // Boost will have a constexpr sqrt

    // https://en.wikipedia.org/wiki/Universal_Transverse_Mercator_coordinate_system#Simplified_formulae
    double const t = sinh(atanh(sin(phi)) - c * atanh(c * sin(phi)));
    if (small_dlambda_) {
      eta_ = shared::fastAtanh(shared::fastSin(dlambda) / sqrt(1 + t * t));
      xi_ = atan(t / shared::fastCos(dlambda));
    } else {
      eta_ = atanh(sin(dlambda) / sqrt(1 + t * t));
      xi_ = atan(t / cos(dlambda));
    }
    /* Calculate scaling. May be used in the future.
    double sigma_sum = 0;
    for (unsigned j = N; j; --j) {
      sigma_sum += 2 * j * alphas[j] * cos(2 * j * xi_) * cosh(2 * j * eta_);
    }
    sigma_ = 1 + sigma_sum;
    for (unsigned j = N; j; --j) {
      tau_ += 2 * j * alphas[j] * sin(2 * j * xi_) * sinh(2 * j * eta_);
    }
    k_ = k0_A *
         sqrt((1 + sq(tan(phi) * (1 - n) / (1 + n))) * (sq(sigma_) + sq(tau_)) /
                   (sq(t) + sq(cos(dlambda)))) /
         a;
    */
  }

Here is the call graph for this function:

Member Function Documentation

DEVICE ALWAYS_INLINE double Transform4326ToUTM::calculateX ( ) const

inline

Definition at line 170 of file Utm.h.

References anonymous_namespace{Utm.h}::alphas, anonymous_namespace{Utm.h}::E0, eta_, shared::fastSinh(), anonymous_namespace{Utm.h}::k0_A, anonymous_namespace{Utm.h}::N, small_dlambda_, and xi_.

Referenced by conv_4326_utm(), and transform_4326_utm_x().

                                                 {
    double sum = 0;  // Sum in reverse to add smallest numbers first
    if (small_dlambda_) {
      for (unsigned j = N; j; --j) {
        sum += alphas[j] * cos(2 * j * xi_) * shared::fastSinh(2 * j * eta_);
      }
    } else {
      for (unsigned j = N; j; --j) {
        sum += alphas[j] * cos(2 * j * xi_) * sinh(2 * j * eta_);
      }
    }
    return E0 + k0_A * (eta_ + sum);
  }

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double Transform4326ToUTM::calculateY ( ) const

inline

Definition at line 195 of file Utm.h.

References anonymous_namespace{Utm.h}::alphas, eta_, shared::fastCosh(), anonymous_namespace{Utm.h}::k0_A, anonymous_namespace{Utm.h}::N, small_dlambda_, srid_, and xi_.

Referenced by conv_4326_utm(), and transform_4326_utm_y().

                                                 {
    double const N0 = (32700 < srid_) * 10e6;  // UTM False northing (m)
    double sum = 0;
    if (small_dlambda_) {
      for (unsigned j = N; j; --j) {
        sum += alphas[j] * sin(2 * j * xi_) * shared::fastCosh(2 * j * eta_);
      }
    } else {
      for (unsigned j = N; j; --j) {
        sum += alphas[j] * sin(2 * j * xi_) * cosh(2 * j * eta_);
      }
    }
    return N0 + k0_A * (xi_ + sum);
  }

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

double Transform4326ToUTM::eta_

private

Definition at line 123 of file Utm.h.

Referenced by calculateX(), calculateY(), and Transform4326ToUTM().

bool Transform4326ToUTM::small_dlambda_

private

Definition at line 125 of file Utm.h.

Referenced by calculateX(), calculateY(), and Transform4326ToUTM().

unsigned const Transform4326ToUTM::srid_

private

Definition at line 122 of file Utm.h.

Referenced by calculateY(), and Transform4326ToUTM().

double Transform4326ToUTM::xi_

private

Definition at line 124 of file Utm.h.

Referenced by calculateX(), calculateY(), and Transform4326ToUTM().

---

The documentation for this class was generated from the following file:

• /home/jenkins-slave/workspace/core-os-doxygen/Geospatial/Utm.h