利用格型递归最小二乘滤波器组的高光谱图像压缩

利用格型递归最小二乘滤波器组的高光谱图像压缩 Lossless compression of hyperspectral images using recursive least square lattice filter group 1 2 first-author 郑　铁（1995-），男，黑龙江伊春人，博士研究生，2017年于东北农业大学获得学士学位。主要从事高光谱图像数据处理方面的研究。E-mail：z_hreo@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332877&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332879&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332878&type= 22.01333618 32.00400543 郑　铁 （1995-），男，黑龙江伊春人，博士研究生，2017年于东北农业大学获得学士学位。主要从事高光谱图像数据处理方面的研究。E-mail： z_hreo@163.com z_hreo@163.com 1 corresp E-mail： xuechangbin@nssc.ac.cn E-mail： xuechangbin@nssc.ac.cn 薛长斌（1972-），男，辽宁锦州人，研究员，博士生导师，1997年获得哈尔滨工业大学硕士学位。2017年于北京理工大学获得博士学位。主要研究方向是空间数据处理。E-mail：xuechang bin@nssc.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332880&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332882&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332881&type= 22.01333618 32.00399780 薛长斌 （1972-），男，辽宁锦州人，研究员，博士生导师，1997年获得哈尔滨工业大学硕士学位。2017年于北京理工大学获得博士学位。主要研究方向是空间数据处理。E-mail：xuechang bin@nssc.ac.cn bin@nssc.ac.cn 3 宋金伟（1983-），男，内蒙古呼和浩特人，高级工程师，主要研究方向图像处理与人工智能。E-mail： ysjk2000@qq.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332883&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332885&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332884&type= 22.01334000 32.00408936 宋金伟 （1983-），男，内蒙古呼和浩特人，高级工程师，主要研究方向图像处理与人工智能。E-mail： ysjk2000@qq.com ysjk2000@qq.com 中国科学院 国家空间科学中心， 北京100190 中国科学院 国家空间科学中心， 北京100190 National Space Science Center， Chinese Academy of Sciences， Beijing 100190， China National Space Science Center， Chinese Academy of Sciences， Beijing 100190， China 中国科学院大学， 北京 100190 中国科学院大学， 北京 100190 The University of Chinese Academy of Sciences， Beijing 100190， China The University of Chinese Academy of Sciences， Beijing 100190， China 国家电网有限公司大数据中心， 北京 100031 国家电网有限公司大数据中心， 北京 100031 Big Data Center， State Grid Corporation of China， Beijing 100031， China Big Data Center， State Grid Corporation of China， Beijing 100031， China 自适应递归最小二乘滤波器具有预测准确、收敛速度快的特点，该滤波器被多种高光谱图像无损压缩方案作为重要组成部分。然而传统递归最小二乘滤波器无法快速找到每个谱带的最优预测长度，其压缩方案的性能有待提升。针对该问题，本文提出基于格型递归最小二乘滤波器组的高光谱图像压缩方案。首先，该方案使用单边高斯预测器对待测像素点做谱带内预测，去除图像的空间相关性。其次，采用格型滤波器组筛选出每个谱带的最优滤波器，获得预测误差。并根据格型滤波器组链式序列更新的特点，简化最优滤波器的筛选过程，大幅度降低计算复杂度。最后对预测误差做算术编码。以AVIRIS 2006高光谱图像为测试数据集，本文算法对16位校准图像、16位未校准图像的平均压缩结果分别为3.34 bits/pixel和5.61 bits/pixel。该算法在获得良好压缩结果的情况下，计算时间低于同类别的其余算法。 The lossless compression algorithms for hyperspectral images based on the recursive least square filter have accurate prediction and fast convergence. However， the traditional compression algorithm is incapable of quickly identifying the optimal prediction length of each band， and the performance of the algorithm needs to be improved. In this paper， we propose a hyperspectral image compression algorithm based on the recursive least square lattice filter group. Considering the spatial correlation of the hyperspectral image， the proposed scheme first uses a single-sided Gaussian predictor for inner-band prediction. The lattice filter group is then used to calculate the optimal filter for each spectral band to obtain the prediction error. Further， the screening process of the optimal filters are simplified according to the characteristics of the chain sequence update of the lattice filter group； this significantly reduces the computational complexity. Finally， the arithmetic code is used to encode the residual data. Experiments on the Airborne Visible/Infrared Imaging Spectrometer （AVIRIS） 2006 hyperspectral image datasets reveal that the average compression results of the algorithm for 16-bit calibrated and uncalibrated images are 3.34 bits/pixel and 5.61 bits/pixel， respectively. The proposed algorithm therefore shows a competitive performance with regard to compression and requires less computation time compared to other outstanding algorithms. 高光谱图像 无损压缩 递归最小二乘法 高斯预测器 格型滤波器组 hyperspectral imagery lossless compression recursive least squares Gaussian predictor lattice filter group 1 引　言 1 高光谱图像被广泛应用于地质勘探、环境分析、农业监测等领域。随高光谱图像分辨率的提升，图像数据量急剧增加，给传输带宽和存储容量带来巨大挑战。因此，高效的高光谱图像压缩技术必不可少 ［ 1 1 - 5 5 1-5 ］ 。高光谱图像的压缩方法分为有损压缩和无损压缩。有损压缩虽然能达到很高的压缩比，但会破坏数据的细节信息，这对需要精确分析的计算场景是难以接受的，只有无损压缩能保证图像信息的完整性 ［ 6 6 - 7 7 6-7 ］ 。 高光谱图像同时具备空间相关性和谱间相关性。传统的二维图像压缩技术只针对空间相关性，无法利用谱间相关性提升压缩比 ［ 8 8 - 10 10 8-10 ］ 。目前，常见的高光谱图像压缩方法分为变换类、矢量量化类和预测类。其中变换类算法是先将图像从空域映射到变换域，使用较少码字表示图像的主要信息，随后再对细节量化编码，因此常用于有损压缩 ［ 11 11 ］ 。矢量量化类算法的计算复杂度会随维数增加呈现指数级的增长，庞大的计算量限制了该类算法的实际应用 ［ 12 12 ］ 。预测类算法先使用预测器计算待测像素点估计值与实际值的差异，去除图像的相关性，再针对预测误差做压缩编码。这类算法在获取良好压缩结果的同时保持较低的算法复杂度，逐渐成为国内外研究热点 ［ 13 13 ］ 。 Mielikainen等人提出基于查找表（Lookup Table， LUT）的预测方案 ［ 14 14 ］ ，该方案将前一波段中与当前像素点具有同位坐标的像素作为索引，使用查找表返回预测结果。H. Bormin等人改进LUT的预测方案，提出局部平均谱间尺度变换的方法（Locally Averaged Inter-Band Scaling-LUT，LAIS-LUT） ［ 15 15 ］ 。该方法先计算待测像素点的LAIS估计值，随后在两个查找表中索引与LAIS估计值最接近的值并作为预测结果。基于LUT的压缩方法虽然具有良好的预测效果，然而它消耗大量的内存资源，不适用于星载图像压缩系统。 聚类差分脉冲编码调制（Clustered Difference Pulse Code Modulation， C-DPCM）是一种将聚类和线性滤波器相结合的预测方案 ［ 16 16 ］ 。该方法首先利用k-means聚类算法将光谱矢量做聚类处理，提升同类矢量间的相似度。随后，采用预测长度固定的线性滤波器做预测压缩。具有自适应预测长度的聚类算法（C-DPCM with Adaptive Prediction Length， C-DPCM-APL）是C-DPCM算法的扩展，它通过搜寻线性滤波器输入的最优预测长度获得良好的压缩效果 ［ 17 17 ］ 。但该方案的计算复杂度较高、消耗资源较大，更适合被应用于离线压缩场景。 空间数据系统协商委员会（Consultative Committee for Space Data Systems， CCSDS）提出基于最小均方误差滤波器的快速无损（Fast lossless， FL）算法 ［ 18 18 ］ ，该算法利用先前的波段预测当前像素，并在线性预测过程中自适应改变预测系数。三阶谱间预测器（The third-order inter-band predictor， IP3） ［ 19 19 ］ 由两阶段预测器和熵编码组成，通过谱间预测器和反向像素搜索（BPS）方案计算最终预测值。宋金伟等人提出基于递归最小二乘法（Recursive Least Squares， RLS）的压缩方案，该方案先将待测像素点做空间滤波，再利用固定长度的递归最小二乘做谱间预测 ［ 20 20 ］ 。由于递归最小二乘滤波器较最小均方误差滤波器的收敛速度快，因此其压缩结果要优于FL算法。高放和郭树旭提出自适应预测波段的传统最小二乘压缩算法（Cnventional Recursive Least-Squares Predictor with Adaptive prediction bands， ACRLS） ［ 21 21 ］ 。该文献指出不同待测波段的最优预测长度是不尽相同的，并通过穷举法搜索每个待测点RLS滤波器的最优预测长度，获得良好的压缩效果。但该算法的计算复杂度过高，消耗大量的计算时间。宋金伟等人提出的快速RLS （Fast Adaptive-Length-Prediction RLS， Fast-RLS-ALP） 预测算法 ［ 22 22 ］ ，相较于传统RLS算法，该算法充分利用已知波段的滤波器做前向预测，降低计算复杂度。但该算法的空间预测方案比较简单，且谱带间滤波器只使用固定长度的预测方案，因此压缩结果并不突出。 本文提出一种基于格型递归最小二乘滤波器组的高光谱图像无损压缩方案（RLS Lattice Filter Group， LFG-RLS）。该压缩方案采用的格型滤波器组充分考虑最优谱带预测长度对预测效果的影响，合理地为每个谱带选择最优滤波器，并利用其链式序列更新的特点，对滤波器组分类、简化，大幅度降低计算时间。 2 压缩方案介绍 1 本文提出的压缩方案分为预测器和编码器，其中预测器包含谱带内预测和谱间预测两个阶段。对于高光谱图像的每个谱带，首先使用单边高斯预测器做谱带内预测，得到谱带内预测误差。其次，采用格型滤波器组筛选每个谱带的最优滤波器，使压缩结果接近最优，并根据格型滤波器组链式更新的特点，对最优滤波器做进一步筛选，大幅度降低计算时间。最后，利用算术编码器对预测误差做压缩编码，得到压缩码流。本文的压缩方案结构图如 图1 图1 所示。 2.1　谱带内预测 2 高光谱图像可被视为一组三维数据，设其尺寸为 <math id="M1"><mi>W</mi><mo>×</mo><mi>H</mi><mo>×</mo><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332891&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332890&type= 19.81199837 2.87866688 （行 <math id="M2"><mo>×</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332893&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332892&type= 3.72533321 2.79399991 列 <math id="M3"><mo>×</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332895&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332894&type= 3.72533321 2.79399991 波段），本文使用 <math id="M4"><msub><mrow><mi>s</mi></mrow><mrow><mi>z</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332897&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332896&type= 12.10733318 4.74133348 表示第 <math id="M5"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332900&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332898&type= 1.69333339 2.79399991 波段中第 <math id="M6"><mi>x</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332905&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332903&type= 2.03200006 2.79399991 行第 <math id="M7"><mi>y</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332909&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332907&type= 1.86266661 3.72533321 列像素点的灰度值，或使用 <math id="M8"><msub><mrow><mi>s</mi></mrow><mrow><mi>z</mi></mrow></msub><mo stretchy="false">(</mo><mi>n</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332914&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332912&type= 7.95866680 4.48733330 表示 <math id="M9"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332919&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332917&type= 1.69333339 2.79399991 波段中第 <math id="M10"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332923&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332921&type= 1.86266661 2.79399991 个像素点的灰度值，其中 <math id="M11"><mi>n</mi><mo>=</mo><mi>W</mi><mo>×</mo><mi>y</mi><mo>+</mo><mi>x</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332928&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332925&type= 24.21466637 3.72533321 。并将具有相同 <math id="M12"><mi>y</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332932&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332930&type= 1.86266661 3.72533321 坐标的平面称为线谱平面（Line-Spectral Plate， LSP），如 图2 图2 所示。 在谱带内预测阶段，通过对待测像素点的上下文像素做预测运算，计算谱带内预测误差，从而去除高光谱图像的空间相关性，其中待测像素点的上下文窗口如 图3 图3 所示。ACRLS ［ 21 21 ］ 和Fast-RLS-ALP ［ 22 22 ］ 算法分别采用24邻域和4邻域均值的预测方式去除空间冗余信息。 因待测像素点与上下文像素之间的相关性随距离的增加而减小，高斯预测器被用于谱带内预测。为保证预测过程的因果性，即使用已知像素点预测当前待测点，采用的单边高斯预测器如 式（1） 式（1） 所示： <math display="block" id="M13"><mtable><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi>r</mi></mrow></msub><mo>=</mo></mtd></mtr><mtr><mtd><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>g</mi></mrow><mrow><mo>-</mo><mi>r</mi><mo>,</mo><mo>-</mo><mi>r</mi></mrow></msub></mtd><mtd><mo>⋯</mo></mtd><mtd><msub><mrow><mi>g</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mo>-</mo><mi>r</mi></mrow></msub></mtd><mtd><msub><mrow><mi>g</mi></mrow><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mo>-</mo><mi>r</mi></mrow></msub></mtd><mtd><msub><mrow><mi>g</mi></mrow><mrow><mn mathvariant="normal">1</mn><mo>,</mo><mo>-</mo><mi>r</mi></mrow></msub></mtd><mtd><mo>⋯</mo></mtd><mtd><msub><mrow><mi>g</mi></mrow><mrow><mi>r</mi><mo>,</mo><mo>-</mo><mi>r</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>g</mi></mrow><mrow><mo>-</mo><mi>r</mi><mo>,</mo><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub></mtd><mtd><mo>⋯</mo></mtd><mtd><msub><mrow><mi>g</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub></mtd><mtd><msub><mrow><mi>g</mi></mrow><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub></mtd><mtd><msub><mrow><mi>g</mi></mrow><mrow><mn mathvariant="normal">1</mn><mo>,</mo><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub></mtd><mtd><mo>⋯</mo></mtd><mtd><msub><mrow><mi>g</mi></mrow><mrow><mi>r</mi><mo>,</mo><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>g</mi></mrow><mrow><mo>-</mo><mi>r</mi><mo>,</mo><mn mathvariant="normal">0</mn></mrow></msub></mtd><mtd><mo>⋯</mo></mtd><mtd><msub><mrow><mi>g</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1,0</mn></mrow></msub></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mo>⋯</mo></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mo>⋯</mo></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mo>⋯</mo></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mo>⋯</mo></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr></mtable></mrow></mfenced></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332953&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332951&type= 71.45866394 41.99466705 ， 其中： <math display="block" id="M14"><msub><mrow><mi>g</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msup><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><msup><mrow><mi mathvariant="normal">e</mi></mrow><mrow><mo>-</mo><mfrac><mrow><mfenced separators="|"><mrow><msup><mrow><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>j</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn><msup><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332958&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332955&type= 28.95599937 9.31333351 . 则当前像素点的估计值 <math id="M15"><msub><mrow><mover accent="true"><mi>s</mi><mo>̂</mo></mover></mrow><mrow><mi>z</mi></mrow></msub><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332962&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332960&type= 11.26066685 4.82599974 如 式（3） 式（3） 所示： <math display="block" id="M16"><msub><mrow><mover accent="true"><mi>s</mi><mo>̂</mo></mover></mrow><mrow><mi>z</mi></mrow></msub><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced><mo>=</mo><mfrac><mrow><munder><mo movablelimits="false">∑</mo><mrow><mn mathvariant="normal">0</mn><mo><</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo><</mo><mi>r</mi></mrow></munder><mrow><msub><mrow><mi>g</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></mrow><msub><mrow><mi>s</mi></mrow><mrow><mi>z</mi></mrow></msub><mfenced separators="|"><mrow><mi>x</mi><mo>-</mo><mi>i</mi><mo>,</mo><mi>y</mi><mo>-</mo><mi>j</mi></mrow></mfenced></mrow><mrow><munder><mo movablelimits="false">∑</mo><mrow><mn mathvariant="normal">0</mn><mo><</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo><</mo><mi>r</mi></mrow></munder><mrow><msub><mrow><mi>g</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></mrow></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332967&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332965&type= 55.37200165 15.32466698 . <math id="M17"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332972&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332970&type= 1.69333339 2.79399991 波段第 <math id="M18"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332977&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332975&type= 1.86266661 2.79399991 个像素点的谱带内预测误差 <math id="M19"><msub><mrow><mi>d</mi></mrow><mrow><mi>z</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332982&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332980&type= 8.38199997 4.06400013 ： <math display="block" id="M20"><msub><mrow><mi>d</mi></mrow><mrow><mi>z</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>=</mo><msub><mrow><mi>s</mi></mrow><mrow><mi>z</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>-</mo><msub><mrow><mover accent="true"><mi>s</mi><mo>̂</mo></mover></mrow><mrow><mi>z</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332986&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332984&type= 33.86666870 4.06400013 . 对于首个波段的预测误差 <math id="M21"><msub><mrow><mi>d</mi></mrow><mrow><mi>z</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332991&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332989&type= 8.38199997 4.06400013 直接使用编码器压缩编码，其余波段则使用格型滤波器组去除谱间冗余信息后再编码。 2.2　谱带间预测 2 高光谱图像在光谱维度上具有较强的相关性，但受噪声影响，各相邻波段谱间相关性的强弱不同。因此，RLS滤波器的预测效果与预测长度（待测像素点滤波器输入向量的长度）并不总是呈现正相关的关系，对于部分谱带，较短的预测长度可以获得良好的预测效果 ［ 17 17 ， 21 21 ］ ，合理地为每个谱带的滤波器选择预测长度，是提升压缩算法性能的关键。 2 . 2 . 1　格型滤波器组 3 对于波段按行交叉格式的高光谱图像，在谱带间预测过程中， <math id="M22"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332995&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332993&type= 1.69333339 2.79399991 波段第 <math id="M23"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333000&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16332998&type= 1.86266661 2.79399991 个像素点共使用 <math id="M24"><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333004&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333002&type= 8.38199997 2.79399991 个不同预测长度的滤波器 <math id="M25"><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo stretchy="false">(</mo><mi>i</mi><mo>=</mo><mo stretchy="false">{</mo><mn mathvariant="normal">1,2</mn><mo>,</mo><mo>⋯</mo><mo>,</mo><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo stretchy="false">}</mo><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333008&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333006&type= 45.29666519 4.06400013 ，其中下标“ <math id="M26"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333013&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333011&type= 1.10066664 2.79399991 ”被用于标识不同预测长度的滤波器。则 <math id="M27"><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333017&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333015&type= 10.32933331 4.06400013 的输入向量 <math id="M28"><msub><mrow><mi mathvariant="bold-italic">x</mi></mrow><mrow><mi mathvariant="normal">z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>=</mo><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><msub><mrow><mi>d</mi></mrow><mrow><mi>z</mi><mo>-</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>,</mo><mo>⋯</mo><mo>,</mo><msub><mrow><mi>d</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></mrow></mfenced></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333022&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333019&type= 51.90066528 5.58799982 ，对应的权重向量 <math id="M29"><msub><mrow><mi mathvariant="bold-italic">w</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>=</mo><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><msub><mrow><mi>w</mi></mrow><mrow><mi mathvariant="normal">z</mi><mo>-</mo><mi mathvariant="normal">i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>,</mo><mo>⋯</mo><mo>,</mo><msub><mrow><mi>w</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></mrow></mfenced></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333026&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333024&type= 54.27133560 5.58799982 ，期望信号为 <math id="M30"><msub><mrow><mi>d</mi></mrow><mrow><mi>z</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333031&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333028&type= 8.38199997 4.06400013 。则 <math id="M31"><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333035&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333033&type= 10.32933331 4.06400013 的计算过程如公式（5）~ 公式（13） 公式（13） 所示。 首先计算 <math id="M32"><msub><mrow><mi mathvariant="normal">F</mi></mrow><mrow><mi mathvariant="normal">z</mi><mo>,</mo><mi mathvariant="normal">i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333040&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333038&type= 10.07533360 4.06400013 的谱带间预测残差： <math display="block" id="M33"><msub><mrow><mi>e</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>=</mo><msub><mrow><mi>d</mi></mrow><mrow><mi>z</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">x</mi></mrow><mrow><mi mathvariant="normal">z</mi><mo>,</mo><mi mathvariant="normal">i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><msub><mrow><mi mathvariant="bold-italic">w</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333045&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333042&type= 56.04933548 4.06400013 . 预测误差方差的累计值 <math id="M34"><msub><mrow><mi>v</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333050&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333047&type= 9.65199947 4.06400013 ： <math display="block" id="M35"><msub><mrow><mi>v</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>=</mo><msub><mrow><msup><mrow><mi>e</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>+</mo><msub><mrow><mi>v</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333055&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333052&type= 47.24399948 4.06400013 . 设： <math display="block" id="M36"><mtable columnalign="left"><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">c</mi></mrow><mrow><mi mathvariant="bold-italic">z</mi><mo>,</mo><mi mathvariant="bold-italic">i</mi></mrow></msub><mfenced separators="|"><mrow><mi mathvariant="bold-italic">n</mi></mrow></mfenced><mo>=</mo></mtd></mtr><mtr><mtd><mfenced separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi mathvariant="bold-italic">c</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi mathvariant="normal">i</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">w</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi mathvariant="normal">i</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><msub><mrow><mi>e</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><msubsup><mrow><mi>v</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></mtd></mtr><mtr><mtd><msub><mrow><mi>e</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><msubsup><mrow><mi>v</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></mtd></mtr></mtable></mrow></mfenced></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333059&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333057&type= 62.73799896 16.93333435 ， <math display="block" id="M37"><msub><mrow><mi>u</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>=</mo><msub><mrow><mi>u</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi mathvariant="normal">i</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>+</mo><msubsup><mrow><mi>e</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><msubsup><mrow><mi>v</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msubsup><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333063&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333061&type= 56.98066711 4.06400013 ， 其中： <math id="M38"><msub><mrow><mi mathvariant="bold-italic">c</mi></mrow><mrow><mi mathvariant="bold-italic">z</mi><mo>,</mo><mi mathvariant="bold-italic">i</mi></mrow></msub><mfenced separators="|"><mrow><mi mathvariant="bold-italic">n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333070&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333067&type= 9.65199947 4.06400013 ， <math id="M39"><msub><mrow><mi>u</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333074&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333072&type= 9.65199947 4.06400013 分别为滤波器的投影向量、投影值。 滤波器的增益为： <math display="block" id="M40"><msub><mrow><mi>k</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mo stretchy="false">(</mo><mi>n</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>c</mi></mrow><mrow><mi mathvariant="normal">z</mi><mo>,</mo><mi mathvariant="normal">i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><msup><mrow><mfenced separators="|"><mrow><mi>λ</mi><mo>+</mo><msub><mrow><mi>u</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></mrow></mfenced></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333079&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333077&type= 46.48200226 5.58799982 ， 则权重向量 <math id="M41"><msub><mrow><mi mathvariant="bold-italic">w</mi></mrow><mrow><mi mathvariant="bold-italic">z</mi><mo>,</mo><mi mathvariant="bold-italic">i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333083&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333081&type= 10.92199993 4.06400013 的递归更新如 式（10） 式（10） 所示： <math display="block" id="M42"><msub><mrow><mi mathvariant="bold-italic">w</mi></mrow><mrow><mi mathvariant="bold-italic">z</mi><mo>,</mo><mi mathvariant="bold-italic">i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">w</mi></mrow><mrow><mi mathvariant="bold-italic">z</mi><mo>,</mo><mi mathvariant="bold-italic">i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced><mo>+</mo><msub><mrow><mi>k</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mo stretchy="false">(</mo><mi>n</mi><mo stretchy="false">)</mo><msub><mrow><mi>e</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333088&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333086&type= 58.84333038 4.06400013 . 为确保滤波器的收敛性及屏蔽一些异常的像素值，当滤波器阶数≤3时，投影向量 <math id="M43"><msub><mrow><mi mathvariant="bold-italic">c</mi></mrow><mrow><mi mathvariant="bold-italic">z</mi><mo>,</mo><mi mathvariant="bold-italic">i</mi></mrow></msub><mfenced separators="|"><mrow><mi mathvariant="bold-italic">n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333092&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333090&type= 9.65199947 4.06400013 及投影值 <math id="M44"><msub><mrow><mi>u</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333097&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333095&type= 9.65199947 4.06400013 采用传统RLS滤波方式更新，如式（11）~ 式（13） 式（13） 所示。 <math display="block" id="M45"><msub><mrow><mi mathvariant="bold-italic">c</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="bold-italic">z</mi><mo>,</mo><mi mathvariant="bold-italic">i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced><msubsup><mrow><mi mathvariant="bold-italic">x</mi></mrow><mrow><mi mathvariant="normal">z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>i</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msubsup><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333105&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333102&type= 49.27600098 4.06400013 ， <math display="block" id="M46"><msub><mrow><mi>u</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">x</mi></mrow><mrow><mi mathvariant="normal">z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>i</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><msub><mrow><mi mathvariant="bold-italic">c</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi mathvariant="normal">i</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333111&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333107&type= 48.76799774 4.06400013 ， <math display="block" id="M47"><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="bold-italic">z</mi><mo>,</mo><mi mathvariant="bold-italic">i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>λ</mi></mrow></mfrac><mfenced open="[" close="]" separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="bold-italic">z</mi><mo>,</mo><mi mathvariant="bold-italic">i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced><mo>-</mo><msub><mrow><mi>k</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mo stretchy="false">(</mo><mi>n</mi><mo stretchy="false">)</mo><msubsup><mrow><mi mathvariant="bold-italic">c</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msubsup><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333117&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333113&type= 63.92333221 7.53533316 ， 其中：根据文献［ 21 21 ］RLS滤波器敏感性测试结果初始化遗忘因子为 <math id="M48"><mi>λ</mi><mo>=</mo><mn mathvariant="normal">0.999</mn><mtext> </mtext><mn mathvariant="normal">5</mn><mo>;</mo><msub><mrow><mi mathvariant="bold-italic">w</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mn mathvariant="normal">0</mn></mrow></mfenced><mo>=</mo><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><mn mathvariant="normal">0,0</mn><mo>,</mo><mo>…</mo><mo>,</mo><mn mathvariant="normal">0</mn></mrow></mfenced></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333123&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333120&type= 56.13400269 4.74133348 。依据上述更新过程可知： <math id="M49"><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>-</mo><mi>i</mi><mo>+</mo><mn mathvariant="normal">1,1</mn></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>,</mo><mo>⋯</mo><mo>,</mo><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333134&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333130&type= 35.72933197 4.06400013 是依次递归更新的，这组具有前向迭代更新关系的滤波器称为链式序列。因高光谱图像共有 <math id="M50"><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333141&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333137&type= 2.87866688 2.87866688 个波段，首个波段的谱带内预测误差直接使用编码器压缩编码，本文提出的格型滤波器组中共包含 <math id="M51"><mi>N</mi><mo>-</mo><mn mathvariant="normal">1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333148&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333145&type= 9.90599918 2.87866688 条链式序列，如 图4 图4 所示，为简洁表示，省略每个滤波器的部分符号 <math id="M52"><mo stretchy="false">(</mo><mi>n</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333153&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333150&type= 4.99533367 3.55599999 。图中每一列表示同一待测点不同预测长度的滤波器；箭头表示同一条链式序列内滤波器的迭代更新方向，将用于迭代更新 <math id="M53"><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333160&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333157&type= 10.32933331 4.06400013 的链式序列记作集合 <math id="M54"><msub><mrow><mi mathvariant="bold">Ω</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333166&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333163&type= 10.83733368 4.06400013 ： <math display="block" id="M55"><msub><mrow><mi mathvariant="bold">Ω</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>=</mo><mfenced open="{" close="}" separators="|"><mrow><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>-</mo><mi>i</mi><mo>+</mo><mn mathvariant="normal">1,1</mn></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>,</mo><mo>⋯</mo><mo>,</mo><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>i</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333172&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333170&type= 62.14533615 4.91066647 . 格型滤波器组在预测过程中， <math id="M56"><msub><mrow><mi>v</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333195&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333190&type= 9.65199947 4.06400013 被用于衡量同一待测点不同预测长度滤波器的预测精度。预测误差被视为服从均值为0的高斯分布，其信息熵 <math id="M57"><mi>H</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333203&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333199&type= 2.96333337 2.79399991 与方差 <math id="M58"><msup><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333212&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333207&type= 3.04799986 3.04799986 呈现正相关关系，如 式（15） 式（15） 所示： <math display="block" id="M59"><mi>H</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msup><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced><mo>+</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333221&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333216&type= 33.02000046 7.45066643 . 因此，在波段 <math id="M60"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333232&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333227&type= 1.69333339 2.79399991 第 <math id="M61"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333245&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333238&type= 1.86266661 2.79399991 个像素点的不同预测长度滤波器中，具有最小预测误差方差累计值 <math id="M62"><msub><mrow><mi>v</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333264&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333255&type= 12.02266693 4.06400013 的滤波器称为最优滤波器 <math id="M63"><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi mathvariant="normal">o</mi><mi mathvariant="normal">p</mi><mi mathvariant="normal">t</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333274&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333269&type= 12.27666664 4.06400013 ，其中 <math id="M64"><msub><mrow><mi>v</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>=</mo><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">{</mo><msub><mrow><mi>v</mi></mrow><mrow><mi>z</mi><mo>,</mo><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><mo>⋯</mo><mo>,</mo><msub><mrow><mi>v</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo stretchy="false">}</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333288&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333280&type= 54.18666840 4.06400013 。最优滤波器的预测长度被称为最优预测长度，记作 <math id="M65"><msub><mrow><mi>L</mi></mrow><mrow><mi>z</mi></mrow></msub><mo stretchy="false">(</mo><mi>n</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333298&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333294&type= 9.14400005 4.48733330 。 若直接使用格型滤波器组顺次计算每个待测点不同预测长度的滤波器，虽然能获得最佳的预测效果，但耗费过长的计算时间。在此，通过对滤波器收敛性的考察，优化了格型滤波器组的更新方式，为衡量 <math id="M66"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333307&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333304&type= 1.69333339 2.79399991 波段第 <math id="M67"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333315&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333311&type= 1.86266661 2.79399991 个与第 <math id="M68"><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333325&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333320&type= 8.55133343 2.79399991 个待测点权值的变化情况，定义相邻滤波器权值的距离（Distance of Weights， DOW） <math display="block" id="M69"><msub><mrow><mi>D</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>=</mo><mfenced open="‖" close="‖" separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">w</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">w</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333336&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333332&type= 53.67866516 4.91066647 . 以AVIRIS 2006高光谱图像集中YellowStone 0校准图像为例，滤波器 <math id="M70"><msub><mrow><mi>F</mi></mrow><mrow><mn mathvariant="normal">223,222</mn></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333345&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333341&type= 14.73199940 4.06400013 的DOW随迭代次数的变化趋势如 图5 图5 所示，随迭代次数的增加，各个滤波器具有明显的收敛趋势。因此，同一谱带不同预测长度滤波器参数 <math id="M71"><msub><mrow><mi>v</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333356&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333350&type= 9.65199947 4.06400013 的排布次序趋于稳定，最优滤波器搜索阶段可以被提前终止，后续谱带的最优滤波器被认为是相同的。 引入限定迭代次数的阈值 <math id="M72"><mi>T</mi><mo>,</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333381&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333376&type= 4.40266657 3.21733332 当 <math id="M73"><mi>n</mi><mo>=</mo><mi>T</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333390&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333386&type= 9.22866726 2.79399991 时，先计算波段 <math id="M74"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333401&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333394&type= 1.69333339 2.79399991 第 <math id="M75"><mi>T</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333412&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333407&type= 2.53999996 2.79399991 个待测点的最优滤波器 <math id="M76"><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>o</mi><mi>p</mi><mi>t</mi></mrow></msub><mo stretchy="false">(</mo><mi>T</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333421&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333417&type= 13.54666710 4.48733330 ，随后更新其链式序列集合 <math id="M77"><msub><mrow><mi mathvariant="bold-italic">Ω</mi></mrow><mrow><mi mathvariant="bold-italic">z</mi><mo>,</mo><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">t</mi></mrow></msub><mfenced separators="|"><mrow><mi>T</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333432&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333425&type= 13.46199894 4.06400013 。为统计当前谱带所有待测点最优滤波器的链式序列，使用最优滤波器集合 <math id="M78"><mi mathvariant="bold-italic">Ω</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333442&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333435&type= 2.96333337 2.79399991 表示每个待测点 <math id="M79"><msub><mrow><mi mathvariant="bold-italic">Ω</mi></mrow><mrow><mi mathvariant="bold-italic">z</mi><mo>,</mo><mi mathvariant="bold-italic">o</mi><mi mathvariant="bold-italic">p</mi><mi mathvariant="bold-italic">t</mi></mrow></msub><mfenced separators="|"><mrow><mi>T</mi></mrow></mfenced><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo><</mo><mi>z</mi><mo>≤</mo><mi>N</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333457&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333451&type= 33.86666870 4.06400013 的并集。当 <math id="M80"><mi>n</mi><mo>></mo><mi>T</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333468&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333463&type= 9.22866726 2.79399991 时，最优谱带搜索阶段将被终止，进入快速预测阶段，只有位于 <math id="M81"><mi mathvariant="bold-italic">Ω</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333480&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333473&type= 2.96333337 2.79399991 中的滤波器需要被持续更新。为直观地描述格型滤波器组的预测过程，算法流程图如 图6 图6 所示。 在上述格型滤波器组的更新过程中，当 <math id="M82"><mi>n</mi><mo>></mo><mi>T</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333504&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333497&type= 9.22866726 2.79399991 时，位于 <math id="M83"><mi mathvariant="bold">Ω</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333513&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333508&type= 3.04799986 2.79399991 中的滤波器需要被持续更新。因此，越多的滤波器位于同一条链式序列中，待更新滤波器的数量越少，消耗的计算时间越短。 2 . 2 . 2　提升策略 3 以校准后YellowStone 0图像的400 th 和500 th LSP为例，各谱带的最优预测长度如 图7 图7 所示，除个别受噪声影响的波段外，大多数谱带的最优预测长度随 <math id="M84"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333521&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333517&type= 1.69333339 2.79399991 的增大而增大。因此，可以经筛选分类，使更多的待更新滤波器处于同一链式序列。本文提出最长更新规则：若多个滤波器的最优预测长度相近时，应将具有最长链式序列的滤波器视为最优滤波器。为简化分类操作，若 <math id="M85"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333533&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333527&type= 1.69333339 2.79399991 波段第 <math id="M86"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333540&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333536&type= 1.86266661 2.79399991 个待测点的最优滤波器预测长度 <math id="M87"><msub><mrow><mi>L</mi></mrow><mrow><mi>z</mi></mrow></msub><mo stretchy="false">(</mo><mi>n</mi><mo stretchy="false">)</mo><mo>≥</mo><mfrac><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333547&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333543&type= 17.10266685 7.53533316 时，则将最长链式序列中的滤波器视为最优滤波器；否则 <math id="M88"><msub><mrow><mi>v</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333554&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333551&type= 9.65199947 4.06400013 最小的滤波器被用于迭代计算。 此外，在文献［ 21 21 ］中，高放等人发现对于同一成像系统拍摄的高光谱图像集，其最优预测长度高度一致，并建议将首个图像的最优预测长度直接应用于预测该场景的其余图像。该加速方案同样应用于本文算法，记为RLS Speed Lattice Filter Group （SLFG-RLS）。 2.3　算术编码 2 本文采用算术编码器做压缩编码。因最终预测误差数据是用16位表示的，故算术编码的码表需要使用65536个符号。根据预测残差的概率分布，大多数符号没有被使用，故采用自适应码表。初始的码表包含0和ESC两个符号。当需要对一个新符号编码时，编码器先使用ESC的概率对其编码，随后再更新码表便于后续编码。 3 实验结果与分析 1 本文压缩方案中的所有程序均在实验平台（Intel i7-7700k CPU 4.20GHZ/16GBRAM）中以C语言进行编写、测试。CCSDS的多光谱和高光谱数据压缩工作组提供的AVIRIS 2006高光谱图像数据集作为本文算法的测试数据源。 测试图像集包括5个16位校准和5个16位未校准的场景，每个场景包含512行，224个谱带的数据，其中校准与非校准场景每行分别有677和680个像素点。在此，将场景名称的首字母及其序列号用作缩写，例如：Yellowstone scene 0记作YS_SC0。 3.1　高斯预测器的参数选择 2 在谱带内预测阶段，为便于计算，本文采用整数型的单边高斯预测器。本文将均值滤波器及不同参数的高斯预测器做测试与比较。以YS_SC0场景为例，实验结果如 表1 表1 所示。 滤波器参数 谱带内预测计算时间/ （ms·band -1 ） 压缩结果/ (bits·pixel -1 ) 2 均值滤波器 r =3 26 3.55 r =4 37 3.55 8 单边高斯 滤波器 r =3， σ =1 39 3.57 r =3， σ =2 39 3.54 r =4， σ =1 54 3.57 r =4， σ =2 54 3.53 r =4， σ =3 54 3.53 r =5， σ =1 85 3.57 r =5， σ =2 85 3.53 r =5， σ =3 85 3.54 在窗口半径相同时，对于每个谱带高斯预测器比均值预测器多消耗十余毫秒的计算时间，但高斯预测器能获得更佳的压缩效果。对于高斯预测器，当 <math id="M89"><mi>σ</mi><mo>≥</mo><mn mathvariant="normal">2</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333583&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333580&type= 8.55133343 3.30199981 时，能获得较好的压缩效果；当窗口半径 <math id="M90"><mi>r</mi><mo>≥</mo><mn mathvariant="normal">4</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333589&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333586&type= 8.38199997 3.30199981 时，预测器的压缩效率没有明显提高，但时间复杂度却骤增。考虑到算法的计算复杂度和实时性，最终选用 <math id="M91"><mi>σ</mi><mo>=</mo><mn mathvariant="normal">2</mn><mo>,</mo><mi>r</mi><mo>=</mo><mn mathvariant="normal">4</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333598&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333594&type= 18.71133423 3.21733332 作为单边高斯预测器的初始化参数。 3.2　迭代阈值的选择 2 随迭代次数的增大，滤波器的最优预测长度趋于稳定。以校准后YS_SC0场景中209，140，94，49波段的最优滤波器为例， 图8 图8 展示不同LSP的最优预测长度 <math id="M92"><msub><mrow><mi>L</mi></mrow><mrow><mi>z</mi></mrow></msub><mo stretchy="false">(</mo><mi>n</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333606&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333602&type= 9.14400005 4.48733330 。 当迭代次数的阈值 <math id="M93"><mi>T</mi><mo>≥</mo><mn mathvariant="normal">50</mn><mtext>  </mtext><mi mathvariant="normal">L</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">P</mi><mi mathvariant="normal">s</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333626&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333620&type= 21.33599854 3.30199981 时，最优预测长度波动很小。根据 公式（15） 公式（15） 和 公式（16） 公式（16） ，随滤波器迭代训练次数的增加，收敛趋势愈加明显，同一谱带不同预测长度滤波器参数 <math id="M94"><msub><mrow><mi>v</mi></mrow><mrow><mi mathvariant="normal">z</mi><mo>,</mo><mi mathvariant="normal">i</mi></mrow></msub><mfenced separators="|"><mrow><mi>n</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333632&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333630&type= 9.39799976 4.06400013 的排布次序及滤波器最优预测长度波动渐小，趋于稳定。为量化描述迭代阈值对算法的压缩结果、计算时间的影响，以校准后YS_SC0场景为例， 表2 表2 展现格型滤波器组在使用提升策略前的压缩结果、计算时间与 <math id="M95"><mi>T</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333637&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333635&type= 2.53999996 2.79399991 的关系。 迭代阈值 /LSPs 压缩结果/ (bits·pixel -1 ) 计算 时间/s 512 3.53 8 016.18 300 3.53 5 364.46 100 3.53 3 209.91 50 3.53 2 468.32 30 3.54 2 361.27 10 3.55 1 913.73 5 3.57 1 550.26 1 3.64 729.98 从 表2 表2 中可知，当 <math id="M96"><mi>T</mi><mo>≥</mo><mn mathvariant="normal">50</mn><mtext>  </mtext><mi mathvariant="normal">L</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">P</mi><mi mathvariant="normal">s</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333646&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333644&type= 21.33599854 3.30199981 时，算法的压缩结果趋于固定，但计算时间急剧上升。考虑到实际压缩结果和时间复杂度的双重影响，本文算法将参数 <math id="M97"><mi>T</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333650&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333648&type= 2.53999996 2.79399991 初始化为 <math id="M98"><mn mathvariant="normal">50</mn><mtext>  </mtext><mi mathvariant="normal">L</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">P</mi><mi mathvariant="normal">s</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333654&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333653&type= 13.63133430 2.79399991 ，从而保证良好的压缩效果以及较低的时间复杂度。 3.3　压缩性能的评估与分析 2 为衡量本文算法的压缩性能，将其与现有的算法（FL， IP3， RLS， ACRLS， and Fast-RLS-ALP）从压缩结果与计算时间两方面作对比。 表3 表3 展示各类算法的压缩结果。其中，ACRLS算法因采用穷举法查找每个波段的最佳预测长度，是已知同类型压缩结果最佳的压缩算法。本文提出的LFG-RLS算法采用格型滤波器组获取图像谱带的最优预测长度，并结合滤波器的收敛性，合理地选择迭代阈值，获取良好的压缩效果。在校准、未校准场景中的平均压缩结果分别是是3.34 bits/pixel和5.61 bits/pixel。该压缩效果与ACRLS非常接近，并优于其他压缩方案。此外，SLFG-RLS的压缩结果证实，当首个场景的预测长度应用于其他场景时，能获得相似的压缩效果。 2 图像名称 FL IP3 RLS Fast-RLS-ALP ACRLS LFG-RLS SLFG-RLS 6 校 准 图 像 YS_SC0 3.87 3.84 3.77 3.62 3.50 3.53 3.53 YS_SC3 3.75 3.74 3.63 3.49 3.39 3.42 3.41 YS_SC10 3.35 3.21 3.24 3.16 3.01 3.07 3.07 YS_SC11 3.56 3.51 3.48 3.34 3.24 3.25 3.26 YS_C18 3.83 3.75 3.67 3.52 3.41 3.43 3.44 平均值 3.67 3.61 3.56 3.43 3.31 3.34 3.34 6 未 校 准 图 像 YS_SC0 6.23 6.08 6.01 5.88 5.79 5.79 5.79 YS_SC3 6.10 5.96 5.86 5.75 5.67 5.70 5.69 YS_SC10 5.57 5.44 5.46 5.40 5.28 5.33 5.33 YS_SC11 5.84 5.74 5.71 5.57 5.50 5.51 5.52 YS_C18 6.29 6.00 5.90 5.78 5.70 5.72 5.72 平均值 6.00 5.84 5.79 5.68 5.59 5.61 5.61 各算法的平均计算时间如 表4 表4 所示，FL，IP3，RLS和Fast-RLS-ALP算法采用固定预测长度，预测方式较为简单，计算时间较短。ACRLS算法采用穷尽式的搜索方案计算所有待测点的最优预测长度，消耗较长的计算时间。 图像名称 FL IP3 RLS Fast - RLS - ALP ACRLS LFG-RLS SLFG-RLS 校准图像 62 121 81 110 11 231 1 050 441 未校准图像 61 121 81 110 11 843 1 064 447 s 本文提出的LFG-RLS算法，对于校准、未校准图像分别仅需1 050 s，1 064 s，相较于使用穷举法搜索最优谱带预测长度的ACRLS算法降低10余倍的运算时间。主要原因是本文合理地限定滤波器的迭代次数，并根据格型滤波器组链式更新的特点，对最优滤波器做分类与简化，大幅度减少待更新滤波器的数量，降低快速预测阶段的计算时间。 此外，谱带间预测器根据迭代阈值 <math id="M99"><mi>T</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333667&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16333664&type= 2.53999996 2.79399991 可分为最优谱带搜索阶段和快速预测阶段，对于校准图像，LFG-RLS压缩方案两阶段平均耗时分别为767.8 s，262.9 s，最优谱带搜索阶段压缩的数据量约占总数据量的10%，但消耗总计算时间的73.1%。因此，若高光谱成像仪采集的图像尺寸越长或图像数量越多，最优谱带搜索阶段消耗计算时间所占比例越小，压缩方案所需总运算时间相对越短，能充分发挥算法的优越性。在SLFG-RLS压缩方案中，系列校准、未校准图像分别是同一系统采集的多幅图像，将首个场景的预测结果推广后，对校准、未校准图像的平均计算时间分别降低至441 s和447 s。 为验证本文压缩算法的无损性，将压缩后的数据流做解压缩处理，得到恢复图像。将原始图像与恢复图像相比较，结果显示二者完全相同，表明高光谱图像信息没有任何损失。 4 结　论 1 为充分挖掘高光谱图像的空间相关性和谱间相关性，本文提出利用格型递归最小二乘法滤波器组的高光谱图像无损压缩方案。该方案中充分考虑最优滤波器对预测效果的影响，采用格型滤波器组计算每个谱带的最优滤波器，并且根据其链式更新的特点，对筛选过程做进一步的分类与简化，大幅度降低压缩方案的计算复杂度。实验结果表明，以CCSDS提供的AVIRIS高光谱图像为测试数据集，本文的压缩方案具有良好的压缩性能，对16位校准图像、未校准图像分别取得3.34 bits/pixel和5.61 bits/pixel的压缩结果，与同类算法的最佳压缩结果相似，且其平均计算时间远低于其余具有相似压缩结果的算法，具有较强的竞争力，为高光谱图像无损压缩提供了一种有效的解决方案。未来将针对压缩方案的并行化处理做研究与探讨，并开发可行的硬件实现方案。 参考文献 陈善学 ， 张燕琪 . 基于自适应波段聚类主成分分析和反向传播神经网络的高光谱图像压缩 ［J］. 电子与信息学报 ， 2018 ， 40 （ 10 ）： 2478 - 2483 . 基于自适应波段聚类主成分分析和反向传播神经网络的高光谱图像压缩 CHEN SH X ， ZHANG Y Q . Hyperspectral image compression based on adaptive band clustering principal component analysis and back propagation neural network ［J］. Journal of Electronics & Information Technology ， 2018 ， 40 （ 10 ）： 2478 - 2483 . （in Chinese） Hyperspectral image compression based on adaptive band clustering principal component analysis and back propagation neural network 张宁 ， 冯书谊 ， 濮建福 ， 等 . 多光谱遥感图像CCSDS动态码率控制近无损压缩 ［J］. 光学 精密工程 ， 2015 ， 23 （ 6 ）： 1783 - 1790 . 多光谱遥感图像CCSDS动态码率控制近无损压缩 ZHANG N ， FENG SH Y ， PU J F ， et al . Dynamic rate control for CCSDS nearly lossless compression of multispectral remote image ［J］. Opt. Precision Eng. ， 2015 ， 23 （ 6 ）： 1783 - 1790 . （in Chinese） Dynamic rate control for CCSDS nearly lossless compression of multispectral remote image 宋金伟 ， 张忠伟 ， 陈晓敏 . 利用线性预测与查表法的高光谱图像压缩 ［J］. 光学 精密工程 ， 2013 ， 21 （ 8 ）： 2201 - 2208 . 利用线性预测与查表法的高光谱图像压缩 SONG J W ， ZHANG ZJ W ， CHEN X M . Hyperspectral imagery compression via linear prediction and lookup tables ［J］. Opt. Precision Eng. ， 2013 ， 21 （ 8 ）： 2201 - 2208 . （in Chinese） Hyperspectral imagery compression via linear prediction and lookup tables XUE J ， ZHAO Y ， LIAO W ， et al . Hyper-Laplacian regularized nonlocal low-rank matrix recovery for hyperspectral image compressive sensing reconstruction ［J］. Information Sciences ， 2019 ， 501 ： 406 - 420 . Hyper-Laplacian regularized nonlocal low-rank matrix recovery for hyperspectral image compressive sensing reconstruction KONG W ， WU J ， HU Z ， et al . Lossless compression codec of aurora spectral data using hybrid spatial-spectral decorrelation with outlier recognition ［J］. Journal of Visual Communication and Image Representation ， 2019 ， 62 ： 174 - 181 . Lossless compression codec of aurora spectral data using hybrid spatial-spectral decorrelation with outlier recognition NAGENDRAN R ， VASUKI A . Hyperspectral image compression using hybrid transform with different wavelet-based transform coding ［J］. International Journal of Wavelets ， Multiresolution and Information Processing， 2019 ， 18 （ 4 ）. Hyperspectral image compression using hybrid transform with different wavelet-based transform coding LUO J ， WU J ， ZHAO S ， et al . Lossless compression for hyperspectral image using deep recurrent neural networks ［J］. International Journal of Machine Learning and Cybernetics ， 2019 ， 10 （ 10 ）： 2619 - 2629 . Lossless compression for hyperspectral image using deep recurrent neural networks 孙蕾 ， 罗建书 ， 谷德峰 . 基于谱间预测和码流预分配的高光谱图像压缩算法 ［J］. 光学 精密工程 ， 2008 （ 4 ）： 752 - 757 . 基于谱间预测和码流预分配的高光谱图像压缩算法 SUN L ， LUO J SH ， GU D F . Hyperspectral image compression algorithm based on prediction between bands and rate pre-allocation ［J］. Opt. Precision Eng. ， 2008 （ 4 ）： 752 - 757 . （in Chinese） Hyperspectral image compression algorithm based on prediction between bands and rate pre-allocation LI R ， PAN Z ， WANG Y . The linear prediction vector quantization for hyperspectral image compression ［J］. Multimedia Tools and Applications ， 2018 ， 78 （ 9 ）： 11701 - 11718 . The linear prediction vector quantization for hyperspectral image compression BAJPAI S ， KIDWAI N R ， SINGH H V ， et al . Low memory block tree coding for hyperspectral images ［J］. Multimedia Tools and Applications ， 2019 ， 78 （ 19 ）： 27193 - 27209 . Low memory block tree coding for hyperspectral images ALVAREZ-CORTES S ， SERRA-SAGRISTA J ， BARTRINA-RAPESTA J ， et al . Regression wavelet analysis for near-lossless remote sensing data compression ［J］. IEEE Transactions on Geoscience and Remote Sensing ， 2020 ， 58 （ 2 ）： 790 - 798 . Regression wavelet analysis for near-lossless remote sensing data compression RYAN M J ， AJ F . The lossless compression of AVIRIS images by vector quantization ［J］. IEEE Transactions on Geoscience & Remote Sensing ， 1997 ， 35 ： 546 - 550 . The lossless compression of AVIRIS images by vector quantization 高放 ， 孙长建 ， 邵庆龙 ， 等 . 基于K-均值聚类和传统递归最小二乘法的高光谱图像无损压缩 ［J］. 电子与信息学报 ， 2016 ， 38 （ 11 ）： 2709 - 2714 . 基于K-均值聚类和传统递归最小二乘法的高光谱图像无损压缩 GAO F ， SUN CH J ， SHAO Q L ， et al . Lossless compression of hyperspectral images using k-means clustering and conventional recursive least-squares predictor ［J］. Opt. Precision Eng. ， 2016 ， 38 （ 11 ）： 2709 - 2714 . （in Chinese） Lossless compression of hyperspectral images using k-means clustering and conventional recursive least-squares predictor MIELIKAINEN J . Lossless compression of hyperspectral images using lookup tables ［J］. IEEE Signal Processing Letters ， 2006 ， 13 （ 3 ）： 157 - 160 . Lossless compression of hyperspectral images using lookup tables BORMIN H ， SRIRAJA Y . Lossless compression of hyperspectral imagery via lookup tables with predictor selection ［J］. Proceedings of the SPIE-The International Society for Optical Engineering ， 2006 ， 6365 . Lossless compression of hyperspectral imagery via lookup tables with predictor selection MIELIKAINEN J ， TOIVANEN P . Clustered DPCM for the lossless compression of hyperspectral images ［J］. IEEE Transactions on Geoscience and Remote Sensing ， 2003 ， 41 （ 12 ）： 2943 - 2946 . Clustered DPCM for the lossless compression of hyperspectral images MIELIKAINEN J ， HUANG B . Lossless compression of hyperspectral images using clustered linear prediction with adaptive prediction length ［J］. IEEE Geoscience and Remote Sensing Letters ， 2012 ， 9 （ 6 ）： 1118 - 1121 . Lossless compression of hyperspectral images using clustered linear prediction with adaptive prediction length CCSDS . Lossless Multispectral & Hyperspectral Image Compression Recommended Standard ［M］. USA ： CCSDS ， 2012 . Lossless Multispectral & Hyperspectral Image Compression Recommended Standard CHENG-CHEN L ， YIN-TSUNG H . An efficient lossless compression scheme for hyperspectral images using two-stage prediction ［J］. IEEE Geoscience and Remote Sensing Letters ， 2010 ， 7 （ 3 ）. An efficient lossless compression scheme for hyperspectral images using two-stage prediction SONG J ， ZHANG Z ， CHEN X . Lossless compression of hyperspectral imagery via RLS filter ［J］. Electronics Letters ， 2013 ， 49 （ 16 ）： 992 - 994 . Lossless compression of hyperspectral imagery via RLS filter GAO F ， GUO S . Lossless compression of hyperspectral images using conventional recursive least-squares predictor with adaptive prediction bands ［J］. Journal of Applied Remote Sensing ， 2016 ， 10 （ 1 ）. Lossless compression of hyperspectral images using conventional recursive least-squares predictor with adaptive prediction bands SONG J ， ZHOU L ， DENG C ， et al . Lossless compression of hyperspectral imagery using a fast adaptive-length-prediction RLS filter ［J］. Remote Sensing Letters ， 2019 ， 10 （ 4 ）： 401 - 410 . Lossless compression of hyperspectral imagery using a fast adaptive-length-prediction RLS filter 10.37188/OPE.20212904.0896 TP751 A 1004-924X（2021）04-0896-10 郑铁,薛长斌,宋金伟.利用格型递归最小二乘滤波器组的高光谱图像压缩[J].光学精密工程,2021,29(04):896-905. 郑铁,薛长斌,宋金伟.利用格型递归最小二乘滤波器组的高光谱图像压缩[J].光学精密工程,2021,29(04):896-905. ZHENG Tie,XUE Chang-bin,SONG Jin-wei.Lossless compression of hyperspectral images using recursive least square lattice filter group[J].Optics and Precision Engineering,2021,29(04):896-905. ZHENG Tie,XUE Chang-bin,SONG Jin-wei.Lossless compression of hyperspectral images using recursive least square lattice filter group[J].Optics and Precision Engineering,2021,29(04):896-905. 中国科学院基金资助项目（No. Y87015A080） 2020-06-29 2020-08-17 2021-04-15 2021-06-03 光学精密工程 4 29