定量高光谱遥感

定量/高光谱遥感之传感器定标我们常用影像的像元值大多是经过量化的、无量纲的DN值，而进行遥感定量化分析时，常用到辐射亮度值、反射率值、温度值等物理量。传感器定标就是要获得这些物理量的过程。本专题叙述的主要是卫星传感器的定标，航空传感器有很多相似地方。本专题包括以下内容：l传感器定标概念l传感器定标类型lENVI下的传感器定标1传感器定标概念传感器定标很多地方又名为辐射定标，严格意义上讲，辐射定标是传感器定标的一部分内容。以下是国内的定义，如赵英时等《遥感应用分析原理与方法》上描述：定标是将遥感器所得的测量值变换为绝对亮度或变换为与地表反射率、表面温度等物理量有关的相对值的处理过程。或者说，遥感器定标就是建立遥感器每个探测器输出值与该探测器对应的实际地物辐射亮度之间的定量关系；建立遥感传感器的数字量化输出值DN与其所对应视场中辐射亮度值之间的定量关系（陈述彭）。辐射亮度的典型的单位为：W/cm2.以m.s（瓦特/平方厘米.微米.球面度）。我们总结以上的定义，通俗的说法：传感器定标就是将图像的数字量化值（DN）转化为辐射亮度值或者反射率或者表面温度等物理量的处理过程。其中反射率又分为大气外层表观反射率和地表实际反射率，后者又属于大气校正的范畴，有的时候也会将大气校正纳入传感器定标的一种途径。2传感器定标类型（定标参数测量）传感器定标可分为绝对定标和相对定标。绝对定标是获取图像上目标物的绝对辐射值等物理量；相对定标是将图像目标物辐射量归一化某个值范围内，比如以其他数据作为基准。传感器定标可分为三个阶段或者说三个方面内容：①发射前的实验室定标；②基于星载定标器的星上定标；③发射后的定标（场地定标）。一、实验室定标在遥感器发射之前对其进行的波长位置、辐射精度、光谱特性等进行精确测量，也就是实验室定标。它一般包含两部分内容：l光谱定标确定遥感传感器每个波段的中心波长和带宽，以及光谱响应函数。l辐射定标在模拟太空环境的实验室中，建立传感器输出的量化值（DN）与传感器入瞳处的辐射亮度之间的模型，一般用线性模型表示。Y=aX+ba和b就是我们常说的偏移和增益。也有些传感器使用其他表达式，如环境小卫星中的元数据说明文件中的这段话明确说明了定标公式、参数值、定标单位信息：（gain1,Fielddata,L=DN/g+L0,W*mA（-2）*srA（-1）*umA（-1））B1:g0.5782,L03.4608,B2:g0.5087,L05.8769,B3:g0.6825,L08.0069,B4:g0.6468,L08.8583二、星上定标有些卫星载有辐射定标源、定标光学系统，在成像时实时、连续的进行定标。三、场地定标场地定标指的是遥感器处于正常运行条件下，选择辐射定标场地。一般选择沙漠地区，它的光谱响应稳定，如利比亚沙漠用于定标AVHRR、北非沙漠定标SPOT影像、敦煌西戈壁沙漠定标CBERS影像、美国的白沙导弹靶场常用于高分辨率图像的定标。通过选择典型的均匀稳定目标，用精密仪器进行地面同步测量感器过顶时的大气环境参量和地物反射率，利用遥感方程，建立图像与实际地物间的数学关系，得到定标参数以完成精确的传感器定标。基本原理：在遥感器飞越辐射定标场地上空时，在定标场地选择偌干个像元区，测量传感器对应的地物各波段光谱反射率和大气环境参量（大气气溶胶光学厚度，大气中水，臭氧含量等）等参量，并利用大气辐射传输模型等手段求解传感器入瞳处各光谱带的辐射亮度，最后确定它与传感器输出的数字量化值之间的数量关系，求解定标系数，并估算定标不确定性。主要方法有：反射率法：在卫星过顶时同步测量地面目标反射率因子和大气光学参量（如大气光学厚度、大气柱水汽含量等）然后利用大气辐射传输模型计算出遥感器入瞳处辐射亮度值。具有较高的精度。辐亮度法：采用经过严格光谱与辐射标定的辐射计，通过航空平台实现与卫星遥感器观测几何相似的同步测量，把机载辐射计测量的辐射度作为已知量，去标定飞行中遥感器的辐射量，从而实现卫星的标定。最后辐射校正系数的误差以辐射计的定标误差为主。辐照度法：又称改进的反射率法，利用地面测量的向下漫射与总辐射度值来确定卫星遥感器高度的表观反射率，进而确定出遥感器入瞳处辐射亮度，。这种方法是使用解析近似方法来计算反射率，从而可大大缩减计算时间和计算复杂性。因此，我们常用的定标参数，有使用实验室定标的结果（如高分辨率传感器QuickBird、WorldView-1等）；也有使用实验室定标与星上定标相结合的参数（如NOAA、MSS等）；由于设备老化，LandsatTM5的定标参数有用实验室定标的（2003年前），也有用经过场地定标的参数（2003年后）；也有些学者为了获得更精确定标数据，使用场地定标的方法完成一些传感器数据的定标过程。定标实现上述传感器定标的方法很多，看起来也非常复杂，其实我们对遥感影像进行定标时候，只需要知道这类影像的定标公式和参数即可，这些公式和参数一般在数据带的元数据文件中，或者固定的公式和参数。现在很多软件提供的特定传感器的定标工具，如果没有也可以使用波段运算的方式实现定标。ENVI下的传感器定标ENVI支持很多数据的定标，包括ASTER、MODIS、AVHRR,MSS,TM,QuickBird,WorldView-1,TIMS等。也可以根据定标参数利用BandMath工具很方便的完成定标。QuickBird,WorldView-1的定标比较简单。或者使用Applygaineandoffset工具可以实现Y=aX+b公式的定标。下面介绍在ENVI下的ASTER、MODIS、AVHRR,Landsat的定标过程。3.1ASTER与MODIS定标对于ASTERL1A/L1B和MODIS02级数据，在打开数据时会自动完成对数据的定标。如图1所示打开ASTERL1B的结果，在波段列表中，自动读取各个波段的中心波长信息，并按照波段范围信息（VNIR、SWIR、TIR）分组波段。其中VNIR、SWIR自动定标为辐射亮度，单位是：W/m2/sr/〃mTIR数据定标为大气表观温度值，单位：开尔文。打开其中一个数据，浏览像元值，可以看到已经定标为浮点型的辐射亮度值。

图1：ASTERL1B数据如图2为打开MODIS02级1km数据，其中250米和500米的波段经过重采样为1km加入这个数据集中。ENVI根据各个波段的中心波长信息定标为三个类型数据：反射率数据（Reflectance）、辐射亮度值数据（Radiance）和发射率数据（Emissive）。其中反射率和发射率为0~1无单位值，辐射亮度值单位是：W/m2/〃m/s。

图2MODIS02级数据如果打开原始的ASTER和MODIS的DN值数据，可以在ENVI主菜单中选择File->preferences,切换到Miscellaneous面板，将Auto-CorrectASTER/MODIS项设置为NO。3.2Landsat数据定标对于Landsat4/5数据可以手动选择以下两种定标公式：Lb=Gain+Bias公式以=LM啊+J掀槌in)(QCAL-QCALMIN)公式2式中：QCAL为原始量化的DN值LMIN入为QCAL=0时的辐射亮度值LMAX入为QCAL=QCALMAX时的辐射亮度值注：LMIN入和LMAX入的值取自Chander,Markham,andHelder(2009)的研究成果。QCALMIN是最小量化定标像素值(与LMIN入类似)。取值如下：1：LPGS产品1：04April2004之后的NLAPS产品0:04April2004之前的NLAPS产品

注：如果没有元数据信息，QCALMIN取默认值1(TM和ETM+))或者0(MSS)。YQCALMAX为最小量化定标像素值(与LMAXX类似)。根据元数据信息取值为127,254,255。当缺少元数据时，QCALMAX取默认值：255(TM和ETM+)或者127(MSS)。作为结果的辐射亮度值的单位为(W/(cm2*sr*pm))。定标参数使用Chander,Markham,andHelder(2009)的研究成果，其中LPGS和NLAPS分别是两种数据处理系统得到的产品：theLevel1ProductGenerationSystem(LPGS)和theNationalLandArchiveProductionSystem(NLAPS),从2008年12月份开始，L7ETM+和L5都是以LPGS系统处理，L4TM和MSS以NLAPS系统处理。具体参数如下所示：图3：NLAPS处理和LPGS处理产品的MSS定标参数rpL5TMILPGS)OLg口侑睡财24口如OJ76-U.9D^1JW1*[g…g2-097-心丁MSIrpL5TMILPGS)OLg口侑睡财24口如OJ76-U.9D^1JW1*[g…g2-097-心丁MSI-05)8O.5K-0.W?Ub24-。.愤』叫而稣U硒”的I0J2-IL佛Zti^7-3tJ47OU452-0JUH.trfKICl砒G-CI.御与QL7?6-ig；L新■1.7K410.43-I?.4?掀17—WH财罪ftM9aH4ii碗sifttfi2222(1铀OU56SQ8矿as4i1..S76ll.tMG122J财部awiaG»eiS4Cii尚祁11435-132kSJ.10曲心相T52mi新息11.17^09S2.W1TO3-IJ7204.500181^^5-IJ7E9・1§】曜X汹-1511皿即37.19anwm«?期3Egw-0.15aosvyau415iiA-ly■质UM77C7■217IW3*1?J3L4心竣h.1.as加泪门冲陀I.DOWOhOS了皿4&125OT0邮5抓(1缸闻-Z.I?・K1539

脉-1521«9•2J91W101765827也口如*-4Iti4J7泌LMJSrTft-211mo-15]221篮亍湖4)033-0-1/3A20ll2l>tS4・心5133032心瓠的方「旧一*|了「!.了门/日4'西巾tUiklL5耕心州电QK山iUML^LMAXtE蹈所UfMltSW；i.wWheiSIVAill2口Hfl)Wl；m-m)图4：NLAPS处理产品的TM4、5定标参数说明：L5TM在(1984-1991)间数据使用(LMAXX=169,333)定标参数。L7ETMt1占皿5Qqzi,■*25%}由M^pK-train缨CenterLM1%原2srUTl)w:irunjEr.:ir■:inr.wiflw甜况u心j<31.452-0514・氏2293.7-TJ8139(72机蝴一6山3D0J，如37茁E3O£3l0U&92G56t-5J0一5甲153?40772-OJi%W5-5.1241.1g宓i5LM7-1748■IjG47.S70.191230■L19必86心1723$IZ04006响-0,07N；A7次掘ISM-035■CL岖S4.WtKi03T5-0SMCLXS6-4.7Mil。.酣5油-5.68g时SOffiS}。或皿Mum-b.21916n.7fay«一g§勒2商斗0641Q攻-e.419&5Q.74KSI<*-74Dg3062啪90痂垢541S3J40.772^0J9eOJ35-51117.4g的础5l.547«1.743g-w碰si如m一"3&JM5Mg体N/A7*3.锵-03510MI找i知CL*-4.715*3炀引732-134图5：LPGS处理产品的ETM+数据定标参数同时ENVI还可以定标为表观大气反射率(pp)：ITT*L/d‘ESUNa*cos0sPp=式中：ESUNa*cos0sL入为辐射亮度值d为天文单位的日地距离ESUN入为太阳表观辐射率均值，对于Landsat7ETM+，ENVI使用“theLandsat7ScienceDataUsersHandbook”上记录的参数;对于LandsatTM4/5,ENVI使用ChanderandMarkham(2003)研究成果Os是以度为单位的太阳高度角。如果定标的Landsat数据带有元数据，ENVI会自动识别定标参数。对于带元数据文件的Landsat7GeoTIFF(*_MTL.txt)，需要通过这种方式打开文件File—OpenExternalFile—Landsat—GeoTIFFwithMetadata。选择BasicTools—Preprocessing—CalibrationUtilities—LandsatCalibration，可以打开定标工具。

图6：Landsat5GeoTIFF（七MTL.txt）定标界面当遇到没有元数据文件的数据时，需要手动输入定标参数，定标参数可以通过NASA网站查询或者是一些研究成果，由于不同处理系统或者数据经过重采样后原始的DN会有所不同，因此定标参数也有很多套标准。如下网址为NASA公布的定标参数。/handbook/handbook_htmls/chapter9/chapter9.html下面表1和表2是Landsat5的其他版本定标参数供参考。表1LandsatTM4/5飞行前辐射定标系数（实验室定标）BandGainoffset10.602-1.5221.170-2.8430.806-1.1740.815-1.5150.108-0.3760.0551.2470.057-0.15表2LandsatTM5辐射定标系数（GyanneshChander等（2003年））从1984年3月1号至2003年5月4号2003年5月5号之后BandLMINLMAXGainoffsetLMINLMAXGainoffset1-1.52152.100.602431-1.52-1.52193.00.762824-1.522-2.84296.811.175100-2.84-2.84365.01.442510-2.843-1.17204.300.805765-1.17-1.17264.01.039880-1.174-1.51206.200.814549-1.51-1.51221.00.872588-1.515-0.3727.190.108078-0.37-0.3730.20.119882-0.3761.237815.3030.0551581.23781.237815.3030.0551581.23787-0.1514.380.056980-0.15-0.1516.50.065294-1.153.3AVHRR数据定标ENVI提供的AVHRR定标工具可以对来自NOAA-6、-7、-8、-9、-10、-11、-12、-14、-15、-16、-17、-18和-19卫星的AVHRR数据进行定标。波段1和2被定标为反射率，波段3、4和5被定标为亮温值(单位：开尔文)。其中NOAA-6、-7、-8、-9、-10、-11选择主模块->BasicTools(或Spectral)->Preprocessing->Data-SpecificUtilities->AVHRR->CalibrateNOAA6/7/8/9/10/11o界面比较简单。在SST技术方面，ENVI使用多通道海面温度(Multi-ChannelSeaSurfaceTemperatureMCSST)计算算法(包括三种：Split-window、dual-window和triple-window),分别用于白天和夜间的数据。每种方法的计算公式如下所示。NOAA-12,-14,-15：lDayMCSSTSplitTs=a0+a1*band4+a2*(band4-band5)+a3*(band4-band5)*(sec(速)-1)lNightMCSSTSplitTs=a0+a1*band4+a2*(band4-band5)+a3*(band4-band5)*(sec(p)-1)lNightMCSSTDualTs=a0+a1*band4+a2*(band3-band4)+a3*(sec(速)-1)lNightMCSSTTripleTs=a0+a1*band4+a2*(band3-band5)+a3*(sec(魂)-1)NOAA-16,-17,-18：lDayMCSSTSplitTs=a0+a1*band4+a2*band5+a3*(band4-band5)*(sec(p)-1)lNightMCSSTSplitTs=a0+a1*band4+a2*band5+a3*(band4-band5)*(sec(p)-1)lNightMCSSTDualTs=a0+a1*band4+a2*(band3-band4)+a3*(sec(速)-1)lNightMCSSTTripleTs=a0+a1*band3+a2*band4+a3*band5+a4*(band3-band5)*(sec(p)-1)NOAA-19lDayMCSSTSplitTs=a0+a1*band4+a2(band4-band5)+a3(band4-band5)(sec(p)-1)lNightMCSSTSplitTs=a0+a1*band4+a2(band4-band5)+a3(band4-band5)(sec(p)-1)lNightMCSSTDualTs=a0+a1*band4+a2(band3-band4)+a3(sec(p)-1)lNightMCSSTTripleTs=a0+a1*band4+a2(band3-band5)+a3(band3-tand5)(sec(p)-1)NOAAMetOp-AFRAC：lDayMCSSTSplit

Ts=a0+a1*band4+a2*(band4-band5)+a3*(band4-band5)*(sec(速)-1)lNightMCSSTSplitTs=a0+a1*band4+a2*(band4-band5)+a3*(band4-band5)*(sec(速)-1)lNightMCSSTDualTs=a0+a1*band4+a2*(band3-band4)+a3*(sec(速)-1)lNightMCSSTTripleTs=a0+a1*band4+a2*(band3-band5)+a3*(band3-band5)*(sec(速)-1)其中Sec(?)是卫星高度角的正切值，Band3、4、5是定标后的亮温值(单位:开尔文)，“a”参数值是根据漂流浮标和热带太平洋固定浮标数据，利用回归模型获得。不同的卫星对应的参数不一样，详细参见表3~6。表3DayMCSSTSplit表4表4NightMCSSTSplitSatellitea0a1a2a3NOAA-12-263.940.9670772.3843760.480788NOAA-14-282.241.0290882.2753850.752567NOAA-15-271.39697240.9938922.75234663690.662999NOAA-16-273.15-2.536553.53160.753291NOAA-17-276.591.01015-2.58151.000541NOAA-18-276.0751.008412.234590.736946NOAA-19-277.713041.014321.917980.72064MetOp-AFRAC-277.4471.013772.523621.03056Satellitea0a1a2a3NOAA-12-263.0060.9635632.5792110.242598NOAA-14-278.431.0173422.1395880.779706NOAA-15-261.0297350.9594562.6635798790.570613NOAA-16-273.773.301267-2.301950.628966NOAA-17-271.2060.992818-2.499160.915103NOAA-18-280.431.024532.100440.0784059NOAA-19-278.745961.019221.722700.80263MetOp-A