CN108761105B - 凝血时间判定方法、电子设备及存储介质 - Google Patents
凝血时间判定方法、电子设备及存储介质 Download PDFInfo
- Publication number
- CN108761105B CN108761105B CN201810554622.7A CN201810554622A CN108761105B CN 108761105 B CN108761105 B CN 108761105B CN 201810554622 A CN201810554622 A CN 201810554622A CN 108761105 B CN108761105 B CN 108761105B
- Authority
- CN
- China
- Prior art keywords
- time
- coagulation
- slope
- blood coagulation
- point search
- 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
- 230000023555 blood coagulation Effects 0.000 title claims abstract description 196
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000005345 coagulation Methods 0.000 claims abstract description 247
- 230000015271 coagulation Effects 0.000 claims abstract description 239
- 238000005070 sampling Methods 0.000 claims abstract description 154
- 238000012360 testing method Methods 0.000 claims abstract description 59
- 230000015654 memory Effects 0.000 claims description 16
- 108090000190 Thrombin Proteins 0.000 claims description 11
- 229960004072 thrombin Drugs 0.000 claims description 11
- 108010094028 Prothrombin Proteins 0.000 claims description 9
- 102100027378 Prothrombin Human genes 0.000 claims description 9
- 229940039716 prothrombin Drugs 0.000 claims description 9
- PGOHTUIFYSHAQG-LJSDBVFPSA-N (2S)-6-amino-2-[[(2S)-5-amino-2-[[(2S)-2-[[(2S)-2-[[(2S)-2-[[(2S)-4-amino-2-[[(2S)-2-[[(2S)-2-[[(2S)-2-[[(2S)-2-[[(2S)-5-amino-2-[[(2S)-5-amino-2-[[(2S)-2-[[(2S)-2-[[(2S)-2-[[(2S,3R)-2-[[(2S)-5-amino-2-[[(2S)-2-[[(2S)-2-[[(2S,3R)-2-[[(2S)-2-[[(2S)-2-[[(2S)-2-[[(2S)-2-[[(2S)-5-amino-2-[[(2S)-1-[(2S,3R)-2-[[(2S)-2-[[(2S)-2-[[(2R)-2-[[(2S)-2-[[(2S)-2-[[2-[[(2S)-2-[[(2S)-2-[[(2S)-2-[[(2S)-1-[(2S)-2-[[(2S)-2-[[(2S)-2-[[(2S)-2-amino-4-methylsulfanylbutanoyl]amino]-3-(1H-indol-3-yl)propanoyl]amino]-5-carbamimidamidopentanoyl]amino]propanoyl]pyrrolidine-2-carbonyl]amino]-3-methylbutanoyl]amino]-4-methylpentanoyl]amino]-4-methylpentanoyl]amino]acetyl]amino]-3-hydroxypropanoyl]amino]-4-methylpentanoyl]amino]-3-sulfanylpropanoyl]amino]-4-methylsulfanylbutanoyl]amino]-5-carbamimidamidopentanoyl]amino]-3-hydroxybutanoyl]pyrrolidine-2-carbonyl]amino]-5-oxopentanoyl]amino]-3-hydroxypropanoyl]amino]-3-hydroxypropanoyl]amino]-3-(1H-imidazol-5-yl)propanoyl]amino]-4-methylpentanoyl]amino]-3-hydroxybutanoyl]amino]-3-(1H-indol-3-yl)propanoyl]amino]-5-carbamimidamidopentanoyl]amino]-5-oxopentanoyl]amino]-3-hydroxybutanoyl]amino]-3-hydroxypropanoyl]amino]-3-carboxypropanoyl]amino]-3-hydroxypropanoyl]amino]-5-oxopentanoyl]amino]-5-oxopentanoyl]amino]-3-phenylpropanoyl]amino]-5-carbamimidamidopentanoyl]amino]-3-methylbutanoyl]amino]-4-methylpentanoyl]amino]-4-oxobutanoyl]amino]-5-carbamimidamidopentanoyl]amino]-3-(1H-indol-3-yl)propanoyl]amino]-4-carboxybutanoyl]amino]-5-oxopentanoyl]amino]hexanoic acid Chemical compound CSCC[C@H](N)C(=O)N[C@@H](Cc1c[nH]c2ccccc12)C(=O)N[C@@H](CCCNC(N)=N)C(=O)N[C@@H](C)C(=O)N1CCC[C@H]1C(=O)N[C@@H](C(C)C)C(=O)N[C@@H](CC(C)C)C(=O)N[C@@H](CC(C)C)C(=O)NCC(=O)N[C@@H](CO)C(=O)N[C@@H](CC(C)C)C(=O)N[C@@H](CS)C(=O)N[C@@H](CCSC)C(=O)N[C@@H](CCCNC(N)=N)C(=O)N[C@@H]([C@@H](C)O)C(=O)N1CCC[C@H]1C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H](CO)C(=O)N[C@@H](CO)C(=O)N[C@@H](Cc1cnc[nH]1)C(=O)N[C@@H](CC(C)C)C(=O)N[C@@H]([C@@H](C)O)C(=O)N[C@@H](Cc1c[nH]c2ccccc12)C(=O)N[C@@H](CCCNC(N)=N)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H]([C@@H](C)O)C(=O)N[C@@H](CO)C(=O)N[C@@H](CC(O)=O)C(=O)N[C@@H](CO)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H](Cc1ccccc1)C(=O)N[C@@H](CCCNC(N)=N)C(=O)N[C@@H](C(C)C)C(=O)N[C@@H](CC(C)C)C(=O)N[C@@H](CC(N)=O)C(=O)N[C@@H](CCCNC(N)=N)C(=O)N[C@@H](Cc1c[nH]c2ccccc12)C(=O)N[C@@H](CCC(O)=O)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H](CCCCN)C(O)=O PGOHTUIFYSHAQG-LJSDBVFPSA-N 0.000 claims description 8
- 108010000499 Thromboplastin Proteins 0.000 claims description 8
- 102000002262 Thromboplastin Human genes 0.000 claims description 8
- 206010053567 Coagulopathies Diseases 0.000 abstract description 21
- 230000035602 clotting Effects 0.000 abstract description 21
- 238000007405 data analysis Methods 0.000 abstract description 7
- 230000008859 change Effects 0.000 description 22
- 239000008280 blood Substances 0.000 description 18
- 210000004369 blood Anatomy 0.000 description 17
- 239000003153 chemical reaction reagent Substances 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 5
- 238000004891 communication Methods 0.000 description 5
- 210000003743 erythrocyte Anatomy 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 108010039209 Blood Coagulation Factors Proteins 0.000 description 3
- 102000015081 Blood Coagulation Factors Human genes 0.000 description 3
- 102000009123 Fibrin Human genes 0.000 description 3
- 108010073385 Fibrin Proteins 0.000 description 3
- BWGVNKXGVNDBDI-UHFFFAOYSA-N Fibrin monomer Chemical compound CNC(=O)CNC(=O)CN BWGVNKXGVNDBDI-UHFFFAOYSA-N 0.000 description 3
- 239000003114 blood coagulation factor Substances 0.000 description 3
- 229950003499 fibrin Drugs 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 229940019700 blood coagulation factors Drugs 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 108010049003 Fibrinogen Proteins 0.000 description 1
- 102000008946 Fibrinogen Human genes 0.000 description 1
- 208000007536 Thrombosis Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000003146 anticoagulant agent Substances 0.000 description 1
- 229940127219 anticoagulant drug Drugs 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 229940012952 fibrinogen Drugs 0.000 description 1
- 208000014951 hematologic disease Diseases 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000007639 printing Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/86—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving blood coagulating time or factors, or their receptors
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Hematology (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Urology & Nephrology (AREA)
- Biotechnology (AREA)
- Microbiology (AREA)
- Cell Biology (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
本发明公开一种凝血时间判定方法、电子设备及存储介质,方法包括:获取凝血测试中随时间采样的多个血液凝固特性采样数据;基于所述血液凝固特性采样数据确定凝血开始时间和凝血结束时间;根据凝血开始时间和凝血结束时间确定凝血时间。本发明通过数据分析来确定凝血时间,有效提高了凝血时间的判断效率。
Description
技术领域
本发明涉及凝血测试相关技术领域,特别是一种凝血时间判定方法、电子设备及存储介质。
背景技术
血液从液态变为胶冻状凝块的过程称为血液凝固过程。此过程中血液中的凝血因子按一定顺序激活,使血液中的纤维蛋白原转变为不溶的纤维蛋白。纤维蛋白交联成网,形成胶冻状凝块。血液凝固的时间称为凝血时间,凝血时间的长短受凝血因子的影响。目前有效的凝血时间测量包括凝血酶原时间(Prothrombin Time,PT),活化部分凝血活酶时间(Activated Partial Thromboplastin Time,APTT),活化凝血时间(Activated ClottingTime,ACT)和凝血酶时间(Thrombin Time,TT)。凝血时间的测定对血栓等血液病的监测,抗凝剂的监测,凝血因子的筛查等至关重要。
现有的基于微流控芯片的检验血液凝固的设备和方法,主要通过在微流控芯片中对血样施加电压差,测量一定持续时间的电信号,并通过对测量结果曲线进行分析,得到与凝血时间相关的测量函数的累计特性。
然而现有技术,存在以下的缺点:
一,对输入数据的要求高。其凝血时间的判断基于几个点之间的面积或电信号强度计算,然而在实际情况下由于噪音等的影响,判断凝血时间的几个点不是很明显;
二,实时鲁棒性差。基于缺点一,当原始数据噪声大,现有算法无法做出正确判断;
三,运算速度慢。其部分后处理过程基于面积作为判断依据,因此运算量大,导致算法速度慢,很难满足实时***要求。
发明内容
基于此,有必要针对现有技术对凝血时间的判断准确性低的技术问题,提供一种凝血时间判定方法、电子设备及存储介质。
本发明提供一种凝血时间判定方法,包括:
获取凝血测试中随时间采样的多个血液凝固特性采样数据;
基于所述血液凝固特性采样数据确定凝血开始时间和凝血结束时间;
根据凝血开始时间和凝血结束时间确定凝血时间。
进一步的,所述基于所述血液凝固特性采样数据确定凝血开始时间和凝血结束时间,具体包括:基于搜索步长对每个所述血液凝固特性采样数据计算对应的斜率;
基于所述斜率确定凝血开始时间和凝血结束时间。
更进一步的,所述搜索步长包括起点搜索步长、以及终点搜索步长;
所述基于搜索步长对每个所述血液凝固特性采样数据计算对应的斜率,具体包括:
基于所述起点搜索步长对每个所述血液凝固特性采样数据计算对应的起点搜索斜率,基于所述终点搜索步长对每个所述血液凝固特性采样数据计算对应的终点搜索斜率;
基于所述斜率确定凝血开始时间和凝血结束时间,具体包括:
基于所述起点搜索斜率确定凝血开始时间,基于所述终点搜索斜率确定凝血结束时间。
再进一步的,所述终点搜索步长采用如下方式确定:
使用多个拟合项数对所述血液凝固特性采样数据进行多次多项式曲线拟合得到多个以拟合项数为参数的拟合曲线;
根据所述凝血测试类型,从多个所述拟合曲线中选择符合预设选择规则的拟合曲线作为参考拟合曲线,将所述参考拟合曲线的拟合项数作为参考拟合项数;
根据所述参考拟合项数计算所述终点搜索步长。
再进一步的,所述选择规则包括:
如果所述凝血测试类型为活化部分凝血活酶时间测试、或活化凝血时间测试,则选择拟合误差最大的拟合曲线作为参考拟合曲线;
如果所述凝血测试类型为凝血酶时间测试、或凝血酶原时间测试,则选择拟合误差最小的拟合曲线作为参考拟合曲线。
再进一步的,所述基于所述起点搜索斜率确定凝血开始时间,具体包括:
选择起点搜索斜率最大的血液凝固特性采样数据对应的采样时间作为凝血开始时间。
再进一步的,所述基于所述终点搜索斜率确定凝血结束时间,具体包括:
在所述凝血开始时间之后预设峰干扰时间段后,基于所述终点搜索斜率确定连续相同斜率时间段、或负斜率时间段;
选择所述连续相同斜率时间段的起始时间、或所述负斜率时间段的起始时间作为所述凝血结束时间。
再进一步的,所述基于所述终点搜索斜率确定连续相同斜率时间段,具体包括:
如果有连续多个终点搜索斜率一致且为正斜率的血液凝固特性采样数据,且数量达到预设第一数量阈值,则将终点搜索斜率一致的血液凝固特性采样数据对应的采样时间段作为连续相同斜率时间段,或者;
如果有连续多个终点搜索斜率为零的血液凝固特性采样数据,且数量达到预设第二数量阈值,则将终点搜索斜率为零的血液凝固特性采样数据对应的采样时间段作为连续相同斜率时间段;
其中,所述终点搜索斜率一致,具体为:终点搜索斜率的差值在预设差值阈值范围内。
再进一步的,所述选择所述连续相同斜率时间段的起始时间、或所述负斜率时间段的起始时间作为所述凝血结束时间,具体包括:
计算所述负斜率时间段的起始时间与所述凝血开始时间的第一差值,如果所述第一差值在预设第一时间范围内,则选择负斜率时间段的起始时间作为所述凝血结束时间,否则;
计算所述连续相同斜率时间段的起始时间与所述凝血开始时间的第二差值,如果所述第二差值在预设第二时间范围外,且所述连续相同斜率时间段的起始时间小于所述负斜率时间段的起始时间,则选择负斜率时间段的起始时间作为所述凝血结束时间,如果所述第二差值在预设第二时间范围外,且所述连续相同斜率时间段的起始时间大于所述负斜率时间段的起始时间,则选择连续相同斜率时间段的起始时间作为所述凝血结束时间,否则;
如果所述第二差值在预设第二时间范围内,则选择连续相同斜率时间段的起始时间作为所述凝血结束时间。
本发明提供一种电子设备,包括:
至少一个处理器;以及,
与所述至少一个处理器通信连接的存储器;其中,
所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够:
获取凝血测试中随时间采样的多个血液凝固特性采样数据;
基于所述血液凝固特性采样数据确定凝血开始时间和凝血结束时间;
根据凝血开始时间和凝血结束时间确定凝血时间。
进一步的,所述基于所述血液凝固特性采样数据确定凝血开始时间和凝血结束时间,具体包括:基于搜索步长对每个所述血液凝固特性采样数据计算对应的斜率;
基于所述斜率确定凝血开始时间和凝血结束时间。
更进一步的,所述搜索步长包括起点搜索步长、以及终点搜索步长;
所述基于搜索步长对每个所述血液凝固特性采样数据计算对应的斜率,具体包括:
基于所述起点搜索步长对每个所述血液凝固特性采样数据计算对应的起点搜索斜率,基于所述终点搜索步长对每个所述血液凝固特性采样数据计算对应的终点搜索斜率;
基于所述斜率确定凝血开始时间和凝血结束时间,具体包括:
基于所述起点搜索斜率确定凝血开始时间,基于所述终点搜索斜率确定凝血结束时间。
再进一步的,所述终点搜索步长采用如下方式确定:
使用多个拟合项数对所述血液凝固特性采样数据进行多次多项式曲线拟合得到多个以拟合项数为参数的拟合曲线;
根据所述凝血测试类型,从多个所述拟合曲线中选择符合预设选择规则的拟合曲线作为参考拟合曲线,将所述参考拟合曲线的拟合项数作为参考拟合项数;
根据所述参考拟合项数计算所述终点搜索步长。
再进一步的,所述选择规则包括:
如果所述凝血测试类型为活化部分凝血活酶时间测试、或活化凝血时间测试,则选择拟合误差最大的拟合曲线作为参考拟合曲线;
如果所述凝血测试类型为凝血酶时间测试、或凝血酶原时间测试,则选择拟合误差最小的拟合曲线作为参考拟合曲线。
再进一步的,所述基于所述起点搜索斜率确定凝血开始时间,具体包括:
选择起点搜索斜率最大的血液凝固特性采样数据对应的采样时间作为凝血开始时间。
再进一步的,所述基于所述终点搜索斜率确定凝血结束时间,具体包括:
在所述凝血开始时间之后预设峰干扰时间段后,基于所述终点搜索斜率确定连续相同斜率时间段、或负斜率时间段;
选择所述连续相同斜率时间段的起始时间、或所述负斜率时间段的起始时间作为所述凝血结束时间。
再进一步的,所述基于所述终点搜索斜率确定连续相同斜率时间段,具体包括:
如果有连续多个终点搜索斜率一致且为正斜率的血液凝固特性采样数据,且数量达到预设第一数量阈值,则将终点搜索斜率一致的血液凝固特性采样数据对应的采样时间段作为连续相同斜率时间段,或者;
如果有连续多个终点搜索斜率为零的血液凝固特性采样数据,且数量达到预设第二数量阈值,则将终点搜索斜率为零的血液凝固特性采样数据对应的采样时间段作为连续相同斜率时间段;
其中,所述终点搜索斜率一致,具体为:终点搜索斜率的差值在预设差值阈值范围内。
再进一步的,所述选择所述连续相同斜率时间段的起始时间、或所述负斜率时间段的起始时间作为所述凝血结束时间,具体包括:
计算所述负斜率时间段的起始时间与所述凝血开始时间的第一差值,如果所述第一差值在预设第一时间范围内,则选择负斜率时间段的起始时间作为所述凝血结束时间,否则;
计算所述连续相同斜率时间段的起始时间与所述凝血开始时间的第二差值,如果所述第二差值在预设第二时间范围外,且所述连续相同斜率时间段的起始时间小于所述负斜率时间段的起始时间,则选择负斜率时间段的起始时间作为所述凝血结束时间,如果所述第二差值在预设第二时间范围外,且所述连续相同斜率时间段的起始时间大于所述负斜率时间段的起始时间,则选择连续相同斜率时间段的起始时间作为所述凝血结束时间,否则;
如果所述第二差值在预设第二时间范围内,则选择连续相同斜率时间段的起始时间作为所述凝血结束时间。
本发明提供一种存储介质,所述存储介质存储计算机指令,当计算机执行所述计算机指令时,用于执行如前所述的凝血时间判定方法的所有步骤。
本发明通过数据分析来确定凝血时间,有效提高了凝血时间的判断效率。
附图说明
图1为本发明一种凝血时间判定方法的工作流程图;
图2a为凝血酶原时间测试的采样数据曲线示意图;
图2b为活化凝血时间测试的采样数据曲线示意图;
图2c为活化部分凝血活酶时间测试的采样数据曲线示意图;
图2d为凝血酶时间测试的采样数据曲线示意图;
图3为凝血曲线示意图;
图4a为凝血开始时间搜索示意图;
图4b为根据连续相同斜率时间段判断凝血结束时间示意图;
图4c为根据连续零斜率时间段判断凝血结束时间示意图;
图4d为根据负斜率时间段判断凝血结束时间示意图;
图5为本发明最佳实施例一种凝血时间判定方法的工作流程图;
图6为本发明最佳实施例基于连续斜率时间和负斜率时间得出凝血时间的工作流程图;
图7为本发明一种电子设备的硬件结构示意图;
图8为本发明最佳实施例的***结构示意图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细的说明。
如图1所示为本发明一种凝血时间判定方法的工作流程图,包括:
步骤S101,获取凝血测试中随时间采样的多个血液凝固特性采样数据;
步骤S102,基于所述血液凝固特性采样数据确定凝血开始时间和凝血结束时间;
步骤S103,根据凝血开始时间和凝血结束时间确定凝血时间。
具体来说,血液凝固特性采样数据,可以是如图2a~2d所示的凝血酶原时间、活化凝血时间、活化部分凝血活酶时间、凝血酶时间等测试中的采样数据曲线示意图。血液凝固特性采样数据可以通过凝血仪表获取,凝血仪表可以对血样施加电信号,并在每个采样时间进行采样得到血液凝固特性采样数据。得到的血液凝固特性采样数据可以以文件形式保存,例如以二进制文件形式保存。
步骤S101可以获取保存血液凝固特性采样数据的文件,对文件进行读取,得到随时间采样的多个血液凝固特性采样数据。血液凝固特性采样数据的文件可以从实时血样特征分析和记录仪器中获取。具体来说,实时血样采集装置根据血液的毛细运动特点将采集的血液样本均匀散播到分析通道中,实时血样特征分析和记录仪器通过测量实时血样采集装置中分析通道的电气特性进行凝血信号记录。电气特性的测量通过间距在3mm-5mm的两个电极进行。两电极间载有凝血实验试剂。。在两电极间施加电势差,记录两电极间电阻,当血液充满分析通道后,电阻开始发生变化。测量分辨率在0.000001Ω,测量频率在每秒3到9次,并且在存储介质中记录这些测量数据。测定一定时间的电气特性后,可以得到电气特性-时间数据流。
步骤S102可以基于,如血液凝固特性采样数据的曲线的第一次变化点、曲线的波峰、或者根据经验设置阈值判断等确定凝血开始时间。
在其中一个实施例中,所述基于所述血液凝固特性采样数据确定凝血开始时间和凝血结束时间,具体包括:基于搜索步长对每个所述血液凝固特性采样数据计算对应的斜率;
基于所述斜率确定凝血开始时间和凝血结束时间。具体来说,可以将所获取的血液凝固特性采样数据进行斜率计算,并基于每个血液凝固特性采样数据对应的斜率来确定凝血开始时间和凝血结束时间。
斜率可以采用如下方式计算得到:
对每个所述血液凝固特性采样数据:
将当前血液凝固特性采样数据作为第一血液凝固特性采样数据,将第一血液凝固特性采样数据对应的采样时间作为当前采样时间,此处的血液凝固特性采样数据为测试得到的电气特性数据;
将当前采样时间增加所述搜索步长后的采样时间作为待计算采样时间,将所述待计算采样时间对应的血液凝固特性采样数据作为第二血液凝固特性采样数据;
计算第一血液凝固特性采样数据对应的搜索斜率为:(第一血液凝固特性采样数据-第二血液凝固特性采样数据)/搜索步长。
图3示意了一个典型的凝血曲线。凝血曲线主要有五个阶段:
第一阶段:血液滴入试剂卡前,由于两电极之间没有液体的联通,因此具有较大的电阻。
第二阶段:血液滴入试剂卡,由于血浆充满试剂槽,造成两电极间电阻的急剧降低。
第三阶段:血液中的红细胞均匀地充满试剂卡,造成两电极间的电阻升高。
第四阶段:试剂槽中后部的红细胞因为毛细作用缓慢向试剂槽前端移动,造成两电极间电阻降低。
第五阶段:由于凝血瀑布反应的发生,血液凝固阻碍了红细胞的向前移动,在此点后两电极间的电阻变化较小。
从图3可以看出,在不同阶段凝血曲线的斜率都有显著变化,例如,当血浆充满试剂槽时,血液凝固特性采样数据对应的斜率会发生特定变化,可以通过该变化来判断凝血开始时间。而当红细胞缓慢向前移动达到一定时间发生凝固时,血液凝固特性采样数据对应的斜率也会发生特定变化,可以通过该变化来判断凝血结束时间。
然后,步骤S103可以采用凝血结束时间与凝血开始时间的差值作为凝血时间。
在其中一个实施例中,所述搜索步长包括起点搜索步长、以及终点搜索步长;
所述基于搜索步长对每个所述血液凝固特性采样数据计算对应的斜率,具体包括:
基于所述起点搜索步长对每个所述血液凝固特性采样数据计算对应的起点搜索斜率,基于所述终点搜索步长对每个所述血液凝固特性采样数据计算对应的终点搜索斜率;
基于所述斜率确定凝血开始时间和凝血结束时间,具体包括:
基于所述起点搜索斜率确定凝血开始时间,基于所述终点搜索斜率确定凝血结束时间。
血液凝固特性采样数据需要进行数据的预处理。预处理包括步长的计算。步长计算可以有很多种方法,包括但不限于通过曲线拟合抽象结果进行估计。步长值也可以通过预设值或预设值和曲线拟合结果的综合加权得出。本实施例对起点和结束采用不同的搜索步长,从而可以采用不同的斜率来分别判断凝血开始时间和凝血结束时间。起点搜索步长可以采用预设值,优选为两点,即每个血液凝固特性采样数据的起点搜索斜率优选为(第一血液凝固特性采样数据-两点后的血液凝固特性采样数据)/2。
在其中一个实施例中,所述起点搜索步长为预设值。
在其中一个实施例中,所述终点搜索步长采用如下方式确定:
使用多个拟合项数对所述血液凝固特性采样数据进行多次多项式曲线拟合得到多个以拟合项数为参数的拟合曲线;
根据所述凝血测试类型,从多个所述拟合曲线中选择符合预设选择规则的拟合曲线作为参考拟合曲线,将所述参考拟合曲线的拟合项数作为参考拟合项数;
根据所述参考拟合项数计算所述终点搜索步长。
可以采用现有的多项式曲线拟合方式对血液凝固特性采样数据进行曲线拟合从而得到拟合曲线。拟合时尝试不同的拟合项数,然后根据选择规则确定出拟合项数,计算得到终点搜索步长。通过拟合项数求步长有多种方法,本实施例认为拟合项数与步长成正比关系,即终点搜索步长=参考拟合项数/预设值。
在其中一个实施例中,所述选择规则包括:
如果所述凝血测试类型为活化部分凝血活酶时间测试、或活化凝血时间测试,则选择拟合误差最大的拟合曲线作为参考拟合曲线;
如果所述凝血测试类型为凝血酶时间测试、或凝血酶原时间测试,则选择拟合误差最小的拟合曲线作为参考拟合曲线。
本实施例对于不同的凝血测试采用不同的拟合项数作为步长计算依据,例如对于数据量较大的曲线如APTT/ACT使用拟合误差最大的拟合曲线的拟合项数计算得到终点搜索步长,因为拟合误差最大的拟合项数能够更大限度地反映曲线的趋势并消除噪音影响;对于数据量较小的曲线如TT/PT使用拟合误差最小的拟合曲线的拟合项数计算得到终点搜索步长,因为拟合误差最小的拟合项数能够更多地保留数据量较小的曲线的细节。
预处理后开始进行数据的二次分析,包括凝血开始时间和凝血结束时间的搜索。
在其中一个实施例中,所述基于所述起点搜索斜率确定凝血开始时间,具体包括:
选择起点搜索斜率最大的血液凝固特性采样数据对应的采样时间作为凝血开始时间。
凝血开始时间搜索的主要依据是凝血曲线在低维空间表达的第一次曲线变化,例如低维空间测量斜率变化是最常用的一种依据,但是寻找凝血开始时间的方法包括但是不限于通过斜率的判断,包括任何其它基于低维空间特征得到的曲线变化点。本实施例对凝血开始时间搜索是基于起点搜索步长计算起点搜索斜率来确定凝血开始时间。起点搜索步长可以采用预设值,优选为两点,即每个血液凝固特性采样数据的起点搜索斜率优选为(第一血液凝固特性采样数据-两点后的血液凝固特性采样数据)/2。从图4a可以看出,在第二阶段血浆填充试剂槽时,其斜率持续增长,然后斜率增长到最高点后反向,该最高点对应的采样时间即为凝血开始时间。因此,本实施例采用起点搜索斜率最大的血液凝固特性采样数据对应的采样时间作为凝血开始时间。
在其中一个实施例中,所述基于所述终点搜索斜率确定凝血结束时间,具体包括:
在所述凝血开始时间之后预设峰干扰时间段后,基于所述终点搜索斜率确定连续相同斜率时间段、或负斜率时间段;
以所述连续相同斜率时间段的起始时间、或所述负斜率时间段的起始时间作为所述凝血结束时间。
本实施例进行凝血结束时间的搜索,在这一步中凝血数据需要进行一次从高维空间到低维空间的特征提取,即通过对低维空间特征曲线斜率进行分析找到凝血结束时间。具体包括:
首先将凝血结束时间搜索起点后移。如图3所示,在第二阶段凝血开始之后,由于第三阶段和第四阶段曲线会形成一个曲线峰。峰的出现会干扰第五阶段凝血结束时间的寻找。因此,需要在峰干扰时间段后开始搜索凝血结束时间,避免峰干扰。峰干扰时间段可以采用不同凝血测试的预设值或预设值和曲线数据量的综合加权得出。接下来基于先前求得的步长在凝血结束时间搜索起点后计算血液凝固特性采样数据点的低维空间特征曲线斜率,即每个血液凝固特性采样数据的终点搜索斜率为(第一血液凝固特性采样数据-步长点后的血液凝固特性采样数据)/步长。
本发明提供的基于低维空间特征曲线斜率估计的血凝结束点有几种方法:1)恒定数据斜率:当数据变化在持续,但是数据斜率是恒定的。这表明数据曲线其实处于一种平稳状态,这个数据段和数据点就可能是凝血结束时间;2)零数据斜率:最低数据变化率是低维度特征变化最慢最低点,通常表明数据在这个时间窗口没有发生变化。这种变化的最常用的测量值是低维度特征的斜率,所以这时候斜率为零并且保持零斜率一段时间窗口的点,就是可能的凝血结束时间;3)基于1,2的平滑窗口的恒定数据变化率:对于噪声相对较大的数据段,虽然数据变化在低维空间呈现锯齿状,但其实这是一种数据稳定无变化的特征。也是一种可能的凝血结束时间;4)负数据斜率:在低维特征空间的数据变化负斜率,这表明数据变化的转折,这在实际中也是可能的凝血结束时间;5)最明显数据变化点:在某些实际凝血数据里,最明显的数据变化点,如斜率极值,可能表明凝血结束。在这一部分,可以根据经验值设置凝血结束点估计方法的优先级,也可以设置不同的算法跳出策略,在找到第一个可行度较高的凝血结束时间时候就跳出后续算法处理,这样以达到快速实时的要求。
因此,本实施例的凝血结束时间的寻找包括连续斜率点的寻找以及负斜率点的寻找。其中连续斜率点的寻找包括连续正斜率点的寻找和连续零斜率点的寻找。连续斜率点可以被视为凝血结束时间,是因为凝固反应发生后红细胞运动停止,导致电阻变化趋于平稳。负斜率点可以被视为凝血结束时间,是因为当凝血发生后纤维蛋白发生交联反应形成网状结构,导致电阻升高,反映到斜率上为负斜率。
在其中一个实施例中,所述基于所述终点搜索斜率确定连续相同斜率时间段,具体包括:
如果有连续多个终点搜索斜率一致且为正斜率的血液凝固特性采样数据,且数量达到预设第一数量阈值,则将终点搜索斜率一致的血液凝固特性采样数据对应的时间段作为连续相同斜率时间段,或者;
如果有连续多个终点搜索斜率为零的血液凝固特性采样数据,且数量达到预设第二数量阈值,则将终点搜索斜率为零的血液凝固特性采样数据对应的时间段作为连续相同斜率时间段;
其中,所述终点搜索斜率一致,具体为:终点搜索斜率的差值在预设差值阈值范围内。
具体来说,如有连续相同终点搜索斜率且为正斜率的点,且数量达到第一数量阈值,例如:50,则记录参数,停止连续相同斜率点搜索,如果出现了连续相同终点搜索斜率的点但数量小于第一数量阈值,则返回检查前面点的终点搜索斜率是否有少量波动,如果有的话可以忽略波动并计数,达到阈值停止搜索;如果没有的话清空计数继续向前搜索。
而第二数量阈值则用于搜索连续多个终点搜索斜率为零的血液凝固特性采样数据,由于终点搜索斜率为零的血液凝固特性采样数据更能反映凝血结束,因此,第二数量阈值可以小于第一数量阈值,优选地,第一数量阈值为50,第二数量阈值为15。
对于出现连续斜率的情况,如图4b和4c所示,凝血结束时间可以参考所记录参数,对应阈值以及步长共同得出。例如凝血结束时间参数=(所记录参数-阈值+步长),凝血结束时间为该参数对应时间,即连续相同斜率时间段的起始时间。
对于负斜率情况,凝血结束时间为记录参数所对应时间,即负斜率点出现的时间。
在其中一个实施例中,所述选择所述连续相同斜率时间段的起始时间、或所述负斜率时间段的起始时间作为所述凝血结束时间,具体包括:
计算所述负斜率时间段的起始时间与所述凝血开始时间的第一差值,如果所述第一差值在预设第一时间范围内,则选择负斜率时间段的起始时间作为所述凝血结束时间,否则;
计算所述连续相同斜率时间段的起始时间与所述凝血开始时间的第二差值,如果所述第二差值在预设第二时间范围外,且所述连续相同斜率时间段的起始时间小于所述负斜率时间段的起始时间,则选择负斜率时间段的起始时间作为所述凝血结束时间,如果所述第二差值在预设第二时间范围外,且所述连续相同斜率时间段的起始时间大于所述负斜率时间段的起始时间,则选择连续相同斜率时间段的起始时间作为所述凝血结束时间,否则;
如果所述第二差值在预设第二时间范围内,则选择连续相同斜率时间段的起始时间作为所述凝血结束时间。
本实施例对负斜率和相同斜率分别判断,以选择合适的凝血结束时间。由于终点搜索斜率为负的血液凝固特性采样数据更能反映凝血结束,因此,如图4d所示,当负斜率点出现时,记录出现的参数。
如图5所示为本发明最佳实施例一种凝血时间判定方法的工作流程图,包括:
步骤S501,读取包括血液凝固特性采样数据的txt文件,并设置参数;
步骤S502,检查是否是新文件,如果是则执行步骤S503,否则给文件编号,重置时间,执行步骤S503;
步骤S503,用多项式拟合功能拟合曲线,尝试不同项数的多项式拟合,寻找拟合误差最大的拟合项数,然后用这个项数计算出终点搜索步长;
步骤S504,通过计算每一血液凝固特性采样数据到两点后的血液凝固特性采样数据的斜率作为起点搜索斜率,起点搜索斜率最大的点为凝血起始点。
步骤S505,设置参数,清空变量为下一步做准备;
步骤S506,将凝血结束时间搜索起始指数调后150,避免峰干扰;
步骤S507,在凝血结束时间搜索起始指数位置开始,计算每一血液凝固特性采样数据到与终点搜索步长距离后的血液凝固特性采样数据的斜率作为终点搜索斜率;
步骤S508,如有连续终点搜索斜率相同的血液凝固特性采样数据,且数量达到50,则记录最后一个终点搜索斜率相同的血液凝固特性采样数据的采样时间为连续斜率时间,停止连续相同斜率点搜索。
步骤S509,如果出现了连续终点搜索斜率相同的血液凝固特性采样数据,但数量在5-50间,返回检查前面血液凝固特性采样数据的终点搜索斜率是否有少量波动,如果有,则可以忽略波动并计数,达到阈值则记录最后一个终点搜索斜率相同的血液凝固特性采样数据的采样时间为连续斜率的结束时间,停止搜索,如果没有,则清空计数继续向前搜索;
步骤S510,如果出现了连续终点搜索斜率为零的血液凝固特性采样数据,且数量达到15,则记录最后一个终点搜索斜率为零的血液凝固特性采样数据的采样时间为连续斜率结束时间,停止搜索;
步骤S511,如果出现了终点搜索斜率为负数的血液凝固特性采样数据,记录第一个出现终点搜索斜率为负数的血液凝固特性采样数据的采样时间为负斜率时间;
步骤S512,基于连续斜率的起始时间和负斜率时间得出各通道的凝血时间;
步骤S513,确定凝血时间值并对本次数据进行错误排查。例如通过上述数据分析筛选决定的最优凝血时间估计在经验允许的范围内,那么这个凝血时间将会显示给用户。但是如果通过上述数据分析筛选决定的最优凝血时间估计处于不正常的经验允许值范围之内,***显示估计失败给用户。
如图6所示为本发明最佳实施例基于连续斜率时间和负斜率时间得出凝血时间的工作流程图,包括:
步骤S601,计算负斜率时间-凝血开始时间,如果在(100,300)内,则以负斜率时间作为凝血结束时间,执行步骤S604,否则执行步骤S602,其中,100、300为数据点的数量,其适用数据一秒采集三个点,换算成时间为33秒、100秒;
步骤S602,计算连续斜率起始时间-凝血开始时间,如果在(100,300)内,则以连续斜率时间作为凝血结束时间,执行步骤S604,否则执行步骤S603;
步骤S603,计算连续斜率时间-负斜率时间,如果小于或等于0,则以负斜率时间作为凝血结束时间,执行步骤S604,否则以连续斜率时间作为凝血结束时间,执行步骤S604;
步骤S604,显示凝血时间为凝血结束时间-凝血开始时间。
如图7所示为本发明一种电子设备1的硬件结构示意图,包括:
至少一个处理器701;以及,
与所述至少一个处理器701通信连接的存储器702;其中,
所述存储器702存储有可被所述一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够:
获取凝血测试中随时间采样的多个血液凝固特性采样数据;
基于所述血液凝固特性采样数据确定凝血开始时间和凝血结束时间;
根据凝血开始时间和凝血结束时间确定凝血时间。
图7中以一个处理器702为例。
电子设备还可以包括:输入装置703和输出装置704。
处理器701、存储器702、输入装置703及显示装置704可以通过总线或者其他方式连接,图中以通过总线连接为例。
处理器701可以是通用处理器(CPU),也可以是专用数字信号处理芯片(DSP),也可以是FPGA。
存储器702作为一种非易失性计算机可读存储介质,可用于存储非易失性软件程序、非易失性计算机可执行程序以及模块,如本申请实施例中的凝血时间判定方法对应的程序指令/模块,例如,图1、图5、图6所示的方法流程。处理器701通过运行存储在存储器702中的非易失性软件程序、指令以及模块,从而执行各种功能应用以及数据处理,即实现上述实施例中的凝血时间判定方法。
存储器702可以包括存储程序区和存储数据区,其中,存储程序区可存储操作***、至少一个功能所需要的应用程序;存储数据区可存储根据凝血时间判定方法的使用所创建的数据,例如***运算的中间结果和最后凝血数据等。此外,存储器702可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实施例中,存储器702可选包括相对于处理器701远程设置的存储器,这些远程存储器可以通过网络连接至执行凝血时间判定方法的装置。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入装置703可接收输入的用户点击,以及产生与凝血时间判定方法的用户设置以及功能控制有关的信号输入。显示装置704可包括显示屏等显示设备。
在所述一个或者多个模块存储在所述存储器702中,当被所述一个或者多个处理器701运行时,执行上述任意方法实施例中的凝血时间判定方法。
在其中一个实施例中,所述基于所述血液凝固特性采样数据确定凝血开始时间和凝血结束时间,具体包括:基于搜索步长对每个所述血液凝固特性采样数据计算对应的斜率;
基于所述斜率确定凝血开始时间和凝血结束时间。
在其中一个实施例中,所述搜索步长包括起点搜索步长、以及终点搜索步长;
所述基于搜索步长对每个所述血液凝固特性采样数据计算对应的斜率,具体包括:
基于所述起点搜索步长对每个所述血液凝固特性采样数据计算对应的起点搜索斜率,基于所述终点搜索步长对每个所述血液凝固特性采样数据计算对应的终点搜索斜率;
基于所述斜率确定凝血开始时间和凝血结束时间,具体包括:
基于所述起点搜索斜率确定凝血开始时间,基于所述终点搜索斜率确定凝血结束时间。
在其中一个实施例中,所述起点搜索步长为预设值。
在其中一个实施例中,所述终点搜索步长采用如下方式确定:
使用多个拟合项数对所述血液凝固特性采样数据进行多次多项式曲线拟合得到多个以拟合项数为参数的拟合曲线;
根据所述凝血测试类型,从多个所述拟合曲线中选择符合预设选择规则的拟合曲线作为参考拟合曲线,将所述参考拟合曲线的拟合项数作为参考拟合项数;
根据所述参考拟合项数计算所述终点搜索步长。
在其中一个实施例中,所述选择规则包括:
如果所述凝血测试类型为活化部分凝血活酶时间测试、或活化凝血时间测试,则选择拟合误差最大的拟合曲线作为参考拟合曲线;
如果所述凝血测试类型为凝血酶时间测试、或凝血酶原时间测试,则选择拟合误差最小的拟合曲线作为参考拟合曲线。
在其中一个实施例中,所述基于所述起点搜索斜率确定凝血开始时间,具体包括:
选择起点搜索斜率最大的血液凝固特性采样数据对应的采样时间作为凝血开始时间。
在其中一个实施例中,所述基于所述终点搜索斜率确定凝血结束时间,具体包括:
在所述凝血开始时间之后预设峰干扰时间段后,基于所述终点搜索斜率确定连续相同斜率时间段、或负斜率时间段;
选择所述连续相同斜率时间段的起始时间、或所述负斜率时间段的起始时间作为所述凝血结束时间。
在其中一个实施例中,所述基于所述终点搜索斜率确定连续相同斜率时间段,具体包括:
如果有连续多个终点搜索斜率一致且为正斜率的血液凝固特性采样数据,且数量达到预设第一数量阈值,则将终点搜索斜率一致的血液凝固特性采样数据对应的采样时间段作为连续相同斜率时间段,或者;
如果有连续多个终点搜索斜率为零的血液凝固特性采样数据,且数量达到预设第二数量阈值,则将终点搜索斜率为零的血液凝固特性采样数据对应的采样时间段作为连续相同斜率时间段;
其中,所述终点搜索斜率一致,具体为:终点搜索斜率的差值在预设差值阈值范围内。
在其中一个实施例中,所述选择所述连续相同斜率时间段的起始时间、或所述负斜率时间段的起始时间作为所述凝血结束时间,具体包括:
计算所述负斜率时间段的起始时间与所述凝血开始时间的第一差值,如果所述第一差值在预设第一时间范围内,则选择负斜率时间段的起始时间作为所述凝血结束时间,否则;
计算所述连续相同斜率时间段的起始时间与所述凝血开始时间的第二差值,如果所述第二差值在预设第二时间范围外,且所述连续相同斜率时间段的起始时间小于所述负斜率时间段的起始时间,则选择负斜率时间段的起始时间作为所述凝血结束时间,如果所述第二差值在预设第二时间范围外,且所述连续相同斜率时间段的起始时间大于所述负斜率时间段的起始时间,则选择连续相同斜率时间段的起始时间作为所述凝血结束时间,否则;
如果所述第二差值在预设第二时间范围内,则选择连续相同斜率时间段的起始时间作为所述凝血结束时间。
作为本发明最佳实施例,图8显示了一种从实时血液采集装置中采集并处理凝血数据的***示意架构。
***包括如图7所示的电子设备1、以及可以记录凝血过程数据的实时血液数据特性分析记录仪2。这种分析记录仪也可以是离线记录数据仪器,***使用移动或固定电源6供电。
图8所示***可以实施一个实时并行历史数据分析模块3,这个模块用来提高具体样本的实时凝血估计的性能。并行历史数据分析模块可以是同一个病人的历史数据,也可以基于大样本的同类病人群体的历史数据。
图8所示***可以实施一个无线物联网通信模块4与其它同类设备进行本地局域网通信,用来实施大规模大数据的凝血估计分析。这个无线物联网通信模块传输的数据可以输入实时并行历史数据分析模块供其综合分析。无线物联网通信模块也可以向上接入广域网的元计算或者云存储平台,由这些云端资源进行后台分析。
图8所示***可以实施连接打印设备5,以打印实时中间数据和任何凝血数据。同样设计可以使用显示器来显示实时中间数据和任何凝血数据。
本发明提供一种存储介质,所述存储介质存储计算机指令,当计算机执行所述计算机指令时,用于执行如前所述的凝血时间判定方法的所有步骤。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (13)
1.一种凝血时间判定方法,其特征在于,包括:
获取凝血测试中随时间采样的多个血液凝固特性采样数据;
基于所述血液凝固特性采样数据确定凝血开始时间和凝血结束时间;
根据凝血开始时间和凝血结束时间确定凝血时间;
所述基于所述血液凝固特性采样数据确定凝血开始时间和凝血结束时间,具体包括:基于搜索步长对每个所述血液凝固特性采样数据计算对应的斜率;
基于所述斜率确定凝血开始时间和凝血结束时间;
所述搜索步长包括起点搜索步长、以及终点搜索步长;
所述基于搜索步长对每个所述血液凝固特性采样数据计算对应的斜率,具体包括:
基于所述起点搜索步长对每个所述血液凝固特性采样数据计算对应的起点搜索斜率,基于所述终点搜索步长对每个所述血液凝固特性采样数据计算对应的终点搜索斜率;
基于所述斜率确定凝血开始时间和凝血结束时间,具体包括:
基于所述起点搜索斜率确定凝血开始时间,基于所述终点搜索斜率确定凝血结束时间;
所述起点搜索斜率区别于所述终点搜索斜率;
所述起点搜索步长和所述终点搜索步长不同;
所述起点搜索步长采用预设值;
所述终点搜索步长采用如下方式确定:
使用多个拟合项数对所述血液凝固特性采样数据进行多次多项式曲线拟合得到多个以拟合项数为参数的拟合曲线;
根据凝血测试类型,从多个所述拟合曲线中选择符合预设选择规则的拟合曲线作为参考拟合曲线,将所述参考拟合曲线的拟合项数作为参考拟合项数;
根据所述参考拟合项数计算所述终点搜索步长。
2.根据权利要求1所述的凝血时间判定方法,其特征在于,所述选择规则包括:
如果所述凝血测试类型为活化部分凝血活酶时间测试、或活化凝血时间测试,则选择拟合误差最大的拟合曲线作为参考拟合曲线;
如果所述凝血测试类型为凝血酶时间测试、或凝血酶原时间测试,则选择拟合误差最小的拟合曲线作为参考拟合曲线。
3.根据权利要求1所述的凝血时间判定方法,其特征在于,所述基于所述起点搜索斜率确定凝血开始时间,具体包括:
选择起点搜索斜率最大的血液凝固特性采样数据对应的采样时间作为凝血开始时间。
4.根据权利要求1所述的凝血时间判定方法,其特征在于,所述基于所述终点搜索斜率确定凝血结束时间,具体包括:
在所述凝血开始时间之后预设峰干扰时间段后,基于所述终点搜索斜率确定连续相同斜率时间段、或负斜率时间段;
选择所述连续相同斜率时间段的起始时间、或所述负斜率时间段的起始时间作为所述凝血结束时间。
5.根据权利要求4所述的凝血时间判定方法,其特征在于,所述基于所述终点搜索斜率确定连续相同斜率时间段,具体包括:
如果有连续多个终点搜索斜率一致且为正斜率的血液凝固特性采样数据,且数量达到预设第一数量阈值,则将终点搜索斜率一致的血液凝固特性采样数据对应的采样时间段作为连续相同斜率时间段,或者;
如果有连续多个终点搜索斜率为零的血液凝固特性采样数据,且数量达到预设第二数量阈值,则将终点搜索斜率为零的血液凝固特性采样数据对应的采样时间段作为连续相同斜率时间段;
其中,所述终点搜索斜率一致,具体为:终点搜索斜率的差值在预设差值阈值范围内。
6.根据权利要求4所述的凝血时间判定方法,其特征在于,所述选择所述连续相同斜率时间段的起始时间、或所述负斜率时间段的起始时间作为所述凝血结束时间,具体包括:
计算所述负斜率时间段的起始时间与所述凝血开始时间的第一差值,如果所述第一差值在预设第一时间范围内,则选择负斜率时间段的起始时间作为所述凝血结束时间,否则;
计算所述连续相同斜率时间段的起始时间与所述凝血开始时间的第二差值,如果所述第二差值在预设第二时间范围外,且所述连续相同斜率时间段的起始时间小于所述负斜率时间段的起始时间,则选择负斜率时间段的起始时间作为所述凝血结束时间,如果所述第二差值在预设第二时间范围外,且所述连续相同斜率时间段的起始时间大于所述负斜率时间段的起始时间,则选择连续相同斜率时间段的起始时间作为所述凝血结束时间,否则;
如果所述第二差值在预设第二时间范围内,则选择连续相同斜率时间段的起始时间作为所述凝血结束时间。
7.一种电子设备,其特征在于,包括:
至少一个处理器;以及,
与所述至少一个处理器通信连接的存储器;其中,
所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够:
获取凝血测试中随时间采样的多个血液凝固特性采样数据;
基于所述血液凝固特性采样数据确定凝血开始时间和凝血结束时间;
根据凝血开始时间和凝血结束时间确定凝血时间;
所述基于所述血液凝固特性采样数据确定凝血开始时间和凝血结束时间,具体包括:基于搜索步长对每个所述血液凝固特性采样数据计算对应的斜率;
基于所述斜率确定凝血开始时间和凝血结束时间;
所述搜索步长包括起点搜索步长、以及终点搜索步长;
所述基于搜索步长对每个所述血液凝固特性采样数据计算对应的斜率,具体包括:
基于所述起点搜索步长对每个所述血液凝固特性采样数据计算对应的起点搜索斜率,基于所述终点搜索步长对每个所述血液凝固特性采样数据计算对应的终点搜索斜率;
基于所述斜率确定凝血开始时间和凝血结束时间,具体包括:
基于所述起点搜索斜率确定凝血开始时间,基于所述终点搜索斜率确定凝血结束时间;
所述起点搜索斜率区别于所述终点搜索斜率;
所述起点搜索步长和所述终点搜索步长不同;
所述起点搜索步长采用预设值;
所述终点搜索步长采用如下方式确定:
使用多个拟合项数对所述血液凝固特性采样数据进行多次多项式曲线拟合得到多个以拟合项数为参数的拟合曲线;
根据凝血测试类型,从多个所述拟合曲线中选择符合预设选择规则的拟合曲线作为参考拟合曲线,将所述参考拟合曲线的拟合项数作为参考拟合项数;
根据所述参考拟合项数计算所述终点搜索步长。
8.根据权利要求7所述的电子设备,其特征在于,所述选择规则包括:
如果所述凝血测试类型为活化部分凝血活酶时间测试、或活化凝血时间测试,则选择拟合误差最大的拟合曲线作为参考拟合曲线;
如果所述凝血测试类型为凝血酶时间测试、或凝血酶原时间测试,则选择拟合误差最小的拟合曲线作为参考拟合曲线。
9.根据权利要求7所述的电子设备,其特征在于,所述基于所述起点搜索斜率确定凝血开始时间,具体包括:
选择起点搜索斜率最大的血液凝固特性采样数据对应的采样时间作为凝血开始时间。
10.根据权利要求7所述的电子设备,其特征在于,所述基于所述终点搜索斜率确定凝血结束时间,具体包括:
在所述凝血开始时间之后预设峰干扰时间段后,基于所述终点搜索斜率确定连续相同斜率时间段、或负斜率时间段;
选择所述连续相同斜率时间段的起始时间、或所述负斜率时间段的起始时间作为所述凝血结束时间。
11.根据权利要求10所述的电子设备,其特征在于,所述基于所述终点搜索斜率确定连续相同斜率时间段,具体包括:
如果有连续多个终点搜索斜率一致且为正斜率的血液凝固特性采样数据,且数量达到预设第一数量阈值,则将终点搜索斜率一致的血液凝固特性采样数据对应的采样时间段作为连续相同斜率时间段,或者;
如果有连续多个终点搜索斜率为零的血液凝固特性采样数据,且数量达到预设第二数量阈值,则将终点搜索斜率为零的血液凝固特性采样数据对应的采样时间段作为连续相同斜率时间段;
其中,所述终点搜索斜率一致,具体为:终点搜索斜率的差值在预设差值阈值范围内。
12.根据权利要求10所述的电子设备,其特征在于,所述选择所述连续相同斜率时间段的起始时间、或所述负斜率时间段的起始时间作为所述凝血结束时间,具体包括:
计算所述负斜率时间段的起始时间与所述凝血开始时间的第一差值,如果所述第一差值在预设第一时间范围内,则选择负斜率时间段的起始时间作为所述凝血结束时间,否则;
计算所述连续相同斜率时间段的起始时间与所述凝血开始时间的第二差值,如果所述第二差值在预设第二时间范围外,且所述连续相同斜率时间段的起始时间小于所述负斜率时间段的起始时间,则选择负斜率时间段的起始时间作为所述凝血结束时间,如果所述第二差值在预设第二时间范围外,且所述连续相同斜率时间段的起始时间大于所述负斜率时间段的起始时间,则选择连续相同斜率时间段的起始时间作为所述凝血结束时间,否则;
如果所述第二差值在预设第二时间范围内,则选择连续相同斜率时间段的起始时间作为所述凝血结束时间。
13.一种存储介质,其特征在于,所述存储介质存储计算机指令,当计算机执行所述计算机指令时,用于执行如权利要求1~6任一项所述的凝血时间判定方法的所有步骤。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810554622.7A CN108761105B (zh) | 2018-06-01 | 2018-06-01 | 凝血时间判定方法、电子设备及存储介质 |
PCT/CN2019/080307 WO2019228048A1 (zh) | 2018-06-01 | 2019-03-29 | 凝血时间判定方法、电子设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810554622.7A CN108761105B (zh) | 2018-06-01 | 2018-06-01 | 凝血时间判定方法、电子设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108761105A CN108761105A (zh) | 2018-11-06 |
CN108761105B true CN108761105B (zh) | 2020-12-01 |
Family
ID=64001785
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810554622.7A Active CN108761105B (zh) | 2018-06-01 | 2018-06-01 | 凝血时间判定方法、电子设备及存储介质 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN108761105B (zh) |
WO (1) | WO2019228048A1 (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108761105B (zh) * | 2018-06-01 | 2020-12-01 | 广州万孚生物技术股份有限公司 | 凝血时间判定方法、电子设备及存储介质 |
CN110967306B (zh) * | 2019-11-14 | 2023-11-24 | 迈克医疗电子有限公司 | 反应稳定起始时间确定方法和装置、分析仪器和存储介质 |
CN112161955A (zh) * | 2020-07-22 | 2021-01-01 | 三诺生物传感股份有限公司 | 一种用于凝血四项的测试方法 |
CN113343899B (zh) * | 2021-06-19 | 2022-10-14 | 南京岚煜生物科技有限公司 | 基于斜率检波的确定凝血突变点的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101484805A (zh) * | 2006-04-07 | 2009-07-15 | 埃葛梅崔克斯股份有限公司 | 用于监测液体样品流动特征改变的方法和器械 |
CN103543278A (zh) * | 2012-07-16 | 2014-01-29 | 微点生物科技有限公司 | 确定血液凝固特性 |
CN104034672A (zh) * | 2013-03-06 | 2014-09-10 | 希森美康株式会社 | 血液凝固分析装置及血液凝固分析方法 |
CN104204787A (zh) * | 2012-04-13 | 2014-12-10 | 索尼公司 | 凝血***分析仪、凝血***分析方法及程序 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105008906B (zh) * | 2013-03-13 | 2020-01-07 | 索尼公司 | 血液状态分析装置、血液状态分析*** |
WO2016049552A1 (en) * | 2014-09-26 | 2016-03-31 | Abbott Point Of Care Inc. | Cartridge device identification for coagulation assays in fluid samples |
CN108761105B (zh) * | 2018-06-01 | 2020-12-01 | 广州万孚生物技术股份有限公司 | 凝血时间判定方法、电子设备及存储介质 |
-
2018
- 2018-06-01 CN CN201810554622.7A patent/CN108761105B/zh active Active
-
2019
- 2019-03-29 WO PCT/CN2019/080307 patent/WO2019228048A1/zh active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101484805A (zh) * | 2006-04-07 | 2009-07-15 | 埃葛梅崔克斯股份有限公司 | 用于监测液体样品流动特征改变的方法和器械 |
CN104204787A (zh) * | 2012-04-13 | 2014-12-10 | 索尼公司 | 凝血***分析仪、凝血***分析方法及程序 |
CN103543278A (zh) * | 2012-07-16 | 2014-01-29 | 微点生物科技有限公司 | 确定血液凝固特性 |
CN104034672A (zh) * | 2013-03-06 | 2014-09-10 | 希森美康株式会社 | 血液凝固分析装置及血液凝固分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108761105A (zh) | 2018-11-06 |
WO2019228048A1 (zh) | 2019-12-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108761105B (zh) | 凝血时间判定方法、电子设备及存储介质 | |
WO2021258657A1 (zh) | 一种电池健康状态获取方法及装置、存储介质 | |
CN108571997B (zh) | 一种检测探头稳定接触被测点的方法和装置 | |
CN110059894A (zh) | 设备状态评估方法、装置、***及存储介质 | |
CN106896273A (zh) | 电池单体的内阻检测方法、检测装置及具有其的车辆 | |
CN109034180A (zh) | 异常检测方法、装置、计算机可读存储介质及电子设备 | |
CN109307795A (zh) | 波形点检测方法及装置 | |
CN109447506A (zh) | 一种电能质量基准水平评价方法及*** | |
CN108982941A (zh) | 一种周期频率的实时测量方法及装置、示波器、存储介质 | |
US11438994B2 (en) | Filament current control method and apparatus | |
WO2021047342A1 (zh) | 一种同源双采样方式的行波测距方法、装置及存储介质 | |
CN117368782A (zh) | 数据分析方法、电池soh估算方法及装置 | |
CN115841049B (zh) | 一种基于多源信息融合的设备寿命可靠性评估方法及*** | |
US6799128B2 (en) | Measurement system for sampling a signal and evaluating data representing the signal | |
CN110888100A (zh) | 一种单相智能电能表在线带载检测***及方法 | |
CN106326098A (zh) | 测试页面感知性能的方法与装置 | |
CN115684711A (zh) | 一种电能计量方法、装置、计算机设备及存储介质 | |
CN114664392A (zh) | 电化学参数预测方法、装置、电子设备及可读存储介质 | |
CN109030934A (zh) | 一种提高峰值功率测量速度的方法 | |
CN112241343B (zh) | 一种慢盘检测方法、装置、电子设备及可读存储介质 | |
CN113435664A (zh) | 一种电费异常数据分析方法、装置、终端设备及介质 | |
CN103713189A (zh) | 一种功率潜动快速检测电路 | |
US10156538B2 (en) | Control method and system for test strip electrodes | |
CN117054879A (zh) | 一种电池寿命的预测方法、装置、电子设备及存储介质 | |
CN115291111B (zh) | 电池静置时间预测模型的训练方法以及静置时间预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |