www.gusucode.com > 《图像配准技术及其Matlab编程实现》--源码程序 > 《图像配准技术及其Matlab编程实现》/match/PointDetect.m

    function PointDetect(filename,number)
%点检测
%参数：fileName：被检测图像的文件名
%     number：检测点的数量
%功能：从被检测图像中检测出指定数量个特征点。
f=imread(filename);
figure,imshow(f);
[width,height]=size(f);
h=zeros(width,height);
if number>width*height
    number=width*height;
end
df=im2double(f);
w=[-1 -1 -1;-1 8 -1;-1 -1 -1];
g=imfilter(df,w);
g=abs(g)./8;
[data index]=sort(g(:));
T=data(width*height - number +1);
for(i=1:width)
    for j=1:height
        if g(i,j)>=T
            h(i,j)=1;
            %h(i,j)=g(i,j);
        end
    end
end
figure,imshow(h)

%PointDetect('ct.bmp',1000)