CN103645491B - Radioactive source localization method, Apparatus and system - Google Patents

Radioactive source localization method, Apparatus and system Download PDF

Info

Publication number
CN103645491B
CN103645491B CN201310602906.6A CN201310602906A CN103645491B CN 103645491 B CN103645491 B CN 103645491B CN 201310602906 A CN201310602906 A CN 201310602906A CN 103645491 B CN103645491 B CN 103645491B
Authority
CN
China
Prior art keywords
crystal
radioactive source
count
counting
angle
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201310602906.6A
Other languages
Chinese (zh)
Other versions
CN103645491A (en
Inventor
魏龙
张译文
李道武
章志明
王宝义
帅磊
王英杰
朱美玲
孙世峰
张玉包
唐浩辉
黄先超
柴培
李婷
庄凯
姜小盼
刘彦韬
周魏
曾凡剑
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of High Energy Physics of CAS
Original Assignee
Institute of High Energy Physics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of High Energy Physics of CAS filed Critical Institute of High Energy Physics of CAS
Priority to CN201310602906.6A priority Critical patent/CN103645491B/en
Publication of CN103645491A publication Critical patent/CN103645491A/en
Application granted granted Critical
Publication of CN103645491B publication Critical patent/CN103645491B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Radiation (AREA)

Abstract

The invention discloses a kind of radioactive source localization method, Apparatus and system.The method includes: according to crystal count in the crystal array of analog data acquisition detector than the relation with radioactive source angle;Measure the actual measurement crystal count of described crystal array during radioactive source emitting radiation to be measured;In the excursion of counting noise, multiple denoising crystal count ratios are obtained according to described actual measurement crystal count and each described counting noise, radioactive source angle, and radioactive source angle that actual measurement radioactive source angle be described to be measured radioactive source when choosing counting noise optimum is respectively surveyed with each described denoising crystal count than corresponding than with the Relation acquisition of radioactive source angle according to described crystal count.This device includes: processing module, measurement module and computing module.This system includes: radioactive source positioner, radioactive source to be measured and detector.The present invention uses crystal count to be used for calculating parameter, uses and effectively deducts noise technique, improves sensitivity and the accuracy of radioactive source location.

Description

Radioactive source localization method, Apparatus and system
Technical field
The present invention relates to radioactive source Detection Techniques, particularly to a kind of radioactive source localization method, Apparatus and system.
Background technology
Radioactive source detection can be widely applied to various radioactive environment, such as laboratory, nuclear power plant, nuclear fuel factory, large-scale heavy Ion accelerator, radioactive source storehouses etc..Varied for detecting the detector of gamma ray, the most more advanced technology is to utilize Particle in gamma ray can be used four pieces of crystal groups by the principle of the absorption of crystals such as cesium iodide (CsI) or bismuth germanium oxide (BGO) Become 2 × 2 crystal arrays, same radioactive source is detected and counts, releases radioactive source according to the count value of each piece of crystal is counter Angle, to determine the orientation of radioactive source.
But, when radioactive source position and crystal array the most at grade time, the radioactive source angle calculated is inadequate Accurately so that positional accuracy is low;Simultaneously as crystal count value is inconspicuous with the change of angle so that existing radioactive source Localization method sensitivity is the highest, and resolution capability is on the low side.
Summary of the invention
Brief overview about the present invention given below, in order to the basic reason about certain aspects of the invention is provided Solve.Should be appreciated that this general introduction is not that the exhaustive about the present invention is summarized.It is not intended to determine the key of the present invention Or pith, nor is it intended to limit the scope of the present invention.Its purpose is only to provide some concept in simplified form, with This is as the preamble in greater detail discussed after a while.
The present invention provides a kind of radioactive source localization method, Apparatus and system, in order to solve radioactive source location in prior art The defect that accuracy is low and sensitivity is low, it is achieved high sensitivity, the radioactive source location of high accuracy.
The invention provides a kind of radioactive source localization method, including:
Step S1, according to reason intend data acquisition detector crystal array in crystal count than the pass with radioactive source angle System;
Step S2, measure the actual measurement crystal count of described crystal array during radioactive source emitting radiation to be measured;
Step S3, counting noise excursion in, according to described actual measurement crystal count and each described counting noise Obtain multiple denoising crystal count ratio, brilliant with each described denoising with the Relation acquisition of radioactive source angle according to described crystal count ratio Body counting is than corresponding each actual measurement radioactive source angle, and actual measurement radioactive source angle when choosing described counting noise optimum is institute State the radioactive source angle of radioactive source to be measured.
Present invention also offers a kind of radioactive source positioner, including:
Processing module, for according to analog data obtain detector crystal array in crystal count than with radioactive source angle Relation;
Measurement module, for measuring the actual measurement crystal count of described crystal array during radioactive source emitting radiation to be measured;
Computing module, in the excursion at counting noise, according to each described counting noise and described measurement mould Described actual measurement crystal count that block measurement obtains obtains multiple denoising crystal count ratios, obtain according to described processing module described in Crystal count respectively surveys radioactive source angle with each described denoising crystal count than corresponding than with the Relation acquisition of radioactive source angle Degree, and the radioactive source angle that actual measurement radioactive source angle is described radioactive source to be measured when choosing described counting noise optimum.
Present invention also offers a kind of radioactive source alignment system, including: above-mentioned radioactive source positioner, radioactive source to be measured and Detector.
The radioactive source localization method of present invention offer, Apparatus and system, use crystal count to be used for calculating parameter, make ginseng Number is obvious with the change of radioactive source angle, uses and effectively deducts noise technique, for radioactive source not in crystal array plane Situation, it is possible to be effectively improved in parameter the changed factor relevant with radioactive source angle, it is to avoid the defect of prior art, improve The sensitivity of radioactive source location and accuracy.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing In having technology to describe, the required accompanying drawing used is briefly described, it should be apparent that, the accompanying drawing in describing below is only this Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, it is also possible to Other accompanying drawing is obtained according to these accompanying drawings.
Fig. 1 is the flow chart of radioactive source localization method first embodiment of the present invention;
Fig. 2 is the flow chart of radioactive source localization method the second embodiment of the present invention;
Fig. 3 is the schematic diagram of 2 × 2 crystal arrays that the embodiment of the present invention uses;
Fig. 4 is the graph of a relation of embodiment of the present invention radioactive source angle and crystal count ratio;
Fig. 5 is the structural representation of radioactive source positioner first embodiment of the present invention;
Fig. 6 is the structural representation of radioactive source positioner the second embodiment of the present invention;
Fig. 7 is the system block diagram of radioactive source alignment system embodiment of the present invention.
Detailed description of the invention
For making the purpose of the embodiment of the present invention, technical scheme and advantage clearer, below in conjunction with the embodiment of the present invention In accompanying drawing, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is The a part of embodiment of the present invention rather than whole embodiments.Described in the accompanying drawing of the present invention or a kind of embodiment Element and feature can combine with the element shown in one or more other accompanying drawing or embodiment and feature.Should Note, for purposes of clarity, accompanying drawing and explanation eliminate unrelated to the invention, known to persons of ordinary skill in the art Parts and the expression of process and description.Based on the embodiment in the present invention, those of ordinary skill in the art are not paying creation Property work on the premise of the every other embodiment that obtained, broadly fall into the scope of protection of the invention.
Fig. 1 is the flow chart of radioactive source localization method first embodiment of the present invention.As it is shown in figure 1, the embodiment of the present invention carries Supply a kind of radioactive source localization method, including:
Step S10, according to analog data obtain detector crystal array in crystal count than the pass with radioactive source angle System, wherein, crystal count than for two pieces of crystal adjacent in crystal array counting and with the counting of crystal array and ratio;
Step S20, the actual measurement crystal count of crystal array when measuring radioactive source emitting radiation to be measured, wherein, survey crystal Counting is the crystal count that the actual measurement of crystal array obtains;
Step S30, in the excursion of counting noise, obtain multiple according to actual measurement crystal count and each counting noise Denoising crystal count ratio, corresponding with each denoising crystal count ratio with the Relation acquisition of radioactive source angle according to crystal count ratio Each actual measurement radioactive source angle, and the radioactive source angle that actual measurement radioactive source angle is radioactive source to be measured when choosing counting noise optimum Degree, wherein, denoising crystal count is than the crystal count ratio calculated for the crystal count after utilizing denoising.
In embodiments of the present invention, first according to analog data acquisition crystal count than the relation with radioactive source angle, should Analog data is situation about being placed on same plane with crystal array by simulation softward simulated emission source and obtains, and is derived from Relation in count noise can be considered zero.Then, for certain environment to be measured, utilize the crystal array of above-mentioned detector to measure actual measurement Crystal count.When radioactive source to be measured and crystal array the most at grade time, the counting on all crystal all will be above with Counting in one plane, the increase (being, count noise) of this counting makes crystal count ratio deviate above-mentioned crystal count ratio Relation with radioactive source angle.Therefore, it can determine the excursion of counting noise according to concrete measured data.To actual measurement crystal Counting carries out denoising, obtains denoising crystal count ratio, and obtains corresponding actual measurement radioactive source angle.Finally, at meter In the excursion of number noise, choose actual measurement radioactive source angle during counting noise optimum, as the radioactive source of radioactive source to be measured Angle.
The embodiment of the present invention provide radioactive source localization method, use crystal count be used for calculate parameter, make parameter with The change of radioactive source angle is obvious, uses and effectively deducts noise technique, for the radioactive source not feelings in crystal array plane Condition, it is possible to be effectively improved in parameter the changed factor relevant with radioactive source angle, it is to avoid the defect of prior art, improves and puts Penetrate sensitivity and the accuracy of source location.
Fig. 2 is the flow chart of radioactive source localization method the second embodiment of the present invention.As in figure 2 it is shown, at said method first On the basis of embodiment, above-mentioned crystal array can be 2 × 2 crystal arrays, and above-mentioned steps S10 can specifically include further:
Step S101, in 2 × 2 crystal arrays optional one piece of crystal as reference crystal;
Step S102, be first crystal group by reference crystal and one piece of crystal combinations being adjacent, by reference crystal and Another block crystal combinations being adjacent is the second crystal group;
Step S103, simulation softward is utilized to measure the first crystal counting of first crystal group than the of AB and the second crystal group Two crystal count are than BC, wherein, AB=(A+B)/(A+B+C+D), BC=(B+C)/(A+B+C+D), A, B, C and D be followed successively by this 2 × The crystal count of four pieces of crystal in 2 crystal arrays;
Step S104, obtain first crystal counting than the discrete relationship of AB and radioactive source angle, θ and the second crystal respectively The counting discrete relationship than BC Yu radioactive source angle, θ;
Step S105, above-mentioned discrete relationship is fitted to linear relationship.
Fig. 3 is the schematic diagram of 2 × 2 crystal arrays that the embodiment of the present invention uses.As it is shown on figure 3, a, b, c and d are 2 × 2 Four pieces of crystal in crystal array.If, selecting crystal b is reference crystal, and crystal a and crystal b is combined as first crystal combination, Crystal b and crystal c are combined as the second crystal combinations.A, B, C and D are followed successively by the counting of crystal a, b, c and d, utilize simulation soft Part is measured first crystal counting and is compared BC than AB and the second crystal count.Obtain respectively the discrete relationship of AB Yu θ and BC Yu θ from The relation of dissipating.Fig. 4 is the graph of a relation of embodiment of the present invention radioactive source angle and crystal count ratio.As shown in Figure 4, the relation of AB Yu θ Can be divided into two parts: θ ∈ (-90 °, 90 °) and θ ∈ (90 °, 270 °), in two intervals, AB Yu θ relation is all close to line Sexual relationship.Equally, BC is also close to linear relationship in θ ∈ (0 °, 180 °) with two intervals of θ ∈ (180 °, 360 °).Therefore, Being fitted it, obtain linear relationship θ=(AB-y1)/x1, θ=(BC-y2)/x2, wherein, x1, x2, y1 and y2 are that straight line is intended Close parameter.
The embodiment of the present invention, when detector works, can directly calculate AB and BC by measuring, therefore can pass through Above-mentioned formula calculates θ, calculates extremely simple, simultaneously figure 4, it is seen that crystal count is than the change with radioactive source angle Substantially, therefore location sensitivity is the highest.
Further, during radioactive source emitting radiation to be measured, the actual measurement crystal count of crystal a, b, c and d can be counted successively For A ', B ', C ' and D ', the excursion of counting noise N can be set to Nmin to Nmax.General N min can be set to 0, and Nmax can To be set in A ', B ', C ' and D ' the 60% of minima.Above-mentioned steps S30 can specifically include:
Step S301, counting noise N is set to Nmin;
Step S302, obtain the first denoising crystal count according to AB '=(A '+B '-2 × N)/(A '+B '+C '+D '-4 × N) Ratio, obtains the second denoising crystal count ratio according to BC '=(B '+C '-2 × N)/(A '+B '+C '+D '-4 × N);
Step S303, obtain, according to linear relationship, the radioactive source angle, θ that radioactive source angle, θ 1 and BC ' corresponding for AB ' is corresponding 2;
Step S304, counting noise N are from adding 1;
Step S305, compare | θ 1-θ 2 | andSize,WithFor parameter preset, if Parameter is then setOtherwise, parameterWithConstant;
Step S306, judge to count noise N whether more than Nmax, the most then perform step S307, otherwise, perform step S302;
Step S307, basisObtain the radioactive source angle of radioactive source to be measured.
When radioactive source to be measured and crystal array the most at grade time, the counting on all crystal all will be above same Counting in plane, the increase (being, count noise) of this counting makes crystal count than shown in AB and BC slip chart 4 Linearity curve.In the embodiment of the present invention, it is assumed that the counting noise on crystal a, b, c and d is N, then obtain the crystal after denoising Counting ratio:
AB '=(A '-N+B '-N)/(A '-N+B '-N+C '-N+D '-N)
=(A '+B '-2 × N)/(A '+B '+C '+D '-4 × N),
BC '=(B '-N+C '-N)/(A '-N+B '-N+C '-N+D '-N)
=(B '+C '-2 × N)/(A '+B '+C '+D '-4 × N).
According to practical situation due to corresponding same radioactive source, angle should be identical.Therefore, the determination side of noise is counted Method is, changes noise size the most continuously so that closest with the angle that curve BC calculates by curve AB, this Time counting noise, be closest to real optimum counting noise.The angle now calculated is closest to putting really Penetrate source actual angle.
Radioactive source to be measured in the embodiment of the present invention can be gamma-ray emission source.
The embodiment of the present invention provide radioactive source localization method, use crystal count be used for calculate parameter, make parameter with The change of radioactive source angle is obvious, uses and effectively deducts noise technique, for the radioactive source not feelings in crystal array plane Condition, removes unhelpful counting noise, it is possible to be effectively improved in parameter the changed factor relevant with radioactive source angle effectively, it is to avoid The defect of prior art, improves sensitivity and the accuracy of radioactive source location.
Fig. 5 is the structural representation of radioactive source positioner first embodiment of the present invention.As it is shown in figure 5, the present invention implements Example provides a kind of radioactive source positioner, including: processing module 51, measurement module 52 and computing module 53.Wherein, mould is processed Block 51 in the crystal array obtaining detector according to analog data crystal count than the relation with radioactive source angle;Measure mould Block 52 is the actual measurement crystal count of crystal array in time measuring radioactive source emitting radiation to be measured;Computing module 53 is for making an uproar at counting In the excursion of sound, measure, according to each counting noise and measurement module 52, the actual measurement crystal count obtained and obtain multiple denoisings Crystal count ratio, the crystal count obtained according to processing module 51 than with the Relation acquisition of radioactive source angle and each denoising crystal meter Number is than corresponding each actual measurement radioactive source angle, and actual measurement radioactive source angle when choosing counting noise optimum is radioactive source to be measured Radioactive source angle.
In embodiments of the present invention, first processing module 51 according to analog data obtain crystal count than with radioactive source angle Relation, wherein, crystal count than for two pieces of crystal adjacent in crystal array counting and with the counting of crystal array and ratio. This analog data is situation about being placed on same plane with crystal array by simulation softward simulated emission source and obtains, and thus obtains The crystal count obtained can be considered zero than with counting noise in the relation of radioactive source angle.Then, for certain environment to be measured, measure mould Block 52 utilizes the crystal array of above-mentioned detector to measure actual measurement crystal count.When radioactive source to be measured and crystal array be not same flat Time on face, the counting on all crystal all will be above counting at grade, and (being, counting is made an uproar in the increase of this counting Sound) make crystal count than deviateing above-mentioned linear relationship.Therefore, computing module 53 can determine counting according to concrete measured data The excursion of noise.Actual measurement crystal count is carried out denoising, obtains denoising crystal count ratio, and obtain corresponding Actual measurement radioactive source angle.Finally, in the excursion of counting noise, choose actual measurement radioactive source angle during counting noise optimum Degree, as the radioactive source angle of radioactive source to be measured.
The embodiment of the present invention provide radioactive source positioner, use crystal count be used for calculate parameter, make parameter with The change of radioactive source angle is obvious, uses and effectively deducts noise technique, for the radioactive source not feelings in crystal array plane Condition, it is possible to be effectively improved in parameter the changed factor relevant with radioactive source angle, it is to avoid the defect of prior art, improves and puts Penetrate sensitivity and the accuracy of source location.
Fig. 6 is the structural representation of radioactive source positioner the second embodiment of the present invention.As shown in Figure 6, at said apparatus On the basis of first embodiment, above-mentioned crystal array can be 2 × 2 crystal arrays, and above-mentioned processing module 51 can be wrapped further Include: select submodule 61, measure submodule the 62, first process submodule 63 and matching submodule 64.Wherein, submodule 61 is selected For one piece of crystal optional in 2 × 2 crystal arrays as reference crystal, by reference crystal and one piece of crystal group being adjacent It is combined into first crystal group, is the second crystal group by reference crystal and another block crystal combinations of being adjacent;Measure submodule 62 First crystal for utilizing simulation softward to measure first crystal group counts than AB and the second crystal count ratio of the second crystal group BC, wherein, AB=(A+B)/(A+B+C+D), BC=(B+C)/(A+B+C+D), A, B, C and D are followed successively by this 2 × 2 crystal array The crystal count of four pieces of crystal;First processes submodule 63 for obtaining first crystal counting respectively than AB and radioactive source angle, θ Discrete relationship and second crystal count discrete relationship than BC Yu radioactive source angle, θ;Matching submodule 64 is for by above-mentioned The discrete relationship that first process submodule 63 gets fits to linear relationship.
In figure 3, four pieces of crystal during a, b, c and d are 2 × 2 crystal arrays.If, select submodule 61 to select the crystal b to be Reference crystal, is combined as crystal a and crystal b first crystal combination, crystal b and crystal c is combined as the second crystal combinations.A、 B, C and D are followed successively by the counting of crystal a, b, c and d, measure submodule 62 utilize simulation softward measure first crystal counting than AB and Second crystal count compares BC.First processes submodule 63 obtains discrete relationship and the discrete relationship of BC Yu θ of AB Yu θ respectively. As shown in Figure 4, the relation of AB Yu θ can be divided into two parts: θ ∈ (-90 °, 90 °) and θ ∈ (90 °, 270 °), in two intervals In, AB Yu θ relation is all close to linear relationship.Equally, BC in θ ∈ (0 °, 180 °) and two intervals of θ ∈ (180 °, 360 °) also It is close to linear relationship.Therefore, it is fitted by matching submodule 64, obtains linear relationship θ=(AB-y1)/x1, θ=(BC- Y2)/x2, wherein, x1, x2, y1 and y2 are Straight Line Fitting Parameters.
The embodiment of the present invention, when detector works, can directly calculate AB and BC by measuring, therefore can pass through Above-mentioned formula calculates θ, calculates extremely simple, simultaneously figure 4, it is seen that crystal count is than the change with radioactive source angle Substantially, therefore location sensitivity is the highest.
Further, during radioactive source emitting radiation to be measured, the actual measurement crystal count of crystal a, b, c and d can be counted successively For A ', B ', C ' and D ', the excursion of counting noise N can be set to Nmin to Nmax.General N min can be set to 0, and Nmax can To be set in A ', B ', C ' and D ' the 60% of minima.
Above-mentioned computing module 53 can specifically include: arranges submodule the 65, first calculating sub module 66, second and processes submodule Block 67, certainly add submodule 68, comparison sub-module 69, judge submodule 71 and the second calculating sub module 72.Wherein, submodule is set 65 for being set to Nmin by counting noise N;First calculating sub module 66 for according to AB '=(A '+B '-2 × N)/(A '+B '+ C '+D '-4 × N) obtain the first denoising crystal count ratio, obtain according to BC '=(B '+C '-2 × N)/(A '+B '+C '+D '-4 × N) Second denoising crystal count ratio;Second processes submodule 67 obtains AB ' for the linear relationship according to matching submodule 64 matching The radioactive source angle, θ 2 that corresponding radioactive source angle, θ 1 and BC ' is corresponding;From adding submodule 68 for counting noise N is added 1 certainly;Ratio Relatively submodule 69 is used for comparing | θ 1-θ 2 | andSize,WithFor parameter preset, if Parameter is then setOtherwise, parameterWithConstant;Judge that submodule 71 is used for judging from adding Whether the counting noise N that submodule 68 draws is more than Nmax;Second calculating sub module 72 is for when counting noise N is more than Nmax Time, according toObtain the radioactive source angle of radioactive source to be measured.
When radioactive source to be measured and crystal array the most at grade time, the counting on all crystal will be above same Counting in plane, the increase (being, count noise) of this counting makes crystal count than shown in AB and BC slip chart 4 Linearity curve.In the embodiment of the present invention, it is assumed that the counting noise on crystal a, b, c and d is N, then obtain the crystal after denoising Counting ratio:
AB '=(A '-N+B '-N)/(A '-N+B '-N+C '-N+D '-N)
=(A '+B '-2 × N)/(A '+B '+C '+D '-4 × N),
BC '=(B '-N+C '-N)/(A '-N+B '-N+C '-N+D '-N)
=(B '+C '-2 × N)/(A '+B '+C '+D '-4 × N).
According to practical situation due to corresponding same radioactive source, angle should be identical.Therefore, the most continuously Change noise size so that closest with the angle that curve BC calculates by curve AB, counting noise now, it is and connects most Nearly real optimum counting noise.The angle now calculated is closest to real radioactive source actual angle.
Radioactive source to be measured in the embodiment of the present invention can be gamma-ray emission source.
The embodiment of the present invention provide radioactive source positioner, use crystal count be used for calculate parameter, make parameter with The change of radioactive source angle is obvious, uses and effectively deducts noise technique, for the radioactive source not feelings in crystal array plane Condition, removes unhelpful counting noise, it is possible to be effectively improved in parameter the changed factor relevant with radioactive source angle effectively, it is to avoid The defect of prior art, improves sensitivity and the accuracy of radioactive source location.
Fig. 7 is the system block diagram of radioactive source alignment system embodiment of the present invention.As it is shown in fig. 7, the embodiment of the present invention provides A kind of radioactive source status system, including: radioactive source positioner 73, radioactive source to be measured 74 and detector 75.Radioactive source location dress Put 73, then to calculate according to this relation than the relation with radioactive source angle according to crystal count in analog data acquisition detector 75 The radioactive source angle of radioactive source 74 to be measured.
In native system embodiment in the function of radioactive source positioner 73 such as above-mentioned Fig. 5 and Fig. 6 shown device embodiment Specifically describe, do not repeat them here.
Radioactive source to be measured 74 in the embodiment of the present invention can be gamma-ray emission source.
The embodiment of the present invention provide radioactive source alignment system, use crystal count be used for calculate parameter, make parameter with The change of radioactive source angle is obvious, uses and effectively deducts noise technique, for the radioactive source not feelings in crystal array plane Condition, removes unhelpful counting noise, it is possible to be effectively improved in parameter the changed factor relevant with radioactive source angle effectively, it is to avoid The defect of prior art, improves sensitivity and the accuracy of radioactive source location.
In the various embodiments described above of the present invention, the sequence number of embodiment is merely convenient of description, does not represent the quality of embodiment.Right The description of each embodiment all emphasizes particularly on different fields, and does not has the part described in detail, may refer to the phase of other embodiments in certain embodiment Close and describe.
One of ordinary skill in the art will appreciate that: all or part of step realizing said method embodiment can be passed through The hardware that programmed instruction is relevant completes, and aforesaid program can be stored in a computer read/write memory medium, this program Upon execution, perform to include the step of said method embodiment;And aforesaid storage medium includes: read only memory (Read- Only Memory, is called for short ROM), random access memory (Random Access Memory, be called for short RAM), magnetic disc or light The various medium that can store program code such as dish.
In the embodiments such as apparatus and method of the present invention, it is clear that each parts or each step be can decompose, combine and/ Or reconfigure after decomposing.These decompose and/or reconfigure the equivalents that should be regarded as the present invention.Meanwhile, the most right In the description of the specific embodiment of the invention, describe for a kind of embodiment and/or the feature that illustrates can be with same or like Mode use in one or more other embodiment, combined with the feature in other embodiment, or substitute it Feature in its embodiment.
It should be emphasized that term " includes/comprises " existence referring to feature, key element, step or assembly herein when using, but also It is not excluded for the existence of one or more further feature, key element, step or assembly or additional.
Although it is last it is noted that the present invention and advantage thereof have been described in detail above it should be appreciated that not Beyond can carry out in the case of the spirit and scope of the present invention that are defined by the claims appended hereto various change, replacement and Conversion.And, the scope of the present invention is not limited only to the concrete reality of the process described by description, equipment, means, method and steps Execute example.One of ordinary skilled in the art will readily appreciate that from the disclosure, can use according to the present invention and hold The row function essentially identical to corresponding embodiment described herein or obtain the result essentially identical with it, existing and future Process, equipment, means, method or step to be developed.Therefore, appended claim is directed in the range of them wrapping Include such process, equipment, means, method or step.

Claims (8)

1. a radioactive source localization method, it is characterised in that including:
Step S1, according to analog data obtain detector crystal array in crystal count than the relation with radioactive source angle;
Step S2, measure the actual measurement crystal count of described crystal array during radioactive source emitting radiation to be measured;
Step S3, in the excursion of counting noise, obtain according to described actual measurement crystal count and each described counting noise Multiple denoising crystal count ratios, according to described crystal count than with the Relation acquisition of radioactive source angle and each described denoising crystal meter Number than corresponding each actual measurement radioactive source angle, and actual measurement radioactive source angle when choosing described counting noise optimum be described in treat Survey the radioactive source angle of radioactive source;Wherein,
Described crystal array is 2 × 2 crystal arrays, and described step S1 includes:
Step S11, in described 2 × 2 crystal arrays optional one piece of crystal as reference crystal;
Step S12, it is first crystal group by described reference crystal and one piece of crystal combinations being adjacent, by described reference crystal And another block crystal combinations being adjacent is the second crystal group;
Step S13, the first crystal utilizing simulation softward to measure described first crystal group count than AB and described second crystal group The second crystal count than BC, wherein, AB=(A+B)/(A+B+C+D), BC=(B+C)/(A+B+C+D), A, B, C and D are successively For the crystal count of four pieces of crystal in described 2 × 2 crystal arrays;
Step S14, obtain described first crystal counting than the discrete relationship of AB and described radioactive source angle, θ and described the respectively Two crystal count discrete relationship than BC Yu described radioactive source angle, θ;
Step S15, described discrete relationship is fitted to linear relationship.
Radioactive source localization method the most according to claim 1, it is characterised in that during described radioactive source emitting radiation to be measured, In described 2 × 2 crystal arrays, the actual measurement crystal count of four pieces of crystal is followed successively by A ', B ', C ' and D ', the change of described counting noise Scope is Nmin to Nmax, and described step S3 includes:
Step S31, described counting noise N is set to Nmin;
Step S32, obtain the first denoising crystal count ratio, root according to AB '=(A '+B '-2 × N)/(A '+B '+C '+D '-4 × N) The second denoising crystal count ratio is obtained according to BC '=(B '+C '-2 × N)/(A '+B '+C '+D '-4 × N);
Step S33, obtain, according to described linear relationship, the radioactive source angle, θ 2 that radioactive source angle, θ 1 and BC ' corresponding for AB ' is corresponding;
Step S34, described counting noise N are from adding 1;
Step S35, compare | θ 1-θ 2 | andSize,WithFor parameter preset, if Parameter is then setOtherwise, parameterWithConstant;
Step S36, judge that described counting noise N whether more than Nmax, the most then performs step S37, otherwise, performs step S32;
Step S37, basisObtain the radioactive source angle of described radioactive source to be measured.
Radioactive source localization method the most according to claim 1 and 2, it is characterised in that described linear relationship is θ=(AB- Y1)/x1, θ=(BC-y2)/x2, wherein, x1, x2, y1 and y2 are Straight Line Fitting Parameters.
4. a radioactive source positioner, it is characterised in that including:
Processing module, in the crystal array obtaining detector according to analog data, crystal count is than the pass with radioactive source angle System;
Measurement module, for measuring the actual measurement crystal count of described crystal array during radioactive source emitting radiation to be measured;
Computing module, in the excursion at counting noise, surveys according to each described counting noise and described measurement module The described actual measurement crystal count measured obtains multiple denoising crystal count ratios, the described crystal obtained according to described processing module Counting than with the Relation acquisition of radioactive source angle and each described denoising crystal count than corresponding each actual measurement radioactive source angle, and The radioactive source angle that actual measurement radioactive source angle is described radioactive source to be measured when choosing described counting noise optimum;Wherein,
Described crystal array is 2 × 2 crystal arrays, and described processing module includes:
Select submodule, for optional one piece of crystal in described 2 × 2 crystal arrays as reference crystal, by described with reference to brilliant Body and the one piece of crystal combinations being adjacent are first crystal group, by described reference crystal and another block crystal group of being adjacent It is combined into the second crystal group;
Measure submodule, count than AB and described second for the first crystal utilizing simulation softward to measure described first crystal group Second crystal count of crystal group is than BC, wherein, AB=(A+B)/(A+B+C+D), BC=(B+C)/(A+B+C+D), A, B, C and D is followed successively by the crystal count of four pieces of crystal in described 2 × 2 crystal arrays;
First processes submodule, for obtaining the described first crystal counting discrete pass than AB Yu described radioactive source angle, θ respectively System and described second crystal count discrete relationship than BC Yu described radioactive source angle, θ;
Matching submodule, the described discrete relationship got for processing submodule by described first fits to linear relationship.
Radioactive source positioner the most according to claim 4, it is characterised in that during described radioactive source emitting radiation to be measured, In described 2 × 2 crystal arrays, the actual measurement crystal count of four pieces of crystal is followed successively by A ', B ', C ' and D ', the change of described counting noise Scope is Nmin to Nmax, and described computing module includes:
Submodule is set, for described counting noise N is set to Nmin;
First calculating sub module, brilliant for obtaining the first denoising according to AB '=(A '+B '-2 × N)/(A '+B '+C '+D '-4 × N) Body counting ratio, obtains the second denoising crystal count ratio according to BC '=(B '+C '-2 × N)/(A '+B '+C '+D '-4 × N);
Second processes submodule, obtains radiation corresponding to AB ' for the described linear relationship according to described matching submodule matching The radioactive source angle, θ 2 that source angle, θ 1 and BC ' is corresponding;
From adding submodule, for described counting noise N is added 1 certainly;
Comparison sub-module, is used for comparing | θ 1-θ 2 | andSize,WithFor parameter preset, if Parameter is then setOtherwise, parameterWithConstant;
Judge submodule, described from whether adding described counting noise N that submodule draws more than Nmax for judging;
Second calculating sub module, is used for when described counting noise N is more than Nmax, according to Obtain described The radioactive source angle of radioactive source to be measured.
6. according to the radioactive source positioner described in claim 4 or 5, it is characterised in that the institute of described matching submodule matching Stating linear relationship is θ=(AB-y1)/x1, θ=(BC-y2)/x2, and wherein, x1, x2, y1 and y2 are Straight Line Fitting Parameters.
7. a radioactive source alignment system, it is characterised in that including: the radioactive source location as described in arbitrary in claim 4-6 Device, radioactive source to be measured and detector.
Radioactive source alignment system the most according to claim 7, it is characterised in that described radioactive source to be measured is that gamma ray is put Penetrate source.
CN201310602906.6A 2013-11-25 2013-11-25 Radioactive source localization method, Apparatus and system Active CN103645491B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310602906.6A CN103645491B (en) 2013-11-25 2013-11-25 Radioactive source localization method, Apparatus and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310602906.6A CN103645491B (en) 2013-11-25 2013-11-25 Radioactive source localization method, Apparatus and system

Publications (2)

Publication Number Publication Date
CN103645491A CN103645491A (en) 2014-03-19
CN103645491B true CN103645491B (en) 2016-12-07

Family

ID=50250739

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310602906.6A Active CN103645491B (en) 2013-11-25 2013-11-25 Radioactive source localization method, Apparatus and system

Country Status (1)

Country Link
CN (1) CN103645491B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104155681A (en) * 2014-07-21 2014-11-19 北京辛耕普华医疗科技有限公司 Locator and method for monitoring radioactive source
CN104155675A (en) * 2014-08-27 2014-11-19 中国科学院高能物理研究所 Radiation source positioning and imaging device
CN104330814B (en) * 2014-10-29 2017-06-20 西南科技大学 A kind of radioactive source localization method and system
CN104597068A (en) * 2015-01-08 2015-05-06 中国科学院高能物理研究所 Collimation structure of space initiative-excitation laser X fluorescent spectrometer
CN104965216B (en) * 2015-06-11 2017-07-07 山东航天电子技术研究所 The method that detector detection efficient based on radiographic source on-orbit calibration system is demarcated
CN105068109B (en) * 2015-07-28 2017-08-01 济南中威仪器有限公司 It is a kind of to be used to search the localization method for losing radioactive source
CN106483545B (en) * 2016-09-14 2021-12-24 中国科学院高能物理研究所 Radioactive source monitoring method and system
CN107064987B (en) * 2017-01-16 2019-10-22 北京科技大学 A kind of radioactive source positioning system and localization method
CN107436446B (en) * 2017-08-07 2020-06-19 西南科技大学 Radioactive source positioning method and positioning system
CN107884805B (en) * 2017-09-28 2019-08-23 苏州瑞派宁科技有限公司 A kind of method and device for penetrating source positioning
CN108445527B (en) * 2018-01-30 2020-04-14 张岚 Method for acquiring position coordinates of ray source
CN114371494B (en) * 2022-03-22 2022-06-24 西南科技大学 Radioactive source scene simulation method for autonomous sourcing robot

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN86106216A (en) * 1986-08-23 1988-03-02 株式会社岛津制作所 The stabilising arrangement of PMT amplification factor
CN1338048A (en) * 1999-01-29 2002-02-27 特罗克斯勒电子实验有限公司 Apparatus and method for calibration of nuclear gauges
WO2003007019A1 (en) * 2001-07-09 2003-01-23 National Institute Of Radiological Sciences Office Of Public Relations Radiation three-dimensional position detector
JP3795041B2 (en) * 2003-12-16 2006-07-12 三菱重工業株式会社 Radioactive substance content measuring method and measuring apparatus
CN201237645Y (en) * 2008-04-30 2009-05-13 宋跃进 On-line positioning and monitoring apparatus for radioactive source
CN101460864A (en) * 2006-09-19 2009-06-17 株式会社岛津制作所 Nuclear medicine diagnosis system
CN101688918A (en) * 2007-04-23 2010-03-31 独立行政法人放射线医学综合研究所 Utilize the radiation detecting method and the device of energy and positional information
CN102759408A (en) * 2011-04-25 2012-10-31 中国科学院空间科学与应用研究中心 Single-photon counting imaging system and method of same
CN102944890A (en) * 2012-11-06 2013-02-27 中国科学院高能物理研究所 PS-PMT (position sensitive-photomultiplier tube) based detector signal readout method and system
CN202815229U (en) * 2012-08-08 2013-03-20 北京辛耕普华医疗科技有限公司 Gamma radioactive source positioning device and gamma camera used for same

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4568818B2 (en) * 2005-09-29 2010-10-27 独立行政法人 日本原子力研究開発機構 Visualization device using gamma ray source
JP5588190B2 (en) * 2010-02-08 2014-09-10 日立アロカメディカル株式会社 Radiation detection module

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN86106216A (en) * 1986-08-23 1988-03-02 株式会社岛津制作所 The stabilising arrangement of PMT amplification factor
CN1338048A (en) * 1999-01-29 2002-02-27 特罗克斯勒电子实验有限公司 Apparatus and method for calibration of nuclear gauges
WO2003007019A1 (en) * 2001-07-09 2003-01-23 National Institute Of Radiological Sciences Office Of Public Relations Radiation three-dimensional position detector
JP3795041B2 (en) * 2003-12-16 2006-07-12 三菱重工業株式会社 Radioactive substance content measuring method and measuring apparatus
CN101460864A (en) * 2006-09-19 2009-06-17 株式会社岛津制作所 Nuclear medicine diagnosis system
CN101688918A (en) * 2007-04-23 2010-03-31 独立行政法人放射线医学综合研究所 Utilize the radiation detecting method and the device of energy and positional information
CN201237645Y (en) * 2008-04-30 2009-05-13 宋跃进 On-line positioning and monitoring apparatus for radioactive source
CN102759408A (en) * 2011-04-25 2012-10-31 中国科学院空间科学与应用研究中心 Single-photon counting imaging system and method of same
CN202815229U (en) * 2012-08-08 2013-03-20 北京辛耕普华医疗科技有限公司 Gamma radioactive source positioning device and gamma camera used for same
CN102944890A (en) * 2012-11-06 2013-02-27 中国科学院高能物理研究所 PS-PMT (position sensitive-photomultiplier tube) based detector signal readout method and system

Also Published As

Publication number Publication date
CN103645491A (en) 2014-03-19

Similar Documents

Publication Publication Date Title
CN103645491B (en) Radioactive source localization method, Apparatus and system
CN103163548B (en) Based on the radioactive substance detection of gamma camera and device thereof and system
CN104536029B (en) Based on multiple NaI(T1)The radioactive substance two-dimensional location method of detector
US20210173096A1 (en) Method and device for multi-dimensional direction measurement of gamma radiation in the far field
CN103901052B (en) A kind of SGS and TGS combined measurement device and collimating apparatus optimization method
Canadas et al. Validation of a small-animal PET simulation using GAMOS: a GEANT4-based framework
CN105496436A (en) Time correction method and device used for PET device
Jan et al. PSF reconstruction for Compton-based prompt gamma imaging
CN100399051C (en) Passive efficient graduating method of radiation detector
CN107195342A (en) A kind of method for measuring nuclear power plant's nucleic activity distribution
JP5523407B2 (en) Radiation detection apparatus and detection method
Xiao et al. A study on scattering correction for γ-photon 3D imaging test method
Wei et al. Calibration of instrument and sample parameters for small angle X-ray scattering
El Bitar et al. A detector response function design in pinhole SPECT including geometrical calibration
CN109581472A (en) The selection of the spectrometer characteristic energy peak road CZT location range and determining method
Li et al. A SPECT system simulator built on the SolidWorksTM 3D-Design package
CN108445029A (en) The method for improving Compton scatter imaging performance
Kowash A rotating modulation imager for the orphan source search problem
Dong et al. Development and preliminary results of a large-pixel two-layer labr3 compton camera prototype
Hetzel et al. Near-field coded-mask technique and its potential for proton therapy monitoring
Shirakawa et al. Remote sensing of nuclear accidents using a direction finding detector
Heath et al. Gantryless Associated-Particle Neutron Radiography
Niknami et al. Design and simulation of a semiconductor detector-based Compton imaging system with efficiency analysis
El Bitar et al. The performance of a hybrid analytical-Monte Carlo system response matrix in pinhole SPECT reconstruction
Toivonen et al. On the symmetry of photon detection arrays: A directionally sensitive 3D model

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant