我的程序中有这个代码:
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应该大得多),即使轮廓是椭圆本身。我怎么解决它?
答案 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)) );
希望这会有所帮助:)