JPH0661328B2 - 再生像の視野外の物体に対するctデータの補償方式 - Google Patents
再生像の視野外の物体に対するctデータの補償方式Info
- Publication number
- JPH0661328B2 JPH0661328B2 JP15010990A JP15010990A JPH0661328B2 JP H0661328 B2 JPH0661328 B2 JP H0661328B2 JP 15010990 A JP15010990 A JP 15010990A JP 15010990 A JP15010990 A JP 15010990A JP H0661328 B2 JPH0661328 B2 JP H0661328B2
- Authority
- JP
- Japan
- Prior art keywords
- processed data
- attenuation
- convolution integral
- attenuation values
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 28
- 230000002093 peripheral effect Effects 0.000 claims description 36
- 230000005540 biological transmission Effects 0.000 claims description 24
- 238000013016 damping Methods 0.000 claims description 5
- 238000002591 computed tomography Methods 0.000 description 17
- 230000010354 integration Effects 0.000 description 13
- 238000012545 processing Methods 0.000 description 11
- 238000005259 measurement Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 10
- 238000001514 detection method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000000007 visual effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 210000003484 anatomy Anatomy 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000000472 traumatic effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/432—Truncation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Medical Treatment And Welfare Office Work (AREA)
Description
【発明の詳細な説明】 発明の背景 この発明は計算機式断層写真法、特にX線減衰測定値か
ら医療用の像を発生するのに使われる計算機式断層写真
法(CT)スキャナに関する。
ら医療用の像を発生するのに使われる計算機式断層写真
法(CT)スキャナに関する。
第1図に示す様に、人間の解剖学的な部分の像を発生す
る為に使われるCTスキャナは、患者テーブル10を持ち、
これはガントリー12の開口11の中に配置することが出来
る。高度にコリメートされたX線の源13が、開口11の片
側で、ガントリー12の中に取付けられ、1つ又は更に多
くの検出器14が開口の反対側に取付けられる。X線源13
及び検出器14は、患者を走査する間、開口11に沿って回
転させ、少なくとも180゜の回転範囲にわたる異なる多
数の角度から見たX線減衰測定値を求める。
る為に使われるCTスキャナは、患者テーブル10を持ち、
これはガントリー12の開口11の中に配置することが出来
る。高度にコリメートされたX線の源13が、開口11の片
側で、ガントリー12の中に取付けられ、1つ又は更に多
くの検出器14が開口の反対側に取付けられる。X線源13
及び検出器14は、患者を走査する間、開口11に沿って回
転させ、少なくとも180゜の回転範囲にわたる異なる多
数の角度から見たX線減衰測定値を求める。
患者の完全な走査は、X線源13及び検出器14の別々の角
度の向きで求められる1組のX線減衰測定値で構成され
る。この各々の1組の測定値が、業界では「ビューP
(View)」と呼ばれ、各組の測定値の結果が透過輪郭で
ある。第2A図に示す様に、各々のビューの1組の測定値
は、X線源13及び検出器14を、矢印15で示す様に、収集
の視野(ビューのフィールド)を横切る様に同時に並進
させることによって求めることが出来る。源13及び検出
器14が並進する時、患者を通抜ける一連のX線減衰の測
定を行ない、その結果得られる1組のデータが1つの角
度の向きに於けける透過輪郭となる。その後、源13及び
検出器14の角度の向きを(例えば1゜)変更し別のビュ
ーを収集する。各々の透過輪郭を収集する別の構成が第
2B図に示されている。この構成では、X線源13が扇形ビ
ームを発生し、それが患者を通過して、検出器配列14に
入射する。この配列にある各々の検出部14が別々の減衰
信号を発生し、全ての検出器14からの信号を別々に収集
して、所定の角度の向きに対する透過輪郭を発生する。
前に述べた構成と同じく、その後、X線源13及び検出器
配列14を異る角度の向きに回転し、次の透過輪郭を収集
する。
度の向きで求められる1組のX線減衰測定値で構成され
る。この各々の1組の測定値が、業界では「ビューP
(View)」と呼ばれ、各組の測定値の結果が透過輪郭で
ある。第2A図に示す様に、各々のビューの1組の測定値
は、X線源13及び検出器14を、矢印15で示す様に、収集
の視野(ビューのフィールド)を横切る様に同時に並進
させることによって求めることが出来る。源13及び検出
器14が並進する時、患者を通抜ける一連のX線減衰の測
定を行ない、その結果得られる1組のデータが1つの角
度の向きに於けける透過輪郭となる。その後、源13及び
検出器14の角度の向きを(例えば1゜)変更し別のビュ
ーを収集する。各々の透過輪郭を収集する別の構成が第
2B図に示されている。この構成では、X線源13が扇形ビ
ームを発生し、それが患者を通過して、検出器配列14に
入射する。この配列にある各々の検出部14が別々の減衰
信号を発生し、全ての検出器14からの信号を別々に収集
して、所定の角度の向きに対する透過輪郭を発生する。
前に述べた構成と同じく、その後、X線源13及び検出器
配列14を異る角度の向きに回転し、次の透過輪郭を収集
する。
各々の透過輪郭用のデータを収集した時、信号をフィル
タにかけ、補正し、コンピュータ・メモリに記憶する為
にディジタル化する。こう云う工程が業界では包括的に
「予備処理」と呼ばれ、データが収集されている時に実
時間行なわれる。その後、収集された透過輪郭を使っ
て、再生視野に於ける各々の容積要素のX線減衰係数を
示す像を再生する。こう云う減衰係数が「CT数」と呼ば
れる整数に変換され、このCT数を使って、CRT表示装置
の対応する画素の輝度を制御する。こうして、患者を通
るスライス中の解剖学的な構造を示す像が発生される。
タにかけ、補正し、コンピュータ・メモリに記憶する為
にディジタル化する。こう云う工程が業界では包括的に
「予備処理」と呼ばれ、データが収集されている時に実
時間行なわれる。その後、収集された透過輪郭を使っ
て、再生視野に於ける各々の容積要素のX線減衰係数を
示す像を再生する。こう云う減衰係数が「CT数」と呼ば
れる整数に変換され、このCT数を使って、CRT表示装置
の対応する画素の輝度を制御する。こうして、患者を通
るスライス中の解剖学的な構造を示す像が発生される。
記憶されている透過輪郭から像を再生するには、かなり
の計算を必要とし、実時間で行なうことは出来ない。像
を再生する為に現在行なわれている方法は、業界ではフ
ィルタ補正逆投影方式と呼ばれ、この方式を使う時に必
要な計算時間は、1つには、各々のビュー又は透過輪郭
の間に収集される減衰データの量によって決定される。
特に、この方式に於けるフィルタ工程はフーリエ変換を
用いて実施され、この変換の計算時間は、透過輪郭の為
に収集されたデータの量の変化によって大幅に影響を受
けることがある。
の計算を必要とし、実時間で行なうことは出来ない。像
を再生する為に現在行なわれている方法は、業界ではフ
ィルタ補正逆投影方式と呼ばれ、この方式を使う時に必
要な計算時間は、1つには、各々のビュー又は透過輪郭
の間に収集される減衰データの量によって決定される。
特に、この方式に於けるフィルタ工程はフーリエ変換を
用いて実施され、この変換の計算時間は、透過輪郭の為
に収集されたデータの量の変化によって大幅に影響を受
けることがある。
第3図について説明すると、像を正しく再生するには、
各々のビューのX線減衰値が、開口11内にある全ての物
体を通ることが必要である。物体が収集された視野より
も大きければ、第3図の垂直向きのビューによって示す
様なある透過輪郭では、支持テーブル10により値が減衰
するが、第3図の水平向きのビューで示す様な他の透過
輪郭では、値が減衰しない。その結果、全ての透過輪郭
を逆投影して、再生された視野に於ける各々の容積要素
のCT数を決定する時、CT数が正確でなくなる。この精度
低下は表示像の中で背景の陰影として見られ、それは解
剖学的細部をぼかす程、輝度又は暗さを強めることがあ
る。
各々のビューのX線減衰値が、開口11内にある全ての物
体を通ることが必要である。物体が収集された視野より
も大きければ、第3図の垂直向きのビューによって示す
様なある透過輪郭では、支持テーブル10により値が減衰
するが、第3図の水平向きのビューで示す様な他の透過
輪郭では、値が減衰しない。その結果、全ての透過輪郭
を逆投影して、再生された視野に於ける各々の容積要素
のCT数を決定する時、CT数が正確でなくなる。この精度
低下は表示像の中で背景の陰影として見られ、それは解
剖学的細部をぼかす程、輝度又は暗さを強めることがあ
る。
勿論、この問題に対する解決策は、開口11内にある物体
全体が収集データの視野の中にある様に保証することで
ある。例えば、患者の頭を作像する時、米国特許第4,40
0,820号に記載されている様な頭保持体15を用い、テー
ブル10の端から伸ばし、頭の輪郭と密接に追述する様に
する。頭保持体15は患者を支持しているが、どの角度で
も視野内に包み込まれるように必要な視野の寸法を目立
って増加させるものではない。
全体が収集データの視野の中にある様に保証することで
ある。例えば、患者の頭を作像する時、米国特許第4,40
0,820号に記載されている様な頭保持体15を用い、テー
ブル10の端から伸ばし、頭の輪郭と密接に追述する様に
する。頭保持体15は患者を支持しているが、どの角度で
も視野内に包み込まれるように必要な視野の寸法を目立
って増加させるものではない。
都合の悪いことに、全ての物体を視野に局限することが
出来ない場合が多くある。例えば、外傷のある患者では
頭保持体を使うことが出来ないことがあり、この場合、
テーブル10を支持の為に使わなければならないし、それ
が開口11の中に入る。こう云う場合、その結果、像の品
質が劣化するのを我慢するか、再生像を作るのに使われ
る収集データの視野を拡大しなければらないが、そうす
れば計算時間が長くなる。
出来ない場合が多くある。例えば、外傷のある患者では
頭保持体を使うことが出来ないことがあり、この場合、
テーブル10を支持の為に使わなければならないし、それ
が開口11の中に入る。こう云う場合、その結果、像の品
質が劣化するのを我慢するか、再生像を作るのに使われ
る収集データの視野を拡大しなければらないが、そうす
れば計算時間が長くなる。
発明の要約 この発明は像を再生するのに使われる透過輪郭用のデー
タを補償する方法、更に具体的に云えば、再生用の視野
の外側にある影響を帳消しにする様に透過輪郭用のデー
タを変更する方法に関する。この方法は、再生用の視野
の外側にある物体を含む視野から透過輪郭用のデータを
収集し、再生用の視野の中にある中心の透過輪郭用のデ
ータを畳込み積分の核と畳込み積分することによって、
透過輪郭用のデータにフィルタ作用を加え、再生用視野
の外側にある周辺の透過輪郭用データを畳込み積分の核
と別々に畳込み積分し、フィルタにかけた周辺の透過輪
郭用データをフィルタにかけた中心の透過輪郭用のデー
タに加算することによって、フィルタにかけた透過輪郭
用データを補償し、補償済みのフィルタにかけた透過輪
郭用データを使って像を再生する工程を含む。
タを補償する方法、更に具体的に云えば、再生用の視野
の外側にある影響を帳消しにする様に透過輪郭用のデー
タを変更する方法に関する。この方法は、再生用の視野
の外側にある物体を含む視野から透過輪郭用のデータを
収集し、再生用の視野の中にある中心の透過輪郭用のデ
ータを畳込み積分の核と畳込み積分することによって、
透過輪郭用のデータにフィルタ作用を加え、再生用視野
の外側にある周辺の透過輪郭用データを畳込み積分の核
と別々に畳込み積分し、フィルタにかけた周辺の透過輪
郭用データをフィルタにかけた中心の透過輪郭用のデー
タに加算することによって、フィルタにかけた透過輪郭
用データを補償し、補償済みのフィルタにかけた透過輪
郭用データを使って像を再生する工程を含む。
この発明の全般的な目的は、像再生の計算時間を目立っ
て長くせずに、再生像のビューの外側にある物体に対す
るCTデータを補償することである。データ収集の視野
は、必要に応じて全ての物体を含む様に増大する。その
結果得られる透過輪郭は、像の再生に必要な減衰値を含
む中心領域と、再生像の視野の外側にある減衰値を含む
周辺領域とを含む。中心領域は、所望の畳込み積分の核
を用いて、普通の様なフィルタ作用にかけ、周辺領域
は、同じ畳込み積分の核を使うが、時間のかゝるフーリ
エ変換を使わないで済む簡略にした手順を用いて、別々
にフィルタ作用にかける。その後、フィルタにかけられ
た中心領域のデータは、フィルタにかけられた周辺領域
のデータを用いて補償し、再生像の視野の外側にある物
体の影響を帳消しにする。計算時間は長くなるが、その
増加は、普通の方法で透過輪郭全体にフィルタ作用をか
けるのに要する時間の増加の僅かな端数である。
て長くせずに、再生像のビューの外側にある物体に対す
るCTデータを補償することである。データ収集の視野
は、必要に応じて全ての物体を含む様に増大する。その
結果得られる透過輪郭は、像の再生に必要な減衰値を含
む中心領域と、再生像の視野の外側にある減衰値を含む
周辺領域とを含む。中心領域は、所望の畳込み積分の核
を用いて、普通の様なフィルタ作用にかけ、周辺領域
は、同じ畳込み積分の核を使うが、時間のかゝるフーリ
エ変換を使わないで済む簡略にした手順を用いて、別々
にフィルタ作用にかける。その後、フィルタにかけられ
た中心領域のデータは、フィルタにかけられた周辺領域
のデータを用いて補償し、再生像の視野の外側にある物
体の影響を帳消しにする。計算時間は長くなるが、その
増加は、普通の方法で透過輪郭全体にフィルタ作用をか
けるのに要する時間の増加の僅かな端数である。
この発明の上記並びにその他の目的及び利点は、以下の
説明から明らかになろう。
説明から明らかになろう。
次にこの発明の好ましい実施例を示した図面について説
明する。然し、この実施例は必ずしもこの発明の範囲全
体を表わすものではなく、この発明の範囲を解釈するに
あたって特許請求の範囲を参照されたい。
明する。然し、この実施例は必ずしもこの発明の範囲全
体を表わすものではなく、この発明の範囲を解釈するに
あたって特許請求の範囲を参照されたい。
好ましい実施例の説明 第5図について説明すると、CT装置の動作が、コンピュ
ータ・プロセッサ26及びディスク・メモリ27を含むプロ
グラム式データを処理装置25によって制御される。ディ
スク・メモリ27が、コンピュータ・プロセッサ26が患者
の走査及び像の再生と表示に使うプログラムを記憶して
いる。更にこれは、収集されたデータと再生像データと
をも短期的に記憶する。コンピュータ・プロセッサが、
図示の様に、装置の他の素子と接続するのに適した入力
及び出力ポートを持つ汎用ミニコンピュータを含む。更
にそれぞれが米国特第4,494,141号に記載される様な配
列プロセッサを含む。
ータ・プロセッサ26及びディスク・メモリ27を含むプロ
グラム式データを処理装置25によって制御される。ディ
スク・メモリ27が、コンピュータ・プロセッサ26が患者
の走査及び像の再生と表示に使うプログラムを記憶して
いる。更にこれは、収集されたデータと再生像データと
をも短期的に記憶する。コンピュータ・プロセッサが、
図示の様に、装置の他の素子と接続するのに適した入力
及び出力ポートを持つ汎用ミニコンピュータを含む。更
にそれぞれが米国特第4,494,141号に記載される様な配
列プロセッサを含む。
コンピュータ・プロセッサ26の出力ポートがX線制御回
路28に接続され、この回路がX線管13を制御する。X線
管13の高圧及びその陰極電流は、正しい線量が得られる
様に制御される。高圧及び陰極電流はオペレータによっ
て選ばれる。このオペレータが、所望の値をオペレータ
・コンソール30から入力し、コンピュータ・プロセッサ
26はその走査計画に従ってX線を発生する様に指示す
る。
路28に接続され、この回路がX線管13を制御する。X線
管13の高圧及びその陰極電流は、正しい線量が得られる
様に制御される。高圧及び陰極電流はオペレータによっ
て選ばれる。このオペレータが、所望の値をオペレータ
・コンソール30から入力し、コンピュータ・プロセッサ
26はその走査計画に従ってX線を発生する様に指示す
る。
前に述べた様にX線は扇形に分散し、ガントリーの開口
11の反対側に取付けられた検出器14の配列がそれを受取
る。個別に852個のセル又は検出素子があり、その各々
が、X線管13を出て、開口11内に配置された患者を通る
直線通路を通る1本のX線を検査する。検出器配線14
は、その両端の各々に、源13からの減衰していないX線
を受取る一群の基準セルを含む。各々の検出素子で得ら
れた電流がアナログ電気信号として集められ、データ収
集装置31にあるA/D変換器によってディジタル数に変
換される。信号は、1番目の検出器から始まり、852番
目の検出器で終わる様に、逐次的にディジタル化され
る。全ての検出器からのディジタル化された測定値が完
全なビューである。米国特許第4,112,303号及び同第4,1
15,596号にはガントリー構造が詳しく説明されており、
米国特許第4,707,607号には検出器配列14が詳しく説明
されており、データ収集装置は米国特許第4,583,240号
に記載されている。ディジタル化された信号がコンピュ
ータ・プロセッサ26に入力される。
11の反対側に取付けられた検出器14の配列がそれを受取
る。個別に852個のセル又は検出素子があり、その各々
が、X線管13を出て、開口11内に配置された患者を通る
直線通路を通る1本のX線を検査する。検出器配線14
は、その両端の各々に、源13からの減衰していないX線
を受取る一群の基準セルを含む。各々の検出素子で得ら
れた電流がアナログ電気信号として集められ、データ収
集装置31にあるA/D変換器によってディジタル数に変
換される。信号は、1番目の検出器から始まり、852番
目の検出器で終わる様に、逐次的にディジタル化され
る。全ての検出器からのディジタル化された測定値が完
全なビューである。米国特許第4,112,303号及び同第4,1
15,596号にはガントリー構造が詳しく説明されており、
米国特許第4,707,607号には検出器配列14が詳しく説明
されており、データ収集装置は米国特許第4,583,240号
に記載されている。ディジタル化された信号がコンピュ
ータ・プロセッサ26に入力される。
データ収集装置31からのディジタル化された減衰測定値
は、「暗」、不均一な検出セルの感度及び利得、並びに
走査全体にわたるX線ビーム強度の変動を補償する様
に、周知の形で予備処理される。その後、ビーム硬化の
補正とデータの対数数形式への変換が行われていて、各
々の測定値がX線ビームの減衰の線積分を表わす様にす
る。この予備処理は、走査を実施している時、実時間で
行なわれ、第6図で示す様に、各々のビューの減衰値32
が、2次元生データ配列33の1行に記憶される。破線34
で示す様に、各々の減衰データが、ある1つの角度で観
察した時の作像すべき物体の透過輪郭となる。
は、「暗」、不均一な検出セルの感度及び利得、並びに
走査全体にわたるX線ビーム強度の変動を補償する様
に、周知の形で予備処理される。その後、ビーム硬化の
補正とデータの対数数形式への変換が行われていて、各
々の測定値がX線ビームの減衰の線積分を表わす様にす
る。この予備処理は、走査を実施している時、実時間で
行なわれ、第6図で示す様に、各々のビューの減衰値32
が、2次元生データ配列33の1行に記憶される。破線34
で示す様に、各々の減衰データが、ある1つの角度で観
察した時の作像すべき物体の透過輪郭となる。
走査が完了した時、生データが配列33は各々の行に、1
つのビューからの透過輪郭34を記憶している。従って、
この配列33の一方の次元は、走査で収集されるビューの
数によって決定され、他方の次元は各々のビューで収集
される検出セル信号の数によって決定される。収集され
る検出セル信号の数が、収集される視野を決定し、この
発明では、この視野(ビューのフィールド)は再生像を
歪める惧れのあるガントリー開口11内の全ての物体を収
める位に大きくすべきある。好ましい実施例では、これ
は852個までの検出セル信号を含むことが出来るる。
つのビューからの透過輪郭34を記憶している。従って、
この配列33の一方の次元は、走査で収集されるビューの
数によって決定され、他方の次元は各々のビューで収集
される検出セル信号の数によって決定される。収集され
る検出セル信号の数が、収集される視野を決定し、この
発明では、この視野(ビューのフィールド)は再生像を
歪める惧れのあるガントリー開口11内の全ての物体を収
める位に大きくすべきある。好ましい実施例では、これ
は852個までの検出セル信号を含むことが出来るる。
生データ配列33にある減衰値から像を再生する前に、第
7図に示す様に、データは畳込み積分の核35のフィルタ
作用を受ける、又それによって修正される。周知の様
に、核35を用いた畳込み積分は、各々の透過輪郭が後で
逆投影されて像を再生する時、各点の周りのぼけがなく
なる様に、減衰データを修正する。畳込み積分の核関数
K(x)はx=0に対して対称的であり、x=0では1
の値を持ち、その両側の値は負である。この畳込み積分
工程は、生データ配列33にある各々の素子の周りの減衰
値と畳込み積分の核35の値との乗算と、その結果を加算
して第6図に示す処理済みデータ配列36を形成すること
ゝを要するから、非常計算量が大きい。この作用を実施
する為に普通使われる方法が、生データ配列33の各行を
フーリエ変換し、畳込み積分の核35をフーリエ変換し、
変換された2つの波形を乗算することである。その後、
逆フーリエ変換を実施して、処理済みデータ配列36を作
る。
7図に示す様に、データは畳込み積分の核35のフィルタ
作用を受ける、又それによって修正される。周知の様
に、核35を用いた畳込み積分は、各々の透過輪郭が後で
逆投影されて像を再生する時、各点の周りのぼけがなく
なる様に、減衰データを修正する。畳込み積分の核関数
K(x)はx=0に対して対称的であり、x=0では1
の値を持ち、その両側の値は負である。この畳込み積分
工程は、生データ配列33にある各々の素子の周りの減衰
値と畳込み積分の核35の値との乗算と、その結果を加算
して第6図に示す処理済みデータ配列36を形成すること
ゝを要するから、非常計算量が大きい。この作用を実施
する為に普通使われる方法が、生データ配列33の各行を
フーリエ変換し、畳込み積分の核35をフーリエ変換し、
変換された2つの波形を乗算することである。その後、
逆フーリエ変換を実施して、処理済みデータ配列36を作
る。
この畳込み積分工程を行なうのに要する処理量は、処理
される視野の中に含まれる検出セル信号の数の関数とし
て、線形に増加するのではない。
される視野の中に含まれる検出セル信号の数の関数とし
て、線形に増加するのではない。
そうではなく、フーリエ変換方法を用いて畳込み積分を
実施するのに要する時間は、Nを処理される検出セル信
号の数より大きいか又はそれに等しい最小の2のべき数
として、大体Nlog2Nに比例する。例えば、512個又は
それより少ない数の検出セル信号を処理する時、N=51
2である。他方、513個又は更に多くの検出セル信号を処
理する時、N=1024であり、畳込み分の処理時間は大体
2倍になる。云い換えれば、検出処理信号の数を増やし
て、2のべき数の境界と交差すると、処理時間は大体2
倍になる。
実施するのに要する時間は、Nを処理される検出セル信
号の数より大きいか又はそれに等しい最小の2のべき数
として、大体Nlog2Nに比例する。例えば、512個又は
それより少ない数の検出セル信号を処理する時、N=51
2である。他方、513個又は更に多くの検出セル信号を処
理する時、N=1024であり、畳込み分の処理時間は大体
2倍になる。云い換えれば、検出処理信号の数を増やし
て、2のべき数の境界と交差すると、処理時間は大体2
倍になる。
この発明は、畳込み積分の積分処理時間を大幅に長くす
る様な2のべき数の境界との交差をせずに、収集用の広
い視野を持つ透過輪郭用データを処理する方法である。
この方法は、畳込み積分過程を実施する為に、コンピュ
ータ・プロセッサ26によって実行行されるプログラムに
よって実施される。次に、この方法を第9図のフローチ
ャート及び第8A図乃至第8D図のグラフについて説明す
る。
る様な2のべき数の境界との交差をせずに、収集用の広
い視野を持つ透過輪郭用データを処理する方法である。
この方法は、畳込み積分過程を実施する為に、コンピュ
ータ・プロセッサ26によって実行行されるプログラムに
よって実施される。次に、この方法を第9図のフローチ
ャート及び第8A図乃至第8D図のグラフについて説明す
る。
畳込み積分は、畳込み積分の核35を用いて、生データ配
列33にある各行のデータ(即ち、透過輪郭)に対して実
施する。第8A図に一番よく示されているが、収集された
データが、618個までの別々の素子又は減衰値を有する
透過輪郭を作る。これらが512個の素子からなる中心領
域50と、夫々53個までの素子からなる2つの周辺領域5
0,52とに分割される。第9図のブロック53で示す様に、
畳込み積分過程の最初の工程は、透過輪郭の512個の素
子からなる中心領域50の普通の畳込み積分を実施するこ
とである。
列33にある各行のデータ(即ち、透過輪郭)に対して実
施する。第8A図に一番よく示されているが、収集された
データが、618個までの別々の素子又は減衰値を有する
透過輪郭を作る。これらが512個の素子からなる中心領
域50と、夫々53個までの素子からなる2つの周辺領域5
0,52とに分割される。第9図のブロック53で示す様に、
畳込み積分過程の最初の工程は、透過輪郭の512個の素
子からなる中心領域50の普通の畳込み積分を実施するこ
とである。
即ち、中心領域50に対してフーリエ変換を行ない、変換
された畳込み積分に結果を乗じ、その和に逆フーリエ変
換を行なう。512個の素子の結果が、処理済みデータ配
列36(第6図)の対応する行に記憶され、この畳込み積
分データ54の見取図が第8B図に示されている。このデー
タを使って像を再生すると、周辺領域51,52に入ってい
る減衰値を考慮しないことによって、陰影効果が生ず
る。過程の残りは、周辺領域51,52を考慮に入れる様
に、処理済みの中心領域のデータ54を補正することであ
る。勿論、618素子の透過輪郭全体の畳込み積分を工程5
3で普通の様に行うことが出来るが、前に述べた様に2
のべき数の境界と交差し、処理時間が2倍近くになる。
された畳込み積分に結果を乗じ、その和に逆フーリエ変
換を行なう。512個の素子の結果が、処理済みデータ配
列36(第6図)の対応する行に記憶され、この畳込み積
分データ54の見取図が第8B図に示されている。このデー
タを使って像を再生すると、周辺領域51,52に入ってい
る減衰値を考慮しないことによって、陰影効果が生ず
る。過程の残りは、周辺領域51,52を考慮に入れる様
に、処理済みの中心領域のデータ54を補正することであ
る。勿論、618素子の透過輪郭全体の畳込み積分を工程5
3で普通の様に行うことが出来るが、前に述べた様に2
のべき数の境界と交差し、処理時間が2倍近くになる。
第8A図に戻って、x1及びX2を周辺領域51の境界と
し、x3及びx4を周辺領域52の境界とする。周辺領域
51の減衰データが任意の点xに於ける透過輪郭全体の畳
込み積分に対して次に示す項の寄与を持つ。
し、x3及びx4を周辺領域52の境界とする。周辺領域
51の減衰データが任意の点xに於ける透過輪郭全体の畳
込み積分に対して次に示す項の寄与を持つ。
こゝでk(x)畳込み積分の核、P(x′)は透過輪
郭、x′は2つの関数の積を、畳込み積分過程で、加算
する点を表わすダミー変数である。
郭、x′は2つの関数の積を、畳込み積分過程で、加算
する点を表わすダミー変数である。
畳込み積分の核K(x−x′)を任意の点x0を中心とす
るテーラー級数に展開すると K(x−x′)=K(x−x0) +(x0−x′)K′(x−x0)+……(2) こゝでK′はKの1次微分である。
るテーラー級数に展開すると K(x−x′)=K(x−x0) +(x0−x′)K′(x−x0)+……(2) こゝでK′はKの1次微分である。
この展開式を式(1)に代入すると 式(3)の第2項はx0を次の様に選ぶことによって除
くことが出来る。
くことが出来る。
式(4)は、周辺領域51の塊の中心の場所を見付ける為
の式である。従って、第1次近似で、畳込み積分過程に
寄与する式(3)の唯一の項は次の通りである。
の式である。従って、第1次近似で、畳込み積分過程に
寄与する式(3)の唯一の項は次の通りである。
この項は、周辺領域51にある全ての減衰値の積分に、周
辺領域51の塊の中心に置いた畳込み積分の核を乗じたも
のである。従って、周辺領域51,52の各々53個の減衰値
の組全体は、この領域の塊の中心にある1個の値に置換
えることが出来る。このことが第8C図で夫々矢印55,56
によって示されている。これらの矢印は、夫々周辺領域
51,52にある全ての減衰値の和に等しい1個の減衰値を
表わす。
辺領域51の塊の中心に置いた畳込み積分の核を乗じたも
のである。従って、周辺領域51,52の各々53個の減衰値
の組全体は、この領域の塊の中心にある1個の値に置換
えることが出来る。このことが第8C図で夫々矢印55,56
によって示されている。これらの矢印は、夫々周辺領域
51,52にある全ての減衰値の和に等しい1個の減衰値を
表わす。
第9図について説明すると、畳込み積分プログラムは、
ブロック57で示す様に、その前にある全ての減衰値を加
算することによって、各々の周辺領域51,52の積分を計
算する。その後、ブロック58で示す様に、式(4)で示
す通りに、各々の周辺領域51,52の塊の中心を計算す
る。次に、ブロック59で示す様に、畳込み積分の核35に
領域51の積分値を乗じて、周辺領域51の畳込み積分の第
1次近似を作る。周辺領域52に対しても同じ作用が行な
われる。ブロック60で示す様に、この様にして畳込み積
分した周辺領域は、その中心を塊の中心(x0)に等し
い分だけその中心から遠ざける向きに割出すことによ
り、畳込み積分された中心領域50に対する位置を占め
る。このことが第8D図に図示されており、畳込み積分し
た周辺領域51の中心インパルス57は中心領域の左側の距
離x0の所にあり、畳込み積分した周辺領域52の中心イ
ンパルス58は中心領域の右側にある。
ブロック57で示す様に、その前にある全ての減衰値を加
算することによって、各々の周辺領域51,52の積分を計
算する。その後、ブロック58で示す様に、式(4)で示
す通りに、各々の周辺領域51,52の塊の中心を計算す
る。次に、ブロック59で示す様に、畳込み積分の核35に
領域51の積分値を乗じて、周辺領域51の畳込み積分の第
1次近似を作る。周辺領域52に対しても同じ作用が行な
われる。ブロック60で示す様に、この様にして畳込み積
分した周辺領域は、その中心を塊の中心(x0)に等し
い分だけその中心から遠ざける向きに割出すことによ
り、畳込み積分された中心領域50に対する位置を占め
る。このことが第8D図に図示されており、畳込み積分し
た周辺領域51の中心インパルス57は中心領域の左側の距
離x0の所にあり、畳込み積分した周辺領域52の中心イ
ンパルス58は中心領域の右側にある。
各々の畳込み積分した周辺領域51,52の内、中心領域50
と重なる部分が、処理済みデータ配列36に記憶されてい
る畳込み積分した中心領域データ54に加算される。この
作用がブロック61で行われ、第8D図には陰影線を施した
領域62,63によって図示されている。勿論、実際には、
畳込み積分した周辺領域データは1次元の配列に記憶さ
れており、プログラムは中心の値から量x0だけ割出
す。この後、この配列内に残っている値が、中心領域の
左側の境界から開始して、処理済みデータ配列36の夫々
の値に加算される。右側に対しても同じ過程を繰返し
て、両方の周辺領域にある物体の効果について像データ
を補償する。
と重なる部分が、処理済みデータ配列36に記憶されてい
る畳込み積分した中心領域データ54に加算される。この
作用がブロック61で行われ、第8D図には陰影線を施した
領域62,63によって図示されている。勿論、実際には、
畳込み積分した周辺領域データは1次元の配列に記憶さ
れており、プログラムは中心の値から量x0だけ割出
す。この後、この配列内に残っている値が、中心領域の
左側の境界から開始して、処理済みデータ配列36の夫々
の値に加算される。右側に対しても同じ過程を繰返し
て、両方の周辺領域にある物体の効果について像データ
を補償する。
判定ブロック65で示す様に、畳込み積分プログラムはル
ープ状に戻って、生データ配列33にある各々の透過輪郭
を処理し、こうして完全な処理済みデータ配列36を作
る。この後、処理済みデータ配列を使って、逆投影方式
を利用して、普通の様に像を再生する。
ープ状に戻って、生データ配列33にある各々の透過輪郭
を処理し、こうして完全な処理済みデータ配列36を作
る。この後、処理済みデータ配列を使って、逆投影方式
を利用して、普通の様に像を再生する。
この発明はCT装置により、再生像を歪みの惧れのある全
ての物体を含む広い視野にわたって、データを収集する
ことが出来る様にする。再生像の視野は一層小さいが、
一層広い収集した視野内にある物体の影響を考慮に入れ
る様に、像データが補償されている。この補償は、フー
リエ変換工程を用いた時の2のべき数の境界との交差に
よる処理時間の不釣合な延長を伴わず達成される。
ての物体を含む広い視野にわたって、データを収集する
ことが出来る様にする。再生像の視野は一層小さいが、
一層広い収集した視野内にある物体の影響を考慮に入れ
る様に、像データが補償されている。この補償は、フー
リエ変換工程を用いた時の2のべき数の境界との交差に
よる処理時間の不釣合な延長を伴わず達成される。
第1図はこの発明を用いることが出来るCT装置の斜視
図、第2A図及び第2B図は第1図のCT装置に使うことが出
来る2種類の走査方式を示す略図、第3図は第2B図の走
査方式の略図で、視野の外側にある物体によって生ずる
問題を例示する図、第4図は第1図のCT装置内での患者
の部分的な斜視図、第5図の第1図のCT装置のブロック
図、第6図は第1図のCT装置によって収集され且つ処理
されるデータを示すグラフ、第7図は第6図のデータの
処理に使われる畳込み積分の核を示すグラフ、第8A図乃
至第8D図は収集されたデータがこの発明によってどの様
に処理されるかを示すグラフ、第9図はこの発明を実施
する為に第1図のCT装置によって実行される畳込み積分
プログラムのフローチャートである。 [主な符号の説明] 34:透過輪郭、 35:畳込み積分の核、 50:中心領域、 51,52:周辺領域、 62,63:推定値。
図、第2A図及び第2B図は第1図のCT装置に使うことが出
来る2種類の走査方式を示す略図、第3図は第2B図の走
査方式の略図で、視野の外側にある物体によって生ずる
問題を例示する図、第4図は第1図のCT装置内での患者
の部分的な斜視図、第5図の第1図のCT装置のブロック
図、第6図は第1図のCT装置によって収集され且つ処理
されるデータを示すグラフ、第7図は第6図のデータの
処理に使われる畳込み積分の核を示すグラフ、第8A図乃
至第8D図は収集されたデータがこの発明によってどの様
に処理されるかを示すグラフ、第9図はこの発明を実施
する為に第1図のCT装置によって実行される畳込み積分
プログラムのフローチャートである。 [主な符号の説明] 34:透過輪郭、 35:畳込み積分の核、 50:中心領域、 51,52:周辺領域、 62,63:推定値。
Claims (6)
- 【請求項1】CT装置を用いて像を再生する方法に於て、
(a)1組の透過輪郭を収集し、各々の透過輪郭は、再
生しようとする像の視野に対応する中心領域に1組の減
衰値を持つと共に、前記中心領域の両側にある周辺領域
に配置された減衰値の組を持っており、(b)収集され
た1組の透過輪郭中の各々の透過輪郭の中心領域を畳込
み積分の核と畳込み積分して処理済みデータを発生し、
(c)各周辺領域にある減衰値の和に等しい1個の減衰
値を畳込積分の核に乗じることにより、周辺領域にある
減衰値の組の畳込み積分を推定し、(d)該推定値を前
記処理済みデータに加算することにより前記処理済みデ
ータを補償して、補償済み処理済みデータを発生し、
(e)該補償済み処理済みデータを逆投影することによ
って像を再生する各段階を含むことを特徴とする方法。 - 【請求項2】周辺領域にある各組の減衰値の畳込み積分
が、周辺領域の塊の中心にあって、該周辺領域にある全
ての減衰値の和に等しい振幅を持つ1個の減衰値に畳込
み積分の核を乗ずることによって推定される請求項1記
載の方法。 - 【請求項3】各々の周辺領域にある減衰値の推定された
畳込み積分が、推定された減衰値の配列であり、推定さ
れた減衰値の選ばれたものが、処理済みデータの要素に
加算される請求項2記載の方法。 - 【請求項4】CT装置を用いて像を再生する装置に於て、
(a)1組の透過輪郭を収集する手段であって、各々の
透過輪郭が、再生しようとする像の視野に対応する中心
領域に1組の減衰値を持つと共に、前記中心領域の両側
にある周辺領域に配置された減衰値の組を持っている手
段と、(b)収集された1組の透過輪郭中の各々の透過
輪郭の中心領域を畳込み積分の核と畳込み積分して処理
済みデータを発生する手段と、(c)各周辺領域にある
減衰値の和に等しい1個の減衰値を畳込み積分の核に乗
じることにより、周辺領域にある減衰値の畳込み積分を
推定する手段と、(d)該推定値を前記処理済みデータ
に加算することにより前記処理済みデータを補償して、
補償済み処理済みデータを発生する手段と、(e)該補
償済み処理済みデータを逆投影することによって像を再
生する手段とを含むことを特徴とする装置。 - 【請求項5】周辺領域にある各組の減衰値の畳込み積分
が、周辺領域の塊の中心にあって、該周辺領域にある全
ての減衰値の和に等しい振幅を持つ1個の減衰値に畳込
み積分の核を乗ずることによって推定される請求項4記
載の装置。 - 【請求項6】各々の周辺領域にある減衰値の推定された
畳込み積分が、推定された減衰値の配列であり、推定さ
れた減衰値の選ばれたものが、処理済みデータの要素に
加算される請求項5記載の装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US364,958 | 1982-04-02 | ||
US07/364,958 US5043890A (en) | 1989-06-12 | 1989-06-12 | Compensation of computed tomography data for objects positioned outside the field of view of the reconstructed image |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH0375045A JPH0375045A (ja) | 1991-03-29 |
JPH0661328B2 true JPH0661328B2 (ja) | 1994-08-17 |
Family
ID=23436870
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP15010990A Expired - Fee Related JPH0661328B2 (ja) | 1989-06-12 | 1990-06-11 | 再生像の視野外の物体に対するctデータの補償方式 |
Country Status (5)
Country | Link |
---|---|
US (1) | US5043890A (ja) |
EP (1) | EP0403209A1 (ja) |
JP (1) | JPH0661328B2 (ja) |
CA (1) | CA2010135A1 (ja) |
IL (1) | IL94627A (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100435732C (zh) * | 2004-04-29 | 2008-11-26 | Ge医疗***环球技术有限公司 | 图像重构方法及x射线ct装置 |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2575523B2 (ja) * | 1990-08-09 | 1997-01-29 | 三菱電機株式会社 | 実時間演算結果表示方式 |
US5611026A (en) * | 1992-12-21 | 1997-03-11 | General Electric Company | Combining a priori data with partial scan data to project three dimensional imaging of arbitrary objects with computerized tomography |
JP3414471B2 (ja) * | 1994-01-24 | 2003-06-09 | 株式会社東芝 | X線ct装置 |
JP4130055B2 (ja) * | 2000-08-31 | 2008-08-06 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 加算断層画像作成方法およびx線ct装置 |
US7254259B2 (en) * | 2002-10-04 | 2007-08-07 | Ge Medical Systems Global Technology Company Llc | Methods and apparatus for truncation compensation |
US6856666B2 (en) * | 2002-10-04 | 2005-02-15 | Ge Medical Systems Global Technology Company, Llc | Multi modality imaging methods and apparatus |
US6810102B2 (en) | 2002-10-04 | 2004-10-26 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for truncation compensation |
KR20030019278A (ko) * | 2002-11-21 | 2003-03-06 | 기태형 | 여과장치가 구비된 활어 수입/수출용 컨테이너 |
US20040120566A1 (en) * | 2002-12-19 | 2004-06-24 | Gines David L. | Compact storage of projection matrix for tomography using separable operators |
US7050528B2 (en) * | 2004-09-30 | 2006-05-23 | Wisconsin Alumni Research Foundation | Correction of CT images for truncated or incomplete projections |
US7782997B2 (en) | 2005-07-05 | 2010-08-24 | Koninklijke Philips Electronics N.V. | Reconstruction algorithm for object point outside the scan-field-of-view |
US7756315B2 (en) * | 2005-11-23 | 2010-07-13 | General Electric Company | Method and apparatus for field-of-view expansion of volumetric CT imaging |
US8213694B2 (en) * | 2008-02-29 | 2012-07-03 | Barco N.V. | Computed tomography reconstruction from truncated scans |
US10706595B2 (en) | 2017-03-13 | 2020-07-07 | General Electric Company | System and method for reconstructing an object via tomography |
CN108596993B (zh) * | 2018-02-26 | 2022-07-12 | 上海奕瑞光电子科技股份有限公司 | 校正图像不饱和伪影的***及校正方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5459856A (en) * | 1977-10-21 | 1979-05-14 | Hitachi Medical Corp | Xxray sectional image reproducer |
JPS5536329A (en) * | 1978-08-31 | 1980-03-13 | Oda Gosen Kogyo Kk | Composite crimped yarn |
US4305127A (en) * | 1979-11-29 | 1981-12-08 | Technicare Corporation | Projection completion method of partial area scan |
JPS58149739A (ja) * | 1982-02-28 | 1983-09-06 | 株式会社島津製作所 | コンピユ−タ断層撮影装置 |
US4550371A (en) * | 1982-09-27 | 1985-10-29 | General Electric Company | Method and apparatus for compensating CT images for truncated projections |
NL8204086A (nl) * | 1982-10-22 | 1984-05-16 | Philips Nv | Werkwijze en inrichting voor het bepalen van een stralingsverzwakkingsverdeling in een vlak van een lichaam. |
US4633398A (en) * | 1983-12-07 | 1986-12-30 | General Electric Company | Attenuation compensated emission reconstruction with simultaneous attenuation factor evaluation |
US4606004A (en) * | 1984-03-21 | 1986-08-12 | General Electric Company | Apparatus for reduction of filtration truncation errors |
JPH0652301B2 (ja) * | 1985-04-11 | 1994-07-06 | 株式会社東芝 | エミツシヨンct装置 |
NL8601004A (nl) * | 1986-04-21 | 1987-11-16 | Philips Nv | Computertomografie-inrichting met gekombineerd recursief- en convolutiefilter. |
CA1274922A (en) * | 1987-09-11 | 1990-10-02 | Terence Taylor | Region of interest tomography employing a differential scanning technique |
-
1989
- 1989-06-12 US US07/364,958 patent/US5043890A/en not_active Expired - Lifetime
-
1990
- 1990-02-15 CA CA002010135A patent/CA2010135A1/en not_active Abandoned
- 1990-06-06 IL IL94627A patent/IL94627A/xx not_active IP Right Cessation
- 1990-06-11 EP EP90306355A patent/EP0403209A1/en not_active Ceased
- 1990-06-11 JP JP15010990A patent/JPH0661328B2/ja not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100435732C (zh) * | 2004-04-29 | 2008-11-26 | Ge医疗***环球技术有限公司 | 图像重构方法及x射线ct装置 |
Also Published As
Publication number | Publication date |
---|---|
JPH0375045A (ja) | 1991-03-29 |
EP0403209A1 (en) | 1990-12-19 |
IL94627A (en) | 1993-05-13 |
CA2010135A1 (en) | 1990-12-12 |
US5043890A (en) | 1991-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2060181C (en) | Method and apparatus for computing tomographic scans | |
US5241576A (en) | Segmented detector containing sub-elements for separate measuring of a fan beam | |
US5293312A (en) | Method and apparatus for computing tomographic scans | |
JP4293303B2 (ja) | イメージング・システムを較正する方法及び装置 | |
US4149081A (en) | Removal of spectral artifacts and utilization of spectral effects in computerized tomography | |
US5225980A (en) | Reduction of image artifacts from support structures in tomographic imaging | |
US6366638B1 (en) | Methods and apparatus for CT scout image processing | |
JPH0661328B2 (ja) | 再生像の視野外の物体に対するctデータの補償方式 | |
US5606585A (en) | Methods and apparatus for multislice helical image reconstruction in a computer tomography system | |
JPH08509408A (ja) | 円すい状ビームデータからの画像の再構成 | |
US5400377A (en) | Artifact reduction method for tomographic image reconstruction using cross-plane rays | |
IL116258A (en) | Image Recovery Device and Spooler Scanning Method | |
JPH0790024B2 (ja) | Ct装置で画像を作成する方法 | |
Tofts et al. | Some sources of artefact in computed tomography | |
US6332013B1 (en) | Methods and apparatus for tilted helical reconstruction multislice CT | |
US5845003A (en) | Detector z-axis gain correction for a CT system | |
JPH0799540B2 (ja) | コンピユータによる写像再構成方法および装置 | |
JP3270153B2 (ja) | コンピュータトモグラフ | |
US6647084B1 (en) | Method and apparatus for filtering projection data of a helical scan | |
JP4582997B2 (ja) | 高速コンピュータ断層撮影方法 | |
EP0608237B1 (en) | Method and apparatus for computing tomographic and panoramic scans | |
JP3678375B2 (ja) | 放射線断層撮影装置 | |
JP4011201B2 (ja) | 測定信号処理方法および装置並びに放射線断層撮影装置 | |
JPH05212026A (ja) | たたみこみ‐逆投影法による医学撮像装置の作動方法 | |
Duerinckx et al. | Non-Linear Smoothing Filters And Noise Structure In Computed Tomography (CT) Scanning: A Preliminary Report |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
LAPS | Cancellation because of no payment of annual fees |