主权项 |
一种基于Surfacelet变换的系数自适应收缩视频去噪方法,其特征在于:包含以下步骤:步骤1.输入含噪视频序列,并对其进行Surfacelet变换,分解层数为4层,每层对应的方向数分别为192、192、48、12;步骤2.对Surfacelet变换域中的系数按如下公式计算最精细层到最粗糙层的各方向子带的阈值:Tl,k=3·δn·El,k式中:Tl,k为第l层k方向子带中ST系数对应的三维矩阵阈值;δn为噪声方差;El,k为高斯噪声在ST域第l层k方向子带中的概率分布特征;步骤3.利用步骤2中得到的阈值对各层各个方向子带的ST系数进行分类得到相应的掩模,计算公式如下: <mrow> <msub> <mi>x</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> <mo>,</mo> <msub> <mi>a</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>></mo> <msub> <mi>T</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> <mo>,</mo> <msub> <mi>a</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>≤</mo> <msub> <mi>T</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>式中:xl,k(m,n,p)为第l层k方向子带中(m,n,p)位置上ST系数对应的掩模;al,k(m,n,p)=|wl,k(m,n,p)|为ST系数的幅值,式中wl,k(m,n,p)为第l层k方向子带在位置(m,n,p)上的ST系数;步骤4.利用步骤2中得到的阈值计算ST系数的能量比: <mrow> <msub> <mi>ξ</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>a</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>2</mn> <mo>·</mo> <msub> <mi>T</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>a</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>式中:ξl,k(m,n,p)反映第l层k方向子带中(m,n,p)位置上ST系数的能量比;步骤5.利用步骤3中得到的掩模计算ST系数的先验比;步骤6.利用步骤4中得到的能量比和步骤5中得到的先验比按以下公式计算ST系数的收缩因子: <mrow> <msub> <mi>q</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mfrac> <mrow> <msub> <mi>ξ</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> <msub> <mi>η</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <msub> <mi>ξ</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> <msub> <mi>η</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <msub> <mi>a</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>≤</mo> <msub> <mrow> <mn>2</mn> <mi>T</mi> </mrow> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> <mo>,</mo> <msub> <mi>a</mi> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>></mo> <msub> <mrow> <mn>2</mn> <mi>T</mi> </mrow> <mrow> <mi>l</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>式中:ql,k(m,n,p)为第l层k方向子带中(m,n,p)位置上ST系数对应的收缩因子;ξl,k(m,n,p)为第l层k方向子带中(m,n,p)位置上ST系数的能量比;ηl,k(m,n,p)为第l层k方向子带中(m,n,p)位置上ST系数的先验比;步骤7.利用收缩因子对ST系数按以下收缩公式进行收缩处理,得到收缩后的ST系数为:yl,k(m,n,p)=wl,k(m,n,p)·ql,k(m,n,p):式中:yl,k(m,n,p)为第l层k方向子带中(m,n,p)位置上收缩后的ST系数;wl,k(m,n,p)为第l层k方向子带中(m,n,p)位置上的ST系数;ql,k(m,n,p)为第l层k方向子带中(m,n,p)位置上ST系数对应的收缩因子;步骤8.对收缩后的ST系数进行Surfacelet逆变换,得到去噪后的视频序列。 |