1. Introduction
Speckle Plethysmography (SPG) and Photoplethysmography (PPG) are fundamentally different biophotonics methods, both based on the measurement of photons to analyse blood dynamics in living tissue. SPG is a technology based on the analysis of Laser Speckle Contrast Imaging (LSCI) [
1,
2,
3,
4,
5,
6]. Spatial changes in speckle contrast are calculated from a sequence of video frames, generating a single value per frame. Changes in speckle contrast are produced by variations, such as movement of the target under analysis (in this case, human tissue cells: skin, blood vessels, blood cells, etc.).
PPG is a technology based on the measurement of light intensity. This technology is used to estimate modulations in the absorption of light by a living tissue, which relates to blood volume [
7]. When cameras are used, intensity changes are calculated throughout a sequence of video frames producing an average value per frame.
Simultaneous monitoring of SPG and PPG could be used to study changes induced by blood circulation modulations [
8]. A well-known technique used to systemically manipulate the peripheral circulation is the Cold Pressor Test (CPT) [
9], which is performed through the application of cold temperatures to a part of the body. To preserve core temperature, sensory nerves trigger a systemic reaction, known as cold-induced vasoconstriction, leading to a decrease in the vessel diameter and increase in tone of the vessel walls. This results in elevated blood pressure (BP), an increase in Pulse Wave Velocity (PWV), and reduced skin blood flow [
10]. Pulse Pressure (PP), which is calculated as the difference between systolic and diastolic BPs, has also been proven to increase [
11]. Stroke volume decreases, which in turn decreases PPG amplitude [
11]. For PPG, this decrease in stroke volume deteriorates the Signal-to-Noise Ratio (SNR) and pulse strength, which sometimes even results in a complete loss of the signal [
12]. Different parts of the body react differently to the CPT. While core body parts do not seem affected, peripheral parts typically react with a decrease in the AC-amplitude of PPG. This has been demonstrated for ear and finger measurements [
13].
The effects of the CPT on peripheral Near-Infrared SPG have been explored by Ghijsen et al. [
1], whose findings show an increase in the time delay between SPG and PPG. SNR changes were also observed but not extensively analysed. Nevertheless, these observations suggested that SPG might be more robust to changes under low perfusion conditions compared to PPG. Therefore, SPG might become a promising alternative to monitor peripheral circulation under low perfusion conditions [
11,
12], although it deserves further investigation.
To our knowledge, the SNR and morphology of dual-wavelength SPG and PPG during the CPT have never been investigated before. This investigation focusses on analysing the changes induced by a lowered perfusion condition on SPG and PPG. The experiments were performed using a custom setup, which can measure simultaneous SPG and PPG from skin at two wavelengths (639 nm and 850 nm) multiplexed in time. Dual-wavelength SPG and PPG were measured at the right index finger location before and during the CPT.
The change produced by the CPT in AC-amplitude of dual-wavelength SPG and PPG signals across 10 subjects was evaluated. In addition, changes in SNR due to the CPT were calculated for SPG and PPG at two wavelengths. Furthermore, using frequency harmonic ratios, SPG and PPG waveform morphology characteristics were analysed for all subjects and compared to the finger Arterial Pressure (fiAP) reference.
2. Materials and Methods
The custom setup used for this investigation incorporated a miniature CMOS camera, and two laser diodes with wavelengths of 639 nm (Thorlabs HL6358MG, Thorlabs, Newton, NJ, USA, 10 mW) and 850 nm (Thorlabs L850P010, 10 mW), embedded in a finger-clip. These wavelengths were chosen because of their difference in penetration depth within human tissue. This, in turn, is caused by the different absorption and scattering coefficient of tissue components at these wavelengths. Especially, the spectral absorption differences between reduced and oxidized haemoglobin have been thoroughly exploited for oxygen saturation measurement (pulse oximetry, based on PPG). Even if absorption is not expected to have a significant effect on SPG performance [
2], the SPG waveform still may be affected by the difference in penetration depth. The setup was configured to allow for the simultaneous measurement of reflective SPG and PPG multiplexing in time between the two wavelengths. The raw camera video stream was recorded at a framerate of 80 Frames Per Second (FPS), while the 639 nm and 850 nm light sources were alternating. We then separated the even and uneven frames into two video streams of 40 FPS each with only one wavelength. The distance between the camera and the finger skin was approximately 1 cm. The subjects positioned their right index finger on the setup.
A Finapres Nova (Finapres Medical Systems, Enschede, The Netherlands) was used as a reference device to monitor fiAP, based on the Peñáz method [
14]. The custom experimental system was synchronized with the reference device using the acquisition platform BIOPAC MP160 (BIOPAC Systems Inc., Goleta, CA, USA). The Finapres finger cuff was applied on the right middle finger. This device was used as reference for heart rate and for harmonic ratios comparison.
The experimental study was performed on 10 subjects without known medical conditions. Both the custom setup and the Finapres recorded continuously during the protocol. The study consisted of two experimental phases: baseline and CPT. During the baseline phase, which lasted 1 min, the subject was asked to sit relatively still. Whereas, during the CPT, the subject was asked to submerge the left hand in ice water (total hand immersion) for 1 min while the reference and custom setup were recording the signals on the right-hand fingers. The duration of 1 min was adapted following the previous literature [
15,
16]. No metabolic activities were considered for this study. The ice water was prepared by mixing 3 litres of cold tap water with 0.5 kg of ice cubes, to be used after waiting for 5 min at room temperature. The experiment took place in a laboratory with an ambient temperature of 20 (+/−2) degrees Celsius.
The subject population consisted of 10 healthy participants ranging from 28 to 53 years old; all subjects were Caucasian, with a 40 percent female population. The experimental protocol was evaluated and executed in accordance with the Declaration of Helsinki 2008 (protocol number: IM-NL-STUDY-2022-0029), and all subjects signed the informed consent prior to the study.
Figure 1 summarizes the experimental protocol and data analysis pipeline.
2.1. Analysis
All the processing and analysis was performed using Python programming language (Python Software Foundation 3.8). Raw videos were recorded at 80 FPS of the index finger, while using synchronized illumination at two alternating wavelengths (40 FPS each). After separating the alternating frames by wavelength, SPG and PPG were extracted for each wavelength at 40 samples per second. The equations used to compute SPG and PPG [
17] are shown below. Let
be a frame of one of the two video streams at time
, with resolution
, corresponding to
pixels. The
PPG was extracted by calculating the average pixel intensity of the frame as
where
and
represent the two indices running through the two frames’ dimensions.
SPG was calculated from the same video as the average spatial variability of each frame as
with
being the spatial variability matrix calculated in the neighbourhood of the pixel at position
. Following the well-known standard deviation definition, the spatial variability can be expressed as
with
being the kernel size, here equal to
pixels based on the literature [
18]. An example of the extracted signals from the videos, both SPG and PPG at the two wavelengths, is shown in
Figure 2.
2.1.1. AC-Amplitude and SNR Calculations
PPG and SPG, obtained from the videos at both wavelengths, were further analysed to monitor the changes in AC and SNR between the baseline and CPT phases. The reference signal, fiAP, was included in the AC analysis; the AC of a BP signal can also be indicated as PP. The reference was not included in the SNR analysis since the signal is not raw but prefiltered by the Finapres system. However, fiAP was used for the estimation of the fundamental frequency in the SNR estimation.
Both AC-amplitude and SNR calculations were performed using a sliding window approach, with a window size of 10 s and an overlap of 9 s. The window size was chosen taking into account the trade-off between time and frequency resolutions. Both phases were included in this calculation. The two phases were separated, and 5 s signals were removed at the start and end of the baseline and CPT phases in order to avoid motion artifacts that could be present, e.g., due to the subject submerging the left hand in the ice water.
Figure 3 summarizes the processing steps for the calculation of AC-amplitude and SNR.
For the AC-amplitude calculations, each window of the signals of interest (SPG, PPG, and fiAP) in both phases and wavelengths was first filtered using a bandpass filter (Butterworth cut-off frequencies: 0.58 Hz–4.17 Hz). Afterwards, peaks and valleys were located, and an interpolation was performed to obtain a line connecting all the peaks and a line connecting all the valleys, similar to an upper and lower envelope. The AC-amplitude of the signal in the current window was then calculated as the median value of the difference between the upper and lower envelopes.
For all subjects, the SNR was analysed by comparing the spectral power contained between certain frequency bands considered as signal and frequency bands of noise. This is carried out using fiAP as reference for the fundamental frequency peak, similar to the method described by Gerard de Haan et al. [
17]. A binary mask was created to split the spectrum between signal and noise. This binary mask was equal to “1” in the frequency bands containing the fundamental frequency peak, calculated from fiAP reference, and in the 4 following harmonics of SPG and PPG. The rest of the binary mask was equal to “0”, representing frequencies corresponding to noise. The spectrum, limited to 800 beats per minute (BPM) was calculated, using zero-padding and a Hanning window, on the high passed filtered signals with a cut-off frequency of 0.58 Hz, in order to remove possible baseline and respiratory contributions. The width of the fundamental band was empirically set to 20 frequency bins and proportionally increased for the harmonic bands [
17]. An example of this procedure is visible in
Figure 4. More in details, let
be the spectrum of the current 10 s window of the signal of interest (SPG, PPG), and
the described binary mask; the
SNR can be calculated as
2.1.2. Harmonic Ratio Calculation
The harmonic content of the signals was analysed to quantify the morphology differences between SPG, PPG, and fiAP. This analysis was based on the ratio between the different frequency harmonics of the signal and the fundamental, inspired by Ghijsen et al. who defined the Third Harmonic Ratio (THR) [
1]. In this work, Second Harmonic Ratio (SHR) and Fourth Harmonic Ratio (FHR) were introduced and included in the analysis together with THR. The harmonics of a signal are influenced by its morphology, and the use of the three harmonic ratios allows for better comparison of the morphological differences and similarities between PPG, SPG, and fiAP. In all three cases, harmonic ratios are calculated as the ratio between the amplitude of the pertaining harmonic peak of interest and the fundamental peak.
The signals chosen for this analysis were produced by the 850 nm laser diode during the baseline test phase since they show the highest quality across subjects. The harmonic content was analysed on a representative beat per signal. This representative beat was selected by first calculating an average beat of each signal per subject and then selecting the real pulse that correlated the most with the average beats. A summary of the steps can be found in
Figure 5.
The representative beats are concatenated 10 times with the purpose of obtaining a Fourier frequency spectrum with reasonable resolution. This results in a frequency spectrum for individual PPG, SPG, and fiAP beats.
Figure 6 shows the representative beat of the collected signals for all the subjects during the baseline phase. It should be noted that the representative beats may not be the same, timewise, for all the three signals because the representative beats are selected based on the correlation coefficient.
2.1.3. Statistical Analysis
For each subject, phase, signal, and wavelength a median value is calculated for both AC and SNR. These values are then analysed with paired comparisons. Four statistical comparisons were performed to analyse the significance of the changes. Two main statistical tests were used, assuming the populations have significant deviations from normality. Firstly, the Wilcoxon signed-rank test was chosen as a non-parametric version of the paired t-test, in order to perform paired comparisons. Secondly, the Levene’s test was chosen to compare variances. Four comparisons were made:
Comparison between phases of the study (baseline and CPT): the Wilcoxon signed-rank test was used inter-subject to analyse if the changes in both AC-amplitude and SNR between the baseline and CPT phases were significant for the signals of interest. This comparison was carried out since it is expected for PPG to have a reduction in AC-amplitude during the CPT and consequently SNR [
15], but it is not clear whether SPG will react in the same way.
Comparison between experimental signals (SPG and PPG): the Wilcoxon signed-rank test was used to test if SNR SPG was significantly higher than SNR PPG for both study phases and wavelengths. This comparison was performed because it has been suggested that SPG has a higher SNR than PPG during the CPT [
1].
Comparison between SNR variance of experimental signals (SPG and PPG): the Levene’s test was used inter-subject to analyse if the SNR variance was significantly different between SPG and PPG. For this test all the SNR estimations were used as populations, and not the median, as for the other tests. This comparison was performed because we expected the SNR of SPG to be more consistent across subjects compared to PPG.
Comparison between harmonic ratios of the signals (fiAP, SPG, and PPG): the Wilcoxon signed-rank test was used inter-subject to test if SHR, THR, and FHR were significantly higher between the different pairs of signals. These comparisons were chosen to verify whether the harmonic content of SPG is higher compared to PPG, and more similar to fiAP.
In addition, the Pearson’s correlation coefficient and the regression line were estimated for the three harmonic ratios to compare fiAP morphology with SPG and PPG.
4. Discussion
In this work, reflective PPG and SPG were obtained at two wavelengths (639 and 850 nm, using time-multiplexing) from the right index finger of 10 healthy volunteers, using fiAP as a reference under lowered perfusion conditions by applying a CPT. The CPT consisted of submerging the subject’s hand in ice water for 1 min. This commonly causes cold-induced vasoconstriction, leading to lowered perfusion conditions also in the non-immersed hand where the measurements took place. The data collection consisted of a baseline phase and a CPT phase. Per wavelength, PPG and SPG were calculated from the same video stream using different processing pipelines.
As seen in the Results Section, both AC-amplitude and SNR of PPG and SPG decrease during CPT-induced lowered perfusion conditions. This only reached statistical significance for the signals at a wavelength of 850 nm and not for 639 nm. Furthermore, SPG has consistently higher SNR values than PPG, with lower intra- and inter-subject variation. In addition, the harmonic analyses revealed that there is a larger resemblance between SPG and fiAP compared to PPG and fiAP.
In healthy subjects, the CPT leads to an increase in BP and a decrease in blood flow and volume. It is expected that PPG will have a reduced AC-amplitude as the blood volume decreases. Since SPG is a speckle-based method, it might be influenced by blood flow [
1] as well as by blood pressure [
8]. From the AC-amplitude analysis, it is shown that the change due to the CPT was significant with 850 nm illumination (see
Figure 9). Both SPG and PPG signals at 850 nm showed a significant decrease, while both signals at 639 nm did not show a significant decrease. This could suggest that the signals at the two wavelengths partly contain different information. A possible explanation might be the difference in penetration depth or changes in chromophores [
19]. However, the non-significance of the AC-amplitude decrease in PPG at 639 nm could also be caused by the lower quality, and hence variation, of the signal overall. fiAP, on the other hand did not show a significant reduction as expected. The PP showed an increase during the CPT, compared to the baseline, of 9% ± 10%.
SPG resulted in having a significantly higher SNR compared to PPG in both the baseline and CPT, and both wavelengths, as shown in
Figure 11 and
Figure 12. During a CPT, the reduction in volumetric flow in the peripheral circulation and the shunting of the capillaries influences the SNR of both PPG and SPG. These lower perfusion conditions depict a challenging situation for PPG measurements. However, the SNR of SPG is significantly higher than PPG during normal conditions but also during the CPT. Therefore, in poor perfusion cases where PPG could be lost, SPG might still preserve enough signal quality for heart rate (HR) measurements. In addition, the variance in the SNR was found to be significantly different between SPG and PPG for all phases and light sources applied, with one exception. This exception was PPG at 639 nm during the CPT, where the SNR was already low before the CPT and became consistently low across subjects during the CPT, reaching a consistency in SNR levels similar to the SPG. This explains how the SNR across subjects during the two phases of the measurement is significantly more consistent for SPG.
The change in SNR was further analysed for each signal between the two phases, as shown in
Figure 13. Similar to the AC-amplitude changes, for 850 nm illumination, there was a statistically significant difference in the SNRs between the baseline and CPT phases. This means that both PPG and SPG have a statistically significant reduction in the SNR during the CPT. This reduction is less present at 639 nm, however, it is worth mentioning that the SNR for PPG at 639 nm was already very poor during baseline measurements, as shown in
Figure 11. This might explain the fact that the decrease in the SNR is not significant.
These results indicate that SPG could be more accurate than PPG for HR measurements for most situations, and especially in conditions with low perfusion. Furthermore, this could even extend to more accurate Heart Rate Variability (HRV) and Pulse Arrival Time (PAT) measurements over long recording periods [
8,
20], which could be applied for more accurate non-invasive blood pressure assessment.
In all subjects, the harmonic ratios measured during baseline on the 850 nm illumination source were found to be significantly different between SPG and PPG, as visible in
Figure 15. fiAP shows the most pronounced harmonic ratios, closely followed by SPG, and PPG shows the lowest harmonic ratios. This means that the secondary peak is usually more pronounced for fiAP and SPG than for PPG. SHR was not significantly different between fiAP and SPG. The correlation between fiAP and SPG is higher than PPG for all harmonic ratios (see
Figure 16), which suggests a bigger relation of SPG with BP than PPG. It could be speculated that differences between harmonic ratios derived from simultaneous SPG and PPG signals might provide information about the state of the circulatory system, but this requires more investigation. However, considering the differences in SPG morphology depend on age, as already shown in [
1], this analysis should be repeated on a larger population. SPG has potential to be used as a method to assess cuffless blood pressure, but care should be taken to directly use AC values, since they are still subjected to perfusion changes independent of BP.
This is the first investigation analysing dual-wavelength SPG and PPG from a reflective contact probe. Measuring in reflective mode is suitable for more body locations than transmission mode; it allows for more compact setups and offers the possibility of integration in wearable devices. In addition, the custom setup designed for this work adds extra information due to the dual-wavelength usage. However, the relatively small population, 10 subjects, forms a limitation of the study. This is especially relevant if the differences in CPT reactions between different ages and sexes are considered [
21].
5. Conclusions
The CPT has an effect on the SNR of both SPG and PPG signals; at 850 nm, most of the participants show a decrease in the SNR under the CPT. However, the SNR of SPG was significantly higher than the PPG one, with significantly lower variability. This suggests that SPG measurements could be more suitable for low perfusion situations. The AC-amplitude of SPG and PPG signals at 850 nm show a significant reduction during the CPT, but this is not seen for 639 nm. This suggests that the two wavelengths might contain different information and might be affected differently by haemodynamics.
The harmonic ratio analysis indicates that fiAP, PPG and SPG are different signals, containing different information about the state of the cardiovascular system and blood dynamics. Harmonic ratio analysis revealed clear differences between SPG and PPG. Furthermore, the correlation between fiAP and SPG is significantly higher than the one between fiAP and PPG.
Thus, to provide more information on the monitoring of blood dynamics, combining SPG and PPG measurements should be considered, instead of using only PPG. Furthermore, while low perfusion conditions that can compromise the pulse amplitude are challenging to monitor with PPG, SPG might offer better results. Some examples of such conditions are Raynaud’s syndrome, diabetes, vasculitis, and advanced age. Studies on a larger population, including healthy subjects and patients with such conditions, should be conducted to further conclude on these points.