首页 / 专利库 / 信号处理 / 稀疏变换 / 一种基于加权L1范数的卷积稀疏编码方法

一种基于加权L1范数的卷积稀疏编码方法

阅读:892发布:2020-05-11

专利汇可以提供一种基于加权L1范数的卷积稀疏编码方法专利检索,专利查询,专利分析的服务。并且一种基于加权L1范数的卷积稀疏编码方法,包括如下步骤,对字典 过滤器 {dm}和目标图像s进行傅里叶变换,并计算L1范数的自适应权重wm;根据交替方向乘子法构建卷积稀疏编码优化模型,将优化问题分解为x子问题和y子问题;根据卷积定律将x子问题转化为频域;通过变量重排列将x子问题的求解分解为若干个独立的线性问题进行求解;通过软 阈值 法对y子问题进行求解;根据获取的x子问题和y子问题的求解更新对偶变量um,计算出原始残差rp和对偶残差rd及计算原始停止条件∈pri和对偶停止条件∈dua,更新惩罚系数ρ,当循环达到最大 迭代 次数或原始残差和对偶残差满足停止条件rp≤∈pri且rd≤∈pri时,输出卷积稀疏特征图,否则循环上述步骤。本方法具有编码效果好,计算速度快的优点。,下面是一种基于加权L1范数的卷积稀疏编码方法专利的具体信息内容。

1.一种基于加权L1范数的卷积稀疏编码方法,其特征在于,包括如下步骤,步骤S1、对字典过滤器{dm}和目标图像s进行傅里叶变换,并计算L1范数的自适应权重wm;
步骤S2、根据交替方向乘子法构建卷积稀疏编码优化模型,将优化问题分解为x子问题和y子问题;
步骤S3、根据卷积定律将x子问题转化为频域;
步骤S4、通过变量重排列将x子问题的求解分解为若干个独立的线性问题进行求解;
步骤S5、通过软阈值法对y子问题进行求解;
步骤S6、根据获取的x子问题和y子问题的求解更新对偶变量um,计算出原始残差rp和对偶残差rd及计算原始停止条件∈pri和对偶停止条件∈dua,更新惩罚系数ρ,当循环达到最大迭代次数或原始残差和对偶残差满足停止条件rp≤∈pri且rd≤∈pri时,输出卷积稀疏特征图,否则调转到步骤S4。
2.根据权利要求1所述的一种基于加权L1范数的卷积稀疏编码方法,其特征在于,所述步骤S1中,对给定的m个字典过滤器{dm}和目标图像s做傅里叶变换,获取 和 其中,目标图像s转化为一个N维向量,N是目标图像s的像素数量;再对 和 的商做傅里叶反变换,获取 其中, 表示傅里叶反变换;则可以计算出自适应权重
其中,γ是给定的权重参数。
3.根据权利要求1所述的一种基于加权L1范数的卷积稀疏编码方法,其特征在于,所述步骤S2中的卷积稀疏编码优化模型为 其中,dm为
一组字典过滤器,s为目标图像,*为卷积操作,xm为一组卷积稀疏特征图,λ为正则化参数,⊙为哈达玛积, 为向量x的L2范数的平方,即向量x各元素的平方和,||x||1为向量x的L1范数。
4.根据权利要求1所述的一种基于加权L1范数的卷积稀疏编码方法,其特征在于,所述步骤S3中,将x子问题转化为频域的具体方法为:首先引入辅助变量zm=ym-um,同时,使用dm的对称托普利兹矩阵与xm的矩阵乘法代替dm和xm的卷积,即Dmxm=dm*xm,通过卷积定理可知函数卷积的傅里叶变换是函数傅里叶变换的乘积,则将dm,Dm,xm,s,zm经傅里叶变换为再利用卷积定理和对称托普利兹矩阵的性质进行变换:
其中, 为傅里叶变换,diag为将N维
向量转化为N×N的对矩阵;变换后可知,x子问题等价于:
其中, 为由M个N×N的对角
矩阵串联而成的对角矩阵,M为字典过滤器的个数,N为输入图像s的维度;则可通过线性系统求解: 其中 为MN×MN的矩阵, 中非零元素在第n列
的行只会与 中非零元素在第n行的列有非零乘积。
5.根据权利要求1所述的一种基于加权L1范数的卷积稀疏编码方法,其特征在于,所述步骤S4的具体方法为:在线性系统 中,右侧第m块为
则有 将向量x的第n个元素作为x(n),则有
向量an为 的第n行的所有非零元素,N个M×M的独立线性系
统为 即
通过Sherman-Morrison公式 可得
将N个独立线性体统全部求解后可得到 再利用傅里叶反变换
获取新的x。
6.根据权利要求1所述的一种基于加权L1范数的卷积稀疏编码方法,其特征在于,所述步骤S5中,通过软阈值法求解为 其中,sign(·)为符号函
数,|·|取绝对值,则Y子问题的解为
7.根据权利要求1所述的一种基于加权L1范数的卷积稀疏编码方法,其特征在于,所述步骤S6中,根据求得的ym和xm,应用公式um=um+xm-ym更新um,原始残差rp=||x-y||2和对偶残差rd=ρ||yprev-y||2,其中yprev等于上一轮的y;原始停止条件
和对偶停止条件 其中,∈abs和
∈rel为给定的绝对停止容忍度和相对停止容忍度;采用自适应策略更新其中,τ和μ为常数,常用的取值为τ=2,μ=10。

说明书全文

一种基于加权L1范数的卷积稀疏编码方法

技术领域

[0001] 本发明属于图像处理技术领域,具体涉及一种基于加权L1范数的卷积稀疏编码方法。

背景技术

[0002] 稀疏编码是一种被广泛应用于信号和图像处理应用的技术,在稀疏编码模型中,输入元素可看作超完备字典中多个原子的线性组合,其中只有部分原子对应的系数非零(当字典中原子的个数大于输入元素的维数时,则称该字典为超完备或冗余字典)。对输入元素的稀疏编码就是要在超完备字典中寻找尽可能少的原子来线性表示输入元素。
[0003] 稀疏编码问题符合人类感知系统中的稀疏编码策略,能够捕捉到数据内在的特性,因此得到了众多研究者的关注以及在计算机视觉领域广泛的应用。但是,在将稀疏表示用于二维图像处理时,考虑到计算的复杂性,大多都是将图像分后对图像块进行单独编码,当将其划分为小块时,图像的潜在结构信息可能会丢失。

发明内容

[0004] 本发明所要解决的技术问题是克服现有技术的不足,提供一种基于加权L1 范数的卷积稀疏编码方法,将加权向量引入卷积稀疏编码中,再将卷积问题转化为频域,通过Sherman-Morrison公式求解减少计算量,具有编码效果好,计算速度快的优点。
[0005] 本发明提供一种基于加权L1范数的卷积稀疏编码方法,包括如下步骤,[0006] 步骤S1、对字典过滤器{dm}和目标图像s进行傅里叶变换,并计算L1范数的自适应权重wm;
[0007] 步骤S2、根据交替方向乘子法构建卷积稀疏编码优化模型,将优化问题分解为x子问题和y子问题;
[0008] 步骤S3、根据卷积定律将x子问题转化为频域;
[0009] 步骤S4、通过变量重排列将x子问题的求解分解为若干个独立的线性问题进行求解;
[0010] 步骤S5、通过软阈值法对y子问题进行求解;
[0011] 步骤S6、根据获取的x子问题和y子问题的求解更新对偶变量um,计算出原始残差rp和对偶残差rd及计算原始停止条件∈pri和对偶停止条件∈dua,更新惩罚系数ρ,当循环达到最大迭代次数或原始残差和对偶残差满足停止条件rp≤∈pri且rd≤∈pri时,输出卷积稀疏特征图,否则调转到步骤S4。
[0012] 作为本发明的进一步技术方案,步骤S1中,对给定的m个字典过滤器{dm} 和目标图像s做傅里叶变换,获取 和 其中,目标图像s转化为一个N 维向量,N是目标图像s的像素数量;再对 和 的商做傅里叶反变换,获取 其中, 表示傅里叶反变换;则可以计算出自适应权重 ,其中,γ是给定的权重参数
[0013] 进一步的,步骤S2中的卷积稀疏编码优化模型为其中,dm为一组字典过滤器,s为目标图像,*为卷积操作,xm为一组卷积稀疏特征图,λ为正则化参数,⊙为哈达玛积, 为向量x的L2范数的平方,即向量x各元素的平方和,||x||1为向量x的L1范数。
[0014] 进一步的,步骤S3中,将x子问题转化为频域的具体方法为:首先引入辅助变量zm=ym-um,同时,使用dm的对称托普利兹矩阵与xm的矩阵乘法代替 dm和xm的卷积,即Dmxm=dm*xm,通过卷积定理可知函数卷积的傅里叶变换是函数傅里叶变换的乘积,则将dm,Dm,xm,s,zm经傅里叶变换为 再利用卷积定理和对称托普利兹矩阵的性质进行变换: 其中, 为傅里叶变换,diag为
将N维向量转化为N×N的对矩阵;变换后可知,x子问题等价于:
[0015] 其中,为由M个N×N的对角矩阵串联而成的块对角矩阵,M 为字典过滤器
的个数,N为输入图像s的维度;则可通过线性系统求解: 其中
为MN×MN的矩阵, 中非零元素在第 n列的行只会与 中非零元素在第n行的列有非零乘积。
[0016] 进一步的,步骤S4的具体方法为:在线性系统 中,右侧第m块为 则有 将向量x的第n个元素作为x(n),则有
向量an为 的第n行的所有非零元素,N个M×M的独立线性系
统为 即
通过Sherman-Morrison公式 可得
将N个独立线性体统全部求解后可得到 再利用傅里叶反变换
获取新的x。
[0017] 进一步的,步骤S5中,通过软阈值法求解为 其中,sign(·)为符号函数,|·|取绝对值,则Y子问题的解为
[0018] 进一步的,步骤S6中,根据求得的ym和xm,应用公式um=um+xm-ym更新um,原始残差rp=||x-y||2和对偶残差rd=ρ||yprev-y||2,其中yprev等于上一轮的y;原始停止条件和对偶停止条件 其中,∈abs和∈rel为给定的绝对停止容忍度和相对停止容忍度;采用自适应策略更新
其中,τ和μ为常数,常用的取值为τ=2,μ=10。
[0019] 本发明将加权向量与卷积稀疏编码算法相结合,提出了一种基于加权L1范数的卷积稀疏编码方法,不仅考虑了图像局部位置不同稀疏度的问题,计算过程中将卷积问题转化为频域,采用Sherman-Morrison公式进行求解,减少了计算量。
[0020] 与现有方案对比,本发明具有编码效果好,计算速度快,图像处理效果好的优点。附图说明
[0021] 图1为本发明的方法流程示意图;
[0022] 图2为本发明的步骤S4中变量重排列示意图;
[0023] 图3为本发明的方法与另外两种方法的图像重建对比图。

具体实施方式

[0024] 请参阅图1,本实施例提供一种基于加权L1范数的卷积稀疏编码方法,包括如下步骤,
[0025] 步骤S1、对字典过滤器{dm}和目标图像s进行傅里叶变换,并计算L1范数的自适应权重wm;
[0026] 步骤S2、根据交替方向乘子法构建卷积稀疏编码优化模型,将优化问题分解为x子问题和y子问题;
[0027] 步骤S3、根据卷积定律将x子问题转化为频域;
[0028] 步骤S4、通过变量重排列将x子问题的求解分解为若干个独立的线性问题进行求解;
[0029] 步骤S5、通过软阈值法对y子问题进行求解;
[0030] 步骤S6、根据获取的x子问题和y子问题的求解更新对偶变量um,计算出原始残差rp和对偶残差rd及计算原始停止条件∈pri和对偶停止条件∈dua,更新惩罚系数ρ,当循环达到最大迭代次数或原始残差和对偶残差满足停止条件 rp≤∈pri且rd≤∈pri时,输出卷积稀疏特征图,否则调转到步骤S4。
[0031] 步骤S1中,对给定的m个字典过滤器{dm}和目标图像s做傅里叶变换,获取 和其中,目标图像s转化为一个N维向量,N是目标图像s的像素数量;再对 和 的商做傅里叶反变换,获取 其中, 表示傅里叶反变换;则可以计算出自适应权重 ,其中,γ是给定的权重参数
[0032] 步骤S2中的卷积稀疏编码优化模型为其中,dm为一组字典过滤器,s为目标图像,*为卷积操作,xm为一组卷积稀疏特征图,λ为正则化参数,⊙为哈达玛积, 为向量x的L2范数的平方,即向量x各元素的平方和,||x||1为向量x的L1范数。
[0033] 步骤S3中,将x子问题转化为频域的具体方法为:首先引入辅助变量 zm=ym-um,同时,使用dm的对称托普利兹矩阵与xm的矩阵乘法代替dm和xm的卷积,即Dmxm=dm*xm,通过卷积定理可知函数卷积的傅里叶变换是函数傅里叶变换的乘积,则将dm,Dm,xm,s,zm经傅里叶变换为 再利用卷积定理和对称托普利兹矩阵的性质进行变换:其中,F为傅里叶变换, diag为将N维
向量转化为N×N的对角矩阵;变换后可知,x子问题等价于:
[0034] 其中,为由M个N×N的对角矩阵串联而成的块对角矩阵,M 为字典过滤器
的个数,N为输入图像s的维度;则可通过线性系统求解: 其中
为MN×MN的矩阵, 中非零元素在第 n列的行只会与 中非零元素在第n行的列有非零乘积。
[0035] 步骤S4的具体方法为:在线性系统 中,右侧第m 块为则有 将向量x的第n个元素作为x(n) ,则有
如图2所示,向量an为 的第n行的所有非零元素;N个M×M的
独立线性系统为
即 通过Sherman-Morrison公式 可得
将N个独立线性体统全部求解后可得到 再利用傅里叶反变换
获取新的x。
[0036] 步骤S5中,通过软阈值法求解为 其中, sign(·)为符号函数,|·|取绝对值,则Y子问题的解为
[0037] 步骤S6中,根据求得的ym和xm,应用公式um=um+xm-ym更新um,原始残差rp=||x-y||2和对偶残差rd=ρ||yprev-y||2,其中yprev等于上一轮的y;原始停止条件和对偶停止条件 其中,∈abs和∈rel为给定的绝对停止容忍度和相对停止容忍度;采用自适应策略更新
其中,τ和μ为常数,常用的取值为τ=2,μ=10。
[0038] 本实施例中,该方法在Lena和barbara上进行图像重建实验。依据本方法步骤将最终求得的卷积稀疏特征图x与过完备字典D做卷积运算,得到重建图像s*,|s-s*|为原始图像和重建图像的差异图像。如图3所示,将本发明的图像重建实验结果的差异图像,与原始图像和另外两种算法CBPDN和OMP进行对比,图例表示各灰度代表的差异数值。可以看出,本发明提出的方法重构图像与原图像最为相似
[0039] 以上显示和描述了本发明的基本原理、主要特征和优点。本领域的技术人员应该了解,本发明不受上述具体实施例的限制,上述具体实施例和说明书中的描述只是为了进一步说明本发明的原理,在不脱离本发明精神范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护的范围由权利要求书及其等效物界定。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈