霍夫变换检测圆和直线

希望变换可以检测所有图形一、实验目的1.掌握MATLAB软件的使用和设计过程。掌握Hopf变换的实现。3.使用MATLAB语言设计基于Hopf变换的随机图形识别。二、实验设备或设备安装MATLAB软件计算机三、总体设计原则和流程图1、编程原理霍夫变换的基本思想是将图像平面上的点映射到参数平面上的曲线上，最后通过统计特性解决问题。具有抗噪音性能好、对部分盖子不敏感等特性。2、流程图程序启动打开图形基于霍夫曼变换的简单图形识别希望变换传感线希望转换检测源标识获得的图像保存图像程序结束四、主要方案代码1、Hough变换检测圆程序Clear allClose allClc% cr，ST=circle _ product；%st(1)-x，st(2)-y，st(3)-rCr=im read(sample . BMP)；%circle1.bmpFigureIMS how(Cr)；row，range=size(Cr)；P=row * rangea=0；Y=zeros(1，p)；Q=zeros(1，p)；For k1=1:1:rowFor k2=1:1:rangeIf cr(k1，k2)=0a=a 1；XXX=k2-1；YYY=row 1-k1；y(a)=YYY；q(a)=XXX；EndEndEndFor k3=1:1:pIf Y(k3)=0Q(k3)=0BreakEndEnd% z=ones(1，5)；% a=1:1:300B=1:136300% 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)*)；% 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% 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:136300B=1:136300Rrr=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).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).2 (b-y (K5)。 2)；r0=(ABS(Z1-z2)=. 1 ABS(Z1=z3)=. 1 ABS(Z1=Z4)=. 1 ABS(Z1=z5)=. 1 ABS()Rr=r0。* z1Aa=r0。* a；Bb=r0。* b；m，n=size(RR)；If rr=zeros(1，n)；as=as 1；ContinueEndpz1，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:pz2RR(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% endEnd% if aa=zeros(1，300)% as=1% endRangex=1:1:rangeRowy=1:1:rowrx，ry=meshgrid (rangex，rowiy)；RT=RXFigure网面(rx、ry、rt)；霍尔德温；Plot3(XX、YY、RR、zh)；% 5rhao，nn=the_max1(XX，YY，RR)；I1=ones(row，range)；RhaoCenter=rhao% t2=1.5T3=1.5For k1=1:1:rowFor k2=1:1:rangeXXX=k2-1；YYY=row 1-k1；if ABS(XXX-center(1)*(XXX-center(1)(yyy-center(2)*(yyy-center(1)I1(k1，k2)=0；EndEndEndImwrite(i1，找到的bmp，BMP)；I2=imread(找到的圆。BMP)；FigureIMS how(I2)；%st2、霍夫变换传感线Clear allClose allClcI=im read(sample . BMP)；row，range=size(I)；P=row * rangea=0；Y=zeros(1，p)；Q=zeros(1，p)；For k1=1:1:rowFor k2=1:1:rangeIf i(k1，k2)=0a=a 1；XXX=k2-1；YYY=row 1-k1；y(a)=YYY；q(a)=YYY-XXX；EndEndEnd% for t=1:1:a% o=Y(t) Q(t)% endx1=0 1；y1=Y(1)Q(1)；Plot(X1，Y1，r)；霍尔德温；YY=zeros(1，a-1)；QQ=zeros(1，a-1)；For t=2:1:ax2=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)% qqqqq=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