劉 杰,陳懷新

(中國西南電子技術(shù)研究所,成都 610036)

1 引 言

圖像能夠提供給人們很多直觀的信息,比較符合人類的認(rèn)識規(guī)律,因而在多媒體通信、計算機系統(tǒng)、醫(yī)學(xué)、物理等眾多領(lǐng)域得到了越來越高的重視。但圖像一般具有龐大的數(shù)據(jù)量,對于通信傳輸網(wǎng)絡(luò)來說,受限的帶寬成為了阻礙其傳輸速度進一步加快的瓶頸;為了解決該問題,人們提出了許多有效的壓縮處理方法,如 JPEG[1]、JPEG2000[2]、壓縮感知[3]等。通過壓縮減少了通信網(wǎng)絡(luò)傳輸?shù)臄?shù)據(jù)量,但這些方法在對分布式壓縮的圖像數(shù)據(jù)進行重建的時候均采用子圖像的分別還原再拼接成完整圖像的技術(shù),在解碼端具有計算量和存儲量大的缺點,不能明顯地提高通信傳輸網(wǎng)絡(luò)的傳輸效率。

本文在基于小波變換的JPEG2000壓縮體制基礎(chǔ)上,提出了基于系數(shù)整合法的分布式圖像壓縮處理方法,本方法通過對原圖像進行分塊小波變換獲得子圖像小波系數(shù)向量,通信傳輸網(wǎng)絡(luò)只分別傳輸圖像壓縮的子圖像小波系數(shù)向量,從而大大降低了壓縮編碼系統(tǒng)的數(shù)據(jù)傳輸量;再對圖像壓縮解碼的子圖像小波系數(shù)進行關(guān)聯(lián)處理,形成原圖像的整體小波變換系數(shù)向量,并通過小波反變換進行原圖像的重構(gòu)。相較于傳統(tǒng)的處理方法,本方法無需直接進行子圖像還原及其圖像的拼接處理,因而能有效緩解解碼端較大數(shù)據(jù)量存儲要求和復(fù)雜的圖像配準(zhǔn)與拼接計算。最后,本文采用MATLAB仿真,對結(jié)果進行了詳細(xì)的數(shù)值分析。

2 基于小波變換的圖像壓縮方法

為了突破通信傳輸網(wǎng)絡(luò)帶寬受限的瓶頸,發(fā)展了許許多多的圖像數(shù)據(jù)壓縮方法,現(xiàn)在運用的最多也是最成熟的技術(shù)就是JPEG2000圖像壓縮標(biāo)準(zhǔn)。

2.1 離散小波變換

離散小波變換[4]在圖像壓縮處理中得到應(yīng)用,采用圖像小波變換系數(shù)的編碼壓縮體制已取代離散余弦變化成為圖像壓縮的新標(biāo)準(zhǔn)。

以此為基函數(shù),所構(gòu)造的小波變換就是離散小波變換:

其逆變換如下:

2.2 JPEG2000圖像壓縮方法

于2000年問世的靜止圖像壓縮標(biāo)準(zhǔn)JPEG2000在編碼端采用n×n的分塊圖像離散小波變換獲得不同分辨率頻段系數(shù),然后進行量化編碼得到了圖像壓縮數(shù)據(jù);在解碼端進行符號解碼后,對分塊壓縮數(shù)據(jù)進行小波反變換得到 n×n的子圖像,再將這些子圖像進行簡單拼接得到完整的大圖像[5],其簡要框圖如圖1所示。

圖1 JPEG2000壓縮框圖Fig.1 The compressed diagram of JPEG2000

圖像經(jīng)過小波變換后被分割成水平、垂直、對角線和低頻4個頻帶,低頻部分還可以繼續(xù)分解。水平、垂直和對角線部分表征了原圖像在3個方向的邊緣信息,具有明顯的方向特性,低頻部分代表著圖像整體的亮度信息。根據(jù)人類視覺的生理和心理特點,對4個部分分別作不同策略的量化和編碼處理,就可以實現(xiàn)圖像的壓縮[6]。

3 基于系數(shù)關(guān)聯(lián)整合的分布式圖像壓縮法

不同于傳統(tǒng)的方法,基于系數(shù)整合(CC)法的分布式圖像壓縮過程在編碼端先將圖像分塊,再進行子圖像小波變換及子圖像小波系數(shù)的量化編碼,解碼端進行小波系數(shù)的符號解碼與關(guān)聯(lián)整合,最后通過小波反變換實現(xiàn)原始圖像重構(gòu),其簡要框圖如圖2所示。

圖2 基于CC法的分布式圖像壓縮框圖Fig.2 The compressed diagram of distributed image based on coefficient conformity

在通信網(wǎng)絡(luò)的解碼端,如果進行子圖像的分別還原再拼接成原圖像反饋給計算機就需要進行復(fù)雜的計算和占用大量的存儲空間,在這種情況下,對系數(shù)進行整合就顯示出了其優(yōu)點,不需要存儲大量的圖像數(shù)據(jù),只需將能重點反映圖像信息的系數(shù)進行整合,最后將這些系數(shù)反饋給計算機終端,將較大的計算量交給功能強大的計算機去完成,這也達(dá)到了資源合理利用的目的。

本文提出的分布式圖像壓縮系數(shù)關(guān)聯(lián)整合處理如下。

首先,將一副N×N的原始圖像f(u,v)進行分塊,子圖像塊大小為n×n,共分成m個子圖像,這m個子圖的灰度值可表示為

其次,每個子圖像進行離散小波變換(CC法基于haar小波)得到4個頻帶信息。

低頻方向為

水平方向為

垂直方向為

對角線方向為

式中 ,l、h=1,2,…,n/4。

量化和編碼以及反量化編碼技術(shù)已經(jīng)相當(dāng)成熟,可以用霍夫曼編解碼[7]、循環(huán)編解碼[8]等技術(shù)進行處理,由于這不是本文研究的內(nèi)容,在這里就不進行詳細(xì)論述了。

考慮到傳輸過程中存在著噪聲影響,所以將接收到的各頻帶信息分別表示為

本文定義 Θq(l,h)為矩陣位置信息,定義Γi(l,h) Θq(l,h)為將 Γi(l,h)的值賦給合成的系數(shù)矩陣的第q塊的對應(yīng)位置。

再次,進行分塊壓縮編碼系數(shù)的關(guān)聯(lián)與組合處理。下面分無重疊和有重疊分塊的情況來分別討論系數(shù)關(guān)聯(lián)整合。

(1)無重疊分塊的情況

低頻方向的系數(shù)關(guān)聯(lián)整合如下:

水平方向的系數(shù)關(guān)聯(lián)整合如下:

垂直方向的系數(shù)關(guān)聯(lián)整合如下:

對角線方向的系數(shù)關(guān)聯(lián)整合如下:

式中 ,r、t=1,2,…,N/4。

(2)有重疊分塊的情況

假設(shè)在分塊的時候相鄰塊均重疊,由于m個子圖像運用不同的處理器進行處理,在進行系數(shù)關(guān)聯(lián)整合的時候重疊部分將進行加權(quán)整合,因為在實際應(yīng)用中各處理器的測量精度可能不同,則對其賦予不同的權(quán)值,將精度高的處理器賦予高的權(quán)值,反之,對精度較低的處理器則賦予較低的權(quán)值。令aj表示第j個處理器的權(quán)值,并且滿足0≤aj≤1,那么可得各個方向的關(guān)聯(lián)整合系數(shù)如下。

低頻方向的關(guān)聯(lián)整合系數(shù):

式中,第一項為無重疊部分的整合,第二、三項分別為兩塊重疊、四塊重疊部分的加權(quán)整合,且滿足 au+av=1和aw+aw+1+ak+ak+1=1。同理可得其它3個方向的關(guān)聯(lián)整合系數(shù)。

在整合分塊的小波系數(shù)基礎(chǔ)上,通過小波反變換進行原圖像的重構(gòu)。

4 仿真及結(jié)果分析

本文以256×256像素的Lena圖像作為仿真的對象,圖像大小為65 kbyte,分別將圖像分成無重疊和有重疊的4塊,無重疊分塊的每塊為128×128像素、16 kbyte,有重疊分塊的每塊為136×136像素、18 kbyte。假設(shè)數(shù)據(jù)在傳輸?shù)倪^程中受到服從標(biāo)準(zhǔn)正態(tài)分布的高斯白噪聲的影響,并根據(jù)本文第2.2節(jié)的方法和第3節(jié)提出的新方法,在MATLAB 7.8.0環(huán)境下進行圖像壓縮的數(shù)值仿真分析。

(1)仿真實驗1:無重疊分塊壓縮的仿真

對圖3所示無重疊分塊壓縮的還原效果如圖4所示。

圖3 無重疊的分塊Fig.3 The division without overlap

圖4 無重疊分塊壓縮的還原效果Fig.4 The reconstructed effect of the division without overlap

為了直觀地反映數(shù)據(jù)處理時效性,本文比較JPEG2000技術(shù)和基于CC法的分布式圖像壓縮處理方法(本文方法)的計算量,并采用圖像數(shù)據(jù)解壓重建時間(TOR)來衡量計算量大小。圖像的還原質(zhì)量通?？衫梅逯敌旁氡?、均方信噪比、均方根誤差以及平均誤差和(Mean Sum of Error,MSE)等來表示,本文采用原始圖像和解壓重建圖像對應(yīng)像素之間灰度值的MSE來評估圖像還原質(zhì)量,其計算公式如下:

式中,t為仿真次數(shù),gi(x,y)為第 i次仿真還原圖像的灰度值。

為了得到準(zhǔn)確、可靠的分析結(jié)果,需要通過大量的重復(fù)實驗來獲得實驗數(shù)據(jù),本文利用5000次的仿真來得到所需的仿真驗證結(jié)果。

基于CC法的無重疊分塊的分布式圖像壓縮和JPEG2000的圖像還原效果對比如表1所示。

表1 無重疊分塊壓縮的還原效果Table 1 The reconstructed effect of the division without overlap

(2)仿真實驗2:有重疊分塊壓縮的仿真

基于CC法的有重疊分塊(圖5)的分布式圖像壓縮和JPEG2000的圖像還原效果對比如表2和圖6所示。

表2 有重疊分塊壓縮的還原效果Table 2 The reconstructed effect of the division with overlap

圖5 有重疊的分塊Fig.5 The division with overlap

圖6 有重疊分塊壓縮的還原效果Fig.6 The reconstructed effect of the division with overlap

從圖4和圖6可以看出,基于CC法的有重疊和無重疊的分布式圖像壓縮還原的圖像效果很好,不論是JPEG2000還是基于CC法的分布式圖像壓縮處理方法,還原出來的圖像和原始圖像之間看不出任何明顯的差異,這說明了基于CC法的分布式圖像壓縮處理通過系數(shù)關(guān)聯(lián)整合的方法是有效和可行的。

對處理時效性利用解壓重建時間進行仿真對比分析,結(jié)果表明基于CC法的有重疊和無重疊的分布式圖像壓縮處理的解壓重建耗時分別為JPEG2000的1/8.8和1/17.6,若解碼端的計算處理能力遠(yuǎn)低于功能強大的計算機,處理時間將按比例的擴大,基于CC法的分布式圖像壓縮處理方法的優(yōu)勢在實際運用中將會體現(xiàn)得更加明顯。

從平均誤差和的結(jié)果可以看出,通過大量的重復(fù)仿真,基于CC法的無重疊分塊的分布式壓縮處理方法的誤差與 JPEG2000相當(dāng),其值略小于JPEG2000;基于CC法的有重疊分塊的分布式壓縮方法的誤差值比JPEG2000少了22.7%,這說明,兩種基于CC法的分布式壓縮方法都比JPEG2000的抗高斯白噪聲影響能力強,基于CC法的有重疊分塊的分布式壓縮方法效果更加明顯。有重疊的分塊方法通過犧牲少量的存儲和傳輸代價換來了較高的抗噪聲性能,對重疊部分的系數(shù)加權(quán)處理有效地減少了噪聲對傳輸信號的影響。

實際應(yīng)用中,在不明顯提高存儲量和計算壓力的前提下,可以通過重疊分塊的方法來達(dá)到較好的圖像還原效果。

5 結(jié) 論

本文在JPEG2000的基礎(chǔ)上,提出了基于CC法的分布式圖像壓縮處理方法。本方法將較大數(shù)據(jù)量的圖像數(shù)據(jù)分解為多個數(shù)據(jù)量較小的子圖像,再對子圖像進行小波分解,可以大大減少單個編碼系統(tǒng)的壓縮編碼數(shù)據(jù);采用分塊小波系數(shù)關(guān)聯(lián)整合處理,避免了傳統(tǒng)方法直接進行子圖像還原及其圖片配準(zhǔn)拼接處理,有效緩解了解碼端的存儲和計算壓力,并且具有很好的圖像重建效果。通過仿真驗證了本算法的有效性,其抗噪聲性能也優(yōu)于傳統(tǒng)的JPEG2000圖像壓縮結(jié)果。這為大數(shù)據(jù)量圖像信息在帶限通信網(wǎng)絡(luò)中高效傳輸提供了一種可行的方案,具有較高的實際應(yīng)用價值。

[1] 羅倩玲.JPEG編碼與矢量量化編碼方法比較[J].電子測量技術(shù),2005(4):35-37.LUO Qian-ling.The compare of coding between JPEG and vector quantization method[J].Electronic Measurement Technique,2005(4):35-37.(in Chinese)

[2] Varma K,Bell A.JPEG2000-choices and trade offs for encod-ers[J].Signal Processing Magazine,2004,21(6):70-75.

[3] David L Donoho.Compressed sensing[J].IEEE Transactions on Information Theory,2006,52(4):1289-1306.

[4] 張賢達(dá).現(xiàn)代信號處理[M].北京:清華大學(xué)出版社,2002:384-386.ZHANG Xian-da.Advanced Signal Processing[M].Beijing:Tsinghua University Press,2002:384-386.(in Chinese)

[5] 李定川.JPEG2000圖像壓縮標(biāo)準(zhǔn)及其關(guān)鍵算法的應(yīng)用[J].影像技術(shù),2010(4):26-31.LI Ding-chuan.The imagecom pressed standard of JPEG2000 and appliance of its pivotal algorithm[J].Video Technique,2010(4):26-31.(in Chinese)

[6] Woods J W,Naveen J.A Filter Based Bit Allocation Scheme for Sub band Compression of HDTV[J].IEEE Transactions on Image Processing,1992,1(7):436-440.

[7] ISO/IEC 10918-1,ITU-T Rec.T.81(JPEG),Digital compression and coding of continuous-tone still images[S].

[8] Jovanovic V M,Budisin S C.On the Coding Gain of Linear Binary Block Codes[J].IEEE Transactions on Communications,1984,32(5):635-638.