本发明专利技术公开了一种基于SPIHT和DWT算法的医学图像压缩方法，包括如下步骤：(1)输入原医学图像；(2)采用小波滤波器对输入图像执行DWT操作；(3)对于经过DWT的小波信号进行熵编码处理，采用SPIHT算法执行压缩处理；(4)对于经过压缩后的比特流SPIHIT解码操作；(5)执行逆DWT处理，生成最终的压缩图像。本发明专利技术的优点在于：基于SPIHT和DWT算法的医学图像压缩方法获得的压缩图像效果要明显优于传统的图像算法。本专利提出的压缩技术在相同压缩比的前提下，具有较低的MSE和较高的PSNR，进一步提高了医学图像压缩算法的质量。图像压缩算法的质量。图像压缩算法的质量。
【技术实现步骤摘要】
一种基于SPIHT和DWT算法的医学图像压缩方法及装置
[0001]本专利技术涉及医学影像领域，特别涉及一种基于SPIHT和DWT算法的医学图像压缩方法。
技术介绍
[0002]由于医学图像的数据量非常大，所需要的存储空间大、传输频带宽，严重地影响了PACS的发展。对于医学图像的压缩技术研发一直是PACS系统的热门研究方向之一，研究在保持高图像质量的前提下提升图像压缩的压缩率在医学图像邻域具有极高的研究价值。JPEG2000是基于小波变换的图像压缩标准，与传统的JPEG算法相比压缩比更高，而且不会产生基于离散余弦变换的JPEG标准产生的块状模糊瑕疵。JPEG2000由于其特有的无损压缩模式常被用于医学图像压缩，然而JPEG2000存在着压缩比低和压缩效率低等缺点。如专利申请号为201210260832.8的一种JPEG图像压缩方法和系统，其就无法解决上述问题，无法应用于医学图像的压缩。
[0003]对于医学图像的来说，图像的细节也可能具有重要的诊断意义，细节部分的失真，可能造成误诊与误判会产生严重的后果。传统的图像压缩方法在高压缩率下容易造成医学图像细节部分存在失真等问题，因此传统的图像压缩方法无法适用于医学图像。
技术实现思路
[0004]本专利技术的目的在于克服现有技术的不足，提供一种基于SPIHT和DWT算法的医学图像压缩方法，解决现有技术无法适用于医学图像的缺陷，提出的图像压缩方法具有提高医学图像重要部分重建质量的能力，可以在保持较高压缩率的前提下保护好图像的细节部分。
[0005]为了实现上述目的，本专利技术采用的技术方案为：一种基于SPIHT和DWT算法的医学图像压缩方法，包括如下步骤：
[0006](1)输入原医学图像；
[0007](2)采用小波滤波器对输入图像执行DWT操作；
[0008](3)对于经过DWT的小波信号进行熵编码处理，采用SPIHT算法执行压缩处理；
[0009](4)对于经过压缩后的比特流SPIHT解码操作；
[0010](5)执行逆DWT处理，生成最终的压缩图像。
[0011]步骤(2)中采用9
‑
7tap小波滤波器对输入图像执行DWT操作。
[0012]SPHIT算法利用空间方向树来制定划分子集的集合分割规则，使用了3个链表来组织空间方向树中的结构：
[0013]重要像素链表LSP、不重要像素链表LIP、不重要集合链表LIS；
[0014]在LSP和LIP中每一个(i,j)都代表一个像素，但是在LIS中(i,j)代表的是一个所有后代D(i,j)集合或所有间接后代L(i,j)集合。
[0015]采用SPIHT算法执行压缩处理包括:
[0016](1)初始化步骤；
[0017](2)排序步骤；
[0018](3)细化步骤；
[0019](4)更新步骤；
[0020](5)对LSP中的值进行编码。
[0021]H为所有空间树结构的根节点，即整个子带LL2。
[0022]初始化步骤包括设置初始阈值T
n
＝2
n
，其中其中c
ij
为任意小波系数；将LSP置空，将H中所有根节点加入LIP，即LIP＝{(i,j)|(i.j)∈H}；然后把根节点H中所有具有子带的节点加入LIS，同时以TypeA对其标记。
[0023]排序步骤是对LIP和LIS表中的重要系数进行判断并将重要系数放入LSP的过程。
[0024]细化步骤是将重要系数输出的过程，最后更新步长。
[0025]对LSP中的值进行编码，将每个值变换成二进制，然后取其第n位为编码值，其中n＝logT式中，T表示阈值；对所有的LIP、LIS以及LSP中的值编码完毕后，将阈值减小一倍，再进行依次编码，直到阈值变为1，则编码完毕。
[0026]一种基于SPIHT和DWT算法的医学图像压缩装置，包括医学图像输入单元、小波滤波器、DWT处理单元、SPIHT压缩处理单元、解码单元、逆DWT单元，医学图像输入单元用于输入原医学图像，其输出端连接小波滤波器，采用小波滤波器对输入图像执行DWT操作；小波滤波器的输出端连接至SPIHT压缩处理单元，其对于经过DWT的小波信号进行熵编码处理，采用SPIHT算法执行压缩处理，其输出端连接至解码单元，所述解码单元对于经过压缩后的比特流SPIHIT解码操作，其输出端连接至逆DWT单元，其用于执行逆DWT处理，生成最终的压缩图像。
[0027]本专利技术的优点在于：基于SPIHT和DWT算法的医学图像压缩方法获得的压缩图像效果要明显优于传统的图像算法。本专利提出的压缩技术在相同压缩比的前提下，具有较低的MSE和较高的PSNR，进一步提高了医学图像压缩算法的质量。
附图说明
[0028]下面对本专利技术说明书各幅附图表达的内容及图中的标记作简要说明：
[0029]图1为本专利技术压缩方法流程图；
[0030]图2为本专利技术处理LIP表中的每一个节点处理方法流程图；
[0031]图3为本专利技术处理LIS表中的TypeA节点的处理方法流程图；
[0032]图4为本专利技术处理LIS表中的TypeB节点的处理方法流程图；
[0033]图5为经过二维离散小波变换的SPIHT算法中的空间树结构。
具体实施方式
[0034]下面对照附图，通过对最优实施例的描述，对本专利技术的具体实施方式作进一步详细的说明。
[0035]本专利技术中提出的图像压缩方法具有提高医学图像重要部分重建质量的能力，可以在保持较高压缩率的前提下保护好图像的细节部分，具体方法如下：
[0036]1.输入原医学图像；
[0037]2.采用9
‑
7tap小波滤波器对输入图像执行DWT操作；
[0038]3.对于经过DWT的小波信号进行熵编码处理，采用本文提出的SPIHT算法替换原有的哈夫曼编码执行压缩处理；
[0039]4.对于经过压缩后的比特流SPIHIT解码操作；
[0040]5.执行逆DWT处理，生成最终的压缩图像。
[0041]步骤3中SPIHT算法过程步骤：
[0042]SPHIT算法利用空间方向树来制定划分子集的集合分割规则，使用了3个链表来组织空间方向树中的结构：
[0043](1)重要像素链表(LSP)(List of Significant Pixels)；
[0044](2)不重要像素链表(LIP)(List of Insignificant Pixels)；
[0045](3)不重要集合链表(LIS)(List of Insignificant Sets)；
[0046]在LSP和LIP中每一个(i,j)都代表一个像素，但是在LIS中(i,j)代表的可能是一个所有后代D(i,j)集合(称之为typeA)或所有间接后代L(i,j)集合(称之为typeB)。...
【技术保护点】
【技术特征摘要】
1.一种基于SPIHT和DWT算法的医学图像压缩方法，其特征在于：包括如下步骤：(1)输入原医学图像；(2)采用小波滤波器对输入图像执行DWT操作；(3)对于经过DWT的小波信号进行熵编码处理，采用SPIHT算法执行压缩处理；(4)对于经过压缩后的比特流SPIHIT解码操作；(5)执行逆DWT处理，生成最终的压缩图像。2.如权利要求1所述的一种基于SPIHT和DWT算法的医学图像压缩方法，其特征在于：步骤(2)中采用9
‑
7tap小波滤波器对输入图像执行DWT操作。3.如权利要求1所述的一种基于SPIHT和DWT算法的医学图像压缩方法，其特征在于：SPHIT算法利用空间方向树来制定划分子集的集合分割规则，使用了3个链表来组织空间方向树中的结构：重要像素链表LSP、不重要像素链表LIP、不重要集合链表；在LSP和LIP中每一个(i,j)都代表一个像素，但是在LIS中(i,j)代表的是一个所有后代D(i,j)集合或所有间接后代L(i,j)集合。4.如权利要求1
‑
3任一所述的一种基于SPIHT和DWT算法的医学图像压缩方法，其特征在于：采用SPIHT算法执行压缩处理包括:(1)初始化步骤；(2)排序步骤；(3)细化步骤；(4)更新步骤；(5)对LSP中的值进行编码。5.如权利要求1或2所述的一种基于SPIHT和DWT算法的医学图像压缩方法，其特征在于：初始化步骤包括设置初始阈值T
n
＝2
n
，其中其中c
ij
【专利技术属性】
技术研发人员：任鹏，刘志恒，赵立军，曹雏清，
申请(专利权)人：哈尔滨工业大学芜湖机器人产业技术研究院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。