CN-116203613-B - Processing method and device of scintillation pulse, digitizing equipment and storage medium

CN116203613BCN 116203613 BCN116203613 BCN 116203613BCN-116203613-B

Abstract

The application discloses a method and a device for processing scintillation pulses, digital equipment and a storage medium. The processing method comprises the steps of obtaining a function model corresponding to a scintillation pulse, wherein the function model is used for representing a pulse waveform of the scintillation pulse, presetting a plurality of thresholds, carrying out multi-threshold sampling on the scintillation pulse based on the thresholds to obtain sampling data, determining a corresponding relation table between pulse energy related to the thresholds and threshold duration, determining whether the scintillation pulse is a stacked pulse based on the sampling data and the corresponding relation table, and if so, determining a composition energy value of a composition pulse forming the scintillation pulse based on the sampling data and the function model. The application can restore the stacked pulse, and accurately calculate the energy and count the pulse after the restoration, thereby avoiding the situation of partial pulse miscounting caused by pulse stacking.

Inventors

FANG LEI
ZHANG BO
CHEN WEICAO
YANG LINGLI
HU YUN
HUANG WENLUE

Assignees

合肥锐世数字科技有限公司

Dates

Publication Date: 20260512
Application Date: 20221230

Claims (20)

1. A method of processing a scintillation pulse, the method comprising: Obtaining a function model corresponding to the scintillation pulse, wherein the function model is used for representing the pulse waveform of the scintillation pulse; Presetting a plurality of thresholds, and carrying out multi-threshold sampling on the scintillation pulse based on the thresholds so as to acquire sampling data; determining a correspondence table between pulse energy associated with the threshold and threshold duration; determining whether the scintillation pulse is a stacked pulse based on the sampling data and the corresponding relation table, wherein the determining comprises determining a comparison threshold duration based on the sampling data, determining whether a first ratio between the comparison threshold duration and a threshold duration in the corresponding relation table is greater than a preset ratio based on the comparison threshold duration and the corresponding relation table; if the scintillation pulse is a stacked pulse, determining a composition energy value of a composition pulse constituting the scintillation pulse based on the sampling data and the function model; If the scintillation pulse is a non-stacked pulse, determining a target energy value of the scintillation pulse based on the sampling data and/or the corresponding relation table, and designating the target energy value to participate in pulse counting and energy spectrum drawing.
2. The processing method according to claim 1, wherein the determining the correspondence table includes: Acquiring a plurality of known pulses with known pulse energy, wherein the known pulses are unstacked pulses; Multi-threshold sampling the known pulse based on the threshold value, determining a first time and a second time at which the known pulse crosses the threshold value; determining the threshold duration based on the first time and the second time; the correspondence table is determined based on a plurality of pulse energies and a plurality of threshold durations.
3. The processing method according to claim 2, wherein the correspondence table reflects correspondence between the number of threshold values crossed by the known pulse, pulse energy of the known pulse, and a first average threshold duration including an average value of threshold durations corresponding to the respective threshold values.
4. The processing method according to claim 2, wherein the correspondence table reflects a correspondence between a threshold duration of a lowest threshold value crossed by the known pulse and a pulse energy of the known pulse.
5. The processing method according to claim 2, wherein the correspondence table reflects a correspondence between a threshold duration of a highest threshold value crossed by the known pulse and a pulse energy of the known pulse.
6. The processing method according to claim 2, wherein the correspondence table reflects correspondence between pulse energies of the known pulses and second average threshold durations including an average value of threshold durations corresponding to a preset one of threshold values traversed by the known pulses.
7. The processing method according to claim 1, wherein the determining of the composition energy values of the composition pulses constituting the scintillation pulse based on the sampling data and the function model includes: determining pulse waveforms of a first constituent pulse and a second constituent pulse constituting the scintillation pulse based on the sampling data and the function model; based on pulse waveforms of the first constituent pulse and the second constituent pulse, constituent energy values of the first constituent pulse and the second constituent pulse are determined.
8. The processing method of claim 7, wherein the sampled data includes a plurality of first threshold-time pairs when rising edges of the scintillation pulses cross a plurality of thresholds and a plurality of second threshold-time pairs when falling edges of the scintillation pulses cross a plurality of thresholds, wherein a fitting operation is performed based on the plurality of first threshold-time pairs to determine a first expression of the functional model to characterize a pulse shape of a first constituent pulse, and wherein a fitting operation is performed based on the plurality of second threshold-time pairs to determine a second expression of the functional model to characterize a pulse shape of a second constituent pulse.
9. The processing method according to claim 8, characterized in that the processing method further comprises: Integrating a first curve corresponding to the first expression to determine a first composition energy value of the first composition pulse; and integrating a second curve corresponding to the second expression to determine a second composition energy value of the second composition pulse.
10. The processing method according to claim 9, characterized in that the processing method further comprises: designating the first and second constituent energy values to participate in pulse counting and energy spectrum drawing.
11. The method of processing of claim 1, wherein determining the target energy value of the scintillation pulse comprises: Performing a fitting operation based on the sampled data to determine a target expression of the functional model; and integrating a target curve corresponding to the target expression to determine the target energy value.
12. The method of processing of claim 1, wherein determining the target energy value of the scintillation pulse comprises: Determining a target threshold duration based on the sampled data; And comparing the target threshold duration with each related threshold duration in the corresponding relation table to determine a second ratio, and when the second ratio is smaller than the threshold duration of the preset ratio, designating pulse energy corresponding to the threshold duration as the target energy value.
13. The method of processing of claim 1, wherein determining the target energy value of the scintillation pulse comprises: Determining a target threshold duration based on the sampled data; Determining a relationship function between pulse energy and threshold duration based on the correspondence table; the target energy value is determined based on the target threshold duration and the relationship function.
14. A processing device for scintillation pulses, the processing device comprising: The acquisition module is configured to acquire a function model corresponding to the scintillation pulse, wherein the function model is used for representing the pulse waveform of the scintillation pulse; The sampling module is configured to preset a plurality of thresholds and perform multi-threshold sampling on the scintillation pulse based on the thresholds so as to acquire sampling data; A determining module configured to determine a correspondence table between pulse energy associated with the threshold and a threshold duration; The judgment module is configured to determine whether the scintillation pulse is a stacked pulse or not based on the sampling data and the corresponding relation table, and comprises a comparison threshold duration time based on the sampling data, a first ratio value between the comparison threshold duration time and the threshold duration time in the corresponding relation table is determined to be larger than a preset ratio value based on the comparison threshold duration time and the corresponding relation table, and if the first ratio value is larger than the preset ratio value, the scintillation pulse is determined to be the stacked pulse; And the calculation module is configured to determine the component energy value of the component pulse forming the scintillation pulse based on the sampling data and the function model when the pulse is a stacked pulse, and determine the target energy value of the scintillation pulse based on the sampling data and/or the corresponding relation table when the scintillation pulse is a non-stacked pulse, and participate in pulse counting and energy spectrum drawing by utilizing the target energy value.
15. The processing apparatus according to claim 14, wherein to determine the correspondence table, the determining module is configured to: Acquiring a plurality of known pulses with known pulse energy, wherein the known pulses are unstacked pulses; Multi-threshold sampling the known pulse based on the threshold value, determining a first time and a second time at which the known pulse crosses the threshold value; determining the threshold duration based on the first time and the second time; the correspondence table is determined based on a plurality of pulse energies and a plurality of threshold durations.
16. The processing apparatus of claim 15, wherein the correspondence table reflects a correspondence between a number of thresholds traversed by the known pulse, a pulse energy of the known pulse, and a first average threshold duration, the first average threshold duration comprising an average of threshold durations corresponding to the respective thresholds.
17. The processing apparatus of claim 15, wherein the correspondence table reflects a correspondence between a threshold duration of a lowest threshold traversed by the known pulse and a pulse energy of the known pulse.
18. The processing apparatus according to claim 15, wherein the correspondence table reflects a correspondence between a threshold duration of a highest threshold value crossed by the known pulse and a pulse energy of the known pulse.
19. The processing apparatus of claim 15, wherein the correspondence table reflects a correspondence between pulse energies of the known pulses and a second average threshold duration comprising an average of threshold durations corresponding to a preset one of the threshold values traversed by the known pulses.
20. The processing apparatus of claim 14, wherein to determine a composition energy value of a composition pulse constituting the scintillation pulse based on the sampling data and the functional model, the calculation module is configured to: determining pulse waveforms of a first constituent pulse and a second constituent pulse constituting the scintillation pulse based on the sampling data and the function model; based on pulse waveforms of the first constituent pulse and the second constituent pulse, constituent energy values of the first constituent pulse and the second constituent pulse are determined.

Description

Processing method and device of scintillation pulse, digitizing equipment and storage medium Technical Field The present application relates to the field of data processing, and in particular, to a method and apparatus for processing scintillation pulse, a digitizing device, and a storage medium. Background Photon counting CT when performing high energy particle/ray detection, the photon counting CT detector can directly convert the detected high energy particles/rays, such as X-rays, into electrical signals. By calculating the energy of each X-ray photon, a count of photons of different energies can be obtained. The higher the X-ray energy, the greater the pulse amplitude. The lower the X-ray energy, the smaller the pulse amplitude. Photon counting CT detectors do photon counting in such a way that the energy of the X-rays is not directly obtained, but the number of photons of each energy segment is obtained. The energy segment count may be obtained by testing the relationship between the pulse amplitude and the pulse energy with a high-speed oscilloscope. The pulse amplitude T corresponding to the energy value E is taken as a comparison threshold, the pulse exceeding the threshold is regarded as energy larger than E, and the pulse smaller than the threshold is regarded as energy smaller than E. In photon counting CT, pulse peaks T0, T1, T2, T3 corresponding to four energies, 25kev,50kev,75kev,100kev, are typically calculated first. Counting the number of photons of T0-T1, T1-T2, T2-T3 and more than T3. Photon counting of the segment energy is achieved. In the prior art, a method of obtaining the segmented energy photon count by presetting a limited threshold value can be utilized, but the accurate energy of each photon can not be obtained, and the real photon count can not be realized. Meanwhile, if one pulse gets a critical position of a peak value and a threshold value, misjudgment on a photon energy segment may occur. In addition, the photon count CT detector detects photons with the number of 10 cps/mm <2 >. The stacking of the electrical pulses of adjacent photon outputs occurs with a high probability, i.e. pileup. For the pileup problem, one approach is to drop both 2 photons that occur in pileup, without recovering the 2 pulse distinction. Because the pulse of 2 photons synthesizes one after pileup occurs, the combined pulse energy will be larger than the conventional energy value, and can be directly discarded without participating in counting. Clearly, in so doing, the photon counts at different energies will deviate. Disclosure of Invention The technical problem to be solved by the embodiments of the present application is how to achieve accurate calculation of pulse energy generated by a single photon. In order to solve the problems, the application discloses a method and a device for processing scintillation pulse, a digitizing device and a storage medium. According to a first aspect of the present application, a method of processing scintillation pulses is provided. The method comprises the steps of obtaining a function model corresponding to a scintillation pulse, wherein the function model is used for representing a pulse waveform of the scintillation pulse, presetting a plurality of thresholds, carrying out multi-threshold sampling on the scintillation pulse based on the thresholds to obtain sampling data, determining a corresponding relation table between pulse energy related to the thresholds and threshold duration, determining whether the scintillation pulse is a stacked pulse based on the sampling data and the corresponding relation table, and if so, determining a composition energy value of a composition pulse forming the scintillation pulse based on the sampling data and the function model. According to some embodiments of the application, the determining the correspondence table includes obtaining a plurality of known pulses having known pulse energies, the known pulses being unstacked pulses, multi-threshold sampling the known pulses based on the threshold, determining a first time and a second time at which the known pulses cross the threshold, determining the threshold duration based on the first time and the second time, and determining the correspondence table based on the plurality of pulse energies and the plurality of threshold durations. According to some embodiments of the application, the correspondence table reflects a correspondence between a number of threshold values traversed by the known pulse, a pulse energy of the known pulse, and a first average threshold duration including an average of threshold durations corresponding to the respective threshold values. According to some embodiments of the application, the correspondence table reflects a correspondence between a threshold duration of a lowest threshold value traversed by the known pulse and a pulse energy of the known pulse. According to some embodiments of the application, the correspondence table reflects a