Computationally efficient noninvasive cardiac activation time imaging.

OBJECTIVE: The computer model-based computation of the cardiac activation sequence in humans has been recently subject of successful clinical validation. This method is of potential interest for guiding ablation therapy of arrhythmogenic substrates. However, computation times of almost an hour are unattractive in a clinical setting. Thus, the objective is the development of a method which performs the computation in a few minutes run time. METHODS: The computationally most expensive part is the product of the lead field matrix with a matrix containing the source pattern on the cardiac surface. The particular biophysical properties of both matrices are used for speeding up this operation by more than an order of magnitude. A conjugate gradient optimizer was developed using C++ for computing the activation map. RESULTS: The software was tested on synthetic and clinical data. The increase in speed with respect to the previously used Fortran 77 implementation was a factor of 30 at a comparable quality of the results. As an additional finding the coupled regularization strategy, originally introduced for saving computation time, also reduced the sensitivity of the method to the choice of the regularization parameter. CONCLUSIONS: As it was shown for data from a WPWpatient the developed software can deliver diagnostically valuable information at a much shorter span of time than current clinical routine methods. Its main application could be the localization of focal arrhythmogenic substrates.

A signal processing pipeline for noninvasive imaging of ventricular preexcitation.

OBJECTIVES: Noninvasive imaging of the cardiac activation sequence in humans could guide interventional curative treatment of cardiac arrhythmias by catheter ablation. Highly automated signal processing tools are desirable for clinical acceptance. The developed signal processing pipeline reduces user interactions to a minimum, which eases the operation by the staff in the catheter laboratory and increases the reproducibility of the results. METHODS: A previously described R-peak detector was modified for automatic detection of all possible targets (beats) using the information of all leads in the ECG map. A direct method was applied for signal classification. The algorithm was tuned for distinguishing beats with an adenosine induced AV-nodal block from baseline morphology in Wolff-Parkinson-White (WPW) patients. Furthermore, an automatic identification of the QRS-interval borders was implemented. RESULTS: The software was tested with data from eight patients having overt ventricular preexcitation. The R-peak detector captured all QRS-complexes with no false positive detection. The automatic classification was verified by demonstrating adenosine-induced prolongation of ventricular activation with statistical significance (p <0.001) in all patients. This also demonstrates the performance of the automatic detection of QRS-interval borders. Furthermore, all ectopic or paced beats were automatically separated from sinus rhythm. Computed activation maps are shown for one patient localizing the accessory pathway with an accuracy of 1 cm. CONCLUSIONS: The implemented signal processing pipeline is a powerful tool for selecting target beats for noninvasive activation imaging in WPW patients. It robustly identifies and classifies beats. The small beat to beat variations in the automatic QRS-interval detection indicate accurate identification of the time window of interest.

Lead field computation for the electrocardiographic inverse problem--finite elements versus boundary elements.

In order to be able to solve the inverse problem of electrocardiography, the lead field matrix (transfer matrix) has to be calculated. The two methods applied for computing this matrix, which are compared in this study, are the boundary element method (BEM) and the finite element method (FEM). The performance of both methods using a spherical model was investigated. For a comparable discretization level, the BEM yields smaller relative errors compared to analytical solutions. The BEM needs less computation time, but a larger amount of memory. Inversely calculated myocardial activation times using either the FEM or BEM computed lead field matrices give similar activation time patterns. The FEM, however, is also capable of considering anisotropic conductivities. This property might have an impact for future development, when also individual myocardial fiber architecture can be considered in the inverse formulation.

Ventricular surface activation time imaging from electrocardiogram mapping data.

Non-invasive imaging of cardiac electrophysiology provides a non-invasive way of obtaining information about electrical excitation. An iterative algorithm based on a general regularisation scheme for non-linear, ill-posed problems in Hilbert scales was applied to the electrocardiographic inverse problem, imaging the ventricular surface activation time (AT) map. This method was applied to electrocardiographic data from a 31-year-old healthy volunteer and a 24-year-old patient suffering from a Wolff-Parkinson-White (WPW) syndrome. The objective was to evaluate non-invasive AT imaging of an autonomous sinus rhythm and to quantify the localisation error of non-invasive AT imaging by localising the accessory pathway of the WPW syndrome and a pacing site for left ventricle pacing. The distances between the invasive and non-invasive localisation of the pacing site and the accessory pathway were 8 mm and 5 mm. The clinical case presented, shows that this non-invasive AT imaging approach may enable the reconstruction of single focal events with sufficient accuracy for potential clinical application.

An iterative algorithm for myocardial activation time imaging.

An iterative algorithm based on a general regularization scheme for nonlinear ill-posed problems in Hilbert scales (method A) is applied to the magnetocardiographic inverse problem imaging the surface myocardial activation time map. This approach is compared to an algorithm using an optimization routine for nonlinear ill-posed problems based on Tikhonov's approach of second order (method B). Method A showed good computational performance and the scheme for determining the proper regularization parameter lambda was found to be easier than in case of method B. The formulation is applied to magnetocardiographic recordings from a patient suffering from idiopathic ventricular tachycardia in which a sinus rhythm sequence was followed by a ventricular extrasystolic beat.

Application of high-order boundary elements to the electrocardiographic inverse problem.

Eight-noded quadrilateral boundary elements are applied to the electrocardiographic inverse problem as an example for high-order boundary elements. It is shown that the choice of the shape functions used for approximation of the potentials has a remarkable influence on the solution obtained if the number of electrodes is smaller than the number of primary source points (under-determined equation system). Three different formulations are investigated considering a concentric spheres problem where an analytic solution is available: (a) the isoparametric formulation; (b) the quasi-first-order formulation; and (c) the pseudo-subparametric formulation as a new method. In a second step the pseudo-subparametric formulation (which provided the best results in the test problem) is applied to real word data. The transmembrane potential pattern of a 40 years old female suffering from severe heart failure and ventricular tachycardia after large anterior wall myocardial infarction is reconstructed for one time instant. Furthermore, an algorithm for the calculation of the transfer matrix is presented which avoids restrictions to the boundary element mesh caused by the placement of the electrodes.

Two-dimensional Fourier representation used in the bioelectric forward problem.

The objective of this paper is the application of two-dimensional discrete Fourier transformation for solving the integral equation of the bioelectric forward problem. Therefore, the potential, the source term, and the integral equation kernel are assumed to be sampled at evenly spaced intervals. Thus the continuous functions of the problem domain can be expressed by their two-dimensional discrete Fourier transform in the spatial frequency domain. The method is applied to compute the surface potential generated by an eccentric dipole in a homogeneous spherical conducting medium. The integral equation for the potential is solved in the spatial frequency domain and the value of the potential at the sampling points is obtained from inverse Fourier transformation. The solution of the presented method is compared to both, an analytic solution and a solution gained from applying the boundary element method. Isoparametric quadrilateral boundary elements are used for modeling the spherical volume conductor in the boundary element solution, while in the two-dimensional Fourier transformation method the volume conductor is represented by a parametric boundary surface approximation.

Innovative telemonitoring system for cardiology: from science to routine operation.

OBJECTIVE: Results of the Austrian MOBITEL (MOBIle phone based TELemonitoring for heart failure patients) trial indicate that home-based telemonitoring improves outcome of chronic heart failure (CHF) patients and reduces both frequency and duration of hospitalizations. Based on lessons learned, we assessed the weak points to clear the way for routine operations. METHODS: We analyzed the system with respect to recommendations of the ESC Guidelines and experiences gained throughout the trial to identify potential improvements. The following components have been identified: a patient terminal with highest usability, integrated way to document drug-intake and well-being, and automated event detection for worsening of CHF. As a consequence the system was extended by Near Field Communication (NFC) technology and by an event management tool. RESULTS: Usability evaluation with 30 adults (14f, median 51y. IQR[45-65]) showed that 21 (8f) were able to immediately operate the system after reading a step-by-step manual. Eight (6f) needed one time demonstration and one man (80y) failed to operate the blood pressure meter. Routine operation of the revised system started in March 2009. Within 9 months, 15 patients (4f, median 74y. IQR[71-83], all NYHA-III) transmitted 17,149 items. 43 events were detected because of body weight gain of more then 2kg within 2 days. 49 therapy adjustments were documented. Three patients stopped using the system, two (1f) because of non-compliance and one (m, 82y) because of death. Overall, the rate of adherence to daily data transfer was 78%. CONCLUSION: First results confirm the applicability of the revised telemonitoring system in routine operation.

A bidomain model based BEM-FEM coupling formulation for anisotropic cardiac tissue.

A hybrid boundary element method (BEM)/finite element method (FEM) approach is proposed in order to properly consider the anisotropic properties of the cardiac muscle in the magneto- and electrocardiographic forward problem. Within the anisotropic myocardium a bidomain model based FEM formulation is applied. In the surrounding isotropic volume conductor the BEM is adopted. Coupling is enabled by requesting continuity of the electric potential and the normal of the current density across the boundary of the heart. Here, the BEM part is coupled as an equivalent finite element to the finite element stiffness matrix, thus preserving in part its sparse property. First, continuous convergence of the coupling scheme is shown for a spherical model comparing the computed results to an analytic reference solution. Then, the method is extended to the depolarization phase in a fibrous model of a dog ventricle. A precomputed activation sequence obtained using a fine mesh of the heart was downsampled and used to calculate body surface potentials and extracorporal magnetic fields considering the anisotropic bidomain conductivities. Results are compared to those obtained by neglecting in part or totally (oblique or uniform dipole layer model) anisotropic properties. The relatively large errors computed indicate that the cardiac muscle is one of the major torso inhomogeneities.