Next Article in Journal
New-Type Shearing Interference Detection System Based on Double-Grating Structure for Suppressing the Skylight Background
Next Article in Special Issue
Low-Power Long-Term Ambulatory Electrocardiography Monitor of Three Leads with Beat-to-Beat Heart Rate Measurement in Real Time
Previous Article in Journal
Efficient Deep Learning Based Hybrid Model to Detect Obstructive Sleep Apnea
Previous Article in Special Issue
Deriving Multiple-Layer Information from a Motion-Sensing Mattress for Precision Care
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

ECG-Free Heartbeat Detection in Seismocardiography Signals via Template Matching

Department of Electrical Engineering and Information Technologies, University of Naples Federico II, Via Claudio 21, 80125 Naples, Italy
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(10), 4684; https://doi.org/10.3390/s23104684
Submission received: 26 April 2023 / Revised: 8 May 2023 / Accepted: 10 May 2023 / Published: 12 May 2023
(This article belongs to the Special Issue Human Signal Processing Based on Wearable Non-invasive Device)

Abstract

:
Cardiac monitoring can be performed by means of an accelerometer attached to a subject’s chest, which produces the Seismocardiography (SCG) signal. Detection of SCG heartbeats is commonly carried out by taking advantage of a simultaneous electrocardiogram (ECG). SCG-based long-term monitoring would certainly be less obtrusive and easier to implement without an ECG. Few studies have addressed this issue using a variety of complex approaches. This study proposes a novel approach to ECG-free heartbeat detection in SCG signals via template matching, based on normalized cross-correlation as heartbeats similarity measure. The algorithm was tested on the SCG signals acquired from 77 patients with valvular heart diseases, available from a public database. The performance of the proposed approach was assessed in terms of sensitivity and positive predictive value (PPV) of the heartbeat detection and accuracy of inter-beat intervals measurement. Sensitivity and PPV of 96% and 97%, respectively, were obtained by considering templates that included both systolic and diastolic complexes. Regression, correlation, and Bland–Altman analyses carried out on inter-beat intervals reported slope and intercept of 0.997 and 2.8 ms (R2 > 0.999), as well as non-significant bias and limits of agreement of ±7.8 ms. The results are comparable or superior to those achieved by far more complex algorithms, also based on artificial intelligence. The low computational burden of the proposed approach makes it suitable for direct implementation in wearable devices.

1. Introduction

Continuous heart rate monitoring is fundamental in critical medical circumstances to assess cardiac function. It is a vital task in intensive care units, surgery, post-anesthesia surveillance, and emergency medicine [1,2,3]. Heart rate monitors are extensively used for analysis of heart rate variability (HRV), which is considered a strong predictor of the outcome in patients with high cardiovascular risk [4,5]. Indeed, people with coronary artery disease, heart failure, valvular heart disease, or after a myocardial infarction, ischemia, or a stroke need to be continuously monitored to provide them with timely therapeutic intervention, thus preventing the deterioration of health conditions and the occurrence of fatal arrhythmic events [1,4,6]. HRV analysis has also been recognized as a valuable non-invasive tool for identifying cardiac autonomic dysfunction in neurological (e.g., multiple sclerosis), renal (e.g., end-stage renal failure), or metabolic (e.g., diabetes, obesity) diseases [6,7,8,9,10]. Moreover, heart rate and HRV indices can be useful: in sleep medicine to assess sleep quality and detect sleep apnea syndrome, or neonatal distress syndrome [1,6,8,11]; in fitness and sport science to monitor exercise training and the performance of athletes [8,12,13]; or even in rehabilitation medicine to track the exercise of chronic patients [1,8]. Furthermore, they can provide valuable information: to evaluate the effect of drugs [8]; to monitor stress levels and obtain clear evidence of emotion disorders (e.g., anxiety, depression, etc.) [1,14]; to promote lifestyle awareness in people and improve their well-being [1,8]; and to enhance road safety by detecting driver fatigue [15].
Electrocardiography (ECG) records the electrical activity produced by the heart on the body’s surface and is currently regarded as the gold standard for heart rate monitoring in clinical routine. Heart rate measurements rely on heartbeat detection, which is usually accomplished by locating the R-wave of each QRS complex in ECG signals [1,16,17,18,19,20,21]. The well-known Pan and Tompkins algorithm is the most popular approach for this purpose [1,22]. However, the ECG has some drawbacks: it is intrinsically affected by electromagnetic interference, needs stable placement of the electrodes over time (electrode slipping and detachment may happen) and good skin–electrode impedance (which worsens with the drying of the electrolytic gel), involves electrical risks. In addition, it must be performed by skilled clinical professionals and is somewhat uncomfortable. Therefore, this technique is not well suited to continuous, long-term monitoring, especially in daily life environments, although some portable instruments, such as the Holter device, have been proposed as possible solutions for home monitoring (the Holter monitor generally allows recordings of no longer than 72 h). Furthermore, the ECG does not provide any information about the cardiac mechanical function [1,16,17,18,19,20,21].
In addition to the ECG, heart rate monitoring can also be pursued by detecting the peripheral effect of cardiac contractions through the Photoplethysmography (PPG) technique. The PPG signal is a non-invasive recording of blood volume changes in a macrovascular bed of the skin due to arterial pulsations. Particularly, blood volume variations are monitored indirectly by measuring changes in the light absorption or transmission of human tissues at characteristic wavelengths. Generally, beat-by-beat heart rate measurement via PPG involves identifying the temporal locations of specific markers on the signal waveform, such as the systolic peak, which are related to important cardiac cycle events. PPG sensors are also integrated into wearable devices (e.g., smartwatches). However, despite the simplicity, low-cost, and unobtrusiveness of this technique, PPG signal quality is very sensitive to motion artifacts, which impair accurate systolic peak detection. Data corrupted by motion artifacts are often unusable, thus jeopardizing continuous heart rate monitoring [1,23,24,25,26,27,28].
Since the second half of the 19th century, a variety of non-invasive, cardio-mechanical monitoring tools have been proposed to record the small vibrations generated by the mechanical activity of the beating heart and to obtain a more comprehensive evaluation of cardiac function. Among these, Ballistocardiography (BCG) [29,30,31,32,33,34], Phonocardiography (PCG) [35,36,37,38], and Seismocardiography (SCG) [39,40,41,42,43,44,45,46] are the subjects of current research. On one hand, BCG records whole-body vibrations due to the blood flowing through the vascular system by means of different instruments, such as weighing scales, systems embedded in a bed or chair, and even wearable sensors, which are particularly appealing for continuous monitoring applications. On the other hand, PCG captures the sonic components of cardiac-induced mechanical vibrations of the thorax, commonly known as heart sounds, by using modern electronic stethoscopes, while SCG measures infrasonic precordial accelerations. Specifically, the availability of lightweight, miniaturized, low-cost accelerometers, manufactured via micro-electromechanical systems (MEMS) technologies, has led SCG to gain particular attention, especially for the development of wearable devices. In recent years, novel, non-invasive cardio-mechanical monitoring techniques have been introduced, namely Gyrocardiography (GCG) [47,48], Kinocardiography (KCG) [49,50], and Forcecardiography (FCG) [51,52,53,54,55,56,57]. GCG records the three-dimensional angular velocities of the precordium via gyroscopes, which are often integrated into inertial measurement units (IMUs) that also contain accelerometers. For this reason, the combined use of GCG and SCG has also been investigated [58]. KCG, instead, results from the simultaneous acquisition of BCG and SCG signals. Finally, FCG measures the local, cardiac-induced mechanical vibrations of the chest wall via piezoresistive or piezoelectric force sensors. Particularly, FCG sensors are characterized by a very broad bandwidth, which allows them to monitor respiration, ventricular volume variations, an SCG-like component, and heart sounds, all simultaneously from a single contact point on the chest.
Measurement of the instantaneous heart rate by cardio-mechanical monitoring techniques is based on beat-by-beat localization of important cardiac cycle events (e.g., the heart valves opening and closing, blood ejection, isovolumic contraction, cardiac filling, etc.), which are marked by specific peaks and valleys. This allows for the estimation not only of the inter-beat intervals (or the instantaneous heart rate), but also further cardiac time intervals of clinical relevance, such as the pre-ejection period, the left ventricular ejection time, and the total systolic time [44,45,54,58,59]. Heartbeat localization incardio-mechanical signals is generally performed by assuming the ECG as a temporal reference, particularly the R-wave in each cardiac cycle. However, research is currently proceeding towards the development of automated, standalone (i.e., ECG-independent) methods, most of which rely on SCG signals. SCG-based heartbeat detectors exploit different approaches, such as signal envelope extraction, often combined with a thresholding operation [60,61,62,63,64,65,66,67,68,69]; continuous or discrete wavelet transform [70,71]; variational mode decomposition [72]; autocorrelated differential algorithm [58]; matched filtering [73]; probabilistic methods, e.g., the hidden Markov model [74]; machine learning [75]; and deep learning [76,77]. Of these, the heartbeat detection methods reported in [58,63] were performed on SCG acquisitions combined with simultaneous GCG recordings. It is important to note that the heartbeat detectors mentioned above, especially those based on artificial intelligence algorithms, have a rather high computational complexity. It is also worth noting that they were tested exclusively on signals from healthy subjects, except for [63], in which 12 patients with coronary artery disease were also considered. This aspect is important because cardiac diseases cause dysfunctions in the mechanical behavior of the beating heart, which result in morphological changes and/or instabilities in the cardio-mechanical signals that may impair the performance of heartbeat detection algorithms. Indeed, in [78], the authors analyzed the waveforms of SCG signals acquired from 90 pathological subjects with myocardial infarction, heart failure, or a transplanted heart, and found that in 38% of subjects the SCG waveform and fiducial points were altered with respect to the traditional references reported in the literature for healthy subjects.
Recently, some authors have presented a template matching approach for ECG-free heartbeat detection in the FCG signal, along with preliminary results on a small cohort of healthy subjects [79]. In particular, it has been applied to the high-frequency component of the FCG signal, which has been shown to share a very high similarity with the SCG signal [51,52,54].
This study was aimed at verifying the feasibility of this approach for SCG signals and achieving an extensive estimation of its performance on a much larger cohort of pathological subjects. Template matching based on normalized cross-correlation is a well-known technique for event detection, and its application to SCG signals could provide a simple, yet effective and robust method for ECG-free heartbeat detection. Indeed, the template matching approach has an intrinsic ability to deal with the inter-subject variability in SCG morphology, because it does not rely on any a priori assumption about the signal shape or specific fiducial points to recognize the signal chunks that match the selected template. Methods based on template matching have previously been applied to BCG [80] and PCG [81] signals, but they have been assessed only on a very small cohort of healthy subjects. The proposed template matching approach, instead, was tested on a publicly available database of SCG signals acquired from patients with valvular heart diseases (VHDs). To the best of our knowledge, no heartbeat detection method reported in the literature has been proven on such a large cohort of pathological subjects. The beat-by-beat performance of the proposed algorithm was assessed by comparing the inter-beat interval estimates obtained from the SCG signals with those obtained from simultaneous ECG recordings. The results of the statistical analyses carried out to compare the SCG- and ECG-based measures demonstrate that the proposed approach can provide very accurate measurements of inter-beat intervals, even on a large cohort of patients affected by one or more VHDs.

2. Materials and Methods

2.1. Dataset

In this study, SCG signals from a publicly available database [82] were analyzed. The database includes SCG recordings from 100 patients (59 males and 41 females, age 68 ± 14 years) with VHDs. The data were collected from two medical centers: 70% in China (patient IDs #CP-01 to #CP-70) and 30% in the United States (patient IDs #UP-01 to #UP-30). All subjects were diagnosed with one or more VHD conditions, namely moderate or greater aortic valve stenosis, aortic valve regurgitation, mitral valve stenosis, mitral valve regurgitation, and tricuspid valve regurgitation. The database also contains ECG signals, which were assumed as the ground truth for heartbeat detection.
Simultaneous SCG and ECG recordings were obtained from a three-dimensional MEMS accelerometer and a biopotential circuit, respectively. Both sensors were embedded in the Shimmer 3 ECG module (Shimmer Sensing, Dublin, Ireland). The device was firmly placed onto the thorax of the subjects, who were asked to lay in a supine position. Four ECG electrodes were attached to the skin of the subjects and connected to the device by cables. The signals from patients #CP-01 to #CP-70 and #UP-01 to #UP-21 were acquired at 256 Hz sampling frequency, while those from patients #UP-22 to #UP-30 were sampled at 512 Hz. All the measurements were carried out during quiet breathing.
The SCG signals from the described database were first explored by visual inspection, and only their dorso-ventral components, corresponding to the z-axis acceleration signals, were taken into account. Particularly, 22 SCG recordings characterized by poor signal quality were excluded from the analysis, as they did not allow a clear identification of the heartbeat waveform. A further SCG signal was also excluded since no ECG recordings had been simultaneously acquired. Therefore, a total of 77 SCG and ECG signals were included in this study. Table 1 reports the patient IDs of discarded SCG signals, along with the reasons for their exclusion.

2.2. Signals Pre-Processing

All the processing operations described in this study were carried out in MATLAB® R2018b (MathWorks, Inc., Natick, MA, USA). To improve the temporal resolution, the SCG and ECG signals were first oversampled at 1 kHz via linear interpolation by using the MATLAB® function “interp1”. Afterwards, the high-frequency component of the SCG signal was extracted via a 4th order zero-lag Butterworth band-pass filter with cut-off frequencies of 7 and 30 Hz. This component has been shown to feature short, oscillatory patterns corresponding to heartbeats [51]. The ECG lead II signal, instead, was first band-pass filtered in the 0.5–40 Hz frequency band via a 4th order zero-lag Butterworth filter. Then, a notch filter was used to remove the 50 Hz powerline interference and its higher harmonics. Finally, the R-peaks were located by means of the well-known Pan and Thompkins algorithm [22], implemented in the “BioSigKit” MATLAB® toolbox [83]. An example of the SCG and ECG signals obtained by the pre-processing operations is depicted in Figure 1.

2.3. Template Matching

After pre-processing, a template matching technique, which had already been used in an FCG study [79], was applied to the SCG signals for heartbeat detection. This technique consists of three steps: (1) selection of a heartbeat template in the SCG signal; (2) computation of the normalized cross-correlation (NCC) function between the selected template and the SCG signal chunks; and (3) localization of the NCC peaks.

2.3.1. Template Selection

The template was defined by selecting a single heartbeat in the SCG signal, which was assumed as a morphological reference. For most of the signals, the typical template waveform included both the systolic and diastolic complexes in the cardiac cycle. This was referred to as “case 1”. Figure 2a shows an example of a template comprising both the systolic and diastolic complexes. However, in a number of SCG signals a clear identification of the diastolic complexes was not possible. In this case, referred to as “case 2”, only the systolic complex of the selected heartbeat was eventually included in the template that was employed for heartbeat detection. An example of a template consisting only of the systolic complex is depicted in Figure 2b.
In detail, 50 SCG recordings fell under “case 1”, while the remaining 27 recordings fell under “case 2”. Table 2 reports the template selection criteria adopted in this study, while also showing the template waveforms recommended for both cases.

2.3.2. Normalized Cross-Correlation and Peaks Localization

After the template had been selected, the NCC function between the SCG signal and the template was computed to evaluate their morphological similarity. The NCC function corresponds to the scalar product between two vectors in the Euclidean space divided by the product of their norms, and is defined by the following mathematical expression:
N C C k = n s n μ s k · t n k μ t n s n μ s k 2 · n t n μ t 2 ,
where s is the SCG signal, t is the selected template, μ t is the mean of the template, and μ s k is the mean of the SCG signal over the shifted template interval. The NCC function has local maxima where the SCG signal reaches the highest similarity with the template. As the template captures the heartbeat pattern in the SCG signal, the highest NCC values occur periodically with the heartbeats. Therefore, the temporal locations of the NCC peaks, which correspond to the heartbeats, were identified via the MATLAB® function “findpeaks”. Specifically, a minimum peak prominence of 0.5 was set for the SCG signals characterized by templates of case 1, while a value in the 0.3–0.7 range was considered for the SCG signals with templates of case 2. These values were selected empirically to fit almost all cases, with the aim to reduce as much as possible the need for parameter optimization by the user. Indeed, the same minimum peak prominence of 0.5 chosen for template 1 was used for 50 out of the 77 subjects. In addition, a minimum peak distance of 500 ms was fixed for all signals. Figure 3 shows some examples of the NCC peak detection. In particular, the NCC signal in panel (a) was obtained by considering a case 1 template, which comprised both the systolic and diastolic complexes in the selected cardiac cycle. The NCC signal in panel (b), instead, was obtained by selecting a case 2 template, which consisted only of the systolic complex.

2.3.3. Inter-Beat Intervals Estimation

Once the NCC peaks had been located, the false positives (FP), false negatives (FN), and detection errors (DE) were identified with respect to the actual heartbeats, provided by the reference ECG signals. The false positives were defined as additional NCC peaks within a single cardiac cycle that were wrongly identified as actual heartbeats (see Figure 4a). The false negatives, instead, were considered as actual heartbeats that were missed by the template matching technique (see Figure 4b). The detection errors were referred to as single NCC peaks occurring in a cardiac cycle, but in a temporal location that was not likely to correspond to a matched heartbeat template (see Figure 4c). For this reason, the detection error represented both a false positive, because it marked the presence of a false heartbeat, and a false negative, because no other peaks were found in the same cardiac cycle that marked the presence of the actual heartbeat. Therefore, the number of detection errors contributed both to the numbers of false negatives and false positives. It is noteworthy to underline that all the NCC peaks corresponding to the ECG heartbeats corrupted by motion artifacts were discarded before FP, FN, and DE recognition. Finally, the inter-beat interval estimates were obtained from the ECG signals as the differences between consecutive R-peaks, and from the SCG signals as the differences between consecutive NCC peaks (see Figure 5).

2.4. Statistical Analyses

To evaluate the performance of the heartbeat detection using the proposed approach, the sensitivity and positive predictive value (PPV) were computed according to the following equations:
S e n s i t i v i t y = T P T P + F N + D E · 100 ,
P P V = T P T P + F P + D E · 100 ,
where TP is the number of true positives, FN is the number of false negatives, FP is the number of false positives, and DE is the number of detection errors. The sensitivity quantifies the ability to correctly identify actual heartbeats in the NCC signals, by considering the ECG as the reference, while the PPV quantifies the ability to identify the actual heartbeats among all the detected NCC peaks. Specifically, sensitivity is defined as the percentage of NCC peaks corresponding to actual heartbeats that were correctly identified over the sum of the true and missed NCC peaks (it can also be expressed as the number of TP over the number of R-peaks, which represent reference heartbeats). On the other hand, PPV is defined as the percentage of NCC peaks corresponding to actual heartbeats that were correctly identified over the number of all the detected NCC peaks.
Furthermore, the inter-beat interval measurements from the NCC signals were compared with those estimated from the reference ECG signals by means of regression, correlation, and Bland–Altman analyses [84,85]. The statistical analyses were carried out via the MATLAB® function “bland-altman-and-correlation-plot” [86]. The inter-beat intervals related to the FN and DE were excluded from the analyses.

3. Results

3.1. NCC Signals with Case 1 Templates

Table 3 outlines the number of R-peaks and NCC peaks detected per subject in the ECG and NCC signals with case 1 templates, respectively. The table also reports the number of FPs, FNs, and DEs identified in the NCC signals and the number of inter-beat intervals considered in the statistical analyses. In detail, a total of 19,496 heartbeats were identified in the ECG signals, while 19,318 NCC peaks were detected in the NCC signals, of which 18,741 were TPs, 156 FPs, 334 FNs, and 421 DEs. Therefore, the proposed heartbeat detection approach achieved a sensitivity of 96% and a PPV of 97%. Furthermore, a total of 18,085 inter-beat intervals were compared via the regression, correlation, and Bland–Altman analyses. A slope of 0.997 and an intercept of 2.80 ms, with a R2 value greater than 0.999 (p-value < 0.0001), were obtained from the regression and correlation analyses (see Figure 6a). Moreover, the Bland–Altman analysis reported a non-significant bias (p-value = 0.36), with limits of agreement (LoA) of ±7.8 ms (see Figure 6b). The results of the statistical analyses are also summarized in Table 5.

3.2. NCC Signals with Case 2 Templates

Table 4 illustrates the number of R-peaks and NCC peaks detected per subject in the ECG and NCC signals with case 2 templates, respectively. The table also reports the number of FPs, FNs, and DEs in the NCC signals, as compared to the reference ECG, along with the number of compared inter-beat intervals. In detail, a total of 13,769 heartbeats were identified in the ECG signals, while 13,241 NCC peaks, of which 10,842 TPs, 524 FPs, 1049 FNs, and 1875 DEs, were detected in the NCC signals. Hence, the NCC signals scored a sensitivity of 79% and a PPV of 82%. Moreover, the statistical analyses were performed on a total of 8955 inter-beat intervals. The regression and correlation analyses (see Figure 7a) reported a slope of 0.990 and an intercept of 7.62 ms, with a R2 value greater than 0.99 (p-value < 0.0001). Furthermore, non-significant bias (p-value = 0.15) and LoA of ±19 ms resulted from the Bland–Altman analysis (see Figure 7b). The results of the statistical analyses are also outlined in Table 5.

4. Discussion

This study described a novel algorithm based on template matching for the detection of heartbeats in SCG signals without the support of concurrent ECG signals. The high inter-subject and intra-subject variability and the noisy nature of SCG signals make the task of ECG-free heartbeat detection rather troublesome. Indeed, the use of ensemble averaging is very common in SCG studies to unveil the underlying waveform corresponding to heartbeats. Unlike previous studies that proposed rather complex algorithms, this study focused on searching for matches between the whole SCG signal and a single heartbeat template selected from the signal itself, by taking advantage of the normalized cross-correlation as a similarity measure with very high robustness to noise.
Two classes of template were used, namely “case 1”, which included both the systolic and diastolic complexes, and “case 2”, which included the systolic complex only. The reason for using the “case 2” template lies in the absence of clear diastolic complexes in the SCG recordings of some patients, which could not be distinguished from the physiological background noise or exhibited significant morphological instabilities that could seriously impair the performance of template matching. According to the results of the statistical analyses, the “case 1” template provided far better performance as compared to the “case 2” template, with a sensitivity of 96% vs. 79%, a PPV of 97% vs. 82%, and limits of agreement of ±7.8 ms vs. ±19 ms, respectively. Indeed, the shapes of individual complexes, both systolic and diastolic, are not very peculiar, so the use of a template that includes only a single complex could lead to detection errors. On the contrary, using a template that includes both systolic and diastolic complexes, which rather represent a pattern that is peculiar to a single heartbeat, is in principle less prone to misdetections, and the experimental results of this study support this hypothesis.
Considering that the SCG signals analyzed in this study had been acquired from pathological subjects affected by one or more valvular heart diseases, ranging from moderate to severe stages, the results achieved with both the “case 1” and “case 2” templates are surprisingly good, all the more so if the low computational burden of the proposed template matching approach is taken into account. Indeed, the proposed approach essentially requires two simple steps, namely NCC computation and peak detection, which make it also reasonably suitable for real-time implementations without significant computational resource requirements. In light of the promising results obtained at the cost of few processing operations, it seems that the use of artificial intelligence approaches for heartbeat detection in SCG signals, which is undoubtedly an emerging trend in modern research, is not worth its price in terms of the computational burden. It is also important to underline that none of the methods described in the literature have been tested on such a large cohort of pathological subjects as the one considered in this study, therefore their performance has yet to be confirmed on actual patients with cardiac diseases, who certainly represent the actual target of cardiac monitoring.
The main limitation to the proposed approach is the need for the manual selection of the heartbeat template, which makes the overall method operator dependent. However, this study provides qualitative guidelines for template selection, so as to foster the reproducibility of the results. Indeed, a fully automated template selection algorithm would improve the feasibility and reliability of the proposed template matching approach. This would also promote its full implementation in wearable devices for SCG-based cardiac monitoring, as opposed to methods based on artificial intelligence that must rely on cloud computing services and telecommunication infrastructures. While the development of an automated template selection algorithm was out of the scope of this study, it could certainly be the subject of future studies on this topic.
Another limitation to this study is that the SCG signals were all recorded on motionless subjects in supine positions from a single specific point on the chest. It is known that SCG signal morphology depends on body posture and sensor location, and that motion artifacts can seriously hinder signal quality. Therefore, it is very important to assess the performance of the proposed template matching approach in these conditions.
An interesting aspect to be investigated is the suitability of the proposed approach for heart rate variability (HRV) analyses based on SCG. Indeed, Sieciński et al. have already shown in [66] that HRV indices estimated from SCG signals are highly correlated with those extracted from ECG, even in patients with VHDs from the same database considered in this study. However, they determined the heartbeat locations via an algorithm that took advantage of the availability of a simultaneously recorded ECG signal. Therefore, it is important to confirm that the measurement errors achieved by the template matching method proposed in the present study are low enough to ensure the reliable estimation of HRV indices.
Finally, the performance of the proposed approach should be assessed also on a study population that includes patients affected by different cardiac diseases, such as branch blocks, myocardial infarction, cardiac dyssynchrony, and heart failure, because the degree to which such pathologies may undermine the effectiveness of the proposed method is currently unknown and cannot be predicted.

Author Contributions

Conceptualization, E.A. and P.B.; methodology, J.C., S.P. and E.A.; software, S.P.; formal analysis, J.C., S.P. and E.A.; investigation, J.C., S.P., E.A., P.B. and D.E.; writing—original draft preparation, J.C., S.P. and E.A.; writing—review and editing, J.C., S.P., E.A., P.B. and D.E.; visualization, J.C. and S.P.; supervision, E.A. and P.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All relevant research data will be made available upon request after the publication of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Galli, A.; Montree, R.J.H.; Que, S.; Peri, E.; Vullings, R. An Overview of the Sensors for Heart Rate Monitoring Used in Extramural Applications. Sensors 2022, 22, 4035. [Google Scholar] [CrossRef] [PubMed]
  2. Hieronymi, U.; Kramme, R. Cardiovascular monitoring. In Springer Handbook of Medical Technology; Kramme, R., Hoffmann, K.P., Pozos, R.S., Eds.; Springer: Berlin/Heidelberg, Germany, 2011; pp. 955–970. [Google Scholar] [CrossRef]
  3. Polley, C.; Jayarathna, T.; Gunawardana, U.; Naik, G.; Hamilton, T.; Andreozzi, E.; Bifulco, P.; Esposito, D.; Centracchio, J.; Gargiulo, G. Wearable Bluetooth Triage Healthcare Monitoring System. Sensors 2021, 21, 7586. [Google Scholar] [CrossRef]
  4. Van Ravenswaaij-Arts, C.M.; Kollée, L.A.; Hopman, J.C.; Stoelinga, G.B.; van Geijn, H.P. Heart rate variability. Ann. Intern. Med. 1993, 118, 436–447. [Google Scholar] [CrossRef] [PubMed]
  5. Perret-Guillaume, C.; Joly, L.; Benetos, A. Heart Rate as a Risk Factor for Cardiovascular Disease. Prog. Cardiovasc. Dis. 2009, 52, 6–10. [Google Scholar] [CrossRef] [PubMed]
  6. Kleiger, R.E.; Stein, P.K.; Bigger, J.T., Jr. Heart rate variability: Measurement and clinical utility. Ann. Noninvas. Electrocardiol. 2005, 10, 88–101. [Google Scholar] [CrossRef]
  7. Cygankiewicz, I.; Zareba, W. Chapter 31—Heart rate variability. In Handbook of Clinical Neurology; Buijs, R.M., Swaab, D.F., Eds.; Elsevier: Amsterdam, The Netherlands, 2013; Volume 117, pp. 379–393. [Google Scholar] [CrossRef]
  8. Faust, O.; Hong, W.; Loh, H.W.; Xu, S.; Tan, R.S.; Chakraborty, S.; Barua, P.D.; Molinari, F.; Acharya, U.R. Heart rate variability for medical decision support systems: A review. Comput. Biol. Med. 2022, 145, 105407. [Google Scholar] [CrossRef]
  9. Hernando, D.; Roca, S.; Sancho, J.; Alesanco, Á.; Bailón, R. Validation of the Apple Watch for Heart Rate Variability Measurements during Relax and Mental Stress in Healthy Subjects. Sensors 2018, 18, 2619. [Google Scholar] [CrossRef]
  10. Böhm, M.; Reil, J.C.; Deedwania, P.; Kim, J.B.; Borer, J.S. Resting heart rate: Risk indicator and emerging risk factor in cardiovascular disease. Am. J. Med. 2015, 128, 219–228. [Google Scholar] [CrossRef]
  11. Stein, P.K.; Pu, Y. Heart rate variability, sleep and sleep disorders. Sleep Med. Rev. 2012, 16, 47–66. [Google Scholar] [CrossRef]
  12. Lundstrom, C.J.; Foreman, N.A.; Biltz, G. Practices and Applications of Heart Rate Variability Monitoring in Endurance Athletes. Int. J. Sport. Med. 2023, 44, 9–19. [Google Scholar] [CrossRef]
  13. Achten, J.; Jeukendrup, A.E. Heart rate monitoring: Applications and limitations. Sport. Med. 2003, 33, 517–538. [Google Scholar] [CrossRef]
  14. Zhu, J.; Ji, L.; Liu, C. Heart rate variability monitoring for emotion and disorders of emotion. Physiol. Meas. 2019, 40, 064004. [Google Scholar] [CrossRef] [PubMed]
  15. Lu, K.; Sjörs Dahlman, A.; Karlsson, J.; Candefjord, S. Detecting driver fatigue using heart rate variability: A systematic review. Accid. Anal. Prev. 2022, 178, 106830. [Google Scholar] [CrossRef] [PubMed]
  16. Clark, J.W.J. The origin of biopotentials. In Medical Instrumentation: Application and Design, 4th ed.; Webster, J.G., Ed.; John Wiley & Sons: Hoboken, NJ, USA, 2010; pp. 147–158. [Google Scholar]
  17. Nazeran, H. Electrocardiography, computers in. In Encyclopedia of Medical Devices and Instrumentation, 2nd ed.; Webster, J.G., Ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2006; pp. 34–51. [Google Scholar]
  18. Vieau, S.; Iaizzo, P.A. Basic ECG theory, 12-lead recordings and their interpretation. In Handbook of Cardiac Anatomy, Physiology and Devices, 3rd ed.; Iaizzo, P.A., Ed.; Springer: Cham, Switzerland, 2015; pp. 321–334. [Google Scholar]
  19. Berbari, E.J. Principles of electrocardiography. In Biomedical Engineering Handbook; Bronzino, J.D., Ed.; CRC Press: Boca Raton, FL, USA, 1999; pp. 231–240. [Google Scholar]
  20. Sattar, Y.; Chhabra, L. Electrocardiogram; StatPearls Publishing: Treasure Island, FL, USA, 2023. Available online: https://www.ncbi.nlm.nih.gov/books/NBK549803/ (accessed on 7 April 2023).
  21. Luz, E.J.; Schwartz, W.R.; Cámara-Chávez, G.; Menotti, D. ECG-based heartbeat classification for arrhythmia detection: A survey. Comput. Methods Programs Biomed. 2016, 127, 144–164. [Google Scholar] [CrossRef]
  22. Pan, J.; Tompkins, W.J. A real-time QRS detection algorithm. IEEE Trans. Biomed. Eng. 1985, 32, 230–236. [Google Scholar] [CrossRef]
  23. Allen, J. Photoplethysmography and its application in clinical physiological measurement. Physiol. Meas. 2007, 28, R1–R39. [Google Scholar] [CrossRef] [PubMed]
  24. Andreozzi, E.; Sabbadini, R.; Centracchio, J.; Bifulco, P.; Irace, A.; Breglio, G.; Riccio, M. Multimodal Finger PulseWave Sensing: Comparison of Forcecardiography and Photoplethysmography Sensors. Sensors 2022, 22, 7566. [Google Scholar] [CrossRef]
  25. Pereira, T.; Tran, N.; Gadhoumi, K.; Pelter, M.M.; Do, D.H.; Lee, R.J.; Colorado, R.; Meisel, K.; Hu, X. Photoplethysmography based atrial fibrillation detection: A review. NPJ Digit. Med. 2020, 3, 3. [Google Scholar] [CrossRef]
  26. Biswas, D.; Simões-Capela, N.; Van Hoof, C.; Van Helleputte, N. Heart Rate Estimation from Wrist-Worn Photoplethysmography: A Review. IEEE Sens. J. 2019, 19, 6560–6570. [Google Scholar] [CrossRef]
  27. Zhang, Z. Photoplethysmography-Based Heart Rate Monitoring in Physical Activities via Joint Sparse Spectrum Reconstruction. IEEE Trans. Biomed. Eng. 2015, 62, 1902–1910. [Google Scholar] [CrossRef]
  28. Ismail, S.; Akram, U.; Siddiqi, I. Heart rate tracking in photoplethysmography signals affected by motion artifacts: A review. EURASIP J. Adv. Signal. Process. 2021, 5, 5. [Google Scholar] [CrossRef]
  29. Gordon, J.W. On certain molar movements of the human body produced by the circulation of blood. J. Anat. Physiol. 1877, 11, 533–536. [Google Scholar] [PubMed]
  30. Knoop, A.A. Experimental investigations on ultra-low frequency displacement ballistocardiography: NASA TT F-269. In NASA Contractor Report—NASA CR; National Aeronautics and Space Administration: Washington, DC, USA, 1965; pp. 1–107. [Google Scholar]
  31. Starr, I. The relation of the ballistocardiogram to cardiac function. Am. J. Cardiol. 1958, 2, 737–747. [Google Scholar] [CrossRef]
  32. Burger, H.C.; Noordergraaf, A. Physical basis of ballistocardiography. III. Am. Heart J. 1956, 51, 179–185. [Google Scholar] [CrossRef]
  33. Sadek, I.; Biswas, J.; Abdulrazak, B. Ballistocardiogram signal processing: A review. Health Inf. Sci. Syst. 2019, 7, 10. [Google Scholar] [CrossRef] [PubMed]
  34. Inan, O.T. Recent advances in cardiovascular monitoring using ballistocardiography. In Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), San Diego, CA, USA, 28 August–1 September 2012; pp. 5038–5041. [Google Scholar] [CrossRef]
  35. Rappaport, M.B.; Sprague, H.B. The graphic registration of the normal heart sounds. Am. Heart J. 1942, 23, 591–623. [Google Scholar] [CrossRef]
  36. Dimond, E.G.; Benchimol, A. Phonocardiography. Calif. Med. 1961, 94, 139–146. [Google Scholar]
  37. Ismail, S.; Siddiqi, I.; Akram, U. Localization and classification of heart beats in phonocardiography signals—A comprehensive review. EURASIP J. Adv. Signal. Process. 2018, 2018, 26. [Google Scholar] [CrossRef]
  38. Giordano, N.; Knaflitz, M. A Novel Method for Measuring the Timing of Heart Sound Components through Digital Phonocardiography. Sensors 2019, 19, 1868. [Google Scholar] [CrossRef]
  39. Zanetti, J.; Salerno, D. Seismocardiography: A new technique for recording cardiac vibrations: Concept, method, and initial observations. J. Cardiovasc. Technol. 1990, 9, 2. [Google Scholar]
  40. Inan, O.T.; Migeotte, P.F.; Park, K.S.; Etemadi, M.; Tavakolian, K.; Casanella, R.; Zanetti, J.; Tank, J.; Funtova, I.; Prisk, G.K.; et al. Ballistocardiography and Seismocardiography: A Review of Recent Advances. IEEE J. Biomed. Health Inform. 2015, 19, 1414–1427. [Google Scholar] [CrossRef] [PubMed]
  41. Taebi, A.; Solar, B.E.; Bomar, A.J.; Sandler, R.H.; Mansy, H.A. Recent Advances in Seismocardiography. Vibration 2019, 2, 64–86. [Google Scholar] [CrossRef] [PubMed]
  42. Zanetti, J.M.; Tavakolian, K. Seismocardiography: Past, present and future. In Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Osaka, Japan, 3–7 July 2013; pp. 7004–7007. [Google Scholar] [CrossRef]
  43. Crow, R.S.; Hannan, P.; Jacobs, D.; Hedquist, L.; Salerno, D.M. Relationship between Seismocardiogram and Echocardiogram for Events in the Cardiac Cycle. Am. J. Noninvas. Cardiol. 1994, 8, 39–46. [Google Scholar] [CrossRef]
  44. Khosrow-Khavar, F.; Tavakolian, K.; Blaber, A.; Menon, C. Automatic and Robust Delineation of the Fiducial Points of the Seismocardiogram Signal for Noninvasive Estimation of Cardiac Time Intervals. IEEE Trans. Biomed. Eng. 2017, 64, 1701–1710. [Google Scholar] [CrossRef]
  45. Di Rienzo, M.; Vaini, E.; Castiglioni, P.; Merati, G.; Meriggi, P.; Parati, G.; Faini, A.; Rizzo, F. Wearable seismocardiography: Towards a beat-by-beat assessment of cardiac mechanics in ambulant subjects. Auton. Neurosci. 2013, 178, 50–59. [Google Scholar] [CrossRef]
  46. Laurin, A.; Blaber, A.; Tavakolian, K. Seismocardiograms return valid heart rate variability indices. Comput. Cardiol. 2013, 2013, 413–416. [Google Scholar]
  47. Jafari Tadi, M.; Lehtonen, E.; Saraste, A.; Tuominen, J.; Koskinen, J.; Teräs, M.; Airaksinen, J.; Pänkäälä, M.; Koivisto, T. Gyrocardiography: A New Non-invasive Monitoring Method for the Assessment of Cardiac Mechanics and the Estimation of Hemodynamic Variables. Sci. Rep. 2017, 7, 6823. [Google Scholar] [CrossRef]
  48. Sieciński, S.; Kostka, P.S.; Tkacz, E.J. Gyrocardiography: A Review of the Definition, History, Waveform Description and Applications. Sensors 2020, 20, 6675. [Google Scholar] [CrossRef]
  49. Hossein, A.; Mirica, D.C.; Rabineau, J.; Del Rio, J.I.J.; Morra, S.; Gorlier, D.; Nonclercq, A.; van de Borne, P.; Migeotte, P.F. Accurate Detection of Dobutamine-induced Haemodynamic Changes by Kino-Cardiography: A Randomised Double-Blind Placebo-Controlled Validation Study. Sci. Rep. 2019, 9, 10479. [Google Scholar] [CrossRef]
  50. Hossein, A.; Rabineau, J.; Gorlier, D.; Del Rio, J.I.J.; van de Borne, P.; Migeotte, P.F.; Nonclercq, A. Kinocardiography Derived from Ballistocardiography and Seismocardiography Shows High Repeatability in Healthy Subjects. Sensors 2021, 21, 815. [Google Scholar] [CrossRef]
  51. Andreozzi, E.; Fratini, A.; Esposito, D.; Naik, G.; Polley, C.; Gargiulo, G.D.; Bifulco, P. Forcecardiography: A Novel Technique to Measure Heart Mechanical Vibrations onto the Chest Wall. Sensors 2020, 20, 3885. [Google Scholar] [CrossRef]
  52. Andreozzi, E.; Gargiulo, G.D.; Esposito, D.; Bifulco, P. A Novel Broadband Forcecardiography Sensor for Simultaneous Monitoring of Respiration, Infrasonic Cardiac Vibrations and Heart Sounds. Front. Physiol. 2021, 12, 725716. [Google Scholar] [CrossRef] [PubMed]
  53. Andreozzi, E.; Centracchio, J.; Punzo, V.; Esposito, D.; Polley, C.; Gargiulo, G.D.; Bifulco, P. Respiration Monitoring via Forcecardiography Sensors. Sensors 2021, 21, 3996. [Google Scholar] [CrossRef]
  54. Centracchio, J.; Andreozzi, E.; Esposito, D.; Gargiulo, G.D.; Bifulco, P. Detection of Aortic Valve Opening and Estimation of Pre-Ejection Period in Forcecardiography Recordings. Bioengineering 2022, 9, 89. [Google Scholar] [CrossRef] [PubMed]
  55. Andreozzi, E.; Centracchio, J.; Esposito, D.; Bifulco, P. A Comparison of Heart Pulsations Provided by Forcecardiography and Double Integration of Seismocardiogram. Bioengineering 2022, 9, 167. [Google Scholar] [CrossRef]
  56. Centracchio, J.; Andreozzi, E.; Esposito, D.; Gargiulo, G.D. Respiratory-Induced Amplitude Modulation of Forcecardiography Signals. Bioengineering 2022, 9, 444. [Google Scholar] [CrossRef]
  57. Centracchio, J.; Esposito, D.; Gargiulo, G.D.; Andreozzi, E. Changes in Forcecardiography Heartbeat Morphology Induced by Cardio-Respiratory Interactions. Sensors 2022, 22, 9339. [Google Scholar] [CrossRef] [PubMed]
  58. D’Mello, Y.; Skoric, J.; Xu, S.; Roche, P.J.R.; Lortie, M.; Gagnon, S.; Plant, D.V. Real-Time Cardiac Beat Detection and Heart Rate Monitoring from Combined Seismocardiography and Gyrocardiography. Sensors 2019, 19, 3472. [Google Scholar] [CrossRef] [PubMed]
  59. Dehkordi, P.; Tavakolian, K.; Tadi, M.J.; Zakeri, V.; Khosrow-Khavar, F. Investigating the estimation of cardiac time intervals using gyrocardiography. Physiol. Meas. 2020, 41, 055004. [Google Scholar] [CrossRef] [PubMed]
  60. Laurin, A.; Khosrow-Khavar, F.; Blaber, A.P.; Tavakolian, K. Accurate and consistent automatic seismocardiogram annotation without concurrent ECG. Physiol. Meas. 2016, 37, 1588–1604. [Google Scholar] [CrossRef] [PubMed]
  61. Jafari Tadi, M.; Lehtonen, E.; Hurnanen, T.; Koskinen, J.; Eriksson, J.; Pänkäälä, M.; Teräs, M.; Koivisto, T. A real-time approach for heart rate monitoring using a Hilbert transform in seismocardiograms. Physiol. Meas. 2016, 37, 1885–1909. [Google Scholar] [CrossRef] [PubMed]
  62. Sieciński, S.; Kostka, P. Determining heart rate beat-to-beat from smartphone seismocardiograms: Preliminary studies. In Proceedings of the Innovations in Biomedical Engineering Conference, Zabrze, Poland, 19–20 October 2017; Advances in Intelligent Systems and Computing. Gzik, M., Tkacz, E., Paszenda, Z., Piętka, E., Eds.; Springer: Cham, Switzerland, 2018; 623, pp. 133–140. [Google Scholar] [CrossRef]
  63. Kaisti, M.; Jafari Tadi, M.; Lahdenoja, O.; Hurnanen, T.; Saraste, A.; Pänkäälä, M.; Koivisto, T. Stand-Alone Heartbeat Detection in Multidimensional Mechanocardiograms. IEEE Sens. J. 2019, 19, 234–242. [Google Scholar] [CrossRef]
  64. Luu, L.; Dinh, A. Using moving average method to recognize systole and diastole on seismocardiogram without ECG signal. In Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Honolulu, HI, USA, 18–21 July 2018; pp. 3796–3799. [Google Scholar] [CrossRef]
  65. Choudhary, T.; Sharma, L.N.; Bhuyan, M.K. Automatic Detection of Aortic Valve Opening Using Seismocardiography in Healthy Individuals. IEEE J. Biomed. Health Inform. 2019, 23, 1032–1040. [Google Scholar] [CrossRef] [PubMed]
  66. Sieciński, S.; Tkacz, E.J.; Kostka, P.S. Comparison of HRV indices obtained from ECG and SCG signals from CEBS database. Biomed. Eng. Online 2019, 18, 69. [Google Scholar] [CrossRef] [PubMed]
  67. Erfianto, B.; Rizal, A.; Suryani, V. Comparison of Seismocardiography Based Heart Rate Measurement Method. J. Southwest. Jiaotong Univ. 2020, 55, 1–15. [Google Scholar] [CrossRef]
  68. Xia, Z.; Shandhi, M.M.H.; Li, Y.; Inan, O.T.; Zhang, Y. The Delineation of Fiducial Points for Non-Contact Radar Seismocardiogram Signals without Concurrent ECG. IEEE J. Biomed. Health Inform. 2021, 25, 1031–1040. [Google Scholar] [CrossRef]
  69. Mora, N.; Cocconcelli, F.; Matrella, G.; Ciampolini, P. Fully Automated Annotation of Seismocardiogram for Noninvasive Vital Sign Measurements. IEEE Trans. Instrum. Meas. 2020, 69, 1241–1250. [Google Scholar] [CrossRef]
  70. García-González, M.A.; Argelagós-Palau, A.; Fernández-Chimeno, M.; Ramos-Castro, J. A comparison of heartbeat detectors for the seismocardiogram. In Proceedings of the Computing in Cardiology 2013, Zaragoza, Spain, 22–25 September 2013; pp. 461–464. [Google Scholar]
  71. Ferdinando, H.; Seppälä, E.; Myllylä, T. Discrete Wavelet Transforms-Based Analysis of Accelerometer Signals for Continuous Human Cardiac Monitoring. Appl. Sci. 2021, 11, 12072. [Google Scholar] [CrossRef]
  72. Choudhary, T.; Das, M.; Sharma, L.N.; Bhuyan, M.K. Analyzing seismocardiographic approach for heart rate variability measurement. Biomed. Signal. Process. Control. 2021, 68, 102793. [Google Scholar] [CrossRef]
  73. Scarpetta, M.; Spadavecchia, M.; Andria, G.; Ragolia, M.A.; Giaquinto, N. Accurate simultaneous measurement of heartbeat and respiratory intervals using a smartphone. J. Instrum. 2022, 17, P07020. [Google Scholar] [CrossRef]
  74. Wahlstrom, J.; Skog, I.; Handel, P.; Khosrow-Khavar, F.; Tavakolian, K.; Stein, P.K.; Nehorai, A. A Hidden Markov Model for Seismocardiography. IEEE Trans. Biomed. Eng. 2017, 64, 2361–2372. [Google Scholar] [CrossRef] [PubMed]
  75. Thakkar, H.K.; Sahoo, P.K. Towards Automatic and Fast Annotation of Seismocardiogram Signals Using Machine Learning. IEEE Sens. J. 2020, 20, 2578–2589. [Google Scholar] [CrossRef]
  76. Chen, Y.; Xu, W.; Zhu, W.; Ma, G.; Chen, X.; Wang, L. Beat-to-beat heart rate detection based on seismocardiogram using BiLSTM network. In Proceedings of the IEEE 20th International Conference on Trust, Security and Privacy in Computing and Communications (TrustCom), Shenyang, China, 20–22 October 2021; pp. 1503–1507. [Google Scholar] [CrossRef]
  77. Duraj, K.M.; Sieciński, S.; Doniec, R.J.; Piaseczna, N.J.; Kostka, P.S.; Tkacz, E.J. Heartbeat detection in seismocardiograms with semantic segmentation. In Proceedings of the 44th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Glasgow, UK, 11–15 July 2022; pp. 662–665. [Google Scholar] [CrossRef]
  78. Işilay Zeybek, Z.M.; Racca, V.; Pezzano, A.; Tavanelli, M.; Di Rienzo, M. Can Seismocardiogram Fiducial Points Be Used for the Routine Estimation of Cardiac Time Intervals in Cardiac Patients? Front. Physiol. 2022, 13, 825918. [Google Scholar] [CrossRef]
  79. Centracchio, J.; Muto, V. Heartbeats localization in forcecardiography signals via template matching. In Proceedings of the E-Health and Bioengineering Conference (EHB), Iasi, Romania, 17–18 November 2022; pp. 1–4. [Google Scholar] [CrossRef]
  80. Shin, J.H.; Choi, B.H.; Lim, Y.G.; Jeong, D.U.; Park, K.S. Automatic ballistocardiogram (BCG) beat detection using a template matching approach. In Proceedings of the 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Vancouver, BC, Canada, 20–25 August 2008; pp. 1144–1446. [Google Scholar] [CrossRef]
  81. Chen, Y.H.; Chen, H.H.; Chen, T.C.; Chen, L.G. Robust heart rate measurement with phonocardiogram by on-line template extraction and matching. In Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Boston, MA, USA, 30 August–3 September 2011; pp. 1957–1960. [Google Scholar] [CrossRef]
  82. Yang, C.; Fan, F.; Aranoff, N.; Green, P.; Li, Y.; Liu, C.; Tavassolian, N. An Open-Access Database for the Evaluation of Cardio-Mechanical Signals from Patients with Valvular Heart Diseases. Front. Physiol. 2021, 12, 750221. [Google Scholar] [CrossRef] [PubMed]
  83. Sedghamiz, H. BioSigKit: A Matlab Toolbox and Interface for Analysis of BioSignals. J. Open. Source Softw. 2018, 3, 671. [Google Scholar] [CrossRef]
  84. Altman, D.G.; Bland, J.M. Measurement in medicine: The analysis of method comparison studies. Statistician 1983, 32, 307–317. [Google Scholar] [CrossRef]
  85. Giavarina, D. Understanding Bland Altman analysis. Biochem. Med. 2015, 25, 141–151. [Google Scholar] [CrossRef] [PubMed]
  86. Ran, K. Bland-Altman and Correlation Plot, MATLAB Central File Exchange. 2020. Available online: https://www.mathworks.com/matlabcentral/fileexchange/45049-bland-altman-and-correlation-plot (accessed on 15 March 2023).
Figure 1. An excerpt ofthe pre-processed SCG (black line) and ECG (blue line) signals from subject #UP-21. The blue points mark the locations of the R-peaks.
Figure 1. An excerpt ofthe pre-processed SCG (black line) and ECG (blue line) signals from subject #UP-21. The blue points mark the locations of the R-peaks.
Sensors 23 04684 g001
Figure 2. Some examples of the ECG signals (blue lines), SCG signals (black lines), and the selected templates (red lines) from: (a) subject #UP-21; (b) subject #UP-11. The template in panel (a) comprises both the systolic and diastolic complexes in the cardiac cycle, while the one in panel (b) includes only the systolic complex.
Figure 2. Some examples of the ECG signals (blue lines), SCG signals (black lines), and the selected templates (red lines) from: (a) subject #UP-21; (b) subject #UP-11. The template in panel (a) comprises both the systolic and diastolic complexes in the cardiac cycle, while the one in panel (b) includes only the systolic complex.
Sensors 23 04684 g002
Figure 3. Some excerpts of the ECG (blue line), SCG (black line), and NCC (orange line) signals from: (a) subject #UP-21; (b) subject #UP-11. The blue and black points mark the locations of the R-peaks and NCC peaks, respectively. The NCC signal in panel (a) was obtained by considering the template depicted in Figure 2a, which comprised both the systolic and diastolic complexes in the selected cardiac cycle. The NCC signal in panel (b), instead, was obtained by considering the template in Figure 2b, which consisted only of the systolic complex of the selected heartbeat.
Figure 3. Some excerpts of the ECG (blue line), SCG (black line), and NCC (orange line) signals from: (a) subject #UP-21; (b) subject #UP-11. The blue and black points mark the locations of the R-peaks and NCC peaks, respectively. The NCC signal in panel (a) was obtained by considering the template depicted in Figure 2a, which comprised both the systolic and diastolic complexes in the selected cardiac cycle. The NCC signal in panel (b), instead, was obtained by considering the template in Figure 2b, which consisted only of the systolic complex of the selected heartbeat.
Sensors 23 04684 g003
Figure 4. An example of the NCC peak detection from subject #CP-36: (a) false positive; (b) false negative; (c) detection error.
Figure 4. An example of the NCC peak detection from subject #CP-36: (a) false positive; (b) false negative; (c) detection error.
Sensors 23 04684 g004aSensors 23 04684 g004b
Figure 5. An example of the ECG (blue line), SCG (black line), and NCC (orange line) signals from subject #UP-21. The blue and black points mark the locations of the R-peaks and NCC peaks, respectively. The temporal distance between two NCC peaks (double arrow green line) is referred to as the inter-beat interval.
Figure 5. An example of the ECG (blue line), SCG (black line), and NCC (orange line) signals from subject #UP-21. The blue and black points mark the locations of the R-peaks and NCC peaks, respectively. The temporal distance between two NCC peaks (double arrow green line) is referred to as the inter-beat interval.
Sensors 23 04684 g005
Figure 6. Statistical analyses on the inter-beat intervals: (a) results of regression and correlation analyses; (b) results of Bland–Altman analysis for the NCC signals with case 1 templates.
Figure 6. Statistical analyses on the inter-beat intervals: (a) results of regression and correlation analyses; (b) results of Bland–Altman analysis for the NCC signals with case 1 templates.
Sensors 23 04684 g006
Figure 7. Statistical analyses on the inter-beat intervals: (a) results of regression and correlation analyses; (b) results of Bland–Altman analysis for the NCC signals with case 2 templates.
Figure 7. Statistical analyses on the inter-beat intervals: (a) results of regression and correlation analyses; (b) results of Bland–Altman analysis for the NCC signals with case 2 templates.
Sensors 23 04684 g007
Table 1. Patient IDs of discarded SCG signals, along with the reasons for exclusion.
Table 1. Patient IDs of discarded SCG signals, along with the reasons for exclusion.
Patient ID #Reason for Exclusion
CP-03, CP-06, CP-17, CP-18, CP-24, CP-25, CP-29, CP-31, CP-35, CP-46, CP-50, CP-51, CP-54, CP-62, CP-67, UP-02, UP-03, UP-05, UP-19, UP-22, UP-25, UP-26Poor SCG signal quality
UP-28Simultaneous ECG signal not acquired
Table 2. Template selection criteria. The templates depicted in the table correspond to those highlighted in red in Figure 2.
Table 2. Template selection criteria. The templates depicted in the table correspond to those highlighted in red in Figure 2.
CaseDescriptionTypical Template Waveform
1The template starts from 2–3 oscillations before the systolic peak (highest local maximum), where the amplitude of the oscillations is significantly reduced with respect to the systolic peak amplitude; the template ends just after the last oscillation of the diastolic complex.Sensors 23 04684 i001
2The template starts and ends at about 1–2 oscillations before and after the highest systolic peak.Sensors 23 04684 i002
Table 3. The number of R-peaks and NCC peaks detected per subject in the ECG and NCC signals by considering case 1 templates, together with the FPs, FNs, and DEs identified in the NCC signals and the number of compared inter-beat intervals.
Table 3. The number of R-peaks and NCC peaks detected per subject in the ECG and NCC signals by considering case 1 templates, together with the FPs, FNs, and DEs identified in the NCC signals and the number of compared inter-beat intervals.
Patient
ID #
R-PeaksNCC PeaksFPFNDECompared Inter-Beat Intervals
CP-01448448000447
CP-05509507358487
CP-085445380618502
CP-1242343314413388
CP-15630623073610
CP-20391391110388
CP-21247247221241
CP-27130130000129
CP-28238237015225
CP-305235100132498
CP-3344945612512417
CP-34462460464444
CP-36386377817101166
CP-39518518008504
CP-41349344162332
CP-44321320230314
CP-4535734301417303
CP-4753752262111476
CP-4945145711547353
CP-53562561010559
CP-57507506015495
CP-58525524011520
CP-59405406100404
CP-60512503091493
CP-61397403601395
CP-63610608020605
CP-64382391901379
CP-65369366031360
CP-66468468110465
CP-68327330232032231
CP-69587585020582
CP-7042241231333338
UP-042582515122233
UP-0645840015928312
UP-07376373142365
UP-08257256011252
UP-09286278080273
UP-101651574122136
UP-1310693013184
UP-14340338134325
UP-152141890250168
UP-162212141818172
UP-18350349011345
UP-20617615022608
UP-21305305002300
UP-23565567200564
UP-243493402119311
UP-2726928316221224
UP-29146142041136
UP-302282441600227
Total19,49619,31815633442118,085
Table 4. The number of R-peaks and NCC peaks detected per subject in the ECG and NCC signals by considering case 2 templates, together with FPs, FNs, and DEs identified in the NCC signals and the number of compared inter-beat intervals.
Table 4. The number of R-peaks and NCC peaks detected per subject in the ECG and NCC signals by considering case 2 templates, together with FPs, FNs, and DEs identified in the NCC signals and the number of compared inter-beat intervals.
Patient
ID #
R-PeaksNCC PeaksFPFNDECompared
Inter-Beat
Intervals
CP-02651643182619571
CP-046616541846556
CP-07451453315441
CP-093643421537105139
CP-105065056364113213
CP-11656662262056525
CP-13837829085814
CP-144724773631163157
CP-16355360151013312
CP-1948449318930414
CP-22610595163151466
CP-2323522501028172
CP-263893815135353
CP-32527486387995225
CP-37342335556290104
CP-3840642018446309
CP-405095052630233114
CP-42346340273341220
CP-43460478371995266
CP-48637626162750500
CP-52728532019638298
CP-5579360215206248127
CP-567426931561145409
UP-01239241201831154
UP-11417414172054305
UP-12339347231564208
UP-176136031116583
Total13,76913,241524104918758955
Table 5. Results of the statistical analyses. Non-significant bias is indicated as “NS”.
Table 5. Results of the statistical analyses. Non-significant bias is indicated as “NS”.
Template CaseCompared
Inter-Beat
Intervals
R2SlopeIntercept
(ms)
Biasp-ValueLoA
(ms)
Sensitivity (%)PPV (%)
118,085>0.9990.9972.80NS0.36±7.89697
28955>0.990.9907.62NS0.15±197982
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Centracchio, J.; Parlato, S.; Esposito, D.; Bifulco, P.; Andreozzi, E. ECG-Free Heartbeat Detection in Seismocardiography Signals via Template Matching. Sensors 2023, 23, 4684. https://doi.org/10.3390/s23104684

AMA Style

Centracchio J, Parlato S, Esposito D, Bifulco P, Andreozzi E. ECG-Free Heartbeat Detection in Seismocardiography Signals via Template Matching. Sensors. 2023; 23(10):4684. https://doi.org/10.3390/s23104684

Chicago/Turabian Style

Centracchio, Jessica, Salvatore Parlato, Daniele Esposito, Paolo Bifulco, and Emilio Andreozzi. 2023. "ECG-Free Heartbeat Detection in Seismocardiography Signals via Template Matching" Sensors 23, no. 10: 4684. https://doi.org/10.3390/s23104684

APA Style

Centracchio, J., Parlato, S., Esposito, D., Bifulco, P., & Andreozzi, E. (2023). ECG-Free Heartbeat Detection in Seismocardiography Signals via Template Matching. Sensors, 23(10), 4684. https://doi.org/10.3390/s23104684

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

Article Metrics

Back to TopTop