一种基于分块压缩感知的光谱图像处理方法技术

本发明专利技术公开了一种基于分块压缩感知的光谱图像处理方法，主要解决传统光谱图像处理中采样成本高以及编码复杂度高的问题，其实现步骤包括：根据谱间相关系数对光谱图像进行分类及差分操作，利用分块压缩感知技术进行采样，改善传统的测量矩阵并利用压缩感知重构算法恢复原始光谱图像。本发明专利技术改善了传统图像处理中采样成本高以及编码复杂高的问题，可以实现在相同采样率的条件下有效提高光谱图像的重构质量。

A spectral image processing method based on block compressed sensing

The invention discloses a block compressed sensing method based on spectral image processing, mainly to solve the high cost and the encoding complexity problems of traditional spectral sampling in the image processing, the method comprises the following steps: according to the correlation coefficient between the spectrum of hyperspectral image classification and differential operation, using block compressed sensing sampling to improve the traditional measurement matrix, and use to restore the original spectral image compression sensing reconstruction algorithm. The invention improves the problems of high sampling cost and high coding complexity in the traditional image processing, and can effectively improve the quality of the spectral image reconstruction under the condition of the same sampling rate.

【技术实现步骤摘要】
一种基于分块压缩感知的光谱图像处理方法
本专利技术属于图像信号处理领域，涉及一种光谱图像处理方法，特别是一种基于分块压缩感知(BlockCompressiveSensing,BCS)的光谱图像处理方法，可以实现在相同采样率的条件下有效提高光谱图像的重构质量。
技术介绍
遥感图像可以根据不同地物的光谱响应特征向人们提供大量的科学数据和动态信息，具有极高的应用价值。然而不同于二维图像，光谱图像呈现出三维特性，数据量十分庞大。在传统的图像处理过程中，压缩发生在数据已经被完整采集到之后，一方面导致了系统复杂度的增加，另一方面，浪费了大量的采样资源。压缩感知作为一种拥有巨大潜力的新兴技术，实现了在采样信号的同时对信号进行压缩，能够大大降低传感器的采样和计算成本。针对海量的光谱数据，传统的压缩感知技术需要占用较高的计算资源，幸运的是，分块压缩感知技术通过将图像数据分块处理可以很大程度上减少能量消耗，降低计算复杂度且容易操作，更适用于在线系统。然而分块压缩感知技术并不能有效利用光谱图像的内在特性，从而进一步降低传感器的采样成本、提高图像的重构质量。
技术实现思路
针对以上问题，本专利技术提出了一种基于分块压缩感知的光谱图像处理方法，通过光谱图像的谱间相关性对光谱图像进行分类以及差分操作，利用分块压缩感知技术进行采样，并将传统的感知矩阵进行加权处理，最终实现在相同采样率的条件下提高光谱图像的重构性能。为实现上述目的，本专利技术所述的基于分块压缩感知的光谱图像处理方法，包括以下步骤：（1）根据光谱图像的谱间相关性对光谱图像进行聚类分组：1a）引入谱间相关系数，并据其确定不同谱段之间的相关性：其中，M为谱段的横向总像素数目,N为谱段的纵向总像素数目,x为像素的横坐标,y为像素的纵坐标,I(i,x,y)表示第i个谱段中坐标(x,y)处的像素值,表示第i个谱段中像素的均值，ri,j表示第i个谱段与第j个谱段之间的相关系数；1b）选择与其它谱段的相关系数和最大的谱段作为初始聚类点；1c）将与步骤1b）选出的聚类点的相关系数大于阈值ε的谱段划分到同一组；1d）重复步骤1b）、1c）直到所有谱段都被归类,经过分类的谱段可以分为三类,即参考谱段,独立重构谱段以及依赖于参考谱段的谱段,分别用R、D和I表示；（2）依次对同一组内的光谱图像进行差分操作：其中，Ψ为稀疏基,Φ为测量矩阵,Θ为Ψ和Φ的乘积,xk为R类谱段，xi为D类谱段，xi可以依赖xk进行重构；（3）利用分块压缩感知技术对经过差分操作的光谱图像进行分块并改进传统的感知矩阵：3a）将经过差分操作的光谱图像分成B×B的子块，边界不足处补零；3b）定义权重矩阵wB＝[wij]B×B，其元素为wij＝1+(i+j-1)×q，并将该矩阵向量化，得到q是一个常数；3c）根据w'B改进传统的感知矩阵：其中，表示原始的感知矩阵，mB表示感知值数目，ΦB表示改进之后的感知矩阵；（4）利用改进后的感知矩阵对分块数据进行采样以及差分重构。本专利技术与现有的技术相比具有以下优点：1、本专利技术根据光谱图像的谱间相关性将光谱图像进行自适应分类，无需事先确定聚类数目，且能够确保相关性较强的谱段划分到同一组。2、本专利技术对经过分类的光谱图像进行差分操作，并利用分块压缩感知技术进行采样，可以有效节约系统资源，编码复杂度低且更容易操作，适用于在线系统。3、本专利技术根据图像中不同频率成分的重要性，对传统的感知矩阵进行加权处理，能够有效提高图像的重构效果，操作简单，易于实现。附图说明图1是本专利技术的总体流程图。图2是本专利技术的基于谱间相关性的光谱图像分类方法流程图。图3是在采样率给定的情况下，不同方法的PSNR性能仿真图。具体实施方式以下参照附图对本专利技术进行详细描述：参照图1，本专利技术基于分块压缩感知的光谱图像处理方法，包括以下步骤：步骤1：根据光谱图像的谱间相关性对光谱图像进行分类。参照图2，本步骤的具体实现如下：1a）引入谱间相关系数，并据其确定不同谱段之间的相关性：其中，M为谱段的横向总像素数目,N为谱段的纵向总像素数目,x为像素的横坐标,y为像素的纵坐标,I(i,x,y)表示第i个谱段中坐标(x,y)处的像素值,表示第i个谱段中像素的均值，ri,j表示第i个谱段与第j个谱段之间的相关系数；1b）选择与其它谱段的相关系数和最大的谱段作为初始聚类点；1c）将与步骤1b）选出的聚类点的相关系数大于阈值ε的谱段划分到同一组；1d）重复步骤1b）、1c）直到所有谱段都被归类，经过分类的谱段可以分为三类,即参考谱段,独立重构谱段以及依赖于参考谱段的谱段,分别用R、D和I表示。步骤2：依次对同一组内的光谱图像进行差分操作。其中，Ψ为稀疏基,Φ为测量矩阵,Θ为Ψ和Φ的乘积,xk为R类谱段，xi为D类谱段，xi可以依赖xk进行重构。步骤3：利用分块压缩感知技术对经过差分操作的光谱图像进行分块并改进传统的感知矩阵。3a）将经过差分操作的光谱图像分成B×B的子块，边界不足处补零；3b）定义权重矩阵wB＝[wij]B×B，其元素为wij＝1+(i+j-1)×q，并将该矩阵向量化，得到q是一个常数；3c）根据w'B改进传统的感知矩阵：其中，表示原始的感知矩阵，mB表示感知值数目，ΦB表示改进之后的感知矩阵。步骤4：利用改进后的感知矩阵对分块数据进行采样以及差分重构。本专利技术的效果可以通过以下仿真结果进一步说明：1、实验条件：本实施例在Matlab语言环境下仿真实现，分别采用DCT矩阵、高斯矩阵作为稀疏基和感知矩阵，分块大小为16，重构方法为OMP，R类、D类和I类谱段的采样率分别为0.4、0.2和0.3,独立重构方法的采样率为0.3。仿真采用的高光谱测试图像来自AVIRIS，分别为LowAltitude的第6景以及第5景，为了方便起见，分别用HS1和HS2标记它们。与此同时，多光谱测试图像为MS1和MS2，它们均为8比特图像且由7个谱段组成，分别是Rio、Paris。为了简便起见，测试图像均使用64×64的子图。根据经验，参数q取为3。2、实验内容A.利用不同感知矩阵对光谱图像进行重构，其结果如图3所示，其中：图3(a)是多光谱图像Pairs的PSNR性能仿真图；图3(b)是高光谱图像LowAltitude第6景中谱段39～45的PSNR性能仿真图；图3表明，对于每个谱段，本专利技术所提方法的PSNR性能均明显优于采用传统未加权处理的感知矩阵的BCS方法。相比该方法，本专利技术所提方法PSNR增益最高可以达到8.4dB，最低也有4.3dB。可见，感知矩阵的好坏对压缩感知的重构性能影响很大，且本专利技术所提方法可以有效提高压缩感知重构性能。B.利用不同方法对HS1中谱段41～44进行重构，这里测试图像的大小为256×256，分块大小为32×32，其结果如表1所示：表1PSNR性能(单位dB)其中，谱段42～44均与谱段41相关性较大，因此它们可以依赖谱段41进行差分重构。表2表明，对于每个谱段，本专利技术方法可以在较低采样率的情况下有效提高图像的重构效果，其中基于权重的独立重构方法性能次之，无权重独立重构方法的性能最差。C.利用不同方法对光谱图像进行分组及差分操作，这里感知矩阵均为本专利技术提出的基于权重的感知矩阵，其PSNR性能如表2所示：表2PSNR性能(单位dB)其本文档来自技高网...

【技术保护点】
一种基于分块压缩感知的光谱图像处理方法，其特征在于，包括如下步骤:1a)引入谱间相关系数，并据其确定不同谱段之间的相关性：

【技术特征摘要】
1.一种基于分块压缩感知的光谱图像处理方法，其特征在于，包括如下步骤:1a)引入谱间相关系数，并据其确定不同谱段之间的相关性：其中，M为谱段的横向总像素数目,N为谱段的纵向总像素数目,x为像素的横坐标,y为像素的纵坐标,I(i,x,y)表示第i个谱段中坐标(x,y)处的像素值,表示第i个谱段中像素的均值，ri,j表示第i个谱段与第j个谱段之间的相关系数；1b)选择与其它谱段的相关系数和最大的谱段作为初始聚类点；1c)将与步骤1b)选出的聚类点的相关系数大于阈值ε的谱段划分到同一组；1d)重复步骤1b)、1c)直到所有谱段都被归类，经过分类的谱段可以分为三类,即参考谱段,独立重构谱段以及依赖...

【专利技术属性】
技术研发人员：肖嵩，牛惠，权磊，杜建超，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61