基于重加权拉普拉斯稀疏先验的高光谱图像压缩感知方法
【技术领域】
[0001] 本发明涉及一种高光谱图像压缩感知方法，特别是涉及一种基于重加权拉普拉斯 稀疏先验的高光谱图像压缩感知方法。
【背景技术】
[0002] 高光谱图像的光谱信息有助于遥感地物的检测定位与分类识别，然而，高光谱图 像的巨大数据量对图像采集、传输和处理中的软硬资源提出了严格的要求，制约了高光谱 图像的应用。因此，高光谱图像压缩算法是高光谱领域的研宄热点之一。目前，大量普通 图像的压缩方法已成功推广应用到高光谱图像上。然而，这类方法只能压缩已经获取的 图像，无法减少成像过程中的巨大资源需求。近年来，压缩感知成像理论证明，仅需要从 场景中采集少量的线性观测值便可以对原始场景的图像进行高精度的重建。ChengBoLi 等人在文献"Acompressivesensingandunmixingschemeforhyperspectraldata processing,IEEETransactionsonImageProcessing, 2012, 21 (3) :1200 _ 1210"中，针对 高光谱图像每个波段使用单像素相机采集少量线性观测值作为压缩数据。重建过程中，引 入少量的端元光谱，结合线性混合模型，在全变差梯度稀疏的约束下，重建空间连续的丰度 值矩阵。最终，利用线性混合模型重建高光谱图像。然而，该方法使用的稀疏约束忽略了稀 疏信号内非零元素之间的关系，重建精度受限；其次，算法性能严重依赖端元的选择，不当 的端元选择可导致重建彻底失败；此外，算法参数需要调节，适应性差。
【发明内容】
[0003] 为了克服现有高光谱图像压缩感知方法重建精度低的不足，本发明提供一种基于 重加权拉普拉斯稀疏先验的高光谱图像压缩感知方法。该方法随机采集每个像素光谱的少 量线性观测作为压缩数据，建立基于重加权拉普拉斯稀疏先验的压缩感知模型和稀疏正则 化的回归模型，对所建立的模型求解。由于随机采集每个像素光谱的少量线性观测作为压 缩数据，减少了图像采集过程中的资源消耗。重建过程中，引入的重加权拉普拉斯稀疏先验 准确刻画了高光谱图像中的强稀疏性，克服了传统拉普拉斯稀疏先验对非零元素的非均匀 约束。通过经验贝叶斯推理，构建噪声鲁棒的非分离稀疏先验约束，实现了高光谱图像的高 精度重建。小波正交基的使用消除了对端元的依赖性，所有的未知参数均可全自动估计。在 真实的卫星图像Urban，PaviaUniversity以及Indiana数据集上的试验结果表明，当采样 率为0. 15,压缩数据中存在信噪比为10db的强噪声时，本发明获得的峰值信噪比相对于背 景技术压缩感知方法提升4db以上。
[0004]本发明解决其技术问题所采用的技术方案是：一种基于重加权拉普拉斯稀疏先验 的高光谱图像压缩感知方法，其特点是包括以下步骤：
[0005] 步骤一、将包含nb个波段，每个波段包含n5个像素的高光谱图像的每一个波段拉 伸成为一个行向量，所有行向量组成一个二维矩阵，。其中，X的每一列表示每一 个像素对应的光谱，称为光谱维；X的每一行对应一个波段的所有像素值，称为空间维。
[0006] 步骤二、采用列归一化的高斯随机观测矩阵采样高光谱图像X的光谱 维，获得压缩数据~，mb表示压缩后的波段长度，P=mb/nb为采样率。
[0007]G=AX+N(1)
[0008] 其中，iVgRW'6X~表示采样中的噪声。
[0009] 步骤三、利用Haar小波基稀疏化高光谱图像的每一个光谱，X=DY，D为小波正交 基，Y为列稀疏的小波系数矩阵。因此，（1)式表示为G=ADY+N。假设采样过程中噪声服 从入IV(〇，!：,,，/)的矩阵正太分布，则对应的似然函数为 _〇]
【主权项】
1. 一种基于重加权拉普拉斯稀疏先验的高光谱图像压缩感知方法，其特征在于包括以 下步骤： 步骤一、将包含nb个波段，每个波段包含η p个像素的高光谱图像的每一个波段拉伸成 为一个行向量，所有行向量组成一个二维矩阵，尤；其中，X的每一列表示每一个像 素对应的光谱，称为光谱维；X的每一行对应一个波段的所有像素值，称为空间维； 步骤二、采用列归一化的高斯随机观测矩阵』￡腿￥"采样高光谱图像X的光谱维，获 得压缩数据Ge ，mb表示压缩后的波段长度，P = mb/nb为采样率； G = AX+N (1) 其中，WeE?4%表示采样中的噪声； 步骤三、利用Haar小波基稀疏化高光谱图像的每一个光谱，X = DY，D为小波正交基， Y为列稀疏的小波系数矩阵；因此，（1)式表示为G = ADY+N;假设采样过程中噪声服从 ΜΛΓ(0,Σ,,,/)的矩阵正太分布，则对应的似然函数为
其中，2 = 2n=diag(A)表示以λ元素为对角线元素的对角矩阵，指 示噪声的强度；||2|k = Σ"β)表示Q矩阵的加权迹范数; 对于稀疏信号Υ，由于拉普拉斯分布不是（2)式中似然函数共轭先验，采用级联分布的 方式构造重加权拉普拉斯稀疏先验；首先假设Y服从如下分布
其中，2y=diag(Y)表示以γ的元素为对角线元素的对角矩阵，，=[汴...，匕^控 制Y中每一行的稀疏度，γ i= O表示Y的第i行为O ;假设F = [Λ，···，JfifJ，则其中任一列 Yi服从高斯分布；假设超参数γ服从以下的伽马分布，
= ;以上两级先验，等价于重加权拉普拉斯分布，因为对于'有
其中，夂=力喂([·/^?"，…，;米用级联先验，λ，γ和κ均为待估计参数； 步骤四、由于λ，γ和κ未知，因此，根据经验贝叶斯框架，先基于压缩数据G利用MP 估计去知参教λ, Y和κ ,加下
通过积分，并引入_21og运算，容易得知（6)式等价于最小化的（7)式
其中，tr( ·)表示迹范数，Sby= Σ n+AD2yDTAT，￡(Α，，，λγ)为代价函数； 通过变形（7)式，得到稀疏信号Y的非分离稀疏约束模型；首先，对（7)式的第一份部 分进行变形
然后，将（8)式带入到（7)式中，得
接着，引入新的代价方程.，如下
显然，￡( y,；l，/,;c) 2 ￡(AU)，而且能够证明，最小化（7)式再对稀疏信号Y进行MAP 估计，与直接最小化（10)式得到的λ，γ和κ相同，关于Y的解仅相差一个常量^;因 此，（10)式看作是关于稀疏信号Y的正则化回归模型，其中+l〇g|2fe,.I为稀疏信号的非 分离稀疏约束；该约束不能拆分成对于Y中每一行的独立约束，因此该约束能同时约束稀 疏信号中非零元素，潜在地考量非零元素之间的相关性；此外，Sby中包含了表征噪声强度 的λ，因此，得到的稀疏约束随着估计的噪声强度自适应的变化，具有噪声鲁棒性； 步骤五、已知压缩后数据G和采样矩阵Α，采用坐标下降法最小化（10)式，每次迭代中 仅优化一个变量而固定剩余的其他变量；具体步骤如下： ① 初始化，λ °，γ °，κ °均初始化为对应长度的全1向量，计数变量t = 0 ; ② 更新中间变量 Ση= diagU，，Sy= diag(Y，，Sby= Σ n+AD2yDTAT; ③ 固定λ \ γ$Ρ κ %根据（10)式得到关于Y的优化形式，如下
求解得到Y的更新规则如下，
④ 固定Yt+1，λ$ρ κ t，得到关于γ的优化形式，如下
求解得到如下的更新形式：
(14) 其中，；κΓ 为 Yt+1的第 i 个元素，2 = + )代表 VT+Yt+1(Yt+1)T的对角线元素组成的向量，Γ = Σ,.-Σ,.Ζ/，Σ以Ζ)Σ,.，~为z的第i个元素； ⑤ 固定Yt+1，Yt+1和κ %得到关于λ的优化形式，如下
求解得到如下的更新形式： A'" = yjdiagJ^Q^yJa (16) 其中，根号运算表示向量每一个元素开方后组成的向量，./运算代表两个向量对应元 素相除后组成的向量，2 α = 代表对角线元素组成的向量； ⑥ 固定Yt+1，Yt+4P λ t+1，得到关于κ的优化形式，如下
求解得到如下的更新形式： /ri+1=」一 (18) γ + d 前式中的加法和除法运算均作用在向量的每一个元素上，得到一个新向量，d = KT6的 引入是为了确保γ中出现〇时，（18)式依然有意义； ⑦ 计算稀疏信号Y更新前后的差异，如下
其中，表示对Yt+1内的每一个元素乘以M · I If表示弗罗贝尼乌斯范 数，如果计数器t >400或者更新差异η < ιοΛ则退出循环；否则，循环执行步骤②至⑦； ⑧假设上述循环结束得到的最优稀疏信号为ym。，则待重建的高光谱图像IJi过如下 方式得到： υ(Κ"·?) (20)。
【专利摘要】本发明公开了一种基于重加权拉普拉斯稀疏先验的高光谱图像压缩感知方法，用于解决现有高光谱图像压缩感知方法重建精度低的技术问题。技术方案是随机采集每个像素光谱的少量线性观测作为压缩数据，建立基于重加权拉普拉斯稀疏先验的压缩感知模型和稀疏正则化的回归模型，对所建立的模型求解。由于随机采集少量线性观测作为压缩数据，减少了图像采集过程中的资源消耗。重加权拉普拉斯稀疏先验准确刻画了高光谱图像中的强稀疏性，克服了传统拉普拉斯稀疏先验对非零元素的非均匀约束，提高了高光谱图像的重建精度。经测试，当采样率为0.15，压缩数据中存在信噪比为10db的强噪声时，本发明的峰值信噪比相对于【背景技术】方法提升4db以上。
【IPC分类】H03M7-30
【公开号】CN104734724
【申请号】CN201510114261
【发明人】魏巍, 张艳宁, 张磊, 严杭琦
【申请人】西北工业大学
【公开日】2015年6月24日
【申请日】2015年3月16日