RS—|下载Landsat8/9数据并进行ENVI大气校正(FLAASH模型）

1、 数据的下载。

下载网址：链接: GloVis (usgs.gov)。下载的数据为2022年8月1日湖南省北部的遥感影像数据。该数据为L1级产品，只经过了几何校正，没有经过辐射定标和大气校正。

图1-1.下载解压后的数据文件


图1-2.用ENVI打开数据（Linear 2%显示）

2、 辐射定标

（辐射定标是用户需要计算地物的光谱反射率或光谱辐射亮度时，或者需要对不同时间、不同传感器获取的图像进行比较时，都必须将图像的亮度灰度值转换为绝对的辐射亮度，这个过程就是辐射定标。）
1） 在ENVI的Toobox中，选择Radiometric Correction>Radiometric Calibration，选择多光谱数据（MultiSpectral多光谱数据（1~7波段）。定标类型选择辐射亮度值（Radiance），输出格式为BIL（因为后续FLAASH大气校正的输入数据类型为BIL），输出数据类型为Float，系数为0.1。

图1-3.进行辐射定标设置

图1-4.完成辐射定标

3、大气校正

（大气校正是指传感器最终测得的地面目标的总辐射亮度并不是地表真实反射率的反映，其中包含了由大气吸收，尤其是散射作用造成的辐射量误差。大气校正就是消除这些由大气影响所造成的辐射误差，反演地物真实的表面反射率的过程。）
1） 计算研究区高程，打开ENVI自带的GMTED2010.jp2全球高程数据，打开工具箱中的Statistics->Compute Statistics，在输入文件对话框中选择GMTED2010.jp2数据，再单击Stats Subset，单击File，选择需要统计高程信息对应的影像区域，取计算出的平均高程结果为0.163562459km；

图1-5.打开GMTED2010.jp2全球高程数据

图1-6.计算影像的区域高程参数设置

图1-7.高程计算结果

2） 选择Toolbox>Radiometric Correction>Atmospheric Correction Module>FLAASH Atmospheric Correction，打开FLAASH大气校正工具。
【Input Radiance Image】 文件输入就是上面辐射定标裁剪后的影像；Lat和Lon（经纬度）和时间将自动获取；传感器类型（Sensor Type）：多光谱数据（MultiSpectral）—Landsat8 OLI（因为其与Landsat9轨道高度等参数接近）

图1-8.打开FLAASH大气校正工具

图1-9.相关参数的设置

3） 大气模型参数选择 Atmospheric Model，根据成像时间和纬度信息依据【help】中的参数对应表；气溶胶模型（Aerosol Model）:Rural郊区 （由图观察；气溶胶反演方法 Aerosol Retrieval：2-band（K-T）（ 初始能见度 Initial Visibility 只有在气溶胶反演方法为 None 时候， 以及 K-T 方法在没有找到黑暗像元的情况下）；

图1-10. 数据经纬度与获取时间的大气模型参数选择对应表

图1-11.相关参数的设置

4） 多光谱数据参数设置 Multispectral Settings，使用的是默认的第一个

图1-12. 多光谱数据参数设置

图1-13.大气校正参数设置完毕

5） 点击Apply，等待处理结果。使用相关的拉伸方式显示并对比大气校正前后的图像，可以看到校正后的图像亮度有所增强。同时对比校正前后同一像素点自然植被的波谱曲线，上方为校正后波谱曲线，下方为校正前，可以看到校正后的波谱曲线趋势符合真实自然植被的波谱，即经过FLAASH校正的影像基本去除了空气中水汽颗粒等因子的影响，植被的波谱曲线趋于正常。

图1-14.大气校正结果处理中

图1-15.大气校正前后图像对比（上图为校正后）

图1-16.对比校正前后波谱曲线

参考文章：
1、ENVI下Landsat8辐射定标与大气校正

2、ENVI软件对Landsat-8数据进行辐射定标、大气校正、提取NDVI、估算植被覆盖度等操作
3、https://blog.csdn.net/betterldd93/article/details/123852327