首页 理论教育二维DCT和DWT全相位技术详解

二维DCT和DWT全相位技术详解

【摘要】:对于ap2 DCT的全相位变换核,有式中符号表示异或运算。基于此结论的ap2 DCT和ap2 DWT的全相位变换基也具有对称性。以图5-20给出的3种卷基窗分别形成的二维apDCT传输特性如图5-21所示。图5-21 ap2nsd DCT特性曲面及曲线图从图中可看出apDCT的平滑性好于apDFT,ap2d DCT的衰减特性比ap2n DCT低20dB,但由于采用与ap2n DCT和ap2s DCT相同的归一化条件,峰值幅度低于1且过渡带宽度增大。

对于二维余弦变换,α=C,β=CT,代入到式(5-64)并简单整理可得:

978-7-111-48233-8-Chapter05-129.jpg

相乘的四个余弦基元素不能进一步合并,所以余弦全相位核没有更具体的表达形式。对于二维沃尔什变换,α=WN,β=α,代入到式(5-64)并简单整理可得:

978-7-111-48233-8-Chapter05-130.jpg

式(5-104)中符号978-7-111-48233-8-Chapter05-131.jpg表示异或运算。对于ap2 DCT的全相位变换核,有

978-7-111-48233-8-Chapter05-132.jpg

N是2的整数幂时,p978-7-111-48233-8-Chapter05-133.jpgu=(N-1-p978-7-111-48233-8-Chapter05-134.jpgN-1-u),所以式(5-105)的结论对二维apDWT的核仍然成立。基于此结论的ap2 DCT和ap2 DWT的全相位变换基也具有对称性。根据式(5-94)的定义,位于第四象限的A—元素Arc中心对称是位于第一象限的A++的元素A-r,-c,两者的关系(假定双基窗均经过归一化处理)推导如下:

978-7-111-48233-8-Chapter05-135.jpg

式(5-106)成立的条件只是要求基窗具有中心对称性,即对于任意传输矩阵H,全相位处理的特性都具有严格零相位。同样可证明,当前后基窗分别具有水平和垂直对称性时,全相位基又满足如下:

978-7-111-48233-8-Chapter05-136.jpg

对于第四象限的A—元素Arc的转置,可推导出如下结论:

978-7-111-48233-8-Chapter05-137.jpg

当前后基窗满足主对称性,且全相位核满足如下条件时,Arc也具有了主对称性。

978-7-111-48233-8-Chapter05-138.jpg

结合式(5-103),若式(5-109)成立,则有

978-7-111-48233-8-Chapter05-139.jpg

对于二维apDWT,则有

978-7-111-48233-8-Chapter05-140.jpg

综合以上,对于ap2d DCT和ap2d DWT,我们有如下结论:

1)当前基窗F和后基窗B具有中心对称性时,全相位基中心对称;

2)当前后基窗水平/垂直对称时,全相位基水平/垂直对称;

3)当前后基窗和传输矩阵H主对称时,全相位基满足主对称性;

在对二维信号进行全相位处理时,基窗的选择和H的设定是必须考虑的问题。如果认为其他像素与中心像素的相关性与距离成反比且各向同性,则基窗应定义为中心各向对称,即同时满足水平、垂直、中心和主对称。根据信号相关性分布方向来自适应设置基窗以达到一定条件下的最优全相位信号处理效果是值得研究的内容。对于H的定义,则与处理的目标任务相关,即由在各个方向上不同列率信号分量的取舍来决定。如果行列上情形相同,则H具有主对称性。下面的讨论中,除非特别指出,前后基窗相等且具有中心各向对称性,H也满足主对称性。假设H定义如下:

978-7-111-48233-8-Chapter05-141.jpg

以图5-20给出的3种卷基窗分别形成的二维apDCT传输特性如图5-21所示。

978-7-111-48233-8-Chapter05-142.jpg

图5-21 ap2nsd DCT特性曲面及曲线图

从图中可看出apDCT的平滑性好于apDFT,ap2d DCT的衰减特性比ap2n DCT低20dB,但由于采用与ap2n DCT和ap2s DCT相同的归一化条件,峰值幅度低于1且过渡带宽度增大。相比而言,ap2s DCT在衰减特性和过渡带宽度都表现出良好特性。生成图5-21所述功能的MATLAB代码(只列出与ap2nsd DFT不同的部分)如下:

978-7-111-48233-8-Chapter05-143.jpg

978-7-111-48233-8-Chapter05-144.jpg

978-7-111-48233-8-Chapter05-145.jpg

与上相同条件形成的二维ap2nsd DWT传输特性如图5-22所示。

978-7-111-48233-8-Chapter05-146.jpg

图5-22 ap2nsd DWT特性曲面及曲线图

生成图5-22结果的MATLAB代码与ap2nsd DCT只有函数getA不同,具体如下:

978-7-111-48233-8-Chapter05-147.jpg

978-7-111-48233-8-Chapter05-148.jpg

上面给出了DCT和DWT变换情况下设计二维全相信号处理基本方法,虽然特性曲面物理意义不如DFT直接,即通过H来定义的传输特性不能与曲面特征建立直接对应关系,但提供了信号处理新思路,且具有严格零相位和优良平滑性等特点。