WO2011145139A1 - ポジトロンｃｔ装置およびタイミング補正方法 - Google Patents

ポジトロンｃｔ装置およびタイミング補正方法 Download PDF

Info

Publication number
WO2011145139A1
WO2011145139A1 PCT/JP2010/003340 JP2010003340W WO2011145139A1 WO 2011145139 A1 WO2011145139 A1 WO 2011145139A1 JP 2010003340 W JP2010003340 W JP 2010003340W WO 2011145139 A1 WO2011145139 A1 WO 2011145139A1
Authority
WO
WIPO (PCT)
Prior art keywords
detector
time difference
scintillator element
detectors
radiation
Prior art date
Application number
PCT/JP2010/003340
Other languages
English (en)
French (fr)
Inventor
津田倫明
佐藤允信
Original Assignee
株式会社島津製作所
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 株式会社島津製作所 filed Critical 株式会社島津製作所
Priority to EP10851704.6A priority Critical patent/EP2573588B1/en
Priority to PCT/JP2010/003340 priority patent/WO2011145139A1/ja
Priority to JP2012515640A priority patent/JP5459397B2/ja
Priority to US13/698,971 priority patent/US9360569B2/en
Priority to CN201080066866.2A priority patent/CN102906595B/zh
Priority to TW100117139A priority patent/TWI444646B/zh
Publication of WO2011145139A1 publication Critical patent/WO2011145139A1/ja
Priority to US15/142,904 priority patent/US9844351B2/en

Links

Images

Classifications

    • 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
    • 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4266—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a plurality of detector units
    • 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/48—Diagnostic techniques
    • A61B6/485—Diagnostic techniques involving fluorescence X-ray imaging
    • 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
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/17—Circuit arrangements not adapted to a particular type of detector
    • G01T1/172—Circuit arrangements not adapted to a particular type of detector with coincidence circuit arrangements
    • 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/20—Measuring radiation intensity with scintillation detectors
    • G01T1/208—Circuits specially adapted for scintillation detectors, e.g. for the photo-multiplier section
    • 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/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
    • G01T1/2914—Measurement of spatial distribution of radiation
    • G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)

Definitions

  • the present invention relates to a positron CT apparatus and timing correction method for detecting radiation emitted from a positron radiopharmaceutical administered into a subject.
  • a positron CT device ie, a PET (Positron EmissionographyTomography) device, detects positrons (Positron), that is, a plurality of ⁇ -rays generated by annihilation of positrons, and a plurality of detectors simultaneously detect ⁇ -rays ( In other words, only when simultaneous counting is performed, the image of the subject is reconstructed.
  • positron positron
  • TOF Time Of Flight
  • TOF uses the fact that the annihilation radiation is at the speed of light, and converts the time difference from the annihilation occurrence point to the detector into the distance difference from the annihilation occurrence point to the light source generation position by the scintillator element of the detector. This is a technology for finding the point of occurrence of annihilation.
  • the detector detects the radiation emitted from the radiation source, and inputs a timing signal indicating the radiation incident timing at which the radiation is incident on the detector to the coincidence counting circuit via the delay adjustment circuit.
  • the output of the coincidence counting circuit is measured, and the sensitivity (ie, counting) of radiation is measured for each signal channel. Thereafter, the sensitivity is measured while changing the delay amount adjusted by the delay adjustment circuit, and a sensitivity distribution with respect to the delay amount change is obtained. By using the delay amount with the highest measured sensitivity in the delay adjustment circuit, the time delay of the signal is corrected.
  • a calibration radiation source (external radiation source) is installed in the field of view (FOV) of the PET device.
  • a plurality of detectors are arranged in a ring (annular).
  • the timing values of a plurality of detectors having a common field of view with the reference detector are averaged, and the averaged timing value is obtained as a time delay value for the reference detector.
  • a time delay value is similarly obtained using a detector adjacent to the reference detector as a reference, and a difference between the time delay value obtained first and the time delay value obtained next is obtained as a reference correction value.
  • Timing correction is performed by aligning the time using the reference correction value. Thereafter, the same calculation is performed on the adjacent detectors in order, so that timing correction for all detectors is performed when the detector makes one turn on the ring.
  • the simulation signal output from the simulation signal generator is input to each of a plurality of signal processors (signal processing units), and calibration data is generated based on the output of each signal processor to correct timing.
  • a DOI detector capable of discriminating a depth-of-interaction light source position (DOI: Depth of Interaction) is incorporated into a TOF type PET apparatus.
  • the DOI detector is constructed by laminating each scintillator element in the depth direction of radiation (here, ⁇ -rays), and the depth direction and the lateral direction (direction parallel to the incident surface) that caused the interaction. ) And the coordinate information with the center of gravity.
  • the detection time correction information corresponding to the coordinate information is written and stored in the table, and the detection time correction information is referred to thereby improving the information accuracy of the flight time difference.
  • Japanese Patent Publication No. 6-19436 Japanese Patent No. 3343122 JP 2006-90827 A JP 2008-51701 A
  • the present invention has been made in view of such circumstances, and it is an object of the present invention to provide a positron CT apparatus and a timing correction method capable of accurately simultaneously counting without repeating many measurements and calculations. To do.
  • the present invention has the following configuration. That is, the positron CT apparatus of the present invention is a positron CT apparatus having a plurality of detectors for detecting radiation emitted from a positron radiopharmaceutical administered into a subject, and each detector for simultaneously counting radiation. With respect to a time difference histogram representing a count value distribution for each pair of time difference changes, two detectors to be simultaneously counted are selected, and one of the two selected detectors is selected, and When a detector different from the other detector is selected and the selection is repeated, the time difference histogram relating to the two detectors selected in the past is used as a reference, and the detector selected this time is selected.
  • the calculation means performs the following calculation on the time difference histogram representing the count value distribution with respect to the time difference change of each pair of detectors that simultaneously count radiation. That is, two detectors to be simultaneously counted are selected, one of the two selected detectors is selected, and another detector is selected from the other detector, When the selection is repeated, the time difference histograms related to the two detectors selected this time are corrected based on the time difference histograms related to the two detectors selected in the past. Then, the operation based on the time difference histogram relating to the two corrected detectors is newly performed.
  • the time difference histograms related to the two detectors selected this time are corrected based on the reference, and the two corrected detectors are corrected.
  • the time difference histogram as a reference for example, as in Patent Document 2 described above, convergence can be obtained when obtaining an optimum time difference histogram, compared to the case where the timing values of a plurality of detectors are averaged. Good. Therefore, an optimal time difference histogram can be obtained without repeating many measurements and calculations.
  • the coincidence counting circuit simultaneously counts the radiation based on the time difference histogram of each pair of the detectors corrected repeatedly as described above, the coincidence can be accurately performed. As a result, simultaneous counting can be performed accurately without repeating many measurements and calculations.
  • the timing correction method of the present invention is a timing correction method used for simultaneously counting radiation emitted from a positron radiopharmaceutical administered into a subject, and is a pair of detectors that simultaneously count radiation.
  • a time difference histogram representing a count value distribution for each time difference change
  • two detectors to be simultaneously counted are selected, one of the two selected detectors is selected, and the other detector is selected.
  • the time difference histograms related to the two detectors selected in the past are used as a reference, and the two selected this time based on the reference. Correct the time difference histogram related to the detector, and repeat the work using the corrected time difference histogram related to the two detectors as a new reference. It is characterized in further comprising a cormorants histogram correction process.
  • the timing correction method of the present invention the following correction is performed in the histogram correction step with respect to the time difference histogram for each detector pair. That is, two detectors to be simultaneously counted are selected, one of the two selected detectors is selected, and another detector is selected from the other detector, When the selection is repeated, the time difference histograms related to the two detectors selected this time are corrected based on the time difference histograms related to the two detectors selected in the past. Then, the operation based on the time difference histogram relating to the two corrected detectors is newly performed. As a result, timing correction can be performed accurately without repeating many measurements and calculations.
  • One example of these positron CT devices and timing correction methods described above uses the time difference at which the total count value in the time difference histogram based on the above reference is the maximum as a reference value, and corrects the time difference of the time difference histogram based on the above reference value. It is.
  • the time point at which the total count value in the time difference histogram is maximum is the timing at which coincidence counting is most likely to occur. Therefore, by correcting the time difference of the time difference histogram based on the reference value that is the timing, it is possible to align at that timing.
  • Another example of the above-described positron CT apparatus and timing correction method is to calculate a time difference that is an intermediate value between the time difference at which the total count value is maximum and the time difference at which the total count value is the second largest in the reference time difference histogram.
  • the reference value is used to correct the time difference of the time difference histogram based on the above-described reference value.
  • the timing at which the scintillator element of the other detector detects the most radiation from the self-radiation of the scintillator element of one detector is the time point where the total count value is maximum, or the total count value is the second largest It is a time point. Therefore, the intermediate value between the two time differences, which are these timings, is the timing at which coincidence counting is most likely to occur. Therefore, by correcting the time difference of the time difference histogram based on the reference value of the intermediate value that is the timing, it is possible to align at that timing.
  • the above-described detector has a structure including a plurality of scintillator elements in addition to a structure including a single scintillator element.
  • correction may be performed for each scintillator element group (that is, detector unit) composed of a plurality of scintillator elements, but correction is performed for each scintillator element unit composed of one scintillator element as described below. To improve accuracy.
  • the following calculation / correction is performed on the time difference histogram for each pair of scintillator elements composed of one scintillator element of each detector that simultaneously counts radiation. That is, two scintillator element units of a detector to be simultaneously counted are selected, and one scintillator element unit of one detector is selected from the two selected scintillator element units.
  • the time difference histogram for the scintillator element units of the two detectors selected in the past is used as a reference, and this time, The time difference histogram relating to the scintillator element unit of the two detectors selected in (1) is corrected.
  • a time difference histogram for each pair of a scintillator element group composed of a plurality of scintillator elements of one detector and a scintillator element unit composed of one scintillator element of the other detector is performed.
  • the scintillator element group of one detector and the scintillator element unit of the other detector selects the scintillator element group of one detector and the scintillator element unit of the other detector, and out of the selected scintillator element group and the scintillator element unit,
  • the scintillator element group and the scintillator element unit selected in the past The time difference histogram relating to the scintillator element group selected this time and the scintillator element unit is corrected based on the reference.
  • the operation based on the time difference histogram relating to the corrected scintillator element group and scintillator element unit is newly performed. Since correction is performed for each scintillator element group and each scintillator element in this way, the accuracy can be further improved compared to when correction is performed for each detector unit. In addition, the calculation time and burden can be reduced compared to when correction is performed for each scintillator element unit.
  • the above-described detector may be a DOI detector configured by laminating each scintillator element in the depth direction of the radiation.
  • positron CT apparatuses and timing correction methods described above include an external radiation source that irradiates the same type of radiation as that of the above-described radiopharmaceutical or a phantom that irradiates the same type of radiation as that of the radiopharmaceutical from the inside.
  • the detector may be acquired based on radiation from a radiation source or a phantom, or the detector may include a scintillator element having self-radiation, and the above-described time difference histogram may be acquired based on radiation from the self-radioactivity. .
  • the positron CT apparatus and the timing correction method according to the present invention, two detectors to be simultaneously counted are selected, one of the two selected detectors is selected, and the other is selected.
  • the time difference histograms related to the two detectors selected in the past are used as a reference, and the 2 selected this time is selected based on the reference. Correct the time difference histogram for two detectors. Then, the operation based on the time difference histogram relating to the two corrected detectors is newly performed. As a result, timing correction can be performed accurately without repeating many measurements and calculations.
  • PET Positron
  • A is a schematic perspective view of a gamma ray detector. It is the front view of the gamma ray detector arranged in the shape of a ring in a PET device, and a block diagram about it. It is a flowchart which shows the flow of a series of timing correction methods.
  • A is a front view of a ⁇ -ray detector when an external radiation source is installed
  • (b) is a front view of the ⁇ -ray detector when a phantom is installed.
  • A), (b) is explanatory drawing of a time difference histogram.
  • FIGS. 8A to 8C are front views showing an embodiment of switching between a reference detector and a correction target detector in another embodiment different from FIG. (A) to (c) are front views showing an embodiment of switching between a reference scintillator element unit and a correction target scintillator element unit.
  • (A) to (c) are front views showing an embodiment of switching between a reference scintillator element group and a correction target scintillator element unit.
  • (A), (b) is explanatory drawing of the time difference histogram at the time of the scintillator element with a self-radioactivity.
  • FIG. 1 is a side view and block diagram of a PET (Positron Emission Tomography) apparatus according to a first embodiment
  • FIG. 2 is a schematic perspective view of a ⁇ -ray detector
  • FIG. 3 is a ring shape of the PET apparatus. It is the front view of the arrange
  • the PET apparatus includes a top plate 1 on which the 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 detector in the present invention.
  • the PET apparatus includes a top board driving unit 4, a controller 5, an input unit 6, an output unit 7, a memory unit 8, a detector signal processing unit 9, a coincidence counting circuit 10, a data collecting / A control unit 11 and a reconstruction processing unit 12 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 coincidence circuit 10 corresponds to the coincidence circuit in the present invention
  • the data collection / control unit 11 corresponds to the computing means in the present invention.
  • the controller 5 comprehensively controls each part constituting the PET apparatus according to the first embodiment.
  • the controller 5 and the data collection / control unit 11 include 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 includes a storage medium represented by ROM (Read-only Memory), RAM (Random-Access Memory), and the like.
  • the count value (count) simultaneously counted by the coincidence circuit 10 data relating to coincidence counting such as a detector pair including two ⁇ -ray detectors 3 and LOR that have been simultaneously counted, and the reconstruction processing unit 12.
  • the image processed in step 1 is written and stored in the RAM, and is read from the RAM as necessary.
  • the ROM stores in advance a program for performing imaging including various types of nuclear medicine diagnosis, and the controller 5 and the data collection / control unit 11 execute the program so that the nuclear medicine corresponding to the program is stored. Make each diagnosis.
  • LOR Line Of Response
  • LOR Line Of Response
  • the reconfiguration processing unit 12 causes the controller 5 to execute, for example, a program stored in a ROM of a storage medium represented by the memory unit 8 or the like described above, or a command input by a pointing device represented by the input unit 6 or the like. It is realized with.
  • the ⁇ -rays generated from the subject M to which the radiopharmaceutical is administered are converted into light by the scintillator block 31 (see FIG. 2) of the ⁇ -ray detector 3, and the converted light is photoelectron of the ⁇ -ray detector 3.
  • a multiplier tube (PMT: Photo Multiplier Tube) 33 (see FIG. 2) multiplies and converts it into an electrical signal. The electric signal is sent to the coincidence counting circuit 10 through the detector signal processing unit 9.
  • the coincidence circuit 10 checks the position of the scintillator block 31 (see FIG. 2) and the incident timing of ⁇ rays, and only when ⁇ rays are simultaneously incident on the two scintillator blocks 31 on both sides of the subject M. The sent event is determined to be appropriate data.
  • the coincidence counting circuit 10 rejects. That is, the coincidence counting circuit 10 detects that ⁇ rays are simultaneously observed in the two ⁇ ray detectors 3 based on the above-described electrical signal.
  • an electrical signal that has been simultaneously observed is determined as image information, and the image information is sent to the reconstruction processing unit 12 via the data collection / control unit 11.
  • the reconstruction processing unit 12 obtains an image of the subject M by performing image reconstruction by forward projection processing or back projection processing.
  • the image is sent to the output unit 7 via the controller 5.
  • nuclear medicine diagnosis is performed based on the image obtained by the reconstruction processing unit 12.
  • a known successive approximation algorithm such as DRAMA method (Dynamic Row-Action Maximum Likelihood Algorithm) is applied. Specific functions of the detector signal processing unit 9 and the data collection / control unit 11 will be described later.
  • the ⁇ -ray detector 3 includes a scintillator block 31, a light guide 32 optically coupled to the scintillator block 31, and photoelectrons optically coupled to the light guide 32.
  • a multiplier (hereinafter simply abbreviated as “PMT”) 33 is provided.
  • 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. Light emitted from the scintillator element is sufficiently diffused by the scintillator block 31 and input to the PMT 33 via the light guide 32.
  • the PMT 33 multiplies the light converted by the scintillator block 31 and converts it into an electric signal.
  • the electric signal is sent to the coincidence counting circuit 10 through the detector signal processing unit 9 as described above.
  • the scintillator block 31 corresponds to the scintillator element in the present invention
  • the PMT 33 corresponds to the photoelectric conversion means in the present invention.
  • the ⁇ -ray detector 3 shown in FIG. 2 includes a plurality of scintillator blocks 31 that fluoresce when ⁇ -rays are incident, and a PMT 33 that detects ⁇ -rays by photoelectrically converting light from each scintillator block 31. It has.
  • the ⁇ -ray detector 3 shown in FIG. 2 is a DOI detector configured by laminating the respective scintillator blocks 31 in the ⁇ -ray depth direction (stacked in four layers in FIG. 2).
  • the ⁇ -ray detector 3 In imaging in a nuclear medicine diagnosis using a normal subject M, the ⁇ -ray detector 3 detects ⁇ -rays generated from the subject M to which a radiopharmaceutical has been administered, and the detector signal processing unit 9 performs simultaneous counting. The data is sent to the reconstruction processing unit 12 via the circuit 10 via the data collection / control unit 11 without timing correction. On the other hand, in data collection using an external source or phantom, the ⁇ -ray detector 3 detects ⁇ -rays from the external source or phantom, and the data is collected through the detector signal processing unit 9 and the coincidence counting circuit 10. Data is sent to the control unit 11, and the data collection / control unit 11 performs feedback correction to the timing correction table 9 b of the detector signal processing unit 9 to perform timing correction, thereby adjusting the delay amount in each coincidence circuit 10. .
  • the detector signal processing unit 9 generates a signal (hereinafter also referred to as “time stamp”) representing the ⁇ -ray incident timing based on the electrical signal output from the PMT 33 of the ⁇ -ray detector 3. And the above-described timing correction table 9b.
  • time stamp a signal representing the ⁇ -ray incident timing based on the electrical signal output from the PMT 33 of the ⁇ -ray detector 3.
  • the above-described timing correction table 9b As described above, each ⁇ -ray detector 3 is arranged in a ring shape, and each ⁇ -ray detector 3 is connected to the timing signal generation circuit 9a of the detector signal processing unit 9 (in FIG. 3, Only two timing signal generation circuits 9a are shown).
  • the timing signal generation circuit 9 a is connected to the timing correction table 9 b, and the timing correction table 9 b is connected to the coincidence counting circuit 10.
  • FIG. 4 is a flowchart showing a flow of a series of timing correction methods
  • FIG. 5 is a front view of a ⁇ -ray detector when an external radiation source or a phantom is installed
  • FIG. 6 is an explanatory diagram of a time difference histogram.
  • FIG. 7 is a front view showing an embodiment of switching between the reference detector and the correction target detector
  • FIG. 8 is an embodiment different from FIG. 7 in the reference detector and the correction target detection. It is a front view which shows one embodiment of switching with a vessel.
  • Step S1 Measurement of Correction Data
  • an external radiation source RI for irradiating a radioactive medicine that is, a radiation of the same kind as a radioisotope (RI) (gamma rays in the first embodiment) is used.
  • the two ⁇ -ray detectors 3 facing each other by about 180 ° can detect the radiation from the external radiation source RI almost simultaneously. Only the time delay of the signal in the signal channel from the line detector 3 to the coincidence circuit 10 is subject to timing correction.
  • a phantom Ph that irradiates the same type of radiation as the radioactive drug from the inside may be installed in the visual field of the PET apparatus. Even when the phantom Ph is installed, it is preferably installed in the central region within the field of view.
  • the external radiation source RI corresponds to the external radiation source in the present invention
  • the phantom Ph corresponds to the phantom in the present invention.
  • the radiation from the external radiation source RI or phantom Ph is detected by the ⁇ -ray detector 3 and acquired as correction data, and sent to the coincidence circuit 10 via the detector signal processing unit 9. As shown in FIG. 5, correction data is measured and acquired by all the ⁇ -ray detectors 3.
  • Step S2 Creation of Time Difference Histogram
  • a time difference histogram is created for each pair of ⁇ -ray detectors 3 that simultaneously count radiation based on the correction data measured in step S1.
  • the time difference histogram is obtained by taking a time stamp difference (that is, time difference) (indicated as “Difference Time” in FIG. 6) for each pair of ⁇ -ray detectors 3 on the horizontal axis, and counting values (FIG. 6).
  • time difference Time that is, time difference
  • FIG. 6 it is a count value distribution with respect to the time difference change when the vertical axis is “Event Counts”.
  • the time difference histogram is obtained as a frequency distribution centered on “0” as shown in FIG.
  • the time difference histogram shifts to the left and right as shown by the solid line in FIG. Therefore, a shift amount that returns the time difference histogram to a frequency distribution centered on “0” shown by the dotted line from the solid line in FIG. 6B is considered as the correction amount.
  • Step S3 Setting of Reference Detector and Correction Target Detector
  • a reference detector and a correction target detector are set.
  • one ⁇ -ray detector 3 is set as a reference
  • the ⁇ -ray detector 3 set as the reference is a reference detector.
  • Each pair of ⁇ -ray detectors 3 in the field of view of the PET apparatus viewed from the reference detector S is set as a correction target
  • the ⁇ -ray detector 3 set as the correction target is set as a correction target detector C. .
  • the ⁇ -ray detector 3 facing the reference detector S by about 180 ° is selected. Assuming that the reference detector S is the “1” -th ⁇ -ray detector 3 and serial numbers are assigned in the clockwise direction, as shown in FIG. This ⁇ -ray detector 3 becomes the ⁇ -ray detector 3 facing the reference detector S by 180 °. In the case of a total of an odd number of ⁇ -ray detectors, the ⁇ -ray detectors 3 facing each other by 180 ° cannot be selected. Therefore, the ⁇ -ray detectors 3 facing each other by about 180 ° may be selected.
  • the ⁇ -ray detector 3 (refer to the gray color in FIG. 7A) facing the reference detector S by about 180 ° is set as the counter detector O.
  • the reference detector S that is the “1” -th ⁇ -ray detector 3 and 180 ° to it.
  • the opposite “61” -th ⁇ -ray detector 3 is selected.
  • the time difference histogram relating to the selected reference detector S, which is the “1” -th ⁇ -ray detector 3, and the “61” -th ⁇ -ray detector 3 (opposite detector O) opposed to the reference detector S by 180 ° is used as a reference.
  • the time difference histogram based on this reference is a frequency distribution centered on “0” as shown in FIG. Accordingly, the time difference at which the total count value in the reference time difference histogram is maximum is set to “0”, and this time difference “0” is set as the reference value.
  • Step S4 Extraction of Time Difference Histogram to be Corrected Among the two ⁇ -ray detectors 3 composed of the reference detector S and the counter detector O selected in step S3, the reference detector is one ⁇ -ray detector 3. While selecting S, a correction target detector C different from the counter detector O which is the other ⁇ -ray detector 3 is selected. That is, the correction target detector C and the reference detector S other than the counter detector O are selected, and the time difference histogram to be corrected is extracted. In the case shown in FIG. 7A, the correction target detector C including the “61” -th ⁇ -ray detector 3 (opposite detector O) in the step S3 has the “49” -th to “73” -th total.
  • time difference histograms relating to the “th” to “60” th and “62” to “73” th correction target detectors C are extracted as the time difference histograms relating to the two ⁇ -ray detectors 3 selected this time.
  • Step S5 Calculation of Correction Amount
  • the time difference histogram relating to the correction target detector C other than the counter detector O and the reference detector S is deviated from the reference time difference histogram.
  • FIG. Shift to the left and right as shown by the solid line. Therefore, the shift amount is corrected so that the time difference histogram relating to the two ⁇ -ray detectors 3 selected at this time in step S4 is returned to the frequency distribution centered on “0” shown by the dotted line from the solid line in FIG.
  • the data collection / control unit 11 obtains the quantity.
  • Step S6 Completion of calculation of correction amounts for all detectors? It is determined whether or not the calculation of the correction amount of all the ⁇ -ray detectors 3 has been completed. If not completed, the process proceeds to step S7. If completed, the process proceeds to step S8.
  • Step S7 Application of correction amount Applying the correction amount obtained in step S5, the time difference histogram relating to the two ⁇ -ray detectors 3 selected this time is a frequency centered on “0” by the correction amount.
  • the time difference histogram is corrected by shifting to return to the distribution. That is, the time difference histogram is corrected based on the time difference “0” that is the reference value.
  • the reference detector S which is the “1” -th ⁇ -ray detector 3, and the “49” th to “60” th and “62” th to “73” th respectively.
  • step S7 If the correction amount is applied and corrected in step S7, the process returns to step S3, and the work based on the two ⁇ -ray detectors 3 corrected in step S7 is newly performed.
  • the correction target detector among the corrected “49” th to “60” th and “62” th to “73” th correction target detectors C.
  • the “49” -th and “73” -th ⁇ -ray detectors 3 at both ends of C are newly used as a reference.
  • step S 3 as shown in FIG. 7B, the “49” -th ⁇ -ray detector 3 and the “73” -th ⁇ -ray detector 3 (see the black painting in FIG. 7B).
  • the reference detector S is assumed.
  • each pair of ⁇ -ray detectors in the field of view of the PET apparatus viewed from each reference detector S is set as a correction target, and the ⁇ -ray detection set as the correction target is set.
  • the device 3 is a correction target detector C.
  • the correction target detector C is a total of 24 ⁇ -ray detectors 3 from the “97” th to “120” th with respect to the reference detector S which is the “49th” ⁇ -ray detector 3.
  • the reference detector S which is the “73” th ⁇ -ray detector 3 there are a total of 24 ⁇ -ray detectors 3 from the “2” th to “25” th.
  • the reference detector S which is the “49” -th ⁇ -ray detector 3, and the “109” -th ⁇ -ray detector 3 (opposite detector O) opposed to the reference detector S by 180 ° (see the gray color in FIG. 7B)
  • the reference detector S which is the “73” th ⁇ -ray detector 3 and the “13th” ⁇ -ray detector 3 (opposite detector O) which is 180 ° opposite to the reference detector S (FIG. 7B ) (See gray).
  • step S4 the reference detector S that is the “49” -th ⁇ -ray detector 3 and “97” other than the “109” -th ⁇ -ray detector 3 (opposite detector O).
  • the time difference histograms related to the “C” to “108” and “110” to “120” correction target detectors C are extracted as the time difference histograms related to the two ⁇ -ray detectors 3 selected this time.
  • the time difference histograms related to the “14” th to “25” th correction target detectors C are extracted as the time difference histograms related to the two ⁇ -ray detectors 3 selected this time.
  • step S5 the correction amount is obtained in step S5
  • step S6 the correction amount calculation completion is determined in step S6, and if not completed, the process proceeds to step S7.
  • step S7 the correction amount is applied to perform correction, and the process returns to step S3, and the work based on the two ⁇ -ray detectors 3 corrected in step S7 is newly performed.
  • a total of four ends of the correction target detector C see the hatched lines in FIG. 7B.
  • the “25” -th and “97” -th ⁇ -ray detectors 3 on the lower side are newly set as a reference.
  • step S3 the “25” -th ⁇ -ray detector 3 and the “97” -th ⁇ -ray detector 3 (see black in FIG. 7C).
  • the reference detector S is assumed.
  • Each ⁇ -ray detector 3 pair within the field of view of the PET apparatus viewed from each reference detector S is set as a correction target, and the ⁇ -ray detector 3 set as the correction target is set as a correction target detector C.
  • the correction target detector C is a total of 23 ⁇ -ray detectors 3 from “74” to “96” with respect to the reference detector S which is the “25” -th ⁇ -ray detector 3.
  • the reference detector S which is the “97” -th ⁇ -ray detector 3 there are a total of 23 ⁇ -ray detectors 3 from the “26” th to “48” th.
  • the reference detector S which is the “25” -th ⁇ -ray detector 3, and the “85” -th ⁇ -ray detector 3 (opposite detector O) opposed to the reference detector S by 180 ° (see the gray color in FIG. 7C)
  • the reference detector S which is the “97” -th ⁇ -ray detector 3, and the “37” -th ⁇ -ray detector 3 (opposite detector O) opposite to the reference detector S by 180 ° (FIG. 7C ) (See gray).
  • step S4 the reference detector S which is the “25” -th ⁇ -ray detector 3 and the “74” -th to “84” -th and “-” other than the “85” -th ⁇ -ray detector 3 (opposite detector O).
  • the time difference histogram relating to each of the 86th to 96th correction target detectors C is extracted as the time difference histogram relating to the two ⁇ -ray detectors 3 selected this time.
  • step S4 the reference detector S which is the “97” -th ⁇ -ray detector 3 and the “26” -th to “36” -th other than the “37” -th ⁇ -ray detector 3 (opposite detector O).
  • the time difference histograms relating to each of the “38” th to “48” th correction target detectors C are extracted as time difference histograms relating to the two ⁇ -ray detectors 3 selected this time.
  • step S5 the correction amount is obtained.
  • step S6 it is determined whether the correction amount has been calculated. If not, the process proceeds to step S7.
  • step S7 the correction amount is applied to perform correction, and the process returns to step S3.
  • the third iteration of steps S6, S7, and S3 to S5 completes the calculation of the correction amounts of all the ⁇ -ray detectors 3 from the “1” to the “120”, and is completed in step S6. to decide. Then, the process proceeds to step S8.
  • Step S8 Setting of Timing Correction Table
  • the correction amount obtained in these steps S3 to S6 (three iteration loops in the case of the embodiment of FIG. 7) is detected from the data collection / control unit 11 by detector signal processing.
  • the timing correction table is set by writing in the timing correction table 9b of the unit 9.
  • the total number of ⁇ -ray detectors 120 is 120, but the number of ⁇ -ray detectors is not particularly limited.
  • the correction target detector C is “1” th.
  • the “36” th and “54th” ⁇ -ray detectors 3 at both ends of the correction target detector C are newly used as a reference.
  • the correction target detector C is “71” to “88” with respect to the reference detector S which is the “36” -th ⁇ -ray detector 3 (the counter detector O is “80”). ”) And a total of 18 ⁇ -ray detectors 3 with respect to the reference detector S that is the“ 54 ”-th ⁇ -ray detector 3, the“ 2 ”-th to“ 19 ”-th (opposite detector) O is the “10th”).
  • the “19” th and “71” th ⁇ -ray detectors 3 on the lower side of the total four ends of the correction target detector C are newly set as a reference.
  • the correction target detector C is “55” th to “70” -th with respect to the reference detector S which is the “19” -th ⁇ -ray detector 3 (the counter detector O is “63”). ”Th) 16 ⁇ -ray detectors 3 in total, and“ 20 ”th to“ 35 ”th (opposite detector) with respect to the reference detector S which is the“ 71 ”th ⁇ -ray detector 3 O is the “27th”).
  • the number of loops in steps S6, S7, and S3 to S5 is three, but the number of loops is not particularly limited.
  • the outline of the setting of the reference detector S and the correction target detector O in the embodiment of FIG. 8 will be described.
  • the first loop as shown in FIG. 3 and the “61” -th ⁇ -ray detector 3 (opposite detector O) facing each other by 180 ° are selected, and the “60” -th and “62” on both sides adjacent to the counter-detector O are selected.
  • the ⁇ -ray detector 3 is referred to as a correction target detector C.
  • a time difference histogram relating to the selected reference detector S which is the “1” -th ⁇ -ray detector 3 and the “61” -th ⁇ -ray detector 3 (opposite detector O) opposed to the reference detector S by 180 ° is obtained.
  • the standard One of the two ⁇ -ray detectors 3 including the reference detector S which is the selected “1” -th ⁇ -ray detector 3 and the “61” -th ⁇ -ray detector 3 (opposite detector O).
  • the reference detector S is selected as the ⁇ -ray detector 3, and the above-mentioned “60” -th and “62” -th correction target detectors C different from the opposing detector O which is the other ⁇ -ray detector 3. Select.
  • the correction target detector C detects the “120” th counter detection with respect to the reference detector S which is the “60th” ⁇ -ray detector 3.
  • the “119” -th ⁇ -ray detector 3 adjacent to the detector O and adjacent to the “2” -th counter detector O with respect to the reference detector S which is the “62” -th ⁇ -ray detector 3 This is the “3” th ⁇ -ray detector 3.
  • the “1” -th ⁇ -ray detector 3 adjacent to the “120” -th counter detector O and adjacent to the “2” -th counter detector O is already selected and thus excluded.
  • the “3” -th and “119” -th ⁇ -ray detectors 3 of the correction target detector C are newly set as a reference.
  • the correction target detector C detects the “63” th counter detection with respect to the reference detector S which is the “3” th ⁇ -ray detector 3.
  • the “64” -th ⁇ -ray detector 3 adjacent to the detector O and adjacent to the “59” -th counter detector O with respect to the reference detector S which is the “119” -th ⁇ -ray detector 3 This is the “58” -th ⁇ -ray detector 3.
  • the “60” -th ⁇ -ray detector 3 adjacent to the “59” -th counter detector O and the “62” -th ⁇ -ray detector 3 adjacent to the “63” -th counter detector O are as follows. Exclude it because it is already selected. After the correction, the “58” th and “64th” ⁇ -ray detectors 3 of the correction target detector C are newly set as a reference.
  • the series of timing corrections shown in FIG. 4 is preferably performed periodically according to the secular change of the time delay of each channel. Further, the correction table is filed and written and stored in the storage medium of the data collection / control unit 11 or the memory unit 8 described above. When the PET apparatus is turned off and turned on again, the storage medium or It is preferable to read from the memory unit 8 and rewrite it in the timing correction table 9b.
  • the data collection / control unit 11 performs the following correction calculation in steps S3 to S6. That is, two ⁇ -ray detectors 3 to be simultaneously counted are selected (in the case of FIG. 7A, the “61” th reference detector S and the “61” th reference detector 180 ° opposite thereto.
  • the counter detector O selects the reference detector S which is one of the two selected ⁇ -ray detectors 3 and a ⁇ -ray detector 3 different from the other counter-detector O.
  • the time difference histograms related to the two ⁇ -ray detectors 3 selected this time are corrected based on the reference.
  • the coincidence circuit 10 simultaneously counts radiation (gamma rays in the first embodiment) based on the time difference histogram for each pair of each ⁇ -ray detector 3 that has been repeatedly corrected as described above. be able to. As a result, it is possible to accurately perform coincidence counting and timing correction without repeating many measurements and calculations.
  • the time difference at which the total count value in the time difference histogram based on the reference is the maximum is used as the reference value, and the time difference in the time difference histogram is corrected based on the reference value.
  • the time point at which the total count value in the time difference histogram is maximum is the timing at which coincidence counting is most likely to occur. Therefore, by correcting the time difference of the time difference histogram based on the reference value that is the timing, it is possible to align at that timing.
  • the ⁇ -ray detector 3 is a DOI detector configured by laminating the scintillator blocks 31 in the ⁇ -ray depth direction.
  • detection time correction information corresponding to coordinate information such as the depth direction in which the interaction occurs as in Patent Document 4 described above is required.
  • the detection time correction information can be obtained. That is, after performing a series of timing corrections shown in FIG. 4, the light source position in the depth direction is specified by the center of gravity calculation, and the detection time correction information corresponding to the coordinate information such as the depth direction is obtained.
  • an external radiation source RI that irradiates the same kind of radiation as the radiopharmaceutical or a phantom Ph that radiates the same kind of radiation as the radiopharmaceutical from the inside is provided, and the above time difference histogram is obtained from the external radiation source RI or phantom Ph. Obtained based on radiation.
  • FIG. 9 is a front view showing one embodiment of switching between the reference scintillator element unit and the correction target scintillator element unit.
  • the PET apparatus according to the second embodiment including the third and fourth embodiments described later is the block diagram shown in FIG. Moreover, the same code
  • one scintillator block 31 of a certain ⁇ -ray detector is set as a reference, and one scintillator block 31 set as the reference is set as a reference scintillator element unit S.
  • One scintillator block 31 of a ⁇ -ray detector different from the reference scintillator element unit S is set as a correction target, and one scintillator block 31 set as the correction target is set as a correction target scintillator element unit C.
  • the specific timing correction is the same method except that it is changed from the detector unit of the first embodiment to the scintillator element unit of the second embodiment, and thus the description thereof is omitted.
  • the corrected scintillator element unit C after correction is a new reference
  • the reference scintillator element unit S is shown in FIG. 9B.
  • One scintillator block 31 of a ⁇ -ray detector different from the reference scintillator element unit S is set as a correction target, and one scintillator block 31 set as the correction target is set as a correction target scintillator element unit C.
  • the corrected scintillator element unit C after correction is set as a reference scintillator element unit S as shown in FIG.
  • One scintillator block 31 of a ⁇ -ray detector different from the reference scintillator element unit S is set as a correction target, and one scintillator block 31 set as the correction target is set as a correction target scintillator element unit C. In this manner, timing correction is performed in the same manner in the scintillator element units of other detectors.
  • two scintillator element units of the ⁇ -ray detector to be simultaneously counted are selected, and the two selected scintillator elements Among the units, a scintillator element unit of one ⁇ -ray detector (reference scintillator element unit S in the case of FIG. 9) is selected, and a scintillator element unit different from the scintillator element unit of the other ⁇ -ray detector ( In the case of FIG.
  • FIG. 10 is a front view showing one embodiment of switching between the reference scintillator element group and the correction target scintillator element unit.
  • the PET apparatus according to the third embodiment including the fourth embodiment which will be described later is the block diagram shown in FIG.
  • the same reference numerals are given to portions common to the above-described first and second embodiments, and the description thereof is omitted.
  • correction is performed for each scintillator element group (that is, in units of detectors) composed of a plurality of scintillator elements (scintillator blocks 31), but in the third embodiment, as in the second embodiment described above.
  • the accuracy is further improved by performing correction for each scintillator element unit composed of one scintillator element as described below.
  • FIG. 9 of the second embodiment only the scintillator block 31 is shown in FIG. 10, and the other structures (light guide 32 and PMT 33) are not shown.
  • a plurality of scintillator blocks 31 of a certain ⁇ -ray detector are set as a reference, and the plurality of scintillator blocks 31 set as the reference are set as a reference scintillator element group S.
  • One scintillator block 31 of a ⁇ -ray detector different from the reference scintillator element group S is set as a correction target, and one scintillator block 31 set as the correction target is set as a correction target scintillator element unit C.
  • the specific timing correction is the same method except that it is changed from the detector unit of the first embodiment to the scintillator element unit of the third embodiment, and thus the description thereof is omitted.
  • the previously selected reference scintillator element group S Based on A scintillator block 31 different from the correction target scintillator element unit C shown in FIG. 10A is set as a correction target, and one scintillator block 31 set as the correction target is set as a correction target scintillator element unit C.
  • the total number of corrected correction target scintillator element units C is integrated into a new standard.
  • a reference scintillator element group S is assumed.
  • a ⁇ -ray detector different from the reference scintillator element unit S (in FIG. 9, a ⁇ -ray detector adjacent to the ⁇ -ray detector of the reference scintillator element group S in FIGS. 9A and 9B).
  • One scintillator block 31 is set as a correction target, and one scintillator block 31 set as the correction target is set as a correction target scintillator element unit C. In this manner, timing correction is performed in the same manner in the scintillator element units of other detectors.
  • the scintillator element group of one ⁇ -ray detector (in the case of FIG. 10).
  • the scintillator element group and the scintillator element group selected this time based on the time difference histogram.
  • the operation based on the time difference histogram relating to the corrected scintillator element group and scintillator element unit is newly performed. Since the correction is performed for each scintillator element group and each scintillator element in this way, the accuracy can be further improved as compared with the case where the correction is performed for each detector as in the first embodiment. In addition, the calculation time and burden can be reduced as compared with the case where correction is performed for each scintillator element as in the second embodiment.
  • FIG. 11 is an explanatory diagram of a time difference histogram for a scintillator element having self-radiation.
  • the PET apparatus according to the fourth embodiment is a block diagram shown in FIG.
  • the same reference numerals are given to portions common to the above-described first to third embodiments, and description thereof is omitted.
  • the scintillator block is a substance to which self-radiation (an element that simultaneously emits a plurality of radiations) represented by Lu-176 or the like is added.
  • the scintillator block is a scintillator element having self-radiation
  • the detector includes a scintillator element having self-radiation.
  • a scintillator block may be configured by attaching a thin film tape made of a substance to which self-radiation is added to a crystal element having no self-radiation (for example, GSO).
  • a scintillator block may be configured by applying a coating agent made of the above substance to a crystal element having no self-radiation.
  • the scintillator element with self-radiation includes nuclides that cause ⁇ or ⁇ -decay and emit ⁇ rays accompanying the ⁇ or ⁇ -decay.
  • the time difference histogram is acquired based on the radiation from the self-radioactivity (here, ⁇ rays).
  • the time difference histogram is a time stamp difference (ie, time difference) for each pair of detectors 3 (indicated as “Difference Time” in FIG. 11). Is a distribution of count values with respect to a change in time difference when a count value (indicated as “Event Counts” in FIG. 11) is plotted on the vertical axis.
  • the time difference histogram shows a time point where the total count value is maximum and a time point where the total count value is the second largest. Therefore, in the fourth embodiment, a time difference that is an intermediate value between the time difference at which the total count value is maximum in the reference time difference histogram and the time difference at which the total count value is the second largest is used as the reference value. This reference value is set to “0”.
  • the time difference histogram to be corrected is shifted from the reference time difference histogram and is shifted left and right as shown by the solid line in FIG. Therefore, a shift amount for returning “0” shown by the dotted line from the solid line in FIG. 11B to the reference value is obtained as a correction amount for the time difference histogram to be corrected. Then, the time difference histogram is corrected by applying this correction amount.
  • an intermediate value between the time difference at which the total count value is maximum in the reference time difference histogram and the time difference at which the total count value is the second largest is obtained.
  • a certain time difference is used as a reference value, and the time difference in the time difference histogram is corrected based on the above-described reference value.
  • the scintillator element of the other detector detects the most radiation from the self-radiation of the scintillator element of one detector.
  • the timing at which the scintillator element of one detector detects the most radiation from the self-radioactivity of the scintillator element of the other detector is the time point where the total count value is maximum, or the total count value is This is the second largest time point. Therefore, the intermediate value between the two time differences, which are these timings, is the timing at which coincidence counting is most likely to occur. Therefore, by correcting the time difference of the time difference histogram based on the reference value of the intermediate value that is the timing, it is possible to align at that timing.
  • the present invention is not limited to the above embodiment, and can be modified as follows.
  • the positron CT apparatus PET apparatus
  • PET apparatus positron CT apparatus
  • ⁇ rays are taken as an example of radiation, but ⁇ rays, ⁇ rays, or the like may be used.
  • the detector is provided with a scintillator element having self-radiation as in the fourth embodiment, in the case where the detector of the scintillator element in which ⁇ or ⁇ -decay has occurred detects ⁇ rays or ⁇ rays.
  • a time difference histogram may be used.
  • the DOI detector is used.
  • the present invention can also be applied to a detector that does not distinguish the depth direction.
  • the present invention can also be applied to a detector having a structure including a single scintillator element.
  • the present invention is applied to detectors arranged in a ring shape, but the present invention can be applied even when a plurality of detectors are provided without being installed in a ring shape. .

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Medical Informatics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (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

ポジトロンＣＴ装置およびタイミング補正方法
　この発明は、被検体内に投与されたポジトロン放射性薬剤から放出される放射線を検出するポジトロンＣＴ装置およびタイミング補正方法に関する。
　ポジトロンＣＴ装置、すなわちＰＥＴ(Positron Emission Tomography)装置は、陽電子（Positron）、すなわちポジトロンの消滅によって発生する複数本のγ線を検出して複数個の検出器でγ線を同時に検出したときのみ（つまり同時計数したときのみ）被検体の画像を再構成するように構成されている。
　同時計数するには同時計数回路を用いる。ところが、各検出器から同時計数回路までの信号チャンネルでの信号の時間遅れが生じる。この時間遅れは、各信号チャンネルでバラツキがある。したがって、各信号チャンネルでの遅れ時間を調整してγ線が同時計数回路に到達するタイミングを一致させるタイミング補正を行う必要がある。
　そこで、かかるタイミング補正を行うのに、校正用の放射線源（外部線源）や模擬信号を用いて校正データを取得し、その校正データに基づいて信号転送の時間的なバラツキを調整する（例えば、特許文献１～３参照）。なお、近年では、ポジトロンの対消滅発生地点を限定する時間差情報（飛行時間）(TOF: Time Of Flight)を利用した技術が提案されている（例えば、特許文献４参照）。ＴＯＦは、消滅放射線が光速であることを用いて、対消滅発生地点から検出器に到達する時間差を、対消滅発生地点から検出器のシンチレータ素子による光源発生位置までの距離差に換算することで、対消滅発生地点を求める技術である。
　上述の特許文献１に記載された信号のタイミング補正方法について述べる。放射線源から照射された放射線を検出器が検出し、放射線が検出器に入射された放射線入射タイミングを表すタイミング信号を、遅延調整回路を経て同時計数回路に入力する。このタイミング信号の入力を受けて同時計数回路の出力を測定して、各信号チャンネル毎に放射線の感度（すなわち計数）を測定する。この後、遅延調整回路で調整される遅延量を変えながら上述の感度を測定して、遅延量変化に対する感度分布を求める。測定された感度が最も高くなる遅延量を遅延調整回路に用いることで、信号の時間遅れを補正する。
　上述の特許文献２に記載された信号のタイミング補正方法について述べる。校正用の放射線源（外部線源）をＰＥＴ装置の視野(FOV: Field of View)内に設置する。ここで、複数の検出器をリング（環状）に並べる。ある検出器を基準とすると、その基準の検出器と共通の視野を有する複数の検出器のタイミング値を平均して、その平均されたタイミング値を基準の検出器に対する時間遅延値として求める。基準の検出器が隣接する検出器を新たに基準として、同様に時間遅延値を求め、最初に求められた時間遅延値と次に求められた時間遅延値との差分を基準補正値として求める。その基準補正値を用いて時間を揃えることでタイミング補正を行う。以下、順に隣接する検出器に同様の演算を行うことで、リング上を１周したときに全ての検出器に関するタイミング補正を行う。
　上述の特許文献３に記載された信号のタイミング補正方法について述べる。模擬信号発生装置から出力された模擬信号を複数の信号処理装置（信号処理ユニット）にそれぞれ入力し、各信号処理装置の出力に基づいて校正データを生成して、タイミング補正を行う。
　また、上述の特許文献４では、相互作用を起こした深さ方向の光源位置(DOI: Depth of Interaction)を弁別することができるＤＯＩ検出器をＴＯＦ型ＰＥＴ装置に組み込む。ＤＯＩ検出器は、各々のシンチレータ素子を放射線（ここではγ線）の深さ方向に積層して構成されたものであり、相互作用を起こした深さ方向と横方向（入射面に平行な方向）との座標情報を重心演算により求める。この座標情報に対応する検出時刻補正情報をテーブルに書き込んで記憶し、その検出時刻補正情報を参照することにより、飛行時間差の情報精度を高めている。
特公平６－１９４３６号公報 特許第３３４３１２２号明細書 特開２００６－９０８２７号公報 特開２００８－５１７０１号公報
　しかしながら、上述の特許文献１の場合、遅延量を求めるために測定を繰り返さなければならない。また、上述の特許文献２の場合、基準の検出器と共通の視野を有する複数の検出器間では時間が揃っていない。したがって、それら検出器のタイミング値を平均して、その平均されたタイミング値を時間遅延値として用いて、リング上を１周しても時間が完全に揃わない。その結果、最適な時間遅延値を求めるまでにリング上を複数周（例えば２～３周）分にわたって廻って上述の演算を繰り返す必要がある。
　また、上述の特許文献３の場合、検出器以外に模擬信号発生装置が必要となる。また、信号処理装置での遅延のみを補正するので、上述の特許文献４のような相互作用を起こした深さ方向の光源位置を弁別するＤＯＩ検出器に対して消滅放射線とシンチレータ素子とが相互作用を起こした深さ方向および横方向の座標情報に対応する検出時刻補正情報が必要な場合に、検出時刻補正情報を求めることができない。
　この発明は、このような事情に鑑みてなされたものであって、測定や演算を多数繰り返さずに、かつ正確に同時計数することができるポジトロンＣＴ装置およびタイミング補正方法を提供することを目的とする。
　この発明は、このような目的を達成するために、次のような構成をとる。
　すなわち、この発明のポジトロンＣＴ装置は、被検体内に投与されたポジトロン放射性薬剤から放出される放射線を検出する複数の検出器を備えたポジトロンＣＴ装置であって、放射線を同時計数する各検出器の対ごとの時間差変化に対する計数値分布を表す時間差ヒストグラムに関して、同時計数する対象の検出器を２つ選択して、当該選択された２つの検出器のうち、一方の検出器を選択するとともに、他方の検出器とは別の検出器を選択し、その選択を繰り返し行った際に、過去に選択された２つの検出器に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つの検出器に関する時間差ヒストグラムを補正し、当該補正された２つの検出器に関する時間差ヒストグラムを新たに基準とする作業を繰り返し行う演算手段と、前記演算手段で繰り返し補正された各検出器の対ごとの前記時間差ヒストグラムに基づいて、放射線を同時計数する同時計数回路とを備えることを特徴とするものである。
　この発明のポジトロンＣＴ装置によれば、放射線を同時計数する各検出器の対ごとの時間差変化に対する計数値分布を表す時間差ヒストグラムに関して、演算手段は以下の演算を行う。すなわち、同時計数する対象の検出器を２つ選択して、当該選択された２つの検出器のうち、一方の検出器を選択するとともに、他方の検出器とは別の検出器を選択し、その選択を繰り返し行った際に、過去に選択された２つの検出器に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つの検出器に関する時間差ヒストグラムを補正する。そして、当該補正された２つの検出器に関する時間差ヒストグラムを新たに基準とする作業を繰り返し行う。このように、過去に選択された２つの検出器に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つの検出器に関する時間差ヒストグラムを補正し、当該補正された２つの検出器に関する時間差ヒストグラムを新たに基準とする作業を繰り返し行うことで、例えば上述した特許文献２のように複数の検出器のタイミング値を平均した場合よりも、最適な時間差ヒストグラムを求める際に収束性がよい。したがって、測定や演算を多数繰り返さずに最適な時間差ヒストグラムを求めることができる。また、上述の繰り返し補正された各検出器の対ごとの時間差ヒストグラムに基づいて、同時計数回路は放射線を同時計数するので、正確に同時計数することができる。その結果、測定や演算を多数繰り返さずに、かつ正確に同時計数することができる。
　また、この発明のタイミング補正方法は、被検体内に投与されたポジトロン放射性薬剤から放出される放射線を同時計数するために用いられるタイミング補正方法であって、放射線を同時計数する各検出器の対ごとの時間差変化に対する計数値分布を表す時間差ヒストグラムに関して、同時計数する対象の検出器を２つ選択して、当該選択された２つの検出器のうち、一方の検出器を選択するとともに、他方の検出器とは別の検出器を選択し、その選択を繰り返し行った際に、過去に選択された２つの検出器に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つの検出器に関する時間差ヒストグラムを補正し、当該補正された２つの検出器に関する時間差ヒストグラムを新たに基準とする作業を繰り返し行うヒストグラム補正工程を備えることを特徴とするものである。
　この発明のタイミング補正方法によれば、各検出器の対ごとの時間差ヒストグラムに関して、ヒストグラム補正工程では以下の補正を行う。すなわち、同時計数する対象の検出器を２つ選択して、当該選択された２つの検出器のうち、一方の検出器を選択するとともに、他方の検出器とは別の検出器を選択し、その選択を繰り返し行った際に、過去に選択された２つの検出器に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つの検出器に関する時間差ヒストグラムを補正する。そして、当該補正された２つの検出器に関する時間差ヒストグラムを新たに基準とする作業を繰り返し行う。その結果、測定や演算を多数繰り返さずに、かつ正確にタイミング補正を行うことができる。
　上述したこれらのポジトロンＣＴ装置およびタイミング補正方法の一例は、上述の基準とした時間差ヒストグラムにおける総計数値が最大となる時間差を基準値とし、上述の基準値に基づいて時間差ヒストグラムの時間差を補正することである。時間差ヒストグラムにおける総計数値が最大となる時間的な箇所は、同時計数が最も起こりうるタイミングである。したがって、そのタイミングである基準値に基づいて時間差ヒストグラムの時間差を補正することで、そのタイミングで揃えることができる。
　上述したこれらのポジトロンＣＴ装置およびタイミング補正方法の他の一例は、上述の基準とした時間差ヒストグラムにおける総計数値が最大となる時間差と、総計数値が二番目に大きい時間差との中間値にある時間差を基準値とし、上述の基準値に基づいて時間差ヒストグラムの時間差を補正することである。後述するように、自己放射能を持つシンチレータ素子を検出器が備えた場合に、一方の検出器のシンチレータ素子の自己放射能からの放射線を他方の検出器のシンチレータ素子が最も多く検出したタイミング、および他方の検出器のシンチレータ素子の自己放射能からの放射線を一方の検出器のシンチレータ素子が最も多く検出したタイミングが、総計数値が最大となる時間的な箇所、あるいは総計数値が二番目に大きい時間的な箇所である。このことから、これらのタイミングである２つの時間差の中間値は、同時計数が最も起こりうるタイミングである。したがって、そのタイミングである中間値の基準値に基づいて時間差ヒストグラムの時間差を補正することで、そのタイミングで揃えることができる。
　上述の検出器は、単数のシンチレータ素子を備えた構造の他に、複数のシンチレータ素子を備えた構造もある。後者の構造の場合には、複数のシンチレータ素子からなるシンチレータ素子群（すなわち検出器単位）ごとに補正を行ってもよいが、下記のように１つのシンチレータ素子からなるシンチレータ素子単位ごとに補正を行って、精度をより向上させてもよい。
　例えば、放射線を同時計数する各検出器の１つのシンチレータ素子からなるシンチレータ素子単位の対ごとの時間差ヒストグラムに関して、以下の演算／補正を行う。すなわち、同時計数する対象の検出器のシンチレータ素子単位を２つ選択して、当該選択された２つのシンチレータ素子単位のうち、一方の検出器のシンチレータ素子単位を選択するとともに、他方の検出器のシンチレータ素子単位とは別のシンチレータ素子単位を選択し、その選択を繰り返し行った際に、過去に選択された２つの検出器のシンチレータ素子単位に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つの検出器のシンチレータ素子単位に関する時間差ヒストグラムを補正する。そして、当該補正された２つの検出器のシンチレータ素子単位に関する時間差ヒストグラムを新たに基準とする作業を繰り返し行う。このようにシンチレータ素子単位ごとに補正を行うので、検出器単位ごとに補正を行うときと比較して精度をより向上させることができる。
　また例えば、放射線を同時計数する各検出器のうち、一方の検出器の複数のシンチレータ素子からなるシンチレータ素子群と他方の検出器の１つのシンチレータ素子からなるシンチレータ素子単位との対ごとの時間差ヒストグラムに関して、以下の演算／補正を行う。すなわち、同時計数する対象の検出器のうち、一方の検出器のシンチレータ素子群と他方の検出器のシンチレータ素子単位とを選択して、当該選択されたシンチレータ素子群とシンチレータ素子単位とのうち、当該シンチレータ素子群を選択するとともに、他方の検出器のシンチレータ素子単位とは別のシンチレータ素子単位を選択し、その選択を繰り返し行った際に、過去に選択されたシンチレータ素子群とシンチレータ素子単位とに関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択されたシンチレータ素子群とシンチレータ素子単位とに関する時間差ヒストグラムを補正する。そして、当該補正されたシンチレータ素子群とシンチレータ素子単位とに関する時間差ヒストグラムを新たに基準とする作業を繰り返し行う。このようにシンチレータ素子群およびシンチレータ素子単位ごとに補正を行うので、検出器単位ごとに補正を行うときと比較して精度をより向上させることができる。また、シンチレータ素子単位ごとに補正を行うときと比較して演算時間や負担を軽減させることができる。
　また、ポジトロンＣＴ装置において、上述の検出器は、各々のシンチレータ素子を放射線の深さ方向に積層して構成されたＤＯＩ検出器であってもよい。この発明をＤＯＩ検出器に適用することで、例えば上述の特許文献４のように相互作用を起こした深さ方向などの座標情報に対応する検出時刻補正情報が必要な場合に、検出時刻補正情報を求めることができる。
　上述したこれらのポジトロンＣＴ装置およびタイミング補正方法では、上述の放射性薬剤と同種の放射線を照射する外部線源、または放射性薬剤と同種の放射線を内部から照射するファントムを備え、上述の時間差ヒストグラムを外部線源またはファントムからの放射線に基づいて取得してもよいし、検出器は、自己放射能を持つシンチレータ素子を備え、上述の時間差ヒストグラムを自己放射能からの放射線に基づいて取得してもよい。
　この発明に係るポジトロンＣＴ装置およびタイミング補正方法によれば、同時計数する対象の検出器を２つ選択して、当該選択された２つの検出器のうち、一方の検出器を選択するとともに、他方の検出器とは別の検出器を選択し、その選択を繰り返し行った際に、過去に選択された２つの検出器に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つの検出器に関する時間差ヒストグラムを補正する。そして、当該補正された２つの検出器に関する時間差ヒストグラムを新たに基準とする作業を繰り返し行う。その結果、測定や演算を多数繰り返さずに、かつ正確にタイミング補正を行うことができる。
各実施例に係るＰＥＴ(Positron Emission Tomography)装置の側面図およびブロック図である。 γ線検出器の概略斜視図である。 ＰＥＴ装置におけるリング状に配置されたγ線検出器の正面図およびそれに関するブロック図である。 一連のタイミング補正方法の流れを示すフローチャートである。 （ａ）は外部線源を設置したときのγ線検出器の正面図、（ｂ）はファントムを設置したときのγ線検出器の正面図である。 （ａ）、（ｂ）は、時間差ヒストグラムの説明図である。 （ａ）～（ｃ）は、基準検出器と補正対象検出器との切り替えの一実施態様を示す正面図である。 （ａ）～（ｃ）は、図７とは別の実施態様で、基準検出器と補正対象検出器との切り替えの一実施態様を示す正面図である。 （ａ）～（ｃ）は、基準シンチレータ素子単位と補正対象シンチレータ素子単位との切り替えの一実施態様を示す正面図である。 （ａ）～（ｃ）は、基準シンチレータ素子群と補正対象シンチレータ素子単位との切り替えの一実施態様を示す正面図である。 （ａ）、（ｂ）は、自己放射能を持つシンチレータ素子のときの時間差ヒストグラムの説明図である。
　１０　…　同時計数回路
　１１　…　データ収集・制御部
　３　…　γ線検出器
　３１　…　シンチレータブロック
　３３　…　光電子増倍管（ＰＭＴ）
　ＲＩ　…　外部線源
　Ｐｈ　…　ファントム
　Ｓ　…　基準検出器、基準シンチレータ素子単位、基準シンチレータ素子群
　Ｃ　…　補正対象検出器、補正対象シンチレータ素子単位
　Ｏ　…　対向検出器
　Ｍ　…　被検体
　以下、図面を参照してこの発明の実施例1を説明する。図１は、実施例１に係るＰＥＴ(Positron Emission Tomography)装置の側面図およびブロック図であり、図２は、γ線検出器の概略斜視図であり、図３は、ＰＥＴ装置におけるリング状に配置されたγ線検出器の正面図およびそれに関するブロック図である。
　後述する実施例２～４も含めて、本実施例１に係るＰＥＴ装置は、図１に示すように、被検体Ｍを載置する天板１を備えている。この天板１は、上下に昇降移動、被検体Ｍの体軸Ｚに沿って平行移動するように構成されている。このように構成することで、天板１に載置された被検体Ｍは、後述するガントリ２の開口部２ａを通って、頭部から順に腹部、足部へと走査されて、被検体Ｍの画像を得る。なお、走査される部位や各部位の走査順序については特に限定されない。
　天板１の他に、本実施例１に係るＰＥＴ装置は、開口部２ａを有したガントリ２と、γ線検出器３とを備えている。γ線検出器３は、被検体Ｍの体軸Ｚ周りを取り囲むようにしてリング状に配置されており、ガントリ２内に埋設されている。γ線検出器３は、この発明における検出器に相当する。
　その他にも、本実施例１に係るＰＥＴ装置は、天板駆動部４とコントローラ５と入力部６と出力部７とメモリ部８と検出器信号処理部９と同時計数回路１０とデータ収集・制御部１１と再構成処理部１２とを備えている。天板駆動部６は、天板１の上述した移動を行うように駆動する機構であって、図示を省略するモータなどで構成されている。同時計数回路１０は、この発明における同時計数回路に相当し、データ収集・制御部１１は、この発明における演算手段に相当する。
　コントローラ５は、本実施例１に係るＰＥＴ装置を構成する各部分を統括制御する。コントローラ５およびデータ収集・制御部１１は、中央演算処理装置（ＣＰＵ）などで構成されている。
　入力部６は、オペレータが入力したデータや命令をコントローラ５に送り込む。入力部６は、マウスやキーボードやジョイスティックやトラックボールやタッチパネルなどに代表されるポインティングデバイスで構成されている。出力部７はモニタなどに代表される表示部やプリンタなどで構成されている。
　メモリ部８は、ＲＯＭ（Read-only Memory）やＲＡＭ（Random-Access Memory）などに代表される記憶媒体で構成されている。本実施例１では、同時計数回路１０で同時計数された計数値（カウント）や同時計数した２つのγ線検出器３からなる検出器対やＬＯＲといった同時計数に関するデータや、再構成処理部１２で処理された画像などについてはＲＡＭに書き込んで記憶し、必要に応じてＲＡＭから読み出す。ＲＯＭには、各種の核医学診断を含めて撮像を行うためのプログラム等を予め記憶しており、そのプログラムをコントローラ５およびデータ収集・制御部１１が実行することでそのプログラムに応じた核医学診断をそれぞれ行う。なお、ＬＯＲ(Line Of Response)とは、同時計数する２つのγ線検出器３間で結ぶ仮想上の直線のことである。
　再構成処理部１２は、例えば上述したメモリ部８などに代表される記憶媒体のＲＯＭに記憶されたプログラムあるいは入力部６などに代表されるポインティングデバイスで入力された命令をコントローラ５が実行することで実現される。
　放射性薬剤が投与された被検体Ｍから発生したγ線をγ線検出器３のシンチレータブロック３１（図２を参照）が光に変換して、変換されたその光をγ線検出器３の光電子増倍管(PMT: Photo Multiplier Tube)３３（図２を参照）は増倍させて電気信号に変換する。その電気信号を検出器信号処理部９を介して同時計数回路１０に送り込む。
　具体的には、被検体Ｍに放射性薬剤を投与すると、ポジトロン放出型のＲＩのポジトロンが消滅することにより、２本のγ線が発生する。同時計数回路１０は、シンチレータブロック３１（図２を参照）の位置とγ線の入射タイミングとをチェックし、被検体Ｍの両側にある２つのシンチレータブロック３１でγ線が同時に入射したときのみ、送り込まれたイベントを適正なデータと判定する。一方のシンチレータブロック３１のみにγ線が入射したときには、同時計数回路１０は棄却する。つまり、同時計数回路１０は、上述した電気信号に基づいて、２つのγ線検出器３においてγ線が同時観測されたことを検出する。
　同時計数回路１０に送り込まれた電気信号のうち、同時観測（すなわち同時計数）された電気信号を画像情報として判定し、その画像情報をデータ収集・制御部１１を介して再構成処理部１２に送り込む。再構成処理部１２は順投影処理や逆投影処理による画像再構成を行って、被検体Ｍの画像を求める。画像を、コントローラ５を介して出力部７に送り込む。このようにして、再構成処理部１２で得られた画像に基づいて核医学診断を行う。再構成処理部１２での画像再構成については、ＤＲＡＭＡ法（Dynamic Row-Action Maximum Likelihood Algorithm）などの公知の逐次近似法アルゴリズムを適用する。検出器信号処理部９およびデータ収集・制御部１１の具体的な機能については後述する。
　γ線検出器３は、図２に示すようにシンチレータブロック３１と、そのシンチレータブロック３１に対して光学的に結合されたライトガイド３２と、そのライトガイド３２に対して光学的に結合された光電子増倍管（以下、単に「ＰＭＴ」と略記する）３３とを備えている。シンチレータブロック３１を構成する各シンチレータ素子は、γ線の入射に伴って発光することでγ線から光に変換する。この変換によってシンチレータ素子はγ線を検出する。シンチレータ素子において発光した光がシンチレータブロック３１で十分に拡散されて、ライトガイド３２を介してＰＭＴ３３に入力される。ＰＭＴ３３は、シンチレータブロック３１で変換された光を増倍させて電気信号に変換する。その電気信号は、上述したように検出器信号処理部９を介して同時計数回路１０に送り込まれる。シンチレータブロック３１は、この発明におけるシンチレータ素子に相当し、ＰＭＴ３３は、この発明における光電変換手段に相当する。
　以上をまとめると、図２に示すγ線検出器３は、γ線の入射により蛍光する複数のシンチレータブロック３１と、各シンチレータブロック３１からの光を光電変換することでγ線を検出するＰＭＴ３３とを備えている。また、図２に示すγ線検出器３は、各々のシンチレータブロック３１をγ線の深さ方向に積層（図２では４層に積層）して構成されたＤＯＩ検出器である。
　なお、同時計数回路１０で同時計数するときには、上述したように各γ線検出器３から同時計数回路１０までの信号チャンネルでの信号の時間遅れが生じて同時計数が正確に行えない。そこで、被検体Ｍを用いた核医学診断よりも前に予め外部線源あるいはファントムを用いて得られた補正用データを、図３に示すように検出器信号処理部９および同時計数回路１０を介して、データ収集・制御部１１に送り込んで、データ収集・制御部１１がデータ収集を行う。そして、その収集された補正用データを検出器信号処理部９のタイミング補正テーブル９ｂに帰還制御（フィードバック制御）することで、データ収集・制御部１１はタイミング補正を行う。
　通常の被検体Ｍを用いた核医学診断での撮像では、放射性薬剤が投与された被検体Ｍから発生したγ線をγ線検出器３が検出して、検出器信号処理部９、同時計数回路１０を経て、データ収集・制御部１１を介して、タイミング補正されることなく、再構成処理部１２に送り込む。一方、外部線源あるいはファントムを用いたデータ収集では、外部線源あるいはファントムからのγ線をγ線検出器３が検出して、検出器信号処理部９、同時計数回路１０を経て、データ収集・制御部１１に送り込んで、データ収集・制御部１１が検出器信号処理部９のタイミング補正テーブル９ｂにフィードバック制御してタイミング補正を行って、各々の同時計数回路１０での遅延量を調整する。
　検出器信号処理部９は、γ線検出器３のＰＭＴ３３から出力された電気信号に基づいて、γ線入射タイミングを表す信号（以下、「タイムスタンプ」とも呼ぶ）を作成するタイミング信号作成回路９ａと、上述のタイミング補正テーブル９ｂとを備えている。上述したように、各々のγ線検出器３はリング状に配置されており、各γ線検出器３は検出器信号処理部９のタイミング信号作成回路９ａに接続されている（図３では、タイミング信号作成回路９ａを２つのみ図示）。タイミング信号作成回路９ａはタイミング補正テーブル９ｂに接続され、タイミング補正テーブル９ｂは同時計数回路１０に接続されている。同時計数回路１０に接続されたデータ収集・制御部１１をタイミング補正テーブル９ｂに接続することで、データ収集・制御部１１はタイミング補正テーブル９ｂにフィードバック制御する。
　次に、タイミング補正について、図４～図８を参照して説明する。図４は、一連のタイミング補正方法の流れを示すフローチャートであり、図５は、外部線源あるいはファントムを設置したときのγ線検出器の正面図であり、図６は、時間差ヒストグラムの説明図であり、図７は、基準検出器と補正対象検出器との切り替えの一実施態様を示す正面図であり、図８は、図７とは別の実施態様で、基準検出器と補正対象検出器との切り替えの一実施態様を示す正面図である。
　（ステップＳ１）補正用データの測定
　図５（ａ）に示すように、放射性薬剤、すなわち放射性同位元素（ＲＩ）と同種の放射線（本実施例１ではγ線）を照射する外部線源ＲＩをＰＥＴ装置の視野内に設置する。好ましくは、外部線源ＲＩを視野内の中央領域に設置することで、約１８０°対向した２つのγ線検出器３が外部線源ＲＩからの放射線をほぼ同時に検出することができ、各γ線検出器３から同時計数回路１０までの信号チャンネルでの信号の時間遅れのみがタイミング補正の対象となる。したがって、タイミング補正を行うのみでγ線が同時計数回路１０に到達するタイミングを一致させることができる。なお、外部線源ＲＩ以外にも、図５（ｂ）に示すように、放射性薬剤と同種の放射線を内部から照射するファントムＰｈをＰＥＴ装置の視野内に設置してもよい。ファントムＰｈを設置する場合においても、好ましくは視野内の中央領域に設置する。外部線源ＲＩは、この発明における外部線源に相当し、ファントムＰｈは、この発明におけるファントムに相当する。
　外部線源ＲＩまたはファントムＰｈからの放射線をγ線検出器３が検出して補正用データとして取得し、検出器信号処理部９を介して同時計数回路１０に送り込む。図５に示すように全てのγ線検出器３で補正用データを測定して取得する。
　（ステップＳ２）時間差ヒストグラムの作成
　ステップＳ１で測定された補正用データに基づいて、放射線を同時計数する各γ線検出器３の対ごとに時間差ヒストグラムを作成する。図６に示すように、時間差ヒストグラムは、各γ線検出器３の対ごとにタイムスタンプの差分（すなわち時間差）（図６では「Difference Time」で表記）を横軸にとり、計数値（図６では「Event Counts」で表記）を縦軸にとったときの時間差変化に対する計数値分布である。
　外部線源ＲＩやファントムＰｈを視野の中央領域に設置すると、時間差ヒストグラムは、図６（ａ）に示すように“０”を中心とした度数分布として得られる。しかし、時間遅れによるズレがある場合には、時間差ヒストグラムは、図６（ｂ）の実線に示すように左右にシフトする。そこで、時間差ヒストグラムを、図６（ｂ）の実線から点線に示す“０”を中心とした度数分布に戻すようなシフト量を補正量と考える。
　（ステップＳ３）基準検出器・補正対象検出器の設定
　次に、基準検出器と補正対象検出器とをそれぞれ設定する。図７（ａ）に示すように、１つのγ線検出器３を基準として設定し、その基準として設定されたγ線検出器３（図７（ａ）の黒塗りを参照）を基準検出器Ｓとする。この基準検出器Ｓから見たＰＥＴ装置の視野内になるそれぞれのγ線検出器３対を補正対象として設定し、その補正対象として設定されたγ線検出器３を補正対象検出器Ｃとする。
　この補正対象検出器Ｃのうち、基準検出器Ｓに対して約１８０°対向したγ線検出器３を選択する。基準検出器Ｓを「１」番目のγ線検出器３として、時計回りに連番を付すると、図７に示すように全部で１２０個のγ線検出器３の場合には「６１」番目のγ線検出器３が、基準検出器Ｓに１８０°対向したγ線検出器３となる。なお、全部で奇数個のγ線検出器の場合には１８０°対向したγ線検出器３が選択できないので、約１８０°対向したγ線検出器３を選択すればよい。基準検出器Ｓに対して約１８０°対向したγ線検出器３（図７（ａ）の灰色を参照）を対向検出器Ｏとする
　このように同時計数する対象のγ線検出器３を２つ選択する場合に、図７に示す実施態様では、「１」番目のγ線検出器３である基準検出器Ｓと、それに１８０°対向した「６１」番目のγ線検出器３（対向検出器Ｏ）とを選択する。当該選択された「１」番目のγ線検出器３である基準検出器Ｓと、それに１８０°対向した「６１」番目のγ線検出器３（対向検出器Ｏ）とに関する時間差ヒストグラムを基準とする。この基準とした時間差ヒストグラムを、図６（ａ）に示すように“０”を中心とした度数分布とする。したがって、基準とした時間差ヒストグラムにおける総計数値が最大となる時間差を“０”として、この時間差“０”を基準値とする。
　（ステップＳ４）補正対象となる時間差ヒストグラムの抽出
　ステップＳ３で選択された基準検出器Ｓおよび対向検出器Ｏからなる２つのγ線検出器３のうち、一方のγ線検出器３として基準検出器Ｓを選択するとともに、他方のγ線検出器３である対向検出器Ｏとは別の補正対象検出器Ｃを選択する。つまり、対向検出器Ｏ以外の補正対象検出器Ｃと基準検出器Ｓとをそれぞれ選択して、補正対象となる時間差ヒストグラムを抽出する。図７（ａ）に示す場合には、ステップＳ３において「６１」番目のγ線検出器３（対向検出器Ｏ）も含めて補正対象検出器Ｃは「４９」番目～「７３」番目の合計２５個のγ線検出器３であって、「１」番目のγ線検出器３である基準検出器Ｓと、「６１」番目のγ線検出器３（対向検出器Ｏ）以外の「４９」番目～「６０」番目および「６２」番目～「７３」番目の各々の補正対象検出器Ｃとに関する時間差ヒストグラムを、今回で選択された２つのγ線検出器３に関する時間差ヒストグラムとして抽出する。
　（ステップＳ５）補正量の算出
　対向検出器Ｏ以外の補正対象検出器Ｃと基準検出器Ｓとに関する時間差ヒストグラムは、基準とした時間差ヒストグラムに対してズレが生じ、上述したように図６（ｂ）の実線に示すように左右にシフトする。そこで、ステップＳ４において今回で選択された２つのγ線検出器３に関する時間差ヒストグラムを、図６（ｂ）の実線から点線に示す“０”を中心とした度数分布に戻すようなシフト量を補正量としてデータ収集・制御部１１は求める。
　（ステップＳ６）全ての検出器の補正量の算出完了？
　全てのγ線検出器３の補正量の算出が完了したか否かを判定する。完了していない場合には、ステップＳ７に進み、完了した場合には、ステップＳ８に進む。
　（ステップＳ７）補正量の適用
　ステップＳ５で求められた補正量を適用して、今回で選択された２つのγ線検出器３に関する時間差ヒストグラムを補正量の分だけ“０”を中心とした度数分布に戻すようにシフトさせて、当該時間差ヒストグラムをそれぞれ補正する。つまり、基準値である時間差“０”に基づいて当該時間差ヒストグラムをそれぞれ補正する。図７（ａ）に示す場合には、「１」番目のγ線検出器３である基準検出器Ｓと、「４９」番目～「６０」番目および「６２」番目～「７３」番目の各々の補正対象検出器Ｃとに関する時間差ヒストグラムを、過去に選択された「１」番目のγ線検出器３である基準検出器Ｓと、それに１８０°対向した「６１」番目のγ線検出器３（対向検出器Ｏ）とに関する時間差ヒストグラムに基づいて補正する。
　ステップＳ７で補正量を適用して補正を行ったら、ステップＳ３に戻って、ステップＳ７で補正された２つのγ線検出器３を新たに基準とする作業を行う。例えば、図７（ａ）に示す場合には、補正された「４９」番目～「６０」番目および「６２」番目～「７３」番目の各々の補正対象検出器Ｃのうち、補正対象検出器Ｃの両端にある「４９」番目および「７３」番目のγ線検出器３（図７（ａ）の斜線を参照）を新たに基準とする。
　そして、ステップＳ３で、図７（ｂ）に示すように、「４９」番目のγ線検出器３および「７３」番目のγ線検出器３（図７（ｂ）の黒塗りを参照）を基準検出器Ｓとする。図７（ａ）と同様に、各々の基準検出器Ｓから見たＰＥＴ装置の視野内になるそれぞれのγ線検出器３対を補正対象として設定し、その補正対象として設定されたγ線検出器３を補正対象検出器Ｃとする。ここでは、補正対象検出器Ｃは、「４９」番目のγ線検出器３である基準検出器Ｓに対しては「９７」番目～「１２０」番目の合計２４個のγ線検出器３であって、「７３」番目のγ線検出器３である基準検出器Ｓに対しては「２」番目～「２５」番目の合計２４個のγ線検出器３である。
　「４９」番目のγ線検出器３である基準検出器Ｓと、それに１８０°対向した「１０９」番目のγ線検出器３（対向検出器Ｏ）（図７（ｂ）の灰色を参照）とを選択するとともに、「７３」番目のγ線検出器３である基準検出器Ｓと、それに１８０°対向した「１３」番目のγ線検出器３（対向検出器Ｏ）（図７（ｂ）の灰色を参照）とを選択する。
　図７（ａ）と同様に、ステップＳ４で「４９」番目のγ線検出器３である基準検出器Ｓと、「１０９」番目のγ線検出器３（対向検出器Ｏ）以外の「９７」番目～「１０８」番目および「１１０」番目～「１２０」番目の各々の補正対象検出器Ｃとに関する時間差ヒストグラムを、今回で選択された２つのγ線検出器３に関する時間差ヒストグラムとして抽出する。また、ステップＳ４で「７３」番目のγ線検出器３である基準検出器Ｓと、「１３」番目のγ線検出器３（対向検出器Ｏ）以外の「２」番目～「１２」番目および「１４」番目～「２５」番目の各々の補正対象検出器Ｃとに関する時間差ヒストグラムを、今回で選択された２つのγ線検出器３に関する時間差ヒストグラムとして抽出する。
　図７（ａ）と同様に、ステップＳ５で補正量を求め、ステップＳ６で補正量の算出完了の判定を行い、完了していない場合には、ステップＳ７に進む。そして、ステップＳ７で補正量を適用して補正を行い、ステップＳ３に戻って、ステップＳ７で補正された２つのγ線検出器３を新たに基準とする作業を行う。例えば、図７（ｂ）に示す場合には、補正された各々の補正対象検出器Ｃのうち、補正対象検出器Ｃの合計４つの両端（図７（ｂ）の斜線を参照）のうちの下側にある「２５」番目および「９７」番目のγ線検出器３を新たに基準とする。
　そして、ステップＳ３で、図７（ｃ）に示すように、「２５」番目のγ線検出器３および「９７」番目のγ線検出器３（図７（ｃ）の黒塗りを参照）を基準検出器Ｓとする。各々の基準検出器Ｓから見たＰＥＴ装置の視野内になるそれぞれのγ線検出器３対を補正対象として設定し、その補正対象として設定されたγ線検出器３を補正対象検出器Ｃとする。ここでは、補正対象検出器Ｃは、「２５」番目のγ線検出器３である基準検出器Ｓに対しては「７４」番目～「９６」番目の合計２３個のγ線検出器３であって、「９７」番目のγ線検出器３である基準検出器Ｓに対しては「２６」番目～「４８」番目の合計２３個のγ線検出器３である。
　「２５」番目のγ線検出器３である基準検出器Ｓと、それに１８０°対向した「８５」番目のγ線検出器３（対向検出器Ｏ）（図７（ｃ）の灰色を参照）とを選択するとともに、「９７」番目のγ線検出器３である基準検出器Ｓと、それに１８０°対向した「３７」番目のγ線検出器３（対向検出器Ｏ）（図７（ｃ）の灰色を参照）とを選択する。
　ステップＳ４で「２５」番目のγ線検出器３である基準検出器Ｓと、「８５」番目のγ線検出器３（対向検出器Ｏ）以外の「７４」番目～「８４」番目および「８６」番目～「９６」番目の各々の補正対象検出器Ｃとに関する時間差ヒストグラムを、今回で選択された２つのγ線検出器３に関する時間差ヒストグラムとして抽出する。また、ステップＳ４で「９７」番目のγ線検出器３である基準検出器Ｓと、「３７」番目のγ線検出器３（対向検出器Ｏ）以外の「２６」番目～「３６」番目および「３８」番目～「４８」番目の各々の補正対象検出器Ｃとに関する時間差ヒストグラムを、今回で選択された２つのγ線検出器３に関する時間差ヒストグラムとして抽出する。
　ステップＳ５で補正量を求め、ステップＳ６で補正量の算出完了の判定を行い、完了していない場合には、ステップＳ７に進む。そして、ステップＳ７で補正量を適用して補正を行い、ステップＳ３に戻る。この３回目のステップＳ６、Ｓ７、Ｓ３～Ｓ５の繰り返しループにより、「１」番目から「１２０」番目までの全てのγ線検出器３の補正量の算出が完了し、ステップＳ６において完了したと判断する。そして、ステップＳ８に進む。上述したこれらのステップＳ３～Ｓ６は、この発明におけるヒストグラム補正工程に相当する。
　（ステップＳ８）タイミング補正テーブルの設定
　これらのステップＳ３～Ｓ６（図７の実施形態の場合には３回の繰り返しループ）で求められた補正量を、データ収集・制御部１１から検出器信号処理部９のタイミング補正テーブル９ｂに書き込むことで、タイミング補正テーブルを設定する。この設定された補正量を、通常の被検体Ｍを用いた核医学診断での撮像に適用することで、タイミング調整を高精度に行うことが可能となり、上述の撮像を行うときに画質の優れた画像を取得することができる。
　図７の実施形態では、全部で１２０個のγ線検出器３の場合であったが、γ線検出器の数については特に限定されない。例えば、全部で８８個のγ線検出器３のときの基準検出器Ｓ・補正対象検出器Ｏの設定について概略を述べると、１回目のループでは、補正対象検出器Ｃは、「１」番目のγ線検出器３である基準検出器Ｓに対しては「３６」番目～「５４」番目（対向検出器Ｏは「４５」番目）の合計１９個のγ線検出器３である。補正を行った後に補正対象検出器Ｃの両端にある「３６」番目および「５４」番目のγ線検出器３を新たに基準とする。
　２回目のループでは、補正対象検出器Ｃは、「３６」番目のγ線検出器３である基準検出器Ｓに対しては「７１」番目～「８８」番目（対向検出器Ｏは「８０」番目）の合計１８個のγ線検出器３であって、「５４」番目のγ線検出器３である基準検出器Ｓに対しては「２」番目～「１９」番目（対向検出器Ｏは「１０」番目）である。補正を行った後に補正対象検出器Ｃの合計４つの両端のうちの下側にある「１９」番目および「７１」番目のγ線検出器３を新たに基準とする。
　３回目のループでは、補正対象検出器Ｃは、「１９」番目のγ線検出器３である基準検出器Ｓに対しては「５５」番目～「７０」番目（対向検出器Ｏは「６３」番目）の合計１６個のγ線検出器３であって、「７１」番目のγ線検出器３である基準検出器Ｓに対しては「２０」番目～「３５」番目（対向検出器Ｏは「２７」番目）である。
　このように、全部で８８個のγ線検出器３においても、３回目のステップＳ６、Ｓ７、Ｓ３～Ｓ５の繰り返しループにより、「１」番目から「８８」番目までの全てのγ線検出器３の補正量の算出が完了する。
　また図７の実施形態では、ステップＳ６、Ｓ７、Ｓ３～Ｓ５のループ回数は３回であったが、ループ回数についても特に限定されない。例えば、図８の実施形態における基準検出器Ｓ・補正対象検出器Ｏの設定について概略を述べると、１回目のループでは、図８（ａ）に示すように「１」番目のγ線検出器３である基準検出器Ｓと、それに１８０°対向した「６１」番目のγ線検出器３（対向検出器Ｏ）を選択し、対向検出器Ｏに隣接した両隣の「６０」番目および「６２」番目のγ線検出器３を補正対象検出器Ｃとする。
　そして、当該選択された「１」番目のγ線検出器３である基準検出器Ｓと、それに１８０°対向した「６１」番目のγ線検出器３（対向検出器Ｏ）とに関する時間差ヒストグラムを基準とする。当該選択された「１」番目のγ線検出器３である基準検出器Ｓおよび「６１」番目のγ線検出器３（対向検出器Ｏ）からなる２つのγ線検出器３のうち、一方のγ線検出器３として基準検出器Ｓを選択するとともに、他方のγ線検出器３である対向検出器Ｏとは別の上述の「６０」番目および「６２」番目の補正対象検出器Ｃを選択する。
　過去に選択された「１」番目のγ線検出器３である基準検出器Ｓと、それに１８０°対向した「６１」番目のγ線検出器３（対向検出器Ｏ）とに関する時間差ヒストグラムに基づいて、今回で選択された「１」番目のγ線検出器３である基準検出器Ｓと「６０」番目の補正対象検出器Ｃとに関する時間差ヒストグラムを補正するとともに、同じく今回で選択された「１」番目のγ線検出器３である基準検出器Ｓと「６２」番目の補正対象検出器Ｃとに関する時間差ヒストグラムを補正する。補正を行った後に補正対象検出器Ｃの「６０」番目のγ線検出器３および「６２」番目のγ線検出器３を新たに基準とする。
　２回目のループでは、図８（ｂ）に示すように、補正対象検出器Ｃは、「６０」番目のγ線検出器３である基準検出器Ｓに対しては「１２０」番目の対向検出器Ｏに隣接した「１１９」番目のγ線検出器３であって、「６２」番目のγ線検出器３である基準検出器Ｓに対しては「２」番目の対向検出器Ｏの隣接した「３」番目のγ線検出器３である。なお、「１２０」番目の対向検出器Ｏに隣接し、「２」番目の対向検出器Ｏの隣接した「１」番目のγ線検出器３については既に選択されているので除外する。補正を行った後に補正対象検出器Ｃの「３」番目および「１１９」番目のγ線検出器３を新たに基準とする。
　３回目のループでは、図８（ｃ）に示すように、補正対象検出器Ｃは、「３」番目のγ線検出器３である基準検出器Ｓに対しては「６３」番目の対向検出器Ｏに隣接した「６４」番目のγ線検出器３であって、「１１９」番目のγ線検出器３である基準検出器Ｓに対しては「５９」番目の対向検出器Ｏに隣接した「５８」番目のγ線検出器３である。なお、「５９」番目の対向検出器Ｏに隣接した「６０」番目のγ線検出器３、および「６３」番目の対向検出器Ｏに隣接した「６２」番目のγ線検出器３については既に選択されているので除外する。補正を行った後に補正対象検出器Ｃの「５８」番目および「６４」番目のγ線検出器３を新たに基準とする。
　以下、「１」番目から「１２０」番目までの全てのγ線検出器３の補正量の算出が完了するまで、ステップＳ６、Ｓ７、Ｓ３～Ｓ５の繰り返しループを行う。
　なお、図４に示す一連のタイミング補正については、各チャンネルの時間遅れの経年変化に応じて定期的に行うのが好ましい。さらに、補正テーブルをファイル化して、データ収集・制御部１１の記憶媒体あるいは上述のメモリ部８に書き込んで記憶し、ＰＥＴ装置の電源をオフして再度オンにして立ち上げるときに、記憶媒体あるいはメモリ部８から読み出してタイミング補正テーブル９ｂに再度書き込むようにするのが好ましい。
　上述の構成を備えた本実施例１に係るＰＥＴ装置およびタイミング補正方法によれば、放射線（本実施例１ではγ線）を同時計数する各γ線検出器３の対ごとの時間差変化に対する計数値分布を表す時間差ヒストグラムに関して、データ収集・制御部１１は、以下の補正演算をステップＳ３～Ｓ６では行う。すなわち、同時計数する対象のγ線検出器３を２つ選択して（図７（ａ）の場合には、「１」番目の基準検出器Ｓと、それに１８０°対向した「６１」番目の対向検出器Ｏ）、当該選択された２つのγ線検出器３のうち、一方の検出器である基準検出器Ｓを選択するとともに、他方の対向検出器Ｏとは別のγ線検出器３（図７（ａ）の場合には、「４９」番目～「６０」番目および「６２」番目～「７３」番目の補正対象検出器Ｃ）を選択し、その選択を繰り返し行った際に、過去に選択された２つのγ線検出器３（図７（ａ）の場合には、「１」番目の基準検出器Ｓと、それに１８０°対向した「６１」番目の対向検出器Ｏ）に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つのγ線検出器３に関する時間差ヒストグラムを補正する。そして、当該補正された２つのγ線検出器３に関する時間差ヒストグラムを新たに基準とする作業を繰り返し行う（図４のステップＳ６、Ｓ７、Ｓ３～Ｓ５の繰り返しループを参照）。
　このように、過去に選択された２つのγ線検出器３に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つのγ線検出器３に関する時間差ヒストグラムを補正し、当該補正された２つのγ線検出器３に関する時間差ヒストグラムを新たに基準とする作業を繰り返し行うことで、例えば上述した特許文献２のように複数の検出器のタイミング値を平均した場合よりも、最適な時間差ヒストグラムを求める際に収束性がよい。したがって、測定や演算を多数繰り返さずに最適な時間差ヒストグラムを求めることができる。
　また、上述の繰り返し補正された各γ線検出器３の対ごとの時間差ヒストグラムに基づいて、同時計数回路１０は放射線（本実施例１ではγ線）を同時計数するので、正確に同時計数することができる。その結果、測定や演算を多数繰り返さずに、かつ正確に同時計数やタイミング補正を行うことができる。
　本実施例１では、上述の基準とした時間差ヒストグラムにおける総計数値が最大となる時間差を基準値とし、上述の基準値に基づいて時間差ヒストグラムの時間差を補正している。時間差ヒストグラムにおける総計数値が最大となる時間的な箇所は、同時計数が最も起こりうるタイミングである。したがって、そのタイミングである基準値に基づいて時間差ヒストグラムの時間差を補正することで、そのタイミングで揃えることができる。
　本実施例１では、複数のシンチレータ素子（シンチレータブロック３１）からなるシンチレータ素子群（すなわち検出器単位）ごとに補正を行っている。また、本実施例１のγ線検出器３は、各々のシンチレータブロック３１をγ線の深さ方向に積層して構成されたＤＯＩ検出器である。本実施例１のようにこの発明をＤＯＩ検出器に適用することで、例えば上述の特許文献４のように相互作用を起こした深さ方向などの座標情報に対応する検出時刻補正情報が必要な場合に、検出時刻補正情報を求めることができる。すなわち、図４に示す一連のタイミング補正を行った後に、深さ方向の光源位置を重心演算により特定して、深さ方向などの座標情報に対応する検出時刻補正情報を求める。
　本実施例１では、放射性薬剤と同種の放射線を照射する外部線源ＲＩ、または放射性薬剤と同種の放射線を内部から照射するファントムＰｈを備え、上述の時間差ヒストグラムを外部線源ＲＩまたはファントムＰｈからの放射線に基づいて取得している。
　次に、図面を参照してこの発明の実施例２を説明する。図９は、基準シンチレータ素子単位と補正対象シンチレータ素子単位との切り替えの一実施態様を示す正面図である。後述する実施例３、４も含めて、本実施例２に係るＰＥＴ装置は、図１に示すブロック図である。また、上述した実施例１と共通する箇所については同じ符号を付して、その説明を省略する。
　上述した実施例１では、複数のシンチレータ素子（シンチレータブロック３１）からなるシンチレータ素子群（すなわち検出器単位）ごとに補正を行っているが、後述する実施例３も含めて、本実施例２では、下記のように１つのシンチレータ素子からなるシンチレータ素子単位ごとに補正を行って、精度をより向上させている。なお、図９ではシンチレータブロック３１のみを図示し、その他の構造（ライトガイド３２やＰＭＴ３３）については図示を省略する。
　同時計数するγ線検出器の基準シンチレータ素子単位と補正対象シンチレータ素子単位とをそれぞれ設定する。図９（ａ）に示すように、あるγ線検出器の１つのシンチレータブロック３１を基準として設定し、その基準として設定された１つのシンチレータブロック３１を基準シンチレータ素子単位Ｓとする。基準シンチレータ素子単位Ｓとは別のγ線検出器の１つのシンチレータブロック３１を補正対象として設定し、その補正対象として設定された１つのシンチレータブロック３１を補正対象シンチレータ素子単位Ｃとする。具体的なタイミング補正については、実施例１の検出器単位から、実施例２ではシンチレータ素子単位に変更されたのを除けば、同じ手法であるので、その説明を省略する。
　補正後の補正対象シンチレータ素子単位Ｃを新たに基準として、図９（ｂ）に示すように、基準シンチレータ素子単位Ｓとする。基準シンチレータ素子単位Ｓとは別のγ線検出器の１つのシンチレータブロック３１を補正対象として設定し、その補正対象として設定された１つのシンチレータブロック３１を補正対象シンチレータ素子単位Ｃとする。
　そして、補正後の補正対象シンチレータ素子単位Ｃを新たに基準として、図９（ｃ）に示すように、基準シンチレータ素子単位Ｓとする。基準シンチレータ素子単位Ｓとは別のγ線検出器の１つのシンチレータブロック３１を補正対象として設定し、その補正対象として設定された１つのシンチレータブロック３１を補正対象シンチレータ素子単位Ｃとする。このように、他の検出器のシンチレータ素子単位においても、タイミング補正を同様に行う。
　上述の構成を備えた本実施例２に係るＰＥＴ装置およびタイミング補正方法によれば、同時計数する対象のγ線検出器のシンチレータ素子単位を２つ選択して、当該選択された２つのシンチレータ素子単位のうち、一方のγ線検出器のシンチレータ素子単位（図９の場合には基準シンチレータ素子単位Ｓ）を選択するとともに、他方のγ線検出器のシンチレータ素子単位とは別のシンチレータ素子単位（図９の場合には補正対象シンチレータ素子単位Ｃ）を選択し、その選択を繰り返し行った際に、過去に選択された２つのγ線検出器のシンチレータ素子単位に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つのγ線検出器のシンチレータ素子単位に関する時間差ヒストグラムを補正する。そして、当該補正された２つのγ線検出器のシンチレータ素子単位に関する時間差ヒストグラムを新たに基準とする作業を繰り返し行う。このようにシンチレータ素子単位ごとに補正を行うので、上述した実施例１のように検出器単位ごとに補正を行うときと比較して精度をより向上させることができる。
　次に、図面を参照してこの発明の実施例３を説明する。図１０は、基準シンチレータ素子群と補正対象シンチレータ素子単位との切り替えの一実施態様を示す正面図である。後述する実施例４も含めて、本実施例３に係るＰＥＴ装置は、図１に示すブロック図である。また、上述した実施例１、２と共通する箇所については同じ符号を付して、その説明を省略する。
　上述した実施例１では、複数のシンチレータ素子（シンチレータブロック３１）からなるシンチレータ素子群（すなわち検出器単位）ごとに補正を行っているが、上述した実施例２と同様に、本実施例３では、下記のように１つのシンチレータ素子からなるシンチレータ素子単位ごとに補正を行って、精度をより向上させている。なお、実施例２の図９と同様に図１０ではシンチレータブロック３１のみを図示し、その他の構造（ライトガイド３２やＰＭＴ３３）については図示を省略する。
　同時計数するγ線検出器の基準シンチレータ素子群と補正対象シンチレータ素子単位とをそれぞれ設定する。図１０（ａ）に示すように、あるγ線検出器の複数のシンチレータブロック３１を基準として設定し、その基準として設定された複数のシンチレータブロック３１を基準シンチレータ素子群Ｓとする。基準シンチレータ素子群Ｓとは別のγ線検出器の１つのシンチレータブロック３１を補正対象として設定し、その補正対象として設定された１つのシンチレータブロック３１を補正対象シンチレータ素子単位Ｃとする。具体的なタイミング補正については、実施例１の検出器単位から、実施例３ではシンチレータ素子単位に変更されたのを除けば、同じ手法であるので、その説明を省略する。
　補正対象シンチレータ素子単位Ｃが属するγ線検出器において、全てのシンチレータ素子単位の補正量の算出が完了するまでは、図１０（ｂ）に示すように、先に選択された基準シンチレータ素子群Ｓを基準とする。図１０（ａ）に示す補正対象シンチレータ素子単位Ｃとは別のシンチレータブロック３１を補正対象として設定し、その補正対象として設定された１つのシンチレータブロック３１を補正対象シンチレータ素子単位Ｃとする。
　補正対象シンチレータ素子単位Ｃが属するγ線検出器において、全てのシンチレータ素子単位の補正量の算出が完了したら、補正後の補正対象シンチレータ素子単位Ｃの全個数を１つにまとめて新たに基準として、図１０（ｃ）に示すように、基準シンチレータ素子群Ｓとする。基準シンチレータ素子単位Ｓとは別のγ線検出器（図９では、図９（ａ）および図９（ｂ）での基準シンチレータ素子群Ｓのγ線検出器に隣接するγ線検出器）の１つのシンチレータブロック３１を補正対象として設定し、その補正対象として設定された１つのシンチレータブロック３１を補正対象シンチレータ素子単位Ｃとする。このように、他の検出器のシンチレータ素子単位においても、タイミング補正を同様に行う。
　上述の構成を備えた本実施例３に係るＰＥＴ装置およびタイミング補正方法によれば、同時計数する対象のγ線検出器のうち、一方のγ線検出器のシンチレータ素子群（図１０の場合には基準シンチレータ素子群Ｓ）と他方のγ線検出器のシンチレータ素子単位（図１０の場合には補正対象シンチレータ素子単位Ｃ）とを選択して、当該選択されたシンチレータ素子群とシンチレータ素子単位とのうち、当該シンチレータ素子群を選択するとともに、他方のγ線検出器のシンチレータ素子単位とは別のシンチレータ素子単位を選択し、その選択を繰り返し行った際に、過去に選択されたシンチレータ素子群とシンチレータ素子単位とに関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択されたシンチレータ素子群とシンチレータ素子単位とに関する時間差ヒストグラムを補正する。そして、当該補正されたシンチレータ素子群とシンチレータ素子単位とに関する時間差ヒストグラムを新たに基準とする作業を繰り返し行う。このようにシンチレータ素子群およびシンチレータ素子単位ごとに補正を行うので、上述した実施例１のように検出器単位ごとに補正を行うときと比較して精度をより向上させることができる。また、上述した実施例２のようにシンチレータ素子単位ごとに補正を行うときと比較して演算時間や負担を軽減させることができる。
　次に、図面を参照してこの発明の実施例４を説明する。図１１は、自己放射能を持つシンチレータ素子のときの時間差ヒストグラムの説明図である。本実施例４に係るＰＥＴ装置は、図１に示すブロック図である。また、上述した実施例１～３と共通する箇所については同じ符号を付して、その説明を省略する。
　上述した実施例１と相違して、本実施例４では、シンチレータブロックは、Ｌｕ－１７６などに代表される自己放射能（複数の放射線を同時に放出する元素）あるいは自己放射能が添加された物質（例えばＬｕを含むＬＹＳＯ）で構成されている。すなわち、シンチレータブロックは、自己放射能を持つシンチレータ素子であって、検出器は、自己放射能を持つシンチレータ素子を備えている。なお、自己放射能が添加された物質からなる薄膜状のテープを、自己放射能を持たない結晶素子（例えばＧＳＯ）に貼り付けてシンチレータブロックを構成してもよいし、自己放射能が添加された物質からなる塗布剤を、自己放射能を持たない結晶素子に貼り付けてシンチレータブロックを構成してもよい。
　自己放射能を持つシンチレータ素子は、αまたはβ－崩壊を起こし、そのαまたはβ－崩壊に伴ったγ線を放出する核種を含んでいる。図１１では、時間差ヒストグラムを自己放射能からの放射線（ここではγ線）に基づいて取得している。
　上述の実施例１の図６と同様に、図１１に示すように、時間差ヒストグラムは、各検出器３の対ごとにタイムスタンプの差分（すなわち時間差）（図１１では「Difference Time」で表記）を横軸にとり、計数値（図１１では「Event Counts」で表記）を縦軸にとったときの時間差変化に対する計数値分布である。
　自己放射能を持つシンチレータ素子のときには、時間差ヒストグラムは、図１１（ａ）に示すように総計数値が最大となる時間的な箇所および総計数値が二番目に大きい時間的な箇所とが現れる。したがって、本実施例４では、基準とした時間差ヒストグラムにおける総計数値が最大となる時間差と、総計数値が二番目に大きい時間差との中間値にある時間差を基準値とする。この基準値を“０”とする。補正の対象となる時間差ヒストグラムは、基準とした時間差ヒストグラムに対してズレが生じ、図１１（ｂ）の実線に示すように左右にシフトする。そこで、補正の対象となる時間差ヒストグラムを、図１１（ｂ）の実線から点線に示す“０”を基準値に戻すようなシフト量を補正量として求める。そして、この補正量を適用することで、時間差ヒストグラムを補正する。
　上述の構成を備えた本実施例４に係るＰＥＴ装置およびタイミング補正方法によれば、基準とした時間差ヒストグラムにおける総計数値が最大となる時間差と、総計数値が二番目に大きい時間差との中間値にある時間差を基準値とし、上述の基準値に基づいて時間差ヒストグラムの時間差を補正している。
　本実施例４のように、自己放射能を持つシンチレータ素子を検出器が備えた場合に、一方の検出器のシンチレータ素子の自己放射能からの放射線を他方の検出器のシンチレータ素子が最も多く検出したタイミング、および他方の検出器のシンチレータ素子の自己放射能からの放射線を一方の検出器のシンチレータ素子が最も多く検出したタイミングが、総計数値が最大となる時間的な箇所、あるいは総計数値が二番目に大きい時間的な箇所である。このことから、これらのタイミングである２つの時間差の中間値は、同時計数が最も起こりうるタイミングである。したがって、そのタイミングである中間値の基準値に基づいて時間差ヒストグラムの時間差を補正することで、そのタイミングで揃えることができる。
　この発明は、上記実施形態に限られることはなく、下記のように変形実施することができる。
　（１）上述した各実施例では、ポジトロンＣＴ装置（ＰＥＴ装置）単独であったが、ＰＥＴ装置とＣＴ装置とを組み合わせたＰＥＴ－ＣＴ装置にも適用することができる。
　（２）上述した各実施例では、放射線としてγ線を例に採って説明したが、α線やβ線などであってもよい。特に、実施例４のように自己放射能を持つシンチレータ素子を検出器が備えた場合には、αまたはβ－崩壊を起こしたシンチレータ素子の検出器自身でα線またはβ線を検出する場合においても時間差ヒストグラムを用いればよい。
　（３）上述した各実施例では、ＤＯＩ検出器であったが、深さ方向を弁別しない検出器にこの発明は適用することもできる。単数のシンチレータ素子を備えた構造の検出器にこの発明は適用することもできる。
　（４）上述した各実施例では、リング状に配置された検出器に適用したが、リング状に設置されずに単に複数の検出器を備えた場合においても、この発明は適用することができる。

Claims (15)

  1. 　被検体内に投与されたポジトロン放射性薬剤から放出される放射線を検出する複数の検出器を備えたポジトロンＣＴ装置であって、
    　放射線を同時計数する各検出器の対ごとの時間差変化に対する計数値分布を表す時間差ヒストグラムに関して、同時計数する対象の検出器を２つ選択して、当該選択された２つの検出器のうち、一方の検出器を選択するとともに、他方の検出器とは別の検出器を選択し、その選択を繰り返し行った際に、過去に選択された２つの検出器に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つの検出器に関する時間差ヒストグラムを補正し、当該補正された２つの検出器に関する時間差ヒストグラムを新たに基準とする作業を繰り返し行う演算手段と、
    　前記演算手段で繰り返し補正された各検出器の対ごとの前記時間差ヒストグラムに基づいて、放射線を同時計数する同時計数回路と
    　を備えることを特徴とするポジトロンＣＴ装置。
  2. 　請求項１に記載のポジトロンＣＴ装置において、
    　前記基準とした前記時間差ヒストグラムにおける総計数値が最大となる時間差を基準値とし、
    　前記演算手段は、前記基準値に基づいて時間差ヒストグラムの時間差を補正することを特徴とするポジトロンＣＴ装置。
  3. 　請求項１に記載のポジトロンＣＴ装置において、
    　前記基準とした前記時間差ヒストグラムにおける総計数値が最大となる時間差と、総計数値が二番目に大きい時間差との中間値にある時間差を基準値とし、
    　前記演算手段は、前記基準値に基づいて時間差ヒストグラムの時間差を補正することを特徴とするポジトロンＣＴ装置。
  4. 　請求項１から請求項３のいずれかに記載のポジトロンＣＴ装置において、
    　前記検出器は、
    　放射線の入射により蛍光する複数のシンチレータ素子と、
    　各シンチレータ素子からの光を光電変換することで放射線を検出する光電変換手段と
    　を備え、
    　放射線を同時計数する各検出器の１つのシンチレータ素子からなるシンチレータ素子単位の対ごとの前記時間差ヒストグラムに関して、同時計数する対象の検出器のシンチレータ素子単位を２つ選択して、当該選択された２つのシンチレータ素子単位のうち、一方の検出器のシンチレータ素子単位を選択するとともに、他方の検出器のシンチレータ素子単位とは別のシンチレータ素子単位を選択し、その選択を繰り返し行った際に、過去に選択された２つの検出器のシンチレータ素子単位に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つの検出器のシンチレータ素子単位に関する時間差ヒストグラムを補正し、当該補正された２つの検出器のシンチレータ素子単位に関する時間差ヒストグラムを新たに基準とする作業を前記演算手段は繰り返し行い、
    　前記演算手段で繰り返し補正された各検出器のシンチレータ素子単位の対ごとの前記時間差ヒストグラムに基づいて、前記同時計数回路は放射線を同時計数することを特徴とするポジトロンＣＴ装置。
  5. 　請求項１から請求項３のいずれかに記載のポジトロンＣＴ装置において、
    　前記検出器は、
    　放射線の入射により蛍光する複数のシンチレータ素子と、
    　各シンチレータ素子からの光を光電変換することで放射線を検出する光電変換手段と
    　を備え、
    　放射線を同時計数する各検出器のうち、一方の検出器の複数のシンチレータ素子からなるシンチレータ素子群と他方の検出器の１つのシンチレータ素子からなるシンチレータ素子単位との対ごとの前記時間差ヒストグラムに関して、同時計数する対象の検出器のうち、一方の検出器のシンチレータ素子群と他方の検出器のシンチレータ素子単位とを選択して、当該選択されたシンチレータ素子群とシンチレータ素子単位とのうち、当該シンチレータ素子群を選択するとともに、前記他方の検出器のシンチレータ素子単位とは別のシンチレータ素子単位を選択し、その選択を繰り返し行った際に、過去に選択されたシンチレータ素子群とシンチレータ素子単位とに関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択されたシンチレータ素子群とシンチレータ素子単位とに関する時間差ヒストグラムを補正し、当該補正されたシンチレータ素子群とシンチレータ素子単位とに関する時間差ヒストグラムを新たに基準とする作業を前記演算手段は繰り返し行い、
    　前記演算手段で繰り返し補正されたシンチレータ素子群とシンチレータ素子単位との対ごとの前記時間差ヒストグラムに基づいて、前記同時計数回路は放射線を同時計数することを特徴とするポジトロンＣＴ装置。
  6. 　請求項１から請求項５のいずれかに記載のポジトロンＣＴ装置において、
    　前記検出器は、
    　放射線の入射により蛍光する複数のシンチレータ素子と、
    　各シンチレータ素子からの光を光電変換することで放射線を検出する光電変換手段と
    　を備え、
    　前記検出器は、各々の前記シンチレータ素子を放射線の深さ方向に積層して構成されたＤＯＩ検出器であることを特徴とするポジトロンＣＴ装置。
  7. 　請求項１から請求項６のいずれかに記載のポジトロンＣＴ装置において、
    　前記放射性薬剤と同種の放射線を照射する外部線源、または前記放射性薬剤と同種の放射線を内部から照射するファントムを備え、
    　前記時間差ヒストグラムを前記外部線源または前記ファントムからの放射線に基づいて取得することを特徴とするポジトロンＣＴ装置。
  8. 　請求項１から請求項６のいずれかに記載のポジトロンＣＴ装置において、
    　前記検出器は、自己放射能を持つシンチレータ素子を備え、
    　前記時間差ヒストグラムを前記自己放射能からの放射線に基づいて取得することを特徴とするポジトロンＣＴ装置。
  9. 　被検体内に投与されたポジトロン放射性薬剤から放出される放射線を同時計数するために用いられるタイミング補正方法であって、
    　放射線を同時計数する各検出器の対ごとの時間差変化に対する計数値分布を表す時間差ヒストグラムに関して、同時計数する対象の検出器を２つ選択して、当該選択された２つの検出器のうち、一方の検出器を選択するとともに、他方の検出器とは別の検出器を選択し、その選択を繰り返し行った際に、過去に選択された２つの検出器に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つの検出器に関する時間差ヒストグラムを補正し、当該補正された２つの検出器に関する時間差ヒストグラムを新たに基準とする作業を繰り返し行うヒストグラム補正工程
    　を備えることを特徴とするタイミング補正方法。
  10. 　請求項９に記載のタイミング補正方法において、
    　前記基準とした前記時間差ヒストグラムにおける総計数値が最大となる時間差を基準値とし、
    　前記ヒストグラム補正工程では、前記基準値に基づいて時間差ヒストグラムの時間差を補正することを特徴とするタイミング補正方法。
  11. 　請求項９に記載のタイミング補正方法において、
    　前記基準とした前記時間差ヒストグラムにおける総計数値が最大となる時間差と、総計数値が二番目に大きい時間差との中間値にある時間差を基準値とし、
    　前記ヒストグラム補正工程では、前記基準値に基づいて時間差ヒストグラムの時間差を補正することを特徴とするタイミング補正方法。
  12. 　請求項９から請求項１１のいずれかに記載のタイミング補正方法において、
    　放射線を同時計数する各検出器の１つのシンチレータ素子からなるシンチレータ素子単位の対ごとの前記時間差ヒストグラムに関して、同時計数する対象の検出器のシンチレータ素子単位を２つ選択して、当該選択された２つのシンチレータ素子単位のうち、一方の検出器のシンチレータ素子単位を選択するとともに、他方の検出器のシンチレータ素子単位とは別のシンチレータ素子単位を選択し、その選択を繰り返し行った際に、前記ヒストグラム補正工程では、過去に選択された２つの検出器のシンチレータ素子単位に関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択された２つの検出器のシンチレータ素子単位に関する時間差ヒストグラムを補正し、当該補正された２つの検出器のシンチレータ素子単位に関する時間差ヒストグラムを新たに基準とする作業を繰り返し行うことで、タイミング補正を行うことを特徴とするタイミング補正方法。
  13. 　請求項９から請求項１１のいずれかに記載のタイミング補正方法において、
    　放射線を同時計数する各検出器のうち、一方の検出器の複数のシンチレータ素子からなるシンチレータ素子群と他方の検出器の１つのシンチレータ素子からなるシンチレータ素子単位との対ごとの前記時間差ヒストグラムに関して、同時計数する対象の検出器のうち、一方の検出器のシンチレータ素子群と他方の検出器のシンチレータ素子単位とを選択して、当該選択されたシンチレータ素子群とシンチレータ素子単位とのうち、当該シンチレータ素子群を選択するとともに、前記他方の検出器のシンチレータ素子単位とは別のシンチレータ素子単位を選択し、その選択を繰り返し行った際に、前記ヒストグラム補正工程では、過去に選択されたシンチレータ素子群とシンチレータ素子単位とに関する時間差ヒストグラムを基準として、当該基準に基づいて、今回で選択されたシンチレータ素子群とシンチレータ素子単位とに関する時間差ヒストグラムを補正し、当該補正されたシンチレータ素子群とシンチレータ素子単位とに関する時間差ヒストグラムを新たに基準とする作業を繰り返し行うことで、タイミング補正を行うことを特徴とするタイミング補正方法。
  14. 　請求項９から請求項１３のいずれかに記載のタイミング補正方法において、
    　前記時間差ヒストグラムは、前記放射性薬剤と同種の放射線を照射する外部線源、または前記放射性薬剤と同種の放射線を内部から照射するファントムからの放射線に基づいて取得されたデータであることを特徴とするタイミング補正方法。
  15. 　請求項９から請求項１３のいずれかに記載のタイミング補正方法において、
    　前記時間差ヒストグラムは、自己放射能からの放射線に基づいて取得されたデータであることを特徴とするタイミング補正方法。
PCT/JP2010/003340 2010-05-18 2010-05-18 ポジトロンｃｔ装置およびタイミング補正方法 WO2011145139A1 (ja)

Priority Applications (7)

Application Number Priority Date Filing Date Title
EP10851704.6A EP2573588B1 (en) 2010-05-18 2010-05-18 Positron ct apparatus and timing correction method
PCT/JP2010/003340 WO2011145139A1 (ja) 2010-05-18 2010-05-18 ポジトロンｃｔ装置およびタイミング補正方法
JP2012515640A JP5459397B2 (ja) 2010-05-18 2010-05-18 ポジトロンｃｔ装置およびタイミング補正方法
US13/698,971 US9360569B2 (en) 2010-05-18 2010-05-18 Positron CT apparatus and a timing correction method
CN201080066866.2A CN102906595B (zh) 2010-05-18 2010-05-18 正电子ct装置和定时校正方法
TW100117139A TWI444646B (zh) 2010-05-18 2011-05-17 正電子ｃｔ裝置及時序修正方法
US15/142,904 US9844351B2 (en) 2010-05-18 2016-04-29 Positron CT apparatus and a timing correction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2010/003340 WO2011145139A1 (ja) 2010-05-18 2010-05-18 ポジトロンｃｔ装置およびタイミング補正方法

Related Child Applications (2)

Application Number Title Priority Date Filing Date
US13/698,971 A-371-Of-International US9360569B2 (en) 2010-05-18 2010-05-18 Positron CT apparatus and a timing correction method
US15/142,904 Division US9844351B2 (en) 2010-05-18 2016-04-29 Positron CT apparatus and a timing correction method

Publications (1)

Publication Number Publication Date
WO2011145139A1 true WO2011145139A1 (ja) 2011-11-24

Family

ID=44991267

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2010/003340 WO2011145139A1 (ja) 2010-05-18 2010-05-18 ポジトロンｃｔ装置およびタイミング補正方法

Country Status (6)

Country Link
US (2) US9360569B2 (ja)
EP (1) EP2573588B1 (ja)
JP (1) JP5459397B2 (ja)
CN (1) CN102906595B (ja)
TW (1) TWI444646B (ja)
WO (1) WO2011145139A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102648856A (zh) * 2012-04-28 2012-08-29 中国科学院高能物理研究所 正电子发射断层扫描仪及其中的符合判选方法
JP2020190446A (ja) * 2019-05-21 2020-11-26 株式会社島津製作所 放射線検出装置の２次元位置マップの校正方法および放射線検出装置

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5526435B2 (ja) * 2011-08-03 2014-06-18 独立行政法人理化学研究所 Ｐｅｔ装置およびそのイメージング方法
CN107076859B (zh) * 2014-10-17 2019-08-20 皇家飞利浦有限公司 具有光共享和相互作用深度估计的pet探测器闪烁体布置
US10004472B2 (en) * 2014-10-27 2018-06-26 Koninklijke Philips N.V. Pet detector timing calibration
US10527741B2 (en) * 2015-04-07 2020-01-07 Siemens Medical Solutions Usa, Inc. Setup of SIPM based PET detector using LSO background radiation
WO2017114267A1 (zh) * 2015-12-29 2017-07-06 上海联影医疗科技有限公司 一种医疗设备的数据采集系统及其配置方法
WO2019000389A1 (en) 2017-06-30 2019-01-03 Shanghai United Imaging Healthcare Co., Ltd. SYSTEM AND METHOD FOR POSITRON EMISSION TOMOGRAPHY
JP2022045163A (ja) * 2020-09-08 2022-03-18 キヤノンメディカルシステムズ株式会社 核医学診断装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004151105A (ja) * 2002-10-31 2004-05-27 Ge Medical Systems Global Technology Co Llc クリスタルベースの一致タイミング較正法
JP2009540339A (ja) * 2006-06-15 2009-11-19 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 飛行時間式ｐｅｔ用の集積複数チャネル時間−デジタル変換器
JP2009281816A (ja) * 2008-05-21 2009-12-03 Shimadzu Corp 断層撮影装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0619436B2 (ja) 1989-03-31 1994-03-16 株式会社島津製作所 ポジトロンｃｔ装置
US5272344A (en) 1992-11-10 1993-12-21 General Electric Company Automated coincidence timing calibration for a pet scanner
JP3717122B2 (ja) * 2003-09-29 2005-11-16 株式会社日立製作所 γ線の検出時刻決定方法、γ線の同時計数方法、及び核医学診断装置
JP2006090827A (ja) 2004-09-24 2006-04-06 Hitachi Ltd 放射線検査装置及びそのタイミング補正方法
JP2006334085A (ja) * 2005-06-01 2006-12-14 Shimadzu Corp 放射線撮像装置および放射線検出信号処理方法
WO2007043137A1 (ja) * 2005-10-04 2007-04-19 Shimadzu Corporation 核医学診断装置
JP4614002B2 (ja) * 2006-06-08 2011-01-19 株式会社島津製作所 クロマトグラフ質量分析用データ処理装置
JP4877766B2 (ja) 2006-08-25 2012-02-15 独立行政法人放射線医学総合研究所 陽電子放射断層撮像装置及び放射線検出器
JP5126049B2 (ja) * 2007-12-28 2013-01-23 株式会社島津製作所 核医学診断装置、形態断層撮影診断装置、核医学用データ演算処理方法および形態断層画像演算処理方法
US9134441B2 (en) * 2008-03-26 2015-09-15 Shimadzu Corporation Tomographic equipment, imaging system provided therewith, and imaging data acquisition method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004151105A (ja) * 2002-10-31 2004-05-27 Ge Medical Systems Global Technology Co Llc クリスタルベースの一致タイミング較正法
JP2009540339A (ja) * 2006-06-15 2009-11-19 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 飛行時間式ｐｅｔ用の集積複数チャネル時間−デジタル変換器
JP2009281816A (ja) * 2008-05-21 2009-12-03 Shimadzu Corp 断層撮影装置

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102648856A (zh) * 2012-04-28 2012-08-29 中国科学院高能物理研究所 正电子发射断层扫描仪及其中的符合判选方法
JP2020190446A (ja) * 2019-05-21 2020-11-26 株式会社島津製作所 放射線検出装置の２次元位置マップの校正方法および放射線検出装置
JP7247745B2 (ja) 2019-05-21 2023-03-29 株式会社島津製作所 放射線検出装置の２次元位置マップの校正方法および放射線検出装置

Also Published As

Publication number Publication date
JPWO2011145139A1 (ja) 2013-07-22
EP2573588B1 (en) 2017-09-27
CN102906595A (zh) 2013-01-30
EP2573588A4 (en) 2015-10-14
TW201200894A (en) 2012-01-01
CN102906595B (zh) 2014-12-31
EP2573588A1 (en) 2013-03-27
US20130062526A1 (en) 2013-03-14
US9360569B2 (en) 2016-06-07
US9844351B2 (en) 2017-12-19
US20160242706A1 (en) 2016-08-25
TWI444646B (zh) 2014-07-11
JP5459397B2 (ja) 2014-04-02

Similar Documents

Publication Publication Date Title
JP5459397B2 (ja) ポジトロンｃｔ装置およびタイミング補正方法
US9029786B2 (en) Nuclear medicine imaging apparatus, and nuclear medicine imaging method
CN102665564B (zh) 图像诊断装置及位置校正方法
EP2287635B1 (en) Two-dimensional position map correcting method, and radiation detecting apparatus
US8278625B2 (en) Two-dimensional position map correcting method
EP3033637B1 (en) A hybrid method based on simulation and experimental data to normalize pet data
JP2007212242A (ja) 核医学装置
KR102026737B1 (ko) 영상 생성 장치 및 방법
JP4993144B2 (ja) Ｐｅｔ検出器の光量出力を補正する方法
US9140804B2 (en) Methods and systems for determining timing recovery information in a positron emission tomography (PET) system
US20070090298A1 (en) Method To Determine The Depth-Of-Interaction Function For PET Detectors
JP2015017907A (ja) 核医学診断装置および画像処理プログラム
US10895651B2 (en) PET device and method of acquiring gamma ray generation position using scattered coincidence with PET device
JP5310642B2 (ja) 放射線断層撮影装置
JP2015017906A (ja) 核医学診断装置および位置推定方法
JP7247745B2 (ja) 放射線検出装置の２次元位置マップの校正方法および放射線検出装置
JP5794765B2 (ja) 核医学イメージング装置および核医学イメージング装置の作動方法
US20220104781A1 (en) Nuclear medicine diagnostic apparatus
JP5454859B2 (ja) 核医学診断装置及び画像処理装置
JP2009098120A (ja) 撮像装置および光源位置算出方法

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 201080066866.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: 10851704

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2012515640

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 13698971

Country of ref document: US

REEP Request for entry into the european phase

Ref document number: 2010851704

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2010851704

Country of ref document: EP