1、a和b就是我们常说的偏移和增益。二、星上定标有些卫星载有辐射定标源、定标光学系统，在成像时实时、连续的进行定标。三、场地定标场地定标指的是遥感器处于正常运行条件下，选择辐射定标场地。一般选择沙漠地区，它的光谱响应稳定，如利比亚沙漠用于定标AVHRR、北非沙漠定标SPOT影像、敦煌西戈壁沙漠定标CBERS影像、美国的白沙导弹靶场常用于高分辨率图像的定标。通过选择典型的均匀稳定目标，用精密仪器进行地面同步测量感器过顶时的大气环境参量和地物反射率，利用遥感方程，建立图像与实际地物间的数学关系，得到定标参数以完成精确的传感器定标。基本原理：在遥感器飞越辐射定标场地上空时，在定标场地选择偌干个像元区，测
2、量传感器对应的地物各波段光谱反射率和大气环境参量（大气气溶胶光学厚度，大气中水，臭氧含量等）等参量，并利用大气辐射传输模型等手段求解传感器入瞳处各光谱带的辐射亮度，最后确定它与传感器输出的数字量化值之间的数量关系，求解定标系数，并估算定标不确定性。主要方法有：反射率法：在卫星过顶时同步测量地面目标反射率因子和大气光学参量（如大气光学厚度、大气柱水汽含量等）然后利用大气辐射传输模型计算出遥感器入瞳处辐射亮度值。具有较高的精度。辐亮度法：采用经过严格光谱与辐射标定的辐射计，通过航空平台实现与卫星遥感器观测几何相似的同步测量，把机载辐射计测量的辐射度作为已知量，去标定飞行中遥感器的辐射量，从而实现卫
3、星的标定。最后辐射校正系数的误差以辐射计的定标误差为主。辐照度法：又称改进的反射率法，利用地面测量的向下漫射与总辐射度值来确定卫星遥感器高度的表观反射率，进而确定出遥感器入瞳处辐射亮度，。这种方法是使用解析近似方法来计算反射率，从而可大大缩减计算时间和计算复杂性。 因此，我们常用的定标参数，有使用实验室定标的结果（如高分辨率传感器QuickBird、WorldView-1等）；也有使用实验室定标与星上定标相结合的参数（如NOAA、MSS等）；由于设备老化，Landsat TM5的定标参数有用实验室定标的（2003年前），也有用经过场地定标的参数（2003年后）；也有些学者为了获得更精确定标数据
4、，使用场地定标的方法完成一些传感器数据的定标过程。3 ENVI下的传感器定标 ENVI支持很多数据的定标，包括ASTER、MODIS、AVHRR, MSS, TM ,QuickBird, WorldView-1,TIMS等。也可以根据定标参数利用BandMath工具很方便的完成定标。QuickBird, WorldView-1的定标比较简单，下面介绍在ENVI下的ASTER、MODIS、AVHRR, Landsat的定标过程。 3.1 ASTER与MODIS定标 在新版本ENVI中（4.5版本以及更新），没有单独设立ASTER和MODIS定标的工具。对于ASTER L1A/L1B和MODIS
5、02级数据，在打开数据时会自动完成对数据的定标。如图1所示打开ASTER L1B的结果，在波段列表中，自动读取各个波段的中心波长信息，并按照波段范围信息（VNIR、SWIR、TIR）分组波段。其中VNIR、SWIR自动定标为辐射亮度，单位是：W/m2/sr/m；TIR数据定标为大气表观温度值，单位：开尔文。打开其中一个数据，浏览像元值，可以看到已经定标为浮点型的辐射亮度值。下载 （95.57 KB）2009-10-28 20:36图1 ASTER L1B数据如图2为打开MODIS 02级1km数据，其中250米和500米的波段经过重采样为1km加入这个数据集中。ENVI根据各个波段的中心波长信
6、息定标为三个类型数据：反射率数据（Reflectance）、辐射亮度值数据（Radiance）和发射率数据（Emissive）。其中反射率和发射率为01无单位值，辐射亮度值单位是：W/m2/m/sr。下载 （88.52 KB）图2 MODIS 02级数据如果打开原始的ASTER和MODIS的DN值数据，可以在ENVI主菜单中选择File-preferences，切换到Miscellaneous面板，将Auto-Correct ASTER/MODIS项设置为NO。3.2 Landsat数据定标ENVI4.7版本改进了Landsat数据定标的功能，对于Landsat4/5数据可以手动选择以下两种定
7、标公式：下载 （12.45 KB）式中:QCAL为原始量化的DN值LMIN为QCAL = 0时的辐射亮度值LMAX为QCAL = QCALMAX时的辐射亮度值注：LMIN和LMAX的值取自Chander, Markham, and Helder （2009）的研究成果。QCALMIN是最小量化定标像素值（与LMIN类似）。取值如下：1：LPGS产品04 April 2004之后的NLAPS产品0：04 April 2004之前的NLAPS产品如果没有元数据信息，QCALMIN取默认值1（TM和ETM+）或者0 （MSS）。QCALMAX为最小量化定标像素值（与LMAX类似）。根据元数据信息取值
8、为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） 和 theNational Land Archive Production System （NLAPS），从2008年12月份开始，L7 ETM+ 和 L5都是
9、以LPGS系统处理，L4 TM和 MSS以NLAPS系统处理。具体参数如下所示：下载 （36.45 KB）NLAPS处理和LPGS处理产品的MSS定标参数下载 （46.78 KB）NLAPS处理产品的TM4、5定标参数说明：L5 TM在（19841991）间数据使用（LMAX=169, 333）定标参数。下载 （31.73 KB）LPGS处理产品的ETM+数据定标参数同时ENVI还可以定标为表观大气反射率（p）：下载 （3.85 KB）公式3 式中：L为辐射亮度值d为天文单位的日地距离ESUN为太阳表观辐射率均值，对于Landsat 7 ETM+，ENVI使用“the Landsat 7 Sc
10、ience 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 Cal
11、ibration，可以打开定标工具。下载 （70.81 KB）图3 Fast格式的L7定标界面下载 （68.32 KB）图4 Landsat 7 GeoTIFF（*_MTL.txt）定标界面当遇到没有元数据文件的数据时，需要手动输入定标参数（默认情况下的参数经常不准），定标参数可以通过NASA网站查询或者是一些研究成果，由于不同处理系统或者数据经过重采样后原始的DN会有所不同，因此定标参数也有很多套标准。如下网址为NASA公布的定标参数。http:/landsathandbook.gsfc.nasa.gov/handbook/handbook_htmls/chapter9/chapter9.h
12、tml注意的是如果定标的数据格式是ENVI标准格式或者TIF格式，一次只能定标一个波段。下载 （63.22 KB）图5 缺少元数据文件时候的定标界面下面表1和表2是Landsat5的其他版本定标参数供参考。表1 Landsat TM5飞行前辐射定标系数（实验室定标）BandGainoffset10.602-1.5221.170-2.8430.806-1.1740.815-1.5150.108-0.3760.0551.2470.057-0.15表2 Landsat TM 5辐射定标系数（Gyannesh Chander等（2003年）从1984年3月1号至2003年5月4号2003年5月5号之后
13、LMINLMAX152.100.602431193.00.762824296.811.175100365.01.442510204.300.805765264.01.039880206.200.814549221.00.87258827.190.10807830.20.1198821.237815.3030.05515814.380.05698016.50.065294-1.153.3 AVHRR数据定标ENVI提供的AVHRR定标工具可以对来自NOAA-6、-7、-8、-9、-10、-11、-12、 -14、 -15、-16、-17、-18和-19卫星的AVHRR数据进行定标。波段1和2被定
14、标为反射率，波段3、4和5被定标为亮温值（单位：开尔文）。其中NOAA-6、-7、-8、-9、-10、-11选择主模块-Basic Tools（或Spectral）-Preprocessing- Data-Specific Utilities-AVHRR-Calibrate NOAA 6/7/8/9/10/11。界面比较简单。在SST技术方面，ENVI使用多通道海面温度（Multi-Channel Sea Surface Temperature MCSST）计算算法（包括三种：Split-window、dual-window和triple-window），分别用于白天和夜间的数据。每种方法的计
15、算公式如下所示。1.NOAA-12,-14,-15：Day MCSST Split Ts = a0 + a1*band4 + a2* （band4 - band5） + a3* （band4 - band5） * （sec（） -1）Night MCSST SplitTs = a0 + a1*band4 + a2* （band4 - band5） + a3* （band4 - band5） * （sec（）-1）Night MCSST DualTs = a0 + a1*band4 + a2* （band3 - band4） + a3* （sec（） -1）Night MCSST TripleT
16、s = a0 + a1*band4 + a2* （band3 - band5） + a3* （sec（）-1）2.NOAA-16,-17,-18： Ts = a0 + a1*band4 + a2*band5 + a3* （band4 - band5） * （sec（） -1）Ts = a0 + a1*band4 + a2*band5 + a3* （band4 - band5） * （sec（） -1）Ts = a0 + a1*band3 + a2*band4 + a3*band5 + a4* （band3 - band5） * （sec（） -1）3.NOAA-19Ts = a0 + a1*b
17、and4 + a2（band4 - band5） + a3（band4 - band5）（sec（） -1）Ts = a0 + a1*band4 + a2（band3 - band4） + a3（sec（） -1）Ts = a0 + a1*band4 + a2（band3 - band5） + a3（band3 - band5）（sec（） -1）4.NOAA MetOp-A FRAC：Ts = a0 + a1*band4 + a2* （band4 - band5） + a3* （band4 - band5） * （sec（） -1）Ts = a0 + a1*band4 + a2* （band
18、3 - band5） + a3*（band3 - band5） * （sec（） -1）其中Sec（）是卫星高度角的正切值，Band3、4、5是定标后的亮温值（单位:开尔文），“a”参数值是根据漂流浮标和热带太平洋固定浮标数据，利用回归模型获得。不同的卫星对应的参数不一样，详细参见表36。表3 Day MCSST SplitSatellite a0 a1 a2 a3 NOAA-12 -263.006 0.963563 2.579211 0.242598 NOAA-14 -278.43 1.017342 2.139588 0.779706 NOAA-15 -261.029735 0.959456
19、 2.663579879 0.570613 NOAA-16 -273.77 3.301267 -2.30195 0.628966 NOAA-17 -271.206 0.992818 -2.49916 0.915103 NOAA-18 -280.43 1.02453 2.10044 0.0784059 NOAA-19 -278.74596 1.01922 1.72270 0.80263 MetOp-A FRAC -273.816 1.00255 2.39451 0.903773 表 4 Night MCSST Split-263.94 0.967077 2.384376 0.480788 -28
20、2.24 1.029088 2.275385 0.752567 -271.3969724 0.993892 2.7523466369 0.662999 -273.15 -2.53655 3.5316 0.753291 -276.59 1.01015 -2.5815 1.000541 -276.075 1.00841 2.23459 0.736946 -277.71304 1.01432 1.91798 0.72064 -277.447 1.01377 2.52362 1.03056 表5 Night MCSST Dual-279.846 1.031355 1.288548 2.265075 -
21、273.914 1.008751 1.409936 1.975581 -283.5117285 1.041037 1.5875819344 1.67743 -273.082 1.50825 1.00413 1.52452 -276.603 1.01805 1.49789 1.96181 -279.755 1.02958 0.0502887 1.78302 -276.61174 1.01873 1.47374 1.88560 -273.235 1.00711 1.49927 1.88373 表6 Night MCSST Triplea4 -271.971 1.000281 0.911173 1.
22、710028 -275.364 1.010037 0.920822 1.760411 -276.7558563 1.015354 1.0635723508 1.294955 -271.763 0.733532 1.01684 0.344182 -0.753123 -274.622 1.00903 0.913248 0.440015 -274.398 1.00820 0.841674 0.377061 -275.24563 1.01084 0.81643 0.43235 -273.044 1.00424 0.894349 0.508159 4 总结 对于一些新上天的传感器，实验室定标参数基本能满足应用需求。当对精度要求较高时候，往往需要通过场地定标的方式提高定标的结果。对于一些服役时间较长的传感器，也需要通过场地定标的方式来校对实验室定标结果，最典型的是Landsat系列卫星。 值得注意的是，为保证精度，定标时的数据源往往是初级别的，因为定标是对DN进行运算，如果经过重采样的数据其DN值会
