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

江北20190411

关注

原创最新推荐文章于 2026-07-21 09:28:03 发布 1.8k 阅读 AI 写同款· GEO 优化 ›

实习一：6S大气校正 ********************************************************************************************* 以下内容只用于记录平时上课实习的实习报告，以便将来查阅。 ********************************************************************************************* 一、实习目的： 大气校正的目的是消除或减少大气对图像的干扰。本次实习要阅读详情

前言

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