专利名称：一种胎动自动识别方法及其装置的制作方法
技术领域：
本发明涉及生物医学信号处理领域，具体的说是涉及一种胎动自动识别方法以及 实现该方法的装置。
背景技术：
胎动是胎儿安危的重要指标，一般从孕观周开始，胎动检测就具有了临床意义。 胎动正常表示着胎儿宫内良好，胎动的减少或消失是胎儿缺氧或死亡的标志，胎动急剧提 示胎儿急性缺氧并可能会死亡。在妊高征、过期妊娠、胎儿IUGR等高危妊娠时，在胎盘功能 影响胎心率之前即可引起胎动减少，一般胎动消失到胎心消失可长达12-48小时。而妊娠 期间，伴随着胎动发生的胎心加速是胎儿健康的表现，尤其在第一产程中，伴随胎动的加速 的出现是胎儿储备良好的重要指标之一。基于超声多普勒技术的测量方法是目前常用的胎动检测手段，超声探头发射的具 有一定频率的超声波束，经过母体腹壁遇到胎儿的运动部位后，形成反射波束，此波束被超 声探头接收，形成Doppler频移信号，通过算法识别出胎动信号，得到胎动的个数以及持续 时间。而实际的Doppler频移信号除了包含胎动信息外，还包括胎儿心跳、母体呼吸以 及母体运动等信息。胎动信号与胎心信号以及其他噪声信号比较，具有频率低、非周期性、 幅度高的特点。而不同孕妇，胎动信号的强度不同，并且在监护过程中，胎心位置改变也会 引起信号强度变化，现有的胎动识别方法中阈值不可调，会造成信号强度变化后无法检测 出胎动或无胎动误检为胎动的情况。目前在临床实时检测中最常用的是阈值检测法，对信号进行简单的滤波后，进行 阈值判断。信号中超过阈值则认为是胎动。但这种方法中阈值多为临床经验值，在临床应 用中无法根据不同孕妇的信号特点来调节，造成无法准确的检测出胎动信息。国内已经公 开一篇关于胎动算法的专利，专利号为CN101194826A，该专利中采用的方法是小波变化与 BP神经网络，这种方法计算的数据量较大，不能满足实时检测的要求。而国外相关检测算法 多基于互相关以及位移检测等理论，这些方法计算过程复杂，不适合临床应用。
发明内容
本发明的目的是为了克服现有的胎动识别方法可调性差，提出的一种自动识别胎 动的方法，其采用自动阈值调节算法，有效检测出胎动的次数以及单次胎动的持续时间。为了实现上述目的，本发明采用以下的技术方案1、一种胎动自动识别方法，其特征在于包括以下步骤A、对包括胎动信号的所有信号进行数据采集；B、将采集后的信号进行滤波预处理，分离出胎动信号；C、根据胎动信号的特点对信号进行实时分析，得到对应的实时阈值；D、根据得到的阈值对分离出的胎动信号进行识别，得到胎动的次数以及持续时间。所述的数据采集是将超声探头接收到的Doppler频移信号经过硬件滤波处理 后分成胎动信号和胎心信号；并对这两路信号进行采样，得到胎动信号m(n)和胎心信号 h(n)。步骤B中的分离胎动信号是采用差分滤波器对胎动信号m(n)进行处理，得到序列 D(η)；并对序列D(n)采用积分滤波器进行处理，得到序列M(η)。步骤C当中还包括a、当所述胎心信号h(n)质量良好即符合倍频关系时，抽取此时所述序列M(n)中 的相应部分，得到新序列Ml (η)，并对Ml (η)进行阈值分析。仏对机⑷加和求平均值，以平均值为参考得到该段时间内信号的基线值B(i)；C、重复步骤a、b，连续计算多次基线值B(i)，对其进行加权平均得到新的基线值 B ；d、以B为阈值对Ml (η)进行分析，统计大于阈值的数据个数M ；如果M > Ni，则认为阈值过低，增大基线值B(i)，重复步骤c、d，重新计算加权平 均值，直到M满足使超过基线值B的噪声信号能量等于A ；如果M < Ni，说明阈值过高，减小基线值B (i)，重复步骤c、d，重新计算加权平均 值；直到M满足使超过基线值B的噪声信号能量小于A ；Nl和A均为临床经验值；满足上述 条件的基线值B即为阈值T ；e、重复步骤a、b、c、d、得到新的阈值T(i)，比较T(i-l)与T(i)的大小，如果T (i) 大，则认为基线值改变，更改阈值为新的阈值T (i)，否则保持不变。步骤D还包括以下步骤根据阈值T对步骤B处理后得到的序列M(n)进行处理；计算连续大于阈值的数据 个数Ll，如果Ll大于N，则认为有胎动，且胎动开始，记录此时定时器的数值Cl ；如果有胎动开始计算连续低于阈值的数据个数L2，如果L2小于M,则认为为同一 次胎动，反之，则认为本次胎动结束，记录此时的定时器的数值C2 ；C2与Cl的差值即为本次 胎动的持续时间，其中N和M均为临床经验值。一种自动识别胎动信息的装置，其特征在于该装置主要包括信号采集模块、预处 理模块、阈值分析模块、胎动识别模块以及显示打印模块；信号采集模块将采集到的信号输 出至预处理模块处理，并经阈值分析模块进行阈值分析；再输送到胎动识别模块进行识别 后显示或打印。所述的信号采集模块包括一超声探头，其接收超声探头采集到的Doppler频移信 号滤波处理后分为两路信号，其中一路为胎动信号，另一路为胎心信号。所述的预处理模块包括差分处理单元与积分处理单元，其用于信号的去噪与预滤 波，得到胎动信号m(n)及胎心信号h(n)。阈值分析模块分为胎心倍频分析单元、阈值计算单元以及阈值质量判断单元；每 隔一定时间通过胎心倍频分析单元对胎心信号进行倍频分析；如果符合倍频，将相应的预 处理后的胎动信号输入阈值计算单元进行基线值计算；求取此段胎动信号Ml (η)的加和平 均值为参考得到基线值；连续计算多次时间的基线值，求加权平均值；阈值质量判断单元 以此平均值为阈值对胎心信号满足倍频时的胎动信号进行分析；如果超过阈值的数据个数> Ni，则认为阈值过低，调节基线值重新求加权平均值，再重新对Ml (η)进行判断；直到使 超过基线值B的噪声信号能量等于A ；如果超过阈值的数据< Ni，则认为阈值过高，调节基 线值，重新计算加权平均值，再对Ml (η)进行判断，直到使超过基线值B的噪声信号能量小 于A ；选取满足上述条件的加权平均值作为胎动识别的阈值。胎动识别模块主要负责采用阈值分析模块得到的阈值，对预处理后的信号进行胎 动识别，得到胎动的次数及每次胎动的持续时间；如果连续超过阈值的数据的个数大于N 则认为是一次胎动，如果连续低于阈值的数据的个数小于Μ，则认为是同一次胎动；通过计 算胎动开始及结束时的定时器的数值可以得到单次胎动的持续时间。采用上述方案，可以有效地滤除胎心跳、母体呼吸等噪声，提取胎动信息；自动阈 值调节算法，可以针对不同孕妇进行阈值调节，并且可以随监护过程中信号强度变化进行 阈值调节，避免了固定阈值引起的有胎动检测不到或误检的情况；可以准确得识别出胎动 信息，计算出监护过程中胎动的次数以及单次胎动的持续时间。
图1是胎动自动识别方法的流程图2是自动阈值调节算法的流程图3是采集的胎动信号示意图4为信号平滑差分处理后的效果图5为信号移动积分处理后的效果图6为胎心信号符合倍频时的处理后的胎动信号；
图7为通过阈值分析后得到的自动调节的阈值效果图
图8为用计算出来的阈值进行胎动识别的效果图9为本发明胎动自动识别装置的控制框图。
具体实施例方式下面结合附图及具体实施例对本发明做进一步的详细说明如图1的一种胎动自动识别方法，主要是处理超声Doppler频移信号中的胎动信 号，超声探头接收到Doppler频移信号中包括胎心信号以及胎动信号，还有母体呼吸、母体 运动等噪声信号，经过前期硬件处理滤波处理后得到胎动信号及胎心信号；其主要包括以 下步骤Α、对需要处理的信号进行数据采集；B、将采集后的信号采用数字滤波的方法进行预处理，分离出胎动信号；C、根据胎动信号的特点对信号进行实时分析，得到对应的实时阈值；D、根据得到的阈值对分离出的胎动信号进行识别，得到胎动的次数以及持续时 间。本发明中超声探头接收到的Doppler频移信号经过硬件滤波处理后分为两路信 号，一路为胎动信号，一路为胎心信号。对这两路信号进行采样，得到胎动信号m(n)及胎心 信号h(n)。如图3和图4、图5，对采集后的胎动信号m(n)进行预处理。采用数字滤波的方法，
6分离出胎动信号。其具体方法是，第一步采用差分滤波器对胎动信号m(n)进行处理，得到序列 D (η)。第二步采用积分滤波器对序列D(η)进行处理得到新的序列Μ(η)。如图2，根据胎动信号的特点对信号进行分析，得到阈值；并根据信号强度的变化 实时更改阈值，从而实现阈值的自动调节，其由以下几个步骤完成；a、胎心信号为Doppler频移信号，当胎心信号质量良好时，即胎心信号符合倍频 关系，一般没有胎动，抽取此时序列M (η)中的相应部分，得到新序列Ml (η)，并对Ml (η)进行 阈值分析。仏对机⑷加和求平均值，以平均值为参考得到该段时间内信号的基线值B (1)；C、重复步a、b，连续计算多次，基线值B (i)，对其进行加权平均得到新的基线值B;d、以B为阈值对Ml (η)进行分析，统计大于阈值的数据个数M ；如果M > Ni，则认为阈值过低，增大基线值B(i)，重复步骤c、d，重新计算加权平 均值，直到M满足使超过基线值B的噪声信号能量等于A ；如果M < Nl，说明阈值过高，减小基线值增大基线值B⑴，重复步骤c、d，重新计 算加权平均值。直到M满足使超过基线值B的噪声信号能量小于A。m和A为临床实验得 到的数据；满足上述条件的基线值B即为阈值T。e、重复上述步骤a、b、c、d、得到新的阈值T(i)，比较T(i-l)与T(i)的大小，如果 T (i)大，则认为基线值改变，更改阈值为新的阈值T (i)，否则保持。其实现方式是，当胎心信号符合倍频时，以此时步骤B处理后的信号M(n)来计算 阈值，可以有效的滤除胎心信号、母体呼吸等噪声信号。每隔一定时间对胎心信号h(η)进 行倍频分析，主要是通过自相关计算、包络分析、峰值检测处理后判断峰值是否符合倍频。 如果满足倍频关系，对步骤B中得到的相应的信号Ml (η)进行分析，否则不进行处理。图6所示，胎心信号质量良好时，没有胎动，此的信号Ml (η)只包括处理后的胎心 跳等噪声信号，没有大的波动。用Ml (η)进行阈值分析，可以有效滤除噪声，识别出胎动信 号。对Ml (η)求积加和平均值，以该平均值为参考得到基线值B(i)。重复第一步、第二 步，计算多次基线值，求多次基线值的加权平均值B，得到基线值B (i)，对其进行加权平均 得到新的基线值B。以基线值B为阈值对Ml (η)进行分析，统计序列中超过阈值的数据个数Μ。如果M > Ni，则认为阈值过低，增大基线值B(i)，重复步骤c、d，重新计算加权平 均值，直到M满足使超过基线值B的噪声信号能量等于A。如果M < Nl，说明阈值过高，减小基线值增大基线值B⑴，重复步骤c、d，重新计 算加权平均值，直到M满足使超过基线值B的噪声信号能量小于A。Nl和A均为临床实验 值；满足上述条件的基线值B作为胎动识别的阈值T。重复上述步骤C，求得新的阈值T⑴，修改胎动识别阈值为新阈值T⑴，并于上次 计算的阈值T(i-1)进行比较，如果T(i)大，则说明基线值改变，修改胎动识别阈值为新阈 值T (i)，否则保持不变。如图7，在tl、t2信号的强度改变了，相应计算得到的阈值也随之改变。对得到的阈值T对B过程中分离出的胎动信号进行识别，得到胎动的个数以及持续时间。其实现方式是根据阈值T对步骤B处理后得到的序列M(n)进行分析。计算连续 大于阈值的数据个数Li，如果Ll大于N，则认为有胎动，且胎动开始，记录此时定时器的数 值Cl。如果有胎动开始计算连续低于阈值的数据个数L2，如果L2小于M,则认为同一次 胎动，反之，则认为本次胎动结束，记录此时的定时器的数值C2 ；C2与Cl的差值即为本次胎 动的持续时间。图8为采用计算得到的阈值处理后的效果图；超过阈值的数据个数和Ml大于N， 本次胎动开始。低于阈值的数据个数和M2、M3都小于M,认为是同一次胎动，直到M4大于 M，本次胎动结束。如图9所示的一种自动识别胎动信息的装置主要包括信号采集模块、预处理模 块、阈值分析模块、胎动识别模块以及显示打印模块；其中信号采集模块将采集到的信号输 出至预处理模块处理，并经阈值分析模块进行阈值分析；再输送到胎动识别模块进行识别 后显示或打印。首先信号采集模块中的超声探头接收到的Doppler频移信号经过硬件滤波处理 后分为两路信号，其中一路为胎动信号，一路为胎心信号，并通过信号采集模块进行采集得 到这两路信号。信号采集后进入预处理模块，预处理模块包括差分处理单元与积分处理单元，主 要用于对采集信号的去噪与预滤波，得到胎动信号m(n)及胎心信号h(n)。胎动信号m(n)及胎心信号h(n)进入阈值分析模块，阈值分析模块分为胎心质量 分析单元、阈值计算单元以及阈值质量判断单元。每隔一定时间通过胎心质量分析单元对 胎心信号进行质量分析。如果质量良好，将相应的预处理后的胎动信号输入阈值计算单元 进行基线值计算。求取此段胎动信号Ml (η)的加和平均值为参考得到基线值，连续计算多 次时间的基线值，求加权平均值。在阈值质量判断单元以此平均值为阈值对胎心信号质量 良好时的胎动信号进行分析。如果超过阈值的数据个数大于W，则认为阈值过低，调节基线 值重新求加权平均值，再重新对Ml (η)进行判断。直到M满足使超过基线值B的噪声信号 能量等于Α。如果超过阈值的数据小于Ni，则认为阈值过高，调节基线值，重新计算加权平 均值，再对Ml (η)进行判断。直到M满足使超过基线值B的噪声信号能量小于Α。选取满足 上述条件的加权平均值作为胎动识别的阈值。胎动识别模块主要负责采用阈值分析模块得到的阈值，并对预处理后的信号进行 胎动识别，得到胎动的次数及每次胎动的持续时间。如果连续超过阈值的数据的个数大于 N则认为是一次胎动，如果连续低于阈值的数据的个数小于Μ,则认为是同一次胎动。通过 计算胎动开始及结束时的定时器的数值可以得到单次胎动的持续时间；得到胎动的次数及 单次胎动的持续时间就送入显示打印模块，对胎动次数及持续时间进行记录或显示。以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定 本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在 不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的 保护范围。
权利要求
1.一种胎动自动识别方法，其特征在于包括以下步骤A、对包括胎动信号的所有信号进行数据采集；B、将采集后的信号进行滤波预处理，分离出胎动信号；C、根据胎动信号的特点对信号进行实时分析，得到对应的实时阈值；D、根据得到的阈值对分离出的胎动信号进行识别，得到胎动的次数以及持续时间。
2.根据权利要求1所述的胎动自动识别方法，其特征在于所述的数据采集是将超声 探头接收到的Doppler频移信号经过硬件滤波处理后分成胎动信号和胎心信号；并对这两 路信号进行采样，得到胎动信号m(n)和胎心信号h(n)。
3.根据权利要求2所述的胎动自动识别方法，其特征在于步骤B中的分离胎动信号 是采用差分滤波器对胎动信号m(n)进行处理，得到序列D(n)；并对序列D(n)采用积分滤 波器进行处理，得到序列M(η)。
4.根据权利要求3所述的胎动自动识别方法，其特征在于步骤C当中还包括a、当所述胎心信号h(n)质量良好即符合倍频关系时，抽取此时所述序列M(n)中的相 应部分，得到新序列Ml (η)，并对Ml (η)进行阈值分析；bjiMl(n)加和求平均值，以平均值为参考得到该段时间内信号的基线值B (i)；c、重复步骤a、b，连续计算多次基线值B(i)，对其进行加权平均得到新的基线值B ；d、以B为阈值对Ml(η)进行分析，统计大于阈值的数据个数M ；如果M > Ni，则认为阈值过低，增大基线值B (i)，重复步骤c、d，重新计算加权平均值， 直到M满足使超过基线值B的噪声信号能量等于A ；如果M < Nl，说明阈值过高，减小基线值B (i)，重复步骤c、d，重新计算加权平均值；直 到M满足使超过基线值B的噪声信号能量小于A ；m和A均为临床经验值；满足上述条件的 基线值B即为阈值T ；e、重复步骤a、b、c、d、得到新的阈值T(i)，比较T(i-l)与T(i)的大小，如果T(i)大， 则认为基线值改变，更改阈值为新的阈值T (i)，否则保持不变。
5.根据权利要求1所述的胎动自动识别方法，其特征在于步骤D还包括以下步骤根据阈值T对步骤B处理后得到的序列M(n)进行处理；计算连续大于阈值的数据个数Li，如果Ll大于N，则认为有胎动，且胎动开始，记录此时定时器的数值Cl ；如果有胎动开始计算连续低于阈值的数据个数L2，如果L2小于M,则认为为同一次胎 动，反之，则认为本次胎动结束，记录此时的定时器的数值C2 ；C2与Cl的差值即为本次胎动 的持续时间，其中N和M均为临床经验值。
6.一种自动识别胎动信息的装置，其特征在于该装置主要包括信号采集模块、预处 理模块、阈值分析模块、胎动识别模块以及显示打印模块；信号采集模块将采集到的信号输 出至预处理模块处理，并经阈值分析模块进行阈值分析；再输送到胎动识别模块进行识别 后显示或打印。
7.根据权利要求6所述的自动识别胎动信息的装置，其特征在于所述的信号采集模 块包括一超声探头，其接收超声探头采集到的Doppler频移信号滤波处理后分为两路信 号，其中一路为胎动信号，另一路为胎心信号。
8.根据权利要求6所述的自动识别胎动信息的装置，其特征在于所述的预处理模块 包括差分处理单元与积分处理单元，其用于信号的去噪与预滤波，得到胎动信号m(n)及胎心信号h(n)。
9.根据权利要求6所述的自动识别胎动信息的装置，其特征在于阈值分析模块分为 胎心倍频分析单元、阈值计算单元以及阈值质量判断单元；每隔一定时间通过胎心倍频分 析单元对胎心信号进行倍频分析；如果符合倍频，将相应的预处理后的胎动信号输入阈值 计算单元进行基线值计算；求取此段胎动信号Ml (η)的加和平均值为参考得到基线值；连 续计算多次时间的基线值，求加权平均值；阈值质量判断单元以此平均值为阈值对胎心信 号满足倍频时的胎动信号进行分析；如果超过阈值的数据个数> Ni，则认为阈值过低，调 节基线值重新求加权平均值，再重新对Ml (η)进行判断；直到使超过基线值B的噪声信号能 量等于A ；如果超过阈值的数据< Ni，则认为阈值过高，调节基线值，重新计算加权平均值， 再对Ml (η)进行判断，直到使超过基线值B的噪声信号能量小于A ；选取满足上述条件的加 权平均值作为胎动识别的阈值。
10.根据权利要求6所述的自动识别胎动信息的装置，其特征在于胎动识别模块主要 负责采用阈值分析模块得到的阈值，对预处理后的信号进行胎动识别，得到胎动的次数及 每次胎动的持续时间；如果连续超过阈值的数据的个数大于N则认为是一次胎动，如果连 续低于阈值的数据的个数小于Μ，则认为是同一次胎动；通过计算胎动开始及结束时的定 时器的数值可以得到单次胎动的持续时间。
全文摘要
本发明公开了一种胎动自动识别方法及其装置，包括以下几个步骤，对包括胎动信号的所有信号进行数据采集；将采集后的信号进行滤波预处理，分离出胎动信号；根据胎动信号的特点对信号进行实时分析，得到对应的实时阈值；根据得到的阈值对分离出的胎动信号进行识别，得到胎动的次数以及持续时间；采用上述方案可以有效地滤除胎心跳、母体呼吸等噪声，提取胎动信息；自动阈值调节算法，可以针对不同孕妇进行阈值调节，并且可以随监护过程中信号强度变化进行阈值调节，避免了固定阈值引起的有胎动检测不到或误检的情况；可以准确得识别出胎动信息，计算出监护过程中胎动的次数以及单次胎动的持续时间。
文档编号A61B5/00GK102090883SQ20101061974
公开日2011年6月15日 申请日期2010年12月30日 优先权日2010年12月30日
发明者刘冰, 孙艳, 陈吴笋, 陈德伟 申请人:深圳市理邦精密仪器股份有限公司