Land Surface Temperature Data Fusion of Landsat ETM and MODISby Combining TsHARP and STITFM Model

Land Surface Temperature (LST) is an important parameter that describes energy balance of substance and energy exchange between the surface and the atmosphere,and LST has widely used in the fields of urban heat island effect,soil moisture and surface radiative flux.Currently,no satellite sensor can deliver thermal infrared data at both high temporal resolution and spatial resolution,which strongly limits the wide application of thermal infrared data.Based on the MODIS land surface temperature product and Landsat ETM+image,a temporal and spatial fusion method is proposed by combining the TsHARP (Thermal sHARPening) model with the STITFM (Spatio\|Temporal Integrated Temperature Fusion Model) algorithm,defined as CTsSTITFM model in this study.The TsHARP method is used to downscale the 1 km MODIS land surface temperature image to LST data at spatial resolution of 250 m.Then the accuracy is verified by the retrieval LST from Landsat ETM+ image at the same time.Land surface temperature image at 30 m spatial scale is predicted by fusing Landsat ETM+ and downscaling MODIS data using STITFM model.The fusion LST image is validated by the estimated LST from Landsat ETM+ data for the same predicted.The results show that the proposed method has a better precision comparing to the STITFM algorithm.Under the default parameter setting,the predicted LST values using CTsSTITFM fusion method have a root mean square error (RMSE) less than 1.33 K.By adjusting the window size of CTsSTITFM fusion method,the fusion results in the selected areas show some regularity with the increasing of the window.In general,a reasonable window size set may slightly improve the effects of LST fusion.The CTsSTITFM fusion method can solve the problem of mixed pixels caused by coarse\|scale MODIS surface temperature images to some degree.     

针对TM6有效波长的计算方法及其对地温反演的影响研究

相对NOAA-AVHRR及MODIS的热波段而言,Landsat TM6有更高的空间分辨率;因而更适合于地面热力场分布特征及微环境研究的准确分析需要。有效波长的确定对反演地表温度有重要影响。本文根据Landsat5及Landsat7 TM6对波长的光谱响应函数(RSR),推导出相应的有效波长值并分析其对地温反演的影响。分析结果表明,一般情况下,就Landsat5而言,用中心波长反演地表温度不会有太大影响;但利用Landsat7反演地表温度时则必须使用有效波长,因为中心波长可能会导致较大的地表温度反演误差。     

Algorithms Comparison of Land Surface Temperature Retrieval from Landsat Series Data：A Case Study in Qiqihar，China

Land Surface Temperature（LST)is considered to be one of the significant indicators of urban environment analysis.Landsat thermal infrared series data is an important data source for retrieving surface temperature.In this paper，the thermal infrared band of the Landsat data in 2002，2008 and 2016 were used to retrieve LST by three different algorithms in municipal area of Qiqihar，China.These algorithms were the Mono-Window algorithm（MW algorithm)，the Single Channel algorithm（SC algorithm) and the Radiation Transport Equation method（RTE algorithm).And the results of the retrieval were compared to each other and verified by MODIS surface temperature products.The LST distribution maps were accomplished according to the retrieval results.The results showed that：(1)The spatial distribution of the LST obtained by the retrieval of the Landsat series by the three algorithms is consistent，and the LSTof the urban center is higher and thetemperature of water is the lowest；(2)Based on ETM+ data，the consistency between SC and RTE algorithm results is good，among which the SC algorithm has the highest precision，and the MW algorithm has large errors in different land cover areas；(3)The retrieval results by MW algorithm based on the TM data has the highest accuracy，RTE algorithm results is second，and the LST form SC algorithm is less consistent with the corresponding MODIS temperature products；(4)Based on the Landsat 8 TIRS data，the SC algorithm has the highest accuracy and the RTE algorithm has a large error.     

Remote sensing of temperature variations around major power plants as point sources of heat

Variations in land surface temperature (LST) around major point sources of heat were studied using the Tampa Bay region as a case study. LST in the Tampa Bay region, FL, USA, was retrieved from Landsat Thematic Mapper (TM) 6 and Enhanced Thematic Mapper Plus (ETM+) 6 high-gain thermal bands. The TM6 image data were obtained on 29 January (winter season) and 3 April 2007 (spring season). The ETM+6 data were obtained on 11 April 2007 (spring season). Spatial profiles of LST around four major fossil-fuelled power plants (FFPPs) were considered in this study. Temperatures were found to be highest at power plants and to decay to an average background temperature within 1.2–2.0 km from the FFPPs. The average background temperatures obtained in January and April were 17°C and 29°C, respectively. Results indicate that LST in close proximity to the FFPPs could be up to 10°C hotter than the surrounding areas. These findings suggest that FFPPs are significant heat sources and populations living within 1–2 km from an FFPP might be at significantly higher risk of heat-related illnesses and mortality.     

使用单窗算法研究北京城区热岛效应

随着全球变暖和城市化进程的加快,大城市城区的热岛效应日益严重。城市下垫面对地表能量交换的影响巨大,引起地表温度分布的不均一性。遥感技术的发展为地表温度的反演提供了可能。近年来人们使用劈窗算法对均一的海面温度的反演很成功,但是受空间分辨率的限制以及陆面的不均一性,陆面温度的反演一直是一个没有解决好的问题。覃志豪提出了一种TM热红外波段单窗算法,可以利用辅助气象资料快速计算出地表温度。本文以北京市城区为研究区,采用LandsetETM第6波段的单窗算法,反演了亮度温度和地表实际温度,分析了城市下垫面情况下NDVI与地表温度的相关关系,并解释了北京城区热岛在空间上的分布及其可能的原因。结果表明:北京市城区热岛效应显著;地表温度与NDVI相关性显著;城区绿地和水体在区域的温度分布中起到重要作用。     

Change analysis of land surface temperature based on robust statistics in the estuarine area of Pearl River (China) from 1990 to 2000 by Landsat TM/ETM+ data

The estuarine area of Pearl River that has taken great changes in land cover since 1990 is a typical area for studying the change of land surface temperature (LST). The LST of the years 1990 and 2000 in this area was estimated from the data of Landsat TM/ETM+ band 6, respectively, and three scales, corresponding to high, normal and low temperature ranges, were divided by a robust statistical method. The results show that the area of high temperature range in 2000 has an increase of 250 km² compared with the year 1990. The urban‐used land and the bare land are the main land cover types constituting the high temperature range area.     

A mono-window algorithm for retrieving land surface temperature from Landsat TM data and its application to the Israel-Egypt border region

Remote sensing of land surface temperature (LST) from the thermal band data of Landsat Thematic Mapper (TM) still remains unused in comparison with the extensive studies of its visible and near-infrared (NIR) bands for various applications. The brightness temperature can be computed from the digital number (DN) of TM6 data using the equation provided by the National Aeronautics and Space Administration (NASA). However, a proper algorithm for retrieving LST from the only one thermal band of the sensor still remains unavailable due to many difficulties in the atmospheric correction. Based on thermal radiance transfer equation, an attempt has been made in the paper to develop a mono-window algorithm for retrieving LST from Landsat TM6 data. Three parameters are required for the algorithm: emissivity, transmittance and effective mean atmospheric temperature. Method about determination of atmospheric transmittance is given in the paper through the simulation of atmospheric conditions with LOWTRAN 7 program. A practicable approach of estimating effective mean atmospheric temperature from local meteorological observation is also proposed in the paper when the in situ atmospheric profile data is unavailable at the satellite pass, which is generally the case in the real world especially for the images in the past. Sensitivity analysis of the algorithm indicates that the possible error of ground emissivity, which is difficult to estimate, has relatively insignificant impact on the probable LST estimation error i T, which is sensible to the possible error of transmittance i 6 and mean atmospheric temperature i T a . Validation of the simulated data for various situations of seven typical atmospheres indicates that the algorithm is able to provide an accurate LST retrieval from TM6 data. The LST difference between the retrieved and the simulated ones is less than 0.4°C for most situations. Application of the algorithm to the sand dunes across the Israel-Egypt border results in a reasonable LST estimation of the region. Based on this LST estimation, spatial variation of the interesting thermal phenomenon has been analysed for comparison of LST difference across the border. The result shows that the Israeli side does have significantly higher surface temperature in spite of its denser vegetation cover than the Egyptian side where bare sand is prevalent.     

波谱响应函数和波宽对地表温度反演的影响

分析了波谱响应函数和波宽对红外数据反演地表温度的影响。结果表明,波谱响应函数对温度反演的影响与传感器波段设置相关,误差随波宽变大而增加。Landsat、CBERS和环境卫星等宽波段红外数据的通道辐亮度和温度的关系,用有效波长比中心波长代替积分效果更好,但都存在一定误差,针对这些传感器数据的单通道算法须考虑波谱响应的差别。窄波段的MDIS地表温度反演,有效波长比中心波长代替积分效果更好,能满足精度要求。Terra和Aqua的MODIS波谱响应函数差别很小,可相互结合提高地表温度的反演精度。     

基于Landsat TM图像的北京城市地表温度遥感反演研究

利用北京地区Landsat TM热红外波段数据，采用单通道算法反演得到北京地区地面温度分布图。从反演结果可以看出，北京城区地面温度比郊区地表温度高，郊区地表温度较低，密云水库、官厅水库等水体的温度最低，总体上北京城市热岛效应显著。地表比辐射率是通过Van经验公式反演得到，通过对比分析，表明该方法对自然地表的比辐射率反演效果较好。     

Disaggregation of land surface temperature over a heterogeneous urban and surrounding suburban area: a case study in Shanghai,China

Many application fields need land surface temperature (LST) with simultaneous high spatial and temporal resolution, which can be achieved through the disaggregation technique. Most published methods built an assumed scale-independent relationship between LST and predictor variables derived from coarse spatial resolution data. However, LST disaggregation in the heterogeneous areas, especially urban areas, is very difficult to achieve and there are few studies on it. In this article, we propose an adjusted stratified stepwise regression method for temperature disaggregation in urban areas. Landsat Enhanced Thematic Plus (ETM+) data from Shanghai, China, were used to construct remote-sensing indices that are related to LST variance and retrieve LST at 60 and 480 m spatial resolution, respectively. Different stepwise regression models at 480 m resolution were built for two stratified regions according to normalized difference vegetation index (NDVI) distribution, and then each independent variable at 60 m resolution was adjusted to calculate disaggregated LST by considering its relationship with the 480 m resolution image. By using LST retrieved directly from ETM+ band 6 at 60 m spatial resolution as the reference, the proposed method comprising resampling disaggregation, the thermal data sharpening model (TsHARP)/disaggregation procedure for radiometric surface temperature (DisTrad) technique, and the LST-principal component analysis (PCA) regression algorithm were verified and compared. The results show that the temperature distribution estimated using the proposed method is most consistent with that of the reference LST in this heterogeneous study area, and that the precision improves significantly, especially for the low vegetation fraction region.     

基于Landsat8的地表温度反演算法研究——以滇池流域为例

2013年2月11日Landsat 8在加州范德堡空军基地发射升空,其携带的热红外传感器为反演地表温度提供了一种新的数据,但目前尚没有针对Landsat 8热红外波段反演地表温度的算法。针对Landsat 8第10波段特征,对现有反演地表温度的单窗算法进行了参数修正,得到了用Landsat 8第10波段反演地表温度的单窗算法系数。为了评价修正后算法的精度,用MODTRAN模拟地表温度为20、30和40℃时大气水汽含量分别为1.0、1.5、2.0和2.5g·cm-2传感器高度处的热辐射值,再将模拟数据用修正后的单窗算法反演地表温度,结果表明:地表温度越低、大气水汽含量越低,误差越小;模拟结果的平均误差为0.74℃。说明基于Landsat 8第10波段用修正后的单窗算法反演地表温度是可行的,该方法可为地表温度反演提供一种途径。最后以滇池流域为例,基于2013年4月20日的Landsat 8热红外数据反演了滇池流域的地表温度,并分析了滇池流域地表温度的分布特征。     

Comparison of impervious surface area and normalized difference vegetation index as indicators of surface urban heat island effects in Landsat imagery

This paper compares the normalized difference vegetation index (NDVI) and percent impervious surface as indicators of surface urban heat island effects in Landsat imagery by investigating the relationships between the land surface temperature (LST), percent impervious surface area (%ISA), and the NDVI. Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) data were used to estimate the LST from four different seasons for the Twin Cities, Minnesota, metropolitan area. A map of percent impervious surface with a standard error of 7.95% was generated using a normalized spectral mixture analysis of July 2002 Landsat TM imagery. Our analysis indicates there is a strong linear relationship between LST and percent impervious surface for all seasons, whereas the relationship between LST and NDVI is much less strong and varies by season. This result suggests percent impervious surface provides a complementary metric to the traditionally applied NDVI for analyzing LST quantitatively over the seasons for surface urban heat island studies using thermal infrared remote sensing in an urbanized environment.     

Improvements in land surface temperature retrieval based on atmospheric water vapour content and atmospheric temperature

Three methods are currently used to retrieve land surface temperatures (LSTs) from thermal infrared data supplied by the Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) sensors: the radiative transfer equation, mono-window, and generalized single-channel algorithms. Most retrieval results obtained using these three methods have an average error of more than 1 K. But if the regional mean atmospheric water vapour content and temperature are supplied by in situ radiosounding observations, the mono-window algorithm is able to provide better results, with a mean error of 0.5 K. However, there are no in situ radiosounding data for most regions. This article provides an improved method to retrieve LST from Landsat TM and ETM+ data using atmospheric water vapour content and atmospheric temperature, which can be obtained from remote-sensing data. The atmospheric water vapour content at the pixel scale was first calculated from Moderate Resolution Imaging Spectroradiometer (MODIS) data. The emissivities of various land covers and uses were then defined by Landsat TM or ETM+ data. In addition, the temperature–vegetation index method was applied to map area-wide instantaneous near-surface air temperatures. The parameters of mean atmospheric water vapour content and temperature and land surface emissivity were finally inputted to the mono-window algorithm to improve the LST retrieval precision. Our results indicate that this improved mono-window algorithm gave a significantly better retrieval of the estimated LST than that using the standard mono-window algorithm, not only in dry and elevated mountain regions but also in humid regions, as shown by the bias, standard deviation (σ), and root mean square deviation (RMSD). In Madoi County, the improved mono-window algorithm validated against the LST values measured in situ produced a bias and RMSD of –0.63 K and 0.91 K, respectively, compared with the mono-window algorithm’s bias and RMSD of –1.08 K and 1.27 K. Validated against the radiance-based method, the improved algorithm shows bias and RMSD values of –1.08 K and 1.27 K, respectively, compared with the initial algorithm’s bias and RMSD –1.65 K and 1.75 K. Additionally, the improved mono-window algorithm also appeared to be more accurate than the mono-window algorithm, with lower error values when validated against in situ measurement and the radiance-based method in the validation area in Zhangye City, Gansu Province, China. Remarkable LST accuracy improvements are shown by the improved mono-window algorithm, with better agreement not only with the in situ measurements but also with the simulated LSTs in the two validation areas, indicating the soundness and suitability of this method.     

作物蒸散发模型研究中的地表温度反演

针对美国加州Merced县2002年8月9日的ETM+影像,利用单窗温度反演算法反演了遥感蒸散发模型S SEBI中的地表温度参数。选用了ETM+热红外的高增益61波段,对热红外波段反射率较低的植被覆盖研究区进行了地表温度反演,并反演地表温度所需要的几个参数：亮度温度、地表比辐射率、大气透射率。最后得出了研究区域地表温度分布结果：水体地表温度低于植被作物,建筑或道路的地表温度最高。不同地物间地温是不同的,作为蒸散发反演的重要参数,这将影响不同地物蒸散发估算。因此精确反演地表温度,将为今后蒸散发的研究打好基础。     

基于Landsat-5 TM的洪河湿地地表温度估算方法对比研究

利用洪河湿地2008年5月15日过境的Landsat/TM图像和实测地面数据以及MODIS 地表发射率数据,分别运用大气辐射传输模型、覃志豪的单窗算法和Jimenez\|Munoz & Sobrino 的单波段算法估算洪河湿地的地表温度,并且对比了大气校正前后的NDVI、LSE以及各种算法估算地表温度的差异。分析估算结果表明,覃志豪的单窗算法与实测地面数据估算结果非常一致。指出在没有实时探空数据的情况下,应用只有一个热红外通道的Landsat/TM数据源,采用覃志豪的单窗算法估算的精度是可以接受的。     

Thermal anomalies associated with faults: a case study of the Jinhua–Quzhou basin of Zhejiang Province,China

Faults provide the path for geothermal natural convection and partially influence the ground surface thermal environment. The land surface temperatures (LSTs) near a fault are higher than in other areas and can indicate the strike trend of an underground fault. However, these anomalies of higher LSTs are not located accurately in the fault centre but near it with some offset, and these LST data may include other thermal information that needs to be eliminated prior to analysis. In this study, LSTs were retrieved from Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) thermal infrared (TIR) images and enhanced with land cover classification and elimination. The spatial patterns of the enhanced images were compared with geophysical prospecting tectonic profiles and with regional geological tectonic maps, revealing the spatial correspondence between the thermal anomalies and the faults. The results indicated that the thermal anomalies are located near the faults and are consistent with the faults' dip planes.     

Evaluation of dynamic linkages between evapotranspiration and land-use/land-cover changes with Landsat TM and ETM+ data

Complexity embedded in coastal management leads to numerous questions as to how inherent spatial and temporal linkages among evapotranspiration (ET), depth to groundwater and land-use/land-cover change (LUCC) could affect the dynamics among these seemingly unrelated events. This article aims to address such unique dynamics in the nexus of physical geography and ecohydrology. To understand such dynamic linkages, a case study was carried out in a fast growing coastal region – the southern Laizhou Bay in Shandong Province, China – by identifying the coastal LUCC at the decadal scale in association with the variations of ET with the aid of Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM+) data. In such a coastal landscape evolutionary assessment, findings show that the major patterns of land use and land cover (LULC) in the study area are farmland, saline-alkali land, developed land, salt land and beach land. Over a 20-year time frame, declining groundwater trends were observed, while ET increased gradually with changing LULC. By using the surface energy balance algorithm for land (SEBAL) with Landsat TM/ETM+ images and additional environmental data, the concomitant response of ET variations due to LUCC becomes lucid among three significantly correlated pairs including fractional vegetation cover (FVC), land surface temperature (LST) and soil heat flux. The dynamic linkages between ET and LULC were finally confirmed with such a pair-wise analysis.     

Land surface temperature retrieval from LANDSAT TM 5

In this paper, three methods to retrieve the land surface temperature (LST) from thermal infrared data supplied by band 6 of the Thematic Mapper (TM) sensor onboard the Landsat 5 satellite are compared. The first of them lies on the estimation of the land surface temperature from the radiative transfer equation using in situ radiosounding data. The others two are the mono-window algorithm developed by Qin et al. [International Journal of Remote Sensing 22 (2001) 3719] and the single-channel algorithm developed by Jiménez-Muñoz and Sobrino [Journal of Geophysical Research 108 (2003)]. The land surface emissivity (LSE) values needed in order to apply these methods have been estimated from a methodology that uses the visible and near infrared bands. Finally, we present a comparison between the LST measured in situ and the retrieved by the algorithms over an agricultural region of Spain (La Plana de Requena-Utiel). The results show a root mean square deviation (rmsd) of 0.009 for emissivity and lower than 1 K for land surface temperature when the Jiménez-Muñoz algorithm is used.     

A method to detect and correct single-band missing pixels in Landsat TM and ETM+ data

Human land use is the main driver of terrestrial ecosystems change, and remote sensing is an important tool to monitor these changes. Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM+) images have been the most important data source to map land cover change, but image artifacts often hinder or even prohibit digital change detection. This paper addresses a group of image distortions that display erroneous values in a single band while leaving the other bands of a spectrum undisturbed. Such artifacts may be due to different phenomena, for instance transmission and ground-processing problems or single event upsets. Automated artifact detection for those phenomena is often difficult, because erroneous band values often lie well within the range of naturally occurring radiance values. We developed IDL-based software that uses edge operators to detect and label affected pixels. Using a least-squares spectral-matching algorithm, the distorted spectrum is compared with undisturbed spectra in the local neighborhood and the undisturbed spectrum of best fit is determined. The erroneous band value is then replaced with the corresponding undisturbed value. This method was tested on seven Landsat TM images and on artificial data. Our results show that the distorted areas are precisely detected and that the correction procedure leads to meaningful spectra. This approach may be useful to minimize the effect of single-band distortions and allows for subsequent image analysis without the need to mask out distorted areas. The software tool includes a user interface and is available online.     

基于劈窗算法的Landsat 8影像地表温度反演

陆地表面温度（LST）是表征地表能量交换和地面特征的重要指标,目前遥感技术逐渐成为区域和全球尺度上LST反演的一种便捷工具,而采样不同算法及不同影像的热红外遥感LST反演研究层出不穷,其中基于Landsat数据的反演成果尤为突出。文章利用劈窗算法对Landsat 8遥感影像进行地表温度反演,对比探讨了根据经验值与借助MODIS热红外数据两种不同方式的LST反演结果,并进行北京市热红外波段辐射亮度温度比较,针对地表温度分级进行统计,分析了当地地表温度分布趋势。结果表明：劈窗算法下Landsat 8数据的反演温度更接近实际温度,精度较高且优于MODIS产品;北京市地表温度空间分布格局受地物结构与反射率所制约,高温区主要集中分布于中东部,中低温区分布与林地及水体分布结构较为吻合。