CCSDS星載圖像壓縮模塊的FPGA設計與實現(xiàn)-

1、收稿日期:2008-11-27.基金項目:國家自然科學基金項目(60805030.光電技術應用CCSDS 星載圖像壓縮模塊的FPGA 設計與實現(xiàn)張學全1,2,顧曉東1,孫輝先1(1.中國科學院空間科學與應用研究中心,北京100080;2.中國科學院研究生院,北京100080摘要:根據(jù)航天應用的需要,采用CCSDS122.02B 21推薦的新一代空間圖像壓縮算法,基于FP GA 設計實現(xiàn)了高速星載圖像壓縮模塊。針對算法的特點,設計了整體架構。圖像采用分塊編碼的方法,對小波變換,采用了邊界補償?shù)娜哂喾謮K算法,提高了信噪比,位平面編碼各模塊采用并行與流水結合的方法,提高了掃描和編碼速率。通過測試,系

2、統(tǒng)性能指標達到了航天任務需求。關鍵詞:CCSDS ;圖像壓縮;小波變換;位平面編碼;FP GA 中圖分類號:V447文獻標識碼:A 文章編號:1001-5868(200906-0935-05Design and Implementation of CCSDS 2based Onboard Image Compression U nit Using FPG AZHAN G Xue 2quan 1,2,GU Xiao 2dong 1,SUN Hui 2xian 1(1.C enter for Space Science and Applied R esearch ,Chinese Academy

3、of Sciences ,B eijing 100080,CHN;2.G radu ate School of Chinese Academy of Sciences ,B eijing 100190,CHN Abstract :Based on CCSDS122.02B 21algorit hm ,a high speed o nboard image compression module is designed and implemented on FP GA.The f unctional f ramework is p roposed after analyzing t he char

4、acteristics of t he algorit hm.The images are processed by blocking and ext ra pixels are compensated for wavelet t ransform ,from which t he compression performance is imp roved.Test result s verify t he applicability of t his designed module in space applications.K ey w ords :CCSDS ;image compress

5、ion ;wavelet ;bit plane encoding ;FP GA0引言星載圖像壓縮歷來是衛(wèi)星在軌數(shù)據(jù)處理的關鍵技術之一。隨著國際國內空間探索活動日益頻繁,星上圖像采集設備分辨率不斷提高,目前高分辨率星載遙感的原始數(shù)據(jù)率達到數(shù)個Gbit/s ,給數(shù)據(jù)緩存和傳輸帶來巨大的壓力。由于大容量存儲器擴展能力和數(shù)據(jù)下行帶寬有限,高速低比特率圖像壓縮就成為解決這一困難的有效途徑。以通用處理器為核心的壓縮方案通常需要較高的時鐘輸入(大于100M Hz 才能滿足高速需求,不符合空間環(huán)境對星上電磁兼容的約束;而圖像壓縮專用ASIC 又不能靈活滿足任務需要,無法進行功能裁剪,且不能保證供貨;采用FP

6、GA 實現(xiàn)壓縮算法雖然較為復雜,但能夠克服上述缺點,且保密性好,成為壓縮系統(tǒng)研制方案的合理選擇。CCSDS (Consultative Committee for Space Data Systems ,空間數(shù)據(jù)系統(tǒng)咨詢委員會2005年正式推出了新一代空間圖像壓縮推薦標準CCSDS122.02B 211,它針對空間圖像處理的特點,協(xié)調了算法性能與硬件實現(xiàn)復雜度的矛盾,兼容CCD 相機的幀格式輸出和推掃格式輸出,適用于從紅外,可見光到紫外等譜段的二維數(shù)字圖像壓縮。本文基于Xilinx FP GA ,采用CCSDS 的推薦算法,設計實現(xiàn)了空間圖像壓縮模塊。文章在第二部分對CCSDS 算法原理和小波

7、分塊算法做簡要介紹,第三部分闡述基于FP GA 的算法整體架構及各個模塊結構的設計和實現(xiàn)方法,第四部分介紹實驗結果,最后給出結論。5391壓縮算法分析CCSDS 算法采用了一種廣泛應用的變換+編碼模型,如圖1所示。輸入圖像數(shù)據(jù)首先利用離散小波變換(DW T ,Discrete Wavelet Transform 進行頻帶分解,以達到能量聚焦和去除相關性的目的,之后模塊對變換得到的小波系數(shù)進行位平面掃描和熵編碼,最后拼接成輸出碼流。由于采用了位平面編碼(B PE ,Bit Plane Encoding ,可對嵌入式碼流隨機截斷,使得圖像可以漸進傳輸,壓縮比/壓縮質量可以連續(xù)調節(jié)。 圖1CCSDS

8、 壓縮算法模型1.1離散小波變換DW T 能夠對數(shù)據(jù)高頻成分在由粗到細漸進的空頻域上正交分解,像顯微鏡一樣處理任意局部特征,實現(xiàn)高質量壓縮。CCSDS 推薦了9/7浮點和整數(shù)兩種雙正交小波變換方法,整數(shù)變換可采用定點運算,節(jié)省資源,運算量小,適合于FP GA 實現(xiàn),且能夠實現(xiàn)完全重構,故本文算法采用整數(shù)小波變換,一維提升公式如式(1、(21。由于二維變換是一維變換的張量積,等效于將圖像沿水平和垂直方向順次進行一維分解。分解得到LL 1、HL 1、L H 1和H H 1子帶,其中LL 1為能量聚集的低子帶,可繼續(xù)進行同樣的分解。根據(jù)協(xié)議規(guī)定,將原始圖像進行三級二維分解,得到10個子帶,如圖2所示

9、,子帶加權后,即可進行位平面編碼。 圖2小波變換系數(shù)位置分布圖D j =x 2j +1-916(x 2j +x 2j +2-116(x 2j -2+x 2j +4+12(1C j =x 2j -D j -1+D j2+12(21.2位平面編碼三級小波變換后的空頻圖像,如圖2,不同頻帶的系數(shù)形成了樹狀關系,每個DC 系數(shù)對應3個四叉樹,共64個成員組成一個家族。位平面編碼模塊首先根據(jù)位置關系,將同一家族成員聚集成8×8的數(shù)據(jù)塊,作為一個單元,并按照相應DC 在LL 3中的排列順序,將各單元依次排列,然后開始掃描和編碼。掃描前首先設定S 個連續(xù)排列的單元為一個段,每次掃描編碼只針對一個段

10、。段與段之間相互獨立掃描編碼,以防止空間信道傳輸中數(shù)據(jù)丟失的情況下,造成解碼端誤碼擴散。編碼分交、直流系數(shù)分別進行。在直流系數(shù)編碼中,DC 系數(shù),AC 深度先按照協(xié)議中給定的量化表進行量化,然后對量化系數(shù)進行Rice 編碼。交流編碼分位平面掃描和熵編碼兩部分。位平面掃描針對交流系數(shù)深度設定一系列閾值,按重要性由高到低,每次得到小波系數(shù)在一個位平面的有效值,對這些有效值按照1中的掃描算法分5階段順序進行,形成轉義字,類型值等中間符號,再對這些符號進行熵編碼,編碼采用查表的方法。1.3圖像分塊算法通常,空間遙感圖像的尺寸都很大(行列通常都大于500像素,如果對整幅像進行變換及編碼處理,從時間延遲和

11、資源占用上考慮都是不合理,甚至難于實現(xiàn)的。面對FP GA 有限的片內邏輯和存儲資源,為了實現(xiàn)高速數(shù)據(jù)壓縮,本文采取分了塊的方法,將整幅圖像分成小塊,對每一塊分別進行處理,塊間進行流水操作,在接收端可以分別解碼重構合成圖像。這種方法,壓縮核只對分塊圖像進行處理,與整幅圖像大小無關,大大提高了壓縮效率和壓縮模塊的適用性。但直接將圖像獨立分塊的方法,一般采用延拓來彌補邊界數(shù)據(jù)的不足,在有損壓縮的情況下,會丟失塊間像素的統(tǒng)計相關性,造成解碼重構圖像的塊效應。本文對9/7整數(shù)小波采用補償邊界的冗余分塊方法2,以較小的代價解決了邊界信息丟失的問題。具體方法如下:以長度為32的一維數(shù)據(jù)為例。由9/7小波變換

12、公式(1、(2可知,每進行一級小波變換,需要在639前一級數(shù)據(jù)邊界左側補充4個數(shù)據(jù),右側補充3個數(shù)據(jù),以此類推,可得三級變換所依賴的原始數(shù)據(jù)范圍如圖3,其中雙豎線為邊界線。由圖可見,最終得到4個直流分量,需要在左側補充28個數(shù)據(jù),在右側補充21個數(shù)據(jù),才能完成運算。因此在本文中,我們每次送入圖像塊時,將其三級變換所依賴范圍的像素數(shù)據(jù)全部送入,則變換所得到的結果就與整幅圖像小波變換所得結果完全一致,避免額外引入塊效應。采用此種方法的代價是增加了約5.4倍的片內存儲,和3.7倍計算量,但因為所分圖像塊較小,并不會對硬件資源造成過大壓力,相對于整幅圖像的計算仍具有較大優(yōu)勢 。圖3三級小波變換邊界補充

13、數(shù)據(jù)范圍2算法的FP GA 設計與實現(xiàn)采用自頂向下的方法,首先設計總體框架,再分別設計功能模塊。2.1算法實現(xiàn)的總體框架如圖4所示,由于壓縮速率和外部數(shù)據(jù)輸入速率存在差異,所以原始圖像數(shù)據(jù)首先輸入緩存,再將緩存內容切分成小塊圖像送入小波變換模塊,變換的中間結果送入小波系數(shù)緩存,當變換完成后,BPE 模塊接到控制指令,從小波系數(shù)緩存中讀取數(shù)據(jù),先進行家族數(shù)據(jù)集合,加權,掃描系數(shù)深度,直流系數(shù)量化等預處理,之后對DC 系數(shù),AC 深度進行Rice 編碼,對交流系數(shù)進行位平面掃描和編碼。直流和交流處理相互獨立,可并行執(zhí)行,以減少壓縮時間,所得到的結果送入碼字拼接模塊,形成連續(xù)的碼流送入輸出緩存。由于

14、DW T 模塊與B PE 模塊是兩個功能獨立的模塊,因此兩模塊間工作為流水模式,如圖5,使得壓縮周期由t DWT +t BPE 變?yōu)閙axt DWT ,t BPE ,因此可以減少時間延遲,提高壓縮效率 。圖4 算法實現(xiàn)整體框架圖5整體流水結構示意圖2.2DWT 模塊的電路結構DW T 模塊包含分塊和小波變換兩部分內容。分塊按照1.3所述方法,每次從緩存讀取67×67大小圖像數(shù)據(jù)塊送入小波變換模塊計算,相鄰讀取的圖像塊會有重疊。對小波變換來說,變換速率和資源占用情況直接影響到整體性能。關于小波變換的硬件實現(xiàn)方法,已有許多相關文獻324作了報道。本文中,因為壓縮速率的瓶頸取決于周期最長的

15、流水環(huán)節(jié),即B PE 環(huán)節(jié),所以從算法整體框架出發(fā),可以適當放寬對小波變換速率的要求,轉而面向資源節(jié)省考慮問題。因此,本文采取復用一維變換的方法來實現(xiàn)三級二維變換。結構如圖6。為降低存儲控制的復雜性和保持數(shù)據(jù)同時讀取和存儲,分配兩塊同等大小的片內緩存。輸入數(shù)據(jù)經(jīng)接口控制模塊進入一維小波變換,行變換后的結果送入緩存2,經(jīng)緩存2讀取數(shù)據(jù)繼續(xù)送入一維小波變換,列變換的結果送入緩存1,此時完成一級變換。接下來由緩存1讀取相應位置的數(shù)據(jù),重復上述過程進行第二級和第三級變換,最終變換結果存放在緩存1中,由接口控制模塊控制結果輸出 。圖6小波變換實現(xiàn)結構一維小波變換采取提升格式,將乘法轉換成加法和移位操作,

16、方法如下:由式(1、(2可以得到其等效表達式(3、(4:D j =x 2j +1-116(x 2j +x 2j +2-x 2j -2-x 2j +4+8(x 2j +x 2j +2+1(3739C j =x 2j -1-D j -1-D j2(4由式(3、(4可得到計算步驟,描述如下:step1:temp1<=x 2j +x 2j +2-x 2j -2-x 2j -4;temp2<=(x 2j +x 2j +2+1<<<3;step2:D j <=x 2j +1-(temp1+temp2>>>4;step3:C j <=x 2j -(

17、2-D j -1-D j >>>2;step 之間采用流水操作,保證每個周期輸出一個結果,以此保持輸出緊湊連續(xù)。2.3BPE 模塊的電路結構BPE 模塊針對一個段的數(shù)據(jù)進行處理,包括預 處理,DC ,ACDept h 編碼,塊內掃描,AC 熵編碼及碼字拼接五個部分。2.3.1預處理首先進行編碼前的預處理。當B PE 模塊從總體控制模塊得到DW T 模塊完成信號時,即開始從小波系數(shù)緩存中讀取數(shù)據(jù),根據(jù)讀取數(shù)據(jù)的地址,可以判斷小波系數(shù)的頻帶位置及家族歸屬情況,以此對小波系數(shù)進行加權和家族聚合,形成掃描單元,同時對讀取的小波系數(shù)測量數(shù)據(jù)深度,最終得到整個段內各掃描單元的DC 深度和

18、AC 深度(ACDept h ,加權表見1。段內全部數(shù)據(jù)讀取之后,暫存于寄存器中,然后對DC 進行量化,量化表見1。2.3.2直流編碼預處理之后,編碼分兩路并行進行。一路是針對DC 量化值和ACDept h 的Rice 編碼,另一路是針對DC 的量化余數(shù)和AC 系數(shù)的位平面編碼。在Rice 編碼中,包括差分,映射,編碼參數(shù)選擇,和變長編碼幾個步驟。其中,編碼參數(shù)的選擇有遍歷尋優(yōu)和次優(yōu)公式兩種方法,為節(jié)省運算,本文采用次優(yōu)公式計算編碼參數(shù)。由于Rice 編碼實現(xiàn)方法較為成熟,本文不作詳細闡述。2.3.3位平面掃描與編碼AC 系數(shù)和DC 量化余數(shù)的位平面編碼,是算法的核心部分,由塊內掃描和AC 熵

19、編碼兩個模塊完成。如圖7所示,bp 表示位平面,B 表示塊,N 為AC 深度,S 為段內所含掃描單元的數(shù)目。輸入數(shù)據(jù)首先存放在寄存器中,AC 系數(shù)按照位平面從高到低逐層掃描,一個位平面包含S 個掃描單元,在一個掃描單元中,系數(shù)掃描的次序如圖8所示,分為五個stage 。stage0為DC 余數(shù)掃描,stage1,stage2,stage3為掃描AC 系數(shù)并形成轉義符號,對上層位平面已掃描中的像素在本平面中由stage4直接記錄剩余數(shù)值而不做符號轉換。掃描的過程中,下一層的掃描結果與上一層的掃描結果緊密相關,因此要記錄掃面歷史信息,在生成新掃描字時進行參考。由于掃描單元間無交叉信息,所以在本文的

20、設計中,采用S 個單元并行掃描以提高編碼效率。掃描中位平面間存在順序關系,但各位平面編碼是相互獨立的,所以多個位平面間掃描與編碼可以采用流水結構處理。掃描模塊對一個位平面進行掃描,掃描過后即送入編碼模塊進行編碼。在編碼的過程中,首先要計算編碼選項,方法是遍歷所有選項的碼字長度,累加比較,挑選使碼長最小的選項,之后,即可根據(jù)選項查表映射碼字。映射表見1。2.3.4碼字拼接與隨機截斷由于碼流以8位或16位字節(jié)為輸出單位,必須將編碼輸出的零散變長碼字拼接成以字或字節(jié)為單位的碼流。本文實現(xiàn)結構如圖9。當讀信號有效,則讀取輸入碼字,截取有效位,與剩余位拼接,如果湊足了新的單位碼流則輸出,并將剩余位右移;

21、若不足以形成新的單位碼流,則拼接后一同作為剩余位暫存。同時,為了對壓縮比進行靈活控制,要求可以隨時截斷碼流,在碼字拼接并輸出的環(huán)節(jié)添加位計數(shù),一旦計數(shù)器滿足了壓縮比,將立刻輸出停止編碼信號。這樣可以省去不必要的計算。為了防止碼流不足的情況發(fā)生,當讀信號無效而補0信號有效時,將不輸入碼字,輸入8位0,用0補充不足的碼流。839 圖9碼字拼接及數(shù)據(jù)截斷功能結構3測試及結果分析本文采用Xilinx 的XC2V6000 器件作為設計驗證平臺,以此器件為核心搭建了驗證系統(tǒng)。上層為驗證板,FP GA 芯片外圍配有4片RAM ,下面有一層基于U SB 的接口板,兩個U SB 口與PC 機相連,分別為上行和下

22、行數(shù)據(jù)傳輸通道,由PC 機作為圖像源和接受方,可對不同大小的圖像以可變的發(fā)送速率進行測試本文選取一系列512×512的8bit 遙感灰度圖進行測試,采用40M Hz 時鐘,進行8倍壓縮,耗時均為35ms 左右,峰值信噪比均大于30dB ,資源占用率約為60%,能夠滿足某些航天應用的需求。圖10是其中三幅圖像的原圖和解壓圖對比,壓縮時間和效果如表1。5為SST (Space Solar telescope 應用的壓縮實現(xiàn)方法,與之相比,本設計更具有適合航天應用的功能優(yōu)勢,對比如表2。圖10所得火星及月球表面圖像的原圖及解壓圖表1實測壓縮性能壓縮比CR峰值信噪比PSNR時間/ms火星表面

23、833.6337.74月面1839.9537.04月面2838.9238.58表2兩種設計壓縮功能對比算法變換實現(xiàn)平臺處理模式壓縮效果5EZW 5/3整數(shù)FP GA +DSP 不分塊,不抗誤碼參考項設計CCSDS 9/7整數(shù)FP GA 分塊,抗誤碼擴散更優(yōu)64結論本文基于CCSDS 圖像壓縮協(xié)議,采用邊界補償?shù)娜哂喾謮K方法,設計實現(xiàn)了高效的星載圖像壓縮模塊。通過采用并行結合流水的功能結構,提高了在軌壓縮的效率,實驗表明,能夠滿足某些航天應用的需求,說明了設計方法的合理性,也為進一步的改進優(yōu)化提供了經(jīng)驗。參考文獻:1CCSDS.Image Data Compression.CCSDS122.02B 21S.2005,11.2Ritter J ,Moli