CN108459344A - Local spectral analysis method and device based on γ cameras - Google Patents
Local spectral analysis method and device based on γ cameras Download PDFInfo
- Publication number
- CN108459344A CN108459344A CN201810136044.5A CN201810136044A CN108459344A CN 108459344 A CN108459344 A CN 108459344A CN 201810136044 A CN201810136044 A CN 201810136044A CN 108459344 A CN108459344 A CN 108459344A
- Authority
- CN
- China
- Prior art keywords
- gamma ray
- visual field
- ray radiator
- photon
- full visual
- 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
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
Abstract
The present invention relates to radioactive source detection technique fields, more particularly to a kind of local spectral analysis method and device based on γ cameras, this method includes obtaining full visual field perspective view, and the photon that the full visual field perspective view is launched by gamma ray radiator is projeced into the detector across the encoding board collimator and is formed.The multiple gamma ray radiator regions for sending out photon are determined in the full visual field perspective view according to algorithm for reconstructing, and calculate weight factor of each described gamma ray radiator region relative to the full visual field perspective view.Each energy section photon counting of actual measurement corresponding to the pixel that the pixel value in each gamma ray radiator region is not zero is simulated to be multiplied with the weight factor, the Energy distribution situation of the photon in each gamma ray radiator region is obtained, the local that each gamma ray radiator is obtained according to the Energy distribution situation of the photon in each gamma ray radiator region is composed.By this programme by establishing local spectrum in self defined area, nucleic classification is analyzed, realizes nucleic early warning etc..
Description
Technical field
The present invention relates to radioactive source detection technique fields, in particular to a kind of local spectrum analysis based on γ cameras
Method and device.
Background technology
Coded imaging γ camera techniques, originate in the fifties in last century.The nineties in last century, encoding board collimate γ phases
Machine technology gradually develops into the known approaches of increased radioactivity imaging and video frequency tracking.Due to encoding board imaging technique overcome it is small
The shortcomings of borescopic imaging sensitivity is low, the parallel aperture collimator visual field is narrow is increasingly becoming the mainstream skill of radio nuclide source terms analysis and tracking
Art.Conventionally, as containing numerous mathematical operators or by iterating to calculate in increased radioactivity algorithm for reconstructing,
It seeks photon projected position and the position after final rebuild is more difficult, can only generally obtain full visual field power spectrum, analyze the full visual field
Interior possible source item is but difficult to realize the energy spectrum analysis of each gamma ray radiator region and hot spot and source item identification in the full visual field.Therefore,
It is very necessary to provide a kind of local spectral analysis method based on γ cameras.
Invention content
The purpose of the present invention is to provide a kind of local spectral analysis methods based on γ cameras, to realize in User Defined
The purpose of identification of nucleic classification and nucleic migrate early warning in region, to adapt to for different scenes.
Another object of the present invention is to provide a kind of local spectrum analysis device based on γ cameras, with realize user from
The purpose of identification of nucleic classification and nucleic migrate early warning in definition region, to adapt to be used for different scenes.
To achieve the goals above, technical solution used in the embodiment of the present invention is as follows：
In a first aspect, an embodiment of the present invention provides a kind of local spectral analysis method based on γ cameras, the γ cameras
Including encoding board collimator and detector, the method includes：
Full visual field perspective view is obtained, the full visual field perspective view is passed through by the photon that all gamma ray radiators are launched in the visual field
The encoding board collimator is projeced into the detector and is formed；
Full visual field source positions and intensity distribution are rebuild according to projection reconstruction algorithm, determine that the multiple γ for sending out photon are put
Penetrate source region；
Calculate each weight factor of the gamma ray radiator region relative to the full visual field；
The corresponding actual measurement of pixel that the pixel value of the independent projections in each gamma ray radiator region is not zero is each
Energy section photon counting is multiplied with the weight factor, obtains the Energy distribution feelings of the photon in each gamma ray radiator region
Condition；
Each described gamma ray radiator is obtained according to the Energy distribution situation of the photon in gamma ray radiator region described in each
The local in region is composed.
Second aspect, the embodiment of the present invention additionally provide a kind of local spectrum analysis device based on γ cameras, the γ phases
Machine includes encoding board collimator and detector, and described device includes：
Acquisition module, for obtaining full visual field perspective view, the full visual field perspective view is sent out by all gamma ray radiators in the visual field
The photon gone out is projeced into the detector across the encoding board collimator and is formed；
Area determination module, for rebuilding full visual field source positions according to projection reconstruction algorithm, determination sends out photon
Multiple gamma ray radiator regions；
First computing module, for calculating power of each described gamma ray radiator region relative to the full visual field perspective view
Repeated factor；
Second computing module, the picture for the pixel value of the independent projections in each gamma ray radiator region to be not zero
The corresponding each energy section photon counting of actual measurement of vegetarian refreshments is multiplied with the weight factor, obtains each gamma ray radiator region
The Energy distribution situation of photon；
Analysis module, for obtaining each according to the Energy distribution situation of the photon in gamma ray radiator region described in each
The local in the gamma ray radiator region is composed.
A kind of local spectral analysis method and device based on γ cameras provided in an embodiment of the present invention, this method include obtaining
Full visual field perspective view, the photon that the full visual field perspective view is launched by gamma ray radiator are projeced into across the encoding board collimator
The detector is formed.The multiple gamma ray radiator areas for sending out photon are determined in the full visual field perspective view according to algorithm for reconstructing
Domain, and calculate weight factor of each described gamma ray radiator region relative to the full visual field perspective view.Simulate each γ
Each energy section photon counting of actual measurement corresponding to the pixel that is not zero of pixel value of radiation source region and the weight factor phase
Multiply, the Energy distribution situation of the photon in each gamma ray radiator region is obtained, according to the photon in each gamma ray radiator region
Energy distribution situation obtains the local spectrum of each gamma ray radiator.By searching gamma ray radiator region in full visual field perspective view,
With respectively to each gamma ray radiator regional analysis, and local is composed and is compared with full visual field power spectrum, and then realization is being made by oneself
Nucleic classification is analyzed in adopted region, realizes nucleic early warning etc..
To enable the above objects, features and advantages of the present invention to be clearer and more comprehensible, preferred embodiment cited below particularly, and coordinate
Appended attached drawing, is described in detail below.
Description of the drawings
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to needed in the embodiment attached
Figure is briefly described, it should be understood that the following drawings illustrates only certain embodiments of the present invention, therefore is not construed as pair
The restriction of range for those of ordinary skill in the art without creative efforts, can also be according to this
A little attached drawings obtain other relevant attached drawings.
Fig. 1 shows a kind of structural schematic diagram of gamma camera provided in an embodiment of the present invention.
Fig. 2 shows a kind of flow signals of the local spectral analysis method based on γ cameras provided in an embodiment of the present invention
Figure.
Fig. 3 shows a kind of schematic diagram of full visual field perspective view provided in an embodiment of the present invention.
Fig. 4 shows a kind of schematic diagram in gamma ray radiator region provided in an embodiment of the present invention.
Fig. 5 shows a kind of stream of the sub-step of local spectral analysis method based on γ cameras provided in an embodiment of the present invention
Journey schematic diagram.
Fig. 6 shows D provided in an embodiment of the present inventionARadiate the schematic diagram of the local spectrum of source region.
Fig. 7 shows D provided in an embodiment of the present inventionBRadiate the schematic diagram of the local spectrum of source region.
Fig. 8 shows the flow signal of another local spectral analysis method based on γ cameras provided in an embodiment of the present invention
Figure.
Fig. 9 shows a kind of schematic diagram of full visual field power spectrum provided in an embodiment of the present invention.
Figure 10 shows a kind of D provided in an embodiment of the present inventionAThe voxel for radiating source region divides schematic diagram.
Figure 11 shows a kind of function module of local spectrum analysis device based on γ cameras provided in an embodiment of the present invention
Schematic diagram.
Diagram：Local spectrum analysis devices of the 310- based on γ cameras；300- gamma cameras；320- memories；330- processing
Device；340- detectors；350- encoding board collimators；311- acquisition modules；312- establishes module；313- category determination modules；
314- area determination modules；The first computing modules of 315-；316- analog modules；The second computing modules of 317-；318- analysis modules；
319- comparison modules.
Specific implementation mode
Below in conjunction with attached drawing in the embodiment of the present invention, technical solution in the embodiment of the present invention carries out clear, complete
Ground describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Usually exist
The component of the embodiment of the present invention described and illustrated in attached drawing can be arranged and be designed with a variety of different configurations herein.Cause
This, the detailed description of the embodiment of the present invention to providing in the accompanying drawings is not intended to limit claimed invention below
Range, but it is merely representative of the selected embodiment of the present invention.Based on the embodiment of the present invention, those skilled in the art are not doing
The every other embodiment obtained under the premise of going out creative work, shall fall within the protection scope of the present invention.
It should be noted that：Similar label and letter indicate similar terms in following attached drawing, therefore, once a certain Xiang Yi
It is defined, then it further need not be defined and explained in subsequent attached drawing in a attached drawing.Meanwhile the present invention's
In description, term " first ", " second " etc. are only used for distinguishing description, are not understood to indicate or imply relative importance.
Local spectral analysis method provided in an embodiment of the present invention based on γ cameras is applied to gamma (γ) camera 300, asks
Referring to Fig.1, it is a kind of structural schematic diagram of gamma camera 300 provided in an embodiment of the present invention, which includes being based on
Local spectrum analysis device 310, memory 320, processor 330, detector 340 and the encoding board collimator 350 of γ cameras,
In,
The memory 320, processor 330, detector 340 and 350 each element of encoding board collimator are direct between each other
Or be electrically connected indirectly, to realize the transmission or interaction of data.For example, these elements can pass through one or more between each other
Communication bus or signal wire, which are realized, to be electrically connected.The local spectrum analysis device 310 based on γ cameras include it is at least one can
The operating system of gamma camera 300 is stored in memory 320 or is solidificated in the form of software or firmware (firmware)
Software function module in (operating system, OS).Processor 330 is used to execute holding of storing in memory 320
Row module, such as the software function module that includes of the local spectrum analysis device 310 based on γ cameras and computer program etc..
Wherein, memory 320 may be, but not limited to, random access memory (Random Access Memory,
RAM), read-only memory (Read Only Memory, ROM), programmable read only memory (Programmable Read-Only
Memory, PROM), erasable read-only memory (Erasable Programmable Read-Only Memory, EPROM),
Electricallyerasable ROM (EEROM) (Electric Erasable Programmable Read-Only Memory, EEPROM) etc..
Wherein, memory 320 is for storing program, and processor 330 executes described program after receiving and executing instruction.
The gamma camera 300 is for detecting gamma ray radiator, and gamma ray radiator correspondence sends out one or more photons, the photon
It is projeced on detector 340 by encoding board collimator 350 and forms full visual field perspective view, and then storage is executed by processor 330
Software function module in memory 320, such as the local spectrum analysis device 310 based on γ cameras are realized and delimit area to user
The nucleic classification of gamma ray radiator in domain is analyzed.
Fig. 2 is please referred to, is that a kind of flow of local spectral analysis method based on γ cameras provided in an embodiment of the present invention is shown
It is intended to, is somebody's turn to do the local spectral analysis method based on γ cameras and is applied to gamma camera 300, it should the local spectrum analysis side based on γ cameras
Method includes：
Step S110, obtains full visual field perspective view, and the full visual field perspective view is passed through by the photon that gamma ray radiator is sent out and compiled
Code plate collimator is projeced into the detector and is formed.
The gamma ray radiator sends out one or more photons, which is projeced by collimator 350 on detector 340, into
And full visual field perspective view as shown in Figure 3 is formed, which is a kind of showing for full visual field perspective view provided in an embodiment of the present invention
It is intended to.
Step S120 rebuilds full visual field source positions and intensity distribution according to projection reconstruction algorithm, and determination sends out photon
Multiple gamma ray radiator regions.
Since a kind of gamma ray radiator sends out one or more photons, a variety of photons collectively form full visual field perspective view, therefore,
The raw bits for the gamma ray radiator for being actually sent out photon are calculated by position of the photon in full visual field perspective view by algorithm for reconstructing
It sets, i.e. gamma ray radiator region, as shown in figure 4, being a kind of schematic diagram in gamma ray radiator region provided in an embodiment of the present invention, at this
In inventive embodiments, is rebuild by distribution of radiation it is found that there are two radioactive source in the full visual field, illustrated in the embodiment of the present invention
Two gamma ray radiator regions in the full visual field perspective view, are respectively defined as D hereinAAnd DB。
Step S130 calculates weight factor of each described gamma ray radiator region relative to the full visual field perspective view.
In embodiments of the present invention, it is to calculate DAAnd DBTwo gamma ray radiator regions relative to full visual field perspective view
Weight factor please refers to Fig. 5, the step of being a kind of local spectral analysis method based on γ cameras provided in an embodiment of the present invention
The flow diagram of the sub-step of S130, step S130 include：
Step S131 calculates all photon counting summations in the full visual field perspective view.
It is counted since the photon of the different-energy in the full visual field perspective view all has, calculate the full visual field projection
The counting summation of all actual measurement photons in figure.
Step S132 calculates photon counting summation in each described gamma ray radiator region；
The photon counting summation in each described gamma ray radiator region is calculated, is difference in embodiments of the present invention
Calculate DAAnd DBThe counting summation of photon in two gamma ray radiator regions.
Step S133, will be in the photon counting summation and the full visual field perspective view in each described gamma ray radiator region
The summations of all photon countings be divided by, obtain the weight factor.
The summation of the photon counting in each gamma ray radiator region is calculated in step S132, i.e., of the invention real
It applies in example, calculates separately out DAIn photon counting summation and DBIn photon counting summation, by DAIn photon
The summation of counting and the counting summation of all photons in full visual field perspective view are divided by, and D is obtainedARadiation source region relative to regarding entirely
The weight factor A of wild perspective view, by DBIn photon counting summation and the counting of all photons in full visual field perspective view it is total
Be divided by, obtain DBRadiate weight factor B of the source region relative to full visual field perspective view.
Step S140, the corresponding each energy of actual measurement of pixel that the pixel value in each gamma ray radiator region is not zero
The counting of amount section photon is multiplied with the weight factor respectively, obtains the energy point of the photon in each gamma ray radiator region
Cloth situation.
In embodiments of the present invention, respectively by DAAnd DBThe pixel value of two radiation source region separate analogue projections is not zero
The counting of each energy section photon of the corresponding actual measurement of pixel be multiplied with corresponding weight factor, in other words, by DARadioactive source
The counting for the corresponding each energy section photon of pixel that pixel value is not zero in region is multiplied with weight factor A, by DBRadioactive source
The counting for the corresponding each energy section photon of pixel that pixel value is not zero in region is multiplied with weight factor B, and then respectively
To the Energy distribution situation of the photon in each gamma ray radiator region.
Step S150 is obtained according to the Energy distribution situation of the photon in gamma ray radiator region described in each described in each
The local in gamma ray radiator region is composed.
In embodiments of the present invention, according to DAThe Energy distribution situation for radiating the photon of source region makes the radiation source region
Local spectrum, as shown in fig. 6, being D provided in an embodiment of the present inventionARadiate the local spectrum of source region.According to DBRadiate source region
The Energy distribution situation of photon makes the local spectrum of the radiation source region, as shown in fig. 7, being D provided in an embodiment of the present inventionBIt puts
Penetrate the local spectrum of source region.The DARadiate source region and DBThe local of radiation source region is composed is with different types of photon energy
Abscissa is counted as ordinate so that this kind of photon energy is corresponding.It is analyzed simultaneously by the local spectrum to gamma ray radiator region
It is compared with full visual field power spectrum, to understand the main nucleic classification in the radiation source region, is provided for nuclear industry associated production
Great convenience, as nucleic is migrated and early warning.
Fig. 8 is please referred to, is the flow of another local spectral analysis method based on γ cameras provided in an embodiment of the present invention
Schematic diagram is somebody's turn to do the local spectral analysis method based on γ cameras and is applied to gamma camera 300, should the local spectrum analysis based on γ cameras
Method includes：
Step S210 obtains full visual field perspective view, what the full visual field perspective view was sent out by all gamma ray radiators in the visual field
Photon is projeced into the detector across the encoding board collimator and is formed.
The gamma ray radiator sends out one or more photons, which is projeced by collimator 350 on detector 340, into
And full visual field perspective view as shown in Figure 3 is formed, which is a kind of showing for full visual field perspective view provided in an embodiment of the present invention
It is intended to
Step S220 is corresponded to according to the energy of the energy of all photons on the full visual field perspective view and each photon
Counting establish full visual field power spectrum.
It is established according to the corresponding counting of the energy of the energy of all photons on full visual field perspective view and each photon
Full visual field power spectrum as shown in Figure 9, as shown in figure 9, being a kind of schematic diagram of full visual field power spectrum provided in an embodiment of the present invention.It should
Full visual field power spectrum is counted as ordinate using the type of photon energy as abscissa with each photon energy is corresponding.
Step S230, according to the nucleic of gamma ray radiator described in the location determination at character pair energy peak in the full visual field power spectrum
Classification.
By the full visual field power spectrum it is found that there are two feature energy outstanding peaks altogether in the full visual field power spectrum, respectively
Corresponding to 0.05MeV and 0.15MeV photon energies, it follows that this is entirely within sweep of the eye, main nucleic classification is
The nucleic of 0.05MeV and 0.15MeV photon energies.
Step S240 rebuilds full visual field source positions and intensity distribution according to projection reconstruction algorithm, and determination sends out photon
Multiple gamma ray radiator regions.
Since different gamma ray radiators send out one or more photons, a variety of photons collectively form full visual field perspective view, therefore,
The raw bits for the gamma ray radiator for being actually sent out photon are calculated by position of the photon in full visual field perspective view by algorithm for reconstructing
It sets, i.e. gamma ray radiator region, as shown in figure 4, being a kind of schematic diagram in gamma ray radiator region provided in an embodiment of the present invention, at this
In inventive embodiments, two gamma ray radiator regions in the full visual field perspective view are illustrated, are respectively defined as D hereinAAnd DB。
Step S250 calculates weight factor of each described gamma ray radiator region relative to the full visual field perspective view.
In embodiments of the present invention, it is to calculate DAAnd DBTwo gamma ray radiator regions relative to full visual field perspective view
Weight factor.Since each photon in the full visual field perspective view has energy, it calculates in the full visual field perspective view
All photons counting summation, meanwhile, calculate separately DAAnd DBThe counting summation of photon in two gamma ray radiator regions.It will
DAIn photon energy summation and the energy summation of all photons in full visual field perspective view be divided by, obtain DARadiate source region
Relative to the weight factor A of full visual field perspective view, by DBIn photon counting summation and full visual field perspective view in all light
The energy summation of son is divided by, and D is obtainedBRadiate weight factor B of the source region relative to full visual field perspective view.
Step S260, using the different gamma ray radiator regions of Monte Carlo simulation, by the gamma ray radiator region point
At multiple regions, the independent projections in each gamma ray radiator region are simulated respectively.
As shown in Figure 10, it is a kind of D provided in an embodiment of the present inventionAThe region segmentation schematic diagram of source region is radiated, is used
Monte Carlo simulation DASource region is radiated, D is obtainedARadiation source region individually projects, it is readily appreciated that, using Monte Carlo simulation
DBSource region is radiated, D is obtainedBRadiation source region is divided into independent projection.
Step S270, the pixel that the pixel value of the independent projections in each gamma ray radiator region is not zero correspond to
Each energy section photon counting of actual measurement be multiplied with the weight factor, obtain the energy of the photon in each gamma ray radiator region
Measure distribution situation.
In embodiments of the present invention, respectively by DAAnd DBThe pixel pair that pixel value in two radiation source regions is not zero
The counting for each energy section photon answered is multiplied with corresponding weight factor, in other words, by DARadiate source region in pixel value not
The counting for the corresponding photon of pixel for being zero is multiplied with weight factor A, by DBThe picture that pixel value is not zero in radiation source region
The counting of the corresponding each energy section photon of vegetarian refreshments is multiplied with weight factor B, and then respectively obtains each gamma ray radiator region
The Energy distribution situation of photon.
Step S280 is obtained according to the Energy distribution situation of the photon in gamma ray radiator region described in each described in each
The local in gamma ray radiator region is composed.
In embodiments of the present invention, according to DAThe Energy distribution situation for radiating the photon of source region makes the radiation source region
Local spectrum, as shown in fig. 6, being D provided in an embodiment of the present inventionARadiate the local spectrum of source region.According to DBRadiate source region
The Energy distribution situation of photon makes the local spectrum of the radiation source region, as shown in fig. 7, being D provided in an embodiment of the present inventionBIt puts
Penetrate the local spectrum of source region.The DARadiate source region and DBThe local of radiation source region is composed is with different types of photon energy
Abscissa is counted as ordinate so that this kind of photon energy is corresponding.
Local spectrum is corresponded to the counting at feature energy peak and corresponding in the full visual field power spectrum by step S290
Counting at the feature energy peak of same position is compared；If the comparative counting at the local spectrum signature energy peak regards entirely more than described
The numerical value corresponded at same position feature energy peak in wild power spectrum, then the local is composed mainly contains in corresponding gamma ray radiator
Corresponding to the category nucleic of the energy of photon.
A kind of full visual field power spectrum provided in an embodiment of the present invention as shown in Figure 9, by the full visual field power spectrum it is found that at this
There are two feature energy outstanding peaks altogether in the power spectrum of the full visual field, respectively corresponding to 0.05MeV and 0.15MeV photon energies, thus
It is found that this is entirely within sweep of the eye, main nucleic classification is the corresponding nucleic of 0.05MeV and 0.15MeV photon energies.
As shown in fig. 6, being D provided in an embodiment of the present inventionARadiate the local spectrum of source region.By DARadiate the office of source region
It composes it is found that in the D in domainAIt is the corresponding energy peak of 0.15MeV photon energies to radiate wave crest more outstanding in the local spectrum of source region,
And DAThe corresponding energy peak of 0.15MeV photon energies is relative to 0.15MeV light in full visual field power spectrum in the local spectrum of radiation source region
The corresponding wave crest of sub- energy increases, therefore, in DARadiate the nucleic of the interior predominantly 0.15MeV photon energies of source region.
As shown in fig. 7, being D provided in an embodiment of the present inventionBRadiate the local spectrum of source region.By DBRadiate the office of source region
It composes it is found that in the D in domainBIt is the corresponding wave crest of 0.05MeV photon energies to radiate wave crest more outstanding in the local spectrum of source region,
And DBThe corresponding wave crest of 0.05MeV photon energies is relative to 0.05MeV light in full visual field power spectrum in the local spectrum of radiation source region
The corresponding wave crest of sub- energy increases, therefore, in DBRadiate the nucleic of the interior predominantly 0.05MeV photon energies of source region.
On the basis of full visual field power spectrum, comparative analysis difference radiates the nucleic classification of source region, to further realize root
The analysis of nucleic classification in the region delimited is needed according to user, while can be widely applied to different application scenarios, such as nucleic
It migrates or nucleic early warning etc., brings great convenience for the working practice of user.
It please refers to shown in Figure 11, is a kind of local spectrum analysis device 310 based on γ cameras provided in an embodiment of the present invention
High-level schematic functional block diagram, should local spectrum analysis device 310 based on γ cameras include acquisition module 311, establish module 312,
Category determination module 313, radioactive source determining module 314, the first computing module 315, analog module 316, the second computing module
317, analysis module 318 and comparison module 319.Wherein,
Acquisition module 311, for obtaining full visual field perspective view, the full visual field perspective view is by all gamma ray radiators in the visual field
The photon sent out is projeced into the detector across the encoding board collimator and is formed.
In embodiments of the present invention, step S110 or step S210 can be executed by acquisition module 311.
Module 312 is established, for according to the energy of all photons on the full visual field perspective view and each photon
Full visual field power spectrum is established in the corresponding counting of energy.
In embodiments of the present invention, step S220 can be executed by establishing module 312.
Category determination module 313, for according to γ described in the location determination at character pair energy peak in the full visual field power spectrum
The nucleic classification of radioactive source.
In embodiments of the present invention, step S230 can be executed by category determination module 313.
Area determination module 314, for rebuilding full visual field source positions according to projection reconstruction algorithm, determination sends out photon
Multiple gamma ray radiator regions.
In embodiments of the present invention, step S120 or step S240 can be executed by radioactive source determining module 314.
First computing module 315, for calculating each described gamma ray radiator region relative to the full visual field perspective view
Weight factor.
In embodiments of the present invention, step S130~step S133 or step S250 can be held by the first computing module 315
Row.
Analog module 316, for using the different gamma ray radiator regions of Monte Carlo simulation, by the gamma ray radiator
It is divided into multiple regions, simulates the independent projections in each gamma ray radiator region respectively.
In embodiments of the present invention, step S260 can be executed by analog module 316.
Second computing module 317, for the pixel value of the independent projections in each gamma ray radiator region to be not zero
Each energy section photon counting of the corresponding actual measurement of pixel be multiplied with the weight factor, obtain each described gamma ray radiator area
The Energy distribution situation of the photon in domain.
In embodiments of the present invention, step S140 or step S270 can be executed by the second computing module 317.
Analysis module 318, for being obtained often according to the Energy distribution situation of the photon in gamma ray radiator region described in each
The local spectrum in one gamma ray radiator region.
In embodiments of the present invention, step S150 or step S280 can be executed by analysis module 318.
Comparison module 319, for corresponding to local spectrum in the counting at feature energy peak and the full visual field power spectrum
The feature energy peak corresponding to same position at counting be compared；If the comparative counting at the local spectrum signature energy peak is more than
The numerical value corresponded at same position feature energy peak in the full visual field power spectrum, then in the corresponding gamma ray radiator of the local spectrum
Mainly contain the category nucleic of the energy corresponding to photon.
In embodiments of the present invention, step S290 can be executed by comparison module 319.
Due to being had been described in the local spectral analysis method part based on γ cameras, details are not described herein.
In conclusion a kind of local spectral analysis method and device based on γ cameras provided in an embodiment of the present invention, the party
Method includes obtaining full visual field perspective view, and the full visual field perspective view is thrown by the photon that gamma ray radiator is launched across the collimator
It penetrates and is formed in the detector.The multiple gamma ray radiators for sending out photon are determined in the full visual field perspective view according to algorithm for reconstructing
Region, and calculate weight factor of each described gamma ray radiator region relative to the full visual field perspective view.By each γ
The energy for the corresponding actual measurement photon of pixel that the pixel value of radiation source region independent projections is not zero and the weight factor phase
Multiply, the Energy distribution situation of the photon in each gamma ray radiator region is obtained, according to the photon in each gamma ray radiator region
Energy distribution situation obtains the local spectrum of each gamma ray radiator.By searching gamma ray radiator region in being imaged in the full visual field, with
It respectively to each gamma ray radiator regional analysis, and then realizes and analyzes nucleic classification in self defined area, realize nucleic early warning
Deng.
It should be noted that herein, relational terms such as first and second and the like are used merely to a reality
Body or operation are distinguished with another entity or operation, are deposited without necessarily requiring or implying between these entities or operation
In any actual relationship or order or sequence.Moreover, the terms "include", "comprise" or its any other variant are intended to
Non-exclusive inclusion, so that the process, method, article or equipment including a series of elements is not only wanted including those
Element, but also include other elements that are not explicitly listed, or further include for this process, method, article or equipment
Intrinsic element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that
There is also other identical elements in process, method, article or equipment including the element.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the invention may be variously modified and varied.All within the spirits and principles of the present invention, any made by repair
Change, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.It should be noted that：Similar label and letter exist
Similar terms are indicated in following attached drawing, therefore, once being defined in a certain Xiang Yi attached drawing, are then not required in subsequent attached drawing
It is further defined and is explained.
Claims (10)
1. a kind of local spectral analysis method based on γ cameras, which is characterized in that the γ cameras include encoding board collimator and
Detector, the method includes：
Full visual field perspective view is obtained, the full visual field perspective view is passed through described by the photon that all gamma ray radiators are launched in the visual field
Encoding board collimator is projeced into the detector and is formed；
Full visual field source positions and intensity distribution are rebuild according to projection reconstruction algorithm, determine the multiple gamma ray radiators for sending out photon
Region；
Calculate weight factor of each described gamma ray radiator region relative to the full visual field perspective view；
The corresponding each energy of actual measurement of pixel that the pixel value of the independent projections in each gamma ray radiator region is not zero
Section photon counting is multiplied with the weight factor, obtains the Energy distribution situation of the photon in each gamma ray radiator region；
Each described gamma ray radiator region is obtained according to the Energy distribution situation of the photon in gamma ray radiator region described in each
Local spectrum.
2. the local spectral analysis method based on γ cameras as described in claim 1, which is characterized in that described each institute of calculating
The step of stating weight factor of the gamma ray radiator region relative to the full visual field perspective view include：
Calculate all photon counting summations in the full visual field perspective view；
Calculate photon counting summation in each described gamma ray radiator region；
By the photon counting summation in each described gamma ray radiator region and all photon meters in the full visual field perspective view
Several summations are divided by, and the weight factor is obtained.
3. the local spectral analysis method based on γ cameras as described in claim 1, which is characterized in that described to obtain full visual field throwing
Further include step after shadow figure：
It is established according to the corresponding counting of the energy of the energy of all photons on the full visual field perspective view and each photon complete
Visual field power spectrum；
According to the nucleic classification of gamma ray radiator described in the location determination at character pair energy peak in the full visual field power spectrum.
4. the local spectral analysis method based on γ cameras as claimed in claim 3, which is characterized in that the method further includes：
Local spectrum is corresponded into the counting at feature energy peak and the spy corresponding to same position in the full visual field power spectrum
Counting at sign energy peak is compared；If the comparative counting at the local spectrum signature energy peak is more than pair in the full visual field power spectrum
Should be in the numerical value at same position feature energy peak, then the local, which is composed, mainly contains in corresponding gamma ray radiator corresponding to photon
The category nucleic of energy.
5. the local spectral analysis method based on γ cameras as described in claim 1, which is characterized in that described each institute of calculating
The weight factor that gamma ray radiator region is stated relative to the full visual field perspective view further includes later step：
Using the different gamma ray radiator regions of Monte Carlo simulation, the gamma ray radiator region is divided into multiple regions, point
The independent projections in each gamma ray radiator region are not simulated.
6. a kind of local spectrum analysis device based on γ cameras, which is characterized in that the γ cameras include encoding board collimator and
Detector, described device include：
Acquisition module, for obtaining full visual field perspective view, what the full visual field perspective view was sent out by all gamma ray radiators in the visual field
Photon is projeced into the detector across the encoding board collimator and is formed；
Area determination module, for rebuilding full visual field source positions according to projection reconstruction algorithm, determination sends out the multiple of photon
Gamma ray radiator region；
First computing module, for calculate each described gamma ray radiator region relative to the full visual field perspective view weight because
Son；
Second computing module, the pixel for the pixel value of the independent projections in each gamma ray radiator region to be not zero
Corresponding each energy section photon counting of actual measurement is multiplied with the weight factor, obtains the photon in each gamma ray radiator region
Energy distribution situation；
Analysis module, for being obtained described in each according to the Energy distribution situation of the photon in gamma ray radiator region described in each
The local in gamma ray radiator region is composed.
7. the local spectrum analysis device based on γ cameras as claimed in claim 6, which is characterized in that first computing module
It is additionally operable to：
Calculate all photon counting summations in the full visual field perspective view；
Calculate the photon counting summation in each described radioactive area γ；
By the photon counting summation in each described gamma ray radiator region and all photon meters in the full visual field perspective view
Number summation is divided by, and the weight factor is obtained.
8. the local spectrum analysis device based on γ cameras as claimed in claim 6, which is characterized in that described device further includes：
Module is established, for being corresponded to according to the energy of all photons on the full visual field perspective view and the energy of each photon
Counting establish full visual field power spectrum；
Category determination module, for according to gamma ray radiator described in the location determination at character pair energy peak in the full visual field power spectrum
Nucleic classification.
9. the local spectrum analysis device based on γ cameras as claimed in claim 8, which is characterized in that described device further includes：
Comparison module, for local spectrum to be corresponded to the counting at feature energy peak and corresponding in the full visual field power spectrum
Counting at the feature energy peak of same position is compared；If the comparative counting at the local spectrum signature energy peak regards entirely more than described
The numerical value corresponded at same position feature energy peak in wild power spectrum, then the local is composed mainly contains in corresponding gamma ray radiator
Corresponding to the category nucleic of the energy of photon.
10. the local spectrum analysis device based on γ cameras as claimed in claim 6, which is characterized in that described device further includes：
Analog module, for using the different gamma ray radiator regions of Monte Carlo simulation, the gamma ray radiator being divided into more
The independent projections in each gamma ray radiator region are simulated in a region respectively.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810136044.5A CN108459344B (en) | 2018-02-09 | 2018-02-09 | Local spectral analysis method and device based on γ camera |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810136044.5A CN108459344B (en) | 2018-02-09 | 2018-02-09 | Local spectral analysis method and device based on γ camera |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108459344A true CN108459344A (en) | 2018-08-28 |
| CN108459344B CN108459344B (en) | 2019-09-06 |
Family
ID=63239973
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810136044.5A Active CN108459344B (en) | 2018-02-09 | 2018-02-09 | Local spectral analysis method and device based on γ camera |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108459344B (en) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3978337A (en) * | 1975-01-29 | 1976-08-31 | Wisconsin Alumni Research Foundation | Three-dimensional time-of-flight gamma camera system |
| JPH04168396A (en) * | 1990-10-31 | 1992-06-16 | Shimadzu Corp | Radiation measurement circuit |
| US20090122958A1 (en) * | 2007-07-09 | 2009-05-14 | Lucian Mihailescu | Hybrid compton camera/coded aperture imaging system |
| CN102540238A (en) * | 2010-12-31 | 2012-07-04 | 同方威视技术股份有限公司 | Gamma camera and method for detecting radiation ray by utilizing same |
| CN202886614U (en) * | 2012-09-07 | 2013-04-17 | 北京中科核安科技有限公司 | Radioactive source position indicator |
| CN103163548A (en) * | 2013-03-07 | 2013-06-19 | 北京辛耕普华医疗科技有限公司 | Radioactive substance detection method, device and system based on gamma camera |
| CN103954987A (en) * | 2014-05-15 | 2014-07-30 | 四川材料与工艺研究所 | Method for measuring mass of radioactive substance through gamma camera and gamma spectrometer |
| US20140301535A1 (en) * | 2013-04-05 | 2014-10-09 | The Arizona Board Of Regents On Behalf Of The University Of Arizona | Gamma Ray Imaging Systems And Methods |
| CN104166153A (en) * | 2014-07-25 | 2014-11-26 | 北京辛耕普华医疗科技有限公司 | Method and device for measuring two-dimensional angle distribution of radiation dose rate of radioactive substance |
| CN104237959A (en) * | 2014-08-21 | 2014-12-24 | 北京辛耕普华医疗科技有限公司 | Method and device for conducting real-time dynamic tracking positioning on radioactive substances |
| CN105974461A (en) * | 2016-05-04 | 2016-09-28 | 中国原子能科学研究院 | Gamma ray scanning imaging and radionuclide identification system and method thereof |
-
2018
- 2018-02-09 CN CN201810136044.5A patent/CN108459344B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3978337A (en) * | 1975-01-29 | 1976-08-31 | Wisconsin Alumni Research Foundation | Three-dimensional time-of-flight gamma camera system |
| JPH04168396A (en) * | 1990-10-31 | 1992-06-16 | Shimadzu Corp | Radiation measurement circuit |
| US20090122958A1 (en) * | 2007-07-09 | 2009-05-14 | Lucian Mihailescu | Hybrid compton camera/coded aperture imaging system |
| CN102540238A (en) * | 2010-12-31 | 2012-07-04 | 同方威视技术股份有限公司 | Gamma camera and method for detecting radiation ray by utilizing same |
| CN202886614U (en) * | 2012-09-07 | 2013-04-17 | 北京中科核安科技有限公司 | Radioactive source position indicator |
| CN103163548A (en) * | 2013-03-07 | 2013-06-19 | 北京辛耕普华医疗科技有限公司 | Radioactive substance detection method, device and system based on gamma camera |
| US20140301535A1 (en) * | 2013-04-05 | 2014-10-09 | The Arizona Board Of Regents On Behalf Of The University Of Arizona | Gamma Ray Imaging Systems And Methods |
| CN103954987A (en) * | 2014-05-15 | 2014-07-30 | 四川材料与工艺研究所 | Method for measuring mass of radioactive substance through gamma camera and gamma spectrometer |
| CN104166153A (en) * | 2014-07-25 | 2014-11-26 | 北京辛耕普华医疗科技有限公司 | Method and device for measuring two-dimensional angle distribution of radiation dose rate of radioactive substance |
| CN104237959A (en) * | 2014-08-21 | 2014-12-24 | 北京辛耕普华医疗科技有限公司 | Method and device for conducting real-time dynamic tracking positioning on radioactive substances |
| CN105974461A (en) * | 2016-05-04 | 2016-09-28 | 中国原子能科学研究院 | Gamma ray scanning imaging and radionuclide identification system and method thereof |
Non-Patent Citations (2)
| Title |
|---|
| 丁长骥等: "康普顿相机的GEANT4模拟与反投影图像重建", 《核技术 》 * |
| 陈立宏等: "基于MURA编码孔径准直器核辐射成像系统设计", 《核技术》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108459344B (en) | 2019-09-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103163548B (en) | Based on the radioactive substance detection of gamma camera and device thereof and system | |
| Hamilton et al. | The origin of the diffuse X-ray background | |
| Marshall et al. | Modelling the Galactic interstellar extinction distribution in three dimensions | |
| Uehara et al. | Measurement of γ γ*→ π 0 transition form factor at Belle | |
| Mohrmann et al. | Validation of open-source science tools and background model construction in γ-ray astronomy | |
| US20160104290A1 (en) | Image based object locator | |
| US8787620B2 (en) | Automated crystal identification achieved via watershed segmentation | |
| CN103745185B (en) | Method and apparatus to recognize the position of detector crystal unit | |
| US20120112087A1 (en) | Apparatus and method for detecting gamma-ray direction | |
| Prieto et al. | Twelve automated thresholding methods for segmentation of PET images: a phantom study | |
| Pierbattista et al. | Young and middle age pulsar light-curve morphology: Comparison of Fermi observations with γ-ray and radio emission geometries | |
| Schulz et al. | A redshift-dependent IRX–β dust attenuation relation for TNG50 galaxies | |
| KR101766294B1 (en) | Device and method for processing 3-dimensional visualization of radiation source distribution using coded-apertures | |
| Dong et al. | Measuring the integrated Sachs–Wolfe effect from the low-density regions of the universe | |
| Riggi et al. | Caesar source finder: recent developments and testing | |
| Jones et al. | Disentangling overlapping astronomical sources using spatial and spectral information | |
| Hewett et al. | Quasar surveys | |
| Wahl et al. | Gamma-ray point-source detection in unknown background using 3D-position-sensitive semiconductor detectors | |
| EP3796049A1 (en) | Substance identification device and method for extracting statistical feature based on cluster analysis | |
| CN108459344B (en) | Local spectral analysis method and device based on γ camera | |
| CN106022217A (en) | Civil airport runway area detection method free from supervision multistage classification | |
| KR102249120B1 (en) | Analysis program and analysis system for radioactive contamination in and outside the body of the subject | |
| CN112269205B (en) | Method for determining parameters of a radiation detector | |
| Kowash | A rotating modulation imager for the orphan source search problem | |
| Capalbo et al. | Morphological analysis of SZ and X-ray maps of galaxy clusters with Zernike polynomials |
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 |