WO2016203822A1 - 放射線検出素子の感度補正方法および放射線断層撮影装置 - Google Patents

放射線検出素子の感度補正方法および放射線断層撮影装置 Download PDF

Info

Publication number
WO2016203822A1
WO2016203822A1 PCT/JP2016/061172 JP2016061172W WO2016203822A1 WO 2016203822 A1 WO2016203822 A1 WO 2016203822A1 JP 2016061172 W JP2016061172 W JP 2016061172W WO 2016203822 A1 WO2016203822 A1 WO 2016203822A1
Authority
WO
WIPO (PCT)
Prior art keywords
sensitivity
radiation
sensitivity coefficient
detector
source
Prior art date
Application number
PCT/JP2016/061172
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 US15/736,168 priority Critical patent/US10126445B2/en
Priority to JP2017524688A priority patent/JP6607256B2/ja
Priority to CN201680035775.XA priority patent/CN107710018B/zh
Publication of WO2016203822A1 publication Critical patent/WO2016203822A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T7/00—Details of radiation-measuring instruments
    • G01T7/005—Details of radiation-measuring instruments calibration techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/17—Circuit arrangements not adapted to a particular type of detector
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/17—Circuit arrangements not adapted to a particular type of detector
    • G01T1/172—Circuit arrangements not adapted to a particular type of detector with coincidence circuit arrangements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
    • G01T1/2914—Measurement of spatial distribution of radiation
    • G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)

Definitions

  • the present invention relates to a method for correcting the sensitivity of a radiation detection element that detects radiation generated from a radiopharmaceutical in a subject, and a radiation tomography apparatus that performs radiation tomography by detecting radiation generated from a radiopharmaceutical in a subject. .
  • sensitivity correction Coefficients for sensitivity correction (hereinafter referred to as “sensitivity coefficients”) are calculated from actual data collected using a calibration source, and there are two methods for calculating sensitivity coefficients.
  • the first is a direct method for obtaining sensitivity coefficients for all detector pairs.
  • the direct method the sensitivity coefficients of all detector pairs can be obtained directly, but when there are many detectors, the number of detector pairs becomes enormous due to the square of the number of detectors. Therefore, the statistical accuracy per pair is lowered, and the number of sensitivity coefficients is enormous.
  • NC uivj ⁇ ui ⁇ ⁇ vj ⁇ b uvk ⁇ d uvrk ⁇ g uvr ⁇ f uv (1)
  • ⁇ ui and ⁇ vj are detector specific sensitivities
  • d uvrk is a crystal interference factor
  • f uv is a ring pair sensitivity
  • b uvk is a block shape factor
  • g uvr is a radial geometric factor.
  • K is the relative crystal position in the block
  • r is the radial position.
  • g uvr , d uvrk and the like are geometrically determined elements (hereinafter referred to as “geometric system factors”). Note that ⁇ ui , ⁇ vj , and b uvk are non-geometric factors that change over time.
  • the sensitivity correction method for each element is represented by the product of sensitivity coefficients, so the number of sensitivity coefficients is smaller than that of the direct method. Further, since the factors obtained by decomposing the sensitivity coefficient are added using geometric symmetry, the statistical accuracy can be improved. However, since the sensitivity coefficient is obtained indirectly by the product of factors, it is an approximate numerical value. In recent PET apparatuses, the element-specific sensitivity correction method has begun to be used for the purpose of high resolution, easy calibration, and improved statistical accuracy.
  • FIG.1 is a calibration flow for calculating geometric factors using a conventional low-scattering calibration source
  • Figure.2 is a calibration flow for calculating non-geometric factors using a conventional uniform cylindrical calibration source. It is. The geometric factor that does not change with time is calculated by the low scattering calibration source, and the non-geometric factor that changes with time is calculated by the data of the uniform cylindrical calibration source corrected by the coefficient of the geometric factor. And the sensitivity coefficient of a specific detector pair is calculated
  • Such a conventional elemental sensitivity correction method has the following problems. That is, there is a problem that the multi-modality compatible flexible PET apparatus does not support the conventional elemental sensitivity correction method.
  • the conventional sensitivity correction method for each element is established when the detector group has a circular arrangement and the positional relationship is fixed.
  • the multi-modality-compatible flexible PET apparatus is composed of two detector heads 2A and 2B in which arc-shaped detector groups are arranged, and each detection is performed as shown by a two-dot chain line in FIG.
  • the instrument heads 2A and 2B are close to each other.
  • the positional relationship becomes variable and becomes elliptical, and the assumption of the elemental sensitivity correction method is broken and cannot be established (see AB in non-proximity and A′B in close proximity shown in FIG. 3).
  • the sensitivity coefficient is unchanged if the detector pair is in the same detector head, but the distance between the detector heads facing each other changes, and ⁇ with respect to the detector changes. As the line incidence angle is different, the sensitivity of the geometric factor at the opposite detector head changes. As a result, artifacts are generated in the reconstructed image.
  • the present invention has been made in view of such circumstances, and a reconstructed image free from artifacts can be obtained even when the arrangement of the plurality of radiation detection elements constituting the detector is variable and not circular. It is an object of the present invention to provide a sensitivity correction method for a radiation detection element and a radiation tomography apparatus.
  • the present invention has the following configuration. That is, the sensitivity correction method according to the present invention is a sensitivity correction method for a radiation detection element that detects radiation generated from a radiopharmaceutical in a subject, and is based on coincidence count data collected by a plurality of the radiation detection elements.
  • the sensitivity coefficient is obtained in the first sensitivity coefficient calculation step, and after the first sensitivity coefficient calculation step, the change in the arrangement of the plurality of radiation detection elements.
  • the sensitivity coefficient is obtained in the second sensitivity coefficient calculation step before, and after the second sensitivity coefficient calculation step, the sensitivity coefficient is obtained in the third sensitivity coefficient calculation step in a state where the arrangement of the plurality of radiation detection elements is changed.
  • the sensitivity coefficient is obtained separately for each sensitivity factor in the first sensitivity coefficient calculation step. Even if the arrangement of the plurality of radiation detection elements does not change, the non-geometric factor changes with time.
  • the sensitivity coefficient of the non-geometric system factor is obtained in the second sensitivity coefficient calculation step based on the collected coincidence data.
  • the third sensitivity coefficient calculation step the sensitivity coefficient derived from the geometric arrangement is obtained based on the coincidence count data collected in a state where the arrangement of the plurality of radiation detection elements is changed.
  • the sensitivity coefficient when the arrangement of a plurality of radiation detection elements is changed is the arrangement for a combination whose positional relationship remains unchanged regardless of the arrangement of the radiation detection elements (for example, the same detector head). It is represented by the product of the sensitivity coefficient of the geometric factor obtained before the change and the sensitivity coefficient of the non-geometric factor. For combinations that are geometrically affected by the positional relationship of the radiation detection elements (for example, detector heads facing each other), the sensitivity factor of the geometric factor obtained after the change in arrangement and the non-geometric system It is expressed as the product of the factor's sensitivity coefficient.
  • a correction coefficient from the ratio of coincidence count data collected at the same time using the same calibration source before and after the change of the arrangement of the plurality of radiation detection elements.
  • the same calibration source the latter: the first: before and after the change of the arrangement of the plurality of radiation detection elements.
  • the difference between the former before the change of the arrangement and the latter after the change of the arrangement is obtained by multiplying the obtained correction coefficient by the right side of the above equation (1). Can also be corrected.
  • the coincidence data collected after the change in the arrangement of the plurality of radiation detection elements is a combination that is geometrically influenced by at least the positional relationship of the radiation detection elements. That is, based on the coincidence data collected after the change of the arrangement, sensitivity coefficient data for a combination (for example, detector heads facing each other) that is geometrically influenced by the positional relationship of the radiation detection elements can be obtained.
  • the calibration rod radiation source can be regarded as a calibration ring radiation source by rotating along the plurality of radiation detection elements. Moreover, you may use the ring source for calibration provided along the some radiation detection element. In addition, when using a calibration plate source that can rotate along a plurality of radiation detection elements, the calibration plate source may be rotated along the plurality of radiation detection elements.
  • the rotated calibration rod source that is, the calibration rod source regarded as the calibration ring source
  • the calibration ring source or the rotated calibration plate source.
  • the sensitivity coefficient is obtained separately for each sensitivity factor.
  • the sensitivity coefficient of the non-geometric factor that changes with time can be obtained.
  • each detector head composed of a plurality of radiation detection elements
  • each detector head is configured such that the distance between the detector heads can be changed, the detector heads are brought close to each other. Therefore, when the arrangement of the plurality of radiation detection elements changes, it is preferable to obtain the sensitivity coefficient using a calibration radiation source as described below. That is, when trying to bring the detector heads close to each other, a calibration rod-shaped radiation source that can be rotated along a plurality of radiation detection elements and a calibration ring radiation source provided along the plurality of radiation detection elements respectively detect each detector head. May interfere with the instrument head.
  • a calibration flat radiation source installed at the center of each detector head is used, or a calibration rod-shaped radiation source capable of reciprocating along the central axis of each detector head is used.
  • the calibration rod-shaped radiation source is reciprocated along the central axis of each detector head to make it a calibration flat plate radiation source.
  • the calibration flat plate source or the reciprocating calibration rod source that is, the calibration rod source assumed to be a calibration flat plate source
  • the sensitivity coefficient derived from the geometry is obtained.
  • a radiation tomography apparatus is a radiation tomography apparatus that performs radiation tomography by detecting radiation generated from a radiopharmaceutical in a subject, and includes a detector head comprising a plurality of radiation detection elements.
  • a plurality of detector heads are configured such that the distances between the detector heads can be changed, and each of the detector heads includes a calculation unit that performs calculation processing.
  • the calculation units are collected by the plurality of radiation detection elements.
  • the first sensitivity coefficient calculation step for obtaining the sensitivity coefficient separately for each sensitivity factor, and after the first sensitivity coefficient calculation step and before the change of the arrangement of the plurality of radiation detection elements
  • a second sensitivity coefficient calculation step for obtaining a sensitivity coefficient of the non-geometric factor, and an arrangement of the plurality of radiation detection elements is changed after the second sensitivity coefficient calculation step.
  • a third sensitivity coefficient calculation step for obtaining a sensitivity coefficient derived from the geometric arrangement based on the coincidence data collected in a state where the first sensitivity coefficient is calculated, the second sensitivity coefficient calculation step, and the second sensitivity coefficient calculation step.
  • the sensitivity correction of the radiation detection element is performed using the sensitivity coefficient obtained in the three sensitivity coefficient calculation step.
  • the first sensitivity coefficient calculation step, the second sensitivity coefficient calculation step, and the third sensitivity coefficient calculation step as described in the sensitivity correction method according to the present invention are performed.
  • Computation means for performing is provided.
  • the sensitivity correction of the radiation detection element is performed using the sensitivity coefficients obtained in the first sensitivity coefficient calculation step, the second sensitivity coefficient calculation step, and the third sensitivity coefficient calculation step.
  • the first sensitivity (determining the sensitivity coefficient for each sensitivity factor based on the coincidence data collected by the plurality of radiation detection elements) Coefficient calculation step (after the first sensitivity coefficient calculation step and before the change of the arrangement of the plurality of radiation detection elements, the sensitivity coefficient of the non-geometric system factor is obtained based on the collected coincidence data)
  • a sensitivity coefficient calculation step and (after the second sensitivity coefficient calculation step, a sensitivity coefficient derived from the geometric arrangement is obtained based on coincidence data collected in a state where the arrangement of the plurality of radiation detection elements has changed)
  • FIG. 1 It is the schematic perspective view and block diagram which show one embodiment of the detector head of the flexible PET apparatus for multi-modality which concerns on an Example.
  • FIG. 1 It is a schematic perspective view of a gamma ray detector.
  • A is the schematic which looked at the detector head from the depth direction
  • (b) is the schematic which looked from the horizontal direction orthogonal to a depth direction.
  • (a) is a schematic diagram of a ring radiation source when not approaching
  • (b) is a schematic diagram of a cylindrical radiation source when not approaching
  • (c) is It is a schematic diagram of the flat plate source at the time of proximity
  • FIG. 1 is a schematic perspective view and a block diagram showing an embodiment of a detector head of a flexible PET apparatus for multi-modality according to an embodiment
  • FIG. 2 is a schematic perspective view of a ⁇ -ray detector.
  • (A) is the schematic which looked at the detector head from the depth direction
  • FIG.3 (b) is the schematic which looked from the horizontal direction orthogonal to a depth direction.
  • a multi-modality compatible PET apparatus (hereinafter simply referred to as “PET apparatus”) 1 includes detector heads 2A and 2B.
  • PET apparatus As shown in FIG. 2, 18 detector blocks 31 ( ⁇ -ray detector 3) in which detector crystals arranged vertically and horizontally in 16 ⁇ 16 are stacked in four layers are arranged in an arc shape in the depth direction z (FIG. 1).
  • a detector head is composed of a total of 18 ⁇ 3 detector blocks 31 arranged in three. A specific configuration of the ⁇ -ray detector 3 other than the detector block 31 will be described later with reference to FIG.
  • the multi-modality flexible PET apparatus (PET apparatus) 1 corresponds to the radiation tomography apparatus in the present invention
  • the detector heads 2A and 2B correspond to the detector head in the present invention
  • the ⁇ -ray detector 3 It corresponds to the radiation detection element in the invention.
  • the detector heads 2A and 2B are configured to face each other.
  • the detector heads 2A and 2B are usually located on a circumference of a 48-gon.
  • the total number of detector blocks 31 ( ⁇ -ray detectors 3) in the cross section of the two detector heads 2A and 2B is 48 squares for a total of 36 pieces (18 x 2 in an arc shape).
  • Each gap region 4 exists.
  • the detector heads 2A and 2B can be brought close to each other with the detector heads 2A and 2B facing each other.
  • each detector head 2A, 2B is comprised so that the distance between each detector head 2A, 2B can be changed.
  • the PET apparatus 1 includes a coincidence counting circuit 5 and an arithmetic circuit 6. Specific functions of the arithmetic circuit 6 will be described later with reference to FIG.
  • the arithmetic circuit 6 corresponds to the arithmetic means in this invention.
  • the coincidence counting circuit 5 checks the position of the detector block 31 (see FIG. 2) and the incident timing of ⁇ rays, and only when ⁇ rays are simultaneously incident on the two detector blocks 31 on both sides of the subject.
  • the sent image information is determined as appropriate data (simultaneous counting data).
  • the coincidence circuit 5 rejects.
  • the coincidence counting circuit 5 detects that ⁇ rays are simultaneously observed (that is, coincidence counting) by the two ⁇ ray detectors 3 based on the above-described electrical signal.
  • Image information (simultaneous counting data) collected by the coincidence circuit 5 is sent to the arithmetic circuit 6.
  • the ⁇ -ray detector 3 includes a detector block 31, a light guide 32 optically coupled to the detector block 31, and an optically coupled to the light guide 32. And a photomultiplier tube (PMT: Photo Multiplier Tube) 33.
  • PMT Photo Multiplier Tube
  • Each detector crystal constituting the detector block 31 emits light with the incidence of ⁇ rays, thereby converting the ⁇ rays into light. By this conversion, the detector crystal detects ⁇ rays.
  • the light emitted from the detector crystal is sufficiently diffused by the detector block 31 and input to a photomultiplier tube (hereinafter simply referred to as “PMT”) 33 through the light guide 32.
  • the PMT 33 multiplies the light converted by the detector block 31 and converts it into an electrical signal.
  • the electric signal is sent to the coincidence counting circuit 5 (see FIG. 1) as image information (simultaneous counting data) as described above.
  • the ⁇ -ray detector 3 is a DOI detector composed of a plurality of layers in the depth direction.
  • a four-layer DOI detector is illustrated, but the number of layers is not particularly limited as long as it is plural. Also, the number of detector crystals arranged vertically and horizontally is not particularly limited as long as it is plural.
  • the DOI detector is constructed by laminating each detector crystal in the depth direction of the radiation, and the interaction depth (DOI: Depth of Interaction) direction and lateral direction (incident) Coordinate information with the direction parallel to the surface) is obtained by calculating the center of gravity.
  • DOI Depth of Interaction
  • lateral direction incident
  • the number of DOI detector layers is the number of detector crystal layers stacked in the depth direction.
  • FIGS. 4 (a) to 4 (c) are schematic diagrams of detector arrangements and types of calibration radiation sources
  • FIG. 4 (a) is a schematic diagram of a ring radiation source when not in proximity.
  • 4 (b) is a schematic diagram of a cylindrical radiation source at the time of non-proximity
  • FIG. 4 (c) is a schematic diagram of a flat-plate radiation source at the time of proximity
  • FIG. 5 is a sensitivity correction for each element according to the embodiment.
  • FIG. 6 is a flowchart showing the flow of the processing procedure of the method
  • FIG. 6 shows a state where a rod-shaped radiation source that can rotate along the ⁇ -ray detector of the detector head is rotated along the ⁇ -ray detector of the detector head.
  • FIG. 4 (a) is a schematic diagram of a ring radiation source when not in proximity.
  • 4 (b) is a schematic diagram of a cylindrical radiation source at the time of non-proximity
  • FIG. 4 (c) is a schematic diagram of
  • FIGS. 4 (a) code S R calibration ring-ray source in (hereinafter, abbreviated as “ring-ray source”) is a
  • FIG. 4 (b) code S C is a calibration cylinder ray source in (hereinafter it is merely referred to as “cylinder-ray source”)
  • FIG. 4 (c) code S P calibration flat plates ray source in (hereinafter, it is abbreviated as "flat-ray source”).
  • a calibration rod-shaped radiation source (rotating along the ⁇ -ray detector 3 (see FIG. 2) (ie, radiation detection element) of the detector heads 2A and 2B ( hereinafter simply using abbreviated as “rod-shaped radiation source”) S L, shown in FIG.
  • a rod-shaped radiation source S L detector heads 2A by rotating along the ⁇ -ray detectors 3 2B Mitateru to ring-ray source S R.
  • the thickness of the ring-ray source S R is about 5 mm
  • the diameter of the cylindrical ray source S C is about 15cm
  • a thickness of the flat plate-ray source S P is approximately 5 mm
  • the diameter of the rod-shaped radiation source S L is about 15cm.
  • the non-illustrated in FIG. 4 rod-shaped radiation source during non-proximity shown in (a) S L ring-ray source likened to the ring is rotated (see Figure 6) S R, 4 (b) A uniform cylindrical source S C installed at the center of the apparatus (detector heads 2A, 2B) when approaching, and a flat plate installed at the center of the apparatus (detector heads 2A, 2B) shown in FIG. It provides three types of sources S P, to collect data, respectively.
  • the calibration radiation source is made of acrylic resin, and a radioactive solution of uniform concentration is sealed inside.
  • the long-axis block shape factor b uvk corresponds to “Transaxial Block Profile Factor” and the crystal interference factor d uvrk corresponds to “Crystal Interference Factor”.
  • Block shape factor b uvk corresponds to “Axial Block Profile Factor”
  • ring pair sensitivity f uv corresponds to “Ring Pair Factor”
  • detector specific sensitivities ⁇ ui , ⁇ vj correspond to “Crystal Efficiency Factor”
  • the crystal interference factor d Plate at the time of the flat plate source corresponds to “Crystal Interference (Plate) Factor”
  • the ring pair sensitivity f Plate at the time of the flat plate source corresponds to “Ring Pair (Plate) Factor”.
  • Step S1 As shown in the flowchart of FIG. 5, first use coincidence data collected by the ring-ray source S R shown in FIG. 4 (a), crystal interfering agents geometric system factors that change over time d uvrk ( “Crystal Calculated as “Interference Factor”. It is preferable to apply path length normalization (denoted as “Apply”) before calculating the crystal interference factor d uvrk (Crystal Interference Factor) of this geometric factor.
  • the path length normalization is to normalize the coincidence count data on a predetermined path using the length that the path passes through the radiation source. By performing path length normalization, a difference in length that intersects the radiation source between predetermined paths is corrected in a pseudo manner.
  • path length standardization refer to Japanese Patent Application Laid-Open No. 2008-170329.
  • crystal interference factor d uvrk Crystal Interference Factor
  • the coefficient of crystal interference factor d uvrk (Crystal Interference Factor) obtained as a result is output as a table.
  • LOR Line Of Response
  • This step S1 corresponds to the first sensitivity coefficient calculation step in the present invention.
  • Step S2 Next, using the coincidence data crystalline interferer d uvrk (Crystal Interference Factor) cylindrical ray source shown in FIG. 4 (b) performing the correction of the geometric-based factors by S C, non-geometric-based time varying
  • the detector's intrinsic sensitivities ⁇ ui and ⁇ vj (expressed as “Crystal Efficiency Factor”) are sequentially calculated.
  • the block shape factor b uvk ( indicated by “Axial Block Profile Factor”) of the body axis section (cross section of the body axis) and the block shape factor b uvk (indicated by “Transaxial Block Profile Factor”) of the long axis calculated, it performs correction on coincidence data of the flat radiation source S P shown in FIG. 4 (c) in step S3.
  • the ring pair sensitivity f uv (expressed as “Ring Pair Factor”) is calculated. Each factor obtained by calculation is output as a table.
  • This step S2 corresponds to the second sensitivity coefficient calculation step in the present invention.
  • Step S3 Finally, using the coincidence counting data of the flat plate-ray source S P shown in FIG. 4 (c), crystal interfering agents d Plate when the flat plate-ray source affected by changes in the geometry of the detector by adjacent (Represented by “Crystal Interference (Plate) Factor”) and ring pair sensitivity f Plate (represented by “Ring Pair (Plate) Factor”) at the time of a flat plate source are calculated.
  • the crystal interference factor d uvrk Crystal Interference Factor
  • Sensitivity f uv Sensitivity f uv (Ring Pair Factor) is not calculated.
  • the combination of the detector heads facing each other is geometrically influenced by the positional relationship of the ⁇ -ray detector 3 (radiation detection element), so that the crystal interference factor d Plate (Crystal Interference ( Plate) Factor) and the ring pair sensitivity f Plate (Ring Pair (Plate) Factor) at the time of a flat plate source.
  • step S3 path length normalization is preferably applied before calculating these factors.
  • the solid angle correction is performed with the coordinates at the time of proximity to the crystal interference factor d Plate (Crystal Interference (Plate) Factor) at the time of the flat plate source.
  • Output the calculated crystal interference factor d Plate (Crystal Interference (Plate) Factor) and flat pair source coefficient f Plate (Ring Pair (Plate) Factor) coefficients for a flat plate source as a table.
  • This step S3 corresponds to the third sensitivity coefficient calculation step in the present invention.
  • Non-Patent Document 1 for specific calculation of each factor.
  • the flexible PET apparatus uses a fan-sum method because there is a gap between the detector heads.
  • Patent Document 2 International Publication No. 2010/013356 Please refer to.
  • N near is the coincidence count data (count value) collected using the flat plate source at the time of proximity
  • N normal is the coincidence count collected using the plate source at the time of non-proximity.
  • Data (count value) L near is the detector distance when close
  • L normal is the detector distance when not close.
  • a sensitivity coefficient is obtained in the first sensitivity coefficient calculation step (step S1), and after the first sensitivity coefficient calculation step (step S1), a plurality of radiation detection elements.
  • the sensitivity coefficient is obtained in the second sensitivity coefficient calculation step (step S2) before the change of the arrangement of the ( ⁇ -ray detector in the present embodiment) (in the case of non-proximity in the present embodiment), and the second sensitivity coefficient calculation step (step After S2), the sensitivity coefficient is obtained in the third sensitivity coefficient calculation step (step S3) in a state where the arrangement of the plurality of radiation detection elements has changed (in this embodiment, the detector heads 2A and 2B are close to each other).
  • the sensitivity coefficient is obtained separately for each sensitivity factor in the first sensitivity coefficient calculation step (step S1). Even if the arrangement of the plurality of radiation detection elements ( ⁇ -ray detectors) does not change, the non-geometric factor changes with time. Therefore, the non-geometric system is used in the second sensitivity coefficient calculation step (step S2) based on the collected coincidence data before the change of the arrangement of the plurality of radiation detection elements ( ⁇ -ray detectors) (when not in proximity). Find the sensitivity factor of the factor. Then, in the third sensitivity coefficient calculation step (step S3), the geometric arrangement is based on the coincidence count data collected in a state where the arrangement of the plurality of radiation detection elements is changed (the detector heads 2A and 2B are close to each other).
  • the radiation detection element (step S1), the second sensitivity coefficient calculation step (step S2), and the third sensitivity coefficient calculation step (step S3) are used to obtain the radiation detection element (
  • a reconstructed image is obtained based on the data after sensitivity correction.
  • a reconstructed image free from artifacts can be obtained even when the arrangement of the plurality of radiation detection elements ( ⁇ -ray detectors) constituting the detector is variable and not circular.
  • the sensitivity coefficient when the arrangement of the plurality of radiation detection elements is changed has the same positional relationship regardless of the arrangement of the radiation detection elements ( ⁇ -ray detectors).
  • the sensitivity factor of the geometric factor obtained before the change of arrangement in the present embodiment in the case of non-proximity
  • the sensitivity of the non-geometric factor It is expressed as a product with a coefficient.
  • the PET apparatus 1 determines whether the detector pair is the same detector head or the detector heads facing each other. )) Determines. If the opposing detector head, switching to the sensitivity coefficients in the nearby obtained from flat plate-ray source S P shown in FIG. 4 (c).
  • CPU central processing unit
  • the same calibration source (same flat plate source in this embodiment) is used before and after the change of the arrangement of the plurality of radiation detection elements ( ⁇ -ray detectors) (when approaching or not approaching) as in the above formula (2). It is preferable to obtain the correction coefficient ⁇ from the ratio of the coincidence count data collected at the same time (the left side of the above equation (2)).
  • a plurality of radiation detection elements are used.
  • the same calibration source (the latter: the second calibration source, in the present embodiment, the flat plate source) is used before and after the change of arrangement (at the time of proximity or non-proximity)
  • the obtained correction coefficient ⁇ is used as the former.
  • the coincidence data collected after the change of the arrangement of multiple radiation detection elements ( ⁇ -ray detectors) (at the time of proximity) has a geometrical influence due to the positional relationship of at least the radiation detection elements ( ⁇ -ray detectors). It is a combination that receives. That is, combinations that are geometrically influenced by the positional relationship of the radiation detection elements ( ⁇ -ray detectors) based on the coincidence data collected after the change in arrangement (at the time of proximity) (in the present embodiment, face each other) Sensitivity coefficient data for the detector head) can be obtained.
  • the first sensitivity coefficient calculation step (step S1) it is preferable to obtain the sensitivity coefficient using a calibration source as described below. That is, as in the present embodiment, when a rod-shaped radiation source S L (see FIG. 6) that can rotate along a plurality of radiation detection elements ( ⁇ -ray detectors) is used, a plurality of rod-shaped radiation sources S L are provided. Mitateru by rotating along a radiation detection element (gamma ray detector) in a ring-ray source S R shown in Figure 4 (a).
  • the first sensitivity coefficient calculating step (Step S1) (in this example ⁇ -rays) radiation from the rotating rod-like radiation source S L (i.e.
  • the rod-shaped radiation source is likened to a ring-ray source S R S L) for detecting the Based on the coincidence data collected in this way, the sensitivity coefficient is obtained separately for each sensitivity factor.
  • the second sensitivity coefficient calculation step (step S2) it is preferable to obtain the sensitivity coefficient using a calibration source as described below. That is, in the second sensitivity coefficient calculating step (step S2), and prior to the change of the arrangement of a plurality of radiation detecting elements (gamma ray detector) (non-proximity state), model and from the cylinder-ray source S C to the actual subject Based on the coincidence data collected by detecting the radiation ( ⁇ rays), the sensitivity coefficient of the non-geometric factor that changes with time can be obtained.
  • a plurality of detector heads (two in this embodiment) including a plurality of radiation detection elements ( ⁇ -ray detectors) are provided, and each detector head is configured such that the distance between the detector heads can be changed.
  • the sensitivity coefficient can be obtained using the following calibration radiation source. preferable. That is, when each detector head is brought close to each other, it is provided along a rod-shaped radiation source that can be rotated along a plurality of radiation detection elements ( ⁇ -ray detectors) or a plurality of radiation detection elements ( ⁇ -ray detectors). There is a possibility that the generated ring source interferes with each detector head.
  • the third sensitivity coefficient calculating step (step S3) and in a state of being close to the detector head one another, based on the coincidence count data collected by detecting radiation (gamma rays) from the flat plate-ray source S P Then, the sensitivity coefficient derived from the geometric arrangement is obtained.
  • Computation means (arithmetic circuit 6 in FIG. 1) that performs the sensitivity coefficient calculation step (step S2) and the third sensitivity coefficient calculation step (step S3) is provided.
  • it was obtained in the first sensitivity coefficient calculation step (step S1), the second sensitivity coefficient calculation step (step S2), and the third sensitivity coefficient calculation step (step S3).
  • a cylindrical radiation source with a uniform radioactivity concentration of 300 mm in diameter is installed at the center of the device (detector heads 2A, 2B) in the proximity of the device (detector heads 2A, 2B) to a diameter of 778mm to 250mm. Collected.
  • the reconstruction result incorporating the conventional elemental sensitivity correction method using the geometric factor of the ring source and the non-geometric system factor of the cylindrical source uses the non-proximate geometric factor. Artifacts occur. Non-geometric factors for cylindrical sources, events in the same detector head are derived from ring sources, events in opposite detector heads are derived from flat sources, and correction factor ⁇ It was confirmed that artifacts disappeared and the image quality improved as a result of the included reconstruction.
  • the present invention is not limited to the above embodiment, and can be modified as follows.
  • the radiation is ⁇ rays, but may be radiation such as ⁇ rays and ⁇ rays.
  • the DOI detector is used.
  • the DOI detector may be applied to a radiation detector that does not distinguish the depth direction.
  • each detector head 2A and 2B has been described, but the number of detector heads is not limited to two. If there is a structure in which a plurality of detector heads each including a plurality of radiation detection elements (gamma ray detectors in the embodiment) are provided and the distance between the detector heads can be changed, each detector head is 3 There may be more than one.
  • the rod-shaped radiation source S L is formed by using the rod-shaped radiation source S L (see FIG. 6) that can rotate along a plurality of radiation detection elements ( ⁇ -ray detectors in the embodiment).
  • S R by rotating along a plurality of radiation detecting elements (gamma ray detector), S R and itself ring-ray source shown in FIG. 4 (a)
  • S P may be rotated along a flat plate-ray source S P a plurality of radiation detecting elements (gamma ray detector) shown in Figure 7 (a). That may be a ring-ray source S R provided along the plurality of radiation detecting elements (gamma ray detector), as shown in Figure 4 (a).
  • FIG. 7 in the case of using a rotatable flat-ray source S P along a plurality of radiation detecting elements as shown in (a) (gamma ray detector)
  • the flat-ray source S P a plurality of radiation detecting elements You may rotate along (gamma ray detector).
  • FIG. 4 A rod-shaped source S L that can reciprocate along the center axis of each detector head shown in FIG. That is, in the case of using a reciprocable rod-shaped radiation source S L along the central axis of each detector head as shown in FIG. 7 (b), along the rod-shaped radiation source S L to the central axis of each detector head Mitateru flat plate-ray source S P by reciprocating Te.
  • the arrangement before the change of the arrangement of the plurality of radiation detection elements is not close, and after the change of the arrangement, the arrangement is close. In other words, before the change of arrangement, it may be close, and after the change of arrangement, it may be non-close. Therefore, the first sensitivity coefficient calculating step (step S1), the based on the coincidence count data collected at proximity with ring-ray source S R oval shown in FIG. 7 (c), the sensitivity coefficient by the sensitivity factor Ask separately. Then, in the third sensitivity coefficient calculation step (step S3), the sensitivity coefficient derived from the geometric arrangement is obtained based on the coincidence data collected at the time of non-proximity.
  • the shape of the calibration radiation source in the first sensitivity coefficient calculation step (step S1), the second sensitivity coefficient calculation step (step S2), and the third sensitivity coefficient calculation step (step S3) is not particularly limited.
  • a calibration radiation source having a shape suitable for the detector arrangement may be used.

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Nuclear Medicine (AREA)
  • Measurement Of Radiation (AREA)

Abstract

回転した校正用棒状線源からのγ線を検出することで収集された同時計数データに基づいて、感度係数を感度要因別に分けて求める第１感度係数算出工程（ステップＳ１）および複数のγ線検出器の配置が変化した状態で収集された同時計数データに基づいて、幾何学配置由来の感度係数を求める第３感度係数算出工程（ステップＳ３）で求められた感度係数を用いて、γ線検出器の感度補正を行うことにより、感度補正後のデータに基づき再構成画像を得る。その結果、検出器を構成する複数のγ線検出器の配置が可変で円形配置でない場合であってもアーティファクトのない再構成画像を得ることができる。

Description

放射線検出素子の感度補正方法および放射線断層撮影装置
　この発明は、被検体内の放射性薬剤から発生した放射線を検出する放射線検出素子の感度補正方法、および被検体内の放射性薬剤から発生した放射線を検出して放射線断層撮影を行う放射線断層撮影装置に関する。
　放射線断層撮影装置のＰＥＴ(Positron Emission Tomography) において、再構成画像を得る際に検出器間の感度差を校正して全検出器の感度を均質にする必要がある。以下、この校正を「感度補正」と呼ぶ。感度補正のための係数（以下、「感度係数」と呼ぶ）は校正用線源を用いて収集された実データから計算され、感度係数の計算方法は２種類ある。
　１つ目が全検出器ペアについて感度係数を求める直接法である。直接法では、全検出器ペアの感度係数が直接的に得られるが、検出器が多い場合に検出器の数の二乗で検出器ペア数が膨大になる。よって、１ペア当たりの統計精度が低くなり、感度係数の個数が膨大になる。
　次に、感度を変動させる要因ごとに分解し、それらの積を検出器ペアの感度係数とする「要素別感度補正法」である（例えば、特許文献１、２、非特許文献１～３参照）。この要素別感度補正法では、同時計数する検出器リング対を（ｕ，ｖ）、同時計数するリング内検出器ペアを（ｉ，ｊ）としたときに、感度補正係数ＮＣｕｉｖｊを下記（１）式のような要素に分解する。
　ＮＣｕｉｖｊ＝εｕｉ×εｖｊ×ｂｕｖｋ×ｄｕｖｒｋ×ｇｕｖｒ×ｆｕｖ　　　　　　　…（１）
　ただし、上記（１）式中のεｕｉ，εｖｊは検出器固有感度、ｄｕｖｒｋは結晶干渉因子、ｆｕｖはリングペア感度、ｂｕｖｋはブロック形状因子、ｇｕｖｒは動径方向幾何学因子、ｋはブロック内の結晶相対位置、ｒは動径方向位置である。これらの要素のうち、ｇｕｖｒ，ｄｕｖｒｋなどは幾何学的に決まる要素（以下、「幾何学系因子」と呼ぶ）である。なお、εｕｉ，εｖｊ，ｂｕｖｋは経時変化する非幾何学系因子である。
　要素別感度補正法では感度係数の積で表すので、直接法に比べて感度係数の個数は少ない。また、感度係数を分解した各因子は幾何学的な対称性を利用して加算するので、統計精度を上げることができる。しかし、感度係数は因子の積で間接的に得られるので近似的な数値となる。近年のＰＥＴ装置は、高分解能や校正の簡便化および統計精度の改善の目的で要素別感度補正法が利用され始めている。
　従来の要素別感度補正法の校正フローは非特許文献１のFigure.1およびFigure.2に示す通りである。Figure.1は、従来の低散乱校正用線源による幾何学系因子の計算の校正フローであり、Figure.2は、従来の均一円筒校正用線源による非幾何学系因子の計算の校正フローである。低散乱校正用線源により経時変化しない幾何学系因子を計算し、幾何学系因子の係数で補正された均一円筒校正用線源のデータにより経時変化する非幾何学系因子を計算する。そして、上記（１）式により、特定の検出器ペアの感度係数を求める。具体的な校正フローについては、非特許文献１を参照されたい。
特開２００８－１７０３２９号公報 国際公開第２０１０／０１３３５６号
Badawi RD, Marsden PK: Developments in component-based normalization for 3D PET. Phys. Med. Biol. 44, 571-594, 1999 Badawi RD, Lodge MA, Marsden PK: Algorithms for calculating detector efficiency normalization coefficients for true coincidences in 3D PET. Phys. Med. Biol. 43, 189-205, 1998 Badawi RD, Ferreira NC, Kohlmyer SG, Dahlbom M, Marsden PK, Lewellen TK: A comparison of normalization effects on three whole-body cylindrical 3D PET systems. Phys. Med. Biol. 45, 3253-3266, 2000
　しかしながら、このような従来の要素別感度補正法の場合には、次のような問題がある。
　すなわち、マルチモダリティ対応フレキシブルＰＥＴ装置は従来の要素別感度補正法に対応していないという問題がある。従来の要素別感度補正法は検出器群が円形配置で、位置関係が固定である場合に成立する。
　マルチモダリティ対応フレキシブルＰＥＴ装置は、図３に示すように、円弧状の検出器群が配置された２つの検出器ヘッド２Ａ，２Ｂで構成され、図３の二点鎖線に示すように、各検出器ヘッド２Ａ，２Ｂが互いに近接する。その結果、位置関係が可変かつ楕円形状になり、要素別感度補正法の仮定が崩れて成立しなくなる（図３に示した非近接時のＡＢおよび近接時のＡ´Ｂを参照）。
　また、２つの検出器ヘッドが近接した場合、同一の検出器ヘッド内の検出器ペアであれば感度係数は不変であるが、互いに対向する検出器ヘッドの間の距離が変わり、検出器に対するγ線の入射角が異なるので、対向する検出器ヘッドでの幾何学系因子の感度は変化する。その結果、再構成画像にアーティファクトが生じる。
　この発明は、このような事情に鑑みてなされたものであって、検出器を構成する複数の放射線検出素子の配置が可変で円形配置でない場合であってもアーティファクトのない再構成画像が得られる放射線検出素子の感度補正方法および放射線断層撮影装置を提供することを目的とする。
　この発明は、このような目的を達成するために、次のような構成をとる。
　すなわち、この発明に係る感度補正方法は、被検体内の放射性薬剤から発生した放射線を検出する放射線検出素子の感度補正方法であって、複数の前記放射線検出素子で収集された同時計数データに基づいて、感度係数を感度要因別に分けて求める第１感度係数算出工程と、前記第１感度係数算出工程の後で、前記複数の放射線検出素子の配置の変化前に、収集された同時計数データに基づいて、非幾何学系因子の感度係数を求める第２感度係数算出工程と、前記第２感度係数算出工程の後で、前記複数の放射線検出素子の配置が変化した状態で収集された同時計数データに基づいて、幾何学配置由来の感度係数を求める第３感度係数算出工程とを備え、前記第１感度係数算出工程，前記第２感度係数算出工程および前記第３感度係数算出工程で求められた感度係数を用いて、放射線検出素子の感度補正を行うことを特徴とするものである。
　［作用・効果］この発明に係る感度補正方法によれば、先ず、第１感度係数算出工程で感度係数を求めて、第１感度係数算出工程の後で、複数の放射線検出素子の配置の変化前に第２感度係数算出工程で感度係数を求めて、第２感度係数算出工程の後で、複数の放射線検出素子の配置が変化した状態で、第３感度係数算出工程で感度係数を求める。具体的には、第１感度係数算出工程で感度係数を感度要因別に分けて求める。複数の放射線検出素子の配置が変化しなくとも、非幾何学系因子は経時変化する。そこで、複数の放射線検出素子の配置の変化前に、収集された同時計数データに基づいて、第２感度係数算出工程で非幾何学系因子の感度係数を求める。そして、第３感度係数算出工程では、複数の放射線検出素子の配置が変化した状態で収集された同時計数データに基づいて、幾何学配置由来の感度係数を求める。このように、第１感度係数算出工程，第２感度係数算出工程および第３感度係数算出工程で求められた感度係数を用いて、放射線検出素子の感度補正を行うことにより、感度補正後のデータに基づき再構成画像を得る。その結果、検出器を構成する複数の放射線検出素子の配置が可変で円形配置でない場合であってもアーティファクトのない再構成画像を得ることができる。
　具体的な一例は、複数の放射線検出素子の配置が変化した場合の感度係数は、放射線検出素子の配置によらず位置関係が不変の組み合わせ（例えば同一の検出器ヘッド）に対しては、配置変化前に求められた幾何学系因子の感度係数と、非幾何学系因子の感度係数との積で表される。そして、放射線検出素子の位置関係により幾何学的な影響を受ける組み合わせ（例えば互いに対向する検出器ヘッド）に対しては、配置変化後に求められた幾何学系因子の感度係数と、非幾何学系因子の感度係数との積で表される。
　また、複数の放射線検出素子の配置の変化前後で同じ校正用線源を用いて同時間にそれぞれ収集された同時計数データの比から補正係数を求めるのが好ましい。これにより、同じ校正用線源を用いた際に配置の変化により生じる差を補正することができる。ひいては、例えば配置の変化前での別の校正用線源（前者：第１校正用線源）を用いた後に、複数の放射線検出素子の配置の変化前後で同じ校正用線源（後者：第２校正用線源）を用いる場合には、求められた補正係数を前者での上記（１）式の右辺に乗算することで、配置の変化前の前者と配置の変化後の後者との差をも補正することができる。
　また、複数の放射線検出素子の配置の変化後に収集された同時計数データは、少なくとも放射線検出素子の位置関係により幾何学的な影響を受ける組み合わせである。すなわち、配置の変化後に収集された同時計数データに基づいて、放射線検出素子の位置関係により幾何学的な影響を受ける組み合わせ（例えば互いに対向する検出器ヘッド）に対する感度係数データを求めることができる。
　また、第１感度係数算出工程では、下記のような校正用線源を用いて感度係数を求めるのが好ましい。すなわち、複数の放射線検出素子に沿って回転可能な校正用棒状線源を用いる場合には、校正用棒状線源を複数の放射線検出素子に沿って回転させることで校正用リング線源に見立てる。また、複数の放射線検出素子に沿って設けられた校正用リング線源を用いてもよい。また、複数の放射線検出素子に沿って回転可能な校正用平板線源を用いる場合には、校正用平板線源を複数の放射線検出素子に沿って回転させてもよい。このように、第１感度係数算出工程では、回転した校正用棒状線源（すなわち校正用リング線源に見立てられた校正用棒状線源），校正用リング線源または回転した校正用平板線源からの放射線を検出することで収集された同時計数データに基づいて、感度係数を感度要因別に分けて求める。
　また、第２感度係数算出工程では、下記のような校正用線源を用いて感度係数を求めるのが好ましい。すなわち、第２感度係数算出工程では、複数の放射線検出素子の配置の変化前に、実際の被検体に模した校正用円筒線源からの放射線を検出することで収集された同時計数データに基づいて、経時変化する非幾何学系因子の感度係数を求めることができる。
　また、複数の放射線検出素子からなる検出器ヘッドを複数備え、各当該検出器ヘッド間の距離が変更可能に各々の検出器ヘッドが構成されている場合で、各検出器ヘッドを互いに近接させることで複数の放射線検出素子の配置が変化するときには、下記のような校正用線源を用いて感度係数を求めるのが好ましい。すなわち、各検出器ヘッドを互いに近接させようとすると、複数の放射線検出素子に沿って回転可能な校正用棒状線源や複数の放射線検出素子に沿って設けられた校正用リング線源が各検出器ヘッドと干渉する恐れがある。そこで、各検出器ヘッドの中心に設置された校正用平板線源を用いる、あるいは各検出器ヘッドの中心軸に沿って往復可能な校正用棒状線源を用いる。当該校正用棒状線源を用いる場合には、校正用棒状線源を各検出器ヘッドの中心軸に沿って往復させることで校正用平板線源に見立てる。第３感度係数算出工程では、各検出器ヘッドを互いに近接させた状態で、校正用平板線源または往復した校正用棒状線源（すなわち校正用平板線源に見立てられた校正用棒状線源）からの放射線を検出することで収集された同時計数データに基づいて、幾何学配置由来の感度係数を求める。このような校正用平板線源または往復した校正用棒状線源（校正用平板線源に見立てられた校正用棒状線源）を用いることで、各検出器ヘッドを近接させても、校正用線源と各検出器ヘッドとの干渉を防止することができる。
　また、この発明に係る放射線断層撮影装置は、被検体内の放射性薬剤から発生した放射線を検出して放射線断層撮影を行う放射線断層撮影装置であって、複数の放射線検出素子からなる検出器ヘッドを複数備え、各当該検出器ヘッド間の距離が変更可能に各々の検出器ヘッドは構成されており、演算処理を行う演算手段を備え、当該演算手段は、複数の前記放射線検出素子で収集された同時計数データに基づいて、感度係数を感度要因別に分けて求める第１感度係数算出工程と、前記第１感度係数算出工程の後で、前記複数の放射線検出素子の配置の変化前に、収集された同時計数データに基づいて、非幾何学系因子の感度係数を求める第２感度係数算出工程と、前記第２感度係数算出工程の後で、前記複数の放射線検出素子の配置が変化した状態で収集された同時計数データに基づいて、幾何学配置由来の感度係数を求める第３感度係数算出工程とを行い、前記第１感度係数算出工程，前記第２感度係数算出工程および前記第３感度係数算出工程で求められた感度係数を用いて、放射線検出素子の感度補正を行うことを特徴とするものである。
　［作用・効果］この発明に係る放射線断層撮影装置によれば、この発明に係る感度補正方法でも述べたような第１感度係数算出工程，第２感度係数算出工程および第３感度係数算出工程を行う演算手段を備えている。この発明に係る感度補正方法でも述べたように、第１感度係数算出工程，第２感度係数算出工程および第３感度係数算出工程で求められた感度係数を用いて、放射線検出素子の感度補正を行うことにより、感度補正後のデータに基づき再構成画像を得る。その結果、検出器を構成する複数の放射線検出素子の配置が可変で円形配置でない場合であってもアーティファクトのない再構成画像を得ることができる。
　この発明に係る放射線検出素子の感度補正方法および放射線断層撮影装置によれば、（複数の放射線検出素子で収集された同時計数データに基づいて、感度係数を感度要因別に分けて求める）第１感度係数算出工程，（第１感度係数算出工程の後で、複数の放射線検出素子の配置の変化前に、収集された同時計数データに基づいて、非幾何学系因子の感度係数を求める）第２感度係数算出工程および（第２感度係数算出工程の後で、複数の放射線検出素子の配置が変化した状態で収集された同時計数データに基づいて、幾何学配置由来の感度係数を求める）第３感度係数算出工程で求められた感度係数を用いて、放射線検出素子の感度補正を行うことにより、感度補正後のデータに基づき再構成画像を得る。その結果、検出器を構成する複数の放射線検出素子の配置が可変で円形配置でない場合であってもアーティファクトのない再構成画像を得ることができる。
実施例に係るマルチモダリティ対応フレキシブルＰＥＴ装置の検出器ヘッドの一実施態様を示す概略斜視図およびブロック図である。 γ線検出器の概略斜視図である。 （ａ）は検出器ヘッドを奥行き方向から見た概略図、（ｂ）は奥行き方向に直交する水平方向から見た概略図である。 検出器配置および校正用線源の種類の模式図であり、（ａ）は非近接時のリング線源の模式図、（ｂ）は非近接時の円筒線源の模式図、（ｃ）は近接時の平板線源の模式図である。 実施例に係る要素別感度補正法の処理手順の流れを示したフローチャートである。 検出器ヘッドのγ線検出器に沿って回転可能な棒状線源を検出器ヘッドのγ線検出器に沿って回転させた時の模式図である。 （ａ）～（ｃ）は変形例に係る検出器配置および校正用線源の種類の模式図である。
　以下、図面を参照してこの発明の実施例を説明する。
　図１は、実施例に係るマルチモダリティ対応フレキシブルＰＥＴ装置の検出器ヘッドの一実施態様を示す概略斜視図およびブロック図であり、図２は、γ線検出器の概略斜視図であり、図３（ａ）は、検出器ヘッドを奥行き方向から見た概略図であり、図３（ｂ）は、奥行き方向に直交する水平方向から見た概略図である。
　図１に示すように、マルチモダリティ対応フレキシブルＰＥＴ装置（以下、単に「ＰＥＴ装置」と略記する）１は、検出器ヘッド２Ａ，２Ｂを備えている。図２に示すように、縦横に16×16に配列された検出器結晶が４層に積層された検出器ブロック３１（γ線検出器３）を円弧状に18個、奥行き方向ｚ（図１を参照）に３個配列された合計18×3個の検出器ブロック３１で１つの検出器ヘッドを構成する。検出器ブロック３１以外のγ線検出器３の具体的な構成については図２で後述する。マルチモダリティ対応フレキシブルＰＥＴ装置（ＰＥＴ装置）１は、この発明における放射線断層撮影装置に相当し、検出器ヘッド２Ａ，２Ｂは、この発明における検出器ヘッドに相当し、γ線検出器３は、この発明における放射線検出素子に相当する。
　図１に示すように、検出器ヘッド２Ａ，２Ｂは、互いに対向して構成されている。検出器ヘッド２Ａ，２Ｂは、通常は48角形の円周上に位置する。２つの検出器ヘッド２Ａ，２Ｂの断面内の検出器ブロック３１（γ線検出器３）数は合計36個（円弧状に18個×2つ）に対して48角形なので、左右に６個分ずつの隙間領域４が存在する。「発明が解決しようとする課題」の欄でも述べたように、図３では検出器ヘッド２Ａ，２Ｂは互いに対向した状態で検出器ヘッド２Ａ，２Ｂを互いに近接させることができる。このように、各検出器ヘッド２Ａ，２Ｂ間の距離が変更可能に各々の検出器ヘッド２Ａ，２Ｂは構成されている。
　その他にも、ＰＥＴ装置１は、同時計数回路５と演算回路６とを備えている。演算回路６の具体的な機能については図４で後述する。演算回路６は、この発明における演算手段に相当する。
　放射性薬剤が投与された被検体（図示省略）に放射性薬剤を投与すると、ポジトロン放出型のＲＩのポジトロンが消滅することにより、２本のγ線が発生する。同時計数回路５は、検出器ブロック３１（図２を参照）の位置とγ線の入射タイミングとをチェックし、被検体の両側にある２つの検出器ブロック３１でγ線が同時に入射したときのみ、送り込まれた画像情報を適正なデータ（同時計数データ）と判定する。一方の検出器ブロック３１のみにγ線が入射したときには、同時計数回路５は棄却する。つまり、同時計数回路５は、上述した電気信号に基づいて、２つのγ線検出器３においてγ線が同時観測（すなわち同時計数）されたことを検出する。同時計数回路５で収集された画像情報（同時計数データ）を、演算回路６に送り込む。
　γ線検出器３は、図２に示すように検出器ブロック３１と、その検出器ブロック３１に対して光学的に結合されたライトガイド３２と、そのライトガイド３２に対して光学的に結合された光電子増倍管(PMT: Photo Multiplier Tube)３３とを備えている。検出器ブロック３１を構成する各検出器結晶は、γ線の入射に伴って発光することでγ線から光に変換する。この変換によって検出器結晶はγ線を検出する。検出器結晶において発光した光が検出器ブロック３１で十分に拡散されて、ライトガイド３２を介して光電子増倍管（以下、単に「ＰＭＴ」と略記する）３３に入力される。ＰＭＴ３３は、検出器ブロック３１で変換された光を増倍させて電気信号に変換する。その電気信号は、上述したように画像情報（同時計数データ）として同時計数回路５（図１を参照）に送り込まれる。
　また、γ線検出器３は、図２に示すように、深さ方向に複数の層からなるＤＯＩ検出器である。図２では、４層のＤＯＩ検出器を図示しているが、層の数については、複数であれば特に限定されない。また、縦横に配列された検出器結晶の数についても、複数であれば特に限定されない。
　ここで、ＤＯＩ検出器は、各々の検出器結晶を放射線の深さ方向に積層して構成されたものであり、相互作用を起こした深さ(DOI: Depth of Interaction)方向と横方向（入射面に平行な方向）との座標情報を重心演算により求める。ＤＯＩ検出器を用いることにより深さ方向の空間分解能をより一層向上させることができる。よって、ＤＯＩ検出器の層の数は、深さ方向に積層された検出器結晶の層の数である。
　次に、演算回路６の具体的な機能について、図４～図６を参照して説明する。図４（ａ）～図４（ｃ）は、検出器配置および校正用線源の種類の模式図であり、図４（ａ）は、非近接時のリング線源の模式図であり、図４（ｂ）は、非近接時の円筒線源の模式図であり、図４（ｃ）は、近接時の平板線源の模式図であり、図５は、実施例に係る要素別感度補正法の処理手順の流れを示したフローチャートであり、図６は、検出器ヘッドのγ線検出器に沿って回転可能な棒状線源を検出器ヘッドのγ線検出器に沿って回転させた時の模式図である。
　図４（ａ）中の符号ＳＲは校正用リング線源（以下、単に「リング線源」と略記する）であり、図４（ｂ）中の符号ＳＣは校正用円筒線源（以下、単に「円筒線源」と略記する）であり、図４（ｃ）中の符号ＳＰは校正用平板線源（以下、単に「平板線源」と略記する）である。なお、本実施例では、図６に示すように、検出器ヘッド２Ａ，２Ｂのγ線検出器３（図２を参照）（すなわち放射線検出素子）に沿って回転可能な校正用棒状線源（以下、単に「棒状線源」と略記する）ＳＬを用いて、棒状線源ＳＬを検出器ヘッド２Ａ，２Ｂのγ線検出器３に沿って回転させることで図４（ａ）に示すリング線源ＳＲに見立てる。以下、図６に示す棒状線源ＳＬを図４（ａ）に示すリング線源ＳＲとして取り扱って説明する。リング線源ＳＲの厚みは約5mm，円筒線源ＳＣの直径は約15cm，平板線源ＳＰの厚みは約5mm，棒状線源ＳＬの直径は約15cmである。
　また、本実施例では、図４（ａ）に示す非近接時に棒状線源ＳＬ（図６を参照）を回転させてリングに見立てたリング線源ＳＲ，図４（ｂ）に示す非近接時に装置（検出器ヘッド２Ａ，２Ｂ）の中心に設置された均一な円筒線源ＳＣ，図４（ｃ）に示す近接時に装置（検出器ヘッド２Ａ，２Ｂ）の中心に設置された平板線源ＳＰの３種類を用意して、それぞれでデータを収集する。校正用線源はアクリル樹脂で製作され、内部に均一濃度の放射能溶液が密封されている。
　また、図５では、長軸断のブロック形状因子ｂｕｖｋが「Transaxial Block Profile Factor」に対応し、結晶干渉因子ｄｕｖｒｋが「Crystal Interference Factor」に対応し、体軸断面（体軸横断面）のブロック形状因子ｂｕｖｋが「Axial Block Profile Factor」に対応し、リングペア感度ｆｕｖが「Ring Pair Factor」に対応し、検出器固有感度εｕｉ，εｖｊが「Crystal Efficiency Factor」に対応し、平板線源の時の結晶干渉因子ｄPlateが「Crystal Interference(Plate) Factor」に対応し、平板線源の時のリングペア感度ｆPlateが「Ring Pair(Plate) Factor」に対応する。
　（ステップＳ１）
　図５のフローチャートに示すように、最初に図４（ａ）に示すリング線源ＳＲで収集された同時計数データを使って、経時変化しない幾何学系因子の結晶干渉因子ｄｕｖｒｋ（「Crystal Interference Factor」で表記）を計算する。この幾何学系因子の結晶干渉因子ｄｕｖｒｋ（Crystal Interference Factor）を計算する前にパス長規格化を適用（「Apply」で表記）するのが好ましい。
　パス長規格化とは、所定のパス上にある同時計数データを、そのパスが線源を通過する長さを用いて規格化することである。パス長規格化を行うことで、所定のパス間の線源と交差する長さの違いが擬似的に補正される。具体的なパス長規格化については、特許文献１：特開２００８－１７０３２９号公報を参照されたい。
　さらに、検出器距離を補正するための立体角補正を行う。具体的には、あるＬＯＲの結晶干渉因子ｄｕｖｒｋ（Crystal Interference Factor）に対して、そのＬＯＲの検出器距離の二乗で割る。その結果で得られた結晶干渉因子ｄｕｖｒｋ（Crystal Interference Factor）の係数をテーブルとして出力する。ここで、ＬＯＲ(Line Of Response)は、同時計数する２つの検出器を結ぶ仮想上の直線である。このステップＳ１は、この発明における第１感度係数算出工程に相当する。
　（ステップＳ２）
　次に、結晶干渉因子ｄｕｖｒｋ（Crystal Interference Factor）による幾何学系因子の補正を行った図４（ｂ）に示す円筒線源ＳＣの同時計数データを使って、経時変化する非幾何学系因子の検出器固有感度εｕｉ，εｖｊ（「Crystal Efficiency Factor」で表記）を順に計算する。このときに、体軸断面（体軸横断面）のブロック形状因子ｂｕｖｋ（「Axial Block Profile Factor」で表記）や長軸断のブロック形状因子ｂｕｖｋ（「Transaxial Block Profile Factor」で表記）も計算し、ステップＳ３で図４（ｃ）に示す平板線源ＳＰの同時計数データに対して補正を行う。
　また、リングペア感度ｆｕｖ（「Ring Pair Factor」で表記）も計算する。計算で得られた各因子をテーブルとして出力する。このステップＳ２は、この発明における第２感度係数算出工程に相当する。
　（ステップＳ３）
　最後に、図４（ｃ）に示す平板線源ＳＰの同時計数データを使って、近接したことによって検出器の幾何学配置の変化による影響を受ける平板線源の時の結晶干渉因子ｄPlate（「Crystal Interference(Plate) Factor」で表記）および平板線源の時のリングペア感度ｆPlate（「Ring Pair(Plate) Factor」で表記）のみを計算する。
　つまり、同一の検出器ヘッドの組み合わせは、γ線検出器３（放射線検出素子）の配置によらず位置関係が不変であるので、ステップＳ３で結晶干渉因子ｄｕｖｒｋ（Crystal Interference Factor）およびリングペア感度ｆｕｖ（Ring Pair Factor）を計算しない。一方、互いに対向する検出器ヘッドの組み合わせは、γ線検出器３（放射線検出素子）の位置関係により幾何学的な影響を受けるので、平板線源の時の結晶干渉因子ｄPlate（Crystal Interference(Plate) Factor）および平板線源の時のリングペア感度ｆPlate（Ring Pair(Plate) Factor）を計算する。
　ステップＳ１と同様に、これらの因子を計算する前にパス長規格化を適用するのが好ましい。平板線源の時の結晶干渉因子ｄPlate（Crystal Interference(Plate) Factor）に対して、近接時の座標で立体角補正を行う。計算で得られた平板線源の時の結晶干渉因子ｄPlate（Crystal Interference(Plate) Factor）および平板線源の時のリングペア感度ｆPlate（Ring Pair(Plate) Factor）の係数をテーブルとして出力する。このステップＳ３は、この発明における第３感度係数算出工程に相当する。
　各因子の具体的な計算については、非特許文献１を参照されたい。なお、フレキシブルＰＥＴ装置は検出器ヘッド間に隙間が存在するので、ファンサム(fan-sum)法を利用する。ファンサム法を利用した長軸断のブロック形状因子ｂｕｖｋ（Transaxial Block Profile Factor）やリングペア感度ｆｕｖ（Ring Pair Factor）の具体的な計算については、特許文献２：国際公開第２０１０／０１３３５６号を参照されたい。
　近接時の平板線源と非近接時のリング線源との差を補正する補正係数αを求める式は、下記（２）式のように表される。
Figure JPOXMLDOC01-appb-M000001
　なお、上記（２）式中のＮｎｅａｒは近接時の平板線源を用いて収集された同時計数データ（カウント値）、Ｎｎｏｒｍａｌは非近接時の平板線源を用いて収集された同時計数データ（カウント値）、Ｌｎｅａｒは近接時の検出器距離、Ｌｎｏｒｍａｌは非近接時の検出器距離である。このように、上記（２）式の右辺である平板線源とリング線源との差は、検出器距離で表す幾何学的な要因（この場合にはＬｎｅａｒ、Ｌｎｏｒｍａｌ）とそれ以外の補正係数αとの積で表される。上記（２）式の左辺は、近接時・非近接時それぞれで同じ平板線源を用いて同時間に同時計数データを収集して実験的に求める。本実施例では、楕円に変形した短辺を検出器間距離Ｌｎｅａｒとして採用している。
　本実施例に係る感度補正方法によれば、先ず、第１感度係数算出工程（ステップＳ１）で感度係数を求めて、第１感度係数算出工程（ステップＳ１）の後で、複数の放射線検出素子（本実施例ではγ線検出器）の配置の変化前（本実施例では非近接時）に第２感度係数算出工程（ステップＳ２）で感度係数を求めて、第２感度係数算出工程（ステップＳ２）の後で、複数の放射線検出素子の配置が変化した（本実施例では検出器ヘッド２Ａ，２Ｂが近接した）状態で、第３感度係数算出工程（ステップＳ３）で感度係数を求める。具体的には、第１感度係数算出工程（ステップＳ１）で感度係数を感度要因別に分けて求める。複数の放射線検出素子（γ線検出器）の配置が変化しなくとも、非幾何学系因子は経時変化する。そこで、複数の放射線検出素子（γ線検出器）の配置の変化前（非近接時）に、収集された同時計数データに基づいて、第２感度係数算出工程（ステップＳ２）で非幾何学系因子の感度係数を求める。そして、第３感度係数算出工程（ステップＳ３）では、複数の放射線検出素子の配置が変化した（検出器ヘッド２Ａ，２Ｂが近接した）状態で収集された同時計数データに基づいて、幾何学配置由来の感度係数を求める。このように、第１感度係数算出工程（ステップＳ１），第２感度係数算出工程（ステップＳ２）および第３感度係数算出工程（ステップＳ３）で求められた感度係数を用いて、放射線検出素子（γ線検出器）の感度補正を行うことにより、感度補正後のデータに基づき再構成画像を得る。その結果、検出器を構成する複数の放射線検出素子（γ線検出器）の配置が可変で円形配置でない場合であってもアーティファクトのない再構成画像を得ることができる。
　本実施例では、複数の放射線検出素子の配置が変化した（検出器ヘッド２Ａ，２Ｂが近接した）場合の感度係数は、放射線検出素子（γ線検出器）の配置によらず位置関係が不変の組み合わせ（本実施例では同一の検出器ヘッド）に対しては、配置変化前（本実施例では非近接時）に求められた幾何学系因子の感度係数と、非幾何学系因子の感度係数との積で表される。そして、放射線検出素子（γ線検出器）の位置関係により幾何学的な影響を受ける組み合わせ（本実施例では互いに対向する検出器ヘッド）に対しては、配置変化後（本実施例では近接時）に求められた幾何学系因子の感度係数と、非幾何学系因子の感度係数との積で表される。
　なお、検出器ヘッド２Ａ，２Ｂが近接した場合の感度係数を求める際には、検出器ペアが同一の検出器ヘッドか互いに対向する検出器ヘッドかをＰＥＴ装置１（例えば中央演算処理装置（ＣＰＵ））が判定する。対向する検出器ヘッドであれば、図４（ｃ）に示す平板線源ＳＰから求められた近接時の感度係数に切り替える。
　上記（２）式のように複数の放射線検出素子（γ線検出器）の配置の変化前後（近接時・非近接時）で同じ校正用線源（本実施例では同じ平板線源）を用いて同時間にそれぞれ収集された同時計数データの比（上記（２）式の左辺）から補正係数αを求めるのが好ましい。ひいては、例えば配置の変化前での別の校正用線源（前者：第１校正用線源、本実施例ではリング線源）を用いた後に、複数の放射線検出素子（γ線検出器）の配置の変化前後（近接時・非近接時）で同じ校正用線源（後者：第２校正用線源、本実施例では平板線源）を用いる場合には、求められた補正係数αを前者（リング線源）での上記（１）式の右辺に乗算することで、配置の変化前の前者（非近接時のリング線源）と配置の変化後の後者（近接時の平板線源）との差をも補正することができる。
　また、複数の放射線検出素子（γ線検出器）の配置の変化後（近接時）に収集された同時計数データは、少なくとも放射線検出素子（γ線検出器）の位置関係により幾何学的な影響を受ける組み合わせである。すなわち、配置の変化後（近接時）に収集された同時計数データに基づいて、放射線検出素子（γ線検出器）の位置関係により幾何学的な影響を受ける組み合わせ（本実施例では互いに対向する検出器ヘッド）に対する感度係数データを求めることができる。
　また、第１感度係数算出工程（ステップＳ１）では、下記のような校正用線源を用いて感度係数を求めるのが好ましい。すなわち、本実施例のように、複数の放射線検出素子（γ線検出器）に沿って回転可能な棒状線源ＳＬ（図６を参照）を用いる場合には、棒状線源ＳＬを複数の放射線検出素子（γ線検出器）に沿って回転させることで図４（ａ）に示すリング線源ＳＲに見立てる。第１感度係数算出工程（ステップＳ１）では、回転した棒状線源ＳＬ（すなわちリング線源ＳＲに見立てられた棒状線源ＳＬ）からの放射線（本実施例ではγ線）を検出することで収集された同時計数データに基づいて、感度係数を感度要因別に分けて求める。
　また、第２感度係数算出工程（ステップＳ２）では、下記のような校正用線源を用いて感度係数を求めるのが好ましい。すなわち、第２感度係数算出工程（ステップＳ２）では、複数の放射線検出素子（γ線検出器）の配置の変化前（非近接時）に、実際の被検体に模した円筒線源ＳＣからの放射線（γ線）を検出することで収集された同時計数データに基づいて、経時変化する非幾何学系因子の感度係数を求めることができる。
　また、複数の放射線検出素子（γ線検出器）からなる検出器ヘッドを複数（本実施例では２つ）備え、各当該検出器ヘッド間の距離が変更可能に各々の検出器ヘッドが構成されている場合で、各検出器ヘッドを互いに近接させることで複数の放射線検出素子（γ線検出器）の配置が変化するときには、下記のような校正用線源を用いて感度係数を求めるのが好ましい。すなわち、各検出器ヘッドを互いに近接させようとすると、複数の放射線検出素子（γ線検出器）に沿って回転可能な棒状線源や複数の放射線検出素子（γ線検出器）に沿って設けられたリング線源が各検出器ヘッドと干渉する恐れがある。そこで、本実施例のように、各検出器ヘッドの中心に設置された図４（ｃ）に示す平板線源ＳＰを用いる。第３感度係数算出工程（ステップＳ３）では、各検出器ヘッドを互いに近接させた状態で、平板線源ＳＰからの放射線（γ線）を検出することで収集された同時計数データに基づいて、幾何学配置由来の感度係数を求める。このような平板線源ＳＰを用いることで、各検出器ヘッドを近接させても、校正用線源と各検出器ヘッドとの干渉を防止することができる。
　また、本実施例に係る放射線断層撮影装置（マルチモダリティ対応フレキシブルＰＥＴ装置１）によれば、本実施例に係る感度補正方法でも述べたような第１感度係数算出工程（ステップＳ１），第２感度係数算出工程（ステップＳ２）および第３感度係数算出工程（ステップＳ３）を行う演算手段（図１では演算回路６）を備えている。本実施例に係る感度補正方法でも述べたように、第１感度係数算出工程（ステップＳ１），第２感度係数算出工程（ステップＳ２）および第３感度係数算出工程（ステップＳ３）で求められた感度係数を用いて、放射線検出素子（γ線検出器）の感度補正を行うことにより、感度補正後のデータに基づき再構成画像を得る。その結果、検出器を構成する複数の放射線検出素子（γ線検出器）の配置が可変で円形配置でない場合であってもアーティファクトのない再構成画像を得ることができる。
　［再構成結果］
　装置（検出器ヘッド２Ａ，２Ｂ）の直径778mmから直径250mmに近接させて、装置（検出器ヘッド２Ａ，２Ｂ）中心に直径300mmの均一放射能濃度の円筒線源を設置して、同時計数データを収集した。リング線源の幾何学因子と円筒線源の非幾何学系因子とを使った従来の要素別感度補正法を組み込んだ再構成結果では、非近接の幾何学因子を使うので、再構成画像にアーティファクトが生じる。円筒線源の非幾何学系因子，同一の検出器ヘッド内のイベントはリング線源由来，互いに対向する検出器ヘッド内のイベントは平板線源由来に切り替える幾何学系因子，さらに補正係数αを入れた再構成結果でアーティファクトがなくなり画質が改善したのが確認された。
　この発明は、上記実施形態に限られることはなく、下記のように変形実施することができる。
　（１）上述した実施例では、放射線はγ線であったが、α線やβ線などの放射線であってもよい。
　（２）上述した実施例では、ＤＯＩ検出器であったが、深さ方向を弁別しない放射線検出器に適用してもよい。
　（３）上述した実施例では、２つの検出器ヘッド２Ａ，２Ｂ（図２～図４および図６を参照）について説明したが、検出器ヘッドの数については２つに限定されない。複数の放射線検出素子（実施例ではγ線検出器）からなる検出器ヘッドを複数備え、各当該検出器ヘッド間の距離が変更可能に各々の検出器ヘッドが構成される構造であれば、３つ以上であってもよい。
　（４）上述した実施例では、複数の放射線検出素子（実施例ではγ線検出器）に沿って回転可能な棒状線源ＳＬ（図６を参照）を用いて、棒状線源ＳＬを複数の放射線検出素子（γ線検出器）に沿って回転させることで図４（ａ）に示すリング線源ＳＲに見立てたが、図４（ａ）に示すリング線源ＳＲそのものや、図７（ａ）に示す平板線源ＳＰを複数の放射線検出素子（γ線検出器）に沿って回転させてもよい。すなわち、図４（ａ）に示すように複数の放射線検出素子（γ線検出器）に沿って設けられたリング線源ＳＲを用いてもよい。また、図７（ａ）に示すように複数の放射線検出素子（γ線検出器）に沿って回転可能な平板線源ＳＰを用いる場合には、平板線源ＳＰを複数の放射線検出素子（γ線検出器）に沿って回転させてもよい。
　（５）上述した実施例では、近接時の校正用平板線源として、各検出器ヘッドの中心に設置された図４（ｃ）に示す平板線源ＳＰを用いたが、図７（ｂ）に示す各検出器ヘッドの中心軸に沿って往復可能な棒状線源ＳＬを用いてもよい。すなわち、図７（ｂ）に示すように各検出器ヘッドの中心軸に沿って往復可能な棒状線源ＳＬを用いる場合には、棒状線源ＳＬを各検出器ヘッドの中心軸に沿って往復させることで平板線源ＳＰに見立てる。
　（６）上述した実施例では、図４（ａ）に示すように真円状のリング線源ＳＲを用いたが、変形した検出器配置に即した形状の校正用リング線源（ここでは図７（ｃ）に示す楕円形のリング線源ＳＲ）でもよい。また、複数の放射線検出素子（実施例ではγ線検出器）に沿って楕円軌道に回転可能な棒状線源を用いて、棒状線源を複数の放射線検出素子（γ線検出器）に沿って楕円軌道に回転させることで図７（ｃ）に示す楕円形のリング線源ＳＲに見立ててもよい。
　（７）上述した実施例では、複数の放射線検出素子（実施例ではγ線検出器）の配置の変化前は非近接時で、配置の変化後は近接時であったが逆でもよい。すなわち、配置の変化前は近接時で、配置の変化後は非近接時でもよい。よって、第１感度係数算出工程（ステップＳ１）では、図７（ｃ）に示す楕円形のリング線源ＳＲを用いて近接時に収集された同時計数データに基づいて、感度係数を感度要因別に分けて求める。そして、第３感度係数算出工程（ステップＳ３）では、非近接時に収集された同時計数データに基づいて、幾何学配置由来の感度係数を求める。
　（８）第１感度係数算出工程（ステップＳ１），第２感度係数算出工程（ステップＳ２）および第３感度係数算出工程（ステップＳ３）での校正用線源の形状については特に限定されない。検出器配置に即した形状の校正用線源を用いればよい。
　１　…　マルチモダリティ対応フレキシブルＰＥＴ装置（ＰＥＴ装置）
　２Ａ，２Ｂ　…　検出器ヘッド
　３　…　γ線検出器
　ＳＲ　…　校正用リング線源（リング線源）
　ＳＣ　…　校正用円筒線源（円筒線源）
　ＳＰ　…　校正用平板線源（平板線源）
　ＳＬ　…　校正用棒状線源（棒状線源）
　α　…　補正係数

Claims (10)

  1. 　被検体内の放射性薬剤から発生した放射線を検出する放射線検出素子の感度補正方法であって、
    　複数の前記放射線検出素子で収集された同時計数データに基づいて、感度係数を感度要因別に分けて求める第１感度係数算出工程と、
    　前記第１感度係数算出工程の後で、前記複数の放射線検出素子の配置の変化前に、収集された同時計数データに基づいて、非幾何学系因子の感度係数を求める第２感度係数算出工程と、
    　前記第２感度係数算出工程の後で、前記複数の放射線検出素子の配置が変化した状態で収集された同時計数データに基づいて、幾何学配置由来の感度係数を求める第３感度係数算出工程と
    　を備え、
    　前記第１感度係数算出工程，前記第２感度係数算出工程および前記第３感度係数算出工程で求められた感度係数を用いて、放射線検出素子の感度補正を行うことを特徴とする感度補正方法。
  2. 　請求項１に記載の感度補正方法において、
    　前記複数の放射線検出素子の配置が変化した場合の感度係数は、放射線検出素子の配置によらず位置関係が不変の組み合わせに対しては、配置変化前に求められた幾何学系因子の感度係数と、非幾何学系因子の感度係数との積で表され、放射線検出素子の位置関係により幾何学的な影響を受ける組み合わせに対しては、配置変化後に求められた幾何学系因子の感度係数と、非幾何学系因子の感度係数との積で表されることを特徴とする感度補正方法。
  3. 　請求項１に記載の感度補正方法において、
    　前記複数の放射線検出素子の配置の変化前後で同じ校正用線源を用いて同時間にそれぞれ収集された同時計数データの比から補正係数を求めることを特徴とする感度補正方法。
  4. 　請求項３に記載の感度補正方法において、
    　前記複数の放射線検出素子の配置の変化前に所定の校正用線源を第１校正用線源として用いた後に、前記第１校正用線源とは別の校正用線源を第２校正用線源としたときに、前記複数の放射線検出素子の配置の変化前後で同じ第２校正用線源を用いて同時間にそれぞれ収集された同時計数データの比から補正係数を求め、
    　前記第１校正用線源での感度補正に関する各因子の積で表された式に、前記補正係数を乗算することを特徴とする感度補正方法。
  5. 　請求項１に記載の感度補正方法において、
    　前記複数の放射線検出素子の配置の変化後に収集された同時計数データは、少なくとも放射線検出素子の位置関係により幾何学的な影響を受ける組み合わせであることを特徴とする感度補正方法。
  6. 　請求項１に記載の感度補正方法において、
    　前記第１感度係数算出工程では、複数の前記放射線検出素子に沿って回転可能な校正用棒状線源，複数の前記放射線検出素子に沿って設けられた校正用リング線源または複数の前記放射線検出素子に沿って回転可能な校正用平板線源からの放射線を検出することで収集された同時計数データに基づいて、感度係数を感度要因別に分けて求めて、
    　前記第２感度係数算出工程では、前記複数の放射線検出素子の配置の変化前に、校正用円筒線源からの放射線を検出することで収集された同時計数データに基づいて、非幾何学系因子の感度係数を求めることを特徴とする感度補正方法。
  7. 　請求項１に記載の感度補正方法において、
    　同時計数する２つの検出器を結ぶ仮想上の直線であるＬＯＲの結晶干渉因子に対して、当該ＬＯＲの検出器距離の二乗で割ることによって、検出器距離を補正するための立体角補正を行うことを特徴とする感度補正方法。
  8. 　請求項１から請求項７のいずれかに記載の感度補正方法において、
    　複数の前記放射線検出素子からなる検出器ヘッドを複数備え、各当該検出器ヘッド間の距離が変更可能に各々の検出器ヘッドは構成されており、
    　各検出器ヘッドの中心に設置された校正用平板線源または各検出器ヘッドの中心軸に沿って往復可能な校正用棒状線源であって、
    　前記第３感度係数算出工程では、各検出器ヘッドを互いに近接させた状態で、当該校正用平板線源または往復した当該校正用棒状線源からの放射線を検出することで収集された同時計数データに基づいて、幾何学配置由来の感度係数を求めることを特徴とする感度補正方法。
  9. 　被検体内の放射性薬剤から発生した放射線を検出して放射線断層撮影を行う放射線断層撮影装置であって、
    　複数の放射線検出素子からなる検出器ヘッドを複数備え、各当該検出器ヘッド間の距離が変更可能に各々の検出器ヘッドは構成されており、
    　演算処理を行う演算手段を備え、
    　当該演算手段は、
    　複数の前記放射線検出素子で収集された同時計数データに基づいて、感度係数を感度要因別に分けて求める第１感度係数算出工程と、
    　前記第１感度係数算出工程の後で、前記複数の放射線検出素子の配置の変化前に、収集された同時計数データに基づいて、非幾何学系因子の感度係数を求める第２感度係数算出工程と、
    　前記第２感度係数算出工程の後で、前記複数の放射線検出素子の配置が変化した状態で収集された同時計数データに基づいて、幾何学配置由来の感度係数を求める第３感度係数算出工程と
    　を行い、
    　前記第１感度係数算出工程，前記第２感度係数算出工程および前記第３感度係数算出工程で求められた感度係数を用いて、放射線検出素子の感度補正を行うことを特徴とする放射線断層撮影装置。
  10. 　請求項９に記載の放射線断層撮影装置において、
    　前記検出器ヘッドは、円弧状に各々の前記放射線検出素子を配置して構成されていることを特徴とする放射線断層撮影装置。
PCT/JP2016/061172 2015-06-18 2016-04-05 放射線検出素子の感度補正方法および放射線断層撮影装置 WO2016203822A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US15/736,168 US10126445B2 (en) 2015-06-18 2016-04-05 Radiation detecting element sensitivity correction method and radiation tomography device
JP2017524688A JP6607256B2 (ja) 2015-06-18 2016-04-05 放射線検出素子の感度補正方法および放射線断層撮影装置
CN201680035775.XA CN107710018B (zh) 2015-06-18 2016-04-05 放射线检测元件的灵敏度校正方法以及放射线断层摄影装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2015122651 2015-06-18
JP2015-122651 2015-06-18

Publications (1)

Publication Number Publication Date
WO2016203822A1 true WO2016203822A1 (ja) 2016-12-22

Family

ID=57545087

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2016/061172 WO2016203822A1 (ja) 2015-06-18 2016-04-05 放射線検出素子の感度補正方法および放射線断層撮影装置

Country Status (4)

Country Link
US (1) US10126445B2 (ja)
JP (1) JP6607256B2 (ja)
CN (1) CN107710018B (ja)
WO (1) WO2016203822A1 (ja)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020530333A (ja) * 2017-08-09 2020-10-22 リフレクション メディカル， インコーポレイテッド 放出誘導放射線療法における異常検出のためのシステムおよび方法
EP3845929A1 (en) * 2020-01-06 2021-07-07 Canon Medical Systems Corporation A pet apparatus and a method
US11627920B2 (en) 2008-03-14 2023-04-18 Reflexion Medical, Inc. Method and apparatus for emission guided radiation therapy
US11675097B2 (en) 2017-07-11 2023-06-13 Reflexion Medical, Inc. Methods for PET detector afterglow management
US11794036B2 (en) 2016-11-15 2023-10-24 Reflexion Medical, Inc. Radiation therapy patient platform
US11904184B2 (en) 2017-03-30 2024-02-20 Reflexion Medical, Inc. Radiation therapy systems and methods with tumor tracking
US11975220B2 (en) 2016-11-15 2024-05-07 Reflexion Medical, Inc. System for emission-guided high-energy photon delivery
US12029921B2 (en) 2017-11-14 2024-07-09 Reflexion Medical, Inc. Systems and methods for patient monitoring for radiotherapy
JP7517966B2 (ja) 2020-01-06 2024-07-17 キヤノンメディカルシステムズ株式会社 Ｐｅｔ装置および算出方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6824133B2 (ja) 2017-09-28 2021-02-03 富士フイルム株式会社 画像処理装置、画像処理方法、及び画像処理プログラム

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006189274A (ja) * 2005-01-04 2006-07-20 Shimadzu Corp 核医学イメージング装置
JP2008170329A (ja) * 2007-01-12 2008-07-24 Shimadzu Corp 放射線同時計数処理方法および断層撮影装置
JP2015099149A (ja) * 2013-11-19 2015-05-28 株式会社東芝 核医学診断装置及び結晶効率補正係数の計算方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3793266B2 (ja) * 1995-10-20 2006-07-05 浜松ホトニクス株式会社 ポジトロンｃｔ装置およびその画像再構成方法
JP4536212B2 (ja) * 2000-05-24 2010-09-01 浜松ホトニクス株式会社 Ｐｅｔ装置
US20070116183A1 (en) * 2003-07-30 2007-05-24 Hironori Ueki Tomograph
JP4983981B2 (ja) * 2008-05-19 2012-07-25 株式会社島津製作所 ２次元位置マップ校正方法および放射線検出装置
JP4968384B2 (ja) 2008-07-31 2012-07-04 株式会社島津製作所 放射線断層撮影装置
US8461518B2 (en) * 2009-01-23 2013-06-11 Shimadzu Corporation Method of collecting calibration data in radiation tomography apparatus
JP5750685B2 (ja) * 2011-01-11 2015-07-22 国立研究開発法人放射線医学総合研究所 Ｐｅｔ装置及びｐｅｔ−ｍｒｉ装置
WO2012127761A1 (ja) * 2011-03-22 2012-09-27 株式会社日立メディコ Ｘ線ｃｔ装置、および方法
US8483992B2 (en) * 2011-04-28 2013-07-09 Kabushiki Kaisha Toshiba Analytical computation of geometric normalization factor for TOF-PET
US8786474B1 (en) * 2013-03-15 2014-07-22 Kabushiki Kaisha Toshiba Apparatus for programmable metastable ring oscillator period for multiple-hit delay-chain based time-to-digital circuits

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006189274A (ja) * 2005-01-04 2006-07-20 Shimadzu Corp 核医学イメージング装置
JP2008170329A (ja) * 2007-01-12 2008-07-24 Shimadzu Corp 放射線同時計数処理方法および断層撮影装置
JP2015099149A (ja) * 2013-11-19 2015-05-28 株式会社東芝 核医学診断装置及び結晶効率補正係数の計算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
KEISHI KITAMURA ET AL.: "Jisedai PET Sochi Kaihatsu no Shinchoku Jokyo (8) jPET-D4 ni Okeru Kakushu Data Hoseiho no Kaihatsu", NIRS-R(NATL INST RADIOL SCI), NO.NIRS-R-54, March 2006 (2006-03-01), pages 47 - 51 *
KEISHI KITAMURA ET AL.: "jPET-D4 ni Okeru Yosobetsu Kando Hoseiho no Kaihatsu", NIRS- M(NATL INST RADIOL SCI, 1 March 2005 (2005-03-01), pages 43 - 48 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11627920B2 (en) 2008-03-14 2023-04-18 Reflexion Medical, Inc. Method and apparatus for emission guided radiation therapy
US11794036B2 (en) 2016-11-15 2023-10-24 Reflexion Medical, Inc. Radiation therapy patient platform
US11975220B2 (en) 2016-11-15 2024-05-07 Reflexion Medical, Inc. System for emission-guided high-energy photon delivery
US11904184B2 (en) 2017-03-30 2024-02-20 Reflexion Medical, Inc. Radiation therapy systems and methods with tumor tracking
US11675097B2 (en) 2017-07-11 2023-06-13 Reflexion Medical, Inc. Methods for PET detector afterglow management
US12032107B2 (en) 2017-07-11 2024-07-09 Reflexion Medical, Inc. Methods for PET detector afterglow management
US11511133B2 (en) 2017-08-09 2022-11-29 Reflexion Medical, Inc. Systems and methods for fault detection in emission-guided radiotherapy
JP7315961B2 (ja) 2017-08-09 2023-07-27 リフレクション メディカル， インコーポレイテッド 放出誘導放射線療法における異常検出のためのシステムおよび方法
JP2020530333A (ja) * 2017-08-09 2020-10-22 リフレクション メディカル， インコーポレイテッド 放出誘導放射線療法における異常検出のためのシステムおよび方法
US12023523B2 (en) 2017-08-09 2024-07-02 Reflexion Medical, Inc. Systems and methods for fault detection in emission-guided radiotherapy
US12029921B2 (en) 2017-11-14 2024-07-09 Reflexion Medical, Inc. Systems and methods for patient monitoring for radiotherapy
EP3845929A1 (en) * 2020-01-06 2021-07-07 Canon Medical Systems Corporation A pet apparatus and a method
JP7517966B2 (ja) 2020-01-06 2024-07-17 キヤノンメディカルシステムズ株式会社 Ｐｅｔ装置および算出方法

Also Published As

Publication number Publication date
US10126445B2 (en) 2018-11-13
JP6607256B2 (ja) 2019-11-20
CN107710018B (zh) 2019-08-13
CN107710018A (zh) 2018-02-16
JPWO2016203822A1 (ja) 2018-03-29
US20180196151A1 (en) 2018-07-12

Similar Documents

Publication Publication Date Title
JP6607256B2 (ja) 放射線検出素子の感度補正方法および放射線断層撮影装置
JP4536212B2 (ja) Ｐｅｔ装置
US9119589B2 (en) Method and system for spectral computed tomography (CT) with sparse photon counting detectors
Tashima et al. Proposed helmet PET geometries with add-on detectors for high sensitivity brain imaging
CN108289647B (zh) 用于稀疏检测器的装置、方法和***
US9709684B2 (en) Systems and methods for scintillators having micro-crack surfaces
JP2005326406A (ja) 陽電子放出断層システムの規格化のための方法及びシステム
JP2013511022A5 (ja)
Jin et al. Experimental evaluation of a 3-D CZT imaging spectrometer for potential use in Compton-Enhanced PET imaging
KR101677715B1 (ko) 방사선 영상화 방법 및 시스템
JP2012173237A (ja) 陽電子放射断層撮像装置
CN107076862B (zh) 射线检测器、射线摄像装置、计算机断层摄影装置以及射线检测方法
US20170108593A1 (en) Device for determining three dimensional locations and energy of gamma incidence events and method for the same
US9390823B2 (en) Radiation image acquiring device
JP2012177555A (ja) 放射線撮像装置
CN111596336A (zh) 基于缝-孔平板准直器的多伽马光子符合成像***及方法
JP2017058191A (ja) 放射線撮像装置
US20200085397A1 (en) Attenuation coefficient image estimation method, attenuation coefficient image estimation program, and positron ct apparatus equipped with the same
Tanaka et al. A 2D/3D hybrid PET scanner with rotating partial slice-septa and its quantitative procedures
Kotasidis et al. Image-based spatially variant and count rate dependent point spread function on the HRRT
US20220334268A1 (en) Compton imaging apparatus and single photon emission and positron emission tomography system comprising same
EP3845929B1 (en) A pet apparatus and a method
McNamara et al. Positron emission tomography coincidence detection with photon polarization correlation
WO2016162962A1 (ja) 放射線撮像装置
Hagiwara et al. A performance study of the next generation PET using a new histogramming method for the DOI detector

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: 16811297

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2017524688

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 15736168

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 16811297

Country of ref document: EP

Kind code of ref document: A1