[发明专利]一种基于3D协同滤波与低秩矩阵重建的视频去噪方法及系统

权利要求书

1.一种基于3D协同滤波与低秩矩阵重建的视频去噪方法，其特征在于，包括如下步骤：

步骤1：对含噪视频进行3D协同滤波预处理，得到基础估计图像；所述的步骤1包括如下步骤；

步骤1.1：对输入视频的每一帧进行分块，对每个图像块分别进行帧内和帧间的相似块匹配，得到每个参考块x_R的相似块分组；

步骤1.2：对每个参考块x_R的相似块分组进行3D协同滤波，得到每个参考块x_R的估计值

步骤1.3：对估计值进行聚集，得到基础估计图像；

步骤2：对得到的基础估计图像再进行低秩矩阵重建进一步去噪，得到最终去噪后的图像；所述的步骤2包括如下步骤；

步骤2.1：对步骤1得到的基础估计图像的每个图像块分别进行帧内和帧间的相似块匹配，得到每个图像块的相似块分组；将分组中的每一个相似块排列成向量，组成一个相似性块矩阵Y；

Y＝(y₁,y₂,…y_m)

其中m为参考块对应相似块组的个数，y₁,y₂,…y_m分别为大小为N₁×N₁的匹配块矩阵；

步骤2.2：将矩阵Y分解为两个矩阵A和E之和，对相似性块矩阵Y进行低秩循环最小化逼近，得到降噪之后的低秩矩阵A；

Y＝A+E

其中矩阵A和E均未知，A为低秩的无噪数据矩阵，E为噪声稀疏矩阵；

所述的2.2中所述的对相似性块矩阵Y进行低秩循环最小化逼近采用非精确拉格朗日乘子法IALM低秩逼近，得到降噪后的低秩矩阵A，包括如下步骤进行：

步骤2.2.1：更新噪声稀疏矩阵E；

Ek+1=D-Ak+Ykμ-PΩ∞λ/μk(D-Ak+Ykμ)]]>

其中，n为参考块大小N₁×N₁×3，σ为噪声水平标准差，m为参考块对应的相似块组的个数，D为相似块矩阵Y的初始值，表示在空间上}的欧几里得投影，λ＝100/(100-σ)，R^m×n表示实数矩阵，k代表当前的迭代循环次数，k为初始值为0的自然数，A_k,Y_k分别表示第k次循环中低秩矩阵A和相似性块矩阵Y的值，A_k,Y_k的初始值为0，E_k+1为新更新的噪声稀疏矩阵；

步骤2.2.2：更新低秩矩阵A，求得更新后低秩矩阵A_k+1，公式如下：

Ak+1=Uk+1Sτ(Σ)Vk+1T]]>

其中U_k+1，∑，通过对矩阵进行奇异值分解SVD求得，其中，U_k+1为左奇异矩阵，∑为奇异值对角矩阵，为右奇异矩阵，V_k+1^T表示矩阵V_k+1的转置，S_τ(∑)表示对奇异值对角矩阵∑进行软阈值收缩；

阈值计算公式为：

τ＝0.5m·n/sum(D)

其中，sum(D)为矩阵D中所有元素的和，矩阵D为相似性块矩阵Y的初始值；

步骤2.2.3：根据更新后低秩矩阵A_k+1，求得更新的矩阵Y_k+1，同时令k＝k+1进行下一次迭代；

Y_k+1＝Y_k+μ(D-A_k+1-E_k+1)

其中A_k+1,E_k+1分别为步骤2.2.1和2.2.2中更新的低秩矩阵A和噪声稀疏矩阵E的值，Y_k为第k次循环中相似性块矩阵的Y值；

步骤2.2.4判断迭代次数k是否达到最大迭代次数，若达到最大迭代次数，则终止迭代，输出去噪低秩矩阵A；否则重复步骤2.2.1到步骤2.2.3，进行下一次迭代。

步骤3.根据步骤2输出的降噪之后的低秩矩阵A，采用均值法对重叠的像素进行聚集，最终输出视频。