GHT 自己编写的广义霍夫变换

自己编写的广义霍夫变换，里面有测试图片，可以作为学习霍夫变换的入门

function acc=ght（inputimageRtable）
[rowscolumns]=size（inputimage）;
BW=edge（inputimage）;
[FxFy]=gradient（inputimage）;
[rowsthcolumnst]=size（Rtable）;
d=pi/rowst;
acc=zeros（rowscolumns）;
for x=1:columns
    for y=1:rows
        if （Fx（yx）~=0）
            ang（yx）=atan（Fy（yx）/Fx（yx））;
        else
            ang（yx）=pi/2;
        end
    end
    
end
for x=1:columns
    for y=1:rows
        if （BW（yx）==1）
            phi=ang（yx）;
            i=round（（phi+（pi/2））/d）;
            if （i==0）
                i=1;
            end
               for j=1:columnst
                   if（Rtable（i1j）==0&Rtable（i2j）==0）
                       j=columnst;
                   else
                        b0=y-Rtable（i2j）;
                        a0=x-Rtable（i1

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2011-03-27 22:36  广义霍夫变换\
     文件        1204  2011-03-27 18:54  广义霍夫变换\1.jpg
     文件       10728  2011-03-26 19:10  广义霍夫变换\2.jpg
     文件       20331  2011-03-27 19:15  广义霍夫变换\3.jpg
     文件        1631  2011-03-27 19:16  广义霍夫变换\4.jpg
     文件        1051  2011-03-27 19:26  广义霍夫变换\5.jpg
     文件        1850  2011-03-27 19:28  广义霍夫变换\6.jpg
     文件       15243  2011-03-20 17:23  广义霍夫变换\7.jpg
     文件        2309  2011-03-20 17:24  广义霍夫变换\8.jpg
     文件        1086  2011-03-26 19:14  广义霍夫变换\ght.m
     文件         833  2011-03-28 23:59  广义霍夫变换\main.m
     文件        1056  2011-03-27 21:14  广义霍夫变换\main1.m
     文件       39831  2011-03-27 16:27  广义霍夫变换\noise.jpg
     文件        6731  2011-03-27 16:29  广义霍夫变换\original.jpg
     文件        1059  2011-03-26 17:36  广义霍夫变换\rtable.m
     文件        2479  2011-03-26 17:48  广义霍夫变换\template1.jpg
     文件        2335  2011-03-27 15:13  广义霍夫变换\template2.jpg
     文件        2265  2011-03-27 15:17  广义霍夫变换\template3.jpg