遥感数据处理之MCD19A2数据的后续处理

以前有篇博客对MCD19A2数据坐了个初步的介绍，但是后续处理没有进行说明，本博客是对MCD19A2后续处理的补充。

利用MCTK处理MCD19A2影像：

由图可知，一般处理出来的数据都有多个波段，但是一般都是拿几个波段有效值合成的数据作为结果数据进行使用的（至少我一般都是这样的），以下附上IDL对MCD19A2几个波段进行有效值均值合成的部分处理代码。

    for b=0,ncols-1 do begin
      for c=0,nrows-1 do begin
        ;先创建一个二维浮点型数组
        data_effect=fltarr(nb)
        ;获取同一位置，每个波段的值。
        for d=0,nb-1 do begin
          data_effect[d]=data[b,c,d]
        endfor       
        ;进行有效值均值合成
        data_effect=data_effect[where(data_effect ge 0 and data_effect lt 50)]
        final_value[b,c]=mean(data_effect)
      endfor
    endfor

结果如下：

以下链接为完整的IDL代码：
https://pan.baidu.com/s/1xw-6YYFt3udqo_moRe4Vmg
提取码：9527

如有错误，欢迎指正，万分感谢！