ENVI5.3保姆级教程:高分二号影像从辐射定标到图像融合,一篇搞定所有预处理
2026/5/5 13:53:30 网站建设 项目流程

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 辐射定标实战步骤

  1. 数据导入

    • 启动ENVI → File → Open As → China Satellites → GF2
    • 选择MSS1.xml文件(不是.tiff文件)
  2. 定标参数设置

    Toolbox → Radiometric Correction → Radiometric Calibration ├── Input File: 自动加载的MSS影像 ├── Calibration Type: Radiance └── Output: output/MSS1_Radiometric.dat

    点击Apply FLAASH Settings自动填充传感器特定参数

  3. 质量检查

    • 使用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参数,极大简化了正射校正流程:

  1. 启动工具:

    Toolbox → Geometric Correction → Orthorectification → RPC Orthorectification Workflow
  2. 关键参数设置:

    • DEM源:默认GMTED2010(900m分辨率)
    • 输出像元大小:4米(保持原始分辨率)
    • 重采样方法:Cubic Convolution
  3. 质量控制:

    • 检查校正前后建筑物边缘对齐情况
    • 对比元数据中的定位精度指标

3. 全色影像专项处理

3.1 全色数据辐射定标

与多光谱定标的主要差异:

参数项多光谱定标全色定标
Calibration TypeRadianceReflectance
Scale Factor110000
Output Data TypeFloatUint
目的转换为辐射亮度匹配多光谱反射率量纲

操作路径:

Toolbox → Radiometric Correction → Radiometric Calibration ├── Input File: PAN1.xml ├── Calibration Type: Reflectance └── Output: output/PAN1_Radiometric.dat

3.2 全色数据正射校正

特别注意分辨率设置:

  • 输出像元大小设为1米(全色原始分辨率)
  • 使用相同DEM数据确保一致性
  • 输出文件建议命名:PAN1_Radiometric_rpcortho.dat

4. 图像融合与成果输出

4.1 配准检查与融合方法选择

在融合前必须检查空间配准情况:

  • 使用Pixel Inspector工具对比同名点
  • 偏移超过2个像元需进行配准处理

ENVI5.3推荐融合方法对比:

方法优点缺点适用场景
NNDiffuse光谱保真度高计算量大多数GF2应用
Gram-Schmidt细节增强明显可能引入噪声城市区域
PCA操作简单光谱失真风险快速处理

4.2 NNDiffuse融合实战

分步操作指南:

  1. 启动工具:

    Toolbox → Extensions → NNDiffuse Pan Sharpening
  2. 参数配置:

    • 低分辨率数据:MSS1_FLAASH_rpcortho.dat
    • 高分辨率数据:PAN1_Radiometric_rpcortho.dat
    • 输出路径:output/GF2_Fusion.tif
  3. 后处理技巧:

    • 使用Data Ignore Value工具处理背景值
    • 建议保存为GeoTIFF格式便于后续使用

4.3 成果质量验证

完整的质量控制 checklist:

  • [ ] 融合图像色彩自然,无明显色偏
  • [ ] 1m分辨率下地物边缘清晰
  • [ ] 波谱曲线形态与原始多光谱一致
  • [ ] 检查元数据中的空间参考信息
  • [ ] 对比融合前后统计值(均值、标准差)

在最近一次的矿区监测项目中,这套流程帮助我们将预处理时间从原来的2天缩短到4小时。特别是在处理山西某地的GF2数据时,正确的MLS大气模型选择使得植被指数计算结果与地面测量值的相关性提高了18%。

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询