WO2006029475A1 - Method and apparatus for resolving individual signals in detector output data. - Google Patents
Method and apparatus for resolving individual signals in detector output data. Download PDFInfo
- Publication number
- WO2006029475A1 WO2006029475A1 PCT/AU2005/001423 AU2005001423W WO2006029475A1 WO 2006029475 A1 WO2006029475 A1 WO 2006029475A1 AU 2005001423 W AU2005001423 W AU 2005001423W WO 2006029475 A1 WO2006029475 A1 WO 2006029475A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- signal
- data
- signals
- detector
- parameters
- Prior art date
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/16—Measuring radiation intensity
- G01T1/17—Circuit arrangements not adapted to a particular type of detector
- G01T1/171—Compensation of dead-time counting losses
-
- 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/16—Measuring radiation intensity
- G01T1/20—Measuring radiation intensity with scintillation detectors
- G01T1/202—Measuring radiation intensity with scintillation detectors the detector being a crystal
Definitions
- the present invention relates generally to the field of the detection and measurement of radiation and in particular, though not exclusively, to a method and apparatus for the recovery, from a radiation detector, of data affected by pulse pile-up.
- the accurate detection and measurement of radiation is employed in many industries including homeland security, scientific instrumentation, medical imaging and the minerals processing industry. These and other industries use the detection and measurement of radiation for the non-invasive analysis of materials or other specimens. Transmission based imaging, spectroscopic analysis or other modalities can be used to perform such analysis.
- Spectroscopy for example, is commonly used to analyze materials. Knowledge about the material is obtained by analysis of radiation emission from elements within the specimen. This emission of radiation can be stimulated emission due to some form of incident radiation or the result of natural emission from the constituent elements.
- Gamma-ray spectroscopy is a form of spectroscopy in which the emitted electromagnetic radiation is in the form of gamma-rays.
- the detection of the resulting radiation is commonly performed with a scintillation crystal (such as thallium-activated sodium iodide, NaI(TI)), though there are a number of other detector types that can also be used.
- NaI(TI) crystals generate ultra-violet photons pursuant to incident gamma-ray radiation. These photons may then be directed to a photomultiplier tube (PMT) which generates a corresponding electrical signal or pulse.
- PMT photomultiplier tube
- pulse-like signals the shape of which is determined by the incident gamma-ray radiation, the detecting crystal and the PMT.
- the fundamental form of these pulse-like signals is referred to as the impulse response of the detector.
- the output from the photomultiplier is an electrical signal representing the summation of input signals, of determined form, generated in response to discrete gamma rays arriving at the scintillation crystal.
- the analysis is affected by a phenomenon known as pulse pile-up [G. F. Knoll, Radiation Detection and Measurement, 3rd edition, Chapter 17, pp. 632-634, 658 and 659, John Wiley and Sons, New York 2000], whereby multiple gamma-rays arriving more or less simultaneously produce signals which sum together and may be counted as a single signal.
- pulse pile-up G. F. Knoll, Radiation Detection and Measurement, 3rd edition, Chapter 17, pp. 632-634, 658 and 659, John Wiley and Sons, New York 2000
- the magnitude of this combined signal is greater than the individual components, leading to errors in later analysis.
- the energy of an incident gamma-ray may be reflected in the amplitude of the pulse- like signal produced by the detector.
- the presence of specific gamma-ray energies within the detector signal is indicative of particular elements in the material from which gamma-rays originate.
- a more sophisticated approach is to make use of prior knowledge about the profile of a single pulse from the detector or to model mathematically the parameters of a signal. It is then possible in principle to distinguish signals or pulses that originate from a single event from those caused by pulse pile-up.
- signals that depart from the simple profile are selected for subsequent analysis. This analysis involves fitting, via an iterative process, two pulses of varying separation and amplitude.
- a method for resolving individual signals in detector output data comprising determining a signal form of signals present in the data (or the impulse response), making parameter estimates of one or more parameters of the signals, wherein the one or more parameters comprise at least signal temporal position, and determining the energy of the signals from at least the signal form and the parameter estimates.
- 'corrupt data' is interchangeable in this context with 'pulse', as it refers to the output corresponding to individual detection events rather than the overall output signal comprising the sum of individual signals.
- 'signal' is interchangeable in this context with 'pulse', as it refers to the output corresponding to individual detection events rather than the overall output signal comprising the sum of individual signals.
- the temporal position (or timing) of a signal can be measured or expressed in various ways, such as according to the time (or position in the time axis) of the maximum of the signal or the leading edge of the signal. Typically this is described as the arrival time ('time of arrival') or detection time.
- 'detector data' refers to data that has originated from a detector, whether processed subsequently by associated or other electronics within or outside the detector.
- the method may include constructing a model of the data from the parameter estimates, and determining the accuracy of the parameter estimates based on a comparison between the detector output data and the model.
- the signal form may be determined by a calibration process that involves measuring the detector's time domain response to one or more single event detections to derive from that data the signal form or impulse response.
- a functional form of this signal form may then be obtained by interpolating the data with (or fitting to the data) a suitable function such as a polynomial, exponential or spline.
- a filter (such as an inverse filter) may then be constructed from this detector signal form.
- An initial estimate of signal parameters may be made by convolution of the output data from the detector with the filter.
- Signal parameters of particular interest include the number of signals and the temporal position (or time of arrival) of each of the signals.
- the particular signal parameters of interest can then be further refined. Firstly, the estimate of the number and arrival times of signals is refined with the application of peak detection and a threshold. Secondly, knowledge of the number of signals and their arrival time, coupled with the detector impulse response (and hence signal form) makes it possible to solve for the energy parameters of the signals.
- the accuracy of the parameter estimation can be determined or 'validated' by comparing a model (in effect, an estimate) of the detector data stream (constructed from the signal parameters and knowledge of the detector impulse response) and the actual detector output. Should this validation process determine that some parameters are insufficiently accurate, these parameters are discarded. In spectroscopic analysis using this method, the energy parameters deemed sufficiently accurate may be represented as a histogram.
- the method may include making the estimates of signal parameters in accordance with the signal form (i.e. the impulse response of the detector used for generating the signal).
- the method may include determining the signal form by a calibration process including measuring the response of the detector to one or more single detections to derive a data based model of the signal form.
- the method may include obtaining a functional form of the model by interpolating the data with a function to generate the expected signal form.
- the function may be a polynomial, exponential or spline function.
- the method may include designing a filter on the basis of the predetermined form of the individual signals produced by the radiation detector.
- the filter may be, for example, of matched filter or inverse filter form.
- the method includes using convolution of the detector output and filter to make an initial estimate of the signal parameters.
- the method may include refining the estimate of the signal parameters.
- the method may include refining the estimate of signal number with a peak detection process.
- the method may include making or refining the estimate of signal temporal position by application of a peak detection process.
- the method may include refining the estimate of signal energy by solving a system of linear equations, by matrix inversion or by iterative techniques.
- the method includes creating a model of the detector output using the signal parameters in combination with the detector impulse response.
- the method may include performing error detection by, for example, comparing the actual detector output data with the model of the detector output, such as by using least-squares or some other measure of the difference between the data and the model.
- the method may include discarding parameters deemed not sufficiently accurately estimated.
- the method includes presenting all sufficiently accurate energy parameters in a histogram.
- the data may include signals of different forms.
- the method may include determining where possible the signal form of each of the signals.
- the signal form depends on the depth within the detector at which the radiation/detector interaction occurs. In other detectors, the signal form may depend on how much time has elapsed since the previous radiation/detector interaction occurred in the same region of the detector.
- the method includes progressively subtracting from the data those signals that acceptably conform to successive signal forms of a plurality of signal forms, and rejecting those signals that do not acceptably conform to any of the plurality of signal forms.
- the invention provides an apparatus forinstalle pile-up recovery from data comprising a plurality of signals output from a radiation detector. The term 'recovery' is used because data that would otherwise be unusable owing to pile-up is 'recovered' and rendered useable.
- the apparatus of this aspect comprises a processor for receiving the data in digitized form, and is programmed to determine the signal form of each of said signals present in the data, to make parameter estimates of one or more parameters of the signals, and to determine the energy of the signals from at least the signal form and the parameter estimates, wherein the one or more parameters comprise at least signal temporal position.
- the apparatus may include an analog to digital converter adapted to receive the data, to convert the data into digitized form, and forward the data in digitized form to the processor. This would be of particular use where the detector outputs analog data.
- the processor may comprise a field programmable gate array (or an array thereof).
- the processor may comprise a digital signal processor (or an array thereof).
- the processor comprises a field programmable gate array (or an array thereof) and a digital signal processor (or an array thereof).
- the apparatus may include an analog front end that includes the analog to digital converter.
- the apparatus may include an electronic computing device in data communication with the processor, for controlling the processor and for displaying an output of the processor.
- the apparatus may include the radiation detector.
- the apparatus may be, for example, a metal detector, a landmine detector, a medical imaging apparatus, a mineral detection apparatus, an oil well logging apparatus, an unexploded ordnance detector, a cargo screening apparatus, an X-ray fluorescence apparatus or an X-ray diffraction apparatus.
- a method for resolving individual signals in detector output data comprising determining the form of signals present in the data, and making parameter estimates of one or more parameters of the signals from at least the form, wherein the one or more parameters comprise at least signal temporal position.
- a method for pulse pile-up recovery from detector output data comprising determining the form of signals present in the data, making parameter estimates of one or more parameters of the signals, wherein the one or more parameters comprise at least signal temporal position, and determining the energy of the signals from at least the form and the parameter estimates.
- Figure 1 is a view of a gamma-ray spectroscopy apparatus according to an embodiment of the present invention
- Figure 2 is a view of a Sodium Iodide NaI(TI) gamma-ray detector of the apparatus of Figure 1 ;
- Figures 3a, 3b and 3c are graphs illustrating pulse pile-up.
- Figure 4 is a diagram illustrating the mathematical modeling of radiation detection used by the signal processing method embodied in the apparatus of Figure
- Figure 5 is a diagram detailing the mathematical model of radiation detection used by the signal processing method embodied in the apparatus of Figure 1 ;
- Figure 6 is a schematic diagram of the apparatus of Figure 1 ;
- Figures 7a, 7b and 7c are plots of unprocessed digitized data collected directly from the output of the detector of Figure 2 over time ranges of 1000 microseconds, 100 microseconds and 10 microseconds respectively;
- Figure 8 is a schematic representation of the signal processing method for pulse pile-up recovery employed by the apparatus of Figure 1 for analyzing spectroscopic data according to this embodiment of the invention
- Figure 9 is a schematic flowchart of the signal processing method for pulse pile-up recovery employed by the apparatus of Figure 1 for analyzing spectroscopic data according to this embodiment of the invention
- Figures 10a, 10b and 10c are plots of the results at different stages of the signal processing method of Figure 9;
- Figure 11 are plots of gamma-ray spectra for a 137Cs source at various input count rates, processed with the method of Figure 9;
- Figure 12 is a plot of the results of a computer simulation of the signal processing method of Figure 9 prepared using a simulated data set produced by a digital nuclear pulse generator;
- Figure 13 is plot of the performance of the simulation of Figure 12 for a gamma-ray source over a range of count rates
- Figures 14a, 14b, 14c and 14d depict the results of applying the signal processing method of Figure 9 to the output of a 76 mm x 76 mm NaI(TI) gamma-ray detector;
- Figures 15a, 15b, 15c and 15d depict the results of applying the signal processing method of Figure 9 to data collected with a HPGe detector;
- Figures 16a, 16b, 16c and 16d depict the results of applying the signal processing method of Figure 9 to the output of a Xenon gas proportional detector.
- Figure 1 is a schematic view of a gamma-ray spectroscopy apparatus adapted to perform pulse pile-up recovery according to an embodiment of the present invention, with an item to be analyzed.
- the apparatus of Figure 1 includes a neutron generator (10) for generating neutrons for interacting with an item under analysis or specimen (12), and a detector unit (14), in the form of a scintillation based gamma-ray radiation detector, for detecting gamma-ray radiation resulting from the interaction of neutrons and the specimen (12).
- the detector unit includes sensors or detector elements (16) that each has a scintillation crystal (in this example, sodium iodide) coupled to a photomultiplier tube (not shown). It will be appreciated that the apparatus could readily be modified, particularly by substituting a different form of detector unit, to detect other forms of radiation.
- a scintillation crystal in this example, sodium iodide
- the apparatus also includes a signal processing unit (18) that comprises two parts: 1 ) an analog to digital converter which produces a digital output corresponding to the analog output of the detector unit, and 2) a processing unit which implements digital signal processing (DSP) routines in accordance with the invention.
- the electrical output signals of the photomultiplier tubes are connected to the signal processing unit.
- the apparatus also includes cables (20) and a computer (22) for display, the former for coupling the output from the signal processing unit to the computer (22).
- FIG 2 is a view of one of the detector elements (16).
- the illustrated detector element is in the form of a NaI(TI) scintillation based gamma-ray detector, and comprises a cylindrical housing in the form of aluminium body (24) with a NaI(TI) crystal (26) located therein at one (forward) end positioned between an aluminium outer end cap (28) (forward of the NaI(TI) crystal (26)) and an inner optical window (30) (rearward of the NaI(TI) crystal (26)).
- the detector includes a photomultiplier tube (32) rearward of the optical window (30).
- Optical coupling fluid (34) may be used between the NaI(TI) crystal (26) and the optical window (30), and between the optical window (30) and the photomultiplier tube (32).
- a further multiplication stage can be provided by a preamplifier (42).
- a preamplifier (42) In this manner an electrical signal, whose amplitude is proportional to the energy of the incident gamma-rays, is present at the detector output terminals (44) of the detector.
- the detector may additionally include a mu metal magnetic shield (46) located about the sides (48) of the photomultiplier tube (32) and extending forwardly of the photomultiplier tube (32) sufficiently far to surround a portion of the NaI(TI) crystal (26).
- Scintillation detectors of the kind last described have high efficiencies, that is, exhibit a high probability of detecting an incident gamma-ray. However, they also exhibit a relatively long detector response time.
- the detector response time is the time required by the detector to detect an incident gamma-ray and return to a state where the next incident gamma-ray can be accurately detected. Radiation detectors with long detector response times are thus prone to pulse pile-up. That is, the output, which ideally consists of completely discrete pulses each corresponding to the incidence of a single gamma-ray, instead exhibits a waveform in which individual pulses can overlap making them difficult to characterize.
- Figures 3a, 3b and 3c illustrate the effect of pulse pile-up, and show illustrative signals or pulses plotted as energy E versus time if (both in arbitrary units).
- Figure 3a illustrates so-called 'tail-end pile-up' where, depending on the type of pulse conditioning employed, the tail (50) of one pulse (51 ) can provide a significant positive or negative bias (positive in the illustrated example) to the amplitude of a subsequent pulse (52).
- the time displacement between the two pulses, N is relatively large (when compared with the overall time interval for which the pulses prevail), the signal envelope or resultant waveform (54) is significantly above zero at the arrival of the second pulse (52).
- FIG 3b illustrates another form of pulse pile-up, 'peak pile-up'.
- two pulses (56) and (58) arrive closely spaced in time (i.e. the time displacement N between the pulses is small compared with the overall time interval over which the pulses prevail).
- the resultant output waveform (60) appears more or less as a single pulse of somewhat greater amplitude than either of the component pulses.
- Figure 3c Multiple signals or pulses (such as those shown at 62) arrive with random time separation ⁇ t and sum to produce a resultant waveform (64) from which the parameters of the component signals are difficult to extract.
- One component of the method of addressing pulse pile-up is the estimation of certain parameters of the signals or pulses; these parameters are the number, time-of-arrival and energy of all gamma-rays in the detector data stream. These parameters are estimated, according to this embodiment, by modeling the signals in the data stream mathematically.
- the model employed in this embodiment includes certain assumptions about the data and the apparatus, as are discussed below.
- Figure 4 is a diagram that illustrates the modeling of the radiation detection process.
- the radiation g(t) (70) is incident on the detector (72) represented by the measurement process m(t), resulting in output data from the detector y(t) (74).
- the addition of a sampling process (76) produces the digital detector data or 'time-series' x[n] (78).
- Figure 5 illustrates a more detailed mathematical mode! of the detection process shown in Figure 4.
- the input g(t) to the detector is characterized by Equation 1 , in which the input g(t) is assumed to be an unknown number [N) of delta-function-like impulses of random amplitude ( ⁇ ) and time of arrival ( ⁇ ).
- N an unknown number of delta-function-like impulses of random amplitude ( ⁇ ) and time of arrival ( ⁇ ).
- An illustrative example of such input is shown at (80).
- the radiation detector is assumed to have a specific response to the incoming radiation, referred to as the detector impulse response d(t) (or, equivalent ⁇ , the signal form of the signals in the data), which is illustrated at (82).
- the digitized version of the detector impulse response i.e. signal form
- the output from the detector is shown at (86) and characterized by Equation 2, in which the detector output y ⁇ t) is the sum of an unknown number of signals of predetermined signal form d(t), with unknown energy ⁇ a) and unknown time of arrival ( ⁇ ). Sources of random noise ⁇ (t) (84) are also considered.
- the digital detector data x[n] (88) is produced by the analog to digital converter (76).
- the digitized signal x[n] (which constitutes a time series of data) at the output of the analog to digital converter (76), as illustrated at (88), is therefore given by
- d[n] is the discrete time form of the signal form Qf(O, ⁇ / is the delay in samples to the /th signal
- ⁇ [n] is the discrete time form of the noise.
- the digitized signal x[n] may also be written in matrix form as
- A an M x N matrix, the entries of which are given by f ⁇ ljn - ⁇ ,] ⁇ , ⁇ n ⁇ min(/W, ⁇ ,. + 7 - 1)
- T is the length of cf[n] in samples
- M is the total number of samples in the digitized signal x[n]
- a is the vector of N signal energies
- ⁇ is the noise vector of length M.
- Matrix A may also be depicted as follows:
- the columns of matrix A contain multiple versions of the signal form.
- the starting point of the signal form is defined by the signal temporal position. For example, if the signals in the data arrive at positions 2, 40, 78 and 125, column 1 of matrix A will have '0' in the first row, the 1st datum point of the signal form in the second row, the 2nd datum point of the signal form in the 3rd row, etc.
- the second column will have '0' up to row 39 followed by the signal form.
- the third column will have '0' up to row 77; the fourth column will have '0' up to row 124 and then the signal form.
- the size of matrix A is determined by the number of identified signals (which becomes the number of columns), while the number of rows depends on the number of samples in the time series.
- the signal processing method of this embodiment thus endeavors to provide an accurate estimate of some unknown parameters of the detector data, including not only the number of component signals [N) in the detector output but also the energy (o) and time-of-arrival ( ⁇ ) of each of the component signals.
- FIG. 6 is a schematic diagram of the functional elements of the gamma-ray spectroscopy apparatus of Figure 1 , and is provided to explain in more detail the signal processing method for pulse pile-up recovery employed by the apparatus of Figure 1.
- the radiation detector unit (14) is connected to a pulse processing board (92) via an analog front end (AFE 94).
- the purpose of the AFE (94) is to digitize the signal produced by the radiation detector unit (14) by performing analog to digital conversion at, in this embodiment, 125 MHz with 12-bit conversion accuracy.
- Figures 7a, 7b and 7c illustrate the waveform resulting from such digitization, over time ranges of 1000 microseconds, 100 microseconds and 10 microseconds respectively.
- the various peaks in these figures correspond to the detection of respective gamma-rays. Some peaks appear as discreet signals or pulses (110, 112) which may indicate the presence of only a single gamma-ray. Other peaks are due to the pile-up either of two peaks (116) or of three or more peaks (114).
- the signal processing method for pulse pile-up recovery is implemented.
- the digital signal produced by the AFE (94) is passed into the pulse processing Field Programmable Gate Array (FPGA) (96).
- the pulse processing FPGA (96) then implements the pulse processing method of this embodiment; a digital signal processing coprocessor (98) may optionally be used to assist the pulse processing FPGA (96) to implement the pulse processing method.
- Variables required by the pulse processing FPGA (96) and data produced at interim steps of the pulse processing method are optionally stored in memory (100).
- the signal processing is controlled via a Data/Control Interface (102) which, in conjunction with a Control Processor (104), can be used to modify the implementation of the signal processing.
- the output data from the signal processing method can be displayed on a display (106) via the Data/Control Interface (102).
- Display (106) is provided in a computer that may, if desired, be used to perform post-processing and system control.
- Figure 8 is a schematic diagram of the signal processing method for pulse pile-up recovery of radiation signals in the detector time series of this embodiment.
- the digitized detector signal (from AFE (94)) forms the input (120) for this signal processing method.
- Offline System Characterization (122) is used to determine the detector impulse response unique to the particular digitized detector signal.
- Characterization data generated in System Characterization phase (122) is use in a Pulse Localization phase (124).
- the Pulse Localization phase (124) estimates, in real-time, the number and temporal position (or time-of-arrival) of radiation pulses within the digitized detector signal.
- a Pulse Identification phase (126) the digitized detector signal, the detector impulse response and the output from the Pulse Localization phase (124) are used to determine the energy of the signals or pulses.
- Validation (128) involves comparing the output of the Pulse Identification phase (126) with the digitized detector signal (120). If this comparison indicates that any of the pulse parameters have been estimated inaccurately, those parameters are rejected so that only valid data is output (130).
- the error signal generated in the Validation phase (128) is also employed in System Characterization (122).
- System Characterization (122) updates the detector impulse response online and adaptively by employing the error signal.
- Such updating of the detector impulse response may be performed with any suitable adaptive method, such as least mean squares adaptation, normalized least mean squares adaptation or recursive least squares adaptation as described, for example, by S. Haykin [Adaptive Filter Theory, 4th Ed, Prentice Hall, 2002].
- FIG. 9 is a flow diagram of the signal processing method of this embodiment.
- calibration is performed. This involves Data Regularization or Conditioning (142), Data Selection and Fitting (144) and Optimal Filter Construction (146).
- Data Regularization 142
- calibration data signals recorded at a low incident radiation flux
- Data Selection and Fitting 144
- selecting only that data corresponding to the detection of single radiation events and constructing a data based model of the detector impulse response.
- a functional form of this model is then obtained by fitting a suitable function to the data, such as a polynomial, exponential or spline function. This results in the expected impulse response of the detector d[n].
- Optimal Filter Construction employs this detector impulse response to construct a suitable filter for the detector, such as an inverse filter or a matched filter.
- step (150) data is acquired, but may be affected by significant pulse pile-up.
- the data may be input (152) either from a file or directly from the detector elements (16).
- step (160) signal processing routines are applied to determine the amplitude and timing parameters of the signals in the time series. Firstly the data is conditioned (162) to remove any bias in the baseline of the data. Next, the detector data is convoluted (164) with the filter derived in step (146) to provide an initial estimate of the time-of-arrival parameters ( ⁇ ) and number of pulses (N). The timing parameters and estimate of the number of pulses are then further refined (166) using a suitable peak detection process, and the energy parameter (a) is determined from ⁇ , N and the detector impulse response d[n] (such as by linear programming, matrix inversion or convolution techniques). Finally, from the number (N), energy (a), timing ( ⁇ ,) and detector impulse response ( ⁇ f[n]), an estimate of the detector data stream (x [n]) is made (168).
- the parameter vector (a) may be determined by linear programming or by solving the system of linear equations defined in Equation 4 using a suitable method for solving such systems of equations, such as one of those described, for example, by G. H. Golub and CF. Van Loan [Matrix Computations, 2nd Ed, Johns Hopkins University Press, 1989].
- the validation phase (128) referred to above is performed, which may be referred to as error checking as, in this embodiment, validation involves determining an error signal e[n], computed successively for the set of samples corresponding to each signal / where 1 ⁇ i ⁇ N (N being the total number of signals in the data stream).
- This error signal is calculated by determining (172) the squares of the differences between the time series data x[n] and the model based data-stream ( x[n] from step (168)); e[n] is thus the square of the difference between x[n] and x[n] , as given in Equation 6.
- e[n] exceeds a predetermined threshold, these parameters are rejected (174) as this condition indicates that the signal parameters do not produce a model of the respective signal that acceptably conforms to that signal (that is, is sufficiently accurate); the relevant signal is deemed to constitute corrupted data and excluded from further spectroscopic analysis.
- the threshold may be varied according to the data and how closely it is desired that the data be modeled; generally, therefore, in any particular specific application, the method of validation and definition of the threshold are chosen to reflect the requirements of that application.
- a threshold is the signal energy a,- multiplied by a suitable factor, such as 0.05.
- a suitable factor such as 0.05.
- Validation may be performed by defining the error signal and threshold in any other suitable way.
- the error signal may be set to the absolute value of the error.
- the threshold may be defined to be a multiple other than 0.05 of the signal amplitude.
- Another threshold comprises a number of noise standard deviations.
- Decreasing the threshold (such as by decreasing the coefficient of ⁇ ,- in Equation 7) enables improved energy resolution at lower throughput, while increasing the threshold enables improved throughput at reduced energy resolution.
- step (180) a decision is made as to whether there is sufficient data. If not, processing continues at step (150). Otherwise, the method proceeds to step (190).
- step (190) a gamma-ray energy spectrum is created.
- the gamma-ray energy parameters determined at step (166), which were deemed to be of sufficient accuracy at step (174), are represented (192) in the form of a histogram. This is the gamma-ray energy spectrum on which spectroscopic analysis may be performed.
- Figures 10a, 10b and 10c are plots of the results at various stages of processing of the digital signal processing method described above by reference to Figures 8 and 9, for digitized data collected with a scintillation gamma-ray detector.
- the detector data stream was digitized by an analog to digital converter at 125 MHz and 12 bit accuracy; the gamma-ray source used was a 137Cs source with a primary gamma- ray emission of 661.7 keV.
- Scintillation detectors employ light generated by the detector/radiation interaction to detect and measure that incident radiation.
- a scintillation detector may comprise organic scintillators or inorganic scintillators.
- Organic scintillators include both organic crystalline scintillators and liquid organic solutions (where the scintillating material has been dissolved to form a liquid scintillator, which can then be plasticized to form a plastic scintillator.
- Inorganic scintillators include crystalline scintillators such as NaI(TI), BGO, CsI(TI) and many others, and photo switch detectors (in which a combination of two or more dissimilar scintillators are optically coupled to a common PMT to exploit the differing decay times of the scintillators to determine where a radiation/detection interaction has occurred).
- the detector comprised a 76 mm x 76 mm NaI(TI) gamma-ray scintillation detector.
- Figure 10a is a plot of a portion of the digitized detector data (200) prior to processing by the signal processing method plotted as energy E(keV) versus time t( ⁇ s), together with the results (for example, at 210) of the signal processing method plotted in terms of the temporal position and energy of the component signals. For example, what may appear to be a single peak (220) in the original digitized detector data (200) at approximately 75.8 ⁇ s has been resolved into two distinct signals (222, 224) at respectively 75.3 and 75.7 ⁇ s.
- Figure 10b is a plot of the resulting data model (230), shown as energy E(keV) versus time t( ⁇ s), of that portion of the digitized detector data stream (200) shown in Figure 10a.
- An inverted error plot (240) comprising a plot of the squares of the differences between the detector data (200) and the data model (230), is also shown, and indicates the error in the model (230).
- the error signal is small where the model (230) has tracked the output of the detector accurately, but the error becomes large when there are inconsistencies between the model (230) of the detector data and the detector data (200) itself. Based on this error signal (240), a decision can be made as to whether to accept or reject the signal parameters estimated by the signal processing method.
- Figure 10c is a gamma-ray energy spectrum (250), shown as a log-linear plot, produced by the signal processing method.
- the energy parameters that have been accepted are plotted as a histogram, where the horizontal axis represents the energy E(keV) of each signal in a respective bin, and the vertical axis represents the number of counts N of that energy determined to have been detected in the collection period (in this example, 1 s).
- Figure 11 is a plot of exemplary gamma-ray energy spectra, collected using a sodium iodide NaI(TI) gamma-ray detector. The gamma-ray energy spectra shown in Figure 11 demonstrate the performance of the signal processing method for pulse pile-up recovery at a range of count rates.
- the experimental data were collected using a 76 mm x 76 mm Canberra brand NaI(TI) gamma-ray detector (model number 802) coupled to a detector base (model number 2007); no preamplifier was used.
- the signal processing hardware was connected to the dynode output of the detector base via a 65 MHz 14-bit analog to digital converter.
- the NaI(TI) crystal was irradiated with a collimated gamma-ray beam, which ensured that the central portion of the detector was illuminated with an essentially parallel beam of gamma-rays; the beam diameter was 50 mm.
- the spectra (260), (262), (264), (266), (268) and (270) were collected at count rates of respectively 529 kHz, 230 kHz, 167 kHz, 124 kHz, 67 kHz and 9 kHz.
- the energy resolution of the data collected with the apparatus and processed with the method of this embodiment deteriorated as the count rate increased.
- the full width at half maximum (FWHM) of the peak was found to be, respectively, 9.6% 7.3%, 7.1 %, 6.9%, 6.7% and 6.7%.
- the energy resolution of the 137Cs gamma-ray energy peak at 661.7 keV remained less than 7.5%; that is, despite more than a 25 fold increase in the count rate from the NaI(TI) detector, the energy resolution at 661.7 keV decreased by less than 0.5%.
- Figure 12 is a log-log plot of the throughput of the signal processing method (i.e. that portion of the input count rate accurately detected) against input count rate from 0.1- 2.5 MHz.
- the theoretical limit i.e. where the throughput equals the input
- This figure demonstrates that, over a very wide range of input count rates, the throughput of the signal processing method remains greater than or equal to 90%.
- Figure 13 is a linear-log plot comparable to Figure 12 but with percentage throughput plotted against input count rate from 0.005-10 MHz.
- Figure 13 includes plots of the energy resolution and peak position performance of the signal processing method of this embodiment. The energy resolution of the 137Cs peak degrades by less than 10% over 0-2.5 MHZ, and the peak position shows very little change over that range.
- Figures 14a, 14b, 14c and 14d also depict the results of applying the signal processing method for pulse pile-up recovery of this embodiment to the output of a 76 mm x 76 mm NaI(TI) gamma-ray detector. Approximately 14 ⁇ s of data was used to generate the data plotted in these figures. The figures are plots of energy E in arbitrary units against time t( ⁇ s).
- Figure 14a is a plot of the output of the AFE (94): an analog to digital conversion rate of 65 MHz and 14 bit resolution was used to covert the time varying voltage output of the detector to digital data.
- Figure 14b is a plot of the results of applying the method. The temporal positions of the signals (depicted as vertical lines) have been resolved, as have the energies of the component signal (depicted as crosses). The temporal position and the energy of the component signal were used as described above, in conjunction with the signal form, to determine a model of the gamma-ray detector output: the resulting model is plotted in Figure 14c.
- the digitized output of the gamma-ray detector was compared with the model of the gamma-ray detector output to derive an estimate of the error made in characterizing the gamma-ray detector output.
- This error signal is plotted in Figure 14d. It is then possible, on the basis of this error signal, to determine thresholds for the exclusion of signal parameter estimates, such as the decision to accept or reject an estimate of signal energy may be determined by the magnitude or the error near the position of a signal peak.
- Figures 15a, 15b, 15c and 15d depict the results of applying the signal processing method for pulse pile-up recovery of this embodiment to data collected with a semiconductor (or solid state) detector.
- a semiconductor or solid state detector.
- Such detectors employ the interaction of incident radiation with the electrons in the crystalline lattice of the semiconductor, forming electron hole pairs. Examples of these detectors include High-Purity Germanium (HPGe) detectors, Silicon Diode detectors, semiconductor drift detectors (such as Silicon Drift detectors), Cadmium Telluride (CdTe) detectors and CZT detectors.
- HPGe High-Purity Germanium
- Silicon Diode detectors Silicon Diode detectors
- semiconductor drift detectors such as Silicon Drift detectors
- Cadmium Telluride (CdTe) detectors Cadmium Telluride
- Figure 1 is a plot of the output of the AFE (94).
- Figure 15b is a plot of the results of applying the method.
- the temporal positions of the signals (depicted as vertical lines) have been resolved, as have the energies of the component signal (depicted as crosses).
- the temporal position, the energy of the component signal and the signal form were used to determine a model of the processed HPGe detector output, which is plotted in Figure 15c.
- Figure 15d is a plot of the error signal, derived from a comparison of the digitized processed output of the HPGe detector and the model of that output. This error signal can again be used to determine thresholds for the exclusion of signal parameter estimates.
- Figures 16a, 16b, 16c and 16d depict the results of applying the signal processing method for pulse pileup recovery of this embodiment to the output of a gas proportional detector used for detecting X-rays.
- Gas proportional detectors are a class of detector whose behavior is similar to that of solid state detectors. Gas proportional detectors rely on the interaction of the radiation with a gas in a chamber.
- An electric field is created in the chamber between an axial wire and the walls of the chamber. Radiation passing through the gas ionizes the gas, which produces electrons that then collect on the wire owing to the electric field, and are output as the detector data.
- the apparatus of Figure 1 was employed, though with a detector unit in the form of a Xenon gas proportional detector substituted for the detector unit (14), and with an X-ray generator from an X-ray diffraction apparatus rather than a neutron source and specimen.
- a detector unit in the form of a Xenon gas proportional detector substituted for the detector unit (14)
- an X-ray generator from an X-ray diffraction apparatus rather than a neutron source and specimen.
- Approximately 300 ⁇ s of data was used to generate the data plotted in Figures 16a, 16b, 16c and 16d, which plot energy E in arbitrary units against time if( ⁇ s).
- a significantly longer data collection period was used compared with that of the previous examples, owing to the relatively long decay time of the xenon gas proportional detector (of the order of 50 ⁇ s or more). For this reason also the sampling rate of the AFE (94) was reduced.
- Figure 16a is a plot of the output of the AFE (94); in this example an analog to digital conversion rate of 15 MHz and 14 bit resolution was used to covert the time varying voltage output of the detector to digital data.
- Figure 16b is a plot of the results of applying the method. The temporal positions of the X-ray signals (depicted as vertical lines) have been resolved, as have the energies of the component signal (depicted as crosses). The temporal position and the energy of the component signal were used as described above, in conjunction with the signal form, to determine a model of the Xenon gas proportional detector output the resulting model is plotted in Figure 16c.
- the digitized output of the Xenon gas proportional detector was compared with the model of the Xenon gas proportional detector output to derive an estimate of the error made in characterizing the Xenon gas proportional detector output.
- This error signal is plotted in Figure 16d. This error signal can then be used to determine thresholds for the exclusion of signal parameter estimates, such as the decision to accept or reject an estimate of signal energy may be determined by the magnitude or the error near the position of a signal peak.
- the form of a given signal may be one of a plurality of possible signal forms. This may be intrinsic to the detector type, or be due to temperature or other measurement-specific factors.
- a CsI(TI) detector is a scintillation detector that, depending on whether a neutron or gamma-ray is being detected, exhibits two distinct signal forms.
- Solid state radiation detectors can exhibit a time-varying signal form, even when detecting only one form of radiation;
- large volume High Purity Germanium (HPGe) detectors for example, can produce an output signal whose form depends on the specific site of interaction between the radiation and the detector.
- HPGe High Purity Germanium
- the electrons and holes travel through the HPGe detector at different velocities, so the charge pulse produced by the electrons generally has a different form from that produced by the holes.
- the pulse produced by the detector (being the sum of the charges carried by both the electrons and holes) has a form dependent on the location of interaction.
- the plurality of signal forms are the result of these varied physical mechanisms.
- the respective signal forms may be denoted di[n], c ⁇ M cfQ[n], where Q is the total number of different signal forms that may be generated by a particular detector type.
- Each of the possible signal forms is characterized in the same way that the signal form of data having a single signal form is characterized.
- the calibration process must be extended for an appropriate length of time to ensure that all of the possible signal forms have been identified and characterized; the estimation of signal parameters, including temporal position and signal energy, can be performed once the form of each signal in the data stream has been identified.
- a number of possible extensions of the method described above may be employed.
- the signal parameters may be estimated for each signal in the data stream by treating all signals in the data stream as having the same form, such as of the first signal, viz. c/i[n].
- the parameters for those signals that do not acceptably conform to signal form c/i[n] are rejected at the validation phase; signals for which the parameters have been estimated successfully and thus acceptably conform to signal form di[n] are subtracted from the data stream.
- This process is repeated successively for cfeM up to CIQ[VI], where at each stage signal parameters are estimated for signals that are of the signal form used at that stage.
- matrix Equation 4 is solved with matrix A constructed repeatedly using, in iteration p, the signal form ⁇ f p [n].
- those signals that have not passed the validation phase for any of the plurality of signal forms are rejected as not acceptably conforming to any of the plurality of signal forms.
- the signal parameters are estimated for each of the signal forms in turn, but the signal estimates are not subtracted at each stage. Instead, the estimated signals are used in a final signal validation stage to determine the signal form and signal parameters that provide the best overall estimate of the data stream. This allows for the possibility that a signal is incorrectly estimated to be of one form, when it is actually of a form that has not yet been used to estimate the signal parameters.
- each of the signal forms d p [n] as a linear combination of two signal forms, termed of ⁇ n] and d 2 [n] for convenience.
- Equation 5 where the sub-matrices A 1 and A 2 are formed from the signal forms di[n] and c/ 2 [n] respectively using Equation 5.
- the new system of linear equations is solved using the same methods as those used to solve the earlier matrix equation, Equation 4.
- this approach eliminates the need for explicitly estimating the unknown constants a and b, and also allows for the possibility that the signal form may be from a continuum of possible signal forms that can be represented as a linear combination of the two signal forms c/- ⁇ [n] and CZ 2 [Yi].
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Chemical & Material Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Measurement Of Radiation (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
Description
Claims
Priority Applications (8)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CA2621736A CA2621736C (en) | 2004-09-16 | 2005-09-16 | Method and apparatus for resolving individual signals in detector output data |
AU2005284602A AU2005284602B2 (en) | 2004-09-16 | 2005-09-16 | Method and apparatus for resolving individual signals in detector output data. |
EP05779842A EP1794884A4 (en) | 2004-09-16 | 2005-09-16 | Method and apparatus for resolving individual signals in detector output data. |
CN2005800392488A CN101124725B (en) | 2004-09-16 | 2005-09-16 | Method and apparatus for resolving individual signals in detector output data |
JP2007531539A JP5022902B2 (en) | 2004-09-16 | 2005-09-16 | Apparatus and method for separating individual signals in detector output data |
US11/685,719 US7383142B2 (en) | 2004-09-16 | 2007-03-13 | Method and apparatus for resolving individual signals in detector output data |
US12/130,671 US20090048795A1 (en) | 2004-09-16 | 2008-05-30 | Method and apparatus for resolving individual signals in detector output data |
US13/283,410 US10310099B2 (en) | 2004-09-16 | 2011-10-27 | Method and apparatus for resolving individual signals in detector output data |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AU2004905364 | 2004-09-16 | ||
AU2004905364A AU2004905364A0 (en) | 2004-09-16 | Method and apparatus for pulse pile-up recovery |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/685,719 Continuation US7383142B2 (en) | 2004-09-16 | 2007-03-13 | Method and apparatus for resolving individual signals in detector output data |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2006029475A1 true WO2006029475A1 (en) | 2006-03-23 |
Family
ID=36059646
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/AU2005/001423 WO2006029475A1 (en) | 2004-09-16 | 2005-09-16 | Method and apparatus for resolving individual signals in detector output data. |
Country Status (6)
Country | Link |
---|---|
US (3) | US7383142B2 (en) |
EP (1) | EP1794884A4 (en) |
JP (1) | JP5022902B2 (en) |
CN (1) | CN101124725B (en) |
CA (1) | CA2621736C (en) |
WO (1) | WO2006029475A1 (en) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007327902A (en) * | 2006-06-09 | 2007-12-20 | Shimadzu Corp | Signal processing device for x ray analysis |
WO2009121130A1 (en) * | 2008-03-31 | 2009-10-08 | Southern Innovation International Pty Ltd | Screening method and apparatus |
WO2009121131A1 (en) * | 2008-03-31 | 2009-10-08 | Southern Innovation International Pty Ltd | Method and apparatus for borehole logging |
WO2009121132A1 (en) * | 2008-03-31 | 2009-10-08 | Southern Innovation International Pty Ltd | Radiation imaging method with individual signal resolution |
GB2460345A (en) * | 2008-05-30 | 2009-12-02 | Precision Energy Services Inc | Borehole logging using LaBr3 scintillation detector and digital spectrometer which rejects pile-up pulses |
WO2010068996A1 (en) | 2008-12-18 | 2010-06-24 | Southern Innovation International Pty Ltd | Method and apparatus for resolving piled-up pulses by using a mathematical transform |
JP2011503550A (en) * | 2007-11-02 | 2011-01-27 | ユニヴァーシティ オブ ワシントン | Data acquisition for positron emission tomography |
JPWO2011039819A1 (en) * | 2009-10-01 | 2013-02-21 | 株式会社島津製作所 | Pulse height analyzer and nuclear medicine diagnostic apparatus provided with the same |
US8581200B2 (en) | 2006-11-17 | 2013-11-12 | Koninklijke Philips N.V. | Radiation detector with multiple electrodes on a sensitive layer |
WO2015085372A1 (en) * | 2013-12-11 | 2015-06-18 | Southern Innovation International Pty Ltd | Method and apparatus for resolving signals in data |
WO2016082006A1 (en) | 2014-11-30 | 2016-06-02 | Southern Innovation International Pty Ltd | Device and method for material characterisation |
JP2017086901A (en) * | 2015-11-02 | 2017-05-25 | 東芝メディカルシステムズ株式会社 | Data collection device, x-ray ct apparatus, and nuclear medicine diagnostic apparatus |
WO2017205914A1 (en) | 2016-05-30 | 2017-12-07 | Southern Innovation International Pty Ltd | Material characterisation system and method |
WO2017205912A1 (en) | 2016-05-30 | 2017-12-07 | Southern Innovation International Pty Ltd | Material characterisation system and method |
US9885674B2 (en) | 2009-12-09 | 2018-02-06 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Method and device for recognition of a material making use of its transmission function |
EP3546926A1 (en) * | 2018-03-28 | 2019-10-02 | PicoQuant Innovations GmbH | Method and device for determining photon arrival times by means of machine learning |
US10485994B2 (en) | 2015-03-30 | 2019-11-26 | Sumitomo Heavy Industries, Ltd. | Charged particle beam treatment apparatus |
Families Citing this family (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8304740B1 (en) * | 2008-05-19 | 2012-11-06 | Innovative American Technology, Inc. | Mobile frame structure with passive/active sensor arrays for non-invasive identification of hazardous materials |
DE10301941B4 (en) * | 2003-01-20 | 2005-11-17 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Camera and method for optical recording of a screen |
US7755015B2 (en) * | 2004-10-22 | 2010-07-13 | Lawrence Livermore National Security, Llc | Real-time multiplicity counter |
US8384004B1 (en) * | 2004-10-22 | 2013-02-26 | Lawrence Livermore National Security | Real-time multi-mode neutron multiplicity counter |
US20100226580A1 (en) * | 2005-12-01 | 2010-09-09 | Innovative American Technology Inc. | System and method for increased gamma/neutron detection |
US8247781B2 (en) | 2005-12-01 | 2012-08-21 | Innovative American Technology, Inc. | Fabrication of a high performance neutron detector with near zero gamma cross talk |
US8330115B2 (en) * | 2005-12-01 | 2012-12-11 | Innovative American Technology, Inc. | High performance neutron detector with near zero gamma cross talk |
US8046766B2 (en) * | 2007-04-26 | 2011-10-25 | Hewlett-Packard Development Company, L.P. | Process assignment to physical processors using minimum and maximum processor shares |
US8173954B2 (en) | 2008-12-30 | 2012-05-08 | Saint-Gobain Ceramics & Plastics, Inc. | Detector for use in well-logging applications |
US9366732B2 (en) * | 2009-09-04 | 2016-06-14 | Board Of Regents, The University Of Texas System | Estimation of state-of-health in batteries |
US8729902B1 (en) * | 2009-10-05 | 2014-05-20 | White's Electronics, Inc. | Metal detector analysis and display methods |
JP5611357B2 (en) * | 2010-09-02 | 2014-10-22 | 株式会社日立製作所 | Radiation measurement equipment |
US20130021166A1 (en) * | 2011-07-20 | 2013-01-24 | Schlumberger Technology Corporation | System and method for borehole communication |
US9268055B2 (en) | 2011-12-30 | 2016-02-23 | Schlumberger Technology Corporation | Well-logging apparatus including azimuthally spaced radiation detectors |
BR112015013392A2 (en) | 2012-12-12 | 2017-07-11 | Koninklijke Philips Nv | imaging system, and method |
US8927920B2 (en) * | 2013-02-07 | 2015-01-06 | Schlumberger Technology Corporation | Correcting gamma-ray energy spectra for pileup degradation |
US11579072B2 (en) | 2013-03-15 | 2023-02-14 | Particles Plus, Inc. | Personal air quality monitoring system |
US12044611B2 (en) | 2013-03-15 | 2024-07-23 | Particles Plus, Inc. | Particle counter with integrated bootloader |
US9677990B2 (en) * | 2014-04-30 | 2017-06-13 | Particles Plus, Inc. | Particle counter with advanced features |
CN104101894B (en) * | 2013-04-07 | 2017-04-12 | 苏州瑞派宁科技有限公司 | Accumulation event detection and reconstruction method and system |
EP3210044A4 (en) * | 2014-10-23 | 2018-05-23 | Bridgeport Instruments, LLC | Performance stabilization for scintillator-based radiation detectors |
CN104382612A (en) * | 2014-11-13 | 2015-03-04 | 沈阳东软医疗系统有限公司 | CT (Computed Tomography) data recovery method and device |
CN105212954B (en) * | 2015-11-05 | 2018-03-16 | 苏州瑞派宁科技有限公司 | A kind of pulse pile-up event real-time processing method and system |
US10261214B2 (en) * | 2016-05-02 | 2019-04-16 | Schlumberger Technology Corporation | Method and apparatus for separating gamma and neutron signals from a radiation detector and for gain-stabilizing the detector |
US10350644B1 (en) | 2016-11-21 | 2019-07-16 | Mss, Inc. | System and method for induction-based metal detection and high resolution sorting |
AU2018202912B1 (en) * | 2018-04-27 | 2019-06-20 | Southern Innovation International Pty Ltd | Input count rate estimation in radiation pulse detectors |
CN113826030A (en) * | 2019-03-22 | 2021-12-21 | 南方创新国际股份有限公司 | Radiation detection with pulse pile-up nonparametric decomposition |
CN111245409B (en) * | 2019-12-30 | 2021-12-24 | 中国科学院高能物理研究所 | Pulse signal processing method and device |
US11988591B2 (en) | 2020-07-01 | 2024-05-21 | Particles Plus, Inc. | Modular optical particle counter sensor and apparatus |
GB2598137A (en) * | 2020-08-20 | 2022-02-23 | Smiths Detection France S A S | Self-aligning array of detector cells |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1413204A (en) * | 1972-02-17 | 1975-11-12 | Atomic Energy Commission | Pulse processing system |
US5225682A (en) * | 1992-01-24 | 1993-07-06 | The United States Of America As Represented By The United States Department Of Energy | Method and apparatus for providing pulse pile-up correction in charge quantizing radiation detection systems |
US5532944A (en) * | 1994-07-28 | 1996-07-02 | Sorrento Electronics, Inc. | Multi-channel analysis system and method using digital signal processing |
WO2000073819A1 (en) * | 1999-05-27 | 2000-12-07 | Analyser Systems Ag | Method and apparatus for multi-parameter digital pulse pile-up rejection |
US6215122B1 (en) * | 1997-04-04 | 2001-04-10 | Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of National Defence | Method for high speed pulse pile-up rejection |
Family Cites Families (56)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4658216A (en) | 1983-07-14 | 1987-04-14 | The United States Of America As Represented By The Department Of Energy | High resolution, high rate X-ray spectrometer |
US5210423A (en) | 1989-09-22 | 1993-05-11 | Siemens Gammasonics, Inc. | Method and apparatus for unpiling pulses generated by piled-up scintillation events |
FR2686982B1 (en) | 1992-01-30 | 1997-12-19 | Schlumberger Services Petrol | METHOD FOR SPECTROSCOPY OF NUCLEAR PULSES WITH A HIGH COUNTING RATE. |
US5349193A (en) | 1993-05-20 | 1994-09-20 | Princeton Gamma Tech, Inc. | Highly sensitive nuclear spectrometer apparatus and method |
EP0880821B1 (en) | 1995-08-14 | 2003-12-17 | William K. Warburton | Method and apparatus for digitally based high speed x-ray spectrometer |
US5884234A (en) * | 1996-08-30 | 1999-03-16 | Schlumberger Technology Corporation | Method for pulse shape regulation and discrimination in a nuclear spectroscopy system |
FR2753279B1 (en) | 1996-09-09 | 1998-11-13 | Smv Int | PROCESS FOR TREATING PULSES DELIVERED BY A GAMMA CAMERA AND GAMMA CAMERA IMPLEMENTING THIS METHOD |
JPH10186041A (en) | 1996-12-21 | 1998-07-14 | Horiba Ltd | Pulse processor of x-ray analyzing device |
US6936822B2 (en) | 1997-05-07 | 2005-08-30 | Board Of Regents, The University Of Texas System | Method and apparatus to prevent signal pile-up |
WO1999058999A1 (en) | 1998-05-09 | 1999-11-18 | Packard Bioscience Company | Digital pulse de-randomization for radiation spectroscopy |
US6347585B1 (en) | 1998-08-04 | 2002-02-19 | Goss Graphic Systems, Inc. | Variable gap stabilizer |
US6207957B1 (en) * | 1998-09-18 | 2001-03-27 | The Regents Of The University Of California | System to quantify gamma-ray radial energy deposition in semiconductor detectors |
US6160259A (en) | 1998-10-23 | 2000-12-12 | Adac Laboratories | Channel-specific control of pulse integration in a gamma camera system in response to pulse pile-up |
US6366872B1 (en) | 1998-12-23 | 2002-04-02 | Varian, Inc. | Data prediction in DSP methods for transient spectroscopies |
JP3836629B2 (en) | 1999-06-22 | 2006-10-25 | 日本電子株式会社 | Digital pulse processor |
EP1221038A4 (en) * | 1999-07-21 | 2004-09-08 | Applera Corp | Luminescence detection workstation |
US7065473B2 (en) * | 1999-08-27 | 2006-06-20 | William K. Warburton | Method and apparatus for improving resolution in spectrometers processing output steps from non-ideal signal sources |
JP3907910B2 (en) | 2000-03-28 | 2007-04-18 | 日本電子株式会社 | Digital pulse processor used in energy dispersive spectrometer |
US20040195512A1 (en) * | 2000-05-16 | 2004-10-07 | Crosetto Dario B. | Method and apparatus for anatomical and functional medical imaging |
JP2001337168A (en) | 2000-05-29 | 2001-12-07 | Jeol Ltd | Ad converter for measuring radiation |
US7133568B2 (en) * | 2000-08-04 | 2006-11-07 | Nikitin Alexei V | Method and apparatus for analysis of variables |
JP4160275B2 (en) * | 2001-05-28 | 2008-10-01 | 浜松ホトニクス株式会社 | Energy measuring method and measuring device |
MXPA04002498A (en) * | 2001-09-19 | 2005-04-11 | Joule Microsystems Canada Inc | A spectrometer incorporating signal matched filtering. |
US20030076914A1 (en) | 2001-10-23 | 2003-04-24 | Tiller Donald E. | Method and apparatus for measuring radiation in a borehole |
US6590957B1 (en) | 2002-03-13 | 2003-07-08 | William K. Warburton | Method and apparatus for producing spectra corrected for deadtime losses in spectroscopy systems operating under variable input rate conditions |
JP3980451B2 (en) | 2002-08-30 | 2007-09-26 | 株式会社東芝 | Waveform discrimination device |
FR2870603B1 (en) | 2004-05-19 | 2006-07-28 | Commissariat Energie Atomique | MEASUREMENT AND PROCESSING OF A SIGNAL COMPRISING ELEMENTARY PULSE STACKS |
US7966155B2 (en) | 2004-06-04 | 2011-06-21 | William K. Warburton | Method and apparatus for improving detection limits in x-ray and nuclear spectroscopy systems |
GB2421317B (en) | 2004-12-15 | 2009-02-11 | Agilent Technologies Inc | A method and apparatus for detecting leading pulse edges |
US7439515B2 (en) | 2005-06-17 | 2008-10-21 | Siemens Medical Solutions Usa, Inc. | Pipeline processing of pulse pile-up correction in a nuclear medicine imaging system |
WO2007049168A2 (en) | 2005-10-28 | 2007-05-03 | Koninklijke Philips Electronics, N.V. | Method and apparatus for spectral computed tomography |
JP4358814B2 (en) | 2005-11-09 | 2009-11-04 | 花王株式会社 | Sample analysis method |
JP4706566B2 (en) | 2006-06-09 | 2011-06-22 | 株式会社島津製作所 | X-ray analysis signal processor |
WO2007146350A2 (en) | 2006-06-14 | 2007-12-21 | Baker Hughes Incorporated | Pileup rejection |
JP4787989B2 (en) | 2006-11-10 | 2011-10-05 | 横河電機株式会社 | Peak detection circuit, multi-channel analyzer and radiation measurement system |
US7725281B2 (en) | 2006-12-21 | 2010-05-25 | Canberra Industries, Inc. | Method for pole-zero adjustment in radiation measurement systems |
US7430481B2 (en) | 2007-01-12 | 2008-09-30 | R. J. Lee Group, Inc. | Methods for detecting and analyzing piled-up X-rays in an X-ray spectrometry system |
CN101680956B (en) | 2007-06-19 | 2013-02-13 | 皇家飞利浦电子股份有限公司 | Digital pulse processing for multi-spectral photon counting readout circuits |
US7741609B2 (en) | 2007-08-03 | 2010-06-22 | Pulsetor, Llc | Adapting a high-performance pulse processor to an existing spectrometry system |
WO2009032452A1 (en) | 2007-08-03 | 2009-03-12 | Pulsetor, Llc | Pileup rejection in an energy-dispersive radiation spectrometry system |
US8415635B2 (en) | 2007-10-18 | 2013-04-09 | Koninklijke Philips Electronics N.V. | Particle-counting apparatus with pulse shortening |
US8003948B2 (en) | 2007-11-02 | 2011-08-23 | University Of Washington | Data acquisition for positron emission tomography |
US9008969B2 (en) | 2007-12-10 | 2015-04-14 | Schlumberger Technology Corporation | Pulsed x-ray signal processing |
JP5043728B2 (en) | 2008-03-19 | 2012-10-10 | 日本電子株式会社 | Pulse processor for radiation measurement |
US7999220B2 (en) | 2008-05-30 | 2011-08-16 | Precision Energy Services, Inc. | Borehole measurements using a fast and high energy resolution gamma ray detector assembly |
US7807973B2 (en) | 2008-08-01 | 2010-10-05 | Pulsetor, Llc | Pileup rejection in an energy-dispersive radiation spectrometry system |
CN101666227B (en) | 2008-09-03 | 2013-02-20 | 中国石油天然气集团公司 | Method for acquiring natural gamma spectra |
US8227761B2 (en) | 2008-09-19 | 2012-07-24 | Canberra Industries, Inc. | True coincidence summing correction and total efficiency computation for radionuclide spectroscopy analysis |
DE102008048306B4 (en) | 2008-09-22 | 2016-06-09 | Siemens Healthcare Gmbh | Method for the detection of X-ray radiation and X-ray system |
DE102009018995B4 (en) | 2009-04-27 | 2017-10-12 | Siemens Healthcare Gmbh | Method and circuit for determining an intensity of ionizing radiation |
WO2011002452A1 (en) | 2009-06-30 | 2011-01-06 | Analogic Corporation | Enhanced photon detection for scanner |
EP2290403A1 (en) | 2009-08-28 | 2011-03-02 | Paul Scherrer Institut | X-ray detector with integrating readout chip for single photon resolution |
FR2951036A1 (en) | 2009-10-01 | 2011-04-08 | Commissariat Energie Atomique | DEVICE FOR PROCESSING A SIGNAL DELIVERED BY A RADIATION DETECTOR |
KR101092648B1 (en) | 2010-04-28 | 2011-12-13 | 서강대학교산학협력단 | Method for detecting arrival time of gamma ray pulse signal in PET detector |
JP5611357B2 (en) | 2010-09-02 | 2014-10-22 | 株式会社日立製作所 | Radiation measurement equipment |
CN102073059B (en) | 2010-12-31 | 2013-05-22 | 华中科技大学 | Digital pileup waveform processing method and system |
-
2005
- 2005-09-16 EP EP05779842A patent/EP1794884A4/en not_active Ceased
- 2005-09-16 CA CA2621736A patent/CA2621736C/en active Active
- 2005-09-16 WO PCT/AU2005/001423 patent/WO2006029475A1/en active Application Filing
- 2005-09-16 CN CN2005800392488A patent/CN101124725B/en active Active
- 2005-09-16 JP JP2007531539A patent/JP5022902B2/en active Active
-
2007
- 2007-03-13 US US11/685,719 patent/US7383142B2/en active Active
-
2008
- 2008-05-30 US US12/130,671 patent/US20090048795A1/en not_active Abandoned
-
2011
- 2011-10-27 US US13/283,410 patent/US10310099B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1413204A (en) * | 1972-02-17 | 1975-11-12 | Atomic Energy Commission | Pulse processing system |
US5225682A (en) * | 1992-01-24 | 1993-07-06 | The United States Of America As Represented By The United States Department Of Energy | Method and apparatus for providing pulse pile-up correction in charge quantizing radiation detection systems |
US5532944A (en) * | 1994-07-28 | 1996-07-02 | Sorrento Electronics, Inc. | Multi-channel analysis system and method using digital signal processing |
US6215122B1 (en) * | 1997-04-04 | 2001-04-10 | Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of National Defence | Method for high speed pulse pile-up rejection |
WO2000073819A1 (en) * | 1999-05-27 | 2000-12-07 | Analyser Systems Ag | Method and apparatus for multi-parameter digital pulse pile-up rejection |
Non-Patent Citations (2)
Title |
---|
R.J. KOMAR; H.-B. MAK: "Digital signal processing for BGO detectors", NUCL. LNST. AND MET . A, vol. 336, 1993, pages 246 - 252 |
See also references of EP1794884A4 * |
Cited By (43)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007327902A (en) * | 2006-06-09 | 2007-12-20 | Shimadzu Corp | Signal processing device for x ray analysis |
US8581200B2 (en) | 2006-11-17 | 2013-11-12 | Koninklijke Philips N.V. | Radiation detector with multiple electrodes on a sensitive layer |
JP2011503550A (en) * | 2007-11-02 | 2011-01-27 | ユニヴァーシティ オブ ワシントン | Data acquisition for positron emission tomography |
AU2009230874B2 (en) * | 2008-03-31 | 2014-06-26 | Southern Innovation International Pty Ltd | Screening method and apparatus |
WO2009121132A1 (en) * | 2008-03-31 | 2009-10-08 | Southern Innovation International Pty Ltd | Radiation imaging method with individual signal resolution |
AU2009230876B2 (en) * | 2008-03-31 | 2014-07-10 | Southern Innovation International Pty Ltd | Radiation imaging method with individual signal resolution |
US9310513B2 (en) | 2008-03-31 | 2016-04-12 | Southern Innovation International Pty Ltd. | Method and apparatus for borehole logging |
EP2260296A4 (en) * | 2008-03-31 | 2016-04-06 | Southern Innovation Internat Pty Ltd | Screening method and apparatus |
EP2260295A4 (en) * | 2008-03-31 | 2015-12-23 | Southern Innovation Internat Pty Ltd | Radiation imaging method with individual signal resolution |
AU2009230875B2 (en) * | 2008-03-31 | 2014-07-10 | Southern Innovation International Pty Ltd | Method and apparatus for borehole logging |
JP2011516838A (en) * | 2008-03-31 | 2011-05-26 | サザン イノヴェーション インターナショナル プロプライアトリー リミテッド | Screening method and apparatus |
JP2011516839A (en) * | 2008-03-31 | 2011-05-26 | サザン イノヴェーション インターナショナル プロプライアトリー リミテッド | Method and apparatus for borehole logging |
JP2011519415A (en) * | 2008-03-31 | 2011-07-07 | サザン イノヴェーション インターナショナル プロプライアトリー リミテッド | Radiation imaging using individual signal resolution |
EP2260334A4 (en) * | 2008-03-31 | 2015-03-11 | Southern Innovation Internat Pty Ltd | Method and apparatus for borehole logging |
US8954300B2 (en) | 2008-03-31 | 2015-02-10 | Southern Innovation International Pty Ltd. | Screening method and apparatus |
JP2014240834A (en) * | 2008-03-31 | 2014-12-25 | サザン イノヴェーション インターナショナル プロプライアトリー リミテッド | Method and apparatus for borehole logging |
US8471210B2 (en) | 2008-03-31 | 2013-06-25 | Southern Innovation International Pty Ltd. | Radiation imaging method with individual signal resolution |
WO2009121131A1 (en) * | 2008-03-31 | 2009-10-08 | Southern Innovation International Pty Ltd | Method and apparatus for borehole logging |
JP2014041147A (en) * | 2008-03-31 | 2014-03-06 | Southern Innovation International Pty Ltd | Screening method and device |
WO2009121130A1 (en) * | 2008-03-31 | 2009-10-08 | Southern Innovation International Pty Ltd | Screening method and apparatus |
US7999220B2 (en) | 2008-05-30 | 2011-08-16 | Precision Energy Services, Inc. | Borehole measurements using a fast and high energy resolution gamma ray detector assembly |
GB2460345A (en) * | 2008-05-30 | 2009-12-02 | Precision Energy Services Inc | Borehole logging using LaBr3 scintillation detector and digital spectrometer which rejects pile-up pulses |
GB2460345B (en) * | 2008-05-30 | 2010-06-09 | Precision Energy Services Inc | Borehole measurements using a fast and high energy resolution gamma ray detector assembly |
AU2009202086B2 (en) * | 2008-05-30 | 2010-09-09 | Precision Energy Services, Inc. | Borehole measurements using a fast and high energy resolution gamma ray detector assembly |
USRE45226E1 (en) | 2008-05-30 | 2014-11-04 | Precision Energy Services, Inc. | Borehole measurements using a fast and high energy resolution gamma ray detector assembly |
AU2009202086B8 (en) * | 2008-05-30 | 2010-09-30 | Precision Energy Services, Inc. | Borehole measurements using a fast and high energy resolution gamma ray detector assembly |
CN102549566A (en) * | 2008-12-18 | 2012-07-04 | 南方创新国际股份有限公司 | Method and apparatus for resolving piled-up pulses by using a mathematical transform |
US8812268B2 (en) | 2008-12-18 | 2014-08-19 | Southern Innovation International Pty. Ltd. | Method and apparatus for resolving piled-up pulses by using a mathematical transform |
EP2359261A4 (en) * | 2008-12-18 | 2015-08-19 | Southern Innovation Internat Pty Ltd | Method and apparatus for resolving piled-up pulses by using a mathematical transform |
WO2010068996A1 (en) | 2008-12-18 | 2010-06-24 | Southern Innovation International Pty Ltd | Method and apparatus for resolving piled-up pulses by using a mathematical transform |
JPWO2011039819A1 (en) * | 2009-10-01 | 2013-02-21 | 株式会社島津製作所 | Pulse height analyzer and nuclear medicine diagnostic apparatus provided with the same |
JP5531021B2 (en) * | 2009-10-01 | 2014-06-25 | 株式会社島津製作所 | Pulse height analyzer and nuclear medicine diagnostic apparatus provided with the same |
US9885674B2 (en) | 2009-12-09 | 2018-02-06 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Method and device for recognition of a material making use of its transmission function |
CN106170919B (en) * | 2013-12-11 | 2019-05-28 | 南方创新国际股份有限公司 | Method and apparatus for parsing the signal in data |
CN106170919A (en) * | 2013-12-11 | 2016-11-30 | 南方创新国际股份有限公司 | For resolving the method and apparatus of the signal in data |
WO2015085372A1 (en) * | 2013-12-11 | 2015-06-18 | Southern Innovation International Pty Ltd | Method and apparatus for resolving signals in data |
US11079502B2 (en) | 2013-12-11 | 2021-08-03 | Southern Innovation International Pty Ltd | Method and apparatus for resolving signals in data |
WO2016082006A1 (en) | 2014-11-30 | 2016-06-02 | Southern Innovation International Pty Ltd | Device and method for material characterisation |
US10485994B2 (en) | 2015-03-30 | 2019-11-26 | Sumitomo Heavy Industries, Ltd. | Charged particle beam treatment apparatus |
JP2017086901A (en) * | 2015-11-02 | 2017-05-25 | 東芝メディカルシステムズ株式会社 | Data collection device, x-ray ct apparatus, and nuclear medicine diagnostic apparatus |
WO2017205914A1 (en) | 2016-05-30 | 2017-12-07 | Southern Innovation International Pty Ltd | Material characterisation system and method |
WO2017205912A1 (en) | 2016-05-30 | 2017-12-07 | Southern Innovation International Pty Ltd | Material characterisation system and method |
EP3546926A1 (en) * | 2018-03-28 | 2019-10-02 | PicoQuant Innovations GmbH | Method and device for determining photon arrival times by means of machine learning |
Also Published As
Publication number | Publication date |
---|---|
CN101124725A (en) | 2008-02-13 |
US20070147702A1 (en) | 2007-06-28 |
CA2621736C (en) | 2018-07-31 |
US7383142B2 (en) | 2008-06-03 |
EP1794884A1 (en) | 2007-06-13 |
CN101124725B (en) | 2012-06-20 |
EP1794884A4 (en) | 2009-03-04 |
US10310099B2 (en) | 2019-06-04 |
JP5022902B2 (en) | 2012-09-12 |
US20120041700A1 (en) | 2012-02-16 |
US20090048795A1 (en) | 2009-02-19 |
CA2621736A1 (en) | 2006-03-23 |
JP2008513740A (en) | 2008-05-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7383142B2 (en) | Method and apparatus for resolving individual signals in detector output data | |
AU2009230875B2 (en) | Method and apparatus for borehole logging | |
AU2009230876B2 (en) | Radiation imaging method with individual signal resolution | |
AU2016200724B2 (en) | Method and apparatus for resolving piled-up pulses by using a mathematical transform | |
RU2681377C1 (en) | Method and apparatus for resolving signals in data | |
AU2009230874B2 (en) | Screening method and apparatus | |
JP5111105B2 (en) | Measurement method and analysis apparatus including signal processing including a chain of main impulses | |
AU2005284602B2 (en) | Method and apparatus for resolving individual signals in detector output data. | |
Coulon et al. | Recent Developments in Count Rate Processing Associated with Radiation Monitoring Systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AK | Designated states |
Kind code of ref document: A1 Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KM KP KR KZ LC LK LR LS LT LU LV LY MA MD MG MK MN MW MX MZ NA NG NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SM SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW |
|
AL | Designated countries for regional patents |
Kind code of ref document: A1 Designated state(s): GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LT LU LV MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
WWE | Wipo information: entry into national phase |
Ref document number: 2005284602 Country of ref document: AU |
|
ENP | Entry into the national phase |
Ref document number: 2005284602 Country of ref document: AU Date of ref document: 20050916 Kind code of ref document: A |
|
WWP | Wipo information: published in national office |
Ref document number: 2005284602 Country of ref document: AU |
|
WWE | Wipo information: entry into national phase |
Ref document number: 11685719 Country of ref document: US |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2007531539 Country of ref document: JP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2005779842 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 200580039248.8 Country of ref document: CN |
|
WWP | Wipo information: published in national office |
Ref document number: 2005779842 Country of ref document: EP |
|
WWP | Wipo information: published in national office |
Ref document number: 11685719 Country of ref document: US |
|
WWG | Wipo information: grant in national office |
Ref document number: 2005284602 Country of ref document: AU |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2621736 Country of ref document: CA |