点检测原理：
function PointDetect(filename,number)
%{
%%%%%%%%%%%%%%%%%%%%%%%%%
%参数：
%filename 被检测图像文件名
%number 检测点数量
%功能：
%从被检测图像中检测出指定数量个特征点
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%}
f=imread(filename);%输入图像
figure,imshow(f)
[width,height]=size(f);%球图像大小
h=zeros(width,height);%构造跟原图一样大小的0矩阵
%如果特征点数目大于图像的像素点，则特征点数就是图像大小
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);%用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)
