JP2013215472A - Radiation tomographic image generating apparatus, radiation tomographic apparatus, and radiation tomographic image generating program - Google Patents
Radiation tomographic image generating apparatus, radiation tomographic apparatus, and radiation tomographic image generating program Download PDFInfo
- Publication number
- JP2013215472A JP2013215472A JP2012090264A JP2012090264A JP2013215472A JP 2013215472 A JP2013215472 A JP 2013215472A JP 2012090264 A JP2012090264 A JP 2012090264A JP 2012090264 A JP2012090264 A JP 2012090264A JP 2013215472 A JP2013215472 A JP 2013215472A
- Authority
- JP
- Japan
- Prior art keywords
- image
- projection data
- measured
- radiation tomographic
- tomographic image
- Prior art date
- Legal status (The legal status 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 status listed.)
- Granted
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 111
- 230000005540 biological transmission Effects 0.000 claims description 120
- 238000010521 absorption reaction Methods 0.000 claims description 115
- 238000012937 correction Methods 0.000 claims description 82
- 238000000034 method Methods 0.000 claims description 57
- 238000013480 data collection Methods 0.000 claims description 34
- 238000003325 tomography Methods 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 7
- 238000005457 optimization Methods 0.000 abstract description 10
- 238000002603 single-photon emission computed tomography Methods 0.000 description 55
- 230000005251 gamma ray Effects 0.000 description 33
- 238000001514 detection method Methods 0.000 description 29
- 238000010586 diagram Methods 0.000 description 17
- 238000006243 chemical reaction Methods 0.000 description 14
- 239000002131 composite material Substances 0.000 description 7
- 238000005259 measurement Methods 0.000 description 7
- 238000012986 modification Methods 0.000 description 7
- 230000004048 modification Effects 0.000 description 7
- 230000000007 visual effect Effects 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 230000001678 irradiating effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 229940121896 radiopharmaceutical Drugs 0.000 description 2
- 239000012217 radiopharmaceutical Substances 0.000 description 2
- 230000002799 radiopharmaceutical effect Effects 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
Images
Landscapes
- Nuclear Medicine (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
本発明は、医療用および異物検査等の産業用に使用され、放射線(X線およびγ線等)を検出して投影データを収集し、収集した投影データを画像再構成して放射線断層画像を生成する放射線断層画像生成装置、放射線断層撮影装置および放射線断層画像生成プログラムに関する。 INDUSTRIAL APPLICABILITY The present invention is used for medical purposes and industrial purposes such as foreign body inspection, detects radiation (X-rays, γ-rays, etc.), collects projection data, and reconstructs the collected projection data to reconstruct a radiation tomographic image. The present invention relates to a radiation tomographic image generation apparatus, a radiation tomography apparatus, and a radiation tomographic image generation program.
従来、この種の放射線断層撮影装置として、PET装置、SPECT装置、X線CT装置、PET/CT装置およびSPECT/CT装置などがある(例えば、特許文献1参照)。これらの装置は、測定によって収集された図14(a)に示す実測投影データ(サイノグラム)を画像再構成することで被検体の断層画像(図14(b))を取得している。 Conventionally, as this type of radiation tomography apparatus, there are a PET apparatus, a SPECT apparatus, an X-ray CT apparatus, a PET / CT apparatus, a SPECT / CT apparatus, and the like (for example, see Patent Document 1). These apparatuses acquire a tomographic image (FIG. 14B) of a subject by reconstructing an image of actual projection data (sinogram) shown in FIG. 14A collected by measurement.
しかしながら、被検体よりも撮影断面の有効視野(FOV)が狭い場合には、被検体が有効視野をはみ出して欠ける、いわゆるトランケーションが生じる。なお、はみ出した部分をトランケーション部分とする。図14(c)は、トランケーションが生じた実測投影データ(サイノグラム)を示す図である。図14(c)中の破線とサイノグラムで囲まれた領域が符号200に示すトランケーション部分である。このようなトランケーションした(被検体が欠けた)実測投影データをそのままの状態で画像再構成(例えばML−EM法)すると、図14(d)に示すように、生成された断層画像(再構成画像)の有効視野内の外端であって被検体がはみ出した領域付近にトランケーションアーチファクトが生じるという問題がある。 However, when the effective field of view (FOV) of the imaging section is narrower than the subject, so-called truncation occurs in which the subject protrudes from the effective field of view and is missing. Note that the protruding part is the truncation part. FIG. 14C is a diagram showing actually measured projection data (sinogram) in which truncation has occurred. A region surrounded by a broken line and a sinogram in FIG. When image reconstruction (for example, the ML-EM method) is performed on the truncated projection data (without the subject) as it is, the generated tomographic image (reconstruction) is obtained as shown in FIG. There is a problem that truncation artifacts occur near the outer edge of the image) in the effective field of view and in the vicinity of the area where the subject protrudes.
さらに、例えば、PET/CT装置やSPECT/CT装置(CTはトランスミッションCTを含む)等の複数の撮影システム(モダリティ)で構成される放射線断層撮影装置においては、一般的にPET装置やSPECT装置よりもCT装置の撮影の有効視野が狭いことが知られている(例えば特許文献1、並びに非特許文献1および2参照)。そのため、CT画像をPET装置やSPECT装置の吸収補正に用いようとする場合には、X線CT画像に生じたトランケーションアーチファクトの影響が吸収補正後に画像再構成されたPET画像にも与えてしまうこととなる。
Further, for example, in a radiation tomography apparatus including a plurality of imaging systems (modalities) such as a PET / CT apparatus and a SPECT / CT apparatus (CT includes a transmission CT), the PET apparatus and the SPECT apparatus are generally used. It is also known that the effective field of view of the CT apparatus is narrow (see, for example,
しかしながら、このような構成を有する従来例の場合には、次のような問題がある。
すなわち、トランケーションアーチファクトの発生は、トランケーション部分を推定することにより防止している。例えば特許文献1では、SPECT画像から被検体の体輪郭を抽出し、抽出した体輪郭を用いて楕円近似によりトランスミッションCT投影データの総和値およびトランスミッションCT像の重心からトランケーション部分を推定し、トランスミッションCT投影データに2次式外挿することでトランケーション補正を行っている。この手法は、実用化されているものの、被検体境界位置の推定などにパラメータの最適化が必要であり、アルゴリズム的に煩雑さを含むものである。
However, the conventional example having such a configuration has the following problems.
That is, the generation of truncation artifacts is prevented by estimating the truncation part. For example, in
本発明は、このような事情に鑑みてなされたものであって、パラメータの最適化をほとんど必要としなく、トランケーション部分の画像推定が行われた放射線断層画像を容易に取得することが可能な放射線断層画像生成装置、放射線断層撮影装置および放射線断層画像生成プログラムを提供することを目的とする。 The present invention has been made in view of such circumstances, and it is possible to easily acquire a radiation tomographic image in which image estimation of a truncation portion is performed with little need for parameter optimization. An object is to provide a tomographic image generation apparatus, a radiation tomography apparatus, and a radiation tomographic image generation program.
本発明は、このような目的を達成するために、次のような構成をとる。すなわち、本発明に係る放射線断層画像生成装置は、被検体に対して異なる複数の方向から収集された実測投影データに基づいて放射線断層画像を生成する放射線断層画像生成装置において、前記実測投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記実測投影データのみを用いて逐次近似画像再構成を行う画像再構成部を備えていることを特徴とするものである。 In order to achieve such an object, the present invention has the following configuration. That is, the radiation tomographic image generation device according to the present invention is a radiation tomographic image generation device that generates a radiation tomographic image based on measured projection data collected from a plurality of different directions with respect to a subject. A reconstructed image area including a circular area set wider than a circular area having a data width as a diameter, and an image reconstructing unit that performs successive approximate image reconstruction using only the measured projection data. It is what.
本発明に係る放射線断層画像生成装置によれば、画像再構成部は、実測投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、実測投影データのみを用いて逐次近似画像再構成を行っている。逐次近似画像再構成の再構成画像領域を広く設定することにより、実測投影データが、被検体Mが検出範囲からはみ出してトランケーションしている場合でも揃っている実測投影データのみでトランケーション部分の画像を推定することができる。揃っている実測投影データは、放射線断層画像を生成するための投影データに推定された部分を有することによる誤差が生じないので、投影データ間の整合性が高くアーチファクトを低減させることができる。したがって、パラメータの最適化などの設定をほとんど必要としなく、トランケーション部分の画像推定が行われた放射線断層画像を容易に取得することができる。 According to the radiation tomographic image generation device according to the present invention, the image reconstruction unit is a reconstructed image area including a circular area set wider than a circular area having the diameter of the data width of the actual measurement projection data. Is used to perform successive approximation image reconstruction. By setting a wide range of reconstructed image areas for successive approximation image reconstruction, the image of the truncation portion can be obtained with only the measured projection data that is obtained even when the measured projection data is truncated even when the subject M is out of the detection range. Can be estimated. Since the measured projection data that has been prepared does not cause an error due to having an estimated portion in the projection data for generating the radiation tomographic image, the consistency between the projection data is high and artifacts can be reduced. Therefore, it is possible to easily obtain a radiation tomographic image in which the image estimation of the truncation portion is performed with little setting such as parameter optimization.
また、本発明に係る放射線断層画像生成装置において、前記画像再構成部は、前記広く設定された円形領域を含む再構成画像領域の初期画像に、天板に対応する画素値を吸収係数とする先見情報を与えることが好ましい。これにより、逐次近似画像再構成の際に、天板に対応する画素値の精度を高めることができ、また、これに伴いトランケーション部分の画素値の精度を高めることができる。 Further, in the radiation tomographic image generation device according to the present invention, the image reconstruction unit uses, as an absorption coefficient, a pixel value corresponding to a top plate in the initial image of the reconstructed image region including the widely set circular region. It is preferable to provide foresight information. Thereby, the accuracy of the pixel value corresponding to the top plate can be increased during the successive approximation image reconstruction, and the accuracy of the pixel value of the truncation portion can be increased accordingly.
また、本発明に係る放射線断層画像生成装置において、前記画像再構成部は、前記広く設定された円形領域を含む再構成画像領域の初期画像に、天板内部に対応する画素値を0とする先見情報を与えることが好ましい。天板内部には被検体が配置されないので、初期画像の天板内部に対応する画素値を0として画像更新を行わないようにしている。これにより、トランケーション部分の画素値の精度を高めることができる。 In the radiation tomographic image generation device according to the present invention, the image reconstruction unit sets a pixel value corresponding to the inside of the top board to 0 in the initial image of the reconstructed image region including the widely set circular region. It is preferable to provide foresight information. Since no subject is arranged inside the top plate, the pixel value corresponding to the inside of the top plate of the initial image is set to 0 so that the image is not updated. Thereby, the precision of the pixel value of the truncation part can be improved.
また、本発明に係る放射線断層画像生成装置において、前記画像再構成部は、前記広く設定された円形領域を含む再構成画像領域の初期画像に、予め設定された天板位置よりも低い画素値を0とする先見情報を与えることが好ましい。天板より低い位置には被検体が配置されないので、初期画像の予め設定された天板位置よりも低い画素値を0として画像更新を行わないようにしている。これにより、トランケーション部分の画素値の精度を高めることができる。 Further, in the radiation tomographic image generation device according to the present invention, the image reconstruction unit includes a pixel value lower than a preset top position in an initial image of the reconstructed image region including the widely set circular region. It is preferable to provide foresight information in which is set to 0. Since the subject is not arranged at a position lower than the top plate, the pixel value lower than the preset top plate position of the initial image is set to 0 so that the image is not updated. Thereby, the precision of the pixel value of the truncation part can be improved.
また、本発明に係る放射線断層画像生成装置において、前記実測投影データは、実測トランスミッション投影データと実測エミッション投影データであり、前記画像再構成部は、実測トランスミッション投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記実測トランスミッション投影データのみを用いて逐次近似画像再構成を行ってトランスミッション画像を生成するトランスミッション画像生成部であり、前記トランスミッション画像に基づいて求められた吸収補正用投影データを用いて実測エミッション投影データに対して吸収補正を行う吸収補正部と、前記吸収補正部で吸収補正された前記実測エミッション投影データに基づいて画像再構成を行ってエミッション画像を生成するエミッション画像生成部と、をさらに備えていることが好ましい。これにより、トランケーション部分の画像が推定されたトランスミッション画像に基づいて、吸収補正部は、実測エミッション投影データを吸収補正している。そのため、トランスミッション画像のトランケーション部分に起因するエミッション画像のアーチファクトを低減させることができる。 In the radiation tomographic image generation device according to the present invention, the measured projection data is measured transmission projection data and measured emission projection data, and the image reconstruction unit is a circle whose diameter is the data width of the measured transmission projection data. A transmission image generation unit that generates a transmission image by performing successive approximation image reconstruction using only the actually measured transmission projection data in a reconstructed image region including a circular region set wider than the region; An absorption correction unit that performs absorption correction on the measured emission projection data using the absorption correction projection data obtained based on the image, and image reconstruction based on the measured emission projection data that has been absorption corrected by the absorption correction unit Go to generate emission images Preferably further includes a transmission image generating unit. Thereby, based on the transmission image in which the image of the truncation portion is estimated, the absorption correction unit absorbs and corrects the actually measured emission projection data. Therefore, it is possible to reduce the artifact of the emission image due to the truncation part of the transmission image.
また、本発明に係る放射線断層画像生成装置において、前記エミッション画像生成部は、前記吸収補正部で吸収補正された実測エミッション投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記吸収補正部で吸収補正された実測エミッション投影データのみを用いて逐次近似画像再構成を行ってエミッション画像を生成することが好ましい。これにより、吸収補正部で吸収補正された実測エミッション投影データが、被検体Mが検出範囲からはみ出してトランケーションしている場合でも揃っている実測エミッション投影データのみでトランケーション部分の画像を推定することができる。 Further, in the radiation tomographic image generation device according to the present invention, the emission image generation unit has a circular area set wider than a circular area whose diameter is the data width of the measured emission projection data subjected to the absorption correction by the absorption correction unit. In the reconstructed image area including the image, it is preferable to perform successive approximate image reconstruction using only the actually measured emission projection data that has been subjected to the absorption correction by the absorption correction unit to generate an emission image. As a result, the image of the truncation portion can be estimated only by the actually measured emission projection data that is obtained even when the subject M is out of the detection range and is truncated when the actual emission projection data subjected to the absorption correction by the absorption correction unit. it can.
また、本発明に係る放射線断層画像生成装置において、前記実測投影データは、実測トランスミッション投影データと実測エミッション投影データであり、前記画像再構成部は、実測トランスミッション投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記実測トランスミッション投影データのみを用いて逐次近似画像再構成を行ってトランスミッション画像を生成するトランスミッション画像生成部であり、前記実測エミッション投影データに基づいて前記トランスミッション画像を利用した逐次近似画像再構成を行ってエミッション画像を生成するエミッション画像生成部と、をさらに備えていることが好ましい。これにより、トランケーション部分の画像が推定されたトランスミッション画像を利用した逐次近似画像再構成を行うことにより、エミッション画像生成部は、吸収補正とエミッション画像の生成を行っている。そのため、トランスミッション画像のトランケーション部分に起因するエミッション画像のアーチファクトを低減させることができる。 In the radiation tomographic image generation device according to the present invention, the measured projection data is measured transmission projection data and measured emission projection data, and the image reconstruction unit is a circle whose diameter is the data width of the measured transmission projection data. A transmission image generation unit configured to generate a transmission image by performing successive approximation image reconstruction using only the measured transmission projection data in a reconstructed image region including a circular region set wider than the region, and the measured emission projection It is preferable to further include an emission image generation unit that generates an emission image by performing successive approximation image reconstruction using the transmission image based on data. Accordingly, the emission image generation unit performs absorption correction and emission image generation by performing successive approximation image reconstruction using the transmission image in which the image of the truncation portion is estimated. Therefore, it is possible to reduce the artifact of the emission image due to the truncation part of the transmission image.
また、本発明に係る放射線断層画像生成装置において、前記エミッション画像生成部は、実測エミッション投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記実測エミッション投影データのみを用いて前記トランスミッション画像を利用した逐次近似画像再構成を行ってエミッション画像を生成することが好ましい。これにより、実測エミッション投影データが、被検体Mが検出範囲からはみ出してトランケーションしている場合でも揃っている実測エミッション投影データのみでトランケーション部分の画像を推定することができる。また、吸収補正も行うことができる。 Further, in the radiation tomographic image generation device according to the present invention, the emission image generation unit is a reconstructed image area including a circular area set wider than a circular area having a data width of measured emission projection data as a diameter. It is preferable to generate an emission image by performing successive approximation image reconstruction using the transmission image using only measured emission projection data. As a result, even when the measured emission projection data is truncated when the subject M is out of the detection range, the image of the truncation portion can be estimated only by the measured emission projection data. Absorption correction can also be performed.
また、本発明に係る放射線断層画像生成装置の一例において、前記放射線断層画像は、エミッション画像またはトランスミッション画像であることである。エミッション画像またはトランスミッション画像の放射線断層画像を生成する際に、実測投影データがトランケーションしている場合でも揃っている実測投影データのみでトランケーション部分の画像を推定することができる。 In the example of the radiation tomographic image generation apparatus according to the present invention, the radiation tomographic image is an emission image or a transmission image. When generating a radiation tomographic image of an emission image or a transmission image, an image of a truncation portion can be estimated using only the measured projection data that is available even when the measured projection data is truncated.
また、本発明に係る放射線断層画像生成装置の一例において、前記実測投影データは、複数の放射線検出器により収集されたものであることである。例えば、検出範囲の大きさが異なる場合でも、揃っている実測投影データのみで画像を推定することができる。 In the example of the radiation tomographic image generation apparatus according to the present invention, the actual projection data is collected by a plurality of radiation detectors. For example, even when the sizes of the detection ranges are different, it is possible to estimate an image using only the measured projection data that is available.
また、本発明に係る放射線断層撮影装置は、放射線断層画像を生成する放射線断層撮影装置であって、被検体に対して異なる複数の方向から実測投影データを収集するデータ収集部と、前記実測投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記実測投影データのみを用いて逐次近似画像再構成を行う画像再構成部と、を備えていることを特徴とするものである。 The radiation tomography apparatus according to the present invention is a radiation tomography apparatus that generates a radiation tomographic image, the data collection unit collecting actual projection data from a plurality of different directions with respect to a subject, and the actual projection An image reconstruction unit that performs successive approximation image reconstruction using only the measured projection data in a reconstructed image region including a circular region set wider than a circular region having a data data diameter as a diameter. It is characterized by being.
本発明に係る放射線断層撮影装置によれば、データ収集部は、被検体に対して異なる複数の方向から実測投影データを収集する。画像再構成部は、実測投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、実測投影データのみを用いて逐次近似画像再構成を行っている。逐次近似画像再構成の再構成画像領域を広く設定することにより、実測投影データが、被検体Mが検出範囲からはみ出してトランケーションしている場合でも揃っている実測投影データのみでトランケーション部分の画像を推定することができる。揃っている実測投影データは、放射線断層画像を生成するための投影データに推定された部分を有することによる誤差が生じないので、投影データ間の整合性が高くアーチファクトを低減させることができる。したがって、パラメータの最適化などの設定をほとんど必要としなく、トランケーション部分の画像推定が行われた放射線断層画像を容易に取得することができる。 According to the radiation tomography apparatus according to the present invention, the data collection unit collects measured projection data from a plurality of different directions with respect to the subject. The image reconstruction unit reconstructs successive approximate images using only the measured projection data in the reconstructed image region including the circular region set wider than the circular region whose diameter is the data width of the measured projection data. . By setting a wide range of reconstructed image areas for successive approximation image reconstruction, the image of the truncation portion can be obtained with only the measured projection data that is obtained even when the measured projection data is truncated even when the subject M is out of the detection range. Can be estimated. Since the measured projection data that has been prepared does not cause an error due to having an estimated portion in the projection data for generating the radiation tomographic image, the consistency between the projection data is high and artifacts can be reduced. Therefore, it is possible to easily obtain a radiation tomographic image in which the image estimation of the truncation portion is performed with little setting such as parameter optimization.
また、本発明に係る放射線断層画像生成プログラムは、被検体に対して異なる複数の方向から収集された実測投影データに基づいて放射線断層画像を生成する処理をコンピュータに実行させるための放射線断層画像生成プログラムであって、前記実測投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記実測投影データのみを用いて逐次近似画像再構成を行う工程を備えていることを特徴とするものである。 Further, the radiation tomographic image generation program according to the present invention is a radiation tomographic image generation for causing a computer to execute a process of generating a radiation tomographic image based on measured projection data collected from a plurality of different directions with respect to a subject. A step of performing successive approximation image reconstruction using only the measured projection data in a reconstructed image area including a circular area set wider than a circular area whose diameter is the data width of the measured projection data. It is characterized by having.
本発明に係る放射線断層画像生成プログラムによれば、実測投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、実測投影データのみを用いて逐次近似画像再構成を行っている。逐次近似画像再構成の再構成画像領域を広く設定することにより、実測投影データが、被検体Mが検出範囲からはみ出してトランケーションしている場合でも揃っている実測投影データのみでトランケーション部分の画像を推定することができる。揃っている実測投影データは、放射線断層画像を生成するための投影データに推定された部分を有することによる誤差が生じないので、投影データ間の整合性が高くアーチファクトを低減させることができる。したがって、パラメータの最適化などの設定をほとんど必要としなく、トランケーション部分の画像推定が行われた放射線断層画像を容易に取得することができる。 According to the radiation tomographic image generation program according to the present invention, a reconstructed image area including a circular area set wider than a circular area whose diameter is the data width of the measured projection data, and sequentially approximates using only the measured projection data. Image reconstruction is performed. By setting a wide range of reconstructed image areas for successive approximation image reconstruction, the image of the truncation portion can be obtained with only the measured projection data that is obtained even when the measured projection data is truncated even when the subject M is out of the detection range. Can be estimated. Since the measured projection data that has been prepared does not cause an error due to having an estimated portion in the projection data for generating the radiation tomographic image, the consistency between the projection data is high and artifacts can be reduced. Therefore, it is possible to easily obtain a radiation tomographic image in which the image estimation of the truncation portion is performed with little setting such as parameter optimization.
本発明に係る放射線断層画像生成装置、放射線断層撮影装置および放射線断層画像生成プログラムによれば、実測投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、実測投影データのみを用いて逐次近似画像再構成を行っている。逐次近似画像再構成の再構成画像領域を広く設定することにより、実測投影データが、被検体Mが検出範囲からはみ出してトランケーションしている場合でも揃っている実測投影データのみでトランケーション部分の画像を推定することができる。揃っている実測投影データは、放射線断層画像を生成するための投影データに推定された部分を有することによる誤差が生じないので、投影データ間の整合性が高くアーチファクトを低減させることができる。したがって、パラメータの最適化などの設定をほとんど必要としなく、トランケーション部分の画像推定が行われた放射線断層画像を容易に取得することができる。 According to the radiation tomographic image generation device, the radiation tomography apparatus, and the radiation tomographic image generation program according to the present invention, a reconstructed image region including a circular region that is set wider than a circular region whose diameter is the data width of measured projection data Thus, successive approximate image reconstruction is performed using only the actually measured projection data. By setting a wide range of reconstructed image areas for successive approximation image reconstruction, the image of the truncation portion can be obtained with only the measured projection data that is obtained even when the measured projection data is truncated even when the subject M is out of the detection range. Can be estimated. Since the measured projection data that has been prepared does not cause an error due to having an estimated portion in the projection data for generating the radiation tomographic image, the consistency between the projection data is high and artifacts can be reduced. Therefore, it is possible to easily obtain a radiation tomographic image in which the image estimation of the truncation portion is performed with little setting such as parameter optimization.
以下、図面を参照して本発明の実施例1を説明する。図1は、実施例1に係るPET/CT装置の構成を示すブロック図であり、図2は、被検体が撮影範囲からはみ出していることの説明に供する図である。図3は、実測投影データ(サイノグラム)の一例を示す図であり、図4は、逐次近似画像再構成の説明に供する図である。図5は、実測投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域の説明に供する図である。本実施例では、放射線断層撮影装置として、PET装置とX線CT装置とを組み合わせたPET/CT装置を一例に説明する。
図1を参照する。PET/CT装置1は、ベッド装置2を備えている。ベッド装置2は、被検体Mを載置する天板2aを上下に昇降移動させ、また、被検体Mの体軸に沿って平行移動させるように構成されている。PET/CT装置1は、天板2aに載置された被検体Mを診断するPET装置3およびX線CT装置4を備えている。
Please refer to FIG. The PET /
PET装置3は、開口部31aを有したガントリ31と、被検体Mから放射したγ線対を検出するγ線検出器32とを備えている。γ線検出器32は、開口部31a、すなわち、被検体Mの体軸周りを取り囲むようにしてリング状に構成され、ガントリ31内に収容されている。
The
γ線検出器32は、例えば、シンチレータブロックとライトガイドと光電子増倍管(PMT)とを備えている(いずれも図示しない)。シンチレータブロックは、複数個のシンチレータからなり、放射性薬剤が投与された被検体Mから放射したγ線を光に変換する。変換された光は、ライトガイドにより案内され、光電子増倍管により光電変換されて電気信号として出力される。
The γ-
PETデータ収集部36は、図示しない同時計数回路等を備えている。PETデータ収集部36は、γ線検出器32で検出されたγ線に基づいてPETデータ(エミッションデータ)、すなわち実際に測定された実測PET投影データを収集する。実測PET投影データは、被検体Mに対して異なる複数の方向から収集される。また、実測PET投影データは、収集した方向の順番に並べられてサイノグラムの形式にされる。
The PET
一方、X線CT装置4は、開口部41aを有したガントリ41を備えている。ガントリ41内には、被検体Mに向けてX線を照射するX線管42と、被検体Mを透過したX線を検出するX線検出器43とが設けられている。X線管42とX線検出器43は対向配置されており、開口部41a内に配置された被検体Mの体軸を中心に、X線管42とX線検出器43とが一体となって回転するように構成されている。X線管42は、ファン状のX線(ファンビーム)を照射する。X線検出器43は、円弧状(図2)の検出面を有して構成されるが、それ以外、例えば、フラットパネル型X線検出器(FPD)であってもよい。
On the other hand, the
X線管42は、X線管制御部44によりX線照射に必要な制御が実行される。X線管制御部44は、X線管42の管電圧や管電流を発生させる高電圧発生部45を有している。X線管制御部44は、管電圧や管電流や照射時間等のX線照射条件に応じてX線管42からX線を照射する。
The
X線CTデータ収集部46は、X線検出器43で検出されたX線に基づいてX線CTデータ(トランスミッションデータ)、すなわち、実際に測定された実測CT投影データを収集する。実測CT投影データは、図2に示すように、被検体Mに対して異なる複数の方向から収集される。また、実測PET投影データは、図3に示すように、収集した方向の順番に並べられてサイノグラム70の形式にされる。なお、ある方向の実測CT投影データを符号70aに示す。また、サイノグラム70の縦軸は、360°に限らず、例えば180°であってもよい。実測CT投影データは、検出されたX線強度と被検体Mが無いときのX線強度の比の対数をとって、吸収係数μ(cm−1)の投影データに処理される。なお、ファンビームの実測CT投影データは、パラレルビームの実測CT投影データに変換される。
The X-ray CT
本実施例では、図2に示すように、被検体Mは、X線検出器43の有効視野、すなわち検出範囲からはみ出して配置されているものとする。そのため、X線CTデータ収集部46には、被検体Mが検出範囲からはみ出して被検体Mが欠けているトランケーションした実測CT投影データが収集される(図3)。
In this embodiment, as shown in FIG. 2, it is assumed that the subject M is disposed so as to protrude from the effective visual field of the
PETデータ収集部36で収集された実測PET投影データ、およびX線CTデータ収集部46で収集された実測CT投影データは、断層画像生成部5に転送される。図1に戻る。断層画像生成部5は、X線CT画像生成部51、吸収係数マップ変換部52、補正用投影データ生成部53、吸収補正部54、PET画像生成部55および重ね合わせ部56を備えている。なお、断層画像生成部5は、本発明の放射線断層画像生成装置に相当し、X線CT画像生成部51は、本発明の画像再構成部に相当する。実測PET投影データは、本発明の実測エミッション投影データおよび実測投影データに相当する。実測CT投影データは、本発明の実測トランスミッション投影データおよび実測投影データに相当する。
The actually measured PET projection data collected by the PET
X線CT画像生成部51は、実測CT投影データを逐次近似法による画像再構成(逐次近似画像再構成)を行ってX線CT画像を生成する。逐次近似法による画像再構成は、例えばML−EM(maximum likelihood - expectation maximization)法により行われる。ML−EM法による画像再構成の計算式を下記の式(1)に示す。但し、式(1)はエミッション画像再構成などの放射型画像再構成を示しており、X線CT画像再構成やトランスミッション画像再構成などの透過型画像再構成においては、実測投影データや検出確率に対数変換等の処理が加わる場合がある。
The X-ray CT
λj (k):k回目の再構成画像
Σj′=1 maij′λj′ (k):順投影
Σi=1 n〔(yiaij)/(Σj′=1 maij′λj′ (k))〕:逆投影
Σi=1 naij:検出確率aijの総和
yi:実測投影データ
i:実測投影データの画素(1〜nまでの通し番号が付され、nはX線検出器の一列のX線検出素子数×投影方向数で算出される)
j:再構成画像の画素(1〜mまでの通し番号が付される)
j: Reconstructed image pixels (serial numbers from 1 to m are given)
ML−EM法による画像再構成を上述の式(1)および図4を参照して説明する。すなわち、式(1)に基づき、次の手順で再構成画像λjの各画素jの更新が行われる。まず、再構成画像λj (k)の各画素jの初期値を仮定する。例えば、全ての画素jの初期値をλj (k)=1とする。次に、再構成画像λj (k)(λj′ (k))からの順投影により投影データを算出する。なお、この算出された投影データを算出投影データと呼ぶこととする。実測投影データyiと算出投影データとの比を計算する。この実測投影データyiと算出投影データとの比を逆投影する。逆投影された値を検出確率の総和Σi=1 naijで規格化し、再構成画像λj (k)を乗算して再構成画像λj (k+1)を計算する。なお、再構成画像λjの各画素jの更新を予め設定された回数を繰り返し行うときは、計算されたλj (k+1)をλj (k)とし、順投影により算出投影データを算出する工程に戻る。
Image reconstruction by the ML-EM method will be described with reference to the above-described equation (1) and FIG. That is, based on equation (1), updating of each pixel j in the reconstructed image lambda j is performed in the following procedure. First, an initial value of each pixel j of the reconstructed image λ j (k) is assumed. For example, the initial value of all the pixels j is set to λ j (k) = 1. Next, projection data is calculated by forward projection from the reconstructed image λ j (k) (λ j ′ (k) ). The calculated projection data is referred to as calculated projection data. The ratio between the measured projection data y i and the calculated projection data is calculated. The ratio between the actually measured projection data y i and the calculated projection data is back projected. Backprojected values normalized by the
X線CT画像生成部51の再構成画像領域は、実測CT投影データのデータ幅(データ画素数)を直径とする円形領域よりも広く設定された円形領域を含んで構成される。すなわち、図5に示すように、実測CT投影データ71のデータ幅72を直径とする円形領域(有効視野)を符号73aで示し、円形領域73aを含む正方形の再構成画像領域を符号73で示す。一方、円形領域73aよりも広く設定された円形領域(有効視野)を符号74aで示し、円形領域74aを含む正方形の再構成画像領域を符号74で示す。なお、再構成画像領域74は、正方形に限定されず、例えば矩形であってもよい。
The reconstructed image area of the X-ray CT
広く設定された円形領域74aは、被検体Mのはみ出し部が無いように設定され、例えばPET装置3の有効視野と同じ大きさ(例えば直径60cm程度)に設定してもよい。再構成画像領域73,74の中心位置75は、ほぼ同じ位置である。また、各方向の実測CT投影データ71のデータ幅72の中心線76は、中心位置75を通るようになっている。
The
X線CT画像生成部51は、実測CT投影データ71のデータ幅72を直径とする円形領域73aよりも広く設定された円形領域74aを含む再構成画像領域74で、実測CT投影データ71のみを用いてML−EM法による画像再構成を行ってX線CT画像を生成する。実測CT投影データ71のみとは、収集されて揃っている実測CT投影データ71のことであり、トランケーション部分を推定した投影データを含まないこととする。
The X-ray CT
図6(a)〜図6(d)は、実測投影データのみを用いたML−EM法による画像再構成の説明に供する図である。図6(a)において、実測CT投影データ71を、そのデータ幅72で投影方向に反映させる。すなわち、実測CT投影データ71を、円形領域73aと、図6(a)中の斜線で示した領域74bとに反映するようにする。具体的には、例えば、まず、実測CT投影データ71と同じ投影方向であって、実測CT投影データ71と同じデータ幅の算出投影データを再構成画像領域74から算出する。そして、上述の式(1)の説明のように、算出投影データと実測CT投影データとの比を求めて比較する。なお、領域74bは、実測CT投影データ71と同じ投影(逆投影)方向であって、実測CT投影データ71と同じデータ幅で限定される領域であり、円形領域73aの外側の領域である。また、算出投影データのデータ幅が広く設定された再構成画像領域74であるときは、実測CT投影データがある部分のみ反映させるようにする。
FIGS. 6A to 6D are diagrams for explaining image reconstruction by the ML-EM method using only measured projection data. In FIG. 6A, the actual
図6(b)および図6(c)は、図6(a)と異なる方向の一例を示す。この場合も同様に、実測CT投影データ71は、円形領域73aと領域74bとに反映される。
FIG. 6B and FIG. 6C show an example of a direction different from that in FIG. In this case as well, the measured
また、図6(d)は、図6(a)〜図6(c)における実測CT投影データが反映される領域を重ね合わせた様子を示した図である。図6(d)中の線77は、投影線である。図6(d)に示すように、広く設定された円形領域74a(を含む再構成画像領域74)のうち有効視野73aの外側の領域では、実測CT投影データ量が制限されている。しかしながら、揃っている実測CT投影データは、推定による誤差を含まないデータであるので、実測CT投影データ間の整合性が高くアーチファクトを低減させることができる。
FIG. 6D is a diagram showing a state in which the regions reflecting the measured CT projection data in FIGS. 6A to 6C are superimposed. A
図1に戻る。吸収係数マップ変換部52は、X線CT画像生成部51で生成されたX線CT画像を511eVの吸収係数μ(cm−1)の吸収係数マップ(吸収係数画像)に変換する。補正用投影データ生成部53は、変換された吸収係数マップを異なる複数の方向に投影して吸収補正用投影データを生成する。吸収補正用投影データは、吸収補正部54に送られる。一方、吸収補正部54には、PETデータ収集部36から吸収補正部54に実測PET投影データが送られる。吸収補正部54は、X線CT画像に基づいて求められた吸収補正用投影データを用いて実測PET投影データに対して吸収補正を行う。
Returning to FIG. The absorption coefficient
PET画像生成部55は、吸収補正部54で吸収補正された実測PET投影データに基づいて画像再構成を行ってPET画像を生成する。画像再構成は、ML−EM法等の逐次近似法でもよく、その他の既存の方法(例えばフィルタ補正逆投影法(FBP法)等の解析的再構成法)でもよい。生成されたPET画像は、重ね合わせ部56または後述する表示部62に表示され、また、後述する記憶部64に記憶される。
The PET
重ね合わせ部56は、PET画像生成部55で生成されたPET画像と、X線CT画像生成部51で生成されたX線CT画像とを重ね合わせて合成画像を生成する。なお、重ね合わせ部56は、後述する記憶部64からPET画像とX線CT画像を取得して合成画像を生成するようにしてもよい。
The superimposing
断層画像生成部5の後段には、主制御部61、表示部62、入力部63および記憶部64が設けられている。主制御部61は、PET/CT装置1の各構成を統括的に制御し、中央演算処理装置(CPU)などで構成される。主制御部61は、例えば、ベッド装置2や、X線管42、X線検出器43、X線管制御部44等の制御を行う。なお、ベッド装置2やX線管42、X線検出器43等は、図示しないラックやピニオンやモータ等の移動機構により移動される。
A
表示部62は、モニタ等で構成され、PET画像やX線CT画像、合成画像を表示する。入力部63は、キーボードやマウス等で構成される。記憶部64は、ROM(Read-only Memory)、RAM(Random-Access Memory)またはハードディスク等、取り外し可能なものを含む記憶媒体で構成される。記憶部64には、PET画像やX線CT画像、合成画像が記憶される。なお、実測PET投影データおよび実測X線CT投影データを記憶部64に記憶させて、記憶部64からX線CT画像生成部51または吸収補正部54に実測PET投影データおよび実測X線CT投影データを送信してもよい。
The
断層画像生成部5は、CPU等で構成される制御部と、表示部と、ROMやRAMやハードディスク等、取り外し可能なものを含む記憶媒体で構成される記憶部とを備えたワークステーションやパーソナルコンピュータで構成してもよい(いずれも図示しない)。記憶部には、X線CT画像生成部51、吸収係数マップ変換部52、補正用投影データ生成部53、吸収補正部54、PET画像生成部55および重ね合わせ部56の各機能(後述するステップS01〜S04,S11〜S14)をプログラムとして記憶する。制御部は、記憶部から必要に応じてそのプログラムを読み出して、実測PET投影データおよび実測CT投影データを処理するようにしてもよい。
The tomographic
また、断層画像生成部5は、CPU等で構成される制御部で構成され、X線CT画像生成部51、吸収係数マップ変換部52、補正用投影データ生成部53、吸収補正部54、PET画像生成部55および重ね合わせ部56の各機能(後述するステップS01〜S04,S11〜S14)を動作プログラムとして記憶部64に記憶させる。そして、断層画像生成部5は、記憶部64から必要に応じてその動作プログラムを読み出して、実測PET投影データおよび実測CT投影データを処理してもよい。この場合において、主制御部61が記憶部64から必要に応じてその動作プログラムを読み出して、実測PET投影データおよび実測CT投影データを処理してもよい。また、その動作プログラムは、LAN等のネットワークシステムでPET/CT装置1と接続されたワークステーションやパーソナルコンピュータ上で実測PET投影データおよび実測CT投影データを処理してもよい。
The tomographic
次に本実施例のPET/CT装置1の動作について説明する。図7は、実施例に係るPET/CT装置(放射線断層撮影装置)の動作を示すフローチャートである。
Next, the operation of the PET /
〔ステップS01〕実測CT投影データの収集(実測トランスミッション投影データの収集)
ガントリ41の開口部41aに被検体Mの関心部位を配置させる。このとき、被検体Mは、図2に示すように、有効視野、すなわちファンビームの検出範囲からはみ出しているものとする。X線管42とX線検出器43とを一体に被検体Mの体軸を中心に回転させる。X線管42は、被検体Mに向けてX線を照射し、被検体Mを透過したX線はX線検出器43で検出される。X線CTデータ収集部46は、X線検出器43で検出されたX線に基づいて、被検体Mが検出範囲からはみ出してトランケーションした実測CT投影データを収集する。実測投影データは、図3に示すように、サイノグラム70形式に順番に並べられる。また、実測CT投影データは、X線CT画像生成部51に送られる。
[Step S01] Collection of measured CT projection data (collection of measured transmission projection data)
The region of interest of the subject M is placed in the
〔ステップS02〕X線CT画像の生成(トランスミッション画像の生成)
X線CT画像生成部51は、実測CT投影データのデータ幅を直径とする円形領域73aよりも広く設定された円形領域74aを含む再構成画像領域74で、実測CT投影データのみを用いてML−EM法(逐次近似法)による画像再構成を行ってX線CT画像を生成する。
[Step S02] X-ray CT image generation (transmission image generation)
The X-ray CT
図6(a)〜図6(d)に示すように、実測CT投影データ71は、広く設定された円形領域74aを含む再構成画像領域74のうち、実測CT投影データ71がある部分のみで反映するように画像再構成される。そのため、揃っている実測CT投影データは、推測による誤差を含まないデータであるので、投影データ間の整合性が高くアーチファクトとなりにくい。X線CT画像生成部51で生成されたX線CT画像は、吸収係数マップ変換部52に送られ、また、重ね合わせ部56や表示部62、記憶部64に送られる。
As shown in FIGS. 6A to 6D, the actual
〔ステップS03〕吸収係数マップ変換
吸収係数マップ変換部52は、X線CT画像のCT値(吸収係数(cm−1))に基づいて511eVの吸収係数μ(cm−1)の吸収係数マップに変換する。吸収係数マップは、補正用投影データ生成部53に送られる。
[Step S03] Absorption Coefficient Map Conversion The absorption coefficient
〔ステップS04〕吸収補正用投影データの生成
吸収補正用投影データ生成部53は、吸収係数マップを投影して吸収補正用投影データを生成する。吸収補正用投影データは、吸収補正部54に送られる。
[Step S04] Generation of Absorption Correction Projection Data The absorption correction projection
〔ステップS11〕実測PET投影データの収集(実測エミッション投影データの収集)
被検体Mには、ポジトロン放射性同位元素(ラジオアイソトープ:RI)で標識した放射性薬剤が予め投与されている。ガントリ31の開口部31aに被検体Mの関心部位を配置させる。被検体Mから180°反対方向に放射する2つのγ線は、γ線検出器32で検出される。PETデータ収集部36は、γ線検出器32で同時計数されたγ線の投影データ、実測PET投影データを収集する。実測PET投影データは、サイノグラム形式に順番に並べられる。なお、被検体Mは、PET装置3の検出範囲からはみ出していないものとする。実測PET投影データは、吸収補正部54に送られる。
[Step S11] Collection of measured PET projection data (collection of measured emission projection data)
The subject M is previously administered with a radiopharmaceutical labeled with a positron radioisotope (radioisotope: RI). The region of interest of the subject M is placed in the
〔ステップS12〕吸収補正
吸収補正部54は、吸収係数マップに基づいて求められた吸収補正用投影データを用いて実測PET投影データに対して吸収補正を行う。ここで、CT装置4の有効視野はPET装置3よりも狭く、実測CT投影データはトランケーションしている。しかしながら、X線CT画像生成部51により、トランケーション部分の画像が推定されたX線CT画像に基づいて、吸収補正部54は、実測PET投影データを吸収補正している。そのため、X線CT画像のトランケーション部分に起因するPET画像のアーチファクトを低減させることができる。すなわち、吸収補正におけるトランケーションアーチファクトの伝播を低減させることができる。
[Step S12] Absorption Correction The
〔ステップS13〕PET画像の生成(エミッション画像の生成)
PET画像生成部55は、吸収補正部54で吸収補正された実測PET投影データに基づいて画像再構成を行ってPET画像を生成する。生成されたPET画像は、重ね合わせ部56や表示部62、記憶部64に送られる。
[Step S13] Generation of PET image (generation of emission image)
The PET
〔ステップS14〕表示・保存
X線CT画像生成部51で生成されたX線CT画像、PET画像生成部55で生成されたPET画像、または重ね合わせ部56によりX線CT画像とPET画像とを重ね合わせた合成画像は、表示部62に表示される。また、X線CT画像、PET画像および合成画像は、記憶部64に記憶される。
[Step S14] Display / Save The X-ray CT image generated by the X-ray CT
本実施例に係るPET/CT装置1によれば、X線CTデータ収集部46は、被検体Mに対して異なる複数の方向から実測CT投影データを収集する。X線CT画像生成部51は、実測CT投影データのデータ幅を直径とする円形領域73aよりも広く設定された円形領域74aを含む再構成画像領域74で、実測CT投影データのみを用いてML−EM法による画像再構成を行っている。ML−EM法による画像再構成の再構成画像領域を広く設定することにより、実測CT投影データが、被検体Mが検出範囲からはみ出してトランケーションしている場合でも揃っている実測CT投影データのみでトランケーション部分の画像を推定することができる。揃っている実測CT投影データは、X線CT画像を生成するための投影データに推定された投影データ部分を有することによる誤差が生じないので、投影データ間の整合性が高くアーチファクトを低減させることができる。したがって、パラメータの最適化などの設定をほとんど必要としなく、トランケーション部分の画像推定が行われた放射線断層画像を容易に取得することができる。トランケーション補正付きの画像再構成をすることができる。
According to the PET /
なお、アーチファクト低減が実現した理由は、少ない投影データであっても投影データ間の整合が取れているので、逐次近似画像再構成の特性でもある、投影データに多少の欠損があっても最も確からしいところに収束する特性が生かされているためと考えられる。 The reason for the reduction in artifacts is that even if there is a small amount of projection data, the projection data is consistent, so this is also the characteristic of successive approximation image reconstruction. This is thought to be due to the fact that the characteristic of convergence at the apparent location is utilized.
また、X線CT画像に基づいて求められた吸収補正用投影データを用いて実測PET投影データに対して吸収補正を行う吸収補正部54と、吸収補正部54で吸収補正された実測PET投影データに基づいて画像再構成を行ってPET画像を生成するPET画像生成部55と、をさらに備えている。これにより、トランケーション部分の画像が推定されたX線CT画像に基づいて、吸収補正部54は、実測PET投影データを吸収補正している。そのため、X線CT画像のトランケーション部分に起因するPET画像のアーチファクトを低減させることができる。
Also, an
次に、図面を参照して本発明の実施例2を説明する。図8(a)〜図8(d)は実施例2に係る逐次近似画像再構成の初期画像の説明に供する図である。なお、図8(a)〜図8(d)において、図示の便宜上、初期画像の画像数を少なく図示している。また、実施例1と重複する説明は省略する。
Next,
実施例2では、実施例1の構成に加え、X線CT画像生成部51は、再構成画像(X線CT画像)領域の初期画像に、予め設定された先見(特徴)情報を与えるようになっている。図8(a)は先見情報を与えていない状態を示す。図8(a)中の円形領域(有効視野)74aには、初期画像の初期値として例えば“1”が与えられる。一方、再構成画像領域74であって円形領域74aの外側には、初期画像の初期値として“0(ゼロ)”が与えられる。これにより、円形領域74aの外側領域では画素値が更新されないようになっている。なお、必要に応じて円形領域74aの外側領域の初期画像の初期値を“1”としてもよい。
In the second embodiment, in addition to the configuration of the first embodiment, the X-ray CT
X線CT画像生成部51は、例えば、ML−EM法による画像再構成の際に、広く設定された円形領域74aを含む再構成画像領域74の初期画像に、天板2aに対応する画素を吸収係数μ(cm−1)とする先見情報を与える。図8(b)に示すように、ML−EM法による画像再構成の初期画像には、天板2aに対応する画素を吸収係数μとする先見情報が与えられる。これにより、ML−EM法による画像再構成の際に、天板2aに対応する画素値の精度を高めることができ、また、これに伴いトランケーション部分の画素値の精度を高めることができる。天板2aに対応する画素は、予め撮影されたX線CT画像の生データであるCT値(μ値)から設定される。このX線CT画像における天板2aの位置は、予め撮影されたX線CT画像の天板2aの位置から設定されている。天板2aの位置を被検体Mの体軸と垂直な面(XY面)で移動すると、その移動に合わせて初期画像内の天板2aに対応する画素位置が移動するようになっている。
The X-ray CT
また、X線CT画像生成部51は、ML−EM法による画像再構成の際に、広く設定された円形領域74aを含む再構成画像領域74の初期画像に、天板2a内部(例えば内部が中空や他材料で構成される)に対応する画素値を“0”とする先見情報を与えるようにしてもよい。すなわち、図8(c)に示すように、天板2a内部に対応する画素値に“0”が与えられる。この場合も同様に、天板2a内部に対応する画素は、予め撮影されたX線CT画像から設定されている。また、天板2aの移動に合わせて初期画像内の天板2a画像も移動する。天板2a内部には被検体Mが配置されないので、初期画像の天板2a内部に対応する画素値を“0”として画像更新を行わないようにしている。これにより、トランケーション部分の画素値の精度を高めることができる。なお、天板2a内部に対応する画素値を吸収係数μとしてもよい。
In addition, the X-ray CT
また、X線CT画像生成部51は、ML−EM法による画像再構成の際に、広く設定された円形領域74aを含む再構成画像74の初期画像に、予め設定された天板2a位置よりも低い位置の画素値を“0”とする先見情報を与えるようにしてもよい。すなわち、図8(d)に示すように、予め設定された天板2a位置によりも低い位置の画素値に“0”が与えられる。予め設定された天板2a位置は、例えば天板2aの載置面2dでもよいし、その反対側の面2eでもよい。また、載置面2dよりも下方の予め設定された位置2fでもよい。天板2aより低い位置には被検体Mが配置されないので、初期画像の予め設定された天板2a位置よりも低い画素値を“0”として画像更新を行わないようにしている。これにより、トランケーション部分の画素値の精度を高めることができる。なお、図8(b)〜図8(d)を組み合わせてもよい。
In addition, the X-ray CT
次に、図面を参照して本発明の実施例3を説明する。図9は、実施例3に係るPET/CT装置の構成を示すブロック図である。なお、各実施例と重複する説明は省略する。
Next,
上述の各実施例では、X線CTデータ収集部46で収集された実測CT投影データに基づいたX線CT画像から吸収補正用投影データを生成し、PETデータ収集部36で収集された実測PET投影データに対して吸収補正を行っていた。しかしながら、図9に示すように、PET/CT装置90は、外部線源91とγ線検出器92とを備え、トランスミッションデータ収集部96で収集された実測トランスミッション投影データに基づいて吸収補正用投影データを生成し、実測PET投影データに対して吸収補正を行うようにしてもよい。
In each of the embodiments described above, absorption correction projection data is generated from an X-ray CT image based on the actual CT projection data collected by the X-ray CT
すなわち、PET/CT装置70は、外部線源91、γ線検出器92、トランスミッションデータ収集部96、およびトランスミッション画像生成部97および吸収係数マップ変換部98をさらに備えている。外部線源91およびγ線検出器92は、例えばPET装置3に設けられる。
That is, the PET /
外部線源91は、例えば68Ge−68Ga、137Csでかつ線状または点状で構成される。外部線源91は、図示しない回転機構により、被検体Mの体軸を中心に被検体Mの周りを回転する。γ線検出器92は、γ線検出器32と同様に、被検体Mの体軸周りを取り囲むようにしてリング状に構成され、ガントリ31内に収容されている。γ線検出器92は、例えばシンチレータブロックとライトガイドと光電子増倍管とを備えている。外部線源91が68Ge−68Gaの場合は、例えば被検体Mの周りを回転する外部線源91の位置に合わせて同時計数を行うγ線検出器32の検出範囲が限定される。また、外部線源91が137Csの場合は、例えばファン状にγ線が照射される。そのため、図2と同様に、PET画像の有効視野に比べ、トランスミッション画像の有効視野が狭くなる。
The
γ線検出器92は、外部線源91から放射されて被検体Mを透過する等したγ線を検出する。トランスミッションデータ収集部96は、γ線検出器92で検出されたγ線に基づいて、トランスミッションデータ、すなわち実測トランスミッション投影データを収集する。なお、トランスミッションデータ収集部96は、外部線源91が68Ge−68Gaの場合、図示しない同時計数回路等を備えている。また、実測トランスミッション投影データは、収集した各方向の順番に並べられてサイノグラム70形式にされる。実測トランスミッション投影データは、例えば予め収集された被検体Mが無いときの実測ブランク投影データ(計数率)を実測トランスミッション投影データ(計数率)で割った値の対数をとって、吸収係数μ(cm−1)の投影データに処理される。
The γ-
トランスミッション画像生成部97は、実測トランスミッション投影データのデータ幅を直径とする円形領域73aよりも広く設定された円形領域74aを含む再構成画像領域74で、実測トランスミッション投影データのみを用いてML−EM法(逐次近似法)による画像再構成を行ってトランスミッション画像を生成する。なお、トランスミッション画像生成部97は、本発明の画像再構成部に相当する。実測トランスミッション投影データは、本発明の実測投影データに相当する。
The transmission
吸収係数マップ変換部98は、トランスミッション画像生成部97で生成されたトランスミッション画像を必要に応じて被検体Mから放射されるγ線のエネルギ(例えば511keV)に変換する。
The absorption coefficient
次に、本実施例のPET/CT装置90の動作について実施例1の図7を参照して説明する。なお、実施例1の図7のフローチャートにおいてステップS04,S11〜S14は、X線CT装置4での説明の一部が外部線源91とγ線検出器92等での説明に置き換わるのみ異なるので説明を省略する。
Next, the operation of the PET /
〔ステップS01a〕トランスミッションデータの収集
ガントリ31の開口部31aに被検体Mの関心部位を配置させる。このとき、被検体Mは、図2と同様に、有効視野、すなわちファンビームの検出範囲からはみ出しているものとする。外部線源(例えば137Cs)91は、リング状のγ線検出器92の内壁を沿うように円の軌跡を描きながら回転移動する。外部線源91は、被検体Mに向けてγ線を照射し、被検体Mを透過したX線はγ線検出器92で検出される。トランスミッションデータ収集部96は、γ線検出器92で検出されたγ線に基づいて、実測トランスミッション投影データを収集する。実測トランスミッション投影データμは、図3に示すように、サイノグラム70形式に順番に並べられる。また、実測トランスミッション投影データは、トランスミッション画像生成部97に送られる。
[Step S01a] Transmission Data Collection The site of interest of the subject M is placed in the
〔ステップS02a〕トランスミッション画像の生成
トランスミッション画像生成部97は、実測トランスミッション投影データのデータ幅を直径とする円形領域73aよりも広く設定された円形領域74aを含む再構成画像領域74で、実測トランスミッション投影データのみを用いてML−EM法(逐次近似法)による画像再構成を行ってトランスミッション画像を生成する。トランスミッション画像は、吸収係数マップ変換部98に送られる。
[Step S02a] Transmission Image Generation The transmission
〔ステップS03a〕吸収係数マップ変換
吸収係数マップ変換部98は、トランスミッション画像生成部97で生成されたトランスミッション画像を必要に応じて被検体Mから放射されるγ線のエネルギ(例えば511eV)の吸収係数μの吸収係数マップに変換する。吸収係数マップは、補正用投影データ生成部53に送られる。
[Step S03a] Absorption Coefficient Map Conversion The absorption coefficient
なお、図9中のX線CT生成部51は、実施例1と同様に、実測CT投影データのデータ幅を直径とする円形領域73aよりも広く設定された円形領域74aを含む再構成画像領域74で、実測CT投影データのみを用いてML−EM法による画像再構成を行ってX線CT画像を生成する。これにより、実測CT投影データが、被検体Mが検出範囲からはみ出してトランケーションしている場合でも揃っている実測CT投影データのみでトランケーション部分の画像を推定することができる。また、重ね合わせ部56は、PET画像生成部55で生成されたPET画像と、X線CT画像生成部51で生成されたX線CT画像とを重ね合わせて合成画像を生成する。
Note that the X-ray
本実施例に係るPET/CT装置90によれば、トランスミッションデータ収集部96は、被検体Mに対して異なる複数の方向から実測トランスミッション投影データを収集する。トランスミッション画像生成部97は、実測トランスミッション投影データのデータ幅を直径とする円形領域73aよりも広く設定された円形領域74aを含む再構成画像領域74で、実測トランスミッション投影データのみを用いてML−EM法による画像再構成を行っている。ML−EM法による画像再構成の再構成画像領域を広く設定することにより、実測トランスミッション投影データが、被検体Mが検出範囲からはみ出してトランケーションしている場合でも揃っている実測トランスミッション投影データのみでトランケーション部分の画像を推定することができる。揃っている実測トランスミッション投影データは、PET画像を生成するための投影データに推定された部分を有することによる誤差が生じないので、投影データ間の整合性が高くアーチファクトを低減させることができる。パラメータの最適化などの設定をほとんど必要としなく、トランケーション部分の画像推定が行われたPET画像を容易に取得することができる。
According to the PET /
また、トランスミッション画像に基づいて求められた吸収補正用投影データを用いて実測PET投影データに対して吸収補正を行う吸収補正部54と、吸収補正部54で吸収補正された実測PET投影データに基づいて画像再構成を行ってPET画像を生成するPET画像生成部36と、をさらに備えている。これにより、トランケーション部分の画像が推定されたトランスミッション画像に基づいて、吸収補正部54は、実測PET投影データを吸収補正している。そのため、トランスミッション画像のトランケーション部分に起因するPET画像のアーチファクトを低減させることができる。
Also, based on the
次に、図面を参照して本発明の実施例4を説明する。図10は、実施例4に係るSPECT/CT装置の構成を示すブロック図であり、図11は、実施例4に係るSPECT装置の構成を示すブロック図である。なお、各実施例と重複する説明は省略する。
Next,
上述した各実施例では、PET装置3とX線CT装置4とを備えたPET/CT装置1,90であった。しかしながら、SPECT装置とCT装置とを備えたSPECT/CT装置120またはSPECT装置130であってもよい。SPECT/CT装置120は、図10に示すように、被検体Mから放射されたγ線を検出するγ線検出器121と、被検体Mにファン状のX線を照射するX線管122と、X線管122と対向配置されて被検体Mを透過したX線を検出するX線検出器123とを備えている。γ線検出器121の入射面側には、パラレルコリメータ121aが設けられている。γ線検出器121とX線管122とX線検出器123は、図示しないガントリに回転自在に取り付けられ、図示しない回転機構により回転される。なお、X線検出器123は、FPD等で構成されている。
In each of the above-described embodiments, the PET /
また、SPECTデータ収集部124は、γ線検出器121で検出されたγ線に基づいて実測SPECT投影データを収集する。一方、トランスミッションデータ収集部125は、X線検出器123で検出されたX線に基づいて実測CT投影データを収集する。
The SPECT
トランスミッション画像生成部126は、実測CT投影データのデータ幅を直径とする円形領域73aよりも広く設定された円形領域74aを含む再構成画像領域74で、実測CT投影データのみを用いてML−EM法による画像再構成を行ってX線CT画像を生成する。吸収係数マップ変換部127は、生成されたX線CT画像を必要に応じて被検体Mから放射されるγ線のエネルギの吸収係数μの吸収係数マップに変換する。なお、トランスミッション画像生成部126は、本発明の画像再構成部に相当する。実測SPECT投影データは、本発明の実測エミッション投影データおよび実測投影データに相当する。
The transmission
SPECT画像生成部128は、実測SPECT投影データと吸収係数マップに基づいて、ML−EM法による画像再構成を行ってSPECT画像を生成する。ML−EM法による画像再構成では、吸収係数マップを利用して吸収補正を行いながら画像再構成を行う。
The SPECT
また、図10に示すX線管122およびX線検出器123に代えて、図11に示すように、外部線源131およびγ線検出器132を備えて、吸収係数マップを生成してもよい。すなわち、SPECT装置130は、外部線源131と、γ線検出器132と、トランスミッションデータ収集部133とを備えている。外部線源131は、例えば99mTcが用いられる。外部線源131の照射口には、コリメータ131aが設けられており、ファンビームを照射するようになっている。γ線検出器132の入射面側には、ファンビームコリメータ132aが設けられている。なお、外部線源131は、被検体Mに向けてファン状のγ線を照射し、被検体Mを透過したγ線は、ファンビームコリメータ132aを経由してγ線検出器132に入射する。
Further, instead of the
トランスミッションデータ収集部133は、γ線検出器132で検出されたγ線に基づいて実測トランスミッション投影データを収集する。また、その他の構成は、図10に示すSPECT/CT装置120と同じであるので説明を省略する。
The transmission
トランスミッション画像生成部126(図10)は、実測トランスミッション投影データのデータ幅を直径とする円形領域73aよりも広く設定された円形領域74aを含む再構成画像領域74で、実測トランスミッション投影データのみを用いてML−EM法による画像再構成を行ってトランスミッション画像を生成する。その後、トランスミッション画像を吸収係数マップに変換し、実測SPECT投影データと吸収係数マップに基づいて、吸収補正と画像再構成を行ってSPECT画像が生成される。
The transmission image generation unit 126 (FIG. 10) uses only the measured transmission projection data in the
なお、SPECT/CT装置120およびSPECT装置130の動作は、図7のフローチャートにおいて、吸収補正用投影データの生成(ステップS04)および吸収補正(ステップS12)の工程を含まない。そして、エミッション画像の生成(ステップS13)において、SPECT画像生成部128は、実測SPECT投影データと吸収係数マップに基づいて、ML−EM法による画像再構成を行ってSPECT画像を生成する。
Note that the operations of the SPECT /
本実施例に係るSPECT/CT装置120およびSPECT装置130によれば、トランスミッションデータ収集部125,133は、被検体Mに対して異なる複数の方向からの実測トランスミッション投影データを収集する。トランスミッション画像生成部126は、実測トランスミッション投影データのデータ幅を直径とする円形領域73aよりも広く設定された円形領域74aを含む再構成画像領域74で、実測トランスミッション投影データのみを用いてML−EM法による画像再構成を行ってトランスミッション画像を生成している。ML−EM法による画像再構成の再構成画像領域を広く設定することにより、実測トランスミッション投影データが、被検体Mが検出範囲からはみ出してトランケーションしている場合でも揃っている実測トランスミッション投影データのみでトランケーション部分の画像を推定することができる。揃っている実測トランスミッション投影データは、トランスミッション画像を生成するための投影データに推定された部分を有することによる誤差が生じないので、投影データ間の整合性が高くアーチファクトを低減させることができる。パラメータの最適化などの設定をほとんど必要としなく、トランケーション部分の画像推定が行われたトランスミッション画像を容易に取得することができる。
According to the SPECT /
また、実測SPECT投影データに基づいてトランスミッション画像を利用したML−EM法による画像再構成を行ってSPECT画像を生成するSPECT画像生成部と、をさらに備えている。これにより、トランケーション部分の画像が推定されたトランスミッション画像を利用したML−EM法による画像再構成を行うことにより、SPECT画像生成部128は、吸収補正とエミッション画像の生成を行っている。そのため、トランスミッション画像のトランケーション部分に起因するSPECT画像のアーチファクトを低減させることができる。
Further, a SPECT image generation unit that generates an SPECT image by performing image reconstruction by the ML-EM method using the transmission image based on the measured SPECT projection data is further provided. Thereby, the SPECT
本発明は、上記実施形態に限られることはなく、下記のように変形実施することができる。 The present invention is not limited to the above embodiment, and can be modified as follows.
(1)上述した各実施例では、トランスミッション画像を生成するX線CT画像作成部51およびトランスミッション画像作成部97,126は、図2に示すファンビームを照射することによるトランケーションした実測投影データに対して画像再構成を行っていた。しかしながら、図12に示すように、エミッション画像を生成するPET画像生成部55およびSPECT画像生成部128は、被検体Mが検出範囲からはみ出すことによるトランケーションした実測投影データに対しても、トランスミッション画像を生成した同様の処理で画像再構成を行ってもよい。
(1) In each of the embodiments described above, the X-ray CT
すなわち、実施例1〜3において、PET画像生成部55は、実測PET投影データのデータ幅を直径とする円形領域73aよりも広く設定された円形領域74aを含む再構成画像領域74で、実測PET投影データのみを用いてML−EM法による画像再構成を行ってPET画像を生成する。これにより、実測投影データが、被検体が検出範囲からはみ出してトランケーションしている場合でも揃っている実測投影データのみでトランケーション部分の画像を推定することができる。なお、この場合、PET画像生成部55は、本発明の画像再構成部に相当する。
That is, in the first to third embodiments, the PET
また、実施例4において、SPECT画像生成部128は、実測SPECT投影データのデータ幅を直径とする円形領域73aよりも広く設定された円形領域74aを含む再構成画像領域74で、実測SPECT投影データのみを用いてトランスミッション画像を利用したML−EM法による画像再構成を行ってSPECT画像を生成する。これにより、実測SPECT投影データが、被検体Mが検出範囲からはみ出してトランケーションしている場合でも揃っている実測投影データのみでトランケーション部分の画像を推定することができる。また、吸収補正も行うことができる。なお、この場合、SPECT画像生成部128は、本発明の画像再構成部に相当する。
In the fourth embodiment, the SPECT
(2)上述した各実施例および各変形例、例えば実施例1では、X線CT装置4は、X線管42とX線検出器43を備え、X線CTデータ収集部46は、1つのX線検出器43で検出されたX線に基づいて実測X線CT投影データを収集していた。しかしながら、図13(a)に示すように、X線CT装置140は、2つ(複数)のX線管42a,42bと2つのX線検出器43a,43bを備え、X線CTデータ収集部46は、2つのX線検出器43a,43bで検出されたX線に基づいて実測X線CT投影データを収集してもよい。
(2) In each of the above-described embodiments and modifications, for example, in
例えば、X線検出器43a,43bは、検出範囲の大きさが異なっていてもよい。図13(b)は、サイノグラム141の一例を示す図であり、符号142は、検出領域の大きさが異なることにより実測CT投影データの無い部分を示す。この場合、X線CT画像生成部は、図13中の実測CT投影データのデータ幅が大きい方を基準として、そのデータ幅を直径とする円形領域73aよりも広く設定された円形領域74aを含む再構成画像領域74で、実測CT投影データのみを用いてML−EM法による画像再構成を行ってX線CT画像を生成する。これにより、符号142の実測CT投影データが無い部分を含めて、被検体Mが検出範囲からはみ出してトランケーションしている場合でも揃っている実測投影データのみでトランケーション部分の画像を推定することができる。なお、X線CT装置4に限定されず、実施例3の外部線源91とγ線検出器92に適用させてもよく、実施例4のSPECT/CT装置120およびSPECT装置130に適用させてもよい。なお、X線管とX線検出器(放射線検出器)は、3つ以上であってもよい。また、符号143は、実測CT投影データを有する部分を示す。
For example, the
(3)上述した各実施例および各変形例では、逐次近似法としてML−EM法を用いていたがこれに限定されない。例えば、OS−EM(ordered subsets - expectation maximization)法、RAMLA(row-action maximum likelihood algorithm)法、DRAMA(dynamic RAMLA)法を用いてもよい。 (3) In each embodiment and each modification described above, the ML-EM method is used as the successive approximation method, but the present invention is not limited to this. For example, an OS-EM (ordered subsets-expectation maximization) method, a RAMLA (row-action maximum likelihood algorithm) method, or a DRAMA (dynamic RAMLA) method may be used.
(4)上述した各実施例および各変形例では、放射線断層撮影装置は、PET/CT装置1,90およびSPECT/CT装置120の複数の撮影装置(モダリティ)で構成されていたが、PET装置、SPECT装置およびX線CT装置等の単体の撮影装置で構成されていてもよい。また、複数の撮影装置の配置順序等が実施例と異なっていてもよい。
(4) In each of the above-described embodiments and modifications, the radiation tomography apparatus is configured by a plurality of imaging apparatuses (modalities) of the PET /
(5)上述した各実施例および各変形例において、ベッド装置2は、被検体Mが載置された天板2aをPET装置3、X線CT装置4、SPECT/CT装置120およびSPECT130等に対して相対的に移動していた。しかしながら、被検体Mが載置された天板2aに対してPET装置3、X線CT装置4、SPECT/CT装置120およびSPECT130等を移動する構成であってもよい。
(5) In each embodiment and each modification described above, the
(6)上述した各実施例および各変形例において、PETデータ収集部36、X線CTデータ収集部46、トランスミッションデータ収集部96,125,133およびSPECTデータ収集部124は、本発明のデータ収集部に相当する。
(6) In each of the above-described embodiments and modifications, the PET
1,90 … PET/CT装置
2a … 天板
3 … PET装置
4 … X線CT装置
36 … PETデータ収集部
46 … X線CTデータ収集部
51 … X線CT画像生成部
53 … 補正用投影データ生成部
54 … 吸収補正部
55 … PET画像生成部
61 … 主制御部
64 … 記憶部
70 … サイノグラム
71 … 実測CT投影データ
72 … データ幅
73 … 再構成画像領域
73a… 円形領域(有効視野)
74 … 再構成画像領域
74a… 円形領域(有効視野)
96,125 … トランスミッションデータ収集部
97,126,133 … トランスミッション画像生成部
120… SPECT/CT装置
124… SPECTデータ収集部
128… SPECT画像生成部
130… SPECT装置
DESCRIPTION OF
74 ...
96, 125 ... Transmission
Claims (13)
前記実測投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記実測投影データのみを用いて逐次近似画像再構成を行う画像再構成部を備えていることを特徴とする放射線断層画像生成装置。 In a radiation tomographic image generation device that generates a radiation tomographic image based on measured projection data collected from a plurality of different directions with respect to a subject,
An image reconstruction unit that performs successive approximation image reconstruction using only the measured projection data in a reconstructed image area including a circular area set wider than a circular area having a diameter of the data width of the measured projection data. A radiation tomographic image generation apparatus characterized by comprising:
前記画像再構成部は、前記広く設定された円形領域を含む再構成画像領域の初期画像に、天板に対応する画素値を吸収係数とする先見情報を与えることを特徴とする放射線断層画像生成装置。 The radiation tomographic image generation apparatus according to claim 1,
The image reconstruction unit is configured to provide foresight information having a pixel value corresponding to a top plate as an absorption coefficient to an initial image of a reconstructed image region including the widely set circular region, and generating a radiation tomographic image apparatus.
前記画像再構成部は、前記広く設定された円形領域を含む再構成画像領域の初期画像に、天板内部に対応する画素値を0とする先見情報を与えることを特徴とする放射線断層画像生成装置。 The radiation tomographic image generating apparatus according to claim 1 or 2,
Radiation tomographic image generation characterized in that the image reconstruction unit gives foreseeing information with a pixel value corresponding to 0 in the top plate to an initial image of the reconstructed image region including the widely set circular region apparatus.
前記画像再構成部は、前記広く設定された円形領域を含む再構成画像領域の初期画像に、予め設定された天板位置よりも低い画素値を0とする先見情報を与えることを特徴とする放射線断層画像生成装置。 The radiation tomographic image generation device according to any one of claims 1 to 3,
The image reconstruction unit provides foreseeing information having a pixel value lower than a preset top position to 0 to an initial image of a reconstructed image region including the widely set circular region. Radiation tomographic image generator.
前記実測投影データは、実測トランスミッション投影データと実測エミッション投影データであり、
前記画像再構成部は、実測トランスミッション投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記実測トランスミッション投影データのみを用いて逐次近似画像再構成を行ってトランスミッション画像を生成するトランスミッション画像生成部であり、
前記トランスミッション画像に基づいて求められた吸収補正用投影データを用いて実測エミッション投影データに対して吸収補正を行う吸収補正部と、
前記吸収補正部で吸収補正された前記実測エミッション投影データに基づいて画像再構成を行ってエミッション画像を生成するエミッション画像生成部と、をさらに備えていることを特徴とする放射線断層画像生成装置。 The radiation tomographic image generation device according to any one of claims 1 to 4,
The measured projection data is measured transmission projection data and measured emission projection data,
The image reconstruction unit is a reconstructed image region including a circular region set wider than a circular region whose diameter is the data width of the measured transmission projection data, and sequentially approximates image reconstruction using only the measured transmission projection data. Is a transmission image generation unit that generates a transmission image by performing
An absorption correction unit that performs absorption correction on the measured emission projection data using the absorption correction projection data obtained based on the transmission image;
A radiation tomographic image generation apparatus, further comprising: an emission image generation unit configured to perform image reconstruction based on the actually measured emission projection data subjected to absorption correction by the absorption correction unit to generate an emission image.
前記エミッション画像生成部は、前記吸収補正部で吸収補正された実測エミッション投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記吸収補正部で吸収補正された実測エミッション投影データのみを用いて逐次近似画像再構成を行ってエミッション画像を生成することを特徴とする放射線断層画像生成装置。 The radiation tomographic image generation device according to claim 5,
The emission image generation unit is a reconstructed image region including a circular region set wider than a circular region having a diameter of the data width of the actually measured emission projection data subjected to absorption correction by the absorption correction unit. A radiation tomographic image generation apparatus characterized by generating an emission image by performing successive approximation image reconstruction using only measured emission projection data subjected to absorption correction.
前記実測投影データは、実測トランスミッション投影データと実測エミッション投影データであり、
前記画像再構成部は、実測トランスミッション投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記実測トランスミッション投影データのみを用いて逐次近似画像再構成を行ってトランスミッション画像を生成するトランスミッション画像生成部であり、
前記実測エミッション投影データに基づいて前記トランスミッション画像を利用した逐次近似画像再構成を行ってエミッション画像を生成するエミッション画像生成部と、をさらに備えていることを特徴とする放射線断層画像生成装置。 The radiation tomographic image generation device according to any one of claims 1 to 4,
The measured projection data is measured transmission projection data and measured emission projection data,
The image reconstruction unit is a reconstructed image region including a circular region set wider than a circular region whose diameter is the data width of the measured transmission projection data, and sequentially approximates image reconstruction using only the measured transmission projection data. Is a transmission image generation unit that generates a transmission image by performing
A radiation tomographic image generation apparatus, further comprising: an emission image generation unit configured to perform successive approximation image reconstruction using the transmission image based on the measured emission projection data to generate an emission image.
前記エミッション画像生成部は、実測エミッション投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記実測エミッション投影データのみを用いて前記トランスミッション画像を利用した逐次近似画像再構成を行ってエミッション画像を生成することを特徴とする放射線断層画像生成装置。 The radiation tomographic image generation device according to claim 7,
The emission image generation unit uses the transmission image by using only the measured emission projection data in a reconstructed image area including a circular area set wider than a circular area whose diameter is the data width of the measured emission projection data. A radiation tomographic image generation apparatus that performs the successive approximation image reconstruction and generates an emission image.
前記放射線断層画像は、エミッション画像であることを特徴とする放射線断層画像生成装置。 The radiation tomographic image generation device according to any one of claims 1 to 4,
The radiation tomographic image generation apparatus, wherein the radiation tomographic image is an emission image.
前記放射線断層画像は、トランスミッション画像であることを特徴とする放射線断層画像生成装置。 The radiation tomographic image generation device according to any one of claims 1 to 4,
The radiation tomographic image generating apparatus, wherein the radiation tomographic image is a transmission image.
前記実測投影データは、複数の放射線検出器により収集されたものであることを特徴とする放射線断層画像生成装置。 The radiation tomographic image generation device according to any one of claims 1 to 4,
3. The radiation tomographic image generation apparatus according to claim 1, wherein the actual projection data is collected by a plurality of radiation detectors.
被検体に対して異なる複数の方向から実測投影データを収集するデータ収集部と、
前記実測投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記実測投影データのみを用いて逐次近似画像再構成を行う画像再構成部と、
を備えていることを特徴とする放射線断層撮影装置。 A radiation tomography apparatus for generating a radiation tomographic image,
A data collection unit for collecting measured projection data from a plurality of different directions with respect to the subject;
An image reconstruction unit that performs successive approximation image reconstruction using only the measured projection data in a reconstructed image region that includes a circular region set wider than a circular region having a data width of the measured projection data as a diameter;
A radiation tomography apparatus comprising:
前記実測投影データのデータ幅を直径とする円形領域よりも広く設定された円形領域を含む再構成画像領域で、前記実測投影データのみを用いて逐次近似画像再構成を行う工程を備えていることを特徴とする放射線断層画像生成プログラム。 A radiation tomographic image generation program for causing a computer to execute a process of generating a radiation tomographic image based on measured projection data collected from a plurality of different directions with respect to a subject,
A step of performing successive approximation image reconstruction using only the measured projection data in a reconstructed image area including a circular area set wider than a circular area whose diameter is the data width of the measured projection data. A radiation tomographic image generation program characterized by
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012090264A JP5849838B2 (en) | 2012-04-11 | 2012-04-11 | Radiation tomographic image generation apparatus, radiation tomography apparatus, and radiation tomographic image generation program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012090264A JP5849838B2 (en) | 2012-04-11 | 2012-04-11 | Radiation tomographic image generation apparatus, radiation tomography apparatus, and radiation tomographic image generation program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2013215472A true JP2013215472A (en) | 2013-10-24 |
JP5849838B2 JP5849838B2 (en) | 2016-02-03 |
Family
ID=49588348
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012090264A Expired - Fee Related JP5849838B2 (en) | 2012-04-11 | 2012-04-11 | Radiation tomographic image generation apparatus, radiation tomography apparatus, and radiation tomographic image generation program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5849838B2 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2018042730A (en) * | 2016-09-14 | 2018-03-22 | 住友重機械工業株式会社 | X-ray CT apparatus |
JP2021032816A (en) * | 2019-08-28 | 2021-03-01 | ヤマハ発動機株式会社 | Inspection device |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004237076A (en) * | 2002-10-04 | 2004-08-26 | Ge Medical Systems Global Technology Co Llc | Method and apparatus for multimodality imaging |
WO2011122613A1 (en) * | 2010-03-30 | 2011-10-06 | 株式会社 日立メディコ | Reconstruction computing device, reconstruction computing method, and x-ray ct device |
JP2013116213A (en) * | 2011-12-02 | 2013-06-13 | Hitachi Medical Corp | X-ray ct apparatus |
-
2012
- 2012-04-11 JP JP2012090264A patent/JP5849838B2/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004237076A (en) * | 2002-10-04 | 2004-08-26 | Ge Medical Systems Global Technology Co Llc | Method and apparatus for multimodality imaging |
WO2011122613A1 (en) * | 2010-03-30 | 2011-10-06 | 株式会社 日立メディコ | Reconstruction computing device, reconstruction computing method, and x-ray ct device |
JP2013116213A (en) * | 2011-12-02 | 2013-06-13 | Hitachi Medical Corp | X-ray ct apparatus |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2018042730A (en) * | 2016-09-14 | 2018-03-22 | 住友重機械工業株式会社 | X-ray CT apparatus |
JP2021032816A (en) * | 2019-08-28 | 2021-03-01 | ヤマハ発動機株式会社 | Inspection device |
JP7303069B2 (en) | 2019-08-28 | 2023-07-04 | ヤマハ発動機株式会社 | inspection equipment |
Also Published As
Publication number | Publication date |
---|---|
JP5849838B2 (en) | 2016-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6747885B2 (en) | Radiation diagnostic apparatus and method for operating radiation diagnostic apparatus | |
EP2041606B1 (en) | Energy spectrum reconstruction | |
JP5860607B2 (en) | System and method for tomographic data collection and image reconstruction | |
JP6192542B2 (en) | Truncation correction for iterative cone beam CT reconstruction for SPECT / CT systems | |
EP3067864B1 (en) | Iterative reconstruction with enhanced noise control filtering | |
US7840052B2 (en) | Restoration of the nuclear medicine 2D planar image by iterative constrained deconvolution | |
JP5611640B2 (en) | Nuclear medicine diagnostic apparatus, medical image processing apparatus, and medical image diagnostic apparatus | |
JP2008062035A5 (en) | ||
JP2005312970A (en) | Reconstruction method of projection data set during dose reduced partial spiral scanning of reduced radiation dosage in computerized tomography | |
US10049465B2 (en) | Systems and methods for multi-modality imaging component alignment | |
US11419566B2 (en) | Systems and methods for improving image quality with three-dimensional scout | |
US20050004452A1 (en) | Systems and methods for analyzing an abnormality of an object | |
US9858688B2 (en) | Methods and systems for computed tomography motion compensation | |
JP2009047602A (en) | Positron emission computerd tomograph, attenuation map creating device, and attenuation map creating program | |
US7853314B2 (en) | Methods and apparatus for improving image quality | |
US10383589B2 (en) | Direct monochromatic image generation for spectral computed tomography | |
JP5849838B2 (en) | Radiation tomographic image generation apparatus, radiation tomography apparatus, and radiation tomographic image generation program | |
US7155047B2 (en) | Methods and apparatus for assessing image quality | |
JP2008267913A (en) | Nuclear medicine diagnostic apparatus and diagnostic system used for same | |
JP6629158B2 (en) | X-ray CT system | |
US10755450B2 (en) | Display of depth location of computed tomography slice images relative to an object to be imaged | |
US20200294281A1 (en) | X-ray ct apparatus, image reconstruction device, and image reconstruction method | |
US20230320685A1 (en) | Nuclear medicine diagnostic apparatus, nuclear medicine diagnostic method, and storage medium | |
JP4371636B2 (en) | Nuclear medicine diagnostic equipment | |
Kim et al. | Performance characterization of microtomography with complementary metal-oxide-semiconductor detectors for computer-aided defect inspection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20140714 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20150213 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150317 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150515 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20151104 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20151117 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 5849838 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
LAPS | Cancellation because of no payment of annual fees |