ENVI 5.3保姆级教程：手把手搞定Landsat 7影像从辐射定标到FLAASH大气校正的全流程
ENVI 5.3实战指南：Landsat 7影像处理全流程精解
遥感影像处理是地理信息科学中的基础技能，而ENVI作为行业标准软件，其功能强大但操作门槛较高。本文将带您从零开始，逐步完成Landsat 7影像从原始数据到最终产品的完整处理流程。不同于零散的教程，我们特别注重步骤间的衔接与常见陷阱的规避，确保您能一次性成功复现整个流程。
1. 环境准备与数据导入
在开始处理前，我们需要做好充分的准备工作。首先确保您的ENVI 5.3已正确安装并激活，建议使用64位版本以获得更好的性能体验。硬件方面，8GB以上内存和SSD硬盘会显著提升处理速度。
数据准备清单：
- Landsat 7 Level-1产品（通常包含多个波段文件和一个MTL元数据文件）
- 确保存储路径为纯英文（如
D:\RS_Data\Landsat7） - 预留至少10GB的临时空间
注意：中文路径或特殊字符可能导致ENVI无法识别文件，这是新手最常见的错误之一
数据导入步骤：
- 启动ENVI 5.3，点击菜单栏的
File > Open - 导航至数据存储位置，选择
*_MTL.txt元数据文件 - ENVI会自动识别并加载所有波段，在
Layer Manager中可查看加载结果
# 快速检查数据完整性的方法 print, '验证波段数量：', n_elements(band_list)2. 辐射定标关键操作
辐射定标是将DN值转换为真实物理量的关键步骤。Landsat 7的定标参数已内置于ENVI，但需要特别注意几个细节：
常见问题排查表：
| 问题现象 | 可能原因 | 解决方案 |
|---|---|---|
| 输出图像全黑 | 未正确选择定标类型 | 确保选择"Radiance"而非"Reflectance" |
| 数值异常大 | 波段顺序错误 | 检查MTL文件中的波段编号 |
| 图像偏色 | 未应用太阳高度角校正 | 勾选"Apply Sun Angle Correction" |
具体操作流程：
- 在工具箱中选择
Radiometric Correction > Radiometric Calibration - 选择输入文件后，设置参数：
- 定标类型：
Radiance - 输出格式：
Floating Point - 勾选太阳高度角校正
- 定标类型：
- 指定输出路径（建议新建
Calibration文件夹）
# 辐射定标后的数据验证代码示例 radiance = ENVISubsetRaster(cal_file, bands=[3,2,1]) # 选取RGB波段 histogram = HISTOGRAM(radiance, MIN=0, MAX=100)3. 波段合成与格式转换
多光谱影像处理需要将单波段文件合并为多波段文件，这个步骤直接影响后续的大气校正质量。
关键参数设置：
- 合成顺序：按波长从短到长排列（B1→B7）
- 输出格式：必须选择
BIP（Band Interleaved by Pixel） - 空间分辨率：保持30米一致
- 坐标系：检查是否与原始数据一致
操作技巧：
- 使用
Layer Stacking工具依次添加各波段 - 在高级选项中设置
Data Interleave为BIP - 输出文件命名建议包含日期和
stacked标识（如LC07_20230520_stacked.dat）
提示：格式转换耗时较长，建议在处理大批量数据时使用ENVI的批处理功能
4. FLAASH大气校正全流程
FLAASH是ENVI中最精确的大气校正模块，但其参数设置复杂，容易出错。以下是经过验证的可靠配置方案：
大气参数设置指南：
| 参数项 | Landsat 7推荐值 | 科学依据 |
|---|---|---|
| 传感器类型 | Landsat 7 ETM+ | 自动匹配波段响应函数 |
| 地面高程 | 从DEM获取或手动输入 | 影响大气厚度计算 |
| 气溶胶模型 | 中纬度夏季/冬季 | 根据影像日期选择 |
| 能见度 | 40km（晴天默认） | 可借助气象数据优化 |
| 水汽反演 | 使用波段比例法 | 适用于多光谱数据 |
逐步操作：
- 打开
FLAASH Atmospheric Correction工具 - 输入辐射定标后的堆叠文件
- 设置输出反射率文件路径
- 在
Advanced Settings中确认：- 波长单位：微米（μm）
- 重采样方法：
Nearest Neighbor - 输出缩放因子：
10000
# FLAASH后处理质量检查代码 ref_img = ENVIOpenRaster('flaash_result.dat') stats = STATISTICS(ref_img, PERCENTILE=[2,98]) print, '反射率合理范围验证：', stats5. 结果验证与常见问题
完成所有步骤后，必须对结果进行质量检查。以下是专业用户常用的验证方法：
多维度验证指标：
- 光谱曲线：植被应呈现典型"双峰"特征
- 数值范围：反射率应在0-1之间（或0-10000若设置了缩放因子）
- 空间一致性：相邻景影像接边处不应有明显色差
- 统计特征：各波段直方图应合理分布
当遇到问题时，可尝试以下排查步骤：
- 检查中间文件的元数据是否完整
- 重新确认各步骤的参数设置
- 在ENVI的
Help文档中搜索特定错误代码 - 对比官方示例数据的结果
实际项目中，我们曾发现当太阳高度角低于30度时，需要额外进行地形校正。这种情况下，建议在FLAASH后接Topographic Correction工具，使用DEM数据消除阴影影响。