利用CE318太阳光度计数据处理反演AOD和水汽含量wv

利用太阳光度计CE318进行基本的数据处理，反演气溶胶光学厚度和水汽含量

#include 
#include 
//#define M 179

#define TMLOGPATH     “E:\\DTF5\\CE318data_software\\x081111.txt“
#define V0RSDPATH     “E:\\DTF5\\CE318data_software\\V0_R_SD.txt“
#define DEPTHPATH     “E:\\DTF5\\CE318data_software\\AODx081111.txt“
#define ALLPATH       “E:\\DTF5\\CE318data_software\\allx081111.txt“

#define ROW 17
void main（）
	{
    	int ijM=0;
		FILE *p;
		char chstr[10]=“x081111“;
	    double timal[ROW]m[ROW]V[ROW][5]mavga[5]b[5]lat=31.9037lon=117.1584;//大营盘lat=40.255lon=115.8;//合肥lat=31.9037lon=117.1584;风电场lat=41.017lon=114.758//经纬度注意是否要修改
		double Vavg[5]V0[5]total=0.0sumV[5]dx=0.0dy[5]dxy[5]R[5]SD[5];
		double bm[ROW]bmavgVm[ROW]Vmavg;
		double tao[ROW][5]taor[5]taoa8[ROW][5]taoa940[ROW]tao940[ROW]taoH2O[ROW]dxy940=0.0ab=0.606ba=0.580;//Wb[M]W[M]
		double taoa5[ROW][4]lnnamta5[4]namta5avgsumtaoa5[ROW]taoa5avg[ROW]dxym[ROW]slope1[ROW]arfa[ROW]beta[ROW];
		double namta8[5]={1.0200.9360.8700.6700.440};
			
		for（i=0;i<5;i++）
		{
		        dxy[i]=0.0;
				dy[i]=0.0;
		        sumV[i]=0.0;
				SD[i]=0.0;
		}

		for（i=0;i		{
			dxym[i]=0.0;
			sumtaoa5[i]=0.0;
		}


		if（（p=fopen（TMLOGPATH“r“））==NULL）
		{
			printf（“can‘t open the file!\n“）;
			return;
		}
		ch=fgetc（p）;
		while（ch!=EOF）
		{
			if（ch==‘\n‘） M++;
			ch=fgetc（p）;
		}
		fclose（p）;

		printf（“%d\n“M）;


		if（（p=fopen（TMLOGPATH“r“））==NULL）
		{
			printf（“can‘t open the file!\n“）;
			return;
		}
		for（i=0;i		{
			fscanf（p“%lf %lf“&timal[i]&m[i]）;
			for（j=0;j<5;j++）
				fscanf（p“%lf“&V[i][j]）;
		}
		fclose（p）;

				
	//计算大气层外仪器测量值V0及定标因子C


	for（i=0;i		total=total+m[i];
	mavg=total/M;
	//printf（“%lf\n“mavg）;
	for（i=0;i		for（j=0;j<5;j++）
	    	sumV[j]=sumV[j]+V[i][j];
	for（j=0;j<5;j++）
	{
		Vavg[j]=sumV[j]/M;
		//printf（“%lf\t“exp（Vavg[i]））;
	}
	for（i=0;i         dx=dx+pow（（m[i]-mavg）2）;
    for（i=0;i<5;i++）
		for（j=0;j		{
			dy[i]=dy[i]+pow（（V[j][i]-Vavg[i]）2）;
			dxy[i]=dxy[i]+（m[j]-mavg）*（V[j][i]-Vavg[i]）;
		}
	for（i=0;i<5;i++）
		{
			a[i]=1.0*dxy[i]/dx;//斜率
			b[i]=Vavg[i]-a[i]*mavg;//截距
			V0[i]=exp（b[i]）;
			R[i]=dxy[i]/sqrt（dx*dy[i]）;
		}
	for（i=0;i<5;i++）
		for（j=0;j			SD[i]=SD[i]+pow（（V[j][i]-a[i]*m[j]-b[i]）2）;
	for（i=0;i<5;i++）
		SD[i]=sqrt（SD[i]/M）;
	for（i=0;i<5;i++）
		printf（“%lf\t%lf\t%lf\n“exp（b[i]）R[i]SD[i]）;

	
	//水汽波段定标	
    //计算大气柱总光学厚度
	for（i=0;i		for（j=0;j<5;j++） 
			V[i][j]=exp（V[i][j]）;
	for（i=0;i		for（j=0;j<5;j++） 
		    tao[i][j]=log（V0[j]/V[i][j]）/m[i];
		
	for（i=0;i<5;i++）
	{
		taor[i]=0.008569*pow（namta8[i]-4）*（1+0.0113*pow（namta8[i]-2）+0.00013*pow（namta8[i]-4））*1000*exp（-0.125017*0.044）/1013.25;  //计算大气分子垂直光学厚度taor[i]=0.0088*pow（namta8[i]-4.05）
		//printf（“%lf\t“taor[i]）;
	}
		
	for（i=0;i		for（j=0;j<5;j++）
			taoa8[i][j]=tao[i][j]-taor[j];  //计算气溶胶光学厚度，含水汽
	
		//利用线性内

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-06-13 14:33  sun_photometer_ce318_data_deal\
     文件        2602  2012-12-27 15:06  sun_photometer_ce318_data_deal\AOD_CMP.CPP
     文件        6540  2012-12-27 15:10  sun_photometer_ce318_data_deal\dingbiao.cpp
     文件        3701  2012-12-27 15:09  sun_photometer_ce318_data_deal\Fit_Ang.cpp
     文件        3742  2018-06-13 14:31  sun_photometer_ce318_data_deal\read_ce318_data.cpp
     文件        2723  2012-12-27 15:05  sun_photometer_ce318_data_deal\WAT_CMP.CPP