CN108763758A - A kind of GATE emulation modes of non-complete ring-type PET rotation sweep patterns - Google Patents
A kind of GATE emulation modes of non-complete ring-type PET rotation sweep patterns Download PDFInfo
- Publication number
- CN108763758A CN108763758A CN201810531389.0A CN201810531389A CN108763758A CN 108763758 A CN108763758 A CN 108763758A CN 201810531389 A CN201810531389 A CN 201810531389A CN 108763758 A CN108763758 A CN 108763758A
- Authority
- CN
- China
- Prior art keywords
- crystal
- detector
- ring
- gate
- emulation
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Nuclear Medicine (AREA)
Abstract
The invention discloses a kind of GATE emulation modes of non-complete ring-type PET rotation sweep patterns, and steps are as follows for this method：Model and parameter prepare before emulation, including establishing non-complete ring-type PET geometrical models under GATE platforms, radiate source model setting, crystal probe parameter setting, detector single rotation angle calculates and is arranged with number of revolutions, and physical action is arranged between particle, analog to digital conversion circuit parameter setting, output file format is arranged, initial sweep time, time slice, the setting of total simulation time；Simulation stage determines detector initial position, completes to carry out rotating and repeating the above scanning to non-complete detector rings again after static scanning under the position, terminates to emulate when overall sweep time reaches preset value.Emulation mode through the invention completes PET device building plan for a small amount of crystal probe and provides simulating, verifying, contributes to optimization and the equipment research and development of rotary PET acquisition pattern, reduces equipment manufacturing costs.
Description
Technical field
The invention belongs to nuclear physics emulation experiment fields, more particularly to a kind of being directed to non-complete ring-type PET rotation sweeps
The GATE emulation modes of pattern.
Background technology
In recent years, radiation detection technology is in industry, agricultural, health care, scientific research extensive use, due to various functions
Powerful, complicated radiation detector assembly and spectrometer are constantly studied, develop, the detector R and D method of early stage with
And it can not adaption demand.Various each scientific research institutions at present and Appliance Development Corp. generally use computer simulation auxiliary detector
Design, Monte Carlo method as a kind of random counting method, to carry can physical process of the particle in detector medium it is similar, because
This this method has become scientific research personnel and engineer into the important means of row detector Computer Aided Design.
Positron emission tomography (Positron Emission Tomography, abbreviation PET) is current medicine neck
The imaging method of domain most advanced level, it is based on positron annihilation principle, can be according to radiotracer in bio-tissue
Distribution and noninvasive, quantitative, dynamically to assess various organs metaboilic level, biochemical reaction and functional activity.And with PET
The continuous development of technology, application field are also no longer limited only to medical domain, are imaged such as the PET of plant research
Instrument, the PET imaging devices of industrial circle particles track all embody the application value that PET is imaged on each field.Positron annihilation
The γ photons of generation are small by material and such environmental effects to there is very strong penetration capacity, and positron radionuclide is tagged to properly
Carrier such as liquid, gas or solid particle medium on, inside injection device, further according to periphery PET device detected
Imaging is, it can be achieved that inside industrial equipment is labeled the positioning of particle.But since medical PET detections ring once builds completion at present,
Layout is difficult change, and domain detector number is more, and hardware is complicated, can not be suitable for complicated changeable industrial equipment；
Simultaneously as PET sensors are expensive, the examination of detector design of many in conceptual phase especially feasibility verification
It tests, can not directly be unfolded to test by PET hardware devices.
Invention content
Goal of the invention：For the above problem of the existing technology, the present invention is intended to provide a kind of non-complete ring-type PET rotations
Turn the GATE emulation modes of scan pattern, proposes that the non-complete PET device of rotation sweep completion, which is added, using crystal probe builds
Scheme, sensor logarithm keep it more flexible in industrial application without limitation；It is simulated by GATE emulation modes true
PET device rotation sweep pattern provides simulating, verifying and optimization reference for rotary PET device building plan, is rotary
PET device research and development provide booster action.
Technical solution：To achieve the purpose of the present invention, the technical solution adopted in the present invention is：A kind of non-complete ring-type
The GATE emulation modes of PET rotation sweep patterns, steps are as follows for the emulation mode：
(1) simulated environment is GATE platforms under Linux Ubuntu systems, and pet detector model is established on platform, is visited
There are many selections for sensor logarithm in survey device, and arrange in cyclic annular at equal intervals；
(2) the two dimensional gray slice map of the radiation source model of preset number is merged to the radioactive source that is converted to GATE and can be read
Model format adjusts moulded dimension, is lived with nucleic according to the distribution of the material of pixel grey scale and GATE material data lab setting models
Spend the configuration documentation of distribution；
(3) crystal probe parameter is set, determines sensor internal crystal coding rule；
(4) PET is solved according to crystal probe size and detects ring rotation parameter, selection monocrystalline is required according to imaging resolution
The integral multiple that body item corresponds to central angle is calculated as rotation step angle by sensor logarithm and the rotation stepping angle of setting
Obtain the number of revolutions that detection ring completes whole circumference scanning；
(5) analog to digital conversion circuit parameter, particle interaction, output file format in pet detector are configured；
(6) detector initial sweep time, time slice and overall sweep time are set；
(7) simulation stage determines detector initial position, and is scanned according to time slice, after completing position scanning
Rotation is carried out to non-complete detector rings to be scanned next time, stops imitating when overall sweep time reaches preset value
Very；
(8) after completing entire scanning process, angle modification is carried out to the event data that meets of output, for counting after amendment
According to recombination, the sinogram after recombination is used for image reconstruction.
Wherein, in step (1), three-dimensional system of coordinate is established in GATE emulation platforms, builds pet detector model, the detection
Device model includes detection ring T, n crystal probe M0~Mn, and the detection that it is D millimeters in internal diameter that sensor, which is spacedly distributed,
On ring T, and face where ring is parallel with xOy planes, and wherein M0 numbers for home sensor, in terms of z-axis negative direction, sensing system
It is incremented by clockwise by M0~Mn.
Wherein, in step (2), the two dimension slicing figure for radiating source model is merged and is converted to " .img ", " .hdr " format,
Setting radioactive source moulded dimension is adjusted by Pixel Dimensions, the different material in different pixel grey scale region setting models,
Different radioactive source activities, final GATE platforms can be directly read with the radiation source model of " .img ", the storage of " .hdr " format.
Wherein, in step (3), sensor parameters, including sensor internal crystal number, array way and ruler are set
Very little, crystal probe M0~Mn local arrays are L*H, and wherein the directions y number is L, and the directions z number is H, crystal bottom width
It is w millimeters, the directions x are the length of crystal item, and crystal item is incremented by again the incremental sequence in the directions z with the first directions y and is numbered,
Middle y-axis forward direction is row, and z-axis forward direction is row.
Wherein, in step (4), the integral multiple β setting rotation sweep stepping angles θ of central angle alpha is corresponded to according to monocrystal itemn,
Pass through the stepping angle θ of settingnWith crystal probe logarithm n, the number of revolutions that detection ring completes whole circumference scanning is calculated
M, computational methods are as follows：It is D millimeters to detect ring diameter, and single crystal width is w millimeters, by approximation method, can be obtained individually brilliant
The corresponding central angle of body itemSingle rotation angle θn=β α °, number of revolutions
Wherein, in step (6), (7), setting emulation initial sweep time, time slice, overall sweep time and detection
Device initial scanning position, emulation initial time are set as Tstart, sweep time of the time slice by detector in each position
TsliceIt determines, overall sweep time TtotalThe number M rotated by detector is determined, meets Ttotal=M*Tslice, determine emulation rank
Section detector initial position R0, according to TsliceIt is scanned, detector is rotated and completed down after completing position scanning
Single pass, multiple scanning stop emulation until overall sweep time reaches setting value.
Wherein, in step (8), after completing scanning, the document that storage meets event data is obtained, it is virtual by establishing one
Complete circle, internal diameter is D millimeter, on annulus connected between crystal item and size also with crystal item on non-complete detector rings T
Size is consistent, and crystal item number in initial position is 0, is incrementally virtually numbered by number clockwise in terms of z-axis negative sense.
Wherein, in step (8), meeting in event data for output includes the energy for receiving γ photons, time, crystal
Number information, line of response are two crystal lines for being nearly simultaneously received γ photons, by detect ring rotate through it is total
Angle is combined with crystal number information, is mapped with virtual ring number, by crystal item weight on crystal position and virtual ring
The number of conjunction is modified to the number on virtual ring, and the event data that meets that detector is rotated to M acquisition all carries out crystal item
Number is corrected, and revised number can reflect position of the line of response on virtual annulus, is then used to meeting event data
Single slice of data reassembly algorithm, is stored in sinogram.
Advantageous effect：Compared with prior art, technical scheme of the present invention has the advantages that：
(1) through the invention, it is provided to a certain degree for various number crystal probes completion PET device building plan imitative
True verification, proposition is the scan pattern for rotating step value with monocrystal integral multiple, reduces line of response to a certain extent
Lack angle.
(2) every detector performance index that emulation obtains will be researched and developed for rotary PET device and provide booster action, help
In rotating acquisition model-based optimization, increase pet detector flexibility.
Description of the drawings
Fig. 1 is the GATE emulation mode implementing procedure figures of non-complete ring-type PET rotation sweep patterns；
Fig. 2 is non-complete ring-type PET sensor distributed model schematic diagrames under GATE emulation platforms；
Fig. 3 is radiation source model and two dimension slicing schematic diagram in GATE emulation platforms；
Fig. 4 is crystal probe crystal inside array and its array code rule schematic diagram in Fig. 2；
Fig. 5 is the simulation process schematic diagram of PET model rotating acquisition patterns；
Fig. 6 is based on crystal code modification principle schematic；
Fig. 7 is that the sinogram that the domain is obtained with data recombination algorithm of the non-complete pet detector Jing Guo angle modification schemes.
Specific implementation mode
Technical scheme of the present invention is further described with reference to the accompanying drawings and examples.
Embodiment 1
The invention discloses a kind of GATE emulation modes of incomplete ring-type PET rotation sweep patterns, method flow such as Fig. 1
Shown, ring-type PET models are arranged under GATE platforms before emulating, radiation source model setting, crystal probe parameter setting, detector
Single rotation angle calculates and is arranged with number of revolutions, and physical action is arranged between particle, analog to digital conversion circuit parameter setting, output
File format is arranged, the setting emulation initial sweep time, when rotation each angle the static scanning time, i.e. time slice, and
Overall sweep time is determined according to number of revolutions and the product of time slice.Simulation stage determines detector initial position, completes
Rotate simultaneously multiple scanning to non-complete detector rings after static scanning under the position, is set in advance when overall sweep time reaches
Definite value terminates to emulate, and acquisition meets event data, and data recombination is realized after carrying out angle modification.The specific step of the emulation mode
It is rapid as follows：
(1) pet detector model is established on emulation platform, there are many selections for sensor logarithm, and are spacedly distributed.Such as
Shown in Fig. 2 (a), non-complete detection ring T1, crystal probe module M0~M11；A three-dimensional seat is established in GATE emulation platforms
Mark system, this coordinate system contribute to the spatial position adjustment of subsequent probe device, radioactive source, generate detector model also in this coordinate system
Middle display.Fig. 2 (a) is the pet detector model of 6 pairs of crystal probe M compositions, and crystal probe M0~M11 is at equal intervals
It being distributed on the ring that an internal diameter is 100mm, and face where ring is parallel with xOy planes, wherein M0 numbers for home sensor, from
Z-axis negative direction sees that sensing system number is incremented by clockwise by M0~M11.
(2) it is superimposed by the radioactive source two dimensional gray slice map of certain amount, establishes radioactive source three-dimensional voxel model, and with
Pixel grey scale is distributed material distribution is arranged with nucleic activity.Radiating source model is built for analog study object Nuclear analysis
Vertical model, as shown in figure 3, S is the radioactive source threedimensional model established, it is characterized in that a cube, contains 4 different shapes
Through-hole, respectively cylindroid, big cylinder, small column and square, GATE platforms can not directly read 3 dimensional format file, L
For the slice map made for this radioactive source threedimensional model, slice map is gray-scale map, and pixel 128*128 gives tacit consent to Pixel Dimensions
For pixel/mm, moulded dimension can be adjusted by medical software Set scale ruler, after the superposition of certain amount slice map, passed through
Medical software conversion generates the model file format " .hdr " and " .img " that GATE platforms can be read.
Cast material, Nuclear analysis are divided into two parts in example：Aluminum shell without nucleic, black portions, gray value
For 0 and the water containing nucleic, i.e. white portion, gray value 255.Radiation source model acquiescence when importing is in three-dimensional system of coordinate
First quartile, therefore to model is placed in other positions in coordinate, then need to be adjusted by coordinate transform.Finally, divide
Not She Zhi radiation source material configuration documentation and nucleic activity distribution configuration documentation complete model foundation, it is flat that material directly invokes GATE
Material database in platform.
Wherein " Alphantom.dat " is material configuration documentation, and " ActivityRange.dat " is activity configuration documentation.
Content example is as follows in configuration documentation：
“Alphantom.dat”
2 # material categories
0~5Aluminium # gray values have been defined in the material of [0,5] in material database
200~255Water # gray values have been defined in the material of [200,255] in material database
“ActivityRange.dat”
2 # nucleic value numbers
0~5 0 # gray values are 0 in the activity of [0,5] pixel
200~255 2000.Bq # gray values are 2000Bq in the activity of [200,225] pixel
(3) sensor parameters, including sensor internal crystal size, array way and number are set.Such as Fig. 4 (a), sheet
Crystal probe M0 local arrays are 6*16 in example, and it is to consider shown in 16, Fig. 4 (b) that wherein the directions y number, which is the directions 6, z number,
Gap crystal bottom width is 2.15mm, and the directions x are the length of crystal item, can use 15mm~20mm, positive vertical with x-axis
And intersect at crystal array center p points.In Fig. 4, crystal probe module M, monocrystal c, crystal System Number c0~
c95。
For sensor M0 crystal inside array coding rules, as shown in Fig. 4 (b), c0~c95 is crystal System Number
0~95, crystal item number increases along y-axis (providing that this direction is row) positive and z-axis (providing that this direction is row) is positive, and first
The directions y be incremented by again the directions z be incremented by, it is similar with M0, M1~M11 crystal inside items number in terms of z-axis negative direction first column direction with
It is incremented by clockwise, then increases toward z-axis forward direction.
(4) detector rotation parameter is solved according to crystal probe size and logarithm, obtains single rotation angle i.e. stepping
Angle.It is less in detector number, in order to reduce the incompleteness of line of response angle, with sensor internal single crystal item
Corresponding central angle alpha is minimum step value, is required according to different imaging resolutions, 1~4 times of α can be selected for step value, α
It can shorten sweep time under the premise of meeting picture quality when larger.Wherein advised according to detection ring size, crystal probe
Lattice calculate rotation angle, if detection ring diameter is D millimeters, single crystal width is w millimeters, then, can by approximate calculation
Obtain the corresponding central angle of single crystal item：
It is β to require the crystal item number for selecting single to rotate by imaging precision and overall sweep time limitation, then single revolves
Corner is
θn=β α ° (2)
Line of response missing is excessive in order to prevent, 1≤β≤4.
Number of revolutions is wherein calculated according to crystal probe logarithm and rotation stepping angle, if existing N is to sensor and waits
It is distributed on ring, rotation angle is above-mentioned θn, then number of revolutions
Special circumstances are when N forms a complete detection ring to sensor, then M=0.
By parameter in step (3) it is found that w=2.15mm in the present embodiment, takes D=100mm, N=6, β=3 are selected, then root
Single rotation angle θ can be obtained according to formula (1), formula (2)n=7.5 °, number of revolutions M=4 is acquired according to formula (3).
(5) three analog to digital conversion circuit parameter, particle interaction, output file format partial parameters need respectively according to reality
Border hardware performance proposes meter detector performance index to simulate, if the problems such as research scatters, data processing means are set
It sets.Certain the crystal probe energy resolution selected in this example is about 15%, hardware circuit temporal resolution 600ps, when meeting
Between window and energy window be respectively 10ns, 350Kev~650Kev, then related simulation parameter can be set accordingly；Hardware requirement rotation is swept
Retouch total time, if research scattering problems, data storage facility etc. determine particle interaction and the output of simulation process
Format is arranged, and positron annihilation, photoelectric effect, compton effect,scattering, Rayleigh scattering etc. are mainly provided in this example, is only exported
It is ROOT formats to meet event and format.
(6) detector initial sweep time, time slice and overall sweep time are set.Emulation initial time is set as
Tstart, sweep time T of the time slice by detector in each positionsliceIt determines, overall sweep time TtotalIt is revolved by detector
The number M turned is determined, meets formula (4)：
Ttotal=M*Tslice (4)
Scanning initial time is traditionally arranged to be 0 in a practical situation, time slice T in the present embodimentsliceIt is set as 10s,
Know that number of revolutions is 4, then T can determine according to formula (4)totalFor 40s, simulation time relevant parameter is configured accordingly.
(7) simulation stage detector initial position is determined, according to TsliceIt is scanned, completes after the position is scanned to detection
Device is rotated and completes to scan next time, stops emulation until overall sweep time reaches setting value.As shown in figure 5, four kinds
Detection ring R0~R3 under rotation position, (a) is simulation model original state R0, after completing R0 emulation, in next timeslice
Initial period, model R0 is rotated to obtain model R1 for 7.5 ° with rotation angle, temporally piece TsliceIt is scanned, and then
It rotates to obtain R2, R3 with same step, until overall sweep time stops emulation when reaching the 40s of preset value.
(8) after completing scanning, the result emulated is largely to meet event data, and it includes to connect to meet in event data
The energy of γ photons, time, crystal number information are received, line of response is two crystal items for being nearly simultaneously received γ photons
Line.Extraction meets the crystal item number data in event, and amendment is numbered according to rotary angle position, finally recombinates
To sinogram.Consider original single slice graph display algorithm, event is determined for compliance with according to crystal probe number M and crystal item number c
Meet position of the event in individual layer annulus, such as Fig. 5, the phase allomeric number of original M0 can be in annulus after rotating acquisition
Upper multiple positions occur, at this time just can not be according to the unique positioning of number.In Fig. 6 (a), crystal probe M0 after rotation, crystal
Item number c0~c5 is repetitively appearing in multiple positions of circumference.
The present invention is D millimeters of complete circle by establishing an internal diameter also, it is closely coupled between crystal item on annulus and
Size is also consistent with the crystal size in non-complete detector rings, and crystal item number in initial position is 0, in terms of z-axis negative sense
It is incrementally virtually numbered by number clockwise.Sensor crystal item number in R1_M0, R2_M0, R3_M0 in 6 (a) is reflected
It is mapped to VR1_M0, VR2_M0, VR3_M0 sensor crystal item in 6 (b) to number, approximation method is used, by the crystal on annulus
Item is numbered to be corresponded with position.In Fig. 6, R0_M0, R1_M0, R2_M0, R3_M0 are respectively under four kinds of rotation positions
Sensor assembly M0, VR0_M0, VR1_M0, R2_M0, R3_M0 are respectively that the crystal number under four kinds of rotation positions is revised
Sensor assembly M0.
(9) similarly, crystal probe M1~M11, crystal c6~c95 can also use same mapping relations one a pair of
It should get up.After obtaining each postrotational mapping relations, multiple rotary scanning is met into event data amendment, superposition recombination is just
It can obtain sinogram.Fig. 7 (a) is the sinogram after single sweep operation recombination under complete circle, and Fig. 7 (b) is that 1 single of embodiment is swept
The sinogram after recombination is retouched, Fig. 7 (c) is the sinogram that embodiment 1 corrects that reassembly algorithm is obtained using number, it is seen that (c) and (a)
Effect due to being incomplete acquisition in figure more closely, (b) be lost many line of response.
Embodiment 2
In specific implementation process, the crystal probe logarithm in PET detection rings can also be other numbers, detector hole
Diameter, crystal size can also according to demand depending on.The present embodiment is compared with embodiment 1, and major embodiment is the emulation mode
Under sensor logarithm, detector rings aperture, crystal size different situations, the calculating of simulation process rotation parameter is specific imitative
True method and step is as follows：
(1) pet detector model is established on emulation platform, there are many selections for sensor logarithm, and are spacedly distributed.Such as
Shown in Fig. 2 (b), non-complete detection ring T2, crystal probe module M0~M11；A three-dimensional seat is established in GATE emulation platforms
Mark system, this coordinate system contribute to the spatial position adjustment of subsequent probe device, radioactive source, generate detector model also in this coordinate system
Middle display.Fig. 2 (b) is the pet detector model of 4 pairs of crystal probe M compositions, and crystal probe M0~M11 is at equal intervals
It being distributed on the ring that an internal diameter is 100mm, and face where ring is parallel with xOy planes, wherein M0 numbers for home sensor, from
Z-axis negative direction sees that sensing system number is incremented by clockwise by M0~M11.
The radiation source model that step (2) is established is in the same manner as in Example 1, and step is consistent.
(3) sensor parameters, including sensor internal crystal size, array way and number are set.Crystal in this example
Sensor M0 local arrays are still 6*16, it is 16 that wherein the directions y number, which is the directions 6, z number, not with Fig. 4 (b) such as Fig. 4 (a)
With, considering that gap crystal bottom width is 3.15mm, the directions x are the length of crystal item, can use 15mm~20mm, with
X-axis is positive vertical and intersects at crystal array center p points.In Fig. 4, crystal probe module M, monocrystal c, crystal system
Number c0~c95；
For sensor M0 crystal inside array coding rules, as shown in Fig. 4 (b), c0~c95 is crystal System Number
0~95, crystal item number increases along y-axis (providing that this direction is row) positive and z-axis (providing that this direction is row) is positive, and first
The directions y directions z again, similar with M0, M1~M7 crystal inside items number is in terms of z-axis negative direction first in column direction to pass clockwise
Increase, then increases toward z-axis forward direction.
(4) detector rotation parameter is solved according to crystal probe size, obtains single rotation angle i.e. stepping angle.PET
Detector model is made of 4 couples of crystal probe M, and detection aperture D is 200mm, and crystal width w is 3.15mm, if wishing to ring
Answer the incompleteness of line smaller, then β should be obtained smaller, and this example β takes 2, can according to the calculation formula of step (4) in embodiment 1
It acquiresThen single rotation angle θnα °=3.6 °=β.
The total rotating acquisition number of detection ring is calculated according to crystal probe logarithm and stepping angle.It is walked according in embodiment 1
Suddenly the calculation formula (3) of (4) acquires number of revolutionsIt is not difficult to obtain such relationship：Number of sensors, it is single
Secondary rotation angle is inversely proportional with rotating acquisition number.
Step (5) is same as Example 1.
(6) time slice T in the present embodimentsliceStill it is set as 10s, it is known that number of revolutions 12 is fallen into a trap according to embodiment 1
Formula (4) can obtain TtotalFor 120s, simulation time relevant parameter is configured accordingly.
(7) simulation stage detector initial position is determined, according to TsliceIt is scanned, completes after the position is scanned to detection
Device is rotated and completes to scan next time, stops emulation until overall sweep time reaches setting value.As shown in Fig. 2 (b), imitate
True mode original state is that x-axis passes through the centers sensor M0, when with TsliceAfter completing initial position scanning, in next timeslice
Initial period is rotated temporally piece T for 3.6 ° with rotation angle to modelsliceIt is scanned, is and then rotated with same step
Until overall sweep time stops emulation when reaching the 120s of preset value.
Step (8), (9) will use identical virtual annulus number in embodiment 1 to carry out angle to the collected data that meet
Degree is corrected, and difference is only that number of revolutions increases, and the crystal number of overlapping becomes 2 from 3, still uses similar in Fig. 6 reflect
The relationship of penetrating is modified, wherein the c2 of R0_M0, and the c0 of R1_M0 after c3 will be rotated with single, c1 are overlapped, and are needed VR1_M0's
Crystal number is modified to c2, c3, c4, c5, c6, c7, remaining crystal item also carries out similar amendment.
Claims (8)
1. a kind of GATE emulation modes of non-complete ring-type PET rotation sweep patterns, which is characterized in that the emulation mode step is such as
Under:
(1) simulated environment is GATE platforms under Linux Ubuntu systems, establishes pet detector model on the platform, is detected
There are many selections for sensor logarithm in device, and in cyclic annular arrangement at equal intervals；
(2) the two dimensional gray slice map of the radiation source model of preset number is merged to the radiation source model that is converted to GATE and can be read
Format adjusts moulded dimension, according to the distribution of the material of pixel grey scale and GATE material data lab setting radioactive source models and nucleic
The configuration documentation of activity distribution；
(3) crystal probe parameter is set, determines sensor internal crystal coding rule；
(4) PET is solved according to crystal probe size and detects ring rotation parameter, selection monocrystal item is required according to imaging resolution
The integral multiple of corresponding central angle is calculated as rotation step angle by sensor logarithm and the rotation step angle of setting
The number of revolutions of whole circumference scanning is completed to detection ring；
(5) analog to digital conversion circuit parameter, particle interaction, output file format in pet detector are configured；
(6) detector initial sweep time, time slice and overall sweep time are set；
(7) simulation stage determines detector initial position, and is scanned according to time slice, completes after the position is scanned to non-
Complete detector rings carry out rotation and are scanned next time, stop emulation when overall sweep time reaches preset value；
(8) after completing entire scanning process, angle modification is carried out to the event data that meets of output, after amendment and is used for data
Recombination, the sinogram after recombination are used for image reconstruction.
2. a kind of GATE emulation modes of non-complete ring-type PET rotation sweep patterns according to claim 1, feature exist
In in step (1), establishing three-dimensional system of coordinate in GATE emulation platforms, build pet detector model, which includes
Ring T, n crystal probe M0~Mn are detected, and sensor is spacedly distributed on the detection ring T that internal diameter is D millimeters, and ring
Place face is parallel with xOy planes, wherein M0 numbers for home sensor, and in terms of z-axis negative direction, sensor number is by M0~Mn
It is incremented by clockwise.
3. a kind of GATE emulation modes of non-complete ring-type PET rotation sweep patterns according to claim 1, feature exist
In, in step (2), by the two dimension slicing figure for radiating source model merge and be converted to " .img ", " .hdr " format, pass through pixel ruler
Radioactive source moulded dimension is arranged in very little adjusting, the different material in different pixel grey scale region setting models, different radiation
Source radiation activity, final GATE platforms can be directly read with the radiation source model of " .img ", the storage of " .hdr " format.
4. a kind of GATE emulation modes of non-complete ring-type PET rotation sweep patterns according to claim 1, feature exist
In, in step (3), setting sensor parameters, including sensor internal crystal number, array way and size, crystal sensing
Device M0~Mn local arrays are L*H, and wherein the directions y number is L, and the directions z number is H, and crystal bottom width is w millimeters, x
Direction is the length of crystal item, and crystal item is incremented by again the incremental sequence in the directions z with the first directions y and is numbered, and wherein y-axis forward direction is
Row, z-axis forward direction are row.
5. a kind of GATE emulation modes of non-complete ring-type PET rotation sweep patterns according to claim 1, feature exist
In in step (4), rotation sweep stepping angle θ is arranged in the integral multiple β that central angle alpha is corresponded to according to monocrystal itemn, pass through setting
Stepping angle θnWith crystal probe logarithm n, the number of revolutions M that detection ring completes whole circumference scanning, computational methods are calculated
It is as follows：It is D millimeters to detect ring diameter, and single crystal width is w millimeters, and by approximation method, it is corresponding can to obtain single crystal item
Central angleSingle rotation angle θn=β α °, number of revolutions
6. a kind of GATE emulation modes of non-complete ring-type PET rotation sweep patterns according to claim 1, feature exist
In, in step (6), (7), setting emulation initial sweep time, time slice, overall sweep time and detector preliminary sweep
Position, emulation initial time are set as Tstart, sweep time T of the time slice by detector in each positionsliceIt determines, always
Swept-volume time TtotalThe number M rotated by detector is determined, meets Ttotal=M*Tslice, at the beginning of determining simulation stage detector
Beginning position R0, according to TsliceIt is scanned, detector is rotated after completing position scanning and completes to scan next time, weight
Multiple scanning stops emulation until overall sweep time reaches setting value.
7. a kind of GATE emulation modes of non-complete ring-type PET rotation sweep patterns according to claim 1, feature exist
In in step (8), after completing scanning, obtaining storage and meet the document of event data, by establishing a virtual complete circle
Ring, internal diameter are D millimeters, are connected between crystal item on annulus and size is also consistent with crystal size on non-complete detector rings T,
Its initial position crystal item number is 0, is incrementally virtually numbered by number clockwise in terms of z-axis negative sense.
8. a kind of GATE emulation modes of non-complete ring-type PET rotation sweep patterns according to claim 7, feature exist
In in step (8), meeting in event data for output includes the energy for receiving γ photons, time, crystal item number letter
Breath, line of response are two crystal lines for being nearly simultaneously received γ photons, by detecting the total angle and crystalline substance that ring rotates through
Body number information combines, and is mapped with virtual ring number, by crystal position and crystal item overlaps on virtual ring number
The number being modified on virtual ring, the event data whole progress crystal item number that meets that detector is rotated to M acquisition are corrected,
Revised number can reflect position of the line of response on virtual annulus, then to meeting event data using single slice of data
Reassembly algorithm is stored in sinogram.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810531389.0A CN108763758B (en) | 2018-05-29 | 2018-05-29 | GATE simulation method of incomplete annular PET (positron emission tomography) rotary scanning mode |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810531389.0A CN108763758B (en) | 2018-05-29 | 2018-05-29 | GATE simulation method of incomplete annular PET (positron emission tomography) rotary scanning mode |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108763758A true CN108763758A (en) | 2018-11-06 |
| CN108763758B CN108763758B (en) | 2022-05-03 |
Family
ID=64003505
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810531389.0A Active CN108763758B (en) | 2018-05-29 | 2018-05-29 | GATE simulation method of incomplete annular PET (positron emission tomography) rotary scanning mode |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108763758B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100074399A1 (en) * | 2008-09-19 | 2010-03-25 | Jefferson Science Associates, Llc | High resolution PET breast imager with improved detection efficiency |
| FR2986080A1 (en) * | 2012-01-23 | 2013-07-26 | Imacisio | Method for calculating position of interaction of gamma photon with scintillator crystal in sensor of positron emission tomography system in medical field, involves calculating interaction position from weighted average of estimates |
| CN104116519A (en) * | 2014-07-30 | 2014-10-29 | 西安电子科技大学 | Rotating plate double-panel PET system based on polygon view symmetry and imaging method thereof |
| JP2015072201A (en) * | 2013-10-03 | 2015-04-16 | 日立アロカメディカル株式会社 | Semiconductor radiation detector, nuclear medicine diagnostic device using semiconductor radiation detector, and semiconductor radiation detector manufacturing method |
| CN106057050A (en) * | 2016-08-16 | 2016-10-26 | 东北大学 | Fluorescent CT system simulation method based on GATE platform |
| CN106137237A (en) * | 2016-08-24 | 2016-11-23 | 沈阳东软医疗系统有限公司 | Portable cord source apparatus and plane source, the bearing calibration of PET device, device |
| CN205880233U (en) * | 2016-02-02 | 2017-01-11 | 武汉数字派特科技有限公司 | Variable structure PET equipment |
| CN107693037A (en) * | 2013-09-18 | 2018-02-16 | 上海联影医疗科技有限公司 | PET scan device and its time migration bearing calibration |
-
2018
- 2018-05-29 CN CN201810531389.0A patent/CN108763758B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100074399A1 (en) * | 2008-09-19 | 2010-03-25 | Jefferson Science Associates, Llc | High resolution PET breast imager with improved detection efficiency |
| FR2986080A1 (en) * | 2012-01-23 | 2013-07-26 | Imacisio | Method for calculating position of interaction of gamma photon with scintillator crystal in sensor of positron emission tomography system in medical field, involves calculating interaction position from weighted average of estimates |
| CN107693037A (en) * | 2013-09-18 | 2018-02-16 | 上海联影医疗科技有限公司 | PET scan device and its time migration bearing calibration |
| JP2015072201A (en) * | 2013-10-03 | 2015-04-16 | 日立アロカメディカル株式会社 | Semiconductor radiation detector, nuclear medicine diagnostic device using semiconductor radiation detector, and semiconductor radiation detector manufacturing method |
| CN104116519A (en) * | 2014-07-30 | 2014-10-29 | 西安电子科技大学 | Rotating plate double-panel PET system based on polygon view symmetry and imaging method thereof |
| CN205880233U (en) * | 2016-02-02 | 2017-01-11 | 武汉数字派特科技有限公司 | Variable structure PET equipment |
| CN106057050A (en) * | 2016-08-16 | 2016-10-26 | 东北大学 | Fluorescent CT system simulation method based on GATE platform |
| CN106137237A (en) * | 2016-08-24 | 2016-11-23 | 沈阳东软医疗系统有限公司 | Portable cord source apparatus and plane source, the bearing calibration of PET device, device |
Non-Patent Citations (4)
| Title |
|---|
| G. PANAYIOTAKIS 等: "Tomographic evaluation of a dual head PET", 《2010 IEEE INTERNATIONAL CONFERENCE ON IMAGING SYSTEMS AND TECHNIQUES》 * |
| WU HEYU 等: "A Feasibility Study of a Prototype PET Insert Device to Convert a General-Purpose Animal PET Scanner to Higher Resolution", 《JOURNAL OF NUCLEAR MEDICINE》 * |
| 刘豪佳 等: "基于GATE的PET系统仿真研究", 《核电子学与探测技术》 * |
| 徐晶华 等: "基于GATE的正电子发射断层成像仿真研究", 《CT理论与应用研究》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108763758B (en) | 2022-05-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN100561520C (en) | Rebuild the method and apparatus of (n+1) dimension image function from the Radon data | |
| CN100502785C (en) | Tomography device and method | |
| CN103908280B (en) | The method of PET scatter corrections | |
| Bequé et al. | Optimization of geometrical calibration in pinhole SPECT | |
| CN105682559B (en) | For estimating the method and system of scattering | |
| CN100587391C (en) | Projection rotating center measurement method suitable for 2D-CT scanning system | |
| Stevens et al. | Alignment of a volumetric tomography system | |
| CN109754445A (en) | Method and apparatus for generating X-ray tomographic image data | |
| CN102789510B (en) | Method for acquiring geometric correction parameter of PET (positron emission tomography) system | |
| CN107837090A (en) | Scattered ray correction based on sinogram in computer tomography | |
| CN101673413B (en) | System and method for modeling three-dimensional configuration of plant branches based on surface data | |
| CN108511043B (en) | X-CT virtual data acquisition and image reconstruction method and system based on numerical simulation | |
| US9858690B2 (en) | Computed tomography (CT) image reconstruction method | |
| CN109561865A (en) | It is a kind of that X-ray is overlapped to rebuild the method and system of 3-D image by space-time | |
| CN103996213A (en) | PET image rebuilding method and system | |
| Gillam et al. | Monte-Carlo simulations and image reconstruction for novel imaging scenarios in emission tomography | |
| CN107796835B (en) | X-ray cylindrical surface three-dimensional cone beam computer tomography method and device | |
| CN105894525A (en) | Detector correction method for positive electron emission tomography system | |
| CN105894550A (en) | Method for synchronously reconstructing dynamic PET image and tracer kinetic parameter on the basis of TV and sparse constraint | |
| CN101825433B (en) | Measuring method of offset of rotating center of rotating table of fan beam 2D-CT scanning system | |
| CN108333197B (en) | Bias industrial CT system turntable rotation center scaling method under scan pattern | |
| CN105832358A (en) | Imaging method based on system calibration for rotating double-panel PET system | |
| CN108763758A (en) | A kind of GATE emulation modes of non-complete ring-type PET rotation sweep patterns | |
| CN110215223A (en) | Scatter correction method, system, readable storage medium storing program for executing and equipment | |
| CN113902823B (en) | Projection method and system for PET image reconstruction |
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 |