搜索
基于EMD的指端光电容积脉搏波中呼吸波提取方法研究李文彪1,陈真诚1,刘福彬2(1.,410083;2.,225003):通过人体指端的脉搏波,提取人体呼吸波信号使用多参数临床生理监护仪,同步采集人体指端光电容积脉搏波信号和胸阻抗法检测的呼吸波信号,对光电容积脉搏波信号做各层的经验模式分解,选择合适频率的本征模函数,与采集的呼吸波信号做相关性分析经验模式分解由脉搏波中所提取的呼吸波与采集的呼吸波有很好的相关性经验模式分解法可有效提取人体指端光电容积脉搏波中所包含的呼吸波成分,对改进医疗监护设备设计实现生理信号的多参数提取和精确分析有重要意义:光电容积脉搏波;经验模式分解;本征模函数;呼吸波:R318:A:10020837(2010)04027904ExtractionofRespiratoryWavefromFingerTipPhotolethysmographySignalsBasedonEMDMethodLIWenbiao,CHENZhencheng,LIUFubin.SpaceMedicine&MedicalEngineering,2010,23(4):279~282Abstract:ObjectiveToextractrespiratorywaveinformationfromphotolethysmography(PPG)signals.MethodsUsingaclinicalmultiparameterphysiologicalmonitor,thefingertipPPGsignalswereacquiredsynchronouslyandtherespiratorywavewasdetectedbythethoracicimpedancemethod.ThenempiricalmodedecompositionwasconductedforthePPGsignals,selectingtheintrinsicmodefunction(IMF)withappropriatefrequencytocarryoutcorrelationanalysistogetherwithrespiratorywavesignals.ResultsTherespiratorywavedirectlyacquiredhadagoodcorrelationwiththatderivedfromphotolethysmographysignalsbyempiricalmodedecompositionmethod.ConclusionItisindicatedthatempiricalmodedecompositionmethodcaneffectivelyextracttherespiratoryinformationcontainedinPPGsignals.Thisallowsforimprovingthedesignofmedicalmonitoringdevicesandisusefulinphysiologicalmultiparameterextractionandaccurateanalysis.Keywords:photolethysmography;empiricalmodedecomposition;intrinsicmodefunction;respiratorywaveAddressreprintrequeststo:LIWenbiao.InstituteofBiomedicalEngineering,SchoolofInfophysics&GeomaticsEngineering,CentralSouthUniversity,ChangshaHunan410083,China,,(),(photolethysmography,PPG),,,PPG[1],,,,(respiratorywave,RW),:20100319:tiger830611@163.com,,PPG,,[23],,,[4]PPG3,,PPG,,,,,[5],,,(empiricalmodedecomposition,EMD)234Vo.l23No.420108SpaceMedicine&MedicalEngineeringAug.2010Huang[6]1998,,,,EMD,,EMD,[78]EMDQRS,[910]EMD,,PPG,,PPG,11.1,,,,:1),1;2),(intrinsicmodefunctions,IMF),IMFIMF,,EMD,1.2xR1,:1)x(t)();2)3,emax(t)emin(t);3):em(t)=emax(t)+emin(t)2(1)4)x(t)em(t):h(t)=x(t)-em(t)(2)5)h(t)x(t),2)~4),,IMF:c1(t)=h1(t)(3)IMFc1(t),,x(t):x1(t)=x(t)-c1(t)(4)x(t)(2)~(4),IMFc2(t),xn(t),,x(t)IMFcj(t)(j=1,2,3n)rn(t):x(t)=!nj=1cj(t)+rn(t)(5)EMD,MATLAB,3:x=cos(10t)+cos(200t)+cos(600t)(6),t,0∀t∀600EMD,1,EMD1EMDFig.1AnexampleofEMDmethod2802322.1(2),423~26,5min,200Hz,10min2.2Huang[7],,,,,,,EMD,,1000(5)PPGIMF,SD(standarddeviation)AR(autoregressive),,IMF32.3MATLABPPG,2min,42Fig.2Dataacquisitionmodule4PPGFig.4OriginalPPGsignalsPPG,2.2,,,5S_RWD_RW,,,[11]AR,6,,,12min,5PPGFig.5RespiratorywavederivedfromthoracicimpedanceandPPGsignals2814EMD1(n=4)Table1Comparisonofmeasuredresults(n=4)subjectsdatasourcerespirationintwominutes(times)averagerespirationrate(times/min)SDvalueofrespirationratecorrelationcoefficient1S_RW2613.40.829993.6%D_RW2513.80.81212S_RW3015.30.926588.9%D_RW3615.60.93623S_RW2814.50.836691.1%D_RW3014.60.81074S_RW2311.20.790496.1%D_RW2310.60.78136ARFig.6PowerspectralanalysisbasedonARmodel,4,,96.1%,,288.9%,,PPG,;,,,SDS_RWD_RW,,3,PPG,,[][1]WANGYan,ZHANGSong,YANGYimin,etal.Developmentofanewlymonitoringmoduleforcardiovascularhemodynamicparametersbasedonthebloodvolumepulse[J].BeijingBiomedicalEngineering,2006,25(2):148150.[2]LIUBaohua,DULiantao,WANGGuoyong.Methodofremovingrespiratoryinterferencefrompulsewave[J].JournalofYanshanUniversity,2008,32(1):5357.[3]GUOJianping,HUMaoq,iZHANGQishan.Themethodofreducingrespiratoryinterferenceinrightchestimpedancewave[J].SpaceMedicine&MedicalEngineering,2000,13(6):431434[4]LeeEM,Kim,NH,TrangNT,etal.Respiratoryratedetectionalgorithmsbyphotoplethysmography[C].30thAnnualInternationalIEEE/EMBSConferenceVancouver,BritishColumbia,Canada,2008:11401143.[5]LIXuan,ZHOUQingfeng,HEChaojin,etal.Comparisonanalysisofempiricalmodedecompositionandwaveletdecompositioninanalogsignals[J].JournalofGeology,2009,33(1):7983.[6]NordenE.Huang,ZhengShen,StevenR.Long,etal.TheempiricalmodedecompositionandtheHilbertspectrumfornonlinearandnonstationarytimeseriesanalysis[C].ProceedingsoftheRoyalSocietyofLondonSERIESA,1998:454,903995.[7]LIXiangjun.AQRSdetectionalgorithmbasedontheempiricalmodedecomposition[J].JournalofUniversityofElectronicScienceandTechnologyofChina,2007,36(4):795797.[8]NimunkarAJ,TompkinsWJ.RpeakdetectionandsignalaveragingforsimulatedstressECGusingEMD[C].Proceedingsofthe29thAnnualInternationalConferenceoftheIEEEEMBSCitInternationale,Lyon,France,2007:12611264.[9]LINShaojie,ZHANGPandeng,WUKa,ietal.AnalgorithmforextractingECGfeaturebasedonempiricalmodedecomposition[J].JournalofBiomedicalEngineeringResearch,2008,26(2):328330.[10]KimJM,ArakawaK,VonLintelT.Useofthepulsewavemonitorasameasurementofdiagnosticsympatheticblockandofsurgicalsympathectomy[J].AnesthAnalg,1995,54:289296.[11]DUJuan,NINGGangmin,ZHANGShaowen,etal.Spectralanalysisofpulsetransittimevariability[J].SpaceMedicine&MedicalEngineering,2007,20(1):3742.[:,,,]28223