基于自适应分块的高光谱图像压缩感知重构方法.pdf

基于自适应分块的高光谱图像压缩感知重构方法王 洋 杨孟宇 赵首博*(哈尔滨理工大学测控技术与通信工程学院 哈尔滨 150080)(广东海洋大学机械与动力工程学院 湛江 524000)摘 要：在对高光谱图像采样重构的研究中，整体采样和固定分块采样没有考虑到高光谱图像复杂的纹理特征分布，使用了相同的测量矩阵导致图像的重构质量较差。针对此问题，该文提出基于2维图像熵自适应分块压缩感知重构方法(ABCS-IE)，该方法以图像2维熵作为高光谱图像纹理细节的度量，根据图像的纹理细节分布自适应改变图像子块的大小，然后为不同的图像块分配特定的采样值，根据分配的采样值设计专有的测量矩阵对图像块进行压缩测量，将采样测量值代入重构算法中进行重构。实验结果表明，与整体采样重构和固定分块采样重构相比，将该方法应用到压缩感知重构算法中对高光谱图像进行采样重构后，重构的图像在视觉效果上有明显的提高，取得的峰值信噪比(PSNR)和结构相似度(SSIM)最大，采样率为0.4时，PSNR提高了24 dB，SSIM最大提高了0.27，均方根误差(RMSE)和信息熵差值(H)也有所降低，说明重构的图像更加接近原始图像。而且运算时间也减少了11.5 s。可见，该方法能充分利用高光谱图像的纹理特征，有效提高图像的重构质量，同时减少重构的运算时间。关键词：高光谱图像重构；自适应分块；压缩感知；图像2维熵中图分类号：TN911.73;TP751文献标识码：A文章编号：1009-5896(2023)07-2605-09DOI:10.11999/JEIT220738Compressed Sensing Reconstruction of Hyperspectral ImagesBased on Adaptive BlockingWANG Yang YANG Mengyu ZHAO Shoubo(School of Measurement and Control Technology and Communication Engineering,Harbin University of Science and Technology,Harbin 150080,China)(School of Mechanical and Power Engineering,Guangdong Ocean University,Zhanjiang 524000,China)Abstract:In the study of sampling and reconstruction of hyperspectral images,global sampling and fixed blocksampling do not take into account the complex texture distribution of hyperspectral images,and the use of thesame measurement matrix results in poor image reconstruction quality.To solve this problem,an AdaptiveBlock Compressed Sensing-Image Entropy(ABCS-IE)is presented.In this method,2-dimensional imageentropy is used as a measure of texture details of hyperspectral images.The size of image blocks is adaptivelychanged according to the texture details distribution of the image.Then,specific sampling values are assignedto different image blocks,and a special measurement matrix is designed to compress the image blocks accordingto the assigned sampling values,and the sampled measurements are brought into the reconstruction algorithmfor reconstruction.The experimental results show that when this method is applied to the compression-awarereconstruction algorithm to sample and reconstruct the hyperspectral image,the visual effect of thereconstructed image is significantly improved.The maximum Peak Signal-to-Noise Ratio(PSNR)andStructural SIMilarity(SSIM)are obtained.When the sampling rate is 0.4,the PSNR is increased by 24 dBand the SSIM is increased by 0.27,the Root Mean Square Error(RMSE)and the information entropydifference(H)are also reduced,indicating that the reconstructed image is closer to the original image.Moreover,the operation time is reduced by 11.5 s.It can be seen that this method can make full use oftexture features of hyperspectral images and improve effectively the quality of image reconstruction,and reduce 收稿日期：2022-06-06；改回日期：2023-01-14；网络出版：2023-02-03*通信作者：赵首博基金项目：国家自然科学基金(61801148)，黑龙江普通本科高等学校青年创新人才培养计划(UNPYSCT-2020187)Foundation Items:The National Natural Science Foundation of China(61801148),The Young Innovative Talents Training Plan ofHeilongjiang Ordinary Undergraduate Colleges and Universities(UNPYSCT-2020187)第45卷第7期电 子 与 信 息 学 报Vol.45No.72023年7月Journal of Electronics&Information TechnologyJul.2023the operation time of reconstruction.Key words:Hyperspectral image reconstruction;Adaptive blocking;Compression Sensing(CS);Image two-dimensional entropy 1 引言高光谱图像是对空间某一区域在几十个甚至数百个波段范围内的同时成像，其每个波段可近似看成一幅自然图像1。由于高光谱图像同时含有目标的空间信息和光谱信息，所以光谱图像已经应用到医学工程2、军事探测3、农业检测4,5、食品安全6以及异常检测7等领域。但高光谱数据是一个庞大的数据集合，图像包含了目标复杂的信息，这对高光谱图像的测量传输和重建造成了一定的难度。因此，在光谱成像技术的研究中，对光谱图像进行有效的压缩测量和重构是十分重要的8。压缩感知(Compressive Sensing,CS)理论可以通过少量的采样值精确地重构出原始目标值，该理论的出现为高光谱图像的压缩重构带来了新思路。压缩感知理论对稀疏或者可压缩的信号提供了一种新的感知框架进行采样，这种采样方式比利用传统的香农-奈奎斯特采样方案更加有效9。CS核心是以远低于奈奎斯特速率的采样率对目标图像进行采样，则原始信号可以通过重构算法由这些少量的测量值重构得到。该理论是利用采样和压缩并行的方式对信号进行压缩测量，打破了采样率对信号采集的限制，减小了后期数据传输和存储的压力10,11。将压缩感知理论应用到光谱数据的采集实验中，可以在数据的采集端去除大量的冗余信息，从而为后期对数据的计算减小压力12。在压缩感知理论基础上发展的光谱成像系统有：单像素相机成像系统13，基于编码孔径的光谱成像系统(Coded,Aperture Snapshot Spectral Imagers,CASSI)14，之后的一些成像系统主要是在编码孔径成像系统中对编码孔径进行的优化15,16。但这些光谱成像系统中，主要是在测量端对光谱图像进行统一单帧测量，然而在以压缩感知理论为基础的图像压缩处理中，将目标图像作为一个整体进行压缩重构时，需要很大的存储空间来存储测量矩阵，也加大了重建过程的难度17。分块压缩感知(BlockCompressed Sensing,BCS)算法18的引入可以有效改善该问题。分块压缩感知算法是对图像进行均匀分块，并使用相同的测量矩阵对各个子块进行单独采样测量，提高了采样效率和重构性能19。但使用分块压缩感知处理图像时，图像主要是被均匀分块的，并使用固定的采样率进行随机采样，但每帧图像各个空间位置的纹理分布都有差异，使用相同的测量矩阵进行采样测量时，尤其在较低的采样率下，会造成纹理细节较少的平滑区域采样点过多，而纹理细节丰富的区域采样点过少的问题，影响后期图像的重构质量。为了解决此问题，许多研究工作将目标聚焦于优化采样测量过程，在采样过程中为不同的图像块动态设置采样率，依据不同的参数对图像块进行分类，然后为不同类别分配不同的采样率20。文献21以图像的纹理显著性作为分配依据，根据显著性信息图的分布给对应的图像块分配采样值；文献22根据图像的空间频率提取图像的纹理特征，以此为不同的图像块分配测量值。但在上述方法中，主要是采用固定尺寸对图像进行均匀分块，在分块时并没有考虑到图像的纹理特征，无法有效地分离出复杂区域和平滑区域，造成图像重构效果不佳。针对上述背景和问题，考虑到高光谱图像数据规模很大，图像内容丰富，对其进行整体测量重构会占用较大的内存，重构的效果也不佳，本文将分块采样测量的方法应用到高光谱图像的采样重建中，提出基于2维图像熵自适应分块压缩感知重构方法(Adaptive Block Compressed Sensing-ImageEntropy,ABCS-IE)。该方法目的是充分利用图像的纹理特征使得分块更加准确，提高各图像块的压缩采样效果，从而改善整体图像的重构效果。本文方法利用高光谱图像复杂的纹理内容，将图像的灰度分布特征和表示灰度空间的特征向量组成图像的2维熵，根据2维熵对图像进行自适应分块，然后将分块后的图像块进行分组，不同组类别根据熵值大小和所占权重自适应分配不同的采样率，分块阈值由各组图像子块2维熵的均值决定。本方法能够充分利用高光谱图像的纹理细节并提高各图像块的压缩采样效果，改善整体图像的重构效果。2 压缩感知理论基础x RN1x=RNN RN1xyMN压缩感知指出，一个N维的离散信号，如果在某变换域中具有稀疏性，则原始信号可以表示为：，其中为稀疏正交基，为稀疏系数，只有少量元素值不为0。则信号 可以由少量的随机观测值 进行重构所得。(MN)为观测矩阵y=x=(1)xyl0式(1)为一个欠定方程，利用压缩感知重构算法求解该方程，将N维的原始信号 从M维的测量信号中重构出来。通常将求解上述方程的过程转换为求解 优化过程2606电 子 与 信 息 学 报第 45 卷=argmin0,s.t.y=x=A(2)其中，常用的压缩感知重构算法主要有贪婪迭代算法、凸优化算法、统计优化方法以及组合算法。B BB RMBB2将压缩感知应用到高光谱图像时，需要使用庞大的观测矩阵对高光谱图像进行测量。为此，结合分块压缩感知算法(BCS)，在分块压缩感知的采样过程中，对大小为NLNW的原图像I进行分块，将其分成互不重叠的大小为的图像块，各子块的测量矩阵为，整个图像的测量矩阵组成为=B0.00B.0.0.0B(3)分块测量的过程为y=y1y2.yb=BBBBx1x2.xb=x1x2.xb(4)即得到的观测向量集合为y=yi|yi=Bxi,i=1,2,.,b(5)xiMBb=(NL NW)/(B B)其中，为第i个子块的灰度值向量，为各图像块的采样数，是图像块的总块数。3 基于2维图像熵自适应分块压缩感知重构方法本方法结合高光谱图像纹理内容丰富的特点，将图像2维熵作为其纹理复杂度的度量，根据图像2维熵自适应选择分块大小，为不同的图像块分配不同的采样值并设计特有的测量矩阵，该分块方式只针对图像的空间维。3.1 图像2维信息熵计算图像的2维熵能够反映图像灰度分布的空间特征，当图像块的2维熵较小时，表明该图像块的像素分布较均匀，即图像块表面相对平滑，细节内容较少，通过少量的采样值可以较容易地重构出原始图像；当图像块的2维熵较大时，表明图像块的像素值分布不均匀，图像块中包含的信息复杂，通过较少的采样值重构出的图像质量较差，需要分配较多的采样值。B Bf(i,j)(i,j=1,2,.,L 1)pij(i,j)假设图像的灰度等级为L(0L1)，图像大小为，用表示图像中灰度值为i、邻域灰度均值为j的像素点出现的频数，为区域灰度均值对发生的概率，则pij=f(i,j)B2(6)图像2维熵的定义为H=L1i=0L1j=0pijlog2pij(7)3.2 自适应分块压缩采样算法以图像2维熵作为依据对图像进行自适应分块，并根据不同块所占权重分配不同的测量次数，图1为基于2维图像熵自适应分块压缩感知采样方法的流程图。自适应分块压缩采样算法主要聚焦在自适应采样测量过程，将该算法与压缩感知重构算法相结合后对图像进行重构，具体重构过程如下：X=NL NW(B/2k1)(B/2k1)Xi(k=1,2,.,k)步骤1将输入图像预分块，分成互不重叠的的图像块，并记录图像块数b，为分块次数。Hi(i=1,2,.,b)Hi=bi=0Hi/b步骤2根据式(7)计算每个图像块的2维熵，并计算各个图像块熵值的平均值。HiTkHiTk32 3232 32步骤3将均值作为每次的分块阈值，比较各子块的熵与阈值的大小，若大于分块阈值，说明纹理细节丰富，继续重复步骤1和步骤2对图像块进行分块处理；若小于分块阈值，停止分块，并保存各个图像子块。经过多次实验表明，当复杂图像块小于时，各区域的图像熵差值已经小于0.1，纹理细节分布无明显差异，无需再继续计算分块。同时为了防止图像块边缘出现严重的块效应，且图像块越小，采样次数越高，失去压缩意义。所以当图像块大小为时，停止分块。Cn(n=1,2,.,n)bnPn=bn/nj=0bj步骤4将步骤3中得到的图像块进行分组处理，相同大小的图像块分为同一组别，每个组别用表示，每一组别内图像块个数记为，计算每组图像块所占的比例。M=SR NL NWRnRn=SR PnMn=Rn(B/2k1)2upper=0.9 (B/2k1)2步骤5假设整幅图像的采样率为SR，测量次数为，各组别内图像块的采样率由表示，根据权重计算每个组别内子块的采样率为，测量次数为，为了防止采样次数过大导致失去压缩意义，设置一个采样上界，若分配的块测量次数超过该上界，则将超出采样上界的次数平均分配给其他组内的图像块，直至所有的块测量次数都小于或等于upper。步骤6根据步骤5获得各组内图像块的测量次第7期王 洋等：基于自适应分块的高光谱图像压缩感知重构方法2607Byi数Mn，构造相应的测量矩阵进行压缩测量，同一组内所有图像块使用相同的测量矩阵，生成各块的测量值向量。y=y1,y2,.,yi步骤7将上述过程得到的观测向量集合利用压缩感知重构算法进行重构，最后整合出原始图像X。4 实验结果与分析 4.1 实验设计512 512 480 480384 256从跨学科计算视觉实验室(InterdisciplinaryComputational Vision Lab,ICVL)提供的自然图像高光谱数据库中选择ulm,BGU,object 3组高光谱测试图像进行实验测试。这3组高光谱图像涵盖了建筑、植物、土壤、布偶、电子器件等内容，在纹理内容分布上也有所差异。所有数据的波长范围为400700 nm，波长间隔为10 nm。实验中，为了验证本文算法，获得更明显的对比展示效果，将3组高光谱图像大小分别截取为,。在实验中，采用不同的测量矩阵进行验证，并将自适应分块采样方法和均匀分块采样方法分别应用在不同的压缩感知重构算法中，对高光谱图像进行重构。选择的重构算法为：基追踪算法(Basis Pursuit,BP)、压缩采样匹配追踪算法(Compressed Sampling Matching Pursuit,CoSaMP)、梯度投影稀疏重建算法(Gradient Pro-jection for Sparse Reconstruction,GPSR)。利用峰值信噪比(Peak Signal-to-Noise Ratio,PSNR)、均方根误差(Root Mean Square Error,RMSE)、结构相似性(Structural SIMilarity,SSIM)、重构图像与原始图像熵的差值(H)作为图像重构的评价指标，并对算法的运行时间(t)进行对比，综合图像的重构效果和算法运行时间来衡量算法的性能。4.2 自适应分块效果分析根据每组图像纹理内容分布的特点，在对图像进行预分块处理时，第1组图像B取256，第2组图像B取240，第3组图像B取128，按照图1所示方法对图像进行分块处理，得到如图2所示的分块结果。观察图2可发现：对于纹理内容特征分布不同的各图像，使用本文提出的分块方法，可以实现纹理复杂区域的块较小，复杂度低的区域块较大的效果。例如图2(a)ulm图像的分块结果中，中间区域内容较复杂，块较小，上半部分细节少、分块结果较大。4.3 不同测量矩阵下的重构结果对比不同的测量矩阵会使重构结果存在差异，为了验证本文所提方法在不同测量矩阵下的适用性，同时选取合适的测量矩阵，选择3种常用的测量矩阵进行测试，分别为高斯矩阵、伯努利矩阵、循环矩阵。表1为使用不同测量矩阵的情况下，采用传统的重构方法和本文提出的方法对图像进行重构得到的结果。通过表1可知，不同的测量矩阵下，通过本文所提方法对图像进行自适应分块重构后，图像的重构质量都得到了明显提高，说明本文提出的方 图 1 基于2维图像熵自适应分块压缩感知采样算法2608电 子 与 信 息 学 报第 45 卷式在不同的测量矩阵下都能适用。同时，由表1可以发现，在高斯矩阵的采样测量下，图像重构的PSNR值最大，而循环矩阵由于在循环构造向量时，元素重复率过高，影响了构造的矩阵向量之间的独立性，导致重构的图像质量低于其他测量矩阵。在以下的对比实验中，选择高斯矩阵对测试图像进行压缩测量。4.4 自适应分块算法应用于不同的重构算法的性能比较将本文ABCS-IE算法和传统的BCS算法分别应用到3种压缩感知重构算法中对测试图像进行重构。图3展示了3种组合方式在不同的采样率下重构ulm图像的PSNR,SSIM。表 1 不同测量矩阵下重构图像的PSNR结果对比(dB)目标图像算法测量矩阵高斯矩阵 伯努利矩阵 循环矩阵ulmCoSaMP31.889 131.446 830.624 0BCS-CoSaMP33.486 433.075 132.543 9ABCS-IE-CoSaMP35.685 734.976 033.480 7BGUCoSaMP30.070 929.689 128.625 4BCS-CoSaMP31.136 830.864 929.795 0ABCS-IE-CoSaMP32.926 832.100 831.543 1objectCoSaMP31.160 530.543 829.267 0BCS-CoSaMP32.456 931.790 230.965 7ABCS-IE-CoSaMP34.596 733.996 132.593 5 图 2 3个测试图像的自适应分块结果 图 3 自适应分块算法重构性能对比结果第7期王 洋等：基于自适应分块的高光谱图像压缩感知重构方法2609从图3可知，在不同的重构算法中，基于本文所提ABCS-IE算法的重构效果比传统CS算法和BCS算法的重构效果好，PSNR和SSIM的数值都有提高。由于采样率较低时，传统的BCS算法在纹理复杂的区域分配的采样值很小，造成重构效果差，而ABCS-IE算法在可以保证整体采样率不变的情况下给复杂区域分配高采样率，使得图像重构效果更好。因此数值在较低的采样率下提高得更加明显。4.5 图像重构效果分析图4展示了3组测试图像重构的视觉效果图，整体采样率为0.4，重构算法选用CoSaMP算法，截取图像的纹理复杂的局部区域于每幅子图的右下角展示。从图4可以看出，基于ABCS-IE算法的重构图像在直观上的效果明显好于传统CS和BCS算法重构图像，而且在纹理复杂的局部区域重构效果更佳。表2和表3为不同采样率下重构图像的RMSE和信息熵差值的结果对比，从中也可以看出，基于ABCS-IE算法的重构图像的均方根误差RMSE也是最小的，重构图像与原图像的信息熵差值(H值)较传统算法也有所降低，说明重构的图像更加接近原始图像。4.6 不同算法的运行时间分析表4为不同重构算法的运行时间对比。可以看到，应用ABCS-IE算法时，重构时间得到了减少，这是因为在传统CS和BCS中，对不同纹理复杂度区域采用相同大小的测量矩阵进行采样，这样不仅使复杂区域获得的采样率较低而导致重构效果差，而且使平滑区域获得的采样率较高而造成计算时间的浪费。但是在本文提出的自适应分块算法中，针对不同块的纹理复杂度进行了自适应采样率的分配，复杂区域采样率高，平滑区域采样率较低，提高图像重构质量的同时，也降低了算法的运行时间。5 结束语本文提出基于2维图像熵自适应分块压缩感知重构方法。本方法将2维图像熵作为高光谱图像纹 图 4 图像重构对比图2610电 子 与 信 息 学 报第 45 卷理复杂度的度量，根据熵值大小和权重，自适应地变化图像块大小，并对不同的图像块分配不同的采样值进行压缩感知重构。通过实验验证，该方法可以根据高光谱图像的纹理特征对图像进行有效的分块，实现纹理复杂区域的图像块较小，平滑区域的图像块较大的效果。在图像重构效果方面，自适应分块压缩采样算法与不同的压缩感知重构算法相结合，可以显著提升传统压缩感知算法的重构性能，重构出的高光谱图像的PSNR和SSIM都有了较大的提高，图像在视觉效果方面也得到了明显的改善，尤其在纹理细节复杂的区域，增强了图像重构的纹理细节。而且本文提出的自适应分块算法在运行时间方面也有明显的改善。因此该方法可以有效地解决传统的算法由于忽略图像的纹理分布特征而造成的重构光谱图像效果不佳的问题。参 考 文 献LIU Lei,SUN Min,REN Xiang,et al.Hyperspectral imagequality based on convolutional network of multi-scaledepthJ.Journal of Visual Communication and ImageRepresentation,2020,71:102721.doi:10.1016/j.jvcir.2019.102721.1FABELO H,ORTEGA S,SZOLNA A,et al.In-vivohyperspectral human brain image database for brain cancerdetectionJ.IEEE Access,2019,7:3909839116.doi:10.1109/ACCESS.2019.2904788.2王建成,朱猛.高光谱侦察技术的发展J.航天电子对抗,3表 2 重构图像的RMSE结果对比目标图像算法采样率0.10.20.30.40.50.60.7ulmCoSaMP0.242 90.205 20.185 00.146 10.125 80.085 20.054 5BCS-CoSaMP0.230 80.185 10.149 00.110 70.105 60.065 30.034 7ABCS-IE-CoSaMP0.211 30.148 70.107 70.071 50.067 00.029 30.009 6BGUCoSaMP0.214 90.194 70.184 40.159 20.121 40.062 10.018 7BCS-CoSaMP0.207 30.174 60.155 40.134 20.092 60.041 70.016 8ABCS-IE-CoSaMP0.196 70.145 30.119 40.097 70.059 40.010 50.007 3objectCoSaMP0.157 40.142 20.122 10.095 60.063 10.023 00.003 1BCS-CoSaMP0.142 80.129 70.102 00.072 50.044 70.012 10.002 6ABCS-IE-CoSaMP0.127 10.106 10.076 30.042 30.020 70.006 30.001 7表 3 重构图像的H结果对比目标图像算法采样率0.10.20.30.40.50.60.7ulmCoSaMP3.288 72.215 71.580 01.003 70.680 80.424 70.265 7BCS-CoSaMP2.889 31.672 81.359 40.612 10.430 10.286 70.220 0ABCS-IE-CoSaMP1.969 51.330 20.529 30.339 30.329 10.133 80.029 9BGUCoSaMP0.951 90.837 70.818 20.628 90.383 20.188 10.060 1BCS-CoSaMP0.878 60.733 40.594 80.474 70.282 80.138 30.007 7ABCS-IE-CoSaMP0.827 30.526 90.437 00.352 60.197 20.077 40.004 5objectCoSaMP1.611 81.479 01.213 00.811 80.398 80.127 60.040 1BCS-CoSaMP1.523 11.247 20.924 80.487 10.291 70.049 00.024 9ABCS-IE-CoSaMP1.285 90.975 70.565 30.455 10.183 80.045 60.009 9表 4 不同重构算法的运行时间t对比(s)算法目标图像ulmBGUobjectBP8.757.547.32BCS-BP8.527.037.25ABCS-IE-BP7.085.876.03CoSaMP6.355.835.80BCS-CoSaMP6.185.795.75ABCS-IE-CoSaMP4.674.344.36GPSR8.428.217.84BCS-GPSR8.368.177.76ABCS-IE-GPSR6.846.736.35第7期王 洋等：基于自适应分块的高光谱图像压缩感知重构方法26112019,35(3):3745.doi:10.16328/j.htdz8511.2019.03.009.WANG Jiancheng and ZHU Meng.Development status ofhyperspectral reconnaissanceJ.Aerospace ElectronicWarfare,2019,35(3):3745.doi:10.16328/j.htdz8511.2019.03.009.LU Bing,HE Yuhong,and DAO P D.Comparing theperformance of multispectral and hyperspectral images forestimating vegetation propertiesJ.IEEE Journal ofSelected Topics in Applied Earth Observations and RemoteSensing,2019,12(6):17841797.doi:10.1109/JSTARS.2019.2910558.4王立国,王丽凤.结合高光谱像素级信息和CNN的玉米种子品种识别模型J.遥感学报,2021,25(11):22342244.doi:10.11834/jrs.20219349.WANG Liguo and WANG Lifeng.Variety identificationmodel for maize seeds using hyperspectral pixel-levelinformation combined with convolutional neural networkJ.National Remote Sensing Bulletin,2021,25(11):22342244.doi:10.11834/jrs.20219349.5邓亚美,王秀娟,杨敏莉,等.成像技术在食品安全与质量控制中的研究进展J.色谱,2020,38(7):741749.doi:10.3724/SP.J.1123.202.DENG Yamei,WANG Xiujuan,YANG Minli,et al.Research advances in imaging technology for food safetyand quality controlJ.Chinese Journal of Chromatography,2020,38(7):741749.doi:10.3724/SP.J.1123.202.6HU Meiqi,WU Chen,ZHANG Liangpei,et al.Hyperspectral anomaly change detection based onautoencoderJ.IEEE Journal of Selected Topics in AppliedEarth Observations and Remote Sensing,2021,14:37503762.doi:10.1109/JSTARS.2021.3066508.7LUO Jiqiang,XU Tingfa,PAN Teng,et al.An efficientcompression method of hyperspectral images based oncompressed sensing and joint OptimizationJ.IntegratedFerroelectrics,2020,208(1):194205.doi:10.1080/10584587.2020.1728625.8SUN Zheng and YAN Xiangyang.Image reconstructionbased on compressed sensing for sparse-data endoscopicphotoacoustic tomographyJ.Computers in Biology andMedicine,2020,116:103587.doi:10.1016/pbiomed.2019.103587.9LI Denghui and WANG Yanhong.Implementation of imageresampling algorithm based on compressed sensingJ.Journal of Physics:Conference Series,2021,1732:012071.doi:10.1088/1742-6596/1732/1/012071.10WANG Rongfang,QIN Yali,WANG Zhenbiao,et al.11Group-based sparse representation for compressed sensingimage reconstruction with joint regularizationJ.Electronics,2022,11(2):182.doi:10.3390/electronics11020182.赵首博,李秀红.基于压缩感知的反射光谱重构算法研究J.光谱学与光谱分析,2021,41(4):10921096.doi:10.3964/j.issn.1000-0593(2021)04-1092-05.ZHAO Shoubo and LI Xiuhong.Research on reflectionspectrum reconstruction algorithm based on compressedsensingJ.Spectroscopy and Spectral Analysis,2021,41(4):10921096.doi:10.3964/j.issn.1000-0593(2021)04-1092-05.12WEI Ziran,ZHANG Jianlin,XU Zhiyong,et al.Optimization methods of compressively sensed imagereconstruction based on single-pixel imagingJ.AppliedSciences,2020,10(9):3288.doi:10.3390/app10093288.13TAO Chenning,ZHU Huanzheng,SUN Peng,et al.Simultaneous coded aperture and dictionary optimization incompressive spectral imaging via coherence minimizationJ.Optics Express,2020,28(18):2658726600.doi:10.1364/OE.396260.14ZHANG Hao,MA Xu,LAU D L,et al.Compressivespectral imaging based on hexagonal blue noise codedaperturesJ.IEEE Transactions on Computational Imaging,2020,6:749763.doi:10.1109/TCI.2020.2979373.15ZHAO Yushi,HE Wenjun,LIU Zhiying,et al.Opticaldesign of an offner coded aperture snapshot spectral imagingsystem based on dual-DMDs in the mid-wave infraredbandJ.Optics Express,2021,29(24):3927139283.doi:10.1364/OE.444460.16WANG Zhongliang,HE Mi,YE Zhen,et al.Reconstructi