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 PDF

Info

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
Application number
CN2012105710582A
Other languages
Chinese (zh)
Other versions
CN103054604B (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.)
SHENYANG NEUSOFT PAISITONG MEDICAL SYSTEM CO Ltd
Original Assignee
SHENYANG NEUSOFT PAISITONG MEDICAL SYSTEM CO Ltd
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 SHENYANG NEUSOFT PAISITONG MEDICAL SYSTEM CO Ltd filed Critical SHENYANG NEUSOFT PAISITONG MEDICAL SYSTEM CO Ltd
Priority to CN201210571058.2A priority Critical patent/CN103054604B/en
Publication of CN103054604A publication Critical patent/CN103054604A/en
Application granted granted Critical
Publication of CN103054604B publication Critical patent/CN103054604B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

Method, the Apparatus and system of a kind of PET system capacity information and positional information conversion
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
Figure BDA00002651237700022
Precision according to making q value and standard value
Figure BDA00002651237700023
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
Figure BDA00002651237700024
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 BDA00002651237700031
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
Figure BDA00002651237700032
T wherein 0Get and make error
Figure BDA00002651237700033
Figure place in default error threshold scope;
Described t 1Specifically equal t 0
Alternatively, described the first intermediate object program
Figure BDA00002651237700034
Concrete according to binary representation
Figure BDA00002651237700035
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
Figure BDA00002651237700036
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 BDA00002651237700037
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
Figure DEST_PATH_GDA00002880332500038
Concrete according to binary representation
Figure DEST_PATH_GDA00002880332500039
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 DEST_PATH_GDA000028803325000310
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
Figure BDA00002651237700041
And
Figure BDA00002651237700042
Precision according to making q value and standard value
Figure BDA00002651237700043
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
Figure BDA00002651237700044
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
Figure BDA00002651237700051
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
Figure BDA00002651237700061
And
Figure BDA00002651237700062
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
Figure BDA00002651237700064
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
Figure BDA00002651237700065
Improvement for The first intermediate object program that the y value is corresponding
Figure BDA00002651237700067
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
Figure BDA00002651237700068
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
Figure BDA00002651237700072
Value, namely the first intermediate object program is specially
Figure BDA00002651237700073
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
q = x × con y = ( H y × x × K ) / 2 t 1 = ( ( 2 t 0 + 1 ) × con y × x ) / 2 t 1 x × con 2 t 0 y × 2 t 1 + x × con y × 2 t 1 ;
Suppose t 1Specifically equal t 0, then error is
Figure BDA00002651237700075
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
Figure BDA00002651237700077
Just can guarantee error
Figure BDA00002651237700078
Less than 1, therefore can get
Figure BDA00002651237700079
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
Figure BDA000026512377000710
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 q = x × con y = ( H y × x × K ) / 2 t 1 = ( ( 2 t 0 + 1 ) × con y × x ) / 2 t 1 x × con 2 t 0 y × 2 t 1 + x × con y × 2 t 1 , 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
Figure BDA00002651237700082
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
Figure BDA00002651237700085
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
Figure BDA00002651237700086
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 BDA00002651237700087
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
Figure BDA00002651237700088
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,
Figure BDA00002651237700089
Result of calculation is 87125 or is shown 10101010001010101 with binary form, then according to q value and standard value
Figure BDA000026512377000810
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
Figure BDA00002651237700091
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
Figure BDA00002651237700092
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
Figure BDA00002651237700093
Concrete according to binary representation
Figure BDA00002651237700094
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
Figure BDA00002651237700096
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
Figure BDA00002651237700097
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
Figure BDA00002651237700101
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
Figure BDA00002651237700103
The restrictive condition of difference in default error threshold calculate through test of many times in advance and obtain, will be with binary representation
Figure BDA00002651237700104
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
Figure BDA00002651237700105
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
Figure BDA00002651237700112
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
Figure BDA00002651237700121
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
Figure BDA00002651237700122
And
Figure BDA00002651237700123
Precision according to making q value and standard value
Figure BDA00002651237700124
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
Figure BDA00002651237700131
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
Figure BDA00002651237700132
And
Figure BDA00002651237700133
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
Figure BDA00002651237700135
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
Figure FDA00002651237600011
And
Figure FDA00002651237600012
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
Figure FDA00002651237600014
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
Figure FDA00002651237600015
Value, wherein con is constant, the t that is determined by the AD figure place of PET system 0Get and make error
Figure FDA00002651237600016
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
Figure FDA00002651237600017
T wherein 0Get and make error
Figure FDA00002651237600018
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
Figure FDA00002651237600019
Concrete according to binary representation
Figure FDA00002651237600021
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
Figure FDA00002651237600022
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 FDA00002651237600023
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
Figure FDA00002651237600024
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 FDA00002651237600026
Figure place in default error threshold scope, t 2Specifically according to this t 2Figure place makes q value and the standard value that obtains
Figure FDA00002651237600027
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
Figure FDA00002651237600031
And Precision according to making q value and standard value
Figure FDA00002651237600033
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
Figure FDA00002651237600034
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
Figure FDA00002651237600035
And
Figure FDA00002651237600036
Precision according to making q value and standard value
Figure FDA00002651237600041
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
Figure FDA00002651237600042
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.
CN201210571058.2A 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 Active CN103054604B (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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