Home Centerline Tracking for Quantification of Reverse Structural Remodeling of the Pulmonary Veins Following Cardiac Ablation Therapy
Post
Cancel

Centerline Tracking for Quantification of Reverse Structural Remodeling of the Pulmonary Veins Following Cardiac Ablation Therapy

Rationale and Objectives

Patients with atrial fibrillation undergo structural remodeling resulting in increased pulmonary vein sizes. Studies have demonstrated that these changes are reversible following successful ablation therapy. To date, analyses of pulmonary vein structure have focused on measurements at the pulmonary vein ostia, and the full extent of reverse remodeling along the length of the pulmonary veins has not yet been fully characterized.

Materials and Methods

An automated, three-dimensional method is proposed that quantifies pulmonary vein geometry starting at the ostia and extending several centimeters into the veins. A centerline is tracked along the length of the pulmonary vein, and orthogonal planes are computed along the curve. The method was validated against manual measurements on each of the four pulmonary veins for 10 subjects. The proposed methodology was used to analyze the pulmonary veins in 21 patients undergoing cardiac ablation therapy with preoperative and postoperative computed tomographic scans.

Results

Validation results demonstrated that the automated measurements closely followed the manual measurements, with an overall mean difference of 11.50 mm 2 . Significant differences in cross-sectional area at the two time points were observed at all pulmonary vein ostia and extending for 2.0 cm (excluding the 0.5-cm interval) into the left inferior pulmonary vein, 3.5 cm into the left superior pulmonary vein, and 2.0 cm into the right superior pulmonary vein.

Conclusions

Quantitative analysis along the length of the pulmonary veins can be accomplished using centerline tracking and measurements from orthogonal planes along the curve. The patient study demonstrated that reverse structural remodeling following ablation therapy occurs not only at the ostia but for several centimeters extending into the pulmonary veins.

Atrial fibrillation is a condition of the heart in which the atria beat rapidly and irregularly, resulting in a lack of synchrony in the beating of the heart. Although the exact etiology of the disease is not yet fully understood, it has been shown that the origin of ectopic beats is frequently located in the pulmonary veins . In catheter ablation therapy, a standard treatment strategy, catheters are guided into the left atrium, and radiofrequency energy is delivered circumferentially to each of the pulmonary veins to create an electrical block to the left atrium. Structural remodeling occurs in conjunction with this disease, with patients with atrial fibrillation having larger left atria and pulmonary veins than normal controls. Following ablation therapy, these structural changes are reversible, a process termed reverse remodeling, with left atrial size and pulmonary vein diameters significantly decreasing in patients who have returned to normal sinus rhythm.

Left atrial size has been quantified using either two-dimensional diameter measurements or three-dimensional (3D) volume measurements from high-resolution computed tomographic and magnetic resonance data sets. Efforts have also been made to quantify pulmonary vein structure , but these techniques primarily include manual line measurements on either volume renderings or within image cross-sections. The number and location of measurements are typically limited to the junction between the pulmonary vein and left atrium and a collection of specified distances within the pulmonary veins. Accurately defining a distance into the pulmonary vein, however, can be difficult because the pulmonary veins traverse a 3D path through the volume. The current lack of 3D quantification tools for pulmonary vein structure limits the ability to assess complex changes that may occur in response to ablation therapy.

Get Radiology Tree app to read full this article<

Open full size image

Figure 1

Back view of left atrium illustrating measurements made along pulmonary vein using successive planes orthogonal to a centerline curve. LIPV, left inferior pulmonary vein; LSPV, left superior pulmonary vein; RIPV, right inferior pulmonary vein; RSPV, right superior pulmonary vein.

Get Radiology Tree app to read full this article<

Materials and methods

Data Preprocessing

Get Radiology Tree app to read full this article<

Figure 2, Three-dimensional software for segmentation of left atrium and pulmonary veins. The pulmonary veins are separated from the left atrium using an interactive cut plane, which is visualized in the three-dimensional display at the top left , and the location of the oblique cut plane is shown in three orthogonal planes along the bottom . The oblique image slice is also displayed at the top right .

Get Radiology Tree app to read full this article<

Overall Method

Get Radiology Tree app to read full this article<

Figure 3, Flowchart demonstrating the overall method for using a centerline curve to compute the cross-sectional area along a pulmonary vein (PV). LA, left atrium.

Get Radiology Tree app to read full this article<

Centerline Calculation

Get Radiology Tree app to read full this article<

Figure 4, Centerlines tracked for (a) left inferior, (b) left superior, (c) right inferior, and (d) right superior pulmonary veins.

Get Radiology Tree app to read full this article<

Cross-sectional Area Computation

Get Radiology Tree app to read full this article<

Figure 5, (a) The centerline curve with orthogonal plane near the ostium. (b) Plane accurately cutting across pulmonary vein (PV) near ostium. LA, left atrium.

Figure 6, In (a) , the orthogonal plane to the centerline curve near the pulmonary vein ostium is shown with (b) its associated image slice. In (c) , the manually adjusted cut plane from the segmentation procedure is shown with (d) its associated image slice.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

v3=α1v1+α2v2, v

3

=

α

1

v

1

+

α

2

v

2

,

where α 1 + α 2 = 1 and α 1 = 1/ c__n . c is a constant, and n is the point number along the curve. c was empirically set using a single data set and held constant for all subsequent data sets. v 3 is the unit vector used to compute the oblique image plane, v 1 is the normal of the segmentation cut plane, and v 2 is defined in equation 2 :

v2=xn−xc|xn+xc|, v

2

=

x

n

x

c

|

x

n

+

x

c

|

,

where x c is the current point and x n is the next point on the centerline. The effect of equations 1 and 2 is that near the pulmonary vein ostia, v 3 is approximately equal to the segmentation plane normal, and further into the vein, v 3 is approximately equal to the vector along the centerline. At each point along the centerline curve, the oblique image plane defined by v 3 is extracted from the volume, and the cross-sectional area for the specified pulmonary vein is computed as

A=(∑iυi)×d2s.t.υi∈PVjwherej={1,…4}, A

=

(

i

υ

i

)

×

d

2

s

.t

i

PV

j

where

j

=

{

1

,

4

}

,

where A is cross-sectional area, υ i are the voxels within the specified pulmonary vein PV j , and d is the voxel resolution.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Figure 7, Cross-sectional area versus distance computed using segmentation cut plane ( blue ), orthogonal to centerline curve ( green ), and orthogonal to combined vector ( magenta ).

Figure 8, Planes orthogonal to the orientation vector at centerline points along the right superior pulmonary vein ( top ) and corresponding oblique images extracted from volume ( bottom ).

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Validation Experiments

Get Radiology Tree app to read full this article<

Figure 9, Interface used to make manual measurements of the cross-sectional area. (a) The initial plane orthogonal to the x axis. (b) The manually adjusted plane.

Get Radiology Tree app to read full this article<

Patient Study

Get Radiology Tree app to read full this article<

Figure 10, Back view of segmented left atrium and pulmonary veins from a patient at (a) baseline and (b) 3-month follow-up: left inferior pulmonary vein ( green ), left superior pulmonary vein ( red ), right inferior pulmonary vein ( blue ), and right superior pulmonary vein ( magenta ).

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Results

Validation Results

Get Radiology Tree app to read full this article<

Table 1

Mean Absolute Differences between Manual and Automatic Measurements for Each of the Four Pulmonary Veins across 10 Subjects

Subject LIPV (mm 2 ) LSPV (mm 2 ) RIPV (mm 2 ) RSPV (mm 2 ) Mean (mm 2 ) 1 5.75 8.35 3.71 9.59 6.85 2 7.84 9.80 6.93 14.40 9.74 3 7.15 21.33 7.70 8.89 11.27 4 10.86 26.67 7.17 12.87 14.39 5 18.75 18.35 17.19 11.68 16.49 6 8.42 9.83 12.38 10.88 10.38 7 16.28 10.62 18.82 8.24 13.49 8 4.11 8.35 5.94 14.87 8.32 9 8.82 18.41 18.06 15.30 15.15 10 13.95 9.17 7.25 5.17 8.89 Mean 10.19 14.09 10.52 11.19 11.50

LIPV, left inferior pulmonary vein; LSPV, left superior pulmonary vein; RIPV, right inferior pulmonary vein; RSPV, right superior pulmonary vein.

Figure 11, Plots of manual ( blue ) and automatic ( red ) cross-sectional area measurements. Circles indicate locations of manual measurements. LIPV, left inferior pulmonary vein; LSPV, left superior pulmonary vein; RIPV, right inferior pulmonary vein; RSPV, right superior pulmonary vein.

Get Radiology Tree app to read full this article<

Patient Study

Get Radiology Tree app to read full this article<

Figure 12, (a) Average cross-sectional area plots at baseline ( blue ) and follow-up ( green ) for all 21 patients. The number of patients included in calculation is shown next to each data point. (b) Mean of paired differences in cross-sectional area between baseline and follow-up scans. Significant differences ( P < .01) are indicated by red stars and error bars indicate confidence intervals. PV, pulmonary vein.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Discussion

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Acknowledgments

Get Radiology Tree app to read full this article<

Appendix

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

References

  • 1. Haïssaguerre M., Jaïs P., Shah D., et. al.: Spontaneous initiation of atrial fibrillation by ectopic beats originating in the pulmonary veins. N Engl J Med 1998; 339: pp. 659-666.

  • 2. Kato R., Lickfett L., Meininger G., et. al.: Pulmonary vein anatomy in patients undergoing catheter ablation of atrial fibrillation. Circulation 2003; 107: pp. 2004-2010.

  • 3. Knackstedt C., Visser L., Plisiene J., et. al.: Dilation of the pulmonary veins in atrial fibrillation: a transesophageal echocardiographic evaluation. Pacing Clin Electrophysiol 2003; 26: pp. 1371-1378.

  • 4. Tsao H.-M., Yu W.-C., Cheng H.-C., et. al.: Pulmonary vein dilation in patients with atrial fibrillation: detection by magnetic resonance imaging. J Cardiovasc Electrophysiol 2001; 12: pp. 809-813.

  • 5. Ueda M., Tada H., Kurosaki K., et. al.: Pulmonary vein morphology before and after segmental isolation in patients with atrial fibrillation. Pacing Clin Electrophysiol 2005; 28: pp. 944-953.

  • 6. Scharf C., Sneider M., Case I., et. al.: Anatomy of the pulmonary veins in patients with atrial fibrillation and effects of segmental ostial ablation analyzed by computed tomography. J Cardiovasc Electrophysiol 2003; 14: pp. 150-155.

  • 7. Tsao H., Wu M., Huang B., et. al.: Morphologic remodeling of pulmonary veins and left atrium after catheter ablation of atrial fibrillation: insight from long-term follow-up of three-dimensional magnetic resonance imaging. J Cardiovasc Electrophysiol 2005; 16: pp. 7-12.

  • 8. Lemola K., Sneider M., Desjardins B., et. al.: Effects of left atrial ablation of atrial fibrillation on size of the left atrium and pulmonary veins. Heart Rhythm 2004; 1: pp. 576-581.

  • 9. Jayam V., Dong J., Vasamreddy C., et. al.: Atrial volume reduction following catheter ablation of atrial fibrillation and relation to reduction in pulmonary vein size: an evaluation using magnetic resonance angiography. J Interv Card Electrophysiol 2005; 13: pp. 107-114.

  • 10. Jeevanantham V., Ntim W., Navaneetham S., et. al.: Meta-analysis of the effect of radiofrequency catheter ablation on left atrial size, volumes and function in patients with atrial fibrillation. Am J Cardiol 2010; 105: pp. 1317-1326.

  • 11. Reant P., Lafitte S., Jais P., et. al.: Reverse remodeling of the left cardiac chambers after catheter ablation after 1 year in a series of patients with isolated atrial fibrillation. Circulation 2005; 112: pp. 2896-2903.

  • 12. denUijl D., Tops L., Delgado V., et. al.: Effect of pulmonary vein anatomy and left atrial dimensions on outcome of circumferential radiofrequency catheter ablation for atrial fibrillation. Am J Cardiol 2011; 107: pp. 243-249.

  • 13. Karim R., Mohiaddin R., Rueckert D.: Left atrium pulmonary veins: segmentation and quantification for planning atrial fibrillation ablation. Proc SPIE 2009; pp. 7261.

  • 14. Rettmann M., Holmes D., Packer D., et. al.: Identification of left pulmonary vein ostia using centerline tracking. Proc SPIE Medl Imaging 2009; pp. 7262.

  • 15. Rettmann M., Gunawan M., Holmes D., et. al.: Quantification of pulmonary vein morphology using centerline tracking. Proc Int Symp Biomed Imaging 2012;

  • 16. Rettmann M., Holmes D., Camp J., et. al.: Validation of semi-automatic segmentation of the left atrium. Proc SPIE 2008; pp. 6916.

  • 17. Rettmann M., Holmes D., Gunawan M., et. al.: Validation of geometric measurements of the left atrium and pulmonary veins for analysis of structural remodeling following ablation therapy. Proc SPIE 2012;

  • 18. Sadleir R., Whelan P.: Fast colon centreline calculation using optimised 3D topological thinning. Comput Med Imaging Graph 2005; 29: pp. 251-258.

  • 19. Dijkstra E.: A note on two problems in connexion with graphs. Numer Math 1959; 1: pp. 269-271.

  • 20. Kim Y.-H., Marom E., Herndon J., et. al.: Pulmonary vein diameter, cross-sectional area, and shape: CT analysis. Radiology 2005; 235: pp. 43-50.

  • 21. Akutsu Y., Kaneko K., Kodama Y., et. al.: Association between left and right atrial remodeling with atrial fibrillation recurrence after pulmonary vein catheter ablation in patients with paroxysmal atrial fibrillation: a pilot study. Circ Cardiovasc Imaging 2011; 4: pp. 524-531.

  • 22. Hof I., Chilukuri K., A-Zadeh A., et. al.: Does left atrial volume and pulmonary venous anatomy predict the outcome of catheter ablation of atrial fibrillation?. J Cardiovasc Electrophysiol 2009; 20: pp. 1005-1010.

  • 23. Nedios S., Tang M., Roser M., et. al.: Characteristic changes of volume and three-dimensional structure of the left atrium in different forms of atrial fibrillation: predictive value after ablative treatment. J Interv Card Electrophysiol 2011; 32: pp. 87-94.

  • 24. Parikh S., Jons C., McNitt S., et. al.: Predictive capability of left atrial size measured by CT, TEE, and TTE for recurrence of atrial fibrillation following radiofrequency catheter ablation. Pacing Clin Electrophysiol 2010; 33: pp. 532-540.

  • 25. Mulder A., Wijffels M., Wever E., et. al.: Pulmonary vein anatomy and long-term outcome after multi-electrode pulmonary vein isolation with phased radiofrequency energy for paroxysmal atrial fibrillation. Europace 2011; 13: pp. 1557-1561.

This post is licensed under CC BY 4.0 by the author.