CN108663707A - A kind of repeatedly two-way S-K smoothing processings system and method - Google Patents
A kind of repeatedly two-way S-K smoothing processings system and method Download PDFInfo
- Publication number
- CN108663707A CN108663707A CN201810285498.9A CN201810285498A CN108663707A CN 108663707 A CN108663707 A CN 108663707A CN 201810285498 A CN201810285498 A CN 201810285498A CN 108663707 A CN108663707 A CN 108663707A
- Authority
- CN
- China
- Prior art keywords
- buffer
- filters
- reversed
- filtering
- repeatedly
- 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
- 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 present invention relates to a kind of repeatedly two-way S-K smoothing processings system and methods, first buffer, forward direction number SK filters, the second buffer and reversed number SK filters are sequentially connected, the reversed number SK filters are connected with the first buffer, and first buffer, forward direction number SK filters, the second buffer and reversed number SK filters are respectively connected with control register and parameter register.It is provided by the invention that repeatedly two-way S-K smoothing processings system and method can be identical as the effect of multiple Gaussian smoothing, compared with Gaussian smoothing, the execution efficiency higher of the present invention, algorithm calculating process is simple, it is easy to use, it is easy to software and hardware realization, present invention can apply to nuclear instrument and actual measurement nuclear energy spectrum analysis, to solve the problem of that characteristics such as Electronics noice and statistic fluctuation in existing nuclear spectrum measurement influence the qualitative and quantitative analysis of low content nucleic.
Description
Technical field
The present invention relates to nuclear spectrum analysis fields, and in particular to a kind of repeatedly two-way S-K smoothing processings system and method.
Background technology
With the progress and development of science, the application field of Nuclear Instruments And Industrial Gauges is also more and more extensive, almost covers agricultural, work
The fields such as industry, health care, geology and mine, environmental protection, space flight, teaching, scientific research, customs.In nuclear decay and nuclear spectrum measurement
Statistic fluctuation characteristic, can make count it is smaller when, real Gaussian Profile desired value can be submerged in noise, though
When counting larger, this statistic fluctuation characteristic also can make measured spectrum data burr phenomena occur, influence later data analysis;For
The presence of reliable recognition low content nucleic characteristic peak, while carrying out qualitative analysis to peak energy amount and position needs research new
Energy spectrum analysis method.
S-K filters propose that low pass S-K filters are in nuclear electronics in nineteen fifty-five by R.P.Sallen and E.L.Key
It is widely used in the filtering forming of core pulse signal, the output of approximate Gaussian waveform can be obtained under less series,
When by nuclear spectrum data back by digital S-K filters, spectral peak can occur to drift about to the left.In live data analysis, for
The good energy spectrum analysis tool of scale, spectral peak drift are unfavorable for data analysis.
Invention content
Technical problem to be solved by the invention is to provide a kind of repeatedly two-way S-K smoothing processings system and methods, solve
The characteristics such as Electronics noice and statistic fluctuation in existing nuclear spectrum measurement, influence asking for the qualitative and quantitative analysis of low content nucleic
Topic.
The technical solution that the present invention solves above-mentioned technical problem is as follows:A kind of repeatedly two-way S-K smoothing processings system, packet
It includes:
First buffer is kept in filtering for the temporary input traffic when filtering beginning through reversed digital
Data flow after SK filters, and output stream after filtering;
Second buffer, for the temporary data flow after positive number SK filters;
Positive number SK filters, for forward direction SK filtering;
Reversed number SK filters, filter for reversed SK;
Control register, for controlling filter times and data flow direction;
Parameter register, for storing filtering parameter;
First buffer, forward direction number SK filters, the second buffer and reversed number SK filters are sequentially connected,
The reversed number SK filters are connected with the first buffer, first buffer, forward direction number SK filters, the second caching
Device and reversed number SK filters are respectively connected with control register and parameter register.
Further, the output signal of second buffer inputs reversed number SK filters after direction is corrected.
A kind of repeatedly two-way S-K smoothing processing methods, include the following steps:
S1, by data stream to be filtered to the first buffer;
S2, by the data flow of the first buffer by being sent to the second buffer after positive number SK filter filterings
Caching;
S3, the data flow cached in the second buffer is corrected into line direction, then passes through reversed number SK filter filterings
After be sent to the first buffer;
S4, filter times are compared with the filtering parameter stored in parameter register, if filter times are less than filtering and join
Number, then enter step S2, if filter times are more than filtering parameter, enters step S5;
S5, the data flow cached in the first buffer is exported.
The beneficial effects of the invention are as follows:It is provided by the invention repeatedly two-way S-K smoothing processings system and method can with it is multiple
The effect of Gaussian smoothing is identical, and compared with Gaussian smoothing, execution efficiency higher of the invention, algorithm calculating process is simple,
It is easy to use, it is easy to software and hardware realization, it is existing to solve present invention can apply to nuclear instrument and actual measurement nuclear energy spectrum analysis
The characteristics such as Electronics noice and statistic fluctuation in nuclear spectrum measurement, the problem of influencing the qualitative and quantitative analysis of low content nucleic.
Description of the drawings
Fig. 1 is principle of the invention figure;
Fig. 2 is flow chart of the present invention;
Fig. 3 is that cadmium rice energy translation measures actual measurement spectral line region of interest enlarged drawing in the embodiment of the present invention;
Fig. 4 is that the multiple Gaussian smoothing of region of interest spectral line and multiple two-way S-K smoothing processings are flat in the embodiment of the present invention
Slip type analysis figure;
Fig. 5 is that the cadmium rice region of interest spectral line after SNIP backgrounds is deducted in the embodiment of the present invention.
Specific implementation mode
The principle and features of the present invention will be described below with reference to the accompanying drawings, and the given examples are served only to explain the present invention, and
It is non-to be used to limit the scope of the present invention.
As shown in Figure 1, a kind of smooth (referred to as MBSK) processing systems of repeatedly two-way S-K, including:
First buffer is kept in filtering for the temporary input traffic when filtering beginning through reversed digital
Data flow after SK filters, and output stream after filtering;
Second buffer, for the temporary data flow after positive number SK filters;
Positive number SK filters, for forward direction SK filtering;
Reversed number SK filters, filter for reversed SK;
Control register, for controlling filter times and data flow direction;
Parameter register, for storing filtering parameter;
First buffer, forward direction number SK filters, the second buffer and reversed number SK filters are sequentially connected,
The reversed number SK filters are connected with the first buffer, first buffer, forward direction number SK filters, the second caching
Device and reversed number SK filters are respectively connected with control register and parameter register.
The output signal of second buffer inputs reversed number SK after the correction of direction first-in last-out by pressing data and filters
Wave device.
As shown in Fig. 2, a kind of repeatedly two-way S-K smoothing processing methods, include the following steps:
S1, by data stream to be filtered to the first buffer;
S2, by the data flow of the first buffer by being sent to the second buffer after positive number SK filter filterings
Caching;
S3, the data flow cached in the second buffer is corrected into line direction, then passes through reversed number SK filter filterings
After be sent to the first buffer;
S4, filter times are compared with the filtering parameter stored in parameter register, if filter times are less than filtering and join
Number, then enter step S2, if filter times are more than filtering parameter, enters step S5;
S5, the data flow cached in the first buffer is exported.
In embodiments of the present invention, using multiple two-way S-K filtering techniques, basic principle is:Make forward filtering first,
Nuclear energy modal data is first passed through into digital S-K filters by sequence from the front to the back;Then inverse filtering is remake, i.e., by previous filter
Gamma-spectrometric data stream after wave is reversed, is allowed to from the front to the back through digital S-K filters.The primary two-way filtered nuclear energy of SK of operation
Modal data, spectrum peak position is constant, facilitates and carries out data analysis with the energy calibration information that proper energy is composed.
Cadmium rice is measured using 2048 X-ray Fluorescence Analyzers, region of interest is intercepted after measurement
(roads 1605-1665), spectral line is as shown in Figure 3.
In embodiments of the present invention, after suitable filtering parameter is arranged, multiple two-way S-K smoothing processings can be with multiple Gauss
The effect of smoothing processing is identical, and compared with Gaussian smoothing, execution efficiency higher of the invention, algorithm calculating process is simple, uses
It is convenient, it is easy to software and hardware realization.
Multiple Gaussian smoothing parameter setting is as follows:Window width is 7, standard deviation 1.25, filter times 55;It is repeatedly double
The filtering parameter smooth to S-K be:K=1.02, a=1, filter times 10.Using MATLAB2016b, record respectively above-mentioned
The Gaussian smoothing of parameter and multiple run time smooth two-way S-K, can obtain:Gaussian smoothing run time 0.0122s, it is repeatedly double
To S-K smoothness run time 0.0022s, it is known that the method for the invention efficiency of code execution higher.Experiment is it is found that two kinds smooth
Data height overlap, the multiple two-way S-K of region of interest is smoothly done into linear fit with data after Gaussian smoothing, such as Fig. 4 institutes
Show, goodness of fit R can be obtained2It is 0.99999, as shown in figure 5, cadmium rice original measurement spectral line and the multiple two-way S-K of warp is smooth
Spectral line afterwards uses SNIP method background corrections, then intercepts the spectral line of region of interest respectively.
The foregoing is merely presently preferred embodiments of the present invention, is not intended to limit the invention, it is all the present invention spirit and
Within principle, any modification, equivalent replacement, improvement and so on should all be included in the protection scope of the present invention.
Claims (3)
1. a kind of repeatedly two-way S-K smoothing processings system, which is characterized in that including:
First buffer is kept in filtering and is filtered through reversed number SK for the temporary input traffic when filtering beginning
Data flow after wave device, and output stream after filtering;
Second buffer, for the temporary data flow after positive number SK filters;
Positive number SK filters, for forward direction SK filtering;
Reversed number SK filters, filter for reversed SK;
Control register, for controlling filter times and data flow direction;
Parameter register, for storing filtering parameter;
First buffer, forward direction number SK filters, the second buffer and reversed number SK filters are sequentially connected, described
Reversed number SK filters are connected with the first buffer, first buffer, forward direction number SK filters, the second buffer and
Reversed number SK filters are respectively connected with control register and parameter register.
2. repeatedly two-way S-K smooth processing units according to claim 1, which is characterized in that second buffer
Output signal inputs reversed number SK filters after direction is corrected.
3. a kind of repeatedly two-way S-K smoothing processing methods, which is characterized in that include the following steps:
S1, by data stream to be filtered to the first buffer;
S2, the data flow of the first buffer is delayed by being sent to the second buffer after positive number SK filter filterings
It deposits;
S3, the data flow cached in the second buffer is corrected into line direction, then by being passed after reversed number SK filter filterings
It send to the first buffer;
S4, filter times are compared with the filtering parameter stored in parameter register, if filter times are less than filtering parameter,
S2 is entered step, if filter times are more than filtering parameter, enters step S5;
S5, the data flow cached in the first buffer is exported.
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 true CN108663707A (en) | 2018-10-16 |
CN108663707B 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) |
Cited By (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 |
CN113189634A (en) * | 2021-03-02 | 2021-07-30 | 四川新先达测控技术有限公司 | 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 |
US20130285746A1 (en) * | 2011-09-23 | 2013-10-31 | 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 |
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 (6)
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 |
US20130285746A1 (en) * | 2011-09-23 | 2013-10-31 | Tensorcom, Inc. | Differential Source Follower having 6dB Gain with Applications to WiGig Baseband Filters |
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 (3)
Title |
---|
C.ARELLANO ET AL.: "Modeling of Opto-Electronics in Complex Photonic Integrated Circuits", 《PROC SPIE》 * |
ZHOU JIANBIN ET AL.: "Improvement of digital S-K filter and its application in nuclear signal processing", 《NUCLEAR SCIENCE AND TECHNIQUES》 * |
纪跃波 等: "零相位数字滤波器", 《重庆大学学报(自然科学版)》 * |
Cited By (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 |
CN113189634A (en) * | 2021-03-02 | 2021-07-30 | 四川新先达测控技术有限公司 | Gaussian-like forming method |
CN114511477A (en) * | 2022-01-20 | 2022-05-17 | 上海盈方微电子有限公司 | NEON optimization method suitable for Gaussian smoothing |
Also Published As
Publication number | Publication date |
---|---|
CN108663707B (en) | 2020-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108663707A (en) | A kind of repeatedly two-way S-K smoothing processings system and method | |
DE112008003554B4 (en) | Multi-channel momentum analyzer for use with multi-channel spectrometers | |
CN106568805A (en) | Highly-integrated Langmuir probe diagnosis system and method | |
CN108983279B (en) | Low-background spectrum resolving method based on sodium iodide detector | |
US5132540A (en) | Dead time correction and regulation method and apparatus for nuclear spectroscopy | |
CN108181648B (en) | The quality control method and device of continuous record aliased seismic acquisition data | |
WO2021174901A1 (en) | Molecular omics data structure implementation method based on data independent acquisition mass spectrum | |
CN108009122B (en) | Improved HHT method | |
CN107219547B (en) | Method for deducting gamma energy spectrum full spectrum background by utilizing successive approximation SNIP (selective non-ip) calculation | |
CN107797154A (en) | A kind of electric imaging logging image scale method and device | |
CN104156565A (en) | System state analysis method and analysis device based on offline CAN bus data | |
CN109634944B (en) | Network loss data cleaning method based on multi-dimensional space-time analysis | |
CN105629290A (en) | Method of forming Mexico hat wavelet by digital nuclear pulse signal | |
Petriska et al. | CDBTools-evaluate positron annihilation coincidence Doppler broadening spectrum | |
CN106656200A (en) | Program counter compression method and hardware circuit thereof | |
CN107220215B (en) | Method for calculating initial parameters of particle spectrum fitting function based on moment estimation method | |
CN111413725B (en) | System and method for realizing gamma-gamma digital coincidence measurement by using virtual instrument technology | |
CN101464420A (en) | Portable high-speed multi-channel energy spectrometer | |
CN109753983A (en) | Image classification method, device and computer readable storage medium | |
CN116319383A (en) | Traffic matrix prediction method, device, computer equipment and storage medium | |
Landsberger et al. | NADA92: An automated, user-friendly program for neutron activation data analysis | |
CN115310472A (en) | Nuclear pulse peak sequence-based one-dimensional convolution neural network nuclide identification method | |
CN109117873A (en) | A kind of user behavior analysis method based on Bayesian Classification Arithmetic | |
CN112035159B (en) | Configuration method, device, equipment and storage medium of audit model | |
EP3676640B1 (en) | Methods and systems for calibration of particle detectors |
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 |