O 引言
隨著遙感技術的發(fā)展，遙感圖像的壓縮受到越來越多的關注。靜態(tài)圖像壓縮標準JPEG2000放棄了JPEG標準采用的以離散余弦變換的區(qū)塊編碼方法，采用以小波變換為主的編碼方式，提高了壓縮性能。本文將分析遙感圖像的特點，在JPEG2000無損壓縮基礎上增加了遙感圖像高頻子帶分級量化方案，提高了遙感圖像壓縮比，實現(xiàn)圖像的高保真壓縮。
1 遙感圖像的特點
遙感圖像分辨率較高，目標物體所占的像素較少，圖像中所含的高頻信息和紋理信息較多，空間頻率的局部變化較快，圖像的不連續(xù)性較強。
為了從空間相關性來分析遙感圖像數(shù)據(jù)之間的特征，本文采用圖像的自相關函數(shù)r(l，k)來衡量：
其中f(x，y)為圖像的灰度函數(shù)。對上式進行歸一化和離散化處理，有：
其中，圖像的大小為M×N，uf為圖像亮度的平均值。
當l=k=1時即表示相鄰兩個像素之間的自相關函數(shù)的值R(1，1)，稱為自相關系數(shù)。
在對遙感圖像進行3層小波分解后，除LL子帶集中了主要的低頻信息，其它所有高頻子帶近似遵循零均值Gam—ma分布。分析表明，除極少數(shù)紋理較強的局部區(qū)域之外，絕大多數(shù)紋理不太強的區(qū)域近似服從廣義Gaussian分布，特別是平坦區(qū)域，均值大多約等于0。遙感圖像的這些特性為我們采用JPEG2000中的算法提高圖像壓縮效率提供了理論依據(jù)。
2 遙感圖像高頻子帶分級量化方案
遙感圖像進行JPEG2000整數(shù)小波變換后，編碼系統(tǒng)需要對分量中的每個拼接塊(tile-component)中的系數(shù)進行量化。量化的關鍵是根據(jù)變換后的圖像的系數(shù)特征、重構圖像質量要求等因素設計合理的量化步長。其中高頻子帶影響較小，只作粗量化編碼；低頻子帶對重建過程影響較大，可進行精細量化編碼。
對于遙感圖像，在經過小波變換后90％的能量集中在低頻子帶LLn中，其中包含了遙感圖像的亮度信息，必須進行無失真壓縮，取量化步長為1，進行無損壓縮；高頻子帶集中了圖像的紋理信息，由于遙感圖像中目標物體數(shù)目多，邊緣信息豐富，紋理信息也占有一部分比重。結合低頻子帶精細量化、高頻子帶粗量化的思想，提出了對圖像高頻子帶劃分為4×4的子塊，采用塊內的方差作為塊的紋理復雜度和重要性度量，對紋理復雜的重要塊采用精細量化，而對于較平坦的區(qū)域采用粗量化，達到在高保真條件下提遙感圖像高壓縮比。這實際上相當于在子帶內的分級量化方案，使得各子塊的恢復誤差大致平衡。
(1)基于紋理特征的高頻子帶分級量化方案設計
在遙感圖像進行JPEG2000無損壓縮過程中，圖像經過n級5／3整數(shù)小波變換后，被分解成3n個高頻子帶HHj，HLi和LHi(j=l，2，3，…，n，其中n為分解級數(shù))和一個低頻子帶LL。其中LL子帶的數(shù)據(jù)量雖然不是很大，但卻承載著整個圖像的亮度信息，HHj，HLj和LHj子帶則包含了圖像的高頻信息。因此，對于LL子帶中的數(shù)據(jù)必須進行量化步長為1的無損壓縮，將該子帶中的小波系數(shù)直接送入EBCOT編碼器并進行位平面編碼和算術編碼；高頻子帶中的小波系數(shù)數(shù)據(jù)量雖然較大，卻只是用來描述圖像的細節(jié)信息，可以將該子帶進一步劃分為4×4的塊，采用塊內方差作為塊的紋理復雜度和重要性的度量，區(qū)分為平坦區(qū)域、次紋理區(qū)域和紋理區(qū)域。其中，平坦區(qū)域為子帶中系數(shù)均勻分布，且變化很小的區(qū)域；次紋理區(qū)域包括部分性質均勻區(qū)域和紋理變化適度的區(qū)域；紋理區(qū)域則為系數(shù)紋理性質較強和結構邊緣或邊界區(qū)域。
高頻子帶中平坦區(qū)域小波系數(shù)數(shù)據(jù)量小，進行粗量化后對圖像恢復質量的影響較小，因此取量化步長△3對該區(qū)域小波系數(shù)進行量化；次紋理區(qū)域為紋理區(qū)域到平坦區(qū)域的過渡部分，數(shù)據(jù)量相對比平坦區(qū)域多，因此取量化步長△2對其小波系數(shù)進行量化；紋理區(qū)往往為重要的目標信息，對圖像質量的影響較大，因此必須進行精細量化，取量化步長為△1。經過量化處理后，將新的高頻子帶小波數(shù)劃分成編碼塊(code-block)后送入EBCOT編碼器和算術編碼器中進行編碼，最終形成JPEG2000碼流格式。
在解碼器端，由于采用本方法壓縮后，碼流組織為JPEG2000碼流結構，可以將其直接進行JPEG2000解碼，獲取恢復圖像。
(2)量化步長的選取
量化步長的選取可以根據(jù)子帶的統(tǒng)計特性和最小失真，計算出量化步長的最優(yōu)值；實際應用中可先選取初始量化器，根據(jù)實際圖像壓縮過程，驗證各量化步長與壓縮性能指標PSNR和CR(峰值信噪比和壓縮比)的關系，經過多次調整量化步長最終確定較優(yōu)值。
3 仿真結果及分析
選取衛(wèi)星遙感圖像中截取的512×512、8bit灰度圖像作為測試圖像(如圖1所示)，其中圖像(a)為航拍遙感圖像，(b)為衛(wèi)星遙感圖像。
小波變換濾波器為Daubechies(5／3)濾波器組，進行3級離散小波變換。表1給出的是上面兩幅遙感圖像采用JPEG2000無損壓縮后的測試數(shù)據(jù)。由于無損壓縮，圖像可以完全重構，PSNR趨近于無窮大。所以采用壓縮重構時間、壓縮比為圖像壓縮的評價標準，具體數(shù)據(jù)如表1。表2給出上面兩幅遙感圖像采用高頻子帶分級量化后的測試數(shù)據(jù)。采用分級量化方案對遙感圖像進行壓縮，由于高頻子帶采用大于l的量化步長．勢必產生一定的失真；另外，分級量化后的參數(shù)整數(shù)化處理也引入了部分失真。因此，在評價此方案的壓縮性能時采用壓縮和重構時間、峰值信噪比、壓縮比為評價標準。
由表l和表2可知，與JPEG2000無損壓縮方法相比較，分級量化方法在圖像壓縮性能上有了明顯提高；盡管引入了部分高頻子帶系數(shù)的失真，峰值信噪比有所下降，但是采用本文方案重建的圖像質量均可達到40db，人眼很難分辨出圖像的失真。因此，本方案能滿足遙感圖像壓縮的要求。
采用分級量化方法對遙感圖像進行高保真壓縮和重構，運算時間均有所增加。由于在編碼過程中，將高頻子帶分塊，計算每個子塊的方差并對其進行量化；解碼過程中，需要判斷高頻子帶子塊的紋理屬性，也通過計算方差并與閾值進行比較。因此，此方法在一定程度上增加了算法的復雜度。
4 結論
本文在JPEG2000無損壓縮的框架上，增加了一種遙感圖像高頻子帶分級量化方案，實驗表明其實現(xiàn)簡單，壓縮效果良好?？梢?，該方法在遙感圖像高保真壓縮方面具有良好的效果。
下載文檔
