在OpenCV中通过轮廓矩查找椭圆参数

时间:2014-04-18 15:36:46

标签: c++ opencv image-processing ellipse

我的程序中有这个代码:

findContours(canny_output, contours, hierarchy, CV_RETR_TREE, CV_CHAIN_APPROX_SIMPLE, Point(0, 0) );
for(int i = 0; i < contours.size(); i++){
    if(hierarchy[i][3] >=0 && arcLength(contours[i], true) > 240){
        Moments mu =  moments(contours[i], false);
        Point2d mc(mu.m10/mu.m00, mu.m01/mu.m00);
        double m11 = mu.m11 / mu.m00;
        double m20 = mu.m20 / mu.m00;
        double m02 = mu.m02 / mu.m00;
        double orientation = 0.5*atan(2*m11 / (m20 - m02));
        double lambda1 = (m20 + m02)/2 + sqrt(4*m11*m11 - (m20 - m02)*(m20 - m02));
        double lambda2 = (m20 + m02)/2 - sqrt(4*m11*m11 - (m20 - m02)*(m20 - m02));
        cout << "labmda1: " << lambda1 << endl << "labmda2: " << lambda2 << endl;
        cout << "orientation: " << orientation << endl << "mc: " << mc << endl << endl;
        //ellipse(drawing_ellipses, mc, Size(lambda2, lambda1), orientation, 0, 360, Scalar(0, 0, 255), 1);
    }
}

如果我理解正确的话,lambda1和lambda2应该是椭圆的半长轴和半短轴,mc--他的中心和方向 - 旋转角度。相反,只有mc正确计算,其他参数是错误的(lambda应该大得多),即使轮廓是椭圆本身。我怎么解决它?

1 个答案:

答案 0 :(得分:0)

我可以问这是为了什么吗?我正在做类似的研究,试图制作一个基准标记。

我的研究问题是:

lambda1 =(m20 + m02)/ 2 + sqrt(4 * m11 * m11 - (m20 - m02)*(m20 - m02));

实际应该是:

lambda1 = srqt((m20 + m02)/ 2 + sqrt(4 * m11 * m11 - (m20 - m02)*(m20 - m02)) );

希望这会有所帮助:)