Validation of non-linear split window algorithm for land surface temperature estimation using Sentinel-3 satellite imagery: Case study; Tehran Province,Iran

In recent years, land surface temperature (LST) has become critical in environmental studies and earth science. Remote sensing technology enables spatiotemporal monitoring of this parameter on large scales. This parameter can be estimated by satellite images with at least one thermal band. Sentinel-3 SLSTR data provide LST products with a spatial resolution of 1 km. In this research, direct and indirect validation procedures were employed to evaluate the Sentinel-3 SLSTR LST products over the study area in different seasons from 2018 to 2019. The validation method was based on the absolute (direct) evaluation of this product with field data and comparison (indirect) evaluation with the MODIS LST product and the estimated LST using the non-linear split-window (NSW) algorithm. Also, two emissivity estimation methods, (1) NDVI thresholding method (NDVI-THM) and (2) classification-based emissivity method (CBEM), were used to estimate the LST using the NSW method according to the two thermal bands of Sentinel-3 images. Then, the accuracy of these methods in estimating LST was evaluated using field data and temporal changes of vegetation, which the NDVI-THM method generated better results. For indirect evaluation between the Sentinel-3 LST product, MODIS LST product, and LST estimated using NSW, four filters based on spatial and temporal separates between pairs of pixels and pixel quality were used to ensure the accuracy and consistency of the compared pairs of a pixel. In general, the accuracy results of the LST products of MODIS and Sentinel-3, and LST estimated using NSW showed a similar trend for LST changes during the seasons. With respect to the two absolute and comparative validations for the Sentinel-3 LST products, summer with the highest values of bias (?1.24 K), standard deviation (StDv = 2.66 K), and RMSE (2.43 K), and winter with the lowest ones (bias of 0.14 K, StDv of 1.13 K, and RMSE of 1.12 K) provided the worst and best results for the seasons in the period of 2018–2019, respectively. According to both absolute and comparative evaluation results, the Sentinel-3 SLSTR LST products provided reliable results for all seasons on a large temporal and spatial scale over our studied area.     

Determination of land surface temperature using precipitable water based Split-Window and Artificial Neural Network in Turkey

Land surface temperature (LST) calculation utilizing satellite thermal images is very difficult due to the great temporal variance of atmospheric water vapor in the atmosphere which strongly affects the thermal radiance incoming to satellite sensors. In this study, Split-Window (SW) and Radial Basis Function (RBF) methods were utilized for prediction of LST using precipitable water for Turkey. Coll 94 Split-Window algorithm was modified using regional precipitable water values estimated from upper-air Radiosond observations for the years 1990–2007 and Local Split-Window (LSW) algorithms were generated for the study area. Using local algorithms and Advanced Very High Resolution Radiometer (AVHRR) data, monthly mean daily sum LST values were calculated. In RBF method latitude, longitude, altitude, surface emissivity, sun shine duration and precipitable water values were used as input variables of the structure. Correlation coefficients between estimated and measured LST values were obtained as 99.23% (for RBF) and 94.48% (for LSW) at 00:00 UTC and 92.77% (for RBF) and 89.98% (for LSW) at 12:00 UTC. These meaningful statistical results suggest that RBF and LSW methods could be used for LST calculation.     

Simulation of the land surface temperature from moon-based Earth observations

The land surface temperature (LST) is a key parameter for the Earth’s energy balance. As a natural satellite of the Earth, the orbital of the moon differs from that of current Earth observation satellites. It is a new way to measure the land surface temperature from the moon and has many advantages compared with artificial satellites. In this paper, we present a new method for simulating the LST measured by moon-based Earth observations. Firstly, a modified land-surface diurnal temperature cycle (DTC) method is applied to obtain the global LST at the same coordinated universal time (UTC) using the Moderate Resolution Imaging Spectroradiometer (MODIS) LST products. The lunar elevation angles calculated using the ephemeris data (DE405) from the Jet Propulsion Laboratory (JPL) were then applied to simulate the Earth coverage observed from the moon. At the same time, the modified DTC model was validated using in situ data, MODIS LST products, and the FengYun-2F (FY-2F) LST, respectively. The results show that the fitting accuracy (root-mean-square error, RMSE) of the modified DTC model is not greater than 0.72?°C for eight in situ stations with different land cover types, and the maximum fitting RMSE of the modified model is smaller than that of current DTC models. By the comparison of the simulated LST with MODIS and FY-2F LST products, the errors of the results were feasible and accredited, and the simulated global LST has a reasonable spatiotemporal distribution and change trend. The simulated LST data can therefore be used as base datasets to simulate the thermal infrared imagery from moon-based Earth observations in future research.     

Surface temperature estimation in Singhbhum Shear Zone of India using Landsat-7 ETM+ thermal infrared data

Land surface temperature (LST) is an important factor in global change studies, heat balance and as control for climate change. A comparative study of LST over parts of the Singhbhum Shear Zone in India was undertaken using various emissivity and temperature retrieval algorithms applied on visible and near infrared (VNIR), and thermal infrared (TIR) bands of high resolution Landsat-7 ETM+ imagery. LST results obtained from satellite data of October 26, 2001 and November 2, 2001 through various algorithms were validated with ground measurements collected during satellite overpass. In addition, LST products of MODIS and ASTER were compared with Landsat-7 ETM+ and ground truth data to explore the possibility of using multi-sensor approach in LST monitoring. An image-based dark object subtraction (DOS3) algorithm, which is yet to be tested for LST retrieval, was applied on VNIR bands to obtain atmospheric corrected surface reflectance images. Normalized difference vegetation index (NDVI) was estimated from VNIR reflectance image. Various surface emissivity retrieval algorithms based on NDVI and vegetation proportion were applied to ascertain emissivities of the various land cover categories in the study area in the spectral range of 10.4–12.5 μm. A minimum emissivity value of about 0.95 was observed over the reflective rock body with a maximum of about 0.99 over dense forest. A strong correlation was established between Landsat ETM+ reflectance band 3 and emissivity. Single channel based algorithms were adopted for surface radiance and brightness temperature. Finally, emissivity correction was applied on ‘brightness temperature’ to obtain LST. Estimated LST values obtained from various algorithms were compared with field ground measurements for different land cover categories. LST values obtained after using Valor’s emissivity and single channel equations were best correlated with ground truth temperature. Minimum LST is observed over dense forest as about 26 °C and maximum LST is observed over rock body of about 38 °C. The estimated LST showed that rock bodies, bare soils and built-up areas exhibit higher surface temperatures, while water bodies, agricultural croplands and dense vegetations have lower surface temperatures during the daytime. The accuracy of the estimated LST was within ±2 °C. LST comparison of ASTER and MODIS with Landsat has a maximum difference of 2 °C. Strong correlation was found between LST and spectral radiance of band 6 of Landsat-7 ETM+. Result corroborates the fact that surface temperatures over land use/land cover types are greatly influenced by the amount of vegetation present.     

Exploring the effect of COVID-19 pandemic lockdowns on urban cooling: A tale of three cities

COVID-19 pandemic has had a major impact on our society, environment and public health, in both positive and negative ways. The main aim of this study is to monitor the effect of COVID-19 pandemic lockdowns on urban cooling. To do so, satellite images of Landsat 8 for Milan and Rome in Italy, and Wuhan in China were used to look at pre-lockdown and during the lockdown. First, the surface biophysical characteristics for the pre-lockdown and within-lockdown dates of COVID-19 were calculated. Then, the land surface temperature (LST) retrieved from Landsat thermal data was normalized based on cold pixels LST and statistical parameters of normalized LST (NLST) were calculated. Thereafter, the correlation coefficient (r) between the NLST and index-based built-up index (IBI) was estimated. Finally, the surface urban heat island intensity (SUHII) of different cities on the lockdown and pre-lockdown periods was compared with each other. The mean NLST of built-up lands in Milan (from 7.71 °C to 2.32 °C), Rome (from 5.05 °C to 3.54 °C) and Wuhan (from 3.57 °C to 1.77 °C) decreased during the lockdown dates compared to pre-lockdown dates. The r (absolute value) between NLST and IBI for Milan, Rome and Wuhan decreased from 0.43, 0.41 and 0.16 in the pre-lockdown dates to 0.25, 0.24, and 0.12 during lockdown dates respectively, which shows a large decrease for all cities. Analysis of SUHI for these cities showed that SUHII during the lockdown dates compared to pre-lockdown dates decreased by 0.89 °C, 1.78 °C, and 1.07 °C respectively. The results indicated a high and substantial impact of anthropogenic activities and anthropogenic heat flux (AHF) on the SUHI due to the substantial reduction of huge anthropogenic pressure in cities. Our conclusions draw attention to the contribution of COVID-19 lockdowns (reducing the anthropogenic activities) to creating cooler cities.     

Influence of land surface parameters on the spatio-seasonal land surface temperature regime in rural West Bengal,India

The magnitude and causes of changes in the land surface temperature of rural areas have not been extensively studied. The thermal band of Landsat imagery is taken to extract winter, summer, and monsoon season land surface temperature (LST) and relate it to surface parameters over a 30-year period. From the extracted parameters constructed a prospective surface temperature (PST) model using Multivariate Adaptive Regression Splines. The Chandrabhaga river basin in West Bengal of the lateritic Rarh Tract at the Chota Nagpur Plateau fringe was chosen as the study area because it is far from urban influences, to avoid the well-known heat island effect. Over the study period, summer and winter average LST increased linearly by 0.085?°C/y and 0.016?°C/y respectively. These results were validated with air temperature (RMSE?=?x and y, respectively). Over time more of the area is in the higher temperature zones, e.g., in April 2011, 4% area exceeded >32°, whereas in 2015 this proportion reached 52%. PST models of all the seasons were moderate to highly correlate (0.57–0.87) with actual LST, showing the value of this model. It also revealed the relative importance of the regional factors. Based on this information factor management is a scientific step to restrict or minimize the temperature rise effect.     

A comparison of different regression models for downscaling Landsat and MODIS land surface temperature images over heterogeneous landscape

Remotely sensed high spatial resolution thermal images are required for various applications in natural resource management. At present, availability of high spatial resolution (<200 m) thermal images are limited. The temporal resolution of such images is also low. Whereas, coarser spatial resolution (∼1000 m) thermal images with high revisiting capability (∼1 day) are freely available. To bridge this gap, present study attempts to downscale coarser spatial resolution thermal image to finer spatial resolution using relationships between land surface temperature (LST) and vegetation indices over a heterogeneous landscape of India. Five regression based models namely (i) Disaggregation of Radiometric Temperature (DisTrad), (ii) Temperature Sharpening (TsHARP), (iii) TsHARP with local variant, (iv) Least median square regression downscaling (LMS_DS) and (v) Pace regression downscaling (PR_DS) are applied to downscale LST of Landsat Thematic Mapper (TM) and Terra MODIS (Moderate Resolution Imaging Spectroradiometer) images. All the five models are first evaluated on Landsat image aggregated to 960 m resolution and downscaled to 480 m and 240 m resolution. The downscale accuracy is achieved using LMS_DS and PR_DS models at 240 m resolution at 0.61 °C and 0.75 °C respectively. MODIS data downscaled from 1000 m to 250 m spatial resolution results root mean square error (RMSE) of 1.43 °C and 1.62 °C for LMS_DS and PR_DS models, respectively. The LMS_DS model is less sensitive to outliers in heterogeneous landscape and provides higher accuracy when compared to other models. Downscaling model is found to be suitable for agricultural and vegetated landscapes up to a spatial resolution of 250 m but not applicable to water bodies, dry river bed sand sandy open areas.     

Landsat 8地表温度产品验证与分析

Landsat系列卫星热波段具有60~120m的空间分辨率，对各种环境监测起到了重要的作用。随着Landsat系列卫星在全球范围内地表温度（land surface temperature，LST）产品的发布，其验证工作也随之展开，然而对于长时间序列的精度验证工作仍然缺乏。以黑河流域中游为研究区，利用研究区内湿地站（SD）、戈壁站（GB）和大满超级站（CJZ）三个气象站的地面测量数据对2013-2016年清晰无云的31景Landsat 8地表温度产品进行了验证与分析，并将Landsat 8地表温度产品与广泛使用的普适性单通道算法（JMS）反演结果进行了对比。结果表明，Landsat 8地表温度产品与普适性单通道算法反演结果精度均较高，在各个站点处R2均优于0.949。基于所有站点分析，Landsat 8地表温度产品精度稍高于普适性单通道算法反演结果。     

复合材料红外无损检测的建模分析和热像处理

红外热成像无损检测由于检测面积大、速度快和非接触从而在复合材料缺陷检测中的应用越来越广泛.对碳纤维增强塑料的红外热成像无损检测从数学建模和热像数据处理2个方面作了研究.针对含有teflon夹片的碳纤维增强塑料试件的红外无损检测实验中出现的缺陷表面温差信号为负的"非传统"现象提出了新的数学模型,作了仿真计算,解决了常规模型预测结果与实验结果之间的矛盾.应用5种国际上流行的热像数据处理方法对实验数据进行了处理,并对各种处理方法用统计方法作了比较,结果表明傅里叶变换和相关分析法具有最高的信噪比.      

ARIMA-SVR组合模型在卫星遥测参数预测中的应用

为辅助卫星在轨运行提供决策分析支持，结合卫星遥测参数的时间序列特性，利用一种ARIMA-SVR组合预测方法，通过对卫星遥测参数进行预测，判定实际遥测数据是否处于正常范围。该组合模型利用ARIMA模型对预处理后的数据进行线性拟合，并利用SVR模型对数据的非线性部分进行补偿。以KX09卫星星敏A的温度遥测数据为基础，分别利用组合模型对短期及中期星敏A温度进行预测，得出短期和中期均方根误差（RMSE）分别为0.768和0.968，相比单一ARIMA模型，短中期RMSE分别提高46.2%和16.4%。此外，对该卫星陀螺B的x轴角速度进行了短中期预测:短期预测中，组合模型比单一ARIMA模型的RMSE提高71.2%；中期预测中，组合模型比单一ARIMA模型的RMSE提高64.2%。实验结果表明，ARIMA-SVR组合模型为保证卫星在轨正常运行提供了有效的决策分析支持。      

Modeling urban heat islands in heterogeneous land surface and its correlation with impervious surface area by using night-time ASTER satellite data in highly urbanizing city,Delhi-India

The present study is an assessment and identification of urban heat island (UHI) in the environment of one of the fastest urbanizing city of India, Delhi Metropolis, employing satellite image of ASTER and Landsat 7 ETM+ in the thermal infrared region 3–14 μm. Temporal (2001 and 2005) ASTER datasets were used to analyze the spatial structure of the thermal urban environment subsequently urban heat island (UHI) in relation to the urban surface characteristics and land use/land cover (LULC). The study involves derivation of parameters governing the surface heat fluxes, constructing statistics of ASTER thermal infrared images along with validation through intensive in situ measurements. The average images reveal spatial and temporal variations of land surface temperature (LST) of night-time and distinct microclimatic patterns. Central Business District (CBD) of Delhi, (Connaught Place, a high density built up area), and commercial/industrial areas display heat islands condition with a temperature greater than 4 °C compared to the suburbs. The small increase in surface temperature at city level is mainly attributed to cumulative impact of human activities, changes in LULC pattern and vegetation density. In this study the methodology takes into account spatially-relative surface temperatures and impervious surface fraction value to measure surface UHI intensity between the urban land cover and rural surroundings. Both the spatial and temporal variation in surface temperature associated with impervious surface area (ISA) has been evaluated to assess the effect of urbanization on the local climate.     

Cramér-rao lower bound for pulsar rotation parameters estimation with X-ray pulsar observation data

In this paper, the Cramér-Rao Lower Bound (CRLB) for estimating the rotation parameters of pulsars using X-ray pulsar observation data is deduced, and the calculation equation is presented. In order to verify the correctness of the deduced equation, we use the X-ray pulsar observation data to estimate pulsar rotation parameters, and obtain the root mean square error (RMSE) of the estimated pulsar rotation parameters through conducting repeated experiments. The experimental results suggest that when the observation time increases, the RMSE gradually approaches the estimated CRLB, and that when the observation time is 2.4 × 10⁶ s, the error between the RMSE of pulsar frequency estimation and the CRLB stays at 10^?11 order of magnitude. This verifies that the CRLB expression deduced in this paper is the theoretical lower bound for estimating pulsar rotation parameters. The deduced CRLB in this paper helps determine the minimum variance estimator for pulsar rotation parameter estimation using X-ray pulsar data, providing a benchmark for the comparison between the minimum variance estimator and other estimators.     

An improved cyclic multi model-eXtreme gradient boosting (CMM-XGBoost) forecasting algorithm on the GNSS vertical time series

The study of GNSS vertical coordinate time series forecasting is helpful for monitoring the crustal plate movement, dam or bridge deformation monitoring, and global or regional coordinate system maintenance. The eXtreme Gradient Boosting (XGBoost) algorithm is a machine learning algorithm that can evaluate features, and it has a great potential and stability for long-span time series forecasting. This study proposes a multi-model combined forecasting method based on the XGBoost algorithm. The method constitutes a new time series as features through the fitting and forecasting results of the forecasting model. The XGBoost model is then used for forecasting. In addition, this method can obtain higher precision forecasting results through circulation. To verify the performance of the forecasting method, 1095 epochs of data in the Up coordinate of 16 GNSS stations are selected for the forecasting test. Compared with the CNN-LSTM model, the experimental results of our forecasting method show that the mean absolute error (MAE) values are reduced by 30.23 %～52.50 % and the root mean square error (RMSE) values are reduced by 31.92 %～54.33 %. The forecasting results have higher accuracy and are highly correlated to the original time series, which can better forecast the vertical movement of the GNSS stations. Therefore, the forecasting method can be applied to the up component of the GNSS coordinate time series.     

纳米酚醛气凝胶材料高温热物性参数辨识方法

临近空间高超声速飞行器大面积区域可能广泛采用纳米酚醛气凝胶(IPC)材料,获取高超声速气动加热作用下IPC材料的高温热物性参数,对于高超声速飞行器热防护系统的精细化设计具有重要的意义。考虑烧蚀效应的材料高温热物性参数辨识方法研究,基于Ablation Workshop烧蚀热响应标准算例对高温热物性参数辨识方法进行验证,计算结果表明：热物性参数辨识分析方法计算精度较高;通过带分层温度/烧蚀传感器的IPC材料电弧风洞试验,得到典型来流状态下不同厚度IPC材料内部的温度分布及热解厚度分布数据,通过辨识获得高温烧蚀条件下IPC材料热导率随温度的变化关系,IPC材料原始层热导率在温度低于800 K时随温度缓慢上升（热导率维持在0.1 W/(m·K)以下）,之后材料热解使得热导率发生突变,碳化层热导率在温度高于800 K时随着温度的上升急剧增大,到1 300 K左右时上升到0.17 W/(m·K)。     

基于DTW的HJ-1B卫星地表温度产品质量检验

HJ-1B卫星热红外数据应用广泛,但其地表温度反演产品的质量检验工作尚不完善。以黑河流域为研究区,利用普适性单通道算法得到HJ 1B地表温度,基于7个地面站点（下垫面为荒漠、沙漠、植被、农作物、雪地和湿地）同步观测资料和MODIS地表温度产品（MOD11A1）,引入动态时间规整方法（DTW）对站点处HJ 1B地表温度进行验证。试验结果表明：HJ 1B地表温度反演产品与地面观测值的偏差值在沙漠和荒漠站点为1K以内,均方根误差在05K左右;对于植被和农作物站点的偏差在2K以内,均方根误差为1~2K;基于DTW的验证对时序不匹配的数据评价结果与现有指标表现一致。HJ-1B地表温度反演产品与地面站点的相对偏差均低于其与MODIS地表温度反演产品的相对偏差。     

影响棱镜扫描式全景相机照相质量的主要因素

简述棱镜扫描式全景相机的基本原理和特点。根据大量的试验数据,详细地分析得出棱镜、象速度和胶片运动速度之间的同步误差以及温度变化是影响相机照相质量的主要因素的结论。文章还讨论了提高速度同步精度的技术途径。     

The effect of spatial truncation error on variance of gravity anomalies derived from inversion of satellite orbital and gradiometric data

The spatial truncation error (STE) is a significant systematic error in the integral inversion of satellite gradiometric and orbital data to gravity anomalies at sea level. In order to reduce the effect of STE, a larger area than the desired one is considered in the inversion process, but the anomalies located in its central part are selected as the final results. The STE influences the variance of the results as well because the residual vector, which is contaminated with STE, is used for its estimation. The situation is even more complicated in variance component estimation because of its iterative nature. In this paper, we present a strategy to reduce the effect of STE on the a posteriori   variance factor and the variance components for inversion of satellite orbital and gradiometric data to gravity anomalies at sea level. The idea is to define two windowing matrices for reducing this error from the estimated residuals and anomalies. Our simulation studies over Fennoscandia show that the differences between the 0.5°×0.5°$0.5°\times 0.5°$ gravity anomalies obtained from orbital data and an existing gravity model have standard deviation (STD) and root mean squared error (RMSE) of 10.9 and 12.1 mGal, respectively, and those obtained from gradiometric data have 7.9 and 10.1 in the same units. In the case that they are combined using windowed variance components the STD and RMSE become 6.1 and 8.4 mGal. Also, the mean value of the estimated RMSE after using the windowed variances is in agreement with the RMSE of the differences between the estimated anomalies and those obtained from the gravity model.     

Modelling of air temperature using remote sensing and artificial neural network in Turkey

The aim of this research was to forecast monthly mean air temperature based on remote sensing and artificial neural network (ANN) data by using twenty cities over Turkey. ANN contained an input layer, hidden layer and an output layer. While city, month, altitude, latitude, longitude, monthly mean land surface temperatures were chosen as inputs, and monthly mean air temperature was chosen as output for network. Levenberg–Marquardt (LM) learning algorithms and tansig, logsig and linear transfer functions were used in the network. The data of Turkish State Meteorological Service (TSMS) and Technological Research Council of Turkey–Bilten for the period from 1995 to 2004 were chosen as training when the data of 2005 year were being used as test. Result of research was evaluated according to statistical rules. The best linear correlation coefficient (R), and root mean squared error (RMSE) between the estimated and measured values for monthly mean air temperature with ANN and remote sensing method were found to be 0.991–1.254 K, respectively.     

随机系统的后验椭球条件滤波方法

针对在随机系统中存在着不确知的控制输入,且系统的状态以及随机干扰的分布和参数未知情况下,基于后验椭球逼近的原理,提出了随机系统的后验椭球条件滤波方法,以保证所求解的椭球序列是包含着系统状态估值域的最小椭球序列.该方法是当先验统计信息不足时,尤其当系统状态及随机干扰的分布形式未知时,是解决随机系统状态估计的有效途径之一.      

PCB红外热像辐射率校准的非线性滤波方法

在印制电路板(PCB)的红外辐射诊断中,为了获得PCB表面温度的真实分布,须对PCB红外热像进行辐射率校准.将PCB表面温度的测量波形看作真实温度分布与辐射率分布的混合波形,提出一种非线性滤波方法,利用二者连续性的差异将其区分开来.讨论了作为该方法关键步骤的红外热像突变点检测的小波变换实现问题.实验结果表明该方法能够准确地估计PCB表面温度的真实分布波形.该方法还可以用于PCB红外热像的热源辨识以及对未知材料的辐射率估计中.