CN111245409B - Pulse signal processing method and device - Google Patents
Pulse signal processing method and device Download PDFInfo
- Publication number
- CN111245409B CN111245409B CN201911394385.3A CN201911394385A CN111245409B CN 111245409 B CN111245409 B CN 111245409B CN 201911394385 A CN201911394385 A CN 201911394385A CN 111245409 B CN111245409 B CN 111245409B
- Authority
- CN
- China
- Prior art keywords
- pulse signal
- waveform
- section
- pulse
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03K—PULSE TECHNIQUE
- H03K5/00—Manipulating of pulses not covered by one of the other main groups of this subclass
- H03K5/13—Arrangements having a single output and transforming input signals into pulses delivered at desired time intervals
- H03K5/135—Arrangements having a single output and transforming input signals into pulses delivered at desired time intervals by the use of time reference signals, e.g. clock signals
Landscapes
- Physics & Mathematics (AREA)
- Nonlinear Science (AREA)
- Measurement Of Radiation (AREA)
Abstract
The application discloses a pulse signal processing method and a device, wherein the method comprises the following steps: carrying out amplitude normalization processing on the pulse signal according to the ideal waveform model; dividing the normalized pulse signal into at least two sections, wherein the first section of pulse signal is a pulse signal generated from the pulse signal to a peak value, and the second section of pulse signal is a pulse signal remained in the pulse signal except the first section of pulse signal; calculating a first variance difference value corresponding to the first section of pulse signals and a second variance difference value corresponding to the second section of pulse signals; judging whether the first variance difference value is smaller than a first threshold value; if yes, judging whether the waveform of the second section of pulse signal needs to be subjected to waveform recovery or not according to the second variance difference value and the second threshold value. The invention reserves and recovers the signals which are not seriously piled up. Thereby enabling the device to detect imaging in a high dose environment at a certain range.
Description
Technical Field
The present invention relates generally to the field of digital signal processing, and more particularly to a method and apparatus for processing a pulse signal.
Background
In the signal processing process of nuclear imaging and nuclear detection, the general method for signal accumulation is as follows: the comparison method searches signal peaks, judges whether the distance between the two peaks exceeds the length of a wavelength, if the distance is greater than the length of the wavelength, the former signal is considered to be a complete pulse signal and can be reserved, otherwise, the former signal is considered to be accumulated, and the two signals are discarded. However, in a high dose rate environment, a large amount of signals are accumulated, the general method is not applicable any more, and a large amount of accumulated waveforms are discarded, so that the equipment cannot respond, and the waveform waste is caused while the imaging cannot be detected.
Disclosure of Invention
In view of the above-mentioned drawbacks and deficiencies in the prior art, it is desirable to provide a pulse signal processing method and apparatus, which can retain and recover signals that are not heavily accumulated, so that the device can detect imaging in a high dose rate environment at a certain range.
In a first aspect, the present invention provides a pulse signal processing method, including:
carrying out amplitude normalization processing on the pulse signal according to the ideal waveform model to obtain a normalized pulse signal;
dividing the normalized pulse signal into at least two sections, wherein the first section of pulse signal is a pulse signal generated from the pulse signal to a peak value, and the second section of pulse signal is a pulse signal remained in the pulse signal except the first section of pulse signal;
calculating a first variance difference value corresponding to the first section of pulse signals and a second variance difference value corresponding to the second section of pulse signals;
judging whether the first variance difference value is smaller than a first threshold value;
if yes, judging whether the waveform of the second section of pulse signal needs to be subjected to waveform recovery or not according to the second variance difference value and the second threshold value.
Further, when the normalized pulse signal is divided into two sections,
judging whether the waveform of the second section of pulse signal needs to be subjected to waveform recovery according to the second variance difference value and the second variance difference threshold value, wherein the judging step comprises the following steps:
if the second variance difference value is smaller than the second threshold value, the pulse signal is an ideal waveform and does not need to be recovered;
and if the second variance difference value is not smaller than the second threshold value, a pile-up signal exists in the second section of pulse signal, and the waveform of the pile-up signal is recovered.
Preferably, when the normalized pulse signal is divided into three or more segments, the second segment of pulse signal includes at least two segments, and further includes:
sequentially judging whether the square difference value corresponding to each section in the second section of pulse signal is smaller than the corresponding square difference threshold value;
if so, the pulse signal is an ideal waveform and does not need to be subjected to waveform recovery;
and if the difference of squares corresponding to the last section in the second section of pulse signals is not smaller than the corresponding threshold value, accumulating signals exist in the last section of the second section of pulse signals, and waveform recovery is carried out on the accumulating signals.
Further, the waveform recovery method includes:
dividing the ideal waveform amplitude corresponding to the accumulation signal by a normalization coefficient to obtain a recovered waveform, wherein the normalization coefficient is a normalization coefficient for normalizing the pulse signal according to an ideal waveform model;
the waveform of the pile-up signal is replaced with a restored waveform.
Preferably, before performing amplitude normalization on the pulse signal according to the ideal waveform model of the pulse signal to obtain a normalized pulse signal, the method further includes:
acquiring a pulse signal;
judging the energy of the pulse signal;
an ideal waveform model matching the pulse signal energy magnitude is selected.
In a second aspect, the present invention provides a pulse signal processing apparatus, including:
the processing unit is used for carrying out amplitude normalization processing on the pulse signals by the ideal waveform model to obtain normalized pulse signals;
the segmentation unit is used for dividing the normalized pulse signal into at least two segments, wherein the first segment of pulse signal is a pulse signal generated from the pulse signal to a peak value, and the second segment of pulse signal is a pulse signal remained in the pulse signal except the first segment of pulse signal;
the calculating unit is used for calculating a first variance difference value corresponding to the first section of pulse signals and a second variance difference value corresponding to the second section of pulse signals;
the first judgment unit is used for judging whether the first variance difference value is smaller than a first threshold value or not, and if so, the second judgment unit is activated;
and the second judging unit judges whether the waveform of the pulse signal needs to be subjected to waveform recovery or not according to the second variance difference value and the second threshold value.
Further, when the normalized pulse signal is divided into three or more segments, the second segment of pulse signal includes at least two segments, and further includes:
at least two judging subunits, which are used for sequentially judging whether the square difference value corresponding to each section in the second section of pulse signal is smaller than the corresponding threshold value; if so, the pulse signal is an ideal waveform and does not need to be subjected to waveform recovery;
and if the difference of squares corresponding to the last section of the second section of pulse signals is not smaller than the corresponding threshold value, accumulating signals exist in the last section of the second section of pulse signals, and waveform recovery is carried out on the accumulating signals.
Further, still include:
and the restoring unit is used for restoring the pile-up signal in the pulse signal.
Further, the recovery unit further includes:
the first processing unit is used for dividing the ideal waveform corresponding to the accumulated signal by a normalization coefficient to obtain a recovered waveform, wherein the normalization coefficient is a normalization coefficient for normalizing the pulse signal according to an ideal waveform model;
and a replacing unit for replacing the waveform of the accumulation signal with the recovery waveform.
Preferably, the method further comprises the following steps:
an acquisition unit configured to acquire a pulse signal;
the judging unit is used for judging the energy of the pulse signal;
and the selection unit is used for selecting the ideal waveform model matched with the pulse signal energy.
The pulse signal processing method and the device have the following beneficial effects:
by selecting the ideal waveform models corresponding to different energy sections of the pulse signals and carrying out normalization processing on the pulse signals, the waveforms and the noises can be more accurately discriminated, and the signal-noise separation performance is improved; through dividing the pulse signal pole after will normalizing, further divide to the different degree of piling up, through square difference value and square difference threshold value comparison, if pile up at first section pulse signal, it is serious to explain the signal piles up, abandons, if pile up at second section pulse signal, keep piling up the not serious signal and resume to make equipment can detect under the high dose rate environment under certain domain.
Drawings
Other features, objects and advantages of the present application will become more apparent upon reading of the following detailed description of non-limiting embodiments thereof, made with reference to the accompanying drawings in which:
FIG. 1 is a schematic flow chart diagram illustrating one embodiment of a pulse signal processing method of the present invention;
FIG. 2 is a flow diagram illustrating one embodiment of the present invention for determining whether a waveform needs to be recovered;
FIG. 3 is a schematic diagram of an ideal waveform model in one embodiment of the invention;
FIG. 4 is a schematic diagram of a piled-up pulse signal in one embodiment of the invention;
FIG. 5 is a schematic flow chart of another embodiment of the present invention for recovering a waveform;
FIG. 6 is a schematic flow chart of the pulse signal selection ideal waveform model of the present invention;
FIG. 7 is a schematic structural diagram of an embodiment of a pulse signal processing apparatus of the present invention;
FIG. 8 is a schematic structural diagram of an embodiment of a second judging unit according to the present invention;
fig. 9 is a schematic structural diagram of another embodiment of the pulse signal processing apparatus of the present invention;
FIG. 10 is a schematic diagram of the recovery unit of the present invention;
fig. 11 is a schematic diagram of the structure of the present invention for selecting an ideal waveform.
Detailed Description
The present application will be described in further detail with reference to the following drawings and examples. It is to be understood that the specific embodiments described herein are merely illustrative of the relevant invention and not restrictive of the invention. It should be noted that, for convenience of description, only the portions related to the present invention are shown in the drawings.
It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict. The present application will be described in detail below with reference to the embodiments with reference to the attached drawings.
Fig. 1 shows a schematic flowchart of an embodiment of the present invention, and a pulse signal processing method according to the present invention includes:
step S101, carrying out amplitude normalization processing on the pulse signal according to an ideal waveform model of the pulse signal to obtain a normalized pulse signal;
step S102, dividing the normalized pulse signal into at least two sections, wherein the first section of pulse signal is the pulse signal from the pulse signal to the peak value, and the second section of pulse signal is the residual pulse signal except the first section of pulse signal in the pulse signal;
step S103, calculating a first variance value corresponding to the first section of pulse signal and a second variance value corresponding to the second section of pulse signal;
step S104, judging whether the first variance difference value is smaller than a first threshold value; if yes, executing step S105;
and step S105, judging whether the waveform of the pulse signal needs to be subjected to waveform recovery or not according to the second variance difference value and the second threshold value.
It should be noted that:
in this embodiment, the ideal waveform model and the standard value of the square difference threshold are written into the chip in advance, and the corresponding ideal waveform model is called according to the energy value of the corresponding pulse signal after the chip is powered on. The square difference threshold table has m values (the m values are calculated in advance on a computer to obtain a square difference value, which is recorded as delta1δ2…δmThe size of m depends on the granularity of the alignment process, and a smaller value of δ indicates that more precise discrimination is desired, and conversely, a lower requirement). According to the detected requirements, the gear is divided into different gears in advance, for example: precision/fuzzy contrast gears, each having a corresponding squared difference threshold table. And selecting a square difference threshold value table according to the detection of the pulse signal energy segment and the estimation of the pile-up degree in advance.
It should be noted here that the ideal waveform model can be obtained by any modeling algorithm, and exemplarily, one implementation of the ideal waveform model establishment is as follows:
collecting a large number of waveform samples through a computer, manually supervising, and selecting m sampling samples of ideal waveforms (under an ideal detection environment and a low counting rate condition, fixing a crystal and an amplification forming circuit of a detector, outputting a signal peak reaching time T1, a signal decay time T2, an amplitude threshold lower limit A1, an amplitude threshold upper limit A2, taking sampling signals with all parameter differences within 10% as ideal waveforms, which can be of various amplitudes), wherein the length of each waveform is n points;
dividing a sampling sample of an ideal waveform into at least one energy segment according to an energy threshold, establishing an ideal waveform model (three are taken as an example, a high energy segment, a middle energy segment and a low energy segment) corresponding to each energy segment by dividing the energy segment, establishing the ideal waveform model by taking the high energy segment as an example, and establishing the ideal waveform model by taking other energy segments as similar;
selecting a waveform X with the amplitude median in the energy segment from m ideal waveform samplesiWave points are noted as (x)i1,xi2,…xin) The other waveform will have an amplitude normalization coefficient A compared with it1,A2…AmWhereinAi＝1；
Converting other m-1 waveforms into amplitude A which is the same as the amplitude median waveform through an amplitude normalization coefficientlxljWherein A islIs an amplitude normalization coefficient, l is 1-m; x is the number ofljWaveform points of the other m-1 waveforms; x is the number ofijWaveform points of an amplitude median waveform; the sum of the squares of the other m-1 waveforms and the magnitude median waveform is calculated using a small two multiplication as follows:
wherein A islIs a discrete fixed value, and L is coupled to AlDerivative that its variables should be continuously variable, where A will belIs marked as AkWhen the square difference sum L is minimum, obtaining the normalization coefficient A of other m-1 waveforms and the waveform of the amplitude mediank：
Multiplying waveform points of m sampled waveforms by an amplitude normalization coefficient AkObtaining m waveforms S with the same amplitude, wherein the same length is n points, and the waveform points are marked as SijI is 1 to m, j is 1 to n; let the waveform to be created be Y and the waveform point be (Y)1,y2…yn) And calculating the sum of the square differences of the normalized m sampling waveforms S and the waveform Y to be established by using a least square method:
due to yjAre discrete fixed values, L2 vs. yjThe derivation is a continuous process that is,thus will yjAnd (5) recording as yk (k is 1-m), and when the sum of squared differences of all waveforms is minimum, obtaining a model of an ideal waveform:
and processing the waveforms of other energy sections according to the same method, calculating an ideal waveform model, and storing the ideal waveform model in a computer for processing the later pulse signals.
In step S101, an amplitude normalization process is performed on the pulse signal according to the selected ideal waveform model to obtain a normalized pulse signal, where the normalization process is to multiply the amplitude of the pulse signal by a normalization coefficient according to the amplitude of the ideal waveform. For example, the following steps are carried out: for an ideal waveform corresponding to a high-energy-band pulse signal, for example, the amplitude range of the ideal waveform is 10V-20V, the median of the amplitude is 15V, the median of the amplitude is 20V-30V, and the median of the amplitude is 25V, so that according to the coefficient 0.6 between the median of the amplitude of the pulse signal and the median of the amplitude of the ideal waveform, the amplitudes of other waveforms of the pulse signal are all multiplied by 0.6 to obtain normalization, and the generated pulse signal is the normalized pulse signal.
In step S102, the normalized pulse signal is divided into at least two segments, and a pulse signal of a photodetector is received, where the pulse signal is a current pulse signal output by each of the current pulse signals formed on the anode after the photoelectrons are multiplied and accelerated by the photomultiplier tube, and the pulse signal starts to rise to a peak value with the passage of time from the beginning of generation, starts to fall back again, rises to another peak value again after falling for a period of time, and then falls again … …. The pulse signal generated from the beginning to the first peak as the first segment of pulse signal can also be called as a rising edge segment, and the rest as the second segment of pulse signal can also be called as a falling edge segment. The purpose is that if the pulse signal is piled up at the rising edge section from the beginning of generation, the pulse signal is serious in pile up, and the pulse signal is definitely piled up at the back, so that the pulse signal needs to be abandoned, the judgment and the processing of the signal are convenient, and the piled-up signal at the falling edge section is generally considered to be recoverable, so that the waste of the pulse signal is avoided, and the device can be imaged under the high dose rate environment. Here, the second-segment signal, that is, the falling edge segment, may be further divided artificially according to the current accumulation degree and the expected value of the signal detection by the tester, and the smaller the number of divided segments, the more accurate the signal detection is, the better the effect after recovery is, but the smaller the number of remaining signals satisfying the condition is, the smaller the detectable imaging probability is, and therefore, the division degree of the second-segment pulse signal is artificially determined according to the accumulation degree. In the embodiment, the pulse signal is divided into at least two sections, so that the balance between the accumulation degree and the expected value of the detection imaging is conveniently realized, and the signal which is not accumulated seriously is reserved and recovered, so that the equipment can detect the imaging under a high-dose environment in a certain range.
The method comprises the following steps:
the total time of a group of pulse signals is 100ns, a peak value is generated for the first time at the position of 20ns, therefore, the 20ns is divided into a rising edge section before and a falling edge section after the 20ns, the falling edge section is divided into one section every 20ns according to the judgment of signal waveforms, the falling edge section is divided into four sections, after the falling edge section is divided into four sections for processing, if the imaging cannot be carried out, the falling edge section can be divided into three sections according to each 25 ns.
It should be noted that, in step 102, the second segment pulse signal is referred to the first segment pulse signal, and the "second" is not limited to one, and is defined to distinguish the first segment pulse signal, and the second variance value in steps 103 to 105 and the "second" in the second threshold value are also defined based on the difference from the "first".
In step S103, a squared difference value, i.e., a variance, corresponding to each segment of the pulse signal is calculated, and a variance between the waveform point of the pulse signal and the waveform point data corresponding to the ideal waveform is calculated.
The method comprises the following steps:
normalization coefficient based on the normalized pulse signal in step S101, here denoted as a*For example, the pulse signal is divided into three segments, which are sequentially recorded as a first segment pulse signal, a second segment pulse signal and a third segment pulse signal, the corresponding squared difference value is recorded as L1, L2 and L3, the calculation method of each squared difference value is the same, taking the calculation of L3 as an example: suppose an input signal Z, n points in length (Z)1,Z2,Z3,…,Zn) The corresponding ideal waveform data is assumed to be Y, and the length is also n points (Y)1,Y2,…Yn) Using the least squares method, then:
for variable A*Derivation is carried out, and when the value of L3 is minimum, A is obtained*A value of (d);
a is to be*Substituting into L3, the squared difference between the normalized input waveform and the ideal waveform is obtained.
According to the corresponding gear of the selected ideal waveform, inquiring the corresponding square difference threshold value to be represented as delta1,δ2…δmIt should be noted that the pulse signal is divided into several segments corresponding to the ideal waveform, and therefore, the square difference threshold, and the process is performed in advanceThe pulse signals are written into the chip first, and when the pulse signals are input, selection is performed according to the accumulation degree of the pulse signals and the expected value of the expected imaging. In step S104, it is determined whether the first variance difference is smaller than a first threshold, where the first threshold refers to a corresponding square difference threshold δ in the ideal waveform model1By analogy, the second threshold refers to a squared difference threshold δ2. It can be said here whether L1 is less than δ1. In step S105, if L1 is less than δ1If the first section of pulse signal is deemed to have no pile-up signal, the pile-up signal should be at a certain position of the second section of pulse signal, the pile-up signal waveform at the second section is deemed to be recoverable, and whether the waveform needs to be recovered is determined according to the second variance difference value and the second threshold value in sequence.
Compared with the prior art, the pulse signal is selected to be the ideal waveform of the corresponding energy section, the problem that the waveforms of various energy sections are slightly different is solved, the selection of the ideal waveform by the energy sections can be expanded to any section in actual work, and the method has better significance for accurately discriminating the waveforms. Compared with the existing method for processing the piled signals by establishing an ideal waveform model comparison, the method is suitable for the situation of serious signal pile-up in the high dose rate environment, and a tester can recover the piled signal waveform in the pulse signals according to the pile-up degree and the expected value of the imaging effect, so that the equipment can image in the high dose rate environment.
Further, in step S105, if the first variance difference is not smaller than the first threshold, there is a pile-up signal in the first segment of the pulse signal, that is, there is pile-up in the rising edge segment, and the pulse signal is discarded.
In some embodiments, as shown in fig. 2, when the normalized pulse signal is divided into two segments, determining whether the waveform of the pulse signal needs to be subjected to waveform restoration according to the second variance difference value and the second variance difference threshold value includes:
step 1052, if the second variance is smaller than the second threshold, the pulse signal is an ideal waveform and does not need to be recovered;
step S1053, if the second variance difference is not smaller than the second threshold, the second segment of pulse signal has a pile-up signal, and the pile-up signal is subjected to waveform restoration.
In a preferred embodiment, when the normalized pulse signal is divided into three or more segments, the second segment of pulse signal includes at least two segments, and in step S103, the squared difference corresponding to each segment in the second segment of pulse signal is simultaneously calculated;
judging whether the waveform of the pulse signal needs to be subjected to waveform recovery or not according to the squared difference value corresponding to each section in the second section of the pulse signal and the corresponding threshold value, and further comprising:
sequentially judging whether the square difference value corresponding to each section in the second section of pulse signal is smaller than the corresponding square difference threshold value;
if so, the pulse signal is an ideal waveform and does not need to be subjected to waveform recovery;
if not, the method is divided into two cases:
firstly, if the variance difference value corresponding to the last section in the second section of pulse signal is not less than the corresponding threshold value, then the last section in the second section of pulse signal has a pile-up signal, and then the pile-up signal is subjected to waveform recovery;
and secondly, if the second section of pulse signals are generated according to the time sequence of the pulse signals in sequence and before the square difference corresponding to the last section is judged, the pulse signals are discarded when the difference of squares of any previous section is not less than the corresponding threshold value.
It should be noted that the waveform recovery of the pulse signal can be understood as the waveform recovery of the section where the pile-up signal is located, when the pulse signal is divided into two sections, the pile-up is located in the second section of the pulse signal, that is, the falling edge section, and the waveform of the signal in the whole falling edge section is recovered; when the pulse signal is divided into three or more sections, the accumulated signal is located in the last section, and the waveform of the pulse signal in the last section is recovered.
For example, the following steps are carried out: after a pulse signal is input, normalization processing is carried out on the pulse signal according to a corresponding ideal waveform model, and the ideal waveformThe model is shown in FIG. 3, corresponding to a squared difference threshold δ1，δ2，δ3，δ4. The normalized pulse signal is divided into four segments, including a first segment pulse signal (rising edge segment) and a second segment pulse signal (falling edge segment), where the falling edge segment pulse signal is divided into three segments as shown in fig. 4, and the corresponding squared differences are respectively denoted as L1, L2, L3, and L4. When L1 is less than δ1When it is, L2 and delta are determined2If L2 is not less than δ2The pulse signal is discarded if L2 is less than δ2Then determine L3 and δ3The size of (d); if L3 is not less than δ3The pulse signal is discarded if L3 is less than δ3Then determine L4 and δ4The size of (d); if L4 is not less than δ4If so, accumulating the last section of the pulse signal, and considering that the waveform of the fourth section needs to be recovered; if L4 is less than δ4Then the waveform is considered to be an ideal waveform.
In this embodiment, the second segment of pulse signal is divided into at least two segments, it is determined that there is no accumulation in the first segment of pulse signal based on the first variance difference and the first threshold, that is, there is no accumulation signal in the rising edge segment, the accumulation signal is located in the second segment of pulse signal, that is, the falling edge segment is divided into multiple segments, and whether the pulse signal needs to be recovered or not is determined according to the variance difference of each segment and the corresponding threshold, so that the accumulation signal can be recovered more accurately and effectively.
Further, another embodiment of the present application provides a method for waveform recovery of a pulse signal, as shown in fig. 5, including:
step S106, dividing the ideal waveform amplitude corresponding to the accumulation signal by a normalization coefficient to obtain a recovered waveform, wherein the normalization coefficient is a normalization coefficient for normalizing the pulse signal according to an ideal waveform model;
in step S107, the accumulated signal waveform is replaced with a restored waveform.
The ideal waveform amplitude corresponding to the pulse signal section where the accumulation signal is located is divided by the normalization coefficient, the processed ideal waveform is moved to the accumulation signal waveform, the pulse signal which is not seriously accumulated is extracted and recovered, and therefore the device achieves detection imaging within a certain range of the domain.
In some embodiments, before performing amplitude normalization on the pulse signal according to the ideal waveform model of the pulse signal to obtain a normalized pulse signal, as shown in fig. 6, the method further includes:
step S1011, acquiring a pulse signal;
step S1012, determining the energy of the pulse signal;
in step S1013, an ideal waveform model matching the pulse signal energy is selected.
The output signals of different photoelectric detectors are different, in the embodiment, the same photoelectric detector inputs the acquired pulse signals of the photoelectric detectors, different energy sections exist in the pulse signals, and the waveforms of the energy sections are different.
For example, a pulse signal is input and judged to belong to a high-energy section, and therefore, an ideal waveform model of the corresponding high-energy section is selected. It should be noted that the signal of the high-energy section may also be a superposition of two or more signals of the low-energy section, and although the signal of the high-energy section is obtained in the determination, the waveform of the high-energy section formed by superposition is obviously different from the waveform of a single signal of the high-energy section, so that such signals of the high-energy section need to be discarded.
In summary, the present application selects a corresponding ideal waveform model for the energy of the pulse signal, and divides the waveform itself into several parts according to the rising edge and the falling edge. The method completely discards the situation that the waveform is piled up at the rising edge part, and considers the situation that each block of the falling edge is piled up respectively, so that the tolerance of the piling up situation and the degree of signal recovery can be flexibly selected completely according to the actual situation of an operator. Some signals which are relatively piled up and not serious can be extracted and recovered, so that the device can achieve detection imaging within a certain range of the domain.
In a second aspect, another embodiment of the present application provides a pulse signal processing apparatus, as shown in fig. 7, including:
the processing unit 110 is configured to perform amplitude normalization processing on the pulse signal by using the ideal waveform model to obtain a normalized pulse signal;
the segmenting unit 120 is configured to segment the normalized pulse signal into at least two segments, where a first segment of the pulse signal is a pulse signal generated from the pulse signal to a peak value, and a second segment of the pulse signal is a pulse signal remaining in the pulse signal except the first segment of the pulse signal;
a calculating unit 130, configured to calculate a first variance difference value corresponding to the first segment of pulse signal and a second variance difference value corresponding to the second segment of pulse signal;
a first determining unit 140, configured to determine whether the first variance difference is smaller than a first threshold, and if so, activate a second determining unit;
the second determining unit 150 determines whether the waveform of the pulse signal needs to be restored according to the second variance difference and the second threshold.
The pulse signal processing apparatus in this embodiment processes a pulse signal output by a front-end device, normalizes the pulse signal through a pre-built-in ideal waveform model to obtain a normalized pulse signal, and a segmentation unit divides the normalized pulse signal into at least two segments, where a first segment of the pulse signal is also called a rising edge segment, a second segment of the pulse signal is also called a falling edge segment, the falling edge segment can be further divided into multiple segments according to a stacking degree and an expected value of detection imaging, a variance difference value of each segment is calculated, and if the variance difference value of the first segment of the pulse signal is not less than a first threshold value, it is indicated that the rising edge segment has a stack, and the pulse signal is discarded; and if the first mean square deviation value is smaller than the first threshold value, judging whether the pulse signal needs waveform recovery according to the second mean square deviation value of the second section of pulse signal.
In some embodiments, when the normalized pulse signal is divided into three or more segments, the second segment pulse signal includes at least two segments, as shown in fig. 8, the second determining unit 105 further includes:
and the at least two judging subunits are used for sequentially judging whether the square difference value corresponding to each section in the second section of pulse signal is smaller than the corresponding threshold value.
If so, the pulse signal is an ideal waveform and does not need to be subjected to waveform recovery;
if not, the method is divided into two cases:
firstly, if the variance difference value corresponding to the last section in the second section of pulse signal is not less than the corresponding threshold value, the last section of pulse signal in the second section of pulse signal has accumulated signals and needs to be subjected to waveform recovery;
and secondly, if the second section of pulse signals are generated according to the time sequence of the pulse signals in sequence and before the square difference corresponding to the last section is judged, the pulse signals are discarded when the difference of squares of any previous section is not less than the corresponding threshold value.
Further, as shown in fig. 9, according to the result of the second determination unit 105, a restoration unit 160 is further included for restoring the pile-up signal waveform in the pulse signal.
Specifically, as shown in fig. 10, the recovery unit 160 includes:
the first processing unit 161 is configured to divide an ideal waveform corresponding to the stacked signal by a normalization coefficient to obtain a recovered waveform, where the normalization coefficient is a normalization coefficient for normalizing the pulse signal according to an ideal waveform model;
and a replacing unit 162 for replacing the pile-up signal with the restored waveform.
In some embodiments, before the processing unit performs normalization according to the ideal waveform model corresponding to the pulse signal, the method further includes, as shown in fig. 11:
an acquisition unit 111 for acquiring a pulse signal;
a judging unit 112, configured to judge an energy level of the pulse signal;
and a selecting unit 113 for selecting an ideal waveform model matched with the pulse signal energy.
The acquisition unit inputs the acquired pulse signals, the judgment unit judges the energy of the pulse signals, and the corresponding ideal waveform model is selected according to the energy of the pulse signals, so that the waveforms and the noises can be more accurately discriminated, and the signal-noise separation performance is improved. In the embodiment, the pulse signals are segmented to judge the size of the square difference value and the threshold value, the position of the accumulated signals is determined and recovered, the pulse signals are divided in energy, the waveform recovery efficiency is improved, and the detection imaging quality of equipment is improved.
The processes described above with reference to fig. 1-2 and 5-6 may be implemented as computer software programs, according to embodiments of the present disclosure. For example, embodiments of the present disclosure include a computer program product comprising a computer program tangibly embodied on a machine-readable medium, the computer program comprising program code for performing the methods of fig. 1-4. In such an embodiment, the computer program may be downloaded and installed from a network via the communication section, and/or installed from a removable medium.
The flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
As another aspect, the present application also provides a computer-readable storage medium, which may be the computer-readable storage medium included in the foregoing device in the foregoing embodiment; or it may be a separate computer readable storage medium not incorporated into the device. The computer readable storage medium stores one or more programs for use by one or more processors in performing the pulse signal processing described herein.
The above description is only a preferred embodiment of the application and is illustrative of the principles of the technology employed. It will be appreciated by a person skilled in the art that the scope of the invention as referred to in the present application is not limited to the embodiments with a specific combination of the above-mentioned features, but also covers other embodiments with any combination of the above-mentioned features or their equivalents without departing from the inventive concept. For example, the above features may be replaced with (but not limited to) features having similar functions disclosed in the present application.
Claims (10)
1. A pulse signal processing method, comprising:
carrying out amplitude normalization processing on the pulse signals according to the ideal waveform models corresponding to different energy sections to obtain normalized pulse signals, wherein the normalization process is as follows: the amplitude of the pulse signal is processed according to the amplitude of the ideal waveform multiplied by a normalization coefficient;
dividing the normalized pulse signal into at least two sections, wherein the first section of pulse signal is a pulse signal generated from the pulse signal to a peak value, and the second section of pulse signal is a pulse signal remained in the pulse signal except the first section of pulse signal;
calculating a first variance difference value corresponding to the first section of pulse signals and a second variance difference value corresponding to the second section of pulse signals;
judging whether the first variance difference value is smaller than a first threshold value;
if yes, judging whether the waveform of the second section of pulse signal needs to be subjected to waveform recovery or not according to the second variance difference value and a second threshold value;
if the waveform needs to be recovered, the recovered waveform is obtained according to the ideal waveform corresponding to the stacked signal, and the stacked signal waveform is replaced by the recovered waveform.
2. The pulse signal processing method according to claim 1,
when the normalized pulse signal is divided into two sections,
the determining whether the waveform of the second segment of pulse signal needs to be waveform restored according to the second variance difference value and the second variance difference threshold value includes:
if the second variance difference value is smaller than the second threshold value, the pulse signal is an ideal waveform and does not need to be recovered;
and if the second variance difference value is not smaller than the second threshold value, a pile-up signal exists in the second section of pulse signal, and the waveform of the pile-up signal is recovered.
3. The method according to claim 1, wherein when the normalized pulse signal is divided into three or more segments, the second segment of the pulse signal includes at least two segments, and further comprising:
sequentially judging whether the square difference value corresponding to each section in the second section of pulse signals is smaller than the corresponding square difference threshold value;
if so, the pulse signal is an ideal waveform and does not need to be subjected to waveform recovery;
and if only the squared difference value corresponding to the last section in the second section of pulse signals is not smaller than the corresponding threshold value, accumulating signals exist in the last section of the second section of pulse signals, and waveform recovery is carried out on the accumulated signals.
4. A pulse signal processing method according to claim 2 or 3, wherein the waveform recovery method comprises:
dividing the ideal waveform amplitude corresponding to the accumulation signal by a normalization coefficient to obtain a recovered waveform, wherein the normalization coefficient is a normalization coefficient for normalizing the pulse signal according to an ideal waveform model;
replacing the waveform of the pile-up signal with the restored waveform.
5. The method according to claim 1, further comprising, before the amplitude normalizing the pulse signal according to the ideal waveform model of the pulse signal to obtain a normalized pulse signal:
acquiring the pulse signal;
judging the energy of the pulse signal;
and selecting an ideal waveform model matched with the pulse signal energy.
6. A pulse signal processing apparatus, characterized by comprising:
the processing unit is used for carrying out amplitude normalization processing on the pulse signals according to the ideal waveform models corresponding to different energy sections to obtain normalized pulse signals, wherein the normalization process is as follows: the amplitude of the pulse signal is processed according to the amplitude of the ideal waveform multiplied by a normalization coefficient;
the segmentation unit is used for dividing the normalized pulse signal into at least two segments, wherein the first segment of pulse signal is a pulse signal generated from the pulse signal to a peak value, and the second segment of pulse signal is a pulse signal remained in the pulse signal except the first segment of pulse signal;
the calculating unit is used for calculating a first variance difference value corresponding to the first section of pulse signals and a second variance difference value corresponding to the second section of pulse signals;
the first judgment unit is used for judging whether the first variance difference value is smaller than a first threshold value or not, and if so, the second judgment unit is activated;
the second judging unit is used for judging whether the waveform of the pulse signal needs to be subjected to waveform recovery or not according to the second variance difference value and a second threshold value; if the waveform needs to be recovered, the recovered waveform is obtained according to the ideal waveform corresponding to the stacked signal, and the stacked signal waveform is replaced by the recovered waveform.
7. The apparatus according to claim 6, wherein when the normalized pulse signal is divided into three or more segments, the second segment pulse signal includes at least two segments, and further comprising:
at least two judging subunits, configured to sequentially judge whether a difference of squares corresponding to each segment in the second segment of pulse signal is smaller than a corresponding threshold; if so, the pulse signal is an ideal waveform and does not need to be subjected to waveform recovery;
and if the difference of squares corresponding to the last section in the second section of pulse signals is not smaller than the corresponding threshold value, accumulating signals exist in the last section of the second section of pulse signals, and waveform recovery is carried out on the accumulating signals.
8. The pulse signal processing device according to claim 6, further comprising:
and the restoring unit is used for restoring the pile-up signal in the pulse signal.
9. The pulse signal processing apparatus according to claim 8, wherein the recovery unit further comprises:
the first processing unit is used for dividing an ideal waveform corresponding to the accumulated signal by a normalization coefficient to obtain a recovered waveform, wherein the normalization coefficient is a normalization coefficient for normalizing the pulse signal according to an ideal waveform model;
and a replacing unit that replaces the waveform of the pile-up signal with the restored waveform.
10. The pulse signal processing device according to claim 6, further comprising:
an acquisition unit configured to acquire the pulse signal;
the judging unit is used for judging the energy of the pulse signal;
and the selection unit is used for selecting an ideal waveform model matched with the pulse signal energy.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911394385.3A CN111245409B (en) | 2019-12-30 | 2019-12-30 | Pulse signal processing method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911394385.3A CN111245409B (en) | 2019-12-30 | 2019-12-30 | Pulse signal processing method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111245409A CN111245409A (en) | 2020-06-05 |
| CN111245409B true CN111245409B (en) | 2021-12-24 |
Family
ID=70879450
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911394385.3A Active CN111245409B (en) | 2019-12-30 | 2019-12-30 | Pulse signal processing method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111245409B (en) |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2621736C (en) * | 2004-09-16 | 2018-07-31 | Southern Innovation International Pty Ltd | Method and apparatus for resolving individual signals in detector output data |
| WO2006034585A1 (en) * | 2004-09-28 | 2006-04-06 | UNIVERSITé DE SHERBROOKE | Method and system for low radiation computed tomography (ct) |
| JP2011519415A (en) * | 2008-03-31 | 2011-07-07 | サザン イノヴェーション インターナショナル プロプライアトリー リミテッド | Radiation imaging using individual signal resolution |
| CN102073059B (en) * | 2010-12-31 | 2013-05-22 | 华中科技大学 | Digital pileup waveform processing method and system |
| CN104101894B (en) * | 2013-04-07 | 2017-04-12 | 苏州瑞派宁科技有限公司 | Accumulation event detection and reconstruction method and system |
| CN103454671B (en) * | 2013-08-21 | 2016-01-06 | 中国人民解放军第二炮兵工程大学 | A kind of nuclear radiation pulse pile-up judgement based on high speed digital sample and bearing calibration |
| CN106170919B (en) * | 2013-12-11 | 2019-05-28 | 南方创新国际股份有限公司 | Method and apparatus for parsing the signal in data |
| CN107767427B (en) * | 2017-09-28 | 2021-05-07 | 东软医疗***股份有限公司 | Signal waveform recovery method and device |
| CN110274921B (en) * | 2019-01-17 | 2022-08-05 | 武汉大学 | Full-digital two-dimensional coincidence Doppler broadening system |
-
2019
- 2019-12-30 CN CN201911394385.3A patent/CN111245409B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN111245409A (en) | 2020-06-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109212597B (en) | Micro seismic source automatic positioning method based on deep belief network and scanning superposition | |
| US8860883B2 (en) | Method and apparatus for providing signatures of audio/video signals and for making use thereof | |
| CN108805269B (en) | Method for picking seismic facies arrival time based on LSTM recurrent neural network | |
| US10078009B2 (en) | Pulse processing circuit with correction means | |
| CN102073059B (en) | Digital pileup waveform processing method and system | |
| US7697749B2 (en) | Stereo image processing device | |
| CN104101894B (en) | Accumulation event detection and reconstruction method and system | |
| CN104698343A (en) | Method and system for judging power grid faults based on historical recording data | |
| KR101666521B1 (en) | Method and apparatus for detecting pitch period of input signal | |
| CN106908655A (en) | Photosignal peak-value detection method and system | |
| CN108549078B (en) | Cross-channel combination and detection method for radar pulse signals | |
| CN107767427B (en) | Signal waveform recovery method and device | |
| CN104699755A (en) | Intelligent multi-target comprehensive identification method based on data mining | |
| KR101932595B1 (en) | Image processing apparatus and method for detecting translucent objects in image | |
| CN109165603B (en) | Ship detection method and device | |
| CN111245409B (en) | Pulse signal processing method and device | |
| CN112327104A (en) | Fault detection and positioning method for power distribution network with distributed power supply | |
| US20040078196A1 (en) | Signal processing method and processor | |
| WO2017187972A1 (en) | Radiation detector and signal processing device for radiation detection | |
| CN102480455B (en) | Detecting method and detecting device for primary synchronization signals in long-term evolution system | |
| CN111122420A (en) | Synchronous identification method and device of multichannel pulse signals and particle detection system | |
| CN112052337A (en) | Target relation detection method, system and storage medium based on time-space correlation | |
| CN115952410A (en) | Landslide hazard detection system based on deep learning | |
| CN111103616A (en) | Gamma energy spectrum measuring method and device considering both energy resolution and detection efficiency | |
| US20060126763A1 (en) | Processing of a signal representing radiation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20230330 Address after: 100049 No. 19, Yuquanlu Road, Beijing, Shijingshan District Patentee after: INSTITUTE OF HIGH ENERGY PHYSICS, CHINESE ACADEMY OF SCIENCES Patentee after: Jinan Zhongke Nuclear Technology Research Institute Address before: 100049 No. 19 (b) hospital, Shijingshan District, Yuquanlu Road, Beijing Patentee before: INSTITUTE OF HIGH ENERGY PHYSICS, CHINESE ACADEMY OF SCIENCES |