发明名称 基于Surfacelet变换的系数自适应收缩视频去噪方法
摘要 本发明公开了一种基于Surfacelet变换的系数自适应收缩视频去噪方法,以改善视频图像边缘模糊和噪声去除不充分等现象,提高视频去噪效果,主要用于视频图像的去噪和三维立体图像的去噪。其实现过程是:(1)输入含噪视频序列,并对其进行Surfacelet变换(ST);(2)对Surfacelet变换域中每一层各个方向子带的ST系数计算阈值;(3)对各方向子带中ST系数进行分类得到相应的掩模;(4)计算各方向子带中ST系数的能量比和先验比,然后计算ST系数的收缩因子;(5)对ST系数进行收缩处理;(6)对收缩后的ST系数进行逆Surfacelet变换,得到去噪后的视频序列。本发明与现有的技术相比,明显地提高了对含噪视频中噪声的抑制能力,同时也能够很好地保留视频图像的纹理细节和边缘信息。
申请公布号 CN102196155A 申请公布日期 2011.09.21
申请号 CN201110081453.8 申请日期 2011.04.01
申请人 西安电子科技大学 发明人 田小林;焦李成;李杰;张小华;王爽;钟桦;于昕;吴建设
分类号 H04N5/21(2006.01)I;H04N5/213(2006.01)I;H04N7/26(2006.01)I;H04N13/00(2006.01)I 主分类号 H04N5/21(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 程晓霞;王品华
主权项 一种基于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>&le;</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>&xi;</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>&CenterDot;</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>&xi;</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>&eta;</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>&xi;</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>&eta;</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>&le;</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逆变换,得到去噪后的视频序列。
地址 710071 陕西省西安市太白南路2号