WO2015056299A1 - 断層画像処理方法およびそれを用いた放射型断層撮影装置 - Google Patents
断層画像処理方法およびそれを用いた放射型断層撮影装置 Download PDFInfo
- 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
Links
- 238000003325 tomography Methods 0.000 title claims description 20
- 238000003672 processing method Methods 0.000 title claims description 14
- 238000012545 processing Methods 0.000 claims abstract description 53
- 238000006243 chemical reaction Methods 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims description 60
- 230000005855 radiation Effects 0.000 claims description 22
- 230000008569 process Effects 0.000 claims description 18
- 230000009467 reduction Effects 0.000 abstract description 3
- 238000002600 positron emission tomography Methods 0.000 description 39
- 230000005251 gamma ray Effects 0.000 description 34
- 238000001514 detection method Methods 0.000 description 16
- 230000006870 function Effects 0.000 description 12
- 238000011156 evaluation Methods 0.000 description 10
- 239000013078 crystal Substances 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000005316 response function Methods 0.000 description 5
- 238000007476 Maximum Likelihood Methods 0.000 description 3
- 210000000481 breast Anatomy 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 229940121896 radiopharmaceutical Drugs 0.000 description 3
- 239000012217 radiopharmaceutical Substances 0.000 description 3
- 230000002799 radiopharmaceutical effect Effects 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000010030 laminating Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
- G01T1/164—Scintigraphy
- G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
- G01T1/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)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Nuclear Medicine (AREA)
Abstract
この発明のPET装置のGPUは、ステップS1~S4の処理を行う。具体的には、リストモードデータを投影データに変換する変換処理工程(ステップS2)と、その変換処理工程(ステップS2)で変換された投影データに基づいて欠損領域データを推定する推定処理工程(ステップS3)とを備え、その推定処理工程(ステップS3)で推定された上述の欠損領域データを取得処理工程(ステップS4)で用いられる投影データとして用いる。変換処理工程(ステップS2)においてリストモードデータを投影データに変換することにより再構成画像(断層画像)の空間分解能(画質)が劣化するが、一部の欠損領域データのみが投影データの形式をとっているだけなのでデータの質自体は劣化しない。よって、投影データによりデータ欠損を補償することができ、従来のリストモード再構成と比較すると、アーティファクトを低減させて画質を改善することができる。
Description
この発明は、放射型の断層画像に関する処理を行う断層画像処理方法およびそれを用いた放射型断層撮影装置に関する。
放射型断層撮影装置として、SPECT装置(Single Photon Emission CT)やPET装置(Positron Emission Tomography)がある。SPECT装置は、単一の放射線(γ線)を検出して被検体の断層画像を再構成する。PET装置は、陽電子(ポジトロン)の消滅によって発生する複数本の放射線(γ線)を検出して複数個の検出器で放射線(γ線)を同時に検出したときのみ(つまり同時計数したときのみ)被検体の断層画像を再構成する。特に、PET装置は「ポジトロン断層撮影装置」とも呼ばれている。以下では、PET装置(ポジトロン断層撮影装置)を例に採って説明する。
図9の下図に示す投影データはサイノグラム(sinogram)であり、横軸は検出器の配置Sを表し、縦軸はγ線の投影角度φである。γ線対を同時計数するので、投影角度φの範囲は0°~360°でなく、その半分の0°~180°となる。
図9(a)に示すように円環状に各々の検出器が並べられた、いわゆる「フルリング型」の検出器ユニットを備えたPET装置では、検出器ユニットが被検体全体を0°~180°でカバーする。よって、それによって得られた投影データは完全な投影データ(完全投影データ)となる。
一方、一部が開口して構成された検出器ユニットにおいて、例えば図9(b)に示すような検出器の一部に抜けがある、いわゆる「部分リング型」の検出器ユニットを備えたPET装置では、その抜けを通過して検出されないγ線が存在する。この場合には、放射性薬剤の体内分布画像を正確に再構成するために必要な最低限のデータ(完全投影データ)を測定することができない。そのように測定されて得られた不完全な投影データ(不完全投影データ)には、図9(b)に示すように欠損領域(図9(b)では符号Dで表記)が生じる。したがって、不完全投影データに対して再構成の計算処理を施すと、再構成画像(断層画像)にデータ欠損に由来した偽像(アーティファクト)が生じる。
なお、図9(b)に示す部分リング型以外にも、***の断層画像を取得するマンモ用PET装置において用いられるC型やコの字型の検出器ユニットや、複数の検出器が互いに対向して平行に配置された検出器ユニット等に例示されるように、一部が開口して構成された検出器ユニットであれば、部分リング型と同様の現象が生じる。すなわち、検出器の一部に抜けがあることによりアーティファクトが生じる。
そこで、アーティファクトを低減させるために、検出時間差(TOF: Time Of Flight)の情報を利用した技術(例えば、特許文献1、非特許文献1参照)や、投影データを完全化する(すなわち欠損領域データを推定する)技術(例えば、非特許文献2参照)がある。検出時間差(TOF)とは、検出されたγ線対の検出時間差である。ポジトロンの対消滅発生地点を、消滅γ線が光速であることを用いて、対消滅発生地点から検出器に到達する時間差を、対消滅発生地点から検出器のシンチレータ素子による光源発生位置までの距離差に換算することで、対消滅発生地点を推定する技術である。なお、非特許文献1は、特許文献1の元となる学術論文である。
特許文献1は、図9(b)に示したような部分リング型TOF-PET装置のシステム構成に関する内容である。検出されたγ線対の検出時間差(TOF)の情報を再構成の計算処理において利用することで、データ欠損による情報量の低下を補償し、再構成画像(断層画像)のアーティファクトを低減させることを目的としている。
特許文献1の元となる学術論文である非特許文献1では、γ線の検出イベント情報(検出器番号、検出時間、γ線のエネルギ等)を時系列で保存した「リストモード」と呼ばれるデータ(「リストモードデータ」または「リストデータ」とも呼ぶ)に対して再構成の計算処理を施している。この再構成方式を「リストモード再構成」という。一方、リストモードデータを検出位置毎にヒストグラム化して得られる投影データに対する再構成方式を「投影データ再構成」という。
非特許文献2は、不完全投影データの欠損部分の投影値を推定するアルゴリズムに関する。このようなアルゴリズムは「投影完全化(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.
しかしながら、上述した特許文献1や非特許文献1や非特許文献2の場合には下記のような問題点がある。
特許文献1や非特許文献1で述べられているように、部分リング型PET装置において、TOFの情報を利用した再構成の計算処理により画像のアーティファクトは低減するが、低減の程度はTOFの情報の精度(検出器の時間分解能)およびデータ欠損量(検出器がない部分の大きさ)に依存する。したがって、TOFの情報を利用することが、満足のいく画質が得られるということには必ずしもならない。つまり、TOFの情報を利用することは、部分リング型PET装置の画質低下の問題に対する根本的な解決策とならない。
欠損領域データを推定して投影データを完全化することができれば、理論上、画像にアーティファクトは発生しない。したがって、部分リング型PET装置における画質低下の問題の根本的な解決策は、非特許文献2にあるような投影完全化アルゴリズムを利用して投影データを完全化することである。
しかし、非特許文献2のような投影データ再構成では、以下のような問題点がある。(1)検出器(シンチレータ結晶)の微細化に伴い結晶数が増加し、結晶数の対(ペア)の分だけ投影データのデータ量が指数関数的に増加してしまう。例えば、結晶数が2倍になると結晶数の対(ペア)の組み合わせや投影データのデータ量は4倍(=22倍)となり、結晶数が3倍になると結晶数の対(ペア)の組み合わせや投影データのデータ量は9倍(=32倍)となる。また、(2)リストモードデータを投影データに変換することが再構成画像(断層画像)の空間分解能の劣化要因となってしまう。一方、リストモードデータのデータ量はシンチレータ結晶数に依存せず、γ線の計数(カウント数)に比例する。
以上の理由により、PET画像の再構成方式ではリストモード再構成が主流になりつつある。その場合には、投影データを作成しないので、当然ながら投影完全化アルゴリズムの出番はない。したがって、従来のリストモード再構成ではデータ欠損は補償されない。
この発明は、このような事情に鑑みてなされたものであって、アーティファクトを低減させて画質を改善することができる断層画像処理方法およびそれを用いた放射型断層撮影装置を提供することを目的とする。
この発明は、このような目的を達成するために、次のような構成をとる。
すなわち、この発明の断層画像処理方法は、放射型の断層画像に関する処理を行う断層画像処理方法であって、放射線を検出することにより得られた事象データから生成されたリストモードデータ、および放射線を検出することにより得られた投影データの両方を用いて、画像を逐次に近似して更新する逐次近似法を行って前記断層画像を取得する処理を行う取得処理工程を備えることを特徴とするものである。
すなわち、この発明の断層画像処理方法は、放射型の断層画像に関する処理を行う断層画像処理方法であって、放射線を検出することにより得られた事象データから生成されたリストモードデータ、および放射線を検出することにより得られた投影データの両方を用いて、画像を逐次に近似して更新する逐次近似法を行って前記断層画像を取得する処理を行う取得処理工程を備えることを特徴とするものである。
この発明の断層画像処理方法によれば、リストモードデータおよび投影データの両方を用いて、画像を逐次に近似して更新する逐次近似法を行っている。すなわち、リストモード再構成および投影データ再構成を組み合わせたハイブリッド方式(以下、「ハイブリッド再構成法」と呼ぶ)を用いて逐次近似法を行っている。よって、投影データによりデータ欠損を補償することができ、従来のリストモード再構成と比較すると、アーティファクトを低減させて画質を改善することができる。
具体的には、上述のリストモードデータを上述の投影データに変換する変換処理工程と、その変換処理工程で変換された投影データに基づいて欠損領域データを推定する推定処理工程とを備え、その推定処理工程で推定された上述の欠損領域データを上述の取得処理工程で用いられる上述の投影データとして用いる。変換処理工程においてリストモードデータを投影データに変換することにより再構成画像(断層画像)の空間分解能(画質)が劣化するが、ハイブリッド再構成法では補償の対象はリストモードデータであってリストモード再構成が主体である。したがって、一部の欠損領域データのみが投影データの形式をとっているだけなのでデータの質自体は劣化しない。よって、完全化された投影データを用いた投影データ再構成よりも、推定処理工程において推定されて完全化された投影データ(欠損領域データ)を用いたハイブリッド再構成法の方が画質は優れる。
また、この発明の放射型断層撮影装置は、上述したこれらの発明の断層画像処理方法を用いた放射型断層撮影装置において、上述のリストモードデータおよび上述の投影データの両方を用いて逐次近似法を行って断層画像を取得する処理を行う断層画像処理手段を備えることを特徴とするものである。
この発明の放射型断層撮影装置によれば、断層画像処理手段は上述のハイブリッド再構成法を用いて逐次近似法を行うので、投影データによりデータ欠損を補償することができ、従来のリストモード再構成と比較すると、アーティファクトを低減させて画質を改善することができる。
この発明に係る断層画像処理方法およびそれを用いた放射型断層撮影装置によれば、(リストモード再構成および投影データ再構成を組み合わせた)ハイブリッド再構成法を用いて逐次近似法を行うので、投影データによりデータ欠損を補償することができ、従来のリストモード再構成と比較すると、アーティファクトを低減させて画質を改善することができる。
以下、図面を参照してこの発明の実施例を説明する。図1は、実施例に係るPET(Positron Emission Tomography)装置の側面図およびブロック図であり、図2は、γ線検出器の概略斜視図であり、図3は、部分リング型の検出器ユニットの概略正面図である。本実施例では、放射型断層撮影装置として、PET装置(ポジトロン断層撮影装置)を例に採って説明するとともに、一部が開口して構成された検出器ユニットとして、図3(図9(b)も参照)に示す部分リング型の検出器ユニットを例に採って説明する。
本実施例に係るPET装置の傍らには、図1に示すように、被検体Mを載置する天板1が配置されている。この天板1は、上下に昇降移動、被検体Mの体軸Zに沿って平行移動するように構成されている。このように構成することで、天板1に載置された被検体Mは、後述するガントリ2の開口部2aを通って、頭部から順に腹部、足部へと走査されて、被検体Mの画像を得る。なお、走査される部位や各部位の走査順序については特に限定されない。なお、本実施例に係るPET装置は、天板1を構成として含ませていないが、天板1を構成として含ませて備えてもよい。
本実施例に係るPET装置は、開口部2aを有したガントリ2と、γ線検出器3とを備えている。γ線検出器3は、被検体Mの体軸Z周りを取り囲むようにして部分リング状に配置されており、ガントリ2内に埋設されている。本実施例では、図3に示すように一部が開口して構成された部分リング型の検出器ユニット30となるように、各々のγ線検出器3が並べられている。γ線検出器3を構成する検出器ユニット30は、この発明における検出器ユニットに相当する。
その他にも、本実施例に係るPET装置は、天板駆動部4とコントローラ5と入力部6と出力部7とメモリ部8と同時計数回路9とGPU(Graphics Processing Unit)10とを備えている。天板駆動部6は、天板1の上述した移動を行うように駆動する機構であって、図示を省略するモータなどで構成されている。GPU10は、この発明における断層画像処理手段に相当する。
コントローラ5は、本実施例に係るPET装置を構成する各部分を統括制御する。コントローラ5は、中央演算処理装置(CPU)などで構成されている。
入力部6は、オペレータが入力したデータや命令をコントローラ5に送り込む。入力部6は、マウスやキーボードやジョイスティックやトラックボールやタッチパネルなどに代表されるポインティングデバイスで構成されている。出力部7はモニタなどに代表される表示部やプリンタなどで構成されている。
メモリ部8は、ROM(Read-only Memory)やRAM(Random-Access Memory)などに代表される記憶媒体で構成されている。本実施例では、同時計数回路9で同時計数された計数値(カウント)や同時計数した2つのγ線検出器3からなる検出器対やLORといった同時計数に関するデータや、GPU10で演算処理された各種のデータなどについてはRAMに書き込んで記憶し、必要に応じてRAMから読み出す。ROMには、各種の制御処理(例えば天板1の駆動制御)や画像処理(例えば断層画像処理)を行うためのプログラム等を予め記憶しており、そのプログラムをコントローラ5およびGPU10が実行することでそのプログラムに応じた制御処理や画像処理をそれぞれ行う。特に、本実施例では、後述する逐次近似法を行って断層画像を取得する処理を行う取得処理,リストモードデータを投影データに変換する変換処理,変換された投影データに基づいて欠損領域データを推定する推定処理に関するプログラムをROMに予め記憶しており、それらの取得処理や変換処理や推定処理に関するプログラムをGPU10が実行することで後述するステップS1~S4の処理を実行する。
放射性薬剤が投与された被検体Mから発生したγ線をγ線検出器3のシンチレータブロック31(図2を参照)が光に変換して、変換されたその光をγ線検出器3の光電子増倍管(PMT: Photo Multiplier Tube)33(図2を参照)は増倍させて電気信号に変換する。その電気信号をイベントとして同時計数回路9に送り込む。
具体的には、被検体Mに放射性薬剤を投与すると、ポジトロン放出型のRIのポジトロンが消滅することにより、2本のγ線が発生する。同時計数回路9は、シンチレータブロック31(図2を参照)の位置とγ線の入射タイミングとをチェックし、被検体Mの両側にある2つのシンチレータブロック31でγ線が同時に入射したときのみ、送り込まれたイベントを適正なデータと判定する。一方のシンチレータブロック31のみにγ線が入射したときには、同時計数回路9は棄却する。つまり、同時計数回路9は、上述した電気信号に基づいて、2つのγ線検出器3においてγ線が同時観測されたことを検出する。
同時計数回路9に送り込まれたイベントを、GPU10に送り込む。GPU10は取得処理や変換処理や推定処理による画像再構成を行って、被検体Mの断層画像を求める。断層画像を、コントローラ5を介して出力部7に送り込む。このようにして、GPU10で得られた断層画像に基づいて断層撮影を行う。GPU10の具体的な機能については後述する。
γ線検出器3は、図2に示すようにシンチレータブロック31と、そのシンチレータブロック31に対して光学的に結合されたライトガイド32と、そのライトガイド32に対して光学的に結合された光電子増倍管(以下、単に「PMT」と略記する)33とを備えている。シンチレータブロック31を構成する各シンチレータ素子は、γ線の入射に伴って発光することでγ線から光に変換する。この変換によってシンチレータ素子はγ線を検出する。シンチレータ素子において発光した光がシンチレータブロック31で十分に拡散されて、ライトガイド32を介してPMT33に入力される。PMT33は、シンチレータブロック31で変換された光を増倍させて電気信号に変換する。その電気信号は、上述したようにイベントとして同時計数回路9(図1を参照)に送り込まれる。
また、図2に示すγ線検出器3は、各々のシンチレータブロック31をγ線の深さ方向に積層(図2では4層に積層)して構成されたDOI検出器である。つまり、DOI検出器は、各々のシンチレータブロック31をγ線の深さ方向に積層して構成されたものであり、相互作用を起こした深さ方向と横方向(入射面に平行な方向)との座標情報を重心演算により求める。これにより、相互作用を起こした深さ方向の光源位置(DOI: Depth of Interaction)を弁別することができる。
図9(b)でも述べたように、図3に示す部分リング型の検出器ユニット30では、γ線検出器3の一部に抜けがあり、その抜けを通過して検出されないγ線が存在する。したがって、完全投影データを測定することができず、部分リング型の検出器ユニット30で得られたリストモードデータを投影データに変換すると、欠損領域が生じた不完全投影データが得られる。
次に、GPU10の具体的な機能について、図4~図5を参照して説明する。図4は、部分リング型PET装置における一連の撮影から再構成画像(断層画像)の出力の流れを示すフローチャートであり、図5は、検出器応答関数(検出確率)の説明に供するγ線検出器での同時計数を示した模式図である。図4ではγ線検出器3として、シンチレータブロック31のみを図示して、ライトガイド32やPMT33については図示を省略する。
(ステップS1)(リストモードデータの)取得処理
図3に示す部分リング型の検出器ユニット30を備えたPET装置(部分リング型PET装置)により被検体M(図1を参照)の撮影を行う。検出器ユニット30のγ線検出器3はγ線を同時計数したときのみ、シンチレータブロック31(図2を参照)の位置とγ線の入射タイミングとを有したイベントデータを同時計数回路9(図1を参照)に送り込み、検出器番号や検出時間やγ線のエネルギ等からなるγ線の検出イベント情報を時系列で保存する。これにより、γ線を検出することにより得られた事象データから生成されたリストモードデータを取得する。
図3に示す部分リング型の検出器ユニット30を備えたPET装置(部分リング型PET装置)により被検体M(図1を参照)の撮影を行う。検出器ユニット30のγ線検出器3はγ線を同時計数したときのみ、シンチレータブロック31(図2を参照)の位置とγ線の入射タイミングとを有したイベントデータを同時計数回路9(図1を参照)に送り込み、検出器番号や検出時間やγ線のエネルギ等からなるγ線の検出イベント情報を時系列で保存する。これにより、γ線を検出することにより得られた事象データから生成されたリストモードデータを取得する。
(ステップS2)変換処理
ステップS1で得られたリストモードデータを投影データに変換する。上述したように、部分リング型の検出器ユニット30(図3を参照)で得られるリストモードデータから変換された投影データは不完全投影データである。このステップS2は、この発明における変換処理工程に相当する。
ステップS1で得られたリストモードデータを投影データに変換する。上述したように、部分リング型の検出器ユニット30(図3を参照)で得られるリストモードデータから変換された投影データは不完全投影データである。このステップS2は、この発明における変換処理工程に相当する。
(ステップS3)推定処理
ステップS2で変換された投影データに基づいて欠損領域データを推定する。そして、投影データを完全化する。欠損領域データを推定して投影データを完全化する投影完全化アルゴリズムの手法については、特に限定されない。
ステップS2で変換された投影データに基づいて欠損領域データを推定する。そして、投影データを完全化する。欠損領域データを推定して投影データを完全化する投影完全化アルゴリズムの手法については、特に限定されない。
例えば、非特許文献2にあるような投影完全化アルゴリズムを利用して投影データを完全化する(非特許文献2のFig.4を参照)。非特許文献2において、不完全投影データをフーリエ変換すると、投影データが完全であれば振幅の小さい信号しか観測されない空間周波数領域に、投影データが不完全であることに起因して振幅の大きい信号が現れる。これを“0”に設定して逆フーリエ変換を行い、欠損領域データを推定する。この演算を繰り返し行うことで、欠損領域データの推定値の精度を向上させる。このステップS3は、この発明における推定処理工程に相当する。
(ステップS4)(断層画像の)取得処理
ステップS3で推定された欠損領域データおよびステップS1で得られたリストモードデータの両方を用いて、画像を逐次に近似して更新する逐次近似法を行って断層画像を取得する。ここでは、逐次近似法として、ML-EM法(Maximum Likelihood Expectation Maximization)を採用する。もちろん、逐次近似法としては、ML-EM法に限定されず、DRAMA法(Dynamic Row-Action Maximum Likelihood Algorithm)でもよいし、スタティックな(つまり静的な)RAMLA法(Row-Action Maximum Likelihood Algorithm)でもよいし、OSEM法(Ordered Subset ML-EM)でもよい。
ステップS3で推定された欠損領域データおよびステップS1で得られたリストモードデータの両方を用いて、画像を逐次に近似して更新する逐次近似法を行って断層画像を取得する。ここでは、逐次近似法として、ML-EM法(Maximum Likelihood Expectation Maximization)を採用する。もちろん、逐次近似法としては、ML-EM法に限定されず、DRAMA法(Dynamic Row-Action Maximum Likelihood Algorithm)でもよいし、スタティックな(つまり静的な)RAMLA法(Row-Action Maximum Likelihood Algorithm)でもよいし、OSEM法(Ordered Subset ML-EM)でもよい。
ML-EM法の更新式は、後述する下記(3)式で表される。つまり、最適化問題(評価関数の最小化問題)の解を再構成画像(断層画像)の画素値として求めている。従来のリストモード再構成の評価関数をL(x)とすると、評価関数L(x)は、下記(1)式で定義される。
一方、ハイブリッド再構成法の評価関数をF(x)とすると、評価関数F(x)は、下記(2)式で定義される。また、評価関数F(x)を構成する評価関数L(x)については、上記(1)式で定義された評価関数L(x)を用いている。
上記(2)式で定義されたハイブリッド再構成法の評価関数F(x)は、上記(1)式で定義されたリストモード再構成の評価関数L(x)に、投影完全化アルゴリズムで推定した欠損領域データpi(i=1,2,…,Mest)に関する項(上記(2)式の右辺の第2項)を追加したものである。
このように、ステップS3で推定された欠損領域データpi(i=1,2,…,Mest)およびステップS1で得られたリストモードデータの両方を用いて、下記(3)式で表されるML-EM法の解を再構成画像(断層画像)の画素値として求める。
なお、上記(3)式中のk(k=1,2,…)は逐次近似回数、xj
(k)は画素(j=1,2,…)の画素値である。図5に示すように、検出器応答関数aijは画素jから放出されたγ線フォトンがLORi(i=1,2,…,Mest)に検出される確率(検出確率)である。ここで、LOR(Line Of Response)とは、同時計数する2つのγ線検出器3を結ぶ仮想上の直線である。よって、検出器応答関数aijは、i番目の結晶ペアで検出される確率となる。また、検出器応答関数bijは、i番目の仮想的なLORで検出される確率となる。
先ず、初期画像であるxj
(1)を適宜に設定する。初期画像xj
(1)については、例えば一様な画素値を有する画像であればよく、xj
(1)>0とする。設定された初期画像xj
(1)を用いて、上記(3)式に繰り返し代入することで、xj
(1),xj
(2),…,xj
(k)が逐次に求められ、xjを順に繰り上げる。反復を表すkの回数については特に限定されず、適宜に設定すればよい。このように最終的に求められたxjをそれに対応する画素jごとに並べることで、画像再構成を行い、被検体Mの再構成画像(断層画像)を求める。このステップS4は、この発明における取得処理工程に相当する。
図4のステップS1~S4の処理を、上述したようにGPU10(図1を参照)が実行する。GPU10で取得された断層画像を、コントローラ5(図1を参照)を介して出力部7(図1を参照)に送り込む。その際に、GPU10で取得された断層画像を、コントローラ5を介してメモリ部8(図1を参照)に書き込んで記憶し、必要に応じてメモリ部8から読み出してもよい。
図4のステップS1~S4の処理に関する断層画像処理方法によれば、リストモードデータおよび投影データの両方を用いて、画像を逐次に近似して更新する逐次近似法を行っている。すなわち、リストモード再構成および投影データ再構成を組み合わせたハイブリッド方式(ハイブリッド再構成法)を用いて逐次近似法を行っている。よって、投影データによりデータ欠損を補償することができ、従来のリストモード再構成と比較すると、アーティファクトを低減させて画質を改善することができる。
具体的には、上述のリストモードデータを上述の投影データに変換する変換処理工程(図4のステップS2)と、その変換処理工程(ステップS2)で変換された投影データに基づいて欠損領域データを推定する推定処理工程(図4のステップS3)とを備え、その推定処理工程(ステップS3)で推定された上述の欠損領域データを取得処理工程(図4のステップS4)で用いられる投影データとして用いる。変換処理工程(ステップS2)においてリストモードデータを投影データに変換することにより再構成画像(断層画像)の空間分解能(画質)が劣化するが、ハイブリッド再構成法では補償の対象はリストモードデータであってリストモード再構成が主体である。したがって、一部の欠損領域データのみが投影データの形式をとっているだけなのでデータの質自体は劣化しない。よって、完全化された投影データを用いた投影データ再構成よりも、推定処理工程(ステップS3)において推定されて完全化された投影データ(欠損領域データ)を用いたハイブリッド再構成法の方が画質は優れる。
上述の構成を備えた本実施例に係るPET装置によれば、断層画像処理手段(本実施例ではGPU10)は上述のハイブリッド再構成法を用いて逐次近似法を行うので、投影データによりデータ欠損を補償することができ、従来のリストモード再構成と比較すると、アーティファクトを低減させて画質を改善することができる。
本実施例では、放射型断層撮影装置において、陽電子(ポジトロン)の消滅によって発生する複数本の放射線(γ線)を検出して複数個の検出器で放射線(γ線)を同時に検出したときのみ(つまり同時計数したときのみ)被検体Mの断層画像を再構成する装置(PET装置)に限定している。また、一部が開口して構成された検出器ユニット30を備えているので、検出器(本実施例ではγ線検出器3)の一部に抜けがあることにより欠損領域が生じたデータ(リストモードデータおよび投影データ)が得られる。
この発明は、上記実施形態に限られることはなく、下記のように変形実施することができる。
(1)上述した実施例では、放射型断層撮影装置として、PET装置(ポジトロン断層撮影装置)を例に採って説明したが、この発明は、単一の放射線を検出して被検体の断層画像を再構成するSPECT(Single Photon Emission CT)装置などにも適用することができる。また、PET装置とCT装置とを組み合わせたPET-CT装置にも適用することができる。
(2)上述した実施例では、被検体の全身の断層画像を取得する放射型断層撮影装置(実施例ではPET装置)について説明したが、撮影対象については全身に限定されない。被検体の頭部の断層画像を取得する頭部用PET装置や被検体の***の断層画像を取得するマンモ用PET装置などに適用してもよい。
(3)上述した実施例では、3次元に配置された複数のシンチレータ素子からなるDOI検出器であったが、2次元あるいは3次元に配置された複数のシンチレータ素子からなる放射線検出器にも適用することができる。
(4)上述した実施例では、固定式の放射型断層撮影装置(実施例ではPET装置)について説明したが、他のモダリティ(modality)装置(例えばCT装置)(図示省略)に隣接して設置するために、図6に示すように移動可能な放射型断層撮影装置(ここではPET装置)に適用してもよい。図6に示すように、検出器ユニット30を埋設した筐体は支持アーム31であり、円弧状にC型で形成されている。アーム保持部32は支持アーム31を保持し、台車33に支持されている。台車33の底部には、後輪34と前輪35とを取り付けており、これらを床上で動かすことで、台車33は搬送可能に構成されている。なお、前輪35にはモータ(図示省略)が駆動軸(図示省略)を介して連結されており、モータの駆動により前輪35を駆動させて、術者が台車33の後ろから任意の方向に押し引きして後輪34を転がすことで、台車33を床面上に任意の方向に動かすことができる。
(5)上述した実施例では、一部が開口して構成された検出器ユニットとして、図3に示す部分リング型の検出器ユニットを例に採って説明したが、図3に示す部分リング型の検出器ユニットに限定されない。上述した図6に示すC型の支持アーム31に埋設されたC型の検出器ユニット30や、マンモ用PET装置において用いられるC型やコの字型の検出器ユニットや、複数の検出器が互いに対向して平行に配置された検出器ユニットであってもよい。その他に、マンモ用PET装置において、図7に示すように両脇でそれぞれ挟む検出器ユニット30に適用してもよい。検出器ユニット30を埋設した検出板41は切り欠きとなっており、この切り欠きに脇で挟むことで***を検査する。また、γ線検出器3(図7では図示省略)は、この切り欠きに合わせて検出板41内に複数に並設されている。
(6)上述した実施例では、図1や図3に示すように被検体Mと検出器ユニット30との距離を変更しなかったが、図8に示すように被検体Mと検出器ユニット30との距離を変更して、検出器ユニット30を被検体Mに近接させてもよい。図8(a)は、複数の検出器(実施例ではγ線検出器3)が互いに対向して平行に配置された検出器ユニット30の構造である。図8(b)は、図3と同様の部分リング型の検出器ユニット30の構造である。図8(a)、図8(b)のいずれの構造であっても、検出器ユニット30の一部が開口して構成され、かつ各々の検出器ユニット30が互いに独立して分離されているので、被検体Mに検出器ユニット30を近接させることができる。また、近接させることで放射線の検出効率(感度)が高くなるという効果をも奏する。
(7)上述した実施例では、断層画像処理手段としてGPUで構成したが、GPUに限定されない。中央演算処理装置(CPU)や、プログラムデータに応じて内部の使用するハードウェア回路(例えば論理回路)が変更可能なプログラマブルデバイス(例えばFPGA(Field Programmable Gate Array))で断層画像処理手段を構成してもよい。
10 … GPU
30 … 検出器ユニット
S2 … 変換処理
S3 … 推定処理
S4 … (断層画像の)取得処理
M … 被検体
30 … 検出器ユニット
S2 … 変換処理
S3 … 推定処理
S4 … (断層画像の)取得処理
M … 被検体
Claims (3)
- 放射型の断層画像に関する処理を行う断層画像処理方法であって、
放射線を検出することにより得られた事象データから生成されたリストモードデータ、および放射線を検出することにより得られた投影データの両方を用いて、画像を逐次に近似して更新する逐次近似法を行って前記断層画像を取得する処理を行う取得処理工程を備えることを特徴とする断層画像処理方法。 - 請求項1に記載の断層画像処理方法において、
前記リストモードデータを前記投影データに変換する変換処理工程と、
その変換処理工程で変換された前記投影データに基づいて欠損領域データを推定する推定処理工程と
を備え、
その推定処理工程で推定された前記欠損領域データを前記取得処理工程で用いられる前記投影データとして用いることを特徴とする断層画像処理方法。 - 請求項1または請求項2に記載の断層画像処理方法を用いた放射型断層撮影装置において、
前記リストモードデータおよび前記投影データの両方を用いて前記逐次近似法を行って前記断層画像を取得する処理を行う断層画像処理手段を備えることを特徴とする放射型断層撮影装置。
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)
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 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Tofpetの再構成を改良するための方法及びシステム |
JP2010151653A (ja) * | 2008-12-25 | 2010-07-08 | Hamamatsu Photonics Kk | 画像処理装置および3次元pet装置 |
WO2011027402A1 (ja) * | 2009-09-04 | 2011-03-10 | 株式会社島津製作所 | 核医学用データ処理方法および核医学診断装置 |
JP2011075549A (ja) * | 2009-10-01 | 2011-04-14 | Toshiba Corp | Pet(陽電子放射断層撮影)においてヘリカルスキャン及びリストモード再構成によりサンプリングを強化するシステム及び方法 |
JP2011153976A (ja) * | 2010-01-28 | 2011-08-11 | Shimadzu Corp | 断層撮影装置 |
-
2013
- 2013-10-15 WO PCT/JP2013/077962 patent/WO2015056299A1/ja active Application Filing
- 2013-10-15 JP JP2015542427A patent/JP6206501B2/ja active Active
Patent Citations (6)
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 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Tofpetの再構成を改良するための方法及びシステム |
JP2010151653A (ja) * | 2008-12-25 | 2010-07-08 | Hamamatsu Photonics Kk | 画像処理装置および3次元pet装置 |
WO2011027402A1 (ja) * | 2009-09-04 | 2011-03-10 | 株式会社島津製作所 | 核医学用データ処理方法および核医学診断装置 |
JP2011075549A (ja) * | 2009-10-01 | 2011-04-14 | Toshiba Corp | Pet(陽電子放射断層撮影)においてヘリカルスキャン及びリストモード再構成によりサンプリングを強化するシステム及び方法 |
JP2011153976A (ja) * | 2010-01-28 | 2011-08-11 | Shimadzu Corp | 断層撮影装置 |
Also Published As
Publication number | Publication date |
---|---|
JPWO2015056299A1 (ja) | 2017-03-09 |
JP6206501B2 (ja) | 2017-10-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20230119427A1 (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 | |
US9645261B2 (en) | Normalization coefficients in PET continuous bed motion acquisition | |
JP5152202B2 (ja) | ポジトロンct装置 | |
Ma et al. | An encoder-decoder network for direct image reconstruction on sinograms of a long axial field of view PET | |
JP2021163493A (ja) | データ処理システム及び訓練済み機械学習ベースシステム作成方法 | |
JP2013192951A (ja) | 画像再構成装置、画像再構成方法およびx線コンピュータ断層撮影装置 | |
US9495770B2 (en) | Practical model based CT construction | |
JP6256608B2 (ja) | 画像再構成処理方法 | |
JP2015102516A (ja) | 散乱成分推定方法 | |
EP2883084B1 (en) | Virtual frames for distributed list-mode time-of-light reconstruction with continuous bed movement | |
US11874411B2 (en) | Estimation of partially missing attenuation in time-of-flight positron emission tomography | |
US11164344B2 (en) | PET image reconstruction using TOF data and neural network | |
JP6206501B2 (ja) | 断層画像処理方法およびそれにおける各工程を実行するための手段を備えた放射線放出型断層撮影装置 | |
US20220148236A1 (en) | Systems and methods for motion estimation in pet imaging using ai image reconstructions | |
CN111080737B (zh) | 图像重建方法、装置及pet扫描*** | |
JP6052425B2 (ja) | 輪郭画像生成装置および核医学診断装置 | |
KR20140130786A (ko) | 콘-빔 기반 반응선 재구성을 이용한 초고해상도 pet 영상 재구성 장치 및 그 방법 | |
KR101958098B1 (ko) | 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법 | |
KR101356881B1 (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 |