International Journal of Biomedical Imaging
Volume 2016 (2016), Article ID 4710842, 8 pages
http://dx.doi.org/10.1155/2016/4710842

Research Article

Appropriate Contrast Enhancement Measures for Brain and Breast Cancer Images

Suneet Gupta¹ and Rabins Porwal²

¹Research Scholar, Mewar University, Chittorgarh 312901, Rajasthan, India
²Department of Information Technology, Institute of Technology & Science (ITS), Ghaziabad 201007, Uttar Pradesh, India

Received 12 November 2015; Accepted 15 March 2016

Academic Editor: Richard H. Bayford

Copyright © 2016 Suneet Gupta and Rabins Porwal. 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.

Medical imaging systems often produce images that require enhancement, such as improving the image contrast as they are poor in contrast. Therefore, they must be enhanced before they are examined by medical professionals. This is necessary for proper diagnosis and subsequent treatment. We do have various enhancement algorithms which enhance the medical images to different extents. We also have various quantitative metrics or measures which evaluate the quality of an image. This paper suggests the most appropriate measures for two of the medical images, namely, brain cancer images and breast cancer images.

1. Introduction

We have different methods of enhancements which enhance the quality of an image [1, 2]. Enhancement can be done in such a way that the quality of image is upgraded visually or suiting to an application [3, 4]. Example of the latter case is medical image where the enhanced image must be more informative rather than being visually more pleasing. Therefore, to assess the quality of image, we have different measures of enhancement. There are some of the measures which suite images having different types of contents and attributes. Contrast based measures of enhancement methods, especially in spatial domain, use the intensities of pixels for calculation purpose. These spatial contrast measures are usually derived from Weber-Fechner law, Michelson contrast measure [5], or Contrast Ratio (CR). The measures PSNR, EME, AMEE [6], SOE, EC, and so forth are examples of such spatial contrast measures [7, 8]. These measures highly depend on image attributes like image content, background (uniform/nonuniform), lighting, texture, patterns which are periodic, single, or multiple targets, randomness, and also noise and distortions. Some measures or metrics are best suited to a set of attributes while other measures are best suited to other set of attributes. If a measure cannot properly tackle a particular image property, the metric or measure will not turn to prove a good measure for that image. Also, by looking at the image attributes, suitability of some measures can be suggested. In this paper, we will suggest the most appropriate or suitable measures of enhancement for brain and breast cancer images.

2. Basic Measures of Image Contrast

2.1. Michelson Contrast

This measure is generally used for images which have equivalent bright and dark features, for example, periodic patterns such as a sinusoidal grating. Michelson contrast is defined as

where and are the highest and lowest intensities.

2.2. Weber-Fechner Law

This measure is generally used for images with a large uniform background and a small test target. It is used for images whose average intensity is approximately equal to background intensity. It is given as

where is the target intensity and is the intensity of the immediate adjacent background.

2.3. Contrast Ratio

This measure is expressed in either linear or logarithmic form:

where is the target intensity and is the intensity of the immediate adjacent background.

2.4. Entropy

Entropy is calculated over entire image from the histogram of that image. It is a measure of randomness characterizing the texture of the image:

where is the histogram count for a segment of image. Increasing the contrast in one portion of an image and decreasing it in another portion may result in similar entropy as the original image. This is because the entropy is calculated over the entire image.

3. Conventional Measures of Contrast

These measures are attractive because they are simple and easy to calculate. MSE and PSNR consider corresponding pair of pixels of two images for calculation purpose. The conventional measures actually calculate the difference between two images and therefore it can be used for both restoration and enhancement. These measures are taken from [5].

3.1. Mean Squared Error (MSE)

Let and be the original and enhanced image, respectively, where represents the size of image. Then, MSE is defined as

3.2. Peak Signal-to-Noise Ratio (PSNR)

Let be an image with gray levels and being the size. Then, PSNR is defined as

3.3. Absolute Mean Brightness Error (AMBE)

Let and be the original and enhanced image, respectively, where represents the size of image. AMBE is calculated from the following equation:

where and are mean gray levels of new and original image, respectively. They are given as

4. Complex Measures of Contrast

These measures are based on the basic contrast measures and are highly sensitive to randomness, periodicity, uniform background, texture, target size, and noise. These measures are taken from [6]. These complex measures are always a good choice as good measures and, by looking at the attributes, one or two measures can be suggested. But doing it quantitatively is always a better option. Figure 1 [6, 9–11] shows images having different attributes.

(a) Large uniform background, small targets, no randomness, and no periodic patterns

(b) Semiperiodic, nonuniform background

(c) Large uniform background, nonperiodic pattern

(d) Periodic, nonuniform background

(e) Random texture, nonperiodic, nonuniform background

(f) Random texture, nonperiodic, nonuniform background

(g) Nonperiodic, nonuniform background

(h) Highly textured, nonperiodic, nonuniform background

(i) Little texture, nonperiodic, nonuniform background

Figure 1: Images with different attributes.

Let be an image which is split into blocks, with being the size of the image. Let each block be represented as , where and refer to th row and th column, respectively, inside each block. Let and be the maximum and minimum intensities in an image block .

4.1. Measure of Enhancement (EME)

The EME is defined as

As defined in (3), the CR of a block is

Therefore,

The EME measure is suitable for images with attributes like noncomplex segments, uniform background in segments, small targets in segments, nonperiodic pattern in segments, and little to no randomness in segments.

4.2. Measure of Enhancement by Entropy (EMEE)

The EMEE is defined as

As defined in (3), the CR of a block is

Therefore,

Therefore, EMEE measure is the entropy of Contrast Ratio for each block , scaled by averaged over the entire image.

The EMEE measure is suitable for images with attributes like noncomplex segments or small size complex segments, nonperiodic patterns in segments, and randomness because of the inclusion of entropy. The more the value of , the more the emphasis on entropy. The term helps to handle more randomness.

4.3. Logarithmic Michelson Contrast Measure (AME)

The AME is defined as

The AME measure is suitable for images with attributes like periodic patterns in segments and no randomness in texture. It cannot be used for images with uniform background and also images with segments of large uniform intensity.

4.4. Logarithmic AME by Entropy (AMEE)

The AMEE is defined as

Michelson contrast is given as

Therefore,

Therefore, AMEE measure is the entropy of Michelson contrast for each block , scaled by averaged over the entire image.

The AMEE measure is suitable for images with the attributes like surplus randomness in texture and periodic patterns in segments. It cannot handle images with large uniform background.

5. Other Measures of Contrast

Here, the first two measures are taken from [7] and the next two from [8].

5.1. Discrete Gray Level Energy (GLE)

The gray level energy measure basically refers to the way of distribution of gray levels. It is given as

where GLE represents the gray level energy with 256 bins and refers to the probability distribution functions of different gray levels, that is, the histogram count. In special case, GLE reaches its maximum value of 1 when the image has a constant value of gray level:

5.2. Relative Entropy (RE)

Let and be two probability functions of processing images; then, RE of with respect to is defined as

Relative entropy is also called the Kullback Leibler distance.

5.3. Second-Order Entropy (SOE)

The first-order entropy is used as enhancement metric as it is less sensitive to noise and has greater response to gray level step. But it has the drawback of ignoring the correlation among pixels. Second-order derivatives enhance noise points much more and they are more sensitive to a line than a step and to a point than to a line.

The SOE as defined in [12] uses a cooccurrence matrix to consider the spatial correlation. Let be an image of size with gray levels; then, element in the cooccurrence matrix of size which contains the information of correlation of the image is given as

where is element at th row and th column of matrix and

The probability of cooccurrence is computed as

and the SOE is computed as

5.4. Edge Content (EC)

Let be an image, where is the row coordinate and is the column coordinate. The gradient vector at any pixel position is computed as

where

Taking absolute value of image gradient as an indicator of contrast,

The measure EC is defined as

where is the size of image and

Contrast changes exist over the entire image and EC takes into account all the contrast changes even for pixels which are very close (adjacent) to each other.

6. Enhancement Procedure

To decide the most appropriate contrast measure, an enhancement procedure is required which enhances the image contrast to different levels. For this purpose, a Matlab function “imadjust()” [13] is used to increase the contrast to different levels. First, an image is read and if required converted to gray image. Let us call this image the original image. From the original image, an adjusted image of reasonably poor contrast is obtained. Let us call this image reference image. The reference image is created using the Matlab command as follows:

 Reference Image, RI = imadjust(Original_Image,  ,  );

where the second parameter, that is, the input range of gray levels, gets mapped to the third parameter, that is, the output range of gray levels. The full range of gray levels for both input and output refers to .

Now, using the reference image, RI, different images of different levels of enhancement are obtained as follows:

  = imadjust(RI, , );  = imadjust(RI, ,  );  = imadjust(RI, ,  );  = imadjust(RI, ,  );  = imadjust(RI, ,  );

where ,  ,  ,  , and  are five enhanced versions of the reference image, RI. The enhancements of images are such that the contrast of . It can be seen that for every next level of enhancement the input range has been narrowed down and the output range has been expanded resulting in contrast stretching and therefore more contrast.

7. Tests and Results

For the purpose of testing, twenty medical images are taken from [14], ten brain cancer images and ten breast cancer images. From each of these images, six versions, namely, one reference image (RI) and five enhanced images ( to ), are created as discussed above. One example of each category, namely, brain cancer and breast cancer, is shown in Figures 2 and 3. Then, the average values of the contrast enhancement measures or metrics as mentioned above are computed and the results are shown in Tables 1 and 2. As far as the complex measures EME, EMEE, AME, and AMEE are concerned, the images are divided into blocks of size 16 × 16 and in case of EMEE and AMEE the value for is taken as 0.1 as the randomness in the images of brain and breast cancer is quite low.

Table 1: Average metrics values for brain cancer images.

RI

MSE
—
8.5517
26.5507
66.3417
109.7868
184.5359

PSNR
—
38.8112
33.8924
29.9151
27.7284
25.4730

AMBE
—
2.6961
4.7250
7.4828
9.5790
12.4371

EME
4.0391
4.3321
4.5948
4.9485
5.2719
5.6941

EMEE
0.0219
0.0236
0.0251
0.0272
0.0292
0.0317

AME
69.8246
68.6713
67.8564
66.5659
65.6326
64.2347

AMEE
0.2166
0.2138
0.2118
0.2087
0.2063
0.2029

GLE
0.5199
0.5226
0.5265
0.5289
0.5331
0.5363

RE
—
12.8020
12.9135
13.0003
13.1236
13.2628

SOE
3.7928
3.7594
3.7108
3.6806
3.6304
3.5931

EC
0.9306
0.9801
1.0253
1.0805
1.1305
1.1884

Table 2: Average metrics values for breast cancer images.

RI

MSE
—
9.0064
30.4972
73.9229
125.6293
207.1439

PSNR
—
38.5890
33.2949
29.4490
27.1487
24.9767

AMBE
—
2.4062
4.3634
6.8167
8.8055
11.3444

EME
6.0779
6.4918
6.8150
7.2779
7.6364
8.1497

EMEE
0.0325
0.0349
0.0368
0.0395
0.0417
0.0448

AME
54.3720
53.5601
53.2630
52.3984
52.0707
51.1420

AMEE
0.1791
0.1766
0.1753
0.1727
0.1714
0.1687

GLE
0.2531
0.2649
0.2825
0.2934
0.3113
0.3240

RE
—
8.9075
9.3810
9.6936
10.0909
10.4761

SOE
6.5488
6.4099
6.2002
6.0690
5.8614
5.7163

EC
1.6252
1.7027
1.7643
1.8450
1.9088
1.9889

(a) Image-RI

(b) Image-

(c) Image-

(d) Image-

(e) Image-

(f) Image-

Figure 2: Different versions of one brain cancer image.

(a) Image-RI

(b) Image-

(c) Image-

(d) Image-

(e) Image-

(f) Image-

Figure 3: Different versions of one breast cancer image.

Defining contrast stretching (CS) quantitatively with respect to the reference image (RI) for the different enhanced versions of the images is as given in the following:

where and are output range and input range of gray levels of the image. The denominator of (31) is 1 since

Therefore, (31) reduces to

Thus, CS with respect to RI for each of the enhanced versions is computed and given in Table 3.

Table 3: CS for different enhanced versions.

CS
1.0690
1.1429
1.2222
1.3077
1.4

Now, to find the most appropriate measure, the values of CS are taken as standard values and correlation between CS and each of the contrast metrics is computed and all are then compared. To calculate correlation, average values of the metrics as given in Tables 1 and 2 are taken and Pearson correlation [15] is used which is given as

where and are two variables between which the correlation, , is to be found and is the number of data. The two variables here would be CS and one enhancement measure. For illustration purpose, one calculation for MSE (brain cancer) using Table 4 has been shown.

Table 4: Illustration for calculation of correlation.

(CS)
(MSE)

1.069
8.5517
9.14177
1.14276
73.1316

1.143
26.5507
30.3448
1.30622
704.940

1.222
66.342
81.0828
1.49377
4401.22

1.308
109.787
143.568
1.71008
12053.1

1.4
184.536
258.350
1.96
34053.5

6.142
395.767
522.4878
7.612834
51285.93

The last row of the table shows the summed up values of each column.

Therefore,

or

The value of ranges from −1.0 to +1.0. When the value of moves away from 0.0 to −1 or +1, the strength of association between the two variables increases; however, they are inversely associated if is negative. Table 5 shows the correlation values of the different contrast enhancement metrics used against brain and breast cancer images. As per Table 5, the most appropriate measure for brain cancer images is Edge Content (EC) and for breast cancer images EC stands as the second choice, with Absolute Mean Brightness Error (AMBE) being the first choice.

Table 5: Correlation values.

Metrics
Brain cancer images
Breast cancer images

MSE
0.9830
0.9846

PSNR
−0.9739
−0.9708

AMBE
0.9990
0.9993

EME
0.9992
0.9989

EMEE
0.9991
0.9989

AME
−0.9983
−0.9987

AMEE
−0.9986
−0.9947

GLE
0.9971
0.9960

RE
0.9987
0.9967

SOE
−0.9966
−0.9955

EC
0.9997
0.9991

8. Conclusion with Discussion

From Table 5, we can see that all measures have quite high correlation values. The conventional measures MSE, PSNR, and AMBE are more used for restoration purpose and also they do not properly reflect the perceived visual quality of the image. Therefore, if we neglect the conventional measures, EC is the most appropriate measure for both brain and breast cancer images. Regarding complex measures EME, EMEE, AME, and AMEE, as they are sensitive to image attributes, they are a good choice for any type of image in general. To find the best among complex measures, we proceed with identifying the attributes of the image. Actually, there are different attributes of images which can be best handled by different measures; one measure cannot tackle all attributes equally well. Therefore, it is important to come to a conclusion of correct set of attributes the image in question has. Then, only a particular measure can be suggested. In this paper, both brain and breast cancer images have the following properties:

(i)Large uniform background.(ii)Random texture in the foreground but no randomness in the background.(iii)Nonperiodic pattern in the foreground but absence of any pattern in the background with the possibility to be considered periodic because the whole background is dark.(iv)Background containing noncomplex segments with the foreground containing complex segments.(v)Background containing large area of uniform intensity which is absent in the foreground.

With these attributes, any recommendation would only be a compromise. Therefore, one way to make things simple is to examine only the foreground of the image as done in mammograms in [6]. Therefore, we forget the dark portion and consider the whole brain or breast as background and cells, tissues, lesions, and so forth as targets in the foreground. Now, the attributes of the image reduce to nonuniform background, random texture, and nonperiodic pattern. The best recommendation among the complex measures for these attributes is AMEE. AMEE cannot handle areas of uniform intensity well and we do have those types of areas but they are small, so we can go ahead with the recommendation.

The other aspect is when we consider the whole image, that is, when we consider the dark portion of the images as background and brain/breast as foreground, then irrespective of the attributes present in the image, the best among the complex measures can be picked up from Table 5 of this paper by seeing the correlation values. EME stands as the best choice whereas EMEE is a strong contender.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

References

Y. Zhou, K. Panetta, and S. Agaian, “Human visual system based mammogram enhancement and analysis,” in Proceedings of the 2nd International Conference on Image Processing Theory, Tools and Applications (IPTA ’10), pp. 229–234, Paris, France, July 2010.

K. Panetta, Y. Zhou, S. Agaian, and H. Jia, “Nonlinear unsharp masking for mammogram enhancement,” IEEE Transactions on Information Technology in Biomedicine, vol. 15, no. 6, pp. 918–928, 2011.

S. S. Agaian, B. Silver, and K. A. Panetta, “Transform coefficient histogram-based image enhancement algorithms using contrast entropy,” IEEE Transactions on Image Processing, vol. 16, no. 3, pp. 741–758, 2007.

K. A. Panetta, E. J. Wharton, and S. S. Agaian, “Human visual system-based image enhancement and logarithmic contrast measure,” IEEE Transactions on Systems, Man, and Cybernetics, Part B: Cybernetics, vol. 38, no. 1, pp. 174–188, 2008.

V. L. Jaya and R. Gopikakumari, “IEM: a new image enhancement metric for contrast and sharpness measurements,” International Journal of Computer Applications, vol. 79, no. 9, 2013.

K. Panetta, A. Samani, and S. Agaian, “Choosing the optimal spatial domain measure of enhancement for mammogram images,” International Journal of Biomedical Imaging, vol. 2014, Article ID 937849, 8 pages, 2014.

Z. Ye, H. Mohamadian, S.-S. Pang, and S. Iyengar, “Image contrast enhancement and quantitative measuringof information flow,” in Proceedings of the 6th International Conference on Information Security and Privacy (WSEAS '07), Tenerife, Spain, December 2007.

A. Saleem, A. Beghdadi, and B. Boashash, “Image fusion-based contrast enhancement,” EURASIP Journal on Image and Video Processing, vol. 2012, article 10, 2012.

http://www.lormed.ru/images/stories/rak.jpg

http://mojo.dailybruin.com/wp-content/uploads/2013/02/viking2mars.gif

S.-M. Guo, C.-B. Li, C.-W. Chen, Y.-C. Liao, and J. S.-H. Tsai, “Enlargement and reduction of image/video via discrete cosine transform pair, part 2: reduction,” Journal of Electronic Imaging, vol. 16, no. 4, Article ID 043007, pp. 1–9, 2007.

N. R. Pal and S. K. Pal, “Entropic thresholding,” Signal Processing, vol. 16, no. 2, pp. 97–108, 1989.

R. Gonzalez, R. Woods, and S. Eddins, Digital Image Processing Using MATLAB, Gatesmark Publishing, 2nd edition, 2010.

https://public.cancerimagingarchive.net/ncia/login.jsf

http://www.statisticshowto.com/what-is-the-pearson-correlation-coefficient/