clc;
clear all;
close all;
%读入图像,并转换为double型
I=imread('D:\Gray Files\4-41.tif');
I_D=im2double(I);
%获得图像的高度和宽度
[M,N]=size(I_D);
%图像中心点
M0=M/2;
N0=N/2;
J=fft2(I_D);
J_shift=fftshift(J);
%%%==================频率域拉普拉斯算子,高通滤波==========================
for x=1:M
for y=1:N
%计算频率域拉普拉斯算子
h_hp=1+4*((x-M0)^2+(y-N0)^2)/(M0*N0);
J_shift(x,y)=J_shift(x,y)*h_hp;
end
end
J=ifftshift(J_shift);
I_D_rep=ifft2(J);
imshow(I_D_rep,[]);