周期噪声频域滤除、散粒噪声中值滤波

对周期噪声，使用频域方式滤波。使用FFT将其变换到频谱域对恒定频率的信号滤除。

对散粒噪声，使用中值滤波。设定5x5的模板统计排序，选取其中值替代模板的中心位置的值。

 

I=imread('D:\文件及下载相关\图片\child.jpg');
G = rgb2gray(I);
[M,N]=size(G);
F = fft2(G);
F1 = fftshift(F);
F2 = log(1+abs(F1));
x = 0:(M-1);
y = 0:(N-1);
[Y,X] = meshgrid(y,x);
Q0 = 10;
y0 = 225;
x0 = 272; 
y1 = 289;
x1 = 272;  
Q1 = sqrt((Y-y0).^2+(X-x0).^2);
Q2 = sqrt((Y-y1).^2+(X-x1).^2);
H = 1./(1+(Q0^2./(Q1.*Q2)).^8);
I =abs(real(ifft2(OUT1.*H)));
figure,imshow(I,[])

 

中值滤波：

clc
I = imread('C:\Users\Desktop\椒盐噪声.jpg');
G = rgb2gray(I);
G1=medfilt2(G,[5,5]);
figure()
subplot(1,2,1);
imshow(G1);
subplot(1,2,2);
imshow(G)