Improvement of Fuzzy Image Contrast Enhancement Using Simulated Ergodic Fuzzy Markov Chains

This paper presents a novel fuzzy enhancement technique using simulated ergodic fuzzy Markov chains for low contrast brain magnetic resonance imaging (MRI). The fuzzy image contrast enhancement is proposed by weighted fuzzy expected value. The membership values are then modified to enhance the image using ergodic fuzzyMarkov chains.The qualitative performance of the proposedmethod is compared to another method in which ergodic fuzzyMarkov chains are not considered.The proposedmethod produces better quality image.


Introduction
Fuzzy set theory is useful in managing various uncertainties in computer vision and image processing applications.Fuzzy image processing is a collection of different fuzzy methods to image processing that can comprehend, characterize, and process the images.It has two main phases, namely, image fuzzification and modification of membership values.
Image enhancement for brightness transformation valid in many practical cases can be position brightness correction and gray-scale transformation.Position brightness correction adjusts the pixel brightness by taking into account the pixel position in the image.Gray-scale transformation changes the pixel brightness, but it does not take into account the position of the pixels in the image.Gray-scale transformation is just a transformation of the gray scale to another scale to increase the contrast.The purpose of the transformation is to improve the visual appearance of an image.
There are many research works on image enhancement [1,2], but in this paper we focus on ergodic fuzzy Markov chains for image enhancement.
Avrachenkov and Sanchez [3] introduced fuzzy Markov chains with a transition possibility measure and a general state space.Also, Kalenatic et al. [4] presented a simulation study on fuzzy Markov chains to identify some characteristics about their behavior, based on matrix analysis.All the aforementioned investigations show that fuzzy Markov chains have a periodic behavior.We improved behavior of fuzzy Markov chains using Halton sequences and simulated ergodic fuzzy Markov chains [5].In this paper, we apply our technique of simulating ergodic fuzzy Markov chains for generating membership values of pixels.Enhancement using ergodic fuzzy Markov chains will improve the quality of the image and provide a clear image to the human observer.
The overall approach of the paper follows.We consider a low quality fuzzy image with  ×  pixels (  ;  = 0, 1 ⋅ ⋅ ⋅  − 1,  = 0, 1 ⋅ ⋅ ⋅  − 1).We then obtain the value of threshold  of the fuzzy image.The main and novel idea to enhance the image contrast is to consider a pixel as a new image and subdivide it to  ×  pixels.To simulate these  ×  pixels we employ the ergodic fuzzy Markov chains and their transition matrix.For a particular pixel, its related  ×  matrix entries represent membership values of gray levels of the pixel.To increase the image contrast of each pixel using subdivided pixel we use the weighted fuzzy expected value approach based on transition matrix entries of ergodic fuzzy Markov chain.We show that proposed method produces better quality images compared to the fuzzy expected value method given in [1].
In Section 2, we discuss the fuzzy image contrast enhancement.In Section 3, we define the similarity measure and distance measure, and then in Section 4 we review ergodic fuzzy Markov chains.The main result of this paper is presented in Section 5 followed by a simulation study in Section 6. Section 7 shows the performance of our new results based on fuzzy image contrast enhancement using simulated ergodic fuzzy Markov chains.Comparison of our result with original MR image is also presented.

Image as a Fuzzy Set, Threshold Technique, and Fuzzy Image Contrast Enhancement
The application of fuzzy set theory in image processing took formal shape only in the 1980s with the pioneering research carried out by Pal et al. [6] and Pal and Rosenfeld [2].
2.1.Image as a Fuzzy Set.The pixel values, which establish an image, may not be accurate and there is basically an intrinsic imprecision or uncertainty embedded in a digital image.While trying to design automated systems for scene analysis and explanation, it may be a good idea to consider the fact that a computer vision system is usually embedded with uncertainty and vagueness, which needs to be taken care of suitably.Proper modeling of this imprecision appearing in a physical phenomenon is an important mission in many applications of image processing.
(1) Equation ( 1) interprets the characteristics of an image with  ×  pixels.The double summations in (1) just refer to a collection of pixels and their membership values not a crisp mathematical summation.

Fuzzy Threshold Technique.
Thresholding is an operation that involves tests against a function  of the form where(, ) is the gray level of pixel (, ) and (, ) denotes some local property of this point.A threshold image (, ) is defined as Thus, pixels label 1 (or any other gray level) correspond to objects.However pixels label 0 (or any other gray level not assigned to objects) correspond to the background.When  depends only on (, ) (that is, only on gray level values) the threshold is called global.If  depends on both (, ) and (, ), the threshold is called local.If, in addition,  depends on the spatial coordinates  and , the threshold is called dynamic or adaptive [7].
Figure 1 captures an MR image, which is not transparent.Its histogram drawn in MATLAB by using pixels properties and their gray levels is also given in this figure.Figures 2 and 3 demonstrate tumor and abdomen images and their related histograms, respectively.As shown in the histograms, in all the cases the ratio of bright pixels to dark ones is low [8, 9].

Fuzzy Image Contrast
Enhancement.Fuzzy image contrast enhancement is established on gray level mapping from a gray plane into a fuzzy plane using a membership value.It uses the principle of contrast stretching where the image gray levels are transformed in such a way that dark pixels appear much darker and bright pixels appear much brighter.The principle of contrast stretching depends on the selection of threshold, , so that the gray levels below the threshold  are reduced and the gray levels above the threshold  are increased in a nonlinear manner.This stretching operation induces saturation at both ends (gray levels): CONT() is an intensifier stretching that transforms a gray level  to an intensified level.,  are the levels between 0 and 1 that decide the percentage stretching of gray level  for a certain threshold "." Fuzzy methods for contrast enhancement employ membership values () to know the degree of brightness or darkness of the pixels in an image.So, ergodic fuzzy Markov chains are used to find the membership values of the pixels in an image that lies in the interval [0, 1] using any membership value.Then, elements of transition possibility matrix of ergodic fuzzy Markov chains are applied on the membership values to generate new membership values of the pixels in the image [1].

Fuzzy Similarity and Distance Measure
In this section, fuzzy distance measure is described using fuzzy set theory.Fuzzy distance measure is studied by many authors [10,11].We consider a universal set  and () to be a fuzzy set.Let () be the class of all crisp sets of .Also, consider two fuzzy sets,  and , such that ,  ∈ ().The properties of distance measure between two fuzzy sets  and  with membership values   () and   () are given as follows [10].

Distance Measure.
A function  : () 2 → [0, ∞] is called a distance measure if it satisfies the following properties.
For example, consider two fuzzy sets  and  in a finite  = { 1 ,  2 , . . .,   } with   () and   () as the membership value of sets  and , respectively.The distance measures are as follow.
Hamming Distance.The Hamming distance between two fuzzy sets  and  is given as Euclidean Distance.The Euclidean distance between two fuzzy sets  and  is given as
We note that it does not need elements of each row of the matrix  to sum up to one.This fuzzy matrix  allows defining all relations among the  states of the fuzzy Markov chain at each time instant , as follows [3].
At each instant ,  = 1, 2 ⋅ ⋅ ⋅ , the state of system is described by the fuzzy set  () .The transition law of a fuzzy Markov chain is given by the fuzzy relational matrix  at instant ,  = 1, 2 ⋅ ⋅ ⋅ , as follows: where  and , ,  = 1, 2, . . ., , are the initial and final states of the transition and  (0) is the initial distribution.Also, Thomason in [12] shows that the powers of a fuzzy matrix are stable over the max-min operator.More information about powers of a fuzzy matrix can be found in [13].Now, a stationary distribution of a fuzzy matrix is defined as follows.
Definition 1.Let the powers of the fuzzy transition matrix  converge in  steps to a nonperiodic solution; then the associated fuzzy Markov chain is called aperiodic fuzzy Markov chain and  * =   is its stationary fuzzy transition matrix.Definition 2. A fuzzy Markov chain is called strong ergodic if it is aperiodic and its stationary transition matrix has identical rows.
A fuzzy Markov chain is called weakly ergodic if it is aperiodic and its stationary transition matrix is stable with no identical rows.
We simulated ergodic fuzzy Markov chains using Halton sequences [5] and make use of them to simulate fuzzy image enhancement.Elements of transition possibility matrix of ergodic fuzzy Markov chains are applied on the membership values to generate new membership values of the pixels in the image.For example  45 = 0.5 in transition possibility matrix of ergodic fuzzy Markov chains corresponding to image means that the membership value of gray pixel (4, 5) is 0.6.
Note.Every pixel   in ( 1) is considered again as an image with  ×  pixels and is simulated by using transition matrix of ergodic fuzzy Markov chains.For examples the matrix of a pixel is shown as follows.

Fuzzy Image Contrast Enhancement Using Ergodic Fuzzy Markov Chains
In recent years, many researchers have applied various fuzzy methods for contrast enhancement [1].In this paper, we discuss fuzzy expected value method.Fuzzy expected value replaces the mean and median value when treating fuzzy sets.Instead of calculating the average value of a set of numbers, we evaluated a more representative value of a set.This value would indicate a typical grade of membership of a fuzzy set.Consider a fuzzy set  in a finite set  = { 1 ,  2 . . .,   } with a membership value   : [0, 1].
Suppose a variational sampling ( ( 1 , 2 ,...,  ) ( 1 , 2 ,...,  ) ) is given,   =   (  ),  = 1, 2, . . .,  are the membership values of some fuzzy set  ⊂ , () is a nonnegative monotonically decreasing function defined over the interval [0, 1], and  > 1 is a real number.Consider the following equation with respect to : The solution of ( 11) is called the weighted fuzzy expected value (WFEV) of order  with the attached weight function  of membership values ( 1 ,  2 . . .  ).The parameter  measures the dependence of frequencies of population pixels on the WFEV.Now, we apply WFEV and simulated ergodic fuzzy Markov chains corresponding to image.Each element of transition possibility matrix of ergodic fuzzy Markov chains is   .
Suppose that the WFEV for (|  − |) =  −|  −| ,  = 2,  0 = FEV is given by the following form: where  = 0, 1, 2, . . .,  − 1 is the gray levels.The equation for "" is in the form  = () and is solved iteratively starting with  +1 = (  )( +1 = (  )), where  0 = FEV.This WFEV is used for image enhancement, and when applied on the gray level represents the most typical gray level.For   , we use ergodic fuzzy Markov chains, which in the last simulated by Halton sequences in [5] and use their fuzzy transition matrix corresponding to each pixel of the image.We note that elements of fuzzy transition possibility matrix are membership values of gray levels.

Algorithm.
The algorithm of our approach is as follows.
(1) Read the original image.
(3) Represent each pixel   as a  ×  matrix.
(4) Simulate the membership values, that is, matrix entries   , by using ergodic fuzzy Markov chains implemented in [5].

Simulation
Consider a low contrast image.We apply FEV and WFEV methods to improve the image contrast.These methods are presented in the last section.In our proposed approach, considering each pixel as a new image and simulating membership values of pixel using the transition matrix of ergodic fuzzy Markov chain and applying WFEV method, we would enhance the image contrast.Figure 5 and its related histogram depict the contrast enhancement of the image using FEV method.Figure 6 also shows the same improvement in image contrast using WFEV and ergodic fuzzy Markov chains.As shown in histograms, the number of dark pixels in Figure 5 is less than those in Figure 2, and the number of dark pixels in Figure 6 is less than those in Figure 5, exhibiting that using ergodic fuzzy Markov chain along with employing WFEV method is a superior approach compared to the others.The same inference results from Figures 7 and 8 as well as Figures 9 and 10.We have employed MATLAB software to implement this approach.

Conclusion
In this paper, using simulated ergodic fuzzy Markov chains given in [5] and WFEV method we not only enhance the (MRI) low image in Figures 6, 8, and 10 but we also obtain better results than those resulted from FEV approach without considering ergodic fuzzy Markov chains for   .This can be seen from simulation study of Section 6.As we subdivide pixels to smaller pixels and treating each new smaller pixel as an image and employing the simulation phase for the image, we would improve the quality of the original image.We hope to employ ergodic fuzzy Markov chains for contrast improvement using an intensification operator (INT approach).

Figure 2 :Figure 3 :
Figure 2: Tumor image and its histogram of gray levels.

Figure 5 :Figure 6 :Figure 7 :
Figure 5: Contrast enhancement of brain MRI using FEV and its histogram.

Figure 8 :Figure 9 :Figure 10 :
Figure 8: Contrast enhancement of tumor image using WFEV and ergodic fuzzy Markov chains and its histogram.

Table 1
ergodic fuzzy Markov chains, which gives most typical value of the membership value , where weights are applied.The weight is calculated as follows.