专利名称：对称延拓双正交小波变换矩阵的构造方法
技术领域：
本发明属于计算机图像处理技术领域。
背景技术：
小波变换，无论分解还是重构都需要逐级计算，而实际所处理的信号和图像都是有限的，现存文献很少研究有限长度信号小波变换的简易算法。如果小波变换的每一级计算都能通过矩阵-向量乘积实现，会非常方便。只要是基于小波变换的方法，如图像压缩、图像去噪、图像增强等，都可以用小波变换矩阵实现。而当滤波器长度大于信号长度时，信号需要被延拓处理。信号延拓方式为周期延拓和对称延拓，发明者曾给出信号周期延拓方式下小波变换矩阵，用于FRIT图像去噪，虽然简化了计算，却表现有明显的边界效应。由于自然图像块间的连续性，信号或图像的对称延拓更有实际意义，所以若能构造出信号对称延拓方式下小波变换矩阵，将更具现实意义和实用价值。
发明内容本发明的目的是克服现有技术存在的不足，提供一种对称延拓双正交小波变换矩阵的构造方法；使信号在对称延拓方式下，有限长信号不需逐级计算而直接得到小波系数的分解矩阵和由这些小波系数重构原信号的重构矩阵的构造方法。
本发明提供的对称延拓双正交小波变换矩阵的构造方法，是基于长度为L＝2k+1的有限长信号，Ck+1＝[ck+1，0，ck+1，1，L，ck+1，L-1]，进行如下方式对称延拓， 小波分解过程可写为COMPk=CkLDk=DEk·Ck+1,]]>其中Ck＝[ck，0，ck，1，…，ck，L/2-1]T为分解后的逼近信号，Dk＝[dk，0，dk，1，…，dk，L/2-1]T为分解后的细节信号，DEk=P~kQ~k]]>为分解矩阵， 均为2k×2k+1的矩阵；求分解矩阵DEk的逆矩阵即得重构矩阵REk=DEk-1.]]>
本发明对于9/7小波变换，其分解矩阵DEk的构造方法如下首先设置一个2k×2k+1的空矩阵；让滤波器系数 以下标递增的方式排成一行放置在空矩阵的第一行，不够则补零，其中令 放置在(0，0)处，其左方是下标为负的滤波器系数，右方是下标为正的滤波器系数；下面各行分别由上一行向右移两位得到，共移位2k-1次；然后，分别以第一列和最后一列的系数为对称轴，将放置在2k×2k+1矩阵外面的系数折叠叠加到矩阵里面即得2k×2k+1的矩阵 用同样的方法在 的下方得到2k×2k+1的矩阵 然后，由DEk=P~kQ~k]]>组成2k+1×2k+1的分解矩阵DEk。
对于长度为L＝2s+1的信号，最多可以有s+1级分解矩阵DEk，k＝0，1，…，s，设分解的深度为m，其中1≤m≤s+1，称m＝s+1为最大分解深度。
本发明给出了分解深度为m的小波分解矩阵DE可通过如下公式(1)得到DE=DEs-m+1OOIL(2m-1-1)2m-1·DEs-m+2OOIL(2m-2-1)2m-2····DEs-2OOIL(22-1)22DEs-1OOIL2[DEs]---(1)]]>相应地，由双正交小波变换理论，若分解时深度为m(1≤m≤s+1)，本发明给出相应的重构矩阵的构造公式(2)如下RE=REsREs-1OOIL2·REs-2OOIL(22-1)22LREs-m+1OOIL(2m-1-1)2m-1---(2).]]>对于9/7小波最大分解深度为3的小波分解矩阵(8×8)为
DE=0.40240.67670.52010.42680.32240.23340.16940.0773-0.5000-0.6158-0.18390.15760.51990.36040.16410.0978-0.4502-0.04430.78890.1941-0.3578-0.14130.0191-0.00850.08740.08860.0213-0.3713-0.7516-0.03490.64290.3175-0.41810.7478-0.3536-0.04070.06450000.0645-0.0407-0.41810.7885-0.4181-0.04070.06450000.0645-0.0407-04181-0.7885-0.3536-0.040700000.1291-0.0814-0.83620.7885]]>相应地，重构矩阵为RE=0.35360.35360.35360.35360.35360.35360.35360.3536-0.4776-0.3254-0.06960.19060.41560.48030.48770.4988-0.6586-0.01860.65910.2463-0.2858-0.2397-0.0955-0.07920.04350.0186-0.0440-0.2463-0.32930.14840.65300.7617-0.75480.7421-0.3536-0.07280.02380.0378000.0477-0.0728-0.37740.8527-0.3774-0.11060.02380.075700.03780.0238-0.1106-0.37740.8905-0.3536-0.22120000.03780.0238-0.1106-0.37740.8527.]]>本发明的优点和积极效果1、本发明给出了信号在对称延拓方式下，有限长信号不需逐级计算而直接得到小波系数的分解矩阵和由这些小波系数重构原信号的重构矩阵的构造方法，并给出常用的相应于9/7小波的分解矩阵和重构矩阵及其基向量，它们可广泛用于基于小波的图像分块处理中，如图像压缩、图像去噪、图像增强等。本发明还给出分解、重构矩阵的互易性，使分解重构矩阵可以互用。2、本发明构造的小波变换矩阵基矢量比周期延拓方式下构造的小波变换矩阵更适于图像分块处理，其构造方法比较直观，容易理解，不仅可简化计算，而且边界效应明显减少。3、由于自然图像块间的连续性，使本发明构造出的信号对称延拓方式下小波变换矩阵，对于信号或图像的对称延拓更具现实意义和实用价值。
图1是9/7小波变换分解基矢量(8×8)图。
具体实施方式实施例1本发明提供的对称延拓双正交小波变换矩阵的构造方法，是基于长度为L＝2k+1的有限长信号，Ck+1＝[ck+1，0，ck+1，1，L，ck+1，L-1]，进行如下方式对称延拓， 小波分解过程可写为COMPk=CkLDk=DEk·Ck+1,]]>其中Ck＝[ck，0，ck，1，…，ck，L/2-1]T为分解后的逼近信号，Dk＝[dk，0，dk，1，…，dk，L/2-1]T为分解后的细节信号，DEk=P~kQ~k]]>为分解矩阵， 均为2k×2k+1的矩阵；求分解矩阵DEk的逆矩阵即得重构矩阵REk=DEk-1.]]>实施例2本发明指出，对于长度为L＝2s+1的信号，最多可以有s+1级分解矩阵DEk(k＝0，1，…，s)，其构造方法由实施例1中给出。设分解的深度为m，其中1≤m≤s+1，称m＝s+1为最大分解深度。
本发明给出了分解深度为m的小波分解矩阵DE可通过如下公式(1)得到DE=DEs-m+1OOIL(2m-1-1)2m-1·DEs-m+2OOIL(2m-2-1)2m-2····DEs-2OOIL(22-1)22DEs-1OOIL2[DEs]---(1)]]>相应地，由双正交小波变换理论，若分解时深度为m(1≤m≤s+1)，本发明给出相应的重构矩阵的构造公式(2)如下RE=REsREs-1OOIL2·REs-2OOIL(22-1)22LREs-m+1OOIL(2m-1-1)2m-1---(2).]]>本发明构造的小波变换矩阵基矢量比周期延拓方式下构造的小波变换矩阵更适于图像分块处理，其构造方法比较直观，容易理解。由图1知，这些分解基矢量具有较好的时—频分析的特点，对于长度为L的数据向量X，设DE代表L×L的小波分解矩阵，RE代表L×L的小波重构矩阵。则其小波分解的系数向量Y可以用矩阵—向量乘法实现Y＝DE·XY的每个分量yi(i＝0，1，L，L-1)是X在基向量Vi上的投影yi=viT·X.]]>小波重构可用如下矩阵乘法实现X＝RE·Y实施例3以9/7小波为例(本发明提出的小波变换矩阵的构造方法同样适用于其他双正交小波，滤波器系数不同而已)，分解矩阵DEk的构造方法如下本发明方法是先设置一个2k×2k+1的空矩阵(如表1中整个阴影部分)，让滤波器系数 以下标递增的方式排成一行放置在空矩阵的第一行(不够则补零)，其中令 放置在(0，0)处，其左方是下标为负的滤波器系数，右方是下标为正的滤波器系数(部分系数可能出现在设置的空矩阵外面)，下面各行分别由上一行向右移两位得到(共移位2k-1次)，表1的上半部分给出了k＝2时的情况。然后，分别以第一列和最后一列的系数为对称轴，将放置在2k×2k+1矩阵外面的系数折叠叠加到矩阵里面即得2k×2k+1的矩阵 如第一行中， 叠加到 上； 叠加到 上；第四行中， 叠加到 上； 叠加到 上； 叠加到 上等。用同样的方法在 的下方得到2k×2k+1的矩阵 然后，由DEk=P~kQ~k]]>组成2k+1×2k+1的分解矩阵DEk。
表1DE2的构造过程 根据以上方法可构造出任意长度为L＝2k+1信号的小波变换矩阵DEk(例如长度L＝8，k＝2，则可构造8×8小波分解矩阵DE2)，求分解矩阵DEk的逆矩阵即得重构矩阵REk=DEk-1,]]>且二者具有互易性。
常用9/7小波最大分解深度为3的小波分解矩阵(8×8)为DE=0.40240.67670.52010.42680.32240.23340.16940.0773-0.5000-0.6158-0.18390.15760.51990.36040.16410.0978-0.4502-0.04430.78890.1941-0.3578-0.14130.0191-0.00850.08740.08860.0213-0.3713-0.7516-0.03490.64290.3175-0.41810.7478-0.3536-0.04070.06450000.0645-0.0407-0.41810.7885-0.4181-0.04070.06450000.0645-0.0407-04181-0.7885-0.3536-0.040700000.1291-0.0814-0.83620.7885]]>相应地，重构矩阵为RE=0.35360.35360.35360.35360.35360.35360.35360.3536-0.4776-0.3254-0.06960.19060.41560.48030.48770.4988-0.6586-0.01860.65910.2463-0.2858-0.2397-0.0955-0.07920.04350.0186-0.0440-0.2463-0.32930.14840.65300.7617-0.75480.7421-0.3536-0.07280.02380.0378000.0477-0.0728-0.37740.8527-0.3774-0.11060.02380.075700.03780.0238-0.1106-0.37740.8905-0.3536-0.22120000.03780.0238-0.1106-0.37740.8527]]>由双正交小波变换性质，分解、重建矩阵具有互易性。一般用带有直流分量的变换矩阵作为分解矩阵，而用其逆矩阵作为重构矩阵。所以选择RE作为9/7小波分解矩阵DE，其分解基矢量图形如下图1所示。
在滤波器长度相对较长时，构造变换矩阵时有可能需要折叠叠加多次，如表2第一行中，元素 以第一列为对称轴折叠到第四列右边，这时要再以第四列为对称轴叠加到 上。
表2DE1的构造过程 得到DE1如下DE1=P~1Q~1=p~0p~1+p~-1p~2+p~-2+p~4+p~-4p~3+p~-3p~-2+p~4p~-1+p~3+p~-3p~0+p~2+p~-4p~1q~0q~1+q~-1q~2+q~-2+q~4q~3q~-2+q~4q~-1+q~3q~0+q~2q~1.]]>
权利要求
1.一种对称延拓双正交小波变换矩阵的构造方法，其特征是本发明基于长度为L＝2k+1的有限长信号，Ck+1＝[Ck+1，0，ck+1，1，L，ck+1，L-1]，进行如下方式对称延拓， 小波分解过程可写为COMPk=CkLDk=DEk·Ck+1,]]>其中Ck＝[Ck，0，Ck，1，…，Ck，L/2-1]T为分解后的逼近信号，Dk＝[dk，0，dk，1，…，dk，L/2-1]T为分解后的细节信号，DEk=P~kQ~k]]>为分解矩阵， 均为2k×2k+1的矩阵；求分解矩阵DEk的逆矩阵即得重构矩阵REk=DEk-1.]]>
2.根据权利要求1所述的对称延拓双正交小波变换矩阵的构造方法，其特征是，本发明对于分解矩阵DEk的构造方法如下a)首先设置一个2k×2k+1的空矩阵； b)让滤波器系数 以下标递增的方式排成一行放置在空矩阵的第一行，不够则补零，其中令 放置在(0，0)处，其左方是下标为负的滤波器系数，右方是下标为正的滤波器系数；c)下面各行分别由上一行向右移两位得到，共移位2k-1次；d)然后，分别以第一列和最后一列的系数为对称轴，将放置在2k×2k+1矩阵外面的系数折叠叠加到矩阵里面即得2k×2k+1的矩阵 e)用同样的方法在 的下方得到2k×2k+1的矩阵 然后，由DEk=P~kQ~k]]>组成2k+1×2k+1的分解矩阵DEk。
3.根据权利要求1、2所述的对称延拓双正交小波变换矩阵的构造方法，其特征是，对于长度为L＝2s+1的信号，最多可以有s+1级分解矩阵DEk，k＝0，1，…，s，设分解的深度为m，其中1≤m≤s+1，称m＝s+1为最大分解深度；本发明给出了分解深度为m的小波分解矩阵DE可通过如下公式(1)得到DE=DEs-m+1OOIL(2m-1-1)2m-1·DEs-m+2OOIL(2m-2-1)2m-2····DEs-2OOIL(22-1)22DEs-1OOIL2[DEs]---(1);]]>相应地，由双正交小波变换理论，若分解时深度为m(1≤m≤s+1)，本发明给出相应的重构矩阵的构造公式(2)如下RE=REsREs-1OOIL2·REs-2OOIL(22-1)22LREs-m+1OOIL(2m-1-1)2m-1---(2);]]>对于长度为L的数据向量X，设DE代表L×L的小波分解矩阵，RF代表L×L的小波重构矩阵，则其小波分解的系数向量Y可以用矩阵—向量乘法实现Y＝DE·X；Y的每个分量yi(i＝0，1，L，L-1)是X在基向量Vi上的投影yi=ViT·X,]]>小波重构可用如下矩阵乘法实现X＝RE·Y。
4.根据权利要求3所述的对称延拓双正交小波变换矩阵的构造方法，其特征是9/7小波最大分解深度为3的小波分解矩阵(8×8)为DE=0.40240.67670.52010.42680.32240.23340.16940.0773-0.5000-0.6158-0.18390.15760.51990.36040.16410.0978-0.4502-0.04430.78890.1941-0.3578-0.14130.0191-0.00850.08740.08860.0213-0.3713-0.7516-0.03490.64290.3175-0.41810.7478-0.3536-0.04070.06450000.0645-0.0407-0.41810.7885-0.4181-0.04070.06450000.0645-0.0407-0.41810.7885-0.3536-0.040700000.1291-0.0814-0.83620.7885]]>相应地，重构矩阵为RE=0.35360.35360.35360.35360.35360.35360.35360.3536-0.4776-0.3254-0.06960.19060.41560.48030.48770.4988-0.6586-0.01860.65910.2463-0.2858-0.2397-0.0955-0.07920.04350.0186-0.0440-0.2463-0.32930.14840.65300.7617-0.75480.7421-0.3536-0.07280.02380.0378000.0477-0.0728-0.37740.8527-0.3774-0.11060.02380.075700.03780.0238-0.1106-0.37740.8905-0.3536-0.22120000.03780.0238-0.1106-0.37740.8527.]]>
全文摘要
一种对称延拓双正交小波变换矩阵的构造方法；本发明方法是基于长度为L＝文档编号G06T9/00GK1828669SQ20061001343
公开日2006年9月6日 申请日期2006年4月4日 优先权日2006年4月4日
发明者侯正信, 杨爱萍 申请人:天津大学