CN105326477B - 一种用于多种医学信号采集的呼吸及心跳同步信号提取方法及装置 - Google Patents
一种用于多种医学信号采集的呼吸及心跳同步信号提取方法及装置 Download PDFInfo
- Publication number
- CN105326477B CN105326477B CN201510720252.6A CN201510720252A CN105326477B CN 105326477 B CN105326477 B CN 105326477B CN 201510720252 A CN201510720252 A CN 201510720252A CN 105326477 B CN105326477 B CN 105326477B
- Authority
- CN
- China
- Prior art keywords
- signal
- value
- peak
- peak value
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明属于医学信号采集领域,具体涉及一种呼吸及心跳同步信号的提取方法及装置。本发明由信号采集模块、信号预处理模块和同步信号提取模块组成。其中信号采集模块完成对呼吸信号和体表心电图(ECG)信号的采集;信号预处理模块对呼吸信号和ECG信号进行放大滤波的预处理;同步信号提取模块通过识别呼吸和ECG信号的特征点,根据要求提取同步信号,即给出其他医学信号采集的时间窗口。本发明能够通过呼吸及心跳同步信号为其他医学信号的采集设定合适的时间窗口,可有效抑制呼吸、心跳对其他医学信号所造成的干扰,方法简单有效。
Description
技术领域
本发明属于医学信号采集技术领域,具体涉及一种呼吸及心跳同步信号的提取方法及装置。
背景技术
在某些医学信号采集时,由于呼吸活动和心脏跳动所造成的躯体振动及阻抗变化,会对被采集信号造成不易消除的干扰。
为了解决上述问题,在目前的相关诊断/治疗中,常用的方法有操作人员通过监视呼吸动态来指导患者在适当时机屏气来消除呼吸运动的影响。它的缺点是操作人员不能够主动控制患者的呼吸运动,造成患者配合不理想,对于儿童、呼吸功能障碍以及依从性差的患者,更加不适用。因此,在医学实践中,需要一种能够同时消除呼吸活动和心脏跳动对医学信号采集的影响的补偿方法。
发明内容
本发明针对多种医学信号采集中呼吸活动和心脏跳动对信号准确采集的影响问题,提出一种呼吸及心跳同步信号提取(即将呼吸/心跳周期中的特定时相,作为其他信号采样的时间窗)的方法,并提出相应装置,能简单有效地减少呼吸及心跳对其他医学信号采集的影响,提高信号采集的准确程度。
本发明提出呼吸及心跳同步信号提取方法,具体步骤包括:
(1)呼吸信号和体表心电图(ECG)信号的采集;
(2)对呼吸信号和ECG信号的预处理;
(3)呼吸和ECG同步信号的提取。
步骤(1)中,呼吸信号使用呼吸传感器采集,如可在胸腔恒压激励回路中串接测量电阻,该电阻的压降的变化与呼吸活动同步,即可得到呼吸信号;ECG信号通过标准II导联标测得到。
步骤(2)中,所述对对呼吸信号和ECG信号的预处理,就是对采集到的呼吸、ECG信号进行放大、滤波的预处理,以滤除信号中的工频干扰、高频噪声以及基线漂移等,同时进行信号的放大,以便后续的进一步处理。
步骤(3)中,呼吸和ECG同步信号的提取,包括:呼吸同步信号的提取,ECG同步信号的提取,时间窗同步选择;其中:
对于呼吸同步信号的提取,本发明提出了一种呼吸同步信号的提取算法,该算法包括:对于经过预处理的呼吸信号求取一阶导数和二阶导数,以此进行峰/谷值判断;同时引入不应期(不应期长度内的数据点不会再被标记为新的峰值或者谷值)避免伪峰干扰,完成对峰/谷值区间的标记,进而得到峰/谷值,并得到当前呼吸速率;最后通过选取合适的可选阈值,得到最终的呼吸同步信号。具体步骤为:
(1)对于经过预处理的呼吸信号求取一阶导数和二阶导数;
(2)选取合适的不应期长度(不应期长度内的数据点不会再被标记为新的峰值或者谷值),用以进行峰/谷值区间标记。一般不应期长度不超过呼吸周期的1/4;
(3)根据信号一阶导数和二阶导数,对信号逐点进行是否满足极值点(包括极大点和极小点)条件的判断;
(4)对于满足极大值条件的点,其前后一个不应期长度均被标记为峰值区间,每个连续峰值区间都被标记上一个唯一的正整数,用以求取每个区间峰值;
(5)对于满足极小值条件的点,其前后一个不应期长度均被标记为谷值区间,每个连续谷值区间都被标记上一个唯一的负整数,用以求取每个区间谷值;
由于当一个点满足极大值或者极小值条件时,其周围一个不应期长度内的数据点不会再被标记为新的峰值或者谷值区间,则伪峰不会被标记为新的峰值或者谷值区间,可以有效防止伪峰干扰,提高最终得到的同步信号的准确程度;
(6)在标记出的各峰值区间中求最大值,得到相应的峰值及其在时间轴上的位置;在标记出的各谷值区间中求最小值,得到相应的谷值及其在时间轴上的位置;
(7)根据得到的最新的4组峰值及谷值在时间轴上的位置,通过取平均的方式计算得到当前的呼吸频率;
(8)根据得到的当前的呼吸频率,选取合适的时间窗宽,用于与峰值进行比较参考。一般时间窗宽选取为峰值点前后各1/4个当前呼吸周期长度时,可以满足下述全部的阈值需要;S212中,阈值选取百分比较高时,也可以根据需要适当减小时间窗宽;
(9)根据需要,由操作者设置阈值(该可选阈值为峰-谷值的百分比,例如95%,90%,70%等),得到合适宽度的同步信号。对一组对应的峰值和谷值,记峰值为Ap,谷值为Av,选取的可选阈值为。对于峰值点周围时间窗宽长度内的点,若对应信号值高于,则同步信号中相应位置取为1;对于谷值点周围时间窗宽长度内的点,若对应信号值低于,则同步信号中相应位置取为-1;其他点在同步信号中相应位置取为0。
经过阈值限定,得到并输出需要的呼吸同步信号。
对于与ECG同步信号的提取,本发明提出了一种心电信号R波峰值检测算法,该算法包括:首先求取预处理后的ECG信号的一阶导数,再进行Hilbert变换,并选取合适的时间窗宽对变换后的数据进行分段。通过计算每个数据段的最大值和均方根值,并进行比较判断,自动得到每一段的阈值,进行峰值区间标记,进而得到峰值,计算当前心率,并引入第二阈值进行漏峰检测。最后通过选取合适的可选阈值,得到需要的心跳同步信号。具体步骤为:
(1)对于预处理后的ECG信号,计算其一阶导数;
(2)对信号一阶导数,作Hilbert变换,用以衰减P波和T波的幅度,并减少其他干扰影响;
(3)选取合适的时间窗宽,用以进行数据分段。时间窗宽度决定了延时的长度,因此时间窗宽不宜选得过大,但至少要保证大于一个心动周期,以防止出现严重漏峰现象。一般选取平静状态下的心动周期的两倍(1.6s)以上;
(4)根据选取的时间窗宽,对得到的信号进行数据分段操作;
(5)计算分段得到的每一个数据段的最大值和均方根值;
(6)根据当前得到的每个数据段的最大值和均方根值,与前一数据段的最大值和均方根值作比较,计算得到相应数据段的划分峰值区间的自动阈值;
(7)用自动阈值标记出每个峰值区间,每个数据段中大于阈值的数据点划入峰值区间;
(8)在标记出的各峰值区间内求最大值,得到相应区间的峰值及其在时间轴上对应的位置;
(9)根据求到的最新的4组峰值在时间轴上的位置,通过取平均的方式计算得到当前的心率,同时得到各组峰值间的R-R间期;
(10)判断当前两个峰值间的R-R间期是否超过了前两个峰值间的R-R间期的1.5倍:
若判断当前两个峰值间的R-R间期超过了前两个峰值间的R-R间期的1.5倍,则采用第二阈值检测当前的R-R间期内是否有遗漏的峰值。第二阈值稍低于对应数据段中得到的自动阈值;
若判断出用第二阈值检测出当前R-R间期内有漏峰,则用第二阈值对该峰值区间进行标记,每个数据段中大于阈值的数据点划入峰值区间。标记完成后,回到步骤(8),再次进行求峰值及检测漏峰步骤;
(11)完成漏峰检测后,根据步骤(9)中得到的当前心率,选取合适的时间窗宽,用于在步骤(12)中与峰值进行比较参考。一般时间窗宽选取为峰值点前后各1/4个当前心动周期长度;
(12)根据需要,由操作者设置阈值(该可选阈值为峰值的百分比,例如95%,90%等),得到合适宽度的同步信号。对一个特定的峰值,记峰值为Ap(Ap的值应为S309中求到的峰值点在时间轴上的位置对应到S301后得到的预处理后的ECG信号的值),选取的可选阈值为。对于峰值点周围时间窗宽长度内的点,若对应信号值高于,则同步信号中相应位置取为1;其他点在同步信号中相应位置取为0。
经过阈值限定,即得到需要的心跳同步信号。
时间窗同步选择,是根据上述提取得到的呼吸同步信号和ECG的同步信号,根据实际需要选择呼吸和心跳同时进行同步,或者仅对呼吸/心跳进行同步,得到最终用于指导其他医学信号采集/处理的时间窗同步信号。
本发明还提供基于上述方法的呼吸及心跳同步信号的提取装置,该装置包括信号采集模块、信号预处理模块、A/D转换器和同步信号提取模块组成。其中,信号采集模块用于采集呼吸信号和体表心电图(ECG)信号,即对应于完成步骤(1)的工作;信号预处理模块用于对呼吸信号和ECG信号进行放大、滤波等的预处理,即对应于完成步骤(2)的工作;A/D转换器把信号预处理模块输出的电信号转换为数字信号,送入同步信号提取模块;同步信号提取模块包括呼吸同步信号的提取算法模块、心电信号R波峰值检测算法模块、时间窗同步选择模块,通过识别呼吸和ECG信号的特征点,根据要求提取同步信号,给出其他医学信号采集的时间窗口,即对应于完成步骤(3)的工作,具体采用微处理器实现。
本发明的有益效果是:能够通过呼吸及/或心跳同步信号为其他医学信号的采集设定合适的时间窗口,可有效抑制呼吸、心跳对其他医学信号所造成的干扰,方法简单有效。且不需要受检者在诊断/治疗过程中进行屏气,也不需要操作者对受检者进行额外的呼吸指导,减少受检者的不适程度,同时有效减少受检者的呼吸活动和心脏跳动对其他医学信号采集的影响,提高了信号采集的准确程度。
附图说明
图1为本发明的一种系统结构框图。
图2为本发明的呼吸同步信号获得的流程图。
图3为本发明的心跳同步信号获得的流程图。
具体实施方式
本发明提出的呼吸及心跳同步方法及装置的总结构示意图如图1所示。系统需要采集的信号包括呼吸信号103及ECG信号104。最终经过本方法的同步选择,得到用以消除呼吸活动和心脏跳动影响的时间窗同步信号117。
呼吸信号采集模块101中,在胸腔前后方向,施加恒压激励,该激励回路中串接一个10Ω的呼吸信号测量电阻,由于呼吸过程中整个胸腔阻抗会发生周期性的变化,将导致呼吸信号测量电阻上的压降发生与呼吸同步的变化,即可测得呼吸信号103。对得到的呼吸信号103,用呼吸信号预处理模块106进行预处理。预处理模块106中包括放大器和一个带通滤波器,放大器的总放大倍数设置为50,带通滤波器的通频带选取为0.1Hz-10Hz,滤除呼吸信号103中的工频干扰、高频噪声等,同时进行信号的放大,得到预处理后的呼吸信号108。
心电图(ECG)信号104,则通过标准II导联标测102得到。通过II导联检测得到的ECG信号104,用ECG的预处理模块107进行预处理。预处理模块107中包括放大器和一个带通滤波器,放大器的总放大倍数设置为1000,带通滤波器的通频带选取为0.5Hz-40Hz,滤除ECG信号104中的工频干扰、高频噪声以及基线漂移等,同时进行信号的放大,得到预处理后的ECG信号109。
预处理后的呼吸信号108及ECG信号109,经过模数转换110,采入微处理器111进行进一步的分析处理。
在微处理器111中,预处理后的呼吸信号108,经过呼吸同步算法112的处理,得到呼吸同步信号113;而预处理后的ECG信号109,经过心电信号的R波峰值检测算法114的处理,得到心跳同步信号115。
如图1所示的呼吸同步算法112的流程图如图2所示。包括以下步骤:
在步骤S201中:微处理器111读入预处理后的呼吸信号108。
在步骤S202中:由于预处理的呼吸信号108中仍然存在一定的基线漂移的情况,采集两个呼吸周期以上的信号后,通过一个下截止频率为0.1Hz的高通滤波器滤除基线漂移。
在步骤S203中:对S202滤除基线漂移后的信号,计算其一阶导数。
在步骤S204中:对S203中得到的信号一阶导数,再次求导,得到预处理后的呼吸信号108的二阶导数。
在步骤S205中:选取合适的不应期长度(不应期长度内的数据点不会再被标记为新的峰值或者谷值),用以进行S207和S208中的峰/谷值区间标记。不应期长度不可太小,但也不能够超过呼吸周期的1/4。当受检者/患者平静呼吸时(即呼吸频率约为12-20次每分钟),可选取不应期长度在0.15-0.25s。
在步骤S206中:根据S203中得到的信号一阶导数,以及S204中得到的信号的二阶导数,对已读入并完成前述步骤处理的信号逐点进行是否满足极值点条件的判断。当该点相应一阶导小于0.001,二阶导小于0,且该点的值大于其前后相邻两点的值,则判断为满足极大值条件;该点相应一阶导小于0.001,二阶导大于0,且该点的值小于其前后相邻两点的值,则判断为满足极小值条件。
在步骤S207中:满足极大值条件的点,其前后一个不应期长度均被标记为峰值区间,每个连续峰值区间都被标记上一个唯一的正整数,用以在S209中求取每个区间峰值。
在步骤S208中:满足极小值条件的点,其前后一个不应期长度均被标记为谷值区间,每个连续谷值区间都被标记上一个唯一的负整数,用以在S209中求取每个区间谷值。由于当一个点满足极大值或者极小值条件时,其周围一个不应期长度内的数据点不会再被标记为新的峰值或者谷值区间,则伪峰不会被标记为新的峰值或者谷值区间,可以有效防止伪峰干扰,提高最终得到的同步信号的准确程度。
在步骤S209中:在S207标记出的各峰值区间中求最大值,得到相应的峰值及其在时间轴上的位置;在S208标记出的各谷值区间中求最小值,得到相应的谷值及其在时间轴上的位置。
在步骤S210中:根据S209中求到的最新的4组峰值及谷值在时间轴上的位置,通过取平均的方式计算得到当前的呼吸频率。并随着数据的不断读入不断更新。
在步骤S211中:根据S210中得到的当前的呼吸频率,选取合适的时间窗宽,用在S212中与峰值进行比较参考。一般时间窗宽选取为峰值点前后各1/4个当前呼吸周期长度时,可以满足S212中全部的阈值需要。S212中,阈值选取百分比较高时,也可以根据需要适当减小时间窗宽。
在步骤S212中:根据需要,由操作者设置阈值(该可选阈值为峰-谷值的百分比,例如95%,90%,70%等,系统默认阈值为峰-谷值的90%),得到合适宽度的同步信号。对一组对应的峰值和谷值,记峰值为Ap,谷值为Av,选取的可选阈值为。对于峰值点周围时间窗宽(S211中选取)长度内的点,若对应信号值(S202后得到的信号)高于,则同步信号中相应位置取为1;对于谷值点周围时间窗宽(S211中选取)长度内的点,若对应信号值(S202后得到的信号)低于,则同步信号中相应位置取为-1;其他点在同步信号中相应位置取为0。
在步骤S213中:经过S212中的阈值限定,得到并输出需要的呼吸同步信号113。
如图1所示的心电信号的R波峰值检测算法114的流程图如图3所示。包括以下步骤:
在步骤S301中:微处理器111读入预处理后的ECG信号109。
在步骤S302中:对S301读入的预处理后的ECG信号109,计算其一阶导数。
在步骤S303中:对S302中计算得到的信号一阶导数,对其作Hilbert变换,用以衰减P波和T波的幅度,并减少其他干扰影响。
在步骤S304中:选取合适的时间窗宽,用以在S305中进行数据分段。时间窗宽度决定了延时的长度,因此时间窗宽不宜选得过大,但至少要保证大于一个心动周期以防止出现严重漏峰现象。一般选取平静状态下的心动周期的两倍(1.6s)以上。本发明中选取了2s的时间窗宽。
在步骤S305中:根据S304中选取的时间窗宽,对S303中得到的信号进行数据分段操作。
在步骤S306中:计算S305中分段得到的每一个数据段的最大值和均方根值。
在步骤S307中:根据S306中得到的每个数据段的最大值和均方根值,与前一数据段的最大值和均方根值作比较,自动计算得到相应数据段的划分峰值区间的阈值。具体计算规则为:当前数据段均方值小于0.18倍的最大值时,阈值取1.6倍的均方值;当前数据段均方值大于0.18倍的最大值,且最大值小于上一数据段最大值的两倍时,阈值取0.39倍的当前数据段最大值;当前数据段均方值大于0.18倍的最大值,且最大值大于上一数据段最大值的两倍时,阈值取0.39倍的上一数据段的最大值。
在步骤S308中:用自动阈值标记出每个峰值区间,每个数据段中大于阈值的数据点划入峰值区间。
在步骤S309中:在S308标记出的各峰值区间内求最大值,得到相应区间的峰值及其在时间轴上对应的位置。
在步骤S310中:根据S309中求到的最新的4组峰值在时间轴上的位置,通过取平均的方式计算得到当前的心率。并随着数据的不断读入不断更新。同时得到各组峰值间的R-R间期。
在步骤S311中:判断当前两个峰值间的R-R间期是否超过了前两个峰值间的R-R间期的1.5倍。
在步骤S312中:若S311中判断当前两个峰值间的R-R间期超过了前两个峰值间的R-R间期的1.5倍,则采用第二阈值检测当前的R-R间期内是否有遗漏的峰值。第二阈值稍低于对应数据段在S307中得到的自动阈值,在本发明中第二阈值取为对应自动阈值的90%。
在步骤S313中:若在S312中,判断出用第二阈值检测出当前R-R间期内有漏峰,则用第二阈值对该峰值区间进行标记,每个数据段中大于阈值的数据点划入峰值区间。标记完成后,回到S309,再次进行求峰值及检测漏峰步骤。
在步骤S314中:完成漏峰检测后,根据S310中得到的当前心率,选取合适的时间窗宽,用在S315中与峰值进行比较参考。一般时间窗宽选取为峰值点前后各1/4个当前心动周期长度时,可以满足S315中全部的阈值需要。S315中,阈值选取百分比较高时,也可以根据需要适当减小时间窗宽。
在步骤S315中:根据需要,由操作者设置阈值(该可选阈值为峰值的百分比,例如95%,90%等,系统默认阈值为峰值的90%),得到合适宽度的同步信号。对一个特定的峰值,记峰值为Ap(Ap的值应为S309中求到的峰值点在时间轴上的位置对应到S301后得到的预处理后的ECG信号的值),选取的可选阈值为。对于峰值点周围时间窗宽(S314中选取)长度内的点,若对应信号值(S301后得到的信号)高于,则同步信号中相应位置取为1;其他点在同步信号中相应位置取为0。
在步骤S316中:经过S315中的阈值限定,得到并输出需要的心跳同步信号115。
如图1所示的微处理器111中,在得到呼吸同步信号113和心跳同步信号115后,将这两路同步信号送入时间窗同步选择模块116中。根据实际需要,操作者可以选择同时对呼吸和心跳进行同步,或者仅对呼吸/心跳进行同步。经过选择,得到最终需要的时间窗同步信号117,用以消除呼吸及/或心跳对其他医学信号造成的干扰。
Claims (4)
1.一种呼吸及心跳同步信号提取方法,具体步骤为:
(1)呼吸信号和ECG信号的采集;
(2)对呼吸信号和ECG信号的预处理;
(3)呼吸和ECG同步信号的提取;
步骤(1)中,呼吸信号使用呼吸传感器采集;ECG信号通过标准II导联标测得到;
步骤(2)中,所述对呼吸信号和ECG信号的预处理,就是对采集到的呼吸、ECG信号进行放大、滤波的预处理,以滤除信号中的工频干扰、高频噪声以及基线漂移,同时进行信号的放大,以便后续的进一步处理;
步骤(3)中,呼吸和ECG同步信号的提取,包括:呼吸同步信号的提取,ECG同步信号的提取,时间窗同步选择;
其特征在于,对于呼吸同步信号的提取,采用呼吸同步信号的提取算法,该算法包括:对于经过预处理的呼吸信号求取一阶导数和二阶导数,以此进行峰/谷值判断;同时引入不应期,避免伪峰干扰,完成对峰/谷值区间的标记,进而得到峰/谷值,并得到当前呼吸速率;最后通过选取合适的可选阈值,得到最终的呼吸同步信号;
对于ECG同步信号的提取,采用心电信号R波峰值检测算法,该算法包括:首先求取预处理后的ECG信号的一阶导数,再进行Hilbert变换,并选取合适的时间窗宽对变换后的数据进行分段;通过计算每个数据段的最大值和均方根值,并进行比较判断,自动得到每一段的阈值,进行峰值区间标记,进而得到峰值,计算当前心率,并引入第二阈值进行漏峰检测;最后通过选取合适的可选阈值,得到需要的心跳同步信号;
时间窗同步选择,是根据上述提取得到的呼吸同步信号和ECG的同步信号,按照实际需要选择呼吸和心跳同时进行同步,或者仅对呼吸/心跳进行同步,得到最终用于指导其他医学信号采集/处理的时间窗同步信号。
2.根据权利要求1所述的呼吸及心跳同步信号提取方法,其特征在于所述的呼吸同步信号的提取算法,其具体步骤为:
(1)对于经过预处理的呼吸信号求取一阶导数和二阶导数;
(2)选取合适的不应期长度,用以进行峰/谷值区间标记;所述不应期长度不超过呼吸周期的1/4;
(3)根据信号一阶导数和二阶导数,对信号逐点进行是否满足极值点条件的判断;
(4)对于满足极大值条件的点,其前后一个不应期长度均被标记为峰值区间,每个连续峰值区间都被标记上一个唯一的正整数,用以求取每个区间峰值;
(5)对于满足极小值条件的点,其前后一个不应期长度均被标记为谷值区间,每个连续谷值区间都被标记上一个唯一的负整数,用以求取每个区间谷值;
(6)在标记出的各峰值区间中求最大值,得到相应的峰值及其在时间轴上的位置;在标记出的各谷值区间中求最小值,得到相应的谷值及其在时间轴上的位置;
(7)根据得到的最新的4组峰值及谷值在时间轴上的位置,通过取平均的方式计算得到当前的呼吸频率;
(8)根据得到的当前的呼吸频率,选取合适的时间窗宽,用于与峰值进行比较参考,时间窗宽选取为峰值点前后各1/4个当前呼吸周期长度;
(9)根据需要,由操作者设置阈值,得到合适宽度的同步信号:对一组对应的峰值和谷值,记峰值为Ap,谷值为Av,选取的可选阈值为;对于峰值点周围时间窗宽长度内的点,若对应信号值高于,则同步信号中相应位置取为1;对于谷值点周围时间窗宽长度内的点,若对应信号值低于,则同步信号中相应位置取为-1;其他点在同步信号中相应位置取为0;
经过阈值限定,得到并输出需要的呼吸同步信号。
3.根据权利要求2所述的呼吸及心跳同步信号提取方法,其特征在于所述的心电信号R波峰值检测算法,其具体步骤为:
(1)对于预处理后的ECG信号,计算其一阶导数;
(2)对信号一阶导数,作Hilbert变换;
(3)选取合适的时间窗宽,用以进行数据分段;时间窗宽度至少大于一个心动周期;
(4)根据选取的时间窗宽,对得到的信号进行数据分段操作;
(5)计算分段得到的每一个数据段的最大值和均方根值;
(6)根据每个数据段的最大值和均方根值,与前一数据段的最大值和均方根值作比较,计算得到相应数据段的划分峰值区间的自动阈值;
(7)用自动阈值标记出每个峰值区间,每个数据段中大于阈值的数据点划入峰值区间;
(8)在标记出的各峰值区间内求最大值,得到相应区间的峰值及其在时间轴上对应的位置;
(9)根据求到的最新的4组峰值在时间轴上的位置,通过取平均的方式计算得到当前的心率,同时得到各组峰值间的R-R间期;
(10)判断当前两个峰值间的R-R间期是否超过了前两个峰值间的R-R间期的1.5倍:
若判断当前两个峰值间的R-R间期超过了前两个峰值间的R-R间期的1.5倍,则采用第二阈值检测当前的R-R间期内是否有遗漏的峰值;第二阈值稍低于对应数据段中得到的自动阈值;
若判断出用第二阈值检测出当前R-R间期内有漏峰,则用第二阈值对该峰值区间进行标记,每个数据段中大于阈值的数据点划入峰值区间;
标记完成后,回到步骤(8),再次进行求峰值及检测漏峰步骤;
(11)完成漏峰检测后,根据步骤(9)中得到的当前心率,选取合适的时间窗宽,用于在步骤(12)中与峰值进行比较参考;时间窗宽选取为峰值点前后各1/4个当前心动周期长度;
(12)根据需要,由操作者设置阈值,得到合适宽度的同步信号:对一个特定的峰值,记峰值为Ap,选取的可选阈值为,对于峰值点周围时间窗宽长度内的点,若对应信号值高于,则同步信号中相应位置取为1;其他点在同步信号中相应位置取为0;
经过阈值限定,即得到需要的心跳同步信号。
4.基于权利要求3所述的方法的呼吸及心跳同步信号提取装置,其特征在于包括:信号采集模块、信号预处理模块、A/D转换器和同步信号提取模块;其中,信号采集模块用于采集呼吸信号和体表ECG信号,即对应于完成步骤(1)的工作;信号预处理模块用于对呼吸信号和ECG信号进行放大、滤波的预处理,即对应于完成步骤(2)的工作;A/D转换器把信号预处理模块输出的电信号转换为数字信号,送入同步信号提取模块;同步信号提取模块包括呼吸同步信号的提取算法子模块、心电信号R波峰值检测算法子模块、时间窗同步选择子模块,对应于完成步骤(3)的工作,具体采用微处理器实现。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510720252.6A CN105326477B (zh) | 2015-10-30 | 2015-10-30 | 一种用于多种医学信号采集的呼吸及心跳同步信号提取方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510720252.6A CN105326477B (zh) | 2015-10-30 | 2015-10-30 | 一种用于多种医学信号采集的呼吸及心跳同步信号提取方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105326477A CN105326477A (zh) | 2016-02-17 |
CN105326477B true CN105326477B (zh) | 2018-10-16 |
Family
ID=55277407
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510720252.6A Expired - Fee Related CN105326477B (zh) | 2015-10-30 | 2015-10-30 | 一种用于多种医学信号采集的呼吸及心跳同步信号提取方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105326477B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108652611B (zh) * | 2018-05-17 | 2021-03-16 | 广东工业大学 | 一种心率的计算方法、装置、设备及计算机可读存储介质 |
US20210169350A1 (en) * | 2018-06-27 | 2021-06-10 | Shenzhen Darma Technology Co., Ltd. | Cardiac physiological parameter measuring method, device, terminal and computer storage medium |
CN110547802B (zh) * | 2019-09-11 | 2022-09-06 | 京东方科技集团股份有限公司 | 识别呼吸状态的装置 |
CN111329481A (zh) * | 2020-03-03 | 2020-06-26 | 中国科学院深圳先进技术研究院 | 生理参数确定方法、装置、生理参数检测设备及介质 |
US20210369176A1 (en) * | 2020-06-02 | 2021-12-02 | Pacesetter, Inc. | Methods, devices and systems for identifying false r-r intervals and false arrhythmia detections due to r-wave undersensing or intermittent av conduction block |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4446872A (en) * | 1977-09-08 | 1984-05-08 | Avl Ag | Method and apparatus for determining systolic time intervals |
CN103876733A (zh) * | 2014-03-12 | 2014-06-25 | 西安交通大学 | 用于心肺系统相位同步分析的系统和方法 |
CN104739418A (zh) * | 2014-11-24 | 2015-07-01 | 中国科学院苏州生物医学工程技术研究所 | 一种基于立体视觉的呼吸门控系统及控制方法 |
CN104799858A (zh) * | 2014-01-26 | 2015-07-29 | 包头市稀宝博为医疗系统有限公司 | 用于磁共振系统的心电呼吸外周门控系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006126020A2 (en) * | 2005-05-24 | 2006-11-30 | Anglia Polytechnic University Higher Education Corporation | Medical data signal processing systems |
-
2015
- 2015-10-30 CN CN201510720252.6A patent/CN105326477B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4446872A (en) * | 1977-09-08 | 1984-05-08 | Avl Ag | Method and apparatus for determining systolic time intervals |
CN104799858A (zh) * | 2014-01-26 | 2015-07-29 | 包头市稀宝博为医疗系统有限公司 | 用于磁共振系统的心电呼吸外周门控系统 |
CN103876733A (zh) * | 2014-03-12 | 2014-06-25 | 西安交通大学 | 用于心肺系统相位同步分析的系统和方法 |
CN104739418A (zh) * | 2014-11-24 | 2015-07-01 | 中国科学院苏州生物医学工程技术研究所 | 一种基于立体视觉的呼吸门控系统及控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105326477A (zh) | 2016-02-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105326477B (zh) | 一种用于多种医学信号采集的呼吸及心跳同步信号提取方法及装置 | |
CN110327036B (zh) | 从穿戴式心电中提取呼吸信号和呼吸频率的方法 | |
CN101732050B (zh) | 一种基于光电容积波的呼吸率监测方法 | |
CN103070683B (zh) | 基于生物电阻抗的睡眠呼吸模式识别方法及装置 | |
CN105147269B (zh) | 一种无创连续血压测量方法 | |
CN104545888A (zh) | 一种基于动态心电与呼吸波采集的睡眠呼吸暂停采集分析系统 | |
CN106037655A (zh) | 一种睡眠周期信号的分离及周期计算方法及装置 | |
CN110236508A (zh) | 一种无创血压连续监测方法 | |
CN107361770A (zh) | 睡眠呼吸暂停事件判别装置 | |
CN106473700A (zh) | 基于hrv频域分析的阻塞性睡眠呼吸暂停综合征的检测系统 | |
CN106308769A (zh) | 一种兼具呼吸波采集的脉搏波诊断系统 | |
Heise et al. | Refinement and evaluation of a hydraulic bed sensor | |
CN109222928B (zh) | 呼吸信号的提取方法、装置、处理设备和系统 | |
CN102860839A (zh) | 一种分布式可穿戴动态心音采集装置 | |
CN110680307A (zh) | 一种运动环境下基于脉搏波传导时间的动态血压监测方法 | |
Zhang et al. | Hear the heart: Daily cardiac health monitoring using Ear-ECG and machine learning | |
CN204765621U (zh) | 一种基于动态心电与呼吸波采集的睡眠呼吸暂停采集分析系统 | |
CN105943040A (zh) | 一种治疗中枢性睡眠呼吸暂停症的方法及装置 | |
CN106108850A (zh) | 心电数据库的干扰数据的识别方法和装置 | |
CN106388788A (zh) | 一种兼具心功能信息采集的脉搏波诊断系统 | |
CN109157220A (zh) | 一种基于多通道动态监测的呼吸指数提取系统及工作方法 | |
KR101755791B1 (ko) | 공진호흡주파수 및 위상 검출을 통한 최적의 호흡패턴 결정 방법, 시스템 및 컴퓨터 프로그램 | |
CN111248900A (zh) | 一种基于单通道心脑信息耦合分析方法及系统 | |
CN115024716B (zh) | 基于心率标签生成的心冲击图信号重建方法 | |
CN209733969U (zh) | 一种脉搏波传播时间的测量设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181016 Termination date: 20211030 |