Computational and Mathematical Methods in Medicine

Volume 2016 (2016), Article ID 5237191, 8 pages

http://dx.doi.org/10.1155/2016/5237191

## New Technique for Automatic Segmentation of Blood Vessels in CT Scan Images of Liver Based on Optimized Fuzzy -Means Method

Department of Computer Engineering, Faculty of Engineering, Arak Branch, Islamic Azad University, ARAK, Markazi, Iran

Received 19 July 2016; Revised 23 October 2016; Accepted 27 October 2016

Academic Editor: Thomas Desaive

Copyright © 2016 Katayoon Ahmadi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Automatic segmentation of medical CT scan images is one of the most challenging fields in digital image processing. The goal of this paper is to discuss the automatic segmentation of CT scan images to detect and separate vessels in the liver. The segmentation of liver vessels is very important in the liver surgery planning and identifying the structure of vessels and their relationship to tumors. Fuzzy -means (FCM) method has already been proposed for segmentation of liver vessels. Due to classical optimization process, this method suffers lack of sensitivity to the initial values of class centers and segmentation of local minima. In this article, a method based on FCM in conjunction with genetic algorithms (GA) is applied for segmentation of liver’s blood vessels. This method was simulated and validated using 20 CT scan images of the liver. The results showed that the accuracy, sensitivity, specificity, and CPU time of new method in comparison with FCM algorithm reaching up to 91%, 83.62, 94.11%, and 27.17 were achieved, respectively. Moreover, selection of optimal and robust parameters in the initial step led to rapid convergence of the proposed method. The outcome of this research assists medical teams in estimating disease progress and selecting proper treatments.

#### 1. Introduction

Liver is a vital organ in the body that cleans the body from toxins and harmful substances. Due to the function of liver, it is vulnerable to cancer. Resection and transplant of the tumor are two main treatments in routine clinical practice while both need planning and quality assessment of CT images for further processing.

Conventionally, detection of malignant tumors in CT scan images is visually applied by experienced physicians but it is time consuming and subject to human errors. Furthermore, CT scan images usually have low contrast and low visibility which increase the probability of false positive and false negative in cancer detection by human observers.

Comprehensive understanding of vascular anatomy of the liver is vital to cut segment-circuit which is derived from Couinaud’s descriptions about the liver segmentation. It is based on distribution of the venous shoots and the position of the liver vessels [1]. Liver surgery based on computer supported planning has a great effect on therapeutic strategies. An important step in the preoperation phase is to visualize the relationship between tumor, liver, and vascular tree of each patient [2].

#### 2. Related Works

So far, multiple image processing techniques have been introduced in the literature for automatic segmentation of liver vessels. Threshold method is perhaps the most important segmentation technique in medical image processing. Selection of proper threshold is the common problem of all threshold-based techniques. In [3], a three-level thresholding method by using fuzzy entropy and genetic algorithm (GA) was presented that degrades problem of threshold selection. Fuzzy -means (FCM) in conjunction with GA was discussed in [4] where the results showed a significant improvement in the accuracy of image segmentation in comparison with other methods. Another study in [5] introduced a modified FCM technique by a regularizing functional and a regularization parameter to balance clustering and smoothing in noisy and incomplete images.

In [6], a region growing vessel segmentation algorithm based on spectrum information was introduced which applied Fourier transform on the Region of Interest (ROI). ROI included vascular structures to gain spectrum information according to extracted primary features. Then combined edge information with primary feature direction computed the center points of the vascular structure as the seeds of region growing segmentation. Finally, improved region growing method with branch-based growth strategy is used to segment the vessels. In order to measure the effectiveness of this algorithm, the results on retinal and abdomen liver vascular CT scan images examined and showed that the proposed method can not only extract the high quality target vessel region, but also reduce effectively the manual intervention.

A novel 4D graph-based method to segment hepatic vasculature and tumors was introduced in [7]. The algorithm uses multiphase CT images to model the differential enhancement of the liver structures and Hessian-based shape likelihoods to avoid the common pitfalls of graph cuts under segmentation and intensity heterogeneity. Veins were tracked using the graph representation and planes fitted to the vessel segments. The method allows the detection of all hepatic tumors and identification of the liver segments with 87.8% accuracy.

Multiple automatic segmentations of liver vessels have been taken into account, for example, adaptive threshold [8], region growing [9], fuzzy -mean clustering, and level set [10]. Adaptive threshold method [8] utilized a two-phased image segmentation including conversion of gray scale to binary image and adapting the threshold value by applying a generated binary mask on the CT scan image. Region growing [9] is performed in the first and the third phases of CT scan by tracing the portal vein and the hepatic vein. Tracing the veins by a 3D labeling operation in the ROI is performed using a threshold. Threshold value helps in separating the blood vessels from the liver’s soft tissue. Erosion and dilation operations remove the adjoining stomach and spleen regions. Morphological dilation and an optimal threshold value are applied to estimate the region of the liver. The analysis showed comparable results of this method in comparison with manual detection of the regions by an expert. In [10], spatial fuzzy -means clustering combined with anatomical prior knowledge is employed to extract liver region, while a distance regularized level set is used for refinement followed by morphological operations. The experiment result shows a high accuracy (0.9986) and specificity (0.9989).

Fast FCM (FFCM) was optimized by Particle Swarm Optimization (PSO) [1] and applied on liver CT images. Evaluation performance was performed in terms of Jaccard Index and Dice Coefficient and based on ANOVA analysis showing higher values than FFCM.

High sensitivity to noisy data and being trapped in local optimum are common disadvantages of FCM and other methods as already mentioned. As the idea of using image energy [11] resolved these problems, in this study, a combination of FCM and GA is presented. Energy-based approach [11] achieved an initial segmentation closed to the liver’s boundary followed by an Active Contour Model (ACM) and GA to produce a proper parameter set closed to the optimal solution. This method had a better ability to segment the liver tissue with respect to the other methods. Furthermore, it showed high accuracy, precision, sensitivity, specificity and low overlap error, MSD, and runtime with fewer ACM iterations.

In this paper, an Innovative heuristic method is introduced based on FCM and GA. By using FCM and GA and selecting suitable features for the liver segmentation in CT images, the sensitivity to noisy data and the heavy dependency on initial data are degraded.

#### 3. Fuzzy Clustering

Fuzzy clustering can be considered as a part of fuzzy data analysis and has two parts: fuzzy data analysis and deterministic data analysis using fuzzy techniques. In the latter, each cluster is assumed as a set of elements. Then by changing the definition of membership in which “an element can only be a member of a cluster (the partition mode)” to the definition that “any element can be a member of multiple clusters with different membership degrees,” classification becomes more compatible with reality [12].

##### 3.1. Fuzzy -Means Algorithm

Fuzzy -means method is one of the most common methods involving feature analysis, clustering, and classifier design. Like the -Means method, FCM is a family member of the clustering algorithms with an objective function where all seek to minimize the function. Using fuzzy membership, FCM algorithm identifies the pixels corresponding to each group.

The first version of FCM algorithm was presented by Duda and Hart (1973) to perform an accurate clustering. Since some data were dependent on multiple clusters, it was impossible to merge them into a single cluster. This algorithm was revised several times. Finally, Bezdek proposed its final version by defining as fuzzifier [13]. The resulting algorithm identifies spherical clouds of points in a P-dimensional space. These clusters are supposed to have roughly the same size. Each cluster is displayed with a centroid. In selecting the centroid, the mean value is used as a representative of all data assigned to the cluster. To calculate the centroid, total degrees of membership of each element are divided by the product of degrees of membership raised to the power . This algorithm does not recognize clusters in different shapes, sizes, and densities. In practice, the final solution can be reached with a few iterations. Therefore, this algorithm reduces the computational time. In other types of FCM, rather than using the identity matrix in determining the distance, matrices such as diagonal might be used for segmentation of the elliptic clusters.

Assume that represents an image with pixels separated in classes, where represents data for shapes. The ISAN classic iterative optimization algorithm that reduces the cost function is defined as shown below.where represents the membership of in the th cluster, . is the center of the th class and is a constant. Parameter controls the fuzziness of the segment. The value of cost function is directly related to the Euclidean distances between the center and pixels, that is, nearness of pixels to the center reduces the cost function and increases the membership values, while pixels far from the center have low member ship values in higher costs. The membership function indicates the probability that a pixel belongs to a specified class. In the FCM algorithm, this probability only depends on the distance between the pixel and the center of each separate class in a specified area. Membership functions and class centers are modified and updated by the following formula: Standard FCM algorithm is called optimized when the values of the pixels are high near the center and have low algebra values in locations far from the center. One of the disadvantages of the standard FCM used for the image separation and segmentation is that it only uses gray-level intensity information with no use of spatial data of pixels. In fact, the probability that gray-level intensity of adjacent pixels belonged to the same class is high.

Chuang et al. proposed a spatial FCM algorithm which merges spatial data with fuzzy membership functions. The spatial function is defined as follows.where represents a square window centered on pixel in the spatial field. Similar to membership function, the spatial function represents the probability that the majority of neighboring pixels belong to the same class. The spatial function is combined with the membership function as follows. where and control the relative importance of and [10].

#### 4. Genetic Algorithm

Genetic algorithm (GA) is a well-known biological-based method taken from the genetics. The GA produces a very large set of possible solutions for a given problem. Each of these solutions is evaluated using a fitness function. Then some of the best solutions are utilized in producing new solutions as the evolved solutions. Hence, the solution is progressed until reaching an optimal one. The effectiveness of GA highly depends on proper selection of parameters [14].

Every solution for a given problem is represented in a list of parameters called chromosome. The most common method of representing chromosomes in the GA is in the form of binary strings, although other representations might be used. Initially, multiple features are randomly created as the first generation. Along with a generation, every single feature is evaluated by a fitness function and its fitness value is measured.

The next step of GA is generating the second generation according to proper selection of features in previous step. This is done by using genetic operators, for example, chromosome join and chromosome modification. Parents are selected for reproduction and are combined using crossover or mutation operators in order to produce new offspring in such a way that the best elements are chosen. Even the weakest elements have the chance of selection to avoid local answer. There are different selection methods in literature including roulette and tournament. This process is repeated until the next generation of population is produced. GA has a probability of connection between 0.6 and 1 showing the probability of generating a child. By connecting two different chromosomes a child is generated and connected to the next generation. This process continues until a good candidate (answer) is found in next generation. This process causes new generations of chromosomes which are different from previous generations. In every step, the population is studied. On the condition that the convergence criteria are met, the process is terminated. Other criteria for termination of GA are runtime and the number of generations.

#### 5. The Proposed Method

Automatic detection of liver cancer based on blood vessels consists of three main steps. They are preprocessing, image segmentation, and ROI classification.

Main problem of the segmentation algorithms is the handling of inhomogeneous or insufficient contrasted images. In line with this, the proposed method initially enhances the quality of CT image by removing noises. This process is known as preprocessing: an input CT image is passed through a Gaussian low-pass filter with relatively large kernel, the noise is eliminated, and a sharper image is obtained. Figure 1 shows the given original image and the image after preprocessing. It is noted that handling of partial volume effects in smaller vessels is not taken into account in this study.