Journal of Medical Engineering

Journal of Medical Engineering / 2016 / Article

Research Article | Open Access

Volume 2016 |Article ID 9823026 |

S. Rangsungnoen, P. Chanbenjapipu, N. Mathuradavong, K. Suwanprasert, "A Hybrid Signal Processing of RR Intervals from QTc Variation Searching Arrhythmia and Improving Heart Rate Variability Assessment in Acute Large Artery Ischemic Stroke", Journal of Medical Engineering, vol. 2016, Article ID 9823026, 6 pages, 2016.

A Hybrid Signal Processing of RR Intervals from QTc Variation Searching Arrhythmia and Improving Heart Rate Variability Assessment in Acute Large Artery Ischemic Stroke

Academic Editor: Nicusor Iftimia
Received30 Jun 2016
Accepted20 Oct 2016
Published14 Nov 2016


Sudden death caused by abnormal QTc and atrial fibrillation (AF) has been reported in stroke. Heart rate variability (HRV) is reduced with missing beats of RRI during arrhythmic episode and abnormal QTc variation during acute stroke. In this study, we develop a hybrid signal processing by Pan Tompkins QRS detection and Kalman filter estimator for meaningful missing beats and searching AF with prolonged QTc. We use this hybrid model to investigate RRIs of Lead II ECG in thirty acute stroke patients with long QTc and AF (LQTc-AF) and normal QTc without AF (NQTc-nonAF) and then assess them by HRV. In LQTc-AF Kalman, higher mean heart rate with lower mean RRIs compared to NQTc-nonAF Kalman was characterized. LQTc-AF Kalman showed significant increase in SDNN, HF, SD2, SD2/SD1, and sample entropy. SDNN and HF associated with high RMSSD, pNN50, and SD1 reflect predominant parasympathetic drive for sympathovagal balance in LQTc-AF Kalman. Greater SD2, SD2/SD1, and sample entropy indicate more scatter of Poincaré plot. Compared with conventional Labchart, fractal scaling exponent of α1 (DFA) is higher in LQTc-AF Kalman. Remarkable complexity with parasympathetic drive in LQTc-AF Kalman suggests an influence of missing beats during stroke.

1. Introduction

Atrial arrhythmia (AF), prolonged QT interval, and tachycardia are the most common cardiac arrhythmia in ischemic stroke with recurrent stroke and sudden death [1, 2]. An inverse correlation between SDNN, parasympathetic HRV parameters (rMSSD and pNN50), and atrial fibrillatory rate (AFR) is reported and it is a possible link between AFR and autonomic modulation [3]. Pathophysiology of silent AF, prolonged or shortened QTc, or even ventricular tachycardia during acute stroke is poorly elucidated. Ventricular repolarization affects QT interval variation and might affect the next P wave on the second ECG cycle. R-R interval is a representative of beat to beat of heart contraction which composes of two cardiac cycles. Hence, abnormal ECG waves behave as any noise or showing missing beat are very interesting which they may influence cardiac autonomic modulation during stroke. Missing beats of R-R interval have been shown contributing abnormal QTc and arrhythmia during acute stroke. Determination of QT correction is essential for R-R interval which can be calculated in Bazett’s method (Figure 1) formula that we called QTc and divided to be short, normal (NQTc), and long QT (LQTc) intervals [4]. During stroke, fluctuation of QTc is obviously found with missing beat of R-R interval. Estimation and computed missing beat of R-R interval by Kalman filter are proposed for better interpretation results with preprocessing to detect QRS complexes by Pan Tomkins algorithm as described in methodology. As we know, brain-heart axis is shifted affecting autonomic modulation in AF patient before and after stroke [3, 58]. These brain-heart interactions will be a new way to explore insights of sudden cardiac event and provide important information for future novel therapies to prevent recurrent arrhythmia and sudden cardiac death [2]. An effective improvement of heart rate variability (HRV) analysis is essential to verify the complexity of autonomic modulation during acute stroke [4, 911]. In this study, we develop a hybrid signal processing by Pan Tompkins QRS detection and Kalman filter estimator for missing beat correction in order to examine proving arrhythmia and improving HRV [1214]. We investigate missing beat behavior in long QTc with AF and normal QTc with nonAF during 24 hours of acute stroke and then assess both groups by HRV analysis. We investigate whether QT interval fluctuation or other silent AF will be examined by hybrid signal processing method and test the difference between LQTc-AF Kalman vs NQTs-nonAF Kalman strokes. Quantitative and qualitative studies of HRV such as sample entropy, Poincare’ plot, and detrended fluctuation analysis (DFA) were analyzed. We hypothesize that developed hybrid signal processing model is fitted with searching silent arrhythmia and improving HRV analysis in QTc variation.

2. Material and Methods

Total of 165 patients were recruited to be studied and 65 cases were characterized to be AF during acute large artery ischemic stroke by visual interpretation from physician. Thirty long QTc with AF (LQTc-AF) and no previous AF history and normal QTc without AF (NQTc-nonAF) stroke patients were investigated in this study (averaged 65 ± 4 years old and matching control). All participants gave informed consent according to local Ethic Committee (MTU-EC-IM-4-018/54). ECG data were digitized directly recording at sampling rate 1000 Hz by PowerLab system monitor. The short-term five-minute HRV was analyzed from ECG signals for linear and nonlinear HRV by using Kubios HRV software (version 2.0, University of Kuopio, Finland) [6]. Thereafter, RR intervals were examined by hybrid signal processing model before they were assessed by autonomic derangement by linear and nonlinear functions of HRV. LQTc-AF group was also examined by conventional Labchart program and compared to a hybrid model. QTc was determined and classified by Bazett’s [4] and Viskin’s methods [5]. QTc calculation was defined by divided distance between the start of the QRS complex and the end of the T wave with the square root of the preceding RR interval.(1)The preceding RR interval can be measured by calculating the distance between the first peak of R wave and the next peak of R wave in beats per minute.(2)The QT interval can be measured by calculating the distance between the start of the QRS complex and the end of the T wave.(3)Divide the QT interval by the square root of the RR interval.

QT Correction (QTc) Measurement. The QT interval was measured at a speed of 25 and 50 mm/s with the PowerLab systems (ADInstruments) and the EKGs were accepted for QT evaluation if the heart rate was between 60 and 100 beats/min. The QT interval was corrected for heart rate using Bazett’s formula (QTc = QT/). Short QTc (SQTc), normal QTc (NQTc), and long QTc (LQTc) were defined as lesser than 360 and 360–390 and greater than 400 milliseconds, respectively [15].

All data in linear and nonlinear functions of HRV were analyzed [9, 10]. For the processing, preprocessing by Pan Tomkins algorithm to detect QRS complexes in EKG was performed [16]. The algorithm is based on analysis of the slope, amplitude, and width of QRS complexes including a series of filters and methods that perform low pass, high pass, derivative, squaring, integration, adaptive threshold, and search procedures. The next step included applying Pan Tomkins for parameter extraction (RR interval) and then applied Kalman filter for RR interval improvement. Kalman filter algorithms have a two-step process [17]: firstly, update step measurement.

Update. The estimate is of valueThe Kalman filter generated estimates of the current state variables, along with their uncertainties. Secondly, the prediction step is as follows: : calculation of Kalman Gain.: estimation from measurement sensor.: the estimate of the error covariance.The outcome of the next measurement is observed; these estimates are updated using a weighted average, with more weight being given to estimates with high certainty as in flow chart (Figure 2).

3. Statistical Analysis

The variables are expressed as mean ± SEM. Data are analyzed by ANOVA. Significant difference is taken at value <0.05.

4. Results and Discussion

Kalman filter algorithm is demonstrated in Figure 2 which determines prediction and correction, described in methodology. Missing beats are evident in LQTc stroke case (dotted line) compared with corrected missing beats by hybrid signal processing method (solid line) as shown in Figure 3. Consistency of RRIs after Kalman estimator suggests an efficient data for further HRV analysis. In Table 1, all values show significant difference between LQTc-AF Kalman stroke and NQTc-nonAF Kalman stroke. Mean RRIs in LQTc-AF Kalman stroke are less than those in NQTc-nonAF Kalman stroke. Significantly greater parasympathetic activity is evident indicating its power drive for sympathovagal balance in LQTc-AF Kalman stroke group as shown by HF, SD2, SD2/SD1, and high value of RMSSD, pNN50, and SD1 as well [1821]. Compared with NQTc-nonAF Kalman stroke, decreasing quantitation of sympathovagal balance is shown by lower LF/HF ratio in LQTc-AF Kalman stroke. An irregularity of signals (SamEn) in LQTc-AF Kalman stroke group is greater than in NQTc-nonAF Kalman stroke as presented by geometry of Poincaré plot in Figure 4 [19]. As in Table 2, significant difference of HRV descriptors between LQTc-AF Kalman and LQTc-AF recorded from conventional Labchart is evident. Significant difference in RMSSD, HF, SD1, SD2/SD1, and fractal scaling exponent of α1 indicates the impact of missing beat and the finding suggests remarkable complexity corresponding with predominant parasympathetic drive (by RMSSD, HF, SD1, and SD2/SD1) [21].

LQTc-AF KalmanNQTc-nonAF Kalman

Mean HR100.30 ± 27.9973.98 ± 12.18
Mean RRI669.3 ± 134.0845.6 ± 154.5
SDNN114.3 ± 68.158.5 ± 41.3
RMSSD93.4 ± 48.334.43 ± 26.68NS
pNN5044.36 ± 26.558.10 ± 11.88NS
VLF (0–0.04 Hz)3001 ± 47291266 ± 1300NS
LF (0.04–0.15 Hz)6404 ± 97151203 ± 1733NS
HF (0.15–0.4 Hz)3893 ± 4022832 ± 1234
LF/HF ratio1.459 ± 0.6182.464 ± 2.501NS
SD166.23 ± 34.2324.45 ± 18.92NS
SD2146.9 ± 91.777.3 ± 54.8
SD2/SD12.249 ± 0.6453.620 ± 2.443
SamEn1.413 ± 0.6550.756 ± 0.460

NS: not significantly different. : significant difference. HR: heart rate (beat/min). RRI: RR interval (msec). SDNN (ms): standard deviation of all NN intervals. RMSSD (ms): square root of the mean squared difference between adjacent NN intervals. PNN50 (%): number of adjacent NN intervals that differ by 50 ms/total number of NN intervals. VLF (ms2, %): very low frequency range (<0.003 Hz). LF (ms2, %): low frequency range (0.04–0.15 Hz). HF (ms2, %): high frequency range (0.15–0.4 Hz). SD1: standard deviation of points perpendicular to the line of identity. SD2: standard deviation of points along to the line of identity. SamEn: sample entropy.

LQTc-AF KalmanLQTc-AF Labchart

Mean HR100.30 ± 27.9985.97 ± 16.74NS
Mean RRI669.3 ± 134.0751 ± 139.3NS
SDNN114.3 ± 68.1148.7 ± 126.3NS
RMSSD93.4 ± 48.3191.9 ± 154.3
pNN5044.36 ± 26.5556.64 ± 34.96NS
VLF (0–0.04 Hz)3001 ± 47298852 ± 27625NS
LF (0.04–0.15 Hz)6404 ± 97154567 ± 3878NS
HF (0.15–0.4 Hz)3893 ± 40228549 ± 7725
LF/HF ratio1.459 ± 0.6180.685 ± 0.412NS
SD166.23 ± 34.23136.2 ± 109.5
SD2146.9 ± 91.7159 ± 144.2NS
SD2/SD12.249 ± 0.6451.470 ± 0.704
SamEn1.413 ± 0.6551.654 ± 0.474NS
α11.0822 ± 0.16910.7714 ± 0.2563
α20.8050 ± 0.2340.6466 ± 0.2871NS

DFA is more fluctuated in LQTc-AF group without improved RR interval (Figure 5(a)) compared with those in LQTc-AF group with improving RR interval by Kalman filter (Figure 5(b)). More complexity as shown by DFA in LQTc-AF Labchart group (Figure 5(a)) may be from other noise or missing beats in ECG tracing and is likely lead to poor result.

AF is a complex arrhythmia related with pulmonary vein (PV) myocardial sleeves which originally generate either reentrant or focal electrical signals. They receive parasympathetic innervation through small branches of vagus nerve [22]. It is well recognized that autonomic control is an important modulator in mediating AF [23]. Vagal stimulation causes several effects such as reducing conduction velocity in atrial tissue, shortening action potential duration and atrial effective refractory periods (ERPs), and facilitating reentrant pathway [24]. Previously, rising atrial fibrillatory rate (AFR) is inversely correlated with declining pNN50 that represents vagal outflow in acute AF stroke patient [3]. Oliveira et al. have reported that vagal activity prolongs interatrial conduction and shortens atrial and PVERP, contributing to induction and maintenance of AF in rabbit [25]. In this study, predominant parasympathetic drive in LQTc-AF Kalman may promote PV myocardial sleeves activity causing AF in arrhythmia episode in stroke. The QT interval reflects two phases of action potential in each cardiac cycle: depolarization of atrium and ventricle and repolarization of ventricle [4]. Any factors that augment depolarization or delay repolarization of the myocardium are able to increase QT interval length [5, 15]. Changes in T wave duration (ventricular repolarization) will affect the next P wave of consecutive cycle. In normal person, vagal activity reduces heart rate per minute and subsequently prolongs RRI [26]. During acute stroke, the modulation of brain-heart axis with predominant parasympathetic activity will facilitate long QTc. In summary, detection of missing beats and correction of RRI by a hybrid processing model are able to improve HRV analysis and determine AF in arrhythmia episode during acute stroke.

5. Conclusions

This finding suggests that a hybrid signal processing that we have developed for better RR detection by estimation of missing beats is essential in QTc variation and arrhythmia episode in stroke. Greater predominant parasympathetic activity in LQTc-AF Kalman stroke plays a vital role in sympathovagal balance of an inconsistency and irregularity signals [13, 14, 27]. High values of sample entropy corresponding with more scatter gram of Poincaré plot indicate signal turbulence and more variation in QTc with AF stroke. Compared with conventional processing, such as Labchart, our developed hybrid model fits searching silent AF and improving HRV, especially nonlinear function. Searching AF during shortened QTc episode will be further studied since it is one major cause of AF episode and acute stroke. Early detection of consistent AF with abnormal QTc variation by our hybrid model will be useful for physician to examine the stroke treatment and recurrent.

Competing Interests

The authors declare that there are no competing interests regarding the publication and the funding that they have received.


This work was supported by National Research University Project of Thailand Office of Higher Education Commission and Faculty of Medicine, Thammasat University.


  1. J. T. Korpelainen, K. A. Sotaniemi, A. Mäkikallio, H. V. Huikuri, and V. V. Myllylä, “Dynamic behavior of heart rate in ischemic stroke,” Stroke, vol. 30, no. 5, pp. 1008–1013, 1999. View at: Publisher Site | Google Scholar
  2. F. Colivicchi, A. Bassi, M. Santini, and C. Caltagirone, “Cardiac autonomic derangement and arrhythmias in right-sided stroke with insular involvement,” Stroke, vol. 35, no. 9, pp. 2094–2098, 2004. View at: Publisher Site | Google Scholar
  3. R. Sungnoon, K. Suwanprasert, and S. Muengtaweepongsa, “Atrial electrophysiological property analysis by sample entropy and atrial fibrillatory rate with cardiac autonomic derangements in acute ischemic stroke with atrial fibrillation,” Neurology Asia, vol. 19, no. 1, pp. 11–18, 2014. View at: Google Scholar
  4. H. C. Bazett, “An analysis of the time—relations of electrocardiograms,” Heart, vol. 7, pp. 353–370, 1920. View at: Google Scholar
  5. S. Viskin, “The QT interval: too long, too short or just right,” Heart Rhythm, vol. 6, no. 5, pp. 711–715, 2009. View at: Publisher Site | Google Scholar
  6. M. P. Tarvainen, J.-P. Niskanen, J. A. Lipponen, P. O. Ranta-aho, and P. A. Karjalainen, “Kubios HRV—heart rate variability analysis software,” Computer Methods and Programs in Biomedicine, vol. 113, no. 1, pp. 210–220, 2014. View at: Publisher Site | Google Scholar
  7. M. Galinier, J.-C. Vialette, J. Fourcade et al., “QT interval dispersion as a predictor of arrhythmic events in congestive heart failure,” European Heart Journal, vol. 19, no. 7, pp. 1054–1062, 1998. View at: Publisher Site | Google Scholar
  8. M. Nagai, S. Hoshide, and K. Kario, “The insular cortex and cardiovascular system: a new insight into the brain-heart axis,” Journal of the American Society of Hypertension, vol. 4, no. 4, pp. 174–182, 2010. View at: Publisher Site | Google Scholar
  9. P. Taggart, H. Critchley, and P. D. Lambiase, “Heartebrain interactions in cardiac arrhythmia,” Heart, vol. 97, no. 9, pp. 698–708, 2011. View at: Publisher Site | Google Scholar
  10. C.-F. Chen, C.-L. Lai, H.-F. Lin, L.-M. Liou, and R.-T. Lin, “Reappraisal of heart rate variability in acute ischemic stroke,” Kaohsiung Journal of Medical Sciences, vol. 27, no. 6, pp. 215–221, 2011. View at: Publisher Site | Google Scholar
  11. K. Intharakham and K. Suwanprasert, “Complexity of autonomic control during cerebrovascular reactivity,” in Proceedings of the 6th Biomedical Engineering International Conference (BMEiCON '13), Krabi, Thailand, October 2013. View at: Publisher Site | Google Scholar
  12. A. E. Aubert, S. Vandeput, F. Beckers, J. Liu, B. Verheyden, and S. Van Huffel, “Complexity of cardiovascular regulation in small animals,” Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 367, no. 1892, pp. 1239–1250, 2009. View at: Publisher Site | Google Scholar
  13. S. De Raedt, A. De Vos, and J. De Keyser, “Autonomic dysfunction in acute ischemic stroke: an underexplored therapeutic area?” Journal of the Neurological Sciences, vol. 348, no. 1-2, pp. 24–34, 2015. View at: Publisher Site | Google Scholar
  14. Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology, “Heart rate variability: standards of measurement, physiological interpretation, and clinical use,” Circulation, vol. 93, no. 5, pp. 1043–1065, 1996. View at: Publisher Site | Google Scholar
  15. W. Zareba and A. B. de Luna, “QT dynamics and variability,” Annals of Noninvasive Electrocardiology, vol. 10, no. 2, pp. 256–262, 2005. View at: Publisher Site | Google Scholar
  16. J. Pan and W. J. Tompkins, “A real-time QRS detection algorithm,” IEEE Transactions on Biomedical Engineering, vol. 32, no. 3, pp. 230–236, 1985. View at: Publisher Site | Google Scholar
  17. M. S. Woolfson, “Study of cardiac arrhythmia using the Kalman filter,” Medical & Biological Engineering & Computing, vol. 29, no. 4, pp. 398–405, 1991. View at: Publisher Site | Google Scholar
  18. F. Beckers, B. Verheyden, and A. E. Aubert, “Aging and nonlinear heart rate control in a healthy population,” American Journal of Physiology—Heart and Circulatory Physiology, vol. 290, no. 6, pp. H2560–H2570, 2006. View at: Publisher Site | Google Scholar
  19. P. K. Stein and A. Reddy, “Non-linear heart rate variability and risk stratification in cardiovascular disease,” Indian Pacing and Electrophysiology Journal, vol. 5, no. 3, pp. 210–220, 2005. View at: Google Scholar
  20. P. Guzik, J. Piskorski, T. Krauze et al., “Correlations between the Poincaré plot and conventional heart rate variability parameters assessed during paced breathing,” Journal of Physiological Sciences, vol. 57, no. 1, pp. 63–71, 2007. View at: Publisher Site | Google Scholar
  21. R.-G. Yeh, G.-Y. Chen, J.-S. Shieh, and C.-D. Kuo, “Parameter investigation of detrended fluctuation analysis for short-term human heart rate variability,” Journal of Medical and Biological Engineering, vol. 30, no. 5, pp. 277–282, 2010. View at: Publisher Site | Google Scholar
  22. R. Khan, “Identifying and understanding the role of pulmonary vein activity in atrial fibrillation,” Cardiovascular Research, vol. 64, no. 3, pp. 387–394, 2004. View at: Publisher Site | Google Scholar
  23. M. J. Shen and D. P. Zipes, “Role of the autonomic nervous system in modulating cardiac arrhythmias,” Circulation Research, vol. 114, no. 6, pp. 1004–1021, 2014. View at: Publisher Site | Google Scholar
  24. J. Chen, S. L. Wasmund, and M. H. Hamdan, “Back to the future: the role of the autonomic nervous system in atrial fibrillation,” Pacing and Clinical Electrophysiology, vol. 29, no. 4, pp. 413–421, 2006. View at: Publisher Site | Google Scholar
  25. M. Oliveira, M. N. da Silva, V. Geraldes et al., “Acute vagal modulation of electrophysiology of the atrial and pulmonary veins increases vulnerability to atrial fibrillation,” Experimental Physiology, vol. 96, no. 2, pp. 125–133, 2011. View at: Publisher Site | Google Scholar
  26. A. R. Magnano, S. Holleran, R. Ramakrishnan, J. A. Reiffel, and D. M. Bloomfield, “Autonomic nervous system influences on QT interval in normal subjects,” Journal of the American College of Cardiology, vol. 39, no. 11, pp. 1820–1826, 2002. View at: Publisher Site | Google Scholar
  27. R. E. Kleiger, P. K. Stein, and J. T. Bigger Jr., “Heart rate variability: measurement and clinical utility,” Annals of Noninvasive Electrocardiology, vol. 10, no. 1, pp. 88–101, 2005. View at: Publisher Site | Google Scholar

Copyright © 2016 S. Rangsungnoen et al. 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.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.