九色国产,午夜在线视频,新黄色网址,九九色综合,天天做夜夜做久久做狠狠,天天躁夜夜躁狠狠躁2021a,久久不卡一区二区三区

打開APP
userphoto
未登錄

開通VIP,暢享免費(fèi)電子書等14項(xiàng)超值服

開通VIP
Python Opencv進(jìn)行識(shí)別相似圖片
在網(wǎng)上看到python做圖像識(shí)別的相關(guān)文章后,真心感覺python的功能實(shí)在太強(qiáng)大,因此將這些文章總結(jié)一下,建立一下自己的知識(shí)體系。

當(dāng)然了,圖像識(shí)別這個(gè)話題作為計(jì)算機(jī)科學(xué)的一個(gè)分支,不可能就在本文簡(jiǎn)單幾句就說清,所以本文只作基本算法的科普向。

看到一篇博客是介紹這個(gè),但他用的是PIL中的Image實(shí)現(xiàn)的,感覺比較麻煩,于是利用Opencv庫(kù)進(jìn)行了更簡(jiǎn)潔化的實(shí)現(xiàn)。

相關(guān)背景

要識(shí)別兩張相似圖像,我們從感性上來談是怎么樣的一個(gè)過程?首先我們會(huì)區(qū)分這兩張相片的類型,例如是風(fēng)景照,還是人物照。風(fēng)景照中,是沙漠還是海洋,人物照中,兩個(gè)人是不是都是國(guó)字臉,還是瓜子臉(還是倒瓜子臉……哈哈……)。

那么從機(jī)器的角度來說也是這樣的,先識(shí)別圖像的特征,然后再相比。

很顯然,在沒有經(jīng)過訓(xùn)練的計(jì)算機(jī)(即建立模型),那么計(jì)算機(jī)很難區(qū)分什么是海洋,什么是沙漠。但是計(jì)算機(jī)很容易識(shí)別到圖像的像素值。

因此,在圖像識(shí)別中,顏色特征是最為常用的。(其余常用的特征還有紋理特征、形狀特征空間關(guān)系特征等)

其中又分為

  • 直方圖

  • 顏色集

  • 顏色矩

  • 聚合向量

  • 相關(guān)圖

直方圖計(jì)算法

這里先用直方圖進(jìn)行簡(jiǎn)單講述。

先借用一下戀花蝶的圖片,



從肉眼來看,這兩張圖片大概也有八成是相似的了。

在python中利用opencv中的calcHist()方法獲取其直方圖數(shù)據(jù),返回的結(jié)果是一個(gè)列表,使用matplotlib,畫出了這兩張圖的直方圖數(shù)據(jù)圖

如下:



是的,我們可以明顯的發(fā)現(xiàn),兩張圖片的直方圖還是比較重合的。所以利用直方圖判斷兩張圖片的是否相似的方法就是,計(jì)算其直方圖的重合程度即可。

計(jì)算方法如下:


其中g(shù)i和si是分別指兩條曲線的第i個(gè)點(diǎn)。

最后計(jì)算得出的結(jié)果就是就是其相似程度。

不過,這種方法有一個(gè)明顯的弱點(diǎn),就是他是按照顏色的全局分布來看的,無法描述顏色的局部分布和色彩所處的位置。

也就是假如一張圖片以藍(lán)色為主,內(nèi)容是一片藍(lán)天,而另外一張圖片也是藍(lán)色為主,但是內(nèi)容卻是妹子穿了藍(lán)色裙子,那么這個(gè)算法也很可能認(rèn)為這兩張圖片的相似的。

緩解這個(gè)弱點(diǎn)有一個(gè)方法就是利用Imagecrop方法把圖片等分,然后再分別計(jì)算其相似度,最后綜合考慮。

圖像指紋與漢明距離

在介紹下面其他判別相似度的方法前,先補(bǔ)充一些概念。第一個(gè)就是圖像指紋

圖像指紋和人的指紋一樣,是身份的象征,而圖像指紋簡(jiǎn)單點(diǎn)來講,就是將圖像按照一定的哈希算法,經(jīng)過運(yùn)算后得出的一組二進(jìn)制數(shù)字。

說到這里,就可以順帶引出漢明距離的概念了。

假如一組二進(jìn)制數(shù)據(jù)為101,另外一組為111,那么顯然把第一組的第二位數(shù)據(jù)0改成1就可以變成第二組數(shù)據(jù)111,所以兩組數(shù)據(jù)的漢明距離就為1

簡(jiǎn)單點(diǎn)說,漢明距離就是一組二進(jìn)制數(shù)據(jù)變成另一組數(shù)據(jù)所需的步驟數(shù),顯然,這個(gè)數(shù)值可以衡量?jī)蓮垐D片的差異,漢明距離越小,則代表相似度越高。漢明距離為0,即代表兩張圖片完全一樣。

如何計(jì)算得到漢明距離,請(qǐng)看下面三種哈希算法

平均哈希法(aHash)

此算法是基于比較灰度圖每個(gè)像素與平均值來實(shí)現(xiàn)的

一般步驟

  • 1.縮放圖片,一般大小為8*8,64個(gè)像素值。

  • 2.轉(zhuǎn)化為灰度圖

  • 3.計(jì)算平均值:計(jì)算進(jìn)行灰度處理后圖片的所有像素點(diǎn)的平均值,直接用numpy中的mean()計(jì)算即可。

  • 4.比較像素灰度值:遍歷灰度圖片每一個(gè)像素,如果大于平均值記錄為1,否則為0.

  • 5.得到信息指紋:組合64個(gè)bit位,順序隨意保持一致性。

最后比對(duì)兩張圖片的指紋,獲得漢明距離即可。

感知哈希算法(pHash)

平均哈希算法過于嚴(yán)格,不夠精確,更適合搜索縮略圖,為了獲得更精確的結(jié)果可以選擇感知哈希算法,它采用的是DCT(離散余弦變換)來降低頻率的方法

一般步驟:

  • 縮小圖片:32 * 32是一個(gè)較好的大小,這樣方便DCT計(jì)算

  • 轉(zhuǎn)化為灰度圖

  • 計(jì)算DCT:利用Opencv中提供的dct()方法,注意輸入的圖像必須是32位浮點(diǎn)型,所以先利用numpy中的float32進(jìn)行轉(zhuǎn)換

  • 縮小DCT:DCT計(jì)算后的矩陣是32 * 32,保留左上角的8 * 8,這些代表的圖片的最低頻率

  • 計(jì)算平均值:計(jì)算縮小DCT后的所有像素點(diǎn)的平均值。

  • 進(jìn)一步減小DCT:大于平均值記錄為1,反之記錄為0.

  • 得到信息指紋:組合64個(gè)信息位,順序隨意保持一致性。

最后比對(duì)兩張圖片的指紋,獲得漢明距離即可。

dHash算法

相比pHash,dHash的速度要快的多,相比aHash,dHash在效率幾乎相同的情況下的效果要更好,它是基于漸變實(shí)現(xiàn)的。

步驟:

  • 縮小圖片:收縮到9*8的大小,以便它有72的像素點(diǎn)

  • 轉(zhuǎn)化為灰度圖

  • 計(jì)算差異值:dHash算法工作在相鄰像素之間,這樣每行9個(gè)像素之間產(chǎn)生了8個(gè)不同的差異,一共8行,則產(chǎn)生了64個(gè)差異值

  • 獲得指紋:如果左邊的像素比右邊的更亮,則記錄為1,否則為0.
    最后比對(duì)兩張圖片的指紋,獲得漢明距離即可。

整個(gè)的代碼實(shí)現(xiàn)如下:

  1. # -*- coding: utf-8 -*-
  2. #feimengjuan
  3. # 利用python實(shí)現(xiàn)多種方法來實(shí)現(xiàn)圖像識(shí)別
  4. import cv2
  5. import numpy as np
  6. from matplotlib import pyplot as plt
  7. # 最簡(jiǎn)單的以灰度直方圖作為相似比較的實(shí)現(xiàn)
  8. def classify_gray_hist(image1,image2,size = (256,256)):
  9. # 先計(jì)算直方圖
  10. # 幾個(gè)參數(shù)必須用方括號(hào)括起來
  11. # 這里直接用灰度圖計(jì)算直方圖,所以是使用第一個(gè)通道,
  12. # 也可以進(jìn)行通道分離后,得到多個(gè)通道的直方圖
  13. # bins 取為16
  14. image1 = cv2.resize(image1,size)
  15. image2 = cv2.resize(image2,size)
  16. hist1 = cv2.calcHist([image1],[0],None,[256],[0.0,255.0])
  17. hist2 = cv2.calcHist([image2],[0],None,[256],[0.0,255.0])
  18. # 可以比較下直方圖
  19. plt.plot(range(256),hist1,'r')
  20. plt.plot(range(256),hist2,'b')
  21. plt.show()
  22. # 計(jì)算直方圖的重合度
  23. degree = 0
  24. for i in range(len(hist1)):
  25. if hist1[i] != hist2[i]:
  26. degree = degree (1 - abs(hist1[i]-hist2[i])/max(hist1[i],hist2[i]))
  27. else:
  28. degree = degree 1
  29. degree = degree/len(hist1)
  30. return degree
  31. # 計(jì)算單通道的直方圖的相似值
  32. def calculate(image1,image2):
  33. hist1 = cv2.calcHist([image1],[0],None,[256],[0.0,255.0])
  34. hist2 = cv2.calcHist([image2],[0],None,[256],[0.0,255.0])
  35. # 計(jì)算直方圖的重合度
  36. degree = 0
  37. for i in range(len(hist1)):
  38. if hist1[i] != hist2[i]:
  39. degree = degree (1 - abs(hist1[i]-hist2[i])/max(hist1[i],hist2[i]))
  40. else:
  41. degree = degree 1
  42. degree = degree/len(hist1)
  43. return degree
  44. # 通過得到每個(gè)通道的直方圖來計(jì)算相似度
  45. def classify_hist_with_split(image1,image2,size = (256,256)):
  46. # 將圖像resize后,分離為三個(gè)通道,再計(jì)算每個(gè)通道的相似值
  47. image1 = cv2.resize(image1,size)
  48. image2 = cv2.resize(image2,size)
  49. sub_image1 = cv2.split(image1)
  50. sub_image2 = cv2.split(image2)
  51. sub_data = 0
  52. for im1,im2 in zip(sub_image1,sub_image2):
  53. sub_data = calculate(im1,im2)
  54. sub_data = sub_data/3
  55. return sub_data
  56. # 平均哈希算法計(jì)算
  57. def classify_aHash(image1,image2):
  58. image1 = cv2.resize(image1,(8,8))
  59. image2 = cv2.resize(image2,(8,8))
  60. gray1 = cv2.cvtColor(image1,cv2.COLOR_BGR2GRAY)
  61. gray2 = cv2.cvtColor(image2,cv2.COLOR_BGR2GRAY)
  62. hash1 = getHash(gray1)
  63. hash2 = getHash(gray2)
  64. return Hamming_distance(hash1,hash2)
  65. def classify_pHash(image1,image2):
  66. image1 = cv2.resize(image1,(32,32))
  67. image2 = cv2.resize(image2,(32,32))
  68. gray1 = cv2.cvtColor(image1,cv2.COLOR_BGR2GRAY)
  69. gray2 = cv2.cvtColor(image2,cv2.COLOR_BGR2GRAY)
  70. # 將灰度圖轉(zhuǎn)為浮點(diǎn)型,再進(jìn)行dct變換
  71. dct1 = cv2.dct(np.float32(gray1))
  72. dct2 = cv2.dct(np.float32(gray2))
  73. # 取左上角的8*8,這些代表圖片的最低頻率
  74. # 這個(gè)操作等價(jià)于c 中利用opencv實(shí)現(xiàn)的掩碼操作
  75. # 在python中進(jìn)行掩碼操作,可以直接這樣取出圖像矩陣的某一部分
  76. dct1_roi = dct1[0:8,0:8]
  77. dct2_roi = dct2[0:8,0:8]
  78. hash1 = getHash(dct1_roi)
  79. hash2 = getHash(dct2_roi)
  80. return Hamming_distance(hash1,hash2)
  81. # 輸入灰度圖,返回hash
  82. def getHash(image):
  83. avreage = np.mean(image)
  84. hash = []
  85. for i in range(image.shape[0]):
  86. for j in range(image.shape[1]):
  87. if image[i,j] > avreage:
  88. hash.append(1)
  89. else:
  90. hash.append(0)
  91. return hash
  92. # 計(jì)算漢明距離
  93. def Hamming_distance(hash1,hash2):
  94. num = 0
  95. for index in range(len(hash1)):
  96. if hash1[index] != hash2[index]:
  97. num = 1
  98. return num
  99. if __name__ == '__main__':
  100. img1 = cv2.imread('10.jpg')
  101. cv2.imshow('img1',img1)
  102. img2 = cv2.imread('11.jpg')
  103. cv2.imshow('img2',img2)
  104. degree = classify_gray_hist(img1,img2)
  105. #degree = classify_hist_with_split(img1,img2)
  106. #degree = classify_aHash(img1,img2)
  107. #degree = classify_pHash(img1,img2)
  108. print degree
  109. cv2.waitKey(0)



本站僅提供存儲(chǔ)服務(wù),所有內(nèi)容均由用戶發(fā)布,如發(fā)現(xiàn)有害或侵權(quán)內(nèi)容,請(qǐng)點(diǎn)擊舉報(bào)。
打開APP,閱讀全文并永久保存 查看更多類似文章
猜你喜歡
類似文章
老外講解什么是直方圖以及如何用opencv繪制直方圖(譯文)
改變灰度圖像直方圖的均值和標(biāo)準(zhǔn)差
特征提取代碼總結(jié)
opencv讀取RGB直方圖
OpenCV_顏色直方圖的計(jì)算、顯示、處理、對(duì)比及反向投影
python+opencv圖像處理(十四)
更多類似文章 >>
生活服務(wù)
熱點(diǎn)新聞
分享 收藏 導(dǎo)長(zhǎng)圖 關(guān)注 下載文章
綁定賬號(hào)成功
后續(xù)可登錄賬號(hào)暢享VIP特權(quán)!
如果VIP功能使用有故障,
可點(diǎn)擊這里聯(lián)系客服!

聯(lián)系客服