在路上走的時候，突然想到前面發布的微博忘了考慮y方向的頻譜，特此更正。
周期噪聲一般產生於圖像采集過程中的電氣或電機的干擾，表現為圖像中周期性的沖擊，如下圖
噪聲的頻率通過截取過頻率零點的橫截線來確定。下圖為上一幅圖的x方向的橫截線，
從圖中可以看出，在x=113和145處有明顯的頻率增強，這就是周期噪聲的頻率。
周期噪聲的去除一般通過頻域濾波實現，如構建一個巴特沃茲陷波器進行頻域濾波。構建的巴特沃茲陷波器如下
圖中最黑處為零，最白處為1。
頻域濾波后，把圖像從頻域轉換到空域即可，下圖就是濾波后圖像。
整個Matlab代碼如下：
[M,N,~] = size(f); %f為含噪圖像
F = fftshift(fft2(f)); % 圖像頻域
figure;plot(abs(F(M/2+1,:))); % x方向橫截線
% 構造n階巴特沃茲陷波器，可以修改D0和n獲得更好的效果。
u = 0:(M-1);
v = 0:(N-1);
[V,U] = meshgrid(v,u);
D0 = 10;
n = 2;
v0 = 113;
u0 = 127; %這張圖y方向沒有周期噪聲
v1 = 145;
u1 = 127;
D1 = sqrt((V-v0).^2+(U-u0).^2);
D2 = sqrt((V-v1).^2+(U-u1).^2);
H = 1./(1+(D0^2./(D1.*D2)).^n);
figure;imshow(H,[])
% 濾波
G = fftshift(F.*H);
g = real(ifft2(G));
figure;
imshow(g,[]);
更多知識可以參考岡薩雷斯的數字圖像處理（matlab版）。