一种基于分层预测的高光谱图像有损压缩方法技术

本发明专利技术从高光谱图像同时存在空间相关性和谱间相关性，而JPEG2000只消除了高光谱图像的空间冗余入手，公开了一种基于分层预测的高光谱图像有损压缩方法。本发明专利技术采用二维小波变换去除空间冗余，利用高光谱图像的谱间相关性，在解码端生成高质量的SI以减少信道编码所需的码率。将高光谱图像分为GOP区域和non_GOP区域，GOP区域中的WZ波段采用分层预测结构生成SI；non_GOP区域的WZ波段将进行特殊处理，只采用前向相邻的一个已解码波段来生成SI。本发明专利技术公开的基于分层预测的高光谱图像有损压缩方法相较于JPEG2000，在全码率段能获得更良好的率失真性能。的率失真性能。的率失真性能。

【技术实现步骤摘要】
一种基于分层预测的高光谱图像有损压缩方法


[0001]本专利技术涉及图像压缩技术，具体涉及一种基于分层预测的高光谱图像有损压缩方法，属于图像通信领域。

技术介绍

[0002]应用搭载在遥感平台上的传感器对目标地物辐射或反射的电磁波信息进行收集、处理，从而形成遥感影像，这一过程通常称作遥感技术。20世纪80年代出现的高光谱遥感技术，能够在电磁波谱的可见光到红外范围内，同时在空间和光谱方向上成像，为每个像元提供数十至数百个近乎连续的窄波段光谱信息。通常高光谱图像的波段宽度小于10nm，具有很高的光谱分辨率。随着成像光谱技术的发展，高光谱成像具有越来越高的光谱分辨率和空间分辨率，这导致获取的高光谱数据量迅速增加，给数据的存储和传输带来了严峻的挑战。因此，开发一种有效的高光谱图像压缩算法一直是遥感领域的热点和难点问题。目前主流的的压缩方案有两个方向：无损压缩和有损压缩。其中，无损压缩是在有限的压缩比下完美地保存了高光谱图像的信息，但由于数据传输的信道带宽有限，仍不能满足星载压缩的要求。然而有损压缩能够获得较大的压缩比，生成不同质量级别的图像，更加满足实现星载压缩的需求。因此，高光谱图像的有损压缩是遥感图像压缩领域的热点研究方向。

技术实现思路

[0003]本专利技术的目的是为了研究高光谱图像的有损压缩，高光谱图像具有较强的空间相关性及谱间相关性，JPEG2000对高光谱图像的每个波段进行独立编解码，仅去除了空间方向的冗余，本专利技术的方法采用二维小波变换去除空间冗余，采用分层预测去除谱间冗余，相比JPEG2000图像编码标准，本专利技术的方法在解码图像具有相同客观质量的条件下，可以节省更多的码率。
[0004]本专利技术的基本思想是将高光谱图像分为GOP区域和non_GOP区域，GOP区域和non_GOP区域中的WZ波段将采用不同的边信息生成方法，且GOP区域中将编码图像分为两种类型，其中K波段采用JPEG2000进行编解码，WZ波段采用分布式编码方案。具体主要包括以下过程步骤：
[0005](1)GOP设置预处理，设高光谱图像总的波段数为M，首先将高光谱图像连续输入的N波段图像作为一个图像组(Group Of Picture,GOP)，则GOP大小为N，每个GOP中的前N
‑
1个波段设为WZ波段，最后一个波段设为K波段(其中第一个波段的图像即波段0作为K波段，但不属于任何GOP)，当M≠n*N+1(n为正整数)时，这些图像能构成的GOP数量为floor((M
‑
1)/N)，而余下的(M
‑1‑
floor((M
‑
1)/N)*N)个图像不足以构成一个GOP，这时WZ波段将处于两种区域即GOP区域和non_GOP区域；
[0006](2)编码，K波段采用JPEG2000进行编码，WZ波段的编码首先减去该波段的均值，进行二维小波变换，然后对得到的小波系数分别进行量化和比特面分解，最后对每个比特面按从高到低的顺序进行LDPC信道编码，得到的校验信息被储存到缓冲器中；
[0007](3)解码，K波段采用JPEG2000进行解码，得到K波段的重建值，对步骤(1)中GOP区域和non_GOP区域的WZ波段采用不同的边信息生成方式，根据多元线性回归模型(Multivariable Linear Regresson Model,MLRM)对当前波段进行预测，生成联合解码的SI生成联合解码的SI，其次SI和编码端传输过来的校验信息一起进行LDPC信道解码，如果解码不成功，则通过反馈信道向编码端请求更多的校验信息，直到解码成功或者达到最大解码次数，然后将解码的比特面合并后得到的重构信息进行反量化与反变换，再加上该波段的均值，得到当前WZ波段的重建值，最后将K波段重建值与WZ波段重建值组合起来，得到解码后的高光谱图像。
[0008]本专利技术的上述技术方案中，所述的一种基于分层预测的高光谱图像有损压缩方法，GOP区域中的WZ波段采用分层预测结构生成边信息，non_GOP区域的WZ波段将进行特殊处理，只采用前向相邻的一个已解码波段来生成边信息。
[0009]本专利技术的上述技术方案中，所述的一种基于分层预测的高光谱图像有损压缩方法，不同于JPEG2000只去除了高光谱图像的空间冗余，本专利技术方法利用高光谱图像的光谱相关性，利用已解码的图像对当前图像加权预测生成高质量的边信息(Side Information,SI)，因此本专利技术方法能减少解码当前图像所需的码率，在解码图像具有相同客观质量的条件下，可以节省更多码率。
[0010]根据本专利技术的上述方法可以编制执行上述所述的一种基于分层预测的高光谱图像有损压缩方法。
[0011]本专利技术是基于以下思路分析而完成的：
[0012]设高光谱图像总的波段数为M，首先将高光谱图像连续输入的N波段图像作为一个图像组(Group Of Picture,GOP)，则GOP大小为N。每个GOP中的前N
‑
1个波段设为WZ波段，最后一个波段设为K波段(其中第一个波段的图像即波段0作为K波段，但不属于任何GOP)。当M≠n*N+1(n为正整数)时，这些图像能构成的GOP数量为floor((M
‑
1)/N)，而余下的(M
‑1‑
floor((M
‑
1)/N)*N)个图像不足以构成一个GOP，这时WZ波段将处于两种区域即GOP区域和non_GOP区域。GOP区域中的WZ波段采用分层预测结构生成边信息；non_GOP区域的WZ波段将进行特殊处理，只采用前向相邻的一个已解码波段来生成SI。当前WZ波段通过MLRM获得预测系数α，SI的具体生成方式如下：
[0013]设当前波段的原始值为f(x,y,k)(x,y∈[1,m],k∈[1,M
‑
1])，平均值为U
k
，波段的重建值为边信息表示为SI(x,y,k)。
[0014]1)当前WZ波段属于GOP区域，则当前WZ波段的参考波段数量P＝2，α＝[α1,α2]T
，且与它的前后向参考波段距离都为l(l＝1,2,3,4)，当前波段的边信息可表示为：
[0015][0016]2)当前WZ波段属于non_GOP区域，则当前WZ波段的参考波段数量为P＝1，α＝α1，且与它的前向参考波段距离都为l(l＝1)，当前波段的边信息可表示为：
[0017][0018]在编码端，K波段采用JPEG2000进行编码。WZ波段的编码首先减去该波段的均值，进行二维小波变换。图像小波变换后的低频子带可以看作原始图像的低分辨率版本，在图像压缩中具有重要意义。在对原始图像进行小波变换之前减去该图像的均值，得到的低频
系数呈现零均值分布和拥有单一尖峰，符合拉普拉斯分布。而比特平面编码是一种对分布在零处有尖锐峰值的图像渐进压缩的有效方法，因此可以减小重建光谱波段的偏差，提高重建图像的质量。然后对得到的小波系数分别进行量化和比特面分解，最后对每个比特面按从高到低的顺序进行LDP本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于分层预测的高光谱图像有损压缩方法，其特征在于主要包括以下过程步骤：(1)GOP设置预处理，设高光谱图像总的波段数为M，首先将高光谱图像连续输入的N波段图像作为一个图像组(Group Of Picture,GOP)，则GOP大小为8，每个GOP中的前7个波段设为WZ波段，最后一个波段设为K波段(其中第一个波段的图像即波段0作为K波段，但不属于任何GOP)，当M≠n*8+1(n为正整数)时，这些图像能构成的GOP数量为floor((M
‑
1)/8)，而余下的(M
‑1‑
floor((M
‑
1)/8)*8)个图像不足以构成一个GOP，这时WZ波段将处于两种区域即GOP区域和non_GOP区域；(2)编码，K波段采用JPEG2000进行编码，WZ波段的编码首先减去该波段的均值，进行二维小波变换，然后对得到的小波系数分别进行量化和比特面分解，最后对每个比特面按从高到低的顺序进行LDPC信道编码，得到的校验信息被储存到缓冲器中；(3)解码，K波段采用JPEG2000进行解码，得到K波段的重建值，对步骤(1)中GOP区域和non_GOP区域的WZ波段采用不同的边信息生成方式，根据多元线性回归模型(Multivaria...

【专利技术属性】
技术研发人员：杨红，王林，卿粼波，何小海，滕奇志，
申请(专利权)人：四川大学，
类型：发明
国别省市：