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

Question

我的程序中有以下代码：

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);
    }
}

如果我理解正确的话，λ1和λ2应该是椭圆的半长轴和半短轴，mc-他的中心和方向-旋转角。相反，只有mc计算正确，其他参数是错误的（lambda比应该的大得多），即使轮廓本身是椭圆。我该怎么修？

Answer 1

我能问一下这是干什么的吗？我正在做类似的研究，试图制作一个基准标记。

我研究的问题是：

lambda1=（m20+m02）/2+sqrt（4*m11*m11-（m20-m02）*（m20-m2））；

实际上应该是：

lambda1=srqt（（m20+m02）/2+sqrt（4*m11*m11-（m20-m02）*（m20-m2））；

希望这能有所帮助：）