Research  Open  Published:
A vessel segmentation method for multimodality angiographic images based on multiscale filtering and statistical models
BioMedical Engineering OnLinevolume 15, Article number: 120 (2016)
Abstract
Background
Accurate segmentation of blood vessels plays an important role in the computeraided diagnosis and interventional treatment of vascular diseases. The statistical method is an important component of effective vessel segmentation; however, several limitations discourage the segmentation effect, i.e., dependence of the image modality, uneven contrast media, bias field, and overlapping intensity distribution of the object and background. In addition, the mixture models of the statistical methods are constructed relaying on the characteristics of the image histograms. Thus, it is a challenging issue for the traditional methods to be available in vessel segmentation from multimodality angiographic images.
Methods
To overcome these limitations, a flexible segmentation method with a fixed mixture model has been proposed for various angiography modalities. Our method mainly consists of three parts. Firstly, multiscale filtering algorithm was used on the original images to enhance vessels and suppress noises. As a result, the filtered data achieved a new statistical characteristic. Secondly, a mixture model formed by three probabilistic distributions (two Exponential distributions and one Gaussian distribution) was built to fit the histogram curve of the filtered data, where the expectation maximization (EM) algorithm was used for parameters estimation. Finally, threedimensional (3D) Markov random field (MRF) were employed to improve the accuracy of pixelwise classification and posterior probability estimation. To quantitatively evaluate the performance of the proposed method, two phantoms simulating blood vessels with different tubular structures and noises have been devised. Meanwhile, four clinical angiographic data sets from different human organs have been used to qualitatively validate the method. To further test the performance, comparison tests between the proposed method and the traditional ones have been conducted on two different brain magnetic resonance angiography (MRA) data sets.
Results
The results of the phantoms were satisfying, e.g., the noise was greatly suppressed, the percentages of the misclassified voxels, i.e., the segmentation error ratios, were no more than 0.3%, and the Dice similarity coefficients (DSCs) were above 94%. According to the opinions of clinical vascular specialists, the vessels in various data sets were extracted with high accuracy since complete vessel trees were extracted while lesser nonvessels and background were falsely classified as vessel. In the comparison experiments, the proposed method showed its superiority in accuracy and robustness for extracting vascular structures from multimodality angiographic images with complicated background noises.
Conclusions
The experimental results demonstrated that our proposed method was available for various angiographic data. The main reason was that the constructed mixture probability model could unitarily classify vessel object from the multiscale filtered data of various angiography images. The advantages of the proposed method lie in the following aspects: firstly, it can extract the vessels with poor angiography quality, since the multiscale filtering algorithm can improve the vessel intensity in the circumstance such as uneven contrast media and bias field; secondly, it performed well for extracting the vessels in multimodality angiographic images despite various signalnoises; and thirdly, it was implemented with better accuracy, and robustness than the traditional methods. Generally, these traits declare that the proposed method would have significant clinical application.
Background
Nowadays, cardio and cerebrovascular diseases have greatly threatened human health. Since the use of imaging techniques such as computed tomography angiography (CTA) and magnetic resonance angiography (MRA) in minimally invasive surgery, high quality image segmentation has become an important area of interest. A detailed review has been made on threshold based, pattern recognition based, and deformable models based segmentation algorithms which were used for medical images [1, 2]. The main tendency of these algorithms with their principle ideas, application field, advantages and disadvantages were discussed. Conclusion has been drawn that each segmentation method with improvement, or in combination with other technique, could provide better performance [3, 4]. Clustering algorithms and supervised classification used for the segmentation of atherosclerotic plaques were analysed in [5].
For the diagnosis and treatment of vascular diseases, it is critical to accurately extract and quantify the blood vessels from the angiographic image. Aimed at the extraction of blood vessels, different kinds of segmentation methods have been proposed, e.g.: multiscale filtering [6–9], deformable models [10], statistical models [11–13], and hybrid methods [14, 15]. Multiscale filtering can enhance the intensity of the vessels while suppresses that of the background [16] in various angiographic data, however, further processing should be conducted in order to label the vessel class out. Deformable models well integrate bottom–up information and top–down priori knowledge, but the segmentation quality mainly depends on the model parameters [17, 18]. For statistical models, the vessels in the angiographic data with a given modality are segmented according to the intensity distributions of anatomical structures, by which voxels with overlapped intensities will inevitably be misclassified [19]. In order to reduce the probability of error classification, hybrid methods have been proposed to take the spatial contextual information into account [20], but limitations still exist for multimodality angiographic images.
So far, statistical models have drawn a lot of attention, and model selection is an important issue in this kind of vessel segmentation techniques. To our known, the earliest statisticalbased segmentation algorithm for extracting 3D cerebral vessels from MRA data was established by Wilson and Noble [19]. They divided the MRA data histogram into three regions, modelled by a mixture model of two Gaussian distributions and one uniform distribution. Hassouna et al. [21] improved the mixture model with a Rayleigh distribution and three Gaussian distributions to provide an accurate fitting of the histogram curve. On the basis of statistical model analysis and improved curve evolution, a fast segmentation algorithm for extracting the 3D cerebral vessels was presented by Gao et al. [22]. To segment fine cerebral vessels with complicated contexts, a statistical method based on maximum a posteriori and Markov random field (MAP–MRF) with multipattern neighbourhood system (MPNBS) was proposed by Zhou et al. [23], where the larger cerebrovascular network can be extracted out automatically. In general, statistical segmentation methods need to construct mixture probability distributions and adjust their parameters to fit a certain image histogram, where the vessels in a particular position are segmented according to the specific intensity distribution of a certain angiography image. To sum up, a given statistical model only deals with a specific imaging modality on a given anatomical part, which is to say that the statistical segmentation method should change its models to fit a certain medical data. This issue has greatly restricted the application of statistical methods.
Being different from the abovementioned statistical vessel segmentation method, in this paper, a flexible method for extracting vessels from various modalities of angiographic images was developed. Up to present, statistical modelling against the data of multiscale filtering response was not been used to vessel segmentation from angiographic data. One reason lies in the complexity of constructing mixture model, and the other lies in the error being not easily controlled with parameters estimation. We have done this attempt, and our method includes the following several steps. In the first step, image data of different modalities were filtered by multiscale vessel enhancement algorithm. The intensity of the vessels was enhanced, whereas the intensity of the background voxels was suppressed, and the filtered data turned to a new distribution. Then, the filtered data was modelled by a new mixture model consisting of reduced distribution classes. Based on the histogram curve analysis, a Gaussian distribution was used to model the high intensity which represents the vessels; while two Exponential distributions were used to model the low and middle intensity regions which represent the background voxels. The parameters of the mixture model were estimated by the EM algorithm. At last, the vessels were marked out by MAP–MRF algorithm.
The rest of this paper is organized as follows. In “Methods” section, the proposed method is described in detail, where the multiscale filtering based vessel enhancement, histogram analysis, the mixture model constructing for histogram curve fitting, parameters estimation, and the MAP–MRF segmentation algorithm are presented in succession. “Experiments” section details the experiments that were conducted on a series of phantoms, the MRA data, and the CTA data, where results are also presented. Summaries and future work are discussed in “Results and discussions” section, and the conclusions are drawing in the last section.
Methods
The framework of the proposed method is shown in Fig. 1. Firstly, to attain a relatively consistent histogram characteristic from multimodality input data, multiscale filtering is employed as the preprocessing step. Secondly, a mixture model composed of three distributions was built to fit the histogram curve of the filtered data. The parameters of the three distributions were automatically estimated through peakpoints detecting of the histogram curve, Kmeans clustering algorithm, and EM algorithm. Finally, the blood vessels were segmented by the MAP–MRF algorithm.
Vessel enhancement algorithm
We use the algorithm of three dimensional multiscale filtering algorithm [6] to strengthen the intensity of the vessels while weakening that of the background. According to the eigenvalues of Hessian, the vessel function with multiscale filtering is expressed as:
where R _{ A }, R _{ B } and S are the three measures that can be respectively controlled by the parameters α, β and c, so that the vessel function be sensitive to the vessel shapes. The relationships between these parameters and the eigenvalues are shown in the following:
Here, h _{1}, h _{2}, and h _{3} are the eigenvalues corresponding to the eigenvectors of three orthonormal directions which are robust to the scaling factor. The resultant ratios R _{ A } and R _{ B } would be sensitive to the shapes of the blood vessel structure (e.g., blob, plate, and linelike structures), while the measure S would be used to suppress image noise.
The abovementioned enhancement algorithm was tested on several angiographic data sets. The histogram curves of the original data sets and those of the corresponding filtered data sets are compared in Fig. 2. Note that the low, middle, and high intensity regions of the filtered data gained from various modality original data exhibit similar distribution features.
After processed using multiscale filtering, the vessels would be enhanced while the background tissues would be weakened, such that the resultant histogram curve of the filtered data would be fitted by reduced model elements.
Mixture model construction
Through the comparison of the histogram curves between the original and the filtered data in Fig. 2, it can be seen obviously that the vessels occupy the high intensities region with a range of gradual change, thus a Gaussian distribution could be used to model the vessel class. For the background class, the Exponential distributions showed the best fitting performance after several probability distribution functions were tested. Therefore, two Exponential distributions were used to model the low and middle intensity regions of background voxels. As a result, the mixture model for the filtered data can be described as a linear combination of the three probability distribution functions:
where the mixture model f(x) is composed of Exponential distributions f _{ E1}(x) and f _{ E2}(x) as well as Gaussian distribution f _{ G }(x). The parameters w _{ E1}, w _{ E2}, and w _{ G } are the proportion coefficients whose sum is unity. The three distributions can be defined as follows:
where parameters λ, μ, σ, and the proportion coefficients w _{ El }, w _{ G } should be firstly estimated before the classification procedure. To begin with, the intensity range of the filtered data was enlarged to be consistent with the original image data. Then, the filtered data was classified into three parts using the Kmeans algorithm. For the initialization of the Kmeans algorithm, the histogram curve of the original image was smoothed, and the peak points were automatically found out through local extremum method. If there were more than three peak points in the original histogram curve, the Kmeans algorithm would be initialized by the first three peak values. Otherwise, it would be initialized according to the following rule:
where I _{max} and I _{min} are the maximum and the minimum intensity values of the filtered data. The variables \(\mu_{k\_1}^{init}\), \(\mu_{k\_2}^{init}\), and \(\mu_{k\_3}^{init}\) are the initial cluster centroids while μ _{ k_l }, σ _{ k_l }, w _{ k_l }, (l = 1, 2, 3) represent the distribution parameters which were all estimated by the Kmeans algorithm.
As the accuracy of theses parameters plays an important role in the final classification, the EM algorithm was used to improve the estimation performances of these parameters. The results of the Kmeans algorithm were used as the initial values for the EM algorithm, as shown in Eq. (7):
Then, the parameters were iteratively updated with the EM algorithm as follows:
where N is the total number of the voxels in the filtered data and x _{ j } is the corresponding intensity of voxel j. The posterior probability f ^{k}(Elx _{ j }) and f ^{k}(Gx _{ j }) are gained by the Bayes formula that are explained in Eq. (10)
MAP–MRF based classification
A series of filtered data was modelled by the proposed mixture model, and the combination of two Exponential distributions and one Gaussian distribution shows a good fitting to the histogram curves of the filtered data from the multimodality angiographic images. According to the MAP algorithm, a voxel belongs to the vessel class only if it meets the following criteria:
Equation (11) gives the simplest case of binary vessel labelling, in which the voxels of the filtered data are labelled as either blood vessels or background according to their posteriori probability. However, misclassification still exists when significant vascular signal losses or the noise has high intensity value because MAP classification only depends on the voxel intensity. Therefore, MRF algorithm is introduced for enhanced segmentation.
Spatial interaction between neighbouring voxels is modelled by MRF. Assume that X is an MRF, observed vector Y is conditionally independent, then MAP–MRF based expressions can be written as
where p(XY), being the low level MRF process, is modelled as a mixture model talked above; the term p(X), being the high level MRF process, is given by the Gibbs distribution [24] as follows
where U(x) is the energy function, Ω denotes the space of all possible labelling, and Z is a normalizing constant to be the sum of the numerator over all possible configurations of x. With the isotropic multilevel logistic (MLL) model [25], U(x) can be expressed as
where η denotes the neighbourhood cube with the size of \(3 \times 3 \times 3\), \(V(x_{s} ,x_{r} )\) is the energy of pairpoints in this MPNBS [23], and the energy function U(x _{ s }) is proportional to the summation of the weighted \(V(x_{s} ,x_{r} )\) with respect to β _{ sr }. The instances of \(V(x_{s} ,x_{r})\) used are given below.
The regularization parameter β _{ sr } is meant for the strength of the interaction between pairwise neighbouring voxels.
Overall, according to the MAP–MRF approach, the posterior probability of the vessel class can be expressed by
Meanwhile, the posterior probability of the background class is given by
Finally, the MAP–MRF model is solved by iteratively using the iterated conditional modes (ICM) [26], and a voxel belongs to the vessel class only if it meets the following rule:
Experiments
The experiment results of the proposed segmentation method are presented in this section. In order to quantitatively and qualitatively validate the proposed method, tests were conducted on phantoms and various modality angiographic images from MRA, and CTA data. Additionally, comparisons with the traditional methods have also been discussed at the end of this section.
Experiment of quantitative validation
Since it was not easy to acquire the golden standard of vessels in real clinical data, phantoms with known ground truth and complicated background noises were designed to quantitatively validate the proposed method. The percentage of the number of misclassified voxels to that of the total voxels of the phantom was calculated as the segmentation error ratio. As a popular metric for evaluating the accuracy of automated or semiautomated segmentation methods, the Dice similarity coefficients (DSCs) was used to estimate the segmentation performance on the phantom data by comparing their results to the ground truth.
Two 3D phantoms were built to simulate the real clinical vessels data, which were shown in Fig. 3.
The simulation processes of mixture distributions for target and background were referred to our approach by Zhou et al. [23]. The main characteristics lied in follows: (1) Considering that the components of the background voxels in the angiographic images are often Gauss signals, different levels of Gaussian noises w.r.t. \((\mu ,\sigma^{2} ,\omega )\) were added into the two phantoms. (2) The proportion of the ground truth is not limited in our method, i.e., the proportions of the target and the background could be randomly set by adjusting the size of their volumes. (3) The vessels occupy the high intensities region, thus, in the setting of parameters \(\text{(}\mu ,\sigma^{2} \text{, }\omega \text{)}\), the value μ of the vessel target should be greater than that of each Gauss noise.
Among the two phantoms one was a single curving tubular target, and the other consists of a series of tubular targets, where the intensity values of the two vessel objects were generated with stochastic ones of 810 and 650 respectively. The diameter of each tubular target smoothly changed, like that of a real vessel structure. Three Gaussian noises with the parameters \(\text{(}\mu ,\sigma^{2} ,\omega \text{)}\) designed as \((148,10^{2} ,0.4)\), \((187,15^{2} ,0.4)\), and \((228,60^{2} ,0.24)\) were added into the first phantom while two Gaussian noises with the parameters \(\text{(}148,10^{2} \text{, 0}\text{.5)}\), and \((187,15^{2} ,0.3)\) were added into the second phantom. The data volumes in the two phantoms were 128 × 128 × 128 and 240 × 240 × 240 voxels respectively. The 3D view of the ground truth and some slices of the phantoms were shown in Fig. 3.
In the first step, the phantoms were preprocessed with the multiscale filtering algorithm. After that, the vessels were marked out using the proposed segmentation method. For the filtered data of the phantoms, the estimated two Exponential distributions and one Gaussian distribution curves were presented in Fig. 4a together with their fitting curves. Slices of the segmentation results and their 3D views were plotted in Fig. 4b. Meanwhile, the segmentation error ratios and the DSCs of the two phantoms were listed in Table 1, where N _{ total } represents the total voxel number of the phantom, N _{ over }, and N _{ under } are the number of voxels in over segmentation and under segmentation respectively. Thus, the sum of N _{ over } and N _{ under } is the number of voxels in misclassification. S _{ g } is the set of pixels assigned to the targets in the ground truth, S _{ m } is the set of pixels assigned to the targets in the segmentation result, and \(N( \cdot )\) represents the number of pixels in the corresponding region. The DSC provided a measure of how accurate a method was, with its value closer to 1 indicating greater accuracy.
Experiment of clinical validation
This section detailed the qualitative validation that was conducted on four clinical data, each of which held a different statistical characteristic. Among the four clinical data, the brain MRA data set was offered by Guangzhou General Hospital of the Chinese People’s Liberation Army (PLA). The placenta MRA, legs MRA, and the thorax CTA data sets were gained from Shenzhen No.2 People’s Hospital. These clinical data sets were shown in Fig. 5, and the sizes and elements spacing of these data were listed in Table 2.
The images of different organs from various imaging equipment exhibited different statistical characteristics. After using the vessel enhancement algorithm, the histograms of the filtered data of these clinical data obtained a new relatively consistent distribution. For the filtered data, the vessel class was modelled by a Gaussian distribution, and the background class was modelled by two Exponential distributions. The results of distribution curves fitting are shown in Fig. 6.
For clinical validation with respect to the real angiographic data, the Zaxis projections as well as the 3D views of the segmentation results were compared to the corresponding maximal intensity projection (MIP) of the original images, as suggested by clinical experts on cardio and cerebrovascular diseases. Parts of the results are given in Fig. 7.
Experiment of comparison tests
As an amendment to abovementioned quantitative and clinical validation experiments, comparison tests were conducted against the traditional methods. Since our study was mainly based on statistical models, three popular modelbased segmentation methods were chose for the tests. The three algorithms were methods of Wilson et al. [20], Gao et al. [22], and Zhou et al. [23]. The mixture models used in the proposed method and the three traditional ones were different from each other, as shown in Table 3. Besides, the proposed mixture model was used to fit the filtered data, and the mixture models of the other three were used to fit the original data. To get the best segmentation results from each of the traditional methods, brain MRA data was used for the tests according to their application ranges. During the validation, the three traditional methods were implemented using models specified in the original papers and parameters as best as possible.
The online brain MRA datasets (http://www.insightjournal.org/midas/community/view/21) have previously been used by many researchers to test their vessel segmentation methods. They were acquired at the element spacing of 0.5 × 0.5 × 0.8 mm^{3} with the data volumes of 448 × 448 × 128 voxels. Comparison tests were conducted on the ten datasets, and the experimental results of a single MRA data from the four segmentation methods are displayed in Fig. 8.
To further verify the effectiveness of the proposed method, comparison tests were also conducted on some data sets offered by the aforementioned Guangzhou General Hospital of the Chinese PLA. The image acquisition parameters of these data together with those of the online brain MRA data sets were listed in Table 4. Because of uneven contrast media and bias field, these MRA data sets exhibited a poor image quality. One brain MRA data set and its histogram curve are shown in Fig. 9, from which we can see this kind of data hold a different statistical property.
Results and discussions
In this paper, multiscale filtering algorithm was firstly used for vessel enhancement and noise suppression before the statistical based classification, where filter parameters settings of α, β, c were same for all the datasets. Suggested by Frangi’s method [6], a compromise value of α and β equal to 0.5 would make multiscale filtering algorithm to differentiate tubularlike structures from plate, and bloblike ones; while the parameter c depended on the greyscale range of the image and was set as half the value of the maximum Hessian norm. The filtered data from multimodality angiographic images turns to a relatively consistent statistical distribution after the vessel enhancement algorithm. The following statistical model used two Exponential distributions and one Gaussian distribution to fitting the relatively consistent statistical distribution, which classified the vessels through adjusting model components to fit a certain histogram curve of the filtered data. To sum up, as a flexible segmentation method, it performs well in vessel enhance and histogram fitting for the filtered data and is robust to various angiographic images from different human organs and different imaging equipment.
For phantom validation, the method performed well for images despite various signal noises. The tubular object was enhanced while background noise was suppressed effectively. The accuracy of the proposed method was tested by several 3D phantoms. The phantoms were built to simulate the real angiographic images where different noises were added to simulate the background voxels. Compared to the known ground truth, the mixture model efficiently fitted the histogram curves of both phantoms with the segmentation error ratios being 0.21 and 0.11%, and the DSCs being 94.74 and 96.94%. See Table 1, which strongly confirmed the segmentation quality.
For clinical validation, the proposed method was tested on multimodality angiographic images from different human organs. The mixture model in our method fitted well in each histogram curve of the filtered data. Comparing with the MIP from the original images and the segmentation results, the complete vessel trees were extracted while lesser nonvessels and background were falsely classified as vessel. The results were satisfactory according to the suggestion of clinical experts. Parts of the results are given in Fig. 7.
In comparison tests, the results from traditional statistical based methods and the proposed method were presented. Experimental results of a single MRA data from the four segmentation methods are displayed in Fig. 8. The mixture model built by Zhou et al. obtained the best fitting performance among the three traditional methods according to the histogram fitting performance and the Zaxis projection of these segmentation results. In their method, a Gaussian distribution was used to model the vessel class, while a Rayleigh distribution and two Gaussian distributions were used to model the background class. At the same time, the mixture model built in this paper fitted well for the histogram of the filtered data, and more small vessels were marked out which can been seen from the Zaxis projection of its segmentation result. From Figs. 8 and 9 we can see the mixture model built in Zhou’s method exhibited a good performance in histogram fitting for the online MRA data, but not for the second MRA data which contains a different statistical characteristic. In contrast, the mixture model generated from our method fitted well for both filtered data whose histograms had a relatively consistent statistical property. As a result, the proposed method demonstrated stronger robustness.
For fitting the histogram curves of the original images, and the mixture models developed by Zhou et al. exhibited the best fitting performance comparing with the traditional models. However, different angiographic images correspond to different histogram curves. The four clinical data sets from different organs that were described in the part of clinical validation experiment were segmented by Zhou et al.’s method, and the model fitting results are shown in Fig. 10. As a result, the traditional method is not available to classify vessels from angiographic images which hold a different statistical characteristic. In other words, the traditional methods are limited in multimodality application.
On the basis of Zhou et al.’s work, our research focuses on vessel segmentation from various angiographic images. In our method, different angiographic images were firstly processed by multiscale filtering algorithm. As shown in Fig. 2, the filtered data sets achieved a relatively consistent distribution. For the vessel class, a Gaussian distribution was built to model its intensity distribution, while two Exponential distributions for the background class. As is shown in Figs. 4, 6, 8, and 9, the above mixture model demonstrated good model fitting performance for the different filtered data sets. As a result, the only known limitation of our approach is that, it is not available for the angiographic images when multiscale filtering is not fitted for their vessel enhancement. However, the proposed method still exhibits important clinical application value towards various angiographic images, and this kind of exploration ideas is very meaningful in other scientific researches. With the aim of contributing to the field of interventional surgery, our future works are primarily focused on vascular path planning and 3D augmented reality.
Conclusions
In this paper, a flexible segmentation method has been proposed for vessel extraction. In the first step, the original image data was filtered with multiscale filtering algorithm. The filtered data obtained a relatively consistent statistical property, and was modelled with a mixture model which consisted of two Exponential distributions and one Gaussian distribution. The EM algorithm was used for parameter estimation, while the MAP–MRF algorithm was used in the Bayesian classification. The proposed method has been tested on phantoms as well as the various modality angiographic images from MRA and CTA data from different imaging equipment. Each histogram of the filtered data from these multimodality angiographic images could be modelled by relatively consistent mixture probability distributions, and the mixture model displayed a good fitting performance. The segmentation error ratios of the phantoms were less than 0.3% and the DSCs were above 94%. The results of the vessel segmentation procedures have indicated that the method in this paper was accurate and robust for vessel segmentation from multimodality angiographic images.
Abbreviations
 EM:

expectation maximization
 3D:

threedimensional
 MAP–MRF:

maximum a posteriori and Markov random field
 MAP:

maximum a posteriori
 MRF:

Markov random field
 CTA:

computed tomography angiography
 MRA:

magnetic resonance angiography
 MLL:

multilevel logistic
 MPNBS:

multipattern neighbourhood system
 ICM:

iterated conditional modes
 PLA:

People’s Liberation Army
 MIP:

maximal intensity projection
References
 1.
Ma Z, Tavares JMRS, Jorge RMN. A review on the current segmentation algorithms for medical images. In: Int. conf. on imaging theory and applications. 2009.
 2.
Ma Z, Tavares JMRS, Jorge RN, Mascarenhas T. A review of algorithms for medical image segmentation and their applications to the female, pelvic cavity. Comput Methods Biomech Biomed Eng. 2010;13(2):235–46.
 3.
Oliveira RB, Filho ME, Ma Z, Papa JP, Pereira AS, Tavares JMRS. Computational methods for the image segmentation of pigmented skin lesions: a review. Comput Methods Programs Biomed. 2016;131:127–41.
 4.
Filho M, Ma Z, Tavares JMRS. A review of the quantification and classification of pigmented skin lesions: from dedicated to handheld devices. J Med Syst. 2015;39(177):1–12.
 5.
Jodas DS, Pereira AS, Tavares JMRS. A review of computational methods applied for identification and quantification of atherosclerotic plaques in images. Expert Syst Appl. 2015;46:1–55.
 6.
Frangi AF, Niessen WJ, Vincken KL, Viergever MA. Multiscale vessel enhancement filtering. In: First international conference on medical image computing and computerassisted intervention. Berlin: Springer; 1998. p. 130–37.
 7.
Zhou C, Chan HP, Sahiner B, Hadjiiski LM, Chughtai A, Patel S, Wei J, Ge J, Cascade PN, Kazerooni EA. Automatic multiscale enhancement and segmentation of pulmonary vessels in ct pulmonary angiography images for cad applications. Med Phys. 2007;34(12):4567–77.
 8.
Ye DH, Kwon D, Yun ID, Sang UL. Fast multiscale vessel enhancement filtering. In: Medical imaging. Bellingham: International Society for Optics and Photonics; 2008. p. 691423.
 9.
Chen C, Zhang JZ. Improved multiscale Retinex enhancement algorithm of medical X ray images. Comput Eng Appl. 2015;9:191–5.
 10.
Mesejo P, Valsecchi A, MarrakchiKacem L, Cagnoni S, Damas S. Biomedical image segmentation using geometric deformable models and metaheuristics. Comput Med Imaging Graph. 2015;43:167–78.
 11.
Gan R, Wong WCK, Chung ACS, Yu SCH. Statistical cerebrovascular segmentation in threedimensional rotational angiography based on maximum intensity projections. Med Phys. 2005;32(9):195–200.
 12.
Xiong TS, Huang YY. Robust gaussian mixture modelling based on spatially constraints for image segmentation. J Inf Hiding Multimed Signal Process. 2015;6(5):857–68.
 13.
Zhang H, Wen T, Zheng Y, Xu D, Wang D, Nguyen TM. Two fast and robust modified Gaussian Mixture models incorporating local spatial information for image segmentation. J Signal Process Syst. 2015;81(1):1–14.
 14.
Hong Q, Li Q, Tian J. Local hybrid levelset method for MRA image segmentation. In: 2010 IEEE 10th international conference on computer and information technology (CIT). New York: IEEE; 2010. p. 1397–402.
 15.
Qian HY, Zhou WP, Bao XD, Luo LM. MRA blood vessel segmentation based on anticoncentration diffusion method and MumfordShah model relating to grayscale histogram. Signal Process. 2007;23(2):273–7.
 16.
Tang S, Lin T, Yang J, Fan J, Ai D, Wang Y. Retinal vessel segmentation using supervised classification based on multiscale vessel filtering and Gabor wavelet. J Med Imaging Health Inform. 2015;5(7):9–17.
 17.
Yan P, Kassim AA. Segmentation of volumetric MRA images by using capillary active contour. Med Image Anal. 2006;10(3):317–29.
 18.
Bauer C, Tsorantin P. Segmentation of interwoven 3d tubular tree structures utilizing shape priors and graph cuts. Med Image Anal. 2010;14(2):172–84.
 19.
Wilson DL, Noble JA. An adaptive segmentation algorithm of timeofflight MRA data. IEEE Trans Med Imaging. 1999;18(10):938–45.
 20.
Yousefi S, Liu T, Wang RK. Segmentation and quantification of blood vessels for OCTbased microangiograms using hybrid shape/intensity compounding. Microvasc Res. 2015;97:37–46.
 21.
Hassouna MS, Farag AA. Cerebrovascular segmentation from TOF using stochastic models. Med Image Anal. 2006;10:2–18.
 22.
Gao X, Uchiyama Y, Zhou X, Hara T, Asano T, Fujita H. A fast and fully automatic method for cerebrovascular segmentation on timeofflight (TOF) MRA image. J Digit Imaging. 2011;24(4):609–25.
 23.
Zhou SJ, Chen WF, Jia FC, Hu QM, Xie YQ, Chen MY, Wu JH. Segmentation of brain magnetic resonance angiography images based on MAP–MRF with multipattern neighborhood system. Med Image Anal. 2013;17:1220–35.
 24.
Besag J. Spatial interaction and the statistical analysis of lattice systems. J R Stat Soc B. 1974;36(2):192–236.
 25.
Derin H, Elliott H. Modeling and segmentation of noisy and textured images using Gibbs random fields. IEEE Trans Pattern Anal Mach Intell. 1987;9(1):39–55.
 26.
Besag J. On the statistical analysis of dirty pictures. J R Stat Soc Ser B. 1986;48:259–302.
Authors’ contributions
PL and SZ determined the research subject, and completed the design of the algorithm. JX is responsible for the acquisition and analysis of the human data. PL, SZ, MC and CW implemented the algorithm, and performed the experiments. ZL, JX, JY and LW provided important suggestions for the paper, and were responsible for the data analysis. All authors read and approved the final manuscript.
Acknowledgements
The study has been financially supported by the following grants: National High Technology Research and Development Program of China (863 Program) (Grant No. 2015AA043203), National Natural Science Foundation of China (Grant No. 61471349), Basic Discipline Layout Project of Shenzhen City (Grant No. JCYJ20150731154850923), and Technology innovation program of Shenzhen City (Grant No. JSGG20160229203812944).
Competing interests
The authors declare that they have no competing interests.
Availability of data and supporting materials
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The human data used in this research has been approved by the ethics committee of Shenzhen No.2 People’s Hospital, and the committee’s reference number is 20160701. We have been informed that the patients were consent to participate in the study.
Funding
The four sources of funding for this research reported are listed as follows: National High Technology Research and Development Program of China (863 Program) (Research and development for cerebral vascular realtime interventional therapy with mobile surgical robot system, Grant No. 2015AA043203), National Natural Science Foundation of China (Computerassisted analysis and clinical application research of cerebral arteriovenous malformations based on MRA, Grant No. 61471349), Basic Discipline Layout Project of Shenzhen City (Research of surgical vessel intervention with immersing multiperception mode, Grant No. JCYJ20150731154850923), and Technology innovation program of Shenzhen City (Research and development of minimally invasive spinal surgery system based on multimodality imaging and Intraoperative navigation, Grant No. JSGG20160229203812944). The first funding provides support in the design of this study, the second one and the third one are responsible for data collection, and the last one for data analysis.
Author information
Additional information
Pei Lu and Jun Xia contributed equally to this work
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Multimodality angiographic images
 Vessel segmentation
 Multiscale filtering
 Statistical mixture model
 Markov random field