CN112120688A - 心电信号处理方法、设备及计算机可读存储介质 - Google Patents
心电信号处理方法、设备及计算机可读存储介质 Download PDFInfo
- Publication number
- CN112120688A CN112120688A CN201910564383.8A CN201910564383A CN112120688A CN 112120688 A CN112120688 A CN 112120688A CN 201910564383 A CN201910564383 A CN 201910564383A CN 112120688 A CN112120688 A CN 112120688A
- Authority
- CN
- China
- Prior art keywords
- maternal
- signal
- electrocardiosignal
- fetal
- target
- 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.)
- Granted
Links
- 238000012545 processing Methods 0.000 title claims abstract description 25
- 238000003672 processing method Methods 0.000 title claims abstract description 22
- 230000008774 maternal effect Effects 0.000 claims abstract description 167
- 230000001605 fetal effect Effects 0.000 claims abstract description 58
- 210000001015 abdomen Anatomy 0.000 claims abstract description 17
- 238000000034 method Methods 0.000 claims description 32
- 238000007781 pre-processing Methods 0.000 claims description 17
- 238000001514 detection method Methods 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 8
- 230000003044 adaptive effect Effects 0.000 claims description 6
- 230000000747 cardiac effect Effects 0.000 claims description 6
- 238000003379 elimination reaction Methods 0.000 claims description 6
- 230000008030 elimination Effects 0.000 claims description 5
- 238000004891 communication Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 210000002458 fetal heart Anatomy 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 210000003754 fetus Anatomy 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000010247 heart contraction Effects 0.000 description 1
- 210000005240 left ventricle Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035935 pregnancy Effects 0.000 description 1
- 238000007639 printing Methods 0.000 description 1
- 210000005241 right ventricle Anatomy 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2503/00—Evaluating a particular growth phase or type of persons or animals
- A61B2503/02—Foetus
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Power Engineering (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种心电信号处理方法,该方法包括:获取胸部心电电极采集到的参考母体心电信号,及腹部心电电极采集到的母胎混合心电信号;确定所述参考母体心电信号中的参考母体QRS信号位置;根据所述参考母体QRS信号位置,确定所述母胎混合心电信号中的目标母体QRS信号位置;根据所述目标母体QRS信号位置,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。本发明还公开了一种心电信号处理设备和一种计算机可读存储介质。本发明能够提高胎儿及母亲心电信号识别的准确性,减少医生误诊。
Description
技术领域
本发明涉及信号处理技术领域,尤其涉及心电信号处理方法、设备及计算机可读存储介质。
背景技术
体表胎儿心电技术是一种无创同时监测母体心率和胎儿心率的技术手段,它将心电电极贴放在母体体表,获取母体、胎儿混合心电信号,然后对母胎混合心电信号进行滤波、识别和分离,进而得到母体心电信号和胎儿心电信号,最后分别对母体心电信号和胎儿心电信号进行分析可以得到波形分析结果和心率监测结果。
现有技术中,在获取体表胎儿心电信号时都是仅仅围绕着胎儿心脏周围(即孕妇腹部)贴放电极,用这种方式能够保证获取到幅度较大的胎儿心电信号,但这也会导致一个严重问题:由于腹部距离母体心脏较远,因此通过腹部电极采集到的母体心电信号的幅度会较小,在孕晚期,随着胎儿心电信号幅度的增大,可能会出现检测到的胎儿心电信号幅度比母体心电信号幅度更大的情况,这时就会导致将母亲的心电信号处理为胎儿的心电信号,将胎儿的心电信号处理为母亲的心电信号,这种误识别会影响医生判断,造成医生误诊。
发明内容
本发明的主要目的在于提出一种心电信号处理方法、设备及计算机可读存储介质,旨在提高胎儿及母亲心电信号识别的准确性,减少医生误诊。
为实现上述目的,本发明提供一种心电信号处理方法,所述心电信号处理方法包括如下步骤:
获取胸部心电电极采集到的参考母体心电信号,及腹部心电电极采集到的母胎混合心电信号;
确定所述参考母体心电信号中的参考母体QRS信号位置;
根据所述参考母体QRS信号位置,确定所述母胎混合心电信号中的目标母体QRS信号位置;
根据所述目标母体QRS信号位置,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。
优选地,所述确定所述参考母体心电信号中的参考母体QRS信号位置的步骤包括:
对所述参考母体心电信号进行预处理,得到参考母体心电包络信号;
获取所述参考母体心电包络信号中大于或等于第一预设阈值的的峰值位置;
根据所述大于或等于第一预设阈值的的峰值位置,确定所述参考母体QRS信号的起始位置和终止位置。
优选地,所述对所述参考母体心电信号进行预处理,得到参考母体心电包络信号的步骤包括:
通过预设的自适应滤波器,对所述参考母体心电信号进行去噪处理;
通过预设的差分滤波器和低通滤波器,对去噪处理后的所述参考母体心电信号进行滤波处理,得到参考母体心电包络信号。
优选地,所述根据所述参考母体QRS信号位置,确定所述母胎混合心电信号中的目标母体QRS信号位置的步骤包括:
获取所述参考母体QRS信号位置的起点和终点;
将所述参考母体QRS信号位置的起点和终点对应作为所述母胎混合心电信号中的目标母体QRS信号的起点和终点;
或者,根据所述参考母体QRS信号位置的起点和终点确定一个检测区间,在所述检测区间内识别所述母胎混合心电信号中的目标母体QRS信号,得到所述目标母体QRS信号位置。
优选地,所述根据所述目标母体QRS信号位置,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号的步骤包括:
根据所述目标母体QRS信号位置,获取所述母胎混合心电信号中的目标母体QRS信号波群;
根据所述目标母体QRS信号波群,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。
优选地,所述根据所述目标母体QRS信号位置,获取所述母胎混合心电信号中的目标母体QRS信号波群的步骤包括:
对所述母胎混合心电信号进行预处理,得到母胎混合心电包络信号;
在所述母胎混合心电包络信号中,在所述目标母体QRS信号位置获取所述目标母体QRS信号波群。
优选地,所述根据所述目标母体QRS信号波群,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号的步骤包括:
对所述母胎混合心电信号进行所述目标母体QRS信号波群消除处理,得到胎儿心电信号;
在所述母胎混合心电信号中消除所述胎儿心电信号,得到目标母体心电信号。
优选地,所述确定所述参考母体心电信号中的参考母体QRS信号位置的步骤之前,还包括:
获取所述参考母体心电信号的信号质量,并判断所述信号质量是否符合预设条件;
当所述信号质量符合预设条件时,执行步骤:确定所述参考母体心电信号中的参考母体QRS信号位置;
当所述信号质量不符合预设条件时,对所述母胎混合心电信号进行预处理,得到母胎混合心电包络信号;
获取所述母胎混合心电包络信号中大于或等于第二预设阈值的峰值位置;
根据所述大于或等于第二预设阈值的峰值位置,获取所述母胎混合心电信号中的目标母体QRS信号波群;
根据所述目标母体QRS信号波群,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。
此外,为实现上述目的,本发明还提供一种心电信号处理设备,所述心电信号处理设备包括:存储器、处理器及存储在所述存储器上并可在所述处理器上运行的心电信号处理程序,所述心电信号处理程序被所述处理器执行时实现如上所述的心电信号处理方法的步骤。
此外,为实现上述目的,本发明还提供一种计算机可读存储介质,所述计算机可读存储介质上存储有心电信号处理程序,所述心电信号处理程序被处理器执行时实现如上所述的心电信号处理方法的步骤。
本发明获取胸部心电电极采集到的参考母体心电信号,及腹部心电电极采集到的母胎混合心电信号;确定所述参考母体心电信号中的参考母体QRS信号位置;根据所述参考母体QRS信号位置,确定所述母胎混合心电信号中的目标母体QRS信号位置;根据所述目标母体QRS信号位置,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。由于参考母体心电信号和母胎混合心电信号是从同一个母体采集而来,因此参考母体心电信号中的参考母体QRS信号位置与母胎混合心电信号中的目标母体QRS信号位置具有同步性,因此本发明通过上述方式,能够对母胎混合心电信号中的目标母体QRS信号进行准确定位,从而保证了胎儿心电信号和目标母体心电信号的准确分离,从而本发明能够提高胎儿及母亲心电信号识别的准确性,减少医生误诊。
附图说明
图1是本发明实施例方案涉及的硬件运行环境的设备结构示意图;
图2为本发明心电信号处理方法第一实施例的流程示意图。
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
如图1所示,图1是本发明实施例方案涉及的硬件运行环境的设备结构示意图。
本发明实施例心电信号处理设备可以是胎儿监护仪。该心电信号处理设备包括一胸部心电电极和一腹部心电电极,其中,胸部心电电极用于贴放在母体心脏区域采集母体心电信号,腹部心电电极用于贴放在母体腹部采集母胎混合心电信号。
如图1所示,该设备可以包括:处理器1001,例如CPU,网络接口1004,用户接口1003,存储器1005,通信总线1002。其中,通信总线1002用于实现这些组件之间的连接通信。用户接口1003可以包括显示屏(Display)、输入单元比如键盘(Keyboard),可选用户接口1003还可以包括标准的有线接口、无线接口。网络接口1004可选的可以包括标准的有线接口、无线接口(如WI-FI接口)。存储器1005可以是高速RAM存储器,也可以是稳定的存储器(non-volatile memory),例如磁盘存储器。存储器1005可选的还可以是独立于前述处理器1001的存储装置。
本领域技术人员可以理解,图1中示出的终端结构并不构成对设备的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件布置。
如图1所示,作为一种计算机存储介质的存储器1005中可以包括操作***、网络通信模块、用户接口模块以及心电信号处理程序。
在图1所示的设备中,网络接口1004主要用于连接后台服务器,与后台服务器进行数据通信;用户接口1003主要用于连接客户端(用户端),与客户端进行数据通信;而处理器1001可以用于调用存储器1005中存储的心电信号处理程序,并执行以下操作:
获取胸部心电电极采集到的参考母体心电信号,及腹部心电电极采集到的母胎混合心电信号;
确定所述参考母体心电信号中的参考母体QRS信号位置;
根据所述参考母体QRS信号位置,确定所述母胎混合心电信号中的目标母体QRS信号位置;
根据所述目标母体QRS信号位置,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。
进一步地,处理器1001可以调用存储器1005中存储的心电信号处理程序,还执行以下操作:
对所述参考母体心电信号进行预处理,得到参考母体心电包络信号;
获取所述参考母体心电包络信号中大于或等于第一预设阈值的的峰值位置;
根据所述大于或等于第一预设阈值的的峰值位置,确定所述参考母体QRS信号的起始位置和终止位置。
进一步地,处理器1001可以调用存储器1005中存储的心电信号处理程序,还执行以下操作:
通过预设的自适应滤波器,对所述参考母体心电信号进行去噪处理;
通过预设的差分滤波器和低通滤波器,对去噪处理后的所述参考母体心电信号进行滤波处理,得到参考母体心电包络信号。
进一步地,处理器1001可以调用存储器1005中存储的心电信号处理程序,还执行以下操作:
获取所述参考母体QRS信号位置的起点和终点;
将所述参考母体QRS信号位置的起点和终点对应作为所述母胎混合心电信号中的目标母体QRS信号的起点和终点;
或者,根据所述参考母体QRS信号位置的起点和终点确定一个检测区间,在所述检测区间内识别所述母胎混合心电信号中的目标母体QRS信号,得到所述目标母体QRS信号位置。
进一步地,处理器1001可以调用存储器1005中存储的心电信号处理程序,还执行以下操作:
根据所述目标母体QRS信号位置,获取所述母胎混合心电信号中的目标母体QRS信号波群;
根据所述目标母体QRS信号波群,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。
进一步地,处理器1001可以调用存储器1005中存储的心电信号处理程序,还执行以下操作:
对所述母胎混合心电信号进行预处理,得到母胎混合心电包络信号;
在所述母胎混合心电包络信号中,在所述目标母体QRS信号位置获取所述目标母体QRS信号波群。
进一步地,处理器1001可以调用存储器1005中存储的心电信号处理程序,还执行以下操作:
对所述母胎混合心电信号进行所述目标母体QRS信号波群消除处理,得到胎儿心电信号;
在所述母胎混合心电信号中消除所述胎儿心电信号,得到目标母体心电信号。
进一步地,处理器1001可以调用存储器1005中存储的心电信号处理程序,还执行以下操作:
获取所述参考母体心电信号的信号质量,并判断所述信号质量是否符合预设条件;
当所述信号质量符合预设条件时,执行步骤:确定所述参考母体心电信号中的参考母体QRS信号位置;
当所述信号质量不符合预设条件时,对所述母胎混合心电信号进行预处理,得到母胎混合心电包络信号;
获取所述母胎混合心电包络信号中大于或等于第二预设阈值的峰值位置;
根据所述大于或等于第二预设阈值的峰值位置,获取所述母胎混合心电信号中的目标母体QRS信号波群;
根据所述目标母体QRS信号波群,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。
本发明心电信号处理设备的实施例与下述心电信号处理方法的各个实施例基本相同,此处不作赘述。
基于上述硬件结构,提出本发明心电信号处理方法实施例。
参照图2,图2为本发明心电信号处理方法第一实施例的流程示意图,所述方法包括:
步骤S10,获取胸部心电电极采集到的参考母体心电信号,及腹部心电电极采集到的母胎混合心电信号;
本实施例心电信号处理方法应用于心电信号处理设备,如胎儿监护仪。该心电信号处理设备包括一胸部心电电极和一腹部心电电极,其中,胸部心电电极用于贴放在母体心脏区域采集母体心电信号,腹部心电电极用于贴放在母体腹部采集母胎混合心电信号。为提高胎儿心电信号和母亲心电信号识别的准确性,本实施例采用胸部心电电极和腹部心电电极同时进行心电信号的采集。
具体地,首先,获取胸部心电电极采集到的参考母体心电信号,及腹部心电电极采集到的母胎混合心电信号,其中,由于参考母体心电信号是从母体心脏周围采集,此处离母体心脏距离近且离胎儿心脏距离远,因此该参考母体心电信号可以视为纯母体心电信号。在采集得到参考母体心电信号和母胎混合心电信号后,可以将采集到的心电信号数据放入指定缓冲区,以供胎儿监护***程序从该缓冲区中读取心电信号数据进行心电分析、心率计算等。
步骤S20,确定所述参考母体心电信号中的参考母体QRS信号位置;
心电信号中包含QRS(QRS wave complex)信号,也即QRS波群,QRS波群反映左、右心室除极电位和时间的变化,一个QRS波群代表一次心脏跳动周期。
在采集到参考母体心电信号后,确定该参考母体心电信号中的参考母体QRS信号位置,该参考母体QRS信号位置即参考母体QRS信号在时序上的位置。
在一实施方式中,上述步骤S20可以包括:对所述参考母体心电信号进行预处理,得到参考母体心电包络信号;获取所述参考母体心电包络信号中大于或等于第一预设阈值的的峰值位置;根据所述大于或等于第一预设阈值的的峰值位置,确定所述参考母体QRS信号的起始位置和终止位置。
上述对所述参考母体心电信号进行预处理,得到参考母体心电包络信号的步骤进一步包括:通过预设的自适应滤波器,对所述参考母体心电信号进行去噪处理;通过预设的差分滤波器和低通滤波器,对去噪处理后的所述参考母体心电信号进行滤波处理,得到参考母体心电包络信号。
具体地,首先通过一个预设的自适应滤波器消除参考母体心电信号中的噪声,如基线漂移干扰、工频干扰和肌电干扰等,然后通过差分滤波器和低通滤波器对消除干扰后的参考母体心电信号进行滤波处理,从而得到参考母体心电包络信号。
通过对参考母体心电信号进行去噪和滤波处理,为后续对该参考母体心电信号进行准确分析提供了前提保证。
在得到参考母体心电包络信号后,获取该参考母体心电包络信号中大于或等于第一预设阈值的峰值位置,该峰值位置对应参考母体QRS信号的峰值,根据该峰值位置结合参考母体心电包络信号波形在该峰值位置附近的变化趋势,可以确定参考母体QRS信号的起始位置和终止位置,即QRS波群在时序上的起点和终点。
由此,实现了对参考母体QRS信号位置的确定。
当然,确定参考母体QRS信号位置也可以采用其他方式,比如通过QRS信号的特征点识别出参考母体QRS信号,进而确定该参考母体QRS信号的位置,具体实施时可灵活设置。
步骤S30,根据所述参考母体QRS信号位置,确定所述母胎混合心电信号中的目标母体QRS信号位置;
该步骤中,根据上述确定的参考母体QRS信号位置,确定母胎混合心电信号中的目标母体QRS信号位置。由于参考母体心电信号和母胎混合心电信号是从同一个母体采集而来,因此参考母体心电信号中的参考母体QRS信号位置与母胎混合心电信号中的目标母体QRS信号位置具有同步性,因此根据参考母体QRS信号可以实现对目标母体QRS信号进行定位。
具体地,上述步骤S30可以进一步包括:获取所述参考母体QRS信号位置的起点和终点;将所述参考母体QRS信号位置的起点和终点对应作为所述母胎混合心电信号中的目标母体QRS信号的起点和终点;或者,根据所述参考母体QRS信号位置的起点和终点确定一个检测区间,在所述检测区间内识别所述母胎混合心电信号中的目标母体QRS信号,得到所述目标母体QRS信号位置。
比如,参考母体QRS信号位置在(t1,t2)区间范围内,则t1,t2可分别作为母胎混合心电信号中的目标母体QRS信号的起点和终点,这种对目标母体QRS信号进行定位的方法简单、方便,易于实现;或者,也可以将(t1-a,t2+b)作为检测区间,其中,a、b可以预先灵活设置,然后在母胎混合心电信号中检测(t1-a,t2+b)区间的信号段,在该信号段内识别目标母体QRS信号,进而得到目标母体QRS信号位置,其中,识别目标母体QRS信号的具体方法与上述识别参考母体QRS信号的方法类似,此处不作赘述,这种方法检测出的目标母体QRS位置准确、可靠。
步骤S40,根据所述目标母体QRS信号位置,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。
本实施例中,根据上述确定的目标母体QRS信号位置,可以识别目标母体QRS信号波群,进而根据目标母体QRS信号波群,可以从母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。
在分离出胎儿心电信号和目标母体心电信号后,可以分别根据胎儿心电信号和目标母体心电信号计算胎儿心率和母体心率,并输出到监护***程序,以完成心率曲线显示、打印以及存储等功能。
在本实施例中,由于参考母体心电信号和母胎混合心电信号是从同一个母体采集而来,因此参考母体心电信号中的参考母体QRS信号位置与母胎混合心电信号中的目标母体QRS信号位置具有同步性,因此本实施例通过上述方式,能够对母胎混合心电信号中的目标母体QRS信号进行准确定位,从而保证了胎儿心电信号和目标母体心电信号的准确分离,从而本发明能够提高胎儿及母亲心电信号识别的准确性,减少医生误诊。
进一步地,基于本发明心电信号处理方法第一实施例,提出本发明心电信号处理方法第二实施例。
在本实施例中,上述步骤S40可以进一步包括:根据所述目标母体QRS信号位置,获取所述母胎混合心电信号中的目标母体QRS信号波群;根据所述目标母体QRS信号波群,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。
具体地,所述根据所述目标母体QRS信号位置,获取所述母胎混合心电信号中的目标母体QRS信号波群的步骤包括:对所述母胎混合心电信号进行预处理,得到母胎混合心电包络信号;在所述母胎混合心电包络信号中,在所述目标母体QRS信号位置获取所述目标母体QRS信号波群。
其中,对母胎混合心电信号进行预处理的方式可以为:通过预设的自适应滤波器,对所述母胎混合心电信号进行去噪处理;通过预设的差分滤波器和低通滤波器,对去噪处理后的所述母胎混合心电信号进行滤波处理,得到母胎混合心电包络信号。之后,在该母胎混合心电包络信号中,截取目标母体QRS信号位置的信号段,该信号段即目标母体QRS信号波群。
进一步地,所述根据所述目标母体QRS信号波群,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号的步骤包括:对所述母胎混合心电信号进行所述目标母体QRS信号波群消除处理,得到胎儿心电信号;在所述母胎混合心电信号中消除所述胎儿心电信号,得到目标母体心电信号。
在本实施例中,在获取到目标母体QRS信号波群后,对母胎混合心电信号进行目标母体QRS信号波群消除处理,即可得到纯净的胎儿心电信号,之后,在该母胎混合心电信号中消除胎儿心电信号,即可得到纯净的目标母体心电信号。
通过上述方式,实现了胎儿心电信号和目标母体心电信号的准确分离。
进一步地,基于本发明心电信号处理方法第一、第二实施例,提出本发明心电信号处理方法第三实施例。
在本实施例中,上述步骤S20之前,还可以包括:获取所述参考母体心电信号的信号质量,并判断所述信号质量是否符合预设条件。
具体地,可以检测参考母体心电信号中参考母体QRS信号的宽度,若宽度在预设宽度范围内,则判定参考母体心电信号的信号质量符合预设条件,反之则判定参考母体心电信号的信号质量不符合预设条件。此外也可以通过检测参考母体心电信号的极值点个数、信号幅值、相邻两个QRS信号之前的间隔等来判断参考母体心电信号的信号质量是否符合预设条件,具体实施时可以灵活设置。
当所述信号质量符合预设条件时,执行步骤S20:确定所述参考母体心电信号中的参考母体QRS信号位置。
当信号质量不符合预设条件时,对所述母胎混合心电信号进行预处理,得到母胎混合心电包络信号;获取所述母胎混合心电包络信号中大于或等于第二预设阈值的峰值位置;根据所述大于或等于第二预设阈值的峰值位置,获取所述母胎混合心电信号中的目标母体QRS信号波群;根据所述目标母体QRS信号波群,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。
在本实施例中,若参考母体心电信号的信号质量不符合预设条件,则不能再使用参考母体心电信号中的参考母体QRS信号位置来定位母胎混合心电信号中的母体QRS位置,而只能利用母胎混合心电信号自身来进行定位,因为此时参考母体心电信号中的记录是不准确的。利用母胎混合心电信号自身来定位目标母体QRS信号的方式为:首先,对母胎混合心电信号进行去噪、滤波等预处理,得到母胎混合心电包络信号,然后获取母胎混合心电包络信号中大于或等于第二预设阈值的峰值位置,从该大于或等于第二预设阈值的峰值位置附近,获取母胎混合心电信号中的目标母体QRS信号波群。
在获取到目标母体QRS信号波群后,对母胎混合心电信号进行目标母体QRS信号波群消除处理,即可得到纯净的胎儿心电信号,之后,在该母胎混合心电信号中消除胎儿心电信号,即可得到纯净的目标母体心电信号。
通过上述方式,可以防止防止当胸部心电电极脱落或者受到干扰时检测到的母体心电信号对母胎混合心电信号的分析造成的影响,从而进一步提高了胎儿及母亲心电信号识别的准确性。
本发明还提供一种计算机可读存储介质。
本发明计算机可读存储介质上存储有心电信号处理程序,所述心电信号处理程序被处理器执行时实现如上所述的心电信号处理方法的步骤。
其中,在所述处理器上运行的心电信号处理程序被执行时所实现的方法可参照本发明心电信号处理方法各个实施例,此处不再赘述。
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者***不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者***所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括该要素的过程、方法、物品或者***中还存在另外的相同要素。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在如上所述的一个存储介质(如ROM/RAM、磁碟、光盘)中,包括若干指令用以使得一台终端设备(可以是手机,计算机,服务器,空调器,或者网络设备等)执行本发明各个实施例所述的方法。
以上仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (10)
1.一种心电信号处理方法,其特征在于,所述心电信号处理方法包括如下步骤:
获取胸部心电电极采集到的参考母体心电信号,及腹部心电电极采集到的母胎混合心电信号;
确定所述参考母体心电信号中的参考母体QRS信号位置;
根据所述参考母体QRS信号位置,确定所述母胎混合心电信号中的目标母体QRS信号位置;
根据所述目标母体QRS信号位置,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。
2.如权利要求1所述的心电信号处理方法,其特征在于,所述确定所述参考母体心电信号中的参考母体QRS信号位置的步骤包括:
对所述参考母体心电信号进行预处理,得到参考母体心电包络信号;
获取所述参考母体心电包络信号中大于或等于第一预设阈值的的峰值位置;
根据所述大于或等于第一预设阈值的的峰值位置,确定所述参考母体QRS信号的起始位置和终止位置。
3.如权利要求2所述的心电信号处理方法,其特征在于,所述对所述参考母体心电信号进行预处理,得到参考母体心电包络信号的步骤包括:
通过预设的自适应滤波器,对所述参考母体心电信号进行去噪处理;
通过预设的差分滤波器和低通滤波器,对去噪处理后的所述参考母体心电信号进行滤波处理,得到参考母体心电包络信号。
4.如权利要求1所述的心电信号处理方法,其特征在于,所述根据所述参考母体QRS信号位置,确定所述母胎混合心电信号中的目标母体QRS信号位置的步骤包括:
获取所述参考母体QRS信号位置的起点和终点;
将所述参考母体QRS信号位置的起点和终点对应作为所述母胎混合心电信号中的目标母体QRS信号的起点和终点;
或者,根据所述参考母体QRS信号位置的起点和终点确定一个检测区间,在所述检测区间内识别所述母胎混合心电信号中的目标母体QRS信号,得到所述目标母体QRS信号位置。
5.如权利要求1所述的心电信号处理方法,其特征在于,所述根据所述目标母体QRS信号位置,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号的步骤包括:
根据所述目标母体QRS信号位置,获取所述母胎混合心电信号中的目标母体QRS信号波群;
根据所述目标母体QRS信号波群,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。
6.如权利要求5所述的心电信号处理方法,其特征在于,所述根据所述目标母体QRS信号位置,获取所述母胎混合心电信号中的目标母体QRS信号波群的步骤包括:
对所述母胎混合心电信号进行预处理,得到母胎混合心电包络信号;
在所述母胎混合心电包络信号中,在所述目标母体QRS信号位置获取所述目标母体QRS信号波群。
7.如权利要求5所述的心电信号处理方法,其特征在于,所述根据所述目标母体QRS信号波群,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号的步骤包括:
对所述母胎混合心电信号进行所述目标母体QRS信号波群消除处理,得到胎儿心电信号;
在所述母胎混合心电信号中消除所述胎儿心电信号,得到目标母体心电信号。
8.如权利要求1至7中任一项所述的心电信号处理方法,其特征在于,所述确定所述参考母体心电信号中的参考母体QRS信号位置的步骤之前,还包括:
获取所述参考母体心电信号的信号质量,并判断所述信号质量是否符合预设条件;
当所述信号质量符合预设条件时,执行步骤:确定所述参考母体心电信号中的参考母体QRS信号位置;
当所述信号质量不符合预设条件时,对所述母胎混合心电信号进行预处理,得到母胎混合心电包络信号;
获取所述母胎混合心电包络信号中大于或等于第二预设阈值的峰值位置;
根据所述大于或等于第二预设阈值的峰值位置,获取所述母胎混合心电信号中的目标母体QRS信号波群;
根据所述目标母体QRS信号波群,从所述母胎混合心电信号中分离出胎儿心电信号和目标母体心电信号。
9.一种心电信号处理设备,其特征在于,所述心电信号处理设备包括:存储器、处理器及存储在所述存储器上并可在所述处理器上运行的心电信号处理程序,所述心电信号处理程序被所述处理器执行时实现如权利要求1至8中任一项所述的心电信号处理方法的步骤。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有心电信号处理程序,所述心电信号处理程序被处理器执行时实现如权利要求1至8中任一项所述的心电信号处理方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910564383.8A CN112120688B (zh) | 2019-06-25 | 2019-06-25 | 心电信号处理方法、设备及计算机可读存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910564383.8A CN112120688B (zh) | 2019-06-25 | 2019-06-25 | 心电信号处理方法、设备及计算机可读存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112120688A true CN112120688A (zh) | 2020-12-25 |
CN112120688B CN112120688B (zh) | 2024-08-09 |
Family
ID=73849496
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910564383.8A Active CN112120688B (zh) | 2019-06-25 | 2019-06-25 | 心电信号处理方法、设备及计算机可读存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112120688B (zh) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5372139A (en) * | 1991-06-24 | 1994-12-13 | Paul Benjamin Crilly | Method for suppressing a maternal electrocardiogram signal from a fetal electrocardiogram signal obtained with invasive and non-invasive techniques using an almost pure maternal electrocardiogram signal as a trigger |
US20050267377A1 (en) * | 2004-05-28 | 2005-12-01 | Dorothee Marossero | Maternal-fetal monitoring system |
CN102204817A (zh) * | 2011-05-31 | 2011-10-05 | 重庆大学 | 两导联的胎儿心电信号提取装置及方法 |
CN103263262A (zh) * | 2013-05-28 | 2013-08-28 | 捷普科技(上海)有限公司 | 测量胎儿心率的***及方法 |
CN104887220A (zh) * | 2015-06-18 | 2015-09-09 | 山东大学 | 一种由腹壁心电信号提取胎儿心电信号的方法及*** |
CN105411577A (zh) * | 2015-12-30 | 2016-03-23 | 深圳先进技术研究院 | 胎儿心电信号分离方法及*** |
GB201618027D0 (en) * | 2016-10-25 | 2016-12-07 | Guangdong University Of Technology | Method for blind extraction of fetal electrocardio based on time-frequency conversion |
CN106667478A (zh) * | 2016-12-07 | 2017-05-17 | 成都亿咖极科技有限公司 | 一种多导联联合检测的智能胎儿心电检测方法及*** |
WO2018023698A1 (zh) * | 2016-08-05 | 2018-02-08 | 深圳先进技术研究院 | 一种胎儿心电分离方法及装置 |
WO2018023697A1 (zh) * | 2016-08-05 | 2018-02-08 | 深圳先进技术研究院 | 一种胎儿心电分离方法及装置 |
CN107684423A (zh) * | 2016-08-05 | 2018-02-13 | 深圳先进技术研究院 | 一种胎儿心电分离方法及装置 |
CN108420425A (zh) * | 2017-08-22 | 2018-08-21 | 索思(苏州)医疗科技有限公司 | 穿戴式动态母体胎儿心电监测仪 |
CN108814590A (zh) * | 2018-03-23 | 2018-11-16 | 江苏华康信息技术有限公司 | 一种心电qrs波群的检测方法及其心电分析方法 |
CN109199375A (zh) * | 2018-11-30 | 2019-01-15 | 东南大学 | 一种无创胎儿心电检测装置及心电信号数据处理方法 |
-
2019
- 2019-06-25 CN CN201910564383.8A patent/CN112120688B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5372139A (en) * | 1991-06-24 | 1994-12-13 | Paul Benjamin Crilly | Method for suppressing a maternal electrocardiogram signal from a fetal electrocardiogram signal obtained with invasive and non-invasive techniques using an almost pure maternal electrocardiogram signal as a trigger |
US20050267377A1 (en) * | 2004-05-28 | 2005-12-01 | Dorothee Marossero | Maternal-fetal monitoring system |
CN102204817A (zh) * | 2011-05-31 | 2011-10-05 | 重庆大学 | 两导联的胎儿心电信号提取装置及方法 |
CN103263262A (zh) * | 2013-05-28 | 2013-08-28 | 捷普科技(上海)有限公司 | 测量胎儿心率的***及方法 |
CN104887220A (zh) * | 2015-06-18 | 2015-09-09 | 山东大学 | 一种由腹壁心电信号提取胎儿心电信号的方法及*** |
CN105411577A (zh) * | 2015-12-30 | 2016-03-23 | 深圳先进技术研究院 | 胎儿心电信号分离方法及*** |
WO2018023697A1 (zh) * | 2016-08-05 | 2018-02-08 | 深圳先进技术研究院 | 一种胎儿心电分离方法及装置 |
WO2018023698A1 (zh) * | 2016-08-05 | 2018-02-08 | 深圳先进技术研究院 | 一种胎儿心电分离方法及装置 |
CN107684423A (zh) * | 2016-08-05 | 2018-02-13 | 深圳先进技术研究院 | 一种胎儿心电分离方法及装置 |
GB201618027D0 (en) * | 2016-10-25 | 2016-12-07 | Guangdong University Of Technology | Method for blind extraction of fetal electrocardio based on time-frequency conversion |
CN106667478A (zh) * | 2016-12-07 | 2017-05-17 | 成都亿咖极科技有限公司 | 一种多导联联合检测的智能胎儿心电检测方法及*** |
CN108420425A (zh) * | 2017-08-22 | 2018-08-21 | 索思(苏州)医疗科技有限公司 | 穿戴式动态母体胎儿心电监测仪 |
CN108814590A (zh) * | 2018-03-23 | 2018-11-16 | 江苏华康信息技术有限公司 | 一种心电qrs波群的检测方法及其心电分析方法 |
CN109199375A (zh) * | 2018-11-30 | 2019-01-15 | 东南大学 | 一种无创胎儿心电检测装置及心电信号数据处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112120688B (zh) | 2024-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3795074B1 (en) | Method and device for automatically determining st segment of electrocardiogram signal based on artificial intelligence technology | |
JP5271718B2 (ja) | 複数の時間区分にわたり胎児及び産婦のecgを特定する方法 | |
CN111067511B (zh) | 一种心电数据中关键点的确定方法、装置及设备 | |
US20180011532A1 (en) | Methods and apparatuses for determining head movement | |
US9314180B2 (en) | Heart electrophysiological signal analysis system | |
CN107233093B (zh) | R波的检测方法、装置及电子设备 | |
CN114052746A (zh) | 运动心电数据处理方法、装置、计算机设备和存储介质 | |
US20090209874A1 (en) | Extrapolating ica knowledge from one epoch to another for improved fetal ecg separation | |
CN113499082A (zh) | Qrs波群检测方法、心电检测装置及可读存储介质 | |
CN113647959B (zh) | 一种心电波形信号的波形识别方法、装置及设备 | |
CN108937920B (zh) | 一种心室纤颤信号检测方法、***及心室纤颤检测装置 | |
CN112120688A (zh) | 心电信号处理方法、设备及计算机可读存储介质 | |
CN110613444B (zh) | 一种母胎混合心电信号的显示方法及装置 | |
CN108968950A (zh) | 一种数字化心电采集方法、装置及*** | |
CN110604564A (zh) | 电极***、胎儿心电信号采集、测量方法及装置 | |
CN104490389A (zh) | 一种用于侦测胎动的可穿戴装置及侦测胎动的方法 | |
JP6706996B2 (ja) | 生体信号処理装置、異常判別方法およびプログラム | |
CN110477904B (zh) | 一种房扑信号检测方法及装置 | |
CN110916649B (zh) | 一种长程心电散点图的处理装置、处理方法及检测装置 | |
CN110613442B (zh) | 一种监测母体子宫收缩活动的方法及装置 | |
Vega-Martínez et al. | Wavelet packet based algorithm for QRS region detection and R/S wave identification | |
US20210052179A1 (en) | Method for determining r peaks of electrocardiogram | |
WO2018023698A1 (zh) | 一种胎儿心电分离方法及装置 | |
CN111449649B (zh) | 心电数据的处理方法、装置及计算机可读存储介质 | |
CN115067967A (zh) | 心拍信号基准点确定方法、心拍类型识别方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |