Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 31 March 2021
Sec. Computational Physiology and Medicine
This article is part of the Research Topic Multi-Scale Computational Cardiology View all 14 articles

A New Method for Detecting Myocardial Ischemia Based on ECG T-Wave Area Curve (TWAC)

\r\nRonghua Li&#x;Ronghua Li1†Xiaoye Zhao&#x;Xiaoye Zhao2†Yinglan GongYinglan Gong1Jucheng ZhangJucheng Zhang3Ruiqing Dong*Ruiqing Dong4*Ling Xia*Ling Xia1*
  • 1Key Laboratory for Biomedical Engineering of Ministry of Education, Institute of Biomedical Engineering, Zhejiang University, Hangzhou, China
  • 2Department of Medical Imaging Technology, North Minzu University, Yinchuan, China
  • 3Department of Clinical Engineering, 2nd Affiliated Hospital, School of Medicine, Zhejiang University, Hangzhou, China
  • 4Dushu Lake Hospital Affiliated to Soochow University, Suzhou, China

In recent years, coronary heart disease (CHD) has become one of the main diseases that endanger human health, with a high mortality and disability rate. Myocardial ischemia (MI) is the main symptom in the development of CHD. Continuous and severe myocardial ischemia will lead to myocardial infarction. The clinical manifestations of MI are mainly the changes of ST-T segment of ECG, that is, ST segment and T wave. Nearly one third of patients with CHD, however, has no obvious ECG changes. In this paper, a new method for detecting MI based on the T-wave area curve (TWAC) was proposed. Through observation and analysis of clinical data, it was found that there exist significant correlation between the morphology of TWAC and MI. The TWAC morphology of normal subject is smooth and gentle, while the TWAC morphology of patients with coronary stenosis is mostly jagged, and the curve becomes more severe with more severe stenosis. The preliminary test results show that the sensitivity, specificity, and accuracy of the proposed method for detecting MI are 84.3, 83.6, and 84%, respectively. This study shows that the TWAC based approach may be an effective method for detecting MI, especially for the CHD patients with no obvious ECG changes.

Introduction

Cardiovascular disease has long been the first cause of death. In China only, the current number of patients with cardiovascular disease is 290 million, among them 11 million are patients with coronary heart disease (CHD) (Ma et al., 2020). For patients with typical CHD, diagnosis can be made based on resting electrocardiogram (ECG), exercise ECG, cardiac ultrasound, coronary angiography, and cardiac magnetic resonance. The resting ECG is widely used in clinic due to its non-invasive, economical and simple reasons. ST-T wave change is the most commonly used ECG feature to judge myocardial ischemia (MI) (Lin and Guo, 2005). According to related research, the positive rate of resting ECG diagnosis for CHD is 71%, and the rest of 29% patients have no ECG changes and some patients with three lesions can be completely normal, indicating that patients with normal ECG cannot rule out CHD (Pijls et al., 1995).

The clinical manifestations of MI are mainly the ST-T segment changes of the electrocardiogram, that is, the ST segment and T wave. The diagnosis of MI by ECG mainly depends on the characteristics of the ST-T segment. However, there are many factors affecting ST segment changes, such as axis shift, heart rate, electrode effects, body position changes, etc., leading to inaccurate detection of feature points, relying on ST segment shift to detect MI has a large false detection rate and missed detection rate. T wave represents the repolarization process of the ventricles. During MI, myocardial repolarization is delayed, resulting in changes in T wave morphology, such as biphasic T wave, inverted T wave, or high-tip T wave. In view of this characteristic of T wave, this paper proposes a method for diagnosing myocardial ischemia based on the T-wave area curve. By locating T-wave characteristic points, calculating T-wave area and drawing a curve with cardiac cycles, based on clinical data, the qualitative relationship between myocardial ischemia and T-wave area curve was analyzed.

Materials and Methods

T-Wave Onset and Offset Detection

T wave is an important part of ECG signal. Accurate localization and morphological recognition of T wave are basic indicators for diagnosis of MI, but the shape of T wave is variable. Hayden et al. (2002) show that when ischemic cardiomyopathy occurs, the T-wave shape will change accordingly, such as inverted, bi-phase, high-point, etc., and the low-frequency components near the end of the T-wave are more abundant than other bands, and more susceptible to noise and baseline drift. Accurate finding the location of the T-wave end point is with some difficulties. Therefore, T-wave detection needs to take into account the changes in T-wave morphology.

At present, T-wave ends can be detected by methods such as area method (Zhang et al., 2006; Vázquez-Seisdedos et al., 2011), wavelet transform (Martínez et al., 2004), pattern recognition (Saini et al., 2013), and artificial neural network (Maglaveras et al., 1998). Due to the noise of the T wave and the baseline drift, the wavelet transform method may contain aliasing part. Therefore, simply using the wavelet method has a higher false detection rate, and some waveforms caused by noise or baseline drift are falsely detected as T waves. The neural network method can adapt to the change of T-wave shape when it is used for feature detection. It has good robustness, but the algorithm is more complicated. For the area-dependent methods, the algorithms can adapt to abnormal changes in T-wave shape, and have strong waveform adaptabilities. Among them, Zhang’s algorithm (Zhang et al., 2006) is a good method for T-wave ends location. It was based on an indicator signal with mathematically proved consistency. It was robust to measurement noise, waveform shape changes and baseline drift, and is suitable for various forms of T waves. The computation burden of the algorithm was very low: its main computation can be implemented as a simple FIR filter. When evaluated with the PhysioNet QT database (Goldberger et al., 2000) in terms of the mean and the standard deviation of the T-wave end location errors, Zhang’s algorithm outperforms the other algorithms evaluated with the same database (Zhang et al., 2006). So, in this study, Zhang’s algorithm is basically used for T-wave ends location.

One key issue of the area-dependent methods is to accurately determine the search boundaries, but the search boundaries are closely related to the RR interval. If the interval of the searching window’s boundaries was set too small which means that two boundary points are near the current R peak, the maximum of sliding area could not be found or the detected onset/offset of T wave are nearer to the R peak. This issue will affect detection accuracy, which results in detection error and vice versa. In order to more accurately model the relationships between RR interval and the searching boundaries, in this study, similar to Shang et al.’s work (Shang et al., 2019), we performed a k-means clustering analysis between RR intervals and RTon (RTon denotes the time interval between the R peak and T wave onset) as well as the relationship between the RR intervals and RToff (RToff the time interval between the R peak and T wave offset), which is implemented by means of the k-means function. The scatter plots with the optimal k-means clustering (k = 3) are showed in Figure 1, and k is determined by combining the results of clustering and the computational complexity of parameters’ settings as well as the adaptiveness of the algorithm. Then, the two relationships (between RR intervals and RTon, and between RR intervals and RToff) are obtained using the following equations:

FIGURE 1
www.frontiersin.org

Figure 1. Clustering results for T wave feature points: (A) clustering information of T wave onsets; (B) clustering information of T wave offsets.

C a s e 1 : R R < 0.67 s , 0.05 s < R T o n < 0.25 s C a s e 2 : 0.67 s R R < 1.23 s , 0.05 s < R T o n < 0.35 s C a s e 3 : R R 1.23 s , 0.05 s < R T o n < 0.45 s C a s e 1 : R R < 0.71 s , 0.2 s < R T o f f < 0.45 s C a s e 2 : 0.71 s R R < 1.1 s , 0.2 s < R T o f f < 0.6 s C a s e 3 : R R 1.1 s , 0.2 s < R T o f f < 0.45 s (1)

Thus, the three piecewise functions for determining the search boundaries for T wave onset and offset detections can be obtained.

Then, the grid search was used to determine the best combination of parameters in Equations (2) and (3), which was implemented by for loop. In a loop, we changed the value of one parameter at a time, kept the other parameters unchanged, and applied the algorithm in the QT database as well as using a fivefold cross-validation. Then, we stored the sensitivity of one loop and started another loop. Through all loops, we traversed all of the combinations of parameters. After comparing the results, the combinations of parameters with the highest sensitivity were chosen. The best parameters’ combinations for T wave onsets are: a = 0.4, b = 0.2, c = 0.4, d = 0.4, e = 0.3, f = 0.0, and for T wave offsets are: a = 0.2, b = 0.1, c = 0.2, d = 0.0, e = 0.0, f = 0.1.

{ t 1 = ( a × R R i + R i + 0.02 ) s t 2 = ( b × ( R R i ) + R i + 0.16 ) s if R R i < 0.67 s t 1 = ( c × R R i + R i + 0.04 ) s t 2 = ( d × ( R R i ) + R i + 0.24 ) s i f 0.67 s R R i < 1.23 s t 1 = ( e × R R i + R i + 0.04 ) s t 2 = ( f × ( R R i ) + R i + 0.4 ) s if R R i 1.23 s (2)
{ t 1 = ( a × R R i + R i + 0.18 ) s t 2 = ( b × ( R R i ) + R i + 0.3 ) s if R R i < 0.71 s t 1 = ( c × R R i + R i + 0.18 ) s t 2 = ( d × ( R R i ) + R i + 0.4 ) s if 0.71 s R R i < 1.1 s t 1 = ( e × R R i + R i + 0.18 ) s t 2 = ( f × ( R R i ) + R i + 0.48 ) s if R R i 1.1 s (3)

T-Wave Area Curve

When the T-wave onset and offset was detected, the T-wave area can be calculated as follows.

Let the T wave onset be Ton and the T wave offset be Toff. Within the fixed window t ∈ [Ton, Toff], calculate the waveform area A_t:

A t = t = T o n T o f f ( s t - s k ¯ ) (4)

where s_t is the amplitude of the t-th sample point, and sk¯ is the local average amplitude (by default, a smoothing window of p = 0.016 s is used), which is defined as:

s ¯ k = 1 2 p + 1 j = t - p t p S j (5)

Calculate the continuous T wave area of each lead, and then draw the connection line with the cardiac cycle number as the abscissa and the T wave area as the ordinate. Based on the T wave morphological variability during MI, it can be inferred that if the line is approximately straight, the ischemia test is negative. If one or more of the leads are serrated, it is positive, and if one or more of the leads are serrated, it may be related to the position of the blocked coronary artery is related.

From the T-wave area curves of 52 healthy samples in the PTB database, it was found that the curves of 43 healthy people were flat and the morphological differences were small. The difference in T-wave area per heartbeat of patients with MI, however, is obviously larger than that of healthy people. The area curve is jagged and irregular. Observing the corresponding T-wave area curve of the 15-s electrocardiogram data of 148 patients with MI in the PTB database, it was found that the curves of 115 MI patients were irregular and chaotic, which was in sharp contrast with the curves of healthy samples.

According to the characteristic that the degree of fluctuation of the T-wave area curve of patients with MI is significantly greater than that of healthy people, we used the TWAC (T-wave area curve) to detect MI. Based on the TWAC form, the gentle TWAC is defined as negative, corresponding to myocardial blood supply was normal in the subject. And the irregular jagged TWAC was defined as positive, corresponding to myocardial ischemia in the subject.

Experimental Data

The sample number in this article is 364 in two groups: the first group contains 148 patients with MI and 52 normal persons from the PTB database (Goldberger et al., 2000); the second group contains 122 patients with suspected MI from Zhejiang Second Hospital and 42 health subjects. The detailed characteristics of the selection patients are shown in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Clinical characteristics of patients.

The clinical data of Zhejiang Second Hospital was conducted from May 2018 to December 2018. Philips TC20 ECG machine was used to record resting ECG data with a sampling rate of 1,000 Hz and 16-bit resolution. We select stable 15 s ECG data (15 s ECG data contains about 20 heart beats) to draw a 12-lead T-wave area curve. This study was approved by Zhejiang Second Hospital Review Board. After signing the informed consent, the patient was placed in the supine position, and the resting ECG was obtained after strictly following the 12-lead ECG collection procedure. Each patient was subjected to coronary angiography after ECG examination. A professional cardiology interventional doctor performed a visual assessment of the degree of coronary artery stenosis. In order to ensure the accuracy of the ECG results of this study, a professional electrocardiologist interpreted the ECG data, but he did not know the coronary angiography examination results of the patients.

Selection criteria were patients with suspected MI who had no MI characteristics on the ECG. The absence of MI characteristics means that the ST segment of any lead with R wave as the main wave is not depressed or depressed < 0.05 mv, T wave is upright and ≥ 1/10R wave. Suspicious MI refers to the clinical manifestations, myocardial enzymes and other tests and coronary CTA and other examinations suggest that the patient may have MI, which is evaluated by a professional cardiologist.

Exclusion criteria were non-sinus electrocardiograms such as atrial fibrillation and atrial flutter; premature beats; atrioventricular block or ventricular block; ventricular pre-excitation patterns; abnormal Q waves or poor R-wave increments in right chest leads; significant sinus Bradycardia (heart rate < 50 beats/min); heart valve disease.

This article used the results of CAG or blood flow reserve fraction (FFR) as the gold standard for MI diagnosis, and defined coronary angiography result as a positive MI if one of the three branches of the coronary artery (anterior descending branch, circumflex branch, right coronary artery) has a stenosis degree of ≥ 70% (Tonino et al., 2009); all without stenosis or a degree of stenosis less than 70% are negative MI. Patients are considered positive MI when FFR value ≥ 0.8, and negative MI if FFR value < 0.8.

Results

As mentioned above, TWAC is a curve of the T-wave area of a conventional 12-lead ECG as a function of the cardiac cycle. By analyzing the TWAC morphology of a lot of clinical data, it is found that the curve of normal people has less fluctuation, and the curve of patients with MI is mostly jagged changes in different cardiac cycles. Figure 2 shows an example of TWAC in two normal people, and Figure 3 shows an example of TWAC in two patients with MI. It shows that the T-wave area curve has a significant correlation with ischemic heart disease. Therefore, in this article, TWAC of one or more leads with jagged fluctuations is identified as positive (with MI), and the fluctuations are small and gentle in TWAC is considered as negative (no MI).

FIGURE 2
www.frontiersin.org

Figure 2. ECG and TWAC Examples of two normal subjects: A005, male, 29 years old; A222, male, 23 years old. (A) ECG of A005; (B) TWAC of A005; (C) ECG of A222; (D) TWAC of A222.

FIGURE 3
www.frontiersin.org

Figure 3. ECG and TWAC examples of two MI patients: T209, female, 66 years of age, with moderate coronary stenosis; T234, male, 73 years, with severe coronary stenosis. (A) ECG of T209; (B) TWAC of T20; (C) ECG of T234; (D) TWAC of T234.

Experiment results are shown in Table 2. The sensitivity of TWAC to the diagnosis of MI was 84.3% (199/236), the specificity was 83.6% (107/128), and the accuracy was 84.0% (306/364).

TABLE 2
www.frontiersin.org

Table 2. Experimental results, where TN, TP, FN, and FP represent true negative, true positive, false negative, and false positive, respectively.

Discussion

Analysis of Misjudgment Data

The total number of patients in the false-negative group was 37, of which 13 patients had a diseased vessel stenosis of 70%. The method used in this study to interpret the degree of coronary stenosis was physician visual assessment (PVA), which is based on the personal experience of the surgeon. Judgment is relatively subjective. Studies have pointed out that the severity of stenosis of coronary lesions in PVA in China is significantly higher than quantitative coronary angiography (QCA) and there is a large difference between hospitals and doctors (Zhang et al., 2018). Therefore, the presence of false-negative patients does not rule out the possibility of overestimation of coronary artery disease caused by PVA. Single-vessel disease is more common in false-negative patients, while multiple-vessel disease is more common in true-positive patients. The range of MI caused by single vessel disease is relatively small, and negative results are easily obtained. In addition, collateral circulation is another factor of false negative results. The collateral circulation rate of patients in the false negative group is higher than that in the true positive group. Professional cardiologists have confirmed that there are 4 cases of coronary stenosis >70% in the false negative group. Collateral circulation can improve the heart blood supply and may make TWAC negative.

A total of 21 false-positive patients was included in this study, 15 of whom had chest tightness and chest pain. Although coronary angiography showed negative results, the possibility of MI could not be completely ruled out. In the FAME study, 35% of the lesions with a degree of coronary stenosis of 50–70% had an FFR of < 0.8, and in Park’s study, 16% of the lesions with a degree of coronary stenosis of < 50% had an FFR of < 0.8. The data indicate that MI is still possible without significant coronary stenosis (Tonino et al., 2009; Park et al., 2012; Ahmadi et al., 2015). The clinical manifestations of all these false positive patients were typical CHD and showed a higher incidence of coronary atherosclerosis than the true negative group.

Analysis of the Relationship Between TWAC Morphology and Coronary Occlusion in Specific Leads

When the TWAC is positive, by observing the TWAC morphology on each lead, we’ve found that in addition to the V1,V3,V4 leads have obvious sharp points or inverted phenomena, the aVL lead and aVR lead curves also fluctuate greatly, as shown in Figure 3D.

T wave can reflect the heterogeneity of ventricular repolarization in patients and predict cardiovascular disease to a certain extent. T wave is a potential wave formed by repolarization of ventricular cells. Repolarization is an active energy-consuming process. When MI occurs, the heart cannot normally deliver blood and nutrients. Therefore, insufficient supply of myocardial energy will cause myocardial contraction and diastolic function is impaired, which may cause T wave changes in patients with ischemic cardiomyopathy.

Studies have shown that the degree of T wave changes in aVL leads reflects the degree of ventricular muscle excitement recovery time, and reflects the heterogeneity of ventricular muscle repolarization. Compared with ST-T changes, aVL lead T-wave changes have a higher sensitivity for the diagnosis of myocardial ischemia, which is significantly related to the number of coronary artery disease vessels and the degree of myocardial ischemic injury (Tepetam et al., 2016).

Studies have also shown that the inversion of T waves on leads V3 and V_4 may indicate a middle obstruction of the left anterior descending coronary artery, which indicates that changes in T waves on specific leads can predict not only extensive myocardial ischemia, but also blockage at specific coronary arteries (De Zwaan et al., 1982). From the azimuth point, leads V3 and V_4 correspond to the proximal anterior descending branch, and the coronary direction of the middle part of the left anterior descending branch is almost parallel to the aVL lead of the ECG. It can be speculated that if the adjacent area of the left anterior descending branch ischemia, the T wave of lead aVL should also change. Therefore, the T wave change of lead aVL also corresponds to the anterior descending branch of the left anterior descending coronary artery.

As shown in Figure 3B, the aVR lead curve fluctuates in a zigzag manner, and the coordinate points in some cardiac cycles exceed the abscissa and become positive values. In a normal 12-lead ECG, the T wave on the aVR lead is inverted. This patient’s CAG showed 70% stenosis of the left circumflex branch. Studies have shown that the morphological changes of the T wave in lead aVR are of great significance in predicting cardiovascular death, and its value is higher than other ECG leads, comparable to the changes in the ST segment of the aVR lead (Tan et al., 2008). If the amplitude of the inverted T wave becomes smaller, it means that the risk of cardiovascular death is gradually increased. When the inverted T wave becomes upright, the risk of cardiovascular death is higher.

The aVR lead has been used only to judge the origin of the heart rhythm, and its role has been seriously underestimated. The aVR lead has a special position on the frontal six-axis system, that is, the aVR lead axis is between the I and II lead axes, the angle between the aVR lead axis and the ventricular depolarization vector is the smallest, and the projection is the largest, which is the most sensitive lead to the changes of the ventricular depolarization vector.

Advantages and Limitations

Coronary atherosclerosis causes a series of electrophysiological changes that affect ventricular repolarization (Downar et al., 1977; Janse and Wit, 1989). During cardiac ischemia, the duration of action potential and conduction velocity decrease, leading to a heterogeneous repolarization process (Janse et al., 1985). Studies have shown that ischemia increases the repolarization dispersion between normal and ischemic fibers, and between the epicardium and the endocardium (Coronel et al., 1988), which refers to the “every other heartbeat” on the ECG. The repolarization pattern has a continuous fluctuation (Arini et al., 2014). This fluctuation refers to the change in the amplitude of the T wave or the change in the ST segment between different cardiac cycles. The amplitude of these bipolar alternations (dispersions) is usually in the microvolt range and cannot be visually recognized. Computer-based signal processing and analysis technology can detect subtle ECG changes. The T wave amplitude and shape on ECG alternately change from beat to beat, which is called T wave electrical alternation (TWA) (Puletti et al., 1980). TWA represents the alternation of cardiac repolarization, is an indicator of ventricular tachycardia and ventricular fibrillation in ischemic myocardium, and can be used as an indicator of risk stratification of acute myocardial infarction. The disadvantage of TWA is that it is susceptible to breathing, electrode and skin interference, wire movement and body position changes, and further research is needed. TWAC analyzes subtle ECG signal changes and amplifies such subtle changes to detect abnormal dynamics of cardiac repolarization, and is robust to acquisition noise, baseline drift and T wave morphology.

TWAC is the T-wave area curve of 12 leads with the cardiac cycle. Different leads of the electrocardiogram record electrical signals at different positions of the heart. Therefore, the degree of fluctuation of TWAC on different leads reflects the degree of stenosis at specific blood vessels to a certain extent. This paper analyzed the TWAC and coronary angiography results of some patients, and found that the T wave changes on leads V3,V4, aVL corresponded to the obstruction of the left anterior descending branch of the coronary artery, and the T wave changes on the aVR lead predicted the stenosis at left trunk coronary arteries.

There are some limitations in this study. First, coronary angiography shows that coronary artery stenosis is not equivalent to MI (Tonino et al., 2009; Park et al., 2012; Ahmadi et al., 2015). Clinically, it can be directly intervened when the diameter of coronary artery is narrower than 90%; if the diameter of coronary artery is not narrower than 90%, it is recommended that only FFR ≤ 0.8, or the disease with corresponding evidence of ischemia, then intervention can be taken. For moderate coronary stenosis, even the experienced cardiologist’s interventional physician’s visual assessment of angiography cannot accurately evaluate its physiological significance (Fischer et al., 2002). The number of patients with mismatched TWAC and CAG results in this study was 58, of which 33 (56.7%) patients had coronary stenosis with a coronary stenosis of 40–70%. It cannot be ruled out that coronary angiography could not determine the myocardial ischemia relatively accurately. The accuracy rate of FFR value ≤ 0.80 for identifying myocardial ischemia caused by coronary artery stenosis is as high as >90%. It has been confirmed by extensive randomized controlled studies that FFR is the gold standard for evaluating the physiological significance of coronary artery stenosis. Changes in hemodynamic factors such as contractile force. In future studies, FFR will be considered as the gold standard for diagnosis of MI, which can further verify the accuracy of TWAC. Secondly, the sample size is small, and more clinical research is needed to further verify the TWAC method.

Conclusion

In this paper, a new method for detecting MI was proposed, unlike the conventional ST-T segment based approaches, it was based on the T-wave area curve (TWAC). Preliminary test results show that the proposed method has good sensitivity, specificity, and accuracy for MI detection, especially for the CHD patients with no obvious ECG changes.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Zhejiang Second Hospital Review Board. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

RL and XZ contributed to the algorithm and the statistical analysis. JZ and RD contributed to the clinical data collection and interpretation. YG, RD, and LX contributed to the funding acquisition, conception, and design of the study. All authors contributed to the writing, critical reading, and approval of the manuscript.

Funding

This work was supported by the Natural Science Foundation of China (NSFC) under grant nos. 61527811 and 61701435, the Key Research and Development Program of Zhejiang Province under grant no. 2020C03016, the Zhejiang Provincial Natural Science Foundation of China under grant no. LY17H180003, and the Medical Health Science and Technology Project of Zhejiang Provincial Health Commission under grant no. 2020RC094.

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.

References

Ahmadi, A., Kini, A., and Narula, J. (2015). Discordance between ischemia and stenosis, or PINSS and NIPSS: are we ready for new vocabulary? JACC Cardiovasc. Imaging 8, 111–114. doi: 10.1016/j.jcmg.2014.11.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Arini, P. D., Baglivo, F. H., Martínez, J. P., and Laguna, P. (2014). Evaluation of ventricular repolarization dispersion during acute myocardial ischemia: spatial and temporal ECG indices. Med. Biol. Eng. Comput. 52, 375–391. doi: 10.1007/s11517-014-1136-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Coronel, R., Fiolet, J. W., Wilms-Schopman, F. J., Schaapherder, A. F., Johnson, T. A., Gettes, L. S., et al. (1988). Distribution of extracellular potassium and its relation to electrophysiologic changes during acute myocardial ischemia in the isolated perfused porcine heart. Circulation 77, 1125–1138.

Google Scholar

De Zwaan, C., Bär, F. W. H. M., and Wellens, H. J. J. (1982). Characteristic electrocardiographic pattern indicating a critical stenosis high in left anterior descending coronary artery in patients admitted because of impending myocardial infarction. Am. Heart J. 103, 730–736.

Google Scholar

Downar, E., Janse, M. J., and Durrer, D. (1977). The effect of acute coronary artery occlusion on subepicardial transmembrane potentials in the intact porcine heart. Circulation 56, 217–224. doi: 10.1161/01.CIR.56.2.217

CrossRef Full Text | Google Scholar

Fischer, J. J., Samady, H., McPherson, J. A., Sarembock, I. J., Powers, E. R., Gimple, L. W., et al. (2002). Comparison between visual assessment and quantitative angiography versus fractional flow reserve for native coronary narrowings of moderate severity. Am. J. Cardiol. 90, 210–215. doi: 10.1016/S0002-9149(02)02456-6

CrossRef Full Text | Google Scholar

Goldberger, A. L., Amaral, L. A. N., Glass, L., Hausdorff, J. M., Ivanov, P. C., Mark, R. G., et al. (2000). PhysioBank. PhysioToolkit, and PhysioNet: components of a new research resource for complex physiologic signals. Circulation 101, e215–e220. doi: 10.1161/01.CIR.101.23.e215

CrossRef Full Text | Google Scholar

Hayden, G. E., Brady, W. J., Perron, A. D., Somers, M. P., and Mattu, A. (2002). Electrocardiographic T-wave inversion: differential diagnosis in the chest pain patient. Am. J. Emerg. Med. 20, 252–262. doi: 10.1053/ajem.2002.32629

PubMed Abstract | CrossRef Full Text | Google Scholar

Janse, M. J., Capucci, A., Coronel, R., and Fabius, M. A. (1985). Variability of recovery of excitability in the normal canine and the ischaemic porcine heart. Eur. Heart J. 6, 41–52. doi: 10.1093/eurheartj/6.suppl_D.41

CrossRef Full Text | Google Scholar

Janse, M. J., and Wit, A. L. (1989). Electrophysiological mechanisms of ventricular arrhythmias resulting from myocardial ischemia and infarction. Physiol. Rev. 69, 1049–1169. doi: 10.1152/physrev.1989.69.4.1049

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, R., and Guo, J. (2005). Modern Perspective of Cardiovascular Disease. Beijing: Peking University Medical Press.

Google Scholar

Ma, L., Chen, W., Gao, R., Liu, L., Zhu, M., Wang, Y., et al. (2020). China cardiovascular diseases report 2018: an updated summary. J. Geriatr. Cardiol. 17, 1–8. doi: 10.11909/j.issn.1671-5411.2020.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Maglaveras, N., Stamkopoulos, T., Pappas, C., and Strintzis, M. G. (1998). An adaptive backpropagation neural network for real-time ischemia episodes detection: development and performance analysis using the European ST-T database. IEEE Trans. Biomed. Eng. 45, 805–813. doi: 10.1109/10.686788

CrossRef Full Text | Google Scholar

Martínez, J. P., Almeida, R., Olmos, S., Rocha, A. P., and Laguna, P. (2004). A wavelet-based ECG delineator: evaluation on standard databases. IEEE Trans. Biomed. Eng. 51, 570–581. doi: 10.1109/TBME.2003.821031

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, S. J., Kang, S. J., Ahn, J. M., Shim, E. B., Kim, Y. T., Yun, S. C., et al. (2012). Visual-functional mismatch between coronary angiography and fractional flow reserve. JACC: Cardiovasc. Interv. 5, 1029–1036. doi: 10.1016/j.jcin.2012.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Pijls, N. H., Van Gelder, B., Van der Voort, P., Peels, K., Bracke, F. A., Bonnier, H. J., et al. (1995). Fractional flow reserve.a useful index to evaluate the influence of an epicardial coronarystenosis on myocardial blood flow. Circulation 92, 3183–3193. doi: 10.1161/01.CIR.92.11.3183

CrossRef Full Text | Google Scholar

Puletti, M., Curione, M., Righetti, G., and Jacobellis, G. (1980). Alternans of the ST segment and T wave in acute myocardial infarction. J. Electrocardiol. 13, 297–300. doi: 10.1016/s0022-0736(80)80035-5

CrossRef Full Text | Google Scholar

Saini, I., Singh, D., and Khosla, A. (2013). Delineation of ECG wave components using K-Nearest neighbor (KNN) algorithm: ECG wave delineation using KNN. In proceedings of the 10th International Conference on Information Technology: New Generations. Piscataway NJ: IEEE 2013, 712–717.

Google Scholar

Shang, H., Wei, S., Liu, F., Wei, D., Chen, L., and Liu, C. (2019). An improved sliding window area method for t wave detection. Comput. Math. Methods Med. 11:3130527. doi: 10.1155/2019/3130527

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, S. Y., Engel, G., Myers, J., Sandri, M., and Froelicher, V. F. (2008). The prognostic value of T wave amplitude in lead aVR in males. Ann. Noninvasive Electrocardiol. 13, 113–119. doi: 10.1111/j.1542-474X.2008.00210.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tepetam, F. M., Dağdeviren, B., Bulut, I., Karabay, C. Y., Barış, S., and Aydıner Karakoç, E. (2016). A patient with mushroom allergy; a new etiological agent of Kounis syndrome. Tuberk. Toraks 64, 171–174. doi: 10.5578/tt.9411

PubMed Abstract | CrossRef Full Text | Google Scholar

Tonino, P. A. L., De Bruyne, B., Pijls, N. H. J., Siebert, U., Ikeno, F., van’ t Veer, M., et al. (2009). Fractional flow reserve versus angiography for guiding percutaneous coronary intervention. N. Engl. J. Med. 360, 213–224. doi: 10.1056/NEJMoa0807611

PubMed Abstract | CrossRef Full Text | Google Scholar

Vázquez-Seisdedos, C. R., Neto, J. E., Marañón Reyes, E. J., Macias, F. S., Rodríguez, V. R., Chacón, L. M., et al. (2011). New approach for T-wave end detection on electrocardiogram: Performance in noisy conditions. Biomed. Eng. Online 10:77. doi: 10.1186/1475-925X-10-77

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Mu, L., Hu, S., Nallamothu, B., Lansky, A. J., Xu, B., et al. (2018). Comparison of physician visual assessment with quantitative coronary angiography in assessment of stenosis severity in China. JAMA Intern. Med. 178, 239–247. doi: 10.1001/jamainternmed.2017.7821

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q., Manriquez, A. I., Médigue, C., Papelier, Y., and Sorine, M. (2006). An algorithm for robust and efficient location of T-wave ends in electrocardiograms. IEEE Trans. Biomed. Eng. 53(12 Pt 1), 2544–2552. doi: 10.1109/TBME.2006.884644

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: myocardial ischemia, ECG, coronary heart disease, heart, electrophysiology

Citation: Li R, Zhao X, Gong Y, Zhang J, Dong R and Xia L (2021) A New Method for Detecting Myocardial Ischemia Based on ECG T-Wave Area Curve (TWAC). Front. Physiol. 12:660232. doi: 10.3389/fphys.2021.660232

Received: 29 January 2021; Accepted: 11 March 2021;
Published: 31 March 2021.

Edited by:

Kuanquan Wang, Harbin Institute of Technology, China

Reviewed by:

Songyun Wang, Renmin Hospital of Wuhan University, China
Xin Li, University of Leicester, United Kingdom

Copyright © 2021 Li, Zhao, Gong, Zhang, Dong and Xia. 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: Ruiqing Dong, ruiqingdong@163.com; Ling Xia, xialing@zju.edu.cn

These authors share first authorship

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.