Skip to content

Instantly share code, notes, and snippets.

@bistaumanga
Created August 22, 2013 16:32
Show Gist options
  • Save bistaumanga/6309599 to your computer and use it in GitHub Desktop.
Save bistaumanga/6309599 to your computer and use it in GitHub Desktop.
Histogram Equalization in python
import numpy as np
def imhist(im):
# calculates normalized histogram of an image
m, n = im.shape
h = [0.0] * 256
for i in range(m):
for j in range(n):
h[im[i, j]]+=1
return np.array(h)/(m*n)
def cumsum(h):
# finds cumulative sum of a numpy array, list
return [sum(h[:i+1]) for i in range(len(h))]
def histeq(im):
#calculate Histogram
h = imhist(im)
cdf = np.array(cumsum(h)) #cumulative distribution function
sk = np.uint8(255 * cdf) #finding transfer function values
s1, s2 = im.shape
Y = np.zeros_like(im)
# applying transfered values for each pixels
for i in range(0, s1):
for j in range(0, s2):
Y[i, j] = sk[im[i, j]]
H = imhist(Y)
#return transformed image, original and new istogram,
# and transform function
return Y , h, H, sk
import pylab as plt
import matplotlib.image as mpimg
import numpy as np
# load image to numpy arrayb
# matplotlib 1.3.1 only supports png images
# use scipy or PIL for other formats
img = np.uint8(mpimg.imread('image.png')*255.0)
# convert to grayscale
# do for individual channels R, G, B, A for nongrayscale images
img = np.uint8((0.2126* img[:,:,0]) + \
np.uint8(0.7152 * img[:,:,1]) +\
np.uint8(0.0722 * img[:,:,2]))
# use hist module from hist.py to perform histogram equalization
from hist import histeq
new_img, h, new_h, sk = histeq(img)
# show old and new image
# show original image
plt.subplot(121)
plt.imshow(img)
plt.title('original image')
plt.set_cmap('gray')
# show original image
plt.subplot(122)
plt.imshow(new_img)
plt.title('hist. equalized image')
plt.set_cmap('gray')
plt.show()
# plot histograms and transfer function
fig = plt.figure()
fig.add_subplot(221)
plt.plot(h)
plt.title('Original histogram') # original histogram
fig.add_subplot(222)
plt.plot(new_h)
plt.title('New histogram') #hist of eqlauized image
fig.add_subplot(223)
plt.plot(sk)
plt.title('Transfer function') #transfer function
plt.show()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment