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

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

ただいまの
回答率

88.63%

ハフ変換によって検出した2つの円の中心間距離

解決済

回答 1

投稿

  • 評価
  • クリップ 0
  • VIEW 2,464

Taka11

score 14

 前提・実現したいこと

ハフ変換で円検出を行っています。
最終的に実現したいことは、検出した2つの円の中心間の距離の算出です。
そのために円の座標情報を得たいです。

print (circles)


で返される値は円の中心座標と円の直径?を示しているという認識でよろしいでしょうか?

また、円の中心間の距離を算出する方法もあわせて教えていただけますと幸いです。

 該当のソースコード

import cv2
import numpy as np

img = cv2.imread('C:/Users/ito/Anaconda3/envs/Sample/pic/testimg/IMG_0634.jpg', 0)
img = cv2.medianBlur(img, 5)
cimg = cv2.cvtColor(img, cv2.COLOR_GRAY2BGR)

circles = cv2.HoughCircles(img, cv2.HOUGH_GRADIENT, 1, 30, param1=50, param2=30, minRadius=0, maxRadius=0)

circles = np.uint16(np.around(circles))
for i in circles[0,:]:
    cv2.circle(cimg, (i[0], i[1]), i[2], (0,255,0), 2)
    cv2.circle(cimg, (i[0], i[1]), 2, (0,0,255), 3)


print(circles)

cv2.imshow('detected circles', cimg)
cv2.waitKey(0)
cv2.destroyAllWindows
cv2.imwrite('C:/Users/ito/Anaconda3/envs/Sample/pic/result/Hough_sample01.jpg', cimg)

 補足情報(FW/ツールのバージョンなど)

windows7 professional
python3.6.5
opencv3.3.1

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

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

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

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 過去に投稿した質問と同じ内容の質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 1

checkベストアンサー

+2

N 個の円が検出された場合、2つの円の中心間の距離ということは、combination(N, 2) 通りの組み合わせが考えられます。
このように「複数の点同士の距離を調べる」という場合、「距離行列」を考えましょう。
距離行列については 過去の回答 を参考にしてください。

cv2.HoughCircles() は (1, NumCircles, 3) の numpy 配列を返し、各要素は (円の中心 x 座標、円の中心 y 座標、円の半径) を表すことに留意すると、以下のようになります。

 サンプルコード

 円を検出する。

イメージ説明
使用した画像

import cv2
import matplotlib.pyplot as plt
import numpy as np

img = cv2.imread('test.png')
img = cv2.medianBlur(img, 5)
gray = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY)

circles = cv2.HoughCircles(
    gray, cv2.HOUGH_GRADIENT, 1, 30, param1=50, param2=30, minRadius=0, maxRadius=0)
print(circles.shape)  # (1, 4, 3)
circles = circles[0]  # (1, NumCircles, 3) -> (NumCircles, 3)

# 検出された円を描画
for x, y, radius in circles:
    cv2.circle(img, (x, y), radius, (0, 255, 0), 2)  # 弧
    cv2.circle(img, (x, y), 2, (0, 0, 255), -1)  # 中心
plt.imshow(img)
plt.show()

イメージ説明

 距離行列を作成する。

# 円の中心同士の距離行列作成
centers = circles[:, :2]
dist_matrix = np.linalg.norm(centers[:, np.newaxis] - centers, axis=-1)
print(dist_matrix)
# [[  0.       127.40487  107.200745 243.3105  ]
#  [127.40487    0.       111.8034   116.51609 ]
#  [107.200745 111.8034     0.       196.86543 ]
#  [243.3105   116.51609  196.86543    0.      ]]

for i, j in zip(*np.triu_indices(len(centers), k=1)):
    print('円 {} の中心と円 {} の中心の距離: {:.2f}'.format(i, j, dist_matrix[i, j]))
 0 の中心と円 1 の中心の距離: 127.40
円 0 の中心と円 2 の中心の距離: 107.20
円 0 の中心と円 3 の中心の距離: 243.31
円 1 の中心と円 2 の中心の距離: 111.80
円 1 の中心と円 3 の中心の距離: 116.52
円 2 の中心と円 3 の中心の距離: 196.87

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2018/10/17 19:15

    ご回答いただきありがとうございます!
    試してみたところ、以下のようなエラーが出てしまうのですが、、


    IndexError Traceback (most recent call last)
    <ipython-input-12-023872d1b40c> in <module>()
    8
    9 for i, j in zip(*np.triu_indices(4, k=1)):
    ---> 10 print('円 {} の中心と円 {} の中心の距離: {:.2f}'.format(i, j, dist_matrix[i, j]))
    11
    12

    IndexError: index 2 is out of bounds for axis 1 with size 2

    キャンセル

  • 2018/10/17 19:16

    これはどういうことなのでしょうか、、?

    キャンセル

  • 2018/10/17 19:19

    すいません。検出された円の個数がハードコーディングされていたので、修正しました。

    ```
    for i, j in zip(*np.triu_indices(4, k=1)):
    ```

    ```
    for i, j in zip(*np.triu_indices(len(centers), k=1)):
    ```

    これでどうでしょうか?

    キャンセル

  • 2018/10/17 19:21

    丁寧にありがとうございます!
    できました!ありがとうございます!

    キャンセル

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

  • ただいまの回答率 88.63%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

同じタグがついた質問を見る