45 lines · 980 B · C
RawDownload
/*  * This file is part of the UCB release of Plan 9. It is subject to the license * terms in the LICENSE file found in the top-level directory of this * distribution and at http://akaros.cs.berkeley.edu/files/Plan9License. No * part of the UCB release of Plan 9, including this file, may be copied, * modified, propagated, or distributed except according to the terms contained * in the LICENSE file. */ #include <u.h>#include <libc.h>#include "map.h" static intXmercator(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);} projmercator(void){	return(Xmercator);} static double ecc = ECC; static intXspmercator(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);} projsp_mercator(void){	return(Xspmercator);}