CN108663707B - Multi-time bidirectional S-K smoothing processing system and method - Google Patents
Multi-time bidirectional S-K smoothing processing system and method Download PDFInfo
- Publication number
- CN108663707B CN108663707B CN201810285498.9A CN201810285498A CN108663707B CN 108663707 B CN108663707 B CN 108663707B CN 201810285498 A CN201810285498 A CN 201810285498A CN 108663707 B CN108663707 B CN 108663707B
- Authority
- CN
- China
- Prior art keywords
- filter
- buffer
- filtering
- digital
- data stream
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
Abstract
The invention relates to a multi-time bidirectional S-K smoothing processing system and a method, wherein a first buffer, a forward digital SK filter, a second buffer and a reverse digital SK filter are sequentially connected, the reverse digital SK filter is connected with the first buffer, and the first buffer, the forward digital SK filter, the second buffer and the reverse digital SK filter are respectively connected with a control register and a parameter register. The multi-time bidirectional S-K smoothing processing system and the method provided by the invention have the same effect as multi-time Gaussian smoothing processing, compared with Gaussian smoothing, the multi-time bidirectional S-K smoothing processing system and the method provided by the invention have the advantages of higher execution efficiency, simple arithmetic operation process, convenience in use and easiness in software and hardware realization, and can be applied to nuclear instruments and actual measurement nuclear spectrum analysis to solve the problem that qualitative and quantitative analysis of low-content nuclides is influenced by characteristics such as electronic noise, statistical fluctuation and the like in the conventional nuclear spectrum measurement.
Description
Technical Field
The invention relates to the field of nuclear energy spectrum analysis, in particular to a multi-time bidirectional S-K smoothing processing system and a method.
Background
With the progress and development of science, the application fields of nuclear instruments and meters are more and more extensive, and almost cover the fields of agriculture, industry, medical treatment and health, geological mines, environmental protection, aerospace, teaching, scientific research, customs and the like. The nuclear decay and the statistical fluctuation characteristic in the nuclear energy spectrum measurement can enable a real Gaussian distribution expected value to be submerged in noise when the count is small, and even when the count is large, the statistical fluctuation characteristic can enable actually-measured energy spectrum data to generate a burr phenomenon to influence later-period data analysis; in order to reliably identify the existence of characteristic peaks of low-content nuclides and simultaneously carry out qualitative analysis on the peak energy and the position, a new energy spectrum analysis method needs to be researched.
An S-K filter was proposed by r.p.sallen and e.l.key in 1955, and a low-pass S-K filter is widely applied to filtering and shaping of nuclear pulse signals in nuclear electronics, and can obtain output of approximate gaussian waveform in fewer stages, and when nuclear energy spectrum data is reversely passed through a digital S-K filter, a spectral peak may drift to the left. In field data analysis, for a well-calibrated energy spectrum analysis tool, the spectrum peak shift is not beneficial to data analysis.
Disclosure of Invention
The invention aims to solve the technical problem of providing a multi-time bidirectional S-K smoothing processing system and a method, and solves the problem that qualitative and quantitative analysis of low-content nuclides is influenced by characteristics such as electronic noise, statistical fluctuation and the like in the conventional nuclear energy spectrum measurement.
The technical scheme for solving the technical problems is as follows: a multi-pass bi-directional S-K smoothing system comprising:
the first buffer is used for temporarily storing an input data stream when the filtering starts, temporarily storing the data stream after the backward digital SK filter in the filtering process and outputting the data stream after the filtering is finished;
the second buffer is used for temporarily storing the data stream after passing through the forward digital SK filter;
a forward digital SK filter for forward SK filtering;
the reverse digital SK filter is used for reverse SK filtering;
the control register is used for controlling the filtering times and the data flow direction;
the parameter register is used for storing the filtering parameters;
the first buffer, the forward digital SK filter, the second buffer and the backward digital SK filter are sequentially connected, the backward digital SK filter is connected with the first buffer, and the first buffer, the forward digital SK filter, the second buffer and the backward digital SK filter are respectively connected with the control register and the parameter register.
Further, the output signal of the second buffer is input to an inverse digital SK filter after being subjected to direction correction.
A multi-time bidirectional S-K smoothing processing method comprises the following steps:
s1, transmitting the data stream to be filtered to a first buffer for buffering;
s2, filtering the data flow buffered by the first buffer through a forward digital SK filter, and then transmitting the data flow to a second buffer for buffering;
s3, correcting the direction of the data flow cached in the second cache, filtering the data flow through a reverse digital SK filter, and transmitting the data flow to the first cache for caching;
s4, comparing the filtering times with the filtering parameters stored in the parameter register, if the filtering times are less than the filtering parameters, entering the step S2, and if the filtering times are more than the filtering parameters, entering the step S5;
and S5, outputting the data stream buffered in the first buffer.
The invention has the beneficial effects that: the multi-time bidirectional S-K smoothing processing system and the method provided by the invention have the same effect as multi-time Gaussian smoothing processing, compared with Gaussian smoothing, the multi-time bidirectional S-K smoothing processing system and the method provided by the invention have the advantages of higher execution efficiency, simple arithmetic operation process, convenience in use and easiness in software and hardware realization, and can be applied to nuclear instruments and actual measurement nuclear spectrum analysis to solve the problem that qualitative and quantitative analysis of low-content nuclides is influenced by characteristics such as electronic noise, statistical fluctuation and the like in the conventional nuclear spectrum measurement.
Drawings
FIG. 1 is a schematic diagram of the present invention;
FIG. 2 is a flow chart of the present invention;
FIG. 3 is an enlarged view of an area of interest of a measured spectral line of cadmium rice energy dispersion X fluorescence measurement in an embodiment of the present invention;
FIG. 4 is a linear analysis graph of multiple Gaussian smoothing and multiple bidirectional S-K smoothing of spectral lines of interest in an embodiment of the present invention;
FIG. 5 shows the spectral lines of interest of cadmium rice after subtraction of SNIP background in an embodiment of the present invention.
Detailed Description
The principles and features of this invention are described below in conjunction with the following drawings, which are set forth by way of illustration only and are not intended to limit the scope of the invention.
As shown in fig. 1, a multi-pass bidirectional S-K smoothing (abbreviated MBSK) processing system includes:
the first buffer is used for temporarily storing an input data stream when the filtering starts, temporarily storing the data stream after the backward digital SK filter in the filtering process and outputting the data stream after the filtering is finished;
the second buffer is used for temporarily storing the data stream after passing through the forward digital SK filter;
a forward digital SK filter for forward SK filtering;
the reverse digital SK filter is used for reverse SK filtering;
the control register is used for controlling the filtering times and the data flow direction;
the parameter register is used for storing the filtering parameters;
the first buffer, the forward digital SK filter, the second buffer and the backward digital SK filter are sequentially connected, the backward digital SK filter is connected with the first buffer, and the first buffer, the forward digital SK filter, the second buffer and the backward digital SK filter are respectively connected with the control register and the parameter register.
And the output signal of the second buffer is input into the inverse digital SK filter after being corrected according to the data input-first and output-second direction.
As shown in fig. 2, a multiple bidirectional S-K smoothing method includes the following steps:
s1, transmitting the data stream to be filtered to a first buffer for buffering;
s2, filtering the data flow buffered by the first buffer through a forward digital SK filter, and then transmitting the data flow to a second buffer for buffering;
s3, correcting the direction of the data flow cached in the second cache, filtering the data flow through a reverse digital SK filter, and transmitting the data flow to the first cache for caching;
s4, comparing the filtering times with the filtering parameters stored in the parameter register, if the filtering times are less than the filtering parameters, entering the step S2, and if the filtering times are more than the filtering parameters, entering the step S5;
and S5, outputting the data stream buffered in the first buffer.
In the embodiment of the invention, a multi-time bidirectional S-K filtering technology is adopted, and the basic principle is as follows: firstly, forward filtering is carried out, namely nuclear energy spectrum data pass through a digital S-K filter from front to back; and then, performing reverse filtering, namely reversing the energy spectrum data flow after the previous filtering so as to enable the energy spectrum data flow to pass through a digital S-K filter from back to front. The nuclear energy spectrum data after one-time bidirectional SK filtering is operated, the position of a spectrum peak is unchanged, and data analysis is conveniently carried out by using energy scale information of an original energy spectrum.
The energy dispersion X-fluorescence analyzer 2048 is used for measuring the cadmium rice, and after the measurement, the region of interest (1605 and 1665) is intercepted, and the spectral lines are shown in FIG. 3.
In the embodiment of the invention, after proper filtering parameters are set, the effect of multiple bidirectional S-K smoothing processes can be the same as that of multiple Gaussian smoothing processes, and compared with Gaussian smoothing, the method has the advantages of higher execution efficiency, simple arithmetic operation process, convenience in use and easiness in software and hardware realization.
Multiple gaussian smoothing parameters are set as follows: the window width is 7, the standard deviation is 1.25, and the filtering frequency is 55; the filter parameters of the multi-time bidirectional S-K smoothing are as follows: k is 1.02, a is 1, and the filtering times are 10. Using MATLAB2016b, recording the running times of gaussian smoothing and multiple bi-directional S-K smoothing of the above parameters, respectively, we can obtain: gaussian smoothing operation time is 0.0122S, and multi-time bidirectional S-K smoothing operation time is 0.0022S, so that the code execution efficiency of the method is higher. Experiments show that the two smooth data are highly overlapped, multiple bidirectional S-K smoothing of the region of interest and Gaussian smoothed data are subjected to linear fitting, and as shown in FIG. 4, the goodness of fit R can be obtained2Is 0.99999, as shown in the figureAnd 5, respectively subtracting the background of the original measuring spectral line of the cadmium rice and the spectral line subjected to multi-time bidirectional S-K smoothing by using an SNIP method, and then intercepting the spectral line of the region of interest.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.
Claims (2)
1. A multi-pass bi-directional S-K smoothing system, comprising:
the first buffer is used for temporarily storing an input data stream when the filtering starts, temporarily storing the data stream after the backward digital SK filter in the filtering process and outputting the data stream after the filtering is finished;
the second buffer is used for temporarily storing the data stream after passing through the forward digital SK filter;
a forward digital SK filter for forward SK filtering;
the reverse digital SK filter is used for reverse SK filtering;
the control register is used for controlling the filtering times and the data flow direction;
the parameter register is used for storing the filtering parameters;
the first buffer, the forward digital SK filter, the second buffer and the backward digital SK filter are sequentially connected, the output end of the backward digital SK filter is connected with the input end of the first buffer, and the first buffer, the forward digital SK filter, the second buffer and the backward digital SK filter are respectively connected with the control register and the parameter register;
the output signal of the second buffer is input into a reverse digital SK filter after being subjected to direction correction;
the data stream is a nuclear spectrum.
2. The multi-pass bidirectional S-K smoothing method of the multi-pass bidirectional S-K smoothing system according to claim 1, comprising the steps of:
s1, transmitting the data stream to be filtered to a first buffer for buffering;
s2, filtering the data flow buffered by the first buffer through a forward digital SK filter, and then transmitting the data flow to a second buffer for buffering;
s3, correcting the direction of the data flow cached in the second cache, filtering the data flow through a reverse digital SK filter, and transmitting the data flow to the first cache for caching;
s4, comparing the filtering times with the filtering parameters stored in the parameter register, if the filtering times are less than the filtering parameters, entering the step S2, and if the filtering times are more than the filtering parameters, entering the step S5;
s5, outputting the data stream buffered in the first buffer;
the data stream is a nuclear spectrum.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810285498.9A CN108663707B (en) | 2018-04-02 | 2018-04-02 | Multi-time bidirectional S-K smoothing processing system and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810285498.9A CN108663707B (en) | 2018-04-02 | 2018-04-02 | Multi-time bidirectional S-K smoothing processing system and method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108663707A CN108663707A (en) | 2018-10-16 |
| CN108663707B true CN108663707B (en) | 2020-12-08 |
Family
ID=63782744
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810285498.9A Active CN108663707B (en) | 2018-04-02 | 2018-04-02 | Multi-time bidirectional S-K smoothing processing system and method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108663707B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110135402A (en) * | 2019-06-04 | 2019-08-16 | 东华理工大学 | A kind of nuclear spectrum data smoothing method based on multi-parameters model |
| CN113189634B (en) * | 2021-03-02 | 2022-10-25 | 四川新先达测控技术有限公司 | Gaussian-like forming method |
| CN114511477A (en) * | 2022-01-20 | 2022-05-17 | 上海盈方微电子有限公司 | NEON optimization method suitable for Gaussian smoothing |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102302363A (en) * | 2011-06-29 | 2012-01-04 | 深圳市理邦精密仪器股份有限公司 | Device and method for automatically identifying fetal heart rate baseline |
| CN103941279A (en) * | 2014-05-08 | 2014-07-23 | 东华理工大学 | Nuclear pulse signal digital Gaussian forming method based on double parameters of pulse width and amplitude width |
| US8803596B2 (en) * | 2011-09-23 | 2014-08-12 | Tensorcom, Inc. | Differential source follower having 6dB gain with applications to WiGig baseband filters |
| CN105806338A (en) * | 2016-03-17 | 2016-07-27 | 孙红星 | GNSS/INS integrated positioning and directioning algorithm based on three-way Kalman filtering smoother |
| CN105989410A (en) * | 2015-03-05 | 2016-10-05 | 成都理工大学 | Overlap kernel pulse separation method |
-
2018
- 2018-04-02 CN CN201810285498.9A patent/CN108663707B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102302363A (en) * | 2011-06-29 | 2012-01-04 | 深圳市理邦精密仪器股份有限公司 | Device and method for automatically identifying fetal heart rate baseline |
| US8803596B2 (en) * | 2011-09-23 | 2014-08-12 | Tensorcom, Inc. | Differential source follower having 6dB gain with applications to WiGig baseband filters |
| CN103941279A (en) * | 2014-05-08 | 2014-07-23 | 东华理工大学 | Nuclear pulse signal digital Gaussian forming method based on double parameters of pulse width and amplitude width |
| CN105989410A (en) * | 2015-03-05 | 2016-10-05 | 成都理工大学 | Overlap kernel pulse separation method |
| CN105806338A (en) * | 2016-03-17 | 2016-07-27 | 孙红星 | GNSS/INS integrated positioning and directioning algorithm based on three-way Kalman filtering smoother |
Non-Patent Citations (2)
| Title |
|---|
| Improvement of digital S-K filter and its application in nuclear signal processing;Zhou Jianbin et al.;《Nuclear Science and Techniques》;20131231;第24卷(第060401期);第060401-1至060401-7页 * |
| Modeling of Opto-Electronics in Complex Photonic Integrated Circuits;C.Arellano et al.;《PRoc SPIE》;20140331;第1-7页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108663707A (en) | 2018-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108663707B (en) | Multi-time bidirectional S-K smoothing processing system and method | |
| Allen | Short term spectral analysis, synthesis, and modification by discrete Fourier transform | |
| Jordanov et al. | Digital synthesis of pulse shapes in real time for high resolution radiation spectroscopy | |
| JP5111105B2 (en) | Measurement method and analysis apparatus including signal processing including a chain of main impulses | |
| CN110609050B (en) | Method and system for eliminating X-ray fluorescence spectrum peak tailing | |
| US20100305873A1 (en) | Method and Apparatus for Spectral Deconvolution of Detector Spectra | |
| CN111881858B (en) | Microseismic signal multi-scale denoising method and device and readable storage medium | |
| CN101627428A (en) | Noise suppression method, device, and program | |
| CN108983279A (en) | A kind of Low background Spectra Unfolding Methods based on Nal detector | |
| CN105301627B (en) | A kind of energy spectrum analysis method, energy spectrum analysis system and gamma-ray detection system | |
| CN105911351A (en) | DSP-based real time ZFFT method | |
| CN105629290B (en) | A kind of number core pulse signal Mexico hat wavelet manufacturing process | |
| CN107407703A (en) | The spectrum analyzer of real-time bandwidth with raising | |
| Selesnick | Sparse deconvolution (an MM algorithm) | |
| CN104597330A (en) | Electromagnetic radiation signal collecting and processing system and method | |
| CN106841238A (en) | The method and system for positron annihilation dopplerbroadening spectrum based on anticoincidence | |
| CN202771001U (en) | Digitalized [gamma] nuclide identification instrument | |
| CN114441897A (en) | Method for identifying arrival time difference of partial discharge pulse of distribution cable line | |
| Zhou et al. | Study of recursive model for pole-zero cancellation circuit | |
| CN111413725B (en) | System and method for realizing gamma-gamma digital coincidence measurement by using virtual instrument technology | |
| CN109584902B (en) | Music rhythm determining method, device, equipment and storage medium | |
| Un et al. | Hardware support of the periodic component detection algorithm based on zero-crossing | |
| CN115508618B (en) | Quasi-synchronous harmonic analysis device and method based on time domain Hermite interpolation | |
| CN111008356B (en) | Gamma energy spectrum set analysis method for deducting background based on WTS VD algorithm | |
| CN104156194B (en) | Generating method and device for quantum random number |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |