Adaptive Noise Reduction Algorithm to Improve R Peak Detection in ECG Measured by Capacitive ECG Sensors.

Electrocardiograms (ECGs) can be conveniently obtained using capacitive ECG sensors. However, motion noise in measured ECGs can degrade R peak detection. To reduce noise, properties of reference signal and ECG measured by the sensors are analyzed and a new method of active noise cancellation (ANC) is proposed in this study. In the proposed algorithm, the original ECG signal at QRS interval is regarded as impulsive noise because the adaptive filter updates its weight as if impulsive noise is added. As the proposed algorithm does not affect impulsive noise, the original signal is not reduced during ANC. Therefore, the proposed algorithm can conserve the power of the original signal within the QRS interval and reduce only the power of noise at other intervals. The proposed algorithm was verified through comparisons with recent research using data from both indoor and outdoor experiments. The proposed algorithm will benefit a noise reduction of noisy biomedical signal measured from sensors.

1. Introduction

The electrocardigram (ECG) is used to diagnose diseases, such as heart disease, by displaying bioinformation, such as R peak magnitude and interval [1,2]. Heart rate, which can be calculated using R peaks, often needs to be monitored for patients with diseases, as well as healthy individuals. To perform an ECG, Ag-AgCl electrodes must be attached to the patient’s body directly, which may be inconvenient in practice.

To overcome measurement limitations, capacitive ECG (cECG) sensors were developed to use capacitive coupling to measure voltage changes. Because ECGs in other intervals than QRS interval are relatively weak in power, the ECGs measured from the cECG sensors are mainly used to effectively detect R peaks, which can be used to derive a lot of medical information. Because of the convenience of cECG sensors, they are widely used in many situations where conventional ECG measurement is difficult but necessary [3,4,5]. For example, ECG measurement cannot be easily performed while driving, but evaluation of a driver may be necessary to prevent emergency situations that can be caused if a driver has a heart attack or falls asleep while driving. Therefore, cECG sensors could be mounted into the driver’s seat to obtain ECGs of drivers [6,7,8].

However, cECG sensors are sensitive to the noise caused by motion. As a result, accurate R peak detection is not possible. It is impossible to limit motion caused by breathing, vibrations, and various other actions of daily life; therefore, noise reduction via signal processing is necessary. Many methods for noise reduction have been suggested using bandpass filters [9], wavelet filters [10,11], independent component analysis [12], injection signals [13], spectro-temporal filtering [14], extended Kalman Filters [15], and active noise cancellation (ANC) [16,17,18]. Among them, adaptive filters, which are used in ANC, can effectively reduce noise that constantly changes by updating the filter weight adaptively as the noise changes. ANC has been studied to find a reference signals that can improve noise-reduction performance, but the commonly used reference signal that is obtained from acceleration sensors was not efficient in removing the noise from cECGs [17,18].

Recently, new reference signals obtained from additional cECG sensors were suggested to reflect the characteristics of the motion noise but not the bioinformation [6]. ANC using the new reference signal was verified to improve R peak detection by comparing the proposed algorithm to previous algorithms. Furthermore, reliability of the reference signal was shown by verifying the assumptions for ANC. However, the assumption for correlation with the reference noise was not clearly verified and it was sometimes wrong according to the subjects. Moreover, the adaptive filter of ANC was updated without considering that the power of the ECG within the QRS interval is much larger than the noise or the intermittence of the QRS interval. These problems caused a reduction in the power of the ECG and degraded R peak detection. Therefore, a new structure is needed for ANC to conserve the power of the ECG while reducing the power of the noise.

This study analyzes the properties of the ECG and the reference signal, and their effect on the update of the filter weight. As the filter weight within the QRS interval is updated as impulsive noise is added, we regard the measured ECG at the QRS interval to be impulsive noise. Accordingly, this paper proposes a new ANC structure using a robust variable step size affine projection sign algorithm (RVSS APA) that does not affect impulsive noise [19]. As a result, the power of the ECG signal within the QRS interval is conserved while the power of the noise is reduced to improve R peak detection.

2. Material

2.1. Measurement System

In recent research [6], four cECG sensors were mounted in a row on the back of a chair and driver’s seat (Figure 1). Conductive fabric was placed on the bottom of the seat and connected to a driven right leg circuit to reduce the common noise component of the measured ECG. cECGs and were obtained from the two sensors outside in a row and additional cECGs and were obtained from two other sensors, which were close to the center. This study used signals that were defined and used in the existing research [6]. The difference between and is filtered through a band pass filter (BPF) with a 0.05–35 Hz pass-band; this was defined as a measured ECG ():

Figure 1

The measurement system.

Two reference signals, and , were obtained from the two sensors of each part as and were selected because they contain motion information, but do not have bioinformation [6]. Then, we used an NI 9205 (National Instruments, Austin, TX, USA) with 16 bits to convert the obtained analog signals to digital signals at 200 Hz.

2.2. Data Acquirement

The data were acquired from both indoor and outdoor experiments. Nineteen subjects participated in the indoor experiment.The participants were between 20 and 30 years old; one was female and the others were male. Ten subjects participated in the outdoor experiment. The participants were between 20 and 30 years old; all were male. All subjects wore cotton T-shirts during the experiments. During both experiments, cECGs were obtained from the cECG sensors. In addition, was simultaneously obtained from conventional contact electrodes to obtain true R peaks.

During the indoor experiment, the subjects sat on the chair in which the measurement system had been installed. Then, they conducted the indoor experiment for 30 min according to following guidelines. For the first 5 min, the subjects remained stable without any movement to obtain an cECG without the noise. Leaning their backs against the cECG sensors, the subjects moved their bodies randomly to obtain an cECG with the noise for 25 min. The size and shape of the movement varied, and short breaks were given between movements. As a result, there was noise with a variety of power and frequency in the measured cECG, as can be assumed to be true for many situations.

During the outdoor experiment, the subjects sat in the driver’s seat, in which the measurement system had been installed (similar to the chair in the indoor experiment). The subjects were asked to drive along a course in the Pohang, Kyungbuk province of Korea. It took an average of 55 min for the participants to drive 31.8 km. The driving course consisted of 2.2 km of on-campus roads, 9.3 km of city roads, 9.2 km of national highway, and 11.1 km of highway. There were road signs and variations along the course, such as signal-induced stops, speed bumps, and unpaved roads. As a result, various types of noise were mixed into the measured cECG. In this experiment, the subjects were required to drive normally and not be concerned with leaning against the cECG sensors.

R peaks in the ECGs obtained during the experiments were detected using an R peak detection algorithm with the Pan and Tompkins algorithm [20] owing to ease of implementation and prior researches [6,16,21]. was used to obtain true R peak positions and to evaluate the accuracy of obtained R peaks from the .

3. Methods

3.1. Active Noise Cancellation

3.1.1. ANC Using Affine Projection Algorithm

ANC is used to remove noise from the measured signal [22,23,24] and its block diagram is depicted in (Figure 2). is a sequence of and equals , where i denotes the iteration. is the original ECG and is the noise from the motion of the subject and the vehicle. The reference signal is filtered through an adaptive filter to produce the output . Then, the error signal is obtained and used as the output of ANC.

Figure 2

The block diagram of ANC.

During ANC, noise can be cancelled when is correlated with , and is not strongly correlated with and [6,22]. The expectation value of squared is obtained as where denotes the expectation. The weight of the adaptive filter is updated to minimize so that . As a result, approaches zero and approximates . The reduction of noise depends on the extent to which the correlation conditions are satisfied. Therefore, the selection of the reference signal is the most important step and has been researched to satisfy these conditions [6,17,18]. Reference noise , determined using the differences between the cECGs from adjacent sensors, has been suggested to improve R peak detection by reducing the noise of [6]: where L denotes half the length of the filter weight and T denotes its transpose. In addition, an affine projection algorithm (APA) is used to update the weight of the adaptive filter. As itself is highly correlated and the correlation reduces the convergence rate of the filter weight [25], the APA can improve the convergence rate. The filter weight of the APA is updated every iteration as follows: where P is the projection order and is the step size of the APA. P determines the number of input vectors of and determines the convergence rate of the filter weight. The convergence rate is important in ANC because it determines the noise reduction performance when noise with different statistical characteristics is added. is the error vector and its first-tap error is the output sequence of ANC.

3.1.2. Properties of Reference Signal and ECG

Although ANC reduces noise, it cannot improve R peak detection because it reduces the power of QRS interval (Figure 3) as well. To investigate this problem, has to be examined. In and , is expected to be removed and only the signal correlated with the is expected to remain only motion information by subtracting each other’s cECG obtained from the adjacent cECG sensors. However, and are still correlated with because is not completely removed (Figure 4c). As a result, the correlation condition is not satisfied and the power of is reduced after filtering. In addition, the QRS interval of is intermittent and the power of at the QRS interval is large enough to ignore . Therefore, within the interval, Equation (4) approximates . As a result, the filter weight is updated to remove from the output signal of the ANC even though the filter weight had been updated to reduce before the QRS interval. When (Figure 4a) is filtered through ANC, the power of the filtered ECG (Figure 4b) decreases at the QRS interval and more incorrect R peaks are detected, where the green o denotes the true R peak position and the red * denotes detected R peaks. Within the QRS interval, the first-tap filter weight dynamically changes (Figure 4d) and the updated quantity of the total filter weight increases (Figure 5a) compared to the S-Q interval (Figure 3). The filter weight has been continuously updated to remove the noise that exists in the S-Q interval, but within the QRS interval is large within the short time frame, causing the filter weight to be updated significantly to remove .

Figure 3

Definition of QRS interval, S-Q interval, , and in the ECG.

Figure 4

ECG filtered by the APA: (a) , (b) ECG filter by the ANC with reduced power at the QRS interval, (c) , , and , and (d) the trajectory of the first-tap filter weight.

Figure 5

Update quantity of the filter weight: (a) update quantity by APA and (b) update quantity by the proposed algorithm using RVSS APA.

To conserve the power of at the QRS interval, a dramatic update of the filter weight should be prevented by slowing the convergence rate. For a slow convergence rate, step size should be small; however, step size adjustment is hard because there is a serious trade-off problem. When the step size is small for slow convergence rate, the power of within QRS interval is conserved, but the power of is not reduced. In contrast, when the step size is large, the power of is reduced. In addition, the step size adjustment depends on the characteristics of each subject’s ECG.

To improve R peak detection with ANC, the filter weight at the QRS interval should be updated only for baseline noise reduction. Further, at the S-Q interval, which is defined as the entire interval except for the QRS interval (Figure 3), the filter weight should converge fast to reduce the noise. Based on the properties of , we regard at the QRS interval as impulsive noise that requires a dramatic update of the filter weight when it is added. Therefore, robust algorithm against impulsive noise should be applied at the QRS interval, whereas the APA works for noise reduction at the other interval.

3.2. Robust Variable Step Size APA

Robust variable step size (RVSS) APA was introduced [19]. The RVSS APA updates its filter weight while switching the APA with the affine projection sign algorithm (APSA) [25]. The cost function of the APSA is -norm of the error vector and the constraint is the update quantity of the filter weight , where is a parameter for limited update quantity [26]. The update equation is obtained by solving the minimization problem: where denotes the sign function. APSA is robust against dynamic changes in filter weight, and its convergence rate is slow, so it should be used when update quantity by APA is large. Furthermore, the APA should be used to update the filter weight when the robustness of APSA is not necessary. In RVSS APA, is updated by the APA when the update quantity of the APA does not exceed threshold . Otherwise, is updated by the APSA with update quantity:

In the equation to update the APSA of the RVSS APA, is replaced with . The is recursively obtained as follows: where is a memory factor within [0,1] and generally close to one.

3.3. Proposed Algorithm with RVSS APA

At the QRS interval, the APSA is used to conserve power and reduce only the baseline noise because is considered impulsive noise. At the S-Q interval, APA is used to reduce the noise where is dominant. To satisfy with both goals, RVSS APA should be used with appropriate parameters, such as threshold and step sizes for the ECG-denoising system: where is the APA step size, is the APSA step size, and is the threshold that determines when APSA is used.

3.3.1. Step Size and

The parameter decides the convergence rate of APA. The noise of the ECG is reduced by the APA especially within the S-Q interval. When becomes large and close to one, the noise is better reduced because the fast convergence rate can quickly reflect the constantly varying characteristic of the noise. However, when is too large, it can remove Q and S waves at the boundary of QRS interval. As a result, R peak detection can be degraded because Q and S waves are important to detect precise R peaks. Therefore, does not have to be too large.

In the APSA, equals [19,25,26]. However, is not small and constantly varies because the filter weight has been constantly updated due to the varying system. As a result, at the QRS interval can be removed with , which is not desirable. Therefore, has to be smaller than , but not zero. If is zero, the filter cannot adjust for the constantly varying system change. Therefore, is a function of to reflect the system change: where the is a scaling factor less than one.

3.3.2. Threshold

During every iteration, the update quantity of the APA is compared to , and APSA is used instead of APA when the update quantity exceeds . The update quantity of the filter weight is constantly changing due to experimental environment factors, such as ambient noise, the subjects, and experimental equipment. If is constant, only the APSA is used so that the noise is not removed, or only the APA is used so that at the QRS interval is removed. The has to be at the boundary of QRS and S-Q intervals. Therefore, has to be a function of the update quantity of the APA as is compared with only the update quantity of the APA as follows:

From this update, can reflect the constantly varying update quantity of the APA. In addition, is a scaling parameter to find only the boundary of the QRS and S-Q intervals. The update quantities of the APA at every iteration are shown in (Figure 5a) when is filtered through the proposed algorithm using RVSS APA. should be within an appropriate range to locate the boundary and detect only the QRS interval. The update quantity is adjusted by the proposed algorithm (Figure 5b).

3.4. Postprocessing

In the research to developethe reference signal, a simple and reasonable postprocessing method is used [6] to prevent incorrect updates owing to exceptional situations that occur in practice. The filtered ECG was verified by two rules using two intervals from three consecutively detected R peaks and the power of ECG. When the filtered ECG is abnormal, its raw ECG is used instead. Owing to the simplicity and reliability of the postprocessing and the fair comparisons in following experiments, postprocessing is used in the proposed algorithm.

4. Experimental Results and Discussion

4.1. Performance Index

To verify the proposed algorithm with the new structure of the ANC, sensitivity () and positive predictivity () were used as performance indices for detected R peaks. and is obtained by

A search window is defined as a 80-ms period centered on each true R peak [6]. True positive (TP) denotes the number of accurately detected R peaks within the search windows, false negative (FN) denotes the number of undetected R peaks within the search windows, and false positive (FP) denotes the number of detected R peaks outside the search windows. decreases when the peaks to be found are not detected and decreases when the wrong peaks are recognized as R peaks. In addition, their average is defined as .

Also, is used as a signal quality index. Because and are difficult to distinguish and obtain in practice, we assumed the original ECG power is the same as the power of the 100 ms centered around a true R peak as we considered the dominance of their power based on the quality index of ECG [27]. In the ECG measured by the cECG sensors, the complete PQRST wave is difficult to obtain, but complete QRS interval can be obtained. Therefore, the S-Q interval is considered the interval including only the noise. In this study, is defined as where is the power of the 100-ms ECGs centered on the detected R peaks and is the power of the ECG outside those 100-ms ECGs (Figure 3). decreases when the power at the QRS interval is not conserved by filtering.

4.2. Experimental Environment

To consider near-real-time noise reduction of the ECG, a 6-s processing window was set to process the length of the ECG, which the previously used 2-s data overlaps. The overlapped data enable R peaks at the boundary to be detected stably.

, , and , as obtained from , , and in the indoor and outdoor experiments, are shown in Table 1, where is the ECG filtered by the ANC using APA [6] and is the ECG filtered by the proposed algorithm using RVSS APA. The same postprocessing was also applied to ANC. For both experiments, the step size of ANC using APA was experimentally set to 0.05, which can lead to the best performance for the comparison. Furthermore, , , and were set as 0.4, 0.01, and 3.5. and were obtained using the entire length of each data set without exception. However, was obtained using only the initial 5-min length to ignore the baseline noise and noise during driving and motion because power is sensitive to large magnitudes of noise even if they do not last long. An increase in compared to that of is defined as and shown in each table to identify performance differences clearly.

Table 1

, , and in the indoor and outdoor experiments.

	Data	ECG_m				ECG_APA					ECG_RVSS
		Se	P^+	Avg		Se	P^+	Avg	d_Avg		Se	P^+	Avg	d_Avg
Indoor	1	80.57	70.00	75.28		87.42	82.50	84.96	9.68		86.26	84.97	85.61	10.33
	2	92.42	87.83	90.13		94.84	92.69	93.77	3.64		95.25	95.35	95.30	5.18
	3	95.83	95.98	95.91		96.38	97.11	96.74	0.84		96.34	97.38	96.86	0.95
	4	91.74	94.02	92.88		93.77	93.44	93.61	0.72		92.58	94.46	93.52	0.64
	5	91.76	82.69	87.23		94.64	89.91	92.28	5.05		94.79	94.79	94.79	7.56
	6	85.10	82.20	83.65		90.81	89.06	89.93	6.29		88.32	89.65	88.99	5.34
	7	91.48	86.93	89.20		95.16	90.53	92.85	3.64		94.20	94.20	94.20	5.00
	8	40.55	38.40	39.47		52.84	48.35	50.59	11.12		49.49	53.38	51.44	11.97
	9	92.11	69.62	80.87		93.55	78.61	86.08	5.21		93.73	88.38	91.05	10.19
	10	58.21	54.84	56.52		62.38	56.54	59.46	2.94		62.48	61.27	61.87	5.35
	11	52.13	39.96	46.05		54.00	40.42	47.21	1.16		63.93	61.86	62.89	16.84
	12	65.76	63.34	64.55		72.21	66.94	69.58	5.03		68.19	69.12	68.65	4.10
	13	63.77	49.24	56.50		69.19	52.69	60.94	4.44		64.82	58.12	61.47	4.97
	14	90.11	72.93	81.52		95.67	85.91	90.79	9.27		96.33	94.63	95.48	13.96
	15	83.87	70.05	76.96		84.48	71.03	77.75	0.79		84.48	82.59	83.54	6.58
	16	74.74	63.59	69.16		83.21	76.32	79.77	10.60		83.13	85.10	84.12	14.95
	17	73.37	84.40	78.88		82.61	87.61	85.11	6.23		82.18	91.18	86.68	7.80
	18	95.44	90.97	93.20		95.87	93.99	94.93	1.73		95.83	96.71	96.27	3.07
	19	81.30	62.92	72.11		92.70	86.77	89.73	17.62		93.13	94.07	93.60	21.49
	Mean	78.96	71.57	75.27		83.77	77.92	80.85	5.58		83.45	83.54	83.49	8.22
Outdoor	1	82.31	72.20	77.26		84.74	76.99	80.86	3.61		86.11	84.18	85.15	7.89
	2	90.17	74.99	82.58		90.09	78.20	84.14	1.56		90.86	84.56	87.71	5.13
	3	86.69	90.38	88.53		87.14	90.30	88.72	0.18		86.48	90.55	88.52	-0.02
	4	80.04	64.79	72.42		86.26	74.47	80.37	7.95		87.37	84.92	86.15	13.73
	5	85.24	70.70	77.97		88.35	77.79	83.07	5.10		88.30	87.43	87.87	9.89
	6	93.52	89.91	91.72		96.47	92.90	94.68	2.96		96.34	95.68	96.01	4.30
	7	81.13	79.93	80.53		82.94	83.15	83.05	2.52		82.35	84.74	83.55	3.02
	8	93.35	84.34	88.84		94.44	89.34	91.89	3.05		95.18	94.09	94.64	5.79
	9	88.32	79.14	83.73		91.99	87.38	89.69	5.96		92.49	92.76	92.63	8.90
	10	80.03	64.97	72.50		80.52	67.62	74.07	1.57		79.42	73.98	76.70	4.20
	Mean	86.08	77.14	81.61		88.29	81.81	85.05	3.45		88.49	87.29	87.89	6.28

The length of the data from the outdoor experiment is so long that the increase seen is small. Therefore, the data were divided into 5 sections based on the roads that were driven. and are shown in Table 2. The 1st section is the initial 5-min section where every subject stayed stable, the 2nd section is the general road section in which many curves, speed bumps, and stops existed, the 3rd section is general road in which speed bumps and stops did not exist and fast driving was possible, the 4th section is a highway in which the noise constantly exited due to vehicle vibrations caused by faster speed, and the 5th section is national highway in which few curves and stops existed and fast driving was sometimes possible.

Table 2

, , and PIR in each section from the outdoor experiment.

Section	ECG_m		ECG_APA			ECG_RVSS		PIR
	Avg		Avg	d_Avg		Avg	d_Avg
1–5	81.61		85.05	3.45		87.89	6.28	1.82
1	87.12		89.83	2.71		91.79	4.67	1.72
2	82.08		85.99	3.91		88.32	6.23	1.59
3	82.62		84.33	1.71		86.52	3.90	2.28
4	78.11		81.40	3.29		84.88	6.77	2.06
5	81.50		84.12	2.62		86.35	4.85	1.85

The of was compared to that of , and , which is the ECG filtered by the ANC using APA with a large step size of 0.4 in Table 3 to show the tradeoff in step size. Furthermore, , , , and are plotted (Figure 6) to show the actual shape of the ECGs after being filtered.

Table 3

comparisons from the indoor and outdoor experiment.

	Data	SNR^
		ECG_m	ECG_APA	ECG_L_APA	ECG_RVSS
Indoor	1	2.00	2.10	1.40	5.22
	2	2.87	3.41	2.21	7.96
	3	6.78	8.55	5.18	19.67
	4	5.91	6.95	4.27	25.82
	5	4.10	4.90	3.41	11.25
	6	2.58	3.93	3.19	10.14
	7	2.51	3.67	2.88	10.39
	8	1.47	0.90	0.44	0.69
	9	2.83	1.23	0.74	1.56
	10	7.41	9.40	6.51	26.95
	11	2.52	1.05	0.43	0.64
	12	1.21	1.08	0.64	1.13
	13	1.23	0.50	0.20	0.25
	14	2.02	2.35	1.38	3.72
	15	2.33	1.45	1.02	2.03
	16	1.85	1.49	0.88	2.09
	17	3.33	5.67	3.14	17.01
	18	2.86	3.57	2.52	10.90
	19	0.87	1.15	0.69	2.69
	Mean	2.98	3.33	2.16	8.43
Outdoor	1	5.15	2.88	2.81	12.58
	2	3.59	4.68	2.41	7.88
	3	6.49	11.13	4.25	29.91
	4	2.09	0.81	0.52	0.62
	5	2.56	2.79	1.96	5.86
	6	4.04	6.85	3.39	14.41
	7	2.54	3.66	2.29	6.43
	8	5.19	7.26	5.38	20.61
	9	3.64	5.69	3.41	14.50
	10	2.04	2.05	1.11	4.86
	Mean	3.73	4.78	2.75	11.77

Figure 6

Comparison of to , and .

4.3. and Comparison

In general, the are used to evaluate the overall R peak detection performance, as both and are important. The for increased by 8.22 (mean) ± 5.55 (standard deviation)% in the indoor experiment and 6.28 ± 3.92% in the outdoor experiment, whereas the for increased by 5.58 ± 4.42% and 3.45 ± 2.32% in the indoor and outdoor experiment (Table 1). For indoor data 19, the increased at most by 21.49%, and for outdoor data 4, the increased at most by 14.11%. Actually, for several of the indoor and outdoor data sets, of was less than that of , but was much larger, causing the to be larger.

In addition, for was 5.62 ± 5.00% larger than of although for was −0.33 ± 2.92% larger than of in the indoor experiments. Likewise, for was 5.48 ± 3.28% larger than that of although for was −0.20 ± 0.83% larger than that of in outdoor experiments. The reason increased is that the proposed algorithm could filter ECG to reduce the power of the noise within the S-Q interval than the ANC using APA did. Thus, the proposed algorithm could prevent the detection of incorrect R peaks outside the search window, decreasing FP.

For more accurate results analysis (Table 2), the performance increase ratio (PIR) of each section were obtained and shown in the last column where

The two largest PIR were shown in Section 3 and Section 4. In those sections, the subjects could drive fast without curves, speed bumps, and stops so that only trivial and moderate noise occurred due to the vibration of the vehicle caused by fast driving and there was no abnormal noise. As a result, the proposed algorithm could reduce the noise and conserve the QRS interval compared to the ANC using APA. In the other sections, the PIR were lower than or similar to the PIR of the total section. There were curves, speed bumps, stops, and unpaved roads, so abnormal noise occurred. The severe vibration of the vehicle (i.e., driving on unpaved roads) resulted in noise which had statistical characteristics similar to QRS interval of original signal so that the QRS interval could not be distinguished (Figure 7a; left plot). Also, the separation of the body and sensor (i.e., steering, driving over speed bumps) resulted in large noise which caused wrong update of the filter weight (Figure 7; right plot). As a result, neither algorithms could reduce the abnormal noise (Figure 7b,c) and PIR did not increase. In conclusion, the performance depended on the type of noise. The proposed algorithm achieved improved results when moderate noise was mixed, but it showed the limitation that abnormal noise was not reduced.

Figure 7

ECGs with abnormal noise filtered by the APA: (a) with abnormal noise, (b) , and (c) .

4.4. ECGs and Comparison

To investigate how much noise was removed and how much the QRS interval was preserved by filtering, ECGs and were compared. , , , and from outdoor data 2 in were plotted in the same chart (Figure 6). Although all R peaks were detected in the , two incorrect R peaks were detected due to motion noise. In , the motion noise was not removed sufficiently using the small step size even though the QRS-interval power was conserved. In contrast, in , the motion noise was removed using a large step size, but the QRS-interval power was not conserved. Because of these problems, step size adjustment does not contribute significantly to improving SNR, which leads to performance degradation. As a result, the R peak detection algorithm detected wrong R peaks. However, the proposed algorithm removed the motion noise and conserved QRS-interval power. Then, SNR of was larger than that of the other ECGs. As a result, only correct R peaks were detected without incorrect R peaks.

On average, of the obtained from both experiments was larger than that of the , , and (Table 3). For indoor data 8, 11, 13, and outdoor data 4, was smaller than . The quality of the to be filtered was abnormally noisy (Table 1), so the was hard to filter. In addition, of the was larger than that of because reduced QRS-interval power more than noise power. This led to performance degradation so that the for the were, on average, 79.35 ± 15.16% for the indoor experiment and 83.26 ± 5.96% for the outdoor experiment. When APA is used for noise reduction in ECG, its step size is normally selected between 0 and 1. Because there is a significant trade off with conservation of QRS interval and noise reduction, appropriate step size must be selected in the range according to experimental environments; in the 29 data of the both experiments, , on average, resulted in the largest among , …, 1. However, in only two data, s were the largest at , but s were the largest at various size of in the other data. The difference between at and the largest in each data was 0.59 ± 0.70%. However, the step size of the proposed algorithm has a significant role for only noise reduction, so the proposed algorithm with step size near 0.4 which is large enough to remove noise achieved high and improved R peak detection; in the 29 data, , on average, resulted in the largest . Furthermore, in 11 data, s were the largest at . Also, in most data, s were the largest at near 0.4. The difference between at and the largest in each data was 0.25 ± 0.29%. In conclusion, the step size of the ANC using APA must be selected differently depending on the experimental environments for the best performance. However, the best performance of the proposed algorithm is achieved when the step size is near 0.4 regardless of the experimental environments.

5. Conclusions

This study proposed a new ANC structure considering the properties of the ECG and the reference signal obtained from cECG sensors. The weight of the adaptive filter is significantly updated at the QRS interval because the correlation condition is not satisfied, and the original signal is regarded as impulsive noise within that interval. Because filter weight is updated if impulsive noise is present in the QRS interval, ANC using RVSS APA with appropriate parameters was suggested to be robust against the QRS-interval original signal. As a result, the raw ECG is filtered out as the power at the QRS interval is conserved and the power of the noise at the S-Q interval is reduced. of the filtered ECG increased, so improved R peak detection was achieved. Motion noise at S-Q interval was reduced; therefore, the inaccurate detection of R peaks in the raw ECG decreased. In the indoor and outdoor experiments, the proposed algorithm could reduce the noise of the raw ECG more than previous research of ANC using APA. As a result, increased. In addition, the proposed algorithm achieves improved R peak detection when its step size is near 0.4 although there is a significant trade off with conservation of QRS interval and noise reduction in ANC using APA. It is because the proposed algorithm focuses on only noise reduction so step size which is large enough to reduce the noise can achieve improved R peak detection under various experimental environments.

However, the proposed algorithm could not reduce abnormal noise which causes wrong update of the filter weight. Then, the noise reduction performance depended on the types of noise. Furthermore, the proposed algorithm has a negative effect of losing other bioinformation (i.e., P wave, T wave) after adaptive filtering because the proposed algorithm focuses on the conservation of only QRS interval and noise reduction in S-Q interval containing P wave, and T wave for improved R peak detection. This negative effect can result in degraded signal quality. However, P wave and T wave are hardly measured by cECG sensors due to the limitation of the sensors and the measurement noise. Therefore, the proposed algorithm aimed at QRS interval and achieved improved R peak detection.

To improve R peak detection in our system, additional channel data can be used for noise reduction by mounting more sensors into the seat or by using the signal from the sensors already installed. As more channels can share the intervals among the R peaks, the QRS interval can be easily found and the power can be better conserved. More channel data can improve R peak detection by using a channel selection method with the proposed algorithm. The proposed algorithm can be applied to many applications measuring noisy biomedical signal.