From 28994509cc11ac6a5443054dfae1fedfb69039bc Mon Sep 17 00:00:00 2001 From: rsc Date: Wed, 21 Apr 2004 22:19:33 +0000 Subject: Why not? --- src/cmd/map/libmap/mercator.c | 36 ++++++++++++++++++++++++++++++++++++ 1 file changed, 36 insertions(+) create mode 100644 src/cmd/map/libmap/mercator.c (limited to 'src/cmd/map/libmap/mercator.c') diff --git a/src/cmd/map/libmap/mercator.c b/src/cmd/map/libmap/mercator.c new file mode 100644 index 00000000..0ab63653 --- /dev/null +++ b/src/cmd/map/libmap/mercator.c @@ -0,0 +1,36 @@ +#include +#include +#include "map.h" + +static int +Xmercator(struct place *place, double *x, double *y) +{ + if(fabs(place->nlat.l) > 80.*RAD) + return(-1); + *x = -place->wlon.l; + *y = 0.5*log((1+place->nlat.s)/(1-place->nlat.s)); + return(1); +} + +proj +mercator(void) +{ + return(Xmercator); +} + +static double ecc = ECC; + +static int +Xspmercator(struct place *place, double *x, double *y) +{ + if(Xmercator(place,x,y) < 0) + return(-1); + *y += 0.5*ecc*log((1-ecc*place->nlat.s)/(1+ecc*place->nlat.s)); + return(1); +} + +proj +sp_mercator(void) +{ + return(Xspmercator); +} -- cgit v1.2.3