6S模型逐像元大气校正MERSI2

前言

      FY3D-MERSI2 6S大气校正
      IDL调用6S模型
     之前完成了6S模型源码中添加MERSI2的光谱响应函数和实现了IDL调用6S模型，然后进行了应用单一参数进行大气校正。
这次实现利用IDL调用6S模型建立的大气校正系数查找表然后对MERSI2进行逐像元大气校正。
     由于MERSI2自带的角度信息是逐像元的数据，而我建立的查找表步长比较大，因此在进行校正前，先对角度信息和MODIS AOD数据进行了处理。实现逐像元角度数据对应查找表中设置的数据，比如20-25之间用20来代替，相当于给原始数据分块(这里这么做,不知道合不合适，但我只能想到这种方案了)。

大气校正结果

单一参数校正时气溶胶浓度输入的是AOD均值，角度信息是中心像元的信息。
取中心像元附近行剖面简单对比下,是有差异的，深入研究还需要定量分析。

对应查找表查找系数计算部分代码。
一块1300*1300的区域校正用了240s左右。

;  y = a*b3 - b
;  acr = y / (1+c*y)
;  lut[3 4 5]分别为a.b.c
  m = size(lut, /dimensions)
  m = m[0] 
  for i = 0,m do begin
      for j=0,n_elements(b4)-1 do begin
        if vz[j] eq lut[0,i] and sz[j] eq lut[1,i] and aod[j] eq lut[2,i] then begin
            b4[j] = (lut[3,i]*b4[j]/100d - lut[4,i]) / $
            (1 + lut[5,i]*(lut[3,i]*b4[j]/100d - lut[4,i]))
        endif
      endfor
  endfor