Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- /*rvi program
- band3:
- http://aitech.ac.jp/~y-yama/data/rvi/band3.raw.zip
- band4:
- http://aitech.ac.jp/~y-yama/data/rvi/band4.raw.zip
- 512*512
- unsigned char
- made by Yoshiyuki Yamamoto
- ver1.0 2004/5/30*/
- #include <stdio.h>
- #include <stdlib.h>
- #include <string.h>
- #include <time.h>
- int main(int argc,char *argv[])
- {
- FILE *fpi1,*fpw,*fpo;
- int i,j,height,width;
- unsigned int dn1,dn2;
- unsigned int sub_dn;
- //float dn3;
- unsigned char line;
- unsigned int *file1;
- //float *file3;
- float start,end;
- start=clock();
- if(argc!=5)
- {
- fprintf(stderr,"Usage: %s\n(1)Input_READ_Binary_filename\n(2)width\n(3)height\n(4)output_filename\n",argv[0]);
- exit(1);
- }
- printf("OPEN R FILE NAME:%s\n",argv[1]);
- //printf("OPEN NIR FILE NAME:\%s\n",argv[2]);
- //printf("OPEN OUTPUT FILE NAME:%s\n",argv[3]);
- width = atoi(argv[2]); printf("W = %d,", width);//1
- height = atoi(argv[3]); printf("H = %d\n", height);
- //header=atoi(argv[6]);printf("header=%d\n",header);
- /*---------ÉÅÉÇÉäÇÃé©ìÆäÑìñ----------*/
- file1= (unsigned int *)malloc(width*height*sizeof(int));
- //file2= (unsigned char *)malloc(width*height*sizeof(char));
- //file3= (float *)malloc(width*height*sizeof(float));
- /*----------------ì¸óÕRâÊëúÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉìäJén----------------------*/
- if((fpi1=fopen(argv[1],"rb"))==NULL)
- { /*ÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉì*/
- fprintf(stderr,"input file open error\n");
- exit(1);
- }
- fread((unsigned int *)file1,sizeof(int),width*height,fpi1);
- /*----------------âÊëúÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉìèIóπ----------------------*/
- /*----------------ì¸óÕNIRâÊëúÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉìäJén----------------------*/
- if((fpw=fopen(argv[4],"w"))==NULL)
- {
- fprintf(stderr,"DSM by GSI データファイル.raw\n");
- exit(1);
- }
- for(i=0;i<height-1;i++)
- {
- for(j=0;j<width;j++)
- {
- dn1=file1[i*width+j];
- dn2=file1[(i+1)*width+j];
- sub_dn=dn2-dn1;
- fprintf(fpw,"%d\n",sub_dn);
- //printf("dn1=%d ",dn1);
- //if(j==0 && dn1==1)
- //{
- //printf("line=%d\n",i);
- //fprintf(fpw,"%d\n",i);
- //}
- }
- }
- free(file1);
- //free(file2);
- //free(file3);
- fclose(fpi1);
- fclose(fpw);
- //fclose(fpi2);
- //fclose(fpo);
- end=clock();
- //printf("äJénéûä‘ÇÕÅC%fïb\n",start/1000);
- //printf("èIóπéûä‘ÇÕÅC%fïb\n",end/1000);
- //printf("åoâfléûä‘ÇÕÅC%fïb\n",(end-start)/1000);
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement