MATLAB实现HARRIS算法

在MATLAB实现Harris角点检测算法。

function T1=coortransf2（T0f）
%%%%图像坐标系到柱面坐标的转换
%对输入图像序列T0按焦距f进行变换
[HWrN]=size（T0）;
w2f=W/2/f;
h2=H/2;
constant2=f*atan（W/（2*f））;
constant1=h2;

for y=1:W      %%%%%%   列
    angle=atan（y/f-w2f）;%%%%atan（（y-W/2）/f）;
     y1=uint16（f*angle+constant2）;
     if  y1==0 
         y1=1; 
     end
    for x=1:H  %%%%%%%%%%%%   行
     x1=uint16（（x-h2）*cos（angle）+constant1）;
     if  x1==0  
         x1=1;
     end
     if r==3   %%%%%%%%%%%%%彩色图像
         for n=1:N  %%%%%%%%%
             if （y==coor_shift（n2）） 
                 coor_shift02（n2）=y1; 
             end%%%对应偏移量
             T1（x1y1:n）=T0（xy:n）;%%%%点的映射
         end
     elseif r==1
          T1（x1y1）=T0（xy）;
     end
    end
end
[hwaN]=size（T1）;
for i=1:60%%%边缘没有映射点的地方填白色
    for j=1:

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     253299  2009-02-21 14:58  harris\X2.jpg

     文件       1016  2009-02-17 21:32  harris\coortransf2.m

     文件       1171  2008-08-08 10:08  harris\harris.m

     文件       1187  2009-02-22 12:44  harris\Harris_detector.m

     文件     241430  2009-02-21 14:58  harris\X1.jpg

     目录          0  2009-05-19 20:41  harris

----------- ---------  ---------- -----  ----

               498103                    6