FPGA implementation of CCSDS compressor for snapshot mosaic hyperspectral images

快照马赛克高光谱图像CCSDS压缩器的FPGA实现 FPGA implementation of CCSDS compressor for snapshot mosaic hyperspectral images 1 first-author 吕德远（1995-），男，山东聊城人，硕士研究生，2018年于长春理工大学获得学士学位，现为深圳大学物理与光电工程学院硕士研究生，主要从事嵌入式系统、高光谱图像处理等方面的研究。E-mail：lvdeyuan@outlook.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051292&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051295&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051329&type= 22.01332855 32.00399780 吕德远 （1995-），男，山东聊城人，硕士研究生，2018年于长春理工大学获得学士学位，现为深圳大学物理与光电工程学院硕士研究生，主要从事嵌入式系统、高光谱图像处理等方面的研究。E-mail： lvdeyuan@outlook.com lvdeyuan@outlook.com 2 2 2 1 2 corresp E-mail：zhaozhigang@sztu.edu.cn E-mail：zhaozhigang@sztu.edu.cn 赵志刚（1977-），男，湖北黄冈人，博士、副教授、硕士生导师，1999年、2005年于军械工程学院（现陆军工程大学）分别获得学士、硕士学位，2010年于华中科技大学获得博士学位，现为深圳技术大学新材料与新能源学院副教授，主要从事微型高光谱成像、实时高光谱物质分类等方面的研究。E-mail： zhaozhigang@sztu.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051297&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051337&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051334&type= 22.01332855 32.00400162 赵志刚 （1977-），男，湖北黄冈人，博士、副教授、硕士生导师，1999年、2005年于军械工程学院（现陆军工程大学）分别获得学士、硕士学位，2010年于华中科技大学获得博士学位，现为深圳技术大学新材料与新能源学院副教授，主要从事微型高光谱成像、实时高光谱物质分类等方面的研究。E-mail： zhaozhigang@sztu.edu.cn zhaozhigang@sztu.edu.cn 深圳技术大学 新材料与新能源学院，广东 深圳 518118 深圳技术大学 新材料与新能源学院，广东 深圳 518118 College of New Materials and New Energies， Shenzhen Technology University， Shenzhen 518118， China College of New Materials and New Energies， Shenzhen Technology University， Shenzhen 518118， China 深圳大学 物理与光电工程学院，广东 深圳 518052 深圳大学 物理与光电工程学院，广东 深圳 518052 College of Physics and Optoelectronic Engineering， Shenzhen University， Shenzhen 518052， China College of Physics and Optoelectronic Engineering， Shenzhen University， Shenzhen 518052， China 为解决快照马赛克高光谱图像数据量大、对传输带宽要求高的问题，采用最新的CCSDS 123.0-B-2多光谱/高光谱无损和近无损压缩国际标准，实现了基于FPGA的快照马赛克高光谱图像无损和近无损压缩。通过改进BIP（Band-Interleaved Pixels）排序算法，使每个时钟周期能够动态处理一个快照马赛克高光谱样本，数据吞吐率约提高至11倍，有效解决了该国际标准在快照马赛克高光谱图像压缩FPGA硬件实现时存在的流水线和并行化处理困难、处理速度慢的问题。实验结果表明，所实现的压缩器在Xilinx XC7Z020CLG-4平台上布局布线后FPGA整体逻辑资源占用率小于10%，100 MHz系统时钟下约22 ms即可完成一张高光谱图像的压缩，无损压缩性能在2.66~4.30 bits/samples之间，近无损压缩性能在1.01~3.70 bits/samples之间，能够满足快照马赛克高光谱成像技术在无线手持及无人机载领域的应用需求。 To solve the problem of large data volume and high transmission bandwidth requirement of snapshot mosaic hyperspectral images， the latest CCSDS 123.0-B-2 multispectral/hyperspectral lossless and near-lossless compression international standard is adopted to realize the FPGA-based lossless and near-lossless compression of snapshot mosaic hyperspectral images. By improving the BIP （Band-Interleaved Pixels） sequencing algorithm， a snapshot mosaic hyperspectral sample can be dynamically processed per clock cycle to 11 times higher of data throughput， and the problems of pipelining and parallelization difficulties and slow processing speed in the FPGA hardware implementation of snapshot mosaic hyperspectral image compression of this international standard have been effectively solved. The experimental results show that the overall logic resource usage of the implemented compressor is less than 10% after the layout and wiring of the FPGA on Xilinx XC7Z020CLG-4 platform， and the compression of a hyperspectral image can be completed in about 22 ms under 100 MHz system clock， with lossless compression performance between 2.66~4.30 bits/samples and the performance near-lossless compression is between 1.01~3.70 bits/samples， which can meet the application requirements of snapshot mosaic hyperspectral imaging technology in wireless handheld and unmanned airborne fields. 高光谱图像压缩 快照马赛克高光谱成像 CCSDS 123.0-B-2 FPGA实现 hyperspectral image compression snapshot mosaic hyperspectral imaging CCSDS 123.0-B-2 FPGA implementation 1 引　言 1 高光谱成像技术结合了数字成像和光谱学的优势 ［ 1 1 ］ ，可以用更高的精度和细节表征场景中的对象，能够应用于遥感 ［ 2 2 ］ 、工业制造 ［ 3 3 ］ 、医学 ［ 4 4 ］ 等领域。近年来，欧洲微电子研究中心（IMEC）推出的新一代快照马赛克高光谱成像技术，通过在CMOS感光面表面沉积不同波段的滤光膜，可同时获取25个不同光谱波段的高光谱数据 ［ 5 5 ， 6 6 ］ ，使得高光谱技术在手持设备或小型无人机遥感领域的应用成为了可能 ［ 7 7 ］ 。 由于包含了多通道的光谱数据信息，高光谱图像通常数据量巨大 ［ 8 8 ］ ，对数据传输带宽要求非常高。国际空间数据系统咨询委员会（The Consultative Committee for Space Data Systems，CCSDS）于2012年5月推出了多光谱/高光谱无损压缩标准CCSDS-123.0-B-1 ［ 9 9 ， 10 10 ］ ，并于2019年2月升级到CCSDS 123.0-B-2标准 ［ 11 11 ］ 。新标准增加了对近无损压缩的支持并向下兼容，拓展了适用范围 ［ 11 11 ］ ，在减少数据带宽、减少缓存和存储的需求以及减少数据传输时间等方面进行了优化 ［ 12 12 ］ ，适合在卫星、无人机或微型无线高光谱采集设备上应用 ［ 13 13 ］ 。 但是CCSDS 123.0-B-2标准新引入的近无损压缩功能会导致算法在FPGA硬件实现时无法有效使用流水线化和并行化技术 ［ 14 14 ， 15 15 ］ ，降低了数据吞吐率，限制了其在快照马赛克高光谱图像压缩中的应用。本文通过改进BIP（Band-Interleaved Pixels）排序算法，尝试在CCSDS 123.0-B-2标准近无损压缩的FPGA硬件实现中引入流水线化和并行化技术，以提高系统的压缩效率。 2 CCSDS 123.0-B-2标准分析 1 CCSDS 123.0-B-2标准适用于由多光谱和高光谱成像光谱仪采集的三维光谱数据集的无损/近无损压缩，支持从2位到32位整数样本值 <math id="M1"><msub><mrow><mi>s</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051390&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051389&type= 6.26533365 3.72533321 所构成的三维矩阵（ <math id="M2"><msub><mrow><mi>N</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>,</mo><msub><mrow><mi>N</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>,</mo><msub><mrow><mi>N</mi></mrow><mrow><mi>z</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051342&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051340&type= 15.83266735 4.14866638 ）数据集。该标准由预测器和编码器两部分组成，原理图如 图1 图1 所示。 预测器（Predictor）使用FLEX（Fast Lossless Extended）算法，通过待测样本的三维邻域获得对应的预测值，计算数据残差以达到减少数据量的目的。预测器典型邻域如 图2 图2 所示，图中 P 为取值范围为［0，15］的用户自定义整数。 预测器的预测过程如下：首先计算同一光谱通道 <math id="M3"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051375&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051374&type= 1.69333339 2.79399991 内与样本 <math id="M4"><msub><mrow><mi>s</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051390&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051389&type= 6.26533365 3.72533321 相邻的样本代表值 <math id="M5"><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051394&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051392&type= 7.53533268 4.14866638 的局部和 <math id="M6"><msub><mrow><mi>σ</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051397&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051384&type= 6.68866634 3.72533321 ，然后通过局部和获得当前光谱通道的方向局部差 <math id="M7"><msubsup><mrow><mi>d</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow><mrow><mi mathvariant="normal">N</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051407&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051400&type= 7.11199999 3.97933316 、 <math id="M8"><msubsup><mrow><mi>d</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow><mrow><mi mathvariant="normal">W</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051411&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051410&type= 7.11199999 3.97933316 和 <math id="M9"><msubsup><mrow><mi>d</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow><mrow><mi mathvariant="normal">N</mi><mi mathvariant="normal">W</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051416&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051425&type= 7.11199999 3.97933316 及先前光谱通道的 P 个中央局部差 <math id="M10"><msub><mrow><mi>d</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051435&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051417&type= 7.11199999 3.72533321 。上述四种局部差共同组成局部差向量 <math id="M11"><msub><mrow><mi mathvariant="bold-italic">U</mi></mrow><mrow><mi>z</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=27051452&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051450&type= 9.05933285 4.48733330 的各个分量。局部差向量 <math id="M12"><msub><mrow><mi mathvariant="bold-italic">U</mi></mrow><mrow><mi>z</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=27051452&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051450&type= 9.05933285 4.48733330 使用与其具有相同维度的权重向量 <math id="M13"><msub><mrow><mi mathvariant="bold-italic">W</mi></mrow><mrow><mi>z</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=27051456&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051454&type= 9.82133293 4.48733330 加权和来计算对应的预测样本值 <math id="M14"><msub><mrow><mover accent="true"><mi>s</mi><mo>̂</mo></mover></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051459&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051457&type= 6.26533365 3.80999994 ，权重向量会随着高光谱图像的特点自适应更新。预测样本值与实际样本值做差得到的预测残差 <math id="M15"><msub><mrow><mi mathvariant="normal">Δ</mi></mrow><mrow><mi>z</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=27051470&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051468&type= 8.63599968 4.48733330 使用步长为 <math id="M16"><mn mathvariant="normal">2</mn><msub><mrow><mi>m</mi></mrow><mrow><mi>z</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>+</mo><mn mathvariant="normal">1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051476&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051473&type= 17.35666656 4.48733330 的均匀量化器进行量化，计算有符号整数量化器索引 <math id="M17"><msub><mrow><mi>q</mi></mrow><mrow><mi>z</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=27051529&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051526&type= 7.95866632 4.82600021 ，使得近无损重建高光谱图像时，可以获得误差不超过 <math id="M18"><msub><mrow><mi>m</mi></mrow><mrow><mi>z</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=27051502&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051483&type= 8.97466660 4.48733330 的样本 <math id="M19"><msub><mrow><mi>s</mi></mrow><mrow><mi>z</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=27051509&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051506&type= 7.36600018 4.48733330 。当 <math id="M20"><msub><mrow><mi>m</mi></mrow><mrow><mi>z</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051524&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051511&type= 15.40933418 4.48733330 时，预测器为无损压缩。最后，有符号量化索引 <math id="M21"><msub><mrow><mi>q</mi></mrow><mrow><mi>z</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=27051529&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051526&type= 7.95866632 4.82600021 被映射为 <math id="M22"><mi>D</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051539&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051521&type= 2.79399991 2.79399991 位无符号映射量化索引 <math id="M23"><msub><mrow><mi>δ</mi></mrow><mrow><mi>z</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=27051562&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051560&type= 8.04333401 4.48733330 作为预测器的输出。 编码器（Encoder）将 <math id="M24"><msub><mrow><mi>δ</mi></mrow><mrow><mi>z</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=27051559&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051557&type= 8.04333401 4.48733330 编码为一种可变长二进制比特流文件，该文件在图像重建时可完整恢复 <math id="M25"><msub><mrow><mi>δ</mi></mrow><mrow><mi>z</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=27051559&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051557&type= 8.04333401 4.48733330 。文件由可变长文件头和文件主体组成，可变长文件头包含压缩算法和被压缩图像的相关参数，文件主体则是映射量化索引的编码。 CCSDS 123.0-B-2标准推荐了三种编码器实现方式 ［ 11 11 ］ ，其中样本自适应编码器的泛用性更强。该编码器是基于指数哥伦布（Golomb-Power-Of-2， GPO2）编码的改进编码算法。编码算法的原理是使用 k 阶长度限制GPO2编码对非负整数 δ z （ t ）进行编码，可以将 <math id="M26"><msub><mrow><mi>δ</mi></mrow><mrow><mi>z</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=27051562&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051560&type= 8.04333401 4.48733330 表示为： <math display="block" id="M27"><msub><mrow><mi>δ</mi></mrow><mrow><mi>z</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mi>u</mi><mo>·</mo><msup><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mi>k</mi></mrow></msup><mo>+</mo><mi>r</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051572&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051563&type= 24.63800049 4.65666676 . 对 δ z （ t ）的编码实际上就是对 u 和 r 的计算，其中 <math id="M28"><mi>u</mi><mo>=</mo><mo>⌊</mo><mi>δ</mi><mo>/</mo><msup><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mi>k</mi></mrow></msup><mi mathvariant="normal">」</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051583&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051574&type= 18.37266731 3.80999994 （向下取整）， <math id="M29"><mi>r</mi><mo>=</mo><mi>δ</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">d</mi><msup><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mi>k</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051581&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051585&type= 17.61066628 3.04800010 （取余）。编码器通过对各光谱通道都设置了单独的累加器 <math id="M30"><msub><mo>∑</mo><mrow><mi>z</mi></mrow></msub><mrow><mo stretchy="false">(</mo></mrow><mi>t</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051600&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051588&type= 9.22866726 4.48733330 和计数器 <math id="M31"><msub><mrow><mi>Γ</mi></mrow><mrow><mi>z</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=27051604&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051592&type= 8.72066593 4.48733330 获得 k 值进而实现对 <math id="M32"><msub><mrow><mi>δ</mi></mrow><mrow><mi>z</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=27051610&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051596&type= 8.04333401 4.48733330 的编码。 CCSDS 123.0-B-2标准最重要的优势在于支持近无损压缩，即一种能够控制重建图像和原始图像之间的最大误差的有损压缩。考虑到新加入的近无损压缩功能会使得图像重建时无法使用原始值恢复高光谱图像，该标准引入了样本代表值的计算，在每个样本量化后计算出样本代表值 <math id="M33"><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051618&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051615&type= 7.53533268 4.14866638 用于下一个样本的预测。但是引入样本代表值后，在FPGA实现的预处理阶段使用推荐的样本输入排序方式时，会因为反馈环路过长、数据依赖不满足等原因，导致对快照马赛克高光谱图像的压缩吞吐率降低。本文将在CCSDS 123.0-B-2标准的FPGA实现中对推荐的样本输入排序算法进行改进，以克服这一不足。 3 CCSDS 123.0-B-2标准的FPGA实现 1 3.1　排序预处理的改进 2 在使用FPGA对高光谱图像进行压缩实现前，首先需要确定各光谱通道数据的样品输入顺序。CCSDS 123.0-B-2标准推荐了BI（Band Interleaved）和BSQ（Band Sequential）两种样本输入排序方式，其中BI排序又可分为BIP（Band-Interleaved Pixels）排序和BIL（Band-Interleaved Lines）排序两类。 图3 图3 为三种推荐样本输入排序方式的原理示意图。 由于CCSDS 123.0-B-2标准包含了 “样本代表值反馈”闭环反馈环路，在使用BSQ排序和BIL排序进行硬件实现时，会有如 图4 图4 所示的时序。从图中可以看到，只有在完成整个反馈回路得到Sample1的样本代表值 <math id="M34"><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</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=27051647&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051627&type= 8.63599968 4.57200003 后，才能开始Sample2的压缩。这就导致吞吐率受限于处理“样本代表值反馈”环路所需时钟周期，使得平均需要十几个时钟才能完成对一个样本的压缩，大幅降低了压缩器的吞吐率。 BIP排序的优势在于硬件实现时不需要等待整个反馈回路完成，就能开始下一个样本的压缩，可以克服BIL和BSQ排序造成的吞吐率低的问题。但是，CCSDS 123.0-B-2标准对样本 <math id="M35"><msub><mrow><mi>s</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051673&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051658&type= 6.26533365 3.72533321 的预测需要使用如 公式（2） 公式（2） 所示的局部差向量 <math id="M36"><msub><mrow><mi mathvariant="bold-italic">U</mi></mrow><mrow><mi>z</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=27051684&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051683&type= 9.05933285 4.48733330 ，而 <math id="M37"><msub><mrow><mi mathvariant="bold-italic">U</mi></mrow><mrow><mi>z</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=27051684&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051683&type= 9.05933285 4.48733330 的计算用到了样本 <math id="M38"><msub><mrow><mi>s</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051668&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051687&type= 9.73666668 3.72533321 的代表值 <math id="M39"><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051726&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051716&type= 11.09133339 4.14866638 。由于BIP排序应用在快照马赛克高光谱图像的流水线压缩时， <math id="M40"><msub><mrow><mi>s</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051696&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051693&type= 9.73666668 3.72533321 必须在 <math id="M41"><msub><mrow><mi>s</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051702&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051709&type= 6.26533365 3.72533321 前一个时钟的开始时刻计算。如 图5 图5 所示，这也就意味着在对样本 <math id="M42"><msub><mrow><mi>s</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051722&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051720&type= 6.26533365 3.72533321 预测时，无法得到样本代表值 <math id="M43"><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051726&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051716&type= 11.09133339 4.14866638 ，根据数据的因果关系，将导致无法进行并行化压缩处理。 <math display="block" id="M44"><msub><mrow><mi mathvariant="bold-italic">U</mi></mrow><mrow><mi>z</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msubsup><mrow><mi>d</mi></mrow><mrow><mi>z</mi></mrow><mrow><mi mathvariant="normal">N</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><msubsup><mrow><mi>d</mi></mrow><mrow><mi>z</mi></mrow><mrow><mi mathvariant="normal">W</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><msubsup><mrow><mi>d</mi></mrow><mrow><mi>z</mi></mrow><mrow><mi mathvariant="normal">N</mi><mi mathvariant="normal">W</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>d</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>d</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>d</mi></mrow><mrow><mi>z</mi><mo>-</mo><mi>P</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mn mathvariant="normal">4</mn><msubsup><mrow><mi>s</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>x</mi></mrow><mrow><mo>″</mo></mrow></msubsup><mo>-</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></mtd></mtr><mtr><mtd><mn mathvariant="normal">4</mn><msubsup><mrow><mi>s</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow><mrow><mo>″</mo></mrow></msubsup><mo>-</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></mtd></mtr><mtr><mtd><mn mathvariant="normal">4</mn><msubsup><mrow><mi>s</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>x</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow><mrow><mo>″</mo></mrow></msubsup><mo>-</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></mtd></mtr><mtr><mtd><mn mathvariant="normal">4</mn><msubsup><mrow><mi>s</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow><mrow><mo>″</mo></mrow></msubsup><mo>-</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></mtd></mtr><mtr><mtd><mn mathvariant="normal">4</mn><msubsup><mrow><mi>s</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">2</mn><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow><mrow><mo>″</mo></mrow></msubsup><mo>-</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>z</mi><mo>-</mo><mn mathvariant="normal">2</mn><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><mn mathvariant="normal">4</mn><msubsup><mrow><mi>s</mi></mrow><mrow><mi>z</mi><mo>-</mo><mi>p</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow><mrow><mo>″</mo></mrow></msubsup><mo>-</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>z</mi><mo>-</mo><mi>p</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051748&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051729&type= 63.24599838 37.59199905 . 本文通过改进的BIP排序方式，将BIP排序在 <math id="M47"><mi>x</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051781&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051779&type= 2.03200006 2.79399991 和 <math id="M48"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051775&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051785&type= 1.69333339 2.79399991 方向上的扫描方式由顺序扫描改为交错扫描，可以有效解决上述问题。如 图6 图6 所示，交错扫描方式通过分离各个光谱通道的压缩顺序，将不同光谱通道相同空间位置样本的处理时间间隔由1个时钟周期改为 N z +1个时钟周期，克服了推荐BIP排序存在的数据依赖问题，同时保留了BIP易于流水线处理、吞吐率高的特点，使得基于新标准实现的FPGA压缩器平均一个时钟可以动态处理一个像素样本。 3.2　压缩器的FPGA实现 2 CCSDS 123.0-B-2标准的硬件实现架构如 图7 图7 所示，架构主要由预测器和编码器两部分组成，使用了多个缓冲器和存储器对数据进行暂存和同步。 在典型的快照马赛克高光谱图像采集系统中，样本会以帧的形式逐个输入到FPGA上部署的压缩器，每个时钟周期都会有一个高光谱图像样本被传输到输入缓存中。预测器按照改进的BIP排序读取样本并开始流水线处理，包含局部和与局部差计算、点乘、预测、量化、映射、权重更新和样本代表计算等一系列操作。预测器还会构建多个存储器，用于储存计算过程中所用到的权重向量、样本代表值和局部差向量等中间数据。编码器则负责样本自适应熵编码。系统的输出缓存将二进制编码重新组织为标准BIP排序，使得通用的解压缩器能够无损或近无损重建原有高光谱图像。 图8 图8 为本文FPGA硬件实现方案对应的流水线时序。 使用Verilog语言将上述压缩器的硬件设计在Vivado2019.2集成开发工具中综合实现，使用Xilinx xc7z020clg400-2 FPGA实现的压缩器顶层结构如 图9 图9 所示。本设计不依赖外部存储资源，共使用了6 731个六输入LUT，2 711个触发器，9.5个4 KB BRAM，以及11个DSP48，占用xz7020逻辑资源少于10%，总的逻辑资源使用情况如 表1 表1 所示。 Resource Utilization Available Utilization LUT 6 731 53 200 12.65% LUTRAM 518 17 400 2.98% FF 2 731 106 400 2.57% BRAM 9.50 140 6.79% DSP 11 220 5.0% BUFG 1 32 3.13% 3.3　快照马赛克高光谱图像输入缓存模块的设计 2 快照马赛克高光谱图像由欧洲微电子研究中心（IMEC）推出的新一代的快照马赛克高光谱传感器采集，它通过在CMOS感光面表面沉积不同的滤光片实现快照马赛克高光谱图像的获取。如 图10 图10 所示，由相邻的25个像素构成一个“成像单元”，可同时获取到25个不同光谱波段的8位光谱数据，传感器像素尺寸为1 088×2 048。 受限于上述独特的“成像单元”，针对普通光谱数据立方体提出的CCSDS 123.0-B-2标准应用于快照马赛克高光谱图像的压缩时有一定限制。因此本文不再使用快照马赛克高光谱图像内相邻像素作为上下文邻域计算局部和与局部差，而是利用“成像单元”每5个像素规律排布的特点的调整邻域，改进邻域后的局部和和局部差计算方式如 式（3） 式（3） 所示： <math display="block" id="M49"><msub><mrow><mi>σ</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub><mo>=</mo><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi><mo>-</mo><mn mathvariant="normal">5</mn></mrow></msub><mo>+</mo><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn mathvariant="normal">5</mn><mo>,</mo><mi>x</mi><mo>-</mo><mn mathvariant="normal">5</mn></mrow></msub><mo>+</mo><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn mathvariant="normal">5</mn><mo>,</mo><mi>x</mi></mrow></msub><mo>+</mo><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn mathvariant="normal">5</mn><mo>,</mo><mi>x</mi><mo>+</mo><mn mathvariant="normal">5</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051868&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051863&type= 72.81333160 4.23333359 <math display="block" id="M50"><mtable><mtr><mtd><msub><mrow><mi>d</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub><mo>=</mo><mn mathvariant="normal">4</mn><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub><mo>+</mo><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn mathvariant="normal">5</mn><mo>,</mo><mi>x</mi></mrow></msub><mo>+</mo><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn mathvariant="normal">5</mn><mo>,</mo><mi>x</mi><mo>+</mo><mn mathvariant="normal">5</mn></mrow></msub></mtd></mtr><mtr><mtd><msubsup><mrow><mi>d</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow><mrow><mi mathvariant="normal">N</mi></mrow></msubsup><mo>=</mo><mn mathvariant="normal">4</mn><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn mathvariant="normal">5</mn><mo>,</mo><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></mtd></mtr><mtr><mtd><msubsup><mrow><mi>d</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow><mrow><mi mathvariant="normal">W</mi></mrow></msubsup><mo>=</mo><mn mathvariant="normal">4</mn><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi><mo>-</mo><mn mathvariant="normal">5</mn></mrow></msub><mo>-</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></mtd></mtr><mtr><mtd><msubsup><mrow><mi>d</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow><mrow><mi mathvariant="normal">N</mi><mi mathvariant="normal">W</mi></mrow></msubsup><mo>=</mo><mn mathvariant="normal">4</mn><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn mathvariant="normal">5</mn><mo>,</mo><mi>x</mi><mo>-</mo><mn mathvariant="normal">5</mn></mrow></msub><mo>-</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi></mrow></msub></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051887&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051870&type= 56.64199829 19.47333336 . 为了计算改进后的局部和与局部差，CCSDS 123.0-B-2标准在FPGA实现过程时使用输入缓存模块获取改进后的上下文邻域。如 图11 图11 所示，该模块使用BRAM开辟了一个深度为2× N Z × N X 的双端口RAM对样本代表值进行暂存和读出，当压缩开始时，写数据写地址先使能，当写入数据量达到 N Z × N X 时，读数据和读地址开始使能。从RAM读出的样本代表值经过FIFO0缓存后得到 <math id="M51"><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi><mo>-</mo><mn mathvariant="normal">5</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051899&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051891&type= 11.09133339 4.14866638 ， <math id="M52"><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>x</mi><mo>-</mo><mn mathvariant="normal">5</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051899&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051891&type= 11.09133339 4.14866638 被传入深度为达 N Z ×（ N X - 2）的FIFO1得到 <math id="M53"><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn mathvariant="normal">5</mn><mo>,</mo><mi>x</mi><mo>+</mo><mn mathvariant="normal">5</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051905&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051902&type= 14.64733315 4.14866638 ；同理使用 N Z 的FIFO2、FIFO3依次获得 <math id="M54"><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn mathvariant="normal">5</mn><mo>,</mo><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051909&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051914&type= 11.09133339 4.14866638 和 <math id="M55"><mi>s</mi><msub><mrow><mo>″</mo></mrow><mrow><mi>z</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn mathvariant="normal">5</mn><mo>,</mo><mi>x</mi><mo>-</mo><mn mathvariant="normal">5</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051928&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27051911&type= 14.64733315 4.14866638 。 输入缓存模块通过四级流水线提高处理带宽，第一级流水线用来生成改进的BIP排序的（ z ， y ， x ）序列，它会随着每个时钟更新；第二级流水线通过计数器生成双端口RAM的写地址，将（ z ， y ， x ）坐标映射为对应位置样本在双端口RAM内的地址；第三级流水线用来对RAM进行写入和读取数据；第四级流水线结合马赛克面阵高光谱图像的特点，获得（ z ， y ， x ）处样本的上下文邻域，输出给局部和与局部差计算模块。 4 压缩器FPGA实现的测试结果 1 目前快照马赛克高光谱图像没有类似AVIRIS、AIRS等的国际通用高光谱遥感图像测试数据集，因此本文使用课题组自行研制的快照马赛克微型高光谱仪进行了测试数据的采集。 图12 图12 为四种不同场景下所采集的快照马赛克高光谱图像。 4.1　时序仿真结果 2 使用Vivado集成的仿真工具Simulator对FPGA设计进行功能仿真和时序仿真， 图13 图13 为激励模块设计。 其中sys_clk和sys_rst_n为系统时钟信号和系统复位信号。复位信号拉高后，Testbench在每个时钟的上升沿都会向压缩器发送一个高光谱样本值。当压缩器获得第一个样本的压缩编码后，二进制编码写入信号code_write会拉高，通过code_string_o信号和code_length_o信号将编码值写入压缩码流文件，所有编码完成后，code_write将拉低以停止码流写入。 以River1的测试结果为例， 图14 图14 是使用CCSDS 123.0-B-2标准推荐的BSQ排序压缩器的时序， 图15 图15 是使用改进BIP排序压缩器的时序。 从时序图可以看到，基于改进BIP排序的压缩器平均一个时钟周期就可以完成一个样本的压缩编码，这是改进排序算法易于流水线和并行化处理的特点所带来的好处。 表2 表2 列出了两种排序方式下压缩器在100 MHz时钟下的性能对比。结果表明，虽然这种改进的排序方式不会改善压缩性能，但是会大大提高压缩吞吐率，压缩一张快照马赛克高光谱图像所需时间仅为22.2 ms，远小于原BSQ排序压缩器所需的242.9 ms，数据吞吐率约提高至11倍。 排序方式 编码时间/ms 压缩编码长度/bits 改进BIP排序 22.2 9 485 844 BSQ排序 242.9 9 485 844 4.2　无损压缩性能 2 表3 表3 列出了本文所实现的压缩器对四种场景下的快照马赛克高光谱图像的无损压缩性能。压缩性能通过bits/sample定义，即压缩后二进制码流长度除以原始图像中的样本数 N X · N Y · N Z ，其数值越小表明压缩效果越好。实验结果显示，最终的压缩性能在2.66~4.30 bits/sample之间，优于原CCSDS 123.0-B-1标准下2.68~4.32 bits/sample的压缩性能。 标准 River1 Grass1 Multiple material1 Red jujube1 CCSDS 123.0-B-1 4.319 6 3.767 4 3.322 2 2.684 4 CCSDS 123.0-B-2 4.295 0 3.744 1 3.283 2 2.661 5 bits/sample 4.3　近无损压缩性能结果分析 2 近无损压缩是一种能够控制重建图像和原始图像之间的最大误差的有损压缩。接下来将介绍压缩器近无损压缩性能，峰值信噪比（PSNR）和光谱角（SA）测试结果。 4 . 3 . 1　近无损压缩性能 3 表4 表4 为当绝对误差极限 a z =16，各场景下快照马赛克高光谱图像的近无损压缩性能。 reError River1 Grass1 Multiple material1 Red jujube1 2 3.696 9 3.608 2.870 8 2.639 5 4 2.724 3 2.878 6 2.117 1 2.047 8 1.828 1.820 7 1.387 4 1.438 2 16 1.361 6 1.340 9 1.089 7 1.084 4 32 1.176 3 1.128 9 1.020 3 1.019 1 64 1.144 8 1.013 5 1.013 5 1.009 6 127 1.143 8 1.012 9 1.012 9 1.008 1 bits/sample 从表中可以看出压缩器的近无损压缩性能约在1~3.70 bits/sample之间，远高于无损压缩性能。但是，近无损压缩性能的提高是以牺牲图像的保真度而获得的。 图16 图16 显示了不同压缩设置下图像细节的损失。 4 . 3 . 2　PSNR测试 3 PSNR是衡量有损压缩图像失真程度的最常用指标，通常来讲PSNR值越小压缩后的图像质量越差。它的计算方式如 公式（4） 公式（4） 所示，其中 MSE 为均方误差， L 为量化等级。 <math display="block" id="M56"><mi>P</mi><mi>S</mi><mi>N</mi><mi>R</mi><mo>=</mo><mn mathvariant="normal">10</mn><mi mathvariant="normal">l</mi><mi mathvariant="normal">g</mi><mfrac><mrow><msup><mrow><mi>L</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mi>M</mi><mi>S</mi><mi>E</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27052011&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27052024&type= 31.58066559 8.21266651 . 表5 表5 给出了3幅快照马赛克高光谱图像在不同的相对误差极限 r z 下压缩处理后的PSNR值。随着相对误差极限 r z 的增加，PSNR值逐渐减小，表明CCSDS 123.0-B-2标准近无损压缩性能的提高是以图像失真为代价。 reError River1 Grass1 Multiple material1 2 21.535 5 26.255 3 28.864 7 6 21.542 5 26.241 7 28.810 0 8 21.544 0 26.227 0 28.802 0 12 21.545 9 26.194 2 28.790 6 16 21.544 8 26.154 3 28.753 4 24 21.500 0 26.056 5 28.785 8 32 21.470 2 25.958 3 28.640 2 48 21.426 2 25.901 0 28.577 8 64 21.407 1 25.901 5 28.515 2 80 21.378 5 25.829 3 28.501 2 96 21.376 9 25.764 1 28.363 2 112 21.374 0 25.745 1 28.399 0 128 21.363 4 25.731 0 28.420 6 4 . 3 . 3　SA测试 3 SA用来表征两个高光谱数据之间光谱的相似性，两光谱向量间的夹角越小，两光谱越相似。定义如 公式（5） 公式（5） 所示。 <math display="block" id="M57"><mi>S</mi><mi>A</mi><mo>=</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><msup><mrow><mi mathvariant="normal">s</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfrac><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi>y</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mrow><mrow><mroot><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msubsup><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mrow><mrow/></mroot><mroot><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msubsup><mrow><mi>y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mrow><mrow/></mroot></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27052034&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27052017&type= 41.82533264 19.89666748 . 图17 图17 给出了快照马赛克高光谱图像SA的变化。由 图17 图17 可知，随着相对误差极限 r z 的增大，SA增大，压缩图像与原始图像的光谱相似程度越差，表明近无损压缩会破坏快照马赛克高光谱图像包含的光谱信息，并且随着压缩性能的提升，这种破坏效果愈加显著。 5 结　论 1 基于多光谱/高光谱无损和近无损压缩标准CCSDS 123.0-B-2，对马赛克高光谱数据的压缩进行了FPGA硬件实现。通过对该标准推荐的BIP排序方式的改进，大幅提高了压缩器的吞吐率，100 MHz的时钟下约22 ms即可完成一帧快照马赛克高光谱图像的压缩。所实现的压缩器的无损压缩性能在2.66~4.30 bits/samples之间，近无损压缩性能可达1.01~3.70 bits/samples，能够满足快照马赛克高光谱成像技术在无线手持及无人机载领域的应用需求。 参考文献 高恒振 . 高光谱遥感图像分类技术研究 ［D］. 长沙 ： 国防科学技术大学 ， 2011 . 高光谱遥感图像分类技术研究 GAO H Z . Research on Classification Techinique for Hyperspectral Remote Sensing Imagery ［D］. Changsha ： National University of Defense Technology ， 2011 . （in Chinese） Research on Classification Techinique for Hyperspectral Remote Sensing Imagery ESA Earth Observation Portal . PRISMA （Hyperspectral Precursor and Application Mission） ［EB/OL］. ［ 2021-01-16 ］ https：//directory.eoportal.org/web/eoportal/satellite-missions/p/prisma-hyperspectral https://directory.eoportal.org/web/eoportal/satellite-missions/p/prisma-hyperspectral . PRISMA （Hyperspectral Precursor and Application Mission） SUN D W . Hyperspectral Imaging for Food Quality Analysis and Control ［M］. London UK ： Academic Press ， 2010 Hyperspectral Imaging for Food Quality Analysis and Control LU G L ， FEI B W . Medical hyperspectral imaging： a review ［J］. Journal of Biomedical Optics ， 2014 ， 19 ： 010901 . doi: 10.1117/1.jbo.19.1.010901 http://dx.doi.org/10.1117/1.jbo.19.1.010901 Medical hyperspectral imaging： a review SHOEIBY M ， PETERSSON L ， ARMIN M A ， et al . Super-resolved chromatic mapping of snapshot mosaic image sensors via a texture sensitive residual network ［C］. 2020 IEEE Winter Conference on Applications of Computer Vision . 15，2020 ， Snowmass ， CO ， USA . IEEE ， 2020 ： 2793 - 2802 . doi: 10.1109/wacv45572.2020.9093518 http://dx.doi.org/10.1109/wacv45572.2020.9093518 Super-resolved chromatic mapping of snapshot mosaic image sensors via a texture sensitive residual network YUAN C ， ZHOU M ， SUN L ， et al . A multi-dimensional hyperspectral image mosaic method and its acquisition system ［C］. 2018 11th International Congress on Image and Signal Processing， BioMedical Engineering and Informatics （CISP-BMEI）. 1315，2018 ， Beijing， China. IEEE ， 2018 ： 1 - 6 . doi: 10.1109/cisp-bmei.2018.8633038 http://dx.doi.org/10.1109/cisp-bmei.2018.8633038 A multi-dimensional hyperspectral image mosaic method and its acquisition system HAGEN N A ， KUDENOV M W . Review of snapshot spectral imaging technologies ［J］. Optical Engineering ， 2013 ， 52 ： 090901 . doi: 10.1117/1.oe.52.9.090901 http://dx.doi.org/10.1117/1.oe.52.9.090901 Review of snapshot spectral imaging technologies DUA Y M ， KUMAR V ， SINGH R S . Comprehensive review of hyperspectral image compression algorithms ［J］. Optical Engineering ， 2020 ， 59 ： 090902 . doi: 10.1117/1.oe.59.9.090902 http://dx.doi.org/10.1117/1.oe.59.9.090902 Comprehensive review of hyperspectral image compression algorithms Consultative Committee for Space Data Systems . Lossless Multispectral and Hyperspectral Image Compression-CCSDS 120.2-G-1 ［S］. Washington DC， USA ： In Green Book ， 2015 . doi: 10.1109/lascas.2015.7250438 http://dx.doi.org/10.1109/lascas.2015.7250438 Consultative Committee for Space Data Systems . LOS SLESS MULTISPECTRAL & HYPERSPECTRAL IMAGE COMPRESSION CCSDS 123.0-B-1 ［S］. Washington DC， USA ： In Blue Book ， 2012 . Consultative Committee for Space Data Systems . Low-Complexity Lossless and Near-Lossless Multispectral and Hyperspectral Image Compression-CCSDS 123.0-B-2 ［S］. Washington DC， USA ： In Blue Book ， 2019 . doi: 10.3390/rs11111390 http://dx.doi.org/10.3390/rs11111390 BLANES I ， KIELY A ， HERNÁNDEZ-CABRONERO M ， et al . Performance impact of parameter tuning on the CCSDS-123.0-B-2 low-complexity lossless and near-lossless multispectral and hyperspectral image compression standard ［J］. Remote Sensing ， 2019 ， 11 （ 11 ）： 1390 . doi: 10.3390/rs11111390 http://dx.doi.org/10.3390/rs11111390 Performance impact of parameter tuning on the CCSDS-123.0-B-2 low-complexity lossless and near-lossless multispectral and hyperspectral image compression standard 张宁 ， 冯书谊 ， 濮建福 ， 等 . 多光谱遥感图像CCSDS动态码率控制近无损压缩 ［J］. 光学 精密工程 ， 2015 ， 23 （ 6 ）： 1783 - 1790 . doi: 10.3788/ope.20152306.1783 http://dx.doi.org/10.3788/ope.20152306.1783 多光谱遥感图像CCSDS动态码率控制近无损压缩 ZHANG N ， FENG S 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） . doi: 10.3788/ope.20152306.1783 http://dx.doi.org/10.3788/ope.20152306.1783 Dynamic rate control for CCSDS nearly lossless compression of multispectral remote image HERNÁNDEZ-CABRONERO M ， KIELY A B ， KLIMESH M ， et al . The CCSDS 123.0-B-2 “Low-Complexity Lossless and Near-Lossless Multispectral and Hyperspectral Image Compression” Standard： a comprehensive review ［J］. IEEE Geoscience and Remote Sensing Magazine ， 2021 ， 9 （ 4 ）： 102 - 119 . doi: 10.1109/mgrs.2020.3048443 http://dx.doi.org/10.1109/mgrs.2020.3048443 The CCSDS 123.0-B-2 “Low-Complexity Lossless and Near-Lossless Multispectral and Hyperspectral Image Compression” Standard： a comprehensive review KIELY A B ， KLIMESH M ， BLANES I ， et al . The new CCSDS Standard for Low-Complexity Lossless and Near-Lossless Multispectral and Hyperspectral Image Compression ［C］. International Workshop on On-Board Payload Data Compression . 2018 （GSFC-E-DAA-TN61202）. The new CCSDS Standard for Low-Complexity Lossless and Near-Lossless Multispectral and Hyperspectral Image Compression 10.37188/OPE.20223008.0883 TP391.4 A 1004-924X（2022）08-0883-11 吕德远,黄浩文,曾珉等.快照马赛克高光谱图像CCSDS压缩器的FPGA实现[J].光学精密工程,2022,30(08):883-893. 吕德远,黄浩文,曾珉等.快照马赛克高光谱图像CCSDS压缩器的FPGA实现[J].光学精密工程,2022,30(08):883-893. LV Deyuan,HUANG Haowen,ZENG Min,et al.FPGA implementation of CCSDS compressor for snapshot mosaic hyperspectral images[J].Optics and Precision Engineering,2022,30(08):883-893. LV Deyuan,HUANG Haowen,ZENG Min,et al.FPGA implementation of CCSDS compressor for snapshot mosaic hyperspectral images[J].Optics and Precision Engineering,2022,30(08):883-893. 横向课题项目（No.20200202002）；校企合作项目（No.2018010802012）；广东省科技创新战略专项资金项目（No.pdjh2020b0526）；深圳技术大学校级科研项目（No.2020028555301020，No.JSZZ2020006）；研究生校企合作项目（No.XQHZ202017，No.XQHZ202008，No.XQHZ202102，No.XQHZ202120） 2021-04-26 2021-06-06 2022-04-25 2022-04-26 光学精密工程 8 30