傅里叶变换离散余弦变换与小波变换

1、数字图像处理实验报告二维离散傅里叶、余弦、小波变换专业班级：10信息安全学生姓名：王猛涛学生学号：20101616310049指导教师：姚孝明完成时间：2013年4月13日14数字图像处理实验三：二维离散傅里叶、余弦、小波变换、实验目的1. 了解图像正变换和逆变换的原理2. 了解图像变换系数的特点。3. 掌握常用图像变换的实现过程。4. 掌握图像的频谱分析方法。5. 了解图像、实验主要仪器设备1. 微型计算机：In tel Pen tium及更高2. MATLAB 软件。三、实验原理二维离散傅里叶变换、余弦变换、小波变换的正逆变换公式，MATLAB中的上述变换的实现函数以及讨论正交变换的应用。

2、1.二维离散傅里叶变换(Discrete Fourier Transform, DFT)对于二维傅立叶变换，其离散形式如式(1)所示；逆变换公式如式(示:F(u,v)二MNMNx £ y =0f (x, y)e_j 2 ： (ux/ M : !vy/N )M 4N f(x, y)二二 F(u,v)ej2：(ux/M ”N) u 30 v 30频谱公式如式(3)所示：F(u,v) =|F(u,v)|ej (u,v) =R(u,v) jI (u,v)|F(u,v)|= ,R2(u,v) I2(u,v)2)所(1)(2)(3)由可傅立叶变换的分离性可知，一个二维傅立叶变换可分解为两步进行，

3、 其中每一步都是一个一维傅立叶变换。先对f(x, y)按列进行傅立叶变换得到F(x,v)，再对F(x, v)按行进行傅立叶变换，便可得到f(x, y)的傅立叶变换结果。显然对f(x, y)先按行进行离散傅立叶变换，再按列进行离散傅立叶变换也是可行的, 这里不再一一赘述此外，在实际工程应用中分析幅度谱较多，习惯上也常把幅度谱称为频谱。 使用DFT变换进行图像处理时，有如下特点：1 M 4N 二(1) 频谱的直流成分为f(o,o)= 2二:二：f (x,y)，说明在频谱原点的傅里N x=0 y=0叶变换F (0,0)等于图像的平均灰度级。(2) 幅度谱|F(u,v) |关于原点对称，即F(u,v)

4、二F(_u,_v)。(3) 图像f(x,y)平移后，幅度谱不发生变化，仅有相位发生变化DFT是一种基本和重要的正交变换。为了提高计算效率，应用时往往采用二维FFT实现。而一般的正交变换图像经过对数变换后便于观察。MATLAB采用fft2和ifft2分别进行二维DFT变换和二维DFT逆变换，采用fftshift将直流 分量移到频谱图的中心以便观察。2.二维离散余弦变换(Discrete Cosine Transform DCT)对于二维余弦变换，其离散形式如式(4)所示，逆变换如式XJ(x,y)coSMNu(T)coSMNv(y1 v = 0F(u,v)=C(u)C叫而式中，C(u)4吉，u=&

5、#176;1,1 Eu EM -1_2 M dN Jf(x, y) 加；二 C(u)C(v)F(u,v)co、MN um心(5)所示：弓(4)C(v) = * J2 '1,1 兰 v 兰 N 12)coSNv(y4)兀SMu(x(5)在MATLAB中，采用dct2和idct2分别进行二维 DCT变换和二维DCT逆变换。二维DCT常用于信号和图像处理，典型应用是对静止图像和运动图像进行性能优良的有损数据压缩。在静止图像编码标准JPEG、运动图像编码标准MJPEG和MPEG等标准中都使用了 8*8块的离散余弦变换，并将结果进行量化 之后进行熵编码。DCT具有很强的能量集中在频谱的低频部分的

6、特性，而且当 信号具有接近马尔科夫过程(Markov processes)的统计特性时，DCT的去相关 性接近于具有最优去相关性的 K-L变换(Karhunen-Loeve变换)的性能。另外，改进的离散余弦变换(Modified Discrete Cosine Transform， MDCT) 对交叠的数据进行 DCT，有助于避免由于区块边界所产生的多余数据，被用在 高级音频编码(Advaneed Audio Coding，AAC )、Ogg Wrbis、AC3 和 MP3 音 频压缩中。3. 二维离散小波变换(D Discrete Space Wavelet Transform DDSWT)

7、对于二维小波变换，其离散形式如式(6)所示；逆变换如式(7)所示：Wf (a,bx,by< f ? a,bxby(x, y) = - . r2 f (x, y *(x b， b)dxdy(6)aa a式中，bx和by分别函数f (x, y)在轴上的x， y平移量。f （x, y 12 R R2Wf （a, bx,by）-a,b（x bxby ）dbxdbyda（7）aa a类似地，可以定义二维离散小波变换逼近，并采用Mallat二维快速算法求解。与DFT类似，可分离二维小波变换最终可转换为两次一维小波变换。对图像进行小波变换的MATLAB常用函数有： 对图像进行一层二维小波分解，常见形式

8、为：CA,CH,CV,CD=dwt2（X, wname'式中，X为图像矩阵；wname'是使用的小波基函数名称，如可选择双正交 样条小波基函数，形式为biorNr.Nd。 查询使用的小波基函数的信息，使用形式为：Waveinfo（ Wname'）式中，小波基名称'wname'可选用'haar'（哈尔小波）、）b）（Daubechies小 波）、bior '（双正交样条小波）等。例如，在命令行状态下键入wavainfo（ bior）进行查询双正交样条小波，可知r表示reconstruction（重建），d表示 decompositi

9、on （分解），N表示相应FIR滤波器的阶数；CA、CH、CV、CD分别是输入矩阵 X小波分解的近似系数矩阵、水平细节系数、垂直细节系数和对角线细节系数。 对二维小波分解的图像进行各种分量的重构，常见函数形式为：Y=upcoef2（O,X, wname',N）式中，X是分解后的细节信号，丫是重构的细节信号分量；N表示对矩阵X 的系数进行重建的步骤数，即重构的层数，默认值为1。O是细节信号的类型。如果O='a',则表示对信号的近似系数进行重建；否则，如果O='h'、v'或''则分别对水平、垂直或对角线细节进行重建。 对应上述的一层二

10、维小波变换 DWT2函数，进行一层二维小波变换逆变 换，常见形式为：X=idwt2（CA,CH,CV ,CD, 'wname'idwt2函数采用wname'所指示的小波、已重建的基于近似矩阵 CA，以及水 平细节CH、垂直细节CV和对角线细节CD计算原图像矩阵X。 对重构的图像进行量化编码，常见函数形式为：Y=wcodemat（X,NBCODES,OPT,ABSOL）式中，X为待进行量化编码的矩阵，丫为编码矩阵。在编码中，把矩阵 X 中元素绝对值最大的作为 NBCODES （整数），绝对值最小的作为1,其他元素 依其绝对值的大小在1与NBCODES中排列。当OPT为&#

11、39;row'时，做行编码；当OPT为'col'时，做列编码；当OPT为'mat'时，做全局编码，即把整个矩阵中元 素绝对值最大的元素作为 NBCODES，最小的作为1.当ABSOL为0时，该函数 返回输入矩阵X的一个编码版本，当ABSOL非0时，返回X的绝对值。四、实验内容1. 在MATLAB环境中，进行图像的离散傅里叶变换和离散余弦变换，观察 图像的频谱并减少DCT系数，观察重建信号和误差信号，理解正交变换在压缩 编码中的应用。2. 在MATLAB环境中，进行图像的离散小波变换，观察图像的近似图像和 各方向的细节图像，观察重建图像，理解小波变换在图像

12、特征检测（如边缘检测、 方向检测等）中的应用。五、实验步骤1. 选择典型图像作为研究对象2. 显示原始图像。3. 进行图像变换4. 对图像进行处理（如选择不同个数的变换系数可以进行压缩、 选择不同方 向的频谱可以进行特征检测等）。5. 对图像进行逆变换复原图像，观察重建图像和误差图像并进行对比六、实验结果及分析1. 二维离散傅里叶变换(Discrete Fourier Transform, DFT)根据实验内容要求，我们通过编写 MATLAB程序（详见附录一），调用函数fft2和ifft2，对所选用的图像（见下图1）进行二维离散傅里叶变换得到如下图1原始图像巨0西附EMl E图2实验1运行结果

13、截图从现代数学的眼光来看，傅里叶变换是一种特殊的积分变换。它能将满足一定条件的某个函数表示成正弦基函数的线性组合或者积分。在不同的研究领域，傅里叶变换具有多种不同的变体形式，如连续傅里叶变换和离散傅里叶变换。当 前傅里叶变换在物理学、数论、组合数学、信号处理、概率论、统计学、密码学、 声学、光学、海洋学、结构动力学等领域都有着广泛的应用（例如在信号处理中， 傅里叶变换的典型用途是将信号分解成幅值分量和频率分量）。2. 二维离散余弦变换（Discrete Cosine Transform，DCT）根据实验内容要求，我们通过编写 MATLAB程序（详见附录二），调用函 数dct2和idct2，对所

14、选用的图像（见下图 3）进行二维离散余弦变换得到如下结果（见图4）:图3原始图像图4实验2运行结果截图从实验结果中可以看出：离散余弦变换的重要特点是能量集中，信号常将其 能量的大部分集中于频率域的一个小范围内， 这样描述不重要的分量只需要很少 的比特数；频率域分解映射了人类感觉系统的处理过程，并允许后继的量化过程 满足其灵敏度的要求。变化后，能量集中的范围可以精细的量化， 其他的范围可 以粗糙量化，这样处理，不会引起太大的精度问题，符合人体的听觉，视觉需要。 这样量化后，可以用小的数据量来保存采集的数据， 对处理音频，视频等数据非 吊有效。此外，离散余弦变换，在当前音频视频编码中起着非常重要的

15、作用。3. 二维离散小波变换（D Discrete Space Wavelet Transform DDSWT）根据实验内容要求，我们通过编写MATLAB程序（详见附录三），调用函数dwt2和idwt2，对所选用的图像（见图5），进行二维离散小波逆变换得到图像（见图6）及一层小波变换后的图像（见图 7）,具体如下所示:逆变損后的團儒1DD1502DO50IM1»找2S0图6逆变换后的图像图5原始图像Approximati on A15010015020025050100150200250Vertical Detail V110015020025050100150200250Horiz

16、ontal Detail H15010015020025050100150200250Diago nal Detail D150100150200250o o O5 0 5o OO 52 2图7 一层小波变换的四个分量经过上述实验，从结果分析并结合查找的资料，我们可以得到小波变换具有 如下优点：（1） 小波分解可以覆盖整个频域（提供了一个数学上完备的描述）；（2）小波变换通过选取合适的滤波器，可以极大的减小或去除所提取得不 同特征之间的相关性；（3）小波变换具有“变焦”特性，在低频段可用高频率分辨率和低时间分辨率（宽分析窗口），在高频段，可用低频率分辨率和高时间分辨率（窄分析窗口）；（4）小波

17、变换实现上有快速算法（Mallat小波分解算法）。由于傅里叶变换分析处理信号的方法也存在着一定的局限性与弱点，傅里叶变换提供了信号在频率域上的详细特征，但却把时间域上的特征完全丢失了。 小 波变换是80年代后期发展起来的新数学分支，它是傅里叶变换的发展与扩充， 在一定程度上克服了傅里叶变换的弱点与局限性。本文从信号分析与处理的角度来介绍小波变换的基本理论与应用，使具有微积分基础的读者通过本文能对这一 新的数学分支有一初步了解。小波变换在函数论、微分方程、数值计算等方面也 有着重要的应用附录：附录一：l=imread('pout.tif);%subplot(2,2,1),imshow(l

18、);%title(' 原始图像');F1=fft2(I);%subplot(2,2,2),imshow(F1); % title('二维离散傅里叶变换后的图像');subplot(2,2,3),imshow(log(abs(F1)+1),0 10);%傅里叶变换读入图像在位置1显示图像计算二维傅立叶变换显示二维傅立叶变换后的图像%显示对数变换后的频谱图title(' 图像的频谱图');将直流分量移到频谱图的中心显示对数变换后中心化的频谱图计算傅立叶变换的逆变换重新显示图像，类似于imread（'pout.tif）读入图像显示原图像F2=f

19、ftshift(F1);%subplot(2,2,4),imshow(log(abs(F2)+1),0 10);%title('中心化的频谱图');F3=ifft2(F1);%figure,imshow(ui nt8(F3);%title('原始图像');附录二：炉维DCT的应用l=imread('kids.tif);%subplot(1,2,1),imshow(l);%title('原始图像');C仁dct2(l);%C2=fftshift(C1);%subplot(1,2,2),imshow(log(abs(C2)+1,0 10);

20、% title('DCT 系数');C3=idct2(C1);%figure,imshow(ui nt8(C3);%title('原始图像');附录三：%matlab图像的连续小波变换load wbarb;%image(X);%colormap(map);%cA1,cH1,cV1,cD1=dwt2(X,'bior3.7'); % 条小波对应的滤波器A1=upcoef2('a',cA1,'bior3.7',1);H仁 upcoef2('h',cH1,'bior3.7',1);V1=upcoef2('v',cV1,'bior3.7',