WO2020103051A1 - 样本吸光度差的测量方法、样本分析仪和存储介质 - Google Patents

样本吸光度差的测量方法、样本分析仪和存储介质

Info

Publication number
WO2020103051A1
WO2020103051A1 PCT/CN2018/116778 CN2018116778W WO2020103051A1 WO 2020103051 A1 WO2020103051 A1 WO 2020103051A1 CN 2018116778 W CN2018116778 W CN 2018116778W WO 2020103051 A1 WO2020103051 A1 WO 2020103051A1
Authority
WO
WIPO (PCT)
Prior art keywords
difference
absorbance
weight
sample
absolute value
Prior art date
Application number
PCT/CN2018/116778
Other languages
English (en)
French (fr)
Inventor
李聪
郭文恒
李坷坷
Original Assignee
北京普利生仪器有限公司
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 北京普利生仪器有限公司 filed Critical 北京普利生仪器有限公司
Priority to CN201880099067.1A priority Critical patent/CN112912712B/zh
Priority to PCT/CN2018/116778 priority patent/WO2020103051A1/zh
Publication of WO2020103051A1 publication Critical patent/WO2020103051A1/zh

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/27—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration

Definitions

  • the automatic coagulation analyzer generally uses immunoturbidimetry to test the concentration of D-dimer (DD), fibrin (pro) degradation products (FDP) and other items, and the most basic biochemical reaction calculation method has three Species: endpoint method (1-Point), 2-point method (2-Point) and rate method (Rate).
  • the endpoint method is: the instrument only detects the absorbance value at a certain time point of the biochemical reaction, which is susceptible to interference, and is rarely used in automatic coagulation analyzers;
  • the two-point method is: the instrument detects the absorbance value at two time points of the biochemical reaction , The absorbance value at the second time point minus the absorbance value at the first time point to obtain the difference in absorbance (absorbance difference).
  • Rate method The instrument continuously monitors the change in absorbance caused by changes in the substrate or product content during the biochemical reaction, obtains the rate of change in absorbance, and determines the difference in absorbance based on the rate of change in absorbance.
  • the key of the rate method is to accurately describe the relationship between the reaction rate and time.
  • the signal-to-noise ratio is relatively low during the detection of low-value samples. The rate of continuous monitoring is easily interfered, and it is difficult to accurately reflect the biochemical reaction process, resulting in low value repeatability bad.
  • the embodiments of the present application provide a method for measuring a sample absorbance difference, a sample analyzer, and a storage medium, to at least solve the problem of measuring repeatability when measuring absorbance difference, when the absorbance difference is low and the absorbance difference is high.
  • Technical issues with linear measurement range are provided.
  • a method for measuring a sample absorbance difference including; acquiring absorbance data of a sample object; calculating the first absorbance difference using the two-point method according to the absorbance data, and calculating the Second absorbance difference; substitute the first absorbance difference and the second absorbance difference into the preset weighted average function used to calculate the final absorbance difference; determine the weight and the second absorbance of the first absorbance difference satisfying the weighted average function according to the preset functional relationship
  • the method further includes: determining the corresponding to the final absorbance difference according to the preset correspondence Sample concentration, where the preset correspondence is the correspondence between absorbance difference and sample concentration.
  • the preset functional relationship includes: within the set definition domain, if the absorbance difference is less than the set threshold, the weight of the first absorbance difference is greater than the weight of the second absorbance difference; if the absorbance difference is greater than the set threshold, the first The weight of one difference in absorbance is less than the weight of the difference in second absorbance.
  • the absolute value of the difference between the weight of the first absorbance difference and the weight of the difference in the second absorbance is negatively related to the difference in absorbance; if the difference in absorbance is greater than the set threshold, the absolute value of the difference in weight Positively correlated with the difference in absorbance.
  • the set threshold is set by the user.
  • the speed of change of the absolute value of the weight difference is set by the user.
  • the preset functional relationship further includes: within the set definition domain, if the absorbance difference is equal to the set threshold, the weight of the first absorbance difference and the second absorbance difference are the same.
  • the method is applied to a sample analyzer.
  • the absolute value of the difference between the weight of the first absorbance difference and the weight of the difference in the second absorbance is negatively related to the difference in absorbance; Positively correlated with the difference in absorbance.
  • the absolute value of the weighted phase difference increases as the absorbance difference increases, and the rate of increase in the absolute value of the weighted phase difference is inversely related to the absorbance difference.
  • the set threshold is set by the user.
  • one or more non-volatile computer-readable storage media are provided on which a computer program is stored, and when the computer program is executed by a processor, the above-mentioned method for measuring the difference in absorbance of a sample is realized .
  • the two-point method is adopted by the low value part to obtain better low value repeatability; the rate method is used to obtain the monotonicity of the calculated value of the difference in absorbance and the higher linear measurement range; the low value and the high value Partially well-connected, there is no technical effect of faults or jumps; thus solving the technical problem of measuring the repeatability of the measurement when the absorbance difference is low, and the linear measurement range when the absorbance difference is high.
  • FIG. 1a is a schematic structural diagram of a sample analyzer according to an embodiment of the present application.
  • FIG. 1b is a schematic flowchart of an optional method for measuring absorbance of a sample according to an embodiment of the present application
  • FIG. 2 is a flowchart of calculating the first absorbance difference by the two-point method provided in the embodiment of the present application;
  • FIG. 3 is a flowchart of calculating the second difference in absorbance by the rate method provided by the embodiment of the present application.
  • FIG. 4 is a flow chart of another rate method for calculating a second difference in absorbance provided by an embodiment of the present application.
  • FIG. 5 is a schematic diagram of weight changes of the two-point method and the rate method determined when the preset function is based on a hyperbolic tangent function according to an embodiment of the present application;
  • FIG. 6 is a schematic diagram of the weight change of the two-point method and the rate method determined when the preset function is based on a first-order system according to an embodiment of the present application;
  • FIG. 9 is a schematic flowchart of a method for measuring a difference in absorbance of a sample according to another embodiment of the present application.
  • Figures 10a-10d are schematic diagrams of the repetitive calculation results of the two-point method and the rate method for calculating the low value part of the difference in absorbance;
  • FIG. 12 is a flowchart of a working method corresponding to a sample analyzer provided by an embodiment of the present application.
  • Rate method Also known as dynamic method and kinetic method, the instrument continuously monitors the absorbance change caused by the change of the substrate or product content during the biochemical reaction, obtains the rate of change in absorbance, and finds the moment when the reaction rate is maximum within the specified sampling time. It is extended to both sides with the center as the center to find an interval that approximately meets the linear reaction rate, and the absorbance difference in this interval is calculated as the absorbance difference of the entire reaction process.
  • Two-point method The absorbance values at two time points of the biochemical reaction are detected by the instrument, which are the sampling start point and the sampling end point, respectively.
  • the absorbance difference of the method is obtained by subtracting the absorbance of the sampling start point from the absorbance of the sampling end point.
  • the two-point method is only suitable for the determination of the linear phase of the reaction rate. It cannot accurately describe the process that the entire reaction rate is continuously decreasing as the substrate is continuously consumed.
  • the main difference between the above two-point method and the rate method in describing the reaction process is that the two-point method is based on that the reaction rate does not change with the consumption of the substrate (expressed as the test time), and the rate method is based on the reaction rate. Changes with consumption.
  • Intensity f (time), where the above Intensity represents the luminous flux data of the scattered light and / or transmitted light after the beam emitted by the light source illuminates the sample, and f (time) represents the luminous flux data corresponding to the time information; it can be based on the luminous flux data and the time information Corresponding relationship, draw the luminous flux curve.
  • Step S206 Determine the first difference in absorbance based on the absorbance curve.
  • dOD Twopoint represents a first difference in absorbance
  • Abs (EndPoint) indicates the termination time point corresponding to the absorbance data
  • Abs (StartPoint) represents an initial time point absorbance data.
  • Step S306 find an interval approximately satisfying the linear reaction process near MaxPoint, and record the start and end points of the interval as Cut_Min and Cut_Max;
  • step S408 the Abs_fit curve is derived, the tangent equation is obtained, and the MaxPoint time point is obtained;
  • Step S414 obtaining the tangent equation Abs_cut and the limit points of the linear range: Cut_Min and Cut_Max;
  • Step S106 substituting the first absorbance difference and the second absorbance difference into a preset weighted average function (which may also be called a weight allocation function or a proportional allocation function) used to calculate the final absorbance difference;
  • a preset weighted average function which may also be called a weight allocation function or a proportional allocation function
  • the preset functional relationship also needs to be considered: within the set definition domain, if the absorbance difference is less than a set threshold, the weight of the first absorbance difference is greater than the first The weight of the second absorbance difference; if the absorbance difference is greater than the set threshold, the weight of the first absorbance difference is less than the weight of the second absorbance difference.
  • the set threshold can be set by the user.
  • the axis is the absorbance difference
  • the vertical axis is the weight (proportion)
  • dOD0 is 5000
  • the absolute value of the difference between the weight of the first absorbance difference and the weight of the second absorbance difference becomes larger gradually, so the absolute value of the difference between the weight of the first absorbance difference and the weight of the second absorbance difference is negatively correlated with the absorbance difference .
  • the absolute value of the weighted phase difference decreases as the absorbance difference increases, and the absolute value of the weighted phase difference The reduction speed is positively correlated with the difference in absorbance.
  • the absorbance difference is greater than the set threshold, the absolute value of the weighted phase difference increases as the absorbance difference increases, and the absolute value of the weighted phase difference The increase speed is inversely related to the difference in absorbance.
  • FIG. 7 is a schematic diagram of changes in the weights of the first difference in absorbance and the second difference in absorbance determined by the preset function.
  • the horizontal axis is the absorbance difference
  • the vertical axis is the ratio of the weight of the first absorbance difference corresponding to the two-point method and the second absorbance difference corresponding to the rate method.
  • the preset function can make the total absorbance difference infinitely close to the calculation result of the two-point method, and make full use of the linear fitting process of the two-point method to resist noise interference and improve the signal-to-noise ratio for better Low value repeatability; in the high value part of the absorbance difference, the preset function can make the total absorbance difference infinitely close to the calculation result of the rate method, make full use of the rate method to actively reflect the change of the reaction rate with time to obtain the difference in absorbance Monotonicity of calculated values and higher linear measurement range.
  • Step S8044 Find an interval that satisfies the linear response process approximately near MaxPoint, and record the start and end points of the interval as Cut_Min and Cut_Max;
  • Step S9066 within the analysis interval, an N-order polynomial fits the absorbance curve (N ⁇ 2), and obtains: Abs fit ;
  • Figure 11 shows that when the difference in absorbance is high, the two-point method, the rate method, and the use of this application A comparison diagram of the difference in absorbance calculated by the measurement method of the difference in absorbance of the sample.
  • the concentration of DD in Figure 11 represents the concentration of DD (D-dimer); the difference in absorbance calculated by the two-point method tends to be saturated or even decreased, and the difference in absorbance calculated by the rate method remains monotonously increasing, effectively improving the linear measurement
  • the linear measurement range of the algorithm used in the embodiment of the present application is basically the same as the rate method.
  • the present application also provides a sample analyzer, and the corresponding working method flowchart is shown in FIG. 12.
  • the working method includes the following steps:
  • the disclosed technical content may be implemented in other ways.
  • the device embodiments described above are only schematic.
  • the division of units may be a division of logical functions.
  • multiple units or components may be combined or integrated into Another system, or some features can be ignored, or not implemented.
  • the displayed or discussed mutual coupling or direct coupling or communication connection may be indirect coupling or communication connection through some interfaces, units or modules, and may be in electrical or other forms.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

一种样本吸光度差的测量方法、样本分析仪和存储介质。该方法包括以下步骤：获取样本对象的吸光度数据；根据吸光度数据应用两点法计算第一吸光度差，根据吸光度数据应用速率法计算第二吸光度差；将第一吸光度差和第二吸光度差代入用于计算最终吸光度差的预设的加权平均函数；根据预设函数关系确定满足所述加权平均函数的第一吸光度差的权重和第二吸光度差的权重以及最终吸光度差。还具有样本分析仪，包括光源（132）、检测器（134）、处理器（136）。

Description

样本吸光度差的测量方法、样本分析仪和存储介质 技术领域
本申请涉及光学测量领域，具体而言，涉及一种样本吸光度差的测量方法、样本分析仪和存储介质。
背景技术
现有技术中，全自动凝血分析仪一般采用免疫比浊法测试D-二聚体(DD)、纤维蛋白(原)降解产物(FDP)等项目的浓度，而最基本的生化反应计算方法有三种：终点法(1-Point assay)、两点法(2-Point assay)和速率法(Rate assay)。其中，终点法为：仪器仅检测生化反应的某一时间点的吸光度值，容易受到干扰、很少用于全自动凝血分析仪；两点法为：仪器检测生化反应两个时间点的吸光度值，第二时间点的吸光度值减去第一时间点的吸光度值得到吸光度的差值即(吸光度差)。两点法仅适用于反应速率的线性期测定，但是随着底物的不断消耗、整个反应速率是不断减小的，因此两点法存在可测量的最高限度，超过限度后会出现吸光度差饱和、甚至降低，进而限制了高值部分的线性能力。速率法：仪器连续监测生化反应过程中底物或产物含量变化导致的吸光度变化，获取吸光度变化速率，依据吸光度变化速率确定吸光度差。速率法的关键在于准确的描述反应速率随时间的变化关系，在低值样本检测过程中信噪比相对较低，连续监测的速率容易受到干扰，难以准确反映生化反应过程，导致低值重复性不良。
当前，全自动凝血分析仪大多都支持两点法和速率法，但是，无论选用了何种计算方法都不能同时满足低值重复性和高值线性。
针对上述问题，目前尚未提出有效的解决方案。
发明内容
本申请实施例提供了一种样本吸光度差的测量方法、样本分析仪和存储介质，以至少解决在测量吸光度差时，无法同时兼顾吸光度差低值时的测量重复性以及吸光度差高值时的线性测量范围的技术问题。
根据本申请实施例的一个方面，提供了一种样本吸光度差的测量方法，包括；获 取样本对象的吸光度数据；根据吸光度数据应用两点法计算第一吸光度差，根据吸光度数据应用速率法计算第二吸光度差；将第一吸光度差和第二吸光度差代入用于计算最终吸光度差的预设的加权平均函数；根据预设函数关系确定满足加权平均函数的第一吸光度差的权重和第二吸光度差的权重以及最终吸光度差，其中，预设函数关系为针对第一吸光度差的权重、第二吸光度差的权重和吸光度差的函数关系。
可选地，根据预设函数关系确定满足加权平均函数的第一吸光度差的权重和第二吸光度差的权重以及最终吸光度差之后，方法还包括：依据预设对应关系确定与最终吸光度差对应的样本浓度，其中，预设对应关系为吸光度差与样本浓度的对应关系。
可选地，预设函数关系包括：在设定的定义域内，如果吸光度差小于设定阈值，则第一吸光度差的权重大于第二吸光度差的权重；如果吸光度差大于设定阈值，则第一吸光度差的权重小于第二吸光度差的权重。
可选地，如果吸光度差小于设定阈值，第一吸光度差的权重与第二吸光度差的权重相差值的绝对值与吸光度差负相关；如果吸光度差大于设定阈值，权重相差值的绝对值与吸光度差正相关。
可选地，如果吸光度差小于设定阈值，权重相差值的绝对值随着吸光度差的增大而减小，权重相差值的绝对值的减小速度和吸光度差正相关。
可选地，如果吸光度差大于设定阈值，权重相差值的绝对值随着吸光度差的增大而增大，权重相差值的绝对值的增大速度和吸光度差负相关。
可选地，设定阈值由用户设定。
可选地，权重相差值的绝对值的变化速度由用户设定。
可选地，预设函数关系还包括：在设定的定义域内，如果吸光度差等于设定阈值时，则第一吸光度差的权重和第二吸光度差的权重是相同的。
可选地，加权平均函数为：F(x)＝a·x+b·(c-x)+d；其中，a为第二吸光度差，x为第二吸光度差的权重，b为第一吸光度差，(c-x)为第一吸光度差的权重，c为大于或等于x的常数，d为大于或等于0的常数，F(x)为最终吸光度差。
可选地，方法应用于样本分析仪。
根据本申请实施例的一个方面，提供了一种样本分析仪，包括：光源、检测器和处理器；其中，光源用于发出用于照射样本的光束；检测器用于检测在光束照射样本后产生的光通量数据；处理器运行程序，其中，程序运行时对于从检测器输出的数据 执行如下处理步骤：根据光通量数据计算吸光度数据；根据吸光度数据应用两点法计算第一吸光度差，根据吸光度数据应用速率法计算第二吸光度差；将第一吸光度差和第二吸光度差代入用于计算最终吸光度差的预设的加权平均函数；根据预设函数关系确定满足加权平均函数的第一吸光度差的权重和第二吸光度差的权重以及最终吸光度差，其中，预设函数关系为针对第一吸光度差的权重、第二吸光度差的权重和吸光度差的函数关系。
可选地，处理器还用于：依据预设对应关系确定与最终吸光度差对应的样本浓度，其中，预设对应关系为吸光度差与样本浓度的对应关系。
可选地，预设函数关系包括：在设定的定义域内，如果吸光度差小于设定阈值，则第一吸光度差的权重大于第二吸光度差的权重；如果吸光度差大于设定阈值，则第一吸光度差的权重小于第二吸光度差的权重。
可选地，如果吸光度差小于设定阈值，第一吸光度差的权重与第二吸光度差的权重相差值的绝对值与吸光度差负相关；如果吸光度差大于设定阈值，权重相差值的绝对值与吸光度差正相关。
可选地，如果吸光度差小于设定阈值，权重相差值的绝对值随着吸光度差的增大而减小，权重相差值的绝对值的减小速度和吸光度差正相关。
可选地，如果吸光度差大于设定阈值，权重相差值的绝对值随着吸光度差的增大而增大，权重相差值的绝对值的增大速度和吸光度差负相关。
可选地，设定阈值由用户设定。
可选地，权重相差值的绝对值的变化速度由用户设定。
可选地，加权平均函数为：F(x)＝a·x+b·(c-x)+d；其中，a为第二吸光度差，x为第二吸光度差的权重，b为第一吸光度差，c-x为第一吸光度差的权重，c为大于或等于x的常数，d为大于或等于0的常数，F(x)为最终吸光度差。
根据本申请实施例的一个方面，提供了一个或多个非易失性计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以上所述的样本吸光度差的测量方法。
在本申请实施例中，根据样本对象的吸光度数据应用两点法计算第一吸光度差，根据吸光度数据应用速率法计算第二吸光度差；并根据预设函数关系确定满足加权平均函数的第一吸光度差的权重和第二吸光度差的权重以及最终吸光度差，其中，预设函数关系为针对第一吸光度差的权重、第二吸光度差的权重和吸光度差的函数关系。 达到了通过低值部分采用两点法，以获得更好的低值重复性；高值部分采用速率法，以获得吸光度差计算值的单调性和更高的线性测量范围；低值和高值部分良好连接，不存在断层或跳变的技术效果；从而解决了在测量吸光度差时，无法同时兼顾吸光度差低值时的测量重复性以及吸光度差高值时的线性测量范围的技术问题。
附图说明
此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：
图1a是根据本申请实施例的一种样本分析仪的结构示意图；
图1b是根据本申请实施例的一种可选的测量样本吸光度方法的流程示意图；
图2为本申请实施例提供的两点法计算第一吸光度差的流程图；
图3为本申请实施例提供的速率法计算第二吸光度差的流程图；
图4为本申请实施例提供的另一种速率法计算第二吸光度差的流程图；
图5为本申请实施例提供的所述预设函数基于双曲正切函数时，确定的两点法与速率法的权重变化的示意图；
图6为本申请实施例提供的所述预设函数基于一阶***时，确定的两点法与速率法的权重变化的示意图；
图7为本申请实施例提供的由所述预设函数确定的第一吸光度差与第二吸光度差的权重的变化情况的示意图；
图8为本申请一实施例提供的样本吸光度差的测量方法流程示意图；
图9为本申请另一实施例提供的样本吸光度差的测量方法的流程示意图；
图10a-10d为采用两点法和速率法计算吸光度差低值部分的重复性计算结果曲线示意图；
图11为吸光度差为高值时，两点法、速率法、与使用本申请所述的样本吸光度差的测量方法计算的吸光度差的对比示意图；
图12为本申请实施例提供的样本分析仪对应的工作方法的流程图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分的实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本申请的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、***、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
为了更好地理解本申请实施例，以下将本申请实施例中涉及的术语简述如下：
速率法：又称之为动态法、动力学方法，通过仪器连续监测生化反应过程中底物或产物含量变化导致的吸光度变化，获取吸光度变化速率，寻找在指定采样时间内反应速率最大的时刻，并以其为中心向两侧扩展、寻找近似满足线性反应速率的区间，计算该区间的吸光度差作为整个反应过程的吸光度差。
速率法的关键在于准确的描述反应速率随时间的变化关系、进而准确的确定吸光度的变化速率。
两点法：通过仪器检测生化反应两个时间点的吸光度值，分别是采样起始点和采样终止点，以采样终止点的吸光度减去采样起始点的吸光度即为本方法获得吸光度差。两点法仅适用于反应速率的线性期测定，不能准确的描述随着底物的不断消耗、整个反应速率是不断减小的过程。
上述两点法与速率法对于反应过程的描述主要差异点在于：两点法基于反应速率不随底物的消耗(表现为测试时间的推移)而变化，速率法基于反应速率会随着底物的消耗而变化。
图1a是根据本申请实施例的一种样本分析仪的结构示意图，如图1a所示，该样本分析包括：光源132、检测器134、处理器136，其中，光源132用于发出用于照射样本的光束；
检测器134用于检测在光束照射样本后产生的光通量数据；
处理器136运行程序，其中，程序运行时对于从检测器输出的数据执行如下步骤：根据光通量数据计算吸光度数据；根据吸光度数据应用两点法计算第一吸光度差，根据吸光度数据应用速率法计算第二吸光度差；将第一吸光度差和第二吸光度差代入用于计算最终吸光度差的预设的加权平均函数；根据预设函数关系确定满足加权平均函数的第一吸光度差的权重和第二吸光度差的权重以及最终吸光度差；其中，所述预设函数关系为针对所述第一吸光度差的权重、第二吸光度差的权重和吸光度差的函数关系。
上述样本分析仪包括但不限于血液分析仪、生化分析仪、免疫分析仪、凝血分析仪等等体外诊断设备。
本申请实施例，提供了一种测量样本吸光度方法，该方法运行于上述样本分析仪。需要说明的是，在以下附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机***中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
图1b是根据本申请实施例的测量样本吸光度方法的流程示意图，如图1b所示，该方法至少包括步骤S102-S108，其中：
步骤S102，获取样本对象的吸光度数据；
在本申请的一些可选的实施例中，可通过执行以下步骤S1022-S1024获取样本对象的吸光度数据：步骤S1022，通过检测光通量的仪器获取光通量数据；此处光通量数据包括透射光和/或散射光的光通量数据；步骤S1024，根据上述光通量数据确定吸光度数据。
获取吸光度数据后，可执行步骤S104。
步骤S104，根据所述吸光度数据应用两点法计算第一吸光度差，根据吸光度数据应用速率法计算第二吸光度差；
在本申请的一些可选的实施例中，图2为本申请实施例提供的两点法计算第一吸光度差的流程图；根据吸光度数据应用两点法计算第一吸光度差可通过以下步骤S202-S206实现：
步骤S202，获取光通量曲线；该光通量曲线即根据所述光通量数据生成的曲线；
在本申请的一些可选的实施例中，可对光通量数据进行采样，并根据采样结果绘制光通量曲线，其中，光通量为光源发出的光束照射样本之后的散射光和/或透射光的光通量；
光通量数据与时间信息对应；例如,光通量数据可以用如下函数表示：
Intensity＝f(time)，其中，上述Intensity表示光源发出的光束照射样本之后的散射光和/或透射光的光通量数据，f(time)表示光通量数据与时间信息对应；可基于光通量数据与时间信息的对应关系，绘制出光通量曲线。
步骤S204，根据光通量曲线确定吸光度曲线；吸光度曲线即根据所述吸光度数据生成的曲线；
在本申请的一些可选的实施例中，吸光度数据与光通量数据具有以下关系:
Abs＝log10(I0/Intensity)；
其中，Abs表示吸光度数据，I0表示光源发出的光束照射样本前的入射光通量数据、参考光通量数据或者其他作为参考的光通量数据，基于光源发出的光束照射样本之后的散射光和/或透射光的光通量数据Intensity以及光源发出的光束照射样本前的入射光通量数据I0确定吸光度数据或吸光度曲线。
步骤S206，基于吸光度曲线确定第一吸光度差。
在本申请的一些可选的实施例中，吸光度差可通过以下公式表示：dOD Twopoint＝Abs(EndPoint)-Abs(StartPoint)；
其中，dOD Twopoint表示第一吸光度差，Abs(EndPoint)表示终止时间点对应的吸光度数据，Abs(StartPoint)表示初始时间点的吸光度数据。
在本申请的一些可选的实施例中，图3为本申请实施例提供的速率法计算第二吸光度差的流程图；根据吸光度数据应用速率法计算第二吸光度差可通过以下步骤S302-S310实现：
步骤S302，采集随反应时间变化的光通量曲线；
步骤S304，寻找反应速率最快的时刻，记作MaxPoint；
步骤S306，在MaxPoint附近寻找近似满足线性反应过程的区间，并将区间的起始点和终止点记作Cut_Min和Cut_Max；
步骤S308，计算该区间的起始点和终止点的吸光度差，分别记作：Abs(Cut_Min)、Abs(Cut_Max)；
步骤S310，计算得到吸光度差：dOD RateMethod＝Abs(Cut_Max)-Abs(Cut_Min)。
在本申请的一些可选的实施例中，图4为本申请实施例提供的另一种速率法计算 第二吸光度差的流程图；根据吸光度数据应用速率法计算第二吸光度差可通过以下步骤S402-S416实现：
步骤S402，仪器采集到光通量曲线：Intensity＝f(time)；其中，上述Intensity表示光源发出的光束照射样本之后的散射光和/或透射光的光通量数据，f(time)表示光通量数据与时间信息对应；可基于光通量数据与时间信息的对应关系，绘制出光通量曲线；
步骤S404，计算获得吸光度曲线：Abs＝log10(I0/Intensity)；Abs表示吸光度数据，I0表示光源发出的光束照射样本前的入射光通量数据、参考光通量数据或者其他作为参考的光通量数据，基于光源发出的光束照射样本之后的散射光和/或透射光的光通量数据Intensity以及光源发出的光束照射样本前的入射光通量数据I0确定吸光度数据或吸光度曲线
步骤S406，在分析区间内，N阶多项式拟合吸光度曲线(N≥2)，获得Abs_fit；
步骤S408，对Abs_fit曲线求导，获得切线方程，求得速率最大的时间点MaxPoint；
步骤S410，从MaxPoint开始、以最短回归时间为单位扩展线性范围；
步骤S412，求得原始吸光度曲线与最快点切线之间的积分面积，求得小于设定阈值的最大线性范围；
步骤S414，获得切线方程Abs_cut和线性范围的极限点：Cut_Min和Cut_Max；
步骤S416，计算得到吸光度差：dOD RateMethod＝Abs_cut(Cut_Max)-Abs_cut(Cut_Min)；
计算得到第一吸光度差和第二吸光度差后，可执行步骤S106。
步骤S106，将第一吸光度差和第二吸光度差代入用于计算最终吸光度差的预设的加权平均函数(也可称权重分配函数或比例分配函数)；
在本申请的一些可选的实施例中，上述加权平均函数为F(x)＝a·x+b·(c-x)+d；其中，a为第二吸光度差，x为第二吸光度差的权重，b为第一吸光度差，c-x为第一吸光度差的权重，c为大于等于x的常数(通常可以取1)，d为大于等于0的常数(通常可以取0)，F(x)为最终吸光度差，·表示乘号。当c＝1且d＝0时，加权平均函数为F(x)＝a·x+b·(1-x)，即第一吸光度差的权重和第二吸光度差的权重之和为1，从而确定了任意一个权重，即可确定另一个权重。在以下说明中，主要以加权平均函数为F(x)＝a·x+b·(1-x)进行说明。
步骤S108，根据预设函数关系确定满足所述加权平均函数的第一吸光度差的权重和第二吸光度差的权重以及最终吸光度差；其中，预设函数关系为针对第一吸光度差的权重、第二吸光度差的权重和吸光度差的函数关系。
可以在数据库中预存“第一吸光度差的权重-第二吸光度差的权重-吸光度差”的函数关系，例如预存迭代方程或大量的“第一吸光度差的权重-第二吸光度差的权重-吸光度差”的数据组合，例如“50％-50％-4000”“98％-2％-2000”“2％-98％-6000”等等数据组合。然后，在已经代入第一吸光度差和第二吸光度差的加权平均函数中，在数据库中迭代寻找满足加权平均函数的数据组合，从而确定第一吸光度差的权重和第二吸光度差的权重以及最终吸光度差。例如，在计算得到第二吸光度差a和第一吸光度差b后，代入加权平均函数为F(x)＝a·x+b·(1-x)，然后再从在数据库中迭代寻找满足加权平均函数的数据组合。
在一个可选实施例中，步骤S108的迭代寻找数据组合的过程可以表现为以下处理过程，但不限于此：依据第一吸光度差的权重和第二吸光度差的权重确定目标吸光度差；判断目标吸光度差、第一吸光度差的权重和第二吸光度差的权重是否满足预设加权平均函数；在判断结果为是时，确定目标吸光度差为最终吸光度差。
在本申请的一些可选的实施例中，上述预设函数关系对应的预设函数主要为了实现：吸光度差值较低的，对应的最终吸光度差的结果无限趋近于第一吸光度差；吸光度差值较高的，对应的最终吸光度差的结果无限趋近于第二吸光度差；吸光度差值处于中间的，过渡区间内尽量平稳转换、无跳变、无重叠、无断层。
在本申请的一些可选的实施例中，所述预设函数关系还需考虑：在设定的定义域内，如果吸光度差小于设定阈值，则所述第一吸光度差的权重大于所述第二吸光度差的权重；如果吸光度差大于所述设定阈值，则所述第一吸光度差的权重小于所述第二吸光度差的权重。在一些实施例中，设定阈值可以由用户自行设置。
由于确定了任意一个权重，即可确定另一个权重，因此所述预设函数关系中的预设函数可以是一个关于其中一个权重的函数，例如关于第二吸光度差的权重的函数。
预设函数可以是Prop＝{tanh[(dOD-dOD0)*k+1]}*0.5；其中，tanh表示双曲正切函数，双曲正切函数的基础函数为：tanh(x)＝[exp(-x)-exp(x)]/[exp(-x)+exp(x)]；其中，Prop为第二吸光度差的权重，dOD表示吸光度差的值，dOD0为上述的预设阈值、k为衰减时间系数，*表示乘号；其中，dOD0与k可以由用户设定。当dOD＝dOD0时，第一吸光度差的权重和第二吸光度差的权重都为0.5。请参阅图5，dOD0为5000，当dOD＝5000时，第一吸光度差的权重和第二吸光度差的权重都为0.5。因此在一些实施 例中，在所述设定的定义域内，如果吸光度差等于所述设定阈值时，则所述第一吸光度差的权重和所述第二吸光度差的权重是相同的。
在本申请的一些可选的实施例中，如果吸光度差小于所述设定阈值，所述第一吸光度差的权重与所述第二吸光度差的权重相差值的绝对值与吸光度差负相关；如果吸光度差大于所述设定阈值，所述权重相差值的绝对值与吸光度差正相关。继续以上述预设函数为例，请参阅图5，图5为本申请实施例提供的所述预设函数基于双曲正切函数时，确定的两点法与速率法的权重变化的示意图；横轴为吸光度差，纵轴为权重(占比)，dOD0为5000，当吸光度差dOD小于5000时，随着dOD的减少第一吸光度差的权重逐渐变大而第二吸光度差的权重逐渐变小，因此第一吸光度差的权重和第二吸光度差的权重的权重相差值的绝对值逐渐变大，因此第一吸光度差的权重与第二吸光度差的权重相差值的绝对值与吸光度差负相关。当dOD大于5000时，随着dOD的增大第一吸光度差的权重逐渐变小而第二吸光度差的权重逐渐变大，因此第一吸光度差的权重和第二吸光度差的权重的权重相差值的绝对值逐渐变大，因此第一吸光度差的权重与第二吸光度差的权重相差值的绝对值与吸光度差正相关。
在本申请的一些可选的实施例中，如果吸光度差小于所述设定阈值，所述权重相差值的绝对值随着吸光度差的增大而减小，所述权重相差值的绝对值的减小速度和所述吸光度差正相关。在本申请的一些可选的实施例中，如果吸光度差大于所述设定阈值，所述权重相差值的绝对值随着吸光度差的增大而增大，所述权重相差值的绝对值的增大速度和所述吸光度差负相关。继续以上述预设函数为例，请参阅图5，dOD0为5000，当吸光度差dOD小于5000时，随着dOD的减少第一吸光度差的权重逐渐变大但是变化速度在下降，而第二吸光度差的权重逐渐变小但是变化速度也在下降，因此随着dOD的变少，第一吸光度差的权重和第二吸光度差的权重的权重相差值的绝对值逐渐变大但是变化速度在下降，即权重相差值的绝对值的减小速度和所述吸光度差正相关。当吸光度差dOD大于5000时，随着dOD的增大第一吸光度差的权重逐渐变小但是变化速度在下降，而第二吸光度差的权重逐渐变大但是变化速度也在下降，因此随着dOD的变大，第一吸光度差的权重和第二吸光度差的权重的权重相差值的绝对值逐渐变大但是变化速度在下降，即权重相差值的绝对值的增大速度和吸光度差负相关。在本申请的一些可选的实施例中，权重相差值的绝对值的减小速度或增大速度由用户设定。例如上述预设函数中的k可由用户进行设置。
当然在其他实施例中，预设函数还可以是一阶***，见图6，例如y＝1-exp(-x/τ)，其中横轴为吸光度差，纵轴为权重(占比)，τ为响应时间，x为吸光度差dOD，y为第二吸光度差的权重。
图7为由所述预设函数确定的第一吸光度差与第二吸光度差的权重的变化情况的示意图。图7中，横轴为吸光度差，纵轴为两点法对应的第一吸光度差与速率法对应的第二吸光度差权重的占比。在吸光度差的低值部分，所述预设函数能够使得总吸光度差无限接近于两点法计算结果，充分利用两点法的线性拟合处理抵抗噪声干扰、提升信噪比，以获得更好的低值重复性；在吸光度差的高值部分，所述预设函数能够使得总吸光度差无限接近于速率法计算结果，充分利用速率法积极反映反应速率随时间的变化关系，以获得吸光度差计算值的单调性和更高的线性测量范围。
例如：当设定阈值为4000时，在dOD≤1000位置(低值所属区间)，两点法计算值占总吸光度差的比重为≥99.75％，必然决定总吸光度差的值；在dOD≥7000位置(中值附近，尚不属于高值部分)，两点法计算值占总吸光度差的比重≤0.25％，总吸光度差就与速率法计算结果保持一致；两点法和速率法在总吸光度差中的分配比例是连续调整的，导数连续、不存在突变，因此整个过程也是连续的。
在本申请的一些可选的实施例中，提供了如图8所示的样本吸光度差的测量方法的流程示意图；该方法包括以下步骤S802、步骤S8022-S8024、步骤S8042-步骤S8048以及步骤S806，其中步骤S8022-S8024用于计算第一吸光度差，步骤S8042-步骤S8048用于计算第二吸光度差。
步骤S802，采集随反应时间变化的光通量曲线；
S8022-S8024具体如下：
步骤S8022，计算采样起始点和采样终止点的吸光度，分别记作：Abs(StartPoint)、Abs(EndPoint)；
步骤S8024，计算得到第一吸光度差：dOD Twopointn＝Abs(EndPoint)-Abs(StartPoint)；
步骤S8042-步骤S8048具体如下：
步骤S8042，寻找反应速率最快的时刻，记作MaxPoint；
步骤S8044，在MaxPoint附近寻找近似满足线性反应过程的区间，并将区间的起始点和终止点记作Cut_Min和Cut_Max；
步骤S8046，计算该区间的起始点和终止点的吸光度差，分别记作：Abs(Cut_Min)、Abs(Cut_Max)；
步骤S8048，计算得到第二吸光度差：dOD RateMethod＝Abs(Cut_Max)-Abs(Cut_Min)；
通过上述步骤得到第一吸光度差dOD Twopoint和第二吸光度差dOD RateMethod后，执行步骤S806，即通过加权平均函数进行分配权重计算最终的吸光度差；
步骤S806，按照一定的配比获得最终的吸光度差dOD＝dOD Twopoint*prop(dOD)+dOD RateMethod*(1-prop(dOD))；
其中，dOD Twopoint为两点法计算的吸光度差(第一吸光度差)，dOD RateMethod为速率法计算的吸光度差(第二吸光度差)，prop(dOD)为两点法计算的吸光度差所占的权重，1-prop(dOD)为速率法计算的吸光度差所占的权重。
在本申请的一些可选的实施例中，提供了如图9所示的样本吸光度差的测量方法的流程示意图；该方法包括以下步骤S902-步骤S904、步骤S9062-步骤S9064、步骤S9066-步骤S90616和步骤S908，其中步骤S9062-步骤S9064用于计算第一吸光度差，S9066-步骤S90616用于计算第二吸光度差。
步骤S902-步骤S904具体如下：
步骤S902，仪器采集到光通量曲线：Intensity＝f(time)；
步骤S904，计算获得吸光度曲线：Abs＝log10(I0/Intensity)；
步骤S9062-步骤S9064具体如下：
步骤S9062，在分析区间内，线性拟合吸光度曲线：Abs fit；
步骤S9064，以两点法计算得到第一吸光度差：dOD Twopoint＝：Abs fit(EndPoint)-：Abs fit(StartPoint)；
步骤S9066-步骤S90616具体如下：
步骤S9066，在分析区间内，N阶多项式拟合吸光度曲线(N≥2)，获得：Abs fit；
步骤S9068，对：Abs fit曲线求导，获得切线方程，求得速率最大点MaxPoint；
步骤S90610，从MaxPoint开始、以最短回归时间为单位扩展线性范围；
步骤S90612，求得原始吸光度曲线与最快点切线之间的积分面积，求得小于设定阈值的最大线性范围；
步骤S90614，获得切线方程Abs cut和线性范围的极限点：Cut Min和Cut Max；
步骤S90616，以速率法计算得到第二吸光度差：dOD RateMethod＝Abs cut(Cut Max)-Abs cut(Cut Min)；
通过上述步骤得到第一吸光度差dOD Twopoint和第二吸光度差dOD RateMethod后，执行步骤S908，即通过加权平均函数进行分配权重计算最终的吸光度差。
步骤S908，按照一定的配比获得最终的吸光度差dOD＝dOD Twopoint*prop(dOD)+dOD RateMethod*(1-prop(dOD))。
其中，dOD Twopoint为两点法计算的吸光度差(第一吸光度差)，dOD RateMethod为速率法计算的吸光度差(第二吸光度差)，prop(dOD)为两点法计算的吸光度差所占的权重，1-prop(dOD)为速率法计算的吸光度差所占的权重。
本申请实施例采用的算法同时兼得两点法和速率法的优点。通过对4个样本进行实验，采用两点法和速率法计算吸光度差低值部分的重复性计算结果如图10中的图10a,图10b,图10c,图10d；可见电路噪声干扰导致吸光度曲线出现波动，两点法可以通过拟合处理获得很好的重复性，速率法则会受到强烈的干扰，导致重复性不良。
如采用本申请实施例的算法，低值重复性与采用两点法计算得到的计算结果基本一致。表1示出了对10个样本进行实验得到的吸光度差的计算结果，可见采用本申请实施例的算法与两点法的计算结果基本一致。
通过对多个不同浓度的样本进行实验计算吸光度差，高值部分的线性计算结果如图11所示，图11为吸光度差为高值时，采用两点法、速率法、与使用本申请所述的样本吸光度差的测量方法计算的吸光度差的对比示意图。图11中DD浓度表示DD(D-二聚体)的浓度；采用两点法计算的吸光度差有饱和、甚至降低的趋势，速率法计算的吸光度差则保持单调递增、有效的提升了线性测量范围，本申请实施例采用的算法的线性测量范围同速率法基本一致。
表1低值重复性的计算结果
Figure PCTCN2018116778-appb-000001
Figure PCTCN2018116778-appb-000002
在一些实施例中，样本分析仪可以提供三种模式，一种模式采用两点法进行计算，一种模式采用速率法进行计算，一种模式采用本申请实施例样本吸光度差的测量方法进行计算，供用户自行选择。
在本申请的一些可选的实施例中，本申请还提供一种样本分析仪，对应的工作方法流程图如图12所示。该工作方法包括以下步骤：
步骤S1202，用户选择是否使用本申请实施例算法，若是，则执行步骤S1204，若否，则执行步骤S1210；
步骤S1204，仪器获取配置参数；
步骤S1206，仪器开展实验，并同时使用两点法和速率法计算吸光度差；
步骤S1208，仪器按照本申请实施例的方法计算吸光度差；
步骤S1210，用户选择其他算法；这里的其他算法可以包括两点法或速率法等等。
步骤S1212，仪器开展实验，按照用户选定的方法计算吸光度差；
步骤S1214，用户获取检测结果。
需要用户选择是否使用本申请实施例的算法，如果不使用，则仪器按照用户选定的计算方法(例如两点法或速率法)开展实验并上报结果。如果使用，则需要进一步获取配置参数。然后，仪器开展实验，并且同时使用两点法和速率法计算吸光度差，按照本申请实施例提供的算法计算最终的吸光度差。最后，仪器上报最终的吸光度差，并根据吸光度差得到检测结果，从而用户获取检测结果，如样本浓度。这里的样本浓度是为了简洁而定义的简称，可以理解为样本中特定物质的浓度，例如可以是D-二聚体(DD)、纤维蛋白(原)降解产物(FDP)等浓度。
在本申请的一些可选的实施例中，根据预设函数关系确定满足所述加权平均函数 的第一吸光度差的权重和第二吸光度差的权重以及最终吸光度差之后，还可以执行以下步骤：依据预设对应关系确定与所述最终吸光度差对应的样本浓度，其中，所述预设对应关系为吸光度差与样本浓度的对应关系。样本浓度是与吸光度差是具有对应关系的，这个对应关系预存在数据库中。确定了吸光度差，即可根据这个对应关系确定样本浓度，由于属于现有技术，在此不再赘述。
在本申请的一些可选的实施例中，上述方法应用于样本分析仪，例如血液分析仪、生化分析仪、免疫分析仪、凝血分析仪等等体外诊断设备。
在本申请实施例中，通过获取样本对象的吸光度数据；根据吸光度数据应用两点法计算第一吸光度差，根据吸光度数据应用速率法计算第二吸光度差；将第一吸光度差和第二吸光度差代入用于计算最终吸光度差的预设的加权平均函数；根据预设函数关系确定满足加权平均函数的第一吸光度差的权重和第二吸光度差的权重以及最终吸光度差，其中，预设函数关系为针对第一吸光度差的权重、第二吸光度差的权重和吸光度差的函数关系。达到了通过低值部分采用两点法，以获得更好的低值重复性；高值部分采用速率法，以获得吸光度差计算值的单调性和更高的线性测量范围；低值和高值部分良好连接，不存在断层或跳变的技术效果；从而解决了在测量吸光度差时，无法同时兼顾吸光度差低值时的测量重复性以及吸光度差高值时的线性测量范围的技术问题。上述样本分析仪的设定阈值可以由用户设定。另外，上述权重相差值的绝对值的减小速度或增大速度也可以由用户设定。
本申请实施例还提供了一个或多个非易失性计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行上述任一实施例样本吸光度差的测量方法中的步骤。
对应上述的样本吸光度差的测量方法，本申请实施例提供了图1a所示的样本分析仪，如图1a所示，该样本分析仪包括：光源132、检测器134、处理器136，其中，光源132用于发出用于照射样本的光束；
检测器134用于检测在光束照射样本后产生的光通量数据；
在本申请的一些可选的实施例中，上述处理器136还用于执行：依据预设对应关系确定与最终吸光度差对应的样本浓度，其中，预设对应关系为吸光度差与样本浓度的对应关系。
在本申请的一些可选的实施例中，上述加权平均函数为F(x)＝a·x+b·(c-x)+d；其中，a为第二吸光度差，x为第二吸光度差的权重，b为第一吸光度差，1-x为第一吸光度差的权重，c为大于等于x的常数，d为大于等于0的常数，F(x)为最终吸光度 差。
在本申请的一些可选的实施例中，上述预设函数关系可以为：在设定的定义域内，如果吸光度差小于设定阈值，则第一吸光度差的权重大于第二吸光度差的权重；如果吸光度差大于设定阈值，则第一吸光度差的权重小于第二吸光度差的权重。
其中，预设函数关系参考参考上述对应图1b的测量样本吸光度方法的相关实施例。
在本申请的一些可选的实施例中，如果吸光度差小于设定阈值，第一吸光度差的权重与第二吸光度差的权重相差值的绝对值权重相差值的绝对值与吸光度差负相关；如果吸光度差大于设定阈值，权重相差值的绝对值与吸光度差正相关。
在本申请的一些可选的实施例中，如果吸光度差小于设定阈值，权重相差值的绝对值随着吸光度差的增大而减小，权重相差值的绝对值的减小速度和吸光度差正相关。如果吸光度差大于设定阈值，权重相差值的绝对值随着吸光度差的增大而增大，权重相差值的绝对值的增大速度和吸光度差负相关。
在本申请的一些可选的实施例中，设定阈值由用户设定。权重相差值的绝对值的减小速度或增大速度由用户设定。
上述实施例中的凝血分析仪，其用于光学测定并分析与血液凝固/纤溶功能有关的特定物质数量及其活性程度，标本为血浆。本实施方式的凝血分析仪用凝固时间法、发色底物法和免疫比浊法对标本进行光学测定。本实施方式所用的凝固时间法是一种将标本凝固过程作为透射光的变化检测的测定方法。测定项目有PT(凝血酶原时间)、APTT(活化部分凝血酶时间)、TT(凝血酶时间)和FIB(纤维蛋白原量)等。发色底物法的测定项目有AT-Ⅲ(抗凝血酶Ⅲ)等，免疫比浊法的测定项目有D-二聚体(D-Dimer)和FDP等。
该凝血分析仪包括：至少一个反应容器，用于为样本与试剂提供反应场所；样本量检测装置，用于对所述样本进行液量检测，得到样本的实测样本量；凝血检测装置，用于对经试剂处理后的样本进行凝血检测，得到反映凝固情况的电信号信息；处理器，用于接收并处理所述凝血检测装置输出的电信号信息，以得到样本的测量参数；其中，所述处理器还用于执行根据上述任一项实施例所述的样本吸光度差的测量方法。
上述本申请实施例序号仅仅为了描述，不代表实施例的优劣。
在本申请的上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。
在本申请所提供的几个实施例中，应该理解到，所揭露的技术内容，可通过其它 的方式实现。其中，以上所描述的装置实施例仅仅是示意性的，例如单元的划分，可以为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个***，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，单元或模块的间接耦合或通信连接，可以是电性或其它的形式。
作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本申请各个实施例方法的全部或部分步骤。而前述的存储介质包括：U盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
以上仅是本申请的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本申请原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本申请的保护范围。

Claims (21)

  1. 一种样本吸光度差的测量方法，其特征在于，包括；
    获取样本对象的吸光度数据；
    根据所述吸光度数据应用两点法计算第一吸光度差，根据吸光度数据应用速率法计算第二吸光度差；
    将第一吸光度差和第二吸光度差代入用于计算最终吸光度差的预设的加权平均函数；
    根据预设函数关系确定满足所述加权平均函数的第一吸光度差的权重和第二吸光度差的权重以及最终吸光度差，其中，所述预设函数关系为针对所述第一吸光度差的权重、第二吸光度差的权重和吸光度差的函数关系。
  2. 根据权利要求1所述的方法，其特征在于，根据预设函数关系确定满足所述加权平均函数的第一吸光度差的权重和第二吸光度差的权重以及最终吸光度差之后，所述方法还包括：
    依据预设对应关系确定与所述最终吸光度差对应的样本浓度，其中，所述预设对应关系为吸光度差与样本浓度的对应关系。
  3. 根据权利要求1所述的方法，其特征在于，所述预设函数关系包括：
    在设定的定义域内，如果吸光度差小于设定阈值，则所述第一吸光度差的权重大于所述第二吸光度差的权重；如果吸光度差大于所述设定阈值，则所述第一吸光度差的权重小于所述第二吸光度差的权重。
  4. 根据权利要求3所述的方法，其特征在于，如果吸光度差小于所述设定阈值，所述第一吸光度差的权重和所述第二吸光度差的权重之间的权重相差值的绝对值与吸光度差负相关；如果吸光度差大于所述设定阈值，所述权重相差值的绝对值与吸光度差正相关。
  5. 根据权利要求4所述的方法，其特征在于，如果吸光度差小于所述设定阈值，所述权重相差值的绝对值随着吸光度差的增大而减小，所述权重相差值的绝对值的减小速度和所述吸光度差正相关。
  6. 根据权利要求4所述的方法，其特征在于，如果吸光度差大于所述设定阈值，所述权重相差值的绝对值随着吸光度差的增大而增大，所述权重相差值的绝对值的增 大速度和所述吸光度差负相关。
  7. 根据权利要求3至6中任一项所述的方法，其特征在于，所述设定阈值由用户设定。
  8. 根据权利要求5或6所述的方法，其特征在于，所述权重相差值的绝对值的变化速度由用户设定。
  9. 根据权利要求3所述的方法，其特征在于，所述预设函数关系还包括：
    在所述设定的定义域内，如果吸光度差等于所述设定阈值时，则所述第一吸光度差的权重和所述第二吸光度差的权重是相同的。
  10. 根据权利要求1所述的方法，其特征在于，所述加权平均函数为：
    F(x)＝a·x+b·(c-x)+d；
    其中，a为所述第二吸光度差，x为所述第二吸光度差的权重，b为所述第一吸光度差，(c-x)为所述第一吸光度差的权重，c为大于或等于x的常数，d为大于或等于0的常数，F(x)为所述最终吸光度差。
  11. 根据权利要求1所述的方法，其特征在于，所述方法应用于样本分析仪。
  12. 一种样本分析仪，其特征在于，包括：光源、检测器和处理器；其中，
    所述光源用于发出用于照射样本的光束；
    所述检测器用于检测在所述光束照射样本后产生的光通量数据；
    所述处理器运行程序，其中，所述程序运行时对于从检测器输出的数据执行如下处理步骤：
    根据所述光通量数据计算吸光度数据；
    根据所述吸光度数据应用两点法计算第一吸光度差，根据吸光度数据应用速率法计算第二吸光度差；
    将第一吸光度差和第二吸光度差代入用于计算最终吸光度差的预设的加权平均函数；
    根据预设函数关系确定满足所述加权平均函数的第一吸光度差的权重和第二吸光度差的权重以及最终吸光度差，其中，所述预设函数关系为针对所述第一吸光度差的权重、第二吸光度差的权重和吸光度差的函数关系。
  13. 根据权利要求12所述的样本分析仪，其特征在于，所述处理器还用于：依据预设 对应关系确定与所述最终吸光度差对应的样本浓度，其中，所述预设对应关系为吸光度差与样本浓度的对应关系。
  14. 根据权利要求12所述的样本分析仪，其特征在于，所述预设函数关系包括：
    在设定的定义域内，如果吸光度差小于设定阈值，则所述第一吸光度差的权重大于所述第二吸光度差的权重；如果吸光度差大于所述设定阈值，则所述第一吸光度差的权重小于所述第二吸光度差的权重。
  15. 根据权利要求14所述的样本分析仪，其特征在于，如果吸光度差小于所述设定阈值，所述第一吸光度差的权重与所述第二吸光度差的权重相差值的绝对值权重相差值的绝对值与吸光度差负相关；如果吸光度差大于所述设定阈值，所述权重相差值的绝对值与吸光度差正相关。
  16. 根据权利要求15所述的样本分析仪，其特征在于，如果吸光度差小于所述设定阈值，所述权重相差值的绝对值随着吸光度差的增大而减小，所述权重相差值的绝对值的减小速度和所述吸光度差正相关。
  17. 根据权利要求15所述的样本分析仪，其特征在于，如果吸光度差大于所述设定阈值，所述权重相差值的绝对值随着吸光度差的增大而增大，所述权重相差值的绝对值的增大速度和所述吸光度差负相关。
  18. 根据权利要求14至17中任一所述的样本分析仪，其特征在于，所述设定阈值由用户设定。
  19. 根据权利要求16或17所述的样本分析仪，其特征在于，所述权重相差值的绝对值的变化速度由用户设定。
  20. 根据权利要求12所述的样本分析仪，其特征在于，所述加权平均函数为：
    F(x)＝a·x+b·(c-x)+d；
    其中，a为第二吸光度差，x为第二吸光度差的权重，b为第一吸光度差，c-x为第一吸光度差的权重，c为大于或等于x的常数，d为大于或等于0的常数，F(x)为最终吸光度差。
  21. 一个或多个非易失性计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1-11任一项所述的方法中的步骤。
PCT/CN2018/116778 2018-11-21 2018-11-21 样本吸光度差的测量方法、样本分析仪和存储介质 WO2020103051A1 (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201880099067.1A CN112912712B (zh) 2018-11-21 2018-11-21 样本吸光度差的测量方法、样本分析仪和存储介质
PCT/CN2018/116778 WO2020103051A1 (zh) 2018-11-21 2018-11-21 样本吸光度差的测量方法、样本分析仪和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2018/116778 WO2020103051A1 (zh) 2018-11-21 2018-11-21 样本吸光度差的测量方法、样本分析仪和存储介质

Publications (1)

Publication Number Publication Date
WO2020103051A1 true WO2020103051A1 (zh) 2020-05-28

Family

ID=70773768

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2018/116778 WO2020103051A1 (zh) 2018-11-21 2018-11-21 样本吸光度差的测量方法、样本分析仪和存储介质

Country Status (2)

Country Link
CN (1) CN112912712B (zh)
WO (1) WO2020103051A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113484244A (zh) * 2021-06-17 2021-10-08 江苏鸿恩医疗器械有限公司 一种基于免疫比浊法的凝血分析仪试剂定标方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001235422A (ja) * 2000-02-22 2001-08-31 Toshiba Corp 自動分析装置
CN101672856A (zh) * 2008-09-11 2010-03-17 深圳迈瑞生物医疗电子股份有限公司 一种半自动生化分析仪及其双波长测试方法
CN102243227A (zh) * 2010-12-16 2011-11-16 浙江亚太药业股份有限公司 非对称二甲基精氨酸浓度的测定方法及测定试剂
WO2017012168A1 (zh) * 2015-07-22 2017-01-26 江苏英诺华医疗技术有限公司 一种优化的适合医学检验的生化检测方法
US9745623B2 (en) * 2011-10-14 2017-08-29 Becton Dickinson And Company Square wave thermal cycling

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1059522A1 (en) * 1999-06-11 2000-12-13 F. Hoffmann-La Roche Ag Method and apparatus for examining fluids of biological origin
WO2006104006A1 (ja) * 2005-03-29 2006-10-05 Sysmex Corporation 検体分析方法および検体分析装置
CN101169439A (zh) * 2006-10-24 2008-04-30 苏州艾杰生物科技有限公司 二氧化碳诊断/测定试剂盒及二氧化碳浓度的测定方法
JP5319696B2 (ja) * 2008-11-17 2013-10-16 株式会社日立ハイテクノロジーズ 自動分析装置
CN105572041B (zh) * 2014-10-16 2020-09-08 北京普利生仪器有限公司 光学检测设备、方法及其光学检测装置和反应杯
CN106290940B (zh) * 2016-08-31 2018-06-29 江苏英诺华医疗技术有限公司 应用于生化分析仪的速率法检测数据智能处理方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001235422A (ja) * 2000-02-22 2001-08-31 Toshiba Corp 自動分析装置
CN101672856A (zh) * 2008-09-11 2010-03-17 深圳迈瑞生物医疗电子股份有限公司 一种半自动生化分析仪及其双波长测试方法
CN102243227A (zh) * 2010-12-16 2011-11-16 浙江亚太药业股份有限公司 非对称二甲基精氨酸浓度的测定方法及测定试剂
US9745623B2 (en) * 2011-10-14 2017-08-29 Becton Dickinson And Company Square wave thermal cycling
WO2017012168A1 (zh) * 2015-07-22 2017-01-26 江苏英诺华医疗技术有限公司 一种优化的适合医学检验的生化检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
MA, FEI ET AL.: "Analysis Method of Biochemistry Analyzer", CHINESE MEDICAL EQUIPMENT JOURNAL, 31 December 2003 (2003-12-31), pages 138, ISSN: 1003-8868 *

Also Published As

Publication number Publication date
CN112912712B (zh) 2023-11-21
CN112912712A (zh) 2021-06-04

Similar Documents

Publication Publication Date Title
JP2500848B2 (ja) 血液凝固時間の決定方法
US8965774B2 (en) Automatic detection of audio compression parameters
CN110082319B (zh) 定标数据修正方法及其电子设备
US7096715B2 (en) Odor discriminating apparatus
JP2012505406A5 (zh)
JP2938302B2 (ja) 血液凝固時間測定方法とその装置
WO2015132861A1 (ja) クロマトグラフ質量分析装置用データ処理装置及びプログラム
CN112345759A (zh) 一种用于荧光强度峰检测的方法
KR102539142B1 (ko) 스펙트럼 분석 장치와 방법, 및 혈당 측정 장치
WO2020103051A1 (zh) 样本吸光度差的测量方法、样本分析仪和存储介质
US20210063362A1 (en) Techniques for exception-based validation of analytical information
CN108761105B (zh) 凝血时间判定方法、电子设备及存储介质
JPH0627115A (ja) 血液凝固時間測定方法とその装置
EP2265965A1 (en) Relative noise of a measured signal
CN110596020A (zh) 物质检测方法及装置、存储介质、通信终端
US20220178865A1 (en) Measuring method for semiconductor gas sensor based on alternating-current impedance
Defrance et al. Finding the onset of a room impulse response: Straightforward?
JP7423831B2 (ja) 流体試験システム
CN113445099A (zh) 槽液分析方法、装置及计算机可读存储介质
US10236167B1 (en) Peak waveform processing device
CN114548189B (zh) 一种凝血酶原时间的检测方法、检测装置及计算机可读存储介质
CN110927068A (zh) 基于非共振光声光谱原理的气体浓度测量装置及方法
Giavarina et al. Method comparison of automated systems for the erythrocyte sedimentation rate
JPH0954072A (ja) クロマトグラフのデータ処理装置
WO2024046267A1 (zh) 基于反应曲线的曲率检测待检测样本的前带现象的方法和装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18941004

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18941004

Country of ref document: EP

Kind code of ref document: A1