基于边缘的高光谱图像自适应无损预测编码系统及方法
【专利摘要】本发明公开了一种基于边缘的高光谱图像自适应无损预测编码系统与方法，主要解决现有技术在去相关时由于未考虑高光谱图像的边缘信息，导致压缩去相关不彻底，使得压缩比受限制的问题。其实现过程是：1)利用基于对角的中值预测算法对高光谱图像进行谱内预测，得到谱内预测残差；2)利用基于边缘的谱间预测算法对高光谱图像进行谱间预测，得到谱间预测残差；3)通过比较两种预测模式的预测残差，为每一谱段选择最终预测模式；4)对最终预测模式对应的预测残差进行熵编码，得到无损压缩码流。本发明充分利用图像边缘信息，使去相关更彻底，获得较高的无损压缩比，适用于高光谱数据处理及传输。
【专利说明】基于边缘的高光谱图像自适应无损预测编码系统及方法
【技术领域】
[0001]本发明涉及图像处理【技术领域】，更进一步涉及一种高光谱图像无损编码系统与方法，可用于各种高光谱数据处理及传输。
【背景技术】
[0002]高光谱图像是由成像光谱仪对同一地物在数百个谱段上成像所获取的同时包含空间信息和光谱信息的三维数据立方体，被广泛应用于资源勘探、目标识别、环境保护等方面。由于高光谱图像数据量非常庞大，需要采用有效的压缩技术以便图像的存储和传输。实际应用中对高光谱图像的重建图像质量要求比较高，应对高光谱图像进行无损压缩。
[0003]在现有压缩方法中，以JPEG-LS (Joint Photographic Experts Group Losslessand near-lossless compression of conituous_tone still image)为代表的基于预测的压缩方法，利用了高光谱图像较强的空间相关性进行自适应预测编码，具有很低的复杂度和良好的无损及近无损压缩性能，已广泛应用于星上高光谱图像压缩。但是该方法只能对高光谱图像的一个谱段图像进行谱内空间预测，没有利用高光谱图像的谱间相关性。
[0004]为了去除高光谱图像的谱间冗余，进一步提高压缩效率，通常采用“谱间去相关+熵编码”的压缩方法，即先对当前谱段图像和参考谱段图像进行差分或线性预测，然后再对残差图像进行熵编码。其中当前谱段是指当前需要被编码的谱段图像，而参考谱段是指与当前谱段相关性较强、用来对当前谱段进行预测的谱段图像，一般情况下，为了减小传输边信息及实现无损编码，参考谱段本身应在当前谱段之前完成编码。于是在JPEG-LS的基础之上提出了差分JPEG-LS (Differential JPEG-LS)，该方法先对高光谱图像的相邻谱段进行差分处理，以去除谱间相关性，再对差分后的残差图像进行JPEG-LS无损压缩。
[0005]为了利用高光谱图像的校正特性，Jarno Mielikainen提出了 LUT, (JarnoMielikainen,"Lossless Co`mpression of Hyperspectral Images Using LookupTables", IEEE SIGNAL PROCESSING LETTERS, VOL.13, N0.3, MARCH2006.),该方法是一种基于查找表的压缩方法；Bormin Huang等人在LUT的基础之上增加一个尺度因子，提出了LAIS-LUT, (Bormin Huang, and Y.Sriraja, "Lossless Compression of HyperspectralImagery via Lookup Tables with Predictor Selection", Image and Signal Processingfor Remote Sensing XII1Proc.0f SPIE, Vol.6365，63650L, 2006.)，该方法提供了两个预测值以供选择。
[0006]在实际应用中，要求压缩方法不仅压缩比高，而且时间复杂度和空间复杂度低。因此，被广泛米用的压缩方法是 SLSQ, (Francesco Rizzo, Bruno Carpentieri, GiovanniMotta, and James A.Storer, ^Low-Complexity Lossless Compression of HyperspectralImagery via Linear Prediction", IEEE SIGNAL PROCESSING LETTERS, VOL.12，N0.2，FEBRUARY2005.)，该方法提供了一种基于最小均方误差的谱间线性预测器，去除高光谱图像的大量谱间冗余信息以获得较高的无损压缩比，且线性运算具有较低的时间复杂度和空间复杂度。在SLSQ的基础上又提出了一种新方法SLSQ-OPT (Francesco Rizzo, BrunoCarpentieri, Giovanni Motta, and James A.Storer,^Low-Complexity LosslessCompression of Hyperspectral Imagery via Linear Prediction",IEEE SIGNALPROCESSING LETTERS, VOL.12，N0.2，FEBRUARY2005.)，该方法在谱内预测和谱间预测中自适应选择预测模式，以增加少量复杂度为代价，进一步提高了压缩比。
[0007]上述现有方法均忽略了一个重要因素，即在去相关时由于没有考虑高光谱图像的边缘信息，导致压缩去相关不彻底，使得压缩比受限制。基于边缘的预测方法可以改善预测的性能，这是因为沿边缘方向的像素点之间的相关性更强。
【发明内容】
[0008]本发明的目的在于针对上述已有技术的不足，提出一种基于边缘的高光谱图像自适应预测无损压缩系统及方法，以在去相关时利用图像的边缘特性，避免因去相关不彻底而导致的图像压缩比不够高的缺点。
[0009]为实现上述目的，本发明的技术方案如下:
[0010]一、本发明提出的高光谱图像自适应预测无损压缩系统包括:
[0011]谱内预测单元，用于对输入的高光谱图像进行谱内预测，得到谱内预测值，并对图像真实值和谱内预测值取差值，该差值作为谱内预测残差输出给预测模式选择单元；
[0012]谱间预测单元，用于对输入的高光谱图像进行谱间预测，通过预测模板，利用参考谱段对当前谱段的当前点进行预测，得到谱间预测值，并对图像真实值和谱间预测值取差值，该差值作为谱间预测残差输出给预测模式选择单元；
[0013]预测模式选择单元，用于比较每一谱段分别由谱内预测单元和谱间预测单元得到的谱内预测残差和谱间预测残差，为每一谱段选择最终的预测模式，并输出相应的预测残差给熵编码单元；
`[0014]熵编码单元，用于对最终的预测残差进行无损熵编码，得到无损压缩码流文件，编码结束；
[0015]其特征在于:谱间预测单元包括:
[0016]四个待选的预测模板:0°方向的预测模板、45°方向的预测模板、135°方向的预测模板和180°方向的预测模板，用于根据当前点所在的边缘方向进行选择，为谱间预测提供相关性较强的上下文。
[0017]边缘检测模块:用于利用边缘检测算子对当前谱段的参考谱段进行边缘检测，用I标记参考谱段中的边缘点，用O标记参考谱段中的非边缘点，得到参考谱段的边缘检测图像，并将边缘检测图像输出给边缘分析子模块；
[0018]边缘分析模块:用于根据参考谱段的边缘检测图像，判断参考谱段当前点是否为边缘点及其所在边缘的方向，从四个待选的预测模板中选择一个预测模板，并将选定的预测模板输出给确定最佳上下文子模块；
[0019]确定最佳上下文模块:用于从确定的预测模板中选择Ne个像素点作为最佳上下文，并输出给线性运算子模块，I ^ Ne ^ 12 ；
[0020]线性运算模块:用于根据最佳上下文，通过线性运算得到谱间预测值和谱间预测残差，并将谱间预测残差输出给预测模式选择单元。
[0021]二.本发明提出的高光谱图像自适应预测无损压缩方法，包括如下步骤:[0022](I)输入一幅高光谱图像{B1，B2，- ,Bk,...，BN}，其中 Bk=Uuk, Ii,2,k, - ,Iijjjk-…，IH,w, J表示第k谱段图像，Iuk表示高光谱图像中第k谱段、第i行、第j列的真实像素灰度值，输入的高光谱图像的谱段数为N，宽度为W，高度为H，i=l, 2，…，H，j=l, 2，…,W, k=l，2，...，N;
[0023](2)对输入的一幅高光谱图像{B\B2，-,Bk,…，BN}进行谱内预测，得到谱内预测值成人式人…^…义/丨和谱内预测残差IeiAei/,…，…，&/}，其中第k谱段图
像Bk的谱内预测值为B1}，第k谱段图像Bk的谱内
了页测残差为:eir — ieir:1, I, k> eir:1, I, k> …，eir:1, j, k> …，eir:H，W，lJ，~ ^i,j,k ^ir.1,j,k } ^iri J,k 和
Amk分别表示高光谱图像中第k谱段第i行第j列的像素点的谱内预测值和谱内预测残差，i=l，2，...，H，j=l，2，...，W，k=l，2，...，N;
[0024](3)对输入的一幅高光谱图像{B\B2，-,Bk,…，BN}进行谱间预测，得到谱间预测值{4Λ4Λ...ΛΛ...,爲和谱间预测残差Iein1, ein2，…，eink，…，einN}，其中第k谱段图
像Bk的谱间预测值为:鸟/}，第k谱段图像Bk的谱间
了页测残差为:ein — {δ?η； Ij Ij k> ?in:1, I, k> ，?in:1, j, k> ，?in:H, W, k-^，^imi, j,k ~ ,,.,.A..，An:;.j',女和
ein:uk分别表示高光谱图像中第k谱段第i行第j列的像素点的谱间预测值和谱间预测残差，i=l，2，...，H，j=l，2，...，W，k=l，2，...，N;
[0025](3a)对高光谱图像第一谱段B1进行谱内预测，得到第一谱段的谱间预测值β 1与第一谱段的谱间预测残差eini ；
[0026](3b)对高光谱图像除第一谱段外的其它谱段{B2，-,Bk, -,BnI进行谱间预测，得到其它谱段的谱间预测值{為?2,...,為Λ...,為/}和其它谱段的谱间预测残差{einV..，eink,…，einN}，其中k=2，…，N ;按如下步骤进行:
[0027](3bl)选择第k-Ι谱段BH作为当前谱段Bk的参考谱段，对参考谱段Blrf进行边缘检测，得到的边缘检测图像为=CBH=ICIut1, CI1,…，CIuw…，CIH，w，k_J，i=l，2，…,H, j=l, 2,…，W，其中，为参考谱段当前点Iulri的边缘检测结果，取值为I或0，1表示参考谱段当前点Imi为边缘点，O表示参考谱段当前点Iun为非边缘点；
[0028](3b2)对参考谱段BH的边缘检测图像CBH进行分析，从给定的四个预测模板中选择一个作为最终的预测模板；
[0029](3b3)从最终预测模板中选择Ne个与参考谱段当前点Ii^1像素灰度值最接近的点作为参考谱段的最佳上下文，记为:
同时将当前谱段中与参考谱段上下文具有相同空间位置的一组点id.,，记为当前谱段的最佳上下文，得到最终的最佳上下文为:
{^ο1 - 1 -..^,ir1 ^...^^ko2^'-Jol^1ONc I， , = 1二 ",Ne ；
[0030](3b4)利用最终的最佳上下文，通过线性运算得到谱间预测值和谱间预测残差；[0031](4)比较每一谱段的谱内预测残差和谱间预测残差eink的熵，选择其预测残差熵较小的预测模式作为最终预测模式，并将与之对应的预测残差作为最终预测残差输出给熵编码单元；
[0032](5)对输入的最终预测残差进行自适应无损区间编码，得到无损压缩码流文件，编码结束。
[0033]本发明与现有技术相比较，具有如下优点:
[0034]第一，本发明充分利用高光谱图像的边缘特性进行谱间和谱内去相关，避免了因去冗余不彻底而造成的压缩比过低的现象。
[0035]第二，本发明将谱内预测和谱间预测结合起来，自适应地为每一谱段选择合适的预测模式进行预测，使压缩去相关更为彻底，有效地提高了压缩比。
【专利附图】
【附图说明】
[0036]图1是本发明的压缩系统框图；
[0037]图2是本发明的压缩方法流程图；
[0038]图3是本发明中谱间预测方法子流程图；
[0039]图4是本发明中谱间预测的边缘分析的子流程图；
[0040]图5是本发明中谱内预测的上下文模板示意图；
[0041]图6是本发明中谱间预测待选的四个方向的预测模板示意图。
【具体实施方式】
[0042]参考图1，本发明基于边缘的高光谱图像自适应预测无损压缩系统，包括谱内预测单元、谱间预测单元、预测模式选择单元和熵编码单元。
[0043]所述的谱内预测单元，对输入的高光谱图像利用基于对角的中值预测算法进行谱内预测，将得到的谱内预测残差输出给预测模式选择单元。该谱内预测单元包括:初始化模块、边缘判断模块和谱内预测计算模块，其中，初始化模块设置非边缘点的判决阈值T1和边缘点的判决阈值T2,并将这两个初始化参数输出给边缘判断子模块；边缘判断模块根据非边缘点的判决阈值T1和边缘点的判决阈值T2判断当前点是否
为边缘点，计算出谱内预测值{BIR1 BIR2…BIRK…BIRN}，并将谱内预测值输出给谱内预测
计算子模块；谱内预测计算模块将当前点的真实值与其谱内预测值做差得到谱内预测残差{eir1, eir2,…，eIRK,…，eIRN},并将谱内预测残差{eit1，eir2,…，eIRK,…，eIRN/}输出给预测模式选择单元，其中，输入的高光谱图像谱段数为N，宽为W，高为B1，B2, BN，…，BN。
[0044]所述的谱间预测单元，对输入的高光谱图像{B1，B2，-,Bk, -,BnI利用基于边缘的谱间预测算法进行谱间预测，将得到的谱间预测残差输出给预测模式选择单元。该谱间预测单元包括:边缘检测模块、边缘分析模块、确定最佳上下文模块和线性运算模块，其中:边缘检测模块对当前谱段Bk的参考谱段Bk-1利用边缘检测算子进行运算，得到边缘检测图像CBk'区分出图像中的边缘点及非边缘点；边缘分析模块根据参考谱段的边缘检测图像CBk-1,判断参考谱段当前点Ii^1是否为边缘点及其所在边缘的方向，从四个待选的预测模板中选择一个作为最终的预测模板；确定最佳上下文模块在选定的最终预测模板中选择Ne个最佳的上下文进行谱间预测；线性运算模块利用选定的最佳上下文计算出线性谱间预测系数α、谱间预测值,/，式,2,…,為Λ...,為/}和谱间预测残差{ein\ ein2,…，eink，…
，einN}，并将谱间预测残差Iein1, ein2,…，eink，…，einN}输出给预测模式选择单元，其中，输入的高光谱图像谱段数为N，宽为W，高为H，k=l, 2，…，N。
[0045]所述的预测模式选择单元，用于比较谱内预测单元输出的谱内预测残差和谱间预测单元输出的谱间预测残差，为每一谱段选择一种预测模式为该谱段的最终预测模式，并将其对应的预测残差作为最终预测残差输出给熵编码单元。该预测模式选择单元包括:计算残差熵模块和比较选择模块，其中:计算残差熵模块分别计算谱内预测残差le。ei/r..，eirk,…，eirN}的熵和谱间预测残差Iein1, ein2,…，eink,…，einN}的熵，并将每一谱段的谱内预测残差的熵和谱间预测残差的熵输出给比较选择子模块；比较选择模块，对每一谱段的谱内预测残差熵和谱间预测残差熵进行比较，为每一谱段选择最终的预测模式，并将最终预测模式对应的残差输出给熵编码单元。
[0046]所述的熵编码单元，对输入的最终预测残差进行自适应区间编码，得到无损压缩码流文件，编码结束。
[0047]参照图2，本发明基于边缘的高光谱图像自适应预测无损压缩方法详细步骤如下:
[0048]步骤1，输入高光谱图像{B1，B2，…，Bk，…，BN}，其中，输入的高光谱图像谱段数为N，宽为W，高为H;
[0049]步骤2，对输入的 高光谱图像{B1，B2，…，Bk，…，BN}利用谱内预测方式进行预测，得到谱内预测值和谱内预测残差。
[0050]所述的对高光谱图像(B11B2,-,Bk,-,BnI进行谱内预测可以采用多种现有方法实现，如JPEG-LS、2D-CALIC等，本实施例采用但不局限于基于对角的中值预测(Eran A.Edirisinghe, Satish Bedi, Christos Grecos, "Improvements to JPEG-LS viadiagonal edge based prediction", Proceedings of SPIE Vol.4671 (2002).),基于对角的中值预测的预测上下文模板示意图参照图5，按如下步骤进行:
[0051](2a)设置初始化参数:!\=60，T2=8，其中，T1为非边缘的判决阈值、T2为边缘的判决阈值；
[0052](2b)对当前谱段当前点Iuk利用基于对角的中值预测方法进行谱内预测，即将当前谱段当前点的左上方点IiHk与当前谱段当前点的上方点I1-Uk和当前谱段当前点的左方点Ii^k进行比较，并利用当前谱段当前点的右上方点I1-Uik进行计算，得到当前谱
段当前点的谱内预测值:
[0053](2bl)对于Ii^k≥max(Iiij^k, IiUk)，分以下两种情况得到当前谱段当前点Iuk的谱内预测值:
[0054]若(Iimk-maxai^k，且(abs(Ii,代,-。,j <T2)，则认为当前谱段当前点Iuk是边缘点，得到其谱内预测值:L J,k =7 U -1i—IjJc + h-\j+\,k '
[0055]否贝IJ，认为当前谱段当前点Iuk不是边缘点，得到其谱内预测值:
=min(Av—U’4U.A-);[0056](2b2)对于Ii^k ≤mindi^k，，分以下两种情况得到当前谱段当前点Iuk的谱内预测值1:,
[0057]若(minduo且(abs (Ii,代,-、J <T2)，则认为当前谱段当前点Iuk是边缘点，得到其谱内预测值: ^k:1j,k = !1.j-lM ~ Ii—LjJc + I?-l.j+l.k '
[0058]否贝IJ，认为当前谱段当前点Iuk不是边缘点，得到其谱内预测值:
1 / =max(J,)
[0059](2b3)对于 mindink, I1-1,J;k) <I1-1,j-ι,k<max(Ii；j-1；k, 1卜丨，j，k)，认为当前谱段当前点Ii, j,k不是边缘点，得到其谱内预测值:Lthk = Iij-U- -1l-UiJ- + Il-Um-;
[0060](2c)利用上述计算得到的当前谱段当前点的谱内预测值，计算当前谱段当
前点的谱内预测残差:-仁#，并输出所有谱段的谱内预测残差le^，eir2,…
，eir，…，eir }，其中，eir _ {eir: 1，l，k，eir: 1，1，k，…，eir:1，j，k，…，eir:H，W，lJ 表不弟 k 谱段 B 的谱内预测残差，i=l, 2，...，H, j=l, 2，...，ff, k=l, 2，...，N。
[0061]步骤3，对输入的一幅高光谱图像{B1，B2，…，Bk,…，BnI进行谱间预测，得到谱间预测值和谱间预测残差。
[0062](3a)对高光谱图像第一谱段B1利用基于对角的中值预测进行谱内预测，得到第一谱段的谱间预测值B与第一谱段的谱间预测残差An1 ；
[0063](3b)对知儿谱图像除第一谱段外的其它谱段{Β2，...，Β\...，Βν}进行谱间预测，得到其它谱段的谱间预测值彳足2,…,^/，…，^/}和其它谱段的谱间预测残差Iein2,...，力上...，力/}，其中1^2，~州；
[0064]所述的对高光谱图像除第一谱段外的其它谱段{B2，…，Bk,…，BnI进行谱间预测可以采用多种现有方法实现，如差分JPEG-LS、LUT、SLSQ等，本实施例提出并采用但不局限于基于边缘的谱间预测，
[0065]参照图3，本步骤的实现如下:
[0066](3bl)选择第k-Ι谱段BH作为当前谱段Bk的参考谱段，对参考谱段Bk4进行边缘检测，得到的边缘检测图像为=CBH=ICIut1, CI1,…，CIuw…，CIH，w，k_J，i=l，2，…,H, j=l, 2,…，W，其中，为参考谱段当前点Iulri的边缘检测结果，取值为I或0，1表示参考谱段当前点Imi为边缘点，O表示参考谱段当前点Iun为非边缘点；
[0067]所述的对参考谱段Bk-1进行边缘检测可以采用多种现有方法实现，如Roberts算子、Sobel算子等，本实施例提出并采用但不局限于Canny算子；
[0068](3b2 )根据参考谱段BH的边缘检测图像CBH，从给定的四个预测模板选择其中的一个作为最终的预测模板:
[0069]给定的四个预测模板如图6所示，其中图(6a)是0°方向的预测模板、图(6b)是45°方向的预测模板、图(6c)是135°方向的预测模板、图(6d)是180°方向的预测模板，
[0070]参照图4，选择规则如下:
[0071]当参考谱段当前点Ii^1的边缘检测结果CIuh=O时，选择0°方向的预测模板作为它的预测模板；
[0072]当参考谱段当前点Ii^1的边缘检测结果CIi^1=I,继续执行以下步骤；
[0073](3b2.1)在以参考谱段当前点的边缘检测结果CIi^1为中心的5*5邻域内，寻找与之相连，且值为I的点；
[0074](3b2.2)以参考谱段当前点的边缘检测结果为中点，作四条直线L1、L2、L3、L4，其方向分别为0°、45°、135°和180°，计算步骤(3b2.1)中找到的点分别到直线L1、L2、L3、L4的距离累加和，分别记为:D1、D2、D3、D4 ；
[0075](3b2.3)求出所述4个距离累加和Dl、D2、D3、D4中的最小值，记为Dmin,求出与最小距离累加和Dmin对应的直线，记为最小直线Lmin ；
[0076](3b2.4)从四个方向预测模板中，选择方向与最小直线Lmin相同的预测模板作为最终预测模板；
[0077](3b3)从最终预测模板中选择Ne个与参考谱段当前点Ii^1像素灰度值最接近的点作为参考谱段的最佳上下文，记为:
同时将当前谱段中与参考谱段上下文具有相同空间位置的一组点(4,4,- A，.，/^}，记为当前谱段的最佳上下文，得到最终的最佳上下文为:
?1,...,/$ V",...,/$,...,/二c}，/ = 1，2，…,iVc，其中 Ne 取值为 5 ；
[0078](3b4)利用最终的最佳上下文，通过线性运算得到谱间预测值和谱间预测残差:
[0079](3b4.1)在最小均方误差的准则下，根据最佳上下文计算线性谱间预测系数为:
[0080]
【权利要求】
1.一种基于边缘的高光谱图像自适应无损预测编码系统，包括: 谱内预测单元，用于对输入的高光谱图像进行谱内预测，得到谱内预测值，并对图像真实值和谱内预测值取差值，该差值作为谱内预测残差输出给预测模式选择单元； 谱间预测单元，用于对输入的高光谱图像进行谱间预测，通过预测模板，利用参考谱段对当前谱段的当前点进行预测，得到谱间预测值，并对图像真实值和谱间预测值取差值，该差值作为谱间预测残差输出给预测模式选择单元； 预测模式选择单元，用于比较每一谱段分别由谱内预测单元和谱间预测单元得到的谱内预测残差和谱间预测残差，为每一谱段选择最终的预测模式，并输出相应的预测残差给熵编码单元； 熵编码单元，用于对最终的预测残差进行无损熵编码，得到无损压缩码流文件，编码结束； 其特征在于:谱间预测单元包括: 四个待选的预测模板:0°方向的预测模板、45°方向的预测模板、135°方向的预测模板和180°方向的预测模板，用于根据当前点所在的边缘方向进行选择，为谱间预测提供相关性较强的上下文。 边缘检测模块:用于利用边缘检测算子对当前谱段的参考谱段进行边缘检测，用I标记参考谱段中的边缘点，用O标记参考谱段中的非边缘点，得到参考谱段的边缘检测图像，并将边缘检测图像输出给边缘分析子模块； 边缘分析模块:用于根据参考谱段的边缘检测图像，判断参考谱段当前点是否为边缘点及其所在边缘的方向，从四个待选的预测模板中选择一个预测模板，并将选定的预测模板输出给确定最佳上下文子模块； 确定最佳上下文模块:用于从确定的预测模板中选择Ne个像素点作为最佳上下文，并输出给线性运算子模块，I ^ Ne ^ 12 ； 线性运算模块:用于根据最佳上下文，通过线性运算得到谱间预测值和谱间预测残差，并将谱间预测残差输出给预测模式选择单元。
2.根据权利要求1所述的基于边缘的高光谱图像自适应无损预测编码系统，其特征在于所述的谱内预测单元，包括: 初始化模块:用于设置非边缘点的判决阈值T1和边缘点的判决阈值T2,并将这两个初始化参数输出给边缘判断子模块； 边缘判断模块:用于根据非边缘点的判决阈值1\和边缘点的判决阈值T2判断当前点是否为边缘点，计算出谱内预测值，并将谱内预测值输出给谱内预测计算子模块； 谱内预测计算模块:用于将当前点的真实值与其谱内预测值做差，将差值作为谱内预测残差输出给预测模式选择单元。
3.根据权利要求1所述的基于边缘的高光谱图像自适应无损预测编码系统，其特征在于所述的预测模式选择单元，包括: 计算残差熵模块:用于分别计算每一谱段谱内预测残差的熵和谱间预测残差的熵，并将谱内预测残差的熵和谱间预测残差的熵输出给比较选择子模块； 比较选择模块:用于比较每一谱段的谱内预测残差熵和谱间预测残差熵，为每一谱段选择最终的预测模式，并将相应的预测残差输出给熵编码单元。
4.一种基于边缘的高光谱图像自适应无损预测编码方法，包括如下步骤: (1)输入一幅高光谱图像IB^B2，-,Bk,…，BN}，其中 #={1“，,，…，…，IH,w,J表示第k谱段图像，Iuk表示高光谱图像中第k谱段、第i行、第j列的真实像素灰度值，输入的高光谱图像的谱段数为N,宽度为W，高度为H，i=l, 2，...，H, j=l, 2，...，W，k=l, 2，…，N ; (2)对输入的一幅高光谱图像(B11B2,-,Bk,-,BnI进行谱内预测，得到谱内预测值{式1,為,Λ...,爲Λ...,為/}和谱内预测残差k1，eir2,…，ej…，e/}，其中第k谱段图像Bk的谱内预测值为:為/}，第k谱段图像Bk的谱内预测残差为:eir — ieir:1, I, k> eir:1, I, k> …，eir:1, j, k> …，eir:H，W，lJ，~ ^1.j,k ^it:1,j,k J ^u:1,j,k 和 eir:1，j，k分别表示高光谱图像中第k谱段第i行第j列的像素点的谱内预测值和谱内预测残差，i=l，2，...，H，j=l，2，...，W，k=l，2，...，N; (3)对输入的一幅高光谱图像{B1，B2，-,Bk,-,BnI进行谱间预测，得到谱间预测值{4ΛΑΛ和谱间预测残差Iein1, ein2，…，eink,…，einN}，其中第k谱段图像Bk的谱间预测值为:Bin — , Αη:ι,2,?； >"'> 5''' > W,k }，弟k谱段图像Bk的谱间预测残差为 Iein-1einllJk, Sin； Ij Ij J5,…，ein:1, j, k> …，ein:H，W，lJ，^in:1,j,k -11.'j 'k, Ima.'/,和 ein:1，j，k分别表示高光谱图像中第k谱段第i行第j列的像素点的谱间预测值和谱间预测残差，i=l，2，...，H，j=l，2，...，W，k=l，2，...，N; (3a)对高光谱图像第一谱段B1进行谱内预测，得到第一谱段的谱间预测值^in1与第一谱段的谱间预测残差An1 ; (3b)对高光谱图像除第一谱段外的其它谱段{B2，一，Bk，一，BN}进行谱间预测，得到其它谱段的谱间预测值{為?2,...,為/，…,戽/)和其它谱段的谱间预测残差Iein2,…，eink，…，einN}，其中k=2，…，N ;按如下步骤进行: (3bl)选择第k-Ι谱段Bk-1作为当前谱段Bk的参考谱段，对参考谱段Bk-1进行边缘检测，得到的边缘检测图像为=CBH=ICIm1, CI1,_，…，CIiH,…，CIH，w，k_J，i=l, 2，…,H, j=l, 2,…，W，其中，为参考谱段当前点Iulri的边缘检测结果，取值为I或O，I表示参考谱段当前点Imi为边缘点，O表示参考谱段当前点Iun为非边缘点； (3b2)对参考谱段Bk4的边缘检测图像CBH进行分析，从给定的四个预测模板中选择一个作为最终的预测模板； (3b3)从最终预测模板中选择Ne个与参考谱段当前点1Uh像素灰度值最接近的点作为参考谱段的最佳上下文，记为=Kir1,...,/:=}，同时将当前谱段中与参考谱段上下文具有相同空间位置的一组点记为当前谱段的最佳上下文，得到最终的最佳上下文为KU, (3b4)利用最终的最佳上下文，通过线性运算得到谱间预测值和谱间预测残差。 (4)比较每一谱段的谱内预测残差ej和谱间预测残差eink的熵，选择其预测残差熵较小的预测模式作为最终预测模式，并将与之对应的预测残差作为最终预测残差输出给熵编码单元； (5)对输入的最终预测残差进行自适应无损区间编码，得到无损压缩码流文件，编码结束。
5.根据权利要求4所述的基于边缘的高光谱图像自适应无损预测编码方法，其特征在于步骤(2)所述的对输入的一幅高光谱图像{B1，B2，…，Bk,…，Bn}，k=l, 2，…，N,进行谱内预测，按如下步骤进行: (2a)设置初始化参数:!\=60，T2=8，其中，T1为非边缘的判决阈值、T2为边缘的判决阈值； (2b)对当前谱段当前点Iuk利用基于对角的中值预测方法进行谱内预测，即将当前谱段当前点的左上方点I1-^k与当前谱段当前点的上方点I1-Uk和当前谱段当前点的左方点IiHk进行比较，并利用当前谱段当前点的右上方点Ii^k进行计算，得到当前谱段当前点的谱内预测值 (2bl)对于Ii^k≥maxdi^k，Ii^k)，分以下两种情况得到当前谱段当前点Iyk的谱内预测值: 若(IiHkIax(IiJoI1-Uk)XTl)且(Bbs(IiHk-1s-UkXT2)，则认为当前谱段当前点I^k是边缘点，得到其谱内预测值
6.根据权利要求4所述的基于边缘的高光谱图像自适应无损预测编码方法，其特征在于步骤(3b2)所述的对参考谱段Bk-1的边缘检测图像CBk-1进行分析，从给定的四个预测模板中选择一个作为最终的预测模板，按如下步骤进行: (3b2.1)当参考谱段当前点Iun的边缘检测结果CIm1=O时,选择0°方向的预测模板作为它的预测模板，将选定的预测模板输出给(3b3)，结束步骤(3b2);当参考谱段当前点Ii^1的边缘检测结果CIi^1=I,继续执行以下步骤；(3b2.2)在以参考谱段当前点的边缘检测结果CIi^1为中心的5*5邻域内，寻找与之相连的且值为I的点； (3b2.3)以参考谱段当前点的边缘检测结果为中点，作四条直线L1、L2、L3、L4，其方向分别为0°、45°、135°和180°，计算步骤(3b2.2)中找到的点到直线L1、L2、L3、L4的距离的累加和，分别记为:D1、D2、D3、D4 ； (3b2.4)求出距离的累加和D1、D2、D3、D4中的最小值,记为Dniin，求出与最小距离累加和Dmin对应的直线Lmin ； (3b25)根据提供的四个方向预测模板，选择方向与直线Lmin相同的预测模板作为最终预测模板，将最终预测模板输出给步骤(3b3 )。
7.根据权利要求4所述的基于边缘的高光谱图像自适应无损预测编码方法，其特征在于步骤(3b4)所述的利用最终的最佳上下文，通过线性运算得到谱间预测值和谱间预测残差，按如下步骤进行: (3b4.1)在最小均方误差的准则下，根据最佳上下文计算线性谱间预测系数为:
【文档编号】H04N19/13GK103780904SQ201410032531
【公开日】2014年5月7日 申请日期:2014年1月23日 优先权日:2014年1月23日
【发明者】王柯俨, 李云松, 廖慧琳, 宋娟, 王丽萍, 胡子帆, 刘凯, 张静, 吴宪云 申请人:西安电子科技大学