辐射定标教程介绍.docVIP

在新版本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*micro;m))。定标参数使用Chander, Markham, and Helder (2009)的研究成果，其中LPGS 和NLAPS分别是两种数据处理系统得到的产品：the Level 1 Product Generation System (LPGS) 和 theNational Land Archive Production System (NLAPS)，从2008年12月份开始，L7 ETM+ 和 L5都是以LPGS系统处理，L4 TM和 MSS以NLAPS系统处理。具体参数如下所示：NLAPS处理和LPGS处理产品的MSS定标参数NLAPS处理产品的TM4、5定标参数说明：L5 TM在(1984–1991)间数据使用(LMAXλ=169, 333)定标参数。LPGS处理产品的ETM+数据定标参数同时ENVI还可以定标为表观大气反射率（ρp）：公式3 式中：Lλ为辐射亮度值d为天文单位的日地距离ESUNλ为太阳表观辐射率均值，对于Landsat 7 ETM+，ENVI使用“the Landsat 7 Science Data Users Handbook”上记录的参数；对于Landsat TM 4/5,ENVI使用Chander and Markham (2003)研究成果θs是以度为单位的太阳高度角。如果定标的Landsat数据带有元数据，ENVI会自动识别定标参数。对于带元数据文件的Landsat 7 GeoTIFF（*_MTL.txt） ，需要通过这种方式打开文件File → Open External File → Landsat → GeoTIFF with Metadata。选择Basic Tools → Preprocessing → Calibration Utilities → Landsat Calibration，可以打开定标工具。 图3 Fast格式的L7定标界面? ?? ? 图4 Landsat 7 GeoTIFF（*_MTL.txt）定标界面当遇到没有元数据文件的数据时，需要手动输入定标参数（默认情况下的参数经常不准），定标参数