CN105193451B - 一种胎儿心率曲线基线的识别方法 - Google Patents

一种胎儿心率曲线基线的识别方法 Download PDF

Info

Publication number
CN105193451B
CN105193451B CN201510694998.4A CN201510694998A CN105193451B CN 105193451 B CN105193451 B CN 105193451B CN 201510694998 A CN201510694998 A CN 201510694998A CN 105193451 B CN105193451 B CN 105193451B
Authority
CN
China
Prior art keywords
mrow
mtd
baseline
fetal heart
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.)
Active
Application number
CN201510694998.4A
Other languages
English (en)
Other versions
CN105193451A (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.)
Beijing University of Technology
Peking University Third Hospital
Original Assignee
Beijing University of Technology
Peking University Third Hospital
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 Beijing University of Technology, Peking University Third Hospital filed Critical Beijing University of Technology
Priority to CN201510694998.4A priority Critical patent/CN105193451B/zh
Publication of CN105193451A publication Critical patent/CN105193451A/zh
Application granted granted Critical
Publication of CN105193451B publication Critical patent/CN105193451B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

一种胎儿心率信号曲线的基线识别方法,包括信号输入检测模块,基线识别提取模块,基线值参数计算模块,基线值参数输出模块。胎儿心率信号输入后通过预处理及有效信号的提取,在基线识别提取模块中进行预计算、平均优化和插值平滑后计算出基线值,之后在参数计算模块计算出基线区间比例和描述性统计参数,最后在参数输出模块输出结果。基线值是利用胎儿心率监护长时间的胎儿心率数据以及胎儿心率基线定义计算获得,具有较明确的生理意义,结合胎儿心率的加速、减速、变异的分类、计数和所占时长比例,能够从胎儿心率信号中获得更为全面的信息。应用提取出的基线值,还可以获取胎儿心率的频域以及熵的相关信息。

Description

一种胎儿心率曲线基线的识别方法
技术领域
本发明涉及生物医学信号处理领域,具体是一种胎儿心率监护中胎儿心率曲线的基线计算识别方法。
背景技术
电子胎儿心率监护作为一种实时获取胎儿心率信号及母体宫缩信号的宫内监护技术,目前在产科临床广泛应用。胎儿心率曲线中蕴含着丰富的有关胎儿的神经与心血管功能信息。由于胎儿位于母体腹部子宫内,胎儿心音信号很容易受到来自母体、传感器以及外界环境的噪声的干扰。目前胎儿心率监护系统大都应用超声多普勒原理并采用自相关算法计算胎儿瞬时心率。胎儿心率曲线所表现出来的形态、振幅强度与一定的周期性等方面的综合信息在相当程度上反映出胎儿在宫内的生理状态。胎儿心率信号的基线作为胎儿电子监护图形中的重要组成部分,是除去加速、减速与显著变异节段的在特定胎儿心率波动范围内的平均胎儿心率,在监护图形上胎儿心率信号围绕基线上下波动。
目前对胎儿心率信号基线的评估主要是通过手工测量与计算和肉眼识别的方法来判定胎儿心率监护指标,对胎儿心率信号的分析与研究不够准确,对基线所含有的信息提取不够充分。利用数字信号处理技术,对胎儿心率信号加以数字化,可获取更为客观的参数值,有助于得到更为全面的胎儿宫内生理信息。
发明内容
现有的技术方案主要是利用多路电极探头采集孕妇与胎儿的心电信号,对心电信号进行滤波等信号处理得到胎儿心率,应用胎儿心率与电子胎儿心率监护中基线的定义对胎儿心率基线进行计算。产前胎儿的监护是一个综合复杂的过程,以手工测量的方法从监护图形中获取的胎儿心率基线信息不够客观,且差异性较大,仅依靠肉眼对胎儿心率基线进行识别,效率较低并且会丢失除时域以外在频域或者由熵提供的更多信息。
为解决上述问题,本发明提出了一种胎儿心率信号曲线的基线识别方法,基线值是利用胎儿心率监护长时间的胎儿心率数据以及胎儿心率基线定义计算获得,具有较明确的生理意义,结合胎儿心率的加速、减速、变异的分类、计数和所占时长比例,能够从胎儿心率信号中获得更为全面的信息。获得的基线通过与特定阈值对比也可以确定一定的生理或病理特征,由此可以改善现有技术对胎儿心率信号中的信息获取较为片面的情况。应用提取出的基线值,还可以获取胎儿心率的频域以及熵的相关信息。
为达到上述目的,本发明所采取的技术方案为:一种胎儿心率信号曲线的基线识别方法,其特征在于:包括信号输入检测模块,基线识别提取模块,基线值参数计算模块,基线值参数输出模块。胎儿心率信号输入后通过预处理及有效信号的提取,在基线识别提取模块中进行预计算、平均优化和插值平滑后计算出基线值,之后在参数计算模块计算出基线区间比例和描述性统计参数,最后在参数输出模块输出结果。
所述信号输入模块,接收胎儿心率监护仪器实测的信号F(t),F(t)中包含母亲心率信号、胎儿心率信号、母亲移动次数和子宫肌电信号,并从信号F(t)中分离出胎儿心率信号f(t);
所述信号检测模块,对输入的时域胎儿心率信号f(t)进行预处理,即将f(t)分成N段时间长度的胎儿心率信号序列f(k),每一段胎儿心率信号序列f(k)中胎儿心率值为0的为丢失信号f(k0),其长度为l(f(k0)),与该段胎儿心率信号序列长度l(f(k))之比为信号丢失率挑选出满足信号丢失率lost小于特定阈值low的胎儿心率信号序列g(k),即其中g(k0)为胎儿心率信号序列g(k)中胎儿心率值为0的丢失信号;有效胎儿心率信号时间总长度为总的信号丢失率为
所述基线识别提取模块接收满足要求的胎儿心率信号并根据基线的定义进行遍历预计算、平均优化和插值平滑,以下为基线具体识别提取方法;
首先进行遍历预计算,对于输入的满足要求的胎儿心率信号bpm(k),将其分为定长为一分钟的信号M段,长度不足一分钟的舍弃,第h段为d(h),并折叠为如下矩阵A:
矩阵A为(M-9)行10列的矩阵,则矩阵A的每行为定长十分钟的胎儿心率信号,对矩阵A的第iA行设为A(iA),其中胎儿心率信号为A(iA(k));特定胎儿心率值B∈[Bmin,Bmax],步长为1,其中Bmax和Bmin分别为胎儿心率最大阈值与最小阈值;对于每段定长十分钟的胎儿心率信号A(iA(k)),计算出满足B-5≤A(iA(k))≤B+5的个数numiA,对区间[Bmin,Bmax]内每一个值B进行个数计算,得出中的最大值此时的特定B值命名为则矩阵A的每一行都有一个和相对应的共计(M-9)个;
对于iA=2,3,4…(M-9),如果满足不变;如果满足的值等于前一个的值;
然后做平均优化处理,计算满足的个数并将满足条件的胎儿心率信号求和平均优化处理后的值为
最后进行插值平滑处理,对相邻两个值做差并求绝对值,对于iA=2,3,4…(M-9),如果满足dif(iA-1)>5,则在iA≤4时,将置为空值;在iA>4时,将置为空值;接下来在空值区域处进行插值,设空值区域前的值为bx,空值区域后的值为by(bx与by都不为0),则此空值区域缺少y-x-1个值;
空值区域的插值项如下:
至此存在值(M-9)个,胎儿心率信号bpm(k)包括时长一分钟的信号M段,设为blj为基线值,其中(j=1,2,3…M)共有M个值;
基线值blj表达式如下:
此时,每分钟有一个基线值,若将基线曲线与胎儿心率信号绘制在一起,需要再次将基线值blj进行插值;对基线值进行三次样条插值得到基线曲线line=spline(blj),其中spline为三次样条插值函数;
所述基线参数计算模块,包括特定区间内基线值所占比例计算与统计分析参数计算;设定正常胎儿心基线范围上限为BLmax,下限为BLmin;将基线值blj与正常胎儿心基线范围上限与下限比较,满足blj>BLmax条件的个数为bnum1,所占比例为Rh;满足BLmin≤blj≤BLmax条件的个数为bnum2,所占比例为Rn;满足blj<BLmin条件的个数为bnum3,所占比例为Rl;特定区间内基线值所占比例计算公式如下:
经统计分析得出基线中位数为median;基线均值公式如下:
所述基线参数输出模块,将上述计算出的胎儿心率监护基线参数结果输出。
胎儿心率曲线的基线识别提取的有益效果在于:
该方法仅需输入胎儿心率监护仪器的实测原始信号,即可分离出胎儿心率信号,并通过遍历、平滑等步骤获取该胎儿心率信号的基线值与相关参数,并绘制基线曲线。步骤简单方便,提出了一种新型的基线曲线提取方法,改善了关于胎儿心率基线的特定参数较难获得的情况。在实际应用过程中可减少通过目测方法带来的误差,有效提高关于胎儿心率监护图的解读效率。
附图说明
图1为胎儿心监护基线特征参数提取结构框图。
图2为胎儿心监护基线特征参数提取操作流程图。
图3为胎儿心率曲线及对应基线示意图。
具体实施方式
以下结合附图对本发明的一种较为典型的具体实施方式进行详细描述。
本发明的一种典型应用场景在于利用胎儿心率监护实测信号获取基线及其特征参数。由此可以利用现有胎儿心率信号采集的成熟技术及基线曲线提取方法获取更为客观全面的关于胎儿心率基线特征的信息。
如图2所示,以一位年龄为32周岁,孕周为37周的被测者为例。
步骤T1中系统开始接收胎儿心率监护仪器实测的信号F(t),并从F(t)中提取处胎儿心率信号f(t)。
步骤T2对输入的胎儿心率信号f(t)进行分段,计算出每段信号序列f(k)的丢失率lost以及整体胎心率信号的丢失率Lost,被测者整体信号丢失率Lost为3.92%。
步骤T3中根据丢失率进行信号有效性判断,对丢失率lost高于特定阈值low的原始信号序列予以舍弃。若出现信号序列f(k)舍弃,重新对下一段信号序列f(k+1)进行计算与判断。
步骤T4对步骤T3中满足要求的有效信号bpm(k)进行基线识别提取,提取出基线值blj
步骤T5中按需求设定正常胎儿心基线范围上限BLmax与正常胎儿心基线范围下限BLmin,这里赋值为BLmax=160,BLmin=110。通过对整体胎儿心率信号的处理,计算出低于正常胎儿心基线范围占比Rl、正常胎儿心基线值占比Rn、高于正常胎儿心基线范围占比Rh,表达式为:
经统计分析得出基线中位数为median。
基线均值公式如下:
经计算可得:
Rh=2.99% Rn=95.28% Rl=1.73% median=129.8 mean=131.9
步骤T6将上述参数按照指定格式输出。

Claims (1)

1.一种胎儿心率信号曲线的基线识别方法,其特征在于:包括信号输入检测模块,基线识别提取模块,基线值参数计算模块,基线值参数输出模块;胎儿心率信号输入后通过预处理及有效信号的提取,在基线识别提取模块中进行预计算、平均优化和插值平滑后计算出基线值,之后在参数计算模块计算出基线区间比例和描述性统计参数,最后在参数输出模块输出结果;
所述信号输入模块,接收胎儿心率监护仪器实测的信号F(t),F(t)中包含母亲心率信号、胎儿心率信号、母亲移动次数和子宫肌电信号,并从信号F(t)中分离出胎儿心率信号f(t);
所述信号检测模块,对输入的时域胎儿心率信号f(t)进行预处理,即将f(t)分成N段时间长度的胎儿心率信号序列f(k),每一段胎儿心率信号序列f(k)中胎儿心率值为0的为丢失信号f(k0),其长度为l(f(k0)),与该段胎儿心率信号序列长度l(f(k))之比为信号丢失率挑选出满足信号丢失率lost小于特定阈值low的胎儿心率信号序列g(k),即其中g(k0)为胎儿心率信号序列g(k)中胎儿心率值为0的丢失信号;有效胎儿心率信号时间总长度为总的信号丢失率为
所述基线识别提取模块接收满足要求的胎儿心率信号并根据基线的定义进行遍历预计算、平均优化和插值平滑,以下为基线具体识别提取方法;
首先进行遍历预计算,对于输入的满足要求的胎儿心率信号bpm(k),将其分为定长为一分钟的信号M段,长度不足一分钟的舍弃,第h段为d(h),并折叠为如下矩阵A:
<mrow> <mi>A</mi> <mo>=</mo> <msub> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>-</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>-</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>-</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>-</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mrow> <mo>(</mo> <mi>M</mi> <mo>-</mo> <mn>9</mn> <mo>)</mo> <mo>&amp;times;</mo> <mn>10</mn> </mrow> </msub> </mrow>
矩阵A为(M-9)行10列的矩阵,则矩阵A的每行为定长十分钟的胎儿心率信号,对矩阵A的第iA行设为A(iA),其中胎儿心率信号为A(iA(k));特定胎儿心率值B∈[Bmin,Bmax],步长为1,其中Bmax和Bmin分别为胎儿心率最大阈值与最小阈值;对于每段定长十分钟的胎儿心率信号A(iA(k)),计算出满足B-5≤A(iA(k))≤B+5的个数对区间[Bmin,Bmax]内每一个值B进行个数计算,得出中的最大值此时的特定B值命名为则矩阵A的每一行都有一个和相对应的共计(M-9)个;
对于iA=2,3,4…(M-9),如果满足不变;如果满足的值等于前一个的值;
然后做平均优化处理,计算满足的个数并将满足条件的胎儿心率信号求和平均优化处理后的值为
最后进行插值平滑处理,对相邻两个值做差并求绝对值,对于iA=2,3,4…(M-9),如果满足dif(iA-1)>5,则在iA≤4时,将置为空值;在iA>4时,将置为空值;接下来在空值区域处进行插值,设空值区域前的值为bx,空值区域后的值为by,bx与by都不为0,则此空值区域缺少y-x-1个值;
空值区域的插值项如下:
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>b</mi> <mrow> <mi>x</mi> <mo>+</mo> <mi>a</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>b</mi> <mi>x</mi> </msub> <mo>+</mo> <mfrac> <mrow> <msub> <mi>b</mi> <mi>y</mi> </msub> <mo>-</mo> <msub> <mi>b</mi> <mi>x</mi> </msub> </mrow> <mrow> <mi>y</mi> <mo>-</mo> <mi>x</mi> </mrow> </mfrac> <mo>*</mo> <mi>a</mi> </mrow> </mtd> <mtd> <mrow> <mo>(</mo> <mi>a</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2...</mn> <mo>(</mo> <mrow> <mi>y</mi> <mo>-</mo> <mi>x</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced>
至此存在值(M-9)个,胎儿心率信号bpm(k)包括时长一分钟的信号M段,设为blj为基线值,其中j=1,2,3…M共有M个值;
基线值blj表达式如下:
此时,每分钟有一个基线值,若将基线曲线与胎儿心率信号绘制在一起,需要再次将基线值blj进行插值;对基线值进行三次样条插值得到基线曲线line=spline(blj),其中spline为三次样条插值函数;
所述基线参数计算模块,包括特定区间内基线值所占比例计算与统计分析参数计算;设定正常胎儿心基线范围上限为BLmax,下限为BLmin;将基线值blj与正常胎儿心基线范围上限与下限比较,满足blj>BLmax条件的个数为bnum1,所占比例为Rh;满足BLmin≤blj≤BLmax条件的个数为bnum2,所占比例为Rn;满足blj<BLmin条件的个数为bnum3,所占比例为Rl;特定区间内基线值所占比例计算公式如下:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <msub> <mi>R</mi> <mi>h</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>b</mi> <mi>n</mi> <mi>u</mi> <mi>m</mi> <mn>1</mn> </mrow> <mi>M</mi> </mfrac> </mtd> </mtr> <mtr> <mtd> <msub> <mi>R</mi> <mi>n</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>b</mi> <mi>n</mi> <mi>u</mi> <mi>m</mi> <mn>2</mn> </mrow> <mi>M</mi> </mfrac> </mtd> </mtr> <mtr> <mtd> <msub> <mi>R</mi> <mi>l</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>b</mi> <mi>n</mi> <mi>u</mi> <mi>m</mi> <mn>3</mn> </mrow> <mi>M</mi> </mfrac> </mtd> </mtr> </mtable> </mfenced>
经统计分析得出基线中位数为median;基线均值公式如下:
<mrow> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>n</mi> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msub> <mi>bl</mi> <mi>j</mi> </msub> </mrow> <mi>M</mi> </mfrac> </mrow>
所述基线参数输出模块,将上述计算出的胎儿心率监护基线参数结果输出。
CN201510694998.4A 2015-10-21 2015-10-21 一种胎儿心率曲线基线的识别方法 Active CN105193451B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510694998.4A CN105193451B (zh) 2015-10-21 2015-10-21 一种胎儿心率曲线基线的识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510694998.4A CN105193451B (zh) 2015-10-21 2015-10-21 一种胎儿心率曲线基线的识别方法

Publications (2)

Publication Number Publication Date
CN105193451A CN105193451A (zh) 2015-12-30
CN105193451B true CN105193451B (zh) 2018-02-27

Family

ID=54941685

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510694998.4A Active CN105193451B (zh) 2015-10-21 2015-10-21 一种胎儿心率曲线基线的识别方法

Country Status (1)

Country Link
CN (1) CN105193451B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106073746B (zh) * 2016-06-15 2018-10-19 北京工业大学 基于胎儿心率曲线的胎儿安静睡眠周期提取方法
CN110866492B (zh) * 2019-11-13 2022-12-13 广州品唯软件有限公司 一种基线分支的识别方法、装置及计算机系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102302363A (zh) * 2011-06-29 2012-01-04 深圳市理邦精密仪器股份有限公司 一种胎心率基线的自动识别装置和方法
CN102319064A (zh) * 2011-10-13 2012-01-18 深圳市理邦精密仪器股份有限公司 一种提高胎心率数据减速识别准确性的装置和方法
CN102319063A (zh) * 2011-07-26 2012-01-18 深圳市理邦精密仪器股份有限公司 一种提高胎心率数据加速识别准确性的装置和方法
CN102940485A (zh) * 2012-10-10 2013-02-27 广州三瑞医疗器械有限公司 一种胎心率基线自动识别算法
CN104287711A (zh) * 2014-09-24 2015-01-21 广州三瑞医疗器械有限公司 胎心率曲线非基线部分和基线的计算方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2447861A1 (en) * 2002-11-01 2004-05-01 Lms Medical Systems Ltd. Method and apparatus for identifying heart rate feature events
US20190298210A1 (en) * 2013-07-01 2019-10-03 Mayo Foundation For Medical Education And Research Algorithms for managing artifact and detecting cardiac events using a patient monitoring system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102302363A (zh) * 2011-06-29 2012-01-04 深圳市理邦精密仪器股份有限公司 一种胎心率基线的自动识别装置和方法
CN102319063A (zh) * 2011-07-26 2012-01-18 深圳市理邦精密仪器股份有限公司 一种提高胎心率数据加速识别准确性的装置和方法
CN102319064A (zh) * 2011-10-13 2012-01-18 深圳市理邦精密仪器股份有限公司 一种提高胎心率数据减速识别准确性的装置和方法
CN102940485A (zh) * 2012-10-10 2013-02-27 广州三瑞医疗器械有限公司 一种胎心率基线自动识别算法
CN104287711A (zh) * 2014-09-24 2015-01-21 广州三瑞医疗器械有限公司 胎心率曲线非基线部分和基线的计算方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A novel cardiotocography fetal heart rate baseline estimation algorithm;Nidhal S, Ali M A M, Najah H;《Scientific Research and Essays》;20101231;第5卷(第24期);第4002-4010页 *
An estimate of fetal autonomic state by time-frequency analysis of fetal heart rate variability;David M, Hirsch M, Karin J等;《Journal of Applied Physiology》;20071231;第102卷(第3期);第1057-1064页 *
Estimation and classification of fetal heart rate baselines using artificial neural networks;de Sa J P M, Reis L P, Lau J N等;《Computers in Cardiology》;19941231;第541-544页 *
利用三次样条差值法抑制脉搏波基线漂移;杨琳,张松等;《北京生物医学工程》;20100430;第29卷(第2期);第198-200页 *
基于单片机和自相关计算的胎心检测仪的实现;朱韶红等;《中国科技信息》;20061231(第23期);第186-187页 *

Also Published As

Publication number Publication date
CN105193451A (zh) 2015-12-30

Similar Documents

Publication Publication Date Title
CN104665849B (zh) 一种基于多生理信号多模型交互的高铁调度员应激检测方法
CN105877766B (zh) 一种基于多生理信号融合的精神状态检测系统及方法
CN102319063B (zh) 一种提高胎心率数据加速识别准确性的装置和方法
CN102319064B (zh) 一种提高胎心率数据减速识别准确性的装置和方法
CN106073746B (zh) 基于胎儿心率曲线的胎儿安静睡眠周期提取方法
CN106934125A (zh) 一种梯形噪声分布的指数模型机械设备剩余寿命预测方法
CN104887220B (zh) 一种由腹壁心电信号提取胎儿心电信号的方法及系统
CN104546007B (zh) 胎动检测抗干扰处理方法及装置
CN107095669B (zh) 一种癫痫患者脑电信号的处理方法及系统
CN102940485B (zh) 一种胎心率基线自动识别算法
CN104504377B (zh) 一种公交车乘客拥挤程度识别系统及方法
CN103083013A (zh) 基于形态学与小波变换的心电信号qrs复波检测方法
Wróbel et al. Improving fetal heart rate signal interpretation by application of myriad filtering
CN105193451B (zh) 一种胎儿心率曲线基线的识别方法
CN104490373B (zh) 脉搏信号的判断方法、判断装置以及生理参数测量设备
CN104055506B (zh) 一种胎监数据处理方法及装置
CN106889987A (zh) 基于分段滤波处理的子宫肌电强度信息提取方法
Gao et al. Fetal heart baseline extraction and classification based on deep learning
CN116725553A (zh) 基于时频域关联特征的睡眠循环交替模式检测方法及装置
CN103942865B (zh) 基于超声传感器的地铁防逃票系统及其设计方法
CN106781397A (zh) 一种牛羊监管系统
CN117354468B (zh) 一种基于大数据的智能化状态感知系统及方法
CN105030231B (zh) 一种胎心监护加速特征参数的提取方法
CN202758414U (zh) 一种公交车智能投币机
CN107977962A (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
TA01 Transfer of patent application right

Effective date of registration: 20170731

Address after: 100124 Chaoyang District, Beijing Ping Park, No. 100

Applicant after: Beijing University of Technology

Applicant after: The Third Affiliated Hospital of Peking University

Address before: 100124 Chaoyang District, Beijing Ping Park, No. 100

Applicant before: Beijing University of Technology

TA01 Transfer of patent application right
CB03 Change of inventor or designer information

Inventor after: Zhang Song

Inventor after: Li Xuwen

Inventor after: Chen Lian

Inventor after: Lu Lu

Inventor after: Zhang Yanhong

Inventor after: Li Guangfei

Inventor after: Zhao Yangyu

Inventor after: Wang Yan

Inventor after: Li Shufang

Inventor after: Zhang Lei

Inventor after: Yang Lin

Inventor after: Hao Dongmei

Inventor after: Yang Yimin

Inventor before: Zhang Song

Inventor before: Li Guangfei

Inventor before: Yang Lin

Inventor before: Yang Yimin

Inventor before: Li Xuwen

Inventor before: Wang Yan

Inventor before: Li Shufang

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant