張亞飛 馮小華

（中南大學(xué) 信息科學(xué)與工程學(xué)院，湖南 長沙 410012）

圖像去噪問題一直都是圖像處理領(lǐng)域的研究重點(diǎn)， 噪聲的干擾往往會(huì)嚴(yán)重影響圖像的質(zhì)量，給后續(xù)的圖像分析及應(yīng)用帶來極大的不便，為此研究者們一直在尋求一種既能有效地減小噪聲，又能很好地保留圖像重要特征的去噪方法。 與傳統(tǒng)的小波變換去噪方法相比， 提升小波變換以其獨(dú)特的算法結(jié)構(gòu)，降低了計(jì)算復(fù)雜度, 易實(shí)現(xiàn)整數(shù)到整數(shù)的變換等特點(diǎn)，成為小波研究與應(yīng)用領(lǐng)域新的熱點(diǎn)。

1 小波變換

小波分析是一種信號的時(shí)頻分析方法，它具有實(shí)時(shí)、多分辨率、自適應(yīng)等特點(diǎn)，在時(shí)頻兩域都具有很強(qiáng)的表征信號局部特征的能力，因而在圖像去噪領(lǐng)域中得到了廣泛的應(yīng)用。

小波一詞由Morlet 和Grossman 在1980 年代早期提出。傳統(tǒng)的小波通常是由定義在L2(R)空間上的容許函數(shù)ψ（t）經(jīng)過伸縮、平移變換生成的，這種小波通常被稱為第一代小波，其基本定義如下：

如果ψ（t）是一滿足條件的基小波或母小波，將其經(jīng)過伸縮和平移變換后，就可得到一簇函數(shù)，形式如下

此函數(shù)被稱之為小波函數(shù)簇， 其中a 為伸縮因子，b 為平移因子。 進(jìn)而定義下式

為關(guān)于基小波ψ 的連續(xù)小波變換。 其中，ψ*（t）表示對ψ（t）的共軛運(yùn)算。

2 提升小波變換

通過以上分析， 傳統(tǒng)的第一代小波變換是在歐氏空間內(nèi)通過基底的平移和伸縮來構(gòu)建小波基的，這種變換雖然克服了傅立葉變換所存在的固定分辨率、難以同時(shí)兼具時(shí)頻分解特性等問題，但其本身也存在著算法復(fù)雜度高、占用內(nèi)存大、浮點(diǎn)運(yùn)算誤差大、不利于硬件實(shí)現(xiàn)等缺點(diǎn)。 為此，1995 年Sweldens 提出了一種不依賴于傅里葉變換的新的小波構(gòu)造方法—提升方法，稱為第二代小波變換。

第二代小波提升算法基本思想是通過一個(gè)預(yù)測算子，確定高頻信息，并初步確定低頻信息；然后通過更新算子，對初步確定的低頻進(jìn)行修正，從而確定低頻信息。 基于提升方案的第二代小波變換過程可分為三個(gè)階段：分解(Split)、預(yù)測(Predict)和更新(Update)。

（1）分解：將輸入信號sj分為2 個(gè)較小的子集sj-1和dj-1。最簡單的辦法是將輸入信號sj根據(jù)奇偶性分為2 組，這種分解方法所產(chǎn)生的小波也稱為“懶小波”，數(shù)學(xué)表達(dá)式描述如下：

（2）預(yù)測：在基于原數(shù)據(jù)相關(guān)性的基礎(chǔ)上，用偶數(shù)序列sj-1的預(yù)測值P（sj-1）去預(yù)測奇數(shù)序列dj-1，然后用dj-1與預(yù)測值P（sj-1）的差來替換dj-1。 用數(shù)學(xué)公式描述如下

式中為預(yù)測算子， 一般可以通過差值細(xì)分方法構(gòu)造。 這樣，產(chǎn)生的dj-1比原來的dj-1包含更少的信息，繼續(xù)用更小的子集sj-1和dj-1來代替原來的信號集sj， 重復(fù)分解和預(yù)測過程，最終原始信號集可以用｛sn，dn，……，s1,d1｝來表示。

（3）更新：更新的思想主要是要構(gòu)造一個(gè)算子U，利用小波子集dj-1對sj-1進(jìn)行更新， 從而使sj-1繼續(xù)保持原來數(shù)據(jù)序列的某些特征。 數(shù)學(xué)公式描述如下

此外，小波提升變換過程是可逆的，只需改變式（3）（4）（5）的計(jì)算次序和符號就可以得到重構(gòu)的數(shù)據(jù)序列。

3 實(shí)例分析

在Matlab（R2010b）環(huán)境下，以woman 為原始圖像進(jìn)行仿真實(shí)驗(yàn)，分別利用傳統(tǒng)小波變換和本文所介紹的提升小波變換對含噪圖像進(jìn)行處理，仿真結(jié)果如圖1 所示。

圖1 提升小波變換圖像去噪效果圖

4 結(jié)論

仿真實(shí)驗(yàn)結(jié)果可見， 基于提升小波變換的圖像去噪效果與傳統(tǒng)小波變換相當(dāng)， 但考慮到提升小波變換運(yùn)算復(fù)雜度低、實(shí)現(xiàn)過程簡單快速等優(yōu)點(diǎn)，以較小的代價(jià)實(shí)現(xiàn)了預(yù)期的處理效果， 這無疑驗(yàn)證了基于提升小波變換的圖像去噪方法的優(yōu)越性。 同時(shí)，針對提升小波變換的實(shí)現(xiàn)過程，我們可以看到數(shù)學(xué)理論方法的改進(jìn)對于提高運(yùn)算效率的重要性，可以預(yù)見在小波分析領(lǐng)域，根據(jù)實(shí)際情況合理選擇先進(jìn)的數(shù)學(xué)變換方法，對于實(shí)現(xiàn)問題最優(yōu)化解的重要意義。

[1]張德豐，等.MATLAB 小波分析[M].北京:機(jī)械工業(yè)出版社，2011

[2]趙慶光.基于提升小波變換的圖像去噪算法研究[D].哈爾濱：哈爾濱工業(yè)大學(xué)，2009

[3]林祥云.基于提升小波變換的圖像去噪[D].武漢：武漢紡織大學(xué)，2010