~ubuntu-branches/ubuntu/vivid/proj/vivid

« back to all changes in this revision

Viewing changes to src/PJ_tmerc.c

Committer: Bazaar Package Importer
Author(s): Peter S Galbraith
Date: 2002-01-11 10:27:12 UTC
Revision ID: james.westby@ubuntu.com-20020111102712-ayi18r8y2eesv0y9

Tags: upstream-4.4.5

Import upstream version 4.4.5

files added:

AUTHORS

COPYING

ChangeLog

INSTALL

Makefile.am

Makefile.in

NEWS

README

aclocal.m4

config.guess

config.sub

configure

configure.in

install-sh

ltconfig

ltmain.sh

man/Makefile.am

man/Makefile.in

man/man1

man/man1/Makefile.am

man/man1/Makefile.in

man/man1/cs2cs.1

man/man1/geod.1

man/man1/nad2nad.1

man/man1/proj.1

man/man3

man/man3/Makefile.am

man/man3/Makefile.in

man/man3/pj_init.3

missing

mkinstalldirs

nad/GL27

nad/Makefile.am

nad/Makefile.in

nad/README

nad/epsg

nad/nad.lst

nad/nad27

nad/nad83

nad/pj_out27.dist

nad/pj_out83.dist

nad/proj_def.dat

nad/test27

nad/test83

nad/world

src/Makefile.am

src/Makefile.in

src/PJ_aea.c

src/PJ_aeqd.c

src/PJ_airy.c

src/PJ_aitoff.c

src/PJ_august.c

src/PJ_bacon.c

src/PJ_bipc.c

src/PJ_boggs.c

src/PJ_bonne.c

src/PJ_cass.c

src/PJ_cc.c

src/PJ_cea.c

src/PJ_chamb.c

src/PJ_collg.c

src/PJ_crast.c

src/PJ_denoy.c

src/PJ_eck1.c

src/PJ_eck2.c

src/PJ_eck3.c

src/PJ_eck4.c

src/PJ_eck5.c

src/PJ_eqc.c

src/PJ_eqdc.c

src/PJ_fahey.c

src/PJ_fouc_s.c

src/PJ_gall.c

src/PJ_gins8.c

src/PJ_gn_sinu.c

src/PJ_gnom.c

src/PJ_goode.c

src/PJ_hammer.c

src/PJ_hatano.c

src/PJ_imw_p.c

src/PJ_labrd.c

src/PJ_laea.c

src/PJ_lagrng.c

src/PJ_larr.c

src/PJ_lask.c

src/PJ_lcc.c

src/PJ_loxim.c

src/PJ_lsat.c

src/PJ_mbt_fps.c

src/PJ_mbtfpp.c

src/PJ_mbtfpq.c

src/PJ_merc.c

src/PJ_mill.c

src/PJ_mod_ster.c

src/PJ_moll.c

src/PJ_mpoly.c

src/PJ_nell.c

src/PJ_nell_h.c

src/PJ_nocol.c

src/PJ_nsper.c

src/PJ_nzmg.c

src/PJ_ob_tran.c

src/PJ_ocea.c

src/PJ_oea.c

src/PJ_omerc.c

src/PJ_ortho.c

src/PJ_poly.c

src/PJ_putp2.c

src/PJ_putp3.c

src/PJ_putp4p.c

src/PJ_putp5.c

src/PJ_putp6.c

src/PJ_robin.c

src/PJ_rpoly.c

src/PJ_sconics.c

src/PJ_somerc.c

src/PJ_stere.c

src/PJ_sts.c

src/PJ_tcc.c

src/PJ_tcea.c

src/PJ_tmerc.c

src/PJ_tpeqd.c

src/PJ_urm5.c

src/PJ_urmfps.c

src/PJ_vandg.c

src/PJ_vandg2.c

src/PJ_vandg4.c

src/PJ_wag2.c

src/PJ_wag3.c

src/PJ_wag7.c

src/PJ_wink1.c

src/PJ_wink2.c

src/aasincos.c

src/adjlon.c

src/bch2bps.c

src/bchgen.c

src/biveval.c

src/cs2cs.c

src/dmstor.c

src/emess.c

src/emess.h

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/mk_cheby.c

src/nad2bin.c

src/nad2nad.c

src/nad_cvt.c

src/nad_init.c

src/nad_intr.c

src/nad_list.h

src/p_series.c

src/pj_apply_gridshift.c

src/pj_auth.c

src/pj_datum_set.c

src/pj_datums.c

src/pj_deriv.c

src/pj_ell_set.c

src/pj_ellps.c

src/pj_errno.c

src/pj_factors.c

src/pj_fwd.c

src/pj_init.c

src/pj_inv.c

src/pj_latlong.c

src/pj_list.c

src/pj_list.h

src/pj_malloc.c

src/pj_mlfn.c

src/pj_msfn.c

src/pj_open_lib.c

src/pj_param.c

src/pj_phi2.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/pj_zpoly1.c

src/proj.c

src/proj.def

src/proj_api.h

src/proj_config.h.in

src/projects.h

src/rtodms.c

src/stamp-h.in

src/strtod.c

src/vector1.c

Show diffs side-by-side

added added

removed removed

src/PJ_tmerc.c

#ifndef lint

static const char SCCSID[]="@(#)PJ_tmerc.c 4.2 94/06/02 GIE REL";

#endif

#define PROJ_PARMS__ \

double esp; \

double ml0; \

double *en;

#define PJ_LIB__

#include <projects.h>

PROJ_HEAD(tmerc, "Transverse Mercator") "\n\tCyl, Sph&Ell";

PROJ_HEAD(utm, "Universal Transverse Mercator (UTM)")

"\n\tCyl, Sph\n\tzone= south";

#define EPS10 1.e-10

#define aks0 P->esp

#define aks5 P->ml0

#define FC1 1.

#define FC2 .5

#define FC3 .16666666666666666666

#define FC4 .08333333333333333333

#define FC5 .05

#define FC6 .03333333333333333333

#define FC7 .02380952380952380952

#define FC8 .01785714285714285714

FORWARD(e_forward); /* ellipse */

double al, als, n, cosphi, sinphi, t;

sinphi = sin(lp.phi); cosphi = cos(lp.phi);

t = fabs(cosphi) > 1e-10 ? sinphi/cosphi : 0.;

t *= t;

al = cosphi * lp.lam;

als = al * al;

al /= sqrt(1. - P->es * sinphi * sinphi);

n = P->esp * cosphi * cosphi;

xy.x = P->k0 * al * (FC1 +

FC3 * als * (1. - t + n +

FC5 * als * (5. + t * (t - 18.) + n * (14. - 58. * t)

+ FC7 * als * (61. + t * ( t * (179. - t) - 479. ) )

)));

xy.y = P->k0 * (pj_mlfn(lp.phi, sinphi, cosphi, P->en) - P->ml0 +

sinphi * al * lp.lam * FC2 * ( 1. +

FC4 * als * (5. - t + n * (9. + 4. * n) +

FC6 * als * (61. + t * (t - 58.) + n * (270. - 330 * t)

+ FC8 * als * (1385. + t * ( t * (543. - t) - 3111.) )

))));

return (xy);

}

FORWARD(s_forward); /* sphere */

double b, cosphi;

b = (cosphi = cos(lp.phi)) * sin(lp.lam);

if (fabs(fabs(b) - 1.) <= EPS10) F_ERROR;

xy.x = aks5 * log((1. + b) / (1. - b));

if ((b = fabs( xy.y = cosphi * cos(lp.lam) / sqrt(1. - b * b) )) >= 1.) {

if ((b - 1.) > EPS10) F_ERROR

else xy.y = 0.;

} else

xy.y = acos(xy.y);

if (lp.phi < 0.) xy.y = -xy.y;

xy.y = aks0 * (xy.y - P->phi0);

return (xy);

}

INVERSE(e_inverse); /* ellipsoid */

double n, con, cosphi, d, ds, sinphi, t;

lp.phi = pj_inv_mlfn(P->ml0 + xy.y / P->k0, P->es, P->en);

if (fabs(lp.phi) >= HALFPI) {

lp.phi = xy.y < 0. ? -HALFPI : HALFPI;

lp.lam = 0.;

} else {

sinphi = sin(lp.phi);

cosphi = cos(lp.phi);

t = fabs(cosphi) > 1e-10 ? sinphi/cosphi : 0.;

n = P->esp * cosphi * cosphi;

d = xy.x * sqrt(con = 1. - P->es * sinphi * sinphi) / P->k0;

con *= t;

t *= t;

ds = d * d;

lp.phi -= (con * ds / (1.-P->es)) * FC2 * (1. -

ds * FC4 * (5. + t * (3. - 9. * n) + n * (1. - 4 * n) -

ds * FC6 * (61. + t * (90. - 252. * n +

45. * t) + 46. * n

- ds * FC8 * (1385. + t * (3633. + t * (4095. + 1574. * t)) )

)));

lp.lam = d*(FC1 -

ds*FC3*( 1. + 2.*t + n -

ds*FC5*(5. + t*(28. + 24.*t + 8.*n) + 6.*n

- ds * FC7 * (61. + t * (662. + t * (1320. + 720. * t)) )

))) / cosphi;

}

return (lp);

}

INVERSE(s_inverse); /* sphere */

double h, g;

h = exp(xy.x / aks0);

g = .5 * (h - 1. / h);

h = cos(P->phi0 + xy.y / aks0);

lp.phi = asin(sqrt((1. - h * h) / (1. + g * g)));

if (xy.y < 0.) lp.phi = -lp.phi;

100

lp.lam = (g || h) ? atan2(g, h) : 0.;

101

return (lp);

102

}

103

FREEUP;

104

if (P) {

105

if (P->en)

106

pj_dalloc(P->en);

107

pj_dalloc(P);

108

}

109

}

110

static PJ *

111

setup(PJ *P) { /* general initialization */

112

if (P->es) {

113

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

114

E_ERROR_0;

115

P->ml0 = pj_mlfn(P->phi0, sin(P->phi0), cos(P->phi0), P->en);

116

P->esp = P->es / (1. - P->es);

117

P->inv = e_inverse;

118

P->fwd = e_forward;

119

} else {

120

aks0 = P->k0;

121

aks5 = .5 * aks0;

122

P->inv = s_inverse;

123

P->fwd = s_forward;

124

}

125

return P;

126

}

127

ENTRY1(tmerc, en)

128

ENDENTRY(setup(P))

129

ENTRY1(utm, en)

130

int zone;

131

132

if (!P->es) E_ERROR(-34);

133

P->y0 = pj_param(P->params, "bsouth").i ? 10000000. : 0.;

134

P->x0 = 500000.;

135

if (pj_param(P->params, "tzone").i) /* zone input ? */

136

if ((zone = pj_param(P->params, "izone").i) > 0 && zone <= 60)

137

--zone;

138

else

139

E_ERROR(-35)

140

else /* nearest central meridian input */

141

if ((zone = floor((adjlon(P->lam0) + PI) * 30. / PI)) < 0)

142

zone = 0;

143

else if (zone >= 60)

144

zone = 59;

145

P->lam0 = (zone + .5) * PI / 30. - PI;

146

P->k0 = 0.9996;

147

P->phi0 = 0.;

148

ENDENTRY(setup(P))

Older »