涂斌斌， 谷麗華， 許 會

(1. 沈陽工業(yè)大學 信息科學與工程學院， 沈陽 110870； 2. 沈陽大學 信息工程學院， 沈陽 110044)

步態(tài)識別的小波去噪質量評價方法*

涂斌斌1，2， 谷麗華1， 許 會1

(1. 沈陽工業(yè)大學 信息科學與工程學院， 沈陽 110870； 2. 沈陽大學 信息工程學院， 沈陽 110044)

針對步態(tài)加速度信號由于采集信號包含大量尖峰或突變等高頻噪聲，導致在步態(tài)分析、周期劃分和特征提取方面存在困難的問題，提出了一種利用小波去噪質量評價來對加速度信號小波去噪的最優(yōu)參數(shù)進行選擇的新方法.根據(jù)質量評價指標的幾何和物理意義，使用評價指標值來指導小波去噪最優(yōu)小波分解與重構層數(shù)的參數(shù)選擇.當評價指標值取最小值時，其對應分解層數(shù)為最優(yōu)分解層數(shù)，采用的小波基函數(shù)為最優(yōu)小波基函數(shù).結果表明，本文方法在步態(tài)識別領域可以為步態(tài)信號小波去噪選擇最優(yōu)小波基函數(shù)和分解層數(shù)提供理論依據(jù).

步態(tài)； 步態(tài)識別； 加速度信號； 小波去噪； 質量評價； 小波基函數(shù)； 分解與重構； 最優(yōu)參數(shù)

近些年，隨著便攜式及可穿戴式智能設備的快速發(fā)展，基于加速度傳感器的步態(tài)生物特征識別技術成為研究熱點[1-3].步態(tài)作為生物特征之一，具有普遍性、可采集性和唯一性，且不易偽造及仿冒，可在無需主動配合的情況下，隨時隨地保護便攜式電子產(chǎn)品的信息安全[4].使用加速度傳感器采集人體步態(tài)信號時噪聲污染難以避免，因此，從混有噪聲的步態(tài)信號中提取有效信息成為研究工作重點.

在步態(tài)識別預處理過程中，信號的去噪方法主要包括：FIR低通濾波、移動均值濾波、加權移動均值濾波及小波去噪等[5-6].由于小波分析具備良好的時域局部性和多分辨率分析能力，借助小波分析與重構技術，小波去噪已經(jīng)在步態(tài)識別預處理中得到了廣泛應用[7].但是采集信號特征的未知性直接導致小波去噪質量的復雜性，很難正確選取去噪所需的最優(yōu)參數(shù).小波去噪原理是基于非統(tǒng)計理論，其特點源于小波分析技術不考慮系統(tǒng)的概率分布，允許系統(tǒng)的概率分布是未知的或非典型的.為此，本文將系統(tǒng)分析現(xiàn)有小波去噪方法在步態(tài)識別預處理中的質量效果，通過完善小波去噪最優(yōu)參數(shù)選擇標準來提高步態(tài)識別預處理的有效性.

1 步態(tài)加速度信號分析

1.1 步態(tài)信號的特點

加速度傳感器采集人體正常行走時三維方向的加速度信號.研究顯示，正?；顒又腥梭w軀干所產(chǎn)生加速度值在6g以下，范圍集中在[-2g，2g][8].一般情況，重力方向的加速度信號幅度為-0.3g～0.8g，身體側向腰部加速度幅度為-0.3g～0.4g.行走時，人體活動產(chǎn)生低頻信號，頻率在20 Hz以下，主要集中在1～5 Hz范圍內.

1.2 信號噪聲的特點

加速度傳感器采集人體步態(tài)信號時，易受多種噪聲源影響，包括身體晃動、環(huán)境噪聲、呼吸狀況、工頻和高頻干擾等.這些產(chǎn)生的噪聲混雜在采集信號中，直接導致有效步態(tài)分析、特征提取的困難，主要表現(xiàn)為步態(tài)信號含有非正常的高頻成分.加速度傳感器步態(tài)采集裝置雖然通過硬件濾波能去除部分噪聲，但高頻信號依然存在，本文以小波去噪軟件濾波的方法來去除噪聲影響.

2 小波去噪質量評價方法

2.1 小波去噪基本原理

小波變換實現(xiàn)信號去噪原理：使用小波基函數(shù)去逼近元函數(shù)，從而完成對元函數(shù)的分解.即利用小波變換的帶通功能，將原信號分解為互不重疊頻率的信號，再擇優(yōu)選取閾值來有效地去除信號中混入的噪聲，小波去噪的基本步驟可概括如下：

1) 擇優(yōu)選擇小波基函數(shù)及分解層數(shù)，計算各層小波分解系數(shù)；

2) 通過對高頻系數(shù)門限閾值量化處理來消除高頻部分的噪聲；

3) 對低頻系數(shù)和去噪后的高頻系數(shù)進行小波重構，得到去噪的真實信號.

2.2 小波去噪質量評價

傳統(tǒng)小波去噪質量評價使用的4個指標為：均方根誤差(RMSE)、信噪比(SNR)、互相關系數(shù)(R)以及平滑度(r)[9]，具體公式和特點如下.

1) 均方根誤差表示為

(1)

2) 信噪比表示為

(2)

式中，psiganl、pnoise為原始信號的功率和噪聲功率.信噪比是指去噪前后信號間功率的比值，其值越大表示去噪效果越好.

3) 互相關系數(shù)表示為

(3)

4) 平滑度表示為

(4)

平滑度是指去噪前后信號一階差分之間方差根的比值，其值越小表示去噪效果越好.

傳統(tǒng)評價指標雖然擁有具體物理意義，但也具有一定局限性.當去噪效果甚微或未達到去噪目的時，將出現(xiàn)均方根誤差值接近0，相關系數(shù)接近1，而信噪比值極大的情況，顯然這種評價結果是無效的.本文利用文獻[10]所提出的復合評價指標來衡量步態(tài)信號小波去噪效果.

復合評價指標公式為

T=WPRMSEPRMSE+WPrPr

(5)

(6)

(7)

式中：PRMSE為標準化的均方根誤差(RMSE)；Pr為標準化的平滑度；σPRMSE為PRMSE的標準差；μPRMSE為PRMSE的均值；σPr為Pr的標準差；μPr為Pr的均值.復合評價指標由均方根誤差、平滑度及其變異系數(shù)定權決定，具有保留細節(jié)信息和逼近信息的物理意義.當復合評價指標值取最小值時，其對應分解層數(shù)為最優(yōu)分解層數(shù)，其采用的小波基函數(shù)可作為最優(yōu)小波基函數(shù)來備選.

3 步態(tài)信號小波去噪

對步態(tài)信號進行分析時，發(fā)現(xiàn)信號包含大量尖峰或突變等噪聲，導致產(chǎn)生極值點、步態(tài)周期劃分和特征提取不一致現(xiàn)象.采用小波變換的方法去噪，能有效去除噪聲，最大限度地提取有用步態(tài)信息.在利用小波變換方法進行去噪的過程中，小波基函數(shù)和分解層數(shù)的選擇十分重要.小波基函數(shù)和分解層數(shù)的差異直接影響真實信號特征提取的差異，Daubechies系列小波基是典型的正交小波基，具有良好的正則性，使得信號重構過程比較光滑.

小波由一組多尺度濾波器迭代而成，小波去噪對消除非平穩(wěn)信號的噪聲具有無可比擬的優(yōu)點.文獻[7]首次將小波去噪方法應用到步態(tài)加速度信號濾波中，以消除因外界環(huán)境或身體晃動等引起的噪聲；文獻[11]使用小波基函數(shù)db13對垂直方向加速度信號進行6層尺度小波分解和重構；文獻[12]使用小波基函數(shù)db8對加速度信號進行2層尺度小波分解和重構；文獻[13]使用小波基函數(shù)db5對加速度信號進行3層尺度小波分解和重構.以上3個文獻都提出了各自使用的小波基函數(shù)和分解尺度層數(shù)，但未給出實際的理論依據(jù)，僅僅對比2種小波基函數(shù)分解和重構后的濾波效果，并不能解釋其為最優(yōu)小波基函數(shù)和最優(yōu)分解層次的根本原因.

4 實驗設計與分析

4.1 實驗數(shù)據(jù)

實驗數(shù)據(jù)來源于第一個公開的步態(tài)加速度數(shù)據(jù)集(ZJU-GaitAcc數(shù)據(jù)集[14])，包括身體5個位置(右手腕、左上臂、髖部、左大腿和右腳踝)的步態(tài)加速度數(shù)據(jù)，所使用的加速度傳感器為ADXL330三軸加速度傳感器，采樣率為100Hz.數(shù)據(jù)集中包含3個子集，每個子集數(shù)據(jù)采集于5個身體位置6次自然行走20m、水平地面的正常行走步態(tài)信號.本文研究僅針對子集1中右髖部位置22名測試者的6次采樣步態(tài)加速度數(shù)據(jù)進行實驗，其加速度信號如圖1所示.

圖1 右髖部位置步態(tài)加速度信號Fig.1 Gait acceleration signals in position of right pelvis

4.2 實驗設計

步態(tài)加速度信號小波去噪質量的復合評價指標實驗設計如下：

1) 確定備選參數(shù)，分解層數(shù)為2～7層，小波基函數(shù)為db5～db13小波，閾值選定為heursure閾值準則及軟閾值處理函數(shù)；

2) 計算均方根誤差及平滑度的值，并歸一化，計算WPRMSE與WPr權重系數(shù)的值；

3) 計算復合評價指標，當其取最小值時，所對應的小波基函數(shù)和分解層數(shù)為最優(yōu)備選參數(shù)；

4) 分析備選參數(shù)小波重構后的濾波效果.

4.3 實驗結果分析

1) 1名測試者1次采集步態(tài)信號的9個小波基函數(shù)與6個分解層數(shù)復合評價指標如表1所示.

表1 多層次多小波基函數(shù)的評價指標Tab.1 Evaluation index of multi-level wavelet basis function

由表1可知，無論采用何種小波基函數(shù)對步態(tài)信號進行分解，當分解層為4層時，得到的復合評價指標值都為最小值，復合指標值小于0.2；當分解層為2層時，得到的復合評價指標值都為最大值，復合指標值大于0.6.復合指標值從2層到3層間遞減幅度較大，3～4層遞減幅度趨于平穩(wěn)，4層時為最小值，4～5層遞增幅度較大，5～7層時增加幅度趨于平穩(wěn)，其中3～7層的復合指標值基本處于[0.1，0.3]范圍內.從數(shù)值大小分析，db8、db11、db12為復合評價指標極小的三種小波基函數(shù)，可以作為最優(yōu)備選參數(shù).

實際上，人正常步行時的步頻包含在0.625～20 Hz頻帶內，其中1～5 Hz頻率范圍內集中了最多正常行走的頻率信息.研究表明，99%以上的運動能量集中在15 Hz以下的信號成分中，正常行走的步頻在1～3 Hz(慢速為1～2 Hz，常速為2～2.5 Hz，快速為2.5～3 Hz)[15].圖2為4層分解中濾波效果對比和頻譜對比圖.從圖2可知，以db5小波基函數(shù)對步態(tài)信號進行4層分解和重構的去噪后，在5～20 Hz頻率范圍內，依然含有頻率信息；以db8小波基函數(shù)對步態(tài)信號進行4層分解和重構的去噪后，步態(tài)信號主要步頻信息集中在1～5 Hz頻率范圍內，滿足對步態(tài)信號研究的濾波需求.圖3為db8濾波效果對比和頻譜對比圖，從圖3可知，使用相同小波基函數(shù)db8對步態(tài)信號進行去噪，分解層數(shù)為3時，在5～10 Hz頻率范圍內含有較多頻率信息，對應濾波后的步態(tài)信號仍然含有毛刺.

圖2 4層分解中濾波效果對比和頻譜對比Fig.2 Comparison in filtering effect and frequency spectrum for 4-level decomposition

圖3 db8濾波效果對比和頻譜對比Fig.3 Comparison in filtering effect and frequency spectrum for db8

2) 1名測試者6次采集步態(tài)信號的9種小波基函數(shù)與6個分解層數(shù)復合評價指標圖如圖4所示.由圖4可見，對同1名測試者進行6次采樣，每一次采樣用9種小波基函數(shù)進行6個層數(shù)分解和重構，所計算的復合評價指標規(guī)律與上述分析結果一致.當分解層為4層時，得到的復合評價指標值都為最小值，且其值小于0.2.最優(yōu)小波基函數(shù)選取方法如表2所示.

圖4 復合評價趨勢圖(1人)Fig.4 Trend graphs of composite evaluation (1 person)

表2 最優(yōu)小波基函數(shù)選取Tab.2 Selection method for optimal wavelet basis function

對比6次評價結果可知，其最小值所對應的小波基函數(shù)無規(guī)律可尋，6次評價結果具有明顯不同，可見每次評價結果只適于每一次采樣數(shù)據(jù).

3) 多名測試者多次采集步態(tài)信號的多小波基函數(shù)、多分解層數(shù)復合評價指標分析.對22名測試者的6次步態(tài)信號進行132×9次復合評價指標分析，評價結果與上述研究基本一致.其中，5～7層層數(shù)分解時，10號、15號、17號測試者有部分或者全部復合指標值處于[0.3，0.4]范圍內；以3層層數(shù)分解時，10號測試者有部分復合指標值大于0.3.圖5為10名測試者采用不同基波函數(shù)的復合評價指標趨勢圖.

5 結 論

圖5 復合評價趨勢圖(10人)Fig.5 Trend graphs of composite evaluation (10 persons)

本文將小波去噪質量評價指標應用到步態(tài)信號預處理過程中，通過評價指標值來指導小波去噪最優(yōu)備選參數(shù)的選取.采用質量評價指標對步態(tài)加速度信號進行分析，得到了滿意的效果，給出如下結論：

1) 以db5～db13小波基函數(shù)進行2層尺度小波分解和小波重構的質量評價指標值最大，根據(jù)判定規(guī)則，不適合用于步態(tài)信號小波去噪處理；

2) 分解層數(shù)為4層時，復合評價指標達到最小值，根據(jù)判定規(guī)則，是最適合的小波分解與重構層數(shù)；

3) 同1名測試者多次采樣的步態(tài)信號進行4層小波分解和小波重構的小波基函數(shù)并不唯一.

[1]Khandelwal S，Wickstr?m N.Gait event detection in real-world environment for long-term applications：incorporating domain knowledge into time-frequency analysis [J].IEEE Transactions on Neural System and Rehabilitation Engineering，2016，99：1-10.

[2]Derlatka M.Biometric system for person recognition using gait [C]//Foundations of Intelligent Systems.Warsaw，Poland，2011：565-574.

[3]Yang C C，Hsu Y L，Shih K S，et al.Real-time gait cycle parameter recognition using a wearable accele-rometry system [J].Sensors，2011，11(8)：7314-7326.

[4]Ren Y，Chen Y，Chuah M C.User verification leveraging gait recognition for smartphone enabled mobile healthcare systems [J].IEEE Transactions on Mobile Computing，2015，14(9)：1961-1974.

[5]Soaz C，Diepold K.Step detection and parameterization for gait assessment using a single waist-worn accelerometer [J].IEEE Transactions on Biomedical Engineering，2015，63(5)：933-942.

[6]Tran K H，Chew M T.Kalman filtering for wearable fitness monitoring [C]//Proceedings of 2015 IEEE International on Instrumentation and Measurement Technology Conference (I2MTC).Ho Chi Minh City，Vietnam，2015：2020-2025.

[7]Liu R，Duan Z，Zhou J Z，et al.Identification of individual walking patterns using gait acceleration [C]//Proceedings of lst International Conference on Bioinformatics and Biomedical Engineering.Wuhan，China，2007：543-546.

[8]谷麗華，崔暢，高松巍，等.基于MPU-6050的步態(tài)信號采集系統(tǒng) [J].沈陽工業(yè)大學學報，2015，37(2)：176-182.

(GU Li-hua，CUI Chang，GAO Song-wei，et al.Gait signal acquisition system based on MPU-6050 [J].Journal of Shenyang University of Technology，2015，37(2)：176-182.)

[9]陶珂，朱建軍.小波去噪質量評價方法的對比研究 [J].大地測量與地球動力學，2012，32(2)：128-133.

(TAO Ke，ZHU Jian-jun.A comparative study on validity assessment of wavelet de-noising [J].Journal of Geodesy and Geodynamics，2012，32(2)：128-133.)

[10]朱建軍，章浙濤，匡翠林，等.一種可靠的小波去噪質量評價指標 [J].武漢大學學報(信息科學版)，2015，40(5)：688-694.

(ZHU Jian-jun，ZHANG Zhe-tao，KUANG Cui-lin，et al.A reliable evaluation indicator of wavelet de-noising [J].Journal of Wuhan University(Geomatics and Information Science)，2015，40(5)：688-694.)

[11]Li Y X，Wang X B，Qiao F.Gait authentication based on acceleration signals of ankle [J].Chinese Journal of Electronics，2011，20(3)：447-451.

[12]侯向鋒，劉蓉，周兆豐.加速度傳感器MMA7260在步態(tài)特征提取中的應用 [J].傳感技術學報，2007，20(3)：507-511.

(HOU Xiang-feng，LIU Rong，ZHOU Zhao-feng.App-lication of accelerometer MMA7260 in the gait feature extraction [J].Chinese Journal of Sensors and Actuators，2007，20(3)：507-511.)

[13]翟文婷.利用加速度傳感器實現(xiàn)人體步態(tài)識別 [D].沈陽：沈陽工業(yè)大學，2013.

(ZHAI Wen-ting.Using acceleration sensor to realize human gait recognition [D].Shenyang：Shenyang University of Technology，2013.)

[14]Zhang Y T，Pan G，Jia K，et al.Accelerometer-based gait recognition by sparse representation of signature points with clusters [J].IEEE Transactions on Cybernetics，2015，45(9)：1864-1875.

[15]Trung N T，Makihara Y，Nagahara H，et al.Performance evaluation of gait recognition using the largest inertial sensor-based gait database [C]//Proceeding of 5th IAPR International Conference on IEEE Biometrics Compendium.Osaka，Japan，2012：360-366.

(責任編輯：景 勇 英文審校：尹淑英)

Quality evaluation method for wavelet de-noising in gait recognition

TU Bin-bin1, 2, GU Li-hua1, XU Hui1

(1. School of Information Science & Engineering, Shenyang University of Technology, Shenyang 110870, China; 2. School of Information Engineering, Shenyang University, Shenyang 110044, China)

In order to overcome the difficulties in gait analysis, cycle division and feature extraction for gait acceleration singals due to the fact that the acquired signals contain a large number of higher frequency noises such as peaks or mutation, a new method, where the quality evaluation for the wavelet de-noising was adopted to select the optimal parameters for the wavelet de-noising of acceleration signals, was proposed. According to the geometric and physical meanings of quality evaluation indexes, the paramenter selection for both decomposition and reconsitution levels of optimal wavelet in the wavelet de-noising were guided by the evaluation index value. When the evaluation index took the minimum value, the correspongding decomposition level was optimal, and the used wavelet basis funtion was also optimal. The results show that the proposed method can provid the theoretical basis for the selection of both optimal wavelet basis funtion and decomposition levels for the wavelet de-noising of gait signals in the field of gait recognition.

gait; gait recognition; acceleration singal; wavelet de-noising; quality evaluation; wavelet basis function; decomposition and reconstruction; optimal parameter

2016-03-14.

國家自然科學基金資助項目(61372154).

涂斌斌(1980-)，女，遼寧沈陽人，講師，博士生，主要從事模式識別、信號處理等方面的研究.

16∶10在中國知網(wǎng)優(yōu)先數(shù)字出版.

http：∥www.cnki.net/kcms/detail/21.1189.T.20160907.1610.056.html

10.7688/j.issn.1000-1646.2017.01.12

TP 212

A

1000-1646(2017)01-0061-06