1*59cc4ca5SDavid du Colombier #include <u.h> 2*59cc4ca5SDavid du Colombier #include <libc.h> 33e12c5d1SDavid du Colombier #include "map.h" 43e12c5d1SDavid du Colombier 5219b2ee8SDavid du Colombier static double r0, a; 63e12c5d1SDavid du Colombier 73e12c5d1SDavid du Colombier static int Xsimpleconic(struct place * place,double * x,double * y)8219b2ee8SDavid du ColombierXsimpleconic(struct place *place, double *x, double *y) 93e12c5d1SDavid du Colombier { 10219b2ee8SDavid du Colombier double r = r0 - place->nlat.l; 11219b2ee8SDavid du Colombier double t = a*place->wlon.l; 123e12c5d1SDavid du Colombier *x = -r*sin(t); 133e12c5d1SDavid du Colombier *y = -r*cos(t); 143e12c5d1SDavid du Colombier return 1; 153e12c5d1SDavid du Colombier } 163e12c5d1SDavid du Colombier 173e12c5d1SDavid du Colombier proj simpleconic(double par0,double par1)18219b2ee8SDavid du Colombiersimpleconic(double par0, double par1) 193e12c5d1SDavid du Colombier { 203e12c5d1SDavid du Colombier struct coord lat0; 213e12c5d1SDavid du Colombier struct coord lat1; 223e12c5d1SDavid du Colombier deg2rad(par0,&lat0); 233e12c5d1SDavid du Colombier deg2rad(par1,&lat1); 243e12c5d1SDavid du Colombier if(fabs(lat0.l+lat1.l)<.01) 253e12c5d1SDavid du Colombier return rectangular(par0); 263e12c5d1SDavid du Colombier if(fabs(lat0.l-lat1.l)<.01) { 273e12c5d1SDavid du Colombier a = lat0.s/lat0.l; 283e12c5d1SDavid du Colombier r0 = lat0.c/lat0.s + lat0.l; 293e12c5d1SDavid du Colombier } else { 303e12c5d1SDavid du Colombier a = (lat1.c-lat0.c)/(lat0.l-lat1.l); 313e12c5d1SDavid du Colombier r0 = ((lat0.c+lat1.c)/a + lat1.l + lat0.l)/2; 323e12c5d1SDavid du Colombier } 333e12c5d1SDavid du Colombier return Xsimpleconic; 343e12c5d1SDavid du Colombier } 35