请教MODIS 250M L1B数据FLAASH大气校正【已解决】

午后阳光

       请教各位大神，用MODIS 250m分辨率L1B数据计算NDVI，需要用反射率数据计算，涉及到FLAASH大气校正，即将辐射定标后的大气反射率校正为真实的地表反射率，从而消除大气影响。而FLAASH校正输入的数据要求是辐射率数据，这样的话，计算NDVI该如何进行大气校正？
（是否是我的理解有误，请求各位大神们指教！不胜感激！）

午后阳光

其实这个困惑我好久的问题，一是自己理解有误，二是没有仔细认真研读书籍。也是给了自己一个深刻的教训，万事不可眼高手低，尤其是零基础门外汉入门系列，更是要好好理解基本概念，不能自己凭一知半解的知识就去搞科研~~~心累。
       我们接收到的卫星数据（L1级数据）通常是地物不同波段的电磁波谱信息，其中的像素值称为亮度值或者灰度值，也就是DN值。DN值需要经过几何定标和辐射定标。这里主要说下辐射定标。在ENVI打开MODIS数据时会自动定标为反射率、辐射率和发射率。此处的反射率为大气表观反射，需要通过大气校正，消除水汽、气溶胶吸收和散射影响，获得真实的地表反射率，再通过波段运算计出NDVI。
       在我以前的误解中，固定思维模式，认为大气表观反射率通过大气校正后就获得真实地表反射率，所以在ENVI预处理数据时，一直尝试用表观反射率数据进行大气校正。发现FLAASH大气校正模式输入数据要求是订正后的辐射率数据，而其他大气校正模式如QUAC，虽然输入数据可以是辐射亮度值。也可以是表观反射率，但是对于波段数有要求。所以对于MODIS 250m的数据进行大气校正就遇到了瓶颈。
      可能是在遥感专业这属于基本知识或者常识，网上查阅并没有找到相关类似问题，一直到后来再次拿起书本，前后好好研读了关于大气校正这部分内容，才发现其实书本上早就说明了。FLAASH大气校正模式对于数据输入要求为订正后的辐射率数据，但是最终通过一系列校正后，输出的数据是地表真实反射率数据（除以10000）。这期间经过了怎样的订正，使得辐射率数据转为反射率数据。我又继续查阅网上相关基础理论知识。
       L=gain*DN+bias   ——根据DN值计算得到各波段辐射亮度值（辐射率）
其中，gain为增益值，bias为偏移值
       ρ=Pi * L * d *d /(ESUN*cosθ)  ——各波段表观反射率计算
其中，ρ为表观反射率或者是行星反射率，L为辐射亮度值，d为日地距离，ESUN为大气层平均入射太阳辐射，θ为太阳天顶角。
      通过上述公式转换可以看出，反射率是由辐射亮度值计算得来，只是没有进行大气校正，计算获得的是表观反射率。

午后阳光

       这个帖子之前已经编辑过，问题已解决。后来家园有几天一直登录不上，再次登陆后就发现第二次编辑的没有被保存，只好再次编辑。话说洋洋洒洒的写了好多，重写一遍，想哭。我就言简意赅吧，和大家分享一下，希望能帮到后面的童孩。