大气校正(ENVI)

《大气校正(ENVI)》由会员分享，可在线阅读，更多相关《大气校正(ENVI)（26页珍藏版）》请在装配图网上搜索。

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

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

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

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

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

6、：l波段最小值(Band Minimum) lROI的平均值(Region Of Interest) l自定义值User Value4在Output Result to中选择File以与相应的输出路径和文件名，单击OK执行操作。下载 (55.67 KB)2010-1-28 16:37图2 Dark Subtraction Parameters面板Band Minimum2.2 基于统计学模型的反射率反演 基于统计学模型的反射率反演的方法主要有平场域法Flat Field，FF、对数残差法Log Residuals、部平均法Internal Average Relative Reflectanc

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

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

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

10、配。也可以用已经存在的系数对数据集进展定标。2.3 不变目标法相对大气校正 相对大气校正按照数学根底可以分为2种，非线性校正法和线性校正法。非线性校正法最典型的是直方图匹配，图像的直方图是图像中所有灰度值的概率分布。即将校正图像的直方图与参考图像的直方图进展匹配，使两幅图像具有一样或相近的灰度值概率分布，达到两幅图像上同名地物具有一样灰度值的目的。利用ENVI中Histogram Matching工具Display中，Enhance-Histogram Matching线性校正法有个前提假设：不同时相的图像灰度值之间满足线性关系，这种假设在近似情况下是成立的。这样就可以通过线性等式来描述不同时

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

12、像、已经经过准确配准、ENVI标准格式的Landsat TM5数据为例2000年和2001年，介绍线性校正法的操作步骤：第一步：PIF选择 选择一幅目视质量较好的图像作为基准图像2001年，另外一幅作为待校正图像2000年。在两个图像上选择一样区域的沥青房顶、砾石面、混凝土停机坪、洁净水体、混凝土、沙地等地物作为PIF，这些地物不会随时间的变化而变化。1在主模块中，选择File-Open Image File，打开两幅图像，并在Display中显示。2在其中一幅影像上点击右键，从快捷菜单中选择Geographic Link，将显示的两幅影像地理。3在显示2000年图像的主图像窗口中，选择Ove

13、rlay- Region of Interest，打开ROI Tool面板。4通过目视方式，从两幅图像找到光谱稳定、一样地物作为样本，用Polygon或者Point类型绘制感兴趣区。5在绘制一定数量感兴趣后(不宜太多，太多后面的回归运算量会很大)，在ROI Tool面板中，选择File- Output ROIs to ASCII。6回到ROI Tool面板中，选择Options- Reconcile ROIs via Map，将前面绘制的ROI转接到2001年的图像上，类似5步的方法将基准图像的ROI对应像素位置和像元值输出为文本文件。分别用记事本打开上面步骤得到两个文本文件，这样我们得到了2

14、000年的图像和2001年的图像相对应伪不变特征要素PIF的像素值。从文本文件中可以看到，两个时相图像中每一个波段的像素值是一一对应关系，刚好对应式1中的x和y。第二步中就是利用这些像素值，根据最小二乘回归分析法获得式中的a和b两个参数。第二步、线性关系式求解使用最小二乘回归的方法来求解线性回归式a和b参数，如表1。表1回归解算的a和b值波段增益a偏移bBand1-34Band2-18Band3-9Band4-16Band54Band73第三步、线性变换利用表1中的a和b值，在ENVI的Band Math工具对待校正图像做线性变换，然后利用Layer Stacking工具将线性变换结果组合成一

15、个多波段文件。不变目标法相对大气校正操作过程已经完成。2.4 热红外大气校正 ENVI提供Thermal Atm Correction工具，可以近似去除热红外辐射数据中的大气影响。在进展大气校正之前，为了得到最好的结果，必须将热红外数据定标为比辐射率数据TIMS的热红外数据必须被转化为辐射亮度数据，并且待校正数据波长在8-14µm之间。下面以ASTER L1A的热红外波段为例，操作过程如下：1打开ASTER，在波段列表中按照波长自动归为4组，2组可见光-近红外VNIR、短波红外SWIR和热红外(IIR)，并根据头文件信息自动定标为辐射亮度值单位W/m2 *µm*sr。2

16、在主菜单中，从以如下表中选择一种方式lBasic Tools - Preprocessing-Calibration Utilities-Thermal Atm CorrectionlBasic Tools-Preprocessing-Data-Specific Utilities-Thermal IR - Thermal Atm CorrectionlBasic Tools- Preprocessing- Data-Specific Utilities -TIMS-Thermal Atm Correction在Thermal Correction Input File对话框中，选择热红外数据

17、。3在Thermal Atm Correction Parameters面板中图3，需要填写以下参数：l数据缩放系数Data Scale Factor：1。将输入数据的单位缩放为W/m2 *µm*sr。l波长单位Wavelength Units：Micrometers。l设定外表温度估算衰退像元Regression Pixels：Allu选择“All，将使用整个输入波长围的亮度温度最大值对每个像元的外表温度进展估算。u选择“Max Hit，仅对那些在特定波长具有最大亮度温度值的像元进展外表温度估算，所说的特定波长是指包含最多具有最大亮度温度的像元的波长围。l散点图拟合技术Fitti

18、ng Technique：Top of Binsu选择“Top of Bins，将会使曲线向辐射率与亮度温度的散点图的上部拟合，散点图的上部对应着发射率接近1的像元，使用该技术的拟合线是通过对散点图上部5%的数据做了标准的最小平方回归得到的。注：该技术易受发生在散点图顶部区域的传感器噪声的影响u选择“Normalized Regression，将会先使用标准的最小平方回归把曲线向辐射率与亮度温度的散点图拟合，然后将拟合线的残差与正态概率分布图相比拟，在正态图中对残差再进展另一个回归，NESR(noise equivalent sensor response)距离的点将被认为是奇异点而被删除，最

19、后利用减少的像元集在散点图上进展一个最终的回归。注：该方法使用散点图中除奇异点以外的所有点，并且不仅仅把曲线向散点图的上部发射率近似为1拟合l设置输出增益与偏移参数为文件Output Gain/Offset File：可选。l设置是否绘制大气透射和上行辐射光谱结果(Plot Transmission/Upwelling?)：Yes。4选择结果输出路径与文件名，单击OK执行校正过程。下载 (77.16 KB)2010-1-28 16:37图3Thermal Atm Correction Parameters面板2.5 QUAC快速大气校正工具 快速大气校正工具简称QUAC自动从图像上收集不同物质

20、的波谱信息，获取经验值完成高光谱和多光谱的快速大气校正图4。它得到结果的精度近似FLAASH或者其他基于辐射传输模型的+/-15%。目前它支持的多光谱和高光谱波谱围是(0.42.5 m)。传感器包括AISA, ASAS, AVIRIS, CAP ARCHER, PASS, HYCAS, HYDICE, HyMap, Hyperion, IKONOS, Landsat TM, LASH, MASTER, MODIS, MTI, QuickBird, RGB, 以与 unknown sensor。下载 (52.06 KB)2010-1-28 16:37图4 快速大气校正流程图QUAC的输入数据可以

21、是辐射亮度值、表观反射率、无单位的raw数据。可以是任何数据储存顺序BIL/BIP/BSQ和储存类型，多光谱和高光谱传感器数据的每个波段必须有中心波长信息。QUAC的操作非常简单，如下：1在ENVI主菜单中，选择以下方式启动lBasic Tools- Preprocessing- Calibration Utilities- QUick Atmospheric CorrectionlSpectral- QUick Atmospheric CorrectionlSpectral- Preprocessing- Calibration Utilities- QUick Atmospheric Co

22、rrection在文件输入对话框中选择校正的图像文件。2打开QUick Atmospheric Correction Parameters面板图5，在Sensor Type中选择相应的传感器类型，选择文件名和路径输出。下载 (43.67 KB)2010-1-28 16:37图5 QUick Atmospheric Correction Parameters面板2.6 FLAASH大气校正 FLAASH是基于MODTRAN4+辐射传输模型，MODTRAN模型是由进展大气校正算法研究的领先者SpectralSpectral Sciences, Inc和美国空军实验室Air Force Resear

23、ch Laboratory共同研发。ITT VIS公司负责集成和GUI设计。特点l支持传感器种类多，包括多光谱的ASTER，AVHRR，GeoEye-1，IKONOS，IRS，Landsat，MODIS，SeaWiFS，SPOT，QuickBird，RapidEye等，高光谱HyMAP、AVIRIS、CASI、HYDICE、HYPERION(EO-1)、AISA等。可以通过自定义波谱响应函数支持更多的传感器。工程化应用价值比拟明显。lFLAASH采用了MODTRAN4+辐射传输模型，该算法精度高。任何有关影像的标准MODTRAN大气模型和气溶胶类型都可以直接使用。l通过影像像素光谱上的特征来估

24、计大气的属性，不依赖遥感成像时同步测量的大气参数数据。l可以有效地去除水蒸气/气溶胶散射效应，同时基于像素级的校正，矫正目标像元和邻近像元交叉辐射的“邻近效应。l对由于人为抑止而导致波谱噪声进展光谱平滑处理。作为结果，除了真实地表反射率外，还可以得到整幅图像的能见度、卷云与薄云的分类影像、水气含量数据。2.使用ENVI大气校正模块ENVI大气校正模块的使用主要又以下7个方面组成：1、输入文件准备，2、根本参数设置，3、多光谱数据参数设置，4、高光谱数据参数设置，5、高级设置，6、输出文件，7、处理结果。下面介绍这7个方面容。(一)输入文件准备1)支持传感器类型高光谱包括：HyMAP、AVIRI

25、S、CASI、HYDICE、HYPERION(EO-1)、AISA等；多光谱包括：ASTER，AVHRR，GeoEye-1，IKONOS，IRS，Landsat，MODIS，SeaWiFS，SPOT，QuickBird，RapidEye等，航空：860nm-1135nm波长围2)数据是经过定标后的辐射亮度辐射率数据，单位是：W/cm2*nm*sr。3)数据带有中心波长wavelenth值，如果是高光谱还必须有波段宽度FWHM,这两个参数都可以通过编辑头文件信息输入Edit Header。4)数据类型支持四种数据类型：浮点型floating、长整型(long integer )、整型intege

26、r和无符号整型 (unsigned int)。数据存储类型： ENVI标准栅格格式文件，且是BIP或者BIL。5)波谱围：flaash能够做的数据光谱围是2500m。(二)根本参数设置1)输入文件与输出路径设置，如图6所示。2)传感器根本信息设置选择传感器类型，成像中心点经纬度，成像时间，高度信息成像区域和传感器飞行高度，传感器这些都可以在数据自带信息文件里获得，需要注意的是南半球和西半球要用负值。下载 (97.86 KB)2010-1-28 16:37图6根本参数设置3)大气模型Atmospheric Model提供6种大气模型，如表2所示，可根据纬度和成像季节对照表3查找对应的大气模型。M

27、odel Atmosphere Water Vapor (std atm-cm) Water Vapor (g/cm2) Surface Air Temperature Sub-Arctic Winter (SAW) 518-16 C or 3 F Mid-Latitude Winter (MLW) 1060-1 C or 30 F U.S. Standard (US) 176215 C or 59 Sub-Arctic Summer (SAS) 258914 C or 57 Mid-Latitude Summer (MLS) 363621 C or 70 Tropical (T) 51192

28、7 C or 80 表2六种标准的大气模型Latitude (N) Jan. March May July Sept. Nov. 80SAWSAWSAWMLWMLWSAW70SAWSAWMLWMLWMLWSAW60MLWMLWMLWSASSASMLW50MLWMLWSASSASSASSAS40SASSASSASMLSMLSSAS30MLSMLSMLSTTMLS20TTTTTT10TTTTTT0TTTTTT-10TTTTTT-20TTTMLSMLST-30MLSMLSMLSMLSMLSMLS-40SASSASSASSASSASSAS-50SASSASSASMLWMLWSAS-60MLWMLWML

29、WMLWMLWMLW-70MLWMLWMLWMLWMLWMLW-80MLWMLWMLWSAWMLWMLW表3数据经纬度与获取时间对应的大气模型4)水气反演Water Retrieval水气反演设置，采用两种方式对水气进展反演：a)利用水气反演模型恢复影像中每个像元的水气量使用水气反演模型，数据必须具有15nm以上波谱分辨率，且至少覆盖以下波谱围之一：l1050-1210 nm (对应 1135 nm)l870-1020 nm (对应940 nm)l770-870 nm (对应820 nm)对于大多数多光谱传感器，水气反演默认显示的是NO，因为大多数传感器没有适当的波段来补偿水气的影响。b)单一

30、的水气因数用于整体影像，默认是1，对于多光谱数据使用水气反演模型，可以在多光谱设置中手动设置水气波段5)气溶胶模型Aerosol Modela)提供四种标准MODTRAN气溶胶模型Rural乡村、Urban城市、Maritime海洋、Tropospheric对流层，能见度在40km以上b)两种气溶胶反演方法2-BandK-T方法(类似模糊减少法)，如果没有找到适应的黑值一般是阴影区或者水体，系统将采用能见度值来计算；所以即使选择了该选项也要给能见度。选择None，采用能见度值参与气溶胶反演，能见度值大约参考值参见表4天气条件 能见度 晴朗 40 to 100 km 中等雾、阴霾 20 to 3

31、0 km 厚雾、阴霾15 km 或者更少 表4天气条件与能见度对照表6)光谱打磨(高光谱) Spectral Polishing对高光谱数据的光谱进展打磨处理，对波谱曲线进展微调，使波谱曲线更加近似于真实地物的波谱曲线。7)重新定标波长Recalibrating the Input Wavelengths 各个波段的中心波长对大气校正结果影像比拟大。(三)多光谱数据参数设置当根本设置里设置了水气反演以与气溶胶反演，相应的要在此设置相关参数单击Multispectral Settings。有两种设置方式：文件方式和图形方式，一般选择图形方式。1)水气反演模型参数Water Retrieval11

32、35 nmabsorption1117 - 1143 nmreference upper wing1184 - 1210 nmreference lower wing1050 - 1067 nm940 nmabsorption935 - 955 nmreference upper wing870 - 890 nmreference lower wing995 - 1020 nm820 nmabsorption810 - 830 nmreference upper wing850 - 870 nmreference lower wing770 - 790 nmAerosol RetrievalK

33、T upper2100 - 2250 nmKT lower640 - 680 nmCloud Maskingcirrus clouds1367 - 1383 nm表5 推荐使用波长与波段对应关系2)气溶胶模型参数设置用气溶胶模型要求数据波段覆盖660nm和2100nm波谱，具体波段选择可以从表5所列中选择。下载 (78.19 KB)2010-1-28 16:37图7 多光谱数据参数设置(四)高光谱数据参数设置3)自动选择通道定义(推荐)4)设置通道定义下载 (33.88 KB)2010-1-28 16:37图8高光谱设置(五)高级设置下载 (75.89 KB)2010-1-28 16:37图9

34、高级参数设置在高级设置面板里可以设置以下几个参数：1)光谱定义文件：置AVIRIS、HYMAP、HYDICE、HYPERION、CASI、AISA。2)气溶胶厚度系数：用于技术邻域效应围。一般值为12km。3)CO2混合比率：默认为390ppm。4)使用领域纠正Use Adjacency Correction。5)使用以前的MODTRAN模型计算结果。6)设置MODTRAN模型的光谱分辨率推荐值5 cm-1)。7)设置MODTRAN多散射模型。提供三种模型供选择Isaacs，DISORT和 Scaled DISORT。默认是Scaled DISORT和streams为8，这种模型对于小于100

35、0nm具有较高的精度；8)天顶角方位角针对非星下点传感器。9)输出反射率缩放系数Output Reflectance Scale Factor：为了降低结果储存空间，默认反射率乘于10000。(六)输出文件经过大气校正模块处理，可以得到以下结果：l外表反射率影像l水气含量数据l云图l日志文件lFLAASH大气校正模板文件(七)处理结果此模块的波谱复原精度高，使用大气校正模块进展大气校正，能高保真地恢复地物波谱信息。下载 (131.98 KB)2010-1-28 16:37图10 校正前后效果图3.常见问题1、如图11所示错误，为因为没有设置输出反射率文件名，解决方法是单击Output Reflectance File按钮，选择反射率数据输出目录与文件名。下载 (33.71 KB)2010-1-28 16:37图11 ENVI Error2、在做modis数据时候，不能使用领域纠正Use Adjacency Correction。3、当数据有地理坐标时候，像元大小要也参数中设置保持一致。如你tm数据是25米的，那么选择Landsat时候，像元大小默认是30米，这个时候如果改成25米会报错。4、得到的结果有负值，这个不好解释，可以理解为误差或者其他因素。或者对太强吸收区域，如水体。26 / 26