Abstract
Plausibly, the first computerized and automated electrocardiogram (ECG) signal processing algorithm was published in the literature in 1961, and since then, the number of algorithms that have been developed to-date for the detection of the QRS-complexes in ECG signals is countless. Both the digital signal processing and artificial intelligence-based techniques have been tested rigorously in many applications to achieve a high accuracy of the detection of the QRS-complexes in ECG signals. However, since the ECG signals are quasi-periodic in nature, a periodicity analysis-based technique would be an apt approach for the detection its QRS-complexes. Ramanujan filter bank (RFB)-based periodicity estimation technique is used in this research for the identification of the QRS-complexes in ECG signals. An added advantage of the proposed algorithm is that, at the instant of detection of a QRS-complex the algorithm can efficiently indicate whether it is a normal or a premature ventricular contraction or an atrial premature contraction QRS-complex. First, the ECG signal is preprocessed using Butterworth low and highpass filters followed by amplitude normalization. The normalized signal is then passed through a set of Ramanujan filters. Filtered signals from all the filters in the bank are then summed up to obtain a holistic time-domain representation of the ECG signal. Next, a Gaussian-weighted moving average filter is used to smooth the time-period-estimation data. Finally, the QRS-complexes are detected from the smoothed data using a peak-detection-based technique, and the abnormal ones are identified using a period thresholding-based technique. Performance of the proposed algorithm is tested on nine ECG databases (totaling a duration of 48.91 days) and is found to be highly competent compared to that of the state-of-the-art algorithms. To the best of our knowledge, such an RFB-based QRS-complex detection algorithm is reported here for the first time. The proposed algorithm can be adapted for the detection of other ECG waves, and also for the processing of other biomedical signals which exhibit periodic or quasi-periodic nature.
1 Introduction
S. K. Yuen discoursed in his Ph.D. thesis (Yuen, 1976) that the cardiovascular diseases were the leading causes of death in the U.S. in 1972, and even in 2021 (Centers for Disease Control and Prevention, 2022) the statistic has not been altered. Perhaps the scenarios are even more alarming in the other parts of the globe. Therefore, it could be said undoubtedly that the research and development in the area of the computerized processing of the ECG signal and automated detection of cardiovascular diseases are to be carried out in higher pace, if the situation is to alter. An electrocardiogram (ECG) signal characterizes the electrical activity of the heart that arises from the depolarization and repolarization activities of the cardiac cells. ECG signals are recorded by means of electrodes that are placed on various standard locations on the body surface. An ECG signal is composed of waves such as P, QRS-complex and T waves, segments such as PR and ST segments and intervals such as R-to-R and P-to-R intervals, and the presence of any cardiac disorder is generally reflected in the shapes, patterns and durations of these waves, segments and intervals, respectively. Among all other waves and segments of an ECG signal, the QRS-complexes are the most salient regions as it contains the highest frequency spectrum and amplitude. Since the publication of the plausibly-first computerized ECG processing algorithm in 1961 (Pipberger et al., 1961), a multitude of algorithms are proposed in the literature to-date on the detection of the QRS-complexes in ECG signals. A methodological review (Berkayaa et al., 2018) suggests that these QRS-complex detection algorithms could be broadly categorized into two groups: 1) digital signal processing (Burguera, 2019; Hossain et al., 2019; Sharma et al., 2019; Zhang et al., 2020; Jorge et al., 2021; Modak et al., 2021; Morshedlou et al., 2021; Rahul et al., 2021; Tueche et al., 2021), and 2) artificial intelligence-based (Mehta and Lingayat, 2008; Merino et al., 2015; Chandra et al., 2019; Goovaerts et al., 2019; Chen and Maharatna, 2020; Jia et al., 2020; He et al., 2021).
Rahul et al. have proposed an amplitude and interval threshold-based QRS-complex detection algorithm in (Rahul et al., 2021). In (Rahul et al., 2021), first, the ECG signal is denoised using bandpass and moving average filters. Next, a series of static and dynamic amplitude-threshold values (7 threshold values in total) are applied on the denoised ECG signal to identify the plausible QRS-complexes. Though the run-time of the algorithm [run-time indicates the amount of time taken by an algorithm to process a certain amount of data (Cormen et al., 2009)] is low, the QRS-complex detection performance of the algorithm is poor compared to that of the state-of-the-art algorithms (Please see Table 6). A similar type of algorithm is proposed in (Tueche et al., 2021) by Tueche et al. In (Tueche et al., 2021), the ECG signal is first denoised, and then the denoised signal is squared to enhance the QRS-complex regions. Next, all the peaks of the enhanced signal are detected, and three different threshold values are applied on those detected peaks to separate out the real QRS-complexes. The algorithm is primarily implemented on software platform, and then it is also implemented on a microcontroller-based system to examine the likelihood of its success in real-time applications. The authors have claimed that they have not found any difference in the results that they obtained on both the software and hardware platforms. However, here also, as in the case of (Rahul et al., 2021), the QRS-complex detection performance of the algorithm (Tueche et al., 2021) is poor. A tunable-Q wavelet transform (TQWT)-based QRS-complex detection algorithm is proposed by Sharma et al. in (Sharma et al., 2019). A QRS-enhanced ECG signal is obtained from a number of suitably chosen wavelet sub-bands. Then, the correntropy of the QRS-enhanced signal is calculated. Finally, a threshold-based peak detection technique is used to identify the QRS complexes. However, neither the QRS-complex detection performance nor the runtime of the algorithm could compete with that of the state-of-the-art algorithms.
Processing of the ECG signals using deep learning models has become one of the most promising areas of research in the last few years. In their review (Hong et al., 2020), Hong et al. have found 109 published articles, which were focused only on the deep learning-based ECG processing algorithms, and all those articles were published after the year 2019. A U-Net and bidirectional long-short term memory-based QRS-complex detection algorithm is proposed in (He et al., 2021) by He et al. In (He et al., 2021), first, an ECG signal is denoised using a median filter and a wavelet transformation-based technique. The denoised signal is then segmented into blocks of length 10 s each. Those blocks are then fed to the deep learning model to train the classifier. However, even after rigorous training and testing, the algorithm is unable to detect the paced QRS-complexes. Moreover, the model demands high computational-memory, and it require a graphics processing unit to get implemented. Chen et al. have proposed a hierarchical clustering-based R-peak detection and discrete wavelet transform-based T-peak detection algorithm in (Chen and Maharatna, 2020). In (Chen and Maharatna, 2020), the ECG signal is first denoised using Butterworth lowpass and highpass filters to expel the high frequency noises and baseline wonder noise out, respectively, and the amplitude of the denoised signal is normalized within [0, 1]. Then the denoised ECG signal is broken into segments of length 1.2 s each and are processed sequentially. The amplitudes and slopes of all the ECG samples in each of the segment are calculated and then they are clustered into two, namely R-cluster and non-R-cluster. The cluster with a smaller number of samples is considered as the R-cluster, and those samples which belong to R-cluster are considered as the R-peaks. However, a window of length 1.2 s restricts the applicability of the algorithm on ECG signals where the heart rate varies between 50 and 100 beats/minute.
The objective and aim of this proposed research work are to design a high-efficient yet fast QRS-complex detection algorithm overriding the shortcomings of other state-of-the-art algorithms as mentioned above. The main contribution of this research work is that the proposed QRS-complex detection algorithm is 1) highly accurate compared to that of the state-of-the art algorithms, 2) less complex, 3) noise tolerant 4) fast and 5) able to identify the premature ventricular contraction and atrial premature contraction QRS-complexes. The novelty of the proposed algorithm lies both in its working principle and performance. It is worth mentioning that a Ramanujan filter bank-based QRS-complex detection algorithm exploiting the time-period estimation of the ECG signal is being proposed here for the first time, to the best of our knowledge.
The paper is organized as follows. The theoretical background of Ramanujan filter bank is presented in Section 2. Section 3 describes the proposed QRS-complex detection algorithm. The performance of the proposed algorithm is investigated in Section 4. Section 5 compares the performance of the proposed algorithm with that of other state-of-the-art algorithms. Limitation and the future direction of the proposed research work are outlined in Section 6. Finally, discussions and conclusions are drawn in Section 7.
2 Ramanujan Filter Bank
In 1918, the Indian mathematician Srinivasa Ramanujan introduced a new set of sequences, which is known as Ramanujan Sums (Ramanujan, 1918). Ramanujan sums (RS) are real sums defined as the powers of the primitive roots of unity. RS is defined as follows.Where, indicates that and are co-prime, i.e., the greatest common divisor (GCD) is unity, and . The summation is obtained for only those values of that are coprime to , and hence, is periodic, i.e., , and the period is equal to . The RS can also be alternatively calculated as given below.where, is the Euler totient function which is a multiplicative arithmetic function defined for the positive integers and is given by the number of positive co-prime integers less than or equal to . (n), which is known as the Mobius function, is also a multiplicative function and is zero for the positive integers which are not square-free (Sugavaneswaran et al., 2012). From Eq. 2 it can be noted that, ; if , and ; if . Another way of deriving the values of is using the Euler’s formula; . A few Ramanujan-sequences for one period is shown below, where .
It can be seen from Eq. 3 that the values of are integer numbers despite the presence of the trigonometric function as shown in Eq. 1. The identification of the periodicity of a signal using RS has been well studied in (Sugavaneswaran et al., 2012; Chen et al., 2013; Tenneti and Vaidyanathan, 2015a; Tenneti and Vaidyanathan, 2015b; Vaidyanathan and Tenneti, 2015; Pei and Chang, 2017; Tenneti, 2018; Yadav et al., 2018; Saatci and Saatci, 2020). It has also been shown in (Mainardi et al., 2007) that a RS-based technique could better capture the periodicity of a signal better than that of a conventional discrete Fourier transform (DFT)-based method even in the presence of a Gaussian types of noise. The likelihood of application of RS in the context of biomedical signal processing is also reported by Mainardi et al. in (Mainardi et al., 2008).
If the RS is regarded as a digital filter having an impulse response , then its frequency response in is as given below.
From Eq. 4 it can be seen that the value of is non-zero only at the co-prime frequencies where is coprime to ; and zero elsewhere. Let us consider a periodic signal with period . Its Fourier transform can be written as below.
If the signal is passed through Ramanujan filter, by comparing Eqs 4, 5, it can be said that the output of the filter will be zero if neither of the Dirac functions in these two expressions coincide. However, for each value of , such that for , the output of Ramanujan filter could be non-zero. As , Ramanujan filter outputs a non-zero value when is a multiple of . The passband of Ramanujan filter is centered around its co-prime frequency. A collection of such filters for periods ranging from 1 to having impulse response ; ( is the length of the signal ), is called the Ramanujan filter bank (RFB) (Tenneti and Vaidyanathan, 2015a), (Vaidyanathan and Tenneti, 2015), and a plot of the outputs of each of the filters for different values of is referred as the time-period plane plot of the signal.
3 The QRS-Complex Detection Algorithm
The proposed RFB-based QRS-complex detection algorithm can be divided into three major parts: 1) Preprocessing, 2) time-period analysis using RFB and the 3) identification of the QRS-complexes. Schematic diagram of the proposed algorithm is shown in Figure 1.
FIGURE 1
3.1 Preprocessing
The clinical bandwidth of an ECG signal ranges between 0.5 and 100 Hz (Tompkins, 1993). In this research work the original ECG signal is denoised using two sets of filters. In the first set, a 4th order zero-phase Butterworth bandpass filter having lower and upper cut-off frequencies 0.5 and 100 Hz, respectively, and a 4th order zero-phase Butterworth notch filter is used to denoise the original ECG signal. The bandpass filter suppresses the out-of-band noises, and the notch filter suppresses the 50/60 Hz powerline interference noise. This version of the denoised signal is referred as in this paper. In the second set, a 4th order zero-phase Butterworth lowpass filter having a cutoff frequency of 20 Hz is used to denoised the original ECG signal, and this version of the denoised signal is referred as in this paper. The roll-off rate of a Butterworth filter (-20dB/decade) is low compared to that of others such as a Chebyshev filter, but the frequency response of a Butterworth is flat within the passband. This is the main reason for what the Butterworth filters have been used in this research work. A Fourier decomposition-based efficient technique for the removal of baseline wander, powerline interference and other high frequency noises from the ECG signals is proposed in (Singhal et al., 2020) and (Singh et al., 2017). Such a technique could also be used for the removal of high and low frequency noises from the ECG signal.
Next, the amplitudes of both the and signals are normalized in between ±1. Undoubtedly, the signal contains less clinical information as compared to that of its counterpart, but it helps better identifying the QRS-complex regions. Figure 2 shows the example of an original ECG signal, its and versions.
FIGURE 2
3.2 Time-Period Analysis Using RFB
In this step, the is fed to the input of the RFB. There are three parameters that are to be optimized in RFB: 1) the number of filters , 2) the number of times the filter coefficients are to be repeated and 3) the length of the average-filter (Tenneti, 2018), (Tenneti and Vaidyanathan, 2015b). Since, each of the filter in the bank captures different periodic components that are present in the input signal, the number of filters that are to be used are event-dependent. As the goal of this research work is to detect the QRS-complexes of the ECG signals, the number of filters that are to be used in the bank depend on the periodic property of the QRS-complexes. As per the cardio-physiologic definition (Goldberger, 2006), the width of a QRS-complex could stretch 120 milliseconds at most, and therefore, the number of the filters that are to be used has to be calculated based on this clinical information. However, as shown in Figure 3, a QRS-complex is not symmetric itself, i.e., the rise and fall times of a QRS-complex are not identical. Therefore, only the rise-time of a QRS-complex (∼60 milliseconds), i.e., the half of the maximum QRS-duration, is used in this research work to optimize the number of number filters that are required to track the periodicity of a QRS-complex in ECG signal and is calculated as given in Eq. 6.
FIGURE 3
In the proposed research work it has been found experimentally, that repeating the coefficients of each of the Ramanujan filter twice , and a length of the average filter equals to six units provides an optimum QRS-complex detection performance. Next, after repetition, each set of the coefficients are normalized with their corresponding Euclidean norm. Eqs 7, 8 summarize the normalization operation, and Table 1 shows the coefficients of the RFB and their corresponding normalized values after repetition. Two different colours; black and red are used in Table 1 for better showing the repetition.where, is the set of a normalized Ramanujan filter coefficients, is the set of a Ramanujan filter coefficients after repetition and is the Euclidean norm of .
TABLE 1
| — | — | — | — | — | — | |
|---|---|---|---|---|---|---|
| 1 | 1 | — | — | — | — | |
| 0.707 | 0.707 | — | — | — | — | |
| 1 | −1 | 1 | −1 | — | — | |
| 0.5 | −0.5 | 0.5 | −0.5 | — | — | |
| 2 | −1 | −1 | 2 | −1 | −1 | |
| 0.577 | −0.289 | −0.289 | 0.577 | −0.289 | −0.2887 |
The coefficients of the RFB after repetition .
Now, the
signal is separately convoluted with each set of
. The convoluted data is denoted as
. This step filters the
signal using the Ramanujan filter coefficients and tries to capture the periodicity of the QRS-complexes. Then, each of the
is again convoluted with the output of the average-filter. The average filter provides a smoothed version of the filtered data (
). The steps, which are followed in performing the convolution operations are shown below.
for q =1 to
ones((i× ), 1)
end
where, the command ‘ones((
i×
), 1)’ creates a vector of length
i×
, and it contains 1.
is the Euclidean norm of
.
denotes the convoluted
signal using the coefficients of the
Ramanujan filter.
is a matrix of dimension
, where
is the length of
.
is called the periodicity-matrix (
Tenneti and Vaidyanathan, 2015b). For an example, if the sampling rate of the ECG signal is 360 Hz, then, the matrix
would have 21 rows. Therefore, following
Eq. 6, 1) the 21st row of the matrix
would represents the periodicity of that ECG event which completes a period in 0.06 s, 2) the 20th row of the matrix
would represents the periodicity of that ECG event which completes a period in ∼0.055 s, and likewise 3) the 1st row of the matrix
would represents the periodicity of that ECG event which completes a period in ∼0.003 s.
Figure 4shows an
signal and the colour map of its corresponding periodicity-matrix. From this figure it can be seen that the RFB is not only able to capture the periodicity of the QRS-complexes, but it also shows a distinct periodic property for an abnormal QRS-complex.
Figure 5shows an
signal and its corresponding filtered outputs from the different filters in the bank.
FIGURE 4
FIGURE 5
Next the outputs of all the Ramanujan filters are summed up so as to obtain a holistic time-period representation () of the signal. Next, the data is smoothed using a Gaussian-weighted moving average filer, and then the amplitude of the smoothed data is normalized to . The normalized smoothed data is denoted as . Figure 6 shows the example of an signal, its corresponding and .
FIGURE 6
3.3 Identification of the QRS-Complexes
Now, all the peaks of the data are detected and are mapped on to the signal. The peaks on the data indicate the probable R-peaks. As per the cardio-physiologic definition (Goldberger, 2006), the width of a QRS-complex could stretch 120 milliseconds at most. Then, the absolute-maximum amplitude is searched around each of the probable-peaks within a window of length milliseconds to locate the true R-peaks. Since, the signal contains the necessary clinical information, the true R-peaks are detected on the signal. Next, the Q and S peaks are detected using the algorithm, which is used in (Mukhopadhyay and Krishnan, 2020). If it is found that the polarity of the detected R-peak is positive, then the index of the minimum amplitude within a fixed time-window to is considered as the Q-peak, and the index of the minimum amplitude within a fixed time-window and is considered as the S-peak. Instead of a minimum amplitude, the index of the maximum amplitude is searched for detecting both the Q and S-peaks, if the polarity of the detected R-peak is found negative. Here, = index of the R-peak—60 milliseconds, = index of the R-peak, and = index of the R-peak + 60 milliseconds. As mentioned above, the duration of a QRS-complex could be 120 milliseconds at most (Goldberger, 2006), and this is the reason for which a half of the maximum QRS-duration (∼60 milliseconds) is used as the length of the window to look for the Q and S-peaks.
A close observation and analysis, which are listed below, reveal that the
data exhibits distinctly disparate patterns for a normal, premature ventricular contraction (PVC) and premature atrial contraction (APC) types of QRS-complexes.
Figure 7exemplifies these observations.
Observation #1: From Figure 7A, it can be seen that the amplitudes of the peaks and the valleys of the do not differ much for a normal sinus rhythm.
Observation #2: Figure 7B shows that the amplitudes of the peaks of the of the APC beats (the 4th QRS-complex) do not differ much compared to that of the normal QRS-complexes. However, amplitudes of the valleys of the of the APC beats significantly differ from that of the normal QRS-complexes.
Observation #3: Figure 7C shows that the amplitudes of both the peaks and the valleys of the of the PVC beats (the 5th QRS-complex) differs significantly compared to that of the normal QRS-complexes.
FIGURE 7
Following these observations, the APC and PVC beats are segregated using the following steps.
Detect the valleys of the(peaks of theare already detected during the detection of the QRS-complexes)
amplitude of thepeak of,amplitude of thepeak of,amplitude of thevalley of,amplitude of thevalley of
if (()≥and ()≥
theQRS-complexPVC-beat
elseif ((and ()≥
theQRS-complexAPC-beat
else
theQRS-complexnormal QRS-complex
end
After testing the proposed algorithm on a large number of ECG data of various pathological and arrhythmic classes, the value of is set to 0.2.
4 Performance Evaluation
The ECG signals are collected from nine publicly available datasets, totaling a duration of 48.91 days, and are used as the performance evaluation testbed of the proposed algorithm. The details of these collected ECG signals are given in Table 2. The algorithm is implemented and tested on MATLAB platform with a computer having 64-bit Windows 10 operating system, 12 GB RAM and Intel Core i7 10510U central processing unit (CPU) 2.30 GHz. Five statistical metrices: sensitivity (), positive predictivity , detection accuracy , processing time and score are used to assess the performance of the algorithm. These metrices are numerically expressed in Eqs 9–13.where, is the total number of QRS-complexes, which the proposed algorithm has detected correctly. is the total number of QRS-complexes, which the algorithm has failed to detect, and is the total number of detections, which are not the real QRS-complexes, nevertheless, the proposed algorithm has detected those as the real QRS-complexes. In the context of the present research work, i.e., the QRS-complex detection, the sensitivity is a measure of probability that an algorithm efficiently detects the true QRS-complexes among a given set of true QRS-complexes, and the positive predictivity is a measure of probability that an algorithm efficiently detects only the true QRS-complexes among a given set of true and false QRS-complexes (Trevethan, 2017). The close the values of and , and to 100%, and the lower the values of , and , the better the performance of the algorithm. Now, the performance of proposed RFB-based QRS-complex detection algorithm on MITDB is shown in Table 3, and the performance of the proposed algorithm in segregating the normal, APC and PVC beats on MITDB is shown in Table 4 in the form of a confusion matrix. The performance of the proposed algorithm is also tested on eight other datasets, and the results are shown in Table 5.
TABLE 2
| Database | No. of leads used | No. of ECG data files | Sampling frequency in Hz | Duration of ECG in each lead | Total duration in minutes |
|---|---|---|---|---|---|
| European ST-T Database (edb) (Physionet, 2022) | 2 | 89 | 250 | 60 min | 10,680 |
| Fantasia database (fantasia) (Physionet, 2022) | 1 | 37 | 250 | 60 min | 2,220 |
| Dreamer (Katsigiannis and Ramzan, 2018) | 2 | 828 | 256 | 1 min | 1,656 |
| The PhysioNet Computing in Cardiology Challenge 2017 (PCCC) (AF Classification from a Short Single Lead ECG Recording, 2022) | 1 | 8,528 | 300 | Not fixed | 4,619 |
| MIT-BIH Arrhythmia database (MITDB) (Physionet, 2022) | 1 | 48 | 360 | 30 min | 1,440 |
| MIT-BIH Noise Stress Test Database (NSTDB) (Physionet, 2022) | 1 | 12 | 360 | 30 min | 360 |
| The China Physiological Signal Challenge 2018 (CPSC) (Liu et al., 2018) | 12 | 6,877 | 500 | Not fixed | 21,935 |
| Chapman University and Shaoxing People’s Hospital database (CUSPH) (Zheng et al., 2020) | 12 | 10,647 | 500 | 10 s | 21,294 |
| PTB Diagnostic ECG database (PTBDB) (Physionet, 2022) | 12 | 519 | 1,000 | 1 min | 6,228 |
| Total | 48.91 days | ||||
Details of the collected ECG signals.
TABLE 3
| Record # | Total no. of QRS | (%) | (%) | (%) | (%) | in seconds | |||
|---|---|---|---|---|---|---|---|---|---|
| 100 | 2,273 | 2,273 | 0 | 0 | 100 | 100 | 100 | 100 | 0.33 |
| 101 | 1866 | 1864 | 1 | 1 | 99.95 | 99.95 | 99.89 | 99.95 | 0.33 |
| 102 | 2,187 | 2,187 | 0 | 0 | 100 | 100 | 100 | 100 | 0.40 |
| 103 | 2084 | 2084 | 0 | 0 | 100 | 100 | 100 | 100 | 0.33 |
| 104 | 2,229 | 2,229 | 0 | 0 | 100 | 100 | 100 | 100 | 0.33 |
| 105 | 2,572 | 2,567 | 5 | 5 | 99.81 | 99.81 | 99.81 | 99.81 | 0.32 |
| 106 | 2027 | 2027 | 0 | 0 | 100 | 100 | 100 | 100 | 0.32 |
| 107 | 2,135 | 2,135 | 0 | 0 | 100 | 100 | 100 | 100 | 0.34 |
| 108 | 1763 | 1763 | 0 | 4 | 100 | 99.77 | 100 | 99.89 | 0.35 |
| 109 | 2,532 | 2,532 | 0 | 0 | 100 | 100 | 100 | 100 | 0.34 |
| 111 | 2,124 | 2,124 | 0 | 0 | 100 | 100 | 100 | 100 | 0.37 |
| 112 | 2,539 | 2,539 | 0 | 0 | 100 | 100 | 100 | 100 | 0.32 |
| 113 | 1794 | 1794 | 0 | 0 | 100 | 100 | 100 | 100 | 0.33 |
| 114 | 1875 | 1874 | 0 | 1 | 100 | 99.95 | 99.95 | 99.97 | 0.32 |
| 115 | 1953 | 1953 | 0 | 0 | 100 | 100 | 100 | 100 | 0.32 |
| 116 | 2,412 | 2,407 | 6 | 0 | 99.75 | 100 | 99.79 | 99.88 | 0.31 |
| 117 | 1,535 | 1,535 | 0 | 0 | 100 | 100 | 100 | 100 | 0.33 |
| 118 | 2,278 | 2,278 | 0 | 0 | 100 | 100 | 100 | 100 | 0.31 |
| 119 | 1987 | 1987 | 0 | 0 | 100 | 100 | 100 | 100 | 0.32 |
| 121 | 1864 | 1864 | 0 | 0 | 100 | 100 | 100 | 100 | 0.31 |
| 122 | 2,476 | 2,476 | 0 | 0 | 100 | 100 | 100 | 100 | 0.31 |
| 123 | 1,519 | 1,519 | 0 | 0 | 100 | 100 | 100 | 100 | 0.31 |
| 124 | 1,619 | 1,619 | 0 | 0 | 100 | 100 | 100 | 100 | 0.33 |
| 200 | 2,599 | 2,598 | 1 | 0 | 99.96 | 100 | 99.9 | 99.98 | 0.32 |
| 201 | 1978 | 1976 | 3 | 0 | 99.85 | 100 | 99.90 | 99.92 | 0.32 |
| 202 | 2,135 | 2,134 | 1 | 0 | 99.95 | 100 | 99.95 | 99.98 | 0.33 |
| 203 | 2,959 | 2,932 | 27 | 9 | 99.09 | 99.69 | 99.09 | 99.39 | 0.33 |
| 205 | 2,656 | 2,654 | 2 | 0 | 99.92 | 100 | 99.92 | 99.96 | 0.33 |
| 207 | 1862 | 1860 | 2 | 4 | 99.89 | 99.79 | 99.89 | 99.84 | 0.31 |
| 208 | 2,948 | 2,945 | 3 | 1 | 99.90 | 99.97 | 99.90 | 99.93 | 0.35 |
| 209 | 3,005 | 3,005 | 0 | 0 | 100 | 100 | 100 | 100 | 0.32 |
| 210 | 2,651 | 2,651 | 0 | 4 | 100 | 99.85 | 100 | 99.92 | 0.32 |
| 212 | 2,748 | 2,748 | 0 | 0 | 100 | 100 | 100 | 100 | 0.34 |
| 213 | 3,248 | 3,242 | 6 | 0 | 99.82 | 100 | 99.82 | 99.91 | 0.32 |
| 214 | 2,262 | 2,259 | 3 | 0 | 99.87 | 100 | 99.87 | 99.93 | 0.31 |
| 215 | 3,362 | 3,362 | 0 | 0 | 100 | 100 | 100 | 100 | 0.32 |
| 217 | 2,206 | 2,206 | 0 | 0 | 100 | 100 | 100 | 100 | 0.32 |
| 219 | 2,154 | 2,154 | 0 | 0 | 100 | 100 | 100 | 100 | 0.32 |
| 220 | 2047 | 2047 | 0 | 0 | 100 | 100 | 100 | 100 | 0.31 |
| 221 | 2,427 | 2,427 | 0 | 0 | 100 | 100 | 100 | 100 | 0.32 |
| 222 | 2,488 | 2,488 | 0 | 11 | 100 | 99.56 | 100 | 99.78 | 0.33 |
| 223 | 2,605 | 2,605 | 0 | 0 | 100 | 100 | 100 | 100 | 0.37 |
| 228 | 2056 | 2053 | 4 | 8 | 99.81 | 99.61 | 99.85 | 99.71 | 0.32 |
| 230 | 2,256 | 2,256 | 0 | 0 | 100 | 100 | 100 | 100 | 0.32 |
| 231 | 1,571 | 1,571 | 0 | 0 | 100 | 100 | 100 | 100 | 0.31 |
| 232 | 1781 | 1781 | 0 | 3 | 100 | 99.83 | 100 | 99.92 | 0.33 |
| 233 | 3,074 | 3,065 | 9 | 0 | 99.71 | 100 | 99.71 | 99.85 | 0.32 |
| 234 | 2,753 | 2,753 | 0 | 0 | 100 | 100 | 100 | 100 | 0.32 |
| Total | Average | ||||||||
| — | 109,474 | 109,402 | 73 | 51 | 99.94 | 99.95 | 99.94 | 99.95 | 0.33 |
Performance of proposed RFB-based QRS-complex detection algorithm on MITDB.
TABLE 4
| Predicted normal | Predicted APC | Predicted PVC | |
|---|---|---|---|
| Actual Normal | 100,045 | 68 | 37 |
| Actual APC | 5 | 2,503 | 34 |
| Actual PVC | 22 | 127 | 6,854 |
Performance of the proposed algorithm in segregating the normal, APC and PVC beats on MITDB.
The bold values in Table 4 only indicate the diagonal elements of the table.
TABLE 5
| Database | Total number of QRS-complexes | Total | Total | Mean | Mean | Mean | Mean | Mean in seconds |
|---|---|---|---|---|---|---|---|---|
| edb | 858,529 | 228 | 736 | 99.97 | 99.999 | 99.97 | 99.94 | 0.42 |
| Fantasia | 283,747 | 3 | 22 | 99.999 | 99.999 | 99.99 | 99.996 | 0.43 |
| Dreamer | 128,181 | 2 | 16 | 99.999 | 99.99 | 99.998 | 99.99 | 0.01 |
| PCCC 2017 | 575,488 | 152 | 84 | 99.97 | 99.99 | 99.97 | 99.98 | 0.01 |
| NSTDB | 25,596 | 822 | 2,722 | 96.79 | 90.10 | 96.79 | 93.32 | 0.32 |
| CPSC | 1,237,860 | 61 | 58 | 99.995 | 99.995 | 99.995 | 99.995 | 0.01 |
| CUSPH | 2,683,044 | 28 | 47 | 99.999 | 99.998 | 99.999 | 99.999 | 0.01 |
| PTBDB | 795,372 | 17 | 554 | 99.998 | 99.93 | 99.998 | 99.96 | 0.16 |
Performance of the proposed algorithm on other ECG datasets.
The proposed algorithm fails to detect a QRS-complex when its amplitude is too small compared to other such as a few of the QRS-complexes following a large motion artefact noise. On the other hand, the algorithm detects a non-QRS-complex as a QRS-complex when its amplitude is high, and the periodicity matches with that of a QRS-complex. From Table 3 it can be seen that the algorithm is unable to detect 27 QRS-complexes of the file number #203 m of MITDB. The reasons behind this are the presence of: 1) QRS-complexes of different morphologies, and 2) intense motion artifact noise. The ECG record #203 m of MITDB contains QRS-complexes of five different morphologies (normal, aberrated APC, PVC, fusion PVC) and of five different rhythms (normal sinus rhythm, atrial flutter, atrial fibrillation, ventricular trigeminy and ventricular tachycardia) (MIT-BIH Arrhythmia Database Directory, 2022). Such a combination intricates the task of interpretation of an ECG signal even for the cardiologists. On the other hand, the proposed algorithm detects 11 non-QRS-complexes in the ECG record #222. In this ECG record the QRS-complexes are annotated as atrial flutter. A few of those atrial flutters are detected as QRS-complexes. Figures 8–10 show the detected QRS-complexes, APC and PVC beats of the ECG signals of different databases. These figures portray the efficiency of the proposed algorithm in detection the QRS-complexes of different morphologies even in the presence of motion artifact noise.
FIGURE 8
Run-time of an algorithm is a function of the size of the input data file, and it is an important figure of merit. Run-time is used to measure the time complexity of an algorithm. Run-time of the proposed algorithm as a function of the length of the signal at different sampling rates are evaluated and shown in Figure 11. From this figure it is to be noted that the time complexity of the proposed algorithm varies almost linearly with the length of the signal and sampling rate. The proposed algorithm is primarily implemented on software platform, and hence, the space complexity and the power requirement are not considered in this research. However, the run-time of proposed algorithm can further be reduced by efficient code-optimization and implementation.
5 Performance Comparison
The performance of the proposed RFB-based QRS-complex detection algorithm on MITDB is compared with that of state-of-the-algorithms in Table 6. The numerical results of the other algorithms are collected from the respective articles.
TABLE 6
| Algorithms | Total number of beats | (%) | (%) | (%) | Mean in seconds | CPU specification | Superiority of the CPU | |||
|---|---|---|---|---|---|---|---|---|---|---|
| U-Net and LSTM (He et al., 2021) | — | — | — | — | 99.06 | 99.22 | — | 1.8 | CPU 3.5 GHz, and NVIDIA Quadro k6000 GPU, 64 GB of memory | ↟ |
| Box-scoring (Hou et al., 2018) | 110,008 | 748 | 610 | 1,358 | 99.32 | 99.45 | 99.38 | — | Intel Core i7-6700, 3.4 GHz | ↟ |
| Envelope-detection (Burguera, 2019) | 109,985 | — | — | — | 99.57 | 99.37 | — | 1.77 | Intel Core i7, 3.1 GHz | ↟ |
| Pan and Tompkins (Pan and Tompkins, 1985) | 109,809 | 277 | 507 | 784 | 99.75 | 99.54 | 99.64 | — | — | — |
| Amplitude difference (Pandit et al., 2017) | 109,432 | 389 | 369 | 758 | 99.65 | 99.66 | 99.65 | — | 2.5 GHz | CNC |
| Embedded algorithm (Tueche et al., 2021) | 108,791 | 398 | 324 | 722 | 99.65 | 99.69 | 99.67 | 0.31 | Intel Core i7-5500, 2.4 GHz | ↟ |
| Dynamic thresholding (Rahul et al., 2021) | 109,494 | 193 | 155 | 384 | 99.82 | 99.85 | 99.84 | 2.35 | Intel Core i7- 7700, 3.60 GHz | ↟ |
| Deterministic automata (Hamdi et al., 2018) | 73,562 | 284 | 92 | 376 | 99.74 | 99.86 | 99.74 | 0.92 | Intel Core i7-2,600, 3.40 GHz | ↟ |
| Time and amplitude threshold (Modak et al., 2021) | 109,494 | 200 | 136 | 336 | 99.82 | 99.88 | 99.85 | 3.5 | Intel Core i7 CPU, 2.6 GHz | ↟ |
| Wavelet (Sharma et al., 2019) | 109,494 | 183 | 131 | 314 | 99.89 | 99.83 | 99.86 | 69.55 | Intel Core i7, 2.5 GHz | ↟ |
| Adaptive weighting (Sharma and Sharma, 2017) | 109,494 | 113 | 136 | 249 | 99.90 | 99.88 | 99.89 | 1.25 | Intel Core i3 2.50 GHz | ↓ |
| Energy (Yazdani et al., 2018) | 110,070 | 103 | 134 | 237 | 99.91 | 99.88 | 99.89 | — | — | — |
| CNN (Chandra et al., 2019) | 109,494 | 172 | 64 | 236 | 99.84 | 99.95 | 99.89 | — | Intel Core i7, 3.40 GHz | ↟ |
| Interval and trigonometric threshold (Mukhopadhyay and Krishnan, 2020) | 109,497 | 74 | 126 | 200 | 99.94 | 99.88 | 99.91 | 2.43 | Intel Xeon E3-1,225 v3, 3.2 GHz | ↓ |
| Hierarchical Clustering (Chen and Maharatna, 2020) | 109,494 | 124 | 63 | 187 | 99.89 | 99.94 | 99.91 | — | — | — |
| Proposed | 109,474 | 73 | 51 | 124 | 99.94 | 99.95 | 99.95 | 0.33 | Intel Core i7 10510U CPU 2.30 GHz | ↟ |
Performance comparison of the proposed algorithm with a few other recently developed algorithms on MITDB.
↟→superior, ↓→inferior compared to that of the CPU, which is used in this research work, CNC→ cannot be compared.
The bold values in Table 6 indicate the best results in the respective columns.
Before going into a detail comparison, it is good to mention that the performance of all the algorithms, which are included in Table 6, are remarkably high. The , and score (except in (He et al., 2021) and (Burguera, 2019)) of all the algorithms are above 99%. However, there are a few more statistical parameters or figures-of-merit, which are also to be taken into consideration such as , and the run-time . Undoubtedly, the value of depends on the configuration of the system the algorithm is being tested on, and therefore, two additional columns, namely, the CPU specification and the superiority of the CPU, have been included in Table 6 so as to justify the comparison.
The values of (Mukhopadhyay and Krishnan, 2020) and the proposed algorithm are identical, but the and score values of the proposed algorithm are better than that of (Mukhopadhyay and Krishnan, 2020). A direct comparison of between the proposed one and (Mukhopadhyay and Krishnan, 2020) is not possible as the system-configurations are different. In order to make a fair comparison, the algorithm which is proposed in (Mukhopadhyay and Krishnan, 2020) is run on the same CPU, which is used in this research. It has been found that the mean of (Mukhopadhyay and Krishnan, 2020) on MITDB is 0.72 s, when run on the same CPU which is used in this research. It suggests that the proposed algorithm is about 2.18 times faster than (Mukhopadhyay and Krishnan, 2020).
As the number of detection of the proposed algorithm and (Chandra et al., 2019) are low; 51 and 64, respectively, the values of both the algorithms are identical (99.95%). However, the number of detection of (Chandra et al., 2019) is about 2.36 times high than that of the proposed one. The run-time of the algorithm is not mentioned.
The number of and detections of the tunable-Q wavelet transform-based technique (Sharma et al., 2019) are about 2.51 and 2.57 times, respectively, higher than that of the proposed algorithm. Even though, both the proposed algorithm and (Sharma et al., 2019) are implemented on Intel Core i7 processors, the proposed algorithm is about 210.76 times faster than that of (Sharma et al., 2019).
The number of and detections of the U-Net and LSTM-based QRS-complex detection algorithm are not mentioned in (He et al., 2021). The algorithm (He et al., 2021) is implemented on GPU-based system as it requires an intense training operation. Nonetheless, the proposed algorithm is about 5.45 times faster than that of (He et al., 2021).
In Table 6, seven metrics (, , , , and ) are considered to compare the performance of the proposed RFB-based QRS-complex detection technique with that of the others on MITDB. It can be seen from Table 6 that the proposed algorithm: 1) solely holds the best results for four metrics (, and ), and 2) jointly holds the best results for two metrics ( and ).
Now, the QRS-complex detection performance of the proposed RFB-based algorithm on other ECG databases is summarized and compared with a few recently developed algorithms in Table 7. It can be seen that the proposed algorithm performs equally well on other ECG datasets.
TABLE 7
| Algorithms | Total number of beats | Total | Total | Mean GPU, 64 GB of memory | Mean | Mean | Mean in seconds | CPU specification | Superiority of CPU | |
|---|---|---|---|---|---|---|---|---|---|---|
| Dataset: edb | Wavelet (Mourad and Fethi, 2016) | 788,772 | 2,760 | 1,258 | 99.65 | 99.84 | 99.75 | 177 | Intel Core 2 Duo | ↓ |
| Adaptive thresholding (Yakut and Bolat, 2018) | 790,565 | 3,051 | 1,371 | 99.61 | 99.83 | 99.72 | 22 | Intel Core i5-3230, 2.60 GHz | ↓ | |
| Interval and trigonometric threshold (Mukhopadhyay and Krishnan, 2020) | 858,529 | 382 | 1,271 | 99.96 | 99.85 | 99.90 | 9.65 | Intel Xeon E3-1,225 v3, 3.20 GHz | ↓ | |
| Proposed | 858,529 | 228 | 736 | 99.97 | 99.999 | 99.94 | 0.42 | Intel Core i7 10510U CPU 2.30 GHz | ↟ | |
| Dataset: Fantasia | Adaptive thresholding (Yakut and Bolat, 2018) | 283,747 | 166 | 58 | 99.94 | 99.98 | 99.96 | 25.5 | Intel Core i5-3230, 2.60 GHz | ↓ |
| Kurtosis (Sharma and Sunkaria, 2016) | 160,844 | 152 | 148 | 99.90 | 99.91 | 99.91 | — | — | — | |
| Interval and trigonometric threshold (Mukhopadhyay and Krishnan, 2020) | 155,792 | 1 | 197 | 99.999 | 99.87 | 99.94 | 7.15 | Intel Xeon E3-1,225 v3, 3.20 GHz | ↓ | |
| Time and amplitude threshold (Modak et al., 2021) | 285,308 | 224 | 284 | 99.92 | 99.90 | 99.91 | - | Intel Core i7 CPU, 2.6 GHz | ↟ | |
| Proposed | 283,747 | 3 | 22 | 99.999 | 99.999 | 99.996 | 0.43 | Intel Core i7 10510U CPU 2.30 GHz | ↟ | |
| Dataset: NSTDB | Adaptive thresholding (Yakut and Bolat, 2018) | 25,590 | 1,633 | 1,389 | 93.62 | 94.52 | 99.07 | 15.39 | Intel Core i5-3230, 2.60 GHz | ↓ |
| Sixth power (Dohare et al., 2014) | 25,590 | - | - | 88.20 | 89.19 | — | — | Intel Core 2 Duo 2.67 GHz | ↓ | |
| Interval and trigonometric threshold (Mukhopadhyay and Krishnan, 2020) | 25,596 | 1,293 | 3,894 | 94.96 | 87.23 | 90.34 | 2.43 | Intel Xeon E3-1,225 v3, 3.20 GHz | ↓ | |
| Proposed | 25,596 | 822 | 2,722 | 96.79 | 90.10 | 93.32 | 0.32 | Intel Core i7 10510U CPU 2.30 GHz | ↟ | |
| Dataset: PTBDB | Interval and trigonometric threshold (Mukhopadhyay and Krishnan, 2020) | 795,372 | 32 | 954 | 99.996 | 99.88 | 99.94 | 0.06 | Intel Xeon E3-1,225 v3, 3.20 GHz | ↓ |
| Proposed | 795,372 | 17 | 554 | 99.998 | 99.93 | 99.96 | 0.16 | Intel Core i7 10510U CPU 2.30 GHz | ↟ | |
| Dataset: Dreamer | Interval and trigonometric threshold (Mukhopadhyay and Krishnan, 2020) | 128,181 | 2 | 37 | 99.998 | 99.97 | 99.98 | 0.02 | Intel Xeon E3-1,225 v3, 3.20 GHz | ↓ |
| Proposed | 128,181 | 2 | 16 | 99.999 | 99.99 | 99.99 | 0.01 | Intel Core i7 10510U CPU 2.30 GHz | ↟ |
Performance comparison of the proposed algorithm with a few other recently developed algorithms on non-MITDB ECG databases.
6 Limitation and Future Direction
It has been observed that sometimes the proposed algorithm fails to detect a few of the QRS-complexes following a large motion artefact noise when the amplitude of the noise is found to be around 10 times or more than that of the QRS amplitude. However, this limitation could be overcome by clipping the amplitude of such noises to a suitable level at the preprocessing stage. From Figure 9B it can also be seen that some of the noise-peaks of very low amplitude values have been wrongly detected by the algorithm as R-peaks. A suitably chosen threshold value could be applied on the data at the time of detection of the QRS-complexes in order to eliminate these noise-peaks. In our future research we will explore the likelihood of detection of other arrhythmic QRS-complexes using their periodic and aperiodic traits. We will also extend the research work in 1) detecting the p and T waves using their corresponding time-period representations, and 2) processing other periodic/quasi-periodic biomedical signals such as photoplethysmogram, in our future research.
FIGURE 9
7 Conclusion and Discussion
A robust, reliable, fast and high-performance QRS-complex detection algorithm is proposed in this paper. It is worth mentioning that Ramanujan filter bank (RFB) is used in this research work for the detection of the QRS-complexes in ECG signals for the first time, to the best of our knowledge. The advantages of using RFB in detecting the QRS-complexes are multifold. First, the number of required filters in the bank and hence the filter-coefficients depends on the periodic property of the QRS-complexes. Since the periodic property of the QRS-complexes are well defined in the domain of cardio physiology, the required number of filters can easily be derived from the domain knowledge. In this research work it has been found that 1) the number of filters, which is equal to the sampling frequency times the half of the QRS duration; as defined in Eq. 6, and 2) the length of the average filter of unit 6; as described in Section 3.2, best suit in detecting the QRS-complexes in ECG signals. Therefore, the task of normalization of the sampling frequency, which is a prime prerequisite in many of the QRS-complex detection algorithms, could be avoided. Second, as described in Section 3.2, the main mathematical computation that is required to perform is the linear convolution operation only, and therefore the use of a mathematically as well as computationally intricated transformation technique could be avoided. Third, the runtime of the proposed algorithm is incredibly low. Fourth, the effect of the motion artifact types of noises on RFB is extremely low. Many of the QRS-complex detection algorithms detects the high amplitude motion artifact types of noises as the plausible QRS-complexes. This is because of the fact that both the frequency spectra and amplitude of the motion artifact noise and QRS-complex overlap. However, a motion artifact noise is not periodic in nature, and its periodic property does not match with that of the QRS-complexes. This is the main reason for what the false positive detection rate of the proposed algorithm low.
As per the cardio physiologic definition, the heights, frequencies, widths and also the periodicities of different arrhythmic QRS-complexes are different. Figures 7–10 manifest and corroborate this notion. A time-period representation clearly portrays and segregates a normal, premature ventricular contraction and atrial premature contraction QRS-complexes with their corresponding unique periodic patterns. In this research work, alongside detecting the normal QRS-complexes, two types of arrhythmic QRS-complexes (premature ventricular contraction and atrial premature contraction) are identified using a period-thresholding-based technique, and this has been done without the intervention of any ad hoc arrhythmic QRS-complex detection technique.
FIGURE 10
FIGURE 11
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://physionet.org/about/database/.
Author contributions
Both the authors; SM and SK, have equal contribution in this proposed research work and also in writing the paper.
Acknowledgments
We would like to thank the funding provided by the Natural Sciences and Engineering Research Council (NSERC) Discovery Grant program of Canada.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AF Classification from a Short Single Lead Ecg Recording (2022). AF Classification from a Short Single Lead ECG Recording - the PhysioNet Computing in Cardiology Challenge 2017 v1.0.0. [Online]. Available at: https://physionet.org/content/challenge-2017/1.0.0/(Accessed April, 2022).
2
BerkayaaS. K.UysalA. K.GunalE. S.ErgincS.GunalS.GulmezogluM. B. (2018). A Survey on ECG Analysis. Biomed. Signal Process. Control43 (5), 216–235.
3
BurgueraA. (2019). Fast QRS Detection and ECG Compression Based on Signal Structural Analysis. IEEE J. Biomed. Health Inf.23 (1), 123–131. 10.1109/jbhi.2018.2792404
4
Centers for Disease Control and Prevention (2022). Centers for Disease Control and Prevention. [Online]. Available at: https://www.cdc.gov/nchs/fastats/leading-causes-of-death.htm (Accessed April, 2022).
5
ChandraB. S.SastryC. S.JanaS. (2019). Robust Heartbeat Detection from Multimodal Data via CNN-Based Generalizable Information Fusion. IEEE Trans. Biomed. Eng.66 (3), 710–717. 10.1109/tbme.2018.2854899
6
ChenG.KrishnanS.BuiT. D. (2013). Matrix-Based Ramanujan-Sums Transforms. IEEE Signal Process. Lett.20 (10), 941–944. 10.1109/lsp.2013.2273973
7
ChenH.MaharatnaK. (2020). An Automatic R and T Peak Detection Method Based on the Combination of Hierarchical Clustering and Discrete Wavelet Transform. IEEE J. Biomed. Health Inf.24 (10), 2825–2832. 10.1109/jbhi.2020.2973982
8
CormenT. H.LeisersonC. E.RivestR. L.SteinC. (2009). Introduction to Algorithms. Third edition July. Cambridge, MA, USA: The MIT Press.
9
DohareA. K.KumarV.KumarR. (2014). An Efficient New Method for the Detection of QRS in Electrocardiogram. Comput. Electr. Eng.40 (5), 1717–1730. 10.1016/j.compeleceng.2013.11.004
10
GoldbergerA. L. (2006). Clinical Electrocardiography, a Simplified Approach. 7th ed.Amsterdam, Netherlands: Elsevier.
11
GoovaertsG.PadhyS.VandenberkB.VaronC.WillemsR.Van HuffelS. (2019). A Machine-Learning Approach for Detection and Quantification of QRS Fragmentation. IEEE J. Biomed. Health Inf.23 (5), 1980–1989. 10.1109/JBHI.2018.2878492
12
HamdiS.AbdallahA. B.BedouiM. H. (2018). A Robust QRS Complex Detection Using Regular Grammar and Deterministic Automata. Biomed. Signal Process. Control40 (2), 263–274. 10.1016/j.bspc.2017.09.032
13
HeR.LiuY.WangK.ZhaoN.YuanY.LiQ.et al (2021). Automatic Detection of QRS Complexes Using Dual Channels Based on U-Net and Bidirectional Long Short-Term Memory. IEEE J. Biomed. Health Inf.25 (4), 1052–1061. 10.1109/jbhi.2020.3018563
14
HongS.ZhouY.ShangJ.XiaoC.SunJ. (2020). Opportunities and Challenges of Deep Learning Methods for Electrocardiogram Data: A Systematic Review. Comput. Biol. Med.122 (7), 103801–103814. 10.1016/j.compbiomed.2020.103801
15
HossainM. B.BasharS. K.WalkeyA. J.McManusD. D.ChonK. H. (2019). An Accurate QRS Complex and P Wave Detection in ECG Signals Using Complete Ensemble Empirical Mode Decomposition with Adaptive Noise Approach. IEEE Access7, 128869–128880. 10.1109/access.2019.2939943
16
HouZ.DongY.XiangJ.LiX.YangB. (2018). A Real-Time QRS Detection Method Based on Phase Portraits and Box-Scoring Calculation. IEEE Sensors J.18 (9), 3694–3702. 10.1109/jsen.2018.2812792
17
JiaM.LiF.WuJ.ChenZ.PuY. (2020). Robust QRS Detection Using High-Resolution Wavelet Packet Decomposition and Time-Attention Convolutional Neural Network. IEEE Access8, 16979–16988.
18
JorgeR. R.DamasI. D. L.BilaJ.ŠkvorJ. (2021). Internet of Things-Assisted Architecture for QRS Complex Detection in Real Time. Internet Things14 (2), 1–15. 10.1016/j.iot.2021.100395
19
KatsigiannisS.RamzanN. (2018). DREAMER: A Database for Emotion Recognition through EEG and ECG Signals from Wireless Low-Cost Off-The-Shelf Devices. IEEE J. Biomed. Health Inf.22 (1), 98–107. 10.1109/jbhi.2017.2688239
20
LiuF.LiuC.ZhaoL.ZhangX.WuX.XuX.et al (2018). An Open Access Database for Evaluating the Algorithms of Electrocardiogram Rhythm and Morphology Abnormality Detection. J. Med. imaging hlth Inf.8 (7), 1368–1373. 10.1166/jmihi.2018.2442
21
MainardiL.BertinelliM.SassiR. (2008). Analysis of T-Wave Alternans Using the Ramanujan Transform. Comput. Cardiol.14-17, 605–608. 10.1109/cic.2008.4749114
22
MainardiL. T.PattiniL.CeruttiS. (2007). Application of the Ramanujan Fourier Transform for the Analysis of Secondary Structure Content in Amino Acid Sequences. Methods Inf. Med.46 (2), 126–129.
23
MehtaS. S.LingayatN. S. (2008). SVM-based Algorithm for Recognition of QRS Complexes in Electrocardiogram. IRBM29 (5), 310–317. 10.1016/j.rbmret.2008.03.006
24
MerinoM.GómezI. M.MolinaA. J. (2015). Envelopment Filter and K-Means for the Detection of QRS Waveforms in Electrocardiogram. Med. Eng. Phys.37 (6), 605–609. 10.1016/j.medengphy.2015.03.019
25
MIT-BIH Arrhythmia Database Directory (2022). MIT-BIH Arrhythmia Database Directory. [Online]. Available at: https://archive.physionet.org/physiobank/database/html/mitdbdir/records.htm#203 (Accessed April, 2022).
26
ModakS.TahaL. Y.Abdel-RaheemE. (2021). A Novel Method of QRS Detection Using Time and Amplitude Thresholds with Statistical False Peak Elimination. IEEE Access9, 46079–46092. 10.1109/access.2021.3067179
27
MorshedlouF.RavanshadN.DehsorkhH. R. (2021). An Ultra-low Power Analog QRS-Detection Circuit for Ambulatory ECG Monitoring. Int. J. Electron. Commun.129 (2), 1–11. 10.1016/j.aeue.2020.153551
28
MouradK.FethiB. R. (2016). Efficient Automatic Detection of QRS Complexes in ECG Signal Based on Reverse Biorthogonal Wavelet Decomposition and Nonlinear Filtering. Measurement94 (18), 663–670. 10.1016/j.measurement.2016.09.014
29
MukhopadhyayS. K.KrishnanS. (2020). Robust Identification of QRS-Complexes in Electrocardiogram Signals Using a Combination of Interval and Trigonometric Threshold Values. Biomed. Signal Process. Control61 (7), 1–12. 10.1016/j.bspc.2020.102007
30
PanJ.TompkinsW. J. (1985). A Real-Time QRS Detection Algorithm. IEEE Trans. Biomed. Eng.BME-32 (3), 230–236. 10.1109/tbme.1985.325532
31
PanditD.ZhangL.LiuC.ChattopadhyayS.AslamN.LimC. P. (2017). A Lightweight QRS Detector for Single Lead ECG Signals Using a Max-Min Difference Algorithm. Comput. Methods Programs Biomed.144 (7), 61–75. 10.1016/j.cmpb.2017.02.028
32
PeiS.-C.ChangK.-W. (2017). Two-Dimensional Period Estimation by Ramanujan's Sum. IEEE Trans. Signal Process.65 (19), 5108–5120. 10.1109/tsp.2017.2726986
33
Physionet (2022). Physionet. [Online]. Available at: https://physionet.org/about/database/(Accessed April, 2022).
34
PipbergerH. V.ArmsR. J.StallmannF. W. (1961). Automatic Screening of Normal and Abnormal Electrocardiograms by Means of a Digital Electronic Computer. Exp. Biol. Med.106 (1), 130–132. 10.3181/00379727-106-26260
35
RahulJ.SoraM.SharmaL. D. (2021). Dynamic Thresholding Based Efficient QRS Complex Detection with Low Computational Overhead. Biomed. Signal Process. Control67 (5), 1–16. 10.1016/j.bspc.2021.102519
36
RamanujanS. (1918). On Certain Trigonometrical Sums and Their Applications in the Theory of Numbers. Trans. Camb. Philosophical Soc.22 (13), 259–276.
37
SaatciE.SaatciE. (2020). Period Determination in Cyclo-Stationary Signals by Autocorrelation and Ramanujan Subspaces. IEEE Signal Process. Lett.27, 266–270. 10.1109/lsp.2020.2966877
38
SharmaA.PatidarS.UpadhyayA.Rajendra AcharyaU. (2019). Accurate Tunable-Q Wavelet Transform Based Method for QRS Complex Detection. Comput. Electr. Eng.75 (3), 101–111. 10.1016/j.compeleceng.2019.01.025
39
SharmaL. D.SunkariaR. K. (2016). A Robust QRS Detection Using Novel Pre-processing Techniques and Kurtosis Based Enhanced Efficiency. Measurement87 (11), 194–204. 10.1016/j.measurement.2016.03.015
40
SharmaT.SharmaK. K. (2017). QRS Complex Detection in ECG Signals Using Locally Adaptive Weighted Total Variation Denoising. Comput. Biol. Med.87 (8), 187–199. 10.1016/j.compbiomed.2017.05.027
41
SinghP.JoshiS. D.PatneyR. K.SahaK. (2017). The Fourier Decomposition Method for Nonlinear and Non-stationary Time Series Analysis. Proc. Math. Phys. Eng. Sci.473 (2199), 20160871. 10.1098/rspa.2016.0871
42
SinghalA.SinghP.FatimahB.PachoriR. B. (2020). An Efficient Removal of Power-Line Interference and Baseline Wander from ECG Signals by Employing Fourier Decomposition Technique. Biomed. Signal Process. Control57 (3), 1–8. 10.1016/j.bspc.2019.101741
43
SugavaneswaranL.XieS.UmapathyK.KrishnanS. (2012). Time-Frequency Analysis via Ramanujan Sums. IEEE Signal Process. Lett.19 (6), 352–355. 10.1109/lsp.2012.2194142
44
TennetiS. V. (2018). The Nested Periodic Subspaces: Extensions of Ramanujan Sums for Period Estimation. Ph.D. dissertation (Pasadena, California: Dept. of Electrical Engineering, California Institute of Technology). 10.7907/1n4t-5876
45
TennetiS. V.VaidyanathanP. P. (2015). Nested Periodic Matrices and Dictionaries: New Signal Representations for Period Estimation. IEEE Trans. Signal Process.63 (14), 3736–3750. 10.1109/tsp.2015.2434318
46
TennetiS. V.VaidyanathanP. P. (2015). “Ramanujan Filter Banks for Estimation and Tracking of Periodicities,” in IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), South Brisbane, QLD, Australia, 19-24 April 2015, 3851–3855. 10.1109/icassp.2015.7178692
47
TompkinsW. J. (1993). Biomedical Digital Signal Processing. New Jersey, USA: Prentice-Hall.
48
TrevethanR. (2017). Sensitivity, Specificity, and Predictive Values: Foundations, Pliabilities, and Pitfalls in Research and Practice. Front. Public Health5, 307–7. 10.3389/fpubh.2017.00307
49
TuecheF.MohamadouY.DjeukamA.KouekeuL. C. N.SeujipR.TonkaM. (2021). Embedded Algorithm for QRS Detection Based on Signal Shape. IEEE Trans. Instrum. Meas.70, 1–12. 10.1109/tim.2021.3051412
50
VaidyanathanP. P.TennetiS. V. (2015). “Properties of Ramanujan Filter Banks,” in 23rd European Signal Processing Conference (EUSIPCO), Nice, France, 31 Aug.- 4 Sept. 2015, 2816–2820. 10.1109/eusipco.2015.7362898
51
YadavD. K.KuldeepG.JoshiS. D. (2018). Orthogonal Ramanujan Sums, its Properties, and Applications in Multiresolution Analysis. IEEE Trans. Signal Process.66 (21), 5789–5798. 10.1109/tsp.2018.2871384
52
YakutÖ.BolatE. D. (2018). An Improved QRS Complex Detection Method Having Low Computational Load. Biomed. Signal Process. Control42 (4), 230–241. 10.1016/j.bspc.2018.02.004
53
YazdaniS.FalletS.VesinJ.-M. (2018). A Novel Short-Term Event Extraction Algorithm for Biomedical Signals. IEEE Trans. Biomed. Eng.65 (4), 754–762. 10.1109/tbme.2017.2718179
54
YuenS. K. (1976). Cardalert: A Portable, Battery-Operated, Real-Time Arrhythmia Detector and Alarm System. Ph.D. dissertation (Champaign, IL, USA: Dept. of Computer Science, University of Illinois at Urbana-Champaign).
55
ZhangZ.YuQ.ZhangQ.NingN.LiJ. (2020). A Kalman Filtering Based Adaptive Threshold Algorithm for QRS Complex Detection. Biomed. Signal Process. Control58 (4), 1–8. 10.1016/j.bspc.2019.101827
56
ZhengJ.ZhangJ.DaniokoS.YaoH.GuoH.RakovskiC. (2020). A 12-lead Electrocardiogram Database for Arrhythmia Research Covering More Than 10,000 Patients. Sci. Data7, 48–8. 10.1038/s41597-020-0386-x
Summary
Keywords
electrocardiogram, period estimation, Ramanujan filter bank, Ramanujan sums, QRS-complex detection
Citation
Mukhopadhyay SK and Krishnan S (2022) Robust Identification of the QRS-Complexes in Electrocardiogram Signals Using Ramanujan Filter Bank-Based Periodicity Estimation Technique. Front. Sig. Proc. 2:921973. doi: 10.3389/frsip.2022.921973
Received
17 April 2022
Accepted
24 May 2022
Published
29 June 2022
Volume
2 - 2022
Edited by
Manjeet Kumar, Delhi Technological University, India
Reviewed by
Pushpendra Singh, National Institute of Technology, Hamirpur, India
Muammar Muhammad Kabir, University Health Network, Canada
Updates
Copyright
© 2022 Mukhopadhyay and Krishnan.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sourav Kumar Mukhopadhyay, sourav.mukhopadhyay@ryerson.ca
This article was submitted to Biomedical Signal Processing, a section of the journal Frontiers in Signal Processing
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.