Next Article in Journal
An Efficient Method for Solving Router Placement Problem in Wireless Mesh Networks Using Multi-Verse Optimizer Algorithm
Next Article in Special Issue
Muscle Engagement Monitoring Using Self-Adhesive Elastic Nanocomposite Fabrics
Previous Article in Journal
A Novel Deep-Learning Method with Channel Attention Mechanism for Underwater Target Recognition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Motion Artifact Reduction in Wearable ECG Measurements Using Impedance Pneumography Signal

1
Academy of Artificial Intelligence, Beijing Institute of Petrochemical Technology, Beijing 102617, China
2
Beijing Academy of Safety Engineering and Technology, Beijing 102617, China
3
Research Institute for Flexible Materials, School of Textiles and Design, Heriot-Watt University, Galashiels TD1 3HF, UK
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(15), 5493; https://doi.org/10.3390/s22155493
Submission received: 30 June 2022 / Revised: 15 July 2022 / Accepted: 19 July 2022 / Published: 23 July 2022
(This article belongs to the Special Issue Advanced Wearable Sensors Technologies for Healthcare Monitoring)

Abstract

:
Noise is a common problem in wearable electrocardiogram (ECG) monitoring systems because the presence of noise can corrupt the ECG waveform causing inaccurate signal interpretation. By comparison with electromagnetic interference and its minimization, the reduction of motion artifact is more difficult and challenging because its time-frequency characteristics are unpredictable. Based on the characteristics of motion artifacts, this work uses adaptive filtering, a specially designed ECG device, and an Impedance Pneumography (IP) data acquisition system to combat motion artifacts. The newly designed ECG-IP acquisition system maximizes signal correlation by measuring both ECG and IP signals simultaneously using the same pair of electrodes. Signal comparison investigations between ECG and IP signals under five different body motions were carried out, and the Pearson Correlation Coefficient |r| was higher than 0.6 in all cases, indicating a good correlation. To optimize the performance of adaptive motion artifact reduction, the IP signal was filtered to a 5 Hz low-pass filter and then fed into a Recursive Least Squares (RLS) adaptive filter as a reference input signal. The performance of the proposed motion artifact reduction method was evaluated subjectively and objectively, and the results proved that the method could suppress the motion artifacts and achieve minimal distortion to the denoised ECG signal.

1. Introduction

In recent years, there has been an increasing interest in wearable electrocardiogram (ECG) monitoring systems. This interest comes from two aspects: societal needs and technology drive. Wearable ECG monitoring systems provide real-time vital sign monitoring, early disease detection, and telemedicine interventions, providing a good solution to the medical challenges posed not only by rapid population aging and healthcare cost inflation, but also by wellbeing requirements of an increasingly more health-conscious society. There is a huge demand for the latter, because as people are becoming more aware of their wellbeing and doing more exercise, these devices help them conduct their own, personal health monitoring independently. This has become possible with the miniaturization of electronics, improvements in the performance of low-power microprocessors, the ability of electronic textiles, and the development of artificial intelligence, rendering wearable healthcare monitoring systems to become increasingly portable, accurate, and powerful.
The ECG signal reflects the electrical activity of the heart, which is a weak physiological signal and hence very susceptible to various kinds of noise [1]. However, unlike the resting ECG measurements, wearable ECG measurements often face electromagnetic interference caused by the complex electromagnetic field environment in daily life, as well as motion artifact interference and EMG interference caused by various movements of the human body. The presence of noise in ECG signals may corrupt the ECG waveform causing inaccurate ECG interpretation, which in turn can lead to misdiagnosis of heart disease, inappropriate treatment decisions, and trigger false alarms. Therefore, to prevent inaccurate ECG waveforms and misdiagnosis of heart disease, reducing the noise in the ECG signal is a very important requirement before ECG signal analysis, especially for wearable ECG monitoring.
Motion artifacts have long been a problem in ECG measurement, the reduction of motion artifacts is difficult and challenging not only because the spectrum of motion artifacts is considerably overlapped with the important spectral components of the ECG signal, but also because the presence of the motion artifacts is inevitable and unpredictable in the wearable ECG measurement data. Motion artifacts are transient (but not step) baseline wander resulting from movement of the wearer, or the electrode [2]. The most widely accepted explanation of motion artifacts is from Tam and Webster [3], who found that the change in skin potential at the skin-electrolyte interface is the main source of motion artifacts. According to Edelberg’s skin model [3], the skin potential is determined by the potential across the sweat duct membrane, the total resistance in the sweat duct, the potential across the epidermis barrier membrane, and the total series resistance of the epidermis. Changes in any of the four variables will cause changes in skin potential.
The skin potential is usually negative, and the magnitude of the skin potential varies depending on the subject, skin site, and the condition of the skin. Floyd and Keele [4] found that increasing gel concentrations of NaCl or KCl increase the magnitude of the skin potential, and that the deformation of skin is the main cause of motion artifacts. Tam and Webster [3] found that the skin potential changes with pressure on the skin and when the skin stretches, the skin potential decreases. “If the thickness of the skin layer is changed by stretching or pressing down on the skin, the skin potential can vary by 5–10 mV compared to the 1 or 2 mV ECG signal” [3], which causes rapid changes in the ECG baseline. Consequently both the wearer’s movement and the respiration can cause deformation of the skin and hence introduce motion artifacts to the ECG signals.
To make wearable ECG monitoring systems practically applicable, motion artifact reduction is one of the most challenging problems that need to be overcome. The performance of different motion artifact reduction methods has been investigated in many studies [5,6,7,8,9,10,11,12,13,14,15,16,17,18], and adaptive filters have been shown to be useful in motion artifact reduction in some of those studies [19,20,21,22,23,24]. Tong et al. [25] used an anisotropic magnetoresistive (AMR) sensor and an accelerometer as the source of reference input, and their results indicate that adaptive filtering can reduce the amount of motion artifact present in the ECG, and the accelerometer-based sensor outperforms the AMR sensor. There are also some other reports that have used an accelerometer as a source of noise reference for adaptive motion artifact reduction; Raya and Sison [26] attached an accelerometer to the wearer’s back, Liu [24] placed an accelerometer next to the left leg electrode, and Lee et al. [27] placed an accelerometer between two electrodes. However, in a study of a wearable monitoring system based on textile electrodes, Zhang et al. [20] reported that motion artifacts were not effectively reduced by adaptive filtering with the accelerometer outputs as reference signals. The skin-electrode impedance signal was also reported to be a useful source as the reference signal for adaptive filtering [28,29,30]. Although the results from most literature prove that the adaptive filter can be used to reduce motion artifacts, there is no conclusion to determine which reference input signal can most effectively reduce motion artifacts caused by different motions [31], and some findings are contradictory, whilst others inconclusive.
In this article, we aimed to use an adaptive filtering technique to solve the problem of motion artifact in wearable ECGs. We have studied the motion artifact reduction performance of adaptive filters that use an impedance pneumography (IP) signal as a reference input signal. To ensure that the IP signal correlated with the motion artifact in ECG signals, we designed hardware that records ECG and IP simultaneously. The correlation between the IP and the motion artifact in the ECG signal was studied in order to test whether the IP signal is suitable as the reference input for adaptive filtering. Finally, we evaluated the performance of the proposed motion artifact reduction method with noisy ECG signals that were recorded while the subjects performed five typical upper body motions.

2. Methodology

2.1. Adaptive Filtering

The adaptive filtering method attracted our attention because it can adaptively track the signal under non-stationary conditions and can adjust its impulse response to filter out the noise in the input with little or no prior knowledge of the signal and noise characteristics [32]. The adaptive filter works generally for the adaptation of signal-changing environments, the spectral overlap between noise and signal, and unknown or time-varying noise [33]. It has the capability of adaptively tracking the signal under non-stationary conditions and can be used for different purposes, such as system identification, prediction, and noise cancellation.
Figure 1 demonstrates the working principle of the adaptive filter in noise cancellation. In adaptive noise cancellation systems, the objective is to produce a system output X’(n) that is the best fit in the least squares sense to the signal X(n) by feeding the system output back to the adaptive filter and adjusting the filter through an adaptive algorithm to minimize total system output power [32]. The adaptive filter contains a digital filter with adjustable coefficients and the adaptive algorithm to modify the values of coefficients for filtering each sample. The coefficients of the digital filter are continuously changed according to the chosen adaptive algorithm so as to minimize the mean squared value of the error signal e(n). The reference signal S(n) is fed into a digital filter to produce an output N’(n), which is as close as possible to the replica of the noise N(n). Subsequently, this filtered signal output y(n) is subtracted from the primary input d(n) to obtain the estimated desired signal X’(n) [5]. As can be seen from Figure 1, the performance of an adaptive filter in noise reduction is mainly determined by the adaptive algorithm and the reference input.
Least mean squares (LMS), Normalised Least Mean Squares (NLMS), and Recursive Least Squares (RLS) are the three most commonly used algorithms in adaptive filtering. In our preliminary trials, we found the three adaptive algorithms LMS, NLMS, and RLS are all capable of removing motion artifacts from the noisy ECG signal, but their noise reduction performance varies with the selection of the step size μ , or the forgetting factor λ . In our proposed adaptive motion artifact reduction method, we adopted the RLS algorithm because it provides a faster convergence speed. Furthermore, the forgetting factor of the RLS algorithm was reduced from 1 to 0.99 to improve the tracking ability of the adaptive filter. However, the reduction of forgetting factor will sacrifice the stability of the algorithm and thus cause a high-frequency noise in the denoised ECG, as demonstrated in Figure 2a, so we added a 5 Hz low-pass filter to the estimated noise N’(n) to reduce the unwanted high-frequency components, the results of the denoised ECG using an RLS adaptive filter and a 5 Hz low-pass filter are shown in Figure 2b. The cut-off frequency of the low-pass filter was set to 5 Hz because the main spectral range of the motion artifacts we wanted to suppress is 0–5 Hz [34].
As for the reference input signal, according to the principles of adaptive filtering, a suitable reference input signal should correlate in some way with the noise to be removed but uncorrelated with the desired signal [32]. However, in some practical cases, the available reference input signal may contain uncorrelated noise components in addition to the correlated noise components or contain low-level of desired signal components. When there is uncorrelated noise presented in the primary and reference input signals, the performance of the adaptive filter in noise cancellation will be degraded. In other words, a reference input signal that is highly correlated with the noise in the primary input is preferred. Moreover, when there are signal components presented in the reference input, some signal distortion will occur and hence the performance of the adaptive filter in noise cancellation will be degraded. The more signal components present in the reference input, the higher the signal distortion in the output. Therefore, in order to have a small signal distortion in the output, a reference input that contains small signal components is desired. In the application of adaptive motion artifact reduction in wearable ECGs, the reference input signal should come from a relatively independent source rather than noisy ECG signals to prevent the suppression and distortion of the ECG waveform.
In summary, a desired reference input signal for motion artifact reduction should have a good correlation with motion artifacts and should not be correlated with the ECG signal itself. Therefore, in our proposed method, we used specially designed hardware to simultaneously record both the ECG signal and the reference input signal in order to maximize the correlation between the primary input and the reference input. Furthermore, we applied a 5 Hz low-pass filter to eliminate uncorrelated noise and unwanted ECG components from the reference input signal. The proposed adaptive motion artifact reduction method is shown in Figure 3.

2.2. Impedance Pneumography

The Impedance Pneumography (IP) signal is a measure of the electrical impedance change of the subject’s thorax caused by respiration, as any increase in the volume of air respiration during each breathing cycle the electrical conductivity of the person’s thorax is reduced. In fact, if we consider the IP signal generated by the change in impedance as noise, it reflects the respiratory motion artifact as well as the motion artifact caused by skin deformation. Thus, from the origin of IP signals, there is a certain correlation between IP signals and the motion artifacts in ECG signals, which makes the IP signal a good option when looking for a reference input signal to adaptive reduce motion artifacts.
In order to suppress motion artifacts from ECG signals, we designed a wearable ECG-IP acquisition system consisting of a microcontroller MSP430F5438A (Texas Instruments, Dallas, TX, USA), an analogue front-end (AFE), a Bluetooth module RN41 (Microchip, Chandler, AZ, USA), and three textile electrodes. The AFE circuit board contains an analogue pre-amplifier ADS1292R (Texas Instruments, Dallas, TX, USA) and its peripheral circuitry to simultaneously record ECG and IP signals. The ADS1292R used in the AFE is of low power, and it is characterized by a high common-mode rejection ratio (−120 dB), 2-channels, and a 24-bit delta-sigma analog-to-digital converter (ADC) with integrated programmable gain amplifiers (PGAs). It is able to sample the ECG and IP signals, convert them from analog to digital and send data via a Serial Port Interface. With the integration of the respiration impedance measurement module, the ADS1292R measures the sum of the skin-electrode impedance and thoracic impedance. The respiration module in the ADS1292R contains a modulation block and a demodulation block. The modulation block produces a 2.42 V, 32 kHz square wave and injects it into the body as a carrier that is amplitude-modulated by a low-frequency signal generated as a result of the breathing action. The generated 32 kHz square wave is fed into the electrodes through the RESP_MODP and RESP_MODN pins, and the amplitude-modulated signal is received by the IN1P and IN1N pins of the ADS1292R. The measurement of the respiration impedance uses the same pair of electrodes for ECG measurement, so the measured IP signal has some correlation with the motion artifacts in the ECG signal. Figure 4 shows the designed analog front-end (AFE) circuit board. The configuration of the peripheral circuit of ADS1292R is referred to the recommendation of the technical datasheet by Texas Instruments [35].
Figure 5 shows examples of the recorded IP signal by the designed ECG -IP acquisition system.
As can be seen in Figure 5a, the measured IP signal contains not only the respiratory impedance signal but also the cardiac-related artifacts, as the transthoracic electrical impedance varies with both respiration and cardiac function [36]. The respiratory impedance signal is caused by the changes in lung volume during respiration because an increase in the amount of air in the lung reduces electrical conductivity. The cardiac-related artifacts [37] are mainly generated by changes in the blood volume during every cardiac cycle and have some correlation with ECG signals. In addition to the cardiac-related artifacts, the IP signal is also susceptible to body movement, as can be seen in Figure 5b. Changes in chest shape and/or in the skin-electrode impedance caused by body movement can produce large motion artifacts in the IP signals. In some cases, the motion artifacts are higher than the respiratory-related signals.
The aim of this work was to reduce the motion artifacts in the ECG signal by adaptive filtering, so the feasibility of using an IP signal as a reference input signal to adaptively reduce motion artifacts from the ECG signal was discussed. Since the ideal reference input for an adaptive filter should be correlated with the motion artifacts to be canceled but uncorrelated with the ECG signal, the motion artifacts in the IP signal should be preserved, but the cardiac-related signal should be eliminated. The frequency of cardiac-related artifacts in the IP signals is usually related to the QRS complex of ECG waves, which has a fundamental frequency of about 10 Hz [38]. As for the frequency of motion artifacts in IP signals, it is usually determined by the frequency of body motion, and the major frequency range of daily activities is 0.3–3.5 Hz [34]. Therefore, in our proposed adaptive motion artifact reduction method, we applied a 5 Hz low-pass FIR filter to the IP signal to suppress its cardiac-related artifacts before feeding it to the adaptive filter, as shown in Figure 3.

3. Experiments and Results

3.1. Experimental Signal Measurement

In the experiments, three textile electrodes were used to measure wearable ECG and IP signals, two of which were connected to the ELA and ERA pin of the AFE circuit board for 1-lead ECG recording and IP recording, and the remaining one textile electrode was connected as a ground electrode to the ERL pin of the AFE circuit board for common mode noise rejection. The fabrication of textile electrodes and their performance for ECG measurement has been discussed in another of our articles [39]. The three textile electrodes were sewn onto a knitted heavy stretch elastic band, as seen in Figure 6. During the measurement, the electrode chest band was placed on the horizontal line of the 6th intercostal space, where the electrodes connected to ELA and ERA were placed on the left and right midclavicular line, respectively, and the ground electrode was placed on the subject’s back.
In order to evaluate the performance of our proposed adaptive motion artifacts reduction method with the actually recorded noisy ECG signals, we simulated five common motions of daily life, and used the designed ECG-IP acquisition system to record the raw ECG and IP signals as disturbed by the different types of motion, as shown in Figure 7. A detailed description of these five types of motion is listed in Table 1. Ten female subjects and ten male healthy subjects between the ages of 21–26 participated in the experiment, each subject repeated each of the five types of motions ten times.
As shown in Table 1, Motion type I is to simulate the interference of external force on the electrodes. Motion types II–IV are to simulate three basic types of upper body movement. Motion type V is to simulate the respiration movement in a static state. In daily activities, various actions of the human body are mainly accomplished through the cooperation of muscles and joints of different body segments (such as upper limbs, upper body, and legs). During ECG measurement, because ECG electrodes are placed on the subject’s chest, only upper-body motions are more likely to produce motion artifacts. Therefore, in this study, three basic upper-body motions were selected.
Figure 8 shows the typical ECG and IP signals recorded in the experiment. It can be seen from Figure 8a that pressing the electrode causes severe baseline drift of the ECG signal, and the IP signal was also severely disturbed. The result indicates that the ECG and IP signals are very sensitive to external interference on the electrodes. As can be seen from Figure 8b–d, the upper body movement causes a certain amount of baseline drift in the ECG and IP signals. Figure 8e shows the typical signals measured when the subject was standing still, in which the respiratory motion artifact can be easily observed.

3.2. Signal Correlations

According to the principles of adaptive filtering, the performance of an adaptive filter in reducing motion artifacts is highly dependent on the characteristics of the reference signal. The higher the correlation between the reference signal and the motion artifacts, the better the performance in reducing motion artifacts. Therefore, the Pearson correlation coefficient was used as an indicator to determine whether the IP signal can be used as a reference input signal for the motion artifacts reduction.
The absolute value of the Pearson correlaion coefficient | r | is described as:
| r | = | COV   ( x , y ) σ x σ y | = | E ( ( x μ x ) ( y μ y ) ) σ x σ y |
where x is the 5 Hz low-pass filtered noisy ECG signal, y is the 5 Hz low-pass filtered IP signal; COV   ( x , y ) is the covariance and E is the expectation; μ x , σ x are the mean and standard deviation of x ; μ y , σ y are the mean and standard deviation of y . Note that in the equation, x should be a motion artifact signal, but in practice, the motion artifact signal is not independent but is included in the noisy ECG signal, so a 5 Hz low-pass filtered noisy ECG signal is used instead of the motion artifact signal to calculate the Pearson correlation coefficient. When | r | > 0.7 , it indicates a strong correlation; when 0.5 < | r | < 0.7 , it indicates a moderate correlation; when 0 < | r | < 0.5 , it indicates a weak correlation.
Figure 9 shows the calculated mean and standard deviation of the absolute correlation coefficient (|r|) under five body motions. The mean and standard deviation were calculated from the 200 measurements of each motion. As can be seen from Figure 9, the IP signal has a good correlation (|r| > 0.6) with noisy ECG signals in most body motions. The results indicate that the low-pass filtered IP signal has a good correlation with the corresponding motion artifacts, so it can be used as a reference input signal for adaptive motion artifact reduction.

3.3. Evaluation Parameters

Parameters such as the improvement of signal-to-noise ratio (SNR), correlation coefficient (c), mean-squared error (MSE), and R-square were often used to evaluate the noise reduction performance of different filters. However, the calculation of these parameters requires prior knowledge of the “clean” signal, which is the main reason that the simulated noisy ECGs were widely used in other relevant studies. But in our work, the actual recorded noisy ECGs were used instead of the simulated noisy ECGs, so there is no independent noise-free “clean” ECG signal that can be used to calculate parameters such as SNR.
Therefore, in this study, the concept “goodness-of-fit” [40] was used to evaluate the performance of the proposed adaptive motion artifact reduction method. According to the principles of adaptive filtering, the denoised ECG signal is the noisy ECG signal minus the estimated noise. When the estimated noise is close enough to the actual noise, the adaptive filter has good noise reduction performance. Thus, the goodness-of-fit between the estimated noise and the actual noise is an alternative measure to evaluate the noise reduction performance of adaptive filtering. To measure goodness-of-fit between the estimated noise and the actual noise, the estimated noise signal and the actual noise signal are essential. The estimated noise signal is the filtered signal output y(n), as described in Figure 3. The actual noise is obtained by low-pass filtering of the noisy ECG signal using an IIR filter with a cut-off frequency of 5 Hz. The 5 Hz low-pass filtered noisy ECG signal contains most motion artifacts but excludes most of the QRS components of the ECG signal. Two parameters were used to measure the “goodness of fit”, as denoted in the equation:
m _ MSE = 1 n i = 1 n ( d lowpass ( i ) y ( i ) ) 2
m _ R 2 = 1 i = 1 n ( d lowpass ( i ) y ( i ) ) 2 i = 1 n ( d lowpass ( i )   d ¯ lowpass ) 2
where d lowpass ( i ) is the 5 Hz low-pass filtered noisy ECG signal, y ( i ) is the estimated noise. The m_MSE value closer to 0 indicates a better fit, and the m_R2 value closer to 1 indicates a better fit.
The performance of the adaptive motion artifact reduction method was also compared with the 0.5 Hz IIR high-pass filter commonly used in practice to verify the efficiency of the proposed method. For this comparison, the standard deviation of the original ECG signal and the standard deviation of the denoised ECG signal were used to measure the reduction of motion artifacts. In statistics, the standard deviation is a measure of the amount of variation in a set of values. Because the motion artifacts are reflected as ECG baseline drift, a low standard deviation indicates less baseline drift (fewer motion artifacts), and a high standard deviation indicates more baseline drift (more motion artifacts). Therefore, the performance of the adaptive motion artifact reduction method and the IIR filter can be evaluated by comparing the STD values of the original noisy ECG and the denoised ECG.
The standard deviation (STD) is defined as:
STD = 1 N 1 i = 1 N | x i   x ¯ | 2
where x i is the value of the sample, x ¯ is the mean value of samples, N is the number of the samples.

3.4. Results and Discussion

Figure 10 shows the examples of the denoised ECG signal (red) compared to the original noisy ECG (blue). As can be seen from the graphs, the proposed adaptive motion artifact reduction method has good performance in reducing motion artifacts and thus preventing ECG waveform distortion. For objective evaluation of the performance of the proposed method, parameters m_MSE and m_R2 were calculated. The m_MSE value closer to 0 and the m_R2 value closer to 1 indicate that the performance of motion artifact reduction is better. Results are shown in Figure 11 and Figure 12.
As can be seen in Figure 11 and Figure 12, the m_MSE is smaller than 0.15, and the m_R2 is larger than 0.7. These results indicate that the proposed adaptive motion artifact reduction method has good performance in reducing motion artifacts. Figure 13 shows the results of the standard deviation (STD) of the original noisy ECG signal, IIR filtered ECG signal and the adaptive filtered ECG signal. It can be seen from the results that the STD of the noisy ECG signal is significantly reduced by both methods. The proposed method has a lower STD than the IIR filter, which proves that the proposed method can effectively reduce the motion artifact in ECG signals. Furthermore, the adaptive motion artifact reduction method has a distinct advantage over the high-pass IIR filter in reducing motion artifacts from ECG signals because it can adaptively track the baseline drift while preserving the ECG waveform.

4. Conclusions

Driven by social needs and technology, the demand for wearable healthcare monitoring devices is growing and is expected to accelerate in the future. Wearable ECG devices have attracted a lot of attention because it makes long-term ambulatory ECG monitoring possible, which helps to identify arrhythmias that are not detected by resting ECG or Holter ECG systems. However, due to the mobility of the human body and the various types of body movement in daily life, motion artifacts are a common and persistent problem for effective wearable ECG monitoring. The occurrence of motion artifacts can corrupt the ECG waveform causing inaccurate signal interpretation, so reducing their presence is necessary.
In this paper, an adaptive motion artifact reduction method based on an improved RLS adaptive filter is proposed to solve the motion artifact problem in wearable ECGs. In this method, the IP signal was used as the reference input signal of the improved RLS adaptive filter. In order to ensure correlation between the reference input signal IP and the motion artifacts in the ECG signal, we have specially designed the ECG-IP synchronous acquisition system discussed, to record both the ECG and the IP signal through the same pair of electrodes. To further improve the signal correlation, we applied a 5 Hz low-pass filter on the IP signal to eliminate the cardiac-related signal in the IP signal while preserving motion artifacts. According to the signal correlation study, we found that the processed IP signal has a good correlation (|r| > 0.6) with the motion artifacts, which proves that the IP signal is suitable as the reference input signal of the adaptive filter for motion artifact reduction.
Supported by the specially designed ECG-IP acquisition system, the performance of the proposed adaptive motion artifact reduction method was investigated by using real recorded noisy ECG signals. By analyzing the denoised ECG signals, it was found that most motion artifacts were removed and little distortion can be noticed from the denoised ECG waveform. By calculating the evaluation parameters, it was found that the mean squared error (m_MSE) is smaller than 0.15 and the R-square (m_R2) is larger than 0.7, which indicates that the estimated noise by the proposed method has a good fitness with the actual motion artifact. Finally, the results prove that the adaptive motion artifact method can suppress most of the motion artifacts and achieve minimal distortion to the denoised ECG signal. In addition, when processing pathological ECG signals, our proposed method should have better performance in preserving pathological characteristics of ECG waveform because it can suppress motion artifacts in the frequency range of 0 to 5 Hz while retaining all frequency components of the pathological ECG waveform. However, this could not be proved in experiments because ECG measurement on wearers who are suffering from heart problem may put wearers at risk, which is contrary to experimental ethics and in any case out of the scope of this research.
This work provides a solution to the noise problem in wearable ECG measurements. Since the presence of motion artifact in the ECG signal may corrupt the ECG waveform resulting in an incorrect diagnosis of heart disease in medical, or performance in sports and well-being, the motion artifact reduction system proposed in this paper can effectively improve the quality of the ECG signal and reduce the possibility of misdiagnosis. In the recent COVID-19 outbreaks, there was a huge demand for wearable health systems that measure heart rate and respiration [41], and these systems also faced motion artifact interference when measuring vital signs. Our method offers a possible solution to these systems, paving the way for truly wearable ECG systems of wider end uses.

Author Contributions

Conceptualization: X.A. and G.K.S.; Data curation: X.A. and Y.L.; Formal analysis: X.A.; Funding acquisition: Y.Z. and Q.L.; Investigation: X.A.; Methodology: X.A. and G.K.S.; Project administration: Q.L.; Resources: Q.L. and G.K.S.; Software: S.L.; Supervision: Q.L. and G.K.S.; Validation: X.A. and Y.Z.; Writing-original draft: X.A.; Writing-review&editing: G.K.S. All authors have read and agreed to the published version of the manuscript.

Funding

(a) This work was supported by the Climbing Program Foundation from Beijing Institute of Petrochemical Technology (Project No. BIPTAAI-2021-002); (b) This work was also supported by the fund of the Beijing Municipal Education Commission, China, under grant number 22019821001.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of Beijing Institute of Petrochemical Technology (protocol code: BIPTAAI002105016 and date of approval: 1 June 2021).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nayak, S.; Soni, M.K.; Bansal, D. Filtering techniques for ECG signal processing. Int. J. Res. Eng. Appl. Sci. 2012, 2, 671–679. [Google Scholar]
  2. Webster, J.G. Medical Instrumentation-Application and Design. J. Clin. Eng. 1978, 3, 306. [Google Scholar] [CrossRef]
  3. Tam, H.; Webster, J.G. Minimizing Electrode Motion Artifact by Skin Abrasion. IEEE Trans. Biomed. Eng. 1977, BME-24, 134–139. [Google Scholar] [CrossRef] [PubMed]
  4. Floyd, W.F.; Keele, C.A. Some observations on skin potentials in human subjects. Trans. Faraday Soc. 1937, 33, 1046–1049. [Google Scholar] [CrossRef]
  5. An, X.; Stylios, G.K. Comparison of Motion Artefact Reduction Methods and the Implementation of Adaptive Motion Artefact Reduction in Wearable Electrocardiogram Monitoring. Sensors 2020, 20, 1468. [Google Scholar] [CrossRef] [Green Version]
  6. Łęski, J.M.; Henzel, N. ECG baseline wander and powerline interference reduction using nonlinear filter bank. Signal Process. 2005, 85, 781–793. [Google Scholar] [CrossRef]
  7. Chouhan, V.S.; Mehta, S.S. Total removal of baseline drift from ECG signal. In Proceedings of the International Conference on Computing: Theory and Applications (ICCTA’07), Kolkata, India, 5–7 March 2007; pp. 512–515. [Google Scholar]
  8. Hao, W.; Chen, Y.; Xin, Y. ECG baseline wander correction by mean-median filter and discrete wavelet transform. In Proceedings of the 2011 Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Boston, MA, USA, 30 August 2011–3 September 2011; pp. 2712–2715. [Google Scholar]
  9. El-Dahshan, E.-S.A. Genetic algorithm and wavelet hybrid scheme for ECG signal denoising. Telecommun. Syst. 2011, 46, 209–215. [Google Scholar] [CrossRef]
  10. Karthikeyan, P.; Murugappan, M.; Yaacob, S. ECG signal denoising using wavelet thresholding techniques in human stress assessment. Int. J. Electr. Eng. Inform. 2012, 4, 306. [Google Scholar] [CrossRef]
  11. Singh, B.N.; Tiwari, A.K. Optimal selection of wavelet basis function applied to ECG signal denoising. Digit. Signal Process. 2006, 16, 275–287. [Google Scholar] [CrossRef]
  12. Alfaouri, M.; Daqrouq, K. ECG signal denoising by wavelet transform thresholding. Am. J. Appl. Sci. 2008, 5, 276–281. [Google Scholar] [CrossRef]
  13. Strasser, F.; Muma, M.; Zoubir, A.M. Motion artifact removal in ECG signals using multi-resolution thresholding. In Proceedings of the 20th European Signal Processing Conference (EUSIPCO), Bucharest, Romania, 27–31 August 2012; pp. 899–903. [Google Scholar]
  14. Zhang, D. Wavelet approach for ECG baseline wander correction and noise reduction. In Proceedings of the 27th Annual Conference on Engineering in Medicine and Biology, Shanghai, China, 17–18 January 2006; pp. 1212–1215. [Google Scholar]
  15. Blanco-Velasco, M.; Weng, B.; Barner, K.E. ECG signal denoising and baseline wander correction based on the empirical mode decomposition. Comput. Biol. Med. 2008, 38, 1–13. [Google Scholar] [CrossRef] [PubMed]
  16. Pan, N.; Mang, V.; Un, M.P.; Hang, P.S. Accurate removal of baseline wander in ECG using empirical mode decomposition. In Proceedings of the 2007 Joint Meeting of the 6th International Symposium on Noninvasive Functional Source Imaging of the Brain and Heart and the International Conference on Functional Biomedical Imaging, Hangzhou, China, 12–14 October 2007; pp. 177–180. [Google Scholar]
  17. Zhao, Z.-D.; Chen, Y.-Q. A new method for removal of baseline wander and power line interference in ECG signals. In Proceedings of the Machine Learning and Cybernetics, Dalian, China, 13–16 August 2006; pp. 4342–4347. [Google Scholar]
  18. Kabir, M.A.; Shahnaz, C. Denoising of ECG signals based on noise reduction algorithms in EMD and wavelet domains. Biomed. Signal Process. Control. 2012, 7, 481–489. [Google Scholar] [CrossRef]
  19. Thakor, N.V.; Zhu, Y.S. Applications of adaptive filtering to ECG analysis: Noise cancellation and arrhythmia detection. IEEE Trans. Biomed. Eng. 1991, 38, 785–794. [Google Scholar] [CrossRef]
  20. Zhang, Z.; Silva, I.; Wu, D.; Zheng, J.; Wu, H.; Wang, W. Adaptive motion artefact reduction in respiration and ECG signals for wearable healthcare monitoring systems. Med. Biol. Eng. Comput. 2014, 52, 1019–1030. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Iyer, V.; Ploysongsang, Y.; Ramamoorthy, P. Adaptive filtering in biological signal processing. Crit. Rev. Biomed. Eng. 1990, 17, 531–584. [Google Scholar]
  22. Correa, A.G.; Laciar, E.; Patino, H.D.; Valentinuzzi, M.E. Artifact removal from EEG signals using adaptive filters in cascade. J. Phys. Conf. Ser. 2007, 90, 012081. [Google Scholar] [CrossRef]
  23. Romero, I.; Geng, D.; Berset, T. Adaptive filtering in ECG denoising: A comparative study. In Proceedings of the 2012 Computing in Cardiology, Krakow, Poland, 9–12 September 2012; pp. 45–48. [Google Scholar]
  24. Liu, S.-H. Motion artifact reduction in electrocardiogram using adaptive filter. J. Med. Biol. Eng. 2011, 31, 67–72. [Google Scholar] [CrossRef]
  25. Tong, D.A.; Bartels, K.A.; Honeyager, K.S. Adaptive reduction of motion artifact in the electrocardiogram. In Proceedings of the Second Joint 24th Annual Conference and the Annual Fall Meeting of the Biomedical Engineering Society, Houston, TX, USA, 23–26 October 2002; pp. 1403–1404. [Google Scholar]
  26. Raya, M.A.D.; Sison, L.G. Adaptive noise cancelling of motion artifact in stress ECG signals using accelerometer. In Proceedings of the Second Joint 24th Annual Conference and the Annual Fall Meeting of the Biomedical Engineering Society EMBS/BMES Conference, Houston, TX, USA, 23–26 October 2002; Volume 2, pp. 1756–1757. [Google Scholar]
  27. Lee, W.-C.; Ou Yang, Y.-S.; Ke, T.-C.; Wei, C.-S.; Lee, H.-C. Adaptive reduction of motion artifact in a portable ECG system. In Proceedings of the Sensors, Waikoloa, HI, USA, 1–4 November 2010; pp. 704–707. [Google Scholar]
  28. Kim, S.; Yazicioglu, R.F.; Torfs, T.; Dilpreet, B.; Julien, P.; Van Hoof, C. A 2.4 µA continuous-time electrode-skin impedance measurement circuit for motion artifact monitoring in ECG acquisition systems. In Proceedings of the Symposium on VLSI Circuits, Honolulu, HI, USA, 16 January 2010; pp. 219–220. [Google Scholar]
  29. Hong, S.; Song, K.; Yan, L.; Yoo, H. A combined method to reduce motion artifact and power line interference for wearable healthcare systems. In Proceedings of the 2010 IEEE Asia Pacific Conference on Circuits and Systems, Kuala Lumpur, Malaysia, 6–9 December 2010; pp. 508–511. [Google Scholar]
  30. Romero, I.; Berset, T.; Buxi, D.; Brown, L.; Penders, J.; Kim, S.; Helleputte, N.V.; Kim, H.; Hoof, C.V.; Yazicioglu, R.F. Motion artifact reduction in ambulatory ECG monitoring: An integrated system approach. In Proceedings of the 2nd Conference on Wireless Health, San Diego, CA, USA, 10 October 2011; pp. 1–8. [Google Scholar]
  31. Ko, B.; Lee, T.; Choi, C.; Kim, Y.; Park, G.; Kang, K.; Bae, S.; Shin, K. Motion artifact reduction in electrocardiogram using adaptive filtering based on half cell potential monitoring. In Proceedings of the 2012 Annual International Conference of the IEEE Engineering in Medicine and Biology Society, San Diego, CA, USA, 28 August 2012–1 September 2012; pp. 1590–1593. [Google Scholar]
  32. Widrow, B.; Glover, J.R.; McCool, J.M.; Kaunitz, J.; Williams, C.; Hearn, R.H.; Zeidler, J.R.; Eugene Dong, J.; Goodlin, R. Adaptive noise cancelling: Principles and applications. Proc. IEEE 1975, 63, 1692–1716. [Google Scholar] [CrossRef]
  33. Tan, L.; Jiang, J. Digital Signal Processing: Fundamentals and Applications; Academic Press: Cambridge, MA, USA, 2018. [Google Scholar]
  34. Sun, M.; Hill, J.O. A method for measuring mechanical work and work efficiency during human activities. J. Biomech. 1993, 26, 229–241. [Google Scholar] [CrossRef]
  35. Texas Instruments. Low-Power, 2-Channel, 24-Bit Analog Front-End for Biopotential Measurements Datasheet; Texas Instruments: Dallas, TX, USA, 2012. [Google Scholar]
  36. Seppä, V.-P.; Hyttinen, J.; Viik, J. A method for suppressing cardiogenic oscillations in impedance pneumography. Physiol. Meas. 2011, 32, 337. [Google Scholar] [CrossRef]
  37. Brown, B.H.; Barber, D.C.; Morice, A.H.; Leathard, A. Cardiac and respiratory related electrical impedance changes in the human thorax. IEEE Trans. Biomed. Eng. 1994, 41, 729–734. [Google Scholar] [CrossRef] [PubMed]
  38. Bailey, J.J.; Berson, A.S.; Garson, A.; Horan, L.G.; Macfarlane, P.W.; Mortara, D.W.; Zywietz, C. Recommendations for standardization and specifications in automated electrocardiography: Bandwidth and digital signal processing. A report for health professionals by an ad hoc writing group of the Committee on Electrocardiography and Cardiac Electrophysiology of the Council on Clinical Cardiology. Am. Heart Association. Circ. 1990, 81, 730–739. [Google Scholar]
  39. An, X.; Stylios, G.K. A Hybrid Textile Electrode for Electrocardiogram (ECG) Measurement and Motion Tracking. Materials 2018, 11, 1887. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Schermelleh-Engel, K.; Moosbrugger, H.; Müller, H. Evaluating the fit of structural equation models: Tests of significance and descriptive goodness-of-fit measures. Methods Psychol. Res. Online 2003, 8, 23–74. [Google Scholar]
  41. Fyntanidou, B.; Zouka, M.; Apostolopoulou, A.; Bamidis, P.D.; Billis, A.; Mitsopoulos, K.; Angelidis, P.A.; Fourlis, A. IoT-based smart triage of Covid-19 suspicious cases in the Emergency Department. In Proceedings of the 2020 IEEE Globecom Workshops, Taipei, Taiwan, 7–11 December 2020; pp. 1–6. [Google Scholar]
Figure 1. Principle of the adaptive filter in noise cancellation.
Figure 1. Principle of the adaptive filter in noise cancellation.
Sensors 22 05493 g001
Figure 2. Samples of denoised ECG waveform: (a) denoised by RLS adaptive filter; (b) denoised by RLS adaptive filter with a 5 Hz low-pass filter.
Figure 2. Samples of denoised ECG waveform: (a) denoised by RLS adaptive filter; (b) denoised by RLS adaptive filter with a 5 Hz low-pass filter.
Sensors 22 05493 g002
Figure 3. The proposed adaptive motion artifact reduction method.
Figure 3. The proposed adaptive motion artifact reduction method.
Sensors 22 05493 g003
Figure 4. The analog front-end (AFE) circuit board: (a) schematic diagram; (b) circuit board.
Figure 4. The analog front-end (AFE) circuit board: (a) schematic diagram; (b) circuit board.
Sensors 22 05493 g004
Figure 5. The IP signal recorded by the ECG -IP acquisition system: (a) the IP signal contains cardiac-related artifacts; (b) the IP signal contains cardiac-related artifacts and motion artifacts.
Figure 5. The IP signal recorded by the ECG -IP acquisition system: (a) the IP signal contains cardiac-related artifacts; (b) the IP signal contains cardiac-related artifacts and motion artifacts.
Sensors 22 05493 g005
Figure 6. The chest band for wearable ECG measurement.
Figure 6. The chest band for wearable ECG measurement.
Sensors 22 05493 g006
Figure 7. Typical motion types to trigger the motion artifacts.
Figure 7. Typical motion types to trigger the motion artifacts.
Sensors 22 05493 g007
Figure 8. The typical ECG and IP signals recorded under different body motions: (a) motion type I, (b) motion type II, (c) motion type III, (d) motion type IV, (e) motion type V.
Figure 8. The typical ECG and IP signals recorded under different body motions: (a) motion type I, (b) motion type II, (c) motion type III, (d) motion type IV, (e) motion type V.
Sensors 22 05493 g008
Figure 9. The correlation between low-pass filtered ECG and IP signal.
Figure 9. The correlation between low-pass filtered ECG and IP signal.
Sensors 22 05493 g009
Figure 10. The denoised ECG signal: (a) motion type I, (b) motion type II, (c) motion type III, (d) motion type IV, (e) motion type V.
Figure 10. The denoised ECG signal: (a) motion type I, (b) motion type II, (c) motion type III, (d) motion type IV, (e) motion type V.
Sensors 22 05493 g010aSensors 22 05493 g010b
Figure 11. The m_MSE of the denoised ECG signal processed by the proposed method.
Figure 11. The m_MSE of the denoised ECG signal processed by the proposed method.
Sensors 22 05493 g011
Figure 12. The m_R2 of the denoised ECG signal processed by the proposed method.
Figure 12. The m_R2 of the denoised ECG signal processed by the proposed method.
Sensors 22 05493 g012
Figure 13. The performance of the IIR filter and the proposed method.
Figure 13. The performance of the IIR filter and the proposed method.
Sensors 22 05493 g013
Table 1. Different motion artifacts.
Table 1. Different motion artifacts.
Motion TypesDescription
IPressing of the electrodeOne electrode is pressed and released
immediately by hand.
IIUpper body anterior flexion and
extension
(upper body front bending)
Bend the spine forward to 45 degrees and then return upright.
Bending speed: about 0.2 cycles/sec
IIIUpper body lateral flexion and
extension
(upper body sideways bending)
Laterally bend the spine to 45 degrees then return to upright.
Bending speed: about 0.2 cycles/sec
IVUpper body rotationTwist the torso to 90 degrees and then
return.
Twisting speed: about 0.2 cycles/sec
VStandingStand still and breath normally
Respiration rate: about 0.3 Hz
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

An, X.; Liu, Y.; Zhao, Y.; Lu, S.; Stylios, G.K.; Liu, Q. Adaptive Motion Artifact Reduction in Wearable ECG Measurements Using Impedance Pneumography Signal. Sensors 2022, 22, 5493. https://doi.org/10.3390/s22155493

AMA Style

An X, Liu Y, Zhao Y, Lu S, Stylios GK, Liu Q. Adaptive Motion Artifact Reduction in Wearable ECG Measurements Using Impedance Pneumography Signal. Sensors. 2022; 22(15):5493. https://doi.org/10.3390/s22155493

Chicago/Turabian Style

An, Xiang, Yanzhong Liu, Yixin Zhao, Sichao Lu, George K. Stylios, and Qiang Liu. 2022. "Adaptive Motion Artifact Reduction in Wearable ECG Measurements Using Impedance Pneumography Signal" Sensors 22, no. 15: 5493. https://doi.org/10.3390/s22155493

APA Style

An, X., Liu, Y., Zhao, Y., Lu, S., Stylios, G. K., & Liu, Q. (2022). Adaptive Motion Artifact Reduction in Wearable ECG Measurements Using Impedance Pneumography Signal. Sensors, 22(15), 5493. https://doi.org/10.3390/s22155493

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