李戰(zhàn)明 楊守祥
基于提升小波的心電信號基線漂移的去除方法*
李戰(zhàn)明①楊守祥①
目的:避免傳統(tǒng)小波變換基于卷積算法中的冗余計算,同時去除心電信號(ECG)在采集中混入其中的基線漂移噪聲。方法:根據(jù)提升小波變換采取雙小波基函數(shù)結(jié)合的方法,經(jīng)分解、含噪聲子帶系數(shù)置零、逆變換形成去噪的心電信號。結(jié)果:運用MATLAB環(huán)境對MIT-BIH數(shù)據(jù)庫提供的心電信號數(shù)據(jù)及基線漂移噪聲信號bw進行去除基線漂移仿真驗證,其基線漂移均被有效去除。結(jié)論:ECG信號經(jīng)該方法處理后其所含有的基線漂移噪聲被準確去除,且原信號中的波形信息被有效保留,可為心電信號特征參數(shù)的檢測提供幫助。
心電信號;基線漂移;提升小波;特征參數(shù)
[First-author’s address]Electrical Engineering and Information Engineering College, Lanzhou University of Technology, Lanzhou 730050, China.
常規(guī)心電信號(electrocardiogram,ECG)是毫伏(mV)級信號,頻帶范圍為0.05~100 Hz[1]。絕大部分的心電信號能量集中在0.25~40 Hz之間,屬于低幅、低頻信號。在心電信號所處的低信噪比的干擾環(huán)境下,其作為微弱信號在采集過程中多種噪聲信號容易被引入。主要的干擾噪聲包括:基線漂移、工頻干擾、肌電干擾等3種[2]。干擾噪聲致使心電信號的準確檢測變得很困難,因此必須引入噪聲去除的預(yù)處理,其中首要的一項是去除基線漂移?;€漂移是由人體的呼吸或者肢體運動所產(chǎn)生的一種上下波動并且緩慢變化的低頻信號,致使心電信號原有的基線位置發(fā)生變化,造成波形的失真。去除基線漂移的方法有很多,目前較常用的方法是利用傳統(tǒng)的小波變換,采用多分辨率分析的方法將基線漂移去除,但是該方法在運行時計算量大,處理時間長,不利于實時應(yīng)用。
提升小波變換與傳統(tǒng)的小波變換相比繼承了傳統(tǒng)小波變換的多分辨率特性,同時避免了傳統(tǒng)的基于卷積算法中的冗余計算[3]。提升小波變換由于其計算速度快,可進行基于原位存儲的計算,因此占用內(nèi)存少,減少了存儲器的開銷,易于硬件的實現(xiàn)。張德平等[4]和翁羽潔等[5]分別指出了用提升小波變換替代傳統(tǒng)的小波變換,計算所需時間約為傳統(tǒng)小波變換的1/3。
基于提升小波變換相比傳統(tǒng)小波變換的一系列優(yōu)勢,本實驗提出了一種在運用提升小波變換的基礎(chǔ)上采用雙小波基函數(shù)結(jié)合的方法,完成對心電信號中基線漂移噪聲的去除,經(jīng)MATLAB對MIT-BIH數(shù)據(jù)庫中的心電信號數(shù)據(jù)進行仿真實驗,證明該算法在去噪精度上要高于基于傳統(tǒng)小波變換的方法。
傳統(tǒng)小波變換作為從傅里葉分析的基礎(chǔ)上發(fā)展而來的一種分析方法,在一定程度上受到其影響。小波分析和多分辨率分析作為小波變換中的兩個核心概念均是建立在二進平移和伸縮思想之上的,這種經(jīng)典的多分辨率分析框架構(gòu)成的小波被稱為第一代小波。1995年Sweldens[6]提出了一種不依賴于傅里葉變換的新的小波構(gòu)造方法—提升格式(Lifting Sweldens),稱之為第二代小波變換。提升格式既保持了第一代小波的特性,同時又克服了其平移和伸縮的不變性,復(fù)雜度只有原來卷積方法的一半左右[7]。提升方案將第一代小波變換過程分為3個階段:分裂(Split)、預(yù)測(Predict)和更新(Update)[8]。其分解與重構(gòu)過程如圖1所示。
圖1 提升小波變換的分解與重構(gòu)
(1)分裂(Split)。將原始信號sj中位于偶下標位置的元素組成一個集合,記為evenj-1;將位于奇下標位置的元素組成一個集合,記為oddj-1,其過程為公式1,懶小波變換即是將信號分解成為奇偶子集序列的過程。
(2)預(yù)測(Predict)。一個局部關(guān)聯(lián)性較強信號的偶子集序列和奇子集序列是密切相關(guān)的。因此,根據(jù)其中的任何一個就能夠使用其在合理的精度范圍內(nèi)去預(yù)測另一個。通常用偶子集evenj-1的預(yù)測值P(evenj-1)去預(yù)測奇子集oddj-1,也即將預(yù)測算子P對偶子集作用以后當(dāng)作奇子集的預(yù)測值使用,將奇子集的實際值與預(yù)測值相減可得到細節(jié)信號dj-1,預(yù)測過程為公式2。
(3)更新(Update)。為了使原始信號的子集能夠繼續(xù)體現(xiàn)其某些全局特性,利用上一步預(yù)測得到的奇子集序列,即細節(jié)信號dj-1對偶子集evenj-1進行更新,更新算子為U,更新的意圖就是要尋找一個更好的子集sj-1,即低頻信號,使其能夠保持原始信號的某一標量特性。更新過程為公式3。
繼續(xù)的分解過程是將sj-1重復(fù)上述公式(1)、(2)、(3)的步驟,這樣即可得到下一層的小波分解系數(shù)。提升格式的小波分解的每一步都是可逆的,其重構(gòu)過程就是分解過程的逆變換,只需改變正變換過程中的加減符號即可得到其逆變換。逆變換的算法為公式4。
式中Merge為合并算子。
2.1 分解尺度的確定
基線漂移的頻率<1 Hz,而心電信號中也同樣含有頻率較低的有用信號,其中ST段的頻率范圍為0.7~2.0 Hz,P、T波頻率范圍為0.7~10 Hz[9-11]。因此,如何在保證有用信號不受損害的同時將基線漂移去除掉是去噪時所考慮的關(guān)鍵。屬低頻段的該噪聲在小波分解的過程中主要分布在高尺度子帶上,在MIT-BIH數(shù)據(jù)庫中所提供的原始信號采樣頻率360 Hz保持不變的情況下,在假定進行8層小波分解時,各尺度上子帶系數(shù)的頻率范圍見表1。
表1 各尺度上子帶系數(shù)的頻率范圍
結(jié)合表1,并根據(jù)基線漂移噪聲及心電信號中低頻波段的頻率范圍,可將基線漂移的頻率范圍認定為0~0.7 Hz,與Jeyarani等[12]和Doudou等[13]認定的范圍相同,即是ca8子帶所覆蓋的頻率范圍,因此將分解尺度確定為j=8即可。
2.2 雙小波基函數(shù)的選取
利用傳統(tǒng)小波進行基線漂移去噪時,對小波基函數(shù)的選取并沒有特殊的要求,但在實驗中發(fā)現(xiàn)如果將提升格式的小波分解直接移植到傳統(tǒng)的小波分解中時,采用單一的小波基函數(shù)并不能將基線漂移去除干凈,有的小波基函數(shù)甚至?xí)π碾妶D波形造成破壞。如采用db2、db4、db5和sym4小波基函數(shù)對含噪信號進行處理時發(fā)現(xiàn),信號尾端中的基線漂移不能去除掉;而如果采用bior4.4小波基函數(shù),經(jīng)處理后所得信號的首段中的基線漂移不能去除掉;而如果采用coif2、db6和db8等數(shù)種在傳統(tǒng)小波變換中常用的小波基函數(shù),不僅不能達到去噪目的,還會破壞原有的心電圖波形?;谏鲜龅某醪綄嶒灲Y(jié)果,為了達到去噪的目的,本文選用了雙小波基函數(shù)結(jié)合的方法,即在選用bior4.4的基礎(chǔ)上,再附加使用db2、db4、db5或sym4小波基函數(shù)。
2.3 去噪過程
整個去噪過程分為5個步驟。
(1)含噪信號xn的讀取。
(2)運用提升小波變換采用bior4.4對整段信號進行8層分解,采用db2、db4、db5或sym4對整段信號中的首2/3段的信號進行8層分解。如圖2所示,即是采用bior4.4對整段信號進行8層分解所得到的分解結(jié)果。
圖2 提升小波(bior4.4)8層分解結(jié)果
(3)對兩次分解中所得到的第8尺度上的近似系數(shù)進行置零處理。
(4)將置零處理后的第8尺度上的近似系數(shù)與1~8尺度上的細節(jié)系數(shù)逐級重構(gòu),得到經(jīng)bior4.4分解的與原信號等長的信號xn1和經(jīng)db2、db4、db5或sym4分解的長度是原信號2/3的信號xn2。
(5)對步驟(4)中所得到的兩段信號在整段信號的0.45處作結(jié)合處理,在0.45處及之前使用xn2中的數(shù)據(jù),在0.45處之后使用xn1中的數(shù)據(jù),在還原與原信號等長心電信號的同時,去除基線漂移噪聲。
由圖2可知,當(dāng)運用提升小波采用bior4.4進行8層分解時,隨著尺度的增加,高頻信號逐級減少而低頻信號逐漸增加,當(dāng)?shù)竭_第8尺度時,子帶ca8幾乎囊括了大部分的基線漂移的信息,因此經(jīng)過步驟(1)至步驟(5)的運算即可達到去噪目的。
3.1 數(shù)據(jù)來源
由于心電信號中波形形態(tài)的多樣性,單一時段的實驗分析不足以充分證明一種算法的有效性,因此本研究在實驗中選用來自于MIT-BIH數(shù)據(jù)庫中的123號數(shù)據(jù)中基本無噪聲的6個不同時段,時段信號長度N=2048,并為了建立含噪模型,在其上疊加了頻率為0.15 Hz的正弦波。為了檢驗本實驗中算法對實際心電信號中基線漂移的去除效果,選取了該數(shù)據(jù)庫中所專門提供的基線漂移的噪聲信號bw與一段含真實基線漂移的心電信號進行實驗對比。
3.2 去噪效果對比分析
實驗采用信噪比SNR作為信號的最終去噪評價標準,其定義為公式5。
式中x(i)為無噪聲的心電信號,y(i)為去噪后的心電信號。
為了確定最佳雙小波基函數(shù)組合并與傳統(tǒng)的小波變換方法進行去噪效果對比,本實驗采用bior4.4+db2、db4、db5或sym4等4種不同的雙小波基函數(shù)搭配,傳統(tǒng)小波變換方法采用bior3.7、db3、db5和sym4等4種文獻[14-17]中在去除基線漂移時所使用的小波基函數(shù)。利用上述兩種方案分別對初次選取的6個不同時段的數(shù)據(jù)進行去噪處理,去噪后的信噪比對比見表2。
表2 兩種去噪方案的實驗信噪比(dB)對比
表2中的信噪比對比結(jié)果顯示,提升小波變換方法完成去噪時平均信噪比較高的小波基函數(shù)組合有bior4.4+db4與bior4.4+db5,而傳統(tǒng)小波變換方法在采用不同小波基函數(shù)完成去噪時的平均信噪比相差不大,約為18.8 dB。就兩種方法中不同小波基函數(shù)的選用而言,提升小波變換采用bior4.4+db4雙小波基函數(shù)組合的方案無論在整體還是在大部分的分段數(shù)據(jù)去噪中都取得最佳的效果,因此最佳的雙小波基函數(shù)的組合即是bior4.4+db4。
利用提升小波采用最佳雙小波基函數(shù)組合進一步對bw中的一段信號及123號數(shù)據(jù)中實際含基線漂移較明顯的一段信號進行去除基線漂移處理,結(jié)果分別如圖3、圖4所示。在該去噪過程中,完成去除所選用的心電信號中的基線漂移時的信噪比SNR=13.5694 dB,高于傳統(tǒng)小波變換在采用sym4小波基函數(shù)完成相同去噪時的信噪比,其信噪比SNR=13.2454 dB。
如圖3所示,通過對比可知去除基線漂移后的信號的基線回到了零點位置,信號中的其他信息得到了很好地保留,基線漂移被有效地分離掉。
圖3 bw的去噪前后對比
如圖4所示,通過對比可知原始信號中含有較明顯的基線漂移噪聲,經(jīng)去噪處理后基線同樣回到零點位置,心電信號中的主要波形信息未被破壞。
基于上述的實驗結(jié)果與分析可知,提升小波變換在采用合理的雙小波基函數(shù)組合的基礎(chǔ)上是能夠?qū)⑿碾娦盘栔械幕€漂移噪聲進行有效去除的,并且去噪效果令人滿意。
圖4 含基線漂移的心電信號去噪前后對比
提升小波變換因其計算量小,易于硬件實現(xiàn)的特點,使得其在心電信號實時處理中有著良好的應(yīng)用前景。本實驗采用提升小波變換利用bior4.4+db4雙小波基函數(shù)結(jié)合的方法完成了心電信號中基線漂移的去除,實驗結(jié)果表明去噪效果理想,較好地保留了原始心電信號中的波形信息,并且噪聲干擾去除后心電信號具有較高的信噪比。
[1]Saritha C,Sukanya V,Murthy YN.ECG signal analysis using wavelet transforms[J].Department of Physics and Electronics,2008,35:68-77.
[2]張徑周,壽國法,戴冠中.基于小波變換的心電信號噪聲處理[J].西北工業(yè)大學(xué)報,2005,23(1):11-14.
[3]葛哲學(xué),沙威.小波分析理論與MATLAB R2007實現(xiàn)[M].北京:電子工業(yè)出版社,2007.
[4]張德平,賈文娜,楊葉青.基于提升方案的心電信號去噪算法[J].計算機仿真,2010,27(10):226-229.
[5]翁羽潔,丁勇,孫立艷,等.一種基于提升小波和中值濾波的心電去噪方法[J].北京生物醫(yī)學(xué)工程,2010,29(5):465-469.
[6]Sweldens W.The lifting scheme:a construction of second generation of wavelets[J].SIAM J Math Anal,1998,29(2):511-546.
[7]張德豐.MATLAB小波分析[M].北京:機械工業(yè)出版社,2009.
[8]孫延奎.小波分析及其應(yīng)用[M].北京:機械工業(yè)出版社,2005.
[9]Mbachu CB,Victor I,Emmanuel I,et al.Filtration of artifacts in ECG signal using rectangular window-based digital filters[J].International Journal of Computer Science Issues,2011,8(5):279-285.
[10]張華,潘冬明,王立會.小波包變換在消除心電圖基線漂移方面的研究[J].中國醫(yī)學(xué)物理學(xué)雜志,2007,24(1):71-73.
[11]李肅義.可穿戴生理參數(shù)監(jiān)測系統(tǒng)的動態(tài)心電信號處理方法研究[D].長春:吉林大學(xué)儀器科學(xué)與電氣工程學(xué)院,2009.
[12]Jeyarani AD,Singh TJ.Analysis of noise reduction techniques on QRS ECG waveformby Applying different filters[J].IEEE Recent Advances in Space Technology Services and Climate Change (RSTSCC),2010,10:149-152.
[13]Doudou Wu,Zhengyao Bai.An improved method for ECG signal feature point detection based on wavelet transform[J].IEEE Conference on Industrial Electronics and Applications (ICIEA),2012,7:1836-1841.
[14]倪原,王曉麗,陸文總.心電信號的小波變換處理算法及仿真[J].西安工業(yè)大學(xué)學(xué)報,2012,32(4):310-314.
[15]任杰,楊麗曉.基于小波變換系數(shù)的心電信號基線漂移噪聲去除方法[J].醫(yī)療衛(wèi)生裝備,2010,31(11):24-26.
[16]張文瓊,劉肖琳,吳濤.一種利用小波變換逼近信號濾除心電圖基線漂移的方法[J].計算機工程與應(yīng)用,2005,41(20):222-224.
[17]張德平.基于提升方案的心電信號去噪與檢測算法研究[D].上海:上海交通大學(xué)電子信息與電氣工程學(xué)院,2010.
Methods of removing the baseline wander in ECG based on the lifting wavelet transform
LI Zhan-ming, YANG Shou-xiang// China Medical Equipment,2014,11(3):16-19.
Objective:To avoid the redundant computation based on the convolution operation in the traditional wavelet transform, and to remove the baseline wander noise existing in the course of collecting the ECG signal.Methods:Use the lifting wavelet transform with two wavelets, and constitute the ECG signal with the noise removed after decomposing, setting the subband coefficient including the noise to zero, and rebuilding.Results:Use MATLAB to remove the baseline wander noise in the ECG signal and bw provided by the MIT-BIH database, and the results show that the baseline wander was removed effectively.Conclusion:The baseline wander noise in the ECG signal can be removed accurately though the method mentioned above, the waveform information in the original ECG signal can be maintained effectively, and subsequently, that can provide help for detecting the characteristic parameters in the ECG signal.
Electrocardiogram; Baseline wander; Lifting wavelet; Characteristic parameters
10.3969/J.ISSN.1672-8270.2014.03.006
1672-8270(2014)03-0016-04
R318.04
A
2013-11-07
教育部博士點基金(20106201110003)“壓縮域檢索與挖掘關(guān)鍵技術(shù)研究”
①蘭州理工大學(xué)電氣工程與信息工程學(xué)院 甘肅 蘭州 730050
李戰(zhàn)明,男,(1962- ),博士,博士生導(dǎo)師,教授。蘭州理工大學(xué)電氣工程與信息工程學(xué)院院長,從事控制理論與控制工程學(xué)科的教學(xué)與科研工作。