Next Article in Journal
EZ-SEP: Extended Z-SEP Routing Protocol with Hierarchical Clustering Approach for Wireless Heterogeneous Sensor Network
Next Article in Special Issue
Efficient Anomaly Detection for Smart Hospital IoT Systems
Previous Article in Journal
Intelligent Decision-Making of Scheduling for Dynamic Permutation Flowshop via Deep Reinforcement Learning
Previous Article in Special Issue
A Framework for Detecting and Analyzing Behavior Changes of Elderly People over Time Using Learning Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Spectral-Based Approach for BCG Signal Content Classification †

1
Higher School of Communication of Tunis, Research Lab. COSIM, Carthage University, Tunis 2088, Tunisia
2
Laboratoire de Recherche en Imagerie et en Orthopédie, CRCHUM, Montreal, QC H2X 0A9, Canada
3
Department of Computer Science, Sherbrooke University, Sherbrooke, QC J1K 2R1, Canada
4
LICEF Institute, TELUQ University, Montreal, QC H2S 3L4, Canada
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in: Ben Nasr, M.C.; Ben Jebara, S.; Otis, S.; Abdulrazak, B.; Mezghani, N. Respiratory Activity Classification Based on Ballistocardiogram Analysis. In Proceedings of the 18th International Conference on Smart Living and Public Health: The Digital Technologies Impact on Public Health in Developed and Developing Countries, ICOST 2020, Hammamet, Tunisia, 24–26 June 2020.
Sensors 2021, 21(3), 1020; https://doi.org/10.3390/s21031020
Submission received: 9 December 2020 / Revised: 27 January 2021 / Accepted: 28 January 2021 / Published: 2 February 2021

Abstract

:
This paper has two objectives: the first is to generate two binary flags to indicate useful frames permitting the measurement of cardiac and respiratory rates from Ballistocardiogram (BCG) signals—in fact, human body activities during measurements can disturb the BCG signal content, leading to difficulties in vital sign measurement; the second objective is to achieve refined BCG signal segmentation according to these activities. The proposed framework makes use of two approaches: an unsupervised classification based on the Gaussian Mixture Model (GMM) and a supervised classification based on K-Nearest Neighbors (KNN). Both of these approaches consider two spectral features, namely the Spectral Flatness Measure (SFM) and Spectral Centroid (SC), determined during the feature extraction step. Unsupervised classification is used to explore the content of the BCG signals, justifying the existence of different classes and permitting the definition of useful hyper-parameters for effective segmentation. In contrast, the considered supervised classification approach aims to determine if the BCG signal content allows the measurement of the heart rate (HR) and the respiratory rate (RR) or not. Furthermore, two levels of supervised classification are used to classify human-body activities into many realistic classes from the BCG signal (e.g., coughing, holding breath, air expiration, movement, et al.). The first one considers frame-by-frame classification, while the second one, aiming to boost the segmentation performance, transforms the frame-by-frame SFM and SC features into temporal series which track the temporal variation of the measures of the BCG signal. The proposed approach constitutes a novelty in this field and represents a powerful method to segment BCG signals according to human body activities, resulting in an accuracy of 94.6%.

1. Introduction

Amidst the expansion in the field of the Internet of Medical Things (IoMT), a significant shift in the paradigm of tele-healthcare has been recorded. In particular, the monitoring of vital signs is gaining attention and has established itself as no longer a commodity but a necessity in this field. Four principal vital signs serve to establish an early warning score, namely the body temperature, blood pressure, heart rate (RR) and respiratory rate (RR). This score is the gold standard when it comes to quantifying the degree of illness of patients [1,2].
The Ballistocardiogram signal (BCG) is a noninvasive, unobstructive measure of the ballistic force generated by the circulation of the blood in the body. Several sensors have been designed to extract the BCG signal using for instance cameras [3], accelerometers [4], weight scales [5,6], ear wearables [7] and piezoelectric sensors [8]. Recently, an interesting microbend Fiber Optic Sensor (FOS) has been placed under mattresses; this noninvasive and unobstructive method is boosting the use of the BCG in telemedicine [9] since it provides more flexibility and eliminates the inconveniences present in invasive or minimally invasive methods [10].
Several works have studied the extraction of vital signs from BCG signals. For instance, in [11], the authors developed different methods for the extraction of the RR such as the Harmonic and Noise (HNM) model and Wavelet Transform (WT), yielding promising results. In [12], the authors proposed a machine learning-based approach for the estimation of the beat-to-beat HR from a highly volatile BCG signal. Features describing the BCG waveform were extracted, and a modified k-means algorithm was used as an inference tool. Other works, such as in [13], implemented a template matching approach to detect the J-peaks of the BCG waveform. In [14], two different methods were implemented for the detection of the HR, namely the Maximum Overlap Discrete Wavelet Transform (MODWT) and Complete Ensemble Empirical Mode Decomposition with Adaptive Noise (CEEMDMAN) methods.
However, the previous works considered signal segments with preset properties and characteristics that describe perfect conditions (a subject laying on the mattress and breathing normally) and are thus distant from the real-life challenges when real data are processed. Indeed, these challenges are present due to the manifestations of spontaneous or illness-related activities of the subject during the recording of the signal. Examples of these instances include movement, coughing, holding breath, etc.
The main focus of this paper is to extend the systems developed in previous works to measure HR and RR by including a preprocessing stage with a set of classification and segmentation tools to decide which portions of the BCG signal, HR and RR are possible to measure. The full approach also allows the precise determination of the start and the endpoints of different kinds of human-body activities that can be extracted from the BCG signal.
We adopted a progressive approach, the flowchart of which is presented in Figure 1. The first step, after the BCG signal pre-processing, of content-based BCG signal exploration and understanding includes the temporal and spectral analysis of the BCG signal during different human body activities (for example, resting, coughing and movement). The analysis permits us to consider the use of two spectral features, namely the Spectral Flatness Measure (SFM) and Spectral Centroid (SC). These are defined according to a frame-by-frame analysis, which means that the BCG signal can be decomposed into pre-defined time duration slices. After features behavior analysis, a clustering task, which is the core of many segmentation methods, is carried out in order to determine the usefulness of selected features for classifying the BCG signal into many classes. The number of classes and the hyper-parameters related to the frame size are optimized thanks to the Gaussian Mixture Model (GMM). Next, in order to establish the link between the number of classes and the nature of human body activity, a supervised classification is proposed. Two types are defined: the first one is a two-class approach that aims to classify a selected frame into one permitting the measurement of HR (binary flag cardiac activity detection (CAD)/NoCAD as an output) and another one permitting the measurement of RR (respiratory activity detection (RAD)/NoRAD). We recall that this kind of binary output is useful to drive HR and RR detectors. In fact, during coughing, there is no need to try to measure the RR, for example, since information on respiratory activity is absent and not relevant. The second supervised classification is a refined one since it gives information about the type of human body activity (resting, coughing, expiration, movement, etc.). We deal with coarse and fine approaches. The coarse approach makes use of features calculated for each frame while the fine approach considers a sample-by-sample classification. To this end, frame-by-frame features are converted to time-series thanks to a novel idea presented in this paper. The segmentation consists of identifying the bounding of each human body activity.

2. Material and Methods

2.1. Data Collection

2.1.1. Acquisition Process

The project was carried out at the Centre Hospitalier de l’Université de Montréal (CHUM) in Montréal, Canada. The protocol was accepted by the institutional ethics committees of the Ecole des Technologies Supérieures (ETS), the TELUQ university, the University of Sherbrooke and CHUM.
The study was conducted on healthy subjects, comprising three males and three females, aged between 20 and 35 years, in good health and without known cardio-respiratory disorders. Once the subjects agreed to participate in the study, information concerning their demographic and anthropomorphic details (age, weight, height, etc. ) were collected.
The system used for collecting data included a small FOS mattress and a module to gather optical data from the mattress. This particular process and more details about the feasibility, conception and performance of the FOS micro-bend sensors are available in the literature (for example, [11] or [15]). The FOS mattress was fixed on the back of a regular office chair, as shown in Figure 2. The system used for collecting data included a small FOS mattress and a module to gather optical data. The raw data were sampled at 50 Hz by the module. During the whole experiment, the subject wore a Hexoskin and was asked to sit still on the chair for a duration of 5 min. During this time, a Raspberry Pi 3B recorded the data streamed on the serial port of the FOS’s module.

2.1.2. Experimental Protocol

To simulate common human body activities, subjects were asked to obey a certain experimental protocol composed of the following steps: (1) normal breathing for 30 s, (2) coughing three times, (3) normal breathing for 300 s, (4) holding breath for 30 s, (5) expiration, (6) normal breathing for 60 s, (7) coughing 10 times with a between-cough interval of 5 s, (8) normal breathing for 30 s, (9) coughing 10 times with a between-cough interval of 2 s, (10) normal breathing for 120 s, (11) standing up, (12) sitting down, (13) steps 10 and 11 repeated 4 times and (14) standing up.
Thus, aside from remaining in a still position, other human body activities that commonly occur and can alter vital sign measurements from BCG signals were introduced. We explored (i) coughing as a frequent activity executed by the human body, which induces a certain internal disturbance in the respiratory and cardiac rhythms; (ii) holding one’s breath as another activity that renders respiration activity impossible to consider while cardiac activity is possible to extract; and (iii) standing up and sitting down as an example of movement types.
Figure 3 illustrates the BCG signal during the experimental protocol activities. This is discussed in detail later.
The ground truth—i.e., the human body activities—was created manually by experts who adjusted the frontiers so that a precise ground truth was created. This label information was used later in frames and then in spectral features.

2.2. Ballistocardiogram Signal Description

2.2.1. Cardiac Information

The BCG signal is a measure of the acceleration of blood through the veins. In broader terms, it is a measure of the ballistic force generated by the heart. For a person placed on a bed with minimal friction and minimal movement artifacts, one cycle of ideal cardiac activity appearing in BCG signal is illustrated in the enlarged part of Figure 4. This contains eight different waves annotated as F, G, H, I, J, K, L and M. They are divided into three groups: pre-ejection ( F , G , H ), ejection ( I , J , K ) and diastolic ( L , M , N ). The F wave (rarely seen in the BCG representation) is closely related to pre-systolic events of the heart. The G wave corresponds to a trough preceding the systolic waves. The H wave is the first one in the systolic cycle; it corresponds to the maximum peak recorded (which is synchronous to the isovolumetric contraction) [17]. I, J and K waves are the most recognized waves in the BCG, occurring during the systole. The I wave simulates the acceleration of blood in the ascending aorta and pulmonary arteries; the J wave is the main positive wave occurring in the systole, simulating the acceleration of the blood when going through the second part of the aorta; and the K wave occurs before the end of the systole. L and M are the diastolic waves that represent headwards deflections following the K wave [18,19].

2.2.2. Respiratory Information

It is worth mentioning that the described waveform depends greatly on the used acquisition system. In fact, multiple sensors using different electrical, optical or mechanical devices are currently available [20]. For example, optical fiber-based systems extract the aforementioned waveform, but it is believed that the resulting signal presents a respiratory component which corresponds to the movement of the thoracic cage of the patient [11,21]. In fact, the BCG signal studied in this work was acquired using a microbend FOS which provides a new way of acquiring the mechanical activity of the human body. More precisely, we used a Juvo’s skin non-contact feature, benefitting from the patented invention of fibre-optic sleep and vital signs sensors [21,22]. The BCG signal formation was achieved through the intensity attenuation of the light passing through an optic fiber in response to a mechanical stimulus on the fiber, as shown in Figure 5.
The FOS sensor, typically placed at the level of the thoracic cage, captures the longitudinal BCG as well as the inspiration and expiration of the body [23]. In fact, during the inspiration phase, the inhalation of air by the patient causes the appearance of a downward trend in the BCG signal, which corresponds to the increase of the force executed on the sensor placed behind the thoracic cage of the body. During the expiration, the exhalation of the air in the body of the patient causes a decrease in the force on the sensor, leading to an upward trend. The coupling of these upward and downward trends constitutes the general respiration waveform of the BCG [21,24]. An example of the respiratory cycle present in the BCG signal acquired with the FOS system is illustrated in the upper part of Figure 4.

2.3. Spectral Features for the Identification of Various Activities

2.3.1. Spectral Features Motivation

The cardiac and respiratory activities were repeated periodically. The term x ( t ) denotes any one of the periodic signals generated in the BCG signal to describe the cardiac or the respiratory activity; according to its periodicity property, its Fourier transform is a sum of Dirac pulses equally spaced in the frequency axis at frequencies that are multiples of f 0 = 1 T 0 :
X ( f ) = k Z X k δ f k f 0 ,
where X k is the k t h Fourier transform coefficient, calculated during one period according to the following formula:
X k = 1 T 0 T 0 2 T 0 2 x ( t ) e j 2 π k T 0 t d t .
Thus, a periodic signal can be represented as a sum of sine waves, and thus the Fourier transformation of this particular signal is spiky.
We used the two features of the Spectral Flatness Measure (SFM) and Spectral Centroid (SC). They are computed as follows. First of all, the discrete-time BCG signal, denoted as x ( n ) , is decomposed into frames of short duration. These frames should be long enough to carry information about the activity but not overly long to avoid an overlap of two or more different activities. In the frequency domain, the short-term Fourier transform is calculated and its amplitude is extracted. Its module is denoted as | X ( m , k ) | , where m is the frame index and k is the discrete frequency.

2.3.2. Spectral Flatness Measure (SFM)

This is also known as Wiener entropy [25]. It is a signal processing measure used to describe the flatness of the spectrum of a certain signal, defined as the ratio of the geometric mean and arithmetic mean of the spectrum:
S F M ( m ) = ( N k = 1 N | X ( m , k ) | ) k = 1 N | X ( m , k ) | N ,
where m is the frame number, k is the frequency bin index and N is the number of frequency bins.
If the spectrum of the frame is flat, S F M ( m ) will be close to 1. In the particular case of white noise, the spectrum is constant, and thus the S F M ( m ) value is equal to one. The SFM is mostly used in audio-related fields, for instance. For example, we relate the detection of voiced and unvoiced speech frames [25], the analysis and recognition of whispered speech [26] and emotional speech synthesis and transformation [27]. However, few works have implemented such features in the framework of vital signals, which makes the subject of this work a novelty. In fact, we believe that the pertinence of this feature is supported by the intuition of their physical significance.

2.3.3. Spectral Centroid (SC)

The Spectral Centroid indicates the location of the center of mass of the spectrum. It is defined as follows:
S C ( m ) = k = 1 N | X ( m , k ) | . f ( k ) k = 1 N | | X ( m , k ) | ,
where f ( k ) is the frequency in Hertz related to the frequency bin k. Once again, because the Spectral Centroid is a good predictor of the “brightness” of a sound, it has been widely used in digital audio and music processing, such as musical genre classification [28]. To the best of our knowledge, its use for a BCG signal is also a novelty.

2.4. Unsupervised Classification for the Exploration of the Number of Activity Classes

We conducted an exploratory data analysis in order to determine the different classes of activities in the signal. First, a general description of the available features and their modeling process in an unsupervised manner was made. No prior knowledge about the available BCG signal was taken into consideration. This justified the existence of different classes and permitted us to define the hyper-parameters that were useful for the effective classification.

2.4.1. Features Histograms and GMM Modeling

The Gaussian Mixture Model (GMM) is used for unsupervised learning algorithms to model clusters of points. Each cluster is assigned a Gaussian function where the mean represents the center of the cluster and the variance describes the spread within it. We denote F ( m ) = [ S F M ( m ) , S C ( m ) ] T as the two-component feature vector of the frame m. The distribution can be modeled using the GMM as follows:
P ( F | μ , Σ ) = k = 1 K π k 2 π | Σ k | 1 / 2 e x p ( 1 2 ( F μ k ) T Σ 1 ( F μ k ) T ) ,
where K is the number of clusters (number of Gaussian components), π k is the weight of the k t h Gaussian (where the sum of all of them is equal to one), μ k is the mean, Σ k is the variance and | Σ k | is its determinant.
Classically, and regardless of the classification problem, many criteria are defined to formally determine the number of Gaussian components; for example, the Akaike Information Criterion ( A I C ) [29] and Bayesian Information Criterion ( B I C ) [30]. They are defined as follows:
A I C ( K ) = 2 K 2 log ( ML ) ,
where ML is the maximum likelihood of the model.
The Bayesian Information Criterion introduces a penalty term,
B I C ( K ) = M log ( M L ) 2 log ( ML ) ,
where M is the number of training examples. Similarly to A I C , low values of B I C indicate a good fit of the model.

2.4.2. Frame Size Optimisation

GMM modeling requires the optimization of the frame size. The used criterion is the silhouette, which provides information about the goodness of fit of a certain model [31]. It takes into consideration the between-cluster distance as well as the within-cluster distances, which makes it particularly interesting for an unsupervised learning context. The expression of the silhouette s ( i ) is the following:
s ( i ) = b ( i ) a ( i ) m a x ( b ( i ) , a ( i ) ) i f c a r d ( C i ) > 1 0 i f c a r d ( C i ) = 1 ,
where a ( i ) is the within-class distance, b ( i ) is the between-class distance, C i is the i t h cluster and c a r d ( . ) is the cardinal operator. a ( i ) and b ( i ) are computed as follows:
a ( i ) = 1 c a r d ( C i ) 1 j C i , i j d ( i , j ) ,
and
b ( i ) = min k i 1 c a r d ( C k ) j C k d ( i , j ) ,
where d ( i , j ) is the distance between the points i and j.
After calculating s ( i ) for each frame under different frame sizes, the mean value is determined for each frame size.

2.5. Supervised Binary Classification for Vital Signs Detection

Supervised classification was performed using the K-Nearest Neighbors (KNN) algorithm thanks to its simplicity and as we knew that other classification techniques have been tested and have shown poorer or similar performances. The KNN algorithm is a non-parametric method that is used for classification. The variant adopted was Fine K N N , which is the finest variation of K N N since it labels the new input with the same label as its nearest neighbor. The BCG signal frames were classified based on their pertinence to detect the heart rate and the respiratory rate. In this regard, all frames were labeled into one of the following activities: cardiac activity detection ( C A D ) or no cardiac activity detection ( N o C A D ) and respiratory activity detection R A D or no respiratory activity detection ( N o R A D ).

2.6. Multi-Class Classification and Segmentation

The BCG signal frames were then classified according to human body activities. Seven classes were used: normal activity in still position, coughing, post-coughing, holding breath, expiration, movement and others. Two different classification approaches were adopted. First, a classification was made for each frame, which is called frame-by-frame classification. The second approach was used for each sample and is called sample-by-sample classification.
Enframing is a classical method in signal processing, and it refers to splitting the BCG signal into temporal segments. For each frame, features are extracted, and decision-making about classification is related to the whole segment.
The approach used for sample-by-sample classification was as follows. Instead of using raw values of SFM and SC for each frame and risking the resulting limitations (coarse decision for long-duration frame, boundary epochs, etc.), we created time-series out of the SFM and SC values computed for each frame. The adopted approach is described below; it is also illustrated in Figure 6.
  •  First, the original BCG signal was decomposed into frames of length 1024 with an overlap of 960 samples (an increase of 64 from frame to frame). The Hamming window was used for this. An example of the original BCG signal, the moving Hamming window and the generated frames are illustrated in the top part of Figure 6.
  •  The features, namely S F M ( m ) and S C ( m ) , were extracted for each frame.
  •  Each scalar feature S F M ( m ) (resp. S C ( m ) ) was used to construct a “sub-time-series”. This latter was a constant vector, whose value was equal to S F M ( m ) (resp. S C ( m ) ) and whose length was equal to the frame size. The constant sub-series are shown in Figure 6 below those of previous steps (horizontal lines labelled F e a t u r e s ( m ) ).
  •  The whole set of sub-series was put together at the input of an “overlap and add module”. This step was equivalent to the inverse of windowing and frames decomposition; the aim of this process was to construct time-series from features that represented the spectral content of the BCG signal. Their length was equal to that of BCG signal, and they are denoted as s f m s i g n a l and s c s i g n a l , respectively, for the Spectral Flatness Measure and Spectral Centroid.
  •  The classification was then conducted on each sample of the new time-series. This process is called sample-by-sample-classification.

2.7. Evaluation Metrics

During this study, some common classification criteria were used. Due to a lack of space, only some of them are given here. In case of binary classification, we mainly retain the True Positive Rate ( T P R ) and the Positive Predicted Value ( P P V ). The T P R , called also sensitivity or recall, measures the proportion of detected positives from the actual positive (true positives and false negatives) while the P P V measures the proportion of detected positives from the real positive.
In the binary classification, the provided confusion matrix contains the T P R and the P P V values on the diagonal, while the complementary metrics define the False Discovery Rates ( F D R ) and the False Omission Rates ( F O R ). The generalized formulas of the evaluation metrics are included below:
T P R = t r u e p o s i t i v e s t r u e p o s i t i v e s + f a l s e n e g a t i v e s .
P P V = t r u e p o s i t i v e s t r u e p o s i t i v e s + f a l s e p o s i t i v e s .
In case of multiple classes (more than two), the confusion matrix in terms of T P R and P P V is used. Each term of the matrix is written as follows:
C o n f T P R ( i , j ) = M i j j = 1 L M i j ,
where M i j is the number of predictions of class i that belong to class j, L is the number of classes. The PPV confusion matrix is given by
C o n f P P V ( i , j ) = M i j i = 1 N M i j .

3. Results and Discussion

3.1. Illustration of Temporal and Spectral Properties

3.1.1. Overall Temporal Evolution

Figure 3 illustrates an example of the BCG signal acquired when following the experimental protocol. The different human body activities are plotted in different colors. The objective of this is to highlight the differences in the BCG signal according to the activity; for example, during normal respiration, the signal periodicity is observable, and so it is thought that we can easily measure the cardiac and respiratory vital signs using classic approaches. During the coughing exercise, it was difficult to detect the periodicity that originated from the respiration rhythm and the cardiac rhythm. In fact, the recording of the BCG shows that the signals attained the maximum value allowed by the acquisition equipment. This can also be noted during the movement phase. The peaks displayed in these different activities were not a result of the physiological activity of the body but of the physical activity. Moreover, in addition to the classes defined in the experimental protocol, a post-cough class has been added. Indeed, careful observation of this step showed particular variations in the BCG signal which did not appear with the other classes. Finally, another class called “others” includes other behaviors not foreseen in the experience.
To conclude, this particular representation, based on manual segmentation, highlights the importance of separately detecting the different activities in the BCG signal. The temporal and spectral contents of the BCG signal during some kinds of activities are detailed below.

3.1.2. Illustration during Quiet Activities

The resulting BCG signal acquired using a microbend FOS during normal activity in still position contained the cardiac waveform enriched with the respiration waveform. The top part of Figure 7a illustrates the temporal evolution of a BCG signal during this kind of activity. One can notice the periodicity of cardiac activity as well as that of respiratory activity. The respiratory periodicity is related to the large period appearing in the BCG signal. The cardiac periodicity is less than the respiration activity. In the frequency domain (bottom part of Figure 7a), one can see the presence of different spikes; they correspond to the respiratory frequency and its harmonics and to the cardiac frequency and its harmonics. In fact, the Fourier transform of a periodic signal is a sum of Dirac pulses equally spaced in the frequency axis at frequencies multiple of the fundamental frequency (see for example [32]). The respiratory frequency appears at the first spike; its first harmonic is the second one rising in the spectrum. However, higher-order harmonics do not appear because they are not relevant. On the other hand, the cardiac harmonics are relevant, and many of them rise significantly in the spectrum.
It is important to note that the values of cardiac and respiratory frequencies are validated by the available ground truth. In fact, in this study, we used the clinically validated Hexoskin intelligent textile. This was embedded with sensors giving accurate measurements [33,34]. It is also worth mentioning that more sophisticated methods can be used to extract these vital signs, such as the Wavelets Transform [35], discrete Fourier transformations [36], empirical mode decomposition [37] and many others (for example, [11]). These methods are, however, beyond the scope of this work.

3.1.3. Illustration during Special Activities

The former illustration describes a perfect environment as the subject remained in a still position. This situation is not always present because of the subject’s movement; thus, it is important to study the BCG signal behavior during various human body activities to determine if the BCG signal continues to manifest cardiac and respiratory vital signs.
The next subject activity studied was coughing. During this test, the subject was asked to cough while the BCG signal was acquired. This experience was studied as we believed that coughing activities reduced the reliability of extracting both vital signs. Furthermore, it was enlightening to capture the coughing frequency to establish information about the medical state of the patient. This application is promising especially in the context of SARS-CoV-2, where the coughing frequency can be an indicator of infectious disease [38].
Figure 7b represents the temporal evolution of the BCG signal (on the top) as well as the frequency representation (on the bottom) during four cough cycles. One can notice the existence of valleys with constant values due to the explosive part of the cough, which led to saturation in the lowest values. These valleys alternated with shapes that can be found in the BCG signal of normal activity. However, the periodicity induced by the respiratory activity was absent. On the other hand, a form of periodicity due to cardiac activity may exist, but it was challenging to accurately extract the corresponding waveform. From frequency representation, the spikes due to normal respiration disappeared. This fact calls into question the possibility of extracting vital signals during coughing.
Another activity performed by the subject was holding their breath. The subject was asked to hold their breath for a fixed amount of time. It was believed that the act of holding breath would nullify the effect of inspiration and expiration, and thus the effect of the respiratory activity. A temporal sequence and its frequency representation are displayed in Figure 7c. The signal portion contains only the cardiac information which allows the measurement of the heart rate.
According to these described situations, and to others that will be described later, it was established that measuring the heart and respiratory signs from BCG signals at all times is challenging and not possible in certain cases. It is thus crucial to carry a study that identifies relevant kinds of body activities in order to determine the possibility of vital signs measurement.

3.2. Illustration of Spectral Features

To illustrate, in a preliminary stage, the effectiveness of SFM and SC to discriminate between different kinds of activities that appear in the BCG signal, the following steps were followed. The signal was decomposed into frames of 20 s duration. This choice was justified as follows: the respiratory cycle is longer than the cardiac cycle, and techniques to measure the longest rate—i.e., the respiratory rate—require the existence of at least three respiratory cycles in the frame. Knowing that the lowest respiratory rate (highest period) is 10 respiration cycles per minute, three cycles were obtained during 3 60 10 = 18 s . We chose a duration of 20 s, giving 1000 samples for a sampling frequency of 50 Hz . It is important to note that this duration seemed to be long enough that boundary epochs (i.e., the time intervals for a transition between one body activity and another occurring) are handled. Unfortunately, these epochs are considered in the spectral features for frame-by-frame analysis, but thanks to the sample-by-sample analysis carried at the end of this study, this limitation was alleviated.
After framing, a Fast Fourier Transform FFT with a length of 1024 (the nearest power of two to 1000 samples) was applied and its module was extracted. Next, SFM and SC features were calculated according to Equations (3) and (4).

3.2.1. Scatter Plot

Figure 8 represents a scatter of these features for the still position, during coughing and during movement. One can notice distinct regions, confirming the ability of these features to discriminate between various activities. Each activity was described by a cluster of points. Moreover, the cluster of the (SFM, SC) pair showed lower values for quiet activity when compared to coughing. Moreover, motion seemed to increase the feature range values. In fact, Figure 8 illustrates the existence of certain decision boundaries to separate between pertinent frames and noisy frames (normal vs coughing and movement) but also a decision boundary to separate between the three classes since movement activities seemed to have a wider range of SFM values.

3.2.2. Features histograms and GMM modeling

Figure 9 shows the histograms of the selected features, namely SFM and SC. A total number of 1500 frames with a length of 1024 were obtained from the dataset. From Figure 9, one can notice the existence of two regions: the first one has low values of SFM and SC and the second one has greater values. The previously established remarks support the hypothesis that the first region corresponds to frames representing quiet physical activity. In these frames, cardiac and respiratory activities are relevant. In opposite, the second region, the higher values of SFM and SC correspond to human activity (movement) or physiological activities (coughing, holding breath, expiration, etc.).

3.2.3. Hyper-Parameter Optimization in Unsupervised Classification

Figure 10 shows the evolution of A I C and B I C as a function of the number of components in the GMM. This number infers the sub-optimal number of classes that can be adopted in the presented BCG signal. The most significant slope (in magnitude) matches the two classes benchmark. This supports the adoption of two classes if the two different classes should be interpreted as those of still activity and moving activity. If more refined classes are needed, Figure 10 shows that the optimal number of classes is around six or seven, as this corresponds to the minimum of the A I C and the B I C representations. This corresponds to the number of activities done by the subject and thus confirms the isolation of the post-cough activity from normal respiration, which corresponds with the theory of multi-labeling [39].

3.2.4. Frame Size Optimisation

In the previous experiments, the size of the frame was chosen equal to 1024 samples (corresponding to 20.48 s). This was chosen as an example to describe the determination of the number of classes. The next experiment aimed to optimize the frame size. Figure 11 shows the evolution of the mean value of s ( i ) versus the frame length. For high frame lengths, the model performs badly since frames can represent different states (different human body activities). The same was observed for very low frame lengths since short frames do not contain enough representative information to characterize the activities. The best trade-off was obtained for frame lengths of 512 or 1024. Note that the tested values are a power of two since the Fast Fourier Transform is applied on frames to compute the features.

3.3. Binary Supervised Classification Results

The first positive class, including quiet activity, post-coughing and holding breath, refers to the possibility of extracting the heart rate from the corresponding frames. The negative class includes movement and coughing. The classification results (in percentage) are presented in Table 1. The obtained results show the very good detection of the C A D class. However, a mis-classification of 33% for N o C A D is observed. This means that roughly a third of N o C A D is classified as relevant to measuring the heart rate. This result supports the idea that cardiac vitals can be extracted from coughing activities, which also explains the results obtained in Figure 8 where the values of SFM and SC occupy different regions for coughing and movement activities. It also indicates the presence of a cardiac component in the coughing activity compared to movement.
Comparing the T P R and P P V values in the confusion matrix, one can observe that the values representing the refinement of the detection of the model show good results for the negative class (85%); this confirms the previous hypothesis stating that certain activities classified as noisy represent cardiac features.
Next, dealing with binary classification based on respiratory rate detection, the first positive class including quiet activity refers to the possibility of extracting the respiratory rate from the corresponding frames. The negative class includes movement, holding breath and coughing. The classification results (in percentages) are presented in Table 2.
Similarly to the case of cardiac classification, good classification results were obtained in terms of the detection of the pertinent frames. However, a mis-classification of 28% in the class of N o R A D was recorded. This suggests the possibility of the extraction of the RR from the activities labeled as N o R A D . A possible explanation of this phenomenon is the possibility of the detection of the RR from other activities since this latter may contain the pertinent I , J , K sub-wave of BCG cycle (illustrated previously in Figure 4). In fact, the process of detecting the RR from the cardiac waveform is well established in the literature with the ECG signal. This latter contains the QRS complex commonly used to measure heart rate [40]. The obtained results support the possibility of the existence of similar extractions in the I , J , K complex.

3.4. Supervised Muti-Classes Results

The multi-class supervised classification corresponded to the seven classes describing major human body activities. They were as follows: normal activity in still position, coughing, post-coughing, holding breath, expiration, movement and others. The results in terms of T P R and P P V are represented in Table 3 and Table 4 for the case of frame-by-frame classification. The obtained classification accuracy was 75.5 % .
The results described in the T P R and P P V confusion matrices show mis-classifications in the class of coughing in terms of the T P R . This mis-classification means that the deployed model classifies all the coughing activities correctly but struggles to detect this class. This is shown by the 100% P P V result compared to only 38% for T P R . This particular interpretation inspired the use of s f m s i g n a l and s c s i g n a l since the process of creating such signals places more emphasis on the temporal evolution of the features and thus serves to detect the coughing class. Mis-classifications between the classes of holding breath and normal respiration present in the P P V matrix are intuitive since both classes present periodicity and represent one of the downsides of the adopted features. Mis-classifications in the class of expiration are mostly due to the availability of samples. We believe that the expiration class is underrepresented in this work, and more sophisticated methods for unbalanced classification such as one class classifiers [41] or model weights [42] can be adopted. However, this is beyond the scope of this research.
In this part, the sample-by-sample classification results are shown. Table 5 and Table 6 show the classification results of the adopted model. First, we can see an improvement in the detection of the coughing activity, as we achieved a T P R of 92% while utilizing the previously described signals. This supports the utility of generating s f m s i g n a l and s c s i g n a l . The mis-classification between the holding breath and normal respiration is still present, and this supports the previously established explanation and the lack of the adopted features in terms of distinguishing between these two classes. The overall classification accuracy using the sample-by-sample approach is 94.6%, which represents a huge improvement compared to the previously recorded 75.5%.

4. Conclusions

In this work, we proposed an approach to generate two binary flags indicating the useful frames permitting the measurement of cardiac and respiratory rates from a BCG signal. We have also presented a refined BCG signal segmentation technique according to the phenomena that occur in the process of recording the BCG signal. These described occurrences represent real-life challenges when it comes to the BCG signal applications since they present different properties than the perfect conditions assumed in the previous works. This fact makes this work a novelty and boosts its adaptability and use case for physical applications.
This work began with a temporal and spectral analysis of the resulting BCG signals. This was followed by unsupervised learning to understand the patterns available in the presented dataset. Furthermore, the different parameters of the segmentation were explored and rigorously determined. Next, two-label classification was carried out to create flags indicating the opportunity of measuring cardiac and respiratory activities. Finally, two different methods of multi-label classification were proposed. The sample-by-sample classification showed more promising results compared with frame-by-frame classification.

Author Contributions

S.B.J. and N.M. conceived the original idea of the paper; S.B.J. devised the idea and the signal processing code; M.C.B.N. developed the machine learning code, created the classification model, and performed data analysis. S.B.J. and M.C.B.N. labeled data. S.O. conceived and planned the experiments. S.O. deployed and managed the data collection. S.B.J. and M.C.B.N. wrote the manuscript. B.A. contributed to the final version of the manuscript; N.M. supervised the project. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Canada Research Chair on Biomedical Data Mining (950-231214).

Institutional Review Board Statement

This study was reviewed and approved by the TELUQ, ETS and CHUM institutional review board (IRB).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AICAkaike Information Criterion
BCGBallistocardiogram
BICBayesian Information Criterion
DFTDiscrete Fourier Transform
EMDEmpirical Mode Decomposition
FOSFiber Optic Sensor
GMMGaussian Mixture Model
IoMTInternet of Medical Things
KNNK-Nearest Neighbors
SFMSpectral Flatness Measure
SCSpectral Centroid
TPRTrue Positive Rate
PPVPositive Predicted Value
WTWavelet Transform

References

  1. Doyle, D. Clinical Early Warning Scores: New Clinical Tools in Evolution. Open Anesth. J. 2018, 12, 26–33. [Google Scholar] [CrossRef]
  2. Alam, N.; Vegting, I.; Houben, E.; Berkel, B.; Louella, V.; Kramer, M.; Nanayakkara, P. Exploring the performance of the National Early Warning Score (NEWS) in a European emergency department. Resuscitation 2015, 90. [Google Scholar] [CrossRef]
  3. Shao, D.; Tsow, F.; Liu, C.; Yang, Y.; Tao, N. Simultaneous monitoring of ballistocardiogram and photoplethysmogram using a camera. IEEE J. Biomed. Health Inform. 2017, 64, 1003–1010. [Google Scholar] [CrossRef] [Green Version]
  4. Mora, N.; Cocconcelli, F.; Matrella, G.; Ciampolini, P. Accurate Heartbeat Detection on Ballistocardiogram Accelerometric Traces. IEEE Trans. Instrum. Meas. 2008, 69, 9000–9009. [Google Scholar] [CrossRef]
  5. Inan, O.T.; Etemadi, M.; Wiard, R.M.; Giovangrandi, L.; Kovacs, G.T. Robust ballistocardiogram acquisition for home monitoring. Physiol. Meas. 2009, 30, 169–185. [Google Scholar] [CrossRef] [PubMed]
  6. Martin, S.L.O.; Carek, A.M.; Kim, C.S.; Ashouri, H.; Inan, O.T.; Hahn, J.O.; Mukkamala, R. Weighing scale-based pulse transit time is a superior marker of blood pressure than conventional pulse arrival time. Sci. Rep. 2016, 6, 39273. [Google Scholar] [CrossRef] [PubMed]
  7. David, H.; Winokur, E.S.; Sodini, C.G. An ear-worn continuous ballistocardiogram (BCG) sensor for cardiovascular monitoring. In Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, San Diego, CA, USA, 28 August–1 September 2012. [Google Scholar]
  8. Alametsä, J.; Viik, J.; Alakare, J.; Värri, A.; Palomäki, A. Ballistocardiography in sitting and horizontal positions. Physiol. Meas. 2008, 29, 1071–1087. [Google Scholar] [CrossRef]
  9. Puentes, J.; Rajeev, B.; Nilmini, W.; Raouf, N. Telemedicine trends and challenges: A technology management perspective. Int. J. Biomed. Eng. Technol. 2007. [Google Scholar] [CrossRef]
  10. GH, K.; Sheryl, N. Invasive electrodes in long-term monitoring. In Textbook of Epilepsy Surgery; Informa: London, UK, 2008. [Google Scholar] [CrossRef]
  11. Samuel, O.; Abdulrazak, B.; Ben Jebara, S.; Tournoux, F.; Mezghani, N. Smart mat for respiratory activity detection: Study in a clinical setting. In International Conference on Smart Homes and Health Telematics; Springer: Cham, Switzerland, 2019; pp. 61–72. [Google Scholar]
  12. Bruser, C.; Stadlthannert, K.; Brauers, A.; Leonhardt, S. Applying machine learning to detect individual heart beats in Ballistocardiograms. IEEE Eng. Med. Biol. Soc. 2010. [Google Scholar] [CrossRef]
  13. Shin, J.; Choi, B.; Lim, Y.H.; Jeong, D.; Park, K.S. Automatic Ballistocardiogram (BCG) beat detection using a template matching approach. IEEE Eng. Med. Biol. Soc. 2008, 1144–1146. [Google Scholar] [CrossRef]
  14. Sadek, I.; Biswas, J.; Abdulrazak, B.; Zhang, H.; Mokhtari, M. Continuous and unconstrained vital signs monitoring with Ballistocardiogram sensors in headrest position. In Proceedings of the IEEE-EMBS International Conferences on Biomedical and Health Informatics, Orlando, FL, USA, 16–19 February 2017. [Google Scholar] [CrossRef] [Green Version]
  15. Ramakrishnan, R.; Rajan, G.; Semenova, Y.; Farrell, G. Overview of fiber optic sensor technologies for strain/temperature sensing applications in composite materials. Sensors 2016, 16, 99. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Otis, S.; Abdulrazak, B.; Mezghani, N. Comparative study of heart rate extraction methods for a novel intelligent mattress. In Proceedings of the International Symposium on Signal, Image, Video and Communications (ISIVC), Rabat, Morocco, 27–29 November 2018. [Google Scholar] [CrossRef] [Green Version]
  17. Yebra, M.A.; Landreani, F.; Casellato, C.; Pavan, E.; Frigo, C.; Migeotte, P.F.; Caiani, E. Studying heart rate variability from ballistocardiography acquired by force platform: Comparison with conventional ECG. Comput. Cardiol. Conf. 2015, 929–932. [Google Scholar] [CrossRef] [Green Version]
  18. Luna-Lozano, P.; Alvarado-Serrano, C. Time and amplitude relationships of the ballistocardiogram in vertical and horizontal direction. In Proceedings of the 9th International Conference on Electrical Engineering, Computing Science and Automatic Control, Mexico City, Mexico, 26–28 September 2012; pp. 1–6. [Google Scholar] [CrossRef]
  19. Noh, Y.-H.; Kew, H.-P.; Jeong, D.-U. BCG monitoring system using unconstrained method with Daubechies Wavelet transform. In Proceedings of the 5th International Conference on Intelligent Manufacturing and Logistics System, Kitakyushu, Japan, 11–13 November 2009. [Google Scholar]
  20. Inan, O.; Migeotte, P.F.; Park, K.; Etemadi, M.; Tavakolian, K.; Casanella, R.; Zanetti, J.; Tank, J.; Funtova, I.; Prisk, K.; et al. Ballistocardiography and Seismocardiography: A review of recent advances. IEEE J. Biomed. Health Inform. 2014, 19, 1414–1427. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Lau, D.; Chen, Z.; Teo, J.T.; Ng, S.; Rumpel, H.; Lian, Y.; Yang, H.; Kei, P. Intensity-Modulated Microbend Fiber Optic Sensor for Respiratory Monitoring and Gating During MRI. IEEE Trans. Biomed. Eng. 2013. [Google Scholar] [CrossRef] [PubMed]
  22. Chen, Z.; Lau, D.; Teo, J.T.; Ng, S.H.; Yang, X.; Kei, P.L. Simultaneous measurement of breathing rate and heart rate using a microbend multimode fiber optic sensor. J. Biomed. Opt. 2014, 19, 057001. [Google Scholar] [CrossRef]
  23. Yoo, W.-J.; Jang, K.-W.; Seo, J.-K.; Heo, J.-Y.; Moon, J.-S.; Park, J.-Y.; Lee, B.-S. Development of respiration sensors using plastic optical fiber for respiratory monitoring inside MRI system. J. Opt. Soc. Korea 2010, 14, 235–239. [Google Scholar]
  24. Giallorenzi, T.; Bucaro, J.; Dandridge, A.; Sigel, G.H.; Cole, J.H.; Rashleigh, S.C.; Priest, R.G. Optical fiber sensor technology. IEEE J. Quantum Electron. 1982. [Google Scholar] [CrossRef] [Green Version]
  25. Gray, A.; Markel, J. A spectral-flatness measure for studying the autocorrelation method of linear prediction of speech analysis. IEEE Trans. Acoust. Speech Signal Process. 1974, 22, 207–217. [Google Scholar] [CrossRef]
  26. Ito, T.; Takeda, K.; Itakura, F. Analysis and recognition of whispered speech. Speech Commun. 2005, 45, 139–152. [Google Scholar] [CrossRef]
  27. Jiri, P.A.; Pribilova, A. Spectral flatness analysis for emotional speech synthesis and transformation. Lect. Notes Comput. Sci. 2009, 5641, 106–115. [Google Scholar] [CrossRef]
  28. Tzanetakis, G.; Cook, P. Musical Genre Classification of Audio Signals. IEEE Trans. Speech Audio Process. 2002, 10, 293–302. [Google Scholar] [CrossRef]
  29. Bozdogan, H. Model selection and Akaike’s Information Criterion (AIC): The general theory and its analytical extensions. Psychometrika 1987, 52, 345–370. [Google Scholar] [CrossRef]
  30. Neath, A.; Cavanaugh, J. The Bayesian information criterion: Background, derivation, and applications. Wiley Interdiscip. Rev. Comput. Stat. 2012, 4. [Google Scholar] [CrossRef]
  31. Rousseeuw, P.J. Silhouettes: A graphical aid to the interpretation and validation of cluster analysis. J. Comput. Appl. Math. 1987, 20, 53–65. [Google Scholar] [CrossRef] [Green Version]
  32. Lyons, R.; Fugal, D. The Essential Guide to Digital Signal Processing; Pearson Education, Inc.: Hoboken, NJ, USA, 2014. [Google Scholar]
  33. Al Sayed, C.; Vinches, L.; Hallé, S. Validation of a wearable biometric System ability to monitor heart rate in two different climate conditions under variable physical activities. IE-Health Telecommun. Syst. Netw. 2017, 6, 19. [Google Scholar] [CrossRef] [Green Version]
  34. Cherif, N.; Mezghani, N.; Gaudreault, N.; Ouakrim, Y.; Mouzoune, I.; Boulay, P. Physiological data validation of the Hexoskin smart textile. In Proceedings of the 11th International Conference on Biomedical Electronics and Devices, Funchal, Portugal, 19–21 January 2018; pp. 150–156. [Google Scholar] [CrossRef]
  35. Alvarado-Serrano, C.; Luna-Lozano, P.; Pallas-Areny, R. An algorithm for beat-to-beat heart rate detection from the BCG based on the continuous spline wavelet transform. Biomed. Signal Process. Control. 2016, 27, 96–102. [Google Scholar] [CrossRef]
  36. Moukadem, A.; Finnaoui, A.; Gassara, H.; Adolphe, D.; Schacher, L.; Dieterlen, A. Time-frequency domain for BCG analysis. In Proceedings of the 2018 International Conference on Computer and Applications (ICCA), Beirut, Lebanon, 25–26 August 2018; pp. 226–230. [Google Scholar] [CrossRef]
  37. Pinheiro, E.; Postolache, O.; Girao, P.S. Assessment of empirical mode decomposition implementation in cardiovascular signals. In Proceedings of the 17th IMEKO TC 4 Symposium, Kosice, Slovakia, 8–10 September 2010. [Google Scholar]
  38. Laguarta, J.; Hueto, F.; Subirana, B. COVID-19 Artificial Intelligence Diagnosis using only Cough Recordings. IEEE Open J. Eng. Med. Biol. 2020, 1, 275–281. [Google Scholar] [CrossRef]
  39. Tsoumakas, G.; Katakis, I. Multi-label classification: An overview. Int. J. Data Warehous. Min. 2009, 3, 1–13. [Google Scholar] [CrossRef] [Green Version]
  40. Sarkar, S.; Bhattacherjee, S.; Saurabh, P. Extraction of respiration signal from ECG for respiratory rate estimation. In Proceedings of the Michael Faraday IET International Summit 2015, Kolkata, India, 12–13 September 2015. [Google Scholar]
  41. Noumir, Z.; Honeine, P.; Richard, C. On simple one-class classification methods. In Proceedings of the IEEE International Symposium on Information Theory—Proceedings, Cambridge, MA, USA, 1–6 July 2012; pp. 2022–2026. [Google Scholar]
  42. Anand, A.; Pugalenthi, G.; Fogel, G.B.; Suganthan, P.N. An approach for classification of highly imbalanced data using weighting and undersampling. Amino Acids 2010, 39, 1385–1391. [Google Scholar] [CrossRef]
Figure 1. Flowchart of adopted methodology. BCG: Ballistocardiogram. CAD: cardiac activity detection; RAD: respiratory activity detection.
Figure 1. Flowchart of adopted methodology. BCG: Ballistocardiogram. CAD: cardiac activity detection; RAD: respiratory activity detection.
Sensors 21 01020 g001
Figure 2. Illustration of the experimental protocol activities [16].
Figure 2. Illustration of the experimental protocol activities [16].
Sensors 21 01020 g002
Figure 3. Illustration of the BCG signal during the experimental protocol activities. Each color corresponds to a human body activity.
Figure 3. Illustration of the BCG signal during the experimental protocol activities. Each color corresponds to a human body activity.
Sensors 21 01020 g003
Figure 4. Theoretical respiratory and cardiac activities present in the BCG waveform.
Figure 4. Theoretical respiratory and cardiac activities present in the BCG waveform.
Sensors 21 01020 g004
Figure 5. Microbend FOS principle. The light passing through the microbend FOS is modulated by the deformations in the optical fiber due to the displacement of the micro-benders [23].
Figure 5. Microbend FOS principle. The light passing through the microbend FOS is modulated by the deformations in the optical fiber due to the displacement of the micro-benders [23].
Sensors 21 01020 g005
Figure 6. Steps of time-series generation.
Figure 6. Steps of time-series generation.
Sensors 21 01020 g006
Figure 7. Illustration of temporal evolution and frequency representation of BCG signals. (a) normal respiration in still position, (b) during coughing activity, (c) while holding breath.
Figure 7. Illustration of temporal evolution and frequency representation of BCG signals. (a) normal respiration in still position, (b) during coughing activity, (c) while holding breath.
Sensors 21 01020 g007
Figure 8. Distribution of the Spectral Flatness Measure (SFM) and Spectral Centroid (SC) values in the dataset.
Figure 8. Distribution of the Spectral Flatness Measure (SFM) and Spectral Centroid (SC) values in the dataset.
Sensors 21 01020 g008
Figure 9. Histograms of SFM and SC values in the dataset.
Figure 9. Histograms of SFM and SC values in the dataset.
Sensors 21 01020 g009
Figure 10. Representation of the Akaike Information Criterion ( A I C ) and Bayes Information Criterion ( B I C ) as a function of the number of Gaussian components.
Figure 10. Representation of the Akaike Information Criterion ( A I C ) and Bayes Information Criterion ( B I C ) as a function of the number of Gaussian components.
Sensors 21 01020 g010
Figure 11. Variation of the mean of silhouette over the frames length.
Figure 11. Variation of the mean of silhouette over the frames length.
Sensors 21 01020 g011
Table 1. True Positive Rate ( T P R ) and Positive Predicted Value ( P P V ) confusion matrix for cardiac activity detection.
Table 1. True Positive Rate ( T P R ) and Positive Predicted Value ( P P V ) confusion matrix for cardiac activity detection.
True Class PPV TPR
CAD NoCAD CAD NoCAD
C A D 95%15%98%2%
N o C A D 5%85%33%67%
Table 2. T P R and P P V confusion matrix for respiratory detection.
Table 2. T P R and P P V confusion matrix for respiratory detection.
Predicted ClassPPV TPR
RAD NoRAD RAD NoRAD
R A D 90%17%95%5%
N o R A D 10%83%28%72%
Table 3. P P V confusion matrix by frame-by-frame classification.
Table 3. P P V confusion matrix by frame-by-frame classification.
ClassNormalCoughingPost-CoughingHolding BreathExpirationMovementOthers
Normal94% 28%30%33%21%45%
Coughing6%100%4%10%11%23%10%
Post-coughing 62%
Holding breath 3%57%
Expiration <1%1%56%
Movement 3%3% 55%
Others <1% 2%45%
Table 4. T P R confusion matrix by frame-by-frame classification.
Table 4. T P R confusion matrix by frame-by-frame classification.
ClassNormalCoughingPost-CoughingHolding BreathExpirationMovementOthers
Normal71% 22%2%1%3%1%
Coughing24%38%16%4%1%17%1%
Post-coughing 100%
Holding breath 35%65%
Expiration 16%4%80%
Movement 21%2% 77%
Others 32% 21%47%
Table 5. T P R confusion matrix for sample-by-sample classification.
Table 5. T P R confusion matrix for sample-by-sample classification.
ClassNormalCoughingPost-CoughingHolding BreathExpirationMovementOthers
Normal98% 1%<1%<1%
Coughing1%92%2%1% 4%
Post-coughing1%2%96% <1%1%
Holding breath31%<1% 68%<1%
Expiration 7%93%
Movement4% 1% 86%9%
Others 7% 93%
Table 6. P P V confusion matrix for sample-by-sample classification.
Table 6. P P V confusion matrix for sample-by-sample classification.
ClassNormalCoughingPost-CoughingHolding BreathExpirationMovementOthers
Normal96% 3%7%13%
Coughing<1%96%1%3% 7%
Post-coughing<1%4%95% 3%4%
Holding breath3%<1% 87%<1%
Expiration 3%84%
Movement1% <1% 90%23%
Others 1% 77%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ben Nasr, M.C.; Ben Jebara, S.; Otis, S.; Abdulrazak, B.; Mezghani, N. A Spectral-Based Approach for BCG Signal Content Classification. Sensors 2021, 21, 1020. https://doi.org/10.3390/s21031020

AMA Style

Ben Nasr MC, Ben Jebara S, Otis S, Abdulrazak B, Mezghani N. A Spectral-Based Approach for BCG Signal Content Classification. Sensors. 2021; 21(3):1020. https://doi.org/10.3390/s21031020

Chicago/Turabian Style

Ben Nasr, Mohamed Chiheb, Sofia Ben Jebara, Samuel Otis, Bessam Abdulrazak, and Neila Mezghani. 2021. "A Spectral-Based Approach for BCG Signal Content Classification" Sensors 21, no. 3: 1020. https://doi.org/10.3390/s21031020

APA Style

Ben Nasr, M. C., Ben Jebara, S., Otis, S., Abdulrazak, B., & Mezghani, N. (2021). A Spectral-Based Approach for BCG Signal Content Classification. Sensors, 21(3), 1020. https://doi.org/10.3390/s21031020

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