基于AVHRR和VIRR数据的改进型Becker“分裂窗”地表温度反演算法

A modified Becker’s split window approach for retrieving land surface temperature from the AVHRR and VIRR data

  • 摘要: 为了将基于NOAA-9/AVHRR数据提出的Becker和Li的“分裂窗”地表温度算法成功地应用于长序列NOAA/AVHRR和FY-3A/VIRR数据的地表温度反演,为气候变化研究提供长序列、高精度、高分辨率的地表温度数据集,从辐射传输方程出发,首先利用MODTRA 4.1模式模拟了多种地表和大气状态下的光谱辐亮度数据,并结合AVHRR和VIRR通道4、5的光谱响应函数建立了温度数据集(TS,T4,T5);然后,基于该数据集采用最小二乘法重新计算了Becker和Li算法中的各参数,提出了一个适用于NOAA/AVHRR和FY-3A/VIRR数据的改进型Becker和Li分裂窗地表温度反演算法;并利用改进型算法对2008年4月27日03时12分(世界时)观测的一景覆盖北京地区的NOAA17/AVHRR数据进行了地表温度的反演,将反演结果与日本东京大学提供的同地区、同时相的MODIS地表温度产品进行了对比分析。结果表明,两种地表温度产品的相关系数为0.88,均方根偏差(RMSD)为2.1 K;在两种地表温度差值图像的频率直方图上有69.6%的像元的值在±2 K之内,37%的像元的值在±1 K之内。

     

    Abstract: In order to successfully apply the Becker and Li’s split window approach, which was proposed based on the NOAA-9 AVHRR data, to estimate the Land Surface Temperature (LST) from the different AVHRRs and VIRR data and further to provide a high-precision, long-time, and high resolution LST dataset for climate change research, a modified Becker and Li’s split widow approach is developed based on the radiative transfer equation in this paper. To begin with, the MODTRAN 4.1 is used to generate the spectral radiance data under a variety of surface and atmosphere conditions. Then, the temperature dataset (TS,T4,T5) is built by convolving the spectral radiance data with the spectral response functions of channels 4 and 5 of the AVHRRs and VIRR. The parameters of the Becker and Li’s split window approach are subsequently recalculated based on the temperature dataset using the least squares regression method. As an example of validation, a single image of the AVHRR-17 over the Beijing area acquired at 03:12 UTC 27 April 2008 was used to retrieve the LST using the modified Becker and Li’s approach, and the comparison between the retrieved LST from AVHRR-17 and that from the MODIS, which was provided by the University of Tokyo in Japan, over the same region at the same time shows that the correlation coefficient is 0.88 and the root mean square deviation (RMSD) is 2.1 K. Furthermore, the frequency histogram of the LST difference image shows that about 69.6% and 37% pixels in the image have the values falling within the scope of ±2 K and ±1 K, respectively.

     

/

返回文章
返回