1*59cc4ca5SDavid du Colombier #include <u.h> 2*59cc4ca5SDavid du Colombier #include <libc.h> 33e12c5d1SDavid du Colombier #include "map.h" 43e12c5d1SDavid du Colombier 53e12c5d1SDavid du Colombier int Xpolyconic(struct place * place,double * x,double * y)6219b2ee8SDavid du ColombierXpolyconic(struct place *place, double *x, double *y) 73e12c5d1SDavid du Colombier { 83e12c5d1SDavid du Colombier double r, alpha; 9219b2ee8SDavid du Colombier double lat2, lon2; 103e12c5d1SDavid du Colombier if(fabs(place->nlat.l) > .01) { 113e12c5d1SDavid du Colombier r = place->nlat.c / place->nlat.s; 123e12c5d1SDavid du Colombier alpha = place->wlon.l * place->nlat.s; 133e12c5d1SDavid du Colombier *y = place->nlat.l + r*(1 - cos(alpha)); 143e12c5d1SDavid du Colombier *x = - r*sin(alpha); 153e12c5d1SDavid du Colombier } else { 163e12c5d1SDavid du Colombier lon2 = place->wlon.l * place->wlon.l; 173e12c5d1SDavid du Colombier lat2 = place->nlat.l * place->nlat.l; 183e12c5d1SDavid du Colombier *y = place->nlat.l * (1+(lon2/2)*(1-(8+lon2)*lat2/12)); 193e12c5d1SDavid du Colombier *x = - place->wlon.l * (1-lat2*(3+lon2)/6); 203e12c5d1SDavid du Colombier } 213e12c5d1SDavid du Colombier return(1); 223e12c5d1SDavid du Colombier } 233e12c5d1SDavid du Colombier 243e12c5d1SDavid du Colombier proj polyconic(void)253e12c5d1SDavid du Colombierpolyconic(void) 263e12c5d1SDavid du Colombier { 273e12c5d1SDavid du Colombier return(Xpolyconic); 283e12c5d1SDavid du Colombier } 29