CN103054604A - Method and device for converting energy information and position information of positron emission tomography (PET) system and PET system - Google Patents
Method and device for converting energy information and position information of positron emission tomography (PET) system and PET system Download PDFInfo
- Publication number
- CN103054604A CN103054604A CN2012105710582A CN201210571058A CN103054604A CN 103054604 A CN103054604 A CN 103054604A CN 2012105710582 A CN2012105710582 A CN 2012105710582A CN 201210571058 A CN201210571058 A CN 201210571058A CN 103054604 A CN103054604 A CN 103054604A
- Authority
- CN
- China
- Prior art keywords
- value
- object program
- intermediate object
- specially
- place
- 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
Images
Landscapes
- Nuclear Medicine (AREA)
Abstract
The invention discloses a method and a device for converting energy information and position information of a positron emission tomography (PET) system and a PET system. The method includes acquiring inputting of a y value, wherein the y value is the sum of output energy of a photomultiplier tube of a detection module; inquiring an association table to obtain a first intermediate result corresponding to the y value, and determining the precision of H/y according to a value which enables the difference of q value and a standard value (x*con)/y in a preset error threshold value, wherein the association table specifically is generated by associating first intermediate results corresponding to various y values in advance; acquiring inputting of an x value, wherein the x value is an energy value of the output energy of the photomultiplier tube of the detection module, and the energy value of the output energy of the photomultiplier tube of the detection module is used for determining an axial coordinate and is smaller than the sum of the output energy of the photomultiplier tube; and multiplying the first intermediate result corresponding to the y value by K times of x to obtain a temporary result, and enabling the temporary result represented by a binary system to perform right shift by t1 to obtain a q value, wherein digits of K and t1 are determined according to the value which enables the difference of the q value and the standard value (x*con)/y in the preset error threshold value, and q value is the axial coordinate.
Description
Technical field
The present invention relates to the PET system, particularly method, the Apparatus and system of the conversion of a kind of PET system capacity information and positional information.
Background technology
Positron emission tomography device (Positron emission tomography, PET) be present high-level nuclear medicine facility, by in human body, injecting with the isotope-labeled chemical compound of positron, the positron that labelling positron isotope discharges is behind motion one segment distance, bury in oblivion with the electronics with negative charge in the surrounding, survey a pair of opposite direction that positron-electron annihilation produces external by detecting module this moment, the γ photon that energy equates, obtain positional information and the energy information of annihilation event, and then obtain the tomography scattergram of positron radionuclide in the body, from image, find the position of focus, thereby reach the purpose of early diagnosis.
Wherein, detecting module is by 4 photomultiplier tubes (PMT, as shown in Figure 4, supposing to demarcate respectively is a, b, c and d), photoconduction and crystal composition.The output of 4 photomultiplier tubes is respectively Ea, Eb, and Ec, Ed supposes that the photomultiplier tube take demarcation as d is initial point, can obtain the position coordinates x1 of annihilation event by the relation between following each formula of Anger logic centroid algorithm, y1, and energy information E, Eac, Eab:
x1=(Eab×con)/E
y1=(Eac×con)/E
(0, con), (0, con), wherein con is constant to y1 ∈ to x1 ∈
E=Ea+Eb+Ec+Ed
Eac=Ea+Ec
Eab=Ea+Eb
The formula of calculating location information in the following formula can be expressed as with general conversion formula: q=(x * con)/y, q ∈ (0, con), wherein q is a wherein axial coordinate of the coordinate system of annihilation event take a certain photomultiplier tube of PET system detecting module as initial point, x is the energy value that is used for determining this axial coordinate in the detecting module photomultiplier tube output energy, for example above-mentioned Eac or, Eab, y are four photomultiplier tube energy summations.
Be converted to optical signal hit the crystal of detecting module when gamma-rays after, photoconduction is reallocated to optical signal, optical signal after the distribution is through behind the PMT, be converted to the signal of telecommunication, produce annihilation event, produce waveform through the RC integrating circuit, as shown in Figure 5, the annihilation event scan period, sampling time and conversion time formed by the integration rise time.Present conversion method is directly utilized above-mentioned conversion formula conversion, the time that equals a multiplication its conversion time adds the time of last division, because present conversion time is long, sampling time from the scan period end that causes the single sweep operation cycle away from, the voltage amplitude value is unstable, the data that collect are not accurate enough, and the single sweep operation overlong time also can cause losing of adjacent events, cause the sensitivity of system.
Summary of the invention
In view of this, main purpose of the present invention is to provide the method for a kind of PET system capacity information and positional information conversion, the purpose that Apparatus and system reduces the conversion time of PET system single sweep operation with realization.
The invention provides the method for a kind of PET system capacity information and positional information conversion, comprising:
Obtain the input of y value, wherein the y value is got detecting module photomultiplier tube output energy sum total;
The inquiry contingency table obtains the first intermediate object program corresponding to this y value, and described the first intermediate object program is
And
Precision according to making q value and standard value
Difference determine that in default error threshold described contingency table is specially the be associated table of generation of the first in advance that various y values are corresponding intermediate object program;
The input of acquisition x, described x value are got in the detecting module photomultiplier tube output energy and be used for determining the wherein energy value of an axial coordinate, and be described for determining that wherein the energy value of an axial coordinate is exported the energy sum total less than photomultiplier tube;
The K that the first intermediate object program that described y value is corresponding multiply by x doubly obtains interim result, so that with the described interim result of the binary representation t that moves to right
1The position obtains q value, wherein K and t
1Figure place is according to making q value and standard value
Difference determine that in default error threshold wherein the q value is described axial coordinate.
Alternatively, described H is specially
Value, wherein con is constant, the t that is determined by the AD figure place of PET system
0Get and make error
Figure place in default error threshold scope;
Described K is specially 1;
Described t
1Specifically equal t
0
Alternatively, described H is specially the value of con, and wherein con is the constant by the AD figure place decision of PET system;
Described K is specially
T wherein
0Get and make error
Figure place in default error threshold scope;
Described t
1Specifically equal t
0
Alternatively, described the first intermediate object program
Concrete according to binary representation
The result of calculation t that moves to right
2The position obtains, wherein said t
2The position is specifically according to this t
2Figure place makes q value and the standard value that obtains
Difference determine that in default error threshold wherein, con is constant, the t that the AD figure place by the PET system determines
0Get and make error
Figure place in default error threshold scope;
Described K is specially 1;
Described t
1Specifically equal t
0-t
2
Alternatively, described the first intermediate object program
Concrete according to binary representation
The result of calculation t that moves to right
2The first acquisition of this left number is removed until left several first place is 1 in position, the n position that moves to left again, and wherein, con is constant, the t by the AD figure place decision of PET system
0Get and make error
Figure place in default error threshold scope, t
2Specifically according to this t
2Figure place makes q value and the standard value that obtains
Difference in default error threshold, determine;
Described K is specially 1;
Described t
1Specifically equal t
0-t
2+ n.
Alternatively, described contingency table is specially in advance the table that the first intermediate object program that various y values are corresponding is associated and generates, and wherein 2
pThe first intermediate object program corresponding to times y value and y value all is associated with same the first intermediate object program, and wherein p is specially the integer more than or equal to 1.
Alternatively, obtain also to comprise after the input of x value that wherein the m position specifically makes t according to this m position with the x value m position that moves to left
0-t
2The value of+n+m is that steady state value is determined;
Described t
1Specifically equal t
0-t
2+ n+m.
Alternatively, described acquisition y value be input to the inquiry contingency table obtain the first intermediate object program corresponding to this y value specifically in the input that obtains x, execution in the m position execution that the x value moved to left.
Alternatively, described contingency table specifically can be stored in SRAM, RAM, Rom or the macroelement.
Alternatively, specifically adopt programmable logic array to realize described method.
Alternatively, described detecting module is specially the detecting module in the PET system.
The present invention also provides the device of a kind of PET system capacity information and positional information conversion, comprising:
The first input block: be used for obtaining the input of y value, wherein the y value is got detecting module photomultiplier tube output energy summation;
Query unit: be used for the inquiry contingency table, obtain the first intermediate object program corresponding to this y value, the first intermediate object program is specially
And
Precision according to making q value and standard value
Difference determine that in default error threshold described contingency table is specially the table that the first in advance that various y values are corresponding intermediate object program is associated and generates, described contingency table is specially in advance the be associated table of generation of the first intermediate object program with various y values correspondences;
The second input block: be used for obtaining the input of x value, described x value is got in the detecting module photomultiplier tube output energy and is used for determining the wherein energy value of an axial coordinate, and is described for determining that wherein the energy value of an axial coordinate is exported the energy sum total less than photomultiplier tube;
Converting unit: be used for the K that the first intermediate object program with described y value correspondence multiply by x and doubly obtain interim result, so that with the interim result of the binary representation t that moves to right
1The position obtains q value, wherein K and t
1Figure place is according to making q value and standard value
Difference determine that in default error threshold wherein the q value is described axial coordinate.
As seen the present invention has following beneficial effect:
Because the present invention is with formula of the prior art
Be improved to
Can obtain by the contingency table that inquiry generates in advance the first intermediate object program corresponding to y value
Multiplication and dextroposition can obtain the q value faster by arithmetic speed again, so that detecting module is when carrying out periodic scan, equal the time of tabling look-up the conversion time that within the single scan period, spends, multiplication time and shift time sum, (the Time Calculation time that x * con)/y adopts the time of a multiplication to add last division will lack so that utilize conversion formula q=when changing than prior art the conversion time of conversion method of the present invention, can be so that move behind the sampling point position, reach more accurately purpose of image data, scan period is shortened, be difficult for losing adjacent events, sensitivity improves.
Description of drawings
Fig. 1 is the method flow diagram of the present invention's a kind of PET system capacity information and positional information conversion;
Fig. 2 is the installation composition figure of the present invention's a kind of PET system capacity information and positional information conversion;
Fig. 3 is the composition diagram of a kind of PET of the present invention system;
Fig. 4 is the photomultiplier tube sketch map of the detecting module of a kind of PET of the present invention system;
Fig. 5 is the scan period sketch map that the PET system surveys annihilation event.
The specific embodiment
For above-mentioned purpose of the present invention, feature and advantage can be become apparent more, below in conjunction with the drawings and specific embodiments the embodiment of the invention is described in further detail.
Referring to Fig. 1, be the method flow diagram of a kind of PET system capacity information provided by the invention and positional information conversion, as shown in the figure, the method comprises:
The input of S101, acquisition y value, wherein the y value is got detecting module photomultiplier tube output energy summation;
For example: the y value is the photomultiplier tube output energy summation E=Ea+Eb+Ec+Ed of detecting module as shown in Figure 4
S102, inquiry contingency table obtain the first intermediate object program corresponding to this y value, and the first intermediate object program is specially
And
Precision according to making q value and standard value
Difference determine that in default error threshold described contingency table is specially the be associated table of generation of the first in advance that various y values are corresponding intermediate object program;
The input of S103, acquisition x value, described x value are got in the detecting module photomultiplier tube output energy and be used for determining the wherein energy value of an axial coordinate, and be described for determining that wherein the energy value of an axial coordinate is exported the energy sum total less than photomultiplier tube;
For example, the x value is photomultiplier tube output ENERGY E ac or the Eab of detecting module as shown in Figure 4, Eac=Ea+Ec wherein, Eab=Ea+Eb.
The K that S104, the first intermediate object program that described y value is corresponding multiply by x doubly obtains interim result, so that with the interim result of the binary representation t that moves to right
1The position obtains q value, wherein K and t
1The position is according to making q value and standard value
Difference determine that in default error threshold wherein the q value is described axial coordinate.
As seen, the method that the energy information that the present invention proposes and positional information are changed is so that detecting module when carrying out periodic scan, equal the time of tabling look-up the conversion time that within the single scan period, spends, multiplication time and shift time sum, (time of the multiplication that x * con)/y calculates adds the Time Calculation time of last division and will lack so that directly utilize conversion formula q=than prior art the conversion time of conversion method of the present invention, can be so that move behind the sampling point position, reach more accurately purpose of image data, scan period is shortened, be difficult for losing adjacent events, sensitivity improves.
Need to prove, in the present invention, with formula
Improvement for
The first intermediate object program that the y value is corresponding
Write in advance in the contingency table, after obtaining the input of y value, can inquire about contingency table and obtain the first intermediate object program
And then again by arithmetic speed faster multiplication and dextroposition can obtain the q value, simultaneously by control H, K and t
1Figure place, can make q value and standard value
Difference in default error threshold.
Below, illustrate H, K and t
1Several modes that can implement of value:
Mode one:
For example, described H is specifically as follows
Value, namely the first intermediate object program is specially
Wherein con is specially 1 by constant, described K that the AD figure place of PET system determines, based on this known conditions, obtains expression formula
Suppose t
1Specifically equal t
0, then error is
t
0Can get and make error
Figure place in default error threshold scope;
In this embodiment, can preset error threshold is 1, because the x value is eternal less than or equal to the y value, therefore, as long as constraint
Just can guarantee error
Less than 1, therefore can get
For example, certain PET system adopts 8 AD, and con=255 supposes that the y value is input as 3, then t
0Can get 10, suppose that the y value is input as 237, then t
0Can get 16.
Need to prove that in this embodiment, t will move to right with the described interim result of binary representation in step S104
1Position, wherein t
1Specifically can equal t
0Make this interim result and standard value
Difference in default error threshold, for example: certain PET system adopts 8 AD, and con=255 supposes that the y value is input as 3, t
0Get 10, the first intermediate object programs and be 87125 or be shown 10101010001010101 with binary form, interim result is 87125 * x, and then will move to right with 87125 * x result of calculation of binary representation 10 obtains the q value.
Mode two:
For example, described H is specifically as follows the value of con, and wherein con is specially 2t by constant, described K that the AD figure place of PET system determines
0+ 1, based on this known conditions, still can obtain expression formula
Therefore, t
1And t
0Value still can be in the manner described above one computational methods obtain, do not repeat them here.
Mode three:
In aforesaid way one of the present invention, described contingency table is specially the first in advance that various y values are corresponding intermediate object program
The value table that is associated and generates, wherein,
The value figure place larger, cause the memory space of contingency table still larger, in order to save memory space, the present invention further proposes:
Described the first intermediate object program
Concrete according to binary representation
The result of calculation t that moves to right
2The position obtains, wherein said t
2The position is specifically according to this t
2Figure place makes q value and the standard value that obtains
Difference determine that in default error threshold wherein, con is constant, the t that the AD figure place by the PET system determines
0Get and make error
Figure place in default error threshold scope;
Described K is specially 1;
T will move to right with the described interim result of binary representation described in the step S104
1The position obtains the q value and is specially the t that moves to right with the described interim result with binary representation
0-t
2Position, i.e. t
1Specifically equal t
0-t
2Make this q value and standard value
Difference in default error threshold.
For example: certain PET system adopts 8 AD, and con=255 supposes that the y value is input as 3, t
0Get 10,
Result of calculation is 87125 or is shown 10101010001010101 with binary form, then according to q value and standard value
Difference in default error threshold, determine t
2Figure place is 8, particularly this t in actual applications
2Figure place can be according to q value and standard value
The restrictive condition of difference in default error threshold calculate through test of many times in advance and obtain, the first intermediate object program is 340 or is shown 101010100 with binary form, interim result is 340 * x, will move to right with 340 * x result of calculation of binary representation 2 and obtain the q value.
As seen, use this embodiment of the invention, according to binary representation
The result of calculation t that moves to right
2The position obtains the first intermediate object program, can reduce the memory space of described contingency table, when calculating the q value, the figure place that moves to right is reduced t simultaneously
2Figure place is calculated the error of the q value that obtains still in default error threshold so that use this embodiment method.
Mode four:
In conjunction with above-described embodiment, can also further shorten the figure place of the second intermediate object program with the purpose of the memory space of the described contingency table of realization minimizing by following steps, specific implementation is:
Described the first intermediate object program
Concrete according to binary representation
The result of calculation t that moves to right
2The first acquisition of this left number is removed until left several first place is 1 in position, the n position that moves to left again, and wherein, con is constant, the t by the AD figure place decision of PET system
0Get and make error
Figure place in default error threshold scope, t
2Specifically according to this t
2Figure place makes q value and the standard value that obtains
Difference in default error threshold, determine;
Described K is specially 1;
T will move to right with the described interim result of binary representation described in the step S104
1The position obtains the q value, is specially the t that moves to right with the described interim result with binary representation
0-t
2+ n position, i.e. t
1Specifically equal t
0-t
2+ n makes this q value and standard value
Difference in default error threshold.
For example: certain PET system adopts 8 AD, and con=255 supposes that the y value is input as 3, t
0Get 10, t
2Get 8, with binary representation
The result of calculation t that moves to right
2The position is 101010100, then according to q value and standard value
Difference determine that in default error threshold the n figure place is 0, this n figure place can be according to q value and standard value in actual applications particularly
The restrictive condition of difference in default error threshold calculate through test of many times in advance and obtain, will be with binary representation
The result of calculation t that moves to right
2Position, the n position that moves to left again are until left several first place is 1, and removing the first result of acquisition of this left number is 01010100, and interim result still is 340 * x, will move to right with 340 * x result of calculation of binary representation 2 and obtain the q value.
As seen, use this embodiment of the invention, by will be with binary representation
The result of calculation t that moves to right
2Move to left the n position until left several first place is 1 behind the position again, remove the first first intermediate object program that obtains of this left number, can further reduce the memory space of described contingency table, simultaneously when calculating the q value, the figure place that moves to right is increased the n position, calculate the error of the q value that obtains still in default error threshold so that use this embodiment method.
Need to prove in addition, for intermediate object program corresponding to different y values in the various embodiments described above, 2
pIntermediate object program corresponding to times y value can be identical with intermediate object program corresponding to y value, in conjunction with the various embodiments described above:
Specifically can be in advance with 2 in the described contingency table
pThe first intermediate object program corresponding to times y value and y value all is associated with same the first intermediate object program, and wherein p is specially the integer more than or equal to 1, although 2
pThe first intermediate object program corresponding to times y value and y value is identical, but finally calculates that to obtain the q value not identical, and reason is that the figure place that is shifted among the step S104 can be with 2 in the various embodiments described above
pDifferent corresponding the changing of times y value and y value still guarantees in the error threshold scope q value of final acquisition.
As seen, use the quantity that this embodiment can reduce intermediate object program, reach the purpose of the memory space of the described contingency table of further minimizing.
From aforesaid way four as seen, conversion method according to the inventive method proposition, so that detecting module is when carrying out periodic scan, equal the time of tabling look-up, multiplication time and shift time sum the conversion time that within the single scan period, spends, in order further to reduce the conversion time that spends in the single scan period, the present invention also propose with time of step S104 a part of displacement with to look into contingency table and multiplication time parallel, to reduce the conversion time that is spent in single scanning week and to be reduced in the complexity that realizes in the programmable array, implementation method is analyzed as follows:
Because the described figure place t that moves to right
1, concrete relevant with the value of y value and con, the variation of y value and con value will cause the figure place t that moves to right
1To change thereupon, to cause in programmable array, realizing that method of the present invention is cumbersome, in order addressing this problem, the present invention proposes: because the x value is necessarily less than or equal to the y value, so the processing that the x value can be moved to left in advance, so that for different y values, the figure place t that in step S104, moves to right
1Can be identical, in programmable array, realize comparatively simply, step and the y value that x value is moved to left in advance process are inputted and the step executed in parallel of tabling look-up, and to reduce the conversion time that is spent in the single scan period, specifically comprise:
Obtain the input of x value at step S103 after, step S103 comprises that also wherein the m position specifically makes t according to this m position with the x value m position that moves to left
0-t
2The value of+n+m is that steady state value is determined;
T moves to right described described interim result with binary representation in step S104
1The position obtains the q value, is specially the t that moves to right with the described interim result with binary representation
0-t
2+ n+m position, i.e. t
1Specifically equal t
0-t
2+ n+m makes this q value and standard value
Difference in default error threshold.
In this embodiment, more preferably, step S101 obtain y value be input to that step S102 inquiry contingency table obtains that the first intermediate object program corresponding to this y value specifically obtains the input of x value at step S103, execution in the m position execution that the x value is moved to left, step and the input of y value of processing tabled look-up and multiplication step executed in parallel so that the x value moves to left in advance, with the conversion time that is spent in the single scan period of further minimizing.
For example: certain PET system adopts 8 AD, and con=255 supposes that the y value is input as 3, t
0Get 10, t
2Getting 8, n figure place is 0, then according to q value and standard value
Difference determine that in default error threshold the m figure place is 6, this m figure place can be according to q value and standard value in actual applications particularly
The restrictive condition of difference in default error threshold calculate through test of many times in advance and obtain, suppose that the x value moves to left to obtain x1 after 6, make t
0-t
2The value of+n+m is steady state value 8, and namely will move to right with 340 * x1 result of calculation of binary representation 8 obtains the q value.
Need to prove H, K and t
1Above-mentioned several value mode only be in order to make more several embodiment of proposing of easy to understand of the present invention, any modification of doing within the spirit and principles in the present invention all should fall within protection scope of the present invention.
Contingency table described in above-described embodiment specifically can be stored in SRAM, RAM, in Rom or the macroelement, if use macroelement, through with the prior art discovery of comparing, prior art PET system need to use general 1200 of macroelement, and the macroelement that the present invention uses probably only has 700, greatly reduced memory space, and because the minimizing of conversion time, two axial coordinates in take a certain photomultiplier tube of described detecting module as the coordinate system of initial point, the time difference that i.e. two axial q values calculating respectively obtain is only by multiplication time difference decision, and time difference reduces.
Need to prove that above-mentioned each embodiment method all can adopt programmable logic array to realize that the detecting module in above-described embodiment is specifically as follows the detecting module in the PET system.
Referring to Fig. 2, be the installation composition figure of a kind of PET system capacity information provided by the invention and positional information conversion, as shown in the figure, this device comprises:
The first input block 201: be used for obtaining the input of y value, wherein the y value is got detecting module photomultiplier tube output energy summation;
Query unit 202: be used for the inquiry contingency table, obtain the first intermediate object program corresponding to this y value, the first intermediate object program is specially
And
Precision according to making q value and standard value
Difference determine that in default error threshold described contingency table is specially the be associated table of generation of the first in advance that various y values are corresponding intermediate object program;
The second input block 203: be used for obtaining the input of x value, described x value is got in the detecting module photomultiplier tube output energy and is used for determining the wherein energy value of an axial coordinate, and is described for determining that wherein the energy value of an axial coordinate is exported the energy sum total less than photomultiplier tube;
Converting unit 204: be used for the K that the first intermediate object program with described y value correspondence multiply by x and doubly obtain interim result, so that with the interim result of the binary representation t that moves to right
1The position obtains q value, wherein K and t
1Figure place is according to making q value and standard value
Difference determine that in default error threshold wherein the q value is described axial coordinate.
Referring to Fig. 3, be a kind of PET provided by the invention system composition diagram, as shown in the figure, this system comprises:
The device 301 of above-mentioned a kind of PET system capacity information and positional information conversion: be used for obtaining the input of y value, wherein the y value is got detecting module photomultiplier tube output energy summation; The inquiry contingency table obtains the first intermediate object program corresponding to this y value, and the first intermediate object program is specially
And
Precision according to making q value and standard value
Difference determine that in default error threshold described contingency table is specially the be associated table of generation of the first in advance that various y values are corresponding intermediate object program; The input of acquisition x value, described x value are got in the detecting module photomultiplier tube output energy and be used for determining the wherein energy value of an axial coordinate, and be described for determining that wherein the energy value of an axial coordinate is exported the energy sum total less than photomultiplier tube; The K that the first intermediate object program that described y value is corresponding multiply by x doubly obtains interim result, so that with the interim result of the binary representation t that moves to right
1The position obtains q value, wherein K and t
1Figure place is according to making q value and standard value
Difference determine that in default error threshold wherein the q value is described axial coordinate;
Detecting module 302: be used for detection of photons, the energy of output photoelectric multiplier tube.
Need to prove, in this article, relational terms such as the first and second grades only is used for an entity or operation are made a distinction with another entity or operation, and not necessarily requires or hint and have the relation of any this reality or sequentially between these entities or the operation.And, term " comprises ", " comprising " or its any other variant are intended to contain comprising of nonexcludability, thereby not only comprise those key elements so that comprise process, method, article or the equipment of a series of key elements, but also comprise other key elements of clearly not listing, or also be included as the intrinsic key element of this process, method, article or equipment.Do not having in the situation of more restrictions, the key element that is limited by statement " comprising ... ", and be not precluded within process, method, article or the equipment that comprises described key element and also have other identical element.
The above is preferred embodiment of the present invention only, is not for limiting protection scope of the present invention.All any modifications of doing within the spirit and principles in the present invention, be equal to replacement, improvement etc., all be included in protection scope of the present invention.
Claims (13)
1. the method for a PET system capacity information and positional information conversion is characterized in that, comprising:
Obtain the input of y value, wherein the y value is got detecting module photomultiplier tube output energy sum total;
The inquiry contingency table obtains the first intermediate object program corresponding to this y value, and described the first intermediate object program is specially
And
Precision according to making q value and standard value
Difference determine that in default error threshold described contingency table is specially the be associated table of generation of the first in advance that various y values are corresponding intermediate object program;
The input of acquisition x value, described x value are got in the detecting module photomultiplier tube output energy and be used for determining the wherein energy value of an axial coordinate, and be described for determining that wherein the energy value of an axial coordinate is exported the energy sum total less than photomultiplier tube;
The K that the first intermediate object program that described y value is corresponding multiply by x doubly obtains interim result, so that with the interim result of the binary representation t that moves to right
1The position obtains q value, wherein K and t
1Figure place is according to making q value and standard value
Difference determine that in default error threshold wherein the q value is described axial coordinate.
2. method according to claim 1 is characterized in that, described H is specially
Value, wherein con is constant, the t that is determined by the AD figure place of PET system
0Get and make error
Figure place in default error threshold scope;
Described K is specially 1;
Described t
1Specifically equal t
0
3. method according to claim 1 is characterized in that, described H is specially con, and wherein con is the constant by the AD figure place decision of PET system;
Described K is specially
T wherein
0Get and make error
Figure place in default error threshold scope;
Described t
1Specifically equal t
0
4. method according to claim 1 is characterized in that, described the first intermediate object program
Concrete according to binary representation
The result of calculation t that moves to right
2The position obtains, wherein said t
2The position is specifically according to this t
2Figure place makes q value and the standard value that obtains
Difference determine that in default error threshold wherein, con is constant, the t that the AD figure place by the PET system determines
0Get and make error
Figure place in default error threshold scope;
Described K is specially 1;
Described t
1Specifically equal t
0-t
2
5. method according to claim 1 is characterized in that, described the first intermediate object program
Concrete according to binary representation
The result of calculation t that moves to right
2The first acquisition of this left number is removed until left several first place is 1 in position, the n position that moves to left again, and wherein, con is constant, the t by the AD figure place decision of PET system
0Get and make error
Figure place in default error threshold scope, t
2Specifically according to this t
2Figure place makes q value and the standard value that obtains
Difference in default error threshold, determine;
Described K is specially 1;
Described t
1Specifically equal t
0-t
2+ n.
6. method according to claim 1 is characterized in that, described contingency table is specially in advance the table that the first intermediate object program that various y values are corresponding is associated and generates, wherein 2
pThe first intermediate object program corresponding to times y value and y value all is associated with same the first intermediate object program, and wherein p is specially the integer more than or equal to 1.
7. method according to claim 5 is characterized in that, obtains also to comprise after the input of x value that wherein the m position specifically makes t according to this m position with the x value m position that moves to left
0-t
2The value of+n+m is that steady state value is determined;
Described t
1Specifically equal t
0-t
2+ n+m.
8. method according to claim 7 is characterized in that, described acquisition y value be input to the inquiry contingency table obtain the first intermediate object program corresponding to this y value specifically in the input that obtains x, execution in the m position execution that the x value moved to left.
9. method according to claim 1, described contingency table specifically can be stored in SRAM, RAM, Rom or the macroelement.
10. method according to claim 1 is characterized in that, specifically adopts programmable logic array to realize described method.
11. method according to claim 1 is characterized in that, described detecting module is specially the detecting module in the PET system.
12. the device of a PET system capacity information and positional information conversion is characterized in that, comprising:
The first input block: be used for obtaining the input of y value, wherein the y value is got detecting module photomultiplier tube output energy summation;
Query unit: be used for the inquiry contingency table, obtain the first intermediate object program corresponding to this y value, the first intermediate object program is specially
And
Precision according to making q value and standard value
Difference determine that in default error threshold described contingency table is specially the be associated table of generation of the first in advance that various y values are corresponding intermediate object program;
The second input block: be used for obtaining the input of x value, described x value is got in the detecting module photomultiplier tube output energy and is used for determining the wherein energy value of an axial coordinate, and is described for determining that wherein the energy value of an axial coordinate is exported the energy sum total less than photomultiplier tube;
Converting unit: be used for the K that the first intermediate object program with described y value correspondence multiply by x and doubly obtain interim result, so that with the interim result of the binary representation t that moves to right
1The position obtains q value, wherein K and t
1Figure place is according to making q value and standard value
Difference determine that in default error threshold wherein the q value is described axial coordinate.
13. a PET system is characterized in that, comprising:
The device of a kind of PET system capacity information and positional information conversion: be used for obtaining the input of y value, wherein the y value is got detecting module photomultiplier tube output energy summation; The inquiry contingency table obtains the first intermediate object program corresponding to this y value, and the first intermediate object program is specially
And
Precision according to making q value and standard value
Difference determine that in default error threshold described contingency table is specially the be associated table of generation of the first in advance that various y values are corresponding intermediate object program; The input of acquisition x value, described x value are got in the detecting module photomultiplier tube output energy and be used for determining the wherein energy value of an axial coordinate, and be described for determining that wherein the energy value of an axial coordinate is exported the energy sum total less than photomultiplier tube; The K that the first intermediate object program that described y value is corresponding multiply by x doubly obtains interim result, so that obtain q value, wherein K and t with the interim result of the binary representation t1 position that moves to right
1Figure place is according to making q value and standard value
Difference determine that in default error threshold wherein the q value is described axial coordinate;
Detecting module: be used for detection of photons, the energy of output photoelectric multiplier tube.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210571058.2A CN103054604B (en) | 2012-12-25 | 2012-12-25 | Method and device for converting energy information and position information of positron emission tomography (PET) system and PET system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210571058.2A CN103054604B (en) | 2012-12-25 | 2012-12-25 | Method and device for converting energy information and position information of positron emission tomography (PET) system and PET system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103054604A true CN103054604A (en) | 2013-04-24 |
CN103054604B CN103054604B (en) | 2015-04-01 |
Family
ID=48097721
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210571058.2A Active CN103054604B (en) | 2012-12-25 | 2012-12-25 | Method and device for converting energy information and position information of positron emission tomography (PET) system and PET system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103054604B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100032574A1 (en) * | 2007-04-23 | 2010-02-11 | National Institute Of Radiological Sciences | Radiation detecting method utilizing energy information and positional information and equipment thereof |
EP2246712A2 (en) * | 2009-04-16 | 2010-11-03 | Kabushiki Kaisha Toshiba | Positron emission tomography apparatus and nuclear medical image generating method |
CN201775635U (en) * | 2010-08-19 | 2011-03-30 | 上海生物医学工程研究中心 | Detecting device for positron emission fault imaging device |
CN102262238A (en) * | 2011-04-19 | 2011-11-30 | 苏州瑞派宁科技有限公司 | Method and device for extracting scintillation pulse information |
-
2012
- 2012-12-25 CN CN201210571058.2A patent/CN103054604B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100032574A1 (en) * | 2007-04-23 | 2010-02-11 | National Institute Of Radiological Sciences | Radiation detecting method utilizing energy information and positional information and equipment thereof |
EP2246712A2 (en) * | 2009-04-16 | 2010-11-03 | Kabushiki Kaisha Toshiba | Positron emission tomography apparatus and nuclear medical image generating method |
CN201775635U (en) * | 2010-08-19 | 2011-03-30 | 上海生物医学工程研究中心 | Detecting device for positron emission fault imaging device |
CN102262238A (en) * | 2011-04-19 | 2011-11-30 | 苏州瑞派宁科技有限公司 | Method and device for extracting scintillation pulse information |
Also Published As
Publication number | Publication date |
---|---|
CN103054604B (en) | 2015-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
von Buddenbrock et al. | The emergence of multi-lepton anomalies at the LHC and their compatibility with new physics at the EW scale | |
Hesterman et al. | Maximum-likelihood estimation with a contracting-grid search algorithm | |
CN101600972B (en) | Time of flight measurements in positron emission tomography | |
RU2431157C2 (en) | Integrated multichannel transducer of time intervals in code for time-of-flight positron-emission tomography | |
CN103745185B (en) | Method and apparatus to recognize the position of detector crystal unit | |
CN102631212B (en) | PET (positron emission tomography) scanner and meeting wherein sentence choosing method | |
CN103460072A (en) | Detector array with time-to-digital conversion having improved temporal accuracy | |
CN101561507B (en) | Ray energy detection method for ray detector | |
CN102981685A (en) | Touch screen apparatus and driving method thereof | |
CN104035122A (en) | Method and device for correcting energy value | |
Erler et al. | Z′ bosons at colliders: a bayesian viewpoint | |
Morrison et al. | Multistage linear feedback shift register counters with reduced decoding logic in 130-nm CMOS for large-scale array applications | |
CN105824817A (en) | Flash pulse digitization method | |
CN105361901A (en) | Method and system for correcting depth effect of positron emission tomography | |
CN102981179B (en) | Position table generating method for scintillation detector | |
CN104184473A (en) | Method and device for carrying out non-linear correction on TDC | |
CN102944890A (en) | PS-PMT (position sensitive-photomultiplier tube) based detector signal readout method and system | |
CN105093262A (en) | Digital neutron spectrometer response matrix obtaining method | |
DeWitt et al. | Design of an FPGA-based algorithm for real-time solutions of statistics-based positioning | |
CN104932000B (en) | The method and device that γ photons positional information obtains in a kind of PET device | |
CN104382611B (en) | Correction method and device of output energy of photomultiplier | |
CN106202865B (en) | A method of calculating arbitrary order coefficient in the discrete locking nub method of neutron transport | |
CN109685866A (en) | Crystal positions look-up table construction method, device, computer equipment and storage medium | |
CN103054604B (en) | Method and device for converting energy information and position information of positron emission tomography (PET) system and PET system | |
CN106859686A (en) | Imaging method and imaging system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |