基于小波变换的医学图像压缩技术的研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
随着现代医院数字化进程的加快,数字化医学图像在医院的诊断和治疗中发挥着越来越重要的作用。数字化医学图像的广泛应用,产生了大量的数据,给图像的存储和传输技术提出了严峻的挑战。因此,研究有针对性的医学图像压缩技术具有重要的实际意义。
     根据医学图像的特点和医学图像对压缩编码方案的特殊要求,本文基于小波变换,嵌入式图像编码,图像分割和感兴趣区域编码的各种理论方法,以从山东省立医院采集到的CT图像,MR图像和超声图像为主要研究对象,结合临床的实际需求,重点解决医学图像高图像质量和高压缩比之间的矛盾,并通过matlab编程对研究过程中涉及的算法进行了实现和分析。论文的主要工作如下:
     (1)研究整数小波变换在医学图像压缩中的应用。由于传统的小波变换需要进行小波系数的舍入运算,无法重建原始无损的图像,本文通过对整数小波变换理论和医学图像特征的分析研究,采用既能实现无损压缩又能实现有损压缩的5/3双正交小波基,通过小波分解后各子带的信息熵值和计算量来综合确定小波变换的级数,采用周期对称延拓方式实现了提升方案下的三级整数小波变换,使能量集中到少数小波系数上。实验结果表明利用整数小波变换进行图像压缩可以避免计算机的舍入误差,有利于医学图像的无损压缩。通过实验方法分析了图像小波系数的分布特点,为选择合适的压缩编码方案奠定了基础。
     (2)研究嵌入式医学图像压缩编码算法。通过对图像小波变换后小波系数分布特点的分析,本文发现它非常适合于嵌入式图像编码算法。在深入分析嵌入式零树小波编码算法和多级树集合分裂(SPIHT)编码算法优缺点的基础上,将SPIHT编码和自适应算术编码相结合,提出了基于整数小波变换的静态医学图像压缩编码方案,并利用此方案对四种不同类型的医学图像(细胞图像,CT图像,MR图像和超声图像)进行了压缩实验,给出了与其它编码方法的比较结果,仿真结果验证了该编码方案对医学图像压缩的有效性。
     (3)研究基于感兴趣区域的医学图像压缩编码算法。对于医学图像,医生往往仅对病变区域感兴趣,为了更好的解决医学图像高图像质量和高压缩比之间的矛盾,本文在提出的静态医学图像压缩编码方案中引进感兴趣区域编码,设计了更具实用性和兼容性的医学图像感兴趣区域编码方案。将一幅医学图像划分为感兴趣区域(病变区域)和背景区域两部分,针对JPEG 2000标准中一般位移法和最大位移法的优缺点,本文提出将图像分割与图像压缩相结合,采用主动轮廓模型算法提取病变区域作为感兴趣区域,并用回溯法生成掩模,解决医学图像压缩技术中感兴趣区域自动提取的问题。采用最大位移法进行了压缩实验,并对实验结果进行了分析。另外,考虑到医学图像周围的相关组织和背景也为正确诊断提供了参照,将能够灵活调整感兴趣区域和背景区域重要性的通用的部分重要位平面偏移方法应用到医学图像压缩中,并将实验结果与最大位移法进行了比较。实验结果表明,利用此方法既可以保证诊断信息有良好的视觉效果,同时又克服了最大位移法的缺点,在提高压缩率的基础上,有效实现了对感兴趣区域和背景区域压缩质量的灵活调整,为医学图像压缩方法提供了一种新思路。
     论文最后对全文进行了总结,并对下一步的研究工作进行了展望。
With the acceleration of digital process in modern hospital, digital medical images are playing an increasingly important role in the diagnostic and treatment of hospital. Large amount of data are produced in the widespread application of medical image and they bring tough challenges to the storage and transmission of medical image. So it is of great practical significance to research targeted medical image compression techniques.
     According to the characteristics of medical image and the special needs of medical image compression scheme, based on theories and methods of wavelet transform, embedded image coding, image segmentation and region of interest coding, CT, MR and ultrasound images obtained from Shandong Provincial Hospital were mainly studied in this paper. Combining with clinical requirement, this paper mainly solved contradiction between high image quality and high compression ratio of medical image and the algorithms involved in research process were realized and analyzed by matlab programming. The main contents are listed below:
     (1) Research on the application of integer wavelet transform in medical image compression. Because traditional wavelet transform need to make rounding operation for wavelet coefficients, it can not reconstruct original lossless images. By analyzing integer wavelet transform theories and medical image characteristics, 5/3 biorthogonal wavelet basis, which can realize either lossless compression or lossy compression, is adopted in this paper. According to the information entropy value of each sub-band and computational complexity, the level of wavelet decomposition is comprehensively determined. Symmetric periodic extension is chose to realize three level integer wavelet transform based on lifting scheme and image energy concentrate on few wavelet coefficients. Experimental results show that the application of integer wavelet transform in image compression can avoid rounding error of computer and it is beneficial to medical image lossless compression. The distribution of wavelet coefficients lays the foundation for suitable coding scheme.
     (2)Research on coding algorithms of embedded medical image. According to the distribution characteristics of wavelet coefficients, it is very suitable for embedded image coding algorithm. On the basis of thorough analysis of advantage and disadvantage for embedded zero tree wavelet coding and set partition in hierachical trees coding, by combing SPIHT coding and adaptive arithmetic coding, a still medical image coding scheme based on integer wavelet transform is proposed in this paper. The compressive experiment to four different types of medical image (cell image, CT image, MR image and ultrasound image) is conducted and comparison results with other coding methods show that the scheme is very effective for medical image compression.
     (3)Research on medical image coding algorithms based on region of interest. Doctors are always only interested in lesion areas of medical image. In order to solve contradiction between high image quality and high compression ratio of medical image better, region of interest coding is introduced in the proposed still medical image coding scheme and ROI coding scheme is designed. Medical image is divided into region of interest area (lesion area) and background area. According to the advantage and disadvantage of Maxshift method and general scaling based method in JPEG 2000 standard, the method of combing image segmentation and image compression is presented in this paper. Active contour model is used to extract ROI region, backtracking is adopted to generate mask to solve the problem of ROI automatic extraction in medical image compression techniques. Maxshift method is adopted and experimental results are analyzed. In addition, considering that related tissue and background of medical image also provide reference for accurate diagnosis, Generalized Partial Bitplanes Shift Method is applied to medical image compression, which can adjust the significance of ROI and background region flexibly. Compared with Maxshift method, experimental results show that this method can guarantee the perfect visual effect of diagnostic information and can overcome the disadvantage of Maxshift. On the basis of improving compression rate, it can realize the flexible adjustment of image quality between ROI area and background area effectively and provide a new method for medical image compression.
     Finally, a conclusion of the main contents is made and the future research directions in this field are proposed.
引文
[1]贾克斌.数字医学图像处理、存档及传输技术[M].北京:科学出版社,2006.
    [2]冯前进.基于小波变换的医学图像编码[A].广州:第一军医大学,2003,5.
    [3]袁志勇,李敏,陈绵云.基于整型小波变换和嵌入式零树编码的医学图像压缩[J].计算机应用研究,2004,(6):109-111.
    [4]吴乐南.数据压缩[M].北京:电子工业出版社,2000.
    [5]薛玲.医学影像压缩技术研究及应用[D]_长春:长春理工大学,2007,6.
    [6]http:www.bioisland.com/ibidb/html/text.asp?key=944,2005,3.
    [7]William A.Pearlman.Lossless volumetric medical image compression[C].Proceedings of SPIE-The International Society for Optical Engineering,1999:305-312.
    [8]杨胜天,童勤业.生物医学信号的近无损压缩[J].中国生物医学工程学报,2003,22(7):235-240.
    [9]顾逶迤,高见,鲍旭东.基于人类视觉特性的医学图像压缩模型[J].中国图像图形学报,2003,8(10):1189-1195.
    [10]王相海,张福炎.静态图像编码研究进展[J].计算机研究与发展,2001,38(11):1315-1326.
    [11]Kunt M,et al.Second-generation image-coding techniques[J].IEEE,1985,73(4):549-574.
    [12]Barnsley M F,Sloan A D.A better way to compress images[J].Byte,1988,13(1):215-223.
    [13]Mallat S G.A theory for multi-resolution signal decomposition:The wavelet representation[J].IEEE Transactions On Pattern Analysis and Machine Intelligence,1989,11(7):674-693.
    [14]Min-Mo Sung,Hee-Joung Kim.Clinical Evaluation of JPEG 2000 Compression Algorithm for Digital Mammography[C].IEEE,2002:1380-1383.
    [15]N.R.Wanigasekara,Shenghao Ding and Zhuangzhi Yan.Quality Evaluation For JPEG Based On Medical Image Compression[C].IEEE Proceedings of the Second Joint EMBS/BMES Conference,Huston,USA,2002,10:1019-1020.
    [16]Kazuhiko Hamamoto.Study On Medical Ultrasonic Echo Image Compression By JPEG 2000-Optimization And The Subjective Assessment Of The Quality [C].Proceedings of the 25~(th) Annual International Conference of the IEEE EMBS,2003,9:833-836.
    [17]郑杰,樊来耀.一种JPEG 2000感兴趣区域编码方法一中央凹移位数字[J].电视技术,2003,5:7-9.
    [18]张立保,王坷.一种基于感兴趣区的小波图像编码算法[J].计算机工程与应用,2003,4:25-28.
    [19]刘杰.医学图像压缩算法研究及进展[J].北京生物医学工程,2006,25(3):330-332.
    [20]Pratt.Digital image processing[M].New York:2~(nd) Edition wiley/Interscience,1991.
    [21]秦安,冯前进,陈武凡.医学图像感兴趣区的内嵌小波编码新方法[J].计算机应用研究,2008,25(2):460-463.
    [22]张海燕,王东木,宋克欧等.系统仿真学报[J].2002,14(7):831-835.
    [23]肖振国,田心.医学图像无损压缩与有损压缩技术的进展[J].国外医学生物医学工程手册,2002,25(3):55-58.
    [24]Mark NelSon.数据压缩技术原理与范例[M].北京:科学出版社,1991.
    [25]张益贞,刘滔.Visual C++实现MPEG/JPEG编解码技术[M].北京:人民邮电出版社,2002.
    [26]钟建国,吴冬梅.基于二值自适应算术编码的图像无损压缩算法研究[J].陕西师范大学学报(自然科学版),2007,35(1):30-33.
    [27]原岛博.图像信息压缩[M].薛培鼎,徐国鼐译.北京:科学出版社,2004.
    [28]DvadiSalomon著.数据压缩原理与应用(第二版)[M].吴乐南等译.北京:电子工业出版社,2003.
    [29]罗军辉,冯平等.Matlab 7.0在图像处理中的应用[M].北京:机械工业出版 社,2005.
    [30]张德丰,马莉,范灵,梁忠宏.基于小波图像压缩技术的算法研究[J].中山大学学报(自然科学版),2008,47(2):42-45.
    [31]张克功.矢量量化和压缩编码技术在语音信号中的应用[J].兰州石化职业技术学院学报,2008,8(2):28-30.
    [32]朱元忠.图像压缩技术的发展研究[J].北京工业职业技术学院学报,2003,2(4):10-12.
    [33]Rabbani M Jones P.Digital Image Compression Techniques[J].SPIE,1991,1450:116-128.
    [34]肖自美.图像信息理论与压缩编码技术[M].广州:中由大学出版社,2000.
    [35]丁广贵,计文平,郭宝龙.Visual C++6.0数字图像编码[M].北京:机械工业出版社,2004.
    [36]王汇源.数字图像通信原理与技术[M].北京:国防工业出版社,2000.
    [37]周成兵.医学图像数据压缩方法研究[D].镇江:江苏大学,2005,4.
    [38]Jacobson R E,Ford AM,Arttidge G G.Evaluation of the effects of compression on the quality of images on a soft display[C].In:Proc SPIE,1997,3016:114-125.
    [39]王纪凤,王昌元,王世刚等.小波变换在医学图像压缩中的应用[J].中国医学物理学杂志,2008,25(4):741-743.
    [40]买买提明·艾海提,吴锡生.小波变换在图像压缩中的应用和研究现状[J].新疆教育学院学报,2008,24(2):128-130.
    [41]M.Boliek,C.Christopoulos and E.Majani.JPEG 2000 Part Final Draft International Standard[S].ISO/IEC JEC1/SC29/WG1N1855,2000.
    [42]李振伟,何继善,刘兵全.小波图像压缩进展[J].微机发展,2004,14(6):29-32.
    [43]David S.Taubman,Michael W.Marcellion.JPEG 2000 Image Compression Fundamentals Standards and Practice[M].Beijing:Publishnig House of Eelectronics Industry,2004.
    [44]赵志英,王景中,于德海.基于整数小波变换和SPIHT的ROI医学图像压 缩[J].北方工业大学学报,2005,17(3):10-14.
    [45]陈佩,张卫东.用lifting方法构造具有线性相位的双正交小波[J].电子与信息学报,2002,24(4):486-491.
    [46]A.R.Calderbank,Ingrid Daubeehies,Wim Sweldens Boon-lock yeo.Wavelet Transforms That Map Integers to Integers[J].Applied and Computational Harmonic Analysis,1998:332-369.
    [47]Adams,M.D.,Kossenili.F.Reversible integer-to-integer wavelet transform for image compression:Performance evaluation and analysis[J].IEEE Trans on Image process,2000,9(6):1010-1024.
    [48]Sweldnes W.The lifting scheme:A custonr design construction of biorthogonal wavelets[J].Appl Comput Harmon Anal,1996,3(2):186-200.
    [49]钱剑.基于小波变换的图像编码技术的研究[D].南京:东南大学,2004,5.
    [50]王良燕.基于小波变换的图像压缩算法研究[D].合肥:合肥工业大学,2007,5.
    [51]周先波,陈文胜.关于提升系数的确定与滤波的构造[J].中山大学学报(自然科学版),2001,40(4):35-37.
    [52]曹丽丹.基于SPIHT的无表零树压缩算法[D].吉林:吉林大学,2005,6.
    [53]董文超,丁振宇.一种基于整数小波变换的图像压缩方法[J].装甲兵工程学院学报,2007,21(2):66-69.
    [54]姜会亮.基于整数小波变换的静态图像压缩技术的研究及实现[D].扬州:扬州大学,2005,4.
    [55]吴宝明,候文生,彭承琳.基于小波变换统计特征的图像压缩算法的研究[J].生物医学工程学杂志,2002,19(3):513-517.
    [56]M.Antonini etal.Image coding using wavelet transformation[J].IEEE Trans on Image Processing,1992,1(2):13-17.
    [57]靳济芳.Visual C++小波变换技术与工程实践[M].北京:人民邮电出版社,2004.
    [58]刘伟章.基于小波变换的图像压缩编码研究[A].西安:西北工业大学.2005,12.
    [59]陈军,李云松,吴成柯.易于硬件实现的无链表小波零树编码算法[J].西安电子科技大学学报,2001,28(2):229-233.
    [60]陈淑珍,刘志军,宋维维等.改进的静止图像小波零树编码[J].武汉大学学报(理学版),2002,48(1):32-37.
    [61]李斌,庄天戈.实现图像感兴趣区域渐进编码的新方法-一种改进的嵌入式零树小波编码算法[J].电视技术,2002,(6):9-12.
    [62]曹怀信,赵建伟.小波分析发展综述[J].咸阳师范学院学报,2002,17(6):5-8.
    [63]王晓芳,杨荣荣.整数小波变换和EZW编码在基于ROI图像渐进传输中应用研究[J].计算机应用软件,2007,24(1):137-139.
    [64]贾志科,崔慧娟,唐昆.改进的SPIHT图像压缩编码方法[J].清华大学学报(自然科学版),2001,41(7):25-28.
    [65]卓力,沈兰荪.基于感兴趣区域的图像编码与渐进传输[J].电子学报,2004,32(3):411-415.
    [66]徐向民,刘伟,邢晓芬等.基于改进SPIHT算法的医学图像ROI编码[J].计算机工程与应用,2006,8:207-209.
    [67]赵志英.基于DICOM标准的医学图像压缩的研究[D].北京:北方工业大学,2005,6.
    [68]李雨,李学明.自适应算术编码及其在图像编码中的应用[J].计算机工程与应用,2003,38(32):73-75.
    [69]Christopoulos C,Askelop J,Larsson M.Efficient method for encoding regions of interest in the upcoming JPEG 2000 still image coding standard[J].IEEE Signal Processing Letters,2000,7(9):247-249.
    [70]胡荣锋,王玲等.基于JPEG 2000标准的ROI技术分析与实现[J].通信技术,2006,7:151-156.
    [71]王飞,郑南宁等.一种适合JPEG 2000编码的动态感兴趣区域提取算法[J].西安交通大学学报,2006,40(2):125-128.
    [72]蔡文君,李玲远,彭义波.基于JPEG 2000标准的感兴趣区域编码[J].电子技术应用,2005,8:44-49.
    [73]李安平,基于活动轮廓模型的轮廓提取[D].昆明:昆明理工大学,2006,5.
    [74]Xu chen-yang,Prince J L.Snakes,shapes and gradient vector flow[J].IEEE Transactions on image processing,1998,7(3):359-369.
    [75]Lu A D,Tang L,Xu Y H,et al.Segmetation of ultrasound images with interactive B-spline snakes and its application[J].Journal of software,2005,12(12):1760-1768.
    [76]Cohen MD,Davey MS.Imaging of gastrointestinal malignancy in childhood[J].Radiol Clin North Am,1996,34(4):717-742.
    [77]Curwen R,Blake A.Dynamic contours:real-time active splines[J].Blake A,Yullie Aeds Active Vision,1992:39-58.
    [78]陶玲,钱志余,陈春晓.主动轮廓模型及其在图像体分割中的应用[J].华南理工大学学报(自然科学版),2008,36(1):67-75.
    [79]马建林,崔志明,张娜敏,吴建,孙晓平.基于一种新的Snake模型的序列图像边缘提取算法[J].计算机应用,2008,28(4):982-1005.
    [80]曹松.JPEG 2000中的感兴趣区域编码技术浅析[J].计算机工程与应用,2003,39(24):67-78.
    [81]张宁.JPEG 2000中的感兴趣区域编码算法[J].计算机与数字工程,2003,31(4):24-27.
    [82]王柯俨.基于JPEG 2000标准的感兴趣区域编码方法[J].航天返回与遥感,2002.2:57-62.
    [83]Keun-hyeong Park,Chul Soo Lee,Hyun Wook Park.A multitiresolutional coding method based on SPIHT[J].Signal Processing:Image Communication,2002,(17):467-476.
    [84]Athanassions Skodras,Charilaos Christopoulos,and Touradj Ebrahimi.The JPEG 2000 still image compression Standard[J].IEEE Signal Processing Magazine,2001,9:325-328.
    [85]朱善安,刘丽芳.支持ROI优先编码策略的自适应码率控制算法[J].光电工程,2006,33(4):105-110.
    [86]Christopoulos C,Skodard A,Ebrehimi T.The JPEG 2000 still image coding system:an overview[J].IEEE Transactions on Consumer Electronics,2000, 46(4):1103-1127.
    [87]Joel Askelof,Mnatias Larsson Carlander,Charilaos Christopoulos.Region of interest coding in JPEG 2000[J].Signal Processing:Image Communication,2002,17:105-111.
    [88]陆波,魏询楷,毕笃彦.JPEG 2000中的感兴趣区域(ROI)编码技术[J].中国有线电视,2004,3:62-64.
    [89]徐磊,胡志忠.基于JPEG 2000的医学图像ROI压缩[J].医疗设备信息,2006,21(2):10-11.
    [90]ISO.ISO / ISC JTC 1/SC 29 //WG 1(ITU-TSG8),JPEG 2000[S].Part Ⅱ Final Committee Draft Version 110,2000(12):135-140.
    [91]朱磊.JPEG 2000中感兴趣区域的编码方法[J].科技信息,2007,(11):296-297.
    [92]蔡文君,李玲远,彭义波.基于JPEG 2000标准的感兴趣区域编码[J].电子技术与应用,2005,(8):44-49.
    [93]ISO.ISO / ISC JTC 1/SC 29 //WG 1(ITU-TSG8),JPEG 2000[S].Part Ⅰ Final Committee Draft Version 110,2000(3):148-151.
    [94]Liu L,Fan G.A new JPEG 2000 Region of interest image coding method:Partial Significant Bitpanes Shift[J].IEEE Signal Processing Letters,2003,10(2):35-38.
    [95]梁燕,刘文耀,郑伟.实现感兴趣区域编码的通用部分位平面偏移法[J].光电子·激光,2004,15(11):1334-1338.
    [96]R.Grosbois,D.S.Cruz,T.Ebrahimi.New approach to JPEG 2000 compliant region-of-interest coding[C].Proceedings of the SPIE 46th Annual Meeting on Applications of Digital Image Processing,2001:267-275.
    [97]Zhou Wang,Alan C.Bitlane-by-Bitplane Shift(BbBShift)-A suggestion for JPEG 2000 Region of Interest Image Coding[J].IEEE Signal Processing Letters,2002,9(5):160-162.
    [98]Han Xian Zhong,Zhang Li Bao.Region-Based Embedded Image Coding using Partial Subbands Maximum BitPlanes Shift[C].Proceedings of the Fifth International Conference on Machine Learning and Cybernetics, 2006:13-16.
    [99] Zhou Wang, Serene Banerjee, Brian L, et al. Generalized Bitplane-by-Bitplane Shift Method for JPEG 2000 ROI Coding[J]. IEEE International Conference on Image Processing, 2002, 8(3): 81-84.
    [100] Lijie Liu, Guoliang Fan. A new JPEG 2000 Region-of-Interest Image Coding Method:Partial Significant Bitplanes Shift[J]. IEEE Signal Processing Letter,2003, 10(2): 35-38.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700