plot_topo_new.c

Go to the documentation of this file.

 
#include "prototypes_new.h"
 
 
void plot_topo(Pixmap pixm,int h, int display_flag)
 
{
 
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;
extern int max_width,max_height;
long zlint1,zlint2;
float zmult;
long xfinal,yfinal;
 
for(i=0;i<max_width;i++) {
    for(j=0;j<max_height;j++) {
        XPutPixel(image, i, j, XBlackPixel(display, 0));
    }
}
 
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;
zlint1=-lint1*zoom;
zlint2=-lint2*zoom;
zmult=cmult*zoom;
xfinal=xmin*zoom+zlint1+xcen;
yfinal=ymin*zoom+zlint2+ycen;
 
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=(float)topo->coord[i][cenj].x*zmult + xfinal;
           
 
            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=(float)topo->coord[ceni][j].y*zmult + yfinal;
 
             if(iflag==0 && nlong < 0)
                        continue;
 
             if(iflag==1 && nlong < YSIZE)
                        continue;
                        
             if(iflag==0)
                     minj=j;
 
 
             if(iflag==1)
                     maxj=j;
 
             iflag++;
 
       }
 
      }
printf("+++++++++++++++++ %d %d %d %d %d %d", max_width, max_height, mini, minj, maxi, maxj);
fflush(stdout);
 
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;
        }
 
        color=(topo->value[i][j]/80)*7 + 3;
               
        if(color > 31)
            color=31;
                                         
     
        tdi=-topo->value[i][j+1]/3+topo->value[i+1][j]/3;
               
        if(tdi > 3)
            tdi=3;
        else if(tdi < -3)
            tdi=-3;
 
        color=color+tdi;
         
        color=tmap[color];          
 
        XSetForeground(display,gc,color);
 
        if(zoom <= 1) {
            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;
            }
        }
        else {
            for(k=0;k<4;k++) {
                nx[k]=(float)topo->coord[i+di[k]][j+dj[k]].x*zmult + xfinal;
                ny[k]=(float)topo->coord[i+di[k]][j+dj[k]].y*zmult + yfinal;
            }
        }
           
              
        if(zoom>=1) {
            for(k=0;k<4;k++) {
                point[k].x=(nx[k]+5)/MY;
                point[k].y=(ny[k]+5)/MY;
                if(point[k].x < 0 || point[k].x >= max_width ||
                    point[k].y < 0 || point[k].y >= max_height)
                    break;      
                }
              
                if(k != 4)
                    continue;
              
              
                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);
                    }
                }
            }
            else  {
 
        }
    }
}
 
XPutImage(display,pixm,gc,image,0,0,0,0,width,height);
 
return;
 
}