FAI计算程序

欢迎共阅欢迎共阅欢迎共阅一、所需数据每日的500米分辨率的HDF文件二、几何校正ENVI下File->OpenExternalFile->EOS->MODIS->打开500米的HDF文件->找到太湖->右键ROIMap->GeoreferenceMODIS->对辐亮度和反射率数据分别进行几何校正->编辑spatialsubset->ROI/EVF->选择区域->OK在GeoreferenceMODISparameters下，投影选择Taihu->OK->输出：日期taihu_radiance_proj日期taihu_reflect_proj三、红光-近红外-短波红外辐亮度值针对MODIS数据，对辐亮度数据，1、2、5波段分别查看每个波段的平均辐亮度值，记录下来，在使用6S模型时使用。也可记录下1、2、5波段的反射率值（那样是话，在6S中RAPP那项前面需加上负号）。如下图四、运行6S模型参数如下：GEOMETRICALCONDITIONS:---------------------igeom(0-7)=0传感器类型YOUROWNCONDITION:---------------------25.44127.9114.0498.6765太阳天顶角与方位角，卫星天顶角与方位角，日期ATMOSPHERICMODEL:---------------------idatm(0-8)=2大气模式。2代表MidlatitudeSummerAEROSOLMODEL(TYPE):---------------------iaer(0-12)=1气溶胶类型。1代表ContinentalAEROSOLMODEL(CONCENTRATION):------------------------------v(visibility)=15.0能见度，需要查阅当地的气象资料xps(inkm)=-0.120当地的海拔，加上-号xpp=-705传感器的高度，加上-号ENTERTHESPECTRALCONDITIONS:------------------------------iwave(-1to34)=42传感器接收波段，如MODIS为42-48GROUNRDREFLECTANCE(TYPE)--------------------------inhomo(0-1)=0是否为均质地物，0均一，1，不均一idirec=0是否考虑二向性，0表示没有，1表示有igroun=4地物类型，如1表示植被，4表示湖泊rapp=48.06356RAPP=gias*DN+bias，大气层顶太阳辐亮度（1、2、5波段分别计算）分别计算1、2、5波段，6S的结果中，在第3页找到计算的瑞利散射反射率值，记录下来五、计算FAI=1\*GB3①ENVI下打开反射率数据=2\*GB3②计算RTOA-R6S，R6S为1、2、5波段分别用6S计算出来的值。BasicTools下找到BandMath输入指定b1，命名为：日期band1_Rrc同样的方法计算band2_Rrc，band5_Rrc=3\*GB3③输入下列，计算band2_Rrc’在这里，b1指定为band1_Rrc，b5指定为band5_Rrc，输出为band2_Rrc’=4\*GB3④计算FAI使用b1-b2，在这里，b1指定为band2_Rrc，b2指定为band2_Rrc’，输出为FAI五、使用FAI提取蓝藻此过程同蓝藻监测中的步骤，先使用模型（ampm），进行分区，再使用决策树设定阈值进行。