CN112710636B - 一种特定蛋白浓度的检测方法及检测装置 - Google Patents

一种特定蛋白浓度的检测方法及检测装置 Download PDF

Info

Publication number
CN112710636B
CN112710636B CN202011450007.5A CN202011450007A CN112710636B CN 112710636 B CN112710636 B CN 112710636B CN 202011450007 A CN202011450007 A CN 202011450007A CN 112710636 B CN112710636 B CN 112710636B
Authority
CN
China
Prior art keywords
specific protein
sample
curve
detected
time
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
CN202011450007.5A
Other languages
English (en)
Other versions
CN112710636A (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.)
Shenzhen Comen Medical Instruments Co Ltd
Original Assignee
Shenzhen Comen Medical Instruments Co Ltd
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 Shenzhen Comen Medical Instruments Co Ltd filed Critical Shenzhen Comen Medical Instruments Co Ltd
Priority to CN202011450007.5A priority Critical patent/CN112710636B/zh
Publication of CN112710636A publication Critical patent/CN112710636A/zh
Application granted granted Critical
Publication of CN112710636B publication Critical patent/CN112710636B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/59—Transmissivity
    • 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/53—Immunoassay; Biospecific binding assay; Materials therefor
    • G01N33/536—Immunoassay; Biospecific binding assay; Materials therefor with immune complex formed in liquid phase
    • 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/68—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving proteins, peptides or amino acids
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2333/00—Assays involving biological materials from specific organisms or of a specific nature
    • G01N2333/435—Assays involving biological materials from specific organisms or of a specific nature from animals; from humans
    • G01N2333/46—Assays involving biological materials from specific organisms or of a specific nature from animals; from humans from vertebrates
    • G01N2333/47—Assays involving proteins of known structure or function as defined in the subgroups
    • G01N2333/4701—Details
    • G01N2333/4737—C-reactive protein

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Immunology (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Molecular Biology (AREA)
  • Biomedical Technology (AREA)
  • Urology & Nephrology (AREA)
  • Hematology (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Microbiology (AREA)
  • Cell Biology (AREA)
  • Biotechnology (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

本发明提供了一种特定蛋白浓度的检测方法及检测装置，检测方法包括在预设的采样总时间内基于透射比浊法对待测样本进行特定蛋白检测，获取待测样本的特定蛋白反应曲线的一阶导数的最大值，根据一阶导数的最大值与预设阈值的关系，确定待测样本的特定蛋白反应曲线是否存在波动，提取待测样本的曲线特征，提取各已知浓度样本的特定蛋白反应曲线在提取时间段对应的曲线特征，建立定标函数关系式，将待测样本的曲线特征输入定标函数关系式中，得到待测样本的特定蛋白浓度。将上述检测方法应用于检测装置中能够在特定蛋白反应曲线发生抖动或鼓包等异常时也能计算出正确的结果，从而提高了检测的准确性和可靠性。

Description

一种特定蛋白浓度的检测方法及检测装置
技术领域
本发明涉及医疗技术领域，具体涉及一种特定蛋白检测方法，以及一种特定蛋白检测装置。
背景技术
C反应蛋白(CRP)是肝脏合成的一种急性炎症正时相反应蛋白，也是一种特定蛋白。正常人血中的CRP浓度很低，在机体突遇紧张、组织创伤和各种炎症刺激时合成快速增加，并从肝细胞中分泌入血液，在感染发生后12-18小时即可检测到高水平的CRP。在感染发生后的12-14天升高的CRP可降至基线水平。因此多年来—直是评价炎症性疾病的指标之一，并且升高幅度与感染的程度呈相关。CRP作为诊断细菌感染的重要标志物之一，已广泛应用于临床。CRP还是一个评估心脏病发生率、复发率、死亡率的临床重要指标。近年来的研究发现，炎症在动脉粥样硬化及肿瘤的发生、发展过程中起重要作用。鉴于血清CRP的重要作用，其测量的准确性受到了广泛的关注。
检测CRP的常用方法多种多样，其中包括散射比浊法、透射比浊法，放射免疫测定、化学发光法、ELISA法及床旁CRP检测(POCT)等。目前临床实验室测定血清中CRP的方法主要是免疫浊度法，包括乳胶增强透射比浊法和速率散射比浊法，这两种方法主要用于自动化分析系统，速率散射比浊法多用于免疫检测领域的封闭检测系统，乳胶增强透射比浊法多用于生化检测领域的开放检测系统。
目前基于透射比浊法对样本的特定蛋白浓度进行检测时，首先得到该样本的特定蛋白反应曲线，然后求得特定蛋白反应曲线的一个曲线特征，接着将该曲线特征代入定标曲线，则可求出样本的特定蛋白浓度。上述方式如果在样本采集(也就是检测样本的特定蛋白浓度)时，如果某一段时间信号存在干扰，特定蛋白反应曲线会出现抖动或鼓包等异常情况，如果选取计算的曲线特征的点正好位于特定蛋白反应曲线抖动或鼓包的位置，就会导致计算结果不准，造成假阴性或假阳性的结果。
发明内容
根据第一方面，一种实施例中提供一种特定蛋白浓度的检测方法,包括：
在预设的采样总时间内基于透射比浊法对待测样本进行特定蛋白检测，以获取待测样本的特定蛋白反应曲线，所述待测样本的特定蛋白反应曲线用于表示采样总时间内由透射比浊法得到的电压的变化；
获取所述待测样本的特定蛋白反应曲线的一阶导数的最大值；
根据所述一阶导数的最大值与预设阈值的关系，确定所述待测样本的特定蛋白反应曲线是否存在波动，如果不存在波动，在所述采样总时间内选取提取时间段，如果存在波动，则获取所述待测样本的特定蛋白反应曲线的二阶导数，并根据所述二阶导数确定采样总时间内中的采样有效时间，在所述采样有效时间内选取提取时间段；
提取所述待测样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征；
根据预先基于透射比浊法得到的多个已知浓度样本的特定蛋白反应曲线，提取各已知浓度样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征；
建立各已知浓度样本的特定蛋白反应曲线的曲线特征与其对应的特定蛋白浓度间的定标函数关系式；
将所述待测样本的曲线特征输入所述定标函数关系式中，得到所述待测样本的特定蛋白浓度。
根据第二方面，一种实施例中提供一种计算机可读存储介质，包括程序，所述程序能够被处理器执行以实现第一方面中的方法。
根据第三方面，一种实施例中提供一种特定蛋白浓度的检测装置，包括：
反应容器，用于容纳待测样本；
光源，用于向所述反应容器内的待测样本提供激光；
光信号接收器，用于采集激光经待测样本形成的透射光，并将透射光的光信号转换为电信号；
存储装置，用于存储各已知浓度样本的特定蛋白浓度和预先基于透射比浊法得到的对应的特定蛋白反应曲线；
数据处理装置，与所述光信号接收器和存储装置分别信号连接，用于：
根据所述电信号获取待测样本的特定蛋白反应曲线，所述特定蛋白反应曲线用于表示采样总时间内得到的电压的变化；
获取所述待测样本的特定蛋白反应曲线的一阶导数的最大值；
根据所述一阶导数的最大值与预设阈值的关系，确定所述待测样本的特定蛋白反应曲线是否存在波动，如果不存在波动，在所述采样总时间内选取提取时间段，如果存在波动，则获取所述待测样本的特定蛋白反应曲线的二阶导数，并根据所述二阶导数确定采样总时间内中的采样有效时间，在所述采样有效时间内选取提取时间段；
提取所述待测样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征；
根据预先基于透射比浊法得到的多个已知浓度样本的特定蛋白反应曲线，提取各已知浓度样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征；
建立各已知浓度样本的特定蛋白反应曲线的曲线特征与其对应的特定蛋白浓度间的定标函数关系式；
将所述待测样本的曲线特征输入所述定标函数关系式中，得到所述待测样本的特定蛋白浓度。
上述实施例中，首先判断特定蛋白反应曲线是否存在波动(例如抖动或鼓包的)，如果不存在波动，直接获取曲线特征，如果存在波动，又识别出特定蛋白反应曲线不存在波动的部分，并由此得到待测样本的曲线特征。该方式降低了样本采集时对特定蛋白反应曲线完整性的依赖，在特定蛋白反应曲线发生抖动或鼓包等异常时也能计算出正确的结果，从而提高了检测的准确性和可靠性。
附图说明
图1为一种实施例的特定蛋白浓度检测装置的结构示意图；
图2为一种实施例的特定蛋白反应曲线的示意图；
图3为另一种实施例的特定蛋白反应曲线的示意图；
图4为又一种实施例的特定蛋白反应曲线的示意图；
图5为一种实施例的特定蛋白浓度检测方法的流程图。
10、反应容器；
20、光源；
30、光信号接收器；
40、存储装置；
50、数据处理装置。
具体实施方式
下面通过具体实施方式结合附图对本发明作进一步详细说明。其中不同实施方式中类似元件采用了相关联的类似的元件标号。在以下的实施方式中，很多细节描述是为了使得本申请能被更好的理解。然而，本领域技术人员可以毫不费力的认识到，其中部分特征在不同情况下是可以省略的，或者可以由其他元件、材料、方法所替代。在某些情况下，本申请相关的一些操作并没有在说明书中显示或者描述，这是为了避免本申请的核心部分被过多的描述所淹没，而对于本领域技术人员而言，详细描述这些相关操作并不是必要的，他们根据说明书中的描述以及本领域的一般技术知识即可完整了解相关操作。
另外，说明书中所描述的特点、操作或者特征可以以任意适当的方式结合形成各种实施方式。同时，方法描述中的各步骤或者动作也可以按照本领域技术人员所能显而易见的方式进行顺序调换或调整。因此，说明书和附图中的各种顺序只是为了清楚描述某一个实施例，并不意味着是必须的顺序，除非另有说明其中某个顺序是必须遵循的。
本文中为部件所编序号本身，例如“第一”、“第二”等，仅用于区分所描述的对象，不具有任何顺序或技术含义。而本申请所说“连接”、“联接”，如无特别说明，均包括直接和间接连接(联接)。
请参考图1，为透射比浊法的原理示意图，包括光源20、反应容器10和光信号接收器30，其中，反应容器10中装有反应液，光源20垂直射入反应容器10，穿过反应容器10壁照射到反应液中的微球上，发生透射，再通过反应容器10壁进入空气，最后到光信号接收器30。
特定蛋白反应检测原理是依据抗原、抗体在特定的电解质溶液中反应，快速形成免疫复合物微粒，使反应液出现浊度，随着时间的推移，复合物微粒聚合得就会越大，浊度也逐渐变大。在透射比浊法中，光信号接收器30接收到光信号的强度，随时间的推移越来越小。C反应蛋白是特定蛋白检测项目的一种。
请参考图2，为一种特定蛋白反应曲线示意图，坐标横轴是采样时间，坐标纵轴是电压，整个曲线是光信号接收器30接收的光信号强度随时间的变化曲线，通常用数学函数公式表示特定蛋白反应曲线，其包括：
V＝F(t)。
其中，0≤t，V是电压，t是采样时间，整个采样过程的持续时间在本发明中定义为采样总时间,特定蛋白反应曲线为一减函数。
请参照图1，本发明提供了一种特定蛋白浓度的检测装置，包括反应容器10、光源20、光信号接收器30、存储装置40和数据处理装置50。
反应容器10用于容纳待测样本，例如可以是透明的反应容器10。
光源20作为透射比浊法使用的光源20，其用于向反应容器10内的待测样本提供激光。
光信号接收器30设置在激光经待测样本形成的透射光的光路上，其用于采集透射光，并将透射光的光信号转换为电信号。
存储装置40用于存储各已知浓度样本的特定蛋白浓度和预先基于透射比浊法得到的对应的特定蛋白反应曲线。
一些实施例中，已知浓度样本的数目为N个，即具有N个按浓度梯度排列的已知浓度样本B1,B2,…，Bi,…,BN；其中，第i个已知浓度样本的特定蛋白浓度为Ci，且0<C1<C2<…<Ci<…<CN，1≤i≤N。
第i个已知浓度样本的定蛋白反应曲线的公式为：
Vbi＝Fi(tb)。
其中，0≤tb≤Tb，tb∈实数，tb为已知浓度样本的采样时间，Vbi为第i个已知浓度样本采样时获取的电压，Tb为已知浓度样本的采样总时间的结束时间点。
一些实施例中，对于任一已知浓度样本的特定蛋白反应曲线，还判断已知浓度样本的特定蛋白反应曲线是否存在波动，如果存在波动，则对存在波动的特定蛋白反应曲线的已知浓度样本重新检测，以再次获取相应的特定蛋白反应曲线，直到每个已知浓度样本的特定蛋白反应曲线均不存在波动。对于第i个已知浓度样本Bi，判断其特定蛋白反应曲线是否存在波动的方式，可以是：获取第i个已知浓度样本的特定蛋白反应曲线的一阶导数Fi’(tb)，然后获取一阶导数Fi’(tb)的最大值max(Fi'(tb))，如果max(Fi'(tb))不大于零，则第i个已知浓度样本的特定蛋白反应曲线不存在波动，否则，第i个已知浓度样本的特定蛋白反应曲线存在波动,那么就再次对该已知浓度样本进行检测，直到该已知浓度样本的特定蛋白反应曲线的一阶导数的最大值不大于0。
数据处理装置50与光信号接收器30和存储装置40分别信号连接，其用于根据电信号获取待测样本的特定蛋白反应曲线，特定蛋白反应曲线用于表示采样总时间内得到的电压的变化。
获取到的待测样本的特定蛋白反应曲线的公式,包括：
Va＝F(ta)。
其中，0≤ta≤Ta，ta∈实数，ta为待测样本的采样时间，Va为采样时获取的电压，Ta为采样总时间的结束时间点。
而后获取待测样本的特定蛋白反应曲线的一阶导数F’(ta)的最大值max(F'(ta))，根据一阶导数的最大值与预设阈值的关系，确定待测样本的特定蛋白反应曲线是否存在波动。一些实施例中，该预设阈值为零，如果max(F'(ta))小于等于零，代表待测样本的特定蛋白反应曲线为一减函数，待测样本的特定蛋白反应曲线不存在波动，如果max(F'(ta))大于零，待测样本的特定蛋白反应曲线存在波动。
如果待测样本的特定蛋白反应曲线不存在波动，在采样总时间内选取提取时间段，如果存在波动，则获取待测样本的特定蛋白反应曲线的二阶导数，并根据二阶导数确定采样总时间内中的采样有效时间，在采样有效时间内选取提取时间段。采样有效时间就是特定蛋白反应曲线不存在波动的时间段，具体的，数据处理装置50获取待测样本的特定蛋白反应曲线的二阶导数F”(ta)，根据二阶导数确定采样有效时间。采样有效时间内待测样本的特定蛋白反应曲线需要满足F”(ta)大于零。并且，0≤t1≤t2≤Ta，t1为采样有效时间的起始时间点，t2为采样有效时间的结束时间点。如图3(箭头指向特定蛋白反应曲线处为波动处)所示，通常，满足上述条件的时间段只有一段。如果该时间段有两段或两段以上，选取时间长度最大的时间段作为采样有效时间，例如如图4所示，一满足上述条件的时间段的长度(t'2-t'1)较小，就不作为采样有效时间。
获取到提取时间段后，数据处理装置50提取待测样本的特定蛋白反应曲线在提取时间段对应的曲线特征。
一些实施例中，如图2所示，待测样本的特定蛋白反应曲线不存在波动，提取待测样本的特定蛋白反应曲线在提取时间段对应的曲线特征，包括：
计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的电压差的绝对值|Da1|，电压差的绝对值|Da1|的计算公式包括：
|Da1|＝|F(ta2)-F(ta1)|。
其中，ta1为提取时间段的起始时间点，ta2为提取时间段的结束时间点，0≤ta1＜ta2≤Ta，Ta为采样总时间的结束时间点，将电压差的绝对值|Da1|作为待测样本的特定蛋白反应曲线的曲线特征。
在另一些实施例中，如图2所示不存在波动，计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的面积Sa1，面积Sa1的计算公式包括：
Figure BDA0002826419530000061
其中，ta1为提取时间段的起始时间点，ta2为提取时间段的结束时间点，0≤ta1＜ta2≤Ta，Ta为采样总时间的结束时间点，将面积Sa1作为待测样本的特定蛋白反应曲线的曲线特征。
一些实施例中，如图3或图4(箭头指向特定蛋白反应曲线处为波动处)所示，待测样本的特定蛋白反应曲线存在波动，提取待测样本的特定蛋白反应曲线在提取时间段对应的曲线特征，包括：
计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的电压差的绝对值|Da2|，电压差的绝对值|Da2|的计算公式包括：
|Da2|＝|F(ta2)-F(ta1)|。
其中，ta1为提取时间段的起始时间点，ta2为提取时间段的结束时间点，t1≤ta1＜ta2≤t2，t1为采样有效时间的起始时间点，t2为采样有效时间的结束时间点，将电压差的绝对值|Da2|作为待测样本的特定蛋白反应曲线的曲线特征。
在另一些实施例中，如图3或图4(箭头指向特定蛋白反应曲线处为波动处)所示存在波动，计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的面积Sa2，面积Sa2的计算公式包括：
Figure BDA0002826419530000071
其中，ta1为提取时间段的起始时间点，ta2为提取时间段的结束时间点，t1≤ta1＜ta2≤t2，t1为采样有效时间的起始时间点，t2为采样有效时间的结束时间点，将面积Sa2作为待测样本的特定蛋白反应曲线的曲线特征。
而后，根据存储装置40中存储的N个已知浓度样本的特定蛋白反应曲线，获取各已知浓度样本的特定蛋白反应曲线在提取时间段对应的曲线特征r1,r2,…，ri,…,rN,1≤i≤N，ri为第i个已知浓度样本的特定蛋白反应曲线在提取时间段对应的曲线特征。
一些实施例中，对于第i个已知浓度样本Bi，计算已知浓度样本的特定蛋白反应曲线上不同两点间的电压差的绝对值，以获取电压差的绝对值|Dbi|，电压差的绝对值|Dbi|的计算公式包括：
|Dbi|＝|Fi(tb2)-Fi(tb1)|。
其中，tb1和tb2是第i个已知浓度样本的特定蛋白反应曲线上两点的采样时间，0≤tb1＜tb2≤Tb，|Dbi|为电压差的绝对值，Tb为已知浓度样本的采样总时间的结束时间点，如果待测样本的曲线特征为|Da1|或|Da2|，则取tb2＝ta2，tb1＝ta1，将电压差的绝对值|Dbi|作为第i个已知浓度样本的特定蛋白反应曲线的曲线特征ri。
一些实施例中，对于第i个已知浓度样本Bi，计算已知浓度样本的特定蛋白反应曲线上不同两点间的面积，以获取面积Sbi，面积Sbi的计算公式包括：
Figure BDA0002826419530000072
其中，Sbi为已知浓度样本的特定蛋白反应曲线上不同两点间的面积，tb1和tb2是已知浓度样本的特定蛋白反应曲线上两点的采样时间，0≤tb1＜tb2≤Tb，Tb为已知浓度样本采样总时间的结束时间点，如果待测样本的曲线特征为Sa1或Sa2，则取tb2＝ta2，tb1＝ta1，将面积Sbi作为第i个已知浓度样本的特定蛋白反应曲线的曲线特征ri。
提取到每个已知浓度样本的特定蛋白反应曲线的特征值后，将每个已知浓度样本的特定蛋白反应曲线的特征值与其对应的特定蛋白浓度建立数学函数对应关系，其中，将第i个已知浓度样本的特定蛋白反应曲线的特征值ri与其对应的特定蛋白浓度Ci建立数学函数对应关系，1≤i≤N，进而获取定标函数关系式：r＝F(C)，其中，C为已知浓度样本的特定蛋白浓度，r为已知浓度样本的特定蛋白反应曲线的特征值。
得到待测样本的特定蛋白反应曲线的曲线特征后，数据处理装置50将待测样本的曲线特征输入定标函数关系式中，从而得到一个特定蛋白浓度，将该特定蛋白浓度作为待测样本的特定蛋白浓度。
本发明还提供了一种特定蛋白浓度的检测方法，如图5所示，包括步骤：
步骤100、在预设的采样总时间内基于透射比浊法对待测样本进行特定蛋白检测，以获取待测样本的特定蛋白反应曲线，待测样本的特定蛋白反应曲线用于表示采样总时间内由透射比浊法得到的电压的变化。
获取到的待测样本的特定蛋白反应曲线的公式,包括：
Va＝F(ta)。
其中，0≤ta≤Ta，ta∈实数，ta为待测样本的采样时间，Va为采样时获取的电压，Ta为采样总时间的结束时间点。
步骤200、获取待测样本的特定蛋白反应曲线的一阶导数的最大值。即获取待测样本的特定蛋白反应曲线的一阶导数F’(ta)的最大值max(F'(ta))。
步骤300、根据一阶导数的最大值与预设阈值的关系，确定待测样本的特定蛋白反应曲线是否存在波动，如果不存在波动，执行步骤400，如果存在波动，执行步骤500。
一些实施例中，该预设阈值为零，如果max(F'(ta))小于等于零，代表待测样本的特定蛋白反应曲线为一减函数，待测样本的特定蛋白反应曲线不存在波动，如果max(F'(ta))大于零，待测样本的特定蛋白反应曲线存在波动。
步骤400、在采样总时间内选取提取时间段后，执行步骤700。
步骤500、获取待测样本的特定蛋白反应曲线的二阶导数，并根据二阶导数确定采样总时间内中的采样有效时间后，执行步骤600。
采样有效时间内待测样本的特定蛋白反应曲线需要满足F”(ta)大于零。并且，0≤t1≤t2≤Ta，t1为采样有效时间的起始时间点，t2为采样有效时间的结束时间点。通常，满足上述条件的时间段只有一段，如果该时间段有两段或两段以上，则选取时间长度最大的时间段作为采样有效时间，例如如图4所示，一满足上述条件的时间段的长度(t'2-t'1)较小，就不作为采样有效时间。
步骤600、在采样有效时间内选取提取时间段后，执行步骤700。
步骤700、提取待测样本的特定蛋白反应曲线在提取时间段对应的曲线特征。
一些实施例中，如图2所示，待测样本的特定蛋白反应曲线不存在波动，提取待测样本的特定蛋白反应曲线在提取时间段对应的曲线特征，包括：
计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的电压差的绝对值|Da1|，电压差的绝对值|Da1|的计算公式包括：
|Da1|＝|F(ta2)-F(ta1)|。
其中，ta1为提取时间段的起始时间点，ta2为提取时间段的结束时间点，0≤ta1＜ta2≤Ta，Ta为采样总时间的结束时间点，将电压差的绝对值|Da1|作为待测样本的特定蛋白反应曲线的曲线特征。
在另一些实施例中，如图2所示不存在波动，计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的面积Sa1，面积Sa1的计算公式包括：
Figure BDA0002826419530000091
其中，ta1为提取时间段的起始时间点，ta2为提取时间段的结束时间点，0≤ta1＜ta2≤Ta，Ta为采样总时间的结束时间点，将面积Sa1作为待测样本的特定蛋白反应曲线的曲线特征。
一些实施例中，如图3或图4(箭头指向特定蛋白反应曲线处为波动处)所示，待测样本的特定蛋白反应曲线存在波动，提取待测样本的特定蛋白反应曲线在提取时间段对应的曲线特征，包括：
计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的电压差的绝对值|Da2|，电压差的绝对值|Da2|的计算公式包括：
|Da2|＝|F(ta2)-F(ta1)|。
其中，ta1为提取时间段的起始时间点，ta2为提取时间段的结束时间点，t1≤ta1＜ta2≤t2，t1为采样有效时间的起始时间点，t2为采样有效时间的结束时间点，将电压差的绝对值|Da2|作为待测样本的特定蛋白反应曲线的曲线特征。
在另一些实施例中，如图3或图4(箭头指向特定蛋白反应曲线处为波动处)所示存在波动，计算待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的面积Sa2，面积Sa2的计算公式包括：
Figure BDA0002826419530000092
其中，ta1为提取时间段的起始时间点，ta2为提取时间段的结束时间点，t1≤ta1＜ta2≤t2，t1为采样有效时间的起始时间点，t2为采样有效时间的结束时间点，将面积Sa2作为待测样本的特定蛋白反应曲线的曲线特征。
步骤800、根据预先基于透射比浊法得到的多个已知浓度样本的特定蛋白反应曲线，提取各已知浓度样本的特定蛋白反应曲线在提取时间段对应的曲线特征。
一些实施例中，已知浓度样本的数目为N个，即具有N个按浓度梯度排列的已知浓度样本B1,B2,…，Bi,…,BN；其中，第i个已知浓度样本的特定蛋白浓度为Ci，且0<C1<C2<…<Ci<…<CN，1≤i≤N。
第i个已知浓度样本的定蛋白反应曲线的公式为：
Vbi＝Fi(tb)。
其中，0≤tb≤Tb，tb∈实数，tb为已知浓度样本的采样时间，Vbi为第i个已知浓度样本采样时获取的电压，Tb为已知浓度样本的采样总时间的结束时间点。
一些实施例中，对于任一已知浓度样本的特定蛋白反应曲线，还判断已知浓度样本的特定蛋白反应曲线是否存在波动，如果存在波动，则对存在波动的特定蛋白反应曲线的已知浓度样本重新检测，以再次获取相应的特定蛋白反应曲线，直到每个已知浓度样本的特定蛋白反应曲线均不存在波动。对于第i个已知浓度样本Bi，判断其特定蛋白反应曲线是否存在波动的方式，可以是：获取第i个已知浓度样本的特定蛋白反应曲线的一阶导数Fi’(tb)，然后获取一阶导数Fi’(tb)的最大值max(Fi'(tb))，如果max(Fi'(tb))不大于零，则第i个已知浓度样本的特定蛋白反应曲线不存在波动，否则，第i个已知浓度样本的特定蛋白反应曲线存在波动,那么就再次对该已知浓度样本进行检测，直到该已知浓度样本的特定蛋白反应曲线的一阶导数的最大值不大于0。
基于上述预先获取的各已知浓度样本的特定蛋白曲线，步骤800中获取各已知浓度样本的特定蛋白反应曲线在提取时间段对应的曲线特征r1,r2,…，ri,…,rN,1≤i≤N，ri为第i个已知浓度样本的特定蛋白反应曲线在提取时间段对应的曲线特征。
一些实施例中，对于第i个已知浓度样本Bi，计算已知浓度样本的特定蛋白反应曲线上不同两点间的电压差的绝对值，以获取电压差的绝对值|Dbi|，电压差的绝对值|Dbi|的计算公式包括：
|Dbi|＝|Fi(tb2)-Fi(tb1)|。
其中，tb1和tb2是第i个已知浓度样本的特定蛋白反应曲线上两点的采样时间，0≤tb1＜tb2≤Tb，|Dbi|为电压差的绝对值，Tb为已知浓度样本的采样总时间的结束时间点，如果待测样本的曲线特征为Da1或Da2，则取tb2＝ta2，tb1＝ta1，将电压差的绝对值|Dbi|作为第i个已知浓度样本的特定蛋白反应曲线的曲线特征ri。
一些实施例中，对于第i个已知浓度样本Bi，计算已知浓度样本的特定蛋白反应曲线上不同两点间的面积，以获取面积Sbi，面积Sbi的计算公式包括：
Figure BDA0002826419530000111
其中，Sbi为已知浓度样本的特定蛋白反应曲线上不同两点间的面积，tb1和tb2是已知浓度样本的特定蛋白反应曲线上两点的采样时间，0≤tb1＜tb2≤Tb，Tb为已知浓度样本采样总时间的结束时间点，如果待测样本的曲线特征为Sa1或Sa2，则取tb2＝ta2，tb1＝ta1，将面积Sbi作为第i个已知浓度样本的特定蛋白反应曲线的曲线特征ri。
步骤900、建立各已知浓度样本的特定蛋白反应曲线的曲线特征与其对应的特定蛋白浓度间的定标函数关系式。
将第i个已知浓度样本的特定蛋白反应曲线的特征值ri与其对应的特定蛋白浓度Ci建立数学函数对应关系，1≤i≤N，进而获取定标函数关系式：r＝F(C)，其中，C为已知浓度样本的特定蛋白浓度，r为已知浓度样本的特定蛋白反应曲线的特征值。
步骤1000、将待测样本的曲线特征输入定标函数关系式中，得到待测样本的特定蛋白浓度。
上述实施例中，首先判断特定蛋白反应曲线是否存在波动(例如抖动或鼓包的)，如果不存在波动，直接获取曲线特征，如果存在波动，又识别出特定蛋白反应曲线不存在波动的部分，并由此得到待测样本的曲线特征。该方式降低了样本采集时对特定蛋白反应曲线完整性的依赖，在特定蛋白反应曲线发生抖动或鼓包等异常时也能计算出正确的结果，从而提高了检测的准确性和可靠性。
本领域技术人员可以理解，上述实施方式中各种方法的全部或部分功能可以通过硬件的方式实现，也可以通过计算机程序的方式实现。当上述实施方式中全部或部分功能通过计算机程序的方式实现时，该程序可以存储于一计算机可读存储介质中，存储介质可以包括：只读存储器、随机存储器、磁盘、光盘、硬盘等，通过计算机执行该程序以实现上述功能。例如，将程序存储在设备的存储器中，当通过处理器执行存储器中程序，即可实现上述全部或部分功能。另外，当上述实施方式中全部或部分功能通过计算机程序的方式实现时，该程序也可以存储在服务器、另一计算机、磁盘、光盘、闪存盘或移动硬盘等存储介质中，通过下载或复制保存到本地设备的存储器中，或对本地设备的系统进行版本更新，当通过处理器执行存储器中的程序时，即可实现上述实施方式中全部或部分功能。
以上应用了具体个例对本发明进行阐述，只是用于帮助理解本发明，并不用以限制本发明。对于本发明所属技术领域的技术人员，依据本发明的思想，还可以做出若干简单推演、变形或替换。

Claims (9)

1.一种特定蛋白浓度的检测方法,所述检测方法用于非诊断和/或治疗目的，其特征在于,包括：
在预设的采样总时间内基于透射比浊法对待测样本进行特定蛋白检测，以获取待测样本的特定蛋白反应曲线，所述待测样本的特定蛋白反应曲线用于表示采样总时间内由透射比浊法得到的电压的变化；
获取所述待测样本的特定蛋白反应曲线的一阶导数的最大值；
根据所述一阶导数的最大值与预设阈值的关系，确定所述待测样本的特定蛋白反应曲线是否存在波动，如果不存在波动，在所述采样总时间内选取提取时间段，如果存在波动，则获取所述待测样本的特定蛋白反应曲线的二阶导数，并根据所述二阶导数确定采样总时间内中的采样有效时间，在所述采样有效时间内选取提取时间段；所述根据所述二阶导数确定采样总时间内中的采样有效时间，包括：获取所述采样总时间内二阶导数大于零的至少一个时间段，将所述至少一个时间段中长度最大的时间段作为采样有效时间；
提取所述待测样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征；
根据预先基于透射比浊法得到的多个已知浓度样本的特定蛋白反应曲线，提取各已知浓度样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征；
建立各已知浓度样本的特定蛋白反应曲线的曲线特征与其对应的特定蛋白浓度间的定标函数关系式；
将所述待测样本的曲线特征输入所述定标函数关系式中，得到所述待测样本的特定蛋白浓度。
2.如权利要求1所述的方法，其特征在于，所述预设阈值的大小为零，所述根据所述一阶导数的最大值与预设阈值的关系，确定所述待测样本的特定蛋白反应曲线是否存在波动，包括：
判断所述一阶导数的最大值是否不大于零，如果不大于零，则所述待测样本的特定蛋白反应曲线不存在波动，如果大于零，则所述待测样本的特定蛋白反应曲线存在波动。
3.如权利要求1所述的方法，其特征在于，获取到的待测样本的特定蛋白反应曲线的公式为：
Va=F(ta)，
其中，0≤ta≤Ta，ta∈实数，ta为待测样本的采样时间，Va为采样时获取的电压，Ta为采样总时间的结束时间点。
4.如权利要求3所述的方法,其特征在于,如果不存在波动，提取所述待测样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征，包括：
计算所述待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的电压差的绝对值|Da1|，所述电压差的绝对值|Da1|的计算公式包括：
|Da1|=|F(ta2)-F（ta1）|，
其中，ta1为所述提取时间段的起始时间点，ta2为所述提取时间段的结束时间点， 0≤ta1＜ta2≤Ta ，Ta为采样总时间的结束时间点，将所述电压差的绝对值|Da1|作为所述待测样本的特定蛋白反应曲线的曲线特征；或者
计算所述待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的面积Sa1，所述面积Sa1的计算公式包括：
Figure 857235DEST_PATH_IMAGE001
，
其中，ta1为所述提取时间段的起始时间点，ta2为所述提取时间段的结束时间点， 0≤ta1＜ta2≤Ta ，Ta为采样总时间的结束时间点，将所述面积Sa1作为所述待测样本的特定蛋白反应曲线的曲线特征。
5.如权利要求3所述的方法，其特征在于，如果存在波动，提取所述待测样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征，包括：
计算所述待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的电压差的绝对值|Da2|，所述电压差的绝对值|Da2|的计算公式包括：
|Da2|=|F(ta2)-F（ta1）|，
其中，ta1为所述提取时间段的起始时间点，ta2为所述提取时间段的结束时间点， t1≤ta1＜ta2≤t2， t1为采样有效时间的起始时间点，t2为采样有效时间的结束时间点，将所述电压差的绝对值|Da2|作为所述待测样本的特定蛋白反应曲线的曲线特征；或者
计算所述待测样本的特定蛋白反应曲线上提取时间段的起始时间点和结束时间点间的面积Sa2，所述面积Sa2的计算公式包括：
Figure DEST_PATH_IMAGE002
，
其中，ta1为所述提取时间段的起始时间点，ta2为所述提取时间段的结束时间点， t1≤ta1＜ta2≤t2， t1为采样有效时间的起始时间点，t2为采样有效时间的结束时间点，将所述面积Sa2作为所述待测样本的特定蛋白反应曲线的曲线特征。
6.如权利要求1所述的方法，其特征在于，基于透射比浊法得到多个已知浓度样本的特定蛋白反应曲线的方式，包括：
对于任一已知浓度样本，基于透射比浊法获取所述已知浓度样本的特定蛋白反应曲线，判断所述已知浓度样本的特定蛋白反应曲线是否存在波动，如果存在波动，则对存在波动的特定蛋白反应曲线的已知浓度样本重新检测，以再次获取相应的特定蛋白反应曲线，直到每个已知浓度样本的特定蛋白反应曲线均不存在波动。
7.如权利要求6所述的方法，其特征在于，所述判断所述已知浓度样本的特定蛋白反应曲线是否存在波动，包括：
获取所述已知浓度样本的特定蛋白反应曲线的一阶导数的最大值；
如果所述已知浓度样本的特定蛋白反应曲线的一阶导数的最大值不大于零，则所述已知浓度样本的特定蛋白反应曲线不存在波动，否则，所述已知浓度样本的特定蛋白反应曲线存在波动。
8.一种计算机可读存储介质，其特征在于，包括程序，所述程序能够被处理器执行以实现如权利要求1-7中任一项所述的方法。
9.一种特定蛋白浓度的检测装置，其特征在于，包括:
反应容器，用于容纳待测样本；
光源，用于向所述反应容器内的待测样本提供激光；
光信号接收器，用于采集激光经待测样本形成的透射光，并将透射光的光信号转换为电信号；
存储装置，用于存储各已知浓度样本的特定蛋白浓度和预先基于透射比浊法得到的对应的特定蛋白反应曲线；
数据处理装置，与所述光信号接收器和存储装置分别信号连接，用于：
根据所述电信号获取待测样本的特定蛋白反应曲线，所述特定蛋白反应曲线用于表示采样总时间内得到的电压的变化；
获取所述待测样本的特定蛋白反应曲线的一阶导数的最大值；
根据所述一阶导数的最大值与预设阈值的关系，确定所述待测样本的特定蛋白反应曲线是否存在波动，如果不存在波动，在所述采样总时间内选取提取时间段，如果存在波动，则获取所述待测样本的特定蛋白反应曲线的二阶导数，并根据所述二阶导数确定采样总时间内中的采样有效时间，在所述采样有效时间内选取提取时间段；所述根据所述二阶导数确定采样总时间内中的采样有效时间，包括：获取所述采样总时间内二阶导数大于零的至少一个时间段，将所述至少一个时间段中长度最大的时间段作为采样有效时间；
提取所述待测样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征；
根据预先基于透射比浊法得到的多个已知浓度样本的特定蛋白反应曲线，提取各已知浓度样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征；
建立各已知浓度样本的特定蛋白反应曲线的曲线特征与其对应的特定蛋白浓度间的定标函数关系式；
将所述待测样本的曲线特征输入所述定标函数关系式中，得到所述待测样本的特定蛋白浓度。
CN202011450007.5A 2020-12-09 2020-12-09 一种特定蛋白浓度的检测方法及检测装置 Active CN112710636B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011450007.5A CN112710636B (zh) 2020-12-09 2020-12-09 一种特定蛋白浓度的检测方法及检测装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011450007.5A CN112710636B (zh) 2020-12-09 2020-12-09 一种特定蛋白浓度的检测方法及检测装置

Publications (2)

Publication Number Publication Date
CN112710636A CN112710636A (zh) 2021-04-27
CN112710636B true CN112710636B (zh) 2022-05-24

Family

ID=75543090

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011450007.5A Active CN112710636B (zh) 2020-12-09 2020-12-09 一种特定蛋白浓度的检测方法及检测装置

Country Status (1)

Country Link
CN (1) CN112710636B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112964673B (zh) * 2021-01-29 2024-02-27 深圳市科曼医疗设备有限公司 特定蛋白反应曲线的异常识别方法及识别装置

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6650719B1 (en) * 2000-07-19 2003-11-18 Tektronix, Inc. MPEG PCR jitter, frequency offset and drift rate measurements
CN1508534A (zh) * 2002-10-21 2004-06-30 ����ɨ�����޹�˾ 缩短终点型反应分布图的分析时间的方法
CN101489035A (zh) * 2008-01-16 2009-07-22 三洋电机株式会社 摄像装置及抖动校正方法
CN105466927A (zh) * 2014-07-02 2016-04-06 深圳迈瑞生物医疗电子股份有限公司 一种比浊法异常反应曲线的识别、修正及报警方法
CN106442682A (zh) * 2016-08-31 2017-02-22 浙江大学 一种钢筋混凝土结构中钢筋析氢临界电流密度的测定方法
JP2018192247A (ja) * 2017-05-12 2018-12-06 キヤノンメディカルシステムズ株式会社 医用画像診断装置
CN109239360A (zh) * 2018-09-14 2019-01-18 深圳开立生物医疗科技股份有限公司 一种反应曲线异常检测方法及装置
CN110967507A (zh) * 2018-09-28 2020-04-07 希森美康株式会社 血液凝固分析方法及装置、及包含可执行程序的存储介质

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6650719B1 (en) * 2000-07-19 2003-11-18 Tektronix, Inc. MPEG PCR jitter, frequency offset and drift rate measurements
CN1508534A (zh) * 2002-10-21 2004-06-30 ����ɨ�����޹�˾ 缩短终点型反应分布图的分析时间的方法
CN101489035A (zh) * 2008-01-16 2009-07-22 三洋电机株式会社 摄像装置及抖动校正方法
CN105466927A (zh) * 2014-07-02 2016-04-06 深圳迈瑞生物医疗电子股份有限公司 一种比浊法异常反应曲线的识别、修正及报警方法
CN106442682A (zh) * 2016-08-31 2017-02-22 浙江大学 一种钢筋混凝土结构中钢筋析氢临界电流密度的测定方法
JP2018192247A (ja) * 2017-05-12 2018-12-06 キヤノンメディカルシステムズ株式会社 医用画像診断装置
CN109239360A (zh) * 2018-09-14 2019-01-18 深圳开立生物医疗科技股份有限公司 一种反应曲线异常检测方法及装置
CN110967507A (zh) * 2018-09-28 2020-04-07 希森美康株式会社 血液凝固分析方法及装置、及包含可执行程序的存储介质

Also Published As

Publication number Publication date
CN112710636A (zh) 2021-04-27

Similar Documents

Publication Publication Date Title
Woimant et al. New tools for Wilson’s disease diagnosis: Exchangeable copper fraction
Hiemann et al. Challenges of automated screening and differentiation of non-organ specific autoantibodies on HEp-2 cells
US11761965B2 (en) Methods for diagnosing systemic lupus erythematosus
CN112710627B (zh) 一种特定蛋白浓度的检测方法及检测装置
WO2022100040A1 (zh) 一种特定蛋白反应检测方法和装置
CN112710636B (zh) 一种特定蛋白浓度的检测方法及检测装置
Nitschke et al. J-edited DIffusional proton nuclear magnetic resonance spectroscopic measurement of glycoprotein and supramolecular phospholipid biomarkers of inflammation in human serum
Srichomkwun et al. Diagnostic dilemma in discordant thyroid function tests due to thyroid hormone autoantibodies
JP6281973B2 (ja) 混合物試料の特性を表現する方法、混合物試料の特性を評価する方法、混合物試料の属性を識別する方法、及び混合物試料に由来する電磁波信号を処理する方法
Trinks et al. Omics-based biomarkers as useful tools in metabolic dysfunction-associated steatotic liver disease clinical practice: How far are we?
CN109030801B (zh) 一种临床样本自动生化分析仪
US11747328B2 (en) Method for measuring the plasma concentration of an analyte directly on a whole blood sample
CN114578040B (zh) 样本检测方法、样本分析仪及其控制装置
CN112485440B (zh) 特定蛋白质反应检测方法、蛋白质检测装置和定标方法
CN112964673B (zh) 特定蛋白反应曲线的异常识别方法及识别装置
Bertin et al. ICARE improves antinuclear antibody detection by overcoming the barriers preventing accreditation
Kopke et al. Response to letter to editor regarding Variability of symmetric dimethylarginine in apparently healthy dogs
CN112485439A (zh) 特定蛋白质反应检测方法、蛋白质检测装置和定标方法
Sturgess et al. Improving the effectiveness of autoantibody testing in the clinic
Ahmed et al. The clinical aspect of overweight on rheumatoid arthritis and disease activity
Nurakhova et al. The experience of application of apparatus “aklides” in Kazakhstan for the diagnosis of systemic autoimmune diseases
CN110082528B (zh) 诊断或预断人类口腔癌的系统和应用
Xiang et al. Clinical value of chemiluminescence method for detection of antinuclear antibody profiles
Wei et al. Method comparison of antibody titre sensitivity using RCPAQAP datasets
Allen Development of a Magnetic Resonance Elastography-Based Model for Diagnosis of Nonalcoholic Steatohepatitis

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