WO2018159548A1 - ベータ線二次元イメージング装置及び方法 - Google Patents

ベータ線二次元イメージング装置及び方法 Download PDF

Info

Publication number
WO2018159548A1
WO2018159548A1 PCT/JP2018/007009 JP2018007009W WO2018159548A1 WO 2018159548 A1 WO2018159548 A1 WO 2018159548A1 JP 2018007009 W JP2018007009 W JP 2018007009W WO 2018159548 A1 WO2018159548 A1 WO 2018159548A1
Authority
WO
WIPO (PCT)
Prior art keywords
beta
nuclide
ray
gamma
gamma ray
Prior art date
Application number
PCT/JP2018/007009
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 US16/490,276 priority Critical patent/US11259768B2/en
Priority to EP18760568.8A priority patent/EP3591438A4/en
Priority to JP2019502985A priority patent/JP7100902B2/ja
Publication of WO2018159548A1 publication Critical patent/WO2018159548A1/ja

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/29Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
    • G01T1/2914Measurement of spatial distribution of radiation
    • G01T1/2985In 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)
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5235Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/06Devices, other than using radiation, for detecting or locating foreign bodies ; determining position of probes within or on the body of the patient
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/161Applications in the field of nuclear medicine, e.g. in vivo counting
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/037Emission tomography

Definitions

  • the present invention relates to a beta ray two-dimensional imaging apparatus and method.
  • a positron emission tomography (PET) apparatus As one of diagnostic apparatuses, a positron emission tomography (PET) apparatus has been put into practical use.
  • the PET device uses positive beta decay ( ⁇ + decay) that emits positrons (also called positron decay), and when the emitted positrons annihilate with electrons in the material.
  • ⁇ + decay positive beta decay
  • positrons also called positron decay
  • the distribution of the radionuclide that causes the beta decay is estimated. For example, if a drug that accumulates in cancer cells is labeled with a positron emission nucleus, and a living body to which the drug is administered is imaged with a PET device as an imaging target, a three-dimensional distribution image of cancer cells in the living body can be obtained. It is possible to obtain.
  • an imaging device in which imaging elements made of semiconductors are arranged in an array on a plane, is also used.
  • Flat panel detectors are those that directly read out ionization in materials caused by radiation, and have become popular in the medical field.
  • an imaging plate that can read out a two-dimensional distribution of chemical changes depending on the radiation dose with a laser beam is widely used in the field of basic research.
  • the imaging plate it is possible to obtain the two-dimensional information of the imaging object with a simple configuration, and the sensitivity is higher than that of the film, and after reading, the information is erased by irradiation with white light and reused. Because it is possible, it is used to examine the distribution of beta-emitting nuclides in tissue sections and cultured cells. However, since these devices all add up the total amount of radiation energy incident on a certain place during the imaging time, the energy of each incident radiation and the incident time are unknown.
  • radiation imaging detectors have developed a technology that can calculate energy and incident time by processing signals for each event (especially called “photon counting” when reading photons). ing.
  • This technology is applied to CT systems that identify substances based on differences in the X-ray absorption rate of substances due to differences in energy by imaging X-rays having different energies for each detection event. ing.
  • the radionuclide distribution in the object can be obtained in two dimensions by imaging with the imaging plate as the object to be imaged.
  • the conventional imaging plate when a plurality of types of nuclides are distributed in the imaging target, a certain nuclide cannot be distinguished from other nuclides, so that it is possible to image the nuclide in an identifiable manner. Can not.
  • Even in such an imaging plate if a multiple molecule simultaneous imaging method can be applied, it is beneficial because more complex analysis is possible than at present.
  • the energy of beta rays emitted by beta decay is not a specific value depending on the nuclide and is continuous, it is difficult to identify the nuclide even if a photon counting type detector is used. .
  • gamma rays high energy electromagnetic waves
  • this gamma ray is generally called deexcitation gamma ray, it has energy intrinsic to the nuclide, so here it is called “intrinsic gamma ray”, and is distinguished from gamma ray due to annihilation of positron and electron.
  • the beta ray two-dimensional imaging apparatus is in the excited state of the daughter nucleus due to beta decay, and the first intrinsic gamma ray when transitioning to the ground state of the daughter nucleus following the beta ray emission due to the beta decay.
  • the first nuclide that emits and the excited state of the daughter nucleus due to beta decay and has a different energy from the first intrinsic gamma ray when transitioning to the ground state of the daughter nucleus following the beta radiation emission due to the beta decay
  • the beta ray is detected by receiving a beta ray based on the first nuclide or the second nuclide from the imaging target including the second nuclide that emits the second intrinsic gamma ray, and the detection position of the beta ray is set to two.
  • a beta ray detector for outputting a beta ray detection signal including position information indicated by a dimension, and a detector for detecting a gamma ray, the first intrinsic gamma ray and the second intrinsic gamma ray
  • a gamma ray detector for distinguishing and detecting a comma line, a beta ray detection timing by the beta ray detector, a gamma ray detection timing by the gamma ray detector, the position information included in the beta ray detection signal, and the gamma ray detector
  • An imaging processing unit capable of distinguishing and generating the distribution image of the first nuclide and the distribution image of the second nuclide based on whether one of the first eigengamma ray and the second eigengamma ray is detected in , Provided.
  • nuclides that perform positive beta decay can be used, but the beta two-dimensional imaging device employs a method of directly imaging beta rays. Therefore, not only nuclides that perform positive beta decay but also nuclides that perform negative beta decay (beta decay with electron emission) are available and useful.
  • a beta ray two-dimensional imaging apparatus is in an excited state of a daughter nucleus due to beta decay, and has a predetermined energy when transitioning to the ground state of the daughter nucleus following emission of beta rays due to the beta decay.
  • a beta ray detector that detects the beta ray by receiving the signal and outputs a beta ray detection signal including position information indicating the detection position of the beta ray in two dimensions, and a gamma ray that detects the gamma ray including the intrinsic gamma ray A detector, a beta ray detection timing by the beta ray detector, a gamma ray detection timing by the gamma ray detector, and the beta ray detection signal
  • the distribution image of the first nuclide and the other distribution image corresponding to the distribution of the second nuclide are distinguished based on the position information to be detected and whether the intrinsic gamma ray is detected by the gamma ray detector.
  • an imaging processing unit that can be generated.
  • nuclides that perform positive beta decay can be used, but the beta two-dimensional imaging device employs a method of directly imaging beta rays. Therefore, not only nuclides that perform positive beta decay but also nuclides that perform negative beta decay (beta decay with electron emission) are available and useful.
  • a beta ray two-dimensional imaging apparatus includes a first nuclide that emits positrons as beta rays by positive beta decay and a second nuclide that emits electrons as beta rays by negative beta decay.
  • the beta ray is detected by receiving a beta ray based on the first nuclide or the second nuclide from the imaging target, and a beta ray detection signal including position information indicating the detection position of the beta ray in two dimensions is obtained.
  • a beta ray detector for outputting, a gamma ray detector for detecting a gamma ray including a pair annihilation gamma ray generated by a pair annihilation between a positron as a beta ray from the first nuclide and an electron in the beta ray detector; Beta ray detection timing by the beta ray detector, gamma ray detection timing by the gamma ray detector, and the positional information included in the beta ray detection signal. Based on whether the pair annihilation gamma ray is detected by the gamma ray detector, imaging that can be generated by distinguishing the distribution image of the first nuclide and the other distribution image according to the distribution of the second nuclide And a processing unit.
  • nuclides that perform positive beta decay can be used, but the beta two-dimensional imaging device employs a method of directly imaging beta rays. Therefore, not only nuclides that perform positive beta decay but also nuclides that perform negative beta decay (beta decay with electron emission) are available and useful.
  • the beta ray two-dimensional imaging apparatus is an intrinsic gamma ray when the daughter nucleus is excited by beta decay and transitions to the ground state of the daughter nucleus following the beta ray emission by the beta decay.
  • an image that can generate a distribution image of the nuclide based on whether the gamma ray is detected by the gamma ray detector. Characterized by comprising a Jingu processing unit.
  • a nuclide distribution image can be generated based on position information when a gamma ray is detected by a gamma ray detector.
  • the position information when the intrinsic gamma ray is detected by the gamma ray detector is the detection position information of the beta ray. It is limited, and a high-resolution distribution image can be generated by generating a nuclide distribution image based on the detected position information. The reason for the high resolution of the generated distribution image will become more apparent from the following description.
  • the beta two-dimensional imaging method provides a first intrinsic gamma ray when the daughter nucleus is excited by beta decay and transitions to the ground state of the daughter nucleus following the beta radiation emission due to the beta decay.
  • the first nuclide that emits and the excited state of the daughter nucleus due to beta decay and has a different energy from the first intrinsic gamma ray when transitioning to the ground state of the daughter nucleus following the beta radiation emission due to the beta decay
  • the beta ray is detected by receiving a beta ray based on the first nuclide or the second nuclide from the imaging target including the second nuclide that emits the second intrinsic gamma ray, and the detection position of the beta ray is set to two.
  • a daughter nucleus is excited by beta decay, and has a predetermined energy when transitioning to the ground state of the daughter nucleus following beta emission by the beta decay.
  • Beta rays based on the first nuclide or the second nuclide from an imaging target including a first nuclide that emits intrinsic gamma rays and a second nuclide that transitions to the ground state of the daughter nucleus due to the release of beta rays due to beta decay.
  • a beta ray detection step for detecting the beta ray by receiving the light and obtaining a beta ray detection signal including position information indicating the detection position of the beta ray in two dimensions; and a gamma ray detection for detecting the gamma ray including the intrinsic gamma ray
  • a beta ray detection timing in the beta ray detection step and a gamma ray detection timing in the gamma ray detection step Others corresponding to the distribution image of the first nuclide and the distribution of the second nuclide based on the position information included in the beta ray detection signal and whether the intrinsic gamma ray was detected in the gamma ray detection step.
  • an imaging processing step capable of generating the distribution images separately.
  • a beta ray two-dimensional imaging method includes a first nuclide that emits positrons as beta rays by positive beta decay and a second nuclide that emits electrons as beta rays by negative beta decay.
  • the beta ray is detected by receiving a beta ray based on the first nuclide or the second nuclide from the imaging target, and a beta ray detection signal including position information indicating the detection position of the beta ray in two dimensions is obtained.
  • a beta ray two-dimensional imaging method provides an intrinsic gamma ray when a daughter nucleus is excited by beta decay and transitions to the ground state of the daughter nucleus following beta emission by the beta decay.
  • Detecting the beta ray by receiving the beta ray from the imaging target including the nuclide that emits the beta ray, and obtaining the beta ray detection signal including the position information indicating the detection position of the beta ray in two dimensions
  • a gamma ray detection step for detecting a gamma ray including the intrinsic gamma ray, a beta ray detection timing in the beta ray detection step, a gamma ray detection timing in the gamma ray detection step, and the position information included in the beta ray detection signal Based on whether the intrinsic gamma ray was detected in the gamma ray detection step, Serial and Imaging Step capable of generating a distribution image of nuclides, characterized by comprising a.
  • the present invention it is possible to provide a beta ray two-dimensional imaging apparatus and method for realizing imaging in a mode in which nuclides can be identified. Moreover, according to this invention, it becomes possible to provide the beta ray two-dimensional imaging apparatus and method which contribute to the resolution improvement in the imaging of a nuclide.
  • FIG. 6 is a diagram showing a configuration example of a beta ray detector in the imaging apparatus of FIG. 5.
  • (A) And (b) is a figure which shows the example of a structure of the scintillator which concerns on embodiment of this invention, respectively, and is a figure which shows the relationship of three positions in connection with a detection of a beta ray.
  • (A) And (b) is a figure which shows the mode of the signal when B type radioactive decay generate
  • FIG. 1 shows the type of A-type radioactive decay
  • FIG. 2 shows the type of B-type radioactive decay.
  • Radionuclides that undergo radioactive decay in the A-type and B-type modes are referred to as A-type and B-type nuclides, respectively.
  • the nuclide before beta decay is called a parent nucleus
  • the nuclide after beta decay is called a daughter nucleus.
  • Beta decay in A-type nuclides may be negative beta decay with the emission of electrons (negatively charged beta particles) or with the release of positrons (positively charged beta particles). It may be a positive beta decay. That is, due to beta decay, beta rays due to electrons or positrons are emitted from the A-type nuclide.
  • the gamma ray having energy specific to the B-type nuclide (hereinafter referred to as gamma decay) (Referred to as intrinsic gamma rays), the transition from the excited state of the daughter nucleus to the ground state of the daughter nucleus.
  • the energy difference between the excited state and the ground state in the daughter nucleus of the B type nuclide becomes the energy of the intrinsic gamma ray.
  • the timing at which intrinsic gamma rays are emitted follows a quantum mechanical probability that depends on the structure of the nucleus.
  • the lifetime index from the daughter nucleus excited state to the daughter nucleus ground state is given as the half-life.
  • the beta decay in the B type nuclide may be a negative beta decay accompanied by the emission of electrons (beta particles having a negative charge) as in the beta decay in the A type nuclide, or may be a positron (positive charge). It may be a positive beta decay accompanied by the release of beta particles). That is, beta rays by electrons or positrons are emitted from the B type nuclide by beta decay.
  • FIG. 3 shows an example of A-type radioactive decay.
  • beta decay occurs in 18 F which is the parent nucleus of the A-type nuclide, it transitions to the ground state of 18 O which is the daughter nucleus with a probability of 100%.
  • the beta decay at this time is positive beta decay, and a positron is emitted by the beta decay.
  • FIG. 4 shows an example of B-type radioactive decay.
  • 22 Na which is the parent nucleus of the B-type nuclide, passes through the excited state of 22 Ne with a probability of about 99.96% before transitioning to the ground state of 22 Ne, which is the daughter nucleus.
  • the gamma decay transitions to the ground state about 3.63fs 22 Ne 22 Ne from the excited state of the at half-life (femtoseconds), when the transition energy 1275keV of (keV) Emits intrinsic gamma rays.
  • the beta decay that causes the transition from 22 Na to the excited state of 22 Ne is a positive beta decay, and a positron is emitted by the beta decay.
  • a plurality of excited states in the daughter nucleus may be passed through in the process from the parent nucleus to the transition to the ground state of the daughter nucleus.
  • FIG. 5 shows a schematic configuration diagram of an imaging apparatus (beta ray two-dimensional imaging apparatus) 1 according to the present embodiment.
  • the imaging apparatus 1 includes a beta ray detector 10 that is a position-resolved beta ray detector, a gamma ray detector 20 that is an energy-resolved gamma ray detector, and an imaging processing unit 30.
  • TG represents an imaging target in the imaging apparatus 1.
  • the X axis, the Y axis, and the Z axis are orthogonal to each other, and a three-dimensional space defined by the X axis, the Y axis, and the Z axis is referred to as an XYZ space.
  • a two-dimensional plane parallel to the X axis and the Y axis is referred to as an XY plane.
  • the X-axis component, Y-axis component, and Z-axis component at a certain position are represented by x, y, and z, respectively, and the position of a certain point in the XYZ space is represented by (x, y, z).
  • a position on the XY plane focusing only on the X-axis component and the Y-axis component is represented by (x, y). Furthermore, for convenience of explanation, the directions of the Z axis toward the positive and negative directions are considered as the upward direction and the downward direction, respectively.
  • the imaging target TG is a two-dimensional object that spreads in the X-axis and Y-axis directions, or an object that can be regarded as a two-dimensional imaging target.
  • the imaging target TG includes two or more types of B type nuclides, or includes one or more types of A type nuclides and one or more types of B type nuclides, and these nuclides are included in the imaging target TG.
  • the drug can be labeled with an A type nuclide or a B type nuclide, and each of the A type nuclide and the B type nuclide may be included in the imaging target TG in a form incorporated in the drug.
  • An agent labeled with an A type nuclide ie, an agent incorporating an A type nuclide
  • an agent labeled with a B type nuclide ie, an agent incorporated with a B type nuclide
  • type B an agent labeled with a B type nuclide
  • Beta-ray detector ⁇ The beta ray detector 10 is disposed adjacent to the imaging target TG on the lower side of the imaging target TG so as to receive beta rays emitted from the radionuclide in the imaging target TG.
  • the beta ray detector 10 detects the beta ray by receiving the beta ray emitted from the radionuclide in the imaging target TG, and outputs a beta ray detection signal every time it detects the incidence of the beta ray on itself.
  • the beta ray detection signal indicates that a beta ray is incident on the beta ray detector 10 and includes position information indicating a beta ray detection position in the beta ray detector 10.
  • the detection position of the beta ray here is a position on the XY plane. That is, the position information in the beta ray detection signal includes only the X axis component and the Y axis component of the beta ray detection position.
  • the imaging target TG is a two-dimensional imaging target and has a low beta ray transmission capability
  • a beta ray is detected at a position (x 1 , y 1 , z 2 ) in the beta ray detector 10.
  • the existence position of the radionuclide that has emitted the beta ray is most likely a position (x 1 , y 1 , z 1 ) in the imaging target TG (see FIG. 6). Therefore, when attention is paid only to the two-dimensional component composed of the X-axis and Y-axis components, the detection position of the beta ray can be regarded as the position of the radionuclide that emits the beta ray.
  • the position information included in the beta-ray detection signal represents the X-axis component and the Y-axis component of the position of the radionuclide that emits beta rays in the imaging target TG.
  • a beta ray detector having such a function is known, and a known position-resolved beta ray detector can be adopted as the beta ray detector 10.
  • the beta ray detector 10 may also detect gamma rays. However, since gamma rays have high penetrating power, they are thick enough to reliably stop beta rays from nuclides used for imaging (minimum required minimum). By preparing a beta ray detector 10 having a thickness (generally within a few millimeters), most gamma rays penetrate the beta ray detector 10.
  • the beta ray detector 10 is composed of a scintillation detector including a scintillator 11 and a photodetection amplifier 12. Any known scintillation detector can be used as the beta ray detector 10, but the configurations and functions of the scintillator 11 and the photodetection amplifier 12 will be briefly described below.
  • the scintillator 11 is arranged adjacent to the imaging target TG on the lower side of the imaging target TG so as to receive beta rays emitted from the radionuclide in the imaging target TG.
  • the scintillator 11 is a plate-like scintillator material that spreads in the X-axis and Y-axis directions, and emits scintillation light by the incidence of beta rays. As long as this function is realizable, the scintillator material which comprises the scintillator 11 is arbitrary.
  • the photodetection amplifier 12 is disposed adjacent to the scintillator 11 on the lower side of the scintillator 11 and optically joined so as to receive the scintillation light from the scintillator 11. Therefore, the scintillator 11 is inserted between the imaging target TG and the photodetection amplifier 12.
  • the photodetection amplifier 12 amplifies the scintillation light from the scintillator 11, converts it into an electrical signal, and outputs it. More specifically, when the scintillation light is incident on the photodetection amplifier 12, information indicating the X-axis component and the Y-axis component of the incident position (detection position) of the scintillation light is used as the position information (x, y).
  • the incident position of the scintillation light in the photodetection amplifier 12 is the detection position of the beta ray in the beta ray detector 10 (in other words, the generation position of the scintillation light in the scintillator 11).
  • the position of the radionuclide that emits beta rays in the imaging target TG is assumed that an MPPC (registered trademark) configured by arranging a high gain avalanche photodiode (hereinafter referred to as APD) operating in Geiger mode in an array is used as the photodetection amplifier 12.
  • APD high gain avalanche photodiode
  • a scintillator 11a (divided) is formed by dividing a plate-like scintillator material along each of the X-axis and Y-axis directions so as to increase the resolution of position detection in the beta ray detector 10. May be used as the scintillator 11.
  • a reflecting material may be provided on the end face orthogonal to the XY plane, or the end face may be fogged (uneven unevenness is provided). This is for reflecting or diffusing the scintillation light reaching the end face.
  • the scintillator material is incident on the scintillator material at the position (x 1 , y 1 , z 2 ), and the scintillation light generated by the scintillator material is directed to the end surface of the scintillator material (the end surface orthogonal to the XY plane).
  • the scintillation light from the scintillator material is incident on the position directly below the position (x 1, y 1, z 2) (x 1, y 1, z 3)
  • the probability of being converted into an electric signal at the position (x 1 , y 1 , z 3 ) is increased (see FIG. 8B). That is, the position (x 1 , y 1 ) on the XY plane of the radionuclide that emits beta rays can be accurately detected.
  • a single and plate-like scintillator material may be used as the scintillator 11, and slits parallel to the X-axis and Y-axis directions may be provided for the scintillator material.
  • the scintillator 11a is used.
  • a reflecting material may be provided on the end surface so that the scintillation light reaching the end surface orthogonal to the XY plane among the end surfaces of each slit is reflected or diffused. It is also possible to provide a rough surface.
  • the optical detection amplifier 12 is formed using MPPC (registered trademark)
  • the optical detection amplifier 12 may be formed using an APD that operates in a normal mode different from the Geiger mode.
  • the photodetection amplifier 12 may be composed of a photomultiplier tube.
  • the beta ray detector 10 may be configured by a semiconductor detector that directly detects beta rays without using the scintillator 11.
  • the gamma ray detector 20 is arranged on the upper side of the imaging target TG and separated from the imaging target TG so as to receive gamma rays emitted from the radionuclide in the imaging target TG.
  • the gamma ray detector 20 detects the gamma ray by receiving the gamma ray (mainly the intrinsic gamma ray from the radionuclide in the imaging target TG), and the gamma ray (mainly the intrinsic gamma ray from the radionuclide in the imaging target TG) enters the gamma ray detector 20.
  • a gamma ray detection signal is output every time.
  • the gamma ray detection signal indicates that the gamma ray was incident on the gamma ray detector 20 and includes energy information E indicating the energy of the incident gamma ray.
  • energy information E indicating the energy of the incident gamma ray.
  • a pulsed voltage signal having a peak value corresponding to the energy of the gamma rays is output as a gamma ray detection signal.
  • the gamma ray detector 20 can be configured by a semiconductor detector made of Ge (germanium). In this case, the gamma ray detector 20 can detect gamma ray energy with an energy resolution of about 0.2 to 0.5%.
  • a semiconductor detector configured using a semiconductor material other than Ge (for example, Si, CdTe, CdZnTe) may be used as the gamma ray detector 20, or the gamma ray detector 20 is configured using a scintillation detector. May be.
  • scintillation detectors are inferior in energy resolution compared to semiconductor detectors, but it is sufficient that they have an energy resolution capable of discriminating the target intrinsic gamma rays.
  • the imaging processing unit 30 generates a distribution image of the target nuclide included in the imaging target TG based on the beta ray detection signal output from the beta ray detector 10 and the gamma ray detection signal output from the gamma ray detector 20.
  • the target nuclide refers to any nuclide that is included in the imaging target TG and is the target of the distribution imaging, for example, may be any specific B type nuclide included in the imaging target TG, When the imaging target TG includes an A type nuclide and a B type nuclide, it may be a combination of the A type nuclide and the B type nuclide.
  • the distribution image of the target nuclide refers to an image showing the distribution of the target nuclide in a plane parallel to the XY plane.
  • an image obtained by projecting the position of the target nuclide in the imaging target TG onto a predetermined projection plane parallel to the XY plane corresponds to the distribution image of the target nuclide.
  • the distribution image of the target nuclide is also the distribution image of the A type probe.
  • the distribution image of the target nuclide is also a distribution image of the B type probe.
  • the imaging processing unit 30 includes a simultaneous measurement determination unit 31 and an image reconstruction unit 32 as shown in FIG.
  • the simultaneous measurement discriminating unit 31 receives the output of the beta ray detector 10 and the output of the gamma ray detector 20.
  • the simultaneous measurement discriminating unit 31 discriminates whether the beta ray detection by the beta ray detector 10 and the gamma ray detection by the gamma ray detector 20 occur at the same time (in other words, whether the beta ray and the gamma ray are measured simultaneously), A simultaneous measurement discrimination signal indicating the discrimination result is output.
  • “simultaneous” is a concept having a predetermined time width, and if the time difference between the beta ray detection timing by the beta ray detector 10 and the gamma ray detection timing by the gamma ray detector 20 is equal to or less than a predetermined value, those It is determined that the detection timings are simultaneous, otherwise, it is determined that the detection timings are not simultaneous. Such determination is hereinafter referred to as simultaneous measurement determination.
  • a timing t beta 1 to beta rays from A type nuclide or B type species within the imaging target TG is released, the timing t beta 2 in which the beta ray is caused to interact with the detector 10 at a beta ray detector 10, beta a timing t beta 3 beta ray detection signal is output from the line detector 10 indicating the detection result, the time difference finite between the timing t beta 4 in which the beta ray detection signal is input to the simultaneous measurement determination section 31 Yes (with the passage of time, timings t ⁇ 1 , t ⁇ 2 , t ⁇ 3 , t ⁇ 4 come in this order).
  • the time from timing t ⁇ 1 to timing t ⁇ 4 is determined by the response of the detector 10, the length of the signal line from the detector 10, and the like, and has a length inherent to the beta ray detector system including the detector 10. Also, due to the difference in time response due to the difference in detection position and the influence of signal noise, from the timing t ⁇ 1 when the beta ray is actually emitted to the timing t ⁇ 4 when the beta ray detection signal is input to the simultaneous measurement determination unit 31. The time difference has a time spread. The spread of the time difference indicates how accurately the radiation emission time can be determined, and is generally called time resolution.
  • timing Tiganma 1 to B type nuclide than gamma rays in the imaging target TG is released
  • the gamma at the gamma ray detector 20 is a timing Tiganma 2 that caused the interaction with the detector 20
  • a timing Tiganma 3 of gamma ray detection signal indicating the detection result from the gamma ray detector 20 is outputted
  • there is a time difference of finite between the timing Tiganma 4 where the gamma ray detection signal is input to the simultaneous measurement determination section 31 (As time passes, the timings t ⁇ 1 , t ⁇ 2 , t ⁇ 3 , t ⁇ 4 come in this order).
  • the time difference between the timings t ⁇ 1 and t ⁇ 4 has a length and a time resolution inherent to the gamma detector system that includes the detector 20.
  • the simultaneous measurement discriminating unit 31 detects the difference in the input timing of the beta ray detection signal and the gamma ray detection signal relative to itself, and detects the beta ray detection timing and the gamma ray detection by the beta ray detector 10. It can be considered that this is a time difference from the gamma ray detection timing by the device 20.
  • the timings t ⁇ 4 and t ⁇ 4 are respectively expressed as the beta ray detection timing by the beta ray detector 10 (beta ray detector system) and the gamma ray by the gamma ray detector 20 (gamma ray detector system). It can also be understood as the detection timing.
  • the simultaneous measurement discriminating unit 31 detects the gamma ray from the gamma ray detector 20 after the input timing and within a predetermined time TTH from the input timing.
  • the beta ray and the gamma ray are measured at the same time and a simultaneous measurement discrimination signal having a logical value of “1” is output. Otherwise, beta rays and gamma rays (including intrinsic gamma rays) are not measured simultaneously.
  • the simultaneous measurement determination signal having the logical value of “0” is output by determination.
  • the time T TH used for simultaneous measurement discrimination is determined mainly depending on the respective time resolutions of the beta ray detector 10 and the gamma ray detection 20.
  • the time spread of the beta ray detection time determination accuracy by the beta ray detector 10 and the time spread of the gamma ray detection time determination accuracy by the gamma ray detector 20 are taken into consideration, so that these time spreads are included.
  • the length of time T TH is determined.
  • the image processing unit 30 creates event data based on the output of the beta ray detector 10 and the output of the gamma ray detector 20 for each event.
  • the event data may be recorded in a memory included in the imaging apparatus 1 (for example, a memory included in the image processing unit 30; not shown). It is efficient to set the conditions (master trigger) for creating and recording the event data only depending on whether or not radiation (beta rays) is detected in the beta ray detector 10. In this case, if radiation (gamma rays) is detected only in the gamma ray detector 20, event data is not recorded. However, imaging cannot be performed only with information from the gamma ray detector 20. It is efficient because it does not take in data.
  • the beta ray detector 10 and the gamma ray detector 20 are operated independently, and a common time stamp is added to the data (data forming the event data) from each detector, and simultaneous measurement discrimination is performed after measurement.
  • radiation detection by each detector serves as a master trigger for recording event data.
  • One event data for one event includes position information (x, y) in the beta ray detection signal for the event and a simultaneous measurement determination signal for the event, and further a gamma ray detection signal for the event. May contain energy information E therein.
  • the position information (x, y) included in the event data indicates the detection position of the beta ray at the beta ray detector 10 in the event, and the beta ray at the event.
  • the position of the nuclide that released are positions in the X-axis direction and the Y-axis direction as described above.
  • energy information E that can be included in the event data indicates the energy of gamma rays detected for the event.
  • FIG. 9 is an explanatory diagram of an event due to beta decay of an A type nuclide (hereinafter referred to as an A type event).
  • a type event when the A type nuclide emits a beta ray, a beta ray detection signal including position information (x, y) indicating the position of the A type nuclide is generated, but an intrinsic gamma ray is not emitted. Therefore, the event data in the A type event includes the position information (x, y) and the simultaneous measurement determination information of “0”, and does not include the energy information E.
  • beta rays due to beta decay of the A type nuclide may not be detected by the beta ray detector 10, but such an event does not form an event (the same applies to the B type nuclide).
  • FIG. 10 (a) and 10 (b) are explanatory diagrams of an event due to beta decay of a B type nuclide (hereinafter referred to as a B type event).
  • a B-type nuclide emits a beta ray to generate a beta-ray detection signal including position information (x, y) indicating the position of the B-type nuclide, while the B-type nuclide A gamma ray detection signal can be generated by emitting a specific gamma ray following the beta ray.
  • the event data is a beta ray detection signal.
  • the event data includes the position information (x, y). It consists of simultaneous measurement discrimination information of “0” and does not include energy information E.
  • the solid angle of detection of the specific gamma ray in the gamma ray detector 20 is increased, the detection probability of the specific gamma ray is increased in the B type event. Therefore, the gamma ray detector 20 is only made to the imaging target TG so that the solid angle is increased. It is advisable to use a gamma ray detector 20 that is arranged close to and has a sufficient size to cover the imaging target TG.
  • the imaging apparatus 1 can easily realize a large solid angle close to 2 ⁇ even when compared to the solid angle of detection of intrinsic gamma rays in a PET apparatus as disclosed in Patent Document 1 (Japanese Patent No. 5526435).
  • the solid angle is formed by configuring the gamma ray detector 20 with a gamma ray detector 20A disposed above the imaging target TG and a gamma ray detector 20B disposed below the imaging target TG. Can be further increased. In this case, the imaging target TG and the beta ray detector 10 are sandwiched between the gamma ray detector 20A and the gamma ray detector 20B. In the example of FIG.
  • the gamma ray detector 20A and the gamma ray detector 20B constitute one gamma ray detector 20, and the gamma ray detector 20 has a logical sum of the output of the gamma ray detector 20A and the output of the gamma ray detector 20B.
  • An output signal is formed. That is, the gamma ray detector 20 in FIG. 11 outputs a gamma ray detection signal when a gamma ray (including a specific gamma ray) is detected at any one of the existence position of the gamma ray detector 20A and the existence position of the gamma ray detector 20B. .
  • the image reconstruction unit 32 generates a distribution image of the target nuclide by performing image reconstruction processing based on a plurality of event data.
  • the intensity distribution of radioactivity is obtained by histogramming the (x, y) beta ray detection position information included in the obtained event data on a pixelated XY plane.
  • the configuration and operation of the imaging apparatus 1 described above will be referred to as a basic embodiment.
  • a specific example of a generated distribution image, a further detailed configuration example, an operation example, and an application example of the imaging apparatus 1 will be described among a plurality of embodiments.
  • the matters described in the basic examples apply to each example described later, and in each example, the description in each example is given priority for items that contradict the basic example. .
  • the matters described in any of the plurality of embodiments described below can be applied to any other embodiment (that is, any two or more of the plurality of embodiments). It is also possible to combine these embodiments).
  • the imaging target TG includes a B type nuclide NC B1 and a B type nuclide NC B2 which are different B type nuclides.
  • the imaging target TG may further include an A type nuclide, but here it is assumed that no A type nuclide is included.
  • the energies of the intrinsic gamma rays of the B type nuclides NC B1 and NC B2 are E 1 and E 2 , respectively, and the energy E 1 and E 2 are different from each other.
  • Gamma ray detector 20 can detect the gamma rays to decompose the energy E 1 and energy E 2. That is, the gamma ray detector 20 can distinguish and detect the intrinsic gamma ray with energy E 1 and the intrinsic gamma ray with energy E 2 .
  • the energy information E indicating that the energy of the detected specific gamma rays the energy E 1 (hereinafter, energy information E 1 Energy indicating that the energy of the detected intrinsic gamma ray is energy E 2 when the intrinsic gamma ray emitted from the B type nuclide NC B2 is detected.
  • a gamma ray detection signal including information E (hereinafter also referred to as energy information E 2 ) is output.
  • the output of the gamma ray detection signal including the energy information E 1 indicates that the intrinsic gamma ray from the B type nuclide NC B1 has been detected by the gamma ray detector 20, and the gamma ray detection signal including the energy information E 2 The output indicates that the gamma ray from the B type nuclide NC B2 is detected by the gamma ray detector 20.
  • FIG. 12 is a conceptual diagram of image reconstruction according to the first embodiment.
  • the image 200 shows the actual distribution of the B type nuclides NC B1 and NC B2 in the imaging target TG.
  • the event data group 210 is a collection of event data for events that do not involve detection of specific gamma ray event data group is a collection of event data for events that involve the detection of specific gamma ray energy E 1 211 and an event data group 212 which is a collection of event data for events involving detection of the intrinsic gamma ray of energy E 2 .
  • the event data regarding the event accompanied by the detection of the intrinsic gamma ray of energy E 1 is event data including a simultaneous measurement determination signal of “1”, and includes position information (x, y) and energy information E 1 .
  • the event data regarding the event accompanied by the detection of the intrinsic gamma ray of energy E 2 is event data including a simultaneous measurement determination signal of “1”, and includes position information (x, y) and energy information E 2 .
  • the event data for an event that does not involve detection of the inherent gamma ray is event data for an event in which the inherent gamma ray is not detected by the gamma ray detector 20, and typically includes an event including a simultaneous measurement determination signal of “0”. Therefore, it includes position information (x, y) but does not include energy information E.
  • position information x, y
  • the simultaneous measurement determination signal in the event data may indicate “1”. If, however, the energy information E in the event data does not indicate the energy E 1 or E 2 , the event data is classified as event data for an event that does not involve detection of intrinsic gamma rays.
  • coincidence coincidence is an event that always occurs in an apparatus that performs simultaneous measurement, although the frequency changes depending on the number of detections per unit time.
  • the coincidence coincidence in the imaging apparatus 1 can be corrected by a technique equivalent to the background removal by the known delayed simultaneous measurement used in PET.
  • the position information (x, y) of each event data in the event data group 211 represents the detected value of the X-axis component and the Y-axis component of the existing position of the B type nuclide NC B1 in the imaging target TG. For this reason, the image reconstruction unit 32 histograms the position information (x, y) of each event data in the event data group 211 on the XY plane, thereby distributing the B type nuclide NC B1 in the imaging target TG. An image 221 can be generated.
  • the imaging processing unit 30 determines that the time difference between the beta ray detection timing by the beta ray detector 10 and the gamma ray detection timing by the gamma ray detector 20 is equal to or less than the predetermined time T TH (that is, the simultaneous measurement determination signal is “1”).
  • the distribution image 221 of the B type nuclide NC B1 can be generated using the position information (x, y) when the energy information E 1 is obtained.
  • the imaging processing unit 30 in the event data indicating that the intrinsic gamma ray of energy E 1 is detected by the gamma ray detector 20 within the predetermined time T TH from the beta ray detection timing by the beta ray detector 10.
  • the distribution image 221 of the B-type nuclide NC B1 can be generated using the position information (x, y).
  • the position information (x, y) of each event data in the event data group 212 represents the detected values of the X-axis component and the Y-axis component of the existence position of the B type nuclide NC B2 in the imaging target TG. Therefore, the image reconstruction unit 32 histograms the position information (x, y) of each event data in the event data group 212 on the XY plane, thereby distributing the B type nuclide NC B2 in the imaging target TG. An image 222 can be generated.
  • the imaging processing unit 30 determines that the time difference between the beta ray detection timing by the beta ray detector 10 and the gamma ray detection timing by the gamma ray detector 20 is equal to or less than the predetermined time T TH (that is, the simultaneous measurement determination signal is “1”).
  • the distribution image 222 of the B type nuclide NC B2 can be generated using the position information (x, y) when the energy information E 2 is obtained.
  • the imaging processing unit 30 in the event data indicating that the intrinsic gamma ray of energy E 2 is detected by the gamma ray detector 20 within the predetermined time T TH from the beta ray detection timing by the beta ray detector 10.
  • the distribution image 222 of the B-type nuclide NC B2 can be generated using the position information (x, y).
  • the position information (x, y) of each event data in the event data group 210 includes the detected values of the X-axis component and the Y-axis component at the position of any of the B type nuclides NC B1 and NC B2 in the imaging target TG. Represents. For this reason, the image reconstruction unit 32 histograms the position information (x, y) of each event data in the event data group 210 on the XY plane, so that all the nuclides included in the imaging target TG (here, A distribution image 220 corresponding to a superposition of the distributions of the B type nuclides NC B1 and NC B2 ) can be generated. That is, the imaging processing unit 30 can generate the distribution image 220 by using the position information (x, y) in the event data for the event that does not involve the detection of the intrinsic gamma ray.
  • the imaging apparatus 1 it is possible to image the distributions of a plurality of types of B-type nuclides separately from each other at a time.
  • the imaging apparatus 1 has a function of generating the distribution images 220, 221 and 222, but generates only one of the distribution images (particularly, for example, the distribution image 221 or 222). It may operate, or may operate so as to generate only two distribution images (particularly, for example, distribution images 221 and 222).
  • pair annihilation occurs between the positrons and nearby electrons after the positrons contribute to the generation of scintillation light by the scintillator 11.
  • a gamma ray having energy of 511 keV generated by pair annihilation may enter the gamma ray detector 20.
  • nuclides that emit positrons due to positive beta decay when the positrons enter the beta-ray detector 10, they finally annihilate a pair, and energy of 511 keV is applied to each of them in directions opposite to each other by 180 °.
  • Two gamma rays (pair annihilation gamma rays) are generated. This pair annihilation gamma ray may be detected by the gamma ray detector 20.
  • characteristic X-rays may be generated, and the characteristic X-rays may enter the gamma ray detector 20.
  • the energy of characteristic X-rays (several keV to hundreds of tens of keV) is lower than the energy of intrinsic gamma rays (for example, about 300 to 2000 keV), and the gamma ray detector 20 is incident on the detector 20 with such characteristic X-rays.
  • the imaging target TG may include three or more types of B-type nuclides that have different specific gamma ray energies. Even in this case, the imaging apparatus 1 generates event data for each event occurrence while setting “0” or “1” to the simultaneous measurement determination signal based on the beta ray detection timing and the gamma ray detection timing. From a collection of data, distributions of three or more types of B-type nuclides can be distinguished from each other and imaged at once (the same applies to other examples described later).
  • the imaging processing unit 30 can generate the distribution image 221 of the B type nuclide NC B1 and the distribution image 222 of the B type nuclide NC B2 as described above, and the beta ray detection timing by the beta ray detector 10 and the gamma ray detector 20.
  • time difference between the gamma ray detection timing is equal to or less than the predetermined time T TH (i.e., they show the simultaneous measurement determination signal is "1") and the position information (x when energy information E of E 3 has been obtained, y) can be used to generate a distribution image (not shown) of the B-type nuclide NC B3 .
  • the imaging target TG includes the A type nuclide NC A1 and the B type nuclide NC B1 described in the first embodiment, and does not include any other nuclides.
  • the gamma ray detector 20 can detect the intrinsic gamma ray emitted from the B type nuclide NC B1 , that is, the intrinsic gamma ray of energy E 1 .
  • energy information E energy information E 1
  • the included gamma ray detection signal is output. That is, the output of the gamma ray detection signal including the energy information E 1 indicates that the gamma ray detector 20 detects the intrinsic gamma ray from the B type nuclide NC B1 .
  • FIG. 13 is a conceptual diagram of image reconstruction according to the second embodiment.
  • the image 300 shows the actual distribution of the A type nuclide NC A1 and the B type nuclide NC B1 in the imaging target TG.
  • the event data group 310 is a collection of event data for events that do not involve detection of specific gamma ray event data group is a collection of event data for events that involve the detection of specific gamma ray energy E 1 311 is obtained.
  • the event data regarding the event accompanied by the detection of the intrinsic gamma ray of energy E 1 is event data including a simultaneous measurement determination signal of “1”, and includes position information (x, y) and energy information E 1 . Event data.
  • the event data for an event that does not involve detection of the inherent gamma ray is event data for an event in which the inherent gamma ray is not detected by the gamma ray detector 20, and typically includes an event including a simultaneous measurement determination signal of “0”. Therefore, it includes position information (x, y) but does not include energy information E.
  • gamma rays other than the intrinsic gamma rays may be incident on the gamma ray detector 20 due to the influence of pair annihilation, and in this case, the simultaneous measurement discrimination signal in the event data is “1”. However, if the energy information E in the event data does not indicate the energy E 1 , the event data is classified as event data for an event that does not involve detection of intrinsic gamma rays.
  • the position information (x, y) of each event data in the event data group 311 represents the detected value of the X-axis component and the Y-axis component of the position where the B type nuclide NC B1 exists in the imaging target TG. For this reason, the image reconstruction unit 32 histograms the position information (x, y) of each event data in the event data group 311 on the XY plane, thereby distributing the B type nuclide NC B1 in the imaging target TG. An image 321 can be generated.
  • the imaging processing unit 30 determines that the time difference between the beta ray detection timing by the beta ray detector 10 and the gamma ray detection timing by the gamma ray detector 20 is equal to or less than the predetermined time T TH (that is, the simultaneous measurement determination signal is “1”).
  • the distribution image 321 of the B type nuclide NC B1 can be generated using the position information (x, y) when the energy information E 1 is obtained.
  • the imaging processing unit 30 in the event data indicating that the intrinsic gamma ray of energy E 1 is detected by the gamma ray detector 20 within the predetermined time T TH from the beta ray detection timing by the beta ray detector 10.
  • the distribution image 321 of the B-type nuclide NC B1 can be generated using the position information (x, y).
  • the gamma ray detector 20 since the nuclide that emits the intrinsic gamma ray is only the B type nuclide NC B1 , the event is limited to two patterns in which the intrinsic gamma ray is detected or not detected. For this reason, the gamma ray detector 20 only needs to detect whether or not a specific gamma ray has been detected.
  • the imaging processing unit 30 determines that the time difference between the beta ray detection timing by the beta ray detector 10 and the gamma ray detection timing by the gamma ray detector 20 is equal to or less than a predetermined time T TH (that is, the simultaneous measurement determination signal is “1”). And a distribution image 321 of the B type nuclide NC B1 can be generated using the position information (x, y) of the event where the specific gamma ray is detected.
  • the position information (x, y) of each event data in the event data group 310 is the X-axis component and the Y-axis component of the existing position of either the A type nuclide NC A1 or the B type nuclide NC B1 in the imaging target TG. Represents the detected value. For this reason, the image reconstruction unit 32 histograms the position information (x, y) of each event data in the event data group 310 on the XY plane, so that all the nuclides included in the imaging target TG (here, A distribution image 320 corresponding to the superposed distribution of the nuclides NC A1 and NC B1 ) can be generated.
  • the imaging processing unit 30 can generate the distribution image 320 by using the position information (x, y) in the event data for the event that does not involve the detection of the intrinsic gamma ray.
  • the imaging processing unit 30 has not detected an intrinsic gamma ray (here, an intrinsic gamma ray of energy E 1 ) by the gamma ray detector 20 within a predetermined time T TH from the beta ray detection timing by the beta ray detector 10.
  • the distribution image 320 can be generated using the position information (x, y) in the event data indicating this.
  • the imaging apparatus 1 in the case where the imaging target TG includes the A type nuclide and the B type nuclide, only the distribution of the B type nuclide can be extracted and imaged. Images obtained by superimposing the distribution and the distribution of B-type nuclides can also be imaged.
  • a distribution image showing the distribution of only the A type nuclide NC A1 (hereinafter referred to as an A type distribution image) can be generated.
  • the pixel value (pixel luminance value) at the target position of the A type distribution image is represented by “P 320 ⁇ ((1 ⁇ ) / ⁇ ) P 321 ”.
  • P 320 represents a pixel value in the distribution image 320 at the target position
  • P 321 represents a pixel value in the distribution image 321 at the target position
  • represents an intrinsic gamma ray emitted from the nuclide in the imaging target TG. Represents the detection efficiency for energy.
  • the imaging apparatus 1 has a function of generating the distribution images 320 and 321 and the A-type distribution image, but generates only one distribution image (particularly, for example, the distribution image 321) among these distribution images. Or may operate so as to generate only two distribution images (particularly, for example, distribution image 321 and A type distribution image).
  • the gamma ray detector 20 in the second embodiment only needs to detect the presence or absence of the intrinsic gamma ray with energy E 1 , it is necessary to distinguish and detect a plurality of intrinsic gamma rays having different energies as in the first embodiment. Compared to a case, the energy resolution required for the gamma ray detector 20 may be low.
  • the detection efficiency ⁇ of intrinsic gamma rays by the gamma ray detector 20 is a value depending on energy. Further, since the detection efficiency ⁇ of the gamma ray detector 20 varies depending on the location, if the thickness of the imaging target TG in the Z-axis direction is ignored and the imaging target TG is considered to be a two-dimensional object, the detection efficiency ⁇ is determined based on the two-dimensional position ( The function depends on x, y). The detection efficiency ⁇ of intrinsic gamma rays at a position (x, y) for a certain energy is represented by ⁇ (x, y).
  • the beta ray detector 10 also has a place-dependent detection efficiency, and when quantifying radioactivity, it is necessary to correct an image obtained by histogramming the number of detections for each position. It is assumed that the image data is corrected with the detection efficiency depending on the position (the same applies to each embodiment other than this embodiment).
  • the actual two-dimensional distributions of the A type nuclide and the B type nuclide in the imaging target TG are respectively expressed as functions of the position (x, y).
  • a (x, y), and B (x, y) an image accompanied by detection of intrinsic gamma rays is represented by “ ⁇ (x, y) B (x, y)”.
  • An image that does not involve detection of intrinsic gamma rays is represented by “A (x, y) + (1 ⁇ (x, y)) B (x, y)”.
  • the beta ray emission rate of each nuclide in the imaging target TG and the probability that the B-type nuclide in the imaging target TG emits intrinsic gamma rays following the beta ray is 1.
  • y may be obtained as a distribution image of B type nuclides (distribution image 321 in this embodiment) and a distribution image of A type nuclides (A type distribution image), respectively.
  • Such a method for deriving a distribution image of each nuclide in consideration of the detection efficiency can be applied to each embodiment other than the present embodiment.
  • the imaging target TG includes the A type nuclide NC A1 and the B type nuclides NC B1 and NC B2 described in the first or second embodiment.
  • the energy of the intrinsic gamma rays of the B type nuclides NC B1 and NC B2 are E 1 and E 2 , respectively, and the energies E 1 and E 2 are different from each other.
  • Gamma ray detector 20 can detect the gamma rays to decompose the energy E 1 and energy E 2. That is, the gamma ray detector 20 can distinguish and detect the intrinsic gamma ray with energy E 1 and the intrinsic gamma ray with energy E 2 .
  • energy information E (energy information E 1 ) indicating that the energy of the detected intrinsic gamma ray is energy E 1 is obtained.
  • energy information E indicating that the energy of the detected intrinsic gamma ray is energy E 1 is obtained.
  • energy information E indicating that the energy of the detected specific gamma rays the energy E 2 (energy information E 2 ) Is output.
  • the output of the gamma ray detection signal including the energy information E 1 indicates that the intrinsic gamma ray from the B type nuclide NC B1 has been detected by the gamma ray detector 20, and the gamma ray detection signal including the energy information E 2 The output indicates that the gamma ray from the B type nuclide NC B2 is detected by the gamma ray detector 20.
  • FIG. 14 is a conceptual diagram of image reconstruction according to the third embodiment.
  • the image 400 shows the actual distribution of the A type nuclide NC A1 and the B type nuclides NC B1 and NC B2 in the imaging target TG.
  • the event data group 410 is a collection of event data for events that do not involve detection of specific gamma ray event data group is a collection of event data for events that involve the detection of specific gamma ray energy E 1 411 and an event data group 412 that is a collection of event data for events involving detection of the intrinsic gamma ray of energy E 2 .
  • the event data regarding the event accompanied by the detection of the intrinsic gamma ray of energy E 1 is event data including a simultaneous measurement determination signal of “1”, and includes position information (x, y) and energy information E 1 .
  • the event data regarding the event accompanied by the detection of the intrinsic gamma ray of energy E 2 is event data including a simultaneous measurement determination signal of “1”, and includes position information (x, y) and energy information E 2 .
  • the event data for an event that does not involve detection of the inherent gamma ray is event data for an event in which the inherent gamma ray is not detected by the gamma ray detector 20, and typically includes an event including a simultaneous measurement determination signal of “0”. Therefore, it includes position information (x, y) but does not include energy information E.
  • gamma rays other than the intrinsic gamma rays may be incident on the gamma ray detector 20 due to the influence of pair annihilation, and in this case, the simultaneous measurement discrimination signal in the event data is “1”. If the energy information E in the event data does not indicate the energy E 1 or E 2 , the event data is classified as event data for an event that does not involve detection of intrinsic gamma rays. .
  • the position information (x, y) of each event data in the event data group 411 represents the detected values of the X-axis component and the Y-axis component of the position where the B type nuclide NC B1 exists in the imaging target TG. For this reason, the image reconstruction unit 32 histograms the position information (x, y) of each event data in the event data group 411 on the XY plane, thereby distributing the B type nuclide NC B1 in the imaging target TG. An image 421 can be generated.
  • the imaging processing unit 30 determines that the time difference between the beta ray detection timing by the beta ray detector 10 and the gamma ray detection timing by the gamma ray detector 20 is equal to or less than the predetermined time T TH (that is, the simultaneous measurement determination signal is “1”).
  • the distribution image 421 of the B type nuclide NC B1 can be generated using the position information (x, y) when the energy information E 1 is obtained.
  • the imaging processing unit 30 in the event data indicating that the intrinsic gamma ray of energy E 1 is detected by the gamma ray detector 20 within the predetermined time T TH from the beta ray detection timing by the beta ray detector 10.
  • the distribution image 421 of the B-type nuclide NC B1 can be generated using the position information (x, y).
  • the position information (x, y) of each event data in the event data group 412 represents the detected values of the X-axis component and the Y-axis component of the position where the B type nuclide NC B2 is present in the imaging target TG. Therefore, the image reconstruction unit 32 histograms the position information (x, y) of each event data in the event data group 412 on the XY plane, thereby distributing the B type nuclide NC B2 in the imaging target TG. An image 422 can be generated.
  • the imaging processing unit 30 determines that the time difference between the beta ray detection timing by the beta ray detector 10 and the gamma ray detection timing by the gamma ray detector 20 is equal to or less than the predetermined time T TH (that is, the simultaneous measurement determination signal is “1”).
  • the distribution image 422 of the B type nuclide NC B2 can be generated using the position information (x, y) when the energy information E 2 is obtained.
  • the imaging processing unit 30 in the event data indicating that the intrinsic gamma ray of energy E 2 is detected by the gamma ray detector 20 within the predetermined time T TH from the beta ray detection timing by the beta ray detector 10.
  • the distribution image 422 of the B-type nuclide NC B2 can be generated using the position information (x, y).
  • the position information (x, y) of each event data in the event data group 410 is the detected value of the X-axis component and the Y-axis component at the position of any of the nuclides NC A1 , NC B1 and NC B2 in the imaging target TG. Represents. For this reason, the image reconstruction unit 32 histograms the position information (x, y) of each event data in the event data group 410 on the XY plane, so that all the nuclides (here, the TG) included in the imaging target TG are included.
  • a distribution image 420 corresponding to the superposed distribution of the nuclides NC A1 , NC B1, and NC B2 ) can be generated.
  • the imaging processing unit 30 can generate the distribution image 420 using the position information (x, y) in the event data for the event that does not involve the detection of the intrinsic gamma ray.
  • the imaging processing unit 30 detects both the intrinsic gamma ray of energy E 1 and the intrinsic gamma ray of energy E 2 by the gamma ray detector 20 within a predetermined time T TH from the beta ray detection timing by the beta ray detector 10.
  • the distribution image 420 can be generated using the position information (x, y) in the event data indicating that it has not been performed.
  • a distribution image (A type distribution image) indicating the distribution of only the A type nuclide NC A1 may be generated. it can.
  • the imaging apparatus 1 has a function of generating the distribution images 420 to 422 and the A type distribution image, but generates only one of the distribution images (particularly, for example, the distribution image 421 or 422). And only two distribution images (particularly, for example, distribution images 421 and 422, or distribution image 421 and A type distribution image, or distribution image 422 and A type distribution image) are generated. Or may operate so as to generate only three distribution images (particularly, for example, distribution images 421 and 422 and an A type distribution image).
  • the imaging target TG may include three or more B-type nuclides with different energy of the corresponding specific gamma rays. Two or more types of A-type nuclides may be included in the imaging target TG.
  • the imaging apparatus 1 can be used for such research.
  • the imaging device 1 is used by using, as the imaging target TG, the one obtained by removing the culture solution 600 from the test container in which the culture solution 600, the first cell group 610, and the second cell group 620 are accommodated. And imaging.
  • a distribution image 631 of the first B type probe that is a distribution image of the B type nuclide NC B1 and a distribution image 632 of the second B type probe that is a distribution image of the B type nuclide NC B2 are obtained. It is possible to know which of the first cell and the second cell is likely to incorporate the first B type probe, and which of the first cell and the second cell is likely to incorporate the second B type probe.
  • the first B type probe is concentrated and distributed at the position of the first cell group 610 in the test container
  • the second B type probe is concentrated at the position of the second cell group 620 in the test container.
  • the first B type probe is easily taken into the first cell
  • the second B type probe is the second cell out of the first cell and the second cell. It turns out that it is easy to be taken in by.
  • the imaging apparatus 1 is also useful at the animal experiment level.
  • a mouse is intravenously injected with a first B type probe labeled with a B type nuclide NC B1, a second B type probe labeled with a B type nuclide NC B2, and an A type probe labeled with an A type nuclide NC A1 or Orally administered, and after a lapse of a certain time, the mouse is dissected to prepare a tissue section of the mouse (for example, a tissue section having a thickness of about 50 ⁇ m). If this tissue section is used as the imaging target TG and the method described in the third embodiment is used, distribution images 420 to 422 in FIG. 14 are obtained for the tissue section.
  • the distribution images 420 to 422 (and also referring to the A type distribution image), it is possible to know which nuclide is easily accumulated in which part of the mouse. For example, in the case where cancer cells are present in a specific part of the mouse liver, if a tissue section of the mouse liver is prepared and imaged, it is only the first type B probe that has accumulated in the mouse liver, It is possible to obtain information such as whether only the second type B probe has accumulated in a specific part of the liver.
  • the living tissue is thin enough to be regarded as a two-dimensional object, there is a possibility that imaging can be performed as the imaging target TG without requiring sectioning.
  • the imaging target TG of the imaging apparatus 1 is not limited to a living body.
  • a plant can be used as the imaging target TG, and an inorganic substance can be used as the imaging target TG.
  • sodium (stable isotope) constituting the battery is used as a pre-process.
  • 23 Na) is replaced with 22 Na which is a radioactive isotope of sodium
  • a part of iron (stable isotope 56 Fe) constituting the battery is replaced with 52 Fe which is a radioactive isotope of iron.
  • 22 Na and 52 Fe function as B type nuclides NC B1 and NC B2 , respectively.
  • distribution images 220 to 222 are obtained by the method described in the first embodiment (see FIG. 12), and the distribution state of sodium and the distribution state of iron in the electrode film of the battery are individually determined based on the distribution images 221 and 222. I can know.
  • trace metals accumulated in thin films such as ion exchange membranes can be activated by beam irradiation from an accelerator, and the thin film containing trace metals after activation can be used as an imaging target TG. It is also possible to set the object subjected to the imaging target TG.
  • the present invention can also be used for particle tracking in which, for example, a probe nuclide having a certain radioactivity is introduced into a single cell and its movement is tracked as a point.
  • the usage examples given in the fourth and fifth embodiments belong to a tracer analysis in which a radioisotope tracer is included in an imaging target TG and the distribution of the tracer in the imaging target TG is analyzed.
  • the imaging apparatus 1 can also be used for activation analysis.
  • high-energy gamma rays are irradiated to an electrode film of a battery that does not contain sodium radioisotope and iron radioisotope using an accelerator.
  • an accelerator for example, in the electrode film of the battery, nuclear transformation occurs in some elements and a radioisotope is generated.
  • the electrode film in which the radioisotope is generated is set as an imaging target TG.
  • 22 Na and 52 Fe in the imaging target TG 22 Na, 52 Fe is respectively function as B type nuclide NC B1, NC B2, described in the first embodiment
  • Distribution images 220 to 222 are obtained by the method (see FIG. 12), and distribution images 221 and 222 can individually know the distribution state of sodium and the distribution state of iron in the electrode film of the battery.
  • a seventh embodiment will be described.
  • the positrons When positrons emitted by positive beta decay enter the beta-ray detector 10, the positrons give energy to the beta-ray detector 10 and eventually become almost in a stopped state (neighboring electrons (beta-ray detector 10)). Annihilation with the electrons inside. By this pair annihilation, two gamma rays (pair annihilation gamma rays) each having energy of 511 keV are generated in directions opposite to each other by 180 °. By simultaneously measuring the annihilation gamma rays together with the positrons at the detectors 10 and 20, it is possible to identify the positive beta decay.
  • positive beta decay nuclides nuclides that perform positive beta decay
  • negative beta decay nuclides nuclides that perform negative beta decay
  • the gamma ray detector 20 is disposed at a position where it can receive pair annihilation gamma rays generated by pair annihilation between a positron from a positive beta decay nuclide and an electron in the beta ray detector 10.
  • the gamma ray detector 20 detects the gamma ray by receiving the gamma ray (here, mainly the annihilation gamma ray), and outputs a gamma ray detection signal every time it detects the incidence of the gamma ray (mainly the annihilation gamma ray) to itself. .
  • the gamma ray detection signal indicates that the gamma ray is incident on the gamma ray detector 20 and includes energy information E indicating the energy of the incident gamma ray.
  • the timings t ⁇ 1 and t ⁇ 2 described above with respect to the response time (time from t ⁇ 1 to t ⁇ 4 ) of the gamma ray detector system are respectively set. It is understood that the timing at which the pair annihilation gamma rays are generated and the timing at which the pair annihilation gamma rays have interacted with the gamma ray detector 20 are obtained.
  • the simultaneous measurement discriminating unit 31 detects the gamma ray from the gamma ray detector 20 after the input timing and within a predetermined time TTH from the input timing. It is determined whether or not a signal is input, and if the determination result is “Yes” (that is, if a gamma ray is detected within a predetermined time T TH from the detection timing of the beta ray), a beta ray and a gamma ray (here In this case, it is determined that the pair annihilation gamma ray) is measured simultaneously and outputs a simultaneous measurement discrimination signal having a logical value of “1”.
  • the beta ray and gamma ray (here mainly the pair annihilation gamma ray) are measured simultaneously It is determined that it has not been performed, and a simultaneous measurement determination signal having a logical value of “0” is output.
  • FIG. 16 is an explanatory diagram of an event due to a negative beta decay (hereinafter referred to as a negative beta decay event).
  • a negative beta decay event a beta beta detection signal including position information (x, y) indicating the position of the beta decay nuclide is generated when the negative beta decay nuclide emits electrons as beta rays. No annihilation gamma rays are generated. Therefore, the event data in the negative beta decay event includes the position information (x, y) and the simultaneous measurement discrimination information of “0”, and does not include the energy information E.
  • beta rays from negative beta decay nuclides may not be detected by the beta detector 10, such an event does not form an event (the same is true for positive beta decay nuclides).
  • FIGS. 17A and 17B are explanatory diagrams of an event due to a positive beta decay (hereinafter referred to as a positive beta decay event).
  • a positive beta decay event two broken line segments with arrows extending in the opposite directions from the vicinity of the beta ray detection position represent a pair of annihilation gamma rays.
  • a positive beta decay nuclide emits a positron as a beta ray to generate a beta ray detection signal including positional information (x, y) indicating the position of the beta decay nuclide
  • the gamma ray detection signal can be generated by detecting the pair annihilation gamma ray generated by the pair annihilation between the positron from the positive beta decay nuclide and the electron in the beta ray detector 10 by the gamma ray detector 20.
  • the positive beta decay event as shown in FIG.
  • the event data includes the position information (x , Y) and simultaneous measurement discrimination information of “0”, energy information E is not included.
  • pair annihilation a pair of pair annihilation gamma rays traveling in opposite directions are generated.
  • the gamma ray detector 20 is disposed only above the imaging target TG, and among the pair of pair annihilation gamma rays, It is assumed that only one pair annihilation gamma ray can be detected by the gamma ray detector 20.
  • the pair annihilation gamma rays are used. You may make it handle as what was detected.
  • energy information E (hereinafter referred to as energy information E 511 ) indicating that the detected gamma ray energy is the energy of the pair annihilation gamma ray (511 keV).
  • a gamma ray detection signal including a signal may be output. That is, the output of the gamma ray detection signal including the energy information E 511 indicates that the annihilation gamma ray is detected by the gamma ray detector 20.
  • FIG. 18 is a conceptual diagram of image reconstruction according to the seventh embodiment.
  • the imaging target TG includes a positive beta decay nuclide NC P1 and a negative beta decay nuclide NC N1 .
  • the detection target of the beta ray detector 10 is a beta ray as a positron emitted from the positive beta decay nuclide NC P1 and a beta ray as an electron emitted from the negative beta decay nuclide NC N1 .
  • the type of radioactive decay of each of the nuclides NC P1 and NC N1 may be A type or B type.
  • the image 700 shows the actual distribution of the positive beta decay nuclide NC P1 and the negative beta decay nuclide NC N1 in the imaging target TG.
  • an event data group 710 that is a collection of event data for events that do not involve detection of pair annihilation gamma rays
  • an event data group 711 that is a collection of event data for events that involve detection of pair annihilation gamma rays, Is obtained.
  • Event data for an event involving detection of a pair of annihilation gamma rays is event data including a simultaneous measurement determination signal of “1”, and event data including position information (x, y) and energy information E 511. is there.
  • the event data for an event that does not involve detection of a pair annihilation gamma ray is event data for an event in which a pair annihilation gamma ray is not detected by the gamma ray detector 20. Therefore, the event data includes the position information (x, y) but does not include the energy information E.
  • a gamma ray other than the pair annihilation gamma ray may enter the gamma ray detector 20, and in this case, the simultaneous measurement determination signal in the event data may indicate “1”, but the energy information E in the event data If does not indicate energy E 511 , the event data is classified as event data for an event that does not involve detection of pair annihilation gamma rays.
  • the position information (x, y) of each event data in the event data group 711 represents the detected values of the X-axis component and the Y-axis component of the position where the positive beta decay nuclide NC P1 is present in the imaging target TG. Therefore, the image reconstruction unit 32 histograms the position information (x, y) of each event data in the event data group 711 on the XY plane, so that the positive beta decay nuclide NC P1 in the imaging target TG is obtained. Distribution image 721 can be generated.
  • the imaging processing unit 30 determines that the time difference between the beta ray detection timing by the beta ray detector 10 and the gamma ray detection timing by the gamma ray detector 20 is equal to or less than the predetermined time T TH (that is, the simultaneous measurement determination signal is “1”).
  • the distribution image 721 of the positive beta decay nuclide NC P1 can be generated using the position information (x, y) when the energy information E 511 is obtained.
  • the imaging processing unit 30 is position information in the event data indicating that a pair annihilation gamma ray is detected by the gamma ray detector 20 within a predetermined time T TH from the beta ray detection timing by the beta ray detector 10. (X, y) can be used to generate a distribution image 721 of the positive beta decay nuclide NC P1 .
  • the pair annihilation gamma rays are limited to two patterns, that is, whether or not they are detected in each event. Therefore, the gamma ray detector 20 only needs to detect whether or not a pair annihilation gamma ray is detected.
  • the position information (x, y) of each event data in the event data group 710 includes the X-axis component of the position of any of the negative beta decay nuclide NC N1 and the positive beta decay nuclide NC P1 in the imaging target TG, and The detected value of the Y-axis component is represented. For this reason, the image reconstruction unit 32 histograms the position information (x, y) of each event data in the event data group 710 on the XY plane, whereby all the nuclides (here, the TGs included in the imaging target TG). A distribution image 720 corresponding to the superposed distribution of the nuclides NC N1 and NC P1 ) can be generated.
  • the imaging processing unit 30 can generate the distribution image 720 using the position information (x, y) in the event data for the event that does not involve detection of the pair annihilation gamma ray.
  • the imaging processing unit 30 detects pair annihilation gamma rays (here, pair annihilation gamma rays of energy E 1 ) by the gamma ray detector 20 within a predetermined time T TH from the beta ray detection timing by the beta ray detector 10.
  • the distribution image 720 can be generated using the position information (x, y) in the event data indicating that the event has not occurred.
  • the imaging apparatus 1 when a negative beta decay nuclide and a positive beta decay nuclide are mixed in the imaging target TG, only the distribution of the positive beta decay nuclide can be extracted and imaged. At the same time, it is possible to image the superposition of the distribution of negative beta decay nuclides and the distribution of positive beta decay nuclides.
  • a distribution image indicating the distribution of only the negative beta decay nuclide NC N1 (hereinafter referred to as a negative beta decay nuclide distribution image) can be generated.
  • the pixel value (pixel luminance value) at the target position of the negative beta decay nuclide distribution image is represented by “P 720 ⁇ ((1 ⁇ ′) / ⁇ ′) P 721 ”.
  • P 720 represents the pixel value in the distribution image 720 at the target position
  • P 721 represents the pixel value in the distribution image 721 at the target position
  • ⁇ ′ represents the detection efficiency for the energy of the pair annihilation gamma rays.
  • This detection efficiency ⁇ ′ is actually a value that depends on the generation position of the pair annihilation gamma ray, but here the dependence on the position is corrected.
  • the method for deriving the distribution image of each nuclide in consideration of the detection efficiency is as described above.
  • the imaging apparatus 1 has a function of generating the distribution images 720 and 721 and the negative beta decay nuclide distribution image, but only one distribution image (particularly, for example, the distribution image 721) among these distribution images. You may operate
  • the B type nuclide NC B1 that emits intrinsic gamma rays of energy E 1 and the B type nuclide NC B2 that emits intrinsic gamma rays of energy E 2 described in the first embodiment are negative beta decay nuclides.
  • the negative beta decay nuclides NC B1 and NC B2 and the positive beta decay nuclide NC P1 are included in the imaging target TG.
  • the beta decay nuclides NC B1 and NC B2 can be identified from the energy of the respective intrinsic gamma rays.
  • the positive beta decay nuclides NC P1 can be identified by the simultaneous measurement of the pair annihilation gamma signal and beta rays by the method described in this example. It becomes.
  • the distribution image of the B type and negative beta decay nuclide NC B1 and the distribution of the B type and negative beta decay nuclide NC B2 respectively .
  • An image, a distribution image of the positive beta decay nuclide NC P1 can be generated individually.
  • the imaging target TG may further include a B-type nuclide other than the B-type nuclide NC B1 (for example, a B-type nuclide NC B2 ) or an A-type nuclide, but here, the imaging target TG other than the B-type nuclide NC B1
  • the nuclides are not included.
  • the energy of the intrinsic gamma ray of the B type nuclide NC B1 is represented by E 1 .
  • the reference imaging apparatus includes a beta ray detector 10 ′ and an imaging processing unit 30 ′ that receives a beta ray detection signal from the beta ray detector 10 ′.
  • the gamma ray detector 20 is not provided.
  • the beta ray detector 10 ′ is the same detector as the beta ray detector 10 described above, but in order to clarify the distinction from the beta ray detector 10 provided in the imaging device 1 of the present embodiment, the beta ray detector 10 ′ is used as a reference imaging device.
  • the provided beta detector is referred to by reference numeral 10 '. Radiation from the B type nuclide NC B1 included in the imaging target TG enters the beta ray detector 10 ′.
  • beta ray detector 10 ' is for detecting beta rays as radiation
  • gamma rays from the B type nuclide NC B1 may interact with the beta ray detector 10'. Radiation is detected by the line detector 10 ′ (the same applies to the beta ray detector 10). However, since the beta ray detector 10 ′ cannot distinguish between two types of radiation (beta ray and gamma ray), the beta ray detection signal output from the beta ray detector 10 ′ actually has a detection beta ray. The position information and the position information about the detected gamma rays are mixed in a state where they cannot be distinguished from each other (the same applies to the beta ray detector 10).
  • the image IM ⁇ corresponds to a superimposition of the image IM ⁇ obtained by detecting beta rays and the image IM ⁇ obtained by detecting gamma.
  • the image IM ⁇ obtained by detection of gamma rays is wider than the distribution of the actual nuclide NC B1 (ie, image IM ⁇ ) (ie, the resolution is poor). ).
  • the thickness of the scintillation detector constituting the beta ray detector 10 ′ is made sufficiently small, most of the gamma rays pass through the scintillation detector (that is, most of the gamma rays are detected by the beta ray detector 10 ′. Therefore, an image close to the image IM ⁇ can be obtained.
  • a decrease in the thickness of the scintillation detector has a limit because it causes a decrease in the detection efficiency of beta rays. Even if the thickness of the scintillation detector is reduced to the limit, the detection of gamma rays in the beta ray detector 10 'cannot be made zero.
  • the imaging apparatus 1 can generate a distribution image of the nuclide NC B1 at a higher resolution than the reference imaging apparatus by using the gamma ray detector 20.
  • An operation of the imaging apparatus 1 according to the eighth example will be described.
  • the gamma ray detector 20 can detect the intrinsic gamma rays emitted from the B type nuclide NC B1 , that is, the intrinsic gamma rays of energy E 1 .
  • energy information E energy information E 1
  • energy information E 1 indicating that the energy of the detected intrinsic gamma ray is energy E 1 is obtained.
  • the included gamma ray detection signal is output. That is, the output of the gamma ray detection signal including the energy information E 1 indicates that the gamma ray detector 20 detects the intrinsic gamma ray from the B type nuclide NC B1 .
  • FIG. 20 is a conceptual diagram of image reconstruction according to the eighth embodiment.
  • An image 800 shows an actual distribution of the B type nuclide NC B1 in the imaging target TG.
  • an event data group 810 is a collection of event data for events that do not involve detection of specific gamma ray event data group is a collection of event data for events that involve the detection of specific gamma ray energy E 1 811 is obtained.
  • the event data regarding the event accompanied by the detection of the intrinsic gamma ray of energy E 1 is event data including a simultaneous measurement determination signal of “1”, and includes position information (x, y) and energy information E 1 . Event data.
  • the event data for an event that does not involve detection of the inherent gamma ray is event data for an event in which the inherent gamma ray is not detected by the gamma ray detector 20, and typically includes an event including a simultaneous measurement determination signal of “0”. Therefore, it includes position information (x, y) but does not include energy information E.
  • gamma rays other than the intrinsic gamma rays may be incident on the gamma ray detector 20 due to the influence of pair annihilation, and in this case, the simultaneous measurement discrimination signal in the event data is “1”. However, if the energy information E in the event data does not indicate the energy E 1 , the event data is classified as event data for an event that does not involve detection of intrinsic gamma rays.
  • the position information (x, y) of each event data in the event data group 811 represents the detected value of the X-axis component and the Y-axis component of the position where the B type nuclide NC B1 exists in the imaging target TG. For this reason, the image reconstruction unit 32 histograms the position information (x, y) of each event data in the event data group 811 on the XY plane, thereby distributing the B type nuclide NC B1 in the imaging target TG. An image 821 can be generated.
  • the imaging processing unit 30 determines that the time difference between the beta ray detection timing by the beta ray detector 10 and the gamma ray detection timing by the gamma ray detector 20 is equal to or less than the predetermined time T TH (that is, the simultaneous measurement determination signal is “1”).
  • the distribution image 821 of the B type nuclide NC B1 can be generated using the position information (x, y) when the energy information E 1 is obtained.
  • the imaging processing unit 30 in the event data indicating that the intrinsic gamma ray of energy E 1 is detected by the gamma ray detector 20 within the predetermined time T TH from the beta ray detection timing by the beta ray detector 10.
  • the distribution image 821 of the B-type nuclide NC B1 can be generated using the position information (x, y).
  • the nuclide that emits intrinsic gamma rays is only the B type nuclide NC B1 , and therefore, in each event, the intrinsic gamma rays are limited to two patterns that are detected or not detected. For this reason, the gamma ray detector 20 only needs to detect whether or not a specific gamma ray has been detected.
  • the imaging processing unit 30 determines that the time difference between the beta ray detection timing by the beta ray detector 10 and the gamma ray detection timing by the gamma ray detector 20 is equal to or less than a predetermined time T TH (that is, the simultaneous measurement determination signal is “1”). And a distribution image 821 of the B-type nuclide NC B1 can be generated using the position information (x, y) of the event where the specific gamma ray is detected.
  • the position information (x, y) of each event data in the event data group 810 represents the detected values of the X-axis component and the Y-axis component of the position where the B type nuclide NC B1 is present in the imaging target TG.
  • a gamma ray (specific gamma ray by the beta ray detector 10 is included).
  • the image reconstruction unit 32 histograms the position information (x, y) of each event data in the event data group 810 on the XY plane, the beta ray detection and the gamma ray detection in the beta ray detector 10 are performed. Thus, an image 820 that is superimposed is generated.
  • the image 820 is an image that spreads over a wider range than the distribution image 821 of the B-type nuclide NC B1 generated based on the event data group 811, and is similar to the image IM ⁇ in FIG.
  • the event data that is the basis of the distribution image 821 is data on an event in which beta rays and intrinsic gamma rays are simultaneously measured, gamma rays (including intrinsic gamma rays) are incident on the beta ray detector 10 in the event. Not done. Therefore, it can be said that the image (that is, the distribution image 821) generated from the position information regarding the event is an extraction of only the beta ray distribution, and has a high resolution. Thus, according to the imaging apparatus 1, it is possible to image a single type of B-type nuclide with high resolution.
  • a thickness sufficient to reliably stop beta rays from nuclides used for imaging (minimum necessary thickness; general It is assumed that most of the gamma rays penetrate the beta ray detector 10 by preparing the beta ray detector 10 within a few millimeters), and based on the detection of the gamma ray based on the beta ray detection signal from the beta ray detector 10 It is assumed that location information is not included.
  • the contents shown in the eighth embodiment can be said to be useful in cases where the detection of gamma rays by the beta ray detector 10 results in noise that cannot be ignored.
  • the imaging apparatus 1 since the detection of the gamma ray by the beta ray detector 10 does not become zero, the imaging apparatus 1 according to the present embodiment is more suitable than the reference imaging apparatus in any case to some extent. It can be said that this contributes to higher resolution of the distribution image.
  • the beta ray detector 10 can be expressed as follows.
  • the beta ray detector 10 detects the radiation by receiving the radiation emitted from the radionuclide in the imaging target TG, and outputs a beta ray detection signal each time the incidence of the radiation on itself is detected.
  • the radiation detected by the beta ray detector 10 is mainly beta rays
  • the beta ray detection from the beta ray detector 10 is performed.
  • the signal can also be read as a radiation detection signal.
  • the radiation detection signal indicates that radiation has been incident on the beta ray detector 10 and position information indicating the radiation detection position in the beta ray detector 10. Including. Further, considering that the radiation detected by the beta ray detector 10 is a gamma ray, the above-described beta ray detection timing can be read as the radiation detection timing by the beta ray detector 10.
  • the process of forming the distribution image of the target nuclide from the beta ray detection signal and the gamma ray detection signal is basically realized by a combination of hardware and software, but is realized by the imaging processing unit 30.
  • Some of the functions may be realized by hardware or may be realized by software. is there.
  • the simultaneous measurement discrimination is realized by software, a common time stamp is added to each of the output signal data of the detector 10 and the output signal data of the detector 20 forming the event data, and the memory (not shown). The data is accumulated on the data, and simultaneous measurement discrimination is performed on the software based on the time stamp.
  • the specific function is described as a program, and the program is executed on a program execution device (for example, a microcomputer constituting the imaging processing unit 30).
  • a program execution device for example, a microcomputer constituting the imaging processing unit 30.
  • a specific function may be realized.
  • the program can be stored and fixed on an arbitrary recording medium.
  • the recording medium for storing and fixing the program may be mounted or connected to a device (server device or the like) different from the imaging processing unit 30.
  • any nuclide that performs A-type radioactive decay can be used as the A-type nuclide
  • any nuclide that performs B-type radioactive decay can be used as the B-type nuclide.
  • any one or more of 11 C, 13 N, and 18 F is set as an A type nuclide to be included in the imaging target TG.
  • the nuclides listed here are typical ones that transition to the ground state of the daughter nucleus with a 100% probability after positive beta decay, and the excited state of the daughter nucleus is highly probable due to the positive beta decay.
  • Other nuclides that transition to the ground state of the daughter nucleus without going through can also be used as A-type nuclides that undergo positive beta decay.
  • any one or more nuclides of Pb and 210 Pb can be A-type nuclides to be included in the imaging target TG.
  • the nuclides listed here are typical ones that transition to the ground state of the daughter nucleus with a negative probability of 100% after the negative beta decay.
  • Other nuclides that transit to the ground state of the daughter nucleus without going through can be used as A-type nuclides that undergo negative beta decay.
  • B-type nuclides that emit positrons by positive beta decay include 14 O, 22 Na, 34 mCl, 38 K, 44 Sc, 48 V, 52 Mn, 52 m Mn, 52 Fe, 60 Cu, 72 As. 76 Br, 82m Rb, 94m Tc, 104m Ag, 110m In, and 124 I, any one or more of the nuclides can be a B-type nuclide to be included in the imaging target TG.
  • the nuclides listed here are representative ones with a high emission rate of intrinsic gamma rays of around 90%, and can be used for other nuclides whose emission rate is lower than these.
  • B type nuclides that perform negative beta decay include 24 Na, 27 Mg, 38 S, 43 K, 46 Sc, 47 Ca, 48 Sc, 55 Mn, 60 Co, 72 Nn, 72 Ga, 78 Ge, 82 Br, 92 Sr, 95 Nb, 96 Nb, 103 Ru, 130 I, 132 I, 134 Cs, 135 Xe, 148 m Pm, 160 Tb, 192 Ir, 198 Au, and 203
  • Any one or more nuclides of Hg can be B-type nuclides to be included in the imaging target TG.
  • the nuclides listed here are representative ones with a high emission rate of intrinsic gamma rays of around 90%, and can be used for other nuclides whose emission rate is lower than these.
  • Patent Document 1 Japanese Patent No. 5526435
  • the imaging apparatus 1 employs a method of directly imaging beta rays, not only nuclides that perform positive beta decay but also nuclides that perform negative beta decay can be used. is there. This brings great benefits.
  • 24 Na which is a radioactive isotope of sodium
  • 47 Ca which is a radioactive isotope of calcium
  • the imaging apparatus 1 can be used.
  • sodium and calcium play an important role in maintaining life activity. Therefore, if the imaging apparatus 1 uses 24 Na and 47 Ca as B-type nuclides to distinguish and perform imaging, it becomes possible to simultaneously know the behavior of sodium and calcium in the living body as the imaging target TG.
  • a living body containing a plurality of types of radionuclides including 47 Ca is imaged as an imaging target TG, the behavior of calcium in the living body can be distinguished from other nuclides (the same applies to sodium).
  • radionuclides do not emit eigengamma rays in the process of transition from the parent nucleus to the daughter nucleus ground state, but do not emit eigengamma rays with a small probability (for example, a probability of less than 1%). Some radionuclides emit. Strictly speaking, such a radionuclide is a B-type nuclide, but such a radionuclide can also be handled as an A-type nuclide in the imaging apparatus 1.
  • a radionuclide treated as an A type nuclide in the imaging apparatus 1 although it may emit intrinsic gamma rays with a small probability is referred to as a pseudo A type nuclide for convenience.
  • the intrinsic gamma rays that can be emitted from the pseudo-A type nuclide are ignored in the image reconstruction process in the imaging processing unit 30.
  • the energy information E in the gamma ray detection signal for an event indicates the energy E 1 or E 2 .
  • the event data of the event may be classified into the event data group 410 by treating the event as an event that does not involve the detection of the inherent gamma ray regardless of the simultaneous measurement determination signal (however, the pseudo-A
  • the energy of the intrinsic gamma rays that can be emitted from the type nuclide is different from the energies E 1 and E 2 and can be distinguished by the gamma ray detector used).
  • nuclides that undergo B-type radioactive decay may transition to an excited state after beta decay, but have a very long half-life as a lifetime of the excited state (called a nuclear isomer or isomer).
  • the half-life of the excited state may range from less than 1 nanosecond to many years, but if the lifetime is longer than the time window of simultaneous measurement, the emission of positrons or electrons and the emission of intrinsic gamma rays Inherent gamma rays cannot be used to identify nuclides because they are not detected as continuous.
  • nuclides that emit positrons due to positive beta decay are used for imaging.
  • positrons are not necessarily emitted. This is a competitive process with the phenomenon that protons absorb the electrons that travel around and change to neutrons. That is, two processes of positron emission and electron capture compete and only one of the processes occurs. The probability of positron emission and the probability of electron capture is determined by the nuclide, but there are nuclides in which positive beta decay occurs only by electron capture and nuclides in which positive beta decay by electron capture is dominant. Not all positive beta decay nuclides are available for PET imaging.
  • negative beta decay nuclides can be used, but as described above, positron emission and electron capture are in a competitive process in positive beta decay. Therefore, the probability of emitting positrons is determined by the nuclide, and some nuclides do not emit positrons. On the other hand, since negative beta decay always emits electrons, negative beta decay nuclides have a larger number of nuclides that can be used in the present invention.
  • Negative beta decay nuclides are nuclides with excess neutrons, so they can also be produced by neutron irradiation in reactors, etc., and protons are excess nuclides, requiring irradiation of charged particles with an accelerator such as cyclotron It is often less expensive to produce than the positive beta decay nuclides.
  • the first B-type nuclide emitting a first intrinsic gamma ray having a relatively high energy and the second B-type nuclide emitting a second intrinsic gamma ray having a relatively low energy are to be imaged.
  • the spectral component due to partial energy absorption of the first intrinsic gamma ray overlaps with the spectral component due to total energy absorption of the second intrinsic gamma ray. Further, these spectral components may overlap with the spectral components due to total energy absorption or partial energy absorption of positron annihilation gamma rays (511 keV gamma rays generated by pair annihilation).
  • the continuous component due to the partial energy absorption existing before and after the peak corresponding to the total energy absorption is assumed to be the same as the continuous component overlapping the peak.
  • the background due to the overlap can be removed.
  • a gamma ray detector having a sufficiently high energy resolution is used as the gamma ray detector 20, and the peak width due to the total energy absorption event (event accompanied by the total energy absorption in the gamma ray detector 20) is narrowed, so The ground can be reduced.
  • the detection efficiency of the intrinsic gamma ray can be improved by the following manner.
  • the B-type nuclide NC B1 that emits intrinsic gamma rays of energy E 1 .
  • the intrinsic gamma ray from the nuclide NC B1 is partially absorbed by the gamma ray detector 20, energy less than the energy E 1 is detected by the gamma ray detector 20, but a predetermined energy smaller than the energy E 1 is obtained. If energy equal to or greater than the threshold is detected, it is determined that intrinsic gamma rays have been detected.
  • This energy threshold value is preferably set in a range where components due to gamma rays different from the intrinsic gamma rays are not mixed in consideration of the energy resolution of the gamma ray detector 20. Thereby, even when the partial energy absorption of the intrinsic gamma ray occurs, the intrinsic gamma ray is detected, so that it is possible to improve the detection efficiency of the intrinsic gamma ray.
  • the imaging apparatus W 1 is a beta ray two-dimensional imaging apparatus, which is in a daughter nucleus excited state due to beta decay,
  • a first nuclide for example, B type nuclide NC B1
  • a first nuclide that emits a first intrinsic gamma ray when transitioning to the ground state of the daughter nucleus following the emission, and the daughter nucleus is in an excited state by beta decay, and beta rays due to the beta decay
  • An imaging object including a second nuclide for example, a B type nuclide NC B2 ) that emits a second intrinsic gamma ray having energy different from that of the first intrinsic gamma ray when transitioning to the ground state of the daughter nucleus following the emission of (TG) detects a beta ray by receiving a beta ray based on the first nuclide or the second nuclide and detects the beta ray detection position in two dimensions.
  • nuclides that perform positive beta decay can be used, which is beneficial.
  • the imaging processing unit uses the positional information when the first intrinsic gamma ray is detected within a predetermined time from the beta ray detection timing to detect the first nuclide.
  • a distribution image may be generated, and the distribution image of the second nuclide may be generated using the position information when the second intrinsic gamma ray is detected within the predetermined time from the beta ray detection timing.
  • An imaging apparatus W 2 is a beta ray two-dimensional imaging apparatus, which is in a daughter nucleus excited state due to beta decay,
  • the first nuclide such as the B type nuclide NC B1
  • the beta ray is detected by receiving a beta ray based on the first nuclide or the second nuclide from the imaging target (TG) including the second nuclide (for example, the A type nuclide NC A1 ).
  • a beta ray detector (10) for outputting a beta ray detection signal including position information indicating two-dimensional detection positions in two dimensions, and a gamma ray detection for detecting a gamma ray including the intrinsic gamma ray.
  • the distribution image of the first nuclide (for example, the distribution image 321) and another distribution image (for example, the distribution image 320 or the A type distribution image) according to the distribution of the second nuclide are divided based on whether And an imaging processing unit (30) that can be separately generated.
  • nuclides that perform positive beta decay can be used, which is beneficial.
  • the imaging processing unit uses the position information when the intrinsic gamma ray is detected within a predetermined time from the beta ray detection timing, and the distribution image of the first nuclide. It is preferable to generate a distribution image in which the distribution of the first nuclide and the distribution of the second nuclide are superimposed using the other position information.
  • a nuclide in which negative beta decay occurs can be used for one or both of the first nuclide and the second nuclide.
  • An imaging apparatus W 3 is a beta-ray two-dimensional imaging apparatus that emits positrons as beta rays by positive beta decay (for example, From the imaging target (TG) including a positive beta decay nuclide NC P1 ) and a second nuclide that emits electrons as beta rays by negative beta decay (for example, the negative beta decay nuclide NC N1 ), the first nuclide or the above
  • a beta ray detector (10) for detecting the beta ray by receiving a beta ray based on the second nuclide, and outputting a beta ray detection signal including position information indicating the detection position of the beta ray in two dimensions;
  • a gamma ray detector (20) for detecting gamma rays including pair annihilation gamma rays generated by pair annihilation between positrons as beta rays from the first nuclide and electrons in the beta ray detector.
  • Beta ray detection timing by the beta ray detector and gamma ray detection timing by the gamma ray detector, the position information included in the beta ray detection signal, and the pair annihilation gamma ray are detected by the gamma ray detector.
  • the distribution image of the first nuclide (for example, the distribution image 721) and the other distribution image (for example, the distribution image 720 or the negative beta decay nuclide distribution image) according to the distribution of the second nuclide are divided.
  • an imaging processing unit (30) that can be separately generated.
  • nuclides that perform positive beta decay can be used, which is beneficial.
  • the imaging processing unit uses the position information when the pair annihilation gamma ray is detected within a predetermined time from the beta ray detection timing to distribute the first nuclide. It is preferable to generate an image and generate a distribution image in which the distribution of the first nuclide and the distribution of the second nuclide are superimposed using the other position information.
  • An imaging method W 1A is a beta-ray two-dimensional imaging method in which a daughter nucleus is excited by beta decay, and the daughter nucleus is brought into a ground state following beta emission by the beta decay.
  • the first nuclide that emits a first intrinsic gamma ray at the time of transition and the excited state of the daughter nucleus due to beta decay, and the transition to the ground state of the daughter nucleus following the release of beta ray by the beta decay.
  • An imaging method W 2A is a beta-ray two-dimensional imaging method, in which a daughter nucleus is excited by beta decay, and the daughter nucleus is brought into a ground state following beta emission by the beta decay.
  • the imaging target including a first nuclide that emits intrinsic gamma rays of a predetermined energy at the time of transition and a second nuclide that transitions to the ground state of the daughter nucleus due to the release of beta rays due to beta decay, the first nuclide or the A beta ray detection step of detecting a beta ray by receiving a beta ray based on the second nuclide and obtaining a beta ray detection signal including position information indicating the detection position of the beta ray in two dimensions; and the specific gamma ray A gamma ray detection step for detecting the included gamma rays, a beta ray detection timing in the beta ray detection step, and a gamma ray detection step.
  • the distribution image of the first nuclide and the second nuclide An imaging processing step for distinguishing and generating other distribution images according to the distribution.
  • An imaging method W 3A includes an imaging target including a first nuclide that emits positrons as beta rays by positive beta decay and a second nuclide that emits electrons as beta rays by negative beta decay. Detecting the beta ray by receiving a beta ray based on the first nuclide or the second nuclide and obtaining a beta ray detection signal including position information indicating the detection position of the beta ray in two dimensions
  • a gamma ray detection step for detecting a gamma ray including a pair annihilation gamma ray generated by a pair annihilation between a positron as a beta ray from the first nuclide and an electron in the beta ray detector, and the beta ray detection step.
  • Beta ray detection timing and gamma ray detection timing in the gamma ray detection step and the beta ray detection signal includes the beta ray detection signal. Based on the position information and whether the pair annihilation gamma ray is detected in the gamma ray detection step, the distribution image of the first nuclide and the other distribution image corresponding to the distribution of the second nuclide can be distinguished and generated. And an imaging process step.
  • An imaging apparatus W 4 is a beta ray two-dimensional imaging apparatus, in which a daughter nucleus is excited by beta decay, and The beta ray is detected by receiving the beta ray from the imaging target (TG) including the nuclide (for example, the B type nuclide NC B1 ) that emits an intrinsic gamma ray when transitioning to the ground state of the daughter nucleus following the emission. And a beta ray detector (10) for outputting a beta ray detection signal including position information indicating the detection position of the beta ray in two dimensions, and a gamma ray detector (20) for detecting a gamma ray including the intrinsic gamma ray.
  • TG the imaging target
  • NC B1 the nuclide
  • a beta ray detector for outputting a beta ray detection signal including position information indicating the detection position of the beta ray in two dimensions
  • a gamma ray detector (20) for detecting a gamma ray including the intrinsic gamma ray.
  • the beta ray detection timing by the beta ray detector and the gamma ray detection timing by the gamma ray detector, and the position included in the beta ray detection signal And an imaging processing unit (30) capable of generating a distribution image (for example, distribution image 821) of the nuclide based on whether the intrinsic gamma ray has been detected by the gamma ray detector.
  • the beta ray detector is for detecting beta rays from the nuclide in the imaging target, but gamma rays from the nuclide may be detected by the beta ray detector.
  • the position information in the beta ray detection signal indicates the position indicating the gamma ray detection position in addition to the position information indicating the beta ray detection position. Information is mixed in a state that is indistinguishable from each other.
  • Position information indicating the detection position of the gamma ray functions as noise, and if a gamma ray detector is not provided, the resolution of the generated nuclide distribution image is degraded. This is because the gamma ray has a longer range than the beta ray, and thus an image of a distribution that is wider than the actual nuclide distribution is generated (see FIG. 19B).
  • a nuclide distribution image can be generated based on position information when a gamma ray is detected by a gamma ray detector.
  • the intrinsic gamma ray is detected by the gamma ray detector, the intrinsic gamma ray is not incident on the beta ray detector. It is limited, and by generating a nuclide distribution image based on the detected position information, a high-resolution distribution image from which the noise is eliminated can be generated.
  • the imaging processing unit generates the nuclide distribution image using the position information when the intrinsic gamma ray is detected within a predetermined time from the beta ray detection timing. Good.
  • the imaging method W 4A includes a nuclide that emits intrinsic gamma rays when the daughter nucleus is excited by beta decay and transitions to the ground state of the daughter nucleus following beta emission by the beta decay.
  • a distribution image of the nuclide based on whether the intrinsic gamma ray was detected in the step.
  • An imaging step can be formed, and further comprising a.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Human Computer Interaction (AREA)
  • Radiology & Medical Imaging (AREA)
  • Measurement Of Radiation (AREA)
  • Nuclear Medicine (AREA)

Abstract

撮像対象(TG)には、ベータ崩壊に続いて第1、第2エネルギーを有する固有ガンマ線を放出する第1、第2核種が分布している。ベータ線検出器(10)は、ベータ線を検出してベータ線の検出位置情報を含むベータ線検出信号を出力する。ガンマ線検出器(20)は、ガンマ線を検出して検出したガンマ線のエネルギー情報を含むガンマ線検出信号を出力する。イメージング処理部(30)は、ベータ線の検出タイミングから所定時間内にガンマ線が検出され且つエネルギー情報が第1、第2エネルギーを示しているときの位置情報から、夫々、第1、第2核種の分布像を生成する。

Description

ベータ線二次元イメージング装置及び方法
 本発明は、ベータ線二次元イメージング装置及び方法に関する。
 診断装置の一つとして陽電子放射断層撮影(Positron emission tomography:PET)装置が実用化されている。PET装置では、正のベータ崩壊(β+崩壊)のうち陽電子(ポジトロン)を放出するもの(陽電子崩壊とも称されうる)を利用し、放出された陽電子が物質中の電子と対消滅した際に放出される2本のガンマ線の検出を通じて、当該ベータ崩壊を起こす放射性核種の分布を推定する。例えば、がん細胞に集積する薬剤を陽電子放出核にて標識し、その薬剤を投与した生体を撮像対象としてPET装置によりイメージングすれば、生体内でのがん細胞の3次元的な分布画像を得ることが可能である。
 ライフサイエンスや医学の進展に伴い、複数の分子動態が複雑にからみ合うことによって生体の機能が発現したり病巣が形成されたりする様子が解明されつつある。そのため、異なる放射性核種によって標識した複数の薬剤を生体に投与し、生体内の核種の分布を、核種を識別可能な態様で一度にイメージングする方法(以下、複数分子同時イメージング法という)の研究が進められており、複数分子同時イメージング法をPET装置において実現する構成が下記特許文献1に示されている。
 一方、X線等の放射線を利用した透過写真撮影において、写真フィルムに代わるものとして、フラットパネルディテクターと呼ばれる、半導体による撮像素子を平面上にアレイ状に配置したイメージング装置も用いられている。フラットパネルディテクターは、放射線による物質内の電離をダイレクトに読み出すものであり、医療の現場で普及してきている。
 上記フラットパネルディテクターに加えて、放射線の照射量による化学的変化の二次元分布をレーザー光などにより読み出し可能なイメージングプレートが基礎研究の分野で普及している。イメージングプレートによれば、簡素な構成にて撮像対象の二次元情報を得ることができ、フィルムと比較して感度が高く、また読み出し後は白色光の照射などにより情報を消去して再利用が可能であることなどから、組織切片や培養細胞中のベータ線放出核種の分布を調べるのに利用されている。しかし、これらの装置は、何れも、撮像時間中において或る場所に入射した放射線エネルギーの総量を積算するものであるため、入射した個々の放射線のエネルギー及び入射時刻は分からない。
 しかし、近年、放射線イメージング検出器において、1イベント毎に信号を処理することで、エネルギーと入射時刻を算出可能な技術(光子を読み出す場合には特に「フォトンカウンティング」と呼ばれる)が登場して来ている。この技術は、エネルギーの異なるX線を検出イベント毎にエネルギーを弁別した上で画像化することにより、エネルギーの違いによる物質のX線吸収率の違いに基づき物質を特定するCTシステムなどに応用されている。
 フォトンカウンティングタイプの計測では、撮像時間内のトータルエネルギーの積算と異なり、各撮像素子における情報をイベント毎に信号処理する必要がある。従って、撮像素子の数に対応するオンラインの1イベント信号処理回路が必要となる。解像度の高い撮像素子では、それだけチャンネル数が増えることになり、高集積且つスループットの高い信号処理装置の開発がフォトンカウンティングイメージングの高分解能化への課題となっている。しかし、近年の集積回路の高密度化などにより、年々取り扱えるチャンネル数は増えて来ており、今後益々フォトンカウンティングによるイメージングの解像度が上がることが期待できる。
 また、解像度向上のための別の手法として、シンチレーション光(可視光)の屈折又は反射を利用したレンズ等により拡大画像を得る方法もある。
特許第5526435号公報
 或る物体にベータ線を放出する放射性核種が含まれているとき、当該物体を撮像対象としてイメージングプレートによりイメージングを行えば、当該物体内における放射性核種の分布を二次元にて得ることができる。但し、従来のイメージングプレートにおいては、複数種類の核種が撮像対象内に分布しているとき、或る核種を他の核種と区別することができないため、核種を識別可能な態様でイメージングすることができない。このようなイメージングプレートにおいても、複数分子同時イメージング法が適用できれば、現在よりも更に複雑な解析が可能となるため有益である。また、ベータ崩壊により放出されるベータ線のエネルギーは核種により特定の値とはならず連続的であるため、フォトンカウンティングタイプの検出器を用いたとしても、核種の識別を行うことは困難である。
 複数種類の核種のイメージングに関わる事情について説明したが、単一種類の核種のイメージングも様々なケースにおいて必要とされ且つ有益である。この際、単一種類の核種のイメージングにおける解像度の向上は有益である。
 そこで本発明は、核種の識別が可能な態様でのイメージングを実現するベータ線二次元イメージング装置及び方法を提供することを目的の1つとする。また本発明は、核種のイメージングにおける解像度向上に寄与するベータ線二次元イメージング装置及び方法を提供することを他の目的の1つとする。
 尚、ベータ崩壊後の娘核において、原子核が励起状態にある場合、主としてガンマ線(高エネルギーの電磁波)としてエネルギーを放出して基底状態に遷移する。このガンマ線は一般的には脱励起ガンマ線と呼ばれるが、核種に固有のエネルギーを持っているため、ここでは「固有ガンマ線」と呼び、陽電子と電子の対消滅によるガンマ線と区別する。
 本発明の一側面に係るベータ線二次元イメージング装置は、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に第1固有ガンマ線を放出する第1核種と、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に前記第1固有ガンマ線とは異なるエネルギーを有する第2固有ガンマ線を放出する第2核種とを含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を出力するベータ線検出器と、ガンマ線を検出する検出器であって、前記第1固有ガンマ線及び前記第2固有ガンマ線を区別して検出するガンマ線検出器と、前記ベータ線検出器によるベータ線検出タイミング及び前記ガンマ線検出器によるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出器にて前記第1固有ガンマ線及び第2固有ガンマ線の何れかが検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布像を区別して生成可能なイメージング処理部と、を備えたことを特徴とする。
 これにより、複数種類の核種の分布を、核種を識別可能な態様で一度にイメージングすることが可能となる。また、対消滅を利用するPET装置では正のベータ崩壊(陽電子の放出を伴うベータ崩壊)を行う核種しか利用できないが、上記ベータ線二次元イメージング装置では、ベータ線を直接イメージングする方法を採用しているため、正のベータ崩壊を行う核種だけでなく負のベータ崩壊(電子の放出を伴うベータ崩壊)を行う核種も利用可能であり、有益である。
 本発明の他の側面に係るベータ線二次元イメージング装置は、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に所定エネルギーの固有ガンマ線を放出する第1核種と、ベータ崩壊によるベータ線の放出により娘核の基底状態に遷移する第2核種とを含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を出力するベータ線検出器と、前記固有ガンマ線を含むガンマ線を検出するガンマ線検出器と、前記ベータ線検出器によるベータ線検出タイミング及び前記ガンマ線検出器によるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出器にて前記固有ガンマ線が検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布に応じた他の分布像を区別して生成可能なイメージング処理部と、を備えたことを特徴とする。
 これにより、複数種類の核種の分布を、核種を識別可能な態様で(少なくとも第1核種を他の核種と識別可能な態様で)一度にイメージングすることが可能となる。また、対消滅を利用するPET装置では正のベータ崩壊(陽電子の放出を伴うベータ崩壊)を行う核種しか利用できないが、上記ベータ線二次元イメージング装置では、ベータ線を直接イメージングする方法を採用しているため、正のベータ崩壊を行う核種だけでなく負のベータ崩壊(電子の放出を伴うベータ崩壊)を行う核種も利用可能であり、有益である。
 本発明の更に他の側面に係るベータ線二次元イメージング装置は、正のベータ崩壊によりベータ線として陽電子を放出する第1核種及び負のベータ崩壊によりベータ線として電子を放出する第2核種を含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を出力するベータ線検出器と、前記第1核種からのベータ線としての陽電子と前記ベータ線検出器内の電子との対消滅により発生する対消滅ガンマ線を含むガンマ線を検出するガンマ線検出器と、前記ベータ線検出器によるベータ線検出タイミング及び前記ガンマ線検出器によるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出器にて前記対消滅ガンマ線が検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布に応じた他の分布像を区別して生成可能なイメージング処理部と、を備えたことを特徴とする。
 これにより、複数種類の核種の分布を、核種を識別可能な態様で(少なくとも第1核種を他の核種と識別可能な態様で)一度にイメージングすることが可能となる。また、対消滅を利用するPET装置では正のベータ崩壊(陽電子の放出を伴うベータ崩壊)を行う核種しか利用できないが、上記ベータ線二次元イメージング装置では、ベータ線を直接イメージングする方法を採用しているため、正のベータ崩壊を行う核種だけでなく負のベータ崩壊(電子の放出を伴うベータ崩壊)を行う核種も利用可能であり、有益である。
 本発明の更に他の側面に係るベータ線二次元イメージング装置は、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に固有ガンマ線を放出する核種を含んだ撮像対象から、前記ベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を出力するベータ線検出器と、前記固有ガンマ線を含むガンマ線を検出するガンマ線検出器と、前記ベータ線検出器によるベータ線検出タイミング及び前記ガンマ線検出器によるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出器にて前記固有ガンマ線が検出されたかと、に基づいて、前記核種の分布像を生成可能なイメージング処理部と、を備えたことを特徴とする。
 上述の如く形成されたベータ線二次元イメージング装置では、ガンマ線検出器にて固有ガンマ線が検出されたときの位置情報に基づいて核種の分布像を生成することができる。ガンマ線検出器にて固有ガンマ線が検出されたとき、ベータ線検出器に固有ガンマ線は入射していないので、ガンマ線検出器にて固有ガンマ線が検出されたときの位置情報はベータ線の検出位置情報に限定され、その検出位置情報に基づき核種の分布像を生成することで、高解像度の分布像を生成することができる。生成分布像の高解像度の理由については後述の説明からより明らかとなる。
 本発明の一側面に係るベータ線二次元イメージング方法は、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に第1固有ガンマ線を放出する第1核種と、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に前記第1固有ガンマ線とは異なるエネルギーを有する第2固有ガンマ線を放出する第2核種とを含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を得るベータ線検出ステップと、ガンマ線を検出するステップであって、前記第1固有ガンマ線及び前記第2固有ガンマ線を区別して検出するガンマ線検出ステップと、前記ベータ線検出ステップにおけるベータ線検出タイミング及び前記ガンマ線検出ステップにおけるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出ステップにて前記第1固有ガンマ線及び第2固有ガンマ線の何れかが検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布像を区別して生成可能なイメージング処理ステップと、を備えたことを特徴とする。
 本発明の他の側面に係るベータ線二次元イメージング方法は、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に所定エネルギーの固有ガンマ線を放出する第1核種と、ベータ崩壊によるベータ線の放出により娘核の基底状態に遷移する第2核種とを含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を得るベータ線検出ステップと、前記固有ガンマ線を含むガンマ線を検出するガンマ線検出ステップと、前記ベータ線検出ステップにおけるベータ線検出タイミング及び前記ガンマ線検出ステップにおけるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出ステップにて前記固有ガンマ線が検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布に応じた他の分布像を区別して生成可能なイメージング処理ステップと、を備えたことを特徴とする。
 本発明の更に他の側面に係るベータ線二次元イメージング方法は、正のベータ崩壊によりベータ線として陽電子を放出する第1核種及び負のベータ崩壊によりベータ線として電子を放出する第2核種を含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を得るベータ線検出ステップと、前記第1核種からのベータ線としての陽電子と前記ベータ線検出器内の電子との対消滅により発生する対消滅ガンマ線を含むガンマ線を検出するガンマ線検出ステップと、前記ベータ線検出ステップにおけるベータ線検出タイミング及び前記ガンマ線検出ステップにおけるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出ステップにて前記対消滅ガンマ線が検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布に応じた他の分布像を区別して生成可能なイメージング処理ステップと、を備えたことを特徴とする。
 本発明の更に他の側面に係るベータ線二次元イメージング方法は、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に固有ガンマ線を放出する核種を含んだ撮像対象から、前記ベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を得るベータ線検出ステップと、前記固有ガンマ線を含むガンマ線を検出するガンマ線検出ステップと、前記ベータ線検出ステップにおけるベータ線検出タイミング及び前記ガンマ線検出ステップにおけるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出ステップにて前記固有ガンマ線が検出されたかと、に基づいて、前記核種の分布像を生成可能なイメージングステップと、を備えたことを特徴とする。
 本発明によれば、核種の識別が可能な態様でのイメージングを実現するベータ線二次元イメージング装置及び方法を提供することが可能となる。また本発明によれば、核種のイメージングにおける解像度向上に寄与するベータ線二次元イメージング装置及び方法を提供することが可能となる。
は、本発明の実施形態で想定されるAタイプの放射性崩壊の様式を示す図である。 は、本発明の実施形態で想定されるBタイプの放射性崩壊の様式を示す図である。 は、本発明の実施形態で想定されるAタイプの放射性崩壊の具体例を示す図である。 は、本発明の実施形態で想定されるBタイプの放射性崩壊の具体例を示す図である。 は、本発明の実施形態に係るイメージング装置の概略構成図である。 は、図5のイメージング装置において、ベータ線を放射する放射性核種の存在位置とベータ線の検出位置との関係を示す図である。 は、図5のイメージング装置におけるベータ線検出器の一構成例を示す図である。 (a)及び(b)は、夫々、本発明の実施形態に係るシンチレータの構成例を示す図、及び、ベータ線の検出に関わる3つの位置の関係を示す図である。 は、本発明の実施形態に係るイメージング装置に関し、Aタイプの放射性崩壊が生じたときの信号の様子を示す図である。 (a)及び(b)は、本発明の実施形態に係るイメージング装置に関し、Bタイプの放射性崩壊が生じたときの信号の様子を示す図である。 は、本発明の実施形態に係るガンマ線検出器の構成例を示す図である。 は、本発明の第1実施例に係る画像再構成の概念図である。 は、本発明の第2実施例に係る画像再構成の概念図である。 は、本発明の第3実施例に係る画像再構成の概念図である。 は、本発明の第4実施例に係り、イメージング装置の利用例の説明図である。 は、本発明の第7実施例に係るイメージング装置に関し、負のベータ崩壊が生じたときの信号の様子を示す図である。 (a)及び(b)は、本発明の第7実施例に係るイメージング装置に関し、正のベータ崩壊が生じたときの信号の様子を示す図である。 は、本発明の第7実施例に係る画像再構成の概念図である。 (a)及び(b)は、本発明の第8実施例に係り、参考イメージング装置の構成と画像再構成の概念図である。 は、本発明の第8実施例に係る画像再構成の概念図である。
 以下、本発明の実施形態の例を、図面を参照して具体的に説明する。参照される各図において、同一の部分には同一の符号を付し、同一の部分に関する重複する説明を原則として省略する。尚、本明細書では、記述の簡略化上、情報、信号、物理量又は部材等を参照する記号又は符号を記すことによって、該記号又は符号に対応する情報、信号、物理量又は部材等の名称を省略又は略記することがある。
 図1及び図2を参照し、まず、放射性核種(放射能を持つ核種)の放射性崩壊における2つの様式を説明する。図1はAタイプの放射性崩壊の様式を示し、図2はBタイプの放射性崩壊の様式を示している。Aタイプ、Bタイプの様式にて放射性崩壊を行う放射性核種を、夫々、Aタイプ核種、Bタイプ核種と称する。また、放射性核種について、ベータ崩壊前の核種を親核と称し、ベータ崩壊後の核種を娘核と称する。
 図1に示す如く、Aタイプ核種では、親核にてベータ崩壊が生じることで娘核の基底状態に遷移する。即ち、Aタイプ核種のエネルギー準位は、ベータ崩壊によって親核のエネルギー準位から娘核の基底状態のエネルギー準位に遷移する。Aタイプ核種でのベータ崩壊は、電子(負の電荷を持ったベータ粒子)の放出を伴う負のベータ崩壊であっても良いし、陽電子(正の電荷を持ったベータ粒子)の放出を伴う正のベータ崩壊であっても良い。即ち、ベータ崩壊により、Aタイプ核種から電子又は陽電子によるベータ線が放出される。Aタイプの放射性崩壊では、後述のBタイプのそれと異なり、ガンマ線は放出されない(或いは、ガンマ線が放出されたとしても当該ガンマ線は、後述される本実施形態のイメージング装置1において無視される)。
 図2に示す如く、Bタイプ核種では、親核にてベータ崩壊が生じることで娘核の励起状態に遷移した後、続けて、ガンマ崩壊により当該Bタイプ核種に固有のエネルギーを有するガンマ線(以下、固有ガンマ線と称する)を放出することで娘核の励起状態から娘核の基底状態に遷移する。Bタイプ核種の娘核における励起状態と基底状態のエネルギー差が、固有ガンマ線のエネルギーとなる。Bタイプ核種において、固有ガンマ線が放出されるタイミングは、核の構造に依存した量子力学的な確率に従う。娘核の励起状態から娘核の基底状態への寿命の指標は半減期として与えられる。Bタイプ核種でのベータ崩壊も、Aタイプ核種でのベータ崩壊と同様、電子(負の電荷を持ったベータ粒子)の放出を伴う負のベータ崩壊であっても良いし、陽電子(正の電荷を持ったベータ粒子)の放出を伴う正のベータ崩壊であっても良い。即ち、ベータ崩壊により、Bタイプ核種から電子又は陽電子によるベータ線が放出される。
 図3に、Aタイプの放射性崩壊の一例を示す。Aタイプ核種の親核である18Fにてベータ崩壊が生じると、100%の確率で、娘核である18Oの基底状態に遷移する。このときのベータ崩壊は正のベータ崩壊であり、当該ベータ崩壊により陽電子が放出される。
 図4に、Bタイプの放射性崩壊の一例を示す。Bタイプ核種の親核である22Naは、娘核である22Neの基底状態に遷移するまでに、約99.96%の確率で22Neの励起状態を経由する。この際、ガンマ崩壊により約3.63fs(フェムト秒)の半減期にて22Neの励起状態から22Neの基底状態へと遷移し、当該遷移のときに、エネルギーが1275keV(キロエレクトロンボルト)の固有ガンマ線を放出する。22Naから22Neの励起状態への遷移をもたらすベータ崩壊は正のベータ崩壊であり当該ベータ崩壊により陽電子が放出される。尚、Bタイプ核種では、親核から娘核の基底状態に遷移するまでの過程において、娘核における複数の励起状態を経由することもある。
[イメージング装置の概略構成]
 図5に、本実施形態に係るイメージング装置(ベータ線二次元イメージング装置)1の概略構成図を示す。イメージング装置1は、位置分解型のベータ線検出器であるベータ線検出器10と、エネルギー分解型のガンマ線検出器であるガンマ線検出器20と、イメージング処理部30と、を備える。TGは、イメージング装置1における撮像対象を表している。
 実空間においてX軸、Y軸及びZ軸は互いに直交しているものとし、X軸、Y軸及びZ軸にて定義される三次元空間をXYZ空間と称する。また、X軸及びY軸に平行な二次元平面をXY平面と称する。或る位置のX軸成分、Y軸成分、Z軸成分を、夫々、x、y、zにて表すと共に、或る点のXYZ空間上での位置を(x,y,z)にて表し、且つ、X軸成分及びY軸成分のみに注目したXY平面上での位置を(x,y)にて表す。更に、説明の便宜上、Z軸の正、負に向かう向きを、夫々、上の向き、下の向きと考える。
 撮像対象TGは、X軸及びY軸方向に広がりを持つ二次元状の物体である、或いは、二次元状の撮像対象とみなすことができる物体である。撮像対象TGには、2種類以上のBタイプ核種が含まれている、或いは、1種類以上のAタイプ核種と1種類以上のBタイプ核種が含まれており、それらの核種は撮像対象TG内においてX軸及びY軸方向に分布している。
 Aタイプ核種又はBタイプ核種によって薬剤を標識することができ、Aタイプ核種及びBタイプ核種は、夫々に、薬剤に取り込まれた形態で撮像対象TGに含まれていても良い。Aタイプ核種によって標識された薬剤(即ち、Aタイプ核種が組み込まれた薬剤)をAタイププローブと称し、Bタイプ核種によって標識された薬剤(即ち、Bタイプ核種が組み込まれた薬剤)をBタイププローブと称する。
―――ベータ線検出器―――
 ベータ線検出器10は、撮像対象TG中の放射性核種から放射されるベータ線を受けることができるよう、撮像対象TGの下側において撮像対象TGに隣接配置される。ベータ線検出器10は、撮像対象TG中の放射性核種から放出されるベータ線を受けることによって当該ベータ線を検出し、自身に対するベータ線の入射を検知する度に、ベータ線検出信号を出力する。ベータ線検出信号は、ベータ線検出器10に対してベータ線の入射があったことを示すと共に、ベータ線検出器10におけるベータ線の検出位置を示す位置情報を含む。ここにおけるベータ線の検出位置はXY平面上の位置である。つまり、ベータ線検出信号中の位置情報は、ベータ線の検出位置のX軸成分及びY軸成分のみを含んでいる。
 撮像対象TGは二次元状の撮像対象であって且つベータ線の透過能力は低いため、ベータ線検出器10内の位置(x1、y1、z2)にてベータ線が検出された場合、当該ベータ線を放出した放射性核種の存在位置は撮像対象TG内の位置(x1、y1、z1)である可能性が最も高い(図6参照)。故に、X軸及びY軸成分から成る二次元成分にのみ注目した場合、ベータ線の検出位置は、ベータ線を放出する放射性核種の位置とみなすことができる。つまり、ベータ線検出信号に含まれる位置情報は撮像対象TG内のベータ線を放出する放射性核種の位置のX軸成分及びY軸成分を表す、と考えることができ、よって、ベータ線検出部10は撮像対象TG内のベータ線を放出する放射性核種の位置を二次元で検出する、とも言える。このような機能を有するベータ線検出器は公知のものであり、公知の位置分解型のベータ線検出器をベータ線検出器10として採用することができる。
 ベータ線検出器10により、ガンマ線も検出することがあるが、ガンマ線は透過力が高いため、イメージングに使用する核種からのベータ線を確実に止めるのに必要十分な厚さ(必要な最低限の厚さ;一般的には数ミリメートル以内)のベータ線検出器10を用意することにより、殆どのガンマ線はベータ線検出器10を突き抜ける。
 ここでは、図7に示す如く、ベータ線検出器10が、シンチレータ11と光検出増幅器12とから成るシンチレーション検出器にて構成されるものとする。公知の任意のシンチレーション検出器をベータ線検出器10として用いることができるが、以下、シンチレータ11及び光検出増幅器12の構成及び機能を簡単に説明する。
 シンチレータ11は、撮像対象TG中の放射性核種から放射されるベータ線を受けることができるよう、撮像対象TGの下側において撮像対象TGに隣接配置される。シンチレータ11は、X軸及びY軸方向に広がる板状のシンチレータ材料であり、ベータ線の入射によりシンチレーション光を発する。この機能を実現できる限り、シンチレータ11を構成するシンチレータ材料は任意である。
 光検出増幅器12は、シンチレータ11からのシンチレーション光を受けることができるよう、シンチレータ11の下側においてシンチレータ11に隣接配置され光学接合されている。従って、撮像対象TGと光検出増幅器12との間にシンチレータ11が挿入されることになる。光検出増幅器12は、シンチレータ11からのシンチレーション光を増幅した上で電気信号に変換して出力する。より具体的には、光検出増幅器12は、シンチレーション光が自身に入射したとき、そのシンチレーション光の入射位置(検出位置)のX軸成分及びY軸成分を示す情報を上記位置情報(x,y)として含んだ電気信号をベータ線検出信号として出力する。位置のX軸成分及びY軸成分にのみ注目したとき、光検出増幅器12におけるシンチレーション光の入射位置は、ベータ線検出器10におけるベータ線の検出位置(換言すればシンチレータ11におけるシンチレーション光の発生位置)とみなされ、撮像対象TGにおけるベータ線を放出する放射性核種の位置とみなされる。ここでは、光検出増幅器12として、ガイガーモードで動作する増倍率の高いアバランシェフォトダイオード(以下、APDと称する)をアレイ状に配置して構成されるMPPC(登録商標)が用いられるものとする。ガイガーモードとは、APDに降伏電圧以上の逆電圧を印加した状態を指す。
 ベータ線検出器10における位置検出の分解能が高まるよう、図8(a)に示す如く、板状のシンチレータ材料をX軸及びY軸方向の夫々に沿って分割して形成されるシンチレータ11a(分割によって得られたアレイ状のシンチレータ群)をシンチレータ11として用いるようにしても良い。分割によって得られた各シンチレータ材料の端面の内、XY平面に直交する端面に反射材を設けるようにしても良いし、当該端面を曇らせる(不均一な凹凸を設ける)ようにしても良い。当該端面に到達したシンチレーション光を反射又は拡散させるためである。図8(a)のシンチレータ11aを用いた場合、撮像対象TG内の或る位置(x1,y1,z1)からシンチレータ11aに向けて放出されたベータ線の多くは当該位置の真下の位置(x1,y1,z2)におけるシンチレータ材料に入射し、当該シンチレータ材料にて発生したシンチレーション光の内、真下に向かうもの以外は、当該シンチレータ材料の端面(XY平面に直交する端面)にて反射又は拡散されるなどし、結果として、当該シンチレータ材料からのシンチレーション光は位置(x1,y1,z2)の真下の位置(x1,y1,z3)に入射して位置(x1,y1,z3)にて電気信号に変換される確率が高くなる(図8(b)参照)。つまり、ベータ線を放出する放射性核種の、XY平面上の位置(x1,y1)を正確に検出することが可能となる。
 或いは、一体且つ板状のシンチレータ材料をシンチレータ11として用いつつ、当該シンチレータ材料に対しX軸及びY軸方向の夫々に平行なスリットを設けるようにしても良く、この場合でも、シンチレータ11aを用いた場合と同様の作用が得られる。この際、各スリットの端面の内、XY平面に直交する端面に到達したシンチレーション光が反射又は拡散されるよう、当該端面に反射材を設けるようにしても良いし、当該端面を曇らせる(不均一な凹凸を設ける)ようにしても良い。
 上記スリットが設けられない、一体且つ板状のシンチレータ材料をシンチレータ11として用いることも可能である。この場合、シンチレータ11内の位置(x1、y1、z2)にてシンチレーション光が発生したとき、そのシンチレーション光の光検出増幅器12での検出位置はXY平面上で相応の広がりを持つことになるが、シンチレーション光の検出強度は位置(x1、y1)において最も大きくなる確率が高いので、信号処理により、本来抽出すべき位置(x1、y1)を大まかには画像化することが可能である。
 また、MPPC(登録商標)を用いて光検出増幅器12を形成することを上述したが、ガイガーモードとは異なるノーマルモードで動作するAPDを用いて光検出増幅器12を形成しても良い。或いは、光電子増倍管にて光検出増幅器12を構成しても良い。更に或いは、シンチレータ11を用いることなくベータ線を直接検出する半導体検出器にてベータ線検出器10を構成するようにしても良い。
 近年の信号処理技術の発達により、ピクセルが細かくチャンネル数の多い半導体イメージング用検出器においても、フォトンカウンティングと同様に放射線を一つずつ数え、エネルギーと検出時刻を決定することができるシステムが登場して来ている。この様なシステムを本発明に応用することにより、将来的には十分に解像度の高い画像が得られて行くと予想される。
―――ガンマ線検出器―――
 ガンマ線検出器20は、撮像対象TG中の放射性核種から放射されるガンマ線を受けることができるよう、撮像対象TGの上側において撮像対象TGから離間して配置される。ガンマ線検出器20は、ガンマ線(主として撮像対象TG中の放射性核種からの固有ガンマ線)を受けることによって当該ガンマ線を検出し、自身に対するガンマ線(主として撮像対象TG中の放射性核種からの固有ガンマ線)の入射を検出する度に、ガンマ線検出信号を出力する。ガンマ線検出信号は、ガンマ線検出器20に対してガンマ線の入射があったことを示すと共に、入射したガンマ線のエネルギーを示すエネルギー情報Eを含む。典型的には例えば、ガンマ線検出器20に対してガンマ線の入射があったとき、当該ガンマ線のエネルギーに応じた波高値を有するパルス状の電圧信号がガンマ線検出信号として出力される。
 例えば、Ge(ゲルマニウム)による半導体検出器にてガンマ線検出器20を構成することができる。この場合、ガンマ線検出器20は、0.2~0.5%程度のエネルギー分解能にて、ガンマ線のエネルギーを検出することが可能である。
 但し、Ge以外の半導体材料(例えば、Si、CdTe、CdZnTe)を用いて構成される半導体検出器をガンマ線検出器20として用いても良いし、シンチレーション検出器を用いてガンマ線検出器20を構成しても良い。一般的にシンチレーション検出器は半導体検出器と比較してエネルギー分解能が劣るが、目的の固有ガンマ線を弁別可能なエネルギー分解能を有していれば良い。
―――イメージング処理部―――
 イメージング処理部30は、ベータ線検出器10から出力されるベータ線検出信号と、ガンマ線検出器20から出力されるガンマ線検出信号とに基づいて、撮像対象TGに含まれる対象核種の分布像を生成する。対象核種とは、撮像対象TGに含まれ且つ分布の画像化の対象となる何れかの核種を指し、例えば、撮像対象TGに含まれる何れかの特定のBタイプ核種であることもあるし、撮像対象TGにAタイプ核種とBタイプ核種が含まれている場合には当該Aタイプ核種と当該Bタイプ核種の組み合わせであることもある。
 対象核種の分布像とは、XY平面に平行な面内における対象核種の分布を示す画像を指す。例えば、XY平面に平行な所定の投影面に対して、撮像対象TG内の対象核種の存在位置を投影した像が、対象核種の分布像に相当する。対象核種としてのAタイプ核種により標識された薬剤(即ちAタイププローブ)が撮像対象TG内に分布している場合、当該対象核種の分布像は、Aタイププローブの分布像でもある。同様に、対象核種としてのBタイプ核種により標識された薬剤(即ちBタイププローブ)が撮像対象TG内に分布している場合、当該対象核種の分布像は、Bタイププローブの分布像でもある。
 イメージング処理部30は、図5に示す如く、同時計測判別部31及び画像再構成部32を備える。
 同時計測判別部31には、ベータ線検出器10の出力とガンマ線検出器20の出力とが入力される。同時計測判別部31は、ベータ線検出器10によるベータ線の検出とガンマ線検出器20によるガンマ線の検出が同時に発生したか(換言すれば、ベータ線とガンマ線が同時に計測されたか)を判別し、その判別結果を示す同時計測判別信号を出力する。ここにおける同時とは、所定の時間幅を有する概念であり、ベータ線検出器10によるベータ線の検出タイミングとガンマ線検出器20によるガンマ線の検出タイミングとの時間差が所定値以下であれば、それらの検出タイミングは同時であると判別され、そうでなければ、それらの検出タイミングは同時でないと判別される。このような判別を、以下、同時計測判別と称する。
 撮像対象TG内のAタイプ核種又はBタイプ核種よりベータ線が放出されるタイミングtβ1と、ベータ線検出器10にて当該ベータ線が検出器10と相互作用を起こしたタイミングtβ2と、ベータ線検出器10から当該検出結果を示すベータ線検出信号が出力されるタイミングtβ3と、当該ベータ線検出信号が同時計測判別部31に入力されるタイミングtβ4との間には有限の時間差がある(時間の経過と共に、タイミングtβ1、tβ2、tβ3、tβ4が、この順番で訪れる)。タイミングtβ1からタイミングtβ4までの時間は、検出器10の応答及び検出器10からの信号線の長さ等により決まり、検出器10を含むベータ線検出器系に固有の長さを持つ。また、検出位置の違いによる時間応答の差や、信号ノイズの影響により、実際にベータ線が放出されたタイミングtβ1から同時計測判別部31にベータ線検出信号が入力されるタイミングtβ4までの時間差は時間的広がりを持つ。この時間差の広がりはどれだけ精確に放射線放出時刻を決定できるかを示しており、一般的に時間分解能と呼ばれる。
 同様に、実際に撮像対象TG内のBタイプ核種よりガンマ線(固有ガンマ線)が放出されるタイミングtγ1と、ガンマ線検出器20にて当該ガンマ線が検出器20と相互作用を起こしたタイミングtγ2と、ガンマ線検出器20から当該検出結果を示すガンマ線検出信号が出力されるタイミングtγ3と、当該ガンマ線検出信号が同時計測判別部31に入力されるタイミングtγ4との間には有限の時間差がある(時間の経過と共に、タイミングtγ1、tγ2、tγ3、tγ4が、この順番で訪れる)。ベータ線検出器系と同様に、タイミングtγ1及びtγ4間の時間差は、検出器20を含むガンマ線検出器系に固有の長さと固有の時間分解能を持つ。
 ベータ線検出器系の応答時間(tβ1からtβ4までの時間)とガンマ線検出器系の応答時間(tγ1からtγ4までの時間)が同じになるように、応答時間が短い方の検出器系に対し信号遅延回路を入れることで、同時計測判別部31は、自身に対するベータ線検出信号及びガンマ線検出信号の入力タイミングの差を、ベータ線検出器10によるベータ線の検出タイミングとガンマ線検出器20によるガンマ線の検出タイミングとの時間差であるとみなすことができる。上記信号遅延回路を考慮すれば、タイミングtβ4、tγ4を、夫々、ベータ線検出器10(ベータ線検出器系)によるベータ線の検出タイミング、ガンマ線検出器20(ガンマ線検出器系)によるガンマ線の検出タイミングと解することもできる。同時計測判別部31は、ベータ線検出器10からベータ線検出信号が入力されたとき、その入力タイミングよりも後であって且つ当該入力タイミングから所定時間TTH内にガンマ線検出器20からガンマ線検出信号が入力されたか否かを判定し、その判定結果が“肯定”であれば(即ち、ベータ線の検出タイミングから所定時間TTH内にガンマ線が検出されれば)、ベータ線とガンマ線(固有ガンマ線を含む)が同時に計測されたと判別して “1”の論理値を有する同時計測判別信号を出力し、そうでなければ、ベータ線とガンマ線(固有ガンマ線を含む)が同時に計測されていないと判別して“0”の論理値を有する同時計測判別信号を出力する。
 ベータ線が放出されてから固有ガンマ線が放出されるまでには時間差が存在するが、娘核の励起状態が特別に長い半減期を持たない限りは、その時間はフェムト秒からピコ秒のオーダーであり、一般的な放射線検出器の時間分解能よりも短い。また、ベータ線の飛行時間、ガンマ線の飛行時間、陽電子が放出されてから対消滅を起こすまでの時間もそれぞれ有限であるが、それらは各検出器の時間分解能よりも十分に短い。従って、同時計測判別に使われる時間TTHは、主としてベータ線検出器10とガンマ線検出20についての各時間分解能に依存して決定される。即ち、ベータ線検出器10によるベータ線の検出時刻決定精度の時間的広がりと、ガンマ線検出器20によるガンマ線の検出時刻決定精度の時間的広がりとを考慮し、それらの時間的広がりを含むように時間TTHの長さが決定される。
 ベータ線検出器10にてベータ線が検出される1つの事象をイベントと称する。イメージ処理部30は、イベントごとに、ベータ線検出器10の出力とガンマ線検出器20の出力とに基づいてイベントデータを作成する。
 イベントデータは、イメージング装置1に内包されるメモリ(例えばイメージ処理部30に内包されるメモリ;不図示)に記録されて良い。イベントデータを作成及び記録するための条件(マスタートリガー)は、ベータ線検出器10において放射線(ベータ線)の検出があったかどうかでのみ設定するのが効率的である。この場合、ガンマ線検出器20にのみ放射線(ガンマ線)の検出があった場合は、イベントデータは記録されないことになるが、ガンマ線検出器20からの情報のみでは画像化を行うことはできないので、無駄なデータを取り込まないことになり効率的である。但し、ベータ線検出器10とガンマ線検出器20を独立に動作させ、それぞれの検出器からのデータ(イベントデータを形成するデータ)に共通のタイムスタンプを付加して、測定後オフラインで同時計測判別を行う場合は、それぞれの検出器による放射線検出がイベントデータ記録のマスタートリガーとなる。
 1つのイベントに対する1つのイベントデータは、当該イベントについてのベータ線検出信号中の位置情報(x,y)と、当該イベントについての同時計測判別信号とを含み、更に、当該イベントについてのガンマ線検出信号中のエネルギー情報Eを含み得る。
 1つのイベントに対する1つのイベントデータに関し、当該イベントデータに含まれる位置情報(x,y)は、当該イベントにおけるベータ線検出器10でのベータ線の検出位置を示すと共に、当該イベントにてベータ線を放出した核種の位置を示す。これらの位置は、上述したようにX軸方向及びY軸方向における位置である。
 1つのイベントに対する1つのイベントデータに関し、当該イベントデータに含まれ得るエネルギー情報Eは、当該イベントについて検出されたガンマ線のエネルギーを示す。
 図9は、Aタイプ核種のベータ崩壊によるイベント(以下、Aタイプイベントと称する)の説明図である。Aタイプイベントでは、Aタイプ核種がベータ線を放出することで当該Aタイプ核種の位置を示す位置情報(x,y)を含んだベータ線検出信号が生成されるが、固有ガンマ線は放出されない。従って、Aタイプイベントにおけるイベントデータは、位置情報(x,y)と“0”の同時計測判別情報とから成り、エネルギー情報Eを含まない。尚、Aタイプ核種のベータ崩壊によるベータ線がベータ線検出器10にて検出されないこともあるが、そのような事象はイベントを形成しない(Bタイプ核種についても同様)。
 図10(a)及び(b)は、Bタイプ核種のベータ崩壊によるイベント(以下、Bタイプイベントと称する)の説明図である。Bタイプイベントでは、Bタイプ核種がベータ線を放出することで当該Bタイプ核種の位置を示す位置情報(x,y)を含んだベータ線検出信号が生成され、一方で、当該Bタイプ核種がベータ線に続いて固有ガンマ線を放出することでガンマ線検出信号が生成され得る。Bタイプイベントにおいて、図10(a)に示す如く、固有ガンマ線がガンマ線検出器20に入射して固有ガンマ線が検出された場合にはガンマ線検出信号が生成されるため、イベントデータはベータ線検出信号中の位置情報(x,y)と、“1”の同時計測判別情報と、ガンマ線検出信号中のエネルギー情報Eと、を含むことになる。一方、Bタイプイベントにおいて、図10(b)に示す如く、固有ガンマ線がガンマ線検出器20に入射しなかった場合には固有ガンマ線が検出されないため、イベントデータは、位置情報(x,y)と“0”の同時計測判別情報とから成り、エネルギー情報Eを含まない。
 ガンマ線検出器20における固有ガンマ線の検出の立体角を高めれば、Bタイプイベントにおいて固有ガンマ線の検出確率が高まるため、当該立体角が高まるように、ガンマ線検出器20を撮像対象TGに対してなるだけ近接して配置すると共に、撮像対象TGを覆うのに十分な大きさ持つガンマ線検出器20を使用すると良い。イメージング装置1では、特許文献1(特許第5526435号公報)に示されたようなPET装置での固有ガンマ線の検出の立体角と比べても、2πに近い大きな立体角を容易に実現できる。
 また複数のガンマ線検出器を使用することもできる。例えば図11に示す如く、ガンマ線検出器20を、撮像対象TGの上方に配置されたガンマ線検出器20Aと撮像対象TGの下方に配置されたガンマ線検出器20Bとで構成することで、上記立体角を更に高めることも可能である。この場合、撮像対象TG及びベータ線検出器10が、ガンマ線検出器20Aとガンマ線検出器20Bとの間に挟まれるようになる。図11の例では、ガンマ線検出器20Aとガンマ線検出器20Bとで1つのガンマ線検出器20が構成され、ガンマ線検出器20Aの出力とガンマ線検出器20Bの出力の論理和にてガンマ線検出器20の出力信号が形成される。即ち、図11のガンマ線検出器20は、ガンマ線検出器20Aの存在位置及びガンマ線検出器20Bの存在位置の何れかにてガンマ線(固有ガンマ線を含む)が検出されたとき、ガンマ線検出信号を出力する。
 画像再構成部32は、複数のイベントデータに基づいて画像の再構成処理を行うことで対象核種の分布像を生成する。実際には、得られたイベントデータに含まれる(x,y)のベータ線検出位置情報を、ピクセル化されたXY平面上にヒストグラミングして行くことで、放射能の強度分布を得る。
 上述したイメージング装置1の構成及び動作等を、便宜上、基本実施例と称する。以下、複数の実施例の中で、生成される分布像の具体例、イメージング装置1の更なる詳細な構成例、動作例、応用例などを説明する。特に記述無き限り且つ矛盾無き限り、基本実施例に記載の事項が後述の各実施例に適用され、各実施例において基本実施例と矛盾する事項については、各実施例での記載が優先される。また矛盾無き限り、以下に述べる複数の実施例の内、任意の実施例に記載した事項を、他の任意の実施例に適用することもできる(即ち複数の実施例の内の任意の2以上の実施例を組み合わせることも可能である)。
<<第1実施例>>
 第1実施例を説明する。第1実施例では、撮像対象TGに、互いに種類の異なるBタイプ核種であるBタイプ核種NCB1及びBタイプ核種NCB2が含まれているものとする。撮像対象TGに、Aタイプ核種が更に含まれていても良いが、ここではAタイプ核種が含まれていないものとする。Bタイプ核種NCB1、NCB2の固有ガンマ線のエネルギーは、夫々、E1、E2であり、エネルギーE1及びE2は互いに異なる。ガンマ線検出器20は、エネルギーE1とエネルギーE2を分解してガンマ線を検出することができる。即ち、ガンマ線検出器20は、エネルギーE1の固有ガンマ線とエネルギーE2の固有ガンマ線を区別して検出することができる。
 ガンマ線検出器20にて、Bタイプ核種NCB1から放出された固有ガンマ線が検出されたとき、検出された固有ガンマ線のエネルギーがエネルギーE1であることを示すエネルギー情報E(以下、エネルギー情報E1と称することもある)を含んだガンマ線検出信号が出力され、Bタイプ核種NCB2から放出された固有ガンマ線が検出されたとき、検出された固有ガンマ線のエネルギーがエネルギーE2であることを示すエネルギー情報E(以下、エネルギー情報E2と称することもある)を含んだガンマ線検出信号が出力される。つまり、エネルギー情報E1を含むガンマ線検出信号の出力は、Bタイプ核種NCB1からの固有ガンマ線がガンマ線検出器20にて検出されたことを示しており、エネルギー情報E2を含むガンマ線検出信号の出力は、Bタイプ核種NCB2からの固有ガンマ線がガンマ線検出器20にて検出されたことを示している。
 図12は、第1実施例に係る画像再構成の概念図である。像200は、撮像対象TG内におけるBタイプ核種NCB1及びNCB2の実際の分布を示している。第1実施例では、固有ガンマ線の検出を伴わないイベントについてのイベントデータの集まりであるイベントデータ群210と、エネルギーE1の固有ガンマ線の検出を伴うイベントについてのイベントデータの集まりであるイベントデータ群211と、エネルギーE2の固有ガンマ線の検出を伴うイベントについてのイベントデータの集まりであるイベントデータ群212と、が得られる。
 エネルギーE1の固有ガンマ線の検出を伴うイベントについてのイベントデータとは、“1”の同時計測判別信号を含むイベントデータであって、且つ、位置情報(x,y)及びエネルギー情報E1を含むイベントデータである。
 エネルギーE2の固有ガンマ線の検出を伴うイベントについてのイベントデータとは、“1”の同時計測判別信号を含むイベントデータであって、且つ、位置情報(x,y)及びエネルギー情報E2を含むイベントデータである。
 固有ガンマ線の検出を伴わないイベントについてのイベントデータとは、固有ガンマ線がガンマ線検出器20にて検出されないイベントについてのイベントデータであって、典型的には“0”の同時計測判別信号を含むイベントデータであり、従って、位置情報(x,y)を含むがエネルギー情報Eは含まない。後述されるように、対消滅等の影響により固有ガンマ線以外のガンマ線がガンマ線検出器20に入射することもあり、その場合には、イベントデータ中の同時計測判別信号が“1”を示すこともあるが、当該イベントデータ中のエネルギー情報EがエネルギーE1又はE2を示していないならば、当該イベントデータは、固有ガンマ線の検出を伴わないイベントについてのイベントデータに分類される。
 尚、固有ガンマ線を放出しない核種からのベータ線が検出されたときと同時に、偶発的に、エネルギーE1又はE2と同じエネルギーのガンマ線が検出されることもある。この偶発同時計数と呼ばれる事象は、同時計測を行う装置においては、単位時間当たりの検出回数に依存して頻度が変化するものの常に起こる事象である。イメージング装置1における偶発同時計数については、PETで用いられる公知の遅延同時計測によるバックグラウンド除去と同等の手法により補正可能である。
 イベントデータ群211中の各イベントデータの位置情報(x,y)は、撮像対象TG中のBタイプ核種NCB1の存在位置のX軸成分及びY軸成分の検出値を表している。このため、画像再構成部32は、イベントデータ群211中の各イベントデータの位置情報(x,y)をXY平面上にヒストグラミングすることにより、撮像対象TG内におけるBタイプ核種NCB1の分布像221を生成することができる。
 つまり、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングとガンマ線検出器20によるガンマ線検出タイミングとの時間差が所定時間TTH以下であって(即ち、同時計測判別信号が“1”を示していて)且つエネルギー情報E1が得られているときの位置情報(x,y)を用いて、Bタイプ核種NCB1の分布像221を生成することができる。換言すれば、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングから所定時間TTH内にガンマ線検出器20にてエネルギーE1の固有ガンマ線が検出されていることを示すイベントデータ中の位置情報(x,y)を用いて、Bタイプ核種NCB1の分布像221を生成することができる。
 イベントデータ群212中の各イベントデータの位置情報(x,y)は、撮像対象TG中のBタイプ核種NCB2の存在位置のX軸成分及びY軸成分の検出値を表している。このため、画像再構成部32は、イベントデータ群212中の各イベントデータの位置情報(x,y)をXY平面上にヒストグラミングすることにより、撮像対象TG内におけるBタイプ核種NCB2の分布像222を生成することができる。
 つまり、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングとガンマ線検出器20によるガンマ線検出タイミングとの時間差が所定時間TTH以下であって(即ち、同時計測判別信号が“1”を示していて)且つエネルギー情報E2が得られているときの位置情報(x,y)を用いて、Bタイプ核種NCB2の分布像222を生成することができる。換言すれば、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングから所定時間TTH内にガンマ線検出器20にてエネルギーE2の固有ガンマ線が検出されていることを示すイベントデータ中の位置情報(x,y)を用いて、Bタイプ核種NCB2の分布像222を生成することができる。
 イベントデータ群210中の各イベントデータの位置情報(x,y)は、撮像対象TG中のBタイプ核種NCB1及びNCB2の何れかの存在位置のX軸成分及びY軸成分の検出値を表している。このため、画像再構成部32は、イベントデータ群210中の各イベントデータの位置情報(x,y)をXY平面上にヒストグラミングすることにより、撮像対象TGに含まれる全ての核種(ここではBタイプ核種NCB1及びNCB2)の分布を重ね合わせたものに相当する分布像220を生成することができる。つまり、イメージング処理部30は、固有ガンマ線の検出を伴わないイベントについてのイベントデータ中の位置情報(x,y)を用いて分布像220を生成することができる。
 このように、イメージング装置1によれば、複数種類のBタイプ核種の分布を、互いに区別して一度にイメージングすることが可能となる。
 但し、イメージング装置1は、分布像220、221及び222を生成する機能を有しつつも、それらの分布像の内、1つの分布像(特に例えば分布像221又は222)のみを生成するように動作しても良いし、2つの分布像(特に例えば分布像221及び222)のみを生成するように動作しても良い。
 尚、Bタイプ核種からベータ線として陽電子が放出される場合、当該陽電子がシンチレータ11にてシンチレーション光の生成に寄与した後、当該陽電子と近傍の電子との間で対消滅が発生する。そして、対消滅により発生した511keVのエネルギーを有するガンマ線がガンマ線検出器20に入射することもある。従って、ベータ崩壊によって陽電子を放出するBタイプ核種を撮像対象TGに含める場合においては、対消滅によるガンマ線と固有ガンマ線との区別を可能とするべく、撮像対象TGに含まれる各Bタイプ核種として、固有ガンマ線のエネルギーが511keVと異なるBタイプ核種を使用すると良い。そして、エネルギーE1とエネルギーE2と511keVとを分解できる程度のエネルギー分解能をガンマ線検出器20に持たせると良い。後述の他の実施例においても同様である。
 上述の如く、正のベータ崩壊により陽電子を放出する核種について、当該陽電子がベータ線検出器10に入射すると、最終的には対消滅を起こし、互いに180°反対向きに、各々に511keVのエネルギーを持った2本のガンマ線(対消滅ガンマ線)を生じる。この対消滅ガンマ線がガンマ線検出器20にて検出されることもある。正のベータ崩壊を行う核種と負のベータ崩壊を行う核種を同時にイメージングする際は、この対消滅ガンマ線の検出の有無により、核種の識別を行うことも可能である。これは後述の他の実施例においても同様であるが、この方法については第7実施例にて詳説する。
 また、ベータ線がシンチレータ11に入射したとき特性X線が発生することがあり、その特性X線がガンマ線検出器20に入射することもある。特性X線のエネルギー(数keV~百数十keV程度)は固有ガンマ線のエネルギー(例えば300~2000keV程度)よりも低く、ガンマ線検出器20は、そのような特性X線の検出器20への入射を阻止するために、低エネルギーであるX線を吸収する吸収体を検出器前に有していると良い。後述の他の実施例においても同様である。
 撮像対象TGに2種類のBタイプ核種が含まれている場合について説明したが、撮像対象TGに、対応する固有ガンマ線のエネルギーが互いに異なる3種類以上のBタイプ核種が含まれていても良い。この場合においても、イメージング装置1は、ベータ線検出タイミング及びガンマ線検出タイミングに基づいて同時計測判別信号に“0”又は“1”を設定しつつ、イベントの発生ごとにイベントデータを生成し、イベントデータの集まりから、3種類以上のBタイプ核種の分布を互いに区別して一度にイメージングすることができる(後述の他の実施例においても同様)。
 より具体的には例えば、Bタイプ核種NCB1及びNCB2に加えて、対応する固有ガンマ線のエネルギーがE3であるBタイプ核種NCB3が撮像対象TGに含まれているとき、イメージング処理部30は、上述の如くBタイプ核種NCB1の分布像221及びBタイプ核種NCB2の分布像222を生成することが可能であると共に、ベータ線検出器10によるベータ線検出タイミングとガンマ線検出器20によるガンマ線検出タイミングとの時間差が所定時間TTH以下であって(即ち、同時計測判別信号が“1”を示していて)且つE3のエネルギー情報Eが得られているときの位置情報(x,y)を用いて、Bタイプ核種NCB3の分布像(不図示)を生成することができる。
<<第2実施例>>
 第2実施例を説明する。第2実施例では、撮像対象TGに、Aタイプ核種NCA1と、第1実施例でも述べたBタイプ核種NCB1が含まれており、それら以外の核種は含まれていないものとする。ガンマ線検出器20は、Bタイプ核種NCB1から放出される固有ガンマ線を検出することができる、即ち、エネルギーE1の固有ガンマ線を検出することができる。ガンマ線検出器20にて、Bタイプ核種NCB1から放出された固有ガンマ線が検出されたとき、検出された固有ガンマ線のエネルギーがエネルギーE1であることを示すエネルギー情報E(エネルギー情報E1)を含んだガンマ線検出信号が出力される。つまり、エネルギー情報E1を含むガンマ線検出信号の出力は、Bタイプ核種NCB1からの固有ガンマ線がガンマ線検出器20にて検出されたことを示している。
 図13は、第2実施例に係る画像再構成の概念図である。像300は、撮像対象TG内におけるAタイプ核種NCA1及びBタイプ核種NCB1の実際の分布を示している。第2実施例では、固有ガンマ線の検出を伴わないイベントについてのイベントデータの集まりであるイベントデータ群310と、エネルギーE1の固有ガンマ線の検出を伴うイベントについてのイベントデータの集まりであるイベントデータ群311と、が得られる。
 エネルギーE1の固有ガンマ線の検出を伴うイベントについてのイベントデータとは、“1”の同時計測判別信号を含むイベントデータであって、且つ、位置情報(x,y)及びエネルギー情報E1を含むイベントデータである。
 固有ガンマ線の検出を伴わないイベントについてのイベントデータとは、固有ガンマ線がガンマ線検出器20にて検出されないイベントについてのイベントデータであって、典型的には“0”の同時計測判別信号を含むイベントデータであり、従って、位置情報(x,y)を含むがエネルギー情報Eは含まない。第1実施例で述べたように、対消滅等の影響により固有ガンマ線以外のガンマ線がガンマ線検出器20に入射することもあり、その場合には、イベントデータ中の同時計測判別信号が“1”を示すこともあるが、当該イベントデータ中のエネルギー情報EがエネルギーE1を示していないならば、当該イベントデータは、固有ガンマ線の検出を伴わないイベントについてのイベントデータに分類される。
 イベントデータ群311中の各イベントデータの位置情報(x,y)は、撮像対象TG中のBタイプ核種NCB1の存在位置のX軸成分及びY軸成分の検出値を表している。このため、画像再構成部32は、イベントデータ群311中の各イベントデータの位置情報(x,y)をXY平面上にヒストグラミングすることにより、撮像対象TG内におけるBタイプ核種NCB1の分布像321を生成することができる。
 つまり、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングとガンマ線検出器20によるガンマ線検出タイミングとの時間差が所定時間TTH以下であって(即ち、同時計測判別信号が“1”を示していて)且つエネルギー情報E1が得られているときの位置情報(x,y)を用いて、Bタイプ核種NCB1の分布像321を生成することができる。換言すれば、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングから所定時間TTH内にガンマ線検出器20にてエネルギーE1の固有ガンマ線が検出されていることを示すイベントデータ中の位置情報(x,y)を用いて、Bタイプ核種NCB1の分布像321を生成することができる。
 第2実施例では、固有ガンマ線を放出する核種はBタイプ核種NCB1のみであるため、各イベントにおいて、固有ガンマ線は検出されるか検出されないかの2パターンに限定される。このため、ガンマ線検出器20は、固有ガンマ線が検出されたか否かのみを検出すれば足る。そして、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングとガンマ線検出器20によるガンマ線検出タイミングとの時間差が所定時間TTH以下であって(即ち、同時計測判別信号が“1”を示していて)且つ固有ガンマ線が検出されたイベントの位置情報(x,y)を用いて、Bタイプ核種NCB1の分布像321を生成できる。
 イベントデータ群310中の各イベントデータの位置情報(x,y)は、撮像対象TG中のAタイプ核種NCA1及びBタイプ核種NCB1の何れかの存在位置のX軸成分及びY軸成分の検出値を表している。このため、画像再構成部32は、イベントデータ群310中の各イベントデータの位置情報(x,y)をXY平面上にヒストグラミングすることにより、撮像対象TGに含まれる全ての核種(ここでは核種NCA1及びNCB1)の分布を重ね合わせたものに相当する分布像320を生成することができる。つまり、イメージング処理部30は、固有ガンマ線の検出を伴わないイベントについてのイベントデータ中の位置情報(x,y)を用いて分布像320を生成することができる。換言すれば、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングから所定時間TTH内にガンマ線検出器20にて固有ガンマ線(ここではエネルギーE1の固有ガンマ線)が検出されていないことを示すイベントデータ中の位置情報(x,y)を用いて、分布像320を生成することができる。
 このように、イメージング装置1によれば、撮像対象TGにAタイプ核種とBタイプ核種が混在している場合において、Bタイプ核種の分布のみを抽出してイメージングでき、それと同時に、Aタイプ核種の分布とBタイプ核種の分布を重ねあわせたものもイメージングできる。
 また、分布像320と分布像321を合成することで、Aタイプ核種NCA1のみの分布を示す分布像(以下、Aタイプ分布像と称する)を生成することもできる。この合成において、Aタイプ分布像の注目位置における画素値(画素の輝度値)は、“P320-((1-ε)/ε)P321 ”によって表される。ここで、P320は注目位置における分布像320中の画素値を表し、P321は注目位置における分布像321中の画素値を表し、εは撮像対象TG中の核種から放出される固有ガンマ線のエネルギーに対する検出効率を表す。
 但し、イメージング装置1は、分布像320及び321並びにAタイプ分布像を生成する機能を有しつつも、それらの分布像の内、1つの分布像(特に例えば分布像321)のみを生成するように動作しても良いし、2つの分布像(特に例えば分布像321及びAタイプ分布像)のみを生成するように動作しても良い。
 尚、第2実施例におけるガンマ線検出器20は、エネルギーE1の固有ガンマ線の有無のみを検出すれば足るため、第1実施例の如くエネルギーが互いに異なる複数の固有ガンマ線を区別して検出する必要がある場合と比べ、ガンマ線検出器20に要求されるエネルギー分解能は低くて良い。
 ここで、検出効率εに関して説明を補足する。ガンマ線検出器20による固有ガンマ線の検出効率εはエネルギーに依存した値となる。また、ガンマ線検出器20の検出効率εは場所により変わるため、Z軸方向における撮像対象TGの厚みを無視して撮像対象TGが2次元物体であると考えると、検出効率εは2次元位置(x,y)に依存した関数となる。あるエネルギーに対する位置(x,y)での固有ガンマ線の検出効率εを、ε(x,y)で表す。
 また、ベータ線検出器10も場所依存の検出効率を持ち、放射能の定量化を行う際は位置毎の検出回数をヒストグラミングして得た画像に対し補正を行う必要となるが、得られた画像データは、位置に依存した検出効率で補正されているものとする(本実施例以外の各実施例においても同様)。
 撮像対象TGにAタイプ核種とBタイプ核種が含まれている場合において、撮像対象TG内のAタイプ核種、Bタイプ核種の実際の二次元分布を、位置(x,y)の関数として、夫々、A(x,y)、B(x,y)で表すと、固有ガンマ線の検出を伴う画像は“ε(x,y)B(x,y)”で表される。また、固有ガンマ線の検出を伴わない画像は“A(x,y)+(1-ε(x,y))B(x,y)”で表される。ここで、撮像対象TG中の各核種のベータ線の放出率、及び、撮像対象TG中のBタイプ核種がベータ線に続けて固有ガンマ線を放出する確率は、1であると仮定している。
 従って、ベータ線と固有ガンマ線の同時検出を伴う画像(本実施例においてイベントデータ群311にて形成される画像)をT(x,y)とし、且つ、固有ガンマ線の検出を伴わない画像(本実施例においてイベントデータ群310にて形成される画像)をD(x,y)とすると、撮像対象TG中のBタイプ核種の実際の分布は“B(x,y)=T(x,y)/ε(x,y)”となる。また、画像の差し引きにより、撮像対象TG中のAタイプ核種の実際の分布を求めることができ、撮像対象TG中のAタイプ核種の実際の分布は“A(x,y)=D(x,y)-T(x,y)(1-ε(x,y))/ε(x,y)”となる。実測によりD(x,y)及びT(x,y)が得られるので、検出効率ε(x,y)を考慮しつつ、上記式を用いて、B(x,y)、A(x,y)を、夫々、Bタイプ核種の分布像(本実施例において分布像321)、Aタイプ核種の分布像(Aタイプ分布像)として求めれば良い。このような検出効率を考慮して各核種の分布像を導出する方法は、本実施例以外の各実施例においても適用可能である。
<<第3実施例>>
 第3実施例を説明する。第3実施例では、撮像対象TGに、第1又は第2実施例で述べたAタイプ核種NCA1並びにBタイプ核種NCB1及びNCB2が含まれているものとする。第1実施例で述べたように、Bタイプ核種NCB1、NCB2の固有ガンマ線のエネルギーは、夫々、E1、E2であり、エネルギーE1及びE2は互いに異なる。ガンマ線検出器20は、エネルギーE1とエネルギーE2を分解してガンマ線を検出することができる。即ち、ガンマ線検出器20は、エネルギーE1の固有ガンマ線とエネルギーE2の固有ガンマ線を区別して検出することができる。
 ガンマ線検出器20にて、Bタイプ核種NCB1から放出された固有ガンマ線が検出されたとき、検出された固有ガンマ線のエネルギーがエネルギーE1であることを示すエネルギー情報E(エネルギー情報E1)を含んだガンマ線検出信号が出力され、Bタイプ核種NCB2から放出された固有ガンマ線が検出されたとき、検出された固有ガンマ線のエネルギーがエネルギーE2であることを示すエネルギー情報E(エネルギー情報E2)を含んだガンマ線検出信号が出力される。つまり、エネルギー情報E1を含むガンマ線検出信号の出力は、Bタイプ核種NCB1からの固有ガンマ線がガンマ線検出器20にて検出されたことを示しており、エネルギー情報E2を含むガンマ線検出信号の出力は、Bタイプ核種NCB2からの固有ガンマ線がガンマ線検出器20にて検出されたことを示している。
 図14は、第3実施例に係る画像再構成の概念図である。像400は、撮像対象TG内におけるAタイプ核種NCA1並びにBタイプ核種NCB1及びNCB2の実際の分布を示している。第3実施例では、固有ガンマ線の検出を伴わないイベントについてのイベントデータの集まりであるイベントデータ群410と、エネルギーE1の固有ガンマ線の検出を伴うイベントについてのイベントデータの集まりであるイベントデータ群411と、エネルギーE2の固有ガンマ線の検出を伴うイベントについてのイベントデータの集まりであるイベントデータ群412と、が得られる。
 エネルギーE1の固有ガンマ線の検出を伴うイベントについてのイベントデータとは、“1”の同時計測判別信号を含むイベントデータであって、且つ、位置情報(x,y)及びエネルギー情報E1を含むイベントデータである。
 エネルギーE2の固有ガンマ線の検出を伴うイベントについてのイベントデータとは、“1”の同時計測判別信号を含むイベントデータであって、且つ、位置情報(x,y)及びエネルギー情報E2を含むイベントデータである。
 固有ガンマ線の検出を伴わないイベントについてのイベントデータとは、固有ガンマ線がガンマ線検出器20にて検出されないイベントについてのイベントデータであって、典型的には“0”の同時計測判別信号を含むイベントデータであり、従って、位置情報(x,y)を含むがエネルギー情報Eは含まない。第1実施例で述べたように、対消滅等の影響により固有ガンマ線以外のガンマ線がガンマ線検出器20に入射することもあり、その場合には、イベントデータ中の同時計測判別信号が“1”を示すこともあるが、当該イベントデータ中のエネルギー情報EがエネルギーE1又はE2を示していないならば、当該イベントデータは、固有ガンマ線の検出を伴わないイベントについてのイベントデータに分類される。
 イベントデータ群411中の各イベントデータの位置情報(x,y)は、撮像対象TG中のBタイプ核種NCB1の存在位置のX軸成分及びY軸成分の検出値を表している。このため、画像再構成部32は、イベントデータ群411中の各イベントデータの位置情報(x,y)をXY平面上にヒストグラミングすることにより、撮像対象TG内におけるBタイプ核種NCB1の分布像421を生成することができる。
 つまり、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングとガンマ線検出器20によるガンマ線検出タイミングとの時間差が所定時間TTH以下であって(即ち、同時計測判別信号が“1”を示していて)且つエネルギー情報E1が得られているときの位置情報(x,y)を用いて、Bタイプ核種NCB1の分布像421を生成することができる。換言すれば、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングから所定時間TTH内にガンマ線検出器20にてエネルギーE1の固有ガンマ線が検出されていることを示すイベントデータ中の位置情報(x,y)を用いて、Bタイプ核種NCB1の分布像421を生成することができる。
 イベントデータ群412中の各イベントデータの位置情報(x,y)は、撮像対象TG中のBタイプ核種NCB2の存在位置のX軸成分及びY軸成分の検出値を表している。このため、画像再構成部32は、イベントデータ群412中の各イベントデータの位置情報(x,y)をXY平面上にヒストグラミングすることにより、撮像対象TG内におけるBタイプ核種NCB2の分布像422を生成することができる。
 つまり、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングとガンマ線検出器20によるガンマ線検出タイミングとの時間差が所定時間TTH以下であって(即ち、同時計測判別信号が“1”を示していて)且つエネルギー情報E2が得られているときの位置情報(x,y)を用いて、Bタイプ核種NCB2の分布像422を生成することができる。換言すれば、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングから所定時間TTH内にガンマ線検出器20にてエネルギーE2の固有ガンマ線が検出されていることを示すイベントデータ中の位置情報(x,y)を用いて、Bタイプ核種NCB2の分布像422を生成することができる。
 イベントデータ群410中の各イベントデータの位置情報(x,y)は、撮像対象TG中の核種NCA1、NCB1及びNCB2の何れかの存在位置のX軸成分及びY軸成分の検出値を表している。このため、画像再構成部32は、イベントデータ群410中の各イベントデータの位置情報(x,y)をXY平面上にヒストグラミングすることにより、撮像対象TGに含まれる全ての核種(ここでは核種NCA1、NCB1及びNCB2)の分布を重ね合わせたものに相当する分布像420を生成することができる。つまり、イメージング処理部30は、固有ガンマ線の検出を伴わないイベントについてのイベントデータ中の位置情報(x,y)を用いて分布像420を生成することができる。換言すれば、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングから所定時間TTH内にガンマ線検出器20にてエネルギーE1の固有ガンマ線及びエネルギーE2の固有ガンマ線の何れも検出されていないことを示すイベントデータ中の位置情報(x,y)を用いて、分布像420を生成することができる。
 このように、イメージング装置1によれば、複数種類のBタイプ核種の分布を、互いに区別して一度にイメージングすることが可能であると共に、それと同時に、Aタイプ核種の分布とBタイプ核種の分布を重ねあわせたものもイメージングできる。また、第2実施例で述べた方法と同様の方法にて分布像420~422を合成することで、Aタイプ核種NCA1のみの分布を示す分布像(Aタイプ分布像)を生成することもできる。
 但し、イメージング装置1は、分布像420~422及びAタイプ分布像を生成する機能を有しつつも、それらの分布像の内、1つの分布像(特に例えば分布像421又は422)のみを生成するように動作しても良いし、2つの分布像(特に例えば分布像421及び422、或いは、分布像421及びAタイプ分布像、或いは、分布像422及びAタイプ分布像)のみを生成するように動作しても良いし、3つの分布像(特に例えば分布像421及び422並びにAタイプ分布像)のみを生成するように動作しても良い。
 これは、特に、後述の放射化分析を行う場合等であって且つ撮像対象TG内に放射化された複数の様々な核種が生成されている状態において、固有ガンマ線を放出する特定のBタイプ核種のみを、固有ガンマ線のエネルギーの検出結果に基づき抽出してイメージングするケース等に有益である。
 尚、第1実施例で述べたように、撮像対象TGに、対応する固有ガンマ線のエネルギーが互いに異なる3種類以上のBタイプ核種が含まれていても良い。また、2種類以上のAタイプ核種が撮像対象TGに含まれていても良い。
<<第4実施例>>
 第4実施例を説明する。第4実施例では、イメージング装置1の利用例を説明する。
 薬の開発過程において、どの薬剤が生体中のどの部位に良く集積するのか等の挙動を調べる際、細胞レベル(ペトリ皿等の平面上に培養した細胞)から研究することが一般的に行われる。このような研究にイメージング装置1を利用できる。
 例えば、図15に示す如く、まず前工程として、Bタイプ核種NCB1により標識された第1Bタイププローブ及びBタイプ核種NCB2により標識された第2Bタイププローブを含んだ培養液600中に、或る特性を有する第1細胞の集まりから成る第1細胞群610と他の特性を有する第2細胞の集まりから成る第2細胞群620を所定時間浸す。これにより、第1Bタイププローブ又は第2Bタイププローブが、第1又は第2細胞の何れかに取り込まれることが期待される。
 上記の前工程の後、培養液600、第1細胞群610及び第2細胞群620が収容されていた試験容器中から培養液600を除去したものを撮像対象TGとして用いて、イメージング装置1にてイメージングを行う。この結果、Bタイプ核種NCB1の分布像である第1Bタイププローブの分布像631及びBタイプ核種NCB2の分布像である第2Bタイププローブの分布像632が得られ、分布像631及び632から、第1細胞及び第2細胞のどちらに第1Bタイププローブが取り込まれ易いのか、第1細胞及び第2細胞のどちらに第2Bタイププローブが取り込まれ易いのかを知ることができる。例えば、試験容器中の第1細胞群610の存在位置に第1Bタイププローブが集中して分布していること及び試験容器中の第2細胞群620の存在位置に第2Bタイププローブが集中して分布していることが分布像631及び632から読み取れる場合には、第1細胞及び第2細胞の内、第1Bタイププローブは第1細胞の方に取り込まれ易く且つ第2Bタイププローブは第2細胞の方に取り込まれ易いということが分かる。
 動物実験レベルにおいてもイメージング装置1は有益である。例えば、マウスに対し、Bタイプ核種NCB1により標識された第1Bタイププローブ、Bタイプ核種NCB2により標識された第2Bタイププローブ及びAタイプ核種NCA1により標識されたAタイププローブを静脈注射又は経口投与し、一定時間の経過後、当該マウスを解剖して当該マウスの組織切片(例えば厚さ50μm程度の組織切片)を作製する。この組織切片を撮像対象TGとして用い、第3実施例で述べた方法を用いれば、当該組織切片に対して図14の分布像420~422が得られる。分布像420~422を参照することで(更にはAタイプ分布像も参照することで)、マウスのどの部位に、どの核種が集積し易いのか等を知ることができる。例えば、マウスの肝臓の特定部分に癌細胞がある場合において、マウスの肝臓の組織切片を作製してイメージングを行えば、マウスの肝臓に集積していたのは第1Bタイププローブのみであるとか、肝臓の特定部位に集積していたのは第2Bタイププローブのみであるか等の情報を得ることができる。
 また、二次元物体とみなせる程度に厚さの小さい生体組織であれば、切片化を必要とすることなく、撮像対象TGとしてイメージングできる可能性もある。
<<第5実施例>>
 第5実施例を説明する。第5実施例では、イメージング装置1の他の利用例を説明する。
 イメージング装置1の撮像対象TGは生体に限定されず、例えば植物を撮像対象TGにすることもできるし、無機物を撮像対象TGにすることもできる。
 例えば、ナトリウムと鉄を含んで構成される電池に関し、当該電池の電極膜におけるナトリウムの分布や鉄の分布を知りたいという要望がある場合、前工程として、当該電池を構成するナトリウム(安定同位体の23Na)の一部をナトリウムの放射性同位体である22Naに置換すると共に、当該電池を構成する鉄(安定同位体の56Fe)の一部を鉄の放射性同位体である52Feに置換する。その後、適宜、当該電池に対して所定の劣化試験を施してから、当該電池の電極膜を撮像対象TGとする。この場合、22Na、52Feが、夫々、Bタイプ核種NCB1、NCB2として機能する。
 そうすると、第1実施例で述べた方法により分布像220~222が得られ(図12参照)、分布像221及び222により、当該電池の電極膜におけるナトリウムの分布状態及び鉄の分布状態を個別に知ることができる。
 この他、イオン交換膜等の薄膜に蓄積した微量金属等を加速器からのビーム照射により放射化させ、該放射化後の微量金属等を含む薄膜を撮像対象TGにすることもできるし、イオンドープが施された物体を撮像対象TGにすることもできる。また、分布画像のみではなく、例えば単一の細胞に対し或る放射能量のプローブ核種を導入し、点としてその動きを追跡するパーティクルトラッキングに本発明を利用することもできる。
<<第6実施例>>
 第6実施例を説明する。第4及び第5実施例で挙げた利用例は、放射性同位体によるトレーサーを撮像対象TGに含ませておいて当該トレーサーが撮像対象TG内でどのように分布したのかを分析するトレーサー分析に属するが、イメージング装置1を放射化分析に利用することもできる。
 例えば、第5実施例で述べた電池の内、ナトリウムの放射性同位体及び鉄の放射性同位体を含んでいない状態の電池の電極膜に対し、加速器を用いて高エネルギーのガンマ線を照射する。これにより、当該電池の電極膜において、幾つかの元素に核変換が生じて放射性同位体が生成される。この放射性同位体が生成された状態の電極膜を撮像対象TGとする。そうすると例えば、撮像対象TG内に22Na及び52Feが含まれている場合には、22Na、52Feが、夫々、Bタイプ核種NCB1、NCB2として機能し、第1実施例で述べた方法により分布像220~222が得られて(図12参照)、分布像221及び222により当該電池の電極膜におけるナトリウムの分布状態及び鉄の分布状態を個別に知ることができる。
<<第7実施例>>
 第7実施例を説明する。正のベータ崩壊で放出された陽電子がベータ線検出器10に入射した場合、陽電子はベータ線検出器10にエネルギーを与え、最終的にはほぼ停止状態になり近隣の電子(ベータ線検出器10内の電子)と対消滅を起こす。この対消滅により、互いに180°反対向きに、各々に511keVのエネルギーを持った2本のガンマ線(対消滅ガンマ線)が発生する。この消滅ガンマ線を陽電子と共に検出器10及び20にて同時計測することにより、正のベータ崩壊であったことを同定することができる。即ち、正のベータ崩壊を行う核種(以下、正のベータ崩壊核種と称する)と負のベータ崩壊を行う核種(以下、負のベータ崩壊核種と称する)が混在している撮像対象TGにおいて、消滅ガンマ線の計測の有無により、これらを区別して同定することが可能である。この方法について説明を加える。
 ガンマ線検出器20は、正のベータ崩壊核種からの陽電子とベータ線検出器10内の電子との対消滅により発生した対消滅ガンマ線を受けることができるような位置に配置されている。ガンマ線検出器20は、ガンマ線(ここでは主として対消滅ガンマ線)を受けることによって当該ガンマ線を検出し、自身に対するガンマ線(ここでは主として対消滅ガンマ線)の入射を検知する度に、ガンマ線検出信号を出力する。上述してきたように、ガンマ線検出信号は、ガンマ線検出器20に対してガンマ線の入射があったことを示すと共に、入射したガンマ線のエネルギーを示すエネルギー情報Eを含む。尚、ガンマ線検出器20にて検出されるべきガンマ線が対消滅ガンマ線である場合、ガンマ線検出器系の応答時間(tγ1からtγ4までの時間)に関して上述したタイミングtγ1、tγ2は、夫々、対消滅ガンマ線が発生したタイミング、当該対消滅ガンマ線がガンマ線検出器20と相互作用を起こしたタイミングであると解される。
 同時計測判別部31は、ベータ線検出器10からベータ線検出信号が入力されたとき、その入力タイミングよりも後であって且つ当該入力タイミングから所定時間TTH内にガンマ線検出器20からガンマ線検出信号が入力されたか否かを判定し、その判定結果が“肯定”であれば(即ち、ベータ線の検出タイミングから所定時間TTH内にガンマ線が検出されれば)、ベータ線とガンマ線(ここでは主として対消滅ガンマ線)が同時に計測されたと判別して “1”の論理値を有する同時計測判別信号を出力し、そうでなければ、ベータ線とガンマ線(ここでは主として対消滅ガンマ線)が同時に計測されていないと判別して“0”の論理値を有する同時計測判別信号を出力する。
 図16は、負のベータ崩壊によるイベント(以下、負のベータ崩壊イベントと称する)の説明図である。負のベータ崩壊イベントでは、負のベータ崩壊核種がベータ線としての電子を放出することで当該ベータ崩壊核種の位置を示す位置情報(x,y)を含んだベータ線検出信号が生成されるが、対消滅ガンマ線は生成されない。従って、負のベータ崩壊イベントにおけるイベントデータは、位置情報(x,y)と“0”の同時計測判別情報とから成り、エネルギー情報Eを含まない。尚、負のベータ崩壊核種からのベータ線がベータ線検出器10にて検出されないこともあるが、そのような事象はイベントを形成しない(正のベータ崩壊核種についても同様)。
 図17(a)及び(b)は、正のベータ崩壊によるイベント(以下、正のベータ崩壊イベントと称する)の説明図である。図17(a)及び(b)において、ベータ線検出位置付近から互いに反対向きに伸びる2本の矢印付き破線線分が一対の対消滅ガンマ線を表している。正のベータ崩壊イベントでは、正のベータ崩壊核種がベータ線としての陽電子を放出することで当該ベータ崩壊核種の位置を示す位置情報(x,y)を含んだベータ線検出信号が生成され、一方で、正のベータ崩壊核種からの陽電子とベータ線検出器10内の電子との対消滅により発生した対消滅ガンマ線がガンマ線検出器20にて検出されることでガンマ線検出信号が生成され得る。正のベータ崩壊イベントにおいて、図17(a)に示す如く、対消滅ガンマ線がガンマ線検出器20に入射して対消滅ガンマ線が検出された場合にはガンマ線検出信号が生成されるため、イベントデータは、ベータ線検出信号中の位置情報(x,y)と、“1”の同時計測判別情報と、ガンマ線検出信号中のエネルギー情報Eと、を含むことになる。一方、正のベータ崩壊イベントにおいて、図17(b)に示す如く、対消滅ガンマ線がガンマ線検出器20に入射しなかった場合には対消滅ガンマ線が検出されないため、イベントデータは、位置情報(x,y)と“0”の同時計測判別情報とから成り、エネルギー情報Eを含まない。
 対消滅により、互いに反対向きに進行する一対の対消滅ガンマ線が発生するが、ここでは、撮像対象TGの上方にのみガンマ線検出器20が配置されているものとし、一対の対消滅ガンマ線の内、一方の対消滅ガンマ線のみがガンマ線検出器20にて検出されうるものとする。但し、図11に示したように撮像対象TGの上方及び下方の夫々にガンマ線検出器を配置し、双方のガンマ線検出器にて一対の対消滅ガンマ線が同時に検出されたときに限り、対消滅ガンマ線が検出されたものとして取り扱うようにしても良い。
 ガンマ線検出器20にて、対消滅ガンマ線がガンマ線として検出されたとき、検出されたガンマ線のエネルギーが対消滅ガンマ線のエネルギー(511keV)であることを示すエネルギー情報E(以下、エネルギー情報E511と称することもある)を含んだガンマ線検出信号が出力される。つまり、エネルギー情報E511を含むガンマ線検出信号の出力は、対消滅ガンマ線がガンマ線検出器20にて検出されたことを示している。
 図18は、第7実施例に係る画像再構成の概念図である。ここでは、撮像対象TGに、正のベータ崩壊核種NCP1及び負のベータ崩壊核種NCN1が含まれているものとする。そうすると、ベータ線検出器10の検出対象は、正のベータ崩壊核種NCP1から放出される陽電子としてのベータ線と、負のベータ崩壊核種NCN1から放出される電子としてのベータ線となる。核種NCP1及びNCN1の夫々の放射性崩壊の様式はAタイプでもBタイプでも良い。
 像700は、撮像対象TG内における正のベータ崩壊核種NCP1及び負のベータ崩壊核種NCN1の実際の分布を示している。第7実施例では、対消滅ガンマ線の検出を伴わないイベントについてのイベントデータの集まりであるイベントデータ群710と、対消滅ガンマ線の検出を伴うイベントについてのイベントデータの集まりであるイベントデータ群711と、が得られる。
 対消滅ガンマ線の検出を伴うイベントについてのイベントデータとは、“1”の同時計測判別信号を含むイベントデータであって、且つ、位置情報(x,y)及びエネルギー情報E511を含むイベントデータである。
 対消滅ガンマ線の検出を伴わないイベントについてのイベントデータとは、対消滅ガンマ線がガンマ線検出器20にて検出されないイベントについてのイベントデータであって、典型的には“0”の同時計測判別信号を含むイベントデータであり、従って、位置情報(x,y)を含むがエネルギー情報Eは含まない。対消滅ガンマ線以外のガンマ線がガンマ線検出器20に入射することもあり、その場合には、イベントデータ中の同時計測判別信号が“1”を示すこともあるが、当該イベントデータ中のエネルギー情報EがエネルギーE511を示していないならば、当該イベントデータは、対消滅ガンマ線の検出を伴わないイベントについてのイベントデータに分類される。
 イベントデータ群711中の各イベントデータの位置情報(x,y)は、撮像対象TG中の正のベータ崩壊核種NCP1の存在位置のX軸成分及びY軸成分の検出値を表している。このため、画像再構成部32は、イベントデータ群711中の各イベントデータの位置情報(x,y)をXY平面上にヒストグラミングすることにより、撮像対象TG内における正のベータ崩壊核種NCP1の分布像721を生成することができる。
 つまり、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングとガンマ線検出器20によるガンマ線検出タイミングとの時間差が所定時間TTH以下であって(即ち、同時計測判別信号が“1”を示していて)且つエネルギー情報E511が得られているときの位置情報(x,y)を用いて、正のベータ崩壊核種NCP1の分布像721を生成することができる。換言すれば、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングから所定時間TTH内にガンマ線検出器20にて対消滅ガンマ線が検出されていることを示すイベントデータ中の位置情報(x,y)を用いて、正のベータ崩壊核種NCP1の分布像721を生成することができる。
 第7実施例では、陽電子を放出する核種は正のベータ崩壊核種NCP1のみであるため、各イベントにおいて、対消滅ガンマ線は検出されるか検出されないかの2パターンに限定される。このため、ガンマ線検出器20は、対消滅ガンマ線が検出されたか否かのみを検出すれば足る。
 イベントデータ群710中の各イベントデータの位置情報(x,y)は、撮像対象TG中の負のベータ崩壊核種NCN1及び正のベータ崩壊核種NCP1の何れかの存在位置のX軸成分及びY軸成分の検出値を表している。このため、画像再構成部32は、イベントデータ群710中の各イベントデータの位置情報(x,y)をXY平面上にヒストグラミングすることにより、撮像対象TGに含まれる全ての核種(ここでは核種NCN1及びNCP1)の分布を重ね合わせたものに相当する分布像720を生成することができる。つまり、イメージング処理部30は、対消滅ガンマ線の検出を伴わないイベントについてのイベントデータ中の位置情報(x,y)を用いて分布像720を生成することができる。換言すれば、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングから所定時間TTH内にガンマ線検出器20にて対消滅ガンマ線(ここではエネルギーE1の対消滅ガンマ線)が検出されていないことを示すイベントデータ中の位置情報(x,y)を用いて、分布像720を生成することができる。
 このように、イメージング装置1によれば、撮像対象TGに負のベータ崩壊核種と正のベータ崩壊核種が混在している場合において、正のベータ崩壊核種の分布のみを抽出してイメージングでき、それと同時に、負のベータ崩壊核種の分布と正のベータ崩壊核種の分布を重ねあわせたものもイメージングできる。
 また、分布像720と分布像721を合成することで、負のベータ崩壊核種NCN1のみの分布を示す分布像(以下、負のベータ崩壊核種分布像と称する)を生成することもできる。この合成において、負のベータ崩壊核種分布像の注目位置における画素値(画素の輝度値)は、“P720-((1-ε’)/ε’)P721 ”によって表される。ここで、P720は注目位置における分布像720中の画素値を表し、P721は注目位置における分布像721中の画素値を表し、ε’は対消滅ガンマ線のエネルギーに対する検出効率を表す。この検出効率ε’は実際には、対消滅ガンマ線の発生位置に依存した値となるが、ここでは位置による依存は補正されているものとしている。検出効率を考慮して各核種の分布像を導出する方法については上述した通りである。
 但し、イメージング装置1は、分布像720及び721並びに負のベータ崩壊核種分布像を生成する機能を有しつつも、それらの分布像の内、1つの分布像(特に例えば分布像721)のみを生成するように動作しても良いし、2つの分布像(特に例えば分布像721及び負のベータ崩壊核種分布像)のみを生成するように動作しても良い。
 本実施例で示した方法を、他の実施例の何れかと組わせることも可能である。例えば、第1実施例で述べた、エネルギーE1の固有ガンマ線を放出するBタイプ核種NCB1及びエネルギーE2の固有ガンマ線を放出するBタイプ核種NCB2が負のベータ崩壊核種であるとした上で、負のベータ崩壊核種NCB1及びNCB2と正のベータ崩壊核種NCP1とを撮像対象TG内に含ませる。そうすると、ベータ崩壊核種NCB1及びNCB2については、夫々の固有ガンマ線のエネルギーから同定可能となる。また、負のベータ崩壊核種NCB1及びNCB2は陽電子を放出しないので、本実施例で述べた方法により対消滅ガンマ信号とベータ線との同時計測にて正のベータ崩壊核種NCP1を同定可能となる。つまり、ベータ線検出器10からのベータ線検出信号及びガンマ線検出器20からのガンマ線検出信号に基づく、エネルギーE1の固有ガンマ線の検出を伴うイベントについてのイベントデータ、エネルギーE2の固有ガンマ線の検出を伴うイベントについてのイベントデータ、対消滅ガンマ線の検出を伴うイベントについてのイベントデータより、夫々、Bタイプ且つ負のベータ崩壊核種NCB1の分布像、Bタイプ且つ負のベータ崩壊核種NCB2の分布像、正のベータ崩壊核種NCP1の分布像を個別に生成することができる。
<<第8実施例>>
 第8実施例を説明する。イメージング装置1を用いて複数種類の核種を識別可能な態様でイメージングする方法を上述したが、イメージング装置1は、単一種類の核種のイメージングにおいても有用である。これについて説明する。
 第8実施例では、撮像対象TGに上述のBタイプ核種NCB1が含まれているものとする。撮像対象TGに、Bタイプ核種NCB1以外のBタイプ核種(例えばBタイプ核種NCB2)やAタイプ核種が更に含まれていても良いが、ここでは、撮像対象TGにBタイプ核種NCB1以外の核種は含まれていないものとする。上述したように、Bタイプ核種NCB1の固有ガンマ線のエネルギーはE1で表される。
 イメージング装置1によるBタイプ核種NCB1のイメージング方法の説明に先立ち、イメージング装置1との対比に供される参考イメージング装置を説明する。図19(a)に示す如く、参考イメージング装置には、ベータ線検出器10’と、ベータ線検出器10’からのベータ線検出信号を受けるイメージング処理部30’とが設けられているが、ガンマ線検出器20は設けられていない。ベータ線検出器10’は上述のベータ線検出器10と同じ検出器であるが、本実施形態のイメージング装置1に設けられるベータ線検出器10との区別を明確にするため、参考イメージング装置に設けられるベータ線検出器を符号10’にて参照する。撮像対象TGに含まれるBタイプ核種NCB1からの放射線がベータ線検出器10’に入射する。
 ベータ線検出器10’はベータ線を放射線として検出するためのものではあるが、Bタイプ核種NCB1からのガンマ線がベータ線検出器10’と相互作用を起こすこともあり、このときにもベータ線検出器10’にて放射線が検出される(ベータ線検出器10でも同様)。但し、ベータ線検出器10’は2種類の放射線(ベータ線とガンマ線)を識別できないため、ベータ線検出器10’から出力されるベータ線検出信号には、実際には、検出ベータ線についての位置情報と検出ガンマ線についての位置情報が互いに区別できない状態で混在する(ベータ線検出器10でも同様)。故に、参考イメージング装置において、ベータ線検出器10’からのベータ線検出信号に基づきイメージング処理部30’にて画像再構成を行った場合、図19(b)に示すような、ベータ線の検出とガンマ線の検出を重ね合わせたような画像IMβγがBタイプ核種NCB1の分布像として生成されることとなる。
 画像IMβγは、ベータ線の検出による画像IMβとガンマの検出による画像IMγとを重ね合わせたものに相当する。一般に、ガンマ線はベータ線よりも飛程が長いため、ガンマ線の検出による画像IMγは、実際の核種NCB1の分布よりも(即ち画像IMβ)よりも広範囲に広がったものとなる(即ち解像度が悪い)。
 ベータ線検出器10’を構成するシンチレーション検出器の厚みを十分に小さくすれば、殆どのガンマ線がシンチレーション検出器を透過するようになるため(即ち殆どのガンマ線がベータ線検出器10’にて検出されなくなるため)、画像IMβに近い画像を得ることができる。しかしながら、シンチレーション検出器の厚みの減少はベータ線の検出効率の低下を招くため限界がある。シンチレーション検出器の厚みを限界まで小さくしたとしても、ベータ線検出器10’におけるガンマ線の検出をゼロにすることはできない。
 一方、本実施形態に係るイメージング装置1は、ガンマ線検出器20を利用することで、参考イメージング装置よりも核種NCB1の分布像を高解像度にて生成することができる。第8実施例に係るイメージング装置1の動作について説明する。既に述べたように、ガンマ線検出器20は、Bタイプ核種NCB1から放出される固有ガンマ線を検出することができる、即ち、エネルギーE1の固有ガンマ線を検出することができる。ガンマ線検出器20にて、Bタイプ核種NCB1から放出された固有ガンマ線が検出されたとき、検出された固有ガンマ線のエネルギーがエネルギーE1であることを示すエネルギー情報E(エネルギー情報E1)を含んだガンマ線検出信号が出力される。つまり、エネルギー情報E1を含むガンマ線検出信号の出力は、Bタイプ核種NCB1からの固有ガンマ線がガンマ線検出器20にて検出されたことを示している。
 図20は、第8実施例に係る画像再構成の概念図である。像800は、撮像対象TG内におけるBタイプ核種NCB1の実際の分布を示している。第8実施例では、固有ガンマ線の検出を伴わないイベントについてのイベントデータの集まりであるイベントデータ群810と、エネルギーE1の固有ガンマ線の検出を伴うイベントについてのイベントデータの集まりであるイベントデータ群811と、が得られる。
 エネルギーE1の固有ガンマ線の検出を伴うイベントについてのイベントデータとは、“1”の同時計測判別信号を含むイベントデータであって、且つ、位置情報(x,y)及びエネルギー情報E1を含むイベントデータである。
 固有ガンマ線の検出を伴わないイベントについてのイベントデータとは、固有ガンマ線がガンマ線検出器20にて検出されないイベントについてのイベントデータであって、典型的には“0”の同時計測判別信号を含むイベントデータであり、従って、位置情報(x,y)を含むがエネルギー情報Eは含まない。第1実施例で述べたように、対消滅等の影響により固有ガンマ線以外のガンマ線がガンマ線検出器20に入射することもあり、その場合には、イベントデータ中の同時計測判別信号が“1”を示すこともあるが、当該イベントデータ中のエネルギー情報EがエネルギーE1を示していないならば、当該イベントデータは、固有ガンマ線の検出を伴わないイベントについてのイベントデータに分類される。
 イベントデータ群811中の各イベントデータの位置情報(x,y)は、撮像対象TG中のBタイプ核種NCB1の存在位置のX軸成分及びY軸成分の検出値を表している。このため、画像再構成部32は、イベントデータ群811中の各イベントデータの位置情報(x,y)をXY平面上にヒストグラミングすることにより、撮像対象TG内におけるBタイプ核種NCB1の分布像821を生成することができる。
 つまり、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングとガンマ線検出器20によるガンマ線検出タイミングとの時間差が所定時間TTH以下であって(即ち、同時計測判別信号が“1”を示していて)且つエネルギー情報E1が得られているときの位置情報(x,y)を用いて、Bタイプ核種NCB1の分布像821を生成することができる。換言すれば、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングから所定時間TTH内にガンマ線検出器20にてエネルギーE1の固有ガンマ線が検出されていることを示すイベントデータ中の位置情報(x,y)を用いて、Bタイプ核種NCB1の分布像821を生成することができる。
 第8実施例では、固有ガンマ線を放出する核種はBタイプ核種NCB1のみであるため、各イベントにおいて、固有ガンマ線は検出されるか検出されないかの2パターンに限定される。このため、ガンマ線検出器20は、固有ガンマ線が検出されたか否かのみを検出すれば足る。そして、イメージング処理部30は、ベータ線検出器10によるベータ線検出タイミングとガンマ線検出器20によるガンマ線検出タイミングとの時間差が所定時間TTH以下であって(即ち、同時計測判別信号が“1”を示していて)且つ固有ガンマ線が検出されたイベントの位置情報(x,y)を用いて、Bタイプ核種NCB1の分布像821を生成できる。
 イベントデータ群810中の各イベントデータの位置情報(x,y)は、撮像対象TG中のBタイプ核種NCB1の存在位置のX軸成分及びY軸成分の検出値を表している。但し、イベントデータ群810中の各イベントデータの位置情報(x,y)には、ベータ線検出器10によるベータ線の検出に基づく位置情報に加えてベータ線検出器10によるガンマ線(固有ガンマ線を含む)の検出に基づく位置情報が含まれる。故に、画像再構成部32にてイベントデータ群810中の各イベントデータの位置情報(x,y)をXY平面上にヒストグラミングしたとき、ベータ線検出器10におけるベータ線の検出とガンマ線の検出を重ね合わせたような画像820が生成されることとなる。画像820は、イベントデータ群811に基づいて生成されるBタイプ核種NCB1の分布像821よりも広い範囲に広がったような画像であり、図19(b)の画像IMβγに類するものである。
 一方、分布像821の元となるイベントデータは、ベータ線と固有ガンマ線とが同時計測されたイベントについてのデータであるので、当該イベントにおいてベータ線検出器10にガンマ線(固有ガンマ線を含む)は入射していない。故に当該イベントについての位置情報から生成された画像(即ち分布像821)は、ベータ線の分布のみを抜き出したものであると言え、高い解像度を有する。このように、イメージング装置1によれば、単一種類のBタイプ核種を高い解像度にてイメージングできる。
 本実施形態に係るイメージング装置1に関し、上述の第7実施例まででは、イメージングに使用する核種からのベータ線を確実に止めるのに必要十分な厚さ(必要な最低限の厚さ;一般的には数ミリメートル以内)のベータ線検出器10を用意することで、殆どのガンマ線がベータ線検出器10を突き抜けると想定し、ベータ線検出器10からのベータ線検出信号にガンマ線の検出に基づく位置情報が含まれていないと仮定している。第8実施例に示した内容は、ベータ線検出器10でのガンマ線の検出が無視できない程度のノイズとなるようなケースにおいて有益であると言える。勿論、ベータ線検出器10でのガンマ線の検出がゼロになることは無いので、程度の多少はあれど任意のケースにおいて、本実施形態に係るイメージング装置1は、参考イメージング装置よりもBタイプ核種の分布像の高解像度化に寄与すると言える。
 尚、ベータ線検出10にてベータ線だけでなくガンマ線も検出されうることを考慮すれば、ベータ線検出器10に関して、以下のように表現することもできる。
 ベータ線検出器10は、撮像対象TG中の放射性核種から放出される放射線を受けることによって当該放射線を検出し、自身に対する放射線の入射を検知する度に、ベータ線検出信号を出力する。ベータ線検出器10にて検出される放射線は主としてベータ線であるが、ベータ線検出器10にて検出される放射線がガンマ線であることをも考慮すると、ベータ線検出器10からのベータ線検出信号を放射線検出信号と読み替えることもでき、当該放射線検出信号は、ベータ線検出器10に対して放射線の入射があったことを示すと共にベータ線検出器10における放射線の検出位置を示す位置情報を含む。また、ベータ線検出器10にて検出される放射線がガンマ線であることをも考慮すると、上述のベータ線検出タイミングをベータ線検出器10による放射線検出タイミングと読み替えることもできる。
<<第9実施例>>
 第9実施例を説明する。イメージング処理部30において、ベータ線検出信号及びガンマ線検出信号から対象核種の分布像を形成する処理は、基本的にハードウェアとソフトウェアの組み合わせにて実現されるが、イメージング処理部30にて実現される機能の一部(例えば、ベータ線検出器20の出力信号から位置情報を導出する処理、同時計測判別の処理)は、ハードウェアで実現される場合もあるし、ソフトウェアで実現される場合もある。特に、同時計測判別をソフトウェアで実現する際は、イベントデータを形成する検出器10の出力信号のデータ及び検出器20の出力信号のデータの夫々に対し共通のタイムスタンプを付加してメモリ(不図示)上に蓄積しておき、蓄積したデータについてタイムスタンプをもとにソフトウェア上で同時計測判別を行う。
 特定の機能をソフトウェアにて実現する場合、その特定の機能をプログラムとして記述しておき、該プログラムをプログラム実行装置(例えば、イメージング処理部30を構成するマイクロコンピュータ)上で実行することによって、その特定の機能を実現するようにしてもよい。上記プログラムは任意の記録媒体に記憶及び固定されうる。上記プログラムを記憶及び固定する記録媒体はイメージング処理部30と異なる機器(サーバ機器等)に搭載又は接続されても良い。
<<第10実施例>>
 第10実施例を説明する。イメージング装置1においては、Aタイプ核種として、Aタイプの放射性崩壊を行う任意の核種を利用することができ、Bタイプ核種として、Bタイプの放射性崩壊を行う任意の核種を利用することができる。
 例えば、正のベータ崩壊により陽電子の放出を行うAタイプ核種として、11C、13N、18F、の内の何れか1つ以上の核種を、撮像対象TGに含めるべきAタイプ核種とすることができる。但し、ここで挙げた核種は、正のベータ崩壊後、100%の確率で娘核の基底状態に遷移する代表的なものであり、正のベータ崩壊により、高い確率で娘核の励起状態を経由せずに娘核の基底状態に遷移する他の核種についても、正のベータ崩壊を行うAタイプ核種とみなして利用可能である。
 例えば、負のベータ崩壊(電子の放出を伴うベータ崩壊)を行うAタイプ核種として、32P、33P、42Ar、66Ni、90Y、106Ru、118Cd、121Sn、143Pr、209Pb、及び、210Pb、の内の何れか1つ以上の核種を、撮像対象TGに含めるべきAタイプ核種とすることができる。但し、ここで挙げた核種は、負のベータ崩壊後、100%の確率で娘核の基底状態に遷移する代表的なものであり、負のベータ崩壊により、高い確率で娘核の励起状態を経由せずに娘核の基底状態に遷移する他の核種についても、負のベータ崩壊を行うAタイプ核種とみなして利用可能である。
 例えば、正のベータ崩壊により陽電子の放出を行うBタイプ核種としては、14O、22Na、34mCl、38K、44Sc、48V、52Mn、52mMn、52Fe、60Cu、72As、76Br、82mRb、94mTc、104mAg、110mIn、及び、124I、の内の何れか1つ以上の核種を、撮像対象TGに含めるべきBタイプ核種とすることができる。但し、ここで挙げた核種は、固有ガンマ線の放出率が9割前後と高い代表的なものであり、放出率がこれらよりも低い他の核種についても利用可能である。
 例えば、負のベータ崩壊(電子の放出を伴うベータ崩壊)を行うBタイプ核種として、24Na、27Mg、38S、43K、46Sc、47Ca、48Sc、55Mn、60Co、72Nn、72Ga、78Ge、82Br、92Sr、95Nb、96Nb、103Ru、130I、132I、134Cs、135Xe、148mPm、160Tb、192Ir、198Au、及び、203Hg、の内の何れか1つ以上の核種を、撮像対象TGに含めるべきBタイプ核種とすることができる。但し、ここで挙げた核種は、固有ガンマ線の放出率が9割前後と高い代表的なものであり、放出率がこれらよりも低い他の核種についても利用可能である。
 特許文献1(特許第5526435号公報)に示された方法は、陽電子放出に伴って生じる対消滅を利用する方法であるため、正のベータ崩壊を行う核種しか利用することができない。これに対し、本実施形態に係るイメージング装置1は、ベータ線を直接イメージングする方法を採用しているため、正のベータ崩壊を行う核種だけでなく、負のベータ崩壊を行う核種も利用可能である。これは、非常に大きなメリットをもたらす。
 例えば、ナトリウムの放射性同位体である24Naやカルシウムの放射性同位体である47Caは、負のベータ崩壊を行うBタイプ核種であるが故に、特許文献1の方法では利用できないが、本実施形態のイメージング装置1では利用できる。周知の如く、ナトリウムやカルシウムは生命活動維持のために重要な役割を担っている。従って、イメージング装置1において、24Na及び47CaをBタイプ核種として用いて区別してイメージングするようにすれば、撮像対象TGとしての生体内のナトリウムとカルシウムの挙動を同時に知ることが可能となる。或いは、47Caを含む複数種類の放射性核種を含んだ生体を撮像対象TGとしてイメージングすれば、生体内のカルシウムの挙動を他の核種と区別して知ることが可能となる(ナトリウムについても同様)。
 尚、放射性核種の中には、親核から娘核の基底状態に遷移する過程において、固有ガンマ線を放出しないことが支配的であるものの、小さい確率(例えば1%未満の確率)で固有ガンマ線を放出する放射性核種もある。このような放射性核種は、厳密にはBタイプ核種と言えるが、そのような放射性核種を、イメージング装置1では、Aタイプ核種として取り扱うことも可能である。
 小さい確率で固有ガンマ線を放出することがあるものの、イメージング装置1においてAタイプ核種として取り扱われる放射性核種を、便宜上、疑似Aタイプ核種と称する。疑似Aタイプ核種から放出され得る固有ガンマ線は、イメージング処理部30での画像の再構成処理において無視される。例えば、図14に対応する第3実施例において、疑似Aタイプ核種がAタイプ核種NCA1として利用される場合、或るイベントについてのガンマ線検出信号中のエネルギー情報EがエネルギーE1又はE2を示していないならば、同時計測判別信号の如何に依らず、当該イベントを固有ガンマ線の検出を伴わないイベントとして取り扱って当該イベントのイベントデータをイベントデータ群410に分類すれば良い(但し、疑似Aタイプ核種から放出され得る固有ガンマ線のエネルギーは、エネルギーE1及びE2と相違し、使用するガンマ線検出器により弁別可能であるものとする)。
 また、Bタイプの放射性崩壊を行う核種によっては、ベータ崩壊後に励起状態に遷移するが、その励起状態の寿命として非常に長い半減期を持つもの(核異性体又はアイソマーと呼ばれる)もある。励起状態の半減期としては、1ナノ秒以下のものから何年にも及ぶものもあるが、同時計測のタイムウインドウ以上に寿命が長い場合は、陽電子又は電子の放出と固有ガンマ線の放出とが連続しているものとして検出されないため、固有ガンマ線を核種の識別に使用することはできない。
 尚、PET装置では、正のベータ崩壊により陽電子を放出する核種をイメージングに利用するが、正のベータ崩壊が起きた際には必ず陽電子が放出されるわけではなく、陽電子放出は電子捕獲(原子核の周りを回る電子を陽子が吸収して中性子に変化する現象)との競合過程となる。即ち、陽電子放出及び電子捕獲の2つの過程が競合して、どちらか一方の過程のみが起こる。陽電子放出が起こる確率及び電子捕獲が起こる確率は、核種によって決まっているが、電子捕獲のみによって正のベータ崩壊が起きる核種や、電子捕獲による正のベータ崩壊が支配的である核種が存在するため、全ての正のベータ崩壊核種がPETイメージングに利用できるわけではない。
 電子捕獲による正のベータ崩壊が起きた場合、陽電子は放出されないが、捕獲された電子の軌道に空きができるため、その軌道に外側の軌道から電子が遷移する際、余ったエネルギーが特性X線として正のベータ崩壊核種から放出される。この特性X線を吸収するのに十分な厚さの検出器を使用すれば、特性X線がイメージング装置1に入射することにより、原理的にはベータ線と同様にイメージングが可能である。
 本発明では、PETと異なり負のベータ崩壊核種も使用可能であるが、上記の如く、正のベータ崩壊において陽電子放出と電子捕獲は競合過程にある。従って、核種によって陽電子を放出する確率が決まっており、陽電子を出さない核種もある。一方で、負のベータ崩壊では、必ず電子を放出するため、負のベータ崩壊核種の方が、本発明で利用可能な核種の数が多い。
 また、負のベータ崩壊核種は中性子が過剰な核種であるため原子炉などによる中性子照射によっても製造が可能であり、陽子が過剰な核種であって、サイクロトロン等の加速器による荷電粒子の照射を必要とする正のベータ崩壊核種よりも製造コストがかからない場合が多い。
<<第11実施例>>
 第11実施例を説明する。上述してきたように、固有ガンマ線は核種に固有のエネルギーを持っているため、これを計測することにより核種を同定することが可能となるが、これは、固有ガンマ線の全エネルギーがガンマ線検出器20にて検出されることを前提としている。しかし、実際には、ガンマ線検出器20では、全エネルギー吸収の他に、コンプトン散乱などによる部分エネルギー吸収も起こる。特にコンプトン散乱により生じる部分エネルギー吸収に関しては散乱角度により連続的に吸収エネルギーが変化するため、ガンマ線検出器20の検出結果から得られるエネルギースペクトルにおいて、全エネルギー吸収によるピークよりエネルギー的に低い部分に、部分エネルギー吸収に対応する成分が連続成分として現れる。このため、相対的に高いエネルギーを有する第1の固有ガンマ線を放出する第1のBタイプ核種と相対的に低いエネルギーを有する第2の固有ガンマ線を放出する第2のBタイプ核種とが撮像対象TGに含まれているとき、第1の固有ガンマ線の部分エネルギー吸収によるスペクトル成分と、第2の固有ガンマ線の全エネルギー吸収によるスペクトル成分とが重なり合う。また、それらのスペクトル成分と、陽電子の対消滅ガンマ線(対消滅により発生した511keVのガンマ線)の全エネルギー吸収又は部分エネルギー吸収によるスペクトル成分とが重なり合うこともある。
 ガンマ線検出器20の検出結果から得られるエネルギースペクトルにおいて、全エネルギー吸収に対応するピークの前後に存在する部分エネルギー吸収による連続成分が、ピークに重なった連続成分と同じであるとみなした上で、ピークの部分から当該連続成分を差し引けば、これらの重なりによるバックグランドを除去可能である。また、ガンマ線検出器20としてエネルギー分解能が十分に高いガンマ線検出器を使い、全エネルギー吸収イベント(ガンマ線検出器20での全エネルギー吸収を伴うイベント)によるピーク幅を狭くすることにより、上記重なりによるバックグランドを減らすことができる。
 また、部分エネルギー吸収イベント(ガンマ線検出器20での部分エネルギー吸収を伴うイベント)については、以下のようにすることで固有ガンマ線の検出効率を向上させることが可能である。説明の具体化のため、エネルギーE1の固有ガンマ線を放出する上述のBタイプ核種NCB1に注目する。核種NCB1からの固有ガンマ線がガンマ線検出器20にて部分エネルギー吸収されたとき、エネルギーE1未満のエネルギーがガンマ線検出器20にて検出されることになるが、エネルギーE1より小さな所定のエネルギー閾値以上のエネルギーが検出されたならば、固有ガンマ線が検出されたと判断するようにする。このエネルギー閾値は、ガンマ線検出器20のエネルギー分解能を考慮した上で、固有ガンマ線と異なるガンマ線による成分が混じり込まない範囲において設定されると良い。これにより、固有ガンマ線の部分エネルギー吸収が生じたときにおいても、固有ガンマ線が検出されることになるため、固有ガンマ線の検出効率の向上が見込める。
 上記の通り、部分エネルギー吸収イベントについては、他の核種の全エネルギー吸収イベントと重なることによる識別精度の低下や、エネルギー閾値の設定による検出効率の向上が考えられるが、上述の各実施例においては、説明の具体化及び簡略化のため、全エネルギー吸収によるイベントを用いた識別方法を前提とした説明を設けている。
 尚、解像度向上のための手法として、シンチレーション光(可視光)の屈折又は反射を利用したレンズ等により拡大画像を得る方法が知られていると上述したが、その方法に属する方式の内、イベント毎にシンチレーション光のカウントが可能な方式であれば、本発明に応用できる。
<<発明の考察>>
 本発明について考察する。
 本発明の一側面に係るイメージング装置W1は(特に例えば図5及び図12参照)、ベータ線二次元イメージング装置であって、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に第1固有ガンマ線を放出する第1核種(例えばBタイプ核種NCB1)と、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に前記第1固有ガンマ線とは異なるエネルギーを有する第2固有ガンマ線を放出する第2核種(例えばBタイプ核種NCB2)とを含んだ撮像対象(TG)から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を出力するベータ線検出器(10)と、ガンマ線を検出する検出器であって、前記第1固有ガンマ線及び前記第2固有ガンマ線を区別して検出するガンマ線検出器(20)と、前記ベータ線検出器によるベータ線検出タイミング及び前記ガンマ線検出器によるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出器にて前記第1固有ガンマ線及び第2固有ガンマ線の何れかが検出されたかと、に基づいて、前記第1核種の分布像(例えば分布像221)及び前記第2核種の分布像(例えば分布像222)を区別して生成可能なイメージング処理部(30)と、を備えたことを特徴とする。
 これにより、複数種類の核種の分布を、核種を識別可能な態様で一度にイメージングすることが可能となる。また、ベータ線を直接イメージングする方法を採用しているため、正のベータ崩壊を行う核種だけでなく、負のベータ崩壊を行う核種も利用可能であり、有益である。
 具体的には例えば、イメージング装置W1において、前記イメージング処理部は、前記ベータ線検出タイミングから所定時間内に前記第1固有ガンマ線が検出されたときの前記位置情報を用いて前記第1核種の分布像を生成し、前記ベータ線検出タイミングから前記所定時間内に前記第2固有ガンマ線が検出されたときの前記位置情報を用いて前記第2核種の分布像を生成すると良い。
 本発明の一側面に係るイメージング装置W2は(特に例えば図5及び図13参照)、ベータ線二次元イメージング装置であって、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に所定エネルギーの固有ガンマ線を放出する第1核種(例えばBタイプ核種NCB1)と、ベータ崩壊によるベータ線の放出により娘核の基底状態に遷移する第2核種(例えばAタイプ核種NCA1)とを含んだ撮像対象(TG)から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を出力するベータ線検出器(10)と、前記固有ガンマ線を含むガンマ線を検出するガンマ線検出器(20)と、前記ベータ線検出器によるベータ線検出タイミング及び前記ガンマ線検出器によるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出器にて前記固有ガンマ線が検出されたかと、に基づいて、前記第1核種の分布像(例えば分布像321)及び前記第2核種の分布に応じた他の分布像(例えば分布像320又はAタイプ分布像)を区別して生成可能なイメージング処理部(30)と、を備えたことを特徴とする。
 これにより、複数種類の核種の分布を、核種を識別可能な態様で一度にイメージングすることが可能となる。また、ベータ線を直接イメージングする方法を採用しているため、正のベータ崩壊を行う核種だけでなく、負のベータ崩壊を行う核種も利用可能であり、有益である。
 具体的には例えば、イメージング装置W2において、前記イメージング処理部は、前記ベータ線検出タイミングから所定時間内に前記固有ガンマ線が検出されたときの前記位置情報を用いて前記第1核種の分布像を生成し、それ以外の前記位置情報を用いて、前記第1核種の分布と前記第2核種の分布を重ねあわせた分布像を生成すると良い。
 また、イメージング装置W1及びW2の夫々において、負のベータ崩壊が生じる核種を、前記第1核種及び前記第2核種の一方又は双方に用いることが可能である。
 本発明の一側面に係るイメージング装置W3は(特に例えば図5及び図18参照)、ベータ線二次元イメージング装置であって、正のベータ崩壊によりベータ線として陽電子を放出する第1核種(例えば正のベータ崩壊核種NCP1)及び負のベータ崩壊によりベータ線として電子を放出する第2核種(例えば負のベータ崩壊核種NCN1)を含んだ撮像対象(TG)から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を出力するベータ線検出器(10)と、前記第1核種からのベータ線としての陽電子と前記ベータ線検出器内の電子との対消滅により発生する対消滅ガンマ線を含むガンマ線を検出するガンマ線検出器(20)と、前記ベータ線検出器によるベータ線検出タイミング及び前記ガンマ線検出器によるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出器にて前記対消滅ガンマ線が検出されたかと、に基づいて、前記第1核種の分布像(例えば分布像721)及び前記第2核種の分布に応じた他の分布像(例えば分布像720又は負のベータ崩壊核種分布像)を区別して生成可能なイメージング処理部(30)と、を備えたことを特徴とする。
 これにより、複数種類の核種の分布を、核種を識別可能な態様で一度にイメージングすることが可能となる。また、ベータ線を直接イメージングする方法を採用しているため、正のベータ崩壊を行う核種だけでなく、負のベータ崩壊を行う核種も利用可能であり、有益である。
 具体的には例えば、イメージング装置W3において、前記イメージング処理部は、前記ベータ線検出タイミングから所定時間内に前記対消滅ガンマ線が検出されたときの前記位置情報を用いて前記第1核種の分布像を生成し、それ以外の前記位置情報を用いて、前記第1核種の分布と前記第2核種の分布を重ねあわせた分布像を生成すると良い。
 本発明の一側面に係るイメージング方法W1Aは、ベータ線二次元イメージング方法であって、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に第1固有ガンマ線を放出する第1核種と、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に前記第1固有ガンマ線とは異なるエネルギーを有する第2固有ガンマ線を放出する第2核種とを含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を得るベータ線検出ステップと、ガンマ線を検出するステップであって、前記第1固有ガンマ線及び前記第2固有ガンマ線を区別して検出するガンマ線検出ステップと、前記ベータ線検出ステップにおけるベータ線検出タイミング及び前記ガンマ線検出ステップにおけるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出ステップにて前記第1固有ガンマ線及び第2固有ガンマ線の何れかが検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布像を区別して生成するイメージング処理ステップと、を備えたことを特徴とする。
 本発明の一側面に係るイメージング方法W2Aは、ベータ線二次元イメージング方法であって、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に所定エネルギーの固有ガンマ線を放出する第1核種と、ベータ崩壊によるベータ線の放出により娘核の基底状態に遷移する第2核種とを含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を得るベータ線検出ステップと、前記固有ガンマ線を含むガンマ線を検出するガンマ線検出ステップと、前記ベータ線検出ステップにおけるベータ線検出タイミング及び前記ガンマ線検出ステップにおけるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出ステップにて前記固有ガンマ線が検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布に応じた他の分布像を区別して生成するイメージング処理ステップと、を備えたことを特徴とする。
 本発明の一側面に係るイメージング方法W3Aは、正のベータ崩壊によりベータ線として陽電子を放出する第1核種及び負のベータ崩壊によりベータ線として電子を放出する第2核種を含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を得るベータ線検出ステップと、前記第1核種からのベータ線としての陽電子と前記ベータ線検出器内の電子との対消滅により発生する対消滅ガンマ線を含むガンマ線を検出するガンマ線検出ステップと、前記ベータ線検出ステップにおけるベータ線検出タイミング及び前記ガンマ線検出ステップにおけるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出ステップにて前記対消滅ガンマ線が検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布に応じた他の分布像を区別して生成可能なイメージング処理ステップと、を備えたことを特徴とする。
 本発明の一側面に係るイメージング装置W4は(特に例えば図5及び図20参照)、ベータ線二次元イメージング装置であって、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に固有ガンマ線を放出する核種(例えばBタイプ核種NCB1)を含んだ撮像対象(TG)から、前記ベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を出力するベータ線検出器(10)と、前記固有ガンマ線を含むガンマ線を検出するガンマ線検出器(20)と、前記ベータ線検出器によるベータ線検出タイミング及び前記ガンマ線検出器によるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出器にて前記固有ガンマ線が検出されたかと、に基づいて、前記核種の分布像(例えば分布像821)を生成可能なイメージング処理部(30)と、を備えたことを特徴とする。
 ベータ線検出器は、撮像対象中の核種からのベータ線を検出するためのものであるが、核種からのガンマ線がベータ線検出器にて検出されることもある。但し、ベータ線検出器はベータ線とガンマ線を区別して検出することができないため、ベータ線検出信号中の位置情報に、ベータ線の検出位置を示す位置情報に加えてガンマ線の検出位置を示す位置情報が互いに区別不能な状態で混在するようになる。ガンマ線の検出位置を示す位置情報がノイズとして機能し、仮にガンマ線検出器が設けられていなかったならば、生成される核種の分布像の解像度を劣化させる。ガンマ線はベータ線よりも飛程が長いため、実際の核種の分布よりも広がった分布の像が生成されるためである(図19(b)参照)。
 他方、上述の如く形成されたイメージング装置W4では、ガンマ線検出器にて固有ガンマ線が検出されたときの位置情報に基づいて核種の分布像を生成することができる。ガンマ線検出器にて固有ガンマ線が検出されたとき、ベータ線検出器に固有ガンマ線は入射していないので、ガンマ線検出器にて固有ガンマ線が検出されたときの位置情報はベータ線の検出位置情報に限定され、その検出位置情報に基づき核種の分布像を生成することで、上記ノイズが排除された高解像度の分布像を生成することができる。
 具体的には例えば、イメージング装置W4において、前記イメージング処理部は、前記ベータ線検出タイミングから所定時間内に前記固有ガンマ線が検出されたときの前記位置情報を用いて前記核種の分布像を生成すると良い。
 本発明の一側面に係るイメージング方法W4Aは、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に固有ガンマ線を放出する核種を含んだ撮像対象から、前記ベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を得るベータ線検出ステップと、前記固有ガンマ線を含むガンマ線を検出するガンマ線検出ステップと、前記ベータ線検出ステップにおけるベータ線検出タイミング及び前記ガンマ線検出ステップにおけるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出ステップにて前記固有ガンマ線が検出されたかと、に基づいて、前記核種の分布像を生成可能なイメージングステップと、を備えたことを特徴とする。
 尚、本発明の実施形態は、特許請求の範囲に示された技術的思想の範囲内において、適宜、種々の変更が可能である。以上の実施形態は、あくまでも、本発明の実施形態の例であって、本発明ないし各構成要件の用語の意義は、以上の実施形態に記載されたものに制限されるものではない。上述の説明文中に示した具体的な数値は、単なる例示であって、当然の如く、それらを様々な数値に変更することができる。
  1 イメージング装置
 10 ベータ線検出器
 20 ガンマ線検出器
 30 イメージング処理部
 31 同時計測判別部
 32 画像再構成部
 TG 撮像対象

Claims (13)

  1.  ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に第1固有ガンマ線を放出する第1核種と、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に前記第1固有ガンマ線とは異なるエネルギーを有する第2固有ガンマ線を放出する第2核種とを含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を出力するベータ線検出器と、
     ガンマ線を検出する検出器であって、前記第1固有ガンマ線及び前記第2固有ガンマ線を区別して検出するガンマ線検出器と、
     前記ベータ線検出器によるベータ線検出タイミング及び前記ガンマ線検出器によるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出器にて前記第1固有ガンマ線及び第2固有ガンマ線の何れかが検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布像を区別して生成可能なイメージング処理部と、を備えた
    ことを特徴とするベータ線二次元イメージング装置。
  2.  前記イメージング処理部は、前記ベータ線検出タイミングから所定時間内に前記第1固有ガンマ線が検出されたときの前記位置情報を用いて前記第1核種の分布像を生成し、前記ベータ線検出タイミングから前記所定時間内に前記第2固有ガンマ線が検出されたときの前記位置情報を用いて前記第2核種の分布像を生成する
    ことを特徴とする請求項1に記載のベータ線二次元イメージング装置。
  3.  ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に所定エネルギーの固有ガンマ線を放出する第1核種と、ベータ崩壊によるベータ線の放出により娘核の基底状態に遷移する第2核種とを含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を出力するベータ線検出器と、
     前記固有ガンマ線を含むガンマ線を検出するガンマ線検出器と、
     前記ベータ線検出器によるベータ線検出タイミング及び前記ガンマ線検出器によるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出器にて前記固有ガンマ線が検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布に応じた他の分布像を区別して生成可能なイメージング処理部と、を備えた
    ことを特徴とするベータ線二次元イメージング装置。
  4.  前記イメージング処理部は、前記ベータ線検出タイミングから所定時間内に前記固有ガンマ線が検出されたときの前記位置情報を用いて前記第1核種の分布像を生成し、それ以外の前記位置情報を用いて、前記第1核種の分布と前記第2核種の分布を重ねあわせた分布像を生成する
    ことを特徴とする請求項3に記載のベータ線二次元イメージング装置。
  5.  負のベータ崩壊が生じる核種を、前記第1核種及び前記第2核種の一方又は双方に用いることが可能である
    ことを特徴とする請求項1~4の何れかに記載のベータ線二次元イメージング装置。
  6.  正のベータ崩壊によりベータ線として陽電子を放出する第1核種及び負のベータ崩壊によりベータ線として電子を放出する第2核種を含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を出力するベータ線検出器と、
     前記第1核種からのベータ線としての陽電子と前記ベータ線検出器内の電子との対消滅により発生する対消滅ガンマ線を含むガンマ線を検出するガンマ線検出器と、
     前記ベータ線検出器によるベータ線検出タイミング及び前記ガンマ線検出器によるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出器にて前記対消滅ガンマ線が検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布に応じた他の分布像を区別して生成可能なイメージング処理部と、を備えた
    ことを特徴とするベータ線二次元イメージング装置。
  7.  前記イメージング処理部は、前記ベータ線検出タイミングから所定時間内に前記対消滅ガンマ線が検出されたときの前記位置情報を用いて前記第1核種の分布像を生成し、それ以外の前記位置情報を用いて、前記第1核種の分布と前記第2核種の分布を重ねあわせた分布像を生成する
    ことを特徴とする請求項6に記載のベータ線二次元イメージング装置。
  8.  ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に第1固有ガンマ線を放出する第1核種と、ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に前記第1固有ガンマ線とは異なるエネルギーを有する第2固有ガンマ線を放出する第2核種とを含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を得るベータ線検出ステップと、
     ガンマ線を検出するステップであって、前記第1固有ガンマ線及び前記第2固有ガンマ線を区別して検出するガンマ線検出ステップと、
     前記ベータ線検出ステップにおけるベータ線検出タイミング及び前記ガンマ線検出ステップにおけるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出ステップにて前記第1固有ガンマ線及び第2固有ガンマ線の何れかが検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布像を区別して生成可能なイメージング処理ステップと、を備えた
    ことを特徴とするベータ線二次元イメージング方法。
  9.  ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に所定エネルギーの固有ガンマ線を放出する第1核種と、ベータ崩壊によるベータ線の放出により娘核の基底状態に遷移する第2核種とを含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を得るベータ線検出ステップと、
     前記固有ガンマ線を含むガンマ線を検出するガンマ線検出ステップと、
     前記ベータ線検出ステップにおけるベータ線検出タイミング及び前記ガンマ線検出ステップにおけるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出ステップにて前記固有ガンマ線が検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布に応じた他の分布像を区別して生成可能なイメージング処理ステップと、を備えた
    ことを特徴とするベータ線二次元イメージング方法。
  10.  正のベータ崩壊によりベータ線として陽電子を放出する第1核種及び負のベータ崩壊によりベータ線として電子を放出する第2核種を含んだ撮像対象から、前記第1核種又は前記第2核種に基づくベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を得るベータ線検出ステップと、
     前記第1核種からのベータ線としての陽電子と前記ベータ線検出器内の電子との対消滅により発生する対消滅ガンマ線を含むガンマ線を検出するガンマ線検出ステップと、
     前記ベータ線検出ステップにおけるベータ線検出タイミング及び前記ガンマ線検出ステップにおけるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出ステップにて前記対消滅ガンマ線が検出されたかと、に基づいて、前記第1核種の分布像及び前記第2核種の分布に応じた他の分布像を区別して生成可能なイメージング処理ステップと、を備えた
    ことを特徴とするベータ線二次元イメージング方法。
  11.  ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に固有ガンマ線を放出する核種を含んだ撮像対象から、前記ベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を出力するベータ線検出器と、
     前記固有ガンマ線を含むガンマ線を検出するガンマ線検出器と、
     前記ベータ線検出器によるベータ線検出タイミング及び前記ガンマ線検出器によるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出器にて前記固有ガンマ線が検出されたかと、に基づいて、前記核種の分布像を生成可能なイメージング処理部と、を備えた
    ことを特徴とするベータ線二次元イメージング装置。
  12.  前記イメージング処理部は、前記ベータ線検出タイミングから所定時間内に前記固有ガンマ線が検出されたときの前記位置情報を用いて前記核種の分布像を生成する
    ことを特徴とする請求項11に記載のベータ線二次元イメージング装置。
  13.  ベータ崩壊によって娘核の励起状態となり、当該ベータ崩壊によるベータ線の放出に続けて娘核の基底状態に遷移する際に固有ガンマ線を放出する核種を含んだ撮像対象から、前記ベータ線を受けることで当該ベータ線を検出し、当該ベータ線の検出位置を二次元で示す位置情報を含んだベータ線検出信号を得るベータ線検出ステップと、
     前記固有ガンマ線を含むガンマ線を検出するガンマ線検出ステップと、
     前記ベータ線検出ステップにおけるベータ線検出タイミング及び前記ガンマ線検出ステップにおけるガンマ線検出タイミングと、前記ベータ線検出信号に含まれる前記位置情報と、前記ガンマ線検出ステップにて前記固有ガンマ線が検出されたかと、に基づいて、前記核種の分布像を生成可能なイメージングステップと、を備えた
    ことを特徴とするベータ線二次元イメージング方法。
PCT/JP2018/007009 2017-03-01 2018-02-26 ベータ線二次元イメージング装置及び方法 WO2018159548A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US16/490,276 US11259768B2 (en) 2017-03-01 2018-02-26 Apparatus and method for beta-emission two-dimensional imaging
EP18760568.8A EP3591438A4 (en) 2017-03-01 2018-02-26 DEVICE AND METHOD FOR BIDIMENSIONAL BETA-RAY IMAGING
JP2019502985A JP7100902B2 (ja) 2017-03-01 2018-02-26 ベータ線二次元イメージング装置及び方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2017-038201 2017-03-01
JP2017038201 2017-03-01

Publications (1)

Publication Number Publication Date
WO2018159548A1 true WO2018159548A1 (ja) 2018-09-07

Family

ID=63370048

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2018/007009 WO2018159548A1 (ja) 2017-03-01 2018-02-26 ベータ線二次元イメージング装置及び方法

Country Status (4)

Country Link
US (1) US11259768B2 (ja)
EP (1) EP3591438A4 (ja)
JP (1) JP7100902B2 (ja)
WO (1) WO2018159548A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020180816A (ja) * 2019-04-24 2020-11-05 国立大学法人千葉大学 反射型断層撮影装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5526435B1 (ja) 1970-11-12 1980-07-12
JP2008209336A (ja) * 2007-02-28 2008-09-11 Hitachi Ltd 核医学診断装置
JP2013033001A (ja) * 2011-08-03 2013-02-14 Institute Of Physical & Chemical Research Pet装置およびそのイメージング方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL83847A (en) * 1987-09-09 1991-06-10 Yehoshua Kalish Scintillation spectrometer
JP4064009B2 (ja) * 1999-07-30 2008-03-19 株式会社東芝 線種弁別型放射線検出装置
US8565860B2 (en) * 2000-08-21 2013-10-22 Biosensors International Group, Ltd. Radioactive emission detector equipped with a position tracking system
ES2258932B1 (es) * 2005-02-25 2007-11-16 Universidad De Barcelona Sensor radioquimico para fluidos.
WO2007005442A2 (en) * 2005-07-01 2007-01-11 Warburton William K Detection of coincident radiations in a single transducer by pulse shape analysis
US7683334B2 (en) * 2007-08-07 2010-03-23 The State Of Oregon Acting By And Through The State Board Of Higher Education On Behalf Of Oregon State University Simultaneous beta and gamma spectroscopy
WO2009093396A1 (ja) 2008-01-25 2009-07-30 Nec Corporation 非線形歪補償回路、送信回路、及び非線形歪補償方法
US8039810B2 (en) * 2008-04-22 2011-10-18 Canberra Industries, Inc. Scintillation-based continuous monitor for beta-emitting radionuclides in a liquid medium
WO2012048274A2 (en) * 2010-10-07 2012-04-12 H. Lee Moffitt Cancer Center & Research Institute Method and apparatus for detection of radioactive isotopes
US8916829B2 (en) * 2011-08-18 2014-12-23 Savannah River Nuclear Solutions, Llc System and method for assaying a radionuclide
US9018586B2 (en) * 2012-06-27 2015-04-28 Batelle Energy Alliance, Llc Apparatuses for large area radiation detection and related method
US9823364B2 (en) 2014-02-03 2017-11-21 Arizona Board Of Regents On Behalf Of The University Of Arizona Beta and alpha emission tomography for three-dimensional autoradiography

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5526435B1 (ja) 1970-11-12 1980-07-12
JP2008209336A (ja) * 2007-02-28 2008-09-11 Hitachi Ltd 核医学診断装置
JP2013033001A (ja) * 2011-08-03 2013-02-14 Institute Of Physical & Chemical Research Pet装置およびそのイメージング方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP3591438A4

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020180816A (ja) * 2019-04-24 2020-11-05 国立大学法人千葉大学 反射型断層撮影装置

Also Published As

Publication number Publication date
JPWO2018159548A1 (ja) 2020-01-09
US20200022667A1 (en) 2020-01-23
US11259768B2 (en) 2022-03-01
JP7100902B2 (ja) 2022-07-14
EP3591438A4 (en) 2020-12-30
EP3591438A1 (en) 2020-01-08

Similar Documents

Publication Publication Date Title
Hajima et al. Proposal of nondestructive radionuclide assay using a high-flux gamma-ray source and nuclear resonance fluorescence
Moskal et al. Novel detector systems for the positron emission tomography
US7253416B2 (en) Radiation detector and detection method, and medical diagnostic apparatus employing same
US8466418B2 (en) Gamma ray detector, radiation diagnostic device, tomography device, and method of analyzing tomography device
JP5526435B2 (ja) Pet装置およびそのイメージング方法
JP5957099B2 (ja) 単一エネルギガンマ線源による、同位体識別、分析およびイメージングのための二重同位体ノッチ観測機
Hausladen et al. Induced-fission imaging of nuclear material
US20230288584A1 (en) Device for the simultaneous deection, identifcation, quantification and/or localization of gamma radiation and neutron sources
CN115485585A (zh) 在重合检测的情况下基于前向散射伽马射线的pet扫描仪中的透射成像
WO2018159548A1 (ja) ベータ線二次元イメージング装置及び方法
Kandzia et al. Development of a liquid scintillator based active fission target for FIPPS
Lang et al. Sub-millimeter nuclear medical imaging with high sensitivity in positron emission tomography using beta-gamma coincidences
Bacon et al. Muon induced fission neutrons in coincidence with muon tomography
US11733183B2 (en) Imaging method and system
Han Advances in Gamma-Ray Imaging with Intensified Quantum-Imaging Detectors
Choghadi et al. Collimator-based Coincidence Imaging for Double Photon Emitting Nuclides
Finocchiaro From nuclear physics to applications: detectors for beam handling, medical diagnostics and radioactive waste monitoring
Belcari et al. Positron emission tomography: its 65 years and beyond
Rego Optimization study of the performance of a LaBr3 monolithic scintillator in a Compton camera system
Borghi From detectors towards systems: enabling clinical TOF-PET with monolithic scintillators
McAreavey Development of the DEPICT system for gamma-ray imaging in molecular radiotherapy
Kalantan Feasibility of using the GRI+ Compton camera for proton therapy range verification
Clarke Development of a scintillation-based CMOS Quantitative Autoradiography imager for Safeguard Applications
Newton Entangled polarisation correlations of annihilation gamma and their applications to PET Imaging
CN116134342A (zh) 信号处理系统、正电子放射断层摄影装置和正电子放射断层摄影方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18760568

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2019502985

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2018760568

Country of ref document: EP

Effective date: 20191001