F. M. Weber, D. U. J. Keller, S. Bauer, O. Dössel, G. Seemann, and C. Lorenz. Predicting tissue conductivity influences on body surface potentials-an efficient approach based on principal component analysis. In IEEE Transactions on Biomedical Engineering, vol. 58(2) , pp. 265-273, 2011
Abstract:
In this paper, we present an efficient method to estimate changes in forward-calculated body surface potential maps (BSPMs) caused by variations in tissue conductivities. For blood, skeletal muscle, lungs, and fat, the influence of conductivity variations was analyzed using the principal component analysis (PCA). For each single tissue, we obtained the first PCA eigenvector from seven sample simulations with conductivities between ±75% of the default value. We showed that this eigenvector was sufficient to estimate the signal over the whole conductivity range of ±75%. By aligning the origins of the different PCA coordinate systems and superimposing the single tissue effects, it was possible to estimate the BSPM for combined conductivity variations in all four tissues. Furthermore, the method can be used to easily calculate confidence intervals for the signal, i.e., the minimal and maximal possible amplitudes for given conductivity uncertainties. In addition to that, it was possible to determine the most probable conductivity values for a given BSPM signal. This was achieved by probing hundreds of different conductivity combinations with a numerical optimization scheme. In conclusion, our method allows to efficiently predict forward-calculated BSPMs over a wide range of conductivity values from few sample simulations.
Conduction velocity (CV) and CV restitution are important substrate parameters for understanding atrial arrhythmias. The aim of this work is to (i) present a simple but feasible method to measure CV restitution in-vivo using standard circular catheters, and (ii) validate its feasibility with data measured during incremental pacing. From five patients undergoing catheter ablation, we analyzed 8 datasets from sinus rhythm and incremental pacing sequences. Every wavefront was measured with a circular catheter and the electrograms were analyzed with a cosine-fit method that calculated the local CV. For each pacing cycle length, the mean local CV was determined. Furthermore, changes in global CV were estimated from the time delay between pacing stimulus and wavefront arrival. Comparing local and global CV between pacing at 500 and 300 ms, we found significant changes in 7 of 8 pacing sequences. On average, local CV decreased by 2015% and global CV by 1713%. The method allows for in-vivo measurements of absolute CV and CV restitution during standard clinical procedures. Such data may provide valuable insights into mechanisms of atrial arrhythmias. This is important both for improving cardiac models and also for clinical applications, such as characterizing arrhythmogenic substrates during sinus rhythm.
Atrial arrhythmias, such as atrial flutter or fibrillation, are frequent indications for catheter ablation. Recorded intracardiac electrograms (EGMs) are, however, mostly evaluated subjectively by the physicians. In this paper, we present a method to quantitatively extract the wave direction and the local conduction velocity from one single beat in a circular mapping catheter signal. We simulated typical clinical EGMs to validate the method. We then showed that even with noise, the average directional error was below 10(°) and the average velocity error was below 5.4 cm/s. In a realistic atrial simulation, the method could clearly distinguish between stimuli from different pulmonary veins. We further analyzed eight clinical data segments from three patients in normal sinus rhythm and with stimulation. We obtained stable wave directions for each segment and conduction velocities between 70 and 115 cm/s. We conclude that the method allows for easy quantitative analysis of single macroscopic wavefronts in intracardiac EGMs, such as during atrial flutter or in typical clinical stimulation procedures after termination of atrial fibrillation. With corresponding simulated data, it can provide an interface to personalize electrophysiological (EP) models. Furthermore, it could be integrated into EP navigation systems to provide quantitative data of high diagnostic value to the physician
Pyramidal GaAs structures on top of GaAs/AlAs distributed Bragg reflectors are investigated as candidates for true three-dimensional cavities with potentially low mode volume and high quality-factor. Different types of single and coupled resonators with base lengths of a few microns are realized using a combination of molecular-beam epitaxy, electron-beam lithography, and wet chemical etching. Embedded InGaAs quantum dots are utilized as light sources to verify the resonator modes. Furthermore, a spatially localized emission through the pyramid facets indicates the future possibility of coupling cavity modes to optical fibers. This could be interesting within the context of single photon emitters.
Aims Atrial cardiomyopathy (ACM) is associated with new-onset atrial fibrillation, arrhythmia recurrence after pulmonary vein isolation (PVI) and increased risk for stroke. At present, diagnosis of ACM is feasible by endocardial contact mapping of left atrial (LA) low-voltage substrate (LVS) or late gadolinium-enhanced magnetic resonance imaging, but their complexity limits a widespread use. The aim of this study was to assess non-invasive body surface electrocardiographic imaging (ECGI) as a novel clinical tool for diagnosis of ACM compared with endocardial mapping. Methods and results Thirty-nine consecutive patients (66 ± 9 years, 85% male) presenting for their first PVI for persistent atrial fibrillation underwent ECGI in sinus rhythm using a 252-electrode-array mapping system. Subsequently, high-density LA voltage and biatrial activation maps (mean 2090 ± 488 sites) were acquired in sinus rhythm prior to PVI. Freedom from arrhythmia recurrence was assessed within 12 months follow-up. Increased duration of total atrial conduction time (TACT) in ECGI was associated with both increased atrial activation time and extent of LA-LVS in endocardial contact mapping (r = 0.77 and r = 0.66, P < 0.0001 respectively). Atrial cardiomyopathy was found in 23 (59%) patients. A TACT value of 148 ms identified ACM with 91.3% sensitivity and 93.7% specificity. Arrhythmia recurrence occurred in 15 (38%) patients during a follow-up of 389 ± 55 days. Freedom from arrhythmia was significantly higher in patients with a TACT <148 ms compared with patients with a TACT ≥148 ms (82.4% vs. 45.5%, P = 0.019). Conclusion Analysis of TACT in non-invasive ECGI allows diagnosis of patients with ACM, which is associated with a significantly increased risk for arrhythmia recurrence following PVI.
Objectives: This study hypothesized that P-wave morphology and timing under left atrial appendage (LAA) pacing change characteristically immediately upon anterior mitral line (AML) block. Background: Perimitral flutter commonly occurs following ablation of atrial fibrillation and can be cured by an AML. However, confirmation of bidirectional block can be challenging, especially in severely fibrotic atria. Methods: The study analyzed 129 consecutive patients (66 ± 8 years, 64% men) who developed perimitral flutter after atrial fibrillation ablation. We designed electrocardiography criteria in a retrospective cohort (n = 76) and analyzed them in a validation cohort (n = 53). Results: Bidirectional AML block was achieved in 110 (85%) patients. For ablation performed during LAA pacing without flutter (n = 52), we found a characteristic immediate V1 jump (increase in LAA stimulus to P-wave peak interval in lead V1) as a real-time marker of AML block (V1 jump ≥30 ms: sensitivity 95%, specificity 100%, positive predictive value 100%, negative predictive value 88%). As V1 jump is not applicable when block coincides with termination of flutter, absolute V1 delay was used as a criterion applicable in all cases (n = 129) with a delay of 203 ms indicating successful block (sensitivity 92%, specificity 84%, positive predictive value 90%, negative predictive value 87%). Furthermore, an initial negative P-wave portion in the inferior leads was observed, which was attenuated in case of additional cavotricuspid isthmus ablation. Computational P-wave simulations provide mechanistic confirmation of these findings for diverse ablation scenarios (pulmonary vein isolation ± AML ± roof line ± cavotricuspid isthmus ablation). Conclusions: V1 jump and V1 delay are novel real-time electrocardiography criteria allowing fast and straightforward assessment of AML block during ablation for perimitral flutter.
Atrial fibrillation (AF) is the most common arrhythmia of the heart in industrialized countries. Its generation and the transitory behavior of paroxysmal AF are still not well understood. In this work we examine the interaction of two activation sources via an isthmus as possible cause for the initiation of fibrillation episodes. For this study, the electrophysiological model of Bueno-Orovio, Cherry and Fenton is adapted to atrial electrophysiology, both for physiological and electrophysiologically remodeled conditions due to AF. We show that the interaction of the pacemakers, combined with the geometrical constraints of the isthmus, can produce fibrillatory-type irregularities, which we quantify by the loss of spatial phase coherence in the transmembrane voltage. Transitions to irregular behavior occur when the frequencies of the pacemakers exceed certain thresholds, suggesting that AF episodes are initiated by frequency changes of the activating sources (sinus node, ectopic focus).
Atrial arrhythmias are frequently treated using catheter ablation during electrophysiological (EP) studies. However, success rates are only moderate and could be improved with the help of personalized simulation models of the atria. In this work, we present a workflow to generate and validate personalized EP simulation models based on routine clinical computed tomography (CT) scans and intracardiac electrograms. From four patient data sets, we created anatomical models from angiographic CT data with an automatic segmentation algorithm. From clinical intracardiac catheter recordings, individual conduction velocities were calculated. In these subject-specific EP models, we simulated different pacing maneuvers and measurements with circular mapping catheters that were applied in the respective patients. This way, normal sinus rhythm and pacing from a coronary sinus catheter were simulated. Wave directions and conduction velocities were quantitatively analyzed in both clinical measurements and simulated data and were compared. On average, the overall difference of wave directions was 15° (8%), and the difference of conduction velocities was 16 cm/s (17%). The method is based on routine clinical measurements and is thus easy to integrate into clinical practice. In the long run, such personalized simulations could therefore assist treatment planning and increase success rates for atrial arrhythmias.
This review article gives a comprehensive survey of the progress made in computa- tional modeling of the human atria during the last 10 years. Modeling the anatomy has emerged from simple peanut-like structures to very detailed models including atrial wall and fiber di- rection. Electrophysiological models started with just two cellular models in 1998. Today, five models exist considering e.g. details of intracellular compartments and atrial heterogeneity. On the pathological side, modeling atrial remodeling and fibrotic tissue are other important aspects. The bridge to data that are measured in the catheter laboratory and on the body surface (ECG) is under construction. Every measurement can be used either for model personalization or for validation. Potential clinical applications are briefly outlined and future research perspectives are suggested.
BACKGROUND: The prevalence of atrial fibrillation is increased in patients with end-stage renal disease. Previous studies suggested that extracellular electrolyte alterations caused by hemodialysis (HD) therapy could be proarrhythmic. METHODS: Multiscale models were used for a consequent analysis of the effects of extracellular ion concentration changes on atrial electrophysiology. Simulations were based on measured electrolyte concentrations from patients with end-stage renal disease. RESULTS: Simulated conduction velocity and effective refractory period are decreased at the end of an HD session, with potassium having the strongest influence. P-wave is prolonged in patients undergoing HD therapy in the simulation as in measurements. CONCLUSIONS: Electrolyte concentration alterations impact atrial electrophysiology from the action potential level to the P-wave and can be proarrhythmic, especially because of induced hypokalemia. Analysis of blood electrolytes enables patient-specific electrophysiology modeling. We are providing a tool to investigate atrial arrhythmias associated with HD therapy, which, in the future, can be used to prevent such complications.
D. U. J. Keller, F. M. Weber, G. Seemann, and O. Dössel. Ranking the Influence of Tissue Conductivities on Forward-Calculated ECGs. In IEEE Transactions on Biomedical Engineering, vol. 57(7) , pp. 1568-1576, 2010
Abstract:
This paper examined the effects that different tissue conductivities had on forward-calculated ECGs. To this end, we ranked the influence of tissues by performing repetitive forward calculations while varying the respective tissue conductivity. The torso model included all major anatomical structures like blood, lungs, fat, anisotropic skeletal muscle, intestine, liver, kidneys, bone, cartilage, and spleen. Cardiac electrical sources were derived from realistic atrial and ventricular simulations. The conductivity rankings were based on one of two methods: First, we considered fixed percental conductivity changes to probe the sensitivity of the ECG regarding conductivity alterations. Second, we set conductivities to the reported minimum and maximum values to evaluate the effects of the existing conductivity uncertainties. The amplitudes of both atrial and ventricular ECGs were most sensitive for blood, skeletal muscle conductivity and anisotropy as well as for heart, fat, and lungs. If signal morphology was considered, fat was more important whereas skeletal muscle was less important. When comparing atria and ventricles, the lungs had a larger effect on the atria yet the heart conductivity had a stronger impact on the ventricles. The effects of conductivity uncertainties were significant. Future studies dealing with electrocardiographic simulations should consider these effects.
Simulations of the electrophysiological behavior of the heart improve the comprehension of the mechanisms of the cardiovascular system. Furthermore, the mathematical modeling will support diagnosis and therapy of patients suffering from heart diseases. In this paper, the chain of modeling of the electrical function in the heart is described. The components are explained briefly, namely modeling of cardiac geometry, reconstructing the cardiac electrophysiology and excitation propagation. Additionally, the mathematical methods allowing to implement and solve these models are outlined. The three recently more investigated cases atrial fibrillation, ischemia and long-QT syndrome are described and show how cardiac modeling can support cardiologists in answering their open questions.
Experimental investigations of the nonlinear properties of superconducting niobium coplanar waveguide resonators are reported. The nonlinearity due to a current dependent kinetic inductance of the center conductor is strong enough to realize bifurcation of the nonlinear oscillator. When driven with two frequencies near the threshold for bifurcation, parametric amplification with a gain of +22.4 dB is observed.
Es wird eine Methode beschrieben, wie medizinische Bilder des Herzens modellbasiert mit EKG-Daten verknüpft werden können, um damit zu einer spezifischen Diagnostik und zu einer besseren Therapieplanung in der Kardiologie zu gelangen. Zunächst wird aus MRT- oder CT-Bildern des Patienten die Geometrie seines Herzens ermittelt. Elektrokardiographische Messungen an der Körperoberfläche (EKG oder Body Surface Potential Mapping) und aus dem Inneren des Herzens (intracardial mapping) werden aufgenommen und die Orte der Messung in den Bilddatensatz eingetragen (registration). Ein elektrophysiologisches Computermodell vom Herzen des Patienten wird mit Hilfe der elektrophysiologischen Messdaten iterativ angepasst. Schließlich entsteht im Computer ein virtuelles Herz des Patienten, welches sowohl die Geometrie als auch die Elektrophysiologie wiedergibt. Ein Modell der Vorhöfe hat beispielsweise das Potenzial, die Ursachen von Vorhofflimmern zu erkennen und die Radiofrequenz-Ablationsstrategie zu optimieren. Ein Modell der Ventrikel des Herzens kann helfen, genetisch bedingte Rhythmusstörungen besser zu verstehen oder auch die Parameter bei der kardialen Resynchronisationstherapie zu optimieren. Die Modellierung des Herzens mit einem Infarktgebiet könnte die elektrophysiologischen Auswirkungen des Infarktes beschreiben und die Risikostratifizierung für gefährliche ventrikuläre Arrhythmien unterstützen oder die Erfolgsrate bei ventrikulären Ablationen erhöhen.
Background: Catheter ablation of complex atrial arrhythmias, such as atrial fibrillation and atypical atrial flutter, is still challenging. Clinically evaluated ablation methods are leading to moderate success rates. Assessments of intracardiac electrograms are often done subjectively by the physician. Automatic algorithms can therefore improve the analysis of complex atrial electrograms (EGMs). In this work, we demonstrate a quantitative analysis of intracardiac EGMs from circular mapping catheters in humans. Both the wave direction and the local conduction velocity (CV) were calculated from individual wave fronts passing the catheter.Methods: Intracardiac EGMs measured with circular mapping catheters in humans were retrospectively analyzed. Five data sets from 3 patients undergoing catheter ablation of atrial fibrillation or flutter were available. Using a nonlinear energy operator, activation times from 9 bipolar catheter signals were calculated for each atrial activity. The resulting activation pattern was fitted to a cosine-shaped data model that has been validated in a previous simulation study. The cosine phase represented the wave direction. From the cosine amplitude and the catheter radius, the conduction velocity was calculated.Results: The wave directions in all five measurements were stable with a standard deviation below 10°. Calculated CVs were in the range of 70 to 110 cm/s, which is in accordance with published values. In one patient, electrograms were recorded during atrial stimulation. Stimulation cycle length was decreased from 500 to 300 milliseconds. Conduction velocity decreased by approximately 10% at a cycle length of 300 milliseconds compared with the CV at 500 milliseconds.Conclusion: The results show the ability to reliably extract wave direction and conduction velocity from intracardiac EGMs recorded with circular mapping catheters. Detected directions were stable, and the CV values were in a physiological range. As individual beats are analyzed, the method will also enable the quantitative study of singular events such as ectopic beats and facilitate the localization of tachycardia origins. Further, it will help to measure substrate parameters such as the CV and even CV restitution behavior. This way, the method can help to identify patient-specific physiological parameters that can be integrated into patient-specific models. Furthermore, it can directly provide quantitative data of high diagnostic value to the examiner and thereby improve clinical success rates.
F. Weber. Towards patient-specific simulations of atrial fibrillation. In VPH Events - 4th Cardiac Physiome Workshop, 2009
Abstract:
Patient-specific cardiac simulations are approaching clinical applications. They could for example improve the treatment of atrial fibrillation (AF). Currently, many patients suffering from AF are treated with minimally-invasive catheter ablation. Using this technique, trigger sources for AF (mainly the pulmonary veins), are electrically isolated from the rest of the atrium. However, a large set of different ablation strategies is currently used in clinical practice. Therefore, the choice of a certain ablation strategy as well as the probability for successful and sustained AF termination are strongly dependent on the experience of the cardiologist. Atrial simulations could assist the cardiologist in the choice of a suitable method for an individual patient. For this, the atrial models have to be adapted to the patient. Besides anatomical modeling, several challenges must be faced in this process. First, an appropriate model of cellular electrophysiology and excitation conduction must be chosen. The model must provide the necessary accuracy and at the same time be fast enough for clinical applications. As a trade-off between accuracy and speed, we propose a minimal model adapted to atrial electrophysiology. Second, a main problem is the adaptation of physiological parameters in the patient-specific model as well as its validation. Therefore, an interface between clinical data and the model is needed. Data collected in standard clinical workflow are mainly intracardiac catheter ECGs. We therefore present techniques to model such catheter measurements. Signals from both circular mapping catheters (such as Lasso or Orbiter) as well as Coronary Sinus catheters can be simulated and compared to clinical signals. These are important steps towards clinical applications of atrial models. The long-term goal then is to assist the cardiologist in the choice of the best treatment for an individual patient.
Patient-specific model adaptation and validation requires a comparison of simulations with measured patient data. For patients suffering from atrial fibrillation, such data is mainly available as intracardiac catheter signals. In this work, we demonstrate the simulation of clinically relevant catheter data as measured using circular mapping catheters (such as Lasso or Orbiter) and coronary sinus catheters using atrial simulations on a realistic geometry. Four circular catheters are modeled using a projection technique for two distinct types of application. We show that in sinus rhythm, the choice of a distinct electrophysiological model does not impair the signal quality. Finally, we compare simulated potentials to a real clinical measurement. In the future, with patient- specific models available, such comparisons can constitute an important interface for personalizing cardiac models.
F. M. Weber, C. Schilling, D. Straub, O. Dössel, G. Seemann, and C. Lorenz. Localizing ectopic foci in the pulmonary veins from intracardiac ECGs a simulation study. In IFMBE Proceedings World Congress on Medical Physics and Biomedical Engineering, vol. 25/4, pp. 645-648, 2009
Abstract:
Intracardiac catheter recordings are available in common clinical practice. They can therefore be employed to adapt and validate atrial computer models of individual patients. Hence, their information content needs to be analyzed quantitatively. During treatment of atrial arrhythmia such as atrial flutter or fibrillation, the location of ectopic foci in the pulmonary veins is of special interest. In this study, virtual catheter signals are extracted from an atrial simulation on a realistic geometry with normal sinus rhythm as well as ectopic stimuli in all four pulmonary veins. Using a simplified Pan-Tompkins algorithm, the activation times are determined. Based on the analysis of the activation sequence in a circular mapping catheter simulated on the posterior left atrial wall, all four ectopic foci can clearly be associated with the pulmonary vein they came from. For a catheter on the anterior wall, this is possible for three of the four ectopic beats. Despite the knowledge gathered for the personalization of patient models, such simulations may help cardiologists to better classify measured signals.
Patient-specific cardiac simulations are approaching clinical applications. They could for example improve the treatment of atrial fibrillation (AF). Currently, many patients suffering from AF are treated with minimally-invasive catheter ablation. Using this technique, trigger sources for AF (mainly the pulmonary veins), are electrically isolated from the rest of the atrium. However, a large set of different ablation strategies is currently used in clinical practice. Therefore, the choice of a certain ablation strategy as well as the probability for successful and sustained AF termination are strongly dependent on the experience of the cardiologist. Atrial simulations could assist the cardiologist in the choice of a suitable method for an individual patient. For this, the atrial models have to be adapted to the patient. Besides anatomical modeling, several challenges must be faced in this process. First, an appropriate model of cellular electrophysiology and excitation conduction must be chosen. The model must provide the necessary accuracy and at the same time be fast enough for clinical applications. As a trade-off between accuracy and speed, we propose a minimal model adapted to atrial electrophysiology. Second, a main problem is the adaptation of physiological parameters in the patient-specific model as well as its validation. Therefore, an interface between clinical data and the model is needed. Data collected in standard clinical workflow are mainly intracardiac catheter ECGs. We therefore present techniques to model such catheter measurements. Signals from both circular mapping catheters (such as Lasso or Orbiter) as well as Coronary Sinus catheters can be simulated and compared to clinical signals. These are important steps towards clinical applications of atrial models. The long-term goal then is to assist the cardiologist in the choice of the best treatment for an individual patient.
Simulation of cardiac excitation is often a trade-off between accuracy and speed. A promising minimal, time-efficient cell model with four state variables has recently been presented together with parametrizations for ventricular cell behaviour. In this work, we adapt the model parameters to reproduce atrial excitation properties as given by the Courtemanche model. The action potential shape is considered as well as the restitution of action potential duration and conduction velocity. Simulation times in a single cell and a tissue patch are compared between the two models. We further present the simulation of a sinus beat on the atria in a realistic 3D geometry using the fitted minimal model in a monodomain simulation.
Background: Perimitral flutter commonly occurs following ablation of atrial fibrillation (AF) and can be cured by an anterior mitral line (AML). However, confirmation of bidirectional block can be challenging. Objective: We hypothesized that P-wave morphology and timing under left atrial appendage (LAA) pacing changes upon AML- block. Methods: We analyzed 129 consecutive patients (66±8 y, 64%male) who developed perimitral flutter after AF ablation. We designed ECG-criteria in a retrospective cohort (n=76) and analyzed them in a validation cohort (n=53). Results: Bidirectional AML-block was achieved in 110 patients (85%). For ablation performed during LAA-pacing without flutter (n=52), we found an immediate V1-jump (increase in LAA- stimulus to P-wave peak in lead V1) as a real-time marker of AML-block (V1-jump ≥30ms: sensitivity 95%, specificity 100%, PPV 100%, NPV 88%). Since V1-jump is not applicable when block coincides with termination of flutter, absolute V1-delay was used as a criterion applicable in all cases (n=129) with a delay of 203ms indicating block (sensitivity 92%, specificity 84%, PPV 90%, NPV 87%). Furthermore, an initial negative P-wave portion in the inferior leads was observed, which was attenuated in case of additional cavotricuspid isthmus (CTI) ablation. Computational P-wave simulations provide mechanistic confirmation of these findings for diverse ablation scenarios (pulmonary vein isolation±AML±roof-line±CTI ablation). Conclusion: V1-jump and V1-delay are novel real-time ECG- criteria allowing fast and straightforward assessment of AML- block during ablation for perimitral flutter.
The arrhythmogenic mechanisms of atrial fibrillation (AF) are still not well understood. Increased atrial fibrosis is a structural hallmark in patients with persistent AF. We assessed the electrogram signature rotational activity and their spatial relationship to low voltage areas in patients with persistent AF. Computer simulations implicating 3- dimensional atrial tissue with different amount of atrial fibrosis were used to assess development and stability of rotational activities during AF. Rotor anchoring occurred at the borderzone between fibrosis and healthy atrial tissue with 12 consecutive rotations prior to rotor extinction. Rotational activity in fibrotic tissue resulted in fractionated signals and were overlapped with large negative electrograms in unipolar recording mode from neighboring healthy tissue impressing as a focal source. Necessary conditions for development and stability of rotational activities around fibrosis were on the one hand a minimum size of atrial fibrosis area equal or larger than 10mm x 10mm and on the other hand the degree of atrial fibrosis of 40%. Clinical data showed that AF termination sites were located within low voltage areas (displaying <0,5mV in AF on the multielectrode mapping catheter) in 80% and at their borderzones in 20% of cases.
A new method to predict changes in a lead-field matrix induced by conductivity variations of a single body tissue is proposed. The approach is based on the princi- ple component analysis (PCA) with three initial lead-field matrices transformed to vectors as input. For each tissue blood, lungs, muscles and fat a PCA was carried out. Further, for each tissue the default conductivity value and the conductivity varied by ±50 % were used to calculate the sample lead-field matrices. The results of the PCAs in- dicate that for every tissue the first principle component suffices to predict the conductivity-induced changes in the samples. With an interpolation of the scores we addition- ally show that the prediction is not bound to the sample ma- trices but moreover every matrix within each conductivity range is possibly estimated and conclusively predicted.
Abstract. Atrial fibrillation (AF) is the most common cardiac arrhyth- mia. Patient-specific computational modeling of the atria can provide a better understanding about mechanisms underlying the arrhythmia and will potentially be used for model-based ablation therapy evaluation and planning. Electrical excitation spreads from the left to the right atrium at discrete locations. The location of the muscular bridges cannot be determined from image data. In the present study, left atrial activation sources were manually identified in local activation time maps of 4 AF patients. This information was used to adjust rule-based placed intera- trial bridges in anatomical atrial models of the patients. Sinus rhythm simulations showed a better qualitative agreement to the measured left atrial activation patterns after the adjustment of the bridges. For one patient, the simulated body surface potential (BSP) pattern after the adjustment correlated better to measured BSP maps. The results show that the fusion of intracardiac electrical measurements of early left atrial activation can be used to refine patient atria models with information of the myocardial structure which cannot be imaged. In future, such personalized atrial models may be used to support EP interventions.
C. Lenk, F. M. Weber, M. Bauer, M. Einax, G. Seemann, and P. Maass. Paroxysmal atrial fibrillation caused by interaction of pacemakerwaves and reduced excitability: Insights from the Bueno-Orovio model adapted to atria. In Computing in Cardiology Conference (CinC), pp. 1079-1082, 2013
Abstract:
As possible cause for atrial fibrillation (AF) we study the influence of a reduced excitability on the interaction of pacemaker waves in the Bueno-Orovio model with parameters adapted to atrial electrophysiology (aBO). One of the two pacemakers represents the sinus node and the other one a self-excitatory source in the left atrium. The pacemakers are spatially separated and their waves get in contact via a small bridge. In previous studies based on the FitzHugh-Nagumo (FHN) model it was shown that three different types of irregular activation patterns can occur in this problem. In the aBO model adapted to physiological conditions only one type is observed because, different from the FHN model, a reduction of excitability due to high-frequency pacing does not occur. If the excitability is reduced in the aBO model, all types of irregularities are recovered and, in addition, a further type is found. Because transitions from regular to irregular behavior depend on the pacing frequency, our findings provide a possible explanation for the phenomenon of paroxysmal AF.
Whole organ scale patient specific biophysical simulations contribute to the understanding, diagnosis and treatment of complex diseases such as cardiac arrhythmia. However, many individual steps are required to bridge the gap from an anatomical scan to a personal- ized biophysical model. In biophysical modeling, differential equations are solved on spatial domains represented by volumetric meshes of high resolution and in model-based segmentation, surface or volume meshes represent the patients geometry. We simplify the personalization pro- cess by representing the simulation mesh and additional relevant struc- tures relative to the segmentation mesh. Using a surface correspondence preserving model-based segmentation algorithm, we facilitate the inte- gration of anatomical information into biophysical models avoiding a complex processing pipeline. In a simulation study, we observe surface correspondence of up to 1.6mm accuracy for the four heart chambers. We compare isotropic and anisotropic atrial excitation propagation in a personalized simulation.
The objective of personalised modelling of the atria is to improve comprehension of the etiology of atrial arrhythmias, to enable specific diagnosis and to optimise therapy. We start with CT or MR datasets and use adapted segmentation procedures to build a patient-specific 3D-model of the atria. Then we include fibre direction based on the rules of atrial anatomy. Work in progress is also considering late enhancement MRI in order to add areas of fibrotic tissue. Next we can use BSPM data of the P-wave and solve the inverse problem of ECG to get a hypothesis about the spread of depolarisation. Finally we use intracardiac catheter signals (e.g. using a circular catheter) to measure direction and conduction velocity of depolarisation waves (sinus rhythm, atrial flutter, or following stimulation). All this is integrated into a personalised model of the atria of an individual patient. Our next goal will be to properly add ablation lines into the model.The research leading to these results has partly received funding from the European Communitys Seventh Framework Programme (FP7/2007-2013) under grant agreement n 224495 (euHeart project).
A framework for step-by-step personalization of a computational model of human atria is presented. Beginning with anatomical modeling based on CT or MRI data, next fiber structure is superimposed using a rule-based method. If available, late-enhancement-MRI images can be considered in order to mark fibrotic tissue. A first estimate of individual electrophysiology is gained from BSPM data solving the inverse problem of ECG. A final adjustment of electrophysiology is realized using intracardiac measurements. The framework is applied using several patient data. First clinical application will be computer assisted planning of RF-ablation for treatment of atrial flutter and atrial fibrillation.
The delineation of anatomical structures in medical images can be achieved in an efficient and robust manner using statistical anatomical organ models, which has been demonstrated for an already considerable set of organs, including the heart. While it is possible to provide models with sufficient shape variability to cope, to a large extent, with inter-patient variability, as long as object topology is conserved, it is a fundamental problem to cope with topological organ variability. We address this by creating a set of model variants and selecting the most appropriate model variant for the patient at hand. We propose a hybrid method combining model-based image analysis with a guided region growing approach for automated anatomical variant selection and apply it to the left atrium in cardiac CT images. Concerning the human heart, the left atrium is the most variable sub-structure with a variable number of pulmonary veins draining into it. It is of large clinical interest in the context of atrial fibrillation and related interventions.
Atrial myofiber orientation is complex and has multiple discrete layers and bundles. A novel robust semi-automatic method to incorporate atrial anisotropy and heterogeneities into patient-specific models is introduced. The user needs to provide 22 distinct seed-points from which a network of auxiliary lines is constructed. These are used to define fiber orientation and myocardial bundles. The method was applied to 14 patient-specific volumetric models derived from CT, MRI and photographic data. Initial electrophysiological simulations show a significant influence of anisotropy and heterogeneity on the excitation pattern and P-wave duration (20.7% shortening). Fiber modeling results show good overall correspondence with anatomical data. Minor modeling errors are observed if more than four pulmonary veins exist in the model. The method is an important step towards creating realistic patient-specific atrial models for clinical applications.
Atrial fibre architecture has complex patterns of bundles and layers and is known to impact on atrial electrophysiology, especially in fast-conducting bundles like Crista Terminalis, Bachmanns bundle and pectinate muscles. Based on a priori knowledge of atrial fibre structure, we incorporated rule-based fibre orientation in seven volumetric models of human atria using a semi-automatic approach. We were able to introduce multiple layers of myofibres and regional heterogeneities of ion channels in the models. We evaluated the influence of complete atrial fibre architecture on multiple modelling scales. First, we simulated atrial excitation in the isotropic and anisotropic models using the model of Courtemanche et al. in combination with the monodomain approach. Second, we computed body surface potentials from the simulated transmembrane voltages and compared these to measured ECGs from the respective patients. Temporal behaviour of the atrial excitation sequences was significantly altered in the anisotropic models compared to the sequences in the isotropic models. Complete atrial activation was achieved approximately 20% faster in the anisotropic models mostly due to fast conducting myofibre bundles. Electrophysiological heterogeneities influenced right atrial transmembrane voltage distribution over time due to a less negative action potential plateau in Crista Terminalis cells. P-wave duration was significantly shorted by the introduction of atrial anisotropy and the error to measured P-wave duration was reduced. Furthermore, a pattern change in body surface potential distribution over time was observed. The anisotropic patterns showed a better match to the measurements. Thus, the modelling error by using generalised fibre architecture for patient-specific models was smaller than by using isotropic models. The results highlight the necessity to incorporate atrial anisotropy in personalised models to produce more realistic simulations. The semi-automatic approach allows the use of these models for future clinical applications.
Background: Patients with end-stage renal disease show an increased prevalence of atrial fibrillation. A combined simulation and electrocardio- gram analysis study revealed a correlation between the changes in plasma electrolytes and intra-atrial conduction velocity related to hemodialysis (HD) session. A recognized limitation of the study is that simulations were performed on single-cell level. We present a computer study to investigate the influence of HD-related electrolyte modifications on atrial electrophys- iology in a volumetric environment.Methods: Based on the Courtemanche-Ramirez-Nattel model and its parameterization for different atrial tissues, we studied action potential, effective refractory period, conduction velocity (CV) restitution, and wave length restitution for common atrial myocardium (CAM) and fast conducting Crista Terminalis (CT). We used isotropic, homogeneous tissue patches. External stimuli were applied with 184 different pacing rates (PRs) from 330 to 1250 milliseconds.Results: The effect of temporary HD- related electrolyte changes on the action potential morphology and effective refractory period showed results consistent with the previous single-cell study. Action potential morphology was not significantly altered both in CAM and CT, but resting potential decreased from ␣82.6 to ␣88.2 mV for CAM and from ␣81.7 to ␣87.3 mV for CT. Effective refractory period decreased from 32 (pre-HD) to 308 milliseconds (end-HD). At a PR of 832 milliseconds, CV dropped by ␣6.3% for both types of tissue (CAM: 741 694 mm/s; CT: 746 699 mm/s). Wave length increased slightly with higher PR, but rapidly fell off below a PR of 450 milliseconds. Wave length was ␣30 mm shorter in the end-HD condition.Conclusions: Conduction velocity decrease and consequent wave length shortening increases vulnerability for atrial fibrillation onset, especially in conjunction with structural dilation often present in atria of end-stage renaldisease patients. Temporary HD-caused electrical remodeling has equal effects on regular and fast-conducting tissue. Although there is no biophysical model for fast interatrial condition pathways (eg, Bachmann dundle) available, the HD influence on them should also be similar and therefore slow down interatrial conduction significantly. It has been suggested that constantly repeating alteration of atrial electrophysiology may lead to a longer lasting electrical atrial remodeling; future studies should therefore investigate the long-term HD effects.
M. W. Krueger, F. M. Weber, G. Seemann, and O. Doessel. Personalizing Anatomical and Electrophysiological Models of the Human Atria. In Biomedizinische Technik / Biomedical Engineering (Proc. BMT 2011), vol. 56(s1) , 2011
Abstract:
IntroductionAtrial fibrillation (AF) is the most common cardiac arrhythmia. Over 4.5 million people in the European Union suffer from AF. The mechanisms leading to AF are still not completely understood although various theories were proposed. Numerical models of the human atria can help to understand these mechanisms. Personalized atrial models may in fu- ture be used to set up patient-specific therapies.MethodsPersonalization of atrial models splits into different tasks. The individual atrial and thorax anatomy are derived from various imaging modalities (CT, MRI). Valuable information is hidden in these data, such as atrial wall thickness and myocardial fiber structure. The missing parts are added to the geometric model using rule-based approaches. Atrial electrophysiology is adapted to different pathologies (e.g. remodeling, genetic defects) and to ECG and intracardiac measurements of the individual patient by tuning model parameters (e.g. conductivity).ResultsPersonalization of atrial anatomy enables a realistic simulation of atrial excitation propagation during sinus rhythm. Ad- justment of the generalized electrophysiology model to the according patient group provides insights into the substrate of the known global effects. Adaption of these model parameters to the individual patient results in a better fit of simu- lated intracardiac and ECG signals to the measurements.ConclusionWith the help of various personalization techniques, generalized atrial models can be adapted to patient data. These models may in future be used for personalized model-based AF-treatment planning.
P. Neher, H. Barschdorf, S. Dries, F. M. Weber, M. W. Krueger, O. Dössel, and C. Lorenz. Automatic segmentation of cardiac CTs - personalized atrial models augmented with electrophysiological structures. In Functional Imaging and Modeling of the Heart 2011, Lecture Notes in Computer Science, vol. 6666, pp. 80-87, 2011
Abstract:
Electrophysiological simulations of the atria could improve diagnosis and treatment of cardiac arrhythmia, like atrial fibrillation or flutter. For this purpose, a precise segmentation of both atria is needed. However, the atrial epicardium and the electrophysiological structures needed for electrophysiological simulations are barely or not at all detectable in CT-images. Therefore, a model based segmentation of only the atrial endocardium was developed as a landmark generator to facilitate the registration of a finite wall thickness model of the right and left atrial myocardium. It further incorporates atlas information about tissue structures relevant for simulation purposes like Bachmanns bundle, terminal crest, sinus node and the pectinate muscles. The correct model based segmentation of the atrial endocardium was achieved with a mean vertex to surface error of 0.53 mm for the left and 0.18 mm for the right atrium respectively. The atlas based myocardium segmentation yields physiologically correct results well suited for electrophysiological simulations.
S. Ponto, C. Schilling, M. W. Krueger, F. M. Weber, S. Seemann, and O. Dössel. Influence of endocardial catheter contact on properties of the atrial signal and comparison with simulated electrograms. In Biomedizinische Technik / Biomedical Engineering (Proc. BMT 2011), vol. 56(s1) , 2011
Abstract:
In spite of the considerable medical and technical progress during the last years, catheter ablation of atrial fibrillation is still challenging. For a successful execution of the ablation and the avoidance of intricacies the catheter must be in contact with the endocardium, which is still difficult to assure with existent techniques. It would be desirable to detect the endocardial catheter contact directly from the signal shape and its properties. In this work, significant signal property changes were detected and investigated, which allow an automatic contact detection. Furthermore, atrial electrograms were simulated and compared with a database of measured and annotated signals. During these simulations, the distance between endocardium and the catheter tip could be chosen discretionary. The simulated signals revealed themselves to be very accurate. Simulations can now be used to analyse intracardiac signals more closely. The exact position of the catheter will hereby always be assured, which is not always granted in clinical practice.
Catheter ablation of complex atrial arrhythmias is a frequently applied procedure, but its success rates are only moderate and highly dependent on the experience of the physician. Personalized atrial simulation models could assist the physician in treatment planning and thus increase success rates. In this work we created a personalized anatomical model for a specific patient from CT image data. Left atrial conduction velocity and local wave directions were determined from intracardiac electrogram (EGM) recordings. We simulated normal sinus rhythm and the clinical pacing protocol using a Cellular Automaton. The incidence direction and conduction velocity were extracted from the simulated data and compared to the results of the clinical EGMs of the same patient. We then showed that the incidence angles differed by less than 15% and that the conduction velocity error was below 12 cm/s. This implies that the model has similar electric properties compared to the real atria. In conclusion, we have presented a workflow for model personalization and validation.
Current models of the human atria represent geometries of single individuals or base on statistical data. We present a work-flow for the creation of patient-specific atrial models. Furthermore we show a framework to compare simulated P- waves and body surface potential maps (BSPMs) of individual patients with measurements. Models of the atrial and thorax anatomy were segmented from MRI data. Volumetric atrial models were semi-automatically enhanced with electrophys- iologically (EP) relevant structures. Simulations were performed on an anisotropic voxel-based mesh and were forward calculated to obtain simulated BSPMs. BSPMs were acquired using a 64 electrode ECG system. Comparison of simulated and measured P-waves in Einthoven leads showed a general agreement of both, although no personalization of the atrial electrophysiology model was performed. P-wave duration was longer in the simulations, highlighting the need for elec- trophysiological model personalization. Simulated and measured BSPMs revealed similar patterns. The presented method enables realistic simulations of atrial activation on patient-specific volumetric atrial models with EP relevant myocardial structures resulting in computed ECGs (P-wave) and BSPMs with show physiological morphologies
A framework for the automatic extraction and generation of patient-specific organ models from different image modalities is presented. These models can be used to extract and represent diagnostic information about the heart and its function. Furthermore, the models can be used for treatment planning and an overlay of the models onto X-ray fluoroscopy images can support navigation when performing an intervention in the CathLab.
Atrial fibrillation is the most common cardiac arrhythmia and often leads to severe complications such as stroke and other embolic incidents. Areas of complex fractionated atrial activity are in the focus of electrophysiologists and have been used as a target for catheter ablation therapy. The underlying mechanisms of complex fractionated atrial electrograms (CFAEs) are not entirely understood. CFAEs may contain concurrent rhythmic episodes of signals with differing characteristic frequencies (CFs). We propose a new algorithm to detect multiple periodicities in atrial signals.First, we preprocess the signal by applying Teager's non-linear energy operator. Next, the first three characteristic frequencies are detected in the frequency spectrum. Information contained in the harmonics is used to recursively detect the exact frequency. Frequency information is then transformed into the time domain, where repeated occurance of signal activity according to the respective cycle length is found. Further, the detection rate and the mean distance to gravity are calculated as key figures to determine more characteristics of the periodicity.The algorithm performs well in detecting the rhythmic components of atrial signals. It has been tested using real patient data acquired during electrophysiological studies in sinus rhythm, atrial flutter and several forms of atrial fibrillation, as well as with simulated data produced by a cellular automaton at our research group.Its application may provide new insights into atrial signals especially CFAEs and the interpretation of characteristic and dominant frequencies. It can be the foundation of displaying rhythmicity and CF information onto the 3D representation of the patient's atrium and give the physician an impression of the organization and regularity of cardiac electrograms.
S. Bauer, D. U. J. Keller, F. M. Weber, P. Tri Dung, O. Dössel, and G. Seemann. How do tissue conductivities impact on forward-calculated ECGs? An efficient prediction based on principal component analysis. In IFMBE Proceedings World Congress on Medical Physics and Biomedical Engineering, vol. 25/4, pp. 641-644, 2009
M. W. Krueger, F. M. Weber, O. Dössel, and G. Seemann. Influence of myocardial structures on electrophysiologic simulations in patient specific atrial models. In The Cardiac Physiome: Multi-scale and Multi-physics Mathematical Modelling Applied to the Heart, 2009
M. W. Krueger, F. M. Weber, O. Dössel, and G. Seemann. Semi-automatic segmentation of sinus node, Bachmann's Bundle and Terminal Crest for patient specific atrial models. In World Congress on Medical Physics and Biomedical Engineering. IFMBE Proceedings, vol. 25/4, pp. 673-676, 2009
Abstract:
The human atria contain fine structures, which can hardly be distinguished with common medical imaging techniques. However, some of these structures play an important role in the electrophysiologic depolarisation sequence of the atria. We present a semi-automatic algorithm to segment the sinus node, Bachmann’s Bundle and the Terminal Crest in given anatomical shape models of the atria. The algorithm bases on anatomical knowledge of the atria and only requires the user to provide few distinct landmarks in the atria as input. Incorporation of these structures into patient individual atrial geometries augments the electrophysiological correctness of the models.
After mathematical modeling of the healthy heart now modeling of diseases comes into the focus of research. Modeling of arrhythmias already shows a large degree of realism. This offers the chance of more detailed diagnosis and computer assisted therapy planning. Options for genetic diseases (channelopathies like Long-QT-syndrome), infarction and infarction-induced ventricular fibrillation, atrial fibrillation (AF) and cardiac resynchronization therapy are demonstrated.
M. Nalbach, J. Nenonen, O. Dössel, and O. Weber. MCG and ECG Source Reconstruction using a 4D-Model of the Human Body. In Biomedizinische Technik, vol. 46-2, pp. 57-59, 2001
Dissertations (1)
F. M. Weber. Personalizing simulations of the human atria: Intracardiac measurements, tissue conductivities, and cellular electrophysiology. KIT Scientific Publishing. Dissertation. 2011
Abstract:
This work addresses major challenges of heart model personalization. Analysis techniques for clinical intracardiac electrograms determine wave direction and conduction velocity from single beats. Electrophysiological measurements are simulated to validate the models. Uncertainties in tissue conductivities impact on simulated ECGs. A minimal model of cardiac myocytes is adapted to the atria. This makes personalized cardiac models a promising technique to improve treatment of atrial arrhythmias.