# MRI Brain Images Healthy and Pathological Tissues Classification with the Aid of Improved Particle Swarm Optimization and Neural Network

Authors: V. Sheejakumari; B. Sankara Gomathi
Journal: Computational and Mathematical Methods in Medicine
 (2015)

---

## Abstract

The advantages of magnetic resonance imaging (MRI) over other diagnostic imaging modalities are its higher spatial resolution and its better discrimination of soft tissue. In the previous tissues classification method, the healthy and pathological tissues are classified from the MRI brain images using HGANN. But the method lacks sensitivity and accuracy measures. The classification method is inadequate in its performance in terms of these two parameters. So, to avoid these drawbacks, a new classification method is proposed in this paper. Here, new tissues classification method is proposed with improved particle swarm optimization (IPSO) technique to classify the healthy and pathological tissues from the given MRI images.  Our proposed classification method includes the same four stages, namely, tissue segmentation, feature extraction, heuristic feature selection, and tissue classification. The method is implemented and the results are analyzed in terms of various statistical performance measures. The results show the effectiveness of the proposed classification method in classifying the tissues and the achieved improvement in sensitivity and accuracy measures. Furthermore, the performance of the proposed technique is evaluated by comparing it with the other segmentation methods.


---

## Body

## 1. Introduction

Normally, in brain tissue segmentation on magnetic resonance (MR) images, the type of tissue present for each pixel or voxel in a 2D or 3D dataset, respectively, is determined based on the information captured from both MR images and prior knowledge of the brain. Segmentation of brain tissue in MRI is a crucial preprocessing step in several medical research and clinical applications, involving measurement of tissue volume, visualization, and analysis of anatomical structures, multimodality fusion and registration, functional brain mapping, pathology recognition, surgical planning and navigation, and brain substructure segmentation [1]. The investigations of medical images for computer-aided diagnosis and treatment often require segmentation as an initial stage. Medical image segmentation is a difficult and challenging task mainly due to the vague nature of the images [2–4]. Magnetic resonance imaging (MRI) is one of the most significant diagnostic imaging techniques, often used for the early detection of anomalous changes in tissues and organs [5, 6], and also it allows a radiologist to produce an image covering the internal features of living tissue because it is a noninvasive imaging technique [7].It is well-known that the brain has a complicated structure; thus, accurate segmentation of brain is very decisive for detecting tumors, edema, necrotic tissues, white matter, gray matter, cerebrospinal fluid (CSF), or vasculature in order to provide proper treatment [8]. A technique to segment tissues into these categories is a vital step in quantitative morphology of brain because most brain structures are defined by boundaries of these tissue classes [9]. Unlike other diagnostic methods, magnetic resonance imaging (MRI) systems create numerous images, where diverse fundamental parameters of internal anatomical structures in the same body section are highlighted by each image with multiple contrasts, based on local variations of spin-spin relaxation time (T2), spin-lattice relaxation time (T1), and proton density (PD) [10]. The major obstruction to segmentation of MR images is the occurrence of noise, flaws in the scanners, and the structural variations of the imaging objects, which can be classified into four types, namely, thermal/electronic noise, magnetic field in homogeneities, biological tissue variations, and partial volume effects [11].Furthermore, manual detection and analysis of lesions from MR brain images are normally time-consuming and expensive and can produce unacceptably high intraobserver and interobserver variability [12]. The efficacy of segmented MR images in the medical diagnostic process depends on the combination of two, often conflicting, requirements, that is, the removal of superfluous information present in the original MR images and the maintenance of significant information in the resulting segmented images [13, 14]. MR image segmentation methods are often analyzed in terms of their potentiality to discriminate (i) between cerebrospinal fluid (CSF), white matter, and gray matter and (ii) between normal tissues and abnormalities [15]. Recently, numerous techniques have been proposed for the segmentation of brain tissues in MR image. Some of them are classical pattern recognition techniques, rule-based systems, image analysis methods, crisp and fuzzy clustering procedures, feed-forward neural networks, fuzzy reasoning [16], geometric models to specify lesion boundaries, connected component analysis, deterministic annealing, atlas based techniques, and contouring approaches [17, 18].


## 2. Related Works

A few recent works related to this existing in the literature are reviewed in the following section.Wang and Chen [19] have proposed a classification technique called vector seeded region growing (VSRG), in which the seed pixel vectors have been selected through standard deviation and relative Euclidean distance. By the VSRG processing, the data dimensionality of MRI has been reduced. A series of experiments has been carried out and compared to the normally used c-means technique for performance evaluation. The results have exhibited the efficacy of the proposed technique in MR image classification.An automatic approach for the segmentation of anatomical 3D brain MR images has been proposed by Cherradi et al. [20] as well. The proposed technique includes many important steps. Initially, noise reduction has been done via median filtering. Secondly, segmentation of brain/nonbrain tissue has been carried out using a threshold morphologic brain extraction technique (TMBE). Subsequently, initial centroids estimation by gray level histogram analysis has been executed, and this stage yields to a modified version of fuzzy C-means algorithm (MFCM) that has been employed for MRI tissue segmentation. At last, 3D visualization of three clusters such as CSF, GM, and WM has been performed. The competency of the proposed technique has been validated by conducting wide segmentation experiments using simulated and real MR images. A confrontation of the technique with similar methods of the literature has been undertaken through diverse performance measures. The MFCM for tissue segmentation has achieved a gain in rapidity of convergence of about 70%.Rajendran and Dhanasekaran [21] have investigated the segmentation of MRI brain image into different tissue types on brain image via possibilistic fuzzy c-means (PFCM) clustering algorithm. Application of this technique to MRI brain image has produced better segmentation result than the fuzzy c-mean (FCM) and fuzzy possibilistic c-means (FPCM) algorithms. The outcomes have been verified quantitatively through similarity metrics, false positive volume function (FPVF), and false negative volume functions (FNVF). These values have proved that the proposed PFCM has segmented the tumor class successfully by employing the membership and possibility (typicality) functions.Mehta et al. [22] have presented an approach to automatically generate fuzzy rules for tissue classification in MRI. The proposed scheme was based on hybrid approach of two popular genetic algorithm based machine learning (GBML) techniques, Michigan and Pittsburg approach. The proposed approach has utilized a training dataset generated from manual segmented images with the aid of an expert in MRI. Features from image histogram and spatial neighborhood of pixels have been employed in fuzzy rules. The approach has been tested for classifying brain T2 weighted 2D axial images obtained by different pulse sequences into three main tissue types, namely, white matter (WM), gray matter (GM), and cerebrospinal fluid (CSF). Also, the experts have matched the results of proposed approach with manual segmentation. Moreover, it has been found that the performance of proposed approach was comparable.Hussain et al. [23] have devised a technique for an accurate segmentation of normal and pathological tissues in the MRI brain images. In the proposed segmentation technique, initially the classification process has been done using fuzzy inference system (FIS) and FFBNN. Both classifiers have used the extracted image features as an input for the classification process. The features have been extorted in two ways from the MRI brain images. The FIS perform the classification process by generating the fuzzy rules using the extracted features. Five features have been extracted from the MRI images: two dynamic statistical features and three 2D wavelet decomposition features. In segmentation, the normal tissues such as WM, GM, and CSF have been segmented from the normal MRI images and pathological tissues such as edema and tumor have been segmented from the anomalous images. In the preprocessing stage, the noncortical tissues in the normal images have been removed. The implementation results have exhibited the potency of the proposed tissue segmentation technique in segmenting the tissues precisely from the MRI images. Moreover, the performance of the segmentation technique has been analyzed using performance measures such as accuracy, specificity, and sensitivity and compared with K-means clustering and fuzzy ANN based segmentation techniques.In the previous tissues classification method, the healthy and pathological tissues are classified from the MRI brain images using HGANN (Figure2) [24]. When compared to other tissues classification methods, this method performance lacks the accuracy and sensitivity measures. The classification method is inadequate in its performance in terms of these two parameters. So, to avoid these drawbacks, a new classification method with improved PSO technique is proposed in this paper. The outline of the paper is as follows. The proposed ISPONN classification process is briefly explained in Section 3. In Section 3.1, tissue segmentation process is performed and feature extraction process is described in Section 3.2. Sections 3.3 and 3.4 explain the heuristic feature selection by ISPO and tissue classification process. The experimental result and conclusion of this paper are given in Sections 4 and 5.


## 3. The Proposed IPSONN Technique for MRI Brain Tissues Classification

In this paper, we proposed an efficient method to classify the normal and pathological tissues in the MRI brain images with IPSO (improved particle swarm optimization) and FFBNN. Four major stages are involved in our proposed methodology:(i)
tissue segmentation;
(ii)
feature extraction;
(iii)
heuristic features selection by IPSO;
(iv)
tissue classification by FFBNN.


### 3.1. Tissue Segmentation

In training images, the tissues are segmented manually by MRI experts. Before the segmentation process, the input images are given to the preprocessing stage. In preprocessing, skull stripping method, described in [23], is applied to the input images in order to remove the distinct dark ring surrounding the brain tissues. After the skull stripping, the brain image tissues are segmented. The normal image obtained after skull stripping is denoted by I
s.


#### 3.1.1. Normal Tissue Segmentation

In MRI brain images, the normal tissues such as WM, GM, and CSF are segmented. The detail description about this segmentation process is given in [23]. The major step is to segment the WM and GM tissues from the image I
s by utilizing gradient method. The smoothing process is performed in the input image I
s by applying Gaussian convolution filter. After that, gradient operation is applied to the image I
G. The gradient of two variables x and y is defined as follows:(1)
∇
I
G
(
x
,
y
)
=
∂
I
G
∂
x
i
∧
+
∂
I
G
∂
y
j
∧
.Then, the binarization process is performed in the edge marked imageE
m in [23]. The MRI brain image WM and GM tissues are segmented based on their intensity values [23]. The detailed process is given in [24]. Consider(2)
I
w
g
=
WM
;
if
I
b
i
=
1
GM
;
if
I
b
i
=
0
.The results of WM and GM segmented images are denoted byI
w, I
g. Another one normal tissue CSF is segmented by the orthogonal polynomial transform (OPT). In orthogonal polynomial transformation, image I
s is computed using the following formula:(3)
I
c
f
=
Sin
I
s
(
i
)
3
100
2
+
0.05
∗
rand
I
s
.In (3), rand
(
) represents a random value. After the polynomial transform, the corresponding CSF region is segmented in the resultant image I
c
f [23].


#### 3.1.2. Abnormal Tissue Segmentation

Edema tissue is segmented from the abnormal imageI
s. Before the edema segmentation process, histogram equalization process is executed over the image I
s. The quality of image I
s is enhanced by the histogram equalization and it is denoted by I
s
′. Then, the enhanced I
s
′ image is converted into indexed image by using multilevel thresholding function. Each pixel in the image is compared with these threshold values to select the pixels(4)
H
⟶
t
3
,
S
⟶
t
4
,
V
⟶
t
5
X
=
p
u
;
t
4
≤
p
u
≤
t
3
,
t
5
0
;
otherwise
.Subsequently, the distance is determined between the coordinates of the center pixels of the regions inX
h
(
c
)
(
x
,
y
) and the tumor centroid coordinate value t
(
x
,
y
)
(5)
O
h
x
,
y
=
X
h
c
x
,
y
-
t
x
,
y
.The resultantO
h
(
x
,
y
) is then verified with threshold value t
6 and edema region coordinate values are obtained(6)
I
e
=
O
h
(
x
,
y
)
≥
t
6
0
;
otherwise
.Then, the morphological dilation and closing operations are performed in imageI
e.The tumor tissue segmentation is performed in the abnormal brain MRI images. The main objective is to segment the tumor tissue in the abnormal imageI
s. Here, we utilize the region growing method (RGM) to segment the tumor tissue. Region growing method is a region based image segmentation method; it selects the initial seed points from the input image I
s. The RGM observes the neighbor pixel values with the initial seed points, which checks whether the neighbor pixels are included in this region or not. The tumor segmentation result is represented as I
t.


### 3.2. Feature Extraction

Training patterns are generated by extracting efficient features from the training MRI images. To accomplish the feature extraction in training and testing phase, the segmented images are utilized. The extracted features from the segmented images are discussed in Section3.2.1.


#### 3.2.1. Feature Analysis

In this paper, we have utilized seven features. Among these seven features, two features are histogram based, two are statistical features, and the other three features are from wavelet. The individual performance of these features in the existing methods has both advantages and disadvantages. Thus, the combination of such features makes our proposed tissue segmentation method more efficient. All features are extracted from these nonzero blocks and computed the mean value for all feature values from each block(7)
I
i
c
=
∑
k
=
1
N
I
(
b
k
)
,
S
i
c
=
∑
k
=
1
N
S
(
b
k
)
,
M
i
c
=
∑
k
=
1
N
M
(
b
k
)
,
E
i
c
=
∑
k
=
1
N
E
(
b
k
)
,
H
i
c
=
∑
k
=
1
N
H
(
b
k
)
,
V
i
c
=
∑
k
=
1
N
V
(
b
k
)
,
D
i
c
=
∑
k
=
1
N
D
(
b
k
)
.The extracted features from the segmented images are represented as(8)
F
=
I
i
c
,
S
i
c
,
M
i
c
,
E
i
c
,
H
i
c
,
V
i
c
,
D
i
c
,
i
=
1,2
,
…
I
,
c
=
1,2
,
…
C
,where i is the number of segmented images and c is the segmented tissue classes. In (8), I denotes histogram intensity, S is the slope along the histogram, M is the mean value, E is the variance value, H denotes horizontal band, V denotes vertical band, and D is the diagonal band.


### 3.3. Heuristic Feature Selection by Improved Particle Swarm Optimization

Here, we utilize an improved PSO (IPSO) technique for improving the accuracy in classification process. The proposed IPSO technique considered both the bad and the good experience components during the new velocity computation. The bad experience component helps the particle to remember its previously visited worst position. To calculate the new velocity, the bad experience of the particle is also taken into consideration. On including the characteristics ofP best and P worst in the velocity updating process along with the difference between the present best particle and current particle, respectively, the convergence towards the solution is found to be faster and an optimal solution is reached in comparison with conventional PSO approaches. Hence, our proposed IPSO selects the most optimal features from more number of features in the feature extraction process. More number of features in the training phase increase the time complexity and reduce the classification accuracy. To shun these drawbacks during the training process, we select the most optimal features from more number of available features. The selected optimal features from the ISPO technique are given to the FFBNN.The process of optimal feature selection by IPSO technique is described as follows (Figure1).Figure 1
Flowchart of IPSO.Figure 2
Proposed HGANN tissue classificationc-FFBNN structure.(i) Initialization. Initially the particles are generated for the class c. The defined particles are composed of the features, represented as G
l
=
{
I
c
,
S
c
,
M
c
,
E
c
,
H
c
,
V
c
,
D
c
}, c
=
1,2
,
…
C.(ii) Parameters. In IPSO, the particles position, velocity, learning parameters, inertia, weight, and maximum number of iterations are defined.(iii) Fitness. Every particle’s fitness value is calculated by using the formula which is given in (9). The particles that have the maximum fitness value are selected as the best particles (9)
f
c
l
a
=
I
c
l
-
I
c
a
2
+
S
c
l
-
S
c
a
2
+
M
c
l
-
M
c
a
n
f
-
1
+
E
c
l
-
E
c
a
+
H
c
l
-
H
c
a
I
c
l
-
I
c
a
2
+
V
c
l
-
V
c
a
+
D
c
l
-
D
c
a
·
n
f
-
1
1
/
2
.In (9),I denotes histogram intensity;
S is the slope along the histogram;
M is the mean value;
E is the variance value;
H denotes the horizontal band;
V denotes the vertical band;
D is the diagonal band;
l is particle in the class c;
a
′ is another particle in the class c;
(
f
c
)
l
a is the best particles having the maximum fitness value.The best particle denotes the maximum fitness value which means all features in the class c have high value. So, the tissue type is more accurately classified by selecting the maximum fitness value particles.Select particles individual best value and particles global best value for each generation. In IPSO, we select the particles individual worst value, that is, the particle too away from the target.(iv) Velocity and Position. Update particle individual best (P best), global best (g best), and particle worst (P worst) in the velocity formula which is given in the following equation and obtain the new velocity:(10)
V
i
=
w
∗
V
i
+
C
1
b
∗
r
1
∗
P
best
i
-
S
i
∗
P
best
i
+
C
1
w
∗
r
2
∗
S
i
-
P
worst
i
∗
P
worst
i
+
C
2
∗
r
3
∗
G
best
i
-
S
i
.In (10), w is the inertia weight;
V
i is the velocity of the particle;
C
1
b is the acceleration coefficient in best position;
C
1
w is acceleration coefficient in worst position;
P
best
i is the best position of the particle i;
S
i is the current position of the particle;
P
worst
i is the worst position of the particle i;
r
1, r
2, r
3 are the uniformly distributed random numbers in the range [0 to 1].Thus, the obtained new velocity value is updated in the original velocity formula given in (11). We obtain the position of the particle by the following: (11)
V
i
=
w
∗
V
i
+
C
1
∗
r
1
∗
P
best
i
-
S
i
+
C
2
∗
r
2
∗
G
best
i
-
S
i
,
(12)
S
i
+
1
=
S
i
+
V
i
.(v) Stopping Criteria. The process is repeated until the maximum number of iterations is reached. The final optimal feature set from IPSO technique is exploited in the tissues classification process.


### 3.4. Tissue Classification by FFBNN

The feature setF
c is given to the FFBNN classifier for training process and this classifier is represented as c-FFBNN. In the training phase, the selected heuristic features are given to the c-FFBNN network. The c-FFBNN network is well trained by these selected features. The network is created with seven input units from feature set F
c, that is, F
c
=
{
I
c
,
S
c
,
M
c
,
E
c
,
H
c
,
V
c
,
D
c
}, H
d hidden units, and one output unit, that is, c. The basic structure of FFBNN network is shown in Figure 3. This FFBNN is used for classifying the input feature set that belongs to which tissue class type.Figure 3
Sample input MRI brain images.The following steps portray the function of the neural network.Step 1.
Assign input weights to all neurons except the neurons in the input layer. The planned bias function and activation function for the neural network are described as follows:(13)
Z
c
=
β
+
∑
n
′
=
0
H
d
-
1
w
c
n
′
I
n
′
c
+
w
c
n
′
S
n
′
c
+
w
c
n
′
M
n
′
c
+
w
c
n
′
E
n
′
c
+
w
c
n
′
H
n
′
c
+
w
c
n
′
V
n
′
c
+
w
c
n
′
D
n
′
c
,
δ
X
c
=
1
1
+
e
-
Z
c
.Equation (10) represents the input layer bias function. Here, (
I
n
′
c
), (
S
n
′
c
), (
M
n
′
c
), (
E
n
′
c
), (
H
n
′
c
), (
V
n
′
c
), and (
D
n
′
c
) are the features of the class c. Equation (11) represents the activation function (δ) for the output layer.Step 2.
Calculate the learning error for the neural network(14)
L
(
e
)
=
1
H
d
∑
n
′
=
0
H
d
-
1
D
n
′
-
Z
n
′
.In (14), L
(
e
) is the learning error rate and D
n
′ and Z
n
′ are the desired and actual outputs, respectively. The error that occurred during the FFBNN training process is minimized by the back propagation algorithm. The error minimization process by back propagation algorithm is briefly explained in [24]. The c-FFBNN network is well trained by these selected features and this network classifies the particular feature values based on the tissue types it belongs or not. Similar FFBNN training process is carried out for the all feature set training.In testing phase, a new set of testing images are taken for the classification. The testing images normal and abnormal tissues are segmented by the method, which is described in Section3.1. From the segmented images, the features M, V, VB, DB, HB, Hist-Inten, and Hist-Slope values are calculated and these extracted feature values are given as input to the well-trained c-FFBNN classifiers. Based on the input features, the specific tissue type c-FFBNN provides any one of the tissue types such as WM, GM, CSF, tumor, or edema as output. In this manner, the MRI brain image tissues are classified.


## 3.1. Tissue Segmentation

In training images, the tissues are segmented manually by MRI experts. Before the segmentation process, the input images are given to the preprocessing stage. In preprocessing, skull stripping method, described in [23], is applied to the input images in order to remove the distinct dark ring surrounding the brain tissues. After the skull stripping, the brain image tissues are segmented. The normal image obtained after skull stripping is denoted by I
s.


### 3.1.1. Normal Tissue Segmentation

In MRI brain images, the normal tissues such as WM, GM, and CSF are segmented. The detail description about this segmentation process is given in [23]. The major step is to segment the WM and GM tissues from the image I
s by utilizing gradient method. The smoothing process is performed in the input image I
s by applying Gaussian convolution filter. After that, gradient operation is applied to the image I
G. The gradient of two variables x and y is defined as follows:(1)
∇
I
G
(
x
,
y
)
=
∂
I
G
∂
x
i
∧
+
∂
I
G
∂
y
j
∧
.Then, the binarization process is performed in the edge marked imageE
m in [23]. The MRI brain image WM and GM tissues are segmented based on their intensity values [23]. The detailed process is given in [24]. Consider(2)
I
w
g
=
WM
;
if
I
b
i
=
1
GM
;
if
I
b
i
=
0
.The results of WM and GM segmented images are denoted byI
w, I
g. Another one normal tissue CSF is segmented by the orthogonal polynomial transform (OPT). In orthogonal polynomial transformation, image I
s is computed using the following formula:(3)
I
c
f
=
Sin
I
s
(
i
)
3
100
2
+
0.05
∗
rand
I
s
.In (3), rand
(
) represents a random value. After the polynomial transform, the corresponding CSF region is segmented in the resultant image I
c
f [23].


### 3.1.2. Abnormal Tissue Segmentation

Edema tissue is segmented from the abnormal imageI
s. Before the edema segmentation process, histogram equalization process is executed over the image I
s. The quality of image I
s is enhanced by the histogram equalization and it is denoted by I
s
′. Then, the enhanced I
s
′ image is converted into indexed image by using multilevel thresholding function. Each pixel in the image is compared with these threshold values to select the pixels(4)
H
⟶
t
3
,
S
⟶
t
4
,
V
⟶
t
5
X
=
p
u
;
t
4
≤
p
u
≤
t
3
,
t
5
0
;
otherwise
.Subsequently, the distance is determined between the coordinates of the center pixels of the regions inX
h
(
c
)
(
x
,
y
) and the tumor centroid coordinate value t
(
x
,
y
)
(5)
O
h
x
,
y
=
X
h
c
x
,
y
-
t
x
,
y
.The resultantO
h
(
x
,
y
) is then verified with threshold value t
6 and edema region coordinate values are obtained(6)
I
e
=
O
h
(
x
,
y
)
≥
t
6
0
;
otherwise
.Then, the morphological dilation and closing operations are performed in imageI
e.The tumor tissue segmentation is performed in the abnormal brain MRI images. The main objective is to segment the tumor tissue in the abnormal imageI
s. Here, we utilize the region growing method (RGM) to segment the tumor tissue. Region growing method is a region based image segmentation method; it selects the initial seed points from the input image I
s. The RGM observes the neighbor pixel values with the initial seed points, which checks whether the neighbor pixels are included in this region or not. The tumor segmentation result is represented as I
t.


## 3.1.1. Normal Tissue Segmentation

In MRI brain images, the normal tissues such as WM, GM, and CSF are segmented. The detail description about this segmentation process is given in [23]. The major step is to segment the WM and GM tissues from the image I
s by utilizing gradient method. The smoothing process is performed in the input image I
s by applying Gaussian convolution filter. After that, gradient operation is applied to the image I
G. The gradient of two variables x and y is defined as follows:(1)
∇
I
G
(
x
,
y
)
=
∂
I
G
∂
x
i
∧
+
∂
I
G
∂
y
j
∧
.Then, the binarization process is performed in the edge marked imageE
m in [23]. The MRI brain image WM and GM tissues are segmented based on their intensity values [23]. The detailed process is given in [24]. Consider(2)
I
w
g
=
WM
;
if
I
b
i
=
1
GM
;
if
I
b
i
=
0
.The results of WM and GM segmented images are denoted byI
w, I
g. Another one normal tissue CSF is segmented by the orthogonal polynomial transform (OPT). In orthogonal polynomial transformation, image I
s is computed using the following formula:(3)
I
c
f
=
Sin
I
s
(
i
)
3
100
2
+
0.05
∗
rand
I
s
.In (3), rand
(
) represents a random value. After the polynomial transform, the corresponding CSF region is segmented in the resultant image I
c
f [23].


## 3.1.2. Abnormal Tissue Segmentation

Edema tissue is segmented from the abnormal imageI
s. Before the edema segmentation process, histogram equalization process is executed over the image I
s. The quality of image I
s is enhanced by the histogram equalization and it is denoted by I
s
′. Then, the enhanced I
s
′ image is converted into indexed image by using multilevel thresholding function. Each pixel in the image is compared with these threshold values to select the pixels(4)
H
⟶
t
3
,
S
⟶
t
4
,
V
⟶
t
5
X
=
p
u
;
t
4
≤
p
u
≤
t
3
,
t
5
0
;
otherwise
.Subsequently, the distance is determined between the coordinates of the center pixels of the regions inX
h
(
c
)
(
x
,
y
) and the tumor centroid coordinate value t
(
x
,
y
)
(5)
O
h
x
,
y
=
X
h
c
x
,
y
-
t
x
,
y
.The resultantO
h
(
x
,
y
) is then verified with threshold value t
6 and edema region coordinate values are obtained(6)
I
e
=
O
h
(
x
,
y
)
≥
t
6
0
;
otherwise
.Then, the morphological dilation and closing operations are performed in imageI
e.The tumor tissue segmentation is performed in the abnormal brain MRI images. The main objective is to segment the tumor tissue in the abnormal imageI
s. Here, we utilize the region growing method (RGM) to segment the tumor tissue. Region growing method is a region based image segmentation method; it selects the initial seed points from the input image I
s. The RGM observes the neighbor pixel values with the initial seed points, which checks whether the neighbor pixels are included in this region or not. The tumor segmentation result is represented as I
t.


## 3.2. Feature Extraction

Training patterns are generated by extracting efficient features from the training MRI images. To accomplish the feature extraction in training and testing phase, the segmented images are utilized. The extracted features from the segmented images are discussed in Section3.2.1.


### 3.2.1. Feature Analysis

In this paper, we have utilized seven features. Among these seven features, two features are histogram based, two are statistical features, and the other three features are from wavelet. The individual performance of these features in the existing methods has both advantages and disadvantages. Thus, the combination of such features makes our proposed tissue segmentation method more efficient. All features are extracted from these nonzero blocks and computed the mean value for all feature values from each block(7)
I
i
c
=
∑
k
=
1
N
I
(
b
k
)
,
S
i
c
=
∑
k
=
1
N
S
(
b
k
)
,
M
i
c
=
∑
k
=
1
N
M
(
b
k
)
,
E
i
c
=
∑
k
=
1
N
E
(
b
k
)
,
H
i
c
=
∑
k
=
1
N
H
(
b
k
)
,
V
i
c
=
∑
k
=
1
N
V
(
b
k
)
,
D
i
c
=
∑
k
=
1
N
D
(
b
k
)
.The extracted features from the segmented images are represented as(8)
F
=
I
i
c
,
S
i
c
,
M
i
c
,
E
i
c
,
H
i
c
,
V
i
c
,
D
i
c
,
i
=
1,2
,
…
I
,
c
=
1,2
,
…
C
,where i is the number of segmented images and c is the segmented tissue classes. In (8), I denotes histogram intensity, S is the slope along the histogram, M is the mean value, E is the variance value, H denotes horizontal band, V denotes vertical band, and D is the diagonal band.


## 3.2.1. Feature Analysis

In this paper, we have utilized seven features. Among these seven features, two features are histogram based, two are statistical features, and the other three features are from wavelet. The individual performance of these features in the existing methods has both advantages and disadvantages. Thus, the combination of such features makes our proposed tissue segmentation method more efficient. All features are extracted from these nonzero blocks and computed the mean value for all feature values from each block(7)
I
i
c
=
∑
k
=
1
N
I
(
b
k
)
,
S
i
c
=
∑
k
=
1
N
S
(
b
k
)
,
M
i
c
=
∑
k
=
1
N
M
(
b
k
)
,
E
i
c
=
∑
k
=
1
N
E
(
b
k
)
,
H
i
c
=
∑
k
=
1
N
H
(
b
k
)
,
V
i
c
=
∑
k
=
1
N
V
(
b
k
)
,
D
i
c
=
∑
k
=
1
N
D
(
b
k
)
.The extracted features from the segmented images are represented as(8)
F
=
I
i
c
,
S
i
c
,
M
i
c
,
E
i
c
,
H
i
c
,
V
i
c
,
D
i
c
,
i
=
1,2
,
…
I
,
c
=
1,2
,
…
C
,where i is the number of segmented images and c is the segmented tissue classes. In (8), I denotes histogram intensity, S is the slope along the histogram, M is the mean value, E is the variance value, H denotes horizontal band, V denotes vertical band, and D is the diagonal band.


## 3.3. Heuristic Feature Selection by Improved Particle Swarm Optimization

Here, we utilize an improved PSO (IPSO) technique for improving the accuracy in classification process. The proposed IPSO technique considered both the bad and the good experience components during the new velocity computation. The bad experience component helps the particle to remember its previously visited worst position. To calculate the new velocity, the bad experience of the particle is also taken into consideration. On including the characteristics ofP best and P worst in the velocity updating process along with the difference between the present best particle and current particle, respectively, the convergence towards the solution is found to be faster and an optimal solution is reached in comparison with conventional PSO approaches. Hence, our proposed IPSO selects the most optimal features from more number of features in the feature extraction process. More number of features in the training phase increase the time complexity and reduce the classification accuracy. To shun these drawbacks during the training process, we select the most optimal features from more number of available features. The selected optimal features from the ISPO technique are given to the FFBNN.The process of optimal feature selection by IPSO technique is described as follows (Figure1).Figure 1
Flowchart of IPSO.Figure 2
Proposed HGANN tissue classificationc-FFBNN structure.(i) Initialization. Initially the particles are generated for the class c. The defined particles are composed of the features, represented as G
l
=
{
I
c
,
S
c
,
M
c
,
E
c
,
H
c
,
V
c
,
D
c
}, c
=
1,2
,
…
C.(ii) Parameters. In IPSO, the particles position, velocity, learning parameters, inertia, weight, and maximum number of iterations are defined.(iii) Fitness. Every particle’s fitness value is calculated by using the formula which is given in (9). The particles that have the maximum fitness value are selected as the best particles (9)
f
c
l
a
=
I
c
l
-
I
c
a
2
+
S
c
l
-
S
c
a
2
+
M
c
l
-
M
c
a
n
f
-
1
+
E
c
l
-
E
c
a
+
H
c
l
-
H
c
a
I
c
l
-
I
c
a
2
+
V
c
l
-
V
c
a
+
D
c
l
-
D
c
a
·
n
f
-
1
1
/
2
.In (9),I denotes histogram intensity;
S is the slope along the histogram;
M is the mean value;
E is the variance value;
H denotes the horizontal band;
V denotes the vertical band;
D is the diagonal band;
l is particle in the class c;
a
′ is another particle in the class c;
(
f
c
)
l
a is the best particles having the maximum fitness value.The best particle denotes the maximum fitness value which means all features in the class c have high value. So, the tissue type is more accurately classified by selecting the maximum fitness value particles.Select particles individual best value and particles global best value for each generation. In IPSO, we select the particles individual worst value, that is, the particle too away from the target.(iv) Velocity and Position. Update particle individual best (P best), global best (g best), and particle worst (P worst) in the velocity formula which is given in the following equation and obtain the new velocity:(10)
V
i
=
w
∗
V
i
+
C
1
b
∗
r
1
∗
P
best
i
-
S
i
∗
P
best
i
+
C
1
w
∗
r
2
∗
S
i
-
P
worst
i
∗
P
worst
i
+
C
2
∗
r
3
∗
G
best
i
-
S
i
.In (10), w is the inertia weight;
V
i is the velocity of the particle;
C
1
b is the acceleration coefficient in best position;
C
1
w is acceleration coefficient in worst position;
P
best
i is the best position of the particle i;
S
i is the current position of the particle;
P
worst
i is the worst position of the particle i;
r
1, r
2, r
3 are the uniformly distributed random numbers in the range [0 to 1].Thus, the obtained new velocity value is updated in the original velocity formula given in (11). We obtain the position of the particle by the following: (11)
V
i
=
w
∗
V
i
+
C
1
∗
r
1
∗
P
best
i
-
S
i
+
C
2
∗
r
2
∗
G
best
i
-
S
i
,
(12)
S
i
+
1
=
S
i
+
V
i
.(v) Stopping Criteria. The process is repeated until the maximum number of iterations is reached. The final optimal feature set from IPSO technique is exploited in the tissues classification process.


## 3.4. Tissue Classification by FFBNN

The feature setF
c is given to the FFBNN classifier for training process and this classifier is represented as c-FFBNN. In the training phase, the selected heuristic features are given to the c-FFBNN network. The c-FFBNN network is well trained by these selected features. The network is created with seven input units from feature set F
c, that is, F
c
=
{
I
c
,
S
c
,
M
c
,
E
c
,
H
c
,
V
c
,
D
c
}, H
d hidden units, and one output unit, that is, c. The basic structure of FFBNN network is shown in Figure 3. This FFBNN is used for classifying the input feature set that belongs to which tissue class type.Figure 3
Sample input MRI brain images.The following steps portray the function of the neural network.Step 1.
Assign input weights to all neurons except the neurons in the input layer. The planned bias function and activation function for the neural network are described as follows:(13)
Z
c
=
β
+
∑
n
′
=
0
H
d
-
1
w
c
n
′
I
n
′
c
+
w
c
n
′
S
n
′
c
+
w
c
n
′
M
n
′
c
+
w
c
n
′
E
n
′
c
+
w
c
n
′
H
n
′
c
+
w
c
n
′
V
n
′
c
+
w
c
n
′
D
n
′
c
,
δ
X
c
=
1
1
+
e
-
Z
c
.Equation (10) represents the input layer bias function. Here, (
I
n
′
c
), (
S
n
′
c
), (
M
n
′
c
), (
E
n
′
c
), (
H
n
′
c
), (
V
n
′
c
), and (
D
n
′
c
) are the features of the class c. Equation (11) represents the activation function (δ) for the output layer.Step 2.
Calculate the learning error for the neural network(14)
L
(
e
)
=
1
H
d
∑
n
′
=
0
H
d
-
1
D
n
′
-
Z
n
′
.In (14), L
(
e
) is the learning error rate and D
n
′ and Z
n
′ are the desired and actual outputs, respectively. The error that occurred during the FFBNN training process is minimized by the back propagation algorithm. The error minimization process by back propagation algorithm is briefly explained in [24]. The c-FFBNN network is well trained by these selected features and this network classifies the particular feature values based on the tissue types it belongs or not. Similar FFBNN training process is carried out for the all feature set training.In testing phase, a new set of testing images are taken for the classification. The testing images normal and abnormal tissues are segmented by the method, which is described in Section3.1. From the segmented images, the features M, V, VB, DB, HB, Hist-Inten, and Hist-Slope values are calculated and these extracted feature values are given as input to the well-trained c-FFBNN classifiers. Based on the input features, the specific tissue type c-FFBNN provides any one of the tissue types such as WM, GM, CSF, tumor, or edema as output. In this manner, the MRI brain image tissues are classified.


## 4. Experimental Results and Discussion

The proposed MRI brain tissues classification technique is implemented in the working platform of MATLAB (version 7.12) with machine configuration as follows:processor: Intel core i5;
OS: Windows 7;
CPU speed: 3.20 GHz;
RAM: 4 GB.The brain MRI images are collected from various medical diagnosis centers; thus, the collected MRI sample images and segmented results are shown in Figures3 and 4.Figure 4
Manually segmented MRI brain tissues: (a) WM, (b) GM, (c) CSF, (d) tumor, and (e) edema.
(a)
(b)
(c)
(d)
(e)From these manually segmented images, seven features such as histogram intensity, slope of histogram, mean, variance, horizontal, vertical, and diagonal functions of 2D wavelet decomposition are extracted, and, using these extracted features, the heuristic features are generated for sample 5 images by IPSO given in Table1.Table 1
Heuristic features from IPSO.
Heuristic features from IPSO
Images
Histogramintensity
Slope ofhistogram
Mean
Variance
Horizontal value of2D wavelet decomposition
Vertical value of 2Dwavelet decomposition
Diagonal value of 2Dwavelet decomposition
1
2
255
1.049225
28.47425
22.17121
8.785642
11.08208
2
2
255
2.263733
66.02611
31.9345
13.95081
6.356591
3
2
255
1.886505
48.81538
23.90411
21.98224
20.34238
4
2
255
1.629013
38.28219
12.74081
18.08271
18.13534
5
2
255
2.083221
63.93179
27.44711
26.96236
19.40522Utilizing these heuristic features, the FFBNN are well trained. In testing process, the testing image’s normal and abnormal tissues are segmented by the methods in [23]. The segmented image results are shown in Figure 5.Figure 5
Segmented MRI brain tissues: (a) WM, (b) GM, (c) CSF, (d) tumor, and (e) edema.
(a)
(b)
(c)
(d)
(e)Performance Analysis. The performance of proposed tissue classification method is analyzed by the statistical measures. The normal and abnormal tissues classification accuracy is calculated by these statistical measures, which are shown in Table 2.Table 2
Performance of the proposed improved particle swarm optimization and neural network (IPSONN)classification method in classifying (a) WM, (b) GM, (c) CSF, (d) edema, and (e) tumor from five different brain MRI images.
(a)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
1
0
4
0
100
0.0
100
100
100
100
0
44.7
2
0
1
4
0
0
20.0
80
80
0
100
100
0
3
1
0
4
0
100
0.0
100
100
100
100
0
44.7
(b)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
1
1
3
0
100
25.0
80
75
50
100
50
43.3
2
1
0
4
0
100
0.0
100
100
100
100
0
44.7
3
1
0
4
0
100
0.0
100
100
100
100
0
44.7
(c)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
1
0
4
0
100
0.0
100
100
100
100
0
44.7
2
1
0
4
0
100
0.0
100
100
100
100
0
44.7
3
0
1
4
0
100
0.0
100
100
100
100
0
44.7
(d)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
1
0
4
0
100
0.0
100
100
100
100
0
44.7
2
1
1
3
0
100
25.0
80
75
50
100
50
43.3
3
1
0
4
0
100
0.0
100
100
100
100
0
44.7
(e)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
0
1
4
0
0
20.0
80
80
0
100
100
0
2
1
0
4
0
100
0.0
100
100
100
100
0
44.7
3
1
0
4
0
100
0.0
100
100
100
100
0
44.7As can be seen from Table2, there are three images utilized in the performance analysis. The classification of normal WM, GM, and CSF tissues has given 93%, 87%, and 93% mean accuracy results, respectively. Thus, this higher accuracy has offered more precise classification in the normal images. Furthermore, the pathological tissues such as edema and tumor have also given 87% and 93% mean accuracy results, respectively. Hence, our proposed tissue classification method has offered more efficient and effective results in both normal and pathological tissues classification processes.Comparative Analysis. The classification performance of proposed IPSONN is analyzed with the existing tissue classification methods [22, 24]. This existing tissue classification method has utilized FRBS and hybrid genetic approach. The existing classification methods results are shown in Tables 3 and 4.Table 3
Performance measures of the existing tissue classification method [22] in classifying (a) WM, (b) GM, (c) CSF, (d) edema, and (e) tumor from five MRI brain images.
(a)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
0
1
3
1
0
25.0
60
75
0
75
100
-
14.4
2
1
2
2
0
100
50.0
60
50
33
100
67
40.8
3
0
1
2
2
0
33.3
40
67
0
50
100
-
28.9
(b)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
1
1
1
2
33
50.0
40
50
50
33
50
-
16.7
2
0
2
3
0
0
40.0
60
60
0
100
100
0
3
1
2
1
1
50
66.7
40
33
33
50
67
-
20.4
(c)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
0
3
2
0
0
60.0
40
40
0
100
100
0
2
1
2
1
2
33
66.7
33
33
33
33
67
-
40.8
3
1
4
0
0
100
100.0
20
0
20
0
80
0
(d)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
0
3
2
0
0
60.0
40
40
0
100
100
0
2
0
3
2
0
0
60.0
40
40
0
100
100
0
3
1
3
1
0
100
75.0
40
25
25
100
75
35.4
(e)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
1
2
1
1
50
66.7
40
33
33
50
67
-
20.4
2
0
3
1
1
0
75.0
20
25
0
50
100
-
106.1
3
0
2
3
0
0
40.0
60
60
0
100
100
0Table 4
Performance measures of the existing HGNN tissue classification method in classifying (a) WM, (b) GM, (c) CSF, (d) edema, and (e) tumor from five MRI brain images.
(a)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
1
0
4
0
100
0.0
100
100
100
100
0
44.7
2
1
0
4
0
100
0.0
100
100
100
100
0
44.7
3
0
1
4
0
0
20.0
80
80
0
100
100
0
(b)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
1
1
3
0
100
25.0
80
75
50
100
50
43.3
2
1
0
3
1
50
0.0
80
100
100
75
0
30.6
3
1
0
4
0
100
0.0
100
100
100
100
0
44.7
(c)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
1
0
4
0
100
0.0
100
100
100
100
0
44.7
2
1
1
3
0
100
25.0
80
75
50
100
50
43.3
3
1
0
4
0
100
0.0
100
100
100
100
0
44.7
(d)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
1
1
3
0
100
25.0
80
75
50
100
50
43.3
2
1
0
4
0
100
0.0
100
100
100
100
0
44.7
3
1
1
3
0
100
25.0
80
75
50
100
50
43.3
(e)
Images
TP
FP
TN
FN
Sensitivity
FPR
ACC
Specificity
PPV
NPV
FDR
MCC
1
1
0
3
1
50
0.0
80
100
100
75
0
30.6
2
1
0
4
0
100
0.0
100
100
100
100
0
44.7
3
1
0
4
0
100
0.0
100
100
100
100
0
44.7The graphical representation of the proposed and existing techniques average performance in tissue classification process is shown in Figure6.Figure 6
Tissues classification result: (a) accuracy, (b) sensitivity, and (c) specificity.
(a)
(b)
(c)Figure6 shows the graphical representation of WM, GM, CSF, tumor, and edema tissue classification performance compared to the existing HGNN and hybrid genetic method. As can be seen from Figure 6(a), the existing hybrid approach has given lower accuracy than the methods HGNN and IPSONN. When compared to HGNN and IPSONN, our proposed IPSONN has given higher accuracy in GM, CSF, and edema than HGNN, but the other two tissues accuracy results are the same for both methods. The tissues classification methods IPSONN, HGNN, and the existing hybrid genetic approach attain an overall mean accuracy of 95%, 91%, and 42%, respectively. Our proposed IPSONN tissues classification method has given higher accuracy result than the existing methods. In Figures 6(b) and 6(c), HGNN attains high sensitivity and specificity in tumor tissues classification, but the other tissues classification our proposed method maintains high and same sensitivity and specificity level than the HGNN and Existing hybrid approach. In sensitivity performance review, our proposed IPSONN and HGNN attain 87% and the existing hybrid genetic approach conquers only 31% sensitivity. In specificity measure, the tissues classification methods like IPSONN, HGNN, and hybrid genetic approach achieve 94%, 92%, and 42% specificity, respectively. Hence, our proposed IPSONN has given high performance in tissues classification than the HGNN and the existing hybrid approach.


## 5. Conclusion

In this paper, we proposed a new tissue classification method called IPSONN to classify the normal and abnormal tissues from the MRI images. The method was implemented and a smaller set of MRI brain images were utilized to analyze the results of the IPSONN classification method. In future works, the results presented in this paper should be extended and verified using larger and diverse image datasets. The performance analysis proved that the IPSONN method offers almost 95% accuracy and 87% and 94% (in average) sensitivity and specificity measures, respectively. The IPSONN method was compared against the existing tissues classification methods to prove the performance. Thus, the results show that the IPSONN achieved more classification performance than the existing hybrid genetic and HGNN methods.



---

*Source: 807826-2015-04-22.xml*