直方图均匀化与直方图的比较
- cv2.equalizeHist函数原型:equalizeHist(src[, dst]) 。函数equalizeHist的作用:直方图均衡化,提高图像质量。
- 直方图均衡化:如果一副图像的像素占有很多的灰度级而且分布均匀,那么这样的图像往往有高对比度和多变的灰度色调。直方图均衡化就是一种能仅靠输入图像直方图信息自动达到这种效果的变换函数。它的基本思想是对图像中像素个数多的灰度级进行展宽,而对图像中像素个数少的灰度进行压缩,从而扩展像元取值的动态范围,提高了对比度和灰度色调的变化,使图像更加清晰。
- 全局直方图均衡化可能得到是一种全局意义上的均衡化,但是有的时候这种操作并不是很好,会把某些不该调整的部分给调整了。Opencv中还有一种直方图均衡化,它是一种局部直方图均衡化,也就是是说把整个图像分成许多小块(比如按10*10作为一个小块),那么对每个小块进行均衡化。
- createCLAHE函数原型:createCLAHE([, clipLimit[, tileGridSize]])
- clipLimit参数表示对比度的大小。
- tileGridSize参数表示每次处理块的大小 。
- clahe = cv.createCLAHE(clipLimit=5, tileGridSize=(8, 8))
- dst = clahe.apply(gray) 把clahe这种局部直方图均衡化应用到灰度图gray
import cv2 as cv
import numpy as np
# 全局
# 直方图的应用
#直方图的均衡化 (即调整图像的对比度) 直方图即统计个像素点的频率
def equalzeHist_demo(img):
gray = cv.cvtColor(img, cv.COLOR_BGR2GRAY)
dst = cv.equalizeHist(gray)
cv.imshow("equalzeHist_demo", dst)
# 局部直方图均衡化
def clahe_demo(img):
gray = cv.cvtColor(img, cv.COLOR_BGR2GRAY)
cla = cv.createCLAHE(clipLimit=5.0, tileGridSize=(8, 8))
dst = cla.apply(gray)
cv.imshow("clahe_demo", dst)
def create_rgb_hist(img):
h, w, c = img.shape
rgbHist = np.zeros([16*16*16, 1], np.float32)
binsize = 256/16
for row in range(h):
for col in range(w):
b = img[row, col, 0]
g = img[row, col, 1]
r = img[row, col, 2]
index = np.int(b/binsize)*16*16 + np.int(g/binsize)*16 + np.int(r/binsize)
rgbHist[np.int(index), 0] = rgbHist[np.int(index), 0] + 1
return rgbHist
def hist_compare(image1, image2):
hise1 = create_rgb_hist(image1)
hise2 = create_rgb_hist(image2)
match1 = cv.compareHist(hise1, hise2, cv.HISTCMP_BHATTACHARYYA)
match2 = cv.compareHist(hise1, hise2, cv.HISTCMP_CORREL)
match3 = cv.compareHist(hise1, hise2,cv.HISTCMP_CHISQR)
print("巴氏距离:{}\n 相关距离:{}\n 卡方距离:{}\n".format(match1, match2, match3))
src = cv.imread("D:/JetBrains/PyCharm 2019.1.3/PycharmProjects/Study/img/Penguins.jpg") # 读取图像
src1 = cv.imread("D:/JetBrains/PyCharm 2019.1.3/PycharmProjects/Study/img/s.jpg")
cv.namedWindow("input img", cv.WINDOW_AUTOSIZE)
cv.imshow("img1", src) # 显示图像
cv.imshow("img2", src1)
equalzeHist_demo(src)
clahe_demo(src)
hist_compare(image1=src, image2=src1)
cv.waitKey(0)
cv.destroyAllWindows()
这两幅图差异很明显 如若导入相同图片,则