opencv如何获取指定圆上所有像素点的坐标？？

[tiger波波]

在opencv里，根据圆心与半径，可以画出一个圆。
请问大家：如何得到这个指定圆上所有像素点的坐标？？？

[赵4老师]

将圆用比如正90边形近似，调用GDI函数LineDDA？

[赵4老师]

或者简单点：
创建一副尺寸够大的黑色位图，在上面画一个白色圆，遍历每个白色像素对应的坐标。

[zgl7903]

数学计算不就可以? x=x0 + r*cosθ; y=y0 + r*sinθ;

[向立天]

你有圆心和半径的不就可以自己算了

[tiger波波]

引用 3 楼 zgl7903 的回复:

数学计算不就可以? x=x0 + r*cosθ; y=y0 + r*sinθ;

我觉得opencv画圆函数，在一张给定的图片上画出的圆，其圆上的像素点是确定的，算出来的和真实画出的圆有差距。另外角度不好取，同样一组角度，大圆和小圆效果就差很多。

[赵4老师]

算出来的值和实际像素位置难免有舍入误差。

[tiger波波]

采用了赵老师的办法弄了一个

	void get_circle_coordinates(int width, int height, cv::Point center, int radius, std::vector<cv::Point> &coordinates)
	{
		cv::Scalar color;

		coordinates.clear();

		cv::Mat image(width, height, CV_8UC3, cv::Scalar(0, 0, 0));
		circle(image, center, radius, cv::Scalar(255, 255, 255), 1);

		for (int i = 0; i < width; i++)
		{
			for (int j = 0; j < height; j++)
			{
				COpenCV::get_color(image, cv::Point(i, j), color);
				if ((255 == color.val[0]) && (255 == color.val[1]) && (255 == color.val[2]))
					coordinates.push_back(cv::Point(i, j));
			}
		}
	}

[向立天]

引用 5 楼 youyingbo 的回复:

[quote=引用 3 楼 zgl7903 的回复:]
数学计算不就可以? x=x0 + r*cosθ; y=y0 + r*sinθ;

我觉得opencv画圆函数，在一张给定的图片上画出的圆，其圆上的像素点是确定的，算出来的和真实画出的圆有差距。另外角度不好取，同样一组角度，大圆和小圆效果就差很多。[/quote]那你就遍历所有像素点，根据颜色值来重新校对一下

[唔枣]

引用 8 楼 youyingbo 的回复:

采用了赵老师的办法弄了一个

	void get_circle_coordinates(int width, int height, cv::Point center, int radius, std::vector<cv::Point> &coordinates)
	{
		cv::Scalar color;

		coordinates.clear();

		cv::Mat image(width, height, CV_8UC3, cv::Scalar(0, 0, 0));
		circle(image, center, radius, cv::Scalar(255, 255, 255), 1);

		for (int i = 0; i < width; i++)
		{
			for (int j = 0; j < height; j++)
			{
				COpenCV::get_color(image, cv::Point(i, j), color);
				if ((255 == color.val[0]) && (255 == color.val[1]) && (255 == color.val[2]))
					coordinates.push_back(cv::Point(i, j));
			}
		}
	}

引用 8 楼 youyingbo 的回复:

采用了赵老师的办法弄了一个

	void get_circle_coordinates(int width, int height, cv::Point center, int radius, std::vector<cv::Point> &coordinates)
	{
		cv::Scalar color;

		coordinates.clear();

		cv::Mat image(width, height, CV_8UC3, cv::Scalar(0, 0, 0));
		circle(image, center, radius, cv::Scalar(255, 255, 255), 1);

		for (int i = 0; i < width; i++)
		{
			for (int j = 0; j < height; j++)
			{
				COpenCV::get_color(image, cv::Point(i, j), color);
				if ((255 == color.val[0]) && (255 == color.val[1]) && (255 == color.val[2]))
					coordinates.push_back(cv::Point(i, j));
			}
		}
	}

画圆函数的最后一个参数要改成“-1”，不然画的是圆不是全白的圆，而是边缘像素为白，中间为黑

[唔枣]

引用 8 楼 youyingbo 的回复:

采用了赵老师的办法弄了一个

	void get_circle_coordinates(int width, int height, cv::Point center, int radius, std::vector<cv::Point> &coordinates)
	{
		cv::Scalar color;

		coordinates.clear();

		cv::Mat image(width, height, CV_8UC3, cv::Scalar(0, 0, 0));
		circle(image, center, radius, cv::Scalar(255, 255, 255), 1);

		for (int i = 0; i < width; i++)
		{
			for (int j = 0; j < height; j++)
			{
				COpenCV::get_color(image, cv::Point(i, j), color);
				if ((255 == color.val[0]) && (255 == color.val[1]) && (255 == color.val[2]))
					coordinates.push_back(cv::Point(i, j));
			}
		}
	}

谢谢大佬的代码，帮我解决了一个问题

[sunmeal]

判断一下， 优化一：一条直线与圆最多只能交与两点，做好这个计算，可以避免过多的遍历。 优化二：直接用sin计算，然后+-1个像素寻找是否存在，一般圆主要是折线时 优化三：注意四个象限的中间折线的那个点的变化规则，可以直接找到的 优化四：直接读OpenCV源码，将画圆的运算方式，找出即可

[GK_2014]

引用 8 楼 tiger波波 的回复:

采用了赵老师的办法弄了一个

	void get_circle_coordinates(int width, int height, cv::Point center, int radius, std::vector<cv::Point> &coordinates)
	{
		cv::Scalar color;

		coordinates.clear();

		cv::Mat image(width, height, CV_8UC3, cv::Scalar(0, 0, 0));
		circle(image, center, radius, cv::Scalar(255, 255, 255), 1);

		for (int i = 0; i < width; i++)
		{
			for (int j = 0; j < height; j++)
			{
				COpenCV::get_color(image, cv::Point(i, j), color);
				if ((255 == color.val[0]) && (255 == color.val[1]) && (255 == color.val[2]))
					coordinates.push_back(cv::Point(i, j));
			}
		}
	}

引用 8 楼 tiger波波 的回复:

采用了赵老师的办法弄了一个

	void get_circle_coordinates(int width, int height, cv::Point center, int radius, std::vector<cv::Point> &coordinates)
	{
		cv::Scalar color;

		coordinates.clear();

		cv::Mat image(width, height, CV_8UC3, cv::Scalar(0, 0, 0));
		circle(image, center, radius, cv::Scalar(255, 255, 255), 1);

		for (int i = 0; i < width; i++)
		{
			for (int j = 0; j < height; j++)
			{
				COpenCV::get_color(image, cv::Point(i, j), color);
				if ((255 == color.val[0]) && (255 == color.val[1]) && (255 == color.val[2]))
					coordinates.push_back(cv::Point(i, j));
			}
		}
	}

既然给定了圆的半径，那么要创建的图片的宽和高用2倍半径就行了，没必要从函数外再传入宽和高。另外创建的图片单通道就行了，为什么还要3通道呢？这种办法对于效率要求不高的场景还是可以用用，更高效的方法就是看opencv在调用circle函数时怎么通过中心点和半径获取圆上点的集合。

[GK_2014]

根据opencv源码改写的函数，根据圆心坐标和半径，输出圆上所有点集。

vector<Point> getCirclePts(Point center, int radius)
{
    vector<Point> ret;
    int err = 0, dx = radius, dy = 0, plus = 1, minus = (radius << 1) - 1;
    while(dx >= dy)
    {
        int mask;
        int y11 = center.y - dy, y12 = center.y + dy, y21 = center.y - dx, y22 = center.y + dx;
        int x11 = center.x - dx, x12 = center.x + dx, x21 = center.x - dy, x22 = center.x + dy;

        ret.push_back(Point(x11,y11));
        ret.push_back(Point(x11,y12));
        ret.push_back(Point(x12,y11));
        ret.push_back(Point(x12,y12));
        ret.push_back(Point(x21,y21));
        ret.push_back(Point(x21,y22));
        ret.push_back(Point(x22,y21));
        ret.push_back(Point(x22,y22));

        dy++;
        err += plus;
        plus += 2;
        mask = (err <= 0) - 1;
        err -= minus & mask;
        dx += mask;
        minus -= mask & 2;
    }
    return ret;
}