ENVI下的大气校正.docx

1、ENVI下的大气校正大气校正是定量遥感中重要的组成部分。本专题包括以下内容：1.大气校正概述；2.ENVI 中的大气校正功能；1.大气校正概述大气校正的目的是消除大气和光照等因素对地物反射的影响，广义上讲获 得地物反射率、辐射率或者地表温度等真实物理模型参数；狭义上是获取地物真 实反射率数据。用来消除大气中水蒸气、氧气、二氧化碳、甲烷和臭氧等物质对 地物反射的影响，消除大气分子和气溶胶散射的影响。 大多数情况下，大气校正 同时也是反演地物真实反射率的过程。图1大气层对成像的影响示意图很多人会有疑问，什么情况下需要做大气校正，我们购买或者其他途径获 取的影像是否做过大气校正。通俗来讲，如果我们需

2、要定量反演或者获取地球信息、精确识别地物等， 需要使用影像上真实反映对太阳光的辐射情况，那么就需要做大气校正。我们购 买的影像，说明文档中会注明是经过辐射校正的， 其实这个辐射校正指的是粗的 辐射校正，只是做了系统大气校正，就跟系统几何校正的意义是一样的。目前，遥感图像的大气校正方法很多。这些校正方法按照校正后的结果可 以分为2种：绝对大气校正方法：将遥感图像的DN(Digital Number值转换为地表反射 率、地表辐射率、地表温度等的方法。相对大气校正方法：校正后得到的图像，相同的DN值表示相同的地物反射 率，其结果不考虑地物的实际反射率。常见的绝对大气校正方法有：基于辐射传输模型MOR

3、TRAN 模型LOWTRAN 模型ATCOR 模型6S 模型等基于简化辐射传输模型的黑暗像元法基于统计学模型的反射率反演；相对大气校正常见的是：基于统计的不变目标法 直方图匹配法等。既然有怎么多的方法，那么又存在方法选择问题。这里有一个总结供参考：1、如果是精细定量研究，那么选择基于基于辐射传输模型的大气校正方法。2、如果是做动态监测，那么可选择相对大气校正或者较简单的方法。3、如果参数缺少，没办法了只能选择较简单的方法了。2.ENVI大气校正功能在ENVI中包含了很多大气校正模型，包括基于辐射传输模型的 MORTRAN模型、黑暗像元法、基于统计学模型的反射率反演。 基于统计的不变目标法可以

4、利用ENVI一些功能实现。其中MORTRA模型集成在ENVI大气校正扩展模块中。 还有直方图匹配等。2.1 简化黑暗像元法大气校正黑暗像元法是一种古老、简单的经典大气校正方法。它的基本原理是在假 设待校正的遥感图像上存在黑暗像元、地表朗伯面反射和大气性质均一，并忽略大气多次散射辐照作用和邻近像元漫反射作用的前提下，反射率很小（近似 0）的黑暗像元由于大气的影响，使得这些像元的反射率相对增加，可以认为这部分 增加的反射率是由于大气影响产生的。 这样，将其他像元减去这些黑暗像元的像元值，就能减少大气（主要是大气散射）对整幅影像的影响，达到大气校正的目 的。整个过程的关键是寻找黑暗像元以及黑暗像元增

5、加的像元值。ENVI下的Dark Subtract工具提供选择波段最小值、ROI的平均值、自定义值三种方式确定黑暗像元的像素值。操作过程如下：（1） 打开待校正图像文件。（2） 在主菜单中，选择 Basic Tools-Preprocessing-General PurposeUtilities-Dark Subtract ，在文件选择对话框中选择待校正图像文件，单击OK按钮，打开 Dark Subtraction Parameters 面板。（3） 在Dark Subtraction Parameters面板中，确定黑暗像素值包括三种方法（Subtraction Method ）:波段最小值

6、（Band Minimum）ROI 的平均值（Region Of Interest）自定义值（User Value）（4） 在Output Result to 中选择File以及相应的输出路径和文件名，单击 OK执行操作。图 2 Dark Subtraction Parameters 面板(Band Minimum)2.2 基于统计学模型的反射率反演基于统计学模型的反射率反演的方法主要有平场域法 (Flat Field , FF)、对数残差法(Log Residuals )、内部平均法(Internal Average Relative Reflectanee ， IARR、经验线性法(Emp

7、irical Line) 。集中在 BasicTools-Preprocessi ng- Calibrati on Utilities 菜单下。1.平场域法(Flat Field)Flat Field 定标工具通过选择图像中一块具有高反射率、光谱变化平坦的区域，利用这个区域的平均光谱值来模拟飞行时的大气条件下的太阳光谱。 将每个像元的DN值除以选择区域的平均光谱值得到相对反射率，以此来消除大气的 影响。在使用这个工具前，需要利用 ENVI提供的感兴趣区绘制工具(ROI Tool) 在被定标图像上选择感兴趣区作为平场域 (Flat Field )，感兴趣区可选择沙漠、大块水泥地、沙地等区域。2.

8、对数残差(Log Residuals )对数残差定标工具将数据除以波段几何均值， 后再除以像元几何均值，可以消除光照、大气传输、仪器系统误差、地形影响和星体反照率对数据辐射的影响。 定标结果的值在1附近。3.内部平均法(Internal Average Relative Reflectanee ， IARR)IAR (In ternal Average Relative) Reflecta nee 定标工具假定整幅图像的平均光谱基本代表了大气影响下的太阳光谱信息。把图像 DN值与整幅图像的平均辐射光谱值相除，得到的结果为相对反射率。该工具特别适用于没有植被的干 旱区域。4.经验线性法(Empi

9、rical Line)Empirical Line 定标方法是假设图像DN值与反射率之间存在线性关系：反射率=增益* DN值+偏移利用两个已知点的地面反射光谱值，再计算图像上对应像元点的平均 DN值，然后利用线性回归求出增益和偏移值，建立 DN值与反射率之间的相互关系 式，进行反射率的定标。消除了太阳辐亮度和大气程辐射。ENVI的Empirical Line 定标工具要求至少需要一个已知区域的地面反射 光谱值(Field Spectra )作为参照波谱，以及图像上对应像元点的波谱曲线(Data Spectra )。它们可以来自波谱剖面或波谱曲线、波谱库、感兴趣区、统计文件和ASCII文件。输入

10、的波谱将自动被重采样，以与选择的数据波长相匹配。也可 以用已经存在的系数对数据集进行定标。2.3不变目标法相对大气校正相对大气校正按照数学基础可以分为 2种，非线性校正法和线性校正法。 非线性校正法最典型的是直方图匹配，图像的直方图是图像中所有灰度值的概率 分布。即将校正图像的直方图与参考图像的直方图进行匹配，使两幅图像具有相 同或相近的灰度值概率分布，达到两幅图像上同名地物具有相同灰度值的目的。利用 ENVI 中 Histogram Matchi ng 工具（Display 中，En ha nce-Histogram Matchi ng）线性校正法有个前提假设：不同时相的图像灰度值之间满足线

11、性关系，这 种假设在近似情况下是成立的。这样就可以通过线性等式来描述不同时相间的灰 度关系，用x表示参考图像，y表示待校正图像，他们之间的线性关系可描述为：y=ax+b （ 式 1）其中：a、b为线性等式中的参数，即为增益和偏移量。根据前述原理，完成线性相对校正需要以下 3个步骤：第一步，在两幅图像中搜寻相对固定目标即光谱稳定的地物样本点，即伪 不变特征要素（PIF： Pseudo-Invariant Features ）；第二步，运用这些伪不变特征点的DN值,利用线性回归的方法求解式（13.6） 中的参数，得到图像间的线性关系；第三步，根据该关系式，通过波段运算，得到与参考图像具有相同或相近

12、 辐射值的结果图像，完成相对大气校正。整个过程的关键是 PIF的选择。下面以两个不同大气环境下成像、 已经经过精确配准、ENVI标准格式的Landsat TM5数据为例（2000年和2001年），介绍线性校正法的操作步骤：第一步：PIF选择选择一幅目视质量较好的图像作为基准图像（2001年），另外一幅作为待 校正图像（2000年）。在两个图像上选择相同区域的沥青房顶、砾石面、混凝 土停机坪、洁净水体、混凝土、沙地等地物作为 PIF，这些地物不会随时间的变化而变化。（1）在主模块中，选择File-Open Image File ，打开两幅图像，并在Display中显示。（2） 在其中一幅影像上点

13、击右键，从快捷菜单中选择 Geographic Link , 将显示的两幅影像地理链接。（3） 在显示2000年图像的主图像窗口中，选择 Overlay-Region of Interest ，打开 ROI Tool 面板。（4） 通过目视方式，从两幅图像找到光谱稳定、相同地物作为样本，用 Polygon或者Point类型绘制感兴趣区。（5） 在绘制一定数量感兴趣后（不宜太多，太多后面的回归运算量会很大）， 在 ROI Tool 面板中，选择 File-Output ROIs to ASCII 。（6） 回到 ROI Tool 面板中，选择 Options- Reconcile ROIs vi

14、a Map ， 将前面绘制的ROI转接到2001年的图像上，类似（5）步的方法将基准图像的 ROI内对应像素位置和像元值输出为文本文件。分别用记事本打开上面步骤得到两个文本文件，这样我们得到了 2000年的 图像和2001年的图像相对应伪不变特征要素 （PIF）的像素值。从文本文件中可 以看到，两个时相图像中每一个波段的像素值是一一对应关系，刚好对应式 1中的x和y。第二步中就是利用这些像素值，根据最小二乘回归分析法获得式13.6 中的a和b两个参数。第二步、线性关系式求解使用最小二乘回归的方法来求解线性回归式 a和b参数，如表1。表1回归解算的a和b值波段增益（a）偏移（b）Band11.0

15、2-34Band21.22-18Band30.92-9Band41.21-16Band50.994Band70.943第三步、线性变换禾用表1中的a和b值，在ENVI的BandMath工具对待校正图像做线性变 换，然后利用Layer Stacking工具将线性变换结果组合成一个多波段文件。不变目标法相对大气校正操作过程已经完成。2.4热红外大气校正ENVI 提供Thermal Atm Correction 工具，可以近似去除热红外辐射数据中的大气影响。在进行大气校正之前，为了得到最好的结果，必须将热红外数据 定标为比辐射率数据(TIMS的热红外数据必须被转化为辐射亮度数据)，并且 待校正数据波

16、长在8-14即之间。下面以ASTER L1A的热红外波段为例，操作过程如下：(1) 打开ASTER在波段列表中按照波长自动归为 4组，2组可见光-近红 外(VNIR)、短波红外(SWIR和热红外(IIR)，并根据头文件信息自动定标为 辐射亮度值(单位 W/ (m2 *卩m*sr )。(2)在主菜单中，从以下列表中选择一种方式Basic Tools - Preprocessing -Calibration Utilities -Thermal AtmCorrectio nBasic Tools -Preprocessing -Data-Specific Utilities -Thermal IR

17、 - Thermal Atm Correcti onBasic Tools - Preprocess ing - Data-Specific Utilities -TIMS-Thermal Atm Correcti on在Thermal Correction In put File 对话框中，选择热红外数据(Wavele ngth : 8.291 to 11.318 )。(3)在 Thermal Atm Correction Parameters 面板中(图 3)，需要填写 以下参数：数据缩放系数(Data Scale Factor ) : 1。将输入数据的单位缩放为 W/(m2 * 卩 m*

18、sr)。波长单位/ Wavele ngth Un its ) : Micrometers。设定表面温度估算衰退像元/ Regression Pixels ) : All选择“ All ”，将使用整个输入波长范围内的亮度温度最大值对每个像元的 表面温度进行估算。选择“ MaxHit ”，仅对那些在特定波长具有最大亮度温度值的像元进行表 面温度估算，所说的特定波长是指包含最多具有最大亮度温度的像元的波长范围。散点图拟合技术/ Fitti ng Tech ni que ) : Top of Bi ns选择“ Top of Bins ”，将会使曲线向辐射率与亮度温度的散点图的上部拟 合，散点图的上部对

19、应着发射率接近1的像元，使用该技术的拟合线是通过对散点图上部5%勺数据做了标准的最小平方回归得到的。（注：该技术易受发生在 散点图顶部区域的传感器噪声的影响）选择“ Normalized Regression ”，将会先使用标准的最小平方回归把曲线 向辐射率与亮度温度的散点图拟合，然后将拟合线的残差与正态概率分布图相比 较，在正态图中对残差再进行另一个回归， NESR（ noise equivale nt sen sorresponse）距离的点将被认为是奇异点而被删除，最后利用减少的像元集在散点 图上进行一个最终的回归。（注：该方法使用散点图中除奇异点以外的所有点， 并且不仅仅把曲线向散点图

20、的上部（发射率近似为 1）拟合）设置输出增益与偏移参数为文件（Output Ga in/Offset File ）:可选。设置是否绘制大气透射和上行辐射光谱结果（PlotTransmission/Upwelling?） : Yes。（4）选择结果输出路径及文件名，单击 0K执行校正过程。2.5QUAC快速大气校正工具快速大气校正工具(简称QUAC自动从图像上收集不同物质的波谱信息， 获取经验值完成高光谱和多光谱的快速大气校正(图 4)。它得到结果的精度近似FLAASH或者其他基于辐射传输模型的+/-15%。目前它支持的多光谱和高光谱波谱范围是(0.4-2.5卩m)。传感器包括AISA,ASAS

21、,AVIRIS, CAFARCHERCOMPASSI YCASHYDICE,HyMap,Hyperio n, IKONOS, Landsat TM, LASH,MASTERMODIS,MTI, QuickBird, RGB,以及 unknown sensor。图4快速大气校正流程图QUAC的输入数据可以是辐射亮度值、表观反射率、无单位的 raw数据。可以是任何数据储存顺序(BIL/BIP/BSQ)和储存类型，多光谱和高光谱传感器数 据的每个波段必须有中心波长信息。QUAC的操作非常简单，如下：(1)在ENVI主菜单中，选择以下方式启动Basic Tools-Preprocess in g- C

22、alibrati on Utilities-QUick Atmospheric Correcti onSpectral-QUick Atmospheric CorrectionS pectral- Preprocess in g- Calibrati on Utilities- QUick Atmospheric Correcti on在文件输入对话框中选择校正的图像文件(2) 打开 QUick Atmospheric Correction Parameters 面板(图 5),在Sensor Type中选择相应的传感器类型，选择文件名和路径输出。图 5 QUick Atmospheric C

23、orrection Parameters 面板2.6FLAASH大气校正FLAASH 是基于MODTRAN辐射传输模型，MODTRA模型是由进行大气校正 算法研究的领先者 Spectral Sciences, Inc 和美国空军实验室(Air ForceResearch Laboratory )共同研发。ITT VIS公司负责集成和 GUI设计。1.FLAASH 特点支持传感器种类多，包括多光谱的ASTERAVHRRGeoEye-1, IKONOSIRS, Landsat, MODISSeaWiFSSPOTQuickBird ,RapidEye 等，高光谱 HyMAPAVIRIS、 CASI、

24、HYDICE HYPERION(EO-1) AISA等。可以通过自定义波谱响应函数支持 更多的传感器。工程化应用价值比较明显。FLAASH采用了 MODTRAN辐射传输模型，该算法精度高。任何有关影像的 标准MODTRAN气模型和气溶胶类型都可以直接使用。通过影像像素光谱上的特征来估计大气的属性，不依赖遥感成像时同步测 量的大气参数数据可以有效地去除水蒸气/气溶胶散射效应，同时基于像素级的校正，矫正目 标像元和邻近像元交叉辐射的“邻近效应”。对由于人为抑止而导致波谱噪声进行光谱平滑处理。作为结果，除了真实 地表反射率外，还可以得到整幅图像内的能见度、卷云与薄云的分类影像、水气 含量数据。2.使

25、用ENVI大气校正模块ENVI大气校正模块的使用主要又以下 7个方面组成：1、输入文件准备，2、 基本参数设置，3、多光谱数据参数设置，4、高光谱数据参数设置，5、高级设 置，6输出文件，7、处理结果。下面介绍这7个方面内容。（一） 输入文件准备1） 支持传感器类型高光谱包括：HyMAPAVIRIS、CASk HYDICE HYPERION（EO-l） AISA等； 多光谱包括：ASTERAVHRRGeoEye-1, IKONOSIRS，Landsat，MODISSeaWiFS SPOT QuickBird，RapidEye 等，航空：860nm-1135nm波长范围2） 数据是经过定标后的辐

26、射亮度（辐射率）数据，单位是：（卩W/（cm2*nm*sr）3） 数据带有中心波长（wavelenth /值，如果是高光谱还必须有波段宽度（FWHM,这两个参数都可以通过编辑头文件信息输入（Edit Header ）。4） 数据类型支持四种数据类型：浮点型（floating ）、长整型（long integer ）、整型（integer / 和无符号整型（unsigned int）。数据存储类型：ENVI标准栅格格式文件，且 是BIP或者BIL。5） 波谱范围：flaash能够做的数据光谱范围是0.4 2500卩（二） 基本参数设置1） 输入文件及输出路径设置，如图6所示。2） 传感器基本信息

27、设置选择传感器类型，成像中心点经纬度，成像时间，高度信息（成像区域和传感器 飞行高度），传感器这些都可以在数据自带信息文件里获得， 需要注意的是南半球和西半球要用负值。图6基本参数设置3)大气模型(Atmospheric Model )4)提供6种大气模型，如表2所示，可根据纬度和成像季节对照表3查找对应 的大气模型。(std atm-cm)Sub-Arctic Win ter (SAW)5180.42-16 C or 3 FMid-Latitude Win ter (MLW)10600.85-1 C or 30 FU.S. Standard (US)17621.4215 Cor 59 Sub

28、-Arctic Summer (SAS)25892.0814 C or 57 Mid-Latitude Summer (MLS)36362.9221 C or 70 Tropical (T)51194.1127 C or 80 表2六种标准的大气模型Latitude ( N) Jan.MarchMayJulySept.Nov.80SAWSAWSAWMLWMLWSAW70SAWSAWMLW MLWMLWSAW60MLWMLWMLW SASSASMLW50MLWMLWSASSASSASSAS40SASSASSASMLSMLSSAS30MLSMLSMLSTTMLS20TTTTTT10TTTTTT0TTTTTT-10TTTTTT-20TTTMLSMLST-30MLSMLSMLSMLSMLSMLS-40 SAS SAS SAS SAS SAS SAS-50 SAS SAS SAS MLW MLW SAS-60 MLW MLW MLW MLW MLW MLW-70 MLW MLW MLW MLW MLW MLW-80 MLW