JP5600742B2 - 投射のデータを処理するための装置及び方法 - Google Patents

投射のデータを処理するための装置及び方法 Download PDF

Info

Publication number
JP5600742B2
JP5600742B2 JP2012527436A JP2012527436A JP5600742B2 JP 5600742 B2 JP5600742 B2 JP 5600742B2 JP 2012527436 A JP2012527436 A JP 2012527436A JP 2012527436 A JP2012527436 A JP 2012527436A JP 5600742 B2 JP5600742 B2 JP 5600742B2
Authority
JP
Japan
Prior art keywords
value
projection
mismatch
image
acquired
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2012527436A
Other languages
English (en)
Other versions
JP2013503683A (ja
Inventor
ニールセン,ティム
ケーラー,トーマス
ヨーハネス ブレンデル,ベルンハルト
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
Koninklijke Philips Electronics NV
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips NV, Koninklijke Philips Electronics NV filed Critical Koninklijke Philips NV
Publication of JP2013503683A publication Critical patent/JP2013503683A/ja
Application granted granted Critical
Publication of JP5600742B2 publication Critical patent/JP5600742B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/412Dynamic

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • General Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Generation (AREA)
  • Processing Or Creating Images (AREA)
  • Image Processing (AREA)

Description

発明の分野
発明は、投射の値を処理する装置、投射の値を処理する方法、及び、獲得された投射の値を処理するための対応するコンピュータープログラムに関係する。発明は、さらに、投射の値を処理する装置を具備する関心のある領域をイメージングするためのイメージングシステムに関係する。
発明の背景
例えば、“Principles of Computerized Tomographic Imaging”,Kak,Slaney,IEEE Press,1999(非特許文献1)においては、コンピュータ断層撮影システムが、X線を発生させるためのX線源及び検出器を具備することは、開示されたことである。X線源によって発生させられたX線が、関心のある領域を通じて横断すると共に投射データを獲得するための検出器によって検出される一方で、X線源及び検出器は、関心のある領域のまわりで回転する。関心のある領域のイメージは、再構築のアルゴリズムを使用することによって再構築されるが、それは、例えば、獲得された投射のデータをフィルター処理すると共に逆投射する。再構築されたイメージは、一般には、イメージのアーチファクトを見せる。
発明の概要
本発明の目的は、獲得された投射の値を処理するための投射の値を処理する装置を提供することであるが、それは、イメージの質を改善するための手段を提供するが、それは、獲得された投射の値から再構築されたものである。本発明のさらなる目的は、対応する投射の値を処理する方法及びコンピュータープログラムを提供することである。本発明のさらなる目的は、関心のある領域をイメージングするための対応するイメージングシステムを提供することであるが、それは、投射の値を処理する装置を具備する。
本発明の第一の態様において、獲得された投射の値を処理するための投射の値を処理する装置は、提供されるが、それにおいて、投射の値を処理する装置は、
− 獲得された投射の値を提供するための獲得された投射の値を提供するユニット、
− 再構築の仮定の考慮の下で獲得された投射の値から第一のイメージを再構築するための再構築ユニット、
− 再構築の仮定の考慮の下で再構築された第一のイメージを通じた投射をシミュレートすることによってシミュレートされた投射の値を決定するためのシミュレートされた投射の値を決定するユニット、
− 獲得された投射の値について不整合の値を決定するための不整合を決定するユニット
:を具備するが、
それにおいて、不整合の値は、獲得された投射の値及びシミュレートされた投射の値を比較することによって、再構築の仮定を備えたそれぞれの獲得された投射の値の不整合の程度を示すものである。
獲得された投射の値からのイメージの再構築は、逆の問題である。獲得された投射の値を予想することを試みる数学的な問題は、与えられたものである。このモデルは、順方向のモデルとみなされることができると共に、シミュレートされた投射の値を決定するためのシミュレートされた投射の値を決定するユニットによって使用されることができる。このモデルは、関心のある領域におけるボクセルの値と同様の変数及び再構築の間に変更されるものではない定数に依存する。イメージを再構築する為には、ボクセルの値と同様の再構築の変数は、獲得された投射の値に基づいて、及び、定数の再構築のパラメータを使用することによって、決定される。定数の再構築のパラメータは、再構築の仮定である。このように、獲得された投射のデータからのイメージの再構築は、再構築の変数及び再構築の仮定に基づいたものであるが、それにおいて、再構築の変数は、再構築の工程の間に最適化されると共に、上に記載されたモデルの部分である他の何かは、再構築の仮定である。
例えば、獲得された投射の値が、コンピュータ断層撮影装置によって獲得されてしまったものであるとすれば、再構築の仮定は、獲得の幾何学である。別の再構築の仮定は、再構築される対象において散乱することの抑制、投射の値を獲得するために使用されてきたものである検出器におけるピクセル間のクロストーク無し、X線ビームの強度の変動の補正、ビーム硬化の補正、その他であることができる、即ち、例えば、X線ビームの強度の変動及び/又はビーム硬化が存在するものではないことは、仮定されることができる。これらの再構築の仮定は、獲得された投射の値が、単一エネルギーの吸収係数によって特徴付けられることができる対象を通じた線積分として表現される、モデルに至る。さらなる再構築の仮定は、投射の値の獲得の間に再構築される対象の時間的な変動が無いというものであることができる。
再構築の仮定が、十分に良好なものではないとすれば、即ち、それらが、現実のものにフィットするものではないとすれば、再構築されたイメージは、獲得された投射の値と整合するものではない。不整合を決定するユニットは、獲得された投射の値及び再構築の仮定の間におけるこの不整合の程度を示すものである不整合の値を決定することに適合させられたものである。これらの不整合の値は、獲得された投射の値からイメージを再構築することの質を改善するために使用されることができる。例えば、再構築の仮定は、変更されることができると共に、不整合の値は、変更された再構築の仮定に関して再度決定されることができる。不整合の値が、変更された再構築の仮定に関してより小さいものであるとすれば、変更された再構築の仮定は、獲得された投射の値及び再構築の仮定の間における不整合によって引き起こされた再構成されたイメージにおけるアーチファクトを低減するための獲得された投射の値からイメージを再構築するために使用されることができる。このように、不整合の値を決定することによって、獲得された投射の値に基づいて再構築されるものである、イメージの質を改善するための手段は、提供される。
獲得された投射の値は、優先的に、コンピューター断層撮影イメージングシステム又はX線C−アーム・システムの投射の値である。獲得された投射の値は、また、核のイメージングシステム、例えば、陽電子放出断層撮影イメージングシステム又は単一光子放出コンピューター断層撮影イメージングシステム、と同様の別のイメージングのモダリティーによって獲得された投射の値であることができる。
獲得された投射の値を提供するユニットは、優先的に、獲得された投射の値が格納される又は格納するユニット又は無線の若しくは有線のデータ接続を介して獲得された投射の値を受容するための及び受容された獲得された投射の値を提供するための獲得された投射の値を受容するユニットである。獲得された投射の値を提供するユニットが、格納するユニット又は獲得された投射の値を受容するユニットであるとすれば、投射の値を処理する装置は、優先的に、少なくとも再構築ユニット、シミュレートされた投射の値を決定するユニット、不整合を決定するユニット、及び獲得された投射の値を提供するユニットを含むコンピュータシステムで構成されたものである。
既に上に述べたように、再構築の仮定は、獲得の幾何学であることができる。獲得の幾何学は、例えば、投射の値が獲得される投射の角度、投射の値を発生させるために使用された放射の強度、投射の値を発生させるために使用された放射の発散の程度によって定義されるが、それにおいて、発散の程度は、錐体の角度及び扇の角度その他のものによってパラメータで表示されることができる。
再構築ユニットは、優先的に、対象の第一のイメージを再構築するためのフィルター処理された逆の投射を行うことに適合させられたものである。
シミュレートされた投射の値を決定するユニットは、優先的に、再構築の仮定に従ってシミュレートされた投射の値を決定するための再構築の仮定の考慮の下で再構築された第一のイメージを通じた前方の投射をシミュレートすること適合させられたものである。
優先的に、各々の獲得された投射の値について、シミュレートされた投射の値は、決定されると共に、各々の獲得された投射の値について、即ち、獲得された投射の値及び対応するシミュレートされた投射の値の各々の対について、不整合の値は、決定される。
不整合を決定するユニットが、獲得された投射の値及びシミュレートされた投射の値の間における差を決定するための差を決定するユニットを具備することは、好適なことであるが、それにおいて、不整合を決定するユニットは、決定された差に基づいて獲得された投射の値についての不整合の値を決定することに適合させられたものである。
差を決定するユニットが、差の高い周波数の成分が差の低い周波数の成分と比べてより多く抑制されたものであるように、差を決定することは、さらに好適なことである。優先的に、高い周波数の成分は、低い周波数の成分ではなく高い周波数の成分を抑制することによって低い周波数の成分と比べてより多く抑制されたものである。
再構築ユニットがイメージを再構築することを可能にするためには、それぞれの再構築のアルゴリズムによって使用されるものである全ての量は、離散的な値として提供されたものである必要がある、即ち、全ての量は、数の有限個のセットによって表されたものである。離散的な値のこの使用は、離散化の効果を生み出す、即ち、たとえ不整合が、ゼロであると思われるとしても、獲得された投射の値及び再構築の仮定の考慮に下でシミュレートされたものであるシミュレートされた投射の値は、離散化の誤差のために異なるものであると思われる。これらの離散化の誤差は、高い空間的な周波数における、即ち、イメージングされる必要があるものである関心のある領域内における構造的な情報を含有する空間的な周波数と比べてより高いものである空間的な周波数における、不整合の値に不利に影響を及ぼす。従って、差の低い周波数の成分が高い周波数の成分と比べてより多く表明されるような差の決定は、決定された差における、及びこのように、決定された不整合の値における、離散化の誤差を抑制することができる。これは、不整合の値の質を改善する。
用語“低い周波数の成分”及び“高い周波数の成分”は、空間的な周波数に関係する。
獲得された投射の値及びシミュレートされた投射の値の間における差は、一般に、ある範囲の空間的な周波数を具備する、特に、それらは、高い周波数の成分の空間的な周波数と比べてより低いものである空間的な周波数を有する低い周波数の成分を具備する。不整合を決定するユニットは、優先的に、高い周波数の成分が抑制されるように、差を決定することに、即ち、差の他の周波数の成分と比べてより大きいものである空間的な周波数を有する周波数の成分が抑制されるように、差を決定することに、適合させられたものである。
離散化の誤差が存在するものである周波数の範囲は、イメージを再構築すると共にシミュレートされた投射の値を決定するための投射をシミュレーションする間に使用された離散化に依存する。この使用された離散化が知られたものであるので、離散化の誤差の周波数の範囲は、決定されることができる。不整合を決定するユニットは、高い周波数の成分、即ち、離散化の誤差に関係付けされたものである決定された周波数の範囲内の周波数の成分が、抑制されるように、優先的に、差を決定することに適合させられたものである。
差を決定するユニットが、獲得された投射の値をローパスフィルター処理することに、並びに、相互のものからそれぞれのローパスフィルター処理された獲得された投射の値及びそれぞれの獲得された投射の値に対応するシミュレートされた投射の値を減算することによって差を決定するために適合させられたものであることは、さらに好適なことである。
離散化の効果は、シミュレートされた投射の値の高い周波数の成分の損失に至る。このように、獲得された投射の値及びシミュレートされた投射の値が、相互のものから減算されるとすれば、減算の結果は、高い周波数の成分における離散化の誤差を示すことになる。獲得された投射の値においてもまた相互のものから獲得された投射の値及びシミュレートされた投射の値を減算する前に、高い周波数の成分が抑制されるとすれば、減算の結果は、高い周波数の成分におけるより少ない離散化の誤差を有する。これは、減算の結果に基づいて決定された不定性の値の質を改善する。
差を決定するユニットが、減算の結果を発生させるために相互のものからそれぞれの獲得された投射の値及びそれぞれの獲得された投射の値に対応するシミュレートされた投射の値を減算することによって、及び、減算の結果をローパスフィルター処理することによって、差を決定することに適合させられたものであることは、さらに好適なことである。
差が、ローパスフィルター処理された減算の結果として決定されたものであるので、離散化の誤差によって引き起こされたものである高い周波数の成分は、ローパスフィルターによって抑制される又は取り除かれる。また、これは、不整合の値の質を改善する。
差を決定するユニットが、減算の結果を発生させるために相互のものからそれぞれの獲得された投射の値及びそれぞれの獲得された投射の値に対応するシミュレートされた投射の値を減算することによって、獲得された投射の値について決定された減算の結果から減算のイメージを再構築することによって、並びに、再構築された減算のイメージを通じた投射をシミュレートすることによって、差を決定することに適合させられたものであることは、さらに好適なことである。
イメージを再構築すると共にシミュレートされた投射の値を決定するための投射をシミュレートする間における離散化は、高い周波数の成分の抑制に至る。このように、減算の結果が、減算のイメージを再構築するために使用されるとすれば、及び、再構築された減算のイメージを通じた投射が、シミュレートされるとすれば、シミュレーションの結果は、高い周波数の成分が、抑制された、特に、取り除かれたものであると共に、それによって減算の結果に基づいて決定されたものである不整合の値における離散化の誤差を抑制する、特に、取り除く、減算の結果である。その上、離散化の誤差は、第一のイメージを再構築すると共に再構築された第一のイメージを通じた投射をシミュレートすることによって投射の値を決定することのプロセスによって発生させられてしまってあるものであるが、それにおいて、高い周波数の成分の抑制、即ち、離散化の誤差の抑制は、同じプロセスによって行われる。従って、高い周波数の成分の抑制は、減算の結果の高い周波数の成分における離散化の誤差に非常に良好に適合させられたものである。これは、さらに、不整合の値の質を改善する。
不整合を決定するユニットが、さらに、獲得された投射の値についてノイズの値を決定するための及び決定されたノイズ値の逆元で決定された差に重みをつけるための、ノイズに重みをつけるユニットを具備することは、さらに好適なことであるが、それにおいて、不整合を決定するユニットは、重みがつけられた差に基づいて獲得された投射の値についての不整合の値を決定することに適合させられたものである。逆元のノイズの値で決定された差に重みをつけることによって、差は、ノイズの値に関して規格化される。これは、さらに、不整合の値の質、特に、異なる獲得された投射の値に属する不整合の値の適合性を改善する。
不整合を決定するユニットが、重みがつけられた差を二乗すること、及び、二乗された重みがつけられた差に基づいて不整合の値を決定ことに適合させられたものであることは、さらに好適なことである。
獲得された投射の値が、獲得の時間に割り当てられてしまってあるものであることは、さらに好適なことであるが、それにおいて、投射の値を処理するユニットは、さらに、それぞれの獲得の時間の獲得された投射の値について決定された不整合の値に獲得の時間を割り当てることによって異なる時間における獲得された投射の値及び再構築の仮定の間における不整合の示すものである時間的な不整合の信号を決定するための時間的な不整合の信号を決定するユニットを具備する。
再構築ユニットが、再構築の仮定として静止の対象を考慮することに適合させられたものであるとすれば、及び、対象が、移動するとすれば、時間的な不整合の信号は、対象の運動を示す時間的な運動の信号とみなされることができる。
優先的に、ユニットを決定する時間的不整合信号は、一次元の時間的不整合信号を生成するために同じ獲得時間に帰属している不整合の値を平均するのに適している。
再構築ユニットが獲得された投射の値及び決定された時間的不整合信号に基づいて第二のイメージを再構築するのに適していることは、更に好まれる。これは、獲得された投射の値及び再構築の仮定間の不整合によって生じるイメージ人工品を還元していた第二のイメージを再構築することを可能にする。
例えば、再構築ユニットは獲得された投射の値だけに基づいてイメージを再構築するのに適していることがありえる。そして、それは時間的不整合信号が閾値以下で不整合を示す時間に得られた。再構築ユニットはまた、第二のイメージを再構築するための獲得された投射の値を使用するのに適していることがありえる。そして、それは時間的不整合信号が所定の不整合範囲の中で不整合を示す時間に得られた。これは、獲得された投射に基づいてイメージを再構築することが類似した不整合を有することを評価すると認める。
投射の値処理ユニットが不整合イメージを出すためのユニットを決定する不整合イメージから更に成ることは更に好まれる。そこにおいて、不整合イメージのイメージ要素のために、それぞれのイメージ要素に貢献する不整合の値は加えられる。不整合イメージを出すために使われる不整合の値は優先的に、獲得された投射の値及びシミュレートされた投射の値を各々から減算することによって決定された差である。そこにおいて、減算の結果から、高い周波数の成分が、特に、取り除かれて抑制されるように、差は決定される。
この種の不整合イメージが、例えば、再構成イメージのどの一部が診断対象のために使われてはならないか指し示すための放射線科医によって使われることができる。
不整合イメージのそれぞれの要素に貢献する不整合の値の和は、フィルター処理されてない背景映写によって優先的に実行される。
例示的な実施形態において、獲得された投射データにフィルターをかけるための適切な低域フィルタ又は不整合が構造(離散化誤差によって生じることは公知である)を示している不整合イメージを評価すると決定するために使われることができる。この種の不整合イメージは、シミュレーションで測定されることができる。この例示的な実施形態において、構造(離散化誤差によって生じることは公知である)が減少するように、低域フィルターは選択されることができる。
それぞれの不整合の値が決定された獲得された投射の値が特定のイメージ要素によって表される対象の範囲内で位置を通り過ぎている光線によって生じる場合、不整合の値は特定のイメージ要素に貢献する。
再構築ユニットが獲得された投射の値及び獲得された投射の値及び再構築の仮定間の不整合によって生じるイメージ人工品を還元するための決定された不整合イメージに基づいて対象の第二のイメージを再構築するのに適していることは、更に好まれる。
投射の値処理ユニットが再構築された第一のイメージ、そして、不整合イメージを覆うための覆っているユニットから更に成ることが、更に好まれる。
これは、直接、ユーザに第一のイメージのどの一部が獲得された投射の値の不整合によって生じるイメージ人工産物に影響を受けるか示すことを可能にする。不整合イメージは、また、前述の第二のイメージのような他のイメージを通じて覆われるために用いることができる。
再構築ユニットが獲得された投射の値及び獲得された投射の値及び再構築の仮定間の不整合によって生じるイメージ人工品を還元するための決定された不整合の値から対象の第二のイメージを再構築するのに適していることは、更に好まれる。
例えば、決定された不整合の値が、特に、最小化されて減少するように、再構築の仮定は変更されることができる。そこにおいて、最終イメージは変更された再構築の仮定で考慮中に獲得された投射の値から再構築される。これは、再構築された最終イメージの質を改善する。
第二のイメージを再構築すると共に、再構築ユニットが所定の閾値以下で不整合の値を有する獲得された投射の値を決定して、所定の閾値以下で不整合の値を有する決定された獲得された投射の値だけを使用するのに適していることは、更に好まれる。
大きい不整合を示している獲得された投射の値がイメージを再構築するために使われないので、再構成イメージの質は改善されることができる。
再構築ユニットが類似した不整合の値を有する獲得された投射の値を決定して、類似した不整合の値を有する獲得された投射の値から対象のイメージを再構築するのに適していることは、更に好まれる。
類似した不整合の値は、対象の類似した様相を示すことができる。このように、類似した不整合の値を有する獲得された投射の値から対象のイメージを再構築することによって、対象のイメージは再構築されることができる。そして、対象の特定の様相に対応する。例えば、再構築の仮定が対象が移動していないということである場合、類似した不整合の値は類似した可動状態を示すことができる。この場合、類似した不整合の値を有する獲得された投射の値から対象のイメージを再構築することによって、特定の可動状態の対象のイメージは、再構築されることができる。
本発明の更なる態様において、興味がある領域を撮像するためのイメージングシステムは示される。そこにおいて、イメージングシステムは成り立つ:興味がある領域の投射の値を得るための投射の値獲得ユニット、そして、投射は、処理ユニットを請求項1に記載の処理獲得された投射の値のために高く評価する。本発明の更なる態様において、処理獲得された投射の値のための投射の値処理方法は提供される。そこにおいて、投射の値処理方法は次の工程を含む:獲得された投射の値を提供すること、獲得された投射の値から再構築の仮定で考慮中に第一のイメージを再構築すること、再構築の仮定で考慮中に再構築された第一のイメージによる投射をシミュレーションすることによってシミュレートされた投射の値を決定すること、獲得された投射の値及びシミュレートされた投射の値を比較することによる獲得された投射の値(不整合の値が再構築の仮定を有するそれぞれの獲得された投射の値のある程度の不整合を表す)のための決定不整合の値。
本発明の更なる態様において、興味がある領域を撮像するイメージング方法は示される。そこにおいて、イメージング方法は次の工程を含む:興味がある領域の投射の値を得ること、そして、処理は、請求項14に記載の投射値を得た。
本発明の更なる態様において、処理獲得された投射の値のためのコンピュータプログラムは提供される。そこにおいて、コンピュータプログラムはプログラムを含む投射を引き起こすためのコード手段は請求項14に記載の投射の値処理方法のステップを実行するために請求項1に記載の処理ユニットを評価する。そのとき、コンピュータプログラムはユニットを処理している投射の値を制御しているコンピュータに動く。
本発明の更なる態様において、興味がある領域を撮像するためのコンピュータプログラムは提供される。そこにおいて、コンピュータプログラムは請求項13に記載のイメージングシステムにイメージング方法のステップを実行させるためのプログラムコード手段を含む。そのとき、コンピュータプログラムはイメージングシステムを制御しているコンピュータに動く。
請求項1、請求項13のイメージングシステム、請求項14の投射の値処理方法、イメージング方法、請求項15の獲得された投射の値を得るためのコンピュータプログラム及び興味がある領域を撮像するためのコンピュータプログラムのユニットを処理している投射の値が従属クレームに記載の類似した及び/又は同一の好ましい実施形態を有すると理解される。
発明の好適な実施形態が、また、それぞれの独立なクレームとの従属のクレームのいずれの組み合わせでもあることができることは、理解されるものとする。
発明のこれらの及び他の態様は、以後に記載された実施形態から明らかなものであると共にそれらを参照して解明されることになる。
図1は、概略的に及び例示的に投射の値を処理する装置を具備する関心のある領域をイメージングするためのイメージングシステムを示す。 図2は、例示的に関心のある領域をイメージングするためのイメージングの方法を図解するフローチャートを示す。 図3は、概略的に及び例示的に獲得された投射の値のシノグラムを示す。 図4は、概略的に及び例示的に再構築された第一のイメージを示すが、それは、獲得された投射の値を使用することによって再構築されてしまってある。 図5は、概略的に及び例示的に不整合のシノグラムを形成する不整合の値を示す。 図6は、概略的に及び例示的に動きの信号である時間的な不整合の信号を示す。 図7は、概略的に及び例示的に不整合のイメージを示す。
図面の手短な記述
後に続く図面において:
図1は、概略的に及び例示的に投射の値を処理する装置を具備する関心のある領域をイメージングするためのイメージングシステムを示す、
図2は、例示的に関心のある領域をイメージングするためのイメージングの方法を図解するフローチャートを示す、
図3は、概略的に及び例示的に獲得された投射の値のシノグラム(sinogram)を示す、
図4は、概略的に及び例示的に再構築された第一のイメージを示すが、それは、獲得された投射の値を使用することによって再構築されてしまってある、
図5は、概略的に及び例示的に不整合のシノグラムを形成する不整合の値を示す、
図6は、概略的に及び例示的に動きの信号である時間的な不整合の信号を示す、並びに、
図7は、概略的に及び例示的に不整合のイメージを示す。
実施形態の詳細な説明
図1は、図式的に、そして、手本となって、コンピューター断層撮影システムであって、興味がある領域を撮像するためのイメージングシステムを示す。コンピューター断層撮影システムは、z方向と平行して伸びる回転軸Rについて回転ができる構台1を含む。この実施形態(X線管)の放射源2(ある)は、構台1に載置する。放射源2はコリメータ3を備えている。そして、それは、本実施形態において、円錐放射ビーム4を放射源2によって発生する放射から形成する。放射はオブジェクト(図示せず)(例えば患者)及び検査ゾーンに対する関心(対象の範囲内で優先的に位置する)の領域を横断する。そして、それは、本実施形態において、円筒状である。検査ゾーン5を横断した後に、放射ビーム4は検出ユニット6上の事件である。そして、それは二次元の検出面から成る。検出ユニット6は、構台1に載置する。
コンピューター断層撮影システムは、2台のモーター7、8を備えている。構台1は、モーター7によって好ましくは恒常的であるが、可調角速度でドライブされる。モーター8は対象(例えば患者)を排気して提供される。そして、その人は、回転軸R又はz軸の方向と平行して、検査ゾーン5の患者テーブルに配置される。例えば、これらのモーター7、8は制御ユニット9によって制御される。そうすると、放射源及び検査は螺旋形の軌道に沿って相対的に各々への検査ゾーン移動の範囲内で5及び、このように、興味がある領域をゾーンに分ける。しかしながら、放射源が対象又は検査ゾーン5と関連して円形の軌道を進めることは、また、もし放射源2だけが回転しなければ、対象又は検査ゾーン5がすなわち移動しない候補者である。さらにまた、もう一つの実施形態では、コリメータ3は他の光線形状(特に扇形ビーム)を形成することに適していることがありえる、そして、検出ユニット6は検出面から成ることができる。そして、それは、扇形ビームに特に、他の光線形状に対応して成形される。
放射源及び検査ゾーン5の相対的な変化の間、検出ユニット6は、検出ユニット6の検出面上の放射事件に従い、投射の値を得る。従って、事項、モーター7、8及び構台1において、放射源2(検査ゾーンと関連して放射源2を移動するための要素)及び検出ユニット6は、興味がある領域の獲得された投射の値を生み出すための投射の値獲得ユニットを形成する。
得された投射の値は、獲得された投射の値を処理する投射の値処理ユニット1に提供される。ユニット1を処理している投射の値には、獲得された投射の値を受信するための、そして、獲得された投射の値を提供するための受信ユニット12が設けられて、再構築ユニット13ように構成されている。興味のある領域は、検査ゾーンの中に位置して、対象の対象又は一部を含む。再構築ユニットは、再構築の仮定で考慮中に獲得された投射の値から対象の第一のイメージを再構築する。本実施形態における、再構築ユニット13が使用することによりぶりのイメージを再構築するのに適していることをフィルターをかける逆投射。再構築の仮定は、例えば、扇形ビーム及び焦点中心の探知器を有する円形のソース軌道、散乱防止格子による散乱放射の抑制、X線ビーム強度バリエーションの変更及び効果を堅くしている光線上でないデータ収集である。測定された投射が評価するその手段は、単一エネルギーのX線吸収係数によって記載されていることができる対象による線積分としてモデル化されることができる。
更なるありうる仮定は、以下の通りである対象は再構築可能な視野に完全に含まれる−対象がデータ収集の間、提議しない−ノイズ統計値がガウス分布によってよく近くされることができるために、X線強度は十分に高い。
二次元のフィルタ処理逆投射は、これらの仮定のための適切な再構築アルゴリズムである。
更なるユニット10を処理している投射の値は、再構築の仮定で考慮中に再構築された第一のイメージによる投射をシミュレーションすることによってシミュレートされた投射の値を決定するためのユニット14を決定するシミュレートされた投射の値を備えている。本実施形態において、シミュレートされた投射は、ユニットが再構築の仮定で考慮中に再構築された第一のイメージによる前方の投射をシミュレーションするのに適していると決定することを評価する。各々の獲得された投射の値のために、シミュレートされた投射の値は、決定される。
投射は、10が更に成る処理ユニットを評価する獲得された投射の値及びシミュレートされた投射の値を比較することによって獲得された投射の値(不整合の値が再構築の仮定を有するそれぞれの獲得された投射の値のある程度の不整合を表す)のための不整合の値を決定するための不整合を決定するユニット15。獲得された投射の値及び対応するシミュレートされた投射の値の各ペアのために、不整合の値は、決定される。
ユニット15を決定する不整合は、獲得された投射の値及びシミュレートされた投射の値の差を決定するためのユニット16を決定する差と、不整合を決定するユニット15がどこで、獲得された投射の値の差に基づいて獲得された投射の値のための不整合の値を決定するのに適しているかと、シミュレートされた投射の値とを備えている。さらに、高い周波数の成分が抑制されるように、ユニット16を決定する差は差を決定することに適している。高い周波数の成分を抑制するために、本実施形態において、ユニット16を決定する差は獲得された投射が評価する低域フィルタに適している、そして、それぞれの低パスを減算することによって差を決定することは各々から獲得された投射の値及び対応するシミュレートされた投射の値にフィルターをかけた。他の実施態様において、差を決定するユニットは、減算の結果を生成するためにそれぞれの獲得された投射の値及び対応するシミュレートされた投射の値を各々から減算することによって、そして、減算の結果にフィルターをかけている低パスによって差を測定するのに適していることがありえる。差を決定するユニットが減算の結果を生成するためにそれぞれの獲得された投射の値及び対応するシミュレートされた投射の値を各々から減算することによって、獲得された投射の値のために決定される減算の結果から減算イメージを再構築することによって、そして、再構築された減算イメージによる投射をシミュレーションすることによって差を測定するのに適していることは、また、可能である。
不整合を決定するユニット15は、獲得された投射の値のためのノイズ値を決定するための、そして、決定されたノイズ値の逆との決定された差を重くするためのノイズに重みをつけるユニット17を備えている。ノイズに重みをつけるユニット17は全ての差を重くするのに適している。そして、それは、それぞれの獲得された投射の値の逆のノイズ値によって、獲得された投射の値のために決定された。
ノイズ加重は、較正データに基づいて探知器ピクセルの検出強度から、優先的に決定される。つまり参照X線源及び反復測定を使用して、システムがそうである既知の事実のためのノイズ特徴は、決定した実験的にX線強度の範囲が治った。
不整合を決定するユニット15は、加重差を正方形にして、二乗加重差として不整合の値を決定するのに更に適している。このように、各々の獲得された投射の値のために、二乗加重差は、不整合の値として決定される。他の実施態様において、正方形の動作及び任意にまた、加重動作は、省略されることができる。それはユニット15が不整合の値として抑制された高い周波数の成分との決定された差を決定するのに適していることがありえると決定する不整合を意味する、又は、不整合の値として抑制された高い周波数の成分を有する加重差を決定するか又は二乗を決定することは不整合の値として抑制された高い周波数の成分との差を重くした。
獲得された投射の値は、獲得時間に割り当てられる。投射の値処理ユニット10は、直説法の時間的不整合信号を決定するためのユニット18を決定する時間的不整合信号から成る獲得された投射の値及び異なる時間の再構築の仮定間の不整合(獲得にそれぞれの獲得時間の獲得された投射の値のために決定される不整合の値の時間を割り当てることによる)。本実施形態において、ユニット18を決定する時間的不整合信号は、一次元の時間的不整合信号を生成するために同じ獲得時間に帰属している不整合の値を平均するのに適している。再構築の仮定がそれである場合、検査ゾーンの対象は移動しない、そして、対象が移動する場合、時間的不整合信号は対象の運動を示している時間的運動信号と考えられることができる。従って、不整合の値が、対象の運動を決定するために使われることができる。
再構築ユニット13は、獲得された投射の値及び決定された時間的不整合信号に基づいて第二のイメージを再構築するのに適していることがありえる。これは、そばに引き起こされるイメージ人工品を還元していた第二のイメージを再構築することを可能にする獲得された投射の値及び再構築の仮定間の不整合。本実施形態において、再構築の仮定(実質的に不整合の値に影響する)は、興味がある領域の対象が移動していないということである。再構築ユニットは従って、第二のイメージを再構築するための獲得された投射の値を使用するのに適している。そして、それは、特に、時間的不整合信号が時間的不整合信号が類似した不整合を示す所定の不整合範囲の中で不整合を示す時間に得られた。本実施形態において、不整合の値が主に物(非移動物体の再構築の仮定と整合していない)の動きによって影響されるので、類似した不整合の値は対象の類似した可動様相を示す。再構築ユニットは、従って、再構築される対象の同じ可動様相に対応する獲得された投射の値に基づいてイメージを再構築するのに適している。
再構築ユニット13はまた、獲得された投射の値だけに基づいて第二のイメージを再構築するのに適していることがありえる。そして、それは時間的不整合信号が閾値以下で不整合を示す時間に得られた。これは大きい不整合(すなわち閾値より上の不整合の値)が画質に減少する影響を及ぼすことを確実にする。それによって、第二のイメージの質を改善する。
時間的不整合信号を使用せずに、再構築ユニット13は、また、獲得された投射の値及び獲得された投射の値及び再構築の仮定間の不整合によって生じるイメージ人工品を還元するための決定された不整合の値から対象の第二のイメージを再構築するのに適していることがありえる。例えば、第一のイメージは、変更された再構築によって再び再構築されることができる同じことに基づく仮定は、投射の値を得た。決定された不整合の値が、特に、最小化されて減少するように、再構築の仮定のこの変更は実行されることができる。そこにおいて、最終イメージは変更された再構築の仮定で考慮中に獲得された投射の値から再構築される。このように、不整合の値が最小化されるまで、再構築の仮定は変更されることができる、そして、これらの変更された再構築の仮定が獲得された投射の値からイメージを再構築するために使われることができる。そして、それは改良された画質を有する。
例えば、変更された再構築の仮定は光線を堅くしている変更を含むことができる、すなわち、堅くなっている光線があると仮定される、そして、再構築アルゴリズムは対応する光線を堅くしている変更を含む。他のありうる変更された再構築の仮定は、投射の値の第一の一部を得ると共に、対象が移動するということで、そして、投射の値の第二部品を得ると共に、対象が移動しないということでありえた。優先的にそれから変更された再構築の仮定を考慮している再構築アルゴリズムは、対象のイメージを再構築するための獲得された投射の値の第二部品を使用するだけである。
時間的不整合信号の上記した用途と同様で、第二のイメージを再構築すると共に、再構築ユニット13はまた、所定の閾値以下で不整合の値を有する獲得された投射の値を決定して、所定の閾値以下で不整合の値を有する獲得された投射の値だけを使用するのに適していることがありえる。さらに、また、上記に対する類似物は、側頭部を使う能力を記載した不整合は信号を送る、再構築ユニット13は類似した不整合の値を有する獲得された投射の値を決定するのに適していることがありえる、そして、獲得された投射から対象のイメージを再構築することが、類似した不整合の値を有することを評価する。
投射の値処理ユニット10は不整合イメージを出すためのユニット19を決定する不整合イメージから更に成る。そこにおいて、不整合イメージのイメージ要素のために、それぞれのイメージ要素に貢献する不整合の値は加えられる。不整合イメージのこの生成は、フィルター処理されてない逆投射によって実行される。このフィルター処理されてない逆投射は、優先的に使用する第一のイメージを再構築するための再構築ユニット13によって使われた再構築の仮定である。
ユニット19が不整合に基づいて不整合イメージを出すのに適していることがありえると決定する不整合イメージは減少する高い周波数の成分を有する減算の結果であるか又は不整合の値に基づいて減少する高い周波数の成分を有するノイズ加重減算の結果であることを評価する、又は、二乗されている不整合の値に基づいて、ノイズは抑制された高周波を有する減算の結果を重くした構成要素。
再構築ユニット13は、獲得された投射の値及び獲得された投射の値及び再構築の仮定間の不整合によって生じるイメージ人工品を還元するための決定された不整合イメージに基づいて対象の第二のイメージを再構築するのに適していることがありえる。例えば、不整合イメージに示される不整合が減少するように、再構築奪取は変更されることができる、特に、不整合イメージによって示される不整合が最小化されるように、再構築奪取は最適化されることができる。興味がある領域のイメージは、それから最適化された再構築の仮定を用いて獲得された投射の値から再構築されることができる。
更なるユニット10を処理している投射の値は、再構築された第一のイメージを覆うための覆っているユニット20と不整合イメージを備えている。これは、直接、ユーザに第一のイメージのどの一部が再構築の仮定に関して獲得された投射の値の不整合によって生じるイメージ人工産物に影響を受けるか示すことを可能にする。不整合イメージは、また、再構築された第二のイメージのような他のイメージを通じて覆われるために用いることができる。
投射の値処理ユニット10は優先的に不整合シノグラムを形成する不整合の値を示すためのディスプレイユニット11、第一のイメージを通じて又は第二のイメージのような他の再構成イメージを通じて覆われる不整合イメージ及び変更イメージ(すなわち第二のイメージ)から更に成る。そして、それは、事項、時間的不整合信号及び/又は不整合イメージで、不整合の値を用いて再構築された。
また、投射の値処理ユニットは、制御ユニット9によって優先的に制御される。
以下において、興味がある領域を撮像するイメージング方法は、図2に示されるフローチャートに関して記載されている。
興味がある領域を撮像するイメージング方法は、本実施形態において、コンピューター断層撮影イメージング方法である。ステップ101において、放射源2は回転軸Rのまわりを回転する、そして、対象又は検査ゾーン5は移動しない(すなわち2が対象又は検査ゾーン5周辺に円形の軌道に沿って進行する放射源)。他の実施形態では、放射源2は、対象と関連して他の軌道(例えば螺旋形の軌道)を進めることができる。放射源2は、少なくとも興味がある領域の対象を横断している放射を発する。放射(対象を横断した)は検出ユニット6によって検出される。
そして、それは投射データである測定されたデータを生成する。獲得された投射の値はユニット10を処理している投射の値へ移される、そして、ステップ102で、提供された獲得された投射の値が考慮中に再構築の仮定で興味がある領域の対象の第一のイメージを再構築するために使われる。ステップ103において、シミュレートされた投射の値は、再構築の仮定で考慮中に再構築された第一のイメージによる投射をシミュレーションすることによって決定される。ステップ104において、獲得された投射の値はフィルターをかけられる低パスである、そして、ステップ105で、低パスは獲得された投射の値にフィルターをかけた、そして、シミュレートされた投射の値は各々から減算される。減算の結果(高い周波数の成分が抑制された)は、不整合の値とすでに誤差を考えることができる。
ステップ106において、各々の獲得された投射の値のために、すなわち各々の不整合の値のために、ノイズ値は決定される、そして、ステップ107で、不整合の値は逆のそれぞれのノイズ値によって傾く。ステップ108において、加重不整合の値は、二乗される。
ステップ109において、時間的不整合信号は、獲得にそれぞれの獲得時間の獲得された投射の値のために決定される不整合の値の時間(獲得された投射の値が得られた)を割り当てることによって決定される。同じ獲得時間に帰属している不整合の値は、一次元の時間的不整合信号を生成するために平均値になる。同じ獲得時間に帰属している不整合の値を平均する代わりに、不整合の値は、また、単に合計されることができる。ステップ110において、第二のイメージは、獲得された投射の値及び獲得された投射データ及び再構築の仮定間の不整合によって生じるイメージ人工品を還元することへの決定された時間的不整合動機に基づいて再構築される。ステップ111において、第二のイメージは、ユーザにディスプレイユニットに表示される。
ステップ102〜111は、投射の値処理方法のステップと考えることができる。投射の値において、処理方法ステップ109〜111は省略されることができる、すなわち、投射の値処理方法は不整合の値を決定するだけであるのに適していることがありえる。また、ノイズの決定、ノイズ加重及び、正方形にすることに関連があるステップ106〜108は、省略されることができる。さらに、投射の値処理方法は、更に詳細で上記の通りに他の方法の不整合価値を使用するための他のステップを含むことができる。
図3は、図式的に、そして、手本となって、獲得された投射の値のシノグラム25を示す。図3において、水平軸は時間tを示す、そして、垂直軸はz’番目がそれぞれの時間tで計量される投射の値を得たことを示す。図3に示されるシノグラムは運動によって損なわれる、すなわち、人工産物は獲得された投射の値(関心上の領域の対象の運動によって影響される)間の不整合を原因として生じるのでシノグラムにおいて誘発される、そして、合図する再構築の仮定は興味がある領域に存在しない。X線が獲得された投射の値を生み出すために横断した対象は、本実施形態において、中心の小さい球体が移動している胸部模型である。この運動は小さいだけである、そして、対応する腐敗は従って、図3に示されるシノグラムにおいて見えない。
図4は、図式的に、そして、手本となっては、再構築された第一のイメージ29(図3に示される獲得された投射の値を用いて再構築された)である。胸部21の範囲内の中心的素子22は移動している。それによって、中心的素子22の近くで再構築された第一のイメージの運動人工産物を引き起こす。
図5は、図式的に、そして、手本となって、獲得された投射の値のために決定された、そして、不整合シノグラム26を形成する不整合の値を示す。図5に示される不整合は、主に中心的素子22の運動によって生じる。このように、図3に示される獲得された投射の値の損なわれたシノグラムが運動によって腐敗を示さないにもかかわらず、結果として生じる不整合(すなわち運動)は不整合シノグラムにおいて検出されることができる。従って、再構築の仮定が運動が再構築される関心の領域に存在しないということである場合、不整合の値が興味がある領域の運動を決定するために使われることができる。第23行は、不整合シノグラムの一部でなくて、ちょうど中心的素子22の運動の期間を示す。それとして図5において見る、不整合の値24の期間は、運動23の期間に対応する。獲得された投射値の構造は、図5に示される不整合シノグラムにおいて、もはや見えない。
図6は、図式的に、そして、手本となって、時間的不整合信号を示す/(t)、それは、図5に示される不整合の値を合計することによって決定された。そして、それは、任意のユニットで、同じ獲得時間に対応する。このように、一次元の時間的不整合信号は、図5に示される不整合シノグラムに基づいて発生した。時間的不整合信号27は、非常によく実際の運動期間と相関して、背景比率に非常に良好な信号を示す。
図7は、図式的に、そして、手本となって、大きい不整合(すなわち低い信頼)が白として描写される不整合イメージ28を示す。低い信頼の領域は、図4に示される再構築された第一のイメージの運動人工産物の位置と一致する。図7に示される不整合イメージは、図5に示される不整合シノグラムのフィルター処理されてない背景映写によって発生した。
投射の値のためのイメージ再構築アルゴリズムは、測定されたデータが静止対象による線積分に対応すると共通に仮定する。例えば、この仮定が運動、対象の範囲内の時間的吸収変化又は金属のために成し遂げられない場合、再構成イメージは人工産物を示す。ユニット及び方法を処理している獲得された投射の値は、獲得された投射の値が再構築アルゴリズムの再構築の仮定と整合しているかどうか分析するために、獲得された投射の値に適用されることができる。再構築の仮定が再構築される対象が移動しないということである場合、獲得された投射の値(すなわち不整合の値)の検出矛盾は運動信号(位置運動がどの時間に起こったかについて示す)を生成するために用いることがありえる、そして、また、獲得された投射の値の不整合に影響を受ける再構成イメージのそれらのボクセルを示す信頼イメージと考えることができる不整合イメージを作成する、いずれがイメージ人工産物の位置を示すかすなわち。獲得された投射の値の不整合に由来する情報が、人工産物を最小化しようとする次のイメージ再構築において使われることができる。
イメージの中で分かれる徴候が診断対象のために使われてはならないように、不整合イメージは放射科医のようなユーザに提出されることができる。それが1ボクセルの強度が正しいという確率を示すように、不整合イメージは構成されることができる。すでに信頼マップより上に又はとして言及されるように、不整合イメージは従って、また、信頼イメージと考えることができる。
優先的に、各々の獲得された投射の値のために、不整合の程度は計算される、すなわち、不整合の値は決定される。不整合の程度は、不整合がどの時間位置で起こったかについて示す時間的信号を生成するために用いることができる。これは、患者でありえる対象の同じ運動様相に帰属する獲得された投射の値の一部を識別するために役立つ。例えば、予想外の呼吸は検出されることができる、又は、サイクルがありえる心臓病患者の位相はいずれが人間の核心の又は動物のイメージを再構築するために使われることができるかについて検出した。
このボクセルのためのカイ2乗値を計算するために、不整合イメージが不整合が評価する、重くなって、正方形にされるノイズであった、そして、所与のボクセルに貢献する全てを合計することによって発生する場合、不整合イメージは基本的に再構成イメージの各々のボクセルのためのカイ2乗検定である。すなわちボクセルが獲得された投射の値及び再構築の仮定間の不整合によって遂行されないために、カイ2乗分布の周知の特性を使用して、それぞれのボクセルの結果として生じる値はこのボクセルの値が正しいという確率を計算するために用いることができる。カイ2乗分布の決定のために、ノイズがガウス分布に続くと優先的に仮定される。計算された確率は、それぞれのボクセルが不整合に影響を受けないという確率を示している変更された不整合イメージを形成することができる。
上記した実施形態においてイメージングシステムがコンピューター断層撮影システムであるにもかかわらず、他の実施形態では、獲得された投射の値は陽電子放出X線断層撮影システム又は単一光子放出コンピューター断層撮影システムのようなC−アーム又は核イメージングシステムのような他のイメージング様相によって得られることができる。
特定の再構築の仮定が使われた、そして、特定の再構築の仮定が変更されたにもかかわらず、上記した実施形態において、上述した実施形態で、更に、他の実施形態でも他の再構築の仮定がイメージを再構築するために使われることができて、変更されることができる。
図面、開示及び添付の請求の範囲の研究から、開示された実施形態に対する他の変化は、理解されることができて、請求された本発明を実施する際の当業者によってもたらされることができる。
請求項において、「成り立っている」語が他の要素又はステップを除外しないこと、そして、不定冠詞「ある」多数を除外する。
単一のユニット又はユニットは、請求項において詳述されるいくつかの部材の機能を成し遂げることができる。特定の計測が相互に異なる従属クレームにおいて詳述されるという単なる事実は、これらの計測の組合せが有効に使われることができないことを示さない。
算出は再構築(シミュレーション又は不整合の値の測定)が1時までに実行されるのを好む、又は、いくつかのユニット又はユニットは他のいかなる数ものユニット又はユニットによって実行されることができる。例えば、ステップ102〜110は、単一ユニットによって又は他のいかなる数ものユニットによって実行されることができる。算出及び/又は判定及び/又はイメージング方法に従うイメージングシステム及び/又は投射の値に従う投射の値処理ユニットの制御の規制処理方法は、コンピュータプログラムのプログラムコード手段として及び/又は専用ハードウエアとして行うことができる。
プログラムがそうすることができるコンピュータは、他のハードウェアと共に又はの一部として供給されて、適切な媒体(例えば光記憶媒体又は固体媒体)に格納されていて/、分散であるまた、インターネット又は他の有線であるかワイヤレス電気通信システムを経てのような、他の形で配布されることができる。
請求項のいかなる参照符号も、範囲を制限するものとして解釈されてはならない。
本発明は処理獲得された投射の値のための投射の値処理ユニットに関する。第一のイメージは、再構築ユニットによって再構築の仮定で考慮中に獲得された投射の値から再構築される。シミュレートされた投射はユニットが再構築奪取で考慮中に再構築された第一のイメージによる投射をシミュレーションすることによってシミュレートされた投射の値を決定すると決定することを評価する、そして、不整合の値は不整合を決定するユニットによって獲得された投射の値のために決定される。そこにおいて、獲得された投射の値及びシミュレートされた投射の値を比較することによって、不整合の値は再構築の仮定を有するそれぞれの獲得された投射の値のある程度の不整合を表す。例えば、不整合の値が異なる対象のために使われることができる。そして、再構成イメージの質を改善するか又は再構築の仮定及び獲得された投射の値間の不整合に影響を受けているイメージ要素を示す。
クレーム:
1. 獲得された投射の値を処理するための投射の値を処理する装置であって、
前記投射の値を処理する装置(10)は、
− 獲得された投射の値を提供するための獲得された投射の値を提供するユニット(12)、
− 再構築の仮定の考慮の下で前記獲得された投射の値から第一のイメージ(29)を再構築するための再構築ユニット(13)、
− 前記再構築の仮定の考慮の下で前記再構築された第一のイメージを通じた投射をシミュレートすることによってシミュレートされた投射の値を決定するためのシミュレートされた投射の値を決定するユニット(14)、
− 前記獲得された投射の値について不整合の値を決定するための不整合を決定するユニット(15)
:を具備する、装置において、
不整合の値は、前記獲得された投射の値及び前記シミュレートされた投射の値を比較することによって、前記再構築の仮定を備えたそれぞれの獲得された投射の値の不整合の程度を示すものである、装置。
2. クレーム1に定義されたような投射の値を処理する装置において、
前記不整合を決定するユニット(15)は、前記獲得された投射の値及び前記シミュレートされた投射の値の間における差を決定するための差を決定するユニット(16)を具備すると共に、
前記不整合を決定するユニットは、前記決定された差に基づいて前記獲得された投射の値についての不整合の値を決定することに適合させられたものである、
装置。
3. クレーム2に定義されたような投射の値を処理する装置において、
前記差を決定するユニット(16)は、前記差の高い周波数の成分が、前記差の低い周波数の成分と比べてより抑制されたものであるように、前記差を決定することに適合させられたものである、装置。
4. クレーム3に定義されたような投射の値を処理する装置において、
前記差を決定するユニット(16)は、前記獲得された投射の値をローパスフィルター処理することに並びに相互のものから前記それぞれのローパスフィルター処理された獲得された投射の値及び前記それぞれの獲得された投射の値に対応する前記シミュレートされた投射の値を減算することによって前記差を決定することに適合させられたものである、装置。
5. クレーム3に定義されたような投射の値を処理する装置において、
前記差を決定するユニットは、減算の結果を発生させるために相互のものから前記それぞれの獲得された投射の値及び前記それぞれの獲得された投射の値に対応する前記シミュレートされた投射の値を減算することによって並びに前記減算の結果をローパスフィルター処理することによって前記差を決定することに適合させられたものである、装置。
6. クレーム3に定義されたような投射の値を処理する装置において、
前記差を決定するユニットは、減算の結果を発生させるために相互のものから前記それぞれの獲得された投射の値及び前記それぞれの獲得された投射の値に対応する前記シミュレートされた投射の値を減算することによって、前記獲得された投射の値について決定された前記減算の結果から減算のイメージを再構築することによって、並びに、前記再構築された減算のイメージを通じた投射をシミュレートすることによって、前記差を決定するために適合させられたものである、装置。
7. クレーム3に定義されたような投射の値を処理する装置において、
前記不整合を決定するユニット(15)は、さらに、獲得された投射の値についてのノイズの値を決定するための及び前記決定されたノイズの値の逆元で前記決定された差に重みをつけるためのノイズに重みをつけるユニット(17)を具備すると共に、
前記不整合を決定するユニット(15)は、前記重みがつけられた差に基づいて獲得された投射の値について前記不整合の値を決定することに適合させられたものである、
装置。
8. クレーム7に定義されたような投射の値を処理する装置において、
前記不整合を決定するユニット(15)は、前記重みがつけられた差を二乗すること及び前記二乗された重みがつけられた差に基づいて不整合の値を決定することに適合させられたものである、装置。
9. クレーム1に定義されたような投射の値を処理する装置において、
前記獲得された投射の値は、獲得の時間へ割り当てられてしまってあると共に、
前記投射の値を処理する装置(10)は、さらに、前記それぞれの獲得の時間の前記獲得された投射の値について決定された前記不整合の値へ前記獲得の時間を割り当てることによって時間的な不整合の信号を決定するユニット(18)を具備する、
装置。
10. クレーム1に定義されたような投射の値を処理する装置において、
前記投射の値を処理する装置は、さらに、不整合のイメージ(28)を発生させるための不整合のイメージを決定するユニット(19)を具備すると共に、
前記不整合のイメージ(28)のイメージの要素について、前記それぞれのイメージの要素に寄与する不整合の値は、加算される、
装置。
11. クレーム10に定義されたような投射の値を処理する装置において、
前記投射の値を処理する装置(10)は、さらに、前記再構築された第一のイメージ及び前記不整合のイメージをオーバーレイするためのオーバーレイするユニット(20)を具備する、装置。
12. クレーム1に定義されたような投射の値を処理する装置において、
前記再構築ユニット(13)は、前記獲得された投射の値及び前記再構築の仮定の間における不整合によって引き起こされたイメージのアーチファクトを低減するために前記獲得された投射の値及び前記決定された不整合の値から対象の第二のイメージを再構築することに適合させられたものである、装置。
13. 関心のある領域をイメージングするためのイメージングシステムであって、
前記イメージングシステム(1)は、
− 前記関心のある領域の投射の値を獲得するための投射の値の獲得ユニット、及び、
− クレーム1に定義されたような獲得された投射の値を処理するための投射の値を処理する装置(10)
:を具備する、イメージングシステム。
14. 獲得された投射の値を処理するための投射の値を処理する方法であって、
前記投射の値を処理する方法は、後に続くステップ:
− 獲得された投射の値を提供すること、
− 再構築の仮定の考慮の下で前記獲得された投射の値から第一のイメージを再構築すること、
− 前記再構築の仮定の考慮の下で前記再構築された第一のイメージを通じた投射をシミュレートすることによってシミュレートされた投射の値を決定すること、
− 前記獲得された投射の値について不整合の値を決定すること
を具備する、方法において、
不整合の値は、前記獲得された投射の値及び前記シミュレートされた投射の値を比較することによって、前記再構築の仮定を備えたそれぞれの獲得された投射の値の不整合の程度を示すものである、方法。
15. 獲得された投射の値を処理するためのコンピュータープログラムであって、
前記コンピュータープログラムは、クレーム1に定義されたような投射の値を処理する装置にクレーム14に定義されたような投射の値を処理する方法のステップを、前記コンピュータープログラムが、前記投射の値を処理する装置を制御するコンピューターにおいて走らせられるとき、実行することを引き起こすためのプログラムコード手段を具備する、コンピュータープログラム。

Claims (15)

  1. 獲得された投射の値を処理するための投射の値を処理する装置であって、
    前記投射の値を処理する装置は、
    獲得された投射の値を提供するための獲得された投射の値を提供するユニット、
    再構築の仮定の考慮の下で前記獲得された投射の値から第一のイメージを再構築するための再構築ユニット、
    前記再構築の仮定の考慮の下で前記再構築された第一のイメージを通じた投射をシミュレートすることによってシミュレートされた投射の値を決定するためのシミュレートされた投射の値を決定するユニット、
    前記獲得された投射の値について不整合の値を決定するための不整合を決定するユニット
    を具備する、装置において、
    不整合の値は、前記獲得された投射の値及び前記シミュレートされた投射の値を比較することによって、前記再構築の仮定を備えたそれぞれの獲得された投射の値の不整合の程度を示すものであると共に、
    前記不整合を決定するユニットは、前記獲得された投射の値及び前記シミュレートされた投射の値の間における差を決定し、前記差の高周波成分が、前記差の低周波成分と比べてより抑制されたものであるようにする、差を決定するユニットを具備すると共に、
    前記不整合を決定するユニットは、前記決定された差に基づいて前記獲得された投射の値についての不整合の値を決定することに適合させられたものである、
    装置。
  2. 請求項1に記載の投射の値を処理する装置において、
    前記差を決定するユニットは、前記獲得された投射の値をローパスフィルター処理することに並びに相互のものから前記それぞれのローパスフィルター処理された獲得された投射の値及び前記それぞれの獲得された投射の値に対応する前記シミュレートされた投射の値を減算することによって前記差を決定することに適合させられたものである、装置。
  3. 請求項1に記載の投射の値を処理する装置において、
    前記差を決定するユニットは、減算の結果を発生させるために相互のものから前記それぞれの獲得された投射の値及び前記それぞれの獲得された投射の値に対応する前記シミュレートされた投射の値を減算することによって並びに前記減算の結果をローパスフィルター処理することによって前記差を決定することに適合させられたものである、装置。
  4. 請求項1に記載の投射の値を処理する装置において、
    前記差を決定するユニットは、減算の結果を発生させるために相互のものから前記それぞれの獲得された投射の値及び前記それぞれの獲得された投射の値に対応する前記シミュレートされた投射の値を減算することによって、前記獲得された投射の値について決定された前記減算の結果から減算のイメージを再構築することによって、並びに、前記再構築された減算のイメージを通じた投射をシミュレートすることによって、前記差を決定するために適合させられたものである、装置。
  5. 請求項1に記載の投射の値を処理する装置において、
    前記不整合を決定するユニットは、さらに、獲得された投射の値についてのノイズの値を決定するための及び前記決定されたノイズの値の逆元で前記決定された差に重みをつけるためのノイズに重みをつけるユニットを具備すると共に、
    前記不整合を決定するユニットは、前記重みがつけられた差に基づいて獲得された投射の値について前記不整合の値を決定することに適合させられたものである、
    装置。
  6. 請求項5に記載の投射の値を処理する装置において、
    前記不整合を決定するユニットは、前記重みがつけられた差を二乗すること及び前記二乗された重みがつけられた差に基づいて不整合の値を決定することに適合させられたものである、装置。
  7. 請求項1に記載の投射の値を処理する装置において、
    前記獲得された投射の値は、獲得の時間へ割り当てられてしまってあると共に、
    前記投射の値を処理する装置は、さらに、前記それぞれの獲得の時間の前記獲得された投射の値について決定された前記不整合の値へ前記獲得の時間を割り当てることによって時間的な不整合の信号を決定するユニットを具備する、
    装置。
  8. 請求項1に記載の投射の値を処理する装置において、
    前記投射の値を処理する装置は、さらに、不整合のイメージを発生させるための不整合のイメージを決定するユニットを具備すると共に、
    前記不整合のイメージのイメージの要素について、前記それぞれのイメージの要素に寄与する不整合の値は、加算される、
    装置。
  9. 請求項8に記載の投射の値を処理する装置において、
    前記投射の値を処理する装置は、さらに、前記再構築された第一のイメージ及び前記不整合のイメージをオーバーレイするためのオーバーレイするユニットを具備する、装置。
  10. 請求項1に記載の投射の値を処理する装置において、
    前記再構築ユニットは、前記獲得された投射の値及び前記再構築の仮定の間における不整合によって引き起こされたイメージのアーチファクトを低減するために前記獲得された投射の値及び前記決定された不整合の値から対象の第二のイメージを再構築することに適合させられたものである、装置。
  11. 関心のある領域をイメージングするためのイメージングシステムであって、
    前記イメージングシステムは、
    前記関心のある領域の投射の値を獲得するための投射の値の獲得ユニット、及び、
    請求項1に記載の獲得された投射の値を処理するための投射の値を処理する装置
    を具備する、イメージングシステム。
  12. 獲得された投射の値を処理するための投射の値を処理する方法であって、
    前記投射の値を処理する方法は、後に続くステップ、
    獲得された投射の値を提供すること、
    再構築の仮定の考慮の下で前記獲得された投射の値から第一のイメージを再構築すること、
    前記再構築の仮定の考慮の下で前記再構築された第一のイメージを通じた投射をシミュレートすることによってシミュレートされた投射の値を決定すること、
    前記獲得された投射の値について不整合の値を決定すること
    を具備する、方法において、
    不整合の値は、前記獲得された投射の値及び前記シミュレートされた投射の値を比較することによって、前記再構築の仮定を備えたそれぞれの獲得された投射の値の不整合の程度を示すものであると共に、
    前記獲得された投射の値及び前記シミュレートされた投射の値の間における差は、前記差の高い周波数の成分が、前記差の低い周波数の成分と比べてより抑制されたものであるように、決定されると共に、
    前記獲得された投射の値についての不整合の値は、前記決定された差に基づいて決定される、
    方法。
  13. 関心のある領域をイメージングするためのイメージング方法において、
    前記イメージング方法は、後に続くステップ、
    前記関心のある領域の投射の値を獲得すること、及び、
    請求項12に定義されたような獲得された投射の値を処理すること
    を具備する、イメージング方法。
  14. 獲得された投射の値を処理するためのコンピュータープログラムであって、
    前記コンピュータープログラムは、請求項1に記載の投射の値を処理する装置に請求項12に記載の投射の値を処理する方法のステップを、前記コンピュータープログラムが、前記投射の値を処理する装置を制御するコンピューターにおいて走らせられるとき、実行することを引き起こすためのプログラムコード手段を具備する、コンピュータープログラム。
  15. 関心のある領域をイメージングするためのコンピュータープログラムにおいて、
    前記コンピュータープログラムは、前記コンピュータープログラムが前記イメージングシステムを制御するコンピューターにおいて走らせられるとき、請求項11に定義されたようなイメージングシステムが請求項13に定義されたようなイメージング方法のステップを実行することを引き起こすためのプログラムコード手段を具備する、コンピュータープログラム。
JP2012527436A 2009-09-07 2010-09-01 投射のデータを処理するための装置及び方法 Expired - Fee Related JP5600742B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP09169583 2009-09-07
EP09169583.3 2009-09-07
PCT/IB2010/053916 WO2011027298A1 (en) 2009-09-07 2010-09-01 Apparatus and method for processing projection data

Publications (2)

Publication Number Publication Date
JP2013503683A JP2013503683A (ja) 2013-02-04
JP5600742B2 true JP5600742B2 (ja) 2014-10-01

Family

ID=43033903

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012527436A Expired - Fee Related JP5600742B2 (ja) 2009-09-07 2010-09-01 投射のデータを処理するための装置及び方法

Country Status (6)

Country Link
US (1) US8666137B2 (ja)
EP (1) EP2476099B1 (ja)
JP (1) JP5600742B2 (ja)
CN (1) CN102483853B (ja)
RU (1) RU2541860C2 (ja)
WO (1) WO2011027298A1 (ja)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2546804A1 (en) * 2011-07-10 2013-01-16 Dürr Dental AG Method and tomography apparatus for reconstruction of a 3D volume
DE102012204019B4 (de) * 2012-03-14 2018-02-08 Siemens Healthcare Gmbh Verfahren zur Reduzierung von Bewegungsartefakten
EP3175790B1 (en) * 2013-11-04 2021-09-08 Ecential Robotics Method for reconstructing a 3d image from 2d x-ray images
EP3123446B1 (en) * 2014-03-26 2020-07-08 Koninklijke Philips N.V. Image generation apparatus
US9986983B2 (en) 2014-10-31 2018-06-05 Covidien Lp Computed tomography enhanced fluoroscopic system, device, and method of utilizing the same
EP3317853B1 (en) * 2015-06-30 2019-08-07 Koninklijke Philips N.V. Registration apparatus for registering images
US10674982B2 (en) 2015-08-06 2020-06-09 Covidien Lp System and method for local three dimensional volume reconstruction using a standard fluoroscope
US10702226B2 (en) 2015-08-06 2020-07-07 Covidien Lp System and method for local three dimensional volume reconstruction using a standard fluoroscope
US10716525B2 (en) 2015-08-06 2020-07-21 Covidien Lp System and method for navigating to target and performing procedure on target utilizing fluoroscopic-based local three dimensional volume reconstruction
US11172895B2 (en) 2015-12-07 2021-11-16 Covidien Lp Visualization, navigation, and planning with electromagnetic navigation bronchoscopy and cone beam computed tomography integrated
US11051886B2 (en) 2016-09-27 2021-07-06 Covidien Lp Systems and methods for performing a surgical navigation procedure
CN106840109A (zh) * 2017-01-20 2017-06-13 长沙全度影像科技有限公司 一种基于差分的单像素成像硬件系统
US10699448B2 (en) 2017-06-29 2020-06-30 Covidien Lp System and method for identifying, marking and navigating to a target using real time two dimensional fluoroscopic data
US10893843B2 (en) 2017-10-10 2021-01-19 Covidien Lp System and method for identifying and marking a target in a fluoroscopic three-dimensional reconstruction
US10905498B2 (en) 2018-02-08 2021-02-02 Covidien Lp System and method for catheter detection in fluoroscopic images and updating displayed position of catheter

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5416815A (en) * 1993-07-02 1995-05-16 General Electric Company Adaptive filter for reducing streaking artifacts in x-ray tomographic images
RU2127075C1 (ru) * 1996-12-11 1999-03-10 Корженевский Александр Владимирович Способ получения томографического изображения тела и электроимпедансный томограф
JP3707347B2 (ja) * 2000-04-07 2005-10-19 株式会社島津製作所 X線ct装置の画像処理方法及びx線ct装置並びにx線ct撮影用記録媒体
US7046831B2 (en) 2001-03-09 2006-05-16 Tomotherapy Incorporated System and method for fusion-aligned reprojection of incomplete data
US6765981B2 (en) * 2002-07-31 2004-07-20 Agilent Technologies, Inc. Computed tomography
US7630528B2 (en) * 2004-03-02 2009-12-08 Koninklijke Philips Electronics N.V. Motion compensation
EP1741062B1 (en) * 2004-04-21 2013-02-13 Philips Intellectual Property & Standards GmbH Cone beam ct apparatus using truncated projections and a previously acquired 3d ct image
US7042974B2 (en) 2004-05-20 2006-05-09 Eastman Kodak Company Method for reconstructing tomographic images
US7468513B2 (en) * 2004-06-18 2008-12-23 The Children's Hospital Of Philadelphia Fast dynamic imaging protocol using a multi-head single photon emission computed tomography system
JP2008515513A (ja) * 2004-10-06 2008-05-15 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ コンピュータ断層撮影方法
EP1804667A4 (en) * 2004-10-12 2009-09-02 Univ Laval METHOD AND APPARATUS FOR REDUCING METAL ARTIFACTS IN COMPUTERIZED TOMOGRAPHY
JP2006212308A (ja) * 2005-02-07 2006-08-17 Ge Medical Systems Global Technology Co Llc 放射線断層撮影装置、放射線画像シミュレーション方法および画像シミュレーション装置
US7570732B2 (en) * 2005-11-09 2009-08-04 Dexela Limited Methods and apparatus for obtaining low-dose imaging
EP2011085A1 (en) 2006-04-25 2009-01-07 Wisconsin Alumni Research Foundation System and method for estimating data missing from ct imaging projections
US8131040B2 (en) * 2006-08-22 2012-03-06 Koninklijke Philips Electronics N.V. Artifact correction for motion artifacted images associated with the pulmonary cycle
US7680240B2 (en) * 2007-03-30 2010-03-16 General Electric Company Iterative reconstruction of tomographic image data method and system
US8472683B2 (en) * 2008-05-09 2013-06-25 General Electric Company Motion correction in tomographic images

Also Published As

Publication number Publication date
US20120148136A1 (en) 2012-06-14
US8666137B2 (en) 2014-03-04
WO2011027298A1 (en) 2011-03-10
CN102483853A (zh) 2012-05-30
EP2476099B1 (en) 2014-06-18
EP2476099A1 (en) 2012-07-18
RU2541860C2 (ru) 2015-02-20
CN102483853B (zh) 2015-08-26
JP2013503683A (ja) 2013-02-04
RU2012113544A (ru) 2013-10-20

Similar Documents

Publication Publication Date Title
JP5600742B2 (ja) 投射のデータを処理するための装置及び方法
US20180330233A1 (en) Machine learning based scatter correction
Wang et al. Scatter correction for cone‐beam computed tomography using moving blocker strips: a preliminary study
JP4532005B2 (ja) X線ct装置及びその画像表示方法
JP6066596B2 (ja) X線撮像における散乱補正の方法及びシステム
US9140803B2 (en) Acquisition protocol assessment apparatus
US9498177B2 (en) Contrast-dependent regularization strength in the iterative reconstruction of CT images
US20130202080A1 (en) System and Method for Denoising Medical Images Adaptive to Local Noise
US10282872B2 (en) Noise reduction in tomograms
CN110602990B (zh) 用于锥形束计算机断层扫描的患者移动校正方法
JP2008505676A (ja) 心臓の円錐ビームct再生における筋状アーチファクト低減
Walker et al. Evaluation of data-driven respiratory gating waveforms for clinical PET imaging
JP2011503570A (ja) 減衰マップを形成するための装置及び方法
US20150247810A1 (en) Method and apparatus for controlling dose in computed x-ray imaging
JP2016152916A (ja) X線コンピュータ断層撮像装置及び医用画像処理装置
RU2599867C2 (ru) Устройство для формирования изображений
JP5514195B2 (ja) 運動する物体の画像を生成するための装置
EP2490180B1 (en) Medical image processing apparatus and medical image imaging method
US10143434B2 (en) Imaging system for generating an image of an object
JP2010269165A (ja) X線ct装置
WO2008099314A2 (en) Apparatus and method for determining high density shadows in projection data
JP4644292B2 (ja) X線ct装置とその画像表示方法
Sun Rigid motion correction for head CT imaging
Buzug Image quality and artifacts
Rajwade Partial-Data Interpolation During Arcing of an X-Ray Tube in a Computed Tomography Scanner

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20130830

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20140228

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140325

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140624

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20140722

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140818

R150 Certificate of patent or registration of utility model

Ref document number: 5600742

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees