基于改进SPIHT的医学图像压缩.docVIP

基于改进SPIHT的医学图像压缩 　　【摘要】为了满足医学数据的存储和传输的需要，在SPIHT算法的基础上，提出了一种基于改进SPIHT的医学图像压缩算法，实现了对医学图像的编码。该算法在减少编码时间的同时，能够得到较高的压缩比，重建图像的峰值信噪比也有明显的提高。选取几幅医学图像进行有效性的验证，试验结果表明：此算法可以更好地保留图像的客观质量，提高了编解码的效率，获得了较好的压缩效果。 　　【关键词】医学图像压缩；改进SPIHT算法；小波变换；峰值信噪比 　　1.引言 　　随着远程医疗系统的逐步成熟和发展，对于数据压缩和通信标准都提出了严格要求。 　　医学图像压缩的目的在于保留原始医学图像数据中蕴含的真实信息前提下，同时用较少的比特数据来表征医学图像，从而减少医学图像数据的存储空间和远程医疗中的传输时间，最终满足海量医学图像数据有效存储和实时传输的实际应用需求。 　　小波变换是近年来发展起来的一种新的时频分析方法，它在时域和频域同时具有良好的局部化特征，特别适合对ECG这类非平稳信号进行分析和数据压缩[2]。1993年，J.M.Shapiro利用小波系数之间的自相似性，提出了嵌入式零树编码算法[3]（Embedded Zerotree Wavelet，EZW）。EZW有很好的均方误差性能、适中的复杂度和嵌入式位流，引起了压缩研究界的兴趣。1996年，由A.Said和W.A.Pearlman提出的多级树集合分裂算法[4]（Set Partitioning in Hierarchical Tree.SPIHT）是对EZW算法的改进，但是具有更高的压缩效率，在相同的比特率下，即使不需要熵编码环节，SPIHT[5]也总能取得优于EZW的量化性能。本文结合小波变换和改进SPIHT算法来实现医学图像的编码，从而达到较优的压缩重建效果。 　　2.改进的SPIHT算法 　　对于SPIHT算法，该算法需要不断重复地对小波系数进行扫描，与其阈值进行比较，才能找到重要的系数，从而提高了访问内存的次数，增加了内存的负荷；同时在进行LSP，LIP分类扫描时，对D型态和L型态进行编码时，存在着对部分重要小波系数进行冗余编码的情况，大大降低了编解码的效率。改进的SPIHT算法可以有效解决这些问题。 　　2.1 调整小波系数的大小顺序 　　将小波系数由大到小进行排序，将排序好的小波系数矩阵存放到一维数组A当中，并且把相对应的小波系数的坐标存放到一维数组B当中，这样的调整使得排序的次数大大减少，从而在进行编码解码时，减少了对小波系数的扫描次数和内存的访问次数，提高编解码效率。 　　2.2 改变集合分布关系 　　SPIHT算法采用的各集合关系为： 　　（1） 　　式中，表示结点和所有后代的坐标集合；表示结点的所有后代的坐标集合；表示结点的直接后代的坐标集合；表示结点除直接后代以外所有后代的坐标集合。 　　由于同一个高频上的第一个分量和其他三个分量存在着很强的相关性，所以在SPIHT算法中，在判定编码是否有效时，可以只对无效集合的第一个分量进行判定，如果其值小于阈值，则此集合判定为无效，这样能减少对无效集合扫描的次数，但也有可能导致误判，最终影响图像的重构效果。改进的SPIHT算法对初始节点后续N代内子孙进行全扫描，而对N代后的子孙只对第1分量进行扫描，（N的取值是由小波系数的能量分布决定的，能量越集中，N值越小），这样能够平衡编码的性能，从而提高编码的效率。 　　改进SPIHT算法的各集合关系为： 　　（2） 　　式中，表示结点的N代以内子孙的集合；表示结点的N代以外子孙的集合；表示子带第1分量上N代以外子孙的集合；表示子带第2，3，4分量上N代以外子孙的集合。 　　2.3 一维数组存放小波系数 　　在编码过程中，SPIHT使用了三个链表：不重要系数表LIP（the list of insignifi-cant pixels），重要系数表LSP（the list of significant pixels），和不重要集合表LIS（the list of insignificant sets）。这三个表与四叉树相对应，当链表过长时，所需内存会增大，这样会使得编解码速率大大降低。 　　改进的SPIHT算法采用一维数组来存放已排序好的小波系数，在扫描小波系数时，小于阈值的小波系数不加以考虑，提高了编解码的速率。 　　2.4 消除部分重要小波系数的冗余编码 　　消除冗余编码的流程图如图1所示。 　　从图1可知：在改进的SPIHT算法中，当判断4个O（a，b）和前三个D（a，b）都不重要时，则不需对L（i，j）和最后一个D（a，b）扫描，就能判断出它的重要性。而在SPIHT算法中，则必须对其进行扫描，这样就会产生