CN102048558A - 一种胎心率信号处理方法及其装置 - Google Patents

一种胎心率信号处理方法及其装置 Download PDF

Info

Publication number
CN102048558A
CN102048558A CN 201110030452 CN201110030452A CN102048558A CN 102048558 A CN102048558 A CN 102048558A CN 201110030452 CN201110030452 CN 201110030452 CN 201110030452 A CN201110030452 A CN 201110030452A CN 102048558 A CN102048558 A CN 102048558A
Authority
CN
China
Prior art keywords
signal
unit
frequency
input
heart rate
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
Application number
CN 201110030452
Other languages
English (en)
Other versions
CN102048558B (zh
Inventor
刘冰
陈吴笋
陈德伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Edan Instruments Inc
Original Assignee
Edan Instruments Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Edan Instruments Inc filed Critical Edan Instruments Inc
Priority to CN2011100304520A priority Critical patent/CN102048558B/zh
Publication of CN102048558A publication Critical patent/CN102048558A/zh
Application granted granted Critical
Publication of CN102048558B publication Critical patent/CN102048558B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

在本发明提供一种胎心率信号处理方法及其装置,将时域处理和频域处理结合起来,进而得到一种即能满足实时性要求,又有较好稳定性的胎心率处理方法,并将其处理方法应用在胎儿监护仪上,使得可以得到的中间监护数据更加准确。本发明能对监护过程中产生的一些特殊的中间监护数据进行校准,使得医生能根据上述经过校准得到的中间监护数据作为进行诊断的参照,极大的减少对胎心跳动错误判断的情况发生,便于医生的诊断。

Description

一种胎心率信号处理方法及其装置
技术领域
本发明涉及胎心信号处理领域,尤其涉及一种胎儿心率信号处理方法及装置。
背景技术
[0002] 胎心率检测在临床诊断中具有重要的意义,尤其是对围产期胎心率的监测,有助于提高胎儿分娩质量,降低胎儿畸形率和死亡率。胎心率检测大致有两种方式,一种是直接法,这种方法是利用螺旋电极勾到胎儿头皮获取心率信号,属于有创的,并且要等到孕妇破膜之后,宫口开到一定程度时才能使用,有一定的局限性;另一种是间接法,目前用的最多的是超声波多谱勒心率检测法,它是通过利用回声和多普勒原理,提取出超声多普勒的频移信号,从而得到胎心包络信号,由此求得胎心的一种方法具有无损伤性和操作简单等优点,得到了广泛的应用。
一个典型的基于超声多普勒获取心率的方法是利用对采集的包络信号进行自相关处理,然后基于自相关结果序列求得信号的周期,进而得到胎儿心率。
自相关函数的公式描述为:
在信息分析中,通常将自相关函数称之为自协方差方程。用来描述信息在不同时间τ的信息函数值的相关性。
自相关具有使能使周期信号得到加强而随机噪声被减弱的特性,在一定程度上可以减少干扰对信号的影响,理论上,自相关结果序列的第一主峰就对应信号周期。但是由于人体是一个复杂的运动环境,尤其是围产期胎儿本身的心跳就比较弱,再加上母体环境的干扰,超声回波信号成分比较复杂,干扰大。因此,采集的信号有时干扰较多,如果只依赖自相关结果序列求得信号周期,有时候很难得到正确的胎心周期,从而有可能导致胎心的减半或加倍现象,给医生的诊断造成困扰。
相对时域的处理方法来说,频域的处理方法实时性不强,但是抗干扰能力比较强,得到的结果要更稳定,常用的求取频域的方法是FFT和CZT。
理论上,FFT就可以得到比较精确的频率,但是由于FFT固有的频率分辨率与计算量之间的矛盾也限制了它的应用,即要提高FFT的频率分辨率,需要增加分析的采样点数,而这将大大增加计算量。因此,我们可以先用FFT对信号频谱进行粗略的估计,获取需要细化分析的窄带频谱范围,再在该窄带范围内进行频域细化分析,以获得满足精度要求的频谱。
CZT是一种典型的频域细化的方法,它可以在Z平面单位圆上取一个自定义的弧段,只在该弧段上进行序列Z变化的均匀取样,而且取样间隔也可以自由确定,如果所取弧段对应于待细化的窄带,则CZT就是窄带中个频率点处的频谱值。因此,CZT适合窄带高分辨率的计算。
从上面的论述可以看出,时域有实时性的优势,频域有稳定性的优势。
目前,现有的方法都是在时域对信号进行自相关处理,进而求得胎心率,虽然在实时性方面满足临床要求,但是在稳定性方面性能不佳,在临床中,当信号不好的情况下,时常会出现减半和加倍的现象。因此,在本专利中,我们将时域处理和频域处理结合起来,应用在胎儿监护仪上,进而得到一种即能满足实时性要求,又有较好稳定性的胎心率处理方法,它可以得到更加准确、可靠的胎儿心率,减少临床上胎心减半或加倍的现象,减少由此导致医生误诊而引发的医疗事故。
发明内容
本发明提供一种更准确地获取胎心率的方法及其实现装置,它能提供更准确、可靠的胎儿心率,特别是在信号比较差的情况下,对于胎心率的加倍和减半有很好的抑制作用。
本发明是通过以下技术方案实现的:一种胎心率信号处理方法,包括以下步骤:
Figure 689923DEST_PATH_IMAGE002
.采集胎心多普勒包络信号;
.对采集的胎心多普勒包络信号进行去噪、滤波处理,得到预处理信号;
Figure 670835DEST_PATH_IMAGE004
.所述预处理信号分为两路信号,一路经过频域分析,另一路经过时域分析;
Figure 84630DEST_PATH_IMAGE005
.对时域分析后确定的胎心信号的初步周期与根据频域分析求得的信号频率进行换算比较,根据比较结果,确认胎心信号的周期;
Figure 547972DEST_PATH_IMAGE006
.对得到的胎心信号的周期进行折算,求取胎儿心率。
进一步地,所述频域分析步骤和时域分析步骤不分先后且可同时进行。
进一步地,所述频域分析步骤为:
Figure 954683DEST_PATH_IMAGE002
′.判断上述预处理信号是否丢失,若丢失,则返回步骤
Figure 187081DEST_PATH_IMAGE002
.重新采集胎心多普勒包络信号,若未丢失,则
Figure 771777DEST_PATH_IMAGE003
′.对预处理信号做FFT运算;
Figure 722416DEST_PATH_IMAGE004
′.对经FFT运算后的信号进行峰值检测,找出测定频率段内的测定极大值点,定位信号的频率范围;
Figure 667238DEST_PATH_IMAGE005
′.在上述定位的信号频率范围内对预处理信号进行CZT运算;
Figure 19722DEST_PATH_IMAGE006
′.对做完CZT运算后的频率数据进行求极值处理,找出极大值,求得信号频率。
进一步地,所述时域分析步骤为:
Figure 772390DEST_PATH_IMAGE007
′.对预处理信号进行自相关处理,得到自相关结果序列;
Figure 210324DEST_PATH_IMAGE008
′.对得到的自相关结果序列进行峰值检测,找到测定峰值;
Figure 958837DEST_PATH_IMAGE009
′.通过上述找到的测定峰值,确定胎心信号的初步周期。
进一步地,所述
Figure 165828DEST_PATH_IMAGE008
′.的步骤包括:
.找出自相关结果序列中所有的峰值;
Figure 17557DEST_PATH_IMAGE011
.通过阈值检测得到检测峰值。
进一步地,还包括:在步骤
Figure 304182DEST_PATH_IMAGE005
之后,当信号好时,经过处理直接得到信号周期;当信号不好时,则将所述自相关结果进行二次分析,得到修正后的信号周期,再将信号周期折算后与经
Figure 444308DEST_PATH_IMAGE006
′步骤得到的信号频率进行比较,根据比较结果,确认胎心信号的周期。
一种获取胎心率的装置,包括数据采集单元201、预处理单元202、信号丢失判断单元203、第一报警单元204、时域分析处理器、后处理单元208、周期计算单元216、心率计算单元217、心率异常判断单元220、第二报警单元221、频域分析处理器、参数更新单元214、控制单元215、显示单元218、打印单元219。
数据采集单元201用于采集胎心多普勒包络信号。
预处理单元202与数据采集单元201连接,用于对数据采集单元采集到的胎心多普勒包络信号进行去噪、滤波处理。
信号丢失判断单元203判断经过预处理的信号是否丢失,若信号丢失,则启动第一报警单元204。
将经预处理后的信号数据分为两路:一路输入频域分析处理器,进行频域分析;另一路输入时域分析处理器,进行时域分析。
将经所述频域分析器分析后的频率信息输入到参数更新单元,保证参数实时性;将更新的频率信息输入到控制单元215,由控制单元对后处理单元208中确认的后处理周期进行判断和控制,保证周期稳定性。
将经所述时域分析处理器时域分析并初步确定信号周期的信号输入到后处理单元208,根据参数更新单元214提供的参数经控制单元215的判断和控制确认后处理周期。
将经过后处理单元208处理后的信号输入到周期计算单元216,得到所需的确定周期。
将得到的确定周期输入到心率计算单元217,得到所需的胎儿心率。
将心率结果输入到显示单元218,完成心率的实时显示。
将显示的心率结果输入到打印单元219,完成心率曲线的实时打印。
心率异常判断单元220对心率是否异常进行判断,并根据不同情况触发第二报警单元221中对应的报警信号。
进一步地,所述时域分析处理器中包括自相关处理单元205、第一峰值检测单元206及初步确定周期单元207;所述经预处理的信号数据输入到自相关处理单元205,对预处理后的信号数据进行自相关处理,得到自相关结果序列;对所述自相关处理单元205得到的自相关结果序列输入到第一峰值检测单元206,查找所需要的测定峰值;将峰值检测处理后的信号输入到初步确定周期单元207,确定信号的初步周期。 
进一步地,所述频域分析处理器中包括FFT处理单元209、第二峰值检测单元210、范围定位单元211、CZT处理单元212及频率确定单元213;经预处理的信号数据输入到FFT处理单元209;信号经FFT处理单元209处理后,得到频域数据,将该频域数据输入到第二峰值检测单元210,找出测定频段内的测定极大值点;根据得到的测定极大值点,在范围定位单元211中得到信号频率的初步范围;将得到信号频率的初步范围的信号输入CZT处理单元212,根据信号频率的初步范围对信号做CZT变换;经CZT变换后的信号输入到频率确定单元213,提取出信号的频率信息。
进一步地,根据所述参数更新单元214提供的参数经控制单元215的判断和控制确认后处理周期;当信号好时,则直接将经过后处理单元208处理后的信号输入到周期计算单元216,得到所需的确定周期;当信号不好时,则对所述自相关结果进行二次分析,得到修正后的后处理周期,然后再将经过后处理单元208处理后的信号输入到周期计算单元216,得到所需的确定周期。
本发明的优点及有益效果。
时域的自相关处理技术具有较好的实时性,但在稳定性方面不如频域得到的结果,而频域的FFT与CZT具有较好的稳定性,但在实时性方面又不如时域的自相关技术,在本发明中,我们将时域处理和频域处理结合起来,进而得到一种即能满足实时性要求,又有较好稳定性的胎心率处理方法,并将其处理方法应用在胎儿监护仪上,使得可以得到的中间监护数据更加准确。本发明能对监护过程中产生的一些特殊的中间监护数据进行校准,使得医生能根据上述经过校准得到的中间监护数据作为进行诊断的参照,极大的减少对胎心跳动错误判断的情况发生,便于医生的诊断。
附图说明
图1.采集的胎心多普勒包络信号图。
图2.自相关结果序列图。
图3.FFT结果,得到粗略的频率信息图。
图4.CZT结果,得到精确的频率信息图。
图5.求取胎心率的信号处理方法流程图。
图6.本发明实施例中胎儿监护仪框图。
具体实施方式
本发明提供一种胎心率信号处理方法及其装置,其中胎心率信号处理方法包括以下步骤:
Figure 525396DEST_PATH_IMAGE002
.采集胎心多普勒包络信号;
Figure 937923DEST_PATH_IMAGE003
.对采集的胎心多普勒包络信号进行去噪、滤波处理,得到预处理信号;
Figure 778971DEST_PATH_IMAGE004
.所述预处理信号分为两路信号,一路经过频域分析,另一路经过时域分析;
Figure 960554DEST_PATH_IMAGE005
.对时域分析后确定的胎心信号的初步周期与根据频域分析求得的信号频率进行换算比较,根据比较结果,确认胎心信号的周期;
.对得到的胎心信号的周期进行折算,求取胎儿心率。
进一步地,所述频域分析步骤和时域分析步骤不分先后且可同时进行。
进一步地,所述频域分析步骤为:
′.判断上述预处理信号是否丢失,若丢失,则返回步骤
Figure 488596DEST_PATH_IMAGE002
.重新采集胎心多普勒包络信号,若未丢失,则
Figure 852582DEST_PATH_IMAGE003
′.对预处理信号做FFT运算;FFT运算,即快速傅氏变换运算,是离散傅氏变换的快速算法,它是根据离散傅氏变换的奇、偶、虚、实等特性,对离散傅立叶变换的算法进行改进获得的。它对傅氏变换的理论并没有新的FFT算法图发现,但是对于在计算机系统或者说数字系统中应用离散傅立叶变换,可以说是进了一大步。
设x(n)为N项的复数序列,由DFT变换,任一X(m)的计算都需要N次复数乘法和N-1次复数加法,而一次复数乘法等于四次实数乘法和两次实数加法,一次复数加法等于两次实数加法,即使把一次复数乘法和一次复数加法定义成一次"运算"(四次实数乘法和四次实数加法),那么求出N项复数序列的X(m),即N点DFT变换大约就需要N2次运算。当N=1024点甚至更多的时候,需要N2=1048576次运算,在FFT中,利用WN的周期性和对称性,把一个N项序列(设N=2k,k为正整数),分为两个N/2项的子序列,每个N/2点DFT变换需要(N/2)2次运算,再用N次运算把两个N/2点的DFT变换组合成一个N点的DFT变换。这样变换以后,总的运算次数就变成N 2(N/2)2=N N2/2。继续上面的例子,N=1024时,总的运算次数就变成了525312次,节省了大约50%的运算量。而如果我们将这种"一分为二"的思想不断进行下去,直到分成两两一组的DFT运算单元,那么N点的DFT变换就只需要Nlog2N次的运算,N在1024点时,运算量仅有10240次,是先前的直接算法的1%,点数越多,运算量的节约就越大,这就是FFT的优越性;
Figure 291784DEST_PATH_IMAGE004
′.对经FFT运算后的信号进行峰值检测,找出测定频率段内的测定极大值点,定位信号的频率范围;
Figure 475641DEST_PATH_IMAGE005
′.在上述定位的信号频率范围内对预处理信号进行CZT(线性调频Z变换)运算;线性调频Z变换(chirp Z transform,CZT)是一种有效的频率细化方法。采用DFT或FFT 计算得到的信号N 点频谱实际上是Z 平面单位圆上N 点等间隔的采样结果;而采用CZT 计算得到的频谱则是Z平面螺旋周线上的等间隔采样结果。(王楠,肖先勇,曾颂崎,顾威.基于插值线性调频Z 变换的间谐波分析方法.电 网 技 术.2007年9月.第31 卷 第18期);
Figure 658492DEST_PATH_IMAGE006
′.对做完CZT运算后的频率数据进行求极值处理,找出极大值,求得信号频率。
进一步地,所述时域分析步骤为:
Figure 876984DEST_PATH_IMAGE007
′.对预处理信号进行自相关处理,得到自相关结果序列;
Figure 674038DEST_PATH_IMAGE008
′.对得到的自相关结果序列进行峰值检测,找到测定峰值;
Figure 95923DEST_PATH_IMAGE009
′.通过上述找到的测定峰值,确定胎心信号的初步周期。
进一步地,所述′.的步骤包括:
.找出自相关结果序列中所有的峰值;
Figure 44791DEST_PATH_IMAGE011
.通过阈值检测得到检测峰值。
进一步地,还包括:在步骤
Figure 974480DEST_PATH_IMAGE005
之后,当信号好时,经过处理直接得到信号周期;当信号不好时,则将所述自相关结果进行二次分析,得到修正后的信号周期,再将信号周期折算后与经
Figure 748401DEST_PATH_IMAGE006
′步骤得到的信号频率进行比较,根据比较结果,确认胎心信号的周期。
一种获取胎心率的装置,包括数据采集单元201、预处理单元202、信号丢失判断单元203、第一报警单元204、时域分析处理器、后处理单元208、周期计算单元216、心率计算单元217、心率异常判断单元220、第二报警单元221、频域分析处理器、参数更新单元214、控制单元215、显示单元218、打印单元219。
所述时域分析处理器中包括自相关处理单元205、第一峰值检测单元206及初步确定周期单元207;所述频域分析处理器中包括FFT处理单元209、第二峰值检测单元210、范围定位单元211、CZT处理单元212及频率确定单元213。
为了更充分理解本发明的技术特征,下面结合附图对本发明的内容进行描述,以下的描述仅是示范性和解释性的,不应对本发明的保护范围有任何的限制作用。
一种胎心计算的信号处理方法,采集和处理的信号是超声回波的胎心包络信号
Figure 613589DEST_PATH_IMAGE012
请参考图5,具体实施的步骤如下:
采集胎心多普勒包络信号101。采集的胎心多普勒信号,如图1所示。
对采集的信号进行预处理102,因为要处理的信号是带限信号,对应的信号频率只能在一定的范围内,所以要对数据进行去噪、滤波等处理。
判断信号是否丢失103,,如果信号丢失,则不对信号进行处理,直接返回数据采集部分;如果信号没有丢失,则对信号进行后面的处理。
将处理后的信号分成两路同样信号,一路用于频域分析(104~108),一路用于时域分析(110~113),,这里先说明频域的处理步骤。
对预处理的信号做FFT运算104,得到信号频域方面的信息,需要注意的是,为了减少计算量,在这一环节不需要求得非常准确的频率信息,只需要粗略估计出信号所在频段即可。图3给出了信号做完FFT后的结果。
对做完FFT之后的信号进行峰值检测105,找出有效频段内的主要极大值点(目前公认的频率段是0.5~4HZ),在这个环节,有可能会出现两个或多个频率幅度相当的点,在这种情况下,要根据之前信号的频率段为参考,来对找到的极大值点进行筛选。从图3可以看出,有效频段内主要的极大值点在3HZ附近。
确定了主要的频率点后,就可以确定信号所在的频率段,即定位了信号的频率范围106,为了防止各方面误差,我们将信号的频率范围稍微加大一点,我们将其信号频率范围定为2~4hz。
对信号进行CZT运算107,CZT运算的范围就是上面定位的信号的频率范围2~4hz,并根据定位的频率范围的大小来适当调整运算的点数。由此就可以得到比较精确的信号频率。图4给出了针对定位后的信号频率范围做CZT得到的结果。
对CZT之后的频率数据进行求极值处理,找出局部最大值,求得信号频率108。这个频率将由于时域处理的后处理中,用于确定信号的周期。从图4可以看出,极大值点对应的频率是2.406HZ,即信号的频率。下面说明时域的处理措施。
将预处理后的信号进行自相关处理109,得到自相关结果序列。图2给出了得到自相关序列。
对得到的自相关结果序列进行峰值检测110,找出测定峰值,即为主要峰值,首先是找到全部的峰值,然后通过阈值检测,得到主要峰值。
通过找出的主要峰值,初步确定信号的周期111。在信号较好的情况下,初步确定的周期是基本正确的,但是在信号较差的情况下,确定的后期主要是通过后处理来求得。
将初步确定的周期和自相关结果进行后处理112,这个步骤主要是确认周期,当信号较好的时候,确认的周期可能就是初步确定的周期,当信号不好时,需要对前面的自相关结果进行二次分析,得到一个或几个较为可能的周期。然后再根据频域提供的频率参数来选择对应的周期。
经过后处理后,很容易就可以获得信号的确定周期113。
对得到的确定周期进行折算,得到所需的胎儿心率114,此处的折算公式为60/确定周期。
如图5所示,包括数据采集单元201、预处理单元202、信号丢失判断单元203、第一报警单元204、时域分析处理器、后处理单元208、周期计算单元216、心率计算单元217、心率异常判断单元220、第二报警单元221、频域分析处理器、参数更新单元214、控制单元215、显示单元218、打印单元219。其中所述时域分析处理器中包括自相关处理单元205、第一峰值检测单元206及初步确定周期单元207;所述频域分析处理器中包括FFT处理单元209、第二峰值检测单元210、范围定位单元211、CZT处理单元212及频率确定单元213。
整个装置的流程如下。
由数据采集单元201负责采集胎心多普勒包络信号。
数据采集完成后,将采集的包络数据输入到预处理单元202;预处理单元202主要负责对数据进行去噪、滤波等处理。
预处理后的信号进入信号丢失判断单元203,判断信号是否丢失,信号如果丢失,就直接返回数据采集部分,并启动第一报警单元204,第一报警单元204发出报警信号,提示医护人员,采取相应的处理措施。
如果信号没有丢失,那么将预处理单元202处理后的数据分成两路,一路输入到自相关处理单元205,完成数据的自相关处理,得到自相关结果序列;另一路输入到FFT处理单元209,对信号进行频域的分析。先介绍频域方面的处理。
将预处理单元202处理后的信号输入到FFT单元209中,对信号进行频域的分析。
信号经过FFT处理单元209后,得到频域数据,将其输入到第二峰值检测单元210,找出有效频段内的主要极大值点(目前公认的有效频率段是0.5~4HZ),即为检测极大值点,在这个环节,有可能会出现两个或多个频率幅度相当的点,在这种情况下,要根据之前信号的频率段为参考,来对找到的极大值点进行筛选。
根据得到的主要极大值点,在范围定位单元211中得到信号频率的大体范围。
将信号输入CZT处理单元212,并根据信号频率的大体范围,对信号做CZT变换,得到比较精确的频率信息。
CZT变换后的信号输入到频率确定单元213,提取出信号的频率信息。
将频率确定单元213中得到的频率信息输入到参数更新单元214,保证参数的实时性。
将更新的参数输入到控制单元215,由控制单元对后处理中周期的确认进行判断和控制,保证周期的稳定性。下面介绍时域方面的处理。
对自相关处理单元205得到的自相关结果序列输入到第一峰值检测单元206,查找所需要的检测峰值。
将峰值检测处理后的信号输入到初步确定周期单元207,初步确定信号的周期。
将信号输入到后处理单元208进行处理,这个步骤主要是确认周期,当信号较好的时候,确认的周期可能就是207单元中初步确定的周期;当信号不好时,确认的周期一般都不是207单元中初步确定的周期。以上两种情况都要根据频域提供的参数频域参数更新单元214提供的参数来判断,当时域得到的周期(即207单元得到的周期)和频域得到的频率(即213单元得到的频率)基本对应时,那么周期就基本是正确的,反之要对前面的自相关结果进行二次分析,得到一个或几个较为可信的周期。然后再根据参数更新单元提供的频域提供的参数来选择对应的周期。
在对胎心率信号进行处理时,会在装置中设置信号的预设值,当信号超过该预设的预设值时,即为前面所述的信号较好;当信号未超过该预设的预设值时,即为前面所述的信号不好。由于硬件设备的精确度不同,所得的结果将有所不同。
将经过后处理单元处理后的信号输入到周期计算单元216,,得到所需的确定周期。
将得到的周期输入到心率计算单元217,得到所需的胎儿心率。
将心率结果输入到显示单元218,完成心率的实时显示。
将显示的心率结果输入到打印单元219,完成心率曲线的实时打印。
心率异常判断单元220完成对心率是否异常的判断,并根据不同的情况触发第二报警单元221中对应的报警信号。

Claims (10)

1.一种胎心率信号处理方法,其特征在于,包括以下步骤:
Figure 402469DEST_PATH_IMAGE001
.采集胎心多普勒包络信号;
.对采集的胎心多普勒包络信号进行去噪、滤波处理,得到预处理信号;
Figure 357972DEST_PATH_IMAGE003
.所述预处理信号分为两路信号,一路经过频域分析,另一路经过时域分析;
Figure 838632DEST_PATH_IMAGE004
.对时域分析后确定的胎心信号的初步周期与根据频域分析求得的信号频率进行换算比较,根据比较结果,确认胎心信号的周期;
Figure 209702DEST_PATH_IMAGE005
.对得到的胎心信号的周期进行折算,求取胎儿心率。
2.根据权利要求1所述的一种胎心率信号处理方法,其特征在于,所述频域分析步骤和时域分析步骤不分先后且可同时进行。
3.根据权利要求1所述的一种胎心率信号处理方法,其特征在于,所述频域分析步骤为:
Figure 801220DEST_PATH_IMAGE001
′.判断上述预处理信号是否丢失,若丢失,则返回步骤
Figure 823403DEST_PATH_IMAGE001
.重新采集胎心多普勒包络信号,若未丢失,则
Figure 474964DEST_PATH_IMAGE002
′.对预处理信号做FFT运算;
Figure 127137DEST_PATH_IMAGE003
′.对经FFT运算后的信号进行峰值检测,找出测定频率段内的测定极大值点,定位信号的频率范围;
Figure 522347DEST_PATH_IMAGE004
′.在上述定位的信号频率范围内对预处理信号进行CZT运算;
′.对做完CZT运算后的频率数据进行求极值处理,找出极大值,求得信号频率。
4.根据权利要求1所述的一种胎心率信号处理方法,其特征在于,所述时域分析步骤为:
Figure 972231DEST_PATH_IMAGE006
′.对预处理信号进行自相关处理,得到自相关结果序列;
′.对得到的自相关结果序列进行峰值检测,找到测定峰值;
Figure 500481DEST_PATH_IMAGE008
′.通过上述找到的测定峰值,确定胎心信号的初步周期。
5.根据权利要求1所述的一种胎心率信号处理方法,其特征在于,所述
Figure 982409DEST_PATH_IMAGE007
′.的步骤包括:
.找出自相关结果序列中所有的峰值;
Figure 792419DEST_PATH_IMAGE010
.通过阈值检测得到检测峰值。
6.根据权利要求4所述的一种胎心率信号处理方法,其特征在于,还包括:在步骤
Figure 795010DEST_PATH_IMAGE004
之后,当信号好时,经过处理直接得到信号周期;当信号不好时,则将所述自相关结果进行二次分析,得到修正后的信号周期,再将信号周期折算后与经′步骤得到的信号频率进行比较,根据比较结果,确认胎心信号的周期。
7.一种获取胎心率的装置,其特征在于,包括数据采集单元(201)、预处理单元(202)、信号丢失判断单元(203)、第一报警单元(204)、时域分析处理器、后处理单元(208)、周期计算单元(216)、心率计算单元(217)、心率异常判断单元(220)、第二报警单元(221)、频域分析处理器、参数更新单元(214)、控制单元(215)、显示单元(218)、打印单元(219);
数据采集单元(201)用于采集胎心多普勒包络信号;
预处理单元(202)与数据采集单元(201)连接,用于对数据采集单元采集到的胎心多普勒包络信号进行去噪、滤波处理;
信号丢失判断单元(203)判断经过预处理的信号是否丢失,若信号丢失,则启动第一报警单元(204);
将经预处理后的信号数据分为两路:一路输入频域分析处理器,进行频域分析;另一路输入时域分析处理器,进行时域分析;
将经所述频域分析器分析后的频率信息输入到参数更新单元,保证参数实时性;将更新的频率信息输入到控制单元(215),由控制单元对后处理单元(208)中确认的后处理周期进行判断和控制,保证周期稳定性;
将经所述时域分析处理器时域分析并初步确定信号周期的信号输入到后处理单元(208),根据参数更新单元(214)提供的参数经控制单元(215)的判断和控制确认后处理周期;
将经过后处理单元(208)处理后的信号输入到周期计算单元(216),得到所需的确定周期;
将得到的确定周期输入到心率计算单元(217),得到所需的胎儿心率;
将心率结果输入到显示单元(218),完成心率的实时显示;
将显示的心率结果输入到打印单元(219),完成心率曲线的实时打印;
心率异常判断单元(220)对心率是否异常进行判断,并根据不同情况触发第二报警单元(221)中对应的报警信号。
8.根据权利要求5所述的一种获取胎心率的装置,其特征在于,所述时域分析处理器中包括自相关处理单元(205)、第一峰值检测单元(206)及初步确定周期单元(207);所述经预处理的信号数据输入到自相关处理单元(205),对预处理后的信号数据进行自相关处理,得到自相关结果序列;对所述自相关处理单元(205)得到的自相关结果序列输入到第一峰值检测单元(206),查找所需要的测定峰值;将峰值检测处理后的信号输入到初步确定周期单元(207),确定信号的初步周期。
9.根据权利要求5所述的一种获取胎心率的装置,其特征在于,所述频域分析处理器中包括FFT处理单元(209)、第二峰值检测单元(210)、范围定位单元(211)、CZT处理单元(212)及频率确定单元(213);所述经预处理的信号数据输入到FFT处理单元(209);信号经FFT处理单元(209)处理后,得到频域数据,将该频域数据输入到第二峰值检测单元(210),找出测定频段内的测定极大值点;根据得到的测定极大值点,在范围定位单元(211)中得到信号频率的初步范围;将得到信号频率的初步范围的信号输入CZT处理单元(212),根据信号频率的初步范围对信号做CZT变换;经CZT变换后的信号输入到频率确定单元(213),提取出信号的频率信息。
10.根据权利要求5所述的一种获取胎心率的装置,其特征在于,根据所述参数更新单元(214)提供的参数经控制单元(215)的判断和控制确认后处理周期;当信号好时,则直接将经过后处理单元(208)处理后的信号输入到周期计算单元(216),得到所需的确定周期;当信号不好时,则对所述自相关结果进行二次分析,得到修正后的后处理周期,然后再将经过后处理单元(208)处理后的信号输入到周期计算单元(216),得到所需的确定周期。
CN2011100304520A 2011-01-28 2011-01-28 一种胎心率信号处理方法及其装置 Active CN102048558B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011100304520A CN102048558B (zh) 2011-01-28 2011-01-28 一种胎心率信号处理方法及其装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011100304520A CN102048558B (zh) 2011-01-28 2011-01-28 一种胎心率信号处理方法及其装置

Publications (2)

Publication Number Publication Date
CN102048558A true CN102048558A (zh) 2011-05-11
CN102048558B CN102048558B (zh) 2012-05-23

Family

ID=43953602

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011100304520A Active CN102048558B (zh) 2011-01-28 2011-01-28 一种胎心率信号处理方法及其装置

Country Status (1)

Country Link
CN (1) CN102048558B (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013013646A3 (zh) * 2011-07-26 2013-03-14 深圳市理邦精密仪器股份有限公司 一种提高胎心率数据加速识别准确性的装置和方法
WO2013053290A1 (zh) * 2011-10-13 2013-04-18 深圳市理邦精密仪器股份有限公司 一种提高胎心率数据减速识别准确性的装置和方法
CN103211624A (zh) * 2013-04-07 2013-07-24 深圳市理邦精密仪器股份有限公司 一种提高多普勒胎心率数据准确性的方法与装置
CN103393435A (zh) * 2013-07-31 2013-11-20 深圳市理邦精密仪器股份有限公司 一种胎心音信号包络的取得方法及装置
CN103800037A (zh) * 2014-01-15 2014-05-21 北京春闱科技有限公司 一种胎心监护系统、胎心监护设备以及监护方法
CN104382618A (zh) * 2014-11-13 2015-03-04 深圳市理邦精密仪器股份有限公司 基于胎心率检测的去噪方法和胎心率检测仪
CN104586379A (zh) * 2015-01-21 2015-05-06 深圳市理邦精密仪器股份有限公司 一种胎心率曲线的参数输出方法及装置
CN104873222A (zh) * 2015-05-08 2015-09-02 浙江铭众生物医学创业投资有限公司 一种基于智能移动终端的胎心率检测方法
CN107510462A (zh) * 2017-09-06 2017-12-26 广州视源电子科技股份有限公司 血氧测量方法和装置
CN107865673A (zh) * 2016-09-23 2018-04-03 深圳市理邦精密仪器股份有限公司 胎心率的获取装置及方法、胎心监测设备
CN108577890A (zh) * 2018-03-12 2018-09-28 深圳市理邦精密仪器股份有限公司 胎心率短变异监测方法及装置
WO2019036908A1 (zh) * 2017-08-22 2019-02-28 深圳市汇顶科技股份有限公司 心率检测方法及装置、电子终端
CN110074816A (zh) * 2019-04-26 2019-08-02 深圳市理邦精密仪器股份有限公司 胎心率确定方法、装置及终端设备
WO2020061887A1 (zh) * 2018-09-27 2020-04-02 深圳市伊欧乐科技有限公司 心率测量方法、装置及计算机可读存储介质
CN111839491A (zh) * 2020-06-29 2020-10-30 北京曙光易通技术有限公司 心搏功能监测方法、心搏功能持续监测方法及电子设备
CN112869724A (zh) * 2021-01-19 2021-06-01 西安交通大学 一种基于多通道被动式采集信号的胎儿健康监测仪

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10285651B2 (en) * 2016-06-30 2019-05-14 Analog Devices, Inc. On-demand heart rate estimation based on optical measurements

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4781200A (en) * 1985-10-04 1988-11-01 Baker Donald A Ambulatory non-invasive automatic fetal monitoring system
CN1195277A (zh) * 1996-04-17 1998-10-07 精工爱普生株式会社 心律失常检测设备
US20010014776A1 (en) * 1994-09-21 2001-08-16 Nancy E. Oriol Fetal data processing system and method employing a time-frequency representation
DE10345717A1 (de) * 2003-10-01 2005-04-28 Trium Analysis Online Gmbh Verfahren und Vorrichtung zur Bestimmung der fötalen Herzfrequenz
CN101689219A (zh) * 2006-11-01 2010-03-31 必安康医疗有限公司 用于监测心肺参数的系统和方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4781200A (en) * 1985-10-04 1988-11-01 Baker Donald A Ambulatory non-invasive automatic fetal monitoring system
US20010014776A1 (en) * 1994-09-21 2001-08-16 Nancy E. Oriol Fetal data processing system and method employing a time-frequency representation
CN1195277A (zh) * 1996-04-17 1998-10-07 精工爱普生株式会社 心律失常检测设备
DE10345717A1 (de) * 2003-10-01 2005-04-28 Trium Analysis Online Gmbh Verfahren und Vorrichtung zur Bestimmung der fötalen Herzfrequenz
CN101689219A (zh) * 2006-11-01 2010-03-31 必安康医疗有限公司 用于监测心肺参数的系统和方法

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9801609B2 (en) 2011-07-26 2017-10-31 Edan Instruments Inc. Device and method for enhancing accuracy of recognizing fetus heart rate acceleration data
WO2013013646A3 (zh) * 2011-07-26 2013-03-14 深圳市理邦精密仪器股份有限公司 一种提高胎心率数据加速识别准确性的装置和方法
US9364197B2 (en) 2011-10-13 2016-06-14 Edan Instruments, Inc Device and method for improving identification accuracy of fetal heart rate deceleration
WO2013053290A1 (zh) * 2011-10-13 2013-04-18 深圳市理邦精密仪器股份有限公司 一种提高胎心率数据减速识别准确性的装置和方法
CN103211624A (zh) * 2013-04-07 2013-07-24 深圳市理邦精密仪器股份有限公司 一种提高多普勒胎心率数据准确性的方法与装置
CN103393435A (zh) * 2013-07-31 2013-11-20 深圳市理邦精密仪器股份有限公司 一种胎心音信号包络的取得方法及装置
CN103800037A (zh) * 2014-01-15 2014-05-21 北京春闱科技有限公司 一种胎心监护系统、胎心监护设备以及监护方法
CN104382618B (zh) * 2014-11-13 2017-09-26 深圳市理邦精密仪器股份有限公司 基于胎心率检测的去噪方法和胎心率检测仪
CN104382618A (zh) * 2014-11-13 2015-03-04 深圳市理邦精密仪器股份有限公司 基于胎心率检测的去噪方法和胎心率检测仪
CN104586379B (zh) * 2015-01-21 2017-05-24 深圳市理邦精密仪器股份有限公司 一种胎心率曲线的参数输出方法及装置
CN104586379A (zh) * 2015-01-21 2015-05-06 深圳市理邦精密仪器股份有限公司 一种胎心率曲线的参数输出方法及装置
CN104873222A (zh) * 2015-05-08 2015-09-02 浙江铭众生物医学创业投资有限公司 一种基于智能移动终端的胎心率检测方法
CN107865673A (zh) * 2016-09-23 2018-04-03 深圳市理邦精密仪器股份有限公司 胎心率的获取装置及方法、胎心监测设备
WO2019036908A1 (zh) * 2017-08-22 2019-02-28 深圳市汇顶科技股份有限公司 心率检测方法及装置、电子终端
US11096628B2 (en) 2017-08-22 2021-08-24 Shenzhen GOODIX Technology Co., Ltd. Heart rate detection method and apparatus, and electronic terminal thereof
CN107510462B (zh) * 2017-09-06 2020-05-05 广州视源电子科技股份有限公司 血氧测量方法和装置
CN107510462A (zh) * 2017-09-06 2017-12-26 广州视源电子科技股份有限公司 血氧测量方法和装置
CN108577890A (zh) * 2018-03-12 2018-09-28 深圳市理邦精密仪器股份有限公司 胎心率短变异监测方法及装置
WO2020061887A1 (zh) * 2018-09-27 2020-04-02 深圳市伊欧乐科技有限公司 心率测量方法、装置及计算机可读存储介质
CN110074816A (zh) * 2019-04-26 2019-08-02 深圳市理邦精密仪器股份有限公司 胎心率确定方法、装置及终端设备
CN111839491A (zh) * 2020-06-29 2020-10-30 北京曙光易通技术有限公司 心搏功能监测方法、心搏功能持续监测方法及电子设备
CN112869724A (zh) * 2021-01-19 2021-06-01 西安交通大学 一种基于多通道被动式采集信号的胎儿健康监测仪
CN112869724B (zh) * 2021-01-19 2022-04-22 西安交通大学 一种基于多通道被动式采集信号的胎儿健康监测仪

Also Published As

Publication number Publication date
CN102048558B (zh) 2012-05-23

Similar Documents

Publication Publication Date Title
CN102048558B (zh) 一种胎心率信号处理方法及其装置
RU2531119C2 (ru) Бесконтактный мониторинг дыхания пациента
US20090062654A1 (en) Method and apparatus for automatic optimization of doppler imaging parameters
CN101766496B (zh) 一种噪声估计方法和一键优化方法及其系统
CN104382618B (zh) 基于胎心率检测的去噪方法和胎心率检测仪
CN106443178A (zh) 一种基于IQuinn‑Rife综合的正弦信号频率估计方法
WO2020072297A1 (en) Direct rf signal processing for heart-rate monitoring using uwb impulse radar
CN104820786A (zh) 一种瞬时加权同步挤压小波双谱分析方法
CN108652640B (zh) 一种基于心电信号的无创血糖检测方法及系统
CN115644840A (zh) 基于毫米波雷达的生命体征检测方法
US8162838B2 (en) Doppler signal spectrum calculating method and apparatus
KR101779018B1 (ko) 초음파 도플러 태아감시 장치의 심박 검출 신호처리 방법
CN106175731A (zh) 非接触式生命体征监测的信号处理系统
CN103494605B (zh) 心率检测方法及装置
CN102809687B (zh) 一种交流电频率的数字化测量方法
CN112965060A (zh) 生命特征参数的检测方法、装置和检测体征点的方法
JP2017127398A (ja) 情報処理装置、情報処理システム、情報処理方法及びプログラム
CN103190913A (zh) 一种胎动识别方法、设备及胎儿监护仪
CN109584902B (zh) 一种音乐节奏确定方法、装置、设备及存储介质
CN104352234A (zh) 一种生理电信号尖峰奇异点检测方法
CN116701840A (zh) 一种机械振动信号的倒频谱优化计算方法和系统
WO2018160890A1 (en) Efficient fetal-maternal ecg signal separation from two maternal abdominal leads via diffusion-based channel selection
CN114052709B (zh) 一种稳健的毫米波雷达生命体征测量方法
CN104873222B (zh) 一种基于智能移动终端的胎心率检测方法
CN113786180A (zh) 基于毫米波雷达信号的人体血压检测方法

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