python OpenCV學(xué)習(xí)筆記之繪制直方圖的方法
本篇文章主要介紹了python OpenCV學(xué)習(xí)筆記之繪制直方圖的方法,小編覺得挺不錯的,現(xiàn)在分享給大家,也給大家做個參考。一起跟隨小編過來看看吧
官方文檔 – https://docs.opencv.org/3.4.0/d1/db7/tutorial_py_histogram_begins.html
直方圖會讓你對圖像的強(qiáng)度分布有一個全面的認(rèn)識。它是一個在x軸上帶有像素值(從0到255,但不總是),在y軸上的圖像中對應(yīng)的像素數(shù)量的圖。
這只是理解圖像的另一種方式。通過觀察圖像的直方圖,你可以直觀地看到圖像的對比度、亮度、強(qiáng)度分布等?,F(xiàn)在幾乎所有的圖像處理工具都提供了直方圖的特性。下面是劍橋彩色網(wǎng)站的圖片,建議去訪問這個網(wǎng)站,了解更多細(xì)節(jié)。

你可以看到圖像和它的直方圖。(這個直方圖是用灰度圖像繪制的,而不是彩色圖像)。直方圖的左邊部分顯示了圖像中較暗像素的數(shù)量,右邊區(qū)域顯示了更明亮的像素。從直方圖中可以看到,深色區(qū)域的像素數(shù)量比亮色區(qū)域更多,而中間色調(diào)的數(shù)量(中值大約在127左右)則少得多。
直方圖
現(xiàn)在我們已經(jīng)知道了什么是直方圖,我們可以看看如何找到它。OpenCV和Numpy都有內(nèi)置的功能。在使用這些函數(shù)之前,我們需要了解一些與直方圖相關(guān)的術(shù)語。
BINS:上面的直方圖顯示了每個像素值的像素數(shù),從0到255。您需要256個值來顯示以上的直方圖。但是,考慮一下,如果您不需要單獨查找所有像素值的像素數(shù)量,而是在一個像素值區(qū)間內(nèi)的像素數(shù)量,該怎么辦?例如,你需要找到介于0到15之間的像素數(shù),然后是16到31……240到255。您只需要16個值來表示這個直方圖。OpenCV Tutorials on histograms中展示了這個例子。
所以你要做的就是把整個直方圖分成16個子部分,每個子部分的值是所有像素數(shù)的和。每個子部分都被稱為“BIN”。在第一種情況下,BINS的數(shù)量是256(每個像素一個),而在第二種情況下,它只有16個。在OpenCV文檔中,用術(shù)語 histSize 表示 BINS。
DIMS:它是我們收集數(shù)據(jù)的參數(shù)的個數(shù)。在這種情況下,我們收集的數(shù)據(jù)只有一件事,強(qiáng)度值。所以這里是1。
RANGE:它是你想測量的強(qiáng)度值的范圍。通常,它是 [ 0,256 ],也就是所有的強(qiáng)度值。
OpenCV中直方圖的計算
現(xiàn)在我們使用cv.calcHist()函數(shù)來找到直方圖。讓我們熟悉一下這個函數(shù)及其參數(shù):
cv.calcHist(images, channels, mask, histSize, ranges[, hist[, accumulate]])
images:它是uint8類型或float32的源圖像。它應(yīng)該用方括號括起來,也就是”[img]”。
channels:它也用方括號括起來。它是我們計算直方圖的信道的索引。例如,如果輸入是灰度圖像,它的值是0。對于顏色圖像,您可以通過0、1或2來分別計算藍(lán)色、綠色或紅色通道的直方圖。
mask:遮罩圖。為了找到完整圖像的直方圖,它被指定為“None”。但如果你想找到圖像的特定區(qū)域的直方圖,你必須為它創(chuàng)建一個遮罩圖,并將其作為遮罩。
histSize:這代表了我們的BINS數(shù)。需要用方括號來表示。在整個范圍內(nèi),我們通過了256。
ranges:強(qiáng)度值范圍,通常是 [ 0,256 ]
讓我們從一個樣本圖像開始。只需在灰度模式下加載圖像并找到其完整的直方圖。
img = cv.imread('home.jpg', 0)
hist = cv.calcHist([img], [0], None, [256], [0,256])
hist是一個256x1陣列,每個值對應(yīng)于該圖像中的像素值及其對應(yīng)的像素值。
Numpy中直方圖的計算
Numpy中提供了np.histogram()方法
hist, bins = np.histogram(img.ravel(), 356, [0,256])
hist和之前計算的一樣。但是bins有257個元素,因為Numpy計算bins是以0-0.99,1-1.99等,所以最后一個是255-255.99。為了表示這一點,他們還在bins的末端添加了256。但我們不需要256。到255就足夠了。
Numpy還有另一個函數(shù),np.bincount(),比np.histograme()要快得多(大約10X)。對于一維直方圖,你可以試一下。不要忘記在np.bincount中設(shè)置minlength=256。例如,hist=np.bincount(img.ravel(),minlength=256)
OpenCV函數(shù)比np.histogram()快(大約40X)。所以堅持用OpenCV函數(shù)。
繪制直方圖
1、使用Matplotlib
Matplotlib有一個繪制直方圖的函數(shù):matplotlib.pyplot.hist()
它直接找到了直方圖并繪制了它。您不需要使用calcHist()或np.histogram()函數(shù)來找到直方圖。看下面的代碼:
import numpy as np
import cv2 as cv
from matplotlib import pyplot as plt
img = cv.imread('home.jpg', 0)
plt.hist(img.ravel(), 256, [0,256])
plt.show()

或者你可以用正常的matplotlib,這對BGR的情節(jié)很有幫助。為此,您需要首先找到直方圖數(shù)據(jù)。試試下面的代碼:
import numpy as np
import cv2 as cv
from matplotlib import pyplot as plt
img = cv.imread('home.jpg')
color = ('b', 'g', 'r')
for i, col in enumerate(color):
histr = cv.calcHist([img], [i], None, [256], [0,256])
plt.plot(histr, color=col)
plt.xlim([0,256])
plt.show()
你可以從上面的圖中扣除,藍(lán)色在圖像中有一些高值區(qū)域(很明顯,它應(yīng)該是由天空引起的)
2、使用OpenCV
這里,你可以調(diào)整直方圖的值和它的bin值,讓它看起來像x,y坐標(biāo),這樣你就可以用cv.line()或cv.polyline()函數(shù)來繪制它,從而生成與上面相同的圖像。這已經(jīng)是OpenCV-Python2官方的樣本了。檢查sampl/python/hist.py的代碼。
應(yīng)用遮罩
我們用cv.calcHist()函數(shù)來找一張完整的圖片的直方圖。但是我們只要圖片的一部分的直方圖呢?在你想要找到的區(qū)域中,創(chuàng)建一個帶有白色的遮罩圖像。然后把它作為遮罩。
img = cv.imread('home.jpg', 0)
# create a mask
mask = np.zeros(img.shape[:2], np.uint8)
mask[100:300, 100:400] = 255
masked_img = cv.bitwise_and(img, img, mask=mask)
#Calculate histogram with mask and without mask
Check third argument for mask
hist_full = cv.calcHist([img], [0], None, [256], [0,256])
hist_mask = cv.calcHist([img], [0], mask, [256], [0,256])
plt.subplot(221), plt.imshow(img, 'gray')
plt.subplot(222), plt.imshow(mask,'gray')
plt.subplot(223), plt.imshow(masked_img, 'gray')
plt.subplot(224), plt.plot(hist_full), plt.plot(hist_mask)
plt.xlim([0,256])
plt.show()
藍(lán)線表示完整圖片的直方圖
綠線表示遮罩之后的直方圖

以上就是本文的全部內(nèi)容,希望對大家的學(xué)習(xí)有所幫助,也希望大家多多支持腳本之家。
相關(guān)文章
python中三種高階函數(shù)(map,reduce,filter)詳解
在Python中,函數(shù)其實也是一種數(shù)據(jù)類型,今天重點給大家介紹python中三種高階函數(shù)(map,reduce,filter)的相關(guān)知識,感興趣的朋友一起看看吧2021-10-10
python通過Seq2Seq實現(xiàn)閑聊機(jī)器人
這篇文章主要介紹了python通過Seq2Seq實現(xiàn)閑聊機(jī)器人,文中有非常詳細(xì)的代碼示例,對正在學(xué)習(xí)python的小伙伴們有很好的幫助,需要的朋友可以參考下2021-04-04
Python統(tǒng)計單詞出現(xiàn)的次數(shù)
最近經(jīng)理交給我一項任務(wù),統(tǒng)計一個文件中每個單詞出現(xiàn)的次數(shù),列出出現(xiàn)頻率最多的5個單詞。本文給大家?guī)砹藀ython 統(tǒng)計單詞次數(shù)的思路解析,需要的朋友參考下吧2018-04-04
利用Python通過商品條形碼查詢商品信息的實現(xiàn)示例
這篇文章主要介紹了利用Python通過商品條形碼查詢商品信息,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2022-07-07

