WO2009141861A1 - 2次元位置マップ校正方法および放射線検出装置 - Google Patents
2次元位置マップ校正方法および放射線検出装置 Download PDFInfo
- Publication number
- WO2009141861A1 WO2009141861A1 PCT/JP2008/059146 JP2008059146W WO2009141861A1 WO 2009141861 A1 WO2009141861 A1 WO 2009141861A1 JP 2008059146 W JP2008059146 W JP 2008059146W WO 2009141861 A1 WO2009141861 A1 WO 2009141861A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- peaks
- dimensional position
- position map
- boundary
- peak
- Prior art date
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims description 21
- 238000000926 separation method Methods 0.000 claims abstract description 49
- 238000001514 detection method Methods 0.000 claims abstract description 20
- 230000035945 sensitivity Effects 0.000 claims description 37
- 230000003287 optical effect Effects 0.000 claims description 11
- 230000005251 gamma ray Effects 0.000 description 26
- 238000002600 positron emission tomography Methods 0.000 description 14
- 239000013078 crystal Substances 0.000 description 10
- 238000009206 nuclear medicine Methods 0.000 description 7
- 238000003745 diagnosis Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 229940121896 radiopharmaceutical Drugs 0.000 description 3
- 239000012217 radiopharmaceutical Substances 0.000 description 3
- 230000002799 radiopharmaceutical effect Effects 0.000 description 3
- 239000000284 extract Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 241000682622 Andracantha sigma Species 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
- G01T1/164—Scintigraphy
- G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
- G01T1/1644—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras using an array of optically separate scintillation elements permitting direct location of scintillations
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
- G01T1/164—Scintigraphy
- G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
- G01T1/1647—Processing of scintigraphic data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T7/00—Details of radiation-measuring instruments
- G01T7/005—Details of radiation-measuring instruments calibration techniques
Definitions
- the present invention relates to a two-dimensional position map calibration method for calibrating a two-dimensional position map used when detecting radiation with a radiation detector comprising a plurality of scintillator elements and an optical sensor optically coupled thereto. And a radiation detection apparatus.
- a PET (Positron Emission Tomography) apparatus As a nuclear medicine diagnostic apparatus, that is, an ECT (Emission Computed Tomography) apparatus, a PET (Positron Emission Tomography) apparatus will be described as an example.
- the PET apparatus detects a plurality of gamma rays generated by annihilation of positrons, that is, positrons, and reconstructs a tomographic image of a subject only when the gamma rays are simultaneously detected by a plurality of detectors. It is configured.
- a radiopharmaceutical containing a positron emitting nuclide is administered into a subject, and a 511 KeV pair annihilation gamma ray released from the administered subject consists of a group of a number of detection elements (for example, scintillators). Detect with a detector. And if ⁇ -rays are detected at the same time by two detectors within a certain period of time, they are counted as a pair of annihilation ⁇ -rays, and the point of occurrence of pair annihilation is on the straight line of the detected detector pair Is identified. By accumulating such coincidence information and performing reconstruction processing, a positron emitting nuclide distribution image (ie, a tomographic image) is obtained.
- a positron emitting nuclide distribution image ie, a tomographic image
- the number of scintillators is increased in order to obtain a more detailed ⁇ -ray incident position at the detector, combined with a photomultiplier tube (PMT: Photo Multiplier Tube) that can detect the position, Discrimination as a scintillator element improves the detection accuracy of ⁇ rays and improves the image resolution of tomographic images. Therefore, the discrimination capability is increased by increasing the number of scintillators.
- DOI detectors have been developed that can discriminate the light source position (DOI: Depth of Interaction) in the depth direction by causing scintillators to be stacked also in the depth direction.
- the two-dimensional position map is a two-dimensional coordinate (X, Y) for an event in which each ⁇ -ray is detected by calculating the center of gravity of an electrical signal obtained by an optical sensor represented by a position detection type photomultiplier tube. It is drawn by performing the work to calculate as. Further, this two-dimensional position map is obtained by irradiating the detector with uniform parallel beam ⁇ -rays and repeating the above-described operation while detecting the ⁇ -rays and integrating the two-dimensional coordinates on the two-dimensional plane. . These are drawn as a distribution having a peak corresponding to each scintillator element (crystal element) position.
- FIG. 9 is a two-dimensional position map in the case of a DOI detector in which four layers of scintillators are stacked in the depth direction, and the position indicated by a white circle (indicated by “ ⁇ ” in FIG. 9) is the first layer (in FIG. 9).
- LUT Look Up Table
- one detector is composed of several hundred to several thousand crystals (scintillator elements), and the same number of peaks appear on the image. Entering the boundary of the two-dimensional position map manually requires a lot of labor. Therefore, it is desirable to define the boundary by automation, but the peak is not sufficiently separated around the two-dimensional position map, and the boundary is often erroneously determined.
- the present invention has been made in view of such circumstances, and an object of the present invention is to provide a two-dimensional position map calibration method and a radiation detection apparatus that can easily determine the boundary.
- the two-dimensional position map calibration method of the present invention comprises a radiation detector composed of a plurality of scintillator elements arranged in one, two, or three dimensions and an optical sensor optically coupled to them.
- a two-dimensional position map that calibrates a two-dimensional position map that is used when detecting the two-dimensional position map in which the signal intensity obtained by the optical sensor is associated with the incident position of the radiation incident on the scintillator element.
- a calibration method wherein a boundary is drawn by the peak of the signal intensity, and each position is separated by the boundary, and a plurality of the peaks are connected to each other, and the separation in the peak separation process has failed.
- a number determination step for determining the number of peaks using the spatial periodicity of the peaks.
- the peak separation step a boundary is drawn by the peak of the signal intensity, and each position is separated by the boundary. If a plurality of peaks are connected and the separation in the above-described peak separation step fails, the number of failed peaks is determined in the number determination step using the spatial periodicity of the peaks. Therefore, the number of failed peaks can be determined by using the peak spatial periodicity, and the boundary can be easily determined.
- each of the scintillator elements that failed to be separated in the above-described peak separation step is determined by setting the boundary so that the sensitivity ratio for each scintillator element matches the sum ratio of the pixels in the peak region. It is preferable to provide a boundary determination step for separating the positions. Utilizing the fact that the sensitivity ratio and the pixel are in a proportional relationship, if the sensitivity ratio for each scintillator element is determined in advance, the scintillator is used to separate each position where the separation in the peak separation process has failed. The boundary can be determined so that the sensitivity ratio for each element matches the sum ratio of pixels in the peak region.
- the signal strengths are compared to obtain respective local minimum values, the positions of the local minimum values are drawn as boundaries, and the respective positions may be separated by the boundaries.
- Each local maximum value may be obtained by comparison, and the position of each local maximum value may be drawn as a boundary, and each position may be separated by the boundary.
- the peak is considered as the maximum, the position of the minimum value almost coincides with the boundary, so that the position of the minimum value is drawn as a boundary as in the former, and each position is separated by the boundary.
- the radiation detection apparatus of the present invention includes a radiation detector including a radiation detector composed of a plurality of scintillator elements arranged one-dimensionally, two-dimensionally or three-dimensionally, and an optical sensor optically coupled thereto.
- a two-dimensional position map representing a two-dimensional position map in which the signal intensity obtained by the optical sensor is two-dimensionally associated with the incident position of the radiation incident on the scintillator element.
- Storage means for storing a table in which each position is associated with each scintillator element; and arithmetic processing means for performing arithmetic processing for calibrating the two-dimensional position map, the calibrated two-dimensional position map;
- the radiation detection position is identified by discriminating the incident position based on the detection result of the radiation, and the arithmetic processing means determines the peak of the signal intensity.
- a peak separation step of drawing a boundary and separating each position by the boundary, and the number of peaks that failed to be separated in the peak separation step by connecting a plurality of the peaks, and the spatial periodicity of the peak And a number determination step that is used and determined, and performing arithmetic processing related to those steps.
- the arithmetic processing means draws a boundary by the peak of the signal intensity and separates each position by the boundary, and the peak separation step by connecting a plurality of peaks. And a number determination step of determining the number of peaks that failed to be separated using the spatial periodicity of the peaks, and performing arithmetic processing related to those steps. If a plurality of peaks are connected and the separation in the above-described peak separation step fails, the number of failed peaks is determined using the peak spatial periodicity in the number determination step. The number of failed peaks can be determined by using periodicity, and the boundary can be easily determined. As a result, the incident position can also be easily distinguished, and the radiation detection position can be easily specified.
- the boundary is drawn by the peak of the signal intensity, and each peak is separated by the boundary, and a plurality of peaks are connected.
- a number determination step of determining the number of peaks that failed to be separated in the peak separation step by using the spatial periodicity of the peaks, and performing arithmetic processing related to those steps. If a plurality of peaks are connected and the separation in the above-described peak separation step fails, the number of failed peaks is determined using the peak spatial periodicity in the number determination step. The number of failed peaks can be determined by using periodicity, and the boundary can be easily determined.
- FIG. 1 is a side view and block diagram of a PET (Positron Emission Tomography) apparatus according to an embodiment. It is a schematic perspective view of a gamma ray detector.
- A) is a top view of a gamma ray detector
- (b) is a side view of a gamma ray detector.
- (A) is a plan view of a two-dimensional position map used for explanation when obtaining a minimum value along the row (x) direction, and (b) is used for explanation when obtaining a minimum value along the column (y) direction. It is a top view of a two-dimensional position map.
- (A) is a plan view of a two-dimensional position map schematically showing a peak
- (b) is a plan view of a two-dimensional position map expressing a sensitivity ratio for each scintillator element representing the spatial periodicity of the peak
- (A) is a plan view of a two-dimensional position map schematically showing a peak when the target region is separated when the separation fails
- (b) is a separation of the target region when the separation fails. It is a top view of the two-dimensional position map which expressed the sensitivity ratio for every scintillator element showing the spatial periodicity of the peak in the case of doing.
- FIG. 1 is a side view and block diagram of a PET (Positron Emission Tomography) apparatus according to an embodiment
- FIG. 2 is a schematic perspective view of a ⁇ -ray detector
- FIG. 3A is a ⁇ -ray detector
- FIG. 3B is a side view of the ⁇ -ray detector.
- the PET apparatus includes a top plate 1 on which a subject M is placed as shown in FIG.
- the top plate 1 is configured to move up and down and translate along the body axis Z of the subject M.
- the subject M placed on the top 1 is scanned from the head to the abdomen and foot sequentially through the opening 2a of the gantry 2, which will be described later. Get the image. Note that there is no particular limitation on the scanned part and the scanning order of each part.
- the PET apparatus includes a gantry 2 having an opening 2a and a ⁇ -ray detector 3.
- the ⁇ -ray detector 3 is arranged in a ring shape so as to surround the body axis Z of the subject M, and is embedded in the gantry 2.
- the ⁇ -ray detector 3 corresponds to the radiation detector in the present invention.
- the PET apparatus includes a table driving unit 4, a controller 5, an input unit 6, an output unit 7, a memory unit 8, a position calculation circuit 9, a look-up table 10, a coincidence circuit 11, A configuration unit 12 and a two-dimensional position map calibration unit 13 are provided.
- the top plate driving unit 6 is a mechanism for driving the top plate 1 so as to perform the above-described movement, and is configured by a motor or the like not shown.
- the look-up table 10 corresponds to the storage means in the present invention
- the two-dimensional position map calibration unit 13 corresponds to the arithmetic processing means in the present invention, and includes the gantry 2, the ⁇ -ray detector 3, the position calculation circuit 9, and the lookup.
- the table 10, the coincidence counting circuit 11, and the two-dimensional position map calibration unit 13 constitute a radiation detection apparatus according to the present invention.
- the controller 5 comprehensively controls each part constituting the PET apparatus according to the present embodiment.
- the controller 5 includes a central processing unit (CPU).
- the input unit 6 sends data and commands input by the operator to the controller 5.
- the input unit 6 includes a pointing device represented by a mouse, a keyboard, a joystick, a trackball, a touch panel, and the like.
- the output unit 7 includes a display unit represented by a monitor, a printer, and the like.
- the memory unit 8 and the lookup table 10 are composed of storage media represented by ROM (Read-only Memory), RAM (Random-Access Memory), and the like.
- the count value (count) simultaneously counted by the coincidence circuit 11 and the image processed by the reconstruction unit 12 are written and stored in the RAM, and are read from the RAM as necessary.
- the center of gravity of an electric signal obtained by a photomultiplier tube 33 (see FIGS. 2 and 3) described later is calculated, and a scintillator of a scintillator block 31 (see FIGS. 2 and 3) described later is used.
- a two-dimensional position map expressed in a two-dimensional manner corresponding to the element position is written and stored in the lookup table 10 as a table in which each position in the two-dimensional position map is associated with each scintillator element.
- the two-dimensional position map is calibrated by the position map calibrating unit 13, it is read from the lookup table 10 and rewritten and calibrated.
- the ROM stores in advance a program for performing various nuclear medicine diagnosis and arithmetic processing by the two-dimensional position map calibration unit 13, and the controller 5 executes the program to execute the nuclear medicine corresponding to the program.
- the calculation process by the diagnosis and the two-dimensional position map calibration unit 13 is performed.
- the reconstruction unit 12 and the two-dimensional position map calibration unit 13 are input by a program stored in a ROM of a storage medium represented by the memory unit 8 or the like described above, or by a pointing device represented by the input unit 6 or the like, for example. This is realized by the controller 5 executing the instruction.
- the ⁇ -rays generated from the subject M to which the radiopharmaceutical is administered are converted into light by the scintillator block 31 (see FIGS. 2 and 3) of the ⁇ -ray detector 3, and the converted light is converted into the ⁇ -ray detector.
- the photomultiplier tube 3 (PMT: Photo Multiplier Tube) 32 (see FIGS. 2 and 3) is multiplied and converted into an electric signal.
- the electric signal is sent to the position calculation circuit 9 as image information (pixel value, that is, a count value simultaneously counted by the ⁇ -ray detector 3).
- the position calculation circuit 9 refers to the look-up table 10 at the time of nuclear medicine diagnosis, and refers to the two-dimensional position map to determine which scintillator element of the scintillator block 31 (see FIGS. 2 and 3) has the counted value. It is determined whether it is incident at. Specifically, the center of gravity is calculated for each incident to obtain the incident position on the scintillator element. The obtained incident position and count value (image information) are sent to the coincidence counting circuit 11.
- the coincidence circuit 11 checks the position of the scintillator block 31 (see FIG. 2 and FIG. 3) (more specifically, the incident position on the scintillator element) and the incident timing of the ⁇ -ray, and 2 on both sides of the subject M. Only when ⁇ rays are simultaneously incident on the two scintillator blocks 31, the sent image information is determined as appropriate data. When ⁇ rays are incident only on one scintillator block 31, the coincidence counting circuit 11 rejects.
- the image information sent to the coincidence circuit 11 is sent to the reconstruction unit 12 as projection data.
- the reconstruction unit 12 reconstructs the projection data to obtain an image of the subject M.
- the image is sent to the output unit 7 via the controller 5. In this manner, nuclear medicine diagnosis is performed based on the image obtained by the reconstruction unit 12.
- the ⁇ -ray detector 3 includes a scintillator block 31 composed of a plurality of scintillator elements, and a photomultiplier tube (hereinafter simply referred to as “PMT”) optically coupled to the scintillator block 31. ”Is abbreviated as“ ”.
- Each scintillator element constituting the scintillator block 31 converts ⁇ rays into light by emitting light with the incidence of ⁇ rays. By this conversion, the scintillator element detects ⁇ rays. The light emitted from the scintillator element is sufficiently diffused by the scintillator block 31 and input to the PMT 32.
- the PMT 32 multiplies the light converted by the scintillator block 31 and converts it into an electrical signal.
- the electric signal is sent as image information (pixel value) to the position calculation circuit 9 (see FIG. 1) and further to the coincidence counting circuit 11 (see FIG. 1) as described above.
- the scintillator elements constituting the scintillator block 31 correspond to the scintillator elements in the present invention
- the photomultiplier tube (PMT) 32 corresponds to the photosensor in the present invention.
- the size of each ⁇ -ray detector 3 is not limited to this.
- a ⁇ -ray detector 3 including a scintillator block 31 in which scintillator elements are arranged in 32 ⁇ 32 ⁇ 4 layers and a PMT 32 having a 16 ⁇ 16 multi-anode is used.
- the number of scintillator elements constituting the scintillator block 31 and the number of multi-anodes of the PMT 32 are not particularly limited.
- FIG. 4A is a plan view of a two-dimensional position map for explanation when obtaining a minimum value along the row (x) direction
- FIG. 4B is a minimum along the column (y) direction
- FIG. 5A is a plan view of a two-dimensional position map schematically showing peaks
- FIG. 5B is a plan view of a two-dimensional position map used for explanation when obtaining values.
- FIG. 6A is a plan view of a two-dimensional position map that represents the sensitivity ratio for each scintillator element that represents the spatial periodicity of FIG. 6, and
- FIG. 6B is a plan view schematically showing a two-dimensional position map, and FIG.
- FIG. 6B shows the sensitivity for each scintillator element representing the spatial periodicity of the peak when the target region is separated when the separation fails.
- FIG. 7 is a plan view of a two-dimensional position map showing a ratio.
- FIG. 8 is a plan view of a two-dimensional position map for explaining the determination of the boundary when separating the target area when defeating.
- FIG. 8 shows the boundary when the target area is separated when the separation fails. It is a top view of the two-dimensional position map which described the sensitivity ratio with which it uses for the description to determine.
- Input a 2D position map (see Fig. 9).
- a Na-22 radiation source is uniformly irradiated from above the ⁇ -ray detector 3 to obtain a two-dimensional position map, and then the region is divided. The entire map is labeled with scintillator element position numbers to create a table, which is written and stored in the lookup table 10.
- the two-dimensional position map is an image of 1024 ⁇ 1024 pixels (pixels)
- the look-up table (LUT) has the position numbers of the scintillator elements in an array of 1024 ⁇ 1024. That is, in the lookup table (LUT), each position (1024 ⁇ 1024 pixels) in the two-dimensional position map is associated with each scintillator element.
- the scintillator block 31 (that is, crystal block) constituting the ⁇ -ray detector 3 is manufactured by assembling a large number (32 ⁇ 32 ⁇ 4 in this embodiment) of scintillator elements (that is, small crystals).
- the number of pixel value peaks is determined by the total number of small crystals constituting the crystal block.
- the total number of peaks among the rows and columns constituting the two-dimensional position map, the number along the row direction is N X and the number along the column direction is N Y. Therefore, the total number of peaks is N X ⁇ N Y. Further, as shown in FIG.
- the density (pixel value) of the peak Peak changes spatially, which is derived from the fact that the sensitivity differs for each layer of the block. It depends on whether it is composed of small crystals.
- the period / number of layers is t.
- pixel values that are signal intensities are compared to determine respective local minimum values or local maximum values, and the positions of these local minimum values or local maximum values are drawn as boundaries, and each position is determined by the boundary. Isolate.
- FIG. 4 (b) along with defining drawing boundaries seeking a minimum value along each line L X each parallel to the line (x) direction as shown in FIG. 4 (a), shown in FIG. 4 (b) determined by drawing a boundary seeking a minimum value along each line L Y of parallel each to the column (y) direction as.
- a graph G X (that is, a pixel value profile in the row direction) is created with the line L X as the horizontal axis and the pixel value as the vertical axis. Then, it extracts a point which is a minimum value in the graph G X (shown in FIG. 4, " ⁇ ").
- Draw boundaries B X at positions corresponding to the minimum value with respect to the line L X of interest. Assume that M X peaks are extracted in the row (x) direction by separating each position by the boundary B X.
- FIG. 4A a graph G X (that is, a pixel value profile in the row direction) is created with the line L X as the horizontal axis and the pixel value as the vertical axis. Then, it extracts a point which is a minimum value in the graph G X (shown in FIG. 4, " ⁇ ").
- Draw boundaries B X at positions corresponding to the minimum value with respect to the line L X of interest. Assume that M X
- a graph G Y (that is, a pixel value profile in the column direction) is created with the line LY as the horizontal axis and the pixel value as the vertical axis. Then, it extracts a point which is a minimum value in the graph G X (shown in FIG. 4, " ⁇ ").
- Draw boundary B Y at a position corresponding to a minimum value with respect to the line L Y of interest.
- M Y-number of peaks in the column (y) direction separates the position of each through the boundary B Y is extracted. From the above, M X ⁇ M Y peaks are extracted (M X ⁇ N X , M Y ⁇ N Y ). This separation operation corresponds to the peak separation step in the present invention.
- the sensitivity ratio and the pixel value are in a proportional relationship, if the sensitivity ratio for each scintillator element is obtained in advance, the pattern of the sensitivity ratio and the shading (pixel value) of the image And the number of peaks can be determined.
- the sensitivity ratio changes spatially as shown in FIG. 5B, and the shading (pixel value) of the image also changes spatially as shown in FIG. 5A.
- t 4
- the sensitivity ratio of the first layer is r1
- the sensitivity ratio of the second layer is r2
- the sensitivity ratio of the third layer is r3
- the sensitivity ratio of the fourth layer is r4.
- the sensitivity ratio and the density of the separated image shown in FIG. 5A correspond to each other. Therefore, the number of peaks that have failed to be separated can be determined using the spatial periodicity of the peaks.
- the sensitivity ratio for each scintillator element may be obtained in advance by experiment, or r1, r2, r3, r4 obtained from the sensitivity ratio at the center of the two-dimensional position map M may be used.
- the left peak is re-separated when the separation fails will be described with reference to FIG.
- the target region is indicated by a thick frame in FIGS.
- two peaks Peak are also crushed due to the connection in the up and down direction. Assume that it has been determined.
- the sensitivity ratio of the target region T when the separation fails is as shown in FIG. 6B, and the sensitivity ratio of the target region T when the separation fails is (r2). , R1, r1, r2, r1, r1).
- the boundary (boundary B X in the case of FIGS. 6 to 8) at the target region T (the left end in FIGS. 6 to 8) when the separation fails is once erased. .
- the left side N L ′ and the right side N R ′ were used.
- boundaries are determined from empirically known sensitivity ratios, and the respective positions where the separation fails are separated, so that N L ′ at the ends are separated. , N R ′ peaks are separated.
- a thick line is inserted so as to divide the target region T into left and right when the separation fails.
- the thick line dividing the region T on the left and right is moved left and right in the direction of the arrow in FIG. 7 to re-separate the thick line where the sensitivity ratio and the sum ratio of the pixels in the peak region match. Determine as the boundary to be.
- the two-dimensional position map calibration unit 13 reads out from the lookup table 10 and rewrites and calibrates the two-dimensional position map to create a lookup table.
- Radiation detection apparatus (gantry 2, ⁇ -ray detector 3, position calculation circuit 9, lookup table 10, coincidence counting circuit 11, and two-dimensional position map calibration provided in the PET apparatus according to the present embodiment having the above-described configuration.
- the two-dimensional position map calibration unit 13 draws a boundary by the peak of the pixel value as the signal intensity, and separates each position by the boundary. If separation fails due to a plurality of peaks connected, the number of failed peaks is determined using the spatial periodicity of the peaks. Therefore, the number of failed peaks can be determined by using the peak spatial periodicity, and the boundary can be easily determined. As a result, the incident position can also be easily distinguished, and the detection position of radiation (gamma rays in the present embodiment) can be easily identified.
- each position where the separation has failed is separated by defining a boundary so that the sensitivity ratio for each scintillator element matches the sum ratio of the pixels in the peak region.
- the sensitivity ratio and the pixel are in a proportional relationship, if the sensitivity ratio for each scintillator element is determined in advance, the sensitivity ratio for each scintillator element is separated in order to separate each position where the separation failed. It is possible to define the boundary so that the sum ratio of the pixels in the peak region matches.
- pixel values that are signal intensities are compared to determine respective local minimum values or local maximum values, and the positions of those local minimum values or local maximum values are drawn as boundaries, and the respective positions are separated by the boundaries. is doing.
- the peak is considered as the maximum, the position of the minimum value almost coincides with the boundary, so that the position of the minimum value is drawn as a boundary as shown in FIG. 4, and each position is separated by the boundary.
- the present invention is not limited to the above embodiment, and can be modified as follows.
- the PET apparatus has been described as an example of the nuclear medicine diagnosis apparatus provided with the radiation detection apparatus.
- the present invention detects tomographic images of a subject by detecting a single gamma ray. It can also be applied to a SPECT (Single-Photon-Emission-CT) apparatus that reconfigures the image.
- SPECT Single-Photon-Emission-CT
- the present invention can also be applied to a PET-CT apparatus that combines a PET apparatus and a CT apparatus. Further, the present invention can also be applied to radiation other than ⁇ (for example, ⁇ rays and ⁇ rays).
- the DOI detector is composed of a plurality of scintillator elements arranged three-dimensionally, but is also applied to a radiation detector composed of a plurality of scintillator elements arranged two-dimensionally or three-dimensionally. can do.
- a photomultiplier tube (PMT) has been described as an example of the photosensor.
- any photosensor optically coupled to the scintillator element may be an avalanche photodiode or silicon photomultiplier. It is not particularly limited as exemplified by the prior.
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Optics & Photonics (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
Description
すなわち、この発明の2次元位置マップ校正方法は、1次元,2次元あるいは3次元に配置された複数のシンチレータ素子と、それらに光学的に結合された光センサとで構成された放射線検出器を検出する際に用いられる、前記光センサで得られた信号強度を前記シンチレータ素子に入射された前記放射線の入射位置に対応させて2次元状に表した2次元位置マップを校正する2次元位置マップ校正方法であって、前記信号強度のピークによって境界を描いて、その境界によって各々の位置を分離するピーク分離工程と、複数の前記ピークが連結することで前記ピーク分離工程での分離に失敗したピークの個数を、前記ピークの空間周期性を用いて判定する個数判定工程とを備えていることを特徴とするものである。
10 … ルックアップテーブル
13 … 2次元位置マップ校正部
31 … シンチレータブロック
32 … 光電子増倍管(PMT)
M … 2次元位置マップ
Claims (5)
- 1次元,2次元あるいは3次元に配置された複数のシンチレータ素子と、それらに光学的に結合された光センサとで構成された放射線検出器を検出する際に用いられる、前記光センサで得られた信号強度を前記シンチレータ素子に入射された前記放射線の入射位置に対応させて2次元状に表した2次元位置マップを校正する2次元位置マップ校正方法であって、前記信号強度のピークによって境界を描いて、その境界によって各々の位置を分離するピーク分離工程と、複数の前記ピークが連結することで前記ピーク分離工程での分離に失敗したピークの個数を、前記ピークの空間周期性を用いて判定する個数判定工程とを備えていることを特徴とする2次元位置マップ校正方法。
- 請求項1に記載の2次元位置マップ校正方法において、前記シンチレータ素子ごとの感度比と前記ピーク領域内の画素の総和比とが一致するように前記境界を定めることで、前記ピーク分離工程での分離に失敗した各々の位置を分離する境界決定工程を備えることを特徴とする2次元位置マップ校正方法。
- 請求項1または請求項2に記載の2次元位置マップ校正方法において、前記分離工程では、前記信号強度を比較して各々の極小値をそれぞれ求めて、それらの極小値の位置を前記境界として描いて、その境界によって各々の位置を分離することを特徴とする2次元位置マップ校正方法。
- 請求項1または請求項2に記載の2次元位置マップ校正方法において、前記分離工程では、前記信号強度を比較して各々の極大値をそれぞれ求めて、それらの極大値の位置を前記境界として描いて、その境界によって各々の位置を分離することを特徴とする2次元位置マップ校正方法。
- 1次元,2次元あるいは3次元に配置された複数のシンチレータ素子と、それらに光学的に結合された光センサとで構成された放射線検出器を備えた放射線検出装置であって、前記光センサで得られた信号強度を前記シンチレータ素子に入射された前記放射線の入射位置に対応させて2次元状に表した2次元位置マップについて、その2次元位置マップにおける各位置と各々のシンチレータ素子とを対応させたテーブルを記憶する記憶手段と、前記2次元位置マップを校正するための演算処理を行う演算処理手段とを備え、前記校正された2次元位置マップと放射線の検出結果とに基づいて前記入射位置を弁別することで放射線の検出位置を特定し、前記演算処理手段は、前記信号強度のピークによって境界を描いて、その境界によって各々の位置を分離するピーク分離工程と、複数の前記ピークが連結することで前記ピーク分離工程での分離に失敗したピークの個数を、前記ピークの空間周期性を用いて判定する個数判定工程とを備え、それらの工程に関する演算処理を行うことを特徴とする放射線検出装置。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP08764354.0A EP2287635B1 (en) | 2008-05-19 | 2008-05-19 | Two-dimensional position map correcting method, and radiation detecting apparatus |
US12/918,595 US8471211B2 (en) | 2008-05-19 | 2008-05-19 | Two-dimensional position map correcting method, and radiation detecting apparatus |
CN2008801281062A CN101978287B (zh) | 2008-05-19 | 2008-05-19 | 二维位置图校正方法及放射线检测装置 |
PCT/JP2008/059146 WO2009141861A1 (ja) | 2008-05-19 | 2008-05-19 | 2次元位置マップ校正方法および放射線検出装置 |
JP2010512869A JP4983981B2 (ja) | 2008-05-19 | 2008-05-19 | 2次元位置マップ校正方法および放射線検出装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2008/059146 WO2009141861A1 (ja) | 2008-05-19 | 2008-05-19 | 2次元位置マップ校正方法および放射線検出装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2009141861A1 true WO2009141861A1 (ja) | 2009-11-26 |
Family
ID=41339836
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2008/059146 WO2009141861A1 (ja) | 2008-05-19 | 2008-05-19 | 2次元位置マップ校正方法および放射線検出装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8471211B2 (ja) |
EP (1) | EP2287635B1 (ja) |
JP (1) | JP4983981B2 (ja) |
CN (1) | CN101978287B (ja) |
WO (1) | WO2009141861A1 (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011155103A1 (ja) * | 2010-06-10 | 2011-12-15 | 株式会社島津製作所 | 2次元位置マップ校正方法 |
JP2017537310A (ja) * | 2014-10-17 | 2017-12-14 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 光共有及び相互作用深さ推定を用いるpet検出器シンチレータアレンジメント |
JP2018105648A (ja) * | 2016-12-22 | 2018-07-05 | 浜松ホトニクス株式会社 | 放射線位置検出方法、放射線位置検出器及びpet装置 |
WO2019008645A1 (ja) * | 2017-07-03 | 2019-01-10 | 株式会社 島津製作所 | 放射線検出装置およびそれを備えた核医学診断装置 |
JP2020190446A (ja) * | 2019-05-21 | 2020-11-26 | 株式会社島津製作所 | 放射線検出装置の2次元位置マップの校正方法および放射線検出装置 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8809793B2 (en) * | 2012-01-27 | 2014-08-19 | General Electric Company | System and method for pixelated detector calibration |
EP3278142B1 (en) * | 2015-03-17 | 2018-11-07 | Koninklijke Philips N.V. | Scintillation event localization in a radiation particle detector |
JP6607256B2 (ja) * | 2015-06-18 | 2019-11-20 | 株式会社島津製作所 | 放射線検出素子の感度補正方法および放射線断層撮影装置 |
US10692188B2 (en) | 2016-12-03 | 2020-06-23 | Tose'e Sanaye Tasvirbardari Parto Negar Persia Company Ltd. | Peak detection in a two dimensional image |
CN107272043B (zh) * | 2017-06-05 | 2019-06-04 | 中派科技(深圳)有限责任公司 | 检测器和具有该检测器的发射成像设备 |
JP7109168B2 (ja) * | 2017-09-13 | 2022-07-29 | 浜松ホトニクス株式会社 | 放射線位置検出方法、放射線位置検出器及びpet装置 |
CN112690818B (zh) * | 2020-12-29 | 2022-06-17 | 赛诺联合医疗科技(北京)有限公司 | Pet探测器的晶体位置查找表的校正方法及pet系统 |
US20220252746A1 (en) * | 2021-02-05 | 2022-08-11 | Canon Medical Systems Corporation | Identifying arrangement errors of detector elements within a gamma-ray detector system |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1123721A (ja) * | 1997-07-07 | 1999-01-29 | Hamamatsu Photonics Kk | 放射線位置検出装置 |
JPH11142524A (ja) * | 1997-11-07 | 1999-05-28 | National Institute Of Radiological Sciences | 放射線入射位置3次元検出器の発光位置特定方法 |
JP2005043104A (ja) | 2003-07-23 | 2005-02-17 | Shimadzu Corp | 放射線位置検出器の校正方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000004403A1 (fr) * | 1998-07-15 | 2000-01-27 | Keiichi Kuroda | Unite d'imagerie a rayonnement numerique |
US6297506B1 (en) * | 2000-03-23 | 2001-10-02 | John W. Young | System and method for reducing pile-up errors in multi-crystal gamma ray detector applications |
US6710349B2 (en) * | 2002-04-16 | 2004-03-23 | General Electric Company | Edge resolved dual scintillator gamma ray detection system and method |
US7449680B2 (en) * | 2006-04-28 | 2008-11-11 | Siemens Medical Solutions Usa, Inc. | Non-uniformity energy correction method and apparatus |
US8063379B2 (en) * | 2006-06-21 | 2011-11-22 | Avraham Suhami | Radiation cameras |
US8278625B2 (en) * | 2008-03-21 | 2012-10-02 | Shimadzu Corporation | Two-dimensional position map correcting method |
-
2008
- 2008-05-19 US US12/918,595 patent/US8471211B2/en not_active Expired - Fee Related
- 2008-05-19 CN CN2008801281062A patent/CN101978287B/zh not_active Expired - Fee Related
- 2008-05-19 WO PCT/JP2008/059146 patent/WO2009141861A1/ja active Application Filing
- 2008-05-19 EP EP08764354.0A patent/EP2287635B1/en not_active Not-in-force
- 2008-05-19 JP JP2010512869A patent/JP4983981B2/ja active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1123721A (ja) * | 1997-07-07 | 1999-01-29 | Hamamatsu Photonics Kk | 放射線位置検出装置 |
JPH11142524A (ja) * | 1997-11-07 | 1999-05-28 | National Institute Of Radiological Sciences | 放射線入射位置3次元検出器の発光位置特定方法 |
JP2005043104A (ja) | 2003-07-23 | 2005-02-17 | Shimadzu Corp | 放射線位置検出器の校正方法 |
Non-Patent Citations (1)
Title |
---|
See also references of EP2287635A4 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011155103A1 (ja) * | 2010-06-10 | 2011-12-15 | 株式会社島津製作所 | 2次元位置マップ校正方法 |
CN102933983A (zh) * | 2010-06-10 | 2013-02-13 | 株式会社岛津制作所 | 二维位置图校正方法 |
JP5664880B2 (ja) * | 2010-06-10 | 2015-02-04 | 株式会社島津製作所 | 2次元位置マップ校正方法 |
EP2581763A4 (en) * | 2010-06-10 | 2015-12-23 | Shimadzu Corp | METHOD FOR CORRECTING A TWO DIMENSIONAL POSITIONS CARD |
JP2017537310A (ja) * | 2014-10-17 | 2017-12-14 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 光共有及び相互作用深さ推定を用いるpet検出器シンチレータアレンジメント |
JP2018105648A (ja) * | 2016-12-22 | 2018-07-05 | 浜松ホトニクス株式会社 | 放射線位置検出方法、放射線位置検出器及びpet装置 |
WO2019008645A1 (ja) * | 2017-07-03 | 2019-01-10 | 株式会社 島津製作所 | 放射線検出装置およびそれを備えた核医学診断装置 |
CN110770604A (zh) * | 2017-07-03 | 2020-02-07 | 株式会社岛津制作所 | 放射线检测装置以及具备该放射线检测装置的核医学诊断装置 |
JPWO2019008645A1 (ja) * | 2017-07-03 | 2020-03-26 | 株式会社島津製作所 | 放射線検出装置およびそれを備えた核医学診断装置 |
CN110770604B (zh) * | 2017-07-03 | 2023-10-24 | 株式会社岛津制作所 | 放射线检测装置以及具备该放射线检测装置的核医学诊断装置 |
JP2020190446A (ja) * | 2019-05-21 | 2020-11-26 | 株式会社島津製作所 | 放射線検出装置の2次元位置マップの校正方法および放射線検出装置 |
JP7247745B2 (ja) | 2019-05-21 | 2023-03-29 | 株式会社島津製作所 | 放射線検出装置の2次元位置マップの校正方法および放射線検出装置 |
Also Published As
Publication number | Publication date |
---|---|
EP2287635A1 (en) | 2011-02-23 |
US20100327168A1 (en) | 2010-12-30 |
US8471211B2 (en) | 2013-06-25 |
JPWO2009141861A1 (ja) | 2011-09-22 |
EP2287635A4 (en) | 2014-03-12 |
CN101978287A (zh) | 2011-02-16 |
JP4983981B2 (ja) | 2012-07-25 |
EP2287635B1 (en) | 2017-09-20 |
CN101978287B (zh) | 2013-05-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4983981B2 (ja) | 2次元位置マップ校正方法および放射線検出装置 | |
JP4835789B2 (ja) | 2次元位置マップ校正方法 | |
US9709684B2 (en) | Systems and methods for scintillators having micro-crack surfaces | |
WO2009101677A1 (ja) | 放射線検出器、およびそれを備えた断層撮影装置 | |
CN110520760A (zh) | 用于辐射检测器的像素设计 | |
JP2011503535A (ja) | 間接放射線検出器 | |
JP5459397B2 (ja) | ポジトロンct装置およびタイミング補正方法 | |
JP5664880B2 (ja) | 2次元位置マップ校正方法 | |
CN102665564A (zh) | 图像诊断装置及方法 | |
JP4993144B2 (ja) | Pet検出器の光量出力を補正する方法 | |
JP4367903B2 (ja) | 放射線位置検出器の校正方法及び放射線位置検出器 | |
JP4933767B2 (ja) | 放射線同時計数処理方法、放射線同時計数処理プログラムおよび放射線同時計数処理記憶媒体、並びに放射線同時計数装置およびそれを用いた核医学診断装置 | |
JP6694213B2 (ja) | 積層型放射線3次元位置検出器 | |
EP3278142B1 (en) | Scintillation event localization in a radiation particle detector | |
JP6827316B2 (ja) | 放射線位置検出方法、放射線位置検出器及びpet装置 | |
JP5082909B2 (ja) | フォトン計数方法、その装置並びにそれを用いた核医学診断装置 | |
JP2017003392A (ja) | 位置マップ作成方法および核医学診断装置 | |
JP7247745B2 (ja) | 放射線検出装置の2次元位置マップの校正方法および放射線検出装置 | |
JP4737104B2 (ja) | 光子検出器の位置情報算出方法及びそれを用いたポジトロンct装置 | |
KR20230095204A (ko) | 섬광 검출기 및 이를 이용한 양전자방출단층촬영장치 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 200880128106.2 Country of ref document: CN |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 08764354 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2010512869 Country of ref document: JP Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 12918595 Country of ref document: US |
|
REEP | Request for entry into the european phase |
Ref document number: 2008764354 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2008764354 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |