Landsat8移除卷云Cirrus波段的大气校正测试

在“ENVI下的Landsat8大气校正（初试）”文章最后提出了一个讨论，这里根据这个设想做出了另外一个结果。结果分析显示，两种方法得到的结果基本一致。
第一步：重新制作波谱响应函数
（1） 启动ENVI classic，选择Window->Start New Plot Window。
（2） 在ENVI Plot Window窗口中，选择File->Input Data –>Spectral Library，打开ldcm_oli.sli波谱响应文件。
（3） 如下图所示选择7个波段的波谱响应函数，选择File->Save plot as-> Spectral Library，按照默认参数保存为.sli文件。

图1：选择7个波段的波谱响应文件
第二步：大气校正
（1） 使用layer stacking工具将辐射定标后的文件保存为7个波段的文件，也就是去除卷云Cirrus波段。
（2） 打开FLAASH工具进行大气校正。
第三步：浏览结果
分别对两种结果进行统计，如下图所示，每个波段的均值和方差相差非常小（个位数以内），折合0~1反射率在10-3范围内，相差甚小。对比单个像素的值也是这个结果。
值得注意的是得到的反射率范围是小于0和大于10000，其实浏览直方图发现，小于0的像素只有不到10000个，占0.02%，大于10000的不到100个像素，属于正常范围内。
因此可以看到，两种方法在精度上相差不大，结果认为是一致的。


图2：统计结果（左-8波段，右-7波段）
文章来源：http://blog.sina.com.cn/s/blog_764b1e9d01018ud7.html