Appropriate Contrast Enhancement Measures for Brain and Breast Cancer Images

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.

Suneet Gupta 1 and Rabins Porwal 2

Academic Editor:Richard H. Bayford

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

Received 12 November 2015; Accepted 15 March 2016

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.

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 [figure omitted; refer to PDF] where [figure omitted; refer to PDF] and [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the target intensity and [figure omitted; refer to PDF] is the intensity of the immediate adjacent background.

2.3. Contrast Ratio

This measure is expressed in either linear or logarithmic form: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the target intensity and [figure omitted; refer to PDF] 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: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] and [figure omitted; refer to PDF] be the original and enhanced image, respectively, where [figure omitted; refer to PDF] represents the size of image. Then, MSE is defined as [figure omitted; refer to PDF]

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

Let [figure omitted; refer to PDF] be an image with [figure omitted; refer to PDF] gray levels and [figure omitted; refer to PDF] being the size. Then, PSNR is defined as [figure omitted; refer to PDF]

3.3. Absolute Mean Brightness Error (AMBE)

Let [figure omitted; refer to PDF] and [figure omitted; refer to PDF] be the original and enhanced image, respectively, where [figure omitted; refer to PDF] represents the size of image. AMBE is calculated from the following equation: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are mean gray levels of new and original image, respectively. They are given as [figure omitted; refer to PDF]

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.

Figure 1: Images with different attributes.

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

[figure omitted; refer to PDF]

(b) Semiperiodic, nonuniform background

[figure omitted; refer to PDF]

(c) Large uniform background, nonperiodic pattern

[figure omitted; refer to PDF]

(d) Periodic, nonuniform background

[figure omitted; refer to PDF]

(e) Random texture, nonperiodic, nonuniform background

[figure omitted; refer to PDF]

(f) Random texture, nonperiodic, nonuniform background

[figure omitted; refer to PDF]

(g) Nonperiodic, nonuniform background

[figure omitted; refer to PDF]

(h) Highly textured, nonperiodic, nonuniform background

[figure omitted; refer to PDF]

(i) Little texture, nonperiodic, nonuniform background

[figure omitted; refer to PDF]

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

4.1. Measure of Enhancement (EME)

The EME is defined as [figure omitted; refer to PDF] As defined in (3), the CR of a block [figure omitted; refer to PDF] is [figure omitted; refer to PDF] Therefore, [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] As defined in (3), the CR of a block [figure omitted; refer to PDF] is [figure omitted; refer to PDF] Therefore, [figure omitted; refer to PDF] Therefore, EMEE measure is the entropy of Contrast Ratio for each block [figure omitted; refer to PDF] , scaled by [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] , the more the emphasis on entropy. The term [figure omitted; refer to PDF] helps to handle more randomness.

4.3. Logarithmic Michelson Contrast Measure (AME)

The AME is defined as [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] Michelson contrast is given as [figure omitted; refer to PDF] Therefore, [figure omitted; refer to PDF] Therefore, AMEE measure is the entropy of Michelson contrast for each block [figure omitted; refer to PDF] , scaled by [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] where GLE represents the gray level energy with 256 bins and [figure omitted; refer to PDF] 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: [figure omitted; refer to PDF]

5.2. Relative Entropy (RE)

Let [figure omitted; refer to PDF] and [figure omitted; refer to PDF] be two probability functions of processing images; then, RE of [figure omitted; refer to PDF] with respect to [figure omitted; refer to PDF] is defined as [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] be an image of size [figure omitted; refer to PDF] with [figure omitted; refer to PDF] gray levels; then, element in the cooccurrence matrix [figure omitted; refer to PDF] of size [figure omitted; refer to PDF] which contains the information of correlation of the image is given as [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is element at [figure omitted; refer to PDF] th row and [figure omitted; refer to PDF] th column of matrix [figure omitted; refer to PDF] and [figure omitted; refer to PDF] The probability of cooccurrence [figure omitted; refer to PDF] is computed as [figure omitted; refer to PDF] and the SOE is computed as [figure omitted; refer to PDF]

5.4. Edge Content (EC)

Let [figure omitted; refer to PDF] be an image, where [figure omitted; refer to PDF] is the row coordinate and [figure omitted; refer to PDF] is the column coordinate. The gradient vector at any pixel position [figure omitted; refer to PDF] is computed as [figure omitted; refer to PDF] where [figure omitted; refer to PDF] Taking absolute value of image gradient as an indicator of contrast, [figure omitted; refer to PDF] The measure EC is defined as [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the size of image and [figure omitted; refer to PDF] 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, [figure omitted; refer to PDF] , [figure omitted; refer to PDF] );

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 [figure omitted; refer to PDF] for both input and output refers to [figure omitted; refer to PDF] .

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

[figure omitted; refer to PDF] = imadjust(RI, [figure omitted; refer to PDF] , [figure omitted; refer to PDF] );

[figure omitted; refer to PDF] = imadjust(RI, [figure omitted; refer to PDF] , [figure omitted; refer to PDF] );

[figure omitted; refer to PDF] = imadjust(RI, [figure omitted; refer to PDF] , [figure omitted; refer to PDF] );

[figure omitted; refer to PDF] = imadjust(RI, [figure omitted; refer to PDF] , [figure omitted; refer to PDF] );

[figure omitted; refer to PDF] = imadjust(RI, [figure omitted; refer to PDF] , [figure omitted; refer to PDF] );

where [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] are five enhanced versions of the reference image, RI. The enhancements of images are such that the contrast of [figure omitted; refer to PDF] . 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 ( [figure omitted; refer to PDF] to [figure omitted; refer to PDF] ), 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 [figure omitted; refer to PDF] 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	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
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	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
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

Figure 2: Different versions of one brain cancer image.

(a) Image-RI

[figure omitted; refer to PDF]

(b) Image- [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(c) Image- [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(d) Image- [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(e) Image- [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(f) Image- [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

Figure 3: Different versions of one breast cancer image.

(a) Image-RI

[figure omitted; refer to PDF]

(b) Image- [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(c) Image- [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(d) Image- [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(e) Image- [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(f) Image- [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

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: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are output range and input range of gray levels of the image. The denominator of (31) is 1 since [figure omitted; refer to PDF] Therefore, (31) reduces to [figure omitted; refer to PDF] 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.

[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
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 [figure omitted; refer to PDF] where [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are two variables between which the correlation, [figure omitted; refer to PDF] , is to be found and [figure omitted; refer to PDF] 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.

[figure omitted; refer to PDF] (CS)	[figure omitted; refer to PDF] (MSE)	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
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, [figure omitted; refer to PDF] or [figure omitted; refer to PDF] The value of [figure omitted; refer to PDF] ranges from -1.0 to +1.0. When the value of [figure omitted; refer to PDF] moves away from 0.0 to -1 or +1, the strength of association between the two variables increases; however, they are inversely associated if [figure omitted; refer to PDF] 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.

References

[1] Y. Zhou, K. Panetta, 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.

[2] K. Panetta, Y. Zhou, S. Agaian, H. Jia, "Nonlinear unsharp masking for mammogram enhancement," IEEE Transactions on Information Technology in Biomedicine , vol. 15, no. 6, pp. 918-928, 2011.

[3] S. S. Agaian, B. Silver, 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.

[4] K. A. Panetta, E. J. Wharton, 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.

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

[6] K. Panetta, A. Samani, S. Agaian, "Choosing the optimal spatial domain measure of enhancement for mammogram images," International Journal of Biomedical Imaging , vol. 2014, 2014.

[7] Z. Ye, H. Mohamadian, S.-S. Pang, 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.

[8] A. Saleem, A. Beghdadi, B. Boashash, "Image fusion-based contrast enhancement," EURASIP Journal on Image and Video Processing , vol. 2012, article 10, 2012.

[9] http://www.lormed.ru/images/stories/rak.jpg

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

[11] S.-M. Guo, C.-B. Li, C.-W. Chen, Y.-C. Liao, 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, pp. 1-9, 2007.

[12] N. R. Pal, S. K. Pal, "Entropic thresholding," Signal Processing , vol. 16, no. 2, pp. 97-108, 1989.

[13] R. Gonzalez, R. Woods, S. Eddins Digital Image Processing Using MATLAB , Gatesmark Publishing, 2010., 2nd.

[14] https://public.cancerimagingarchive.net/ncia/login.jsf

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