points-polygons2

/*Points-Polygon2 program
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,*fpi2,*fpo;
        int i,j,width;
    int k,l;
    int height1,height2;
    //unsigned char dn1,dn2;
    float dn1;
    int dn2;
    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];
            fwrite(&dn1_x,sizeof(float),1,fpo); /*míPà ïWçÇÉfÅ[É^ÇÃèëÇ´çûÇ›*/
            dn1_y=file1[dn2*width+1];
            fwrite(&dn1_y,sizeof(float),1,fpo); /*míPà ïWçÇÉfÅ[É^ÇÃèëÇ´çûÇ›*/
                dn1_z=file1[dn2*width+2];
                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);

}