Harris角点检测matlab代码

本文格式为Word版，下载可任意编辑——Harris角点检测matlab代码Harris角点检测代码

%%%PrewittOperatorCornerDetection.m%%%时间优化--相邻像素用取差的方法%%clear;

fornfigure=1:6

t=input('inputyourfigure’sname(includingitsexternname):','s');

%t1=tic;%测算时间

FileInfo=imfinfo(t);%保存图像的所有信息Image=imread(t);%读取图像%转为灰度值图像(IntensityImage)

if(strcmp('truecolor',FileInfo.ColorType)==1)%转为灰度值图像Image=im2uint8(rgb2gray(Image));end

dx=[-101;-101;-101];%dx：横向Prewitt差分模版Ix2=filter2(dx,Image).^2;

Iy2=filter2(dx',Image).^2;Ixy=filter2(dx,Image).*filter2(dx',Image);

%生成9*9高斯窗口。窗口越大，探测到的角点越少。h=fspecial('gaussian',9,2);

A=filter2(h,Ix2);%用高斯窗口差分Ix2得到AB=filter2(h,Iy2);C=filter2(h,Ixy);nrow=size(Image,1);ncol=size(Image,2);

Corner=zeros(nrow,ncol);%矩阵Corner用来保存候选角点位置,初值全零，值为1的点是角点

%真正的角点在137和138行由(row_ave,column_ave)得到%参数t:点(i,j)八邻域的“相像度〞参数，只有中心点与邻域其他八个点的像素值之差在%（-t,+t）之间，才确认它们为相像点，相像点不在候选角点之列

t=20;

%我并没有全部检测图像每个点，而是除去了边界上boundary个像素，%由于我们感兴趣的角点并不出现在边界上boundary=8;

fori=boundary:nrow-boundary+1forj=boundary:ncol-boundary+1nlike=0;%相像点个数

ifImage(i-1,j-1)>Image(i,j)-tend

ifImage(i-1,j)>Image(i,j)-tend

ifImage(i-1,j+1)>Image(i,j)-tend

ifImage(i,j-1)>Image(i,j)-tend

ifImage(i,j+1)>Image(i,j)-tend

ifImage(i+1,j-1)>Image(i,j)-tend

ifImage(i+1,j)>Image(i,j)-tend

ifImage(i+1,j+1)>Image(i,j)-tend

ifnlike>=2%假使周边有0，1，7，8个相像与中心的（i,j）%那(i,j)就不是角点，所以，直接忽略end;

end;end;

CRF=zeros(nrow,ncol);%CRF用来保存角点响应函数值,初值全零CRFmax=0;%图像中角点响应函数的最大值，作阈值之用t=0.05;%计算CRF

%工程上常用CRF(i,j)=det(M)/trace(M)计算CRF，那么此时应当将下面第105行的%比例系数t设置大一些，t=0.1对采集的这几幅图像来说是一个比较合理的经验值fori=boundary:nrow-boundary+1forj=boundary:ncol-boundary+1ifCorner(i,j)==1%只关注候选点M=[A(i,j)C(i,j);C(i,j)B(i,j)];

CRF(i,j)=det(M)-t*(trace(M))^2;ifCRF(i,j)>CRFmaxCRFmax=CRF(i,j);end;endend;end;%CRFmax

count=0;%用来记录角点的个数t=0.01;

%下面通过一个3*3的窗口来判断当前位置是否为角点fori=boundary:nrow-boundary+1forj=boundary:ncol-boundary+1

ifCorner(i,j)==1%只关注候选点的八邻域

ifCRF(i,j)>t*CRFmax%这个是角点，count加1

else%假使当前位置（i,j）不是角点，则在Corner(i,j)中删除对该候选角点的记录Corner(i,j)=0;end;end;end;end;

%disp('角点个数');%disp(count)

figure,imshow(Image);%displayIntensityImageholdon;%toc(t1)

fori=boundary:nrow-boundary+1forj=boundary:ncol-boundary+1column_ave=0;row_ave=0;k=0;

ifCorner(i,j)==1

forx=i-3:i+3%7*7邻域fory=j-3:j+3ifCorner(x,y)==1

%用算数平均数作为角点坐标，假使改用几何平均数求点的平均坐标，对角点的提取意义不大

row_ave=row_ave+x;column_ave=column_ave+y;k=k+1;