An Efficient Teager Energy Operator-Based Automated QRS Complex Detection

Beyramienanlou, Hamed; Lotfivand, Nasser

doi:https://doi.org/10.1155/2018/8360475

Journal of Healthcare Engineering

On this page

Abstract Introduction Results Discussion Conclusion Data Availability Conflicts of Interest References Copyright Related Articles

Research Article | Open Access

Volume 2018 | Article ID 8360475 | https://doi.org/10.1155/2018/8360475

An Efficient Teager Energy Operator-Based Automated QRS Complex Detection

Hamed Beyramienanlou¹and Nasser Lotfivand¹

Academic Editor: Norio Iriguchi

Received03 May 2018

Revised13 Jul 2018

Accepted06 Aug 2018

Published18 Sept 2018

Abstract

Database. The efficiency and robustness of the proposed method has been tested on Fantasia Database (FTD), MIT-BIH Arrhythmia Database (MIT-AD), and MIT-BIH Normal Sinus Rhythm Database (MIT-NSD). Aim. Because of the importance of QRS complex in the diagnosis of cardiovascular diseases, improvement in accuracy of its measurement has been set as a target. The present study provides an algorithm for automatic detection of QRS complex on the ECG signal, with the benefit of energy and reduced impact of noise on the ECG signal. Method. The method is basically based on the Teager energy operator (TEO), which facilitates the detection of the baseline threshold and extracts QRS complex from the ECG signal. Results. The testing of the undertaken method on the Fanatasia Database showed the following results: sensitivity (Se) = 99.971%, positive prediction (P+) = 99.973%, detection error rate (DER) = 0.056%, and accuracy (Acc) = 99.944%. On MIT-AD involvement, Se = 99.74%, P+ = 99.97%, DER = 0.291%, and Acc = 99.71%. On MIT-NSD involvement, Se = 99.878%, P+ = 99.989%, DER = 0.134%, and Acc = 99.867%. Conclusion. Despite the closeness of the recorded peaks which inflicts a constraint in detection of the two consecutive QRS complexes, the proposed method, by applying 4 simple and quick steps, has effectively and reliably detected the QRS complexes which make it suitable for practical purposes and applications.

1. Introduction

Cardiovascular disease is the primary global cause of death. According to the World Health Organization, about 17.3 million people died of cardiovascular disease in 2008, which represented 30 percent of all global deaths. This number is predicted to grow to more than 23.6 million by 2030 [1, 2]. Heart diseases like cardiovascular disease, sudden death, ischemic heart disease, and cardiac arrhythmias are all diagnosed by analyzing the heart’s signal [3–5]. The electrocardiogram is a noninvasive method of recording signals of heart muscle contractions over a period of time. Therefore, the accurate analysis of these signals will result in a more accurate diagnosis of cardiovascular diseases [6, 7]. An ECG signal is a combination of QRS complex, P and T peaks, and sometimes includes U peak. The detection of the QRS complex also helps to detect and determine P and T peaks, QT interval, ST interval, and the respiratory rate, which are considered as the human’s vital signs. Therefore, the accurate recognition of QRS complex has a significant role in the accurate diagnosis of heart disease [8].

During the past decades, a variety of QRS complex detection methods have been developed [9] such as Pan–Tompkins method of R-wave detection [10], support vector machine [11], and the wavelet method, which is an analytical technique based on time-frequency chromatography. The wavelet transform is widely used in medical signal analysis such as EEG or ECG. However, this has a drawback because by applying a fixed scale, the signal characteristics are ignored [3, 12, 13]. Kalman filters use a dynamic model derived from a dynamic system to predict the hidden state in a nonlinear approach [14]. Artificial neural networks are an ideal self-correcting nonlinear process used in a wide range of tasks [15]. Shannon energy computes the average signal energy in a signal spectrum. In other words, it reduces the high intensity to balance out with the low intensity [16–18]. The Adaptive Double Threshold Filter (ADTF) and Discrete Wavelet Transform (DWT) are used to reduce the noise in the ECG signal to improve the ECG signal filtering [19]. Hermit transformation, which is used as an alternative to the Fourier transformation, may by optimization, shows an improved performance [20]. Teager Energy Operator (TEO) mainly shows the frequency and instantaneous changes of the signal amplitude that is very sensitive to subtle changes. Although TEO was first proposed for modeling nonlinear speech signals, it was later widely applied in the audio signal processing. Using TEO can minimize the effects of P and T waves on QRS complex detection [21]. Remarkable research efforts have been developed to analyze the sensitive point of the ECG signal based on TEO [22, 23].

The aim of this study is to propose a new approach based on an innovative viewpoint using TEO to detect QRS complex in the ECG signal. The recorded ECG signal may be affected by noise interference, such as power line interference, which must be eradicated for more accuracy. The Section 2 of the study consists of a series of preprocessing measures to minimize the noises before QRS complex detection on the ECG signal. This includes low-pass filter which removes noises such as power line interference. Since P or T peaks may interfere with the TEO computation, the moving average technique is used to smoothen and envelope the spikes in the signal. Sensitivity, positive prediction, and accuracy of the proposed algorithm from Fantasia Database, MIT-AD, and MIT-NSD are evaluated in the Section 3 of this article. Finally, in Sections 4 and 5 of the study, a discussion and conclusions are presented.

2. Methodology

The details of the proposed method are illustrated in Figure 1. The QRS complex detection procedure involves four steps.(i)Most of the time, the recorded ECG signal is afflicted by noises [24]. The noise frequencies generated by the power lines’ interferences are in the range of 50 to 60 Hz. The noises generated by the muscle contractions and the electrodes placed on the body skin are in the range of 38 to 48 Hz. These greatly impede on the ECG signals. However, a notch filter is very effective in removing these noises. The maximum density of the QRS complex is between 5 to 20 Hz [17, 25]. Therefore, the IIR Butterworth digital filter is the best compromise for phase response and signal attenuation. It has no ripple in the band-pass and is more efficient than the FIR filter [26, 27]. To reduce the noises from the electrical (device) components in order to make a peak clearer for detection, a Butterworth low-pass digital filter, with order 4 and cut-off frequency of 15 Hz, was used.(ii)TEO has various applications, especially in AM and FM signal processing such as speech signals. TEO can be driven from a second-order differential equation [28]. The total energy of oscillation (i.e., the sum of kinetic and potential energies) can be obtained from the following equation: where is the mass of the oscillating body and is the spring constant. Using the formula in (1), a periodic harmonic formula can be obtained: where is the phase shift, is the oscillation frequency, is the oscillation amplitude, and denotes the position of the oscillating body with respect to time. Using (1) and (2), the essential harmonic energy to generate signals can be calculated: The following is a simplified form of TEO: Substituting for we will get the following equation: where is the energy operator for continuous time , is the signal component, and are the first and second derivatives of , respectively, is the sample period, and is the sample size [28, 29]. The dynamicity of the heart beats creates an intermittent and nonlinear pattern for TEO. Since TEO itself is a nonlinear operator, it nonlinearly captivates the intermittent characteristics.(iii)After computing TEO, in some signals, spikes of energies are observable and are attributed to P and T peaks in QRS complex. Although not wide in range, they hamper the accurate detection of the QRS complex. To resolve the situation, these spikes should be converted into energy envelopes. There are several methods for this such as Hilbert transform [6] or averaging method [30]. In the present study, moving average the following equation is used: Here defines a rectangle with length, is a constant and is equal to 1, and defines Teager energy from previous steps. To increase the small amplitude, square root is used: Here is the moving average obtained from the previous step. To decrease the baseline signal’s value below zero, the following formula is used [16]:where is the standard deviation and is defined as signal average.(iv)The process of peak detection includes the following step:where baseline (0) is the threshold level for peak detection and R peaks are found in an ECG signal by searching the maximum peak within ±50 samples (length of window = min (RR interval)) of the recognized location of the candidate R peak in the previous step (Equation (9)).

3. Results

Under the supervision of the National Research Center, the PhysioBank database was developed by the National Institute of Health in order to do a clinical diagnosis and conduct research on complex cardiovascular physiologic signals [31]. The proposed method was tested on three different ECG databases [32] including Fantasia Database (FTD), MIT-BIH Arrhythmia Database (MIT-AD), and MIT-BIH Normal Sinus Rhythm Database (MIT-NSD).

The suggested peak detection method based on TEO has been implemented with MATLAB R2016a on a minimum laptop with a 4 GB of memory and Intel core i3-4000M 2.4 GHz CPU on Windows 10. This algorithm takes less than 0.026 second.

The following formulae were used to determine the performance, sensitivity, error rate detection, positive prediction, and the accuracy of the proposed method:

is the number of R peaks, is the number of missed R peaks, and is the false positive prediction of R peak due to the existing noise with dispositioned true R peak.

3.1. MIT-BIH Arrhythmia Database

MIT-AD contains slightly over 30 minutes of recordings in 48 records. The sampling frequency was set to 360 Hz with 11-bit ADC resolution. The subjects who were chosen for this study were 22 women aged 23 to 89 years and 25 men aged 32 to 89 years [31, 33]. Table 1 depicts the details of detected QRS complex in channel 1. The results showed that sensitivity was at 99.74% with a 0.391% (detection) error and 99.97% positive prediction with an accuracy of 99.71%. Table 2 compares the proposed algorithm in this study with those of other studies. All the stages and the process of QRS detections are illustrated in Figures 2–4. The MIT-AD is available on [36].

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

In Figures 2 and 3, (a) reveals wandering signals. (b) shows that after calculating Teager energy, the amplitudes of the signals are very low and close to zero. Therefore, small values with low energy are reduced to zero, and the wandering signals (drift) are canceled.

3.2. Fantasia Database

Fantasia Database (FTD) contains 40 cases in both groups: the young group aged 21 to 34 years (f1y01 … f2y10 and f2y01 … f2y10) and the elderly group aged 68 to 85 years (f2o01 … f2o10 and f2y01 … f2y10), with an average of 5 men and 5 women in each group. The members of each group underwent 120 minutes of continuous supine resting with complete care. The sampling frequency was set at 250 Hz, with a 16- and 12-bit resolutions for ADC. The records included 2 or 3 channels, such as respiration, ECG signal, and blood pressure [31, 37]. Fantasia Database is available on [38].

The QRS complex detection details for the channel 2 in Fantasia Database are presented in Table 3. Here too, the results showed 99.971% sensitivity with 0.056% detection error, and 99.973% positive prediction with an accuracy of 99.944%. Table 4 shows the comparison of the proposed method with the other studies. Figures 5 and 6 illustrate another example of detection: QRS Complex in the Fantasia Database with both elderly and young subjects. As shown in this figure, the proposed method can remove drift noise and detect correct location beat.

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

3.3. MIT-BIH Normal Sinus Rhythm Database

MIT-NSD contains 18 long-term two-channel ECG recordings. This database includes 5 men, aged 26 to 45 and 13 women, aged 20 to 50. Frequency sampling equals to 128 Hz with 12-bit ADC resolution [31]. The details of QRS complex detection of channel 1 of MIT-NSD is presented in Table 5. The obtained values showed that sensitivity was equal to 99.878%, with an error equal to 0.134, positive prediction was equal to 99.989%, and accuracy was equal to 99.867%. Table 6 includes a comparison of the proposed algorithm with the other studies. Figure 7 illustrates the QRS detection in record with Gaussian white noise. As shown in the figure, the proposed method removed Gaussian white noise, but T peak was detected as a beat. MIT-NSD Database is available on [39].

4. Discussion

The aim of the present research is to use a novel algorithm based on the Teager energy operator in ECG signal to detect QRS complex. The main findings of the study indicated the high reliability and accuracy of this method in QRS complex detection. In spite of applying zero-phase digital filter to maintain QRS complex location, the zero-phase filter is anticausal, and the results showed that the present method had faced a little lag which was less than 0.026 second. Only a detection shift of less than 0.05 second is acceptable [40].

In testing the present method on MIT-AD, some records such as 203 and 210 are main sources of error. The error rate is higher than 1%, which is equal to 0.291. Record 203 has a great number of QRS complexes with multiform ventricular arrhythmia. The TEO phase revealed that the amplitudes are very low and close to zero. Due to this fact, the present method indicated quite a weak performance about records: 203, 19090, and 19830. Records 230, 114, 113, 107, and 106 contain high and sharp T peaks. Record number 207 includes some ventricular flutter (VF) intervals. Those intervals are not interpreted and they go out of studies. One of the constraints of the proposed method is when the QRS complex locations are very close to each other. The length (L) of moving average of phase 2 is assumed to be about . The recorded signals of Fantasia Database included a variety of cardiac morphology, heart failures, and noises from sources like power lines, white Gaussian noise, and flicker noise. Lowering the baseline is the main factor contributing to R-peak losses in the MIT-BIH Normal Sinus Rhythm Database.

The advantages of the proposed method are a reduced number of steps to implement, no need for an excessive memory capacity or learning stage, a fast method of detection, a set baseline threshold, and no complex mathematical relationships.

5. Conclusion

The present study detects QRS complex based on Teager energy, which was tested on four databases. It is a novel algorithm with an acceptable accuracy for ECG baseline prediction. The obtained results from testing the presented method on the Fantasia Database involved: sensitivity (Se) = 99.971%, positive prediction (+P) = 99.973%, detection error rate (DER) = 0.056%, and accuracy (Acc) = 99.944%. On MIT-AD involvement, Se = 99.74%, +P = 99.97%, DER = 0.291%, and Acc = 99.71%. On MIT-NSD involvement, Se = 99.878%, +P = 99.989%, DER = 0.134%, and Acc = 99.867%. The provided results indicate that the presented method is reliable to detect QRS complex, and because the relationships are simple, the proposed method has a better performance than other sophisticated techniques such as neural networks. The results show that the proposed method is simple, effective, accurate, and suitable for practical application. To avoid the lag from zero-phase filter, a low-pass filter and a moving average were used, but still, the signal faced a shift that was about 0.026 s.

Data Availability

The data used to support the findings of this study are available in [36, 38, 39].

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

P. Balakumar, K. Maung-U, and G. Jagadeesh, “Prevalence and prevention of cardiovascular disease and diabetes mellitus,” Pharmacological Research, vol. 113, pp. 600–609, 2016.
View at: Publisher Site | Google Scholar
J. M. de Jesus, S. Kahan, and R. H. Eckel, “Nutrition interventions for cardiovascular disease,” Medical Clinics of North America, vol. 100, no. 6, pp. 1251–1264, 2016.
View at: Publisher Site | Google Scholar
M. Yochum, C. Renaud, and S. Jacquir, “Automatic detection of P, QRS and T patterns in 12 leads ECG signal based on CWT,” Biomedical Signal Processing and Control, vol. 25, pp. 46–52, 2016.
View at: Publisher Site | Google Scholar
F. Jin, L. Sugavaneswaran, S. Krishnan, and V. S. Chauhan, “Quantification of fragmented QRS complex using intrinsic time-scale decomposition,” Biomedical Signal Processing and Control, vol. 31, pp. 513–523, 2017.
View at: Publisher Site | Google Scholar
S. Westerman, N. Engberding, and N. K. Wenger, “Pathophysiology and lifetime risk factors for atherosclerosis and coronary artery disease in women and in the elderly,” in Pathophysiology and Pharmacotherapy of Cardiovascular Disease, G. Jagadeesh, P. Balakumar, and K. Maung-U, Eds., pp. 425–441, Springer International Publishing, Cham, Switzerland, 2015.
View at: Publisher Site | Google Scholar
S. Sahoo, P. Biswal, T. Das, and S. Sabut, “De-noising of ECG signal and QRS detection using hilbert transform and adaptive thresholding,” Procedia Technology, vol. 25, pp. 68–75, 2016.
View at: Publisher Site | Google Scholar
D. Ricciardi, I. Cavallari, A. Creta et al., “Impact of the high-frequency cutoff of bandpass filtering on ECG quality and clinical interpretation: a comparison between 40 Hz and 150 Hz cutoff in a surgical preoperative adult outpatient population,” Journal of Electrocardiology, vol. 49, no. 5, pp. 691–695, 2016.
View at: Publisher Site | Google Scholar
L. D. Sharma and R. K. Sunkaria, “A robust QRS detection using novel pre-processing techniques and kurtosis based enhanced efficiency,” Measurement, vol. 87, pp. 194–204, 2016.
View at: Publisher Site | Google Scholar
A. A. Fedotov, A. S. Akulova, and S. A. Akulov, “Effective QRS-detector based on Hilbert transform and adaptive thresholding,” in IFMBE Proceedings, E. Kyriacou, S. Christofides, and C. S. Pattichis, Eds., vol. 57, pp. 140–144, Springer International Publishing, Cham, Switzerland, 2016.
View at: Publisher Site | Google Scholar
J. Pan and W. J. Tompkins, “A real-time QRS detection algorithm,” IEEE Transactions on Biomedical Engineering, vol. BME-32, no. 3, pp. 230–236, 1985.
View at: Publisher Site | Google Scholar
S. S. Mehta and N. S. Lingayat, “SVM-based algorithm for recognition of QRS complexes in electrocardiogram,” IRBM, vol. 29, no. 5, pp. 310–317, 2008.
View at: Publisher Site | Google Scholar
J. Yan and L. Lu, “Improved Hilbert-Huang transform based weak signal detection methodology and its application on incipient fault diagnosis and ECG signal analysis,” Signal Processing, vol. 98, pp. 74–87, 2014.
View at: Publisher Site | Google Scholar
S. Rekik and N. Ellouze, “Enhanced and optimal algorithm for QRS detection,” IRBM, vol. 38, no. 1, pp. 56–61, 2017.
View at: Publisher Site | Google Scholar
D. Panigrahy and P. K. Sahu, “Extended Kalman smoother with differential evolution technique for denoising of ECG signal,” Australasian Physical & Engineering Sciences in Medicine, vol. 39, no. 3, pp. 783–795, 2016.
View at: Publisher Site | Google Scholar
J. Mateo, A. M. Torres, M. A. García, and J. L. Santos, “Noise removal in electroencephalogram signals using an artificial neural network based on the simultaneous perturbation method,” Neural Computing and Applications, vol. 27, no. 7, pp. 1941–1957, 2016.
View at: Publisher Site | Google Scholar
H. Beyramienanlou and N. Lotfivand, “Shannon’s energy based algorithm in ECG signal processing,” Computational and Mathematical Methods in Medicine, vol. 2017, Article ID 8081361, 16 pages, 2017.
View at: Publisher Site | Google Scholar
M. Rakshit, D. Panigrahy, and P. K. Sahu, “An improved method for R-peak detection by using Shannon energy envelope,” Sadhana—Academy Proceedings in Engineering Sciences, vol. 41, no. 5, pp. 469–477, 2016.
View at: Google Scholar
J. S. Park, S. W. Lee, and U. Park, “R peak detection method using wavelet transform and modified Shannon energy envelope,” Journal of Healthcare Engineering, vol. 2017, Article ID 4901017, 14 pages, 2017.
View at: Publisher Site | Google Scholar
W. Jenkal, R. Latif, A. Toumanari, A. Dliou, O. El B’Charri, and F. M. R. Maoulainine, “An efficient algorithm of ECG signal denoising using the adaptive dual threshold filter and the discrete wavelet transform,” Biocybernetics and Biomedical Engineering, vol. 36, no. 3, pp. 499–508, 2016.
View at: Publisher Site | Google Scholar
M. Brajović, I. Orović, M. Daković, and S. Stanković, “On the parameterization of Hermite transform with application to the compression of QRS complexes,” Signal Processing, vol. 131, pp. 113–119, 2017.
View at: Publisher Site | Google Scholar
A. Zhang, L. Chai, and H. Dong, “QRS complex detection of ECG signal by using teager energy operator,” in Proceedings of 2nd International Conference on Bioinformatics and Biomedical Engineering (iCBBE 2008), pp. 2095–2098, Shanghai, China, May 2008.
View at: Publisher Site | Google Scholar
V. Sharmila and K. A. Reddy, “Identification of premature ventricular cycles of electrocardiogram using discrete cosine transform-Teager energy operator model,” Journal of Medical Engineering, vol. 2015, Article ID 438569, 9 pages, 2015.
View at: Publisher Site | Google Scholar
C. Kamath, “ECG beat classification using features extracted from Teager energy functions in time and frequency domains,” IET Signal Processing, vol. 5, no. 6, p. 575, 2011.
View at: Publisher Site | Google Scholar
M. M. Tantawi, K. Revett, A. B. Salem, and M. F. Tolba, “A wavelet feature extraction method for electrocardiogram (ECG)-based biometric recognition,” Signal, Image and Video Processing, vol. 9, no. 6, pp. 1271–1280, 2015.
View at: Publisher Site | Google Scholar
A. R. Verma and Y. Singh, “Adaptive tunable notch filter for ECG Signal enhancement,” Procedia Computer Science, vol. 57, pp. 332–337, 2015.
View at: Publisher Site | Google Scholar
S. Winder, “Time and frequency response,” in Analog and Digital Filter Design, pp. 41–82, Elsevier, New York City, NY, USA, 2002.
View at: Publisher Site | Google Scholar
S. Winder, “Iir filter design,” in Analog and Digital Filter Design, pp. 395–408, Elsevier, New York, NY, USA, 2002.
View at: Publisher Site | Google Scholar
R. S. Holambe and M. S. Deshpande, “Nonlinear measurement and modeling using Teager energy operator,” in Advances in Non-Linear Modeling for Speech Processing. SpringerBriefs in Electrical and Computer Engineering, pp. 45–59, Springer Science & Business Media, Berlin, Germany, 2012.
View at: Publisher Site | Google Scholar
C. Kamath, “A new approach to detect congestive heart failure using Teager energy nonlinear scatter plot of R-R interval series,” Medical Engineering & Physics, vol. 34, no. 7, pp. 841–848, 2012.
View at: Publisher Site | Google Scholar
A. Felinger, “7 Signal enhancement,” in Data Handling in Science and Technology, vol. 21, pp. 143–181, Elsevier, New York, NY, USA, 1998.
View at: Publisher Site | Google Scholar
A. L. Goldberger, L. A. N. Amaral, L. Glass et al., “PhysioBank, PhysioToolkit, and PhysioNet : components of a new research resource for complex physiologic signals,” Circulation, vol. 101, no. 23, pp. e215–e220, 2000.
View at: Publisher Site | Google Scholar
Physiobank Archieve, “MIT-BIH arrhythmia database,” 2015, https://www.physionet.org/physiobank/database/mitdb/.
View at: Google Scholar
G. B. Moody and R. G. Mark, “The impact of the MIT-BIH arrhythmia database,” IEEE Engineering in Medicine and Biology Magazine, vol. 20, no. 3, pp. 45–50, 2001.
View at: Publisher Site | Google Scholar
F. Zhang and Y. Lian, “QRS detection based on multiscale mathematical morphology for wearable ECG devices in body area networks,” IEEE Transactions on Biomedical Circuits and Systems, vol. 3, no. 4, pp. 220–228, 2009.
View at: Publisher Site | Google Scholar
M. S. Manikandan and K. P. Soman, “A novel method for detecting R-peaks in electrocardiogram (ECG) signal,” Biomedical Signal Processing and Control, vol. 7, no. 2, pp. 118–128, 2012.
View at: Publisher Site | Google Scholar
PhysioNet, “MIT-BIH arrhythmia database,” https://physionet.org/physiobank/database/mitdb/.
View at: Google Scholar
N. Iyengar, C. K. Peng, R. Morin, A. L. Goldberger, and L. A. Lipsitz, “Age-related alterations in the fractal scaling of cardiac interbeat interval dynamics,” American Journal of Physiology-Regulatory, Integrative and Comparative Physiology, vol. 271, no. 4, pp. R1078–R1084, 1996.
View at: Publisher Site | Google Scholar
PhysioNet, “Fantasia database,” https://physionet.org/physiobank/database/fantasia/.
View at: Google Scholar
PhysioNet, “The MIT-BIH normal sinus rhythm database,” https://physionet.org/physiobank/database/nsrdb/.
View at: Google Scholar
A. Amirou, D. Ould-Abdeslam, Z. Zidelmal, M. Aidene, and J. Merckle, “Using S-transform and Shannon energy for electrical disturbances detection,” in Proceedings of IECON 2014-40th Annual Conference of the IEEE Industrial Electronics Society, vol. 116, pp. 2452–2457, Dallas, TX, USA, October 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2018 Hamed Beyramienanlou and Nasser Lotfivand. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

4388

Downloads

1274

Citations