[发明专利]基于非局部小波系数收缩的图像去噪方法

权利要求书

1.一种基于非局部小波系数收缩的图像去噪方法，包括如下步骤：

(1)在输入含噪图像中以步长3取参考图像块，根据欧式距离公式计算该参考块与其邻域内所有图像块的距离d(Z_i,Z_i,j)，选取距离最小的N₂个图像块构成该参考块的相似组：S_i＝{Z_i,j:min(d(Z_i,Z_i,j)),j＝1...N₂,i表示相似组序号}，其中，Z_i为参考图像块，Z_i,j为Z_i邻域中的图像块，N₂为相似组中相似块个数；

(2)对相似组中的相似块进行二维小波变换,得到相似块的小波系数：

α_i,j＝Τ_2D(Z_i,j),Z_i,j∈S_i

其中，Τ_2D表示二维小波变换，α_i,j为第i个相似组中第j个相似块的小波系数；

(3)根据非局部均值公式，计算每个相似组小波系数的非局部均值μ_i：

μi=Σj∈Siωi,jαi,j,j=1...N2]]>

ω_i,j＝exp(-d(Z_i,Z_i,j))/h)/W

其中，ω_i,j为相似块对应的权值，h为12σ_n，σ_n为噪声标准差，W表示归一化操作；

(4)使用双L₁范数模型的收缩函数，计算相似块的小波系数估计值

α^i,j=Sτ1,τ2,μi(αi,j)μi≥0-Sτ1,τ2,-μi(αi,j)μi<0]]>

其中，为收缩函数，其定义为：

Sτ1,τ2,b(t)=t+τ1+τ2t<-τ1-τ20-τ1-τ2≤t≤τ1-τ2t-τ1+τ2τ1-τ2≤t≤τ1-τ2+bbτ1-τ2+b≤t≤τ1+τ2+bt-τ1-τ2t>τ1+τ2+b]]>

其中，t为待收缩系数，b是待收缩系数对应的非局部均值，τ₁和τ₂为两个不同的收缩阈值，c₁，c₂为两个不同数值的常数，c₁取值为0.1，c₂取值为0.9，σ_n为噪声标准差，σ_i是由α_i估计得到的标准差，α_i＝{α_i,1...α_i,j}，δ_i由α_i-μ_i估计得到的标准差；

(5)对小波系数估计值进行小波反变换，得到相似块估计值，整合所有估计值，得到一次估计图像Y：

Y^i,j=T2D-1(α^i,j)]]>

其中，X为含噪图像，x为X中的像素点，表示二维小波反变换，χ_j是图像块的特征函数，时，χ_j取值为1，否则为0；

(6)对得到的一次估计图像Y进行残差补回，在残差补回图上重复步骤(1)-(5)，得到基础估计图像Y'；当σ_n＞50时，对基础估计图像Y'再执行一次相同的残差补回操作。

(7)对得到的基础估计图像Y'，以步长N_s选取参考块，根据欧式距离公式构建其相似组，并记录相似组中相似块的坐标信息：

Si1={Y^i,j':||Y^i'-Y^i,j'||22(N1')2<τ}]]>

其中，和分别是基础估计图像中的参考块和候选块，表示的L₂范数平方，N₁'×N₁'是图像块大小，τ是判断两个图像块是否相似的阈值，σ_n≤40时，τ取值为400，σ_n＞40时，τ取值为3500；

(8)根据相似块坐标信息，从含噪图像中提取相应的图像块，构成相似组对基础估计图像的相似组和含噪图像的相似组均进行三维变换，得到各自的变换系数，根据相似组的变换系数计算维纳收缩系数W_i，再根据相似组的变换系数与维纳收缩系数W_i，得到相似组中相似块估计值：

Wi=|T3D(Si1)|2|T3D(Si1)|2+σn2]]>

Y^Si2'=T3D-1(WiT3D(Si2))]]>

其中，Τ_3D表示由二维DCT变换和块间一维小波变换组成的三维变换，表示三维反变换，表示对相似组的变换系数取绝对值平方，W_i为相似组对应的维纳收缩系数；

(9)对得到的相似组中相似块估计值加权平均，得到去噪图像Y''：

其中，ω_i是相似组对应的权值，表示W_i的L₂范数平方。

2.根据权利要求1所述的去噪方法，其中步骤(1)所述的根据欧式距离公式计算参考块与其邻域内所有图像块的距离d(Z_i,Z_i,j)，其公式如下：

d(Zi,Zi,j)=||Zi-Zi,j||22(N1)2]]>

其中，Z_i为参考图像块，Z_i,j为Z_i邻域中的图像块，表示Z_i-Z_i,j的L₂范数平方，N₁×N₁为图像块大小。