当前位置:首页 » 图片资讯 » 大气校正后的图片是什么格式
扩展阅读
美女健身跳河视频 2023-08-31 22:08:21
西方贵族美女照片真人 2023-08-31 22:08:15

大气校正后的图片是什么格式

发布时间: 2022-06-16 09:17:35

A. envi大气校正的时候图像格式转换之前为什么要将波段运算b1/10.0

这个是因为物理量的单位换算需要
经过定标后的辐射亮度(辐射率)数据,单位是:(μW)/(cm2*nm*sr)
根据波段波长之间的关系,换算过程是需要对单位进行换算,这里差了10倍,所以要÷10.

B. ENVI下FLASSH大气校正后的数据范围正常是少,输出格式是float还是什么格式的,谢谢

大气校正后默认乘以10000,因此范围为0-10000,输出的应该是整型的

C. 大气校正

遥感器接收目标辐射或反射的电磁波所形成的遥感原始图像与目标相比是失真的,这是因为在太阳-大气-目标-大气-遥感器的光线传播路径中,许多因素的影响造成接收的信号不能准确的反映地表物理特征。这些因素归结为以下几个方面:

(1)大气内容物的影响

大气主要由大气分子和气溶胶组成,这两者的影响行为是不相同的。大气分子瑞利散射、气溶胶的Mie散射;大气分子与气溶胶的吸收及两者的耦合作用。一方面,大气的吸收导致消光,减少了辐射量,降低了图像对比度,使图像变得暗淡;另一方面,大气散射导致的程辐射,增加了辐射量。

(2)表面因素的贡献

在一般的应用中,为了简化计算,假定地表为朗伯体,反射与方向无关。事实上任何表面在物理特性与物质结构上都不是理想朗伯体,因此认为地面是朗伯体会带来误差,而当地表方向反射特性突出时,假设地面是朗伯体的大气纠正方法精度受到限制。另一个因素是由于大气散射的存在,邻近像元的反射光也会进入目标视场从而影响辐射量,即交叉辐射。

(3)地形因素的影响

目标高度与坡向会对辐射造成影响。

(4)太阳辐射光谱的影响

太阳本身是一个黑体,其光谱辐射按照普朗克定律有一定的形状,这个因素在反射率反演中需要予以考虑。

由以上可知,大气对光学遥感的影响是十分复杂的。为此,学者们尝试着提出不同的大气纠正模型来模拟大气的影响。但是对于任一幅图像,其对应的大气数据几乎是永远变化的,且难以获得,因而应用完整的模型纠正每个像元是不可能的。最早的大气纠正方法是从图像本身来估计大气参数,反复运用大气模拟模型进行纠正。结合地面实况数据进行大气校正是另一种方法,其包括两种类型:一种是通过地面测定大气参数(如可见光近红外的气溶胶的密度及红外区域的水汽浓度),再结合辐射传输方程作近似求解;另一种是测得地面目标物的反射率,再与图像数据进行比较来消除大气的影响。地面同步测量有助于提高精度,但是却需要人力物力,且应用区域也有限。此外还有一些大气纠正的方法。例如在同一平台上,除了安装获取目标图像的遥感器以外,也安装上专门测量大气参数的遥感器,利用这些数据进行大气校正。

3.4.1 基于影像特征的校正模型

基于图像特征的相对校正法是在没有条件进行地面同步测量的情况下,借用统计方法进行图像相对反射率转换。从理论上来讲,基于图像特征的大气校正方法都不需要进行实际地面光谱及大气环境参数的测量,而是直接从图像特征本身出发消除大气影响,进行反射率反演,基本属于数据归一化的范畴。精确的大气校正需要精确的测量大气参数和复杂的运算,这些在许多遥感应用中,往往很难满足。并且在某些应用中不一定需要绝对的辐射校正。此时,这种基于图像的相对校正就能满足其要求。

基于图像特征的相对校正法主要有内部平均法、平场域法、对数残差法等。

(1)内部平均法

假定一幅图像内部的地物充分混杂,整幅图像的平均光谱基本代表了大气影响下的太阳光谱信息。因而,把图像DN值与整幅图像的平均辐射光谱值的比值确定为相对反射率光谱,即

ρλ = Rλ /Fλ (3.21)

式中:Rλ为像元在该波段的辐射值;Fλ为整幅图像的平均辐射光谱值;ρλ为该像元的相对反射率。

(2)平场域法

平场域法是选择图像中一块面积大且亮度高而光谱响应曲线变化平缓的区域,利用其平均光谱辐射值来模拟飞行时大气条件下的太阳光谱。将每个像元的DN值与该平均光谱辐射值的比值作为地表反射率,以此来消除大气的影响。

ρλ = Rλ /Fλ (3.22)

式中:Rλ为像元在该波段的辐射值;Fλ为平场域的平均辐射光谱值;ρλ为该像元的相对反射率。

利用平场域消除大气影响并建立反射率光谱图像有两个重要的假设条件:一个是平场域自身的平均光谱没有明显的吸收特征;另一个是平场域辐射光谱主要反映的是当时大气条件下的太阳光谱。

平场域模型已广泛应用于遥感数据处理中,它是在内部平均法模型基础上发展起来的,这种模型克服了内部平均法模型易受像幅内吸收特征影响而出现假反射峰的弱点,而且计算量更小,其不足之处在于选取光谱地理平台单元时,会引入人为的误差,而且需要对研究区内地物光谱有一定的先验了解,当选取具有不同反射率等级的地理平台单元时,会引出不同处理结果。当研究区位于山区或其他地形起伏较大的复杂地区时,选择地理平台单元较为困难。

(3)对数残差法

对数残差法的意义是为了消除光照及地形因子的影响。按照一定的规则调节每个像元值,使其在每一个被选定的波段上的值等于整个图幅的最大值,然后对每一个波段减去其归一化后的平均值。假设有:

DNij = TiRijIj (3.23)

式中:DNij为像元i的j波段的灰度值;Ti 为像元i处表征表面变化的地貌因子,对确定的像元所有的波段该值都相同;Rij为像元i波段j的反射率;Ij为波段j的光照因子。

由表3.12我们可以看出,以上三种方法中,只有残差图像法是真正意义上的辐射校正。

表3.12 高光谱基于图像特征的相对校正法对各种影响辐射的物理因素的补偿能力比较

3.4.2 地面线性回归经验模型

基于地面线性回归经验模型法是一个比较简便的定标算法,国内外已多次成功地利用该模型进行遥感定标实验。它首先假设地面目标的反射率与遥感器探测的信号之间具有线性关系,通过获取遥感影像上特定地物的灰度值及其成像时相应的地面目标反射光谱的测量值,建立两者之间的线性回归方程式,在此基础上对整幅遥感影像进行辐射校正。该方法数学和物理意义明确,计算简单,但必须以大量野外光谱测量为前提,因此成本较高,对野外工作依赖性强,且对地面定标点的要求比较严格。这种方法仅适用于地面实况数据特定的地区及时间。

3.4.3 利用波段特性进行大气校正

在利用卫星遥感中,有相当部分的大气散射光未经过地物反射,通过大气吸收后,直接进入传感器。我们叫这种辐射为程辐射。严格地说,程辐射的大小与像元位置有关,随大气条件、太阳照射方向和时间变化而变化,但因其变化量微小而忽略。可以认为,程辐射度在同一幅图像的有限面积内是一个常数,其值的大小只与波段有关。一般来说,程辐射度主要来自米氏散射,即散射主要发生在短波波段,其散射强度随波长的增大而减小,到红外波段基本接近于零。可以把近红外波段作为无散射影响的标准图像,通过对不同波段图像的对比分析来计算大气影响。根据这个原理主要有三种方法:单影像直方图调整法、单影像回归分析法和多时相影像归一化分析法。

(1)单影像直方图调整法

采用单影像直方图调整方法的前提条件是在一幅影像中存在某种地物如深海水体、高山背阴处等,其辐亮度值或反射率接近于0,这时其图像直方图的最小值就应该为0,如果不为0,就认为是大气散射导致的。

(2)单影像回归分析法

假定某红外波段,程辐射影响接近于零,设为波段a,现需要找到其他波段相应的亮度最小值,这个值一定比a波段的亮度最小值大一些,设为波段b。分别以a,b波段的像元亮度值为坐标,作二维光谱空间,两个波段中对应像元在坐标系内用一个点表示。由于波段之间的相关性,通过回归分析在众多点中一定能找到一条直线与波段b的亮度Lb 轴相交,即用最小二乘法拟合出一条直线,回归方程为

Lb = KLa + c (3.24)

式中:c为拟合的直线在Lb轴的截距;K为拟合直线的斜率。

高光谱遥感技术原理及矿产与能源勘查应用

高光谱遥感技术原理及矿产与能源勘查应用

式中:La为假定波段亮度最小值;

为所选黑区的像元均值。C是波段a中的亮度值为0时,在波段b中所具有的亮度值。可认为C就是波段b的程辐射值。然后将波段b中的所有像元值都减去这个截距值C,来去掉波段b中的程辐射值。

(3)多时相影像归一化分析法

多时相影像归一化首先要选取基准影像(设为b),然后对不同时相的所有其他影像的光谱特征进行转换,使它们具有与基准影像基本相同的辐射量级。多时相影像归一化分析方法的一个重要步骤是选取伪不变特征(Pseudo-Invariant Features,PIFs),也称为辐射地面控制点。伪不变特征具有如下特点:尽管某些变化是不可避免的,伪不变特征的光谱特性应该随时间变化很小,如深海水体、裸土。大屋顶或其他同质地物都是不错的选择;在一景影像中,伪不变特征与其他地方的高程应该大致相同,山顶伪不变特征在估计近海面大气条件中的作用不大,因为大气中的多数气溶胶都出现在低于1000 m的大气中;伪不变特征包含的植被应尽可能少,由于环境胁迫和气候周期的影响,植被光谱反射率会随时间变化;伪不变特征应该选在相对平坦的区域,使太阳高度角的逐日变化与所有归一化目标的太阳光直射光束之间具有增加或减小的比例。

利用基准影像与其他时相影像的PIFs光谱特性之间的联系进行回归分析。该方法是假定时相b-1或b+1的影像像元与基准影像b相同位置上的像元是线性相关的。这意味着,采样像元的光谱反射特性在这一段时间内没有发生变化,所以多时相影像回归分析的关键是选取伪不变特征。

地面覆盖的遥感分类能力依赖于遥感亮度值(BV)和实际地表条件的稳定联系。然而,太阳高度角、日地距离、各种不同传感器系统的探测器定标差异、大气条件和太阳-目标传感器的几何关系等因素会影响像元亮度值。影像归一化减少了由非地表因素引起的像元亮度值变化,使不同时相的像元亮度值变化与地表条件的实际变化相联系。归一化处理使得从基准影像中得到的像元分类可用于其他的归一化影像上。

3.4.4 大气辐射传输模型理论方法

1972年,Turner与Spencer提出的通过模拟大气-地表系统来评估大气影响的方法,可作为最早的大气辐射传输模型之一,当时研究的重点在于消除大气对影像对比度的影响。20世纪80年代,许多学者对卫星影像的大气校正研究做了大量工作,在模拟地-气过程的能力上有了很大提高,发展了一系列辐射传输模型,例如我们熟知的LOWTRAN系列模型和5S模型。

自1990年以来,许多的辐射传输模型被用于大气校正算法中,涌现出一大批新的大气校正模型,其中有的方法使用一些先进的数学算法提高计算速度(如6S),试图寻找精度与速度的最佳平衡点。

基于大气辐射传输理论的大气纠正模型主要考虑的问题有以下几个方面:

1)构成大气的气体分子和气溶胶的散射和吸收特性及两者耦合效果的研究。其中,各模型主要考虑的是吸收及气溶胶散射。大气辐射传输模型中用到的大气参数包括气温、气压、水汽含量、臭氧含量、能见度、水平气象视距、灰尘颗粒度等,这些参数用于计算辐射传输方程中大气的吸收透过率与散射透过率,以及气溶胶光学厚度,因此输入大气参数的精度直接影响大气校正的最终结果。同步实地观测可以为大气校正提供所需的大气参数,但同步实地观测需耗费大量人力物力,且对历史数据无能为力。为此,6S和MODTRAN中提供了一系列既定参数供用户选择,这些参数是对大量观测数据统计分析得到的,旨在模拟遥感器过境时的大气状况,但这与实际的大气状况存在差距。

2)地表特性的假设。高精度的大气校正必须考虑地表非朗伯体特性。在6S中可以选择均一非朗伯体模型。

3)模型中算法的选择。更精确的算法往往会伴随巨大的计算量,以往大气纠正的过程中,学者往往会牺牲一定的精度来满足计算速度的需求,现在随着计算机科技的发展,越来越多的模型选择了复杂而更精确的算法来满足高精度的需求。

在已有的模型中,最着名的辐射传输模型是MODTRAN和6S。分别是在对LOWTRAN与5S改进的基础上发展而来。由于高光谱相机波段范围是400~2500nm,拟分别采用6S 辐射传输模型和MODTRAN辐射传输模型进行大气纠正。6S源代码开放,可以很方便地进行修改和移植;MODTRAN可供自定义的参数多,均适合于相应的地表反射率反演系统开发。利用若干典型区域的长期地基观测数据(如AERONET观测站点所在区域),建立起局地气溶胶模式和类型,结合6S和MODTRAN分别建立针对这些区域的大气纠正模式。同时,在实验验证的基础上对重点区域大气辐射传输方程进行简化,在不降低反演精度的前提下,减少运算次数,提高高光谱图像在这些区域的地表反射率反演效率。

6S模型是1997年由Vermote和Tanre等人用Fortran语言编写的适用于太阳反射波段(0.25~4μm)的大气辐射传输模式。由于计算机水平和其他相关知识的发展,6S模型对5S模型提出了一系列的改进。主要改进如下:

1)在5S模型中,瑞利散射的大气函数ρ,T,S被制成表,给使用带来不便。在6S中,用满足精度的解析表达式代替。

2)6S模型选用高精度的SOS模型代替原有方法处理分子和气溶胶散射。SOS模型可以精确模拟机载遥感,并且提供处理非朗伯体(BRDF)临近问题所需的输入参数。

3)5S模型假定吸收作用与散射作用可以耦合,就像吸收粒子位于散射层的上面一样。6S假设散射和吸收互不影响,主要考虑水汽吸收和气溶胶散射的三种极端耦合情况:水汽吸收粒子覆在气溶胶层之上;水汽吸收粒子在气溶胶层之下;有一半水汽吸收粒子与气溶胶辐射路径混在一起。

4)5S中,气体吸收传输用的是随机波段模型。这个模型有两个主要问题:首先,使用的是AFGL在1982年公布的大气吸收线性参数,并没有考虑太阳反射光谱段的一些其他吸收气体;其次,采用20 cm-1的波段间隔(过大)模拟宽波段辐射计(如1000 cm-1)的吸收,这个较宽的波段间隔不适用于模拟更高光谱分辨率(如100 cm-1)光谱仪器的吸收情况。在6S中,不仅考虑新的吸收分子种类的影响,并且气体的吸收以10 cm-1的光谱间隔来计算。

5)为了兼顾计算效率,5S代码仅模拟海平面上均匀朗伯体目标的反射率。在6S中,目标高程Zt 可作为一项输入:可依据Zt 去除目标高度以下的大气层,计算新的大气廓线;由于Zt对主要分布在低层大气中的H2 O产生很大的影响,故可根据Zt 重新计算H2 O含量,同理,可根据Zt 重新计算气溶胶的含量;6S将光学厚度视为目标高度处压强的比例函数,从而很高精度的计算了Zt 对分子光学厚度的影响。

6S模型定义了地表的反射率模型,包括均一地表和非均一地表两种情况,在均一地表中又考虑了无方向性反射问题,在考虑方向性问题时用了九种不同的模型。利用较高精度的新模型解释BRDF作用和临近效应。

6)6S对5S数据库的改进:

光谱积分步长达到了2.5nm(相比于原来的5nm)。

增加了新的吸收气体(如CO2,N2O,CH4),6S模型用HITRAN数据库以10cm-1分辨率计算波段吸收参数。

IRC定义的四种基本气溶胶微粒以更好的步长重新计算一次。

且6S中新加了5S中难以计算的气溶胶模型(平流层型、沙漠型,以及生物燃烧产生的气溶胶类型)。

6S给定了九种比较成熟的BRDF供用户选择,也可以自定义BRDF函数,作为参数输入到6S,验证研究反射率与地表BRDF的关系(表3.13 ,表3.14)。

表3.13 6S模型输入参数

表3.14 6S模型输出参数

LOWTRAN和MODTRAN模型是由美国空军地球物理实验室研制的大气辐射模拟计算程序,在遥感领域被广泛应用于图像的大气校正。

LOWTRAN是一个光谱分辨率为20 cm-1的大气辐射传输实用软件,它提供了六种参考大气模式的温度、气压、密度的垂直廓线;水汽、臭氧、甲烷、一氧化碳、一氧化二氮的混合比垂直廓线;其他13种微量气体的垂直廓线;城乡大气气溶胶、雾、沙尘、火山喷发物、云、雨的廓线;辐射参量(如消光系数、吸收系数、非对称因子的光谱分布)及地外太阳光谱。它可以根据用户的需要,设置水平、倾斜及垂直路径,地对空、空对地等各种探测几何形式,适用对象广泛。LOWTRAN的基本算法包括透过率计算方法,多次散射处理和几何路径计算。

(1)多次散射处理

LOWTRAN采用改进的累加法,自海平面开始向上直至大气的上界,全面考虑整层大气和地表、云层的反射贡献,逐层确定大气分层每一界面上的综合透过率、吸收率、反射率和辐射通量。再用得到的通量计算散射源函数,用二流近似解求辐射传输方程。

(2)透过率计算

该模型在单纯计算透过率或仅考虑单次散射时,采用参数化经验方法计算带平均透过率,在计算多次散射时,采用k分布法。

(3)光线几何路径计算

考虑了地球曲率和大气折射效应,将大气看作球面分层,逐层考虑大气折射效应。

MODTRAN模型可以计算0到50000 cm-1的大气透过率和辐亮度,它在440 nm到无限大的波长范围精度是2cm-1 ,在22680cm-1到50000cm-1紫外波(200~440nm)范围的精度是20cm-1 ,在给定辐射传输驱动、气溶胶和云参数、光源与遥感器的几何立体对和地面光谱信息的基础上,根据辐射传输方程来计算大气的透过率以及辐亮度。

开发MODTRAN是为了改进 LOWTRAN 的光谱分辨率,MODTRAN 将光谱的半高全宽度(full width half maximum,FWHM)由LOWTRAN的20cm-1减少到2cm-1 ,目前的MODTRAN4.0它的光谱分辨率已经达到2 cm-1 ,改进了瑞利散射和复折射系数的计算精度,增加了DISORT计算太阳散射贡献的方位角相关选项,并将七种BRDF模型引进到模型中,使地表的参数化输入成为可能。

MODTRAN以卡片的形式来进行参数设置,操作起来清晰简洁,可以在文本格式下直接对其输出输入参数文件进行修改。

D. ENVI在大气校正是FLAASH中的Lat 和Lon的值要怎么在元数据文件中找到啊

1、 打开tm原始影像数据
ENVI > file > open image file > ‘L5142032_03220090725_MTL.txt’ > 打开文件,如下:

图2 自动读取数据列表展示
2、 TM数据辐射定标
ENVI > basic tools > preprocessing > calibration utilities > Landsat calibration,弹出如下对话框,图3:

图3 辐射定标参数设置对话框

3、 储存顺序调整
Flassh大气校正对于波段存储的要求为:BIL,BIP格式,上述计算得到的存储方式为BSQ,在此进行波段存储顺序的转化,具体操作如下:
ENVI > basic tools > convert data (BSQ ,BIL ,BIP)

图 4 存放顺序转换
4、 Flaash校正参数设置
大气校正的前期准备工作完毕,现在进行校正参数的设置:
ENVI > basic tools > preprocessing > calibration utilities > FLAASH,弹出对话框:

图 5 FLAASH参数界面设置

图 6 多光谱设置对对话框
根据上述图中的参数设置,然后点击ok,运行flaash大气校正。

E. 用ENVI做大气校正必须用fst格式的文件吗

不用,该有的信息有就可以。

F. 什么叫大气校正

这是一个专业名词,为您解答,请参考:
中文名称:大气校正
英文名称:atmospheric correction
应用学科:测绘学(一级学科);摄影测量与遥感学(二级学科)
概念:消除遥感图像中由大气散射引起的辐射误差的处理过程。
遥感所利用的各种辐射能均要与地球大气层发生相互作用-或 散射、或吸收,而使能量衰减,并使光谱分布发生变化。大气的衰减作用对不同波长的光是有选择性的,因而大气对不同波段的图像的影像是不同的。另外,太阳-目标-遥感器之间的几何关系不同,则所穿越的大气路径长度不同,使图像中不同地区地物的像元灰度值所受大气影响程度不同,且同一地物的像元灰度值在不同获取时间所受大气影响程度也不同。消除这些大气影响的处理,称为大气校正。
即使遥感系统工作正常,获取的数据仍然带有辐射误差。两种最重要的环境衰减是1)由大气散射和吸收引起的大气衰减;2)地形衰减。然而,在所有的遥感应用中都进行大气校正可能没有必要。是否进行大气校正,取决于问题本身、可以得到的遥感数据的类型取的历史与当前实测大气信息的数量和从遥感数据中提取生物物理信息所要求的精度。

G. 老师好,我做完大气校正后的landsat8图像无法打开,大气校正过程应该没错,请老师指导下,谢谢您!

这个是envi文件系统,不用管它,可以打开另存为tif等格式的

H. 大气校正是什么

即使遥感系统工作正常,获取的数据仍然带有辐射误差。两种最重要的环境衰减是1)由大气散射和吸收引起的大气衰减;2)地形衰减。然而,在所有的遥感应用中都进行大气校正可能没有必要。是否进行大气校正,取决于问题本身、可以得到的遥感数据的类型取的历史与当前实测大气信息的数量和从遥感数据中提取生物物理信息所要求的精度。

大气校正的目的是消除大气和光照等因素对地物反射的影响,广义上讲获得地物反射率地表温度等真实物理模型参数;狭义上是获取地物真实反射率数据。用来消除大气中水蒸气二氧化碳、甲烷和臭氧等物质对地物反射的影响,消除大气分子和气溶胶散射的影响。大气校正同时也是反演地物真实反射率的过程。

I. 遥感中BIP、BSQ、BIL优缺点是什么

BSQ在图像显示速度上更快,而BIL、BIP在图像处理上又明显优于BSQ格式;另外像FLAASH大气校正是基于像元处理的,BSQ是不满足的,另外像图像融合以及图像的分类等都是基于像元进行处理的。

BIL,BSQ和BIP文件是栅格数据记录的方式不同。

BSQ:波段顺序格式。

BIP:波段按像元交叉格式。

BIL:波段按行交叉格式。

各种可以互相转换 在ENVI中就可以实现转换,ENVI主菜单--basic tools--convert data(/BSQ/BIP/BIL)。

遥感(remote sensing)是指非接触的,远距离的探测技术。一般指运用传感器/遥感器对物体的电磁波的辐射、反射特性的探测。遥感是通过遥感器这类对电磁波敏感的仪器,在远离目标和非接触目标物体条件下探测目标地物。

获取其反射、辐射或散射的电磁波信息(如电场、磁场、电磁波、地震波等信息),并进行提取、判定、加工处理、分析与应用的一门科学和技术。

J. 实验十九 遥感图像辐射校正

一、实验目的

通过对ENVI的FLAASH 功能运用,掌握对Landsat卫星遥感影像做辐射定标与FLAASH大气校正的技术操作,加深对辐射定标和大气校正原理的理解。

二、实验内容

①桂林市TM 遥感影像辐射定标;②桂林市TM 遥感影像FLAASH大气校正处理。

三、实验要求

①明确ENVI的Lndsat辐射定标各项参数的意义及作用;②明确ENVI FLAASH大气校正处理各步操作的作用;③对FLAASH大气校正前后同名点的植被波谱曲线差异进行比较分析。编写实验报告。

四、技术条件

①微型计算机;②国际分幅127-43 TM 影像;③ENVI软件;④Photoshop软件(ver.6.0以上)和ACDSee软件(ver.4.0以上)。

五、实验步骤

辐射校正包括辐射定标与大气校正两部分。辐射定标是将传感器记录的电压或数字量化值(DN)转换成绝对辐射亮度值(辐射率)。大气校正就是消除大气和光照等因素对地物反射的影响,获得地物反射率、辐射率和地表温度等物理参数真实数值的物理数学计算。

1.辐射定标

ENVI提供Calibration Utilities(校正工程)工具,可以利用定标系数完成ASTER、MSS、TM、ETM+、QuickBird等传感器的辐射定标。本次实验介绍Landsat传感器定标的具体操作过程:

(1)在ENVI主菜单中,选择“File>Open External File>Iandsat>GeoTIFF with Metadata”,打开Landsat5 TM数据L5127043_04320051009 MTL.txt。

(2)在 ENVI主菜单中,选择“Basic Tools> Prepr coessing > Calibration Utilities >Landsat Calibration”,在打开的“Landsat Calibration Input File”对话框中,选择Land ast5 TM 文件,单击【OK】按钮,出现“ENVI Landsat Calibration”对话框(图19-1)。

(3)在“ENVI Landsat Calibration”对话框中,ENVI将自动从元数据中获取下列参数:Landsat卫星类型(Landsat Satellite Sensor);图像成像时间(Data Acquisition Month/Day/Year);太阳高度角(Sun Elevation)。

(4)定标类型(Calibration Type)包括“Radiance”(辐射亮度值)和“Reflectance”(反射率),选择“Radiance”(辐射亮度值)。

(5)点击【Edit Calibration Parameters】按钮,可以打开定标参数对话框,可以自己修改定标参数。

如果定标的数据格式是 ENVI标准或者TIF格式,需要手动输入“ENVI Landsat Calibration”对话框中的参数(Landsat卫星类型、图像成像时间和太阳高度角),并且每次只能定标一个波段。

图19-1 ENVI陆地卫星定标对话框

(6)选择输出路径及文件名,单击【OK】按钮,执行定标过程。

2.大气校正

目的是消除大气中的水蒸气、氧气、二氧化碳、甲烷和臭氧等对地物反射的影响,消除大气分子和气溶胶散射的影响,目前遥感图像的大气校正方法很多,本次实验选择FLAASH大气校正工具进行大气校正,FLAASH 对大气校正的输入图像作了一些要求,具体如下:

◎卫星图像波段范围为400~2500nm;航空图像为860~1135nm。

◎文件类型为ENVI标准栅格格式,BIP或BIL存储格式。

◎数据头文件中包含中心波长(Wavelength)值(表19-1),如果是高光谱还需要有波段宽度(FWHM),这两个参数都可以通过编辑头文件信息输入(Edit Header)。

表19-1 Ladnst中心波长

续表

图19-2 转换文件参数设置

FLAASH大气校正具体操作步骤如下:

(1)数据储存顺序调整。由于FLAASH大气校正工具要求数据存储顺序为BIP或者BIL格式,而ENVI默认的数据存储格式为BSQ格式,我们需要进行数据储存格式调整,在ENVI主菜单中,选择“Basic Tools>Convert Data(BSQ、BIP、BIL)”,在“Convert File Input”对话框中选择上一步辐射定标的结果,单击【OK】按钮打开“Convert File Parameters”对话框(图19-2)。

◎选择“Output Interleave”:BIL,

◎Convert In Place: Yes

单击【OK】按钮,执行处理。

(2)文件输入与输出信息。在ENVl主菜单中,选择“Spectral>FLAASH”,打开“FLAASH”功能,如图19-3所示。

点击【Input Radiance lmage】按钮,选择准备好的辐射亮度值数据,由于辐射定标时得到的辐射亮度单位为W/(m2 ·μm ·sr),FLA ASH要求的辐射亮度单位为μW/(cm2 nm ·sr),两者相差10倍,所以点击【Input Radiance Image】按钮会出现“Radiance Scale Factors”对话框(图19-4),选择“Us esingle scale factorfor all bands”(对所有波段使用统一比例因子),“single scale factor”值输入10,点击【OK】按钮,这样就获得符合FLAASH要求的辐射亮度值。

图19-3 FLAASH模块参数对话框

图19-4 辐射亮度比例因子设置

单击图19-3中的【Output Reflectance File】按钮,选择输出文件名和路径。

(3)传感器与图像目标信息。

Scene Center Location:从元数据文件中获取,当图像位于西半球时,经度为负值;位于南半球时,纬度为负值。

Sensor Type:选择辐射亮度图像对应的传感器类型,本实验使用数据为Landsat TM5。

Ground Elevation(km):可从相应区域的DEM 获取平均值。

Flight Date:图像成像日期,可以从元数据文件中获取。

Flight Time GMT:图像成像格林尼治时间,可以从元数据文件中获取。

(4)大气模型(Atmospheric Model)。ENVI提供了六种大气模型:亚极地冬季(Sub-Arctic Winter)、中纬度冬季(Mid-Latitude Winter)、美国标准大气模型(U.S.Standard)、亚极地夏季(Sub-Arctic Summer)、中纬度夏季(Mid-Latitude Summer)、热带(TroPical)。可以通过季节-纬度信息选择大气模型。

(5)水汽反演(Water Retrieval)。多光谱数据由于缺少相应波段和光谱分辨率太低不执行水汽反演。选择No。

(6)气溶胶模型(Aerosol Model)。ENVI提供五种气溶胶模型:无气溶胶(No Aerosol)、乡村(Rural)、城市(Urban)、海面(Maritime)、对流层(Tropospheric)。

(7)气溶胶反演(Aerosol Retrieval):选择2-Band(K-T), K-T气溶胶反演方法。

(8)初始能见度(Initial Visibility)(km):见表19-2。

表19-2 天气条件与估算能见度

(9)多光谱设置(MultispectraI Settings)。在“Sensor Type”选项中选择多光谱传感器时,出现【Multispectral Settings】按钮,单击此按钮可以打开多光谱设置对话框(图19-5),有两种设置方式:文件方式(File)和图形方式(GUI),一般选择图形方式。由于多光谱数据一般不用于水汽,因此多光谱设置对话框中主要的参数见“Kaufman Taner Aerosol Retrieval”对话框中,选择“Defaults”下拉框:Over-Land Retrieval Standard(660:2100),即上行通道660nm,下行通道2100nm,点击【OK】按钮后完成设置,回到FLAASH模块参数对话框。

图19-5 多光谱设置对话框

(10)高级设置(Advanced Settings)。按照默认设置。

上述设置完成后,点击【Apply】按钮,执行FLAASH大气校正。

3.大气校正前后同名点植被波谱曲线对比

(1)分别在“Display”窗口中显示原始数据图像和FLAASH 校正后的反射率图像。

(2)在其中一个主窗口中选择“Link>Geographic Link”,将两个“Display”窗口地理链接。

(3)将大气校正前后影像窗口移动到植被区域,分别在两个“Display”窗口中单击鼠标右键选择“Z Pro-file(Spectrmu)”,获取两个图像上的植被波谱曲线(图19-6),可以看到经过大气校正后的植被波谱曲线更加接近真实植被波谱曲线。

图19-6 大气校正前(左图)后(右图)植被波谱曲线对比

六、实验报告

(1)简述实验过程。

(2)回答问题:①什么是辐射定标?ENVl软件采用何种功能命令进行辐射定标?②什么是大气校正?ENVI软件采用何种功能命令进行大气校正?③辐射校正前后地物波谱特征对比分析:分别提取辐射校正前、后的漓江水体、岩溶石山、碎屑岩土山、农田、桂林市区和飞机场六种典型地物的波谱曲线。比较分析它们的差异。为何它们会存在差异?

实验报告格式见附录一。