plot_topo.c

Go to the documentation of this file.

 
#include "prototypes.h"
 
 
void plot_topo(Pixmap pixm,int h, int display_flag)
 
{
 
extern float topo_scale;
/*
extern XImage *image;
*/
extern Display *display;
extern int npoint[2],zoom;
extern GC gc;
extern struct topo *topo;
extern unsigned long tmap[50];
extern Widget drawing_area;
extern struct display_set display_set[4];
int i,j,k;
int color,MY=10;
float xrat,yrat,mult,dmult,cmult;
long nx[4],ny[4];
Dimension width,height;
Arg args[10];
unsigned long XSIZE,YSIZE;
XPoint point[5];
long xmin,xmax,ymin,ymax,lint1,lint2;
signed long xcen,ycen;
int iflag,maxi,maxj,mini,minj,nlong,dx,dy,di[4],dj[4];
int tdi,ceni,cenj;
int xx,yy;
 
  dx=1;dy=1;
 
  di[0]=0;dj[0]=0;
  di[1]=1;dj[1]=0; 
  di[2]=1;dj[2]=1;
  di[3]=0;dj[3]=1;
 
dmult=display_set[display_flag].dmult[h];
xmin= display_set[display_flag].xmin[h];
xmax= display_set[display_flag].xmax[h];
ymin= display_set[display_flag].ymin[h];
ymax= display_set[display_flag].ymax[h];
xcen= display_set[display_flag].xcen[h];
ycen= display_set[display_flag].ycen[h];
 
XtSetArg(args[0],XtNwidth,&width);
XtSetArg(args[1],XtNheight,&height);
XtGetValues(drawing_area,args,2);
 
XSIZE=(long)width*10L;
YSIZE=(long)height*10L;
 
xrat=(float)XSIZE/(float)12800;
yrat=(float)YSIZE/(float)9600;
 
if(xrat < yrat)
             mult=xrat;
 
else
             mult=yrat;
 
cmult=dmult*mult;
lint1=(long)npoint[0]*(long)MY;
lint2=(long)npoint[1]*(long)MY;
 
XSetFillRule(display,gc,WindingRule);
 
mini=minj=0;
maxi=topo->maxi;
maxj=topo->maxj;
ceni=maxi/2;
cenj=maxj/2;
 
if(zoom > 1) {
 
        iflag=0;
 
        for(i=0;i<maxi-dx;i++) {
 
            nlong=(xmin+((float)topo->coord[i][cenj].x*cmult-lint1))*zoom + xcen;
           
 
            if(iflag==0 && nlong < xmin)
                       continue;
 
            if(iflag==1 && nlong < xmax)
                       continue;
 
            if(iflag==0)
                   mini=i;
            
            if(iflag==1)
                   maxi=i;
 
            iflag++;
 
 
                          }
 
            iflag=0;
 
         for(j=0;j<maxj-dy;j++) {
            
 
             nlong=(ymin+((float)topo->coord[ceni][j].y*cmult-lint2))*zoom + ycen;
            
 
             if(iflag==0 && nlong < 0)
                        continue;
 
             if(iflag==1 && nlong < YSIZE)
                        continue;
 
             if(iflag==0)
                     minj=j;
 
 
             if(iflag==1)
                     maxj=j;
 
             iflag++;
 
       }
 
      }
 
 
for(i=mini;i<maxi-dx;i=i+dx) {
for(j=minj;j<maxj-dy;j=j+dy) {
 
                 if(topo->value[i][j]==0)
                                     continue;
 
                 if(topo->value[i][j] < 80*topo_scale)
                                     color=3;
 
                 else if(topo->value[i][j] < 160*topo_scale)
                                        color=10;
 
                 else if(topo->value[i][j] < 240*topo_scale)
                                        color=17;
 
                 else if(topo->value[i][j] < 320*topo_scale)
                                        color=24; 
 
                 else
                                        color=31;
     
               tdi=-topo->value[i][j+1]/3+topo->value[i+1][j]/3;
               if(tdi > 3)
                        tdi=3;
 
                 if(tdi < -3)
                        tdi=-3;
 
                 color=color+tdi;
         
                 color=tmap[color];          
/*     
                 if(color==0)
                        continue;
*/
              XSetForeground(display,gc,color);
 
              if(zoom <= 1) {
 
/* fix for wider gaps */
 
              for(k=0;k<4;k++) {               
 
                  nx[k]=xmin+(float)topo->coord[i+di[k]][j+dj[k]].x*cmult;                
                  ny[k]=ymin+(float)topo->coord[i+di[k]][j+dj[k]].y*cmult;
 
        }
/*
 
                  nx[0]=xmin+(float)topo->coord[i][j].x*cmult;                
                  ny[0]=ymin+(float)topo->coord[i][j].y*cmult;
 
                  nx[1]=xmin+(float)topo->coord[i+1][j].x*cmult;                
                  ny[1]=ymin+(float)topo->coord[i+1][j].y*cmult;
 
                  nx[2]=xmin+(float)topo->coord[i+1][j+1].x*cmult;                
                  ny[2]=ymin+(float)topo->coord[i+1][j+1].y*cmult;
 
                  nx[3]=xmin+(float)topo->coord[i][j+1].x*cmult;                
                  ny[3]=ymin+(float)topo->coord[i][j+1].y*cmult;
 
 
*/
 
        }
 
              else {
  
                     for(k=0;k<4;k++) {
 
                     nx[k]=(xmin+((float)topo->coord[i+di[k]][j+dj[k]].x*cmult-
                           lint1))*zoom + xcen;
                     ny[k]=(ymin+((float)topo->coord[i+di[k]][j+dj[k]].y*cmult-
                           lint2))*zoom + ycen;
           }
                                      
 
 
           }
 
          if(zoom>=1) {
   
                 for(k=0;k<4;k++) {
 
                        point[k].x=(nx[k]+5)/MY;
                        point[k].y=(ny[k]+5)/MY;
 
 
              }
              
                 point[4].x=point[0].x;
                 point[4].y=point[0].y;
/*
                 for(xx=point[0].x;xx<=point[1].x;xx++){
 
                 for(yy=point[0].y;yy<=point[3].y;yy++){
 
                    XPutPixel(image,xx,yy,color);
 
                  }
           }
*/
 
 
                 XFillPolygon(display,pixm,gc,point,5,Convex,CoordModeOrigin);
 
           }
                         
 
              else  {
 
                                                            
 
 
 
          }
 
            
                 }
 
 
}
 
return;
 
}