JP2008519634A - エネルギー分解コンピュータトモグラフィ - Google Patents
エネルギー分解コンピュータトモグラフィ Download PDFInfo
- Publication number
- JP2008519634A JP2008519634A JP2007540763A JP2007540763A JP2008519634A JP 2008519634 A JP2008519634 A JP 2008519634A JP 2007540763 A JP2007540763 A JP 2007540763A JP 2007540763 A JP2007540763 A JP 2007540763A JP 2008519634 A JP2008519634 A JP 2008519634A
- Authority
- JP
- Japan
- Prior art keywords
- energy
- detector
- distance
- interest
- electromagnetic 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.)
- Withdrawn
Links
- 238000002591 computed tomography Methods 0.000 title claims abstract description 39
- 238000005070 sampling Methods 0.000 claims abstract description 23
- 238000007689 inspection Methods 0.000 claims abstract description 10
- 238000012546 transfer Methods 0.000 claims abstract description 6
- 230000005855 radiation Effects 0.000 claims description 59
- 238000001514 detection method Methods 0.000 claims description 54
- 230000005670 electromagnetic radiation Effects 0.000 claims description 37
- 238000000034 method Methods 0.000 claims description 15
- 230000001427 coherent effect Effects 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 10
- 239000000463 material Substances 0.000 claims description 6
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 239000013056 hazardous product Substances 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000004154 testing of material Methods 0.000 description 1
- 230000003936 working memory Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
- G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
-
- 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/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4241—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
-
- 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/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- 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/48—Diagnostic techniques
- A61B6/483—Diagnostic techniques involving scattered radiation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Toxicology (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
コンピュータトモグラフィ手荷物検査または医療アプリケーションにおけるフィルタ後方投影再構成の前に運動量移動の補間をしても、必ずしも画像品質は最高にならず、計算コストも最低にならない。本発明の一実施形態では、エネルギー分解単一列ディテクタの非線形エネルギー・ビニングを提供する。これにより、自動的に、パラレル・ビニングされたディテクタのカーテシアンqサンプリングができる。これにより、フィルタ後方投影再構成の前のq補間を回避でき、空間解像度が高くなり、計算量が減少し、画像品質がよくなる。
Description
本発明は、医療用アプリケーション等におけるコンピュータトモグラフィの分野に関する。特に、本発明は、コンピュータトモグラフィ装置、放射線ディテクタ、コンピュータトモグラフィ装置における関心対象(object of interest)の検査方法、コンピュータトモグラフィ装置における関心対象の検査を実行するコンピュータプログラムに関する。
干渉性散乱(CS)コンピュータトモグラフィ(CT)(coherent scatter computed tomography)は、干渉性散乱X線光子に基づく新しい画像化方法である。干渉性散乱CTシステムは、対象の1つのスライスを照射するX線管と、検出システムとより構成され、両方とも患者またはその他の観測対象の周りを回転する。検出システムは、2次元ディテクタであってオフプレーン(off−plane)の散乱光子を測定するもの、または単一列ディテクタであって主ビームのプレーンから距離Hのところにあるものであって、これらは散乱光子のエネルギー分解測定を行う。測定した投影データから、主放射面の2つの空間次元(x,y)により画成される3次元ボリュームを再構成する。第3の次元は散乱光子の運動量移動qにより数値化される。
しかし、フィルタ後方投影再構成(filtered back−projection reconstruction)の前に、運動量移動値qを補間しなければならない。これにより画像品質が低下し、追加的な計算が必要となる。
関心対象の改良された検査が望まれている。
本発明の一実施形態によると、上記課題は、関心対象を検査するためのコンピュータトモグラフィ装置により解決される。該コンピュータトモグラフィ装置は、回転面中で回転し、関心対象に電磁放射ビームを照射する回転放射源と、第1のエネルギー区間内の第1のエネルギーを有する、前記関心対象の第1の対象点から干渉性散乱された電磁放射を検出するように構成された第1の検出要素と、第2のエネルギー区間内の第2のエネルギーを有する、前記関心対象の第2の対象点から干渉性散乱された電磁放射を検出するように構成された第2の検出要素と、前記第1の検出要素は前記回転平面から第1の距離に配置され、前記第2の検出要素は前記回転平面から第2の距離に配置され、前記第1の距離は前記第2の距離と基本的に同じであり、前記第1のエネルギー区間は前記第2のエネルギー区間と異なる。
言い換えると、放射源の回転平面に平行な線上に位置し、異なるエネルギー範囲に感度を有する第1の検出要素と第2の検出要素とを有するディテクタを提供する。それゆえ、第1の検出要素は、第2の検出要素とエネルギーが異なる散乱光子を検出できる。
有利にも、これにより、手荷物検査やカーテシアンqサンプリングができる医療アプリケーションのためのエネルギー分解干渉性散乱コンピュータトモグラフィ装置が得られる。
本発明の他の実施形態によると、第1の対象点と第2の対象点は、電磁放射のビームの中央光線と垂直な線上に配置される。
本発明の他の一実施形態によると、前記第1のエネルギー区間は、前記中央光線と、前記放射源から前記第1の対象点に放射された光線との間の第1のファン角度の所定の関数であり、前記第2のエネルギー区間は、前記中央光線と、前記放射源から前記第2の対象点に放射された光線との間の第2のファン角度の所定の関数である。
有利にも、第1と第2のファン角度の第1と第2のエネルギー区間の依存性は関数により予め決められるので、関心対象の検査が始まる前に、第1と第2の検出要素をそれぞれ第1と第2のエネルギー区間内の放射エネルギーに非常に感度よいように構成できる。
本発明の他の一実施形態によると、本コンピュータトモグラフィ装置はデータプロセッサをさらに有し、前記データプロセッサは、各検出要素のエネルギーを線形サンプリングし、ビームをパラレルビーム配置(parallel−beam geometry)にパラレルビーム・リビニング(rebinning)して、補間をしなくても、各検出要素の検出放射の運動量移動の等距離サンプリングとなる段階を実行するように構成されている。
このように、有利にも、パラレル・リビニングされたディテクタ平面上でqのカーテシアンサンプリングとなる。これにより、再構成ボリューム中の解像度がよくなる。
本発明の他の一実施形態によると、前記第1の検出要素と前記第2の検出要素は放射ディテクタの一部であり、前記放射ディテクタは、回転平面から一定の距離にあるフォーカス中心単一列エネルギー分解ディテクタと前記回転平面から一定の距離にある平面単一列エネルギー分解ディテクタとのうちの一方である。
有利にも、これにより、コンピュータトモグラフィ装置の計算効率が向上する。
本発明の他の一実施形態によると、前記電磁放射源は多色X線源であり、前記放射源は前記関心対象のまわりのらせん経路に沿って動き、前記ビームはファンビーム配置を有する。
多色X線源の使用は有利である。多色X線は容易に発生でき、光子束が大きいからである。
本発明の他の実施形態は、干渉性散乱コンピュータトモグラフィ装置として構成されたコンピュータトモグラフィ装置を提供する。
本コンピュータトモグラフィ装置は、手荷物検査装置、医療アプリケーション装置、材料テスト装置、及び材料科学分析装置よりなる群のうちの1つとして構成されてもよい。しかし、本発明を応用する最も好ましい分野は、手荷物検査と医療アプリケーションである。本発明により空間解像度が高くなり、計算量が減少し、画像品質が向上するからである。本発明により、一定タイプの材料を自動認識し、必要に応じて、危険物があると警告を発する高品質の自動化システムができる。かかる検査システムは、本発明のコンピュータトモグラフィ装置を使用し、X線源でX線を放射し、そのX線は検査する手荷物を透過、または手荷物で散乱する。それによりエネルギー分解して干渉性散乱された放射を検出できる。
本発明は、さらに放射ディテクタにも関する。該放射ディテクタは、第1のエネルギー区間内の第1のエネルギーを有する、電磁放射源から放射され、前記関心対象の第1の対象点から干渉性散乱された電磁放射を検出するように構成された第1の検出要素と、第2のエネルギー区間内の第2のエネルギーを有する、電磁放射源から放射され、前記関心対象の第2の対象点から干渉性散乱された電磁放射を検出するように構成された第2の検出要素と、第1の検出要素と第2の検出要素は、放射源の回転平面から等距離には位置され、第1のエネルギー区間は第2のエネルギー区間と異なる。
有利にも、これにより改良された放射ディテクタが提供され、空間解像度が向上し、画像品質がよくなる。
本発明の他の一実施形態によると、コンピュータトモグラフィにおける関心対象の検査方法が開示されている。該方法は以下の段階を有する:回転平面内で電磁放射源を回転する段階と、前記放射源から関心対象に電磁放射ビームを放射する段階と、第1の検出要素により、第1のエネルギー区間内の第1のエネルギーを有する、前記関心対象の第1の対象点から干渉性散乱された電磁放射を検出する段階と、第2の検出要素により、第2のエネルギー区間内の第2のエネルギーを有する、前記関心対象の第2の対象点から干渉性散乱された電磁放射を検出する段階。第1の検出要素は回転平面から第1の距離に配置され、第2の検出要素は回転平面から第2の距離に配置され、第1の距離は基本的に第2の距離と等しく、第1のエネルギー区間は第2のエネルギー区間と異なる。
本発明は、例えば、プロセッサ(画像プロセッサ等)上で実行されるコンピュータプログラムにも関する。このようなコンピュータプログラムは、例えば、CTスキャナシステムの一部である。本発明の一実施形態によるコンピュータプログラムは、好ましくはデータプロセッサのワーキングメモリにロードされる。このように、データプロセッサは本発明の方法の一実施形態例を実行する。本コンピュータプログラムは、C++等のどの好適なプログラミング言語で記述してもよく、CD−ROM等のコンピュータ読み取り可能媒体に格納してもよい。また、コンピュータプログラムは、ワールドワイドウェブ等のネットワークから入手することもできる。本コンピュータプログラムをネットワークから画像処理ユニットまたはプロセッサ、またはその他の好適なコンピュータにダウンロードしてもよい。
本発明の一態様は、エネルギー分解単一列ディテクタの非線形エネルギー・ビニングであり、これにより、自動的に、パラレル・ビニングされたディテクタのカーテシアンqサンプリングができる。これにより、フィルタ後方投影再構成の前のq補間を回避でき、空間解像度が高くなり、計算量が減少し、画像品質がよくなる。
本発明の上記その他の態様は、以下に説明する実施形態から明らかとなり、この実施形態を参照して詳しく説明される。
添付図面を参照して、発明の実施形態例を次に説明する。
以下、図1を参照して、エネルギー分解CSCTを実施するコンピュータトモグラフィ装置を説明する。別の図面においても、同様または同一の要素には同じ参照数字を付した。
この実施形態例を参照して、本発明を医療画像化に応用する場合について説明する。しかし、留意すべきことは、本発明は、医療が増加の分野における応用に限定されず、手荷物等に関する危険物(爆発物等)検出のための手荷物検査や、材料テスト等のその他の産業応用に応用することができることである。
図1に示したスキャナはファンビーム(fan−beam)CTスキャナである。図1に示したCTスキャナは、回転軸2の回りに回転可能なガントリ1を有する。ガントリ1はモータ3により駆動される。参照数字4は、X線光源等の放射源を示す。このX線光源は、本発明の一態様では、多色放射ビームを放射する。
参照数字5は、放射光源から放射された放射光ビームをコーン形状の放射光ビーム6に形成するアパチャーシステム(aperture system)を指す。コーン形状の放射ビーム6を放射した後、そのビームはスリットコリメータ(図1には図示せず)を通して導かれ、主ファンビームを形成して、対象領域にある対象7にあたる。
ファンビーム6は、(図1では誇張して示した;その経路で散乱されなければ、実際には検出要素の中心列だけにあたる、)ガントリ1の中心、すなわちCSCTスキャナの検査領域に配置された関心対象7を透過して、ディテクタ8にあたるように方向付けられる。図1から分かるように、ディテクタ8は、放射源4に対向するガントリ1上に配置されて、ディテクタ8の表面はファンビーム6でカバーされる。図1に示したディテクタ8は複数のディテクタ要素を有する。
関心対象7をスキャンする際、放射源4、アパチャーシステム5、及びディテクタ8は、矢印16で示された方向にガントリ1の回りを回転する。放射光源4、アパチャーシステム5、及びディテクタ8を有するガントリ1の回転のため、モータ3がモータ制御部17に接続されている。そのモータ制御部17は計算部18に接続されている。
スキャンの際、放射ディテクタ8は所定の時間間隔でサンプリングされる。放射ディテクタ8から読み出したサンプリング結果は電気信号、すなわち電気的データであり、以下の説明では投影(projection)と呼ぶ。それゆえ、関心対象のスキャンで得られたデータセットは、複数の投影により構成されている。投影の数は、放射ディテクタ8がサンプリングされる時間間隔に対応する。また、複数の投影を併せてボリューメトリックデータと呼ぶ。さらに、ボリューメトリックデータは心電図データも含み得る。
図1では、関心対象はコンベヤベルト19上に配置される。関心対象7のスキャン中、ガントリ1は患者7の回りを回転するが、コンベヤベルト19は関心対象7をガントリ1の回転軸2と平行な方向に沿って並進移動する。これにより、関心対象17はヘリカルスキャンパスに沿ってスキャンされる。コンベヤベルト19は、スキャン中は停止してもよい。コンベヤベルト19を設ける替わりに、例えば、医療への応用では、関心対象7が患者の場合には可動テーブルを使用する。しかし、留意すべきことは、説明する全ての場合において、回転スキャンを実行することが可能である。この場合、回転軸2に平行な方向での並進移動はなく、回転軸2に回りにガントリ1が回転するだけである。
ディテクタ8は計算部18に接続されている。計算部18は、検出結果、すなわち、ディテクタ8のディテクタ要素からのリードアウトを受け取り、そのリードアウトに基づいてスキャン結果を決定する。ディテクタ8のディテクタ要素は、関心対象7によりファンビーム6に生じた減衰、またはエネルギーが一定のエネルギーインターバルにある関心対象7の対象点から干渉性散乱されたX線のエネルギーと強さを測定するように構成されている。さらに、計算部18は、ガントリ1の動きをモータ3、20またはコンベヤベルト19と調和させるために、モータ制御部17と通信する。
計算部18は、ディテクタ8のリードアウトから、画像を再構成するように構成されている。計算部18により生成された画像は、インターフェイス22を介してディスプレイ(図1には示さず)に出力される。
データプロセッサにより実現される計算部18は、各ディテクタ要素のエネルギーのリニアサンプリング(linear sampling)を実行し、パラレルビーム配置(parallel−beam geometry)へのビームのパラレルビーム・リビニング(parallel−beam rebinning)を行う。本発明の一態様によると、これにより、補間されていない各ディテクタ要素の検出された放射の運動量移動における等距離サンプリング(equidistant sampling)となる。
さらに、図1から分かるように、計算部18はラウドスピーカ21に接続され、例えば自動的に警告を発するようになっている。
図2は、本発明の一実施形態によるCSCTの取得配置(acquisition geometry)を示す概略図である。図2に示したCSCT配置は、単一ラインのフォーカス中心(focus−centered)ディテクタ23を示す。このディテクタ23は、本発明の一実施形態によるCSCTコンピュータトモグラフィ装置の一部である。この装置において、本発明の一実施形態による方法が実行される。
主ビーム平面から一定距離27にあるエネルギー分解フォーカス中心単一列検出システム23を使用する場合、(中心が放射源4と回転中心33の中間にある)円弧50上の対象点241、242、243、244、245からの散乱角Θは、ファン角度βの非線形関数として変化する。
本発明によると、ファン角度βの関数としてエネルギーサンプリング範囲の上下限を調整することにより、この効果を補正することができる。これらの2つの値はディテクタ列の要素毎に異なるが、それらの間のエネルギーは線形にサンプリングできる。その後のパラレルビーム配置へのファンビーム・パラレルビーム・リビニング(fan−beam to parallel−beam rebinning)により、ディテクタ平面は長方形となり等距離のq個の列が含まれる。
この方法により、CSCT再構成の空間的解像度が改善される。1つの補間ステップがディテクタ依存のエネルギー・ビニング(binning)により行われるからである。また、計算効率も向上する。再構成の前処理において1つの補間ステップを実行しなくてよいからである。最後に、等距離のq個の列を含む長方形のディテクタ形状により、CSCT再構成の結果の画像品質が、標準的なフィルタ後方投影(standard filtered back−projection)アプローチと比較して、向上する。
干渉性散乱コンピュータトモグラフィのフィルタ後方投影再構成の基本的方法は、U. van Stevendaal、J. P. Schlomka、 A. Harding、 M. Grass著「フィルタ後方投影に基づく干渉性散乱コンピュータトモグラフィの再構成アルゴリズム(A reconstruction algorithm for coherent scatter computed tomography based on filtered back−projection)」(Med. Phys. 30 (9) (2003) pp. 2465−2474)に記載されている。この文献はここに参照援用する。
本発明は、主ビーム30の平面から一定の距離27にある単一列のディテクタを有するCSCTシステムを開示している。このシステムは、干渉性散乱のエネルギー分解測定を実行するものである。測定した投影データから、主放射面30の2つの空間次元(x,y)により画成される3次元ボリュームを再構成し、第3の次元は散乱光子の運動量移動qにより数値化される。検出要素36と37を含むエネルギー分解ディテクタ列23の非線形エネルギー・ビニング(binning)を実行する。これにより、パラレル・リビニングされた(parallel rebinned)ディテクタ31(図3参照)でカーテシアンqサンプリング(Cartesian q sampling)が行われる。これにより再構成の際のq補間(q interpolation)が回避される。
第1の放射光は、放射源4から第1の関心点244に向けて中央光線26と放射源4から点244への光線との間の第1のファン角度で放射され、第1の関心点244において第1の散乱角度で第1の検出要素36に散乱される。さらに、第2の放射光線は、放射源4から第2の対象点245に(第2のファン角度で)放射され、第2の対象点245から第2の検出要素37に向けて干渉性散乱される。
回転平面30から一定の距離27にあるエネルギー分解フォーカス中心単一列検出システムを使用する場合、対象点244、245からの散乱角度Θは、(区間[−β0;+β0]にある)それぞれのファン角度βにより、次式(1)に従って変化する:
ここで、GとSは、それぞれ放射源からディテクタまで、及び放射源から回転中心33までの距離である。運動量移動qは、次式(2)に従って、散乱角度と光子エネルギーとに関係する:
ここで、Eは光子のエネルギーであり、hとcはそれぞれプランク定数と光速度である。
回転平面30から一定の距離27にあるエネルギー分解フォーカス中心単一列検出システムを使用する場合、対象点244、245からの散乱角度Θは、(区間[−β0;+β0]にある)それぞれのファン角度βにより、次式(1)に従って変化する:
q0を一定距離Hに対して+/−β0で測定した散乱角度であるとする。EmaxとEminを最大エネルギーと最小エネルギーであるとする。これらはエネルギー分解ディテクタ23を用いて検出できる。最大ファン角度β0において、最大検出可能運動量移動は次式(3)である:
これは、その列のすべてのディテクタ要素の最大検出可能値である。
ファン角度がβc=0°である中央光線では、すべてのディテクタ要素の最小検出可能運動量移動は次式(4)である:
ファン角度βが変化するすべての異なるディテクタ要素位置のqにおいて一定のサンプリングを実現するために、EmaxとEminを次式(5)及び(6)に従ってβの関数として決めなければならない:
及び、
ディテクタ要素毎の最大及び最小サンプリングエネルギーと、その間のエネルギー範囲の線形サンプリングとを用いて、q方向で補間をしなくても、パラレル・リビニングされた(parallel−rebinned)ディテクタ平面上のqにおけるカーテシアンサンプリング(Cartesian sampling)をすることができる。これにより、再構成ボリュームの解像度が向上する。
平面単一ラインエネルギー分解ディテクタ(planar single−line energy−resolved detector)の場合、散乱角度Θとファン角度βの間の関係は、次式(7)である
図3は、本発明の一実施形態によるファンビーム・パラレルビーム・リビニング(fan−beam to parallel−beam rebinning)CSCTの取得配置(acquisition geometry)を示す概略図である。図3の左側の図は、回転軸29に沿った断面図であり、右側の図は、回転軸29に垂直な方向の断面図である。ディテクタ要素毎のエネルギー・ビニング(energy binning)を延長した放射源4、41、42、43、44により表したが、これらは可変なので、長方形ディテクタ31のq個の列は等距離になる。
電磁放射源4は平面30内で回転し、電磁放射ビームを対象点241、242、243、244、245に放射する。これらの対象点は、回転軸29と垂直な線上に配置されている。電磁放射は対象点により散乱され、線形ディテクタ31の一部である検出要素により検出される。各検出要素のエネルギーを線形サンプリングし、ビームをパラレルビーム配置(parallel−beam geometry)にパラレルビーム・リビニング(rebinning)すると、追加的に補間をしなくても、各検出要素の検出放射の運動量移動の等距離サンプリングとなる。
図4は、本発明による方法の一実施形態例を示すフローチャートである。本方法は、ステップS1において始まり、投影データセットが取得される。これは、例えば、好適なCSCTスキャナシステムを使用することにより、または記憶装置から投影データを読み出すことにより行われる。その後、ステップS2において、第1のエネルギー区間内のエネルギーを有する、関心対象の第1の対象点から干渉性散乱された電磁放射は、第1の検出要素で検出される。同時に、またはその前に、またはその後に、第2のエネルギー区間内のエネルギーを有する、関心対象の第2の対象点から干渉性散乱された電磁放射が、第2の検出要素で検出される。第1と第2の検出要素は、回転放射源の回転平面と、フォーカス中心単一列エネルギー分解ディテクタまたは平面単一列エネルギー分解ディテクタの一方から同じ距離に配置される。第1のエネルギー区間は前記第2のエネルギー区間と異なる。第1と第2のエネルギー区分は、それぞれのファン角度に対応し、中央光線と放射源からそれぞれの対象点に放射された光線との間のそれぞれのファン角度の所定の関数により決定される。
第1の対象点と第2の対象点は、電磁放射のビームの中央光線と垂直な線上に配置される。
されに別のステップでは、各ディテクタ要素のエネルギーを線形サンプリングし、ビームをパラレルビーム配置(parallel−beam geometry)にパラレルビーム・リビニング(rebinning)すると、補間をしなくても、各ディテクタ要素の検出された放射の運動量移動(q)の等距離サンプリングとなる。
図5は、本発明による方法の一実施形態例を実行する、本発明によるデータ処理装置の一実施形態例を示す図である。図5に示したデータ処理装置は、関心対象を示す画像を格納するメモリ152に接続された中央処理装置(CPU)または画像プロセッサ151を含む。データプロセッサ151は、複数の入出力ネットワークまたは診断装置(CSCT装置等)に接続されてもよい。データプロセッサは、さらに、データプロセッサ151で計算、または受け入れられた情報または画像を表示するコンピュータモニター等の表示装置154に接続されている。オペレータまたはユーザは、キーボード155及び/またはその他の出力装置(図5には図示せず)を介してデータプロセッサ151を操作する。
さらに、バスシステム153を介して、画像処理及び制御プロセッサ151を例えば動きモニターに接続することも可能である。この動きモニターは関心対象の動きを監視する。例えば、患者の肺を画像化する場合、動きセンサーは例えば呼気センサーである。心臓を画像化する場合、動きセンサーは心電計である。
留意すべきことは、「有する」という用語は他の要素やステップを排除するものではなく、「1つの」という用語は複数の場合を排除するものではないことである。また、請求項に記載されたいくつかの手段の機能を単一のプロセッサまたはシステムが満たすことができることである。また、異なる実施形態に関して説明した要素を組み合わせてもよい。
請求項中の参照符号は、その請求項の範囲を限定するものと解釈してはならないことにも留意すべきである。
Claims (15)
- 関心対象を検査するためのコンピュータトモグラフィ装置であって、
回転面中で回転し、関心対象に電磁放射ビームを照射する回転放射源と、
第1のエネルギー区間内の第1のエネルギーを有する、前記関心対象の第1の対象点から干渉性散乱された電磁放射を検出するように構成された第1の検出要素と、
第2のエネルギー区間内の第2のエネルギーを有する、前記関心対象の第2の対象点から干渉性散乱された電磁放射を検出するように構成された第2の検出要素と、
前記第1の検出要素は前記回転平面から第1の距離に配置され、前記第2の検出要素は前記回転平面から第2の距離に配置され、前記第1の距離は前記第2の距離と基本的に同じであり、
前記第1のエネルギー区間は前記第2のエネルギー区間と異なる装置。 - 第1の対象点と第2の対象点は電磁放射ビームの中央光線と垂直な線上に位置する、請求項1に記載のコンピュータトモグラフィ装置。
- 前記第1のエネルギー区間は、前記中央光線と、前記放射源から前記第1の対象点に放射された光線との間の第1のファン角度の所定の関数であり、
前記第2のエネルギー区間は、前記中央光線と、前記放射源から前記第2の対象点に放射された光線との間の第2のファン角度の所定の関数である、請求項1に記載のコンピュータトモグラフィ装置。 - データプロセッサをさらに有し、
前記データプロセッサは、
各検出要素のエネルギーの線形サンプリングと、
前記ビームのパラレルビーム配置へのパラレルビーム・リビニングを適用する段階であって、補間されていない各検出要素の検出された放射の運動量移動の等距離サンプリングをする段階とを実行するように構成された請求項1に記載のコンピュータトモグラフィ装置。 - 前記第1の検出要素と前記第2の検出要素は放射ディテクタの一部であり、
前記放射ディテクタは、回転平面から一定の距離にあるフォーカス中心単一列エネルギー分解ディテクタと前記回転平面から一定の距離にある平面単一列エネルギー分解ディテクタとのうちの一方である、請求項1に記載のコンピュータトモグラフィ装置。 - 前記電磁放射源は多色X線源であり、
前記放射源は前記関心対象のまわりのらせん経路に沿って動き、前記ビームはファンビーム配置を有する、請求項1に記載のコンピュータトモグラフィ装置。 - 干渉性散乱コンピュータトモグラフィ装置として構成されている、請求項1に記載のコンピュータトモグラフィ装置。
- 手荷物検査装置、医療アプリケーション装置、材料テスト装置、及び材料科学分析装置よりなる群のうちの1つとして構成されている、請求項1に記載のコンピュータトモグラフィ装置。
- 放射ディテクタであって、
第1のエネルギー区間内の第1のエネルギーを有する、電磁放射源から放射され、前記関心対象の第1の対象点から干渉性散乱された電磁放射を検出するように構成された第1の検出要素と、
第2のエネルギー区間内の第2のエネルギーを有する、電磁放射源から放射され、前記関心対象の第2の対象点から干渉性散乱された電磁放射を検出するように構成された第2の検出要素と、
前記第1の検出要素は前記回転平面から第1の距離に配置され、前記第2の検出要素は前記回転平面から第2の距離に配置され、前記第1の距離は前記第2の距離と基本的に同じであり、
前記第1のエネルギー区間は前記第2のエネルギー区間と異なる放射ディテクタ。 - 第1の対象点と第2の対象点は電磁放射ビームの中央光線と垂直な線上に位置する、
前記第1のエネルギー区間は、前記中央光線と、前記放射源から前記第1の対象点に放射された光線との間の第1のファン角度の所定の関数であり、
前記第2のエネルギー区間は、前記中央光線と、前記放射源から前記第2の対象点に放射された光線との間の第2のファン角度の所定の関数である、請求項9に記載の放射ディテクタ。 - 各ディテクタ要素のエネルギーの線形サンプリングと、前記ビームをパラレルビーム配置へのパラレルビーム・リビニングにより、補間をしなくても、各検出要素の検出放射の運動量移動の等距離サンプリングとなる、請求項9に記載の放射ディテクタ。
- 前記放射ディテクタは、回転平面から一定の距離にあるフォーカス中心単一列エネルギー分解ディテクタと、平面単一列エネルギー分解ディテクタとのうちの一方である、請求項9に記載の放射ディテクタ。
- コンピュータトモグラフィにおける関心対象の検査方法であって、
回転平面内で電磁放射源を回転する段階と、
前記放射源から関心対象に電磁放射ビームを放射する段階と、
第1の検出要素により、第1のエネルギー区間内の第1のエネルギーを有する、前記関心対象の第1の対象点から干渉性散乱された電磁放射を検出する段階と、
第2の検出要素により、第2のエネルギー区間内の第2のエネルギーを有する、前記関心対象の第2の対象点から干渉性散乱された電磁放射を検出する段階と、
前記第1の検出要素は前記回転平面から第1の距離に配置され、前記第2の検出要素は前記回転平面から第2の距離に配置され、前記第1の距離は前記第2の距離と基本的に同じであり、
前記第1のエネルギー区間は前記第2のエネルギー区間と異なる方法。 - 各ディテクタ要素のエネルギーを線形サンプリングする段階と、
前記ビームのパラレルビーム配置へのパラレルビーム・リビニングを適用する段階であって、補間されていない各ディテクタ要素の検出された放射の運動量移動の等距離サンプリングをする段階とをさらに有し、
第1の対象点と第2の対象点は電磁放射ビームの中央光線と垂直な線上に位置し、
前記第1のエネルギー区間は、前記中央光線と、前記放射源から前記第1の対象点に放射された光線との間の第1のファン角度の所定の関数であり、
前記第2のエネルギー区間は、前記中央光線と、前記放射源から前記第2の対象点に放射された光線との間の第2のファン角度の所定の関数である、請求項13に記載の方法。 - コンピュータトモグラフィ装置における関心対象の検査を実行するコンピュータプログラムであって、前記コンピュータプログラムは、プロセッサにより実行されると、前記プロセッサに以下の動作を実行させるコンピュータプログラム:
回転面中で回転し、関心対象に電磁放射ビームを照射する回転放射源によりデータセットをロードし、前記データセットは、
第1の検出要素により検出され、第1のエネルギー区間内の第1のエネルギーを有する、前記関心対象の第1の対象点から干渉性散乱された電磁放射に対応するデータと、
第2の検出要素により検出され、第2のエネルギー区間内の第2のエネルギーを有する、前記関心対象の第2の対象点から干渉性散乱された電磁放射に対応するデータとを有し、
前記第1の検出要素は前記回転平面から第1の距離に配置され、前記第2の検出要素は前記回転平面から第2の距離に配置され、前記第1の距離は前記第2の距離と基本的に同じであり、
前記第1のエネルギー区間は前記第2のエネルギー区間と異なるコンピュータプログラム。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| GB0424877A GB0424877D0 (en) | 2004-11-11 | 2004-11-11 | Energy resolved computer tomography |
| PCT/IB2005/053565 WO2006051443A1 (en) | 2004-11-11 | 2005-11-02 | Energy resolved computer tomography |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2008519634A true JP2008519634A (ja) | 2008-06-12 |
Family
ID=33523527
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2007540763A Withdrawn JP2008519634A (ja) | 2004-11-11 | 2005-11-02 | エネルギー分解コンピュータトモグラフィ |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US20090060124A1 (ja) |
| EP (1) | EP1814460A1 (ja) |
| JP (1) | JP2008519634A (ja) |
| CN (1) | CN101056583A (ja) |
| GB (1) | GB0424877D0 (ja) |
| WO (1) | WO2006051443A1 (ja) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0423707D0 (en) * | 2004-10-26 | 2004-11-24 | Koninkl Philips Electronics Nv | Computer tomography apparatus and method of examining an object of interest with a computer tomography apparatus |
| US7646850B2 (en) * | 2007-01-18 | 2010-01-12 | The Research Foundation Of State University Of New York | Wide-field, coherent scatter imaging for radiography using a divergent beam |
| US9086366B2 (en) * | 2012-02-15 | 2015-07-21 | L-3 Communications Security And Detection Systems, Inc. | Determining a material property based on scattered radiation |
| US9489752B2 (en) | 2012-11-21 | 2016-11-08 | General Electric Company | Ordered subsets with momentum for X-ray CT image reconstruction |
| US11058369B2 (en) | 2019-11-15 | 2021-07-13 | GE Precision Healthcare LLC | Systems and methods for coherent scatter imaging using a segmented photon-counting detector for computed tomography |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2145485C1 (ru) * | 1998-03-12 | 2000-02-20 | Кванта Вижн, Инк. | Ультрамалоугловая рентгеновская томография |
| WO2004066215A1 (en) * | 2003-01-21 | 2004-08-05 | Philips Intellectual Property & Standards Gmbh | Computed tomography method with coherent scattered rays, and computed tomograph |
| WO2004105610A1 (en) * | 2003-05-28 | 2004-12-09 | Philips Intellectual Property & Standards Gmbh | Fan-beam coherent-scatter computer tomography |
-
2004
- 2004-11-11 GB GB0424877A patent/GB0424877D0/en not_active Ceased
-
2005
- 2005-11-02 CN CNA2005800385709A patent/CN101056583A/zh active Pending
- 2005-11-02 US US11/718,718 patent/US20090060124A1/en not_active Abandoned
- 2005-11-02 JP JP2007540763A patent/JP2008519634A/ja not_active Withdrawn
- 2005-11-02 EP EP05805180A patent/EP1814460A1/en not_active Withdrawn
- 2005-11-02 WO PCT/IB2005/053565 patent/WO2006051443A1/en not_active Ceased
Also Published As
| Publication number | Publication date |
|---|---|
| US20090060124A1 (en) | 2009-03-05 |
| CN101056583A (zh) | 2007-10-17 |
| GB0424877D0 (en) | 2004-12-15 |
| EP1814460A1 (en) | 2007-08-08 |
| WO2006051443A1 (en) | 2006-05-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7453974B2 (en) | Beam-hardening and attenuation correction for coherent-scatter CT | |
| US7778383B2 (en) | Effective dual-energy x-ray attenuation measurement | |
| US20090304145A1 (en) | Computer tomography apparatus | |
| US7590215B2 (en) | Coherent-scatter computer tomograph | |
| US7580499B2 (en) | Coherent-scatter computed tomography | |
| US20070019782A1 (en) | Fan-beam coherent-scatter computed tomography | |
| JP2008519976A (ja) | コンピュータトモグラフィ装置及び関心対象の検査方法 | |
| WO2006109227A2 (en) | Energy distribution reconstruction in ct | |
| CN101218607B (zh) | 用于在扫描视场外的对象点的重建算法 | |
| JP2008519634A (ja) | エネルギー分解コンピュータトモグラフィ | |
| CN101501478A (zh) | 用于获取图像数据的系统和方法 | |
| US20080095304A1 (en) | Energy-Resolved Computer Tomography | |
| US20080205598A1 (en) | Coherent Scatter Computer Tomography Material Identification | |
| US20090245458A1 (en) | Extension of the q-range in csct |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20081031 |
|
| A761 | Written withdrawal of application |
Free format text: JAPANESE INTERMEDIATE CODE: A761 Effective date: 20100201 |
