caicai.txt

对图像添加盐噪

function [N,M] = getimage(u)
clear, close all;
I = imread('sun.Jpg');
figure,imshow(I);
I2 = imnoise(I,'salt');%加入盐噪声
figure,imshow(I2);
f = double(I2);
g = fft2(f);%傅立叶变换  
g = fftshift(g);%转换数据矩阵
[N,M] = fun(g);
 result = ifftshift(result);
    X2 = ifft2(result);
    J = uint8(real(X2));
    figure,imshow(J);

function [N,M] = fun(binImage)
[N1,N2]=size(binImage);
n = 2;
d0 = 80;%d0为截止频率
n1 = fix(N1/2);
n2 = fix(N2/2);
for i = 1:N1;%对频域中循环滤波
    for j = 1:N2;
        d = sqrt((i-n1)^2+(j-n2)^2);
        h = 1/(1+0.414*(d/d0)^(2*n));
        result(i,j) = h*binImage(i,j);
                end
end