CN112710627B - 一种特定蛋白浓度的检测方法及检测装置 - Google Patents
一种特定蛋白浓度的检测方法及检测装置 Download PDFInfo
- Publication number
- CN112710627B CN112710627B CN202011449118.4A CN202011449118A CN112710627B CN 112710627 B CN112710627 B CN 112710627B CN 202011449118 A CN202011449118 A CN 202011449118A CN 112710627 B CN112710627 B CN 112710627B
- 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
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/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/39—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using tunable lasers
-
- 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/75—Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated
- G01N21/77—Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated by observing the effect on a chemical indicator
- G01N21/82—Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated by observing the effect on a chemical indicator producing a precipitate or turbidity
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Engineering & Computer Science (AREA)
- Optics & Photonics (AREA)
- Plasma & Fusion (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壁照射到反应液中的微球上,发生散射,散射光从非90度(与反应容器10侧壁)的角度进入反应容器10壁,再通过反应容器10壁进入空气,最后到光信号接收器30。整个过程散射光分别经过了反应液、反应容器10和空气三种不同的介质,因为散射光从非90度的方向经过这三种介质,所以散射光在传播过程中会发生折射。
特定蛋白反应检测原理是依据抗原、抗体在特定的电解质溶液中反应,快速形成免疫复合物微粒,使反应液出现浊度,随着时间的推移,复合物微粒聚合得就会越大,浊度也逐渐变大。在散射比浊法中,光信号接收器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)的最小值minFi’(tb),如果minFi’(tb)不小于零,则第i个已知浓度样本的特定蛋白反应曲线不存在波动,否则,第i个已知浓度样本的特定蛋白反应曲线存在波动,那么就再次对该已知浓度样本进行检测,直到该已知浓度样本的特定蛋白反应曲线的一阶导数不小于0。
数据处理装置50与光信号接收器30和存储装置40分别信号连接,其用于根据电信号获取待测样本的特定蛋白反应曲线,特定蛋白反应曲线用于表示采样总时间内得到的电压的变化。
获取到的待测样本的特定蛋白反应曲线的公式,包括:
Va=F(ta)。
其中,0≤ta≤Ta,ta∈实数,ta为待测样本的采样时间,Va为采样时获取的电压,Ta为采样总时间的结束时间点。
而后获取待测样本的特定蛋白反应曲线的一阶导数F’(ta)的最小值minF’(ta),根据一阶导数的最小值与预设阈值的关系,确定待测样本的特定蛋白反应曲线是否存在波动。一些实施例中,该预设阈值为零,如果minF’(ta)大于等于零,代表待测样本的特定蛋白反应曲线为一增函数,待测样本的特定蛋白反应曲线不存在波动,如果minF’(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的计算公式包括:
其中,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的计算公式包括:
其中,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的计算公式包括:
其中,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)的最小值minF’(ta)。
步骤300、根据一阶导数的最小值与预设阈值的关系,确定待测样本的特定蛋白反应曲线是否存在波动,如果不存在波动,执行步骤400,如果存在波动,执行步骤500。
一些实施例中,该预设阈值为零,如果minF’(ta)大于等于零,代表待测样本的特定蛋白反应曲线为一增函数,待测样本的特定蛋白反应曲线不存在波动,如果minF’(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的计算公式包括:
其中,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的计算公式包括:
其中,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)的最小值minFi’(tb),如果minFi’(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的计算公式包括:
其中,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 (8)
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的计算公式包括:
其中,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的计算公式包括:
其中,ta1为所述提取时间段的起始时间点,ta2为所述提取时间段的结束时间点, t1≤ta1<ta2≤t2, t1为采样有效时间的起始时间点,t2为采样有效时间的结束时间点,将所述面积Sa2作为所述待测样本的特定蛋白反应曲线的曲线特征。
6.如权利要求1所述的方法,其特征在于,基于散射比浊法得到多个已知浓度样本的特定蛋白反应曲线的方式,包括:
对于任一已知浓度样本,基于散射比浊法获取所述已知浓度样本的特定蛋白反应曲线,判断所述已知浓度样本的特定蛋白反应曲线是否存在波动,如果存在波动,则对存在波动的特定蛋白反应曲线的已知浓度样本重新检测,以再次获取相应的特定蛋白反应曲线,直到每个已知浓度样本的特定蛋白反应曲线均不存在波动。
7.如权利要求6所述的方法,其特征在于,所述判断所述已知浓度样本的特定蛋白反应曲线是否存在波动,包括:
获取所述已知浓度样本的特定蛋白反应曲线的一阶导数的最小值;
如果所述已知浓度样本的特定蛋白反应曲线的一阶导数的最小值不小于零,则所述已知浓度样本的特定蛋白反应曲线不存在波动,否则,所述已知浓度样本的特定蛋白反应曲线存在波动。
8.一种特定蛋白浓度的检测装置,其特征在于,包括:
反应容器,用于容纳待测样本;
光源,用于向所述反应容器内的待测样本提供激光;
光信号接收器,用于采集激光经待测样本形成的散射光,并将散射光的光信号转换为电信号;
存储装置,用于存储各已知浓度样本的特定蛋白浓度和预先基于散射比浊法得到的对应的特定蛋白反应曲线;
数据处理装置,与所述光信号接收器和存储装置分别信号连接,用于:
根据所述电信号获取待测样本的特定蛋白反应曲线,所述特定蛋白反应曲线用于表示采样总时间内得到的电压的变化;
获取所述待测样本的特定蛋白反应曲线的一阶导数的最小值;
根据所述一阶导数的最小值与预设阈值的关系,确定所述待测样本的特定蛋白反应曲线是否存在波动,如果不存在波动,在所述采样总时间内选取提取时间段,如果存在波动,则获取所述待测样本的特定蛋白反应曲线的二阶导数,并根据所述二阶导数确定采样总时间内中的采样有效时间,在所述采样有效时间内选取提取时间段;其中, 根据二阶导数确定采样有效时间,包括:获取所述采样总时间内二阶导数小于零的至少一个时间段;将所述至少一个时间段中长度最大的时间段作为采样有效时间;
提取所述待测样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征;
根据预先基于散射比浊法得到的多个已知浓度样本的特定蛋白反应曲线,提取各已知浓度样本的特定蛋白反应曲线在所述提取时间段对应的曲线特征;
建立各已知浓度样本的特定蛋白反应曲线的曲线特征与其对应的特定蛋白浓度间的定标函数关系式;
将所述待测样本的曲线特征输入所述定标函数关系式中,得到所述待测样本的特定蛋白浓度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011449118.4A CN112710627B (zh) | 2020-12-09 | 2020-12-09 | 一种特定蛋白浓度的检测方法及检测装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011449118.4A CN112710627B (zh) | 2020-12-09 | 2020-12-09 | 一种特定蛋白浓度的检测方法及检测装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112710627A CN112710627A (zh) | 2021-04-27 |
CN112710627B true CN112710627B (zh) | 2022-05-27 |
Family
ID=75543087
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011449118.4A Active CN112710627B (zh) | 2020-12-09 | 2020-12-09 | 一种特定蛋白浓度的检测方法及检测装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112710627B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112964673B (zh) * | 2021-01-29 | 2024-02-27 | 深圳市科曼医疗设备有限公司 | 特定蛋白反应曲线的异常识别方法及识别装置 |
CN113466180B (zh) * | 2021-09-02 | 2021-11-16 | 天津迈科隆生物科技有限公司 | 一种特定蛋白检测方法、电子设备及计算机可读存储介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102428373A (zh) * | 2009-05-20 | 2012-04-25 | 株式会社日立高新技术 | 自动分析装置以及分析方法 |
CN102866256A (zh) * | 2012-04-23 | 2013-01-09 | 深圳市希莱恒医用电子有限公司 | 一种超敏c反应蛋白检测方法及检测试剂 |
CN105466927A (zh) * | 2014-07-02 | 2016-04-06 | 深圳迈瑞生物医疗电子股份有限公司 | 一种比浊法异常反应曲线的识别、修正及报警方法 |
CN106442682A (zh) * | 2016-08-31 | 2017-02-22 | 浙江大学 | 一种钢筋混凝土结构中钢筋析氢临界电流密度的测定方法 |
CN107703132A (zh) * | 2017-09-30 | 2018-02-16 | 深圳迈瑞生物医疗电子股份有限公司 | 反应曲线异常处理方法及装置、生化分析仪、存储介质 |
CN109239360A (zh) * | 2018-09-14 | 2019-01-18 | 深圳开立生物医疗科技股份有限公司 | 一种反应曲线异常检测方法及装置 |
CN109977916A (zh) * | 2019-04-09 | 2019-07-05 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种基于嵌入式平台的变压器油中溶解气体在线监测系统 |
CN112014335A (zh) * | 2019-05-30 | 2020-12-01 | 深圳迈瑞生物医疗电子股份有限公司 | 基于比浊法测量分析物的方法、样本分析仪及存储介质 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5897323B2 (ja) * | 2011-12-26 | 2016-03-30 | 株式会社日立ハイテクノロジーズ | 自動分析装置および測定値異常検出方法 |
-
2020
- 2020-12-09 CN CN202011449118.4A patent/CN112710627B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102428373A (zh) * | 2009-05-20 | 2012-04-25 | 株式会社日立高新技术 | 自动分析装置以及分析方法 |
CN102866256A (zh) * | 2012-04-23 | 2013-01-09 | 深圳市希莱恒医用电子有限公司 | 一种超敏c反应蛋白检测方法及检测试剂 |
CN105466927A (zh) * | 2014-07-02 | 2016-04-06 | 深圳迈瑞生物医疗电子股份有限公司 | 一种比浊法异常反应曲线的识别、修正及报警方法 |
CN106442682A (zh) * | 2016-08-31 | 2017-02-22 | 浙江大学 | 一种钢筋混凝土结构中钢筋析氢临界电流密度的测定方法 |
CN107703132A (zh) * | 2017-09-30 | 2018-02-16 | 深圳迈瑞生物医疗电子股份有限公司 | 反应曲线异常处理方法及装置、生化分析仪、存储介质 |
CN109239360A (zh) * | 2018-09-14 | 2019-01-18 | 深圳开立生物医疗科技股份有限公司 | 一种反应曲线异常检测方法及装置 |
CN109977916A (zh) * | 2019-04-09 | 2019-07-05 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种基于嵌入式平台的变压器油中溶解气体在线监测系统 |
CN112014335A (zh) * | 2019-05-30 | 2020-12-01 | 深圳迈瑞生物医疗电子股份有限公司 | 基于比浊法测量分析物的方法、样本分析仪及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112710627A (zh) | 2021-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Woimant et al. | New tools for Wilson’s disease diagnosis: Exchangeable copper fraction | |
CN109239360B (zh) | 一种反应曲线异常检测方法及装置 | |
Willis et al. | Immune biomarkers: the promises and pitfalls of personalized medicine | |
CN112710627B (zh) | 一种特定蛋白浓度的检测方法及检测装置 | |
US20190025326A1 (en) | Methods for Diagnosing Systemic Lupus Erythematosus | |
WO2022100040A1 (zh) | 一种特定蛋白反应检测方法和装置 | |
Jasensky et al. | Evaluation of three different point‐of‐care tests for quantitative measurement of canine C‐reactive protein | |
Yoo et al. | Performance of an automated fluorescence antinuclear antibody image analyzer | |
CN112710636B (zh) | 一种特定蛋白浓度的检测方法及检测装置 | |
Craig et al. | Serum C‐reactive protein and S100A12 concentrations in dogs with hepatic disease | |
Srichomkwun et al. | Diagnostic dilemma in discordant thyroid function tests due to thyroid hormone autoantibodies | |
Sahn et al. | Histological features of ileitis differentiating pediatric Crohn disease from ulcerative colitis with backwash ileitis | |
CN112014335A (zh) | 基于比浊法测量分析物的方法、样本分析仪及存储介质 | |
CN109030801B (zh) | 一种临床样本自动生化分析仪 | |
CN106872379B (zh) | 一种尿液成分测量分析仪 | |
CN112485440B (zh) | 特定蛋白质反应检测方法、蛋白质检测装置和定标方法 | |
CN114578040B (zh) | 样本检测方法、样本分析仪及其控制装置 | |
CN112964673B (zh) | 特定蛋白反应曲线的异常识别方法及识别装置 | |
Macias‐Muñoz et al. | Analytical and clinical evaluation of DiaSorin Liaison® Calprotectin fecal assay adapted for serum samples | |
KR102593221B1 (ko) | 분석 장치 및 분석 방법 | |
Nagy et al. | Anti-neutrophil cytoplasmic antibody testing by indirect immunofluorescence: computer-aided versus conventional microscopic evaluation of routine diagnostic samples from patients with vasculitis or other inflammatory diseases | |
CN114839140A (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 | |
CN106290817B (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 |