MATLAB图像处理编程与应用 作者张涛 第11章 图像正交变换.pptVIP

第11章 图像正交变换 问题的提出： 视觉所感受到的是在空间域和时间域的信号。但是，往往许多问题在频域中讨论时，有其非常方便分析的一面。 图像变换的目的： 使图像处理问题简化 有利于图像特征提取 有助于从概念上增强对图像信息的理解　　　 变换问题的引入 什么是图像变换 将图像看成是线性叠加系统 图像在空域上相关性很强 图像变换是将图像从空域变换到其它域如频域的数学变换 常用的变换： 傅立叶变换、沃尔什变换、哈达玛变换、离散余弦变换、离散K-L变换、小波变换 先对行做变换： Fourier 变换示意图 Fourier变换的频率特性 Fourier变换的低通滤波 Fourier变换的高通滤波 Fourier变换的压缩原理 Fourier变换的压缩原理 快速傅里叶变换 问题的提出： 离散傅里叶变换已成为数字信号处理的重要工具。然而，它的计算量较大，运算时间长，在某种程度上却限制了它的使用范围。 11.2 离散余弦变换（DCT） Fourier变换的一个最大的问题是：它的参数都是复数，在数据的描述上相当于实数的两倍。为此，我们希望有一种能够达到相同功能但数据量又不大的变换。 在此期望下，产生了DCT变换。 一维离散余弦变换 二维离散余弦变换 离散余弦变换的计算 DCT变换的应用： 余弦变换实际上是傅立叶变换的实数部分。 余弦变换主要用于图像的压缩，如目前的国际压缩标准的JPEG格式中就用到了DCT变换。具体的做法与DFT 相似。给高频系数大间隔量化，低频部分小间隔量化。 d=zeros(32,32); %图象大小32*32 d(13:20,13:20)=1; %中心白色方块大小为8*8 subplot(221); imshow(d,notruesize); title(原始图像) D=fft2(d); %计算图象d的傅立叶变换 subplot(222); %显示图象d的傅立叶变换谱 imshow(abs(D),[-1 5],notruesize);title(傅立叶变换谱) subplot(223); %显示图象d的傅立叶变换对数谱 imshow(log(abs(D)),[-1 5],notruesize);title(傅立叶变换对数谱) subplot(224); DF=fftshift(D); %显示图象d的傅立叶变换中心谱 imshow(log(abs(DF)),[-1 5],notruesize);title(傅立叶变换中心谱) 程序生成一个矩形函数，区域内像素值为1，区域外为0。然后对矩形做二维傅立叶变换，由于图像的傅立叶变换矩阵元素一般是复数，不能直接显示，需要调用abs函数对变换后的结果求模，图中下面两幅图分别是傅立叶变换的对数谱和中心谱。 一维离散余弦反变换 二维离散反余弦变换 如果令N=4，由一维解析式定义可得如下展开式： 写成矩阵形式： [F(u)]=[A][f(x)] 同理，可得到反变换展开形式： 写成矩阵形式： [f(x)]=[A]T[F(u)] 二维离散余弦变换为： [F(u,v)]=[A][f(x,y)][A]T [f(x,y)]=[A]T[F(u,v)] [A]T 与傅立叶变换一样，离散余弦变换可以由定义出发进行计算，但这样的计算量太大，在实际应用中很不方便，寻找快速算法 首先，从定义出发，作如下推导 取实部的意思 如果把时域数据向量作系列延拓，即 (9)平均值性质 平均值定义如下 平均值性质如下： 即： 结论：二维离散函数的平均值等于其傅立叶变换在频率原点处值的1/MN。 (10)卷积定理： f(x,y)*h(x,y) = F(u,v)H(u,v) f(x,y)h(x,y) = F(u,v)*H(u,v) 二维傅立叶变换(幅值及相位)意义 左边一列: 上方为原始图像,下方为本图的相关说明说明; 中间一列: 上图幅值谱,下图为根据幅值谱的傅立叶逆变换(忽略相位信息,设相位为0); 右边一列: 上图相位谱,下图为根据相位谱的傅立叶逆变换(忽略幅值信息,设幅值为某一常数); 图像的说明 返回 返回 返回 另一幅图像效果 压缩率为：1.7：1 压缩率为：2.24：1 压缩率为：3.3：1 返回 压缩率为：8