Matlab牛刀小试

读入matlab 自带图像cameraman.tif，添加椒盐噪声后，使用fspecial 函数获取高斯空间滤波器，使用freqz2 函数将获得空间滤波器转换为频率域滤波器，然后使用dftfilt 函数进行频率域滤波，要求结果显示读入的原图像、空间滤波后结果及频率域滤波后结果。

% 读入 cameraman.tif 图像

img = imread('cameraman.tif');

% 显示原始图像

figure(1)

imshow(img);

title('Original Image')

% 添加椒盐噪声

noisy_img = imnoise(img, 'salt & pepper', 0.1);

% 显示添加噪声后的图像

figure(2)

imshow(noisy_img);

title('Noisy Image')

% 使用 fspecial 函数获取高斯空间滤波器

gaussian_filter = fspecial('gaussian', [5 5], 2);

% 显示高斯空间滤波器

figure(3)

surf(gaussian_filter);