辐射定标教程 - 图文 下载本文

在新版本ENVI中(4.5版本以及更新),没有单独设立ASTER和MODIS定标的工具。对于

ASTER L1A/L1B和MODIS 02级数据,在打开数据时会自动完成对数据的定标。 如图1所示打开ASTER L1B的结果,在波段列表中,自动读取各个波段的中心波长信息,并按照波段范围信息(VNIR、SWIR、TIR)分组波段。其中VNIR、SWIR自动定标为辐射亮度,单位是:W/m2/sr/μm;TIR数据定标为大气表观温度值,单位:开尔文。打开其中一个数

据,浏览像元值,可以看到已经定标为浮点型的辐射亮度值。

图1 ASTER L1B数据

如图2为打开MODIS 02级1km数据,其中250米和500米的波段经过重采样为1km加入这个数据集中。ENVI根据各个波段的中心波长信息定标为三个类型数据:反射率数据(Reflectance)、辐射亮度值数据(Radiance)和发射率数据(Emissive)。其中反射率和发

射率为0~1无单位值,辐射亮度值单位是:W/m2/μm/sr。

图2 MODIS 02级数据

如果打开原始的ASTER和MODIS的DN值数据,可以在ENVI主菜单中选择File->preferences,

切换到Miscellaneous面板,将Auto-Correct ASTER/MODIS项设置为NO。

3.2 Landsat数据定标

ENVI4.7版本改进了Landsat数据定标的功能,对于Landsat4/5数据可以手动选择以下两种

定标公式:

式中:

?QCAL为原始量化的DN值 ?LMINλ为QCAL = 0时的辐射亮度值 ?LMAXλ为QCAL = QCALMAX时的辐射亮度值

注:LMINλ和LMAXλ的值取自Chander, Markham, and Helder (2009)的研究成果。

?QCALMIN是最小量化定标像素值(与LMINλ类似)。取值如下:

1:LPGS产品

1:04 April 2004之后的NLAPS产品

0:04 April 2004之前的NLAPS产品

注:如果没有元数据信息,QCALMIN取默认值1(TM和ETM+))或者0 (MSS)。 ?QCALMAX为最小量化定标像素值(与LMAXλ类似)。根据元数据信息取值为127, 254, 255。当缺少元数据时,QCALMAX 取默认值:255 (TM 和 ETM+)或者127 (MSS)。作为结果的辐

射亮度值的单位为(W/(cm2*sr*µm))。

定标参数使用Chander, Markham, and Helder (2009)的研究成果,其中LPGS 和NLAPS分别是

两种数据处理系统得到的产品:the Level 1 Product Generation System (LPGS) 和 the National Land Archive Production System (NLAPS),从2008年12月份开始,L7 ETM+ 和 L5都

是以LPGS系统处理,L4 TM和 MSS以NLAPS系统处理。具体参数如下所示:

NLAPS处理和LPGS处理产品的MSS定标参数