地球信息科学学报 ›› 2014, Vol. 16 ›› Issue (1): 45-53.doi: 10.3724/SP.J.1047.2014.00045

• 地球信息科学理论与方法 • 上一篇    下一篇

改进的DisTrad模型在地形起伏区地表温度空间分辨率提升的应用

庞庆非, 权凌   

  1. 电子科技大学资源与环境学院, 成都 611731
  • 收稿日期:2013-04-02 修回日期:2013-05-20 出版日期:2014-01-05 发布日期:2014-01-05
  • 作者简介:庞庆非(1987-),男,山东聊城人,硕士生,研究方向为热红外遥感。E-mail:pcrs2006@gmail.com
  • 基金资助:

    国家自然科学基金项目(41101380);国家重点基础研究发展计划(2013CB733406);中国博士后科学基金面上项目(2011M500145);江苏省资源环境信息工程重点实验室(中国矿业大学)开放基金资助项目(JS201103)。

Disaggregation of MODIS Land Surface Temperature Using Stepwise Regression:A Case Study over Sichuan Basin

PANG Qingfei, QUAN Ling   

  1. School of Resources and Environment, University of Electronic Science and Technology of China, Chengdu 611731, China
  • Received:2013-04-02 Revised:2013-05-20 Online:2014-01-05 Published:2014-01-05

摘要:

DisTrad(Disaggregation procedure for radiometric surface temperature)模型是常用于遥感地表温度空间分辨率提升的主要模型之一。DisTrad模型常面向空间范围有限、地形相对平坦的研究区域,且常选用植被参数(如植被指数或植被覆盖度等)作为关键参数。然而在空间范围较大、地形起伏地区,地表温度的空间变异可能无法完全通过植被参数解释。本研究选取四川盆地及毗邻地区为研究区,通过模拟数据研究DisTrad模型在地形起伏区地表温度空间分辨率提升中的适用性。数字高程模型(Digital Elevation Model,DEM)、归一化差值植被指数(Normalized Difference Vegetation Index,NDVI)等参数,采用滑动窗口逐步回归,将空间分辨率为6km的地表温度提升至空间分辨率为1km。研究结果表明,改进的模型在平原及海拔较低的高原地区提升获得的地表温度空间分辨率具有较高精度,均方根误差(Root Mean Square Error,RMSE)为0.5K左右;在地形起伏较大的地区,RMSE为4K,验证了改进的模型提升的可行性。

关键词: 逐步回归, 滑动窗口, DisTrad模型, 地表参量, 地表温度

Abstract:

Accurate temporal and spatial estimation of land surface temperature (LST) is important for evaluating climate change, global hydrological cycle and monitoring urban heat islands (UHI). LSTs with high quality can be routine by using satellite remote sensing. However, characters of both high spatial and temporal resolutions have been difficult. Cloud cover further reduces the useable observations of surface conditions. Monthly LST product (MOD11C3) composited and averaged temperature values at 0.05 degree latitude/longitude grids (CMG) have coarse spatial resolution (~5.5 km). An alternative to the lack of high-resolution observations is to disaggregate LST data using other products of MODIS of 1 km observations. Historically, disaggregation of LST at high resolutions (1 km) has relied on vegetation index, e.g. NDVI (Normalized Difference Vegetation Index). However, this downscaling method is not adequate for areas encompass basin and upland. We applied Digital Elevation Model (DEM), NDVI, Enhanced Vegetation Index (EVI), Albedo, and slope to resolve this drawback by utilizing stepwise regression method with a moving window. The following is the algorithm. Land surface parameter (LSP) data are sampled to the coarser thermal resolution. A stepwise regression is performed between the monthly temperature product and sampled land surface parameters, then a function f (LSP) framed. The parameters of the regression function are applied to LSP data at high, target resolution. Coarse-scale residual field represent variability in temperature driven by other factors other than vegetation and DEM is added back into the high-resolution base map. So, we utilize LSP to sharpen original images. A reasonable rectangle box that making certain pixel be center is outlined for stepwise regression. Function is obtained by stepwise between LST and LSP. Loop and downscale the other pixels until image processed. Coefficients and intercept are saved as images. The disaggregation LST is achieved by substituting images at target resolution to function. The size of the box flowed over the image in this paper is 19 by 19. Stepwise disaggregation algorithm is applied to the resample MOD13A1 and DEM data. The fitting parameters vary with different window scenes. In contrast, the number of DEM entered function is much larger than NDVI. That indicated DEM is more significant than NDVI, EVI, albedo and slope in most fields of the study area, especially in mountain area. The RMSE of downscaling LST is 4.93K. Image sharpening is therefore not a replacement for high-resolution thermal imaging sensors. Nevertheless, in the absence of thermal imagery because of cloud in Sichuan, DisTrad seems to be able to enhance the resolution of MOD11C3 product.

Key words: DisTrad, moving window, stepwise, land surface temperature, land surface parameter