/* 
  convert CCT count to real data (for WV, IR1&2) 
  for: gcc (no need options) 
  made by: Atsushi Higuchi, CEReS, Chiba Univ. Japan

Version history:
ver.0.1 31 mar 2007: start make program
*/
#include <stdio.h>
#include <stdlib.h>
#define	 data 2096704
#define	 bit 256 
FILE	*fpi, *fpo, *fpc;
/* open file sub-rutine */
int	open_file(inname, outname, calibname)
char	*inname, *outname, *calibname;
{
	if((fpi = fopen(inname, "rb")) == NULL){
		printf("don't open %s!\n", inname);
		return(-1);
	}
	if((fpo = fopen(outname, "wb")) == NULL){
		printf("don't write %s!\n", outname);
		return(-1);
	}
        if((fpc = fopen(calibname, "rt")) == NULL){
                printf("don't read %s!\n", calibname);
                return(-1);
	}
return(0);
}
/* main */
int	main(argc, argv)
int	argc;
char	*argv[];
{
int	i,dummy;
long	j;
unsigned char input[data];
float	output[data];
float	calib[bit];
	if(argc !=4){
	printf("Usage: program [input] [output] [calibration file]\n", argv[0]);
	exit(-1);
	}
	open_file(argv[1], argv[2], argv[3]);
	/* read original CCT count data */
	fread (input, sizeof(unsigned char), data, fpi);
	/* read calibration data */
	for(i=0;i<bit;i++){
		fscanf(fpc, "%d %f", &dummy, &calib[i]);
	}
	/* conversion from CCT to real data */
	for(j=0;j<data;j++){
		for(i=0;i<bit;i++){
			if(input[j]==i){
			output[j]=calib[i];
			}
		}
	}
	fwrite(output, sizeof(float), data, fpo);
fclose(fpi);fclose(fpo);fclose(fpc);
return(0);
}