魏玉芬,野金花,丁艷清
(黑龍江八一農墾大學文理學院,黑龍江 大慶 163319)
嵌入式零樹小波編碼(Embedded Zerotree Wavelets Encoding,EZW)是基于比特連續逼近的圖像編碼方法:按位平面分層進行孤立系數和零樹的判決和熵編碼,而判決閾值則逐層折半遞減,又稱之為多層(或位平面)零樹編碼方法,這種編碼方法十分有效。近年來,在EZW算法的基礎上,出現了許多新的改進算法,如集合分裂嵌入塊編碼、可逆嵌入小波壓縮算法、最優截斷嵌入式塊編碼等,其中最優截斷嵌入式塊編碼算法,更是成為新圖像編碼標準JPEG2000中所采用的算法[1~3]。
EZW算法簡單,無需任何訓練,支持多碼率,具有較高的信噪比和較好的圖像復原品質。然而實驗結果表明,該算法由于不能完全利用小波系數的特點和充分考慮人眼視覺特性(HSV),導致編碼時間過長,編碼效率低,影響圖像復原品質。本文提出視覺加權小波零樹編碼方法,通過對高頻子帶進行視覺加權,修改高頻子帶系數掃描策略和量化閾值等方法進行編碼,對低頻子帶采用無損預測編碼方法,進一步提高圖像復原的品質和編碼效率。
EZW算法利用小波系數的特點,較好地實現了圖像編碼的嵌入功能,主要包括以下3個過程:零樹預測,掃描方法,逐次逼近量化[4]。
1.1 零樹預測
在變換編碼中,變換系數矩陣經過量化后,產生大量零符號。編碼的后續過程,就是有效地表示那些非零符號,包括非零符號的位置和大小。表示量化后非零值位置的過程,也就是表示有效值位置的過程,稱為有效值映射(significance map)。
為了構成一個完整的有效值映射,除了定義零樹根(ZTR)外,還需定義其他3個符號:
(1)孤零(IZ),表示當前系數值是無效值,但它的子孫系數中至少有一個是有效值;
(2)正有效值(POS),表示當前系數是一個正的有效值;
(3)負有效值(NEG),表示當前系數是一個負的有效值。小波樹結構示意圖如圖1所示。
1.2 掃描方法
EZW算法對小波系數進行編碼的次序叫做掃描。有了上述4個符號,可以按一定的順序掃描小波變換系數矩陣,從而形成一個符號表,它就是要得到的有效值映射。掃描過程中,各子帶按圖2所示的次序掃描。
圖1 三級小波分解和樹結構示意
圖2 子帶掃描順序
1.3逐次逼近量化
為了使零樹表示構成一個有效的嵌入式碼流,結合逐次逼近量化技術,來依次確定有效值和有效值映射,門限值之間滿足依次減半。
通過對人眼某些視覺現象的觀察,并結合視覺生理、心理學等方面的研究,發現恢復圖像的視覺失真,不僅取決于整體均方誤差,而且還取決于這個失真空間頻率上的分布。由于人眼具有視覺閾值效應和掩蓋效應,因此人的視覺分辨率力有限,可以容忍一定的圖像失真,但是人眼對不同頻率視覺信號的敏感度是不一樣的:
(1)人眼對圖像邊緣區信息的失真很敏感;
(2)人眼對圖像平滑區信息的失真比較感敏;
(3)人對圖像紋理區信息的失真不敏感。
長期以來,人們試圖建立人類視覺系統的數學模型,本文采用Ngan[5]提出的模型,即
其中,
a,b,c為常數;
f為視頻率,單位為周/像素(cycle per pixel,簡稱CPD)。
根據Ngan的研究,人眼對3~6 CPD的空間頻率分量最為敏感。利用上述視覺模型對小波系數進行加權,可把空間頻率f的單位變換為小波系數域(付氏域)的單位周/像素(cycle per pixel)。表1給出三級小波高頻系數的視覺加權值。
如圖3所示編碼流程圖,是將最低頻子圖和其他子圖分開處理的方法,用無損預測編碼對最低頻子圖單獨編碼,而對高頻圖像的編碼采用如下的改進方法:
(1)小波分解高頻系數進行視覺加權,對傳統的掃描順序進行改變,掃描不是從最后一級低頻系數開始,而是從最后一級高頻系數最大值開始,且按先水平方向高頻、垂直方向低頻,然后水平方向低頻、垂直方向高頻,再對角線方向高頻;
(2)選擇閾值的時候,是先在LHi最后一級分解的高頻系數中選擇,然后在LHi最后一級的高頻系數中選擇,最后是HHi中選擇;依然按照2的整數冪從高到低排列量化閾值。根據閾值選取的子帶順序不同,產生的零樹編碼碼流分別采用自適應算術編碼、Huffman編碼。
(3)如果第j層(三級小波變換,2≤j≤3)高頻子帶中的一個系數及其所有后代為零,定義該節點為j-1級ZTR(零樹根),然后依此類推。如果同一關系樹中有數個ZTR,則只對最高級ZTR進行編碼。解碼時,根據ZTR的級數,確定開始補零的層以及具體的起點;
(4)根據給定的閾值T做出一輪掃描后,將高頻子帶系數分成了3類:非零點、ZTR(零樹根))、IZ(孤立零)。顯然,要表現這3種情況都需要消耗額外的比特數,由于IZ的產生是因為第j層子帶的一個系數無意義,但其后代中有非零點。將IZ的量化值重新恢復進來,與非零點一并進行編碼,這時所有的高頻子帶系數就只有兩類了:非零點、ZTR(零樹根)。這樣在有效減小圖像失真的同時,提高了編碼效率。
(5)對除LL3子帶外的每一個高頻子帶,都建立一個編碼信息表,用來表明對應系數的分類情況。如果值為1,則表示對應系數為非零點,參與量化編碼;如果值為0,則表示對應系數為ZTR,不需要參與量化編碼。
圖3 HVEZW編碼流程圖
為了驗證上述算法的正確性,對MPEG-4標準測試序列圖像Claire第15幀灰度圖像進行EZW編碼壓縮和HVEZW編碼壓縮。如表1得出在不同的壓縮比下,得到的重構圖像的PSNR及MSE。
表1 EZW和改DFEZW對Claire圖像編碼結果
通過對比能夠得出下面的結論:在相同的壓縮比下,應用HVEZW算法,能夠提高重構圖像的PSNR、降低圖像的MSE,明顯地提高了圖像編碼的效率,即在較低碼率的情況下仍能很好的恢復圖像。這主要是因為HVEZW算法在進行零樹量化之前,對最低頻子帶進行單獨處理,對高頻子帶進行視覺加權,這樣在不影響圖像整體壓縮比的情況下,提高了編碼效率和恢復圖像的品質。如圖4,從視覺效果上來講,采用HVEZW算法重構圖像的視覺效果,要比采用EZW算法的重構圖像效果稍好。
圖4EZW編碼與DFEZW編碼結果比較
基于視覺特性的小波零樹圖像編碼壓縮方法,充分利用了小波系數的特點和人眼視覺特性,重建圖像的主觀效果比較好,達到了壓縮比高而失真度減少的目的。
[1]張旭東,盧國棟,馮 健.圖像編碼基礎和壓縮技術-原理、算法和標準[M].北京:清華大學出版社,2004.
[2]沃焱,韓國強,張 波.基于視覺特性的靜態圖象壓縮編碼[J].中國圖象圖形學報,2003,10(10):248-302.
[3]王 娜,賈傳熒.一種改進的無表零樹編碼的圖像壓縮方法[J].電子學報,2003,31(6):879-881.
[4]沈蘭蓀,卓 力,田 棟,汪孔橋.視頻編碼與低速率傳輸[M].北京:電子工業出版社,2001.
[5]Kunt M,Kocher M.Second generation image coding techniques[J].Proceedings of the IEEE,1985,73(4):549-575.