一种基于改进ezw的声纳图像数据压缩方法
【专利摘要】本发明公开了一种基于改进EZW的声纳图像数据压缩方法，包括步骤：1)采集声纳图像信号，分析其统计特性和能量分布特点；2)对声纳图像进行整数提升小波变换，建立小波系数树结构；3)将小波变换后图像的低频子带经过调整后直接存储；4)将高频部分应用改进的嵌入式零树编码(EZW)；5)合并低高频数据，实现声纳数据的压缩。本发明的优点在于：采用小波变换对声纳图像信号进行处理，使其幅值分布相对集中，更利于压缩编码；采用整数提升小波变换进行声纳图像处理，可有效地提高运算效率，实现图像的无损压缩；压缩声纳数据其压缩速度和PSNR值均优于传统的EZW算法。
【专利说明】一种基于改进EZW的声纳图像数据压缩方法
【技术领域】
[0001]本发明涉及涉及图像压缩方法，尤其是涉及一种基于改进EZW的声纳图像数据压缩方法。
【背景技术】
[0002]声纳是水下目标探测的重要工具。获取声纳数据的主要方法是:将声纳探测到的数据通过无线或有线通信介质实时传输回水面。这种方式的灵活性和实时性都非常好，但是，由于信道容量，尤其是无线信道容量非常有限，而声纳数据量通常又较大，因此，对声纳数据进行压缩就非常必要。
[0003]目前，声纳数据压缩的方法有很多，主要方法有:基于DCT变换的压缩方法和基于小波变换的压缩方法。传统的基于DCT变换的图像压缩虽然在较高码率下能够提供较好的图像质量，但在码率低于0.25bpp时，重构图像存在严重的方块效应；同时，基于DCT变换的图像压缩不能在单一码流中实现图像的有损和无损压缩，从而不能实现从有损到无损的累进式传输。基于小波变换的压缩方法的提出有效地克服了以上缺点。由Wim Sweldens等人提出的整数提升的小波变换通过简单地分裂、预测和更新等步骤来完成对一列数字信号的变换，是整数到整数的小波变换，有效地提高运算效率。
[0004]Shapr i ο根据小波分解后同方向子带中的小波系数存在的相似性，利用一种称为小波树的树形结构来组织这些小波系数，设计了嵌入式零树小波编码方法(EmbeddedZero-tree Wavelet, EZff)，有效地利用了小波系数的特性，是至今最有效的小波编码方法之一。但该方法在压缩声纳数据上存在一些缺点，主要表现在:①经过小波变换后，声纳图像的大部分能量集中在低频子带，低频子带的编码几乎不影响整个图像的压缩比，但其较小损失就可能对恢复图像的质量造成较大的影响。而EZW将低频数据和高频数据用同样的方法进行编码，在较低码率的情况下，低频子带的信息损失较多，难以保证恢复图像的质量。②在量化过程中，门限T的值按2的负幂级逐步减小，每减小一次T，就需要对整个高频子带重新做零树分类。对于每一次扫描，得到的重要系数的绝对误差上限为T/4，对于重要系数不是很多的声纳图像来说，每一轮的扫描输出较多是表示零树根和孤立零的码字，而表示重要系数的码字则相对较少，因此在输出的码流中，大部分的复原数据为0，这些O值对复原图像是无效的。③在每次主扫描后，要将重要系数的相关信息存储在辅表中，在辅扫描的过程中对辅表中的元素逐个细化编码。若能在主扫描中对重要系数直接细化编码，则将有效地节省时间内存。
[0005]以上研究虽然对经典小波压缩算法进行了改进，但由于压缩结果与图像本身特点关系密切，不同类型图像其统计特性不同，即使同一类图像统计特性也有差别，故其压缩算法不仅要能针对特定类型的图像进行，而且还需要能够在一定范围内适应统计特性的变化。针对声纳图像数据的特点，如何用声纳获得图像数据并进行高效率的压缩，是本专利介绍的内容。
[0006]综上所述，由于现有技术存在不足，就需一种高效率的声纳图像数据压缩方法。
【发明内容】
[0007]本发明主要是解决现有技术所存在的技术问题；提供了一种采用小波变换对声纳图像信号进行处理，使其幅值分布相对集中，更利于压缩编码的一种基于改进EZW的声纳图像数据压缩方法。
[0008]本发明还有一目的是解决现有技术所存在的技术问题；提供了一种采用整数提升小波变换进行声纳图像处理，是整数到整数的小波变换，可有效地提高运算效率，实现图像的无损压的一种基于改进EZW的声纳图像数据压缩方法。
[0009]本发明的上述技术问题主要是通过下述技术方案得以解决的:
[0010]一种基于改进EZW的声纳图像数据压缩方法，其特征在于，包括以下步骤:
[0011]步骤1、采集声纳图像信号，并采用提升格式的LeGall5，3小波对声纳图像进行小波变换，建立小波系数树结构；提升小波变换分为分裂、预测和更新三个步骤，LeGall小波整数实现形式如下:
[0012]正变换:
[0013]
【权利要求】
1.一种基于改进EZW的声纳图像数据压缩方法，其特征在于，包括以下步骤: 步骤1、采集声纳图像信号，并采用提升格式的LeGall5，3小波对声纳图像进行小波变换，建立小波系数树结构；提升小波变换分为分裂、预测和更新三个步骤，LeGall小波整数实现形式如下: 正变换:
式中符号L」表示取整运算，Xext表示周期对称延拓后的信号； 按照LeGall小波变换算法，完成一次对图像水平方向上的提升小波变换，得到水平方向上低频L和高频H两个部分； 接下来用同样的方法分别再对这两个部分进行垂直方向上的提升小波变换小波，得到LL, LH, HL, HH ;这样整个2维的提升小波变换就完成了 ；到这里，完成的是一级的2维提升小波变换，对小波变换后的低频部分再做小波变换，如此循环N次，就得到N级的2维小波变换； 经过N级小波变换的小波图像，对于低频子图中的某一系数而言，与其对应的具有相同空间定位的高频子图中的系数称为是它的子孙，从图像的低频层开始依照子孙关系延伸，得到树形结构； 步骤2、将步骤I中得到的高频与低频数据分开编码；若为低频子带的数据，按公式
将低频数据映射到[0，255]之间后直接存储，存储后转至步骤3.5 ;若为高频子带的数据，则进行步骤3 ;式中，c为小波系数，Min为小波系数的最小值，Max为小波系数的最大值，f(c)为映射后的值； 步骤3、将高频部分应用改进的嵌入式零树编码，具体包括如下子步骤: 步骤3.1、初始化阈值
ci，j为小波系数； 步骤3.2、改进的主扫描过程:按“Z”字形扫描，对于不同的系数类型，做不同的处理:选择处理一:如果为正重要系数:输出符号POS，再根据其绝对值输出幅值码:如果在区间[T，T+T/4)则输出OO ;如果在区间[Τ+Τ/4，Τ+Τ/2)输出01 ;如果在区间[Τ+Τ/2，Τ+3Τ/4)输出10 ;如果在区间[Τ+3Τ/4，2Τ)输出11 ； 选择处理二:如果为负重要系数:输出符号NEG，再根据其绝对值输出幅值码:如果在区间[Τ，Τ+Τ/4)输出00 ;如果在区间[Τ+Τ/4，Τ+Τ/2)输出01 ;如果在区间[Τ+Τ/2，Τ+3Τ/4)输出10 ;如果在区间[T+3T/4，2T)输出11 ； 选择处理三:如果为孤立零点:输出符号IZ ； 选择处理四:如果为零树根:输出符号ZTR ； 步骤3.3、重新设置阈值T = Τ/2，如果T= I或者达到码率要求，算法终止，否则转至步骤3.2 ； 步骤3.4、对主扫描输出的类型码和幅值码进行自适应二进制编码； 步骤3.5、输出高频数据码流，编码过程结束； 步骤4、将低高频数据合并，输出码流，实现声纳数据的编码压缩； 步骤5、解码为编码的逆过程，首先对得到的码流进行算术解码；对于最低频子带的编码码流按公式
行调整，就得到相应的最低频子带小波变换重建系数；对于其他高频子带的码流，解码器利用接收到的编码器发送过来的相关信息，设置相应的阈值，进行主扫描解码:在解码出POS、NEG、ZTR、IZ类型码之后，如果是重要系数POS或者NEG，则直接进行幅值码的解码，即根据区间码进行重要系数的重构；“00”重构为.1.1251\，“01”重构为1.375?\，“10”重构为1.725Ti; “ 11”重构为1.875^ ;阈值减半，重复主扫描过程，直到阈值为I或者达到码率要求，结束解码；最后进行(5，3)逆整数提升小波变换，即可转换为原始声纳图像了。
【文档编号】H04N19/96GK104079947SQ201410291450
【公开日】2014年10月1日 申请日期:2014年6月25日 优先权日:2014年6月25日
【发明者】饶云华, 曾敏 申请人:武汉大学