python opencv 圆形roi

拟合圆：

#! /usr/bin/env python
# -*- coding: utf-8 -*-
import cv2
from numpy import *
from scipy import optimize
import functools

def countcalls(fn):
    @functools.wraps(fn)
    def wrapped(*args):
        wrapped.ncalls +=1
        return fn(*args)

    wrapped.ncalls = 0
    return wrapped

def calc_R(xc, yc):
    return sqrt((xx - xc) ** 2 + (yy - yc) ** 2)
@countcalls
def f_2(c):
    Ri = calc_R(*c)
    return Ri - Ri.mean()

def nihe(xx,yy):
    # 质心坐标
    x_m = mean(xx)
    y_m = mean(yy)
    #圆心估计
    center_estimate = x_m, y_m
    center_2, _ = optimize.leastsq(f_2, center_estimate)

    xc_2, yc_2 = center_2
    Ri_2  = calc_R(xc_2, yc_2)
    #拟合圆的半径
    R_2 = Ri_2.mea