JPH0458336B2 - - Google Patents
Info
- Publication number
- JPH0458336B2 JPH0458336B2 JP59156976A JP15697684A JPH0458336B2 JP H0458336 B2 JPH0458336 B2 JP H0458336B2 JP 59156976 A JP59156976 A JP 59156976A JP 15697684 A JP15697684 A JP 15697684A JP H0458336 B2 JPH0458336 B2 JP H0458336B2
- Authority
- JP
- Japan
- Prior art keywords
- view
- data
- view data
- scan
- channel
- 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.)
- Expired - Lifetime
Links
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
【発明の詳細な説明】
(産業上の利用分野)
本発明は、スキヤンの開始点近傍のビユー範囲
内及び終了点近傍のビユー範囲内の各方向のビユ
ーデータを得るに際して、当該方向に関する実測
によるスキヤン・データと、当該方向とは180゜±
βずれた範囲内の方向のビユーデータ群より算出
した当該方向に関する算出ビユーデータとを用
い、体動によるアーテイフアクト(Artifact)を
軽減化した各方向の補正ビユーデータを生成し、
画像再構成時に、この補正ビユーデータを、前記
スキヤンの開始点近傍のビユー範囲内及び終了点
近傍のビユー範囲内の当該方向のビユーデータと
して使用するCT(Computer Tomography)装
置に関する。
内及び終了点近傍のビユー範囲内の各方向のビユ
ーデータを得るに際して、当該方向に関する実測
によるスキヤン・データと、当該方向とは180゜±
βずれた範囲内の方向のビユーデータ群より算出
した当該方向に関する算出ビユーデータとを用
い、体動によるアーテイフアクト(Artifact)を
軽減化した各方向の補正ビユーデータを生成し、
画像再構成時に、この補正ビユーデータを、前記
スキヤンの開始点近傍のビユー範囲内及び終了点
近傍のビユー範囲内の当該方向のビユーデータと
して使用するCT(Computer Tomography)装
置に関する。
(従来の技術)
CTでのスキヤン・タイムはかなり長く(例え
ば第3世代型で数秒〜10数秒かかる)、その間撮
影収集されるデータは理想的静止状態の各方向の
投影データと異なり、動きのある(体動を伴つ
た)撮影データとなる。特に開始と終了点でのデ
ータは時間的に最も隔たりのあるデータとなる。
即ち、360゜方向のスキヤン・データを使つてイメ
ージ再構成を行う従来の方式では、その開始点
(角度θ0)と終了点(角度θ0±(360゜−Δθ))で
最
も時間的ずれが大きいため、この方向に体動によ
るアーテイフアクトが発生しやすい。これを軽減
するために、オーバーラツプ・スキヤンなしで
360゜以内(180゜+α以上)の各角度方向のスキヤ
ン・データをもとにして体動の補正を行い、体動
によるアーテイフアクトの軽減を行える画像再構
成方式がある。具体的には、スキヤンの開始点近
傍のビユー範囲内及び終了点近傍のビユー範囲内
の各方向のビユーデータを得るに際して、当該方
向に関する実測によるスキヤン・データと、当該
方向とは180゜±βずれた範囲内の方向のビユーデ
ータ群より算出した当該方向に関する算出ビユー
データとを重み付け加算して、体動によるアーテ
イフアクトを軽減化した各方向の補正ビユーデー
タを生成し、画像再構成時に、この補正ビユーデ
ータを、前記スキヤンの開始点近傍のビユー範囲
内及び終了点近傍のビユー範囲内の当該方向のビ
ユーデータとして使用するCT装置がある。
ば第3世代型で数秒〜10数秒かかる)、その間撮
影収集されるデータは理想的静止状態の各方向の
投影データと異なり、動きのある(体動を伴つ
た)撮影データとなる。特に開始と終了点でのデ
ータは時間的に最も隔たりのあるデータとなる。
即ち、360゜方向のスキヤン・データを使つてイメ
ージ再構成を行う従来の方式では、その開始点
(角度θ0)と終了点(角度θ0±(360゜−Δθ))で
最
も時間的ずれが大きいため、この方向に体動によ
るアーテイフアクトが発生しやすい。これを軽減
するために、オーバーラツプ・スキヤンなしで
360゜以内(180゜+α以上)の各角度方向のスキヤ
ン・データをもとにして体動の補正を行い、体動
によるアーテイフアクトの軽減を行える画像再構
成方式がある。具体的には、スキヤンの開始点近
傍のビユー範囲内及び終了点近傍のビユー範囲内
の各方向のビユーデータを得るに際して、当該方
向に関する実測によるスキヤン・データと、当該
方向とは180゜±βずれた範囲内の方向のビユーデ
ータ群より算出した当該方向に関する算出ビユー
データとを重み付け加算して、体動によるアーテ
イフアクトを軽減化した各方向の補正ビユーデー
タを生成し、画像再構成時に、この補正ビユーデ
ータを、前記スキヤンの開始点近傍のビユー範囲
内及び終了点近傍のビユー範囲内の当該方向のビ
ユーデータとして使用するCT装置がある。
(発明が解決しようとする問題点)
この方式は先頭スキヤン・データ群と最終スキ
ヤン・データ群のモーシヨンの影警を軽減する優
れた手法であるが、 (1) 第2図イに示すように、先頭ビユー、最終ビ
ユー方向に多数の直線的なストリーク状のアー
テイフアクトA1が発生する。
ヤン・データ群のモーシヨンの影警を軽減する優
れた手法であるが、 (1) 第2図イに示すように、先頭ビユー、最終ビ
ユー方向に多数の直線的なストリーク状のアー
テイフアクトA1が発生する。
(2) (1)の軽減の目的で同一ビユーデータについて
チヤネル方向のスムージングを加えると、高コ
ントラスト変化部分で、 (a) 第2図ロのように接線方向のアーテイフア
クトA2が発生する。
チヤネル方向のスムージングを加えると、高コ
ントラスト変化部分で、 (a) 第2図ロのように接線方向のアーテイフア
クトA2が発生する。
(b) 第2図ハのように長い直線部分を持つ被撮
影体で、直線状のアーテイフアクトA3が発
生する。
影体で、直線状のアーテイフアクトA3が発
生する。
(c) スムージングに伴うイメージの不自然さが
ある。
ある。
等の問題があつた。
本発明はこのような点に鑑みてなされたもの
で、その目的は、先頭ビユー、最終ビユー方向の
直線的なストリーク状のアーテイフアクトだけで
なく、高コントラスト変化部分での接線方向のア
ーテイフアクトや、長い直線部分を持つ被検体で
の直線状のアーテイフアクトをも除去できると共
に、自然なイメージを生成できるCT装置を提供
することにある。
で、その目的は、先頭ビユー、最終ビユー方向の
直線的なストリーク状のアーテイフアクトだけで
なく、高コントラスト変化部分での接線方向のア
ーテイフアクトや、長い直線部分を持つ被検体で
の直線状のアーテイフアクトをも除去できると共
に、自然なイメージを生成できるCT装置を提供
することにある。
(問題点を解決するための手段)
上記問題点を解決する本発明は、スキヤンの開
始点近傍のビユー範囲内及び終了点近傍のビユー
範囲内の各方向のビユーデータを得るに際して、
当該方向に関する実測によるスキヤン・データ
と、当該方向とは180゜±βずれた範囲内の方向の
ビユーデータ群より算出した当該方向に関する算
出ビユーデータとを用い、体動によるアーテイフ
アクトを軽減化した各方向の補正ビユーデータを
生成し、画像再構成時に、この補正ビユーデータ
を、前記スキヤンの開始点近傍のビユー範囲内及
び終了点近傍のビユー範囲内の当該方向のビユー
データとして使用するCT装置において、第Vビ
ユーデータの第Vチヤネルの実測による前記スキ
ヤン・データをDv(I)、前記180゜±βずれた範囲
内の方向のビユーデータ群より算出した第Vビユ
ーの第Iチヤネルの前記算出ビユーデータをHv
(I)で表わす時、第V−kビユーの第Iチヤネルの
算出ビユーデータHv−k(I)と第V+kビユーの
第Iチヤネルの算出ビユーデータHv+kv(I)(た
だし、=1,2,……)との一組又は異なる複数
個のkに対応する複数組のビユーデータ群の加算
平均であるCvv(I)と、前記第ビユーの第Iチヤ
ネルの実測によるスキヤン・データDvv(I)とを重
み付け加算したものを求める手段を設け、該手段
により求めたものを、第Vビユーの第Iチヤネル
の前記補正ビユーデータとして用いることを特徴
とするものである。
始点近傍のビユー範囲内及び終了点近傍のビユー
範囲内の各方向のビユーデータを得るに際して、
当該方向に関する実測によるスキヤン・データ
と、当該方向とは180゜±βずれた範囲内の方向の
ビユーデータ群より算出した当該方向に関する算
出ビユーデータとを用い、体動によるアーテイフ
アクトを軽減化した各方向の補正ビユーデータを
生成し、画像再構成時に、この補正ビユーデータ
を、前記スキヤンの開始点近傍のビユー範囲内及
び終了点近傍のビユー範囲内の当該方向のビユー
データとして使用するCT装置において、第Vビ
ユーデータの第Vチヤネルの実測による前記スキ
ヤン・データをDv(I)、前記180゜±βずれた範囲
内の方向のビユーデータ群より算出した第Vビユ
ーの第Iチヤネルの前記算出ビユーデータをHv
(I)で表わす時、第V−kビユーの第Iチヤネルの
算出ビユーデータHv−k(I)と第V+kビユーの
第Iチヤネルの算出ビユーデータHv+kv(I)(た
だし、=1,2,……)との一組又は異なる複数
個のkに対応する複数組のビユーデータ群の加算
平均であるCvv(I)と、前記第ビユーの第Iチヤ
ネルの実測によるスキヤン・データDvv(I)とを重
み付け加算したものを求める手段を設け、該手段
により求めたものを、第Vビユーの第Iチヤネル
の前記補正ビユーデータとして用いることを特徴
とするものである。
本発明におけるフアンビーム方式でのビユーデ
ータ生成アルゴリズムは次の通りである。
ータ生成アルゴリズムは次の通りである。
Cv(I)を第Vビユー(View)の第Iチヤネルに
おける新しい対向ビユーデータとすると、新しい
対向ビユーデータ生成のアルゴリズムを、 Cvv(I) =Hv(I)+[(1/2){Hv−1(I) +Hv+1(I)}−Hv(I)] =(1/2){Hv−1(I)+Hv+1(I)} ……(1) とする。
おける新しい対向ビユーデータとすると、新しい
対向ビユーデータ生成のアルゴリズムを、 Cvv(I) =Hv(I)+[(1/2){Hv−1(I) +Hv+1(I)}−Hv(I)] =(1/2){Hv−1(I)+Hv+1(I)} ……(1) とする。
ここで、
Hv(I)
=(1−v′+V′){(1−i′+I′)
・Dv′(I′)+(i′+I′)
・Dv′(I′+1)}
+(v′−V′){(1−i′+I′)
・Dv′+1((I′)+(i′−I′)
・Dv′+1(I′+1)}
但し、
v′=V−(π/Δθ)
−2(I−cc)・Δγ/Δθ
(V=1,2,…,VM)
i′=2・cc−I
(I=1,2,…,Im)
V′=trunc(v′)
I′=trunc(i′)
尚、Dv(I)は第Vビユーの第Iチヤネルにおけ
るスキヤン・データ(ビユーデータ)、Hv(I)は第
Vビユーの第Iチヤネルにおける対向ビユーデー
タ(スキヤン・データより演算により求められ
る)、Δθは隣接2ビユーの角度差(ラジアン)、
Δγは検出器が等角度間隔に配置されるとしたと
きの隣接2チヤネル間の中心角誤差(ラジアン)、
ccは中心投影データのチヤネル番号である(ccw
でΔθ>0、cwでΔθ<0)。
るスキヤン・データ(ビユーデータ)、Hv(I)は第
Vビユーの第Iチヤネルにおける対向ビユーデー
タ(スキヤン・データより演算により求められ
る)、Δθは隣接2ビユーの角度差(ラジアン)、
Δγは検出器が等角度間隔に配置されるとしたと
きの隣接2チヤネル間の中心角誤差(ラジアン)、
ccは中心投影データのチヤネル番号である(ccw
でΔθ>0、cwでΔθ<0)。
次にこのアルゴリズムの対向ビユーデータ生成
を用いた体動補正(画像再構成)を説明する。
を用いた体動補正(画像再構成)を説明する。
本発明によるスキヤン及びデータ収集、アーテ
イフアクト補正用データ生成方法(アルゴリズ
ム)を示す。第3図イ,ロは反時計回り方向
(CCW)と時計回り方向(CW)におけるビユー
の構成を示す。V1s,V1E,V2s,V2Eは各方向で
の投影データ群(各ビユー)を示す。Cは再構成
の中心点(Iso−center)でV1SC―――→は投影データ
の中心の方向とみなすことができる(V1EC―――→,
…も同様)。V1Sは第1ビユーで、時間の経過
(増加)と共に、順次V1E,V2s,V2E,……のよ
うにビユーが変化する。第3図ハは、各ビユーデ
ータの第1ビユーV1sよりの経過時間(又はスキ
ヤン時間)又は、第1ビユーからの相対回転角度
を表わす。第3図ニは、各ビユーデータに対応す
る重み関数Weightを示す。
イフアクト補正用データ生成方法(アルゴリズ
ム)を示す。第3図イ,ロは反時計回り方向
(CCW)と時計回り方向(CW)におけるビユー
の構成を示す。V1s,V1E,V2s,V2Eは各方向で
の投影データ群(各ビユー)を示す。Cは再構成
の中心点(Iso−center)でV1SC―――→は投影データ
の中心の方向とみなすことができる(V1EC―――→,
…も同様)。V1Sは第1ビユーで、時間の経過
(増加)と共に、順次V1E,V2s,V2E,……のよ
うにビユーが変化する。第3図ハは、各ビユーデ
ータの第1ビユーV1sよりの経過時間(又はスキ
ヤン時間)又は、第1ビユーからの相対回転角度
を表わす。第3図ニは、各ビユーデータに対応す
る重み関数Weightを示す。
(1) スキヤン及びデータ収集;第1ビユー〜第
V2Eビユーまで行う。
V2Eビユーまで行う。
(2) アーテイフアクト補正用データの生成;スキ
ヤンで収集したビユーデータの第1ビユーデー
タ群をD1(I)、第2ビユーデータ群をD2(I),…
…で表わし、補正されたビユーデータの第1ビ
ユーデータ群をR1(I)、第2ビユーデータ群を
R2(I),……とし、ビユー番号をVで表わせば
(I=1,2,3,……) (I) V1s≦V≦V1Eに対し Rv(I) =w11・Dv(I)+w21・Cv(I) ……(2) w11+w21=1 ……(2−1) () V1E<V′<V2sに対し Rv(I)=Dv(I) (3) () V2s≦V≦V2Eに対し Rv(I) =w12・Dv(I)+w22・Cv(I) ……(4) w12+w22=1 ……(4−1) () V2E<V≦VMに対し Rv(I)=Cv(I) ……(5) に従つて補正データを生成する。
ヤンで収集したビユーデータの第1ビユーデー
タ群をD1(I)、第2ビユーデータ群をD2(I),…
…で表わし、補正されたビユーデータの第1ビ
ユーデータ群をR1(I)、第2ビユーデータ群を
R2(I),……とし、ビユー番号をVで表わせば
(I=1,2,3,……) (I) V1s≦V≦V1Eに対し Rv(I) =w11・Dv(I)+w21・Cv(I) ……(2) w11+w21=1 ……(2−1) () V1E<V′<V2sに対し Rv(I)=Dv(I) (3) () V2s≦V≦V2Eに対し Rv(I) =w12・Dv(I)+w22・Cv(I) ……(4) w12+w22=1 ……(4−1) () V2E<V≦VMに対し Rv(I)=Cv(I) ……(5) に従つて補正データを生成する。
V1s≦V≦V1Eでは原ビユーデータと対向ビユ
ーデータとの間でオーバーラツプを行い、V1E<
V<V2sでは原ビユーデータをそのまま使用す
る。V2s≦V≦V2Eでは原ビユーデータと対向ビ
ユーデータとの間でオーバラツプを行い、V2E<
V≦VMでは対向ビユーデータを使用する。V1s〜
V1E,V2s〜V2E従つてV1E〜V2s,V2E〜VMは、若
干の制限を除いて、ビユーの位置と大きさは自由
に選ぶことができる。第4図には、画像再構成処
理における手順の例を示してある。
ーデータとの間でオーバーラツプを行い、V1E<
V<V2sでは原ビユーデータをそのまま使用す
る。V2s≦V≦V2Eでは原ビユーデータと対向ビ
ユーデータとの間でオーバラツプを行い、V2E<
V≦VMでは対向ビユーデータを使用する。V1s〜
V1E,V2s〜V2E従つてV1E〜V2s,V2E〜VMは、若
干の制限を除いて、ビユーの位置と大きさは自由
に選ぶことができる。第4図には、画像再構成処
理における手順の例を示してある。
(実施例)
第1図は本発明によるX線CTの構成例である。
図において、テーブル・ガントリ−TGには、X
線発生部XGとX線検出部Sとが装備され、被検
体PA(患者)を寝かせるように構成されている。
この被検体PAにX線発生部XGから公知の手法
により、X線を発生し被検体PAを透過させ、X
線検出部Sで電気信号に変換される。この電気信
号はデータ収集装置DASで増幅、積分され、
A/D変換器によりデイジタル信号に変換して透
視像データを得る。この透視像データはデータ記
憶装置DS1に一時記憶され、前処理装置PPCによ
り、断面像再構成のための前処理を施す。その前
処理の結果得られたデータはデータ記憶装置DS
2に一時記憶される。このデータ記憶装置DS2
に記憶されたデータは、オーバーラツプ処理装置
OVLPで処理され、データ記憶装置DS3に一時
記憶される。
図において、テーブル・ガントリ−TGには、X
線発生部XGとX線検出部Sとが装備され、被検
体PA(患者)を寝かせるように構成されている。
この被検体PAにX線発生部XGから公知の手法
により、X線を発生し被検体PAを透過させ、X
線検出部Sで電気信号に変換される。この電気信
号はデータ収集装置DASで増幅、積分され、
A/D変換器によりデイジタル信号に変換して透
視像データを得る。この透視像データはデータ記
憶装置DS1に一時記憶され、前処理装置PPCによ
り、断面像再構成のための前処理を施す。その前
処理の結果得られたデータはデータ記憶装置DS
2に一時記憶される。このデータ記憶装置DS2
に記憶されたデータは、オーバーラツプ処理装置
OVLPで処理され、データ記憶装置DS3に一時
記憶される。
このデータに対して、フーリエ変換装置FFT
との間でフーリエ変換の処理が施され、その結果
はデータ記憶装置DS4に一時記憶される。その
データ記憶装置DS4に記憶されたデータについ
て、フイルタ装置FILTによりフイルタ処理が施
され、更に、逆フーリエ変換装置IFFTにより、
逆フーリエ変換処理が施され、データ記憶装置
DS5に一時記憶される。このデータは逆投影制
御装置BPCにより処理されて映像データとなり、
映像データ記憶装置IMに記憶される。この映像
データは像写真撮影装置MFCで写真撮影され、
画像表示装置GDCに表示される。
との間でフーリエ変換の処理が施され、その結果
はデータ記憶装置DS4に一時記憶される。その
データ記憶装置DS4に記憶されたデータについ
て、フイルタ装置FILTによりフイルタ処理が施
され、更に、逆フーリエ変換装置IFFTにより、
逆フーリエ変換処理が施され、データ記憶装置
DS5に一時記憶される。このデータは逆投影制
御装置BPCにより処理されて映像データとなり、
映像データ記憶装置IMに記憶される。この映像
データは像写真撮影装置MFCで写真撮影され、
画像表示装置GDCに表示される。
上記テーブル・ガントリ−TGにはテーブル・
ガントリ−制御装置TGCが接続され、X線発生
部XGにはX線発生部制御装置XGCが接続され、
それぞれX線発生及び照射の制御を行うように構
成されている。
ガントリ−制御装置TGCが接続され、X線発生
部XGにはX線発生部制御装置XGCが接続され、
それぞれX線発生及び照射の制御を行うように構
成されている。
これらのテーブル・ガントリ−制御装置TGC
及びX線発生部制御装置XGC、更に、各データ
記憶装置DS1〜DS5、各処理装置、制御装置及
び変換装置は、撮影制御装置SCCにより制御され
るように構成されている。
及びX線発生部制御装置XGC、更に、各データ
記憶装置DS1〜DS5、各処理装置、制御装置及
び変換装置は、撮影制御装置SCCにより制御され
るように構成されている。
このように構成されたX線CTでは、被検体PA
をテーブル・ガントリ−TGに置き、被検体PA
の断面部分に関してそれぞれ複数N個のサンプリ
ングデータとして得られる透視像データを多数の
角度方向に収集してデータ記憶装置DS1に記憶
し、これに前処理を施した後に、オーバーラツプ
処理装置OVLPで次のような処理を実行する。
OVLPは、必要に応じてDS2に演算により生成
後格納された対向ビユーデータより式(1)のCv(I)
を計算し、DS2に格納されてあるDv(I)との間で
ウエイト処理2又は4の演算を行い、その結果
Rv(I)をDS3に格納する(Rv(I)は、式(2)又は(3)
又は(4)又は(5))。そのデータについて、フーリエ
変換装置FFTでフーリエ変換により上記断面部
分に関するフーリエ像データを得る。このデータ
はデータ記憶装置DS4に記憶され、フイルタ装
置FILTによるフイルタ処理及び逆フーリエ変換
装置IFFTによる逆フーリエ変換装処理が施され
て、放射線吸収係数に対応するデータとなり、逆
投影制御装置BPCにより逆投影されて、各部の
放射線吸収係数の分布をなす画像データとなる。
これは、画像表示装置GDCに表示すると共に、
像写真撮影装置MFCで写真撮影を行う。
をテーブル・ガントリ−TGに置き、被検体PA
の断面部分に関してそれぞれ複数N個のサンプリ
ングデータとして得られる透視像データを多数の
角度方向に収集してデータ記憶装置DS1に記憶
し、これに前処理を施した後に、オーバーラツプ
処理装置OVLPで次のような処理を実行する。
OVLPは、必要に応じてDS2に演算により生成
後格納された対向ビユーデータより式(1)のCv(I)
を計算し、DS2に格納されてあるDv(I)との間で
ウエイト処理2又は4の演算を行い、その結果
Rv(I)をDS3に格納する(Rv(I)は、式(2)又は(3)
又は(4)又は(5))。そのデータについて、フーリエ
変換装置FFTでフーリエ変換により上記断面部
分に関するフーリエ像データを得る。このデータ
はデータ記憶装置DS4に記憶され、フイルタ装
置FILTによるフイルタ処理及び逆フーリエ変換
装置IFFTによる逆フーリエ変換装処理が施され
て、放射線吸収係数に対応するデータとなり、逆
投影制御装置BPCにより逆投影されて、各部の
放射線吸収係数の分布をなす画像データとなる。
これは、画像表示装置GDCに表示すると共に、
像写真撮影装置MFCで写真撮影を行う。
尚、本発明は前述の実施例に限定されることは
なく、例えば、次のようにすることもできる。
なく、例えば、次のようにすることもできる。
(1) 本発明は、360゜以下のスキヤンデータをもと
に体動等のモーシヨン・アテイフアクト
(Motion Artifact)を軽減するCT装置の他
に、180゜+η(ηは放射線源より見た被検体領
域を含む最小のフアン・アングル)(Fan
Angle)のスキヤンデータより360゜方向の画像
再構成を行う画像再構成方式(ハーフスキヤン
(Half Scan)イメージ再構成やセグメント・
イメージ再構成等)等への適用も可能である。
に体動等のモーシヨン・アテイフアクト
(Motion Artifact)を軽減するCT装置の他
に、180゜+η(ηは放射線源より見た被検体領
域を含む最小のフアン・アングル)(Fan
Angle)のスキヤンデータより360゜方向の画像
再構成を行う画像再構成方式(ハーフスキヤン
(Half Scan)イメージ再構成やセグメント・
イメージ再構成等)等への適用も可能である。
(2) フアン・ビームCTのみならず平行ビーム
(Parallel Beam)を用いるCT装置にも実施で
きる。
(Parallel Beam)を用いるCT装置にも実施で
きる。
(3) CT装置への適用では、その構成等において
様々の変形が可能である。
様々の変形が可能である。
(4) アルゴリズムの変形
アルゴリズムは式(1)に限定されない。
例えば、下記のようにしてもよい。
Cv(I)
=(1/2){HvN(I)
+Hv+N(I)}
(N=1,2,……)
Cv(I)
=(1/4){Hv−2(I)
+Hv−1(I)+Hv+1(I)
+Hv+2(I)}
(発明の効果)
以上説明したように、本発明によれば次のよう
な効果がある。
な効果がある。
(1) 体動補正のためのオーバーラツプの方向で発
生する多数の直線的なストリーク状のアーテイ
フアクトがなくなる。
生する多数の直線的なストリーク状のアーテイ
フアクトがなくなる。
(2) 高コントラスト変化部での接線方向のアーテ
イフアクトがなくなる。
イフアクトがなくなる。
(3) 長い直線部を持つ被検体で、直線状のアーテ
イフアクトがなくなる。
イフアクトがなくなる。
(4) スムージングを全く加えることなく、上記(1)
〜(3)の実現が可能なので、生成されるイメージ
はCT画像再構成の原理に非常に忠実な、従つ
て自然なシメージとなる。
〜(3)の実現が可能なので、生成されるイメージ
はCT画像再構成の原理に非常に忠実な、従つ
て自然なシメージとなる。
以上のように、本質的には放射線源〜検出器系
のノイズ等に起因すると考えられるアーテイフア
クトを軽減・除去し、より再構成原理に忠実なイ
メージの生成を可能とすることができる。
のノイズ等に起因すると考えられるアーテイフア
クトを軽減・除去し、より再構成原理に忠実なイ
メージの生成を可能とすることができる。
第1図は本発明に係るX線CTの構成図、第2
図は従来の方式におけるアーテイフアクトの発生
例を示す図、第3図は本発明を説明するための
図、第4図は画像再構成処理における手順の例を
示すフローチヤートである。 PA…被検体、S…X線検出部、DAS…データ
収集装置、DS1,DS2,DS3,DS4,DS5…
データ記憶装置、PPC…前処理装置、OVLP…オ
ーバーラツプ処理装置、FFT…フーリエ変換装
置、FILT…フイルタ装置、IFFT…逆フーリエ
変換装置、BPC…逆投影制御装置、IM…映像デ
ータ記憶装置、MFC…像写真撮影装置、GDC…
画像表示装置、TGC…テーブル・ガントリー制
御装置、XGC…X線発生部制御装置、XG…X線
発生部、TG…テーブル・ガントリー、SCC…撮
影制御装置。
図は従来の方式におけるアーテイフアクトの発生
例を示す図、第3図は本発明を説明するための
図、第4図は画像再構成処理における手順の例を
示すフローチヤートである。 PA…被検体、S…X線検出部、DAS…データ
収集装置、DS1,DS2,DS3,DS4,DS5…
データ記憶装置、PPC…前処理装置、OVLP…オ
ーバーラツプ処理装置、FFT…フーリエ変換装
置、FILT…フイルタ装置、IFFT…逆フーリエ
変換装置、BPC…逆投影制御装置、IM…映像デ
ータ記憶装置、MFC…像写真撮影装置、GDC…
画像表示装置、TGC…テーブル・ガントリー制
御装置、XGC…X線発生部制御装置、XG…X線
発生部、TG…テーブル・ガントリー、SCC…撮
影制御装置。
Claims (1)
- 【特許請求の範囲】 1 スキヤンの開始点近傍のビユー範囲内及び終
了点近傍のビユー範囲内の各方向のビユーデータ
を得るに際して、当該方向に関する実測によるス
キヤン・データと、当該方向とは180゜±βずれた
範囲内の方向のビユーデータ群より算出した当該
方向に関する算出ビユーデータとを用い、体動に
よるアーテイフアクトを軽減化した各方向の補正
ビユーデータを生成し、画像再構成時に、この補
正ビユーデータを、前記スキヤンの開始点近傍の
ビユー範囲内及び終了点近傍のビユー範囲内の当
該方向のビユーデータとして使用するCT装置に
おいて、 第Vビユーの第Iチヤネルの実測による前記ス
キヤン・データをDv(I)、前記180゜±βずれた範
囲内の方向のビユーデータ群より算出した第Vビ
ユーの第Iチヤネルの前記算出ビユーデータを
Hv(I)で表わす時、第V−kビユーの第Iチヤネ
ルの算出ビユーデータHv−k(I)と第V+kビユ
ーの第Iチヤネルの算出ビユーデータHv+k(I)
ただし、k=1,2,……)との一組又は異なる
複数個のkに対応する複数組のビユーデータ群の
加算平均であるCv(I)と、前記第Vビユーの第
チヤネルの実測によるスキヤン・データDv(I)と
を重み付け加算したものを求める手段を設け、 該手段により求めたものを、第Vビユーの第I
チヤネルの前記補正ビユーデータとして用いるこ
とを特徴とするCT装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP15697684A JPS6137139A (ja) | 1984-07-27 | 1984-07-27 | Ct装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP15697684A JPS6137139A (ja) | 1984-07-27 | 1984-07-27 | Ct装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS6137139A JPS6137139A (ja) | 1986-02-22 |
| JPH0458336B2 true JPH0458336B2 (ja) | 1992-09-17 |
Family
ID=15639443
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP15697684A Granted JPS6137139A (ja) | 1984-07-27 | 1984-07-27 | Ct装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS6137139A (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2839486B2 (ja) * | 1986-12-12 | 1998-12-16 | 株式会社東芝 | 断層撮影装置 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE3132674C2 (de) * | 1981-08-19 | 1983-12-08 | Degussa Ag, 6000 Frankfurt | Verfahren zur Herstellung von Preßlingen |
| JPS5832746A (ja) * | 1981-08-19 | 1983-02-25 | 株式会社東芝 | 放射線断層撮影装置 |
| JPS58118735A (ja) * | 1982-01-06 | 1983-07-14 | 株式会社東芝 | コンピユ−タ断層撮影装置 |
-
1984
- 1984-07-27 JP JP15697684A patent/JPS6137139A/ja active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| JPS6137139A (ja) | 1986-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cho et al. | Cone-beam CT for radiotherapy applications | |
| US6470070B2 (en) | Image reconstruction using multiple X-ray projections | |
| US5663995A (en) | Systems and methods for reconstructing an image in a CT system performing a cone beam helical scan | |
| US7949088B2 (en) | X-ray CT apparatus and method for processing X-ray projection data | |
| US6292530B1 (en) | Method and apparatus for reconstructing image data acquired by a tomosynthesis x-ray imaging system | |
| JP3547455B2 (ja) | 転頭運動をするスライスのct画像の再構成 | |
| US7315605B2 (en) | Method and device for reconstructing a 3D image data set of a moving object | |
| JPH04332538A (ja) | 物体の画像を作成する方法及び装置 | |
| WO1995022115A1 (en) | Reconstruction of images from cone beam data | |
| JP2011031070A (ja) | 半撮影域のみをカバーする縮小サイズ検出器を利用するコンピュータ断層撮影システムに用いる装置及び方法 | |
| JP2001511385A (ja) | 非対称検出器によるハーフスキャンct復元 | |
| CN108876730B (zh) | 校正运动伪影的方法、装置及设备和存储介质 | |
| JP2001161678A (ja) | コーン・ビーム・マルチスライス式ct補正方法及び装置 | |
| JPH11253441A (ja) | スペクトル・ア―ティファクトを有する画像の補正方法 | |
| US6332013B1 (en) | Methods and apparatus for tilted helical reconstruction multislice CT | |
| JPH0767445B2 (ja) | X線ct装置 | |
| US8213565B2 (en) | Method for correcting truncated projection data | |
| US6061423A (en) | Fluoroscopy image reconstruction | |
| JP4676641B2 (ja) | マルチ・スライスct走査の螺旋再構成の方法及び装置 | |
| JP3006722B2 (ja) | コンピュータ断層撮影装置 | |
| JPH119582A (ja) | X線コンピュータ断層撮影装置 | |
| JPH0966051A (ja) | X線断層撮影装置 | |
| JPH0414008B2 (ja) | ||
| JPH0458336B2 (ja) | ||
| JPH09164133A (ja) | X線撮影装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EXPY | Cancellation because of completion of term |