![数字图像处理及应用:使用MATLAB分析与实现](https://wfqqreader-1252317822.image.myqcloud.com/cover/293/25111293/b_25111293.jpg)
3.4 离散沃尔什变换
DFT或DCT变换都是以正弦或余弦三角函数为基本的正交函数基,而DFT的快速运算是复数范围内进行运算,不仅运算量庞大,而且运算复杂,DCT变换虽然避免了复数运算,但需要进行三角函数运算,运算复杂程度依然很高。因此,无论DFT,还是DCT运算,占用时间仍然较多,而在实际应用中,有时需要更为有效和便利的变换方法。沃尔什变换就是其中一种,它有效地避免了因数据本身原因所产生的运算复杂性,它由只有+1和—1两个数值所构成的完备正交基组成。由于沃尔什函数基是二值正交基,与数字逻辑的两个状态相对应,因而更加适用于计算机技术、数字信号处理等应用领域。
此外,与DFT相比,沃尔什变换减少了存储空间并提高了运算速度,这一点对图像处理来说是至关重要的,特别是在大量数据需要进行实时处理时,沃尔什变换更加显示出其优越性。
3.4.1 一维离散沃尔什变换
沃尔什(Walsh)变换与沃尔什函数密切相关,沃尔什函数是1923年由美国数学家沃尔什首先提出的。沃尔什函数系是一个完备正交函数系,其值只能取+1和—1。从排列次序上可将沃尔什函数分为三种定义方法:一是按照佩利排列来定义(按自然排序);二是按照沃尔什排列来定义(按列率排序);三是按照哈达玛排列来定义,又称为哈达玛变换。
1. 沃尔什正变换
设f(x)表示N点的一维离散序列,则一维沃尔什变换定义为
![](https://epubservercos.yuewen.com/ED2C5A/13467200603426206/epubprivate/OEBPS/Images/Figure-P88_15677.jpg?sign=1739623613-lDhix6Ac5wTaKFU8z83STpqkqXtHkR5i-0-66823457456cc990cd57245dfd3b14f0)
其中,一维沃尔什变换核为
![](https://epubservercos.yuewen.com/ED2C5A/13467200603426206/epubprivate/OEBPS/Images/Figure-P88_15678.jpg?sign=1739623613-g8DfRjfnbAQ9tZHIamIAvCoCtHLVAp0y-0-a9f6e351c0e082b5ea9d8a45d90ed1bf)
式中,u=0,1,2,3,…,N—1;x=0,1,2,3,…,N—1。N是沃尔什变换的阶数,N=2n。bi(z)是z的二进制数的第i位数值,取值为0或1。如i=6,由于6的二进制表示为110,因此b0(z)=0,b1(z)=1,b2(z)=1。
2. 沃尔什逆变换
一维离散沃尔什逆变换定义为
![](https://epubservercos.yuewen.com/ED2C5A/13467200603426206/epubprivate/OEBPS/Images/Figure-P88_15682.jpg?sign=1739623613-kXAZnRwlDx3HaPnjnId8YUYXLW2uotGH-0-13c48d2f58f8705a812309ec28bd6a75)
逆变换的核为
h(x,u)=g(x,u)
一维沃尔什正、反变换核相同,沃尔什变换核是一个对称阵列,其行和列是正交的。沃尔什正、反变换形式本质上相同,因此,计算沃尔什变换的算法可直接用来求其逆变换。一维沃尔什变换也具有快速算法,简称为FWT,在形式上和FFT算法类似。当N=8时,其变换核用矩阵表示为
![](https://epubservercos.yuewen.com/ED2C5A/13467200603426206/epubprivate/OEBPS/Images/Figure-P89_15684.jpg?sign=1739623613-FjhIzad3EXNUqpEMds1uDGpHIQzADUmd-0-b827be0252d9739ad2bf7b48d05aef31)
3.4.2 二维离散沃尔什变换
1. 二维沃尔什正变换
设f(x,y)表示M×N的二维离散序列,则二维沃尔什变换定义为
![](https://epubservercos.yuewen.com/ED2C5A/13467200603426206/epubprivate/OEBPS/Images/Figure-P89_15686.jpg?sign=1739623613-8hdkXZn0uOVs7eAwCF3IKnKSi6id7o5V-0-d2e4b9aeef7bf0d12c821326f0a94db1)
其中,二维离散沃尔什变换核为
![](https://epubservercos.yuewen.com/ED2C5A/13467200603426206/epubprivate/OEBPS/Images/Figure-P89_15688.jpg?sign=1739623613-VypTq9Jp7Bl6sU6xhPTTb81HWC1fFTjq-0-c114c646c8389da034ddd71fddff4a2a)
式中,M=2m,N=2n。
二维离散沃尔什变换的变换核是可分离的,即
![](https://epubservercos.yuewen.com/ED2C5A/13467200603426206/epubprivate/OEBPS/Images/Figure-P89_15690.jpg?sign=1739623613-G2pShg7vnwmBtZsjFitBBK6OIQNDzixd-0-48857b1613fea50337fac540a3c4dfb9)
根据沃尔什变换的定义形式可以得出,二维沃尔什变换具有可分离特性,即一次二维沃尔什变换可以通过两次一维沃尔什变换来实现。
2. 二维沃尔什逆变换
二维沃尔什逆变换定义为
![](https://epubservercos.yuewen.com/ED2C5A/13467200603426206/epubprivate/OEBPS/Images/Figure-P89_15692.jpg?sign=1739623613-nmsqievjdNWYETzhueqXJGBS1PvJwnhY-0-39c6373fef8a34234e9aec0e8ce87e50)
二维离散沃尔什逆变换核为
![](https://epubservercos.yuewen.com/ED2C5A/13467200603426206/epubprivate/OEBPS/Images/Figure-P89_15694.jpg?sign=1739623613-3TCEqVNaT5080KoslVzOQSaIkF4tMRoh-0-42d118189cdb2b8ef98af2c5f191081d)
式中,M=2m,N=2n。
二维沃尔什逆变换的核为
h(x,u,y,v)=g(x,u,y,v)
同样,二维逆变换具有可分离性。二维沃尔什变换也可以表示为矩阵形式,即
![](https://epubservercos.yuewen.com/ED2C5A/13467200603426206/epubprivate/OEBPS/Images/Figure-P90_15701.jpg?sign=1739623613-7DK8II2ssXrMe0vKCwYdF1iUYXcVbmcO-0-38eb31f9be5febd4487a7cc71dfb6259)
式中,G1为M×M变换核方阵;G2为N×N变换核方阵。
3.4.3 沃尔什变换的频谱
和DCT变换相同,二维沃尔什变换也具有能量集中的性质,原始图像数据越是均匀分布,沃尔什变换后的数据越集中于矩阵的边角上,因此,二维沃尔什变换也常用于压缩图像信息。
例3-5:沃尔什变换实例
如图3-8所示是沃尔什变换的实例,图3-8(a)为对应于图3-7(a)的沃尔什变换结果,图3-8(b)为对应于图3-7(c)的沃尔什变换结果。
![](https://epubservercos.yuewen.com/ED2C5A/13467200603426206/epubprivate/OEBPS/Images/Figure-P90_4883.jpg?sign=1739623613-CQLAl5P9fsXEMBmIURHfrqbtXQyQPEIq-0-bc296fd6a3b3a14a7bc98cc2e771acc3)
图3-8 沃尔什变换的频谱分布
沃尔什变换是将一个函数变换成取值为+l或—1的基本函数构成的级数,用它来逼近数字脉冲信号时要比DFT有利。因此,它在图像传输、通信技术和数据压缩中获得了广泛的使用。同时,沃尔什变换是实数,所以对工程应用问题,沃尔什变换的存储量比DFT少,而且运算速度非常快。
例3-6:有两个二维数字图像信号,即
(1)
(2)
分别求f1和f2的二维沃尔什变换。
离散二维沃尔什变换的求解过程如下:
根据理论,对于(1)和(2),均有M=N=4,因此,其二维沃尔什变换的核为
![](https://epubservercos.yuewen.com/ED2C5A/13467200603426206/epubprivate/OEBPS/Images/Figure-P91_15707.jpg?sign=1739623613-TJS71hfcEZy4C2pWsgtPFqGzQYwwjHA8-0-47b0e5c09e3f8f2830086b0e7cae86ae)
由此可得
![](https://epubservercos.yuewen.com/ED2C5A/13467200603426206/epubprivate/OEBPS/Images/Figure-P91_15708.jpg?sign=1739623613-sggbeCrlv1R3vK1AbDJQ2xrXM35y4Ut9-0-2e670ca8101ed9c5e9fadb287a321a50)