13e12c5d1SDavid du Colombier #include "map.h" 23e12c5d1SDavid du Colombier 33e12c5d1SDavid du Colombier int 4*219b2ee8SDavid du Colombier Xpolyconic(struct place *place, double *x, double *y) 53e12c5d1SDavid du Colombier { 63e12c5d1SDavid du Colombier double r, alpha; 7*219b2ee8SDavid du Colombier double lat2, lon2; 83e12c5d1SDavid du Colombier if(fabs(place->nlat.l) > .01) { 93e12c5d1SDavid du Colombier r = place->nlat.c / place->nlat.s; 103e12c5d1SDavid du Colombier alpha = place->wlon.l * place->nlat.s; 113e12c5d1SDavid du Colombier *y = place->nlat.l + r*(1 - cos(alpha)); 123e12c5d1SDavid du Colombier *x = - r*sin(alpha); 133e12c5d1SDavid du Colombier } else { 143e12c5d1SDavid du Colombier lon2 = place->wlon.l * place->wlon.l; 153e12c5d1SDavid du Colombier lat2 = place->nlat.l * place->nlat.l; 163e12c5d1SDavid du Colombier *y = place->nlat.l * (1+(lon2/2)*(1-(8+lon2)*lat2/12)); 173e12c5d1SDavid du Colombier *x = - place->wlon.l * (1-lat2*(3+lon2)/6); 183e12c5d1SDavid du Colombier } 193e12c5d1SDavid du Colombier return(1); 203e12c5d1SDavid du Colombier } 213e12c5d1SDavid du Colombier 223e12c5d1SDavid du Colombier proj 233e12c5d1SDavid du Colombier polyconic(void) 243e12c5d1SDavid du Colombier { 253e12c5d1SDavid du Colombier return(Xpolyconic); 263e12c5d1SDavid du Colombier } 27