[0002] 压缩感知(Compressed Sensing:CS)理论指出:只要
信号是可压缩的或在某个变换域是稀疏的,那么就可以用一个与变换基不相关的
观测矩阵将变换所得高维信号投影到一个低维空间上,然后通过求解一个优化问题就可以从这些少量的投影中以高概率重构出原始信号。其理论
框架如图1所示,图中CS信息矩阵A=ΦΨ*,其中Φ为测量矩阵,Ψ为稀疏矩阵,上标“*”表示共轭转置。由图可知:
[0014] ij=-ji=k,jk=-kj=i,ki=-ik=j. (5)
[0022] θ=arccos(S(q)/|q|).
[0027] 目前压缩感知对彩色图像的处理均采用传统的RGB三个分量单独处理的方式,而没有考虑三个通道数据之间的相关性,使得恢复算法的约束信息不够完善,恢复出来的图像效果没有达到理想的效果。
[0029] 为了将彩色图像三个通道数据之间的相关性相结合,将彩色图像的RGB三个分量分别放在四元数的一个实部和两个虚部,这样将彩色图像中的信息转化为四元数矩阵。同时我们将四元数信号表示成欧拉形式,即幅度和相位的形式,这样我们就保留了三个通道数据之间的相关性,从而能更好的恢复彩色图像。
[0030] 技术方案:一种基于伪牛顿法的四元数域彩色图像压缩感知恢复方法,包括编码端对彩色图像进行稀疏表示及测量编码,解码端对接收的数据进行解码重构,得到恢复图像。具体包括以下步骤:
[0032] Q=R(T)+G(T)i+B(T)j+0k;i,j,k是三个虚数单位,R(T),G(T)和B(T)分别表示彩色图像R分量、G分量和B分量的数据值;
[0034] 步骤C、用大小为M×N的观测矩阵Φ对稀疏后的四元信号矩阵Q0进行观测,得到M×N的观测值Y1,观测矩阵Φ为随机高斯矩阵;(注:对四元信号矩阵Q0进行观测相当于按每列进行观测处理)其中,M<<N;
[0035] 步骤D、用观测矩阵Φ对稀疏后的四元信号矩阵Q0的转置进行观测,得到M×N的观测值Y2;(注:对四元信号矩阵Q0的转置进行观测相当于对四元信号矩阵Q0按每行进行观测处理)
[0037] 步骤F、用观测值Y1和Y2,以及稀疏矩阵W和观测矩阵Φ,结合四元数信号压缩感知恢复算法,从观测值中高概率重构出四元数信号 和 其中 是通过观测值Y1重构得到的信号, 是通过观测值Y2重构得到的信号;
[0038] 步骤G、将得到的四元数信号 和 相加求平均,得到四元数矩阵 将 的实部赋值给彩色图像的R通道,将 的两个虚部分别赋值给彩色图像的G通道和B通道,得到恢复的彩色图像。
[0039] 将彩色图像T的RGB三个通道的数据分别赋值给大小为N×N四元数矩阵Q的一个实部和两个虚部过程中,两个虚部的选择可以是任意的,但一定要赋值给实部,保证相位有意义。
[0043] 改写为lp(其中0
[0044]
[0045]-6
[0046] 其中ε是一个很小的常数,常取1.0×10 ,D代表梯度算子,W′表示W的转置;
[0047] 步骤2)、对目标函数J(Q)关于四元数信号Q求导,得到目标函数的梯度ΔJ(Q),[0048] 通过测量梯度的变化构造一个近似的海森矩阵H(Q);
[0049] 步骤3)、设置迭代的步长因子γ,γ取值为0.7;
[0050] 步骤4)、用梯 度ΔJ(Qt),海森 矩阵H(Qt)及步 长因 子γ,通 过公 式-1Q1=Qt-γ[H(Qt)] ΔJ(Qt)求得经过第一次迭代后的值Q1;
[0051] 步骤5)、将Qt更新成Q1,重新计算梯度ΔJ(Qt)和海森矩阵H(Qt),重复步骤I和J,直到迭代过程收敛;
[0052] 步骤6)、如果观测值为Y1时,收敛后的Q1不再赋给Qt,把Q1赋给 作为Y1的处理结果;
[0053] 如果观测值为Y2时,收敛后的Q1不再赋给Qt;把Q1赋给 作为Y2的处理结果。
[0054] 有益效果:本发明的彩色图像压缩感知恢复算法,在编码时,将彩色图像的RGB三个分量分别赋值给四元数信号的一个实部和两个虚部,将彩色图像的压缩感知恢复问题转换到四元数域的压缩感知恢复问题,对稀疏后四元数矩阵信号Q0按行进行观测得到观测值Y1,对Q0按列进行观测得到观测值Y2;进行解码时,通过观测值Y1和Y2,以及稀疏矩阵W和观测矩阵Φ,得到解码重构恢复的四元数矩阵 和 将四元数矩阵 和 相加求平均,得到结果四元数矩阵 将四元数矩阵 的实部及两个虚部,分别赋值给彩色图像的RGB三个通道得到结果图像。本发明的方法与现有的将彩色图像RGB三个分量数据单独处理的方法相比,只需进行一次压缩感知就能恢复彩色图像,比传统的进行三次压缩感知更节约时间。同时本发明利用了四元数矩阵的幅度和相位信息,结合了三个通道之间的相关性,使得彩色图像的恢复更为精确,这是传统的彩色图像处理技术难以企及的。
附图说明
[0055] 图1为压缩感知的理论框架图;
[0056] 图2为本发明方法与现有方法的比较
流程图;
[0057] 图3为
采样率从60%到85%时,按列进行处理时lena图像的恢复效果比较图;
[0058] 图4为采样率从60%到85%时,按行进行处理时lena图像的恢复效果比较图;
[0059] 图5为采样率从60%到85%时,按行处理与按列处理相结合时lena图像的恢复效果比较图;
[0060] 图6为采样率从60%到85%时,单独按行处理、单独按列处理、按行处理与按列处理相结合的方法所对应的图像的不同恢复效果比较图。
具体实施方式
[0061] 下面结合具体
实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的
修改均落于本
申请所附
权利要求所限定的范围。
[0062] 如图2所示,基于伪牛顿法的四元数域彩色图像压缩感知恢复方法,按照前述步骤进行:
[0063] 步骤1、将彩色图像表示成四元数矩阵的形式,一幅64×64彩色lena图像可以用一个64×64的四元数矩阵来存储其三个通道的值。
[0064] 例如,彩色图像右上角3×3的
像素矩阵是:
[0065]
[0066] 则相对应的四元数矩阵表示为:
[0067]
[0068] 因此,我们可以得到大小为64×64的四元数矩阵Q。
[0069] 步骤2、将四元数矩阵Q在冗余字典下进行稀疏分解,实验中我们取的是大小为64×64
离散小波变换(DWT)矩阵W,根据以下公式进行稀疏分解:
[0070] Q0=W×Q×W′,W′表示矩阵W的转置. (10)
[0071] 例如,彩色图像右上角3×3的像素矩阵稀疏后得到:
[0072]
[0073] 对得到的3×3四元数矩阵中四元数幅度的值接近0的元素,我们将其置零,得到如下更为稀疏的四元数矩阵:
[0074]
[0075] 因此,我们可以得到大小为64×64的稀疏四元数矩阵Q0,并且Q0中的非零元素远远少于N×N。
[0076] 步骤3、用大小为M×N的随机高斯矩阵Φ对稀疏后的四元数信号Q0的按行进行观测,得到M×N的观测值Y1,接着对Q0按列进行观测,得到M×N的观测值Y2。
[0077] 步骤401、对原始四元数矩阵Q进行压缩感知编码以后,我们需要通过解决如下优化问题来解码恢复彩色图像:
[0078]
[0079]
[0080] 其中D代表梯度算子,形式如下:
[0081]
[0082] 步骤402、为了避免lp范数在初始点的不可微,我们将lp范数改写为如下形式:
[0083]-6
[0084] 其中ξ是一个很小的常数,常取1.0×10 。这样(11)式就可以改写为:
[0085]
[0086]
[0087] 其中ε是一个很小的常数,常取1.0×10-6,D代表梯度算子,W′表示W的转置;
[0088] 步骤403、对目标函数J(Q)求导,得到目标函数的梯度ΔJ(Q),通过测量梯度的变化构造一个近似的海森矩阵H(Q)。
[0089] 步骤404、设置迭代的初始值Qt=Φ*×Y(这里的Y分别指Y1和Y2),通过如下公式计算Qt的幅度|Qt|和相位θt:
[0090]
[0091]
[0092]
[0093] θ=arccos(S(Q)/|Q|)
[0094] 例如将得到的64×64的四元数矩阵Qt的右上角3×3的元素改写为幅度和相位的形式:
[0095]
[0096]
[0097] 步骤405、设置迭代的步长因子γ,γ取值常设为(0,1]之间,本发明γ取值为0.7。
[0098] 步骤406、结合梯度ΔJ(Q)和海森矩阵H(Q)用如下伪牛顿迭代式求得经过第一次迭代的值Q1。
[0099] Q1=Qt—γ[H(Qt)]-1ΔJ(Qt) (16)
[0100] 步骤407、将Qt更新成Q1,重新计算梯度ΔJ(Q)和海森矩阵H(Q),重复步骤406到步骤407,直到迭代过程收敛。
[0101] 步骤5、将通过观测值Y1和Y2重构得到的四元数矩阵 和 相加求平均,得到四元数矩阵 将四元数矩阵 的实部及两个虚部,分别赋值给彩色图像的RGB三个通道得到结果图像T。
[0102] 为了验证本发明方法的效果,进行了以下实验:
[0103] 1、实验条件:
[0104] 在一台计算机上进行验证实验,该计算机配置为AMD单核处理器(3200兆赫兹)和4096兆字节随机存取
存储器(RAM),编程语言用的是Matlab(8.0版本)。
[0105] 2、实验方法:
[0106] 使用大小为64×64的彩色图像lena作为实验对象,采用以下步骤完成lena图像的压缩感知恢复。
[0107] 1)采用具体实施方式中的步骤1,将彩色图像二维矩阵转化为四元数二维矩阵Q;
[0108] 2)采用具体实施方式中的步骤2,对四元数矩阵Q进行稀疏,得到稀疏后的四元数矩阵Q0;
[0109] 3)设置不同的采样率,对四元数矩阵Q0的按行进行观测,得到观测值Y1。将公式(14)中的三个约束项的权重设置为
[0110] λ1=0.01,λ2=0.0114,λ3=0.0028,p=0.9,迭代的步长因子r=0.7,用具体实施方式中的步骤401到步骤406,得到恢复矩阵
[0111] 4)设置不同的采样率,对四元数矩阵Q0的按列进行观测,得到观测值Y2。将公式(14)中的三个约束项的权重设置为
[0112] λ1=0.005,λ2=0.0124,λ3=0.0048,p=0.9,迭代的步长因子r=0.7,用具体实施方式中的步骤401到步骤406,得到恢复矩阵
[0113] 5)采用具体实施方式中的步骤5,得到结果图像T。
[0114] 3、实验结果的评价指标:
[0115] 实验中采用如下公式计算彩色图像的峰值
信噪比:
[0116] 两个M×N的彩色图像I和K,它们RGB各个通道的均方差定义为:
[0117]
[0118]
[0119]
[0121]
[0122]
[0123]
[0124] 其中,MAXr,MAXg,MAXb表示图像I和K的RGB各个通道图像点
颜色的最大值,如果每个采样点用8位表示,那么它们就是255。
[0125] 彩色图像I和K的峰值信噪比定义为:
[0126] PSNR=(PSNRr+PSNRg+PSNRb)/3 (19)
[0127] 图3显示了采样率从60%到85%时,按列进行处理时lena图像的恢复效果。表1为采样率与峰值信噪比。
[0128] 表1采样率与峰值信噪比
[0129]采样率 60% 65% 70% 75% 80% 85%
信噪比 18.9260 21.4313 21.7546 22.0619 25.4913 28.6617
[0130] 单独按列进行处理时,原始图像的行和行之间的相关性被割裂了。从图3的恢复结果可以看到,恢复图像的列与列之间存在明显的裂痕。当采样率提高时,这种裂痕的数量变少了,说明在提高采样率时,每列的恢复效果更为精确,可以近似地保持行与行之间相关性,但裂痕仍然会存在。
[0131] 从表1的结果可以看到,当采样率提高时,恢复图像与原始图像的峰值信噪比提高了,这与压缩感知的理论知识是一致的。
[0132] 图4显示了采样率从60%到85%时,按行处理时lena图像的恢复效果。表2为采样率与峰值信噪比。
[0133] 表2采样率与峰值信噪比
[0134]采样率 60% 65% 70% 75% 80% 85%
信噪比 24.0742 24.5024 25.8496 26.0599 29.3833 33.3818
[0135] 单独按行进行处理时,原始图像的列和列之间的相关性被割裂了,从图5的恢复结果可以看到,恢复图像的行与行之间存在明显的裂痕。同单独按行处理一样,当采样率提高时,这种裂痕的数量变少了,说明在提高采样率时,每行的恢复效果更为精确,可以近似地保持列与列之间相关性,但裂痕仍然会存在。
[0136] 从表2的结果可以看到,当采样率提高时,恢复图像与原始图像的峰值信噪比提高了,这与压缩感知的理论知识是一致的。对比表1和表2我们发现,单独按行进行处理比单独按列进行处理得到的恢复图像与原始图像的峰值信噪比要高,这说明按行进行处理时,每行的稀疏度比每列的稀疏度要大,这一点与压缩感知的理论知识也是一致的。
[0137] 图5显示了采样率从60%到85%时,按行处理与按列处理相结合图像的恢复效果。表3为采样率与峰值信噪比。
[0138] 使用按行处理与按列处理相结合的方法,从图5的恢复结果可以看到,恢复图像行与行之间的裂痕、恢复图像列与列之间的裂痕得到了明显的改善,在采样率提高时,恢复图像更为平滑,取得了较理想的结果。
[0139] 从表3的结果可以看到,当采样率提高时,恢复图像与原始图像的峰值信噪比提高了,这与压缩感知的理论知识是一致的。对比表2和表3我们发现,当采用按行处理与按列处理相结合这种方法时,恢复图像与原始图像的峰值信噪比得到了小幅度的改善,这与我们预计的结果是一致的。
[0140] 为了更好的对比不同方法的恢复结果,图6显示了采样率从60%到85%时,单独按行处理、单独按列处理、按行处理与按列处理相结合的方法所对应的图像的不同恢复效果。