/src/PROJ/src/projections/urmfps.cpp
Line | Count | Source |
1 | | |
2 | | |
3 | | #include <errno.h> |
4 | | #include <math.h> |
5 | | |
6 | | #include "proj.h" |
7 | | #include "proj_internal.h" |
8 | | |
9 | | PROJ_HEAD(urmfps, "Urmaev Flat-Polar Sinusoidal") "\n\tPCyl, Sph\n\tn="; |
10 | | PROJ_HEAD(wag1, "Wagner I (Kavrayskiy VI)") "\n\tPCyl, Sph"; |
11 | | |
12 | | namespace { // anonymous namespace |
13 | | struct pj_urmfps { |
14 | | double n, C_y; |
15 | | }; |
16 | | } // anonymous namespace |
17 | | |
18 | 336 | #define C_x 0.8773826753 |
19 | 43 | #define Cy 1.139753528477 |
20 | | |
21 | 336 | static PJ_XY urmfps_s_forward(PJ_LP lp, PJ *P) { /* Spheroidal, forward */ |
22 | 336 | PJ_XY xy = {0.0, 0.0}; |
23 | 336 | lp.phi = aasin(P->ctx, |
24 | 336 | static_cast<struct pj_urmfps *>(P->opaque)->n * sin(lp.phi)); |
25 | 336 | xy.x = C_x * lp.lam * cos(lp.phi); |
26 | 336 | xy.y = static_cast<struct pj_urmfps *>(P->opaque)->C_y * lp.phi; |
27 | 336 | return xy; |
28 | 336 | } |
29 | | |
30 | 0 | static PJ_LP urmfps_s_inverse(PJ_XY xy, PJ *P) { /* Spheroidal, inverse */ |
31 | 0 | PJ_LP lp = {0.0, 0.0}; |
32 | 0 | xy.y /= static_cast<struct pj_urmfps *>(P->opaque)->C_y; |
33 | 0 | lp.phi = aasin(P->ctx, |
34 | 0 | sin(xy.y) / static_cast<struct pj_urmfps *>(P->opaque)->n); |
35 | 0 | lp.lam = xy.x / (C_x * cos(xy.y)); |
36 | 0 | return lp; |
37 | 0 | } |
38 | | |
39 | 43 | static PJ *urmfps_setup(PJ *P) { |
40 | 43 | static_cast<struct pj_urmfps *>(P->opaque)->C_y = |
41 | 43 | Cy / static_cast<struct pj_urmfps *>(P->opaque)->n; |
42 | 43 | P->es = 0.; |
43 | 43 | P->inv = urmfps_s_inverse; |
44 | 43 | P->fwd = urmfps_s_forward; |
45 | 43 | return P; |
46 | 43 | } |
47 | | |
48 | 3 | PJ *PJ_PROJECTION(urmfps) { |
49 | 3 | struct pj_urmfps *Q = |
50 | 3 | static_cast<struct pj_urmfps *>(calloc(1, sizeof(struct pj_urmfps))); |
51 | 3 | if (nullptr == Q) |
52 | 0 | return pj_default_destructor(P, PROJ_ERR_OTHER /*ENOMEM*/); |
53 | | |
54 | 3 | P->opaque = Q; |
55 | | |
56 | 3 | if (!pj_param(P->ctx, P->params, "tn").i) { |
57 | 3 | proj_log_error(P, _("Missing parameter n.")); |
58 | 3 | return pj_default_destructor(P, PROJ_ERR_INVALID_OP_MISSING_ARG); |
59 | 3 | } |
60 | | |
61 | 0 | Q->n = pj_param(P->ctx, P->params, "dn").f; |
62 | 0 | if (Q->n <= 0. || Q->n > 1.) { |
63 | 0 | proj_log_error(P, |
64 | 0 | _("Invalid value for n: it should be in ]0,1] range.")); |
65 | 0 | return pj_default_destructor(P, PROJ_ERR_INVALID_OP_ILLEGAL_ARG_VALUE); |
66 | 0 | } |
67 | | |
68 | 0 | return urmfps_setup(P); |
69 | 0 | } |
70 | | |
71 | 43 | PJ *PJ_PROJECTION(wag1) { |
72 | 43 | struct pj_urmfps *Q = |
73 | 43 | static_cast<struct pj_urmfps *>(calloc(1, sizeof(struct pj_urmfps))); |
74 | 43 | if (nullptr == Q) |
75 | 0 | return pj_default_destructor(P, PROJ_ERR_OTHER /*ENOMEM*/); |
76 | 43 | P->opaque = Q; |
77 | | |
78 | 43 | static_cast<struct pj_urmfps *>(P->opaque)->n = |
79 | 43 | 0.8660254037844386467637231707; |
80 | 43 | return urmfps_setup(P); |
81 | 43 | } |
82 | | |
83 | | #undef C_x |
84 | | #undef Cy |