Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (209)

Search Parameters:
Keywords = split window

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
26 pages, 658 KiB  
Article
Realistic Data Delays and Alternative Inactivity Definitions in Telecom Churn: Investigating Concept Drift Using a Sliding-Window Approach
by Andrej Bugajev, Rima Kriauzienė and Viktoras Chadyšas
Appl. Sci. 2025, 15(3), 1599; https://doi.org/10.3390/app15031599 - 5 Feb 2025
Viewed by 488
Abstract
Predicting customer churn is essential for telecommunications companies to maintain profitability. However, training models on historical models leads to performance degradation when they are applied to future conditions—a phenomenon known as concept drift. We employ a sliding-window approach that separates the training and [...] Read more.
Predicting customer churn is essential for telecommunications companies to maintain profitability. However, training models on historical models leads to performance degradation when they are applied to future conditions—a phenomenon known as concept drift. We employ a sliding-window approach that separates the training and testing time windows, creating a future-based “true test”. Using unique real data, we show that a CatBoost classifier model trained on older data can remain relevant when new, unseen intervals are used. A key innovation of our work is the use of 40-day “partial churn” labels; a model trained on these labels accurately predicts 90-day churn by simply adjusting the decision threshold. Out of the six modeled scenarios, in the main realistic scenario, CatBoost retained an accuracy above 0.798 and an F1 of near 0.704, reflecting its robustness even under real-world delays and potential drift. Overall, our findings emphasize that models do not necessarily “expire” with time; rather, their performance varies according to when they are tested. This research underscores the importance of a truly future-based evaluation (instead of artificial splits) and offers practical guidance for earlier churn detection when facing real-world data delays. Full article
Show Figures

Figure 1

Figure 1
<p>The number of articles by year obtained using the CatBoost method. This result was obtained from the Web Of Science platform with an advanced search query, “TS = (CatBoost) OR AK = (CatBoost)”.</p>
Full article ">Figure 2
<p>Simulation modeling. The red dotted line represents the moment in time that was simulated.</p>
Full article ">Figure 3
<p>Churn rate.</p>
Full article ">Figure 4
<p>F1 score.</p>
Full article ">Figure 5
<p>AUC (Area under curve).</p>
Full article ">Figure 6
<p>Accuracy.</p>
Full article ">Figure 7
<p>Recall.</p>
Full article ">Figure 8
<p>Precision.</p>
Full article ">Figure 9
<p>F1 score for different testing time windows vs. modeling moment time windows.</p>
Full article ">Figure 10
<p>F1 score differences using the L2 norm between scenarios.</p>
Full article ">Figure 11
<p>AUC differences using the L2 norm between scenarios.</p>
Full article ">Figure 12
<p>F1 performance with different model thresholds.</p>
Full article ">Figure 13
<p>F1 scores with the default and corrected thresholds.</p>
Full article ">Figure 14
<p>Confusion matrices with the default and corrected thresholds.</p>
Full article ">Figure A1
<p>Accuracy differences using L2 norm between scenarios.</p>
Full article ">Figure A2
<p>Accuracy differences using L2 norm between scenarios.</p>
Full article ">Figure A3
<p>Accuracy differences using L2 norm between scenarios.</p>
Full article ">Figure A4
<p>Definition of the hyperparameter space using the <tt>Hyperopt</tt> library.</p>
Full article ">
8 pages, 1424 KiB  
Proceeding Paper
A Convolutional Neural Network for Early Supraventricular Arrhythmia Identification
by Emilio J. Ochoa and Luis C. Revilla
Eng. Proc. 2025, 83(1), 8; https://doi.org/10.3390/engproc2025083008 - 8 Jan 2025
Viewed by 358
Abstract
Supraventricular arrhythmias (SVAs), including the often-asymptomatic supraventricular extrasystole (SVE), pose significant challenges in early detection and precise diagnosis. These challenges are of paramount importance, as recurrent SVEs may elevate the risk of developing severe SVAs, potentially resulting in cardiac weakening and subsequent heart [...] Read more.
Supraventricular arrhythmias (SVAs), including the often-asymptomatic supraventricular extrasystole (SVE), pose significant challenges in early detection and precise diagnosis. These challenges are of paramount importance, as recurrent SVEs may elevate the risk of developing severe SVAs, potentially resulting in cardiac weakening and subsequent heart failure. In the study conducted, an innovative approach was introduced that combined a convolutional neural network (CNN) architecture to enable the early identification and characterization of SVEs within electrocardiogram (ECG) signals. The analysis leveraged a dataset comprising 78 half-hour recordings from the highly regarded MIT-BIH Arrhythmia Database, which included annotation headers serving as labels for each recording. Signals were down-sampled by a factor of 2 and split into windows of 512 samples, with 12,288 observations for training. Following the methodology, classic signal preprocessing techniques (filtering and data normalization) were used. The proposed model was based on the UNET 1D model. A binary cross-entropy loss function, Adam optimizer, and a batch size of 128 were obtained after a hyperparameter tuning. As a training-validation methodology, a 50-fold cross-validation technique was used. The approach demonstrated a Dice coefficient of 79.01%, a precision of 80.96%, and a recall rate of 86.60% in detecting SVE events. These findings were corroborated through meticulous comparison with the annotations provided by the MIT-BIH database. The results underscore the immense potential of CNN and deep learning techniques in the early detection of supraventricular arrhythmias. This approach not only offers a valuable tool for healthcare professionals engaged in telemonitoring and early intervention strategies but also represents a significant contribution to the field of cardiac health monitoring. By facilitating efficient and precise identification of SVEs, our research sets the stage for improved patient outcomes and the prevention of severe SVAs, marking substantial advancements in this critical domain. Full article
Show Figures

Figure 1

Figure 1
<p>Architecture of the proposed convolutional neural network model.</p>
Full article ">Figure 2
<p>Project flowchart. It starts with the choice of databases, a preprocessing of the dataset, the preparation of the neural network, its evaluation, and, finally, its validation.</p>
Full article ">Figure 3
<p>Results of neural network segmentation for the detection of supraventricular extrasystole in an electrocardiogram signal. Interval detected.</p>
Full article ">Figure 4
<p>Results of neural network segmentation for the detection of supraventricular extrasystole in an electrocardiogram signal. ECG signal superposed.</p>
Full article ">Figure 5
<p>Probability level of recognition of an SVE (<b>left</b>) and its respective segmented signal in an ECG: red (training), green (validation), and brown (coincidence) (<b>right</b>).</p>
Full article ">
14 pages, 2052 KiB  
Article
Human Activity Recognition: A Comparative Study of Validation Methods and Impact of Feature Extraction in Wearable Sensors
by Saeed Ur Rehman, Anwar Ali, Adil Mehmood Khan and Cynthia Okpala
Algorithms 2024, 17(12), 556; https://doi.org/10.3390/a17120556 - 5 Dec 2024
Viewed by 728
Abstract
With the increasing availability of wearable devices for data collection, studies in human activity recognition have gained significant popularity. These studies report high accuracies on k-fold cross validation, which is not reflective of their generalization performance but is a result of the inappropriate [...] Read more.
With the increasing availability of wearable devices for data collection, studies in human activity recognition have gained significant popularity. These studies report high accuracies on k-fold cross validation, which is not reflective of their generalization performance but is a result of the inappropriate split of testing and training datasets, causing these models to evaluate the same subjects that they were trained on, making them subject-dependent. This study comparatively discusses this validation approach with a universal approach, Leave-One-Subject-Out (LOSO) cross-validation which is not subject-dependent and ensures that an entirely new subject is used for evaluation in each fold, validated on four different machine learning models trained on windowed data and select hand-crafted features. The random forest model, with the highest accuracy of 76% when evaluated on LOSO, achieved an accuracy of 89% on k-fold cross-validation, demonstrating data leakage. Additionally, this experiment underscores the significance of hand-crafted features by contrasting their accuracy with that of raw sensor models. The feature models demonstrate a remarkable 30% higher accuracy, underscoring the importance of feature engineering in enhancing the robustness and precision of HAR systems. Full article
(This article belongs to the Section Algorithms for Multidisciplinary Applications)
Show Figures

Figure 1

Figure 1
<p>The process of the activity recognition chain.</p>
Full article ">Figure 2
<p>Distribution of Activities in PAMAP2 dataset.</p>
Full article ">Figure 3
<p>Box plot of hand features.</p>
Full article ">Figure 4
<p>Variation in activities based on sensor readings.</p>
Full article ">Figure 5
<p>Plot of all models’ training vs. testing accuracies for k-fold and LOSO.</p>
Full article ">Figure 6
<p>Testing accuracy comparison of raw sensor models and feature-trained models validated with k-fold and LOSO on random forest and logistic regression classifiers.</p>
Full article ">Figure 7
<p>Training vs. testing accuracy in each subject in LOSO.</p>
Full article ">
15 pages, 1726 KiB  
Article
Forecasting Wind–Photovoltaic Energy Production and Income with Traditional and ML Techniques
by Giovanni Masala and Amelie Schischke
Econometrics 2024, 12(4), 34; https://doi.org/10.3390/econometrics12040034 - 12 Nov 2024
Viewed by 963
Abstract
Hybrid production plants harness diverse climatic sources for electricity generation, playing a crucial role in the transition to renewable energies. This study aims to forecast the profitability of a combined wind–photovoltaic energy system. Here, we develop a model that integrates predicted spot prices [...] Read more.
Hybrid production plants harness diverse climatic sources for electricity generation, playing a crucial role in the transition to renewable energies. This study aims to forecast the profitability of a combined wind–photovoltaic energy system. Here, we develop a model that integrates predicted spot prices and electricity output forecasts, incorporating relevant climatic variables to enhance accuracy. The jointly modeled climatic variables and the spot price constitute one of the innovative aspects of this work. Regarding practical application, we considered a hypothetical wind–photovoltaic plant located in Italy and used the relevant climate series to determine the quantity of energy produced. We forecast the quantity of energy as well as income through machine learning techniques and more traditional statistical and econometric models. We evaluate the results by splitting the dataset into estimation windows and test windows, and using a backtesting technique. In particular, we found evidence that ML regression techniques outperform results obtained with traditional econometric models. Regarding the models used to achieve this goal, the objective is not to propose original models but to verify the effectiveness of the most recent machine learning models for this important application, and to compare them with more classic linear regression techniques. Full article
Show Figures

Figure 1

Figure 1
<p>Power plant location.</p>
Full article ">Figure 2
<p>Power curve of the turbine.</p>
Full article ">Figure 3
<p>Backtesting and rolling windows.</p>
Full article ">Figure 4
<p>(<b>a</b>) Autocorrelation functions for the main variables. (<b>b</b>) Partial autocorrelation functions for the main variables.</p>
Full article ">
22 pages, 2149 KiB  
Article
Robust Biometric Verification Using Phonocardiogram Fingerprinting and a Multilayer-Perceptron-Based Classifier
by Roberta Avanzato, Francesco Beritelli and Salvatore Serrano
Electronics 2024, 13(22), 4377; https://doi.org/10.3390/electronics13224377 - 8 Nov 2024
Viewed by 807
Abstract
Recently, a new set of biometric traits, called medical biometrics, have been explored for human identity verification. This study introduces a novel framework for recognizing human identity through heart sound signals, commonly referred to as phonocardiograms (PCGs). The framework is built on extracting [...] Read more.
Recently, a new set of biometric traits, called medical biometrics, have been explored for human identity verification. This study introduces a novel framework for recognizing human identity through heart sound signals, commonly referred to as phonocardiograms (PCGs). The framework is built on extracting and suitably processing Mel-Frequency Cepstral Coefficients (MFCCs) from PCGs and on a classifier based on a Multilayer Perceptron (MLP) network. A large dataset containing heart sounds acquired from 206 people has been used to perform the experiments. The classifier was tuned to obtain the same false positive and false negative misclassification rates (equal error rate: EER = FPR = FNR) on chunks of audio lasting 2 s. This target has been reached, splitting the dataset into 70% and 30% training and testing non-overlapped subsets, respectively. A recurrence filter has been applied to also improve the performance of the system in the presence of noisy recordings. After the application of the filter on chunks of audio signal lasting from 2 to 22 s, the performance of the system has been evaluated in terms of recall, specificity, precision, negative predictive value, accuracy, and F1-score. All the performance metrics are higher than 97.86% with the recurrence filter applied on a window lasting 22 s and in different noise conditions. Full article
Show Figures

Figure 1

Figure 1
<p>Block diagram of biometric authentication modes: (<b>a</b>) identification mode; (<b>b</b>) verification mode.</p>
Full article ">Figure 2
<p>Examples of PCG recordings for 2 subjects: (<b>a</b>) 1st PCG recording of a female subject <math display="inline"><semantics> <msub> <mi>F</mi> <mrow> <mi>x</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </semantics></math>, (<b>b</b>) 2nd PCG recording of a female subject <math display="inline"><semantics> <msub> <mi>F</mi> <mrow> <mi>x</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </semantics></math>, (<b>c</b>) 1st PCG recording of a male subject <math display="inline"><semantics> <msub> <mi>M</mi> <mrow> <mi>y</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </semantics></math>, and (<b>d</b>) 2nd PCG recording of a male subject <math display="inline"><semantics> <msub> <mi>M</mi> <mrow> <mi>y</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </semantics></math>.</p>
Full article ">Figure 2 Cont.
<p>Examples of PCG recordings for 2 subjects: (<b>a</b>) 1st PCG recording of a female subject <math display="inline"><semantics> <msub> <mi>F</mi> <mrow> <mi>x</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </semantics></math>, (<b>b</b>) 2nd PCG recording of a female subject <math display="inline"><semantics> <msub> <mi>F</mi> <mrow> <mi>x</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </semantics></math>, (<b>c</b>) 1st PCG recording of a male subject <math display="inline"><semantics> <msub> <mi>M</mi> <mrow> <mi>y</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </semantics></math>, and (<b>d</b>) 2nd PCG recording of a male subject <math display="inline"><semantics> <msub> <mi>M</mi> <mrow> <mi>y</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </semantics></math>.</p>
Full article ">Figure 3
<p>Details (extracted at t = 20 s and lasting 2 s) of PCG recordings for 2 subjects: (<b>a</b>) 1st PCG recording of a female subject <math display="inline"><semantics> <msub> <mi>F</mi> <mrow> <mi>x</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </semantics></math>, (<b>b</b>) 2nd PCG recording of a female subject <math display="inline"><semantics> <msub> <mi>F</mi> <mrow> <mi>x</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </semantics></math>, (<b>c</b>) 1st PCG recording of a male subject <math display="inline"><semantics> <msub> <mi>M</mi> <mrow> <mi>y</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </semantics></math>, and (<b>d</b>) 2nd PCG recording of a male subject <math display="inline"><semantics> <msub> <mi>M</mi> <mrow> <mi>y</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </semantics></math>.</p>
Full article ">Figure 4
<p>Mel-frequency spectrograms and image representation of the related averaged vector <math display="inline"><semantics> <mover accent="true"> <mi>c</mi> <mo>→</mo> </mover> </semantics></math> evaluated for (<b>a</b>) 1st PCG chunk for female subject <math display="inline"><semantics> <msub> <mi>F</mi> <mi>x</mi> </msub> </semantics></math>, (<b>b</b>) 2nd PCG chunk for female subject <math display="inline"><semantics> <msub> <mi>F</mi> <mi>x</mi> </msub> </semantics></math>, (<b>c</b>) 1st PCG chunk for male subject <math display="inline"><semantics> <msub> <mi>M</mi> <mi>y</mi> </msub> </semantics></math>, and (<b>d</b>) 2nd PCG chunk for male subject <math display="inline"><semantics> <msub> <mi>M</mi> <mi>y</mi> </msub> </semantics></math>.</p>
Full article ">Figure 5
<p>Architecture of MLP binary classifier implemented to perform human identity verification by means of PCG features extracted from a segment of audio recording lasting 2 s; the output <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> </semantics></math> is true for identity verified, false for identity not verified.</p>
Full article ">Figure 6
<p>Architecture of the recurrence filter.</p>
Full article ">Figure 7
<p>Example of chunk affected by noises (2nd PCG recording of a female subject <math display="inline"><semantics> <msub> <mi>F</mi> <mrow> <mi>x</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </semantics></math>): (<b>a</b>) signal plus office noise, SNR = 15 dB (in red), (<b>b</b>) signal plus babble noise, SNR = 20 dB (in red), (<b>c</b>) signal plus babble noise, SNR = 30 dB (in red), (<b>d</b>) Mel spectrogram and features for chunk in (<b>a</b>), (<b>e</b>) Mel spectrogram and features for chunk in (<b>b</b>), and (<b>f</b>) Mel spectrogram and features for chunk in (<b>c</b>).</p>
Full article ">Figure 8
<p>Output of the loss function during the training phase for each epoch.</p>
Full article ">Figure 9
<p>ROC curve for HSCT-11 database without application of the “recurrence filter”.</p>
Full article ">Figure 10
<p>DET curve for HSCT-11 database without application of the “recurrence filter”. Red dashed line corresponds to the EER.</p>
Full article ">Figure 11
<p>Confusion matrices varying “recurrence filter” length obtained with the testing set portion of the system database: (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>L</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math>: no filtering, response time 2 s; (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>L</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>: response time 6 s; (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>L</mi> <mo>=</mo> <mn>4</mn> </mrow> </semantics></math>: response time 10 s; (<b>d</b>) <math display="inline"><semantics> <mrow> <mi>L</mi> <mo>=</mo> <mn>6</mn> </mrow> </semantics></math>: response time 14 s.</p>
Full article ">Figure 12
<p>Performance indexes varying “recurrence filter” length (response time) for the four datasets: (<b>a</b>) recall, (<b>b</b>) specificity, (<b>c</b>) precision, (<b>d</b>) negative predictive value, (<b>e</b>) accuracy, and (<b>f</b>) F1-score.</p>
Full article ">
23 pages, 32897 KiB  
Article
On the Suitability of Different Satellite Land Surface Temperature Products to Study Surface Urban Heat Islands
by Alexandra Hurduc, Sofia L. Ermida and Carlos C. DaCamara
Remote Sens. 2024, 16(20), 3765; https://doi.org/10.3390/rs16203765 - 10 Oct 2024
Cited by 2 | Viewed by 1451
Abstract
Remote sensing satellite data have been a crucial tool in understanding urban climates. The variety of sensors with different spatiotemporal characteristics and retrieval methodologies gave rise to a multitude of approaches when analyzing the surface urban heat island effect (SUHI). Although there are [...] Read more.
Remote sensing satellite data have been a crucial tool in understanding urban climates. The variety of sensors with different spatiotemporal characteristics and retrieval methodologies gave rise to a multitude of approaches when analyzing the surface urban heat island effect (SUHI). Although there are considerable advantages that arise from these different characteristics (spatiotemporal resolution, time of observation, etc.), it also means that there is a need for understanding the ability of sensors in capturing spatial and temporal SUHI patterns. For this, several land surface temperature products are compared for the cities of Madrid and Paris, retrieved from five sensors: the Spinning Enhanced Visible and InfraRed Imager onboard Meteosat Second Generation, the Advanced Very-High-Resolution Radiometer onboard Metop, the Moderate-resolution Imaging Spectroradiometer onboard both Aqua and Terra, and the Thermal Infrared Sensor onboard Landsat 8 and 9. These products span a wide range of LST algorithms, including split-window, single-channel, and temperature–emissivity separation methods. Results show that the diurnal amplitude of SUHI may not be well represented when considering daytime and nighttime polar orbiting platforms. Also, significant differences arise in SUHI intensity and spatial and temporal variability due to the different methods implemented for LST retrieval. Full article
(This article belongs to the Section AI Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Land cover resampled for the three projections of LST products (<b>a</b>–<b>f</b>) along with the percentage of urban pixels (<b>g</b>–<b>l</b>).</p>
Full article ">Figure 2
<p>Time of observation of each sensor: (<b>a</b>) for Madrid during daytime and time of minimum SUHI (SUHI<sub>min</sub>), (<b>b</b>) for Madrid during nighttime and time of maximum SUHI (SUHI<sub>max</sub>), (<b>c</b>) for Paris during daytime and SUHI<sub>max</sub>, (<b>d</b>) for Paris during nighttime and SUHI<sub>min</sub>. Colored bins are sampled every 15 min.</p>
Full article ">Figure 3
<p>Mean DJF (December, January, February) LST for all products considered. (<b>a1</b>–<b>a9</b>) The spatial pattern during daytime, and in the case of the MLST, the most frequent hour of the LST maximum and SUHI minimum are shown; (<b>b1</b>–<b>b9</b>) histograms of urban and rural LST shown in (<b>a1</b>–<b>a9</b>); (<b>c1</b>–<b>c7</b>) as in the first line but for nighttime and for the LST minimum and SUHI maximum; (<b>d1</b>–<b>d7</b>) as in (<b>b1</b>–<b>b9</b>) but for nighttime. Please note that color bars are different amongst the different products to allow a better visualization of patterns, but value ranges of the histograms are the same.</p>
Full article ">Figure 4
<p>As <a href="#remotesensing-16-03765-f003" class="html-fig">Figure 3</a> but for JJA (June, July, and August). (<b>a1</b>–<b>a9</b>) The spatial pattern during daytime, and in the case of the MLST, the most frequent hour of the LST maximum and SUHI minimum are shown; (<b>b1</b>–<b>b9</b>) histograms of urban and rural LST shown in (<b>a1</b>–<b>a9</b>); (<b>c1</b>–<b>c7</b>) as in the first line but for nighttime and for the LST minimum and SUHI maximum; (<b>d1</b>–<b>d7</b>) as in (<b>b1</b>–<b>b9</b>) but for nighttime. Please note that color bars are different amongst the different products to allow a better visualization of patterns, but value ranges of the histograms are the same.</p>
Full article ">Figure 5
<p>As <a href="#remotesensing-16-03765-f003" class="html-fig">Figure 3</a> but for Paris. (<b>a1</b>–<b>a9</b>) The spatial pattern during daytime, and in the case of the MLST, the most frequent hour of the LST and SUHI maximum; (<b>b1</b>–<b>b9</b>) histograms of urban and rural LST shown in (<b>a1</b>–<b>a9</b>); (<b>c1</b>–<b>c7</b>) as in the first line but for nighttime and for the LST and SUHI minimum; (<b>d1</b>–<b>d7</b>) as in (<b>b1</b>–<b>b9</b>) but for nighttime. Please note that color bars are different amongst the different products to allow a better visualization of patterns, but value ranges of the histograms are the same.</p>
Full article ">Figure 6
<p>As <a href="#remotesensing-16-03765-f003" class="html-fig">Figure 3</a> but for Paris and DJF; (<b>a1</b>–<b>a9</b>) The spatial pattern during daytime, and in the case of the MLST, the most frequent hour of the LST and SUHI maximum; (<b>b1</b>–<b>b9</b>) histograms of urban and rural LST shown in (<b>a1</b>–<b>a9</b>); (<b>c1</b>–<b>c7</b>) as in the first line but for nighttime and for the LST and SUHI minimum; (<b>d1</b>–<b>d7</b>) as in (<b>b1</b>–<b>b9</b>) but for nighttime, an extension of the histogram in (<b>d6</b>) is seen in (<b>d8</b>). Please note that color bars are different amongst the different products to allow a better visualization of patterns, but value ranges of the histograms are the same.</p>
Full article ">Figure 7
<p>Diurnal cycle of SUHI for Madrid: (<b>a</b>) DJF, (<b>b</b>) MAM, (<b>c</b>) JJA, (<b>d</b>) SON.</p>
Full article ">Figure 8
<p>As <a href="#remotesensing-16-03765-f007" class="html-fig">Figure 7</a> but for Paris.</p>
Full article ">Figure 9
<p>Correlation of monthly SUHI anomalies between all products considered: (<b>a</b>) daytime, (<b>b</b>) nighttime. Blank spaces correspond to pairs of products with no significant correlation (<span class="html-italic">p</span>-value &gt; 0.05).</p>
Full article ">Figure 10
<p>As <a href="#remotesensing-16-03765-f009" class="html-fig">Figure 9</a> but for Paris. (<b>a</b>) daytime, (<b>b</b>) nighttime. Blank spaces correspond to pairs of products with no significant correlation (<span class="html-italic">p</span>-value &gt; 0.05).</p>
Full article ">
16 pages, 13318 KiB  
Article
Investigation and Validation of Split-Window Algorithms for Estimating Land Surface Temperature from Landsat 9 TIRS-2 Data
by Qinghua Su, Xiangchen Meng and Lin Sun
Remote Sens. 2024, 16(19), 3633; https://doi.org/10.3390/rs16193633 - 29 Sep 2024
Viewed by 1252
Abstract
Land surface temperature (LST) is important in a variety of applications, such as urban thermal environment monitoring and water resource management. In this paper, eleven candidate split-window (SW) algorithms were adapted to Thermal Infrared Sensor-2 (TIRS-2) data of the Landsat 9 satellite for [...] Read more.
Land surface temperature (LST) is important in a variety of applications, such as urban thermal environment monitoring and water resource management. In this paper, eleven candidate split-window (SW) algorithms were adapted to Thermal Infrared Sensor-2 (TIRS-2) data of the Landsat 9 satellite for estimating the LST. The simulated dataset produced by extensive radiative transfer modeling and five global atmospheric profile databases was used to determine the SW algorithm coefficients. Ground measurements gathered at Surface Radiation Budget Network sites were used to confirm the efficiency of the SW algorithms after their performance was initially examined using the independent simulation dataset. Five atmospheric profile databases perform similarly in training accuracy under various subranges of total water vapor. The candidate SW algorithms demonstrate superior performance compared to the radiative transfer equation algorithm, exhibiting a reduction in overall bias and RMSE by 1.30 K and 1.0 K, respectively. It is expected to provide guidance for the generation of the Landsat 9 LST using the SW algorithms. Full article
(This article belongs to the Special Issue Surface Radiative Transfer: Modeling, Inversion, and Applications)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Performance of the eleven candidate SW algorithms under different total water vapor (TWV) subranges: (<b>a</b>) 0 &lt; TWV &lt; 1.5 g/cm<sup>2</sup>, (<b>b</b>) 1.5 &lt; TWV &lt; 3.0 g/cm<sup>2</sup>, (<b>c</b>) 3.0 &lt; TWV &lt; 4.5 g/cm<sup>2</sup>, (<b>d</b>) 4.5 &lt; TWV &lt; 10 g/cm<sup>2</sup>, and (<b>e</b>) 0 &lt; TWV &lt; 10 g/cm<sup>2</sup>.</p>
Full article ">Figure 1 Cont.
<p>Performance of the eleven candidate SW algorithms under different total water vapor (TWV) subranges: (<b>a</b>) 0 &lt; TWV &lt; 1.5 g/cm<sup>2</sup>, (<b>b</b>) 1.5 &lt; TWV &lt; 3.0 g/cm<sup>2</sup>, (<b>c</b>) 3.0 &lt; TWV &lt; 4.5 g/cm<sup>2</sup>, (<b>d</b>) 4.5 &lt; TWV &lt; 10 g/cm<sup>2</sup>, and (<b>e</b>) 0 &lt; TWV &lt; 10 g/cm<sup>2</sup>.</p>
Full article ">Figure 2
<p>Performance of the eleven candidate SW algorithms before and after adding random noise under different total water vapor (TWV) subranges: (<b>a</b>–<b>c</b>) 0 &lt; TWV &lt; 1.5 g/cm<sup>2</sup>, (<b>d</b>–<b>f</b>) 4.5 &lt; TWV &lt; 10 g/cm<sup>2</sup>. The random noise is added in the emissivity (<b>a</b>,<b>d</b>), brightness temperature (<b>b</b>,<b>e</b>), and TWV (<b>c</b>,<b>f</b>).</p>
Full article ">Figure 3
<p>The bias (<b>a</b>) and RMSE (<b>b</b>) between the retrieved Landsat 9 LST and the in situ LST for all matchups at the SURFRAD sites.</p>
Full article ">Figure 4
<p>The bias (<b>a</b>) and RMSE (<b>b</b>) between the retrieved Landsat 9 LST and the MODIS LST.</p>
Full article ">Figure 5
<p>Spatial distribution of the difference between the retrieved Landsat 9 LST and the MODIS LST. (<b>a</b>,<b>e</b>) LST retrieved from RTE algorithm, (<b>b</b>,<b>f</b>) LST retrieved from SW1 algorithm, (<b>c</b>,<b>g</b>) LST retrieved from SW3 algorithm, and (<b>d</b>,<b>h</b>) LST retrieved from SW5 algorithm.</p>
Full article ">Figure 5 Cont.
<p>Spatial distribution of the difference between the retrieved Landsat 9 LST and the MODIS LST. (<b>a</b>,<b>e</b>) LST retrieved from RTE algorithm, (<b>b</b>,<b>f</b>) LST retrieved from SW1 algorithm, (<b>c</b>,<b>g</b>) LST retrieved from SW3 algorithm, and (<b>d</b>,<b>h</b>) LST retrieved from SW5 algorithm.</p>
Full article ">Figure A1
<p>Scatterplots between the air temperature at the first layer and the total water vapor of the six atmospheric profile databases. (<b>a</b>) CLAR, (<b>b</b>) ERA5, (<b>c</b>) GAPRI (<b>d</b>) SeeBor, (<b>e</b>) TIGR and (<b>f</b>) WYO.</p>
Full article ">Figure A1 Cont.
<p>Scatterplots between the air temperature at the first layer and the total water vapor of the six atmospheric profile databases. (<b>a</b>) CLAR, (<b>b</b>) ERA5, (<b>c</b>) GAPRI (<b>d</b>) SeeBor, (<b>e</b>) TIGR and (<b>f</b>) WYO.</p>
Full article ">
17 pages, 6764 KiB  
Article
Fault Diagnosis Method for Vacuum Contactor Based on Time-Frequency Graph Optimization Technique and ShuffleNetV2
by Haiying Li, Qinyang Wang and Jiancheng Song
Sensors 2024, 24(19), 6274; https://doi.org/10.3390/s24196274 - 27 Sep 2024
Viewed by 672
Abstract
This paper presents a fault diagnosis method for a vacuum contactor using the generalized Stockwell transform (GST) of vibration signals. The objective is to solve the problem of low diagnostic performance efficiency caused by the inadequate feature extraction capability and the redundant pixels [...] Read more.
This paper presents a fault diagnosis method for a vacuum contactor using the generalized Stockwell transform (GST) of vibration signals. The objective is to solve the problem of low diagnostic performance efficiency caused by the inadequate feature extraction capability and the redundant pixels in the graph background. The proposed method is based on the time-frequency graph optimization technique and ShuffleNetV2 network. Firstly, vibration signals in different states are collected and converted into GST time-frequency graphs. Secondly, multi-resolution GST time-frequency graphs are generated to cover signal characteristics in all frequency bands by adjusting the GST Gaussian window width factor λ. The OTSU algorithm is then combined to crop the energy concentration area, and the size of these time-frequency graphs is optimized by 68.86%. Finally, considering the advantages of the channel split and channel shuffle methods, the ShuffleNetV2 network is adopted to improve the feature learning ability and identify fault categories. In this paper, the CKJ5-400/1140 vacuum contactor is taken as the test object. The fault recognition accuracy reaches 99.74%, and the single iteration time of model training is reduced by 19.42%. Full article
Show Figures

Figure 1

Figure 1
<p>Acceleration sensor mounting position diagram.</p>
Full article ">Figure 2
<p>Closing vibration signal waveforms in different states.</p>
Full article ">Figure 3
<p>GST time-frequency graphs of vibration signals in different states: (<b>a</b>) normal; (<b>b</b>) iron core rusting; (<b>c</b>) closing spring fatigue; (<b>d</b>) base screw loosening.</p>
Full article ">Figure 4
<p>Vibration signal GST time-frequency graph data augmentation.</p>
Full article ">Figure 5
<p>Cropping optimization of the GST time-frequency graph.</p>
Full article ">Figure 6
<p>ShuffleNetV2 network architecture.</p>
Full article ">Figure 7
<p>Two basic modules for ShuffleNetV2: (<b>a</b>) basic unit; (<b>b</b>) down-sampling unit.</p>
Full article ">Figure 8
<p>The framework of vacuum contactor fault diagnosis.</p>
Full article ">Figure 9
<p>Multi-resolution GST time-frequency graphs.</p>
Full article ">Figure 10
<p>The values of SSIM between time-frequency graph augmented data and original data.</p>
Full article ">Figure 11
<p>Optimization results of cropped GST time-frequency graphs: (<b>a</b>) normal; (<b>b</b>) iron core rusting; (<b>c</b>) closing spring fatigue; (<b>d</b>) base screw loosening.</p>
Full article ">Figure 12
<p>State recognition confusion matrices for GST time-frequency graphs: (<b>a</b>) without cropping optimization; (<b>b</b>) with cropping optimization.</p>
Full article ">Figure 13
<p>The accuracy of the test set for GST time-frequency graphs in four networks. AlexNet [<a href="#B28-sensors-24-06274" class="html-bibr">28</a>], ResNet50 [<a href="#B29-sensors-24-06274" class="html-bibr">29</a>], ResNeXt50 [<a href="#B30-sensors-24-06274" class="html-bibr">30</a>], ShuffleNetV2 [<a href="#B33-sensors-24-06274" class="html-bibr">33</a>].</p>
Full article ">
16 pages, 8514 KiB  
Article
Sinus Floor Augmentation with Synthetic Hydroxyapatite (NanoBone®) in Combination with Platelet-Rich Fibrin: A Case Series
by Luís Francisco, Manuel Francisco, Rosana Costa, Miguel Nunes Vasques, Marta Relvas, António Rajão, Luís Monteiro, Paulo Rompante, Fernando Guerra and Marco Infante da Câmara
Biomedicines 2024, 12(8), 1661; https://doi.org/10.3390/biomedicines12081661 - 25 Jul 2024
Viewed by 1489
Abstract
Several techniques have been described for maxillary sinus graft augmentation, including the lateral window technique and crestal approach with osteotomes or osseodensification. Platelet-rich fibrin has been used in maxillary sinus lift procedures due to its ability to accelerate soft and hard tissue healing. [...] Read more.
Several techniques have been described for maxillary sinus graft augmentation, including the lateral window technique and crestal approach with osteotomes or osseodensification. Platelet-rich fibrin has been used in maxillary sinus lift procedures due to its ability to accelerate soft and hard tissue healing. The aim of this study was to evaluate the potential of PRF in combination with the synthetic hydroxyapatite NanoBone® to enhance bone regeneration in sinus floor elevation with the lateral window technique. Out of the 50 individuals screened in a preoperative assessment visit from the CESPU—Famalicão clinical unit and intervened upon between January 2023 and December 2023, only 6 patients who met the study’s inclusion criteria consented to participate. In a split-mouth study, twelve sinus graft surgeries were carried out. Our observations reveal that for the test group (NanoBone®/PRF), there is a 27.5 ± 4.9% increase new vital bone, 23.0 ± 3.7% increase in inert bone particles, and 49.4 ± 2.8% increase in connective tissue. Meanwhile, for the control group (NanoBone®), there is a 19.5 ± 3.0% increase in new vital bone, 23.4 ± 5.7% increase in inert bone particles, and 57.0 ± 3.5% increase in connective tissue. The results strongly indicate that mixing liquid PRF with NanoBone® does not have a negative influence on the amount of viable bone formation, and it seems to slightly increase the amount of new bone formation and revascularization in sinus bone graft procedures with the lateral window technique compared to the single use of NanoBone®. Full article
(This article belongs to the Section Nanomedicine and Nanobiology)
Show Figures

Figure 1

Figure 1
<p>A visual representation of the surgery. (<b>a</b>) Bony window with Acteon Satelec<sup>®</sup> piezoelectric device; (<b>b</b>) Schneiderian membrane elevated; (<b>c</b>) Aggregation of NanoBone<sup>®</sup> with liquid fibrin; (<b>d</b>) Biomaterial insertion in the sinus cavity; (<b>e</b>) Straumann<sup>®</sup> Fex collagen membrane over the bony window.</p>
Full article ">Figure 2
<p>Harvesting of the bone specimen with a 2.5 mm diameter trephine bur. (<b>a</b>) Surgical site; (<b>b</b>) Trephine drill with collected bone; (<b>c</b>) Trephine drill with bone specimen in a 10% formaldehyde sterile vial.</p>
Full article ">Figure 3
<p>CONSORT flowchart.</p>
Full article ">Figure 4
<p>Newly formed bone, residual bone graft material, and connective tissue for both control (NanoBone<sup>®</sup> alone) and test groups (NanoBone<sup>®</sup>/liquid fibrin). Bars represent the mean and standard deviation of individuals results. (A) New vital bone, (B) particles, (C) new vital bone + particles, and (D) connective tissue. (*** <span class="html-italic">p</span> &lt; 0.0005).</p>
Full article ">Figure 5
<p>Histologic results in both groups: (<b>a</b>) test group; (<b>b</b>) control group. 10 × 0.5 magnification and toluidine blue staining.</p>
Full article ">Figure 6
<p>Histologic results in both groups: (<b>a</b>) test group; (<b>b</b>) control group. 10 × 0.5 magnification and toluidine blue staining.</p>
Full article ">Figure 7
<p>Histologic results in both groups: (<b>a</b>) test group; (<b>b</b>) control group. 10 × 0.5 magnification and toluidine blue staining.</p>
Full article ">Figure 8
<p>Histologic results in both groups: (<b>a</b>) test group; (<b>b</b>) control group. 10 × 0.5 magnification and toluidine blue staining.</p>
Full article ">Figure 9
<p>Histologic results in the test group. 10 × 0.5 magnification and toluidine blue staining.</p>
Full article ">Figure 10
<p>Histologic results in both groups: (<b>a</b>,<b>b</b>) test group; (<b>c</b>) control group. 10 × 0.5 magnification and toluidine blue staining.</p>
Full article ">Figure 11
<p>Histologic results in the control group. 10 × 0.5 magnification and toluidine blue staining.</p>
Full article ">
10 pages, 959 KiB  
Communication
One Algorithm to Rule Them All? Defining Best Strategy for Land Surface Temperature Retrieval from NOAA-AVHRR Afternoon Satellites
by Yves Julien, José A. Sobrino and Juan-Carlos Jiménez-Muñoz
Remote Sens. 2024, 16(15), 2720; https://doi.org/10.3390/rs16152720 - 25 Jul 2024
Viewed by 743
Abstract
The NOAA-AVHRR (National Oceanographic and Atmospheric Administration–Advanced Very High-Resolution Radiometer) archive includes data from 1981 onwards, which allow for estimating land surface temperature (LST), a key parameter for the study of global warming as well as surface characterization. However, algorithms for LST retrieval [...] Read more.
The NOAA-AVHRR (National Oceanographic and Atmospheric Administration–Advanced Very High-Resolution Radiometer) archive includes data from 1981 onwards, which allow for estimating land surface temperature (LST), a key parameter for the study of global warming as well as surface characterization. However, algorithms for LST retrieval were developed before the latest sensors and were based on more reduced atmospheric datasets. Here, we present 50 novel sets of coefficients for an LST retrieval algorithm from NOAA-AVHRR sensors, to which we added one historical methodology, which we validate against historical in situ as well as independent satellite data. This validation shows that the historical algorithm performs surprisingly well, with an in situ RMSE below 1.5 K and a quasi-null bias when compared with independent satellite data. A couple of the novel algorithms also perform within expectations (errors below 1.5 K), so any of these could be used for the complete processing of the AVHRR dataset. In our case, considering consistency with previous works, we opt for the use of the historical algorithm, now also tested for more recent periods. Full article
Show Figures

Figure 1

Figure 1
<p>Distribution of the satellite validation data through 1981–2022. The number of validation data points for each satellite are in parentheses.</p>
Full article ">
9 pages, 2884 KiB  
Comment
Comment on Yu et al. Land Surface Temperature Retrieval from Landsat 8 TIRS—Comparison between Radiative Transfer Equation-Based Method, Split Window Algorithm and Single Channel Method. Remote Sens. 2014, 6, 9829–9852
by Almustafa Abd Elkader Ayek and Bilel Zerouali
Remote Sens. 2024, 16(14), 2514; https://doi.org/10.3390/rs16142514 - 9 Jul 2024
Viewed by 1275
Abstract
Accurate land surface temperature (LST) retrieval from satellite data is pivotal in environmental monitoring and scientific research. This study addresses the impact of variability in the effective wavelengths used for LST retrieval from the Thermal Infrared Sensor (TIRS) data of Landsat 8. We [...] Read more.
Accurate land surface temperature (LST) retrieval from satellite data is pivotal in environmental monitoring and scientific research. This study addresses the impact of variability in the effective wavelengths used for LST retrieval from the Thermal Infrared Sensor (TIRS) data of Landsat 8. We conduct a detailed analysis comparing the effective wavelengths reported by Yu et al. (2014) and those derived from data provided by the USGS. Our analysis reveals significant variability in the effective wavelengths for bands 10 and 11 of Landsat 8. By applying Planck’s Law and utilizing the K1 and K2 coefficients available in the metadata of Landsat 8 products, we derive the effective wavelengths for bands 10 and 11. We also rederive the effective wavelength by integrating the spectral response function of the TIRS1 sensor. Our findings indicate that the effective wavelength for band 10 is 10.814 μm, aligning with the USGS data, while the effective wavelength for band 11 is 12.013 μm. We discuss the implications of these corrected effective wavelengths on the accuracy of LST retrieval algorithms, particularly the single channel algorithm (SC) and the radiative transfer equation (RT) employed by Yu et al. The importance of using precise effective wavelengths in satellite-based temperature retrieval is emphasized, to ensure the reliability and consistency of results. This analysis underscores the critical role of accurate spectral calibration parameters in remote sensing studies and provides valuable insights in the field of land surface temperature retrieval from Landsat 8 TIRS data. Full article
Show Figures

Figure 1

Figure 1
<p>TIRS relative spectral response.</p>
Full article ">Figure 2
<p>Spectral response function for Landsat 8 product (B10).</p>
Full article ">Figure 3
<p>Standard deviation for Landsat 8 product (B10).</p>
Full article ">Figure 4
<p>Spectral response function for Landsat 8 product (B11).</p>
Full article ">Figure 5
<p>Standard deviation for Landsat 8 product (B11).</p>
Full article ">Figure 6
<p>The relative spectral response function for band 10 using cubic interpolation.</p>
Full article ">Figure 7
<p>The relative spectral response function for band 11 using cubic interpolation.</p>
Full article ">
24 pages, 8649 KiB  
Article
Assessing the Impact of Land Use and Land Cover Changes on Surface Temperature Dynamics Using Google Earth Engine: A Case Study of Tlemcen Municipality, Northwestern Algeria (1989–2019)
by Imene Selka, Abderahemane Medjdoub Mokhtari, Kheira Anissa Tabet Aoul, Djamal Bengusmia, Malika KACEMI and Khadidja El-Bahdja Djebbar
ISPRS Int. J. Geo-Inf. 2024, 13(7), 237; https://doi.org/10.3390/ijgi13070237 - 2 Jul 2024
Cited by 3 | Viewed by 3270
Abstract
Changes in land use and land cover (LULC) have a significant impact on urban planning and environmental dynamics, especially in regions experiencing rapid urbanization. In this context, by leveraging the Google Earth Engine (GEE), this study evaluates the effects of land use and [...] Read more.
Changes in land use and land cover (LULC) have a significant impact on urban planning and environmental dynamics, especially in regions experiencing rapid urbanization. In this context, by leveraging the Google Earth Engine (GEE), this study evaluates the effects of land use and land cover modifications on surface temperature in a semi-arid zone of northwestern Algeria between 1989 and 2019. Through the analysis of Landsat images on GEE, indices such as normalized difference vegetation index (NDVI), normalized difference built-up index (NDBI), and normalized difference latent heat index (NDLI) were extracted, and the random forest and split window algorithms were used for supervised classification and surface temperature estimation. The multi-index approach combining the Normalized Difference Tillage Index (NDTI), NDBI, and NDVI resulted in kappa coefficients ranging from 0.96 to 0.98. The spatial and temporal analysis of surface temperature revealed an increase of 4 to 6 degrees across the four classes (urban, barren land, vegetation, and forest). The Google Earth Engine approach facilitated detailed spatial and temporal analysis, aiding in understanding surface temperature evolution at various scales. This ability to conduct large-scale and long-term analysis is essential for understanding trends and impacts of land use changes at regional and global levels. Full article
Show Figures

Figure 1

Figure 1
<p>Location map of Tlemcen, Algeria.</p>
Full article ">Figure 2
<p>Conceptual framework of the study.</p>
Full article ">Figure 3
<p>Classified LULC maps of Tlemcen: (<b>a</b>) 2019, (<b>b</b>) 2009, (<b>c</b>) 1999, and (<b>d</b>) 1989.</p>
Full article ">Figure 4
<p>Description of LULC statistics of Tlemcen from 1989 to 2019.</p>
Full article ">Figure 5
<p>LST distribution for Tlemcen in (<b>a</b>) 2019, (<b>b</b>) 2009, (<b>c</b>) 1999, and (<b>d</b>) 1989.</p>
Full article ">Figure 6
<p>Spatial normal difference vegetation index (NDVI) in Tlemcen. NDVI maps are shown for (<b>a</b>) 2019, (<b>b</b>) 2009, (<b>c</b>) 1999, and (<b>d</b>) 1989.</p>
Full article ">Figure 7
<p>Spatial NDBI in Tlemcen. NDBI maps are shown for (<b>a</b>) 2019, (<b>b</b>) 2009, (<b>c</b>) 1999, and (<b>d</b>) 1989.</p>
Full article ">Figure 8
<p>Spatial NDLI in Tlemcen. NDLI maps are shown for (<b>a</b>) 2019, (<b>b</b>) 2009, (<b>c</b>) 1999, and (<b>d</b>) 1989.</p>
Full article ">
22 pages, 6217 KiB  
Article
A Mathematical Model for Integrated Disaster Relief Operations in Early-Stage Flood Scenarios
by Nur Insani, Sona Taheri and Mali Abdollahian
Mathematics 2024, 12(13), 1978; https://doi.org/10.3390/math12131978 - 26 Jun 2024
Cited by 2 | Viewed by 2445
Abstract
When a flood strikes, the two most critical tasks are evacuation and relief distribution. It is essential to integrate these tasks, particularly before the floodwater reaches the vulnerable area, to minimize loss and damage. This paper presents a mathematical model of vehicle routing [...] Read more.
When a flood strikes, the two most critical tasks are evacuation and relief distribution. It is essential to integrate these tasks, particularly before the floodwater reaches the vulnerable area, to minimize loss and damage. This paper presents a mathematical model of vehicle routing problems to optimize an integrated disaster relief operation. The model addresses routing for both the evacuation and relief distribution tasks in the early stages of a flood, aiming to identify a minimal number of vehicles required with their corresponding routes to transport vulnerable individuals and simultaneously distribute emergency relief. The new model incorporates several features, including vehicle reuse, multi-trip and split delivery scenarios for evacuees and emergency relief items, uncertainty in evacuation demands, and closing time windows at evacuation points. Due to the complexity of vehicle routing problems, particularly in large-scale scenarios, the exact approach for obtaining optimal solutions is time-consuming. Therefore, we propose the use of a metaheuristic algorithm, specifically a modified genetic algorithm, to find an approximate solution for the proposed model. We apply the developed model and modified algorithm to various simulated flood scenarios and a real-life case study from Indonesia. The experimental results demonstrate that our approach requires fewer vehicles compared to standard models for similar scenarios. Moreover, while the exact approach fails to find optimal solutions within a reasonable timeframe for large-scale scenarios, our new approach provides near-optimal solutions in a much shorter time. In smaller simulated scenarios, the modified genetic algorithm obtains optimal or near-optimal solutions approximately 92.5% faster than the exact approach. Full article
(This article belongs to the Section E: Applied Mathematics)
Show Figures

Figure 1

Figure 1
<p>Relief disaster operations with evacuation and relief distribution tasks.</p>
Full article ">Figure 2
<p>The proposed modified genetic algorithm (MGA) utilising a modified best cost route as the crossover operator, a modified local search as the mutation operator and a constructive heuristic for generating an initial population and repairing it.</p>
Full article ">Figure 3
<p>Network traveling times.</p>
Full article ">Figure 4
<p>Evacuation and relief distribution plans of the integrated (<b>left</b>) and hierarchical (<b>right</b>) models.</p>
Full article ">Figure 5
<p>Bontoala district in Gowa regency, South Sulawesi, Indonesia.</p>
Full article ">Figure 6
<p>Distribution of temporary shelters in Bontoala, South Sulawesi, Indonesia.</p>
Full article ">
25 pages, 2095 KiB  
Article
Operational Angular Track Reconstruction in Space Surveillance Radars through an Adaptive Beamforming Approach
by Marco Felice Montaruli, Maria Alessandra De Luca, Mauro Massari, Germano Bianchi and Alessio Magro
Aerospace 2024, 11(6), 451; https://doi.org/10.3390/aerospace11060451 - 1 Jun 2024
Cited by 2 | Viewed by 1418
Abstract
In the last few years, many space surveillance initiatives have started to consider the problem represented by resident space object overpopulation. In particular, the European Space Surveillance and Tracking (EUSST) consortium is in charge of providing services like collision avoidance, fragmentation analysis, and [...] Read more.
In the last few years, many space surveillance initiatives have started to consider the problem represented by resident space object overpopulation. In particular, the European Space Surveillance and Tracking (EUSST) consortium is in charge of providing services like collision avoidance, fragmentation analysis, and re-entry, which rely on measurements obtained through ground-based sensors. BIRALES is an Italian survey radar belonging to the EUSST framework and is capable of providing measurements including Doppler shift, slant range, and angular profile. In recent years, the Music Approach for Track Estimate and Refinement (MATER) algorithm has been developed to retrieve angular tracks through an adaptive beamforming technique, guaranteeing the generation of more accurate and robust measurements with respect to the previous static beamforming approach. This work presents the design of a new data processing chain to be used by BIRALES to compute the angular track. The signal acquired by the BIRALES receiver array is down-converted and the receiver bandwidth is split into multiple channels, in order to maximize the signal-to-noise ratio of the measurements. Then, the signal passes through a detection block, where an isolation procedure creates, for each epoch, signal correlation matrices (CMs) related to the channels involved in the detection and then processes them to isolate the data stream related to a single detected source. Consequently, for each epoch and for each detected source, just the CM featuring the largest signal contribution is kept, allowing deriving the Doppler shift measurement from the channel illumination sequence. The MATER algorithm is applied to each CM stream, first estimating the signal directions of arrival, then grouping them in the observation time window, and eventually returning the target angular track. Ambiguous estimates may be present due to the configuration of the receiver array, which cause spatial aliasing phenomena. This problem can be addressed by either exploiting transit prediction (in the case of cataloged objects), or by applying tailored criteria (for uncatalogued objects). The performance of the new architecture was assessed in real operational scenarios, demonstrating the enhancement represented by the implementation of the channelization strategy, as well as the angular measurement accuracy returned by MATER, in both nominal and off-nominal scenarios. Full article
(This article belongs to the Special Issue Track Detection of Resident Space Objects)
Show Figures

Figure 1

Figure 1
<p>BIRALES receivers (red dots) installed along the antenna focal line of the eight cylindrical parabolic concentrators called <span class="html-italic">1N-section</span>. Each receiver contains 16 dipoles.</p>
Full article ">Figure 2
<p>MUSIC array responses, generated through simulated data based on BIRALES characteristics, for a single (<b>a</b>) and two sources (<b>b</b>) simultaneously detected.</p>
Full article ">Figure 3
<p>Track estimate (<b>a</b>) and track refinement (<b>b</b>) results in the cataloged case.</p>
Full article ">Figure 4
<p>Track estimate (<b>a</b>) and track refinement (<b>b</b>) results in the catalogued case.</p>
Full article ">Figure 5
<p>Track estimate (<b>a</b>) and track refinement (<b>b</b>) results in the cataloged case for a multiple source scenario.</p>
Full article ">Figure 6
<p>Track estimate result (<b>a</b>) and the DOA estimation selection based on <span class="html-italic">delta-k</span> technique-based criterion (<b>b</b>) in the uncatalogued case for the multiple sources scenario.</p>
Full article ">Figure 7
<p>Flow diagram of the BIRALES processing chain for the MATER algorithm.</p>
Full article ">Figure 8
<p>CM eigenvalues related to the acquisition on 18 April, by function of the observation time window and illuminated channels, before and after removing noise. Two detection trends can be noticed: the largest CM eigenvalue trend represents the satellite Sentinel-3B, for which the observation was scheduled, while the smaller one concerns the satellite Starlink-4114, which was also detected during the acquisition.</p>
Full article ">Figure 9
<p>Time evolution of the detected signal, in terms of illuminated channels and Doppler shift. Two detection trends can be noticed: the longer one (Detection A) is related to the detection of the satellite Sentinel-3B, while the shorter one is related (Detection B) to the satellite Starlink-4114.</p>
Full article ">Figure 10
<p>Clustered CM eigenvalues related to the acquisition on 18 April. Two detection trends can be noticed: the larger eigenvalue one (Detection A) is related to the detection of the satellite Sentinel-3B, while the smaller eigenvalue is related (Detection B) to the satellite Starlink-4114.</p>
Full article ">Figure 11
<p>Doppler shift time profile returned by the pipeline. Detection A refers to the satellite Sentinel-3B, while Detection B refers to the satellite Starlink-4114.</p>
Full article ">Figure 12
<p>DOA estimations and angular tracks returned by the MATER algorithm. Detection A refers to the satellite Sentinel-3B, while Detection B refers to the satellite Starlink-4114. The angular track (continuous black line) is longer than the path covered by the estimated signal DOAs, as the former is related to all detection epochs, while the latter is related to a reduced time window, guaranteeing a more accurate DOA estimation in the track estimate phase of the MATER algorithm.</p>
Full article ">Figure 13
<p>SARAL satellite observation: eigenvalue time trend. By processing the data with the original architecture (<b>a</b>), that is without the channelization strategy, the two largest CM eigenvalue trends could be identified, related to a Cassiopea-A radio-astronomical source (the larger) and to the target satellite (the smaller). With the new architecture pipeline (<b>b</b>), just the latter trend remained.</p>
Full article ">Figure 14
<p>SARAL satellite observation: array response at a generic acquisition epoch. By processing the data with the original architecture (<b>a</b>), that is, without the channelization strategy, both the DOA estimation of the Cassiopea-A radio-astronomical source and that of the target satellite can be noticed, together with the related ambiguities. With the new architecture pipeline (<b>b</b>), at the same generic epoch, only the DOA estimation of the target satellite remains, together with the ambiguous solutions.</p>
Full article ">Figure 15
<p>SARAL satellite observation: MATER angular track. By processing the data with the original architecture (<b>a</b>), that is, without the channelization strategy, the DOA estimations of both the Cassiopea-A radio-astronomical source (circle clusters) and that of the target satellite (which constitutes angular sequences) can be noticed. With the new architecture pipeline (<b>b</b>), only the DOA estimations related to the target satellite remain. It can be also noticed that the DOA estimations in the new architecture are denser than in the original one, because of the shorter integration time, and the angular track returned by MATER is longer.</p>
Full article ">Figure 16
<p>Observation of the Aeolus satellite on 27 July 2023, during its assisted re-entry: DOA estimations and angular track returned by the MATER algorithm in the uncatalogued case. The ambiguity was solved thanks to the geometrical filters.</p>
Full article ">Figure 17
<p>Observation of the Aeolus satellite on 27 July 2023, during its assisted re-entry: transit prediction delay with respect to the MATER estimation represented through the angular coordinate evolution throughout the entire acquisition (<b>a</b>) and using the BIRALES array response at a generic acquisition epoch (<b>b</b>).</p>
Full article ">Figure 18
<p>Calibrator data set: cumulative distribution functions describing the orbits of the observed objects at the detection epoch. The orbits are represented in terms of the semi-major axis, eccentricity, inclination, and right ascension of the ascending node.</p>
Full article ">Figure 19
<p>Calibrator data set: cumulative distribution functions of the angular RMSE, for both cataloged and uncatalogued cases, for the two angular coordinates separately. Overall, the cataloged and uncatalogued cases showed similar results, and the angular track turned out to be generally more accurate along the <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>γ</mi> <mn>1</mn> </msub> </mrow> </semantics></math> direction than along the <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>γ</mi> <mn>2</mn> </msub> </mrow> </semantics></math> one, because of the length of the angular path traveled during the integration time.</p>
Full article ">Figure 20
<p>Calibrators data set: angular track and DOA estimations returned by the MATER algorithm for the observation of the satellite COSMOS 2517 (norad ID 41579) conducted on 7 April 2022. An incorrect angular track was returned by the exploited ambiguity solving criterion, as the assumption that the source spent a remarkable portion of the transit in the central region of the receiver FoV was not matched.</p>
Full article ">
30 pages, 1637 KiB  
Article
Enhancing Monthly Streamflow Prediction Using Meteorological Factors and Machine Learning Models in the Upper Colorado River Basin
by Saichand Thota, Ayman Nassar, Soukaina Filali Boubrahimi, Shah Muhammad Hamdi and Pouya Hosseinzadeh
Hydrology 2024, 11(5), 66; https://doi.org/10.3390/hydrology11050066 - 1 May 2024
Cited by 1 | Viewed by 4076
Abstract
Streamflow prediction is crucial for planning future developments and safety measures along river basins, especially in the face of changing climate patterns. In this study, we utilized monthly streamflow data from the United States Bureau of Reclamation and meteorological data (snow water equivalent, [...] Read more.
Streamflow prediction is crucial for planning future developments and safety measures along river basins, especially in the face of changing climate patterns. In this study, we utilized monthly streamflow data from the United States Bureau of Reclamation and meteorological data (snow water equivalent, temperature, and precipitation) from the various weather monitoring stations of the Snow Telemetry Network within the Upper Colorado River Basin to forecast monthly streamflow at Lees Ferry, a specific location along the Colorado River in the basin. Four machine learning models—Random Forest Regression, Long short-term memory, Gated Recurrent Unit, and Seasonal AutoRegresive Integrated Moving Average—were trained using 30 years of monthly data (1991–2020), split into 80% for training (1991–2014) and 20% for testing (2015–2020). Initially, only historical streamflow data were used for predictions, followed by including meteorological factors to assess their impact on streamflow. Subsequently, sequence analysis was conducted to explore various input-output sequence window combinations. We then evaluated the influence of each factor on streamflow by testing all possible combinations to identify the optimal feature combination for prediction. Our results indicate that the Random Forest Regression model consistently outperformed others, especially after integrating all meteorological factors with historical streamflow data. The best performance was achieved with a 24-month look-back period to predict 12 months of streamflow, yielding a Root Mean Square Error of 2.25 and R-squared (R2) of 0.80. Finally, to assess model generalizability, we tested the best model at other locations—Greenwood Springs (Colorado River), Maybell (Yampa River), and Archuleta (San Juan) in the basin. Full article
Show Figures

Figure 1

Figure 1
<p>Visual Representation of Analyzed Factors (Temperature, Snow Water Equivalent, Precipitation, and Streamflow) in Streamflow Prediction.</p>
Full article ">Figure 2
<p>Colorado River Basin by United States Geological Survey (USGS) (<a href="https://www.usgs.gov/media/images/colorado-river-basin-map" target="_blank">https://www.usgs.gov/media/images/colorado-river-basin-map</a> (accessed on 10 February 2024)).</p>
Full article ">Figure 3
<p>Spatial Distribution of SNOTEL Weather Stations and USGS Gauges within the Upper Colorado River Basin.</p>
Full article ">Figure 4
<p>Time Series Plots of SNOTEL Data - Snow Water Equivalent (SWE) in inches, Precipitation Accumulation (Prcp_Acc) in inches, and Average Temperature (Avg_Temp) in degrees Fahrenheit.</p>
Full article ">Figure 5
<p>Time Series Plot of Streamflow Data from the US Bureau of Reclamation.</p>
Full article ">Figure 6
<p>Temporal Variation of Streamflow in Relation to Meteorological Variables.</p>
Full article ">Figure 7
<p>Correlation Analysis between Streamflow and Meteorological Variables.</p>
Full article ">Figure 8
<p>A Simplified Structure of RFR.</p>
Full article ">Figure 9
<p>The Structure of LSTM Memory Unit.</p>
Full article ">Figure 10
<p>Comparison of RMSE Distribution between Univariate and Multivariate Models.</p>
Full article ">Figure 11
<p>Comparison of Median <math display="inline"><semantics> <msup> <mi>R</mi> <mn>2</mn> </msup> </semantics></math> between Univariate and Multivariate Time-Series Models.</p>
Full article ">Figure 12
<p>MAE, RMSE, SMAPE, MAPE, and <math display="inline"><semantics> <msup> <mi>R</mi> <mn>2</mn> </msup> </semantics></math> Results of Univariate Time-Series ML Models.</p>
Full article ">Figure 13
<p>RMSE Distribution of Univariate Time-Series Models.</p>
Full article ">Figure 14
<p>RMSE Distribution of Multivariate Time-Series Models.</p>
Full article ">Figure 15
<p>RMSE, SMAPE, MAPE, and <math display="inline"><semantics> <msup> <mi>R</mi> <mn>2</mn> </msup> </semantics></math> Results of Multivariate Time-Series ML Models.</p>
Full article ">Figure 16
<p>Comparison of Predicted and Observed USBR Streamflow Using Multivariate Time-Series Models (RFR, LSTM, GRU) with Input-Output Combinations of 12-12 and 24-12, for the Period from May 2015 to April 2016 over the Test Set.</p>
Full article ">Figure 17
<p>Comparison of Predicted vs. Observed Streamflow from Multivariate Time-Series Models (RFR, LSTM, GRU) over the Test Set (May 2015–April 2016).</p>
Full article ">Figure 18
<p>Heatmap Visualization of Evaluation Metrics (MAPE, SMAPE, <math display="inline"><semantics> <msup> <mi>R</mi> <mn>2</mn> </msup> </semantics></math>, MAE, and RMSE) for Univariant RFR Model.</p>
Full article ">Figure 19
<p>Predictions of Univariate Time-Series RFR Models with Different Look-back Span Windows against Ground Truth (May 2018–April 2019).</p>
Full article ">Figure 20
<p>Heatmap Visualization of Evaluation Metrics (MAPE, SMAPE, <math display="inline"><semantics> <msup> <mi>R</mi> <mn>2</mn> </msup> </semantics></math>, MAE, and RMSE) for Multivariate RFR Model.</p>
Full article ">Figure 21
<p>Predictions of Multivariate Time-Series RFR Models with Different Look-back Span Windows against Ground Truth (May 2018–April 2019).</p>
Full article ">Figure 22
<p>Top Features’ Combinations Count.</p>
Full article ">Figure 23
<p>Feature Importance of Meteorological Factors.</p>
Full article ">Figure 24
<p>Boxplot of <math display="inline"><semantics> <msup> <mi>R</mi> <mn>2</mn> </msup> </semantics></math> for RFR Model Across Different Streamflows.</p>
Full article ">
Back to TopTop