本发明涉及医学信号处理技术领域,尤其涉及QRS波群提取方法及装置。
背景技术:
当前监测人体健康,通过可穿戴式设备提取心电信号中QRS波群(反应左右心室的电激动过程)的方法基于两种算法。第一种是PT,Pan-Tompkins算法,主要是将ECG,Electrocardiograph心电信号通过带通滤波器,差分器和平方模块,基于一阶微分的希尔伯特变换方法,在这种方法中,阈值可变,并且不受人为干扰;第二种是FDD算法,主要使用一个FIR带通滤波器,滤波器的系数至于分式的阶数有关。上面两种提取方法皆是基于一阶微分信号的检测方法,后续其他算法皆是对其二者的改进方法。考虑到可穿戴式设备由于电路的设计不能过于复杂以及不会选用成本昂贵的传感器,往往降低了QRS波群检测的准确性。
上述内容仅用于辅助理解本发明的技术方案,并不代表承认上述内容是现有技术。
技术实现要素:
本发明的主要目的在于提供一种QRS波群提取方法及装置,旨在解决现有技术中可穿戴式设备由于电路的设计不能过于复杂以及不会选用成本昂贵的传感器,往往降低了QRS波群检测的准确性的问题。
为实现上述目的,本发明提供的一种QRS波群提取方法,包括步骤:
获取原始信号;
结构元素化所述原始信号;
运用形态学操作结构元素化后的信号得到原始信号的特征信号;
从所述特征信号中提取QRS波群。
优选地,所述从所述特征信号中提取QRS波群的步骤包括:
从特征信号中确定QRS波群的起点和终点;
从特征信号中确定最大值点作为R1点;
从特征信号中的起点到R1点范围内确定最小值点为Q1点;
从特征信号中的R1点到终点范围内确定最小值点为S1点;
根据QRS波群的起点、R1点、Q1点、S1点和终点确定QRS波群。
优选地,所述从所述特征信号中提取QRS波群的步骤包括:
从特征信号中确定最小值点作为R2点;
从特征信号中的起点到R2点范围内确定最大值点为Q2点;
从特征信号中的R2点到终点范围内确定最大值点为S2点;
根据QRS波群的起点、R2点、Q2点、S2点和终点确定QRS波群。
优选地,所述方法还包括:
计算正极点的值为-Q1+(R1-Q1)+(R1-S1)-S1;
计算负极点的值为-Q2+(Q2-Q1)+(R1-S1)-S1;
判断正极点的值是否大于负极点的值;
在正极点的值小于负极点的值时,R-peak=R1,正极点下的为原始信号的QRS波群;
在正极点的值大于或等于负极点的值时,R-peak=R2,负极点下的为原始信号的QRS波群;其中peak为QRS波群的波峰值。
优选地,所述从所述特征信号中提取QRS波群的步骤之后,还包括:
获取结构元素的更新方式,更新结构元素,并按照更新后的结构元素结构化所述原始信号。
此外,为实现上述目的,本发明还提供一种QRS波群提取装置,包括:
获取模块,用于获取原始信号;
处理模块,用于结构元素化所述原始信号;
计算模块,用于运用形态学操作结构元素化后的信号得到原始信号的特征信号;
提取模块,用于从所述特征信号中提取QRS波群。
优选地,所述提取模块包括:
确定单元,用于从特征信号中确定QRS波群的起点和终点;
确定单元还用于从特征信号中确定最大值点作为R1点;确
定单元还用于从特征信号中的起点到R1点范围内确定最小值点为Q1点;
确定单元还用于从特征信号中的R1点到终点范围内确定最小值点为S1点;
提取单元,用于根据QRS波群的起点、R1点、Q1点、S1点和终点确定QRS波群。
优选地,所述确定单元,还用于从特征信号中确定最小值点作为R2点;
确定单元还用于从特征信号中的起点到R2点范围内确定最大值点为Q2点;
确定单元还用于从特征信号中的R2点到终点范围内确定最大值点为S2点;
所述提取单元,还用于根据QRS波群的起点、R2点、Q2点、S2点和终点确定QRS波群。
优选地,所述装置还包括:判断模块,
所述计算模块,还用于计算正极点的值为-Q1+(R1-Q1)+(R1-S1)-S1;
计算模块还用于计算负极点的值为-Q2+(Q2-Q1)+(R1-S1)-S1;
所述判断模块,用于判断正极点的值是否大于负极点的值;
提取单元,还用于在正极点的值小于负极点的值时,R-peak=R1,正极点下的为原始信号的QRS波群;
提取单元还用于在正极点的值大于或等于负极点的值时,R-peak=R2,负极点下的为原始信号的QRS波群;其中peak为QRS波群的波峰值。
优选地,所述装置还包括:更新模块,用于获取结构元素的更新方式,更新结构元素;
所述处理模块,还用于按照更新后的结构元素结构化所述原始信号。
本发明通过运用形态学的方法可以起到突出QRS波群,抑制心电信号其他部分的作用,可以很好的从原始采集的心电信号中提取QRS波群,计算量小,且提高了提取的准确度。
附图说明
图1为本发明QRS波群提取方法的第一实施例的流程示意图;
图2为本发明一实施例中确定QRS波群的位置的流程示意图;
图3为本发明另一实施例中确定QRS波群的位置的流程示意图;
图4为本发明一实施例中提取的QRS波群的示意图;
图5为本发明QRS波群提取方法的第二实施例的流程示意图;
图6为本发明QRS波群提取装置的第一实施例的功能模块示意图;
图7为图6中提取模块一实施例的细化功能模块示意图;
图8为本发明QRS波群提取装置的第二实施例的功能模块示意图。
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明提供一种QRS波群提取方法。
参照图1,图1为本发明QRS波群提取方法的第一实施例的流程示意图。
在一实施例中,所述QRS波群提取方法包括:
步骤S10,获取原始信号;
在本实施例中,所述原始信号可以是ECG心电图仪采集到的心电信号,也可以是采集到的其他生理信号。在需要对采集的信号进行处理时,从心电图仪获取原始信号。在本发明一较佳实施例中,在获取到原始信号后,还可以对原始信号进行低通滤波处理,去掉显著的工频干扰,低通滤波处理的过程可以根据设置选择处理或不处理。
步骤S20,结构元素化所述原始信号;
在获取到原始信号后,初始化结构元素,结构元素化所述原始信号。结构元素的初始化:结构元素由5个基准点构成:起点、Q、R、S和终点。结构元素的长度为90ms,结构元素的幅度为心电信号初始2s信号段中最大值与最小值之差。在本发明一较佳实施例中,窗口滑动的规则是:如果没有检测出候选区域,则下一次窗口的起点定位当前窗口起点加上半秒窗口的位置;如果已经检测出了QRS波群,则下一次窗口的起点定位QRS波群的终点加上20ms的位置。
步骤S30,运用形态学操作结构元素化后的信号得到原始信号的特征信号;
在结构元素初始化后,进行形态学操作,运用形态学操作结构元素化后的信号得到原始信号的特征信号。具体的,形态学操作:按照下面公式1进行运算,FS=f-(f0g+f·g)/2------公式1;形态学的基本操作包括:
Top_Hat:That(f)=f-fοg。
其中,g为结构元素,f为窗口长度为1s的心电信号数据,FS为特征信号。通过上述形态学操作,可以起到突出QRS波群,抑制心电信号其他部分的作用。
步骤S40,从所述特征信号中提取QRS波群。
在得到特征信号后,从所述特征信号中提取QRS波群。找出特征信号FS中非零的部分,然后选择其中持续时间大于70ms的部分作为QRS的候选部分,然后从70ms候选的部分中按照以下方式确定QRS的位置。具体的,参考图2,QRS的位置确定的方式可以包括:
步骤S41,从特征信号中确定QRS波群的起点和终点;
步骤S42,从特征信号中确定最大值点作为R1点;
步骤S43,从特征信号中的起点到R1点范围内确定最小值点为Q1点;
步骤S44,从特征信号中的R1点到终点范围内确定最小值点为S1点;
步骤S45,根据QRS波群的起点、R1点、Q1点、S1点和终点确定QRS波群。
参考图3,QRS的位置确定的方式还可以包括:
步骤S46,从特征信号中确定最小值点作为R2点;
步骤S47,从特征信号中的起点到R2点范围内确定最大值点为Q2点;
步骤S48,从特征信号中的R2点到终点范围内确定最大值点为S2点;
步骤S49,根据QRS波群的起点、R2点、Q2点、S2点和终点确定QRS波群。
进一步地,在得到QRS波群的起点、R1、Q1、S1、终点、R2、Q2和S2后,计算正极点的值为-Q1+(R1-Q1)+(R1-S1)-S1;计算负极点的值为-Q2+(Q2-Q1)+(R1-S1)-S1;比较正极点的值与负极点的值,判断正极点的值是否大于负极点的值,在正极点的值小于负极点的值时,R-peak=R1,原始信号的QRS波群为正极点下确定的QRS波群;在正极点的值大于或等于负极点的值时,R-peak=R2,原始信号的QRS波群为负极点下确定的QRS波群。根据正负极点的判断结果确定QRS点的位置,得到所述原始信号的QRS波群,完成提取QRS波群的操作。其中peak为波峰的值,R为QRS波群的R的值。参考图4,为提取的QRS波群的示意图,其中圈出的为QRS波群的R点。
本实施例通过运用形态学的方法可以起到突出QRS波群,抑制心电信号其他部分的作用,可以很好的从原始采集的心电信号中提取QRS波群,计算量小,且提高了提取的准确度。
参照图5,图5为本发明QRS波群提取方法的第二实施例的流程示意图。基于上述QRS波群提取方法的第一实施例,所述方法还包括:
步骤S50,获取结构元素的更新方式,更新结构元素,并按照更新后的结构元素结构化所述原始信号。根据提取的QRS波群的信息对结构元素进行更新。
结构元素更新:Q R S终点都以起点为基准,按照下面的公式进行更新:
Newloc=(1-a)*Curr_loc+a*extractedloc;
Newam==(1-a)*Curr_am+a*extractedam;其中,Curr_loc和Curr_am分别指的是在现有结构元素中相对于起始点的位置和幅度;extractedloc和extractedam分别是指在提取的QRS波群中相对于起点的位置和幅度;而a按照下面公式更新:a=a-b,NPA>OPA*1.1;a=a+b,NPA<OPA*0.9;a=0.3其他条件,a的初始值为0.9,b为0.05,NPA和OPA分别指当前一次和前一次的PA值,PA=∑ie候选部分|FS(i)|。本实施例通过对结构元素的不断更新,使得获取QRS波群的位置的准确度进一步提高,提高医疗信号检测的准确性。
本发明进一步提供一种QRS波群提取装置。
参照图6,图6为本发明QRS波群提取装置的第一实施例的功能模块示意图。
在一实施例中,所述QRS波群提取装置包括:获取模块10、处理模块20、计算模块30、提取模块40和判断模块50。
所述获取模块10,用于获取原始信号;
在本实施例中,所述原始信号可以是ECG心电图仪采集到的心电信号,也可以是采集到的其他生理信号。在需要对采集的信号进行处理时,从心电图仪获取原始信号。在本发明一较佳实施例中,在获取到原始信号后,还可以对原始信号进行低通滤波处理,去掉显著的工频干扰,低通滤波处理的过程可以根据设置选择处理或不处理。
所述处理模块20,用于结构元素化所述原始信号;
在获取到原始信号后,初始化结构元素,结构元素化所述原始信号。结构元素的初始化:结构元素由5个基准点构成:起点、Q、R、S和终点。结构元素的长度为90ms,结构元素的幅度为心电信号初始2s信号段中最大值与最小值之差。在本发明一较佳实施例中,窗口滑动的规则是:如果没有检测出候选区域,则下一次窗口的起点定位当前窗口起点加上半秒窗口的位置;如果已经检测出了QRS波群,则下一次窗口的起点定位QRS波群的终点加上20ms的位置。
所述计算模块30,用于运用形态学操作结构元素化后的信号得到原始信号的特征信号;
在结构元素初始化后,进行形态学操作,运用形态学操作结构元素化后的信号得到原始信号的特征信号。具体的,形态学操作:按照下面公式1进行运算,FS=f-(f0g+f·g)/2------公式1;形态学的基本操作包括:
Top_Hat:That(f)=f-fοg。
其中,g为结构元素,f为窗口长度为1s的心电信号数据,FS为特征信号。通过上述形态学操作,可以起到突出QRS波群,抑制心电信号其他部分的作用。
所述提取模块40,用于从所述特征信号中提取QRS波群。
在得到特征信号后,从所述特征信号中提取QRS波群。找出特征信号FS中非零的部分,然后选择其中持续时间大于70ms的部分作为QRS的候选部分,然后从70ms候选的部分中按照以下方式确定QRS的位置。具体的,参考图7,所述提取模块40包括:确定单元41和提取单元42,
所述确定单元41,用于从特征信号中确定QRS波群的起点和终点;
确定单元41还用于从特征信号中确定最大值点作为R1点;
确定单元41还用于从特征信号中的起点到R1点范围内确定最小值点为Q1点;
确定单元41还用于从特征信号中的R1点到终点范围内确定最小值点为S1点;
所述提取单元42,用于根据QRS波群的起点、R1点、Q1点、S1点和终点确定QRS波群。
进一步地,所述确定单元41,还用于从特征信号中确定最小值点作为R2点;
确定单元41还用于从特征信号中的起点到R2点范围内确定最大值点为Q2点;
确定单元41还用于从特征信号中的R2点到终点范围内确定最大值点为S2点;
所述提取单元42,还用于根据QRS波群的起点、R2点、Q2点、S2点和终点确定QRS波群。
进一步地,在得到QRS波群的起点、R1、Q1、S1、终点、R2、Q2和S2后,计算正极点的值为-Q1+(R1-Q1)+(R1-S1)-S1;计算模块30,还用于计算负极点的值为-Q2+(Q2-Q1)+(R1-S1)-S1;判断模块50,用于比较正极点的值与负极点的值,判断正极点的值是否大于负极点的值,计算模块30还用于在正极点的值小于负极点的值时,R-peak=R1,原始信号的QRS波群为正极点下确定的QRS波群;在正极点的值大于或等于负极点的值时,R-peak=R2,原始信号的QRS波群为负极点下确定的QRS波群。根据正负极点的判断结果确定QRS点的位置,得到所述原始信号的QRS波群,完成提取QRS波群的操作。其中peak为波峰的值,R为QRS波群的R的值。
本实施例通过运用形态学的方法可以起到突出QRS波群,抑制心电信号其他部分的作用,可以很好的从原始采集的心电信号中提取QRS波群,计算量小,且提高了提取的准确度。
参照图8,图8为本发明QRS波群提取装置的第二实施例的功能示意图。所述装置还包括:更新模块60,
所述更新模块60,用于获取结构元素的更新方式,更新结构元素;
所述处理模块20,还用于按照更新后的结构元素结构化所述原始信号。根据提取的QRS波群的信息对结构元素进行更新。
结构元素更新:Q R S终点都以起点为基准,按照下面的公式进行更新:
Newloc=(1-a)*Curr_loc+a*extractedloc;
Newam==(1-a)*Curr_am+a*extractedam;其中,Curr_loc和Curr_am分别指的是在现有结构元素中相对于起始点的位置和幅度;extractedloc和extractedam分别是指在提取的QRS波群中相对于起点的位置和幅度;而a按照下面公式更新:a=a-b,NPA>OPA*1.1;a=a+b,NPA<OPA*0.9;a=0.3其他条件,a的初始值为0.9,b为0.05,NPA和OPA分别指当前一次和前一次的PA值,PA=∑ie候选部分|FS(i)|。本实施例通过对结构元素的不断更新,使得获取QRS波群的位置的准确度进一步提高,提高医疗信号检测的准确性。
以上仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。