CN115048613B - 一种指标同质化换算方法、装置、电子设备及存储介质 - Google Patents
一种指标同质化换算方法、装置、电子设备及存储介质 Download PDFInfo
- Publication number
- CN115048613B CN115048613B CN202210978315.8A CN202210978315A CN115048613B CN 115048613 B CN115048613 B CN 115048613B CN 202210978315 A CN202210978315 A CN 202210978315A CN 115048613 B CN115048613 B CN 115048613B
- Authority
- CN
- China
- Prior art keywords
- standard curve
- end point
- concentration
- arc segment
- target
- 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
- 238000000265 homogenisation Methods 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 24
- 238000012360 testing method Methods 0.000 claims abstract description 88
- 238000005259 measurement Methods 0.000 claims abstract description 13
- 230000010354 integration Effects 0.000 claims description 32
- 230000006870 function Effects 0.000 claims description 18
- 239000002131 composite material Substances 0.000 claims description 16
- 150000001875 compounds Chemical class 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims 3
- 238000001514 detection method Methods 0.000 abstract description 52
- 238000012545 processing Methods 0.000 description 16
- 230000008569 process Effects 0.000 description 7
- 230000009471 action Effects 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000013329 compounding Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000002835 absorbance Methods 0.000 description 1
- 238000011088 calibration curve Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004393 prognosis Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Indication And Recording Devices For Special Purposes And Tariff Metering Devices (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
本申请提供一种指标同质化换算方法、装置、电子设备及存储介质，该方法包括：获取实测信号值；确定实测信号值在第一标准曲线的浓度；根据第一标准曲线的第一端点和第一标准曲线的第二端点，获得第一标准曲线的第一目标弧线段，根据第一标准曲线的第一端点和第一标准曲线测试样本对应的端点，获得第一标准曲线的第二目标弧线段；根据第二标准曲线的第一端点和第二标准曲线的第二端点，获得第二标准曲线的第三目标弧线段；根据第一目标弧线段、第二目标弧线段和第三目标弧线段，确定第四目标弧线段的弧长方程；根据第四目标弧线段的弧长方程，确定实测信号值对应的同质化浓度。该方法用以实现不同检测设备的检测结果的同质化。
Description
技术领域
本申请涉及检测技术领域，具体而言，涉及一种指标同质化换算方法、指标同质化换算装置、电子设备及计算机可读存储介质。
背景技术
在医疗检测场景中，一种场景下：同一实验室采用两台设备对同一样本进行同一指标检测，获得的检测结果差别较大；另一种场景下：不同的实验室采用各自的设备对同一样本进行同一指标检测，检测结果差别也很大。
对于两台检测设备来说，在上述两种情况下，都存在不同设备检测同一样本，检测结果存在差异的问题。从临床医生的角度而言，希望同一个样本在各个设备上进行同一指标的检测，获得的检测结果保持一致，即检测结果在较小的范围内波动，而这种波动不影响临床医生对疾病诊断和预后判断。
然而，目前市面上缺少实现不同设备间对同一样本的指标检测换算方法，即缺少使两台检测设备针对同一样本的同一检测指标的检测结果趋于一致的方法。
发明内容
有鉴于此，本申请实施例的目的在于提供一种指标同质化换算方法、指标同质化换算装置、电子设备、存储介质，用以实现不同检测设备的检测结果的同质化。
第一方面，本申请实施例提供一种指标同质化换算方法，包括：获取实测信号值，所述实测信号值为测试设备对测试样本进行检测获得的信号值；确定所述实测信号值在第一标准曲线的对应浓度；所述第一标准曲线为所述测试设备预设的浓度-信号值关系曲线；根据所述第一标准曲线的第一端点和所述第一标准曲线的第二端点，获得所述第一标准曲线的第一目标弧线段，根据所述第一标准曲线的第一端点和所述第一标准曲线所述测试样本对应的端点，获得所述第一标准曲线的第二目标弧线段，所述第一标准曲线的第一端点和第二端点为不同比对样本对应的浓度-信号值；根据第二标准曲线的第一端点和所述第二标准曲线的第二端点，获得所述第二标准曲线的第三目标弧线段，所述第二标准曲线为参比设备预设的浓度-信号值关系曲线，所述第二标准曲线的第一端点与所述第一标准曲线的第一端点为同一比对样本对应的浓度-信号值，所述第二标准曲线的第二端点与所述第一标准曲线的第二端点为同一比对样本对应的浓度-信号值；根据所述第一目标弧线段、所述第二目标弧线段和所述第三目标弧线段，确定第四目标弧线段的弧长方程，所述第四目标弧线段为所述第二标准曲线的第一端点到所述测试样本在所述第二标准曲线对应的端点；根据所述第四目标弧线段的弧长方程，确定所述实测信号值对应的同质化浓度。
在本申请实施例中，不同的检测设备分别为测试设备和参比设备，通过检测设备对测试样本进行检测，获得测试样本对应的信号值和浓度，从而确定测试样本在第一标准曲线的坐标，再根据第一标准曲线的第一端点和第二端点，获得第一标准曲线上的第一端点和第二端点连接的第一目标弧线段，根据第一标准曲线的第一端点与测试样本的坐标，获得第一标准曲线上的第一端点和测试样本坐标的第二目标弧线段，同理，获得在第二标准曲线上的第三目标弧线段。进一步地，由于第一标准曲线上的第一端点和第二标准曲线上的第一端点都是通过同一对比样本获得的不同检测结果，两点具有对应关系，同理，第一标准曲线上的第二端点与第二标准曲线上的第二端点对应。
在获得第一目标弧线段、第二目标弧线段、第三目标弧线段后，根据类似原理，建立第四目标弧线段的弧长方程，其中，第四目标弧线段为第二标准曲线的第一端点到测试样本在第二标准曲线对应的端点的弧线。再对第四目标弧线段的弧长方程进行求解，即可确定测试样本在第二标准曲线上的浓度，即，确定实测信号值对应的同质化浓度，最终实现不同检测设备的检测结果的同质化。
一实施例中，所述根据所述第一目标弧线段、所述第二目标弧线段和所述第三目标弧线段，确定第四目标弧线段的弧长方程，包括：根据所述第一目标弧线段、所述第二目标弧线段、第三目标弧线段各自的弧长，确定出所述第四目标弧线段的弧长与各目标弧线段的弧长的长度关系；基于弧长计算公式和所述长度关系，确定所述第四目标弧线段的弧长方程。
一实施例中，所述根据所述第四目标弧线段的弧长方程，确定所述实测信号值对应的同质化浓度，包括：通过复合辛普森求积分公式和二分法，对所述弧长方程进行处理，确定所述实测信号值对应的同质化浓度。
在本申请实施例中，通过复合辛普森求积分公式和二分法对弧长方程进行处理，可以减小计算量，进而提高数据处理的效率。
一实施例中，所述通过复合辛普森求积分公式和二分法，对所述弧长方程进行处理，确定所述实测信号值对应的同质化浓度，包括：通过复合辛普森求积分公式对所述弧长方程进行转换，获得目标积分方程；通过二分法对所述目标积分方程进行求解，确定所述实测信号值对应的同质化浓度。
一实施例中，在所述基于弧长计算公式和所述长度关系，确定所述第四目标弧线段的弧长方程之前，所述方法还包括：根据所述第一标准曲线的第一端点、所述第一标准曲线的第二端点、所述第一标准曲线中所述测试样本对应的端点、所述第二标准曲线的第一端点和所述第二标准曲线的第二端点，确定所述复合辛普森求积分公式的步长值。
在本申请实施例中，通过第一标准曲线的第一端点、第一标准曲线的第二端点、第一标准曲线中测试样本对应的端点、第二标准曲线的第一端点和第二标准曲线的第二端点，对步长值进行确定，以减少最终获取同质化浓度的误差。
一实施例中，在所述根据所述第四目标弧线段的弧长方程，确定所述实测信号值对应的同质化浓度之后，所述方法还包括：判断所述同质化浓度与所述第二标准曲线的第一端点或第二端点的浓度之差是否小于预设阈值，若是，则根据公式：或，对所述同质化浓度进行校准；其中，为校准后的同质化浓度，为所述第二目标弧线段的第一端点的浓度，为所述第二目标弧线段的第二端点的浓度，为所述同质化浓度。
在本申请实施例中，通过将同质化浓度与第二标准曲线的第一端点或第二端点的浓度之差进行判断，以在浓度之差小于预设阈值时，通过公式或对同质化浓度进行校准，使得最终的同质化浓度的数值可控，提高数据运算处理的效率，避免计算结果为无限小数或者迭代次数过多等情况出现，导致设备运算宕机。
第二方面，本申请实施例提供一种同质化换算装置，包括：获取模块，用于获取实测信号值，所述实测信号值为测试设备对测试样本进行检测获得的信号值；确定模块，用于确定所述实测信号值在第一标准曲线的浓度；所述第一标准曲线为所述测试设备预设的浓度-信号值关系曲线；所述获取模块，还用于根据所述第一标准曲线的第一端点和所述第一标准曲线的第二端点，获得所述第一标准曲线的第一目标弧线段，根据所述第一标准曲线的第一端点和所述第一标准曲线所述测试样本对应的浓度-信号值，获得所述第一标准曲线的第二目标弧线段，所述第一标准曲线的第一端点和第二端点为不同比对样本对应的端点；根据第二标准曲线的第一端点和所述第二标准曲线的第二端点，获得所述第二标准曲线的第三目标弧线段，所述第二标准曲线为参比设备预设的浓度-信号值关系曲线，所述第二标准曲线的第一端点与所述第一标准曲线的第一端点为同一比对样本对应的浓度-信号值，所述第二标准曲线的第二端点与所述第一标准曲线的第二端点为同一比对样本对应的浓度-信号值；所述确定模块，还用于根据所述第一目标弧线段、所述第二目标弧线段和所述第三目标弧线段，确定第四目标弧线段的弧长方程，所述第四目标弧线段为所述第二标准曲线的第一端点到所述测试样本在所述第二标准曲线对应的端点；根据所述第四目标弧线段的弧长方程，确定所述实测信号值对应的同质化浓度。
一实施例中，所述确定模块还用于根据所述第一目标弧线段、所述第二目标弧线段、第三目标弧线段各自的弧长，确定出所述第四目标弧线段的弧长与各目标弧线段的弧长的长度关系；基于弧长计算公式和所述长度关系，确定所述第四目标弧线段的弧长方程。
第三方面，本申请实施例提供一种电子设备，包括存储器和处理器，所述存储器中存储有计算机可读指令，所述计算机可读指令被所述处理器执行时，使得所述处理器执行上述第一方面的指标同质化换算方法或实现上述指标同质化换算装置的功能。
第四方面，本申请实施例提供一种存储有计算机可读指令的非易失性可读存储介质，所述计算机可读指令被处理器执行时，使得所述处理器执行第一方面的指标同质化换算方法或实现上述指标同质化换算装置的功能。
本申请的一个或多个实施例的细节在下面的附图和描述中提出。本申请的其它特征、目的和优点将从说明书、附图以及权利要求书变得明显。
附图说明
为了更清楚地说明本申请实施例的技术方案，下面将对本申请实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本申请的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1为本申请实施例提供的信号值与浓度的标准曲线的第一示例图；
图2为本申请实施例提供的信号值与浓度的标准曲线的第二示例图；
图3为本申请实施例提供的信号值与浓度的标准曲线的第三示例图；
图4为本申请实施例提供的信号值与浓度的标准曲线的第四示例图；
图5为本申请实施例提供的指标同质化换算方法的流程图；
图6为本申请实施例提供的目标弧线段示意图；
图7为本申请实施例提供的指标同质化换算装置的结构图；
图8为本申请实施例提供的电子设备的结构图。
图标：指标同质化换算装置700；获取模块710；确定模块720；校准模块730。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。
本申请实施例提供的技术方案可以应用于医疗场景中，即对医疗场景中涉及到的检测设备检测出的结果进行同质化。其中，检测设备可以为穿戴式检测设备，或其他检测设备。
检测设备检测样本获得的参数包括信号值和浓度，应用过程中，检测设备对样本进行检测，获得信号值和浓度，其中，浓度是最终需要的检测结果。浓度可以通过信号值与浓度之间的标准曲线进行计算，而标准曲线是各检测设备预先设置的，各检测设备中的标准曲线都存在一定差异。
示例性地，针对同一样本，在不同检测设备上进行检测，获得的浓度和信号值不同，但都有一定的对应关系。示例性地，测试设备1检测样本A，获得的信号值为5，浓度为6，测试设备2检测样本A，获得信号值为8，浓度值为15。尽管是同个样本，但通过不同检测设备检测出的信号值与浓度不同，造成该测试结果差异的是各检测设备之间检测信号敏感度、度量单位值及内置的标准曲线不同等众多种因素所造成。
在一些实施例中，该技术方案也可以应用于其他应用场景的检测设备的同质化，并不限定医疗场景。相应地，检测设备对样本检测，获得的参数也不限定信号值和浓度。也就是说，可以将本案的发明构思应用到其他应用场景中，也在本申请的保护范围内。
另外，本申请实施例提供的技术方案的硬件运行环境可以是检测设备自身，即检测设备自身对样本的检测结果进行同质化处理。作为另外一种实施方式还可以是单独的数据处理端，该设备处理端用于对不同的检测设备的检测结果进行同质化处理，该数据处理端可以是各类电子产品，在此不作限定。
示例性地，检测设备两台，一台为检测设备，另一台为参比设备。针对同一个样本，检测同一项指标，测试设备检测的信号值为
y C ，对应的浓度为
x C ，参比设备检测的信号值为
y B ，对应的浓度为
x B 。测试设备和参比设备均有两种检测方法可供选择，非竞争法（坐标轴两参数
x和
y为正相关关系曲线）和竞争法（坐标轴两参数
x和
y为负相关关系曲线）。两台检测设备对同一个样品进行同一项指标检测，共有4种情况：
情况1：参比设备与测试设备均采用非竞争法所得信号值与浓度标准曲线，如图1所示。
情况2：参比设备与测试设备均采用竞争法所得信号值与浓度标准曲线，如图2所示。
情况3：参比设备采用非竞争法、测试设备采用竞争法所得信号值与浓度标准曲线，如图3所示。
情况4：参比设备采用竞争法、测试设备采用非竞争法所得信号值与浓度标准曲线，如图4所示。
其中，在图1-4中，
L C 为测试设备的信号值与浓度对应的标准曲线，
L B 为参比设备的信号值与浓度对应的标准曲线。
如上述4种情况，测试设备和参比设备对同一个样品进行同一个指标检测，无论是选择相同的检测方法，还是选择不同的检测方法，信号值和浓度均不相同。若想实现测试设备与参比设备对相同样品、同一指标的检测结果同质化，关键是确定测试设备标准曲线上的点，然后在参比设备标准曲线上找到对应点，两个检测结果经过转换即可实现结果的同质化。
基于上述发明构思，请参照图5，本申请实施例提供的指标同质化换算方法的流程图，该方法包括：
S510，获取实测信号值。
其中，实测信号值为测试设备对测试样本进行检测获得的信号值。
一实施例中，信号值可以为发光值或吸光度等，具体可以根据检测项目决定，本申请并不以此为限。
S520，确定实测信号值在第一标准曲线的对应浓度。
其中，第一标准曲线为测试设备预设的浓度-信号值关系曲线，具体可以参照上文对标准曲线的说明，此处不再赘述。
需要说明的是，实测信号值为测试设备对测试样本进行检测所获得的信号值，基于该信号值和第一标准曲线，可获得与该信号值对应的浓度。
S530，根据第一标准曲线的第一端点和第一标准曲线的第二端点，获得第一标准曲线的第一目标弧线段，根据第一标准曲线的第一端点和第一标准曲线测试样本对应的端点，获得第一标准曲线的第二目标弧线段。
其中，第一标准曲线的第一端点和第二端点为不同比对样本对应的浓度-信号值。
S540，根据第二标准曲线的第一端点和所述第二标准曲线的第二端点，获得所述第二标准曲线的第三目标弧线段。
其中，第二标准曲线为参比设备预设的浓度-信号值关系曲线，第二标准曲线的第一端点与第一标准曲线的第一端点为同一比对样本对应的浓度-信号值，第二标准曲线的第二端点与第一标准曲线的第二端点为同一比对样本对应的浓度-信号值。
在本实施例中，S530和S540中第一标准区线的第一端点和第二端点以及第二标准曲线的第一端点和第二端点可以通过以下方法确定，该方法可以包括：
预先通过参比设备和测试设备对n个比对样本的信号值和浓度进行记录，获得第一标准曲线上的n个端点的浓度-信号值以及第二标准曲线上的n个端点的浓度-信号值。
可以理解，第一标准曲线和第二标准曲线都是由多个不同的曲线函数及各曲线函数对应的区间构成。各比对样本的选取可以是这些曲线函数的区间端点对应的浓度-信号值。
其中，可以根据厂家提供的检测设备（参比设备和测试设备）中对标准曲线的各区间段的说明获取标准曲线中各曲线函数以及对应区间，或者通过现有的比对样本选取方法获取标准曲线中各曲线函数以及对应区间，例如：申请号CN202210701442.3一种比对样本的确定方法及装置、存储介质、电子设备。
在确定测试样本在第一标准曲线的实测信号值以及对应浓度后，确定与该测试样本所处的区间，根据区间确定第一标准曲线的第一端点（区间端点）和第二端点（区间端点）。
进一步地，根据第一标准曲线与第二标准曲线上各端点的映射关系，也就可以通过第一标准曲线的第一端点确定第二标准曲线的第一端点，以及根据第一标准曲线的第二端点确定第二标准曲线的第二端点。
示例性地，请参照图6，S530的实现过程可以为：测试样本E的实测信号值及对应浓度落在对比样本
j和对比样本
j-1在第一标准曲线
L C 对应的
C j-1 端点和
C j 端点区间内，因此确定
C j-1 为第一标准曲线
L C 的第一端点，
C j 为第一标准曲线
L C 的第二端点，根据第一端点与第二端点在第一标准曲线
L C 上的位置，对第一标准曲线
L C 进行截取，获得第一目标弧线段。根据
C j-1 与E在第一标准曲线
L C 上的位置，对第一标准曲线
L C 进行截取，获得第二目标弧线段。
S540的实现过程可以为：根据对比样本
j和对比样本
j-1在第二标准曲线
L B 的浓度-信号值，可以确定
B j-1 为第二标准曲线
L B 的第一端点，
B j 为第二标准曲线
L B 的第二端点，通过对
B j-1 端点与
B j 端点在第二标准曲线
L B 上的部分进行截取，获得第三目标弧线段。
S550，根据第一目标弧线段、第二目标弧线段和第三目标弧线段，确定第四目标弧线段的弧长方程。
其中，第四目标弧线段为第二标准曲线的第一端点到测试样本在第二标准曲线对应的端点。
一实施例中，S550的具体实施方式可以是：根据第一目标弧线段、第二目标弧线段、第三目标弧线段各自的弧长，确定出第四目标弧线段的弧长与各目标弧线段的弧长的长度关系。基于弧长计算公式和长度关系，确定第四目标弧线段的弧长方程。
示例性地，请再次参考图6，由于第二标准曲线的第一端点(
B j-1 )与第一标准曲线的第一端点(
C j-1 )为同一比对样本（
j-1）获得的检测结果，第二标准曲线的第二端点（
B j ）与第一标准曲线的第二端点（
C j ）为同一比对样本（
j）获得的检测结果。此时，可以认为第一目标弧线段（）上的各端点与第二目标弧线段（）上的各端点是一一对应的关系且各端点在目标弧线段上的变化趋势均匀。
因此，可以得到以下各目标弧线段的弧长长度关系等式：
其中，E为第一标准曲线上测试样本对应的端点，为第二标准曲线上测试样本同质化后的端点。
根据弧长计算公式，对各弧长进行计算可以得到：
通过改变积分上限，得到弧线的长度为
其中，测试设备的第一标准曲线在的函数记为
f(x)，其导数记为
f'(x)。
同理，对于参比设备的第二标准曲线方程在的函数记为
g(x)，其导数记为
g'(x)。可以得到弧线的长度为
弧线的长度为
带入公式（1）换算，可得
其中，公式（2）为第四目标弧线段的弧长方程。
需要说明的是，为第一标准曲线的第一端点的横坐标，为第一标准曲线的第二端点的横坐标，为第二标准曲线的第一端点的横坐标，为第二标准曲线的第二端点的横坐标，为测试样本在第一标准曲线的横坐标，为测试样本在第二标准曲线的横坐标（实测信号值对应的同质化浓度）。
S560，根据第四目标弧线段的弧长方程，确定实测信号值对应的同质化浓度。
一实施例中，S560的具体实施过程，可以为：通过复合辛普森求积分公式和二分法，对第四目标弧线段的弧长方程进行处理，确定实测信号值对应的同质化浓度。
具体地，通过复合辛普森求积分公式对弧长方程进行转换，获得目标积分方程，通过二分法对目标积分方程进行求解，确定实测信号值对应的同质化浓度。
可以理解，考虑到寻找弧长公式中的积分的原函数以对该积分进行积分运算的计算可行性较低，因此，可以通过数值积分中的复合辛普森求积分公式对弧长方程进行转换。
示例性地，复合辛普森求积公式对弧长方程转换的过程如下：
将记为
F(x)，记为
G(x)，
m C ，
m B ，
m E 分别为，，区间内的小区间的划分个数，，，分别表示对应区间的步长值，从而可有对应的积分表达式：
其中，，
k C =0,1,2...
m C ，表示的中点。同理，对应的积分表达式为：
其中，，
k B =0,1,2...
m B ，表示的中点。
类似地，对应的积分表达式为：
其中，，
k E =0,1,2...
m E ，表示的中点。
对应地，可有对应的积分表达式为：
经过上述积分表达式进行转换后，除了未知，其他参数已知，因此可求得，因此，可以构建目标积分方程：
一实施例中，在上述获取的过程中，复合辛普森求积分公式的步长值的确定方式可以为：根据第一标准曲线的第一端点、第一标准曲线的第二端点、第一标准曲线中测试样本对应的端点、第二标准曲线的第一端点和第二标准曲线的第二端点，确定复合辛普森求积分公式的步长值。
示例性地，对的计算误差为：
其中，公式（3）由复合辛普森求积分公式的误差经数值计算中的误差传播计算所得，公式（3）中的各字符在前文已进行解释，在此不在过多赘述。
为便于计算，可以将公式（3）等同于。
为了数据计算处理效率和准确性，可以通过使用直线代替弧线的办法简单估算，和的值，即：
其中，为通过直线代替弧线估算的值，为通过直线代替弧线估算的值，为通过直线代替弧线估算的值。
在通过估算的方式得到A，B，C的值后，可以调整设置，，值，使得公式（3）得到的结果尽可能的小，公式（3）的计算结果越小，则表示后续计算结果的误差越小。
进一步地，若想继续减少误差，可比较A，B，C三个值的大小，取其最大的所对应的项进行增大。示例性地，若A>B>C，则我们首先可以增大，且保持和不变，以最有效地提高精度。
在调整好，，的值后，分别代入，，，以确定步长值。
可以理解，通过第一标准曲线的第一端点、第一标准曲线的第二端点、第一标准曲线测试样本对应的端点、第一目标弧线段的测试样本的端点、第二标准曲线的第一端点和第二标准曲线的第二端点，对步长值进行确定，以减少最终获取同质化浓度的误差。
进一步地，在获得目标积分方程后，通过二分法对目标积分方程进行求解，即可获得实测信号值对应的同质化浓度（即）。
需要说明的是，适用于二分法求解的函数为单调函数，而目标积分方程满足该要求。
示例性地，可以参照以下伪代码进行通过二分法对函数进行求解。
输入:点的横坐标;点的横坐标;计算得到的比值 复合辛普森求积分公式；待积函数精度eps; 迭代次数上限N。
1.,
2.
3.for do
4.
5.if then
6.
7.
8.if then
9.
10.else
11.
12.else
13.return
14.end for
以下是对上述伪代码的说明，设定区间【a，b】，其中a为点的横坐标，b为点的横坐标，d为同质化浓度与第二标准曲线的第二端点的浓度之差，如果浓度之差大于精度eps，则取a和b的中点值c，并将该中点值c和a和待积函数（即弧长方程）输入至复合辛普森求积分公式L中进行计算，得到中间同质化浓度std，若中间同质化浓度std小于，则重新对b赋值为c，若中间同质化浓度std大于或等于时，则重新对a赋值为c，一直重复此过程，直至最终循环结束得到的最终数值为同质化浓度。
可以理解，不同的检测设备分别为测试设备和参比设备，通过检测设备对测试样本，获得测试样本对应的信号值和浓度，从而确定测试样本在第一标准曲线的坐标，再根据第一标准曲线的第一端点和第二端点，获得第一标准曲线上的第一端点和第二端点连接的第一目标弧线段，根据第一标准曲线的第一端点与测试样本的坐标，获得第一标准曲线上的第一端点和测试样本坐标的第二目标弧线段，同理，获得在第二标准曲线上的第三目标弧线段。进一步地，由于第一标准曲线上的第一端点和第二标准曲线上的第一端点都是通过同一对比样本获得的不同检测结果，两点具有对应关系，同理，第一标准曲线上的第二端点与第二标准曲线上的第二端点对应。
在获得第一目标弧线段、第二目标弧线段、第三目标弧线段后，根据相似原理，建立第四目标弧线段的弧长方程，其中，第四目标弧线段为第二标准曲线的第一端点到测试样本在第二标准曲线对应的端点的弧线。再对第四目标弧线段的弧长方程进行求解，即可确定实测信号值在第二标准曲线上的浓度，即，确定实测信号值对应的同质化浓度，最终实现不同检测设备的检测结果的同质化。
一实施例中，在S560之后，指标同质化换算方法还可以包括：判断同质化浓度与第二标准曲线的第一端点或第二端点的浓度之差是否小于预设阈值，若是，则根据公式：或，对同质化浓度进行校准；其中，为校准后的同质化浓度，为第二目标弧线段的第一端点的浓度，为第二目标弧线段的第二端点的浓度，为同质化浓度。
可以理解，为了避免计算结果为无限小数或者迭代次数过多等情况，导致设备运算宕机，以及提高数据运算处理的效率，因此可以设定预设阈值，提高设备的数据处理效率。即，设定精度eps。在同质化浓度与第二标准曲线的第二端点的浓度之差小于精度eps时，可近似认为，使得最终的同质化浓度的数值可控，减少设备对同质化浓度的计算量，使得计算具有可行性，提高数据处理的效率。
可以理解，本领域技术人员可以根据实际应用场景进行灵活调整预设阈值。
基于同一发明构思，请参照图7，本申请实施例中还提供一种指标同质化换算装置700。该装置700，包括：获取模块710和确定模块720。
获取模块710用于获取实测信号值，实测信号值为测试设备对测试样本进行检测获得的信号值。确定模块720用于确定实测信号值在第一标准曲线的浓度；第一标准曲线为测试设备预设的浓度-信号值关系曲线。获取模块710还用于根据第一标准曲线的第一端点和第一标准曲线的第二端点，获得第一标准曲线的第一目标弧线段，根据第一标准曲线的第一端点和第一标准曲线测试样本对应的端点，获得第一标准曲线的第二目标弧线段，第一标准曲线的第一端点和第二端点为不同比对样本对应的端点。根据第二标准曲线的第一端点和第二标准曲线的第二端点，获得第二标准曲线的第三目标弧线段，第二标准曲线为参比设备预设的浓度-信号值关系曲线，第二标准曲线的第一端点与第一标准曲线的第一端点为同一比对样本对应的浓度-信号值，第二标准曲线的第二端点与第一标准曲线的第二端点为同一比对样本对应的浓度-信号值。确定模块720还用于根据第一目标弧线段、第二目标弧线段和第三目标弧线段，确定第四目标弧线段的弧长方程，第四目标弧线段为第二标准曲线的第一端点到测试样本在第二标准曲线对应的端点根据第四目标弧线段的弧长方程，确定实测信号值对应的同质化浓度。
一实施例中，确定模块720还用于根据第一目标弧线段、第二目标弧线段、第三目标弧线段各自的弧长，确定出第四目标弧线段的弧长与各目标弧线段的弧长的长度关系，基于弧长计算公式和长度关系，确定第四目标弧线段的弧长方程。
一实施例中，确定模块720还用于通过复合辛普森求积分公式和二分法，对弧长方程进行处理，确定实测信号值对应的同质化浓度。
一实施例中，确定模块720还用于通过复合辛普森求积分公式对弧长方程进行转换，获得目标积分方程，通过二分法对目标积分方程进行求解，确定实测信号值对应的同质化浓度。
一实施例中，确定模块720还用于根据第一标准曲线的第一端点、第一标准曲线的第二端点、第一标准曲线测试样本对应的端点、第一目标弧线段的测试样本的端点、第二标准曲线的第一端点和第二标准曲线的第二端点，确定复合辛普森求积分公式的步长值。
一实施例中，指标同质化换算装置700还包括校准模块730。
校准模块730用于判断同质化浓度与第二标准曲线的第一端点或第二端点的浓度之差是否小于预设阈值，若是，则根据公式：或，对同质化浓度进行校准；其中，为校准后的同质化浓度，为第二目标弧线段的第一端点的浓度，为第二目标弧线段的第二端点的浓度，为同质化浓度。
可以理解，本申请提供的指标同质化装置700与本申请提供的指标同质化方法对应，为使说明书简洁，相同或相似部分可以参照指标同质化方法部分的内容，在此不再赘述。
上述指标同质化装置700中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于服务器中的处理器中，也可以以软件形式存储于服务器中的存储器中，以便于处理器调用执行以上各个模块对应的操作。该处理器可以为中央处理单元（CPU)、微处理器、单片机等。
上述指标同质化方法和/或指标同质化装置可以实现为一种计算机可读指令的形式，计算机可读指令可以在如图8所示的电子设备上运行。
本申请实施例还提供的一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机可读指令，该处理器执行该程序时实现上述的指标同质化方法。
本实施例还提供了一种计算机可读存储介质，如软盘、光盘、硬盘、闪存、U盘、SD（SecureDigitalMemoryCard，安全数码卡）卡、MMC（MultimediaCard，多媒体卡）卡等，在该计算机可读存储介质中存储有实现上述各个步骤的一个或者多个程序，这一个或者多个程序可被一个或者多个处理器执行，以实现上述实施例一和/或实施例二中的编程界面显示控制方法。在此不再赘述。
在本申请所提供的实施例中，应该理解到，所揭露装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。
在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。
在本文中，多个是指两个或两个以上。
以上所述仅为本申请的实施例而已，并不用于限制本申请的保护范围，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。
Claims (5)
1.一种指标同质化换算方法，其特征在于，包括：
获取实测信号值，所述实测信号值为测试设备对测试样本进行检测获得的信号值；
确定所述实测信号值在第一标准曲线的对应浓度；所述第一标准曲线为所述测试设备预设的浓度-信号值关系曲线；
根据所述第一标准曲线的第一端点和所述第一标准曲线的第二端点，获得所述第一标准曲线的第一目标弧线段，根据所述第一标准曲线的第一端点和所述第一标准曲线所述测试样本对应的端点，获得所述第一标准曲线的第二目标弧线段，所述第一标准曲线的第一端点和第二端点为不同比对样本对应的浓度-信号值；
根据第二标准曲线的第一端点和所述第二标准曲线的第二端点，获得所述第二标准曲线的第三目标弧线段，所述第二标准曲线为参比设备预设的浓度-信号值关系曲线，所述第二标准曲线的第一端点与所述第一标准曲线的第一端点为同一比对样本对应的浓度-信号值，所述第二标准曲线的第二端点与所述第一标准曲线的第二端点为同一比对样本对应的浓度-信号值；
根据所述第一目标弧线段、所述第二目标弧线段和所述第三目标弧线段，确定第四目标弧线段的弧长方程，所述第四目标弧线段为所述第二标准曲线的第一端点到所述测试样本在所述第二标准曲线对应的端点；所述弧长方程为：
其中，，为所述第一目标弧线段的弧长，为所述第二目标弧线段的弧长，为所述第三目标弧线段的弧长，为所述第四目标弧线段的弧长；为所述第二标准曲线方程在所述第三目标弧线段上的函数的导数，为所述第二标准曲线的第一端点的横坐标，为所述测试样本在所述第二标准曲线的横坐标，即所述实测信号值对应的同质化浓度；
根据所述第四目标弧线段的弧长方程，确定所述实测信号值对应的同质化浓度，包括：
根据所述第一标准曲线的第一端点、所述第一标准曲线的第二端点、所述第一标准曲线中所述测试样本对应的端点、所述第二标准曲线的第一端点和所述第二标准曲线的第二端点，确定复合辛普森求积分公式的步长值；
通过所述复合辛普森求积分公式对所述弧长方程进行转换，获得目标积分方程；
通过二分法对所述目标积分方程进行求解，确定所述实测信号值对应的同质化浓度；
其中，所述根据所述第一标准曲线的第一端点、所述第一标准曲线的第二端点、所述第一标准曲线中所述测试样本对应的端点、所述第二标准曲线的第一端点和所述第二标准曲线的第二端点，确定复合辛普森求积分公式的步长值，包括：
通过如下误差公式计算所述复合辛普森求积分公式计算的误差，
，其中，A=，B=，C=，分别为区间内的小区间的划分个数，为所述第一标准曲线的第一端点的横坐标，为所述第一标准曲线的第二端点的横坐标，为所述测试样本在所述第一标准曲线的横坐标；
通过使用直线代替弧线的办法估算, 和的值；
通过估算的值分别估算所述A、所述B、所述C；
调整的值使得所述误差公式的值尽可能的小；
在调整好的值后，分别代入, ,，以确定所述复合辛普森求积分公式的步长值，其中，分别为区间对应的步长值。
2.如权利要求1所述的方法，其特征在于，在所述根据所述第四目标弧线段的弧长方程，确定所述实测信号值对应的同质化浓度之后，所述方法还包括：
判断所述同质化浓度与所述第二标准曲线的第一端点或第二端点的浓度之差是否小于预设阈值，
若是，则根据公式：或，对所述同质化浓度进行校准；
其中，为校准后的同质化浓度，为所述第二目标弧线段的第一端点的浓度，为所述第二目标弧线段的第二端点的浓度，为所述同质化浓度。
3.一种指标同质化换算装置，其特征在于，包括：
获取模块，用于获取实测信号值，所述实测信号值为测试设备对测试样本进行检测获得的信号值；
确定模块，用于确定所述实测信号值在第一标准曲线的对应浓度；所述第一标准曲线为所述测试设备预设的浓度-信号值关系曲线；
所述获取模块，还用于根据所述第一标准曲线的第一端点和所述第一标准曲线的第二端点，获得所述第一标准曲线的第一目标弧线段，根据所述第一标准曲线的第一端点和所述第一标准曲线所述测试样本对应的端点，获得所述第一标准曲线的第二目标弧线段，所述第一标准曲线的第一端点和第二端点为不同比对样本对应的浓度-信号值；根据第二标准曲线的第一端点和所述第二标准曲线的第二端点，获得所述第二标准曲线的第三目标弧线段，所述第二标准曲线为参比设备预设的浓度-信号值关系曲线，所述第二标准曲线的第一端点与所述第一标准曲线的第一端点为同一比对样本对应的浓度-信号值，所述第二标准曲线的第二端点与所述第一标准曲线的第二端点为同一比对样本对应的浓度-信号值；
所述确定模块，还用于根据所述第一目标弧线段、所述第二目标弧线段和所述第三目标弧线段，确定第四目标弧线段的弧长方程，所述第四目标弧线段为所述第二标准曲线的第一端点到所述测试样本在所述第二标准曲线对应的端点；根据所述第四目标弧线段的弧长方程，确定所述实测信号值对应的同质化浓度；所述弧长方程为：
其中，，为所述第一目标弧线段的弧长，为所述第二目标弧线段的弧长，为所述第三目标弧线段的弧长，为所述第四目标弧线段的弧长；为所述第二标准曲线方程在所述第三目标弧线段上的函数的导数，为所述第二标准曲线的第一端点的横坐标，为所述测试样本在所述第二标准曲线的横坐标，即所述实测信号值对应的同质化浓度；
所述确定模块具体用于：根据所述第一标准曲线的第一端点、所述第一标准曲线的第二端点、所述第一标准曲线中所述测试样本对应的端点、所述第二标准曲线的第一端点和所述第二标准曲线的第二端点，确定复合辛普森求积分公式的步长值；通过所述复合辛普森求积分公式对所述弧长方程进行转换，获得目标积分方程；通过二分法对所述目标积分方程进行求解，确定所述实测信号值对应的同质化浓度；
其中，所述根据所述第一标准曲线的第一端点、所述第一标准曲线的第二端点、所述第一标准曲线中所述测试样本对应的端点、所述第二标准曲线的第一端点和所述第二标准曲线的第二端点，确定复合辛普森求积分公式的步长值，包括：
通过如下误差公式计算所述复合辛普森求积分公式计算的误差，
，其中，A=，B=，C=，分别为区间内的小区间的划分个数，为所述第一标准曲线的第一端点的横坐标，为所述第一标准曲线的第二端点的横坐标，为所述测试样本在所述第一标准曲线的横坐标；
通过使用直线代替弧线的办法估算, 和的值；
通过估算的值分别估算所述A、所述B、所述C；
调整的值使得所述误差公式的值尽可能的小；
在调整好的值后，分别代入, ,，以确定所述复合辛普森求积分公式的步长值，其中，分别为区间对应的步长值。
4.一种电子设备，其特征在于，包括存储器以及处理器，所述存储器中存储有计算机程序指令，所述计算机程序指令被所述处理器读取并运行时，执行如权利要求1或2所述的指标同质化换算方法。
5.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序指令，所述计算机程序指令被计算机读取并运行时，执行如权利要求1或2所述的指标同质化换算方法。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210978315.8A CN115048613B (zh) | 2022-08-16 | 2022-08-16 | 一种指标同质化换算方法、装置、电子设备及存储介质 |
| PCT/CN2022/122661 WO2023151286A1 (zh) | 2022-02-14 | 2022-09-29 | 不同设备检测同一指标同质化换算方法、装置、电子设备 |
| EP22925644.1A EP4310701A1 (en) | 2022-02-14 | 2022-09-29 | Homogenization conversion method and apparatus for detection of same index by different devices, and electronic device |
| US18/003,654 US20230273973A1 (en) | 2022-02-14 | 2022-09-29 | Method and device for homogenization conversion of same index detected by different equipment, and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210978315.8A CN115048613B (zh) | 2022-08-16 | 2022-08-16 | 一种指标同质化换算方法、装置、电子设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115048613A CN115048613A (zh) | 2022-09-13 |
| CN115048613B true CN115048613B (zh) | 2023-05-12 |
Family
ID=83166885
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210978315.8A Active CN115048613B (zh) | 2022-02-14 | 2022-08-16 | 一种指标同质化换算方法、装置、电子设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115048613B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP4310701A1 (en) * | 2022-02-14 | 2024-01-24 | West China Hospital, Sichuan University | Homogenization conversion method and apparatus for detection of same index by different devices, and electronic device |
| CN117076938B (zh) * | 2023-10-17 | 2024-01-26 | 四川大学华西医院 | 最大样本数量确定方法、装置及设备 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101059335A (zh) * | 2007-05-10 | 2007-10-24 | 上海交通大学 | 三维空间连杆曲线匹配度的检测方法 |
| CN107657134A (zh) * | 2017-10-27 | 2018-02-02 | 天津百利机械装备集团有限公司中央研究院 | 一种基于Creo的弧面分度凸轮建模方法 |
| CN111368420A (zh) * | 2020-03-02 | 2020-07-03 | 厦门大学 | 一种gpu加速的参数曲线弧长计算及弧长参数化方法 |
| CN114166768A (zh) * | 2022-02-14 | 2022-03-11 | 四川大学华西医院 | 不同设备检测同一指标同质化换算方法、装置、电子设备 |
| CN114265367A (zh) * | 2021-12-31 | 2022-04-01 | 杭州电子科技大学 | 一种高精度参数路径插补方法及系统 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113130008A (zh) * | 2019-12-31 | 2021-07-16 | 科美诊断技术股份有限公司 | 测试结果确定方法及装置 |
| CN113283677B (zh) * | 2021-06-29 | 2022-08-26 | 中国平安人寿保险股份有限公司 | 指标数据处理方法、装置、设备及存储介质 |
-
2022
- 2022-08-16 CN CN202210978315.8A patent/CN115048613B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101059335A (zh) * | 2007-05-10 | 2007-10-24 | 上海交通大学 | 三维空间连杆曲线匹配度的检测方法 |
| CN107657134A (zh) * | 2017-10-27 | 2018-02-02 | 天津百利机械装备集团有限公司中央研究院 | 一种基于Creo的弧面分度凸轮建模方法 |
| CN111368420A (zh) * | 2020-03-02 | 2020-07-03 | 厦门大学 | 一种gpu加速的参数曲线弧长计算及弧长参数化方法 |
| CN114265367A (zh) * | 2021-12-31 | 2022-04-01 | 杭州电子科技大学 | 一种高精度参数路径插补方法及系统 |
| CN114166768A (zh) * | 2022-02-14 | 2022-03-11 | 四川大学华西医院 | 不同设备检测同一指标同质化换算方法、装置、电子设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115048613A (zh) | 2022-09-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115048613B (zh) | 一种指标同质化换算方法、装置、电子设备及存储介质 | |
| CN114166768B (zh) | 不同设备检测同一指标同质化换算方法、装置、电子设备 | |
| WO2006014509A2 (en) | Quantitative pcr data analysis system (qdas) | |
| WO2006014464A2 (en) | Method for quantitative pcr data analysis system (qdas) | |
| EP4310701A1 (en) | Homogenization conversion method and apparatus for detection of same index by different devices, and electronic device | |
| WO2022116673A1 (zh) | 色谱峰检测方法、装置、计算机设备和存储介质 | |
| WO2018161245A1 (zh) | 一种染色体变异的检测方法及装置 | |
| CN103234647A (zh) | 一种嵌入式系统的温度校准方法及系统 | |
| CN107924390A (zh) | 借助于校准曲线来确定样本的浓度范围的方法和系统 | |
| CN114780905B (zh) | 一种比对样本的确定方法及装置、存储介质、电子设备 | |
| CN112182067A (zh) | 个体身高预测方法、装置、电子设备及存储介质 | |
| CN114972339A (zh) | 用于推土机结构件生产异常检测的数据增强系统 | |
| CN107312850A (zh) | 一种pcr无效扩增的检测方法 | |
| US20110130989A1 (en) | System and method for identifying a peripheral component interconnect express signal | |
| CN108109675B (zh) | 一种实验室质量控制数据管理系统 | |
| CN115452101A (zh) | 一种仪表检定方法、装置、设备及介质 | |
| CN109920474B (zh) | 绝对定量方法、装置、计算机设备和存储介质 | |
| CN113761456A (zh) | 免疫荧光层析曲线的分析方法、装置及电子设备 | |
| CN113445099A (zh) | 槽液分析方法、装置及计算机可读存储介质 | |
| Wagner | The XL-mHG test for gene set enrichment | |
| CN110134265A (zh) | 滑鼠及其设定参数取得方法 | |
| CN113435058A (zh) | 配网自愈测试模型的数据降维方法、系统、终端及介质 | |
| CN114885231B (zh) | 一种通讯协议自适应的信号采集方法、系统、终端及介质 | |
| CN115062273B (zh) | 工业互联网的光电传感器精度控制方法及系统 | |
| JP2010193832A (ja) | 遺伝子解析方法および遺伝子解析システム |
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 |