We propose a region growing vessel segmentation algorithm based on spectrum information. First, the algorithm does Fourier transform on the region of interest containing vascular structures to obtain its spectrum information, according to which its primary feature direction will be extracted. Then combined edge information with primary feature direction computes the vascular structure’s center points as the seed points of region growing segmentation. At last, the improved region growing method with branch-based growth strategy is used to segment the vessels. To prove the effectiveness of our algorithm, we use the retinal and abdomen liver vascular CT images to do experiments. The results show that the proposed vessel segmentation algorithm can not only extract the high quality target vessel region, but also can effectively reduce the manual intervention. 1. Introduction Nowadays, the vessel segmentation technique [1–3] is still a bottleneck of medical image processing. As vascular has an extremely complex topology structure, making most of the conventional image segmentation methods difficult to segment vascular structures accurately, so how to fast, accurately and effectively segment vascular structures from medical images becomes an important issue. However, just using simple segmentation method cannot achieve expectant purpose. And effective vascular structure enhancement is also a key point for vascular segmentation. Hessian matrix [4] is often used in enhancement filter on tubular structures. Traditional single scale-based enhancement filter is not well adapted for vascular structures with large-scale changes. Some multiscale filters such as “cole” [5] and steerable filters [6], as well as using the Hessian matrix to determine the local direction of the target [7], however, do not have good and automatic scale selection approach. Region growing algorithm [8, 9] has small calculation complexity and high speed and is widely used in vascular image segmentation. The basic idea of the traditional growth region is to collect pixels that have similar properties together to form a region. Its performance depends largely on the position of seed points and growth conditions. On the other hand, the traditional method needs to select seed points manually in order to ensure the stability. Because the blood vessels have a very wide gray level distribution, making the traditional region growing algorithm and binary segmentation method difficult to accurately segment the vascular area, for the region growing conditions also, depends on range of the image gray
References
[1]
B. D. Thackray and A. C. Nelson, “Semi-automatic segmentation of vascular network images using a rotating structuring element (ROSE) with mathematical morphology and dual feature thresholding,” IEEE Transactions on Medical Imaging, vol. 12, no. 3, pp. 385–392, 1993.
[2]
S. R. Aylward and E. Bullitt, “Initialization, noise, singularities, and scale in height ridge traversal for tubular object centerline extraction,” IEEE Transactions on Medical Imaging, vol. 21, no. 2, pp. 61–75, 2002.
[3]
D. J. Stevenson, L. D. R. Smith, and G. Robinson, “Working towards the automatic detection of blood vessels in X-ray angiograms,” Pattern Recognition Letters, vol. 6, no. 2, pp. 107–112, 1987.
[4]
I. Jafar and H. Ying, “New algorithms for contrast enhancement in grayscale images based on the variational definition of histogram equalization,” Integrated Computer-Aided Engineering, vol. 15, no. 2, pp. 131–147, 2008.
[5]
S. R. Aylward and E. Bullitt, “Initialization, noise, singularities, and scale in height ridge traversal for tubular object centerline extraction,” IEEE Transactions on Medical Imaging, vol. 21, no. 2, pp. 61–75, 2002.
[6]
T. M. Koller, G. Gerig, G. Szekely, and D. Dettwiler, “Multiscale detection of curvilinear structures in 2-D and 3-D image data,” in Proceedings of the 5th International Conference on Computer Vision, pp. 864–869, June 1995.
[7]
Y. Sato, S. Nakajima, H. Atsumi et al., “3D multi-scale line filter for segmentation and visualization of curvilinear structures in medical images,” in Proceedings of the Lecture Notes in Computer Science, pp. 213–222, Springer, Berlin, Germany, 1997.
[8]
R. Adams and L. Bischof, “Seeded region growing,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 16, no. 6, pp. 641–647, 1994.
[9]
Preetha, “Image segmentation using seeded region growing,” in Proceedings of the International Conference on Computing, Electronics and Electrical Technologies (ICCEET '12), pp. 576–583, March 2012.
[10]
Z. Xizhi, “The application of wavelet transform in digital image processing,” in Proceedings of the International Conference on MultiMedia and Information Technology (MMIT '08), pp. 326–329, December 2008.
[11]
S. K. Alam, J. Mamou, E. J. Feleppa, A. Kalisz, and S. Ramachandran, “Comparison of template-matching and singular-spectrum-analysis methods for imaging implanted brachytherapy seeds,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 58, no. 11, pp. 2484–2491, 2011.
[12]
F. Chadiri and T. Akbarzadeh, “Vessel segmentation based on Sobel operator and fuzzy reasoning,” in Proceedings of the 1st International Conference on Computer and Knowledge Engineering (ICCKE '11), pp. 189–194, 2011.