CN108287361B - A kind of detection method and device in single event dead time - Google Patents
A kind of detection method and device in single event dead time Download PDFInfo
- Publication number
- CN108287361B CN108287361B CN201810005487.0A CN201810005487A CN108287361B CN 108287361 B CN108287361 B CN 108287361B CN 201810005487 A CN201810005487 A CN 201810005487A CN 108287361 B CN108287361 B CN 108287361B
- Authority
- CN
- China
- Prior art keywords
- single event
- data
- dead time
- decay
- counting rate
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 17
- 230000002285 radioactive effect Effects 0.000 claims abstract description 56
- 239000013078 crystal Substances 0.000 claims abstract description 47
- 238000000034 method Methods 0.000 claims abstract description 21
- 206010033799 Paralysis Diseases 0.000 claims description 16
- 238000010408 sweeping Methods 0.000 claims description 8
- 230000005253 gamme decay Effects 0.000 description 18
- 230000005255 beta decay Effects 0.000 description 9
- 238000001228 spectrum Methods 0.000 description 8
- 230000005855 radiation Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000005264 electron capture Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
- G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
The embodiment of the invention discloses a kind of detection method and device in single event dead time, it include the PET system of intrinsic radionuclide suitable for scintillation crystal, first obtain the first data after sky is swept, the second data after placing radioactive source scanning are obtained again, data after sky is swept are the corresponding data of intrinsic radionuclide decay, and the energy area of Radioactive Source Decay is different from the intrinsic radionuclide area the first data Zhong Neng, in this way, the energy area for counting of intrinsic radionuclide can be distinguished in the data after radioactive source scanning, to, the corresponding counting rate of intrinsic radionuclide decay after placing radioactive source scanning can be obtained, and the photon count rate that identical decay generates after the twice sweep of front and back is identical, pass through first data and second data, the single event dead time can be obtained.This method times of collection is few and required time is short, and easy to operate, improves the utilization efficiency of PET system.
Description
Technical field
The present invention relates to medical image imaging and processing technology field more particularly to a kind of detection sides in single event dead time
Method and device.
Background technique
Positron emission tomography (Positron Emission Computed Tomography, abbreviation PET) system
System is one of current state-of-the-art medical diagnosis imaging device, using the nucleic that β decay occurs, generates just in decay
Electronics is buried in oblivion effect and is generated γ photon by the electron capture of surrounding, counts to these photons, and then be implemented as
Picture.
It is made of in PET system multiple subsystems, one minimum of each subsystem demands when handling continuous two events
The time is handled, which is referred to as the dead time.The PET system dead time is the minimum treat that detector processing module requires
Time, within the minimum treat time, it can be ensured that the processing for completing an annihilation event, in this way, the latter is avoided to bury in oblivion thing
Part has little time the previous annihilation event that is disposed before occurring, and leads to the loss of previous annihilation event, causes the dead time
Loss, the loss in dead time will cause photon counting loss, and then underestimate the pixel value of image, influence image quality.Currently, main
If determining the PET system dead time using decayed source method, this method needs the radioactive source of high activity to be measured for a long time, surveys
It tries difficulty height and the time is long.The PET system dead time includes the single event dead time and meets the event dead time, and PET system is dead
Time is mainly contributed by the single event dead time.
Summary of the invention
It is easily operated and high-efficient the present invention provides a kind of detection method and device in single event dead time.
The present invention provides a kind of detection methods in single event dead time, are applied to PET system, the spy of the PET system
The scintillation crystal surveyed in device includes intrinsic radionuclide, comprising:
Execute sky sweep, obtain sky sweep after the first data, first data include the first score single event counting rate and
The first single event counting rate in energy area, the energy area correspond to the full energy peak region of the intrinsic radionuclide decay；
The scanning for placing radioactive source is executed, the second data are obtained, second data are counted including the second score single event
Second single event counting rate of rate and the energy area, wherein the full energy peak region of the Radioactive Source Decay is different from the energy area
Corresponding full energy peak region；
By first data and second data, the single event dead time is obtained.
Optionally, the energy area is one or more.
Optionally, the intrinsic radionuclide of the scintillation crystal includes Lu176 element.
Optionally, the energy area is one, described by first data and second data, and it is dead to obtain single event
Time, comprising:
According to the corresponding true single event counting rate of decay same after the non-paralysis formula and front and back twice sweep in dead time
It is identical, by first data and second data, obtain the single event dead time.
Optionally, the specific energy area is multiple, described to pass through first data and second data, obtains single thing
The part dead time, comprising:
According to respectively decay after the non-paralysis formula and front and back twice sweep in dead time corresponding true single event counting rate it
With it is constant, pass through first data and the second data, obtain the single event dead time.
A kind of detection device in single event dead time is applied to PET system, the flashing in the detector of the PET system
Crystal includes intrinsic radionuclide, comprising:
Scanning element after executing sky and sweeping, obtains the first data, institute for executing the empty scanning for sweeping and placing radioactive source
The first single event counting rate that the first data include the first score single event counting rate and energy area is stated, the energy area corresponds to described
The full energy peak region of intrinsic radionuclide decay；After executing the scanning for placing radioactive source, the second data of acquisition, described second
Data include the second single event counting rate of the second score single event counting rate and the energy area, wherein the Radioactive Source Decay
Full energy peak region be different from it is described can the corresponding full energy peak region in area；
Single event dead time determination unit, for it is dead to obtain single event by first data and second data
Time.
Optionally, the energy area is one or more.
Optionally, the intrinsic radionuclide of the scintillation crystal includes Lu176 element.
Optionally, the energy area is one, according to decay pair same after the non-paralysis formula and twice sweep in dead time
The true single event counting rate answered is identical, by first data and second data, obtains the single event dead time.
Optionally, the energy area is multiple, public according to the non-paralysis in dead time in the single event dead time determination unit
The sum of the corresponding true single event counting rate that respectively decays after formula and twice sweep is constant, passes through first data and the second number
According to the acquisition single event dead time.
The detection method and device in single event dead time provided in an embodiment of the present invention, it includes intrinsic for being suitable for scintillation crystal
The PET system of radionuclide first obtains the first data after sky is swept, then obtains the second data after placing radioactive source scanning,
Data after sky is swept are the corresponding data of intrinsic radionuclide decay, and the energy area of Radioactive Source Decay and intrinsic radioactivity
The area nucleic the first data Zhong Neng is different, in this way, can distinguish intrinsic radionuclide in the data after radioactive source scanning
For counting energy area, it is thus possible to obtain place radioactive source scanning after the corresponding counting of intrinsic radionuclide decay
Rate, and the photon count rate that identical decay generates after the twice sweep of front and back is identical, passes through first data and second number
According to the single event dead time can be obtained.This method can by reasonably selecting radioactive source, by the enumeration data of twice sweep,
It can determine the single event dead time, times of collection is few and required time is short, and easy to operate, and that improves PET system utilizes effect
Rate.
Detailed description of the invention
In order to more clearly explain the technical solutions in the embodiments of the present application, make required in being described below to embodiment
Attached drawing is briefly described, it should be apparent that, the drawings in the following description are only some examples of the present application, for
For those of ordinary skill in the art, without any creative labor, it can also be obtained according to these attached drawings
His attached drawing.
Fig. 1 is the Energy distribution schematic diagram of the discretization of LYSO crystal；
Fig. 2 is the flow chart according to the detection method in the single event dead time of the embodiment of the present invention；
Fig. 3 is the structural schematic diagram according to the determining device in the single event dead time of the embodiment of the present invention.
Specific embodiment
Below in conjunction with the attached drawing in the embodiment of the present application, technical solutions in the embodiments of the present application carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of embodiments of the present application, instead of all the embodiments.It is based on
Embodiment in the application, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall in the protection scope of this application.
The present invention provides a kind of detection methods in single event dead time, are applied to PET system.In order to better understand originally
The technical solution of invention is below first described the scintillation crystal of PET system and its characteristic.The single event dead time of the application
The PET system that is applicable in of detection method in, the scintillation crystal of detector includes intrinsic radionuclide, that is to say, that flashing
Crystal will naturally decay, and the decay period of the intrinsic radionuclide of scintillation crystal is long and stablizes, it is therefore contemplated that
The intrinsic radionuclide decay and the photon counting generated is invariable.In this way, without external radioactive source and having
In the case where external source, by the counting rate for scintillation crystal itself decay that detector detects, and then consolidated according to scintillation crystal
There is radionuclide decay to count constant characteristic, it can be by twice sweep data, to obtain the single event dead time.
The intrinsic radionuclide decay of scintillation crystal generally includes a variety of different decays, the energy of decay or ray
Type can be different, and for different decays, different peak regions is showed in the power spectrum of Energy distribution, namely
Full energy peak region, different peak regions correspond to different decays, after power spectrum is carried out discretization, different decays pair
Should be in different energy areas, energy area is the energy section that continuous road location is constituted in entire score.
In a specific embodiment of the present invention, scintillation crystal can be the scintillation crystal comprising Lu176, that is, flashing crystalline substance
The intrinsic radionuclide of body includes Lu176, such as can be LSO or LYSO crystal or LFS crystal etc., is existed in these crystal
Natural radionuclide can introduce the Lu176 element of synthesized element in crystal growing process, and Lu176 element can occur to decline certainly
Become, half-life period is extremely long, and related decay energy and activity stabilization are outstanding.LSO, LFS or LYSO crystal, radioactivity show as β
It decays and with corresponding gamma decay, the energy that β decay occurs is 589keV, is respectively with the energy that gamma decay occurs
307keV、202keV、88keV。
Refering to what is shown in Fig. 1, the Energy distribution schematic diagram of the discretization for LYSO crystal, wherein horizontal axis is road location, and the longitudinal axis is
Counting rate, it is each can area correspond to several segments in entire score, the Lu176 element in LYSO crystal occur β decay with
And with occur gamma decay, β decay energy be 589keV, with occur gamma decay energy be respectively 307keV, 202keV,
88keV, it can be seen that, the road location range entirely to decay is { 0,1500 }, that is, the road of score from Energy distribution schematic diagram
Location range, four energy peak dislocation positions arrangement of decay, in specific application, the region where can choose full energy peak is institute
The energy area to be counted, for example, can area A, B, C, D, the location Qi Dao section be respectively { 200,300 }, { 600,800 }, { 1000,1100 },
{ 1150,1300 }, these can areas correspond respectively to the gamma decay of 88keV, the gamma decay of 202keV, 307keV gamma decay and
The β of 589Kev decays.It is understood that the range in energy area herein is merely illustrative, energy area used by the embodiment of the present invention
Road location range be not limited to this.
It include the PET system of intrinsic radionuclide based on the scintillation crystal in detector, below with reference to specific reality
Example is applied the detection method in single event dead time is described in detail.
Embodiment one
Refering to what is shown in Fig. 2, executing sky in step S101 and sweeping, the first data after sky is swept are obtained, the first data include first
First single event counting rate of score single event counting rate and energy area, the energy area correspond to the intrinsic radionuclide decay
Full energy peak region.
First data are that PET system carries out the enumeration data obtained after single pass, this time scanning is swept for sky, and sky, which is swept, is
Refer to and do not place any radiation source in the visual field of pet detector and the scanning that carries out, in this way, scanning it for the first time at this
Afterwards, the number of photons whole scintillation crystal itself generated in detector decays and generates.
After sky is swept, PET system acquisition counter data can obtain the counting of the first score single event from enumeration data
The first single event counting rate in the area Shuai Heneng, single event counting rate, that is, photon counting of collected unit time, score single event
Counting rate is the unit time photon counting within the scope of the road location of entire power spectrum, and the single event counting rate in energy area is power spectrum intermediate energy region
Road location within the scope of unit time photon counting.
In embodiments of the present invention, the energy area counted can for one or more, that is, can obtain one or
The single event counting rate in multiple energy areas.In some embodiments, the scintillation crystal of detector is LYSO crystal, specifically, being used for
The energy area of counting can be the γ corresponding to the β decay of 589keV, the gamma decay of 307keV, the gamma decay of 202keV or 88keV
One or more of full energy peak region of decay then in this step counts the single event for obtaining one or more energy areas
Rate.
In the present embodiment, it is illustrated with the first single event counting rate for obtaining an area Ge Neng.Specifically show at one
In example, the scintillation crystal of detector is LYSO crystal, in the step, is obtained after sky sweeps, the first score list of the crystal decay
Event count rate ma1And the first single event counting rate m in the corresponding energy area in full energy peak region of the β decay of 589Kevp1。
In step S102, the scanning for placing radioactive source is executed, obtains the second data, the second data include the second score list thing
Second single event counting rate of part counting rate and the energy area, wherein the full energy peak region of the Radioactive Source Decay is different from institute
State the corresponding full energy peak region in specific energy area.
Second data are to place radioactive source to carry out scanning the enumeration data obtained later again, in this time scanning, PET detection
It placed external radioactive source in the visual field of device.
After placing radioactive source scanning, external radioactive source is contained in the enumeration data of acquisition and scintillation crystal is inherently put
The counting that the decay of penetrating property nucleic generates needs to select suitable external radiation according to the energy area for counting in previous step
Source, so that external radioactive source will not generate interference to the counting in the specific energy area of intrinsic radionuclide decay.
In embodiments of the present invention, full energy peak region and above-mentioned intrinsic radioactive nucleus in the power spectrum of the external Radioactive Source Decay
The first area data Zhong Neng obtained in element decay is different, that is, any full energy peak in the power spectrum of external Radioactive Source Decay is not
There is coincidence in full energy peak region corresponding with above-mentioned energy area, in this way, making external radioactive source will not be to intrinsic radionuclide decay
Energy area counting generate interference.
Score single event counting rate in second data is the unit time photon counting within the scope of the road location of entire power spectrum,
The power spectrum contains the decay of external radioactive source and the intrinsic radionuclide decay of scintillation crystal, also includes in the second data
The single event counting rate in the area You Neng, the energy area are identical energy area in previous step, that is, the intrinsic radioactivity of scintillation crystal
The energy area of the counting of nucleic decay.
In a specific example, the scintillation crystal of detector is LYSO crystal, and the first area data Zhong Neng is
The full energy peak region of the β decay of 589Kev select external radioactive source of the full energy peak not for the β decay of 589Kev, example in the step
Such as can be54Mn (835keV) in the step, is obtained after placing external radioactive source scanning, the external radioactive source and crystal
Second score single event counting rate m of decaya2And the second of the corresponding energy area in full energy peak region of the β decay of 589Kev is single
Event count rate mp2。
It, can be with it should be noted that the application is to the execution of above-mentioned steps S101 and step S102 sequence and without limitation
It executes step S101 and executes step S102 later, execute step S101 after step S102 can also be executed.
In step S103, the photon count rate generated according to decay identical after the twice sweep of front and back is identical, passes through described the
One data and second data obtain the single event dead time.
Counting rate in above-mentioned the first data and the second data is all the photon counting of collected unit time, also
It is the counting rate that detector detects, the single event counting loss that detector detects is the photon counting generated with practical decay
Rate is directly proportional, and the photon count rate which generates usually is denoted as true single event counting rate, is placing radioactive source scanning
The photon count rate of generation of decaying afterwards is higher than after sky is swept the photon count rate of generation of decaying, and places single thing after radioactive source scanning
Part counting loss is equally higher than the single event counting loss after empty sweep, also identical for the counting in energy area, and due to scintillation crystal
The stability of intrinsic radionuclide decay, the identical corresponding true single event that decays of intrinsic radionuclide after this twice sweep
Counting rate should remain unchanged.
In the present embodiment, the photon count rate generated according to decay same after the twice sweep of front and back is identical, specific real
In existing, the photon count rate for the generation that decays can be embodied by the corresponding true single event counting rate of the decay, be obtained by above-mentioned
The first data and the second data obtained, available single event dead time.In specific application, the non-of dead time can be used
Paralysis formula is calculated.In other application, when counting rate is lower, the close of the non-paralysis formula in dead time also can use
It is calculated like formula.
Specifically, sky sweeps the corresponding true single event counting rate n in rear a certain energy area according to the non-paralysis formula in dead timep1Are as follows:
Wherein, ma1The first score single event counting rate after being swept for sky, mp1For first single event in the area Kong Sao Hou Gaineng
Counting rate, τ are the single event dead time.
Place the corresponding true single event counting rate n in the radiation source scanning area Hou Nengp2Are as follows:
Wherein, ma2To place the second score single event counting rate after radioactive source scanning, mp2After placing radioactive source scanning
The second single event counting rate in the energy area, τ are the single event dead time.
It is recognised that np1=np2, by formula (1) and formula (2) and the first data and the second data, then can obtain
Obtain single event dead time τ are as follows:
In this way, just obtaining the single event dead time in PET system, it includes intrinsic radiation that this method, which is suitable for scintillation crystal,
Property nucleic PET system, first obtain sky sweep after the first data, then obtain place radioactive source scanning after the second data, in sky
Data after sweeping are the corresponding data of intrinsic radionuclide decay, and the energy area of Radioactive Source Decay and intrinsic radionuclide
The energy Qu Butong for counting in first data, in this way, intrinsic radiation can be distinguished in the data after radioactive source scanning
Property nucleic for counting energy area, it is thus possible to obtain place radioactive source scanning after the specific decay of intrinsic radionuclide
Corresponding counting rate, and the photon count rate that same decay generates after the twice sweep of front and back is identical, by first data and
Second data, can obtain the single event dead time.Pass through the enumeration data of twice sweep, it can when determining that single event is dead
Between, times of collection is few and required time is short, and easy to operate, improves the utilization efficiency of PET system.
Embodiment two
What is different from the first embodiment is that in the present embodiment for counting energy area be it is multiple, pass through it is multiple can areas count
Rate determines single event dead time, the part different from embodiment one described in detail below, and same section will not be described in great detail.
In step S201, the first data after sky is swept are obtained, the first data include the first score single event counting rate and more
The first single event counting rate in the area Ge Neng, the energy area correspond to the full energy peak region of the intrinsic radionuclide decay.
In the present embodiment, it is possible to which area is multiple, that is, obtain the single event counting rate in multiple energy areas.It is specific at one
Example in, the scintillation crystal of detector is LYSO crystal, in the step, is obtained after sky sweeps, obtaining the first data includes: this
First score single event counting rate m of crystal decaya1And the full energy peak region corresponding first energy area of the gamma decay of 307keV
First single event counting rate m1p1, the first single event counting rate in the full energy peak region of the gamma decay of 202keV corresponding second energy area
m2p1And the full energy peak region of the gamma decay of 88keV corresponds to the first single event counting rate m3 in third energy areap1。
In step S202, the second data after placing radioactive source scanning are obtained, the second data include the second score single event
Second single event counting rate of counting rate and multiple energy areas, wherein the full energy peak region of the Radioactive Source Decay is different from
It is described can the corresponding full energy peak region in area.
In the present embodiment, in the full energy peak region of Radioactive Source Decay and above-mentioned first data all energy areas for counting
All different, external radioactive source will not generate interference to the counting in energy all these in the first data area.
Still by taking above-mentioned specific example as an example, full energy peak is selected in the step not is the gamma decay of 88keV, 202keV
Gamma decay, 307keV gamma decay external radioactive source, such as can be137Cs (662keV) or54Mn (835keV), the step
In, it obtains after placing the scanning of external radioactive source, the second data of acquisition include: the of the external radioactive source and crystal decay
Two score single event counting rate ma2And the second single event meter in the full energy peak region corresponding first energy area of the gamma decay of 307keV
Digit rate m1p2, the second single event counting rate m2 in the full energy peak region of the gamma decay of 202keV corresponding second energy areap2, and
The full energy peak region of the gamma decay of 88keV corresponds to the second single event counting rate m3 in third energy areap2。
In step S203, the photon count rate generated according to decay identical after twice sweep is identical, passes through first number
According to second data, obtain the single event dead time.
In the embodiment, the photon count rate according to the generation respectively to decay after the twice sweep of front and back is identical, then each decay pair
The sum of the true single event counting rate answered is constant, and by the first data and the second data of above-mentioned acquisition, available single event is dead
Time.In specific application, it can be calculated using the non-paralysis formula in dead time.In other application, work as counting rate
When lower, the approximate formula that also can use the non-paralysis formula in dead time is calculated.
By taking above-mentioned specific example as an example, according to the non-paralysis formula in dead time, sky sweeps the corresponding very single thing of rear each energy region
Part counting rate is respectively as follows:
Wherein, ma1The first score single event counting rate after being swept for sky, m1p1、m2p1、m3p1Respectively it is empty sweep after first
The first single event counting rate in energy area, the second energy area, third energy area, τ are the single event dead time.
The corresponding true single event counting rate of each energy region is respectively as follows: after placing radiation source scanning
Wherein, ma2To place the second score single event counting rate after radioactive source scanning, m1p2、m2p2、m3p2Respectively put
The second single event counting rate in the first energy area, the second energy area, third energy area after setting radioactive source scanning, τ is the single event dead time.
It is recognised that n1p1+n2p1+n3p1=n1p2+n2p2+n3p2, pass through above-mentioned formula (4)-(9), and and first
Data and the second data, then can obtain single event dead time τ.
In this embodiment, the counting rate for obtaining the multiple energy areas of scintillation crystal passes through the sum of multiple energy areas counting rate and protects
It holds constant, calculate and obtain the single event dead time, can increase counting statistics, shorten acquisition time, reduce statistical error.
The detection method in the single event dead time of the embodiment of the present invention is described in detail above, in addition, this hair
It is bright to additionally provide the detection device for realizing the single event dead time of the above method, it is applied to PET system, the spy of the PET system
The scintillation crystal surveyed in device includes intrinsic radionuclide, refering to what is shown in Fig. 3, including:
Scanning element 300 after executing sky and sweeping, obtains the first data for executing the empty scanning for sweeping and placing radioactive source,
First data include the first single event counting rate of the first score single event counting rate and energy area, and the energy area corresponds to institute
State the full energy peak region of intrinsic radionuclide decay；After executing the scanning for placing radioactive source, the second data are obtained, described the
Two data include the second single event counting rate of the second score single event counting rate and the energy area, wherein the radioactive source declines
The full energy peak region of change is different from the corresponding full energy peak region in the energy area；
Single event dead time determination unit 320, the photon count rate phase for being generated according to decay identical after twice sweep
Together, by first data and second data, the single event dead time is obtained.
Further, the energy area is one or more.
Further, the intrinsic radionuclide of the scintillation crystal includes Lu176 element.
Further, it is described can area be one, in the single event dead time determination unit 320, it is non-according to the dead time
The same corresponding true single event counting rate that decays is identical after paralysis formula and twice sweep, passes through first data and institute
The second data are stated, the single event dead time is obtained.
Further, it is described can area be it is multiple, it is non-according to the dead time in the single event dead time determination unit 320
The sum of the corresponding true single event counting rate that respectively decays is constant after paralysis formula and twice sweep, by first data and
Second data obtain the single event dead time.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality
For applying example, since it is substantially similar to the method embodiment, so describing fairly simple, related place is referring to embodiment of the method
Part explanation.System embodiment described above is only schematical, wherein described be used as separate part description
Module or unit may or may not be physically separated, the component shown as module or unit can be or
Person may not be physical unit, it can and it is in one place, or may be distributed over multiple network units.It can root
According to actual need that some or all of the modules therein is selected to achieve the purpose of the solution of this embodiment.Ordinary skill
Personnel can understand and implement without creative efforts.
It should be noted that, in this document, relational terms such as first and second and the like are used merely to a reality
Body or operation are distinguished with another entity or operation, are deposited without necessarily requiring or implying between these entities or operation
In any actual relationship or order or sequence.Moreover, the terms "include", "comprise" or its any other variant are intended to
Non-exclusive inclusion, so that the process, method, article or equipment including a series of elements is not only wanted including those
Element, but also including other elements that are not explicitly listed, or further include for this process, method, article or equipment
Intrinsic element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that
There is also other identical elements in process, method, article or equipment including the element.
Claims (10)
1. a kind of detection method in single event dead time is applied to PET system, the flashing in the detector of the PET system is brilliant
Body includes intrinsic radionuclide characterized by comprising
It executes sky to sweep, obtains the first data after sky is swept, first data include the first score single event counting rate and can area
The first single event counting rate, it is described can area correspond to the intrinsic radionuclide decay full energy peak region；
Execute place radioactive source scanning, obtain the second data, second data include the second score single event counting rate and
The second single event counting rate in the energy area, wherein it is corresponding that the full energy peak region of the Radioactive Source Decay is different from the energy area
Full energy peak region；
According to respectively decaying after the non-paralysis formula and twice sweep in dead time, the sum of corresponding true single event counting rate is constant,
By first data and the second data, the single event dead time is obtained.
2. the method according to claim 1, wherein the energy area is one.
3. the method according to claim 1, wherein the intrinsic radionuclide of the scintillation crystal includes
Lu176 element.
4. the method according to claim 1, wherein it is described can area be one, it is described to pass through first data
With second data, the single event dead time is obtained, comprising:
It is identical according to the corresponding true single event counting rate of decay same after the non-paralysis formula and twice sweep in dead time, lead to
First data and second data are crossed, the single event dead time is obtained.
5. the method according to claim 1, wherein the energy area is multiple.
6. a kind of detection device in single event dead time is applied to PET system, the flashing in the detector of the PET system is brilliant
Body includes intrinsic radionuclide characterized by comprising
Scanning element, for executing the empty scanning for sweeping and placing radioactive source, after executing sky and sweeping, the first data of acquisition, described the
One data include the first single event counting rate of the first score single event counting rate and energy area, and the energy area corresponds to described intrinsic
The full energy peak region of radionuclide decay；After executing the scanning for placing radioactive source, the second data, second data are obtained
Including the second score single event counting rate and it is described can area the second single event counting rate, wherein the Radioactive Source Decay it is complete
Energy peak region is different from the corresponding full energy peak region in the energy area；
Single event dead time determination unit, for respectively decay to correspond to after the non-paralysis formula and twice sweep according to the dead time
The sum of true single event counting rate it is constant, pass through first data and the second data, obtain the single event dead time.
7. device according to claim 6, which is characterized in that the energy area is one.
8. device according to claim 6, which is characterized in that the intrinsic radionuclide of the scintillation crystal includes
Lu176 element.
9. device according to claim 6, which is characterized in that the energy area is one, and the single event dead time determines
It is identical according to the corresponding true single event counting rate of decay same after the non-paralysis formula and twice sweep in dead time in unit,
By first data and second data, the single event dead time is obtained.
10. device according to claim 6, which is characterized in that the energy area is multiple.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810005487.0A CN108287361B (en) | 2018-01-03 | 2018-01-03 | A kind of detection method and device in single event dead time |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810005487.0A CN108287361B (en) | 2018-01-03 | 2018-01-03 | A kind of detection method and device in single event dead time |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108287361A CN108287361A (en) | 2018-07-17 |
| CN108287361B true CN108287361B (en) | 2019-08-27 |
Family
ID=62834929
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810005487.0A Active CN108287361B (en) | 2018-01-03 | 2018-01-03 | A kind of detection method and device in single event dead time |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108287361B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3845928A1 (en) * | 2019-12-30 | 2021-07-07 | Canon Medical Systems Corporation | Pet apparatus and method of performing timing calibration |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110327065B (en) * | 2019-07-18 | 2023-01-03 | 上海联影医疗科技股份有限公司 | Single-event counting rate correction method and device, computer equipment and storage medium |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100508891C (en) * | 2007-09-04 | 2009-07-08 | 陈武凡 | Maximum posteriori optimizing image rebuilding method in PET imaging |
| CN101856236B (en) * | 2010-06-13 | 2012-07-04 | 苏州瑞派宁科技有限公司 | Positron emission tomography (PET) method and device with application adaptability |
| CN104224219B (en) * | 2014-08-11 | 2017-05-03 | 沈阳东软医疗系统有限公司 | Conforming pair judging method and equipment |
| US10101474B2 (en) * | 2014-12-15 | 2018-10-16 | Koninklijke Philips N.V. | Pixel based dead time correction |
| CN106251313B (en) * | 2016-08-15 | 2020-06-26 | 上海联影医疗科技有限公司 | Medical imaging method and system |
| CN106296768B (en) * | 2016-08-16 | 2019-02-19 | 赛诺联合医疗科技（北京）有限公司 | A kind of method and apparatus of Data correction |
| CN106353786B (en) * | 2016-09-30 | 2019-06-28 | 上海联影医疗科技有限公司 | Positron emission tomography time-of-flight capability detection method and device |
| CN106344060B (en) * | 2016-09-30 | 2019-06-28 | 上海联影医疗科技有限公司 | Coincidence correction method and dead time detection method for pet detector |
-
2018
- 2018-01-03 CN CN201810005487.0A patent/CN108287361B/en active Active
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3845928A1 (en) * | 2019-12-30 | 2021-07-07 | Canon Medical Systems Corporation | Pet apparatus and method of performing timing calibration |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108287361A (en) | 2018-07-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sarrut et al. | Advanced Monte Carlo simulations of emission tomography imaging systems with GATE | |
| US7238946B2 (en) | Nuclear imaging system using scintillation bar detectors and method for event position calculation using the same | |
| CN104412127B (en) | Digital positron emission tomography（DPET）Energy calibration system and method | |
| EP2047297B1 (en) | Method and system for improved tof pet reconstruction | |
| Gonias et al. | Validation of a GATE model for the simulation of the Siemens biograph™ 6 PET scanner | |
| CN102007430B (en) | Radiation tomographic equipment | |
| US20180114345A1 (en) | Correcting count loss | |
| US7470897B2 (en) | Nuclear medical apparatus | |
| CN107250844A (en) | Method and apparatus for handling PET data | |
| CN105637386B (en) | Histogram in positron emission tomography (PET) energy histogram is smooth | |
| CN108287361B (en) | A kind of detection method and device in single event dead time | |
| Aklan et al. | GATE Monte Carlo simulations for variations of an integrated PET/MR hybrid imaging system based on the Biograph mMR model | |
| Wang et al. | Development and characterization of modular readout design for two-panel head-and-neck dedicated PET system based on CZT detectors | |
| RU2705758C2 (en) | Low-cost design of digital pet (positron emission tomograph) | |
| JP5736566B2 (en) | Radioactivity absolute measurement method by liquid scintillation and calibration method of radioactivity measurement device | |
| CN107462917A (en) | A kind of ionising radiation detection device and method | |
| US9140804B2 (en) | Methods and systems for determining timing recovery information in a positron emission tomography (PET) system | |
| CN107440734A (en) | A kind of PET scan method and apparatus | |
| Ahmed et al. | A validated Geant4 model of a whole-body PET scanner with four-layer DOI detectors | |
| US8957384B2 (en) | Gamma ray detector linearity calibration | |
| Efthimiou et al. | Data-driven, energy-based method for estimation of scattered events in positron emission tomography | |
| DE112015002766B4 (en) | Calibration in single-photon emission computed tomography with multi-emission energies | |
| Brown | In-silico optimisation of tileable philips digital SiPM based thin monolithic scintillator detectors for SPECT applications | |
| US10517559B2 (en) | Nuclear medicine diagnostic apparatus, diagnostic imaging apparatus, and image processing method | |
| US5347455A (en) | Scintillation camera |
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 | ||
| CB02 | Change of applicant information |
Address after: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Applicant after: Shenyang Neusoft Medical Systems Co.,Ltd. Address before: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Applicant before: SHENYANG NEUSOFT MEDICAL SYSTEMS Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20230411 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. |
|
| TR01 | Transfer of patent right |