遥测地表温度与比辐射率的迭代反演方法──理论推导与数值模拟

首先在地表比辐射率为已知的条件下,提出一个非线性迭代温度反演模型,我们对不同的地表和大气条件进行了模拟计算,结果表明当大气温度廓线误差－２Ｋ－２Ｋ,水汽廓线误差±２０％时的温度均方根误差为０．４８Ｋ。当大气模式误差一个模式时反演的温度均方根误差为０．８５Ｋ。在此基础上,引人相邻像元的概念,相邻像元的大气状况可以认为是相同的,应用两个时相的遥感影像数据,假定在两个相近时相之间地表比辐射率值不变,建立地表比辐射率与温度的反演模型。我们对不同的地表和大气条件进行了模拟计算,结果表明当大气温度廓线误差－２Ｋ－２Ｋ,水汽廓线误差±２０％时地表温度均方根误差小于１．５Ｋ,地表比辐射率均方根误差小于０．０２,地表辐射均方根误差为１％;当大气温度廓线误差－２Ｋ－２Ｋ,水汽廓线误差±１０％时,地表温度均方根误差小于１．０Ｋ,地表辐射均方根误差小于０．６％。     

陆地卫星TM6波段范围内地表比辐射率的估计

地表比辐射率是用热红外波段遥感数据反演地表温度的关键参数。目前,应用陆地卫星TM6波段数据反演地表温度共有3种算法,即大气校正法、单窗算法和单通道算法。这3种算法都需要TM6波段范围内的地表比辐射率作为地表温度反演参数。本文首先简介这3种反演算法;然后着重探讨TM6波段地表比辐射率估计方法;最后,利用这一方法对山东省陵县附近农田地区进行地表比辐射率估计和地表温度反演。结果表明,该方法能获得较合理的地表温度反演结果。     

Contribution of Landsat 8 data for the estimation of land surface temperature in Batna city,Eastern Algeria

In this study, we presented a mono-window (MW) algorithm for land surface temperature retrieval from Landsat 8 TIRS. MW needs spectral radiance and emissivity of thermal infrared bands as input for deriving LST. The spectral radiance was estimated using band 10, and the surface emissivity value was derived with the help of NDVI and vegetation proportion parameters for which OLI bands 5 and 4 were used. The results in comparison with MODIS (MOD11A1) products indicated that the proposed algorithm is capable of retrieving accurate LST values, with a correlation coefficient of 0.850. The industrial area, public facilities and military area show higher surface temperature (more than 37 °C) in comparison with adjoining areas, while the green spaces in urban areas (34 °C) and forests (29 °C) were the cooler part of the city. These successful results obtained in the study could be used as an efficient method for the environmental impact assessment.     

Effects of prediction accuracy of the proportion of vegetation cover on land surface emissivity and temperature using the NDVI threshold method

Predicting land surface energy budgets requires precise information of land surface emissivity (LSE) and land surface temperature (LST). LST is one of the essential climate variables as well as an important parameter in the physics of land surface processes at local and global scales, while LSE is an indicator of the material composition. Despite the fact that there are numerous publications on methods and algorithms for computing LST and LSE using remotely sensed data, accurate prediction of these variables is still a challenging task. Among the existing approaches for calculating LSE and LST, particular attention has been paid to the normalised difference vegetation index threshold method (NDVI^THM), especially for agriculture and forest ecosystems. To apply NDVI^THM, knowledge of the proportion of vegetation cover (P_V) is essential. The objective of this study is to investigate the effect of the prediction accuracy of the P_V on the estimation of LSE and LST when using NDVI^THM. In August 2015, a field campaign was carried out in mixed temperate forest of the Bavarian Forest National Park, in southeastern Germany, coinciding with a Landsat-8 overpass. The P_V was measured in the field for 37 plots. Four different vegetation indices, as well as artificial neural network approaches, were used to estimate P_V and to compute LSE and LST. The results showed that the prediction accuracy of P_V improved using an artificial neural network (R²_CV = 0.64, RMSE_CV = 0.05) over classic vegetation indices (R²_CV = 0.42, RMSE_CV = 0.06). The results of this study also revealed that variation in the accuracy of the estimated P_V affected calculation results of the LSE. In addition, our findings revealed that, though LST depends on LSE, other parameters should also be taken into account when predicting LST, as more accurate LSE results did not increase the prediction accuracy of LST.     

基于大气温度和水汽含量改进HJ-1B地表温度反演

为了提高地面气象站稀少地区地表温度遥感反演的精度,本文基于多源遥感数据的优势,首先利用MODIS影像获取研究区像元尺度上平均大气水汽含量;然后利用同时相的HJ-1B影像估算区域地表比辐射率,再采用温度-植被指数法获取近地表大气温度;最后将以上3个参数输入单窗体算法,改进其地表温度反演的精度。研究结果表明,改进单窗体算法反演地表温度与地面实测温度的偏差小于1 K,为地面气象站点稀少的植被覆盖区域提供了一种可行的精确遥感反演地表温度方法。     

Investigation and validation of algorithms for estimating land surface temperature from Sentinel-3 SLSTR data

Land surface temperature (LST) is an important indicator of global ecological environment and climate change. The Sea and Land Surface Temperature Radiometer (SLSTR) onboard the recently launched Sentinel-3 satellites provides high-quality observations for estimating global LST. The algorithm of the official SLSTR LST product is a split-window algorithm (SWA) that implicitly assumes and utilizes knowledge of land surface emissivity (LSE). The main objective of this study is to investigate alternative SLSTR LST retrieval algorithms with an explicit use of LSE. Seventeen widely accepted SWAs, which explicitly utilize LSE, were selected as candidate algorithms. First, the SWAs were trained using a comprehensive global simulation dataset. Then, using simulation data as well as in-situ LST, the SWAs were evaluated according to their sensitivity and accuracy: eleven algorithms showed good training accuracy and nine of them exhibited low sensitivity to uncertainties in LSE and column water vapor content. Evaluation based on two global simulation datasets and a regional simulation dataset showed that these nine SWAs had similar accuracy with negligible systematic errors and RMSEs lower than 1.0 K. Validation based on in-situ LST obtained for six sites further confirmed the similar accuracies of the SWAs, with the lowest RMSE ranges of 1.57–1.62 K and 0.49−0.61 K for Gobabeb and Lake Constance, respectively. While the best two SWAs usually yielded good accuracy, the official SLSTR LST generally had lower accuracy. The SWAs identified and described in this study may serve as alternative algorithms for retrieving LST products from SLSTR data.     

基于ASTER GED产品的地表发射率估算

地表发射率是地表温度反演的重要输入参数,为了解决现有地表发射率估算方法在裸露地表精度较差的问题,本文基于最新的ASTER全球地表发射率产品(ASTER GED)和基于植被覆盖度的方法(VCM),提出了一个改进的地表发射率估算方法。首先,利用ASTER GED产品求解裸土发射率,然后,利用ASTER波谱库中的植被发射率和植被覆盖度结合VCM方法计算地表发射率。利用张掖地区2012年11景ASTER TES算法反演的地表发射率产品和实测地表发射率数据进行了验证,同时利用一景Landsat 8 TIRS数据分析了对地表温度反演精度的影响。结果表明该方法估算的地表发射率整体精度较高,可以有效改进裸露地表的发射率估算精度,用于支持利用多种热红外传感器数据生产高精度的地表温度产品。     

Land surface temperature and its impact factors in Western Sichuan Plateau,China

The understanding influence of multiple factors variations on land surface temperature (LST) remains elusive. LST was retrieved by the atmospheric correction algorithms. Based on the correlation coefficients, stepwise regression analysis was developed to examine how multiple factors variability led to LST variations. The differences in LST between impact factors vary depending on time in a day. The elevation and land use types significantly affect the LST in sunny slope or shadow areas has a significantly quadratic curve correlation or a negative linear correlation with it, the influence of slope and aspect is not very significant. LST for forestland, grassland and bare land in the sunny slope and shadow area was the cubic polynomial related to its elevation. Normalized difference vegetation index (NDVI) and normalized difference moisture index (NDMI) effectively express LST in mountainous. LST and NDMI or NDVI have a significantly negative correlation, NDMI is more effective and more applicable for the expression of LST.     

AAC算法的大气校正复合改进算法

利用模拟数据,评价Autonomous Atmospheric Compensation(AAC)算法的抗噪性,认为AAC算法的抗噪性较弱。基于TASI实测数据,利用AAC算法开展反演计算时,计算结果呈现出多样性问题。结合In-scene Atmospheric Compensation(ISAC)算法中黑体像元的标定方法,提出了一种复合改进算法。首先,利用ISAC算法反演的大气透过率和路径辐射,重新计算AAC算法中大气透过率之比(Tr)和相邻两强弱吸收通道的路径辐射之差(Pd),再次,运用经验公式获得稳定的大气反演结果(大气透过率和路径辐射),有效解决了计算结果多样性的问题。利用复合改进算法,开展的温度与发射率分离实验,证明反演得到的发射率波谱更接近野外实测波谱。     

单窗算法的大气参数估计方法

根据地表热辐射传导方程,提出了一个简单可行且精度较高的专门用于从TM6数据中演算地表温度的方法——单窗算法.这一算法把大气和地表状态对地表热传导的影响直接包括在演算公式中.该方法需要两个大气参数进行地表温度的演算,即大气平均作用温度和大气透射率.本文论述这两个大气参数的估计方法:根据大气水分含量或地表附近空气湿度来估计大气透射率;通过分析标准大气剖面资料,尤其是大气水分和气温随高程的变化规律,根据地表温度推算大气平均作用温度.     

Landsat卫星热红外数据地表温度遥感反演研究进展

作为驱动地表与大气之间能量交换的关键物理量,地表温度在众多领域中都发挥着重要作用,包括气候变化、环境监测、蒸散发估算以及地热异常勘探等。Landsat热红外数据因其时间连续性和高空间分辨率等特点被广泛应用于地表温度反演中。本文详细地介绍了Landsat热红外传感器及其可用的数据与产品的现状,梳理了2001年—2020年20年间基于Landsat热红外数据的地表温度遥感反演与应用的相关文献发表及互引情况,系统地综述了基于Landsat热红外数据的地表温度反演算法,包括基于辐射传输方程的算法、单窗算法、普适性单通道算法、实用单通道算法和分裂窗算法等。在此基础上,进一步介绍了每种算法的参数化方案,包括地表比辐射率和大气参数的估算方法。最后针对Landsat热红外数据地表温度遥感反演提出了未来可能的发展趋势与研究方向。     

An approach of surface coal fire detection from ASTER and Landsat-8 thermal data: Jharia coal field,India

Radiant temperature images from thermal remote sensing sensors are used to delineate surface coal fires, by deriving a cut-off temperature to separate coal-fire from non-fire pixels. Temperature contrast of coal fire and background elements (rocks and vegetation etc.) controls this cut-off temperature. This contrast varies across the coal field, as it is influenced by variability of associated rock types, proportion of vegetation cover and intensity of coal fires etc. We have delineated coal fires from background, based on separation in data clusters in maximum v/s mean radiant temperature (13th band of ASTER and 10th band of Landsat-8) scatter-plot, derived using randomly distributed homogeneous pixel-blocks (9 × 9 pixels for ASTER and 27 × 27 pixels for Landsat-8), covering the entire coal bearing geological formation. It is seen that, for both the datasets, overall temperature variability of background and fires can be addressed using this regional cut-off. However, the summer time ASTER data could not delineate fire pixels for one specific mine (Bhulanbararee) as opposed to the winter time Landsat-8 data. The contrast of radiant temperature of fire and background terrain elements, specific to this mine, is different from the regional contrast of fire and background, during summer. This is due to the higher solar heating of background rocky outcrops, thus, reducing their temperature contrast with fire. The specific cut-off temperature determined for this mine, to extract this fire, differs from the regional cut-off. This is derived by reducing the pixel-block size of the temperature data. It is seen that, summer-time ASTER image is useful for fire detection but required additional processing to determine a local threshold, along with the regional threshold to capture all the fires. However, the winter Landsat-8 data was better for fire detection with a regional threshold.     

针对Terra/MODIS数据的改进分裂窗地表温度反演算法

针对Terra/MODIS数据提出改进的分裂窗地表温度反演算法。充分考虑了传感器观测角度(VZA)的影响,并对地表和有效大气辐射按照不同的亮度温度区间分别进行Planck函数简化。利用TIGR3大气廓线库中的875条晴空大气廓线,ASTER波谱库中的106条地物发射率波谱,结合MODTRAN4大气辐射传输模型模拟得到分裂窗算法系数。利用MODTRAN4模拟数据对算法精度进行验证,结果表明本文的改进算法和原算法的均方根误差RMSE分别为0.34K和0.65K。敏感性分析表明,在中等湿润的大气条件下,算法对大气水汽含量并不敏感。该算法降低了传感器观测角度带来的地表温度反演误差。利用2009年6月美国SURFRAD辐射观测网6个站点的实测数据对改进算法、原算法以及MOD11_L2地表温度产品进行了对比验证,RMSE分别是0.93K、1.49K和1.0K,表明本文算法可以提高反演精度。     

热红外遥感地表温度与发射率地面验证进展

地表温度与发射率是地表—大气系统长波辐射和潜热通量交换的直接驱动力,是描述区域和全球尺度上地表能量平衡与水平衡的重要参数,其时空变化信息在气象预测、气候变化、水循环、地质勘探、农林监测和城市热环境等诸多领域具有广泛的应用。热红外遥感作为当前获取区域或全球尺度上地表温度和发射率的最有效手段之一,相较于传统的地面点位测量方法,具有空间覆盖范围大和重复观测等优势。对热红外遥感定量反演的地表温度与发射率产品进行地表真实性验证,有利于发现遥感数据自身或其反演算法的缺陷,确定产品的精度与不确定度,便于遥感产品的应用与推广。本文首先回顾了地表温度和发射率的定义,阐述了热红外遥感可反演、地面可测量的地表温度和发射率的科学内涵,并对利用热红外遥感数据反演地表温度和发射率的理论和方法作了概述;对地表温度和发射率地面验证的框架体系、验证指标进行总结,建立了基于精度、精确度、不确定度、完整性和稳定性的验证评价指标体系;总结了地表温度和地表发射率的地面验证方法、地面测量方法、辅助数据的获取方法、地表温度地面测量的采样方法,以及在验证异质非同温地表时从点到像元尺度的地表温度尺度转换方法等,分析了地面验证过程的主要误差来源;归纳了目前地表温度和地表发射率主要验证站点、观测网络及其空间分布特征;最后,本文讨论了地表温度与发射率地面验证存在的若干问题,并对地表温度与发射率验证工作的发展前景和趋势进行了相关展望。     

Land surface temperature and emissivity estimation for Urban Heat Island assessment using medium- and low-resolution space-borne sensors: A review

The Urban Heat Island (UHI) phenomenon, a typical characteristic on urban landscapes, has been recognised as a key driver to the transformation of local climate. Reliable retrieval of urban and intra-urban thermal characteristics using satellite thermal data depends on accurate removal of the effects of atmospheric attenuations, angular and land surface emissivity. Several techniques have been proposed to retrieve land surface temperature (LST) from coarse resolution sensors. Medium spatial resolution sensors like the Advanced Space-borne Thermal Emission and Reflection Radiometer and the Landsat series offer a viable option for assessing LST within urban landscapes. This paper reviews the theoretical background of LST estimates from the thermal infrared part of the electromagnetic spectrum, LST retrieval algorithms applicable to each of the commonly used medium-resolution sensors and required variables for each algorithm. The paper also highlights LST validation techniques and concludes by stipulating the requirements for LST temporal and spatial configuration.     

基于ASTER数据反演我国南方山地陆表温度

 以贵州省黎平县山地植被覆盖区为例，基于ASTER遥感数据进行15 m分辨率的归一化植被指数制图和地表发射率制图，在利用 MODTRAN 4大气辐射传输模型进行大气订正的基础上，基于普朗克辐射方程的推导反演陆表温度，取得了较为理想的结果。     

Angular effect of MODIS emissivity products and its application to the split-window algorithm

The angular effects of emissivity are ignored in current land surface temperature (LST) products. As a result, the directionality of these LST products limits their further application in many fields. Accurate correction of the angular problem of LST products requires explicit understanding of the angular effects of emissivity at the pixel scale. Currently, nearly ten years of global emissivity products of MODIS are available. However, the pixel-scale directionality of emissivity has never been analyzed. By performing a statistical analysis of 5-year MODIS emissivity products over most of East Asia, we generated the empirical relationships between the directional emissivity, land cover, and seasonal variations. Two look-up tables (LUTs) of directional emissivity were created for typical land cover types and applied to the generalized split-window algorithm to modify the MODIS LST. The results showed that the angular effect of emissivity could introduce a significant bias of −1-3 K to the 1 km resolution LST. Finally, the spatial scale effects of emissivity were analyzed, and it was found that the temperature differences caused by scale effects fell within +/−0.5 K for most pixels if 5 km emissivity was used in 1 km LST retrieval. Therefore, wide use of the LUTs can be expected.     

Satellite-based analysis of the role of land use/land cover and vegetation density on surface temperature regime of Delhi,India

The knowledge of the surface temperature is important to a range of issues and themes in earth sciences central to urban climatology, global environmental change and human-environment interactions. The study analyses land surface temperature (LST) estimation using temporal ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) datasets (day time and night time) over National Capital Territory Delhi using the surface emissivity information at pixel level. The spatial variations of LST over different land use/land cover (LU/LC) at day time and night time were analysed and relationship between the spatial distribution of LU/LC and vegetation density with LST was developed. Minimum noise fraction (MNF) was used for LU/LC classification which gave better accuracy than classification with original bands. The satellite derived emissivity values were found to be in good agreement with literature and field measured values. It was observed that fallow land, waste land/bare soil, commercial/industrial and high dense built-up area have high surface temperature values during day time, compared to those over water bodies, agricultural cropland, and dense vegetation. During night time high surface temperature values are found over high dense built-up, water bodies, commercial/industrial and low dense built-up than over fallow land, dense vegetation and agricultural cropland. It was found that there is a strong negative correlation between surface temperature and NDVI over dense vegetation, sparse vegetation and low dense built-up area while with fraction vegetation cover, it indicates a moderate negative correlation. The results suggest that the methodology is feasible to estimate NDVI, surface emissivity and surface temperature with reasonable accuracy over heterogeneous urban area. The analysis also indicates that the relationship between the spatial distribution of LU/LC and vegetation density is closely related to the development of urban heat islands (UHI).     

比辐射率封闭法测定技术及误差分析

劈窗算法是目前由热红外遥感图像数据获取陆面温度最主要的方法。由于进行地面像元尺度实时陆面温度同步测量的困难,尚无法直接对现有各劈窗算法进行评判。该文借助于辐射传输模型ＬＯＷＴＲＡＮ７及其提供的６ 种标准大气模式,进行模拟计算,分析了６ 种主要劈窗算法对大气廓线误差和比辐射率的敏感度,作为对劈窗算法适用性的一间接判据。     

基于ASTER数据的城市热环境遥感监测研究

以ASTER数据为数据源,采用同一颗卫星上的MODIS数据得到大气透过率;利用可见光和近红外波段对下垫面类型进行分类和利用JPL(Jet Propulsion Lab)提供的光谱库计算地表比辐射率,进而采用劈窗算法进行地表温度(Land Surface Temperature,LST)的反演。在此基础上,利用反演的LST、分类结果和归一化差值植被指数(NDVI),对沧州地区的城市热环境进行了定量分析,研究结果可为进一步深入探讨城市热岛的发生发展规律以及城市热环境的模拟调控、优化布局提供一定的科学依据。