Landsat 5和IRS P6数据辐射定标说明.doc

Landsat 5和IRS P6数据辐射定标说明Landsat 5和IRS P6数据是两种应用很广泛的遥感数据，在利用两种数据来 进行定量遥感研究过程中，不少用户会对它们的辐射定标问题产生疑问，现就中 国科学院对地观测与数字地球科学中心（中国遥感卫星地面站）生产的 Landsat 5 TM 和 IRS P6 数据的辐射定标问题做一说明（针对系统级纠正的 FAST 格式产 品）1. 辐射亮度计算辐射亮度计算采用如下的通用计算公式：-(QCAL - QCAL ) + Lmin min1）L - L -厶.L max mtn 九 QCAL - QCALmax min其中，L）为光谱辐射亮度（w/（m2卩msr））, QCAL为像元的DN值，QCAL 和QCALm.n分别为像元可以取得的最大和最小DN值（对于中国科学院对地观 测与数字地球科学中心（中国遥感卫星地面站）生产的Landsat 5和IRS P6数据, QCAL . =0）, L 和 L .分别为 QCAL= QCAL 和 QCAL= QCAL .时的光 min max min max min谱辐射亮度值对于 Landsat 5 TM 数据，QCAL = 255max对于 IRS P6 LISS 3 和 LISS 4 数据，QCAL =255，对于 AWIFS 数据,maxQCAL = 1023max1.1 TM数据Lmax和Lmin获取max minTM数据头文件中“RAD GAINS/BIASES”之后的部分给出的是各波段的L 和L .（注意这里给出的并非GAIN和BIAS），单位为mw/（cm2sr），没有考 max min虑波段宽度，需要除以波段宽度将单位变为mw/（cm2pmsr），然后乘以10将辐射 亮度单位变为标准单位（w/（m2pmsr）），再参与有关的定量遥感计算。

从头文件计算TM数据各波段的辐射亮度要考虑波段宽度并进行单位换算， 比较复杂，另外中国科学院对地观测与数字地球科学中心（中国遥感卫星地面站） TM数据处理系统采用的辐射定标算法与USGS TM数据处理系统采用的辐射定 标算法接近，而采用 USGS 给出的定标参数来计算辐射亮度不用考虑波段宽度和 单位换算，比较方便和精确，因此建议直接利用USGS的定标参数来进行TM数 据的辐射定标［1］，具体如下：L^=gain*DN+bias （2）其中，L九为光谱辐射亮度（w/（m2pmsr）， gain和bias为各波段的增益和偏置 （w/（m2ymsr）USGS在2003年5月4日和2007年4月1日分别对TM数据辐射 定标算法进行了两次更新，gain和bias参数的变化情况如表1：表1 USGS各波段的gain和bias值处理1984年3月1日至2003年5月5日至2007年4月1日至今时间2003年5月4日2007年4月1日接收1984年3月1日至1984年3月1日至1984年3月1日至1992年1月1日至时间2003年5月4日2007年4月1日1991年12月31日今波段gainbiasgainbiasgainbiasgainbias10.602431-1.520.762824-1.520.668706-1.520.762824-1.5221.175100-2.841.442510-2.841.317020-2.841.442510-2.8430.805765-1.171.039880-1.171.039880-1.171.039880-1.1740.814549-1.510.872588-1.510.872588-1.510.872588-1.5150.108078-0.370.119882-0.370.119882-0.370.119882-0.3760.0551581.23780.0551581.23780.0551581.23780.0551581.237870.056980-0.150.065294-0.150.065294-0.150.065294-0.15中国科学院对地观测与数字地球科学中心(中国遥感卫星地面站)于 2004 年5月对Landsat 5 TM数据辐射定标算法进行了更新，这次更新与USGS 2003 年5月 4日的更新对应，因此中国科学院对地观测与数字地球科学中心(中国遥 感卫星地面站)2004年5月以后处理的1级产品应使用USGS 2003年5月给出 的定标参数计算辐射亮度。

Landsat 5 TM第6波段迄今一直采用基于星上内部定标灯的辐射定标算法， 定标系数一直没有变化从1级产品DN值计算辐射亮度L (w/m2pmsr)的公式 为：L=0.055158xDN+1.2378 (3)但是近年来研究发现1999年4月以后接收的TM数据第6波段的辐射亮度 被低估了 0.092w/m2pmsr,温度被低估了 0.68K(300K 时)[2]2 IRS P6数据Lmax和Lmn获取P6 数据头文件中 “REV CBIASES AND GAINS IN THE BAND ORDER ASON THIS TAPE”后给出的是各波段的L 和L i ,与TM数据不同的是，此处 max min的单位是mw/(cm2pmsr)，已经考虑了波段宽度，只需乘以10将单位转换为 w/(m2pmsr)就可以利用(1)式得到光谱辐射亮度2. 表观反射率计算在得到光谱辐射亮度后有时还需要将其转换为表观反射率，转换公式如下：p=冗今2 (4)E cos(0 )0zP为表观反射率，d为日地距离，0为太阳天顶角(与头文件中给出的太阳 z高度角互为余角)，E0是大气层外相应波长的太阳光谱辐照度，TM和P6数据各 波段的E0值见表2和表3。

表2 TM数据各波段E0值(w/(m2ym))⑶波段「E0值Band11957Band21826Band31554Band41036Band5215.0Band780.67表3 P6数据各波段E0值(w/(m2卩m))[4]传感器Band2Band3Band4Band5LISS-31849.51553.01092.0239.52LISS-41853.61581.61114.3AWIFS1854.71556.41082.4239.84参考文献[1] 王文,李子扬，唐伶俐.中国遥感卫星地面站Landsat 5 TM的辐射定标.中国遥感卫星地面 站用户简讯,2004,(4):3.[2][3] Gyanesh Chander, Brian Markham. Revised Landsat 5 TM radiometric calibration procedures and postcalibration dynamic ranges[J]. IEEE Transactions on Geoscience and remote sensing, 2003, 41(11): 2674-2677.[4]。