hough.m

来自「是一个用MATLAB编的一个系统」· M 代码 · 共 27 行

%tic 
function res=hough(im,RHO_MAX,THETA_MAX) 
[X,Y]=size(im); 
d_rho=X/RHO_MAX; 
d_theta=pi/THETA_MAX; 

theta=-pi/4:d_theta:pi-pi/4-d_theta; 

smat=sin(theta); 
cmat=cos(theta); 

[x,y]=find(im); 
% translation by a pixel so that low left pixel has (0,0) 
%coordinates 
x=x-1; 
y=y-1; 

h= (y * smat + x * cmat )/d_rho; 
h=round(h); 

res=zeros(RHO_MAX,THETA_MAX); 
for j=0:RHO_MAX-1 
    temp=(h==j); 
    res(j+1,:)=sum(temp); 
end 
%toc