从原始数据到精准分析：ENVI5.3驱动下的高分二号影像全流程预处理实战

1. 高分二号数据与ENVI5.3环境准备

高分二号（GF-2）作为我国自主研制的亚米级高分辨率遥感卫星，其多光谱（4米）与全色（1米）数据组合在农业监测、城市规划等领域应用广泛。但在使用前，必须完成三个关键动作：

1. 安装ENVI5.3主程序：建议选择默认安装路径（如C:\Program Files\Harris\ENVI53），避免后续插件识别问题。实测发现安装在中文路径下会导致国产卫星支持工具加载异常。

2. 配置国产卫星插件：安装完成后需额外操作：

   bash复制# 示例插件安装路径（需替换实际路径）
   cd "C:\Program Files\Harris\ENVI53\extensions"
   git clone https://github.com/ChinaSatelliteTools/GF2_Support.git
   

   安装后重启ENVI，在File > Open As菜单中会出现China Satellites选项，这是处理高分数据的关键入口。

3. 数据目录规范：新建项目文件夹时，建议采用纯英文路径结构。例如：

   code复制D:\GF2_Processing
   ├── raw_data      # 存放原始.xml和.tiff文件
   ├── output        # 处理结果输出目录
   └── temp          # 临时文件
   

   我曾遇到因路径包含中文导致辐射定标生成的.dat文件丢失的情况，这点要特别注意。

提示：首次打开GF2数据时，建议通过View Metadata检查RPC参数是否完整。若缺失需联系数据提供商补充，否则无法进行后续正射校正。

2. 多光谱影像的辐射定标实战

辐射定标是将原始DN值转换为具有物理意义的辐射亮度值的关键步骤。在ENVI5.3中操作时，有几个易错点需要特别注意：

操作流程分解：

1. 通过File > Open As > China Satellites > GF2加载MSS1.xml文件
2. 在Toolbox中选择Radiometric Correction > Radiometric Calibration
3. 参数设置技巧：
   • Calibration Type选择Radiance（必须）
   • 勾选Apply FLAASH Settings自动填充光谱参数
   • 输出文件命名建议添加_Rad后缀（如MSS1_Rad.dat）

避坑指南：

• 当遇到"Failed to generate output"报错时，90%的情况是输出路径包含中文或特殊符号
• 定标后的波谱曲线验证：在Display > Profiles > Spectral中对比原始DN值与辐射亮度曲线，正常情况后者数值范围应在0-100之间

实测案例：处理一景陕西榆林地区的GF2数据时，未勾选Apply FLAASH Settings导致后续大气校正失败。正确的辐射亮度值范围应是：

python复制# 典型GF2辐射亮度值范围（单位：W/(m²·sr·μm)）
band1 = 0-85.32 
band2 = 0-74.23
band3 = 0-59.34
band4 = 0-46.25

3. 大气校正的精细化操作

大气校正环节对最终反射率精度影响极大。GF2数据需重点关注以下参数：

3.1 高程数据准备

使用GMTED2010全球DEM计算区域平均高程时，有个效率技巧：

1. 用ROI工具绘制区域时，先在全幅影像上标记四个角点
2. 在Statistics面板中查看Mean值，直接除以1000得到千米单位值
3. 复杂地形建议分区块计算后取加权平均

3.2 FLAASH参数详解

在Atmospheric Correction Module中，关键参数设置逻辑：

markdown复制| 参数项               | GF2典型值          | 选择依据                     |
|----------------------|-------------------|----------------------------|
| Sensor Type          | MSI              | 国产多光谱传感器             |
| Sensor Altitude      | 631 km           | GF2轨道高度固定值            |
| Ground Elevation     | 1.2 km           | 需实际计算                  |
| Atmospheric Model    | MLS              | 中纬度夏季(7月拍摄)          |
| Aerosol Model        | Rural            | 非城市区域                  |
| Aerosol Retrieval    | None             | 无实测能见度数据时选择       |

光谱响应函数配置：
找到ENVI安装目录下的filterfuncs文件夹，GF2-PMS1应选择gf2_pms1_mss.sli文件。常见错误是选错PMS1/PMS2版本，会导致反射率曲线异常。

校正效果验证：同一地物在校正前后波谱曲线应保持形状一致，仅数值范围变化。建议选择水体或裸土等稳定地物进行验证。

4. 正射校正与影像融合

4.1 多光谱数据正射校正

使用RPC Orthorectification Workflow工具时：

1. DEM源选择：优先使用30米分辨率DEM（如ASTER GDEM）
2. 输出像元大小严格设置为4米（保持原始分辨率）
3. 重采样方法选择Cubic Convolution能在清晰度和平滑度间取得平衡

4.2 全色数据处理差异点

• 辐射定标时选择Reflectance类型
• Scale Factor设为10000（与大气校正结果匹配）
• 正射校正输出分辨率设为1米

4.3 NNDiffuse融合技巧

在Extensions > NNDiffuse Pan Sharpening中：

1. 输入数据必须满足分辨率整数倍关系（4米多光谱+1米全色）
2. 融合后图像发白的处理方法：

   python复制# 通过头文件修改背景值
   with open('NNDiffusePanSharpening.hdr', 'a') as f:
       f.write('data ignore value = 0\n')
   

3. 质量评估：检查道路边缘、农田边界等高频特征是否保持清晰

最终成果建议保存为GeoTIFF格式，便于在GIS软件中使用。一套完整的GF2处理流程通常需要2-3小时（16GB内存配置），其中大气校正步骤耗时占比超过50%。对于批量处理，可以编写ENVI IDL脚本自动化流程。