CN103607205B - Signal processing method, apparatus and equipment - Google Patents
Signal processing method, apparatus and equipment Download PDFInfo
- Publication number
- CN103607205B CN103607205B CN201310566812.8A CN201310566812A CN103607205B CN 103607205 B CN103607205 B CN 103607205B CN 201310566812 A CN201310566812 A CN 201310566812A CN 103607205 B CN103607205 B CN 103607205B
- Authority
- CN
- China
- Prior art keywords
- signal
- decay period
- sampled result
- threshold value
- predetermined threshold
- 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
Abstract
The invention discloses a signal processing method, apparatus and equipment. The signal processing method comprises the following steps that: an amplified electric pulse signal is obtained; analog-to-digital conversion is carried out on the amplified electric pulse signal so as to obtain a digital signal; sampling is carried out on the digital signal according to a preset time interval; whether a sampling result is larger than a preset threshold value is determined; and if so, the sampling result is stored; and two sampling results that are larger than the preset threshold value and are first to be stored in a cycle of decay are utilized to identify time information of the cycle of decay. According to the invention, the electric pulse signal is converted into the digital signal; time information identification processing is carried out on the sampling result of the digital signal; and the first two sampling results that are larger than the preset threshold value in the cycle of decay are used for identifying time information of the cycle of decay. Therefore, there is no need to use a time processing circuit for time identification process; the structure of a positron emission tomography (PET) system is simplified; and the system cost is lowered. Moreover, the digital signals can be transmitted in parallel; the transmission speed is fast; the time delay error is small; and the precision is high.
Description
Technical field
The present invention relates to field of computer technology, more particularly to a kind of signal processing method, device and equipment.
Background technology
Positron emission tomography（Positron Emission Tomography, PET）Equipment is that medical treatment is led now
A kind of more advanced molecular imaging equipment in domain.In fluorodeoxyglucose, labelling radionuclide, as tracer, is put
The positron that penetrating property nucleic discharges, after motion one segment distance, is fallen into oblivion with the electronics with negative charge in surrounding
Go out, produce a pair of energy γ photon equal, in opposite direction.
Traditional PET device typically adopts closed-loop panel detector structure, and detector rings are by the assembly of multiple detector modules
Become.Each detector module is by multiple scintillation crystals and photomultiplier tube（Photomultiplier Tube, PMT）Composition
's.Bury in oblivion produced γ photon to be received by detector ring, decay after γ photon incides scintillation crystal produces light pulse letter
Number, light pulse signal is converted into electric impulse signal by PMT.
After the electric impulse signal of PMT output is filtered and amplifies, electric impulse signal is carried out signal processing, obtains γ light
The temporal information of son.Temporal information is mainly used in determining the moment receiving γ photon, and then carries out meeting the judgement of event.When one
The time difference that γ photon is reached with detector is less than time window set in advance（It is usually 8 to 12ns）Then it is assumed that this is to γ light
Son stems from same positron annihilation events, that is, produce one and meet event, this is then remembered to information produced by γ photon
Record；Otherwise, this is considered as then two single incidents to γ photon, and produced information is dropped.By the electricity that PMT is exported
Pulse signal is analyzed, the moment value that a pair of γ photon of differentiation is received, and judges whether to meet event, and judges to receive
Produce the position of the scintillation crystal of γ photon meeting event, and then obtain radionuclide concentration distribution in vivo.
At present, the identification to the temporal information of γ photon mainly adopts time-triggered protocol circuit, and time-triggered protocol circuit receives
The γ light period of the day from 11 p.m. to 1 a.m, the high-speed comparator in circuit produces a time pulse, is produced with the measurement of high-precision measure of time chip
Time pulse and the time difference between the nominal time of master clock, so that it is determined that receiving the temporal information of γ photon.
Those skilled in the art when electric impulse signal PMT being exported using said method is processed, find just like
Lower shortcoming：
Time-triggered protocol circuit needs using high-speed comparator and high-precision measure of time chip, and circuit cost is high, often
The passage of individual PET system is required for determining the temporal information of γ photon using time-triggered protocol circuit, leads to whole PET system structure
Complexity, cost is high.
Content of the invention
In view of this, the invention provides a kind of signal processing method, device are equipment, by the electric impulse signal after amplifying
It is converted into digital signal, digital signal samples result is identified process, need not be using multiple complicated time-triggered protocol circuit
Realize temporal information identifying processing respectively.
A kind of signal processing method, methods described includes：
Obtain the electric impulse signal after amplifying, described electric impulse signal is carried out obtaining after opto-electronic conversion by light pulse signal；
Electric impulse signal after amplifying is carried out analog digital conversion and obtains digital signal；
Described digital signal is sampled according to prefixed time interval；
Judge whether sampled result is more than predetermined threshold value, if it is, storing described sampled result；
Identify this decay using the sampled result that two storing at first in a decay period are more than predetermined threshold value
One temporal information in cycle.
Optionally, methods described also includes：
Identify an energy information of this decay period using all sampled result being stored in a decay period.
Optionally, described all sampled result being stored in a decay period are utilized to identify one of this decay period
Energy information includes：
All sampled result in one decay period are sued for peace summation of must sampling, will sampling summation and time interval product
Energy information as this decay period；
Or,
The product of each sampled result being stored in a decay period and time interval as a sub- energy value,
Every sub- energy value is sued for peace as an energy information in a decay period.
Optionally, described using two stored at first in a decay period be more than predetermined threshold value sampled result
Identify that a temporal information of this decay period includes：
Calculate two storing at first and be more than the double sampling moment of predetermined threshold value and being compared to of the difference of temporal information
For the first ratio；
Calculate the two sub- energy values of double sampling being more than predetermined threshold value storing at first is used for the second ratio；
The difference that two storing at first are more than between the double sampling moment of predetermined threshold value is prefixed time interval；
Equal with the second ratio using the first ratio, first sampling instant is calculated according to prefixed time interval and believes with the time
The difference of breath is as initial time section；
The sampling instant of first sampled result of calculating and the difference between the demarcation moment are as total period；
Calculate the difference between total period and initial time section as a temporal information of this decay period.
Optionally, the described electric impulse signal obtaining after amplifying includes：
Obtain the electric impulse signal after preposition amplification；
Or,
Obtain the preposition electric impulse signal amplifying after amplifying with variable gain.
A kind of signal processing apparatus, described device includes：
Signal acquiring unit, for obtaining the electric impulse signal after amplifying, described electric impulse signal is entered by light pulse signal
Obtain after row opto-electronic conversion；
AD conversion unit, obtains digital signal for the electric impulse signal after amplifying is carried out analog digital conversion；
Sampling unit, for sampling according to prefixed time interval to described digital signal；
Identifying unit, for judging whether sampled result is more than predetermined threshold value, if it is, store described sampled result；
Time recognition unit, for being more than adopting of predetermined threshold value using two storing at first in a decay period
Sample result identifies a temporal information of this decay period.
Optionally, described device also includes：
Thin interbed unit, for identifying this decay period using all sampled result being stored in a decay period
An energy information.
Optionally, described Thin interbed unit includes：
First Thin interbed subelement, for summation of must sampling that all sampled result in a decay period are sued for peace, will
The product of sampling summation and time interval is as the energy information of this decay period；
Or,
Second Thin interbed subelement, for each sampled result of being stored in a decay period and time interval
Product as a sub- energy value, every sub- energy value is sued for peace as an energy information in a decay period.
Optionally, described time recognition unit includes：
First computation subunit, for calculate store at first the two double sampling moment being more than predetermined threshold value and when
Between information difference be used for the first ratio；
Second computation subunit, for calculating the two sub- energy values of double sampling being more than predetermined threshold value storing at first
Be used for the second ratio；
3rd computation subunit, for two differences being more than between the double sampling moment of predetermined threshold value storing at first
It is worth for prefixed time interval；
4th computation subunit, for equal with the second ratio using the first ratio, calculates the according to prefixed time interval
The difference of one sampling instant and temporal information is as initial time section；
5th computation subunit, for calculating the sampling instant of first sampled result and demarcating the difference work between the moment
For total period；
6th computation subunit, for calculating the difference between total period and initial time section as this decay period
One temporal information.
Optionally, described signal acquiring unit includes：
First signal acquisition subelement, for obtaining the electric impulse signal after preposition amplification；
Or,
Secondary signal obtains subelement, for obtaining the preposition electric impulse signal amplifying after amplifying with variable gain.
A kind of signal handling equipment, described equipment includes：
Photomultiplier tube, amplifying circuit, modulus conversion chip and programmable logic device are sequentially connected；
Photomultiplier tube, for being converted into electric impulse signal by the light pulse signal being received；
Amplifying circuit, for being amplified to described electric impulse signal；
Modulus conversion chip, obtains digital signal for the electric impulse signal after amplifying is carried out analog digital conversion；
Programmable logic device, for sampling according to prefixed time interval to described digital signal；Judge sampled result
Whether it is more than predetermined threshold value, if it is, storing described sampled result；To the sampled result being stored in a decay period
The time of carrying out identifying processing and/or Thin interbed are processed.
As shown in the above, the method have the advantages that：
The invention discloses a kind of signal processing method, device and equipment, obtain the electric impulse signal after amplifying, described electricity
Pulse signal is carried out obtaining after opto-electronic conversion by light pulse signal；Electric impulse signal after amplifying is carried out analog digital conversion number
Word signal；Described digital signal is sampled according to prefixed time interval；Judge whether sampled result is more than predetermined threshold value,
If it is, storing described sampled result；It is more than adopting of predetermined threshold value using two storing at first in a decay period
Sample result identifies a temporal information of this decay period, the light pulse signal being received is converted into electric impulse signal, then turns
Change digital signal into, temporal information identifying processing is carried out to the sampled result of digital signal, utilizes in a decay period
Two first storing identify a temporal information of this decay period more than the sampled result of predetermined threshold value, need not adopt the time
Process circuit carries out time identifying processing, simplifies the structure of PET system, reduces system cost, and digital signal can be parallel
Transmission, transmission speed is fast, and time identifying processing realizes process simply, and delay time error is little, high precision.
Brief description
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
Have technology description in required use accompanying drawing be briefly described it should be apparent that, drawings in the following description be only this
Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, acceptable
Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is a kind of present invention signal processing method embodiment one flow chart；
Fig. 2 is digital signal samples result schematic diagram of the present invention；
Fig. 3 is a kind of present invention signal processing method embodiment two flow chart；
Fig. 4 obtains decay period temporal information t method schematic diagram for the present invention；
Fig. 5 is a kind of present invention signal processing apparatus embodiment three structural representation；
Fig. 6 is a kind of present invention signal processing apparatus example IV structural representation；
Fig. 7 is a kind of present invention signal handling equipment embodiment five structural representation.
Specific embodiment
The invention discloses a kind of signal processing method, device are equipment, the electric impulse signal after amplifying is converted into counting
Word signal, carries out temporal information identifying processing, simplied system structure, reduces system cost to digital signal samples result.
Below in conjunction with the accompanying drawings the embodiment of the present invention is described in detail.
Embodiment one
Fig. 1 is a kind of present invention signal processing method embodiment one flow chart, and methods described includes：
Step 101：Obtain the electric impulse signal after amplifying, described electric impulse signal carries out opto-electronic conversion by light pulse signal
After obtain.
Produced by positron annihilation, γ photon is received by the crystal on detector, and γ photon decays generation in crystal can
See the light pulse signal of light, light pulse signal is through photomultiplier tube（Photomultiplier Tube, PMT）It is converted into electric pulse
Signal output., in different crystal, decay period is different, for example for γ photon：In bismuth germanium oxide（Bi4Ge3O12, BGO）In crystal,
The decay period of γ photon is 300ns.
Due to visible ray produced by γ Photon Decay very faint it is impossible between detection although through PMT receive after carry out
Amplify, but the signal after amplifying still is not enough to directly be detected the electric pulse in addition it is also necessary to after PMT is amplified and changed
Signal amplifies further and just can be detected.
Optionally, the described electric impulse signal obtaining after amplifying includes two kinds of possible embodiments：
The first possible embodiment：
Obtain the electric impulse signal after preposition amplification.
Electric impulse signal can be amplified to the scope that can be directly detected by preposition amplification.
The possible implementation of second：
Obtain the preposition electric impulse signal amplifying after amplifying with variable gain.
Not only preposition amplification is carried out to electric impulse signal, also variable gain amplification is carried out to electric impulse signal, can be effective
Reduce electric impulse signal in noise.
Step 102：Electric impulse signal after amplifying is carried out analog digital conversion and obtains digital signal.
Electric impulse signal after amplifying is adopted analog-digital converter（Analog-to-Digital Converter ADC）Or
Person's modulus conversion chip is converted into digital signal.
Step 103：Described digital signal is sampled according to prefixed time interval.
Using programmable gate array（Field-Programmable Gate Array, FPGA）To the numeral letter after conversion
Number sampled according to prefixed time interval, sample frequency can not be less than 40MHz.As a example BGO crystal, the decay week of γ photon
Phase is 300ns, and the sample frequency of FPGA is 100M, often receives a γ photon and can sample 30 times.
Step 104：Judge whether sampled result is more than predetermined threshold value, if it is, execution step 105；If not, execution
Step 107.
Step 105：Store described sampled result.
Preset a predetermined threshold value in FPGA, when sampled value is more than predetermined threshold value, show to receive γ
Light pulse signal produced by photon, that is, receive an event, stores described sampled result.
Step 106：Known using the sampled result that two storing at first in a decay period are more than predetermined threshold value
One temporal information of this decay period other.
After light pulse signal produced by one γ Photon Decay cycle is converted into electric impulse signal, through analog digital conversion
The digital signal obtaining, when being sampled according to prefixed time interval, can collect multiple sampled result.Because γ photon declines
During change, the process releasing energy is one and first increases the process reducing again, then the value of sampled result is also first to increase to reduce again.From
First sampled result being more than predetermined threshold value, the sampled result being more than pre-determined threshold to last, owns in the meantime
Continuous sampled result be all the sampled result that a decay period is stored, as shown in Fig. 2 a decay period one coexists
23 sampled result are stored up, the value of sampled result is according to the sequencing sampled, first increases and then decreases.
Time identifying processing is carried out to sampled result：
Obtain γ photon more than two sampled result of predetermined threshold value and decline according to storing at first in a decay period
The concrete temporal information becoming.The temporal information obtaining is exported to meeting datatron, meeting datatron can will be most recently received
Two are carried out meeting time judgement by the temporal information that different signal processing apparatus send, when the time difference of two temporal informations
Less than time window set in advance（It is usually 8 to 12ns）Then it is assumed that this two temporal informations produce one and meet event, to institute
State and meet event and recorded；If the time difference of two temporal informations is more than time window set in advance, two time letters
Breath is respectively two single events, abandons described temporal information.The number of times meeting event according to meeting datatron record is
Follow-up process provides effective data foundation.
Step 107：Abandon described sampled result.
When sampled value is less than predetermined threshold value it is believed that current sampled value is noise, abandon described sampled result.
As shown in the above, the method have the advantages that：
The invention discloses a kind of signal processing method, device and equipment, obtain the electric impulse signal after amplifying, described electricity
Pulse signal is carried out obtaining after opto-electronic conversion by light pulse signal；Electric impulse signal after amplifying is carried out analog digital conversion number
Word signal；Described digital signal is sampled according to prefixed time interval；Judge whether sampled result is more than predetermined threshold value,
If it is, storing described sampled result；It is more than adopting of predetermined threshold value using two storing at first in a decay period
Sample result identifies a temporal information of this decay period, the light pulse signal being received is converted into electric impulse signal, then turns
Change digital signal into, temporal information identifying processing is carried out to the sampled result of digital signal, utilizes in a decay period
Two first storing identify a temporal information of this decay period more than the sampled result of predetermined threshold value, need not adopt the time
Process circuit carries out time identifying processing, simplifies the structure of PET system, reduces system cost, and digital signal can be parallel
Transmission, transmission speed is fast, and time identifying processing realizes process simply, and delay time error is little, high precision.
Embodiment two
Fig. 3 is a kind of present invention signal processing method embodiment two flow chart, and compared with embodiment one, methods described is also wrapped
Include Thin interbed to process, methods described includes：
Step 301：Obtain the electric impulse signal after amplifying, described electric impulse signal carries out opto-electronic conversion by light pulse signal
After obtain.
Step 302：Electric impulse signal after amplifying is carried out analog digital conversion and obtains digital signal.
Step 303：Described digital signal is sampled according to prefixed time interval.
Step 304：Judge whether sampled result is more than predetermined threshold value, if it is, execution step 305；If not, execution
Step 308.
Step 305：Store described sampled result.
Step 301- step 305 is similar with embodiment one, the description of reference implementation example one, repeats no more here.
Step 306：Known using the sampled result that two storing at first in a decay period are more than predetermined threshold value
One temporal information of this decay period other.
Optionally, described using two stored at first in a decay period be more than predetermined threshold value sampled result
Identify that a temporal information of this decay period includes：
Calculate two storing at first and be more than the double sampling moment of predetermined threshold value and being compared to of the difference of temporal information
For the first ratio；
Calculate the two sub- energy values of double sampling being more than predetermined threshold value storing at first is used for the second ratio；
The difference that two storing at first are more than between the double sampling moment of predetermined threshold value is prefixed time interval；
Equal with the second ratio using the first ratio, first sampling instant is calculated according to prefixed time interval and believes with the time
The difference of breath is as initial time section；
The sampling instant of first sampled result of calculating and the difference between the demarcation moment are as total period；
Calculate the difference between total period and initial time section as a temporal information of this decay period.
Optionally, as shown in figure 4, formula can also be adopted（1）Calculate a temporal information identifying this decay period：
Wherein, m by the sampling instant of first sampled result and demarcates the prefixed time interval being comprised between the moment
Number, E1For the sub- energy value of first sampled result, E2For the sub- energy value of second sampled result, t0For between Preset Time
Every t is temporal information.
Step 307：Identify an energy of this decay period using all sampled result being stored in a decay period
Amount information.
23 sampled result that a decay period as shown in Figure 2 is gathered, between the time of each sampled result and sampling
Every between product as a sub- energy value, all 23 sub- energy values summations are obtained the energy letter of a decay period
Breath.
Optionally, described all sampled result being stored in a decay period are utilized to identify one of this decay period
Energy information has two kinds of possible embodiments：
The first possible embodiment：Each sampled result being stored in a decay period and time interval
Product, as a sub- energy value, every sub- energy value is sued for peace as an energy information in a decay period.
According to formula（2）Obtain energy information E：
Wherein, hiFor the height of each sampled result, t0For prefixed time interval, n is the number of sampled result.
The possible embodiment of second：All sampled result in one decay period are sued for peace summation of must sampling, will adopt
The product of sample summation and time interval is as the energy information of this decay period.
According to formula（3）Obtain energy information E：
Wherein, hiFor the height of each sampled result, t0For prefixed time interval, n is the number of sampled result.
The sampled result being stored in one decay period is carried out after Thin interbed process, obtained energy information can
To provide reliable data foundation for later process：
The energy letter that one γ Photon Decay cycle produced is calculated according to the sampled result that stored in a decay period
Breath, this is the energy information produced by γ Photon Decay that a PMT is received.
The matrix that one detector is made up of 11 crystal taking advantage of 11, takes advantage of 2 PMT to receive by 2 after each detector
Visible ray produced by γ Photon Decay in crystal.Four PMT form a four-quadrant receiving area, and each PMT receives γ
The number of visible ray produced by Photon Decay is different, and in each PMT, produced electric impulse signal is also different.
In a decay period, the electric impulse signal that each PMT is changed in γ Photon Decay carries out above-mentioned process
The energy information obtaining afterwards is analyzed, and the energy information obtaining after the electric impulse signal subsequent treatment being exported by PMT is big
Little difference, release energy when can be obtained and being decayed by analysis position the strongest, that is, the crystalline substance receiving γ photon on detector
The particular location of body, and then can learn that positron decay produces the position of γ photon, to obtain radionuclide in biology
Internal concentration distribution.
Need exist for illustrating, the execution sequence of step 306 and step 307 is not defined it is also possible to first carry out step
Rapid 307, in execution step 306, or only carry out one of step.Step 308：Abandon described sampled result.
Embodiment three
Fig. 5 is a kind of present invention signal processing apparatus embodiment three structural representation, is and the method described in embodiment one
Corresponding device, described device includes：
Signal acquiring unit 501, for obtaining the electric impulse signal after amplifying, described electric impulse signal is by light pulse signal
Obtain after carrying out opto-electronic conversion.
Optionally, described signal acquiring unit 501 includes：
First signal acquisition subelement, for obtaining the electric impulse signal after preposition amplification；
Or,
Secondary signal obtains subelement, for obtaining the preposition electric impulse signal amplifying after amplifying with variable gain.
AD conversion unit 502, obtains digital signal for the electric impulse signal after amplifying is carried out analog digital conversion.
Sampling unit 503, for sampling according to prefixed time interval to described digital signal.
Identifying unit 504, for judging whether sampled result is more than predetermined threshold value, if it is, storage described sampling knot
Really.
Time recognition unit 505, for being more than predetermined threshold value using two storing at first in a decay period
Sampled result identify a temporal information of this decay period.
Example IV
Fig. 6 is a kind of present invention signal processing apparatus example IV structural representation, is and the method described in embodiment two
Corresponding device, described device includes：
Signal acquiring unit 501, for obtaining the electric impulse signal after amplifying, described electric impulse signal is by light pulse signal
Obtain after carrying out opto-electronic conversion.
AD conversion unit 502, obtains digital signal for the electric impulse signal after amplifying is carried out analog digital conversion.
Sampling unit 503, for sampling according to prefixed time interval to described digital signal.
Identifying unit 504, for judging whether sampled result is more than predetermined threshold value, if it is, storage described sampling knot
Really.
Time recognition unit 505, for being more than predetermined threshold value using two storing at first in a decay period
Sampled result identify a temporal information of this decay period.
Optionally, described time recognition unit 505 includes：
First computation subunit, for calculate store at first the two double sampling moment being more than predetermined threshold value and when
Between information difference be used for the first ratio；
Second computation subunit, for calculating the two sub- energy values of double sampling being more than predetermined threshold value storing at first
Be used for the second ratio；
3rd computation subunit, for two differences being more than between the double sampling moment of predetermined threshold value storing at first
It is worth for prefixed time interval；
4th computation subunit, for equal with the second ratio using the first ratio, calculates the according to prefixed time interval
The difference of one sampling instant and temporal information is as initial time section；
5th computation subunit, for calculating the sampling instant of first sampled result and demarcating the difference work between the moment
For total period；
6th computation subunit, for calculating the difference between total period and initial time section as this decay period
One temporal information.
Optionally, described time recognition unit 505 includes：
Time identifies subelement, for being more than predetermined threshold value using two storing at first in a decay period
Sampled result, calculates a temporal information identifying this decay period using following formula：
Wherein, m by the sampling instant of first sampled result and demarcates the prefixed time interval being comprised between the moment
Number, E1For the sub- energy value of first sampled result, E2For the sub- energy value of second sampled result, t0For between Preset Time
Every t is temporal information.
Thin interbed unit 601, for identifying this decay using all sampled result being stored in a decay period
One energy information in cycle.
Optionally, described Thin interbed unit 601 includes：
First Thin interbed subelement, for summation of must sampling that all sampled result in a decay period are sued for peace, will
The product of sampling summation and time interval is as the energy information of this decay period；
Or,
Second Thin interbed subelement, for each sampled result of being stored in a decay period and time interval
Product as a sub- energy value, every sub- energy value is sued for peace as an energy information in a decay period.
Embodiment five
Fig. 7 is a kind of present invention signal handling equipment embodiment five structural representation, and described equipment includes：
Photomultiplier tube 701, amplifying circuit 702, modulus conversion chip 703 and programmable logic device 704 are sequentially connected.
Photomultiplier tube 701, for being converted into electric impulse signal by the light pulse signal being received.
Amplifying circuit 702, for being amplified to described electric impulse signal.
Modulus conversion chip 703, obtains digital signal for the electric impulse signal after amplifying is carried out analog digital conversion.
Programmable logic device 704, for sampling according to prefixed time interval to described digital signal；Judge sampling knot
Whether fruit is more than predetermined threshold value, if it is, storing described sampled result；To the sampling knot being stored in a decay period
Fruit carries out time identifying processing and/or Thin interbed is processed.
The foregoing is only presently preferred embodiments of the present invention, not in order to limit the present invention, all essences in the present invention
Within god and principle, any modification, equivalent substitution and improvement made etc., should be included within the scope of the present invention.
Claims (9)
1. a kind of signal processing method is it is characterised in that methods described includes：
Obtain the electric impulse signal after amplifying, described electric impulse signal is carried out obtaining after opto-electronic conversion by light pulse signal；
Electric impulse signal after amplifying is carried out analog digital conversion and obtains digital signal；
Described digital signal is sampled according to prefixed time interval；
Judge whether sampled result is more than predetermined threshold value, if it is, storing described sampled result；
Identify this decay period using the sampled result that two storing at first in a decay period are more than predetermined threshold value
A temporal information；
Described using two stored at first in a decay period be more than predetermined threshold value sampled result identify this decay
One temporal information in cycle includes：
Calculate store at first two be more than the double sampling moment of predetermined threshold value and the difference of temporal information be used for the
One ratio；
Calculate the two sub- energy values of double sampling being more than predetermined threshold value storing at first is used for the second ratio；
The difference that two storing at first are more than between the double sampling moment of predetermined threshold value is prefixed time interval；
Equal with the second ratio using the first ratio, first sampling instant and temporal information are calculated according to prefixed time interval
Difference is as initial time section；
The sampling instant of first sampled result of calculating and the difference between the demarcation moment are as total period；
Calculate the difference between total period and initial time section as a temporal information of this decay period.
2. method according to claim 1 is it is characterised in that methods described also includes：
Identify an energy information of this decay period using all sampled result being stored in a decay period.
3. method according to claim 2 is it is characterised in that described all adopted using store in a decay period
Sample result identifies that an energy information of this decay period includes：
All sampled result in one decay period are sued for peace summation of must sampling, using the product of sampling summation and time interval as
The energy information of this decay period；
Or,
The product of each sampled result being stored in a decay period and time interval, will be every used as a sub- energy value
Individual sub- energy value is sued for peace as an energy information in a decay period.
4. the method according to claim 1-3 any one it is characterised in that described acquisition amplify after electric impulse signal
Including：
Obtain the electric impulse signal after preposition amplification；
Or,
Obtain the preposition electric impulse signal amplifying after amplifying with variable gain.
5. a kind of signal processing apparatus are it is characterised in that described device includes：
Signal acquiring unit, for obtaining the electric impulse signal after amplifying, described electric impulse signal carries out light by light pulse signal
Obtain after electricity conversion；
AD conversion unit, obtains digital signal for the electric impulse signal after amplifying is carried out analog digital conversion；
Sampling unit, for sampling according to prefixed time interval to described digital signal；
Identifying unit, for judging whether sampled result is more than predetermined threshold value, if it is, store described sampled result；
Time recognition unit, for being more than the sampling knot of predetermined threshold value using two storing at first in a decay period
Fruit identifies a temporal information of this decay period；
Described time recognition unit includes：
First computation subunit, for calculating, two storing at first are more than the double sampling moment of predetermined threshold value and the time believes
Breath difference be used for the first ratio；
Second computation subunit, for calculating the ratio of the two sub- energy values of double sampling being more than predetermined threshold value storing at first
As the second ratio；
3rd computation subunit, for store at first two differences being more than between the double sampling moment of predetermined threshold value be
Prefixed time interval；
4th computation subunit, for equal with the second ratio using the first ratio, calculates first according to prefixed time interval
The difference of sampling instant and temporal information is as initial time section；
5th computation subunit, for calculating the difference between the sampling instant of first sampled result and demarcation moment as total
Time period；
6th computation subunit, for calculating the difference between total period and initial time section as one of this decay period
Temporal information.
6. device according to claim 5 is it is characterised in that described device also includes：
Thin interbed unit, for identifying the one of this decay period using all sampled result being stored in a decay period
Individual energy information.
7. device according to claim 6 is it is characterised in that described Thin interbed unit includes：
First Thin interbed subelement, for summation of must sampling that all sampled result in a decay period are sued for peace, will sample
The product of summation and time interval is as the energy information of this decay period；
Or,
Second Thin interbed subelement, for taking advantage of of each sampled result of being stored in a decay period and time interval
Amass as a sub- energy value, every sub- energy value is sued for peace as an energy information in a decay period.
8. the device according to claim 5-7 any one is it is characterised in that described signal acquiring unit includes：
First signal acquisition subelement, for obtaining the electric impulse signal after preposition amplification；
Or,
Secondary signal obtains subelement, for obtaining the preposition electric impulse signal amplifying after amplifying with variable gain.
9. a kind of signal handling equipment is it is characterised in that described equipment includes：
Photomultiplier tube, amplifying circuit, modulus conversion chip and programmable logic device are sequentially connected；
Photomultiplier tube, for being converted into electric impulse signal by the light pulse signal being received；
Amplifying circuit, for being amplified to described electric impulse signal；
Modulus conversion chip, obtains digital signal for the electric impulse signal after amplifying is carried out analog digital conversion；
Programmable logic device, for sampling according to prefixed time interval to described digital signal；Whether judge sampled result
More than predetermined threshold value, if it is, storing described sampled result；The sampled result being stored in a decay period is carried out
Time identifying processing and/or Thin interbed are processed；
One temporal information of one decay period of identification includes：
Calculate store at first two be more than the double sampling moment of predetermined threshold value and the difference of temporal information be used for the
One ratio；
Calculate the two sub- energy values of double sampling being more than predetermined threshold value storing at first is used for the second ratio；
The difference that two storing at first are more than between the double sampling moment of predetermined threshold value is prefixed time interval；
Equal with the second ratio using the first ratio, first sampling instant and temporal information are calculated according to prefixed time interval
Difference is as initial time section；
The sampling instant of first sampled result of calculating and the difference between the demarcation moment are as total period；
Calculate the difference between total period and initial time section as a temporal information of this decay period.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310566812.8A CN103607205B (en) | 2013-11-12 | 2013-11-12 | Signal processing method, apparatus and equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310566812.8A CN103607205B (en) | 2013-11-12 | 2013-11-12 | Signal processing method, apparatus and equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103607205A CN103607205A (en) | 2014-02-26 |
| CN103607205B true CN103607205B (en) | 2017-02-08 |
Family
ID=50125409
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310566812.8A Active CN103607205B (en) | 2013-11-12 | 2013-11-12 | Signal processing method, apparatus and equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103607205B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105490681B (en) * | 2015-11-20 | 2020-01-14 | 东软医疗系统股份有限公司 | Signal processing method and device |
| CN107981877B (en) * | 2016-10-27 | 2021-01-05 | 北京东软医疗设备有限公司 | Sampling method and sampling device of scanning equipment |
| CN109270490A (en) * | 2018-09-30 | 2019-01-25 | 成都精位科技有限公司 | Positioning signal processing method and processing device |
| CN109330598A (en) * | 2018-10-31 | 2019-02-15 | 上海健康医学院 | It is a kind of based on the positron radionuclide real-time counting detector for meeting detection principle |
| CN109884684B (en) * | 2019-01-10 | 2020-09-15 | 中国科学院高能物理研究所 | Electronic data processing method and system for performance test of PET scintillation crystal |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101285715A (en) * | 2008-05-06 | 2008-10-15 | 中国科学院上海光学精密机械研究所 | Weak pulse luminous signs detection apparatus and method |
| CN101452080A (en) * | 2007-12-03 | 2009-06-10 | 通用电气公司 | Integrated neutron-gamma radiation detector with adaptively selected gamma threshold |
| CN101903798A (en) * | 2007-11-02 | 2010-12-01 | 华盛顿大学 | Be used for the data acquisition of positron emission tomography |
| CN102262238A (en) * | 2011-04-19 | 2011-11-30 | 苏州瑞派宁科技有限公司 | Method and device for extracting scintillation pulse information |
-
2013
- 2013-11-12 CN CN201310566812.8A patent/CN103607205B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101903798A (en) * | 2007-11-02 | 2010-12-01 | 华盛顿大学 | Be used for the data acquisition of positron emission tomography |
| CN101452080A (en) * | 2007-12-03 | 2009-06-10 | 通用电气公司 | Integrated neutron-gamma radiation detector with adaptively selected gamma threshold |
| CN101285715A (en) * | 2008-05-06 | 2008-10-15 | 中国科学院上海光学精密机械研究所 | Weak pulse luminous signs detection apparatus and method |
| CN102262238A (en) * | 2011-04-19 | 2011-11-30 | 苏州瑞派宁科技有限公司 | Method and device for extracting scintillation pulse information |
Non-Patent Citations (2)
| Title |
|---|
| "PET系统中数字化精确定时的设计与实现";杜垚垚;《中国硕士学位论文全文数据库信息科技辑》;20130215(第2期);I137-52 * |
| "基于线性光电器件的PET闪烁脉冲数字化方法研究";王希;《中国博士学位论文全文数据库医药卫生科技辑》;20130815(第8期);E076-11 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103607205A (en) | 2014-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103607205B (en) | Signal processing method, apparatus and equipment | |
| CN105212954B (en) | A kind of pulse pile-up event real-time processing method and system | |
| US6936822B2 (en) | Method and apparatus to prevent signal pile-up | |
| CN102262238B (en) | Method and device for extracting scintillation pulse information | |
| EP2867702B1 (en) | Spectral photon counting detector | |
| US6297506B1 (en) | System and method for reducing pile-up errors in multi-crystal gamma ray detector applications | |
| CN102073059B (en) | Digital pileup waveform processing method and system | |
| WO2007043137A1 (en) | Nuclear medical diagnosis device | |
| US10027340B1 (en) | Method and apparatus to digitize pulse shapes from radiation detectors | |
| WO2014121548A1 (en) | Method and system for correcting baseline of digitized scintillation pulse | |
| US20130087710A1 (en) | Trigger Methods in Nuclear Medical Imaging | |
| CN104360376A (en) | Gamma camera having function of identifying radioactive source, namely nuclide, and nuclide identification method | |
| CN105182401B (en) | A kind of nuclear detector crystal positions identification device | |
| WO2021135337A1 (en) | Signal sampling circuit, detection apparatus, and imaging system | |
| CN109581460A (en) | Composite detecting device | |
| CN111751865A (en) | Method for constructing double-time-step PET scanner | |
| CN106842277B (en) | Stacking event processing method and device | |
| CN104932000B (en) | The method and device that γ photons positional information obtains in a kind of PET device | |
| CN107462917A (en) | A kind of ionising radiation detection device and method | |
| KR101174485B1 (en) | Positron Emission Tomogrphy and Signal Process Method of Positron Emission Tomogrphy | |
| US10007011B2 (en) | System for acquisition of tomographic measurement data | |
| CN211698223U (en) | Signal sampling circuit, detection device and imaging system | |
| CN106821409B (en) | Stacking event processing method and device | |
| KR101352771B1 (en) | Coincidence circuit using shift register and PET data aquisition system, radiation counting system, medical diagnosis device using the same | |
| CN108008438A (en) | The measuring device and method of a kind of ray energy |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Patentee after: Shenyang Neusoft Medical Systems Co.,Ltd. Address before: Hunnan New Century Road 110179 Shenyang city of Liaoning Province, No. 16 Patentee before: SHENYANG NEUSOFT MEDICAL SYSTEMS Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20230512 Address after: Room 308, No. 177-2 Chuangxin Road, Hunnan District, Shenyang City, Liaoning Province, 110167 Patentee after: Shenyang Zhihe Medical Technology Co.,Ltd. Address before: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Patentee before: Shenyang Neusoft Medical Systems Co.,Ltd. |