CN104777508A - Digital pulse overlapping peak separation algorithm based on model base - Google Patents
Digital pulse overlapping peak separation algorithm based on model base Download PDFInfo
- Publication number
- CN104777508A CN104777508A CN201510143272.1A CN201510143272A CN104777508A CN 104777508 A CN104777508 A CN 104777508A CN 201510143272 A CN201510143272 A CN 201510143272A CN 104777508 A CN104777508 A CN 104777508A
- Authority
- CN
- China
- Prior art keywords
- pulse
- pulses
- data
- pending
- overlapping
- 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
Abstract
The invention discloses a digital pulse overlapping peak separation algorithm based on a model base. The digital pulse overlapping peak separation algorithm comprises establishment of a pulse model base and processing of digital pulse overlapping peaks. The algorithm comprises basic steps as follows: firstly, a peak height is taken as a classification index in digital pulse data, pulses which have remarkable pulse peaks and are not overlapped are selected to be template pulses, data such as peak heights, pulse widths and the like are extracted, and the pulse model base is established; secondarily, corresponding pulse templates are matched according to to-be-separated overlapping pulse peak heights, and overlapping peaks are separated.
Description
Technical field
The present invention relates to and belong to nuclear spectrum measurement research field, be specifically related to a kind of digit pulse overlap peak separation algorithm based on model bank.
Background technology
In recent years, line Measurement Technique develop rapidly, nuclear spectrum be measured as a kind of conventional line Measurement Technique.Utilize the spectrometer of high count rate to measure nuclear spectrum, need in its process to process the digit pulse after high speed acquisition.In these digit pulses, overlap peak phenomenon clearly, therefore needs to process these overlap peaks, to improve the resolution of surveying instrument.The overlap peak separation algorithm of current existence, the methods such as the matching of overlap peak, iteration, least square that mainly have employed are separated overlap peak, algorithm is comparatively complicated, calculated amount is large, and these algorithms are separated mainly for the overlap peak in spectral line data, and be separated not too applicable for the overlap peak in digit pulse.Invent a kind of digit pulse overlap peak separation algorithm based on model bank for this reason, overlap peak has been processed.
Summary of the invention
For the deficiencies in the prior art, the present invention proposes a kind of digit pulse overlap peak separation algorithm based on model bank, the method is simple, laminating is actual, result and legitimate reading error less, and algorithm automaticity is high, can be effective to the online process of digit pulse overlap peak.
To achieve these goals, the technical solution used in the present invention is:
The present invention is based on the digit pulse overlap peak separation algorithm of model bank, the method comprises following step:
Step one: pending digit pulse data according to condition screened, obtains multiple complete independent pulse, and it can be used as multiple template pulses, and wherein, pending digit pulse data are the digit pulse data comprising overlapping pulses and not overlapping pulses.
Step 2: extract each point data of the pulse height in multiple template pulses, pulse width and negative edge, set up impulse model storehouse.
Step 3: with the pulse height of negative edge overlap in one group of overlapping pulses to be separated in pending digit pulse data for criterion, search in impulse model storehouse, obtains template pulses highly corresponding with it.
Step 4: the pulse reappearing the negative edge overlap described in pending digit pulse by the template pulses that the height obtained in step 3 is corresponding, thus the pulse of isolating described negative edge overlap; Then the rising edge each point ordinate data of the rising edge overlapping pulses overlapped with it and the negative edge each point ordinate data of correspondence position in the negative edge overlapping pulses be separated are made difference to process, and the rising edge of the pulse of the rear data reproduction rising edge overlap obtained of difference process will be made, thus isolate the pulse of rising edge overlap.
Step 5: repeat step 3 to step 4, successively overlapping pulses separation is carried out to overlapping pulses to be separated in pending digit pulse data.
Further, the screening of described step one is specially: the pulse width that the pulse width do not overlapped in given and pending digit pulse data is identical is judge index, using pulses different from given pulse width phase co-altitude in pending digit pulse data as independent pulse, wherein, be not highly all interpulse height relative deviation and be greater than 2%.
Further, in described step 3, obtain highly corresponding being with it specially: the pulse height of the template pulses in the impulse model storehouse pulse height relative deviation overlapping with negative edge is within 2%.
Further, in described step 4, the concrete formula doing difference process employing is:
Y
2’(x)=Y
2(x)-Y
1 (x)（1）
Wherein, certain any horizontal ordinate in x indicating impulse, Y(x) represent that horizontal ordinate is the ordinate of the point of x; Y
2' the ordinate data of (x) to be the rising edge horizontal ordinate of rising edge overlapping pulses obtained after doing difference process be x point, Y
2x rising edge horizontal ordinate that () is the front rising edge overlapping pulses of separation is the ordinate data of x point, Y
1x in the negative edge overlapping pulses after () separation, negative edge horizontal ordinate is the ordinate data of x point; Wherein x represents the time, and unit is 5x10
-9/ s; Y indicating impulse intensity.
Compared with prior art, the invention has the beneficial effects as follows:
1) pulse width that the present invention is based on all independent pulses in the digit pulse overlap peak separation algorithm of model bank is equal, tallies with the actual situation and simplifies calculating process.
2) the present invention is based on the digit pulse overlap peak separation algorithm of model bank by setting up impulse model storehouse, impulse model storehouse can be used for the digit pulse processing procedure of the same type of identical sample frequency, only need continue to expand model bank, wide accommodation.
3) method that the digit pulse overlap peak separation algorithm that the present invention is based on model bank obtains impulse model by search in impulse model storehouse processes overlapping pulses, and method calculating process is simple, and only need subtraction, therefore the degree of accuracy of separating resulting is very high.
4) the present invention is separated overlapping pulses according to impulse model storehouse, and method is simple, and laminating is actual, result and legitimate reading error less, and algorithm automaticity is high, can be effective to the online process of digit pulse overlap peak.
Accompanying drawing explanation
Fig. 1 is the digit pulse overlap peak separation algorithm process flow diagram that the present invention is based on model bank;
Fig. 2 is the pending digit pulse data of the present invention;
Fig. 3 is the present invention's template pulses;
Fig. 4 is one group of overlapping pulses in Fig. 2 of the present invention in pending digit pulse data;
Fig. 5 is the present invention's template pulses corresponding with negative edge overlapping pulses height in Fig. 4 overlapping pulses;
Fig. 6 is two independently pulses that in Fig. 4 of the present invention, overlapping pulses obtain after being separated.
Embodiment
Below in conjunction with accompanying drawing, embodiments of the invention are elaborated:
As shown in Figure 1, the present invention is based on the digit pulse overlap peak separation algorithm of model bank, comprise following step:
Step one: pending digit pulse data according to condition screened, obtains multiple complete independent pulse, and it can be used as multiple template pulses, and wherein, pending digit pulse data are the digit pulse data comprising overlapping pulses and not overlapping pulses.
Step 2: extract each point data of the pulse height in multiple template pulses, pulse width and negative edge, set up impulse model storehouse.
Step 3: in pending digit pulse data with the pulse height of negative edge overlap in one group of overlapping pulses to be separated for criterion, search in impulse model storehouse, obtain template pulses highly corresponding with it.
Step 4: the pulse reappearing the negative edge overlap described in pending digit pulse by the template pulses that the height obtained in step 3 is corresponding, thus the pulse of isolating described negative edge overlap; Then the rising edge each point ordinate data of the rising edge overlapping pulses overlapped with it and the negative edge each point ordinate data of correspondence position in the negative edge overlapping pulses be separated are made difference to process, and the rising edge of the pulse of the data reproduction rising edge overlap obtained after being used as difference process, thus isolate the pulse of rising edge overlap.
Step 5: repeat step 3 to step 4, successively overlapping pulses separation is carried out to overlapping pulses to be separated in pending digit pulse data.
Further, the screening of described step one is specially: the pulse width that the pulse width do not overlapped in given and pending digit pulse data is identical is judge index, using pulses different from given pulse width phase co-altitude in pending digit pulse data as independent pulse, wherein, be not highly all interpulse height relative deviation and be greater than 2%.
Further, in described step 3, obtain highly corresponding being with it specially: the pulse height of the template pulses in the impulse model storehouse pulse height relative deviation overlapping with negative edge is within 2%.
Further, in described step 4, the concrete formula doing difference process employing is:
Y
2’(x)=Y
2(x)-Y
1 (x) （1）
Wherein, certain any horizontal ordinate in x indicating impulse, Y(x) represent that horizontal ordinate is the ordinate of the point of x; Y
2' the ordinate data of (x) to be the rising edge horizontal ordinate of rising edge overlapping pulses obtained after doing difference process be x point, Y
2x rising edge horizontal ordinate that () is the front rising edge overlapping pulses of separation is the ordinate data of x point, Y
1x in the negative edge overlapping pulses after () separation, negative edge horizontal ordinate is the ordinate data of x point; X represents the time, and unit is 5x10
-9/ s; Y shows pulse strength.
Embodiment 1:
The present embodiment take sample frequency as 200MHz, baseline coordinate is 4050 as pending digit pulse data instance, and be separated wherein one group of overlapping pulses, as shown in figures 2-6, concrete steps are as follows:
Step one: as shown in Figure 2, sampling instant is at 3243x5x10
-9s is to 3601 x5x10
-9a pulse in s is independent pulse, and its pulse width is 358x5x10
-9s; By above-mentioned pending digit pulse data with pulse width 358x5x10
-9s is that criterion is screened, and to obtain multiple complete independent pulse, and it can be used as multiple template pulses, Figure 3 shows that one of them template pulses.
Step 2: extract each point data of pulse height, pulse width and negative edge in above-mentioned multiple template pulses, set up impulse model storehouse.
Step 3: as shown in Figure 4, to sampling instant in pending digit pulse partial data at 1x5x10
-9s is to 462 x5x10
-9overlapping pulses in s carry out overlap peak separation, in pending digit pulse data with overlapping pulses to be separated in the height 6033 of pulse of negative edge overlap be criterion, search in impulse model storehouse, obtain template pulses highly corresponding with it, as shown in Figure 5, the pulse height of this highly corresponding template pulses is 6005, be 0.46% with the relative deviation of the height of negative edge overlapping pulses, be less than 2%, think highly equal, the sampling instant therefore in this model pulse and pending digit pulse is at 1x5x10
-9s is to 462 x5x10
-9negative edge overlapping pulses in overlapping pulses in s match.
Step 4: adopt in step 3 the pulse height obtained to be the template pulses of 6005, reappears sampling instant in pending digit pulse at 1x5x10 with it
-9s is to 462 x5x10
-9the pulse of the negative edge overlap in the overlapping pulses in s, thus isolated the pulse of described negative edge overlap.
Then by the sampling instant in pending digit pulse at 1x5x10
-9s is to 462 x5x10
-9the rising edge each point ordinate data of the pulse of the rising edge overlap in the overlapping pulses in s and the negative edge each point ordinate data of correspondence position in the negative edge overlapping pulses be separated are made difference and are processed, and the rising edge of the pulse of the rear data reproduction rising edge overlap obtained of difference process will be made, thus isolate the pulse of rising edge overlap; So far, complete one group of overlap peak in pending digit pulse data and be separated, two independent pulses that these group overlapping pulses obtain after being separated, as shown in Figure 6.
If also need, other group overlapping pulses in pending digit pulse data are separated, repeat step 3 to step 4, complete the separation of overlapping pulses in pending digit pulse data.
Claims (3)
1., based on a digit pulse overlap peak separation algorithm for model bank, it is characterized in that, the method comprises following step:
Step one: pending digit pulse data screened, obtains multiple complete independent pulse, and it can be used as multiple template pulses, and wherein, pending digit pulse data are the digit pulse data comprising overlapping pulses and not overlapping pulses;
Step 2: extract each point data of the pulse height in multiple template pulses, pulse width and negative edge, set up impulse model storehouse;
Step 3: with the pulse height of negative edge overlap in one group of overlapping pulses to be separated in pending digit pulse data for criterion, search in impulse model storehouse, obtains template pulses highly corresponding with it;
Step 4: the template pulses adopting in step 3 the height that obtains corresponding reappears the pulse of the negative edge overlap described in pending digit pulse, thus isolates the pulse of described negative edge overlap; Then the rising edge each point ordinate data of the rising edge overlapping pulses overlapped with it and the negative edge each point ordinate data of correspondence position in the negative edge overlapping pulses be separated are made difference to process, and the rising edge of the pulse of the data reproduction rising edge overlap obtained after being used as difference process, thus isolate the pulse of rising edge overlap;
Step 5: repeat step 3 to step 4, successively overlapping pulses separation is carried out to overlapping pulses to be separated in pending digit pulse data.
2. the digit pulse overlap peak separation algorithm based on model bank according to claim 2, it is characterized in that, the screening of described step one is specially: the pulse width that a pulse width underlapped in given and pending digit pulse data is identical is judge index, using pulses different from given pulse width phase co-altitude in pending digit pulse data as independent pulse, wherein, be not highly all interpulse height relative deviation and be greater than 2%.
3. the digit pulse overlap peak separation algorithm based on model bank according to claim 1, it is characterized in that, in described step 3, obtain highly corresponding being with it specially: the pulse height of the template pulses in the impulse model storehouse pulse height relative deviation overlapping with negative edge is within 2%.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510143272.1A CN104777508B (en) | 2015-03-30 | 2015-03-30 | A kind of digit pulse overlap peak separation algorithm based on model library |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510143272.1A CN104777508B (en) | 2015-03-30 | 2015-03-30 | A kind of digit pulse overlap peak separation algorithm based on model library |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104777508A true CN104777508A (en) | 2015-07-15 |
| CN104777508B CN104777508B (en) | 2017-06-13 |
Family
ID=53619085
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510143272.1A Expired - Fee Related CN104777508B (en) | 2015-03-30 | 2015-03-30 | A kind of digit pulse overlap peak separation algorithm based on model library |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104777508B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105005068A (en) * | 2015-06-25 | 2015-10-28 | 南京瑞派宁信息科技有限公司 | Pulse classification method and pulse classification system |
| CN107884346A (en) * | 2017-10-12 | 2018-04-06 | 浙江全世科技有限公司 | A kind of Overlapped spectral line separation method based on MPT spectroscopic datas |
| CN108398711A (en) * | 2018-01-31 | 2018-08-14 | 成都理工大学 | A kind of pulse recognition method based on the double-deck parameter model |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62180273A (en) * | 1986-02-04 | 1987-08-07 | Kuria Parusu Kk | Pulse height analyzer |
| CN1049562A (en) * | 1989-07-26 | 1991-02-27 | 法国国家电力企业 | Measurement has the method and the application in the γ spectroscopic assay thereof of physical quantity pulsed or that be convertible into pulsed of random character |
| CN1058470A (en) * | 1991-07-02 | 1992-02-05 | 中国建筑材料科学研究院水泥科学研究所 | Novel multi-element analyzer with isotope and X fluorescence |
| CN1427956A (en) * | 2000-03-31 | 2003-07-02 | 爱克斯康特公司 | Spectrally resolved detection of ionizing radiation |
| US20070114427A1 (en) * | 2005-03-11 | 2007-05-24 | Sumitomo Heavy Industries, Ltd. | Radiation detection circuit and apparatus for radiographic examination |
| CN102338880A (en) * | 2011-06-02 | 2012-02-01 | 中国科学技术大学 | Nuclear pulse amplitude digitizing method and system |
| CN103713310A (en) * | 2014-01-02 | 2014-04-09 | 四川新先达测控技术有限公司 | Nuclear pulse signal high speed digital separation method and system and digital spectrometer |
-
2015
- 2015-03-30 CN CN201510143272.1A patent/CN104777508B/en not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62180273A (en) * | 1986-02-04 | 1987-08-07 | Kuria Parusu Kk | Pulse height analyzer |
| CN1049562A (en) * | 1989-07-26 | 1991-02-27 | 法国国家电力企业 | Measurement has the method and the application in the γ spectroscopic assay thereof of physical quantity pulsed or that be convertible into pulsed of random character |
| CN1058470A (en) * | 1991-07-02 | 1992-02-05 | 中国建筑材料科学研究院水泥科学研究所 | Novel multi-element analyzer with isotope and X fluorescence |
| CN1427956A (en) * | 2000-03-31 | 2003-07-02 | 爱克斯康特公司 | Spectrally resolved detection of ionizing radiation |
| US20070114427A1 (en) * | 2005-03-11 | 2007-05-24 | Sumitomo Heavy Industries, Ltd. | Radiation detection circuit and apparatus for radiographic examination |
| CN102338880A (en) * | 2011-06-02 | 2012-02-01 | 中国科学技术大学 | Nuclear pulse amplitude digitizing method and system |
| CN103713310A (en) * | 2014-01-02 | 2014-04-09 | 四川新先达测控技术有限公司 | Nuclear pulse signal high speed digital separation method and system and digital spectrometer |
Non-Patent Citations (1)
| Title |
|---|
| 周星风等: "谱图解析中重叠峰分离的计算方法研究(Ⅰ)", 《计算机与应用化学》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105005068A (en) * | 2015-06-25 | 2015-10-28 | 南京瑞派宁信息科技有限公司 | Pulse classification method and pulse classification system |
| CN105005068B (en) * | 2015-06-25 | 2018-04-06 | 南京瑞派宁信息科技有限公司 | A kind of method and system of pulse classification |
| CN107884346A (en) * | 2017-10-12 | 2018-04-06 | 浙江全世科技有限公司 | A kind of Overlapped spectral line separation method based on MPT spectroscopic datas |
| CN108398711A (en) * | 2018-01-31 | 2018-08-14 | 成都理工大学 | A kind of pulse recognition method based on the double-deck parameter model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104777508B (en) | 2017-06-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103955518A (en) | Method for matching a test object spectrum with database spectrum | |
| CN102565857B (en) | Automatic remaining dynamic correction method | |
| CN104545887B (en) | The recognition methods of artifact ecg wave form and device | |
| WO2013140132A3 (en) | Multi-dimensional survey scans for improved data dependent acquisitions | |
| WO2016141598A1 (en) | Surface consistent amplitude compensation method and device in common attitude gather | |
| CN104155691B (en) | Converted wave anisotropic velocity analysis method and device | |
| CN104777508A (en) | Digital pulse overlapping peak separation algorithm based on model base | |
| WO2012035412A3 (en) | Data independent acquisition of production spectra and reference spectra library matching | |
| CN105334279A (en) | High-resolution mass spectrum data processing method | |
| CN104766341B (en) | A kind of remote sensing image building change detecting method based on DSM and Density Estimator | |
| CN105114067A (en) | Lithology electrofacies method | |
| CN109507709A (en) | Core pulse signal processing method and processing device | |
| Wang et al. | Waveform timing performance of a 5 GS/s fast pulse sampling module with DRS4 | |
| CN104375181B (en) | A kind of method quickly generating three-dimensional common receiver first arrival road collection | |
| CN104182658A (en) | Tandem mass spectrogram identification method | |
| CN110244361A (en) | The seismic data separation method and system handled simultaneously based on interpolation and anti-aliasing | |
| CN103512877B (en) | A kind of Quick lookup method for Raman spectrum substance detection sample | |
| Li et al. | New method of horizon recognition in seismic data | |
| CN106680875A (en) | Analysis method and apparatus of wide-azimuth seismic migration speed | |
| CN105301655B (en) | A kind of common imaging gather line noise minimizing technology and device | |
| CN103364023A (en) | Total least square-based extraction method for central point of instrument image | |
| CN107728168B (en) | Cycle slip detection method based on morphological filtering and singular value decomposition | |
| CN106127154A (en) | The Wavelet Package Denoising Method of partial discharge radio-frequency signal is extracted from white noise | |
| CN102636718B (en) | Mutual capacitance touch screen detection method | |
| CN103837890A (en) | Seismic data acquiring method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170613 |