CN109613593A - The gammaphoton active position methods of sampling, device, electronic equipment and storage medium - Google Patents
The gammaphoton active position methods of sampling, device, electronic equipment and storage medium Download PDFInfo
- Publication number
- CN109613593A CN109613593A CN201811417389.4A CN201811417389A CN109613593A CN 109613593 A CN109613593 A CN 109613593A CN 201811417389 A CN201811417389 A CN 201811417389A CN 109613593 A CN109613593 A CN 109613593A
- Authority
- CN
- China
- Prior art keywords
- active position
- crystal
- gammaphoton
- target
- effect
- 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
- 238000005070 sampling Methods 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000013078 crystal Substances 0.000 claims abstract description 317
- 230000000694 effects Effects 0.000 claims abstract description 93
- 239000000463 material Substances 0.000 claims description 24
- 238000004891 communication Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 5
- 239000004744 fabric Substances 0.000 claims description 3
- 230000006870 function Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 6
- 238000004590 computer program Methods 0.000 description 5
- 239000002356 single layer Substances 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 238000002600 positron emission tomography Methods 0.000 description 3
- 239000008710 crystal-8 Substances 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 210000003733 optic disk Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/20—Measuring radiation intensity with scintillation detectors
- G01T1/202—Measuring radiation intensity with scintillation detectors the detector being a crystal
Landscapes
- Chemical & Material Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of Radiation (AREA)
Abstract
The application proposes the gammaphoton active position methods of sampling, device, electronic equipment and storage medium in a kind of scintillation crystal detectors.Wherein, this method comprises: determining the crystal identification that gammaphoton is had an effect according to the output signal of scintillation crystal detectors, and obtain the probability distribution of crystals active position corresponding with crystal identification, and the probability distribution according to crystals active position, the position that gammaphoton is had an effect inside target crystal is sampled, to obtain the target effect site that gammaphoton is had an effect inside target crystal.Thus, probability distribution based on crystals active position, it is accurately sampled to the position that gammaphoton is had an effect inside target crystal, has accurately obtained the target effect site that crystals are had an effect, be the subsequent basis established based on target effect site accurate reconstruction image.
Description
Technical field
This application involves gammaphotons in terminal device technical field more particularly to a kind of scintillation crystal detectors to act on position
Set the gammaphoton active position methods of sampling, device, electronic equipment and storage medium in methods of sampling scintillation crystal detectors.
Background technique
Scintillation crystal detectors are widely used in gamma photon detection, the basic principle is that being stopped simultaneously using scintillation crystal
The energy of gammaphoton is converted into scintillation photons, detects gamma light indirectly using photon detector detection scintillation photons
Spatial position, time and the energy information of son effect.Scintillation crystal detectors are widely used in all kinds of imaging systems, especially
It is positron emission tomography (Positron Emission Tomography, PET) and single photon emission computed
Computed tomography (SPECT) system (Single Photon Emission Computed Tomography, SPECT).PET and SPECT pass through
Detection gammaphoton reconstructs the distribution of Radioactive isotope agent in vivo, to obtain the intracorporal metabolism of biology
Etc. information, be widely used in the diagnosis of medicine, drug research and clinical disease.
In the discrete crystal array detector design being widely adopted at present, detector only can determine that gammaphoton occurs
In which root scintillation crystal, and it can not determine the specific location that gammaphoton is had an effect inside the scintillation crystal.In image
During reconstruction, it is thus necessary to determine that the line of response or projection line of each gammaphoton event, traditional method are often simple
It is assumed that gammaphoton active position some determination point inside the central point or scintillation crystal of scintillation crystal incidence surface,
However, above-mentioned this mode will lead to, in the position for deviateing system imaging central region, there are parallax effects, so as to cause space
The decline of resolution ratio；Simultaneously this sparse discrete sampling will also result in image reconstruction process system transmission matrix calculate,
Line of response or projection line position calculate and project the inaccurate of back projection's process calculating process, to affect reconstruction image
The picture quality of spatial resolution and uniformity etc..Therefore, based on above-mentioned analysis, it can be seen that accurate to determine gamma light
The specific location that son is had an effect inside the scintillation crystal can have a great impact for subsequent image procossing.Therefore, such as
The problem of what accurate specific location for determining that gammaphoton is had an effect inside the scintillation crystal is current urgent need to resolve.
Summary of the invention
The application is intended to solve at least some of the technical problems in related technologies.For this purpose, the application
First purpose is to propose that the gammaphoton active position methods of sampling in a kind of scintillation crystal detectors, this method are based on crystal
The probability distribution of internal action position accurately samples to the position that gammaphoton is had an effect inside target crystal, quasi-
The target effect site that crystals are had an effect really has been obtained, has been established to be subsequent based on target effect site accurate reconstruction image
Basis.
Second purpose of the application is to propose gammaphoton active position sampling dress in a kind of scintillation crystal detectors
It sets.
The third purpose of the application is to propose a kind of electronic equipment.
The 4th purpose of the application is to propose a kind of storage medium.
The 5th purpose of the application is to propose a kind of computer program product.
To achieve the above object, the application first aspect embodiment proposes gammaphoton in a kind of scintillation crystal detectors
The active position methods of sampling, comprising: determine what gammaphoton was had an effect according to the output signal of the scintillation crystal detectors
Crystal identification；Obtain the probability distribution of crystals active position corresponding with the crystal identification；According to the crystals
The probability distribution of active position is sampled the position that gammaphoton is had an effect inside target crystal, to obtain gamma
The target effect site that photon is had an effect inside the target crystal, wherein the target crystal and the crystal identification
It is corresponding.
Optionally, further includes: obtain the probability value of each active position in the target crystal；According to each active position
Probability value, pre-establish the probability distribution of the crystals active position.
Optionally, the probability value for obtaining each active position in the target crystal, comprising: according to each effect position
The space coordinate of the spatial coordinated information, default incidence point set and each active position are respectively vertically mapped in target crystal table
The space coordinate of location point on face determines each active position respectively corresponding gammaphoton incident angle respectively, wherein
The location point and the default incidence point are in approximately the same plane；According to the corresponding gammaphoton of each active position
It is every inside the target crystal to determine that gammaphoton reaches for the crystal arrangement information of incident angle and the scintillation crystal detectors
The path length of the scintillator crystal materials passed through at a active position；Decaying according to gammaphoton in scintillator crystal materials is advised
Rule and each corresponding path length of active position, obtain the probability of each active position in the target crystal
Value.
Optionally, the probability value for obtaining each active position in the target crystal, comprising: S21, obtain the mesh
Mark the active position sequence in crystal, wherein the active position sequence includes multiple active positions arranged in sequence；
S22, for the current active position in the target crystal, according to the spatial coordinated information of current active position, default incidence
The space coordinate and current active position of point are vertically mapped in the space coordinate of the location point on target crystal surface, determine institute
State gammaphoton incident angle corresponding to current location；S23, according to the gammaphoton incident angle and the scintillation crystal
The crystal arrangement information of detector determines gammaphoton to reach inside the target crystal and is passed through at the current active position
Scintillator crystal materials path length；S24, according to gammaphoton scintillator crystal materials attenuation law and the path length
Degree, obtains the probability value of current active position described in the target crystal；S25, according to the active position sequence, obtain with
The adjacent next active position of the current active position, and using next active position as the current effect position
It sets, and repeats the step S22 to S24, until traversing all active positions in the active position sequence.
Optionally, the probability distribution according to the crystals active position, to gammaphoton in target crystal
The position that portion has an effect is sampled, to obtain the interacting goals position that gammaphoton is had an effect inside the target crystal
It sets, comprising: according to the probability distribution of sampling prescription and the crystals active position, to gammaphoton inside target crystal
The position having an effect carries out multiple sampling, and obtains multiple target effect sites according to sampling results.
To achieve the above object, the application second aspect embodiment proposes gammaphoton in a kind of scintillation crystal detectors
Active position sampling apparatus, comprising: determining module, for determining gamma light according to the output signal of the scintillation crystal detectors
The crystal identification that son is had an effect；First obtains module, acts on position for obtaining crystals corresponding with the crystal identification
The probability distribution set；Processing module, for the probability distribution according to the crystals active position, to gammaphoton in target
The position that crystals are had an effect is sampled, to obtain the target that gammaphoton is had an effect inside the target crystal
Active position, wherein the target crystal is corresponding with the crystal identification.
Optionally, second module is obtained, for obtaining the probability value of each active position in the target crystal；Establish mould
Block pre-establishes the probability distribution of the crystals active position for the probability value according to each active position.
Optionally, described second module is obtained, be specifically used for: according to the spatial coordinated information of each active position, preset
The space that the space coordinate of incidence point and each active position are respectively vertically mapped in the location point on target crystal surface is sat
Mark determines each active position respectively corresponding gammaphoton incident angle respectively, wherein the location point and described presets
Incidence point is in approximately the same plane；It is brilliant according to the corresponding gammaphoton incident angle of each active position and the flashing
The crystal arrangement information of bulk detector determines that gammaphoton reaches and is passed through at each active position inside the target crystal
The path length of scintillator crystal materials；According to gammaphoton scintillator crystal materials attenuation law and each active position
Corresponding path length obtains the probability value of each active position in the target crystal.
To achieve the above object, the application third aspect embodiment proposes a kind of electronic equipment, comprising: processor is deposited
Reservoir, communication interface and bus；The processor, the memory and the communication interface are connected and are completed by the bus
Mutual communication；The memory stores executable program code；The processor is stored by reading in the memory
Executable program code run program corresponding with the executable program code, for executing the application first aspect
The display methods at the application programs interface of embodiment.
To achieve the above object, the application fourth aspect embodiment proposes a kind of storage medium, wherein the storage is situated between
Matter is for storing application program, and the application program for executing gal in scintillation crystal detectors described herein at runtime
The horse photon active position methods of sampling.
To achieve the above object, the 5th aspect embodiment of the application proposes a kind of computer program product, wherein described
Computer program product is taken out for executing gammaphoton active position in scintillation crystal detectors described herein at runtime
Quadrat method.
In this application, crystal mark that gammaphoton is had an effect is determined according to the output signal of scintillation crystal detectors
Know, and obtains the probability distribution of crystals active position corresponding with crystal identification, and according to crystals active position
Probability distribution, the position that gammaphoton is had an effect inside target crystal is sampled, to obtain gammaphoton in mesh
The target effect site that mark crystals are had an effect.Probability distribution as a result, based on crystals active position, to gamma light
The position that son is had an effect inside target crystal is accurately sampled, and has accurately been obtained the target that crystals are had an effect and has been made
It is the subsequent basis established based on target effect site accurate reconstruction image with position.
Detailed description of the invention
Fig. 1 is the gammaphoton active position methods of sampling in the scintillation crystal detectors according to the application one embodiment
Flow chart；
Fig. 2 is the gammaphoton active position methods of sampling in the scintillation crystal detectors according to the application one embodiment
Flow chart；
Fig. 3 is the refined flow chart according to the application one embodiment step 201；
Fig. 4 is the refined flow chart according to the step 201 of the application another embodiment；
Fig. 5 is active position A in target crystal in the application one embodiment, default incidence point B, mapping point C and enters
Penetrate the schematic diagram of angle, θ；
Fig. 6 is gammaphoton active position sampling apparatus in the scintillation crystal detectors according to the application one embodiment
Structural schematic diagram；
Fig. 7 is gammaphoton active position sampling apparatus in the scintillation crystal detectors according to another embodiment of the application
Structural schematic diagram；
Fig. 8 is the structural schematic diagram according to the electronic equipment of the application one embodiment.
Specific embodiment
Embodiments herein is described below in detail, examples of the embodiments are shown in the accompanying drawings, wherein from beginning to end
Same or similar label indicates same or similar element or element with the same or similar functions.Below with reference to attached
The embodiment of figure description is exemplary, it is intended to for explaining the application, and should not be understood as the limitation to the application.
Below with reference to the accompanying drawings gammaphoton active position sampling side in the scintillation crystal detectors of the embodiment of the present application is described
Method, device, electronic equipment and storage medium.
Fig. 1 is the gammaphoton active position methods of sampling in the scintillation crystal detectors according to the application one embodiment
Flow chart.
As shown in Figure 1, the gammaphoton active position methods of sampling includes: in the scintillation crystal detectors
S101 determines the crystal identification that gammaphoton is had an effect according to the output signal of scintillation crystal detectors.
Specifically, the gammaphoton active position methods of sampling in scintillation crystal detectors provided by the embodiments of the present application, it can
To be executed by gammaphoton active position sampling apparatus in scintillation crystal detectors provided by the embodiments of the present application, which can be with
It is configured in the imaging system with scintillation crystal detectors, for example, single photon emission computerized tomography system is (i.e.
SPECT), positron emission tomography (i.e. PET) etc..
Wherein, crystal identification unique identification crystal can be indicated by a variety of expression crystal, for example, can be compiled with user's crystal
Number indicate crystal identification, in another example, can be used letter indicate crystal identification.
Wherein, which indicates crystal identification with crystal number.
Wherein, it should be noted that the scintillation crystal detectors of the embodiment can be single layer scintillation crystal detectors,
It can be the double-deck scintillation crystal detectors.
Wherein, single layer scintillation crystal detectors are that the flickering position in detector is arranged to make up by single layer scintillation crystal.
Wherein, the double-deck scintillation crystal detectors are that the flickering position in detector is arranged to make up by the double-deck scintillation crystal.
Wherein, which is described so that scintillation crystal detectors are single layer scintillation crystal detectors as an example.
As a kind of possible implementation, in the output signal for obtaining scintillation crystal detectors, can according to output signal,
The co-ordinate position information of pixel hit is determined, then, according to the region for the pixel that each crystal prestored is hit
Information is compared, and according to comparison result, determines the crystal identification that gammaphoton is had an effect.
S102 obtains the probability distribution of crystals active position corresponding with crystal identification.
It specifically, can be according to the probability of the crystal identification and crystals active position that prestore after obtaining crystal identification
The corresponding relationship of distribution obtains the probability distribution of the corresponding crystals active position of the crystal identification.
S103 has an effect inside target crystal to gammaphoton according to the probability distribution of crystals active position
Position be sampled, to obtain the target effect site that gammaphoton is had an effect inside target crystal, wherein target is brilliant
Body is corresponding with crystal identification.
Wherein, active position is the position that gammaphoton is had an effect in target crystal.
The embodiment of property as an example can obtain in order to subsequent according to the active position reconstruction image extracted
The probability distribution of crystals active position corresponding with crystal identification is taken, it can be according to the probability of the crystals active position point
Cloth carries out random sampling to the position that gammaphoton is had an effect inside target crystal.
Calculating and the image reconstruction of system transmission matrix are carried out in order to be more accurate, as an example the embodiment party of property
Formula, the position that can be had an effect inside target crystal according to the probability distribution of the crystals active position to gammaphoton into
The multiple random sampling of row, and multiple target effect sites are obtained according to multiple random sampling result.Then, it is taken out by repeatedly random
Multiple target effect sites of sample result carry out calculating and the image reconstruction process of subsequent transmission matrix.That is, should
Target effect site in embodiment can be multiple.
Wherein, it should be noted that for each active position in target crystal, the probability of active position is not
It together, can be that each active position be set according to the probability of active position in order to improve the accuracy of extracted target effect site
Set Sampling weights, wherein the probability of active position is bigger, and corresponding Sampling weights are bigger.
The gammaphoton active position methods of sampling, is visited according to scintillation crystal in the scintillation crystal detectors of the embodiment of the present application
The output signal for surveying device determines the crystal identification that gammaphoton is had an effect, and obtains crystals corresponding with crystal identification and make
With the probability distribution of position, and according to the probability distribution of crystals active position, to gammaphoton inside target crystal
The position having an effect is sampled, to obtain the target effect site that gammaphoton is had an effect inside target crystal.By
This carries out the position that gammaphoton is had an effect inside target crystal based on the probability distribution of crystals active position
Accurate sampling, has accurately obtained the target effect site that crystals are had an effect, and is subsequent accurate based on target effect site
The basis that reconstruction image is established.
On the basis of based on the above embodiment, in order to can be according to the probability distribution of crystals active position, to gamma light
The position that son is had an effect inside target crystal is sampled, and is determining gamma according to the output signal of scintillation crystal detectors
Before the crystal identification that photon is had an effect, as shown in Fig. 2, this method may include:
Step 201, the probability value of each active position in target crystal is obtained.
Wherein, it should be noted that according to the difference of application scenarios, can obtain in several ways each in target crystal
The probability value of active position, is illustrated below:
The first example
In this example, can according to the device identification of scintillation crystal detectors and the corresponding crystal identification of target crystal,
From the probability value for obtaining each active position in target crystal in server.
After the probability value for obtaining each active position, place can also be normalized to the probability value of each active position
Reason, so that the probability summation of each active position of crystals is 1, and according to the normalization probability of each active position
Value, determines the probability distribution of the active position of crystals.
Second of example
In this example, as shown in figure 3, above-mentioned steps 201, may include:
Step 211, according to the spatial coordinated information of each active position, the space coordinate of default incidence point and each work
It is respectively vertically mapped in the space coordinate of the location point on target crystal surface with position, determines each active position respectively respectively
Corresponding gammaphoton incident angle, wherein location point and default incidence point are in approximately the same plane.
Step 212, according to the corresponding gammaphoton incident angle of each active position and scintillation crystal detectors
Crystal arrangement information determines that gammaphoton reaches the scintillator crystal materials passed through at each active position inside target crystal
Path length.
Wherein, the crystal arrangement information of scintillation crystal detectors can include but is not limited to crystal arrangement mode, crystal
The information such as length, width and height.
Step 213, corresponding in the attenuation law of scintillator crystal materials and each active position according to gammaphoton
Path length obtains the probability value of each active position in target crystal.
The embodiment of property as an example, according to gammaphoton scintillator crystal materials attenuation law index of coincidence
Decaying.
Wherein, the attenuation law according to gammaphoton in scintillator crystal materials and the corresponding path of each active position
Length, the formula for calculating the probability value of each active position in target crystal are as follows:
P (A)=e-μD (1)
Wherein, the P in formula (A) indicates the probability of crystals active position A, and μ is preparatory according to scintillator crystal materials
Determining attenuation parameter, D indicate path length.
Step 202, according to the probability value of each active position, the probability distribution of crystals active position is pre-established.
The embodiment of property as an example can also be to each work after the probability value for obtaining each active position
It is normalized with the probability value of position, so that the probability summation of each active position of crystals is 1, and according to
The normalization probability value of each active position, determines the probability distribution of the active position of crystals.
The third example
In this example, as shown in figure 4, above-mentioned steps 201 may include:
S21, active position sequence in target crystal is obtained, wherein active position sequence includes multiple arranging in sequence
The active position of column；
S22, for the current active position in target crystal, according to the spatial coordinated information of current active position, default
The space coordinate of incidence point and current active position are vertically mapped in the space coordinate of the location point on target crystal surface, really
Gammaphoton incident angle corresponding to settled front position.
S23, according to the crystal arrangement information of gammaphoton incident angle and scintillation crystal detectors, determine that gammaphoton arrives
The path length of the scintillator crystal materials passed through at current active position inside up to target crystal.
S24, according to gammaphoton in the attenuation law and path length of scintillator crystal materials, obtain in target crystal current
The probability value of active position.
S25 obtains the next active position adjacent with current active position according to active position sequence, and will be next
A active position repeats step S22 to S24 as current active position, until the institute in ergodic action position sequence
There is active position.
Specifically, according to active position sequence, obtain the next active position adjacent with current active position, and will under
One active position judges whether current active position is the last one work in active position sequence as current active position
With position, if so, determining all active positions traversed in active position sequence.If it is not, then repeating step S22
To S23.
The probability value for determining active position A in target crystal is described below with reference to Fig. 2, it is assumed that according to scintillation crystal
The crystal number that the output signal of detector determines that gammaphoton is had an effect is crystal 8.
Firstly, the information such as space coordinate of crystal 8 being had an effect according to gammaphoton, calculate gammaphoton entrance crystal
Angle, θ；
Specifically, it is sat according to the space of the space coordinate of default incidence point B, the space coordinate of active position A and C point
Mark.
Wherein, C point is that active position A is vertically mapped to corresponding location point on 8 surface of crystal.
Wherein, it should be noted that location point C and default incidence point A are in approximately the same plane.
Then, it according to crystal arrangement mode in gammaphoton entrance crystal angle and detector, calculates gammaphoton and reaches
The path length D for the scintillator crystal materials that the place crystals active position A (x, y, z) passes through；
Later, according to above-mentioned formula (1), the probability value of active position A is calculated.
Wherein, it is to be understood that, can be by identical with active position A for other active positions in the crystal
Mode calculates the probability value of other active positions, which repeats no more.
In order to realize above-described embodiment, the present invention also proposes that gammaphoton active position is taken out in a kind of scintillation crystal detectors
Sampling device.
Fig. 6 is gammaphoton active position sampling apparatus in scintillation crystal detectors according to an embodiment of the invention
Structural schematic diagram.
As shown in fig. 6, in the scintillation crystal detectors gammaphoton active position sampling apparatus include determining module 110,
First obtains module 120 and processing module 130, in which:
Determining module 110 determines the crystalline substance that gammaphoton is had an effect for the output signal according to scintillation crystal detectors
Body mark.
First obtains module 120, for obtaining the probability distribution of crystals active position corresponding with crystal identification.
Processing module 130, for the probability distribution according to crystals active position, to gammaphoton in target crystal
The position that portion has an effect is sampled, to obtain the target effect site that gammaphoton is had an effect inside target crystal,
Wherein, target crystal is corresponding with crystal identification.
Wherein, it should be noted that aforementioned that the gammaphoton active position methods of sampling in scintillation crystal detectors is implemented
The explanation of example is also applied for gammaphoton active position sampling apparatus in the scintillation crystal detectors of the embodiment, herein not
It repeats again.
Gammaphoton active position sampling apparatus, is visited according to scintillation crystal in the scintillation crystal detectors of the embodiment of the present application
The output signal for surveying device determines the crystal identification that gammaphoton is had an effect, and obtains crystals corresponding with crystal identification and make
With the probability distribution of position, and according to the probability distribution of crystals active position, to gammaphoton inside target crystal
The position having an effect is sampled, to obtain the target effect site that gammaphoton is had an effect inside target crystal.By
This carries out the position that gammaphoton is had an effect inside target crystal based on the probability distribution of crystals active position
Accurate sampling, has accurately obtained the target effect site that crystals are had an effect, and is subsequent accurate based on target effect site
The basis that reconstruction image is established.
In one embodiment of the application, on the basis of shown in Fig. 6, as shown in fig. 7, the device can also include:
Second obtains module 140, for obtaining the probability value of each active position in target crystal.
Module 150 is established, for the probability value according to each active position, pre-establishes the general of crystals active position
Rate distribution.
Wherein, it should be noted that according to the difference of application scenarios, can obtain in several ways each in target crystal
The probability value of active position, is illustrated below:
The first example
In this example, second obtain module 140 can device identification according to scintillation crystal detectors and target crystal
Corresponding crystal identification, from the probability value for obtaining each active position in target crystal in server.
After the second acquisition module 140 obtains the probability value of each active position, establishing module can be to each active position
Probability value be normalized so that the probability summation of each active position of crystals is 1, and according to each work
With the normalization probability value of position, the probability distribution of the active position of crystals is determined.
Second of example
In this example, second module is obtained, be specifically used for: according to the spatial coordinated information of each active position, preset
The space that the space coordinate of incidence point and each active position are respectively vertically mapped in the location point on target crystal surface is sat
Mark determines each active position respectively corresponding gammaphoton incident angle respectively, wherein location point at default incidence point
In in approximately the same plane；According to the crystalline substance of each active position corresponding gammaphoton incident angle and scintillation crystal detectors
Body arrangement information determines that gammaphoton reaches the road of the scintillator crystal materials passed through at each active position inside target crystal
Electrical path length；According to gammaphoton scintillator crystal materials attenuation law and the corresponding path length of each active position,
Obtain the probability value of each active position in target crystal.
Wherein, second the detailed process that module 140 obtains the probability value of each active position in target crystal is obtained, such as schemed
Shown in 3, details are not described herein again.
Second of example, the second acquisition module 140 can obtain each effect in target crystal by mode as shown in Figure 4
The probability value of position.
The detailed process for obtaining the probability value of each active position in 140 target crystal of module about second, reference can be made on
The corresponding portion of method description is stated, details are not described herein again.
In order to realize above-described embodiment, the application also proposes a kind of electronic equipment.
Fig. 8 is the structural schematic diagram according to the electronic equipment of the application one embodiment.
As shown in figure 8, the electronic equipment includes processor 81, memory 82, communication interface 83 and bus 84, in which:
Processor 81, memory 82 and communication interface 83 connect by bus 84 and complete mutual communication；Memory
82 storage executable program codes；Processor 81 is run by reading the executable program code stored in memory 82 and can
The corresponding program of program code is executed, for executing gammaphoton active position in the scintillation crystal detectors in upper embodiment
The methods of sampling.
The electronic equipment of the embodiment of the present application determines that gammaphoton is made according to the output signal of scintillation crystal detectors
Crystal identification, and the probability distribution of crystals active position corresponding with crystal identification is obtained, and according in crystal
The probability distribution of portion's active position is sampled the position that gammaphoton is had an effect inside target crystal, to obtain gal
The target effect site that horse photon is had an effect inside target crystal.Probability point as a result, based on crystals active position
Cloth accurately samples to the position that gammaphoton is had an effect inside target crystal, has accurately obtained crystals generation
The target effect site of effect is the subsequent basis established based on target effect site accurate reconstruction image.
In order to realize above-described embodiment, the application also proposes a kind of imaging system, which includes above-described embodiment
Scintillation crystal detectors in gammaphoton active position sampling apparatus.
Wherein, imaging system can be disconnected for single photon emission computerized tomography system (i.e. SPECT), positron emission
Layer imaging system (i.e. PET) etc..
In order to realize above-described embodiment, the application also proposes a kind of storage medium, wherein storage medium is for storing application
Program, application program is for gammaphoton active position sampling side in the scintillation crystal detectors of execution the application at runtime
Method.
In order to realize above-described embodiment, the application also proposes a kind of computer program product, wherein computer program product
The gammaphoton active position methods of sampling in scintillation crystal detectors for executing the application at runtime.
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show
The description of example " or " some examples " etc. means specific features, structure, material or spy described in conjunction with this embodiment or example
Point is contained at least one embodiment or example of the application.In the present specification, schematic expression of the above terms are not
It must be directed to identical embodiment or example.Moreover, particular features, structures, materials, or characteristics described can be in office
It can be combined in any suitable manner in one or more embodiment or examples.In addition, without conflicting with each other, the skill of this field
Art personnel can tie the feature of different embodiments or examples described in this specification and different embodiments or examples
It closes and combines.
In addition, term " first ", " second " are used for descriptive purposes only and cannot be understood as indicating or suggesting relative importance
Or implicitly indicate the quantity of indicated technical characteristic.Define " first " as a result, the feature of " second " can be expressed or
Implicitly include at least one this feature.In the description of the present application, the meaning of " plurality " is at least two, such as two, three
It is a etc., unless otherwise specifically defined.
Any process described otherwise above or method description are construed as in flow chart or herein, and expression includes
It is one or more for realizing specific logical function or process the step of executable instruction code module, segment or portion
Point, and the range of the preferred embodiment of the application includes other realization, wherein can not press shown or discussed suitable
Sequence, including according to related function by it is basic simultaneously in the way of or in the opposite order, Lai Zhihang function, this should be by the application
Embodiment person of ordinary skill in the field understood.
Expression or logic and/or step described otherwise above herein in flow charts, for example, being considered use
In the order list for the executable instruction for realizing logic function, may be embodied in any computer-readable medium, for
Instruction execution system, device or equipment (such as computer based system, including the system of processor or other can be held from instruction
The instruction fetch of row system, device or equipment and the system executed instruction) it uses, or combine these instruction execution systems, device or set
It is standby and use.For the purpose of this specification, " computer-readable medium ", which can be, any may include, stores, communicates, propagates or pass
Defeated program is for instruction execution system, device or equipment or the dress used in conjunction with these instruction execution systems, device or equipment
It sets.The more specific example (non-exhaustive list) of computer-readable medium include the following: there is the electricity of one or more wirings
Interconnecting piece (electronic device), portable computer diskette box (magnetic device), random access memory (RAM), read-only memory
(ROM), erasable edit read-only storage (EPROM or flash memory), fiber device and portable optic disk is read-only deposits
Reservoir (CDROM).In addition, computer-readable medium can even is that the paper that can print described program on it or other are suitable
Medium, because can then be edited, be interpreted or when necessary with it for example by carrying out optical scanner to paper or other media
His suitable method is handled electronically to obtain described program, is then stored in computer storage.
It should be appreciated that each section of the application can be realized with hardware, software, firmware or their combination.Above-mentioned
In embodiment, software that multiple steps or method can be executed in memory and by suitable instruction execution system with storage
Or firmware is realized.It, and in another embodiment, can be under well known in the art for example, if realized with hardware
Any one of column technology or their combination are realized: having a logic gates for realizing logic function to data-signal
Discrete logic, with suitable combinational logic gate circuit specific integrated circuit, programmable gate array (PGA), scene
Programmable gate array (FPGA) etc..
Those skilled in the art are understood that realize all or part of step that above-described embodiment method carries
It suddenly is that relevant hardware can be instructed to complete by program, the program can store in a kind of computer-readable storage medium
In matter, which when being executed, includes the steps that one or a combination set of embodiment of the method.
It, can also be in addition, can integrate in a processing module in each functional unit in each embodiment of the application
It is that each unit physically exists alone, can also be integrated in two or more units in a module.Above-mentioned integrated mould
Block both can take the form of hardware realization, can also be realized in the form of software function module.The integrated module is such as
Fruit is realized and when sold or used as an independent product in the form of software function module, also can store in a computer
In read/write memory medium.
Storage medium mentioned above can be read-only memory, disk or CD etc..Although having been shown and retouching above
Embodiments herein is stated, it is to be understood that above-described embodiment is exemplary, and should not be understood as the limit to the application
System, those skilled in the art can be changed above-described embodiment, modify, replace and become within the scope of application
Type.
Claims (10)
1. the gammaphoton active position methods of sampling in a kind of scintillation crystal detectors characterized by comprising
The crystal identification that gammaphoton is had an effect is determined according to the output signal of the scintillation crystal detectors；
Obtain the probability distribution of crystals active position corresponding with the crystal identification；
According to the probability distribution of the crystals active position, the position having an effect inside target crystal to gammaphoton
It is sampled, to obtain the target effect site that gammaphoton is had an effect inside the target crystal, wherein the target
Crystal is corresponding with the crystal identification.
2. the method as described in claim 1, which is characterized in that further include:
Obtain the probability value of each active position in the target crystal；
According to the probability value of each active position, the probability distribution of the crystals active position is pre-established.
3. method according to claim 2, which is characterized in that each active position is general in the acquisition target crystal
Rate value, comprising:
It is respectively hung down according to the spatial coordinated information of each active position, the space coordinate of default incidence point and each active position
Directly be mapped in the space coordinate of location point on target crystal surface, determine each active position respectively respectively corresponding to gamma
Photon incident angle, wherein the location point and the default incidence point are in approximately the same plane；
According to the crystal arrangement of each active position corresponding gammaphoton incident angle and the scintillation crystal detectors
Information determines that gammaphoton reaches the path of the scintillator crystal materials passed through at each active position inside the target crystal
Length；
Attenuation law and each corresponding path length of the active position according to gammaphoton in scintillator crystal materials,
Obtain the probability value of each active position in the target crystal.
4. method according to claim 2, which is characterized in that each active position is general in the acquisition target crystal
Rate value, comprising:
Active position sequence in S21, the acquisition target crystal, wherein the active position sequence includes multiple according to suitable
The active position of sequence arrangement；
S22, for the current active position in the target crystal, according to the spatial coordinated information of current active position, default
The space coordinate of incidence point and current active position are vertically mapped in the space coordinate of the location point on target crystal surface, really
Gammaphoton incident angle corresponding to the fixed current location；
S23, according to the crystal arrangement information of the gammaphoton incident angle and the scintillation crystal detectors, determine gamma light
Son reaches the path length of the scintillator crystal materials passed through at the current active position inside the target crystal；
S24, the attenuation law according to gammaphoton in scintillator crystal materials and the path length, obtain in the target crystal
The probability value of the current active position；
S25, according to the active position sequence, obtain the next active position adjacent with the current active position, and general
Next active position repeats the step S22 to S24 as the current active position, until traversal institute
State all active positions in active position sequence.
5. the method as described in claim 1, which is characterized in that the probability according to the crystals active position point
Cloth is sampled the position that gammaphoton is had an effect inside target crystal, brilliant in the target to obtain gammaphoton
The target effect site that internal portion has an effect, comprising:
According to the probability distribution of sampling prescription and the crystals active position, gammaphoton is occurred inside target crystal
The position of effect carries out multiple sampling, and obtains multiple target effect sites according to sampling results.
6. gammaphoton active position sampling apparatus in a kind of scintillation crystal detectors characterized by comprising
Determining module, for determining crystal mark that gammaphoton is had an effect according to the output signal of the scintillation crystal detectors
Know；
First obtains module, for obtaining the probability distribution of crystals active position corresponding with the crystal identification；
Processing module, for the probability distribution according to the crystals active position, to gammaphoton inside target crystal
The position having an effect is sampled, to obtain the interacting goals position that gammaphoton is had an effect inside the target crystal
It sets, wherein the target crystal is corresponding with the crystal identification.
7. device as claimed in claim 6, which is characterized in that further include:
Second obtains module, for obtaining the probability value of each active position in the target crystal；
Module is established, for the probability value according to each active position, pre-establishes the probability of the crystals active position
Distribution.
8. device as claimed in claim 7, which is characterized in that described second obtains module, is specifically used for:
It is respectively hung down according to the spatial coordinated information of each active position, the space coordinate of default incidence point and each active position
Directly be mapped in the space coordinate of location point on target crystal surface, determine each active position respectively respectively corresponding to gamma
Photon incident angle, wherein the location point and the default incidence point are in approximately the same plane；
According to the crystal arrangement of each active position corresponding gammaphoton incident angle and the scintillation crystal detectors
Information determines that gammaphoton reaches the path of the scintillator crystal materials passed through at each active position inside the target crystal
Length；
Attenuation law and each corresponding path length of the active position according to gammaphoton in scintillator crystal materials,
Obtain the probability value of each active position in the target crystal.
9. a kind of electronic equipment, which is characterized in that including reason device, memory, communication interface and bus；
The processor, the memory are connected by the bus with the communication interface and complete mutual communication；
The memory stores executable program code；
The processor is run and the executable program by reading the executable program code stored in the memory
The corresponding program of code, to make for executing gammaphoton in scintillation crystal detectors as described in any one in claim 1-5
With the position methods of sampling.
10. a kind of storage medium, which is characterized in that for storing application program, the application program is used for the storage medium
Gammaphoton active position sampling side in scintillation crystal detectors as described in any one in claim 1-5 is executed at runtime
Method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811417389.4A CN109613593B (en) | 2018-11-26 | 2018-11-26 | Gamma photon action position sampling method, gamma photon action position sampling device, electronic equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811417389.4A CN109613593B (en) | 2018-11-26 | 2018-11-26 | Gamma photon action position sampling method, gamma photon action position sampling device, electronic equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109613593A true CN109613593A (en) | 2019-04-12 |
| CN109613593B CN109613593B (en) | 2020-10-20 |
Family
ID=66005082
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811417389.4A Active CN109613593B (en) | 2018-11-26 | 2018-11-26 | Gamma photon action position sampling method, gamma photon action position sampling device, electronic equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109613593B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110664423A (en) * | 2019-09-12 | 2020-01-10 | 东软医疗系统股份有限公司 | Imaging method, imaging device, detector, terminal equipment and PET system |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2001071382A1 (en) * | 2000-03-23 | 2001-09-27 | Cti Pet Systems, Inc. | Pile-up reduction in multi-crystal gamma detectors |
| CN101163989A (en) * | 2005-04-22 | 2008-04-16 | 皇家飞利浦电子股份有限公司 | Pet/mr scanner with time-of-flight capability |
| CN102727237A (en) * | 2011-04-06 | 2012-10-17 | 株式会社东芝 | Medical image diagnostic apparatus and control method |
| US20130056639A1 (en) * | 2008-09-11 | 2013-03-07 | Siemens Medical Solutions Usa, Inc. | Timing Response Improvement in Light-Sharing Detectors |
| CN103547942A (en) * | 2012-05-16 | 2014-01-29 | 株式会社东芝 | Random event reduction method, random event reduction device, and non-temporary computer-readable recording medium |
| CN104155677A (en) * | 2014-07-21 | 2014-11-19 | 北京辛耕普华医疗科技有限公司 | Coincidence response line positioning method and device for PET detector |
| CN104337531A (en) * | 2013-07-25 | 2015-02-11 | 苏州瑞派宁科技有限公司 | Online energy conforming method and system for full-digital PET system |
| US20150057535A1 (en) * | 2013-05-08 | 2015-02-26 | Arkadiusz Sitek | Systems and methods for attenuation correction in time-of-flight positron emission tomography |
| CN104732541A (en) * | 2015-03-27 | 2015-06-24 | 北京永新医疗设备有限公司 | Method and device for acquiring PET crystal position lookup table |
| US20160170040A1 (en) * | 2014-12-15 | 2016-06-16 | General Electric Company | Systems and methods for scintillators having micro-crack surfaces |
| US9645258B1 (en) * | 2014-11-25 | 2017-05-09 | Sandia Corporation | Segmented scintillation antineutrino detector |
| CN107684436A (en) * | 2017-06-30 | 2018-02-13 | 上海联影医疗科技有限公司 | The method that the position of incident photon is determined in PET system |
-
2018
- 2018-11-26 CN CN201811417389.4A patent/CN109613593B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2001071382A1 (en) * | 2000-03-23 | 2001-09-27 | Cti Pet Systems, Inc. | Pile-up reduction in multi-crystal gamma detectors |
| CN101163989A (en) * | 2005-04-22 | 2008-04-16 | 皇家飞利浦电子股份有限公司 | Pet/mr scanner with time-of-flight capability |
| US20130056639A1 (en) * | 2008-09-11 | 2013-03-07 | Siemens Medical Solutions Usa, Inc. | Timing Response Improvement in Light-Sharing Detectors |
| CN102727237A (en) * | 2011-04-06 | 2012-10-17 | 株式会社东芝 | Medical image diagnostic apparatus and control method |
| CN103547942A (en) * | 2012-05-16 | 2014-01-29 | 株式会社东芝 | Random event reduction method, random event reduction device, and non-temporary computer-readable recording medium |
| US20150057535A1 (en) * | 2013-05-08 | 2015-02-26 | Arkadiusz Sitek | Systems and methods for attenuation correction in time-of-flight positron emission tomography |
| CN104337531A (en) * | 2013-07-25 | 2015-02-11 | 苏州瑞派宁科技有限公司 | Online energy conforming method and system for full-digital PET system |
| CN104155677A (en) * | 2014-07-21 | 2014-11-19 | 北京辛耕普华医疗科技有限公司 | Coincidence response line positioning method and device for PET detector |
| US9645258B1 (en) * | 2014-11-25 | 2017-05-09 | Sandia Corporation | Segmented scintillation antineutrino detector |
| US20160170040A1 (en) * | 2014-12-15 | 2016-06-16 | General Electric Company | Systems and methods for scintillators having micro-crack surfaces |
| CN104732541A (en) * | 2015-03-27 | 2015-06-24 | 北京永新医疗设备有限公司 | Method and device for acquiring PET crystal position lookup table |
| CN107684436A (en) * | 2017-06-30 | 2018-02-13 | 上海联影医疗科技有限公司 | The method that the position of incident photon is determined in PET system |
Non-Patent Citations (1)
| Title |
|---|
| 牛珂楠等: "LYSO闪烁晶体PET探测器设计及位置映射算法研究", 《北京生物医学工程》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110664423A (en) * | 2019-09-12 | 2020-01-10 | 东软医疗系统股份有限公司 | Imaging method, imaging device, detector, terminal equipment and PET system |
| CN110664423B (en) * | 2019-09-12 | 2023-05-09 | 沈阳智核医疗科技有限公司 | Imaging method, imaging device, detector, terminal equipment and PET system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109613593B (en) | 2020-10-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3793266B2 (en) | Positron CT apparatus and image reconstruction method thereof | |
| US5793045A (en) | Nuclear imaging using variable weighting | |
| CN105125231B (en) | A kind of minimizing technology and device of PET image ring artifact | |
| CN100506154C (en) | Nuclear medical apparatus | |
| CN101842806B (en) | Dirty isotope PET reconstruction | |
| CN101490581A (en) | Method and system for improved TOF PET reconstruction | |
| CN107976706B (en) | A kind of counting loss bearing calibration of PET system and device | |
| DE102010037037A1 (en) | Method for providing motion correction in tomographic imaging such as positron emission tomography (PET) imaging, involves combining registered gated emission tomography images to generate motion corrected images | |
| CN108523916A (en) | A kind of more gammaphotons meet imaging system and method | |
| CN111788499B (en) | Dead time correction method in quantitative Positron Emission Tomography (PET) reconstruction of various objects and radioactivity distributions | |
| JP2000271089A (en) | Quality control system for medical diagnostic device | |
| US10245002B2 (en) | Isotope specific calibration of a dose calibrator for quantitative functional imaging | |
| CN114666495A (en) | High-resolution Compton camera imaging method and device, electronic equipment and medium | |
| CN108209958B (en) | Method for determining and acquiring normalization correction factor and medical imaging method | |
| JP2012088300A (en) | Positron radiation tomography (pet) imaging system, method and nuclear medicine imaging system | |
| CN109077748A (en) | A kind of accurately PET normalization method | |
| CN108186043A (en) | The method and apparatus for correcting PET device sensitivity | |
| CN102487607B (en) | Nuclear medicine imaging device and control method | |
| US9910162B2 (en) | Calibrating in single photon emission computed tomography with multi-emission energies | |
| CN109613593A (en) | The gammaphoton active position methods of sampling, device, electronic equipment and storage medium | |
| US7564035B2 (en) | Generating detector efficiency estimates for a pet scanner | |
| CN109864754A (en) | The system and method being imaged for using improved dosage | |
| CN103842847B (en) | Nuclear medicine diagnostic apparatus, image diagnosing system and image processing method | |
| CN107456235A (en) | Method for correcting position and system | |
| CN109875592A (en) | A kind of methods, devices and systems that PET and SPECT is imaged simultaneously |
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 | ||
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: 610219 Chengdu Tianfu International Biological City, Chengdu, Sichuan Province (No. 618 Fenghuang Road, Shuangliu District) Patentee after: Chengdu Yongxin Medical Equipment Co.,Ltd. Country or region after: China Address before: 102206 1st floor, building 1, yard 20, shengshengyuan Road, Changping District, Beijing Patentee before: BEIJING NOVEL MEDICAL EQUIPMENT Ltd. Country or region before: China |