CC BY
76
A TIME-DOMAIN HYBRID ANALYSIS METHOD FOR DETECTING AND QUANTIFYING T-WAVE ALTERNANS
1 11 2 Xiangkui Wan , Kanghui Yan , Jun Zhang , Yanjun Zeng
T-wave alternans (TWA) in surface electrocardiograph (ECG) signals has been recognized as a marker of cardiac electrical instability, and is hypothesized to be related with patients at increased risk for ventricular arrhythmias. A novel timedomain TWA hybrid analysis method (HAM) utilizing the correlation method and least squares regression technique is described in this paper. Simulated ECGs containing artificial TWA (cases of absence of TWA and presence of stationary or time-varying or phase-reversal TWA) under different baseline wanderings are used to test the method, and the results shows the HAM has a better ability to quantifying TWA amplitude compared with the Correlation Method (CM) and Adapting Match Filter Method (AMFM). The HAM is subsequently used to analyze the clinical ECGs, and results produced by the HAM have, in general, demonstrated the consistency with those produced by the CM and the AMFM, while the quantifying TWA amplitudes by the HAM are universally higher than those by the other two methods.
Russ J Cardiol 2014, 4 (108), Engl.: 46-53
Key words: electrocardiograph, T-wave alternans, hybrid analysis method, correlation method, adapting match filter method, least squares regression.
1School of Information Engineering, Guangdong University of Technology, Guangzhou; 2Biomedical Engineering Center, Beijing University of Technology, Beijing, China.
Corresponding author. Yanjun Zeng, Professor of Biomechanics & Medical Information Institute, Beijing University of Technology, Beijing 100022, China, Tel: +8610 67391809, Fax: +8610 67391975, e-mail: yjzeng@bjut.edu.cn
AMFM — Adapting Match Filter Method, CM — Correlation Method, ECG — electrocardiograph, HAM — hybrid analysis method TWA — T-wave alternans.
Received September 11, 2013. Revision received September 18, 2013. Accepted September 25, 2013.
ГИБРИДНЫЙ МЕТОД АНАЛИЗА ВРЕМЕННОГО ДОМЕНА ДЛЯ ВЫЯВЛЕНИЯ И КОЛИЧЕСТВЕННОЙ ОЦЕНКИ ИЗМЕНЕНИЙ Т-ВОЛН
11 12 Xiangkui Wan , Kanghui Yan , Июнь Zhang , Yanjun Zeng
Изменения Т-волн (TWA) в сигналах электрокардиографии (ЭКГ) были признаны маркером электрической нестабильности сердца, и гипотетически связаны с пациентами, у которых повышен риск развития желудочковых аритмий. Новый гибридный метод анализа временного домена TWA (HAM) использующий метод корреляции и методики регрессии наименьших квадратов, описан в этой статье. Симулируемые ЭКГ, содержащие искусственные TWA (в случаях отсутствия TWA и присутствия стационарных или изменяющихся во времени или фаза-восстановление TWA) при различных блужданиях исходных показателей используются для тестирования метода, и результаты показывают, что HAM имеет лучшую возможность количественной оценки амплитуды TWA, по сравнению с корреляционным методом (СМ) и методом адапта-
ции матч-фильтрации (AMFM). HAM в дальнейшем используется для анализа клинических ЭКГ, и результаты, полученные НАМ, в целом, показали согласованность с теми, что были получены при проведенных СМ и AMFM, а количественная оценка амплитуд TWA при НАМ является повсеместно выше, чем при использовании двух других методов.
Российский кардиологический журнал 2014, 4 (108), Англ.: 46-53
Ключевые слова: электрокардиограф, изменения T-волн, гибридные методы анализа, метод корреляции, метод адаптации матч-фильтрации, метод регрессии наименьших квадратов.
Introduction
The T-wave alternans (TWA) has been considered as one of the most promising markers of sudden cardiac death (SCD) over past 10 years. TWA is a phenomenon appearing in the surface electrocardiograph (ECG) as a consistent fluctuation in the repolarisation morphology on an "every-other-beat" basis (2:1 behaviour). This fluctuation refers to a beat-to-beat variability in the amplitude, morphology and/or polarity of the T-wave. Numerous clinical studies have demonstrated TWA is associated with ventricular arrhythmias. Nowadays TWA has been considered an independent predictor for cardiac arrhythmias.
Several signal processing methods have been proposed to detect and estimate TWA in the ECG on a single-lead or multi-lead basis [1—8]. And a comprehensive and systematic discussion of methods for TWA detection and analysis is reported in [9]. Most widely used TWA detection methods work in two different domains: time and frequency.
The disadvantage of the frequency based methods is that they treat the alternans signal as a stationary wave with the constant amplitude and phase, which is not true in
general. They cannot detect non-stationary characteristics of the signal.
The time domain methods can detect TWA in short-time, non-stationary ECG signal, and they have also been used on Holter data. Correlation Method (CM) [6, 7] as a well known time-domain method performs well under different conditions, but it is sensitive to noise, especially the baseline wandering. In the presence of baseline oscillations at TWA frequency, a strong overestimation of TWA mean amplitude, and even TWA detection from TWA-free ECG tracings, is produced by the CM. And in the presence of higher frequency baseline fluctuations, the CM is not able to identify TWA [10]. An adapting match filter method (AMFM) was proposed by the same authors of the CM to overcome the CM limitations [11]. The AMFM yielded a significant improvement in algorithm-based identification of duration and amplitude of TWA from ECG tracings with frequency of baseline oscillations both lower and higher than that of TWA. Nevertheless, in the presence of baseline fluctuations at the TWA frequency, it produced erroneous TWA detection from ECG tracing
with no TWA, and even strong overestimation of TWA amplitude, when present.
Based on above background, we propose a hybrid approach for the TWA detection, which is based on correlation method and least squares regression technique. The study aims to develop a novel TWA detector to overcome the CM limitations, which can detect and measure transient TWA with more accuracy in the time domain, even in the presence of higher frequency baseline fluctuations.
The rest of the paper is organized as follows. In Section II, we present a novel method of TWA detection, simulated cases and clinical cases are also prepared. Then, in Section III, we report the results of its validation on the simulation database and clinical database, and compare the results to that of the CM and the AMFM. Next, in Section IV, we give the discussion. Finally, we summarize the conclusions of this work in Section V.
Material and methods
1. The hybrid analysis method (HAM) using correlation method and least squares regression technique. The hybrid analysis method consists of three different blocks: pre-processing, TWA detection and TWA evaluation. The whole TWA analysis process is described as below:
1) Data pre-process.
Before detecting TWA, the clinical ECG used here are required to be submitted to a preliminary pre-processing stage. This consists of various steps, which are: baseline wandering suppression, QRS complex detection and segmentation of the T-wave.
♦ Baseline wandering suppression: this is performed using a cubic spline interpolation technique [12].
♦ QRS complex detection: it is determined using a wavelet-based algorithm [13].
♦ T-wave segmentation: it is done by selecting intervals of 300 ms, beginning at a distance from the QRS fiducial point dependent on the RR interval. The interval onset for the ith beat, bt, is given by the expression:
b, = 4® +1.3 R V2 (ms)
♦ T-wave alignment: after T-wave segmentation, 128 consecutive T-waves present in the ECG are used to compute the median T-wave (Tm , which has each sample point given by the median value of the corresponding sample points of the 128 available T-waves), which is used as a template. Synchronization of the ith T-wave is performed according to a recursive procedure that keeps the segmented T-wave window length constant but varies its position ±30ms from the original position, with a time increment of one sample point. For each position of the T-wave window, the windowed ith T-wave is cross-correlated against the template. Optimal alignment occurs when maximum correlation is reached.
2) Qualitative detection of TWA.
After the data pre-processing, TWA is detected by looking for an alternating trend in the T-wave morphology quantified by a correlation index. To this aim, an alternans
correlation index (ACI) is computed to measure morphological changes of each of the consecutive T, waves in com-[12], which is as shown in equation (1)
XT, ( j)Tm ( j)
parison to T^
ACL = j
i = 1,2, —,128 (1),
N
X[Tm ( j]2
j =1
where Tm is the median T wave computed using 128 T waves available in each ECG tracing. N is the number of samples in each T wave.
ACI, is defined as the ratio of the maximum value of the cross-correlation function of Ti and Tm over the maximum value of the auto-correlation function of Tm . T, is classified as alternating.
The presence of TWA is considered when the value of ACI strictly oscillates (not necessary around one) in the case of monophasic TWA at least 7 consecutive beats. Figure 1 shows an example of alternating values of ACI,, indicating the presence of TWA.
To limit false detections caused by noise, a local threshold criterion, with ThACI equal to 0.06 [6], is considered, such that ACI values alternations have to exceed 0.12 for at least seven consecutive beats to be detected as TWA.
3) Quantitative estimation of TWA.
The odd and even beats of above detected consecutive beats are labelled as A and B , respectively. The odd T waves are obtained from A series and the even T waves are obtained from B series. The odd T waves constitute a matrix:
"(TA„
,TA„ ) =
TA A
Ta
ta A
ta
Where TA is the nth point of the mth odd T wave.
Anal-
gously the even T-wave
matrix TB can be constituted.
m,n
The amplitude corrections of odd and even T waves are performed using the first-degree polynomial as shown in equation (2).
f (i) = ® hk + b (2),
where Tt k is the ith row and kth column point of odd (or even) T-wave matrix. And the coefficients a, b are estimated by the linear least squares fitting process.
Each column vector TAk of T^is divided into 7-point epochs, and the equation (2) is recursively applied to each epoch throughout the entire TAk . Denote di as the ith deviation point of TAk from the fitting line:
- f (i)| (3).
^ = T
Then the mean deviation value of TA can be expressed
as equation (4):
(4).
A
0 , 0
A
1,0
1,1
1,n
T
A
mxn
T
T
A
A
A
m,0
m,1
o
m
If™ x(Oi) > 3 xO , then the Oi is considered to be corrected and replaced by the ith column mean value T(k) (as shown in equation (5)) of odd (or even) T-wave matrix.
T
T (k) =
i=0
(5).
And the amplitude correction of entire TAk is recalculated, until the m x(Oi) < 3 xO or n x(O1) ^ 2mV .
A specific example of amplitude correction of odd T wave matrix using linear fitting function is shown as Figure 2. The Figure 2 (a) represents the uncorrected T waves, and the Figure 2 (b) represents the corrected T waves.
Measure TWAk as the maximum absolute value of the
difference between TA and TB :
Ak Bk
TWA(k) = max' JTojrset
i T onset
TAl (i) - TBk (i) (6),
where TWA(k) denotes the kth local TWA (i.e. relative to a single odd (or even) beat), k = 1,2—,m.
The TWA of the analyzed consecutive ECG segment (segment TWA) is measured as the mean value of measured local TWAs:
m
X TWA(k)
TWA = -k=1-
1 rrjrLseg
(7).
And the global TWA (i.e., relative to the entire ECG tracing analyzed) is measured as the mean value of segment TWAs:
X TWAseg (i)
TWA = J
(8)
I
Above process can be described as the block diagram (Fig. 3).
2. Simulated cases. There is no generally accepted TWA-measuring criterion to be used as a gold-standard. Therefore, a simulation approach was used in the present study in different controlled cases.
A realistic, clean simulated ECG was obtained as a K-fold repetition of a single beat extracted from a real ECG [14].This guarantees that all the T waves of the simulated ECG are identical, so no TWA can be present in the original signal. In particular, we used a 0.7-s beat sampled at 500 samples per second. Length of each simulated ECG tracing was assumed to count 128 consecutive heart beats. Our choice relies on the fact that 128 consecutive beats were originally used for SM applications, and later on, this became the standard ECG length [5] for TWA detection and quantification. A constant RR interval of 0.7 s was assumed, so that TWA fundamental frequency was 0.71 Hz (that is, 1/ (0.7 x 2s) or 0.5 cycles per beat). TWA was simulated by varying T-wave amplitude (10, 50 and 100 mV) in a time window of 160 ms centered around the T-wave apex.
Four different sets of ECG simulation were considered, respectively reproducing the cases relative to the absence
of TWA, the presence of stationary TWA, the presence of time-varying TWA and the phase-reversal TWA, which are described below.
1) Case 1: simulated ECG tracing with no TWA.
The simulated ECG tracing with no TWA (N_TWA) is assumed not to be affected by any kind of noise. This simulated signal is thought to test the ability of recognizing the absence of TWA, which is represented in Figure 4 (a).
2) Case 2: simulated ECG tracings with stationary TWA.
The simulated ECG tracings with stationary TWA (S_TWA) are designed to test the ability of quantifying TWA amplitude in the presence of stationary alternating T-wave profiles. Three kinds of simulated ECG tracings were considered; namely, a tracing with a 10 mV TWA (S_TWA10), a tracing with 50 mVTWA (S_TWA50), and a tracing with 100 mVTWA (S_TWA100). An example of a tracing with 50 mV TWA is represented in Figure 4 (b).
3) Case 3: simulated ECG tracings with time-varying TWA.
ECG with visible TWA clearly shows the non-stationary nature of this phenomenon, whose variability often shows on-off or cyclic trends. Evaluation of dynamic aspects of TWA is important in clinics since transient TWA has been observed during acute ischemia [15]. To test the ability of the HAM in detecting non-stationary TWA, two simulated ECG tracings were considered, each one incorporating a specific beat-to-beat varying (and then, time-varying) A (n) sequence. A sinusoidal A (n) sequences, with 128 beats period, were affecting the first (TV_TWA1) ECG tracing, while An A (n) varying from 50 mV to 20 mV, following a smoothed (24 beats transition) step pattern, was affecting the second ECG tracing (cascaded TWA, TV_TWA2). The two simulated tracings were characterized by a uniform profile of TWA, which are represented in Figure 5 (a) and (b). The examples of TV_ TWA1 and TV_TWA2 are represented in Figure 6 (a) and (b), respectively.
4) Case 4: simulated ECG tracing with phase-reversal TWA.
Arrhythmias can sometimes trigger a phase reversal so that the alternans pattern changes from ABABAB to BABABA [5]. The simulated ECG tracing with phase-reversal TWA (PR_TWA) is designed to test the ability of the method in detecting phase-reversal TWA. PR_TWA tracing incorporates a stationary 10 mV TWA, which changes phase twice, at beats 40 and 80, respectively. This simulated case may also be used to help the interpretation of realistic cases in which a beat is missed (false negative QRS detection) or wrongly inserted (false positive QRS detection). An example is represented in Figure 6 (c).
Finally the noise is also considered to be added to above simulated ECG tracings in this study. In clinical settings, power line interfere is generally eliminated by hardware filter. When computing the ACI indexes (equation (1)) the white noise is already taken into account. So baseline wandering is
k
m
considered to the present simulated cases which might cause erroneous detection of TWA. Baseline wandering can be eliminated by pre-processing stage, whereas ECG amplitude modulation may survive, and elimination by pre-processing of ECG modulation related to T-wave variability should be prevented because TWA is a specific case of it [4]. Based on these considerations, ECG simulations with baseline wandering are considered. Baseline wanderings are simulated with a sinusoid of 0.1 mV amplitude and various frequencies: 0.30, 0.71 and 1.50 Hz, respectively, which we denote as bw030, bw071 and bw150. These frequencies are, respectively, lower, equal and greater than TWA frequency. The frequency of 0.30 Hz relates to usual breathing pattern in patients. And the baseline fluctuations are simply added to each simulated ECG tracing. Two representative examples of our simulated ECG tracings, with and without baseline fluctuations, are displayed in Figure 7.
3. Clinical cases. Two clinical data sets are considered in this study: ECG tracings from healthy subjects (^subjects) and that from patients.
ECG tracings from H_subjects belong to the Digital Electrocardiology Study databases of Liuhuaqiao Hospital, Guangzhou, which includes 320 Holter ECG tracings from H_subjects. The study was approved by the institutional research ethics committee of Guangzhou Medical College, and it was conducted following required rules for human subjects' research principles, according to the Declaration of Helsinki, as well as to Title 45, U. S. Code of Federal Regulations, Part 46, Protection of Human Subjects, Revised November 13, 2001, effective December 13, 2001. Each subject underwent 10-min ECG recording in resting conditions. Nine standard leads (V1-V6, I, II, and III) were recorded using equipment by Siemens-El-ema AB and digitized at a sampling rate of 500 Hz with amplitude resolution of 0.6 mV. Leads aVF, aVR, and aVL were derived from leads I. II, and III.
ECG tracings from patients belong to the T-Wave Alternans Challenge Database (TWACD) [16], which contains 100 multichannel ECG records sampled at 500 Hz with 16 bit resolution over a ± 32 mV range. The subjects include patients with myocardial infarctions, transient ischemia, ventricular tachyarrhythmia, and other risk factors for sudden cardiac death, as well as healthy controls and synthetic cases with calibrated amounts of T-wave alternans. the databases is chosen for two reasons: one is that previous studies found T-wave alternans episodes, some of them related to annotated ischemic episodes. Another is that the databases are well-known and available by many research groups.
In the specific, a group of fourteen healthy subjects was compared with a group of fourteen patients. A subject was classified as belonging to the H-group when fulfilling the following criteria [17]:
1). no overt cardiovascular disease or history of cardiovascular disorders (including stroke, TIA, and peripheral vascular disease);
Figure 1. An example of the presence of TWA.
Figure 2. A specific example of amplitude correction of odd T wave matrix. The corrected matrixes for odd and even T waves are known as TA and TB , respectively.
Figure 3. The block diagram of T-waves amplitude correction and TWA estimation.
Figure 4. An example of simulated ECG tracings with/without TWA. (a) simulated ECG tracings not affected by TWA (N_TWA); (b) simulated ECG tracings with stationary 50mV TWA (S_TWA50).
Figure 5. The different cases of nonstationary TWA. (a) sinusoidal trend of TWA amplitude signals; (b) cascaded TWA amplitude signals; (c) phase-reversal TWA.
main methods, which are the CM and the AMFM, are used here for comparison.
In our simulation study, the root mean square error (RMSE) in the estimate of TWA amplitudes are computed
[13]: i-
RMSEm Jx^JTWAM^A^ (9),
Figure 6. The simulated ECG tracings with different nonstationary TWA. (a) TV TWA1; (b) TV_TWA2; (c) PR_TWA.
Figure 7. Two examples of the simulated tracings affected by baseline wanderings. (a) The simulated tracings without baseline; (b) The simulated tracings with 0.3Hz baseline.
Figure 8. The time-varying TWA measurement results under different baseline conditions.
2). no history of high blood pressure (>150/90 mmHg);
3). not taking medication;
4). no other chronic illness (e.g. diabetes, asthma, chronic obstructive pulmonary disease, etc.);
5). diagnosed as being healthy if evaluated by a physician for cardiovascular-related syndrome (chest pain, palpitation, syncope);
6). normal physical examination;
7). sinus rhythm in 12-lead ECG without any suspicious abnormalities (e.g. signs ofventricular hypertrophy, inverted T-wave, intraventricular conduction disturbances);
8). normal echo and normal ECG exercise testing in presence of suspicious ECG changes;
9). no pregnancy.
4. Statistics. To evaluate the ability of the presented method to quantify TWA, the other two related time-do-
where N is the total number of beats in an ECG tracing, and A(n) (relative to the nth beat) is assumed equal to the absolute value of the maximum difference between the nth and the (n+1)th T-wave sample amplitude, and TWAM(n) is the estimated local TWA (relative to the nth beat) by the three competing methods. Subscript M is for either the HAM or the CM or the AMFM. In this study, the resolution of RMSE is 0.1mV and the predefined A(n) are considered as constitutive reference TWA-amplitude signals (gold-standard).
When analyzing clinical data, the Lilliefors test was used to evaluate the hypothesis that estimated TWA had a normal distribution (significance was set at 5% level) over a population. Comparisons between normal distributions were performed using Student's test, whereas distributions that could not be considered normal would be compared using the Wilcoxon rank sum test. Statistical significant differences were assumed for P<0.05.
Results
For the simulated data and clinical data set, ECG segments of 128 consecutive beats were randomly extracted and directly submitted to the AMFM, which does not require pre-processing [11]. Rather, a data preprocessing stage, described in 2.1 section, was performed prior to submit the CM and the HAM.
1. Simulated cases. For the simulated cases, results obtained from TWA analysis, by applying the CM, the AMFM and the HAM respectively, are reported in below tables.
In Table 1 are reported the results obtained from the simulated ECG tracing with no TWA (N_TWA). These three methods applying to tracings with no baseline yielded accurate identification of TWA amplitude. In the presence of 0.30 and 1.50 Hz baseline wandering, a slight overestimation of TWA amplitude was produced by the CM. In the presence of baseline fluctuations with frequency equal (0.71 Hz) to that of TWA, the strong overes-timations of TWA amplitude were produced by the three methods.
In Table 2 are reported the results obtained from the simulated ECG tracing with stationary TWA (S_TWA). In the presence of 0.30 and 1.50 Hz baseline wandering, the CM and the AMFM produced underestimation of TWA amplitude for different stationary TWA, while the HAM yielded an accurate identification of TWA (RMSE-IAM=0mV). In the presence of baseline fluctuations at the TWA frequency of 0.71 Hz, the CM produced underestimation of TWA amplitude, and the AMFM still produced
Tablel
TWA amplitude measurements and errors evaluation for N_TWA case
TWACM RMSECM RMSEAual AMFM RMSEmu HAM
N_TWA (no bw) 0 0 0 0 0 0
N_TWA (bw030) 2 0 0 4.6 0 0
N_TWA (bw071) 34 200 23 34.4 385.9 22.6
N_TWA (bw150) 1 0 0 5.2 0 0
Table 2
TWA amplitude measurements and errors evaluation for S_TWA case
™AoU twaamfm twaham rmsecm RMSEamfm RMSEham
S_TWA10 (no bw) 6.5 5.7 10.0 3.2 4.23 0
S TWA10 (bw030) 6.5 5.8 10.0 4.7 4.2 0
S_TWA10 (bw071) 12 219.7 10.1 21.6 209.7 0.2
S TWA10 (bw150) 7.2 5.9 10.0 l|26 4.2 0
S_TWA50 (no bw) 36 28.6 50.0 14.3 21.4 0
S TWA50 (bw030) 36 28.8 50.0 l|14.8 28.9 0
S_TWA50 (bw071) 17.3 198.8 32.0 32.4 198.8 18.4
S TWA50 (bw150) 38.2 29 50.0 11.7 29 0
S_TWA100 (no bw) 74.9 57.7 100.0 25.3 42.3 0
S TWA100 (bw030) 74.9 57.5 70.0 ||25.6 42.4 0
S_TWA100 (bw071) 54.7 176.2 100.0 45.0 76.0 29.5
S_TWA100 (bw150) 79.2 57.6 100.0 21.1 42.4 0
strong overestimation of TWA amplitude. While the HAM, even produced a slight underestimation of TWA amplitude, showed a better ability to quantifying TWA amplitude with more accuracy in this case, and RMSETT„,
J ' HAM
obtained 0mVfor S_TWA10, 18.4mVfor S_TWA50, 29.5mV for S_TWA100, respectively.
For TV_TWA1 and TV_TWA2 cases, the local TWA comparisons are considered because of time-varying amplitudes. A graphical representation of the results obtained from ECG simulations with presence of time-varying TWA (TV_TWA1 and TV_TWA2) is depicted in Figure 8. The columns of panels from left to right display simulated TWA-amplitude signals (128 beats), and detected TWA-amplitude signals provided by the CM, the AMFM and the HAM, respectively. For the cases of the simulated ECG tracing with 0.30 and 0.71 Hz baseline wandering, analogous results are obtained.
The root mean square errors obtained are reported in Table 3. The three methods were able to track the time course of TWA. But the local TWA-amplitude signals provided by the CM showed vigorous amplitude fluctuation, and RMSE are higher then RMSErtr and RMSE_
CM c AMFM HAM
uniformly. The CM and the AMFM produced underestimation of TWA amplitude, which are same to the above mentioned cases, while the HAM provided a good estimate of TWA (RMSEHM<1.5mV, except the case of frequency of baseline equal to that of TWA)
In Table 4 are reported the results obtained from the simulated ECG tracing with phase-reversal TWA
Table 3
TWA errors evaluation for TV_TWA1 and TV_TWA2 cases
RMSECU RMSEamfm Rmseham
TV_TWA1 (no bw) 13.8 12.7 0.9
TV TWA1 (bw030) 16.1 II 10.9 1.2
TV_TWA1 (bw071) 25.5 296.9 14.9
TV TWA1 (bw150) 16.1 II 12.6 1
TV_TWA2 (no bw) 6.9 2.8 0.5
TV TWA2 (bw030) 10.1 2.2 1.3
TV_TWA2 (bw071) 19.9 342.1 15.1
TV_TWA2 (bw150) 7.5 II 3.0 0.5
(PR_TWA). The AMFM produced underestimation of TWA amplitude (40%) in the presence of 0.30 and 1.50 Hz baseline wandering, while the CM and the HAM produced a good results (RMSECM=0mV, RMSEham = 0m V). In the presence of baseline fluctuations at the TWA frequency of 0.71 Hz, the three methods produced strong overestimation of TWA amplitude, but obviously the results provided by the HAM are more close to the simulated TWA (TWA„a =19.7u.V, RMSE-
HAM
HAM=13.W
2. Clinical cases. TWA levels quantified by the three competing methods in the H-subjects and patients data are reported in Table 5, The CM, the AMFM, and the HAM detected various levels of TWA in same H-subjects and all patients. TWA was detected in two H-subjects by the CM and the AMFM, while only one H-subjects was affected by TWA according to the HAM (Table 5). And the three
methods detected the presence of TWA in all patients. TWA showed a normal distribution over patients' populations. Mean TWA values estimated by the HAM in H-sub-jects (0.5±1.9mV) and patients (10.8±3.7^F) were higher than the corresponding mean TWA estimates provided by the AMFM (H-subjects: 0.4±1.9mV; patients: 8.6±3.3mV) and the CM (H-subjects: 0.5±1.3mF; patients: 9.5±3.5mF). All these methods provided mean TWA estimates which showed significant differences between H-subject and patient groups.
The CM, the AMFM and the HAM detected the presence of TWA in all patients, and provided similar TWA estimates. The CM and the AMFM tend to underestimate TWA (Fig. 9 and simulation study results), and this finding is confirmed by our clinical result.
Discussion
In this study four simulated cases were generated with characters of: absence of TWA; presence of different kinds of stationary TWA; presence of two kinds of non-stationary (time-varying) TWA; and presence of phase-reversal TWA. The other two time-domain methods, namely the CM and the AMFM, are compared with the HAM in detection TWA. Results of our simulation study indicate
that the HAM allows detection and quantification of TWA better than the CM and the AMFM.
The CM was found to underestimate TWA amplitude in the simulated ECG tracing, since it assumed TWA being distributed along the entire length of the T wave [18]. And in the case of ECG simulations with presence of time-varying TWA, the CM produced the worst results compared with other method (Figure 9 and Table. 3).
The AMFM showed a good performance of time-varying TWA detection, due to its heart-rate adaptive-match-filter yielded the suppression of all ECG and interferences frequency components, while it produced strong overesti-mation of TWA amplitude in the presence of baseline fluctuations at the TWA frequency, and the reason and a potential solution were given in literature [10].
We can find that the HAM yielded, in general, a more accuracy TWA estimation in the simulated cases, although in the presence of baseline fluctuations with frequency equal to that of TWA the deviation from TWA amplitude was produced which are also produced by the CM, and the reason is that the accuracy of isoelectric line estimation by the cubic spline interpolation technique reduces. And all simulation cases showed RMSETT„, were
HAM
systematically smaller than TWA„„ and TWA„r„ ,„ even in
J J CM AMFM"
TWA amplitude measurements and errors evaluation for PR_TWA cases
Table 4
TWA0M RMSE0M RMSEMlal AMFM rmseuau HAM
PR_TWA (no bw) 10 6 10 0 2.8 0
PR_TWA (bw030) 10 6 10 5.1 3.1 1.3
PR_TWA (bw071) 32 230 19.7 22.5 348.3 13.1
PR_TWA (bw150) 9 6 10 0 3.2 0
Table 5
TWA amplitude measurements of clinical data applying the CM, the AMFM and the HAM
H-subjects ™acu twaamfm TWAham TWACD twacu TWAamfm twahaam
1 0 0 0 Twa06 6.55 5.73 7.05
2 0 0 0 Twa09 || 6.91 1 6.21 1 8.11
3 0 0 0 Twa10 7.05 5.93 7.65
4 3.02 2.11 0 Twa18 3.75 3.23 4.25
5 0 0 0 Twa22 14.17 12.56 15.10
6 0 0 0 Twa23 12.91 | 10.83 | 13.47
7 0 0 0 Twa41 12.08 11.10 14.48
8 0 0 0 Twa46 5.95 4.91 7.55
9 4.12 4.08 7.11 Twa61 11.83 11.06 13.06
10 0 0 0 Twa71 7.12 7.10 | 8.80
11 0 0 0 Twa85 10.88 10.06 13.16
12 0 0 0 Twa92 13.10 12.55 14.21
13 0 0 0 Twa94 13.79 12.94 15.41
14 0 0 0 ||Twa99 H 7.16 | 6.38 | 8.23
0.5±1.3 0.4±1.9 0.5±1.9 9.5±3.5* 8.6±3.3* 10.8±3.7*
P < 0.05 when comparing H-subjects vs. patients with the t-test for normal distributions.
the presence of baseline fluctuations at the TWA frequency of 0.71 Hz.
The HAM performs a amplitude corrections procedure based on the linear least squares fitting technique before calculating the local TWA, which further suppresses the interferences, and the local threshold criterion, integrated in the HAM, appears to help improving detecting accuracy. The limitation of the CM is that then computing the ACI, the exact location of the maximum amplitude difference between the two waves is lost, so that a mean (over T wave) TWA amplitude value is provided (assumption of uniformly distributed TWA), while in our method TWA is measured by the maximum absolute value of the difference between the corrected matrixes for odd and even T waves, which also improves the accuracy of TWA estimation. The baselines with various frequencies are considered in the simulated cases, and the test results also show the HAM is robust to the noise.
Our results relative to the clinical data highlighted consistency in the detection and quantification of TWA by the three different methods, and significant differences between H-subject and patient groups are manifested, as shown in Table 5. While the TWA amplitudes measured by the CM and the AMFM are slight lower than that by the HAM. The results of our simulation test help interpretation of TWA data obtained from clinical cases.
Conclusions
A novel time-domain TWA detector is presented in this paper based on the correlation method and linear least squares fitting technique. Although the method is simple, it was validated using simulated ECG test signals with artificial TWA of various amplitudes and baseline wanderings, and achieved a good performance under reasonable levels
References
1. Adam DR, Akselrod S, Cohen RJ. Estimation of ventricular vulnerability to fibrillation through T-wave time series analysis, Comput. Cardiol. 1981; 8: 307.
2. Smith JM, Clancy EA, Valeri CR, et al. Electrical alternans and cardiac electrical instability, Circulation 1988;77:110.
3. Nearing BD, Huang AH, Verrier RL. Dynamic tracking of cardiac vulnerability by complex demodulation of the T wave, Science 1991; 252: 437.
4. Laguna P, Ruiz M, Moody GB. Repolarization alternans detection using the KL transform and the beatquency spectrum, Computing in Cardiology 1996; 23: 673.
5. Nearing BD, Verrier RL. Modified moving average analysis of T-wave alternans to predict ventricular fibrillation with high accuracy, J Appl Physiol 2002; 92: 541.
6. Burattini L, Zareba W, Moss AJ. Correlation method for detection of transient T-wave alternans in digital Holter ECG recordings, ANE 1999; 44:416.
7. Burattini L, Bini S, Burattini R. Comparative analysis of methods for automatic detection and quantification of microvolt T-wave alternans, Medical Engineering & Physics 2009; 31:1290.
8. Monasterio V, Laguna P, Martinez JP, Multilead analysis of T-wave alternans in the ECG using principal component analysis, IEEE Trans. Biomed. Eng. 2009; 56: 1880.
9. Mart'inez JP, Olmos S. Methodological principles of T wave alternans analysis: a unified framework, IEEE Trans Biomed Eng 2005;52: 599.
10. Burattini L, Zareba W, Burattini R. Automatic detection of microvolt T-wave alternans in Holter recordings: Effect of baseline wandering, Biomedical Signal Processing and Control 2006; 1: 162.
of noise. The results of our simulation study indicate that the HAM provides a more accurate TWA estimation than the CM and the AMFM.
Results of TWA detection produced by the three methods in real clinical ECG records show high consistency, which confirms the TWA detection power of the hybrid method for clinical data, although the quantifying TWA amplitudes by the HAM are universally higher than that by the CM and the AMFM.
Declarations
Funding: This work was supported by the National Nature Science Foundation of China (No. 60901027).
Competing Interests: All authors of the present work exclude any financial and personal relationships with other people or organizations that could inappropriately influence this job.
Ethical approval: The study was approved by the institutional research ethics committee of Liuhuaqiao Hospital, Guangzhou, and the Approved No. of ethic committee is No. 20111212. Informed consent was obtained from each subject.
The study protocol conforms to the ethical guidelines of the World Medical Association, Declaration of Helsinki — Ethical Principles for Medical Research Involving Human Subjects adopted by the 18 th WMA General Assembly, Helsinki, Finland, June1964, as revised in Tokyo 2004, as reflected in a priori approval by the appropriate institutional review committee. Declaration of Helsinki-Ethical Principles for Medical Research Involving Human Subjects adopted by the 18th WMA General Assembly, Helsinki, Finland, June1964, as revised in Tokyo 2004, as reflected in a priori approval by the appropriate institutional review committee.
11. Burattini L, Zareba W, Burattine R. Adaptive Match Filter Based Method for Time vs. Amplitude Characterization of Microvolt ECG T-Wave Alternans, Annals of Biomedical Engineering 2008; 36: 1558.
12. Meyer C, Keiser H. Electrocardiogram baseline noise estimation and removal using cubic splines and statespace computation techniques, Comput. Biamed. Res. 1977; 10:459.
13. Martinez JP, Almeida R, Olmos S, et al. A wavelet-based ECG delineator: Evaluation on standard databases, IEEE Trans. Biomed. Eng. 2004; 51: 570.
14. Narayan SM, Smith JM. Spectral analysis of periodic fluctuations in electrocardiographic repolarization, IEEE Trans Biomed Eng 1999; 46: 203.
15. Martinez JP, Olmos S, Wagner G, et al. Characterization of repolarization alternans during ischemia: time-course and spatial analysis, IEEE Trans. Biomed. Eng. 2006; 53:701.
16. Moody GB. The PhysioNet/Computers in Cardiology Challenge 2008: T-Wave Alternans. Computers in Cardiology, 2008; 35:505-508.
17. Burattini L, Bini S, Burattini R. Correlation method versus enhanced modified moving average method for automatic detection of T-wave alternans, Computer methods and programs in biomedicine 2010; 98: 94.
18. Burattini L, Bini S, Burattini R. Automatic microvolt T-wave alternans identification in relation to ECG interferences surviving preprocessing, Medical Engineering & Physics 2011; 33:17.
