JPH0446578B2 - - Google Patents
Info
- Publication number
- JPH0446578B2 JPH0446578B2 JP17717783A JP17717783A JPH0446578B2 JP H0446578 B2 JPH0446578 B2 JP H0446578B2 JP 17717783 A JP17717783 A JP 17717783A JP 17717783 A JP17717783 A JP 17717783A JP H0446578 B2 JPH0446578 B2 JP H0446578B2
- Authority
- JP
- Japan
- Prior art keywords
- projection
- projections
- truncated
- moments
- determining
- 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
- 238000003384 imaging method Methods 0.000 claims description 23
- 238000000034 method Methods 0.000 claims description 15
- 230000005855 radiation Effects 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 230000000149 penetrating effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 208000031872 Body Remains Diseases 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
-
- 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
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Description
【発明の詳細な説明】
(産業上の利用分野)
この発明は計算機式断層写真法、更に具体的に
云えば、物体が走査器の撮影領域(視野)を越え
ることによつてCT像に生ずるアーテイフアクト
(擬似画像)を少なくする方法並びに装置に関す
る。
云えば、物体が走査器の撮影領域(視野)を越え
ることによつてCT像に生ずるアーテイフアクト
(擬似画像)を少なくする方法並びに装置に関す
る。
(従来技術)
CT走査器が正確な再生像を発生する為には、
検査する身体が走査器の撮影領域内に限定される
こと、並びに正確な基準チヤンネル・データを利
用し得ることが必要である。典型的には、減衰デ
ータを測定するチヤンネルに隣接している基準チ
ヤンネルは、入射放射ビーム強度の監視装置とし
て作用する。撮影領域は、身体を露出し、身体の
周りの複数個の角度から放射を検出する為に、透
過性放射源か、或いは透過性放射源と検出器がそ
の周りに回転する区域であり、典型的には円形で
ある。身体が撮影領域を越えて伸びると、測定さ
れる幾つかの投影(幾つかのビユー・データ)
は、(イ)基準検出器が遮蔽され、こうして幾つかの
ビユーに対する基準チヤンネル情報がくずれるこ
と、並びに(ロ)身体の内、撮影領域を越えて伸びる
部分に減衰情報を収集することが出来なくて、截
頭したビユー(身体の一部が撮影領域の外側に出
ているビユー)になる為に、誤りがある。前者の
場合、最終的な画像は不正な基準レベルのシフト
により、低い周波数の陰影を持つ。後者の場合、
アーテイフアクトが典型的には低い空間周波数の
陰影又はカツピングとなつて現われ、時には縞と
なつて現われる。
検査する身体が走査器の撮影領域内に限定される
こと、並びに正確な基準チヤンネル・データを利
用し得ることが必要である。典型的には、減衰デ
ータを測定するチヤンネルに隣接している基準チ
ヤンネルは、入射放射ビーム強度の監視装置とし
て作用する。撮影領域は、身体を露出し、身体の
周りの複数個の角度から放射を検出する為に、透
過性放射源か、或いは透過性放射源と検出器がそ
の周りに回転する区域であり、典型的には円形で
ある。身体が撮影領域を越えて伸びると、測定さ
れる幾つかの投影(幾つかのビユー・データ)
は、(イ)基準検出器が遮蔽され、こうして幾つかの
ビユーに対する基準チヤンネル情報がくずれるこ
と、並びに(ロ)身体の内、撮影領域を越えて伸びる
部分に減衰情報を収集することが出来なくて、截
頭したビユー(身体の一部が撮影領域の外側に出
ているビユー)になる為に、誤りがある。前者の
場合、最終的な画像は不正な基準レベルのシフト
により、低い周波数の陰影を持つ。後者の場合、
アーテイフアクトが典型的には低い空間周波数の
陰影又はカツピングとなつて現われ、時には縞と
なつて現われる。
或る種の走査では、身体が撮影領域の中にとゞ
まる様に保証しながら、所望のスライスを得るこ
とが簡単には出来ない。例えば、大きな人の胴体
の上側部分を走査するる時、肩が撮影領域の中に
とゞまる様に保証することが出来ないことがあ
る。
まる様に保証しながら、所望のスライスを得るこ
とが簡単には出来ない。例えば、大きな人の胴体
の上側部分を走査するる時、肩が撮影領域の中に
とゞまる様に保証することが出来ないことがあ
る。
(発明の目的)
以上に述べたことから、この発明の目的は、截
頭した投影の組から作られた画像のアーテイフア
クトを最小限に抑えることである。
頭した投影の組から作られた画像のアーテイフア
クトを最小限に抑えることである。
この目的を達成する為、この発明の目的は、不
良の基準情報を持つ投影を特定(識別)し、正規
化(基準情報を用いて投影データを相対比し、レ
ベル合せをする処理)する前に、この基準情報を
改善することである。
良の基準情報を持つ投影を特定(識別)し、正規
化(基準情報を用いて投影データを相対比し、レ
ベル合せをする処理)する前に、この基準情報を
改善することである。
この発明の別の目的は、截頭した投影データの
組を識別すると共に、截頭しなかつた組から取出
した最善の評価情報を用いてこう云う組を完全に
することである。最後に、この発明の目的は、
CT像を発生するのに必要な時間を許容し難い程
延ばさない様な能率の良い形で前述の目的を達成
することである。
組を識別すると共に、截頭しなかつた組から取出
した最善の評価情報を用いてこう云う組を完全に
することである。最後に、この発明の目的は、
CT像を発生するのに必要な時間を許容し難い程
延ばさない様な能率の良い形で前述の目的を達成
することである。
この他の目的並びに利点は、以下図面について
詳しく説明する所から明らかになろう。
詳しく説明する所から明らかになろう。
(実施例の記載)
この発明の或る好ましい実施例を説明するが、
この発明をこう云う実施例に制限するつもりはな
いことをことわつておきたい。むしろ以下の説明
の意図する所は、特許請求の範囲の記載によつて
定められたこの発明の範囲内に含まれる全ての変
更並びに均等物を包括することである。
この発明をこう云う実施例に制限するつもりはな
いことをことわつておきたい。むしろ以下の説明
の意図する所は、特許請求の範囲の記載によつて
定められたこの発明の範囲内に含まれる全ての変
更並びに均等物を包括することである。
第1図にCT走査器の主な素子が概略的に示さ
れている。例として、「回転専用」又は「第3世
代」走査器が示されているが、この発明の考えが
この様な形状に制限されないことに注意された
い。走査器が透過性放射源10を持つており、こ
れは回転陽極形X線管の形にする場合が非常に多
い。X線管10によつて発生された放射が11で
コリメートされ、薄い扇形の放射ビーム12とし
て発生し、これが患者用開口13を介してX線検
出器の配列14に向けて投射される。患者15の
様な検査しようとする身体を扇形X線ビーム12
の通路内で患者用開口13の中に位置ぎめし、身
体を通過するビームが、出会つた物体の密度に応
じて減衰する様にする。その結果、各々の検出セ
ル14a,14b等が、そのセルが受取つた放射
の強度に依存する電気信号を発生する。従つて、
こうして発生される信号は、X線ビームが通過し
た身体の部分によるX線ビームの減衰の測定値で
ある。
れている。例として、「回転専用」又は「第3世
代」走査器が示されているが、この発明の考えが
この様な形状に制限されないことに注意された
い。走査器が透過性放射源10を持つており、こ
れは回転陽極形X線管の形にする場合が非常に多
い。X線管10によつて発生された放射が11で
コリメートされ、薄い扇形の放射ビーム12とし
て発生し、これが患者用開口13を介してX線検
出器の配列14に向けて投射される。患者15の
様な検査しようとする身体を扇形X線ビーム12
の通路内で患者用開口13の中に位置ぎめし、身
体を通過するビームが、出会つた物体の密度に応
じて減衰する様にする。その結果、各々の検出セ
ル14a,14b等が、そのセルが受取つた放射
の強度に依存する電気信号を発生する。従つて、
こうして発生される信号は、X線ビームが通過し
た身体の部分によるX線ビームの減衰の測定値で
ある。
動作について説明すると、透過性放射源及び検
出器の配列を患者用開口の周りに回転する時、患
者に対する複数個の角度位置で、各々のセルから
のX線の指示値を求める。特定の角度に於ける
夫々一組の指示値は、投影データ又はビユーデー
タ(図データ)と呼ぶ場合が多い。利用し得る1
つの走査器に対する典型的なビユーは、512個の
個別の検出器の指示値で構成される。各々のビユ
ーに対する指示値をデイジタル化し、再構成計算
機16に送る。計算機16は、利用し得る多数の
アルゴリズムの内の1つを用いて、扇形ビームが
横切つた断面の画像を発生することが出来る。こ
の画像はCRT17に表示することが出来るし、
或いはその代わりに診断医が更に検討する為のフ
イルムを作るのに使うことが出来る。
出器の配列を患者用開口の周りに回転する時、患
者に対する複数個の角度位置で、各々のセルから
のX線の指示値を求める。特定の角度に於ける
夫々一組の指示値は、投影データ又はビユーデー
タ(図データ)と呼ぶ場合が多い。利用し得る1
つの走査器に対する典型的なビユーは、512個の
個別の検出器の指示値で構成される。各々のビユ
ーに対する指示値をデイジタル化し、再構成計算
機16に送る。計算機16は、利用し得る多数の
アルゴリズムの内の1つを用いて、扇形ビームが
横切つた断面の画像を発生することが出来る。こ
の画像はCRT17に表示することが出来るし、
或いはその代わりに診断医が更に検討する為のフ
イルムを作るのに使うことが出来る。
患者用開口13の中にあつて、透過性放射源及
び検出器の回転軸線を中心とする円20は、第1
図に示した形状に対する撮影領域を表わす。正確
な再構成を行なうとする場合、何よりも信頼性の
ある基準チヤンネル情報を得る為に、そして2番
目に截頭しない投影を発生する為に、患者15は
撮影領域20内に配置されていなければならな
い。
び検出器の回転軸線を中心とする円20は、第1
図に示した形状に対する撮影領域を表わす。正確
な再構成を行なうとする場合、何よりも信頼性の
ある基準チヤンネル情報を得る為に、そして2番
目に截頭しない投影を発生する為に、患者15は
撮影領域20内に配置されていなければならな
い。
図示の実施例では、検出器配列の末端の検出セ
ル14x,14yは患者の減衰を測定する為に使
わず、その代わりにX線源の減衰しない強度に特
有な基準の指示値を発生する為に使われる。従つ
て、患者用開口13を通過するが、撮影領域20
の外側を通過するX線21,22によつて表わさ
れる放射が、何の干渉も受けずに、セル14x,
14yに入射することが重要である。こう云う指
示値が全体を23で示した扇形の残りの部分から
の指示値を正規化する為に再構成計算機16で使
われる。
ル14x,14yは患者の減衰を測定する為に使
わず、その代わりにX線源の減衰しない強度に特
有な基準の指示値を発生する為に使われる。従つ
て、患者用開口13を通過するが、撮影領域20
の外側を通過するX線21,22によつて表わさ
れる放射が、何の干渉も受けずに、セル14x,
14yに入射することが重要である。こう云う指
示値が全体を23で示した扇形の残りの部分から
の指示値を正規化する為に再構成計算機16で使
われる。
第2図に問題を図式的に示してある。撮影領域
20が、撮影領域の中にある部分15aと撮影領
域を越える部分15b,15cを持つ被検体15
と共に示されている。1対の投影を25,26に
示してある。投影25は基準X線21,22が撮
影領域20を妨げなしに通過して、関連した検出
セルによつて受取られると云う点で、「正常」で
あることが判る。同様に、身体全体が検出される
放射を通過させる点で、身体15の全体が走査さ
れる。
20が、撮影領域の中にある部分15aと撮影領
域を越える部分15b,15cを持つ被検体15
と共に示されている。1対の投影を25,26に
示してある。投影25は基準X線21,22が撮
影領域20を妨げなしに通過して、関連した検出
セルによつて受取られると云う点で、「正常」で
あることが判る。同様に、身体全体が検出される
放射を通過させる点で、身体15の全体が走査さ
れる。
然し、ビユー26は2通りに截頭される投影を
示している。先づ、両方の基準ビーム21,22
が身体を通過するが、基準チヤンネル14x,1
4yの指示値は、このビユーに対するX線源の減
衰しない強度の真の指示値ではない。第2に、身
体全体の減衰がビユー26では測定されない。こ
れは身体の部分15b,15cが撮影領域の中に
なく、それらがX線に露出したとしても、こうい
うX線は検出されないからである。従つて、ビユ
ー26に対しては、走査器は身体の一部分15a
だけで構成された別の身体を測定していた場合の
様に作用する。
示している。先づ、両方の基準ビーム21,22
が身体を通過するが、基準チヤンネル14x,1
4yの指示値は、このビユーに対するX線源の減
衰しない強度の真の指示値ではない。第2に、身
体全体の減衰がビユー26では測定されない。こ
れは身体の部分15b,15cが撮影領域の中に
なく、それらがX線に露出したとしても、こうい
うX線は検出されないからである。従つて、ビユ
ー26に対しては、走査器は身体の一部分15a
だけで構成された別の身体を測定していた場合の
様に作用する。
上に述べた2つの状態は、CT像を普通の様に
単に処理した場合、CT像に重大な影響を持つこ
とがある。先づ、幾つかの投影で基準チヤンネル
がぼけることにより、この様にぼけた各々の投影
は一様な直流シフトを含むことになる。この直流
シフトはは、ぼけた投影を用いて得られた再構成
に低い周波数の陰影を発生する傾向がある。この
発明では、単に基準チヤンネルがぼけたことによ
つて生じた直流シフトは、典型的な臨床では
70H.U.程度になることが判つたが、これはかな
りの量である。
単に処理した場合、CT像に重大な影響を持つこ
とがある。先づ、幾つかの投影で基準チヤンネル
がぼけることにより、この様にぼけた各々の投影
は一様な直流シフトを含むことになる。この直流
シフトはは、ぼけた投影を用いて得られた再構成
に低い周波数の陰影を発生する傾向がある。この
発明では、単に基準チヤンネルがぼけたことによ
つて生じた直流シフトは、典型的な臨床では
70H.U.程度になることが判つたが、これはかな
りの量である。
この発明を実施する時、基準チヤンネルがぼけ
た截頭図に対する対数投影の直流シフトは、基準
チヤンネル情報がぼけた各々のビユーを識別し、
截頭しなかつた投影から取出した一層よい基準チ
ヤンネル情報でそれを置き換えることにより、実
質的に減少する。このやり方は、第3図について
説明するのが一番判り易い。第3図は第1図の装
置を用いて実施される一連の工程を示している。
た截頭図に対する対数投影の直流シフトは、基準
チヤンネル情報がぼけた各々のビユーを識別し、
截頭しなかつた投影から取出した一層よい基準チ
ヤンネル情報でそれを置き換えることにより、実
質的に減少する。このやり方は、第3図について
説明するのが一番判り易い。第3図は第1図の装
置を用いて実施される一連の工程を示している。
この手順を工程30で開始した後、最初に装置
は工程31で、透過性放射源及び検出器が患者の
周りを或る範囲のθにわたつて回転する時、各々
の角度θに対する投影データ又はビユー・データ
を発生するサイクルを進む。各々の投影に対して
約511個の指示値(投影データ)を発生すると共
に、透過性放射源及び検出器が身体の周りを回転
する時、約576回の投影をとることの出来る走査
器を利用することが出来る。この各々の投影には
第3図でIrと示した基準セル強度情報が関連して
いる。工程31で全ての投影を発生したら、装置
は工程32に進む。この工程32は、全ての基準
チヤンネル情報の内の最大強度Inaxを探して、
各々の投影データを調べる。この最大強度情報が
得られたら、装置は最大強度Inaxの或る端数であ
る閾値強度ITを決定する。次に各々の基準強度Ir
を閾値ITと比較して、何れかのビユーの基準チヤ
ンネルがぼけた惧れがあるかどうかの判定を下す
ことが出来る。工程34に示す様に、最初のビユ
ーを選択し、工程35の所で、そのビユーに対す
る基準チヤンネル情報Irを求める。工程36で、
基準チヤンネル情報を閾値強度ITと比較する。基
準チヤンネルが閾値より低いことが判れば、工程
37により、後で走査する為に、この投影にフラ
グをつける。これに反して、基準チヤンネル情報
が閾値より高ければ、このビユーはぼやけていな
い基準チヤンネルによつて撮影されたと結論す
る。その場合、工程38で強度情報Irを貯蔵し、
工程39で指数を増数して、後で使う為の典型的
な基準情報を発生する。工程40を使つて、次の
ビユーを選択し、この手順を繰返す。各々の場
合、閾値強度ITに対する関係に応じて、ビユーは
工程37でフラグをつけるか、或いは工程38で
和に加算される。全てのビユーを処理した後、工
程41で典型的な強度IAVGを計算する。図示の実
施例では、よいビユーに対する平均強度であり、
これは加算した強度の値をこうして加算した数値
の数で除すことによつて得られる。次にプログラ
ムは工程42で、工程37でフラグを付けた各々
のビユー・データ又は投影データを検索し、そう
いう各々のビユーに対し、工程42で、基準強度
情報Irを工程41で計算した平均情報IAVGに置き
換える。
は工程31で、透過性放射源及び検出器が患者の
周りを或る範囲のθにわたつて回転する時、各々
の角度θに対する投影データ又はビユー・データ
を発生するサイクルを進む。各々の投影に対して
約511個の指示値(投影データ)を発生すると共
に、透過性放射源及び検出器が身体の周りを回転
する時、約576回の投影をとることの出来る走査
器を利用することが出来る。この各々の投影には
第3図でIrと示した基準セル強度情報が関連して
いる。工程31で全ての投影を発生したら、装置
は工程32に進む。この工程32は、全ての基準
チヤンネル情報の内の最大強度Inaxを探して、
各々の投影データを調べる。この最大強度情報が
得られたら、装置は最大強度Inaxの或る端数であ
る閾値強度ITを決定する。次に各々の基準強度Ir
を閾値ITと比較して、何れかのビユーの基準チヤ
ンネルがぼけた惧れがあるかどうかの判定を下す
ことが出来る。工程34に示す様に、最初のビユ
ーを選択し、工程35の所で、そのビユーに対す
る基準チヤンネル情報Irを求める。工程36で、
基準チヤンネル情報を閾値強度ITと比較する。基
準チヤンネルが閾値より低いことが判れば、工程
37により、後で走査する為に、この投影にフラ
グをつける。これに反して、基準チヤンネル情報
が閾値より高ければ、このビユーはぼやけていな
い基準チヤンネルによつて撮影されたと結論す
る。その場合、工程38で強度情報Irを貯蔵し、
工程39で指数を増数して、後で使う為の典型的
な基準情報を発生する。工程40を使つて、次の
ビユーを選択し、この手順を繰返す。各々の場
合、閾値強度ITに対する関係に応じて、ビユーは
工程37でフラグをつけるか、或いは工程38で
和に加算される。全てのビユーを処理した後、工
程41で典型的な強度IAVGを計算する。図示の実
施例では、よいビユーに対する平均強度であり、
これは加算した強度の値をこうして加算した数値
の数で除すことによつて得られる。次にプログラ
ムは工程42で、工程37でフラグを付けた各々
のビユー・データ又は投影データを検索し、そう
いう各々のビユーに対し、工程42で、基準強度
情報Irを工程41で計算した平均情報IAVGに置き
換える。
全てのビユーに対してこういうことを行なう
と、最良の基準チヤンネル情報が利用可能にな
り、工程44を実施して、普通の様に各々のビユ
ーを正規化することが出来る。
と、最良の基準チヤンネル情報が利用可能にな
り、工程44を実施して、普通の様に各々のビユ
ーを正規化することが出来る。
基準チヤンネル情報をこの様に処理すると、約
70H.U.の誤差を50H.U.以上補正することになり、
従つて貴重な利点になることが判つた。残りの誤
差の大部分は、前に述べた2番目の形式の収差、
即ち、身体の一部分が一組の投影に寄与しない様
なビユーが存在することによるものである。問題
のこういう面は、第4図について説明するのが一
番判り易い。
70H.U.の誤差を50H.U.以上補正することになり、
従つて貴重な利点になることが判つた。残りの誤
差の大部分は、前に述べた2番目の形式の収差、
即ち、身体の一部分が一組の投影に寄与しない様
なビユーが存在することによるものである。問題
のこういう面は、第4図について説明するのが一
番判り易い。
第4図は或るビユー又は投影に対して測定され
た対数減衰値のグラフである。縦軸に振幅をと
り、扇形の中での検出器の位置を横軸に沿つてφ
として示す。中心の検出器は配列の中心のφ=0
に示す。このグラフは、検出器自体がφ1からφ2
までの扇形の角度だけを受入れるという点で、截
頭したビユーを示している。Pn(φ)は検出器の
指示値によつて得られた値を表わす。一組の投影
がφ1及びφ2を越えて夫々φ1′及びφ2′まで拡がつて
いることが示されている。関数g1(φ)及びg2
(φ)は截頭した一組の投影を示す。第2図の身
体の部分15b,15cが撮影領域より突出して
いることを前に述べたが、一組の投影のg1(φ)
及びg2(φ)の部分が失われていることが理解さ
れよう。従つて、問題の投影を収集する際に得ら
れた全ての情報がPn(φ)で表わされており、こ
れを普通の様にたゝみ込めば(重畳積分処理すれ
ば)、g1(φ)及びg2(φ)情報が無視されている
ことにより、この結果得られる画像にアーテイフ
アクトが生ずる。
た対数減衰値のグラフである。縦軸に振幅をと
り、扇形の中での検出器の位置を横軸に沿つてφ
として示す。中心の検出器は配列の中心のφ=0
に示す。このグラフは、検出器自体がφ1からφ2
までの扇形の角度だけを受入れるという点で、截
頭したビユーを示している。Pn(φ)は検出器の
指示値によつて得られた値を表わす。一組の投影
がφ1及びφ2を越えて夫々φ1′及びφ2′まで拡がつて
いることが示されている。関数g1(φ)及びg2
(φ)は截頭した一組の投影を示す。第2図の身
体の部分15b,15cが撮影領域より突出して
いることを前に述べたが、一組の投影のg1(φ)
及びg2(φ)の部分が失われていることが理解さ
れよう。従つて、問題の投影を収集する際に得ら
れた全ての情報がPn(φ)で表わされており、こ
れを普通の様にたゝみ込めば(重畳積分処理すれ
ば)、g1(φ)及びg2(φ)情報が無視されている
ことにより、この結果得られる画像にアーテイフ
アクトが生ずる。
第4図に示す様に投影が截頭になつている時、
即ち、Pn(φ)情報しか利用することが出来ない
時、物体の一意的な再構成は不可能であることが
知られている。然し、妥当な数の投影(例えば少
なくとも半分の撮影)が截頭されていなければ、
「完備関数」(complection functions)の低い空
間周波数の特徴の妥当な評価をし、更に処理する
為にこの最善の評価の完備投影を使うことが可能
である。
即ち、Pn(φ)情報しか利用することが出来ない
時、物体の一意的な再構成は不可能であることが
知られている。然し、妥当な数の投影(例えば少
なくとも半分の撮影)が截頭されていなければ、
「完備関数」(complection functions)の低い空
間周波数の特徴の妥当な評価をし、更に処理する
為にこの最善の評価の完備投影を使うことが可能
である。
この発明を実施する時、截頭したものも、截頭
しなかつたものも、投影の複数個のモーメント、
実施例では0次及び1次モーメントを使つて、截
頭の程度を決定すると共に、各々のビユーに対
し、モーメントの他に、問題のビユー内にある情
報だけを使つて、截頭情報を補正関数に割振る。
しなかつたものも、投影の複数個のモーメント、
実施例では0次及び1次モーメントを使つて、截
頭の程度を決定すると共に、各々のビユーに対
し、モーメントの他に、問題のビユー内にある情
報だけを使つて、截頭情報を補正関数に割振る。
第6A図で、装置が70で動作を開始し、最初
に工程71で、全ての角度θを含む一組の投影を
前に述べた様に発生することが判る。工程72を
使つて、ぼやけた基準チヤンネルの影響を取去る
為に、第3図の手順全体を行なう。次に工程73
を実施して、この一組内にある各々の投影に対
し、0次及び1次モーメントQ0及びQ1を決定す
る。次に第5図について、一組の投影のモーメン
トの性格を簡単に説明する。モーメントを計算す
る式は次の通りである。
に工程71で、全ての角度θを含む一組の投影を
前に述べた様に発生することが判る。工程72を
使つて、ぼやけた基準チヤンネルの影響を取去る
為に、第3図の手順全体を行なう。次に工程73
を実施して、この一組内にある各々の投影に対
し、0次及び1次モーメントQ0及びQ1を決定す
る。次に第5図について、一組の投影のモーメン
トの性格を簡単に説明する。モーメントを計算す
る式は次の通りである。
Qk=∫〓2′〓1′P(φ)φkdφ,k=0,1 (1)
こゝでQはモーメントであり、kはその次数を
示す指数であり、残りの変数は前に説明した通り
である。従つて、0次モーメントは、投影P(φ)
の振幅の積分と云う性質を持ち、投影内にある合
計「質量」に関するものとみなすことが出来る。
この質量は、物体全体が撮影領域の中にとゞまつ
ていれば、あまり変化しないから、0次モーメン
トは全てのビユーに対してかなり一定のまゝに
とゞまる筈である。ビームが平行な場合、0次モ
ーメントは実際に定数であるが、第1図に示す様
にビームが扇形である場合、0次モーメント55
は正弦状の関数になることがあり、一組全体にわ
たつて1サイクル以上にはならず、振幅は比較的
小さい。1次モーメントを56に示してあり、こ
れは0次モーメントによつて表わされる質量のビ
ユー心の場所に関するものと考えることが出来
る。0次モーメントと同じく、これは角度と共に
変化するが、一組の投影の360゜あたり、完全な2
周期しか持たない。
示す指数であり、残りの変数は前に説明した通り
である。従つて、0次モーメントは、投影P(φ)
の振幅の積分と云う性質を持ち、投影内にある合
計「質量」に関するものとみなすことが出来る。
この質量は、物体全体が撮影領域の中にとゞまつ
ていれば、あまり変化しないから、0次モーメン
トは全てのビユーに対してかなり一定のまゝに
とゞまる筈である。ビームが平行な場合、0次モ
ーメントは実際に定数であるが、第1図に示す様
にビームが扇形である場合、0次モーメント55
は正弦状の関数になることがあり、一組全体にわ
たつて1サイクル以上にはならず、振幅は比較的
小さい。1次モーメントを56に示してあり、こ
れは0次モーメントによつて表わされる質量のビ
ユー心の場所に関するものと考えることが出来
る。0次モーメントと同じく、これは角度と共に
変化するが、一組の投影の360゜あたり、完全な2
周期しか持たない。
この発明では、上に述べたモーメント情報を調
べて、予想される挙動からの変化を検出し、この
様な変化を生じたビユーを截頭したものと認識
し、完備したビユーで考えられる更によいモーメ
ントを評価し、その情報を利用して截頭した投影
を完成する。第5図はモーメントの予想しない様
な挙動が起る2つの区域を示してある。これらの
区域は57及び58の括弧で示してある。区域5
7では、0次モーメントQ0の振幅が予想しない
形で増加し、最初は挙動のよかつた1次モーメン
トQ1が予想しない起伏を持つことが判る。この
為、区域57内に入る々のビユーは截頭したもの
と云うことが出来、ビユーを完備したものにする
為に、実際のモーメントと比較すべき予想モーメ
ントを決定すべきである。この為、0次モーメン
トでは、多少とも真直ぐな線60を、完備したデ
ータの組のモーメントがそうなると思われる値の
よりよい評価として決定し、真直ぐな線分61を
一組の投影が完全であれば、その区域に対して考
えられる更によい1次モーメントになると評価す
る。同様な判断基準により、58の括弧に囲まれ
たビユーは截頭したものと認識され、真直ぐな線
62が0次モーメントになり得るものと評価さ
れ、真直ぐな線63が、ビユーが完備していれ
ば、この区域の1次モーメントとして考えられる
ものと評価される。真直ぐな線を使うのが便利で
あるが、希望によつては他の曲線を用いてもよ
い。
べて、予想される挙動からの変化を検出し、この
様な変化を生じたビユーを截頭したものと認識
し、完備したビユーで考えられる更によいモーメ
ントを評価し、その情報を利用して截頭した投影
を完成する。第5図はモーメントの予想しない様
な挙動が起る2つの区域を示してある。これらの
区域は57及び58の括弧で示してある。区域5
7では、0次モーメントQ0の振幅が予想しない
形で増加し、最初は挙動のよかつた1次モーメン
トQ1が予想しない起伏を持つことが判る。この
為、区域57内に入る々のビユーは截頭したもの
と云うことが出来、ビユーを完備したものにする
為に、実際のモーメントと比較すべき予想モーメ
ントを決定すべきである。この為、0次モーメン
トでは、多少とも真直ぐな線60を、完備したデ
ータの組のモーメントがそうなると思われる値の
よりよい評価として決定し、真直ぐな線分61を
一組の投影が完全であれば、その区域に対して考
えられる更によい1次モーメントになると評価す
る。同様な判断基準により、58の括弧に囲まれ
たビユーは截頭したものと認識され、真直ぐな線
62が0次モーメントになり得るものと評価さ
れ、真直ぐな線63が、ビユーが完備していれ
ば、この区域の1次モーメントとして考えられる
ものと評価される。真直ぐな線を使うのが便利で
あるが、希望によつては他の曲線を用いてもよ
い。
第6図に戻つて、工程74は区域57及び58
の同定を行ない、この様な括弧の区内にある各々
のビユーにフラグを付ける。次の工程75が、一
組の投影が前に述べた様に完備していた場合、モ
ーメントがどの様になつたかを決定する。この様
な決定をすると、装置はフラグを付した各々のビ
ユーに戻り、前に述べたモーメント情報をそのビ
ユーの情報と共に用いて、延長関数(第4図のg
(φ))を決定することにより、このビユーを完成
する。これは、利用し得る情報に低次多項式をは
めることによつて行なわれる。gi(φ)を表わす
のに選んだ多項式の次数が、表わさなければなら
ない制約の数を決定する。次にこの数が表示で計
算しなければならないモーメントの数を決定す
る。こゝで説明した特定の実施例は1次多項式
(即ち直線)を用いてお、この為にはQ0及びQ1だ
けを考慮すればよい。
の同定を行ない、この様な括弧の区内にある各々
のビユーにフラグを付ける。次の工程75が、一
組の投影が前に述べた様に完備していた場合、モ
ーメントがどの様になつたかを決定する。この様
な決定をすると、装置はフラグを付した各々のビ
ユーに戻り、前に述べたモーメント情報をそのビ
ユーの情報と共に用いて、延長関数(第4図のg
(φ))を決定することにより、このビユーを完成
する。これは、利用し得る情報に低次多項式をは
めることによつて行なわれる。gi(φ)を表わす
のに選んだ多項式の次数が、表わさなければなら
ない制約の数を決定する。次にこの数が表示で計
算しなければならないモーメントの数を決定す
る。こゝで説明した特定の実施例は1次多項式
(即ち直線)を用いてお、この為にはQ0及びQ1だ
けを考慮すればよい。
特定の低次多項式の係数を決定する為、一組の
延長投影に第1の制約が加えられ、その振幅が或
る点でゼロになることを必要とする。これは身体
が或る点で終るという肉体的な条件に類似する。
第4図は延長した投影の組がゼロになる点φ1′及
びφ2′を示している。φ1′及びφ2′の実際の位置も
やはり決定しなければならない。
延長投影に第1の制約が加えられ、その振幅が或
る点でゼロになることを必要とする。これは身体
が或る点で終るという肉体的な条件に類似する。
第4図は延長した投影の組がゼロになる点φ1′及
びφ2′を示している。φ1′及びφ2′の実際の位置も
やはり決定しなければならない。
0次及び1次モーメントを使う関数には別の2
つの制約が加えられる。上に述べたモーメントに
対して最善の評価値を決定したら、截頭した投影
の組に対する実際のモーメントが判つていれば、
予想のモーメントと同形になる様に測定したモー
メントを変えることが条件である。
つの制約が加えられる。上に述べたモーメントに
対して最善の評価値を決定したら、截頭した投影
の組に対する実際のモーメントが判つていれば、
予想のモーメントと同形になる様に測定したモー
メントを変えることが条件である。
有用であることが判つた最後の制約は、截頭し
た投影の組と延長した投影の組の間の境界にわた
つて、投影並びにその微係数の連続性があること
である。これはφに対する或る次数のgiの微係数
がφ1及びφ2で、φに対するPiの微係数と等しいと
云つてもよい。この制約に要求される微係数の次
数は、多項式の次数によつて決定される。1次多
項式、即ち直線を選んだ場合、完備した投影が境
界で連続的でであることを保証しさえすればよ
い。単に截頭した組の最後の点の振幅P1又はP2
を、夫々延長した組g1又はg2の最初の点に割当て
ることによつて、連続性を達成することが出来
る。
た投影の組と延長した投影の組の間の境界にわた
つて、投影並びにその微係数の連続性があること
である。これはφに対する或る次数のgiの微係数
がφ1及びφ2で、φに対するPiの微係数と等しいと
云つてもよい。この制約に要求される微係数の次
数は、多項式の次数によつて決定される。1次多
項式、即ち直線を選んだ場合、完備した投影が境
界で連続的でであることを保証しさえすればよ
い。単に截頭した組の最後の点の振幅P1又はP2
を、夫々延長した組g1又はg2の最初の点に割当て
ることによつて、連続性を達成することが出来
る。
第6A図及び第6b図で、工程76,77が第
1の截頭図(第4図のPn(φ)に対応する)を得
る様に作用する。次に工程78が前に説明した様
に行なわれ、延長した投影の組の最初の点の値を
この組の各々の端に於ける截頭しなかつた投影の
組の最後の点の値に設定する。この動作が行なわ
れる点が第4図の79,80に示してある。
1の截頭図(第4図のPn(φ)に対応する)を得
る様に作用する。次に工程78が前に説明した様
に行なわれ、延長した投影の組の最初の点の値を
この組の各々の端に於ける截頭しなかつた投影の
組の最後の点の値に設定する。この動作が行なわ
れる点が第4図の79,80に示してある。
次に装置は工程82に進み、この工程で0次モ
ーメントの情報を使うことによつて、截頭された
質量の量を決定する。更に具体的に云うと、第5
図で述べたことを念頭において、装置が括弧を付
した部分57内のビユーに対して作用しているこ
とを仮定すれば、装置は評価した0次モーメント
60と截頭データからの対応する測定モーメント
との間の差を決定する。
ーメントの情報を使うことによつて、截頭された
質量の量を決定する。更に具体的に云うと、第5
図で述べたことを念頭において、装置が括弧を付
した部分57内のビユーに対して作用しているこ
とを仮定すれば、装置は評価した0次モーメント
60と截頭データからの対応する測定モーメント
との間の差を決定する。
次に工程83を実施して、工程75からの1次
モーメントの情報を用いて、延長した投影の組に
於ける質量の合計量のビユー心を決定する。この
ことの物理的な意味は第4図から理解されよう。
点80,81が判つているから、工程82が曲線
g1(φ)及びg2(φ)の下にある合計面積を決定す
る。こういう関数を更に限定する為、工程83が
投影の組全体の質量の中心が何処にあるべきかを
決定し、この情報から、曲線g1(φ)及びg2(φ)
の下にある面積を割振ることが出来る。これは工
程85で、延長した投影の組に低次多項式をはめ
合せることによつて行なわれ、こうして図示の実
施例では、g1及びg2を直線であるとして、その勾
配及び交点を決定することにより、g1(φ)及び
g2(φ)を限定する。次に工程86が次の截頭し
た投影を選択し、この手順全体を繰返す。最終結
果として、再構成アルゴリズムの少なくとも低周
波数の連続性の条件を一番よく充たすと思われる
評価データで、各々の截頭した投影が完成され
る。
モーメントの情報を用いて、延長した投影の組に
於ける質量の合計量のビユー心を決定する。この
ことの物理的な意味は第4図から理解されよう。
点80,81が判つているから、工程82が曲線
g1(φ)及びg2(φ)の下にある合計面積を決定す
る。こういう関数を更に限定する為、工程83が
投影の組全体の質量の中心が何処にあるべきかを
決定し、この情報から、曲線g1(φ)及びg2(φ)
の下にある面積を割振ることが出来る。これは工
程85で、延長した投影の組に低次多項式をはめ
合せることによつて行なわれ、こうして図示の実
施例では、g1及びg2を直線であるとして、その勾
配及び交点を決定することにより、g1(φ)及び
g2(φ)を限定する。次に工程86が次の截頭し
た投影を選択し、この手順全体を繰返す。最終結
果として、再構成アルゴリズムの少なくとも低周
波数の連続性の条件を一番よく充たすと思われる
評価データで、各々の截頭した投影が完成され
る。
第6A図及び第6B図について説明した操作工
程は数学的に説明することも出来る。次に、この
発明の方法を実施する様に、再構成用計算機をプ
ログラムする助けとして、この数学的な説明をす
る。最初に、延長関数g1(φ)又はg2(φ)、又は
iを1又は2として一般的にgi(φ)に対して定
めるべき直線は、次の様に定義することが出来
る。
程は数学的に説明することも出来る。次に、この
発明の方法を実施する様に、再構成用計算機をプ
ログラムする助けとして、この数学的な説明をす
る。最初に、延長関数g1(φ)又はg2(φ)、又は
iを1又は2として一般的にgi(φ)に対して定
めるべき直線は、次の様に定義することが出来
る。
gi(φ)=aiφ+bi (2)
関数を限定する為に、係数ai及びbiを決定しな
ければならない。投影の連続性から gi(φi)=Pn(φi)≡Pi (3) これは数学的には、截頭曲線が第4図の点80,
81で延長曲線と滑らかに接続することを述べて
いる。この条件から、振幅がゼロの時のφiの値が
次の式で定められることが判る。
ければならない。投影の連続性から gi(φi)=Pn(φi)≡Pi (3) これは数学的には、截頭曲線が第4図の点80,
81で延長曲線と滑らかに接続することを述べて
いる。この条件から、振幅がゼロの時のφiの値が
次の式で定められることが判る。
φi′=−bi/ai (4)
これを念頭におけば、式(3)は次の様に書き直す
ことが出来る。
ことが出来る。
gi(φ)=ai(φ−φi)+Pi (5)
更に、延長関数(φi′−φi)の水平座標に関係
する関数δiを定義することが出来る。これは式(4)
に述べた交点及び勾配に等しいとおくことが出来
る。従つて δi=φi′−φi=−Pi/ai (6) 第6図の工程82及び83は数学的に次の式に
よつて表わすことが出来る。
する関数δiを定義することが出来る。これは式(4)
に述べた交点及び勾配に等しいとおくことが出来
る。従つて δi=φi′−φi=−Pi/ai (6) 第6図の工程82及び83は数学的に次の式に
よつて表わすことが出来る。
Rk≡Qk−∫〓2〓1Pn(φ)φkdφ
=∫〓1〓1g2(φ)φkdφ+∫〓2′〓2g2(φ)φkd
φ(7) この式は全ての次数のモーメントに一般的に適
用されるが、指数kを0及び1に等しいとおくこ
とにより、0次及び1次モーメントを用いる。量
Rkは截頭した組のモーメントとその組が完全で
あつた場合の予想モーメントとの間の差である。
量Qkがこの後者を表わし、φ1及びφ2の限界の間
のPnの積分が截頭データの組からの測定モーメ
ントを表わす。この差は、延長関数に対するg1
(φ)及びg2(φ)の積分について述べた様に、延
長関数の間で分割しなければならない量である。
φ(7) この式は全ての次数のモーメントに一般的に適
用されるが、指数kを0及び1に等しいとおくこ
とにより、0次及び1次モーメントを用いる。量
Rkは截頭した組のモーメントとその組が完全で
あつた場合の予想モーメントとの間の差である。
量Qkがこの後者を表わし、φ1及びφ2の限界の間
のPnの積分が截頭データの組からの測定モーメ
ントを表わす。この差は、延長関数に対するg1
(φ)及びg2(φ)の積分について述べた様に、延
長関数の間で分割しなければならない量である。
式(18)は、延長関数の係数を決定する際にモーメ
ント及びビユーの情報を更に利用し易い形になる
様に操作することが出来る。或る操作の後、0次
及び1次モーメントの両方を使うと、式(7)は次の
様になる。
ント及びビユーの情報を更に利用し易い形になる
様に操作することが出来る。或る操作の後、0次
及び1次モーメントの両方を使うと、式(7)は次の
様になる。
P1 2/2ξ1−P2 2/2ξ2=R0 (8)
φ1P1 2/2ξ1−P1 3/6ξ1 2−φ2P2 2/2ξ2+P2 3/6
ξ2 2=R1(9) ξ≡1/ai (10) ξを消去すると V1ξ2 2+V2ξ2+V3=0 (11) こゝで V1=P2 3/6(1−P2/P1) V2=P2 2/2(φ1−φ2−4/3 R0/P1 (12) V3=−R1+φ1R0−2/3 R0 2/P1 2次式の解を求める普通の式を使い、式(11)をξ2
について解くと ξ2={−V2−√2 2−41 3}/2V1 (13) 最後に、ξ2を使つてξ1、即ち延長関数の勾配の逆
数を解くと、 ξ2=P2 2/P1 2ξ2+2R0/P1 2 (14) 最後に実際の截頭モーメントと、問題の投影を
截頭しない場合の予想モーメントを決定すれば、
単にモーメント情報及び問題のビユーのデータを
使うことにより、任意のビユーをアクセスして完
成することが出来ることは明らかである。最初
に、P1及びP2で表わした点80及び81(第4
図)の振幅、0次及び1次モーメントR0,R1及
び切断点φ1及びφ2に於けるチヤンネルの角度を
使つて、式(12)に示した係数を計算することが必要
である。次にこういう係数を式(13)で使つて、ξ2を
計算し、その後式(14)でξ2を使つてξ1を計算する。
交点は式(4)によつて与えられ、そのビユーの中に
ある情報だけを使つて、関数g1(φ)及びg2(φ)
に対する係数が決定され、截頭データの組が完成
する。
ξ2 2=R1(9) ξ≡1/ai (10) ξを消去すると V1ξ2 2+V2ξ2+V3=0 (11) こゝで V1=P2 3/6(1−P2/P1) V2=P2 2/2(φ1−φ2−4/3 R0/P1 (12) V3=−R1+φ1R0−2/3 R0 2/P1 2次式の解を求める普通の式を使い、式(11)をξ2
について解くと ξ2={−V2−√2 2−41 3}/2V1 (13) 最後に、ξ2を使つてξ1、即ち延長関数の勾配の逆
数を解くと、 ξ2=P2 2/P1 2ξ2+2R0/P1 2 (14) 最後に実際の截頭モーメントと、問題の投影を
截頭しない場合の予想モーメントを決定すれば、
単にモーメント情報及び問題のビユーのデータを
使うことにより、任意のビユーをアクセスして完
成することが出来ることは明らかである。最初
に、P1及びP2で表わした点80及び81(第4
図)の振幅、0次及び1次モーメントR0,R1及
び切断点φ1及びφ2に於けるチヤンネルの角度を
使つて、式(12)に示した係数を計算することが必要
である。次にこういう係数を式(13)で使つて、ξ2を
計算し、その後式(14)でξ2を使つてξ1を計算する。
交点は式(4)によつて与えられ、そのビユーの中に
ある情報だけを使つて、関数g1(φ)及びg2(φ)
に対する係数が決定され、截頭データの組が完成
する。
投影の片側だけが截頭される特別な場合が存在
し得ることに注意されたい。この場合、勾配を見
出すのに、式(3)と1次モーメントの式(7)だけを使
えばよい。投影の左側の端が脱落している場合、
勾配の逆数は次式によつて決定される。
し得ることに注意されたい。この場合、勾配を見
出すのに、式(3)と1次モーメントの式(7)だけを使
えばよい。投影の左側の端が脱落している場合、
勾配の逆数は次式によつて決定される。
ξ1=2R0/P1 2 (15)
これに対して右側の端では、
ξ2=−2R0/P2 2 (16)
従つて、2つの延長関数の係数が完全に決定さ
れ、元の截頭した投影と組合さつて、再構成アル
ゴリズムを少なくとも低い周波数の連続性の条件
を充たす様な、延長した投影が限定される。
れ、元の截頭した投影と組合さつて、再構成アル
ゴリズムを少なくとも低い周波数の連続性の条件
を充たす様な、延長した投影が限定される。
この後、この延長関数を使い、それを標準型の
たゝみ込み波器でたゝみ込んで、画像を再構成
する為の逆向き投影に使われる波データの組を
発生することが出来る。然し、多くの今日のCT
走査器では、たゝみ込みはフーリエ領域で行なわ
れ、最初に速いフーリエ変換を使つた後、乗算及
び逆変換をしてから逆向き投影を行なう。
たゝみ込み波器でたゝみ込んで、画像を再構成
する為の逆向き投影に使われる波データの組を
発生することが出来る。然し、多くの今日のCT
走査器では、たゝみ込みはフーリエ領域で行なわ
れ、最初に速いフーリエ変換を使つた後、乗算及
び逆変換をしてから逆向き投影を行なう。
特定のCT走査器では、循環を避ける為に、約
511個の要素を持つ非截頭データの組に1024点の
FFTを使つた。従つて、上に述べた様に投影の
組を延長することにより、データの組が拡大し
て、1024点のFFTで十分でないと想定するのは
考えられないことでもない。従つて、この延長デ
ータの組に標準型のたゝみ込み方法を使うには、
少なくとも2048点のFFTが必要であり、変換を
行なうのに必要な時間は大体2倍になる。
511個の要素を持つ非截頭データの組に1024点の
FFTを使つた。従つて、上に述べた様に投影の
組を延長することにより、データの組が拡大し
て、1024点のFFTで十分でないと想定するのは
考えられないことでもない。従つて、この延長デ
ータの組に標準型のたゝみ込み方法を使うには、
少なくとも2048点のFFTが必要であり、変換を
行なうのに必要な時間は大体2倍になる。
この発明の特定の1面として、データの組の非
截頭部分に対して普通のFFTによるたゝみ込み
を使い、延長部分に対して近似的な閉じた形の解
を用い、その後2つのたゝみ込みを組合せて、逆
向き投影の為の変更投影の組を発生することによ
り、この難点を避ける。即ち、第6B図で、工程
90が、截頭したものであつてもなくても、各々
の測定した投影の組をたゝみ込む。次に、工程9
1が各々の截頭した組を選択して、延長のたゝみ
込みの為、閉じた形の解を求める。次に工程92
が、各々の截頭した投影の組に対し、閉じた形の
たゝみ込みを工程90のFFTにたゝみ込みと組
合せ、身体の一部分が撮影領域を越えたことを考
慮した更に正確な一組の波データを発生する。
次に工程93が画像を発生するのに使われる通常
の手順、典型的には逆向き投影を遂行する。
截頭部分に対して普通のFFTによるたゝみ込み
を使い、延長部分に対して近似的な閉じた形の解
を用い、その後2つのたゝみ込みを組合せて、逆
向き投影の為の変更投影の組を発生することによ
り、この難点を避ける。即ち、第6B図で、工程
90が、截頭したものであつてもなくても、各々
の測定した投影の組をたゝみ込む。次に、工程9
1が各々の截頭した組を選択して、延長のたゝみ
込みの為、閉じた形の解を求める。次に工程92
が、各々の截頭した投影の組に対し、閉じた形の
たゝみ込みを工程90のFFTにたゝみ込みと組
合せ、身体の一部分が撮影領域を越えたことを考
慮した更に正確な一組の波データを発生する。
次に工程93が画像を発生するのに使われる通常
の手順、典型的には逆向き投影を遂行する。
評価しようとする延長部分の挟み込みは次の様
に表わすことが出来る。
に表わすことが出来る。
qi(φ)=[gi(φ)cosφ]*K(φ),
φ1<φ<φ2 (17)
こゝでK(φ)は核であり、この式は撮影領域
の中、即ちφ1とφ2の間だけで評価しさえすれば
よい。特定の核が選ばれた時、今の例では、
LAKS核と知られるものとを選んだ時、式(17)は閉
じた形で評価することが出来る。最初に、これは
次の様に近似することが出来る。
の中、即ちφ1とφ2の間だけで評価しさえすれば
よい。特定の核が選ばれた時、今の例では、
LAKS核と知られるものとを選んだ時、式(17)は閉
じた形で評価することが出来る。最初に、これは
次の様に近似することが出来る。
qi(φ)=cos[gi(φ)*K(φ)] (18)
次に
g2(φ)=P2(1−φ−φ2/δ2) (19)
に対するフーリエ変換は
G2(k)=P2e−ikφ2
G2(k)=P2e−ikφ2
・{1−coskδ2/R2δ2+i(sinkδ2/k2δ2−1/
k)}(20) 従つて、LAKS核に対して、たゝみ込みは次の
様に表わすことが出来る。
k)}(20) 従つて、LAKS核に対して、たゝみ込みは次の
様に表わすことが出来る。
q2′(φ)≡q2(φ)/P2cosφ2=1/2π・∞ -∞|
k|G2(k)eik〓dk(21) その解は次の様になる。
k|G2(k)eik〓dk(21) その解は次の様になる。
q2(φ)=1/πδ2log[δ2+φ2−φ/φ2−φ]
−1/π(φ2−φ1) (22) これは撮影領域の限界内に対するものである。
−1/π(φ2−φ1) (22) これは撮影領域の限界内に対するものである。
q1に対して同じ様な式を計算することが出来、
撮影領域の外側にある区域のたゝみ込みの、撮影
領域の内部にある要素に対する影響を評価するこ
とが出来る。
撮影領域の外側にある区域のたゝみ込みの、撮影
領域の内部にある要素に対する影響を評価するこ
とが出来る。
扇形のビームの場合に“LAKS”核を使うこと
は厳密には有効ではない。この発明では、誤差が
生じても、この延長過程で回復するのは低い空間
周波数だけであるから、無視し得ることが判つ
た。
は厳密には有効ではない。この発明では、誤差が
生じても、この延長過程で回復するのは低い空間
周波数だけであるから、無視し得ることが判つ
た。
延長投影の最初の2つのモーメントを他のビユ
ーと合う様に作る上述の手順を用いた実験で、截
頭ビユーが原因で起る陰影のアーキテクチユアが
目立つて減少することが判つた。更に、投影の組
の延長のたゝみ込みの評価に閉じた形の解を使う
ことにより、処理時間を許容し難い程長くせず
に、この結果が得られる。
ーと合う様に作る上述の手順を用いた実験で、截
頭ビユーが原因で起る陰影のアーキテクチユアが
目立つて減少することが判つた。更に、投影の組
の延長のたゝみ込みの評価に閉じた形の解を使う
ことにより、処理時間を許容し難い程長くせず
に、この結果が得られる。
第1図はCT走査器の主な素子を示す略図、第
2図は物体が撮影領域を越えることを示す略図、
第3図はこの発明に従つて基準チヤンネル情報を
改善する様子を示すフローチヤート、第4図は截
頭投影の組を示す図、第5図は截頭投影の組に対
する0次及び1次モーメントを示す図、第6図は
この発明を実施して截頭投影の組を完備する場合
を示すフローチヤートである。 主な符号の説明、14:検出器配列、15:患
者、20:撮影領域、21,22:基準チヤンネ
ルX線。
2図は物体が撮影領域を越えることを示す略図、
第3図はこの発明に従つて基準チヤンネル情報を
改善する様子を示すフローチヤート、第4図は截
頭投影の組を示す図、第5図は截頭投影の組に対
する0次及び1次モーメントを示す図、第6図は
この発明を実施して截頭投影の組を完備する場合
を示すフローチヤートである。 主な符号の説明、14:検出器配列、15:患
者、20:撮影領域、21,22:基準チヤンネ
ルX線。
Claims (1)
- 【特許請求の範囲】 1 身体15を該身体の周りの複数個の角度の放
射に対して露出し、前記複数個の角度で身体を通
過した放射を検出する手段14を有して、その投
影の幾つかは身体が撮影領域を越えることによつ
て影響を受ける投影の組を作り、各々の投影には
身体に関連する基準チヤンネル情報21,22が
有るCT像を補償する装置に於て、 予定のレベルより低い基準チヤンネル情報を検
出して、基準チヤンネル情報が身体の影響を受け
た投影を識別する手段32,33,36と、 影響を受けなかつた投影から典型的な基準チヤ
ンネル・レベルを決定する手段41と、 前記典型的な基準チヤンネル・レベルを前記識
別された投影の基準チヤンネルに割当てる手段4
3と、 前記基準レベルを割当てて各々の投影を正規化
する手段44と、 一組の正規化された投影から画像を発生する手
段16,17と、 前記一組の投影の0次及び1次モーメントを決
定する手段73と、 該0次モーメントから截頭した投影から脱落し
ている質量の量を決定する手段82と、 1次モーメントから各々の截頭した投影に関係
する延長した投影の図心を決定する手段83と、 各々の截頭した投影に対する脱落した質量及び
図心の情報を用いて前記投影を延長する手段8
4,85と、 各々の投影を重畳積分処理して一組のたゝみ込
み投影を発生する手段90,91,92と、 該一組のたゝみ込み投影から画像を発生する手
段93を含むCT像を補償する装置。 2 特許請求の範囲1に記載した装置に於て、重
畳積分処理の工程が、 速いフーリエ変換を用いて各々の截頭しない投
影並びに各々の截頭した投影の截頭しなかつた部
分のたゝみ込む手段90と、 閉じた形の解を用いて截頭した投影の延長部分
をたゝみ込む手段91と、 各々の截頭した投影の組に対し、閉じた形の解
をたゝみ込みと組合せる手段92を含んでいる装
置。 3 特許請求の範囲1乃至2のいずれか一項に記
載した装置に於て、 前記一組の投影の複数個のモーメントを決定す
る手段73と、 該モーメントから各々の截頭した投影から脱落
した質量の量並びに各々の截頭した投影の脱落し
た質量の分布を決定する手段82,83と、 低次多項式にモーメント及び境界条件をはめて
各々の截頭した投影を延長する手段85と、 各々の投影を核とたゝみ込んで一組のたゝみ込
み投影を発生する手段90,91,92と、 該一組のたゝみ込み投影から画像を発生する手
段93を含むCT像を補償する装置。 4 特許請求の範囲3に記載した装置に於て、複
数個のモーメントを決定する手段が、0次及び1
次モーメントを決定する手段73を含み、 脱落した質量の量を決定する手段82が0次モ
ーメントを使い、 脱落している質量の分布を決定する手段83が
1次モーメントを使う装置。 5 特許請求の範囲3に記載した装置に於て、
たゝみ込む手段が、 速いフーリエ変換を用いて各々の截頭しない投
影並びに各々の截頭した投影の截頭しなかつた部
分をたゝみ込む手段90と、 閉じた形の解を使つて截頭した投影の延長部分
をたゝみ込む手段91と、 各々の截頭した投影に対し閉じた形の解をたゝ
み込みと組合せる手段92を含む装置。 6 特許請求の範囲4乃至5のいずれか一項に記
載した装置に於て、低次多項式をはめる手段85
が、 投影の截頭しなかつた部分と滑らかに合流する
点から延長部分を開始し、 0次及び1次モーメントを用いて前記開始する
点と延長した投影がゼロになる点との間で低次多
項式をはめることを含む装置。 7 特許請求の範囲6に記載した装置に於て、低
次多項式が1次多項式である装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US424501 | 1982-09-27 | ||
US06/424,501 US4550371A (en) | 1982-09-27 | 1982-09-27 | Method and apparatus for compensating CT images for truncated projections |
Publications (2)
Publication Number | Publication Date |
---|---|
JPS5980234A JPS5980234A (ja) | 1984-05-09 |
JPH0446578B2 true JPH0446578B2 (ja) | 1992-07-30 |
Family
ID=23682857
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP58177177A Granted JPS5980234A (ja) | 1982-09-27 | 1983-09-27 | Ct像を補償する装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US4550371A (ja) |
EP (1) | EP0106964B1 (ja) |
JP (1) | JPS5980234A (ja) |
DE (1) | DE3381384D1 (ja) |
IL (1) | IL69092A0 (ja) |
Families Citing this family (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IL69327A (en) * | 1983-07-26 | 1986-11-30 | Elscint Ltd | Automatic misregistration correction |
US4606004A (en) * | 1984-03-21 | 1986-08-12 | General Electric Company | Apparatus for reduction of filtration truncation errors |
JPS60220050A (ja) * | 1984-04-17 | 1985-11-02 | 横河メディカルシステム株式会社 | 計算機トモグラフィ装置 |
US4674045A (en) * | 1984-12-31 | 1987-06-16 | Picker International Inc. | Filter for data processing |
US4809172A (en) * | 1987-03-11 | 1989-02-28 | Kabushiki Kaisha Toshiba | Method for obtaining image data with a tomographic apparatus |
CA1274922A (en) * | 1987-09-11 | 1990-10-02 | Terence Taylor | Region of interest tomography employing a differential scanning technique |
DE8713524U1 (de) * | 1987-10-08 | 1989-02-02 | Siemens AG, 1000 Berlin und 8000 München | Computertomograph der dritten Generation |
US5043890A (en) * | 1989-06-12 | 1991-08-27 | General Electric | Compensation of computed tomography data for objects positioned outside the field of view of the reconstructed image |
US5225980A (en) * | 1991-11-27 | 1993-07-06 | General Electric Company | Reduction of image artifacts from support structures in tomographic imaging |
US5241576A (en) * | 1991-12-23 | 1993-08-31 | General Electric Company | Segmented detector containing sub-elements for separate measuring of a fan beam |
US5485492A (en) * | 1992-03-31 | 1996-01-16 | Lunar Corporation | Reduced field-of-view CT system for imaging compact embedded structures |
US6701000B1 (en) * | 1999-04-30 | 2004-03-02 | General Electric Company | Solution to detector lag problem in a solid state detector |
US6879656B2 (en) * | 2002-07-23 | 2005-04-12 | General Electric Company | Method and apparatus for deriving motion information from projection data |
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 |
WO2005104038A1 (en) * | 2004-04-21 | 2005-11-03 | Philips Intellectual Property & Standards Gmbh | Cone beam ct apparatus using truncated projections and a previously acquired 3d ct image |
US7050528B2 (en) * | 2004-09-30 | 2006-05-23 | Wisconsin Alumni Research Foundation | Correction of CT images for truncated or incomplete projections |
JP2006102299A (ja) * | 2004-10-07 | 2006-04-20 | Ge Medical Systems Global Technology Co Llc | X線線量補正方法およびx線ct装置 |
US7359478B2 (en) * | 2004-11-18 | 2008-04-15 | Toshiba Medical Systems Corporation | Method for restoring truncated helical cone-beam computed tomography data |
JP2008525142A (ja) * | 2004-12-29 | 2008-07-17 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | X線投影のアーチファクト補正のための装置及び方法 |
KR100687846B1 (ko) * | 2005-01-21 | 2007-02-27 | 경희대학교 산학협력단 | 국부 고해상도 엑스선 단층 영상 재구성 방법 및 국부고해상도 엑스선 단층 영상 재구성 장치 |
US7502440B2 (en) * | 2005-04-05 | 2009-03-10 | Kabushiki Toshiba | Radiodiagnostic apparatus |
JP2006312027A (ja) * | 2005-04-05 | 2006-11-16 | Toshiba Corp | 放射線診断装置 |
US20070076933A1 (en) * | 2005-09-30 | 2007-04-05 | Jared Starman | Estimating the 0th and 1st moments in C-arm CT data for extrapolating truncated projections |
US7515678B2 (en) * | 2005-11-23 | 2009-04-07 | General Electric Company | Method and system for performing CT image reconstruction with motion artifact correction |
US7756315B2 (en) * | 2005-11-23 | 2010-07-13 | General Electric Company | Method and apparatus for field-of-view expansion of volumetric CT imaging |
EP1825811B1 (en) * | 2006-02-27 | 2015-08-26 | Kabushiki Kaisha Toshiba | Image display apparatus and x-ray computed tomography apparatus |
DE102006014629A1 (de) * | 2006-03-29 | 2007-10-04 | Siemens Ag | Verfahren zur Korrektur von Trunkierungsartefakten |
JP2008220653A (ja) * | 2007-03-13 | 2008-09-25 | Toshiba Corp | X線ct装置、被検体外形推定方法、画像再構成方法 |
DE102007054079A1 (de) * | 2007-11-13 | 2009-05-20 | Siemens Ag | Verfahren zur Extrapolation abgeschnittener, unvollständiger Projektionen für die Computertomographie |
DE102008003173B4 (de) * | 2008-01-04 | 2016-05-19 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung für die Computertomographie zur Be- |
KR100923094B1 (ko) | 2008-01-15 | 2009-10-22 | (주)이우테크놀로지 | 트렁케이션 아티팩트를 보정하는 방법 |
US8213694B2 (en) * | 2008-02-29 | 2012-07-03 | Barco N.V. | Computed tomography reconstruction from truncated scans |
RU2568321C2 (ru) * | 2009-08-20 | 2015-11-20 | Конинклейке Филипс Электроникс Н.В. | Реконструкция изображения изучаемой области |
DE102009047867B4 (de) * | 2009-09-30 | 2016-10-06 | Siemens Healthcare Gmbh | Verfahren und Vorrichtung zur Korrektur von trunkierten Projektionsdaten |
US8238518B2 (en) * | 2010-06-23 | 2012-08-07 | The Institute Of Cancer Research | Radiotherapy system |
DE102012207629B4 (de) * | 2012-05-08 | 2023-01-05 | Siemens Healthcare Gmbh | CT-Bildrekonstruktion im erweiterten Messfeld |
US9224216B2 (en) * | 2013-07-31 | 2015-12-29 | Kabushiki Kaisha Toshiba | High density forward projector for spatial resolution improvement for medical imaging systems including computed tomography |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5672849A (en) * | 1979-11-21 | 1981-06-17 | Hitachi Medical Corp | Method of regenerating tomographing image |
JPS5689243A (en) * | 1979-11-29 | 1981-07-20 | Technicare Corp | Method for compensating projection of partial region scanning |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1584954A (en) * | 1976-08-17 | 1981-02-18 | Emi Ltd | Radiography |
DE2753260A1 (de) * | 1977-11-30 | 1979-05-31 | Philips Patentverwaltung | Verfahren zur ermittlung der raeumlichen verteilung der absorption einer strahlung |
DE2926456A1 (de) * | 1979-06-30 | 1981-01-15 | Philips Patentverwaltung | Verfahren zur ermittlung des randes eines koerpers mittels am koerper gestreuter strahlung |
DE3012648A1 (de) * | 1980-04-01 | 1981-10-08 | Philips Patentverwaltung Gmbh, 2000 Hamburg | Computer-tomographiegeraet |
FR2519772B1 (fr) * | 1982-01-08 | 1986-01-31 | Thomson Csf | Dispositif de detection de reference pour tomodensitometre multidetecteur et tomodensitometre comportant un tel dispositif |
US4472823A (en) * | 1982-03-17 | 1984-09-18 | U.S. Philips Corporation | Computed tomography apparatus with detector sensitivity correction |
-
1982
- 1982-09-27 US US06/424,501 patent/US4550371A/en not_active Expired - Lifetime
-
1983
- 1983-06-28 IL IL69092A patent/IL69092A0/xx not_active IP Right Cessation
- 1983-08-18 EP EP19830108168 patent/EP0106964B1/en not_active Expired - Lifetime
- 1983-08-18 DE DE8383108168T patent/DE3381384D1/de not_active Expired - Fee Related
- 1983-09-27 JP JP58177177A patent/JPS5980234A/ja active Granted
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5672849A (en) * | 1979-11-21 | 1981-06-17 | Hitachi Medical Corp | Method of regenerating tomographing image |
JPS5689243A (en) * | 1979-11-29 | 1981-07-20 | Technicare Corp | Method for compensating projection of partial region scanning |
Also Published As
Publication number | Publication date |
---|---|
EP0106964B1 (en) | 1990-03-28 |
IL69092A0 (en) | 1983-10-31 |
US4550371A (en) | 1985-10-29 |
JPS5980234A (ja) | 1984-05-09 |
EP0106964A3 (en) | 1986-02-12 |
EP0106964A2 (en) | 1984-05-02 |
DE3381384D1 (de) | 1990-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPH0446578B2 (ja) | ||
US4590558A (en) | Method and apparatus for removing objects from CT images | |
US6094467A (en) | Method for improving CT images having high attenuation objects | |
CA1074460A (en) | Device for measuring the spatial distribution of radiation absorption in a slice of a body | |
US6415013B1 (en) | Backprojection methods and apparatus for computed tomography imaging systems | |
US6219441B1 (en) | Reconstruction of images from three-dimensional cone beam data | |
US4703424A (en) | Method and means of correcting for a shift in the center of rotation of a rotating fan beam CT system | |
US20150124927A1 (en) | Dark field computed tomography imaging | |
US7280630B2 (en) | Calculation of additional projection data from projection data acquired with a divergent beam | |
US20040073584A1 (en) | Methods and apparatus for truncation compensation | |
US6628744B1 (en) | Off-focal radiation correction in CT | |
Hsieh | Nonstationary noise characteristics of the helical scan and its impact on image quality and artifacts | |
US20040240604A1 (en) | Methods and devices for CT reconstruction using a grangeat approach | |
US4228505A (en) | Method for computed tomography | |
US7379527B2 (en) | Methods and apparatus for CT calibration | |
Zou et al. | Partial volume and aliasing artefacts in helical cone-beam CT | |
JPS582695B2 (ja) | 透過性放射線による検査装置 | |
US6980681B1 (en) | Methods and apparatus for helical reconstruction for multislice CT scan | |
US6084936A (en) | Almost-everywhere extrapolation from cone-beam data | |
US20210295574A1 (en) | Motion compensated reconstruction for helical computer tomography | |
US6411671B2 (en) | Method for reducing line artifacts in a CT image and device for implementing the method | |
JPH0966051A (ja) | X線断層撮影装置 | |
NL193556C (nl) | Computer-tomogram-scanner. | |
Rathee et al. | Image restoration in computed tomography: Estimation of the spatially variant point spread function | |
JP2003135450A (ja) | X線ct再構成画像におけるアーチファクトの低減方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
LAPS | Cancellation because of no payment of annual fees |