胡锡坤 金 添
(国防科学技术大学电子科学与工程学院 长沙 410073)
基于自适应小波尺度选择的生物雷达呼吸与心跳分离方法
胡锡坤 金 添*
(国防科学技术大学电子科学与工程学院 长沙 410073)
当使用生物雷达进行生命体征参数提取时,由于心脏跳动产生的位移形变很小,回波较为微弱,而呼吸带动的胸腔起伏回波强度较大,基于简单傅里叶变换的周期性信息检测,往往无法有效提取心跳信号。采用小波变换的方法可以较好地分离出含有心跳、呼吸运动信息的信号分量,但小波尺度的选择对于不同的场景存在细微差异,影响到了分离的效果。针对这一问题,该文采用Morlet二进小波变换,提出了一种基于信噪比阈值定标的自适应小波尺度选择方法,有效解决了不同场景的呼吸心跳分离问题。最后通过实测结果验证了算法的准确性和可行性。
生物雷达;自适应小波选择;呼吸和心跳信号分离
生物雷达[1]通过发射电磁波信号穿透衣物等非金属障碍物,分析接收到的回波信号,对被测主体进行非接触式的体征探测,而且可以给出生命体心肺活动的有关医学信息,具有探测距离较远,不易受到天气、光线(如昼夜)、温度等实际因素影响的优点。生物雷达在医疗监护领域譬如非接触式呼吸心跳监测[2]、心律失常检测[3]、病态呼吸模式检测[4],尤其是在婴儿猝死综合症和睡眠呼吸中止症方面有广阔的应用[5],同时在野战救护、反恐处突和应急搜索等领域有广泛的应用前景[6],成为了近年来雷达领域的研究热点。早期的生物雷达以单频连续波雷达为主,多用于生命检测[7],随后发展出的超宽带技术由于具有良好的穿透能力,辅以杂波相消,信号分离技术可实现对生命体的呼吸频率和心率的检测,从而实现诸如呼吸模式研究、心律失常检测、胸腔形变检测等更为复杂的应用[8]。
在生命特征提取的应用中,呼吸和心跳通常被建模成正弦信号进行表征,认为心率、呼吸频率等同于它们频谱函数中的基频值,需要从有限长信号中对其进行分离。常规的生物雷达信号处理方法是首先使用带通滤波器将呼吸和心跳信号进行分离,之后利用离散傅里叶变换进行心率和呼吸频率的参数提取[9]。基于多分辨率的特性,小波变换可用作滤波器,通过对小波尺度等参数的选取实现对呼吸和心跳信号的分离。文献[10]和文献[11]采用小波变换处理雷达回波信号,提取生命体征的瞬时频率和微多普勒信息。文献[12]通过计算雷达回波信号各小波系数的自相关,再求得峰值比最大的自相关序列,以此峰值计算心率。该算法能有效提取周期性信息,得到平均心率信息。由于受外界环境噪声以及监测条件变化的影响,呼吸和心跳信号所对应的小波尺度无法预先设定,小波变换并不能一直很好地用作生命体征信号的提取。因此,需要针对小波尺度的选择方法进行相应的研究。本文提出了一种以信噪比为准则的自适应小波尺度选择方法,实现了不同场景下的非接触式生命体征信号分离。
2.1 频率分离与小波尺度对应关系
小波变换具有多分辨率特性,可以由粗到精的逐步观察信号,通过选择适当的尺度因子,确定一组伸缩窗,同时选择恰当的母小波,则可以在时域和频域表征信号的局部特征。信号f(t)的连续小波变换(Continuous-Wavelet Transform, CWT):
其中,ψ为母小波,a为尺度因子,b为位移因子,分别决定小波的时频窗在频域和时域的位置。小波分析中的尺度因子a的变化不仅改变小波的频谱结构,而且改变其窗口的大小和形状,较大的尺度因子对应低频段,此时重构信号频率分辨率高,时间分辨率低;反之,较小的尺度因子对应高频段,此时重构信号频率分辨率低,时间分辨率高,这就是它的多分辨率特性。对于处于不同频率带上的呼吸和心跳信号,可以利用不同的小波尺度选择进行分离。
采用小波技术提取生命体征信号需要对信号时频域都有较高的刻画,在众多母小波中,Morlet小波有以下3个特点:一是对于低频率振动信号具有很好的时间和频率定位功能,可以较好地揭示信号的振动模式;二是其傅里叶变换为Gauss函数,是一个特性较理想的带通滤波器;三是具有良好的时频局部化特性,并且对称性较好,其小波变换系数的相位信息比较规则[13–15]。所以本文选择Morlet母小波来进行时频分析。
Morlet基本小波的时域和频域表达式为:
a对应的Morlet小波函数的时域和频域表达式为[14]:
其中,fb是母小波的带宽参数,可预先设为定值,k为常数,与信号采样时间和采样点数有关。f0为小波中心频率,表征了小波的总体特性。值得注意的是实际应用中的CWT需要事先定义尺度范围a,因此,需要采用合适的方法来确定需要估计的频率分量的尺度范围。根据式(5),对于给定的a,当akf=f0时,Ga(f)的值最大,使得呼吸与心跳回波在该小波尺度下的能量达到最大即为所选尺度范围。因此,对于呼吸分量fr和心率分量fh,分别存在如下关系[13–18]:
其中,k表示信号采样间隔,为慢时间维的脉冲重复频率fs的导数。由式(6)可得,不同频率下的尺度选择与母小波的选择、中心频率选取及采样频率有关。在母小波中心频率f0和采样频率fs已知条件下,呼吸频率fr和心率fh对应着不同的小波尺度ar和ah。
2.2 小波尺度粗集合
在式(6)的基础上,依据呼吸和心跳的频率范围可以初步确定呼吸信号和心跳信号小波尺度范围,首先假设呼吸频率范围和心率范围分别为:
为了方便计算机处理,在CWT的基础上通常将尺度参数a离散化,这种小波变换被称为二进小波变换(Dyadic Wavelet Transform, DWT),即将小波尺度特殊离散化a=2m。将式(7)代入式(6)可确定包含有整个呼吸频率和心率范围的整数形式的小波尺度粗集合:
2.3 小波尺度精集合
式(9)所示集合内每一个小波尺度均可以重构出一组信号:
其中,si表示小波尺度为时的重构信号,sj 表示小波尺度为时的重构信号。每一个重构信号经快速傅里叶变换后可得其频谱峰值,即为对应呼吸频率和心率估计值:
同样可以更新该尺度下的重构信号频率估计值集合为:
2.4 信噪比阈值确定
首先给出信噪比的计算公式,假设小波尺度对应重构信号的频谱峰值频率为则重构信号的信噪比为[19]:
其中,B代表信号周期图频率分辨率,由具体FFT点数和采样率决定。P(f)表示在频率为f 时的信号输出功率。则小波尺度范围对应重构信号的信噪比集合为:
其中,SNRi表示小波尺度时重构信号si的信噪比,SNRj表示小波尺度时重构信号sj的信噪比。
(1) 理想测量条件
理想条件一般指近距离非接触测量,大约在0.1 m左右进行探测。若已知参考呼吸频率为fr,参考心率为fh,分别与式(14)中重构频率估计值进行比较,计算其准确率:
其中,φi和θj分别为小波尺度时重构信号si的准确率和小波尺度时重构信号sj的准确率。理想情况下,通过DWT可以100%重构恢复呼吸和心跳信号,因此使得φi=1和θj=1时的小波尺度即为最终选择尺度:
(2) 边界测量条件
在给定准确度σ下,雷达系统具有一个测量边界,在该边界测量的呼吸频率和心率准确度恰好均不低于σ。已知参考呼吸频率为fr,参考心率为fh,依据式(17)计算其准确率。因为心跳比较呼吸更加难以测量,根据测量前规定的准确度σ,在保证心跳具有σ的准确度前提下,可以确信呼吸准确度也不低于σ,因此使得φi最大和θj=σ时的小波尺度即为最终选择尺度:
综合两种特殊测量条件下的检测,得到适合绝大部分中间场景的呼吸尺度集合下的信噪比阈值和心跳尺度集合下的信噪比阈值SNRTh:
值得注意的是,式(23)中针对两个下限取最大值的过程将会舍弃少部分边界条件下测量场景,致使该场景下的检测得不到结果,实际操作过程中这一现象的出现取决于准确度σ的预先设定。
2.5 算法流程
通过理想条件测量和边界条件测量,确定了信噪比阈值后,算法可以进行自适应的小波尺度选择。最终尺度集合内尺度i和j从最小尺度值开始递增并依次得到重构信号当重构信号的信噪比SNRi>SNRTr,则输出呼吸频率当SNRj>SNRTh,则输出心率如果尺度集合中的重构信号的信噪比不符合阈值要求,则需要改变测量环境重新测量。依据前面的理论推导,自适应小波尺度选择算法流程如图1所示。
3.1 雷达系统与实验场景
实验采用的冲激超宽带雷达系统如图2所示,天线采用Vivaldi天线,发射信号中心频率6.8 GHz,带宽2.3 GHz。
实验设置和测试设备如图3所示,选用了综合医疗器械DynaMap Suite-SA7925输出的呼吸数据和心跳数据作为生物雷达检测的呼吸和心跳的参考信号。
3.2 信噪比阈值定标
依据本文第2部分,在进行体征信号检测前,我们需要通过理想条件和边界条件的测量来确定信噪比阈值。依据算法流程,首先需要根据重构呼吸频率和心率范围,初步确定小波尺度范围Ar和Ah;然后计算落于呼吸频率和心率范围内的频率峰值,缩小尺度范围至最后根据该小波尺度下的信噪比分析得到阈值
图1 自适应小波尺度选择算法流程图Fig. 1 The flow chart of the adaptive wavelet scales selection algorithm
图2 冲激超宽带雷达系统Fig. 2 Impulse ultra-wideband radar system
图3 实验设置和测试设备Fig. 3 Experimental setup and measuring devices
(1) 理想测量条件
理想条件下人体距离天线约为0.1 m,被测目标为一名健康男性。
(a) 确定小波尺度粗集合
采用二进小波变换,选择Morlet母小波中心频率参数值为8,即f0=28Hz。雷达接收回波结果脉冲积累后以65 Hz的实际慢时间维采样频率进行采样存储,即fs=65 Hz,对于呼吸频率fr,其范围一般为0.13~0.65 Hz,心率fh一般为0.83~1.80 Hz,根据式(6),计算得:
可以确定小波尺度粗集合Ar={14, 15, 16, 17}和Ah={12, 13, 14, 15}。
(b) 确定小波尺度精集合
计算小波尺度范围内各尺度对应重构信号的峰值频率和信噪比,如表1所示。Ar集合中所有尺度对应重构信号的峰值频率均落于0.13~0.65 Hz呼吸频率范围内,而Ah集合中尺度14和15对应重构信号的峰值频率落于0.83~1.80 Hz心率范围外,则可更新得到小波尺度精集合以及
表1 理想测量条件下不同小波尺度下的峰值频率与信噪比Tab. 1 The peak frequency and SNR for each different wavelet scale under ideal measuring condition
(c) 确定信噪比阈值范围
已知呼吸频率的参考值为0.22 Hz,心跳频率参考值为1.25 Hz,比较呼吸频率最终尺度范围内的峰值频率,其中尺度为17时的峰值频率与参考呼吸频率达到了100%的准确率;比较心率最终尺度范围在尺度为12时的峰值频率与参考心率达到了100%的准确度,具体对比如图4所示。从而确定以及依据式(19),根据表1可以得到理想条件下呼吸信号和心跳信号的信噪比阈值为:
图4 理想条件下的重构生命体征信号与参考信号的比较Fig. 4 The comparison between reconstructed vital sign signals and reference signals
(2) 边界测量条件
首先设定边界条件下准确率σ=0.95,经过多次实验,对不同体型、性别和距离远近的人体进行检测,与参考信号进行比较,在保证呼吸频率和心率检测准确率均不小于95%的限制条件下,确定边界测量条件。与理想条件下的步骤类似,可以给出所选小波尺度所对应峰值频率符合限制条件下的信噪比,如表2所示。
表2 边界测量条件下求得小波尺度的峰值频率与信噪比Tab. 2 The peak frequency and SNR for each different wavelet scale under complex measuring condition
由式(21)和式(22)可得边界测量条件下的阈值范围为:
依据式(23),由式(25)和式(26)得到测量准确度不低于95%条件下的定标阈值范围:
在此范围内我们一般选取较为接近下限的值为信噪比阈值,即呼吸信号分离信噪比阈值为以及心跳信号分离信噪比阈值为SNRTh=–9.1,这样可以保证系统适用于较多的实验测量场景,如果阈值选取越大,则相应分离准确度越高,同时无结果输出的可能性也更高,因为适用的场景变小,实际中参考阈值范围的下限进行选取就可以较好地达到准确度高于σ的设定要求。
3.3 实测场景下呼吸和心跳信号分离结果
改变被测场景,选取一名健康女性作为被测目标,在人体相距天线1.0 m条件下进行生物雷达呼吸心跳检测。当信噪比阈值确定后,采用提出的自适应小波尺度选择方法对测试目标的呼吸和心跳信号进行分离,呼吸频率和心率的测量准确度均不低于95%。图5给出了实测场景下的自适应小波尺度选择后的呼吸和心跳信号分离结果,其中心跳信号选择的尺度为12,呼吸信号选择的尺度为16,测试结果证明算法较好地实现了呼吸与心跳信号小波尺度的选取,同时实现了较高准确度的呼吸和心跳检测,具有可操作性与实用性。
图5 自适应小波尺度选择下的呼吸与心跳分离结果Fig. 5 Separation results using adaptive wavelet scales selection algorithm
本文针对生物雷达呼吸与心跳分离方法中小波变换尺度选择不固定这一问题,从理论分析到实验验证,对实际条件下的信噪比阈值进行了确定,实现了呼吸与心跳信号的有用分离,该方法可以灵活设定准确度范围,使生物雷达系统达到较高精度的非接触测量。
[1]胡巍. 基于多普勒雷达的非接触式生命体征检测技术研究[D]. [博士论文], 中国科学技术大学, 2014: 1–5. Hu Wei. Non-contact vital sign detection based on Doppler radar[D]. [Ph.D. dissertation], University of Science and Technology of China, 2014: 1–5.
[2]Li Changzhi, Lin Jenshan, and Xiao Yanming. Robust overnight monitoring of human vital signs by a non-contact respiration and heartbeat detector[C]. IEEE EMBS Annual International Conference, New York, NY, USA, 2006: 2235–2238.
[3]Hu Wei, Zhao Zhangyan, Wang Yunfeng, et al.. Noncontact accurate measurement of cardiopulmonary activity using a compact quadrature Doppler radar sensor[J]. IEEE Transactions on Biomedical Engineering, 2014, 61(3): 725–735.
[4]Gu Changzhan and Li Changzhi. Assessment of human respiration patterns via noncontact sensing using Doppler multi-radar system[J]. Sensors, 2015, 15(3): 6383–6398.
[5]Gu Changzhan. Short-Range noncontact sensors for healthcare and other emerging applications: A review[J]. Sensors, 2016, 16(8). Pii: E1169. DOI: 10.3390/s16081169.
[6]Chen Kun-Mu, Huang Yong, Zhang Jianping, et al.. Microwave life-detection systems for searching humansubjects under earthquake rubble or behind barrier[J]. IEEE Transactions on Biomedical Engineering, 2000, 47(1): 105–114.
[7]Gu Changzhan and Li Changzhi. From tumor targeting to speech monitoring: accurate respiratory monitoring using medical continuous-wave radar sensors[J]. IEEE Microwave Magazine, 2014, 15(4): 66–76.
[8]Fontana R J. Recent system applications of short-pulse ultra-wideband (UWB) technology[J]. IEEE Transactions on Microwave Theory and Techniques, 2004, 52(9): 2087–2104.
[9]Lazaro A, Girbau D, and Villarino R. Analysis of vital signs monitoring using an IR-UWB radar[J]. Progress In Electromagnetics Research, 2010, 100: 265–284.
[10]Tariq A and Ghafouri-Shiraz H. Vital signs detection using Doppler radar and continuous wavelet transform[J]. Journal of Nursing Science, 2011, 10(5): 285–288.
[11]Tariq A and Ghafouri-Shiraz H. Doppler radar vital signs monitoring using wavelet transform[C]. 2010 Antennas & Propagation Conference, Loughborough, 2010: 293–296.
[12]Morgan D and Zierdt M. Novel signal processing techniques for Doppler radar cardiopulmonary sensing[J]. Signal Processing, 2009, 89(1): 45–66.
[13]易丽, 何怡刚, 方葛丰, 等. 基于Morlet小波在雷达信号脉内特征提取中的应用[J]. 计算机应用研究, 2013, 30(1): 172–175. Yi Li, He Yigang, Fang Gefeng, et al.. Application of modified Morlet wavelet in feature extraction of radar signal to intra-pulse[J]. Application Research of Computers, 2013, 30(1): 172–175.
[14]薛蕙, 杨仁刚. 基于连续小波变换的非整数次谐波测量方法[J].电力系统自动化, 2003, 27(5): 49–53. Xue Hui and Yang Rengang. A novel method for noninteger harmonics measurement using continuous wavelet transform[J]. Automation of Electric Power Systems, 2003, 27(5): 49–53.
[15]卢新城, 龚沈光, 周骏, 等. 基于Morlet小波的高分辨信号频谱估计[J]. 武汉理工大学学报(交通科学与工程版), 2002, 27(5): 61–74. Lu Xincheng, Gong Shenguang, Zhou Jun, et al.. High resolution signal spectrum estimation based on Morletwavelet[J]. Journal of Wuhan University of Technology (Transportation Science & Engineering), 2002, 27(5): 61–74.
[16]Yang Yushan and Li Yuanyuan. Spectral analysis of geomagnetic data with continuous wavelet transform[C]. International Congress on Image and Signal Processing, 2009: 1–4.
[17]Rueda J L, Juarez C A, and Erlich I. Wavelet-based analysis of power system low-frequency electromechanical oscillations[J]. IEEE Transactions on Power Systems, 2011, 26(3): 1733–1743.
[18]Ball John E and Tolley Alan. Low SNR radar signal detection using the continuous wavelet transform (CWT) and a Morlet wavelet[C]. 2008 IEEE Radar Conference, 2008: 1–6.
[19]Lazaro Antonio, Girbau David, and Villarino Ramon. Techniques for clutter suppression in the presence of body movements during the detection of respiratory activity through UWB radars[J]. Sensors, 2014, 14(2): 2595–2618.
Adaptive Wavelet Scale Selection-based Method for Separating Respiration and Heartbeat in Bio-radars
Hu Xikun Jin Tian
(College of Electronic Science and Engineering, National University of Defense Technology, Changsha 410073, China)
Extracting periodic heartbeat signals based on the traditional Fourier transform using a noncontact bio-radar is difficult because chest displacements caused by the heart are much smaller than those caused by respiration. Normally, they can be separated using the continuous wavelet transform; however, the miniscule difference of wavelet scale selection under different conditions may influence the separation performance to some extent. To solve this problem, this study proposes a method based on signal-to-noise ratio calibration to adaptively select the Morletdyadic wavelet scales and then separate the heartbeat signal from the respiration one using the selected scales, which can be applied to detect vital signs of different conditions. The experimental results have exhibited the accuracy and feasibility of the proposed method.
Bio-radar; Adaptive wavelet selection; Respiration and heartbeat separation
TN959
A
2095-283X(2016)05-0462-08
10.12000/JR16103
胡锡坤, 金添. 基于自适应小波尺度选择的生物雷达呼吸与心跳分离方法[J]. 雷达学报, 2016, 5(5): 462–469.
10.12000/JR16103.
Reference format: Hu Xikun and Jin Tian. Adaptive wavelet scale selection-based method for separating respiration and heartbeat in bio-radars[J]. Journal of Radars, 2016, 5(5): 462–469. DOI: 10.12000/JR16103.
胡锡坤(1994–),男,湖北天门人,国防科学技术大学硕士研究生,主要研究方向为超宽带雷达生命信号检测技术。
E-mail: xikung_hu@126.com
金 添(1980–),男,国防科学技术大学教授,博士生导师,主要从事隐蔽目标雷达成像与检测识别、新型微波传感器机理与系统实现等方面的研究工作。2009年获全国优秀博士学位论文奖,2010年入选教育部“新世纪优秀人才支持计划”,2014年获国际无线电科学联盟青年科学家奖。承担国家自然科学基金、武器装备探索等多项课题,获省部级科技进步一等奖1项、二等奖2项。“信号处理与系统”国家精品课程和资源共享课主讲教师,信号处理系列课程国家级教学团队主要成员。已发表论文100余篇,获授权国家发明专利5项,出版专著3部、译著1部、教材1部。
E-mail: tianjin@nudt.edu.cn
2016-09-13;改回日期:2016-10-26;
2016-11-11
*通信作者:金添 tianjin@nudt.edu.cn
国家自然科学基金(61271441)
Foundation Item: The National Natural Science Foundation of China (61271441)