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

作为驱动地表与大气之间能量交换的关键物理量,地表温度在众多领域中都发挥着重要作用,包括气候变化、环境监测、蒸散发估算以及地热异常勘探等.Landsat热红外数据因其时间连续性和高空间分辨率等特点被广泛应用于地表温度反演中.本文详细地介绍了Landsat热红外传感器及其可用的数据与产品的现状,梳理了2001年-2020年...     

从ASTER和Landsat／TM卫星数据反演日本千叶地区地表反射率和气溶胶光学厚度分布

在卫星遥感大气研究中，已有的精确反演海洋表面反射率和其上空气溶胶光学厚度分布的算法，由于陆地地表像元反射率的不均一性，使得这些方法在应用于陆地表面反射率的反演中具有一定的局限性。而利用三步校正法可以去除地表邻近像元的影响，从而消除这种局限性。文中介绍了三步校正的算法，并对日本千叶地区ASTER卫星数据进行了大气和地表邻近像元的影响校正，精确地反演了该地区的地表反射率，通过地表反射率和气溶胶光学厚度之间的关系计算得到了气溶胶光学厚度的分布。同时证明了在洁净天反演的地表反射率可以应用于反演同一季节中非洁净天的气溶胶光学厚度分布，这样不但减小了气溶胶模式选择的影响，而且实现了在缺少太阳辐射计数据情况下获得气溶胶光学厚度分布的目的。     

Characterizing bi-temporal patterns of land surface temperature using landscape metrics based on sub-pixel classifications from Landsat TM/ETM+

Landscape patterns in a region have different sizes, shapes and spatial arrangements, which contribute to the spatial heterogeneity of the landscape and are linked to the distinct behavior of thermal environments. There is a lack of research generating landscape metrics from discretized percent impervious surface area data (ISA), which can be used as an indicator of urban spatial structure and level of development, and quantitatively characterizing the spatial patterns of landscapes and land surface temperatures (LST). In this study, linear spectral mixture analysis (LSMA) is used to derive sub-pixel ISA. Continuous fractional cover thresholds are used to discretize percent ISA into different categories related to urban land cover patterns. Landscape metrics are calculated based on different ISA categories and used to quantify urban landscape patterns and LST configurations. The characteristics of LST and percent ISA are quantified by landscape metrics such as indices of patch density, aggregation, connectedness, shape and shape complexity. The urban thermal intensity is also analyzed based on percent ISA. The results indicate that landscape metrics are sensitive to the variation of pixel values of fractional ISA, and the integration of LST, LSMA. Landscape metrics provide a quantitative method for describing the spatial distribution and seasonal variation in urban thermal patterns in response to associated urban land cover patterns.     

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.     

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.     

两种不同尺度下地表温度之间的转换——以艾比湖流域精河县为例

针对Landsat数据收集易受到传感器及外部环境的影响,数据成像周期长,不能获取所需要的数据;而中分辨率成像光谱仪数据具有较高的波谱分辨率、时间分辨率和免费获取,适合于进行大尺度、周期性的研究等优点,该文提出了建立Landsat反演的温度数据与中分辨率成像光谱仪温度产品两种数据间的关系,以便两种数据间的相互转换及验证.采用单窗算法对2011年9月13日Landsat热红外数据进行了地表温度反演,并将其与同一天MOD11A1数据中的LST_day_1km数据进行了比较,选取典型地表类型进行讨论.结果表明在所选的典型地物中,绿地(耕地、林地)温度的相关性最高,地表覆盖单一,像元的混合度较低;盐漠(包括盐渍地、荒漠等)温度的相关性较低,其地处绿洲与沙漠过渡带,地表覆盖相较于绿地复杂,像元的混合程度较高;湖面温度的相关性最低,湖面像元受到湖水深度、影像分辨率等的影响,其温度受到很大影响.该文所建立的函数关系,可为今后利用中分辨率成像光谱仪温度产品数据推导出同一时间上热混合影像中的地表温度提供参考.     

中国区域植树造林对地表温度的影响

近些年随着中国植树造林计划的持续进行,植被覆盖不断增长,进而影响中国地表环境.地表温度是表征地表物理过程的重要参数,是地表—大气能量交换中的驱动因子,广泛应用于气候、水文、生态及气象等研究中,是众多基础学科研究的关键参数之一.为了分析植树造林战略对局地和区域尺度的影响,以指导相应政策的制定,本文采用IBM(Intrin...     

Change detection of surface mining activity and reclamation based on a machine learning approach of multi-temporal Landsat TM imagery

Being able to quantify land cover changes due to mining and reclamation at a watershed scale is of critical importance in managing and assessing their potential impacts to the Earth system. In this study, a remote sensing-based methodology is proposed for quantifying the impact of surface mining activity and reclamation from a watershed to local scale. The method is based on a Support Vector Machines (SVMs) classifier combined with multi-temporal change detection of Landsat TM imagery. The performance of the technique was evaluated at selected open mining sites located in the island of Milos in Greece. Assessment of the mining impact in the studied areas was based on the confusion matrix statistics, supported by co-orbital QuickBird-2 very high spatial resolution imagery. Overall classification accuracy of the thematic land cover maps produced was reported over 90%. Our analysis showed expansion of mining activity throughout the whole 23-year study period, while the transition of mining areas to soil and vegetation was evident in varying rates. Our results evidenced the ability of the method under investigation in deriving highly and accurate land cover change maps, able to identify the mining areas as well as those in which excavation was replaced by natural vegetation. All in all, the proposed technique showed considerable promise towards the support of a sustainable environmental development and prudent resource management.     

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

地表温度与发射率是地表—大气系统长波辐射和潜热通量交换的直接驱动力,是描述区域和全球尺度上地表能量平衡与水平衡的重要参数,其时空变化信息在气象预测、气候变化、水循环、地质勘探、农林监测和城市热环境等诸多领域具有广泛的应用.热红外遥感作为当前获取区域或全球尺度上地表温度和发射率的最有效手段之一,相较于传统的地面点位测量方...     

Landsat 8热红外数据定标参数的变化及其对地表温度反演的影响

Landsat系列卫星上的TIRS热红外传感器数据已被大量应用,针对TIRS数据的地表温度反演也相继开发出一些算法,并有一些研究对TIRS数据的定标及其地表温度反演算法的精度进行了对比。本文主要就TIRS热红外传感器定标参数的变化,结合这些定标参数变化的时间点对有关地表温度反演算法的适用性和有效性进行分析,特别是对劈窗算法是否适合当前的TIRS数据进行了讨论,以使用户能够对Landsat 8 TIRS热红外数据的正确使用有进一步的认识。总的看来,由于视域外杂散光的影响,TIRS数据的定标精度仍达不到设计目标,TIRS第11波段的不确定性仍成倍大于TIRS 10波段。因此,在Landsat团队未彻底解决这一问题之前,同时用TIRS第10、第11这两个差距较大的波段构成的劈窗算法来反演地表温度,其精度存在较大的不确定性,US6-S团队仍在致力于改进第11波段的精度,改进后的波段可以用劈窗策法。目前应以TIRS第10单波段的方式来反演地表温度为宜。     

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.     

Landsat 8地表温度反演及验证—以黑河流域为例

地表温度是区域和全球尺度地表物理过程的一个重要参数,目前已有的地表温度产品空间分辨率较低,缺乏高空间分辨率的地表温度产品。Landsat系列卫星提供了大量免费的高空间分辨率遥感数据,然而对应的高空间分辨率地表温度产品还未见到,为了获取长时间序列的高空间分辨率地表温度数据,针对Landsat 8 TIRS数据提出了一个物理单通道地表温度反演算法。该算法首先利用ASTER全球地表发射率产品(ASTER GED)结合Landsat 8地表反射率产品计算Landsat 8影像的地表发射率,然后利用快速辐射传输模型RTTOV结合MERRA大气廓线数据对热红外影像进行大气校正,最后利用物理单通道地表温度反演算法得到地表温度。利用黑河流域HiWATER试验2013年—2015年15个站点的实测地表温度数据对本文方法和普适性单通道算法进行了验证,同时对验证站点的空间异质性进行了分析。结果表明,本文方法和普适性单通道算法估算的地表温度整体精度均较高,能够获取高精度、高空间分辨率的地表温度数据,可以服务于城市热岛效应、地表蒸散发估算等相关研究。     

应用天宫一号高光谱红外谱段数据反演地表温度产品及验证

本文利用天宫一号获取的具有高空间分辨率的红外谱段数据,选择黑河中游作为反演区域,基于辐射传输模型和Planck定律进行地表温度反演.影像的反演结果表明使用天宫一号高空间分辨率的红外谱段反演得到的地表温度能更细致合理的刻画地表温度分布,这对地表温度异质性的研究具有重要意义.为了进一步验证反演精度,本文依托黑河流域已有观测平台和系统,选取黑河中游8个典型下垫面作为验证场地,开展地面同步观测实验.实验结果表明通过反演获得的地表温度与地面实测数据的平均偏差为-0.375℃,这也说明本文的反演过程能够较为精确的使用天宫一号红外谱段反演地表温度.     

城市天空可视因子对地表热环境的影响分析

针对城市热岛效应的各类影响因素,该文主要关注城市三维形态与城市热环境的关系,并研究这种关系因建筑高度、密度、季节而变化的原因.作为重要的城市结构参数,天空可视因子能够较好地描述建筑物的三维形态,并可以表征有效太阳辐射吸收强度.因此,该文结合建筑物高度数据集和半球模型计算天空可视因子,并利用Landsat-8热红外波段和辐射传输模型反演地表温度,然后分析了天空可视因子与地表温度在不同季节的相关性.选择地表异质性较高的北京核心建成区作为研究区,并按照不同建筑高度和密度选择4个实验区.实验结果表明,天空可视因子与地表温度的关系存在季节差异,冬季两者的相关性更强.通过对4个实验区的定量分析可知,实验区a中天空可视因子与地表温度夏季是正相关,冬季是非线性关系;实验区b和c中是两者的关系正相关;而实验区d中无明显相关趋势.所以,建筑物的高度与密度对天空可视因子与地表温度的关系有显著影响,并且这种关系随季节变化而不同.     

The impact of the dynamic sea surface topography on the quasi-geoid in shallow coastal waters

In this study, we examine the impact of instantaneous dynamic sea surface topography (DT) corrections to be applied to altimeter-derived sea surface slopes on the quasi-geoid in the shallow and coastal waters of the North Sea. In particular, we investigate the added value of DT corrections obtained from a shallow-water hydrodynamic model. These corrections comprise the contributions of ocean tides, wind- and pressure-driven (surge), and density-driven (baroclinic) water-level variations including the interactions between them. As a reference, we used tidal corrections derived from the global ocean tide model GOT4.7, surge corrections derived from the MOG2D model, and corrections for the time-averaged baroclinic contribution computed as differences between the DTU10 mean sea surface model and the EGG08 quasi-geoid. From a spectral analysis, we found that the baroclinic and surge parts of the DT mainly contribute to improvements in the signal-to-noise ratio (SNR) at longer wavelengths down to $100{-}200~\hbox {km}$ and that the improvements increase towards the southern North Sea. We also found that the shallow-water hydrodynamic model provides better tidal corrections compared to the GOT4.7 global ocean tide model, which are most pronounced in the southern North Sea and affect almost the entire spectrum. Very small differences (mostly below ${\pm } 2~\hbox {cm}$ ) are observed between the quasi-geoid solutions obtained using the different sets of DT corrections. We showed that the variance component estimation provides too optimistic variance factors for the shipboard data set relative to the altimeter-derived quasi-geoid slopes. Hence, the limited impact of DT corrections is due to the fact that altimeter-derived quasi-geoid slopes hardly contribute to the estimated quasi-geoid if shipboard gravity data are included. When computing quasi-geoid solutions without shipboard gravity data, we found that less accurate or incomplete DT corrections may cause errors in the quasi-geoid with systematic spatial patterns. These systematic patterns disappear or are reduced significantly when using the DT corrections provided by the shallow-water hydrodynamic model. The main contributor to this improvement is the better tidal correction provided by the shallow-water hydrodynamic model compared to the GOT4.7 global ocean tide model. Seen the improvements of the global ocean tide models over the last two decades, we expect that in the near future global ocean tide models perform as well as dedicated regional models such as DCSM. Critical issue is, however, access to high-quality local bathymetric data.     

HJ-1B热红外辐射定标对地表温度反演的影响

辐射定标是地表温度反演的关键环节。在实际应用中,定标参数的选取直接影响辐射定标结果。由于卫星在轨定标系统的灵敏度会随时间发生改变,这使得实际的辐射定标存在不确定性。本文基于单波段热红外辐射传输方程给出了辐射定标影响地表温度反演的敏感性分析模型,并重点对HJ-1B热红外波段的辐射定标问题及其对地表温度反演的影响进行讨论和分析。敏感性分析模型显示,在一般的情形下,辐射定标偏差与地表温度反演误差的比值在数值上约为1∶11,即0.1个单位(W·m–2·sr–1·μm–1)的辐射定标偏差可能引起1.1 K左右的温度反演误差。由于对HJ-1B热红外波段定标参数的更新存在滞后,即使仅考虑相邻年份间的定标偏差,相应的地表温度平均绝对误差也在1.0—2.0 K左右。在更严重的情形下,因辐射定标偏差引起的地表温度误差甚至高达5.0 K以上。同时,分析结果也表明对HJ-1B的存档影像,无论使用其头文件中的定标参数还是当年公布的定标参数,定标后辐射亮度值的偏差均存在一定的不稳定性。针对当前在HJ-1B定标参数更新方面的不足,本文假定传感器辐射性能变化是匀速稳定的,相应地提出了线性内插修正法和线性外推法两个定标参数估计/修正方案。通过线性内插的方法进行辐射定标修正后,在一定程度上减小了误差的不稳定性,与MODIS B31波段的亮温相比,偏差在1.0 K以内(约为0.5 K)。当定标参数尚未公布时,由线性外推法修正后的辐射定标结果所得亮温的偏差约为1.0 K。与在大亚湾核电站附近海域进行的一次卫星同步观测(2011年12月18日)相比,经两个方案修正后的海表温度偏差分别约为0.2 K和0.7 K。案例分析的结果表明,在现阶段,线性内插修正法和线性外推法是简单有效的,可供用户在实际应用中参考和使用。     

Precise surface topography of the grounded ice ridges at the Ekströmisen, Antarctica, based on several geophysical data sets

As part of the CryoSat Cal/Val activities and the pre-site survey for an ice core drilling contributing to the International Partnerships in Ice Core Sciences (IPICS), ground-based kinematic GPS measurements were conducted in early 2007 in the hinterland of the German overwintering station Neumayer (8.25 W,70.65 S). The investigated area comprises the regions of the ice ridges Halvfarryggen and Søråsen, which rise from the Ekströmisen to a maximum of about 760 m surface elevation, and have an areal extent of about 100 km×50 km each. Available digital elevation models (DEMs) from radar altimetry and the Antarctic Digital Database show elevation differences of up to hundreds of meters in this region, which necessitated an accurate survey of the conditions on-site. An improved DEM of the Ekströmisen surroundings is derived by a combination of highly accurate ground-based GPS measurements, satellite derived laser altimetry data (ICESat), airborne radar altimetry (ARA), and radio echo sounding (RES). The DEM presented here achieves a vertical accuracy of about 1.3 m and can be used for improved ice dynamic modelling and mass balance studies.     

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蚀变异常信息提取

遥感蚀变异常信息在一定程度上可以反映出围岩蚀变情况.为了利用遥感数据提取围岩蚀变异常信息,通过分析蚀变矿物的波谱曲线,首先得出其在ASTER各波段的吸收、反射特征;然后选择特征明显的波段做主成分分析,并依据其光谱特征与特征向量的对应关系,确定出主要包含围岩蚀变信息的主分量;最后对该分量用分形模型计算其异常灰度的阚值,来...