霍夫變換檢測圓和直線

1、霍夫變換檢測任意形狀一、 實驗目的1 掌握MATLAB軟件的使用，以及其設計流程；2 掌握霍夫變換的實現(xiàn)方法；3 用MATLAB語言設計基于霍夫變換的任意圖形的識別。二、 實驗儀器或設備 裝MATLAB軟件的微機一臺三、 總體設計原理及流程圖1、程序設計的原理霍夫變換的基本思想就是把圖像平面上的點對應到參數(shù)平面上的曲線，最后通過統(tǒng)計特性來解決問題。具有良好的抗噪聲性能和對部分遮蓋的不敏感等特性 。2、程序流程圖程序開始打開圖形基于霍夫曼變換的簡單圖形識別霍夫變換檢測直線霍夫變換檢測圓識別得到的圖像保存圖像程序結束四、 主要程序代碼1、 霍夫變換檢測圓程序clear all; close all

2、; clc;% cr,st=circle_product;%st(1)-x,st(2)-y,st(3)-rcr=imread(sample.bmp);%circle1.bmpfigure;imshow(cr);row,range=size(cr);p=row*range;a=0;Y=zeros(1,p);Q=zeros(1,p);for k1=1:1:row for k2=1:1:range if cr(k1,k2)=0 a=a+1; XXX=k2-1;YYY=row+1-k1; Y(a)=YYY;Q(a)=XXX; end endendfor k3=1:1:p if Y(k3)=0&Q(k3

3、)=0 break; endend% z=ones(1,5);% a=1:1:300;b=1:1:300;% z1=sqrt(a-Q(1).2+(b-Y(1).2);% z2=sqrt(a-Q(2).2+(b-Y(2).2);% z3=sqrt(a-Q(3).2+(b-Y(3).2);% z4=sqrt(a-Q(4).2+(b-Y(4).2);% z5=sqrt(a-Q(5).2+(b-Y(5).2);% % z1=sqrt(a-Q(1)*(a-Q(1)+(b-Y(1)*(b-Y(1);% % Z1=(a-Q(1).2+(b-Y(1).2;r0=(abs(z1-z2)=0.1&abs(z1=z

4、3)=0.1&abs(z1=z4)=0.1&abs(z1=z5)=0.1&abs(z2=z3)=0.1&abs(z2=z4)=0.1&abs(z2=z5)=0.1);% r=r0.*z1;% aa=r0.*a;% bb=r0.*b;% rrr=round(r(r0=0)% aaa=aa(r0=0)% bbb=bb(r0=0)k7=floor(k3/6);k5=2;a=1:1:300;b=1:1:300;rrr=zeros(1,p);aaa=zeros(1,p);bb=zeros(1,p);k6=0;as=0;k11=1;for k5=1:1:k7% k6=0;%z1=sqrt(a-Q(k5).

5、2+(b-Y(k5).2);z2=sqrt(a-Q(k5+1).2+(b-Y(k5+1).2);z3=sqrt(a-Q(k5+2).2+(b-Y(k5+2).2);z4=sqrt(a-Q(k5+3).2+(b-Y(k5+3).2);z5=sqrt(a-Q(k5+4).2+(b-Y(k5+4).2);z6=sqrt(a-Q(k5+5).2+(b-Y(k5+5).2);r0=(abs(z1-z2)=.1&abs(z1=z3)=.1&abs(z1=z4)=.1&abs(z1=z5)=.1&abs(z1=z6)=.1&abs(z2=z3)=1&abs(z2=z4)=.1&abs(z2=z5)=.1&a

6、bs(z2=z6)=.1);rr=r0.*z1;aa=r0.*a;bb=r0.*b;m,n=size(rr);if rr=zeros(1,n); as=as+1; continue;endpz1,pz=size(rr(r0=0);rrr=round(rr(r0=0)aaa=aa(r0=0)bbb=bb(r0=0)pz2=pz+k11-1;k9=1;for k8=k11:pz2 RR(k8)=rrr(k9); XX(k8)=aaa(k9); YY(k8)=bbb(k9); k9=k9+1;endk11=k8+1;% if rr=zeros(1,n);% as=as+1; continue;% e

7、ndend% if aa=zeros(1,300)% as=1% endrangex=1:1:range;rowy=1:1:row;RX,RY=meshgrid(rangex,rowy);RT=RX;figure;mesh(RX,RY,RT);hold on;plot3(XX,YY,RR,ko);%5rhao,nn=the_max1(XX,YY,RR); i1=ones(row,range);rhaocenter=rhao;% t2=1.5;t3=1.5; for k1=1:1:row for k2=1:1:range XXX=k2-1;YYY=row+1-k1; if abs(XXX-cen

8、ter(1)*(XXX-center(1)+(YYY-center(2)*(YYY-center(2)-center(3)*center(3)=101 i1(k1,k2)=0; end end end imwrite(i1,找出的圓.bmp,bmp); i2=imread(找出的圓.bmp); figure;imshow(i2);%st2、 霍夫變換檢測直線clear all; close all; clc;i=imread(sample.bmp);row,range=size(i);p=row*range;a=0;Y=zeros(1,p);Q=zeros(1,p);for k1=1:1:ro

9、w for k2=1:1:range if i(k1,k2)=0 a=a+1; XXX=k2-1;YYY=row+1-k1; Y(a)=YYY;Q(a)=YYY-XXX; end endend% for t=1:1:a% o=Y(t) Q(t)% endX1=0 1;Y1=Y(1) Q(1);plot(X1,Y1,r);hold on;YY=zeros(1,a-1);QQ=zeros(1,a-1);for t=2:1:a X2=0 1;Y2=Y(t) Q(t); plot(X2,Y2); YY(t-1),QQ(t-1)=pll(X1,Y1,X2,Y2);% yyy=YY(t-1)% qqq=Q

10、Q(t-1)%看交點 plot(YY(t-1),QQ(t-1),ko);end % res=gather_point(YY,QQ,QQ);nn=0;res% m,n3=size(YY);% for no1=1:1:n3% n=0;% for no2=(no1+1):1:n3% if (abs(res(no2,1)-res(no1,1)=0.0001&abs(res(no2,2)-res(no1,2)=0.0001&abs(res(no2,2)-res(no1,2)=n3/5% rhao=hao;break;% elseif nnn% rhao=hao;nn=n;% end% endrhao=t

11、he_max1(YY,QQ,QQ); % yyy=YY% qqq=QQ% t2=0;t3=0;% for t1=2:1:t% t2=YY(t1)+t2;% t3=t3+QQ(t1);% end% t2=t2/(t1-1)% t3=t3/(t1-1)%t2=m,t3=b% i1=ones(row,range);% t2=2;t3=0; i1=ones(row,range);rhaot2=rhao(1);t3=rhao(2);t2=round(t2);t3=round(t3);% t2=1.5;t3=1.5; for k1=1:1:row for k2=1:1:range XXX=k2-1;YYY=row+1-k1; if abs(YYY-(t2*XXX+t3)=1%YYY=t2*XXX+t3 i1(k1,k2)=0; end end en