diff options
Diffstat (limited to 'src/projections/tcc.cpp')
| -rw-r--r-- | src/projections/tcc.cpp | 34 |
1 files changed, 34 insertions, 0 deletions
diff --git a/src/projections/tcc.cpp b/src/projections/tcc.cpp new file mode 100644 index 00000000..64fdc182 --- /dev/null +++ b/src/projections/tcc.cpp @@ -0,0 +1,34 @@ +#define PJ_LIB__ + +#include <math.h> + +#include "proj.h" +#include "projects.h" + +PROJ_HEAD(tcc, "Transverse Central Cylindrical") "\n\tCyl, Sph, no inv"; + +#define EPS10 1.e-10 + + +static XY s_forward (LP lp, PJ *P) { /* Spheroidal, forward */ + XY xy = {0.0, 0.0}; + double b, bt; + + b = cos (lp.phi) * sin (lp.lam); + if ((bt = 1. - b * b) < EPS10) { + proj_errno_set(P, PJD_ERR_TOLERANCE_CONDITION); + return xy; + } + xy.x = b / sqrt(bt); + xy.y = atan2 (tan (lp.phi) , cos (lp.lam)); + return xy; +} + + +PJ *PROJECTION(tcc) { + P->es = 0.; + P->fwd = s_forward; + P->inv = nullptr; + + return P; +} |
