霍夫變換檢測(cè)圓和直線

1、霍夫變換檢測(cè)任意形狀一、 實(shí)驗(yàn)?zāi)康? 掌握MATLAB軟件的使用，以及其設(shè)計(jì)流程；2 掌握霍夫變換的實(shí)現(xiàn)方法；3 用MATLAB語言設(shè)計(jì)基于霍夫變換的任意圖形的識(shí)別。二、 實(shí)驗(yàn)儀器或設(shè)備 裝MATLAB軟件的微機(jī)一臺(tái)三、 總體設(shè)計(jì)原理及流程圖1、程序設(shè)計(jì)的原理霍夫變換的基本思想就是把圖像平面上的點(diǎn)對(duì)應(yīng)到參數(shù)平面上的曲線，最后通過統(tǒng)計(jì)特性來解決問題。具有良好的抗噪聲性能和對(duì)部分遮蓋的不敏感等特性 。2、程序流程圖程序開始打開圖形基于霍夫曼變換的簡單圖形識(shí)別霍夫變換檢測(cè)直線霍夫變換檢測(cè)圓識(shí)別得到的圖像保存圖像程序結(jié)束四、 主要程序代碼1、 霍夫變換檢測(cè)圓程序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、 霍夫變換檢測(cè)直線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)%看交點(diǎn) 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