一种凝血酶原时间的检测方法、检测装置及计算机可读存储介质与流程

1.本发明涉及电化学检测技术领域，具体涉及一种凝血酶原时间的检测方法、检测装置及计算机可读存储介质。


背景技术：

2.凝血酶原时间是检查外源凝血系统诸因子及相关抑制物的重要筛选试验，也是口服抗凝剂治疗的主要手段。凝血酶原时间的常规检测方法有光学法、磁路法、电化学法等。但是，光学法易受特异血浆干扰，磁路法易受磁珠质量和杯壁光滑程度的影响；电化学法因方便、快捷、简单、易用等优点，在外科手术和抗凝患者日常血凝监控等方面发挥着重要作用。
3.基于电化学法，凝血酶原时间检测卡的基本原理是通过检测酶反应过程中产生的电信号来测定凝血酶原时间。因此，在采集到参数数据后，需要对其进行平滑和寻峰处理。但是在现有技术中，在进行平滑处理时很难避免脉冲噪声、椒盐噪声的干扰，从而影响凝血酶原时间的检测精度；另外，通过单个特征点计算凝血酶原时间，也会存在较大的误差。
4.综上所述，急需一种凝血酶原时间的检测方法，能够有效避免电化学方法中噪声对检测精度的影响，提高检测的准确性和可靠性。


技术实现要素：

5.本发明目的在于提供一种凝血酶原时间的检测方法，能够避免噪声对凝血酶原时间检测精度的干扰，提高检测的准确性和可靠性，具体技术方案如下：一种凝血酶原时间的检测方法，具体包括如下步骤：数据采集：将待测样本滴入凝血酶原时间电化学试纸的进样口，按照预设时间间隔对原始数据进行采集，得到参数数据m；所述参数数据m为阻抗-时间曲线或者电流-时间曲线上的数据；数据处理：通过中值滤波方法对参数数据m进行平滑处理获得平滑的曲线；对平滑的曲线进行多高斯拟合后，根据最小二乘法获得n组高斯函数特征参数的特征值；通过多条件判断，对n组高斯函数特征参数的特征值进行筛选，筛选出n组高斯函数特征值；根据n组高斯函数特征值结合如下表达式计算凝血酶原时间：；其中：pt代表凝血酶原时间，为经过多条件判断筛选后的第1组至第n组高斯函数特征值。
6.优选的，所述中值滤波方法具体包括如下步骤：对参数数据m进行处理，具体是：步骤
①
、对参数数据m中的第i个点进行处理，获取第i个点前后各v个点，共得到2v
1个点；其中：i为大于等于2且小于等于m的自然数，v为大于等于1的自然数；步骤
②
、对2v 1个点的数值大小进行排序，形成单调上升或者单调下降序列，取序号为中心点的数值作为滤波输出；步骤
③
、取i=i 1，若i小于m，则返回步骤
①
；否则结束。
7.优选的，中值滤波方法的数学表达式为：；其中：zi为滤波输出的数值，zi为第i个点的数值，n
*
代表正整数集。
8.优选的，多高斯拟合采用的函数表达式为：；其中：aj、bj以及cj均为多高斯拟合函数的特征参数，t为时间轴上的参数，j为正整数，e为自然对数。
9.优选的，多条件判断具体包括如下步骤：步骤
①
、对n组高斯函数特征参数的特征值中的第k组特征值（ak，bk，ck）进行判断，剔除假峰和无效峰的特征值，其中：k为小于等于n的自然数，判断条件具体包括：条件一：对所有特征值aj的组合进行冒泡排序，获得最大值a
max
；当所述特征值时，将此第k组特征值删除；条件二：对所有特征值bj的组合进行均值计算，得到平均值b
mid
；当时，将此第k组特征值删除；条件三：当特征值bk小于5s或者大于40s时，将此第k组特征值删除；步骤
②
、取k=k 1，若k小于n，则返回步骤
①
；否则进入下一步；步骤
③
、将筛选出的n组高斯函数特征值保留，形成新的数组。
10.应用本发明提供的检测方法，具有的有益效果为：本发明提供的凝血酶原检测方法中，通过中值滤波方法对参数数据进行了平滑处理，可以去除电信号中的异常点、噪声点以及干扰点，能够有效避免脉冲噪声、椒盐噪声等噪声干扰对凝血酶原时间检测精度的影响，进而提高检测结果的信噪比，从而提高检测的准确性和可靠性；除此之外，本发明通过多条件判断，对高斯函数特征参数的特征值进行了筛选，剔除了假峰和无效峰的特征值，再进行凝血酶原时间的计算，能够有效减少不对称性对高斯拟合结果的影响，同时又减轻了指数修正不对称高斯拟合函数的计算难度。
11.本发明还提供一种凝血酶原时间的检测装置，包括数据采集模块和数据处理模块：所述数据采集模块用于获取利用电化学方法采集的阻抗-时间曲线或者电流-时间曲线的原始数据，获得参数数据m；所述数据处理模块用于对数据采集模块获得的参数数据m进行中值滤波、多高斯拟合、特征值计算、多条件判断以及凝血酶原时间计算的数据分析、存储以及传输。
12.优选的，还包括中心控制模块和用户界面模块：所述中心控制模块用于控制数据采集模块、数据处理模块以及用户界面模块的操
作；所述用户界面模块包括按键和显示屏，所述按键用于设置仪器参数和回放测试结果；所述显示屏用于显示输入的检测信息、存储的病人信息和操作者信息以及检测卡编码的相关指令。
13.此外，为实现上述目的，本发明还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现上述凝血酶原时间的检测方法。
14.除了上面所描述的目的、特征和优点之外，本发明还有其它的目的、特征和优点。下面将参照图，对本发明作进一步详细的说明。
附图说明
15.构成本技术的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：图1是本发明优选实施例1提供的凝血酶原时间检测方法的流程示意图；图2是本发明优选实施例2提供的凝血酶原时间检测装置的模块示意图。
具体实施方式
16.以下结合附图对本发明的实施例进行详细说明，但是本发明可以根据权利要求限定和覆盖的多种不同方式实施。
17.实施例1：参见图1，本实施例公开一种优选的凝血酶原时间检测方法的流程示意图，包括：原始数据采集：将待测样本滴入凝血酶原时间电化学试纸的进样口，按照预设时间间隔对原始数据进行采集，得到参数数据m；所述参数数据m为阻抗-时间曲线或者电流-时间曲线上的数据；中值滤波：对参数数据m进行中值滤波平滑处理获得平滑的曲线；中值滤波的数学表达式为：；其中：zi为滤波输出的数值，zi为第i个点的数值，n
*
代表正整数集；处理的具体步如下：步骤
①
、对参数数据m中的第i个点进行处理，获取第i个点前后各v个点，共得到2v 1个点；其中：i为大于等于2且小于等于m的自然数，v为大于等于1的自然数；步骤
②
、对2v 1个点的数值大小进行排序，形成单调上升或者单调下降序列，取序号为中心点的数值作为滤波输出；步骤
③
、取i=i 1，若i小于m，则返回步骤
①
；否则结束。
18.此时，通过中值滤波方法对参数数据m进行处理，可以去除电信号中的异常点、噪声点以及干扰点，能够有效避免脉冲噪声和椒盐噪声对检测精度的干扰。
19.多高斯拟合及特征值计算：考虑到参数数据m所在的阻抗-时间曲线或者电流-时间曲线均为非对称性的多峰峰谱，进一步地，对中值滤波平滑后的曲线进行多高斯拟合，再根据最小二乘法计算n组高斯函数特征参数的特征值；处理的具体步骤如下：
多高斯拟合的函数表达式为：；其中：aj、bj以及cj均为多高斯拟合函数的特征参数，t为时间轴上的参数，j为正整数，e为自然对数。
20.优选的，多高斯拟合的函数表达式中，参数j为1-9范围内的正整数，当参数j的取值太小时，拟合曲线的准确度较低；当参数j的取值太大时，多高斯拟合处理时间较长，因此，参数j的优选范围为3-7，本实施例优选参数j=4。
21.此时，再将中值滤波的分析结果进行多高斯拟合，根据最小二乘法获得n组高斯函数特征参数的特征值（a1，b1，c1；a2，b2，c2；a3，b3，c3；.......；an，bn，cn）。
22.多条件判断：对n组高斯函数特征参数的特征值进行多条件判断，剔除假峰和无效峰的特征值，筛选出n组高斯函数特征值；判断的具体步骤为：步骤
①
、对n组高斯函数特征参数的特征值中的第k组特征值（ak，bk，ck）进行判断，剔除假峰和无效峰的特征值，其中：k为小于等于n的自然数，判断条件具体包括：条件一：对所有特征值aj的组合进行冒泡排序，获得最大值a
max
；当所述特征值时，将此第k组特征值删除；条件二：对所有特征值bj的组合进行均值计算，得到平均值b
mid
；当时，将此第k组特征值删除；条件三：当特征值bk小于5s或者大于40s时，将此第k组特征值删除；步骤
②
、取k=k 1，若k小于n，则返回步骤
①
；否则进入下一步；步骤
③
、将筛选出的n组高斯函数特征值保留，形成新的数组。
23.构建公式，计算凝血酶原时间：根据筛选出的n组高斯函数特征值结合如下表达式计算凝血酶原时间：；其中：pt代表凝血酶原时间，为经过多条件判断筛选后的第1组至第n组高斯函数特征值。
24.对比例1：在对比例1中，将中值滤波方法替换为均值滤波方法；所述均值滤波方法的具体步骤如下：步骤
①
、对参数数据m的第p个点进行处理，获取第p个点前后各q个点，共得到2q 1个点；其中：p为大于等于2且小于等于m的自然数，q为大于等于1的自然数；步骤
②
、计算所述2q 1个点的平均值，取平均值作为滤波输出；步骤
③
、取p=p 1，若p小于m，则返回步骤
①
；否则结束。
25.其他步骤、参数以及条件与实施例1中的凝血酶原时间的检测方法一致。
26.对比例2：在对比例2中，省略实施例1中多高斯拟合及之后的步骤，直接在经过中值滤波平
滑处理后的曲线上寻找最大值所对应的时间，即为凝血酶原时间值。
27.其他步骤、参数以及条件与实施例1中的凝血酶原时间的检测方法一致。
28.对比例3：选用市售的便携式血凝仪及其检测卡。
29.对实施例1、对比例1、对比例2以及对比例3进行检验，具体的检验过程如下：将实施例1、对比例1、对比例2以及对比例3的检测卡与配套电化学检验仪配合使用。
30.四种检测卡分别对同样的血浆样品进行凝血酶原时间测试，其重复性结果见表1。
31.表1 凝血酶原时间测定结果表1表明，本实施例1、对比例1、对比例2以及对比例3相比，实施例1的变异系数值＜对比例3的变异系数值≦对比例2的变异系数值＜对比例1的变异系数。其中：实施例1的变异系数最小，且测试得到的凝血酶原时间重复性好。因此，通过本发明提供的检测方法来测定凝血酶原时间，能够有效避免噪声干扰对凝血酶原时间检测精度的影响，从而提高检测的准确性和可靠性。
32.实施例2：参见图2，本实施例公开一种优选的凝血酶原时间的检测装置，包括中心控制模块21、数据采集模块22、数据处理模块23以及用户界面模块24。
33.所述中心控制模块21运用mcu芯片，负责对数据采集模块、数据处理模块以及用户界面模块进行总控制。
34.所述数据采集模块22包含电化学试纸以及进样口，将待测样本滴入进样口后，能够通过电化学试纸检测酶反应过程中产生的电信号，对电信号等原始数据进行采集，获得参数数据，再通过中心控制模块将参数数据传输到数据处理模块中。
35.所述数据处理模块23主要用于实现对数据采集模块获得的参数数据的处理、存储以及传输功能，具体包括数据分析模块、数据存储模块以及数据传输模块。
36.所述数据分析模块用于对参数数据进行中值滤波处理，得到平滑的曲线；对平滑的曲线进行多高斯拟合，获取多个高斯函数特征参数的特征值；对高斯函数特征参数的特征值进行多条件判断筛选，并根据筛选出的特征值，构造公式计算凝血酶原时间。（参照实施例1）所述数据存储模块用于存储参数数据、检测结果以及病人信息和操作者信息等，本实施例不加以具体限定。
37.所述数据传输模块用于接收数据采集模块22传输的参数数据，以及向用户界面模块24输送检测结果、病人信息以及操作者信息等，本实施例不加以具体限定。
38.所述用户界面模块24包括按键和显示屏，所述按键用于设置仪器参数和回放测试结果；所述显示屏用于显示检测信息、病人信息和操作者信息等，本实施例不加以具体限定。
39.以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。