JPH0630927A - 放射線ｃｔの高速逆投影演算を行う画像再構成方法と画像再構成装置 - Google Patents
放射線ｃｔの高速逆投影演算を行う画像再構成方法と画像再構成装置Info
- Publication number
- JPH0630927A JPH0630927A JP4190762A JP19076292A JPH0630927A JP H0630927 A JPH0630927 A JP H0630927A JP 4190762 A JP4190762 A JP 4190762A JP 19076292 A JP19076292 A JP 19076292A JP H0630927 A JPH0630927 A JP H0630927A
- Authority
- JP
- Japan
- Prior art keywords
- data
- backprojection
- image
- memory
- data memory
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 13
- 230000005855 radiation Effects 0.000 title claims description 6
- 230000015654 memory Effects 0.000 claims abstract description 47
- 230000003111 delayed effect Effects 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 25
- 238000013480 data collection Methods 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 16
- 239000002872 contrast media Substances 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
(57)【要約】
【目的】 連続的にスキャンを行って多くのイメージを
得る場合にも費用の増大を招くことはなく、多くのイメ
ージを短時間で得るための高速逆投影画像再構成方法と
その装置を実現することである。 【構成】 収集データメモリ２１のデータを第１逆投影
演算器２２でスキャン開始時からのデータを逐次逆投影
演算を行って第１逆投影データメモリ２３に逐次格納す
る。第２逆投影演算器２４は３６０°遅れて収集データ
メモリ２１に格納されているスキャン開始時点からのデ
ータを逐次逆投影演算を行って第２逆投影データメモリ
２５に格納する。上記各逆投影データメモリ２３，２５
からのデータは減算器２６で減算されて任意のスライス
位置における逆投影データを得ることができる。
得る場合にも費用の増大を招くことはなく、多くのイメ
ージを短時間で得るための高速逆投影画像再構成方法と
その装置を実現することである。 【構成】 収集データメモリ２１のデータを第１逆投影
演算器２２でスキャン開始時からのデータを逐次逆投影
演算を行って第１逆投影データメモリ２３に逐次格納す
る。第２逆投影演算器２４は３６０°遅れて収集データ
メモリ２１に格納されているスキャン開始時点からのデ
ータを逐次逆投影演算を行って第２逆投影データメモリ
２５に格納する。上記各逆投影データメモリ２３，２５
からのデータは減算器２６で減算されて任意のスライス
位置における逆投影データを得ることができる。
Description
【０００１】
【産業上の利用分野】本発明は放射線ＣＴの高速逆投影
演算を行う画像再構成装置とその方法に関する。
演算を行う画像再構成装置とその方法に関する。
【０００２】
【従来の技術】ＣＴにおいて、被検体の周囲３６０°に
亘ってデータを収集し、これらのデータから画像再構成
することにより、被検体の断層像を得て表示することが
行われている。この画像再構成のアルゴリズムは大きく
分類して、逆投影法、逐次近似法、解析的方法とに分け
られる。
亘ってデータを収集し、これらのデータから画像再構成
することにより、被検体の断層像を得て表示することが
行われている。この画像再構成のアルゴリズムは大きく
分類して、逆投影法、逐次近似法、解析的方法とに分け
られる。
【０００３】ここで、従来の逆投影法について説明す
る。この逆投影演算では通常３６０°毎に計算すること
によって画像データを得ている（３６０°毎に画像再構
成することをフルリコンといっており、１８０°毎に行
う場合をハーフリコンといっている）。
る。この逆投影演算では通常３６０°毎に計算すること
によって画像データを得ている（３６０°毎に画像再構
成することをフルリコンといっており、１８０°毎に行
う場合をハーフリコンといっている）。
【０００４】
【発明が解決しようとする課題】今、ヘリカルスキャン
を行う場合に、或るデーブル速度によって１秒間のスキ
ャンに対しテーブルの進む距離よりも小さいイメージピ
ッチを設定すると、重なった３６０°のデータが作成さ
れ、それを用いて逆投影演算によってイメージが得られ
る。
を行う場合に、或るデーブル速度によって１秒間のスキ
ャンに対しテーブルの進む距離よりも小さいイメージピ
ッチを設定すると、重なった３６０°のデータが作成さ
れ、それを用いて逆投影演算によってイメージが得られ
る。
【０００５】図６はテーブル速度を１０mm/sとし、イメ
ージピッチを２．５mmとした時の１回転のスキャン時間
が１秒のヘリカルスキャンによって得られるデータの図
である。図において、（イ）図は連続回転により得られ
るデータと、画像再構成しようとするスライス位置にお
けるデータを示す図である。この例では、１回転のスキ
ャン時間が１秒で、テーブル速度が１０mm/sなので、
２．５mm間隔のスライス位置は９０°間隔になってい
る。（ロ）図は各スライス位置における収集データを示
している。即ち、スライス位置０mmにおいては、スキャ
ンの０°から３６０°までの３６０°の範囲のデータセ
ット１によりイメージNo.1の画像データが得られる。ス
ライス位置２．５mmでは、（イ）図のスキャンの９０°
から４５０°に至る３６０°の範囲のデータを０°〜３
６０°のデータセット２としてイメージNo.2の画像デー
タを得ている。同様にしてデータセット５までのデータ
を画像再構成して、イメージNo.3，イメージNo.4，イメ
ージNo.5の画像データを得ている。このように、２回転
のヘリカルスキャンによるデータによって、５個のデー
タセットが得られ、５個のイメージが５回の逆投影演算
を行うことによって得られる。
ージピッチを２．５mmとした時の１回転のスキャン時間
が１秒のヘリカルスキャンによって得られるデータの図
である。図において、（イ）図は連続回転により得られ
るデータと、画像再構成しようとするスライス位置にお
けるデータを示す図である。この例では、１回転のスキ
ャン時間が１秒で、テーブル速度が１０mm/sなので、
２．５mm間隔のスライス位置は９０°間隔になってい
る。（ロ）図は各スライス位置における収集データを示
している。即ち、スライス位置０mmにおいては、スキャ
ンの０°から３６０°までの３６０°の範囲のデータセ
ット１によりイメージNo.1の画像データが得られる。ス
ライス位置２．５mmでは、（イ）図のスキャンの９０°
から４５０°に至る３６０°の範囲のデータを０°〜３
６０°のデータセット２としてイメージNo.2の画像デー
タを得ている。同様にしてデータセット５までのデータ
を画像再構成して、イメージNo.3，イメージNo.4，イメ
ージNo.5の画像データを得ている。このように、２回転
のヘリカルスキャンによるデータによって、５個のデー
タセットが得られ、５個のイメージが５回の逆投影演算
を行うことによって得られる。
【０００６】一般に１回の画像再構成に要する時間をｔ
とし、スキャン時間を１秒，テーブル速度をｖmm/s，ピ
クセルサイズをｐとして、ｎ回転のデータを用いて画像
再構成を行うと、得られるイメージの数Ｎは次式のよう
になる。
とし、スキャン時間を１秒，テーブル速度をｖmm/s，ピ
クセルサイズをｐとして、ｎ回転のデータを用いて画像
再構成を行うと、得られるイメージの数Ｎは次式のよう
になる。
【０００７】Ｎ＝（ｖ／ｐ）×（ｎ−１）＋１ 従って画像再構成に要する時間Ｔは次式で得られる。 Ｔ＝Ｎｔ＝｛（ｖ／ｐ）×（ｎ−１）＋１｝ｔ ｖ＝１０mm／ｓ，ｐ＝０．５，ｎ＝１１回転，ｔ＝１０
秒とすると、 Ｔ＝｛（１０／０．５）×（１１−１）＋１｝×１０ ＝２０１０ｓ＝３３．５min となる。即ち、ピクセルサイズと同じイメージ間ピッチ
で、１０cmの間のスライスデータにより３Ｄ用のイメー
ジを得るためには３３．５分の時間が必要となり、実用
的でない。これはヘリカルスキャンに限らず、同一位置
における連続スキャンにより、例えば、造影剤の動きを
見るような場合にも同様である。
秒とすると、 Ｔ＝｛（１０／０．５）×（１１−１）＋１｝×１０ ＝２０１０ｓ＝３３．５min となる。即ち、ピクセルサイズと同じイメージ間ピッチ
で、１０cmの間のスライスデータにより３Ｄ用のイメー
ジを得るためには３３．５分の時間が必要となり、実用
的でない。これはヘリカルスキャンに限らず、同一位置
における連続スキャンにより、例えば、造影剤の動きを
見るような場合にも同様である。
【０００８】この所要時間を減少するためには、各デー
タセット毎に逆投影演算用のＣＰＵを並列に設置するこ
とが考えられるが、コストが増大して実用には限界があ
る。本発明は上記の点に鑑みてなされたもので、その目
的は、連続的にスキャンを行って多くのイメージを得る
場合においても、費用の増大を制限しながら多くのイメ
ージを短時間で得るための高速逆投影画像再構成方法と
その装置を実現することにある。
タセット毎に逆投影演算用のＣＰＵを並列に設置するこ
とが考えられるが、コストが増大して実用には限界があ
る。本発明は上記の点に鑑みてなされたもので、その目
的は、連続的にスキャンを行って多くのイメージを得る
場合においても、費用の増大を制限しながら多くのイメ
ージを短時間で得るための高速逆投影画像再構成方法と
その装置を実現することにある。
【０００９】
【課題を解決するための手段】前記の課題を解決する第
１の発明は、連続スキャンにより得たデータを収集して
前処理後収集データメモリに格納し、このデータを逐次
読み出して第１の逆投影演算を行い、次に、１画像再構
成期間だけ遅れて前記収集データメモリに格納されてい
るデータ収集開始時からのデータについて第２の逆投影
演算を行い、前記第１の逆投影演算の結果のデータから
第２の逆投影演算の結果のデータを減じて必要な１画像
再構成期間の範囲の画像データを得ることを特徴とする
ものである。
１の発明は、連続スキャンにより得たデータを収集して
前処理後収集データメモリに格納し、このデータを逐次
読み出して第１の逆投影演算を行い、次に、１画像再構
成期間だけ遅れて前記収集データメモリに格納されてい
るデータ収集開始時からのデータについて第２の逆投影
演算を行い、前記第１の逆投影演算の結果のデータから
第２の逆投影演算の結果のデータを減じて必要な１画像
再構成期間の範囲の画像データを得ることを特徴とする
ものである。
【００１０】第２の発明は、連続回転により収集したデ
ータに逆投影法を用いて画像再構成を行う放射ＣＴの画
像再構成装置において、多回転スキャンにより得られ前
処理を受けた収集データを逐次格納する収集データメモ
リと、該収集データメモリから読み出されたデータを逐
次加算する第１の逆投影演算を行い、更に１画像再構成
期間遅れて前記収集データメモリに格納されているスキ
ャン開始時点からのデータを読み出して第２の逆投影演
算を行う逆投影演算手段と、前記第１の逆投影演算結果
のデータを格納する第１の逆投影データメモリと、前記
第２の逆投影演算の結果のデータを格納する第２の逆投
影データメモリと、前記第１の逆投影データメモリから
読み出されたデータと、前記第２の逆投影データメモリ
から読み出されたデータとの差を求める減算器とを具備
することを特徴とするものである。
ータに逆投影法を用いて画像再構成を行う放射ＣＴの画
像再構成装置において、多回転スキャンにより得られ前
処理を受けた収集データを逐次格納する収集データメモ
リと、該収集データメモリから読み出されたデータを逐
次加算する第１の逆投影演算を行い、更に１画像再構成
期間遅れて前記収集データメモリに格納されているスキ
ャン開始時点からのデータを読み出して第２の逆投影演
算を行う逆投影演算手段と、前記第１の逆投影演算結果
のデータを格納する第１の逆投影データメモリと、前記
第２の逆投影演算の結果のデータを格納する第２の逆投
影データメモリと、前記第１の逆投影データメモリから
読み出されたデータと、前記第２の逆投影データメモリ
から読み出されたデータとの差を求める減算器とを具備
することを特徴とするものである。
【００１１】
【作用】連続スキャンによるデータをスキャン開始時点
から収集データメモリに格納し、これをスキャン開始時
点から読み出して逆投影演算手段で第１の逆投影演算を
行って第１の逆投影データメモリに格納し、次に１画像
再構成期間だけ遅れて、前記収集データメモリに格納さ
れているデータ収集開始時点からのデータを逐次逆投影
演算して第２逆投影データメモリに格納する。格納され
た第１の逆投影データメモリのイメージから第２の逆投
影データメモリのイメージデータが減算器で減算され、
任意のスライス位置の逆投影データが出力される。
から収集データメモリに格納し、これをスキャン開始時
点から読み出して逆投影演算手段で第１の逆投影演算を
行って第１の逆投影データメモリに格納し、次に１画像
再構成期間だけ遅れて、前記収集データメモリに格納さ
れているデータ収集開始時点からのデータを逐次逆投影
演算して第２逆投影データメモリに格納する。格納され
た第１の逆投影データメモリのイメージから第２の逆投
影データメモリのイメージデータが減算器で減算され、
任意のスライス位置の逆投影データが出力される。
【００１２】
【実施例】以下、図面を参照して本発明の実施例をＸ線
ＣＴについて詳細に説明する。図１は本発明の一実施例
の画像再構成装置のブロック図、図２はこの画像再構成
装置を用いたＸ線ＣＴのブロック図である。
ＣＴについて詳細に説明する。図１は本発明の一実施例
の画像再構成装置のブロック図、図２はこの画像再構成
装置を用いたＸ線ＣＴのブロック図である。
【００１３】図２において、１はＸ線ＣＴを構成する構
成部分であるガントリで、Ｘ線管２と検出器３を被検体
の周囲で各種のスキャン方式に応じた動作をさせる。４
は被検体５を載置してガントリ１の内部に送り込むため
のテーブルである。ガントリ１のティルト，テーブル４
の移動等はテーブル・ガントリ制御装置６によって制御
される。
成部分であるガントリで、Ｘ線管２と検出器３を被検体
の周囲で各種のスキャン方式に応じた動作をさせる。４
は被検体５を載置してガントリ１の内部に送り込むため
のテーブルである。ガントリ１のティルト，テーブル４
の移動等はテーブル・ガントリ制御装置６によって制御
される。
【００１４】Ｘ線管２はＸ線管駆動発生制御装置７の制
御により回転停止及びＸ線の発生休止を行う。検出器３
は検出器駆動装置９の制御により被検体５の周囲を回転
する。Ｘ線管２の照射により被検体５を透過したＸ線は
検出器３で検出され、データ採集装置１０でデータが採
集される。採集されたデータは画像再構成装置１１に転
送される。 画像再構成装置１１は入力されたデータに
逆投影演算を施こし、イメージを再構成する。再構成さ
れたイメージは表示装置１２に表示され、或るいは、デ
ータ保存装置１３に保存される。
御により回転停止及びＸ線の発生休止を行う。検出器３
は検出器駆動装置９の制御により被検体５の周囲を回転
する。Ｘ線管２の照射により被検体５を透過したＸ線は
検出器３で検出され、データ採集装置１０でデータが採
集される。採集されたデータは画像再構成装置１１に転
送される。 画像再構成装置１１は入力されたデータに
逆投影演算を施こし、イメージを再構成する。再構成さ
れたイメージは表示装置１２に表示され、或るいは、デ
ータ保存装置１３に保存される。
【００１５】データ保存装置１３に保存されたイメージ
データは読み出されて３Ｄ構成装置１４で集団検診に適
した３Ｄ像を形成され、その投影イメージは表示装置１
２で表示される。システム制御装置１５は上記各装置の
動作及びそのタイミングを制御する。
データは読み出されて３Ｄ構成装置１４で集団検診に適
した３Ｄ像を形成され、その投影イメージは表示装置１
２で表示される。システム制御装置１５は上記各装置の
動作及びそのタイミングを制御する。
【００１６】次に、本実施例の画像再構成装置１１の行
う逆投影法による画像再構成の方法を図１の画像再構成
装置１１のブロック図により説明する。図において、図
２と同一の部分には同一の符号を付してある。２１はデ
ータ採集装置１０が採集した複数回転によるスキャンデ
ータを、すべて生のデータのまま格納する収集データメ
モリである。
う逆投影法による画像再構成の方法を図１の画像再構成
装置１１のブロック図により説明する。図において、図
２と同一の部分には同一の符号を付してある。２１はデ
ータ採集装置１０が採集した複数回転によるスキャンデ
ータを、すべて生のデータのまま格納する収集データメ
モリである。
【００１７】２２はデータ採集装置１０が採集したスキ
ャン開始からの生データを収集データメモリ２１から読
み出して次々と加算する第１逆投影演算器で、そのデー
タは第１逆投影データメモリ２３に格納される。
ャン開始からの生データを収集データメモリ２１から読
み出して次々と加算する第１逆投影演算器で、そのデー
タは第１逆投影データメモリ２３に格納される。
【００１８】２４は収集データメモリ２１に格納されて
いるデータの内、１回転目が終った後、第１逆投影演算
器２２で演算されるデータから３６０°（フルリコンの
場合）遅れたデータを逐次加算する第２逆投影演算器
で、この演算結果は第２逆投影データメモリ２５に格納
される。２６は第１逆投影データメモリ２３に格納され
ているデータと第２逆投影データメモリ２５に格納され
ているデータとが入力されて、前者から後者を各同一ス
キャン角度毎に減ずる減算器である。
いるデータの内、１回転目が終った後、第１逆投影演算
器２２で演算されるデータから３６０°（フルリコンの
場合）遅れたデータを逐次加算する第２逆投影演算器
で、この演算結果は第２逆投影データメモリ２５に格納
される。２６は第１逆投影データメモリ２３に格納され
ているデータと第２逆投影データメモリ２５に格納され
ているデータとが入力されて、前者から後者を各同一ス
キャン角度毎に減ずる減算器である。
【００１９】次に上記のように構成された画像再構成装
置１１の動作を図３のタイムチャートを参照しながら説
明する。図３において、（イ）図は、スキャンによりデ
ータ採集装置１０が採集し、収集データメモリ２１に格
納される生データの格納のタイミングを示す図で、時刻
ｔ1 でスキャンを開始し、時刻ｔ2 で２回転目、時刻ｔ
3 で３回転目に入っていることを示している。
置１１の動作を図３のタイムチャートを参照しながら説
明する。図３において、（イ）図は、スキャンによりデ
ータ採集装置１０が採集し、収集データメモリ２１に格
納される生データの格納のタイミングを示す図で、時刻
ｔ1 でスキャンを開始し、時刻ｔ2 で２回転目、時刻ｔ
3 で３回転目に入っていることを示している。
【００２０】（ロ）図は、第１逆投影演算器２２で演算
され、第１逆投影データメモリ２３に格納されるイメー
ジデータＩ1 のタイミングを示す図で、スキャン角度０
°から３６０°に至る１回転目では生データＡのイメー
ジデータＡ1 、３６０°から７２０°に至る２回転目で
は生データＡと生データＢのイメージデータＡ1 とＢ 1
が逐次加算されたＡ1 ＋Ｂ1 ，７２０°か１０８０°で
は生データＡ，Ｂ，ＣのイメージデータＡ1 ，Ｂ1 ，Ｃ
1 が逐次加算されたイメージデータＡ1 ＋Ｂ1＋Ｃ1 が
示されている。
され、第１逆投影データメモリ２３に格納されるイメー
ジデータＩ1 のタイミングを示す図で、スキャン角度０
°から３６０°に至る１回転目では生データＡのイメー
ジデータＡ1 、３６０°から７２０°に至る２回転目で
は生データＡと生データＢのイメージデータＡ1 とＢ 1
が逐次加算されたＡ1 ＋Ｂ1 ，７２０°か１０８０°で
は生データＡ，Ｂ，ＣのイメージデータＡ1 ，Ｂ1 ，Ｃ
1 が逐次加算されたイメージデータＡ1 ＋Ｂ1＋Ｃ1 が
示されている。
【００２１】（ハ）図は、第２逆投影演算器２４で演算
され、第２逆投影データメモリ２５に格納されるイメー
ジデータＩ2 のタイミングを示す図で、第１逆投影演算
器２２による演算開始時刻ｔ1 から１回転分即ち３６０
°遅れて収集データメモリ２１から生データＡを読み出
して逆投影演算を行って第２逆投影データメモリ２５に
格納するイメージデータＩ2 のタイミングを示す図で、
２回転目の始まる時刻ｔ2 から生データＡをイメージデ
ータＡ1 に演算し、時刻ｔ3 から生データＢをイメージ
データＢ1 に演算しながらイメージデータＡ1 に加算し
てイメージデータＡ1 ＋Ｂ1 を算出している状態を示し
ている。
され、第２逆投影データメモリ２５に格納されるイメー
ジデータＩ2 のタイミングを示す図で、第１逆投影演算
器２２による演算開始時刻ｔ1 から１回転分即ち３６０
°遅れて収集データメモリ２１から生データＡを読み出
して逆投影演算を行って第２逆投影データメモリ２５に
格納するイメージデータＩ2 のタイミングを示す図で、
２回転目の始まる時刻ｔ2 から生データＡをイメージデ
ータＡ1 に演算し、時刻ｔ3 から生データＢをイメージ
データＢ1 に演算しながらイメージデータＡ1 に加算し
てイメージデータＡ1 ＋Ｂ1 を算出している状態を示し
ている。
【００２２】（ニ）図は、（ロ）図のイメージデータＩ
1 から（ハ）図のイメージデータＩ 2 を減じたイメージ
データＩ3 の図で、各回転角度におけるデータから、そ
のまま逆投影演算されたイメージデータが得られる状態
を示している。
1 から（ハ）図のイメージデータＩ 2 を減じたイメージ
データＩ3 の図で、各回転角度におけるデータから、そ
のまま逆投影演算されたイメージデータが得られる状態
を示している。
【００２３】図４は図３に示した各データをスキャン角
度を基準として描いた図で、（イ）図、（ロ）図は図３
と同様である。（ハ）図は第１逆投影演算器２２の演算
開始から３６０°遅れてスキャン角度０°の生データＡ
を演算してイメージデータＡ 1 を算出し、スキャン角度
３６０°からスキャン角度α2 迄演算してイメージデー
タＡ1 ＋Ｂ1 を算出する状態を示している。
度を基準として描いた図で、（イ）図、（ロ）図は図３
と同様である。（ハ）図は第１逆投影演算器２２の演算
開始から３６０°遅れてスキャン角度０°の生データＡ
を演算してイメージデータＡ 1 を算出し、スキャン角度
３６０°からスキャン角度α2 迄演算してイメージデー
タＡ1 ＋Ｂ1 を算出する状態を示している。
【００２４】（ニ）図は、減算器２６でイメージデータ
Ｉ1 からイメージデータＩ2 を減算し、スキャン角度α
2 からスキャン角度α1 迄の３６０°に亘るデータＩ3
が出力される状態を示している。
Ｉ1 からイメージデータＩ2 を減算し、スキャン角度α
2 からスキャン角度α1 迄の３６０°に亘るデータＩ3
が出力される状態を示している。
【００２５】図５は図６に示す従来の場合との比較によ
り、同じイメージデータを得る場合の説明図である。図
において、（イ）図は図４の（ロ）図のα1 を４５０°
とし、（ロ）図は図（４）の（ハ）図のα2 を９０°と
した図である。減算器２６から得られる出力は（ハ）図
のイメージデータＩ3 で、これは図６の（ロ）図のデー
タセット２から成るイメージNo.2と同じデータであっ
て、画像再構成するスライス位置を２．５mm位置の逆投
影演算角度は、この例のようにα1 ＝α2 ＋３６０°＝
４５０°となる。
り、同じイメージデータを得る場合の説明図である。図
において、（イ）図は図４の（ロ）図のα1 を４５０°
とし、（ロ）図は図（４）の（ハ）図のα2 を９０°と
した図である。減算器２６から得られる出力は（ハ）図
のイメージデータＩ3 で、これは図６の（ロ）図のデー
タセット２から成るイメージNo.2と同じデータであっ
て、画像再構成するスライス位置を２．５mm位置の逆投
影演算角度は、この例のようにα1 ＝α2 ＋３６０°＝
４５０°となる。
【００２６】以上説明したように本実施例によれば、２
台の逆投影演算器と２個の逆投影データメモリを用いる
ことによって、任意の時刻において、各回転のスキャン
によって得られるデータを逆投影演算して得たイメージ
データが常に求められ、この再構成に要する時間はイメ
ージの枚数とは無関係であって、所要時間はスキャンの
回転の回数で決まるものである。つまり、従来の通常の
スキャンと同じ再構成時間で、極めて小さいピクセル単
位の極めて多数のイメージデータを得ることができ、こ
れを用いて短時間で良い画質の３次元のヘリカルスキャ
ンによるイメージを得ることができる。又、連続回転の
ダイナミックＣＴにおいて、実時間で造影剤の造影効果
を観測することが可能になる。
台の逆投影演算器と２個の逆投影データメモリを用いる
ことによって、任意の時刻において、各回転のスキャン
によって得られるデータを逆投影演算して得たイメージ
データが常に求められ、この再構成に要する時間はイメ
ージの枚数とは無関係であって、所要時間はスキャンの
回転の回数で決まるものである。つまり、従来の通常の
スキャンと同じ再構成時間で、極めて小さいピクセル単
位の極めて多数のイメージデータを得ることができ、こ
れを用いて短時間で良い画質の３次元のヘリカルスキャ
ンによるイメージを得ることができる。又、連続回転の
ダイナミックＣＴにおいて、実時間で造影剤の造影効果
を観測することが可能になる。
【００２７】尚、本発明は上記実施例に限定されるもの
ではない。 実施例では３６０°のフルリコンで説明したが、１
８０°のハーフリコンでも差し支えない。 逆投影演算器は２台の場合で説明したが、１台の逆
投影演算器を用いて時分割で演算し、第１，第２の逆投
影データメモリにそれぞれ格納するようにしてもよい。 上記の説明はＸ線ＣＴについて行ったが、アイソト
ープ等放射性元素を用いたエミッションＣＴにも用いる
ことができる。
ではない。 実施例では３６０°のフルリコンで説明したが、１
８０°のハーフリコンでも差し支えない。 逆投影演算器は２台の場合で説明したが、１台の逆
投影演算器を用いて時分割で演算し、第１，第２の逆投
影データメモリにそれぞれ格納するようにしてもよい。 上記の説明はＸ線ＣＴについて行ったが、アイソト
ープ等放射性元素を用いたエミッションＣＴにも用いる
ことができる。
【００２８】
【発明の効果】以上詳細に説明したように本発明によれ
ば、逆投影演算手段の数を制限して費用の増大を抑えな
がら多くのイメージを短時間で得ることができるように
なり、実用上の効果は大きい。
ば、逆投影演算手段の数を制限して費用の増大を抑えな
がら多くのイメージを短時間で得ることができるように
なり、実用上の効果は大きい。
【図１】本発明の一実施例の画像再構成装置の構成のブ
ロック図である。
ロック図である。
【図２】本発明の画像再構成装置を用いたＸ線ＣＴのブ
ロック図である。
ロック図である。
【図３】画像再構成装置の逆投影演算のタイムチャート
である。
である。
【図４】画像再構成装置のスキャン角度を基準とした逆
投影演算のタイムチャートである。
投影演算のタイムチャートである。
【図５】スライス位置２．５mmのイメージデータを得る
場合の動作のタイムチャートである。
場合の動作のタイムチャートである。
【図６】従来のＸ線ＣＴのヘリカルスキャンによるスラ
イス間隔の２．５mmの場合のデータの図である。
イス間隔の２．５mmの場合のデータの図である。
１１ 画像再構成装置 ２１ 収集データメモリ ２２ 第１逆投影演算器 ２３ 第１逆投影データメモリ ２４ 第２逆投影演算器 ２５ 第２逆投影データメモリ ２６ 減算器
Claims (2)
- 【請求項１】 連続スキャンにより得たデータを収集し
て前処理後収集データメモリに格納し、このデータを逐
次読み出して第１の逆投影演算を行い、次に、１画像再
構成期間だけ遅れて前記収集データメモリに格納されて
いるデータ収集開始時からのデータについて第２の逆投
影演算を行い、前記第１の逆投影演算の結果のデータか
ら第２の逆投影演算の結果のデータを減じて必要な１画
像再構成期間の範囲の画像データを得ることを特徴とす
る放射線ＣＴの高速逆投影演算を行う画像再構成方法。 - 【請求項２】 連続回転により収集したデータに逆投影
法を用いて画像再構成を行う放射ＣＴの画像再構成装置
において、 多回転スキャンにより得られ前処理を受けた収集データ
を逐次格納する収集データメモリ（２１）と、 該収集データメモリ（２１）から読み出されたデータを
逐次加算する第１の逆投影演算を行い、更に１画像再構
成期間遅れて前記収集データメモリ（２１）に格納され
ているスキャン開始時点からのデータを読み出して第２
の逆投影演算を行う逆投影演算手段（２２，２４）と、 前記第１の逆投影演算結果のデータを格納する第１の逆
投影データメモリ（２３）と、 前記第２の逆投影演算の結果のデータを格納する第２の
逆投影データメモリ（２５）と前記第１の逆投影データ
メモリ（２３）から読み出されたデータと、前記第２の
逆投影データメモリ（２５）から読み出されたデータと
の差を求める減算器（２６）とを具備することを特徴と
する放射線ＣＴの高速逆投影演算を行う画像再構成装
置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4190762A JPH0630927A (ja) | 1992-07-17 | 1992-07-17 | 放射線ｃｔの高速逆投影演算を行う画像再構成方法と画像再構成装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4190762A JPH0630927A (ja) | 1992-07-17 | 1992-07-17 | 放射線ｃｔの高速逆投影演算を行う画像再構成方法と画像再構成装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH0630927A true JPH0630927A (ja) | 1994-02-08 |
Family
ID=16263307
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP4190762A Pending JPH0630927A (ja) | 1992-07-17 | 1992-07-17 | 放射線ｃｔの高速逆投影演算を行う画像再構成方法と画像再構成装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0630927A (ja) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003528683A (ja) * | 2000-03-30 | 2003-09-30 | ザ、ボード、オブ、トラスティーズ、オブ、ザ、ユニバシティー、オブ、イリノイ | ３ｄラドン変換のための高速階層的逆投影方法 |
| JP2009500114A (ja) * | 2005-07-08 | 2009-01-08 | ウイスコンシン アラムナイ リサーチ フオンデーシヨン | 高度に限定されたイメージの再構成法 |
| JP2009500115A (ja) * | 2005-07-08 | 2009-01-08 | ウイスコンシン アラムナイ リサーチ フオンデーシヨン | Ｃｔ画像の逆投影再構成法 |
| JP2013166018A (ja) * | 2007-07-25 | 2013-08-29 | Toshiba Corp | Ｘ線コンピュータ断層撮影装置 |
-
1992
- 1992-07-17 JP JP4190762A patent/JPH0630927A/ja active Pending
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003528683A (ja) * | 2000-03-30 | 2003-09-30 | ザ、ボード、オブ、トラスティーズ、オブ、ザ、ユニバシティー、オブ、イリノイ | ３ｄラドン変換のための高速階層的逆投影方法 |
| JP2009500114A (ja) * | 2005-07-08 | 2009-01-08 | ウイスコンシン アラムナイ リサーチ フオンデーシヨン | 高度に限定されたイメージの再構成法 |
| JP2009500115A (ja) * | 2005-07-08 | 2009-01-08 | ウイスコンシン アラムナイ リサーチ フオンデーシヨン | Ｃｔ画像の逆投影再構成法 |
| JP2013166018A (ja) * | 2007-07-25 | 2013-08-29 | Toshiba Corp | Ｘ線コンピュータ断層撮影装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6256370B1 (en) | Method and apparatus for performing tomosynthesis | |
| US8731269B2 (en) | Method and system for substantially reducing artifacts in circular cone beam computer tomography (CT) | |
| RU2510080C2 (ru) | Устройство для обработки изображения, способ обработки изображения и среда долговременного хранения информации | |
| CN101454804B (zh) | 动态计算机断层摄影成像 | |
| JPH08509408A (ja) | 円すい状ビームデータからの画像の再構成 | |
| JPH03188832A (ja) | 断層撮影像作成方法および装置 | |
| US7209535B2 (en) | Fourier space tomographic image reconstruction method | |
| US4495645A (en) | Computed tomography device | |
| JPH08263638A (ja) | 物体の断層写真像を作成するシステム及び方法 | |
| JPH08308824A (ja) | 物体の断層写真像を作成するシステム | |
| US5438602A (en) | Correction of CT attenuation data using fan beam reprojections | |
| CN105118039A (zh) | 实现锥束ct图像重建的方法及系统 | |
| JP3917684B2 (ja) | 物体の断層写真像を作成する方法及び装置 | |
| JP3006722B2 (ja) | コンピュータ断層撮影装置 | |
| JP3878788B2 (ja) | トモシンセシスを行う方法及び装置 | |
| JPH0798039B2 (ja) | コンピュータ断層撮影装置 | |
| US6999550B2 (en) | Method and apparatus for obtaining data for reconstructing images of an object | |
| JP2000083948A (ja) | 断層写真画像を作る装置と方法 | |
| JP3285944B2 (ja) | コンピュータ断層撮影装置 | |
| JP3529874B2 (ja) | Ｘ線コンピュータ断層撮影装置 | |
| JPH11113892A (ja) | 画像再構成装置 | |
| JP3419116B2 (ja) | Ｘ線ｃｔ装置 | |
| JP3332087B2 (ja) | Ｘ線ｃｔ装置 | |
| JPH0531105A (ja) | Ｘ線ｃｔ装置 | |
| JPH0824253A (ja) | Ｘ線ｃｔ装置 |