转载请注明来自我的CSDN博客:黄朝辉的博客
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]);
加入高斯躁声
ga