Donate Help Contact The AHA Sign In Home
American Heart Association
Circulation
Search: search_blue_button Advanced Search
Circulation. 2005;112:3769-3776
Published online before print December 5, 2005, doi: 10.1161/CIRCULATIONAHA.105.565002
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
112/24/3769    most recent
CIRCULATIONAHA.105.565002v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowRequest Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ector, J.
Right arrow Articles by Heidbüchel, H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ector, J.
Right arrow Articles by Heidbüchel, H.
Right arrowPubmed/NCBI databases
Medline Plus Health Information
*Diagnostic Imaging
Related Collections
Right arrow Electrophysiology
Right arrow Ablation/ICD/surgery
Right arrow CT and MRI
Right arrowRelated Article

(Circulation. 2005;112:3769-3776.)
© 2005 American Heart Association, Inc.


Imaging

Cardiac Three-Dimensional Magnetic Resonance Imaging and Fluoroscopy Merging

A New Approach for Electroanatomic Mapping to Assist Catheter Ablation

Joris Ector, MD; Stijn De Buck, MSc; Jef Adams, MS; Steven Dymarkowski, MD, PhD; Jan Bogaert, MD, PhD; Frederik Maes, MSc, PhD; Hein Heidbüchel, MD, PhD

From the Departments of Cardiology (J.E., J.A., H.H.), Radiology (S.D., J.B.), and Electrical Engineering (S.D.B., F.M.), University Hospital Gasthuisberg, University of Leuven, Leuven, Belgium.

Correspondence to Prof Dr Hein Heidbüchel, MD, PhD, Cardiology-Electrophysiology, University Hospital Gasthuisberg, Herestraat 49, B-3000, Leuven, Belgium. E-mail Hein.Heidbuchel{at}uz.kuleuven.ac.be

Received May 30, 2005; revision received August 24, 2005; accepted September 12, 2005.


*    Abstract
up arrowTop
*Abstract
down arrowIntroduction
down arrowMethods
down arrowResults
down arrowDiscussion
down arrowReferences
 
Background— Modern nonfluoroscopic mapping systems construct 3D electroanatomic maps by tracking intracardiac catheters. They require specialized catheters and/or dedicated hardware. We developed a new method for electroanatomic mapping by merging detailed 3D models of the endocardial cavities with fluoroscopic images without the need for specialized hardware. This developmental work focused on the right atrium because of the difficulties in visualizing its anatomic landmarks in 3D with current approaches.

Methods and Results— Cardiac MRI images were acquired in 39 patients referred for radiofrequency catheter ablation using balanced steady state free-precession sequences. We optimized acquisition and developed software for construction of detailed 3D models, after contouring of endocardial cavities with cross-checking of different imaging planes. 3D models were then merged with biplane fluoroscopic images by methods for image calibration and registration implemented in a custom software application. The feasibility and accuracy of this merging process were determined in heart-cast experiments and electroanatomic mapping in patients. Right atrial dimensions and relevant anatomic landmarks could be identified and measured in all 3D models. Cephalocaudal, posteroanterior, and lateroseptal diameters were, respectively, 65±11, 54±11, and 57±9 mm; posterior isthmus length was 26±6 mm; Eustachian valve height was 5±5 mm; and coronary sinus ostium height and width were 16±3 and 12±3 mm, respectively (n=39). The average alignment error was 0.2±0.3 mm in heart casts (n=40) and 1.9 to 2.5 mm in patient experiments (n=9), ie, acceptable for clinical use. In 11 patients, reliable catheter positioning and projection of activation times resulted in 3D electroanatomic maps with an unprecedented level of anatomic detail, which assisted ablation.

Conclusions— This new approach allows activation visualization in a highly detailed 3D anatomic environment without the need for a specialized nonfluoroscopic mapping system.


Key Words: catheter ablation • atrium • mapping • magnetic resonance imaging • imaging


*    Introduction
up arrowTop
up arrowAbstract
*Introduction
down arrowMethods
down arrowResults
down arrowDiscussion
down arrowReferences
 
Radiofrequency catheter ablation techniques have had a dramatic impact on the treatment of various forms of arrhythmias. Fluoroscopy-guided catheter navigation and ablation of accessory pathways, AV nodal slow pathways, and atrial flutter are accepted as first-line therapy with excellent outcomes.1–3 In the late 1990s, catheter ablation of more complex arrhythmias, such as postsurgical intra-atrial reentrant tachycardia, ventricular tachycardia, and atrial fibrillation, encouraged the development of new nonfluoroscopic mapping techniques in which electrical data and anatomic information are combined into 3D electroanatomic models.4–6 These systems continuously record the location of a roving mapping catheter and create a map of electrical activity in 3D space. However, they require a dedicated system and often special catheters, which makes these map- ping systems expensive. Moreover, by limiting the anatomic information to points collected by a roving catheter, the resulting model is only a rough geometric approximation of the endocardial cavity, because points that are difficult to reach are not included. The latest technological developments allow merging of these electroanatomic maps with MRI- or CT-based 3D models.7 However, this approach requires registration of two 3D models with often very differing appearances to project the electrical activation map on the anatomic model. Most research work has focused on 3D modeling of the left atrium.8–11 The right atrium is less suitable for 3D modeling because of its anatomic complexity and inhomogeneous filling with gadolinium and iodine contrast agents. Nevertheless, it is the site of most arrhythmias treated in the electrophysiology laboratory.

Editorial p 3677

We developed an approach to provide detailed, accurate anatomic localization of complex arrhythmogenic substrates that can be used in a classic biplane fluoroscopy setup without the need for a specialized nonfluoroscopic mapping system. Anatomic structures relevant to electrophysiological procedures (eg, eustachian valve, coronary sinus, terminal crest, and pulmonary veins) were visualized in a detailed 3D model of the endocardial cavity that was based on cardiac magnetic resonance images. This detailed 3D model was then directly integrated with the electrical data during the electrophysiological procedure, with biplane fluoroscopic and angiographic images used for image calibration and registration. In this report, we describe our new method of electroanatomic mapping and our initial experience in a group of 39 patients referred for catheter ablation. This developmental work concentrated on the right atrium, given the challenge it poses for detailed electroanatomic mapping, but the methodology can easily be extrapolated to the left atrium.


*    Methods
up arrowTop
up arrowAbstract
up arrowIntroduction
*Methods
down arrowResults
down arrowDiscussion
down arrowReferences
 
Image Acquisition
Cardiac magnetic resonance images were acquired on a 1.5-T Intera CV MRI unit (Philips Medical Systems) at least 1 day before the radiofrequency catheter ablation procedure (n=39). A balanced fast-field echo (bFFE) sequence with a slice thickness of 6 mm was performed in 3 orthogonal planes during different breath holds. In 11 of 39 patients, we performed an additional 3D balanced turbo-field echo (bTFE) sequence during free breathing, using prospective respiratory navigator echo gating. This more time-consuming sequence produced a high-spatial-resolution data set with a slice thickness of 1.5 mm, which could be resliced to any imaging plane during postacquisition processing. In all sequences, ECG gating was performed at end systole, to obtain image acquisition during atrial diastole for optimizing later image alignment with atrial angiography. We did not administer intravenous gadolinium contrast agent, because all peripherally injected contrast agents are inhomogeneously distributed in the right atrium and therefore do not facilitate right atrial 3D modeling (Figure 1).



View larger version (40K):
[in this window]
[in a new window]
 
Figure 1. Inhomogeneous contrast filling of the right atrium illustrated on transverse 64-slice CT image. A, Iodine contrast agent is heterogeneously distributed in the right atrium (asterisk and small arrows), which causes dropout of the posterior part of the right atrium on an intensity-based 3D model, as shown by the asterisk in B. C, 3D model of the same right atrium based on manual contouring, showing a prominent eustachian valve (EV) that was not demonstrated in B because of inhomogeneous contrast distribution. RA indicates right atrium; RAA, right atrial appendage; SVC, superior vena cava; RV, right ventricle; CS, coronary sinus; LV, left ventricle; RAO, right anterior oblique view; and LAO, left anterior oblique view.

Image Processing to Generate Detailed 3D Anatomic Models
The endocardial contours of the heart chambers were manually delineated on the MRI images with in-house developed software that allowed for cross-checking of the contours with those in other imaging planes (Figure 2). The software also compensated for positional differences in orthogonal planes that were acquired during different breath holds in the bFFE sequence, with automated 3D image registration by maximization of mutual information.12 All manually delineated contours in different 2D imaging planes were then combined into a single 3D surface model with a scattered data interpolation algorithm based on radial basis functions.13 Right atrial dimensions were measured directly on the 3D model along the cephalocaudal, posteroanterior and lateroseptal axes. In addition, posterior isthmus length, eustachian valve height, and coronary sinus ostium height and width were measured on the 3D models.



View larger version (70K):
[in this window]
[in a new window]
 
Figure 2. Cardiac 3D bTFE MRI imaging sequence. A, Manual contour of right atrium in axial imaging plane (yellow line), taking into account metal artifact (Ar) due to sternotomy clips. B, Cross-checking of contours in vertical long-axis plane (purple dots). C, Resulting 3D model of right atrium in right anterior oblique view. RA indicates right atrium; RV, right ventricle; LV, left ventricle; CS, coronary sinus; EV, eustachian valve; SVC, superior vena cava; and IVC, inferior vena cava.

Integration of 3D Model in Biplane Fluoroscopic Framework
All radiofrequency catheter ablation procedures were performed with a conventional biplane fluoroscopy system (Bicor, Siemens) and standard electrophysiological mapping and ablation catheters. A custom software application was developed to allow integration of the 3D anatomic models into the biplane fluoroscopy framework. This integration process consisted of 2 steps:

  1. Fluoroscopy geometry calibration: To render the 3D model in the biplane fluoroscopic framework in its proper aspect ratio, a calibration procedure was performed based on the fluoroscopic view angles and a calibration object. The fluoroscopic view angles were determined automatically from the DICOM (Digital Imaging and Communications in Medicine) header information of the fluoroscopic images. As a calibration object, we used intracardiac catheters with known interelectrode distances, which were measured on fluoroscopic images in both viewing planes (a minimum of 7 electrode positions were marked in each plane).
  2. Image registration: After the calibration procedure, the 3D model could be introduced in the fluoroscopic framework in its proper aspect ratio at a standard offset. Initial model rotation was calculated from the fluoroscopic view angles. Subsequently, the transparent model was carefully manually translated and rotated by means of a 3D mouse (SpaceMouse, 3Dconnexion Inc) to properly align it with angiographic images of the heart in 2 projections (at end systole) by a modified version of the visual matching technique (Figure 3).14 The acquisition of angiographic images was performed with the same end-systolic triggering interval as for the MRI acquisition to avoid synchronicity mismatch.



View larger version (138K):
[in this window]
[in a new window]
 
   Figure 3. Integration of 3D model of right atrium and right ventricle in fluoroscopy framework. Alignment of the model to the angiography image is performed by visual matching. RA indicates right atrium; RV, right ventricle; CS, coronary sinus; RAO, right anterior oblique view; LAO, left anterior oblique view; SVC: superior vena cava; and IVC, inferior vena cava.

The feasibility and accuracy of this integration process were assessed in 2 ways:

  1. Phantom experiments: Registration by the proposed visual matching technique was tested on a set of simulated fluoroscopy images of a plastic endocardial negative heart cast. These fluoroscopic images were rendered in a 3D environment together with the 3D model of the cast, which was positioned at a random translational and rotational offset. Subsequently, the 3D model was visually aligned with the simulated fluoroscopic images from 20 different offsets by 2 independent operators (1 cardiologist familiar with cardiac anatomy and 1 clinical engineer). Alignment error in the x, y, and z axes, angular deviation, and the time to complete the task were recorded.
  2. Patient validation: After registration of right atrial 3D models with right atrial angiography in 9 patients, we measured the distance between heart boundaries visible on angiography and projected boundaries of the 3D model rendered in the fluoroscopy framework. Distance between boundaries was defined as the area between the 2 manually outlined boundaries divided by boundary length (Figure 4). This alignment error was evaluated separately for posterior, inferior, and anterior borders in the right anterior oblique view and for lateral and septal borders in the left anterior oblique view. If no clear border could be visualized on angiography (eg, due to incomplete filling with contrast agent), this region was not used for measurement.



View larger version (73K):
[in this window]
[in a new window]
 
   Figure 4. Right anterior oblique view of right atrium. A, Clearly visible right atrial boundaries are marked on angiography image (black line) and measured to determine boundary length. B, 3D model integrated in fluoroscopy framework. C, The area between boundaries marked on angiography and the boundaries of the 3D model is labeled in red. Mean difference in distance for each boundary is calculated as area divided by boundary length. SVC indicates superior vena cava; IVC, inferior vena cava; and RA, right atrium.

Projection of Catheter Positions From Fluoroscopic Images to 3D Model
After proper image calibration and registration, the position of intracardiac catheters could be projected from the fluoroscopy images to the 3D model. For this purpose, 3 anchor points of the catheter were manually marked on the fluoroscopy images and automatically projected onto the 3D model as a virtual catheter. Further modifications of the catheter curve were possible in the 3D environment using 3D crosshairs. The position of the virtual coronary sinus catheter within the coronary sinus ostium of the 3D model was used for additional evaluation of correct fluoroscopy model alignment in 11 patients (Figure 5).



View larger version (146K):
[in this window]
[in a new window]
 
Figure 5. A, Projection of coronary sinus (CS) catheter from fluoroscopic image to 3D model (left anterior oblique view). B, Resulting CS catheter position in relation to 4-chamber 3D model of same patient. C, Endocardial view of right atrial 3D model in same patient, with virtual CS catheter visible in CS ostium. RV indicates right ventricle; LV, left ventricle; Ao, aorta; LA, left atrium; CS, coronary sinus catheter; IVC, inferior vena cava; ER, eustachian ridge; Isthmus, posterior isthmus; and TA, tricuspid annulus.

Electroanatomic Mapping on 3D Model
During conventional electrophysiological mapping of arrhythmias, multiple positions from different catheters were projected onto the 3D model by the method described above. Their electrograms were recorded, and bipolar activation times on different electrode pairs (2 mm spacing) were graphically represented on the 3D surface, which resulted in a color-coded electrical activation map of the arrhythmia (Figure 6). During this developmental work, activation times were manually marked on the catheter electrograms. Fractionated and double potentials could also be represented on the 3D model, as well as scar tissue and orifices. Multielectrode mapping catheters were used for rapid electrical data acquisition at different catheter positions to accelerate the mapping process.15



View larger version (86K):
[in this window]
[in a new window]
 
Figure 6. Color-coded activation map of right atrium in right anterior oblique view, showing a counterclockwise reentrant tachycardia around an old atriotomy scar (SCAR) in a patient with previous mitral valve surgery. Left panels show the fluoroscopic framework used for creation of the activation map. Yellow dots mark recorded activation times. SVC indicates superior vena cava; IVC, inferior vena cava; and TC, terminal crest.

The investigational review board approved this study, and informed consent was obtained from all patients.

Statistical Analysis
Data are presented as mean±SD unless otherwise stated.


*    Results
up arrowTop
up arrowAbstract
up arrowIntroduction
up arrowMethods
*Results
down arrowDiscussion
down arrowReferences
 
Image Acquisition and Processing to 3D Model
Imaging time required for acquisition of the 3D bTFE sequence exceeded imaging time for the bFFE sequence by more than 10-fold (15±2.5 versus 1±0.2 minutes). Time required for manual contouring of the right atrium was not different between the 2 sequences (bFFE 1.5±0.4 hours, bTFE 1.6±0.4 hours, P=NS), even though more slices had to be contoured in the bTFE sequence (bTFE 120±5 axial slices, bFFE 20 slices in the axial, sagittal, and coronal planes). This was explained by increased complexity of contouring and cross-checking contours in different imaging planes in the 3 orthogonal bFFE acquisitions.

The accuracy of distance measurements was evaluated on an MRI-based 3D model of a static plastic endocardial negative heart cast. Distance measurements of structures on the 3D model corresponded exactly to distances on the 3D model (differences <0.5 mm). The endocavitary volume calculated from the 3D model was 126.9 cm3, which corresponded well to the 127.5 cm3 observed by water immersion of the phantom.

The inclusion of all 3 imaging planes for contouring in the bFFE sequence allowed for identification of fine anatomic structures, such as the eustachian valve and coronary sinus, with the same accuracy as the 3D bTFE sequence. This was assessed by measuring the eustachian valve height in its middle portion and coronary sinus ostium height and width on 3D models of the same patients based on the bFFE and 3D bTFE sequences (n=11). The mean eustachian valve height was 8±5 mm for both sequences, with a mean measurement difference of 0±1 mm. The coronary sinus ostium height and width were 16±3 and 12±3 mm, respectively, in the bFFE sequence versus 15±4 and 10±3 mm in the 3D bTFE sequence, with a mean measurement difference of 0±2 mm for coronary sinus ostium height and 2±2 mm for coronary sinus ostium width. Right atrial dimensions, measured along different axes (n=39, bFFE sequence), were 65±11 mm cephalocaudal, 54±11 mm posteroanterior, and 57±9 mm lateroseptal. Mean posterior isthmus length was 26±6 mm, coronary sinus ostium height and width were 16±3 and 12±3 mm, respectively, and eustachian valve height was 5±5 mm (median 5 mm, interquartile range 0 to 9 mm).

Integration of 3D Model in Biplane Fluoroscopy Framework
The registration procedure in the heart-cast experiment could be completed within 5.0 minutes in 95% of cases by the cardiologist and in 70% of cases by the engineer. The initial random alignment errors were 0±15 mm translational and 0±5° rotational for the x, y, and z axes separately. After the registration procedure, the mean translational alignment error in the x, y, and z axes (ie, horizontal and vertical axes in the right anterior oblique view and the horizontal axis in the left anterior oblique view, respectively) for 20 registration procedures of the heart-cast model was 0.0±0.2, 0.4±0.3, and 0.2±0.3 mm, respectively, for the cardiologist (overall mean 0.2±0.3 mm) and 0.3±0.3, 0.2±0.3, and 0.4±0.3 mm, respectively, for the engineer (overall mean 0.2±0.4, difference=NS; Figure 7A). The mean angular deviations around the same x, y, and z axes were –1.0±1.6°, –0.7±0.6°, and 0.1±1.3°, respectively, for the cardiologist and –0.3±1.2°, –0.4±0.6°, and –0.3±0.8° for the engineer (difference=NS; Figure 7B).



View larger version (20K):
[in this window]
[in a new window]
 
Figure 7. A and B, Translational and rotational error as measured in 40 visual matching registration procedures of the plastic heart cast and simulated fluoroscopy images. X represents horizontal axis in right anterior oblique view; Y, vertical axis in right anterior oblique view; Z, horizontal axis in left anterior oblique view. C, Distance between borders visualized on angiography and 3D model in 9 patients.

The median distances between angiographic and 3D model boundaries in 9 patients were 1.9 mm (interquartile range 1.0 to 2.7 mm) for posterior, 2.4 mm (1.9 to 3.3 mm) for inferior, and 2.5 mm (1.5 to 9.9 mm) for anterior borders in the right anterior oblique view. In the left anterior oblique view, there was a median distance of 2.1 mm (interquartile range 1.7 to 3.1 mm) between lateral boundaries and 2.1 mm (2 to 5.1 mm) between septal boundaries (Figure 7C). When this methodology was used to calculate alignment errors between the boundaries of the MRI-based 3D model of the heart cast and its static fluoroscopy projection, alignment errors for all borders were negligible (0.37 to 0.6 mm). This suggests that most of the observed alignment errors in patients are due to cardiac motion or changes of the endocardial cavity dimension between imaging and the ablation procedure.

Coronary sinus catheters correctly followed the course of the coronary sinus as visualized on the 3D models in all 11 patients. Their most proximal segment projected correctly within the coronary sinus ostium in 7 of 11 patients and within 5 mm of the ostium in all 11 patients. This small error at the level of the coronary sinus ostium is most probably due to catheter-related pressure exerted in this region (see below).

Projection of Catheter Positions and Electroanatomic Mapping
The construction of a color-coded activation map during tachycardia or localized pacing was possible in all patients in whom it was attempted (n=11). They included patients with clockwise or counterclockwise atrial flutter (n=4), accessory pathways (n=2), AV nodal reentrant tachycardia (n=2), and complex arrhythmogenic substrates such as periatriotomy or congenital reentrant tachycardia (n=3) (Figure 6). The average time needed for completion of these activation maps (based on current manual annotation of activation times) was 30±6 minutes for an average of 62±25 mapped points. There were no significant differences in total procedural times, fluoroscopy times, and number of radiofrequency applications compared with patients who were treated conventionally in the same time period for the same type of arrhythmias (Table). The mean number of radiofrequency applications in the group with complex arrhythmogenic substrates was higher (although not statistically significant) in the electroanatomic mapping group owing to 1 patient requiring a total of 17 radiofrequency applications.


View this table:
[in this window]
[in a new window]
 
Procedure Characteristics


*    Discussion
up arrowTop
up arrowAbstract
up arrowIntroduction
up arrowMethods
up arrowResults
*Discussion
down arrowReferences
 
The main results of our research are the development of a mapping system that (1) is primarily based on highly detailed anatomic information, (2) does not require special mapping catheters, and (3) can be applied in a regular biplane fluoroscopy setup. Moreover, by using cardiac MRI for image acquisition, we did not invoke extra radiation for the patients (although our methodology could also be applied to cardiac CT).

Although MRI- or CT-based 3D volume rendering of the left atrium and pulmonary veins is frequently performed to assess left atrial anatomy (such as for catheter ablation of atrial fibrillation9,10), a similar approach for the right atrium has not yet been developed. Intravenously administered contrast agents are not homogeneously distributed in the right atrium (in contrast to the left atrium) and therefore do not allow easy segmentation and 3D visualization of the right atrium. Nevertheless, the majority of radiofrequency catheter ablation procedures are right-sided, and exact knowledge of right atrial anatomy could certainly benefit mapping procedures in patients with complex arrhythmogenic substrates, such as those with congenital heart disease or previous cardiac surgery. Therefore, we developed a methodology to construct right atrial 3D models without contrast agents by manually delineating and cross-checking endocardial contours in different imaging planes. This resulted in surface-rendered 3D models of the right atrium with an unprecedented level of detail, which also represented fine anatomic structures relevant to electrophysiology procedures. The time required for developing such a detailed 3D model of the right atrium is still substantial (1.5 hours), but it may be reduced by automating certain steps in the delineation process with model-based segmentation approaches.16 The preprocedural availability of a detailed and anatomically correct model may also allow preprocedural planning and selection of specialized guiding sheaths and ablation catheters. This is an important advantage compared with existing nonfluoroscopic mapping systems.

The methodology of manual contouring in differently oriented imaging planes allowed the use of the bFFE MRI sequence with 6-mm slice thickness to construct 3D models without losing the ability to include fine anatomic structures such as the eustachian valve and coronary sinus. The possibility of using the bFFE sequence is an advantage, because the image acquisition time is substantially lower than for the 3D bTFE sequence, and the time needed for manual delineation is comparable. The accuracy of MRI, including these steady state free-precession sequences, with regard to quantitative assessment of cardiac dimensions and function has been validated in many reports.17–20 The present study data confirm the possibility of using these MRI sequences for visualization of detailed atrial structures and their relation to surrounding anatomy. This level of anatomic detail cannot be achieved when using a roving catheter for construction of 3D geometry, even when efforts are made to map a sufficient number of points.21 This limitation encouraged the new technologies for merging electroanatomic maps with CT/MRI-based 3D models. As these merging approaches are further developed, a more direct comparison between MRI/CT-based and roving catheter-based models will be available.

The integration of detailed cardiac 3D models in the biplane fluoroscopic framework is a new concept. Its main advantage is the wide availability and general use of biplane fluoroscopy in electrophysiology laboratories around the world. The registration of the 3D models with angiography images was performed by a visual matching technique. Theoretically, automated registration of magnetic resonance images with fluoroscopy is possible with intensity-based registration measures,22 but we found these methods not to be reliable enough for use in clinical practice. The visual matching procedure could be performed with high accuracy and reproducibility by 2 independent operators. The maximal translational alignment errors in the heart-cast experiments were <1 mm, and the median distance between angiographic and model borders in patients ranged from 1.9 to 2.5 mm. The large error for the anterior border (median 2.5 mm, interquartile range 1.5 to 9.9 mm) is explained by the fact that this border represents the tricuspid valve in the right anterior oblique view. Very small differences in timing of the angiography image acquisition can have a profound impact on valve position and are the most likely explanation for the observed difference. The observed alignment error is certainly acceptable for use in clinical practice. This was confirmed by the finding of correct catheter projection with respect to anatomic landmarks (coronary sinus ostium experiments) and consistent electroanatomic maps of complex arrhythmias.

To further diminish alignment errors of different borders and enhance the anatomic accuracy of the mapping process, we are investigating a periprocedural correction of the 3D model to angiographic images. Such an approach for periprocedural fine-tuning of 3D models is desirable because changes in atrial dimensions might occur owing to (1) differences in heart rhythm and volume-loading conditions between MRI imaging and ablation procedure and (2) pressure exerted from sometimes highly curved catheters on the endocardial wall. The projection of the coronary sinus catheter (<5 mm) outside the coronary sinus ostium of the 3D model in 4 of 11 patients might have been caused in part by the latter effect. The quantification of these changes and the need for periprocedural correction in different clinical settings will be the subject of future investigations. Acquisition of MRI images and angiographic images was triggered to the same end-systolic interval to avoid synchronicity mismatch. This not only allows imaging during atrial diastole but also causes less error when changes in heart rate occur between MRI acquisition and the ablation procedure, because systolic duration is less influenced by heart rate than diastolic duration.

The current delay in completion of the activation map (30±6 minutes) did not result in a significant increase in total procedure times, although more patients are needed for further evaluation of procedure and fluoroscopy times, number of radiofrequency applications required, and clinical outcome within different patient groups. Moreover, by automating the marking of catheter positions on the fluoroscopic images and the determination of activation times, the time required to complete the activation map will be shortened so that this electroanatomic approach can be used efficiently in clinical practice without the need for specialized catheters or hardware. The lower cost of this software and fluoroscopy-based approach may bring the benefits of electroanatomic mapping to a broader array of electrophysiology laboratories and to more patients.

Limitations
In its present form, our new concept for electroanatomic mapping has some shortcomings compared with available electroanatomic mapping systems. Currently, real-time tracking of the position of the ablation/mapping catheter is not possible. However, the ability to indicate an electrode in biplane fluoroscopic images in both the right and left anterior oblique views by simply clicking on it allows for its localization on the 3D model (in an epicardial or endocardial view), a maneuver that only requires a few seconds. This is a limitation compared with other merge technologies, but we think that it might suffice in most procedures. In the future, we intend to combine real-time catheter tracking within the 3D angiography framework.

MRI acquisition was performed at least 1 day before the ablation procedure, with possibly different heart rate, rhythm, and volume-loading conditions. The effect of these factors on model accuracy must be determined in future work.

The accuracy of the MRI-based 3D models was assessed by calculating the alignment error between heart boundaries of the 3D model and angiography images. A more direct evaluation of model accuracy is not possible. In our opinion, this is the most realistic method to evaluate model accuracy in vivo, despite the absence of a "gold standard." Measurements on a phantom, although not subject to cardiac motion, showed the high accuracy of our 3D angiography merge process.


*    Acknowledgments
 
Dr Ector is a Research Assistant and Dr Heidbüchel is a Fundamental Clinical Investigator of the Fund for Scientific Research-Flanders. Dr Heidbüchel is holder of the AstraZeneca Chair in Cardiac Electrophysiology, University of Leuven.

Disclosures

Dr Heidbüchel is supported by an unconditional research grant from Medtronic, Inc. He is a member of the advisory board of Biosense Webster, Inc, and St. Jude Medical, Inc. The remaining authors report no conflicts of interest.


*    References
up arrowTop
up arrowAbstract
up arrowIntroduction
up arrowMethods
up arrowResults
up arrowDiscussion
*References
 
1. Jackman WM, Wang XZ, Friday KJ, Roman CA, Moulton KP, Beckman KJ, McClelland JH, Twidale N, Hazlitt HA, Prior MI. Catheter ablation of accessory atrioventricular pathways (Wolff-Parkinson-White syndrome) by radiofrequency current. N Engl J Med. 1991; 324: 1605–1611.[Abstract]

2. Haissaguerre M, Warin JF, Lemetayer P, Saoudi N, Guillem JP, Blanchot P. Closed-chest ablation of retrograde conduction in patients with atrioventricular nodal reentrant tachycardia. N Engl J Med. 1989; 320: 426–433.[Abstract]

3. Heidbuchel H, Willems R, van RH, Adams J, Ector H, Van de WF. Right atrial angiographic evaluation of the posterior isthmus: relevance for ablation of typical atrial flutter. Circulation. 2000; 101: 2178–2184.[Abstract/Free Full Text]

4. Ben-Haim SA, Osadchy D, Schuster I, Gepstein L, Hayam G, Josephson ME. Nonfluoroscopic, in vivo navigation and mapping technology. Nat Med. 1996; 2: 1393–1395.[CrossRef][Medline] [Order article via Infotrieve]

5. Smeets JL, Ben-Haim SA, Rodriguez LM, Timmermans C, Wellens HJ. New method for nonfluoroscopic endocardial mapping in humans: accuracy assessment and first clinical results. Circulation. 1998; 97: 2426–2432.[Abstract/Free Full Text]

6. Gornick CC, Adler SW, Pederson B, Hauck J, Budd J, Schweitzer J. Validation of a new noncontact catheter system for electroanatomic mapping of left ventricular endocardium. Circulation. 1999; 99: 829–835.[Abstract/Free Full Text]

7. Verma A, Marrouche N, Natale A. Novel method to integrate three-dimensional computed tomographic images of the left atrium with real-time electroanatomic mapping. J Cardiovasc Electrophysiol. 2004; 15: 968.[Medline] [Order article via Infotrieve]

8. Cronin P, Sneider MB, Kazerooni EA, Kelly AM, Scharf C, Oral H, Morady F. MDCT of the left atrium and pulmonary veins in planning radiofrequency ablation for atrial fibrillation: a how-to guide. Am J Roentgenol. 2004; 183: 767–778.[Free Full Text]

9. Wittkampf FHM, Vonken EJ, Derksen R, Loh P, Velthuis B, Wever EFD, Boersma LVA, Rensing BJ, Cramer MJ. Pulmonary vein ostium geometry: analysis by magnetic resonance angiography. Circulation. 2003; 107: 21–23.[Abstract/Free Full Text]

10. Kato R, Lickfett L, Meininger G, Dickfeld T, Wu R, Juang G, Angkeow P, LaCorte J, Bluemke D, Berger R, Halperin HR, Calkins H. Pulmonary vein anatomy in patients undergoing catheter ablation of atrial fibrillation: lessons learned by use of magnetic resonance imaging. Circulation. 2003; 107: 2004–2010.[Abstract/Free Full Text]

11. Lin WS, Prakash VS, Tai CT, Hsieh MH, Tsai CF, Yu WC, Lin YK, Ding YA, Chang MS, Chen SA. Pulmonary vein morphology in patients with paroxysmal atrial fibrillation initiated by ectopic beats originating from the pulmonary veins: implications for catheter ablation. Circulation. 2000; 101: 1274–1281.[Abstract/Free Full Text]

12. Maes F, Collignon A, Vandermeulen D, Marchal G, Suetens P. Multimodality image registration by maximization of mutual information. IEEE Trans Med Imaging. 1997; 16: 187–198.[CrossRef][Medline] [Order article via Infotrieve]

13. Carr JC, Fright WR, Beatson RK. Surface interpolation with radial basis functions for medical imaging. IEEE Trans Med Imaging. 1997; 16: 96–107.[CrossRef][Medline] [Order article via Infotrieve]

14. De Buck S, Van Cleynenbreugel J, Marchal G, Suetens P. Evaluation of a visual matching registration technique for intraoperative surgery support. In: Mun SK, ed. Medical Imaging, Visualization, Image-Guided Procedures, and Display. Bellingham, Wash: SPIE; 2002: 492–499.

15. Anne W, van RH, Adams J, Ector H, Van de WF, Heidbuchel H. Ablation of post-surgical intra-atrial reentrant tachycardia: predilection target sites and mapping approach. Eur Heart J. 2002; 23: 1609–1616.[Abstract/Free Full Text]

16. Frangi AF, Rueckert D, Schnabel JA, Niessen WJ. Automatic construction of multiple-object three-dimensional statistical shape models: application to cardiac modeling. IEEE Trans Med Imaging. 2002; 21: 1151–1166.[CrossRef][Medline] [Order article via Infotrieve]

17. Miquel ME, Hill DLG, Baker EJ, Qureshi SA, Simon RDB, Keevil SF, Razavi RS. Three- and four-dimensional reconstruction of intra-cardiac anatomy from two-dimensional magnetic resonance images. Int J Cardiovasc Imaging. 2003; 19: 239–254.[Medline] [Order article via Infotrieve]

18. Sorensen TS, Korperich H, Greil GF, Eichhorn J, Barth P, Meyer H, Pedersen EM, Beerbaum P. Operator-independent isotropic three-dimensional magnetic resonance imaging for morphology in congenital heart disease: a validation study. Circulation. 2004; 110: 163–169.[Abstract/Free Full Text]

19. Shors SM, Fung CW, Francois CJ, Finn JP, Fieno DS. Accurate quantification of right ventricular mass at MR imaging by using cine true fast imaging with steady-state precession: study in dogs. Radiology. 2004; 230: 383–388.[Abstract/Free Full Text]

20. Beygui F, Furber A, Délépine S, Helft G, Metzger JP, Geslin P, Jeune JJL. Routine breath-hold gradient echo MRI-derived right ventricular mass, volumes and function: accuracy, reproducibility and coherence study. Int J Cardiovasc Imaging. 2004; 20: 509–516.[CrossRef][Medline] [Order article via Infotrieve]

21. Kubo R, Shoda M, Fuda Y, Sugiura R, Kasanuki H. Anatomical structure of the isthmus between the inferior vena cava and tricuspid annulus investigated with a three-dimensional electroanatomical mapping system. Heart Vessels. 2005; 20: 50–55.[Medline] [Order article via Infotrieve]

22. Rohlfing T, Maurer CR Jr. A novel image similarity measure for registration of 3-D MR images and x-ray projection images. In: Dohi T, Kikinis R, eds. In: Medical Image Computing and Computer-Assisted Intervention: MICCAI 2002: 5th International Conference, Tokyo, Japan, September 25-28, 2002, Proceedings, Part II. New York, NY: Springer-Verlag; 2002:469–476.


Related Article:

Issue Highlights
Circulation 2005 112: 3667. [Extract] [Full Text]



This article has been cited by other articles:


Home page
Eur Heart JHome page
B. A. Hoffmann, A. Koops, T. Rostock, K. Mullerleile, D. Steven, R. Karst, M. U. Steinke, I. Drewitz, G. Lund, S. Koops, et al.
Interactive real-time mapping and catheter ablation of the cavotricuspid isthmus guided by magnetic resonance imaging in a porcine model
Eur. Heart J., November 5, 2009; (2009) ehp460v1.
[Abstract] [Full Text] [PDF]


Home page
EuropaceHome page
J. Sra
Cardiac registration: going further than atrial fibrillation ablation
Europace, November 1, 2009; 11(11): 1417 - 1418.
[Full Text] [PDF]


Home page
EuropaceHome page
P. Kirchhof, J. Bax, C. Blomstrom-Lundquist, H. Calkins, A. J. Camm, R. Cappato, F. Cosio, H. Crijns, H.-C. Diener, A. Goette, et al.
Early and comprehensive management of atrial fibrillation: Proceedings from the 2nd AFNET/EHRA consensus conference on atrial fibrillation entitled 'research perspectives in atrial fibrillation'
Europace, July 1, 2009; 11(7): 860 - 885.
[Full Text] [PDF]


Home page
EuropaceHome page
J.-A. Cabrera, S. Y. Ho, and D. Sanchez-Quintana
How anatomy can guide ablation in isthmic atrial flutter
Europace, January 1, 2009; 11(1): 4 - 6.
[Full Text] [PDF]


Home page
EuropaceHome page
P. Kirchhof, M. Ozgun, S. Zellerhoff, G. Monnig, L. Eckardt, K. Wasmer, W. Heindel, G. Breithardt, and D. Maintz
Diastolic isthmus length and 'vertical' isthmus angulation identify patients with difficult catheter ablation of typical atrial flutter: a pre-procedural MRI study
Europace, January 1, 2009; 11(1): 42 - 47.
[Abstract] [Full Text] [PDF]


Home page
EuropaceHome page
L. F. Tops, M. J. Schalij, D. W. den Uijl, T. P. Abraham, H. Calkins, and J. J. Bax
Image integration in catheter ablation of atrial fibrillation
Europace, November 1, 2008; 10(suppl_3): iii48 - iii56.
[Abstract] [Full Text] [PDF]


Home page
CirculationHome page
S. R. Dukkipati, R. Mallozzi, E. J. Schmidt, G. Holmvang, A. d'Avila, R. Guhde, R. D. Darrow, G. Slavin, M. Fung, Z. Malchano, et al.
Electroanatomic Mapping of the Left Ventricle in a Porcine Model of Chronic Myocardial Infarction With Magnetic Resonance-Based Catheter Tracking
Circulation, August 19, 2008; 118(8): 853 - 862.
[Abstract] [Full Text] [PDF]


Home page
EuropaceHome page
S. Knecht, H. Skali, M. D. O'Neill, M. Wright, S. Matsuo, G. M. Chaudhry, C. I. Haffajee, I. Nault, G. H.M. Gijsbers, F. Sacher, et al.
Computed tomography-fluoroscopy overlay evaluation during catheter ablation of left atrial arrhythmia
Europace, August 1, 2008; 10(8): 931 - 938.
[Abstract] [Full Text] [PDF]


Home page
RadiologyHome page
F. Saremi, L. Pourzand, S. Krishnan, O. Ashikyan, S. V. Gurudevan, J. Narula, K. Kaushal, and A. Raney
Right Atrial Cavotricuspid Isthmus: Anatomic Characterization with Multi-Detector Row CT
Radiology, June 1, 2008; 247(3): 658 - 668.
[Abstract] [Full Text] [PDF]


Home page
RadioGraphicsHome page
J. Rademaker, F. Saremi, and S. Krishnan
Invited Commentary * Authors' Response
RadioGraphics, November 1, 2007; 27(6): 1566 - 1567.
[Full Text] [PDF]


Home page
CirculationHome page
J. Ector, D. Loeckx, W. Coudijzer, S. De Buck, F. Maes, S. Dymarkowski, J. Bogaert, and H. Heidbuchel
Changes in Left Atrial and Pulmonary Venous Anatomy During Respiration: A 4-Dimensional Computed Tomography Based Assessment and Implications for Atrial Fibrillation Ablation
Circulation, June 12, 2007; 115(23): e617 - e619.
[Full Text] [PDF]


Home page
CirculationHome page
R. de Silva, L. F. Gutierrez, A. N. Raval, E. R. McVeigh, C. Ozturk, and R. J. Lederman
X-Ray Fused With Magnetic Resonance Imaging (XFM) to Target Endomyocardial Injections: Validation in a Swine Model of Myocardial Infarction
Circulation, November 28, 2006; 114(22): 2342 - 2350.
[Abstract] [Full Text] [PDF]


Home page
CirculationHome page
J. Triedman
Virtual Reality in Interventional Electrophysiology
Circulation, December 13, 2005; 112(24): 3677 - 3679.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
112/24/3769    most recent
CIRCULATIONAHA.105.565002v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowRequest Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ector, J.
Right arrow Articles by Heidbüchel, H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ector, J.
Right arrow Articles by Heidbüchel, H.
Right arrowPubmed/NCBI databases
Medline Plus Health Information
*Diagnostic Imaging
Related Collections
Right arrow Electrophysiology
Right arrow Ablation/ICD/surgery
Right arrow CT and MRI
Right arrowRelated Article