Advertisement
yamaji14

how many points of polygon

Jan 26th, 2015
234
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 2.56 KB | None | 0 0
  1. /*rvi program
  2.  
  3. band3:
  4. http://aitech.ac.jp/~y-yama/data/rvi/band3.raw.zip
  5. band4:
  6. http://aitech.ac.jp/~y-yama/data/rvi/band4.raw.zip
  7. 512*512
  8. unsigned char
  9.  
  10. made by Yoshiyuki Yamamoto
  11. ver1.0 2004/5/30*/
  12.  
  13.  
  14. #include <stdio.h>
  15. #include <stdlib.h>
  16. #include <string.h>
  17. #include <time.h>
  18.  
  19. int main(int argc,char *argv[])
  20. {
  21.     FILE *fpi1,*fpw,*fpo;              
  22.         int i,j,height,width;  
  23.     unsigned int dn1,dn2;
  24.     unsigned int sub_dn;
  25.     //float dn3;
  26.     unsigned char line;
  27.     unsigned int *file1;
  28.     //float *file3;          
  29.     float start,end;
  30.    
  31.     start=clock();
  32.  
  33.     if(argc!=5)
  34.     {
  35.         fprintf(stderr,"Usage: %s\n(1)Input_READ_Binary_filename\n(2)width\n(3)height\n(4)output_filename\n",argv[0]);
  36.         exit(1);
  37.     }
  38.  
  39.         printf("OPEN R FILE NAME:%s\n",argv[1]);
  40.     //printf("OPEN NIR FILE NAME:\%s\n",argv[2]);
  41.     //printf("OPEN OUTPUT FILE NAME:%s\n",argv[3]);
  42.  
  43.     width = atoi(argv[2]); printf("W = %d,", width);//1
  44.         height = atoi(argv[3]); printf("H = %d\n", height);
  45.     //header=atoi(argv[6]);printf("header=%d\n",header);
  46.  
  47.     /*---------ÉÅÉÇÉäÇÃé©ìÆäÑìñ----------*/
  48.         file1= (unsigned int *)malloc(width*height*sizeof(int));
  49.     //file2= (unsigned char *)malloc(width*height*sizeof(char));
  50.     //file3= (float *)malloc(width*height*sizeof(float));
  51.  
  52.     /*----------------ì¸óÕRâÊëúÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉìäJén----------------------*/
  53.         if((fpi1=fopen(argv[1],"rb"))==NULL)
  54.     {       /*ÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉì*/
  55.             fprintf(stderr,"input file open error\n");
  56.         exit(1);
  57.         }
  58.  
  59.         fread((unsigned int *)file1,sizeof(int),width*height,fpi1);
  60.     /*----------------âÊëúÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉìèIóπ----------------------*/
  61.     /*----------------ì¸óÕNIRâÊëúÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉìäJén----------------------*/
  62.  
  63.     if((fpw=fopen(argv[4],"w"))==NULL)
  64.         {
  65.                 fprintf(stderr,"DSM by GSI データファイル.raw\n");
  66.                 exit(1);
  67.         }
  68.  
  69.     for(i=0;i<height-1;i++)
  70.     {
  71.         for(j=0;j<width;j++)
  72.         {
  73.             dn1=file1[i*width+j];
  74.             dn2=file1[(i+1)*width+j];
  75.             sub_dn=dn2-dn1;
  76.             fprintf(fpw,"%d\n",sub_dn);
  77.             //printf("dn1=%d ",dn1);
  78.            
  79.             //if(j==0 && dn1==1)
  80.             //{
  81.                 //printf("line=%d\n",i);
  82.                 //fprintf(fpw,"%d\n",i);
  83.             //}
  84.         }
  85.     }
  86.  
  87.     free(file1);
  88.     //free(file2);
  89.     //free(file3);
  90.     fclose(fpi1);
  91.     fclose(fpw);
  92.     //fclose(fpi2);
  93.     //fclose(fpo);
  94.     end=clock();
  95.  
  96.     //printf("äJénéûä‘ÇÕÅC%fïb\n",start/1000);
  97.     //printf("èIóπéûä‘ÇÕÅC%fïb\n",end/1000);
  98.     //printf("åoâfléûä‘ÇÕÅC%fïb\n",(end-start)/1000);
  99.  
  100. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement