Emotion Explorer

• points - 2Dポイントセットを入力します
  
• retval - 返却値 ellipse / rotatedRectデータ
  返却値データは下記のタプル形式をとります。
  
  

  ((cx, cy), (h, w), deg)

  
  
  • cx - 中心X
    
  • cy - 中心Y
    
  • h - 楕円縦方向の長さ
    
  • w - 楕円横方向の長さ
    
  • deg - 傾き角度
    
  
  この関数を実行する前処理として輪郭を求める必要があります。
  ちなみに外接矩形の角4点を求める場合、cv2.boxPoints()を使います。
  

contours,hierarchy =  cv2.findContours(binimg,cv2.RETR_TREE,cv2.CHAIN_APPROX_SIMPLE)
for i, cnt in enumerate(contours):
    ellipse = cv2.fitEllipse(cnt)
    resimg = cv2.ellipse(resimg,ellipse,(255,0,0),2)

    #外接矩形4点を求める
    box = cv2.boxPoints(ellipse)
    box = np.int0(box)
    im = cv2.drawContours(im,[box],0,(0,0,255),2)

C/C++の場合のfitEllipse()返却値について

 

今回のテスト用画像はこの画像を使います。それぞれの図を楕円フィッテングしてみて
どうなるか？試します。

• 一般的なグレイスケール・二値化
  三角形が二値化できなかった例をしめします。
  
  一般的なグレイスケール化サンプル
  

  gray = cv2.cvtColor(img, cv2.COLOR_RGB2GRAY)
  _, binimg = cv2.threshold(gray, 0, 255, cv2.THRESH_BINARY | cv2.THRESH_OTSU)
  

  
  
  二値化結果
   
  
  手動でしきい値を設定するのも有なんですけれど、それもなんかやだな。
  



• カスタムなグレイスケール・二値化
  RGB３プレーンをすべてAND合成してグレイスケールを作ります。
  
  カスタムなグレイスケール化
  

  gray1 = cv2.bitwise_and(img[:,:,0], img[:,:,1])
  gray1 = cv2.bitwise_and(gray1, img[:,:,2])
  _, binimg = cv2.threshold(gray1 , 0, 255, cv2.THRESH_BINARY | cv2.THRESH_OTSU)
  

  
  
  二値化結果
   
  
  出来ていますね。
  

import cv2
import math
import numpy as np

def main():
    img = cv2.imread('ellipse2.png', cv2.IMREAD_COLOR)

    # グレイスケール化
    gray1 = cv2.bitwise_and(img[:,:,0], img[:,:,1])
    gray1 = cv2.bitwise_and(gray1, img[:,:,2])

    # 二値化
    _, binimg = cv2.threshold(gray1, 0, 255, cv2.THRESH_BINARY | cv2.THRESH_OTSU)
    binimg = cv2.bitwise_not(binimg)

    # 結果画像の黒の部分を灰色にする。
    bimg = binimg // 4 + 255 * 3 //4
    resimg = cv2.merge((bimg,bimg,bimg)) 

    # 輪郭取得
    contours,hierarchy =  cv2.findContours(binimg,cv2.RETR_TREE,cv2.CHAIN_APPROX_SIMPLE)
    for i, cnt in enumerate(contours):
        # 楕円フィッティング
        ellipse = cv2.fitEllipse(cnt)
        print(ellipse)

        cx = int(ellipse[0][0])
        cy = int(ellipse[0][1])

        # 楕円描画
        resimg = cv2.ellipse(resimg,ellipse,(255,0,0),2)
        cv2.drawMarker(resimg, (cx,cy), (0,0,255), markerType=cv2.MARKER_CROSS, markerSize=10, thickness=1)
        cv2.putText(resimg, str(i+1), (cx+3,cy+3), cv2.FONT_HERSHEY_SIMPLEX, 0.8, (0,80,255), 1,cv2.LINE_AA)

    cv2.imshow('resimg',resimg)
    cv2.waitKey()

if __name__ == '__main__':
    main()

• 2020/3/13
  OpenCV4でfindContours()の返却値が変わりましたので修正しました。
  
  

  contours, hierarchy = cv.findContours(...)

  
  ソースコードも修正しておきます。
  また、動作環境も更新しました。
  


• 2020/10/27
  cv2.fitEllipse()の返却値とC++の場合について、少し加筆しました。


• 2021/3/14
  cv2.fitEllipse()の関数の説明を加筆。


• 2021/3/23
  cv2.fitEllipse()の関数の説明をさらに加筆。


• 2021/05/02 修正
  IntersectionObserver's による画像のオフスクリーン遅延読み込み処理に変更。
  IE11未対応。