get_topo_coord_new.c

Go to the documentation of this file.

#include "prototypes.h"
 
void get_topo_coord(char *fname)
 
{
 
 extern double max_lat,min_lat,cen_lon;
 double max_lon;
 int i,j,ier,maxj,maxi,max_value,k;
 double lat,lon;
 extern struct topo *topo;
 FILE *fp;
 char *p,kbuf[10000];
 double mlat,mlon,delta_lat,delta_lon,total_lat,total_lon;
 int ii,jj;
 double xlat,xlon;
 int taxi,taxj;
 
 topo=(struct topo *) calloc(1,sizeof(struct topo));
 if(topo==NULL) {
 
                   printf("could not allocate topo space\n");
                   exit(1);
 
                  }          
 
 max_value=0;
 
 fp=fopen(fname,"r");
 
 if(fp==NULL) {
 
              printf("could not open %s\n",fname);
              exit(1);
 
              }
 
 p=fgets(kbuf,80,fp);
 
 ier=sscanf(kbuf,"%lf %lf %lf %lf %lf %lf",&mlat,&mlon,&total_lat,&total_lon,&delta_lat,&delta_lon);
 
 delta_lat=delta_lat/60.;
 delta_lon=delta_lon/60.;
 
 taxi=total_lon/delta_lon;
 taxj=total_lat/delta_lat;
 
 maxi=10/delta_lon+1;
 maxj=(max_lat-min_lat)/delta_lat+1;
 
 max_lon=cen_lon+5;
 
 topo->maxi=maxi;
 topo->maxj=maxj;
 topo->max_lat=max_lat;
 topo->max_lon=cen_lon+5;
 topo->total_lat=max_lat-min_lat;
 topo->total_lon=10;
 topo->delta_lat=delta_lat;
 topo->delta_lon=delta_lon;
 
 topo->coord=(struct coord**) calloc(maxi,sizeof(struct coord *));
 topo->value=(short int**) calloc(maxi,sizeof(short int *));
 
 if(topo->coord==NULL || topo->value==NULL) {
 
                            printf("no memory for topo array\n");
                            exit(1);
 
              }
 
 for(i=0;i<maxi;i++) {
 
      topo->coord[i]=(struct coord *) calloc(maxj,sizeof(struct coord));
 
      topo->value[i]=(short int *) calloc(maxj,sizeof(short int));
 
      if(topo->coord[i]==NULL || topo->value[i]==NULL) {
 
                printf("no memory for topo array\n");
                exit(1);
 
          }
 
    }
 
 for(i=0;i<maxi;i++) {
 for(j=0;j<maxj;j++) {
 
topo->value[i][j]=-9999; 
 
}
}
 
 for(j=0;j<taxj;j++) {
 
 ier=fread(kbuf,sizeof(char),taxi*4,fp);
 
 xlat=mlat-(double)j*delta_lat;
 
 jj=((max_lat-xlat)/delta_lat+.0000001);
 
 if(jj < 0 || jj > maxj-1)
                    continue; 
 
 k=0;
 
 
 
 for(i=0;i<taxi;i++) {
 
         xlon=mlon-(double)i*delta_lon;
         
         ii=((max_lon-xlon)/delta_lon+.0000001);
         
         if(ii < 0 || ii > maxi-1) {
         
                            k=k+4;
         
                            continue; 
         
                            } 
                           
         topo->value[ii][jj]=atoi(&kbuf[k]);
 
         lat=mlat-(float)j*delta_lat;
         lon=mlon-(float)i*delta_lon;
 
         get_topo_pixels(ii,jj,lat,lon);
 
         k=k+4;         
 
    }
 
   }
 
fclose(fp);
 
return;
 
}
void get_topo_pixels(int i, int j, float lat,float lon)
 
  {
 
     extern struct dval dval;
     long r,s;
     float conv=.0174;
     extern struct topo *topo;
     
     r=dval.a * cos(lat*conv)/(1+sin(lat*conv))
       * cos((lon-dval.lo-90)*conv) + dval.xo +.5;
 
     topo->coord[i][j].x=r;
 
     s=dval.a * cos(lat*conv)/(1+sin(lat*conv))
       * sin((lon-dval.lo-90)*conv) + dval.yo + .5;
       
     topo->coord[i][j].y=s;
 
     return;
 
   }