Next Article in Journal
An Efficient Framework for Securing the Smart City Communication Networks
Next Article in Special Issue
Single-Element and MIMO Circularly Polarized Microstrip Antennas with Negligible Back Radiation for 5G Mid-Band Handsets
Previous Article in Journal
Clinical–Functional Evaluation and Test–Retest Reliability of the G-WALK Sensor in Subjects with Bimalleolar Ankle Fractures 6 Months after Surgery
Previous Article in Special Issue
What Can 5G Do for Public Safety? Structural Health Monitoring and Earthquake Early Warning Scenarios
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Contribution of Singular Spectral Analysis to Forecasting and Anomalies Detection of Indoors Air Quality

by
Felipe Espinosa
1,
Ana B. Bartolomé
1,
Pablo Villoria Hernández
2 and
M. C. Rodriguez-Sanchez
2,*
1
Electronics Department, University of Alcala, E-28801 Alcalá de Henares, Spain
2
Electronics Technology Department, University Rey Juan Carlos, E-28933 Móstoles, Spain
*
Author to whom correspondence should be addressed.
Submission received: 24 March 2022 / Revised: 11 April 2022 / Accepted: 11 April 2022 / Published: 15 April 2022
(This article belongs to the Special Issue Sensors in 2022)

Abstract

:
The high impact of air quality on environmental and human health justifies the increasing research activity regarding its measurement, modelling, forecasting and anomaly detection. Raw data offered by sensors usually makes the mentioned time series disciplines difficult. This is why the application of techniques to improve time series processing is a challenge. In this work, Singular Spectral Analysis (SSA) is applied to air quality analysis from real recorded data as part of the Help Responder research project. Authors evaluate the benefits of working with SSA processed data instead of raw data for modelling and estimation of the resulting time series. However, what is more relevant is the proposal to detect indoor air quality anomalies based on the analysis of the time derivative SSA signal when the time derivative of the noisy original data is useless. A dual methodology, evaluating level and dynamics of the SSA signal variation, contributes to identifying risk situations derived from air quality degradation.

1. Introduction

Air quality, indoors and outdoors, has a direct impact on the environment and, fundamentally, on the health and comfort of citizens [1]. Improving air quality is therefore a challenge for both public administrations and private entities. The tendency of populations to concentrate in large cities has increased the interest in this topic, hence the continuous references in publications related to Smart Cities [2,3,4,5,6].
To assess air pollution, different pollutants are considered: carbon monoxide and dioxide, ozone, nitrogen dioxide and sulphur dioxide. Although these are not detectable with the naked eye, there are commercial devices that include the necessary algorithm to report standardised air quality indices (AQ Index) [7,8,9]. Besides, there are visible air pollutants in the form of solid particles emitted to the atmosphere, particulate matter PM10 and PM2.5 [10].
The starting point for monitoring air quality is to have sensor devices [11], generally forming wireless sensor network nodes supported by IoT technology [5,12,13]. Subsequent data processing is essential, both to have models that anticipate the evolution of this variable, and to alert researchers to anomalous variations detected in its recording.
The record provided by air quality sensors usually incorporates significant noise components. Part of this noise is due to the nature of the variable itself, conditioned by climatological phenomena (temperature, humidity, atmospheric pressure), and part to the characteristic noise of the sensor itself. Therefore, the processing of the original signal is essential to facilitate its subsequent analysis, modelling and prediction. In the case of univariate time series, techniques such as autocorrelation and periodogram are traditionally used to reduce the noise components inherent to the recording of time series; however, these techniques have significant limitations in defining the number of regressors in the modelling of the series [14,15].
The contribution of this work focuses on the application of Singular Spectral Analysis (SSA) to the processing of data from air quality sensors. From the principal components of the signal, a double objective is proposed: to obtain a model that allows predicting its evolution with a defined horizon and to analyse the temporal evolution changes in order to detect risk situations due to the spontaneous or forced introduction of pollutants. Figure 1 shows a schematic presentation of this contribution.
The authors’ contribution to smart node Yongzhi sensing can be extended to a spatially distributed wireless sensing network as described in [13,16] in the Smart City context.

2. SSA Fundamentals

Singular Spectrum Analysis (SSA) is a powerful technique in time series analysis with a diverse application field. The first publications date from 1986, but theoretical and methodology aspects can also be found in several recent references [16,17,18]; this is why a mathematical review is included in this chapter.
The aim of SSA is to make a decomposition of the original series into the sum of a small number of independent and interpretable components and a structureless noise. This way, slowly varying trends, seasonality components or cycles with small and large periods can be isolated without needing to know the model of the original signal.
We consider a real-valued nonzero time series of length N, X N = ( x 1 ,     x N ) . The main purpose of SSA is to decompose the original series into a sum of components (identified as dominant/principal components) and noise. This is followed by a reconstruction the original series.
Given a window length L   ( 1 < L < N ) , we construct the L-lagged vectors:
X i = ( x i ,     x i + L 1 ) T ,   i = 1 ,   2 ,   ,   K = N L + 1
and compose these vectors into the X matrix:
X = ( x i + j 1 ) i , j = 1 L , K = [ X 1 : : X K ]
It is a Hankel (square and symmetric) matrix with size L x K and is often called `trajectory matrix’. The singular-value decomposition of the matrix C = X X T yields a collection of L eigenvalues λ i and eigenvectors v i :
C = E Λ E T
where Λ is a diagonal matrix containing the singular values ( σ i = λ i , ) of C and E = [ v 1 . v L ] is the orthogonal matrix formed with the eigenvectors v i .
Operating with E and C matrices, the diagonal (eigenvalues) Λ matrix is obtained as shown in Figure 2. A particular combination of a certain number l < L ( l or principal components) of these eigenvectors, determines an l -dimensional subspace in L . The L-dimensional data X = { X 1 ,   ,   X K } is then projected onto this l -dimensional subspace E ¯ to obtain the X ¯ submatrix (sized K x l ) that allows to reconstruct the Hankel matrix X ˜ , as symbolically described in Figure 3. The X ˜ is considered as an approximation to X , where the parameters L and l (Principal Components—PC) are chosen by the designer for a specific application.

3. HelpResponder Project—Case Study: Air Quality Evaluation

There are several alternatives and devices for monitoring air quality. In [19], an MQ135 sensor is used for its wide range for gases detection (NH3, NOx, CO2, etc.). In [8], a low-cost device to measure Indoor Air Quality (IAQ) is used: BME680. The gas sensor BME680 is a digital device capable to measure gas (in a unit of Ohm), humidity (RH), pressure (hPa) and temperature (°C) that complies with ISO16000-29 standard “Test methods for Volatile Organic Compounds (VOC) detectors”. The operating range of BME680 is −40/+85 °C, 0/100% RH and 300/1100 hPa, which suits general conditions. The IAQ is calculated by a smart algorithm using these measurement parameters, developed by Bosch Software Environmental Cluster (BSEC) based on the German Federal Environmental Agency’s guidelines.
The interpretation of computed IAQ is shown in Table 1, where lower IAQ infers better air quality [20]. The combination of both gas and chemoresistive sensors allows the calculation of an Air Quality Index based on the concentration of volatile organic compounds (VOCs), CO and NOx.
The authors belong to the research team involved in the HelpResponder project. The global goal of the project is to provide Emergency Professionals with innovative technology making their emergency activity easier and guaranteeing their own safety as well as that of the people to be evacuated. The final platform deals with the integration of electronics, communication and services that increase the performance of portable equipment and its interaction with Management and Supervision Staff. Among the tasks to be carried out is the detection of level and degradation of indoor air quality.
In this context, Figure 4 shows an actual record obtained at the Alcorcón Fire Department (Madrid). The test was carried out between 27 October and 3 November 2020; recording IAQ values provided by a BME680 sensor with a 1 min sampling. Over some time intervals, smoke was introduced to force a degradation of the air quality. As can be seen, between minute 3000 and 8800 there are intervals with very bad IAQ levels (see Table 1). The records in blue show results obtained under normal condition and the data in red those obtained under forced disturbance condition.

4. SSA Contribution to Air Quality Estimation

The next measurement estimation or the variable prediction with a given time horizon requires a model of time series under study. In the case of non-linear time series such as air quality measurements, intelligent different solutions have been recently published: based on shallow and deep-learning predictors [21], neuro-fuzzy network approach [22] or alternatives based on support vector machines [23]. However, to highlight the SSA contribution to this goal, a first approach combining linear and non-linear modelling techniques is discussed, as shown in Figure 5, from a set of j regressors, and adjusting the corresponding offset.
In this work we have used a solution in which the linear component responds to an autoregressive AR model; for the non-linear component we have compared three solutions: a sigmoidal function, another based on the wavelet transform and another based on the treepartition technique. Matlab has specific functions for each (sigmoidnet, wavenet and treepartition) [24,25,26]. The main configuration parameters are the number of regressors (order) and the number of non-linear functions (units), which can be chosen by the user or set automatically.
To evaluate the SSA contribution to air quality modelling, the following methodology is proposed: (a) modelling and prediction of the time series from the original IAQ data provided by the sensor; (b) modelling and prediction of the time series after SSA data processing.
For the first phase of the study, the original 3000 data recorded under normal (smokeless) condition were selected. The first 2000 data were used for identification and the last 1000 for validation [27]. A minimum modelling fit (FIT) of 80% was set as an acceptance criterion, and automatic identification was carried out for the three alternatives previously mentioned. The validation results for an order of 50 are shown in Table 2.
After the obtained results, the treepartition model of 50 order and 31 units was chosen, as it presented the lowest MSE (Mean Squared Error) with the shortest computation load. The relative value t = 1 corresponds to 1.447 s using an HP Pavilion x360 Convertible PC, with an Intel Core i7-7500U processor, 2 cores, 12 GB DDR4 RAM and a NVIDIO GeForce 940 MX graphics card.
Considering this model, an estimation was made with a half-hour horizon (30 samples) and the model results were compared with the recorded results (last 1100 data in Figure 4). In Figure 6, the grey trace represents part of the original signal data_zn_val(y1) and the blue trace is the estimation provided by the model with the indicated time horizon. The percentage of fit is slightly higher than 70% mainly due to the limitations of the model to fit the original signal with significant noise content.
One of the objectives of time series modelling is to predict the variable under study beyond the last recorded sample. For this purpose, the Matlab forecast function [25,26] was used in this work. Figure 7 shows the IAQ prediction results for half an hour after the last recorded data.
In the following, the SSA contribution to air quality modelling and prediction is evaluated. Taking as reference the works of [17,28], a time window LN/12 (250 samples) was chosen, where N is the total of the recorded data for the analysis in the time interval without forced disturbance. The covariance matrix was calculated and the eigenvalues shown in Figure 8 were obtained. The signal reconstruction obtained with L = 250 and the dominant eigenvalues (RCs = 4) is shown in Figure 9; the percentage of fit is 80.27%. As indicated in Table 3, increasing the number of principal components (eigenvalues) does not significantly improve the reconstruction percentage of the time series.
Using the same nonlinear model structure (treepartition, 50 order and 31 units), but with the SSA processed data, the model validation results are: Fit = 99.12%, MSE = 0.0035 and computation time 0.71 s. Figure 10 shows the graphical modelling results. Compared to the modelling results in Table 2, obtained using the raw data (original signal), the MSE was divided by 433, the Fit was improved by 16.88% and the computation time was reduced by 51%. The same model is applicable to the prediction of IAQ evolution. Figure 11 shows the one-day prediction (1440 samples, with a period of one minute). Compared to the prediction in Figure 6, the one obtained from the SSA processed signal shows a smoother dynamic.

5. Air Quality Anomalies Detection

Air quality degradation can have different causes. These can be natural, such as the increase of COx, NOx, Sox or O3 levels, generally with slow dynamics, or forced by emergency situations such as fire, criminal actions with chemical components, etc. In any case, these anomalies result in very bad IAQ levels, above 60 (see Table 1) and fast dynamics. This is the case of the forced disturbances shown in red in Figure 12.
Knowing the critical consequences of IAQ degradation, the authors propose a dual approach to characterise areas of anomalous behaviour: level and dynamic detection, as shown in Figure 13. A level alarm will be generated when the IAQ value exceeds the threshold I A Q t h = 60% (very bad IAQ according to Table 1). In the case of dynamic detection, it is unfeasible to draw conclusions from the time derivative of the original signal (raw data), as shown in Figure 14(left). However, the time derivative of the SSA processed signal (Figure 14(right)) provides relevant information about anomaly detection.
Focused on SSA IAQ time derivative, Figure 15(top) compares the two evaluated scenarios: without (magenta line) and with (blue line) forced disturbance. As can be seen, in the experimental recording performed in the HelpResponder Project, the smoke introduced in the scenario under test generated a time derivative (DZ derivative) with a higher level and duration. In order to better discriminate between the Normal Zone (NZ) and the Disturbance Zone (DZ), the integral was applied to both signals, for | d I A Q d t | > 0.1 % / min values, obtaining the representation shown in Figure 15(down). The interpretation of the area generated from the time derivatives, in both cases, facilitates the detection of air quality anomalies through a threshold for the area value.

6. Conclusions

Singular Spectral Analysis is a well-known noisy signals processing technique, that can also be used for time series. It facilitates the study and analysis of a reconstructed signal from its dominant components. In the context of the HelpResponder Project, the authors have highlighted the contribution of SSA to the processing of the air quality signal recorded with commercial sensors at 1 min sampling.
SSA contributes to the modelling and estimation of non-linear time series such as the one addressed in this paper. After a comparative study of Matlab functions for non-linear time series modelling, treepartition was chosen. The estimation result of the IAQ signal processed with SSA, with respect to the original one (raw data), provides a higher fit level, lower MSE error and less computational time. In addition, it facilitates the prediction of future values, such as those shown for air quality with a prediction horizon of one day.
However, where the contribution of SSA stands out most is in the analysis of the time derivative of the series, admitting that the IAQ sensor does not saturate. The derivative of the original signal, with an important high-frequency component, is of little use, as shown in this work. However, after applying SSA, the time derivative of the series provides information of interest, both in terms of reached level (%/s) and time duration. In short, the application of SSA to air quality monitoring offers a different alternative of signal treatment for the rapid detection of degradation, whether due to natural or forced causes (emergencies).

Author Contributions

Data curation, F.E. and P.V.H.; Formal analysis, F.E.; Funding acquisition, M.C.R.-S.; Investigation, F.E. and M.C.R.-S.; Methodology, F.E. and M.C.R.-S.; Project administration, M.C.R.-S.; Software, A.B.B.; Supervision, F.E. and M.C.R.-S.; Validation, F.E. and P.V.H.; Writing—original draft, F.E., A.B.B. and P.V.H.; Writing—review & editing, F.E. and M.C.R.-S.. All authors have read and agreed to the published version of the manuscript.

Funding

This research and APC are funded by [Comunidad de Madrid] grant number [M2184].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization (WHO). Monitoring Health for the SDGs: Sustainable Development Goals; World Health Organization (WHO): Geneva, Switzerland, 2019. [Google Scholar]
  2. Rodríguez-Molano, J.I.; Obando-Bobadilla, L.M.; Ruiz-Nieto, M.P. Of cities traditional to smart cities. In Proceedings of the 2018 13th Iberian Conference on Information Systems and Technologies (CISTI), Caceres, Spain, 13–16 June 2018; IEEE: Piscataway, NJ, USA, 2018. [Google Scholar]
  3. Sanchez, L.; Galache, J.A.; Gutierrez, V.; Hernandez, J.M.; Bernat, J.; Gluhak, A.; Garcia, T. SmartSantander: The meeting point between Future Internet research and experimentation and the smart cities. In Proceedings of the 2011 Future Network Mobile Summit, Warsaw, Poland, 15–17 June 2011; pp. 1–8. Available online: http://0-ieeexplore-ieee-org.brum.beds.ac.uk/abstracl/documenl/6095264 (accessed on 10 April 2022).
  4. Santos, C.; Jiménez, J.A.; Espinosa, F. Effect of event-based sensing on IoT node power efficiency. Case study: Air quality monitoring in smart cities. IEEE Access 2019, 7, 132577–132586. [Google Scholar] [CrossRef]
  5. Kalajdjieski, J.; Stojkoska, B.R.; Trivodaliev, K. IoT Based Framework for Air Pollution Monitoring in Smart Cities. In Proceedings of the 28th Telecommunications forum TELFOR 2020, Belgrade, Serbia, 24–25 November 2020. [Google Scholar]
  6. González, E.; Casanova-Chafer, J.; Romero, A.; Vilanova, X.; Mitrovics, J.; Llobet, E. LoRa Sensor Network Development for Air Quality Monitoring or Detecting Gas Leakage Events. Sensors 2020, 20, 6225. [Google Scholar] [CrossRef] [PubMed]
  7. Popoola, O.A.M.; Carruthers, D.; Lad, C.; Bright, V.B.; Mead, M.I.; Stettler, M.E.J.; Saell, J.R.; Jones, R.L. Use of networks of low cost air quality sensors to quantify air quality in urban settings. Atmos. Environ. 2018, 194, 58–70. [Google Scholar] [CrossRef]
  8. Lasomsri, P.; Yanbuaban, P.; Kerdpoca, O.L.; Ouypornkochagorn, T. A Development of Low-Cost Devices for Monitoring Indoor Air Quality in a Large-Scale Hospital. In Proceedings of the 2018 15th International Conference on Electrical Engineering/Electronics, Computer, Telecommunications and Information Technology, Chiang Rai, Thailand, 18–21 July 2018. [Google Scholar] [CrossRef]
  9. Polichetti, T.; Miglietta, M.L.; Alfano, B.; Massera, E.; de Vito, S.; di Francia, G.; Faucon, A.; Saoutie, E.; Boisseau, S.; Marchand, N.; et al. A Networked Wearable Device for Chemical Multisensing”. Lecture Notes in Electrical Engineering; Springer International Publishing: Cham, Germany, 2019; pp. 17–24. [Google Scholar]
  10. Okigbo, C.A.; Seeam, A.; Guness, S.P.; Bellekens, X.; Bekaroo, G.; Ramsurrun, V. Low Cost Air Quality Monitoring: Comparing the Energy Consumption of an Arduino against a Raspberry Pi Based System; 2020 International Conference on Intelligent and Innovative Computing Applications (ICONIC); ACM: New York, NY, USA, 2020; 8p. [Google Scholar] [CrossRef]
  11. Kim, J.Y.; Chu, C.H.; Shin, S.M. ISSAQ: An integrated sensing systems for real-time indoor air quality monitoring. IEEE Sens. J. 2014, 14, 4230–4244. [Google Scholar] [CrossRef]
  12. Sendra, S.; Garcia-Navas, J.L.; Romero-Diaz, P.; Lloret, J. Collaborative LoRa-Based Sensor Network for Pollution Monitoring in Smart Cities. In Proceedings of the 2019 Fourth International Conference on Fog and Mobile Edge Computing (FMEC), Rome, Italy, 10–13 June 2019; IEEE: Piscataway, NJ, USA, 2019. [Google Scholar]
  13. Khalifeh, A.; Darabkh, K.A.; Khasawneh, A.M.; Alqaisieh, I.; Salameh, M.; AlAbdala, A.; Alrubaye, S.; Alassaf, A. Wireless Sensor Networks for Smart Cities: Network Design, Implementation and Performance Evaluation. Electronics 2021, 10, 218. [Google Scholar] [CrossRef]
  14. Lütkepohl, H. New Introduction to Multiple Time Series Analysis; Springer: Berlin/Heidelberg, Germany, 2005; ISBN 3-540-40172-5. [Google Scholar]
  15. Brockwell, P.J.; Davis, R.A. Intoduction to Time Series and Forecasting; Springer: Berlin/Heidelberg, Germany, 2016; ISBN 978-3-319-29852-8. [Google Scholar] [CrossRef]
  16. Yongzhi, J. Application of fault detection using distributed sensors in smart cities. Phys. Commun. 2021, 46, 101182. [Google Scholar] [CrossRef]
  17. Golyandina, N.E.; Zhigljavsky, A. Singular Spectrum Analysis for Time Series; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar] [CrossRef]
  18. Hassani, H. Singular Spectrum Analysis: Methodology and Comparison. J. Data Sci. 2007, 5, 239–257. [Google Scholar] [CrossRef]
  19. Ali-Kazmi, S.N.; Ulasyar, A.; Nadeem-Khan, M.F. IoT based Energy Efficient Smart Street Lighting Technique with Air Quality Monitoring. In Proceedings of the 2020 14th International Conference on Open Source Systems and Technologies (ICOSST), Lahore, Pakistan, 16–17 December 2020. [Google Scholar] [CrossRef]
  20. EPA. Indoor Air Quality|EPA’s Report on the Environment (ROE)|US EPA. 2021. Available online: https://www.epa.gov/report-environment/indoor-air-quality (accessed on 25 November 2021).
  21. Liu, H.; Yan, G.; Duan, Z.; Chen, C. Intelligent modeling strategies for forecasting air quality time series: A review. Appl. Soft Comput. 2021, 102, 106957. [Google Scholar] [CrossRef]
  22. Lin, Y.-C.; Lee, S.-J.; Ouyang, C.-S.; Wu, C.-H. Air quality prediction by neuro fuzzy modeling approach. Appl. Soft Comput. 2020, 86, 105898. [Google Scholar] [CrossRef]
  23. Espinosa, R.; Palma, J.; Jiménez, F.; Kaminska, J.; Sciavicco, G.; Lucena-Sánchez, E. A time series forecasting based multi-criteria methodology for air quality prediction. Appl. Soft Comput. 2021, 113, 107850. [Google Scholar] [CrossRef]
  24. Ljung, L. System Identification Toolbox. User’s Guide; Matlab & Simulink, MathWorks: Natick, MA, USA, 2020. [Google Scholar]
  25. Ljung, L. System Identification Toolbox. Getting Started Guide; Matlab & Simulink, MathWorks: Natick, MA, USA, 2020. [Google Scholar]
  26. Ljung, L. System Identification Toolbox. Reference; Matlab & Simulink, MathWorks: Natick, MA, USA, 2020. [Google Scholar]
  27. Ljung, L. System Identification. Theory for the User, 2nd ed.; Prentice Hall: Hoboken, NJ, USA, 1999; ISBN 1-13-656695-2. [Google Scholar]
  28. Marques, C.A.F.; Ferreira, J.A.; Rocha, A.; Castanheira, J.M.; Melo-Gonçalves, P.; Vaz, N.; Dias, J.M. Singular spectrum analysis and forecasting of hydrological time series. Phys. Chem. Earth Parts A/B/C 2006, 31, 1172–1179. [Google Scholar] [CrossRef]
Figure 1. Contribution scheme of SSA to air quality forecasting and anomalous detection.
Figure 1. Contribution scheme of SSA to air quality forecasting and anomalous detection.
Sensors 22 03054 g001
Figure 2. Steps involved in the SSA decomposition stage.
Figure 2. Steps involved in the SSA decomposition stage.
Sensors 22 03054 g002
Figure 3. SSA reconstruction stage from the principal components.
Figure 3. SSA reconstruction stage from the principal components.
Sensors 22 03054 g003
Figure 4. Recorded IAQ data under normal condition and incorporating smoke int the indoor area under test.
Figure 4. Recorded IAQ data under normal condition and incorporating smoke int the indoor area under test.
Sensors 22 03054 g004
Figure 5. Scheme of non-linear time series modelling.
Figure 5. Scheme of non-linear time series modelling.
Sensors 22 03054 g005
Figure 6. Comparison between original signal (do not included at the identification stage) and estimation with the treepartition model proposed by authors.
Figure 6. Comparison between original signal (do not included at the identification stage) and estimation with the treepartition model proposed by authors.
Sensors 22 03054 g006
Figure 7. IAQ prediction with a half-hour horizon, from the treepartition time series model.
Figure 7. IAQ prediction with a half-hour horizon, from the treepartition time series model.
Sensors 22 03054 g007
Figure 8. Eigenvalue representation applying SSA to the IAQ time series, being the time window L = 250.
Figure 8. Eigenvalue representation applying SSA to the IAQ time series, being the time window L = 250.
Sensors 22 03054 g008
Figure 9. IAQ time series (without disturbances) reconstruction with L = 250 and four dominant eigenvalues.
Figure 9. IAQ time series (without disturbances) reconstruction with L = 250 and four dominant eigenvalues.
Sensors 22 03054 g009
Figure 10. Estimation result using the treepartition model (50 order and 31 units) from the SSA reconstructed signal instead of the original one.
Figure 10. Estimation result using the treepartition model (50 order and 31 units) from the SSA reconstructed signal instead of the original one.
Sensors 22 03054 g010
Figure 11. One day IAQ prediction (1440 min) based on the SSA processing of the original measurements.
Figure 11. One day IAQ prediction (1440 min) based on the SSA processing of the original measurements.
Sensors 22 03054 g011
Figure 12. Original (left) and SSA processed (right) IAQ signal. The blue line belongs to registered data in normal conditions, while the red line is related to forced disturbances in the scenario under test. According to Table 1, 60% is the threshold to consider very bad air quality.
Figure 12. Original (left) and SSA processed (right) IAQ signal. The blue line belongs to registered data in normal conditions, while the red line is related to forced disturbances in the scenario under test. According to Table 1, 60% is the threshold to consider very bad air quality.
Sensors 22 03054 g012
Figure 13. Dual methodology for air quality anomalies detection.
Figure 13. Dual methodology for air quality anomalies detection.
Sensors 22 03054 g013
Figure 14. Original (left) and SSA processed temporal derivative (dIAQ/dt) signal (right).
Figure 14. Original (left) and SSA processed temporal derivative (dIAQ/dt) signal (right).
Sensors 22 03054 g014
Figure 15. Top: Time derivative of the processed IAQ signal in both cases: without disturbance (NZ) and with disturbance (DZ). Down: Area calculated from the derivative signals with higher level to 0.1%/min; the area value is also referred.
Figure 15. Top: Time derivative of the processed IAQ signal in both cases: without disturbance (NZ) and with disturbance (DZ). Down: Area calculated from the derivative signals with higher level to 0.1%/min; the area value is also referred.
Sensors 22 03054 g015
Table 1. Interpretation of IAQ based on the Bosch algorithm.
Table 1. Interpretation of IAQ based on the Bosch algorithm.
IAQ Index (%) (*)Air Quality
0–10Good
10–20Average
20–30Little Bad
30–40Bad
40–60Worse
60–100Very Bad
(*) The values are converted to the percentage of the maximum value of 500.
Table 2. Results of the validation stage for different nonlinear time series model (order = 50), working with the IAQ original data.
Table 2. Results of the validation stage for different nonlinear time series model (order = 50), working with the IAQ original data.
Matlab FunctionFIT (%)MSERelative Time
Wavenet81.671.6142.684
Sigmoidnet81.771.59715.923
Treepartition82.241.5161
Table 3. Effect of the eigenvalues number on the IAQ time series reconstruction (fit percentage) applying SSA.
Table 3. Effect of the eigenvalues number on the IAQ time series reconstruction (fit percentage) applying SSA.
SSA: Reconstruction Fit of the IAQ Time Series, L = 250
Eigenvalues4681012
Fit (%)80.2781.3381.7782.0982.42
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Espinosa, F.; Bartolomé, A.B.; Hernández, P.V.; Rodriguez-Sanchez, M.C. Contribution of Singular Spectral Analysis to Forecasting and Anomalies Detection of Indoors Air Quality. Sensors 2022, 22, 3054. https://0-doi-org.brum.beds.ac.uk/10.3390/s22083054

AMA Style

Espinosa F, Bartolomé AB, Hernández PV, Rodriguez-Sanchez MC. Contribution of Singular Spectral Analysis to Forecasting and Anomalies Detection of Indoors Air Quality. Sensors. 2022; 22(8):3054. https://0-doi-org.brum.beds.ac.uk/10.3390/s22083054

Chicago/Turabian Style

Espinosa, Felipe, Ana B. Bartolomé, Pablo Villoria Hernández, and M. C. Rodriguez-Sanchez. 2022. "Contribution of Singular Spectral Analysis to Forecasting and Anomalies Detection of Indoors Air Quality" Sensors 22, no. 8: 3054. https://0-doi-org.brum.beds.ac.uk/10.3390/s22083054

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop