中国陆地区域陆表温度业务化遥感反演算法及产品运行系统

地表温度反演的裂窗算法已成功应用于NOAA系列卫星热红外遥感数据。目前，裂窗算法中应用较为广泛的一种是Becker等人于1990年提出的局地裂窗算法，主要是通过辐射传输模型模拟不同地表条件和大气状况下，地表温度和发射率对红外辐射亮温的影响，从而发展出一个利用AVHRR4，5通道亮温数据反演地表温度的线性模型。在晴空无云和地表比辐射率能精确估算的情况下，Becker算法反演地表温度的精度在1K以内。Becker算法用Lowtran程序模拟计算地表辐射量，且模型中参数主要针对NOAA-9传感器特性得到。本文在Becker算法的基础上，针对NOAA-16／17传感器热红外通道光谱响应函数特性，利用最新的、计算光谱分辨率更高的MODTRAN程序模拟不同大气状况下，不同地表温度和发射率对NOAAAVHRR4，5通道辐射亮温响应特性的影响，改进Becker算法中模型参数，使之能适用于NOAA-16／17热红外数据。同时，本文利用植被指数NDVI，在中国陆地区域lkm分辨率最新地表分类数据的基础上，得到模型中需要的地表比辐射率参数，将改进的模型应用于1km分辨率NOAA17数据，得到了旬合成中国陆地区域范围地表温度，通过地面气象台站实测数据对比验证．取得了较好的结果。     

Retrieval of land surface temperature (LST) from landsat TM6 and TIRS data by single channel radiative transfer algorithm using satellite and ground-based inputs

The present study proposes land surface temperature (LST) retrieval from satellite-based thermal IR data by single channel radiative transfer algorithm using atmospheric correction parameters derived from satellite-based and in-situ data and land surface emissivity (LSE) derived by a hybrid LSE model. For example, atmospheric transmittance (τ) was derived from Terra MODIS spectral radiance in atmospheric window and absorption bands, whereas the atmospheric path radiance and sky radiance were estimated using satellite- and ground-based in-situ solar radiation, geographic location and observation conditions. The hybrid LSE model which is coupled with ground-based emissivity measurements is more versatile than the previous LSE models and yields improved emissivity values by knowledge-based approach. It uses NDVI-based and NDVI Threshold method (NDVI^THM) based algorithms and field-measured emissivity values. The model is applicable for dense vegetation cover, mixed vegetation cover, bare earth including coal mining related land surface classes. The study was conducted in a coalfield of India badly affected by coal fire for decades. In a coal fire affected coalfield, LST would provide precise temperature difference between thermally anomalous coal fire pixels and background pixels to facilitate coal fire detection and monitoring. The derived LST products of the present study were compared with radiant temperature images across some of the prominent coal fire locations in the study area by graphical means and by some standard mathematical dispersion coefficients such as coefficient of variation, coefficient of quartile deviation, coefficient of quartile deviation for 3rd quartile vs. maximum temperature, coefficient of mean deviation (about median) indicating significant increase in the temperature difference among the pixels. The average temperature slope between adjacent pixels, which increases the potential of coal fire pixel detection from background pixels, is significantly larger in the derived LST products than the corresponding radiant temperature images.     

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

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

Land surface temperature retrieval and its validation using NOAA AVHRR thermal data

The retrieval of land (soil-vegetation complex) surface temperature (LST) was carried out over semi-arid mixed agriculture landscape of Gujarat using thermal bands (channel 4 and 5) and ground emissivity from atmospherically corrected NDVI of NOAA AVHRR LAC images. The atmospheric correction of Visible and NIR band reflectance was done using SMAC model. The LST computed from split-window method and subsequently corrected with fractional vegetation cover were then compared with near synchronous ground observations of soil and air temperatures made during 13–17 January and April, 1997 at five Land Surface Processes Experiment (LASPEX) sites of Anand, Sanand, Derol, Arnej and Khandha covering 100 km x 100 km. The fractional vegetation cover corrected LST at noon hrs. varied from 301.6 – 311.9K in January and from 315.8 – 325.6K in April. The LST_corr were found to lie in the mid way between AT and ST during January. But in April, LST were found to be more close to ST which may be due to relatively poor vegetation growth as indicated by lower NDVI values in April indicating more contribution to LST from exposed soil surface.     

城市地表温度空间降尺度研究——以北京市为例

地表温度LST (Land Surface Temperature)是城市热环境研究的重要参数之一,城市下垫面极为复杂,LST空间差异性较高.高空间分辨率LST对精细化城市热环境监测和缓解具有重要意义.目前大部分城市遥感LST降尺度研究仍以二维角度为主,缺乏建筑三维结构的考虑.本研究同时考虑地表二维和三维指标,构建基于...     

Retrieval of land surface temperature from the Kalpana-1 VHRR data using a single-channel algorithm and its validation over western India

Indian geostationary satellite Kalpana-1 (K1) offers a potential to capture the diurnal cycle of land surface temperature (LST) through thermal infrared channel (10.5–12.5 μm) observations of the Very High Resolution Radiometer (VHRR) sensor. A study was carried out to retrieve LST by adapting a generalized single-channel (SC) algorithm (Jiménez-Muñoz and Sobrino, 2003) for the VHRR sensor over India. The basis of SC algorithm depends on the concept of Atmospheric Functions (AFs) that are dependent on transmissivity, upwelling and downwelling radiances of the atmosphere. In the present study AFs were computed for the VHRR sensor through the MODTRAN simulations based upon varying atmospheric and surface inputs. The AFs were fitted with the atmospheric columnar water vapour content and a set of coefficients was derived for LST retrieval. The K1-LST derived with the SC algorithm was validated with (a) in situ measurements at two sites located in western parts of India and (b) the MODIS LST products. Comparison of K1-LST with the in situ measurements demonstrated that SC algorithm was successful in capturing the prominent diurnal variations of 283–332 K in the LST at desert and agriculture experimental sites with a rmse of 1.6 K and 2.7 K, respectively. Inter comparison of K1-LST and MODIS LST showed a reasonable agreement between these two retrievals up to LST of 300 K, however a cold bias up to 7.9 K was observed in MODIS LST for higher LST values (310–330 K) over the hot desert region.     

用被动微波AMSR数据反演地表温度及发射率的方法研究

 针对对地观测卫星多传感器的特点，提出了借助MODIS地表温度产品从被动微波数据中反演地表温度的方法。即利用MODIS地表温度产品和AMSR不同通道之间的亮度温度，建立地表温度的反演方程。该方法克服了以往需要测量同步数据的困难，为不同传感器之间的参数反演相互校正和综合利用多传感器的数据提供实际应用和理论依据。文中以MODIS地表温度产品作为评价标准，对方法进行检验，其平均误差为2~3℃。另外，微波的发射率是土壤水分反演的关键参数，在对微波地表温度反演的基础上，进一步对发射率进行了研究。     

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

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

A case study for intercomparison of land surface temperature retrieved from GOES and MODIS

In recent years, algorithms have been developed to derive land surface temperature (LST) from geostationary and polar satellite systems. However, few works have addressed the intercomparison between Geostationary Operational Environmental Satellites (GOES) and the available suite of polar sensors. In this study, differences in LSTs between GOES and MODerate resolution Imaging Spectroradiometer (MODIS) have been compared and also evaluated against ground observations. Due to the lack of split-window (SW) channels in the GOES M (12)-Q era, a dual-window algorithm using a mid-infrared 3.9 µm channel is compared with traditional SW algorithm. It is found that the differences in LST between different platforms are bigger during daytime than those during nighttime. During daytime, LSTs from GOES with the dual-window algorithm are warmer than MODIS LSTs, while LSTs from the SW algorithm are close to MODIS LSTs. The difference during daytime is found to be related to anisotropy in satellite viewing geometry, and land surface properties, such as vegetation cover and especially surface emissivity at middle infrared (MIR) channel. When evaluated against ground observations, the standard deviation (precision) error (2.35 K) from the dual window algorithm is worse than that (1.83 K) from the SW algorithm, indicating the lack of split-window channel in the GOES M(12)-Q era may degrade the performance of LST retrievals.     

青藏高原那曲地区MODIS 地表温度估算

地表温度是区域和全球尺度陆面过程研究中的一个关键参数,利用遥感卫星资料反演得到的地表温度数据在气象、水文和生态领域研究中有重要作用.本文基于改进后的针对MODIS 数据的分裂窗口算法,对MODIS L1B 卫星数据进行实用而简便的云检测处理,并根据青藏高原陆地、水体和冰雪等常见下垫面状况的遥感影像分类结果,反演得到了2007-01-03 、04-18 、06-12 和10-02 四日的无云下垫面地表温度.最后,将Sobrino 结果在青藏高原那曲地区与MODIS 日地表温度产品及CAMP/Tibet 观测站地表温度数据进行了对比验证分析.结果表明,该方法得到的地表温度结果与MODIS 数据产品具有较好的一致性,并且地表温度结果与地面观测数据(去除可疑点后)的平均误差仅为1.435 K .     

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

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

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

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

Atacama Field Campaign: laboratory and in-situ measurements for remote sensing applications

This work presents the preliminary results of the first field calibration campaign performed in the Atacama Desert, Chile, between the 18 and 22 August 2014, called the Atacama Field Campaign (ATAFIC 2014). In situ measurements were performed in order to spectrally characterize the surface reflectance spectra between 0.3 and 2.5?µm, radiometric temperature (8.0–14.0?µm) and atmospheric measurements. A soil sample was collected and analyzed using Fourier Transform Infrared Spectroscopy and X-Ray Diffraction techniques to characterize the surface reflectance spectra and mineralogical composition, respectively. ASTER land surface emissivity in addition to GOES, MODIS and Landsat-8 land surface temperature (LST) were also used. Results showed that the spectral features of the Atacama soil and the characteristics of this geographical zone, which is featured as the most hyper-arid and cloudless place in the world, make this area a potential target for surface reflectance characterization. Day and night LST comparison between field and remote sensing data are lower than 2?K and the Root Mean Square Error for land surface emissivity is close to 2%. This work opens the possibilities to consider the Atacama Desert as a reference target for calibration and validation activities for earth observation missions’ purposes.     

针对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,表明本文算法可以提高反演精度。     

MISR和MODIS二向性反射数据产品的对比分析

利用4种地表类型共16个地面站点的MISR 9个角度反射率数据产品和MODIS反射率数据产品(MOD09A1), 对比两种传感器获取的地表方向反射观测数据和用BRDF模型参数产品计算的方向反射数据, 分析了这两种BRDF模型参数的外推能力。研究结果表明: (1) MISR和MODIS BRDF模型参数数据产品都具有一定的对观测数据以外方向反射的外推能力, 相比用MISR BRDF模型参数数据外推到相应MODIS观测角度方向反射的结果一致性较强。(2) 采用两种BRDF模型参数推算的方向反射率与观测数据的接近程度有随观测天顶角的增大而减弱的趋势, 在观测天顶角较小时一致性较好。(3) 所用数据处理结果显示, MODIS BRDF模型参数数据产品在近垂直主平面方向反射率推算结果相比近主平面方向的推算结果更接近观测数据。     

Developing a temporally land cover-based look-up table (TL-LUT) method for estimating land surface temperature based on AMSR-E data over the Chinese landmass

The land surface temperature (LST) is an important parameter when studying the interface between the atmosphere and the Earth's surface. Compared to satellite thermal infrared (TIR) remote sensing, passive microwave (PMW) remote sensing is better able to overcome atmospheric influences and to estimate the LST, especially in cloudy regions. However, methods for estimating PMW LSTs at the country and continental scales are still rare. The necessity of training such methods from a temporally dynamic perspective also needs further investigations. Here, a temporally land cover based look-up table (TL-LUT) method is proposed to estimate the LSTs from AMSR-E data over the Chinese landmass. In this method, the synergies between observations from MODIS (Moderate Resolution Imaging Spectroradiometer) and AMSR-E (Advanced Microwave Scanning Radiometer for EOS), which are onboard the same Aqua satellite, are explored. Validation with the synchronous MODIS LSTs demonstrates that the TL-LUT method has better performances in retrieving LSTs with AMSR-E data than the method that uses a single brightness temperature in 36.5 GHz vertical polarization channel. The accuracy of the TL-LUT method is better than 2.7 K for forest and 3.2 K for cropland. Its accuracy varies according to land cover type, time of day, and season. When compared with the in-situ measured LSTs at four sites without urban warming in the Tibet Plateau, the standard errors of estimation between the estimated AMSR-E LST and in-situ measured LST are from 5.1 K to 6.0 K in the daytime and 3.1 K to 4.5 K in the nighttime. Further comparison with the in-situ measured air temperatures at 24 meteorological stations confirms the good performance of the TL-LUT method. The feasibility of PMW remote sensing in estimating the LST for China can complement the TIR data and can, therefore, aid in the generation of daily LST maps for the entire country. Further study of the penetration of PMW radiation would benefit the LST estimations in barren and other sparsely vegetated environments.     

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.     

光谱指数趋势面的城市地表温度降尺度转换

针对城市及周边区域建造区和自然地表交织分布的特点,探讨了利用归一化植被指数(NDVI)和归一化建造指数(NDBI)构造趋势面的地表温度(LST)降尺度方法,以北京市市区及周边较平坦区域为例实现了LST自960 m向120 m的降尺度转换。分析了LST空间分布特征及NDVI、NDBI对地物的指示性特征;以北京市四至六环为界分析NDVI、NDBI趋势面对地表温度的拟合程度及各自的适用区域;在120 m、240 m、480 m和960 m 4个尺度上评价了NDVI、NDBI和NDVI+NDBI趋势面对LST的拟合程度和趋势面转换函数的尺度效应;对NDVI、NDBI和NDVI NDBI等3种方法的降尺度结果分覆盖类型、分区域对比评价。实验结果表明结合两种光谱指数的NDVI NDBI方法降尺度转换精度有所改善,改善程度取决于地表覆盖类型组合。     

利用HJ-1B模拟数据反演地表温度的两种单通道算法

针对HJ-1B热红外波段特点,采用修正型QK & B算法,反演广州市2013-01-14的地表温度（land surface temperature,LST）。建立偏微分方程得出,当辐射率误差为0.01时,引起的LST误差约为0.6 K,LST误差与大气透过率成反比,与大气透过率误差成正比,0.1的透过率误差引起LST误差约1 K。大气水汽含量w误差与LST误差成线性关系,当大气水汽含量误差为0.1 g/cm²时,引起LST误差约为0.2 K。LST反演误差与近地表气温误差和大气平均作用温度误差均成正比,1 K的近地表气温误差引起LST反演误差约1K。总的来说,LST反演误差与区间比值和大气平均作用温度误差和近地表气温误差相关。用算法反演出来的广州市地表温度与MOD11_L2温度产品具有较强的空间一致性,温度差值曲线呈正态分布,主要集中在-0.9~0.9℃区域,选取广州市6个观测点,得出修正型QK & B算法和实测地温平均值相差约为0.31 K,MOD11_L2与实测地温的温度平均值相差0.65 K,误差均小于1 K。通过对修正型QK & B算法偏微分方程的推导,可对HJ-1B/IRS中的LST反演进行更细致和精确的分析,为其他针对环境卫星热红外波段类似反演LST的算法提供一定的借鉴,也为后续提高LST反演精度提供科学依据。     

Assessment of Vegetation Response to Drought in Nebraska Using Terra-MODIS Land Surface Temperature and Normalized Difference Vegetation Index

Eight-day composite Terra-MODIS cumulative LST and NDVI timeseries data were used to analyze the responses of crop and grassland cover types to drought in Nebraska. Four hundred ninety 1 km pixels that included irrigated and non-irrigated corn and soybeans and three grassland cover types were selected across the state of Nebraska. Statistical analyses revealed that the majority of the land cover pixels experienced significantly higher daytime and nighttime LSTs and lower NDVI during the drought-year growing season (p < 0.01). Among the land cover types analyzed, grassland experienced the highest increase in daytime LST and decrease in NDVI.