求关于椭圆检测的程序

junren

请教
谁有检测椭圆的程序。
具体说来就是：在一幅图像中，有一幅已经边缘检测好的图像。首先要验证图像中存在椭圆，（这个椭圆不一定是标准椭圆！）其次，确定椭圆的中心，长半轴，短半轴，（倾斜角度）。

[ 本帖最后由 eight 于 2008-5-14 17:06 编辑 ]

frogfish

hough变换检测椭圆的MATLAB程序

1. [row col]=size(fedge);
   
2. minofa=a;
   
3. maxofa=round(row/2);
   
4. minofy0=round(col/2)-30;
   
5. maxofy0=round(col/2)+30;
   
6. minofb=round(col/2)-60;
   
7. maxofb=round(col/2);
   
8. maxofx=round(row/2);
   
9. scalor=4;
   
10. H=zeros(floor((maxofa-minofa)/scalor)+1,floor((maxofa-minofa)/scalor)+1,...
    
11. floor((maxofy0-minofy0)/scalor)+1,floor((maxofb-minofb)/scalor)+1);
    
12. for x=1:maxofx
    
13.    for y=1:col
    
14.      temp=fedge(x,y);
    
15.        if temp==255
    
16.          for a=minofa:scalor:maxofa
    
17.            for x0=a:scalor:maxofa
    
18.              for b=minofb:scalor:maxofb
    
19.                for y0=minofy0:scalor:maxofy0
    
20.                   temp=((y-y0)/b)^2+((x-x0)/a)^2;
    
21.                     if abs(temp-1)<=0.01
    
22.                          xtemp=floor((x0-minofa)/scalor)+1;
    
23.                          atemp=floor((a-minofa)/scalor)+1;
    
24.                          ytemp=floor((y0-minofy0)/scalor)+1;
    
25.                          btemp=floor((b-minofb)/scalor)+1;
    
26.                          H(xtemp,atemp,ytemp,btemp)=H(xtemp,atemp,ytemp,btemp)+1;
    
27.                      end
    
28.                 end
    
29.               end
    
30.             end
    
31.           end
    
32.         end
    
33.    end
    
34. end
    
35. 
    
36. maxofH=max(max(max(max(H))));
    
37. for i=1:floor((maxofa-minofa)/scalor)+1
    
38.     for j=1:floor((maxofa-minofa)/scalor)+1
    
39.         for m=1:floor((maxofy0-minofy0)/scalor)+1
    
40.            for n=1:floor((maxofb-minofb)/scalor)+1
    
41.                  temp=H(i,j,m,n);
    
42.                    if temp==maxofH
    
43.                        xtemp=i;
    
44.                        atemp=j;
    
45.                        ytemp=m;
    
46.                        btemp=n;
    
47.                        break;
    
48.                     end
    
49.             end
    
50.          end
    
51.      end
    
52. end
    
53. 
    
54. x0=(xtemp-1)*scalor+minofa;
    
55. a=(atemp-1)*scalor+minofa;
    
56. y0=(ytemp-1)*scalor+minofy0;
    
57. b=(btemp-1)*scalor+minofb;
    
58. figure;
    
59. imshow(fedge,[]);
    
60. hold on;
    
61. plot(y0,x0,'r+');
    
62. t=0:pi/180:2*pi;
    
63. plot(y0+b*sin(t),x0+a*cos(t),'r-');
    
64. hold off;

复制代码

看看这个好不好使

ivylife

请问下a是什么？
fedge是什么？正好要用到这个程序。多谢了