【分享】基于图像处理的圆拟合

addoil

根据已经拍摄的圆照片,来对照片中的圆进行拟合,求取其圆心坐标和半径.为了简化问题,自己用phototshop画了一张圆照片.下面是我的处理方法,希望对需要这方面的东西的人有点帮助,更希望大家谈谈自己的意见,提出更好的方法.

addoil

拟合:
  
a=imread('c1.bmp'); %读入图像数据
b=im2bw(a); %图像二值化
c=edge(b); %边缘提取
  
n=1;
for i=1:70 %提取数据点坐标,70和100分别为图像的水平和竖直分辨率
     for j=1:100
         if c(i,j)==1  
             x(n)=i;
             y(n)=j;
             n=n+1;
         end;
     end;
end;
x=x';
y=y';
  
plot(x,y,'bo'), hold on;  
[xc,yc,r,f]=tlscirc(x,y) %对数据进行圆拟合   
rectangle('Position', [xc-r,yc-r,2*r,2*r],'Curvature', [1,1],'EdgeColor','b')

addoil

tlscirc是圆拟合的算法(感谢Genial兄提供的算法).
请把下面附件中的c1.bmp和tlscirc.m文件拷贝到matlab工作目录下,然后执行circle.m中的所有命令.

addoil

计算结果:
  
xc = 27.1470
yc = 51.1470
r =   13.4848(和我用逆向工程软件imageware拟合的结果很接近)

addoil

欢迎大家来讨论,提出自己的意见,提出更好的方法!:I

cwit

请问：如果你的圆是椭圆，或者一段不规则的曲线，可以很好的拟合吗？

addoil

算法的代码我还没有仔细看(这段时间太忙了,很久不见了哈,cwit版主),不敢肯定地回答你这个问题.不过我想应该是没有问题的

Genial

椭圆的拟合我见过,
这个我贴出来的这个好像不是的
一段不规则的曲线这样直接拟合好像更加不行吧
  
对于不规则曲线的拟合,我粗略地看了下svm的工具箱
它是那个做到的,在选择适当的核的时候
不过,它好像又无法正确地拟合诸如圆,椭圆之类的标准图形
  
不过,记得有人说过,其实,万能方法只能是无能方法
或许多少有些道理

addoil

”万能方法只能是无能方法“
说的好啊,Genial，在曲线拟合上要想找一个通用的方法很难啊。

njwy1978

楼上，你好，你的算法看过，但是我找不到tlscirc.m文件，能否帮我发一份。
e-mail:njwy1978@sohu.com

hf_lwb

没看到tlscirc.m文件啊，能再传一下吗？

paradiseboy

呵呵，我作过椭圆拟合的程序，比较有意思。

xjtueecjun

你好可否分享下圆拟合的算法，[email=sunboy@stu.xjtu.edu.cn谢谢sunboy@stu.xjtu.edu.cn谢谢[b[/email]] 3# addoil

xjtueecjun

你好你有椭圆拟合的程序吗？ 13# paradiseboy

messenger

见http://forum.simwe.com/thread-835047-1-1.html，http://forum.simwe.com/thread-948512-1-1.html。

你好你有椭圆拟合的程序吗？ 13# paradiseboy
xjtueecjun 发表于 2011-1-19 16:03

qzh001

dddddddddddd