/*------------------------------------------------------------------------------
	NDVIの各種パラメータを求める  [GIMMSバージョン]
        0:年最大値
        1:積算（0.1以上)
  ------------------------------------------------------------------------------*/

/*　格納されているデータの値域（要確認）

        NDVI 0-1 : 0-1000（ただし、負の値もあり）
        Water -0.1 : -10000
        Null -0.05 : （未確認）
*/

#include        <stdio.h>
#include        <stdlib.h>
#include        <fcntl.h>
#include        <sys\types.h>
#include        <sys\stat.h>
#include        <process.h>
#include        <dos.h>
#include        <io.h>
#include        <share.h>
#include        <math.h>
#include        <string.h>
#include        <ctype.h>
#include        <conio.h>
#include        <errno.h>

#define         NUMBER_OF_LINES		2091
#define         NUMBER_OF_PIXELS        4950
#define         NUMBER_OF_CHANNELS	24
#define         NUMBER_OF_PRM           2

signed short    ndvi[NUMBER_OF_CHANNELS][NUMBER_OF_PIXELS];
signed short    nprm[NUMBER_OF_PRM][NUMBER_OF_PIXELS];
unsigned char   buffer[NUMBER_OF_PIXELS*2];

void    main(argc,argv)
int     argc;
char    *argv[];
{
	int	fi,fo;
	int	i,j,ch;
	int	isw;    /* 0:最大値、1:積算 */
        signed short  max,count;
        int     sum;

	if(argc!=4){
                printf("Extraction of NDVI Parameter - GIMMS version\n");
		printf("USAGE: getviprm [input] [output] [switch]\n");
                printf("(1) input  : GIMM NDVI - annual 24 bi-monthly data. BIL format.\n");
                printf("(2) output : 0:max, 1:sum(>0.1)\n");
                printf("note) output values are NDVI*1000\n");
                exit(-1);
	}
 
        if( (fi=open(argv[1],O_RDONLY | O_BINARY))==-1){
                       printf("Can't open %s.\n",argv[1]);
                       exit(-1);
        }

        if( (fo=open(argv[2],O_RDWR | O_TRUNC | O_CREAT | O_BINARY,
                        S_IREAD | S_IWRITE))==-1){
                       printf("Can't open %s.\n",argv[2]);
                       exit(-1);
        }

	isw=atoi(argv[3]);
	if(isw<0 || isw>1){
		printf("Switch error %d\n",isw);
		exit(-1);
	}

/* 計算 */

        for(i=0;i<NUMBER_OF_LINES;i++){
                if( !(i%100) ) printf("*");

/* 24半月読みとり [BIL] */

                for(ch=0;ch<NUMBER_OF_CHANNELS;ch++){
                        read(fi,ndvi[ch],NUMBER_OF_PIXELS*sizeof(short));
                }

/* isw ... 0:最大値、1:積算(>0.1) */

                for(j=0;j<NUMBER_OF_PIXELS;j++){

                        if(ndvi[0][j]<-9998){     /* 領域外(水域):-9999 */
                                nprm[0][j]=-9999;
                                nprm[1][j]=-9999;
                                continue;
                        }

                        sum=0;
                        max=-9999;
                        count=0;
                        for(ch=0;ch<NUMBER_OF_CHANNELS;ch++){ /* 半月スキャン */

                                if(ndvi[ch][j]>max){
                                        max=ndvi[ch][j];
                                }

                                if(ndvi[ch][j]>100 ){ /* NDVI 0.1=100 植生域のみ積算 */
                                        sum += ndvi[ch][j]; /* 1000倍の値をそのまま積算 */
                                        count++;
                                }
                        }/* next ch */

                     /* 最大値格納 isw=0 */

                        if(isw==0){
                                nprm[0][j]=max;
                        }

                     /* 積算値格納 isw=1*/

                        if(isw==1){
                                nprm[1][j]=sum;
                        }

                }/* next j */

/* 出力 NDVI*1000 */

                write(fo,nprm[isw],NUMBER_OF_PIXELS*sizeof(short));

	}/* next line */

        printf("\n");

}/* end of main */