摘 要：針對(duì)醫(yī)學(xué)圖像信息量大的特征，提出基于視覺特性的多興趣區(qū)域圖像壓縮方案，通過仿真實(shí)驗(yàn)表明，該方法得到較高的壓縮比和壓縮視覺效果，為醫(yī)學(xué)圖像遠(yuǎn)程傳輸和存儲(chǔ)提供更為有效的方案。

關(guān)鍵詞：曲波；視覺冗余；感興趣區(qū)域；醫(yī)學(xué)圖像壓縮

1 概述

醫(yī)學(xué)圖像壓縮是醫(yī)學(xué)圖像研究領(lǐng)域中的一門重要課題[1]，近年來備受研究人員關(guān)注，特別在遠(yuǎn)程醫(yī)療，PACS（Picture Archiving and Communication Systems）醫(yī)學(xué)影像存檔與通訊系統(tǒng)中顯得至關(guān)重要。曲波（Curvelet）是一種新的多尺度分析方法，無論在分析二維圖像中的曲線（直線狀邊緣）特征或者是稀疏表達(dá)能力方面都比小波更勝一籌，文章充分利用曲波這一特性并結(jié)合醫(yī)學(xué)圖像特點(diǎn)，提出基于視覺特性的多感興趣區(qū)域圖像壓縮方案，即先根據(jù)視覺特性對(duì)醫(yī)學(xué)圖像提取人眼感興趣區(qū)域，再對(duì)不重要的區(qū)域進(jìn)行JPEG2000有損壓縮，消除不重要的視覺冗余，為醫(yī)學(xué)圖像遠(yuǎn)程傳輸和存儲(chǔ)提供更為有效的方案。

2 基于視覺冗余的多感興趣區(qū)域醫(yī)學(xué)圖像壓縮

2.1 基于視覺冗余的多感興趣區(qū)域壓縮具體過程

（1）根據(jù)需求人為選擇一個(gè)或多個(gè)感興趣區(qū)域。（2）醫(yī)學(xué)圖像數(shù)據(jù)分離。對(duì)圖像經(jīng)過人為選擇感興趣區(qū)域后，醫(yī)學(xué)圖像被分為R0I區(qū)域和非R0I區(qū)域，一般而言，R0I區(qū)域是一個(gè)連續(xù)的矩形區(qū)域，對(duì)這個(gè)連續(xù)的區(qū)域可直接送入JPEG2000編碼器，而對(duì)不連續(xù)的矩形區(qū)域非R0I圖像區(qū)域則是采用曲波迭代量化編碼。（3）圖像壓縮。分離后的圖像數(shù)據(jù)，將分別根據(jù)圖像數(shù)據(jù)系數(shù)矩陣送入Curvelet和JPEG2000編碼器最后得到各自分離部分的碼流。（4）碼流綜合。經(jīng)過編碼后的碼流包含了Curvelet碼流和JPEG2000碼流。文章利用曲波變換理論，根據(jù)人眼視覺冗余特性，人為劃分ROI和非ROI區(qū)域，在對(duì)非ROI系數(shù)進(jìn)行曲波迭代量化編碼，而對(duì)ROI區(qū)域則采用無損JPEG2000，從而達(dá)到在保證圖像重建后的人眼視覺效果的同時(shí)，也提高了圖像的壓縮效率。

2.2 第二代Curvelet變換理論

Donoho于1999年提出曲波變換理論[3]，即第一代Curvelet變換，第一代Curvelet的數(shù)字實(shí)現(xiàn)比較復(fù)雜，需要子帶分解平滑分塊、正規(guī)化和Ridgelet分析等一系列步驟，而且其中的Curvelet金字塔的分解帶來了巨大的數(shù)據(jù)冗余量，因此E.J.Candes等又提出了實(shí)現(xiàn)更簡(jiǎn)單、更便于理解的快速Curvelet變換算法，也稱第二代Curvelet變換。Curvelet變換是在研究小波變換的基礎(chǔ)上發(fā)展起來的，它克服了小波變換在應(yīng)用中的不足，顯示出了許多獨(dú)到之處。

2.3 根據(jù)視覺冗余特性，提取感興趣區(qū)域

醫(yī)學(xué)圖像感興趣區(qū)域的提取是圖像壓縮領(lǐng)域的熱點(diǎn)問題，文章采用人工干預(yù)的方法，利用人眼視覺冗余[4]的特性，即去掉那些人眼看不到或者可有可無的圖像數(shù)據(jù)。在圖像壓縮或量化過程中，人們根據(jù)視覺特性對(duì)圖像數(shù)據(jù)進(jìn)行取舍，如果取舍的圖像數(shù)據(jù)在接收端不能被視覺所感覺到，都被認(rèn)為圖像質(zhì)量是完好的或是足夠好的，那么這部分圖像數(shù)據(jù)就可以被舍棄。通過人工干預(yù)確定出ROI區(qū)域。選擇好ROI區(qū)域后，在程序中需要保留以下信息：即保留醫(yī)學(xué)圖像中ROI區(qū)域的相對(duì)位置信息，且將圖像平面分成N-1個(gè)感興趣區(qū)域ROI（1=1，…，N-l）和一個(gè)背景子區(qū)域，N-1個(gè)不同的感興趣區(qū)域需要采用不同的優(yōu)先級(jí)進(jìn)行編碼，一般情況下，感興趣區(qū)域是不規(guī)則形狀，一般為了計(jì)算方便，通常采用了矩形區(qū)域來表示ROI，對(duì)感興趣區(qū)域部分采用低壓縮率的無損壓縮編碼，而對(duì)非感興趣區(qū)域采用較高的有損壓縮，也降低了算法的復(fù)雜性。

3 實(shí)驗(yàn)結(jié)果與分析

文章采在MatLab環(huán)境中下實(shí)驗(yàn)上述算法，采用256×256像素的8bit深度的灰度醫(yī)學(xué)圖像子宮多發(fā)性平滑肌瘤CT為例，根據(jù)圖片的特點(diǎn)，方便觀看重構(gòu)效果，本案例劃分了單個(gè)的ROI，圖1（a）當(dāng)CR=40時(shí)，采用Huffman/EZW面向任務(wù)的醫(yī)學(xué)圖像編碼算法重建圖；圖1（b）當(dāng)CR=40時(shí)，采用文章所提出的結(jié)合Curvelet/JPEG2000編碼壓縮重建圖。

文章采用峰值信噪比（PSNR）和壓縮比（CR）兩項(xiàng)指標(biāo)對(duì)圖1進(jìn)行客觀評(píng)價(jià)，如圖2所示：

從圖2中可得出，由于圖像經(jīng)Curvelet變換后，存在大量冗余性，即使進(jìn)行迭代量化有效的降低冗余性，但在相同壓縮比的前提下，Curvelet/JPEG2000的PSNR是沒有JPEG和JPEG/JPEG2000高，明顯優(yōu)于EZW/Huffman算法，由于Curvelet對(duì)邊緣具有較強(qiáng)的表達(dá)的能力，在同壓縮比下，壓縮的視覺效果要明顯好于JPEG，很好的解決了JPEG的“蚊式效應(yīng)”，且隨著壓縮比的提高，壓縮效果將逼近JEPG2000/JPEG。JPEG/JPEG2000算法對(duì)感興趣區(qū)域采用JPEG2000的壓縮，因此壓縮效果要優(yōu)于JPEG，即圖像分別對(duì)ROI和非ROI區(qū)域進(jìn)行不同壓縮算法，其重建后質(zhì)量遠(yuǎn)遠(yuǎn)優(yōu)于對(duì)整副圖像采用單一的壓縮算法效果。

4 結(jié)束語

本章提出的基于Curvelet/JPEG2000相結(jié)合的多ROI醫(yī)學(xué)圖像壓縮方法，人為定義醫(yī)學(xué)圖像感興趣區(qū)域的基礎(chǔ)上，用戶可結(jié)合自己的需求與所需處理醫(yī)學(xué)圖像本身特點(diǎn)，結(jié)合曲波自身的變換特性和JPEG2000壓縮特點(diǎn)，選擇單個(gè)或多個(gè)感興趣區(qū)域，較好地解決了醫(yī)學(xué)圖像的高壓縮比和高質(zhì)量之間的矛盾，讓算法更具適應(yīng)性，不僅實(shí)現(xiàn)簡(jiǎn)單，運(yùn)行時(shí)間較短，在遠(yuǎn)程醫(yī)學(xué)圖像壓縮中具有很高的實(shí)踐性和應(yīng)用價(jià)值。

參考文獻(xiàn)

[1]康曉東.醫(yī)學(xué)影像圖像處理[M].人民衛(wèi)生出版社，2009.

[2]LIU Wei，YANGSheng，MediealImageROICompressionBasedon JPEG/JPEG2000 [J〕Measurement and Control Teehnology，2006.

[3]J.L.StarekE.J.CandesandD.L.Dohono.TheDenoising.IEEETrans.on

ImageProeessing[J].CurveletTransformforImage2002，11（6）：670-683.

[4]張帆.基于心理視覺冗余和PCA的圖像壓縮算法[J].科學(xué)技術(shù)與工程，2013，13（26）：7688-7691.

[5]胡敏，張長(zhǎng)江，魏晗，等.基于混合編碼和IWT的醫(yī)學(xué)圖像的近無損壓縮[J].光電工程，2008（05）：114-118（EI檢索）.

*通訊作者（文稿內(nèi)容責(zé)任人；在校生為其指導(dǎo)教師）：胡敏，女，講師。