CN102715897B - 一种基于胎心电频谱图确定胎儿平均心率的方法 - Google Patents
一种基于胎心电频谱图确定胎儿平均心率的方法 Download PDFInfo
- Publication number
- CN102715897B CN102715897B CN201210245589.2A CN201210245589A CN102715897B CN 102715897 B CN102715897 B CN 102715897B CN 201210245589 A CN201210245589 A CN 201210245589A CN 102715897 B CN102715897 B CN 102715897B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- average heart
- fetus
- heart rate
- frequency
- 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.)
- Active
Links
- 210000003754 fetus Anatomy 0.000 title claims abstract description 38
- 230000001605 fetal effect Effects 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000001228 spectrum Methods 0.000 claims abstract description 102
- 230000003993 interaction Effects 0.000 claims abstract description 12
- 210000002784 stomach Anatomy 0.000 claims description 9
- 208000007502 anemia Diseases 0.000 claims description 8
- 230000003595 spectral effect Effects 0.000 abstract 2
- 210000003815 abdominal wall Anatomy 0.000 abstract 1
- 210000002458 fetal heart Anatomy 0.000 description 7
- 210000004556 brain Anatomy 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 230000002452 interceptive effect Effects 0.000 description 3
- 230000033764 rhythmic process Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000003187 abdominal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008774 maternal effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003340 mental effect Effects 0.000 description 1
- 230000005622 photoelectricity Effects 0.000 description 1
- 230000035935 pregnancy Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明涉及一种基于胎儿心电频谱图确定胎儿平均心率的方法，它包括两个步骤：(1)读取从孕妇腹壁检测的胎儿心电信号、制作幅值谱或功率谱图；(2)在制作的幅值谱或功率谱图上通过人机交互确定胎儿平均心率。其特征是，在第(2)步人机交互中，在幅值谱或功率谱图上设置一组游标(须两个以上)，它们的频率坐标始终成倍数关系，操作者拖动这组游标中其中一个时，组内所有游标同时游动，并且各游标频率坐标总保持着始终不变的倍数关系。该方法将人工测算一组谱峰频率是否成倍数关系的问题变成了频率坐标成倍数关系的一组游标是否都经过谱峰的认定问题，可提高胎儿平均心率的确定效率5倍以上。
Description
技术领域
本申请涉及一种基于胎儿心电频谱图确定胎儿平均心率的方法。
胎儿在整个孕期处在非常艰难的生存环境中，胎心率是胎儿健康状况的一个重要指标。多普勒超声方法检测胎心率时，探头要对中胎儿心脏，在实际操作中多有不便，尤其在早期胎心比较小的时候。在孕妇腹壁体表安放电极采集胎儿心电信号，据此确定胎心率的方法，不要求电极对中胎心，操作较为方便，但也存在信号信噪比低、稳定地获取胎心率不容易的问题。
文献中介绍的基于胎儿心电信号确定平均胎心率的方法主要是时域的方法：(1)首先从孕妇腹壁检测的胎儿心电信号中分离去主要的干扰成分——母体心电成分，并尽量多地去掉基线漂移等其他干扰成分；(2)去除干扰后的保留信号波形，被认为主要是胎儿心电成分，基于差分阈值法或小波奇异性检测等方法，进行R波波峰检测，利用检测的各R波波峰，计算瞬时胎心率，进而计算胎儿平均心率。这样的时域方法，实现比较复杂；有时，保留信号会残留有较大母体心电分量，导致胎儿R波波峰检测失败，方法失效。
本申请涉及的是基于胎儿心电信号频谱图确定胎儿平均心率的问题，不同于以上所述的时域方法。
背景技术
利用胎儿心电信号频谱确定胎儿平均心率的方法，其步骤可概括为：(1)读取从孕妇腹壁体表检测的胎儿心电信号，直接制作其幅值谱或功率谱，或者对其进行诸如去干扰、降噪、变换等处理后再制作其幅值谱或功率谱；(2)在幅值谱或功率谱上，从左到右扫寻谱峰，每扫到一个谱峰，检查该谱峰频率的倍频处是否也存在谱峰，如果存在，并且该谱峰频率又不对应母体平均心率，则判定该谱峰对应频率为平均胎心率，否则为伪差谱峰。
在频谱上确定胎儿平均心率方法的原理和在频谱上确定母体平均心率方法的原理相同，都是利用谱峰的倍频关系。但是，确定母体平均心率比较容易，确定胎儿平均心率比较困难，特别是要稳定可靠地确定胎儿平均心率很困难，因为胎儿分量往往很小、在频谱上胎儿分量对应的一组谱峰往往不是很突出，存在很多伪差谱峰。
相对而言，从频谱上确定胎儿平均心率的比较放心可靠的方法，不是让计算机程序从左到右自动扫寻谱峰并自动判定一组谱峰频率是否有倍频关系，而是把幅值谱或功率谱图画出来、展现在人机界面上，由操作者通过人机交互来确定胎儿平均心率。比较容易想到的模式是：用眼睛和大脑来扫寻频谱图上的各谱峰，目测出其频率成倍数关系的一组谱峰，若其肯定不对应母体分量，则选中这组谱峰中最左边的一个谱峰，将其对应频率确定为胎儿平均心率。
胎儿心电技术发展严重滞后于需求，以上提到的比较容易想到的在频谱上通过人机交互确定胎儿平均心率的方法，虽然简单，也未见文献提到。另一方面，该简单方法一旦实施也会发现新的问题：由于很多时候胎儿心电分量很弱，在频谱上对应胎儿的谱峰不是很突出，频谱上会存在很多伪差谱峰；此时，单纯用眼睛和大脑来扫寻频谱图上的各谱峰，并且目测出一组谱峰的频率是否存在倍数关系，操作者会很累，须仔细分辨和心算，操作起来很慢，效率低。
参考文献：
[1]Karvounis，E.C.，C.Papaloukas，et al.″Fetal heart rate extraction from composite maternal ECG using complex continuous wavelet transform.″Computers in Cardiology，2004，Vol31：737-740.
[2]Karvounis，E.C.，M.G.Tsipouras，et al.″An automated methodology for fetal heart rate extraction from the abdominal electrocardiogram.″IEEE Transactions on Information Technology in Biomedicine，2007，11(6)：628-638.
[3]S.Kadambe，R.Murray，and G.Faye Boudreaux-Bartels，“Wavelet Transform-Based QRS Complex Detector，”IEEE Trans.on Biomed.Eng.，vol.46，no.7，pp.838-848，July.1999.
发明内容
发明目的。
提出一种基于胎儿心电频谱图确定胎儿平均心率的方法，改进容易想到的简单人机交互模式，以提高胎儿平均心率的估计效率，减轻操作者的劳动强度。
技术方案。
一种基于胎儿心电频谱图确定胎儿平均心率的方法，它包括两个步骤：(1)读取从孕妇腹壁检测的胎儿心电信号，制作幅值谱或功率谱图；(2)在制作的幅值谱或功率谱图上通过人机交互确定胎儿平均心率，其特征是，所述步骤(2)中的人机交互中，在幅值谱或功率谱图上设置一组联动游标，拖动时这组联动游标的频率坐标始终保持着倍数关系，当操作者将该组联动游标中最左边的游标移至频谱图上的一个谱峰时，观察位于该谱峰所在频率的倍频处的其他游标是否也对应谱峰，当发现一组谱峰正好具有倍数关系，且又不对应于母体的心电分量时，则该组谱峰最左边谱峰的频率确定为胎儿平均心率。
以上方法第(1)步骤中，读取孕妇腹壁检测胎心电信号后，可直接制作幅值谱或功率谱，也可以对其进行诸如去干扰、降噪、变换等各种处理后再制作其幅值谱或功率谱，此处不构成本方案的技术特征。
以上方法第(2)步骤中，由于频谱图上设置了一组始终具有倍频关系的游标，当操作者将最左边的游标移至一个谱峰时，则其他游标正好位于该谱峰频率的倍频处，非常容易观察倍频处是否也有谱峰；当发现一组谱峰的频率正好具有倍数关系，又不对应母体心电分量时，可将该组谱峰最左边谱峰的频率确定为胎儿平均心率。
人的眼睛和大脑擅长发现识别一个频谱图上的谱峰，但很不擅长测算谱峰频率之间是否存在倍数关系。前面提到的容易想到的简单人机交互模式，需要用人的眼睛和大脑来测算谱峰频率之间是否存在倍数关系，这是方法低效的根源；但是，本技术方案采用了新的交互模式，不再需要人的眼睛和大脑来测算谱峰频率之间是否是倍数关系，把人工测算一组谱峰频率是否成倍数关系的问题变成了频率坐标成倍数关系的一组游标是否都经过谱峰的认定问题，简言之，把一个计算问题变成了一个眼睛和大脑擅长的单纯的识别问题。本技术方案反映了人机交互确定平均心率的本质需要，是客观规律的自觉利用。
有益效果。
目前，尚未见文献报道基于胎儿心电频谱图人机交互确定胎儿平均心率的方法；当然，基于在胎儿心电频谱图上设置一组有倍频关系的游标来人机交互确定胎儿平均心率的方案更是未见报道。因此，本技术方案具有新颖性。
经大量测试，本技术方案在胎儿心电频谱图上设置一组有倍频关系的游标这一做法，与容易想到的简单人机交互方法相比，可提高胎儿平均心率的确定效率5倍以上。因此，提出的方法有实用性，有实质进展。
附图说明
图1，本发明基于胎儿心电频谱图确定胎儿平均心率的步骤及特征。
图2，一腹壁实测胎儿心电信号的波形图。
图3，幅值谱图及经过母体平均心率处谱峰的一组游标示意图。
图4，幅值谱图及经过一伪差谱峰的一组游标示意图。
图5，幅值谱图及经过胎儿平均心率处谱峰的一组游标示意图。
具体实施方式(实施例)
在医院对某一孕妇用日本光电标准12导心电图机——ECG-1350P将电极置于腹壁测了一组数据，采样率500Hz。
首先，读取从孕妇腹壁检测的一个导联信号的胎儿心电信号，选择以0.5Hz到100Hz为通带的带通滤波器去除基线漂移作为预处理，去除基线漂移后的信号波形图见图2所示，图中横轴为时间轴，时间刻度0、1、…、8的单位为秒；然后制作信号幅值谱，见图3、图4、图5所示。在图3、图4、图5幅值谱中，横坐标均为频率轴，频率刻度0、0.5、1、1.5、…、10的单位为Hz。
在制作的幅值谱图上设置了一组游标(此处为三个)，见图3、图4、图5中的三根竖线、它们的频率坐标始终成倍数关系。操作者从幅值谱图左边开始向右拖动游标，首先将游标落到第一个谱峰及其倍频位置处，见图3所示，发现正好对应3个谱峰，则易断定此组谱峰中的第1个谱峰对应母体平均心率；继续向右移动游标，将游标落到第2个谱峰及其倍频位置处，见图4所示，发现后面的两个游标没有对上谱峰，可断定此谱峰为1个伪差谱峰；继续向右移动游标，将游标落到第3个谱峰及其倍频位置处，见图5所示，可发现后面的两个游标也对上了谱峰，说明该组游标所指谱峰是对应胎儿心电的，则可将此时左边第1个游标指向的频率确定为胎儿平均心率。
Claims (1)
1.一种基于胎儿心电频谱图确定胎儿平均心率的方法，它包括两个步骤：(1)读取从孕妇腹壁检测的胎儿心电信号，制作幅值谱或功率谱图；(2)在制作的幅值谱或功率谱图上通过人机交互确定胎儿平均心率，其特征是，所述步骤(2)中的人机交互中，在幅值谱或功率谱图上设置了一组联动游标，拖动时这组联动游标的频率坐标始终保持着倍数关系，当操作者将该组联动游标中最左边的游标移至频谱图上的一个谱峰时，观察位于该谱峰所在频率的倍频处的其他游标是否也对应谱峰，当发现一组谱峰正好具有倍数关系，且又不对应于母体的心电分量时，则最左边谱峰的频率确定为胎儿平均心率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210245589.2A CN102715897B (zh) | 2012-07-17 | 2012-07-17 | 一种基于胎心电频谱图确定胎儿平均心率的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210245589.2A CN102715897B (zh) | 2012-07-17 | 2012-07-17 | 一种基于胎心电频谱图确定胎儿平均心率的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102715897A CN102715897A (zh) | 2012-10-10 |
| CN102715897B true CN102715897B (zh) | 2014-06-25 |
Family
ID=46941872
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210245589.2A Active CN102715897B (zh) | 2012-07-17 | 2012-07-17 | 一种基于胎心电频谱图确定胎儿平均心率的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102715897B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103767701A (zh) * | 2014-01-10 | 2014-05-07 | 中南民族大学 | 一种基于智能手机的可穿戴式胎儿心电监护系统及方法 |
| DE112015007217T5 (de) | 2015-12-23 | 2018-09-13 | Intel Corporation | HMM-basiertes adaptives Spektrogrammnachverfolgungsverfahren |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5666959A (en) * | 1995-08-30 | 1997-09-16 | British Technology Group Limited | Fetal heart rate monitoring |
| GB9923955D0 (en) * | 1999-10-08 | 1999-12-08 | Univ Nottingham | Fetal heart beat |
| CN101554325B (zh) * | 2009-05-15 | 2011-04-13 | 南京大学 | 胎儿心电图的分离提取方法 |
| CN102068249B (zh) * | 2009-11-23 | 2013-05-01 | 财团法人工业技术研究院 | 母体胎儿监视装置与方法 |
| CN102160787B (zh) * | 2011-05-31 | 2013-06-26 | 广东工业大学 | 一种基于时频变换的胎儿心电盲提取方法 |
-
2012
- 2012-07-17 CN CN201210245589.2A patent/CN102715897B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN102715897A (zh) | 2012-10-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102178522B (zh) | 一种母亲及胎儿心电信号qrs波中r波的检测定位方法 | |
| Jezewski et al. | A novel technique for fetal heart rate estimation from Doppler ultrasound signal | |
| CN102440768B (zh) | 脉搏波形特征点提取方法 | |
| CN101972145B (zh) | 一种基于源信号时域相对稀疏性的胎儿心电盲分离方法 | |
| JP2013532570A (ja) | 胎児モニタリング中に超音波信号の不明確性を削減する方法 | |
| US11253204B2 (en) | Method for assessing electrocardiogram signal quality | |
| Desai et al. | A real-time fetal ECG feature extraction using multiscale discrete wavelet transform | |
| Al-Angari et al. | A hybrid EMD-kurtosis method for estimating fetal heart rate from continuous Doppler signals | |
| CN103845079A (zh) | 一种基于盲分离的多普勒胎心音瞬时心率的检测方法 | |
| CN104473631A (zh) | 一种基于非负盲分离胎儿心电瞬时心率识别方法及系统 | |
| CN103006207A (zh) | 基于形态学的t波交替散点图法的心电信号分析方法 | |
| CN101627907B (zh) | 心电信号分析处理装置及方法 | |
| US7340289B2 (en) | Biomagnetic field measuring apparatus | |
| CN111358464A (zh) | 一种针对卧床病人的非接触式生命体征监测方法 | |
| CN104305992A (zh) | 一种交互式胎儿心电快速自动提取方法 | |
| CN106137185A (zh) | 一种基于模板构造小波的癫痫特征波检测方法 | |
| Ghaffari et al. | Robust fetal QRS detection from noninvasive abdominal electrocardiogram based on channel selection and simultaneous multichannel processing | |
| CN102715897B (zh) | 一种基于胎心电频谱图确定胎儿平均心率的方法 | |
| Desai et al. | A comparison and quantification of fetal heart rate variability using Doppler ultrasound and direct electrocardiography acquisition techniques | |
| CN207996202U (zh) | 孕妇体检用的胎心监护仪 | |
| CN2675047Y (zh) | 非线性胎儿心率检测仪 | |
| Lee et al. | Advances in fetal heart rate monitoring using smart phones | |
| CN201481411U (zh) | 心电信号分析处理装置 | |
| CN201743707U (zh) | 胎心率和胎心波形同时显示的胎儿心率仪 | |
| Jamaluddin et al. | Wavelet analysis on FECG detection using two electrodes system device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |