Hough变换重点.docVIP

%Hough变换?clc;clear;close?all?f=imread(D:\dj.jpg);? %若是彩色图片转灰度化?if?length(size(f))2?????f=rgb2gray(f);?end?%? figure(1)? subplot(121);imshow(f);? %利用edge函数进行边缘检测?j=edge(f,Sobel);? subplot(122);imshow(j);?[row,col]=size(j);?pinfang=round((row*row+col*col)^0.5);?A=zeros(2*pinfang,180);?? for?m=1:row?????for?n=1:col?????????if?j(m,n)0? ????????????for?thera=1:180? ????????????????r=thera/180*pi;??%角度转弧度? ????????????????rho=round(m*cos(r)+n*sin(r));%ρ=cosθ+sinθ?????????????????rho=rho+pinfang+1;%-l:l转换到1：2l+1?????????????????A(rho,thera)=A(rho,thera)+1;?????????????end?????????end?????end?end?? [r