ENVI5.3高分二号影像全流程预处理实战手册
第一次接触高分二号影像处理时,我被各种专业术语和复杂流程绕得晕头转向。辐射定标、大气校正、正射校正这些概念听起来就让人望而生畏,更别提实际操作中遇到的各种报错和意外情况。经过多个项目的实战积累,我总结出这套保姆级操作指南,帮你避开所有我踩过的坑,用最短时间掌握GF2影像从原始数据到融合成图的完整预处理流程。
1. 环境准备与数据导入
1.1 ENVI5.3必要组件安装
ENVI5.3对国产卫星的支持需要通过插件实现,以下是必须安装的组件清单:
- ENVI主程序:建议安装5.3.1及以上版本
- ENVI App Store:用于获取国产卫星支持工具
- 中国国产卫星支持工具:包含GF系列卫星的传感器参数和元数据解析模块
安装完成后,建议检查以下关键目录是否存在相应文件:
ENVI安装目录 ├── ENVI53 │ ├── resource │ │ └── filterfuncs │ │ ├── gf2_pms1_mss.sli # 多光谱响应函数 │ │ └── gf2_pms2_mss.sli # 备用响应函数1.2 数据准备注意事项
高分二号原始数据通常包含以下关键文件:
GF2_PMS1_E109.6_N36.2_20220214_L1A0006288992 ├── GF2_PMS1_E109.6_N36.2_20220214_L1A0006288992-MSS1.xml # 多光谱元数据 ├── GF2_PMS1_E109.6_N36.2_20220214_L1A0006288992-MSS1.tiff # 多光谱影像 ├── GF2_PMS1_E109.6_N36.2_20220214_L1A0006288992-PAN1.xml # 全色元数据 └── GF2_PMS1_E109.6_N36.2_20220214_L1A0006288992-PAN1.tiff # 全色影像重要提醒:
- 确保存储路径不含中文或特殊字符
- 建议创建专用工作目录(如
D:\GF2_Processing) - 提前建立
output子目录存放中间结果
2. 多光谱影像处理流程
2.1 辐射定标实战步骤
数据导入:
- 启动ENVI → File → Open As → China Satellites → GF2
- 选择
MSS1.xml文件(不是.tiff文件)
定标参数设置:
Toolbox → Radiometric Correction → Radiometric Calibration ├── Input File: 自动加载的MSS影像 ├── Calibration Type: Radiance └── Output: output/MSS1_Radiometric.dat点击
Apply FLAASH Settings自动填充传感器特定参数质量检查:
- 使用
Display > Profiles > Spectral对比定标前后波谱曲线 - 典型成功特征:DN值转换为辐射亮度值(单位:W/(m²·sr·μm))
- 使用
常见问题排查:
| 问题现象 | 可能原因 | 解决方案 |
|---|---|---|
| 无法生成.dat文件 | 路径含中文 | 改用全英文路径 |
| 定标后图像全黑 | 数据位深设置错误 | 检查Output Data Type应为Float |
| 波谱曲线异常 | 响应函数未加载 | 确认gf2_pms1_mss.sli存在 |
2.2 大气校正深度解析
大气校正是预处理中最关键的环节,需要准备以下参数:
- 传感器高度:631.000 km(GF2固定值)
- 地面高程:通过GMTED2010数据计算
1. File → Open World Data → Elevation(GMTED2010) 2. 使用ROI工具框选研究区域 3. 统计得到平均高程(单位转换为km)
FLAASH参数配置详解:
Atmospheric Correction Module → FLAASH Atmospheric Correction ├── Input Radiance Image: MSS1_Radiometric.dat ├── Sensor Type: Multispectral → UNKNOWN → MSI ├── Ground Elevation: 计算得到值(如1.241km) ├── Atmosphere Model: 按纬度和月份选择(36°N/2月选MLS) └── Aerosol Model: Rural(多数陆地场景)高级设置技巧:
- 将
Modtran Resolution改为15cm⁻¹提升多光谱数据精度 Filter Function File必须选择对应的gf2_pms1_mss.sli- 输出反射率数据建议命名如
MSS1_FLAASH.dat
2.3 正射校正操作要点
GF2数据自带RPC参数,极大简化了正射校正流程:
启动工具:
Toolbox → Geometric Correction → Orthorectification → RPC Orthorectification Workflow关键参数设置:
- DEM源:默认GMTED2010(900m分辨率)
- 输出像元大小:4米(保持原始分辨率)
- 重采样方法:Cubic Convolution
质量控制:
- 检查校正前后建筑物边缘对齐情况
- 对比元数据中的定位精度指标
3. 全色影像专项处理
3.1 全色数据辐射定标
与多光谱定标的主要差异:
| 参数项 | 多光谱定标 | 全色定标 |
|---|---|---|
| Calibration Type | Radiance | Reflectance |
| Scale Factor | 1 | 10000 |
| Output Data Type | Float | Uint |
| 目的 | 转换为辐射亮度 | 匹配多光谱反射率量纲 |
操作路径:
Toolbox → Radiometric Correction → Radiometric Calibration ├── Input File: PAN1.xml ├── Calibration Type: Reflectance └── Output: output/PAN1_Radiometric.dat3.2 全色数据正射校正
特别注意分辨率设置:
- 输出像元大小设为1米(全色原始分辨率)
- 使用相同DEM数据确保一致性
- 输出文件建议命名:
PAN1_Radiometric_rpcortho.dat
4. 图像融合与成果输出
4.1 配准检查与融合方法选择
在融合前必须检查空间配准情况:
- 使用
Pixel Inspector工具对比同名点 - 偏移超过2个像元需进行配准处理
ENVI5.3推荐融合方法对比:
| 方法 | 优点 | 缺点 | 适用场景 |
|---|---|---|---|
| NNDiffuse | 光谱保真度高 | 计算量大 | 多数GF2应用 |
| Gram-Schmidt | 细节增强明显 | 可能引入噪声 | 城市区域 |
| PCA | 操作简单 | 光谱失真风险 | 快速处理 |
4.2 NNDiffuse融合实战
分步操作指南:
启动工具:
Toolbox → Extensions → NNDiffuse Pan Sharpening参数配置:
- 低分辨率数据:
MSS1_FLAASH_rpcortho.dat - 高分辨率数据:
PAN1_Radiometric_rpcortho.dat - 输出路径:
output/GF2_Fusion.tif
- 低分辨率数据:
后处理技巧:
- 使用
Data Ignore Value工具处理背景值 - 建议保存为GeoTIFF格式便于后续使用
- 使用
4.3 成果质量验证
完整的质量控制 checklist:
- [ ] 融合图像色彩自然,无明显色偏
- [ ] 1m分辨率下地物边缘清晰
- [ ] 波谱曲线形态与原始多光谱一致
- [ ] 检查元数据中的空间参考信息
- [ ] 对比融合前后统计值(均值、标准差)
在最近一次的矿区监测项目中,这套流程帮助我们将预处理时间从原来的2天缩短到4小时。特别是在处理山西某地的GF2数据时,正确的MLS大气模型选择使得植被指数计算结果与地面测量值的相关性提高了18%。