数字图像处理及应用:使用MATLAB分析与实现
上QQ阅读APP看书,第一时间看更新

3.4 离散沃尔什变换

DFT或DCT变换都是以正弦或余弦三角函数为基本的正交函数基,而DFT的快速运算是复数范围内进行运算,不仅运算量庞大,而且运算复杂,DCT变换虽然避免了复数运算,但需要进行三角函数运算,运算复杂程度依然很高。因此,无论DFT,还是DCT运算,占用时间仍然较多,而在实际应用中,有时需要更为有效和便利的变换方法。沃尔什变换就是其中一种,它有效地避免了因数据本身原因所产生的运算复杂性,它由只有+1和—1两个数值所构成的完备正交基组成。由于沃尔什函数基是二值正交基,与数字逻辑的两个状态相对应,因而更加适用于计算机技术、数字信号处理等应用领域。

此外,与DFT相比,沃尔什变换减少了存储空间并提高了运算速度,这一点对图像处理来说是至关重要的,特别是在大量数据需要进行实时处理时,沃尔什变换更加显示出其优越性。

3.4.1 一维离散沃尔什变换

沃尔什(Walsh)变换与沃尔什函数密切相关,沃尔什函数是1923年由美国数学家沃尔什首先提出的。沃尔什函数系是一个完备正交函数系,其值只能取+1和—1。从排列次序上可将沃尔什函数分为三种定义方法:一是按照佩利排列来定义(按自然排序);二是按照沃尔什排列来定义(按列率排序);三是按照哈达玛排列来定义,又称为哈达玛变换。

1. 沃尔什正变换

fx)表示N点的一维离散序列,则一维沃尔什变换定义为

其中,一维沃尔什变换核为

式中,u=0,1,2,3,…,N—1;x=0,1,2,3,…,N—1。N是沃尔什变换的阶数,N=2nbiz)是z的二进制数的第i位数值,取值为0或1。如i=6,由于6的二进制表示为110,因此b0z)=0,b1z)=1,b2z)=1。

2. 沃尔什逆变换

一维离散沃尔什逆变换定义为

逆变换的核为

hx,u=gx,u

一维沃尔什正、反变换核相同,沃尔什变换核是一个对称阵列,其行和列是正交的。沃尔什正、反变换形式本质上相同,因此,计算沃尔什变换的算法可直接用来求其逆变换。一维沃尔什变换也具有快速算法,简称为FWT,在形式上和FFT算法类似。当N=8时,其变换核用矩阵表示为

3.4.2 二维离散沃尔什变换

1. 二维沃尔什正变换

fx,y)表示M×N的二维离散序列,则二维沃尔什变换定义为

其中,二维离散沃尔什变换核为

式中,M=2m,N=2n

二维离散沃尔什变换的变换核是可分离的,即

根据沃尔什变换的定义形式可以得出,二维沃尔什变换具有可分离特性,即一次二维沃尔什变换可以通过两次一维沃尔什变换来实现。

2. 二维沃尔什逆变换

二维沃尔什逆变换定义为

二维离散沃尔什逆变换核为

式中,M=2m,N=2n

二维沃尔什逆变换的核为

hx,u,y,v=gx,u,y,v

同样,二维逆变换具有可分离性。二维沃尔什变换也可以表示为矩阵形式,即

式中,G1M×M变换核方阵;G2N×N变换核方阵。

3.4.3 沃尔什变换的频谱

和DCT变换相同,二维沃尔什变换也具有能量集中的性质,原始图像数据越是均匀分布,沃尔什变换后的数据越集中于矩阵的边角上,因此,二维沃尔什变换也常用于压缩图像信息。

例3-5:沃尔什变换实例

如图3-8所示是沃尔什变换的实例,图3-8(a)为对应于图3-7(a)的沃尔什变换结果,图3-8(b)为对应于图3-7(c)的沃尔什变换结果。

图3-8 沃尔什变换的频谱分布

沃尔什变换是将一个函数变换成取值为+l或—1的基本函数构成的级数,用它来逼近数字脉冲信号时要比DFT有利。因此,它在图像传输、通信技术和数据压缩中获得了广泛的使用。同时,沃尔什变换是实数,所以对工程应用问题,沃尔什变换的存储量比DFT少,而且运算速度非常快。

例3-6:有两个二维数字图像信号,即

(1)

(2)

分别求f1f2的二维沃尔什变换。

离散二维沃尔什变换的求解过程如下:

根据理论,对于(1)和(2),均有M=N=4,因此,其二维沃尔什变换的核为

由此可得