WO2015056299A1 - 断層画像処理方法およびそれを用いた放射型断層撮影装置 - Google Patents

断層画像処理方法およびそれを用いた放射型断層撮影装置 Download PDF

Info

Publication number
WO2015056299A1
WO2015056299A1 PCT/JP2013/077962 JP2013077962W WO2015056299A1 WO 2015056299 A1 WO2015056299 A1 WO 2015056299A1 JP 2013077962 W JP2013077962 W JP 2013077962W WO 2015056299 A1 WO2015056299 A1 WO 2015056299A1
Authority
WO
WIPO (PCT)
Prior art keywords
data
projection data
tomographic image
projection
detector
Prior art date
Application number
PCT/JP2013/077962
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 PCT/JP2013/077962 priority Critical patent/WO2015056299A1/ja
Priority to JP2015542427A priority patent/JP6206501B2/ja
Publication of WO2015056299A1 publication Critical patent/WO2015056299A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
    • G01T1/164—Scintigraphy
    • G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
    • G01T1/1647—Processing of scintigraphic data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/037—Emission tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image

Definitions

  • the present invention relates to a tomographic image processing method for performing processing related to a radial tomographic image and a radial tomographic apparatus using the same.
  • the SPECT apparatus reconstructs a tomographic image of a subject by detecting a single radiation ( ⁇ ray).
  • the PET device detects only multiple rays ( ⁇ rays) generated by the annihilation of positrons (positrons) and detects the radiation ( ⁇ rays) with multiple detectors simultaneously (that is, only when they are counted simultaneously). Reconstruct a tomographic image of the subject.
  • the PET apparatus is also called a “positron tomography apparatus”.
  • a PET apparatus positron tomography apparatus
  • the projection data shown in the lower diagram of FIG. 9 is a sinogram, the horizontal axis represents the detector arrangement S, and the vertical axis represents the projection angle ⁇ of the ⁇ -ray. Since the gamma ray pairs are counted simultaneously, the range of the projection angle ⁇ is not 0 ° to 360 °, but half of that is 0 ° to 180 °.
  • the detector unit moves the entire subject from 0 ° to 180 °. Cover with °. Therefore, the projection data obtained thereby becomes complete projection data (complete projection data).
  • a detector unit configured to be partially opened, for example, a PET having a so-called “partial ring type” detector unit in which a part of the detector as shown in FIG. 9B is missing.
  • the minimum data (complete projection data) necessary for accurately reconstructing the biodistribution image of the radiopharmaceutical cannot be measured.
  • Incomplete projection data (incomplete projection data) obtained by such a measurement has a defect area (denoted by reference sign D in FIG. 9B) as shown in FIG. 9B. Therefore, when reconstruction calculation processing is performed on incomplete projection data, a false image (artifact) derived from data loss occurs in the reconstructed image (tomographic image).
  • a C-shaped or U-shaped detector unit used in a mammographic PET apparatus for acquiring a tomographic image of a breast, or a plurality of detectors face each other. Then, as exemplified by the detector units arranged in parallel, the same phenomenon as in the partial ring type occurs if the detector unit is configured to be partially opened. That is, an artifact is caused by a missing part of the detector.
  • the detection time difference is a detection time difference between detected ⁇ -ray pairs.
  • the time difference from the pair annihilation occurrence point to the detector is the distance from the pair annihilation occurrence point to the light source generation position by the scintillator element of the detector. It is a technique that estimates the point of occurrence of annihilation by converting to a difference.
  • Non-Patent Document 1 is an academic paper that is the basis of Patent Document 1.
  • Patent Document 1 relates to the system configuration of a partial ring type TOF-PET apparatus as shown in FIG. 9B.
  • TOF detected time difference
  • Non-Patent Document 1 which is an academic paper that is the basis of Patent Document 1, data called “list mode” that stores ⁇ -ray detection event information (detector number, detection time, ⁇ -ray energy, etc.) in time series Reconfiguration calculation processing is performed on (also referred to as “list mode data” or “list data”).
  • This reconstruction method is called “list mode reconstruction”.
  • a reconstruction method for projection data obtained by converting the list mode data into a histogram for each detection position is referred to as “projection data reconstruction”.
  • Non-Patent Document 2 relates to an algorithm for estimating a projection value of a missing portion of incomplete projection data. Such an algorithm is called a “projection completion algorithm”.
  • Patent Document 1 Non-Patent Document 1
  • Non-Patent Document 2 described above have the following problems.
  • Non-Patent Document 2 a fundamental solution to the problem of image quality degradation in the partial ring type PET apparatus is to complete projection data using a projection perfection algorithm as described in Non-Patent Document 2.
  • the projection data reconstruction as in Non-Patent Document 2 has the following problems.
  • (2) conversion of list mode data into projection data becomes a cause of deterioration of the spatial resolution of the reconstructed image (tomographic image).
  • the data amount of the list mode data does not depend on the number of scintillator crystals but is proportional to the count of ⁇ rays (count number).
  • list mode reconstruction is becoming mainstream in the PET image reconstruction method.
  • the projection data is not created, there is of course no turn of the projection perfection algorithm. Therefore, data loss is not compensated for by the conventional list mode reconstruction.
  • the present invention has been made in view of such circumstances, and an object of the present invention is to provide a tomographic image processing method capable of improving image quality by reducing artifacts and a radiation tomography apparatus using the same.
  • the tomographic image processing method of the present invention is a tomographic image processing method for performing processing related to a radiation type tomographic image, and includes list mode data generated from event data obtained by detecting radiation, and radiation.
  • An acquisition processing step for performing processing for acquiring the tomographic image by performing a successive approximation method that sequentially approximates and updates an image using both of the projection data obtained by detection is provided. It is.
  • the successive approximation method is performed in which images are sequentially approximated and updated using both list mode data and projection data. That is, the successive approximation method is performed using a hybrid method (hereinafter referred to as “hybrid reconstruction method”) that combines list mode reconstruction and projection data reconstruction. Therefore, data loss can be compensated by projection data, and artifacts can be reduced and image quality can be improved compared to conventional list mode reconstruction.
  • hybrid reconstruction method combines list mode reconstruction and projection data reconstruction. Therefore, data loss can be compensated by projection data, and artifacts can be reduced and image quality can be improved compared to conventional list mode reconstruction.
  • the above-described missing area data estimated in the estimation processing step is used as the above-described projection data used in the above-described acquisition processing step.
  • the object of compensation is list mode data in the hybrid reconstruction method. Reconstruction is the main subject. Therefore, since only a part of the missing area data takes the form of projection data, the data quality itself does not deteriorate. Therefore, the image quality of the hybrid reconstruction method using the projection data (defect region data) estimated and completed in the estimation processing step is superior to the projection data reconstruction using the complete projection data.
  • the radiation tomography apparatus of the present invention is a successive approximation method using both the above-described list mode data and the above-mentioned projection data in the above-described radiation tomography apparatus using the tomographic image processing methods of these inventions.
  • a tomographic image processing means for performing processing for acquiring a tomographic image is provided.
  • the tomographic image processing means performs the successive approximation method using the hybrid reconstruction method described above, so that the data loss can be compensated for by the projection data, and the conventional list mode reconstruction can be performed. Compared to the configuration, the image quality can be improved by reducing artifacts.
  • the successive approximation method is performed using the hybrid reconstruction method (combining the list mode reconstruction and the projection data reconstruction).
  • Data loss can be compensated by projection data, and artifacts can be reduced and image quality can be improved compared to conventional list mode reconstruction.
  • FIG. 1 is a side view and block diagram of a PET (Positron Emission Tomography) apparatus according to an embodiment. It is a schematic perspective view of a gamma ray detector. It is a schematic front view of a partial ring type detector unit. It is a flowchart which shows the flow of the output of a reconstruction image (tomographic image) from a series of imaging
  • (A), (b) is a schematic front view of the detector unit brought close to the subject.
  • (A) is a schematic front view of a full ring type detector unit and projection data at that time
  • (b) is a schematic front view of a partial ring type detector unit and projection data at that time.
  • FIG. 1 is a side view and a block diagram of a PET (Positron Emission Tomography) apparatus according to an embodiment
  • FIG. 2 is a schematic perspective view of a ⁇ -ray detector
  • FIG. 3 is a partial ring type detector unit.
  • a PET apparatus positron tomography apparatus
  • FIG. 3 FIG. 9B
  • a top plate 1 on which the subject M is placed is disposed beside the PET apparatus according to the present embodiment.
  • the top plate 1 is configured to move up and down and translate along the body axis Z of the subject M.
  • the subject M placed on the top 1 is scanned from the head to the abdomen and foot sequentially through the opening 2a of the gantry 2, which will be described later. Get the image.
  • the PET apparatus which concerns on a present Example does not include the top plate 1 as a structure, you may provide the top plate 1 as a structure.
  • the PET apparatus includes a gantry 2 having an opening 2a and a ⁇ -ray detector 3.
  • the ⁇ -ray detector 3 is arranged in a partial ring shape so as to surround the body axis Z of the subject M, and is embedded in the gantry 2.
  • the respective ⁇ -ray detectors 3 are arranged so as to form a partial ring type detector unit 30 that is partially opened.
  • the detector unit 30 constituting the ⁇ -ray detector 3 corresponds to the detector unit in the present invention.
  • the PET apparatus includes a top board driving unit 4, a controller 5, an input unit 6, an output unit 7, a memory unit 8, a coincidence circuit 9, and a GPU (Graphics Processing Unit) 10. Yes.
  • the top plate driving unit 6 is a mechanism for driving the top plate 1 so as to perform the above-described movement, and is configured by a motor or the like not shown.
  • the GPU 10 corresponds to the tomographic image processing means in this invention.
  • the controller 5 comprehensively controls each part constituting the PET apparatus according to the present embodiment.
  • the controller 5 includes a central processing unit (CPU).
  • the input unit 6 sends data and commands input by the operator to the controller 5.
  • the input unit 6 includes a pointing device represented by a mouse, a keyboard, a joystick, a trackball, a touch panel, and the like.
  • the output unit 7 includes a display unit represented by a monitor, a printer, and the like.
  • the memory unit 8 includes a storage medium represented by ROM (Read-only Memory), RAM (Random-Access Memory), and the like.
  • ROM Read-only Memory
  • RAM Random-Access Memory
  • the count value (count) simultaneously counted by the coincidence circuit 9 the data relating to the coincidence counting such as the detector pair consisting of the two ⁇ -ray detectors 3 and the LOR, and the calculation processing performed by the GPU 10.
  • Various data and the like are written and stored in the RAM, and read from the RAM as necessary.
  • the ROM stores in advance programs for performing various control processes (for example, drive control of the top board 1) and image processes (for example, tomographic image processes), and the controller 5 and the GPU 10 execute the programs. Then, control processing and image processing corresponding to the program are performed.
  • a program related to the estimation process to be estimated is stored in the ROM in advance, and the GPU 10 executes a program related to the acquisition process, the conversion process, and the estimation process, thereby executing processes in steps S1 to S4 described later.
  • the ⁇ -rays generated from the subject M to which the radiopharmaceutical is administered are converted into light by the scintillator block 31 (see FIG. 2) of the ⁇ -ray detector 3, and the converted light is photoelectron of the ⁇ -ray detector 3.
  • a multiplier tube (PMT: Photo Multiplier Tube) 33 (see FIG. 2) multiplies and converts it into an electrical signal. The electric signal is sent to the coincidence circuit 9 as an event.
  • the coincidence circuit 9 checks the position of the scintillator block 31 (see FIG. 2) and the incident timing of the ⁇ rays, and only when the ⁇ rays are simultaneously incident on the two scintillator blocks 31 on both sides of the subject M. The sent event is determined to be appropriate data.
  • the coincidence counting circuit 9 rejects. That is, the coincidence counting circuit 9 detects that ⁇ rays are simultaneously observed in the two ⁇ ray detectors 3 based on the above-described electrical signal.
  • the event sent to the coincidence circuit 9 is sent to the GPU 10.
  • the GPU 10 obtains a tomographic image of the subject M by performing image reconstruction through acquisition processing, conversion processing, and estimation processing.
  • the tomographic image is sent to the output unit 7 via the controller 5. In this way, tomography is performed based on the tomographic image obtained by the GPU 10. Specific functions of the GPU 10 will be described later.
  • the ⁇ -ray detector 3 includes a scintillator block 31, a light guide 32 optically coupled to the scintillator block 31, and photoelectrons optically coupled to the light guide 32.
  • a multiplier (hereinafter simply abbreviated as “PMT”) 33 is provided.
  • Each scintillator element constituting the scintillator block 31 converts ⁇ rays into light by emitting light with the incidence of ⁇ rays. By this conversion, the scintillator element detects ⁇ rays.
  • Light emitted from the scintillator element is sufficiently diffused by the scintillator block 31 and input to the PMT 33 via the light guide 32.
  • the PMT 33 multiplies the light converted by the scintillator block 31 and converts it into an electric signal. The electric signal is sent to the coincidence circuit 9 (see FIG. 1) as an event as described above.
  • the ⁇ -ray detector 3 shown in FIG. 2 is a DOI detector configured by laminating each scintillator block 31 in the depth direction of the ⁇ -ray (in FIG. 2, four layers).
  • the DOI detector is configured by laminating the respective scintillator blocks 31 in the depth direction of the ⁇ -ray, and the depth direction and the lateral direction (direction parallel to the incident surface) in which the interaction has occurred. Is obtained by the center of gravity calculation. This makes it possible to discriminate the light source position (DOI: Depth of Interaction) in the depth direction where the interaction has occurred.
  • FIG. 4 is a flowchart showing a flow of output of a reconstructed image (tomographic image) from a series of imaging in the partial ring type PET apparatus
  • FIG. 5 is a gamma ray detection for explaining a detector response function (detection probability). It is the schematic diagram which showed the coincidence count in a vessel.
  • the scintillator block 31 is illustrated as the ⁇ -ray detector 3, and the light guide 32 and the PMT 33 are not illustrated.
  • Step S1 (List Mode Data) Acquisition Processing
  • the subject M (see FIG. 1) is imaged by a PET apparatus (partial ring type PET apparatus) provided with the partial ring type detector unit 30 shown in FIG. .
  • the ⁇ -ray detector 3 of the detector unit 30 Only when the ⁇ -ray detector 3 of the detector unit 30 simultaneously counts ⁇ -rays, the event data having the position of the scintillator block 31 (see FIG. 2) and the incident timing of the ⁇ -rays is output to the coincidence circuit 9 (FIG. 1), and the ⁇ -ray detection event information including the detector number, detection time, ⁇ -ray energy, etc. is stored in time series.
  • list mode data generated from event data obtained by detecting ⁇ -rays is acquired.
  • Step S2 Conversion Process
  • the list mode data obtained in step S1 is converted into projection data.
  • the projection data converted from the list mode data obtained by the partial ring type detector unit 30 is incomplete projection data.
  • This step S2 corresponds to the conversion processing step in the present invention.
  • Step S3 Estimation Processing The missing area data is estimated based on the projection data converted in step S2. Then, the projection data is completed.
  • the method of the projection perfection algorithm that estimates the missing area data and completes the projection data is not particularly limited.
  • the projection data is completed using a projection perfection algorithm as described in Non-Patent Document 2 (see Fig. 4 of Non-Patent Document 2).
  • Non-Patent Document 2 when incomplete projection data is Fourier-transformed, a signal having a large amplitude due to the fact that the projection data is incomplete in a spatial frequency region where only a signal having a small amplitude is observed if the projection data is complete. Appears. This is set to “0” and inverse Fourier transform is performed to estimate missing area data. By repeating this calculation, the accuracy of the estimated value of the missing area data is improved.
  • This step S3 corresponds to the estimation processing step in this invention.
  • Step S4 Tomographic Image Acquisition Process
  • a sequential approximation method is performed in which the image is sequentially approximated and updated.
  • the ML-EM method Maximum Likelihood Expectation Maximization
  • the successive approximation method is not limited to the ML-EM method, but may be a DRAMA method (Dynamic Row-Action Maximum Likelihood Algorithm) or a static (that is, static) RAMLA method (Row-Action Maximum Likelihood Algorithm).
  • the OSEM method Ordered Subset ML-EM
  • the update formula of the ML-EM method is expressed by the following formula (3) described later. That is, the solution of the optimization problem (evaluation function minimization problem) is obtained as the pixel value of the reconstructed image (tomographic image).
  • the evaluation function for the conventional list mode reconstruction is L (x)
  • the evaluation function L (x) is defined by the following equation (1).
  • the evaluation function F (x) is defined by the following equation (2). Further, for the evaluation function L (x) that constitutes the evaluation function F (x), the evaluation function L (x) defined by the above equation (1) is used.
  • the evaluation function F (x) of the hybrid reconstruction method defined by the above equation (2) is estimated by the projection perfection algorithm to the evaluation function L (x) of the list mode reconstruction defined by the above equation (1).
  • LOR Line Of Response
  • the detector response function a ij is a probability of being detected by the i-th crystal pair.
  • the detector response function b ij is a probability of being detected by the i-th virtual LOR.
  • the initial image x j (1) is set appropriately.
  • the initial image x j (1) may be an image having a uniform pixel value, for example, and x j (1) > 0.
  • X j (k) are sequentially obtained by repeatedly substituting into the above equation (3).
  • X j in order.
  • the number of times k representing repetition is not particularly limited, and may be set as appropriate.
  • x j finally obtained is arranged for each pixel j corresponding to it, thereby performing image reconstruction and obtaining a reconstructed image (tomographic image) of the subject M.
  • This step S4 corresponds to the acquisition processing step in this invention.
  • the GPU 10 executes the processing of steps S1 to S4 in FIG. 4 as described above.
  • the tomographic image acquired by the GPU 10 is sent to the output unit 7 (see FIG. 1) via the controller 5 (see FIG. 1).
  • the tomographic image acquired by the GPU 10 may be written and stored in the memory unit 8 (see FIG. 1) via the controller 5 and read from the memory unit 8 as necessary.
  • a sequential approximation method is performed in which images are approximated and updated sequentially using both list mode data and projection data. That is, the successive approximation method is performed using a hybrid method (hybrid reconstruction method) that combines list mode reconstruction and projection data reconstruction. Therefore, data loss can be compensated by projection data, and artifacts can be reduced and image quality can be improved compared to conventional list mode reconstruction.
  • the conversion process step (step S2 in FIG. 4) for converting the list mode data into the projection data described above, and the missing region data based on the projection data converted in the conversion processing step (step S2).
  • Projection processing step (step S3 in FIG. 4) for estimating the above-mentioned missing area data estimated in the estimation processing step (step S3) and projection data used in the acquisition processing step (step S4 in FIG. 4) Used as Although the spatial resolution (image quality) of the reconstructed image (tomographic image) is deteriorated by converting the list mode data into projection data in the conversion processing step (step S2), the object of compensation is list mode data in the hybrid reconstruction method. Therefore, the list mode reconstruction is the main subject.
  • the hybrid reconstruction method using the projection data (missing area data) estimated and completed in the estimation processing step (step S3) is more suitable than the projection data reconstruction using the completed projection data.
  • the image quality is excellent.
  • the tomographic image processing means (GPU 10 in the present embodiment) performs the successive approximation method using the above-described hybrid reconstruction method. Can be compensated, and compared to conventional list mode reconstruction, artifacts can be reduced and image quality can be improved.
  • the present invention is not limited to the above embodiment, and can be modified as follows.
  • the PET apparatus positron tomography apparatus
  • the present invention detects tomographic images of a subject by detecting a single radiation. It can also be applied to a SPECT (Single-Photon-Emission-CT) apparatus that reconfigures the image. The present invention can also be applied to a PET-CT apparatus that combines a PET apparatus and a CT apparatus.
  • SPECT Single-Photon-Emission-CT
  • the radiation tomography apparatus (PET apparatus in the embodiment) that acquires a tomographic image of the whole body of the subject has been described.
  • the imaging target is not limited to the whole body. You may apply to the PET apparatus for heads which acquires the tomographic image of the head of a subject, the PET apparatus for mammons which acquires the tomographic image of the breast of a subject.
  • the DOI detector is composed of a plurality of scintillator elements arranged in three dimensions.
  • the present invention is also applicable to a radiation detector composed of a plurality of scintillator elements arranged in two dimensions or three dimensions. can do.
  • the fixed type radiation tomography apparatus (PET apparatus in the embodiment) has been described. However, it is installed adjacent to another modality apparatus (for example, a CT apparatus) (not shown). Therefore, the present invention may be applied to a movable radiation tomography apparatus (here, a PET apparatus) as shown in FIG.
  • a movable radiation tomography apparatus here, a PET apparatus
  • FIG. 6 the housing in which the detector unit 30 is embedded is a support arm 31 and is formed in a C shape in an arc shape.
  • the arm holding part 32 holds the support arm 31 and is supported by the carriage 33.
  • a rear wheel 34 and a front wheel 35 are attached to the bottom of the carriage 33, and the carriage 33 is configured to be transportable by moving these on the floor.
  • a motor (not shown) is connected to the front wheel 35 via a drive shaft (not shown).
  • the front wheel 35 is driven by driving the motor, and the surgeon pushes in an arbitrary direction from behind the carriage 33.
  • the carriage 33 By pulling and rolling the rear wheel 34, the carriage 33 can be moved in any direction on the floor surface.
  • the partial ring type detector unit shown in FIG. 3 has been described as an example of the detector unit that is partially opened.
  • the partial ring type shown in FIG. It is not limited to the detector unit.
  • the detector units may be arranged opposite to each other in parallel.
  • the present invention may be applied to a detector unit 30 that is sandwiched between both sides as shown in FIG.
  • the detection plate 41 in which the detector unit 30 is embedded is a notch, and the breast is inspected by being sandwiched by this notch.
  • a plurality of ⁇ -ray detectors 3 (not shown in FIG. 7) are arranged in parallel in the detection plate 41 in accordance with the notches.
  • FIG. 8A shows a structure of a detector unit 30 in which a plurality of detectors ( ⁇ -ray detector 3 in the embodiment) are arranged in parallel to face each other.
  • FIG. 8B shows the structure of a partial ring type detector unit 30 similar to FIG. 8A and 8B, a part of the detector unit 30 is configured to be open, and the detector units 30 are separated from each other independently. Therefore, the detector unit 30 can be brought close to the subject M. Moreover, the effect that the detection efficiency (sensitivity) of a radiation becomes high by making it adjoin is also show
  • the tomographic image processing unit is configured by the GPU, but is not limited to the GPU.
  • the tomographic image processing means is composed of a central processing unit (CPU) and programmable devices (for example, FPGA (Field Programmable Gate Gate Array)) that can change hardware circuits (for example, logic circuits) used in accordance with program data. May be.
  • CPU central processing unit
  • programmable devices for example, FPGA (Field Programmable Gate Gate Array)
  • hardware circuits for example, logic circuits

Landscapes

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

Abstract

　この発明のＰＥＴ装置のＧＰＵは、ステップＳ１～Ｓ４の処理を行う。具体的には、リストモードデータを投影データに変換する変換処理工程（ステップＳ２）と、その変換処理工程（ステップＳ２）で変換された投影データに基づいて欠損領域データを推定する推定処理工程（ステップＳ３）とを備え、その推定処理工程（ステップＳ３）で推定された上述の欠損領域データを取得処理工程（ステップＳ４）で用いられる投影データとして用いる。変換処理工程（ステップＳ２）においてリストモードデータを投影データに変換することにより再構成画像（断層画像）の空間分解能（画質）が劣化するが、一部の欠損領域データのみが投影データの形式をとっているだけなのでデータの質自体は劣化しない。よって、投影データによりデータ欠損を補償することができ、従来のリストモード再構成と比較すると、アーティファクトを低減させて画質を改善することができる。

Description

断層画像処理方法およびそれを用いた放射型断層撮影装置
　この発明は、放射型の断層画像に関する処理を行う断層画像処理方法およびそれを用いた放射型断層撮影装置に関する。
　放射型断層撮影装置として、ＳＰＥＣＴ装置（Single Photon Emission CT）やＰＥＴ装置(Positron Emission Tomography)がある。ＳＰＥＣＴ装置は、単一の放射線（γ線）を検出して被検体の断層画像を再構成する。ＰＥＴ装置は、陽電子（ポジトロン）の消滅によって発生する複数本の放射線（γ線）を検出して複数個の検出器で放射線（γ線）を同時に検出したときのみ（つまり同時計数したときのみ）被検体の断層画像を再構成する。特に、ＰＥＴ装置は「ポジトロン断層撮影装置」とも呼ばれている。以下では、ＰＥＴ装置（ポジトロン断層撮影装置）を例に採って説明する。
　図９の下図に示す投影データはサイノグラム(sinogram)であり、横軸は検出器の配置Ｓを表し、縦軸はγ線の投影角度φである。γ線対を同時計数するので、投影角度φの範囲は０°～３６０°でなく、その半分の０°～１８０°となる。
　図９（ａ）に示すように円環状に各々の検出器が並べられた、いわゆる「フルリング型」の検出器ユニットを備えたＰＥＴ装置では、検出器ユニットが被検体全体を０°～１８０°でカバーする。よって、それによって得られた投影データは完全な投影データ（完全投影データ）となる。
　一方、一部が開口して構成された検出器ユニットにおいて、例えば図９（ｂ）に示すような検出器の一部に抜けがある、いわゆる「部分リング型」の検出器ユニットを備えたＰＥＴ装置では、その抜けを通過して検出されないγ線が存在する。この場合には、放射性薬剤の体内分布画像を正確に再構成するために必要な最低限のデータ（完全投影データ）を測定することができない。そのように測定されて得られた不完全な投影データ（不完全投影データ）には、図９（ｂ）に示すように欠損領域（図９（ｂ）では符号Ｄで表記）が生じる。したがって、不完全投影データに対して再構成の計算処理を施すと、再構成画像（断層画像）にデータ欠損に由来した偽像（アーティファクト）が生じる。
　なお、図９（ｂ）に示す部分リング型以外にも、乳房の断層画像を取得するマンモ用ＰＥＴ装置において用いられるＣ型やコの字型の検出器ユニットや、複数の検出器が互いに対向して平行に配置された検出器ユニット等に例示されるように、一部が開口して構成された検出器ユニットであれば、部分リング型と同様の現象が生じる。すなわち、検出器の一部に抜けがあることによりアーティファクトが生じる。
　そこで、アーティファクトを低減させるために、検出時間差(TOF: Time Of Flight)の情報を利用した技術（例えば、特許文献１、非特許文献１参照）や、投影データを完全化する（すなわち欠損領域データを推定する）技術（例えば、非特許文献２参照）がある。検出時間差（ＴＯＦ）とは、検出されたγ線対の検出時間差である。ポジトロンの対消滅発生地点を、消滅γ線が光速であることを用いて、対消滅発生地点から検出器に到達する時間差を、対消滅発生地点から検出器のシンチレータ素子による光源発生位置までの距離差に換算することで、対消滅発生地点を推定する技術である。なお、非特許文献１は、特許文献１の元となる学術論文である。
　特許文献１は、図９（ｂ）に示したような部分リング型ＴＯＦ－ＰＥＴ装置のシステム構成に関する内容である。検出されたγ線対の検出時間差（ＴＯＦ）の情報を再構成の計算処理において利用することで、データ欠損による情報量の低下を補償し、再構成画像（断層画像）のアーティファクトを低減させることを目的としている。
　特許文献１の元となる学術論文である非特許文献１では、γ線の検出イベント情報（検出器番号、検出時間、γ線のエネルギ等）を時系列で保存した「リストモード」と呼ばれるデータ（「リストモードデータ」または「リストデータ」とも呼ぶ）に対して再構成の計算処理を施している。この再構成方式を「リストモード再構成」という。一方、リストモードデータを検出位置毎にヒストグラム化して得られる投影データに対する再構成方式を「投影データ再構成」という。
　非特許文献２は、不完全投影データの欠損部分の投影値を推定するアルゴリズムに関する。このようなアルゴリズムは「投影完全化（projection completion）アルゴリズム」と呼ばれる。
米国特許出願公開第２０１０／０１０８８９６号明細書
S Surti and JS Karp: Design considerations for a limited-angle, dedicated breast, TOF PET scanner. Phys Med Biol, vol. 53, 2911-2921, 2008. J. S. Karp, G. Muehllehner and R. M. Lewitt: Constrained Fourier space method for compensation of missing data in emission computed tomography. IEEE Trans Med Imag, vol. 7, 22-25, 1998.
　しかしながら、上述した特許文献１や非特許文献１や非特許文献２の場合には下記のような問題点がある。
　特許文献１や非特許文献１で述べられているように、部分リング型ＰＥＴ装置において、ＴＯＦの情報を利用した再構成の計算処理により画像のアーティファクトは低減するが、低減の程度はＴＯＦの情報の精度（検出器の時間分解能）およびデータ欠損量（検出器がない部分の大きさ）に依存する。したがって、ＴＯＦの情報を利用することが、満足のいく画質が得られるということには必ずしもならない。つまり、ＴＯＦの情報を利用することは、部分リング型ＰＥＴ装置の画質低下の問題に対する根本的な解決策とならない。
　欠損領域データを推定して投影データを完全化することができれば、理論上、画像にアーティファクトは発生しない。したがって、部分リング型ＰＥＴ装置における画質低下の問題の根本的な解決策は、非特許文献２にあるような投影完全化アルゴリズムを利用して投影データを完全化することである。
　しかし、非特許文献２のような投影データ再構成では、以下のような問題点がある。（１）検出器（シンチレータ結晶）の微細化に伴い結晶数が増加し、結晶数の対（ペア）の分だけ投影データのデータ量が指数関数的に増加してしまう。例えば、結晶数が２倍になると結晶数の対（ペア）の組み合わせや投影データのデータ量は４倍（＝２２倍）となり、結晶数が３倍になると結晶数の対（ペア）の組み合わせや投影データのデータ量は９倍（＝３２倍）となる。また、（２）リストモードデータを投影データに変換することが再構成画像（断層画像）の空間分解能の劣化要因となってしまう。一方、リストモードデータのデータ量はシンチレータ結晶数に依存せず、γ線の計数（カウント数）に比例する。
　以上の理由により、ＰＥＴ画像の再構成方式ではリストモード再構成が主流になりつつある。その場合には、投影データを作成しないので、当然ながら投影完全化アルゴリズムの出番はない。したがって、従来のリストモード再構成ではデータ欠損は補償されない。
　この発明は、このような事情に鑑みてなされたものであって、アーティファクトを低減させて画質を改善することができる断層画像処理方法およびそれを用いた放射型断層撮影装置を提供することを目的とする。
　この発明は、このような目的を達成するために、次のような構成をとる。
　すなわち、この発明の断層画像処理方法は、放射型の断層画像に関する処理を行う断層画像処理方法であって、放射線を検出することにより得られた事象データから生成されたリストモードデータ、および放射線を検出することにより得られた投影データの両方を用いて、画像を逐次に近似して更新する逐次近似法を行って前記断層画像を取得する処理を行う取得処理工程を備えることを特徴とするものである。
　この発明の断層画像処理方法によれば、リストモードデータおよび投影データの両方を用いて、画像を逐次に近似して更新する逐次近似法を行っている。すなわち、リストモード再構成および投影データ再構成を組み合わせたハイブリッド方式（以下、「ハイブリッド再構成法」と呼ぶ）を用いて逐次近似法を行っている。よって、投影データによりデータ欠損を補償することができ、従来のリストモード再構成と比較すると、アーティファクトを低減させて画質を改善することができる。
　具体的には、上述のリストモードデータを上述の投影データに変換する変換処理工程と、その変換処理工程で変換された投影データに基づいて欠損領域データを推定する推定処理工程とを備え、その推定処理工程で推定された上述の欠損領域データを上述の取得処理工程で用いられる上述の投影データとして用いる。変換処理工程においてリストモードデータを投影データに変換することにより再構成画像（断層画像）の空間分解能（画質）が劣化するが、ハイブリッド再構成法では補償の対象はリストモードデータであってリストモード再構成が主体である。したがって、一部の欠損領域データのみが投影データの形式をとっているだけなのでデータの質自体は劣化しない。よって、完全化された投影データを用いた投影データ再構成よりも、推定処理工程において推定されて完全化された投影データ（欠損領域データ）を用いたハイブリッド再構成法の方が画質は優れる。
　また、この発明の放射型断層撮影装置は、上述したこれらの発明の断層画像処理方法を用いた放射型断層撮影装置において、上述のリストモードデータおよび上述の投影データの両方を用いて逐次近似法を行って断層画像を取得する処理を行う断層画像処理手段を備えることを特徴とするものである。
　この発明の放射型断層撮影装置によれば、断層画像処理手段は上述のハイブリッド再構成法を用いて逐次近似法を行うので、投影データによりデータ欠損を補償することができ、従来のリストモード再構成と比較すると、アーティファクトを低減させて画質を改善することができる。
　この発明に係る断層画像処理方法およびそれを用いた放射型断層撮影装置によれば、（リストモード再構成および投影データ再構成を組み合わせた）ハイブリッド再構成法を用いて逐次近似法を行うので、投影データによりデータ欠損を補償することができ、従来のリストモード再構成と比較すると、アーティファクトを低減させて画質を改善することができる。
実施例に係るＰＥＴ(Positron Emission Tomography)装置の側面図およびブロック図である。 γ線検出器の概略斜視図である。 部分リング型の検出器ユニットの概略正面図である。 部分リング型ＰＥＴ装置における一連の撮影から再構成画像（断層画像）の出力の流れを示すフローチャートである。 検出器応答関数（検出確率）の説明に供するγ線検出器での同時計数を示した模式図である。 変形例に係るＰＥＴ装置の側面図である。 さらなる変形例に係るＰＥＴ装置の側面図およびブロック図である。 （ａ）、（ｂ）は、被検体に対して近接させた検出器ユニットの概略正面図である。 （ａ）はフルリング型の検出器ユニットの概略正面図およびそのときの投影データ、（ｂ）は部分リング型の検出器ユニットの概略正面図およびそのときの投影データである。
　以下、図面を参照してこの発明の実施例を説明する。図１は、実施例に係るＰＥＴ(Positron Emission Tomography)装置の側面図およびブロック図であり、図２は、γ線検出器の概略斜視図であり、図３は、部分リング型の検出器ユニットの概略正面図である。本実施例では、放射型断層撮影装置として、ＰＥＴ装置（ポジトロン断層撮影装置）を例に採って説明するとともに、一部が開口して構成された検出器ユニットとして、図３（図９（ｂ）も参照）に示す部分リング型の検出器ユニットを例に採って説明する。
　本実施例に係るＰＥＴ装置の傍らには、図１に示すように、被検体Ｍを載置する天板１が配置されている。この天板１は、上下に昇降移動、被検体Ｍの体軸Ｚに沿って平行移動するように構成されている。このように構成することで、天板１に載置された被検体Ｍは、後述するガントリ２の開口部２ａを通って、頭部から順に腹部、足部へと走査されて、被検体Ｍの画像を得る。なお、走査される部位や各部位の走査順序については特に限定されない。なお、本実施例に係るＰＥＴ装置は、天板１を構成として含ませていないが、天板１を構成として含ませて備えてもよい。
　本実施例に係るＰＥＴ装置は、開口部２ａを有したガントリ２と、γ線検出器３とを備えている。γ線検出器３は、被検体Ｍの体軸Ｚ周りを取り囲むようにして部分リング状に配置されており、ガントリ２内に埋設されている。本実施例では、図３に示すように一部が開口して構成された部分リング型の検出器ユニット３０となるように、各々のγ線検出器３が並べられている。γ線検出器３を構成する検出器ユニット３０は、この発明における検出器ユニットに相当する。
　その他にも、本実施例に係るＰＥＴ装置は、天板駆動部４とコントローラ５と入力部６と出力部７とメモリ部８と同時計数回路９とＧＰＵ（Graphics Processing Unit）１０とを備えている。天板駆動部６は、天板１の上述した移動を行うように駆動する機構であって、図示を省略するモータなどで構成されている。ＧＰＵ１０は、この発明における断層画像処理手段に相当する。
　コントローラ５は、本実施例に係るＰＥＴ装置を構成する各部分を統括制御する。コントローラ５は、中央演算処理装置（ＣＰＵ）などで構成されている。
　入力部６は、オペレータが入力したデータや命令をコントローラ５に送り込む。入力部６は、マウスやキーボードやジョイスティックやトラックボールやタッチパネルなどに代表されるポインティングデバイスで構成されている。出力部７はモニタなどに代表される表示部やプリンタなどで構成されている。
　メモリ部８は、ＲＯＭ（Read-only Memory）やＲＡＭ（Random-Access Memory）などに代表される記憶媒体で構成されている。本実施例では、同時計数回路９で同時計数された計数値（カウント）や同時計数した２つのγ線検出器３からなる検出器対やＬＯＲといった同時計数に関するデータや、ＧＰＵ１０で演算処理された各種のデータなどについてはＲＡＭに書き込んで記憶し、必要に応じてＲＡＭから読み出す。ＲＯＭには、各種の制御処理（例えば天板１の駆動制御）や画像処理（例えば断層画像処理）を行うためのプログラム等を予め記憶しており、そのプログラムをコントローラ５およびＧＰＵ１０が実行することでそのプログラムに応じた制御処理や画像処理をそれぞれ行う。特に、本実施例では、後述する逐次近似法を行って断層画像を取得する処理を行う取得処理，リストモードデータを投影データに変換する変換処理，変換された投影データに基づいて欠損領域データを推定する推定処理に関するプログラムをＲＯＭに予め記憶しており、それらの取得処理や変換処理や推定処理に関するプログラムをＧＰＵ１０が実行することで後述するステップＳ１～Ｓ４の処理を実行する。
　放射性薬剤が投与された被検体Ｍから発生したγ線をγ線検出器３のシンチレータブロック３１（図２を参照）が光に変換して、変換されたその光をγ線検出器３の光電子増倍管(PMT: Photo Multiplier Tube)３３（図２を参照）は増倍させて電気信号に変換する。その電気信号をイベントとして同時計数回路９に送り込む。
　具体的には、被検体Ｍに放射性薬剤を投与すると、ポジトロン放出型のＲＩのポジトロンが消滅することにより、２本のγ線が発生する。同時計数回路９は、シンチレータブロック３１（図２を参照）の位置とγ線の入射タイミングとをチェックし、被検体Ｍの両側にある２つのシンチレータブロック３１でγ線が同時に入射したときのみ、送り込まれたイベントを適正なデータと判定する。一方のシンチレータブロック３１のみにγ線が入射したときには、同時計数回路９は棄却する。つまり、同時計数回路９は、上述した電気信号に基づいて、２つのγ線検出器３においてγ線が同時観測されたことを検出する。
　同時計数回路９に送り込まれたイベントを、ＧＰＵ１０に送り込む。ＧＰＵ１０は取得処理や変換処理や推定処理による画像再構成を行って、被検体Ｍの断層画像を求める。断層画像を、コントローラ５を介して出力部７に送り込む。このようにして、ＧＰＵ１０で得られた断層画像に基づいて断層撮影を行う。ＧＰＵ１０の具体的な機能については後述する。
　γ線検出器３は、図２に示すようにシンチレータブロック３１と、そのシンチレータブロック３１に対して光学的に結合されたライトガイド３２と、そのライトガイド３２に対して光学的に結合された光電子増倍管（以下、単に「ＰＭＴ」と略記する）３３とを備えている。シンチレータブロック３１を構成する各シンチレータ素子は、γ線の入射に伴って発光することでγ線から光に変換する。この変換によってシンチレータ素子はγ線を検出する。シンチレータ素子において発光した光がシンチレータブロック３１で十分に拡散されて、ライトガイド３２を介してＰＭＴ３３に入力される。ＰＭＴ３３は、シンチレータブロック３１で変換された光を増倍させて電気信号に変換する。その電気信号は、上述したようにイベントとして同時計数回路９（図１を参照）に送り込まれる。
　また、図２に示すγ線検出器３は、各々のシンチレータブロック３１をγ線の深さ方向に積層（図２では４層に積層）して構成されたＤＯＩ検出器である。つまり、ＤＯＩ検出器は、各々のシンチレータブロック３１をγ線の深さ方向に積層して構成されたものであり、相互作用を起こした深さ方向と横方向（入射面に平行な方向）との座標情報を重心演算により求める。これにより、相互作用を起こした深さ方向の光源位置(DOI: Depth of Interaction)を弁別することができる。
　図９（ｂ）でも述べたように、図３に示す部分リング型の検出器ユニット３０では、γ線検出器３の一部に抜けがあり、その抜けを通過して検出されないγ線が存在する。したがって、完全投影データを測定することができず、部分リング型の検出器ユニット３０で得られたリストモードデータを投影データに変換すると、欠損領域が生じた不完全投影データが得られる。
　次に、ＧＰＵ１０の具体的な機能について、図４～図５を参照して説明する。図４は、部分リング型ＰＥＴ装置における一連の撮影から再構成画像（断層画像）の出力の流れを示すフローチャートであり、図５は、検出器応答関数（検出確率）の説明に供するγ線検出器での同時計数を示した模式図である。図４ではγ線検出器３として、シンチレータブロック３１のみを図示して、ライトガイド３２やＰＭＴ３３については図示を省略する。
　（ステップＳ１）（リストモードデータの）取得処理
　図３に示す部分リング型の検出器ユニット３０を備えたＰＥＴ装置（部分リング型ＰＥＴ装置）により被検体Ｍ（図１を参照）の撮影を行う。検出器ユニット３０のγ線検出器３はγ線を同時計数したときのみ、シンチレータブロック３１（図２を参照）の位置とγ線の入射タイミングとを有したイベントデータを同時計数回路９（図１を参照）に送り込み、検出器番号や検出時間やγ線のエネルギ等からなるγ線の検出イベント情報を時系列で保存する。これにより、γ線を検出することにより得られた事象データから生成されたリストモードデータを取得する。
　（ステップＳ２）変換処理
　ステップＳ１で得られたリストモードデータを投影データに変換する。上述したように、部分リング型の検出器ユニット３０（図３を参照）で得られるリストモードデータから変換された投影データは不完全投影データである。このステップＳ２は、この発明における変換処理工程に相当する。
　（ステップＳ３）推定処理
　ステップＳ２で変換された投影データに基づいて欠損領域データを推定する。そして、投影データを完全化する。欠損領域データを推定して投影データを完全化する投影完全化アルゴリズムの手法については、特に限定されない。
　例えば、非特許文献２にあるような投影完全化アルゴリズムを利用して投影データを完全化する（非特許文献２のFig.4を参照）。非特許文献２において、不完全投影データをフーリエ変換すると、投影データが完全であれば振幅の小さい信号しか観測されない空間周波数領域に、投影データが不完全であることに起因して振幅の大きい信号が現れる。これを“０”に設定して逆フーリエ変換を行い、欠損領域データを推定する。この演算を繰り返し行うことで、欠損領域データの推定値の精度を向上させる。このステップＳ３は、この発明における推定処理工程に相当する。
　（ステップＳ４）（断層画像の）取得処理
　ステップＳ３で推定された欠損領域データおよびステップＳ１で得られたリストモードデータの両方を用いて、画像を逐次に近似して更新する逐次近似法を行って断層画像を取得する。ここでは、逐次近似法として、ＭＬ－ＥＭ法(Maximum Likelihood Expectation Maximization)を採用する。もちろん、逐次近似法としては、ＭＬ－ＥＭ法に限定されず、ＤＲＡＭＡ法(Dynamic Row-Action Maximum Likelihood Algorithm)でもよいし、スタティックな（つまり静的な）ＲＡＭＬＡ法(Row-Action Maximum Likelihood Algorithm)でもよいし、ＯＳＥＭ法(Ordered Subset ML-EM)でもよい。
　ＭＬ－ＥＭ法の更新式は、後述する下記（３）式で表される。つまり、最適化問題（評価関数の最小化問題）の解を再構成画像（断層画像）の画素値として求めている。従来のリストモード再構成の評価関数をＬ（ｘ）とすると、評価関数Ｌ（ｘ）は、下記（１）式で定義される。
Figure JPOXMLDOC01-appb-M000001
　一方、ハイブリッド再構成法の評価関数をＦ（ｘ）とすると、評価関数Ｆ（ｘ）は、下記（２）式で定義される。また、評価関数Ｆ（ｘ）を構成する評価関数Ｌ（ｘ）については、上記（１）式で定義された評価関数Ｌ（ｘ）を用いている。
Figure JPOXMLDOC01-appb-M000002
　上記（２）式で定義されたハイブリッド再構成法の評価関数Ｆ（ｘ）は、上記（１）式で定義されたリストモード再構成の評価関数Ｌ（ｘ）に、投影完全化アルゴリズムで推定した欠損領域データｐｉ（ｉ＝１，２，…，Ｍｅｓｔ）に関する項（上記（２）式の右辺の第２項）を追加したものである。
　このように、ステップＳ３で推定された欠損領域データｐｉ（ｉ＝１，２，…，Ｍｅｓｔ）およびステップＳ１で得られたリストモードデータの両方を用いて、下記（３）式で表されるＭＬ－ＥＭ法の解を再構成画像（断層画像）の画素値として求める。
Figure JPOXMLDOC01-appb-M000003
　なお、上記（３）式中のｋ（ｋ＝１，２，…）は逐次近似回数、ｘｊ （ｋ）は画素（ｊ＝１，２，…）の画素値である。図５に示すように、検出器応答関数ａｉｊは画素ｊから放出されたγ線フォトンがＬＯＲｉ（ｉ＝１，２，…，Ｍｅｓｔ）に検出される確率（検出確率）である。ここで、ＬＯＲ(Line Of Response)とは、同時計数する２つのγ線検出器３を結ぶ仮想上の直線である。よって、検出器応答関数ａｉｊは、ｉ番目の結晶ペアで検出される確率となる。また、検出器応答関数ｂｉｊは、ｉ番目の仮想的なＬＯＲで検出される確率となる。
　先ず、初期画像であるｘｊ （１）を適宜に設定する。初期画像ｘｊ （１）については、例えば一様な画素値を有する画像であればよく、ｘｊ （１）＞０とする。設定された初期画像ｘｊ （１）を用いて、上記（３）式に繰り返し代入することで、ｘｊ （１），ｘｊ （２），…，ｘｊ （ｋ）が逐次に求められ、ｘｊを順に繰り上げる。反復を表すｋの回数については特に限定されず、適宜に設定すればよい。このように最終的に求められたｘｊをそれに対応する画素ｊごとに並べることで、画像再構成を行い、被検体Ｍの再構成画像（断層画像）を求める。このステップＳ４は、この発明における取得処理工程に相当する。
　図４のステップＳ１～Ｓ４の処理を、上述したようにＧＰＵ１０（図１を参照）が実行する。ＧＰＵ１０で取得された断層画像を、コントローラ５（図１を参照）を介して出力部７（図１を参照）に送り込む。その際に、ＧＰＵ１０で取得された断層画像を、コントローラ５を介してメモリ部８（図１を参照）に書き込んで記憶し、必要に応じてメモリ部８から読み出してもよい。
　図４のステップＳ１～Ｓ４の処理に関する断層画像処理方法によれば、リストモードデータおよび投影データの両方を用いて、画像を逐次に近似して更新する逐次近似法を行っている。すなわち、リストモード再構成および投影データ再構成を組み合わせたハイブリッド方式（ハイブリッド再構成法）を用いて逐次近似法を行っている。よって、投影データによりデータ欠損を補償することができ、従来のリストモード再構成と比較すると、アーティファクトを低減させて画質を改善することができる。
　具体的には、上述のリストモードデータを上述の投影データに変換する変換処理工程（図４のステップＳ２）と、その変換処理工程（ステップＳ２）で変換された投影データに基づいて欠損領域データを推定する推定処理工程（図４のステップＳ３）とを備え、その推定処理工程（ステップＳ３）で推定された上述の欠損領域データを取得処理工程（図４のステップＳ４）で用いられる投影データとして用いる。変換処理工程（ステップＳ２）においてリストモードデータを投影データに変換することにより再構成画像（断層画像）の空間分解能（画質）が劣化するが、ハイブリッド再構成法では補償の対象はリストモードデータであってリストモード再構成が主体である。したがって、一部の欠損領域データのみが投影データの形式をとっているだけなのでデータの質自体は劣化しない。よって、完全化された投影データを用いた投影データ再構成よりも、推定処理工程（ステップＳ３）において推定されて完全化された投影データ（欠損領域データ）を用いたハイブリッド再構成法の方が画質は優れる。
　上述の構成を備えた本実施例に係るＰＥＴ装置によれば、断層画像処理手段（本実施例ではＧＰＵ１０）は上述のハイブリッド再構成法を用いて逐次近似法を行うので、投影データによりデータ欠損を補償することができ、従来のリストモード再構成と比較すると、アーティファクトを低減させて画質を改善することができる。
　本実施例では、放射型断層撮影装置において、陽電子（ポジトロン）の消滅によって発生する複数本の放射線（γ線）を検出して複数個の検出器で放射線（γ線）を同時に検出したときのみ（つまり同時計数したときのみ）被検体Ｍの断層画像を再構成する装置（ＰＥＴ装置）に限定している。また、一部が開口して構成された検出器ユニット３０を備えているので、検出器（本実施例ではγ線検出器３）の一部に抜けがあることにより欠損領域が生じたデータ（リストモードデータおよび投影データ）が得られる。
　この発明は、上記実施形態に限られることはなく、下記のように変形実施することができる。
　（１）上述した実施例では、放射型断層撮影装置として、ＰＥＴ装置（ポジトロン断層撮影装置）を例に採って説明したが、この発明は、単一の放射線を検出して被検体の断層画像を再構成するＳＰＥＣＴ（Single Photon Emission CT）装置などにも適用することができる。また、ＰＥＴ装置とＣＴ装置とを組み合わせたＰＥＴ－ＣＴ装置にも適用することができる。
　（２）上述した実施例では、被検体の全身の断層画像を取得する放射型断層撮影装置（実施例ではＰＥＴ装置）について説明したが、撮影対象については全身に限定されない。被検体の頭部の断層画像を取得する頭部用ＰＥＴ装置や被検体の乳房の断層画像を取得するマンモ用ＰＥＴ装置などに適用してもよい。
　（３）上述した実施例では、３次元に配置された複数のシンチレータ素子からなるＤＯＩ検出器であったが、２次元あるいは３次元に配置された複数のシンチレータ素子からなる放射線検出器にも適用することができる。
　（４）上述した実施例では、固定式の放射型断層撮影装置（実施例ではＰＥＴ装置）について説明したが、他のモダリティ（modality）装置（例えばＣＴ装置）（図示省略）に隣接して設置するために、図６に示すように移動可能な放射型断層撮影装置（ここではＰＥＴ装置）に適用してもよい。図６に示すように、検出器ユニット３０を埋設した筐体は支持アーム３１であり、円弧状にＣ型で形成されている。アーム保持部３２は支持アーム３１を保持し、台車３３に支持されている。台車３３の底部には、後輪３４と前輪３５とを取り付けており、これらを床上で動かすことで、台車３３は搬送可能に構成されている。なお、前輪３５にはモータ（図示省略）が駆動軸（図示省略）を介して連結されており、モータの駆動により前輪３５を駆動させて、術者が台車３３の後ろから任意の方向に押し引きして後輪３４を転がすことで、台車３３を床面上に任意の方向に動かすことができる。
　（５）上述した実施例では、一部が開口して構成された検出器ユニットとして、図３に示す部分リング型の検出器ユニットを例に採って説明したが、図３に示す部分リング型の検出器ユニットに限定されない。上述した図６に示すＣ型の支持アーム３１に埋設されたＣ型の検出器ユニット３０や、マンモ用ＰＥＴ装置において用いられるＣ型やコの字型の検出器ユニットや、複数の検出器が互いに対向して平行に配置された検出器ユニットであってもよい。その他に、マンモ用ＰＥＴ装置において、図７に示すように両脇でそれぞれ挟む検出器ユニット３０に適用してもよい。検出器ユニット３０を埋設した検出板４１は切り欠きとなっており、この切り欠きに脇で挟むことで乳房を検査する。また、γ線検出器３（図７では図示省略）は、この切り欠きに合わせて検出板４１内に複数に並設されている。
　（６）上述した実施例では、図１や図３に示すように被検体Ｍと検出器ユニット３０との距離を変更しなかったが、図８に示すように被検体Ｍと検出器ユニット３０との距離を変更して、検出器ユニット３０を被検体Ｍに近接させてもよい。図８（ａ）は、複数の検出器（実施例ではγ線検出器３）が互いに対向して平行に配置された検出器ユニット３０の構造である。図８（ｂ）は、図３と同様の部分リング型の検出器ユニット３０の構造である。図８（ａ）、図８（ｂ）のいずれの構造であっても、検出器ユニット３０の一部が開口して構成され、かつ各々の検出器ユニット３０が互いに独立して分離されているので、被検体Ｍに検出器ユニット３０を近接させることができる。また、近接させることで放射線の検出効率（感度）が高くなるという効果をも奏する。
　（７）上述した実施例では、断層画像処理手段としてＧＰＵで構成したが、ＧＰＵに限定されない。中央演算処理装置（ＣＰＵ）や、プログラムデータに応じて内部の使用するハードウェア回路（例えば論理回路）が変更可能なプログラマブルデバイス（例えばＦＰＧＡ(Field Programmable Gate Array)）で断層画像処理手段を構成してもよい。
　１０　…　ＧＰＵ
　３０　…　検出器ユニット
　Ｓ２　…　変換処理
　Ｓ３　…　推定処理
　Ｓ４　…　（断層画像の）取得処理
　Ｍ　…　被検体

Claims (3)

  1. 　放射型の断層画像に関する処理を行う断層画像処理方法であって、
    　放射線を検出することにより得られた事象データから生成されたリストモードデータ、および放射線を検出することにより得られた投影データの両方を用いて、画像を逐次に近似して更新する逐次近似法を行って前記断層画像を取得する処理を行う取得処理工程を備えることを特徴とする断層画像処理方法。
  2. 　請求項１に記載の断層画像処理方法において、
    　前記リストモードデータを前記投影データに変換する変換処理工程と、
    　その変換処理工程で変換された前記投影データに基づいて欠損領域データを推定する推定処理工程と
    　を備え、
    　その推定処理工程で推定された前記欠損領域データを前記取得処理工程で用いられる前記投影データとして用いることを特徴とする断層画像処理方法。
  3. 　請求項１または請求項２に記載の断層画像処理方法を用いた放射型断層撮影装置において、
    　前記リストモードデータおよび前記投影データの両方を用いて前記逐次近似法を行って前記断層画像を取得する処理を行う断層画像処理手段を備えることを特徴とする放射型断層撮影装置。
PCT/JP2013/077962 2013-10-15 2013-10-15 断層画像処理方法およびそれを用いた放射型断層撮影装置 WO2015056299A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
PCT/JP2013/077962 WO2015056299A1 (ja) 2013-10-15 2013-10-15 断層画像処理方法およびそれを用いた放射型断層撮影装置
JP2015542427A JP6206501B2 (ja) 2013-10-15 2013-10-15 断層画像処理方法およびそれにおける各工程を実行するための手段を備えた放射線放出型断層撮影装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2013/077962 WO2015056299A1 (ja) 2013-10-15 2013-10-15 断層画像処理方法およびそれを用いた放射型断層撮影装置

Publications (1)

Publication Number Publication Date
WO2015056299A1 true WO2015056299A1 (ja) 2015-04-23

Family

ID=52827773

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2013/077962 WO2015056299A1 (ja) 2013-10-15 2013-10-15 断層画像処理方法およびそれを用いた放射型断層撮影装置

Country Status (2)

Country Link
JP (1) JP6206501B2 (ja)
WO (1) WO2015056299A1 (ja)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004237076A (ja) * 2002-10-04 2004-08-26 Ge Medical Systems Global Technology Co Llc マルチ・モダリティ・イメージング方法及び装置
JP2009544944A (ja) * 2006-07-21 2009-12-17 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Ｔｏｆｐｅｔの再構成を改良するための方法及びシステム
JP2010151653A (ja) * 2008-12-25 2010-07-08 Hamamatsu Photonics Kk 画像処理装置および３次元ｐｅｔ装置
WO2011027402A1 (ja) * 2009-09-04 2011-03-10 株式会社島津製作所 核医学用データ処理方法および核医学診断装置
JP2011075549A (ja) * 2009-10-01 2011-04-14 Toshiba Corp Ｐｅｔ（陽電子放射断層撮影）においてヘリカルスキャン及びリストモード再構成によりサンプリングを強化するシステム及び方法
JP2011153976A (ja) * 2010-01-28 2011-08-11 Shimadzu Corp 断層撮影装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004237076A (ja) * 2002-10-04 2004-08-26 Ge Medical Systems Global Technology Co Llc マルチ・モダリティ・イメージング方法及び装置
JP2009544944A (ja) * 2006-07-21 2009-12-17 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Ｔｏｆｐｅｔの再構成を改良するための方法及びシステム
JP2010151653A (ja) * 2008-12-25 2010-07-08 Hamamatsu Photonics Kk 画像処理装置および３次元ｐｅｔ装置
WO2011027402A1 (ja) * 2009-09-04 2011-03-10 株式会社島津製作所 核医学用データ処理方法および核医学診断装置
JP2011075549A (ja) * 2009-10-01 2011-04-14 Toshiba Corp Ｐｅｔ（陽電子放射断層撮影）においてヘリカルスキャン及びリストモード再構成によりサンプリングを強化するシステム及び方法
JP2011153976A (ja) * 2010-01-28 2011-08-11 Shimadzu Corp 断層撮影装置

Also Published As

Publication number Publication date
JP6206501B2 (ja) 2017-10-04
JPWO2015056299A1 (ja) 2017-03-09

Similar Documents

Publication Publication Date Title
US11517197B2 (en) Apparatus and method for medical image reconstruction using deep learning for computed tomography (CT) image noise and artifacts reduction
US9990741B2 (en) Motion correction in a projection domain in time of flight positron emission tomography
Tong et al. Image reconstruction for PET/CT scanners: past achievements and future challenges
JP2023159080A (ja) 医用画像処理装置及び医用画像処理システム
US9031300B1 (en) System and method reconstructing a nuclear medicine image using deformed attenuation image
EP3067864B1 (en) Iterative reconstruction with enhanced noise control filtering
US8903152B2 (en) Methods and systems for enhanced tomographic imaging
US8457380B2 (en) PET local tomography
US9645261B2 (en) Normalization coefficients in PET continuous bed motion acquisition
JP5152202B2 (ja) ポジトロンｃｔ装置
HU231302B1 (hu) Eljárás és rendszer, valamint tárolóeszköz kép rekonstrukció végrehajtására egy térfogathoz projekciós adathalmaz alapján
Ma et al. An encoder-decoder network for direct image reconstruction on sinograms of a long axial field of view PET
JP2013192951A (ja) 画像再構成装置、画像再構成方法およびｘ線コンピュータ断層撮影装置
US9495770B2 (en) Practical model based CT construction
WO2012120791A1 (ja) ポジトロンｃｔ装置および再構成方法
JP6256608B2 (ja) 画像再構成処理方法
JP2015102516A (ja) 散乱成分推定方法
EP2883084B1 (en) Virtual frames for distributed list-mode time-of-light reconstruction with continuous bed movement
US20220091286A1 (en) Estimation of partially missing attenuation in time-of-flight positron emission tomography
JP2021163493A (ja) データ処理システム及び訓練済み機械学習ベースシステム作成方法
US11164344B2 (en) PET image reconstruction using TOF data and neural network
JP6206501B2 (ja) 断層画像処理方法およびそれにおける各工程を実行するための手段を備えた放射線放出型断層撮影装置
CN111080737B (zh) 图像重建方法、装置及pet扫描系统
JP6052425B2 (ja) 輪郭画像生成装置および核医学診断装置
KR20140130786A (ko) 콘-빔 기반 반응선 재구성을 이용한 초고해상도 ｐｅｔ 영상 재구성 장치 및 그 방법

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

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2015542427

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 13895472

Country of ref document: EP

Kind code of ref document: A1