aboutsummaryrefslogtreecommitdiff
path: root/src/PJ_crast.c
diff options
context:
space:
mode:
authorFrank Warmerdam <warmerdam@pobox.com>1999-03-18 16:34:52 +0000
committerFrank Warmerdam <warmerdam@pobox.com>1999-03-18 16:34:52 +0000
commit565a4bd035b9d4a83955808efef20f1d8dfa24cf (patch)
tree75785fc897708023f1ccdaf40079afcbaaf0fd3a /src/PJ_crast.c
downloadPROJ-565a4bd035b9d4a83955808efef20f1d8dfa24cf.tar.gz
PROJ-565a4bd035b9d4a83955808efef20f1d8dfa24cf.zip
New
git-svn-id: http://svn.osgeo.org/metacrs/proj/trunk@776 4e78687f-474d-0410-85f9-8d5e500ac6b2
Diffstat (limited to 'src/PJ_crast.c')
-rw-r--r--src/PJ_crast.c25
1 files changed, 25 insertions, 0 deletions
diff --git a/src/PJ_crast.c b/src/PJ_crast.c
new file mode 100644
index 00000000..a72203e8
--- /dev/null
+++ b/src/PJ_crast.c
@@ -0,0 +1,25 @@
+#ifndef lint
+static const char SCCSID[]="@(#)PJ_crast.c 4.1 94/02/15 GIE REL";
+#endif
+#define PJ_LIB__
+# include <projects.h>
+PROJ_HEAD(crast, "Craster Parabolic (Putnins P4)")
+"\n\tPCyl., Sph.";
+#define XM 0.97720502380583984317
+#define RXM 1.02332670794648848847
+#define YM 3.06998012383946546542
+#define RYM 0.32573500793527994772
+#define THIRD 0.333333333333333333
+FORWARD(s_forward); /* spheroid */
+ lp.phi *= THIRD;
+ xy.x = XM * lp.lam * (2. * cos(lp.phi + lp.phi) - 1.);
+ xy.y = YM * sin(lp.phi);
+ return (xy);
+}
+INVERSE(s_inverse); /* spheroid */
+ lp.phi = 3. * asin(xy.y * RYM);
+ lp.lam = xy.x * RXM / (2. * cos((lp.phi + lp.phi) * THIRD) - 1);
+ return (lp);
+}
+FREEUP; if (P) pj_dalloc(P); }
+ENTRY0(crast) P->es = 0.; P->inv = s_inverse; P->fwd = s_forward; ENDENTRY(P)