Harris角点检测matlab代码.doc

Harris角点检测代码 %%%Prewitt Operator Corner Detection.m%%%时间优化--相邻像素用取差的方法%%clear;for nfigure=1:6????t=input(input your figure’s name(including its extern name):,s); % t1 = tic;??????????????????????? %测算时间FileInfo = imfinfo(t);???????????? % 保存图像的所有信息Image = imread(t);???????????????? % 读取图像% 转为灰度值图像(Intensity Image)if(strcmp(truecolor,FileInfo.ColorType) == 1) %转为灰度值图像Image = im2uint8(rgb2gray(Image));???end????dx = [-1 0 1;-1 0 1;-1 0 1];? %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得到A?B = 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;for i=boundary:nrow-boundary+1???? for j=boundary:ncol-boundary+1??????? nlike=0; %相似点个数??????? if Image(i-1,j-1)Image(i,j)-t Image(i-1,j-1)Image(i,j)+t???????????? nlike=nlike+1;??????? end??????? if Image(i-1,j)Image(i,j)-t Image(i-1,j)Image(i,j)+t????????????? nlike=nlike+1;??????? end??????? if Image(i-1,j+1)Image(i,j)-t Image(i-1,j+1)Image(i,j)+t????????????? nlike=nlike+1;??????? end????????? if Image(i,j-1)Image(i,j)-t Image(i,j-1)Image(i,j)+t????????????? nlike=nlike+1;??????? end??????? if Image(i,j+1)Image(i,j)-t Image(i,j+1)Image(i,j)+t????????????? nlike=nlike+1;??????? end??????? if Image(i+1,j-1)Image(i,j)-t Image(i+1,j-1)Image(i,j)+t????????????? nlike=nlike+1;??????? end??????? if Image(i+1,j)Image(i,j)-t Image(i+1,j)Image(i,j)+t???????