灰階圖像的頻率域濾波——拉普拉斯——高頻提升（Matlab）

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);

%%%%%%%================高頻提升（拉普拉斯算子）============================
%參數A>=1,當其等于1時，為普通的高通濾波器
A=2
for x=1:M
    for y=1:N
        %計算頻率域拉普拉斯算子
        h_hp=1+4*((x-M0)^2+(y-N0)^2)/(M0*N0);
        h_bp=(A-1)+h_hp;
        J_shift(x,y)=J_shift(x,y)*h_bp;
    end
end

J=ifftshift(J_shift);
I_D_rep=ifft2(J);
imshow(I_D_rep)