 Open Access
 Authors : Dr. B. P. Santhosh Kumar , T. Vishnu Vardhan , Y. Naga Jagadeesh , N. Shiva Ravi Kumar, C Lakshmi Devi
 Paper ID : IJERTV9IS100062
 Volume & Issue : Volume 09, Issue 10 (October 2020)
 Published (First Online): 14102020
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
A Combinational Histogram Equalization Related to Contrast Image Enhancement
Dr. B. P. Santosh Kumar
Assistant Professor, B.Tech, M.Tech, Ph.D., MISTE, Department of ECE, YSR Engineering College of Yogi Vemana University, Proddatur, Kadapa, AP,
T. Vishnu Vardhan; Y. Naga Jagadeesh;
N.Shiva Ravi Kumar; C.Lakshmi Devi.

ech Student, Department of ECE, YSR Engineering College of Yogi Vemana University, Proddatur, Kadapa ,
AbstractThe limitation to the most commonly used histogram equalization (HE) technique is the inconsideration of the neighbourhood info near each pixel for contrast enhancement. This gives rise to noise in the output image. Image enhancement tools are often classified into point operations and spatial operations. To overcome this effect, a combinational images of histogram equalization (CHE) based technique is suggested. The focus is to utilize the information among each pixel and its neighbours, which improves the contrast of an image. This method is used in an only twodimensional images. This histogram is using the original image and its intermediate image. Further, it does not require a target uniform distribution for generating the output. The 2D cumulative distribution function (CDF) is utilized as a mapping function to get the output pixel intensity. The experimental analysis indicates that the procedure produces better results than the HE based contrast enhancement algorithms. More significantly, it produces the best results even for images having a narrow dynamic range. The implementation simplicity of the proposed algorithm may attract researchers to explore the idea for new applications in image processing.
Keywords Images; intermediate image; histogram equalization; cumulative distribution.

INTRODUCTION
Presently, highquality digital cameras are indeed the most widely used devices to acquire images. They are extensively used in cell phones, personal digital assistants, robots, medical systems and surveillance, and home security systems. Over the years, the quality of the images acquired has significantly improved due to the development of technology. Still, there are varieties of problems that need to be addressed regarding the quality of the images. Some of the problems include contrast defects, chromatic anomalies, noise, and geometrical distortions, focus defects, etc. Many image processing techniques are reported in the literature to address such problems (Gonzalez and Woods, 2009). In this paper, we are primarily concerned with image enhancement or contrast enhancement to be more specific. Image enhancement entails the modification of an image such that the output image is either more pleasing to the human eye, or contains more information and less noise, useful for further processing. Image enhancement techniques are used either as preprocessing steps or as post processing steps to generate a visually desirable image. This includes various contrast enhancement techniques to enhance the edges of the image. Image enhancement finds application in numerous areas, for instance, medical imaging, remote sensing, television, microscopic imaging, etc.
One of the main quality impairment of a digital image is low contrast. Low contrast is caused by many factors such as uneven illumination, the addition of noise during transmission, analog to digital transformations, etc. Several algorithms for contrast enhancement are proposed in the literature. Based on the information available in the literature (Tang et al., 2003), contrast enhancement algorithms are classified into two categories: transform domain and spatial domain. The transform domain algorithms partition an input image into various subbands of frequency components. The image is then enhanced by modifying the magnitude of frequency components locally or globally. Such calculations are computationally complex and timeconsuming. Moreover, to produce an output image free from distortions and visual artefacts, such algorithms require proper settings of the related parameters. The strategy accomplishes both global and local contrast improvement with an appropriate parameter choice.
In general, it is impossible to design a contrast enhancement method that produces a visual artefact free output. Choosing an appropriate contrast enhancement algorithm is difficult because of the absence of trustworthy and dependable measures to evaluate the quality of the output image. Moreover, enhancement algorithms generally depend on the legitimate parameter choice, which also experiences the absence of the tried and true measures.
The proposed method uses the correlation between the intensity of a pixel and the average intensity value of its neighbourhood to improve the contrast of the image. The correlation is achieved by building the histogram.
Fig 1: Schematic block diagram of the suggested technique.
A histogram is created by choosing a group of local pixel attributes and building a multidimensional histogram. The individual cells in the histogram matrix represent the number of pixels in the image expressed by a certain combination of attribute values. For instance, consider a histogram that combines the pixel intensity information of an image with the pixel neighbourhood average intensity information from another image.
A given pixel in an image has intensity levels
{0, L1} and its neighbourhood average intensity levels {0, L1}. The histogram will contain L x L entries. Each value corresponds to a particular pixel intensity value and its corresponding neighbourhood average intensity value at the same location. The value stored here is the number of times the intensity pair occurs in the two images. More precisely, we can develop a histogram with a given set of k attributes, where the Ith attribute has 711 possible values. A histogram is a kdimensional matrix so that each element represents the number of pixels in an image expressed by a ktuple of attribute values. The dimension of the histogram matrix is, therefore nm, the number of feasible permutations of the
y) = j. Now utilizing (i, j) and h (i, j), the histogram is formed.
B. calculate the Probability Distribution Function of an image:
The conventional approach for image enhancement is to replace individual pixel intensities with the required intensity values by forming a onedimensional (1D) histogram of the input image. Then the required intensity values are computed from the corresponding probability distribution. The 1D histogram equalization technique uses a discrete transformation function as defined below to generate the output pixel intensity.
= ( ) = ( 1) =0 (), = 0,1,2, . . 1
(2)
Where represents the equalized grey level in the output image, , is the grey level of the pixel in the input image. T () represents the transformation operator. ( ) Represents the probability of occurrence of grey level . It is to be noted that the probability of occurrence of grey level
is represented as,
( ) = 1 , = 0,1,2, . , 1 (3)
values of each attribute.
=0
A. calculating mean value of an image:
Let I represent a grey scale image of size M x N with L intensity levels G= {0, 1 L1}.Let f(x, y) be the grey value of the pixel 3 at the location (x, y) where x {1, 2 M}, y {1, 2 N}. the total number of pixels is M x N.
Let represent the average image derived from I. The size
of is the same as I with L intensity levels. Let g(x, y) be the grey value of the pixel at the corresponding location (x,
y)in . The g(x, y) is computed in a 2 neighbouring window, expressed as,
Note that this equation is not used in any of our calculations. However, it is presented here for the sake of completeness. Such an approach does not consider the information contained around each pixel while constructing the histogram. The histogram takes into account the local neighbourhood information around each pixel. In order to construct the histogram, the average intensity value around each pixel is calculated. This histogram takes into consideration the correlation existing between intensity values in a small neighbourhood of the image, which the conventional approaches ignore. The histogram is expressed
(, ) = [ 1
=
=
( + , + )] (1)
as:
Where k=/2, Note that w < min (M, N). is typically set as an odd value. Here, it is taken as three. However, researchers can use any other values of w for computing the average image. Now f(x, y) (from image I) and g(x, y) (from image ) are taken as a feature to construct the histogram. Let h (i, j) be the count of the number of times the pair (i, j) appears, where f (x, y) = i and g(x,
= {(, ) 0 1,0 1}. (4)
The term h (i, j) is the number of occurrences of the grey level pair f (x, y) and g(x, y) at the same spatial location (x, y) of the images I and respectively. It represents the
count function. As i and j can take any possible integer value between 0 and L1, the number of pixel pair combinations possible are L x L. Therefore, the histogram H will contain L x L entries.
The combined histogram for an example island image is shown in Fig. 1.
Original Image Mean image ID histogram
Histogram Matrix representation of the histogram
Fig 2: Representation of histogram concept.
Fig 2 represents the histogram concept. The original island image is taken. Its average image is obtained using (1). Here the window size is taken three. The ID histogram is shown for comparison. The histogram is formed using (4). The histogram is shown as a matrix. The horizontal axis #1 (f(x, y) = i) represents the grey levels of the original image and the horizontal axis #2 (g(x, y) = j) represents the grey levels of the average image in the 3D plot. The entries of the histogram h (i, j) (count function) show the number of times the pixel pair (i, j) appears and are represented by the vertical axis of the 3D plot.
C.Caluculate the Cumulative Distribution Function of an image:
The twodimensional (2D) cumulative distribution function (CDF) is obtained from the count function as given below:
have multiple occurrences in the original subimage. Separately for each g(x, y) = j, the original values of f(x, y) = i are replaced by the equalized values. This is further explained in the example given below. It gives the enhanced image whose dynamic range is {0 LI}. A schematic block diagram of our suggested technique is shown in Fig 2.

ALGORITHM OF THE CONTRAST HISTOGRAM EQUALISATION
The algorithm for the suggested technique consists of the following processing steps.
Step1: Consider the input image I and compute the average image using (1) where each pixel intensity value is substituted by the average intensity value of its neighbouring
=0
=0
(, ) =
=0
(, )
(5)
pixels.
Note that the computation of CDF does not depend on the size (M, N) of the images.
In this paper, we use this 2D CDF value to generate the contrast enhanced output pixel intensity. The equalized value of the intensity pairs (i, j) in the output image using the proposed method is obtained as:
Step2: Compute the histogram by comparing the input image I and the average image using (4).
Step3: The histogram count function is then used to obtain the twodimensional cumulative distribution function using (5).
Step4: The output pixel intensity is calculated using (6) and
= ( 1
1
((, ) (, )
)). (6)
the original intensity values F(x, y) =i are mapped to the
equalized ones at all occurrences of pair (i, j) only. The final
Where (, )) the minimum nonzero value of the CDF is computed using (5). Here, the size of the image
(MN) appears in the denominator so that the equalized intensity levels remain within the range {0… L1}. The equalized histogram matrix is now represented as:
= { (, ) 0 1,0 1}. (7)
It is to be noted that the dynamic range of entries of the equalized histogram matrix is extended. Then the original intensity values f(x, y) =i are replaced by (, ) at all the occurrences of i with j only. Note that f(x, y) = i may
mapping produces the output image with a dynamic range that spans a wider range of greylevel scale.
To explain the proposed method, a small subimage matrix is considered and the results are shown below in Table 1 Let an 8bit grey scale subimage I of size 6 x 6 have intensity values as shown in Table.1 (a). The average sub
image obtained using (1) is shown in Table. 1 (b).
102
123
153
196
174
43
63
96
99
100
98
64
58
109
167
211
189
65
93
143
149
152
148
96
29
87
167
240
204
94
91
141
149
153
151
98
14
72
182
247
218
123
89
139
149
152
148
96
0
131
225
145
80
21
93
144
150
147
141
92
138
255
233
51
7
36
63
98
100
95
90
58

(b)
149
Intensity level pairs
Count
CDF
(, )
Intensity level pairs
Count
CDF
(, )
I
I
112
93
1
1
0
146
144
1
19
131
122
90
1
2
7
148
63
1
20
138
130
89
1
3
14
148
147
1
21
145
132
92
1
4
21
149
99
1
22
153
133
91
1
5
29
149
2
24
167
134
58
1
6
36
150
98
1
25
174
135
64
1
7
43
150
149
1
26
182
136
95
1
8
51
152
148
1
27
189
137
93
1
9
58
154
100
1
28
196
139
96
1
10
65
154
151
1
29
204
140
139
1
11
72
155
152
1
30
211
140
141
1
12
80
156
148
1
31
218
141
141
1
13
87
156
150
1
32
225
142
98
1
14
94
158
100
1
33
233
143
63
1
15
102
158
153
1
34
240
143
143
1
16
109
160
152
1
35
247
145
96
2
18
123
164
98
1
36
255
(c)
102
123
153
196
174
43
58
109
167
211
189
65
29
87
167
240
204
94
14
72
182
247
218
123
0
131
225
145
80
21
138
255
233
51
7
36
(d)
Table 1: Example subimage to demonstrate histogram equalization. (a) Sample sub image matrix,

Average sub image matrix, (c) histogram equalization, (d) Equalized sub image matrix.
An example cameraman grey scale image and its contrast enhanced results using different methods are presented in table 1. The cameraman image shown has bright and dark areas. The background contains a light coloured sky and building. The unequalised histogram confirms the concentration of intensity values near the origin. The output
images found with HE, MWCVMHE and FHSABP methods look similar. However, the sky and cameraman's face looks degraded. HMF method gives a better contrast with minor deformations on the sky area and somehow darkens the whole image.
Original histogram HE
MWCVMHE FHSABP HMF
CEBGA 2DHE CHE
Fig 3: Results of contrast enhancement for cameraman image using different methods.
Note that the sample images used for comparison are represented as follows.
Image 1 Plane image Image 5Cessna image
Image 2 Tank image Image 6 Lighthouse image
Image 3Cameraman image Image 7_ Beach image
Image 4 Baboon image Image 8 Island image
The theoretical time complexity of the proposed algorithm is computed as in Cao et al. (2018). The
the absolute difference between the average values of the input and the output images as Wang and Ward (2007).
computation time required to construct the joint histogram for a 'bbit' grey scale image of size M x N is O (2MN). The
(, ) = 1 1+()()
(8)
mapping function using two dimensional CDF requires time O (2b ). The mapping requires O (MN) time. Hence, the total time complexity is O (3MN + 2b ). It is to be noted that the time complexity of the proposed technique may be slightly higher than the conventional HE technique. However, it is very effective while enhancing narrow dynamic range images and reduces the noise in the enhanced image.
Existing parameters:
The objective measurement utilized to evaluate the performance of the approaches in preserving the original brightness of the input image is AMBE. This is expressed as
Here, I and O represent the input and output image,
respectively. Note that E (.) represents the statistical mean value. A higher AMBE value shows improved brightness preservation. Discrete Entropy is a statistical quantity of randomness which is utilized to describe the textural characteristics of the input image. It determines the content in an image. A higher value indicates richer details.
=0
=0
=0
=0
It is expressed as () = 1 ( ) log ( ) where ( ) is the probability of pixel intensity computed from the normalized histogram of the input image I. Note that () = 1 ( ) log ( ) where
( )the probability of pixel intensity is computed from
the normalized histogram of the output image O. Then the normalized DE for the input image I and the corresponding output image O is defined as
The concept is employed to develop the suggested technique. The suggested joint histogram equalization scheme utilizes the intensity distribution surrounding each
1
log(256)()
=
=
1+ log(256)()
(9)
Table 2
pixel in an image to improvise the contrast. It is implemented for both the grey and coloured images.
Comparison of AMBE for different methods.
/table>
Table 3
Comparison of for different methods.
Image
HE
MWCVMHE
FHSABP
HMF
CEBGA
2DHE
CHE
Im N1
Im N2
Im N3
Im N4
Im N5
Im N6
Im N7
Im N8
Average
0.0260
0.4928
0.0944
0.3618
0.0197
0.1141
0.0198
0.2495
0.1723
0.0824
0.1039
0.1471
0.4936
0.2926
0.7862
0.1603
0.5470
0.3266
0.1246
0.2140
0.5470
0.7338
0.2308
0.5283
0.6743
0.5625
0.4519
0.2328
0.0661
0.2601
0.0938
0.0392
0.1587
0.0467
0.0932
0.1238
0.0330
0.0239
0.0782
0.0564
0.2225
0.4376
0.0371
0.1524
0.1301
0.3269
0.0360
0.0473
0.1276
0.5973
0.0973
0.0245
0.1873
0.1805
0.3291
0.2639
0.1653
0.2487
0.6108
0.2099
0.1712
0.2461
0.2806
Image
HE
MWCVMHE
FHSABP
HMF
CEBGA
2DHE
CHE
Im N1 Im N2 Im N3
Irn N4 Im N5 Im N6 Im N7 Im N8
Average
0.4920
0.4880
0.4458
0.4572
0.4623
0.4502
0.4528
0.4532
0.4627
0.4909
0.4867
0.4516
0.4630
0.4690
0.4597
0.4630
0.4573
0.4677
0.4973
0.4895
0.4432
0.4588
0.4558
0.4538
0.4501
0.4517
0.4625
0.4935
0.4945
0.4585
0.4785
0.4522
0.4847
0.4619
0.4843
0.4760
0.4975
0.4525
0.3876
0.3197
0.4560
0.3625
0.4389
0.3787
0.4117
0.4990
0.4951
0.4812
0.4802
0.4815
0.4866
0.4767
0.4670
0.4830
0.5117
0.5211
0.4951
0.4992
0.4923
0.4961
0.4931
0.5009
0.5012
House image Beach image Enhanced image Enhanced image Fig. 4: Results of contrast enhancement for narrow dynamic range images using the JHE method.


CONCLUSION:

In this work, a truly two dimensional (domain) concept is employed to develop the suggested technique. The suggested combinational histogram equalization
scheme utilizes the intensity distribution surrounding each
pixel in an image to improvise the contrast. It is implemented for both the grey and coloured images. The experimental analysis shows that the algorithm produces better or
competitive results with respect to most of the stateofthe art algorithms. The histogram equalization produces the best results for images having a narrow dynamic range. Different quantitative assessment metrics are used to validate the algorithm.
The histogram concept used in the proposed method may pave the way for investigation of more sophisticated contrast enhancement algorithms. Further, this method can also be used along with other techniques to achieve both global and local contrast enhancement. The suggested technique is very simple but useful for contrast enhancement. It does not require a target uniform distribution and uses the CDF as the mapping function to generate the output image.
REFERENCE:

Arici, T., Dikbas, S., Altunbasak, Y., 2009. A histogram modification framework and its application for image contrast enhancement. IEEE Trans. Image Process. 18 (9), 19211935.

Beghdadi, A, Negrate, L, 1989. Contrast enhancement technique based on local detection of edges. Computer Vision, Graphics, Image Process. 46 (2), 162174.

Cao, G., Tian, H., Yu, L, Huang, X., Wang, Y., 2018.Acceleration of histogrambased contrast enhancement via selective down sampling. IET Image Proc. 12 (3), 447452.

Celik, T., 2012. Twodimensional histogram equalization and contrast enhancement. Pattern Recognition 45 (10), 38103824.

Celik, T., 2016. Spatial mutual information and PageRankbased contrast enhancement and qualityaware relative contrast measure. IEEE Trans. Image Process 25 (10), 47194728.

Chen, S.D., Ramli, A, 2003. Minimum mean brightness error bi histogram equalization in contrast enhancement. IEEE Trans. Consum. Electron. 49 (4), 13101319.

Chen, S.D., Ramli, AR., 2003. Contrast enhancement using recursive mean separate histogram equalization for scalable brightness preservation. IEEE Trans. Consum. Electron. 49 (4), 13011309.

DaleJones, R., Tjahjadi, T., 1993. A study and modification of the local histogram equalization algorithm. Pattern Recogn. 26 (9), 13731381.

Gonzalez, R.C., Woods, R.E., 2009. Digital Image Processing. Pearson, New Delhi, India.

Hashemi, S., Kiani, S., Noroozi, N., Moghaddam, M.E., 2010. An image contrast enhancement method based on genetic.
ABOUT AUTHORS:
Dr. B. P. Santosh Kumar
Presently working as Assistant Professor, Department of ECE, YSR Engineering College of Yogi Vemana University, Proddatur, India. He received the B.Tech. Degree from Jawaharlal Nehru Technological University, Hyderabad, India and the M.Tech. Degree from Kerala University, Thiruvananthapuram, India. He received the Ph.D. degree from Yogi Vemana University, Kadapa, India. His current research interests include image processing.
T.Vishnu Vardhan
Presently pursuing bachelor of engineering degree in Electronics and Communication Engineering at Y.S.R.Engineering College of Yogi Vemana University proddutur, AP, India.
Y. Naga Jagadeesh
Presently pursuing bachelor of engineering degree in Electronics and Communication Engineering at Y.S.R.Engineering College of Yogi Vemana University proddutur, AP, India.
N. Shiva Ravi Kumar
Presently pursuing bachelor of engineering degree in Electronics and Communication Engineering at Y.S.R.Engineering College of Yogi Vemana University proddutur, AP, India.
C. Lakshmi Devi
Presently pursuing bachelor of engineering degree in Electronics and Communication Engineering at Y.S.R.Engineering College of Yogi Vemana University proddutur, AP, India.