Characterization of Temperature Gradients According to Height in a Baroque Church by Means of Wireless Sensors
Abstract
:1. Introduction
1.1. Microclimatic Monitoring for the Preservation of Cultural Heritage
1.2. Microclimatic Studies with Sensors Located at Different Heights
2. Materials and Methods
2.1. Description of the Monitoring System
2.2. Experiment for the Calibration of Temperature Sensors
2.3. Installation of Wireless Nodes
- : nodes I and J were located at the upper part of the retable decorating the presbytery.
- : nodes A, F, P, and Q were placed at the upper position, close to the ceiling vaults.
- : nodes L and M were also located at the retable.
- : nodes G, H and O were placed near to the main altar.
- : it corresponds to node U, which was located near the main entrance.
- : nodes K and N were also installed at the retable.
- : node D was located close to the altarpiece of Saint Joseph.
- : nodes B and T were positioned near to the main entrance.
- : nodes C, E, R, and S were located, as indicated in Figure 4, at the lowest level.
2.4. Data Pretreatment
2.5. Statistical Methods
2.5.1. Identification of Stages in the Time Series
2.5.2. Estimation of the Vertical Gradient of Temperature for Each Month
2.5.3. Calculation of Classification Variables
- Method 1: Using Time Series FunctionsThis method consists of computing features from the observed time series , in some cases, and from the time series after applying the logarithm transformation and regular differencing to . The goal of using this transformation and differencing was to stabilize the variance and remove the trend of the series in order to extract information about the seasonal component. Features were calculated by means of values of sample Auto Correlation Function (ACF), sample Partial Auto Correlation Function (PACF), periodogram, Moving Range (MR) [71,72], as well as features defined using quantiles [50]. Each variable was computed for each month and sensor. These correspond to estimates of the following parameters:
- (a)
- mean.ts: Mean of recorded in the month. This parameter allows to compare the level of the different time series.
- (b)
- sd.ts: Standard deviation of , which provides information about the variability of the recorded values.
- (c)
- range.ts: Range of (i.e., by subtracting the minimum to the maximum). It reflects the amplitude of the time series of T and gives information about the dispersion.
- (d)
- mean.mr: Mean of MR values with order 24 of . MR computes the moving range for all sequences of 24 consecutive observations.
- (e)
- median.mr: Median of MR values with order 24 of . This parameter and the previous one are helpful for capturing the daily variability of the different time series of T.
- (f)
- mean.acf: Mean of the first 72 lags () of sample ACF applied to time series. Each value of ACF for at lag l () is the correlation coefficient between the observations that are lagged for a time gap l. It is given by , i.e., Pearson’s correlation coefficient between the time series and the lagged values (i.e., the time gap which is considered). The value 72 was used because sample ACF values computed for were comprehended within the limits of a 95% confidence interval in the correlogram. This parameter provides information about the dynamic structure of the time series.
- (g)
- median.acf: Median of the first 72 lags of sample ACF applied to . As in the previous case, this parameter can be useful for comparing the dynamic structure of the time series.
- (h)
- sd.acf: Standard deviation of the first 72 lags of sample ACF of .
- (i)
- pacf: First 4 lags () of sample PACF applied to . A value of PACF at lag l measures the autocorrelation between the observation and , which is not accounted for by lags 1 to . The first four values of PACF are usually the most important ones for capturing the most significant autocorrelation information. These four values were computed trying to differentiate the dynamic structure of the different time series.
- (j)
- maximum.I: Maximum value from the periodogram (I), which is employed for identifying the dominant periods or frequencies of time series of T. This parameter is helpful for recognizing the dominant cyclical behavior in a series.
- (k)
- range.I: Range of values of the periodogram. This parameter can be useful to compare the impact of the dominant cyclical pattern in the different series.
- (l)
- maximum.slps: Maximum increase of T in one hour found in the month (i.e., ). This parameter allows the comparison of the maximum changes of for two consecutive hours, and it is intended to capture the information of abnormal peaks or sudden increases due to occasional events.
- (m)
- median.abs.sd: Median of absolute values of the deviation between the values of and the median of . It is given by . This parameter is somewhat related to the variance (i.e., average of the squared deviations with respect to the mean) and, hence, it is another measurement of data dispersion.
- (n)
- t.p.r.m20: It is computed as , being the percentile a of values in the month. Thus, it is the ratio of percentiles (60th–40th) over (95th–5th) of . The numerator is the range of variability corresponding to 20% of the central part of the original time series. The denominator is basically the range of the original time series after removing the lowest 5% and highest 5%. An equivalent interpretation corresponds to the parameters t.p.r.m35, t.p.r.m50, and t.p.r.m80 described next.
- (o)
- t.p.r.m35: It is computed as , which is the ratio of percentiles (67.5th–32.5th) over (95th–5th) of .
- (p)
- t.p.r.m50: Ratio of percentiles (75th–25th) over (95th–5th) of .
- (q)
- t.p.r.m80: Ratio of percentiles (90th–10th) over (95th–5th) of .
- (r)
- p.d.f.p: Ratio of percentiles (95th–5th) over the median of . This parameter divides the amplitude (range) of the time series, after removing the lowest 5% and highest 5% of observations, by the median of .
This list comprises a set of 21 variables that were computed for each one of the seven months, which implies 147 variables in total. They were arranged in a matrix denoted as comprised of 21 rows (one per node) and 147 columns (one per variable). - Method 2: Additive Seasonal Holt-Winters Method (SH-W)This approach calculates features from time series of T, by using the Holt-Winters method (SH-W) [73], which is an extension of the Holt’s method [74]. It captures the level, trend, and seasonality of the different time series and is comprised of the forecast equation and three smoothing equations (i.e., one for the level , one for the trend or slope , and one for the seasonal component ) with corresponding smoothing parameters , , and [75]. According to the additive SH-W, the forecast equation for a time series of T with period length p is given by Equation (2) (in this study, p is 24), where k is the integer part of , and is the forecast at step [75].Slope, level and seasonal components at step i are estimated by using the three smoothing equations (i.e., for , , and ), respectively. If the algorithm converges, a, b and to are the estimations for the level, trend or slope and seasonal components. This algorithm was run by using the function HoltWinters of the stats package [76] of R software.The flow diagram for the additive SH-W method is displayed in Figure 6. In this diagram, all the steps are repeated with each observation of time series , . However, in step (1), the initial values of level (), trend, () and seasonal coefficients () are only used once to start up the algorithm. The initial conditions are estimated through a simple decomposition in trend and seasonal component by using moving averages. After initialization, steps from (2) to (4) perform the forecast task internally, these values were updated and stored for the next step [76]. In step (2), the estimation of slope requires knowledge of the level at steps i, , and so on until , as well as slope at steps , and so on until . In step (3), as in step (2), the equation is solved recursively. Estimation of the level requires knowledge of the level, slope, seasonal components at different steps starting at (for ), (for ), (for ), and finishing when the values are , , and . It also requires values of T at steps i, and so on until , where is just the oldest data point in the training data set (i.e., a set of observations starting from until the current observation ). Note that the weighting coefficients , and need to be computed for running steps (2), (3) and (4). Such coefficients are calculated by minimizing the squared one-step prediction error [76]. Now that the level, trend and seasonal component at time step i have been estimated, the forecast at step with can be estimated by using the three values of components together.According to this method, the level, trend, and seasonal components are updated over a historical period. For example, when the method is applied per month, the components are updated every hour over each month. If the algorithm converges, a, b and to are the estimated values for the level, trend and seasonal components at the last instant of time in the month.The level at a time t corresponds to a weighted average between the seasonally adjusted temperature and the level forecast, based on the level and slope at the previous instance of time . This component gives an estimate of the local mean (i.e., mean per hour in this study). Regarding the slope component, it expresses the linear increment of the level, over an hour. Finally, the seasonality component estimates the deviation from the local mean, due to seasonality.The features calculated per sensor are the following:
- (a)
- a: Estimated value for the level for each month of the time series.
- (b)
- b: Estimated value for the trend (slope) for each month.
- (c)
- s1,s2,…,s24: Estimated values for the seasonal components for each month.
- (d)
- sse: Sum of squared estimate of errors per month.
- (e)
- maximum.I: Maximum value of the periodogram computed with the residuals of SH-W for each month.
- (f)
- mean.acf: Mean of sample ACF of residuals at lags 1 to 72 per month.
- (g)
- median.acf: Median of sample ACF of residuals at lags 1 to 72 for each month.
- (h)
- range.acf: Range of sample ACF of residuals at lags 1 to 72 per month.
- (i)
- Dn: Statistic of the Kolgomorov–Smirnov () normality test [77] of the residuals derived from SH-W, per month of the time series. The normality test was employed to compare the empirical distribution function of the residuals with the cumulative distribution function of the normal model.
- (j)
- Wn: Statistic of the Shapiro–Wilk test () [78] of the residuals per month. This test was used to detect deviations from normality, because of either kurtosis or skewness, or both. The Dn and Wn statistics were also used as classification variables, because they provide information about deviation from normality for the residuals derived from the SH-W method.
- (k)
- fcast: 24 forecasts of T (i.e., , ) for a unique additive SH-W model that was fitted using the complete time series without splitting it in different months.
Features calculated from (a) to (j) imply a set of 33 variables computed for each month. By including the 24 forecasts as explained in (k), the total number of variables was , which were organized as a matrix denoted as , comprised of 21 rows (one per sensor) and 255 columns (one per variable).
2.5.4. sPLS
2.5.5. Linear Discriminant Analysis (LDA)
3. Results and Discussion
3.1. Vertical Gradients of Temperature
3.2. Ventilation of the Church of Saint Thomas and Saint Philip Neri
3.3. Application of sPLS to Identify Key Features Correlated with Height
3.4. Discrimination of Sensors in Three Categories by Means of LDA
4. Conclusions
- This research reports a microclimatic study in the church of Saint Thomas and Saint Philip Neri in Valencia for the first time, which is of relevant interest because inappropriate conditions of temperature can affect the valuable artworks. The results suggest that temperature gradients in this church were comparable to those estimated at the Duomo in Milan and Santa Maria Maggiore in Rome, Italy. Moreover, it turned out that the identification of such gradients was restricted to a very limited period (August–September) during summertime. Furthermore, the results found in this study might provide guidelines for establishing a plan for thermal monitoring and preventive conservation in similar churches.
- The first methodology is based on Pearson’s correlation coefficient and linear regression. This methodology, which could help to determine reference thermal gradients for art conservation, could be improved using smoothing techniques and nonparametric regression. Furthermore, taking into account that datasets about indoor air conditions in historical buildings in Mediterranean climates are scarce, the confidence interval (95%) of the vertical gradient found in summer (0.030 C/m, 0.057 C/m), could be considered as a reference for further similar studies. Results obtained can be extrapolated to similar scenarios, whether in a heritage building or others, such as an industrial building, warehouse or farm of similar volume and height, with little ventilation, in a similar climate, according to some climate classification criteria (e.g., Köppen [114] and Trewartha [115]).
- The second methodology proposed here combines sPLS [85] and LDA. Furthermore, it employs variables computed from the seasonal H-W method, or functions that are applied to time series. This methodology helped to obtain parsimonious models with a small subset of variables, leading to satisfactory discrimination and easy interpretation of the different clusters of the time series. Furthermore, it was useful for identifying the most important variables for classifying time series. The variables computed from the seasonal H-W method yielded better results. In other studies, SH-W has also been shown to provide efficient results. This method was more flexible for fitting the distinct time series and obtaining low values of the classification error rate. The new methodology proposed allowed an efficient characterization of T at high, medium and low altitude levels. This approach had the best results according to the classification error rate and number of selected variables, when compared to results from SPLSDA [87] and sPLS-DA [35]. When using variables from seasonal H-W as input for either sPLS with LDA, sPLS-DA, or SPLSDA, both the error rate and the number of selected variables were better.
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Acknowledgments
Conflicts of Interest
References
- Perles, A.; Pérez-Marín, E.; Mercado, R.; Segrelles, J.D.; Blanquer, I.; Zarzo, M.; Garcia-Diego, F.J. An energy-efficient internet of things (IoT) architecture for preventive conservation of cultural heritage. Future Gener. Comput. Syst. 2018, 81, 566–581. [Google Scholar] [CrossRef]
- Varas-Muriel, M.J.; Fort, R.; Martínez-Garrido, M.I.; Zornoza-Indart, A.; López-Arce, P. Fluctuations in the indoor environment in Spanish rural churches and their effects on heritage conservation: Hygro-thermal and CO2 conditions monitoring. Build. Environ. 2014, 82, 97–109. [Google Scholar] [CrossRef] [Green Version]
- Pavlogeorgatos, G. Environmental parameters in museums. Build. Environ. 2003, 38, 1457–1462. [Google Scholar] [CrossRef]
- Camuffo, D.; Sturaro, G.; Valentino, A. Thermodynamic exchanges between the external boundary layer and the indoor microclimate at the basilica of Santa Maria Maggiore, Rome, Italy: The problem of conservation of ancient works of art. Bound. Layer Meteorol. 1999, 92, 243–262. [Google Scholar] [CrossRef]
- Tabunschikov, Y.; Brodatch, M. Indoor air climate requirements for Russian churches and cathedrals. Indoor Air J. 2004, 14 (Suppl. 7), 168–174. [Google Scholar] [CrossRef]
- Vuerich, E.; Malaspina, F.; Barazutti, M.; Georgiadis, T.; Nardino, M. Indoor measurements of microclimate variables and ozone in the church of San Vincenzo (Monastery of Bassano Romano Italy): A pilot study. Microchem. J. 2008, 88, 218–223. [Google Scholar] [CrossRef]
- Loupa, G.; Charpantidou, E.; Kioutsioukis, I.; Rapsomanikis, S. Indoor microclimate, ozone and nitrogen oxides in two medieval churches in Cyprus. Atmos. Environ. 2006, 40, 7457–7466. [Google Scholar] [CrossRef]
- Camuffo, D. Microclimate for Cultural Heritage, 1st ed.; Elsevier: Amsterdam, The Netherlands, 1998. [Google Scholar]
- Aste, N.; Adhikari, R.; Buzzetti, M.; Della Torre, S.; Del Pero, C.; Huerto, C.H.; Leonforte, F. Microclimatic monitoring of the duomo (Milan cathedral): Risks-based analysis for the conservation of its cultural heritage. Build. Environ. 2019, 148, 240–257. [Google Scholar] [CrossRef]
- Hall, C.; Hamilton, A.; Hoff, W.D.; Viles, H.A.; Eklund, J.A. Moisture dynamics in walls: Response to micro-environment and climate change. Proc. R. Soc. A Math. Phys. Eng. Sci. 2011, 467, 194–211. [Google Scholar] [CrossRef] [Green Version]
- Lucchi, E.; Roberti, F. Diagnosi, Simulazione e Ottimizzazione Energetica Degli Edifici Storici Energy Diagnosis, Simulation and Optimization of Historic Buildings. 2016. Available online: https://www.researchgate.net/profile/Elena_Lucchi3/publication/311846765_Diagnosi_simulazione_e_ottimizzazione_energetica_degli_edifici_storici_Energy_diagnosis_simulation_and_optimization_of_historic_buildings/links/585d308708ae6eb8719ff392/Diagnosi-simu (accessed on 3 October 2021).
- Adán, A.; Pérez, V.; Vivancos, J.L.; Aparicio-Fernández, C.; Prieto, S.A. Proposing 3D thermal technology for heritage building energy monitoring. Remote Sens. 2021, 3, 1537. [Google Scholar] [CrossRef]
- Camuffo, D.; Pagan, E.; Rissanen, S.; Bratasz, Ł.; Kozłowski, R.; Camuffo, M.; della Valle, A. An advanced church heating system favourable to artworks: A contribution to European standardisation. J. Cult. Herit. 2010, 11, 205–219. [Google Scholar] [CrossRef]
- Muñoz-González, C.; León-Rodríguez, A.L.; Navarro-Casas, J. Air conditioning and passive environmental techniques in historic churches in Mediterranean climate. A proposed method to assess damage risk and thermal comfort pre-intervention, simulation-based. Energy Build. 2016, 130, 567–577. [Google Scholar] [CrossRef]
- EN15757. Conservation of Cultural Property. Specifications for Temperature and Relative Humidity to Limit Climate-Induced Mechanical Damage in Organic Hygroscopic Materials. 2010. Available online: https://standards.iteh.ai/catalog/standards/cen/ad03d50b-22dc-4c57-b198-2321863f3870/en-15757-2010 (accessed on 3 October 2021).
- Mecklenburg, M.F.; Tumosa, C.S.; Erhardt, W.D. Structural Response of Painted Wood Surfaces to Changes in Ambient Relative Humidity in Painted Wood: History and Conservation; Dorge, V., Howlett, F.C., Eds.; The Getty Conservation Institute: Los Angeles, CA, USA, 1998. [Google Scholar]
- Bratasz, L.; Koztowski, R.; Camuffo, D.; Pagan, E. Impact of indoor heating on painted wood—Monitoring the altarpiece in the church of Santa Maria Maddalena in Rocca Pietore, Italy. Stud. Conserv. 2007, 52, 199–210. [Google Scholar] [CrossRef]
- Bratasz, L.; Camuffo, D.; Kozlowski, R. Target microclimate for preservation derived from past indoor conditions. In Contributions to the Museum Microclimates Conference; The National Museum of Denmark: Copenhagen, Denmark, 2007; pp. 129–134. [Google Scholar]
- Silva, H.E.; Henriques, F.M. Microclimatic analysis of historic buildings: A new methodology for temperate climates. Build. Environ. 2014, 82, 381–387. [Google Scholar] [CrossRef]
- Weatherbase. Weatherbase. 2021. Available online: https://www.weatherbase.com/weather/weather-summary.php3?s=48280&cityname=Valencia,+Spain (accessed on 3 October 2021).
- Collectioncare.eu. European Horizon 2020 Project Collectioncare: Innovative and Affordable Service for the Preventive Conservation Monitoring of Individual Cultural Artefacts during Display, Storage, Handling and Transport. 2020. Available online: https://www.collectioncare.eu/ (accessed on 3 October 2021).
- Corgnati, S.P.; Filippi, M. Assessment of thermo-hygrometric quality in museums: Method and in-field application to the Duccio di Buoninsegna exhibition at Santa Maria della Scala (Siena, Italy). J. Cult. Herit. 2010, 11, 345–349. [Google Scholar] [CrossRef]
- García-Diego, F.J.; Verticchio, E.; Beltrán, P.; Siani, A.M. Assessment of the minimum sampling frequency to avoid measurement redundancy in microclimate field surveys in museum buildings. Sensors 2016, 16, 1291. [Google Scholar] [CrossRef]
- Herráez, J.; Enríquez de Salamanca, M.; Pastor Arenas, T.; Gil Muñoz, M. Manual de Seguimiento y Análisis de Condiciones Ambientales. Plan Nacional de Conservación Preventiva. 2014. Available online: https://www.libreria.culturaydeporte.gob.es/libro/manual-de-seguimiento-y-analisis-de-condiciones-ambientales_2654/ (accessed on 3 October 2021).
- Sesana, E.; Gagnon, A.S.; Bertolin, C.; Hughes, J. Adapting cultural heritage to climate change risks: Perspectives of cultural heritage experts in Europe. Geosciences 2018, 8, 305. [Google Scholar] [CrossRef] [Green Version]
- Angelini, E.; Grassini, S.; Corbellini, S.; Parvis, M.; Piantanida, M. A multidisciplinary approach for the conservation of a building of the seventeenth century. Appl. Phys. A 2010, 100, 763–769. [Google Scholar] [CrossRef]
- Lourenço, P.B.; Luso, E.; Almeida, M.G. Defects and moisture problems in buildings from historical city centres: A case study in Portugal. Build. Environ. 2006, 41, 223–234. [Google Scholar] [CrossRef] [Green Version]
- Ramírez, S.; Zarzo, M.; Perles, A.; Garcia-Diego, F.J. Methodology for discriminant time series analysis applied to microclimate monitoring of fresco paintings. Sensors 2021, 21, 436. [Google Scholar] [CrossRef]
- García-Diego, F.J.; Zarzo, M. Microclimate monitoring by multivariate statistical control: The renaissance frescoes of the cathedral of Valencia (Spain). J. Cult. Herit. 2010, 11, 339–344. [Google Scholar] [CrossRef]
- Zarzo, M.; Fernández-Navajas, A.; García-Diego, F.J. Long-term monitoring of fresco paintings in the cathedral of Valencia (Spain) through humidity and temperature sensors in various locations for preventive conservation. Sensors 2011, 11, 8685–8710. [Google Scholar] [CrossRef] [PubMed]
- Merello, P.; García-Diego, F.J.; Zarzo, M. Microclimate monitoring of Ariadne’s house (Pompeii, Italy) for preventive conservation of fresco paintings. Chem. Cent. J. 2012, 6, 145. [Google Scholar] [CrossRef] [Green Version]
- Muñoz-González, C.M.; León-Rodríguez, Á.L.; Suárez Medina, R.C.; Teeling, C. Hygrothermal performance of worship spaces: Preservation, comfort, and energy consumption. Sustainability 2018, 10, 3838. [Google Scholar] [CrossRef] [Green Version]
- Ramírez, S.; Zarzo, M.; García-Diego, F.J. Multivariate time series analysis of temperatures in the archaeological museum of L’Almoina (Valencia, Spain). Sensors 2021, 21, 4377. [Google Scholar] [CrossRef] [PubMed]
- Klein, L.J.; Bermudez, S.A.; Schrott, A.G.; Tsukada, M.; Dionisi-Vici, P.; Kargere, L.; Marianno, F.; Hamann, H.F.; López, V.; Leona, M. Wireless sensor platform for cultural heritage monitoring and modeling system. Sensors 2017, 17, 1998. [Google Scholar] [CrossRef] [Green Version]
- Lê Cao, K.A.; Boitard, S.; Besse, P. Sparse PLS discriminant analysis: Biologically relevant feature selection and graphical displays for multiclass problems. J. Clin. Bioinform. 2011, 12, 253. [Google Scholar] [CrossRef] [Green Version]
- Camuffo, D. The friendly heating project and the conservation of the cultural heritage preserved in churches. In Proceedings of the Developments in Climate Control of Historic Buildings. Proceedings from the International Conference “Climatization of Historic Buildings, State of the Art”, Linderhof Palace, Ettal, Germany, 2 December 2010; pp. 699–710. [Google Scholar]
- García-Diego, F.; Navajas, Á.F.; Beltrán, P.; Merello, P. Study of the effect of the strategy of heating on the Mudejar Church of Santa Maria in Ateca (Spain) for preventive conservation of the Altarpiece Surroundings. Sensors 2013, 13, 11407–11423. [Google Scholar] [CrossRef] [Green Version]
- Al-Omari, A.; Brunetaud, X.; Beck, K.; Muzahim, A. Effect of thermal stress, condensation and freezing–thawing action on the degradation of stones on the Castle of Chambord, France. Environ. Earth Sci. 2014, 71, 3977–3989. [Google Scholar] [CrossRef]
- Brunetaud, X.; De Luca, L.; Janvier-Badosa, S.; Beck, K.; Al-Mukhtar, M. Application of digital techniques in monument preservation. Eur. J. Environ. Civ. Eng. 2012, 16, 543–556. [Google Scholar] [CrossRef]
- USA, O. Hobo Data Loggers. 2021. Available online: http://www.onsetcomp.com/ (accessed on 3 October 2021).
- Visco, G.; Plattner, S.; Fortini, P.; Di Giovanni, S.; Sammartino, M.P. Microclimate monitoring in the Carcer Tullianum: Temporal and spatial correlation and gradients evidenced by multivariate analysis; first campaign. Chem. Cent. J. 2012, 6, 3977–3989. [Google Scholar] [CrossRef] [PubMed]
- USA, O. Temperature Logger iButton with 8KB Data-Log Memory. 2015. Available online: https://datasheets.maximintegrated.com/en/ds/DS1922L-DS1922T.pdf (accessed on 3 October 2021).
- Valero, M.A.; Merello, P.; Navajas, A.F.; Garcia-Diego, F.J. Statistical tools applied in the characterisation and evaluation of a thermo-hygrometric corrective action carried out at the Noheda archaeological site (Noheda, Spain). Sensors 2014, 14, 1665–1679. [Google Scholar] [CrossRef] [PubMed]
- iButton. Hygrochron Temperature/Humidity Logger iButton with 8KB Data-Log Memory. 2021. Available online: https://ibutton.cl/product/datalogger-temperatura-humedad-hygrochron-ds1923/ (accessed on 3 October 2021).
- Merello, P.; Fernandez-Navajas, A.; Curiel-Esparza, J.; Zarzo, M.; García-Diego, F.J. Characterisation of thermo-hygrometric conditions of an archaeological site affected by unlike boundary weather conditions. Build Environ. 2014, 76, 125–133. [Google Scholar] [CrossRef]
- Merello, P.; García-Diego, F.J.; Zarzo, M. Diagnosis of abnormal patterns in multivariate microclimate monitoring: A case study of an open-air archaeological site in Pompeii (Italy). Sci. Total Environ. 2014, 488–489, 14–25. [Google Scholar] [CrossRef]
- García-Diego, F.J.; Borja, E.; Merello, P. Design of a hybrid (Wired/Wireless) acquisition data system for monitoring of cultural heritage physical parameters in smart cities. Sensors 2015, 15, 7246–7266. [Google Scholar] [CrossRef]
- Best, D.J.; Roberts, D.E. Algorithm AS 89: The upper tail probabilities of Spearman’s rho. J. R. Stat. Soc. Ser. C Appl. Stat. 1975, 24, 377–379. [Google Scholar] [CrossRef]
- Hastie, T. Statistical Models in S; Chapter 4; Routledge: Oxfordshire, UK, 1992. [Google Scholar]
- Richards, J.W.; Starr, D.L.; Butler, N.R.; Bloom, J.S.; Brewer, J.M.; Crellin-Quick, A.; Higgins, J.; Kennedy, R.; Rischard, M. On machine-learned clasification of variable stars with sparse and noisy time-series data. Astrophys. J. 2011, 733, 10. [Google Scholar] [CrossRef]
- Chung, D.; Keles, S. Sparse partial least squares classification for high dimensional data. Stat. Appl. Genet. Mol. Biol. 2010, 9, 1–30. [Google Scholar] [CrossRef]
- Vilar, J.; Montero, P. TSclust: An R package for time series clustering. J. Stat. Softw. 2014, 62, 1–43. [Google Scholar] [CrossRef] [Green Version]
- Batista, E.; Wang, X.; Keogh, E. A complexity-invariant distance measure for time series. In Proceedings of the Eleventh SIAM International Conference on Data Mining, SDM11, Mesa, AZ, USA, 28–30 April 2011; pp. 699–710. [Google Scholar] [CrossRef] [Green Version]
- Piccolo, D. A distance measure for classifying ARIMA models. J. Time Ser. Anal. 1990, 11, 153–164. [Google Scholar] [CrossRef]
- Elorrieta, F.; Eyheramendy, S.; Jordán, A.; Dékány, I.; Catelan, M.; Angeloni, R.; Alonso-García, J.; Contreras-Ramos, R.; Gran, F.; Hajdu, G.; et al. A machine learned classifier for RR Lyrae in the VVV survey. A&A 2016, 595, A82. [Google Scholar] [CrossRef]
- Elorrieta, F. Classification and Modeling of Time Series of Astronomical Data. Ph.D. Thesis, Pontificia Universidad Catolica de Chile, Santiago, Chile, 2018. [Google Scholar]
- The Sensor Company. Sensirion. 2011. Available online: https://www.sensirion.com/fileadmin/user_upload/customers/sensirion/Dokumente/2_Humidity_Sensors/Datasheets/Sensirion_Humidity_Sensors_SHT1x_Datasheet.pdf (accessed on 3 October 2021).
- Perles, A.; Mercado, R.; Capella, J.V.; Serrano, J.J. Ultra-low power optical sensor for xylophagous insect detection in wood. Sensors 2016, 16, 1977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Standard, O. MQTT: The Standard for IoT Messaging OASIS. 2019. Available online: https://docs.oasis-open.org/mqtt/mqtt/v5.0/mqtt-v5.0.html (accessed on 3 October 2021).
- Redash. 2021. Available online: https://redash.io (accessed on 3 October 2021).
- Visco, G.; Plattner, S.; Fortini, P.; Sammartino, M.P. A multivariate approach for a comparison of big data matrices. Case study: Thermo-hygrometric monitoring inside the Carcer Tullianum (Rome) in the absence and in the presence of visitors. Environ. Sci. Pollut. Res. 2017, 24, 3977–3989. [Google Scholar] [CrossRef] [PubMed]
- Little, R.J.; Rubin, D.B. Statistical Analysis with Missing Data, 2nd ed.; Wiley: Hoboken, NJ, USA, 2002. [Google Scholar]
- Moritz, S.; Bartz-Beielstein, T. imputeTS: Time Series Missing Value Imputation in R. R J. 2017, 9, 207–218. [Google Scholar] [CrossRef] [Green Version]
- R Core Team. R: A Language and Environment for Statistical Computing. 2014. Available online: https://www.r-project.org/about.html (accessed on 3 October 2021).
- Rohart, F.; Gautier, B.; Singh, A.; Lê Cao, K.A. mixOmics: An R package for omics feature selection and multiple data integration. PLoS Comput. Biol. 2017, 13, e1005752. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Roever, C.; Raabe, N.; Luebke, K.; Ligges, U.; Szepannek, G.; Zentgraf, M.; Meyer, D. Package `klaR: Classification and Visualization’. 2020. Available online: https://cran.r-project.org/web/packages/klaR/index.html (accessed on 3 October 2021).
- Chung, D.; Chun, H.; Keles, S. Package Spls: Sparse Partial Least Squares (SPLS) Regression and Classification. 2019. Available online: https://cran.r-project.org/web/packages/spls/index.html (accessed on 3 October 2021).
- Leisch, F.; Hornik, K.; Kuan, C.M. Monitoring structural changes with the generalized fluctuation test. J. Econ. Theory 2000, 16, 835–854. [Google Scholar] [CrossRef] [Green Version]
- Zeileis, A.; Leisch, F.; Hornik, K.; Kleiber, C. Strucchange: An R package for testing for structural change in linear regression models. J. Stat. Softw. 2002, 7, 1–38. [Google Scholar] [CrossRef] [Green Version]
- Hamilton, J.D. Time Series Analysis, 5th ed.; Princeton University Press: Princeton, NJ, USA, 1994. [Google Scholar]
- Venables, W.N.; Ripley, B.D. Modern Applied Statistics with S, 4th ed.; Springer: New York, NY, USA, 2002. [Google Scholar]
- Brockwell, P.J.; Davis, R.A. Time Series: Theory and Methods, 2nd ed.; Springer: New York, NY, USA, 1987. [Google Scholar]
- Winters, P.R. Forecasting sales by exponentially weighted moving averages. Manag. Sci. 1960, 6, 324–342. [Google Scholar] [CrossRef]
- Holt, C.C. Forecasting seasonals and trends by exponentially weighted moving averages. Int. J. Forecast. 2004, 20, 5–10. [Google Scholar] [CrossRef]
- Hyndman, R.; Athanasopoulos, G. Forecasting: Principles and Practice. OTexts. 2013. Available online: http://otexts.org/fpp/ (accessed on 3 October 2021).
- R Core Team and Contributors Worldwide. Package `Stats’. 2020. Available online: https://www.rdocumentation.org/packages/stats/versions/3.6.2 (accessed on 3 October 2021).
- Conover, W.J. Practical Nonparametric Statistics, 3rd ed.; John Wiley and Sons, INC: Hoboken, NJ, USA, 1999. [Google Scholar]
- Royston, J.P. Algorithm AS 181: The W test for normality. J. R. Stat. Soc. Ser. C Appl. Stat. 1982, 31, 176–180. [Google Scholar] [CrossRef]
- Box, G.E.P.; Cox, D.R. An analysis of transformations. J. R. Stat. Soc. Ser. B 1964, 26, 211–243. [Google Scholar] [CrossRef]
- Bacon, C. Practical Portfolio Performance Measurement and Attribution, 2nd ed.; John Wiley and Sons: New Delhi, India, 2008. [Google Scholar]
- Wold, H. Multivariate Analysis; Wiley: New York, NY, USA, 1966. [Google Scholar]
- Tenenhaus, M. La Régression PLS: Thórie et Pratique; Editions Technip: Paris, France, 1998. [Google Scholar]
- Dray, S.; Pettorelli, N.; Chessel, D. Multivariate analysis of incomplete mapped data. Trans. GIS 2003, 7, 411–422. [Google Scholar] [CrossRef]
- Jolliffe, I.T.; Trendafilov, N.T.; Uddin, M. A modified principal component technique based on the LASSO. J. Comput. Graph. Stat. 2003, 12, 531–547. [Google Scholar] [CrossRef] [Green Version]
- Lê Cao, K.A.; Rossouw, D.; Robert-Granié, C.; Besse, P. A sparse PLS for variable selection when integrating omics data. Stat. Appl. Genet. Mol. Biol. 2008, 7, 1–29. [Google Scholar] [CrossRef] [PubMed]
- Lê Cao, K.A.; Martin, P.G.P.; Robert-Granie, C.; Besse, P. Sparse canonical methods for biological data integration: Application to a cross-platform study. J. Clin. Bioinform. 2009, 10, 34. [Google Scholar] [CrossRef]
- Chun, H.; Keleş, S. Sparse partial least squares regression for simultaneous dimension reduction and variable selection. J. R. Stat. Soc. Ser. B 2010, 72, 3–25. [Google Scholar] [CrossRef] [Green Version]
- Tibshirani, R. Regression shrinkage and selection via the Lasso. J. R. Stat. Soc. Ser. B 1996, 58, 267–288. [Google Scholar] [CrossRef]
- Hastie, T.; Tibshirani, R.; Friedman, J.; Freedman, J. The Elements of Statistical Learning: Data Mining, Inference and Prediction, 2nd ed.; Springer: Standford, CA, USA, 2009. [Google Scholar]
- Shen, H.; Huang, J.Z. Sparse principal component analysis via regularized low rank matrix approximation. J. Multivar. Anal. 2008, 99, 1015–1034. [Google Scholar] [CrossRef] [Green Version]
- Lê Cao, K.A.; Dejean, S. mixOmics: Omics Data Integration Project. 2020. Available online: https://bioconductor.org/packages/release/bioc/html/mixOmics.html (accessed on 3 October 2021). [CrossRef]
- Mevik, B.H.; Wehrens, R. The pls package: Principal component and partial least squares regression in R. J. Stat. Softw. Artic. 2007, 18, 1–23. [Google Scholar] [CrossRef] [Green Version]
- Sartorius Stedim Data Analytics AB. SIMCA 15 Multivariate Data Analysis Solution User Guide; Sartorius Stedim Data Analytics AB: Umeå, Sweden, 2017. [Google Scholar]
- Ramírez, S.; Zarzo, M.; Perles, A.; Garcia-Diego, F.J. sPLS-DA to discriminate time series. In JSM Proceedings, Statistics and the Environment JSM2020; American Statistical Association: Alexandria, VA, USA, 2021; pp. 107–135. [Google Scholar]
- Kuhn, M.; Wing, J.; Weston, S.; Williams, A. Package `Caret’. 2021. Available online: https://cran.r-project.org/web/packages/caret/caret.pdf (accessed on 3 October 2021).
- Camuffo, D. Microclimate for Cultural Heritage Conservation, Restoration, and Maintenance of Indoor and Outdoor Monuments, 2nd ed.; Elsevier Science: Amsterdam, The Netherlands, 2014. [Google Scholar]
- Jiang, Y.; Lu, L. A study of dust accumulating process on solar photovoltaic modules with different surface temperatures. Energy Procedia 2015, 75, 337–342. [Google Scholar] [CrossRef] [Green Version]
- Hastie, T.; Tibshirani, R. Generalized Additive Models; Chapman & Hall/CRC Monographs on Statistics & Applied Probability; Taylor & Francis: London, UK; Boca Raton, FL, USA; New York, NY, USA; Washington, DC, USA, 1990. [Google Scholar]
- EN16883. Conservation of Cultural Heritage. Guidelines for Improving the Energy Performance of Historic Buildings. 2017. Available online: https://standards.iteh.ai/catalog/standards/cen/189eac8d-14e1-4810-8ebd-1e852b3effa3/en-16883-2017 (accessed on 3 October 2021).
- EN16141. Conservation of Cultural Heritage. Guidelines for Management of Environmental Conditions. Open Storage Facilities: Definitions and Characteristics of Collection Centres Dedicated to the Preservation and Management of Cultural Heritage. 2012. Available online: https://standards.iteh.ai/catalog/standards/cen/452b6461-3ba0-4f70-8007-f6702a7e804d/en-16141-2012 (accessed on 3 October 2021).
- EN16242. Conservation of Cultural Heritage. Procedures and Instruments for Measuring Humidity in the Air and Moisture Exchanges between Air and Cultural Property. 2012. Available online: https://standards.iteh.ai/catalog/standards/cen/12d5d8fb-f047-4d49-a750-9349dff5f7f4/en-16242-2012 (accessed on 3 October 2021).
- EN15898. Conservation of Cultural Property. Main General Terms and Definitions. 2019. Available online: https://standards.iteh.ai/catalog/standards/cen/d8ac0bbb-7dfb-4c68-a51b-9c0b66ae8f63/en-15898-2011 (accessed on 3 October 2021).
- EN15758. Conservation of Cultural Property. Procedures and Instruments for Measuring Temperatures of the Air and the Surfaces of Objects. 2010. Available online: https://standards.iteh.ai/catalog/standards/cen/8e796a23-395f-4ed2-8154-33ef479b3231/en-15758-2010 (accessed on 3 October 2021).
- EN16893. Conservation of Cultural Heritage. Specifications for Location, Construction and Modification of Buildings or Rooms Intended for the Storage or Use of Heritage Collections. 2018. Available online: https://standards.iteh.ai/catalog/standards/cen/b53c7f39-b809-4998-be91-31f362a8a2cb/en-16893-2018 (accessed on 3 October 2021).
- UNE-EN-15759-1. Conservación del Patrimonio Cultural. Clima Interior. Parte 1: Recomendaciones Para la Calefacción de Iglesias, Capillas y Otros Lugares de Culto. 2012. Available online: https://www.une.org/encuentra-tu-norma/busca-tu-norma/norma?c=N0050451 (accessed on 3 October 2021).
- UNE-EN 15757. Conservación del Patrimonio Cultural. Especificaciones de Temperatura y Humedad Relativa para Limitar los Daños Mecánicos Causados por el Clima a Los Materiales Orgánicos Higroscópicos. 2011. Available online: https://tienda.aenor.com/norma-une-en-15757-2011-n0046628 (accessed on 3 October 2021).
- ASHRAE Standard 55-2007, A.S. Thermal Environment Conditions for Human Occupancy. 2020. Available online: https://www.ashrae.org/technical-resources/bookstore/standard-55-thermal-environmental-conditions-for-human-occupancy (accessed on 3 October 2021).
- UNI 10829. Works of Art of Historical Importance. Ambient Conditions for the Conservation. Measurement and Analysis. UNI Ente Nazionale Italiano di Unificazione. 1999. Available online: http://www.brescianisrl.it/newsite/ita/xprodotti.php?sub_content=1&sub_prodotto=47&hash=db99bf602ca201c16ad187434828d334 (accessed on 13 October 2021).
- Andretta, M.; Coppola, F.; Modelli, A.; Santopuoli, N.; Seccia, L. Proposal for a new environmental risk assessment methodology in cultural heritage protection. J. Cult. Herit. 2017, 23, 22–32. [Google Scholar] [CrossRef]
- Fabbri, K.; Bonora, A. Two new indices for preventive conservation of the cultural heritage: Predicted risk of damage and heritage microclimate risk. J. Cult. Herit. 2021, 47, 208–217. [Google Scholar] [CrossRef]
- García-Diego, J.F.; Scatigno, C.; Merello, P.; Bustamante, E. Preliminary data of CFD modeling to assess the ventilation in an archaelogical building. In Proceedings of the 8th International Congress on Archaeology, Computer Graphics, Cultural Heritage and Innovation ‘ARQUEOLÓGICA 2.0’, Valencia, Spain, 5–7 September 2016; pp. 504–507. [Google Scholar]
- Albero, S.; Giavarini, C.; Santarelli, M.L.; Vodret, A. CFD modeling for the conservation of the Gilded Vault Hall in the Domus Aurea. J. Cult. Herit. 2004, 5, 197–203. [Google Scholar] [CrossRef]
- Oetelaar, T. CFD, thermal environments, and cultural heritage: Two case studies of Roman baths. In Proceedings of the 2016 IEEE 16th International Conference on Environment and Electrical Engineering (EEEIC), Florence, Italy, 7–10 June 2016; pp. 1–6. [Google Scholar] [CrossRef]
- Handbuch der Klimatologie; Chapter Das Geographisca System der Klimate; Köppen, W. (Ed.) Verlag Von Gebrüder Borntraeger: Berlin, Germany, 1936; pp. 1–44. [Google Scholar]
- Trewartha, G.; Horn, L. An Introduction to Climate, 5th ed.; McGraw-Hill: New York, NY, USA, 1980. [Google Scholar]
Node | B | T | U | S | R | C | D | G | E | O | K |
---|---|---|---|---|---|---|---|---|---|---|---|
Bias | −0.280 | 0.097 | 0.160 | −0.003 | 0.069 | −0.088 | 0.077 | 0.009 | −0.019 | −0.089 | −0.036 |
Node | N | L | M | I | J | Q | A | F | P | H | |
Bias | −0.046 | −0.249 | 0.000 | 0.150 | 0.189 | −0.277 | −0.098 | 0.276 | 0.335 | 0.175 |
(a) | (b) | |||||||
---|---|---|---|---|---|---|---|---|
Stage | V | r | p-Value | Stage | V | r | p-Value | |
1 | Stg1 | mean.ts | 0.86 | 0.000 | Stg4 | a | −0.33 | 0.138 |
2 | Stg5 | pacf4 | 0.28 | 0.217 | Stg1 | s7 | 0.83 | 0.000 |
3 | Stg4 | pacf4 | 0.36 | 0.108 | Stg1 | s8 | 0.79 | 0.000 |
4 | Stg7 | pacf3 | 0.51 | 0.019 | Stg1 | s6 | 0.80 | 0.000 |
5 | Stg3 | pacf4 | 0.43 | 0.054 | Stg1 | s19 | −0.77 | 0.000 |
6 | Stg1 | mean.mr | −0.65 | 0.001 | Stg5 | a | −0.31 | 0.178 |
7 | Stg5 | pacf3 | 0.31 | 0.167 | Stg4 | s12 | 0.74 | 0.000 |
8 | Stg1 | s18 | −0.75 | 0.000 | ||||
9 | Stg4 | s6 | −0.23 | 0.319 | ||||
10 | Stg1 | a | 0.77 | 0.000 | ||||
11 | Stg7 | s20 | −0.41 | 0.065 | ||||
12 | Stg7 | s16 | 0.70 | 0.000 | ||||
13 | Stg4 | s23 | −0.69 | 0.001 |
(a) | (b) | |||
---|---|---|---|---|
Classification Method | Error Rate (%) | N | Error Rate (%) | N |
sPLS-DA | 35.06 | 10 | 18.75 | 15 |
SPLSDA | 19.04 | 42 | 19.04 | 11 |
sPLS [85] with LDA | 14.28 | 15 | 4.76 | 15 |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Ramírez, S.; Zarzo, M.; Perles, A.; García-Diego, F.-J. Characterization of Temperature Gradients According to Height in a Baroque Church by Means of Wireless Sensors. Sensors 2021, 21, 6921. https://doi.org/10.3390/s21206921
Ramírez S, Zarzo M, Perles A, García-Diego F-J. Characterization of Temperature Gradients According to Height in a Baroque Church by Means of Wireless Sensors. Sensors. 2021; 21(20):6921. https://doi.org/10.3390/s21206921
Chicago/Turabian StyleRamírez, Sandra, Manuel Zarzo, Angel Perles, and Fernando-Juan García-Diego. 2021. "Characterization of Temperature Gradients According to Height in a Baroque Church by Means of Wireless Sensors" Sensors 21, no. 20: 6921. https://doi.org/10.3390/s21206921
APA StyleRamírez, S., Zarzo, M., Perles, A., & García-Diego, F. -J. (2021). Characterization of Temperature Gradients According to Height in a Baroque Church by Means of Wireless Sensors. Sensors, 21(20), 6921. https://doi.org/10.3390/s21206921