matlab 实现数字图像的傅立叶变换及滤波锐化

最新推荐文章于 2025-10-16 14:17:04 发布

原创

最新推荐文章于 2025-10-16 14:17:04 发布 · 9.9k 阅读

55 ·

CC 4.0 BY-SA版权

文章标签：

#matlab #傅立叶变换 #锐化滤波

本文详细介绍了如何使用MATLAB进行数字图像处理，包括启动MATLAB程序读取图像并进行傅立叶变换，展示原始图像及其频谱图。接着，文章探讨了如何在图像上添加椒盐和高斯噪声，并通过FFT观察噪声影响。最后，通过高斯低通、巴特沃兹低通、高斯高通和巴特沃兹高通滤波器对比分析了图像的锐化和平滑效果。

转载请注明来自我的优快云博客：黄朝辉的博客

1. 启动MATLAB程序，读入一幅图像；对图像做FFT。使用’subplot’命令，同时显示原始图像其频谱图；

IenaImg=imread('lena.jpg');   %读入原图像文件
fftI=fft2(double(IenaImg));        %二维离散傅立叶变换
sfftI=fftshift(fftI);   %直流分量移到频谱中心
RR=real(sfftI);     %取傅立叶变换的实部
fftImg = log(RR);
subplot(1,2,1)
imshow(IenaImg);        %显示原图像
subplot(1,2,2)
imshow(fftImg, [8,10]);       %显示原图像的频谱

2. 读入一幅图像，分别为图像添加椒盐、高斯噪声，做FFT变换。使用’subplot’命令，将原始图像、原始图像频谱图、添加噪声后的图像，以及噪声图像的频谱图同时显示出来。

lenaImg=imread('lena.jpg');   %读入原图像文件
fftI=fft2(double(lenaImg));        %二维离散傅立叶变换
sfftI=fftshift(fftI);   %直流分量移到频谱中心
RR=real(sfftI);     %取傅立叶变换的实部
fftLenaImg = log(RR);

加入椒盐躁声

saltImg=imnoise(lenaImg,'salt & pepper',0.02); %加入椒盐躁声
fftI=fft2(double(saltImg));        %二维离散傅立叶变换
sfftI=fftshift(fftI);   %直流分量移到频谱中心
RR=real(sfftI);     %取傅立叶变换的实部
fftSaltImg = log(RR);

subplot(2,2,1);
imshow(lenaImg);
subplot(2,2,2);
imshow(fftLenaImg , [8,10]);
subplot(2,2,3);
imshow(saltImg);
subplot(2,2,4);
imshow(fftSaltImg , [8,10]);

加入高斯躁声

最低0.47元/天解锁文章