# Improvement of the Non-Destructive Testing of Heritage Mural Paintings Using Stimulated Infrared Thermography and Frequency Image Processing

Kamel Mouhoubi, Vincent Detalle, Jean-Marc Vallet, Jean-Luc Bodnar

► To cite this version:

Kamel Mouhoubi, Vincent Detalle, Jean-Marc Vallet, Jean-Luc Bodnar. Improvement of the Non-Destructive Testing of Heritage Mural Paintings Using Stimulated Infrared Thermography and Frequency Image Processing. Journal of Imaging, 2019, 5 (9), pp.72. 10.3390/jimaging5090072. hal-03164311

HAL Id: hal-03164311

https://hal.science/hal-03164311

Submitted on 29 Aug 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Improvement of the Non-Destructive Testing of Heritage Mural Paintings Using Stimulated Infrared Thermography and Frequency Image Processing

Kamel Mouhoubi <fnref n="1" />,* and Vincent Detalle <fnref n="2" />, Jean-Marc Vallet <fnref n="3" /> and Jean-Luc Bodnar <fnref n="1" />

Received: 29 April 2019; Accepted: 18 August 2019; Published: 29 August 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Within the framework of conservation and assistance for the restoration of cultural property, a method of analysis assistance has been developed to help in the restoration of cultural heritage. Several collaborations have already demonstrated the possibility of defects detection (delamination, salts) in murals paintings using stimulated infrared thermography. One of the difficulties encountered with infrared thermography applied to the analysis of works of art is the remanence of the pictorial layer. This difficulty can sometimes induce detection artifacts and false positives. A method of thermograms post-processing called PPT (pulse phase thermography) is described. The possibilities offered by the PPT in terms of reducing the optical effects associated with the pictorial layer are highlighted first with a simulation, and then through experiments. This approach can significantly improve the study of painted works of art such as wall paintings.

**Keywords:** stimulated infrared thermography; work of art; post-processing; Fourier transform; image enhancement; NDT

## 1. Introduction

Heritage conservation is an important societal challenge. The uniqueness and value of cultural heritage involves as much effort as possible to implement non-destructive testing techniques for diagnosis purposes. A lot of them are tested, including stimulated infrared thermography [1–15]. “Laboratoire de Recherche des Monuments Historiques” (LRMH), “Centre Interdisciplinaire de Conservation et Restauration du Patrimoine” (CICRP) and “Institut de Thermique, Mécanique, Matériaux” (ITheMM) of the Reims University (URCA) are thus collaborating in order to develop this last one. First results showed that stimulated infrared thermography allows detecting detachments within wall paintings [7,8]. However, the sensitivity of this method to the thermal response of pigments can lead to induce artifacts in the detection of defects. This impact can be reduced in our point of view by the way of Pulse Phase Thermography analysis.

This method is based on discrete Fourier transform analysis of the thermographic film obtained. It gives access to two characteristic parameters: The amplitude and the phase of the photothermal signal. The amplitude is potentially sensitive to heterogeneities of energy deposition on it. The phase measures the delay of the photothermal response in relation to the visible excitation as a function of frequency and it is thus potentially less sensitive to surface radiative issues.

Our work has been focused on this last property according to two research axes: We first made a simulation based on finite element modelling of the photothermal experiment which has been developed on an academic sample. We also developed an experimental study on two laboratory samples. The obtained results are presented according to the following four sections:

<footnote n="1">¹ ITHEMM, UFR Sciences Exactes et Naturelles, BP 1039, 51687 Reims CEDEX 02, France</footnote>

<footnote n="2">² C2RMF, 14 Quai François Mitterrand, 75001 Paris, France</footnote>

<footnote n="3">³ CICRP, 21 rue Guibal, 13003 Marseille, France</footnote>

<footnote>* Correspondence: kmouhoubi@gmail.com; Tel.: +33-6-669-722-90</footnote>
First, the principle of the implemented photothermal model is developed. In a second step, we present the obtained simulation results and we show that the phase parameter is interesting in terms of rejection of the optical effects induced by the different colors of the pictorial layer. In a third step, the experimental study is described: We present the samples and we describe the experimental system implemented. Finally, we discuss the obtained experimental results and their match to the theoretical conclusions.

## 2. The Photothermal Model Developed for the Experiment

### 2.1. Principle of the Pulse Phase Thermography (PPT)

The principle of stimulated infrared thermography is divided into three steps: First, the sample under analysis is excited by a luminous flux. This energy input causes an increase of the temperature on the surface of the sample. The heat is then diffused into the analyzed sample. The presence of a defect disturbs the thermal diffusion, and leads to a slower cooling of the analyzed surface in the case of a displacement. This last change can be detected by an infrared thermography camera.

In the case of the PPT in infrared thermography, the excitation signal is of the crenel type. The thermographic film is then subjected to a mathematical post-processing, such as discrete Fourier transform (DFT), that is defined by the following expression:

$$F(u) = \frac{1}{N} \sum_{n=0}^{N-1} h(x) \exp \left[ -\frac{j2\pi ux}{N} \right] = R(u) + jI(u) \quad (1)$$

When $T(t)$ is a measured thermal signal at a pixel, after sampling it gives a succession of measurement points of index $x$, $h(x)$ represents the temperature for this pixel at the index $x$, $N$ is the number of measurement points resulting from the sampling and $u$ is the frequency calculated by applying Expression (1). Maldague et al. [16–18] introduced this mode of analysis for the thermography. It allows, while being quick (crenel excitation), access to the characteristic parameters provided by a sinusoidal excitation: the amplitude and the phase of the photothermal signal from the real and imaginary components $R(u)$ and $I(u)$. They are obtained from Equation (1). Refer to the books [19,20] on the discrete Fourier transform for more information.

The sequential processing of the thermogram is made pixel by pixel in a fixed time interval. It makes it possible to extract an image of amplitude and phase for each frequency component $u$.

Initially, this variant had the objective of a good noise rejection and consequently a detection of low photothermal signatures. In our case, our interest is about the temporal connotation of the phase. Because it is less sensitive to energy fluctuations, such as heterogeneities of energy deposition. These heterogeneities are induced by different radiative absorptivities coming from the pictorial layer that are present in the raw thermographic film. So, it enables us to detect the defects better.

### 2.2. Simulation Developed for the Study

The simulation developed for the study is based on solving the heat equation by the finite element method. The process is performed in two steps: first, the calculation engine determines numerically the spatial and temporal distribution of temperature, and second, we collect these data and perform a discrete Fourier transform (DFT) on the surface temperatures of the sample. This finally results in amplitude and phase maps of photothermal signal as a function of frequency.

The sample we analyzed theoretically (Figure 1) is a block of plaster which has thermophysical properties very similar to those of a mural painting. A Figure 2 shows the meshing of the sample studied, a fine mesh has been chosen.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** The studied sample: Overview (a) and view from the side (b) of the theoretically studied sample.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Meshing of the sample studied.</figcaption></figure>

It is 160 mm by 120 mm in size and 20 mm thick. The retained thermophysical properties are as follows: a thermal conductivity of 0.4 W/mK, a density of 1100 kg/m³, a calorific capacity of 830 J/kg K, which gives a thermal diffusivity of $4.38 \times 10^{-7}$ m²/s. Six parallelepiped spaces containing air (detachments) were arranged in the volume of this sample. Their dimensions are 20 mm in length and in width and 4 mm in thickness. Their depths vary from 2 mm to 12 mm in 2 mm increments.

The thermophysical properties taken into account for these defects are those of the air at 20 °C: a thermal conductivity of 0.026 W/mK, a density of 1.17 kg/m³, a calorific capacity of 1006 J/kg K or a thermal diffusivity of $2.22 \times 10^{-5}$ m²/s.

Finally, two areas were drawn on the surface of the sample in order to simulate the optical effects induced by the pictorial layer: The first area has an absorptivity 1.5 times higher than the second area and covers defects located at 4 mm, 8 mm and 12 mm deep. The second zone covers defects located at 2 mm, 6 mm and 10 mm in depth. The type of excitation signal to the sample is a crenel of 2 seconds duration (Figure 3). The analysis time is 200 s. The acquisition frequency is fixed at 1 Hertz. The power of energy used is 1500 W. The energy distribution was adjusted as follows: the zone of lower absorptivity received a total power of 1000 W while the zone of higher absorptivity received a total power of 1500 W. We considered a no-loss thermal model.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Excitation signal.</figcaption></figure>

**3. Results Obtained**

In Figure 4, we present the obtained raw thermograms at t = 28 s, t = 91 s, t = 177 s and t = 200 s.
They show a higher photothermal heating at the location of the defects.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Example of theoretical thermograms obtained (t = 28 s, t = 91 s, t = 177 s and t = 200 s).</figcaption></figure>

The intensity of the background is variable and comparable defects in terms of depth do not appear with the same contrast. The Figure 4 shows that differences in surface absorptivity can lead to the creation of detection artifacts. We are trying to reduce this phenomenon.

In Figure 5, we then present the results obtained after a PPT-type frequency analysis. It shows the obtained module images, corresponding respectively to frequencies of 5 mHz, 10 mHz, 20 mHz, 30 mHz, 40 mHz and 80 mHz.

This figure shows that frequency images of amplitude, like raw signals, can detect all 6 defects. It also shows that this detection is biased by the inhomogeneities of energetic deposition due to the pictorial layer; defects at near depths do not appear with the same intensity. As expected, the amplitude analysis does not allow us to reduce the parasitic effects due to the pictorial layer.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Amplitude images of the calculated harmonic response.</figcaption></figure>

Figure 6 shows the phase images of the calculated harmonic response at the same frequencies. The 6 defects located in the analyzed sample are well detected. This phase analysis shows that it has been possible to homogenize the background, even in bands with different optical absorptivity. This result is encouraging, as it seems to reduce the parasitic optical effects induced by the different colors of the pictorial layer. It is now necessary to confirm this experimentally.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Phase images of the calculated harmonic response.</figcaption></figure>

# 4. Experimental Study

## 4.1. *Samples Description*

Based on the results obtained by the theoretical study, an experimental validation of the PPT should be implemented in order to ensure the demonstration. Two samples have been built with different internal defects.
The first one dedicated to our research was built in our laboratory, and is a block of plaster whose thermal properties are close to those of a wall paint. It is covered with a pictorial layer (acrylic paint) composed of 20 rectangular areas of various colors covering the visible spectrum (from purple to red), moreover this choice of pigments allows representing the different cases of figure in terms of absorptivities: The black color has a high absorptivity, the white color is highly diffuse, the silver color is very reflective, and so on. These different colors will have very different behaviors in terms of energy absorption.

The geometric dimensions of this block are 68 cm wide, 50 cm high and 5 cm thick in order to be able to study in one recorded area with a sufficient resolution.

In order to simulate delaminations, 32 extruded polystyrene inserts covering the entire surface were introduced into this sample (Figures 7 and 8). These are discs with a diameter of 20 mm and a thickness of 2 mm placed at a depth of 5 mm.

<figure><img src="image_10.png" /><figcaption>Figure 7. The sample analyzed: Front panel.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>Figure 8. The sample analyzed: Rear side.</figcaption></figure>

The second built sample is a partial copy of the mural painting "Saint Maurille en Evêque" of the Angers Cathedral painted in 1270-1280 (Figure 9). This sample was built and painted according to the
painting technics studied before by the research laboratory of historical monuments that discovered one of the first oil paintings of mural painting. It therefore has thermal and optical properties close to the work of art itself.

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Mural painting "Saint Maurille en Evèque".</figcaption></figure>

This copy is shown in Figure 10. It represents the character located to the right of the "Saint Maurille" (white border of Figure 9). We choose this part of the wall painting because it offered different radiative behaviors. The black and dark parts easily absorb the exciter light flux. The lighter parts diffuse it partially. Finally, the gold insert (clover-shaped) behaves like a diffusing mirror. The geometric dimensions of this sample copy are a width and length of 15 cm. Its thickness is equal to 15 mm. To simulate local deterioration of the support stone, 5 extruded polystyrene inserts were introduced into this sample, see Figure 10a. These are 30 mm diameter and 2 mm thick discs. They were placed at variable depths of 2 mm, 3 mm, 4 mm, 5 mm and 6 mm. They were spread as shown in Figure 11.

<figure><img src="image_13.png" /><figcaption>**Figure 10.** The copy of the mural painting "Saint Maurille en Evèque" studied: Front panel (a), back side (b).</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>**Figure 11.** The copy of the mural painting "Saint Maurille en Evèque" studied: Positioning of defects.</figcaption></figure>

## 4.2. *The Experimental Set-Up Implemented for the Study*

Both samples were analyzed using the SAMMTHIR system of the laboratory. This device first includes an excitation by two halogen lamps. It then includes an infrared thermography camera. It is the FLIR SC655 micro-bolometer camera, used for the acquisition of the photothermal signal. It includes an electronic control and data acquisition system, in addition to a data analysis and processing software developed by FLIR ©.

In order to study the first sample, the experimental conditions selected are as follows: the camera is at a distance of about 150 cm from the sample. The 500 W halogen light sources are 45° oriented in the direction of the sample, and they are placed on either side of the camera in order to homogenize the heating flux. The excitation is of the crenel type with a duration time equal to 1 min. The analysis time is equal to 194 s with an acquisition frequency of 1 Hz. Figure 12 shows the first sample during the experiment.

<figure><img src="image_15.png" /><figcaption>**Figure 12.** The sample under analysis.</figcaption></figure>

For the second sample, we have chosen experimental conditions close to the previous ones: the infrared thermography camera type SC655 was placed roughly 50 cm from the work of art. The halogen excitation sources were placed on either side of this camera, 50 cm from the sample and inclined about 45° from its surface with an average power of 2 × 500 Watts. The excitation time was equal to 2 min at 1 Hz of frequency (see Figure 13 of the sample under analysis).
<figure><img src="image_16.png" /><figcaption>**Figure 13.** The copy of the mural painting "Saint Maurille en Evèque" under analysis.</figcaption></figure>

### 4.3. Results and Discussion

The first experiments were performed on the plaster block described above. In Figure 14, we first present the raw thermograms obtained at times t = 10 s, t = 60 s and t = 120 s. They clearly show the possibility of detecting, after 120 s, most of the inserts present in the sample studied. At the same time, the beginning of the recording time, they also show very clearly the important influence of the pictorial layer emissivity on the collected photothermal signal.

<figure><img src="image_17.png" /><figcaption>**Figure 14.** Example of experimental thermograms obtained (at t = 10 s, t = 60 s and t = 120 s).</figcaption></figure>

In Figures 15 and 16, we present the results after a Fourier transform treatment and analysis. Figure 16 shows the obtained module images corresponding respectively to frequencies of 5.15 mHz, 10.30 mHz, and 0.9 Hz. As expected, it shows that these signatures allow a partial detection of the defects, but also that they are very sensitive to inhomogeneities of energy deposition.

<figure><img src="image_18.png" /><figcaption>**Figure 15.** Module images of the calculated harmonic response.</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>**Figure 16.** Phase images of the calculated harmonic response.</figcaption></figure>

Figure 16 shows the obtained phase images at the same frequencies. On one side, we get the possibility of detecting almost all the defects located in the sample studied. On the other side, they show a significant rejection of the optical effects induced by the pictorial layer, which was the desired goal. This method, which is based on the PPT, seems to be really appropriate to reduce the optical effects induced by a pictorial layer.

With this first encouraging experimental result, we were focused on the copy of the "Saint Maurille en Evèque". In Figure 17, we present the raw thermograms obtained at times t = 1 s, t = 41 s, and t = 108 s. On the one hand, they show a particular photo thermal signature at the place of the 5 defects since it allows their detection. On the other hand, they show the important influence of the pictorial layer on the photothermal signal collected. We made the same conclusions previously.

<figure><img src="image_20.png" /><figcaption>**Figure 17.** Example of experimental thermograms obtained (at t = 1 s, t = 41 s and t = 108 s).</figcaption></figure>

The results obtained after post-treatment are shown in Figures 18 and 19. Figure 19 shows the obtained module pictures of the photothermal signal. They correspond respectively to 16 mHz, 24 mHz, and 32 mHz frequencies. It shows a better detection of defects, but also a sensitivity to the patterns of the pictorial layer.

Figure 19 shows the obtained phase images at the same frequencies. On the one hand, they show the possibility of detecting all the defects located in the sample studied. On the other hand, they show a significant rejection of the optical effects induced by the pictorial layer, which was the desired objective. This confirms on the first part the previous theoretical and experimental results. This demonstrates the high interest of a PPT-type analysis in assisting in the conservation of wall paintings.
<figure><img src="image_21.png" /><figcaption>**Figure 18.** Amplitude images of the calculated harmonic response (at f = 16m Hz, f = 24 mHz and f = 32 mHz).</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>**Figure 19.** Phase images of the calculated harmonic response (at f = 16 mHz, f = 24 mHz and f = 32 mHz).</figcaption></figure>

**5. Conclusions**

In this work, we approached the contribution of a Pulse Phase Thermography type to improve the application during the restoration and conservation process of works of art, and more specifically in mural painting applications.

First, we underlined that stimulated infrared thermography was already very effective for assistance in the restoration and conservation, allowing the detection of delamination or chemical contamination located in murals paintings (salts for example). The difficulty is that the effect of the pictorial layer does not absorb the excitatory flux in homogeneous ways. This difficulty induces the detection of artifacts signal that conduct to a false identification of defects. This ghost defect should be removed in order to ensure the delivery of the right state of conservation.

To solve this problem, we have suggested studying the contribution of a Pulse Phase Thermography (PPT) analysis.

Firstly, a simulation approach of the photothermal experiment has been developed, using finite element methods to solve the heat equation. We then implemented an experimental study of two representative samples of heritage mural paintings to demonstrate the efficiency of the approach.

In all three cases, we studied the raw thermograms and then the amplitude and phase maps given by a Fourier Transform. These cases, the raw thermograms and the amplitude maps, allow a good detection of localized failures in the samples studied. However, they remain very sensitive to the different surface radiative properties.

On the other terms, phase maps show, both theoretically and experimentally, a good detection of defects, but also a very good rejection of these optical problems.
These results are very encouraging. They seem to improve the non-destructive testing of cultural
heritage using stimulated infrared thermography. Further studies in order to generalize these results
and to confirm them must be carried by in situ analyses.

**Author Contributions:** Writing—original draft, K.M.; software, K.M., J.-L.B.; writing—review and editing, K.M., V.D., J.-M.V. and J.-L.B.

**Funding:** This research received no external funding.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Ambrosini, D.; Daffara, C.; Di Biase, R.; Paoletti, D.; Pezzati, L.; Bellucci, R.; Bettini, F. Integrated reflectography and thermography for wooden paintings diagnostics. *J. Cult. Herit.* **2010**, *11*, 196–204. [CrossRef]

2. Avdelidis, N.; Moropoulou, A. Applications of infrared thermography for the investigation of historic structures. *J. Cult. Herit.* **2004**, *5*, 119–127. [CrossRef]

3. Barreira, E.; de Freitas, V.P. Evaluation of building materials using infrared thermography. *Constr. Build. Mater.* **2007**, *21*, 218–224. [CrossRef]

4. Bison, P.; Bortolin, A.; Cadelano, G.; Ferrarini, G.; Finesso, L.; Lopez, F.; Maldague, X. Evaluation of frescoes detachments by partial least square thermography. In Proceedings of the 12th International Conference on Quantitative Infrared Thermography, Bordeaux, France, 7–11 July 2014.

5. Bison, P.; Bortolin, A.; Cadelano, G.; Ferrarini, G.; López, F.; Maldague, X. Comparison of image processing techniques for the on-site evaluation of damaged frescoes. In Proceedings of the SPIE Sensing Technology+Applications, Baltimore, MD, USA, 21 May 2014; p. 91050E.

6. Blessley, K.; Young, C.; Nunn, J.; Coddington, J.; Shepard, S. The feasibility of flash thermography for the examination and conservation of works of art. *Stud. Conserv.* **2010**, *55*, 107–120. [CrossRef]

7. Bodnar, J.L.; Nicolas, J.L.; Mouhoubi, K.; Candore, J.C.; Detalle, V. Characterization of an Inclusion of Plastazote Located in an Academic Fresco by Photothermal Thermography. *Int. J. Thermophys.* **2013**, *34*, 1633–1637. [CrossRef]

8. Candore, J.C.; Bodnar, J.L.; Detalle, V.; Grossel, P. Characterization of defects situated in a fresco by stimulated infrared thermography. *Eur. Phys. J. Appl. Phys.* **2012**, *57*, 11002. [CrossRef]

9. Carlomagno, G.M.; Meola, C. Comparison between thermographic techniques for frescoes NDT. *NDT E Int.* **2002**, *35*, 559–565. [CrossRef]

10. Chambard, J.P.; Roche, A. Defect detection of wall paintings in the Chateau de Versailles using TV-holography and IR thermography. In Proceedings of the Optical Metrology, Munich, Germany, 16 July 2007; p. 66180Y.

11. Daffara, C.; Pezzati, L.; Ambrosini, D.; Paoletti, D.; Di Biase, R.; Mariotti, P.I.; Frosinini, C. Wide-band IR imaging in the NIR-MIR-FIR regions for in situ analysis of frescoes. In Proceedings of the SPIE Optical Metrology, Munich, Germany, 20 June 2011.

12. Bruder, R.; Detalle, V.; Coupry, C. An example of the complementarity of laser-induced breakdown spectroscopy and Raman microscopy for wall painting pigments analysis. *J. Raman Spectrosc.* **2007**, *38*, 909–915. [CrossRef]

13. Anglos, D.; Detalle, V. Cultural Heritage Applications of LIBS. In *Laser-Induced Breakdown Spectroscopy*; Springer: Berlin/Heidelberg, Germany, 2014; pp. 531–554.

14. Tornari, V.; Bonarou, A.; Castellini, P.; Esposito, E.; Osten, W.; Kalms, M.K.; Stasinopulos, S. Laser-based systems for the structural diagnostic of artwork: An application to XVII-century Byzantine icons. In Proceedings of the Laser Techniques and Systems in Art Conservation, Munich, Germany, 22 October 2001; Volume 4402, pp. 172–184.

15. Blumich, B.; Casanova, F.; Perlo, J.; Presciutti, F.; Anselmi, C.; Doherty, B. Noninvasive testing of art and cultural heritage by mobile NMR. *Acc. Chem. Res.* **2010**, *43*, 761–770. [CrossRef] [PubMed]

16. Maldague, X.; Marinetti, S. Pulse phase infrared thermography. *J. Appl. Phys.* **1996**, *79*, 2694–2698. [CrossRef]

17. Maldague, X.; Galmiche, F.; Ziadi, A. Advances in pulsed phase thermography. *Infrared Phys. Technol.* **2002**, *43*, 175–181. [CrossRef]

18. Ibarra-Castanedo, C.; Maldague, X. Pulsed phase thermography reviewed. *Quant. Infrared Thermogr.* **J. 2004**, *1*, 47–70. [CrossRef]
19. Briggs, W.L. *The DFT: An Owners' Manual for the Discrete Fourier Transform*; Siam: Philadelphia, PA, USA, 1995; Volume 45.

20. Rao, K.R.; Yip, P. *Discrete Cosine Transform: Algorithms, Advantages, Applications*; Academic Press: Cambridge, MA, USA, 2014.

<figure><img src="image_23.png" /></figure>

© 2019 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 (http://creativecommons.org/licenses/by/4.0/).