JPH04370784A - Ｓｐｅｃｔ画像の再構成方法 - Google Patents
Ｓｐｅｃｔ画像の再構成方法Info
- Publication number
- JPH04370784A JPH04370784A JP3147596A JP14759691A JPH04370784A JP H04370784 A JPH04370784 A JP H04370784A JP 3147596 A JP3147596 A JP 3147596A JP 14759691 A JP14759691 A JP 14759691A JP H04370784 A JPH04370784 A JP H04370784A
- Authority
- JP
- Japan
- Prior art keywords
- image
- collimator
- fan beam
- spect
- data
- 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
- 238000000034 method Methods 0.000 title claims description 13
- 238000002603 single-photon emission computed tomography Methods 0.000 claims abstract description 25
- 238000012545 processing Methods 0.000 claims abstract description 5
- 239000000941 radioactive substance Substances 0.000 claims abstract description 4
- 238000001514 detection method Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 abstract description 2
- 230000000007 visual effect Effects 0.000 abstract 3
- 230000004304 visual acuity Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 9
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000013480 data collection Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005251 gamma ray Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
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/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
- G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Nuclear Medicine (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【０００１】
【産業上の利用分野】本発明は、被写体内に分布した放
射性物質（以下、「ＲＩ」という）から放出されるγ線
を検出して、その３次元分布像を得るＳＰＥＣＴ画像の
再構成方法に関する。
射性物質（以下、「ＲＩ」という）から放出されるγ線
を検出して、その３次元分布像を得るＳＰＥＣＴ画像の
再構成方法に関する。
【０００２】
【従来の技術】従来、放射線を検出した画像再構成する
装置として、Ｘ線ＣＴ装置やＳＰＥＣＴ装置などがある
。そして、Ｘ線ＣＴ装置は、Ｘ線源がＸ線ファンビーム
の焦点位置となる幾何学的位置関係でＸ線源と検出素子
群とを設け、被写体の体軸を回転中心としてそのＸ線源
及び検出素子群を回転させながらＸ線源より直線的に被
写体を透過してきたＸ線を検出素子群で検出して投影デ
ータを収集し、この収集した投影データを基に被写体に
おけるＸ線の吸収分布を画像再構成するものである。
装置として、Ｘ線ＣＴ装置やＳＰＥＣＴ装置などがある
。そして、Ｘ線ＣＴ装置は、Ｘ線源がＸ線ファンビーム
の焦点位置となる幾何学的位置関係でＸ線源と検出素子
群とを設け、被写体の体軸を回転中心としてそのＸ線源
及び検出素子群を回転させながらＸ線源より直線的に被
写体を透過してきたＸ線を検出素子群で検出して投影デ
ータを収集し、この収集した投影データを基に被写体に
おけるＸ線の吸収分布を画像再構成するものである。
【０００３】他方、ＳＰＥＣＴ装置は、被写体内に分布
したＲＩから放出されるγ線を検出するため、仮想焦点
を構成するファンビームコリメータと、このファンビー
ムコリメータにより選択通過したγ線を光に変換して検
出するシンチレータ検出器とからなるガンマカメラを設
け、このガンマカメラを用いて収集した投影データを基
に被写体におけるＲＩそのものの分布をγ線の吸収を無
視して画像再構成するものである。
したＲＩから放出されるγ線を検出するため、仮想焦点
を構成するファンビームコリメータと、このファンビー
ムコリメータにより選択通過したγ線を光に変換して検
出するシンチレータ検出器とからなるガンマカメラを設
け、このガンマカメラを用いて収集した投影データを基
に被写体におけるＲＩそのものの分布をγ線の吸収を無
視して画像再構成するものである。
【０００４】いずれの装置においても、画像再構成する
場合に、フィルター補正逆投影法（ｆｉｌｔｅｒｄ ｂ
ａｃｋ ｐｒｏｊｅｃｔｉｏｎ）が適用されることが多
い。このフィルター補正逆投影法を適用した一例として
「資料１」に示すものがある。
場合に、フィルター補正逆投影法（ｆｉｌｔｅｒｄ ｂ
ａｃｋ ｐｒｏｊｅｃｔｉｏｎ）が適用されることが多
い。このフィルター補正逆投影法を適用した一例として
「資料１」に示すものがある。
【０００５】資料１：ＣＯＮＶＯＬＵＴＩＯＮ ＲＥＣ
ＯＮＳＴＲＵＣＴＩＯＮ ＴＥＣＨＮＩＱＵＥＳ ＦＯ
Ｒ ＤＩＶＥＲＧＥＮＴ ＢＥＡＭＳ．ｃｏｍｐｕｔ．
Ｂｉｏｌ，Ｍｅｄ，１９７６，ＶＯＬ．６，ＰＰ２５９
−２７１ この資料１には、被写体の周囲を回転するＸ線源から被
写体に向けてＸ線ファンビームを爆射し、この被写体を
透過したＸ線を多数の検出素子を直線状あるいは円弧上
に配列した検出器で検出し、この検出した投影データを
フィルタ処理により補正後、Ｘ線ビームに沿って有効視
野内全部に逆投影してＸ線ＣＴ画像を再構成する旨が示
されている。この場合には、Ｘ線源が焦点にあるため、
画像上のある点における３６０°のあらゆる方向からの
情報量は相対的に等価の関係で取扱える。
ＯＮＳＴＲＵＣＴＩＯＮ ＴＥＣＨＮＩＱＵＥＳ ＦＯ
Ｒ ＤＩＶＥＲＧＥＮＴ ＢＥＡＭＳ．ｃｏｍｐｕｔ．
Ｂｉｏｌ，Ｍｅｄ，１９７６，ＶＯＬ．６，ＰＰ２５９
−２７１ この資料１には、被写体の周囲を回転するＸ線源から被
写体に向けてＸ線ファンビームを爆射し、この被写体を
透過したＸ線を多数の検出素子を直線状あるいは円弧上
に配列した検出器で検出し、この検出した投影データを
フィルタ処理により補正後、Ｘ線ビームに沿って有効視
野内全部に逆投影してＸ線ＣＴ画像を再構成する旨が示
されている。この場合には、Ｘ線源が焦点にあるため、
画像上のある点における３６０°のあらゆる方向からの
情報量は相対的に等価の関係で取扱える。
【０００６】
【発明が解決しようとする課題】しかしながら、ファン
ビームコリメータを用いたＳＰＥＣＴ装置にフィルター
補正逆投影法を適用して画像再構成した場合においては
、コリメータに近い所の空間分解能は良いがそのコリメ
ータから離れるに従って空間分解能が低下してゆくため
に、コリメータから離れた所に於けるＲＩ分布は投影デ
ータの情報に充分反映されない。換言すれば、再構成さ
れた画像の解像力が著しく劣化することになる。
ビームコリメータを用いたＳＰＥＣＴ装置にフィルター
補正逆投影法を適用して画像再構成した場合においては
、コリメータに近い所の空間分解能は良いがそのコリメ
ータから離れるに従って空間分解能が低下してゆくため
に、コリメータから離れた所に於けるＲＩ分布は投影デ
ータの情報に充分反映されない。換言すれば、再構成さ
れた画像の解像力が著しく劣化することになる。
【０００７】従って、従来は、ＳＰＥＣＴ装置により画
像再構成したＳＰＥＣＴ画像を得た場合、そのＳＰＥＣ
Ｔ画像の解像力が不充分のものであった。
像再構成したＳＰＥＣＴ画像を得た場合、そのＳＰＥＣ
Ｔ画像の解像力が不充分のものであった。
【０００８】本発明は、上記事情に着目してなされたも
ので、その目的とするところは、ＳＰＥＣＴ画像の解像
度を向上させることができるＳＰＥＣＴ画像の再構成方
法を提供することにある。
ので、その目的とするところは、ＳＰＥＣＴ画像の解像
度を向上させることができるＳＰＥＣＴ画像の再構成方
法を提供することにある。
【０００９】
【課題を解決するための手段】本発明は上記の目的を達
成するために、被写体内に分布した放射性物質が放出さ
れるγ線をファンビーム状に異なる方向から検出し、こ
の検出により収集した投影データをフィルタ処理により
補正後、逆投影してＳＰＥＣＴ画像を再構成する際、逆
投影する範囲を、再構成画像上の個々の画素についてフ
ァンビームコリメータに近い側の有効視野内とすること
を特徴とする。
成するために、被写体内に分布した放射性物質が放出さ
れるγ線をファンビーム状に異なる方向から検出し、こ
の検出により収集した投影データをフィルタ処理により
補正後、逆投影してＳＰＥＣＴ画像を再構成する際、逆
投影する範囲を、再構成画像上の個々の画素についてフ
ァンビームコリメータに近い側の有効視野内とすること
を特徴とする。
【００１０】
【作用】本発明にあっては、有効視野内において、コリ
メータに近い側の投影データの方がコリメータから遠い
側のそれよりも多くの情報量をになっている点、及びコ
リメータに近いほど被写体によるγ線の吸収が少い点の
各観点に従って逆投影する範囲を、再構成画像上の個々
の画素についてコリメータに近い側の有効視野内とした
。このように逆投影する範囲を設定して投影データを逆
投影すると、ＳＰＥＣＴ画像の解像力が従来に比べて大
幅に向上したものとなる。
メータに近い側の投影データの方がコリメータから遠い
側のそれよりも多くの情報量をになっている点、及びコ
リメータに近いほど被写体によるγ線の吸収が少い点の
各観点に従って逆投影する範囲を、再構成画像上の個々
の画素についてコリメータに近い側の有効視野内とした
。このように逆投影する範囲を設定して投影データを逆
投影すると、ＳＰＥＣＴ画像の解像力が従来に比べて大
幅に向上したものとなる。
【００１１】
【実施例】図１は、本発明のＳＰＥＣＴ画像の再構成方
法が適用されたＳＰＥＣＴ装置の概略を示す図である。
法が適用されたＳＰＥＣＴ装置の概略を示す図である。
【００１２】このＳＰＥＣＴ装置は、架台１、データ収
集部２、画像再構成部３、表示部４からなり、画像再構
成部３で実施する画像再構成処理に本発明方法を以下説
明するように用いる。
集部２、画像再構成部３、表示部４からなり、画像再構
成部３で実施する画像再構成処理に本発明方法を以下説
明するように用いる。
【００１３】架台１は、被検体Ｐに分布したＲＩから放
出されるγ線を３６０°方向から検出するため、ファン
ビームコリメータ５が取付けられたガンマカメラ６を被
検体Ｐを中心として３６０°回転させることができる。 あるいは、２台以上のガンマカメラを用い、それぞれ分
割した角度だけ回転させ、全体として３６０°方向から
の投影データを収集できるようにしてもよい。この架台
１のガンマカメラ６からの投影データは、データ収集部
２で収集される。これにより画像再構成部３では、デー
タ収集部２より受けた投影データをフィルター処理によ
り補正後、逆投影してＳＰＥＣＴ画像を画像再構成する
際、逆投影する範囲を、再構成画像上の個々の画素につ
いてファンビームコリメータ５に近い側の有効視野内と
する。
出されるγ線を３６０°方向から検出するため、ファン
ビームコリメータ５が取付けられたガンマカメラ６を被
検体Ｐを中心として３６０°回転させることができる。 あるいは、２台以上のガンマカメラを用い、それぞれ分
割した角度だけ回転させ、全体として３６０°方向から
の投影データを収集できるようにしてもよい。この架台
１のガンマカメラ６からの投影データは、データ収集部
２で収集される。これにより画像再構成部３では、デー
タ収集部２より受けた投影データをフィルター処理によ
り補正後、逆投影してＳＰＥＣＴ画像を画像再構成する
際、逆投影する範囲を、再構成画像上の個々の画素につ
いてファンビームコリメータ５に近い側の有効視野内と
する。
【００１４】このように逆投影する範囲をファンビーム
コリメータ５に近い側の有効視野内とするのは以下の理
由による。
コリメータ５に近い側の有効視野内とするのは以下の理
由による。
【００１５】ファンビームコリメータを用いてＳＰＥＣ
Ｔ画像を再構成する場合には、図２で示すような有効視
野（ＦＯＶ：Ｆｉｅｌｄ ｏｆ Ｖｉｅｗ ）とファン
ビームコリメータ５及び焦点Ｓの幾何学的な位置関係を
採る。なお、図２中、βはガンマカメラの角度、Ｐ（β
，Ｓ）はｌ−ｌ′方向投影像である。
Ｔ画像を再構成する場合には、図２で示すような有効視
野（ＦＯＶ：Ｆｉｅｌｄ ｏｆ Ｖｉｅｗ ）とファン
ビームコリメータ５及び焦点Ｓの幾何学的な位置関係を
採る。なお、図２中、βはガンマカメラの角度、Ｐ（β
，Ｓ）はｌ−ｌ′方向投影像である。
【００１６】この図２の関係にあるとき、画像再構成が
された結果であるｆ（ｘ，ｙ）は、次の（１）式で表現
される。なお、（１） 式中、−πからπまでの積分は
逆投影の式である。
された結果であるｆ（ｘ，ｙ）は、次の（１）式で表現
される。なお、（１） 式中、−πからπまでの積分は
逆投影の式である。
【００１７】
【数１】
【００１８】このとき、ガンマカメラからの投影データ
：投影像Ｐ（β，Ｓ）は、再構成フィルターによりコン
ボリューションされて、−πからπまで３６０°方向の
データについて有効視野内全ての範囲ｌ−ｌ′に沿って
逆投影されることを意味している。
：投影像Ｐ（β，Ｓ）は、再構成フィルターによりコン
ボリューションされて、−πからπまで３６０°方向の
データについて有効視野内全ての範囲ｌ−ｌ′に沿って
逆投影されることを意味している。
【００１９】しかし、逆投影する範囲を、従来のように
有効視野内全てとした場合には、図３のように再構成画
像上の点Ｘに注目すると、この点Ｘを通る３６０°のあ
らゆる方向からデータ収集が行われたことになる。この
場合、点Ｘを通る任意の直線をひいた関係において、点
Ｘに近いコリメータ位置Ａのデータと点ＸのＲＩ濃度と
の相関、点Ｘに遠いコリメータ位置Ｂのデータと点Ｘの
ＲＩ濃度との相関と実際上相違することになる。
有効視野内全てとした場合には、図３のように再構成画
像上の点Ｘに注目すると、この点Ｘを通る３６０°のあ
らゆる方向からデータ収集が行われたことになる。この
場合、点Ｘを通る任意の直線をひいた関係において、点
Ｘに近いコリメータ位置Ａのデータと点ＸのＲＩ濃度と
の相関、点Ｘに遠いコリメータ位置Ｂのデータと点Ｘの
ＲＩ濃度との相関と実際上相違することになる。
【００２０】その第１の理由は、点ＸのＲＩ濃度とファ
ンビームコリメータにより選択通過したγ線を検出する
ガンマカメラ６でのカウントとの相関は、図４のように
ファンビームコリメータ５表面からの距離ｄが大きいほ
ど小さくなる。即ち、解像限界（ＦＷＨＭ）が距離ｄに
よらず一定となるような理想的なファンビームコリメー
タであれば、図３においてコリメータ位置Ａのデータと
点ＸのＲＩ濃度との相関は、コリメータ位置Ｂのデータ
と点ＸのＲＩ濃度との相関と同じになる。実際は、図４
のようにコリメータ位置Ｂよりコリメータ位置Ａのデー
タの方が点ＸのＲＩ濃度との相関関係が高い。つまり、
より多くの情報量とになっている。
ンビームコリメータにより選択通過したγ線を検出する
ガンマカメラ６でのカウントとの相関は、図４のように
ファンビームコリメータ５表面からの距離ｄが大きいほ
ど小さくなる。即ち、解像限界（ＦＷＨＭ）が距離ｄに
よらず一定となるような理想的なファンビームコリメー
タであれば、図３においてコリメータ位置Ａのデータと
点ＸのＲＩ濃度との相関は、コリメータ位置Ｂのデータ
と点ＸのＲＩ濃度との相関と同じになる。実際は、図４
のようにコリメータ位置Ｂよりコリメータ位置Ａのデー
タの方が点ＸのＲＩ濃度との相関関係が高い。つまり、
より多くの情報量とになっている。
【００２１】一方、コリメータ位置Ｂのデータからは点
線の周辺の点のＲＩ濃度の情報を点ＸでのＲＩ濃度の情
報と区別できない。これはコリメータ位置Ｂから点Ｘの
距離が長いのでＦＷＨＭがそれだけ広いからである。
線の周辺の点のＲＩ濃度の情報を点ＸでのＲＩ濃度の情
報と区別できない。これはコリメータ位置Ｂから点Ｘの
距離が長いのでＦＷＨＭがそれだけ広いからである。
【００２２】第２の理由は、コリメータ位置ＡＢが点Ｘ
の間の距離が、コリメータ位置Ｂ及び点Ｘの間の距離よ
り短いということは、点Ｘからコリメータ位置Ａまでγ
線が移動するとき、点Ｘからコリメータ位置Ｂまでγ線
が移動するときよりも、点Ｘから出たガンマ像の受ける
被写体自体による吸収が小さいことを意味する。このた
め、コリメータ位置Ａのデータはコリメータ位置Ｂのデ
ータより正確なカウントで、且つノイズが少い。
の間の距離が、コリメータ位置Ｂ及び点Ｘの間の距離よ
り短いということは、点Ｘからコリメータ位置Ａまでγ
線が移動するとき、点Ｘからコリメータ位置Ｂまでγ線
が移動するときよりも、点Ｘから出たガンマ像の受ける
被写体自体による吸収が小さいことを意味する。このた
め、コリメータ位置Ａのデータはコリメータ位置Ｂのデ
ータより正確なカウントで、且つノイズが少い。
【００２３】以上の２つの理由により、本発明の一実施
例では、ガンマカメラの角度βで、投影像のポイントＰ
（β，Ｓ）が上記（１） に従ってコンボリューション
バックプロジェクションされる条件で、逆投影する範囲
を、図５で示す（ｇ−ａ）／２を半径とする円：無円の
外側のみｌ−ｌ″間のみの範囲、換言すれば再構成画像
上の個々の画素についてコリメータに近い側のＦＯＶ内
とした。
例では、ガンマカメラの角度βで、投影像のポイントＰ
（β，Ｓ）が上記（１） に従ってコンボリューション
バックプロジェクションされる条件で、逆投影する範囲
を、図５で示す（ｇ−ａ）／２を半径とする円：無円の
外側のみｌ−ｌ″間のみの範囲、換言すれば再構成画像
上の個々の画素についてコリメータに近い側のＦＯＶ内
とした。
【００２４】ところが、図３の関係においてコリメータ
位置Ａ，Ｂどちらかも点Ｘが同じ位の距離の場合、少し
でも近い方のデータを使ったからといって、空間分解能
が少し遠い方のデータを使ったときよりもそれほど向上
するわけではない。むしろ。コリメータ位置Ａ，Ｂの各
データがほぼ同等の情報量であることから、両方のデー
タを使った方がノイズの平均化をなし得るので有利であ
る。そして、その両方のデータを使った場合には、点Ｘ
に於けるノイズが約１／２１／２ になる効果がある。
位置Ａ，Ｂどちらかも点Ｘが同じ位の距離の場合、少し
でも近い方のデータを使ったからといって、空間分解能
が少し遠い方のデータを使ったときよりもそれほど向上
するわけではない。むしろ。コリメータ位置Ａ，Ｂの各
データがほぼ同等の情報量であることから、両方のデー
タを使った方がノイズの平均化をなし得るので有利であ
る。そして、その両方のデータを使った場合には、点Ｘ
に於けるノイズが約１／２１／２ になる効果がある。
【００２５】もし、このような漸進的な接続をせず、図
６に示すような不連続をつなぎ方（ある点Ｘではコリメ
ータ位置Ａのデータだけ、点Ｘのすぐ近く点ではＢのデ
ータだけとした関係でバックプロジェクション）をする
と、そのつなぎ目のところは段差が生じることがおこり
得る。このような段差をアーチファクトの原因となり画
質を劣化させるおそれがある。
６に示すような不連続をつなぎ方（ある点Ｘではコリメ
ータ位置Ａのデータだけ、点Ｘのすぐ近く点ではＢのデ
ータだけとした関係でバックプロジェクション）をする
と、そのつなぎ目のところは段差が生じることがおこり
得る。このような段差をアーチファクトの原因となり画
質を劣化させるおそれがある。
【００２６】次に、上記した漸進的な接続をする場合の
一例を図７，図８に従って説明する。
一例を図７，図８に従って説明する。
【００２７】図７に示すように無円ＦＯを中心として直
径がＤｍになるように縮小した円ＦＤｍと直径ＦＤｊに
なるように拡大した円ＦＤｊを想定する。但し、｜ＯＤ
ｍ｜＝｜ＤｊＯ｜＝ｋとする。
径がＤｍになるように縮小した円ＦＤｍと直径ＦＤｊに
なるように拡大した円ＦＤｊを想定する。但し、｜ＯＤ
ｍ｜＝｜ＤｊＯ｜＝ｋとする。
【００２８】そして、ＦのビームＬがこれらの円と交差
する点をそれぞれＸ，Ｘｍ，Ｘｊとする。
する点をそれぞれＸ，Ｘｍ，Ｘｊとする。
【００２９】次に、線Ｌに沿って、図８に示すように重
みｗをつける。この重みは、ファン逆投影における重み
（Ｆからの距離の逆２乗）とは別個であり、ファン逆投
影における重みに更にｗを掛算するのである。ここで、
ｗはＸｍ〜Ｘｊにおいて、０から１へ滑らかにつながる
ようになっていて、しかも、 ｗ（ｄ）＋ｗ（１ＦＦ′−ｄ）＝２ でなければならない。これを満たすためには、ｄ＝Ｘｍ
〜Ｘｊにおいて例えば次の３次式を使うとよい。
みｗをつける。この重みは、ファン逆投影における重み
（Ｆからの距離の逆２乗）とは別個であり、ファン逆投
影における重みに更にｗを掛算するのである。ここで、
ｗはＸｍ〜Ｘｊにおいて、０から１へ滑らかにつながる
ようになっていて、しかも、 ｗ（ｄ）＋ｗ（１ＦＦ′−ｄ）＝２ でなければならない。これを満たすためには、ｄ＝Ｘｍ
〜Ｘｊにおいて例えば次の３次式を使うとよい。
【００３０】つまり、ｗ（ｄ）は、
ｄ＜Ｄ／２−ｕ→０
ｄ＞Ｄ／２＋ｕ→１
ｏｔｈｅｒ ｗｉｓｅ→１／２［−（ｄ−Ｄ／２）
３ ＋３（ｄ−Ｄ／２）］＋１の各条件を満すものとす
る。
３ ＋３（ｄ−Ｄ／２）］＋１の各条件を満すものとす
る。
【００３１】これによって図６の点Ａでの不連続が回避
され、アーチファクトが無くなるのみならず画像中心部
分でのノイズが下る。
され、アーチファクトが無くなるのみならず画像中心部
分でのノイズが下る。
【００３２】
【発明の効果】以上説明したように本発明のＳＰＥＣＴ
画像の再構成方法によれば、逆投影する範囲をファンビ
ームコリメータに近い側の有効視野内としてＳＰＥＣＴ
画像を再構成するから、ＳＰＥＣＴ画像の解像力が従来
に比較して大幅に向上される。
画像の再構成方法によれば、逆投影する範囲をファンビ
ームコリメータに近い側の有効視野内としてＳＰＥＣＴ
画像を再構成するから、ＳＰＥＣＴ画像の解像力が従来
に比較して大幅に向上される。
【図１】本発明のＳＰＥＣＴ画像の再構成方法が適用さ
れたＳＰＥＣＴ装置の概略を示す図である。
れたＳＰＥＣＴ装置の概略を示す図である。
【図２】ファンビームコリメータを用いてＳＰＥＣＴ画
像を再構成する原理を示す図である。
像を再構成する原理を示す図である。
【図３】再構成画像上の点とコリメータ位置とを示す図
である。
である。
【図４】再構成画像上の任意の点のＲＩ濃度とコリメー
タ位置のデータとの相関を説明するために用いた図であ
る。
タ位置のデータとの相関を説明するために用いた図であ
る。
【図５】逆投影する範囲を再構成画像上の個々の画素に
ついてファンビームコリメータに近い側の有効視野内と
してＳＰＥＣＴ画像を再構成する概略を示す図である。
ついてファンビームコリメータに近い側の有効視野内と
してＳＰＥＣＴ画像を再構成する概略を示す図である。
【図６】逆投影する範囲を再構成画像上の個々の画素に
ついてファンビームコリメータに近い側の有効視野内分
とした際に生じるデータの不連続について説明するため
に用いた図である。
ついてファンビームコリメータに近い側の有効視野内分
とした際に生じるデータの不連続について説明するため
に用いた図である。
【図７】データの不連続を補間する方法を説明するため
に用いた図である。
に用いた図である。
【図８】データの不連続を補間する具体例を説明するた
めに用いた図である。
めに用いた図である。
Claims (1)
- 【請求項１】 被写体内に分布した放射性物質が放出
されるγ線をファンビーム状に異なる方向から検出し、
この検出により収集した投影データをフィルタ処理によ
り補正後、逆投影してＳＰＥＣＴ画像を再構成する際、
逆投影する範囲を、再構成画像上の個々の画素について
ファンビームコリメータに近い側の有効視野内とするこ
とを特徴とするＳＰＥＣＴ画像の再構成方法。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP03147596A JP3089050B2 (ja) | 1991-06-19 | 1991-06-19 | Ｓｐｅｃｔ画像の再構成方法 |
| US07/899,693 US5324946A (en) | 1991-06-19 | 1992-06-17 | Method and apparatus for reconstructing spect image with maintaining high spatial resolution |
| EP92110261A EP0519445B1 (en) | 1991-06-19 | 1992-06-17 | Method and apparatus for reconstructing spect image with maintaining high spatial resolution |
| DE69212723T DE69212723T2 (de) | 1991-06-19 | 1992-06-17 | Verfahren und Vorrichtung zur Spect-Bildzusammensetzung unter Aufrechterhaltung einer hohen Raumauflösung |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP03147596A JP3089050B2 (ja) | 1991-06-19 | 1991-06-19 | Ｓｐｅｃｔ画像の再構成方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH04370784A true JPH04370784A (ja) | 1992-12-24 |
| JP3089050B2 JP3089050B2 (ja) | 2000-09-18 |
Family
ID=15433926
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP03147596A Expired - Fee Related JP3089050B2 (ja) | 1991-06-19 | 1991-06-19 | Ｓｐｅｃｔ画像の再構成方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5324946A (ja) |
| EP (1) | EP0519445B1 (ja) |
| JP (1) | JP3089050B2 (ja) |
| DE (1) | DE69212723T2 (ja) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5420803A (en) * | 1992-11-18 | 1995-05-30 | Technology Development Group, Inc. | Enhanced resolution wafer thickness measurement system |
| US5434416A (en) * | 1993-03-24 | 1995-07-18 | Kabushiki Kaisha Toshiba | Method and apparatus for reconstructing SPECT image by utilizing two separate reconstruction filter functions |
| US5528042A (en) * | 1995-06-14 | 1996-06-18 | Siemens Medical Systems, Inc. | Retrospectively determining the center of rotation of a scintillation camera detector from SPECT data acquired during a nuclear medicine study |
| US6452183B1 (en) | 1999-09-13 | 2002-09-17 | Kabushiki Kaisha Toshiba | Nuclear medical diagnosis apparatus and image reconstruction method therefor |
| JP3909048B2 (ja) * | 2003-09-05 | 2007-04-25 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｘ線ｃｔ装置およびｘ線管 |
| ES2292327B1 (es) * | 2005-12-26 | 2009-04-01 | Consejo Superior Investigaciones Cientificas | Mini camara gamma autonoma y con sistema de localizacion, para uso intraquirurgico. |
| US20080008372A1 (en) * | 2006-07-07 | 2008-01-10 | General Electric Company | A method and system for reducing artifacts in a tomosynthesis imaging system |
| CN103163548B (zh) * | 2013-03-07 | 2016-01-20 | 北京永新医疗设备有限公司 | 基于伽马相机的放射性物质探测方法及其装置和系统 |
| WO2017096609A1 (en) | 2015-12-11 | 2017-06-15 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image reconstruction |
| US10255696B2 (en) | 2015-12-11 | 2019-04-09 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image reconstruction |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4580219A (en) * | 1983-05-02 | 1986-04-01 | General Electric Company | Method for reducing image artifacts due to projection measurement inconsistencies |
| US4593355A (en) * | 1983-11-21 | 1986-06-03 | American Science And Engineering, Inc. | Method of quick back projection for computed tomography and improved CT machine employing the method |
| US4729100A (en) * | 1984-08-28 | 1988-03-01 | Kabushiki Kaisha Toshiba | CT System which convolutes projection data with a frequency varying filter function |
| US4831261A (en) * | 1986-06-20 | 1989-05-16 | Digital Scintigraphics, Inc. | Compound collimator and tomography camera using same |
| JP2957623B2 (ja) * | 1990-01-31 | 1999-10-06 | 株式会社東芝 | Ｓｐｅｃｔ装置用のアーチファクト乃至感度補正装置 |
-
1991
- 1991-06-19 JP JP03147596A patent/JP3089050B2/ja not_active Expired - Fee Related
-
1992
- 1992-06-17 EP EP92110261A patent/EP0519445B1/en not_active Expired - Lifetime
- 1992-06-17 DE DE69212723T patent/DE69212723T2/de not_active Expired - Fee Related
- 1992-06-17 US US07/899,693 patent/US5324946A/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP3089050B2 (ja) | 2000-09-18 |
| DE69212723D1 (de) | 1996-09-19 |
| DE69212723T2 (de) | 1997-03-13 |
| US5324946A (en) | 1994-06-28 |
| EP0519445B1 (en) | 1996-08-14 |
| EP0519445A1 (en) | 1992-12-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2394494C2 (ru) | Система сканирования х-кт | |
| US8787646B2 (en) | Method and apparatus for computed tomography image reconstruction | |
| US7424089B2 (en) | System and method for reconstructing image by using straight-line trajectory scan | |
| US7903779B2 (en) | Apparatus and method for reconstruction of volumetric images in a divergent scanning computed tomography system | |
| JP5019193B2 (ja) | スキャン対象に関するｃｔ画像内の再構成点における画像データ値を決定する再構成方法及びｘ線コンピュータ断層撮影装置 | |
| US6907100B2 (en) | Cone beam type of X-ray CT system for three-dimensional reconstruction | |
| JPH0919425A (ja) | Ｘ線コンピュータ断層撮影装置 | |
| JP2010501856A (ja) | 動きアーティファクト画像についてのアーティファクト補正 | |
| JP2008520257A (ja) | ボクセルに依存する補間を用いる画像再構成 | |
| JP3089050B2 (ja) | Ｓｐｅｃｔ画像の再構成方法 | |
| JPH0720245A (ja) | ポジトロンｃｔ装置 | |
| US4692624A (en) | Emission computed tomography apparatus | |
| US5434416A (en) | Method and apparatus for reconstructing SPECT image by utilizing two separate reconstruction filter functions | |
| EP1921467B1 (en) | Tomograph, tomography, tomography program, and computer-readable recording medium where the program is recorded | |
| JP3583567B2 (ja) | Ｘ線断層撮像装置 | |
| JPH119582A (ja) | Ｘ線コンピュータ断層撮影装置 | |
| EP1386289A2 (en) | Fast transform for reconstruction of rotating-slat data | |
| JP2006133204A (ja) | Ｓｐｅｃｔ装置 | |
| JP4387758B2 (ja) | Ｓｐｅｃｔ装置及びｓｐｅｃｔ画像再構成方法 | |
| JP3352215B2 (ja) | Ｓｐｅｃｔ装置及びｓｐｅｃｔ画像の再構成方法 | |
| JP2003010169A (ja) | コンピュータ断層撮影装置 | |
| JP2002204796A (ja) | ３次元ｘ線ｃｔ装置 | |
| JP3490699B2 (ja) | Ｘ線コンピュータ断層撮影装置 | |
| JP2004524943A (ja) | フィルターを伴ったコンピューター断層撮影システム | |
| JP3101308B2 (ja) | Ｓｐｅｃｔ装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |