Shan Gao1, Ronald van't Klooster1, Anne Brandts2, Stijntje D. Roes2, Reza Alizadeh Dehnavi3, Albert de Roos2, Jos J.M. Westenberg1, and Rob J. van der Geest1
1Division of Image Processing, Department of Radiology, Leiden University Medical Center, Leiden, Netherlands, 2Department of Radiology, Leiden University Medical Center, Leiden, Netherlands, 3Department of Cardiology, Leiden University Medical Center, Leiden, Netherlands
Synopsis
Atherosclerosis is characterized by thickening
of the arterial wall. To monitor disease progression and regression, vessel
wall segmentation is required for wall thickness assessment. In this study, we developed a fully automatic
method to detect the vessel wall boundaries and quantify the wall thickness for
both the common carotid artery and the descending aorta in MR images. The
results demonstrated that our method has good agreement with manual
segmentation, excellent scan-rescan reproducibility and capability of detecting
vessel wall thickening in hypertensive patients. Therefore, the presented
method is promising for application in future cohort studies.INTRODUCTION:
Aortic and carotid vessel wall MRI have been applied in
clinical trials to monitor the progression and regression of atherosclerosis by
assessing the wall thickness (WT). The analysis of these data requires vessel
wall segmentation and several studies have reported automated methods
[1,2,3].
However, most of the methods require user interaction as an input to the
automated process. In addition, all these segmentation algorithms are dedicated
to either the carotid artery or the aorta.
Aim: To develop a method that can fully automatically
segment the vessel wall and calculate WT for both the common carotid artery and
the descending aorta in MR images, which is especially suitable for cohort studies
where the detection of temporal WT changes is needed in a large dataset.
METHODS
AND MATERIALS:
MRI: In 10 healthy volunteers repeated MRI scans of the carotid artery were
acquired and in another 10 healthy volunteers repeated MRI scans of aorta. In 10
patients with essential hypertension and 10 age-matched and sex-matched healthy
volunteers one MRI scan of both the carotid artery and the aorta was acquired. All vessel wall images were acquired using a
3.0T scanner with an ECG triggered DIR (black-blood) segmented k-space
gradient-echo sequence with fat suppression. Cross-sectional images of the left common carotid were acquired using
the following parameters: TR/TE = 12/3.6 ms,
FA= 45◦, FOV = 140×140 mm
2, resolution = 0.27×0.27×2 mm
3. Cross-sectional images of descending aorta were
acquired with: TR/TE = 4.9/2.5 ms, FA = 20◦, FOV = 270×270 mm
2,
resolution = 0.53×0.53×2 mm
3.
Automated segmentation: The proposed
algorithm includes the following four steps which are all performed
automatically:
Step 1: localization
of the vessel of interest based on the circular Hough Transform (CHT) (Figure 1
b, e);
Step 2: estimation of the
lumen center and radius using oriented-edge based Hough Transform (OEHT)
(Figure 1 c, f);
Step 3: initialization
of the tube model for lumen followed by detection and refinement of the lumen
boundary based on a tube-fitting method
[1] and dynamic programing
[2];
Step 4: initialization of the tube
model for outer vessel wall based on the dilation of lumen segmentation
followed by segmentation of the outer wall using a similar manner as for lumen.
Manual
segmentation: Two trained radiologists (observer 1 and observer 2)
delineated the vessel wall contours using VesselMASS software
[2]. The
delineations performed by observer 1 served as the gold standard, and those
made by observer 2 were used to quantify inter-observer variability.
Evaluation: WT was calculated by taking the median length of
100 chords connecting lumen and outer wall boundaries
[1]. 1) To
directly compare the manual and automatic generated contours, the mean contour
distance (MCD) was calculated by taking the average of the distances between
the two sets of contour points. 2) To evaluate the performance and inter-scan
reproducibility of the automated WT assessment, Bland-Altman analysis and
paired t-test were performed. 3) To compare the WT between the healthy
volunteers and hypertensive patients, unpaired T-test was computed. For all statistical
tests, a P < 0.05 was considered to be statistically significant.
RESULTS:
Generally, all the mean differences between the manual
and automated segmentation results were smaller than one pixel (Table 1). The
vessel wall thickness was slightly, but significantly underestimated by the
automated segmentation for carotid and aorta; however, the difference between
observer 1 and observer 2 was also statistically significant (Table 1). For
both carotid and aorta, no significant difference was observed between the
inter-scan wall thickness assessments using automated segmentation (Table 1). Automated
segmentation detected significantly higher carotid and aortic wall thickness in
the hypertensive patients (Table 2).
CONCLUSION:
We presented a fully automated contour detection and
vessel wall thickness quantification method. The results demonstrated that our
method has good agreement with manual segmentation, excellent scan-rescan
reproducibility and capability of detecting vessel wall thickening in
hypertensive patients. Therefore, the presented method is promising for application
in future cohort studies.
Acknowledgements
No acknowledgement found.References
[1]van’t Klooster et al., J Magn Reson Imaging 2012. [2]Adame et al., J Magn Reson
Imaging 2006. [3] Underhill et al., J Magn Reson Imaging 2006.