CN105100810A - 一种成像声纳实时处理系统中的图像压缩解压方法及系统 - Google Patents

Abstract

本发明提供一种成像声纳实时处理系统中的图像压缩解压方法及系统，所述压缩方法包含：步骤101)对实时的声纳图像的行数据进行DCT变换，得到DCT系数；步骤102)将得到的DCT系数划分为S段，截取第一段DCT系数；步骤103)将截取的第一段DCT系数再划分为若干子段，然后对各子段包含的DCT系数分别进行量化压缩处理，完成图像压缩。其中，所述的量化压缩处理的原则为：针对较低频段数据采用较多位数进行整型量化，针对较高频段数据采用较少位数进行整型量化。本发明利用DCT变换的能量集中特性进行DCT系数截取，有很高的压缩率；利用分段量化进一步提高压缩率；对声纳行数据进行连续处理，实现声纳图像压缩传输的实时性。

Description

一种成像声纳实时处理系统中的图像压缩解压方法及系统

技术领域

本发明涉及声纳信号处理的技术领域，特别涉及一种成像声纳实时处理系统中的图像压缩解压方法及系统。

背景技术

成像声纳是利用水下声波对目标成像，进而进行探测和定位的设备。随着声纳技术的发展，以及计算机处理能力的不断进步，军民两用领域均对成像声纳的处理速度和处理效率提出了更高的要求。处理速度需求的提高驱使产生了声纳实时处理系统，处理效率需求的提高使得声纳的技术指标，尤其是测绘带宽度有了很大的提升。

上述背景决定了成像声纳的两个性质：实时性和大数据量。实时性带来的特点是，在声纳系统中，数据以流的方式进行传输和交互，理论上可以无限期运行，从而忽略了文件的概念。这就表明，原始数据在系统中是以每次一帧或数帧的方式周期传输的，图像是以每次一行或者若干行的方式周期传输的。因此原始数据和图像数据在传输时均为一维的概念，而不是从文件角度考虑的二维概念。

随着声纳平台的多样化发展，部分脱离母船的成像声纳平台要求图像数据在空中做无线传输，以供母船实时监视和操控。典型的有半潜式航行器，也称为遥控多功能航行器(RemoteMulti-MissionVehicle,RMMV)，使用这种平台的如洛克希德·马丁公司推出的搭载有AN/AQS-20A声纳的遥控猎雷系统(RemoteMinehuntingSystem，RMS)。当前无线传输的最大问题是带宽不足，当声纳在湖泊或海洋的恶劣环境中使用时，带宽更加受到限制。而应用需求又带来数据量大的特点，为此亟需改进现有的声纳图像压缩传输的技术。

发明内容

本发明的目的在于，为克服采用现有技术的多种平台的声纳实时系统中，声纳处理子系统所在平台与母船显控子系统之间通信时可能带宽不足的问题，本发明提供一种成像声纳实时处理系统中的图像压缩解压方法及系统。

为实现上述目的，本发明提供了一种成像声纳实时处理系统中的图像压缩方法，所述压缩方法包含：

步骤101)对实时的声纳图像的行数据进行DCT变换，得到DCT系数；

步骤102)将得到的DCT系数划分为S段，截取第一段DCT系数；

步骤103)将截取的第一段DCT系数再划分为若干子段，然后对各子段包含的DCT系数分别进行量化压缩处理，完成图像压缩；

其中，所述的量化压缩处理的原则为：针对较低频段数据采用较多位数进行整型量化，针对较高频段数据采用较少位数进行整型量化。

可选的，DCT变换的公式为：

$y\left(k\right)=\mathrm{\ω}\left(k\right)\underset{n=1}{\overset{N}{\mathrm{\Σ}}}x\left(n\right)\cos \left(\frac{\mathrm{\π}(2n-1)(k-1)}{2N}\right),k=\mathrm{1,2},...N$

其中，x(n)表示输入信号序列；N为序列的点数；y(k)为DCT变换后得到的系数。

可选的，上述步骤103)进一步包含：

步骤103-1)将截取的第一段DCT系数划分为4个子段，其中第一子段占截取段落的1/8，第二子段占截取段落的1/8，第三子段占截取段落的1/4，第四子段占截取段落的1/2；

步骤103-2)对第一子段的浮点型数据使用16位整型量化，具体的量化过程为：

设第一子段的DCT系数对应的数据为y₁(k)，首先求得y₁(k)序列的绝对值的最大值，再将整个y₁(k)序列归一化处理，其后乘以16位有符号整型量化范围的一半，再就近取整，计算公式为：

${y_1}^{\′}\left(k\right)=\mathrm{round}(\frac{y_1\left(k\right)}{\max \left(\mathrm{abs}\left(y_1\left(k\right)\right)\right)}\×2^{15})$

其中，函数abs()表示求绝对值，函数round()表示就近取整，函数max()表示求序列的最大值；

步骤103-3)采用如下三个公式分别对第二、三、四子段的数据进行量化处理，其中第二子段使用12位整型量化，第三子段使用8位整型数据量化，第四子段使用4位整型数据量化：

${y_2}^{\′}\left(k\right)=\mathrm{round}(\frac{y_2\left(k\right)}{\max \left(\mathrm{abs}\left(y_2\left(k\right)\right)\right)}\×2^{11})$

${y_3}^{\′}\left(k\right)=\mathrm{round}(\frac{y_3\left(k\right)}{\max \left(\mathrm{abs}\left(y_3\left(k\right)\right)\right)}\×2^7)$

${y_4}^{\′}\left(k\right)=\mathrm{round}(\frac{y_4\left(k\right)}{\max \left(\mathrm{abs}\left(y_4\left(k\right)\right)\right)}\×2^3)$

其中，y_i(k)表示划分得到的第i个子段，其中i＝1,2,3,4，y_i'(k)表示对第i子段量化后得到的序列。

可选的，上述的段数S应满足如下公式：

$\mathrm{\ρ}\left(S\right)=\frac{\underset{k=1}{\overset{k=\frac{N}{S}}{\mathrm{\Σ}}}{y}^2\left(k\right)}{\underset{k=1}{\overset{k=N}{\mathrm{\Σ}}}{y}^2\left(k\right)}$

ρ(S)＞ρ₀

其中，ρ(S)为前1/S段DCT系数的均方和与整个DCT系数序列的均方和的比值，ρ₀为设置的能量含量阈值。

针对上述压缩方法，本发明还提供了一种成像声纳实时处理系统中的图像解压方法，所述解压方法包含：

步骤201)根据反量化表达式对接收的数据进行反量化，对第一子段，第二子段、第三子段和第四个子段的数据的反量化公式分别如下：

$\left\{\begin{array}{c}{\tilde{y}}_1\left(k\right)=\frac{{y_1}^{\′}\left(k\right)}{2^{15}}\×\max \left(\mathrm{abs}\left(y_1\left(k\right)\right)\right)\\ {\tilde{y}}_2\left(k\right)=\frac{{y_2}^{\′}\left(k\right)}{2^{11}}\×\max \left(\mathrm{abs}\left(y_2\left(k\right)\right)\right)\\ {\tilde{y}}_3\left(k\right)=\frac{{y_3}^{\′}\left(k\right)}{2^7}\×\max \left(\mathrm{abs}\left(y_3\left(k\right)\right)\right)\\ {\tilde{y}}_4\left(k\right)=\frac{{y_4}^{\′}\left(k\right)}{2^3}\×\max \left(\mathrm{abs}\left(y_4\left(k\right)\right)\right)\end{array}\right.$

将反量化后得到的四段数据连接为一段数据，公式为：

$\tilde{y}\left(k\right)=\left[\begin{array}{cccc}{\tilde{y}}_1\left(k\right)& {\tilde{y}}_2\left(k\right)& {\tilde{y}}_3\left(4\right)& {\tilde{y}}_4\left(k\right)\end{array}\right]$

步骤201)将反量化后的数据进行DCT反变换，得到实际的图像数据；所述的DCT反变换为：

$x\left(n\right)=\underset{k=1}{\overset{N}{\mathrm{\Σ}}}\mathrm{\ω}\left(k\right)y\left(k\right)\cos \left(\frac{\mathrm{\π}(2n-1)(k-1)}{2N}\right),n=\mathrm{1,2},...N$

其中

$\mathrm{\ω}\left(k\right)=\left\{\begin{array}{c}\frac{1}{\sqrt{N}},k=1\\ \sqrt{\frac{2}{N}},2\≤k\≤N\end{array}\right.$

其中，N为对实时的声纳图像的行数据进行DCT变换得到的DCT系数的总长度。

此外，本发明还提供了一种用于成像声纳的图像压缩及解压系统，所述系统压缩子系统和解压子系统，所述压缩子系统包含：

DCT变换模块，用于对实时的声纳图像的行数据进行DCT变换，得到DCT系数；

截取模块，用于将得到的DCT系数划分为S段，截取第一段DCT系数；

分段量化处理模块，用于将截取的第一段DCT系数再划分为若干子段，然后对各子段包含的DCT系数分别进行量化压缩处理，完成图像压缩；

所述解缩子系统包含：

反量化处理模块，用于根据反量化表达式对接收的数据进行反量化，对第一子段，第二子段、第三子段和第四个子段的数据的反量化公式分别如下：

$\left\{\begin{array}{c}{\tilde{y}}_1\left(k\right)=\frac{{y_1}^{\′}\left(k\right)}{2^{15}}\×\max \left(\mathrm{abs}\left(y_1\left(k\right)\right)\right)\\ {\tilde{y}}_2\left(k\right)=\frac{{y_2}^{\′}\left(k\right)}{2^{11}}\×\max \left(\mathrm{abs}\left(y_2\left(k\right)\right)\right)\\ {\tilde{y}}_3\left(k\right)=\frac{{y_3}^{\′}\left(k\right)}{2^7}\×\max \left(\mathrm{abs}\left(y_3\left(k\right)\right)\right)\\ {\tilde{y}}_4\left(k\right)=\frac{{y_4}^{\′}\left(k\right)}{2^3}\×\max \left(\mathrm{abs}\left(y_4\left(k\right)\right)\right)\end{array}\right.$

将反量化后得到的四段数据连接为一段数据，公式为：

$\tilde{y}\left(k\right)=\left[\begin{array}{cccc}{\tilde{y}}_1\left(k\right)& {\tilde{y}}_2\left(k\right)& {\tilde{y}}_3\left(4\right)& {\tilde{y}}_4\left(k\right)\end{array}\right]$

DCT反变换，用于将反量化后的数据进行DCT反变换，得到实际的图像数据；所述的DCT反变换为：

$x\left(n\right)=\underset{k=1}{\overset{N}{\mathrm{\Σ}}}\mathrm{\ω}\left(k\right)y\left(k\right)\cos \left(\frac{\mathrm{\π}(2n-1)(k-1)}{2N}\right),n=\mathrm{1,2},...N$

其中

$\mathrm{\ω}\left(k\right)=\left\{\begin{array}{c}\frac{1}{\sqrt{N}},k=1\\ \sqrt{\frac{2}{N}},2\≤k\≤N\end{array}\right.$

其中，N为对实时的声纳图像的行数据进行DCT变换得到的DCT系数的总长度。

可选的，上述分段量化处理模块进一步包含：

分段子模块，根据DCT变换的能量集中特性，进一步将截取的第一段落划分为四个子段，第一子段占第一段落总长度的1/8，第二子段占第一段落总长度1/8，第三子段占第一段落总长度的1/4，第四子段占第一段落总长度的1/2；

量化子模块，使用不同精度的整型数据分别对各子段进行量化处理，其中第一子段使用16位整型量化，第二子段使用12位整型量化，第三子段使用8位整型量化，第四子段使用4位整型量化。

可选的，上述的段数S应满足下式：

$\mathrm{\ρ}\left(S\right)=\frac{\underset{k=1}{\overset{k=\frac{N}{S}}{\mathrm{\Σ}}}{y}^2\left(k\right)}{\underset{k=1}{\overset{k=N}{\mathrm{\Σ}}}{y}^2\left(k\right)}$

ρ(S)＞ρ₀

其中，ρ(S)为前1/S段DCT系数的均方和与整个DCT系数序列的均方和的比值，ρ₀为设置的能量含量阈值。

综上所述，本发明提供一种在实时处理系统中对图像数据进行压缩的技术。在实时处理端利用本方法进行压缩，在显控端利用逆向的方法进行解压缩。所述的压缩方法主要包括三步骤：(1)对实时的声纳图像的行数据进行DCT变换；(2)变换后利用DCT变换的能量集中特性，对DCT系数进行截断处理，截取其前1/S段系数，其中常数S可以根据实际情况灵活设定；(3)对截断后的DCT系数，再次利用能量集中特性，分段对DCT系数进行不同程度的量化压缩处理。

在显控端，主要的解压缩也分为三个步骤：(1)对接收到的数据，根据反量化表对数据进行反量化；(2)进行DCT反变换，得到实际的图像数据；(3)实时地输出图像至显示界面。

与现有技术相比，本发明的技术优势在于：

(1)利用DCT变换的能量集中特性进行DCT系数截取，有很高的压缩率；

(2)利用FFTW实现DCT变换，运算速度快；

(3)利用分段量化进一步提高压缩率；

(4)对声纳行数据进行连续处理，实现声纳图像压缩传输的实时性。

附图说明

图1是本发明提供的成像声纳实时处理系统中的图像压缩方法的处理流程图；

图2是本发明提供的分段量化的示意图；

图3是本发明的人工目标压缩处理结果(范围：200m×35m)；其中，图3(a)是原图，图3(b)是IR＝2的目标恢复图，图3(c)是IR＝4的目标恢复图，图3(d)是IR＝6的目标恢复图，图3(e)是IR＝8的目标恢复图，图3(f)是IR＝16的目标恢复图；

图4是人工目标压缩处理结果局部图(范围：10m×7.5m)，其中，图4(a)是原图，图4(b)是IR＝2的目标恢复图，图4(c)是IR＝4的目标恢复图，图4(d)是IR＝6的目标恢复图，图4(e)是IR＝8的目标恢复图，图4(f)是IR＝16的目标恢复图。

具体实施方式

下面结合附图对本发明的技术方案进行详细阐述。

本发明的关键在于两个部分：

(1)压缩阶段的DCT变换和解压阶段的DCT反变换；

(2)压缩阶段的量化处理和解压阶段的反量化处理。

以下分别进行陈述。

1DCT变换和DCT反变换

DCT变换常用于信号处理和图像处理，尤其多用于对信号或图像进行有损数据压缩。其特点是大多数自然信号经过变换后其能量主要集中在低频部分，而高频部分能量很少，这就是DCT变换的“能量集中”特性。利用该特性，仅使用少量的DCT系数就可以重建信号，并且信号失真不大。

DCT变换如下式所示。

$y\left(k\right)=\mathrm{\ω}\left(k\right)\underset{n=1}{\overset{N}{\mathrm{\Σ}}}x\left(n\right)\cos \left(\frac{\mathrm{\π}(2n-1)(k-1)}{2N}\right),k=\mathrm{1,2},...N$

其中

$\mathrm{\ω}\left(k\right)=\left\{\begin{array}{c}\frac{1}{\sqrt{N}},k=1\\ \sqrt{\frac{2}{N}},2\≤k\≤N\end{array}\right.$

相应的DCT反变换(也称IDCT变换)如下式所示。

$x\left(n\right)=\underset{k=1}{\overset{N}{\mathrm{\Σ}}}\mathrm{\ω}\left(k\right)y\left(k\right)\cos \left(\frac{\mathrm{\π}(2n-1)(k-1)}{2N}\right),n=\mathrm{1,2},...N$

其中，x(n)表示输入信号序列，N为序列的点数；y(k)为DCT变换后得到的系数，系数ω(k)与DCT变换的系数相同。

在具体实现中，可以使用FFTW函数库来实现DCT和IDCT的快速执行。FFTW函数库支持数据结构复杂的DCT和IDCT变换，是目前已知免费的计算FFT、DCT最快的函数库。

基于上述公式和说明，对实时的声纳图像的行数据进行DCT变换，得到DCT系数；将得到的DCT系数划分为S段，截取第一段DCT系数(第一段DCT系数对应低频段数据)。

2分段量化

在截取段落中，DCT系数仍然符合低频能量较高，高频能量较低的特点，因此可以进行分段量化，以进一步提高压缩效率。

提出的量化规则为，设待处理的信号长度为M，则对截取的第一段落的前M/8数据使用16位整型(int16)量化，其后的M/8用12位整型(int12)量化，再之后的M/4用8位整型(int8)量化，最后的M/2段落数据用4位整型(int4)量化。如图2所示。

3压缩效率

压缩效率是表征压缩效果的重要指标。本发明所涉及压缩方法的压缩效率关键有两点。一是DCT变换后的数据截取，其压缩倍数(InterceptionRatio,IR)为可设常数S。

二是分段量化带来的压缩。声纳成像处理子系统获得的图像数据通常是浮点型，以32位浮点型(float32)为例进行计算，则量化带来的压缩倍数为

$\frac{M\·\mathrm{sizeof}\left(\mathrm{float}32\right)}{\frac{M}{8}\·\mathrm{sizeof}\left(\mathrm{int}16\right)+\frac{M}{8}\·\mathrm{sizeof}\left(\mathrm{int}12\right)+\frac{M}{4}\·\mathrm{sizeof}\left(\mathrm{int}8\right)+\frac{M}{2}\·\mathrm{sizeof}\left(\mathrm{int}4\right)}=4.267$

其中sizeof()函数表示数据类型的字节数。

总的压缩倍数(CompressionRatio,CR)为

CR＝4.267·IR

例如DCT截取段落为总长的1/8时，压缩倍数为34.14.

4处理实例

图3为圆柱形目标的处理结果，图像垂直航迹方向(图中横轴方向)的长度为200m，沿航迹方向(图中竖轴方向)的长度为35m。

从处理结果可以看出，对于截断倍数IR依次为2,4,6,8,16的情形，从大尺度来看，目标在大范围背景下解压后仍能清晰显示。将目标所在局部截取出来进行对照，如图4所示。其垂直航迹方向为10m，沿航迹方向为7.5m。从局部对照可以看出，随着截断倍数的增加，目标的清晰度也在不断下降。

为了进一步对图像压缩效果进行判断，采用图像压缩客观评价中的几种标准。如下所述。

均方误差：设原始图像行信号为f(n)，恢复图像信号为g(n)，且信号长度均为N。则均方误差(MeanSquareError,MSE)定义为

$\mathrm{MSE}=\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{[g\left(n\right)-f\left(n\right)]}^2$

峰值信噪比：峰值信噪比(PeakSignaltoNoiseRatio,PSNR)定义为

$\mathrm{PSNR}=10\lg \left(\frac{\underset{n=1}{\overset{N}{\mathrm{\Σ}}}f{\left(n\right)}^2}{\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{[g\left(n\right)-f\left(n\right)]}^2}\right)=10\lg \left(\frac{\underset{n=1}{\overset{N}{\mathrm{\Σ}}}f{\left(n\right)}^2}{N\·\mathrm{MSE}}\right)$

相关系数：相关系数(CorrelationCoefficient,CC)定义为

$\mathrm{CC}=\frac{\underset{n=1}{\overset{N}{\mathrm{\Σ}}}[f\left(n\right)\·g\left(n\right)]}{\underset{n=1}{\overset{N}{\mathrm{\Σ}}}f{\left(n\right)}^2}$

平均差异：平均差异(AverageDifference,AD)定义为

$\mathrm{AD}=\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}[f\left(n\right)-g\left(n\right)]$

对图3处理结果的客观评价如下表所示。

表1图像压缩效果客观评价表

最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (8)

1.一种成像声纳实时处理系统中的图像压缩方法，所述压缩方法包含：

步骤101)对实时的声纳图像的行数据进行DCT变换，得到DCT系数；

步骤102)将得到的DCT系数划分为S段，截取第一段DCT系数；

步骤103)将截取的第一段DCT系数再划分为若干子段，然后对各子段包含的DCT系数分别进行量化压缩处理，完成图像压缩；

其中，所述的量化压缩处理的原则为：针对较低频段数据采用较多位数进行整型量化，针对较高频段数据采用较少位数进行整型量化。

2.根据权利要求1所述的成像声纳实时处理系统中的图像压缩方法，其特征在于，DCT变换的公式为：

$y\left(k\right)=\mathrm{\ω}\left(k\right)\underset{n=1}{\overset{N}{\mathrm{\Σ}}}x\left(n\right)\cos \left(\frac{\mathrm{\π}(2n-1)(k-1)}{2N}\right),k=\mathrm{1,2},...N$

其中，x(n)表示输入信号序列；N为序列的点数；ω(k)为系数；y(k)为DCT变换后得到的系数。

3.根据权利要求1所述的成像声纳实时处理系统中的图像压缩方法，其特征在于，所述步骤103)进一步包含：

步骤103-1)将截取的第一段DCT系数划分为4个子段，其中第一子段占截取段落的1/8，第二子段占截取段落的1/8，第三子段占截取段落的1/4，第四子段占截取段落的1/2；

步骤103-2)对第一子段的浮点型数据使用16位整型量化，具体的量化过程为：

设第一子段的DCT系数对应的数据为y₁(k)，首先求得y₁(k)序列的绝对值的最大值，再将整个y₁(k)序列归一化处理，其后乘以16位有符号整型量化范围的一半，再就近取整，计算公式为：

${y_1}^{\′}\left(k\right)=\mathrm{round}(\frac{y_1\left(k\right)}{\max \left(\mathrm{abs}\left(y_1\left(k\right)\right)\right)}\×2^{15})$

其中，函数abs()表示求绝对值，函数round()表示就近取整，函数max()表示求序列的最大值；

步骤103-3)采用如下三个公式分别对第二、三、四子段的数据进行量化处理，其中第二子段使用12位整型量化，第三子段使用8位整型数据量化，第四子段使用4位整型数据量化：

${y_2}^{\′}\left(k\right)=\mathrm{round}(\frac{y_2\left(k\right)}{\max \left(\mathrm{abs}\left(y_2\left(k\right)\right)\right)}\×2^{11})$

${y_3}^{\′}\left(k\right)=\mathrm{round}(\frac{y_3\left(k\right)}{\max \left(\mathrm{abs}\left(y_3\left(k\right)\right)\right)}\×2^7)$

${y_4}^{\′}\left(k\right)=\mathrm{round}(\frac{y_4\left(k\right)}{\max \left(\mathrm{abs}\left(y_4\left(k\right)\right)\right)}\×2^3)$

其中，y_i(k)表示划分得到的第i个子段，其中i＝1,2,3,4，y_i'(k)表示对第i子段量化后得到的序列。

4.根据权利要求1所述的成像声纳实时处理系统中的图像压缩方法，其特征在于，所述的段数S应满足如下公式：

$\mathrm{\ρ}\left(S\right)=\frac{\underset{k=1}{\overset{k=\frac{N}{S}}{\mathrm{\Σ}}}{y}^2\left(k\right)}{\underset{k=1}{\overset{k=N}{\mathrm{\Σ}}}{y}^2\left(k\right)}$

ρ(S)＞ρ₀

其中，ρ(S)为前1/S段DCT系数的均方和与整个DCT系数序列的均方和的比值，ρ₀为设置的能量含量阈值。

5.一种成像声纳实时处理系统中的图像解压方法，所述的解压方法用于解压采用权利要求3记载的压缩方法的数据，所述解压方法包含：

步骤201)根据反量化表达式对接收的数据进行反量化，对第一子段，第二子段、第三子段和第四个子段的数据的反量化公式分别如下：

$\left\{\begin{array}{c}{\tilde{y}}_1\left(k\right)=\frac{{y_1}^{\′}\left(k\right)}{2^{15}}\×\max \left(\mathrm{abs}\left(y_1\left(k\right)\right)\right)\\ {\tilde{y}}_2\left(k\right)=\frac{{y_2}^{\′}\left(k\right)}{2^{11}}\×\max \left(\mathrm{abs}\left(y_2\left(k\right)\right)\right)\\ {\tilde{y}}_3\left(k\right)=\frac{{y_3}^{\′}\left(k\right)}{2^7}\×\max \left(\mathrm{abs}\left(y_3\left(k\right)\right)\right)\\ {\tilde{y}}_4\left(k\right)=\frac{{y_4}^{\′}\left(k\right)}{2^3}\×\max \left(\mathrm{abs}\left(y_4\left(k\right)\right)\right)\end{array}\right.$

将反量化后得到的四段数据连接为一段数据，公式为：

$\tilde{y}\left(k\right)=\left[\begin{array}{cccc}{\tilde{y}}_1\left(k\right)& {\tilde{y}}_2\left(k\right)& {\tilde{y}}_3\left(4\right)& {\tilde{y}}_4\left(k\right)\end{array}\right]$

步骤201)将反量化后的数据进行DCT反变换，得到实际的图像数据；所述的DCT反变换为：

$x\left(n\right)=\underset{k=1}{\overset{N}{\mathrm{\Σ}}}\mathrm{\ω}\left(k\right)y\left(k\right)\cos \left(\frac{\mathrm{\π}(2n-1)(k-1)}{2N}\right),n=\mathrm{1,2},...N$

其中

$\mathrm{\ω}\left(k\right)=\left\{\begin{array}{c}\frac{1}{\sqrt{N}},k=1\\ \sqrt{\frac{2}{N}},2\≤k\≤N\end{array}\right.$

其中，N为对实时的声纳图像的行数据进行DCT变换得到的DCT系数的总长度。

6.一种用于成像声纳的图像压缩及解压系统，其特征在于，所述系统压缩子系统和解压子系统，所述压缩子系统包含：

DCT变换模块，用于对实时的声纳图像的行数据进行DCT变换，得到DCT系数；

截取模块，用于将得到的DCT系数划分为S段，截取第一段DCT系数；

分段量化处理模块，用于将截取的第一段DCT系数再划分为若干子段，然后对各子段包含的DCT系数分别进行量化压缩处理，完成图像压缩；

所述解缩子系统包含：

反量化处理模块，用于根据反量化表达式对接收的数据进行反量化，对第一子段，第二子段、第三子段和第四个子段的数据的反量化公式分别如下：

$\left\{\begin{array}{c}{\tilde{y}}_1\left(k\right)=\frac{{y_1}^{\′}\left(k\right)}{2^{15}}\×\max \left(\mathrm{abs}\left(y_1\left(k\right)\right)\right)\\ {\tilde{y}}_2\left(k\right)=\frac{{y_2}^{\′}\left(k\right)}{2^{11}}\×\max \left(\mathrm{abs}\left(y_2\left(k\right)\right)\right)\\ {\tilde{y}}_3\left(k\right)=\frac{{y_3}^{\′}\left(k\right)}{2^7}\×\max \left(\mathrm{abs}\left(y_3\left(k\right)\right)\right)\\ {\tilde{y}}_4\left(k\right)=\frac{{y_4}^{\′}\left(k\right)}{2^3}\×\max \left(\mathrm{abs}\left(y_4\left(k\right)\right)\right)\end{array}\right.$

将反量化后得到的四段数据连接为一段数据，公式为：

$\tilde{y}\left(k\right)=\left[\begin{array}{cccc}{\tilde{y}}_1\left(k\right)& {\tilde{y}}_2\left(k\right)& {\tilde{y}}_3\left(4\right)& {\tilde{y}}_4\left(k\right)\end{array}\right]$

DCT反变换，用于将反量化后的数据进行DCT反变换，得到实际的图像数据；所述的DCT反变换为：

$x\left(n\right)=\underset{k=1}{\overset{N}{\mathrm{\Σ}}}\mathrm{\ω}\left(k\right)y\left(k\right)\cos \left(\frac{\mathrm{\π}(2n-1)(k-1)}{2N}\right),n=\mathrm{1,2},...N$

其中

$\mathrm{\ω}\left(k\right)=\left\{\begin{array}{c}\frac{1}{\sqrt{N}},k=1\\ \sqrt{\frac{2}{N}},2\≤k\≤N\end{array}\right.$

其中，N为对实时的声纳图像的行数据进行DCT变换得到的DCT系数的总长度。

7.根据权利要求6所述的用于成像声纳的图像压缩解压系统，其特征在于，所述分段量化处理模块进一步包含：

分段子模块，根据DCT变换的能量集中特性，进一步将截取的第一段落划分为四个子段，第一子段占第一段落总长度的1/8，第二子段占第一段落总长度1/8，第三子段占第一段落总长度的1/4，第四子段占第一段落总长度的1/2；

量化子模块，使用不同精度的整型数据分别对各子段进行量化处理，其中第一子段使用16位整型量化，第二子段使用12位整型量化，第三子段使用8位整型量化，第四子段使用4位整型量化。

8.根据权利要求6所述的成像声纳实时处理系统中的图像压缩解压系统，其特征在于，所述的段数S应满足下式：

$\mathrm{\ρ}\left(S\right)=\frac{\underset{k=1}{\overset{k=\frac{N}{S}}{\mathrm{\Σ}}}{y}^2\left(k\right)}{\underset{k=1}{\overset{k=N}{\mathrm{\Σ}}}{y}^2\left(k\right)}$

ρ(S)＞ρ₀

其中，ρ(S)为前1/S段DCT系数的均方和与整个DCT系数序列的均方和的比值，ρ₀为设置的能量含量阈值。