Doctor of Philosophy

dissertationInverse Electrocardiography (ECG) aims to noninvasively estimate the electrophysiological activity of the heart from the voltages measured at the body surface, with promising clinical applications in diagnosis and therapy. The main challenge of this emerging technique lies in its mathematical foundation: an inverse source problem governed by partial differential equations (PDEs) which is severely ill-conditioned. Essential to the success of inverse ECG are computational methods that reliably achieve accurate inverse solutions while harnessing the ever-growing complexity and realism of the bioelectric simulation. This dissertation focuses on the formulation, optimization, and solution of the inverse ECG problem based on finite element methods, consisting of two research thrusts. The first thrust explores the optimal finite element discretization specifically oriented towards the inverse ECG problem. In contrast, most existing discretization strategies are designed for forward problems and may become inappropriate for the corresponding inverse problems. Based on a Fourier analysis of how discretization relates to ill-conditioning, this work proposes refinement strategies that optimize approximation accuracy o f the inverse ECG problem while mitigating its ill-conditioning. To fulfill these strategies, two refinement techniques are developed: one uses hybrid-shaped finite elements whereas the other adapts high-order finite elements. The second research thrust involves a new methodology for inverse ECG solutions called PDE-constrained optimization, an optimization framework that flexibly allows convex objectives and various physically-based constraints. This work features three contributions: (1) fulfilling optimization in the continuous space, (2) formulating rigorous finite element solutions, and (3) fulfilling subsequent numerical optimization by a primal-dual interiorpoint method tailored to the given optimization problem's specific algebraic structure. The efficacy o f this new method is shown by its application to localization o f cardiac ischemic disease, in which the method, under realistic settings, achieves promising solutions to a previously intractable inverse ECG problem involving the bidomain heart model. In summary, this dissertation advances the computational research of inverse ECG, making it evolve toward an image-based, patient-specific modality for biomedical research

Assessment of the equivalent dipole layer source model in the reconstruction of cardiac activation times on the basis of BSPMs produced by an anisotropic model of the heart

Contains fulltext : 191955.pdf (Publisher’s version ) (Open Access

Reduced-order modeling for cardiac electrophysiology. Application to parameter identification

A reduced-order model based on Proper Orthogonal Decomposition (POD) is proposed for the bidomain equations of cardiac electrophysiology. Its accuracy is assessed through electrocardiograms in various configurations, including myocardium infarctions and long-time simulations. We show in particular that a restitution curve can efficiently be approximated by this approach. The reduced-order model is then used in an inverse problem solved by an evolutionary algorithm. Some attempts are presented to identify ionic parameters and infarction locations from synthetic ECGs.Comment: No. RR-7811 (2011

DICOM for EIT

With EIT starting to be used in routine clinical practice [1], it important that the clinically relevant information is portable between hospital data management systems. DICOM formats are widely used clinically and cover many imaging modalities, though not specifically EIT. We describe how existing DICOM specifications, can be repurposed as an interim solution, and basis from which a consensus EIT DICOM ‘Supplement’ (an extension to the standard) can be writte

Bayesian Inference with Combined Dynamic and Sparsity Models: Application in 3D Electrophysiological Imaging

Data-driven inference is widely encountered in various scientific domains to convert the observed measurements into information that cannot be directly observed about a system. Despite the quickly-developing sensor and imaging technologies, in many domains, data collection remains an expensive endeavor due to financial and physical constraints. To overcome the limits in data and to reduce the demand on expensive data collection, it is important to incorporate prior information in order to place the data-driven inference in a domain-relevant context and to improve its accuracy. Two sources of assumptions have been used successfully in many inverse problem applications. One is the temporal dynamics of the system (dynamic structure). The other is the low-dimensional structure of a system (sparsity structure). In existing work, these two structures have often been explored separately, while in most high-dimensional dynamic system they are commonly co-existing and contain complementary information. In this work, our main focus is to build a robustness inference framework to combine dynamic and sparsity constraints. The driving application in this work is a biomedical inverse problem of electrophysiological (EP) imaging, which noninvasively and quantitatively reconstruct transmural action potentials from body-surface voltage data with the goal to improve cardiac disease prevention, diagnosis, and treatment. The general framework can be extended to a variety of applications that deal with the inference of high-dimensional dynamic systems

Doctor of Philosophy

dissertationComputational simulation has become an indispensable tool in the study of both basic mechanisms and pathophysiology of all forms of cardiac electrical activity. Because the heart is comprised of approximately 4 billion electrically active cells, it is not possible to geometrically model or computationally simulate each individual cell. As a result computational models of the heart are, of necessity, abstractions that approximate electrical behavior at the cell, tissue, and whole body level. The goal of this PhD dissertation was to evaluate several aspects of these abstractions by exploring a set of modeling approaches in the field of cardiac electrophysiology and to develop means to evaluate both the amplitude of these errors from a purely technical perspective as well as the impacts of those errors in terms of physiological parameters. The first project used subject specific models and experiments with acute myocardial ischemia to show that one common simplification used to model myocardial ischemia-the simplest form of the border zone between healthy and ischemic tissue-was not supported by the experimental results. We propose a alternative approximation of the border zone that better simulates the experimental results. The second study examined the impact of simplifications in geometric models on simulations of cardiac electrophysiology. Such models consist of a connected mesh of polygonal elements and must often capture complex external and internal boundaries. A conforming mesh contains elements that follow closely the shapes of boundaries; nonconforming meshes fit the boundaries only approximately and are easier to construct but their impact on simulation accuracy has, to our knowledge, remained unknown. We evaluated the impact of this simplification on a set of three different forms of bioelectric field simulations. The third project evaluated the impact of an additional geometric modeling error; positional uncertainty of the heart in simulations of the ECG. We applied a relatively novel and highly efficient statistical approach, the generalized Polynomial Chaos-Stochastic Collocation method (gPC-SC), to a boundary element formulation of the electrocardiographic forward problem to carry out the necessary comprehensive sensitivity analysis. We found variations large enough to mask or to mimic signs of ischemia in the ECG

Nanoparticle electrical impedance tomography

We have developed a new approach to imaging with electrical impedance tomography (EIT) using gold nanoparticles (AuNPs) to enhance impedance changes at targeted tissue sites. This is achieved using radio frequency (RF) to heat nanoparticles while applying EIT imaging. The initial results using 5-nm citrate coated AuNPs show that heating can enhance the impedance in a solution containing AuNPs due to the application of an RF field at 2.60 GHz

Estimation of thorax shape for forward modelling in lungs EIT

The thorax models for pre-term babies are developed based on the CT scans from new-borns and their effect on image reconstruction is evaluated in comparison with other available models