摘 要：在分析了JPEG2000中的感興趣區(qū)域編碼算法的基礎(chǔ)上，針對(duì)其不足提出了一種基于區(qū)域重要性的嵌入式零樹小波編碼方法。算法無(wú)須修改小波系數(shù)值而用域值控制感興趣區(qū)的恢復(fù)質(zhì)量，并基于此修改了零樹編碼中的幅值編碼。為了使得符號(hào)流更適合游程編碼修改了零樹編碼中的掃描順序，獲得了較好的實(shí)驗(yàn)結(jié)果。

關(guān)鍵詞：感興趣區(qū);零樹編碼;幅值編碼;游程編碼

中圖分類號(hào)：TN919.8 文獻(xiàn)標(biāo)識(shí)碼：B 文章編號(hào)：1004373X(2008)1716503

Embedded Zero Tree Wavelet Coding Based on Area Significance

LIU Chunjuan

(Shaanxi Institute of Education，Xi′an，710061，China)

Abstract：Thorough analysis of the region of interest coding in JPEG2000，aiming at the detects，an improved embedding zerotree wavelet algorithm based on area significiance are proposed.This method uses threshold to control the reconstructed quality of the region of interest，and yet doesn′t need change the wavelet coefficient.On the basis of this idea，magnitude coding in EZW is modified.To make symbol stream adapt to runlenth coding，it uses a different scanning order with EZW.The experimental results are appreciated.

Keywords：interest；zero-tree code；amplitude coding；run-length coding

1 引 言

根據(jù)重構(gòu)圖像與原圖像是否完全一致，圖像壓縮可以分為有損壓縮和無(wú)損壓縮兩種。其中，有損壓縮可以獲得較高的壓縮比，但重建圖像的質(zhì)量相對(duì)較差，尤其在高壓縮比下，難以滿足需求；然而，無(wú)損壓縮圖像質(zhì)量雖然很好但由于其壓縮比很底，隨圖像數(shù)據(jù)量的增長(zhǎng)，其實(shí)用性也越來(lái)越小。實(shí)際應(yīng)用中，用戶往往只對(duì)圖像中某一特定區(qū)域的內(nèi)容感興趣，而且由人眼的視覺特性知觀察者對(duì)于圖像中的突變邊緣具有相當(dāng)?shù)囊曈X敏感度，這些區(qū)域統(tǒng)稱為感興趣區(qū)(ROI)，別的區(qū)域相應(yīng)稱為背景區(qū)(BG)?？梢酝ㄟ^向用戶傳送局部無(wú)損或近無(wú)損的高壓縮比圖像來(lái)解決以上的矛盾。另外也可通過漸進(jìn)編碼的方式，使呈現(xiàn)在用戶面前的圖像清晰度逐漸提高，用戶可以根據(jù)圖像質(zhì)量的狀況隨時(shí)切斷碼流。

小波變換的時(shí)頻局部化特性使上述想法成為可能。JPEG2000是一個(gè)基于小波變換的編碼系統(tǒng)，它主要采用的編碼方式為嵌入式零樹編碼(EZW)，而且其一大特色就是支持對(duì)感興趣區(qū)域的編碼，在高壓縮比下，感興趣區(qū)可以獲得比背景區(qū)高得多的圖像壓縮質(zhì)量。針對(duì)JPEG2000給出的基于感興趣區(qū)編碼方案的不足，筆者提出了一種新的基于區(qū)域重要性的零樹小波編碼方案，并在實(shí)驗(yàn)中取得了較好的結(jié)果。

2 JPEG2000中的感興趣區(qū)域編碼

2.1 JPEG2000中的感興趣區(qū)域編碼

JPEG2000支持的感興趣區(qū)域算法是基于對(duì)小波系數(shù)的放大scale。該算法的原則就是按比例增大與感興趣區(qū)域有關(guān)的小波系數(shù)(或按比例縮小背景區(qū)域的小波系數(shù))，使他們能夠位于較高的比特平面上。在接下來(lái)進(jìn)行的嵌入式編碼的編碼比特流里，這些被上移的系數(shù)被放在背景區(qū)的小波變換系數(shù)前。這樣感興趣區(qū)就會(huì)先于背景區(qū)進(jìn)行編碼和細(xì)化。即使編碼比特流被截?cái)嗔?，或者在?duì)整幅圖像完全編碼之前終止了編碼過程，感興趣區(qū)域也會(huì)擁有比圖像的其他區(qū)域高的逼真度。對(duì)于任意形狀的感興趣區(qū)域，如果使用普通的scale-based算法，其縮放值是任意的，因此需要編碼感興趣區(qū)形狀，以備解碼時(shí)將感興趣區(qū)系數(shù)按比例縮小，這樣增加了比特率。而JPEG2000采用的是MAXSHIFT算法，該算法是普通的scale-based的一種推廣和修訂。MAXSHIFT算法并不事先定義出期望的縮放值，而是首先由編碼器掃描量化系數(shù)，并自適應(yīng)的選出一個(gè)縮放值，使得感興趣區(qū)域系數(shù)里的最小值比所有的背景區(qū)系數(shù)都大。即選擇使s×min(ROI(i，j))≥max(BG(i，j))成立的最小整數(shù)s，記錄值s×min(ROI(i，j))為m。解碼器收到比特流后開始解碼每一個(gè)大于或等于m的系數(shù)一定是屬于感興趣區(qū)域的，應(yīng)該按比例縮小。但由于有限的設(shè)備精度且感興趣系數(shù)的最小值和背景區(qū)系數(shù)最大值差別往往很大，放大小波系數(shù)可能會(huì)導(dǎo)致溢出。

JPEG中，基于小波系數(shù)的按比例放大的方法的主要步驟如下：

計(jì)算小波變換;若已選定了感興趣區(qū)，計(jì)算感興趣掩模以確定感興趣系數(shù);按一定比例放大(縮小)感興趣區(qū)內(nèi)(背景區(qū))小波系數(shù);EZW編碼;熵編碼。算法比較如圖1所示。

2.2 小波零樹編碼概述

小波分解后信號(hào)的低頻部分對(duì)應(yīng)于原信號(hào)的平滑版本，而高頻部分對(duì)應(yīng)于這兩個(gè)信號(hào)的差別信息。一幅圖像經(jīng)過若干級(jí)小波分解后，在不同的子帶相同的空間位置的變換系數(shù)是相關(guān)的。在大的尺度上的小波系數(shù)稱為父節(jié)點(diǎn)，在較小的尺度上的相應(yīng)位置上的小波系數(shù)稱為子節(jié)點(diǎn)。對(duì)于小波3個(gè)高頻分解方向的子圖像的小波系數(shù)的所有的父節(jié)點(diǎn)都有4個(gè)子節(jié)點(diǎn)，對(duì)于最大尺度上的低頻子帶，每個(gè)父節(jié)點(diǎn)都有3個(gè)子節(jié)點(diǎn)。

由于圖像的小波系數(shù)經(jīng)過量化后，會(huì)出現(xiàn)大量的零系數(shù)，因此如何用最少的符號(hào)來(lái)表示這些零系數(shù)的位置，則是圖像提高壓縮比的關(guān)鍵。Shapiro提出的嵌入式零樹小波編碼提供了一個(gè)絕好的方法，小波零樹算法基于的重要假設(shè)是：首先，對(duì)于一個(gè)給定的域值T，一個(gè)系數(shù)x當(dāng)它的幅值不大于此門限時(shí)，即|x|

POS(正的重要系數(shù));NEG(負(fù)的重要系數(shù));IZ(孤立零);ZTR(零樹根)。一次掃描結(jié)束之后，輸出符號(hào)流并將POS和NEG的系數(shù)置為零，域值減半，繼續(xù)編碼輸出直到達(dá)到期望壓縮率或圖像質(zhì)量。

3 算法描述

通過以上對(duì)scale_based 及Maxshift的優(yōu)劣分析，本文給出一種新的基于區(qū)域重要性的小波圖像編碼方案。

3.1 本算法中的掃描方式

從小波分析的角度講，零樹代表圖像的平滑部分，當(dāng)圖像足夠平滑時(shí)就會(huì)產(chǎn)生零樹。零樹出現(xiàn)的位置取決于圖像的平滑程度和門限T的選擇。而且零樹編碼輸出是一系列符號(hào)流如果將這些符號(hào)流排序整理為連續(xù)出現(xiàn)的POS，NEG，IZ，ZTR，就可利用游程編碼極大的提高編碼效率。在前敘之零樹編碼中已經(jīng)知道，T給定時(shí)每一個(gè)小波系數(shù)屬于哪一種符號(hào)是確定的，但它何時(shí)出現(xiàn)是由掃描順序決定的，那么適當(dāng)?shù)匦薷脑搾呙桧樞蚓涂蓾M足上述愿望。

我們將低頻部分劃分為n*n個(gè)小塊，每塊具有相近的像素值表示塊內(nèi)平滑。n的大小視圖像的平滑程度不同而定，對(duì)于復(fù)雜度大的相關(guān)性較小的圖像，n選擇小些。本文采用的是4*4的塊。人為選擇感興趣區(qū)或?qū)γ恳粔K計(jì)算其方差那些大于低頻子圖像總方差的小塊為客觀感興趣區(qū)(這類小塊對(duì)應(yīng)圖像中邊緣部分)。由于我們將圖像粗略按平滑程度分塊，而且考慮到在二維小波變換中，低通濾波方向比高通濾波方向更能體現(xiàn)出圖像的趨勢(shì)，因此對(duì)每個(gè)HL，LH帶內(nèi)塊內(nèi)塊間均沿L方向進(jìn)行，如圖4所示。

圖4 本文算法中的掃描順序按這樣的順序掃描后，所得連續(xù)出現(xiàn)某類符號(hào)平均長(zhǎng)度根據(jù)所實(shí)驗(yàn)的圖像的平滑度不同而不同，當(dāng)用標(biāo)準(zhǔn)圖lena為實(shí)驗(yàn)圖左眼為感興趣區(qū)，碼率為0.125 b/pixel時(shí)，連續(xù)出現(xiàn)零樹根的平均長(zhǎng)度提高了9%，連續(xù)出現(xiàn)正重要系數(shù)平均長(zhǎng)度提高了7%，孤立零和負(fù)重要系數(shù)連續(xù)出現(xiàn)的平均長(zhǎng)度也均有10%和2%的提高。平均長(zhǎng)度按下式計(jì)算：ι^-=∑w_i×n_i/N(連續(xù)出現(xiàn)w_i個(gè)某類符號(hào)有n_i個(gè)，N是連續(xù)出現(xiàn)某類符號(hào)的總個(gè)數(shù))。但對(duì)標(biāo)準(zhǔn)圖barbara上述各平均長(zhǎng)度同原零樹編碼基本相同。

3.2 感興趣掩模的生成及雙域值零樹編碼

感興趣掩模的生成采用的是塊零樹的方法，即設(shè)某一大尺度上(低頻部分除外)感興趣塊的起始坐標(biāo)為(i，j)，結(jié)束坐標(biāo)為(s，t)，則其對(duì)應(yīng)的相同空間位置的小尺度上感興趣子塊的起始坐標(biāo)為(2×i－1，2×j－1)，結(jié)束坐標(biāo)為(2×s，2×t)，其實(shí)質(zhì)依然是四叉樹結(jié)構(gòu)(見圖3)。當(dāng)然這樣的構(gòu)造感興趣掩模的方法會(huì)造成小塊的邊界部分失真，但本文著眼的是塊內(nèi)圖像的高質(zhì)量恢復(fù)，而且JPEG所支持之感興趣掩模的確定依賴于所選小波基的不同而不同^［3］，具有較高的計(jì)算復(fù)雜度。

感興趣系數(shù)確定后本算法并不改變ROI或BG區(qū)小波系數(shù)的值，為在高碼率壓縮時(shí)，使ROI獲得比BG區(qū)質(zhì)量高的恢復(fù)圖像，在小波零樹編碼時(shí)對(duì)ROI或BG區(qū)使用不同的域值控制重要系數(shù)的輸出。BG區(qū)初始域值的選擇同零樹編碼即取T=2^logmax(N(i，j))，其中N(i，j)為(i，j)處的小波系數(shù)，ROI的初始域值RT小于T，它的值依賴于碼率或恢復(fù)圖像質(zhì)量的控制。實(shí)驗(yàn)中我們選擇RT=T/23。

EZW編碼中幅值編碼時(shí)需根據(jù)域值T的改變而相應(yīng)增加一個(gè)比特描述其精度，由于本文著眼于低比特的圖像壓縮，域值不會(huì)達(dá)到很低則圖像的小波系數(shù)在編碼后無(wú)須精確重構(gòu)。本文采用的BG幅值編碼是將區(qū)間［T，2T］分成以8為區(qū)間長(zhǎng)度的T/8個(gè)區(qū)間，判斷待編碼幅值屬于第幾個(gè)區(qū)間，用此區(qū)間的編號(hào)n編碼該幅值。解碼時(shí)用所對(duì)應(yīng)區(qū)間的中間值進(jìn)行重構(gòu)，即1/2×［T+(n－1)×8+T+n×8］，這樣每個(gè)幅值的重構(gòu)誤差不大于4。這樣的編碼方式保證了重要系數(shù)的重構(gòu)質(zhì)量。ROI重要系數(shù)幅值第一次編碼時(shí)先將［RT，2T)等分為以RT為區(qū)間長(zhǎng)度的(2T－RT)/RT=m個(gè)區(qū)間，判斷該值屬于哪個(gè)區(qū)間輸出區(qū)間代碼，然后再將所對(duì)應(yīng)區(qū)間如上所述BG區(qū)幅值編碼方法來(lái)編碼(見圖5)。這就需要增加log m個(gè)比特來(lái)描述，不過因?yàn)榫幋a比特中符號(hào)流編碼占絕大部分(即重要系數(shù)位置比特開銷要遠(yuǎn)大于重要系數(shù)值的開銷)而且此增加的比特僅在第一次掃描時(shí)增加，RT往后也如T一樣減半進(jìn)行下一次編碼掃描，無(wú)須增加比特，這樣所增加的比特就微不足道。

4 實(shí)驗(yàn)結(jié)果

以標(biāo)準(zhǔn)測(cè)試圖lena，Barbara(大小為256*256，每像素8 b)為例，給出了實(shí)驗(yàn)結(jié)果，在編碼之前，圖像經(jīng)過4級(jí)小波分解。對(duì)lena圖選擇左眼所處的方塊為感興趣區(qū)，對(duì)barbara圖選擇面部所在的方塊為感興趣區(qū)。

5 結(jié) 語(yǔ)

本文分析了JPEG2000中感興趣編碼的方法的優(yōu)劣，提出了一種基于視覺重要性的小波零樹編碼算法。該算法不對(duì)感興趣區(qū)域小波系數(shù)值進(jìn)行放大避免了溢出，而用不同的域值來(lái)控制感興趣區(qū)和背景區(qū)的重要系數(shù)的輸出，基于此修改了EZW中的幅值編碼方案，獲得了較好的實(shí)驗(yàn)結(jié)果。本文中感興趣區(qū)峰值信噪比的提高是以背景區(qū)的失真為代價(jià)的使得在低比特壓縮時(shí)感興趣區(qū)可以獲得比背景區(qū)高的恢復(fù)圖像質(zhì)量。而且在EZW編碼時(shí)運(yùn)用了不同的掃描順序，改善了對(duì)符號(hào)流的編碼使之更有利于游程編碼，有較好的應(yīng)用前景。

參 考 文 獻(xiàn)

［1］Shapiro J M.Embedded Image Coding Using Zerotrees of Wavelet Coefficients.IEEE Trans.on S.P.，1993，41(12):3 445-3 462.

［2］Christopoulos C，Skodras A，Ebrahimi T.The JPEG2000 Still Image Coding System:An Overview.IEEE Transactions on Consumer Electronics，2000，46(4):1 103-1 127.

［3］陳軍，吳成柯，李云松.基于零樹結(jié)構(gòu)的感興趣區(qū)圖像內(nèi)嵌編碼算法［J］.西安電子科技大學(xué)學(xué)報(bào)，2002，29(3):343-347.

［4］黃卓君，馬爭(zhēng)鳴.一種零樹與游程相結(jié)合的小波圖像編碼方法［J］.中國(guó)圖像圖形學(xué)報(bào)，2001，6(11):1 118-1 124.

［5］張孝杰，張專成，李研.一種改進(jìn)的嵌入式零樹小波圖像編碼算法.計(jì)算機(jī)應(yīng)用，2004，24(7):52-53.

［6］李津平，芮小平，楊崇俊.一種基于改進(jìn)EZW算法的圖像壓縮方法 ［J］.計(jì)算機(jī)應(yīng)用研究，2006，23(8):174-196.

［7］馬社祥，劉鐵根.變深度嵌入小波零樹自然圖像壓縮算法 ［J］.計(jì)算機(jī)工程，2006，32(20):1-3，61.

［8］陳偉，吳謹(jǐn).基于EZW的圖像編碼改進(jìn)算法 ［J］.武漢科技大學(xué)學(xué)報(bào):自然科學(xué)版，2006，29(2):142-144，147.

［9］盧輝斌，耿玉靜，趙華.一種基于小波數(shù)據(jù)形態(tài)學(xué)描述的圖像編碼算法 ［A］.計(jì)算機(jī)技術(shù)與應(yīng)用進(jìn)展——全國(guó)第17屆計(jì)算機(jī)科學(xué)與技術(shù)應(yīng)用(CACIS)學(xué)術(shù)會(huì)議論文集(上冊(cè)) ［C］.2006.

［10］王向陽(yáng)，楊紅穎.基于改進(jìn)零樹集合的新圖像編碼算法 ［J］.測(cè)繪學(xué)報(bào)，2003，32(1):58-62.