CN102063729A - 基于二维稀疏性的压缩感知图像重建方法 - Google Patents
基于二维稀疏性的压缩感知图像重建方法 Download PDFInfo
- Publication number
- CN102063729A CN102063729A CN 201010615683 CN201010615683A CN102063729A CN 102063729 A CN102063729 A CN 102063729A CN 201010615683 CN201010615683 CN 201010615683 CN 201010615683 A CN201010615683 A CN 201010615683A CN 102063729 A CN102063729 A CN 102063729A
- Authority
- CN
- China
- Prior art keywords
- mrow
- vec
- dimensional
- formula
- compressed sensing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 239000011159 matrix material Substances 0.000 claims description 22
- 238000005259 measurement Methods 0.000 claims description 21
- 238000005457 optimization Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 abstract 1
- 238000007796 conventional method Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
基于二维稀疏性的压缩感知图像重建方法。它涉及一种压缩感知系统中的图像重建方法。它解决了现有压缩感知图像重建方法中所存在的只利用一维稀疏性的不足。步骤如下:一、对二维图像的测量值Y2执行矢量化操作vec;二、利用Kronecker积将公式五变换为如下形式:Y=ΦΨu=Θu公式十一;三、利用传统的压缩感知重建方法求解公式十一获得系数向量u;四、再对系数向量u执行逆矢量化操作ivec,获得二维图像X2的稀疏域表示系数S2:S2=ivec(u)再利用X2=Ψ2S2Ψ3重建出原始二维图像X2。本发明方法的重建图像峰值信噪比均高于传统方式的重建图像峰值信噪比,也即本发明的方法的重建质量均高于传统方式的重建图像质量。
Description
技术领域
本发明涉及一种压缩感知系统中的图像重建方法。
背景技术
压缩感知理论是一种新的信号获取、处理方面的理论。与在数据获取领域通常的经验相反,该理论认为许多自然信号和人造信号都是稀疏的,利用随机测量矩阵可把一个具有稀疏特性(或可压缩)的高维信号投影到低维的空间上,并证明了这样的随机投影包含了重建信号的足够信息,然后通过求解一个优化问题就可以从这较少的投影值中完全重建出原始信号。
现有的一维信号X(N×1维)和二维图像X2(N1×N2维)的压缩感知测量和重建过程描述如下:
根据稀疏矩阵Ψ(N×N维),得到一维信号X的稀疏域表示:
X=ΨS 公式一
其中,S是一维信号X的稀疏域表示系数(N×1维),具有K个非零系数和(N-K)个零;
对于一维信号X的压缩感知测量过程如下:
利用测量矩阵Φ(M×N维),将一维信号X投影到M(M<<N)个测量值Y(M×1维)上;
Y=ΦX=ΦΨS=ΘS 公式二
其中Θ=ΦΨ,Θ是个M×N维矩阵;
在测量过程后,先利用公式二从测量值Y中求取出一维信号X的稀疏域表示系数S,然后再利用公式一重建出原始一维信号X,在重建的过程中利用了一维信号X的稀疏性;
对于二维图像X2(N1×N2维),利用二维图像的行稀疏矩阵Ψ2(N1×N1维)和二维图像的列稀疏矩阵Ψ3(N2×N2维),得到二维图像X2的稀疏域表示:
X2=Ψ2S2Ψ3 公式三
其中,S2是二维图像X2的稀疏域表示系数(N1×N2维);
利用二维图像的测量矩阵Φ2(M×N1×N2维),获得M个二维图像的测量值Y2(M×1维),上述过程用如下公式表示:
Y2=Φ2*X2=Φ2*(Ψ2S2Ψ3) 公式四
在测量过程后,需要根据公式四从二维图像的测量值Y2中求取出二维图像X2的稀疏域表示系数S2,然后再利用公式三重建出原始二维图像X2,重建过程就利用到了原始二维图像X2的稀疏特性,如公式三所示。
然而,对于公式四的求解,由于涉及到二维矩阵的求解难度较大;另一方面,现有的压缩感知重建方法大都是针对一维信号的,即使是二维图像获取和重建也是将二维图像逐行(或者列)将二维图像看成一维信号操作的,只用到了一维行(或者列)的稀疏特性,相邻的行(或者列)间的相关性没有用到。这样的求解会影响重建图像的质量,造成重建图像质量的下降。考虑到这样的问题,我们希望在重建图像的过程中应用图像的二维稀疏性,而不是只应用一维行(或者列)的稀疏性,从而提高图像重建的质量。
发明内容
本发明为了解决现有压缩感知图像重建方法中所存在的只利用一维稀疏性的不足,提出一种基于二维稀疏性的压缩感知图像重建方法。
基于二维稀疏性的压缩感知图像重建方法的步骤如下:
步骤一:对二维图像的测量值Y2执行矢量化操作vec:
步骤二:利用Kronecker积将公式五变换为如下形式:
Y=ΦΨu=Θu 公式十一
步骤三:利用传统的压缩感知重建方法求解公式十一获得系数向量u;
步骤四:再对系数向量u执行逆矢量化操作ivec,获得二维图像X2的稀疏域表示系数S2:
S2=ivec(u) 公式十二
再利用X2=Ψ2S2Ψ3重建出原始二维图像X2。
本发明中提到的一维信号包括具有稀疏特性的语音信号、雷达信号、通信信号等,二维图像包括具有稀疏特性的遥感图像、自然图像、医学图像(CT、核磁)等。为了验证所提出的算法的性能,进行了计算机仿真实验。针对一幅遥感图像(图1)和一幅自然图像(图4),对图中取出的32×32的子图(图2和图5)进行了实验。采用随机测量矩阵测量,实验中分别采用了不同的获取率,重建使用的稀疏基是傅里叶基,使用的重建方法是典型的稀疏重建的梯度追踪算法(Gradient Pursuit for Sparse Reconstruction,GPSR)方法,对比了传统的只利用一维稀疏性的重建结果和本发明利用二维稀疏性的重建结果,分别标记为传统方法和本发明方法,重建质量用重建图像的峰值信噪比(Peak Signal to Noise Ratio,PSNR)评价。实验结果(图3和图6),从图3和图6中可见,在多种获取率下,本发明方法的重建图像峰值信噪比均高于传统方式的重建图像峰值信噪比,也即本发明的方法的重建质量均高于传统方式的重建图像质量。
附图说明
图1是一幅遥感图像;图2是图1中取出的32×32的子图;图3是本发明方法的重建图像峰值信噪比与传统方式的重建图像峰值信噪比的对比图;图4是一幅自然图像;图5是图4中取出的32×32的子图;图6是本发明方法的重建图像峰值信噪比与传统方式的重建图像峰值信噪比的对比图。
具体实施方式
具体实施方式一:本实施方式的步骤如下:
步骤一:对二维图像的测量值Y2执行矢量化操作vec:
如果考虑到压缩感知过程中测量矩阵中系数的随机性,只要二维图像的测量矩阵Φ2不是非结构性,那么对二维图像的测量值Y2执行矢量化操作vec:
公式五中,矢量化操作vec的含义如下:
设A=(ai,j)M×N,称M×N维的列向量vec(A)为矩阵A的按行展开,或矩阵A按行拉直的列向量,如下式:
vec(A)=(a11,…,a21,…,am1,…,amn)T 公式六
公式五中,是二维图像的测量矩阵Φ2的变换形式(维的矩阵,),如果令中的每一个元素为1≤i≤M,令二维图像的测量矩阵Φ2中的每一个元素为Φ2(i,j,k),1≤i≤M,1≤j≤N1,1≤k≤N2,那么是由二维图像的测量矩阵Φ2确定的,关系如下:
步骤二:利用Kronecker积将公式五变换为如下形式:
Kronecker积的具体解释如下:
根据矩阵分析中的理论,设A∈Cm×n,X∈Cn×p,B∈Cp×q,那么则有下式成立
其定义如公式十,其中设A=(ai,j)∈Cm×n,B=(bi,j)∈Cp×q。
Y=ΦΨu=Θu 公式十一
步骤三:利用传统的压缩感知重建方法求解公式十一获得系数向量u,传统的压缩感知重建方法包括凸优化算法中的包括基追踪(Basis Pursuit,BP)、稀疏重建的梯度追踪算法(Gradient Pursuit for Sparse Reconstruction,GPSR)、Bregman迭代;贪婪方法算法中的匹配追踪(Matching Pursuit,MP)、正交匹配追踪(Orthogonal Matching Pursuit,OMP)、梯度追踪(Gradient Pursuit,GP)和逐步正交匹配追踪(Stagewise Orthogonal Matching Pursuit,StOMP);
步骤四:再对系数向量u执行逆矢量化操作ivec,所述的逆矢量化操作ivec操作是矢量化操作vec的逆过程,获得二维图像X2的稀疏域表示系数S2:
S2=ivec(u) 公式十二
再利用X2=Ψ2S2Ψ3重建出原始二维图像X2。这种求解方式利用到了二维图像的稀疏特性,而不是按照传统的方法只利用其一维稀疏特性求解。
Claims (2)
2.根据权利要求1所述的基于二维稀疏性的压缩感知图像重建方法,其特征在于步骤三中的传统的压缩感知重建方法包括凸优化算法中的包括基追踪、稀疏重建的梯度追踪算法、Bregman迭代;贪婪方法算法中的匹配追踪、正交匹配追踪、梯度追踪和逐步正交匹配追踪。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010615683 CN102063729A (zh) | 2010-12-30 | 2010-12-30 | 基于二维稀疏性的压缩感知图像重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010615683 CN102063729A (zh) | 2010-12-30 | 2010-12-30 | 基于二维稀疏性的压缩感知图像重建方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102063729A true CN102063729A (zh) | 2011-05-18 |
Family
ID=43998993
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010615683 Pending CN102063729A (zh) | 2010-12-30 | 2010-12-30 | 基于二维稀疏性的压缩感知图像重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102063729A (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102438102A (zh) * | 2011-10-26 | 2012-05-02 | 西安电子科技大学 | 基于压缩编码孔径的超分辨率成像系统及成像方法 |
CN102609920A (zh) * | 2012-02-17 | 2012-07-25 | 上海交通大学 | 一种基于压缩感知的彩色数字图像修复方法 |
CN102801428A (zh) * | 2012-08-13 | 2012-11-28 | 程涛 | 一种0-1稀疏循环矩阵的近似优化和信号采集重构方法 |
CN102890275A (zh) * | 2011-07-22 | 2013-01-23 | 中国科学院电子学研究所 | 一种基于交叉熵最小化的压缩感知雷达成像方法 |
CN103237204A (zh) * | 2013-03-21 | 2013-08-07 | 上海交通大学 | 基于高维压缩感知的视频信号采集与重构系统 |
CN103584835A (zh) * | 2013-09-24 | 2014-02-19 | 南京大学 | 一种基于压缩感知的光声图像重建方法 |
CN104933685A (zh) * | 2015-06-16 | 2015-09-23 | 西安电子科技大学 | 基于三维张量压缩感知的高光谱压缩成像方法 |
CN106204501A (zh) * | 2016-07-29 | 2016-12-07 | 上海科技大学 | 一种压缩感知恢复方法 |
CN107037429A (zh) * | 2017-04-17 | 2017-08-11 | 电子科技大学 | 基于门限梯度追踪算法的线阵sar三维成像方法 |
CN107689067A (zh) * | 2017-08-30 | 2018-02-13 | 中国科学院光电技术研究所 | 一种基于压缩感知在图像重建中投影方式的优化方法 |
CN109089123A (zh) * | 2018-08-23 | 2018-12-25 | 江苏大学 | 基于1比特矢量量化的压缩感知多描述编解码方法 |
CN109509229A (zh) * | 2018-11-12 | 2019-03-22 | 凌云光技术集团有限责任公司 | 一种基于二维线性变换的模板重建装置及方法 |
CN111314000A (zh) * | 2020-02-18 | 2020-06-19 | 五邑大学 | 一种基于低秩数据驱动权重模型的射频层析成像方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009087493A1 (en) * | 2008-01-11 | 2009-07-16 | Zoran (France) | Sparse geometry for super resolution video processing |
CN101739666A (zh) * | 2009-12-15 | 2010-06-16 | 西南交通大学 | 基于一维快速哈特莱变换和匹配追踪的图像稀疏分解快速算法 |
-
2010
- 2010-12-30 CN CN 201010615683 patent/CN102063729A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009087493A1 (en) * | 2008-01-11 | 2009-07-16 | Zoran (France) | Sparse geometry for super resolution video processing |
CN101739666A (zh) * | 2009-12-15 | 2010-06-16 | 西南交通大学 | 基于一维快速哈特莱变换和匹配追踪的图像稀疏分解快速算法 |
Non-Patent Citations (5)
Title |
---|
《光电子·激光》 20071130 李林静,冯勇,刘洪臣 基于投影的单幅图像超分辨率重构算法 第1379页 1-2 第18卷, 第11期 * |
《光电子·激光》 20071130 李林静,冯勇,刘洪臣 基于投影的单幅图像超分辨率重构算法 第1379页 1-2 第18卷, 第11期 2 * |
《电子学报》 20090430 刘洪臣,冯勇,赵永平 基于Kronecker积的图像超分辨率快速算法 第840页至843页 1-2 第37卷, 第4期 2 * |
《计算机应用》 20080930 傅迎华 可压缩传感重构算法与近似QR分解 第2300页 1-2 第28卷, 第9期 * |
《计算机科学》 20101031 李小燕,尹忠科 图像1DFFT2MP 稀疏分解算法研究 第246页到第247页 1-2 第37卷, 第10期 2 * |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102890275A (zh) * | 2011-07-22 | 2013-01-23 | 中国科学院电子学研究所 | 一种基于交叉熵最小化的压缩感知雷达成像方法 |
CN102890275B (zh) * | 2011-07-22 | 2014-04-16 | 中国科学院电子学研究所 | 一种基于交叉熵最小化的压缩感知雷达成像方法 |
CN102438102A (zh) * | 2011-10-26 | 2012-05-02 | 西安电子科技大学 | 基于压缩编码孔径的超分辨率成像系统及成像方法 |
CN102438102B (zh) * | 2011-10-26 | 2014-06-18 | 西安电子科技大学 | 基于压缩编码孔径的超分辨率成像方法 |
CN102609920B (zh) * | 2012-02-17 | 2014-09-17 | 上海交通大学 | 一种基于压缩感知的彩色数字图像修复方法 |
CN102609920A (zh) * | 2012-02-17 | 2012-07-25 | 上海交通大学 | 一种基于压缩感知的彩色数字图像修复方法 |
CN102801428A (zh) * | 2012-08-13 | 2012-11-28 | 程涛 | 一种0-1稀疏循环矩阵的近似优化和信号采集重构方法 |
CN102801428B (zh) * | 2012-08-13 | 2015-02-18 | 程涛 | 一种0-1稀疏循环矩阵的近似优化和信号采集重构方法 |
CN103237204B (zh) * | 2013-03-21 | 2016-03-02 | 上海交通大学 | 基于高维压缩感知的视频信号采集与重构系统 |
CN103237204A (zh) * | 2013-03-21 | 2013-08-07 | 上海交通大学 | 基于高维压缩感知的视频信号采集与重构系统 |
CN103584835B (zh) * | 2013-09-24 | 2015-05-13 | 南京大学 | 一种基于压缩感知的光声图像重建方法 |
CN103584835A (zh) * | 2013-09-24 | 2014-02-19 | 南京大学 | 一种基于压缩感知的光声图像重建方法 |
CN104933685A (zh) * | 2015-06-16 | 2015-09-23 | 西安电子科技大学 | 基于三维张量压缩感知的高光谱压缩成像方法 |
CN106204501B (zh) * | 2016-07-29 | 2019-03-19 | 上海科技大学 | 一种压缩感知恢复方法 |
CN106204501A (zh) * | 2016-07-29 | 2016-12-07 | 上海科技大学 | 一种压缩感知恢复方法 |
CN107037429A (zh) * | 2017-04-17 | 2017-08-11 | 电子科技大学 | 基于门限梯度追踪算法的线阵sar三维成像方法 |
CN107037429B (zh) * | 2017-04-17 | 2020-06-16 | 电子科技大学 | 基于门限梯度追踪算法的线阵sar三维成像方法 |
CN107689067A (zh) * | 2017-08-30 | 2018-02-13 | 中国科学院光电技术研究所 | 一种基于压缩感知在图像重建中投影方式的优化方法 |
CN107689067B (zh) * | 2017-08-30 | 2020-06-19 | 中国科学院光电技术研究所 | 一种基于压缩感知在图像重建中投影方式的优化方法 |
CN109089123A (zh) * | 2018-08-23 | 2018-12-25 | 江苏大学 | 基于1比特矢量量化的压缩感知多描述编解码方法 |
CN109089123B (zh) * | 2018-08-23 | 2021-08-03 | 江苏大学 | 基于1比特矢量量化的压缩感知多描述编解码方法 |
CN109509229A (zh) * | 2018-11-12 | 2019-03-22 | 凌云光技术集团有限责任公司 | 一种基于二维线性变换的模板重建装置及方法 |
CN109509229B (zh) * | 2018-11-12 | 2020-12-15 | 凌云光技术集团有限责任公司 | 一种基于二维线性变换的模板重建装置及方法 |
CN111314000A (zh) * | 2020-02-18 | 2020-06-19 | 五邑大学 | 一种基于低秩数据驱动权重模型的射频层析成像方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102063729A (zh) | 基于二维稀疏性的压缩感知图像重建方法 | |
CN103400402B (zh) | 基于低秩结构稀疏的压缩感知mri图像重建方法 | |
CN103955904B (zh) | 一种基于离散分数阶傅里叶变换相位信息的信号重建方法 | |
Yang et al. | Seismic data reconstruction via matrix completion | |
CN111047515A (zh) | 一种基于注意力机制的空洞卷积神经网络图像超分辨率重建方法 | |
CN103745487B (zh) | 基于结构化稀疏先验的贝叶斯高光谱解混压缩感知方法 | |
CN107016656B (zh) | 一种基于压缩感知在图像重建中的小波稀疏基优化方法 | |
CN110596645B (zh) | 二维免求逆稀疏贝叶斯学习快速稀疏重构方法 | |
CN103279932B (zh) | 一种二维合成稀疏模型和基于该模型的字典训练方法 | |
CN112884851A (zh) | 一种展开迭代优化算法的深度压缩感知网络 | |
CN103077510B (zh) | 基于小波hmt模型的多变量压缩感知重构方法 | |
CN103279959B (zh) | 一种二维分析稀疏模型、其字典训练方法和图像去噪方法 | |
CN104779960A (zh) | 一种基于分块压缩感知的信号重构方法 | |
CN104506752B (zh) | 一种基于残差压缩感知的相似图像压缩方法 | |
CN108419083B (zh) | 一种图像多级小波全子带压缩感知编码方法 | |
CN111562545A (zh) | 基于pd-alm算法的稀疏阵列doa估计方法 | |
CN103558606A (zh) | 一种基于压缩感知的条件部分测量关联成像方法 | |
CN104915935B (zh) | 基于非线性压缩感知与字典学习的压缩光谱成像方法 | |
CN102163338B (zh) | 一种压缩感知系统中的高效重建方法 | |
CN104933685A (zh) | 基于三维张量压缩感知的高光谱压缩成像方法 | |
CN104008536A (zh) | 基于cs-chmt和idpcnn的多聚焦噪声图像融合方法 | |
CN104732566A (zh) | 基于非分离稀疏先验的高光谱图像压缩感知方法 | |
CN109887050A (zh) | 一种基于自适应字典学习的编码孔径光谱成像方法 | |
CN104185026B (zh) | 随机投影域下相位编码的红外高分辨率成像方法及其装置 | |
CN104881846A (zh) | 基于双密度双树复小波的结构化图像压缩感知还原方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20110518 |