pythonで取得した全ての座標が入る最小の円を描画し、その中心点を求める方法を教えてください。

Question

上記の画像の青色の三角形の座標を取得するソフトウェアを作りました。その導線を追うために、取得した座標が全て入る最小の円を描画し、その中心点を求める機能を追加したいですが、やり方がわかりません。

画像の青色の三角形の座標を取得するプログラムは

# -*- coding:utf-8 -*-
import cv2
import numpy as np

# 青色の検出
def detect_blue_color(img):
    # HSV色空間に変換
    hsv = cv2.cvtColor(img, cv2.COLOR_BGR2HSV)

    # 青色のHSVの値域1
    hsv_min = np.array([40,140,120])
    hsv_max = np.array([150,255,255])

    # 青色領域のマスク（255：赤色、0：赤色以外）    
    mask = cv2.inRange(hsv, hsv_min, hsv_max)

    # マスキング処理
    masked_img = cv2.bitwise_and(img, img, mask=mask)

    return mask, masked_img

for i in range(392):
    img = cv2.imread('CoD/filename_{0:05d}.png'.format(i))

    # 色検出（赤、緑、青）
    mask, masked_img = detect_blue_color(img)

    # 輪郭を抽出する。
    contours, hierarchy = cv2.findContours(
        mask, cv2.RETR_EXTERNAL, cv2.CHAIN_APPROX_SIMPLE
    )

    # 小さい輪郭は誤検出として削除する
    contours = list(filter(lambda x: cv2.contourArea(x) > 10, contours))


    def get_center(contour):
        """輪郭の中心を取得する。
        """
        # 輪郭のモーメントを計算する。
        M = cv2.moments(contour)
        # モーメントから重心を計算する。
        cx = M["m10"] / M["m00"]
        cy = M["m01"] / M["m00"]

        return cx, 127-cy

    # 輪郭の中心を取り出す。
    centers = [get_center(x) for x in contours]

    # 結果を CSV で保存する。
    with open('centers2.csv','a') as f_handle: 
     np.savetxt(f_handle, centers, delimiter=",", fmt="%.0f")

で、cxと127-cyが座標を表しています。その座標をcsvファイルにて保存しているのがこのプログラムです。
詳しい方よろしくお願いします。

Answer 1

精度はintですがOpenCVのminEnclosingCircleで求めることができます。

    # 略

    def get_center(contour):

        # 略

        #return cx, 127-cy
        return cx, cy # 修正

    # 輪郭の中心を取り出す。
    centers = [get_center(x) for x in contours]

    # 確認用に各中心を描画
    centers = np.array(centers, dtype=int)
    for c in centers:
        img = cv2.circle( img, tuple(c), 3, (0,0,255), 1)

    # 最小外接円を求める
    cs = centers.reshape(len(centers), 1, 2).astype(int) # カタチを整える
    (x,y),radius = cv2.minEnclosingCircle(cs)
    print(x,y,radius) # 中心座標と半径

    # 描画
    img = cv2.circle( img, (int(x), int(y)), int(radius), (255,0,0), 1)
    cv2.imwrite( 'ret.png', img)

Comments

[gigagiga64]

ありがとうございました！参考にさせていただきます！