《课062.12.2二维卷积与DFT》由会员分享，可在线阅读，更多相关《课062.12.2二维卷积与DFT（36页珍藏版）》请在金锄头文库上搜索。
1、第二章第二章 数字图象的线性处理数字图象的线性处理2. .1 离散卷积与离散相关离散卷积与离散相关2. .2 二维离散傅里叶变换（二维离散傅里叶变换（2D-DFT）2. .3 离散沃尔什变换（离散沃尔什变换（DWT）2.1.1 二维离散卷积二维离散卷积 2.1.2 二维离散卷积定理二维离散卷积定理 2.1.3 二维离散相关二维离散相关 2.2.1 定义与讨论定义与讨论 2.2.2 矢量矩阵表示矢量矩阵表示 2.2.3 傅里叶变换是酉变换傅里叶变换是酉变换 2.2.4 2.3.1 离散沃尔什函数离散沃尔什函数 2.3.2 一维一维DWT2.3.3 二维二维DWT8/26/20241第二章第二章
2、数字图象的线性处理数字图象的线性处理 2.6.1 奇数点的奇数点的CosineCosine变换变换2. .4 离散哈达玛变换（离散哈达玛变换（DHT）2. .5 离散卡离散卡- -洛变换（洛变换（DKLT） 2.4.1 哈达玛矩阵哈达玛矩阵 2.4.2 一维一维DHT 2.4.3 二维二维DHT 2.5.1 一维连续一维连续K-L展开展开 2.5.2 一维离散一维离散KLT 2.5.3 一维离散一维离散KLT T2. .6 离散余弦变换（离散余弦变换（DCT） 2.6.1 偶数点的偶数点的CosineCosine变换变换 2.6.3 DCT的性能的性能8/26/202422. .1 离散卷积与
3、离散相关离散卷积与离散相关 2.1.1 二维离散卷积二维离散卷积 2.1.2 离散卷积定理离散卷积定理v 定义定义v 对卷积矩阵对卷积矩阵Te的讨论的讨论v 小结小结 2.1.3 二维离散相关二维离散相关v 定义定义v 二维离散相关定理与性质二维离散相关定理与性质返回返回8/26/20243v 定义定义2. .1. .1 二维二维离散卷积离散卷积定义二维离散卷积定义二维离散卷积：( ( 个样本值个样本值) )( ( 个样本值个样本值) )设两个二维离散函数设两个二维离散函数：式中式中， 与与 分别是分别是 、 的周期化函数的周期化函数。上页上页8/26/202442. .1. .1 二维二维离
4、散卷积离散卷积上页上页即即： 和和 的周期为的周期为：定义所给出的定义所给出的 阶函数阵列阶函数阵列，是二维离散卷积的一个周期是二维离散卷积的一个周期。8/26/202452. .1. .1 二维二维离散卷积离散卷积 例题例题 求两个求两个22 阶二维离散函数的卷积阶二维离散函数的卷积： 解法一（解析法）解法一（解析法）解法一（解析法）解法一（解析法） (2) 求求 F和和G 的列矢量（按行扫描的列矢量（按行扫描- -堆叠方式）堆叠方式）：上页上页 (1) F和和G周期化周期化： 周期周期 M = N =A+B-1=38/26/20246卷积矩阵卷积矩阵 2. .1. .1 二维二维离散卷积离
5、散卷积(3) 按一维离散卷积方法计算卷积按一维离散卷积方法计算卷积： 上页上页8/26/202472. .1. .1 二维二维离散卷积离散卷积(4) 卷积结果卷积结果( (矩阵形式矩阵形式) )：l Te为为 阶方阵阶方阵；l 有有 个分块子阵每个子阵为个分块子阵每个子阵为 阶阶；上页上页l 共有共有 M 组相同的子阵组相同的子阵。 其中其中卷积矩阵卷积矩阵：8/26/202482. .1. .1 二维二维离散卷积离散卷积其中，分块子阵其中，分块子阵：10 02 22 210 00 012 2v 对卷积矩阵的进一步讨论对卷积矩阵的进一步讨论上页上页 箭头所指各列与箭头所指各列与 f fe e
6、的零元素相乘，因此，可改的零元素相乘，因此，可改动相关元素值，使每个分块子阵均成为循环矩阵。动相关元素值，使每个分块子阵均成为循环矩阵。T Te e1 1T Te e1 1T Te e1 1T Te e2 2T Te e2 2T Te e2 2T Te e3 3T Te e3 3T Te e3 38/26/202492. .1. .1 二维二维离散卷积离散卷积小结：小结：Te 为分块循环矩阵，其中各列（行）是分块子阵的循环移位为分块循环矩阵，其中各列（行）是分块子阵的循环移位； 各分块子阵是由各分块子阵是由 g e子矢量诸元素构成的循环矩阵。子矢量诸元素构成的循环矩阵。上页上页本例本例：其中各
7、子矢量为其中各子矢量为：8/26/2024102. .1. .1 二维二维离散卷积离散卷积返回返回1-22-12-11-2-21-121-22-1xy0(2) 垂直垂直/ /水平折叠水平折叠1234(1) 将将F和和G表示成表示成“掩模掩模”形式形式 解法二（图解法）解法二（图解法）解法二（图解法）解法二（图解法）8/26/2024112. .1. .1 二维二维离散卷积离散卷积上页上页123-2142-12342-211-1 (3) 水平水平/ /垂直平移折叠后的掩模，求乘积和垂直平移折叠后的掩模，求乘积和24-21-1231342-212-11422-23-11142-2-121318/2
8、6/2024122. .1. .2 离散卷积定理离散卷积定理上页上页 二维二维 若二维离散函数的傅里叶变换对为若二维离散函数的傅里叶变换对为则有则有则有则有 一维一维 若离散函数若离散函数 、 的傅里叶变换分别为的傅里叶变换分别为 、 ，即即8/26/2024132. .1. .3 二维二维离散相关离散相关上页上页定义二维离散相关定义二维离散相关：以及以及v 定义定义 对两个二维离散实函数对两个二维离散实函数： 8/26/2024142. .1. .3 二维二维离散相关离散相关上页上页则则v 二维相关定理二维相关定理若若自相关函数自相关函数：v 自相关函数及其性质自相关函数及其性质8/26/2
9、024152. .1. .3 二维二维离散相关离散相关上页上页当当 时时，自相关函数的傅氏变换自相关函数的傅氏变换：8/26/2024162. .2 二维离散傅里叶变换二维离散傅里叶变换 2.2.1 定义与讨论定义与讨论 2.2.2 矩阵矢量表示矩阵矢量表示 2.2.3 傅里叶变换是酉变换傅里叶变换是酉变换下页下页 2.2.4 常用性质常用性质8/26/2024172. .2. .1 定义与讨论定义与讨论 设二维数据阵列设二维数据阵列 f (x,y)为为 NN 方阵方阵, , 定义定义2D- -DFT： 反变换：反变换：F(u, v)也是也是 NN 方阵方阵。 正变换：正变换：v 2 2D-D
10、FTD-DFT定义定义上页上页8/26/2024182. .2. .1 定义与讨论定义与讨论上页上页 列N -1012LN -112Muv行 列N -1012LN -112Mxy行变换对变换对：v 讨论讨论 变量说明变量说明：x, y 离散图像像素点在空间的行、列位置标号；离散图像像素点在空间的行、列位置标号；u, v 表示变换域中样点行、列位置标号表示变换域中样点行、列位置标号。 8/26/2024192. .2. .1 定义与讨论定义与讨论上页上页记记 其中，正变换核其中，正变换核： 变换核变换核： 可分离性可分离性：行分量行分量列分量列分量 令令，则，则DFT的变换核为的变换核为；反变换
11、核反变换核： 其中其中：8/26/2024202. .2. .1 定义与讨论定义与讨论改写正变换公式改写正变换公式：两次一维变换两次一维变换先作行运算 xu 一个二维一个二维DFT，可连续运用两次一维，可连续运用两次一维DFT来实现，因而，来实现，因而，可采用一维可采用一维FFT进行快速运算进行快速运算。示意如下：示意如下： 结论结论 再作列运算 yv返回目录返回目录8/26/2024212. .2. .2 矢量矢量- -矩阵表示矩阵表示则则 DFT的矢量形式为的矢量形式为：v DFT的矢量表示的矢量表示 正变换正变换A变换矩阵变换矩阵其中其中：频率域矩阵为频率域矩阵为 Q ( (列列) )矢
12、量为矢量为 q设图像空间域矩阵为设图像空间域矩阵为 F ( (列列) )矢量为矢量为 f对应对应 反变换反变换B反变换矩阵反变换矩阵其中其中：上页上页8/26/202422 举例举例 设空域矩阵设空域矩阵F与频域矩阵与频域矩阵Q分别为分别为：2. .2. .2 矢量矢量- -矩阵表示矩阵表示相应的列矢量为相应的列矢量为：v 对变换矩阵对变换矩阵A的讨论的讨论 （举例举例 推广）推广），上页上页问题问题：A=? 为此，可根据定义式分别求相应的变换系数为此，可根据定义式分别求相应的变换系数。由由DFT变换式变换式8/26/2024232. .2. .2 矢量矢量- -矩阵表示矩阵表示由正变换定义式
13、可得由正变换定义式可得：式中式中：上页上页8/26/202424变换矩阵变换矩阵变换矩阵变换矩阵A A2. .2. .2 矢量矢量- -矩阵表示矩阵表示上页上页将以上方程组写成矩阵形式将以上方程组写成矩阵形式： 进一步进一步改写变换矩阵改写变换矩阵：8/26/202425即即即即A A的的的的行、列是可分离的行、列是可分离的行、列是可分离的行、列是可分离的2. .2. .2 矢量矢量- -矩阵表示矩阵表示上页上页令令 ， 推广推广 一般地，变换矩阵一般地，变换矩阵A的行、列分量为的行、列分量为显然存在关系显然存在关系： 直积直积：则则A可表为可表为AC与与AR的的8/26/2024262. .
14、2. .2 矢量矢量- -矩阵表示矩阵表示l 正变换正变换：l 反变换反变换：等式两边左乘等式两边左乘 ，右乘，右乘 v DFT的矩阵表示的矩阵表示2 2阶阶DFT的矩阵形式可表示为的矩阵形式可表示为： 得到得到：返回目录返回目录8/26/202427 对于变换矩阵对于变换矩阵T，酉变换酉变换满足以下条件：满足以下条件：T-1 = ( T* ) T即即T的逆为的逆为T的共轭转置，这时称的共轭转置，这时称T为为“酉矩阵酉矩阵”。当。当T为实酉矩阵时，为实酉矩阵时，则为正交矩阵。正交矩阵中各行（列）矢量是正交归一的，即则为正交矩阵。正交矩阵中各行（列）矢量是正交归一的，即 T-1 T = ( T*
15、 ) T T = I TT-1 = T ( T* ) T = I2. .2. .3 傅里叶变换是傅里叶变换是酉酉变换变换上页上页酉变换是普遍意义上的线性正交变换，而且是可逆的。酉变换是普遍意义上的线性正交变换，而且是可逆的。v 何谓酉变换？何谓酉变换？对于对于DFT，若能证明变换矩阵，若能证明变换矩阵A满足满足 则则A是酉矩阵是酉矩阵，进而，进而DFT是酉变换。是酉变换。8/26/2024282. .2. .3 傅里叶变换是酉变换傅里叶变换是酉变换上页上页其逆矩阵为其逆矩阵为 已知已知DFT变换矩阵变换矩阵A的行、列分量为的行、列分量为比较显见比较显见：8/26/2024292. .2. .3
16、 傅里叶变换是酉变换傅里叶变换是酉变换 由于由于 ,所以所以 即即DFT正变换矩阵正变换矩阵A是酉矩阵是酉矩阵。 结论结论 傅里叶变换是酉变换（正交变换）傅里叶变换是酉变换（正交变换）。因而逆变换矩阵因而逆变换矩阵B也是酉矩阵也是酉矩阵。根据根据 有有故逆变换核同样是可分离的。故逆变换核同样是可分离的。且因且因进一步进一步返回目录返回目录8/26/2024302. .2. .4 傅里叶变换的常用性质傅里叶变换的常用性质上页上页可得可得 因此，因此，F (u,v)和和 f (x, y)均具有周期性，周期为均具有周期性，周期为NN。根据根据DFT定义和复指数函数的周期性定义和复指数函数的周期性：同
17、理，反变换同理，反变换v 周期性周期性8/26/2024312. .2. .4 傅里叶变换的常用性质傅里叶变换的常用性质上页上页与一维与一维与一维与一维DFTDFT类似，在求类似，在求类似，在求类似，在求 f f ( (x x, , y y) )二维二维二维二维DFTDFT时，意味着在空间域和频时，意味着在空间域和频时，意味着在空间域和频时，意味着在空间域和频率域两方面都周期化了率域两方面都周期化了率域两方面都周期化了率域两方面都周期化了。f f ( (x x, , y y) )和和和和F F ( (u, vu, v) )的的的的 N NN N 个样本，均表个样本，均表个样本，均表个样本，均表
18、示二维周期函数的一个周期。因此，欲了解示二维周期函数的一个周期。因此，欲了解示二维周期函数的一个周期。因此，欲了解示二维周期函数的一个周期。因此，欲了解F F ( (u, vu, v) )的全貌，只需的全貌，只需的全貌，只需的全貌，只需一个完整的周期即可。一个完整的周期即可。一个完整的周期即可。一个完整的周期即可。二维二维二维二维DFTDFT也具有也具有也具有也具有“循环循环循环循环”特性，如循环空特性，如循环空特性，如循环空特性，如循环空间位移，循环卷积等。间位移，循环卷积等。间位移，循环卷积等。间位移，循环卷积等。 结论结论 8/26/2024322. .2. .4 傅里叶变换的常用性质傅
19、里叶变换的常用性质上页上页 由周期性由周期性：作变量代换作变量代换 u u， v v，并取，并取共轭共轭： 可得可得： 等式左边等式左边： 等式右边（观察指数部分）等式右边（观察指数部分）:v 共轭对称性共轭对称性8/26/2024332. .2. .4 傅里叶变换的常用性质傅里叶变换的常用性质上页上页 令令 m , n = 1，进一步有，进一步有 m , n = 0， 解释共轭对称性解释共轭对称性 假设假设 f (x, y) 为为正正、实实函数函数(1)(1) F (u, v) 实部实部8/26/2024342. .2. .4 傅里叶变换的常用性质傅里叶变换的常用性质(2)(2) F (u, v) 虚部虚部返回目录返回目录8/26/202435矩阵直积矩阵直积返回返回 矩阵的直积运算，又称为矩阵的直积运算，又称为“克罗内克克罗内克”积，符号为积，符号为，规则如下规则如下：若若 则则 8/26/202436
