python opencv实现图像配准与比较

╰+哭是因爲堅強的太久メ 2022-10-29 07:51 258阅读 0赞

本文实例为大家分享了python opencv实现图像配准与比较的具体代码，供大家参考，具体内容如下

代码

1  
2  
3  
4  
5  
6  
7  
8  
9  
10  
11  
12  
13  
14  
15  
16  
17  
18  
19  
20  
21  
22  
23  
24  
25  
26  
27  
28  
29  
30  
31  
32  
33  
34  
35  
36  
37  
38  
39  
40  
41  
42  
43  
44  
45  
46  
47  
48  
49  
50  
51  
52  
53  
54  
55  
56  
57  
58  
59  
60  
61  
62  
63  
64  
65  
66  
67  
68  
69  
70  
71  
72  
73  
74  
75  
from skimage import io  
import cv2 as cv  
import numpy as np  
import matplotlib.pyplot as plt  
  
img\_path1 = '2\_HE\_maxarea.png'  
img\_path2 = '2\_IHC\_maxarea.png'  
  
img1 = io.imread(img\_path1)  
img2 = io.imread(img\_path2)  
img1 = np.uint8(img1)  
img2 = np.uint8(img2)  
  
\# find the keypoints and descriptors with ORB  
orb = cv.ORB\_create()  
kp1, des1 = orb.detectAndCompute(img1,None)  
kp2, des2 = orb.detectAndCompute(img2,None)  
  
\# def get\_good\_match(des1,des2):  
\# bf = cv.BFMatcher()  
\# matches = bf.knnMatch(des1, des2, k=2)  
\# good = \[\]  
\# for m, n in matches:  
\# if m.distance < 0.75 \* n.distance:  
\# good.append(m)  
\# return good,matches  
\# goodMatch,matches = get\_good\_match(des1,des2)  
\# img3 = cv.drawMatchesKnn(img1,kp1,img2,kp2,matches\[:20\],None,flags=2)  
  
\# create BFMatcher object  
bf = cv.BFMatcher(cv.NORM\_HAMMING, crossCheck=True)  
\# Match descriptors.  
matches = bf.match(des1,des2)  
\# Sort them in the order of their distance.  
matches = sorted(matches, key = lambda x:x.distance)  
\# Draw first 20 matches.  
img3 = cv.drawMatches(img1,kp1,img2,kp2,matches\[:20\],None, flags=2)  
  
  
goodMatch = matches\[:20\]  
if len(goodMatch) > 4:  
ptsA= np.float32(\[kp1\[m.queryIdx\].pt for m in goodMatch\]).reshape(-1, 1, 2)  
ptsB = np.float32(\[kp2\[m.trainIdx\].pt for m in goodMatch\]).reshape(-1, 1, 2)  
ransacReprojThreshold = 4  
H, status =cv.findHomography(ptsA,ptsB,cv.RANSAC,ransacReprojThreshold);  
\#其中H为求得的单应性矩阵矩阵  
\#status则返回一个列表来表征匹配成功的特征点。  
\#ptsA,ptsB为关键点  
\#cv2.RANSAC, ransacReprojThreshold这两个参数与RANSAC有关  
imgOut = cv.warpPerspective(img2, H, (img1.shape\[1\],img1.shape\[0\]),flags=cv.INTER\_LINEAR + cv.WARP\_INVERSE\_MAP)  
  
\# 叠加配准变换图与基准图  
rate = 0.5  
overlapping = cv.addWeighted(img1, rate, imgOut, 1-rate, 0)  
io.imsave('HE\_2\_IHC.png', overlapping)  
err = cv.absdiff(img1,imgOut)  
  
\# 显示对比  
plt.subplot(221)  
plt.title('orb')  
plt.imshow(img3)  
  
plt.subplot(222)  
plt.title('imgOut')  
plt.imshow(imgOut)  
  
plt.subplot(223)  
plt.title('overlapping')  
plt.imshow(overlapping)  
  
plt.subplot(224)  
plt.title('diff')  
plt.imshow(err)  
  
plt.show()  
结果：

以上就是本文的全部内容，希望对大家的学习有所帮助