#include #include #include #include "map.h" #include "iplot.h" #define NSYMBOL 20 enum flag { POINT,ENDSEG,ENDSYM }; struct symb { double x, y; char name[10+1]; enum flag flag; } *symbol[NSYMBOL]; static int nsymbol; static double halfrange = 1; extern int halfwidth; extern int vflag; static int getrange(FILE *); static int getsymbol(FILE *, int); static void setrot(struct place *, double, int); static void dorot(struct symb *, double *, double *); void getsyms(char *file) { FILE *sf = fopen(file,"r"); if(sf==0) filerror("cannot open", file); while(nsymbolname)!=1) return 0; i = 0; while(fscanf(sf,"%1s",s)!=EOF) { switch(s[0]) { case 'r': if(!getrange(sf)) break; continue; case 'm': if(i>0) symbol[n][i-1].flag = ENDSEG; continue; case ':': ungetc(s[0],sf); break; default: ungetc(s[0],sf); case 'v': if(fscanf(sf,"%lf %lf",&x,&y)!=2) break; sp[i].x = x*halfwidth/halfrange; sp[i].y = y*halfwidth/halfrange; sp[i].flag = POINT; i++; sp = symbol[n] = (struct symb*)realloc(symbol[n], (i+1)*sizeof(struct symb)); continue; } break; } if(i>0) symbol[n][i-1].flag = ENDSYM; else symbol[n] = 0; return 1; } static int getrange(FILE *sf) { double x,y,xmin,ymin; if(fscanf(sf,"%*s %lf %lf %lf %lf", &xmin,&ymin,&x,&y)!=4) return 0; x -= xmin; y -= ymin; halfrange = (x>y? x: y)/2; if(halfrange<=0) error("bad ra command in -y file"); return 1; } /* r=0 upright;=1 normal;=-1 reverse*/ int putsym(struct place *p, char *name, double s, int r) { int x,y,n; struct symb *sp; double dx,dy; int conn = 0; for(n=0; symbol[n]; n++) if(strcmp(name,symbol[n]->name)==0) break; sp = symbol[n]; if(sp==0) return 0; if(doproj(p,&x,&y)*vflag <= 0) return 1; setrot(p,s,r); for(;;) { dorot(sp,&dx,&dy); conn = cpoint(x+(int)dx,y+(int)dy,conn); switch(sp->flag) { case ENDSEG: conn = 0; case POINT: sp++; continue; case ENDSYM: break; } break; } return 1; } static double rot[2][2]; static void setrot(struct place *p, double s, int r) { double x0,y0,x1,y1; struct place up; up = *p; up.nlat.l += .5*RAD; sincos(&up.nlat); if(r&&(*projection)(p,&x0,&y0)) { if((*projection)(&up,&x1,&y1)<=0) { up.nlat.l -= RAD; sincos(&up.nlat); if((*projection)(&up,&x1,&y1)<=0) goto unit; x1 = x0 - x1; y1 = y0 - y1; } else { x1 -= x0; y1 -= y0; } x1 = r*x1; s /= hypot(x1,y1); rot[0][0] = y1*s; rot[0][1] = x1*s; rot[1][0] = -x1*s; rot[1][1] = y1*s; } else { unit: rot[0][0] = rot[1][1] = s; rot[0][1] = rot[1][0] = 0; } } static void dorot(struct symb *sp, double *px, double *py) { *px = rot[0][0]*sp->x + rot[0][1]*sp->y; *py = rot[1][0]*sp->x + rot[1][1]*sp->y; }