基于MODTRAN高光谱快速大气校正方法.docVIP

基于MODTRAN高光谱快速大气校正方法 　　摘 要： 高光谱遥感成像过程中一般会受大气成分的影响，为了利用高光谱数据进行有效的地物识别，这些影响必须消除。详细介绍大气校正的一个过程：首先采用大气预处理微分吸收法（APDA）和暗像元法，通过并行调用MODTRAN先后建立水汽和能见度的查找表，从而计算出图像中水汽含量和能见度，并利用得出的这些参数建立反射率查找表，最终对高光谱图像进行逐像元校正。该过程结合辐射传输模型法和基于图像信息的校正方法的优点，不需要实时的大气参数也能进行绝对大气校正。同时，加入并行处理使运行速度显著提高。通过对校正后的数据进行对比分析，结果表明校正后的光谱曲线与地面实物和其他软件校正后的光谱特征相吻合，该方法能够满足大气辐射校正的基本要求。 　　关键词： 大气校正； MODTRAN； APDA； 暗像元法 　　中图分类号： TN911?34 文献标识码： A 文章编号： 1004?373X（2015）10?0023?05 　　0 引 言 　　20世纪80年代初期成像光谱（高光谱）技术的出现，使光学遥感进入了崭新的阶段。高光谱遥感能在电磁波谱的可见光和红外区域， 获取许多非常窄的光谱图像数据，为每个像元提供数十至数百个窄波段（一般小于10 nm）的光谱信息，产生一条完整而连续的光谱曲线，以便于识别出地物特征，具有波段多、光谱分辨率高、图谱合一等特点。大气校正是遥感信息定量化过程中不可缺少的环节。在高光谱图像的拍摄过程中，由于大气的存在，进入大气的太阳辐射会发生反射、折射、吸收、散射[1]，其中对传感器接受影响较大的是吸收和散射。由于大气的吸收和散射，使得进入到传感器中的信号减弱[2]。同时，大气的散射光也有一部分直接或经过背景地物反射进入到传感器时却又使得原型号增强了。为了消除这些因素的影响，必须对高光谱数据进行大气辐射校正。 　　目前大气校正的方法主要有基于图像特征的经验线性法、相对校正方法和基于物理的辐射传输模型法[3]。经验线性法需要做大量的地物反射率测量工作，通过统计的方法建立传感器接收到的辐亮度与地物反射率的线性回归方程。基于图像特征的相对校正方法得到的不是地物的真实反射率，而是相对反射率。它不需要实际地面光谱和大气参数的测量，可以应用于历史数据。它校正后的结果不是具体的反射率或辐亮度，相同的结果只表示相同的反射率，至于多少不能给定[4]。辐射传输模型法建立在大气辐射传输方程的基础上，在进行大气辐射校正时需要输入各种大气参数（例如：卫星过境时的地物反射率、大气的能见度、太阳天顶角和卫星传感器的标定参数等），获得这些参数的代价昂贵，又不易实现，尤其对于历史存档数据，获取这些参数更是不可能的[5]。陶东兴，赵慧洁通过调用MODTRAN分别建立水汽吸收通道和2个参考通道的2个查找表[6]，搜索这2个查找表使用三次样条插值方法计算水汽含量，实验结果表明生成的水汽含量分布图精度与MODTRAN计算结果只有0.1%的误差。另外还有暗像元法，暗目标方法通过确定暗目标在可见光波段的地表反射率来估计气溶胶参数。赵祥等选取小麦作为新的暗目标[4]，使用波段间的线性回归与插值的方法对常规暗目标方法进行改进，从而将暗目标方法扩展应用于高光谱遥感数据的大气校正。彭妮娜等提出了浓密植物像元660 nm与830 nm波段处的线性关系进行气溶胶含量反演[7]，为获取气溶胶含量提供了一种新途径。 　　1 概 述 　　为了保证大气校正的精确性和程序运行速度，结合前人的工作经验，总结出一套大气辐射校正方法，形成一整套大气辐射校正的具体流程，在流程中加入并行计算提高运行效率。 　　该方法具体思路是通过大气预处理微分吸收法（APDA法）并调用MODTRAN建立的水汽查找表来获得水汽含量，然后采用暗像元法并调用MODTRAN建立的能见度查找表获得能见度，最后利用水汽含量和能见度参数调用MODTRAN建立反射率查找表，计算FLAASH模型中的大气辐射校正的参数从而实现大气辐射校正的目标。 　　运用该方法开发了用于大气辐射校正的软件。并对调用MODTRAN建立查找表进行并行处理，提高运行效率。该软件包括水汽反演、能见度反演、大气辐射校正这三大模块。为了验证这三大模块的精确性，对东天山地区高光谱数据（已完成辐射定标和几何校正）进行大气辐射校正，将校正后的光谱和地面采集的光谱进行对应像元对比，分析校正结果的准确性。本文方法的大体思路流程见图1。 　　 　　图1 大气校正流程图 　　2 MODTRAN快速查找表并行处理 　　2.1 调用MODTRAN计算大气参数 　　在调用MODTRAN后会生成tape7文件，该文件得到一系列参数（FREQ，TOTAL RAD，GRND RFLT，DRCT RFLT等）。 　　利用这