获取多个圆的圆心

 

原始图片：

                                                    图A

                                                          图B

 

两种方法：

 

第一种方法：扫描整个图像，对于每个圆形区域，最先扫描到的点应为其顶点A，如图1所示，获取顶点A和底部点B的坐标，AB的中点即为圆心坐标。同时将该圆形区域的灰度值置为其他值，这样在扫描下一行时就不会在扫描到该圆。

                                     图1

 

实现代码：

int label(byte [,] f, int x, int y, int L)//将连通区域（圆）的灰度值标记为L { if(f[x,y]==255) { f[x,y] = (byte)L; label(f, x-1, y, L); label(f, x+1, y, L); label(f, x, y-1, L); label(f, x, y+1, L); label(f, x-1, y-1, L); label(f, x-1, y+1, L); label(f, x+1, y-1, L); label(f, x+1, y+1, L); return 1; } return 0; } void GetCen(byte [,] f)//获取圆心 { int w = f.GetLength(0); int h = f.GetLength(1); int cenx = 0;//圆心横坐标 int ceny_up = 0;//圆顶点的纵坐标 int ceny_bottom = 0;//圆底部点的纵坐标 for(int y = 0; y < h; y++) for(int x = 0; x < w; x++) { if( f[x,y] == 255 ) { label(f, x, y, 128);//将该连通区域（圆）灰度值标记为128 cenx = x;//获取圆顶点横坐标，圆顶点横坐标与圆心横坐标相同 ceny_up = ceny_bottom = y;//获取圆顶点纵坐标 while( f[cenx, ceny_bottom] == 128 )//获取圆底部点的纵坐标 ceny_bottom++; f[cenx,(ceny_up + ceny_bottom)/2] = 254;//将圆心灰度值标记为254，(cenx, (ceny_up + ceny_bottom)/2)即为圆心坐标 } } } void main() { byte [,] f = LoadImg(); ShowImg("",f); GetCen(f); ShowImg("Label", f); }

 

对A图的运算效果如图2：效果可以接受

                                图4

                                  图2

 

对B图的运算效果如图3：由于B图中的圆有些不是很光滑，顶部有些地方是直线，所以第一个扫描的点不是A，而是A前面的某一点，所以有些圆并未取到圆心。

 

                                图3

 

针对B图这种情况，采用另一种算法，求质心法。质心的求法为：

，

其中，f(x,y)是(x,y)点处的灰度值。

仍然扫描图像，对每个圆形区域求质心，并将区域的灰度值重新标定，以免该区域被再次扫描到，代码：

int label(byte [,] f, int x, int y, int L, ref int sum_x, ref int sum_y, ref int sum)//将连通区域的灰度值标记为L，并求连通区域的质心,ref表示引用用法 { if(f[x,y]==255) { f[x,y] = (byte)L; sum_x += x; sum_y += y; sum++; label(f, x-1, y, L, ref sum_x, ref sum_y, ref sum); label(f, x+1, y, L, ref sum_x, ref sum_y, ref sum); label(f, x, y-1, L, ref sum_x, ref sum_y, ref sum); label(f, x, y+1, L, ref sum_x, ref sum_y, ref sum); label(f, x-1, y-1, L, ref sum_x, ref sum_y, ref sum); label(f, x-1, y+1, L, ref sum_x, ref sum_y, ref sum); label(f, x+1, y-1, L, ref sum_x, ref sum_y, ref sum); label(f, x+1, y+1, L, ref sum_x, ref sum_y, ref sum); return 1; } return 0; } void GetCen(byte [,] f) { int w = f.GetLength(0); int h = f.GetLength(1); int sum_x = 0; int sum_y = 0; int sum = 0; for(int x = 0; x < w; x++) for(int y = 0; y < h; y++) { if( f[x,y] == 255 ) { label(f, x, y, 128, ref sum_x, ref sum_y, ref sum); f[sum_x/sum, sum_y/sum] = 254;//圆心坐标即质心坐标：(sum_x/sum, sum_y/sum) } sum_x = sum_y = sum = 0; } } void main() { byte [,] f = LoadImg(); ShowImg("",f); GetCen(f); ShowImg("Label", f); }

对B图的效果见图4：效果可以接受

                          图4