~ubuntu-branches/ubuntu/breezy/proj/breezy

« back to all changes in this revision

Viewing changes to src/PJ_lcca.c

Committer: Bazaar Package Importer
Author(s): Peter S Galbraith
Date: 2004-11-06 19:44:53 UTC
mto: This revision was merged to the branch mainline in revision 3.
Revision ID: james.westby@ubuntu.com-20041106194453-axnsmkh1zplal8mz

Tags: upstream-4.4.9

Import upstream version 4.4.9

files added:
depcomp

jniwrap

jniwrap/Makefile.am

jniwrap/Makefile.in

jniwrap/README

jniwrap/build.xml

jniwrap/doxygen.cfg

jniwrap/org

jniwrap/org/Makefile.am

jniwrap/org/Makefile.in

jniwrap/org/proj4

jniwrap/org/proj4/LatLong.java

jniwrap/org/proj4/Makefile.am

jniwrap/org/proj4/Makefile.in

jniwrap/org/proj4/Others.java

jniwrap/org/proj4/Proj4.java

jniwrap/org/proj4/Proj4Factory.java

jniwrap/org/proj4/ProjectionData.java

jniwrap/org/proj4/Projections.java

nad/esri

nad/td_out.dist

nad/testntv2

nad/testvarious

src/PJ_geos.c

src/PJ_krovak.c

src/PJ_lcca.c

src/PJ_sterea.c

src/jniproj.c

src/org_proj4_Projections.h

src/pj_gauss.c

src/pj_geocent.c

src/pj_gridinfo.c

src/pj_gridlist.c

files removed:
src/stamp-h.in

src/strtod.c

files modified:
ChangeLog

Makefile.am

Makefile.in

NEWS

README

aclocal.m4

config.guess

config.sub

configure

configure.in

ltmain.sh

man/Makefile.in

man/man1/Makefile.in

man/man1/cs2cs.1

man/man1/geod.1

man/man1/proj.1

man/man3/Makefile.in

missing

nad/Makefile.am

nad/Makefile.in

nad/README

nad/epsg

src/Makefile.am

src/Makefile.in

src/PJ_aea.c

src/PJ_aeqd.c

src/PJ_airy.c

src/PJ_aitoff.c

src/PJ_goode.c

src/PJ_nzmg.c

src/PJ_ob_tran.c

src/PJ_oea.c

src/adjlon.c

src/bchgen.c

src/biveval.c

src/cs2cs.c

src/dmstor.c

src/emess.c

src/gen_cheb.c

src/geocent.c

src/geocent.h

src/geod.c

src/geod_for.c

src/geod_inv.c

src/geod_set.c

src/geodesic.h

src/makefile.vc

src/nad2nad.c

src/nad_cvt.c

src/nad_init.c

src/pj_apply_gridshift.c

src/pj_datums.c

src/pj_ellps.c

src/pj_errno.c

src/pj_init.c

src/pj_list.c

src/pj_list.h

src/pj_msfn.c

src/pj_open_lib.c

src/pj_param.c

src/pj_pr_list.c

src/pj_qsfn.c

src/pj_release.c

src/pj_strerrno.c

src/pj_transform.c

src/pj_tsfn.c

src/pj_units.c

src/pj_utils.c

src/proj.c

src/proj.def

src/proj_api.h

src/proj_config.h.in

src/projects.h

Show diffs side-by-side

added added

removed removed

src/PJ_lcca.c

static const char RCS_ID[] = "$Id: PJ_lcca.c,v 1.1 2003/03/04 02:59:41 warmerda Exp $";

/* PROJ.4 Cartographic Projection System -- Revision Log:

**$Log: PJ_lcca.c,v $

**Revision 1.1 2003/03/04 02:59:41 warmerda

**New

#define MAX_ITER 10

#define DEL_TOL 1e-12

#define PROJ_PARMS__ \

double *en; \

double r0, l, M0; \

double C;

#define PJ_LIB__

#include <projects.h>

PROJ_HEAD(lcca, "Lambert Conformal Conic Alternative")

"\n\tConic, Sph&Ell\n\tlat_0=";

static double /* func to compute dr */

fS(double S, double C) {

return(S * ( 1. + S * S * C));

}

static double /* deriv of fs */

fSp(double S, double C) {

return(1. + 3.* S * S * C);

}

FORWARD(e_forward); /* ellipsoid */

double S, S3, r, dr;

S = pj_mlfn(lp.phi, sin(lp.phi), cos(lp.phi), P->en) - P->M0;

dr = fS(S, P->C);

r = P->r0 - dr;

xy.x = P->k0 * (r * sin( lp.lam *= P->l ) );

xy.y = P->k0 * (P->r0 - r * cos(lp.lam) );

return (xy);

}

INVERSE(e_inverse); /* ellipsoid & spheroid */

double theta, dr, S, dif;

int i;

xy.x /= P->k0;

xy.y /= P->k0;

theta = atan2(xy.x , P->r0 - xy.y);

dr = xy.y - xy.x * tan(0.5 * theta);

lp.lam = theta / P->l;

S = dr;

for (i = MAX_ITER; i ; --i) {

S -= (dif = (fS(S, P->C) - dr) / fSp(S, P->C));

if (fabs(dif) < DEL_TOL) break;

}

if (!i) I_ERROR

lp.phi = pj_inv_mlfn(S + P->M0, P->es, P->en);

return (lp);

}

FREEUP; if (P) { if (P->en) pj_dalloc(P->en); pj_dalloc(P); } }

ENTRY0(lcca)

double s2p0, N0, R0, tan0, tan20;

if (!(P->en = pj_enfn(P->es))) E_ERROR_0;

if (!pj_param(P->params, "tlat_0").i) E_ERROR(50);

if (P->phi0 == 0.) E_ERROR(51);

P->l = sin(P->phi0);

P->M0 = pj_mlfn(P->phi0, P->l, cos(P->phi0), P->en);

s2p0 = P->l * P->l;

R0 = 1. / (1. - P->es * s2p0);

N0 = sqrt(R0);

R0 *= P->one_es * N0;

tan0 = tan(P->phi0);

tan20 = tan0 * tan0;

P->r0 = N0 / tan0;

P->C = 1. / (6. * R0 * N0);

P->inv = e_inverse;

P->fwd = e_forward;

ENDENTRY(P)

Older »