JP2002224099A - 螺旋状の相対運動を含むコンピュータ断層撮影方法 - Google Patents
螺旋状の相対運動を含むコンピュータ断層撮影方法Info
- Publication number
- JP2002224099A JP2002224099A JP2001373109A JP2001373109A JP2002224099A JP 2002224099 A JP2002224099 A JP 2002224099A JP 2001373109 A JP2001373109 A JP 2001373109A JP 2001373109 A JP2001373109 A JP 2001373109A JP 2002224099 A JP2002224099 A JP 2002224099A
- Authority
- JP
- Japan
- Prior art keywords
- link
- attenuation
- values
- ray
- radiation
- 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.)
- Ceased
Links
Classifications
-
- G06T12/20—
-
- 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (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)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
(57)【要約】
【課題】 円錐状の放射線ビームが回転軸に垂直な方向
及び平行な方向に大きい開口角を有する場合でも少ない
画像アーティファクトで画像が再現されるようコンピュ
ータ断層撮影を行うことを目的とする。 【解決手段】 本発明はフィルタリングされた逆投影に
必要とされる加算演算が２段階で行われ、フィルタリン
グが第１の加算段階の後にのみ行われるコンピュータ断
層撮影方法に関する。この方法は、平行な光線によって
決められる表面上の平行な平面上に位置し螺旋の対向し
て位置するセグメントをリンクする点について行われる
ことが望ましい。
及び平行な方向に大きい開口角を有する場合でも少ない
画像アーティファクトで画像が再現されるようコンピュ
ータ断層撮影を行うことを目的とする。 【解決手段】 本発明はフィルタリングされた逆投影に
必要とされる加算演算が２段階で行われ、フィルタリン
グが第１の加算段階の後にのみ行われるコンピュータ断
層撮影方法に関する。この方法は、平行な光線によって
決められる表面上の平行な平面上に位置し螺旋の対向し
て位置するセグメントをリンクする点について行われる
ことが望ましい。
Description
【０００１】
【発明の属する技術分野】本発明は、円錐状放射線ビー
ムを放出する放射線源が検査ゾーンに対して螺旋状の相
対運動を行うコンピュータ断層撮影方法に関する。ま
た、本発明は、かかる方法を実行するためのコンピュー
タ断層撮影装置及びコンピュータプログラムに関する。
ムを放出する放射線源が検査ゾーンに対して螺旋状の相
対運動を行うコンピュータ断層撮影方法に関する。ま
た、本発明は、かかる方法を実行するためのコンピュー
タ断層撮影装置及びコンピュータプログラムに関する。
【０００２】
【従来の技術】上述の種類のコンピュータ断層撮影方法
は、ＰＣＴ／ＩＢ９９／０００２７より公知である。Ｃ
Ｔ画像は、測定値がいわゆる平行リビニング（rebinnin
g）演算の後に１次元の傾斜状の（ramp-like）フィルタ
リングを受けるフィルタ補正逆投影法を適用することに
より、捕捉された測定値から再構成される。続いて、検
査ゾーン中の所与の点を通過した光線に関連するフィル
タデータが加算され、当該の点の放射線ビームの光線に
ついての減衰係数が与えられる。
は、ＰＣＴ／ＩＢ９９／０００２７より公知である。Ｃ
Ｔ画像は、測定値がいわゆる平行リビニング（rebinnin
g）演算の後に１次元の傾斜状の（ramp-like）フィルタ
リングを受けるフィルタ補正逆投影法を適用することに
より、捕捉された測定値から再構成される。続いて、検
査ゾーン中の所与の点を通過した光線に関連するフィル
タデータが加算され、当該の点の放射線ビームの光線に
ついての減衰係数が与えられる。
【０００３】
【発明が解決しようとする課題】このようにして形成さ
れたＣＴ（ＣＴ＝コンピュータ断層撮影）画像は非常に
よい画質であるが、特に円錐状の放射線ビームが回転軸
に垂直な方向及び平行な方向に大きい開口角を有する場
合は、やはり画像アーティファクトを含みうる。
れたＣＴ（ＣＴ＝コンピュータ断層撮影）画像は非常に
よい画質であるが、特に円錐状の放射線ビームが回転軸
に垂直な方向及び平行な方向に大きい開口角を有する場
合は、やはり画像アーティファクトを含みうる。
【０００４】
【課題を解決するための手段】本発明は、上述の種類の
方法によって与えられる画質を更に向上させることを目
的とする。この目的は、（ａ）放射線源から発せられ検
査ゾーン又はその中にある対象を横切る複数の光線を含
む円錐状放射ビームを発生する段階と、（ｂ）放射線源
と検査ゾーン又は対象との間で回転軸回りの回転及び回
転軸に平行な変位を含み螺旋の形状とされる相対的な運
動を発生する段階と、（ｃ）検出器ユニットを用いて、
相対的な運動中の検査ゾーン中の光線の減衰に依存する
測定値を捕捉する段階と、（ｄ）光線の位置及び向きを
表わす３次元パラメータ空間中の網のリンクに沿って測
定値を加算することによりリンク値を計算する段階と、
（ｅ）検査ゾーンの所与の面を通過する光線に関連付け
られるリンクについてのフィルタデータを生成するため
にリンク値をフィルタリングする段階と、（ｆ）面上の
画素における放射線の減衰を、画素を通過する光線によ
ってパラメータ空間中で決められる軌跡を近似するリン
クのフィルタデータを加算することにより計算する段階
と、（ｇ）回転軸の方向上互いにずらされた他の面につ
いて少なくとも段階（ｅ）及び（ｆ）を繰り返す段階と
を含む、コンピュータ断層撮影方法によって達成され
る。
方法によって与えられる画質を更に向上させることを目
的とする。この目的は、（ａ）放射線源から発せられ検
査ゾーン又はその中にある対象を横切る複数の光線を含
む円錐状放射ビームを発生する段階と、（ｂ）放射線源
と検査ゾーン又は対象との間で回転軸回りの回転及び回
転軸に平行な変位を含み螺旋の形状とされる相対的な運
動を発生する段階と、（ｃ）検出器ユニットを用いて、
相対的な運動中の検査ゾーン中の光線の減衰に依存する
測定値を捕捉する段階と、（ｄ）光線の位置及び向きを
表わす３次元パラメータ空間中の網のリンクに沿って測
定値を加算することによりリンク値を計算する段階と、
（ｅ）検査ゾーンの所与の面を通過する光線に関連付け
られるリンクについてのフィルタデータを生成するため
にリンク値をフィルタリングする段階と、（ｆ）面上の
画素における放射線の減衰を、画素を通過する光線によ
ってパラメータ空間中で決められる軌跡を近似するリン
クのフィルタデータを加算することにより計算する段階
と、（ｇ）回転軸の方向上互いにずらされた他の面につ
いて少なくとも段階（ｅ）及び（ｆ）を繰り返す段階と
を含む、コンピュータ断層撮影方法によって達成され
る。
【０００５】本発明は、公知の方法によって生ずるアー
ティファクトは、共通な（一次元であり傾斜状の）フィ
ルタ演算を受ける測定値の組合せがフィルタ演算毎に異
なるという認識によるものである。しかしながら、本発
明によれば種々のフィルタ演算は、検査ゾーン中の同一
の面を少なくとも近似的に通過する光線に関連付けられ
る測定値に対してのみ行われる。
ティファクトは、共通な（一次元であり傾斜状の）フィ
ルタ演算を受ける測定値の組合せがフィルタ演算毎に異
なるという認識によるものである。しかしながら、本発
明によれば種々のフィルタ演算は、検査ゾーン中の同一
の面を少なくとも近似的に通過する光線に関連付けられ
る測定値に対してのみ行われる。
【０００６】フィルタ演算をこれらの測定値に限ること
が可能であるのは、段階（ｄ）及び（ｆ）に従って２ス
テップ加算演算が実行されるため、また、これらの２つ
の段階の間に段階（ｅ）に従ってフィルタ演算が実行さ
れるためである。
が可能であるのは、段階（ｄ）及び（ｆ）に従って２ス
テップ加算演算が実行されるため、また、これらの２つ
の段階の間に段階（ｅ）に従ってフィルタ演算が実行さ
れるためである。
【０００７】請求項２は、本発明の望ましい変形例を開
示する。請求項２に定義される面を通過する全ての光線
は、回転軸に平行に延び共に正確に１８０°の角度範囲
を決める面上に位置する。従って、この面での減衰は、
冗長な測定データを使用することなく再構成されうる。
他の全ての面について、より大きい角度範囲から光線を
考慮する必要がある。
示する。請求項２に定義される面を通過する全ての光線
は、回転軸に平行に延び共に正確に１８０°の角度範囲
を決める面上に位置する。従って、この面での減衰は、
冗長な測定データを使用することなく再構成されうる。
他の全ての面について、より大きい角度範囲から光線を
考慮する必要がある。
【０００８】本願の出願人は雑誌"International Journ
al of Imaging Systems and Technology", Vol. 11, 20
00(pp.91-100)に、（逆投影よりも先に行われる）１次
元フィルタリングが傾斜線に沿った回転軸を含む平面を
通る光線の測定値を使用するフィルタリングされた逆投
影の方法について発表している。これらの線は、平面上
にＰｉ（パイ）面（請求項２に定義される面と同じであ
る）と称される面の投影を近似する。しかしながら、投
影は幾つかの投影方向についてのみ線形である。従っ
て、一連のかかる面について毎回２次元に減衰が再構成
される場合は、最適な画質は得られない。更に、かかる
方法ではフィルタリング及び逆投影は従来と同様の順序
で行われ、一方、本発明ではこれらの２つの段階は上述
の３つの段階（ｄ）、（ｅ）、（ｆ）を含むものであり
間に挟まれている。
al of Imaging Systems and Technology", Vol. 11, 20
00(pp.91-100)に、（逆投影よりも先に行われる）１次
元フィルタリングが傾斜線に沿った回転軸を含む平面を
通る光線の測定値を使用するフィルタリングされた逆投
影の方法について発表している。これらの線は、平面上
にＰｉ（パイ）面（請求項２に定義される面と同じであ
る）と称される面の投影を近似する。しかしながら、投
影は幾つかの投影方向についてのみ線形である。従っ
て、一連のかかる面について毎回２次元に減衰が再構成
される場合は、最適な画質は得られない。更に、かかる
方法ではフィルタリング及び逆投影は従来と同様の順序
で行われ、一方、本発明ではこれらの２つの段階は上述
の３つの段階（ｄ）、（ｅ）、（ｆ）を含むものであり
間に挟まれている。
【０００９】請求項３に記載の更なる変形例は、検出器
ユニットが回転軸の方向に比較的大きい寸法を有する場
合に特に有利である。小さい寸法の場合は、全ての光線
について余弦関数が適切に値１を近似するため、なくさ
れてもよい。
ユニットが回転軸の方向に比較的大きい寸法を有する場
合に特に有利である。小さい寸法の場合は、全ての光線
について余弦関数が適切に値１を近似するため、なくさ
れてもよい。
【００１０】請求項４に従って正則デカルトグリッドの
グリッド点において減衰値が得られる場合、検査ゾーン
中の任意の面のＣＴ画像は単に後の段階において形成さ
れる。
グリッド点において減衰値が得られる場合、検査ゾーン
中の任意の面のＣＴ画像は単に後の段階において形成さ
れる。
【００１１】外見上、フィルタリングの前にリンク値の
計算のために測定値が得られる角度範囲が大きいほど、
データはより特定的に１つのＰｉ面にのみ関連する。一
方、元の投影データを従来の方法でフィルタリングしな
いための要件はリンクが直線の線分であることである。
角度的な範囲が大きいほど、測定空間中の画素軌跡の片
ごとの線形的な近似の精度は低くなる。これらの矛盾し
た要件の間の適切な妥協策については、請求項５に記載
の変形によって提供される。
計算のために測定値が得られる角度範囲が大きいほど、
データはより特定的に１つのＰｉ面にのみ関連する。一
方、元の投影データを従来の方法でフィルタリングしな
いための要件はリンクが直線の線分であることである。
角度的な範囲が大きいほど、測定空間中の画素軌跡の片
ごとの線形的な近似の精度は低くなる。これらの矛盾し
た要件の間の適切な妥協策については、請求項５に記載
の変形によって提供される。
【００１２】光線は、請求項６に従っていわゆる平行リ
ビニングを受けることが有利である。
ビニングを受けることが有利である。
【００１３】請求項７は本発明の方法を実行するのに適
したコンピュータ断層撮影装置について記載したもので
あり、請求項８は本発明による再構成方法を実行するた
めのコンピュータプログラムについて記載したものであ
る。
したコンピュータ断層撮影装置について記載したもので
あり、請求項８は本発明による再構成方法を実行するた
めのコンピュータプログラムについて記載したものであ
る。
【００１４】
【発明の実施の形態】以下、添付の図面を参照して本発
明について詳述する。図１に示されるコンピュータ断層
撮影装置は、回転軸１４回りに回転可能なガントリー１
の形状の走査ユニットを含む。このために、ガントリー
はモータ２によって望ましくは一定であるが調整可能な
角速度で駆動される。例えばＸ線源である放射線源Ｓ
は、ガントリー１上に取り付けられる。源には、放射線
源Ｓによって生成される放射線から円錐状の放射線ビー
ム４を、即ち、回転軸に垂直な平面上に且つ回転軸の方
向上にゼロ以外の有限の寸法を有する放射線ビームを形
成するコリメータ装置３が設けられる。
明について詳述する。図１に示されるコンピュータ断層
撮影装置は、回転軸１４回りに回転可能なガントリー１
の形状の走査ユニットを含む。このために、ガントリー
はモータ２によって望ましくは一定であるが調整可能な
角速度で駆動される。例えばＸ線源である放射線源Ｓ
は、ガントリー１上に取り付けられる。源には、放射線
源Ｓによって生成される放射線から円錐状の放射線ビー
ム４を、即ち、回転軸に垂直な平面上に且つ回転軸の方
向上にゼロ以外の有限の寸法を有する放射線ビームを形
成するコリメータ装置３が設けられる。
【００１５】放射線ビーム４は、患者が患者台上に（い
ずれも図示せず）載せられうる検査ゾーン１３を通過す
る。検査ゾーン１３は円筒形状とされる。Ｘ線ビーム４
は、円筒状の検査ゾーンを通った後、ガントリー１に接
続されマトリックス状に配置される多数の検出器素子を
含む２次元検出器ユニット１６に入射する。各検出器素
子は、放射線源の各位置の放射線ビーム４の光線につい
ての測定値を生じさせうる。検出器素子は、行及び列に
配置される。検出器の列は回転軸に平行に延びる。検出
器の行は、回転軸に垂直に延びる平面上、例えば放射線
源Ｓの回りの円の円弧の上に配置されうる。検出器の行
は、通常は検出器の行に含まれる検出器素子の数（例え
ば１６個）よりもかなり多くの検出器素子（例えば１０
００個）を含みうる。
ずれも図示せず）載せられうる検査ゾーン１３を通過す
る。検査ゾーン１３は円筒形状とされる。Ｘ線ビーム４
は、円筒状の検査ゾーンを通った後、ガントリー１に接
続されマトリックス状に配置される多数の検出器素子を
含む２次元検出器ユニット１６に入射する。各検出器素
子は、放射線源の各位置の放射線ビーム４の光線につい
ての測定値を生じさせうる。検出器素子は、行及び列に
配置される。検出器の列は回転軸に平行に延びる。検出
器の行は、回転軸に垂直に延びる平面上、例えば放射線
源Ｓの回りの円の円弧の上に配置されうる。検出器の行
は、通常は検出器の行に含まれる検出器素子の数（例え
ば１６個）よりもかなり多くの検出器素子（例えば１０
００個）を含みうる。
【００１６】放射線ビーム４の開口角αmax（開口角
は、放射線ビーム４の縁に位置するビーム４の光線が放
射線源Ｓ及び回転軸１４によって決まる平面に対して成
す角度と定義される）は、測定値の捕捉中に検査される
べき対象を配置すべき円筒の直径を決める。検査される
べき対象、又は患者台もまた、モータ５によって回転軸
１４に平行に変位されうる。この変位の速度は望ましく
は一定であるが、調整可能である。２つのモータ５及び
２を同時に作動させると、放射線源Ｓ及び検出器ユニッ
ト１６は螺旋走査運動を行う。
は、放射線ビーム４の縁に位置するビーム４の光線が放
射線源Ｓ及び回転軸１４によって決まる平面に対して成
す角度と定義される）は、測定値の捕捉中に検査される
べき対象を配置すべき円筒の直径を決める。検査される
べき対象、又は患者台もまた、モータ５によって回転軸
１４に平行に変位されうる。この変位の速度は望ましく
は一定であるが、調整可能である。２つのモータ５及び
２を同時に作動させると、放射線源Ｓ及び検出器ユニッ
ト１６は螺旋走査運動を行う。
【００１７】回転するガントリー１上の検出器ユニット
１６によって捕捉される測定値は、通常は空間中の定点
に配置され非接触で動作するデータループ（図示せず）
を介して検出器ユニットに接続される画像処理コンピュ
ータ１０に与えられる。画像処理コンピュータ１０は、
種々の画像処理演算を実行することが可能である。特
に、捕捉された測定値から検査ゾーン１３中のＸ線の減
衰を再構成し、それにより３Ｄデータセットを生成する
ことが可能である。
１６によって捕捉される測定値は、通常は空間中の定点
に配置され非接触で動作するデータループ（図示せず）
を介して検出器ユニットに接続される画像処理コンピュ
ータ１０に与えられる。画像処理コンピュータ１０は、
種々の画像処理演算を実行することが可能である。特
に、捕捉された測定値から検査ゾーン１３中のＸ線の減
衰を再構成し、それにより３Ｄデータセットを生成する
ことが可能である。
【００１８】以下、図２に示すフローチャートを参照し
て本発明について詳述する。初期化（ブロック１０１）
の後、放射線源Ｓが検査されるべき対象に対して螺旋状
の路に沿って動くよう、モータ２及び５を同時に作動さ
せる。ステップ１０２では、検出器ユニット１６は、た
ぶん平滑化し対数をとった後は、測定値を測定した光線
に沿った減衰の線積分に対応する測定値を捕捉する。３
次元（β，γ，ｓ）パラメータ空間では、かかる測定値
ｐ（β，γ，ｓ）又は関連する光線は、当該の放射線源
の位置から回転軸１４への垂線の方向β（βは１回転以
上した後は２πよりも大きくなりうる）と、当該の測定
値に関連する光線が上記の垂線を含む平面に対して成す
角度γと、（上記の垂線に対して垂直であり、回転軸を
含む平面上の）光線の高さ座標ｓとによって特徴付けら
れる。各光線は、この３次元パラメータ空間中の点によ
って特徴付けられる。
て本発明について詳述する。初期化（ブロック１０１）
の後、放射線源Ｓが検査されるべき対象に対して螺旋状
の路に沿って動くよう、モータ２及び５を同時に作動さ
せる。ステップ１０２では、検出器ユニット１６は、た
ぶん平滑化し対数をとった後は、測定値を測定した光線
に沿った減衰の線積分に対応する測定値を捕捉する。３
次元（β，γ，ｓ）パラメータ空間では、かかる測定値
ｐ（β，γ，ｓ）又は関連する光線は、当該の放射線源
の位置から回転軸１４への垂線の方向β（βは１回転以
上した後は２πよりも大きくなりうる）と、当該の測定
値に関連する光線が上記の垂線を含む平面に対して成す
角度γと、（上記の垂線に対して垂直であり、回転軸を
含む平面上の）光線の高さ座標ｓとによって特徴付けら
れる。各光線は、この３次元パラメータ空間中の点によ
って特徴付けられる。
【００１９】必要であれば、この段階において、全ての
測定値は、関連する光線が回転軸に垂直に延びる平面に
対して成す角度の余弦で重み付けされうる。しかしなが
ら、余弦が全ての光線について略１の値を有する場合、
この段階は省かれうる。これは、検出器１６の寸法が放
射線源からの距離と比較してもはや無視できるほど小さ
くない場合に必要となる。
測定値は、関連する光線が回転軸に垂直に延びる平面に
対して成す角度の余弦で重み付けされうる。しかしなが
ら、余弦が全ての光線について略１の値を有する場合、
この段階は省かれうる。これは、検出器１６の寸法が放
射線源からの距離と比較してもはや無視できるほど小さ
くない場合に必要となる。
【００２０】次の段階（１０３）において、以下の式、
【００２１】
【数１】 に従っていわゆる平行リビニング演算が実行されうる。
【００２２】このようにして生成された新しい組の測定
値ｐＰ（θ，ｔ，ｓ）は回転軸１４に垂直な平面上の当
該の測定値に関連付けられる光線の（投影）方向θと、
この光線と回転軸１４との間の距離ｔと、上述の座標ｓ
とによって特徴付けられる。このようにして再ソート及
び再補間によって得られる新しい測定値は元の測定値ｐ
（β，γ，ｓ）と同等であるが、以下、明瞭性のため投
影値と称するものとし、上付きの添え字Ｐを付す。この
測定値に関連付けられる光線は、（θ，ｔ，ｓ）パラメ
ータ空間中で正則グリッドを定義する。
値ｐＰ（θ，ｔ，ｓ）は回転軸１４に垂直な平面上の当
該の測定値に関連付けられる光線の（投影）方向θと、
この光線と回転軸１４との間の距離ｔと、上述の座標ｓ
とによって特徴付けられる。このようにして再ソート及
び再補間によって得られる新しい測定値は元の測定値ｐ
（β，γ，ｓ）と同等であるが、以下、明瞭性のため投
影値と称するものとし、上付きの添え字Ｐを付す。この
測定値に関連付けられる光線は、（θ，ｔ，ｓ）パラメ
ータ空間中で正則グリッドを定義する。
【００２３】上述の文献により公知の再構成方法及び本
発明による再構成方法について、以下、簡単化された例
に基づいて詳述する。簡単化は、検出器ユニットがｓ＝
０であるような検出器の単一の行のみを含むとの仮定に
基づく。公知の方法に従って、リビニング演算の後に、
同じ投影角度θを有し異なるｔを有する全ての測定値が
共通のフィルタリング演算を受けるフィルタリング段階
が行われる。
発明による再構成方法について、以下、簡単化された例
に基づいて詳述する。簡単化は、検出器ユニットがｓ＝
０であるような検出器の単一の行のみを含むとの仮定に
基づく。公知の方法に従って、リビニング演算の後に、
同じ投影角度θを有し異なるｔを有する全ての測定値が
共通のフィルタリング演算を受けるフィルタリング段階
が行われる。
【００２４】フィルタリングされたデータは、その後に
のみ空間へ逆投影され、検査ゾーン中の点について減衰
値を計算するために、当該の点を通過した全ての光線か
らの寄与はその後に加算される。この段階について、投
影角度θを横座標とし、縦座標ｔを有する２次元（θ，
ｔ）パラメータ空間を示す図３を参照する。パラメータ
空間中の太い破線は、０乃至πの投影角度範囲θに亘っ
て当該の点を通過する光線を相互接続する軌跡３０を示
す。これらの光線は、通常は、投影値ｐＰが置かれる
（θ，ｔ）パラメータ空間中のグリッド点のうちの１つ
には一致しない。従って、逆投影のために使用される値
は補間によって決定されねばならない。軌跡３０が正弦
状の進路であるため、図３の表現はシノグラムと称され
ることが多い。
のみ空間へ逆投影され、検査ゾーン中の点について減衰
値を計算するために、当該の点を通過した全ての光線か
らの寄与はその後に加算される。この段階について、投
影角度θを横座標とし、縦座標ｔを有する２次元（θ，
ｔ）パラメータ空間を示す図３を参照する。パラメータ
空間中の太い破線は、０乃至πの投影角度範囲θに亘っ
て当該の点を通過する光線を相互接続する軌跡３０を示
す。これらの光線は、通常は、投影値ｐＰが置かれる
（θ，ｔ）パラメータ空間中のグリッド点のうちの１つ
には一致しない。従って、逆投影のために使用される値
は補間によって決定されねばならない。軌跡３０が正弦
状の進路であるため、図３の表現はシノグラムと称され
ることが多い。
【００２５】本発明は異なったアプローチをとる。０乃
至πの投影角度範囲に亘って逆投影中に行われる加算演
算は２つのサブステップで行われ、フィルタリング段階
は逆投影の前ではなく逆投影の２つのサブステップの間
に行われる。このため、パラメータ空間中に網４０が形
成される。この網のノードは投影値が決められる立方グ
リッドのグリッド点の距離よりもかなり大きいθ方向及
びｔ方向上の距離に配置される。隣接するノードの列
は、リンク５１によって相互接続される。これらのリン
ク５１のうちの幾つかは、シノグラム中の軌跡３０の少
なくとも最大の勾配に等しい勾配を有する。このように
全ての軌跡は一組のリンクによって近似されうる。図３
中の太い線は軌跡３０を近似するリンクを示す。
至πの投影角度範囲に亘って逆投影中に行われる加算演
算は２つのサブステップで行われ、フィルタリング段階
は逆投影の前ではなく逆投影の２つのサブステップの間
に行われる。このため、パラメータ空間中に網４０が形
成される。この網のノードは投影値が決められる立方グ
リッドのグリッド点の距離よりもかなり大きいθ方向及
びｔ方向上の距離に配置される。隣接するノードの列
は、リンク５１によって相互接続される。これらのリン
ク５１のうちの幾つかは、シノグラム中の軌跡３０の少
なくとも最大の勾配に等しい勾配を有する。このように
全ての軌跡は一組のリンクによって近似されうる。図３
中の太い線は軌跡３０を近似するリンクを示す。
【００２６】第１の部分加算演算は、網の各リンクに沿
った投影値を計算し加算することによってリンク値を計
算することによって行われる。これについて、図３のシ
ノグラムの一部と、点（θｉ１，ｔ１及びθｉ２，
ｔ２）をリンクする単一のリンク５１とを示す図４を参
照して説明する。値ｐＰが定義される幾つかのグリッド
点は、図４中の十字によって示される。図４に示される
ように、グリッド点は各投影角度増分について得られる
リンク上の点に一致しない。リンク５１の端点もまたグ
リッド点と一致する必要はない。点（θｉ１，ｔ１）と
点（θｉ２，ｔ２）の間のリンク５１についてのリンク
値Ｉ（θｉ１，ｔ１；θｉ２，ｔ２）の計算は、以下の
式、
った投影値を計算し加算することによってリンク値を計
算することによって行われる。これについて、図３のシ
ノグラムの一部と、点（θｉ１，ｔ１及びθｉ２，
ｔ２）をリンクする単一のリンク５１とを示す図４を参
照して説明する。値ｐＰが定義される幾つかのグリッド
点は、図４中の十字によって示される。図４に示される
ように、グリッド点は各投影角度増分について得られる
リンク上の点に一致しない。リンク５１の端点もまたグ
リッド点と一致する必要はない。点（θｉ１，ｔ１）と
点（θｉ２，ｔ２）の間のリンク５１についてのリンク
値Ｉ（θｉ１，ｔ１；θｉ２，ｔ２）の計算は、以下の
式、
【００２７】
【数２】 に従う。
【００２８】まだ説明されておらず、続いて実行される
フィルタリング段階は、リンク値からフィルタリングさ
れたデータを生成する。図５に従った第２の加算段階
中、軌跡３０を近似するリンクのフィルタデータ
Ａ１...Ｄ１及びＡ２...Ｄ２が加算される。図５は、軌
跡の一部及びリンクのみを示す。以下画素と称され軌跡
３０に関連付けられる点（ｘ，ｙ）の減衰値ｆ（ｘ，
ｙ）は、以下の式、
フィルタリング段階は、リンク値からフィルタリングさ
れたデータを生成する。図５に従った第２の加算段階
中、軌跡３０を近似するリンクのフィルタデータ
Ａ１...Ｄ１及びＡ２...Ｄ２が加算される。図５は、軌
跡の一部及びリンクのみを示す。以下画素と称され軌跡
３０に関連付けられる点（ｘ，ｙ）の減衰値ｆ（ｘ，
ｙ）は、以下の式、
【００２９】
【数３】 に従って計算される。
【００３０】式中、Ａj，Ｂj，...Ｄjは軌跡３０が近似
されるリンクについてのフィルタリングされたデータで
あり、ｋは網４０によってシノグラムが分割される投影
間隔の数である。値ｗj及びｗj+1はリンクと軌跡３０の
間の距離に依存する重み付け係数であり、その意味は図
５から明らかである。投影間隔の数ｋは、投影間隔の長
さに反比例して変化する。ｋ値が大きい場合は、網のメ
ッシュが細かいこと及びリンクにより軌跡が適切に近似
されることに対応する。しかしながら、これらはかなり
の計算努力を必要とする。精度と計算努力との間の適切
な妥協策は、以下の式、
されるリンクについてのフィルタリングされたデータで
あり、ｋは網４０によってシノグラムが分割される投影
間隔の数である。値ｗj及びｗj+1はリンクと軌跡３０の
間の距離に依存する重み付け係数であり、その意味は図
５から明らかである。投影間隔の数ｋは、投影間隔の長
さに反比例して変化する。ｋ値が大きい場合は、網のメ
ッシュが細かいこと及びリンクにより軌跡が適切に近似
されることに対応する。しかしながら、これらはかなり
の計算努力を必要とする。精度と計算努力との間の適切
な妥協策は、以下の式、
【００３１】
【数４】 によって与えられる。式中、Ｎθは、１８０°に亘って
規則的に分布され、それについて３次元パラメータ空間
が決められうる投影値の数である。
規則的に分布され、それについて３次元パラメータ空間
が決められうる投影値の数である。
【００３２】個々のリンクもまた他の軌跡による軌跡の
近似のために使用されうるため、１回捕捉された関連す
るフィルタデータは数回使用され、それにより必要な計
算努力を減少させる。以下、この方法を３次元データセ
ットの場合について説明するが、３次元の場合は、図３
の２次元の場合の測定によってなんらかの方法で強調さ
れる面はない。それでもなお、ステップ１０４では、放
射線の減衰の画素が再構成されるべき画素について面が
与えられる。
近似のために使用されうるため、１回捕捉された関連す
るフィルタデータは数回使用され、それにより必要な計
算努力を減少させる。以下、この方法を３次元データセ
ットの場合について説明するが、３次元の場合は、図３
の２次元の場合の測定によってなんらかの方法で強調さ
れる面はない。それでもなお、ステップ１０４では、放
射線の減衰の画素が再構成されるべき画素について面が
与えられる。
【００３３】これについて、放射線源が下から上へ動く
と想定して、検査ゾーンに対する放射線源の螺旋状の軌
跡１７を示す図６を参照して説明する。例えば、所与の
投影角度について、相互に平行であり、回転軸１４に平
行な平面上で検査ゾーン（図示せず）を横切る幾つかの
扇形ビーム４１，４２．．．４５が示される。これにつ
いて、上側エッジ及び下側エッジは螺旋１７の対向して
位置するセグメントを通ると想定する。検査ユニット
が、螺旋の隣接した巻きの間の領域中の光線についての
み測定値が得られるよう構成されるとき（又はこれらの
光線の測定値のみが評価されるとき）、（引用された文
献から公知であるように）検査ゾーン中の各点が正確に
１８０°の投影角度範囲から照射されることが確実とさ
れる。
と想定して、検査ゾーンに対する放射線源の螺旋状の軌
跡１７を示す図６を参照して説明する。例えば、所与の
投影角度について、相互に平行であり、回転軸１４に平
行な平面上で検査ゾーン（図示せず）を横切る幾つかの
扇形ビーム４１，４２．．．４５が示される。これにつ
いて、上側エッジ及び下側エッジは螺旋１７の対向して
位置するセグメントを通ると想定する。検査ユニット
が、螺旋の隣接した巻きの間の領域中の光線についての
み測定値が得られるよう構成されるとき（又はこれらの
光線の測定値のみが評価されるとき）、（引用された文
献から公知であるように）検査ゾーン中の各点が正確に
１８０°の投影角度範囲から照射されることが確実とさ
れる。
【００３４】図は、回転軸１４を含み、投影方向に垂直
に延びる矩形の表面１８０も示す。矩形の上辺及び下辺
は平行な扇形ビームの上側エッジ光線及び下側エッジ光
線に一致する。左側の辺及び右側の辺は、検査ゾーン１
３の周囲によって決められる（図１参照）。
に延びる矩形の表面１８０も示す。矩形の上辺及び下辺
は平行な扇形ビームの上側エッジ光線及び下側エッジ光
線に一致する。左側の辺及び右側の辺は、検査ゾーン１
３の周囲によって決められる（図１参照）。
【００３５】全ての光線はこのように所与の方向につい
て矩形１８０を通るため、この矩形は以下検出器ウィン
ドウとも称される。他の投影方向については、この検出
器窓の空間中の位置及び向きはそれを通る扇形ビームと
同じ程度まで変化する。放射線源が螺旋１７に沿って上
向きに動くと想定すると、検査ゾーン中の全ての点はま
ず検出器ウィンドウの上側のエッジへ投影され、１８０
°の投影角度の変化後は下側エッジを通る。図６の表現
において検出器ウィンドウ６の上側エッジを通る光線、
例えば光線４４及び４５は、１８０°に亘る投影方向の
増加の後、検出器ウィンドウ１６０の下側エッジ（その
高さにより上向きに移動されている）を通る。
て矩形１８０を通るため、この矩形は以下検出器ウィン
ドウとも称される。他の投影方向については、この検出
器窓の空間中の位置及び向きはそれを通る扇形ビームと
同じ程度まで変化する。放射線源が螺旋１７に沿って上
向きに動くと想定すると、検査ゾーン中の全ての点はま
ず検出器ウィンドウの上側のエッジへ投影され、１８０
°の投影角度の変化後は下側エッジを通る。図６の表現
において検出器ウィンドウ６の上側エッジを通る光線、
例えば光線４４及び４５は、１８０°に亘る投影方向の
増加の後、検出器ウィンドウ１６０の下側エッジ（その
高さにより上向きに移動されている）を通る。
【００３６】ステップ１０４では、螺旋の対向して位置
するセグメントを相互接続し、平行な平面上に延びる光
線によって形成される面７０が予め選択される。この
（平坦でない）面は、（θ，ｔ，ｓ）パラメータ空間で
はこの面を通過するかこの面に位置する光線に関連する
全ての投影値ｐＰは正確に１８０°の投影角度範囲にあ
るため、Ｐｉ面とも称される。
するセグメントを相互接続し、平行な平面上に延びる光
線によって形成される面７０が予め選択される。この
（平坦でない）面は、（θ，ｔ，ｓ）パラメータ空間で
はこの面を通過するかこの面に位置する光線に関連する
全ての投影値ｐＰは正確に１８０°の投影角度範囲にあ
るため、Ｐｉ面とも称される。
【００３７】この面に位置するかこの面を通過する光線
のみが共通フィルタリング演算を受け、最初は減衰の２
次元再構成はこの面に対してのみ実行される。例えば平
坦な面といった他の面もまた選択されうる。しかしなが
ら、その場合、１８０°以上の投影角度範囲から投影値
ｐＰを考慮に入れることが必要であり、これはフィルタ
リング及び再構成をより困難なものとする。
のみが共通フィルタリング演算を受け、最初は減衰の２
次元再構成はこの面に対してのみ実行される。例えば平
坦な面といった他の面もまた選択されうる。しかしなが
ら、その場合、１８０°以上の投影角度範囲から投影値
ｐＰを考慮に入れることが必要であり、これはフィルタ
リング及び再構成をより困難なものとする。
【００３８】図７は、夫々２２．５°だけ相互にずらさ
れた投影方向について、Ｐｉ面中の画素を検出器ウィン
ドウ１６０に投影した場合について示す図である。上述
のように、Ｐｉ面上の画素の投影は、放射線ビームへ入
るときは検出ウィンドウ１６０の上側エッジａに一致
し、出るときは下側エッジｂに一致する。検出器ウィン
ドウの上側エッジａは式ｓ／ｈ＝０．２５に従って水平
に延び、下側エッジｂが式ｓ／ｈ＝−０．２５に従うの
と同様であり、式中、ｈは螺旋の２つの隣接する巻きの
間の距離である。投影のために、Ｐｉ面７０の画素の間
に位置する方向は直線上に投影されるのではなく、水平
な直線ａ及びｂに対して多少傾斜した狭い片（ｃ，ｄ，
ｅ，ｆ等）へ投影される。隣接の投影方向の片は重なり
合う場合がある。
れた投影方向について、Ｐｉ面中の画素を検出器ウィン
ドウ１６０に投影した場合について示す図である。上述
のように、Ｐｉ面上の画素の投影は、放射線ビームへ入
るときは検出ウィンドウ１６０の上側エッジａに一致
し、出るときは下側エッジｂに一致する。検出器ウィン
ドウの上側エッジａは式ｓ／ｈ＝０．２５に従って水平
に延び、下側エッジｂが式ｓ／ｈ＝−０．２５に従うの
と同様であり、式中、ｈは螺旋の２つの隣接する巻きの
間の距離である。投影のために、Ｐｉ面７０の画素の間
に位置する方向は直線上に投影されるのではなく、水平
な直線ａ及びｂに対して多少傾斜した狭い片（ｃ，ｄ，
ｅ，ｆ等）へ投影される。隣接の投影方向の片は重なり
合う場合がある。
【００３９】従って、Ｐｉ面を通過する（θ，ｔ，ｓ）
パラメータ空間中の光線を含む網４０は、もはや（２次
元の場合のように）１つの平面上に位置するのではな
く、曲面上に位置する。明らかに、網のこの部分の
（θ，ｔ）平面上への投影は図３の表現に対応するが、
この網のノードは異なるｓ座標を有する。座標θ１，ｔ
１及びθ２，ｔ２を有する２つの光線（又は網４０のノ
ード）をリンクするリンク５１は、座標ｓ１及びｓ２を
有する。これらの２つの光線は、座標ｘ，ｙを有するＰ
ｉ面上の画素で交わる。投影方向θ１及びθ２からこの
画素（ｘ，ｙ）を検出ウィンドウ１６０へ投影すること
により、リンク５１に関連する座標ｓ１及びｓ ２が生ず
る。図８は、（θ，ｔ，ｓ）空間中の網の上記の部分の
斜視図である。
パラメータ空間中の光線を含む網４０は、もはや（２次
元の場合のように）１つの平面上に位置するのではな
く、曲面上に位置する。明らかに、網のこの部分の
（θ，ｔ）平面上への投影は図３の表現に対応するが、
この網のノードは異なるｓ座標を有する。座標θ１，ｔ
１及びθ２，ｔ２を有する２つの光線（又は網４０のノ
ード）をリンクするリンク５１は、座標ｓ１及びｓ２を
有する。これらの２つの光線は、座標ｘ，ｙを有するＰ
ｉ面上の画素で交わる。投影方向θ１及びθ２からこの
画素（ｘ，ｙ）を検出ウィンドウ１６０へ投影すること
により、リンク５１に関連する座標ｓ１及びｓ ２が生ず
る。図８は、（θ，ｔ，ｓ）空間中の網の上記の部分の
斜視図である。
【００４０】この網についてのリンクのリンク値は、ス
テップ１０５において計算される。式２は、
テップ１０５において計算される。式２は、
【００４１】
【数５】 となる。
【００４２】ステップ１０６において、図８に示される
網の部分に関連するリンクのリンク値は、以下の式、
網の部分に関連するリンクのリンク値は、以下の式、
【００４３】
【数６】 に従ってフィルタリングされ、式中、
【００４４】
【外１】 は畳み込み演算を示し、ｇ（ｔ）は畳み込みにより所望
の傾斜状のフィルタリングを生じさせるフィルタリング
関数である。式５は、同一の投影角度（θi1）から発せ
られ同一の勾配Δｔ＝ｔ2−ｔ1を有する全てのリンクが
（共通の）フィルタリング演算を受ける１次元の畳み込
みを表わす。これは、所与の投影角度θから出るリンク
を示す図９に示され、図中、同じ勾配を有するリンクは
同じ破線又は点線で示される。フィルタリング演算は、
例えば参照番号５１によって示されるリンクのリンク値
を含む。
の傾斜状のフィルタリングを生じさせるフィルタリング
関数である。式５は、同一の投影角度（θi1）から発せ
られ同一の勾配Δｔ＝ｔ2−ｔ1を有する全てのリンクが
（共通の）フィルタリング演算を受ける１次元の畳み込
みを表わす。これは、所与の投影角度θから出るリンク
を示す図９に示され、図中、同じ勾配を有するリンクは
同じ破線又は点線で示される。フィルタリング演算は、
例えば参照番号５１によって示されるリンクのリンク値
を含む。
【００４５】続く段階１０７において、Ｐｉ面７０（図
６）の個々の画素（ｘ，ｙ）の軌跡３０（図８参照）を
近似するリンクのフィルタリングされたデータが加算さ
れる。このＰｉ面の画素の全ての軌跡は直線ｓ／ｈ＝
０．２５及びθ＝θi（図８のθi＝０）から開始する。
全ての軌跡は、座標ｓ／ｈ＝−０．２５及びθ＝θi＋
πによって決められる共通の線上で終端する。計算は式
３に従って行なわれ、値Ａj...Ｄjはステップ１０６に
おいて計算され軌跡３０を近似する
６）の個々の画素（ｘ，ｙ）の軌跡３０（図８参照）を
近似するリンクのフィルタリングされたデータが加算さ
れる。このＰｉ面の画素の全ての軌跡は直線ｓ／ｈ＝
０．２５及びθ＝θi（図８のθi＝０）から開始する。
全ての軌跡は、座標ｓ／ｈ＝−０．２５及びθ＝θi＋
πによって決められる共通の線上で終端する。計算は式
３に従って行なわれ、値Ａj...Ｄjはステップ１０６に
おいて計算され軌跡３０を近似する
【００４６】
【外２】 によって生成される。
【００４７】ステップ１０７は、ステップ１０４で与え
られるＰｉ面７０上の減衰値ｆ（ｘ，ｙ）を生じさせ
る。続くステップ１０９において、異なるＰｉ面が選択
される。この面は例えば投影方向をインクリメントさせ
ることにより先行するＰｉ面とは異なり、回転軸１４の
方向にずれている。ステップ１０５及び１０７はこの新
しいＰｉ面のために再び行われ、続いて更なるＰｉ面の
ために繰り返される。最後に、ｚ方向に互いにずれてお
り、投影角度のインクリメントにより互いに対して回転
している多数のＰｉ面（図１０に図示）について減衰値
が得られる。
られるＰｉ面７０上の減衰値ｆ（ｘ，ｙ）を生じさせ
る。続くステップ１０９において、異なるＰｉ面が選択
される。この面は例えば投影方向をインクリメントさせ
ることにより先行するＰｉ面とは異なり、回転軸１４の
方向にずれている。ステップ１０５及び１０７はこの新
しいＰｉ面のために再び行われ、続いて更なるＰｉ面の
ために繰り返される。最後に、ｚ方向に互いにずれてお
り、投影角度のインクリメントにより互いに対して回転
している多数のＰｉ面（図１０に図示）について減衰値
が得られる。
【００４８】例えば立方グリッドといった正則デカルト
グリッドの点において２つの外側Ｐｉ面によって決めら
れる３次元ゾーン中の減衰ｆ（ｘ，ｙ，ｚ）の空間的な
分布の表現を可能とするため、ステップ１０８において
補間が行われる。同じｘ及びｙ座標について個々のＰｉ
面のために減衰値が既に計算されているため、ｚ方向の
補間のみが必要とされる。続いて方法は終了する（ステ
ップ１１０）。
グリッドの点において２つの外側Ｐｉ面によって決めら
れる３次元ゾーン中の減衰ｆ（ｘ，ｙ，ｚ）の空間的な
分布の表現を可能とするため、ステップ１０８において
補間が行われる。同じｘ及びｙ座標について個々のＰｉ
面のために減衰値が既に計算されているため、ｚ方向の
補間のみが必要とされる。続いて方法は終了する（ステ
ップ１１０）。
【００４９】ステップ１０３における式１に従ってリビ
ニングがステップ１０４における式４に従ったリンク値
の計算にシフトしたとき、リンク値はリビニングを必要
とすることなく、測定値ｐ（β，γ，ｓ）から直接計算
されうる。ステップ１０３におけるリビニングは、この
ように省かれうる。
ニングがステップ１０４における式４に従ったリンク値
の計算にシフトしたとき、リンク値はリビニングを必要
とすることなく、測定値ｐ（β，γ，ｓ）から直接計算
されうる。ステップ１０３におけるリビニングは、この
ように省かれうる。
【００５０】更に、回転軸（又はｚ軸）に垂直に延びる
平面に配置される検出器行を含む検出器の場合、検出器
行が例えば放射線源の周りの縁の円弧上に配置されてい
る場合、値ｓは検出器行については一定ではない。その
場合、３次元パラメータ空間において一定の値ｓを生じ
させるリビニングもまた必要ではない。その代わりに、
式４の各加数に対して補間が行われうる。エイリアシン
グ効果を回避するため、（式４が基づいている）投影角
度の増分は減少されるべきである。
平面に配置される検出器行を含む検出器の場合、検出器
行が例えば放射線源の周りの縁の円弧上に配置されてい
る場合、値ｓは検出器行については一定ではない。その
場合、３次元パラメータ空間において一定の値ｓを生じ
させるリビニングもまた必要ではない。その代わりに、
式４の各加数に対して補間が行われうる。エイリアシン
グ効果を回避するため、（式４が基づいている）投影角
度の増分は減少されるべきである。
【００５１】ステップ１０５におけるリンク値の計算
は、図２に従った方法について想定されるように夫々の
面に限られる必要はない。しかしながら、その場合、必
要な全てのリンク値は（Ｐｉ）面が選択される前に計算
されねばならず、関連するリンク値はフィルタリングさ
れる。
は、図２に従った方法について想定されるように夫々の
面に限られる必要はない。しかしながら、その場合、必
要な全てのリンク値は（Ｐｉ）面が選択される前に計算
されねばならず、関連するリンク値はフィルタリングさ
れる。
【図１】本発明による方法を実行するのに適したコンピ
ュータ断層撮影装置を示す図である。
ュータ断層撮影装置を示す図である。
【図２】本発明による方法を示すフローチャートであ
る。
る。
【図３】光線の位置及び向きを表わす２次元パラメータ
空間を示す図である。
空間を示す図である。
【図４】図３のパラメータ空間におけるリンクの位置を
示す図である。
示す図である。
【図５】画素を通過する光線を接続する軌跡のセグメン
トを近似する複数のかかるリンクを示す図である。
トを近似する複数のかかるリンクを示す図である。
【図６】本発明の前提となる幾何学的条件について示す
図である。
図である。
【図７】回転軸を含む平面上の種々の画素の投影を示す
図である。
図である。
【図８】３次元パラメータ空間をリンクの網と共に示す
図である。
図である。
【図９】共通のフィルタリング演算を受けるリンクを示
す図である。
す図である。
【図１０】放射線源によって表わされる螺旋パスに対す
る再構成された面の位置を示す図である。
る再構成された面の位置を示す図である。
１ ガントリー ２ モータ ３ コリメータ装置 ４ 放射線ビーム ５ モータ １０ 画像処理コンピュータ １３ 検査ゾーン １４ 回転軸 １６ ２次元検出器ユニット
───────────────────────────────────────────────────── フロントページの続き (71)出願人 590000248 Ｇｒｏｅｎｅｗｏｕｄｓｅｗｅｇ １, 5621 ＢＡ Ｅｉｎｄｈｏｖｅｎ， Ｔｈ ｅ Ｎｅｔｈｅｒｌａｎｄｓ (72)発明者 ヘンリク ヴァルデマール トゥールベル スウェーデン国，58212 リンチェピング, フリール・ガタン 25 Ｆターム(参考） 4C093 AA22 BA03 BA10 CA13 FD05 FD12 FE12 FE15
Claims (8)
- 【請求項１】 （ａ）放射線源から発せられ検査ゾーン
又はその中にある対象を横切る複数の光線を含む円錐状
放射ビームを発生する段階と、 （ｂ）上記放射線源と上記検査ゾーン又は上記対象との
間で回転軸回りの回転及び上記回転軸に平行な変位を含
み螺旋の形状とされる相対的な運動を発生する段階と、 （ｃ）検出器ユニットを用いて、上記相対的な運動中の
上記検査ゾーン中の光線の減衰に依存する測定値を捕捉
する段階と、 （ｄ）上記光線の位置及び向きを表わす３次元パラメー
タ空間中の網のリンクに沿って上記測定値を加算するこ
とによりリンク値を計算する段階と、 （ｅ）上記検査ゾーンの所与の面を通過する光線に関連
付けられるリンクについてのフィルタデータを生成する
ために上記リンク値をフィルタリングする段階と、 （ｆ）上記面上の画素における放射線の減衰を、上記画
素を通過する光線によって上記パラメータ空間中で決め
られる軌跡を近似するリンクのフィルタデータを加算す
ることにより計算する段階と、 （ｇ）上記回転軸の方向上互いにずらされた他の面につ
いて少なくとも上記段階（ｅ）及び（ｆ）を繰り返す段
階とを含む、コンピュータ断層撮影方法。 - 【請求項２】 上記各面は、上記回転軸に平行に延びる
面に位置し上記螺旋の対向して位置するセグメントをリ
ンクする夫々の光線の組によって決められることを特徴
とする、請求項１記載のコンピュータ断層撮影方法。 - 【請求項３】 上記測定値は上記夫々の測定値に関連す
る光線が上記回転軸に垂直に延びる平面に対して成す角
度の余弦で重み付けされることを特徴とする、請求項１
記載のコンピュータ断層撮影方法。 - 【請求項４】 上記減衰値は、上記面上の画素の減衰値
からの正則３次元グリッド、望ましくは立方グリッドの
グリッド点において計算されることを特徴とする、請求
項１記載のコンピュータ断層撮影方法。 - 【請求項５】 上記測定値は各リンク値を計算するため
の１８０°／Ｎ１の角度範囲から取られ、但し、Ｎ１は
少なくとも値√Ｎθを有し、Ｎθは１８０°の角度範囲
に亘って分布され測定値が得られる放射線源の位置の数
である、請求項１記載のコンピュータ断層撮影方法。 - 【請求項６】 上記捕捉された測定値は、回転軸を含む
平面上での光線の方向及び位置によって決められるパラ
メータ空間中のデカルトグリッド上に位置する光線に関
連付けられる新しい測定値の組によって形成されるよう
リビニングされる、請求項１記載のコンピュータ断層撮
影方法。 - 【請求項７】 円錐状放射線ビームを発生する放射線源
と、 上記放射線源と上記検査ゾーンとの間に螺旋の形状の相
対運動を生じさせる駆動装置と、 上記検査ゾーン中の光線の減衰に依存する測定値の捕捉
のための検出器ユニットと、 上記検査ゾーン中の減衰値の空間的な分布を決める再構
成ユニットとを含み、 上記再構成ユニットは、 （ｄ）上記光線の位置及び向きを表わす３次元パラメー
タ空間中の網のリンクに沿って上記測定値を加算するこ
とによりリンク値を計算する段階と、 （ｅ）上記検査ゾーンの所与の面を通過する光線に関連
付けられるリンクについてのフィルタデータを生成する
ために上記リンク値をフィルタリングする段階と、 （ｆ）上記面上の画素における放射線の減衰を、上記画
素を通過する光線によって上記パラメータ空間中で決め
られる軌跡を近似するリンクのフィルタデータを加算す
ることにより計算する段階と、 （ｇ）上記回転軸の方向上互いにずらされた他の面につ
いて少なくとも上記段階（ｅ）及び（ｆ）を繰り返す段
階とを含む、請求項１記載の方法を行うためのコンピュ
ータ断層撮影装置。 - 【請求項８】 請求項７記載のコンピュータ断層撮影装
置の測定値を処理するためのコンピュータプログラムで
あって、 （ｄ）上記光線の位置及び向きを表わす３次元パラメー
タ空間中の網のリンクに沿って上記測定値を加算するこ
とによりリンク値を計算する段階と、 （ｅ）上記検査ゾーンの所与の面を通過する光線に関連
付けられるリンクについてのフィルタデータを生成する
ために上記リンク値をフィルタリングする段階と、 （ｆ）上記面上の画素における放射線の減衰を、上記画
素を通過する光線によって上記パラメータ空間中で決め
られる軌跡を近似するリンクのフィルタデータを加算す
ることにより計算する段階と、 （ｇ）上記回転軸の方向上互いにずらされた他の面につ
いて少なくとも上記段階（ｅ）及び（ｆ）を繰り返す段
階とを含む、請求項７記載のコンピュータ断層撮影装置
の測定値を処理するためのコンピュータプログラム。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| DE10061120A DE10061120A1 (de) | 2000-12-07 | 2000-12-07 | Computertomographie-Verfahren mit helixförmiger Relativbewegung |
| DE10061120.6 | 2000-12-07 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2002224099A true JP2002224099A (ja) | 2002-08-13 |
Family
ID=7666316
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001373109A Ceased JP2002224099A (ja) | 2000-12-07 | 2001-12-06 | 螺旋状の相対運動を含むコンピュータ断層撮影方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US6542572B2 (ja) |
| EP (1) | EP1213685B1 (ja) |
| JP (1) | JP2002224099A (ja) |
| DE (2) | DE10061120A1 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005296653A (ja) * | 2004-04-08 | 2005-10-27 | Siemens Ag | コンピュータ断層撮影装置によるコンピュータ断層撮影画像形成方法およびコンピュータ断層撮影装置 |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4090970B2 (ja) * | 2003-09-09 | 2008-05-28 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 放射線断層画像撮影装置と放射線断層画像撮影方法および画像生成装置と画像生成方法 |
| WO2005121836A1 (en) * | 2004-06-09 | 2005-12-22 | Philips Intellectual Property & Standards Gmbh | Computerized tomography method with helical relative movement and conical beam |
| US7573973B2 (en) * | 2005-05-17 | 2009-08-11 | General Electric Company | Methods and systems to facilitate reducing cone beam artifacts in images |
| JP4495109B2 (ja) * | 2006-04-06 | 2010-06-30 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｘ線ｃｔ装置 |
| US20090010518A1 (en) * | 2007-07-02 | 2009-01-08 | Harald Schoendube | Method for reconstruction of computed tomography representations from x-ray ct data sets of an examination subject with spiral scanning |
| GB2483831B (en) * | 2009-07-14 | 2015-07-22 | Rapiscan Systems Inc | Method for image reconstruction by using multi-sheet surface rebinning |
| US8811707B2 (en) | 2011-08-25 | 2014-08-19 | General Electric Company | System and method for distributed processing of tomographic images |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3637074B2 (ja) * | 1992-12-15 | 2005-04-06 | 株式会社東芝 | ヘリカルスキャン方式のコンピュータ断層撮影装置 |
| SE9602594D0 (sv) * | 1996-07-01 | 1996-07-01 | Stefan Nilsson | Förfarande och anordning vid datortomografi |
| US5802134A (en) * | 1997-04-09 | 1998-09-01 | Analogic Corporation | Nutating slice CT image reconstruction apparatus and method |
-
2000
- 2000-12-07 DE DE10061120A patent/DE10061120A1/de not_active Withdrawn
-
2001
- 2001-11-30 US US10/015,840 patent/US6542572B2/en not_active Expired - Fee Related
- 2001-12-05 DE DE50101837T patent/DE50101837D1/de not_active Expired - Fee Related
- 2001-12-05 EP EP01000712A patent/EP1213685B1/de not_active Expired - Lifetime
- 2001-12-06 JP JP2001373109A patent/JP2002224099A/ja not_active Ceased
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005296653A (ja) * | 2004-04-08 | 2005-10-27 | Siemens Ag | コンピュータ断層撮影装置によるコンピュータ断層撮影画像形成方法およびコンピュータ断層撮影装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| EP1213685A1 (de) | 2002-06-12 |
| EP1213685B1 (de) | 2004-03-31 |
| DE50101837D1 (de) | 2004-05-06 |
| US20020113215A1 (en) | 2002-08-22 |
| DE10061120A1 (de) | 2002-06-13 |
| US6542572B2 (en) | 2003-04-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8094910B2 (en) | Method of reconstructing an image function from Radon data | |
| US20040179643A1 (en) | Apparatus and method for reconstruction of volumetric images in a divergent scanning computed tomography system | |
| US6130930A (en) | Exact region of interest cone beam imaging without circle scans | |
| JP2000081318A (ja) | ３次元コンピュ―タトモグラフィイメ―ジングのためのスキャニングおよびデ―タ収集方法およびイメ―ジング装置 | |
| JP4342164B2 (ja) | コンピュータ断層撮影装置 | |
| CN1947154A (zh) | 使用截短的投影和在先采集的3d ct图像的锥形束ct设备 | |
| JP3787374B2 (ja) | 断層撮影画像内のアーチファクトを低減するための方法および装置 | |
| JP2000107168A (ja) | 被検体内の関心領域の３次元コンピュ―タ断層撮影方法及び装置 | |
| JP2005504571A (ja) | 多機能コーンビーム結像装置とその方法 | |
| JP2000107167A (ja) | 円錐状放射線ビ―ムを使用するコンピュ―タ断層撮影方法 | |
| JPH11326243A (ja) | Ｃｔイメ―ジング装置を動作する方法及び走査及びデ―タ収集装置 | |
| JP4813681B2 (ja) | コンピュータ断層撮影方法 | |
| JP2004188163A (ja) | 断層撮影装置 | |
| JPH10146331A (ja) | 対象物の関心ある領域の３次元コンピュータトモグラフ撮影のための走査及びデータ取得方法および装置 | |
| JP2000102530A (ja) | 被検体の中の関心領域の３次元コンピュ―タ断層画像形成のための方法および装置 | |
| JP4087547B2 (ja) | コンピュータ断層撮影装置 | |
| JP3290726B2 (ja) | 透過型三次元断層撮影装置 | |
| JP4347061B2 (ja) | 逐次コンピュータ断層撮影方法 | |
| US5748697A (en) | Method and apparatus for eliminating boundary errors in cone beam imaging | |
| US6333960B1 (en) | Exact region of interest cone beam imaging without circle scans | |
| JP3980696B2 (ja) | 画像再構成処理装置 | |
| CN1936958B (zh) | 用于从二维投影图像中重建三维图像体积的方法和装置 | |
| JP3788846B2 (ja) | Ｘ線コンピュータ断層撮影装置 | |
| JP4129572B2 (ja) | 傾斜三次元ｘ線ｃｔ画像の再構成方法 | |
| JP2002224099A (ja) | 螺旋状の相対運動を含むコンピュータ断層撮影方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20041202 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20061030 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080304 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080603 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080819 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20081118 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20090707 |
|
| A045 | Written measure of dismissal of application [lapsed due to lack of payment] |
Free format text: JAPANESE INTERMEDIATE CODE: A045 Effective date: 20091124 |