- Research
- Open Access
- Published:

# Retrieving hidden atrial repolarization waves from standard surface ECGs

*BioMedical Engineering OnLine***volume 17**, Article number: 146 (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.

## 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,4,5,6,7,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,15,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,18,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,25,26,27,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,37,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}}}\).

To model equivalent current density, the entire myocardium is divided into grid meshes. Following the suggestion in [39], boundary element methods are applied. The potential \(\varPhi\) at the body surface is maintained as \(\varPhi\), and TMP is denoted as \(\varvec{u}\). By tessellating and vectorizing all cardiac and thorax surfaces, a discrete matrix Eq. (1) is obtained as suggested in [32] and [40].

where \(\varvec{L}\) is the discretized transfer matrix that converts TMP \(\varvec{u}\) to surface potential \(\phi_{8}\). When the vectorized body surface potentials are only sampled at eight electrode positions for the standard 12-lead ECG signals, the potentials are denoted as \(\varPhi_{8}\) for clarity.

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

Electrical propagation between myocardia is typically modeled differently in terms of complexity level—from the simplest Eikonal model at the tissue level, through bidomain/monodomain models and phenomenological models, to the most complicated ionic models at the cellular level. Phenomenological models focus at the macroscopic level and ranges from 2-variable equations [14, 37] to the complicated 15-variable Luo–Rudy model [45]. Resolution is not a concern in extracting *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:

where \(\varvec{u}\) and \(\varvec{v}\) are the column vectors of TMPs and recovery current, respectively; and the operator \(< , >\) represents a component-wise multiplication. \(D\) is the diffusion tensor; and \(k\), \(a\), and \(e\) are the parameters. By converting the equation into finite element meshes [47], the reaction–diffusion system can then be used as an effective constraint in solving the inverse problem. Let \(\varvec{x} = [\varvec{u},\,\varvec{v}]\). The system can then be written as \(\dot{\varvec{x}} = F_{d} (\varvec{x})\), where \(F_{d} (\varvec{x}) = \left[ {\left( {\nabla \left( {\varvec{D}\nabla \varvec{u}} \right) + k\varvec{u},\varvec{u} - \varvec{a},1 - \varvec{u} - < \varvec{u},\varvec{v} > } \right)\;and \; - e(\varvec{v} + k < \varvec{u},\varvec{u} - \varvec{a} - 1 > } \right]\).

### Hierarchical estimation

Our problem contains a large number of uncertainties, and thus, advanced Bayesian statistics can be a viable approach [44]. The basic idea is to estimate the posterior probability of the unknown cardiac source \(P(\varvec{x}_{k} |\phi_{1:k} )\) based on an a priori distribution of the sources \(P(\varvec{x})\) and a group of affecting parameters. When (1) and (2) are combined, we obtain the data model as follows (3):

where \(\varvec{H} = [\varvec{L} + \Delta \varvec{L}\;;\;0]\) is the output matrix with uncertainty \(\Delta \varvec{L}\), and \(\varvec{w}\) and \(\varvec{z}\) are two i.i.d. error processes with zero means and covariances \(\varvec{\xi}_{w}\) and \(\varvec{\xi}_{z}\). Given that the model does not rely on the accuracy of the heart and torso geometries, the error terms in the elements of the transfer matrix \(L\) are embedded into the matrix with random variables \(\Delta \varvec{L}\). Let \(\theta = (k,a,e)\) to incorporate the parameters in the reaction–diffusion function \(F_{d} ( \cdot )\). Therefore, the parameters for the process comprise \(\Delta \varvec{L}\) and \(\theta = (k,a,e)\).

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}\).

To deal with a large number of parameters, the guideline in [46] and [47] indicates that the complicated joint distribution in data model (3) can be formulated as a hierarchical model and factorized into a series of conditional distributions. The guideline suggests that the random variables to be estimated can be factored into three stages, such that \(p({\text{process}},{\text{parameters}}|{\text{data}}) \propto\) \(p({\text{data}}|{\text{process}},{\text{parameters}})\) \(p({\text{process}}|{\text{parameters}})\) \(p({\text{parameters}})\). Therefore, the joint posterior distribution can be written in a hierarchical form as follows:

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

To conduct the following experiments, 3D geometric models of a complete heart and torso are necessary. Cardiac geometric data were adopted from the ECGSim data set, which described a healthy normal young male using complete atria and ventricles (Fig. 1, with 1634 nodes for atria and 1500 nodes for ventricles) [50]. Given that a 3D imaging will not be constructed on the epicardial surface, the requirement for grid size is low. Resolution is further reduced to prevent the introduction of excessive numerical difficulties from the source of the standard 12-lead ECG.

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,52,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].

This study develops a model that retrieves hidden atrial repolarization waves by solving an inverse problem from surface ECG to cardiac TMPs (Fig. 2), where an ill-posed problem is constrained by temporal and spatial electrophysio-relations. The modeling approach can only be maintained at a coarse level because the source data are limited by the number of channels in the standard lead ECG. By contrast, cardiac electrical signals can be estimated by being modeled as a stochastic process with unknown excitation parameters and continuous acquisition of signals. In the solving process, several issues are encountered and need to discuss further.

The experiment presents good results. As shown in Fig. 3, the top panel presents the inverse solution for TMPs in the atrial part of the myocardium. The figure reflects the correct excitation sequence starting from the atrium to the end of the apex. When we multiply the entire TMPs to the transfer matrix, the forward problem restores the original ECG, as shown in the third panel. The figure exhibits good approximation of the original ECG (second panel), except for several ripples near the end of the cycle. This result is considered good because the resolution is under 14 nodes on the body surface and 20 nodes in the myocardium. The bottom panel shows the extracted atrial electric activities. Each line in the graph corresponds to one of the 14 nodes that constitute the standard 12-lead ECG.

## 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.

## References

- 1.
Ihara Z, van Oosterom A, Hoekema R. Atrial repolarization as observable during the PQ interval. J Electrocardiol. 2006;39(3):290–7.

- 2.
Holmqvist F, Carlson J, Platonov PG. Detailed ecg analysis of atrial repolarization in humans. Ann Noninvasive Electrocardiol. 2009;14(1):13–8.

- 3.
Stridh M, Sörnmo L. Spatiotemporal qrst cancellation techniques for analysis of atrial fibrillation. IEEE Trans Biomed Eng. 2001;48(1):105–11.

- 4.
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.

- 5.
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.

- 6.
Rieta JJ, Alcaraz R. Applications of signal analysis to atrial fibrillation. In: Atrial fibrillation-mechanisms and treatment. New York: InTech; 2013. p. 155–180.

- 7.
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.

- 8.
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.

- 9.
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.

- 10.
van Dam PM, van Oosterom A. Volume conductor effects involved in the genesis of the p wave. Europace. 2005;7(s2):S30–8.

- 11.
Childers R. Atrial repolarization: its impact on electrocardiography. J Electrocardiol. 2011;44(6):635–40.

- 12.
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.

- 13.
Jacquemet V. Modeling left and right atrial contributions to the ECG: a dipole-current source approach. Comput Biol Med. 2015;65:192–9.

- 14.
FitzHugh R. Impulses and physiological states in theoretical models of nerve membrane. Biophys J. 1961;1(6):445.

- 15.
Miller WT, Geselowitz DB. Simulation studies of the electrocardiogram. I. The normal heart. Circ Res. 1978;43(2):301–15.

- 16.
Geselowitz DB. On the theory of the electrocardiogram. Proc IEEE. 1989;77(6):857–76.

- 17.
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.

- 18.
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.

- 19.
Franzone PC, Pavarino LF, Scacchi S. Mathematical cardiac electrophysiology. 13th ed. Berlin: Springer; 2014.

- 20.
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.

- 21.
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.

- 22.
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.

- 23.
Rudy Y. Noninvasive electrocardiographic imaging of arrhythmogenic substrates in humans. Circ Res. 2013;112(5):863–74.

- 24.
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.

- 25.
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.

- 26.
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.

- 27.
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.

- 28.
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.

- 29.
van Oosterom A. The use of the spatial covariance in computing pericardial potentials. IEEE Trans Biomed Eng. 1999;46(7):778–87.

- 30.
Greensite F. The temporal prior in bioelectromagnetic source imaging problems. IEEE Trans Biomed Eng. 2003;50(10):1152–9.

- 31.
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.

- 32.
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.

- 33.
Villaverde AF, Banga JR. Reverse engineering and identification in systems biology: strategies, perspectives and challenges. J R Soc Interface. 2014;11(91):20130505.

- 34.
Lopez-Perez A, Sebastian R, Ferrero JM. Three-dimensional cardiac computational modelling: methods, features and applications. Biomed Eng Online. 2015;14(1):1.

- 35.
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.

- 36.
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.

- 37.
Aliev RR, Panfilov AV. A simple two-variable model of cardiac excitation. Chaos Solitons Fractals. 1996;7(3):293–301.

- 38.
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.

- 39.
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.

- 40.
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.

- 41.
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.

- 42.
Serinagaoglu Y, Brooks DH, MacLeod RS. Bayesian solutions and performance analysis in bioelectric inverse problems. IEEE Trans Biomed Eng. 2005;52(6):1009–20.

- 43.
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.

- 44.
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.

- 45.
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.

- 46.
Roth BJ. Electrical conductivity values used with the bidomain model of cardiac tissue. IEEE Trans Biomed Eng. 1997;44(4):326–8.

- 47.
Rogers JM, McCulloch AD. A collocation-Galerkin finite element model of cardiac action potential propagation. IEEE Trans Biomed Eng. 1994;41(8):743–57.

- 48.
Berliner LM. Hierarchical bayesian time series models. Maximum entropy and Bayesian methods. Berlin: Springer; 1996. p. 15–22.

- 49.
Wikle CK, Berliner LM. A bayesian tutorial for data assimilation. Physica D. 2007;230(1):1–16.

- 50.
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.

- 51.
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.

- 52.
Van Oosterom A, Oostendorp T. Ecgsim: an interactive tool for studying the genesis of QRST waveforms. Heart. 2004;90(2):165–8.

- 53.
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.

- 54.
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.

- 55.
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.

## 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.

## Author information

## Rights and permissions

**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.

## About this article

#### Published

#### DOI

### Keywords

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