一种基于解析稀疏表示的图像盲分离方法
【技术领域】
[0001] 本发明属于图像处理技术领域。涉及一种图像盲分离方法。
【背景技术】
[0002] 在实际生活，我们得到的信号往往是多个信号的混合物，为了有效的提取信号中 的信息，对这些混合物进行分离是很有必要的。盲分离是指在源信号及其混合方式未知的 情况下，仅仅利用观测信号来恢复出源信号的过程。盲分离问题可以用以下模型来描述：
[0003] Y = AX+V (1)
[0004] 上式的含义是η维源信号X= [Xl，x2，...，xn] e RnXN通过mXn维混合矩阵A 混合，再叠加噪声VeRmXN得到m维混合信号γ= [yi，y2，...，yJ eRmXN。近年来信号 盲分离已经取得了很多的研究成果，提出了很多算法来解决源信号的恢复问题，如独立 分量分析（Independent Component Analysis, ICA)，稀疏分量分析（Sparse Component Analysis, SCA)和非负矩阵分解（Nonnegative Matrix Factorization, NMF)等方法。
[0005] 稀疏分量分析方法的先验知识通常是源信号的稀疏性。而现实生活中的信号大多 不是稀疏，因此对信号的稀疏表示是很有必要的。信号的稀疏表示大多采用综合稀疏模型 (Synthesis Sparse Model)，即在满足一定的近似条件下，用少量过完备字典中的原子的 线性组合来表示信号。在综合稀疏模型中，信号由少数原子构成的子空间来表示，使得信号 稀疏表示受个别原子影响较大。解析稀疏模型（Analysis Sparse Model)克服了综合稀疏 模型稀疏表示性能较差的缺点。假设输入信号X e Rk，解析字典为Ω e RpXk，也称之为解 析算子（Analysis Operator)，与综合字典不同，其行向量ω/为解析字典的原子（（·）1'表 示转置运算）。Ω -般也是过完备字典，但与综合字典不同，它的行数大于列数，即p > k。 信号X的解析稀疏模型定义为
[0006] 1: = p_ | | Ω X | |。 (2)
[0007] 上式中Ωχ为X的解析稀疏表示系数，1为共稀疏度，它是向量Ωχ中零的个数， 也就是Ω中与X正交的原子的个数，其数值越大，解析稀疏表示系数就越稀疏。在解析稀 疏模型中信号用所有与之正交的原子构成的子集来表示，所以信号稀疏表示受个别原子影 响小，具有较好的稀疏表示性能，能很好地应用于源信号的恢复问题。然而目前大多数方法 采用计算量很大的贪婪类追踪算法来估计源信号，很显然这些方法并不是最优的。
【发明内容】
[0008] 本发明的目的在于针对现有盲分离方法存在的不足，提出了一种基于解析稀疏表 示的图像盲分离方法，该方法分别利用子集追踪、分裂Bregman算法和最小二乘法学习解 析字典、恢复源信号和估计混合系统，实现混合图像的盲分离。
[0009] 为了达到上述目的，本发明采用下述技术方案。
[0010] -种基于解析稀疏表示的图像盲分离方法，其特征在于首先利用图像信号的解析 稀疏先验，采用子集追踪算法，通过学习得到源图像的解析字典，然后利用Bregman距离作 为目标函数，采用分裂Bregman算法进行恢复源信号，再利用最小二乘法来估计混合系统， 循环进行以上三个步骤得到最终源图像。
[0011] 上述图像分离方法的具体步骤如下：
[0012] ①根据已知的混合图像Y= [yi，y2，...，yJ e 假设已知的混合矩阵 A e RmXn得到初始的源图像X = [X χ2,...，XJ e RnXN，从初始源图像x]中重叠抽取出K 个#xW大小的图像块，将这些图像块按列排列得到训练数据矩阵尸L.，》。
[0013] ②对\利用子集追踪算法训练得到解析字典Ω ]ε RpXd。
[0014] ③先定义第j个源信号的第i个图像块&的剩余量
，其中，
为混合图像的图像块形式： 为源图像的 图像块形式。
[0015] 采用分裂Bregman算法恢复源图像，该算法的优化函数为：
[0017] 其中心为源图像的估计，2]1为&的解析稀疏表示，λ和μ为可调节的系数。
[0018] ④将第③步得到的所有％按列顺序存储可得到?;?，然后将& 中每列排列 为乂7X 大小的图像块，则恢复图像X]由这些图像块拼接而成。
[0019] ⑤利用这个拼接得到的源图像X]，用最小二乘法求解混合矩阵的第j列a、首
[0020] 先定义第j个源图像\的剩余量
Μ古计混合矩阵的目标函数是
[0022] 本发明方法与现有技术相比较，具有如下显著的特点和优点：
[0023] 本发明提供的基于解析稀疏的图像盲分离方法是首先利用图像信号的解析稀疏 先验，采用子集追踪算法，通过学习得到源图像的解析字典，然后利用Bregman距离作为目 标函数，采用分裂Bregman算法恢复源信号，再利用最小二乘法来估计混合系统，实现混合 图像的盲分离。具体特点和优点为：
[0024] (1)针对现有的基于稀疏表示的图像盲分离算法大多采用综合稀疏表示这一情 况，提出了基于解析稀疏表示的混合图像盲分离算法。
[0025] (2)对于学习字典存在的不足，将子集追踪应用到解析字典学习中，为学习字典提 供了比较快速有效的方法。
[0026] (3)对现有解析稀疏方法中采用的计算量大的贪婪类追踪算法的不足，将分裂 Bregman算法应用于源信号的恢复中，加快了算法的迭代速度。
[0027] (4)通过对解析字典、源信号和混合矩阵的迭代估计过程，得到最优的分离结果， 这个迭代过程提高了分离的有效性。
[0028] 本发明提供的基于解析稀疏表示的图像盲分离方法能提高分离信号的质量。在无 线电通信、雷达与声纳信号处理、医学图像分析、图像信号处理以及语音识别等领域均有广 泛的应用。
【附图说明】
[0029] 图1为在适定情况下本发明和FastICA方法以及基于综合稀疏表示的图像 盲分离方法的分离结果的比较。其中，（a)为本发明（Analysis-based)的恢复源信 号，（b)为FastICA方法的恢复源信号，（c)为基于综合稀疏表示的图像盲分离方法 (Synthesis-based)的恢复源信号。
[0030] 图2为过定情况下分别利用本发明和基于综合稀疏表示的图像盲分离方法得到 的恢复源信号。其中，左起第一和第二幅图为本发明的恢复源信号，第三和第四幅图为基于 综合稀疏表示的图像盲分离的恢复源信号。σ表示的是添加噪声的标准差，其中（a)为〇=〇，（b)为 〇 = 5，（c)为 0 = 10，（d)为 0 = 15。
【具体实施方式】
[0031] 本发明将通过以下实施例作进一步说明。
[0032] 本发明方法首先利用图像信号的解析稀疏先验，采用子集追踪算法，通过学习得 到源图像的解析字典，然后利用Bregman距离作为目标函数，采用分裂Bregman算法恢复源 信号，再利用最小二乘法来估计混合系统，通过对以上三个步骤的迭代求解过程，得到最终 的分离源图像，达到盲分离的目的。
[0033] 具体步骤为：
[0034] ①从由混合图像估计得到的源图像中重叠抽取出K个忍x#大小的图像块，将 各图像块按列排列得到训练数据矩阵
[0035] ②对&利用子集追踪算法训练得到解析字典Ω , e R pXd。
[0036] ③采用分裂Bregman算法估计源图像，该算法的优化函数为：
[0038] 其中
是第j个源图像的第i个图像块％的剩余量，％为源图 像的估计，2]1为的解析稀疏表示，λ为可调节的系数。
[0039] ④将第③步得到的所有1^按列顺序存储可得到\ 。
[0040] ⑤将$中每列排列为vQv W大小的图像块，估计源图像由这些图像块拼接而 成，重叠的区域由平均得到。
[0041] ⑥利用这个拼接得到的源图像X]，用最小二乘法求解混合矩阵的第j列a、首
[0042] 先定义第j个源图像\的剩余量°估计混合矩阵的目标函数 为：
[0044] ⑦对以上的②③④⑤⑥步骤循环迭代，得到最优的分离源图像。
[0045] 为评价盲分离方法的性能，采用了峰值信噪比（PSNR)衡量估计源信号的质量。
[0046] 在适定情况下，对三幅源图像（Paraty、Lena和Peppers)的混合图像进行分离，得 到表1中的结果。从表1中可以看出，本盲分离方法取得了很好的分离效果，有效地从混合 图像中分离出源信号。从表1可知，本发明的在适定情况下的分离效果优于FastICA方法 和基于综合稀疏表示的图像盲分离方法。
[0047] 在过定情况下，对两幅源图像（Paraty和Lena)的混合图像进行分离，得到表2中 的结果。由表2可以看出，本发明对于含有噪声的图像盲分离也能取得很好的分离效果。
[0048] 表1适定情况下估计源信号PSNR比较（dB)
[0050] 表2过定情况下估计源信号PSNR比较（dB)
【主权项】
1. 一种基于解析稀疏表示的图像盲分离方法，其特征是首先利用图像信号的解析稀疏 先验，采用子集追踪算法，通过学习得到源图像的解析字典；然后利用Bregman距离作为目 标函数，采用分裂Bregman算法进行恢复源信号；再利用最小二乘法来估计混合系统，循环 进行以上三个步骤得到最终源图像。2. 根据权利要求1所述的基于解析稀疏表示的图像盲分离方法，其特征是按如下步 骤： ① 根据已知的混合图像和假设已知的混合矩阵A e 到初始的源图像从初始源图像Xj中重叠抽取出K个大小的图像块，将这些图像块按列排列得到训练数据矩阵② 对$利用子集追踪算法训练得到解析字典③ 先定义第j个源信号的第i个图像块&的剩余量其中，为混合图像的图像块形式，为源图像的 图像块形式； 采用分裂Bregman算法恢复源图像，其优化函数为：其中％·为源图像的估计，2]1为的解析稀疏表示，λ和μ为可调节的系数； ④ 将第③步得到的所有' ?,χ按列顺序存储可得到\,，然后将中每列排列为大小的图像块，则恢复图像Χ]由这些图像块拼接而成； ⑤ 利用这个拼接得到的源图像\，用最小二乘法求解混合矩阵的第j列a、];首先定义 第j个源图像\的剩余量估计混合矩阵的目标函数是：
【专利摘要】一种基于解析稀疏表示的图像盲分离方法，首先利用图像信号的解析稀疏先验，采用子集追踪算法，通过学习得到源图像的解析字典；然后利用Bregman距离作为目标函数，采用分裂Bregman算法进行恢复源信号；再利用最小二乘法来估计混合系统，循环进行以上三个步骤得到最终源图像。本发明为学习字典提供了比较快速有效的方法；将分裂Bregman算法应用于源信号的恢复中，加快了算法的迭代速度；通过对解析字典、源信号和混合矩阵的迭代估计过程，得到最优的分离结果，这个迭代过程提高了分离的有效性。在无线电通信、雷达与声纳信号处理、医学图像分析、图像信号处理以及语音识别等领域均有广泛的应用。
【IPC分类】G06T5/00
【公开号】CN105354807
【申请号】CN201510815145
【发明人】张烨, 方婉婷, 徐标, 张文全
【申请人】南昌大学
【公开日】2016年2月24日
【申请日】2015年11月23日