广义霍夫变换的MATLAB源码

function [score  y x ] = Generalized_hough_transform（IsItm） 
%Find template/shape Itm in greyscale image Is using generalize hough trasform

%Use generalized hough transform to find Template/shape binary image given in binary image Itm in grayscale image Is （greyscale image）

%Is is greyscale  picture were the template Itm should be found 
%Itm is bool edge image of the template with edges markedd ones
%Return the xy location  cordniates  which gave the best match 
%Also return the score of each this point （number of point matching）

%The xy are the cordinates in image Is in which the  the top left edge of image Itm （11） should be positioned in order to give the best match

%Is=imread（‘‘）;
%Itm=imread（‘‘）;
%if nargin<3 thresh=1;end;

%--------------------------create edge and system edge images------------------------------------------------------------------------------------------------------------------------

%Is=rgb2gray（Is）;

Iedg=edge（Is‘canny‘）; % Take canny edge images of Is with automatic threshold
%}
%--------------------------------------------------------------------------------------------------------------------------------------
[y x]=find（Itm>0）; % find all yx cordinates of all points equal 1 inbinary template image Itm
nvs=size（x）;% number of points in the  template image
%-------------------Define Yc and Xc ----------------------------------------------
Cy=1;%round（mean（y））;% find object y center note that any reference point will do so the origin of axis hence 1 could be used just as well
Cx=1;%round（mean（x））;% find object z center note that any reference point will do so the origin of axis hence 1 could be used just as well

%------------------------------create gradient map of Itm distrobotion between zero to pi %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
GradientMap = gradient_direction（ Itm ）;

%%%%%%%%%%%%%%%%%%%%%%%Create an R-Table of Itm gradients to  parameter space in parameter space.%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%---------------------------create template descriptor array------------------------------------
MaxAngelsBins=30;% devide the angel space to MaxAngelsBins uniformed space bins
MaxPointsPerangel=nvs（1）;% maximal amount of points corresponding to specific angel

PointCounter=zeros（MaxAngelsBins）;% counter for the amount of edge points associate with each angel gradient
Rtable=zeros（MaxAngelsBinsMaxPointsPerangel2）; % assume maximum of 100 points per angle with MaxAngelsBins angles bins between zero and pi and xy for the vector to the center of each point
% the third adimension are vector from the point to the center of the vessel

%------------------fill the  angel bins with points in the Rtable---------------------------------------------------------
for f=1:1:nvs（1）
    bin=round（（GradientMap（y（f） x（f））/pi）*（MaxAngelsBins-1））+1; % transform from continues gradient angles to discrete angle bins and 
    PointCounter（bin）=PointCounter