/src/proj/src/projections/denoy.cpp
Line | Count | Source (jump to first uncovered line) |
1 | | |
2 | | #include <math.h> |
3 | | |
4 | | #include "proj.h" |
5 | | #include "proj_internal.h" |
6 | | |
7 | | PROJ_HEAD(denoy, "Denoyer Semi-Elliptical") "\n\tPCyl, no inv, Sph"; |
8 | | |
9 | 0 | #define C0 0.95 |
10 | 0 | #define C1 -0.08333333333333333333 |
11 | 0 | #define C3 0.00166666666666666666 |
12 | 0 | #define D1 0.9 |
13 | 0 | #define D5 0.03 |
14 | | |
15 | 0 | static PJ_XY denoy_s_forward(PJ_LP lp, PJ *P) { /* Spheroidal, forward */ |
16 | 0 | PJ_XY xy = {0.0, 0.0}; |
17 | 0 | (void)P; |
18 | 0 | xy.y = lp.phi; |
19 | 0 | xy.x = lp.lam; |
20 | 0 | lp.lam = fabs(lp.lam); |
21 | 0 | xy.x *= cos((C0 + lp.lam * (C1 + lp.lam * lp.lam * C3)) * |
22 | 0 | (lp.phi * (D1 + D5 * lp.phi * lp.phi * lp.phi * lp.phi))); |
23 | 0 | return xy; |
24 | 0 | } |
25 | | |
26 | 0 | PJ *PJ_PROJECTION(denoy) { |
27 | 0 | P->es = 0.0; |
28 | 0 | P->fwd = denoy_s_forward; |
29 | |
|
30 | 0 | return P; |
31 | 0 | } |
32 | | |
33 | | #undef C0 |
34 | | #undef C1 |
35 | | #undef C3 |
36 | | #undef D1 |
37 | | #undef D5 |