在“ENVI下的Landsat8大气校正(初试)”文章最后提出了一个讨论,这里根据这个设想做出了另外一个结果。结果分析显示,两种方法得到的结果基本一致。
第一步:重新制作波谱响应函数
(1)启动ENVIclassic,选择Window->Start New PlotWindow。
(2)在ENVI PlotWindow窗口中,选择File->Input Data–>Spectral Library,打开ldcm_oli.sli波谱响应文件。
(3)如下图所示选择7个波段的波谱响应函数,选择File->Save plot as->Spectral Library,按照默认参数保存为.sli文件。
图1:选择7个波段的波谱响应文件
第二步:大气校正
(1)使用layerstacking工具将辐射定标后的文件保存为7个波段的文件,也就是去除卷云Cirrus波段。
(2)打开FLAASH工具进行大气校正。
第三步:浏览结果
分别对两种结果进行统计,如下图所示,每个波段的均值和方差相差非常小(个位数以内),折合0~1反射率在10-3范围内,相差甚小。对比单个像素的值也是这个结果。
值得注意的是得到的反射率范围是小于0和大于10000,其实浏览直方图发现,小于0的像素只有不到10000个,占0.02%,大于10000的不到100个像素,属于正常范围内。
因此可以看到,两种方法在精度上相差不大,结果认为是一致的。
图2:统计结果(左-8波段,右-7波段)