diff options
| author | Frank Warmerdam <warmerdam@pobox.com> | 1999-03-18 16:34:52 +0000 |
|---|---|---|
| committer | Frank Warmerdam <warmerdam@pobox.com> | 1999-03-18 16:34:52 +0000 |
| commit | 565a4bd035b9d4a83955808efef20f1d8dfa24cf (patch) | |
| tree | 75785fc897708023f1ccdaf40079afcbaaf0fd3a /src/PJ_rpoly.c | |
| download | PROJ-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_rpoly.c')
| -rw-r--r-- | src/PJ_rpoly.c | 38 |
1 files changed, 38 insertions, 0 deletions
diff --git a/src/PJ_rpoly.c b/src/PJ_rpoly.c new file mode 100644 index 00000000..7873b8c0 --- /dev/null +++ b/src/PJ_rpoly.c @@ -0,0 +1,38 @@ +#ifndef lint +static const char SCCSID[]="@(#)PJ_rpoly.c 4.1 94/02/15 GIE REL"; +#endif +#define PROJ_PARMS__ \ + double phi1; \ + double fxa; \ + double fxb; \ + int mode; +#define EPS 1e-9 +#define PJ_LIB__ +#include <projects.h> +PROJ_HEAD(rpoly, "Rectangular Polyconic") + "\n\tConic, Sph., no inv.\n\tlat_ts="; +FORWARD(s_forward); /* spheroid */ + double fa; + + if (P->mode) + fa = tan(lp.lam * P->fxb) * P->fxa; + else + fa = 0.5 * lp.lam; + if (fabs(lp.phi) < EPS) { + xy.x = fa + fa; + xy.y = - P->phi0; + } else { + xy.y = 1. / tan(lp.phi); + xy.x = sin(fa = 2. * atan(fa * sin(lp.phi))) * xy.y; + xy.y = lp.phi - P->phi0 + (1. - cos(fa)) * xy.y; + } + return (xy); +} +FREEUP; if (P) pj_dalloc(P); } +ENTRY0(rpoly) + if ((P->mode = (P->phi1 = fabs(pj_param(P->params, "rlat_ts").f)) > EPS)) { + P->fxb = 0.5 * sin(P->phi1); + P->fxa = 0.5 / P->fxb; + } + P->es = 0.; P->fwd = s_forward; +ENDENTRY(P) |