数字图象处理实验内容.docVIP

实验一 数字图像的傅立叶变换 一、 实验目的 （一）用快速傅立叶变换FFT实现数字图像的傅立叶变换，进一步加深对DFT算法原理和基本性质的理解（因为FFT只是DFT的一种快速算法，所以FFT的算法结果必然满足DFT的基本性质）。 （二）熟悉FFT算法原理和FFT子程序的应用。 （三）学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分布误差及其原因，以便在实际中正确应用FFT。 二、 实验原理和方法 （一）傅立叶（Fourier）变换的定义 对于二维信号，二维连续傅立叶变换定义为： 正变换： 反变换： 二维离散傅立叶变换为： 正变换： 反变换： 图像的傅立叶变换与一维信号的傅立叶变换变换一样，有快速算法。实际上，现在有实现傅立叶变换的芯片，可以实时实现傅立叶变换。 （二） 其逆变换为 离散余弦变换是图像压缩中常用的一个变换编码方法，任何实对称函数的傅立叶变换中只含余弦项，就成为余弦变换，因此余弦变换是傅立叶变换的特例。余弦变换与傅立叶变换一样有明确的物理意义，是简化傅立叶变换的重要方法。 （三）矩阵形式的傅立叶变换的算法如下： 数字图像f的傅立叶正变换： 数字图像f的傅立叶反变换： 变换矩阵，其中，N 为图像的维数。 三、 实验结果及分析 （一）对原图像进行傅立叶变换 分析： （二）输出彩色图像greens.jpg的傅立叶频谱： 分析： （三）对彩色图像football.jpg进行二维DCT变换， 分析： 四、 思考题 （一）在数字图像傅立叶变换的频谱图上原点附近出现的亮点的来源是什么？ （二）如何在数字地图（或普通景物的数字图像）的频谱图上识别地物（或类别）的延伸方向？ 利用Matlab语言编写的数字图像处理的例程如下： 傅立叶变换Matlab图像的DFT clc; figure 1 ; load imdemos saturn2; imshow saturn2 ; title 原图像 ; figure 2 ; S fftshift fft2 saturn2 ; imshow log abs S ,[] ; title 原图像傅立叶频谱 ; 彩色图像的傅立叶频谱 figure 1 ; A imread greens.jpg ; B rgb2gray A ; imshow B ; title 原图像 ; S fftshift fft2 B ; figure 2 ; imshow log abs S ,[] ; title 彩色图像的傅立叶频谱 ; 二维DCT变换 RGB imread football.jpg ; figure 1 ; imshow RGB ; title 彩色图像 ; GRAY rgb2gray RGB ; figure 2 ; imshow GRAY ; title 灰色图像 ; DCT dct2 GRAY ; figure 3 ; imshow log abs DCT ,[] ; title 二维DCT变换 ; fftshift 作用：将零频点移到频谱的中间 用法： Y fftshift X Y fftshift X，dim 描述：fftshift移动零频点到频谱中间，重新排列fft,fft2和fftn的输出结果。将零频点放到频谱的中间对于观察傅立叶变换是有用的。 示例： clf; fs 100;N 256;?? %采样频率和数据点数n 0:N-1;t n/fs;?? %时间序列x 0.5*sin 2*pi*15*t +2*sin 2*pi*40*t ; %信号 y1 fft x,N ;??? %对信号进行快速Fourier变换y2 fftshift y1 ; mag1 abs y1 ; %求得Fourier变换后的振幅mag2 abs y2 ; f1 n*fs/N;??? %频率序列f2 n*fs/N-fs/2; subplot 3,1,1 ,plot f1,mag1,r ;?? %绘出随频率变化的振幅xlabel 频率/Hz ;ylabel 振幅 ;title 图1：usual FFT,color,r ;grid on; subplot 3,1,2 ,plot f2,mag1,b ;?? %绘出随频率变化的振幅xlabel 频率/Hz ;ylabel 振幅 ;title 图2：FFT without fftshift,color,b ;grid on; subplot 3,1,3 ,plot f2,mag2,c ;?? %绘出随频率变化的振幅xlabel 频率/Hz ;ylabel 振幅 ;title 图3：FFT after fftshift,color,c ;grid on; 结论： 1）如果期望绘制的幅频图的频率范围为0~fs