数字图像处理dc.docxVIP

《数字图像处理》实验 报 告 学院：信息与通信工程学院 班级: 学号： 姓名： 实验一 图像变换 实验目的 掌握图像变换的方法； 掌握离散二维余弦变换的原理和性质； 掌握离散二维傅里叶变换的原理和性质。 实验原理 1. 1-D FFT 序列的DFT定义为： （1） 设，，，则上式变为： （2） （3） （4） 将序列按序号的奇偶分成两组，即： （5） 因此，的傅里叶变换可以写成： （6） 由此可得： （7） 式中： （8） （9） 它们分别是和的点DFT。上面的推导表明：一个点的DFT被分解为两个点的DFT，这两个点的DFT又可合成一个点的DFT。但上面的公式仅能得到的前点的值，运用权系数的周期性与对称性，即： ， （10） 因此，的后点的值可表示为： （11） 通过上面的推导可以看出，一个点的DFT可以分解为两个点的DFT，每个点的DFT又可以分解为两个点的DFT。以此类推，当为2的整数次幂时（），由于每分解一次降低一阶幂次，所以通过次分解，最后全部成为一系列2点的DFT运算，以上就是按时间抽取的FFT算法。 2. 2-D FFT 以图形为例，将式（1）推广，可得2-D FFT式： （12） 将式（12）写成如下的分离形式： （13） 由上述分离形式可知一个2-D DFT可由连续2次运用1-D DFT来实现。如式（13）可分成下列两式： （14） （15） 对每个值，式（14）方括号中是一个1-D DFT。所以可由沿的每一列求变换得到。在此基础上，再对的每一行求变换就可得到。由上述可知，对于2-D FFT同样可以采取分别进行2次1-D FFT来得到。 3. 数字图像FFT流程图 根据图像阵列的特性，对其按下列步骤进行FFT变换： （1）将图像数据阵列变换为按列存储，即从下到上，从左到右； （2）对每一列图像数据进行1-D FFT； （3）将按列处理后的数据结果存储，并对它们按原图像阵列的形式，即按行重新存储； （4）对重新排列的数据逐行进行1-D FFT； （5）将每行处理后的数据存储起来，即得到2-D 数字图像FFT结果。 对一幅图像进行FFT其流程图如下所示： 图1 2-D FFT流程图 其中的1-D FFT需要倒位序和递推计算，它们的流程图如下： 图2 1-D FFT流程图 1-D FFT中整个L级递推过程由三个嵌套循环构成。外层的一个循环控制L(L=)级的顺序循环；内层的两个循环控制同一级各蝶形结的运算，其中最内一层循环控制同一种（指的同一个r）蝶形结的运算，而中间一层循环则是对不同种（即中的r不同）蝶形结做运算。 4. 2-D离散傅里叶反变换 （16） 5. 1-D离散余弦变换 （17） （18） 其中由下式定义： 6. 2-D离散余弦变换 （19） （20） 实验内容和步骤 实验内容： 对标准图像进行离散傅里叶变换并在计算机屏幕观测其频谱，验证二维傅里叶变换的常用性质。 对标准图像进行离散余弦变换并在计算机屏幕观测其频谱，验证二维余弦变换的常用性质，了解二维余弦变换用在图像压缩中的原因。 对标准图像离散傅里叶变换和离散余弦变换的频谱进行比较。 实验步骤： 学习Matlab语言的常用和图像处理函数库，会使用其函数进行编程。 用Matlab语言、C语言或C++语言编制标准图像离散傅里叶变换和离散余弦变换的程序。 将处理结果存盘。 实验设备 计算机、Matlab软件、Photoshop软件。 程序源代码及所得对应图像 %DFT频谱和DCT频谱比较 clc clear %创建矩形函数 f=zeros(256,256); f(180:240,130:150)=1; figure; imshow(f); %DFT2 F=fft2(f,256,256); figure; imshow(log(abs(F)),[]); %将直流分量移到中心 F2=fftshift(F); figure; imshow(log(abs(F2)),[]); %DCT2 C=d