hgtrans.m

神经网络VC++代码 人工神经网络原理及仿真实例.

function r=hgtrans(bw,r_num,theta_num)
%Hough变换
%bw 为二值化后的图像，r_num 为Hough 直线截距的个数，
%theta_num 为Hough空间角度的个数
[m,n]=size(bw);
r_max=sqrt((m-1)^2+(n/2)^2);r_min=-n/2+1;
dr=(r_max-r_min)/(r_num-1); dtheta=pi/(theta_num-1);
%dr截距分辨率，dtheta角度分辨率

theta=8*dtheta:dtheta:(r_num-1)*dtheta;
theta=[theta,0:dtheta:7*dtheta];
s=sin(theta); c=cos(theta);
[i,j]=find(bw);
x=j-n/2;y=m-i;  %x,y为二值图像在图像空间实际的坐标
h=(x*c+y*s-r_min)/dr;
h=round(h);

r=zeros(r_num,theta_num);
for i=0:r_num-1
   y=(h==i);
   r(i+1,:)=sum(y);
end