matlab数字图像的傅立叶变换实验

实验三  图像的傅立叶变换


1.启动MATLAB程序，读入一幅图像；对图像做FFT。使用’subplot’命令，同时显示原始图像其频谱图；
1.1实验过程：
首先读取一幅图像，然后将这幅图像归一化到0~1之间，然后对图像做二维离散傅立叶变换，然后做快速傅里叶变换，即直流分量移到频谱中心，让正半轴部分和负半轴部分的图像分别关于各自的中心对称，然后取傅立叶变换的实部，然后做频谱对数变换，就得到原图像的频谱图，最后将原图像和对应的频谱图显示出来。
1.2代码及图像






I=imread('f:\lena.jpg'); %读入图像
F=fft2(im2double(I)); %FFT
F=fftshift(F); %FFT频谱平移
F=real(F);
T=log(F+1); %频谱对数变换
subplot(1,2,1),imshow(I),title('原始图像');
subplot(1,2,2),imshow(T,[]),title('原始图像其频谱图');
















 






2. 读入一幅图像，分别为图像添加椒盐、高斯噪声，做FFT变换。使用’subplot’命令，将原始图像、原始图像频谱图、添加噪声后的图像，以及噪声图像的频谱图同时显示出来。
2.1实验过程
得到原始图像的频谱图的方法同上，然后对原图像模拟叠加密度为0.04的椒盐噪声，然后添加模拟均值为0方差为0.02的高斯噪声，然后得到噪声图像的频谱图，方法同上，最后显示所有图像。
2.2代码及图像
I=imread('f:\lena.jpg&#