CN107101980A - The fluorescence detector virtual value compensation method of microscratch amount and system - Google Patents
The fluorescence detector virtual value compensation method of microscratch amount and system Download PDFInfo
- Publication number
- CN107101980A CN107101980A CN201710109577.XA CN201710109577A CN107101980A CN 107101980 A CN107101980 A CN 107101980A CN 201710109577 A CN201710109577 A CN 201710109577A CN 107101980 A CN107101980 A CN 107101980A
- Authority
- CN
- China
- Prior art keywords
- value
- data
- peak value
- diaphragm
- mentioned
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/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
- 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"
- G01N21/643—Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes" non-biological material
-
- 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"
- G01N2021/6432—Quenching
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2201/00—Features of devices classified in G01N21/00
- G01N2201/12—Circuits of general importance; Signal processing
- G01N2201/127—Calibration; base line adjustment; drift compensation
Landscapes
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Optics & Photonics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Investigating Or Analysing Materials By The Use Of Chemical Reactions (AREA)
Abstract
The present invention relates to a kind of microscratch amount fluorescence detector virtual value compensation method, including：A. collection obtains the initial data of diaphragm fluorescent value；B. smothing filtering is carried out to above-mentioned initial data；C. according to the data after above-mentioned smothing filtering, calculating obtains detecting peak value；D. above-mentioned detection peak value is modified, to be compensated to fluorescent material detection virtual value.The invention further relates to a kind of microscratch amount fluorescence detector virtual value compensation system.The present invention can carry out dynamic compensation data, the sensitivity of stabilizer instrument and rate of false alarm according to diaphragm use state, increase the service life of diaphragm.
Description
Technical field
The present invention relates to a kind of microscratch amount fluorescence detector virtual value compensation method and system.
Background technology
Microscratch amount fluorescence detector compares traditional explosive detection instrument, the characteristics of having high sensitivity, fast response time,
Have wide practical use in the important safety check place such as subway, airport.The core of this quasi-instrument is to be coated with the glass substrate
Machine thin-film material, this material has the attribute of fluorescent quenching to explosive molecule.When there is explosive molecule contacts thin-film material,
Fluorescence volume can be made to diminish, instrument determines whether there is the presence of explosive molecule by the change of fluorescence intensity.
Inevitably, the same with most of organic fluorescence materials, this kind of organic film material is due to photobleaching, diaphragm
The factors such as pollution, as the fluorescence intensity using diaphragm can slowly die down, diaphragm has certain life-span limitation.Entirely using
Time range in, the variable quantity of fluorescent quenching also can dying down and diminish with diaphragm integral fluorescence intensity, if with fixation
Threshold value come judge whether detection explosive when, such situation just occurs：Some materials that new diaphragm can be detected
And can't detect using old diaphragm, the state of the detection sensitivity heavy dependence diaphragm of instrument.Because new diaphragm is to finite concentration
The change in fluorescence that material is produced exceedes alarm threshold value, and in the later stage used, with the consumption of diaphragm, same amount of material is produced
Change in fluorescence amount again be less than alarm threshold value, it is impossible to alarm.
Can all there are two parameters of rate of false alarm and sensitivity for most explosive detection instrument, and the two are crucial
Parameter mutually restriction again, one of parameter is improved simply can cause the reduction of another parameter, for fluoroscopic examination class instrument
Device is particularly true, because fluorescence diaphragm is running stores, it is necessary to be regularly replaced according to service condition, the diaphragm newly changed is sensitiveer
With actively, data dithering is relatively large, easily wrong report, in order to not reduce rate of false alarm, it is necessary to which alarm threshold value is set relatively
Height, such rate of false alarm and sensitivity obtain an active balance, but with using, more and more blunt, the data dithering that diaphragm becomes
Also it is relatively small, but in the case that alarm threshold value is constant, some sensitivity can be caused to decline rapidly, cause what some should be measured
Thing can't detect.
The content of the invention
In view of this, it is necessary to which a kind of microscratch amount fluorescence detector virtual value compensation method and system are provided.
The present invention provides a kind of microscratch amount fluorescence detector virtual value compensation method, and this method comprises the following steps：A. adopt
Collection obtains the initial data of diaphragm fluorescent value；B. smothing filtering is carried out to above-mentioned initial data；C. according to above-mentioned smothing filtering after
Data, calculating obtain detect peak value；D. above-mentioned detection peak value is modified, to be mended to fluorescent material detection virtual value
Repay.
Wherein, described step a is specifically included：The initial data A for obtaining diaphragm fluorescent value is gathered using analog-digital converter
(A1,A2,A3,,,,,,,An-1,An,An+1)。
Described step b is specifically included：Wn=(An-j,+An-j+1,+An,+An+j-1,+An+j)/(2j-1), obtain data
(W1,W2,W3,,,,,,,Wn-1,Wn,Wn+1), wherein j is constant.
Described step c is specifically included：Detection peak value H is asked for using differential methodi=Wi-n-Wi, obtain data (H1,H2,
H3,,,,,,,Hn-1,Hn,Hn+1)。
Described step d is specifically included：Using formula：Jn=(Q/Wn)b* Hn is modified to above-mentioned detection peak value, its
In, Q and b are constant.
The present invention also provides a kind of microscratch amount fluorescence detector virtual value compensation system, and the system includes acquisition module, filter
Ripple module, computing module and compensating module, wherein：The acquisition module is used to gather the initial data for obtaining diaphragm fluorescent value；
The filtration module is used to carry out smothing filtering to above-mentioned initial data；After the computing module is used for according to above-mentioned smothing filtering
Data, calculating obtain detect peak value；The compensating module is used to be modified above-mentioned detection peak value, to be examined to fluorescent material
Virtual value is surveyed to compensate.
Wherein, described acquisition module specifically for：The original number for obtaining diaphragm fluorescent value is gathered using analog-digital converter
According to A (A1,A2,A3,,,,,,,An-1,An,An+1)。。
Described filtration module specifically for：Wn=(An-j,+An-j+1,+An,+An+j-1,+An+j)/(2j-1), obtain data
(W1,W2,W3,,,,,,,Wn-1,Wn,Wn+1), wherein j is constant.
Described computing module specifically for：Detection peak value H is asked for using differential methodi=Wi-n-Wi, obtain data (H1,
H2,H3,,,,,,,Hn-1,Hn,Hn+1)。
Described compensating module specifically for：Using formula：Jn=(Q/Wn)b* Hn is modified to above-mentioned detection peak value,
Wherein, Q and b is constant.
Microscratch amount fluorescence detector virtual value compensation method of the present invention and system, can enter action according to diaphragm use state
The compensation data of state, the sensitivity of stabilizer instrument and rate of false alarm, increase the service life of diaphragm.
Brief description of the drawings
Fig. 1 is a kind of flow chart of microscratch amount fluorescence detector virtual value compensation method of the invention；
Fig. 2 is a kind of hardware architecture diagram of microscratch amount fluorescence detector virtual value compensation system of the invention.
Embodiment
Below in conjunction with the accompanying drawings and specific embodiment the present invention is further detailed explanation.
As shown in fig.1, being a kind of work of microscratch amount fluorescence detector virtual value compensation method preferred embodiment of the invention
Industry flow chart.
Step S1, collection obtains the initial data of diaphragm fluorescent value.Specifically：
The original for obtaining diaphragm fluorescent value is gathered using ADC (Analog-to-Digital Converter, analog-digital converter)
Beginning data A (A1,A2,A3,,,,,,,An-1,An,An+1)。
Step S2, smothing filtering is carried out to above-mentioned initial data.It is specific as follows：
Initial data A is subjected to smothing filtering, to eliminate big burr, specific method is Wn=(An-j,+An-j+1,+
An,+An+j-1,+An+j)/(2j-1), obtain data (W1,W2,W3,,,,,,,Wn-1,Wn,Wn+1), wherein j is constant.It is appreciated that
, according to the actual conditions of filtering, if noise is very serious, j can be larger, if noise is not serious, and j can take small by one
Point.J=3 is taken in the present embodiment.
Step S3, according to the data after above-mentioned smothing filtering, calculating obtains detecting peak value.Specifically：
Detection peak value H is asked for using differential methodi=Wi-n-Wi, obtain data (H1,H2,H3,,,,,,,Hn-1,Hn,Hn+1)。
Wherein, time span n is set, such as the present embodiment sampling 10 per second according to sampling rate and the difference of substance detection time
Data, and single sample detection maximum time is 10 seconds, then the present embodiment just takes N=10*10.
Step S4, is modified to above-mentioned detection peak value, to be compensated to fluorescent material detection virtual value.It is specific and
Speech：
Correction formula is：Jn=(Q/Wn)b*Hn。
Wherein, Q and b is constant, and Q is relevant with the maximum of sampled data, and b is relevant with diaphragm sensitivity and fluorescence decay.
In the present embodiment, the digital-to-analogue converter of 16 is used, therefore changes into maximum after the decimal system and is
65535, therefore Q is set to Q=65535.B is relevant with diaphragm sensitivity and fluorescence decay, can be drawn from many experiments sensitive
Degree and fluorescence decay by b into being set to b=2 in quadratic relationship, therefore the present embodiment.
As shown in fig.2, being a kind of hardware architecture diagram of microscratch amount fluorescence detector virtual value compensation system of the invention.Should
System includes acquisition module, filtration module, computing module and compensating module.
The acquisition module is used to gather the initial data for obtaining diaphragm fluorescent value.Specifically：
The acquisition module obtains film using ADC (Analog-to-Digital Converter, analog-digital converter) collections
Initial data A (the A of piece fluorescent value1,A2,A3,,,,,,,An-1,An,An+1)。
The filtration module is used to carry out smothing filtering to above-mentioned initial data.It is specific as follows：
Initial data A is carried out smothing filtering by the filtration module, and to eliminate big burr, specific method is Wn=
(An-j,+An-j+1,+An,+An+j-1,+An+j)/(2j-1), obtain data (W1,W2,W3,,,,,,,Wn-1,Wn,Wn+1), wherein j is
Constant.It is understood that according to the actual conditions of filtering, if noise is very serious, j can be larger, if noise is not tight
Weight, j can take a little bit smaller.J=3 is taken in the present embodiment.
The computing module is used for according to the data after above-mentioned smothing filtering, and calculating obtains detecting peak value.Specifically：
The computing module asks for detection peak value H using differential methodi=Wi-n-Wi, obtain data (H1,H2,H3,,,,,,,
Hn-1,Hn,Hn+1).Wherein, time span n is set, such as the present embodiment according to sampling rate and the difference of substance detection time
10 data of sampling per second, and single sample detection maximum time is 10 seconds, then the present embodiment just takes N=10*10.
The compensating module is used to be modified above-mentioned detection peak value, to be mended to fluorescent material detection virtual value
Repay.Specifically：
Correction formula is：Jn=(Q/Wn)b*Hn。
Wherein, Q and b is constant, and Q is relevant with the maximum of sampled data, and b is relevant with diaphragm sensitivity and fluorescence decay.
In the present embodiment, the digital-to-analogue converter of 16 is used, therefore changes into maximum after the decimal system and is
65535, therefore Q is set to Q=65535.B is relevant with diaphragm sensitivity and fluorescence decay, can be drawn from many experiments sensitive
Degree and fluorescence decay by b into being set to b=2 in quadratic relationship, therefore the present embodiment.
The present invention can be derived that effective penalty coefficient based on a kind of normalization backoff algorithm, the algorithm.Make diaphragm whole
Similar responsiveness is shown in individual life span, so that the sensitivity of instrument will not drastically decline with the use of diaphragm,
Meanwhile, the present invention also increases the Acceptable life of diaphragm.
Although the present invention is described with reference to current better embodiment, those skilled in the art should be able to manage
Solution, above-mentioned better embodiment is only used for illustrating the present invention, any in the present invention not for limiting protection scope of the present invention
Spirit and spirit within, any modification, equivalence replacement, improvements for being done etc., should be included in the present invention right guarantor
Within the scope of shield.
Claims (10)
1. a kind of microscratch amount fluorescence detector virtual value compensation method, it is characterised in that this method comprises the following steps：
A. collection obtains the initial data of diaphragm fluorescent value；
B. smothing filtering is carried out to above-mentioned initial data；
C. according to the data after above-mentioned smothing filtering, calculating obtains detecting peak value；
D. above-mentioned detection peak value is modified, to be compensated to fluorescent material detection virtual value.
2. the method as described in claim 1, it is characterised in that described step a is specifically included：
Initial data A (the A for obtaining diaphragm fluorescent value are gathered using analog-digital converter1,A2,A3,,,,,,,An-1,An,An+1)。
3. method as claimed in claim 2, it is characterised in that described step b is specifically included：
Wn=(An-j,+An-j+1,+An,+An+j-1,+An+j)/(2j-1), obtain data (W1,W2,W3,,,,,,,Wn-1,Wn,Wn+1),
Wherein j is constant.
4. method as claimed in claim 3, it is characterised in that described step c is specifically included：
Detection peak value H is asked for using differential methodi=Wi-n-Wi, obtain data (H1,H2,H3,,,,,,,Hn-1,Hn,Hn+1)。
5. method as claimed in claim 4, it is characterised in that described step d is specifically included：
Using formula：Jn=(Q/Wn)b* Hn is modified to above-mentioned detection peak value, wherein, Q and b are constant.
6. a kind of microscratch amount fluorescence detector virtual value compensation system, it is characterised in that the system includes acquisition module, filtering mould
Block, computing module and compensating module, wherein：
The acquisition module is used to gather the initial data for obtaining diaphragm fluorescent value；
The filtration module is used to carry out smothing filtering to above-mentioned initial data；
The computing module is used for according to the data after above-mentioned smothing filtering, and calculating obtains detecting peak value；
The compensating module is used to be modified above-mentioned detection peak value, to be compensated to fluorescent material detection virtual value.
7. system as claimed in claim 6, it is characterised in that described acquisition module specifically for：
Initial data A (the A for obtaining diaphragm fluorescent value are gathered using analog-digital converter1,A2,A3,,,,,,,An-1,An,An+1)。
8. system as claimed in claim 7, it is characterised in that described filtration module specifically for：
Wn=(An-j,+An-j+1,+An,+An+j-1,+An+j)/(2j-1), obtain data (W1,W2,W3,,,,,,,Wn-1,Wn,Wn+1),
Wherein j is constant.
9. system as claimed in claim 8, it is characterised in that described computing module specifically for：
Detection peak value H is asked for using differential methodi=Wi-n-Wi, obtain data (H1,H2,H3,,,,,,,Hn-1,Hn,Hn+1)。
10. system as claimed in claim 9, it is characterised in that described compensating module specifically for：
Using formula：Jn=(Q/Wn)b* Hn is modified to above-mentioned detection peak value, wherein, Q and b are constant.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710109577.XA CN107101980B (en) | 2017-02-27 | 2017-02-27 | Effective value compensation method and system for micro trace fluorescence detector |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710109577.XA CN107101980B (en) | 2017-02-27 | 2017-02-27 | Effective value compensation method and system for micro trace fluorescence detector |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107101980A true CN107101980A (en) | 2017-08-29 |
| CN107101980B CN107101980B (en) | 2020-09-25 |
Family
ID=59675543
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710109577.XA Active CN107101980B (en) | 2017-02-27 | 2017-02-27 | Effective value compensation method and system for micro trace fluorescence detector |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107101980B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100032583A1 (en) * | 2005-04-25 | 2010-02-11 | Polestar Technologies, Inc. | Optical Chemical Sensor Feedback Control System |
| CN202057569U (en) * | 2011-05-18 | 2011-11-30 | 北京汇丰隆经济技术开发有限公司 | Automatic correction system for aerosol particle analyzer |
| CN103140753A (en) * | 2010-10-08 | 2013-06-05 | 梅特勒-托利多公开股份有限公司 | Method for operating an optochemical sensor and measuring instrument |
| CN203405428U (en) * | 2013-04-26 | 2014-01-22 | 上海精谱科技有限公司 | Energy dispersion x fluorescence analyzer |
| CN104641219A (en) * | 2012-08-21 | 2015-05-20 | 奥里巴Abx股份有限公司 | Method for compensating for the aging of a reagent during fluorescence measurements carried out on particles, and biological analysis device implementing the method |
| JP2015222179A (en) * | 2014-05-22 | 2015-12-10 | 凸版印刷株式会社 | Fluorescent chip for correction |
-
2017
- 2017-02-27 CN CN201710109577.XA patent/CN107101980B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100032583A1 (en) * | 2005-04-25 | 2010-02-11 | Polestar Technologies, Inc. | Optical Chemical Sensor Feedback Control System |
| CN103140753A (en) * | 2010-10-08 | 2013-06-05 | 梅特勒-托利多公开股份有限公司 | Method for operating an optochemical sensor and measuring instrument |
| CN202057569U (en) * | 2011-05-18 | 2011-11-30 | 北京汇丰隆经济技术开发有限公司 | Automatic correction system for aerosol particle analyzer |
| CN104641219A (en) * | 2012-08-21 | 2015-05-20 | 奥里巴Abx股份有限公司 | Method for compensating for the aging of a reagent during fluorescence measurements carried out on particles, and biological analysis device implementing the method |
| CN203405428U (en) * | 2013-04-26 | 2014-01-22 | 上海精谱科技有限公司 | Energy dispersion x fluorescence analyzer |
| JP2015222179A (en) * | 2014-05-22 | 2015-12-10 | 凸版印刷株式会社 | Fluorescent chip for correction |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107101980B (en) | 2020-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hansen et al. | tweezercalib 2.0: Faster version of MatLab package for precise calibration of optical tweezers | |
| JP4397692B2 (en) | System and method for multi-photon counting with time correlation | |
| CN109073536B (en) | Information processing device, information processing method, program, and information processing system | |
| US11781974B2 (en) | Method for detecting gas concentration in glass bottle with dynamical threshold adjustment | |
| CN101069093A (en) | Self-adjusting electrochemical sensor | |
| CN110441374B (en) | Dissolved oxygen sensor compensation method based on SA-BP neural network | |
| CN112345759B (en) | Method for detecting fluorescence intensity peak | |
| CN108630330A (en) | Pressurized-water reactor nuclear power plant instrument system detector test process method, apparatus and system | |
| CN104749156A (en) | Raman spectrum detection method | |
| Nakhostin | Recursive algorithms for digital implementation of neutron/gamma discrimination in liquid scintillation detectors | |
| Halter et al. | An automated protocol for performance benchmarking a widefield fluorescence microscope | |
| TW202033939A (en) | Cavity ring-down spectroscopy having interleaved data acquisition for interference mitigation | |
| CN111521883A (en) | Method and system for obtaining electric field measurement value of high-voltage direct-current transmission line | |
| CN108268025B (en) | Elasticity evaluation method for networked control system under random disturbance | |
| CN107101980A (en) | The fluorescence detector virtual value compensation method of microscratch amount and system | |
| Casciati et al. | Structural health monitoring by Lyapunov exponents of non‐linear time series | |
| CN117556327A (en) | Method for detecting diseases and insect pests of red-tassel sorghum | |
| CN108169215A (en) | A kind of computational methods of emission spectrometer time of integration upper limit setting | |
| US11525781B2 (en) | Method for measuring oxygen and apparatus for measuring oxygen | |
| CN106248653A (en) | A kind of method improving LIBS quantitative analysis long-time stability | |
| CN109933883A (en) | A kind of current in resistance property Growth Rate Calculation method based on multiple linear regression | |
| CN115980230A (en) | Liquid chromatogram-mass spectrum combined instrument calibration device and calibration method thereof | |
| Ahmadzadeh et al. | Automated analysis of individual particles using a commercial capillary electrophoresis system | |
| CN117132504B (en) | Method for removing interference based on single photon spectrum counting technology imaging | |
| CN110333148B (en) | Soil dynamic shear modulus testing method based on fine analysis of vibration attenuation curve |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 518000 six layer C unit of A building, Shenfu Bao science and Technology Industrial Park, Futian District Free Trade Zone, Shenzhen, Guangdong Applicant after: Shenzhen sword Defense Technology Co. Ltd. Address before: 518000 six layer C unit of A building, Shenfu Bao science and Technology Industrial Park, Futian District Free Trade Zone, Shenzhen, Guangdong Applicant before: Shaanxi Normal University |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |