aboutsummaryrefslogtreecommitdiff
path: root/src/projections/denoy.cpp
diff options
context:
space:
mode:
authorKristian Evers <kristianevers@gmail.com>2018-12-26 10:53:01 +0100
committerGitHub <noreply@github.com>2018-12-26 10:53:01 +0100
commit81ec8c0b58d66333fccd3703dab62a11132a0331 (patch)
tree41b4675d169f849516649700ab4f50d4cffa4c02 /src/projections/denoy.cpp
parent355d681ed88019e97742344bd642c2fd97e700a1 (diff)
parent80dad6ef2bed4a83008db06986dc168918d48476 (diff)
downloadPROJ-81ec8c0b58d66333fccd3703dab62a11132a0331.tar.gz
PROJ-81ec8c0b58d66333fccd3703dab62a11132a0331.zip
Merge pull request #1203 from rouault/cpp_conversion
Conversion of most C files to C++
Diffstat (limited to 'src/projections/denoy.cpp')
-rw-r--r--src/projections/denoy.cpp32
1 files changed, 32 insertions, 0 deletions
diff --git a/src/projections/denoy.cpp b/src/projections/denoy.cpp
new file mode 100644
index 00000000..5c337c45
--- /dev/null
+++ b/src/projections/denoy.cpp
@@ -0,0 +1,32 @@
+#define PJ_LIB__
+#include <math.h>
+
+#include "projects.h"
+
+PROJ_HEAD(denoy, "Denoyer Semi-Elliptical") "\n\tPCyl, no inv, Sph";
+
+#define C0 0.95
+#define C1 -0.08333333333333333333
+#define C3 0.00166666666666666666
+#define D1 0.9
+#define D5 0.03
+
+
+static XY s_forward (LP lp, PJ *P) { /* Spheroidal, forward */
+ XY xy = {0.0, 0.0};
+ (void) P;
+ xy.y = lp.phi;
+ xy.x = lp.lam;
+ lp.lam = fabs(lp.lam);
+ xy.x *= cos((C0 + lp.lam * (C1 + lp.lam * lp.lam * C3)) *
+ (lp.phi * (D1 + D5 * lp.phi * lp.phi * lp.phi * lp.phi)));
+ return xy;
+}
+
+
+PJ *PROJECTION(denoy) {
+ P->es = 0.0;
+ P->fwd = s_forward;
+
+ return P;
+}