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

doi:10.3724/SP.J.1047.2014.00045

摘要/Abstract

摘要：

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

庞庆非, 权凌. 改进的DisTrad模型在地形起伏区地表温度空间分辨率提升的应用[J]. 地球信息科学学报, 2014, 16(1): 45-53.DOI:10.3724/SP.J.1047.2014.00045

PANG Qingfei, QUAN Ling. Disaggregation of MODIS Land Surface Temperature Using Stepwise Regression：A Case Study over Sichuan Basin[J]. , 2014, 16(1): 45-53.DOI:10.3724/SP.J.1047.2014.00045

参考文献

[1] Kustas W P, Norman J M, Anderson M C, et al. Estimating subpixel surface temperature and energy fluxes from the vegetation index-radiometric temperature relationship[J]. Remote Sensing of Environment, 2003(85):429-440.

[2] Agam N, Kustas W P, Anderson M, et al. A vegetation index based technique for spatial sharpening of thermal imagery[J]. Remote Sensing of Environment, 2007(107):545-558.

[3] Zakšek K, Oštir K. Downscaling land surface temperature for urban heat island diurnal cycle analysis[J]. Remote Sensing of Environment, 2012(117):114-124.

[4] 吕京国, 张小咏, 蒋玲梅, 等.MODIS地表产品数据的相关算法及处理过程[J].遥感应用, 2009(4):25-29.

[5] 梁顺林著, 范闻捷, 等译.定量遥感[M].北京:科学出版社, 2009, 240.

[6] Hais M, Kucera T. The influence to topography on the forest surface temperature retrieved from Landsat TM, ETM+ and ASTER thermal channels[J]. Journal of Photogrammetry and Remote Sensing, 2009(64):585-591.

[7] Markus N. Estimating daily land surface temperature in mountainous environments by reconstructed MODIS LST data[J]. Remote Sensing, 2010(2):333-351.

[8] 李粉玲, 李京忠, 张琦翔.DEM提取坡度坡向算法的对比研究[J].安徽农业科学, 2008, 36(17):7355-7357.

[9] 卢毅敏, 岳天祥, 陈传法, 等.中国太阳总辐射的多元逐步回归模拟[J].遥感学报, 2010, 14(5):858-864.

[10] Land Processes DAAC USGS Earth Resources Observation and Science(EROS) Center. MODIS Reprojection Tool User's Manual. 2008. Release 4.0.

[11] Agam N, Kustas W P, Anderson M C, et al. Utility of thermal image sharpening for monitoring field-scale evapotranspiration over rain-fed and irrigated agricultural regions[J]. Geophysical Research Letters, 2008, 35(L02402):10.1029/2007GL032195.

[12] 杨静学, 苏华, 王云鹏. DisTrad热像元分解模型运用于高植被覆盖区的问题及改进[J].遥感技术与应用, 2010, 25(3):346-352.

[13] Agam N, Kustas W P, Anderson M C, et al. Utility of thermal sharpening over Texas high plains irrigated agricultural fields[J]. Journal of Geophysical Research, 2007, 112(19):1-10.

[14] Zhan W F, Chen Y H, Zhou J, et al. Sharpening thermal imageries: a generalized theoretical framework from an assimilation perspective[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(2):773-789.

[15] Nichol J. An emissivity modulation method for spatial enhancement of thermal satellite images in urban heat island analysis[J]. Photogrammetric Engineering and Remote Sensing, 2009(75):1-10.

[16] Inamdar A K, French A. Disaggregation of GOES land surface temperatures using surface emissivity[J]. Geophysical Research Letters, 2009, 36(2):1-5.

[17] Mao K, Qin Z, Shi J, et al. A practical split-window algorithm for retrieving land-surface temperature from MODIS data[J]. International Journal of Remote Sensing, 2005, 26(15):3181-3204.

[18] Jimenez-Munoz J C, Sobrino J A. A generalized single-channel method for retrieving land surface temperature from remote sensing data[J]. Journal of Geophysical Research, 2003, 108(22):1-9.

[19] Zhou J, Li J, Zhang L X, et al. Intercomparison of methods for estimating land surface temperature from a Landsat-5 TM image in an arid region with low water vapour in the atmosphere[J]. International Journal of Remote Sensing, 2012, 33(3):2582-2602.

[20] Inamdar A K, French A, Hook S, et al. Land surface temperature retrieval at high spatial and temporal resolutions over the sourthwestern United States[J]. Journal of Geophysical Research, 2008, 113(7):1-18.

[21] Menenti M, Bastiaanssen W, van Eick D, et al. Linear relationships between surface reflectance and temperature and their application to map actual evaporation of groundwater[J]. Advances in Space Research, 1989, 9(1):165-176.

[22] Merlin O, Duchemin B, Hagolle O, et al. Disaggregation of MODIS surface temperature over an agricultural area using a time series of Formosat-2 images[J]. Remote Sensing of Environment, 2010, 114(11):2500-2512.

[23] Merlin O, Chehbouni A, Walker J P, et al. A simple method to disaggregate passive microwave based soil moisture[J]. IEEE Transactions Geoscience Remote Sensing (SMOS Special Issue), 2008, 46(3):786-796.

[24] Merlin O, Bitar A A, Walker J P, et al. A sequential model for disaggregating near-surface soil moisture observations using multi-resolution thermal sensors[J]. Remote Sensing of Environment, 2009(113):2275-2284.

[25] Price J C. Estimating surface temperatures from satellite thermal infra-red data: A simple formulation for the atmospheric effect[J]. Remote Sensing of Environment, 1983(13):353-361.

[26] Jia S F, Zhu W B, Lu A F, et al. A statistical spatial downscaling algorithm of TRMM precipitation based on NDVI and DEM in the Qaidam Basin of China[J]. Remote Sensing of Environment, 2011, 115(12):3069-3079.

[27] Zhan W F, Chen Y H, Wang J F, et al. Downscaling land surface temperatures with multi-spectral and multi-resolution images[J]. International Journal of Applied Earth Observation and Geoinformation, 2012(18):23-36.

[28] Zhan W F, Chen Y H, Zhou J, et al. Disaggregation of remotely sensed land surface temperature: Literature survey, taxonomy, issues, and caveats[J]. Remote Sensing of Environment, 2013(131):119-139.