2022年FAI计算流程.doc

一、所需数据每日的500米分辨率的 HDF文件二、几何校正ENVI下 File -> Open External File ->EOS -> MODIS->打开500米的HDF文件->找到太湖 ->右键ROI Map-> Georeference MODIS ->对辐亮度和反射率数据分别进展几何校正 ->编辑spatial subset->ROI/EVF->选择区域->OK 在Georeference MODIS parameters下，投影选择Taihu->OK->输出：日期taihu_radiance_proj日期 taihu_reflect_proj三、记录红光-近红外-短波红外辐亮度值 针对MODIS数据，对辐亮度数据，1、2、5波段分别查看每个波段的平均辐亮度值，记录下来，在使用6S模型时使用。也可记录下1、2、5波段的反射率值（那样是话，在6S中RAPP那项前面需加上负号）。如以下图 四、运转6S模型参数如下：GEOMETRICAL CONDITIONS:-igeom (0-7)= 0 传感器类型YOUR OWN CONDITION:-25.44 127.91 14.04 98.67 6 5 太阳天顶角与方位角，卫星天顶角 与方位角，日期 ATMOSPHERIC MODEL:-idatm (0-8) = 2 大气形式。2代表Midlatitude SummerAEROSOL MODEL (TYPE):- iaer (0-12) = 1 气溶胶类型。1代表ContinentalAEROSOL MODEL (CONCENTRATION):-v (visibility) = 15.0 能见度，需要查阅当地的气象材料xps (in km) = -0.120 当地的海拔，加上-号xpp =-705 传感器的高度，加上-号ENTER THE SPECTRAL CONDITIONS:-iwave (-1 to 34) = 42 传感器接收波段，如MODIS为42-48 GROUNRD REFLECTANCE (TYPE)-inhomo (0-1) = 0 是否为均质地物，0均一，1，不均一idirec = 0 是否考虑二向性，0表示没有，1表示有 igroun = 4 地物类型，如1表示植被，4表示湖泊rapp = 48.06356 RAPP=gias*DN+bias，大气层顶太阳辐亮度 （1、2、5波段分别计算）分别计算1、2、5波段，6S的结果中，在第3页找到计算的瑞利散射反射率值，记录下来五、计算FAI ENVI下 打开反射率数据 计算RTOA-R6S，R6S为1、2、5波段分别用6S计算出来的值。Basic Tools下找到Band Math 输入指定b1，命名为：日期band1_Rrc同样的方法计算band2_Rrc，band5_Rrc 输入以下公式，计算band2_Rrc在这里，b1指定为band1_Rrc ，b5指定为band5_Rrc，输出为band2_Rrc 计算FAI 使用b1-b2， 在这里，b1指定为band2_Rrc ，b2指定为band2_Rrc，输出为FAI五、使用FAI提取蓝藻此过程同蓝藻监测报告中的步骤，先使用模型（am pm），进展分区，再使用决策树设定阈值进展。