高分一号PMS数据预处理实战：从L1A级数据到融合影像的ENVI全流程解析

1. 高分一号PMS数据预处理入门指南第一次接触高分一号PMS数据时我也被那一堆专业术语搞得晕头转向。L1A级数据、全色多光谱、辐射定标...这些词听着就让人头大。但别担心今天我就用最直白的语言带你一步步搞定这个看似复杂的预处理流程。高分一号是我国重要的对地观测卫星它的PMS传感器能同时获取2米分辨率的全色影像和8米分辨率的多光谱影像。简单来说全色影像就像黑白照片细节丰富多光谱影像则像彩色照片包含多个波段的光谱信息。L1A级是卫星下传的原始数据级别相当于相机的RAW格式需要经过一系列处理才能用于分析。ENVI作为遥感领域的瑞士军刀是我们处理这些数据的得力助手。最新版的ENVI5.6已经内置了对国产卫星的良好支持这让我们省去了很多麻烦。整个预处理流程可以概括为数据加载→辐射定标→大气校正→正射校正→影像融合。听起来步骤不少但跟着我的操作走保证你能轻松掌握。2. 数据准备与环境配置2.1 安装必备工具首先确保你的ENVI版本是5.6或更新。打开ENVI后第一件事就是安装国产卫星支持工具。这个工具就像是个翻译官能让ENVI正确读取高分一号的数据格式。具体操作点击ENVI顶部菜单的Extensions→Extension Manager在商店里搜索China Satellite Support安装后重启ENVI。这个步骤我当初花了半小时才搞明白现在你只需要30秒就能完成。2.2 加载原始数据安装好支持工具后就可以加载数据了。高分一号的PMS数据通常会包含两个文件一个是8米分辨率的多光谱影像MSS另一个是2米分辨率的全色影像PAN。在ENVI中点击File→Open找到你的数据存放位置。这里有个小技巧如果数据是压缩包格式建议先解压到单独的文件夹避免路径过长导致读取失败。我第一次处理时就因为路径问题卡了半天后来发现把文件夹改名为简单的GF1就解决了。3. 辐射定标实战操作3.1 多光谱影像定标辐射定标就是把原始的DN值数字量化值转换为更有物理意义的辐射亮度值。在ENVI中这个操作非常简单打开Toolbox→Radiometric Correction→Radiometric Calibration选择你的多光谱影像在参数设置界面点击Apply FLAASH Settings这个按钮简直是懒人福音设置输出路径和文件名建议命名为MS_RadCal方便后续识别这里有个重要细节定标系数通常会自动读取但如果遇到错误提示可能需要手动输入。高分一号的定标系数可以在卫星数据说明文档中找到一般格式为Gain和Bias两个参数。3.2 全色影像特殊处理全色影像的定标稍有不同因为它后续要参与影像融合。我们需要把它转换为与多光谱数据匹配的表观反射率同样使用辐射定标工具选择全色影像在输出类型中选择Reflectance关键步骤勾选Scale Factor并设置为10000输出命名为PAN_RadCal这个缩放操作是为了让全色影像的数值范围与经过大气校正的多光谱影像保持一致。我第一次做的时候忽略了这点结果融合后的影像色彩严重失真不得不返工重做。4. 大气校正关键步骤4.1 FLAASH参数设置大气校正是预处理中最复杂的环节之一但ENVI的FLAASH模块让它变得相对简单打开Toolbox→Radiometric Correction→FLAASH Atmospheric Correction输入刚才辐射定标后的多光谱影像海拔设置点击File→Open World Data→Elevation(GMTED2010)使用Compute Band Statistics工具获取影像区域的平均海拔这里有个实用技巧如果处理区域地形起伏不大直接输入该地区的平均海拔即可。我曾经为了精确试图用每个像元的海拔结果不仅耗时翻倍最终效果提升却微乎其微。4.2 高级参数调整在FLAASH的高级设置中有几个关键参数需要注意大气模型根据影像拍摄时间和地理位置选择。中国大部分地区适合选择Mid-Latitude Summer或Mid-Latitude Winter气溶胶模型城市区域选Urban乡村选Rural内存设置8G以下电脑建议开启分块计算Tile Size设为100-200MB特别提醒全色影像不需要进行大气校正这是很多新手容易犯的错误。全色波段单一无法完成大气校正所需的计算强行处理只会得到错误结果。5. 几何精校正技巧5.1 多光谱正射校正高分一号的影像已经带有RPC有理多项式系数文件我们可以利用这个信息进行正射校正打开Toolbox→Geometric Correction→RPC Orthorectification Workflow选择大气校正后的多光谱影像DEM选择优先使用影像自带的DEM如果没有再用GMTED2010在Advanced面板中设置输出像元大小为8米输出命名为MS_Ortho这里有个常见问题校正后的影像边缘可能出现变形。如果研究区域不在边缘可以直接忽略若在边缘可以考虑适当扩大处理范围。5.2 全色影像校正全色影像的校正流程与多光谱类似但有两点不同输入的是辐射定标后的全色影像输出像元大小设为2米命名为PAN_Ortho特别注意全色和多光谱影像必须使用相同的DEM和校正参数否则后续融合时会出现配准问题。我曾经因为DEM不一致导致融合影像出现重影浪费了一整天时间排查问题。6. 影像融合终极操作6.1 数据格式转换ENVI中影像存储顺序有BSQ、BIL、BIP三种。为了提高融合速度建议将多光谱影像转为BIP格式打开Toolbox→Raster Management→Convert Interleave选择正射校正后的多光谱影像输出格式选BIP命名为MS_Ortho_BIP6.2 黑边处理技巧高分一号影像常有黑色边框会影响融合效果。处理方法是双击图层打开属性窗口点击Edit Metadata添加Data Ignore Value并设为0保存修改这个步骤看似简单却能显著提升融合质量。我对比过处理前后的结果忽略黑边后的融合影像边缘更加自然。6.3 Gram-Schmidt融合终于到了最激动人心的融合环节打开Toolbox→Image Sharpening→Gram-Schmidt Pan Sharpening先选择低分辨率的多光谱影像MS_Ortho_BIP再选择高分辨率的全色影像PAN_Ortho设置输出路径和文件名建议Fusion_ResultGram-Schmidt方法能很好地保持光谱特性是高分一号数据融合的首选算法。融合完成后你会得到兼具2米空间分辨率和多光谱特性的完美影像。第一次看到自己亲手处理的融合结果时那种成就感至今难忘。