CN103150709A - 一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法 - Google Patents

一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法 Download PDF

Info

Publication number
CN103150709A
CN103150709A CN2013100579842A CN201310057984A CN103150709A CN 103150709 A CN103150709 A CN 103150709A CN 2013100579842 A CN2013100579842 A CN 2013100579842A CN 201310057984 A CN201310057984 A CN 201310057984A CN 103150709 A CN103150709 A CN 103150709A
Authority
CN
China
Prior art keywords
matrix
quaternion
compressed sensing
image
observed reading
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.)
Granted
Application number
CN2013100579842A
Other languages
English (en)
Other versions
CN103150709B (zh
Inventor
严路
伍家松
杨淳沨
沈傲东
舒华忠
韩旭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201310057984.2A priority Critical patent/CN103150709B/zh
Publication of CN103150709A publication Critical patent/CN103150709A/zh
Application granted granted Critical
Publication of CN103150709B publication Critical patent/CN103150709B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Processing (AREA)

Abstract

本发明公开了一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法,将彩色图像二维矩阵转化到四元数域的二维矩阵,对四元数域的二维矩阵只进行一次压缩感知就能恢复原始的彩色图像,比传统的对彩色图像的RGB三个分量分别进行压缩感知要节约运算时间。本发明将四元数矩阵信号写成四元数的欧拉形式,用幅度和相位作为压缩感知优化问题新的约束项,比传统将彩色图像RGB三个通道的数据转化为三个实数二维矩阵分别处理的恢复结果更好。本发明的结果图像是按行进行压缩传感解码得到的图像与按列进行压缩传感解码得到的图像的均值,这样比单独按行处理或单独按列处理恢复的图像更为平滑。

Description

一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法
技术领域
本发明涉及一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法,属于数字图像处理技术领域。
背景技术
压缩感知(Compressed Sensing:CS)理论指出:只要信号是可压缩的或在某个变换域是稀疏的,那么就可以用一个与变换基不相关的观测矩阵将变换所得高维信号投影到一个低维空间上,然后通过求解一个优化问题就可以从这些少量的投影中以高概率重构出原始信号。其理论框架如图1所示,图中CS信息矩阵A=ΦΨ*,其中Φ为测量矩阵,Ψ为稀疏矩阵,上标“*”表示共轭转置。由图可知:
(1)CS的最终目的是用尽可能少的测量数据来恢复原始信号。
(2)CS理论框架主要分为三步:信号的稀疏表示、信号测量和信号恢复。其中稀疏矩阵Ψ可以根据信号自身的特点在正交基或字典中灵活选取。因此关于CS的研究主要集中在后两个步骤:测量矩阵Φ的设计和快速鲁棒的信号恢复算法。
对于图1中第三步,信号的恢复算法可以等价于求解如下优化问题:
x ^ = arg min J ( x ) s . t . J ( x ) = | | y - Ax | | 2 2 + λ 1 | | ψ * x | | 0 , - - - ( 1 )
上式求解是个NP-Hard问题,由于l1范数在一定条件下和l0范数具有等价性,可得到相同的解。那么上式可以转化为l1范数下的最优化问题:
x ^ = arg min J ( x ) s . t . J ( x ) = | | y - Ax | | 2 2 + λ 1 | | ψ * x | | 1 , - - - ( 2 )
l1范数下最优化问题又称为基追踪(Basis Pursuit:BP),其常用实现算法有:伪牛顿法、内点法和梯度投影法。
四元数q是一种特殊的超复数,是复数的扩展形式。四元数q由1个实部和3个虚部组成:
q=S(q)+X(q)i+Y(q)j+Z(q)k,    (3)
其中S(q),X(q),Y(q),
Figure BDA00002852066500015
表示实数域,i,j,k是三个虚数单位,满足如下性质:
i2=j2=k2=ijk=-1,    (4)
ij=-ji=k,jk=-kj=i,ki=-ik=j.    (5)
四元数的共轭和模分别定义为:
q ‾ = S ( q ) - X ( q ) i - Y ( q ) j - Z ( q ) k , - - - ( 6 )
| q | = S 2 ( q ) + X 2 ( q ) + Y 2 ( q ) + Z 2 ( q ) . - - - ( 7 )
四元数q也可以用欧拉形式表达为幅度和相位的形式:
q = S ( q ) + X ( q ) i + Y ( q ) j + Z ( q ) k = | q | e η θ ,
| q | = S 2 ( q ) + X 2 ( q ) + Y 2 ( q ) + Z 2 ( q )
η = X ( q ) i + Y ( q ) j + Z ( q ) k X 2 ( q ) + Y 2 ( q ) + Z 2 ( q ) , - - - ( 8 )
θ=arccos(S(q)/|q|).
η表示归一化的纯四元数,θ表示四元数相位角。
类似地,(3)和(8)的定义形式也适用于四元数矩阵x。
四元数矩阵x的lp范数为:
| | x | | p = ( &Sigma; r = 1 m | x r | p ) 1 / p 0 < p < 2 - - - ( 9 )
目前压缩感知对彩色图像的处理均采用传统的RGB三个分量单独处理的方式,而没有考虑三个通道数据之间的相关性,使得恢复算法的约束信息不够完善,恢复出来的图像效果没有达到理想的效果。
发明内容
发明目的:针对传统彩色图像压缩感知恢复算法中的将RGB三个通道数据之间的相关性割裂的问题与不足,本发明结合压缩感知理论和四元数信号处理理论,提供一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法。将彩色图像转化为四元数矩阵信号,采用四元数矩阵信号的欧拉表示形式,在压缩感知恢复算法中加入了四元数矩阵信号的幅度和相位作为新的约束项,使恢复的图像更为平滑。
为了将彩色图像三个通道数据之间的相关性相结合,将彩色图像的RGB三个分量分别放在四元数的一个实部和两个虚部,这样将彩色图像中的信息转化为四元数矩阵。同时我们将四元数信号表示成欧拉形式,即幅度和相位的形式,这样我们就保留了三个通道数据之间的相关性,从而能更好的恢复彩色图像。
技术方案:一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法,包括编码端对彩色图像进行稀疏表示及测量编码,解码端对接收的数据进行解码重构,得到恢复图像。具体包括以下步骤:
步骤A、将大小为N×N的彩色图像T的RGB三个通道的数据分别赋值给N×N的四元数矩阵Q的一个实部和两个虚部,并将第三个虚部置零,即
Q=R(T)+G(T)i+B(T)j+0k;i,j,k是三个虚数单位,R(T),G(T)和B(T)分别表示彩色图像R分量、G分量和B分量的数据值;
步骤B、用大小为N×N的稀疏矩阵W对四元数矩阵Q进行稀疏,得到稀疏后的四元信号矩阵Q0,稀疏矩阵W为离散小波变换(DWT)矩阵;
步骤C、用大小为M×N的观测矩阵Φ对稀疏后的四元信号矩阵Q0进行观测,得到M×N的观测值Y1,观测矩阵Φ为随机高斯矩阵;(注:对四元信号矩阵Q0进行观测相当于按每列进行观测处理)其中,M<<N;
步骤D、用观测矩阵Φ对稀疏后的四元信号矩阵Q0的转置进行观测,得到M×N的观测值Y2;(注:对四元信号矩阵Q0的转置进行观测相当于对四元信号矩阵Q0按每行进行观测处理)
步骤E、设置压缩感知优化方程迭代的初始值Qt,Qt=Φ*×Y(这里的Y分别指Y1和Y2),上标“*”表示共轭转置,并且将得到的四元数信号矩阵Qt改写为幅度和相位的形式:η表示归一化的纯四元数,θt表示四元数相位角;
Figure BDA00002852066500032
表示两个矩阵的点乘积,即矩阵对应位置的元素相乘;
步骤F、用观测值Y1和Y2,以及稀疏矩阵W和观测矩阵Φ,结合四元数信号压缩感知恢复算法,从观测值中高概率重构出四元数信号
Figure BDA00002852066500034
其中是通过观测值Y1重构得到的信号,
Figure BDA00002852066500036
是通过观测值Y2重构得到的信号;
步骤G、将得到的四元数信号
Figure BDA00002852066500037
Figure BDA00002852066500038
相加求平均,得到四元数矩阵
Figure BDA00002852066500039
Figure BDA000028520665000310
的实部赋值给彩色图像的R通道,将的两个虚部分别赋值给彩色图像的G通道和B通道,得到恢复的彩色图像。
将彩色图像T的RGB三个通道的数据分别赋值给大小为N×N四元数矩阵Q的一个实部和两个虚部过程中,两个虚部的选择可以是任意的,但一定要赋值给实部,保证相位有意义。
所述从观测值中高概率重构出四元数信号
Figure BDA00002852066500042
的步骤为:
步骤1)、将压缩感知恢复算法中经典的l1范数最优化问题
x ^ = arg min J ( x ) s . t . J ( x ) = | | y - Ax | | 2 2 + &lambda; 1 | | &psi; * x | | 1
改写为lp(其中0<p<2)范数下的最优化问题,同时用幅度和相位作为压缩感知优化问题新的约束项,得到目标函数J(Q);
Q ^ = arg min J ( Q )
s . t . J ( Q ) = | | Y - &Phi; WQW ' | 2 2 + &lambda; 1 &Sigma; i = 1 N ( ( WQW ' ) i 2 + &epsiv; ) p / 2 + &lambda; 2 &Sigma; i = 1 N ( ( D | Q | ) i 2 + &epsiv; ) p / 2 + &lambda; 3 &Sigma; i = 1 N ( ( D&theta; ) i 2 + &epsiv; ) p / 2
其中ε是一个很小的常数,常取1.0×10-6,D代表梯度算子,W′表示W的转置;
步骤2)、对目标函数J(Q)关于四元数信号Q求导,得到目标函数的梯度ΔJ(Q),
通过测量梯度的变化构造一个近似的海森矩阵H(Q);
步骤3)、设置迭代的步长因子γ,γ取值为0.7;
步骤4)、用梯度ΔJ(Qt),海森矩阵H(Qt)及步长因子γ,通过公式Q1=Qt-γ[H(Qt)]-1ΔJ(Qt)求得经过第一次迭代后的值Q1
步骤5)、将Qt更新成Q1,重新计算梯度ΔJ(Qt)和海森矩阵H(Qt),重复步骤I和J,直到迭代过程收敛;
步骤6)、如果观测值为Y1时,收敛后的Q1不再赋给Qt,把Q1赋给
Figure BDA00002852066500047
作为Y1的处理结果;
如果观测值为Y2时,收敛后的Q1不再赋给Qt;把Q1赋给
Figure BDA00002852066500048
作为Y2的处理结果。
有益效果:本发明的彩色图像压缩感知恢复算法,在编码时,将彩色图像的RGB三个分量分别赋值给四元数信号的一个实部和两个虚部,将彩色图像的压缩感知恢复问题转换到四元数域的压缩感知恢复问题,对稀疏后四元数矩阵信号Q0按行进行观测得到观测值Y1,对Q0按列进行观测得到观测值Y2;进行解码时,通过观测值Y1和Y2,以及稀疏矩阵W和观测矩阵Φ,得到解码重构恢复的四元数矩阵
Figure BDA00002852066500051
Figure BDA00002852066500052
将四元数矩阵
Figure BDA00002852066500053
Figure BDA00002852066500054
相加求平均,得到结果四元数矩阵
Figure BDA00002852066500055
将四元数矩阵
Figure BDA00002852066500056
的实部及两个虚部,分别赋值给彩色图像的RGB三个通道得到结果图像。本发明的方法与现有的将彩色图像RGB三个分量数据单独处理的方法相比,只需进行一次压缩感知就能恢复彩色图像,比传统的进行三次压缩感知更节约时间。同时本发明利用了四元数矩阵的幅度和相位信息,结合了三个通道之间的相关性,使得彩色图像的恢复更为精确,这是传统的彩色图像处理技术难以企及的。
附图说明
图1为压缩感知的理论框架图;
图2为本发明方法与现有方法的比较流程图;
图3为采样率从60%到85%时,按列进行处理时lena图像的恢复效果比较图;
图4为采样率从60%到85%时,按行进行处理时lena图像的恢复效果比较图;
图5为采样率从60%到85%时,按行处理与按列处理相结合时lena图像的恢复效果比较图;
图6为采样率从60%到85%时,单独按行处理、单独按列处理、按行处理与按列处理相结合的方法所对应的图像的不同恢复效果比较图。
具体实施方式
下面结合具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
如图2所示,基于伪牛顿法的四元数域彩色图像压缩感知恢复方法,按照前述步骤进行:
步骤1、将彩色图像表示成四元数矩阵的形式,一幅64×64彩色lena图像可以用一个64×64的四元数矩阵来存储其三个通道的值。
例如,彩色图像右上角3×3的像素矩阵是:
( 151,111,108 ) ( 84,37,34 ) ( 105,65,66 ) ( 130,102,98 ) ( 59,28,24 ) ( 80,62,57 ) ( 93,99 , 161 ) ( 20,29,65 ) ( 55,62,78 )
则相对应的四元数矩阵表示为:
( 151 + 111 i + 108 j + 0 k ) ( 84 + 37 i + 34 j + 0 k ) ( 105 + 65 i + 66 j + 0 k ) ( 130 + 102 i + 98 j + 0 k ) ( 59 + 28 i + 24 j + 0 k ) ( 80 + 62 i + 57 j + 0 k ) ( 93 + 99 i + 161 j + 0 k ) ( 20 + 29 i + 65 j + 0 k ) ( 55 + 62 i + 78 j + 0 k )
因此,我们可以得到大小为64×64的四元数矩阵Q。
步骤2、将四元数矩阵Q在冗余字典下进行稀疏分解,实验中我们取的是大小为64×64离散小波变换(DWT)矩阵W,根据以下公式进行稀疏分解:
Q0=W×Q×W′,W′表示矩阵W的转置.    (10)
例如,彩色图像右上角3×3的像素矩阵稀疏后得到:
( 4 - 2 i - j + 0 k ) ( 6 - i - 6 j + 0 k ) ( 4 - i - 3 j + 0 k ) ( - 2 + 0 i - j + 0 k ) ( - 6 + i + 0 j + 0 k ) ( - 3 - i + 2 j + 0 k ) ( - 2 - i - j + 0 k ) ( 5 - i + 0 j + 0 k ) ( 2 - i - j + 0 k )
对得到的3×3四元数矩阵中四元数幅度的值接近0的元素,我们将其置零,得到如下更为稀疏的四元数矩阵:
( 4 - 2 i - j + 0 k ) ( 6 - i - 6 j + 0 k ) ( 4 - i - 3 j + 0 k ) ( 0 + 0 i + 0 j + 0 k ) ( - 6 + i + 0 j + 0 k ) ( 0 + 0 i + 0 j + 0 k ) ( 0 + 0 i + 0 j + 0 k ) ( 5 - i - 0 j + 0 k ) ( 0 + 0 i + 0 j + 0 k )
因此,我们可以得到大小为64×64的稀疏四元数矩阵Q0,并且Q0中的非零元素远远少于N×N。
步骤3、用大小为M×N的随机高斯矩阵Φ对稀疏后的四元数信号Q0的按行进行观测,得到M×N的观测值Y1,接着对Q0按列进行观测,得到M×N的观测值Y2
步骤401、对原始四元数矩阵Q进行压缩感知编码以后,我们需要通过解决如下优化问题来解码恢复彩色图像:
Q ^ = arg min J ( Q )
s . t . J ( Q ) = | | Y - &Phi; WQW ' | | 2 2 + &lambda; 1 | | WQW ' | | p p + &lambda; 2 | | D | Q | | p p + &lambda; 3 | | D&theta; | | p p - - - ( 11 )
其中D代表梯度算子,形式如下:
D = - 1 1 0 0 . . 0 0 - 1 1 0 . . 0 . . . . . . . 0 0 0 0 0 - 1 1 0 0 0 0 0 0 0 - - - ( 12 )
步骤402、为了避免lp范数在初始点的不可微,我们将lp范数改写为如下形式:
| | f | | p p &ap; &Sigma; i = 1 N ( | ( f ) i | 2 + &xi; ) p / 2 - - - ( 13 )
其中ξ是一个很小的常数,常取1.0×10-6。这样(11)式就可以改写为:
Q ^ = arg min J ( Q )
s . t . J ( Q ) = | | Y - &Phi; WQW ' | 2 2 + &lambda; 1 &Sigma; i = 1 N ( ( WQW ' ) i 2 + &epsiv; ) p / 2 + &lambda; 2 &Sigma; i = 1 N ( ( D | Q | ) i 2 + &epsiv; ) p / 2 + &lambda; 3 &Sigma; i = 1 N ( ( D&theta; ) i 2 + &epsiv; ) p / 2 - - - ( 14 )
其中ε是一个很小的常数,常取1.0×10-6,D代表梯度算子,W′表示W的转置;
步骤403、对目标函数J(Q)求导,得到目标函数的梯度ΔJ(Q),通过测量梯度的变化构造一个近似的海森矩阵H(Q)。
步骤404、设置迭代的初始值Qt=Φ*×Y(这里的Y分别指Y1和Y2),通过如下公式计算Qt的幅度|Qt|和相位θt
Q = S ( Q ) + X ( Q ) i + Y ( Q ) j + 0 k = | Q | e &eta; &theta;
| Q | = S 2 ( Q ) + X 2 ( Q ) + Y 2 ( Q )
&eta; = I ( Q ) i + J ( Q ) j I 2 ( Q ) + J 2 ( Q ) - - - ( 15 )
θ=arccos(S(Q)/|Q|)
例如将得到的64×64的四元数矩阵Qt的右上角3×3的元素改写为幅度和相位的形式:
| Q t | = 202.1930 133.8469 131.0572 163.6490 122.6051 115.8836 109.8818 120.3578 190.3418
&theta; t = 0.7275 0.5929 0.6022 0.6528 0.5882 0.5630 0.5617 0.6049 0.5626
步骤405、设置迭代的步长因子γ,γ取值常设为(0,1]之间,本发明γ取值为0.7。
步骤406、结合梯度ΔJ(Q)和海森矩阵H(Q)用如下伪牛顿迭代式求得经过第一次迭代的值Q1
Q1=Qt—γ[H(Qt)]-1ΔJ(Qt)    (16)
步骤407、将Qt更新成Q1,重新计算梯度ΔJ(Q)和海森矩阵H(Q),重复步骤406到步骤407,直到迭代过程收敛。
步骤5、将通过观测值Y1和Y2重构得到的四元数矩阵
Figure BDA00002852066500084
相加求平均,得到四元数矩阵
Figure BDA00002852066500085
将四元数矩阵
Figure BDA00002852066500086
的实部及两个虚部,分别赋值给彩色图像的RGB三个通道得到结果图像T。
为了验证本发明方法的效果,进行了以下实验:
1、实验条件:
在一台计算机上进行验证实验,该计算机配置为AMD单核处理器(3200兆赫兹)和4096兆字节随机存取存储器(RAM),编程语言用的是Matlab(8.0版本)。
2、实验方法:
使用大小为64×64的彩色图像lena作为实验对象,采用以下步骤完成lena图像的压缩感知恢复。
1)采用具体实施方式中的步骤1,将彩色图像二维矩阵转化为四元数二维矩阵Q;
2)采用具体实施方式中的步骤2,对四元数矩阵Q进行稀疏,得到稀疏后的四元数矩阵Q0
3)设置不同的采样率,对四元数矩阵Q0的按行进行观测,得到观测值Y1。将公式(14)中的三个约束项的权重设置为
λ1=0.01,λ2=0.0114,λ3=0.0028,p=0.9,迭代的步长因子r=0.7,用具体实施方式中的步骤401到步骤406,得到恢复矩阵
Figure BDA00002852066500091
4)设置不同的采样率,对四元数矩阵Q0的按列进行观测,得到观测值Y2。将公式(14)中的三个约束项的权重设置为
λ1=0.005,λ2=0.0124,λ3=0.0048,p=0.9,迭代的步长因子r=0.7,用具体实施方式中的步骤401到步骤406,得到恢复矩阵
Figure BDA00002852066500092
5)采用具体实施方式中的步骤5,得到结果图像T。
3、实验结果的评价指标:
实验中采用如下公式计算彩色图像的峰值信噪比:
两个M×N的彩色图像I和K,它们RGB各个通道的均方差定义为:
MSE r = 1 MN &Sigma; i = 0 M - 1 &Sigma; j = 0 N - 1 | | I ( i , j , 1 ) - K ( i , j , 1 ) | | 2
MSE g = 1 MN &Sigma; i = 0 M - 1 &Sigma; j = 0 N - 1 | | ( i , j , 2 ) - K ( i , j , 2 ) | | 2 - - - ( 17 )
MSE b = 1 MN &Sigma; i = 0 M - 1 &Sigma; j = 0 N - 1 | | ( i , j , 3 ) - K ( i , j , 3 ) | | 2
RGB各个通道的峰值信噪比定义为:
PSNR r = 10 &CenterDot; log 10 ( MAX r 2 MSE r )
PSNR g = 10 &CenterDot; log 10 ( MAX g 2 MSE g ) - - - ( 18 )
PSNR b = 10 &CenterDot; log 10 ( MAX b 2 MSE b )
其中,MAXr,MAXg,MAXb表示图像I和K的RGB各个通道图像点颜色的最大值,如果每个采样点用8位表示,那么它们就是255。
彩色图像I和K的峰值信噪比定义为:
PSNR=(PSNRr+PSNRg+PSNRb)/3    (19)
图3显示了采样率从60%到85%时,按列进行处理时lena图像的恢复效果。表1为采样率与峰值信噪比。
表1采样率与峰值信噪比
采样率 60% 65% 70% 75% 80% 85%
信噪比 18.9260 21.4313 21.7546 22.0619 25.4913 28.6617
单独按列进行处理时,原始图像的行和行之间的相关性被割裂了。从图3的恢复结果可以看到,恢复图像的列与列之间存在明显的裂痕。当采样率提高时,这种裂痕的数量变少了,说明在提高采样率时,每列的恢复效果更为精确,可以近似地保持行与行之间相关性,但裂痕仍然会存在。
从表1的结果可以看到,当采样率提高时,恢复图像与原始图像的峰值信噪比提高了,这与压缩感知的理论知识是一致的。
图4显示了采样率从60%到85%时,按行处理时lena图像的恢复效果。表2为采样率与峰值信噪比。
表2采样率与峰值信噪比
采样率 60% 65% 70% 75% 80% 85%
信噪比 24.0742 24.5024 25.8496 26.0599 29.3833 33.3818
单独按行进行处理时,原始图像的列和列之间的相关性被割裂了,从图5的恢复结果可以看到,恢复图像的行与行之间存在明显的裂痕。同单独按行处理一样,当采样率提高时,这种裂痕的数量变少了,说明在提高采样率时,每行的恢复效果更为精确,可以近似地保持列与列之间相关性,但裂痕仍然会存在。
从表2的结果可以看到,当采样率提高时,恢复图像与原始图像的峰值信噪比提高了,这与压缩感知的理论知识是一致的。对比表1和表2我们发现,单独按行进行处理比单独按列进行处理得到的恢复图像与原始图像的峰值信噪比要高,这说明按行进行处理时,每行的稀疏度比每列的稀疏度要大,这一点与压缩感知的理论知识也是一致的。
图5显示了采样率从60%到85%时,按行处理与按列处理相结合图像的恢复效果。表3为采样率与峰值信噪比。
使用按行处理与按列处理相结合的方法,从图5的恢复结果可以看到,恢复图像行与行之间的裂痕、恢复图像列与列之间的裂痕得到了明显的改善,在采样率提高时,恢复图像更为平滑,取得了较理想的结果。
从表3的结果可以看到,当采样率提高时,恢复图像与原始图像的峰值信噪比提高了,这与压缩感知的理论知识是一致的。对比表2和表3我们发现,当采用按行处理与按列处理相结合这种方法时,恢复图像与原始图像的峰值信噪比得到了小幅度的改善,这与我们预计的结果是一致的。
为了更好的对比不同方法的恢复结果,图6显示了采样率从60%到85%时,单独按行处理、单独按列处理、按行处理与按列处理相结合的方法所对应的图像的不同恢复效果。

Claims (3)

1.一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法,其特征在于,包括以下步骤:
步骤A、将大小为N×N的彩色图像T的RGB三个通道的数据分别赋值给N×N的四元数矩阵Q的一个实部和两个虚部,并将第三个虚部置零,即Q=R(T)+G(T)i+B(T)j+0k;i,j,k是三个虚数单位,R(T),G(T)和B(T)分别表示彩色图像R分量、G分量和B分量的数据值;
步骤B、用大小为N×N的稀疏矩阵W对四元数矩阵Q进行稀疏,得到稀疏后的四元信号矩阵Q0,稀疏矩阵W为离散小波变换矩阵;
步骤C、用大小为M×N的观测矩阵Φ对稀疏后的四元信号矩阵Q0进行观测,得到M×N的观测值Y1,观测矩阵Φ为随机高斯矩阵;其中,M<<N;
步骤D、用观测矩阵Φ对稀疏后的四元信号矩阵Q0的转置进行观测,得到M×N的观测值Y2
步骤E、设置压缩感知优化方程迭代的初始值Qt,Qt=Φ*×Y,Φ*表示Φ的共轭转置,Y为观测值Y1或Y2。并且将得到的四元数信号矩阵Qt改写为幅度和相位的形式:
Figure FDA00002852066400011
η表示归一化的纯四元数,θt表示四元数相位角;ηθt表示两个矩阵的点乘积,即矩阵对应位置的元素相乘;
步骤F、用观测值Y1和Y2,以及稀疏矩阵W和观测矩阵Φ,结合四元数信号压缩感知恢复算法,从观测值中高概率重构出四元数信号
Figure FDA00002852066400012
其中
Figure FDA00002852066400014
是通过观测值Y1重构得到的信号,是通过观测值Y2重构得到的信号;
步骤G、将得到的四元数信号
Figure FDA00002852066400016
Figure FDA00002852066400017
相加求平均,得到四元数矩阵
Figure FDA00002852066400018
Figure FDA00002852066400019
的实部赋值给彩色图像的R通道,将
Figure FDA000028520664000110
的两个虚部分别赋值给彩色图像的G通道和B通道,得到恢复的彩色图像。
2.如权利要求1所述的基于伪牛顿法的四元数域彩色图像压缩感知恢复方法,其特征在于:将彩色图像T的RGB三个通道的数据分别赋值给大小为N×N四元数矩阵Q的一个实部和两个虚部过程中,两个虚部的选择可以是任意的,但要赋值给实部。
3.如权利要求1所述的基于伪牛顿法的四元数域彩色图像压缩感知恢复方法,其特征在于,所述从观测值中高概率重构出四元数信号
Figure FDA00002852066400021
Figure FDA00002852066400022
的步骤为:
步骤1)、将压缩感知恢复算法中经典的l1范数最优化问题
x ^ = arg min J ( x ) s . t . J ( x ) = | | y - Ax | | 2 2 + &lambda; 1 | | &psi; * x | | 1
改写为lp范数下的最优化问题,其中0<P<2,同时用幅度和相位作为压缩感知优化问题新的约束项,得到目标函数J(Q);
Q ^ = arg min J ( Q )
s . t . J ( Q ) = | | Y - &Phi;WQ W &prime; | | 2 2 + &lambda; 1 &Sigma; i = 1 N ( ( WQW &prime; ) i 2 + &epsiv; ) p / 2 + &lambda; 2 &Sigma; i = 1 N ( ( D | Q | ) i 2 + &epsiv; ) p / 2 + &lambda; 3 &Sigma; i = 1 N ( ( D&theta; ) i 2 + &epsiv; ) p / 2
其中ε是一个很小的常数,常取1.0×10-6,D代表梯度算子,W′表示W的转置;
步骤2)、对目标函数J(Q)关于四元数信号Q求导,得到目标函数的梯度ΔJ(Q),通过测量梯度的变化构造一个近似的海森矩阵H(Q);
步骤3)、设置迭代的步长因子γ,γ取值为0.7;
步骤4)、用梯度ΔJ(Qt),海森矩阵H(Qt)及步长因子γ,通过公式Q1=Qt-γ[H(Qt)]-1ΔJ(Qt)求得经过第一次迭代后的值Q1
步骤5)、将Qt更新成Q1,重新计算梯度ΔJ(Qt)和海森矩阵H(Qt),重复步骤4)和5),直到迭代过程收敛;
步骤6)、如果观测值为Y1时,收敛后的Q1不再赋给Qt,把Q1赋给
Figure FDA00002852066400027
作为Y1的处理结果;
如果观测值为Y2时,收敛后的Q1不再赋给Qt;把Q1赋给
Figure FDA00002852066400028
作为Y2的处理结果。
CN201310057984.2A 2013-02-22 2013-02-22 一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法 Expired - Fee Related CN103150709B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310057984.2A CN103150709B (zh) 2013-02-22 2013-02-22 一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310057984.2A CN103150709B (zh) 2013-02-22 2013-02-22 一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法

Publications (2)

Publication Number Publication Date
CN103150709A true CN103150709A (zh) 2013-06-12
CN103150709B CN103150709B (zh) 2016-10-26

Family

ID=48548763

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310057984.2A Expired - Fee Related CN103150709B (zh) 2013-02-22 2013-02-22 一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法

Country Status (1)

Country Link
CN (1) CN103150709B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104200421A (zh) * 2014-08-01 2014-12-10 南京信息工程大学 基于四元数正交变换的彩色图像加密方法及解密方法
CN106952315A (zh) * 2017-03-22 2017-07-14 广东工业大学 一种基于bfgs的对太赫兹复值数据进行图像快速重构的方法
CN109409383A (zh) * 2018-09-26 2019-03-01 重庆邮电大学 基于扩展四元局部二值模式的彩色图像特征提取方法
CN110430270A (zh) * 2019-08-08 2019-11-08 网易(杭州)网络有限公司 载具数据的同步方法及装置
CN111292284A (zh) * 2020-02-04 2020-06-16 淮阴师范学院 基于双树-四元数小波变换的彩色图像融合方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090274356A1 (en) * 2008-05-02 2009-11-05 The Uwm Research Foundation, Inc. Toeplitz random encoding for reduced acquisition in compressed sensing magnetic resonance imaging
CN101908889A (zh) * 2010-07-30 2010-12-08 哈尔滨工业大学 一种块稀疏度未知的稀疏信号压缩感知重构方法
CN102135618A (zh) * 2010-12-23 2011-07-27 北京航空航天大学 一种针对空域稀疏sar图像的频域压缩感知方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090274356A1 (en) * 2008-05-02 2009-11-05 The Uwm Research Foundation, Inc. Toeplitz random encoding for reduced acquisition in compressed sensing magnetic resonance imaging
CN101908889A (zh) * 2010-07-30 2010-12-08 哈尔滨工业大学 一种块稀疏度未知的稀疏信号压缩感知重构方法
CN102135618A (zh) * 2010-12-23 2011-07-27 北京航空航天大学 一种针对空域稀疏sar图像的频域压缩感知方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
EMMANUEL J.CANDÈS ET AL.: "Robust Uncertainty Principles: Exact Signal Reconstruction From Highly Incomplete Frequency Information", 《IEEE TRANSACTIONS ON INFORMATION THEORY》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104200421A (zh) * 2014-08-01 2014-12-10 南京信息工程大学 基于四元数正交变换的彩色图像加密方法及解密方法
CN104200421B (zh) * 2014-08-01 2017-07-28 南京信息工程大学 基于四元数正交变换的彩色图像加密方法及解密方法
CN106952315A (zh) * 2017-03-22 2017-07-14 广东工业大学 一种基于bfgs的对太赫兹复值数据进行图像快速重构的方法
CN109409383A (zh) * 2018-09-26 2019-03-01 重庆邮电大学 基于扩展四元局部二值模式的彩色图像特征提取方法
CN110430270A (zh) * 2019-08-08 2019-11-08 网易(杭州)网络有限公司 载具数据的同步方法及装置
CN110430270B (zh) * 2019-08-08 2022-03-25 网易(杭州)网络有限公司 载具数据的同步方法及装置
CN111292284A (zh) * 2020-02-04 2020-06-16 淮阴师范学院 基于双树-四元数小波变换的彩色图像融合方法
CN111292284B (zh) * 2020-02-04 2024-03-01 淮阴师范学院 基于双树-四元数小波变换的彩色图像融合方法

Also Published As

Publication number Publication date
CN103150709B (zh) 2016-10-26

Similar Documents

Publication Publication Date Title
Brunton et al. Compressed sensing and dynamic mode decomposition
CN103150709A (zh) 一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法
CN102708576A (zh) 基于结构字典的分块图像压缩感知重建方法
CN105791189B (zh) 一种提高重构精度的稀疏系数分解方法
Hong et al. An efficient algorithm for designing projection matrix in compressive sensing based on alternating optimization
Balasubramanee et al. Twitter bootstrap and AngularJS: Frontend frameworks to expedite science gateway development
CN104199627B (zh) 基于多尺度在线字典学习的可分级视频编码系统
Hartley et al. Reconstruction from projections using Grassmann tensors
CN102075749B (zh) 压缩感知框架下基于非凸模型的图像压缩重构方法
CN104732535A (zh) 一种约束稀疏的非负矩阵分解方法
CN104933685A (zh) 基于三维张量压缩感知的高光谱压缩成像方法
Li et al. On joint optimization of sensing matrix and sparsifying dictionary for robust compressed sensing systems
Dhumal et al. Survey on comparative analysis of various image compression algorithms with singular value decomposition
CN103957011A (zh) 基于阈值收缩迭代的压缩感知含噪信号的恢复方法
CN103955956B (zh) 一种面向压缩感知的图像联合重构方法
CN108288295A (zh) 基于结构信息的红外小目标图像的快速重构方法及系统
Wang et al. Image reconstruction from patch compressive sensing measurements
CN104243986B (zh) 基于数据驱动张量子空间的压缩视频采集与重构系统
CN104867167B (zh) 一种基于压缩感知的影像两步重构法
CN107818325A (zh) 基于集成字典学习的图像稀疏表示方法
Thepade et al. Varying proportions of constituent transforms to generate hybrid wavelet transform for image compression
CN104952046B (zh) 基于统计的图像压缩感知低复杂度还原方法
Elamparuthi Implementation of a hybrid color image compression technique using principal component analysis and discrete tchebichef transform
Narayanankutty et al. Understanding theory behind compressed sensing
Li et al. Mobile 3D visual search using the Helmert transformation of stereo features

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20161026

Termination date: 20190222

CF01 Termination of patent right due to non-payment of annual fee