专利名称：一种基于胎心电频谱图确定胎儿平均心率的方法
技术领域：
本申请涉及一种基于胎儿心电频谱图确定胎儿平均心率的方法。胎儿在整个孕期处在非常艰难的生存环境中，胎心率是胎儿健康状况的一个重要指标。多普勒超声方法检测胎心率时，探头要对中胎儿心脏，在实际操作中多有不便，尤其在早期胎心比较小的时候。在孕妇腹壁体表安放电极采集胎儿心电信号，据此确定胎心率的方法，不要求电极对中胎心，操作较为方便，但也存在信号信噪比低、稳定地获取胎心率不容易的问题。文献中介绍的基于胎儿心电信号确定平均胎心率的方法主要是时域的方法(I)首先从孕妇腹壁检测的胎儿心电信号中分离去主要的干扰成分——母体心电成分，并尽量多地去掉基线漂移等其他干扰成分；(2)去除干扰后的保留信号波形，被认为主要是胎儿心电成分，基于差分阈值法或小波奇异性检测等方法，进行R波波峰检测，利用检测的各R 波波峰，计算瞬时胎心率，进而计算胎儿平均心率。这样的时域方法，实现比较复杂；有时，保留信号会残留有较大母体心电分量，导致胎儿R波波峰检测失败，方法失效。本申请涉及的是基于胎儿心电信号频谱图确定胎儿平均心率的问题，不同于以上所述的时域方法。
背景技术：
利用胎儿心电信号频谱确定胎儿平均心率的方法，其步骤可概括为(I)读取从孕妇腹壁体表检测的胎儿心电信号，直接制作其幅值谱或功率谱，或者对其进行诸如去干扰、降噪、变换等处理后再制作其幅值谱或功率谱；(2)在幅值谱或功率谱上，从左到右扫寻谱峰，每扫到一个谱峰，检查该谱峰频率的倍频处是否也存在谱峰，如果存在，并且该谱峰频率又不对应母体平均心率，则判定该谱峰对应频率为平均胎心率，否则为伪差谱峰。在频谱上确定胎儿平均心率方法的原理和在频谱上确定母体平均心率方法的原理相同，都是利用谱峰的倍频关系。但是，确定母体平均心率比较容易，确定胎儿平均心率比较困难，特别是要稳定可靠地确定胎儿平均心率很困难，因为胎儿分量往往很小、在频谱上胎儿分量对应的一组谱峰往往不是很突出，存在很多伪差谱峰。相对而言，从频谱上确定胎儿平均心率的比较放心可靠的方法，不是让计算机程序从左到右自动扫寻谱峰并自动判定一组谱峰频率是否有倍频关系，而是把幅值谱或功率谱图画出来、展现在人机界面上，由操作者通过人机交互来确定胎儿平均心率。比较容易想到的模式是用眼睛和大脑来扫寻频谱图上的各谱峰，目测出其频率成倍数关系的一组谱峰，若其肯定不对应母体分量，则选中这组谱峰中最左边的一个谱峰，将其对应频率确定为胎儿平均心率。胎儿心电技术发展严重滞后于需求，以上提到的比较容易想到的在频谱上通过人机交互确定胎儿平均心率的方法，虽然简单，也未见文献提到。另一方面，该简单方法一旦实施也会发现新的问题由于很多时候胎儿心电分量很弱，在频谱上对应胎儿的谱峰不是很突出，频谱上会存在很多伪差谱峰；此时，单纯用眼睛和大脑来扫寻频谱图上的各谱峰，并且目测出一组谱峰的频率是否存在倍数关系，操作者会很累，须仔细分辨和心算，操作起来很慢，效率低。参考文献[ I] Karvouni s, E. C. , C. Papaloukas , et al. " Fetal heart rateextraction from composite maternal ECG using complex continuous wavelettransform." Computers in Cardiology，2004，Vol 31:737-740.[2]Karvounis, E. C.，M. G. Tsipouras, et al. " An automated methodologyfor fetal heart rate extraction from the abdominal electrocardiogram. " IEEETransactions on Information Technolo gy in Biomedicine，2007，11(6) :628-638.[3] S. Kadambe，R. Murray, and G. Faye Boudreaux-Bartels, ^WaveletTransform-Based QRS Complex Detector，，，IEEE Trans, on Biomed. Eng.，vol. 46, no. 7，pp. 838-848，July. 1999.
发明内容
发明目的。提出一种基于胎儿心电频谱图确定胎儿平均心率的方法，改进容易想到的简单人机交互模式，以提高胎儿平均心率的估计效率，减轻操作者的劳动强度。技术方案。本发明涉及一种基于胎儿心电频谱图确定胎儿平均心率的方法，它包括两个步骤——⑴读取从孕妇腹壁检测的胎儿心电信号、制作幅值谱或功率谱图，⑵在制作的幅值谱或功率谱图上通过人机交互确定胎儿平均心率，其特征是，在第(2)步人机交互中，在幅值谱或功率谱图上设置一组游标(须两个以上)，它们的频率坐标始终成倍数关系，操作者拖动这组游标中其中一个时，组内所有游标同时游动，并且各游标频率坐标保持着始终不变的倍数关系。以上方法第(I)步骤中，读取孕妇腹壁检测胎心电信号后，可直接制作幅值谱或功率谱，也可以对其进行诸如去干扰、降噪、变换等各种处理后再制作其幅值谱或功率谱，此处不构成本方案的技术特征。以上方法第(2)步骤中，由于频谱图上设置了一组始终具有倍频关系的游标，当操作者将最左边的游标移至一个谱峰时，则其他游标正好位于该谱峰频率的倍频处，非常容易观察倍频处是否也有谱峰；当发现一组谱峰的频率正好具有倍数关系，又不对应母体心电分量时，可将最左边谱峰的频率确定为胎儿平均心率。人的眼睛和大脑擅长发现识别一个频谱图上的谱峰，但很不擅长测算谱峰频率之间是否存在倍数关系。前面提到的容易想到的简单人机交互模式，需要用人的眼睛和大脑来测算谱峰频率之间是否存在倍数关系，这是方法低效的根源；但是，本技术方案采用了新的交互模式，不再需要人的眼睛和大脑来测算谱峰频率之间是否是倍数关系，把人工测算一组谱峰频率是否成倍数关系的问题变成了频率坐标成倍数关系的一组游标是否都经过谱峰的认定问题，简言之，把一个计算问题变成了一个眼睛和大脑擅长的单纯的识别问题。本技术方案反映了人机交互确定平均心率的本质需要，是客观规律的自觉利用。有益效果。
目前，尚未见文献报道基于胎儿心电频谱图人机交互确定胎儿平均心率的方法；当然，基于在胎儿心电频谱图上设置一组有倍频关系的游标来人机交互确定胎儿平均心率的方案更是未见报道。因此，本技术方案具有新颖性。经大量测试，本技术方案在胎儿心电频谱图上设置一组有倍频关系的游标这一做法，与容易想到的简单人机交互方法相比，可提高胎儿平均心率的确定效率5倍以上。因此，提出的方法有实用性，有实质进展。
图1，本发明基于胎儿心电频谱图确定胎儿平均心率的步骤及特征。图2，一腹壁实测胎儿心电信号的波形图。
图3，幅值谱图及经过母体平均心率处谱峰的一组游标示意图。图4，幅值谱图及经过一伪差谱峰的一组游标示意图。图5，幅值谱图及经过胎儿平均心率处谱峰的一组游标示意图。
具体实施方式
(实施例)在医院对某一孕妇用日本光电标准12导心电图机——ECG-1350P将电极置于腹壁测了一组数据，采样率500Hz。首先，读取从孕妇腹壁检测的一个导联信号的胎儿心电信号，选择以0.5Hz到IOOHz为通带的带通滤波器去除基线漂移作为预处理，去除基线漂移后的信号波形图见图2所示，图中横轴为时间轴，时间刻度0、1、…、8的单位为秒；然后制作信号幅值谱，见图3、图4、图5所示。在图3、图4、图5幅值谱中，横坐标均为频率轴，频率刻度0、0. 5、I、I. 5、…、10的单位为Hz。在制作的幅值谱图上设置了一组游标(此处为三个)，见图3、图4、图5中的三根竖线、它们的频率坐标始终成倍数关系。操作者从幅值谱图左边开始向右拖动游标，首先将游标落到第一个谱峰及其倍频位置处，见图3所示，发现正好对应3个谱峰，则易断定此组谱峰中的第I个谱峰对应母体平均心率；继续向右移动游标，将游标落到第2个谱峰及其倍频位置处，见图4所示，发现后面的两个游标没有对上谱峰，可断定此谱峰为I个伪差谱峰；继续向右移动游标，将游标落到第3个谱峰及其倍频位置处，见图5所示，可发现后面的两个游标也对上了谱峰，说明该组游标所指谱峰是对应胎儿心电的，则可将此时左边第I个游标指向的频率确定为胎儿平均心率。
权利要求
1.本发明涉及一种基于胎儿心电频谱图确定胎儿平均心率的方法，它包括两个步骤一⑴读取从孕妇腹壁检测的胎儿心电信号、制作幅值谱或功率谱图，⑵在制作的幅值谱或功率谱图上通过人机交互确定胎儿平均心率，其特征是，在第(2)步人机交互中，在幅值谱或功率谱图上设置一组游标(须两个以上)，它们的频率坐标始终成倍数关系，操作者拖动这组游标中其中一个时，组内所有游标同时游动，并且各游标频率坐标保持着始终不变的倍数关系。
全文摘要
本发明涉及一种基于胎儿心电频谱图确定胎儿平均心率的方法，它包括两个步骤(1)读取从孕妇腹壁检测的胎儿心电信号、制作幅值谱或功率谱图，(2)在制作的幅值谱或功率谱图上通过人机交互确定胎儿平均心率，其特征是，在第(2)步人机交互中，在幅值谱或功率谱图上设置一组游标(须两个以上)，它们的频率坐标始终成倍数关系，操作者拖动这组游标中其中一个时，组内所有游标同时游动，并且各游标频率坐标总保持着始终不变的倍数关系。该方法将人工测算一组谱峰频率是否成倍数关系的问题变成了频率坐标成倍数关系的一组游标是否都经过谱峰的认定问题，可提高胎儿平均心率的确定效率5倍以上。
文档编号A61B5/0444GK102715897SQ201210245589
公开日2012年10月10日 申请日期2012年7月17日 优先权日2012年7月17日
发明者刘红星, 刘铁兵, 宁新宝, 张洁敏, 李萍, 汤黎明, 石立娟, 管群, 肇莹, 闫华文 申请人:南京大学