python 简单图像处理（7） 计算直方图

好啦。现在我们来介绍下图像的直方图

说白了，直方图就是统计图像中像素点为某值的个数

对于8位灰度图来说，颜色有256级，统计每级的个数

然后把结果用图表表示出来就可以了

看看程序吧

import cv

def Hist(image):
    a = [0]*256
    w = image.width
    h = image.height
    iHist = cv.CreateImage((256,256),8,3)
    for i in range(h):
        for j in range(w):
            iGray = int(image[i,j])
            a[iGray] = a[iGray] + 1
    
    S = max(a)
    c = cv.RGB(200,150,255)
    
    for k in range(256):
        a[k] = a[k]*200/S
        x = (k,255)
        y = (k,255-a[k])
        cv.Line(iHist,x,y,c)
    
    return iHist

image = cv.LoadImage('lena.jpg',0)
iHist = Hist(image)
cv.ShowImage('image',image)
cv.ShowImage('iHist',iHist)
cv.WaitKey(0)

 

我们先统计各灰度值像素点个数

然后我们把个数值映射到（0，200）的范围

方便我们画图

我们来看看结果吧。

 

还可以啦

其实OpenCV提供了直方图函数

我简单看了下，没有丝毫思路，以后看懂了再来对比一下

 

现在我们把一幅图的三个通道放开来看看吧

 

import cv
def Hist(image,color = cv.RGB(255,255,255)):
    a = [0]*256
    w = image.width
    h = image.height
    iHist = cv.CreateImage((256,256),8,3)
    for i in range(h):
        for j in range(w):
            iGray = int(image[i,j])
            a[iGray] = a[iGray] + 1
    
    S = max(a)
    
    for k in range(256):
        a[k] = a[k]*200/S
        x = (k,255)
        y = (k,255-a[k])
        cv.Line(iHist,x,y,color)
    
    return iHist
def GetChannel(image,num):
    w = image.width
    h = image.height
    size = (w,h)
    iGetChannel = cv.CreateImage(size,image.depth,1)
    for i in range(h):
        for j in range(w):
            iGetChannel[i,j] = image[i,j][num]
    return iGetChannel
    
image = cv.LoadImage('lena.jpg',0)
image2 = cv.LoadImage('lena.jpg',1)
iHist = Hist(image)
iRed = GetChannel(image2,2)
iGreen = GetChannel(image2,1)
iBlue = GetChannel(image2,0)
iHistRed = Hist(iRed,cv.RGB(255,0,0))
iHistGreen = Hist(iGreen,cv.RGB(0,255,0))
iHistBlue = Hist(iBlue,cv.RGB(0,0,255))
cv.ShowImage('image',image)
cv.ShowImage('iRed',iRed)
cv.ShowImage('iGreen',iGreen)
cv.ShowImage('iBlue',iBlue)
cv.ShowImage('iHist',iHist)
cv.ShowImage('iHistRed',iHistRed)
cv.ShowImage('iHistGreen',iHistGreen)
cv.ShowImage('iHistBlue',iHistBlue)
cv.WaitKey(0)  

 

有一点要说明的是OpenCV中每个像素三色是按照BGR排列的

运行效果

 

 

这些直方图有什么作用呢？

以后再来说明