霍夫变换试验.doc

实验报告册数字图形图像处理霍夫变换试验学 院： 人民武装学院学院 专 业： 计算机科学与技术班 级： 11级计科班 学 号： 1120070544 学生姓名： 苏靖 指导教师： 实验目的：1掌握MATLAB软件的使用，以及其设计流程； 2掌握霍夫变换的实现方法；3用MATLAB语言设计基于霍夫变换的任意图形的识别。4.掌握图像分割的基本原理及应用。实验内容：1. 通过全局阈值及局部阈值实现灰度图像二值化；2. 通过Hough检测图像中的直线。二、实验仪器或设备 装MATLAB软件的微机一台 三、总体设计原理及流程图 1、程序设计的原理 霍夫变换的基本思想就是把图像平面上的点对应到参数平面上的曲线，最后通过统计特性来解决问题。具有良好的抗噪声性能和对部分遮盖的不敏感等特性 。2、程序流程图四、实验要求1. 程序结构清晰，运行结果正确。2. 描述上述程序的设计、实现和结果，并对结果进行分析。五、主要程序代码 1、霍夫变换检测圆程序 clear all; close all; clc; % cr,st=circle_product;%st(1)-x,st(2)-y,st(3)-r cr=imread(sample.bmp);%circle1.bmp figure; 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 end end for k3=1:1:p if Y(k3)=0&Q(k3)=0 break; end end % 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=z3)=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:k7k6=0;% z1=sqrt(a-Q(k5).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&abs(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; end pz1,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; end k11=k8+1; % if rr=zeros(1,n); % as=as+1; continue; % end end % if aa=zeros(1,300) % as=1 % end rangex=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);%5 rhao,nn=the_max1(XX,YY,RR); i1=ones(row,range);rhao center=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-center(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:row 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 end end % for t=1:1:a % o=Y(t) Q(t) % end X1=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=QQ(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 % end rhao=the_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); rhao t2=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: