本发明专利技术公开了一种基于重加权拉普拉斯稀疏先验的高光谱图像压缩感知方法，用于解决现有高光谱图像压缩感知方法重建精度低的技术问题。技术方案是随机采集每个像素光谱的少量线性观测作为压缩数据，建立基于重加权拉普拉斯稀疏先验的压缩感知模型和稀疏正则化的回归模型，对所建立的模型求解。由于随机采集少量线性观测作为压缩数据，减少了图像采集过程中的资源消耗。重加权拉普拉斯稀疏先验准确刻画了高光谱图像中的强稀疏性，克服了传统拉普拉斯稀疏先验对非零元素的非均匀约束，提高了高光谱图像的重建精度。经测试，当采样率为0.15，压缩数据中存在信噪比为10db的强噪声时，本发明专利技术的峰值信噪比相对于背景技术方法提升4db以上。
【技术实现步骤摘要】
本专利技术涉及一种高光谱图像压缩感知方法，特别是涉及一种基于重加权拉普拉斯 稀疏先验的高光谱图像压缩感知方法。
技术介绍
高光谱图像的光谱信息有助于遥感地物的检测定位与分类识别，然而，高光谱图 像的巨大数据量对图像采集、传输和处理中的软硬资源提出了严格的要求，制约了高光谱 图像的应用。因此，高光谱图像压缩算法是高光谱领域的研宄热点之一。目前，大量普通 图像的压缩方法已成功推广应用到高光谱图像上。然而，这类方法只能压缩已经获取的 图像，无法减少成像过程中的巨大资源需求。近年来，压缩感知成像理论证明，仅需要从 场景中采集少量的线性观测值便可以对原始场景的图像进行高精度的重建。ChengBoLi 等人在文献"Acompressivesensingandunmixingschemeforhyperspectraldata processing,IEEETransactionsonImageProcessing, 2012, 21 (3) :1200 _ 1210"中，针对 高光谱图像每个波段使用单像素相机采集少量线性观测值作为压缩数据。重建过程中，引 入少量的端元光谱，结合线性混合模型，在全变差梯度稀疏的约束下，重建空间连续的丰度 值矩阵。最终，利用线性混合模型重建高光谱图像。然而，该方法使用的稀疏约束忽略了稀 疏信号内非零元素之间的关系，重建精度受限；其次，算法性能严重依赖端元的选择，不当 的端元选择可导致重建彻底失败；此外，算法参数需要调节，适应性差。
技术实现思路
为了克服现有高光谱图像压缩感知方法重建精度低的不足，本专利技术提供一种基于 重加权拉普拉斯稀疏先验的高光谱图像压缩感知方法。该方法随机采集每个像素光谱的少 量线性观测作为压缩数据，建立基于重加权拉普拉斯稀疏先验的压缩感知模型和稀疏正则 化的回归模型，对所建立的模型求解。由于随机采集每个像素光谱的少量线性观测作为压 缩数据，减少了图像采集过程中的资源消耗。重建过程中，引入的重加权拉普拉斯稀疏先验 准确刻画了高光谱图像中的强稀疏性，克服了传统拉普拉斯稀疏先验对非零元素的非均匀 约束。通过经验贝叶斯推理，构建噪声鲁棒的非分离稀疏先验约束，实现了高光谱图像的高 精度重建。小波正交基的使用消除了对端元的依赖性，所有的未知参数均可全自动估计。在 真实的卫星图像Urban，PaviaUniversity以及Indiana数据集上的试验结果表明，当采样 率为0. 15,压缩数据中存在信噪比为10db的强噪声时，本专利技术获得的峰值信噪比相对于背 景技术压缩感知方法提升4db以上。本专利技术解决其技术问题所采用的技术方案是：一种基于重加权拉普拉斯稀疏先验 的高光谱图像压缩感知方法，其特点是包括以下步骤： 步骤一、将包含nb个波段，每个波段包含n5个像素的高光谱图像的每一个波段拉 伸成为一个行向量，所有行向量组成一个二维矩阵，。其中，X的每一列表示每一 个像素对应的光谱，称为光谱维；X的每一行对应一个波段的所有像素值，称为空间维。 步骤二、采用列归一化的高斯随机观测矩阵采样高光谱图像X的光谱 维，获得压缩数据~，mb表示压缩后的波段长度，P=mb/nb为采样率。G=AX+N(1) 其中，iVgRW'6X~表示采样中的噪声。 步骤三、利用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的对角线元素组成的向量，Γ = Σ,.-Σ,.Ζ/，Σ以Ζ本文档来自技高网...
【技术保护点】
一种基于重加权拉普拉斯稀疏先验的高光谱图像压缩感知方法，其特征在于包括以下步骤：步骤一、将包含nb个波段，每个波段包含np个像素的高光谱图像的每一个波段拉伸成为一个行向量，所有行向量组成一个二维矩阵，其中，X的每一列表示每一个像素对应的光谱，称为光谱维；X的每一行对应一个波段的所有像素值，称为空间维；步骤二、采用列归一化的高斯随机观测矩阵采样高光谱图像X的光谱维，获得压缩数据mb表示压缩后的波段长度，ρ＝mb/nb为采样率；G＝AX+N (1)其中，表示采样中的噪声；步骤三、利用Haar小波基稀疏化高光谱图像的每一个光谱，X＝DY，D为小波正交基，Y为列稀疏的小波系数矩阵；因此，(1)式表示为G＝ADY+N；假设采样过程中噪声服从的矩阵正太分布，则对应的似然函数为p(G|Y,λ)=exp{-12||ADY-G||Σn2}(2π)mbnp/2|Σn|np/2,Σn=diag(λ),---(2)]]>其中，Σn＝diag(λ)表示以λ元素为对角线元素的对角矩阵，指示噪声的强度；表示Q矩阵的加权迹范数；对于稀疏信号Y，由于拉普拉斯分布不是(2)式中似然函数共轭先验，采用级联分布的方式构造重加权拉普拉斯稀疏先验；首先假设Y服从如下分布p(Y|γ)=exp{-12||Y||Σy2}(2π)nbnp/2|Σy|np/2,Σy=diag(γ),---(3)]]>其中，Σy＝diag(γ)表示以γ的元素为对角线元素的对角矩阵，控制Y中每一行的稀疏度，γi＝0表示Y的第i行为0；假设则其中任一列yi服从高斯分布；假设超参数γ服从以下的伽马分布，p(γ|κ)=Πi=1nbGamma(1,2κi)=Πi=1nbκi2exp(-κiγi2)---(4)]]>以上两级先验，等价于重加权拉普拉斯分布，因为对于yi有p(yi)=∫p(yi|γ)p(γ|κ)dγ=12Πj=1nbκj1/2exp(-||Kyi||1)---(5)]]>其中，采用级联先验，λ，γ和κ均为待估计参数；步骤四、由于λ，γ和κ未知，因此，根据经验贝叶斯框架，先基于压缩数据G利用MAP估计未知参数λ，γ和κ，如下{λopt,γopt,κopt}=argmaxλ,γ,κp(λ,γ,κ|G)∝argmaxλ,γ,κ∫p(G|Y,λ)p(Y|γ)p(γ|κ)dY---(6)]]>通过积分，并引入‑2log运算，容易得知(6)式等价于最小化的(7)式其中，tr(·)表示迹范数，Σby＝Σn+ADΣyDTAT，为代价函数；通过变形(7)式，得到稀疏信号Y的非分离稀疏约束模型；首先，对(7)式的第一份部分进行变形tr(nnp-1GTΣby-1G)=tr[(GTnp)Σby-1(Gnp)]=minY||ADY-Gnp||Σn2+||Y||Σy2.---(8)]]>然后，将(8)式带入到(7)式中，得接着，引入新的代价方程如下显然，而且能够证明，最小化(7)式再对稀疏信号Y进行MAP估计，与直接最小化(10)式得到的λ，γ和κ相同，关于Y的解仅相差一个常量因此，(10)式看作是关于稀疏信号Y的正则化回归模型，其中为稀疏信号的非分离稀疏约束；该约束不能拆分成对于Y中每一行的独立约束，因此该约束能同时约束稀疏信号中非零元素，潜在地考量非零元素之间的相关性；此外，Σby中包含了表征噪声强度的λ，因此，得到的稀疏约束随着估计的噪声强度自适应的变化，具有噪声鲁棒性；步骤五、已知压缩后数据G和采样矩阵A，采用坐标下降法最小化(10)式，每次迭代中仅优化一个变量而固定剩余的其他变量；具体步骤如下：①初始化，λ0，γ0，κ0均初始化为对应长度的全1向量，计数变量t＝0；②更新中间变量Σn＝diag(λt)，Σy＝diag(γt)，Σby＝Σn+ADΣyDTAT；③固定λt，γt和κt，根据(10)式得到关于Y的优化形式，如下求解得到Y的更新规则如下，Yt+1=ΣyDTATΣby-1Gnp---(12)]]>④固定Yt+1，λt和κt，得到关于γ的优化形式，如下求解得到如下的更新形式：&...
【技术特征摘要】
【专利技术属性】
技术研发人员：魏巍，张艳宁，张磊，严杭琦，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西;61
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。