- Research
- Open Access

# Retrieving hidden atrial repolarization waves from standard surface ECGs

- Wei-Hua Tang
^{1}, - Wen-Hsien Ho
^{2, 3}Email author and - Yenming J. Chen
^{4}Email author

**17 (Suppl 2)**:146

https://doi.org/10.1186/s12938-018-0576-3

© The Author(s) 2018

**Published:**6 November 2018

## Abstract

### Background

This study estimates atrial repolarization activities (*Ta* waves), which are typically hidden most of the time from body surface electrocardiography when diagnosing cardiovascular diseases. The morphology of *Ta* waves has been proven to be an important marker for the early sign of inferior injury, such as acute atrial infarction, or arrhythmia, such as atrial fibrillation. However, *Ta* waves are usually unseen except during conduction system malfunction, such as long QT interval or atrioventricular block. Therefore, justifying heart diseases based on atrial repolarization becomes impossible in sinus rhythm.

### Methods

We obtain TMPs in the atrial part of the myocardium which reflects the correct excitation sequence starting from the atrium to the end of the apex.

### Results

The resulting TMP shows the hidden atrial part of ECG waves.

### Conclusions

This extraction makes many diseases, such as acute atrial infarction or arrhythmia, become easily diagnosed.

## Keywords

- Atrial repolarization waves
- Reaction–diffusion system
- Electrophysiology model
- Kalman filter
- Monte Carlo Markov chain

## Background

Atrial repolarization waves (*Ta* waves) are equally important as ventricular repolarization waves (*T* waves) and can exhibit significant potential as an effective biomarker for clinic diagnosis [1, 2]. However, depolarization is stronger than repolarization, and ventricular activities are even considerably stronger than atrial activities. To separate relatively weak activities in atria, existing studies have invested significant efforts in this aspect to help in myocardial diagnosis for either sinus rhythms and arrhythmia [3–8]. In the literature, no study has yet retrieved atrial repolarization activities because these activities can only be normally observed in patients under certain special ectopic pacing [2].

In this study, we develop a method to reveal hidden *Ta* waves in the majority of patients. Unlike in the direct processing of electrocardiogram (ECG) signals, an independent cycle of atria electrical activities is obtained by extracting the atrial part of transmembrane potentials (TMPs) of whole-heart activities and mapping partial TMPs back onto the ECG of a body surface. Despite the simple back-forward process, obtaining *Ta* waves with acceptable quality is challenging in computational design because the accuracy of inverse transformation highly relies on the fidelity of organ geometries and conductivities [9, 10]. Constructing a 3D mesh model for each patient can be expensive, and thus, our solution for the inverse problem assumes that all the physical measures of an individual are unreliable, and these parameters should be re-estimated during computation. The algorithm only requires long observation periods to perform self-adaptive and reliable statistical inferences for unsure parameters.

The morphologies of *P* waves and hidden *Ta* waves are typically considered important among healthcare professionals. In [1] and [11], *Ta* waves were suggested to be an early sign of inferior injury, such as acute atrial infarction or arrhythmia. One example is atrial fibrillation if the waves can be observed during conduction system malfunction, such as long QT interval or atrioventricular block.

*P* waves alone are important in justifying heart diseases in both sinus rhythm and fibrillation [12]. However, *P* waves are relatively weak, and thus, current understanding on atria activities are generally insufficient. Few studies have investigated activities from the atria part of myocytes. In [13], the problem was initially addressed via computer simulation. Without involving inverse problems, [13] facilitated a forward model, which mapped current dipoles onto atrial mid-myocardium to surface ECGs under a set of predetermined parameters to understand the contributions of right and left atria activities to the observed *P* waves during sinus rhythm and atrial fibrillation.

In clinical situations, *Ta* waves are not only hidden, but *P* waves are also tangled with the adjacent QRS complex. In [7] and [8], statistical and signal processing methods were, respectively used to single out *P* waves from the rest of the QRS part. Surface signals are complex combinations of current stimuli from millions of cardiomyocytes; thus, the signal separation task must be performed at the level of myocardium cells, and solving an ill-posed inverse problem is inevitable. Pioneer source models, such as those in [14–16], have been integrated into and advanced for contemporary computational electrocardiology models that establish models for various processes ranging from cellular bioelectrical activities to body surface potential distribution [17–19]. The forward problems involve mapping from inter- to intra-cellular currents onto body surface potential distributions [20, 21]. Given the nature of complexity in the biological field, the large degree of freedom poses a challenge in evaluating inverse problems.

When dealing with the inverse problem, most regulation methods can only condition numerical difficulty from the mathematical point of view; however, the problem of multiple solutions must be addressed by restoring possible missing constraints. In [22], this difficulty was solved by reducing the inverse problem in limited mapping from epicardial to body surfaces. In [23], the inverse problem was addressed by proposing a special equipment that could collect thousands of potentials on the body surface instead of the standard 12-lead ECG. Activation time sequencing or imaging can be evaluated for primitive diagnosis without retrieving detailed cellular activities [24–28]. However, electric current constraints from the ionic behavior of individual cells, such as in [17], impose additional computational challenge given the millions of myocytes. Moreover, additional constraints introduce other unknown parameters, and the degree of freedom remains high. Therefore, advanced statistical methods should be applied to obtain reliable solutions.

Existing studies continuously contribute to addressing the problem of physiological models, and specific body surface electrical data from patients are always being corrupted by noise and the incorrect construction of organ geometries. In [29], spatial covariance in a volume conductor was facilitated for maximum a posteriori (MAP) equation. In [30], temporal and spatial covariances were estimated under certain mathematical assumptions based on structures that were inherent in the space–time correlation matrix. In [31], the facilitation of multiple information sources to improve the efficiency of Bayesian MAP formulation was suggested. In [32], TMPs were constrained using a diffusion–reaction model from cellular activation dynamics, which limited the inverse problem in both spatial and temporal dimensions. This work suggested relying on a statistical method to address both model and data errors in terms of prior knowledge on cell current dynamics and evidence for surface potential data. In [33], the progress of statistical identification from the perspective of systems biology was reviewed.

## Methods

As mentioned earlier, the extraction of *P* waves should be conducted at the electric current level in myocardial sources. The model for the cardiac computational system comprises two parts according to the component guideline in [34]. The first part involves mapping between body surface potentials and intra-cellular TMPs. Evaluating TMPs is considered a difficult inverse problem given a potential map of a body surface [35, 36]. The second part aims to constrain the inverse problem, in which the constraint describes changes in TMPs in terms of electrical propagation between myocardia. Most electrophysiological models are diffusion–reaction systems [17, 36–38].

### Inverse problem

We first consider the forward problem from equivalent current–dipole sources to body surface potentials. The sources of bioelectric currents across cell membranes excite the movement of cardiomyocytes and induce potential fields, which can be detected via surface electrodes. The total current density is presented as \(\varvec{J}(\varvec{r}) = \varvec{J}_{s} (\varvec{r}) + \sigma \varvec{E}(\varvec{r})\), where \(\varvec{J}_{s}\) is the net source current density (\(A/m^{2}\)); \(\sigma\) is conductivity in homogenous dielectric media; and \(\varvec{E}\) is the electric field, which exhibits the relation \(\varvec{E} = - \nabla \varPhi\) for potential function \(\varPhi (\varvec{r})\). Vector fields are denoted as bold face symbols, such as current density \(\varvec{J}(\varvec{r})\), which is a vector field at location \(\varvec{r}\). The total current \(\nabla \cdot \varvec{J} = 0\) diverges without external current under quasi-static conditions. Thus, \(\nabla \cdot (\sigma \nabla \varPhi ) = \nabla \cdot \varvec{J}_{s}\), and the relation between measured potentials and heart sources is transformed into a Poisson equation. For cardiac volume \(V_{H}\), the potentials are primitively expressed as \(\varPhi (\varvec{r}) = \frac{1}{4\pi \sigma }\iiint_{{V_{H} }} {\varvec{J}_{s} (\varvec{r^{\prime}}) \cdot \nabla \left( {\frac{1}{{|\varvec{r} - \varvec{r^{\prime}}|}}} \right)d^{3} \varvec{r^{\prime}}}\).

The transfer matrix \(\varvec{L}\) is synthesized with the geometries and conductivities of the organs inside the thorax. The geometric coordinates are segmented and discretized via magnetic resonance imaging (MRI) or computed tomography for a specific patient. Given numerical sensitivity and unavoidable movement, the forward model may suffer from geometric errors and should be incorporated as a part of modeling [9, 41]. In [42], geometric errors were suggested to be overcome by using Bayesian MAP estimation or Kalman filtering with Gaussian geometric errors. In the present study, we do not rely on the accuracy of geometry and conductivity. We estimate the parameters along with the process of estimating TMPs [43, 44]. Bayesian estimation in error covariance enables performance analysis to statistically characterize solutions.

### Reaction–diffusion systems

*P*waves. Electrical propagation is captured using the reaction–diffusion system [37] with the same setting as that in [46, 47]. Considering the balance between precision and computation, a simple system is sufficient to constrain the ill-posed inverse problem. Therefore, we adopt the system from [37] as follows:

### Hierarchical estimation

The recursive estimation for the posterior probability density \(P(\varvec{x}_{k} |\phi_{1:k} )\) can be conceptually achieved in two steps. The forecast term \(P(\varvec{x}_{k} |\phi_{1:k - 1} )\) can be obtained through Chapman–Kolmogorov integration \(\mathop \smallint \nolimits P(\varvec{x}_{k} |\varvec{x}_{k - 1} )P(\varvec{x}_{k - 1} |\phi_{1:k - 1} )d\varvec{x}_{k - 1}\), given that the posterior \(P(\varvec{x}_{k - 1} |\phi_{1:k - 1} )\) is known from time \(k - 1\), and \(P(\varvec{x}_{k} |\varvec{x}_{k - 1} )\) is determined from the system equation. The current time posterior \(P(\varvec{x}_{k} |\phi_{1:k} )\) is updated using the Bayes rule \(\frac{{P\left( {\phi_{k} |\varvec{x}_{k} } \right)P\left( {\varvec{x}_{k} |\phi_{1:k - 1} } \right)}}{{P\left( {\phi_{k} |\phi_{1:k - 1} } \right)}}\), where \(P(\phi_{k} |\phi_{1:k - 1} ) = \mathop \smallint \nolimits P(\phi_{k} |\varvec{x}_{k} )P(\varvec{x}_{k} |\phi_{1:k - 1} )d\varvec{x}_{k}\).

Following the suggestion in [47], a Monte Carlo Markov chain (MCMC) slice sampler [48] is applied in the Bayesian computation model because of the high dimension in our complex problem. A full Bayesian analysis of this problem is achieved by sampling the joint posterior distribution (13) using an MCMC technique called slice sampling [49]. Another potential solution for reducing the constraining effects of prior knowledge is the simultaneous estimation of the TMP dynamics and electrophysiological properties of the myocardium. This method has the advantage that the constraining models can be modified according to the collected data of patients with filtering of unknown parameters.

### Experiment setup

The geometry of a torso was adopted from the PhysioNet data archive, which also originated from the body surface mapping data of Dalhousie University [51–53]. Although accuracy is not a concern, mapping between surface nodes to the electrode positions of standard leads should be specified. Given the well-prepared recording and documentation in the data set, the detailed mapping from the surface nodes to the 15 standard leads was elaborated.

The ECG data were also adopted from PhysioNet: ptbdb and incartdb [51]. The signals were preprocessed to eliminate electromagnetic interference, baseline wandering (e.g., electromyographic noise), and various artifacts (e.g., electrode motion) [54].

The implementation programs for the experiments were developed in MATLAB and R. The transfer matrix was produced using the open source SCIRun/BioPSE from the Scientific Computing and Imaging Institute of the University of Utah [18, 55].

## Discussions

The error caused by spatial digitization may overwhelm other sources of error terms because of the resolution limitation in mesh grids. The errors caused by modeling, geometric, and measuring uncertainties are integrated into a reaction–diffusion system, which is a stochastic state-space system with a powerful capability to estimate true values among noise background. The two-variable reaction–diffusion system (2) is linearized for convenient computation experiments. All the parameters are embedded into each individual matrix element, which will be estimated in the subsequent MCMC step. Although sophisticated and accurate models for cardiac electrophysiological dynamics can be used, the estimation results may not improve in our case.

One of the error sources originates from the inhomogeneity and anisotropy of conductivity among different torso tissues. To obtain an acceptable imaging of the heart surface, the anisotropy of the intracellular conductivity tensor along the fiber structure of myocardia must be considered. In this study, we adopt a real-time estimation approach and leave all the parameters to be estimated during the online training step.

As mentioned earlier, given the lack of spatial resolution, we do not expect the complete accuracy of the inverse and forward simulations. Although an ordinary Moore–Penrose pseudo-inverse can obtain a perfectly low root-mean-square error after performing backward and forward projections, the cardiac potential maps obtained from this pseudo-inverse are meaningless. Therefore, the computation objective is to reconstruct a TMP distribution that closely mimics a real TMP distribution [17]. We determine that the sequence and heterogeneity of TMPs throughout the myocardium significantly influence the final waveform of surface ECG.

The selection of prior distributions and initial conditions can influence the final results of the inverse solution when the Bayesian approach is used. During our MCMC estimation step, we perform our experiments by assuming that the noise and error terms are Gaussian distributed. Accordingly, we also select Gaussian as their conjugate prior for assigning prior distribution. The Gaussian assumption is reasonable for the coarse computation level. If further investigation is conducted, then non-Gaussian or empirical distribution can be used to estimate unknown parameters under the hierarchical Bayesian estimation framework.

By contrast, initial conditions are considerably sensitive to the algorithm. If parameters or matrices are available in published studies, then we will use their well-known values as initial amounts. Patient-specific variation should be accommodated by the random sampling process. A good initial condition can make convergence fast and accurate. In certain situations, convergence may not reach each selected initial value. For example, numerical difficulty may occur for a certain prior distribution. Several procedures can improve numerical stability, including increasing the condition numbers of covariance matrices. However, their details are considered state-of-art and require further elaboration.

For computational convenience, the reaction–diffusion system (2) is linearized into a matrix equation because matrix operations are more efficient than the *n*th order Runge–Kutta methods with nonlinear cross-terms. However, the order of the system will increase significantly from nonlinear to linear. We artificially augment a large number of co-state variables in system (2); that is, we add ten times more dummy variables following variable \(\varvec{v}\). The number of dummy variables significantly influences the quality of TMP reconstruction.

## Conclusions

In this study, we consider only signals in the standard 12-lead ECG measurement. Unlike most studies, our approach relies on neither accurate geometries nor high-resolution body surface potential maps. We only use a rough model as an initial guess and then simultaneously estimate the inverse solutions and model parameters. In this setup, we eliminate the requirement to construct a geometric model of the heart using off-line MRI pictures. Our method simultaneously estimates inverse potential solutions and geometric details.

We separate the atrial part of ECG waves from the rest of the signals based on a standard clinical ECG recording. Traditional approaches [7, 8] work only on signal morphology and attempt to cancel out the large QRS-T complex from several standard templates. We proceed directly through TMPs in the myocardial cells and extract the electricity generated by atria only. The obtained surface ECG contains only *P* waves and is more robust than those obtained via waveform approach are.

Future development can focus on improving prior distribution specifications and initial conditions. We will additionally construct empirical distribution for all error terms to be a basis of the prior distribution. The algorithm will be continuously improved and parallelized for better computational efficiency.

## Declarations

### Declarations

**Authors’ contributions** YJC designed the study, performed the experiments, analyzed the data, and wrote the manuscript. WHT and WHH participated in designing the study, analyzing the data, and editing the manuscript. All authors read and approved the final manuscript.

### Authors’ information

Dr. Wei-Hua Tang is with Division of Cardiology, Department of Internal Medicine, National Yang-Ming University Hospital, Yilan, Taiwan.

Wen-Hsien Ho is currently a Professor with Healthcare Administration and Medical Informatics, Kaohsiung Medical University, Taiwan. His research interests include intelligent systems control, computational intelligence, robust control, and quality engineering.

Yenming J. Chen received his Ph.D. in Systems Science and Mathematics in 1998 at Washington University in St. Louis. He is currently a professor of Logistics Management, National Kaohsiung University of Science & Technology in Kaohsiung, Taiwan. His research interests include industrial data analytics, IoT/AI, agriculture supply chain, and game theory.

### Acknowledgements

We thank the anonymous reviewers and editors for their careful reading of our manuscript and their many insightful comments and suggestions.

### Competing interests

The authors declare that they have no competing interests.

### Availability of data and materials

The datasets generated and/or analysed during the current study are available in the PhysioNet repository, http://www.physionet.org/physiobank/database/.

### Consent for publication

Not applicable.

### Ethics approval and consent to participate

Not applicable.

### Funding

Publication of this article was funded by the Ministry of Science and Technology, Taiwan, under Grant MOST 105-2221-E-037-004, 104-2410-H-327-018-MY3, 106-2221-E-037-001, 106-2622-E-037-005-CC3, 106-2218-E-327-001, 107-2221-E-037-006, 107-2218-E-992-308, and the “Intelligent Manufacturing Research Center” (iMRC) from the Featured Areas Research Center Program within the framework of the Higher Education Sprout Project by the Ministry of Education (MOE) in Taiwan.

### About this supplement

This article has been published as part of BioMedical Engineering OnLine Volume 17 Supplement 2, 2018: Proceedings of the International Conference on Biomedical Engineering Innovation (ICBEI) 2016. The full contents of the supplement are available online at https://biomedical-engineering-online.biomedcentral.com/articles/supplements/volume-17-supplement-2.

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access**This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Ihara Z, van Oosterom A, Hoekema R. Atrial repolarization as observable during the PQ interval. J Electrocardiol. 2006;39(3):290–7.View ArticleGoogle Scholar
- Holmqvist F, Carlson J, Platonov PG. Detailed ecg analysis of atrial repolarization in humans. Ann Noninvasive Electrocardiol. 2009;14(1):13–8.View ArticleGoogle Scholar
- Stridh M, Sörnmo L. Spatiotemporal qrst cancellation techniques for analysis of atrial fibrillation. IEEE Trans Biomed Eng. 2001;48(1):105–11.View ArticleGoogle Scholar
- Rieta JJ, Castells F, Sánchez C, Zarzoso V, Millet J. Atrial activity extraction for atrial fibrillation analysis using blind source separation. IEEE Trans Biomed Eng. 2004;51(7):1176–86.View ArticleGoogle Scholar
- Lemay M, Vesin J-M, Van Oosterom A, Jacquemet V, Kappenberger L. Cancellation of ventricular activity in the ECG: evaluation of novel and existing methods. IEEE Trans Biomed Eng. 2007;54(3):542–6.View ArticleGoogle Scholar
- Rieta JJ, Alcaraz R. Applications of signal analysis to atrial fibrillation. In: Atrial fibrillation-mechanisms and treatment. New York: InTech; 2013. p. 155–180.Google Scholar
- Donoso FI, Figueroa RL, Lecannelier EA, Pino EJ, Rojas AJ. Atrial activity selection for atrial fibrillation ECG recordings. Comput Biol Med. 2013;43(10):1628–36.View ArticleGoogle Scholar
- Perlman O, Katz A, Weissman N, Amit G, Zigel Y. Atrial electrical activity detection using linear combination of 12-lead ECG signals. IEEE Trans Biomed Eng. 2014;61(4):1034–43.View ArticleGoogle Scholar
- Weixue L, Ling X. Computer simulation of epicardial potentials using a heart-torso model with realistic geometry. IEEE Trans Biomed Eng. 1996;43(2):211–7.View ArticleGoogle Scholar
- van Dam PM, van Oosterom A. Volume conductor effects involved in the genesis of the p wave. Europace. 2005;7(s2):S30–8.Google Scholar
- Childers R. Atrial repolarization: its impact on electrocardiography. J Electrocardiol. 2011;44(6):635–40.View ArticleGoogle Scholar
- Nielsen JB, Kühl JT, Pietersen A, Graff C, Lind B, Struijk JJ, Olesen MS, Sinner MF, Bachmann TN, Haunsø S, et al. P-wave duration and the risk of atrial fibrillation: results from the Copenhagen ECG study. Heart Rhythm. 2015;12(9):1887–95.View ArticleGoogle Scholar
- Jacquemet V. Modeling left and right atrial contributions to the ECG: a dipole-current source approach. Comput Biol Med. 2015;65:192–9.View ArticleGoogle Scholar
- FitzHugh R. Impulses and physiological states in theoretical models of nerve membrane. Biophys J. 1961;1(6):445.View ArticleGoogle Scholar
- Miller WT, Geselowitz DB. Simulation studies of the electrocardiogram. I. The normal heart. Circ Res. 1978;43(2):301–15.View ArticleGoogle Scholar
- Geselowitz DB. On the theory of the electrocardiogram. Proc IEEE. 1989;77(6):857–76.View ArticleGoogle Scholar
- Courtemanche M, Ramirez RJ, Nattel S. Ionic mechanisms underlying human atrial action potential properties: insights from a mathematical model. Am J Physiol Heart Circ Physiol. 1998;275(1):H301–21.View ArticleGoogle Scholar
- Nordsletten D, Niederer S, Nash M, Hunter P, Smith N. Coupling multi-physics models to cardiac mechanics. Prog Biophys Mol Biol. 2011;104(1):77–88.View ArticleGoogle Scholar
- Franzone PC, Pavarino LF, Scacchi S. Mathematical cardiac electrophysiology. 13th ed. Berlin: Springer; 2014.MATHGoogle Scholar
- Barr RC, Ramsey M III, Spach MS. Relating epicardial to body surface potential distributions by means of transfer coefficients based on geometry measurements. IEEE Trans Biomed Eng. 1977;24(1):1–11.View ArticleGoogle Scholar
- Gulrajani RM, Mailloux GE. A simulation study of the effects of torso inhomogeneities on electrocardiographic potentials, using realistic heart and torso models. Circ Res. 1983;52(1):45–56.View ArticleGoogle Scholar
- Horáček BM, Clements JC. The inverse problem of electrocardiography: a solution in terms of single-and double-layer sources on the epicardial surface. Math Biosci. 1997;144(2):119–54.MathSciNetMATHView ArticleGoogle Scholar
- Rudy Y. Noninvasive electrocardiographic imaging of arrhythmogenic substrates in humans. Circ Res. 2013;112(5):863–74.View ArticleGoogle Scholar
- Modre R, Tilg B, Fischer G, Wach P. Noninvasive myocardial activation time imaging: a novel inverse algorithm applied to clinical ECG mapping data. IEEE Trans Biomed Eng. 2002;49(10):1153–61.View ArticleGoogle Scholar
- He B, Li G, Zhang X. Noninvasive imaging of cardiac transmembrane potentials within three-dimensional myocardium by means of a realistic geometry anisotropic heart model. IEEE Trans Biomed Eng. 2003;50(10):1190–202.View ArticleGoogle Scholar
- Messnarz B, Tilg B, Modre R, Fischer G, Hanser F. A new spatiotemporal regularization approach for reconstruction of cardiac transmembrane potential patterns. IEEE Trans Biomed Eng. 2004;51(2):273–81.View ArticleGoogle Scholar
- Franzone PC, Pavarino L, Taccardi B. Simulating patterns of excitation, repolarization and action potential duration with cardiac bidomain and monodomain models. Math Biosci. 2005;197(1):35–66.MathSciNetMATHView ArticleGoogle Scholar
- Villongco CT, Krummen DE, Stark P, Omens JH, McCulloch AD. Patient-specific modeling of ventricular activation pattern using surface ECG-derived vectorcardiogram in bundle branch block. Prog Biophys Mol Biol. 2014;115(2):305–13.View ArticleGoogle Scholar
- van Oosterom A. The use of the spatial covariance in computing pericardial potentials. IEEE Trans Biomed Eng. 1999;46(7):778–87.View ArticleGoogle Scholar
- Greensite F. The temporal prior in bioelectromagnetic source imaging problems. IEEE Trans Biomed Eng. 2003;50(10):1152–9.View ArticleGoogle Scholar
- Serinagaoglu Y, Brooks DH, MacLeod RS. Improved performance of bayesian solutions for inverse electrocardiography using multiple information sources. IEEE Trans Biomed Eng. 2006;53(10):2024–34.View ArticleGoogle Scholar
- Wang L, Zhang H, Wong KC, Liu H, Shi P. Physiological-model-constrained noninvasive reconstruction of volumetric myocardial transmembrane potentials. IEEE Trans Biomed Eng. 2010;57(2):296–315.View ArticleGoogle Scholar
- Villaverde AF, Banga JR. Reverse engineering and identification in systems biology: strategies, perspectives and challenges. J R Soc Interface. 2014;11(91):20130505.View ArticleGoogle Scholar
- Lopez-Perez A, Sebastian R, Ferrero JM. Three-dimensional cardiac computational modelling: methods, features and applications. Biomed Eng Online. 2015;14(1):1.View ArticleGoogle Scholar
- Barr R, Pilkington T, Boineau J, Spach M. Determining surface potentials from current dipoles, with application to electrocardiography. IEEE Trans Biomed Eng. 1996;2:88–92.Google Scholar
- Lopez-Rincon A, Bendahmane M, Ainseba B. On 3D numerical inverse problems for the bidomain model in electrocardiology. Comput Math Appl. 2015;69(4):255–74.MathSciNetMATHView ArticleGoogle Scholar
- Aliev RR, Panfilov AV. A simple two-variable model of cardiac excitation. Chaos Solitons Fractals. 1996;7(3):293–301.View ArticleGoogle Scholar
- Hurtado DE, Kuhl E. Computational modelling of electrocardiograms: repolarisation and T-wave polarity in the human heart. Comput Methods Biomech Biomed Eng. 2014;17(9):986–96.View ArticleGoogle Scholar
- Seger M, Fischer G, Modre R, Messnarz B, Hanser F, Tilg B. Lead field computation for the electrocardiographic inverse problem—finite elements versus boundary elements. Comput Methods Programs Biomed. 2005;77(3):241–52.View ArticleGoogle Scholar
- Wang D, Kirby RM, Johnson CR. Resolution strategies for the finite-element-based solution of the ECG inverse problem. IEEE Trans Biomed Eng. 2010;57(2):220–37.View ArticleGoogle Scholar
- Jiang M, Xia L, Shou G, Wei Q, Liu F, Crozier S. Effect of cardiac motion on solution of the electrocardiography inverse problem. IEEE Trans Biomed Eng. 2009;56(4):923–31.View ArticleGoogle Scholar
- Serinagaoglu Y, Brooks DH, MacLeod RS. Bayesian solutions and performance analysis in bioelectric inverse problems. IEEE Trans Biomed Eng. 2005;52(6):1009–20.View ArticleGoogle Scholar
- Erem B, Coll-Font J, Orellana RM, Št P, Brooks DH, et al. Using transmural regularization and dynamic modeling for noninvasive cardiac potential imaging of endocardial pacing with imprecise thoracic geometry. IEEE Trans Med Imaging. 2014;33(3):726–38.View ArticleGoogle Scholar
- Rahimi A, Sapp J, Xu J, Bajorski P, Horacek M, Wang L. Examining the impact of prior models in transmural electrophysiological imaging: a hierarchical multiple-model bayesian approach. IEEE Trans Med Imaging. 2016;35(1):229–43.View ArticleGoogle Scholar
- Luo C-H, Rudy Y. A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes. Circ Res. 1994;74(6):1071–96.View ArticleGoogle Scholar
- Roth BJ. Electrical conductivity values used with the bidomain model of cardiac tissue. IEEE Trans Biomed Eng. 1997;44(4):326–8.View ArticleGoogle Scholar
- Rogers JM, McCulloch AD. A collocation-Galerkin finite element model of cardiac action potential propagation. IEEE Trans Biomed Eng. 1994;41(8):743–57.View ArticleGoogle Scholar
- Berliner LM. Hierarchical bayesian time series models. Maximum entropy and Bayesian methods. Berlin: Springer; 1996. p. 15–22.MATHView ArticleGoogle Scholar
- Wikle CK, Berliner LM. A bayesian tutorial for data assimilation. Physica D. 2007;230(1):1–16.MathSciNetMATHView ArticleGoogle Scholar
- Hadfield JD, et al. Mcmc methods for multi-response generalized linear mixed models: the MCMC glmm R package. J Stat Softw. 2010;33(2):1–22.View ArticleGoogle Scholar
- Parker SG, Johnson CR. Scirun: a scientific programming environment for computational steering. In: Proceedings of the 1995 ACM/IEEE conference on supercomputing. New York: ACM; 1995. p. 52.Google Scholar
- Van Oosterom A, Oostendorp T. Ecgsim: an interactive tool for studying the genesis of QRST waveforms. Heart. 2004;90(2):165–8.View ArticleGoogle Scholar
- Goldberger AL, Amaral LA, Glass L, Hausdorff JM, Ivanov PC, Mark RG, Mietus JE, Moody GB, Peng C-K, Stanley HE. Physiobank, physiotoolkit, and physionet components of a new research resource for complex physiologic signals. Circulation. 2000;101(23):e215–20.View ArticleGoogle Scholar
- Dawoud F, Wagner GS, Moody G, Horáček BM. Using inverse electrocardiography to image myocardial infarction-reflecting on the 2007 physionet/computers in cardiology challenge. J Electrocardiol. 2008;41(6):630–5.View ArticleGoogle Scholar
- Sapp JL, Dawoud F, Clements JC, Horáček BM. Inverse solution mapping of epicardial potentials: quantitative comparison to epicardial contact mapping. Circulation Arrhythmia Electrophysiol. 2012;5:1001–9.View ArticleGoogle Scholar