章浙濤，朱建軍，匡翠林，周 璀

中南大學(xué)地球科學(xué)與信息物理學(xué)院，湖南長沙 410083

小波包多閾值去噪法及其在形變分析中的應(yīng)用

章浙濤，朱建軍，匡翠林，周 璀

中南大學(xué)地球科學(xué)與信息物理學(xué)院，湖南長沙 410083

在形變分析中，傳統(tǒng)的小波去噪只保留低頻上的有用信息，很容易去掉中頻以及高頻上的有用信息。小波包分析方法則同時考慮了各個頻段上的有用信息，因此是一種更為精細(xì)的去噪方法，小波包去噪的關(guān)鍵是對小波包分解系數(shù)選取合適的閾值準(zhǔn)則并進(jìn)行閾值處理，但傳統(tǒng)的小波包去噪并沒有對此進(jìn)行充分的研究。針對傳統(tǒng)小波、小波包分析的不足，本文提出一種基于頻率順序并依據(jù)信息類型分段的多閾值準(zhǔn)則小波包去噪法。通過理論分析與實際應(yīng)用，結(jié)果表明該方法能夠高效剔除各頻段的噪聲，同時當(dāng)采樣頻率較低時能有效保留去噪信號中頻率較高的有用信息，其去噪能力優(yōu)于傳統(tǒng)的小波、小波包等其他去噪方法，因此可以廣泛應(yīng)用于高精度變形監(jiān)測領(lǐng)域中。

小波包分析；變形監(jiān)測；數(shù)據(jù)去噪；頻率順序；多閾值準(zhǔn)則

1 引 言

小波分析因其具有良好的時頻特性，已被廣泛應(yīng)用于形變分析中，并取得了不錯的效果［1］。例如對高層建筑的動態(tài)監(jiān)測［2-3］，礦層的移動分析［4］以及大壩的變形預(yù)報［5］等。但小波分析主要是提取低頻有用信息，容易忽略中頻以及高頻有用信息，這會對高精度GPS變形監(jiān)測的分析結(jié)果產(chǎn)生不良影響。例如對高層建筑進(jìn)行變形監(jiān)測，不僅要監(jiān)測出其主模態(tài)的振動頻率，有時還需提取其風(fēng)致振動變形、臺風(fēng)暴雨等極端天氣引起的變形信息，這些變形信息的振動頻率已高達(dá)甚至超過10－1Hz［6］，當(dāng)采樣頻率較低時傳統(tǒng)的小波分析就容易產(chǎn)生過度去噪的現(xiàn)象，無法準(zhǔn)確提取有用信息。小波包分析是一種新的小波分析方法，具有更好的去噪能力，它不僅像小波分析一樣對觀測數(shù)據(jù)的低頻部分進(jìn)行分解，還能對觀測數(shù)據(jù)的高頻部分進(jìn)行分解，能相對有效地提取各頻段的有用信息［7］。目前，小波包分析已在數(shù)據(jù)去噪［8］、振動信息提?。?］等變形監(jiān)測領(lǐng)域廣泛應(yīng)用。其中小波包應(yīng)用的關(guān)鍵是對小波包分解系數(shù)選取合適的閾值準(zhǔn)則并進(jìn)行閾值處理，現(xiàn)今已有了諸如改進(jìn)的單閾值處理［10］和普通的多閾值處理［11-12］的方法，較傳統(tǒng)小波包去噪的效果有一定改善，但是上述方法并沒有完全充分考慮信號及其噪聲的分布，去噪精度有待進(jìn)一步提高。

如何對小波包分解系數(shù)排序重組并進(jìn)行多閾值處理是高精度去噪的關(guān)鍵問題。文獻(xiàn)［11］的方法是對小波包樹的第一個節(jié)點和最后一個節(jié)點進(jìn)行了特殊化處理，雖然提高了去噪精度，但是事實上最后一個節(jié)點并不是頻率最高部分，如果不按照小波包樹節(jié)點的頻率順序，而利用小波包樹分解的自然順序進(jìn)行多閾值處理，雖然能夠?qū)Ω哳l部分的噪聲進(jìn)行相對有效地處理，但卻是一種比較粗糙的多閾值處理方法。文獻(xiàn)［12］的方法是依據(jù)小波包分解系數(shù)的能量來進(jìn)行分組，同樣提高了去噪精度，并且比文獻(xiàn)［11］的方法更為精確，但是當(dāng)各個頻段中噪聲污染不穩(wěn)定時，即噪聲污染比較復(fù)雜，已嚴(yán)重影響了信號的能量分布，這時候根據(jù)能量來進(jìn)行分組會對有用信息的提取造成困難。

基于這些問題，本文提出了一種小波包去噪方法，即根據(jù)不同信號及其噪聲的分布，對小波包分解系數(shù)按照頻率大小的順序進(jìn)行排列，依據(jù)信息類型分段，對每個頻段選取合適的閾值準(zhǔn)則并進(jìn)行閾值處理。它是一種基于頻率順序的多閾值準(zhǔn)則（the multi-threshold criteria based on frequency order）小波包去噪法，稱為FMC法。理論和實例證明，該方法去噪能力優(yōu)于傳統(tǒng)方法，并且能夠有效提取信號中各頻段的有用信息，可以高效處理變形監(jiān)測數(shù)據(jù)。

2 小波包多閾值去噪法的關(guān)鍵技術(shù)

2.1 小波包分解系數(shù)的排序重組

小波包理論的分解算法和重構(gòu)算法的數(shù)學(xué)表達(dá)式如下［9］。

式中，h、g為濾波器系數(shù)；d為小波包分解系數(shù)；p、t為分解層數(shù)；j、k為小波包節(jié)點號。

由于噪聲在頻率域上的分布主要集中在頻率較高的部分［13-14］。因此在利用小波包去噪時，假如不同頻段信息采用同一種閾值處理方法，或者同一個頻段信息采用不同閾值處理方法勢必會影響去噪精度。因此如何準(zhǔn)確劃分頻段是進(jìn)行多閾值去噪的關(guān)鍵問題。這里提出了一種按照頻率大小順序排列，并依據(jù)信息類型重組的劃分方法。根據(jù)小波包變換理論，它可以將隱含在低頻、中頻、高頻中的有用信息提取出來，以分3層為例，信號分解如圖1（a）所示。以采樣頻率為2 Hz的GPS變形監(jiān)測數(shù)據(jù)為例，同時根據(jù)那奎斯特（Nyquist）采樣定理，將信號用小波包分解如圖1（b）所示。然后再將第3層每一個小波包樹節(jié)點與圖1（a）進(jìn)行逐一對應(yīng)，得表1。

圖1 小波包3層分解節(jié)點號及頻率分解示意圖Fig．1 Three layers of decomposition node numbers and frequency decompositions with wavelet packet

通過表1可以發(fā)現(xiàn)，小波包樹節(jié)點的自然順序與頻率順序存在不一致現(xiàn)象。最低頻部分對應(yīng)的是第1個節(jié)點，最高頻部分對應(yīng)的是第5個節(jié)點。出現(xiàn)這種現(xiàn)象的原因是由于小波包分解時，高通濾波器會進(jìn)行一次“翻轉(zhuǎn)”操作?？梢酝ㄟ^小波包的性質(zhì)證明，任何小波包的分解都會產(chǎn)生自然順序與頻率順序不一致的現(xiàn)象，并且不一致的現(xiàn)象的情況是一樣的［15］。因此，小波包每分解一層時，低頻分解部分按頻率從小到大排列，高頻部分按頻率從大到小排列。在變形監(jiān)測領(lǐng)域，因為噪聲分布和頻率順序緊密相關(guān)，因此需要的是頻率順序而不是自然順序。

表1 第3層小波包樹節(jié)點號與小波包分解頻率對應(yīng)表Tab.1 The third layer of wavelet packet tree nodes corresponding to the frequency of wavelet packet decompositions

在形變分析中，往往有需要集中保留的頻段以及集中去噪的頻段。例如圖1（b）情況的分解，如果需要監(jiān)測0．15 Hz的建筑物主模態(tài)振動頻率信息，根據(jù)計算，［3，0］、［3，1］為有用信息集中頻段，即為低頻段；而噪聲污染情況一般，［3，4］為無用信息集中頻段，即為高頻段；其余節(jié)點為過渡頻段，即為中頻段。因此，可以根據(jù)需要監(jiān)測的頻段以及噪聲污染情況，對小波包分解系數(shù)進(jìn)行分組，即為依據(jù)信息類型重組的方法。

具體而言，根據(jù)信號及其噪聲的分布，以分3層為例，設(shè)α為低頻段比例系數(shù)，β為高頻段比例系數(shù)，n為信號分解層數(shù)。則有小波包樹節(jié)點分組如（3）所示。需要指出的是，在某些特殊情況下，也可以對頻段進(jìn)一步細(xì)分，分為低頻、次低頻、中頻、次高頻以及高頻，但基本原理不變。

式中，小波包樹節(jié)點［n，m0］、［n，m1］、…、［n，m2n－1］按頻率順序從小到大進(jìn)行排列。

2.2 小波包分解系數(shù)的多閾值處理

單一閾值準(zhǔn)則由于沒有充分考慮信號及其噪聲的分布，因此去噪時很容易出現(xiàn)去噪效果不明顯或是去噪過度等現(xiàn)象，而多閾值準(zhǔn)則就很好地克服了上述缺點。多閾值準(zhǔn)則的基本思想就是根據(jù)一定法則，對每個小波包分解系數(shù)靈活的選用不同閾值準(zhǔn)則，從而最大程度地保留有用信息，去掉無用信息。小波包分析中常用的四種閾值準(zhǔn)則包括固定形式閾值準(zhǔn)則（sqtwolog）、自適應(yīng)閾值準(zhǔn)則（rigrsure）、啟發(fā)式閾值準(zhǔn)則（heursure）以及極小化極大閾值準(zhǔn)則（minimaxi）［16］，由于各自的選取規(guī)則不一，去噪效果不盡相同。因此，每一個小波包分解系數(shù)都有它最合適的閾值準(zhǔn)則。

（1）固定形式閾值準(zhǔn)則。它所采用的是固定形式的閾值，其公式為

式中，σ為噪聲信號的標(biāo)準(zhǔn)差；N為信號長度。

（2）自適應(yīng)閾值準(zhǔn)則。它是基于Stein的無偏似然估計原理的一種閾值，其公式為

式中，σ為噪聲信號的標(biāo)準(zhǔn)差；Qa為依據(jù)無偏似然估計原理而得到的小波包分解系數(shù)的平方。

（3）啟發(fā)式閾值準(zhǔn)則。它是一種無偏風(fēng)險閾值，是前兩種閾值準(zhǔn)則的綜合，選取規(guī)則為

（4）極小化極大閾值準(zhǔn)則。也是一種固定形式的閾值，是在最壞條件下求最大均方誤差的最小值，這在統(tǒng)計學(xué)上叫做“悲觀法”，選取規(guī)則為

式中，σ為噪聲信號的標(biāo)準(zhǔn)差；N為信號長度。

4種閾值準(zhǔn)則都有其各自特點。其中，sqtwolog閾值準(zhǔn)則以及heursure閾值準(zhǔn)則兩種方法類似，都是將全部系數(shù)進(jìn)行了處理［17-18］，因此可以較強(qiáng)地去除噪聲，與之相對應(yīng)的，容易過度去噪，因而可以稱之為“激進(jìn)”的去噪準(zhǔn)則，適合處理GPS信號的高頻部分。rigrsure閾值準(zhǔn)則以及minimaxi閾值準(zhǔn)則是將部分系數(shù)進(jìn)行了處理［17-18］，是一種較為折中的處理方法，因此，它可以防止過度去噪，與之相對應(yīng)的容易出現(xiàn)去噪不明顯的現(xiàn)象，因而可以稱之為“保守”的去噪準(zhǔn)則，適合處理低頻部分。其中，中頻段更適合采用rigrsure閾值準(zhǔn)則以及minimaxi閾值準(zhǔn)則［11，15，19］。

為此，建立如表2所示的多閾值準(zhǔn)則選取表，以用于對小波包分解系數(shù)的閾值處理。

表2 多閾值準(zhǔn)則選取表Tab.2 The selections of multi-threshold criteria

2.3 算法步驟

基于頻率順序并依據(jù)信息類型進(jìn)行分段的小波包多閾值準(zhǔn)則處理方法是一種精度較高的去噪方法。下面給出其算法步驟，以分高中低3個頻段為例。

（1）小波包分解。根據(jù)所給信號選定一個合適的小波基函數(shù)φ以及決定所需的分解層次n，對信號進(jìn)行小波包分解，得到小波包樹。

（2）優(yōu)化小波包樹（可選）。確定熵標(biāo)準(zhǔn)從而選定最優(yōu)小波基，得到最優(yōu)樹。

（3）根據(jù)信號及其噪聲的分布，選定低頻段比例系數(shù)α和高頻段比例系數(shù)β。其中，有用信息集中的頻段越多，α越大；無用信息集中的頻段越多，β越大。通過排序重組，得到各頻段的小波包樹節(jié)點信息。

（4）根據(jù)表2的選取標(biāo)準(zhǔn)，決定低頻段閾值準(zhǔn)則A、中頻段閾值準(zhǔn)則B、高頻段閾值準(zhǔn)則C，對各頻段的小波包分解系數(shù)分別進(jìn)行閾值處理。

（5）小波包重構(gòu)。對處理后的小波包分解系數(shù)進(jìn)行小波包重構(gòu)，最終得到去噪信號。

3 仿真試驗

為了驗證此法的有效性和優(yōu)越性，進(jìn)行了一系列的對比分析。以Matlab為平臺，使用評價小波去噪效果質(zhì)量的兩種專用測試信號：blocks信號和doppler信號，在這些信號上加上高斯白噪聲，信噪比為10 d B，信號長度為1024個采樣點，采樣頻率為1 Hz。這兩種信號都為非平穩(wěn)信號，前者具有突變性，后者具有漸變性，因此基本涵蓋了變形監(jiān)測數(shù)據(jù)的特點。其加噪信號如圖2所示。選用均方根誤差（RMSE）、信噪比（SNR）以及與真實信號的相關(guān)系數(shù)（r）為評價指標(biāo)［20］，用來定量比較去噪效果，其中均方根誤差越小、信噪比越大、相關(guān)系數(shù)越高則去噪效果越好。

圖2 兩種加噪信號Fig．2 Two noisy signals

為了較為全面地認(rèn)識基于頻率順序的多閾值準(zhǔn)則（FMC）小波包去噪法的去噪效果，將上述測試信號分別用小波通用去噪，小波包強(qiáng)制去噪，小波包單一全局閾值準(zhǔn)則去噪（sqtwolog，rigrsure，heursure，minimaxi），文獻(xiàn)［11］的方法以及FMC法進(jìn)行處理，每種方法得到對應(yīng)的RMSE、SNR以及r。上述各方法若涉及選取小波基的，均使用工程實際中被廣泛應(yīng)用的db4小波基，db4小波基具有正交性，可以任意分解，另外還有保證高頻信息不丟失的高度緊支撐性［21］。根據(jù)信噪比，分解層數(shù)定為3層。涉及閾值類型的選取時都采用軟閾值，這樣可以更好地保護(hù)有用信息。參數(shù)的統(tǒng)一還可以保證去噪效果僅與分解方法（小波與小波包）、閾值準(zhǔn)則以及頻率順序有關(guān)。使用小波通用去噪的目的是作為一個基準(zhǔn)，因為此法具有簡單、快速、穩(wěn)定的特點，在實際工程應(yīng)用中被廣泛使用。使用小波包強(qiáng)制去噪的目的是為了比較沒有帶閾值準(zhǔn)則的處理與帶閾值準(zhǔn)則的處理的差別。分別利用4種閾值準(zhǔn)則進(jìn)行小波包單一全局閾值準(zhǔn)則去噪是為了比較單一閾值準(zhǔn)則的處理與多閾值準(zhǔn)則的處理的差別。使用文獻(xiàn)［11］的方法去噪是為了驗證基于頻率順序并依據(jù)信息類型排序重組的有效性，其中，節(jié)點［3，0］采用minimaxi準(zhǔn)則，節(jié)點［3，7］采用sqtwolog準(zhǔn)則，其余節(jié)點采用rigrsure準(zhǔn)則。

圖3是利用FMC法對兩種加噪信號進(jìn)行去噪的結(jié)果。具體過程如下：選定db4為小波基函數(shù)，分解層數(shù)為3，對信號進(jìn)行小波包分解，得到小波包樹；熵標(biāo)準(zhǔn)定為Shannon，優(yōu)化小波基，得到最優(yōu)樹；根據(jù)信號及其噪聲的分布，設(shè)低頻段比例系數(shù)α＝1／8，高頻段比例系數(shù)β＝1／8，重組小波包樹節(jié)點，具體而言，低頻段節(jié)點為［3，0］，高頻段節(jié)點為［3，4］，其余節(jié)點分在中頻段；依據(jù)表2，選低頻段閾值準(zhǔn)則A為minimaxi、中頻段閾值準(zhǔn)則B為rigrsure、高頻段閾值準(zhǔn)則C為sqtwolog，分別對小波包樹節(jié)點系數(shù)進(jìn)行閾值處理，具體而言，節(jié)點［3，0］采用minimaxi準(zhǔn)則，節(jié)點［3，4］采用sqtwolog準(zhǔn)則的，其余節(jié)點采用rigrsure準(zhǔn)則；對處理后的小波包分解系數(shù)進(jìn)行小波包重構(gòu)，得到去噪信號。

圖3 兩種去噪信號（FMC法）Fig．3 Two de-noising signals（FMC）

將各種去噪方法后得到的信號的均方根誤差、信噪比以及相關(guān)系數(shù)列表如表3所示。

表3 各種方法去噪效果比較Tab.3 The result comparison between different de-noising methods

由表3可知，F(xiàn)MC法去噪后的均方根誤差最小，信噪比最大，相關(guān)系數(shù)最高，因此去噪效果最好。通過分析可以得到以下幾點結(jié)論：所有方法都可以不同程度地去噪，因此可以應(yīng)用于GPS變形監(jiān)測領(lǐng)域；去噪效果最差的是小波包強(qiáng)制去噪，原因是因為沒有利用到閾值準(zhǔn)則，這在調(diào)整系數(shù)時，會產(chǎn)生很多誤差，因此必要的閾值準(zhǔn)則是提高去噪效果的有效方法；小波包單一全局閾值準(zhǔn)則去噪效果優(yōu)于小波通用去噪，說明小波包去噪的精度整體上是優(yōu)于小波去噪的；4種閾值準(zhǔn)則的去噪效果相當(dāng)，對于不同信號，各有優(yōu)劣，如何正確認(rèn)識信號以及閾值準(zhǔn)則的特點并靈活使用，是一個值得注意的問題；FMC法和文獻(xiàn)［11］的方法去噪效果優(yōu)于小波包單一全局閾值準(zhǔn)則去噪，說明多閾值準(zhǔn)則處理能夠提高去噪精度；FMC法去噪效果優(yōu)于文獻(xiàn)［11］的方法，證明基于頻率順序并依據(jù)信息類型排序重組可以切實提高精度。

4 工程實例

為了驗證該法在實際工程中的去噪效果，進(jìn)行了一組試驗。數(shù)據(jù)是從某教學(xué)樓頂連續(xù)觀測所得的GPS數(shù)據(jù)，基線長度約為13 m，數(shù)據(jù)采樣率為1 Hz，連續(xù)觀測7 d，經(jīng)雙差解算、多路徑重復(fù)性周期對齊之后等預(yù)處理而得到y(tǒng)方向的坐標(biāo)殘差數(shù)據(jù)。這7 d數(shù)據(jù)具有高度相關(guān)性，其中存在著多路徑效應(yīng)、建筑物主模態(tài)的振動頻率，以及風(fēng)致振動變形等信息，所有這些有用信息的振動頻率達(dá)到10－1Hz數(shù)量級。

利用FMC法對該GPS變形監(jiān)測數(shù)據(jù)進(jìn)行去噪。根據(jù)原始信號特點，選定具有正交性以及保證高頻信息不丟失的高度緊支撐性的db4為小波基函數(shù)，熵標(biāo)準(zhǔn)定為shannon；根據(jù)小波包分解性質(zhì)以及那奎斯特采樣定理，采樣頻率為1 Hz，所需要的小波包分解系數(shù)分辨率為0．01 Hz，分解層數(shù)為5時分辨率為0．015 625 Hz，符合要求，確定分解層數(shù)為5；根據(jù)原理，0．015 625為每一個小波包分解系數(shù)的頻段大小，而有用信息高達(dá)10－1Hz，因此低頻段和中頻段的分界點為0．1 Hz，即低頻段比例系數(shù)α＝8／32，又噪聲污染情況一般，即高頻段比例系數(shù)β＝4／32；依據(jù)表2，選定閾值準(zhǔn)則A為minimaxi、B為rigrsure、C為sqtwolog。具體而言，節(jié)點［5，0］、［5，1］、［5，3］、［5，2］、［5，7］、［5，6］、［5，4］、［5，5］采用minimaxi準(zhǔn)則，節(jié)點［5，23］、［5，22］［5，20］、［5，21］采用sqtwolog準(zhǔn)則，其余節(jié)點采用rigrsure準(zhǔn)則。以第1天為例，原始信號和去噪信號如圖4所示，從結(jié)果可以看出去噪后信號的噪聲得到了有效剔除。

此外，還使用小波通用去噪、文獻(xiàn)［11］的方法對該數(shù)據(jù)進(jìn)行了去噪。并計算了第一天與后面各天坐標(biāo)殘差間的原始相關(guān)系數(shù)與各種方法得到的新的相關(guān)系數(shù)r，如表4所示；表5同時計算了各天坐標(biāo)殘差的標(biāo)準(zhǔn)偏差。通過分析能夠發(fā)現(xiàn)，所有經(jīng)去噪后的信號中，F(xiàn)MC法處理后的相關(guān)系數(shù)提高最明顯，標(biāo)準(zhǔn)偏差最小，由此最大程度地去除了噪聲的污染，去噪效果最好。

表4 第1天與后面各天坐標(biāo)殘差間的相關(guān)系數(shù)表Tab.4 Correlation coefficients between the first coordinate residual and others

表5 各天坐標(biāo)殘差的標(biāo)準(zhǔn)偏差表Tab.5 Standard deviations of coordinate residuals mm

為更直觀地說明本方法去噪效果好的原理，圖5給出了利用welch平均周期圖法得到的原始信號、FMC法去噪信號以及傳統(tǒng)小波去噪信號的功率譜密度估計圖。由圖5能夠看出，兩種去噪信號的功率譜密度在10－2Hz以下的信息都得到了很好的保留，已基本滿足去噪需要；而在大于10－1Hz的信息部分，F(xiàn)MC法的去噪之后的功率譜從形狀以及減弱情況來看也比小波去噪要好。通過圖5放大圖可以發(fā)現(xiàn)，在高精度GPS變形監(jiān)測所需要保留的10－2～10－1Hz的信息，傳統(tǒng)小波去噪之后的功率譜密度與原信號的功率譜密度相差較大，而FMC法則更接近原信號，保留了有用信息。因此，在高精度GPS變形監(jiān)測中，F(xiàn)MC法去噪能夠防止去噪不明顯或去噪過度等現(xiàn)象，是一種穩(wěn)定的去噪方法。

圖5 功率譜密度圖Fig．5 Power spectrum density

5 結(jié) 論

通過試驗分析與工程應(yīng)用，可以發(fā)現(xiàn)本文方法去噪效果優(yōu)于傳統(tǒng)的小波、小波包去噪方法以及其他改進(jìn)的小波包去噪方法。本文的創(chuàng)新點如下：①傳統(tǒng)小波、小波包去噪并沒有充分考慮信號及其噪聲的分布，而基于頻率順序并依據(jù)信息類型劃分的多閾值準(zhǔn)則處理方法，可以在剔除噪聲的同時有效保護(hù)各頻段的有用信息；②基于頻率順序并依據(jù)信息類型劃分的觀點用于小波去噪，這為小波、小波包的多閾值準(zhǔn)則去噪提供了一種思路。

［1］ HUANG Shengxiang，LIU Jingnan．A Novel Method for Reducing Noises in GPS Deformation Monitoring System［J］．Acta Geodaetica et Cartographica Sinica，2002，31（2）：104-107．（黃聲享，劉經(jīng)南．GPS變形監(jiān)測系統(tǒng)中消除噪聲的一種有效方法［J］．測繪學(xué)報，2002，31（2）：104-107．）

［2］ H U A NG Shengxiang，LIU Jingnan，LIU Xianglin．Deformation Analysis Based on Wavelet and Its Application in Dynamic Monitoring for High-rise Buildings［J］．Acta Geodaetica et Cartographica Sinica，2003，32（2）：153-157．（黃聲享，劉經(jīng)南，柳響林．小波分析在高層建筑動態(tài)監(jiān)測中的應(yīng)用［J］．測繪學(xué)報，2003，32（2）：153-157．）

［3］ KUANG Cuilin，DAI Wujiao．Measurement of Windinduced Vibration of Tall Buildings Using GPS and Wavelet Application［J］．Geomatics and Information Science of Wuhan University，2010，35（9）：1024-1028．（匡翠林，戴吾蛟．GPS監(jiān)測高層建筑風(fēng)致振動變形及小波應(yīng)用［J］．武漢大學(xué)學(xué)報：信息科學(xué)版，2010，35（9）：1024-1028．）

［4］ GU Qinghua，LU Caiwu，GUO Jinping，et al．Dynamic Management System of Ore Blending in an Open Pit Mine Based on GIS／GPS／GPRS［J］．Mining Science and Technology，2010，20（1）：132-137．

［5］ MA Pan，MENG Lingkui，WEN Hongyan．Kalman Filtering Model of Dynamic Deformation Based on Wavelet Analysis［J］．Geomatics and Information Science of Wuhan University，2004，29（4）：349-353．（馬攀，孟令奎，文鴻雁．基于小波分析的Kalman濾波動態(tài)變形模型研究［J］．武漢大學(xué)學(xué)報：信息科學(xué)版，2004，29（4）：349-353．）

［6］ H U ANG Dingfa，CH EN Yongqi，DING Xiaoli，et al．Wavelet Based Analysis Technique for the Monitoring of Tall Stracture under Normal Loading Using GPS［J］．Journal of Vibration and Shock，2001，20（1）：12-17．（黃丁發(fā)，陳永奇，丁曉利，等．GPS高層建筑物常荷載振動測試的小波分析［J］．振動與沖擊，2001，20（1）：12-15．）

［7］ DAUBECHIES I．Ten Lectures on Wavelets［M］．SIAM：Society for Industrial and Applied Mathematics，1992：312-313．

［8］ YI Ting Hua，LI Hongnan，GU Ming．Experimental Assessment of High-rate GPS Receivers for Deformation Monitoring of Bridge［J］．Measurement，2012，46：420-432．

［9］ WU Jizhong，HUA Xianghong，GAO Junqiang．Feature Extraction of Structure Natural Vibration and Multipath Separation Based on Wavelet Packet Decomposition［J］．Geomatics and Information Science of Wuhan University，2010，35（4）：486-490．（吳繼忠，花向紅，高俊強(qiáng)．基于小波包分解的結(jié)構(gòu)自振特征提取及多路徑誤差分離［J］．武漢大學(xué)學(xué)報：信息科學(xué)版，2010，35（4）：486-490．）

［10］ LI Qingwu，HE Chunyuan．A New Thresholding Method in Wavelet Packet Analysis for Image De-noising［C］∥Proceedings of IEEE International Conference on Mechatronics and Automation．Piscatawan：IEEE，2006：2074-2078．

［11］ GUO Xiaoxia，YANG Huizhong．Wavelet Packet De-noising Based on Several Thresholds［C］∥Proceedings of the 27thChinese Control Conference．Beijing：Beihang University Press，2008：169-172．（郭曉霞，楊慧中．基于多閾值的小波包去噪［C］∥第27屆中國控制會議．北京：北京航空航天大學(xué)出版社，2008：169-172．）

［12］ TAN Wencai，ZHANG Qiuju．An Improved De-noising Method by Wavelet Packet Multi-threshold［J］．Journal of Jiangnan University：Natural Science Edition，2012，11（2）：178-181．（譚文才，張秋菊．小波包多閾值去噪的一種方法改進(jìn)［J］．江南大學(xué)學(xué)報：自然科學(xué)版，2012，11（2）：178-181．）

［13］ WU Zhaocai，LIU Tianyou．Wavelet Transform Methods in Seismic Data Noise Attenuation［J］．Progress in Geophysics，2008，23（2）：493-499．（吳招才，劉天佑．地震數(shù)據(jù)去噪中的小波方法［J］．地球物理學(xué)進(jìn)展，2008，23（2）：493-499．）

［14］ OGAJA C，RIZOS C，WANG Jingling，et al．A Dynamic GPS System for On-line Structural Monitoring［C］∥International Symposium on Kinematic Systems in Geodesy，Geomatics＆Navigation（KIS 2001）．Banff：［s．n．］，2001：5-8．

［15］ JI Yuebo．Frequency Order of Wavelet Packet［J］．Journal of Vibration and Shock，2005，24（3）：97-99．（紀(jì)躍波．小波包的頻率順序［J］．振動與沖擊，2005，24（3）：96-98．）

［16］ GE Zhexue，SHA wei．Wavelet Analysis Theory and MATLAB R2007 Implementation［M］．Beijing：Publishing House of Electronics Industry，2007：340-342．（葛哲學(xué)，沙威．小波分析理論與MATLAB R2007實現(xiàn)［M］．北京：電子工業(yè)出版社，2007：340-342．）

［17］ QU Guoqing，DANG Yamin，ZHANG Chuanyin，et al．Analysis and Improvement of De-noising Method for Wavelet Packet［J］．Journal of Geodesy and Geodynamics，2008，28（4）：102-110．（曲國慶，黨亞民，章傳銀，等．小波包消噪方法分析及改進(jìn)［J］．大地測量與地球動力學(xué)，2008，28（4）：102-110．）

［18］ ZHANG Defeng．Matlab Wavelet Analysis［M］．Beijing：China Machine Press，2009：93-97．（張德豐．Matlab小波分析［M］．北京：機(jī)械工業(yè)出版社，2009：93-97．）

［19］ VIRMANI J，KUMAR V，KALRA N，et al．SVM-Based Characterization of Liver Ultrasound Images Using Wavelet Packet Texture Descriptors［J］．Journal of Digital Imaging，2013：1-14．

［20］ TAO Ke，ZHU Jianjun．A Hybrid Indicator for Determining the Best Decomposition Scale of Wavelet Denoising［J］．Acta Geodaetica et Cartographica Sinica，2012，41（5）：749-755．（陶珂，朱建軍．多指標(biāo)融合的小波去噪最佳分解尺度選擇方法［J］．測繪學(xué)報，2012，41（5）：749-755．）

［21］ WU Fumei，YANG Yuanxi．GPS／INS Integrated Navigation by Adaptive Filtering Based on Wavelet Threshold De-noising［J］．Acta Geodaetica et Cartographica Sinica，2007，36（2）：124-128．（吳富梅，楊元喜．基于小波閾值消噪自適應(yīng)濾波的GPS／INS組合導(dǎo)航［J］．測繪學(xué)報，2007，36（2）：124-128．）

（責(zé)任編輯：宋啟凡）

Multi-threshold Wavelet Packet De-noising Method and lts Application in Deformation Analysis

ZHANG Zhetao，ZHU Jianjun，KUANG Cuilin，ZHOU Cui
School of Geosciences and lnfo-Physics，Central South University，Changsha 410083，China

ln the field of deformation monitoring，the traditional wavelet de-noising method retains only the low frequency of useful information．lt is easy to get rid of intermediate frequency and high frequency useful information．The wavelet packet analysis is a new kind of wavelet analysis method developed in recent years，which is a more subtle de-noising method for considering useful information of the various bands．The key of the wavelet packet de-noising is to select the appropriate threshold criteria and to process the wavelet packet decomposition coefficients by the threshold，but the researches using traditional wavelet packet de-noising method are not sufficient．This article is for the lack of traditional wavelet and wavelet packet analysis．According to the distribution of different signals and their noise，wavelet packet decomposition coefficients are arranged by the frequency order，and segmented in accordance with information type，to select the appropriate threshold criteria for each band and to perform threshold processing．lt is the method of wavelet packet de-noising with multi-threshold criteria based on frequency order．The results show that this method can effectively remove the noise of each band through theoretical analysis and practical applications．The de-noising ability of this method is better than the other methods such as traditional wavelet de-noising or wavelet packet de-noising．Studies have shown that this method can preserve the useful information from the de-noising signal after de-noising when the sampling frequency is low．Therefore，it can be widely used in the field of high-precision deformation monitoring．

wavelet packet analysis；deformation monitoring；data de-noising；frequency order；multithreshold criteria

ZHANG Zhetao（1988—），male，postgraduate，majors in modern surveying data processing．

P207

A

1001-1595（2014）01-0013-08

國家973計劃子課題（2013CB733303）；國家863計劃（2012AA121301）；國家自然科學(xué)基金（41274010；40974007）；中南大學(xué)中央高?；究蒲袠I(yè)務(wù)費專項資金（2013zzts254）

2012-12-27

章浙濤（1988—），男，碩士生，主要從事現(xiàn)代測量數(shù)據(jù)處理理論方面的研究。

E-mail：zhetaozhang＠y eah．net

ZHANG Zhetao，ZHU Jianjun，KUANG Cuilin，et al．Multi-threshold Wavelet Packet De-noising Method and Its Application in Deformation Analysis［J］．Acta Geodaetica et Cartographica Sinica，2014，43（1）：13-20．（章浙濤，朱建軍，匡翠林，等小波包多閾值去噪法及其在形變分析中的應(yīng)用［J］．測繪學(xué)報，2014，43（1）：13-20．）

10．13485／j．cnki．11-2089．2014．0003

修回日期：2013-06-25