Construction of boundary element models in bioelectromagnetism

Multisensor electro- and magnetoencephalographic (EEG and MEG) as well as electro- and magnetocardiographic (ECG and MCG) recordings have been proved useful in noninvasively extracting information on bioelectric excitation. The anatomy of the patient needs to be taken into account, when excitation sites are localized by solving the inverse problem. In this work, a methodology has been developed to construct patient specific boundary element models for bioelectromagnetic inverse problems from magnetic resonance (MR) data volumes as well as from two orthogonal X-ray projections. The process consists of three main steps: reconstruction of 3-D geometry, triangulation of reconstructed geometry, and registration of the model with a bioelectromagnetic measurement system. The 3-D geometry is reconstructed from MR data by matching a 3-D deformable boundary element template to images. The deformation is accomplished as an energy minimization process consisting of image and model based terms. The robustness of the matching is improved by multi-resolution and global-to-local approaches as well as using oriented distance maps. A boundary element template is also used when 3-D geometry is reconstructed from X-ray projections. The deformation is first accomplished in 2-D for the contours of simulated, built from the template, and real X-ray projections. The produced 2-D vector field is back-projected and interpolated on the 3-D template surface. A marching cube triangulation is computed for the reconstructed 3-D geometry. Thereafter, a non-iterative mesh-simplification method is applied. The method is based on the Voronoi-Delaunay duality on a 3-D surface with discrete distance measures. Finally, the triangulated surfaces are registered with a bioelectromagnetic measurement utilizing markers. More than fifty boundary element models have been successfully constructed from MR images using the methods developed in this work. A simulation demonstrated the feasibility of X-ray reconstruction; some practical problems of X-ray imaging need to be solved to begin tests with real data.reviewe

Non-invasive Localization of the Ventricular Excitation Origin Without Patient-specific Geometries Using Deep Learning

Ventricular tachycardia (VT) can be one cause of sudden cardiac death affecting 4.25 million persons per year worldwide. A curative treatment is catheter ablation in order to inactivate the abnormally triggering regions. To facilitate and expedite the localization during the ablation procedure, we present two novel localization techniques based on convolutional neural networks (CNNs). In contrast to existing methods, e.g. using ECG imaging, our approaches were designed to be independent of the patient-specific geometries and directly applicable to surface ECG signals, while also delivering a binary transmural position. One method outputs ranked alternative solutions. Results can be visualized either on a generic or patient geometry. The CNNs were trained on a data set containing only simulated data and evaluated both on simulated and clinical test data. On simulated data, the median test error was below 3mm. The median localization error on the clinical data was as low as 32mm. The transmural position was correctly detected in up to 82% of all clinical cases. Using the ranked alternative solutions, the top-3 median error dropped to 20mm on clinical data. These results demonstrate a proof of principle to utilize CNNs to localize the activation source without the intrinsic need of patient-specific geometrical information. Furthermore, delivering multiple solutions can help the physician to find the real activation source amongst more than one possible locations. With further optimization, these methods have a high potential to speed up clinical interventions. Consequently they could decrease procedural risk and improve VT patients' outcomes.Comment: 14 pages, 9 figures. Abstract was shortened for arXi

Non-Invasive Electrocardiographic Imaging of Ventricular Activities: Data-Driven and Model-Based Approaches

Die vorliegende Arbeit beleuchtet ausgewählte Aspekte der Vorwärtsmodellierung, so zum Beispiel die Simulation von Elektro- und Magnetokardiogrammen im Falle einer elektrisch stillen Ischämie sowie die Anpassung der elektrischen Potentiale unter Variation der Leitfähigkeiten. Besonderer Fokus liegt auf der Entwicklung neuer Regularisierungsalgorithmen sowie der Anwendung und Bewertung aktuell verwendeter Methoden in realistischen in silico bzw. klinischen Studien

Electrocardiographic Imaging Using a Spatio-Temporal Basis of Body Surface Potentials - Application to Atrial Ectopic Activity

Electrocardiographic imaging (ECGI) strongly relies on a priori assumptions and additional information to overcome ill-posedness. The major challenge of obtaining good reconstructions consists in finding ways to add information that effectively restricts the solution space without violating properties of the sought solution. In this work, we attempt to address this problem by constructing a spatio-temporal basis of body surface potentials (BSP) from simulations of many focal excitations. Measured BSPs are projected onto this basis and reconstructions are expressed as linear combinations of corresponding transmembrane voltage (TMV) basis vectors. The novel method was applied to simulations of 100 atrial ectopic foci with three different conduction velocities. Three signal-to-noise ratios (SNR) and bases of six different temporal lengths were considered. Reconstruction quality was evaluated using the spatial correlation coefficient of TMVs as well as estimated local activation times (LAT). The focus localization error was assessed by computing the geodesic distance between true and reconstructed foci. Compared with an optimally parameterized Tikhonov-Greensite method, the BSP basis reconstruction increased the mean TMV correlation by up to 22, 24, and 32% for an SNR of 40, 20, and 0 dB, respectively. Mean LAT correlation could be improved by up to 5, 7, and 19% for the three SNRs. For 0 dB, the average localization error could be halved from 15.8 to 7.9 mm. For the largest basis length, the localization error was always below 34 mm. In conclusion, the new method improved reconstructions of atrial ectopic activity especially for low SNRs. Localization of ectopic foci turned out to be more robust and more accurate. Preliminary experiments indicate that the basis generalizes to some extent from the training data and may even be applied for reconstruction of non-ectopic activity

Electrocardiogram Signal Analysis and Simulations for Non-Invasive Diagnosis - Model-Based and Data-Driven Approaches for the Estimation of Ionic Concentrations and Localization of Excitation Origins

Das Elektrokardiogramm (EKG) ist die Standardtechnik zur Messung der elektrischen Aktivität des Herzens. EKG-Geräte sind verfügbar, kostengünstig und erlauben zudem eine nichtinvasive Messung. Das ist insbesondere wichtig für die Diagnose von kardiovaskulären Erkrankungen (KVE). Letztere sind mit verursachten Kosten von 210 Milliarden Euro eine der Hauptbelastungen für das Gesundheitssystem in Europa und dort der Grund für 3,9 Millionen Todesfälle – dies entspricht 45% aller Todesfälle. Neben weiteren Risikofaktoren spielen chronische Nierenerkrankungen und strukturelle Veränderungen des Herzgewebes eine entscheidende Rolle für das Auftreten von KVE. Deshalb werden in dieser Arbeit zwei Pathologien, die in Verbindung zu KVE stehen, betrachtet: Elektrolytkonzentrationsveränderungen bei chronisch Nierenkranken und ektope Foki, die autonom Erregungen iniitieren. In beiden Projekten ist die Entwicklung von Methoden mithilfe von simulierten Signalen zur Diagnoseunterstützung das übergeordnete Ziel. Im ersten Projekt helfen simulierte EKGs die Signalverarbeitungskette zur EKG-basierten Schätzung der Ionenkonzentrationen von Kalium und Calcium zu optimieren. Die Erkenntnisse dieser Optimierung fließen in zwei patienten-spezifische Methoden zur Kaliumkonzentrationsschätzung ein, die wiederum mithilfe von Patientendaten ausgewertet werden. Die Methoden lieferten im Mittel einen absoluten Fehler von 0,37 mmol/l für einen patienten-spezifischen Ansatz und 0,48 mmol/l für einen globalen Ansatz mit zusätzlicher patienten-spezifischer Korrektur. Die Vorteile der Schätzmethoden werden gegenüber bereits existierender Ansätze dargelegt. Alle entwickelten Algorithmen sind ferner unter einer Open-Source-Lizenz veröffentlicht. Das zweite Projekt zielte auf die Lokalisierung von ektopen Foki mithilfe des EKGs ohne die Nutzung der individuellen Patientengeometrie. 1.766.406 simulierte EKG-Signale (Body Surface Potential Maps (BSPMs)) wurden zum Trainieren von zwei Convolutional Neural Networks (CNNs) erzeugt. Das erste CNN sorgt für die Schätzung von Anfang und Ende der Depolarisation der Ventrikel. Das zweite CNN nutzt die Information der Depolarisation im BSPM zur Schätzung des Erregungsurpsrungs. Der spezielle Aufbau des CNNs ermöglicht die Darstellung mehrerer Lösungen, wie sie durch Mehrdeutigkeiten im BSPM vorliegen können. Der kleinste Median des Lokalisierungsfehlers lag bei 1,54 mm für den Test-Datensatz der simulierten Signale, bzw. bei 37 mm für Patientensignale. Somit erlaubt die Kombination beider CNNs die verlässliche Lokalisierung von ektopen Foki auch anhand von Patientendaten, obwohl Patientendaten vorher nicht im Training genutzt wurden. Die Resultate dieser zwei Projekte demonstrieren, wie EKG-Simulationen zur Entwicklung und Verbesserung von EKG-Signalverarbeitungsmethoden eingesetzt werden und bei der Diagnosefindung helfen können. Zudem zeigt sich das Potential der Kombination von Simulationen und CNNs, um einerseits die zumeist raren klinischen Signale zu ersetzen und andererseits Modelle zu finden, die für mehrere Patienten/-innen gültig sind. Die vorgestellten Methoden bergen die Möglichkeit, die Diagnosestellungen zu beschleunigen und mit hoher Wahrscheinlichkeit den Therapieerfolg der Patienten zu verbessern

Spatiotemporal Sparse Bayesian Learning with Applications to Compressed Sensing of Multichannel Physiological Signals

Energy consumption is an important issue in continuous wireless telemonitoring of physiological signals. Compressed sensing (CS) is a promising framework to address it, due to its energy-efficient data compression procedure. However, most CS algorithms have difficulty in data recovery due to non-sparsity characteristic of many physiological signals. Block sparse Bayesian learning (BSBL) is an effective approach to recover such signals with satisfactory recovery quality. However, it is time-consuming in recovering multichannel signals, since its computational load almost linearly increases with the number of channels. This work proposes a spatiotemporal sparse Bayesian learning algorithm to recover multichannel signals simultaneously. It not only exploits temporal correlation within each channel signal, but also exploits inter-channel correlation among different channel signals. Furthermore, its computational load is not significantly affected by the number of channels. The proposed algorithm was applied to brain computer interface (BCI) and EEG-based driver's drowsiness estimation. Results showed that the algorithm had both better recovery performance and much higher speed than BSBL. Particularly, the proposed algorithm ensured that the BCI classification and the drowsiness estimation had little degradation even when data were compressed by 80%, making it very suitable for continuous wireless telemonitoring of multichannel signals.Comment: Codes are available at: https://sites.google.com/site/researchbyzhang/stsb

Proceedings of the second "international Traveling Workshop on Interactions between Sparse models and Technology" (iTWIST'14)

The implicit objective of the biennial "international - Traveling Workshop on Interactions between Sparse models and Technology" (iTWIST) is to foster collaboration between international scientific teams by disseminating ideas through both specific oral/poster presentations and free discussions. For its second edition, the iTWIST workshop took place in the medieval and picturesque town of Namur in Belgium, from Wednesday August 27th till Friday August 29th, 2014. The workshop was conveniently located in "The Arsenal" building within walking distance of both hotels and town center. iTWIST'14 has gathered about 70 international participants and has featured 9 invited talks, 10 oral presentations, and 14 posters on the following themes, all related to the theory, application and generalization of the "sparsity paradigm": Sparsity-driven data sensing and processing; Union of low dimensional subspaces; Beyond linear and convex inverse problem; Matrix/manifold/graph sensing/processing; Blind inverse problems and dictionary learning; Sparsity and computational neuroscience; Information theory, geometry and randomness; Complexity/accuracy tradeoffs in numerical methods; Sparsity? What's next?; Sparse machine learning and inference.Comment: 69 pages, 24 extended abstracts, iTWIST'14 website: http://sites.google.com/site/itwist1

The Application of Computer Techniques to ECG Interpretation

This book presents some of the latest available information on automated ECG analysis written by many of the leading researchers in the field. It contains a historical introduction, an outline of the latest international standards for signal processing and communications and then an exciting variety of studies on electrophysiological modelling, ECG Imaging, artificial intelligence applied to resting and ambulatory ECGs, body surface mapping, big data in ECG based prediction, enhanced reliability of patient monitoring, and atrial abnormalities on the ECG. It provides an extremely valuable contribution to the field

General rational approximation of Gaussian wavelet series and continuous-time gm-C filter implementation

© 2020 John Wiley & Sons, Ltd. This is the accepted version of the following article: Li, M, Sun, Y. General rational approximation of Gaussian wavelet series and continuous‐time g m ‐C filter implementation. Int J Circ Theor Appl. 2020; 1– 17., which has been published in final form at https://doi.org/10.1002/cta.2834.A general method of rational approximation for Gaussian wavelet series and Gaussian wavelet filter circuit design with simple gm-C integrators is presented in this work. Firstly, the multi-order derivatives of Gaussian function are analysed and proved as wavelet base functions. Then a high accuracy general approximation model of Gaussian wavelet series is constructed and the transfer function of first order derivative of Gaussian wavelet filter is obtained using quantum differential evolution (QDE) algorithm. Thirdly, as an example, a 5th order continuous-time analogue first order derivative of Gaussian wavelet filter circuit is designed based on multiple loop feedback structure with simple gm-C integrator as the basic blocks. Finally, simulation results demonstrate the proposed method is an excellent way for the wavelet transform implementation. The designed first order derivative of Gaussian wavelet filter circuit operates from a 0.53V supply voltage and a bias current 2.5nA. The power dissipation of the wavelet filter circuit at the basic scale is 41.1nW. Moreover, the high accuracy QRS detection based on the designed wavelet filter has been validated in application analysis.Peer reviewe