CN106539591A - PET flight time state quality detection methods and PET scan device - Google Patents
PET flight time state quality detection methods and PET scan device Download PDFInfo
- Publication number
- CN106539591A CN106539591A CN201510603207.2A CN201510603207A CN106539591A CN 106539591 A CN106539591 A CN 106539591A CN 201510603207 A CN201510603207 A CN 201510603207A CN 106539591 A CN106539591 A CN 106539591A
- Authority
- CN
- China
- Prior art keywords
- flight
- time
- flight time
- pet
- response
- 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.)
- Granted
Links
Landscapes
- Nuclear Medicine (AREA)
Abstract
The invention discloses a kind of PET flight time state quality detection method, including：Die body is placed in into PET scan within sweep of the eye, the die body axis and scan vision scope diameter parallel；Detection meets event, obtains meeting data, and meets data acquisition motif position according to described；According to motif position, the differential time of flight determined by position for meeting in line of response that two photon of event reaches a pair of detector crystals is calculated；According to the flight time sinogram for meeting data and time bias table to be detected, the flight time aberration measurements of line of response are obtained；By comparing, whether the differential time of flight that determined by position is equal with flight time aberration measurements to judge whether flight time state quality is qualified.The die body of the inventive method accurately detects flight time state quality by need not being placed in the center for scanning cavity, simple to operate.Additionally, the present invention also provides a kind of PET scan device.
Description
【Technical field】
The present invention relates to the detection of medical imaging device quality state, more particularly to a kind of PET flight time shape
State quality determining method and PET scan device.
【Background technology】
Positron emission computerized tomography technology (Positron Emission Tomography, PET) conduct is worked as
Modern high-level nuclear medicine technology, and the state-of-the-art larger medical diagnosing image that Medical circle is generally acknowledged sets
It is one of standby, it has also become tumor, the indispensable important method of cardio-cerebral diseases diagnosis.Its operation principle is main
For：The radiopharmaceutical that positron can be launched is injected into biological internal, the positron electricity internal with biology
Son is buried in oblivion, according to law of conservation of energy and the law of conservation of momentum, produce there is 511KeV, to contrary
A pair of gammaphotons that direction flies out；The detector array being looped around around measurand is transported to these opposite directions
Dynamic photon carry out coincidence measurement (ray that two detector crystals are produced to two gammaphotons while detection,
Also referred to as meet event), and produce coincidence measurement line of response (Line of Response, LOR)；Using calculating
Machine process these data for projection carry out image reconstruction can solve radionuclide in measurand body point
Cloth；The multiple tomographies for reconstructing sequentially are stacked to form into 3-D view.
Common PET is required for injecting radioactive indicator before checking, in above process, measures two
Time difference of the photon to detector, according to known diameter detector and the light velocity, it is possible to determine that photon occurs
Position (the transmitting position of positron), so as to obtain the position of tracer decays, above-mentioned technology is referred to as
Flight time (Time of Flight, TOF).TOF technologies have the advantage that：(1) it can be found that very
To less than 2 millimeters of tumor focus, improve picture quality, substantially increase minimal disease Detection capability and really
Rate is examined, the tumor focus of patient's investigation early stage is helped and for minimal neoplastic metastasis Accurate Diagnosis；(2)
The conventional general physical checkup time can be completed in more than ten minutes, and patient radiation's property drug dose also substantially subtracts
It is few, radiation dose is reduced to greatest extent, the time required to substantially reducing scanning；(3) hypersensitivity,
Can complete the preoperative positioning to epileptic, apply also for presenile dementia diagnosis, mental sickness and examine
It is disconnected.
Existing PET scan TOF quality testings are usually to place radioactive source in PET scan chamber central,
The annihilation photon of radioactive source is collected, detector crystal obtains data and calculated, so as to provide PET scan
The time state of instrument.However, this method common problem is radioactive source need to be positioned at PET scan
Chamber central, cannot be processed during placement eccentric for radioactive source, be positioned at scanning chamber central for behaviour
Need to take a substantial amount of time for author.Based on this, it is necessary to which existing TOF quality determining methods are entered
Row is improved.
【The content of the invention】
The technical problem to be solved is to propose a kind of time of flight mass detection side simple to operate
Method.
The present invention solve the technical scheme that adopted of above-mentioned technical problem for：PET flight time state quality is examined
Survey method, including：
Die body is placed in the range of the scan vision of PET scan cavity, the die body axis and scan vision
Scope diameter parallel；
Detection meets event, obtains meeting data, and meets data acquisition motif position according to described；
According to motif position, calculate meet in line of response two photon of event reach a pair of detector crystals by position
The differential time of flight of decision is put, the line of response is the flight path for meeting two photon of event；
According to the flight time sinogram for meeting data and time bias table to be detected, the flight of line of response is obtained
Time difference measurements value；
Flown by comparing the differential time of flight judgement whether equal with flight time aberration measurements determined by position
Whether row time state quality is qualified.
Further, the die body is linear radiation source, uniform tubbiness radioactive source, uniform hollow tubbiness radiation
Any one in source.
Further, according to motif position, calculate and in line of response, meet two photon of event, a pair of detectors of arrival
The differential time of flight detailed process determined by position of crystal is：
The range difference of two detector distance motif positions in the line of response, root are calculated according to motif position coordinate
Differential time of flight is obtained with the ratio of the speed of photon according to the range difference of two detector distance motif positions.
Further, the range difference of motif position coordinate and two detector distance motif positions in the line of response
Relation is：Wherein, Δ s represent in the line of response two detectors away from
From the range difference of motif position, x0Represent abscissa of the die body in XOY plane, y0Represent die body in XOY
The vertical coordinate of plane,Line of response angle in the counterclockwise direction with coordinate axess y positive directions is represented,
Further, according to the flight time sinogram for meeting data and time bias table to be detected, obtain and ring
The flight time aberration measurements of line, concrete formula are answered to be：Wherein Δ t ' represents line of response
Flight time aberration measurements, i represents and the flight time sinogram for meeting data made every fixed interval
Histogrammic numbering, niRepresent the frequency that numbering is light pulse in the rectangular histogram of i, Δ τiRepresent that numbering is i's
In rectangular histogram, the time of measuring of light pulse arrival detector is poor.
Further, being determined by position for detector crystal is reached according to meeting event two photon in line of response
Differential time of flight generates position time difference curve, when being generated according to the flight time aberration measurements of every bar response line
Between aberration measurements curve, by judging whether position time difference curve and time difference measurements value curve match determination
Whether flight time state quality is qualified.
Further, each data point difference on comparison position time difference curve and time difference measurements value curve
Whether in the range of given threshold, if position time difference curve and the every number on time difference measurements value curve
Strong point difference judges that flight time state quality is qualified in the range of given threshold, then；Otherwise, then judge to fly
Row time state is off quality.
Further, the given threshold is less than 25ps.
The present invention also provides a kind of PET scan device, including frame, detector, time migration correction
Die body and flight time state quality detecting system, the flight time state quality detecting system are included：
Motif position computing unit, meets detector detection data for basis and is calculated motif position；
The differential time of flight computing unit determined by position, meets detector institute for calculating according to motif position
What is detected meets the differential time of flight determined by position of two photons, a pair of detector crystals of arrival in event；
Flight time aberration measurements computing unit, for according to the flight time sinogram for meeting detector detection data
The flight time aberration measurements of line of response are obtained with time bias table to be detected；
Matching unit, the differential time of flight determined by position for basis and flight time aberration measurements whether phase
Deng judgement flight time state quality.
Further, also including flight time correction module, for detecting when the flight time state quality
When system judges that flight time state quality is unqualified, the flight time for correcting the PET scan device is inclined
Move.
Compared with prior art, it is an advantage of the current invention that：Event is met by detection and is calculated die body position
Put, can be calculated according to motif position and on every bar response line, meet two photon of event, a pair of detector crystals of arrival
Differential time of flight in XOY coordinate systems, by the flight time sinogram for meeting event and time to be detected
Compensation table obtains the time difference measurements value in every bar response line TOS coordinate systems, when differential time of flight and time
Aberration measurements are matched, then judge that the time bias meter accuracy of PET scan device meets and require, the flight time
State quality detection is qualified, when differential time of flight and time difference measurements value are mismatched, then judges PET scan
The time bias meter accuracy of device is unsatisfactory for requiring that the detection of flight time state quality is unqualified；As die body
Radioactive source need not be placed in scanning chamber center, as long as in the scan vision model of PET scan cavity
It is in enclosing, simple to operate.
【Description of the drawings】
Fig. 1 detector detection gammaphoton schematic diagrams used by the present invention；
Fig. 2 is the flight time state quality detecting system structure chart of the present invention；
Fig. 3 is PET flight time state quality detection method flow chart of the present invention；
Fig. 4 is present invention PET scan cavity used in XOY coordinate system schematic diagrams；
Fig. 5 is present invention PET scan cavity used in TOS coordinate system schematic diagrams；
Fig. 6 is that one embodiment of the invention judges flight time state quality pass result figure；
Fig. 7 is that another embodiment of the present invention judges the unqualified result figure of flight time state quality.
【Specific embodiment】
It is understandable to enable the above objects, features and advantages of the present invention to become apparent from, below in conjunction with the accompanying drawings and
Embodiment is described in detail to the specific embodiment of the present invention.
In PET imagings, in the positron and detected object body that the nucleic carried in tracer drug is launched
Electronics produces after burying in oblivion that both direction is contrary, energy for 511KeV gammaphoton (such once
Bury in oblivion and be referred to as once meeting event).The probe of PET is rearranged by several detector rings, detection
Device ring number determine PET axial direction the visual field size and fault plane number, with detector rings plane hang down
The straight axial visual field for PET.Projection line in coincidence detection is referred to as line of response (Line of as shown in Figure 1
Response, LOR), obtain it is substantial amounts of meet event after, by image reconstruction can obtain radioactive nucleus
Distribution of the element in detected object body.The imaging of flying time technology (TOF-PET) and tradition PET
The composition of process and system is essentially identical, but their time measurement system has essential distinction：In tradition
In PET, the function of time measurement system is to meet whether determination true coincidence in time window occurs, to obtain
Obtain the projection of organ and organizational structure；And in TOF-PET, the function of time measurement system is when meeting
Between position and the intensity of Radio-nuclide distribution are determined in window.TOF-PET systems, for temporal resolution
Require higher.In traditional PET, due to cannot coordinate detection is arrived on LOR annihilation event, therefore,
The annihilation event that every LOR is detected in image reconstruction is distributed in whole roads of the LOR by equal weight
On footpath.In contrast to this, in TOF-PET, gammaphoton can be detected and reaches two crystal bars
Time difference, only can determine the position at vanishing point place of falling into oblivion in theory completely by time difference.By TOF information
For image reconstruction, signal to noise ratio and contrast recovery coefficient of image etc. can be significantly improved, and then improves figure
As resolution capability, the detector efficiency of small lesion is improved, the early diagnosiss for tumor and cardiovascular disease provide heavily fortified point
Real technical guarantee.Advantage simultaneously for TOF-PET in terms of picture quality, can be relatively low
Drug dose and satisfied picture quality under the conditions of shorter sweep time, is obtained, substantially reduces radiation risk,
Significantly improve inspection efficiency.However, the measurement of time difference is determined by time difference with certain uncertainty
The position of annihilation event also have uncertainty.The accuracy measurement of flight time can improve PET scan
The performance of instrument, therefore, before PET scan imaging, it is necessary to PET flight time state-detection, i.e.,
Detection time compensates the precision of table, to improve scanning imagery quality.
PET scan device used in the present invention by detector, electronic system front end and meet system,
The cylindric frame and the die body composition as time migration correction of reconstructing system and support detector.Machine
Multiple detectors are usually provided with frame, detector is by radiosensitive scintillator crystal materials and photoelectric device group
Into its function is to detect the detection that a pair of gammaphotons being converted into during positron annihilation are hit respectively
On device ring, (the 511KeV gamma rays of positron annihilation events are linearly turned by crystal for the position of scintillation crystal bar
It is changed to flash of light), and these position signallings are converted into the signal of telecommunication (photoelectric device generation), related gammaphoton
Energy signal and the temporal information of due in be sent in follow-up front end electronics line system together.Need
Illustrate, in the present embodiment, scintillator crystal materials are lutetium yttrium orthosilicate (LYSO), naturally it is also possible to selected
Select other scintillator crystal materials such as silicic acid lutecium (LSO), gadolinium siliate (GSO), bismuth germanium oxide (BGO) with
And mix the crystalline material or semi-conducting material CdTe, CdZnTe of cerium, i.e., for crystalline material is not specific
Restriction.Photoelectric device is photomultiplier tube and silicon photomultiplier or other solid state photodetector compositions
Photodetector.Front end electronics circuit is by amplifier, flash A/D, location tables, energy scale etc.
Composition, is sent together by the location data of front end electronics circuit output, timing data and " example " efficacy data
What is entered below meets in process circuit.It should be noted that in PET detections, having three kinds of physics symbols
Conjunction form, including true coincidence, random coincidence and scattering meet.Wherein, random coincidence is one in human body
It is uniformly distributed, this homogeneous radiation distribution can increase the noise number of whole visual field, so as to reduce picture quality；
Scattering meets relevant with the form size of tested internal organs, can make pattern distortion, therefore, both meet all
It is not that we are desired.Two gammaphotons are only met simultaneously to occur simultaneously, are mutually 180 °, and two
When the ability of gammaphoton is 511KeV, true coincidence is just produced.The function of meeting system is to determine and meets,
It is whether to judge a pair of gammaphotons by once burying in oblivion what example was sent, so as to remove random signals example,
Select.Hereafter, the two basic points hit by example are really met by what Jing was selected
The coordinate Jing computer interfaces of detector scintillation crystal bar, are sent to connected computer.Computer and work station
And the function of corresponding each software kit is data acquisition, system monitoring, image reconstruction and image procossing,
Realize clinically various operations and the requirement for diagnosing.When the PET scan device of the present invention also includes a flight
Between state quality detecting system, as shown in Fig. 2 the flight time state quality detecting system includes die body position
The flight time computing unit 200 that put computing unit 100, determines by position, flight time aberration measurements meter
Calculate unit 300 and matching unit 400.
Motif position computing unit 100, meets detector for basis and (meets one of photon hit in event
To detector crystal) the data that meet be calculated motif position, it should be noted that in the present embodiment
Die body (radioactive source) can be linear radiation source, uniform barreled radioactive source, in uniform hollow tubbiness radioactive source
Any one, the motif position needs to meet die body and cylindrical shape frame concentric co-axial in being different from prior art,
As long as in the range of scan vision and parallel with cylindrical shape frame central axle.
The differential time of flight computing unit 200 determined by position, is connected with motif position computing unit 100,
A pair of spies are reached for two photons in event that meet detected by calculating and meet detector according to motif position
Survey the differential time of flight determined by position of device crystal；And by positron source, connection meets two light of event
The flight path of son forms line of response；
Flight time aberration measurements computing unit 300, with the differential time of flight computing unit determined by position
200 connections, for according to the flight time sinogram and time bias to be detected for meeting detector detection data
Table obtains the flight time aberration measurements per bar response line；
Matching unit 400, is connected with flight time aberration measurements computing unit 300, for basis by position
The conjunction for judging flight time state quality whether equal with flight time aberration measurements of the differential time of flight of decision
Lattice situation.
PET flight time state quality detection methods are as shown in figure 3, comprise the following steps：
S10, die body is placed in the range of the scan vision of PET scan cavity, the die body axis and scanning
Field range diameter parallel.It should be noted that in this particular embodiment, the position of die body is different from now
Need to meet die body and cylindrical shape frame concentric co-axial in having technology, as long as in the range of scan vision, i.e.,
Do not have strict restriction to the position of die body, such operator need not take a significant amount of time accurately adjustment mould
Body has saved the operating time to center.
S20, detection meet event, obtain meeting data, and are calculated mould according to the data that meet for obtaining
Body position.The positron-emitting radionuclides that the die body of PET scan intracavity contains can launch positron in decay,
And the proton decay in atomic nucleus is neutron and discharges positron and neutrino, quality and the electronics phase of positron
Deng electricity is identical with the electricity of electronics, and symbol is contrary.Positron is advanced before burying in oblivion in tissue
1-3mm, while producing the gammaphoton of the 511keV for being mutually 180 °.Two relative gammaphotons for pointing to
Ray two crystal counters substantially simultaneously in collision detector ring, and produce and meet trigger.Symbol
That closes the transmission of detector pick-up probe meets trigger, and determines whether in a scheduled time window
Two events are same annihilation event.If same annihilation event, then meet detector and start line of response
(LOR) processor is determining line of response that annihilation event occurs along which.If Fig. 4 is the embodiment of the present invention
PET scan cavity coordinate schematic diagram (to scan the center of cavity as zero, scanning cavity longitudinal direction tangent plane
The radius of circumference is for r), scanning cavity surrounding is distributed with multiple annular detectors, can detect and meet in event
Two photons for sending, such as detector A (xa, ya) and detector B (xb, yb) really meet event and hit
Two crystal counters, the line of two detectors is line of response, and starts flight time computer to be based on two
The individual relative reception time for meeting trigger determines that annihilation event occurs the position that point occurs along LOR, i.e.,
Radioactive source motif position (x0, y0), 0 ＜ x0＜ r, and 0 ＜ y0＜ r.
S30, according to motif position (x0, y0), calculate two photons in meeting event on every bar response line and reach one
To the differential time of flight on detector crystal, wherein line of response is by the detector crystal in detector rings and institute
State the relative detector crystal line of detector crystal to produce, motif position and Jing are obtained according to motif position
Cross the distance distance at the line of response center of the die body.As shown in figure 4, detector is arranged closely in detector rings
Zhou Shang, transverse section of the X-Y plane for PET are parallel with detection plane of a loop；Z axis are the major axis of PET,
With detector loop plane normal.With the axial visual field for PET of detector loop plane normal, and PET length
The greatest length of detectable true coincidence time in the range of axle.Axial tomography number=(number of rings × 2) -1 of PET.
Through die body C and two ends correspond to the line of response of detector A and B respectively, according to geometric operation, this
Line of response midpoint D with the distance of motif position C isIt is line of response along inverse
The angle of clockwise and coordinate axess y positive directions,According to symmetry, then in the response
On line, two detector distance motif positions apart from length difference for Δ s=2 Δs 1 isCombine and can meet two photon of event in the line of response using above-mentioned relation
Differential time of flight be：(fly) of the wherein c for the light velocity, i.e. photon
Speed.The time difference that two photons in all line of response reach correspondence detector crystal can be obtained by that analogy,
And form the Δ t time difference curves that correspondence is determined by position.
(d), according to the flight time sinogram and time bias table to be detected for meeting data, obtain per bar response
The flight time aberration measurements of line.The detection of coincidence detection module meets event, calculates each and meets two in the time
Individual photon is to the differential time of flight to a pair of detector crystals, and electronic system front end also includes that the time is inclined
Correcting module is moved, which prestores time bias table, when meeting the flight of two photons in the time in theory for a pair
Between difference reach time difference for two rightabout photons of detector record, and actual optical signal is converted into the signal of telecommunication
Transmitting procedure in, often have certain time-delay, detector can be improved by time migration correcting module
Precision.In such as Fig. 4, the actual flying time of two photons arrival detector A and detector B is poor (during measurement
Between it is poor) Δ τAB=(TA-OTA)-(TB-OTB)=(TA-TB)-(OTA-OTB), wherein OTAIt is for spy
Survey the photon compensation time of device A, OTBBe for detector B photon compensation the time, OTA-OTBCan root
Obtain according to time bias table, TAAnd TBThe logging of flight time value of detector is reached for photon.But PET
In scanning means, the precision of the time migration offset (correction value) of each detector crystal might not meet
Actual demand, therefore, it is the accuracy for definitely measuring, needs to detect the time bias table (correction chart)
Accuracy, i.e., detect to the flight time state quality of PET scan device.
In this particular embodiment, detector array is listed in collection edge in the set timeAll line of response in direction,
Obtain two-dimensional bodies to correspond toThe One Dimensional Projection distribution in direction, can be obtained by that analogy around two-dimensional bodies
Different angles projective distribution measured value, the PET measurement count form thus expressed be TOF it is sinusoidal
Figure, using time bias table (Time Offset) to be detected as input, and calculates each using sinogram
The time difference of bar response line, concrete formula is：Wherein to meeting flight time of data
Sinogram (light pulse that each detects reaches the time of detector crystal), 5-15ps takes a sub-cell (bin)
Make rectangular histogram, in the present embodiment, be preferably 10ps, i represents histogrammic numbering (i=1,2,3 ...), niTable
Show the frequency (meeting event number) that numbering is light pulse in the rectangular histogram of i, Δ τiRepresent the Nogata that numbering is i
It (is that the original time value recorded during each meets event is gone that in figure, the time of measuring of light pulse arrival detector is poor
Except the time of measuring of gained after the time migration correction value in time bias table to be measured it is poor), according to above-mentioned formula
Flight time aberration measurements Δ t ' is obtained.
For each bar response line LOR in XOY plane, all there is an Adjoint Coordinates system TOS.As schemed
Shown in 5, it is that, parallel to line of response LOR, s axles are perpendicular to line of response for t axles the characteristics of TOS coordinate systems
LOR, and through the center of the line of response, the sinusoidal index of this line of response is usedRepresent,
Wherein, for the numbering of line of response, (scope of j can determine j according to practical situation, in the present embodiment
J=1,2,3 ... 576),Represent that the corresponding coordinate axess t axles of line of response that numbering is j are square with coordinate axess y-axis
To the angle angle of-π to π scopes (will 576 deciles are divided into), and as t axles are parallel to the response
Line LOR, then hereinEqual to the angle of the line of response and coordinate axess y-axis positive direction that numbering is jsj
Represent and number the distance for being the line of response of j with the center of circle (or scanning chamber center), raAnd rbRepresent that numbering is j
Line of response two-end-point detector crystal numbering.According to above-mentioned relation, line of response LOR and coordinate system XOY
Relation translate into the relation between coordinate system TOS and coordinate system XOY.It is through die body and corresponding same
The line of response of detector is had it is a plurality of, therefore also need by Δ t ' flight time aberration measurements eachSj
Dimension is averaged (i.e.For constant or definite value, Δ t ' is sought with regard to sjRelation).The like can obtain Δ t '
Time difference measurements value curve.
(e), by the comparison position time difference curve judgement flight time whether equal with time difference measurements value curve
Whether state quality is qualified.Detailed process is：Comparison position time difference curve and time difference measurements value curve
Whether each data point difference is in the range of given threshold, if position time difference curve and time difference measurements value
Each data point difference on curve then judges that the time of PET scan device mends in the range of given threshold
Repay meter accuracy and meet requirement, flight time state quality is qualified；Otherwise, then judge PET scan device
Time bias meter accuracy is unsatisfactory for requiring that flight time state quality is unqualified.It is illustrated in figure 6 an enforcement
Flight time state quality pass result figure is judged using the inventive method in example, wherein abscissa represents response
The numbering of line, that is, represent each numbering it is corresponding with y-axis positive direction intoThe line of response of angle, vertical coordinate table
Show time value, smooth curve is the differential time of flight curve determined by position, and the curve of edge roughness is winged
Row time difference measurements value curve, for eachAll meet | Δ t- Δ t ' |≤25ps, " | | " to take absolute value,
When the time difference numerical value for being determined by position in circumferencial direction each angle value corresponding line of response is with flight
Between aberration measurements difference be both less than 25ps, two heights of curve matching, be now judged to flight time state quality
It is qualified.Judge that flight time state quality does not conform to using the inventive method in being illustrated in figure 7 another embodiment
Lattice result figure, wherein smooth curve P1It is the differential time of flight curve determined by position, edge roughness
Curve P2For flight time aberration measurements curve, two curves are determined by position multiple angle values are corresponding
There is relatively large deviation between time difference numerical value and flight time aberration measurements, be now judged to flight time state
It is off quality, then need to be corrected flight time compensation table.
In PET scan device, the time migration correction value of each detector is calculated using die body, it is existing
In technology, die body and cylindrical shape frame concentric co-axial need to be ensured, in embodiments of the present invention, at die body center
In the case that axle is non-overlapped with cylindrical shape frame central axle, i.e., die body symmetrical centre is located at PET scan device
During the middle position of field range, above-mentioned flight time matter is adopted using flight time state quality detecting system
Quantity measuring method can carry out quality testing to the flight time state of scanning means, and without the need for strict guarantee die body
Central shaft is Chong Die with cylindrical shape frame central axle, saves the central shaft overlapping correction before detection.In addition, this
The PET scan device of invention also includes flight time correction module, when the detection of flight time state quality is
When system judges that flight time state quality is unqualified, flight time bias correction is obtained, correction PET sweeps
The flight time skew of imaging apparatus.
Presently preferred embodiments of the present invention is the foregoing is only, it is not to limit the present invention, all at this
Within bright spirit and principle, any modification, equivalent substitution and improvements done etc. should be included in this
Within bright protection domain.
Claims (10)
1. a kind of PET flight time state quality detection method, it is characterised in that include：
Die body is placed in the range of the scan vision of PET scan cavity, the die body axis and scan vision
Scope diameter parallel；
Detection meets event, obtains meeting data, and meets data acquisition motif position according to described；
According to motif position, calculate meet in line of response two photon of event reach a pair of detector crystals by position
The differential time of flight of decision is put, the line of response is the flight path for meeting two photon of event；
According to the flight time sinogram for meeting data and time bias table to be detected, the flight of line of response is obtained
Time difference measurements value；
Flown by comparing the differential time of flight judgement whether equal with flight time aberration measurements determined by position
Whether row time state quality is qualified.
2. PET flight time state quality detection method as claimed in claim 1, it is characterised in that
The die body is linear radiation source, uniform tubbiness radioactive source, any one in uniform hollow tubbiness radioactive source.
3. PET flight time state quality detection method as claimed in claim 1, it is characterised in that
According to motif position, calculate and in line of response, meet being determined by position for two photon of event, a pair of detector crystals of arrival
Fixed differential time of flight detailed process is：
The range difference of two detector distance motif positions in the line of response, root are calculated according to motif position coordinate
Differential time of flight is obtained with the ratio of the speed of photon according to the range difference of two detector distance motif positions.
4. PET flight time state quality detection method as claimed in claim 3, it is characterised in that
The motif position coordinate with the range difference relation of two detector distance motif positions in the line of response is：Wherein, Δ s represents two detector distance die body position in the line of response
The range difference put, x0Represent abscissa of the die body in XOY plane, y0Represent die body in XOY plane
Vertical coordinate,Line of response angle in the counterclockwise direction with coordinate axess y positive directions is represented,
5. PET flight time state quality detection method as claimed in claim 1, according to meeting data
Flight time sinogram and time bias table to be detected, obtain line of response flight time aberration measurements, tool
Body formula is：Wherein Δ t ' represents the flight time aberration measurements of line of response, and i represents right
Meet the flight time sinogram of data every the made histogrammic numbering of fixed interval, niRepresent that numbering is i
Rectangular histogram in frequency in light arteries and veins, Δ τiRepresent that numbering is that light pulse reaches detector in the rectangular histogram of i
Time of measuring is poor.
6. PET flight time state quality detection method as claimed in claim 1, it is characterised in that
According to the flight time bad student determined by position for meeting two photon of event arrival detector crystal in line of response
Into position time difference curve, time difference measurements value is generated according to the flight time aberration measurements of every bar response line bent
Line, by judging whether position time difference curve and time difference measurements value curve match determination flight time state
Whether quality is qualified.
7. PET flight time state quality detection method as claimed in claim 6, it is characterised in that
Whether each the data point difference on comparison position time difference curve and time difference measurements value curve is in setting threshold
In the range of value, if position time difference curve is being set with each the data point difference on time difference measurements value curve
Determine in threshold range, then judge that flight time state quality is qualified；Otherwise, then judge flight time state matter
Amount is unqualified.
8. PET flight time state quality detection method as claimed in claim 7, it is characterised in that
The given threshold is less than 25ps.
9. die body and fly that a kind of PET scan device, including frame, detector, time migration are corrected
Row time state quality detecting system, the flight time state quality detecting system are included：
Motif position computing unit, meets detector detection data for basis and is calculated motif position；
The differential time of flight computing unit determined by position, meets detector institute for calculating according to motif position
What is detected meets the differential time of flight determined by position of two photons, a pair of detector crystals of arrival in event；
Flight time aberration measurements computing unit, for according to meet detector detection data flight time just
String figure and time bias table to be detected obtain the flight time aberration measurements of line of response；
Matching unit, the differential time of flight determined by position for basis and flight time aberration measurements whether phase
Deng judgement flight time state quality.
10. PET scan device as claimed in claim 9, it is characterised in that also including the flight time
Correction module, for judging that flight time state quality does not conform to when the flight time state quality detecting system
During lattice, the flight time skew of the PET scan device is corrected.
Priority Applications (8)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510603207.2A CN106539591B (en) | 2015-09-21 | 2015-09-21 | PET flight time state quality detection method and PET scanning device |
| US15/171,877 US10663608B2 (en) | 2015-09-21 | 2016-06-02 | System and method for calibrating a PET scanner |
| KR1020187010475A KR102709104B1 (en) | 2015-09-21 | 2016-09-14 | System and method for calibrating a ET scanner |
| PCT/CN2016/099078 WO2017050180A1 (en) | 2015-09-21 | 2016-09-14 | System and method for calibrating a pet scanner |
| EP16848062.2A EP3352671B1 (en) | 2015-09-21 | 2016-09-14 | System and method for calibrating a pet scanner |
| JP2018514820A JP6828024B2 (en) | 2015-09-21 | 2016-09-14 | Systems and methods for calibrating PET scanners |
| US16/882,556 US11619755B2 (en) | 2015-09-21 | 2020-05-25 | System and method for calibrating a PET scanner |
| US18/295,239 US20230243988A1 (en) | 2015-09-21 | 2023-04-03 | System and method for calibrating a pet scanner |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510603207.2A CN106539591B (en) | 2015-09-21 | 2015-09-21 | PET flight time state quality detection method and PET scanning device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106539591A true CN106539591A (en) | 2017-03-29 |
| CN106539591B CN106539591B (en) | 2020-06-26 |
Family
ID=58364255
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510603207.2A Active CN106539591B (en) | 2015-09-21 | 2015-09-21 | PET flight time state quality detection method and PET scanning device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106539591B (en) |
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106539593A (en) * | 2016-11-30 | 2017-03-29 | 上海联影医疗科技有限公司 | Medical image system and its method for testing performance |
| CN108542412A (en) * | 2017-04-24 | 2018-09-18 | 沈阳东软医疗系统有限公司 | Image rebuilding method and PET device |
| CN109363709A (en) * | 2018-12-19 | 2019-02-22 | 上海联影医疗科技有限公司 | PET data acquisition method and PET system |
| CN109567853A (en) * | 2019-01-17 | 2019-04-05 | 上海联影医疗科技有限公司 | Method and device, the PET scan equipment of acquisition time scale value |
| CN109632837A (en) * | 2018-12-28 | 2019-04-16 | 中国原子能科学研究院 | The automatic time scaling method of associated particle neutron detection |
| CN110432923A (en) * | 2019-08-14 | 2019-11-12 | 上海联影医疗科技有限公司 | Image rebuilding method, device, computer equipment and storage medium |
| CN110584698A (en) * | 2019-09-23 | 2019-12-20 | 上海联影医疗科技有限公司 | Detector quality control validation method, detector quality control validation device, computer equipment and storage medium |
| CN110988974A (en) * | 2019-12-12 | 2020-04-10 | 明峰医疗系统股份有限公司 | Calibration system for flight time of PET (positron emission tomography) detector of curved radiation source |
| CN111080734A (en) * | 2019-11-25 | 2020-04-28 | 中国科学院深圳先进技术研究院 | Method and terminal for processing Positron Emission Tomography (PET) data |
| CN111568452A (en) * | 2020-05-25 | 2020-08-25 | 上海联影医疗科技有限公司 | PET system state detection method, PET system state detection device, computer equipment and storage medium |
| CN111839569A (en) * | 2020-07-29 | 2020-10-30 | 明峰医疗系统股份有限公司 | PET time resolution measuring, calculating and reconstructing system and method |
| CN112102426A (en) * | 2020-08-28 | 2020-12-18 | 上海联影医疗科技股份有限公司 | Background coincidence event judging and selecting method, device, equipment and readable storage medium |
| CN112450959A (en) * | 2020-12-08 | 2021-03-09 | 鲁东大学 | Single-mode PET imaging system and method capable of achieving accurate positioning |
| CN112807009A (en) * | 2021-02-02 | 2021-05-18 | 上海联影医疗科技股份有限公司 | Radioactive source positioning method, system, equipment, electronic device and storage medium |
| WO2021102614A1 (en) * | 2019-11-25 | 2021-06-03 | 中国科学院深圳先进技术研究院 | Method and terminal for processing positron emission tomography (pet) data |
| CN112971824A (en) * | 2021-02-08 | 2021-06-18 | 上海联影医疗科技股份有限公司 | PET dynamic image scanning method, device and computer equipment |
| CN112998736A (en) * | 2021-02-26 | 2021-06-22 | 中派科技（深圳）有限责任公司 | Time correction system and time correction method of scanning device |
| CN112998735A (en) * | 2021-02-26 | 2021-06-22 | 中派科技（深圳）有限责任公司 | System and method for reconstructing image by scanning device |
| CN113116371A (en) * | 2019-12-30 | 2021-07-16 | 佳能医疗系统株式会社 | PET device and calibration method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070131857A1 (en) * | 2002-11-07 | 2007-06-14 | Thompson Christopher J | Instrument and method to facilitate and improve the timing alignment of a pet scanner |
| US20070205368A1 (en) * | 2005-09-14 | 2007-09-06 | General Electric Company | Method and system for calibrating a time of flight positron emission tomography system |
| US20110127413A1 (en) * | 2009-12-01 | 2011-06-02 | Siemens Medical Solutions Usa, Inc. | Systems and Methods for Calibrating Time Alignment For A Positron Emission Tomography (PET) System |
| CN102283665A (en) * | 2010-06-17 | 2011-12-21 | 株式会社东芝 | Nuclear medicine imaging apparatus, and nuclear medicine imaging method |
| US8796637B1 (en) * | 2013-05-24 | 2014-08-05 | Kabushiki Kaisha Toshiba | Timing calibration for time-of-flight (TOF) PET using positron-emitting isotopes and annihilation targets |
| WO2015028602A1 (en) * | 2013-08-30 | 2015-03-05 | Uniwersytet Jagiellonski | A method for calibration of tof-pet detectors using cosmic radiation |
-
2015
- 2015-09-21 CN CN201510603207.2A patent/CN106539591B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070131857A1 (en) * | 2002-11-07 | 2007-06-14 | Thompson Christopher J | Instrument and method to facilitate and improve the timing alignment of a pet scanner |
| US20070205368A1 (en) * | 2005-09-14 | 2007-09-06 | General Electric Company | Method and system for calibrating a time of flight positron emission tomography system |
| US20110127413A1 (en) * | 2009-12-01 | 2011-06-02 | Siemens Medical Solutions Usa, Inc. | Systems and Methods for Calibrating Time Alignment For A Positron Emission Tomography (PET) System |
| CN102283665A (en) * | 2010-06-17 | 2011-12-21 | 株式会社东芝 | Nuclear medicine imaging apparatus, and nuclear medicine imaging method |
| US8796637B1 (en) * | 2013-05-24 | 2014-08-05 | Kabushiki Kaisha Toshiba | Timing calibration for time-of-flight (TOF) PET using positron-emitting isotopes and annihilation targets |
| WO2015028602A1 (en) * | 2013-08-30 | 2015-03-05 | Uniwersytet Jagiellonski | A method for calibration of tof-pet detectors using cosmic radiation |
Cited By (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106539593A (en) * | 2016-11-30 | 2017-03-29 | 上海联影医疗科技有限公司 | Medical image system and its method for testing performance |
| CN108542412A (en) * | 2017-04-24 | 2018-09-18 | 沈阳东软医疗系统有限公司 | Image rebuilding method and PET device |
| CN108542412B (en) * | 2017-04-24 | 2021-11-09 | 东软医疗系统股份有限公司 | Image reconstruction method and PET device |
| CN109363709B (en) * | 2018-12-19 | 2020-11-24 | 上海联影医疗科技有限公司 | PET data acquisition method and PET system |
| CN109363709A (en) * | 2018-12-19 | 2019-02-22 | 上海联影医疗科技有限公司 | PET data acquisition method and PET system |
| CN109632837A (en) * | 2018-12-28 | 2019-04-16 | 中国原子能科学研究院 | The automatic time scaling method of associated particle neutron detection |
| CN109567853A (en) * | 2019-01-17 | 2019-04-05 | 上海联影医疗科技有限公司 | Method and device, the PET scan equipment of acquisition time scale value |
| CN110432923A (en) * | 2019-08-14 | 2019-11-12 | 上海联影医疗科技有限公司 | Image rebuilding method, device, computer equipment and storage medium |
| CN110584698A (en) * | 2019-09-23 | 2019-12-20 | 上海联影医疗科技有限公司 | Detector quality control validation method, detector quality control validation device, computer equipment and storage medium |
| WO2021102614A1 (en) * | 2019-11-25 | 2021-06-03 | 中国科学院深圳先进技术研究院 | Method and terminal for processing positron emission tomography (pet) data |
| CN111080734A (en) * | 2019-11-25 | 2020-04-28 | 中国科学院深圳先进技术研究院 | Method and terminal for processing Positron Emission Tomography (PET) data |
| US20220283325A1 (en) * | 2019-11-25 | 2022-09-08 | Shenzhen Institutes Of Advanced Technology Chinese Academy Of Sciences | Method and terminal device for processing positron emission tomography data |
| CN111080734B (en) * | 2019-11-25 | 2023-10-20 | 中国科学院深圳先进技术研究院 | Method and terminal for processing Positron Emission Tomography (PET) data |
| CN110988974A (en) * | 2019-12-12 | 2020-04-10 | 明峰医疗系统股份有限公司 | Calibration system for flight time of PET (positron emission tomography) detector of curved radiation source |
| CN113116371A (en) * | 2019-12-30 | 2021-07-16 | 佳能医疗系统株式会社 | PET device and calibration method |
| CN113116371B (en) * | 2019-12-30 | 2024-04-19 | 佳能医疗系统株式会社 | PET device and calibration method |
| CN111568452A (en) * | 2020-05-25 | 2020-08-25 | 上海联影医疗科技有限公司 | PET system state detection method, PET system state detection device, computer equipment and storage medium |
| CN111839569B (en) * | 2020-07-29 | 2022-11-01 | 明峰医疗系统股份有限公司 | PET time resolution measuring, calculating and reconstructing system and method |
| CN111839569A (en) * | 2020-07-29 | 2020-10-30 | 明峰医疗系统股份有限公司 | PET time resolution measuring, calculating and reconstructing system and method |
| CN112102426B (en) * | 2020-08-28 | 2024-03-26 | 上海联影医疗科技股份有限公司 | Background coincidence event judging and selecting method, device and equipment and readable storage medium |
| CN112102426A (en) * | 2020-08-28 | 2020-12-18 | 上海联影医疗科技股份有限公司 | Background coincidence event judging and selecting method, device, equipment and readable storage medium |
| CN112450959A (en) * | 2020-12-08 | 2021-03-09 | 鲁东大学 | Single-mode PET imaging system and method capable of achieving accurate positioning |
| CN112450959B (en) * | 2020-12-08 | 2023-08-15 | 鲁东大学 | Single-mode PET imaging system and method capable of realizing accurate positioning |
| CN112807009A (en) * | 2021-02-02 | 2021-05-18 | 上海联影医疗科技股份有限公司 | Radioactive source positioning method, system, equipment, electronic device and storage medium |
| CN112971824A (en) * | 2021-02-08 | 2021-06-18 | 上海联影医疗科技股份有限公司 | PET dynamic image scanning method, device and computer equipment |
| CN112998735B (en) * | 2021-02-26 | 2022-09-02 | 中派科技（深圳）有限责任公司 | System and method for reconstructing image by scanning device |
| CN112998735A (en) * | 2021-02-26 | 2021-06-22 | 中派科技（深圳）有限责任公司 | System and method for reconstructing image by scanning device |
| CN112998736A (en) * | 2021-02-26 | 2021-06-22 | 中派科技（深圳）有限责任公司 | Time correction system and time correction method of scanning device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106539591B (en) | 2020-06-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106539591A (en) | PET flight time state quality detection methods and PET scan device | |
| US7402807B2 (en) | Method for reducing an electronic time coincidence window in positron emission tomography | |
| CN103908280B (en) | The method of PET scatter corrections | |
| CN1307431C (en) | Resolution enhancement for dual head gamma camera | |
| CN107110983A (en) | Coincidence correction based on pixel | |
| EP3033637B1 (en) | A hybrid method based on simulation and experimental data to normalize pet data | |
| CN105496436B (en) | For the time-correcting method and device of PET device | |
| US20110127436A1 (en) | Radiation tomography apparatus | |
| EP2748641A2 (en) | Data-driven optimization of event acceptance/rejection logic | |
| CN106556857A (en) | Response location localization method and device | |
| CN108338805B (en) | Time correction device for PET system | |
| US6403960B1 (en) | Correction for spatial variations of deadtime of a monolithic scintillator based detector in a medical imaging system | |
| CN205826876U (en) | Positron emission tomography | |
| CN106353786A (en) | Method and device for detecting flight time performance of positron emission tomography system | |
| CN101530330B (en) | Nuclear medicine diagnosis device and method, form tomography diagnosis device and method | |
| EP3039442B1 (en) | A hybrid tof-pet/mri tomograph | |
| CN108152849A (en) | Radiation detecting apparatus | |
| Stolin et al. | Preclinical positron emission tomography scanner based on a monolithic annulus of scintillator: initial design study | |
| US9279892B2 (en) | Systems and methods for scintillators having polished and roughened surfaces | |
| CN112998737B (en) | Time offset correction system and method for scanning device | |
| Domingo-Pardo et al. | A Position Sensitive $\gamma $-Ray Scintillator Detector With Enhanced Spatial Resolution, Linearity, and Field of View | |
| CN103126701A (en) | Positron emission computed tomography apparatus and image processing apparatus | |
| Xi et al. | Optimization of the SiPM pixel size for a monolithic PET detector | |
| Stefaan et al. | State of the art in total body PET | |
| Sattler | 22 Clinical Imaging Molecular PET/MRI Hybrid |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai. Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai. Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |