图像处理单一频点周期噪声matlab

以下是使用MATLAB进行图像处理单一频点周期噪声的方法：

1.读取图像并显示

img = imread('swanNoise.jpg');
imshow(img);

2.计算傅里叶变换并将零频分量移到中心

f = fft2(img);
fshift = fftshift(f);

3.生成滤波器并将其应用于频域图像

[M, N] = size(f);
D0 = 50; % 滤波器半径
H = zeros(M, N);
for i = 1:M
    for j = 1:N
        D = sqrt((i-M/2)^2 + (j-N/2)^2);
        if D <= D0
            H(i, j) = 1;
        end
    end
end
g = fshift .* H;

4.计算逆傅里叶变换并显示结果

gshift = ifftshift(g);
result = ifft2(gshift);
result = uint8(real(result));
imshow(result);