A calibration method, device, storage medium and computer device

Abstract

Embodiments of the present application disclose a kind of scaling method, device, storage medium and computer equipment, wherein the method comprises: using Spline interpolation algorithm, the function relationship between concentration value and reaction degree is obtained, based on the function relationship and curve fitting data, the calibration sample of reaction degree measurement error is judged, so as to the reaction degree of the calibration sample of measurement error is re-measured, until all judging conditions meet preset condition, end judgment, according to the concentration of calibration sample and the standard curve obtained by corresponding measurement error-free reaction degree, by improving the accuracy of calibration sample reaction degree measurement, to improve the accuracy of the determined standard curve, and then improve the accuracy of scaling.

Description

Technical Field

[0001] The present invention relates to the technical field of spline interpolation curve fitting calibration, and in particular to a calibration method, device, storage medium and computer equipment. Background Art

[0002] At present, due to the increase in the number of calibration points during the measurement of calibrants by the instrument, as well as the occasional factors of the instrument itself during the testing process, the calibration method used for immune response detection has certain errors, making its accuracy unable to be guaranteed. As a result, the results of immune reaction detection are not accurate enough, affecting clinical diagnosis.

[0003] In clinical immunology testing, the calibration method is usually to measure a dose-response curve (i.e., a standard curve) using a series of concentration calibrators, and use this to infer the concentration of the unknown sample to be tested. The standard curve is generally obtained by using interpolation and fitting methods. Both methods are based on known discrete data to find a suitable function expression. The difference is that the function obtained by the interpolation method can pass through known points, while the fitting method only requires the function graph to be similar. Although the interpolation method can more accurately reflect the functional relationship of discrete data, it is more dependent on the precision and accuracy of the data. Therefore, there is an urgent need for a calibration method that can ensure the accuracy of the standard curve. Summary of the Invention

[0004] The main purpose of the present invention is to provide a calibration method, device, computer equipment and storage medium, which can solve the problem of low accuracy of standard curve in the prior art.

[0005] To achieve the above object, the present invention provides a calibration method in a first aspect, the method comprising:

[0006] Obtaining the reactivity of the calibrator, taking the corresponding relationship between the concentration value and the reactivity of the calibrator as a set of data, arranging each set of data according to the size of the concentration, and grouping the concentration values ​​of each two adjacent sets of data after the arrangement into a concentration interval, obtaining N-1 concentration intervals, where N is the total number of calibrators;

[0007] For a target concentration interval, the functional relationship between the concentration values ​​and the reactivity in two adjacent groups of data corresponding to the target concentration interval is calculated based on the Spline interpolation algorithm to obtain the functional relationship corresponding to the target concentration interval; wherein the target concentration interval is any one of the N-1 concentration intervals;

[0008] According to the function relationship corresponding to the target concentration interval, extreme points of the function relationship in the target concentration interval are calculated, if there is an extreme point in the target concentration interval, whether there is a target calibrant with incorrect reaction rate in the adjacent two groups of data corresponding to the target concentration interval is determined according to the comparison result of the extreme point and the maximum concentration value and the minimum concentration value constituting the target concentration interval and the number of extreme points, if the target calibrant exists, the reaction rate of the target calibrant is re-measured, the reaction rate re-measured is used to replace the reaction rate of the target calibrant, and the step of arranging each group of data according to the size of the concentration is re-executed until there is no extreme point in the target concentration interval;

[0009] If there is no extreme point in the target concentration interval, curve fitting data are obtained according to the function relationship of each concentration interval, a first-order difference array of the curve fitting data is calculated, a second-order difference array is calculated according to the first-order difference array, whether there is a calibrant with incorrect reaction rate is determined according to the second-order difference array, if there is a calibrant with incorrect reaction rate, the reaction rate of the calibrant with incorrect reaction rate is re-measured, if there is no calibrant with incorrect reaction rate, the judgment is ended, and a standard curve is obtained based on a Spline interpolation algorithm according to the concentration of the calibrant and the reaction rate measured correctly.

[0010] In combination with the first aspect, in a possible implementation manner, before the concentration values in each adjacent two groups of data are grouped into a concentration interval to obtain N-1 concentration intervals, the following is included: the reaction rate in the i+1th group of data and the reaction rate in the ith group of data are compared, when each group of data is sorted from small to large according to the concentration value of the calibrant, if the reaction rate in the i+1th group of data is less than the reaction rate in the ith group of data, it is determined that the reaction rate in the i+1th group of data is measured with error, and the reaction rate of the corresponding calibrant in the i+1th group of data is re-measured; the reaction rate re-measured for the calibrant corresponding to the i+1th group of data is used to replace the original reaction rate in the i+1th group of data, and the step of comparing the reaction rate in the i+1th group of data and the reaction rate in the ith group of data is re-executed until the reaction rate in the i+1th group of data is greater than the reaction rate in the ith group of data; wherein, the value of i is from 1 to N-1; when each group of data is sorted from large to small according to the concentration value of the calibrant, if the reaction rate in the i+1th group of data is greater than the reaction rate in the ith group of data, it is determined that the reaction rate in the ith group of data is measured with error, and the reaction rate of the corresponding calibrant in the ith group of data is re-measured; the reaction rate re-measured for the calibrant corresponding to the ith group of data is used to replace the original reaction rate in the ith group of data, and the step of comparing the reaction rate in the i+1th group of data and the reaction rate in the ith group of data is re-executed until the reaction rate in the i+1th group of data is less than the reaction rate in the ith group of data.

[0011] With reference to the first aspect, in a possible implementation manner, the determining, according to the comparison result of the extreme value point and the maximum concentration value and the minimum concentration value constituting the target concentration interval and the number of extreme value points, whether there is a target calibrator with a wrong reaction degree measurement in the adjacent two groups of data corresponding to the target concentration interval, comprises: comparing the extreme value point with the maximum concentration value and the minimum concentration value, if the difference between the extreme value point and the minimum concentration value is greater than the difference between the maximum concentration value and the extreme value point, and the number of extreme value points is greater than a first threshold, it is determined that the reaction degree measurement of the calibrator with the maximum concentration value in the adjacent two groups of data corresponding to the target concentration interval is wrong; if the difference between the extreme value point and the minimum concentration value is less than the difference between the maximum concentration value and the extreme value point, and the number of extreme value points is greater than the first threshold, it is determined that the reaction degree measurement of the calibrator with the minimum concentration value in the adjacent two groups of data corresponding to the target concentration interval is wrong.

[0012] With reference to the first aspect, in a possible implementation manner, the determining, according to the comparison result of the extreme value point and the maximum concentration value and the minimum concentration value constituting the target concentration interval and the number of extreme value points, whether there is a target calibrator with a wrong reaction degree measurement in the adjacent two groups of data corresponding to the target concentration interval, comprises: comparing the extreme value point with the maximum concentration value and the minimum concentration value, if the difference between the extreme value point and the minimum concentration value is greater than the difference between the maximum concentration value and the extreme value point, and the number of extreme value points is greater than a first threshold, it is determined that the reaction degree measurement of the calibrator with the maximum concentration value in the adjacent two groups of data corresponding to the target concentration interval is wrong; if the difference between the extreme value point and the minimum concentration value is less than the difference between the maximum concentration value and the extreme value point, and the number of extreme value points is greater than the first threshold, it is determined that the reaction degree measurement of the calibrator with the minimum concentration value in the adjacent two groups of data corresponding to the target concentration interval is wrong.

[0013] In a possible implementation of the first aspect, the function relationship according to the target concentration interval is calculated by deriving the function relationship to obtain a derivative function of the function relationship, and the extreme point of the target concentration interval is calculated according to the derivative function.

[0014] To achieve the above object, the second aspect of the present application provides a calibration device, which comprises:

[0015] The interval division module is configured to obtain the reaction degree of the calibrators, take the corresponding relationship between the concentration values of the calibrators and the reaction degrees as a group of data, arrange each group of data according to the concentration values, group the concentration values in each adjacent two groups of data after the arrangement as a concentration interval, and obtain N-1 concentration intervals, where N is the total number of the calibrators.

[0016] The function establishment module is configured to calculate, for a target concentration interval, a function relationship between the concentration values and the reaction degrees in the adjacent two groups of data corresponding to the target concentration interval based on a Spline interpolation algorithm, and obtain the function relationship corresponding to the target concentration interval, where the target concentration interval is any one of the N-1 concentration intervals.

[0017] The first judgment module is configured to calculate an extreme point of the function relationship in the target concentration interval according to the function relationship corresponding to the target concentration interval, and if the extreme point exists in the target concentration interval, judge whether there is a target calibrator with a wrong reaction degree in the adjacent two groups of data corresponding to the target concentration interval according to the comparison result of the extreme point and the maximum and minimum concentration values constituting the target concentration interval and the number of the extreme points, if the target calibrator exists, re-measure the reaction degree of the target calibrator, replace the reaction degree of the target calibrator with the re-measured reaction degree, and re-perform the step of arranging each group of data according to the concentration values until the target concentration interval does not have the extreme point.

[0018] The second judgment module is configured to, if the target concentration interval does not have the extreme point, obtain curve fitting data according to the function relationships of the concentration intervals, calculate a first-order difference array of the curve fitting data, calculate a second-order difference array according to the first-order difference array, judge whether there is a calibrator with a wrong reaction degree according to the second-order difference array, if the calibrator with the wrong reaction degree exists, re-measure the reaction degree of the calibrator with the wrong reaction degree, if the calibrator with the wrong reaction degree does not exist, end the judgment, and obtain a standard curve based on the Spline interpolation algorithm according to the concentration of the calibrators and the reaction degrees without errors.

[0019] To achieve the above object, the third aspect of the present application provides a computer readable storage medium storing a computer program, where the computer program is executed by a processor to make the processor perform the following steps:

[0020] obtaining the reaction degree of the calibration samples, taking the corresponding relationship between the concentration values and the reaction degrees of the calibration samples as a group of data, arranging each group of data according to the size of the concentration, grouping the concentration values in each adjacent two groups of data after the arrangement into a concentration interval, and obtaining N-1 concentration intervals, wherein N is the total number of the calibration samples;

[0021] for a target concentration interval, calculating the functional relationship between the concentration values and the reaction degrees in the adjacent two groups of data corresponding to the target concentration interval based on a Spline interpolation algorithm, and obtaining the functional relationship corresponding to the target concentration interval; wherein the target concentration interval is any one of the N-1 concentration intervals;

[0022] calculating the extreme points of the functional relationship in the target concentration interval according to the functional relationship corresponding to the target concentration interval, if there is an extreme point in the target concentration interval, judging whether there is a target calibration sample with a wrong reaction degree measurement in the adjacent two groups of data corresponding to the target concentration interval according to the comparison result of the extreme point and the maximum and minimum concentration values constituting the target concentration interval and the number of extreme points, if there is the target calibration sample, re-measuring the reaction degree of the target calibration sample, replacing the reaction degree of the target calibration sample with the re-measured reaction degree, and re-executing the step of arranging each group of data according to the size of the concentration until there is no extreme point in the target concentration interval;

[0023] if there is no extreme point in the target concentration interval, obtaining the curve fitting data according to the functional relationship of each concentration interval, calculating the first-order difference array of the curve fitting data, calculating the second-order difference array according to the first-order difference array, judging whether there is a calibration sample with a wrong reaction degree measurement according to the second-order difference array, if there is a calibration sample with a wrong reaction degree measurement, re-measuring the reaction degree of the calibration sample with the wrong reaction degree measurement, if there is no calibration sample with a wrong reaction degree measurement, ending the judgment, and obtaining the standard curve based on the Spline interpolation algorithm according to the concentration of the calibration samples and the reaction degrees without wrong measurement.

[0024] To achieve the above object, the fourth aspect of the present application provides a computer device, comprising a memory and a processor, the memory stores a computer program, and the computer program is executed by the processor to make the processor execute the following steps:

[0025] obtaining the reaction degree of the calibration samples, taking the corresponding relationship between the concentration values and the reaction degrees of the calibration samples as a group of data, arranging each group of data according to the size of the concentration, grouping the concentration values in each adjacent two groups of data after the arrangement into a concentration interval, and obtaining N-1 concentration intervals, wherein N is the total number of the calibration samples;

[0026] For the target concentration interval, based on the Spline interpolation algorithm, a function relationship between the concentration values and the reaction degrees in the adjacent two groups of data corresponding to the target concentration interval is calculated to obtain a function relationship corresponding to the target concentration interval; wherein the target concentration interval is any one of the N-1 concentration intervals;

[0027] According to the function relationship corresponding to the target concentration interval, extreme points of the function relationship in the target concentration interval are calculated, if there is an extreme point in the target concentration interval, then according to the comparison result of the extreme point and the maximum concentration value and the minimum concentration value constituting the target concentration interval and the number of extreme points, it is judged whether there is a target calibrator with incorrect reaction degree measurement in the adjacent two groups of data corresponding to the target concentration interval, if there is the target calibrator, then the reaction degree of the target calibrator is re-measured, the reaction degree obtained by re-measuring is used to replace the reaction degree of the target calibrator, and the step of arranging each group of data according to the size of the concentration is re-executed until there is no extreme point in the target concentration interval;

[0028] If there is no extreme point in the target concentration interval, then according to the function relationship of each concentration interval, curve fitting data is obtained, a first-order difference array of the curve fitting data is calculated, a second-order difference array is calculated according to the first-order difference array, whether there is a calibrator with incorrect reaction degree measurement is judged according to the second-order difference array, if there is a calibrator with incorrect reaction degree measurement, then the reaction degree of the incorrect calibrator is re-measured, if there is no calibrator with incorrect reaction degree measurement, then the judgment is ended, and based on the Spline interpolation algorithm, a standard curve is obtained according to the concentration of the calibrator and the correct reaction degree.

[0029] By adopting the embodiment of the present application, the following beneficial effects are achieved:

[0030] The application provides a calibration method, which comprises the following steps: obtaining the reaction degree of a calibration sample, taking the corresponding relationship between the concentration value and the reaction degree of the calibration sample as a group of data, arranging each group of data according to the concentration, grouping the concentration values in each adjacent two groups of data after the arrangement into a concentration interval, and obtaining N-1 concentration intervals, wherein N is the total number of the calibration samples; for a target concentration interval, calculating the functional relationship between the concentration value and the reaction degree in the adjacent two groups of data corresponding to the target concentration interval based on a Spline interpolation algorithm, and obtaining the functional relationship corresponding to the target concentration interval; wherein the target concentration interval is any one of the N-1 concentration intervals; calculating the extreme points of the functional relationship in the target concentration interval according to the functional relationship corresponding to the target concentration interval, and if there is an extreme point in the target concentration interval, judging whether there is a target calibration sample with a wrong reaction degree in the adjacent two groups of data corresponding to the target concentration interval according to the comparison result of the extreme point and the maximum concentration value and the minimum concentration value constituting the target concentration interval and the number of the extreme points, if there is a target calibration sample, re-measuring the reaction degree of the target calibration sample, replacing the reaction degree of the target calibration sample with the re-measured reaction degree, and re-executing the step of calculating the functional relationship between the concentration value and the reaction degree of the target concentration interval based on the Spline interpolation algorithm until there is no extreme point in the target concentration interval; if there is no extreme point in the target concentration interval, obtaining curve fitting data according to the functional relationship of each concentration interval, judging whether there is a calibration sample with a wrong reaction degree according to the curve fitting data, if there is a calibration sample with a wrong reaction degree, re-measuring the reaction degree of the calibration sample with the wrong reaction degree, if there is no calibration sample with a wrong reaction degree, ending the judgment, and obtaining a standard curve according to the concentration of the calibration sample and the correct reaction degree when there is no calibration sample with a wrong reaction degree. Based on the technical scheme, the Spline interpolation algorithm is used to obtain the functional relationship between the concentration value and the reaction degree, the calibration sample with a wrong reaction degree is judged based on the functional relationship and the curve fitting data, the reaction degree of the calibration sample with a wrong reaction degree is re-measured until all the judgment conditions meet the preset conditions, the judgment is ended, the accuracy of the reaction degree measurement of the calibration sample is improved, the accuracy of the obtained standard curve is improved, and the calibration accuracy is improved. BRIEF DESCRIPTION OF DRAWINGS

[0031] In order to more clearly illustrate the technical solutions of the embodiments of the present application or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or the prior art description. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without any creative effort.

[0032] In the formula, the concentration of the calibration sample is represented by C, the reaction degree of the calibration sample is represented by R, and the function is represented by F.

[0033] Figure 1 a flowchart of a calibration method according to an embodiment of the present application;

[0034] Figure 2 a structural block diagram of a calibration device according to an embodiment of the present application;

[0035] Figure 3 a structural block diagram of a computer device according to an embodiment of the present application. DETAILED DESCRIPTION

[0036] The technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all the other embodiments obtained by a person of ordinary skill in the art without creative work fall within the protection scope of the present application.

[0037] The embodiments of the present application provide a calibration method, referring to Figure 1 , Figure 1 A flowchart of a calibration method according to an embodiment of the present application is shown in FIG. 1, and the specific steps of the method are as follows: Figure 1

[0038] In step S101, the reaction degree of the calibration sample is obtained, and the corresponding relationship between the concentration value and the reaction degree of the calibration sample is taken as a group of data. Each group of data is arranged according to the size of the concentration, and the concentration values in each adjacent two groups of arranged data are grouped into a concentration interval, thereby obtaining N-1 concentration intervals.

[0039] Wherein, N is the total number of calibration samples.

[0040] The reaction degree of the calibration sample is measured on the instrument. The concentration values of the calibration samples are not the same, that is, the reaction degrees of calibration samples with different concentrations are measured on the instrument. In addition, due to the limitation of the measurement accuracy of the instrument itself and the interference of some unknown factors, errors may exist in the measurement. In order to improve the accuracy of the reaction degree of the calibration sample, the reaction degree of each concentration calibration sample is measured multiple times on the instrument, and then the average value of the measured multiple reaction degrees is calculated as the final reaction degree of the calibration sample. In this embodiment, the reaction degree of each concentration calibration sample is measured at least three times on the instrument. For example, for a calibration sample with a concentration value of C1, the measured reaction degrees are R1 1 , R1 2 , R1 3 , the average value R1 1 , R1 2 , R1 3 is calculated, and R1 is taken as the reaction degree of the calibration sample. ​

[0041] In the embodiment, after all the calibration samples to be measured are measured, the host computer stores the corresponding relationship between the concentration values of the calibration samples and the reaction degrees as a set of data in the first array. For example, the calibration samples to be measured include a calibration sample 1 with a concentration C1, a calibration sample 2 with a concentration C2, and a calibration sample 3 with a concentration C3. The measured reaction degree of the calibration sample 1 is R1, the reaction degree of the calibration sample 2 is R2, and the reaction degree of the calibration sample 3 is R3. Then, C1 and R1 of the calibration sample 1 are taken as a set of data, C2 and R2 of the calibration sample 3 are taken as a set of data, and C3 and R3 of the calibration sample 3 are taken as a set of data. The three sets of data are stored in the first array.

[0042] The host computer sends the first array to the algorithm module. After receiving the first array from the host computer, the algorithm module sorts each set of data in the first array according to the concentration values of the calibration samples. The sorting can be from small to large according to the concentration values of the calibration samples, or from large to small according to the concentration values of the calibration samples.

[0043] Step S1011, the reaction degree in the i+1th set of data is compared with the reaction degree in the ith set of data. When each set of data is sorted from small to large according to the concentration values of the calibration samples, if the reaction degree in the i+1th set of data is less than the reaction degree in the ith set of data, it is determined that the measurement of the reaction degree in the i+1th set of data has an error, and the reaction degree of the corresponding calibration sample in the i+1th set of data is re-measured. The reaction degree re-measured for the calibration sample corresponding to the i+1th set of data replaces the original reaction degree in the i+1th set of data, and the step of comparing the reaction degree in the i+1th set of data with the reaction degree in the ith set of data is re-executed until the reaction degree in the i+1th set of data is greater than the reaction degree in the ith set of data.

[0044] Step S1012, when each set of data is sorted from large to small according to the concentration values of the calibration samples, if the reaction degree in the i+1th set of data is greater than the reaction degree in the ith set of data, it is determined that the measurement of the reaction degree in the ith set of data has an error, and the reaction degree of the corresponding calibration sample in the ith set of data is re-measured. The reaction degree re-measured for the calibration sample corresponding to the ith set of data replaces the original reaction degree in the ith set of data, and the step of comparing the reaction degree in the i+1th set of data with the reaction degree in the ith set of data is re-executed until the reaction degree in the i+1th set of data is less than the reaction degree in the ith set of data.

[0045] After the sorting, the reaction rates in each adjacent two groups of data are compared. Specifically, when each group of data is sorted in ascending order according to the concentration values of the calibrators, if the reaction rate in the i+1th group of data is less than the reaction rate in the ith group of data, it indicates that an abnormality of the calibration point occurs, that is, there is a calibrator with a wrong reaction rate measurement, and it is determined that the reaction rate measurement of the ith group of data has a large deviation and needs to be re-measured. At the same time, the determination information is reported to the host, such as reporting to the host that the reaction rate measurement of the calibrator corresponding to the ith group of data is wrong and needs to be re-measured. When each group of data is sorted in descending order according to the concentration values of the calibrators, if the reaction rate in the i+1th group of data is greater than the reaction rate in the ith group of data, it indicates that an abnormality of the calibration point occurs, that is, there is a calibrator with a wrong reaction rate measurement, and it is determined that the reaction rate measurement of the ith group of data has a large deviation and needs to be re-measured. At the same time, the determination information is reported to the host, such as reporting to the host that the reaction rate measurement of the calibrator corresponding to the ith group of data is wrong and needs to be re-measured.

[0046] After the instrument re-measures the reaction rate of the calibrator with wrong measurement, the host re-provides the re-measured reaction rate to the algorithm module, the algorithm module updates the reaction rate in the data corresponding to the calibrator with wrong measurement, that is, replaces the original reaction rate in the data with the re-measured reaction rate, and re-executes the comparison of the reaction rates in each adjacent two groups of data. If each group of data is sorted in ascending order according to the concentration values of the calibrators, until the reaction rate in the i+1th group of data is greater than the reaction rate in the ith group of data. If each group of data is sorted in descending order according to the concentration values of the calibrators, until the reaction rate in the i+1th group of data is less than the reaction rate in the ith group of data.

[0047] When the comparison results of the reaction rates in each adjacent group of data all meet the conditions, the concentration values in each adjacent two groups of data are grouped into a concentration interval, and N-1 concentration intervals are obtained, that is, the concentration values C i (i∈1,2,3,L,N) in the two sorted groups of data are grouped into a small interval, so that N-1 concentration intervals are divided as (C1, C2), (C2, C3), …, (C N-1 , C N ), where each interval can be represented by x j (j∈1,2,3,L,N-1). It should be noted that the representation symbol is not limited to the letter symbol x.

[0048] For example, if the concentration values in each group of data after arrangement are in the order of 0.1, 0.2, 0.3, 0.4, and 0.5, the four concentration intervals divided are (0.1, 0.2), (0.2, 0.3), (0.3, 0.4), and (0.4, 0.5), respectively. If the concentration values in each group of data after arrangement are in the order of 0.5, 0.4, 0.3, 0.2, and 0.1, the four concentration intervals divided are also (0.1, 0.2), (0.2, 0.3), (0.3, 0.4), and (0.4, 0.5), respectively.

[0049] In step S102, for the target concentration interval, a function relationship between the concentration values and the reaction degrees in the adjacent two groups of data corresponding to the target concentration interval is calculated based on the Spline interpolation algorithm, to obtain a function relationship corresponding to the target concentration interval.

[0050] In this embodiment, the function relationship S(x j ) between the concentration values and the reaction degrees in the adjacent two groups of data corresponding to each concentration interval is calculated. j (j∈1,2,3,L,N-1), j represents the jth concentration interval in the concentration intervals, Y j (j∈1,2,3,L,N-1) represents the function relationship corresponding to the jth concentration interval. It can be understood that, the function relationship between the concentration values and the reaction degrees in the adjacent two groups of data corresponding to the target concentration interval is calculated based on the Spline interpolation algorithm, to obtain a function relationship corresponding to the target concentration interval, wherein the target concentration interval is any one of the N-1 concentration intervals.

[0051] In step S103, an extreme point of the function relationship in the target concentration interval is calculated according to the function relationship corresponding to the target concentration interval. If there is an extreme point in the target concentration interval, whether there is a target calibrator with a measurement error in the reaction degree in the adjacent two groups of data corresponding to the target concentration interval is determined according to the comparison result of the extreme point and the maximum and minimum concentration values constituting the target concentration interval and the number of extreme points. If there is the target calibrator, the reaction degree of the target calibrator is re-measured, the reaction degree re-measured is used to replace the reaction degree of the target calibrator, and the step of arranging each group of data according to the concentration is re-executed until there is no extreme point in the target concentration interval.

[0052] Further, the derivative of the function relationship S(x j ) corresponding to the target concentration interval is obtained, to obtain the derivative function The extreme point C site, determine whether there is a target calibrator with incorrect reactivity measurement in two adjacent groups of data corresponding to the target concentration range according to the comparison result of the extreme point with the maximum concentration value and the minimum concentration value constituting the target concentration range and the number of extreme points. Among them, if the target concentration range is (a, b), the maximum concentration value is b and the minimum concentration value is a.

[0053] Step S1031: Compare the extreme point with the maximum concentration value and the minimum concentration value. If the difference between the extreme point and the minimum concentration value is greater than the difference between the maximum concentration value and the extreme point, and the number of extreme points is greater than the first threshold, then determine that the reactivity measurement of the calibrator with the maximum concentration value in two adjacent groups of data corresponding to the target concentration range is incorrect.

[0054] Step S1032: If the difference between the extreme point and the minimum concentration value is less than the difference between the maximum concentration value and the extreme point, and the number of extreme points is greater than the first threshold, then determine that the reactivity measurement of the calibrator with the minimum concentration value in two adjacent groups of data corresponding to the target concentration range is incorrect.

[0055] Specifically, if there is an extreme point C site , it indicates that there is an abnormal calibration point situation, that is, there is a target calibrator with incorrect reactivity measurement. Then record the total number of extreme points C number , and compare the extreme point C site with the minimum concentration value C L of the target concentration range, and the maximum concentration value C H . If |C site - C L | > |C H - C site |, then store C H in the second array, otherwise store C L in the second array.

[0056] If C number > TH1, the algorithm module reports a signal of calibration failure to the host, such as reporting to the host "Calibration failed, please calibrate again!". After receiving this signal, the host clears all reactivities in the first array, re-measures the reactivities of all calibrators and stores them in the first array for calibration, that is, re-executes step S101 until C number < TH1 or there is no extreme point C site .

[0057] If C numberIf TH1, the host reports that the reaction degree measurement of the target calibration sample with the concentration of the value in the second array is incorrect, and the reaction degree of these target calibration samples needs to be re-measured. Wherein, TH1 is a first threshold value, which can be determined according to experience. For example, when C H is stored in the second array, it is determined that the reaction degree measurement of the calibration sample with the maximum concentration value C H is incorrect, when C L is stored in the second array, it is determined that the reaction degree measurement of the calibration sample with the minimum concentration value C L is incorrect.

[0058] After the instrument re-measures the reaction degree of the calibration sample with incorrect measurement, the host re-provides the re-measured reaction degree to the algorithm module, the algorithm module updates the reaction degree in the data corresponding to the calibration sample with incorrect measurement in the first array, that is, replaces the original reaction degree in the data with the re-measured reaction degree, and re-executes the step of arranging each group of data according to the size of the concentration until there is no extreme point in the target concentration interval. For example, the target concentration interval is (a, b), the maximum concentration value is b, the calibration sample corresponding to the maximum concentration value b is C2, the minimum concentration value is a, and the calibration sample corresponding to the minimum concentration value a is C1. When it is determined that the reaction degree measurement of the calibration sample with the maximum concentration value b is incorrect, that is, the reaction degree measurement of C2 is incorrect, the reaction degree of C2 is re-measured to obtain a new reaction degree of C2, and the original reaction degree in the data corresponding to C2 in the first array is replaced with the new reaction degree of C2.

[0059] In step S104, if there is no extreme point in the target concentration interval, the curve fitting data is obtained according to the function relationship of each concentration interval, the first-order difference array of the curve fitting data is calculated, the second-order difference array is calculated according to the first-order difference array, and it is judged whether there is a calibration sample with incorrect reaction degree measurement according to the second-order difference array. If there is a calibration sample with incorrect reaction degree measurement, the reaction degree of the calibration sample with incorrect measurement is re-measured, if there is no calibration sample with incorrect reaction degree measurement, the judgment is ended, and the standard curve is obtained based on the Spline interpolation algorithm according to the concentration of the calibration sample and the correct reaction degree.

[0060] If there is no extreme point C site , that is, C number =0 is satisfied, the curve fitting data is obtained according to the function relationship S(x j ) = Y j (j∈1,2,3,L,N-1) of each concentration interval, the first-order difference array R of the curve fitting data is calculated, and the second-order difference array R1 is calculated according to the first-order difference array R. Whether there is a calibration sample with incorrect reaction degree measurement is judged according to the second-order difference array R1.

[0061] Wherein, according to the second-order difference array R1, whether there is a calibration sample with incorrect reactivity measurement is determined by the following specific steps:

[0062] Step S1041, the product between the first second-order difference array and the second second-order difference array is calculated.

[0063] Step S1042, if the product between the first second-order difference array and the second second-order difference array is not less than 0, it is determined that there is no calibration sample with incorrect reactivity measurement, and the judgment is ended.

[0064] Step S1043, if the product between the first second-order difference array and the second second-order difference array is less than 0, the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is calculated, if the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is less than the second threshold value, it is determined that there is no calibration sample with incorrect reactivity measurement, and the judgment is ended; if the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is greater than the second threshold value, the absolute value of the difference between the concentration value and k in each concentration interval is calculated, and it is determined that the reactivity measurement of the calibration sample corresponding to the minimum value in the absolute value of the difference between the concentration value and k is incorrect, the reactivity of the calibration sample with incorrect reactivity measurement is re-measured, the reactivity of the calibration sample with incorrect reactivity measurement re-measured is replaced with the original reactivity in the data corresponding to the calibration sample with incorrect reactivity measurement, and the step of arranging each group of data according to the size of the concentration is re-executed until the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is less than the second threshold value, or the product between the first second-order difference array and the second second-order difference array is not less than 0.

[0065] Wherein, the first second-order difference array refers to the second-order difference array under k+1, the second second-order difference array refers to the second-order difference array under k, and k is a preset gradient value.

[0066] Specifically, if R1[k+1]xR1[k]≥0 is met in the R1 array, it means that there is no abnormal calibration point, i.e., there is no calibration sample with incorrect reactivity measurement, the judgment is ended, and finally, the standard curve is obtained according to all the concentrations and the corresponding reactivity stored in the first array after the judgment is ended.

[0067] If R1[k+1]xR1[k]<0 is satisfied in the R1 array, indicating that a calibration sample with a wrong reaction rate is present, the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is calculated, where R1[k+1] is the first second-order difference array, R1[k] is the second second-order difference array, and k represents a preset concentration gradient divided according to the maximum concentration value and the minimum concentration value in the concentration values of all calibration samples. For example, all calibration samples include a calibration sample with a concentration value of 0.1, a calibration sample with a concentration value of 0.4, and a calibration sample with a concentration value of 0.6. The maximum concentration value is 0.6, and the minimum concentration value is 0.1. The concentration gradient from 0.1 to 0.6 can be divided into 0.1, 0.2, 0.3, 0.4, 0.5, and 0.6, that is, k can be 0.1, 0.2, 0.3, 0.4, 0.5, and 0.6.

[0068] If |R1[k+1]|+|R1[k]|>TH2 is present, R j j -k|(j∈1,2,3,L,N) is calculated, and R j The minimum value of C j is stored in the second array, and a signal indicating that the reaction rate of the calibration sample with the concentration value in the second array is wrong is reported to the host, such as reporting to the host that the reaction rate of the calibration sample with the concentration value in the second array is wrong, and the reaction rate of these calibration samples needs to be measured again.

[0069] After the instrument re-measures the reaction rate of the calibration sample with the wrong reaction rate, the host re-provides the re-measured reaction rate to the algorithm module, the algorithm module updates the reaction rate in the data corresponding to the calibration sample with the wrong reaction rate, that is, replaces the original reaction rate in the data with the re-measured reaction rate, and re-executes the step of arranging each group of data according to the size of the concentration until |R1[k+1]|+|R1[k]|<TH2.

[0070] If |R1[k+1]|+|R1[k]|<TH2 is satisfied, the judgment is ended, and at this time, all the concentrations corresponding to the reaction rates stored in the first array are valid. Therefore, a standard curve can be obtained based on the Spline interpolation algorithm according to the concentration of the calibration sample and the measured reaction rate, that is, a standard curve is obtained according to all the concentrations and the corresponding reaction rates stored in the first array after the end of the judgment, to complete the calibration.

[0071] ​Based on the above method, by acquiring the reaction degree of the calibration samples, taking the corresponding relationship between the concentration value of the calibration samples and the reaction degree as a set of data, arranging each set of data according to the size of the concentration, grouping the concentration values in each adjacent two sets of data after the arrangement into a concentration interval, N-1 concentration intervals are obtained, wherein N is the total number of calibration samples; for the target concentration interval, based on the Spline interpolation algorithm, the function relationship between the concentration value and the reaction degree in the adjacent two sets of data corresponding to the target concentration interval is calculated, and the function relationship corresponding to the target concentration interval is obtained; wherein the target concentration interval is any one of the N-1 concentration intervals; according to the function relationship corresponding to the target concentration interval, the extreme point of the function relationship in the target concentration interval is calculated, if there is an extreme point in the target concentration interval, then according to the comparison result of the extreme point and the maximum concentration value and the minimum concentration value constituting the target concentration interval and the number of extreme points, it is judged whether there is a target calibration sample with incorrect reaction degree measurement in the adjacent two sets of data corresponding to the target concentration interval, if there is a target calibration sample, then the reaction degree of the target calibration sample is re-measured, the reaction degree obtained by re-measuring is used to replace the reaction degree of the target calibration sample, and the step of calculating the function relationship between the concentration value and the reaction degree of the target concentration interval based on the Spline interpolation algorithm is re-executed until there is no extreme point in the target concentration interval; if there is no extreme point in the target concentration interval, then the curve fitting data is obtained according to the function relationship of each concentration interval, whether there is a calibration sample with incorrect reaction degree measurement is judged according to the curve fitting data, if there is a calibration sample with incorrect reaction degree measurement, then the reaction degree of the incorrect calibration sample is re-measured, if there is no calibration sample with incorrect reaction degree measurement, then the judgment is ended, and when there is no calibration sample with incorrect reaction degree measurement, the standard curve is obtained based on the Spline interpolation algorithm according to the concentration of the calibration samples and the correct reaction degree. Based on the technical scheme, the function relationship between the concentration value and the reaction degree is obtained by using the Spline interpolation algorithm, the calibration sample with incorrect reaction degree measurement is judged based on the function relationship and the curve fitting data, so that the reaction degree of the calibration sample with incorrect measurement is re-measured until all the judgment conditions meet the preset conditions, the judgment is ended, the accuracy of the reaction degree measurement of the calibration sample is improved, the accuracy of the obtained standard curve is improved, and the accuracy of the calibration is improved.

[0072] In order to better realize the above method, the embodiment provides a calibration device, which refers to Figure 2 , Figure 2 The structural block diagram of the calibration device provided by the embodiment is shown in Figure 2 The device 20 comprises:

[0073] The interval division module 201 is configured to obtain the reaction degree of the calibrators, take the corresponding relationship between the concentration value and the reaction degree of the calibrators as a group of data, arrange each group of data according to the concentration, group the concentration values in each adjacent two groups of data after the arrangement into a concentration interval, and obtain N-1 concentration intervals, wherein N is the total number of the calibrators.

[0074] The function establishment module 202 is configured to, for a target concentration interval, calculate the function relationship between the concentration value and the reaction degree in the adjacent two groups of data corresponding to the target concentration interval based on a Spline interpolation algorithm, and obtain the function relationship corresponding to the target concentration interval; wherein the target concentration interval is any one of the N-1 concentration intervals.

[0075] The first judgment module 203 is configured to calculate the extreme points of the function relationship in the target concentration interval according to the function relationship corresponding to the target concentration interval, and if there is an extreme point in the target concentration interval, judge whether there is a target calibrator with a wrong reaction degree in the adjacent two groups of data corresponding to the target concentration interval according to the comparison result of the extreme point and the maximum concentration value and the minimum concentration value constituting the target concentration interval and the number of the extreme points, if there is the target calibrator, re-measure the reaction degree of the target calibrator, replace the reaction degree of the target calibrator with the re-measured reaction degree, and re-perform the step of arranging each group of data according to the concentration until there is no extreme point in the target concentration interval.

[0076] The second judgment module 204 is configured to, if there is no extreme point in the target concentration interval, obtain the curve fitting data according to the function relationship of each concentration interval, calculate the first-order difference array of the curve fitting data, calculate the second-order difference array according to the first-order difference array, judge whether there is a calibrator with a wrong reaction degree according to the second-order difference array, if there is a calibrator with a wrong reaction degree, re-measure the reaction degree of the calibrator with the wrong reaction degree, if there is no calibrator with a wrong reaction degree, end the judgment, and obtain the standard curve based on the Spline interpolation algorithm according to the concentration of the calibrators and the reaction degree without errors.

[0077] In a possible design, the interval division module 201 further includes a comparison and judgment module, which is configured to compare the reaction degrees in the i+1th group of data and the reaction degrees in the ith group of data, when each group of data is sorted in ascending order according to the concentration values of the calibrators, if the reaction degree in the i+1th group of data is less than the reaction degree in the ith group of data, it is determined that the reaction degree measurement in the i+1th group of data is erroneous, and the reaction degree of the corresponding calibrator in the i+1th group of data is re-measured; the reaction degree re-measured for the calibrator corresponding to the i+1th group of data replaces the original reaction degree in the i+1th group of data, and the step of comparing the reaction degree in the i+1th group of data and the reaction degree in the ith group of data is re-executed until the reaction degree in the i+1th group of data is greater than the reaction degree in the ith group of data; wherein the value of i is from 1 to N-1; when each group of data is sorted in descending order according to the concentration values of the calibrators, if the reaction degree in the i+1th group of data is greater than the reaction degree in the ith group of data, it is determined that the reaction degree measurement in the ith group of data is erroneous, and the reaction degree of the corresponding calibrator in the ith group of data is re-measured; the reaction degree re-measured for the calibrator corresponding to the ith group of data replaces the original reaction degree in the ith group of data, and the step of comparing the reaction degree in the i+1th group of data and the reaction degree in the ith group of data is re-executed until the reaction degree in the i+1th group of data is less than the reaction degree in the ith group of data.

[0078] In a possible design, the one-time judgment module 203 further includes a derivative calculation module, which is configured to derive the function relationship to obtain a derivative function of the function relationship, and calculate extreme points according to the derivative function.

[0079] In a possible design, the one-time judgment module 203 is specifically configured to: compare the extreme points with the maximum concentration value and the minimum concentration value, if a difference between the extreme points and the minimum concentration value is greater than a difference between the maximum concentration value and the extreme points, and the number of the extreme points is greater than a first threshold value, it is determined that the reaction degree measurement of the calibrator with the maximum concentration value in the adjacent two groups of data corresponding to the target concentration interval is erroneous; if the difference between the extreme points and the minimum concentration value is less than the difference between the maximum concentration value and the extreme points, and the number of the extreme points is greater than the first threshold value, it is determined that the reaction degree measurement of the calibrator with the minimum concentration value in the adjacent two groups of data corresponding to the target concentration interval is erroneous.

[0080] In a possible design, the secondary judgment module 204 is specifically configured to: calculate a product between a first second-order difference array and a second second-order difference array; the first second-order difference array refers to a second-order difference array at k+1, and the second second-order difference array refers to a second-order difference array at k, where k is a preset gradient value; if the product between the first second-order difference array and the second second-order difference array is not less than 0, it is determined that there is no calibration sample with incorrect reaction rate measurement, and the judgment is ended; if the product between the first second-order difference array and the second second-order difference array is less than 0, the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is calculated; if the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is less than a second threshold, it is determined that there is no calibration sample with incorrect reaction rate measurement, and the judgment is ended; if the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is greater than the second threshold, the absolute value of the difference between the concentration value in each concentration interval and k is calculated, and it is determined that the calibration sample corresponding to the minimum value in the absolute value of the difference between the concentration value and k has incorrect reaction rate measurement, the reaction rate of the calibration sample with incorrect reaction rate measurement is re-measured, the reaction rate of the calibration sample with incorrect reaction rate measurement re-measured is used to replace the original reaction rate in the data corresponding to the calibration sample with incorrect reaction rate measurement, the step of arranging each group of data according to the size of the concentration is re-executed until the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is less than the second threshold, or the product between the first second-order difference array and the second second-order difference array is not less than 0.

[0081] Based on the above device, by acquiring the reaction degree of the calibration sample, taking the corresponding relationship between the concentration value of the calibration sample and the reaction degree as a group of data, arranging each group of data according to the size of the concentration, grouping the concentration values in each adjacent two groups of data after arrangement into a concentration interval, N-1 concentration intervals are obtained, wherein N is the total number of calibration samples; for the target concentration interval, based on the Spline interpolation algorithm, the function relationship between the concentration value and the reaction degree in the adjacent two groups of data corresponding to the target concentration interval is calculated, and the function relationship corresponding to the target concentration interval is obtained; wherein the target concentration interval is any one of the N-1 concentration intervals; according to the function relationship corresponding to the target concentration interval, the extreme point of the function relationship in the target concentration interval is calculated, if there is an extreme point in the target concentration interval, then according to the comparison result of the extreme point and the maximum concentration value and the minimum concentration value constituting the target concentration interval and the number of extreme points, it is judged whether there is a target calibration sample with incorrect reaction degree measurement in the adjacent two groups of data corresponding to the target concentration interval, if there is a target calibration sample, then the reaction degree of the target calibration sample is re-measured, the reaction degree obtained by re-measuring is used to replace the reaction degree of the target calibration sample, and the step of calculating the function relationship between the concentration value and the reaction degree of the target concentration interval based on the Spline interpolation algorithm is re-executed until there is no extreme point in the target concentration interval; if there is no extreme point in the target concentration interval, then the curve fitting data is obtained according to the function relationship of each concentration interval, whether there is a calibration sample with incorrect reaction degree measurement is judged according to the curve fitting data, if there is a calibration sample with incorrect reaction degree measurement, then the reaction degree of the incorrect calibration sample is re-measured, if there is no calibration sample with incorrect reaction degree measurement, then the judgment is ended, and when there is no calibration sample with incorrect reaction degree measurement, the standard curve is obtained based on the Spline interpolation algorithm according to the calibration sample concentration and the correct reaction degree. Based on the technical scheme, the function relationship between the concentration value and the reaction degree is obtained by using the Spline interpolation algorithm, the calibration sample with incorrect reaction degree measurement is judged based on the function relationship and the curve fitting data, so that the reaction degree of the calibration sample with incorrect measurement is re-measured until all the judgment conditions meet the preset conditions, and the judgment is ended, which improves the accuracy of the reaction degree measurement of the calibration sample, thereby improving the accuracy of the obtained standard curve, and further improving the accuracy of the calibration.

[0082] Figure 3 The internal structure diagram of the computer device in one embodiment is shown. The computer device can be a terminal or a server. As shown in FIG. 1, the computer device includes a central processing unit (CPU), a memory, a storage device, a keyboard, a mouse, a display screen, and a network interface. Figure 3As shown, the computer device includes a processor, a memory and a network interface connected through a system bus. Among them, the memory includes a non-volatile storage medium and an internal memory. The non-volatile storage medium of the computer device stores an operating system, and can also store a computer program, which, when executed by the processor, can enable the processor to execute each step of the above method. The internal memory can also store a computer program, which, when executed by the processor, can enable the processor to execute each step of the above method. Those skilled in the art can understand that Figure 3 The structure shown in the figure is only a block diagram of part of the structure related to the scheme of the present application, and does not constitute a limitation on the computer device to which the scheme of the present application is applied. The specific computer device can include more or fewer components than those shown in the figure, or combine certain components, or have a different component arrangement.

[0083] In one embodiment, a computer device is provided, including a memory and a processor, the memory stores a computer program, the computer program is executed by the processor to enable the processor to execute all steps of the above method.

[0084] In one embodiment, a computer readable storage medium is provided, which stores a computer program, the computer program is executed by a processor to enable the processor to execute all steps of the above method.

[0085] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiments can be completed by a computer program instructing related hardware. The program can be stored in a non-volatile computer readable storage medium, and the program can include the processes of the above-mentioned embodiments when executed. Any reference to memory, storage, database or other medium used in each embodiment provided by the present application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration but not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDR SDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0086] Any combination of the technical features in the above embodiments can be made. For the sake of brevity, the foregoing description has not described all possible combinations of the technical features in the above embodiments. However, as long as the combination of the technical features does not contradict, it should be considered within the scope of the present disclosure.

[0087] The above embodiments only express several implementation manners of the present application, and the description is relatively specific and detailed, but it should not be understood as a limitation on the patent scope of the present application. It should be pointed out that, for ordinary skilled persons in the art, several modifications and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the patent protection scope of the present application should be subject to the appended claims.

Claims

1. A method of calibrating, characterized by, The method comprises: obtaining the reaction degree of the calibrators, taking the corresponding relationship between the concentration value and the reaction degree of the calibrators as a group of data, arranging each group of data according to the concentration, grouping the concentration values in each adjacent two groups of data after the arrangement into a concentration interval, and obtaining N-1 concentration intervals, wherein N is the total number of the calibrators; for a target concentration interval, based on a Spline interpolation algorithm, calculating the functional relationship between the concentration values and the reaction degrees in the adjacent two groups of data corresponding to the target concentration interval, and obtaining the functional relationship corresponding to the target concentration interval; wherein the target concentration interval is any one of the N-1 concentration intervals; according to the functional relationship corresponding to the target concentration interval, calculating the extreme points of the functional relationship in the target concentration interval, if there is an extreme point in the target concentration interval, then according to the comparison result of the extreme point and the maximum concentration value and the minimum concentration value constituting the target concentration interval and the number of extreme points, judging whether there is a target calibrator with a measurement error of the reaction degree in the adjacent two groups of data corresponding to the target concentration interval, if there is the target calibrator, then re-measuring the reaction degree of the target calibrator, replacing the reaction degree of the target calibrator with the re-measured reaction degree, and re-executing the step of arranging each group of data according to the concentration until there is no extreme point in the target concentration interval; if there is no extreme point in the target concentration interval, then according to the functional relationship of each concentration interval, obtaining curve fitting data, calculating a first-order difference array of the curve fitting data, calculating a second-order difference array according to the first-order difference array, judging whether there is a calibrator with a measurement error of the reaction degree according to the second-order difference array, if there is a calibrator with a measurement error of the reaction degree, then re-measuring the reaction degree of the calibrator with the measurement error, if there is no calibrator with a measurement error of the reaction degree, then ending the judgment, and based on the Spline interpolation algorithm, obtaining a standard curve according to the concentration of the calibrators and the measurement error-free reaction degree; wherein the judgment of whether there is a target calibrator with a measurement error of the reaction degree in the adjacent two groups of data corresponding to the target concentration interval according to the comparison result of the extreme point and the maximum concentration value and the minimum concentration value constituting the target concentration interval and the number of extreme points comprises: comparing the extreme point with the maximum concentration value and the minimum concentration value, if the difference between the extreme point and the minimum concentration value is greater than the difference between the maximum concentration value and the extreme point, and the number of extreme points is greater than a first threshold, then it is determined that the reaction degree of the calibrator with the maximum concentration value in the adjacent two groups of data corresponding to the target concentration interval is measured with an error; if the difference between the extreme point and the minimum concentration value is less than the difference between the maximum concentration value and the extreme point, and the number of extreme points is greater than the first threshold, then it is determined that the reaction degree of the calibrator with the minimum concentration value in the adjacent two groups of data corresponding to the target concentration interval is measured with an error. If the second-order difference array indicates that there is a calibration sample with a measurement error, the reactivity of the calibration sample with the measurement error is re-measured; if there is no calibration sample with a measurement error, the judgment is ended, comprising: calculating the product between the first second-order difference array and the second second-order difference array; wherein the first second-order difference array refers to the second-order difference array at k+1, the second second-order difference array refers to the second-order difference array at k, and k is a preset gradient value; if the product between the first second-order difference array and the second second-order difference array is not less than 0, it is determined that there is no calibration sample with a measurement error, and the judgment is ended; if the product between the first second-order difference array and the second second-order difference array is less than 0, the absolute value of the first second-order difference array and the absolute value of the second second-order difference array are calculated, if the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is less than a second threshold value, it is determined that there is no calibration sample with a measurement error, and the judgment is ended; if the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is greater than the second threshold value, the absolute value of the difference between the concentration value in each concentration interval and k is calculated, and it is judged that the reactivity of the calibration sample corresponding to the minimum value in the absolute value of the difference between the concentration value and k is measured with an error, the reactivity of the calibration sample measured with an error is re-measured, the reactivity of the calibration sample measured with an error is replaced by the original reactivity in the data corresponding to the calibration sample, and the step of arranging each group of data according to the size of the concentration is re-executed until the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is less than the second threshold value, or the product between the first second-order difference array and the second second-order difference array is not less than 0.

2. The method of claim 1, wherein, Before the concentration values in every adjacent two groups of arranged data are grouped into a concentration interval to obtain N-1 concentration intervals, comprising: comparing the reactivity in the i+1th group of data and the reactivity in the ith group of data, when each group of data is sorted from small to large according to the concentration value of the calibration sample, if the reactivity in the i+1th group of data is less than the reactivity in the ith group of data, it is determined that the reactivity in the i+1th group of data has a measurement error, and the reactivity of the corresponding calibration sample in the i+1th group of data is re-measured; the reactivity of the calibration sample corresponding to the i+1th group of data re-measured is replaced by the original reactivity in the i+1th group of data, and the step of comparing the reactivity in the i+1th group of data and the reactivity in the ith group of data is re-executed until the reactivity in the i+1th group of data is greater than the reactivity in the ith group of data; wherein the value of i is from 1 to N-1; When each group of data is ordered from large to small according to the concentration values of the calibrators, if the reaction degree in the i+1th group of data is greater than the reaction degree in the ith group of data, it is determined that the reaction degree measurement in the ith group of data has an error, and the reaction degree of the corresponding calibrator in the ith group of data is re-measured; the reaction degree re-measured for the calibrator corresponding to the ith group of data replaces the original reaction degree in the ith group of data, and the step of comparing the reaction degree in the i+1th group of data with the reaction degree in the ith group of data is re-executed until the reaction degree in the i+1th group of data is less than the reaction degree in the ith group of data.

3. The method of claim 1, wherein, The function relationship corresponding to the target concentration interval is calculated, and an extreme point of the function relationship in the target concentration interval is obtained, including: The derivative of the function relationship is obtained, and the extreme point is calculated according to the derivative.

4. A calibration device, characterized by The device comprises: An interval division module is configured to obtain the reaction degrees of the calibrators, take the corresponding relationship between the concentration values and the reaction degrees of the calibrators as a group of data, arrange each group of data according to the concentration, group the concentration values in each adjacent two groups of data arranged, and obtain N-1 concentration intervals, where N is the total number of the calibrators; A function establishment module is configured to, for a target concentration interval, calculate a function relationship between the concentration values and the reaction degrees in the adjacent two groups of data corresponding to the target concentration interval based on a Spline interpolation algorithm, and obtain the function relationship corresponding to the target concentration interval; the target concentration interval is any one of the N-1 concentration intervals; A first judgment module is configured to calculate an extreme point of the function relationship in the target concentration interval according to the function relationship corresponding to the target concentration interval, and if the extreme point exists in the target concentration interval, judge whether there is a target calibrator with an erroneous reaction degree measurement in the adjacent two groups of data corresponding to the target concentration interval according to the comparison result of the extreme point and the maximum and minimum concentration values constituting the target concentration interval and the number of the extreme points, re-measure the reaction degree of the target calibrator if the target calibrator exists, replace the reaction degree of the target calibrator with the re-measured reaction degree, and re-execute the step of arranging each group of data according to the concentration until the target concentration interval does not have the extreme point; A second judgment module is configured to, if the target concentration interval does not have the extreme point, obtain curve fitting data according to the function relationship of each concentration interval, calculate a first-order difference array of the curve fitting data, calculate a second-order difference array according to the first-order difference array, judge whether there is a calibrator with an erroneous reaction degree measurement according to the second-order difference array, re-measure the reaction degree of the calibrator with the erroneous reaction degree measurement if the calibrator with the erroneous reaction degree measurement exists, and end the judgment if the calibrator with the erroneous reaction degree measurement does not exist, and obtain a standard curve based on the Spline interpolation algorithm according to the concentration of the calibrators and the reaction degrees without errors. The first determining module is specifically configured to compare the extreme value points with the maximum concentration value and the minimum concentration value, and if a difference between the extreme value points and the minimum concentration value is greater than a difference between the maximum concentration value and the extreme value points, and the number of the extreme value points is greater than the first threshold value, it is determined that the reaction degree measurement of the calibrant with the maximum concentration value in the adjacent two groups of data corresponding to the target concentration interval is incorrect; if the difference between the extreme value points and the minimum concentration value is less than the difference between the maximum concentration value and the extreme value points, and the number of the extreme value points is greater than the first threshold value, it is determined that the reaction degree measurement of the calibrant with the minimum concentration value in the adjacent two groups of data corresponding to the target concentration interval is incorrect; The second determining module is specifically configured to calculate a product between a first second-order difference array and a second second-order difference array; wherein the first second-order difference array refers to a second-order difference array under k+1, and the second second-order difference array refers to a second-order difference array under k, and k is a preset gradient value; if the product between the first second-order difference array and the second second-order difference array is not less than 0, it is determined that there is no calibrant with incorrect reaction degree measurement, and the judgment is ended; if the product between the first second-order difference array and the second second-order difference array is less than 0, the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is calculated, if the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is less than a second threshold value, it is determined that there is no calibrant with incorrect reaction degree measurement, and the judgment is ended; if the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is greater than the second threshold value, the absolute value of the difference between the concentration value in each concentration interval and k is calculated, and it is determined that the reaction degree measurement of the calibrant corresponding to the minimum value in the absolute value of the difference between the concentration value and k is incorrect, the reaction degree of the calibrant with incorrect reaction degree measurement is re-measured, the reaction degree of the calibrant with incorrect reaction degree measurement re-measured is replaced with the original reaction degree in the data corresponding to the calibrant with incorrect reaction degree measurement, and the step of arranging each group of data according to the size of the concentration is re-executed until the sum of the absolute value of the first second-order difference array and the absolute value of the second second-order difference array is less than the second threshold value, or the product between the first second-order difference array and the second second-order difference array is not less than 0.

5. The apparatus of claim 4, wherein, The interval dividing module comprises: The comparison and judgment module is configured to compare the reaction degree in the i+1th group of data with the reaction degree in the ith group of data, when each group of data is sorted in ascending order according to the concentration values of the calibration samples, if the reaction degree in the i+1th group of data is less than the reaction degree in the ith group of data, it is determined that the reaction degree measurement in the i+1th group of data has an error, and the reaction degree of the corresponding calibration sample in the i+1th group of data is re-measured; the reaction degree of the corresponding calibration sample in the i+1th group of data is re-measured, and the original reaction degree in the i+1th group of data is replaced, and the step of comparing the reaction degree in the i+1th group of data with the reaction degree in the ith group of data is re-executed until the reaction degree in the i+1th group of data is greater than the reaction degree in the ith group of data; wherein, the value of i is from 1 to N-1; when each group of data is sorted in descending order according to the concentration values of the calibration samples, if the reaction degree in the i+1th group of data is greater than the reaction degree in the ith group of data, it is determined that the reaction degree measurement in the ith group of data has an error, and the reaction degree of the corresponding calibration sample in the ith group of data is re-measured; the reaction degree of the corresponding calibration sample in the ith group of data is re-measured, and the original reaction degree in the ith group of data is replaced, and the step of comparing the reaction degree in the i+1th group of data with the reaction degree in the ith group of data is re-executed until the reaction degree in the i+1th group of data is less than the reaction degree in the ith group of data.

6. The apparatus of claim 5, wherein, The first judgment module comprises: The derivative calculation module is configured to derive the function relationship to obtain a derivative function of the function relationship, and calculate an extreme point according to the derivative function.

7. A computer readable storage medium storing a computer program, wherein the computer program comprises program instructions configured to cause a processor to perform the method according to any one of claims 1 to 6. The computer program is executed by the processor to make the processor execute the steps of the method of any one of claims 1 to 3. 8.A computer device, comprising a memory and a processor, and characterized in that, The computer program is executed by the processor to make the processor execute the steps of the method of any one of claims 1 to 3.

Images