
==== Front
Digit Health
Digit Health
DHJ
spdhj
Digital Health
2055-2076
SAGE Publications Sage UK: London, England

10.1177/20552076241277746
10.1177_20552076241277746
Original Research Article
Solving variability: Accurately extracting feature components from ballistocardiograms
https://orcid.org/0000-0002-0711-919X
Yang Tianyi 1
Yuan Haihang 1
Yang Junqi 1
Zhou Zhongchao 1
Abe Masayuki 2
Nakayama Yoshitake 3
Huang Shao Ying 4
Yu Wenwei 15
1 Department of Medical Engineering, 12737 Chiba University , Chiba City, Chiba Prefecture, Japan
2 Nanayume Co. Ltd, Chiba City, Chiba Prefecture, Japan
3 Center for Preventive Medical Sciences, 12737 Chiba University , Chiba City, Chiba Prefecture, Japan
4 Engineering Product Development, 233793 Singapore University of Technology and Design , Singapore, Singapore
5 Center for Frontier Medical Engineering, 12737 Chiba University , Chiba City, Chiba Prefecture, Japan
Wenwei Yu, Department of Medical Engineering, Natural Sciences Research Building 1-5F, Room 510, Chiba University, Nishi-Chiba Campus, 1-33 Yayoi-cho, Inage-ku, Chiba City, Chiba Prefecture 263-8522, Japan. Email: yuwill@faculty.chiba-u.jp
5 9 2024
Jan-Dec 2024
10 2055207624127774611 3 2024
8 8 2024
© The Author(s) 2024
2024
SAGE Publications Ltd, unless otherwise noted. Manuscript content on this site is licensed under Creative Commons Licenses
https://creativecommons.org/licenses/by-nc/4.0/ This article is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License (https://creativecommons.org/licenses/by-nc/4.0/) which permits non-commercial use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access page (https://us.sagepub.com/en-us/nam/open-access-at-sage).
Objective

A ballistocardiogram (BCG) is a vibration signal generated by the ejection of the blood in each cardiac cycle. The BCG has significant variability in amplitude, temporal aspects, and the deficiency of waveform components, attributed to individual differences, instantaneous heart rate, and the posture of the person being measured. This variability may make methods of extracting J-waves, the most distinct components of BCG less generalizable so that the J-waves could not be precisely localized, and further analysis is difficult. This study is dedicated to solving the variability of BCG to achieve accurate feature extraction.

Methods

Inspired by the generation mechanism of the BCG, we proposed an original method based on a profile of second-order derivative of BCG waveform (2ndD-P) to capture the nature of vibration and solve the variability, thereby accurately localizing the components especially when the J-wave is not prominent.

Results

In this study, 51 recordings of resting state and 11 recordings of high-heart-rate from 24 participants were used to validate the algorithm. Each recording lasts about 3 min. For resting state data, the sensitivity and positive predictivity of proposed method are: 98.29% and 98.64%, respectively. For high-heart-rate data, the proposed method achieved a performance comparable to those of low-heart-rate: 97.14% and 99.01% for sensitivity and positive predictivity, respectively.

Conclusion

Our proposed method can detect the peaks of the J-wave more accurately than conventional extraction methods, under the presence of different types of variability. Higher performance was achieved for BCG with non-prominent J-waves, in both low- and high-heart-rate cases.

Ballistocardiogram
second-order derivate
dynamic time warping
component extraction
J-wave
typesetterts19
cover-dateJanuary-December 2024
==== Body
pmcIntroduction

A ballistocardiogram (BCG) is a vibration signal generated by the ejection of blood in each cardiac cycle 1 and can be used to monitor the cardiovascular system. 2 Several studies have shown that BCGs can be used for monitoring the cardiac system, 3 cardiovascular diagnosis,4,5 and sleep monitoring. 6 BCGs can be collected using mechanical sensors, such as piezoelectric sheet sensor 7 and accelerometers. 8 Compared to piezoelectric sheet sensors, the measurement of BCG using accelerometer mostly requires the user to wear the sensor, for instance, in the form of a wristwatch. 9 The piezoelectric sheet sensor generates voltage by vibration or bending, therefore the BCG signal collected is vibration displacement. Since this type of BCG sensors can be integrated with furniture, such as chair cushions and mattresses, to monitor BCG signals in daily life, such that users are not required to wear any sensors, it does not impose additional burden on users, which makes it easier to use in non-clinical daily life environments. Therefore, it is very practical to use piezoelectric sheet sensor to measure BCG to realize some applications and this study focuses on BCG collected by piezoelectric sheet sensor.

Many applications based on the BCG require the precise extraction of its various components. As shown in Figure 1, a standard BCG waveform is characterized by distinct components labeled as H, I, J, K, and L waves. The IJK wave group are considered to originate during the systolic period of the heart. These waves, in conjunction with H and L waves, create a pattern notably resembling the letter “W” representing the most identifiable group of features within the BCG waveform. 2 In some applications of BCGs, information regarding the feature components is critical. For example, in Ref. 10 , the difference in amplitude between the I-J and J-K segments was an important parameter for estimating blood pressure.

Figure 1. Example of a standard BCG from piezoelectric sheet sensor using in this research.

In this field, the J-wave tends to be the most prominent component in the BCG, and after obtaining the location of the peak of J-wave, information on other BCG components is easy to obtain. 11 Thus, researchers have focused on methods of accurately extracting J-waves. 12 However, owing to the significant variability of the BCG waveforms, many J-wave-extraction methods have limitations.

Because the J-wave is the most prominent component generally, several methods have been proposed. The simplest method of detecting J-waves is to use a threshold that is computed automatically using a heart rate (HR) algorithm. 13 Some studies have used the length transform on the BCG and then search for local maxima to detect J-waves. 14 Additionally, the short-time energy (STE) of the BCG can be used to detect J-waves. 15 However, owing to waveform variability, the J-wave is not necessarily the most prominent component. Therefore, these methods are less effective in certain scenarios. Researchers have also used pattern recognition-based methods to seek similar components in each heartbeat cycle to accurately determine the HR. 16 These methods attempt to circumvent the discussion of BCG waveform variability; therefore, it is not certain whether the observed component is the J-wave. Recently, some researchers have used neural networks (NNs).17–19 However, NNs require training with a large amount of data, and even if there are enough data, there may still be the problem of imbalanced dataset. This may lead to biased learning of NNs. Moreover, the BCG variability has not been the focus of those studies. Even if the variability problem could be partially resolved relying on the generality of NNs, however, their explainability is the problem.

According to existing literatures and experimental observations, there are three main types of BCG waveform variability. The first is the variability in the relative amplitude (Va).20–22 As shown in Figure 2, the relative amplitude of the different components is likely to be different for different individuals. Even for the same individual, posture can affect the waveform. In some cases, the J-wave is not the most prominent component. The second is the variability in the time of component occurrence (Vt). 23 Each heartbeat is not exactly the same; therefore, differences would also occur in the timing of the occurrence of the components of the BCG. The third is the variability in the deficiency of the waveform components (Vd). As the HR increases, the components of the BCG become aliased and occasionally indistinguishable.16,17 Moreover, some cardiovascular diseases may also cause changes in the heartbeat, resulting in this variability. 24

Figure 2. Variability of BCG waveforms. Significant difference between the waveforms of two participants, a and b. Waveforms collected from the same participant, b, are different depending on sitting posture, (2) in sitting posture 1 and (3) in sitting posture 2. (4) shows the waveform of participant b at 140 bpm. Compared with (2) and (3), Vt and Vd are apparent.

Although the Va is apparent, the relative amplitudes of the BCG components are relatively consistent for the same measurement of the same participant. 25 Therefore, the use of a template-based approach is more feasible.22,26,27 Creating a template is crucial for template-based methods. Using predefined templates is unlikely to handle a full range of scenarios, and it is better to extract templates from the BCG signal itself. This requires a method that can solve Va to determine the location of the main features of the BCG better than the traditional method represented by STE.

In terms of physiological signal processing, using second-order derivative is the classical operation to obtain additional information. 28 In recent years, studies have been conducted on the physiological mechanism of BCG generation. A research group used invasive intravascular pressure data to build a mathematical model and successfully predicted BCG waveforms. 29 As disclosed in, 29 a J-wave occurs when a maximum blood pressure gap exists between the apex and outlet of the aorta. The second-order derivative of the vibration displacement BCG at the measured site contains significant information about the pressure. Therefore, it can better capture the vibration feature of the cardiovascular system to control the Va to achieve a lock on the main features of the BCG.

By locking the J-wave in the template, only the point-to-point correspondence between the template and each signal segment containing the main components need be determined to extract all J-waves. Vt must be addressed to obtain accurate point-to-point correspondence.

Dynamic time warping (DTW) is an effective method of computing the optimal alignment between two time sequences.30,31 It can be applied to signals that change in the temporal domain but have consistent signal components. The occurrence time of the feature components of the BCG also change, but these components have some consistency in their relative magnitude. 25 Therefore, DTW is considered effective for solving Vt.

In the experimental data acquired in this study, Vd is clearly observed in a high-heart-rate BCG,16,17 where the IJK wave group occurring in systole is noticeable, but the other components may be aliased owing to the small heartbeat interval, resulting in a deficiency in the waveform components. Very few studies have tested high-heart-rate BCG signals. In ref. 17, three existing methods were tested; however, these methods performed poorly because of Vd. In this study, high-heart-rate cases were carefully investigated and compared with low-heart-rate cases to identify the differences and classify the high- and low-heart-rate cases for further analysis. In many physiological signal analysis methods, adjusting the parameters of the method according to the different conditions of the signals under investigation may also play a key role. 32 By experiments, it is clear that 2ndD-Ps of different rates of heart beats have obvious difference in frequency spectrum. Therefore, in this study, based on the frequency spectrum of the 2ndD-P, the Low- and High-Heart-Rate classification is developed to adjust the parameters.

Based on the above considerations, this paper proposes a novel method for identifying the peak of the J-wave, which is a salient component of BCGs. The 2ndD-P, template, DTW, and Low- and High-Heart-Rate Classifications based on the 2ndD-P frequency spectrum were combined to address the BCG waveform variability. Resting-state and post-exercise signals (high-heart-rate cases) were collected to validate the proposed method. Through experiments, it can be clearly seen that 2ndD-P, template, DTW, and low- and high-heart-rate classification play indispensable and important roles. Compared with the traditional method represented by STE, this method has more advantages in the sensitivity and positive predictivity and calculation of HRV parameters. In addition, the filtering strategy in signal pre-processing also has an impact on feature extraction. 33 To make sure that the proposed method does not depend on filtering strategies, preprocessing using different filtering strategies is also performed.

The remainder of this paper is organized as follows. The proposed method is described in detail in Methods section. Experiments section presents the experimental design of this study. Results section presents the results of the study. A discussion of the results is presented in Discussion section. Finally, Conclusion section concludes the paper.

Methods

The nature of this study is the processing and feature extraction of physiological signals with different types of variability.

Process of the proposed method

The schematic diagram of the proposed method is shown in Figure 3. For each BCG section, a small portion of the signals was first extracted as a trial segment (Seg-T). The 2ndD-P is then obtained for the Seg-T with default parameters, and used as the input for the Low- and High-Heart-Rate Classification, to adjust the parameters for the further analysis of BCG signals of different rates of heart beats. The 2ndD-P of the whole BCG section is acquired with the adjusted parameters. After that, the BCG section then is divided to pulse segments, each of which is aligned to be centered at its corresponding local maximum point in the 2ndD-P of the BCG section. Each pulse segment has same length, supposed to contain IJK wave group of a BCG pulse. A template is further generated for the BCG section by averaging the pulse segments, and used to localize the J-waves through template matching with DTW, which could deal with the Vt.

Figure 3. Schematic diagram of the proposed method.

As shown, the 2ndD-P, low- and high-heart-rate classification, template matching with DTW work cohesively in the proposed method for dealing with different types of BCG variability. The details of the proposed method are as follows.

Pre-processing

In this study, signals were collected at a sampling rate of 1000 Hz. To reduce the computational burden, we first decreased the sampling rate to 250 Hz. The spectrum energy of the BCG is considered to be concentrated at 10 Hz.15,34 Wavelet decomposition and reconstruction methods were used to remove noise. In this method, to avoid losing the original signal features, we selected the Sym8 wavelet, which is close to the shape of the BCG, to perform a 7-layer wavelet decomposition. 13 Subsequently, the detailed coefficients containing high-frequency noise and the approximation containing the baseline shift and respiration were removed, and signals ranging from 1 to 15.6 Hz were preserved. Owing to the individual differences among the participants, a significant difference was observed in the amplitude of the original signal. Therefore, the signal was standardized to eliminate this difference.

2ndD-P

For the pre-processed BCG signal, the Savitzky-Golay (SG) filter is first used to obtain the second-order derivative of the BCG signal. 35 In this study, to preserve the features of the BCG as much as possible and eliminate noise with sharp shapes, we set the window width and poly-order of the SG filter to 0.18 s and 3, respectively. As shown in Figure 4, the result revealed many components with high second-order derivative value near the J-wave position, even though the J-wave of the signal is not the most apparent component. Therefore, the positions of the IJK wave group in the BCG were effectively highlighted by the second-order derivative waveform. Then, the second-order derivative waveform must be simplified to a profile.

Figure 4. BCG and the absolute value of its second-order derivative waveform (au: arbitrary units). The absolute values of second-order derivative are used for highlighting the IJK wave group, which has significantly large values than other parts of BCG, even in the case that the J-wave was not the most prominent component.

For the second-order derivative waveform, referring to the method in, 15 the STE transformation was applied to obtain its profile using a sliding window. For waveform xi , in the sliding window, the short-time energy Ei is defined as follows: (1) Ei=∑n=0NA−1xi(n)2

where NA is a suitable segment size to ensure that the IJK wave group is included in segment xi . In the implementation, a sliding window of length NA with a step size of one point was used to scan all second-order derivative waveform to obtain the corresponding profile, 2ndD-P, as shown in Figure 5.

Figure 5. BCG and its 2ndD-P. The IJK wave group position can be located through the LIJK .

Subsequently, the local maximum points of the 2ndD-P, LIJK , considered to be the position with the highest energy of the pulse cycle, are most likely to correspond to the center of IJK wave group. LIJK are determined as:

in∈N exists such that for all i with |i−in|<WL , Ei<Ei0 . (2) LIJK={i1,i2,…,in}

As mentioned in the introduction, the parameters NA and WL must be determined according to the classification of the HR to achieve better results, which is described in the next section.

Low- and high-heart-rate classification

As described in the introduction, the parameters ( NA,WL ) must be adjusted according to the HR to have a better effect. In this study, a HR below 90 bpm was considered a low HR. When the HR is low, the interval between heartbeats is at least 0.66 s, WL was set to 0.6 s. The NA was set to 0.24 s based on the statistics of the I-K interval data (mean = 158 ms, standard deviation = 35 ms) in. 29 The parameters for a low HR were set as the default parameter.

The 2ndD-P of the BCG signals created using the default parameters for different HRs were significantly different, as shown in Figure 6. From the spectrum obtained by the Fourier transform, we observed that when the HR was higher, the 2ndD-P had a higher direct current (DC) component, and the lower frequency band (< 0.5 Hz) had a significant energy distribution. However, when the HR is low, the higher frequency band (> 0.5 Hz) has apparent energy, particularly where it corresponds to the HR.

Figure 6. 2ndD-P and its spectrum. Left: Mean HR: 86. Right: Mean HR: 149. At high HRs, there was a significant energy within 0.5 Hz.

Based on these observations, a method is proposed for classifying based on the 2ndD-P. Typically, it is exceedingly rare for an individual with normal physiological conditions to exhibit a HR below 30 bpm (0.5 Hz). 36 In this study, a HR greater than 90 bpm (1.5 Hz) was defined as a high HR. If the HR is low (30 < HR < 90), it has high energy in the frequency band between 0.5 and 1.5 Hz. However, in a high-heart-rate state, the frequency band within 0.5 Hz has an energy distribution and no apparent energy distribution in the frequency band between 0.5 and 1.5 Hz. To highlight the difference in spectral energy distribution between the low and high HRs, we used the ratio of spectral energy integrals, and 0.5 and 1.5 Hz were set as the cutoff. This ratio is defined as follows: (3) Pl=∑f=0f=0.5P(f)

(4) Ph=∑f=0.5f=1.5P(f)

(5) Ratio=PlPh

where f is the frequency, and P(f) is the spectral energy when the frequency is f. Ideally, Ratio should be larger when the HR is higher. The parameters ( NA,WL ) of the proposed method vary with Ratio . The rationale of parameter setting is described as follows:

When the HR is high, the width ( WL ) of the window used to detect the local maximum should be reduced. Otherwise, the method likely misses heartbeats at small adjacent intervals. If the HR increases further, WL should be further reduced to ensure that no heartbeat is missed. Meanwhile, components with high second-order derivative occur frequently, and components with low second-order derivative become fewer. If the NA which equal to 0.24 s is used, the window has many high second-order derivative components, and the low second-order derivative components occupy an extremely small portion of the window. Consequently, 2ndD-P has difficulty distinguishing the high second-order derivative in the signal, which is the location of the IJK wave group. Therefore, NA should also be reduced. The values of the parameters were determined through careful analysis and experimental verification and are listed in Table 1.

Table 1. The relationship between HR level and Ratio, and corresponding parameter settings.

HR Level	Ratio	NA,WL	Note	
Low	Ratio  ≤ 1.7	0.24 s, 0.6 s	Default parameter	
High	1.7 ≤  Ratio  < 4.5	0.16 s, 0.32 s		
Much high	4.5 <  Ratio	0.16 s, 0.08 s		
Note. HR Level: “Low”, HR ≤ 90; “High”, 90 < HR ≤ 140; “Much High”, HR > 140. Ratio is the index calculated from the spectral information of 2ndD-P, as shown in equation (5), used for classifying HR level. NA,WL are parameters of the proposed method, set based on Ratio.

In this study, the first 16 s segment of the signal used for analysis, before creating the template of this signal, were used for Low- and High-Heart-Rate Classification. The segment length was 16 s to ensure a sufficient length of the signal to be Fourier transformed to avoid distortion of the spectral information, whereas a long signal length increased the computational burden.

Template creation

LIJK often deviates from the J-wave. Therefore, to ensure that the pulse segment, SIJK , to contain the IJK wave group, 0.48 s (sampling rate 250 Hz, 120 points) BCG signal segment was divided with LIJK as the center: (6) SIJK={x(LIJK−60),x(LIJK−59),x(LIJK−58),…,x(LIJK+60)}

where x is the BCG signal. After all the SIJK values are obtained in the signal, the template TIJK can be obtained using the following formula: (7) TIJK=1N∑i=1NSIJK

where N is the number of SIJK in a signal. Creating templates based on these means of calculating the averages enables some noise with spectral energy ranges very similar to those of the BCG to be removed and retains the features of the IJK wave group to the greatest extent.

Determining the J-wave in the template

After obtaining TIJK , the peak of the J-wave in TIJK must be determined. By carefully analyzing the BCG signals obtained in this study, we observed that the I–J segment was the most apparent feature of the BCG waveform. We decided to divide TIJK into two types based on the difference in the amplitudes of the I- and J-waves in the template, as shown in Figure 7.

Figure 7. Two types of TIJK . Upper Point (red): The J-wave peak. Lower Point (green): The trough of the I-wave. Left: The amplitude of the J-wave is larger. Right: The I-wave is larger than the J-wave. The horizontal midline represents the threshold value.

In the first type of TIJK , the amplitude of the J-wave is larger, and I-wave is not apparent in comparison. Thus, the positive extreme point with the highest amplitude in the template was determined to be the peak of the J-wave.

In the second type of TIJK , the I-wave is more prominent, and its amplitude is larger. Therefore, the I-wave was detected using a threshold-based method. The threshold was set to half the min negative amplitude in the TIJK . Therefore, the center of the TIJK should be within the range of the IJK wave group. If the I-wave occurs at 0.28 s (sampling rate 250 Hz, point 70), considering the average I-K interval of 158 ms, 29 the K-wave would occur at 0.438 s in TIJK , which is very close to the end of TIJK (0.48 s, point 120). This indicates that the template's center is significantly misaligned with the IJK wave group. This misalignment contradicts the proposed method, suggesting that the template does not belong to the second type and should be discussed separately. Therefore, for the second type of TIJK , the I-wave needs to be constrained to occur within the first 0.28 s of TIJK . In chronological order, the first minimal point that satisfies the above condition is considered to be trough of the I-wave. The first maximum point after the I-wave is determined to be the peak of the J-wave.

Considering the robustness of the algorithm and the complexity and diversity of the data collected in practice, if a TIJK does not conform to the above two types, the point with the maximum second-order derivative value is determined to be the peak of the J-wave.

J-wave identification based on DTW

In this study, DTW was used to implement template matching to obtain point-to-point correspondence 37 between TIJK and each SIJK , as shown in Figure 8.

Figure 8. Point-to-point correspondence between TIJK and SIJK . Upper: TIJK ; Lower: SIJK .

The point-to-point correspondence obtained using DTW is not a single point-to-single point, but rather a point-to-many points. In this case, among all points corresponding to the peak of the J-wave in the TIJK , the point with the largest signal amplitude is determined as the recognition result of the peak of the J-wave.

Reference methods

Reference method 1 (RM1): method based on the short time energy

In ref. 38 , a comparative study of the BCG featuring component extraction algorithms, the method based on the short time energy, which was proposed in ref. 15 , exhibited the highest efficiency and was easy to implement. Therefore, it was selected as a representative of the traditional method for comparison with the proposed method. The steps and specific implementation details of this method are described below.

In this method, a short-time energy transformation is performed on the BCG using a sliding window. For segment xi , in the sliding window, the short-time energy Ei is defined as shown in equation (1). However, the segment size NE was set to 0.32 s (sampling rate 250 Hz, 80 points) to ensure that the IJK wave group was included in segment xi . In the implementation, the same sliding window and step size were used to obtain the STE profile (STE-P), as shown in Figure 9. The local maximum points of STE-P were then determined. In most practical applications, the point of the local maximum does not correspond well with the J-wave of the BCG. Therefore, the positive extreme point in the BCG signal closest to the local maximal point in STE-P was determined as the peak of the J-wave.

Figure 9. BCG and its 2ndD-P and STE-P.

Reference method 2 (RM2): using the STE-P to create the template

For a fair comparison and to better demonstrate the advantages of the proposed method, we propose an improved method RM1, called RM2, in which a template is created using STE-P, and then DTW is used to determine all the peaks of the J-wave, as in the proposed method. Based on this concept, we propose RM2 to clarify the necessity and advantages of 2ndD-P.

Experiments

Experimental design and setup

Two types of experiments were conducted. One of them was the resting-state case collection data for component extraction. A high HR or hypertensive state can lead to variation in the BCG waveform. After sports activities, the HR gradually slows, blood pressure gradually returns to the normal level, and Vt and Vd occur in this process. To verify the validity of the proposed method for Vt and Vd, we performed a second experiment after the participants completed an exercise.

Participants and instruments

This study has been conducted for about one year at Center for Frontier Medical Engineering, Chiba University.

Twenty-four participants were invited to participate in the experiment. Seventeen of them were male and seven were female. The demographics of the participants are shown in Table 2. The age range of the participants was 22–55 years, and their body mass index (BMI) ranged from 17.22 to 29.41. None of the participants had a known history of cardiovascular disease. Before conducting the experiments, participants were informed of all details regarding the experiment, and all participants provided written consent to participate in the experiment.

Table 2. Demographics of the participants.

Characteristic	Mean ± SD	
Age (years)	26.92 ± 6.66	
Weight (kg)	67.66 ± 14.56	
Height (cm)	172.29 ± 8.62	
BMI	22.57 ± 3.49	

For each participant, the BCG and finger pulse (FPu) signal which could assist in labeling the reference value of the J-wave location of the BCG signal were collected.

The BCG sensor (Vital Beats from Future Ink) is a piezoelectric sheet sensor that generates voltage from vibrations or bending. 39 The sensor was placed on a chair with soft cushions to ensure that the vibration of the sensor was not limited. The participants were asked to sit directly on the sensor to collect the BCG signal. The FPu sensor was a Pulse Transducer TN1012/ST (ADInstruments). The participants were asked to wear the sensor with their left thumb. The signals from both sensors were input into Powerlab 2/26 (AD Instruments) to perform analog-to-digital conversion. The experimental setup is shown in Figure 10.

Figure 10. Experimental setup. Left: Top view; Right: Side view.

Resting state

In this experiment, at least two recordings were made for each participant and each recording lasted approximately 3 min. The participants were instructed to remain as still as possible and maintain the resting state during the experiment. To verify that the proposed method could adapt to BCG waveforms in different sitting postures, after completing one recording collection, the participants were asked to adjust their sitting postures, such as moving their hands from the tabletop to their thighs, for a new recording collection. A total of 51 recordings were obtained. For the naming of the dataset, “participant number-data number” format was used, for example, “003-2,” to denote the second recording for the third participant.

For this experiment, not only is the proposed method, RM1 and RM2, proposed method without DTW and proposed method without Low- and High-Heart-Rate Classification are also performed to show the importance of DTW and Low- and High-Heart-Rate Classification.

After exercise

In this experiment, participants were requested to perform ergometer cycling (ERGOMETER 232CXL) to increase their HR. The pedaling power was set according to the individual weight of each participant to ensure that there was no health risk. Participants were instructed to maintain a speed of approximately 90 revolutions per minute for 5 min, but without a target HR to reach. Also, in this experiment, if one participant experienced physical discomfort or believed they could not continue the exercise or measurement, the experiment is to be immediately terminated to ensure safety. After the exercise, the participants were requested to sit and remain as still as possible to collect the BCG and FPu signal for approximately three minutes. As shown in Table 3, 11 participants participated in the experiment.

Table 3. Participants in the experiment “after exercise.”

No.	Gender	Age (years)	Height (cm)	Weight (kg)	BMI	
002	M	25	179	88	27.46	
005	F	27	172	69	23.32	
006	F	24	155	50	20.81	
007	M	23	172	80	27.04	
008	M	27	178	79	24.93	
009	M	24	176	62	20.02	
010	M	24	185	68	19.87	
017	M	27	184	87	25.70	
018	F	27	160	46.4	18.13	
021	F	24	165	52	19.10	
022	M	25	176	82	26.47	

For these recordings, to distinguish them from the resting state recording, we named them in the format “participant number-H”. To prove that the proposed method can be adapted to high HR recordings compared with the reference method, and to show the importance of Low- and High-Heart-Rate Classification. Proposed method, RM1 and RM2, proposed method without Low- and High-Heart-Rate Classification are performed. Moreover, for a fair comparison, the parameters of RM2 were manually adjusted to accommodate a high-heart-rate state as well. Here, NE was adjusted to 0.24 s, and to recording 008-H whose HR is highest, the length of the window for finding the local maxima was set to 0.08 s.

Filtering strategy

In this study, the preprocessing filter is based on the wavelet decomposition and reconstruction method. Some studies have shown that filtering strategies may have a key impact on the feature extraction of physiological signals. 33 Therefore, it is also necessary to compare different filtering strategies. Moreover, it is important to make sure that the proposed method does not rely on pre-processing techniques.

In this experiment, three filtering techniques are compared: Wavelet decomposition and reconstruction: the specific implementation details have been described in 2.2.

Widely used IIR Butterworth filter: A 4-order high-pass filter with a cutoff frequency of 1 Hz and a 4-order low-pass filter with a cutoff frequency of 15.6 Hz were used to remove the baseline drift and high-frequency noise.

Elliptic IIR filter which has effective performance in processing of PPG 33 : A 4-order bandpass elliptic IIR filter is used with a lower and upper cutoff frequency of 1 Hz 15.6 Hz.

Reference signal and manual labeling

In this study, the FPu signal was used as a reference signal to assist with the labeling of the J-waves. As shown in Figure 11, the peak of the FPu signal was generally extremely close to and occurred before the peak of the J-wave of the BCG collected in the sitting state.

Figure 11. Preprocessed BCG and FPu. The peak of the J-wave can be determined intuitively through the FPu.

Based on this observation, all peaks of the J-wave were manually labeled to build a dataset of reference values for measuring the effectiveness of the method. In the actual recordings, we observed that a few segments of the BCG signal were severely disturbed by body motion and the J-wave could not be determined. In this case, the peaks of the J-wave were labeled according to the interval between the other peaks of FPu signals and J-waves in the same recording.

And with this reference value dataset, instantaneous HR of each recording was calculated. In addition, a 16-s width sliding window with a 15-s step length was used to calculate the mean HR and Ratio within the window to study the relationship between them.

Performance metrics

Sensitivity and positive predictivity

Using the reference value dataset, the performance of the method was evaluated using precise metrics. Metrics widely used in the field of ECG analysis, 40 sensitivity (Se) and positive predictivity (+P) were employed to evaluate the performance of the method. Therefore, the numbers of false positives ( FP ), true positives ( TP ), and false negatives ( FN ) for each set of data must be counted. The formulae for calculating Se and +P are as follows: (13) Se=TPTP+FN

(14) +P=TPTP+FP

Standard deviation of normal-normal intervals (SDNN)

Heart rate variability (HRV) is of great significance for the monitoring of daily life.41,42 Therefore, for signals related to cardiac activity, the calculation of HRV parameters is very important. Ref. 43 proved that in the calculation of the SDNN which is an important HRV parameter, the BCG and pulse wave signals measured by the pressure sensor were highly consistent with the ECG. Therefore, in this study, the SDNN was selected as the metric, and FPu was used as standard. The SDNN is calculated as follows: 44 (15) SDNN=1N−1∑i=1N−1(JJi−JJ¯)2

where N is the number of heartbeat intervals. JJ is the interval between the detected peaks of the J-wave and the FPu signal.

Results

Low- and high-heart-rate classification

Ratio is calculated using the default parameter ( NA=0.24s,60points ). As described in 3.2, a total of 780 data were obtained by using the sliding window and the results are shown in Figure 12. A clear positive relationship was observed between Ratio and HR (correlation value: 0.80). This result demonstrates that the proposed method is effective in Low- and High-Heart-Rate Classification. The experimental results also prove the role of the classification.

Figure 12. Ratio and HR. Blue dash-dot line: HR = 90, Ratio  = 1.7; Red dashed line: HR = 140, Ratio  = 4.5.

As for comparison, instead of 2ndD-P, the HR-Ratio relationship in terms of STE-P is shown in the lower part of Figure 12. Although the correlation value (0.75) was similar to result obtained using 2ndD-P, it is difficult to set thresholds to distinguish between low- and high-heart-rate.

J-wave extraction

Resting state

The results of the methods for 51 recordings of the resting-state experiments are listed in Table 4. As shown in Table 4, the proposed method had the highest Se and +P . RM1 performed the worst. The proposed method and RM2, which used both the template and DTW, yielded good results. It is worth noting that the proposed method has apparent advantages over RM2 on recording 003-1 and recording 024-2. Excluding these two recordings manually, the results were as follows: (Proposed: Se : 98.37%, +P : 98.71%; RM2: Se : 98.24%, +P : 97.99%). The proposed approach has several advantages.

Table 4. J-wave identification results.

	Proposed	Proposed without low- and high-heart-rate classification	Proposed without DTW	RM1	RM2	
Recording	Se (%)	+P (%)	Se (%)	+P (%)	Se (%)	+P (%)	Se (%)	+P (%)	Se (%)	+P (%)	
001-1	100.00	100.00	100.00	100.00	100.00	100.00	100.00	99.32	100.00	99.32	
…	…	…	…	…			…	…	…	…	
003-1	99.27	99.51	99.27	99.51	97.83	98.30	83.29	84.94	88.38	90.12	
010-1	100.00	100.00	100.00	100.00	100.00	100.00	100.00	100.00	100.00	100.00	
022-1	99.42	99.71	98.83	99.71	99.42	99.42	99.71	99.12	99.13	99.42	
022-2	98.22	99.10	95.55	99.38	97.35	97.63	99.38	98.80	98.22	99.40	
024-2	92.27	94.86	92.27	94.86	51.80	53.74	42.53	42.34	66.06	65.77	
…	…	…	…	…	…	…	…	…	…	…	
Mean	98.29	98.64	98.22	98.65	90.95	91.21	86.78	86.57	97.42	97.20	
SD	2.41	2.09	2.43	2.09	16.44	16.34	19.01	18.84	5.41	5.32	
Note. Bold: 1.7 ≤  Ratio  < 4.5. Only the results for a subset of the resting state recordings are illustrated here; full results can be found in the supplementary materials Section 2.

As shown, if DTW is not used, compared with the full implementation of the proposed method, the accuracy is significantly reduced. Moreover, the Low- and High-Heart-Rate Classification didn’t make difference for the resting state BCGs.

Because the waveform of each heartbeat of BCG collected in the same sitting position in the resting state has high consistency, recordings can be further divided according to the morphology of the template. If the J-wave amplitude was smaller than that of the other components in the template, the recording was considered a J-wave non-prominent recording, as shown in Figure 13. For these recordings, the results are shown in Table 5.

Figure 13. The J-wave amplitude was smaller than that of the other components. The point on the waveform: The J-wave peak.

Table 5. Results for the recordings with non-prominent J-wave.

	Proposed	Proposed without DTW	RM1	
Recording	Se (%)	+P (%)	Se (%)	+P (%)	Se (%)	+P (%)	
001-1	100.00	100.00	100.00	100.00	100.00	99.32	
001-2	97.96	97.96	93.20	93.20	92.57	93.20	
001-3	91.50	91.50	81.05	81.05	81.05	81.05	
003-1	99.27	99.51	97.83	98.30	83.29	84.94	
010-2	100.00	100.00	64.76	65.04	37.00	37.00	
014-1	96.19	96.60	93.22	93.62	95.36	95.36	
014-2	98.03	98.42	96.85	97.23	96.06	96.44	
015-1	97.74	98.48	66.04	66.54	24.44	24.34	
015-2	87.96	88.60	60.22	60.66	36.50	36.63	
024-1	98.23	98.23	50.44	50.44	36.73	36.73	
Mean	96.69	96.93	80.36	80.61	68.30	68.50	
SD	3.73	3.63	17.44	17.44	28.98	29.08	

As can be seen from Table 5, the proposed method has apparent advantages for BCG waveforms with non-prominent J-waves.

Meanwhile, the advantages brought by the combined use of DTW cannot be ignored. As shown in Figure. 14, the accuracy of RM2 is significantly improved compared with that of RM1. However, due to the Va, RM2 is less effective than the proposed method in some cases.

Figure 14. Results of proposed method and reference methods. FP FN TP Left: For recording 015-1, compared with RM1, the accuracy of RM2 qualitatively improved. Right: For recording 003-1, if the J-wave is not prominent, RM2 may cause many errors.

The different details of the use of STE-P and 2ndD-P are shown in Figure 15. The local maximum point of the STE-P deviated more from the true J-wave position than that of the 2ndD-P. The local maximum point of the 2ndD-P also showed some deviation from the J-wave. However, owing to the use of DTW, the point-to-point relationship between the template and signal segment allows the peaks of the J-wave to be correctly labeled.

Figure 15. Details of the BCG segmentation results. Left: Labeling result of a BCG segment; Right: Details of DTW results for the fifth heartbeat of the left signal segment. Because of the variability of BCG, 2ndD-P is more suitable for determining the approximate position of the J-waves than STE-P. Owing to the use of DTW, the peak of J-wave can be correctly marked even if LIJK deviate from the real J-wave.

About the SDNN calculation, the proposed method was closer to the results calculated using the FPu signal than reference methods. As shown in Figure 16, the proposed method exhibited a significantly better correlation.

Figure 16. Bland–Altman plot between FPu-SDNN and BCG-SDNN obtained by (a) the proposed method, (b) RM1, and (c) RM2. SDNN calculation using the positions of the J-wave peak obtained using the proposed method has advantages over reference methods.

Additionally, the three evaluative metrics associated with the Bland–Altman plot (the mean bias, standard deviation, and limits of agreement) are also presented, as shown in Table 6. The proposed method exhibited distinct advantages compared with reference methods.

Table 6. Evaluation metrics of Bland–Altman plot.

	Proposed	RM1	RM2	
Bias	−0.021	−0.053	0.029	
SD	0.027	0.041	0.036	
Limits of Agreement	−0.073, 0.032	−0.133, 0.027	−0.099, 0.041	

After exercise

The results for 11 recordings collected after exercise are presented in Table 7. Compared with the reference methods, the proposed approach has obvious advantages. On the other hand, the results also reflect the importance of parameter adjustment based on the Low- and High-Heart-Rate Classification.

Table 7. Results for “after exercise.”

Method	Proposed	Proposed without Low- and high-heart-rate Classification	RM1	RM2	RM2 with adjustment	
Recording	Mean HR
(Max, Min)	Se (%)	+P (%)	Se (%)	+P (%)	Se (%)	+P (%)	Se (%)	+P (%)	Se (%)	+P (%)	
002-H	109
(126, 98)	98.28	100.00	91.07	100.00	95.88	100.00	95.88	100.00	97.94	100.00	
005-H	107
(123, 96)	98.36	100.00	90.13	100.00	77.45	84.95	91.45	99.64	96.05	99.66	
006-H	96
(116, 84)	95.36	97.16	93.50	98.05	14.24	14.70	94.74	97.76	95.05	96.85	
007-H	98
(110, 92)	100.00	100.00	99.25	100.00	97.74	97.74	100.00	100.00	100.00	100.00	
008-H	130
(149, 120)	97.27	98.74	53.69	100.00	66.16	83.87	78.88	100.00	96.74	99.48	
009-H	121
(136, 107)	94.44	98.74	74.15	99.35	75.60	86.23	87.44	99.72	93.72	99.74	
010-H	106
(115, 101)	100.00	100.00	99.30	100.00	100.00	100.00	100.00	100.00	100.00	100.00	
017-H	115
(137, 87)	90.45	96.41	71.63	96.96	82.02	95.11	83.43	96.74	89.61	96.67	
018-H	101
(109, 96)	100.00	100.00	98.77	100.00	99.69	99.69	100.00	100.00	100.00	100.00	
021-H	105
(112, 101)	99.10	99.10	96.69	99.07	29.52	29.52	99.10	99.10	99.10	99.10	
022-H	121
(126, 119)	95.24	99.01	73.97	99.57	92.06	98.31	92.38	98.64	95.87	98.69	
Mean		97.14	99.01	85.64	99.36	75.49	80.92	92.97	99.24	96.73	99.11	
SD		2.85	1.18	14.26	0.96	27.58	28.52	6.85	1.05	3.06	1.18	
Note. “Mean HR” represents the mean HR across the recordings. “Max” and “Min” represent, respectively, the maximum and minimum heart rates obtained for this recording using the sliding window method described in section Reference signal and manual labeling.

HR decreases over time and the change of Ratio is also consistent with the change of HR, as shown in Figure 17.

Figure 17. Ratio and instantaneous HR for recording 017-H. Red points: The mean HR; blue diamonds: The Ratio value. Both were calculated by using a 16 s sliding window.

Compared with low-heart-rate signals, FNs were more likely to occur in the detection results of high-heart-rate signals, as shown in Figure 18.

Figure 18. Segment from the recording 008-H. There are two FNs. The periodicity of the BCG is clear.

Filtering strategy

After preprocessing with three different filtering techniques, the proposed method is used to extract J-wave from 51 recordings of resting state and 11 recordings collected after exercise. The results (mean of the metrics values) are shown in Table 8.

It should be noted that the elliptic filter can cause severe distortion in the waveforms of some recordings. If these distorted recordings are excluded, the results are as follows: (resting state: Se: 98.49%, +P: 98.87%; After Exercise: Se: 97.03%, +P: 98.97%). The results indicate that different filtering techniques do not have a significant impact on the proposed method.

Discussion

To solve variability of BCG

The experiments showed that the proposed method performs satisfactorily in J-wave extraction and HR variability parameter calculations, as shown in Tables 4–7. This is because the waveform variability of the BCG is fully considered in this method. 2ndD-P, the template and DTW, and Low- and High-Heart-Rate Classification, were used to address the three types of variability (Va, Vt, and Vd) in the BCG.

2ndD-P is the core of the proposed method. It is the basis for creating templates for template matching and implementing Low- and High-Heart-Rate Classification. Using the 2ndD-P, the positions of the IJK wave groups could be found more accurately. Specially, as shown in Tables 4 and 5, even without the DTW, the 2ndD-P based method has obvious advantages over RM1. These results show that the proposed method can deal with Va using the 2ndD-P, particularly while the J-wave is not the most prominent component. Although the J-wave is not prominent, it and its surrounding components still have higher second-order derivative values in one cardiac cycle. Therefore, the 2ndD-P proposed in this study can better capture the vibrational nature of BCG.

Another key factor that addresses waveform variability is the combined use of the template and DTW. Owing to variability, determining apparent commonalities in BCG among different participants as the basis for feature extraction is difficult. By creating a participant- and moment-dependent template and then using the template to match the feature segments, detecting all the peaks of the J-wave can be effectively reduced to determining the peak in the template. The validity of this concept was demonstrated through experimental results. As shown in Table. 4 and Figure 14, RM2 using the template method had clear advantages over RM1.

After determining the peak of the J-wave in the template, using this information to identify each peak of the J-wave should be considered, which must use DTW owing to Vt. As shown in Table. 4, the proposed method exhibited a significant improvement compared with the proposed method without DTW. This is because the occurrence times of the BCG components is not uniform. Comparing RM1 with the proposed method without DTW, 2ndD-P has apparent advantages in finding the segment containing the IJK wave group. However, because of Vt, this segment was not always consistent with the J-wave, as shown in Figure 15, which requires DTW to adapt to Vt.

In this study, Vd primarily occurred in recordings of a high HR. As shown in Figures 2 and 17, for the signal of participant 008 in the resting state, a clear interval was observed between each heartbeat. For a signal with a high HR, the interval was very short; thus, the components with high second-order derivative value appeared at shorter intervals, and some other components of the BCG waveform were aliased and could not be distinguished. To solve Vd, this paper proposes a Low- and High-Heart-Rate Classification. For the resting-state recordings, the presence or absence of Low- and high-heart-rate classification did not exhibit a clear difference in the results. As shown in Table 4, in the resting-state recordings, only the “010-1”, “022-1”, and “022-2” recordings had an average HR higher than 90 and a Ratio value greater than 1.7. The effect was apparent in 022-1 and 022-2, and the sensitivity significantly improved. In terms of high-heart-rate recordings, Low- and High-Heart-Rate Classification can greatly improve the sensitivity, as shown in Table 7. If the parameters are not modified, as shown in Figure 18, many J-waves are missed. Further supporting evidence is that after manually modifying the parameters of RM2, the accuracy of RM2 for high-heart-rate recordings is also greatly improved. The internal reason for the success of the Low- and High-Heart-Rate Classification is that the amplitude of the BCG with a high HR showed a clear periodic variation, which can be found in the spectrum of the 2ndD-P, as shown in Figure 6. Additionally, for recordings exhibiting a gradual decline in HR, as illustrated in Fig. 17, satisfactory results are achieved by calculating the Ratio using only the signal from the first 16 s, without further updating the Ratio . This is due to Ratio and its corresponding parameter design incorporate a certain level of tolerance (Section 3: The Tolerance of the Parameters in the supplementary materials).

It is worth noting, as shown in Table 8, that variations in the filtering strategy do not significantly affect the performance of the proposed method. The results from three different filtering techniques (though the results from the elliptic filter require the exclusion of recordings with obvious distortion) show no significant differences. This indicates that the proposed method does not rely on the preprocessing filtering techniques. And through the ablation studies on the parameter WL (Section 4: Ablation Studies in the supplementary material), the proposed method does not exhibit significant changes with regard to slight modification of the values of the parameter.

Table 8. The result of different filtering strategies.

	Wavelet	Butterworth	Elliptic IIR	
	Se (%)	+P (%)	Se (%)	+P (%)	Se (%)	+P (%)	
Resting state	98.29	98.64	98.26	98.61	94.64	95.01	
After exercise	97.11	99.01	97.26	99.14	70.65	72.06	

In summary, the best performance was achieved by proposed method whose components work cohesively with each other to solve three types of variability. 2ndD-P plays the most important role in this work. Based on the frequency spectrum of 2ndD-P, the parameters can be adjusted according to the heart rate to obtain the locations of IJK wave group more accurately. These locations are the basis for subsequent template creation and template matching with DTW. By solving variability, the proposed method has obvious advantages over the reference methods. Compared with the NN based methods17–19 reported in recent years, the proposed method is simple and reliable, and fully explainable.

Contribution and limitations

Contribution

In this work, for the first time, the three types of variability of BCG are addressed systematically, by a 2ndD-P based method, which integrates a template matching, DTW, and Low- and High-Heart-Rate Classification based on 2ndD-P frequency spectrum. Moreover, this is the first time for the 2ndD-P to be used to extract feature components for BCG. It can better capture the nature of the BCG as a vibration signal and has great help in dealing with the Va of the BCG. Especially in the case of BCG with non-prominent J-wave, 2ndD-P can help capture the feature efficiently. The low- and high- heart-rate classification based on the spectral information of 2ndD-P contributes to the BCG in different heart rate. Based on this classification, the method can be adapted to high-heart-rate recordings, which have not been carefully analyzed in most related studies. In this study, BCG is one step closer to its application in daily life.

Limitations

The largest limitation of this study is that all data were measured in a state in which the participants were required to remain as still as possible. Data on body-motion artifacts have not been fully discussed. On the other hand, for template-based methods, the signal period length used to create the template needs to be discussed. However, in this study, there was no significant difference in the results of different template creating methods. This is most likely because the participants were instructed to remain still during the experiment. When body movements, especially changes in participant posture, are included, template updates have to be considered.

And the data in this study included only healthy individuals and did not explore cases of cardiovascular diseases. Theoretically, cardiovascular diseases are likely to have a strong impact on BCG signal. 45 For example, mitral valve membrane insufficiency may cause blood flow vibrations in the heart to be less consistent than those observed in a healthy person. 46 It is likely that there was significant Vd in these patients. Because of Vd, the template of BCG may also be more types than the two presented here. Therefore, BCG signals in people with cardiovascular diseases must be carefully studied.

Finally, since the sensor used in this study is a piezoelectric sheet sensor, the proposed method may be suitable only for the vibration displacement BCG collected by this sensor. The applicability of the proposed method for BCG collected by other kinds of sensors needs to be further discussed.

Future work

The current data were collected when the participants were required to maintain a resting state that deviated from their actual life situation. Signals with body motion or other artifacts must be considered in future studies. We proposed a BCG signal quality evaluation method that can remove the poor-quality BCG segments that are heavily affected by body movement. 47 By combining this method, the proposed method can be applied to a signal with body-motion interference. Accordingly, we plan to conduct additional tests in the future. In addition, to further study how to solve Vd and whether the proposed method can be applied to the elderly or people with cardiovascular diseases, BCG signals from elderly and non-healthy individuals must be collected and tested.

Conclusion

This paper presents a method of extracting the BCG J-wave using 2ndD-P, a template, DTW, and Low- and High-Heart-Rate Classification, with full consideration of BCG waveform variability. Fifty-six recordings from 24 participants were tested, and the method performed well with satisfactory Se and +P . Compared with the direct use of the original BCG waveform, the use of the second-order derivative waveform can be useful in solving Va and determining the location of the main features of the BCG. Using a template and DTW, the proposed method can solve Vt well; therefore, J-waves can be determined more accurately. Through Low- and High-Heart-Rate Classification based on the spectrum information of the 2ndD-P, adaptive parameters are also achieved to solve Vd to some extent, such that the proposed method also has a good effect on high HR recordings.

Supplemental Material

sj-docx-1-dhj-10.1177_20552076241277746 - Supplemental material for Solving variability: Accurately extracting feature components from ballistocardiograms

Supplemental material, sj-docx-1-dhj-10.1177_20552076241277746 for Solving variability: Accurately extracting feature components from ballistocardiograms by Tianyi Yang, Haihang Yuan, Junqi Yang, Zhongchao Zhou, Masayuki Abe, Yoshitake Nakayama, Shao Ying Huang and Wenwei Yu in DIGITAL HEALTH

Acknowledgements

We would like to thank all the participants in this research.

Contributorship: TY was responsible for the methodology, management of the experimental data of the study, and wrote the manuscript. HY was responsible for conducting the experiment and participated in the management of the experimental data. JY and ZZ contributed to the experimental design and method development. MA, YN, and SH contributed to the discussion and review of the paper. WY was responsible for the supervision of the study and for the review and editing of the manuscript.

All authors reviewed and edited the manuscript and approved the final version of the manuscript.

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Ethics approval and consent to participate: This research was conducted in accordance with the Declaration of Helsinki. And the Joint Institutional Review Board of the Graduate School of Science and Engineering at the Design Research Institute, Center of Frontier Medical Engineering (R5-01) approved this research.

Before starting the experiment in this research, participants were informed of the specific details. After obtaining informed consent from the participants, the experiment was conducted.

Funding: The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by JST SPRING, Grant Number JPMJSP2109.

Guarantor: WY is the guarantor and corresponding author of this research.

ORCID iD: Tianyi Yang https://orcid.org/0000-0002-0711-919X

Supplemental material: Supplemental material for this article is available online.
==== Refs
References

1 Giovangrandi L Inan OT Wiard RM , et al. Ballistocardiography – a method worth revisiting. Conf Proc Annu Int Conf IEEE Eng Med Biol Soc IEEE Eng Med Biol Soc Annu Conf 2011; 2011 : 4279–4282.
2 Pinheiro E Postolache O Girão P . Theory and developments in an unobtrusive cardiovascular system representation: ballistocardiography. Open Biomed Eng J 2010; 4 : 201–216.21673836
3 Bruser C Diesel J Zink MDH , et al. Automatic detection of atrial fibrillation in cardiac vibration signals. IEEE J Biomed Health Inform 2013; 17 : 162–171.23086532
4 Ozcelik STA Uyanık H Deniz E , et al. Automated hypertension detection using ConvMixer and spectrogram techniques with ballistocardiograph signals. Diagnostics 2023; 13 : 82.
5 Feng S Zhang H Bao A , et al. Diagnosis of heart failure using high quality ballistocardiography and respiratory effort signals: a pilot study. In: 2022 14th Biomedical Engineering International Conference (BMEiCON 2022), Held 10–13 November 2022, Songkhla, Thailand, pp.1–5. 10.1109/BMEiCON56653.2022.10012098
6 Qi P Gong S Jiang N , et al. Mattress-Based non-influencing sleep apnea monitoring system. Sensors 2023; 23 : 3675.37050735
7 Inan OT Migeotte P-F Park K-S , et al. Ballistocardiography and seismocardiography: a review of recent advances. IEEE J Biomed Health Inform 2015; 19 : 1414–1427.25312966
8 Clausen T Jura T Jähne-Raden N , et al. A precise, parallel and scalable measurement system for ballistocardiographic research. Smart Health 2021; 19 : 100169.
9 Xing X Li H Chen Q , et al. Blood pressure monitoring with piezoelectric bed sensor systems. Biomed Signal Process Control 2024; 87 : 105479.
10 Su BY Enayati M Ho KC , et al. Monitoring the relative blood pressure using a hydraulic bed sensor system. IEEE Trans Biomed Eng 2019; 66 : 740–748.30010544
11 Feng J Huang W Jiang J , et al. Non-invasive monitoring of cardiac function through Ballistocardiogram: an algorithm integrating short-time Fourier transform and ensemble empirical mode decomposition. Front Physiol 2023; 14 . https://www.frontiersin.org/articles/10.3389/fphys.2023.1201722 (accessed November 2, 2023).
12 Pino EJ Chávez JAP Aqueveque P . BCG algorithm for unobtrusive heart rate monitoring. In: 2017 IEEE Healthcare Innovations and Point of Care Technologies (HI-POCT, 2017), pp.180–183. 10.1109/HIC.2017.8227614 .
13 Jin J Wang X Li S , et al. A novel heart rate detection algorithm in ballistocardiogram based on wavelet transform. In: 2009 Second International Workshop on Knowledge Discovery and Data Mining, 2009: pp.76–79. 10.1109/WKDD.2009.98 .
14 Pino EJ Chávez JAP Aqueveque P . Noninvasive ambulatory measurement system of cardiac activity. In: 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC, 2015), pp.7622–7625. 10.1109/EMBC.2015.7320157 .
15 Lydon K Su BY Rosales L , et al. Robust heartbeat detection from in-home ballistocardiogram signals of older adults using a bed sensor. In: 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC, 2015). pp.7175–7179. 10.1109/EMBC.2015.7320047 .
16 Albukhari A Lima F Mescheder U . Bed-embedded heart and respiration rates detection by longitudinal ballistocardiography and pattern recognition. Sensors 2019; 19 : 1451.30934577
17 Mai Y Chen Z Yu B , et al. Non-Contact heartbeat detection based on ballistocardiogram using UNet and bidirectional long short-term memory. IEEE J Biomed Health Inform 2022; 26 : 3720–3730.35333727
18 Huang Y Jin T Sun C , et al. Efficient J peak detection from ballistocardiogram using Lightweight convolutional neural network. In: 2021 43rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC, 2021), pp.269–272. 10.1109/EMBC46164.2021.9630255 .
19 Zhou T Men S Liang J , et al. 1d u-net++: an effective method for ballistocardiogram j-peak detection. J Mech Med Biol 2021; 21 : 2140058.
20 Javaid AQ Töreyin H Inan OT . Quantifying the effects of blood pressure changes on ballistocardiogram signals. In: 2017 25th Signal Processing and Communications Applications Conference (SIU, 2017), pp.1–4. 10.1109/SIU.2017.7960423 .
21 Etemadi M T O . Inan, wearable ballistocardiogram and seismocardiogram systems for health and performance. J Appl Physiol 2018; 124 : 452–461.28798198
22 Cathelain G Rivet B Achard S , et al. Dynamic Time warping for heartbeat detection in ballistocardiography. In: 2019 Computing in Cardiology (CinC, 2019), pp.1–4. 10.22489/CinC.2019.145 .
23 Suliman A Mowla MR Alivar A , et al. Effects of ballistocardiogram peak detection jitters on the quality of heart rate variability features: a simulation-based case study in the context of sleep staging. Sensors 2023; 23 : 2693.36904896
24 Jiang F Hong C Cheng T , et al. Attention-based multi-scale features fusion for unobtrusive atrial fibrillation detection using ballistocardiogram signal. Biomed Eng Online 2021; 20 : 12.33509212
25 Alametsä J Palomäki A Viik J . Short and longer term repeatability of ballistocardiography in a sitting position with EMFi sensor. Med Biol Eng Comput 2011; 49 : 881–889.21336678
26 Xie Q Wang M Zhao Y , et al. A personalized beat-to-beat heart rate detection system from ballistocardiogram for smart home applications. IEEE Trans Biomed Circuits Syst 2019; 13 : 1593–1602.31804943
27 Jiao C Lyons P Zare A , et al. Heart beat characterization from ballistocardiogram signals using extended functions of multiple instances. In: 2016 38th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC, 2016), pp.756–760. 10.1109/EMBC.2016.7590812 .
28 Park J Seok HS Kim S-S , et al. Photoplethysmogram Analysis and Applications: An Integrative Review. Front Physiol 2022; 12 , https://www.frontiersin.org/journals/physiology/articles/10.3389/fphys.2021.808451 (accessed February 23, 2024).
29 Kim C-S Ober SL McMurtry MS , et al. Ballistocardiogram: mechanism and potential for unobtrusive cardiovascular health monitoring. Sci Rep 2016; 6 : 31297.27503664
30 Müller M (eds). Dynamic time warping in information retrieval music motion. Berlin, Heidelberg: Springer, 2007, pp.69–84. 10.1007/978-3-540-74048-3_4 .
31 Shen J Bao S-D Yang L-C , et al. The PLR-DTW method for ECG based biometric identification. In: 2011 Annual International Conference of the IEEE Engineering in Medicine and Biology Society, pp.5248–5251. 10.1109/IEMBS.2011.6091298 22255521
32 Ghassemi M Lehman L Snoek J , et al. Global optimization approaches for parameter tuning in biomedical signal processing: a focus on multi-scale entropy. Comput Cardiol 2014; 2014 : 993–996. https://ieeexplore.ieee.org/document/7043212 (accessed November 2, 2023).
33 Mejía-Mejía E Kyriacou PA . Effects of noise and filtering strategies on the extraction of pulse rate variability from photoplethysmograms. Biomed Signal Process Control 2023; 80 : 104291.
34 Xie Q Li Y Wang G , et al. An unobtrusive system for heart rate monitoring based on ballistocardiogram using Hilbert transform and Viterbi decoding. IEEE J Emerg Sel Top Circuits Syst 2019; 9 : 635–644.
35 Savitzky A Golay MJE . Smoothing and differentiation of data by simplified least squares procedures. Anal Chem 1964; 36 : 1627–1639.
36 Opthof T . The normal range and determinants of the intrinsic heart rate in man. Cardiovasc Res 2000; 45 : 177–184.10841627
37 Keogh E Ratanamahatana CA . Exact indexing of dynamic time warping. Knowl Inf Syst 2005; 7 : 358–386.
38 Suliman A Carlson C Ade CJ , et al. Performance comparison for ballistocardiogram peak detection methods. IEEE Access 2019; 7 : 53945–53955.
39 PRODUCTS | 株式会社Future Ink, 株式会社Future Ink (2019). http://www.futureink.co.jp/eng/products/ (accessed December 30, 2023).
40 Arzeno NM Deng Z-D Poon C-S . Analysis of first-derivative based QRS detection algorithms. IEEE Trans Biomed Eng 2008; 55 : 478–484.18269982
41 Malik M Camm AJ . Heart rate variability. Clin Cardiol 1990; 13 : 570–576.2204508
42 Rajendra Acharya U Paul Joseph K Kannathal N , et al. Heart rate variability: a review. Med Biol Eng Comput 2006; 44 : 1031–1051.17111118
43 Cui H Wang Z Yu B , et al. Statistical analysis of the consistency of HRV analysis using BCG or pulse wave signals. Sensors 2022; 22 : 2423.35336592
44 Wang H-M Huang S-C . SDNN/RMSSD as a surrogate for LF/HF: a revised investigation. Model Simul Eng 2012; 2012 : 16:16.
45 Zaid M Sala L Ivey JR , et al. Mechanism-Driven Modeling to Aid Non-invasive Monitoring of Cardiac Function via Ballistocardiography. Front Med Technol 2022; 4 , https://www.frontiersin.org/articles/10.3389/fmedt.2022.788264 (accessed November 2, 2023).
46 Nishimura RA Vahanian A Eleid MF , et al. Mitral valve disease—current management and future challenges. Lancet 2016; 387 : 1324–1334.27025438
47 Yang T Yuan H Zhou Z , et al. Evaluation of ballistocardiogram based on short-time energy and dynamic time warping. In: Lee S-G An J Chong NY Strand M Kim JH (eds) Intell. Auton. Syst. 18. Cham, Switzerland: Springer Nature, 2024, pp.333–346. 10.1007/978-3-031-44981-9_28 .
