【matlab】图像去噪的代码测试

%% 自己设置频域的滤波窗口
girl=imread('F:\Users\*****\Pictures\CGS_stripe1.bmp');
girl=rgb2gray(girl);
girl=im2double(girl);
subplot(,,);imshow(girl);%显示原始图像
title('原始图像');
G=fft2(double(girl));% 二维傅里叶变换
G=fftshift(G);%象限转换
subplot(,,);imshow(log(abs(G)+),[ ]);%显示原始图像频谱图
title('原始图像频谱图');figure;
%构造理想低通滤波器，并用它滤波
[height  width]=size(G);
H(: height,: width)=0.9;
x0= height/; y0= width/;%确定原点
for x=:height
    for y=:width
        if(sqrt((x- x0)*(x- x0)+(y-y0)*(y-y0))>) %确定滤波半径
            H(x,y)=;
        end
        F(x,y)=G(x,y)*H(x,y);
    end
end
imshow(log(abs(F)+),[ ]);%显示中间图像频谱图
title('中间图像频谱图');
g=ifft2(F);% 傅里叶反变换
figure;imshow(abs(real(g)),[ ]);%显示处理后的图像
title('处理后的图像');