基于Context模型和双树复小波变换的水下声纳图像的去噪方法技术

本发明专利技术的目的在于提供基于Context模型和双树复小波变换的水下声纳图像的去噪方法，包括如下步骤：对水下声纳图像进行双树复小波分解，图像经四层双树复小波分解后获得的低频近似分量保持不变，对图像的高频分量进行去噪处理，对处理后的复小波系数进行双树复小波反变换，获得最终去噪后的图像。本发明专利技术采用Context模型衡量双树复小波系数之间能量的相似性，将能量相接近的系数归类，针对每一类系数确定不同的阈值，并结合软阈值函数实现噪声的去除。优化了阈值的选取，在去除噪声的同时保留了更多的图像细节，抑制了系数被过扼杀的现象。

【技术实现步骤摘要】
本专利技术涉及的是一种图像处理
的去噪方法。
技术介绍
图像在采集和传输的过程中往往会受到外界环境的干扰，造成图像质量退化，这严重影响了后续的图像处理工作。因而图像的去噪预处理也成为了图像处理的研究热点，它为后续高层次的图像处理提供了更好的研究基础。图像去噪的方法的种类很多，主要分为空间域方法和变换域方法。空间域方法是 针对图像的灰度值直接进行运算处理，常用的方法有中值滤波、均值滤波、非局部均值滤波等等。中值滤波和均值滤波都会存在一定程度的边缘模糊现象，非局部均值滤波的计算量较大。变换域方法中常用的方法是基于小波的萎缩阈值法，传统的全局阈值对每一个变换系数并不是最优的，导致图像的细节信息丢失。
技术实现思路
本专利技术的目的在于提供具备良好的平移不变性和方向选择性的基于Context模型和双树复小波变换的水下声纳图像的去噪方法。本专利技术的目的是这样实现的本专利技术基于Context模型和双树复小波变换的水下声纳图像的去噪方法，其特征是(I)对水下声纳图像进行双树复小波分解，分解层数为四层图像经过四层双树复小波变换，获得一个低频近似分量和四个尺度下的高频分量，yl代表经过4层分解后图像的低频近似分量的系数矩阵，yhk分别表示第一、二、三、四层的高频分量，k = I, 2，3，4，每一层高频分量包含6个复系数矩阵，分别代表该尺度下双树复小波变换分解的6个方向的高频细节信息；(2)保留图像分解后的低频分量yl ；(3)对图像分解的四层高频分量yhp yh2、yh3和yh4进行去噪处理I)图像的第一层高频分量yhi和第二层高频分量yh2，采用非线性软阈值方法进行去噪处理；2)图像的第三层高频分量yh3和第四层高频分量yh4，采用Context模型对系数进行分类，针对每一类系数确定各自的阈值，结合软阈值函数进行去噪处理①获取图像的第三层高频分量yh3中六个复系数矩阵的实部和虚部系数；图像的第三层高频分量yh3同样由六个复系数矩阵组成，每一个系数矩阵代表着该尺度下图像6个方向的高频信息，yh3 (d) (d= 1，2，...，6)表示该尺度下任意一个方向的复系数矩阵，它的实部和虚部系数分别为a3 (d) = real (yh3 (d)) (d = I, 2, . . . , 6)b3 (d) = imag (yh3 (d)) (d = 1,2,...,6)a3(d)和133((1)表示第三层高频分量的任意一个方向复系数的实数和虚数部分，实现后续对获得的实部和虚部分别进行处理；②利用Context模型对系数进行分类；用Context模型来衡量变换系数的能量的接近程度，首先计算出每个系数的Context 值 Z(i, j)z (i, j) = 2 X (IY (i-i, j-i) I +1Y (i-i, j) I +1Y (i-i, j+i) I +1Y (i, j-i) I +1Y (i, j)+ |Y(i, j+1) +|Y(i+l，j-1) +|Y(i+l，j) +|Y(i+l，j+1) 1))/9式中i，j表示变换系数在矩阵中的位置，Y表示实数部分或虚数部分，利用上式计算所有方向的复系数的实部和虚部的Context值,将Context值分类,在Context值小于20时，每隔5分一级，20以后增大分级间隔；③确定阈值选取经Context模型分类后a3(4)系数进行阈值估计，确定每一类系数的阈值首先估计噪声的标准差σ n，计算噪声的能量σ η2，然后求取每一类系数的方差，再根据下式获得信号的能量σχ2:σχ2 = var (y) - σ n2式中y是每一类系数，var (y)是求解y的方差，最后采用Bayes阈值估计方法确定每一类系数的阈值T :T = Txafl /σχ式中r是可调节参数；④对任意方向经Context模型分类后的变换系数的实部和虚部分别采用软阈值函数式进行滤波处理；每一类系数所选取的阈值由式/’ = Γχσ; /σχ确定，经过处理后的实部和虚部系数表示为 al3 (d) (d = I, 2, . . , 6)和 bl3 (d) (d = 1,2,...,6)；⑤将处理后的实部和虚部系数重新组合，获得处理后的复系数，复系数表达式为yh3 (d) = al3 (d) +bl3 (d) X i (d = I, 2, . . . , 6)式中i表示虚数单位；⑥对图像分解的第四层高频分量按照步骤①-⑤进行处理；(4)对处理后的复小波系数进行双树复小波反变换，获得最终去噪后的图像。本专利技术还可以包括I、所述的非线性软阈值方法为(I)获取图像的第一层高频分量Yh1中六个复系数矩阵的实部和虚部系数图像的第一层高频分量Yh1由六个复系数矩阵组成，分别代表着该尺度下图像6个方向的高频信息，yhjd) (d= 1,2,...,6)表示该尺度下任意一个方向的复系数矩阵，它的实部和虚部系数分别为H1 (d) = real (Yh1 (d)) (d = 1,2, , 6) h1 (d) = imag Cyh1 (d)) (d = I, 2, . . . , 6)式中real和imag表示求取实数和虚数部分,S1 (d)和Id1 (d)表示第一层高频分量的任意一个方向复系数的实数和虚数部分；(2)确定阈值T :选取ai (4)系数进行阈值估计，阈值通过非线性软阈值获得Γ = σ XIog(AZ)式中= η= median(\a{ (4)|) / 0.6745 <是 σ η 的估计；Ν 是 ⑷中含有的系数个数；(3)采用软阈值函数对变换系数的实部ai(d) (d= 1,2,...,6)和虚部bJdMd = 1，2, · · ·，6)分别进行去噪处理软阈值函数表达式为 f_(x)x(|x|-r)，Ixl >ΓJ= Μ 'τ[O,\χ\<τ式中X是需要进行阈值处理的变换系数，对任意方向的变换系数的实部和虚部分别采用软阈值函数进行滤波处理，经过处理后的实部和虚部系数表示为aljd) (d =1，2，· · ·，6)和 M1(Cl) (d = 1，2，· · ·，6)；(4)将处理后的实部和虚部系数重新组合，获得处理后的复系数复系数表达式为Yh1 (d) = ali (d) +M1 (d) X i (d = I, 2, . . . , 6)式中i表示虚数单位；(5)对图像分解的第二层高频分量按照步骤(I) - (4)进行处理。本专利技术的优势在于本专利技术采用Context模型衡量双树复小波系数之间能量的相似性，将能量相接近的系数归类，针对每一类系数确定不同的阈值，并结合软阈值函数实现噪声的去除。优化了阈值的选取，在去除噪声的同时保留了更多的图像细节，抑制了系数被过扼杀的现象。附图说明图I为本专利技术的流程图；图2为一维双树复小波分解滤波器结构图；图3a为噪声方差σ = 10时的含噪图像；图3b为噪声方差σ = 10时采用本专利技术方法的去噪后图像；图4a为噪声方差σ = 15时的含噪图像；图4b为噪声方差σ = 15时采用本专利技术方法的去噪后图像；图5a为噪声方差σ = 20时的含噪图像；图5b为噪声方差σ = 20时采用本专利技术方法的去噪后图像；图6不同噪声方差下去噪前后图像的峰值信噪比(PSNR)对比图。具体实施例方式下面结合附本文档来自技高网...

【技术保护点】
基于Context模型和双树复小波变换的水下声纳图像的去噪方法，其特征是：（1）对水下声纳图像进行双树复小波分解，分解层数为四层：图像经过四层双树复小波变换，获得一个低频近似分量和四个尺度下的高频分量，yl代表经过4层分解后图像的低频近似分量的系数矩阵，yhk分别表示第一、二、三、四层的高频分量，k＝1,2,3，4，每一层高频分量包含6个复系数矩阵，分别代表该尺度下双树复小波变换分解的6个方向的高频细节信息；（2）保留图像分解后的低频分量yl；（3）对图像分解的四层高频分量yh1、yh2、yh3和yh4进行去噪处理：1）图像的第一层高频分量yh1和第二层高频分量yh2，采用非线性软阈值方法进行去噪处理；2）图像的第三层高频分量yh3和第四层高频分量yh4，采用Context模型对系数进行分类，针对每一类系数确定各自的阈值，结合软阈值函数进行去噪处理：①获取图像的第三层高频分量yh3中六个复系数矩阵的实部和虚部系数；图像的第三层高频分量yh3同样由六个复系数矩阵组成，每一个系数矩阵代表着该尺度下图像6个方向的高频信息，yh3(d)(d＝1,2,..,6)表示该尺度下任意一个方向的复系数矩阵，它的实部和虚部系数分别为：a3(d)＝real(yh3(d))(d＝1,2,...,6)b3(d)＝imag(yh3(d))(d＝1,2,...,6)a3(d)和b3(d)表示第三层高频分量的任意一个方向复系数的实数和虚数部分，实现后续对获得的实部和虚部分别进行处理；②利用Context模型对系数进行分类；用Context模型来衡量变换系数的能量的接近程度，首先计算出每个系数的Context值Z(i,j)：Z(i,j)＝2×(|Y(i？1,j？1)|+|Y(i？1,j)|+|Y(i？1,j+1)|+|Y(i,j？1)|+|Y(i,j)|+|Y(i,j+1)|+|Y(i+1,j？1)|+|Y(i+1,j)|+|Y(i+1,j+1)|))/9式中i,j表示变换系数在矩阵中的位置，Y表示实数部分或虚数部分，利用 上式计算所有方向的复系数的实部和虚部的Context值，将Context值分类，在Context值小于20时，每隔5分一级，20以后增大分级间隔；③确定阈值：选取经Context模型分类后a3(4)系数进行阈值估计，确定每一类系数的阈值：首先估计噪声的标准差σn，计算噪声的能量σn2，然后求取每一类系数的方差，再根据下式获得信号的能量σx2＝var(y)？σn2式中y是每一类系数，var(y)是求解y的方差，最后采用Bayes阈值估计方法确定每一类系数的阈值T：T=r×σn2/σx式中r是可调节参数；④对任意方向经Context模型分类后的变换系数的实部和虚部分别采用软阈值函数式进行滤波处理；每一类系数所选取的阈值由式确定，经过处理后的实部和虚部系数表示为a13(d)(d＝1,2，...,6)和b13(d)(d＝1,2，...,6)；⑤将处理后的实部和虚部系数重新组合，获得处理后的复系数，复系数表达式为：yh3(d)＝a13(d)+b13(d)×i(d＝1,2,...,6)式中i表示虚数单位；⑥对图像分解的第四层高频分量按照步骤①？⑤进行处理；（4）对处理后的复小波系数进行双树复小波反变换，获得最终去噪后的图像。FDA00002186697400021.jpg,FDA00002186697400023.jpg...

【技术特征摘要】
1.基于Context模型和双树复小波变换的水下声纳图像的去噪方法，其特征是 (1)对水下声纳图像进行双树复小波分解，分解层数为四层 图像经过四层双树复小波变换，获得一个低频近似分量和四个尺度下的高频分量，yl代表经过4层分解后图像的低频近似分量的系数矩阵，yhk分别表示第一、二、三、四层的高频分量，k = 1，2，3，4，每一层高频分量包含6个复系数矩阵，分别代表该尺度下双树复小波变换分解的6个方向的高频细节信息； (2)保留图像分解后的低频分量yl； (3)对图像分解的四层高频分量yhpyh2、yh3和yh4进行去噪处理 .1)图像的第一层高频分量Yh1和第二层高频分量yh2，采用非线性软阈值方法进行去噪处理； .2)图像的第三层高频分量yh3和第四层高频分量yh4，采用Context模型对系数进行分类，针对每一类系数确定各自的阈值，结合软阈值函数进行去噪处理 ①获取图像的第三层高频分量yh3中六个复系数矩阵的实部和虚部系数； 图像的第三层高频分量yh3同样由六个复系数矩阵组成，每一个系数矩阵代表着该尺度下图像6个方向的高频信息，yh3 (d) (d= 1,2,..,6)表示该尺度下任意一个方向的复系数矩阵，它的实部和虚部系数分别为a3 (d) = real (yh3 (d)) (d = 1,2,..., 6)b3 (d) = imag (yh3 (d)) (d = 1,2,..., 6) a3(d)和133((1)表示第三层高频分量的任意一个方向复系数的实数和虚数部分，实现后续对获得的实部和虚部分别进行处理； ②利用Context模型对系数进行分类； 用Context模型来衡量变换系数的能量的接近程度，首先计算出每个系数的Context值 Z(i，j)Z(i, j) = 2X (|Y(i-l, j-1) | + |Y(i_l，j) | + |Y(i_l，j+1) | + |Y(i，j_l) | + |Y(i，j) | + |Y(i,j+1) +|Y(i+l，j-1...

【专利技术属性】
技术研发人员：李一兵，张静，叶方，汤春瑞，付强，李一晨，李骜，田雪怡，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：