点检测

%点检测
f=imread('single_pixel.tif');%读取图像
subplot(1,2,1),imshow(f),title('原图像')

w=[-1,-1,-1;-1,8,-1;-1,-1,-1];%点检测掩模
g=abs(imfilter(double(f),w));%计算滤波图像
T=max(g(:));%使用滤波图像的数据找到T值
g=g>=T;%将已滤波图像与T比较
subplot(1,2,2),imshow(g),title('检测到的点')