参考图5-1的信号移动规律,二维全相位信号处理可以按照行或者列优先的方式来实现,下面以行优先为例给出实现原理图。式表明,二维全相位信号处理对应xi,j的输出等于在行列方向上以xi,j为中心的共2个数据以A为加权的总和。图5-8 二维全相位行信号处理框图图5-8中,输入是按照行方向移动的大小为N×N的数据块,假设待处理图像大小为M×L,且进行边界扩展成×以保证输出图像大小不变。......
2023-06-23
对于二维余弦变换,α=C,β=CT,代入到式(5-64)并简单整理可得:
相乘的四个余弦基元素不能进一步合并,所以余弦全相位核没有更具体的表达形式。对于二维沃尔什变换,α=WN,β=α,代入到式(5-64)并简单整理可得:
式(5-104)中符号表示异或运算。对于ap2 DCT的全相位变换核,有
当N是2的整数幂时,pu=(N-1-p)(N-1-u),所以式(5-105)的结论对二维apDWT的核仍然成立。基于此结论的ap2 DCT和ap2 DWT的全相位变换基也具有对称性。根据式(5-94)的定义,位于第四象限的A—元素Ar,c的中心对称是位于第一象限的A++的元素A-r,-c,两者的关系(假定双基窗均经过归一化处理)推导如下:
式(5-106)成立的条件只是要求基窗具有中心对称性,即对于任意传输矩阵H,全相位处理的特性都具有严格零相位。同样可证明,当前后基窗分别具有水平和垂直对称性时,全相位基又满足如下:
对于第四象限的A—元素Ar,c的转置,可推导出如下结论:
当前后基窗满足主对称性,且全相位核满足如下条件时,Ar,c也具有了主对称性。
结合式(5-103),若式(5-109)成立,则有
对于二维apDWT,则有
综合以上,对于ap2d DCT和ap2d DWT,我们有如下结论:
1)当前基窗F和后基窗B具有中心对称性时,全相位基中心对称;
2)当前后基窗水平/垂直对称时,全相位基水平/垂直对称;
3)当前后基窗和传输矩阵H主对称时,全相位基满足主对称性;
在对二维信号进行全相位处理时,基窗的选择和H的设定是必须考虑的问题。如果认为其他像素与中心像素的相关性与距离成反比且各向同性,则基窗应定义为中心各向对称,即同时满足水平、垂直、中心和主对称。根据信号相关性分布方向来自适应设置基窗以达到一定条件下的最优全相位信号处理效果是值得研究的内容。对于H的定义,则与处理的目标任务相关,即由在各个方向上不同列率信号分量的取舍来决定。如果行列上情形相同,则H具有主对称性。下面的讨论中,除非特别指出,前后基窗相等且具有中心各向对称性,H也满足主对称性。假设H定义如下:
以图5-20给出的3种卷基窗分别形成的二维apDCT传输特性如图5-21所示。
图5-21 ap2nsd DCT特性曲面及曲线图
从图中可看出apDCT的平滑性好于apDFT,ap2d DCT的衰减特性比ap2n DCT低20dB,但由于采用与ap2n DCT和ap2s DCT相同的归一化条件,峰值幅度低于1且过渡带宽度增大。相比而言,ap2s DCT在衰减特性和过渡带宽度都表现出良好特性。生成图5-21所述功能的MATLAB代码(只列出与ap2nsd DFT不同的部分)如下:
与上相同条件形成的二维ap2nsd DWT传输特性如图5-22所示。
图5-22 ap2nsd DWT特性曲面及曲线图
生成图5-22结果的MATLAB代码与ap2nsd DCT只有函数getA不同,具体如下:
上面给出了DCT和DWT变换情况下设计二维全相信号处理基本方法,虽然特性曲面物理意义不如DFT直接,即通过H来定义的传输特性不能与曲面特征建立直接对应关系,但提供了信号处理新思路,且具有严格零相位和优良平滑性等特点。
有关全相位数字信号处理方法及MATLAB实现的文章
参考图5-1的信号移动规律,二维全相位信号处理可以按照行或者列优先的方式来实现,下面以行优先为例给出实现原理图。式表明,二维全相位信号处理对应xi,j的输出等于在行列方向上以xi,j为中心的共2个数据以A为加权的总和。图5-8 二维全相位行信号处理框图图5-8中,输入是按照行方向移动的大小为N×N的数据块,假设待处理图像大小为M×L,且进行边界扩展成×以保证输出图像大小不变。......
2023-06-23
系统原理同图2-26所示,对应的输出为式中式中Tk的表达式与apDCT时相同。此外,apDWT具有与apDCT相同的传输特性函数形式,因此也具有严格零相位的特性。图4-11 基于apDWT子带分解信号图实际上,不论是何种正交变换,按照全相位信号处理方法得到的ap都可以对信号进行去除分段效应的滤波和处理,它是建立在新概念上的一种滤波系统。......
2023-06-23
国际标准建议CCITT H.261、JPEG、MPEG中都采用8×8的块做DCT,但基于DCT的数字滤波尚不普遍。为使DCT所具有的优良特性能在数字滤波中得到充分发挥,结合全相滤波的概念提出了DCT域ap系统设计。TN-2TN-1],其传输特性为显然,apDCT具有严格零相位。实现无窗、单窗和双窗apDCT与传统滤波器特性的MATLAB代码如下:显然,apDCT具有严格零相位。利用apDCT的各子带滤波器对图像的一行进行子带滤波的结果如图4-8所示。......
2023-06-23
二维全相位在各种正交变换的基础上进行了全相位信号处理,一方面扩展了输入信号信息,另一方面对于不同相位信号进行不同权值的加权,从而使二维全相位信号处理在保留传统分块固有性质的同时也具备了一些特殊性质。可把基窗分解成若干个子基窗,从而得到信号不同加权情况下的全相位处理输出,其特性取决于子基窗的形状。以式中低通特性HLL为例,为保持全相位处理的严格零相位特性,仍使子基窗满足中心对称。......
2023-06-23
在一维信号处理中,用来抑制噪声干扰并保留原有有用信号的滤波方法以及获得各种信号变换的数字滤波器,有着极其重要和广泛的应用。所以,统计结果表明二维图像信号谱能量90%以上集中在二维低频区域,如果进行2∶1或更高亚取样,谱也未发生重叠。为了克服边缘效应,一种有效的方法是对图像进行重叠滤波,在重叠滤波中的子图像在行列方向上逐点进行重叠。......
2023-06-23
由第5章内容可知,二维DCT运算可以转化成DFT运算,然后取结果的实部作为DCT变换的结果。DCT基可以设计任何大小的内插模板,而不限于2的整数次幂。以上面傅里叶基构造的结果直接得到DCT基滤波器H和对应的非重叠模板。图8-1c对应的定制H1/3如下,其反变换即为非重叠模板T1/3。......
2023-06-23
长度为N的单边带H0按上面3种方法设计生成的ap半带滤波器H的长度为2N-1。在满足条件1、2、3下,ap1 nsd半带滤波器H、单边带H0和卷积窗C的零点分布如图7-22所示。ap1ns半带滤波器的H零点由单位圆上和内及对应的倒数组成,而ap1 d半带滤波器的H零点只由单位圆内及对应的倒数组成。即把半带滤波器H进行分解后找到所有圆内的零点就可以确定分析低通滤波器h0,再按照PR条件即式关系依次求出分析高通滤波器h1,综合滤波器f1和f2。......
2023-06-23
图4-12 变换域信号处理图对于N阶系统,时刻n输入xn=[x,x(n-1),…其中X1是1维N阶全相位输入信号,即式就是全相位数字信号处理得到的输出与输入的数学关系表达。把式称为全相位变换,矩阵A称为全相位变换基,矩阵T则为全相位变换核。综合以上,变换域全相位数字信号处理一种可能实现如图4-16所示。......
2023-06-23
相关推荐