??? 摘 要: 分形圖像壓縮編碼" title="壓縮編碼">壓縮編碼是近年來產生的新的圖像壓縮編碼技術，由于其具有極高的壓縮比而獲得廣泛的關注。主要討論了圖像小波" title="小波">小波域的去噪" title="去噪">去噪問題以及如何將小波域的去噪與分形圖像壓縮方法結合起來，以獲得良好的編碼效率和圖像質量。
??? 關鍵詞: 分形編碼? 圖像壓縮? 噪聲
?
??? 分形的概念是由數學家B.Mandelbrot于1975年提出的，他把分形定義為“一種由許多個與整體有某種相似性的局部所構成的形體”。分形概念的提出及分形幾何學的創立為描述客觀世界提供了更準確的數學模型。圖形學是幾何學的延伸與發展，分形模型研究成果的積累形成了新的圖像學分支——分形圖像學。而基于分形的圖像編碼方法實質是對圖像中一個或多個相對大的部分施行壓縮變換來逼近圖像的每一部分。1990年，A.Jacquin提出了全自動的可行的分形壓縮編碼方法，由于其可以獲得極高的壓縮比而得到廣泛關注。在實際的圖像編碼過程中，原始圖像經常被噪聲（最常見的是高斯白噪聲）污染。由于噪聲的存在，一方面使得圖像編碼的時間延長，另一方面，降低了圖像的信噪比，圖像質量明顯下降。因此，筆者試圖尋找一種方法，在不影響圖像壓縮比的前提下去噪，從而提高編碼效率和圖像質量。
1 分形圖像壓縮編碼方法
1．1拼接定理（Collage Theorem）
??? 拼接定理是分形圖像壓縮技術的核心[1]。
??? 設{RT: wi,i=1,2,…，p}是T維的收縮仿射變換集合，
收縮因子s∈(0,1),且滿足：
??? 
??? 其中，A為IFS的吸引子，h(A,B)為Hausdorff距離。拼接定理給出了數集V與吸引子IFS之間逼近程度的一個上限值，即拼接誤差的上限值。
??? 拼接定理提供了用IFS進行圖像壓縮的理論依據。對于一般的灰度圖像，可認為是一張原始灰度曲面（R3空間上的一個緊子集）進行抽樣和量化得到的。盡管無法使原始圖像(V)成為某一個迭代函數系統IFS的吸引子，但是如果能找到一組收縮仿射變換wi,i=1,2,…，p，
該IFS{RT:wi,i=1,2,…，p}對應吸引子的良好逼近。
??? 在利用收縮仿射變換{RT:wi,i=1,2,…p}對圖像進行解壓縮時，迭代過程與初始條件無關。也就是說，對任意給定的初始圖像數據進行多次迭代，就可以完成對原始圖像的重構" title="重構">重構。
1．2 分形圖像壓縮編碼的實現
??? 所謂局部IFS（LIFS）是指其變換的定義域" title="定義域">定義域由原來的整個區域放寬為全部區域的某些子集。通過將IFS理論從全局擴展到局部，可以得到一種全自動的分形壓縮方案，Jacquin方法的實質是固定方塊定義域塊的大小并限制仿射變換為一定的形式，然后搜索圖像尋找與定義域塊相匹配的值域塊（它的大小是定義域塊的4倍），搜索中要配合Jacquin提出的八種對稱變換算子對值域塊進行變換。
??? 編碼壓縮過程：把原始圖像分成互不重疊的定義域塊(Range塊)，這些塊能將原始圖像全部覆蓋，每一個Range塊大小均為B×B；再把原始圖像分成互相重疊的值域塊（Domain塊），每一Domain塊的大小均為D×D,且通常D=2B。為使壓縮后重構圖像的質量更好，相鄰的Domain塊之間在水平及垂直方向均有重疊，水平及垂直方向上位移量為B。依次對每一個Domain塊中相鄰的4個灰度值求平均，于是每一個大小為D×D(2B×2B) Domain塊就變成了大小為B×B的Sub_ Domain塊。之后利用最小二乘法，并配合八種對稱變換算子，將收縮后的Domain塊（Sub_ Domain塊）與Range塊進行匹配運算：
??? 
??? 若計算出的MSE小于給定誤差，則認為匹配成功，否則繼續匹配，在找到最佳值域塊及仿射變換后，需要存儲其參數以便傳輸。這些參數包括：匹配成功的Domain塊塊號、Jacquin變換算子的編號以及比例因子S、偏移量O。
??? 解碼重構過程：從任意的初始圖像開始（當然該圖的大小要與原圖一致），和原始圖像（如圖1所示）一樣，對任意圖像也分為相同大小的定義域塊和值域塊，根據傳輸的參數，用值域塊對定義域塊進行逼近，當全部的定義域塊都被逼近一次后，稱為完成一次迭代，所得圖像作為下一次迭代的初始圖像，一直迭代下去直至重構圖像不再隨迭代發生顯著變化為止，重構結束。
2 對噪聲圖像進行分形壓縮編碼的結果
??? 以256×256的Lena灰度圖像為例（噪聲均為零均值的高斯白噪聲）,在對圖像進行分割時，取B=4,D=8。表1給出了用分形壓縮算法對噪聲圖像編碼的結果。圖2和圖3給出了噪聲方差為50時的噪聲圖像和分形壓縮圖像。
?
?
??? 從表1中可以明顯看出，當圖像被高斯白噪聲污染后，分形編碼算法的性能明顯降低，一是編碼時間大大增加，而且噪聲方差越大，編碼時間越長，這主要是因為圖像被污染的程度越大，圖像信息損害的就越厲害，圖像的分形特性遭到破壞，在編碼過程中尋找匹配塊的難度就越大、時間越久，分形算法的編碼效率大大降低。二是重構圖像信噪比降低，圖像質量嚴重惡化，尤其是方差增大時，圖像質量惡化更嚴重，如方差為200時，重構圖像的信噪比降到了24.36（如圖3所示），圖像的一些信息特征已經不存在了。因此要提高含噪聲圖像的編碼效率和信噪比，就必須對圖像先進行去噪處理，同時要保證圖像的分形特點。由于小波已被廣泛地用于圖像處理并獲得良好的效果，所以利用小波方法對圖像進行去噪處理。
3? 圖像小波域的去噪方法
??? 小波分析用于圖像去噪處理，主要是針對圖像信號與噪聲信號經小波變換后在不同的分辨率下呈現不同的規律，在不同的分辨率下，設定閾值門限，調整小波系數，達到圖像去噪目的[2]。
??? 在小波系數進行取舍之前，實際上按照一定準則（或者閾值化）將小波系數劃分為兩類：一類是重要的、規則的小波系數；另一類是被看作非重要的或者受噪聲干擾較大的小波系數。通常以小波系數的絕對值作為小波系數的分類單元。小波系數絕對值趨向于零，意味著小波系數所包含的信息量少并且強烈地受噪聲干擾。最常用的閾值化去噪方法：一是默認閾值消噪處理，即在消噪處理過程中采用程序中設定的閾值，對分解信號進行分類處理，以求消除噪聲；二是給定軟（或硬）閾值消噪處理，閾值通過某一個經驗公式獲得，該閾值比默認的閾值去噪效果更有說服力。
??? 對于“軟閾值化”，絕對值小于閾值?啄的小波系數數值用零代替；絕對值大于閾值?啄的小波系數數值用?啄來縮減。如下式所示:
??? 
??? 式中，W表示小波系數的數值；sgn(·)是符號函數，當數值大于零，符號為正，反之符號為負。對于“硬閾值化”，僅僅保留絕對值大于閾值?啄的小波系數，并且被保留系數與原始系數相同（沒有被縮減）,如下式所示:
??? 
??? 兩種方法各有差異，前者具有連續性，在數學上容易處理，后者更接近實際應用。
??? 閾值化處理的關鍵在于選擇合適的閾值?啄。如果閾值門限太小，處理后的信號仍有噪聲存在；閾值太大，重要的圖像特征將被濾掉，引起偏差。大多數閾值的選取過程是針對一組小波系數，即根據本組小波系數的統計特性，計算出一個閾值?啄。Donoho[3]等提出了一種典型的閾值選取方法，從理論上證明閾值與噪聲的方差成正比，為：
??? 
??? 其中,Nj表示第j層子帶上小波系數的個數。
??? 通過對信號的收縮閾值處理，能較好去除噪聲，提高重構圖像的信噪比和圖像的編碼效率。
4 去噪聲圖像的分形壓縮編碼結果及討論
??? 為了便于對比，這里所采用的圖像仍然是256×256的Lena灰度圖像，分解所用的濾波器是雙正交B97濾波器,分解級數為4。實驗均在Pentium 4/1.5G、RAM128MB的機器上編譯完成。表2給出了用分形壓縮算法對消噪圖像編碼的結果。圖4和圖5是消噪圖像（原噪聲圖像方差=100）和消噪圖像分形壓縮結果。
?
?
??? 從表2和表1的對比中可以明顯看出，經小波域去噪的圖像分形壓縮編碼性能明顯提高，一是體現在編碼時間上，去噪圖像的編碼時間與有噪圖像相比大大縮短，效率大大提高，分形壓縮方法本身特點就是壓縮時間長，因此這種效率的提高是非常可觀的；二是圖像質量明顯好轉，信噪比有較大的提高，保證了分形壓縮編碼的質量（如圖5所示）。以上實驗和結論充分說明了圖像小波域的去噪方法能夠對噪聲圖像進行良好的消噪處理，保證了分形壓縮編碼方法的編碼效率和圖像質量。
參考文獻
1 杜世培,楊 木.分形圖像壓縮方法的改進與實現[J].貴州工業大學學報(自然科學版),2000；（8）：56~59
2 陳武凡.小波分析及其在圖像處理中的應用[M].北京：科學出版社,2002
3 David L.Donoho .De-noising by soft-thresholding[J].IEEE?Trans.On Information Theory.1995；41(3)：617~627