画像の輪郭を楕円でフィッティングする方法2

Question

図の中央の黒い輪の外側の輪郭を楕円でフィッティングしたいのです。
周囲の黒いつぶつぶを除去できれば、OpeningやClosingを用いて、楕円でのフィッティングができそうです。
周囲の黒いつぶつぶはどう除去するといいでしょうか。

fourteenlength様のご回答を少し変更して、次のように楕円のフィッティングができました。
しかし少し大きめに評価しており、困っています。
やはり、周囲の黒いつぶつぶが邪魔です。
消す方法にいいアイディアあれば、教えていただきたいです。

---追記---
画面に接するくらいの大きさの円のマスクを使って、周りをすべて取り除いたところ、うまくいきました。

Answer 1

import cv2
import numpy as np

# 白黒で読み込んでボカす
img = cv2.imread('img.png',0)
img_mod =  cv2.GaussianBlur(img,(11,11),0)

# 極端な値で二値化する
# ※ただし、白部分を検索するため、背景を黒にするように注意
ret,thresh = cv2.threshold(img_mod,240,255,cv2.THRESH_BINARY_INV)
cv2.imshow("thresh",thresh)

# 領域に分ける
contours,hierarchy = cv2.findContours(thresh, cv2.RETR_EXTERNAL, cv2.CHAIN_APPROX_TC89_L1)

# 最大値を探す
area = []
for contour in contours:
    area.append(cv2.contourArea(contour))
idx = area.index(max(area))

# マスクを作る
img_mask = np.zeros_like(img)
img_mask = cv2.fillPoly(img_mask, [contours[idx]],(255))
cv2.imshow("img_mask",img_mask)

# マスクに従ってノイズを消す
img_mod[img_mask==0] = 0
cv2.imshow("img_mod",img_mod)


# とりあえず線を引く
# img_mod = cv2.drawContours(img, contours, -1, (128), 3)

# 最大値に線を引く
# img_mod = cv2.polylines(img, [contours[idx]],True,(128), thickness=3)
# cv2.imshow("img_mod",img_mod)
# cv2.waitKey(0)

cv2.imwrite("img_mod.png",img_mod)
cv2.imwrite("img_mask.png",img_mask)

マスク画像

ノイズ塗りつぶし

Comments

[退会済みユーザー]

OpenCVのHough検出だと(真)円に近似した検出しかできなさそうですが、scikit-imageのHough検出であれば楕円も検出できそうですね。
https://sabopy.com/py/scikit-image-22/

[surphy]

ご回答ありがとうございます。
初心者なもので、いただいたコードを解読中です。

ひとまず、詳細を理解する前にコード内のcontours[idx]を楕円でフィッティング
ellipse = cv2.fitEllipse(cnt)
することで、やりたいことができました。

丁寧にどうもありがとうございます。

[退会済みユーザー]

> cv2.fitEllipse
こんな関数があったんですね、こちらこそ勉強になりました。
今日の午前中に気になってScikit-imageの楕円フィッティングを試したのですがフリーズして困ってました。

[surphy]

私もまずはScikit-imageの楕円フィッティングを走らせてみました。
同じくフリーズしました。
深追いしなかったですが、楕円が複数検出されてしまうためでしょうかね。

こちらも少しお役に立てたみたいで、うれしいです。