基于opencv的卡尺找圆

一个没有素质的人

已于 2023-11-07 16:24:25 修改

阅读量4.7k

点赞数 1

分类专栏： Opencv 文章标签： opencv

于 2022-07-02 16:07:55 首次发布

本文链接：https://blog.csdn.net/weixin_43493903/article/details/125573775

版权

Opencv 专栏收录该内容

5 篇文章 22 订阅

订阅专栏

与卡尺找线的原理一样，只是把直线模型换成了圆模型。直线的在这里基于opencv的卡尺找线
边缘点的筛选还是用RANSNC筛一遍，

void CCaliperGraphics::RansacCircleFiler(const vector<Point2d>& points, vector<Point2d>& vpdExceptPoints, double sigma)
{
	unsigned int n = points.size();

	if (n < 3)
	{
		return;
	}

	RNG random;
	double bestScore = -1.;
	vector<Point2d>vpdTemp;
	int iterations = log(1 - 0.99) / (log(1 - (1.00 / n)));

	for (int k = 0; k < iterations; k++)
	{
		int i1 = 0, i2 = 0, i3 = 0;
		Point2d p1(0, 0), p2(0, 0), p3(0, 0);
		while (true)
		{
			i1 = random(n);
			i2 = random(n);
			i3 = random(n);
			if ((i1 != i2 && i1 != i3 && i2 != i3))
			{
				if ((points[i1].y != points[i2].y) && (points[i1].y != points[i3].y))
				{
					break;
				}
			}
		}
		p1 = points[i1];
		p2 = points[i2];
		p3 = points[i3];

		//use three points to caculate a circle
		Point2d pdP12 = GetPPCenter(p1, p2);
		double dK1 = -1 / GetLineSlope(p1, p2);
		double dB1 = pdP12.y - dK1 * pdP12.x;
		Point2d pdP13 = GetPPCenter(p1, p3);
		double dK2 = -1 / GetLineSlope(p1, p3);
		double dB2 = pdP13.y - dK2 * pdP13.x;
		Point2d pdCenter(0, 0);
		pdCenter.x = (dB2 - dB1) / (dK1 - dK2);
		pdCenter.y = dK1 * pdCenter.x + dB1;
		double dR = GetPPDistance(pdCenter, p1);
		double score = 0;
		vpdTemp.clear();
		for (int i = 0; i < n; i++)
		{
			double d = dR - GetPPDistance(points[i], pdCenter);
			if (fabs(d) < sigma)
			{
				score += 1;
			}
			else
			{
				vpdTemp.push_back(points[i]);
			}
		}
		if (score > bestScore)
		{
			bestScore = score;
			vpdExceptPoints = vpdTemp;
		}
	}
}