一个从ASTER数据中反演地表温度的劈窗算法

根据EOS/Terra多传感器的特点,提出了一个适合于ASTER数据的劈窗算法,该算法包括两个必要的参数大气透过率和比辐射率。大气透过率是通过利用MODIS的3个近红外波段反演大气水汽含量并根据大气水汽含量与热红外波段的统计关系计算得到。由于MODIS和ASTER是在同一颗星上,这种大气透过率估计方法保证了地表温度反演过程中所需大气参数的同步获取。对于比辐射率则是通过分类和JPL提高的光谱库获得。最后用大气模拟校正法对算法进行了验证,在比辐射率已知的情况下,当使用大气模型模拟得到的大气透过率时,对Planck函数优化简化后的平均精度为0.56℃;当大气透过率是从大气水汽含量计算得到时,优化平均精度为0.58℃,表明该算法可行。     

MODIS数据反演地表温度的参数敏感性分析

在利用MODIS卫星遥感数据进行地表温度反演过程中，有两个基本参数需要确定，即地表比辐射率和大气透过率，尽管采用了比较合理的参数估计方法，但仍会有一些不可避免的因素导致误差的产生。为了进一步研究可能的参数误差对地表温度反演精度的影响，我们对该算法的两个参数进行敏感性分析。结果表明，当31、32两个波段的参数估计都有中等误差时，可能的地表温度误差对大气透过率和地表比辐射率都不敏感，所引起的地表温度误差大约为0．6～0．8℃，算法能够得到较高精度的地表温度反演结果。     

陆地资源卫星数据地表温度反演

针对目前陆地资源卫星(Landsat-8)地表温度反演过程中,地表比辐射率估计和敏感度分析中存在的不足,对这两方面进行改进,提出了一种基于Landsat-8数据的地表温度反演算法。该文主要从劈窗算法的推导、参数的估计、敏感度分析等方面进行研究。对于大气透过率的计算,首先用与其有相邻过境时间的MODIS数据反演大气水汽含量,然后通过中分辨率的大气传输模型(Moderate Resolution Atmospheric Transmission,MODTRAN)模拟大气水汽含量与透过率的关系,最后得到大气透过率。对于发射率的计算,通过分类和ASTER提供的光谱库获得。将大气辐射传输方程模拟的地表温度与此劈窗算法反演的地表温度做比较,结果表明平均精度达到0.82K。最后研究了大气水汽含量对地表温度的影响。结果显示,当大气水汽含量误差为0.1g/cm2,其对温度反演精度的影响最大不超过0.3K;当大气水汽含量的反演误差较大的时候,其对温度反演精度的影响较大。     

基于MODIS 影像数据的劈窗算法研究及其参数确定

劈窗算法是目前由热红外遥感数据获取陆面温度的主要方法。在介绍劈窗算法的一般表现形式的基础上, 我们推导出适合于MOD IS 影像数据的劈窗算法。大气透过率和地表比辐射率是求解地表温度的两个关键参数。由于MOD IS 图像分辨率较低,MOD IS 像元主要由水面、植被和裸土3种地物类型构成, 故可依据这3 种地物的构成比例确定地表比辐射率。从遥感影像上反演大气的水汽含量, 再根据大气水汽含量与大气透过率的关系计算出大气透过率。最后将文中推导的劈窗算法用于江苏省地表温度的反演。反演出来的地表温度图显示出明显的地表温度空间差异、城市热岛效应和不同的地物类型。     

基于MODIS 数据的南京市夏季城市热岛分析

城市热岛效应是当前城市环境与气候主要研究内容之一。地表温度与气温之间有紧密的联系, 通过遥感反演地表温度已成为研究城市热岛的有效手段。利用MODIS 数据, 获取地表比辐射率与大气透过率2 个基本参数, 运用劈窗算法反演南京市夏季地表温度。基于不同时相的MODIS数据, 对4 幅南京市地表温度反演图像作对比分析, 较好地显示了南京市城市热岛的空间分布、热岛范围和城市热岛强度, 结果表明南京市夏季热岛问题较为严重。     

MODIS数据反演地表温度的传感器视角校正研究

提出了MODIS数据传感器视角的计算方法,并分析了视角的变化对大气透过率以及地表温度反演的影响。选取渤海地区一景影像进行实验研究,结果表明MODIS数据边缘的传感器视角可达55.02°,由此引起的大气透过率降低近0.086,引起的地表温度误差最高可达3.64℃。     

基于分组分裂窗算法的MTSAT-1R地表温度反演

以黑河流域上游和中游为研究区,针对MTSAT-1R卫星数据,运用MODTRAN 4.0及晴空状态下的TIGR大气廓线数据,发展了根据地表比辐射率、大气水汽含量、传感器观测角度分组模拟的分裂窗算法,进行地表温度反演。分析了传感器噪声、地表比辐射率和大气水汽含量3个参数对该算法的影响,并结合模拟数据、地面观测数据及MODIS地表温度产品,对反演结果进行分析评价。结果表明:当传感器垂直观测或大气水汽含量小于2.5g/cm2时,反演精度在1K以内;反演结果与地面观测数据对比差异较小,在阿柔站RMSE为3.7 K(日)/1.4 K(夜),在盈科站RMSE为2.4K(日)/2.0K(夜);与MODIS地表温度产品比较,空间分布呈现出一致性。总之,分组分裂窗算法能较好地用于MTSAT-1R卫星数据进行地表温度反演。     

一个从MODIS数据同时反演地表温度和发射率的神经网络算法

MODIS的三个热红外波段29、31、32建立了三个辐射传输方程,这三个方程包含了5个未知数(大气平均作用温度、地表温度和三个波段的发射率)。用JPL提供的大约160种地物的波谱数据对MODIS三个波段(29/31/32)发射率之间的关系和用MODTRAN4对大气透过率和大气水汽含量之间关系进行模拟分析。分析结果表明地球物理参数之间存在着大量的潜在信息。由于潜在的信息难以严格地用数学表达式来描述,因此神经网络是非常适合被用来解这种病态反演问题。利用辐射传输模型(RM)和神经网络(NN)反演分析表明神经网络能够被用来精确地同时从MODIS数据中反演地表温度和发射率。地表温度的平均反演误差在0.4°C以下;波段29/31/32发射率平均反演误差都在0.008以下。     

中国江淮、黄淮地区陆面微波比辐射率的变化特征

陆面微波比辐射率较高且易变,造成陆面上反演降水以及其它大气参数较为困难。对于地表特征复杂的中国,陆面微波比辐射率的研究还很有限。通过利用Tropical Rainfall Measuring Mission（TRMM）卫星上同步扫描的VIRS（红外和可见光）与TMI（微波）资料以及微波辐射传输模式反演了中国江淮、黄淮地区陆面微波比辐射率。然后,结合MODIS提供的地表类型数据,分析了江淮、黄淮地区不同地表微波比辐射率的时空变化特征。 结果表明该地区的农作物地表比辐射率最小,垂直与水平比辐射率极化差最大;而森林地表比辐射率最大,极化差最小。此外,不同地表的微波比辐射率昼夜变化明显,季节变化不明显。比辐射率估算误差中,地表温度、微波亮温和大气相对湿度3因子的准确计算对22GHz和85GHz的影响较为明显,对其它通道影响较小。对于小于85GHZ的通道,比辐射率估算精度受微波亮温的影响最为明显,地表温度其次,相对湿度最小;对于高频85OHz,相对湿度的影响最明显,其次是微波亮温,最后是地表温度。     

中国江淮、黄淮地区陆面微波比辐射率的变化特征

陆面微波比辐射率较高且易变,造成陆面上反演降水以及其它大气参数较为困难。对于地表特征复杂的中国,陆面微波比辐射率的研究还很有限。通过利用Tropical Rainfall Measuring Mission (TRMM)卫星上同步扫描的VIRS（红外和可见光）与TMI（微波）资料以及微波辐射传输模式反演了中国江淮、黄淮地区陆面微波比辐射率。然后,结合MODIS提供的地表类型数据,分析了江淮、黄淮地区不同地表微波比辐射率的时空变化特征。 结果表明该地区的农作物地表比辐射率最小,垂直与水平比辐射率极化差最大;而森林地表比辐射率最大,极化差最小。此外,不同地表的微波比辐射率昼夜变化明显,季节变化不明显。比辐射率估算误差中,地表温度、微波亮温和大气相对湿度3因子的准确计算对22 GHz和85 GHz的影响较为明显,对其它通道影响较小。对于小于85 GHz的通道,比辐射率估算精度受微波亮温的影响最为明显,地表温度其次,相对湿度最小;对于高频85 GHz,相对湿度的影响最明显,其次是微波亮温,最后是地表温度。     

A practical split‐window algorithm for retrieving land‐surface temperature from MODIS data

This paper presents a practical split‐window algorithm utilized to retrieve land‐surface temperature (LST) from Moderate‐resolution Imaging Spectroradiometer (MODIS) data, which involves two essential parameters (transmittance and emissivity), and a new method to simplify Planck function has been proposed. The method for linearization of Planck function, how to obtain atmosphere transmittance from MODIS near‐infrared (NIR) bands and the method for estimating of emissivity of ground are discussed with details. Sensitivity analysis of the algorithm has been performed for the evaluation of probable LST estimation error due to the possible errors in water content and emissivity. Analysis indicates that the algorithm is not sensitive to these two parameters. Especially, the average LST error is changed between 0.19–1.1°C when the water content error in the simulation standard atmosphere changes between ?80 and 130%. We confirm the conclusion by retrieving LST from MODIS image data through changing retrieval water content error. Two methods have been used to validate the proposed algorithm. Results from validation and comparison using the standard atmospheric simulation and the comparison with the MODIS LST product demonstrate the applicability of the algorithm. Validation with standard atmospheric simulation indicates that this algorithm can achieve the average accuracy of this algorithm is about 0.32°C in LST retrieval for the case without error in both transmittance and emissivity estimations. The accuracy of this algorithm is about 0.37°C and 0.49°C respectively when the transmittance is computed from the simulation water content by exponent fit and linear fit respectively.     

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.     

Evaluating NDVI-based emissivities of MODIS bands 31 and 32 using emissivities derived by Day/Night LST algorithm

Surface emissivity estimation is a significant factor for the land surface temperature estimation from remotely sensed data. For fully vegetated surfaces, the emissivity estimation is performed in a simple manner since the emissivity is relatively uniform. However, for arid land with sparse vegetation, the estimation is more complicated since the emissivity of the exposed soil and rock is highly variable. In this study, mean and difference emissivity for bands 31 and 32 of MODIS sensor have been derived based on NDVI values. First, the NDVI thresholds have been determined to separate bare soil, partially vegetated soil and fully vegetated land. Then regression relations have been derived to estimate mean and difference emissivity of the bare soil samples and partially vegetated surfaces. A constant emissivity is also used for fully vegetated area. Along with the correlations, standard deviations of the regression relations have been examined for a set of representative soil types. Standard deviations smaller than 0.003 in mean emissivity and smaller than 0.004 in difference emissivity are resulted in regression linear relations. Evaluation of the NDVI derived regression relations has been performed using the results of MODIS Day/Night Land Surface Temperature (LST) algorithm on a pair of MODIS images. Using around 45,500 pixels with different soil and land cover types, emissivity of each pixel in bands 31 and 32 have been estimated. The calculated emissivities have been compared with emissivities calculated by MODIS Day/Night LST algorithm. Biases and standard deviations of NDVI-based relations show relatively high agreement for mean and difference emissivity relations with Day/Night method results. It may be concluded that the proposed algorithm can be used as a rather simple alternative to complex emissivity estimation algorithms.     

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

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

HJ-1B卫星遥感数据在重庆市城市热环境监测的应用分析

本文以MODIS反演大气透射率,以HJ-1B/CCD分类结果反演地表比辐射率,并基于单窗算法,利用HJ-1B/IRS4数据反演地表温度.在此基础上,提取研究区的热场变异指数来分析重庆热岛空间分布特征,并就NDVI与NDBI对热岛效应的影响进行了分析.其结果如下:1)重庆城市热岛大致位于中梁山、铜锣山之间,呈东北、西南走向分布;2)热岛中心不在市中心,而是集中在大渡口工业园区、江北机场这些能耗大、人口密集区域,热岛强度范围在5?C-10?C之间;3)接近长江、嘉陵江水域的建筑用地密集区域,其热岛效应并不明显;4)NDVI与热岛强度呈负相关关系,NDBI与热岛强度呈现较为明显的正相关关系,二者对热岛都有重要影响,而NDBI的影响更大.因此,利用HJ-1B数据监测城市热环境,能较好地揭示重庆城市热岛空间分布特征,为城市环境监测与改善提供参考.     

A generalized split-window algorithm for land surface temperature estimation from MSG-2/SEVIRI data

This paper aims to determine land surface temperature (LST) using data from a spinning enhanced visible and infrared imager (SEVIRI) on board Meteosat Second Generation 2 (MSG-2) by using the generalized split-window (GSW) algorithm. Coefficients in the GSW algorithm are pre-determined for several overlapping sub-ranges of the LST, land surface emissivity (LSE), and atmospheric water vapour content (WVC) using the data simulated with the atmospheric radiative transfer model MODTRAN 4.0 under various surface and atmospheric conditions for 11 view zenith angles (VZAs) ranging from 0° to 67°. The results show that the root mean square error (RMSE) varies with VZA and atmospheric WVC and that the RMSEs are within 1.0 K for the sub-ranges in which the VZA is less than 30° and the atmospheric WVC is less than 4.25 g cm^?2. A sensitivity analysis of LSE uncertainty, atmospheric WVC uncertainty, and instrumental noise (NEΔT) is also performed, and the results demonstrate that LSE uncertainty can result in a larger LST error than other uncertainties and that the total error for the LST is approximately 1.21 and 1.45 K for dry atmosphere and 0.86 and 2.91 K for wet atmosphere at VZA = 0° and at VZA = 67°, respectively, if the uncertainty in the LSE is 1% and that in the WVC is 20%. The GSW algorithm is then applied to the MSG-2 – SEVIRI data with the LSE determined using the temperature-independent spectral indices method and the WVC either determined using the measurements in two split-window channels or interpolated temporally and spatially using European Centre for Medium Range Weather Forecasting (ECMWF) data. Finally, the SEVIRI LST derived in this paper (SEVIRI LST1) is evaluated through comparisons with the SEVIRI LST provided by the land surface analysis satellite applications facility (LSA SAF) (SEVIRI LST2) and the Moderate Resolution Imaging Spectroradiometer (MODIS) LST product (MOD11B1 LST product). The results show that more than 80% of the differences between SEVIRI LST1 and SEVIRI LST2 are within 2 K, and approximately 70% of the differences between SEVIRI LST1 and MODIS LST are within 4 K. Furthermore, compared to MODIS LST, for four specific areas with different land surfaces, our GSW algorithm overestimates the LST by up to 1.0 K for vegetated surfaces and by 1.3 K for bare soil.     

Land surface temperature and emissivity retrieval by integrating MODIS data onboard Terra and Aqua satellites

This paper presents an algorithm to retrieve land surface temperature (LST) and emissivity by integrating MODIS (Moderate Resolution Imaging Spectroradiometer) data onboard Terra and Aqua satellites. For a study area, there will be four pairs of day and night observations by MODIS onboard two satellites every day. Solar zenith angle, view zenith angle, and atmospheric water vapour have first been taken as independent variables to analyse their sensitivities to the same infrared channel measurements of MODIS on both Terra and Aqua satellites. Owing to their similar influences on the same MODIS band from Terra and Aqua satellites, four pairs of MODIS data from Terra and Aqua satellites can be thought of as MODIS measurement on a satellite at different viewing angles and viewing time. Comparisons between the retrieved results and in-situ measurements at three test sites (Qinghai Lake, Poyang Lake and Luancheng in China) indicate that the root mean square (rms) error is 0.66 K, except for the sand in Poyang Lake area. The rms error is less than 0.7 K when the retrieved results are compared with Earth Observing System (EOS) MODIS LST data products using the physics-based day/night algorithm. Emissivities retrieved by this algorithm are well compared to EOS MODIS emissivity data products (V5). The proposed algorithm can therefore be regarded as complementary and an extension to the EOS physics-based day/night algorithm.