The present disclosure relates to detection, and more particularly, to a method for reducing noise in AES detection.
ICP-AES has become a main member among laboratory analysis instruments due to advantages of a low detection limit, a wide dynamic range, a wide detection waveband and the like, and is widely used in the fields of environments, semiconductors, medical treatment, food, metallurgy, nuclear industry and the like. In recent years, with the innovation of technologies and the increase in low content demands of detection substances, the improvement of ICP detection performance has become one of hotspots of manufacturers. The detection performance is directly derived from two factors, signal response and noise, and how to improve the signal-to-noise ratio is a main issue of the industry.
Some manufacturers start from detectors and have customized and developed low noise detectors. However, as restricted by the current domestic process level, this solution has a large cost for domestic instrument manufacturers and cannot ensure stable supply.
Some manufacturers focus on improving the response and improving the detection limit by increasing a light-transmitting aperture and the like. However, this changes the original optical system at the cost of degrading other parameters.
Meanwhile, a scientific research institution proposed and implemented to obtain better photo response by increasing the excitation efficiency by means of solution pretreatment, for example, special treatment near an atomizer to reduce surface adhesion of aerosols. Although this manner provides a certain improvement, reliability, stability, and applicability for all elements still need to be verified, and this manner cannot be industrialized at present.
To solve the defects in the solutions of the related art, the present disclosure provides a method for reducing noise in AES detection to improve the detection limit.
An object of the present disclosure is implemented by the following technical solutions.
A method for reducing noise in AES detection is provided. The method for reducing noise in AES detection includes steps of:
A1, obtaining G based on a sub-array {tilde over (Z)} of detection data:
G=(Sx{tilde over (Z)})2+(Sy{tilde over (Z)})2+(Sr{tilde over (Z)})2+(Sl{tilde over (Z)})2,
where Sx, Sy, Sr, and Sl represents a sobel horizontal gradient operator, a sobel vertical gradient operator, a sobel right diagonal gradient operator, and a sobel left diagonal gradient operator, respectively, and SxZ˜, SyZ˜, SrZ˜, and SrZ˜ each represent convolution processing of a sobel operator on Z˜ with a step of 1, a convolution edge being processed by filling 0 or by symmetrization, an influence of the edge on the result being negligible, and
where
and a size of {tilde over (Z)} is (m, n);
A2, adding one column of zero values to both a left end and a right end of G for only a calculation purpose in this step without changing an original size of G, and recording the added columns as a 0-th column and a (n+1)-th column; and for each element in G, forming a set of data using three elements including the element and a previous element and a next element of the element in a column direction, and sorting the set of data in a descending order: Gi,d≥Gi,e≥Gi,f, to obtain an array {tilde over (D)}:
where α=0.618, β=(1−α)*α, and γ=(1−α)2; and d, e, and f are different from each other and are each taken from a set {j, (j−1), (j+1)}, where 1i≤m and 1≤j≤n;
A3, performing normalization processing on the array {tilde over (D)} to obtain an array D;
A4, adding one column of zero values to both a left end and a right end of {tilde over (Z)} for only a calculation purpose in this step without changing an original size of {tilde over (Z)}, and recording the added columns as a 0-th column and a (n+1)-th column; and for each element in {tilde over (Z)}, forming a set of data using three elements including the element in {tilde over (Z)} and a previous element and a next element of the element in {tilde over (Z)} in the column direction, and sorting the set of data in a descending order: {tilde over (Z)}i,g≥{tilde over (Z)}i,h≥{tilde over (Z)}i,k, to obtain an array U of m rows by n columns:
where g, h, and k are different from each other and are each taken from a set {j, (j−1), (j+1)}; and
where 1≤i≤m and 1≤j≤n;
A5, calculating a noise difference value in the column direction, the noise difference value being an array C of m rows by n−1 columns:
where {tilde over (Z)}j˜max is a maximum value of a column in which the element is located, 1≤i≤m, and 1≤j≤n;
A6, formulating a noise array N of m rows by n columns as:
where x is a column coordinate of a maximum value in the array {tilde over (Z)}; and
A7, constructing a new sub-array PN:
P
i,j
N
={tilde over (Z)}
i,j
−N
i,j,
where 1≤i≤m and 1≤j≤n.
Compared with the related art, the present disclosure has the following beneficial effects.
The noise caused by non-uniform response of a detector device, differences in attributes, and light distribution is similar in each frame, and the influence of the noise is reduced by processing each image frame through an algorithm, thereby reducing the fluctuation in a time domain, reducing data variance, and improving the detection limit.
The disclosed content of the present disclosure will become more understandable by referring to the accompanying drawings. It will be readily understood by those skilled in the art that these drawings are merely used for illustrating the technical solutions of the present disclosure, and are not intended to limit the scope of the present disclosure. In the drawings:
At step A1, a multi-gradient quadratic sum G is obtained based on a sub-array {tilde over (Z)} of detection data. A specific manner is as follows: since an image of a real object generated by an optical system is a Gaussian distribution two-dimensional light spot, and data read by a detector is also a two-dimensional array, introducing Sobel operators including a vertical operator, a horizontal operator, and diagonal operators in two directions, which are recorded as Sx, Sy, Sr, and Sl; and
G=(Sx{tilde over (Z)})2+(Sy{tilde over (Z)})2+(Sr{tilde over (Z)})2+(Sl{tilde over (Z)})2,
where
and a size of {tilde over (Z)} is (m, n).
At step A2, one column of zero values is added to both a left end and a right end of G for only a calculation purpose in this step without changing an original size of G, and the added columns are recorded as a 0-th column and a (n+1)-th column; and for each element in G, a set of data is formed by using three elements comprising the element and a previous element and a next element of the element in a column direction, and the set of data is sorted in a descending order: Gi,d≥Gi,e≥Gi,f, to obtain an array {tilde over (D)}:
where α=0.618, β=(1−α)*α, and γ=(1−α)2; and d, e, and f are different from each other and are each taken from a set {j, (j−1), (j+1)}, where 1≤i≤m and 1≤j≤n.
At step A3, normalization processing is performed on the array {tilde over (D)}, i.e., each element in the array is divided by a sum of values in a column where the element is located to obtain an array D.
At step A4, one column of zero values is added to both a left end and a right end of {tilde over (Z)} for only a calculation purpose in this step without changing an original size of {tilde over (Z)}, and the added columns are recorded as a 0-th column and a (n+1)-th column; and for each element in {tilde over (Z)}, a set of data is formed by using three elements including the element in {tilde over (Z)} and a previous element and a next element of the element in {tilde over (Z)} in the column direction, and the set of data is sorted in a descending order: {tilde over (Z)}i,g≥{tilde over (Z)}i,h≥{tilde over (Z)}i,k, to obtain an array U of in rows by n columns:
where g, h, and k are different from each other and are each taken from a set {j, (j−1), (j+1)}; and
where 1≤i≤m and 1≤j≤n.
At step A5, a noise difference value in the column direction is calculated, and the noise difference value is an array C of m rows by n−1 columns:
where {tilde over (Z)}j˜max is a maximum value of a column in which the element is located, 1≤i≤m, and 1≤j≤n.
At step A6, a noise array N of m rows by n columns is formulated as:
where x is a column coordinate of a maximum value in the array {tilde over (Z)}.
At step A7, a new sub-array PN is constructed:
P
i,j
N
={tilde over (Z)}
i,j
−N
i,j,
where 1≤i≤m and 1≤j≤n.
The method for reducing noise in AES detection according to Example 1 of the present disclosure was applied in ICP-AES.
In this application example, Ba233.525, Cu324.754, Cr283.563, and Mn257.610 were used for testing:
1. A standard curve for Ba233.525, a standard curve for Cu324.754, a standard curve for Cr283.563, and a standard curve for Mn257.610 were constructed and are as illustrated in
2. The instrument detection limit, i.e., 3 times a standard deviation, was detected by loading a blank water sample without adding an algorithm, and the data is as follows (two groups in the middle are taken as controls).
Test Group 1
Test Group 2
1. Example 1 was used to process original sub-array data corresponding to the above test groups, and the results are shown below.
Test Group 1 after Algorithm Processing
Test Group 2 after Algorithm Processing
Number | Date | Country | Kind |
---|---|---|---|
202010539508.4 | Jun 2020 | CN | national |
This application is a continuation of International Application NO. PCT/CN2021/099738, filed on Jun. 11, 2021, which claims priority to Chinese Patent Application No. 202010539508.4, entitled “METHOD FOR REDUCING NOISE IN AES DETECTION” and filed by HANGZHOU PUYU TECHNOLOGY DEVELOPMENT CO., LTD. on Jun. 14, 2020, the entire disclosures of which are incorporated here by reference.
Number | Date | Country | |
---|---|---|---|
Parent | PCT/CN2021/099738 | Jun 2021 | US |
Child | 18080773 | US |