質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.48%
Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Q&A

解決済

1回答

3666閲覧

python opencvで特徴量マッチング

soepu

総合スコア41

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

0グッド

0クリップ

投稿2019/06/09 08:07

python opencv で2枚の画像(比較される画像、パターン画像)を比較して、
特徴量マッチングをおこない、一致した場所を四角形でくくりたいと考えています。
いろいろ検索してインターネット上から、上記内容ができそうな
参考プログラムを見つけたため、試してみました。
作ったコード

python

1import numpy as np 2import cv2 3import sys 4 5img1 = cv2.imread("src/star.png",0) 6img2 = cv2.imread("src/starPtn.png",0) 7 8akaze = cv2.AKAZE_create() 9 10kp1, des1 = akaze.detectAndCompute(img1,None) 11kp2, des2 = akaze.detectAndCompute(img2,None) 12 13bf = cv2.BFMatcher() 14matches = bf.knnMatch(des1,des2,k = 2) 15 16good = [] 17for m,n in matches: 18 if m.distance < 0.7*n.distance: 19 good.append(m) 20 21if len(good)>MIN_MATCH_COUNT: 22 src_pts = np.float32([ kp1[m.queryIdx].pt for m in good ]).reshape(-1,1,2) 23 dst_pts = np.float32([ kp2[m.trainIdx].pt for m in good ]).reshape(-1,1,2) 24 25 M, mask = cv2.findHomography(src_pts, dst_pts, cv2.RANSAC,5.0) 26 matchesMask = mask.ravel().tolist() 27 28 h,w = img1.shape 29 pts = np.float32([ [0,0],[0,h-1],[w-1,h-1],[w-1,0] ]).reshape(-1,1,2) 30 dst = cv2.perspectiveTransform(pts,M) 31 32 print('true') 33 print(np.int32(dst)[0][0][0],np.int32(dst)[0][0][1]) 34 print(np.int32(dst)[1][0][0],np.int32(dst)[1][0][1]) 35 print(np.int32(dst)[2][0][0],np.int32(dst)[2][0][1]) 36 print(np.int32(dst)[3][0][0],np.int32(dst)[3][0][1]) 37 38else: 39 print('false') 40

print 出力結果:
true
-204 -82
-204 297
488 296
488 -81

となります。
得られた結果は、なにを表しているのでしょうか?
(座標を表しているとは思いますが、マイナス[-]がでるのは?)

使っている画像は、
上:比較画像、下:パターン画像
star.png
starPtn.png

以上、宜しくお願い致します。

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

提示コードの動作よく分かっていませんが、比較される画像とパターン画像を逆にするとそれらしい結果が得られました。
参考:Python3+OpenCV3でパッケージ認識してみましょう

Python

1import numpy as np 2import cv2 3import sys 4 5img1 = cv2.imread("starPtn.png",1) 6img2 = cv2.imread("star.png",1) 7 8akaze = cv2.AKAZE_create() 9kp1, des1 = akaze.detectAndCompute(img1,None) 10kp2, des2 = akaze.detectAndCompute(img2,None) 11 12bf = cv2.BFMatcher() 13matches = bf.knnMatch(des1,des2,k = 2) 14 15good = [] 16for m,n in matches: 17 if m.distance < 0.7*n.distance: 18 good.append(m) 19 20if len(good): 21 src_pts = np.float32([ kp1[m.queryIdx].pt for m in good ]).reshape(-1,1,2) 22 dst_pts = np.float32([ kp2[m.trainIdx].pt for m in good ]).reshape(-1,1,2) 23 24 M, mask = cv2.findHomography(src_pts, dst_pts, cv2.RANSAC,5.0) 25 matchesMask = mask.ravel().tolist() 26 27 h,w,c = img1.shape 28 pts = np.float32([ [0,0],[0,h-1],[w-1,h-1],[w-1,0] ]).reshape(-1,1,2) 29 dst = cv2.perspectiveTransform(pts,M) 30 31 dst_pts = [ (np.int32(dst)[0][0][0],np.int32(dst)[0][0][1]), 32 (np.int32(dst)[1][0][0],np.int32(dst)[1][0][1]), 33 (np.int32(dst)[2][0][0],np.int32(dst)[2][0][1]), 34 (np.int32(dst)[3][0][0],np.int32(dst)[3][0][1])] 35 36 for idx in range(len(dst_pts)): 37 cv2.line(img2, dst_pts[idx], dst_pts[(idx+1) % len(dst_pts)], (0,0,255),2) 38 39 cv2.imshow("img2", img2) 40 k = cv2.waitKey(0)

イメージ説明

投稿2019/06/09 12:03

can110

総合スコア38262

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.48%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問