36 lines · 567 B · C
RawDownload
#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);}