基于隐空间样例学习的非线性压缩光谱成像方法
【专利摘要】本发明公开了一种基隐空间样例学习的非线性压缩光谱成像方法，主要解决现有技术在利用核函数进行非线性空间下的字典学习中时间复杂度和空间复杂度过高的问题。其实现步骤为：1.对训练样本进行预处理，获得虚拟训练样本；2.通过线性字典学习的方法对虚拟训练样本进行训练，获得稀疏字典；3.随机初始化观测矩阵，通过核压缩感知的方法实现非线性压缩感知光谱成像；4.利用pre?image方法恢复出原信号。实验结果表明：在相同的采样率下，本发明方法同现有KPCA的字典学习方法相比，其重构效果较好，且时间复杂度大大降低，可用于高光谱图像的低速率采样和恢复。
【专利说明】
基于隐空间样例学习的非线性压缩光谱成像方法
技术领域
[0001] 本发明属于图像处理技术领域，特别涉及一种压缩光谱成像方法，可用于高光谱 图像的低成本的获取。
【背景技术】
[0002] 压缩感知是近年来信号处理技术领域中发展起来的一种新的采样理论，通过利用 信号的稀疏可压缩特性，可在远小于传统奈奎斯特采样率的条件下，实现信息的精确恢复。 目前现存的压缩感知的方法均是基于显式的线性稀疏表示模型的。线性稀疏表示模型具有 简单直观，容易理解，容易操作等优势。但是，实际场景信息比较复杂，很难在线性稀疏表示 模型下获得足够稀疏的表示。如果使用线性稀疏表示模型，则获得的稀疏系数的稀疏度较 低，这就使得恢复出原信号所需要的测量值更多。目前，非线性稀疏表示模型下的字典学习 方法，均是通过对线性模型下的字典学习方法的"核化"来进行的。但是，引入了核函数后所 构造的核矩阵的大小依据样本数量而定，往往比较大，在运算处理过程中且需要一直保留， 因此基于'核化'的字典学习有很高的时间复杂度和空间复杂度。
【发明内容】
[0003] 本发明的目的在于针对上述已有技术的不足，提出一种基于隐空间样例学习的非 线性压缩光谱成像方法，以降低成像时间复杂度及空间复杂度。
[0004] 本发明的技术方案是：引入隐空间样例学习的方法，对训练样本进行预处理，得到 隐空间样本，从而运用线性字典学习的方法获得稀疏字典;随机生成观测矩阵，通过核压缩 感知的方法实现非线性压缩感知光谱成像。其实现步骤如下：
[0005] (1)选取三组Ii1Xn2Xn3的高光谱图像，除第10谱段外随机选择η个谱段构造大小 为￥父11训练样本矩阵:￥=[71，72，...，7」，...，711]，其中，'? = 111\112，111\112表示高光谱图像 的大小，η3表示总谱带的数目，yj表示第j个谱段拉成的列向量，j = l，2,. . .，η，η为训练样本 的个数；
[0006] (2)对训练样本矩阵Y进行预处理，获得隐空间样本F;
[0007] 2a)从η个训练样本中随机抽取t个样本，组成一个子样本矩阵，记为YkEMw'其 中 t = 20%n;
[0008] 2b)选定核函数为多项式核函数1^(1，7)=€15(〈1，7>) = (〈1，7>+8)(1，求核函数简化 矩阵:R = k(Y，YR)，其中，C&]T'S为核函数的截距参数，其值S s = 0.5，d为指数参数，其 值为d = 5;
[0009] 2c)求子样本矩阵Yr的核矩阵:B = k(YRT，YRT)，其中FeKixi，T表示转置；
[0010] 2d)对子样本核矩阵B进行特征值分解，即B = U Λ Ut，其中Λ是对角矩阵，U是特征 值所对应的特征向量矩阵；
[0011] 2e)设贡献度ζ = 0.9,计算对角矩阵Σ中每一个正的特征值在所有正的特征值总 和中的贡献率，根据降序排列的特征值的顺序，依次累加正的特征值的贡献率，当累加后的 总和大于贡献度Z时停止，把所累加的特征值的数量记为P;
[0012] 2f)根据上述累加的p个特征值，得出相对应的对角矩阵Λ p以及特征向量矩阵Up;
[0013] 2g)根据上述对角矩阵Ap、特征向量矩阵叫和核函数简化矩阵R计算隐空间样本： F = (Λ；^ ，其中 F e IT' t 表示伪逆；
[0014] (3)利用隐空间训练样本F，采用主成份分析PCA的方法求出稀疏字典Ψ，其中 ψ6：Ε ΛΧΒ；
[0015] (4)设测量值的点数为m，并随机生成大小为mXw的高斯随机矩阵，作为观测矩阵 Φ ;
[0016] (5)把三组高光谱图像的第10个谱段的图像作为测试图像，并把其拉成列向量，分 别记为ei，e2，e3;
[0017] ⑷根据⑷所生成的观测矩阵Φ，对(5)中的三幅测试图像ei，e2， e3进行非线性压 缩成像得到测量值M = fk(< Φ，ei> )，其中ei表示第i组测试样本，i从1到3，fk为选定的核 函数；
[0018] (7)据观测矩阵Φ、样本矩阵Y和稀疏字典Ψ，结合核函数，计算感知矩阵G = fk(< Φ ,Υ>)Ψ ；
[0019] ⑶根据测量值M和感知矩阵G利用最小二乘法计算稀疏系数於=GiM
[0020] (9)根据稀疏系数β和稀疏字典Ψ，利用pre-image方法重构出原图像I，i = l，2, 3,其中每为第i组恢复的图像。
[0021 ]本发明与现有的技术相比有以下优点：
[0022] 1.本发明由于对训练样本进行预处理，使用隐空间样本训练字典，大大降低了字 典学习过程中的时间复杂度以及空间复杂度。
[0023] 2.本发明由于对训练样本进行预处理，可使得一切线性字典学习方法均能应用于 非线性空间下的字典学习中。
【附图说明】
[0024] 图1是本发明的实现流程图；
[0025] 图2是用本发明和现有KPCA字典学习方法在采样率为10 %时，对测试图像 IndianPines的重构效果对比图；
[0026]图3是用本发明和现有KPCA字典学习方法在采样率为10%时对测试图像Moffet的 重构效果对比图；
[0027]图4是用本发明和现有KPCA字典学习方法在采样率为10 %时对测试图像 WashtonDC的重构效果对比图。
[0028]具体实施方法
[0029] 参照图1，本发明的具体实现步骤如下：
[0030] 步骤1.构建训练样本矩阵。
[0031] 选取三组mXn2 Xn3的高光谱图像，除第10谱段外随机选择η个谱段构造大小为W Xn训练样本矩阵Y= [yi，y2,. . .，yj,. . .，yn]，其中，w = niXn2，mXn2表示高光谱图像的大 小，M表示总谱带的数目，yj表示第j个谱段拉成的列向量，j = l，2,. . .，n，n为训练样本的个 数；
[0032] 步骤2.构建隐空间训练样本。
[0033] 2a)从η个训练样本中随机抽取t个样本，组成一个子样本矩阵:其中t = 20%n；
[0034] 2b)选定核函数为多项式核函数1^(1，7)=€1<(〈1，7>) = (〈1，7>+8)(1，求核函数简化 矩阵:R=k(Y，Yr)，其中i? e RW，s为核函数的截距参数，其值为s = 0.5，d为指数参数，其值 为 d = 5;
[0035] 2c)求子样本矩阵Yr的核矩阵:B = k(YRT，Yrt)，其中丑e Kwi，T表示转置；
[0036] 2d)对子样本核矩阵B进行特征值分解，即B = UAUT，Λ是对角矩阵，U是特征值所 对应的特征向量矩阵；
[0037] 2e)设贡献度ζ = 0.9,计算对角矩阵Σ中每一个正的特征值在所有正的特征值总 和中的贡献率，根据降序排列的特征值的顺序，依次累加正的特征值的贡献率，当累加后的 总和大于贡献率z = 0.9时停止。把所累加的特征值的数量记为p;
[0038] 2f)根据上述累加的p个特征值，得出相对应的对角矩阵Ap以及特征向量矩阵U p;
[0039] 2g)根据对角矩阵Ap、特征向量矩阵1^和核函数简化矩阵R计算隐空间训练样本：
表示伪逆。
[0040] 步骤3.利用隐空间训练样本F，求出稀疏字典Ψ，其中Ψ e TTxb。
[0041 ] 线性字典学习的方法有很多，例如?04、100、1(-3￥0、？001^等，考虑实现的复杂性，本 发明采用PCA的方法，其实现步骤如下：
[0042] 3a)计算隐空间训练样本的协方差矩阵：Γ =FTF，其中Te ITyrt ;
[0043] 3b)对协方差矩阵Γ进行特征值分解，即Γ = VΣ Vt，其中Σ是对角矩阵，V是特征 值所对应的特征向量矩阵；
[0044] 3c)将特征向量矩阵U作为稀疏字典，即Ψ =V。
[0045] 步骤4.构造观测矩阵。
[0046] 设测量值的点数为m，随机生成大小为mXw的高斯随机矩阵，将该高斯随机矩阵作
为观测矩P 是观测矩阵Φ的第i行行向量，i = l，2,...，m。
[0047] 步骤b.犾I测讯囹像。
[0048]把三组高光谱图像的第10个谱段的图像作为测试图像，并把其拉成列向量，分别 记为 ei，e2，e3〇
[0049] 步骤6.根据上述步骤(3)所求的字典Ψ和步骤(4)生成的观测矩阵Φ，利用核压缩 感知的方法对步骤(5)中的三幅测试图像 ei，e2，e3进行非线性压缩成像。
[0050] 6a)结合核函数，计算测量值向量M和感知矩阵G:
[0051]
[0052]
[0053] 6b)根据压缩观测方程M = GP和6a)中得到的测量值向量M和感知矩阵G，采用最小 二乘法计算稀疏系数=(:ΓΜ a
[0054] 步骤7.根据稀疏系数β和稀疏字典Ψ重构出原图像。
[0055] 利用pre-image方法，通过如下公式重构出原图像4，i = l ,2,3:
[0056]
[0057] 其中，m表示单位正交基的第1列，1 = 1，2，...，w，w为高光谱图像像素点的个数， W0=[C1，C2, . . .，Cj, . . .，Cn]T，Cj表示ψβ的第j个元素，fk为选定的多项式核函数，flT1是fk 的逆函数。
[0058] 本发明的效果可以通过以下实验进一步说明：
[0059] 1)实验条件
[0060]本实验所用的三组高光谱图像为典型的AVIRIS高光谱数据= IndianPines和 Moffet，以及WashtonDC JndianPines数据是1992年由AVIRIS传感器对印第安那州西北农 业区成像所得，Moffet图像是由1992年8月由AVIRIS传感器对加利福尼亚州的Moffett地区 成像所得;波长范围为〇. 4um~2.5um，共224个光谱带，去掉所有像素为0和不透明的波段后 有200个谱段，空间分辨率为SOmc3WashtonDC图像由HYDICE光谱仪对Washtington DC Mali 地区成像而来，波长范围为0 · 4um~0 · 25um，共210个谱段，预处理后选191个波段，空间分辨 率2.8m。IndianPines图像的大小为 145 X 145 X 200，Moffet图像的大小为 145 X 145 X 200， WashtonDC 图像的大小为 145 X 145 X 191。
[0061 ] 实验仿真环境:采用软件MATLAB 2012R作为仿真工具，CPU是AMD A8-5550M，主频 为2 · IOGHz，内存16G，操作系统为Windows 7旗舰版。
[0062]从每组高光谱图像中随机抽取第10个谱段作为测试图像，其他谱段的图像作为训 练样本。
[0063] 2)仿真内容
[0064]仿真1:在10%的采样率下，分别采用本发明与现有KPCA的字典学习方法对测试图 像IndianPines进行压缩感知仿真实验，实验结果如图2所示，其中：
[0065]图2(a)是测试图像IndianPines第10谱段的原始图像；
[0066]图2(b)是采用现有KPCA的字典学习方法的重构图像，其PSNR为34.8322dB;
[0067] 图2(c)是采用本发明的重构图像，其PSNR为35.361dB;
[0068] 仿真2:在10%的采样率下，分别采用本发明与现有KPCA的字典学习方法对测试图 像Moffet进行压缩感知仿真实验，实验结果如图3所示，其中：
[0069]图3(a)是测试图像Moffet第10谱段的原始图像；
[0070]图3(b)是采用现有KPCA的字典学习方法的重构图像，其PSNR为47.2954dB;
[0071] 图3 (c)是采用本发明的重构图像，其PSNR为48.5724dB;
[0072] 仿真3:在10%的采样率下，分别采用本发明与现有KPCA的字典学习方法对测试图 像WashtonDC进行压缩感知仿真实验，实验结果如图4所示，其中：
[0073]图4(a)是测试图像WashtonDC第10谱段的原始图像；
[0074]图4(b)是采用现有KPCA的字典学习方法的重构图像，其PSNR为44.1374dB;
[0075] 图4(c)是采用本发明的重构图像，其PSNR为45.2161dB;
[0076] 从图2、图3和图4的实验结果可以看出，KPCA的字典学习方法和本发明都能够精确 的重构出原图像;在相同的采样率下，同KPCA的字典学习方法相比，本发明重构图像的PSNR 略高一点，说明本发明重构效果较好。
[0077]仿真4:在1%和10%的不同采样率下，分别采用本发明与KPCA的字典学习方法对 测试图像进行非线性压缩感知仿真实验，给出实验运行时间，实验结果如表1所示。
[0078] 表1两种方法的运行时间
[0080]从表1可以看出，同现有KPCA的字典学习方法相比，本发明的方法运行时间更短， 说明本发明方法能够降低时间复杂度。
【主权项】
1. 一种基于隐空间样例学习的非线性压缩光谱成像方法，包括如下步骤： (1) 选取Ξ组mXmXm的高光谱图像，除第10谱段外随机选择η个谱段构造大小为wX η训练样本矩阵Υ= [yi，y2,. ..，yj,...，yn]，其中，w=niXn2，ni Χπ2表示高光谱图像的大小， Π 3表示总谱带的数目，yj表示第j个谱段拉成的列向量，j = l，2，. . .，η，η为训练样本的个数； (2) 对训练样本矩阵Υ进行预处理，获得隐空间训练样本F; 2a)从η个训练样本中随机抽取t个样本，组成一个子样本矩阵，记为政"W:，其中t = 0.2η； 2b)选定核函数为多项式核函数k(x，y)=fk(<x，y〉）= (<x，y〉+s)d，求核函数简化矩阵： R = k(Y，YR)，其中，W e脏"X'，S为核函数的截距参数，其值为s = 0.5，d为指数参数，其值为d =5； 2c)求子样本矩阵Yr的核矩阵:B = k(YrT，YrT)，其中FT e吸Μ，T表示转置； 2d)对子样本核矩阵W进行特征值分解，即B = UAUT，其中Λ是对角矩阵，U是特征值所对 应的特征向量矩阵； 2e)设贡献度z = 0.9,计算对角矩阵Σ中每一个正的特征值在所有正的特征值总和中 的贡献率，根据降序排列的特征值的顺序，依次累加正的特征值的贡献率，当累加后的总和 大于贡献度Z时停止，把所累加的特征值的数量记为P; 2f)根据上述累加的P个特征值，得出相对应的对角矩阵ApW及特征向量矩阵Up; 2g)根据上述对角矩阵Λ P、特征向量矩阵Up和核函数简化矩阵R计算隐空间训练样本：其中护e联PXH，t表示伪逆. (3) 利用隐空间训练样本F，采用主成份分析PCA的方法求出稀疏字典Ψ，其中Ψε盼""； (4) 设测量值的点数为m，并随机生成大小为mXw的高斯随机矩阵，作为观测矩阵Φ; (5) 把Ξ组高光谱图像的第10个谱段的图像作为测试图像，并把其拉成列向量，分别记 为61心，63; (6) 根据(4)所生成的观测矩阵Φ，对(5)中的Ξ幅测试图像ei，e2，e3进行非线性压缩成 像得到测量值M=fk(<〇，ei>)，其中e康示第i组测试样本，i从巧Ij3，fk为选定的核函数； (7) 根据观测矩阵Φ、样本矩阵Y和稀疏字典Ψ，结合核函数，计算感知矩阵G = fk(<〇, Υ>)Ψ; (8) 根据测量值Μ和感知矩阵G利用最小二乘法，计算稀疏系数公= (9) 根据稀疏系数β和稀疏字典Ψ，利用pre-image方法重构出原图像' i = l ,2,3。2. 根据权利要求1所述的方法，其中步骤（3)利用隐空间训练样本F，采用主成份分析 PCA的方法求出稀疏字典Ψ，按如下步骤进行： 2a)计算隐空间训练样本的协方差矩阵：Γ =fTf，其中Γ e胶胃； 2b)对协方差矩阵Γ进行特征值分解，即Γ =vsyT，其中Σ是对角矩阵，V是特征值所对 应的特征向量矩阵； 2c)将特征向量矩阵U作为稀疏字典，即Ψ = V。3. 根据权利要求1所述的方法，其中步骤(9)中根据稀疏系数β和稀疏字典Ψ，利用pre- image方法重构出原图像?,.。由于选择的多项式核函数是可逆的，则重构图像可通过将原图 像投影到单位正交基上经过核函数映射后再利用核函数的逆函数计算原图像。
【文档编号】G06T5/00GK105844591SQ201610015219
【公开日】2016年8月10日
【申请日】2016年1月11日
【发明人】杨淑媛, 蔡朝东, 焦李成, 刘芳, 马晶晶, 马文萍, 熊涛, 刘红英, 李斌, 金莉
【申请人】西安电子科技大学昆山创新研究院