WO2022100236A1 - Method for fluorescence intensity peak detection - Google Patents
Method for fluorescence intensity peak detection Download PDFInfo
- Publication number
- WO2022100236A1 WO2022100236A1 PCT/CN2021/116363 CN2021116363W WO2022100236A1 WO 2022100236 A1 WO2022100236 A1 WO 2022100236A1 CN 2021116363 W CN2021116363 W CN 2021116363W WO 2022100236 A1 WO2022100236 A1 WO 2022100236A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- value
- fluorescence intensity
- slope
- peak
- data
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000001514 detection method Methods 0.000 title claims abstract description 25
- 238000003018 immunoassay Methods 0.000 claims abstract description 11
- 238000012360 testing method Methods 0.000 claims description 18
- 230000035772 mutation Effects 0.000 claims description 14
- 238000001914 filtration Methods 0.000 claims description 8
- 239000003153 chemical reaction reagent Substances 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract description 16
- 230000008859 change Effects 0.000 abstract description 9
- 230000000694 effects Effects 0.000 abstract description 3
- 238000012545 processing Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 102000004190 Enzymes Human genes 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 229940088597 hormone Drugs 0.000 description 1
- 239000005556 hormone Substances 0.000 description 1
- 238000003317 immunochromatography Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 239000012488 sample solution Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/58—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving labelled substances
- G01N33/582—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving labelled substances with fluorescent label
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
- G01N21/6428—Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes"
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/53—Immunoassay; Biospecific binding assay; Materials therefor
- G01N33/558—Immunoassay; Biospecific binding assay; Materials therefor using diffusion or migration of antigen or antibody
-
- 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
Definitions
- the invention relates to the technical field of fluorescence immunodetection, and more particularly, to a method for fluorescence intensity peak detection.
- Fluorescence immunoassay technology has the advantages of strong specificity, high sensitivity and good practicability, and can be used to detect biologically active compounds with very low content, such as proteins (enzymes, receptors, antibodies), hormones, drugs and microorganisms.
- the working process of fluorescence immunoassay by the fluorescence immunoassay analyzer is as follows: the reagent card passes through the detection area, the LED light source generates excitation fluorescence, and gathers on the target detection object of the reagent card. After the target detection object is excited, fluorescence is generated, which is detected by the analyzer. , to obtain a fluorescence intensity data curve, in which there will be a peak corresponding to the contained target. The appearance position of the peak depends on the type of the component, and the size of the peak (ie, its height or area) depends on the amount or concentration of the component corresponding to the peak.
- the wave peak value is generally determined by the wave peak value calculation method, which directly selects the maximum fluorescence intensity value as the wave peak value; or selects several fluorescence intensity values with relatively large values, and takes their average value as the wave peak value. The position corresponding to the peak value of the obtained wave is then the position of the wave peak.
- the fluorescence intensity value may change abruptly, so that the fluorescence intensity data curve has some maximum or larger values similar to the wave peak, but these maximum values are different from those to be measured. There is no direct correspondence between the concentrations of the substances. Using the above wave peak value calculation method, the maximum or larger value of this sudden change will also be included in the calculation, which will lead to deviations in the wave peak value calculation, resulting in a decrease in the accuracy of the final detection concentration.
- the Chinese Patent Publication No. CN106645708A provides a quantitative detection and calculation method based on fluorescence immunochromatography technology, which includes the following steps: inserting the fluorescence immunoassay test strip with the sample solution to be tested dropwise into the fluorescence immunoassay analyzer, and sequentially passing through the LED light source Irradiation, filter processing, photodetector detection, electrical signal processing, AD conversion processing, filtering algorithm processing, baseline fitting, wave peak processing, calculation of T/C area ratio, calculation of concentration process, so as to pass the calculated Concentration results are judged.
- the above calculation method filters out all kinds of noise and signal interference through filtering algorithm processing and baseline fitting, and provides more accurate detection results, but there is still a possibility that the abnormal mutation signal cannot be excluded, so the abnormal mutation signal is included in the calculation of the peak value and affects the detection result. Condition.
- the purpose of the present invention is to provide a method for detecting the fluorescence intensity peak in order to make up for the deficiencies of the prior art.
- a method for fluorescence intensity peak detection comprising the following steps:
- Peak detection find out the peak position from the fluorescence intensity scattergram, including the following steps:
- N Determine whether the number of data N is greater than M: if N ⁇ M, exclude the data points and endpoints between i min and i max in the slope array obtained by the closest step, and the remaining slope values form a new , repeat the above steps bd in the new slope array until N>M is satisfied; if N>M, the corresponding output numerical sequence number in the slope array obtained in the closest step is between i min and i max Find the value with the smallest absolute value of the slope value among the slope values between the two, and determine the output value number i corresponding to the value with the smallest absolute value of the slope value.
- the fluorescence intensity value corresponding to the output value number i is the peak value, and its position is the peak value. Location.
- the slope value K i is used to represent the data change trend between the fluorescence intensity point corresponding to the output numerical number i and the two adjacent fluorescence intensity points; the slope value K i can be calculated according to the following formula :
- any positive integer value of m can be selected according to actual needs.
- the theoretically calculated slope The value should be more accurate; however, it is sufficient to select 1 or 2 for m in actual operation, so m is preferably 1 or 2 in this case.
- the A is preferably 3 or 4.
- the fluorescence immunoassay analyzer in step (1) adopts ADUCM360 as the main chip for data acquisition, the set test frequency f is 500Hz, and the test time t is 2s.
- the ADUCM360 is a low-power precision analog microcontroller integrating dual-channel ⁇ - ⁇ ADC and ARMCortex-M3.
- the data register of the AD converter has a total of 24 bits, and the system selects 16-bit valid data, which is conducive to improving accuracy of the test.
- the test card is driven by the step frequency of 500Hz, and the data acquisition module synchronously collects the fluorescence intensity data through the frequency of 500Hz, and obtains 1000 fluorescence intensity data sequentially in a total of 2S time.
- the step (2) before the step (2), it further includes performing a smoothing filtering process on the fluorescence intensity points in the step (1), so as to obtain a processed fluorescence intensity scattergram; the method for the smoothing filtering processing adopts the prior art , which will not be described in detail here, the data is smoothed by the smoothing filtering process, and the noise interference is reduced.
- the present invention calculates the slope value of the fluorescence intensity point by setting the method for fluorescence intensity peak detection, and adopts the slope value calculation method to calculate the data number between the maximum value and the minimum value of the slope value.
- the detection of the fluorescence intensity peak can improve the accuracy of the determination of the peak value and its position, and effectively solve the existing method for determining the peak value.
- the occurrence of the mutation point affects the determination of the peak value and its position. The accuracy of the technical problem, especially when the mutation point is too large, the effect is more significant.
- FIG. 1 is a schematic flowchart of the method for detecting fluorescence intensity peaks according to the present invention
- Fig. 2 is the fluorescence intensity scatter diagram obtained after smooth filtering processing in the present invention
- FIG. 3 is a scatter diagram corresponding to the output numerical serial number i and the slope K i according to the present invention.
- this embodiment discloses a method for fluorescence intensity peak detection, which includes the following steps:
- the fluorescent immune analyzer adopts ADUCM360 as the main chip for data acquisition, the set test frequency f is 500Hz, and the test time t is 2s, so that the output numerical serial number i is the abscissa and the fluorescence intensity value. 1000 fluorescence intensity points are obtained on the coordinate system where F is the ordinate, and the corresponding fluorescence intensity scattergram is obtained; the fluorescence intensity points are subjected to smooth filtering to obtain the processed fluorescence intensity scattergram; Figure see Figure 2;
- Peak detection find out the peak position from the fluorescence intensity scattergram, including the following steps:
- +1
- This example describes in detail how to specifically determine the correct wave peak and wave peak position using the method for fluorescence intensity peak detection described in this application when there is only one mutation point in the fluorescence intensity scattergram obtained from the test;
- the method for detecting fluorescence intensity peaks if there are 2 or more mutation points, multiple mutation points can be eliminated through multiple cycles of steps b-d, which will not be described in detail here.
- the existing wave peak value is determined by directly selecting the maximum value of the fluorescence intensity value as the wave peak value, or selecting several fluorescence intensity values with relatively large values, and taking their average value as the wave peak value. It can be seen that, if the existing wave peak determination method is used to determine the wave peak value in the fluorescence intensity scattergram shown in Figure 2, since the fluorescence intensity value corresponding to the abrupt change point is larger or even greater than the wave peak value, the sudden change point will be It will be included in the calculation, resulting in a large wave peak value, which leads to the inability to accurately determine the wave peak value and its corresponding position, which affects the detection accuracy.
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Immunology (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Urology & Nephrology (AREA)
- Hematology (AREA)
- Pathology (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Cell Biology (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Microbiology (AREA)
- Mathematical Analysis (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Biotechnology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Optics & Photonics (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
The present invention relates to the technical field of fluorescence immunoassay, and more particularly to a method for fluorescence intensity peak detection. In the present invention, a method for performing fluorescence intensity peak detection is provided, the method comprising: calculating slope values of fluorescence intensity points using a slope value calculation method, so as to determine the number of pieces of data between the maximum slope value and the minimum slope value, and comparing the number of pieces of data with a set value, so as to exclude change points. The fluorescence intensity peak detection can improve the accuracy of the determination of peak values and the positions thereof, and can effectively solve the technical problem of the existing peak value determination method affecting the accuracy of the determination of peak values and the positions thereof due to the presence of change points, and has a more remarkable effect especially when the change points are relatively large.
Description
本发明涉及荧光免疫检测技术领域，更具体地说，涉及一种用于荧光强度峰检测的方法。The invention relates to the technical field of fluorescence immunodetection, and more particularly, to a method for fluorescence intensity peak detection.
荧光免疫检测技术具有专一性强、灵敏度高、实用性好等优点，可被用于检测含量很低的生物活性化合物，如蛋白质（酶、接受体、抗体）、激素、药物及微生物等。荧光免疫分析仪进行荧光免疫检测的工作过程是：试剂卡通过检测区，LED光源产生激发荧光，并聚集在试剂卡的目标检测物上，目标检测物受激发后产生荧光，由分析仪检测获取，从而获得的荧光强度数据曲线，曲线中会出现与所含的目标物相对应的波峰。所述波峰的出现位置取决于成分的种类，波峰的大小(即其高度或面积)取决于与所述波峰相对应的成分的量或浓度。Fluorescence immunoassay technology has the advantages of strong specificity, high sensitivity and good practicability, and can be used to detect biologically active compounds with very low content, such as proteins (enzymes, receptors, antibodies), hormones, drugs and microorganisms. The working process of fluorescence immunoassay by the fluorescence immunoassay analyzer is as follows: the reagent card passes through the detection area, the LED light source generates excitation fluorescence, and gathers on the target detection object of the reagent card. After the target detection object is excited, fluorescence is generated, which is detected by the analyzer. , to obtain a fluorescence intensity data curve, in which there will be a peak corresponding to the contained target. The appearance position of the peak depends on the type of the component, and the size of the peak (ie, its height or area) depends on the amount or concentration of the component corresponding to the peak.
目前一般采用波峰值计算方法确定波峰的大小，其采用直接选取荧光强度值最大值作为波峰值；或者选择若干个数值相对较大的荧光强度值，取它们的平均值作为波峰值。所得波峰值对应的位置则为波峰的位置。然而，在荧光测量过程中，由于一些干扰因素的存在，荧光强度值可能会出现突变的情况，使得荧光强度数据曲线出现一些与波峰类似的最大值或较大值，但这些最大值与待测物浓度并没有直接对应关系。采用上述波峰值计算方法，会将这种突变的最大值或较大值也纳入计算，这样便会导致波峰值计算出现偏差，导致最终的检测浓度准确度降低。At present, the wave peak value is generally determined by the wave peak value calculation method, which directly selects the maximum fluorescence intensity value as the wave peak value; or selects several fluorescence intensity values with relatively large values, and takes their average value as the wave peak value. The position corresponding to the peak value of the obtained wave is then the position of the wave peak. However, in the process of fluorescence measurement, due to the existence of some interference factors, the fluorescence intensity value may change abruptly, so that the fluorescence intensity data curve has some maximum or larger values similar to the wave peak, but these maximum values are different from those to be measured. There is no direct correspondence between the concentrations of the substances. Using the above wave peak value calculation method, the maximum or larger value of this sudden change will also be included in the calculation, which will lead to deviations in the wave peak value calculation, resulting in a decrease in the accuracy of the final detection concentration.
中国专利公开号为CN106645708A提供一种基于荧光免疫层析技术的定量检测计算方法，包括以下步骤：将滴加有待测样本溶液的荧光免疫试纸条插入荧光免疫分析仪中，依次经过LED光源照射、滤光片滤光处理、光电探测器检测、电信号处理、AD转换处理、滤波算法处理、基线拟合、寻波峰处理、计算T/C面积比值、计算浓度的过程，从而通过计算的浓度结果进行判断。上述计算方法通过滤波算法处理及基线拟合过滤掉各种噪声、信号干扰，提供较为精确的检测结果，但是仍存在无法排除突变异常信号，从而将突变异常信号纳入波峰值计算从而影响检测结果的情况。The Chinese Patent Publication No. CN106645708A provides a quantitative detection and calculation method based on fluorescence immunochromatography technology, which includes the following steps: inserting the fluorescence immunoassay test strip with the sample solution to be tested dropwise into the fluorescence immunoassay analyzer, and sequentially passing through the LED light source Irradiation, filter processing, photodetector detection, electrical signal processing, AD conversion processing, filtering algorithm processing, baseline fitting, wave peak processing, calculation of T/C area ratio, calculation of concentration process, so as to pass the calculated Concentration results are judged. The above calculation method filters out all kinds of noise and signal interference through filtering algorithm processing and baseline fitting, and provides more accurate detection results, but there is still a possibility that the abnormal mutation signal cannot be excluded, so the abnormal mutation signal is included in the calculation of the peak value and affects the detection result. Condition.
故，现有技术具有较大的改进空间。Therefore, the prior art has a large room for improvement.
本发明的目的是为了弥补现有技术的不足，提出一种用于荧光强度峰检测的方法。The purpose of the present invention is to provide a method for detecting the fluorescence intensity peak in order to make up for the deficiencies of the prior art.
为了达到上述目的，本发明通过以下技术方案实现：In order to achieve the above object, the present invention realizes through the following technical solutions:
一种用于荧光强度峰检测的方法，包括以下步骤：A method for fluorescence intensity peak detection, comprising the following steps:
（1）采用荧光免疫分析仪以设定的测试频率f、测试时间t对试剂卡进行测试，按序输出荧光强度值，从而在以输出数值序号i为横坐标、荧光强度值F为纵坐标的坐标系上得到（f×t）个荧光强度点，得到对应的荧光强度散点图；(1) Use a fluorescence immunoassay analyzer to test the reagent card with the set test frequency f and test time t, and output the fluorescence intensity values in sequence, so that the output numerical serial number i is the abscissa and the fluorescence intensity value F is the ordinate. (f×t) fluorescence intensity points are obtained on the coordinate system of , and the corresponding fluorescence intensity scattergram is obtained;
（2）峰检测：从荧光强度散点图中找出波峰位置，包括以下步骤：(2) Peak detection: find out the peak position from the fluorescence intensity scattergram, including the following steps:
a.计算荧光强度点的斜率值K
i，将斜率值以其对应的输出数值序号i从小到大作为排列顺序进行排序得到对应的斜率数组；设定斜率数组中斜率值最大值与最小值间数据个数的最小阈值M ；
a. Calculate the slope value K i of the fluorescence intensity point, sort the slope value with its corresponding output value number i from small to large as the arrangement order to obtain the corresponding slope array; set the slope value between the maximum value and the minimum value in the slope array The minimum threshold M of the number of data;
b.遍历上一步骤所得的斜率数组找出斜率值的最大值和最小值，根据斜率值的最大值、最小值找到其对应的荧光强度数据点对应的输出数值序号，分别记为i
min、i
max；
b. Traverse the slope array obtained in the previous step to find the maximum value and the minimum value of the slope value, and find the output numerical sequence number corresponding to the corresponding fluorescence intensity data point according to the maximum value and the minimum value of the slope value, and denote them as i min , respectively. i max ;
c.计算i
min、i
max之间的数据个数N=|i
max-i
min|+1；
c. Calculate the number of data between i min and i max N=|i max -i min |+1;
d.判断所述数据个数N是否大于M：若N≤M，则在最接近的一个步骤所得斜率数组中将i
min、i
max之间的数据点及端点排除，剩余的斜率值组成新的斜率数组，在新的斜率数组中重复进行上述步骤b-d直至满足N＞M；若N＞M，则在最接近的一个步骤所得斜率数组中其对应的输出数值序号在i
min至i
max之间的斜率值中找到斜率值绝对值最小的值，确定该斜率值绝对值最小的值对应的输出数值序号i，该输出数值序号i对应的荧光强度值则为波峰值，其位置则为波峰位置。
d. Determine whether the number of data N is greater than M: if N≤M, exclude the data points and endpoints between i min and i max in the slope array obtained by the closest step, and the remaining slope values form a new , repeat the above steps bd in the new slope array until N>M is satisfied; if N>M, the corresponding output numerical sequence number in the slope array obtained in the closest step is between i min and i max Find the value with the smallest absolute value of the slope value among the slope values between the two, and determine the output value number i corresponding to the value with the smallest absolute value of the slope value. The fluorescence intensity value corresponding to the output value number i is the peak value, and its position is the peak value. Location.
根据以上方案，所述斜率值K
i用于表示输出数值序号i对应的荧光强度点与前后相邻的两个荧光强度点之间的数据变化趋势；所述斜率值K
i可按以下公式计算：
According to the above scheme, the slope value K i is used to represent the data change trend between the fluorescence intensity point corresponding to the output numerical number i and the two adjacent fluorescence intensity points; the slope value K i can be calculated according to the following formula :
根据以上方案，所述m可根据实际需要选取任意正整数值，其取值越大，在计算斜率值公式时，选取的连续数据点越多，计算工作量也越大，理论上计算的斜率值应该越准确；然而在现实操作时所述m选取1或2已足够，因此本案所述m优选为1或2。According to the above scheme, any positive integer value of m can be selected according to actual needs. The larger the value is, the more continuous data points are selected when the slope value formula is calculated, and the calculation workload is also larger. The theoretically calculated slope The value should be more accurate; however, it is sufficient to select 1 or 2 for m in actual operation, so m is preferably 1 or 2 in this case.
根据以上方案，所述最小阈值M为斜率数组中斜率值最大值与最小值间数据个数，本领域技术人员可以根据经验进行设定；也可以按M= A×（2m+1）计算得出；其中，A=希望排除突变点的个数，可以根据仪器、试剂等硬件的情况进行确定，一般来说A取值为1、2、3、4、5的数值。如果突变点过多，则需要先验证硬件是否故障。According to the above scheme, the minimum threshold M is the number of data between the maximum value and the minimum value of the slope value in the slope array, which can be set by those skilled in the art according to experience; it can also be calculated according to M=A×(2m+1) Among them, A = the number of mutation points that you want to exclude, which can be determined according to the hardware conditions such as instruments and reagents. Generally speaking, A is a value of 1, 2, 3, 4, and 5. If there are too many mutation points, you need to verify whether the hardware is faulty first.
根据以上方案，所述A优选为3或4。According to the above scheme, the A is preferably 3 or 4.
根据以上方案，步骤（1）中所述荧光免疫分析仪采用ADUCM360作为数据采集主芯片，所述设定的测试频率f为500Hz，所述测试时间t为2s。所述ADUCM360是一种集成双通道Σ-Δ型ADC和ARMCortex-M3的低功耗精密模拟微控制器，其AD转换器的数据寄存器总共有24位，系统选用16位有效数据，有利于提高测试的精确度。通过500Hz的步进频率带动测试卡，数据采集模块同步通过500Hz的频率采集荧光强度数据，总共2S时间按序获得1000个荧光强度数据。According to the above scheme, the fluorescence immunoassay analyzer in step (1) adopts ADUCM360 as the main chip for data acquisition, the set test frequency f is 500Hz, and the test time t is 2s. The ADUCM360 is a low-power precision analog microcontroller integrating dual-channel Σ-Δ ADC and ARMCortex-M3. The data register of the AD converter has a total of 24 bits, and the system selects 16-bit valid data, which is conducive to improving accuracy of the test. The test card is driven by the step frequency of 500Hz, and the data acquisition module synchronously collects the fluorescence intensity data through the frequency of 500Hz, and obtains 1000 fluorescence intensity data sequentially in a total of 2S time.
根据以上方案，所述步骤（2）前还包括将步骤（1）所述荧光强度点进行平滑滤波处理，从而得到处理后的荧光强度散点图；所述平滑滤波处理的方法采用现有技术，在此不再详述，通过平滑滤波处理对数据起到平滑作用，减小噪音干扰。According to the above solution, before the step (2), it further includes performing a smoothing filtering process on the fluorescence intensity points in the step (1), so as to obtain a processed fluorescence intensity scattergram; the method for the smoothing filtering processing adopts the prior art , which will not be described in detail here, the data is smoothed by the smoothing filtering process, and the noise interference is reduced.
本发明的有益效果在于：The beneficial effects of the present invention are:
本发明通过设置用于荧光强度峰检测的方法，采用斜率值计算方法对荧光强度点的斜率值进行计算，确定斜率值最大值、最小值之间的数据个数，通过数据个数与设定值进行大小比较，从而能够排除突变点，所述荧光强度峰检测能够提高波峰值及其位置确定的准确性，有效解决现有波峰值的确定方法由于突变点的出现影响波峰值及其位置确定的准确性的技术问题，尤其在突变点偏大时其效果更为显著。The present invention calculates the slope value of the fluorescence intensity point by setting the method for fluorescence intensity peak detection, and adopts the slope value calculation method to calculate the data number between the maximum value and the minimum value of the slope value. The detection of the fluorescence intensity peak can improve the accuracy of the determination of the peak value and its position, and effectively solve the existing method for determining the peak value. The occurrence of the mutation point affects the determination of the peak value and its position. The accuracy of the technical problem, especially when the mutation point is too large, the effect is more significant.
图1为本发明所述用于荧光强度峰检测的方法的流程示意图；1 is a schematic flowchart of the method for detecting fluorescence intensity peaks according to the present invention;
图2为本发明中平滑滤波处理后所得的荧光强度散点图；Fig. 2 is the fluorescence intensity scatter diagram obtained after smooth filtering processing in the present invention;
图3为本发明所述输出数值序号i与斜率K
i对应的散点图。
FIG. 3 is a scatter diagram corresponding to the output numerical serial number i and the slope K i according to the present invention.
为了更好地理解本发明，下面结合实施例进一步阐明本发明的内容，但本发明不仅仅局限于下面的实施例。In order to better understand the present invention, the content of the present invention is further illustrated below in conjunction with the examples, but the present invention is not limited to the following examples.
如图1所示，本实施例公开了一种用于荧光强度峰检测的方法，包括以下步骤：As shown in FIG. 1 , this embodiment discloses a method for fluorescence intensity peak detection, which includes the following steps:
（1）采用荧光免疫分析仪以设定的测试频率f、测试时间t对试剂卡进行测试，按序输出荧光强度值，从而在以输出数值序号i为横坐标、荧光强度值F为纵坐标的坐标系上得到（f×t）个荧光强度点，得到对应的荧光强度散点图；(1) Use a fluorescence immunoassay analyzer to test the reagent card with the set test frequency f and test time t, and output the fluorescence intensity values in sequence, so that the output numerical serial number i is the abscissa and the fluorescence intensity value F is the ordinate. (f×t) fluorescence intensity points are obtained on the coordinate system of , and the corresponding fluorescence intensity scattergram is obtained;
进一步地，所述荧光免疫分析仪采用ADUCM360作为数据采集主芯片，所述设定的测试频率f为500Hz，所述测试时间t为2s，从而在以输出数值序号i为横坐标、荧光强度值F为纵坐标的坐标系上得到1000个荧光强度点，得到对应的荧光强度散点图；将所述荧光强度点进行平滑滤波处理，从而得到处理后的荧光强度散点图；具体的散点图见图2；Further, the fluorescent immune analyzer adopts ADUCM360 as the main chip for data acquisition, the set test frequency f is 500Hz, and the test time t is 2s, so that the output numerical serial number i is the abscissa and the fluorescence intensity value. 1000 fluorescence intensity points are obtained on the coordinate system where F is the ordinate, and the corresponding fluorescence intensity scattergram is obtained; the fluorescence intensity points are subjected to smooth filtering to obtain the processed fluorescence intensity scattergram; Figure see Figure 2;
（2）峰检测：从荧光强度散点图中找出波峰位置，包括以下步骤：(2) Peak detection: find out the peak position from the fluorescence intensity scattergram, including the following steps:
a.计算荧光强度点的斜率K
i，将斜率值以其对应的输出数值序号i从小到大作为排列顺序进行排序得到含有998个斜率值的斜率数组；所述斜率值K
i按以下公式计算：
a. Calculate the slope K i of the fluorescence intensity point, sort the slope values with their corresponding output numerical serial numbers i from small to large as the arrangement order to obtain a slope array containing 998 slope values; the slope value K i is calculated according to the following formula :
b.遍历上一步骤所得的斜率数组找出斜率值的最大值和最小值，根据斜率值的最大值、最小值找到其对应的荧光强度数据点对应的输出数值序号，分别记为i
min、i
max；结合图3可见，i
min=501,i
max=499;
b. Traverse the slope array obtained in the previous step to find the maximum value and the minimum value of the slope value, and find the output numerical sequence number corresponding to the corresponding fluorescence intensity data point according to the maximum value and the minimum value of the slope value, and denote them as i min , respectively. i max ; in conjunction with Fig. 3, i min =501, i max =499;
c.计算i
min、i
max之间的数据个数N=|i
max-i
min|+1=|499-501|+1=3；
c. Calculate the number of data between i min and i max N=|i max -i min |+1=|499-501|+1=3;
d.判断所述数据个数N是否大于M：从步骤a、c的计算可知，N=3、M=9，即N＜M，则在最接近的一个步骤即步骤a所述斜率数组中将i
min、i
max之间的数据点及端点排除，剩余的斜率值组成新的斜率数组，在新的斜率数组中重复进行上述步骤b-d直至满足N＞M，具体循环计算过程如下：
d. Judging whether the number of data N is greater than M: From the calculation of steps a and c, it can be seen that N=3, M=9, that is, N<M, then in the closest step, that is, the slope array in step a The data points and endpoints between i min and i max are excluded, the remaining slope values form a new slope array, and the above steps bd are repeated in the new slope array until N>M is satisfied. The specific cyclic calculation process is as follows:
将步骤a所述斜率数组中输出数值序号为499-501对应的3个斜率值删除，剩余的995个斜率值组成新的斜率数组；遍历含有剩余的995个斜率值的斜率数组找到斜率值的最大值和最小值，根据斜率值的最大值、最小值找到其对应的荧光强度数据点对应的输出数值序号，分别记为i
min、i
max；结合图3可见，i
min=592,i
max=486；计算i
min、i
max之间的数据个数N=|i
max-i
min|+1=|486-592|+1=104；判断所述数据个数N是否大于M：由于M=9，而上述该循环计算中的N=104，N＞M，那么含有剩余的995个斜率值的斜率数组中其对应的输出数值序号在486至592之间的斜率值中斜率值绝对值最小的值，确定该斜率值绝对值最小的值对应的输出数值序号i=536，该输出数值序号i=536对应的荧光强度值则为波峰值，其位置则为波峰位置。
Delete the 3 slope values corresponding to the output numerical serial numbers 499-501 in the slope array described in step a, and the remaining 995 slope values form a new slope array; traverse the slope array containing the remaining 995 slope values to find the slope value. The maximum value and the minimum value, according to the maximum value and the minimum value of the slope value, find the corresponding output numerical serial number of the corresponding fluorescence intensity data point, which are respectively recorded as i min and i max ; it can be seen from Figure 3 that i min =592, i max =486; calculate the number of data between i min and i max N=|i max -i min |+1=|486-592|+1=104; determine whether the number of data N is greater than M: because M =9, and N=104 in the above loop calculation, N>M, then in the slope array containing the remaining 995 slope values, its corresponding output numerical serial number is between 486 and 592. The absolute value of the slope value The smallest value, determine the output value number i=536 corresponding to the value with the smallest absolute value of the slope value, the fluorescence intensity value corresponding to the output value number i=536 is the peak value, and its position is the peak position.
本实施例详细说明了测试所得荧光强度散点图仅存在1个突变点时，采用本申请所述用于荧光强度峰检测的方法如何具体确定正确的波峰值、波峰位置；按照本申请所述用于荧光强度峰检测的方法处理，如果有2个及以上的突变点，则可以通过多次循环步骤b-d从而排除掉多个突变点，在此不再详述。This example describes in detail how to specifically determine the correct wave peak and wave peak position using the method for fluorescence intensity peak detection described in this application when there is only one mutation point in the fluorescence intensity scattergram obtained from the test; In the method for detecting fluorescence intensity peaks, if there are 2 or more mutation points, multiple mutation points can be eliminated through multiple cycles of steps b-d, which will not be described in detail here.
现有波峰值的确定方法是通过直接选取荧光强度值最大值作为波峰值，或选择若干个数值相对较大的荧光强度值，取它们的平均值作为波峰值。由此可见，若采用现有波峰值的确定方法在图2所述荧光强度散点图中确定波峰值，由于出现的突变点对应的荧光强度值较大甚至大于波峰值，因此该突变点将会被纳入计算导致波峰值偏大，从而导致不能准确地确定波峰值及其对应的位置，影响检测准确性。The existing wave peak value is determined by directly selecting the maximum value of the fluorescence intensity value as the wave peak value, or selecting several fluorescence intensity values with relatively large values, and taking their average value as the wave peak value. It can be seen that, if the existing wave peak determination method is used to determine the wave peak value in the fluorescence intensity scattergram shown in Figure 2, since the fluorescence intensity value corresponding to the abrupt change point is larger or even greater than the wave peak value, the sudden change point will be It will be included in the calculation, resulting in a large wave peak value, which leads to the inability to accurately determine the wave peak value and its corresponding position, which affects the detection accuracy.
而采用本申请所述用于荧光强度峰检测的方法，由于突变点的出现，导致该突变点的前后输出数值序号对应的斜率值出现较大的变化（变化趋势过大），因此在波峰值确定的过程中会被排除，从而有效解决由于突变点的出现影响波峰值及其对应的位置确定的准确性，尤其在突变点偏大时其效果更为显著。However, using the method for detecting the fluorescence intensity peak described in this application, due to the appearance of the mutation point, the slope value corresponding to the output numerical sequence number before and after the mutation point changes greatly (the change trend is too large). In the process of determination, it will be excluded, so as to effectively solve the problem that the appearance of the mutation point affects the accuracy of the determination of the wave peak and its corresponding position, especially when the mutation point is too large, the effect is more significant.
以上所述仅是本发明的较佳实施方式，故凡依本发明专利申请范围所述的构造、特征及原理所做的等效变化或修饰，均包括于本发明专利申请范围内。The above descriptions are only the preferred embodiments of the present invention, so all equivalent changes or modifications made according to the structures, features and principles described in the scope of the patent application of the present invention are included in the scope of the patent application of the present invention.
Claims (7)
- 一种用于荧光强度峰检测的方法，其特征在于，包括以下步骤：A method for detecting a fluorescence intensity peak, comprising the following steps:（1）采用荧光免疫分析仪以设定的测试频率f、测试时间t对试剂卡进行测试，按序输出荧光强度值，从而在以输出数值序号i为横坐标、荧光强度值F为纵坐标的坐标系上得到（f×t）个荧光强度点，得到对应的荧光强度散点图；(1) Use a fluorescence immunoassay analyzer to test the reagent card with the set test frequency f and test time t, and output the fluorescence intensity values in sequence, so that the output numerical serial number i is the abscissa and the fluorescence intensity value F is the ordinate. (f×t) fluorescence intensity points are obtained on the coordinate system of , and the corresponding fluorescence intensity scattergram is obtained;（2）峰检测：从荧光强度散点图中找出波峰位置，包括以下步骤：(2) Peak detection: find out the peak position from the fluorescence intensity scattergram, including the following steps:a.计算荧光强度点的斜率值K i，将斜率值以其对应的输出数值序号i从小到大作为排列顺序进行排序得到对应的斜率数组；设定斜率数组中斜率值最大值与最小值间数据个数的最小阈值M ； a. Calculate the slope value K i of the fluorescence intensity point, sort the slope value with its corresponding output value number i from small to large as the arrangement order to obtain the corresponding slope array; set the slope value between the maximum value and the minimum value in the slope array The minimum threshold M of the number of data;b.遍历上一步骤所得的斜率数组找出斜率值的最大值和最小值，根据斜率值的最大值、最小值找到其对应的荧光强度数据点对应的输出数值序号，分别记为i min、i max； b. Traverse the slope array obtained in the previous step to find the maximum value and the minimum value of the slope value, and find the output numerical sequence number corresponding to the corresponding fluorescence intensity data point according to the maximum value and the minimum value of the slope value, and denote them as i min , respectively. i max ;c.计算i min、i max之间的数据个数N=|i max-i min|+1； c. Calculate the number of data between i min and i max N=|i max -i min |+1;d.判断所述数据个数N是否大于M：若N≤M，则在最接近的一个步骤所得斜率数组中将i min、i max之间的数据点及端点排除，剩余的斜率值组成新的斜率数组，在新的斜率数组中重复进行上述步骤b-d直至满足N＞M；若N＞M，则在最接近的一个步骤所得斜率数组中其对应的输出数值序号在i min至i max之间的斜率值中找到斜率值绝对值最小的值，确定该斜率值绝对值最小的值对应的输出数值序号i，该输出数值序号i对应的荧光强度值则为波峰值，其位置则为波峰位置。 d. Determine whether the number of data N is greater than M: if N≤M, exclude the data points and endpoints between i min and i max in the slope array obtained by the closest step, and the remaining slope values form a new , repeat the above steps bd in the new slope array until N>M is satisfied; if N>M, the corresponding output numerical sequence number in the slope array obtained in the closest step is between i min and i max Find the value with the smallest absolute value of the slope value among the slope values between the two, and determine the output value number i corresponding to the value with the smallest absolute value of the slope value. The fluorescence intensity value corresponding to the output value number i is the peak value, and its position is the peak value. Location.
- 根据权利要求要求1所述用于荧光强度峰检测的方法，其特征在于，步骤a中所述斜率值K i按以下公式计算： The method for fluorescence intensity peak detection according to claim 1, wherein the slope value K i in step a is calculated according to the following formula:；其中，x为输出数值序号；y为输出数值序号对应的荧光强度值，i=m+1，m+2…, f×t-m；m为预设的已知整数，m≥1。
; where, x is the serial number of the output value; y is the fluorescence intensity value corresponding to the serial number of the output value, i=m+1, m+2…, f×tm; m is a preset known integer, m≥1. - 根据权利要求要求2所述用于荧光强度峰检测的方法，其特征在于，所述m为1或2。The method for fluorescence intensity peak detection according to claim 2, wherein the m is 1 or 2.
- 根据权利要求要求2所述用于荧光强度峰检测的方法，其特征在于，所述最小阈值M =A×（2m+1），A=希望排除突变点的个数。The method for detecting fluorescence intensity peaks according to claim 2, characterized in that, the minimum threshold M=A×(2m+1), A=the number of mutation points desired to be excluded.
- 根据权利要求要求4所述用于荧光强度峰检测的方法，其特征在于，所述A为3或4。The method for fluorescence intensity peak detection according to claim 4, wherein the A is 3 or 4.
- 根据权利要求要求1所述用于荧光强度峰检测的方法，其特征在于，步骤（1）中所述荧光免疫分析仪采用ADUCM360作为数据采集主芯片，所述设定的测试频率f为500Hz，所述测试时间t为2s。The method for detecting fluorescence intensity peaks according to claim 1, wherein in step (1), the fluorescence immunoassay analyzer adopts ADUCM360 as the main chip for data acquisition, and the set test frequency f is 500 Hz, The test time t is 2s.
- 根据权利要求要求1所述用于荧光强度峰检测的方法，其特征在于，所述步骤（2）前还包括将步骤（1）所述荧光强度点进行平滑滤波处理，从而得到处理后的荧光强度散点图。The method for detecting fluorescence intensity peaks according to claim 1, characterized in that before the step (2), the method further comprises smoothing and filtering the fluorescence intensity points in the step (1), so as to obtain the processed fluorescence Intensity scatter plot.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011276183.1A CN112345759B (en) | 2020-11-16 | 2020-11-16 | Method for detecting fluorescence intensity peak |
| CN202011276183.1 | 2020-11-16 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2022100236A1 true WO2022100236A1 (en) | 2022-05-19 |
Family
ID=74363924
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/CN2021/116363 WO2022100236A1 (en) | 2020-11-16 | 2021-09-03 | Method for fluorescence intensity peak detection |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112345759B (en) |
| WO (1) | WO2022100236A1 (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117786281A (en) * | 2024-02-23 | 2024-03-29 | 中国海洋大学 | Optimization calculation method for deposition rate and error of deposit columnar sample |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112345759B (en) * | 2020-11-16 | 2022-07-12 | 三诺生物传感股份有限公司 | Method for detecting fluorescence intensity peak |
| CN113376384A (en) * | 2021-05-19 | 2021-09-10 | 棒糖科技(杭州)股份有限公司 | Method for calculating hormone concentration of fluorescent test paper |
| CN113761456A (en) * | 2021-09-07 | 2021-12-07 | 杭州凯曼健康科技有限公司 | Immunofluorescence chromatography curve analysis method and device and electronic equipment |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103969379A (en) * | 2013-01-29 | 2014-08-06 | 深圳普门科技有限公司 | An acquisition method for obtaining the measured value of a liquid chromatography test |
| CN104458981A (en) * | 2013-09-18 | 2015-03-25 | 株式会社岛津制作所 | Waveform processing assistance method and system |
| CN104730181A (en) * | 2013-12-18 | 2015-06-24 | 北京普源精电科技有限公司 | Chromatographic peak end point adjusting method and chromatographic work station having chromatographic peak end point adjusting function |
| TW201811261A (en) * | 2016-09-20 | 2018-04-01 | 捷騰光電股份有限公司 | Signal detection method for accurately detecting the peak-to-peak interval of the photoplethysmography signal |
| CN111248895A (en) * | 2018-11-30 | 2020-06-09 | 无锡祥生医疗科技股份有限公司 | Electrocardiosignal characteristic detection method and system |
| CN112345759A (en) * | 2020-11-16 | 2021-02-09 | 三诺生物传感股份有限公司 | Method for detecting fluorescence intensity peak |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6070097A (en) * | 1998-12-30 | 2000-05-30 | General Electric Company | Method for generating a gating signal for cardiac MRI |
| CN105030233B (en) * | 2015-07-08 | 2017-10-13 | 上海师范大学 | A kind of electrocardiosignal ST sections of recognition methods |
| CN106974617A (en) * | 2016-01-19 | 2017-07-25 | 深圳市卡迪赛克科技有限公司 | The Signal Pre-Processing Method and signal wave crest detection method of a kind of efficiently and accurately |
| CA3065208A1 (en) * | 2018-12-20 | 2020-06-20 | Queen's University At Kingston | Long qt syndrome diagnosis and classification |
| CN109859188B (en) * | 2019-01-31 | 2021-04-06 | 领航基因科技（杭州）有限公司 | Fluorescence crosstalk correction method based on mean shift algorithm and application thereof |
| CN111259311B (en) * | 2020-01-14 | 2023-03-24 | 西安应用光学研究所 | Peak noise processing method |
| CN111879744A (en) * | 2020-08-06 | 2020-11-03 | 深圳市锦瑞生物科技有限公司 | Method for detecting concentration of substance to be detected, fluorescence immunoassay analyzer and storage medium |
-
2020
- 2020-11-16 CN CN202011276183.1A patent/CN112345759B/en active Active
-
2021
- 2021-09-03 WO PCT/CN2021/116363 patent/WO2022100236A1/en active Application Filing
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103969379A (en) * | 2013-01-29 | 2014-08-06 | 深圳普门科技有限公司 | An acquisition method for obtaining the measured value of a liquid chromatography test |
| CN104458981A (en) * | 2013-09-18 | 2015-03-25 | 株式会社岛津制作所 | Waveform processing assistance method and system |
| CN104730181A (en) * | 2013-12-18 | 2015-06-24 | 北京普源精电科技有限公司 | Chromatographic peak end point adjusting method and chromatographic work station having chromatographic peak end point adjusting function |
| TW201811261A (en) * | 2016-09-20 | 2018-04-01 | 捷騰光電股份有限公司 | Signal detection method for accurately detecting the peak-to-peak interval of the photoplethysmography signal |
| CN111248895A (en) * | 2018-11-30 | 2020-06-09 | 无锡祥生医疗科技股份有限公司 | Electrocardiosignal characteristic detection method and system |
| CN112345759A (en) * | 2020-11-16 | 2021-02-09 | 三诺生物传感股份有限公司 | Method for detecting fluorescence intensity peak |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117786281A (en) * | 2024-02-23 | 2024-03-29 | 中国海洋大学 | Optimization calculation method for deposition rate and error of deposit columnar sample |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112345759A (en) | 2021-02-09 |
| CN112345759B (en) | 2022-07-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2022100236A1 (en) | Method for fluorescence intensity peak detection | |
| WO2022100234A1 (en) | Method for fluorescence intensity basis calculation | |
| CN104515725B (en) | A kind of method and system and its cytoanalyze for recognizing abnormal particle | |
| CN106645708A (en) | Quantitative detection calculation method based on fluorescent immuno-chromatographic technology | |
| CN105334243A (en) | Blood analyzer | |
| CN105044332B (en) | A kind of method for monitoring HOOK effects in immune colloid gold experiment | |
| CN107312850A (en) | A kind of detection method of the invalid amplifications of PCR | |
| CN110231328B (en) | Raman spectrum quantitative analysis method based on half-peak height distance method | |
| CN105403712A (en) | High performance detection kit for human urine alpha 1 acidoglycoprotein | |
| CN105334317A (en) | Anti-LKM (liver-kidney microsomal) 1 antibody detection kit and detection method | |
| Li et al. | Moving rate of positive patient results as a quality control tool for high-sensitivity cardiac troponin T assays | |
| CN108614017A (en) | A kind of prothrombin time detection method and device | |
| CN111366692A (en) | Gas environment parameter monitoring system and method | |
| DE60102797D1 (en) | FLUORESCENCE INTENSITY ANALYSIS USING A VARIETY OF DISTRIBUTIONS: COMPETITIVE DETERMINATION OF DIFFUSION TIME AND MOLECULAR BRIGHTNESS | |
| CN103616522A (en) | Immune chromatography result identification method based on envelope area and twice compensation | |
| CN1220043C (en) | Ozone concentration detection method and device used for ozone box | |
| CN108931516B (en) | System parameter optimization method capable of saving sample introduction amount and serum element quantitative analysis method | |
| CN211978750U (en) | NO in air2Monitoring system | |
| CN108226238A (en) | A kind of method that nitrite concentration is measured using interdigital electrode | |
| CN113640253A (en) | Turbidity detection method | |
| CN113219015A (en) | Method and device for detecting HIV-1P24 antigen based on silicon nitride (SiNx) solid nano-pores | |
| CN206609868U (en) | A kind of ammonia nitrogen water quality analyzer | |
| Tian et al. | Construction of high sensitivity electrochemiluminescence sensor and its application in nt-probnp detection | |
| CN111339498B (en) | Rapid correction method and system for fluorescence dissolved organic matter data | |
| CN116539831B (en) | Water environment data monitoring processing method based on big data analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 21890764 Country of ref document: EP Kind code of ref document: A1 |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 21890764 Country of ref document: EP Kind code of ref document: A1 |