13e12c5d1SDavid du Colombier #include "map.h" 23e12c5d1SDavid du Colombier 33e12c5d1SDavid du Colombier static struct coord stdpar0, stdpar1; 4*219b2ee8SDavid du Colombier static double k; 5*219b2ee8SDavid du Colombier static double yeq; 63e12c5d1SDavid du Colombier 73e12c5d1SDavid du Colombier static int 8*219b2ee8SDavid du Colombier Xtrapezoidal(struct place *place, double *x, double *y) 93e12c5d1SDavid du Colombier { 103e12c5d1SDavid du Colombier *y = yeq + place->nlat.l; 113e12c5d1SDavid du Colombier *x = *y*k*place->wlon.l; 123e12c5d1SDavid du Colombier return 1; 133e12c5d1SDavid du Colombier } 143e12c5d1SDavid du Colombier 153e12c5d1SDavid du Colombier proj 16*219b2ee8SDavid du Colombier trapezoidal(double par0, double par1) 173e12c5d1SDavid du Colombier { 183e12c5d1SDavid du Colombier if(fabs(fabs(par0)-fabs(par1))<.1) 193e12c5d1SDavid du Colombier return rectangular(par0); 203e12c5d1SDavid du Colombier deg2rad(par0,&stdpar0); 213e12c5d1SDavid du Colombier deg2rad(par1,&stdpar1); 223e12c5d1SDavid du Colombier if(fabs(par1-par0) < .1) 233e12c5d1SDavid du Colombier k = stdpar1.s; 243e12c5d1SDavid du Colombier else 253e12c5d1SDavid du Colombier k = (stdpar1.c-stdpar0.c)/(stdpar0.l-stdpar1.l); 263e12c5d1SDavid du Colombier yeq = -stdpar1.l - stdpar1.c/k; 273e12c5d1SDavid du Colombier return Xtrapezoidal; 283e12c5d1SDavid du Colombier } 29