Return to route.c CVS log

Up to [Research Unix] / researchv10no / cmd / map / export

Annotation of researchv10no/cmd/map/export/route.c, revision 1.1

1.1     ! root        1: #include "map.h"
        !             2: 
        !             3: /* Given two lat-lon pairs, find an orientation for the
        !             4:    -o option of "map" that will place those two points
        !             5:    on the equator of a standard projection, equally spaced
        !             6:    about the prime meridian.   Both an erect and an inverted
        !             7:    orientation are found.
        !             8: 
        !             9:    Option -t prints out a series of
        !            10:    coordinates that follows the (great circle) track
        !            11:    in the original coordinate system,
        !            12:    followed by ".
        !            13:    This data is just right for map -t.
        !            14: */
        !            15: struct place pole;
        !            16: struct coord twist;
        !            17: int track;
        !            18: 
        !            19: extern void doroute(double, double, double, double, double);
        !            20: 
        !            21: void
        !            22: dorot(double a, double b, double *x, double *y, void (*f)(struct place *))
        !            23: {
        !            24:        struct place g;
        !            25:        deg2rad(a,&g.nlat);
        !            26:        deg2rad(b,&g.wlon);
        !            27:        (*f)(&g);
        !            28:        *x = g.nlat.l/RAD;
        !            29:        *y = g.wlon.l/RAD;
        !            30: }
        !            31: 
        !            32: void
        !            33: rotate(double a, double b, double *x, double *y)
        !            34: {
        !            35:        dorot(a,b,x,y,normalize);
        !            36: }
        !            37: 
        !            38: void
        !            39: rinvert(double a, double b, double *x, double *y)
        !            40: {
        !            41:        dorot(a,b,x,y,invert);
        !            42: }
        !            43: 
        !            44: main(int argc, char **argv)
        !            45: {
        !            46: #pragma ref argv
        !            47:        double an,aw,bn,bw;
        !            48:        track = argc>1;
        !            49:        for(;;){
        !            50:                if(scanf("%lf%lf%lf%lf",&an,&aw,&bn,&bw)!=4)
        !            51:                        break;
        !            52:                doroute(-90.,an,aw,bn,bw);
        !            53:                if(!track)
        !            54:                        doroute(90.,an,aw,bn,bw);
        !            55:        }
        !            56:        return 0;
        !            57: }
        !            58: 
        !            59: void
        !            60: doroute(double dir, double an, double aw, double bn, double bw)
        !            61: {
        !            62:        double an1,aw1,bn1,bw1,pn,pw;
        !            63:        double theta;
        !            64:        double cn,cw,cn1,cw1;
        !            65:        int i,n;
        !            66:        orient(an,aw,0.);
        !            67:        rotate(bn,bw,&bn1,&bw1);
        !            68: /*     printf("b %f %f\n",bn1,bw1);*/
        !            69:        orient(an,aw,bw1);
        !            70:        rinvert(0.,dir,&pn,&pw);
        !            71: /*     printf("p %f %f\n",pn,pw);*/
        !            72:        orient(pn,pw,0.);
        !            73:        rotate(an,aw,&an1,&aw1);
        !            74:        rotate(bn,bw,&bn1,&bw1);
        !            75:        theta = (aw1+bw1)/2;
        !            76: /*     printf("a %f %f \n",an1,aw1);*/
        !            77:        orient(pn,pw,theta);
        !            78:        rotate(an,aw,&an1,&aw1);
        !            79:        rotate(bn,bw,&bn1,&bw1);
        !            80:        if(fabs(aw1-bw1)>180)
        !            81:                if(theta<0.) theta+=180;
        !            82:                else theta -= 180;
        !            83:        orient(pn,pw,theta);
        !            84:        rotate(an,aw,&an1,&aw1);
        !            85:        rotate(bn,bw,&bn1,&bw1);
        !            86:        if(!track) {
        !            87:                printf("-o %.4f %.4f %.4f\n",pn,pw,theta);
        !            88:                printf("A %.4f %.4f\n",an1,aw1);
        !            89:                printf("B %.4f %.4f\n",bn1,bw1);
        !            90:        } else {
        !            91:                cn1 = 0;
        !            92:                n = 1 + fabs(bw1-aw1)/.2;
        !            93:                for(i=0;i<=n;i++) {
        !            94:                        cw1 = aw1 + i*(bw1-aw1)/n;
        !            95:                        rinvert(cn1,cw1,&cn,&cw);
        !            96:                        printf("%f %f\n",cn,cw);
        !            97:                }
        !            98:                printf("\"\n");
        !            99:        }
        !           100: }

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.