plot_isoh.c

Go to the documentation of this file.

 
#include "prototypes.h"
 
 
void plot_isoh(Pixmap pixm,int h, int display_flag)
 
{
 
extern int isom;
extern int iscale;
extern Display *display;
extern int npoint[2],zoom;
extern GC gc;
extern struct isoh *isoh;
extern unsigned long imap[16];
extern Widget drawing_area;
extern struct display_set display_set[4];
int i,j,k,m;
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 ceni,cenj;
float ddi;
float fm;
 
  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;
 
 if(zoom==1 && (isoh->maxi > 500 || isoh->maxj > 500)){
 
      dx=2;dy=2;
      di[0]=0;dj[0]=0;
      di[1]=2;dj[1]=0; 
      di[2]=2;dj[2]=2;
      di[3]=0;dj[3]=2;
 
    }
 
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=isoh->maxi;
maxj=isoh->maxj;
ceni=maxi/2;
cenj=maxj/2;
 
if(zoom > 1) {
 
        iflag=0;
 
        for(i=0;i<maxi-dx;i++) {
 
            nlong=(xmin+((float)isoh->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)isoh->coord[ceni][j].y*cmult-lint2))*zoom + ycen;
            
 
             if(iflag==0 && nlong < ymin)
                        continue;
 
             if(iflag==1 && nlong < ymax)
                        continue;
 
             if(iflag==0)
                     minj=j;
 
 
             if(iflag==1)
                     maxj=j;
 
             iflag++;
 
       }
 
      }
 
 
if(iscale==0)
      fm=.5;
else
      fm=pow(2,(double)(iscale-1));
 
 
printf("isoh %d %d %f\n",xmin,ymin,cmult);
   
for(i=mini;i<maxi-dx;i=i+dx) {
 
for(j=minj;j<maxj-dy;j=j+dy) {
 
                  if(isoh->value[isom][i+di[0]][j+dj[0]] < 0 ||
                     isoh->value[isom][i+di[1]][j+dj[1]] < 0 ||
                     isoh->value[isom][i+di[2]][j+dj[2]] < 0 ||
                     isoh->value[isom][i+di[3]][j+dj[3]] < 0) 
                                                    continue;
 
                  ddi=(float)(isoh->value[isom][i+di[0]][j+dj[0]]+
                              isoh->value[isom][i+di[1]][j+dj[1]] +
                              isoh->value[isom][i+di[2]][j+dj[2]]+
                              isoh->value[isom][i+di[3]][j+dj[3]])/(4*25.3);
 
                    if(ddi <= 0)
                              continue;
 
                    for(h=1;h <= 16;h++) {
 
                           if(ddi < fm*h ) {
 
                                        color=h;
                                        break;
 
                      }
 
                                  }        
 
                    if(color > 16)
                                color=16;        
 
                    color=imap[color-1];
 
              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)isoh->coord[i+di[k]][j+dj[k]].x*cmult;                
                  ny[k]=ymin+(float)isoh->coord[i+di[k]][j+dj[k]].y*cmult;
 
        }
        }
 
              else {
  
                     for(k=0;k<4;k++) {
 
                     nx[k]=(xmin+((float)isoh->coord[i+di[k]][j+dj[k]].x*cmult-
                           lint1))*zoom + xcen;
                     ny[k]=(ymin+((float)isoh->coord[i+di[k]][j+dj[k]].y*cmult-
                           lint2))*zoom + ycen;
 
                                      } 
 
    }      
              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;
 
              XFillPolygon(display,pixm,gc,point,5,Convex,CoordModeOrigin);
                                     
                 }
 
 
}
 
return;
}