1、。SUSAN角点检测算子的MATLAB实现filename,pathname,=uigetfile(cnt=0;%计数专用，统*.jpg, 选择 JPG格式图片 );计圆形邻域内满足条件的像素点个数if ischar(filename)for p=1:7returnfor q=1:7endifstr=pathname filename;(p-4)2+(q-4)2=12% 半径pic=imread(str);一般在 34 之间if length(size(pic)=3ifimg=rgb2gray(pic);abs(img(i,j)-tmp(p,q)tendM,N=size(img);cnt=cn
2、t+1;timg=zeros(M+6,N+6);endtimg(4:end-3,4:end-3)=img;end%扩展图像边缘 3个像素endimg=timg;endt=45;%阈值USAN=USAN cnt;USAN=;%用于存放 USANendfor i= 4:M+3endfor j=4:N+3tmp=img(i-3:i+3,j-3:j+3);g=max(USAN)/2;% 给定的阈-可编辑修改 -。值for i=2:M-1for k=1:length(USAN)for j=2:N-1if USAN(k)sq%USAN(k)=0;局部非极大值抑制endloc=loc;j,i;endendendimgn=reshape(USAN,M,N);%endUSAN 向量张成二维图像imshowpair(pic,pic,montage);imgn=fliplr(imrotate(imgn,-90);hold on%调整图像plot(loc(:,1)+size(pic,2),loc(:,2),*);hold offloc=;-可编辑修改 -。运行结果图如下：-可编辑修改 -
