Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include <stdio.h>
- #include <stdlib.h>
- #include <string.h>
- #include <time.h>
- int main(int argc,char *argv[])
- {
- FILE *fpi1,*fpi2,*fpo;
- int i,j,width;
- int k,l;
- int height1,height2;
- //unsigned char dn1,dn2;
- float dn1;
- int dn2;
- float dn1_x,dn1_y,dn1_z;
- float dn3;
- //unsigned char *file1,*file2;
- float *file1;
- int *file2;
- //float *file3;
- //float start,end;
- //start=clock();
- if(argc!=7)
- {
- fprintf(stderr,"Usage: %s\n(2)Input_R_READ_filename\n(3)Input_NIR_READ_filename\n(4)Output_RVI_filename\n(5)width\n(6)height\n(7)header_size\n",argv[0]);
- exit(1);
- }
- printf("OPEN Points FILE NAME:%s\n",argv[1]);
- printf("OPEN Verteices FILE NAME:\%s\n",argv[2]);
- printf("OPEN OUTPUT FILE NAME:%s\n",argv[3]);
- width = atoi(argv[4]); printf("W = %d,", width);
- height1 = atoi(argv[5]); printf("H1 = %d\n", height1);
- height2 = atoi(argv[6]); printf("H2 = %d\n", height2);
- //header=atoi(argv[6]);printf("header=%d\n",header);
- /*---------ÉÅÉÇÉäÇÃé©ìÆäÑìñ----------*/
- file1= (float *)malloc(width*height1*sizeof(float));
- file2= (int *)malloc(width*height2*sizeof(int));
- //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((float *)file1,sizeof(float),width*height1,fpi1);
- /*----------------âÊëúÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉìèIóπ----------------------*/
- /*----------------ì¸óÕNIRâÊëúÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉìäJén----------------------*/
- if((fpi2=fopen(argv[2],"rb"))==NULL)
- { /*ÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉì*/
- fprintf(stderr,"output file open error\n");
- exit(1);
- }
- fread((int *)file2,sizeof(int),width*height2,fpi2);
- /*----------------âÊëúÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉìèIóπ----------------------*/
- /*----------------èoóÕRVIâÊëúÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉìäJén----------------------*/
- if((fpo=fopen(argv[3],"wb"))==NULL)
- { /*ÉtÉ@ÉCÉãÇÃÉIÅ[ÉvÉì*/
- fprintf(stderr,"output file open error\n");
- exit(1);
- }
- //fread((float *)file3,sizeof(float),width*height,fpo);
- for(i=0;i<height2;i++)
- {
- for(j=0;j<width;j++)
- {
- dn2=file2[i*width+j];//vertex No
- dn1_x=file1[dn2*width];
- printf("%f\n",dn1_x);
- fwrite(&dn1_x,sizeof(float),1,fpo); /*míPà ïWçÇÉfÅ[É^ÇÃèëÇ´çûÇ›*/
- dn1_y=file1[dn2*width+1];
- printf("%f\n",dn1_y);
- fwrite(&dn1_y,sizeof(float),1,fpo); /*míPà ïWçÇÉfÅ[É^ÇÃèëÇ´çûÇ›*/
- dn1_z=file1[dn2*width+2];
- printf("%f\n",dn1_z);
- fwrite(&dn1_z,sizeof(float),1,fpo); /*míPà ïWçÇÉfÅ[É^ÇÃèëÇ´çûÇ›*/
- }
- }
- free(file1);
- free(file2);
- //free(file3);
- fclose(fpi1);
- 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