JP2002177254A - 造影剤に覆われた被検体のダイナミック３次元画像を再構成するための方法及びデバイス - Google Patents

造影剤に覆われた被検体のダイナミック３次元画像を再構成するための方法及びデバイス

Info

Publication number
JP2002177254A
JP2002177254A JP2001232982A JP2001232982A JP2002177254A JP 2002177254 A JP2002177254 A JP 2002177254A JP 2001232982 A JP2001232982 A JP 2001232982A JP 2001232982 A JP2001232982 A JP 2001232982A JP 2002177254 A JP2002177254 A JP 2002177254A
Authority
JP
Japan
Prior art keywords
subject
dimensional image
dimensional
image
contrast agent
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.)
Granted
Application number
JP2001232982A
Other languages
English (en)
Other versions
JP4745551B2 (ja
Inventor
Regis Vaillant
レジス・ヴァイアン
Dorothee Freymann
ドロシー・フレイマン
Laurent Launay
ローラン・ローナイ
Erwan Kerrien
エルワン・ケリエン
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.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
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 GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Publication of JP2002177254A publication Critical patent/JP2002177254A/ja
Application granted granted Critical
Publication of JP4745551B2 publication Critical patent/JP4745551B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • 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/58—Testing, adjusting or calibrating thereof
    • A61B6/582—Calibration
    • A61B6/583—Calibration using calibration phantoms
    • A61B6/584—Calibration using calibration phantoms determining position of components of the apparatus or device using images of the phantom
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/404—Angiography
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/412—Dynamic

Landscapes

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

Abstract

(57)【要約】 【課題】 造影剤に覆われた被検体のダイナミック３次
元画像を再構成するための方法及びデバイスを提供す
る。 【解決手段】 本方法は、被検体（１）の周りの様々な
カメラ位置に対してそれぞれ取得される第１組の２次元
ディジタル投影画像から被検体の静止３次元画像ＩＦを
再構成するフェーズと、同じ第１のカメラ位置に対して
ｎ回の連続する造影剤伝播時刻において少なくとも第２
組のｎ個の静止２次元投影画像ＩＢｉを収集するフェー
ズ５０と、第２組の各静止２次元画像ＩＢｉ及び再構成
した静止３次元画像ＩＦから被検体のダイナミック３次
元画像ＩＦＤを再構成するフェーズ５１と、を含む。

Description

【発明の詳細な説明】
【０００１】
【発明の属する技術分野】本発明は、造影剤に覆われた
被検体の３次元画像の再構成に関する。
【０００２】検査中の患者の内部構造の再構成、特に血
管造影画像の再構成（すなわち、造影剤の注入により不
透明にした脈管構造の画像の取得）を実施しているよう
な医療分野において、本発明の応用は特に重要である。
【０００３】
【発明の背景】しかしながら、本発明は、その他の分
野、とりわけ、医学的検査と同じタイプの検査が実行さ
れる非破壊的工業制御の分野にも応用できる。
【０００４】医療分野では、被検体（例えば、患者の頭
部）の２次元投影画像は、一般に、Ｘ線カメラを被検体
の周りに回転させることにより取得している。Ｘ線イメ
ージングでの再構成アルゴリズムには本質的に２つのタ
イプがある。第１のタイプでは、逆投影及びフィルタ処
理による計算や、さらには数次元でのフーリエ変換再構
成を規定している。第２のタイプは、代数法とも呼ば
れ、反復再構成法に関するものである。こうした代数的
アルゴリズムの原理は、例えば、フランス特許出願第８
９０３６０６号、同第８９１６９０６号または同第９８
０７３７１号に開示されており、これらの出願には、多
重分解能ボリュームに対する画像の代数的再構成での反
復アルゴリズムの応用について記載されている。
【０００５】要約して言えば、装置較正を使用して基本
的ボリューム要素（すなわち、ボクセル）に分解された
観察ボリュームの収集画像の投影面における投影パラメ
ータ（これらの較正パラメータにより投影マトリックス
が形成される）を決定した後、代数的画像再構成アルゴ
リズムを使用してこれらの２次元投影画像から３次元ボ
リュームを再構成させる。このアルゴリズムの基本原理
は、ボリュームの各ボクセルを所定の初期値（例えば、
ゼロ値）に初期化すること、並びに、仮想画像を取得す
るために収集した各画像の面にボクセルを投影する演
算、投影ボリューム（仮想画像）と対応する収集画像と
の差を決定する演算、さらにボリュームのこの差の逆投
影演算という３つの演算を多数回反復することである。
多数回反復した後、Ｘ線撮影した血管内に注入された造
影剤の密度を表す推定値が各ボクセルがごとに取得さ
れ、これによりこれらＸ線撮影した血管の地図を３次元
で視覚化することが可能となる。
【０００６】これらの３次元画像は、診断のためである
か、治療処置計画のためであるか、あるいは検査対象の
形状及び大きさの評価のためであるかを問わず、神経科
医や外科医に対して貴重な支援となる。
【０００７】一方、こうして再構成したボリューム画像
は１つの大きな欠点を有する。現実に、こうした画像で
は、３次元画像を再構成させる元になる収集２次元画像
がほとんど静止状態の造影剤に対して収集されたもので
あるため、動脈に注入される造影剤の伝播を視覚化する
ことができない。実際に、造影剤がほとんど静止の特性
であるため本明細書において「静止(static)」という語
を先頭に付けて記載するような３次元画像の再構成で
は、むしろ、同じ程度の造影剤伝播に対応する画像の組
を取得することが追求される。
【０００８】その結果、こうした再構成では、単に３次
元情報及び血流情報（造影剤の伝播）によっては患者の
脈管系を解析することができない。
【０００９】
【発明の概要】本発明の実施の一形態は、上述した静止
３次元画像とは異なり、本明細書において「ダイナミッ
ク」と記載するような、被検体内の造影剤伝播を３次元
で視覚化することが可能な３次元画像の再構成を対象と
する。
【００１０】本発明の実施の一形態は、造影剤に覆われ
た被検体のダイナミック３次元画像を再構成する方法で
あって、［ａ］被検体の周りの異なるイメージング位置
に対してそれぞれ取得される、被検体の２次元ディジタ
ル化投影画像の第１の組から被検体の静止３次元画像を
再構成（すなわち、再形成）するステップと、［ｂ］同
じイメージング位置に対してｎ回の連続する造影剤伝播
時刻においてそれぞれ取得される、ｎ個の静止２次元投
影画像からなる少なくとも第２の組を収集するステップ
と、［ｃ］第２の組の静止２次元画像の各々と再構成
（すなわち、再形成）した静止３次元画像とから、被検
体のダイナミック３次元画像を再構成するステップと、
を含む方法を対象とする。
【００１１】本発明はまた、造影剤に覆われた被検体の
ダイナミック３次元画像を再構成させるためのデバイス
であって、［ａ］被検体の周りの異なるイメージング位
置に対してそれぞれ取得された被検体の２次元ディジタ
ル化投影画像の第１の組から、被検体の静止３次元画像
を再構成する第１の手段と、［ｂ］ｎ回の連続する造影
剤伝播時刻に対応して、同じ第１のカメラ位置からｎ個
の静止２次元投影画像からなる少なくとも第２の組を収
集する手段と、［ｃ］第２の組の各静止２次元画像から
及び再構成した静止３次元画像から、被検体のダイナミ
ック３次元画像を再構成するための第２の手段と、を備
えるデバイスを提案する。
【００１２】本発明はまた、以下で規定する方法の使用
に適した、造影剤に覆われた被検体のダイナミック３次
元画像を再構成させるためのデバイスに関するものであ
る。
【００１３】本方法の使用を可能にするような別の手段
は、マイクロプロセッサにより実行させるソフトウェア
手段である。このソフトウェア手段は、読出し専用メモ
リ及び／またはディスケットなどのサポートで提供する
ことができる。したがって、本発明はさらに、プロセッ
サで実行させる際に以下に規定する方法を利用している
プログラムコード手段を備えたコンピュータ・プログラ
ムを提案する。
【００１４】本発明はさらに、プロセッサにより読み取
りが可能であって、プログラムをプロセッサで実行させ
る際に本方法を適用できるプログラムコード手段を含ん
でいるような、読出し専用メモリやディスケットなどの
サポートを提案する。
【００１５】本発明のその他の利点及び特徴は、非限定
的な使用法及び実施形態に関する詳細な記述、並びに添
付の図面を検討することにより明らかとなろう。
【００１６】
【発明の実施の形態】本方法の実施の一形態によれば、
本方法は、被検体を取り囲みかつボクセルに分解される
仮想ボリュームを生成させるような装置較正を含む。次
いで再構成した静止３次元画像は、その各値が仮想ボリ
ュームのボクセルにそれぞれ関連しており、被検体に注
入された造影剤の密度に関する対応するボクセルを表し
ている静止推定値と呼ぶ推定値から構成されている。ダ
イナミック３次元画像は、実際に、ｎ回の造影剤伝播時
刻に対応するｎ個の基本３次元画像から構成されてい
る。現在の伝播時刻に対応する現在の基本３次元画像
は、（１）被検体に注入された造影剤の密度の静止推定
値から、（２）各ボクセルごとに、ボクセルをその上に
投影させている現在の伝播時刻に対応する静止２次元画
像のピクセル強度から、並びに（３）ピクセルの各々に
対して当該ピクセルに関連する観察線に沿って位置する
仮想ボリュームのボクセルの密度推定値の積算値から、
再構成される。
【００１７】幾つかの場合では、数本の血管が同じピク
セル上に投影されることがあり得る。こうしたケース
は、一般に、ボクセル全体の１０％に相当する。これに
より、再構成したダイナミック３次元画像の推定は曖昧
となる可能性がある。こうした曖昧さは、第２の組のｎ
個の２次元画像と異なる観測点から収集した、別系列の
静止２次元画像を用いることにより解決できる。
【００１８】換言すると、別の実施形態では、第１の位
置（第２の組の２次元投影画像の収集フェーズに対応す
る）とは別の同じ第２のカメラ位置に対して、同じ連続
する造影剤伝播時刻でそれぞれ収集される第３の組の２
次元投影画像の収集フェーズが設けられる。次いで、再
形成した静止３次元画像から、第２の組の各画像から、
並びに第３の組の各画像から、被検体のダイナミック３
次元画像が再構成される。
【００１９】本発明を限定するものではないが、患者
（詳細には、患者の頭部）のダイナミック３次元血管造
影画像の再構成に対する本方法の応用について記載す
る。
【００２０】特に、図１及び２を参照すると、本発明を
適用するために使用可能なイメージング・システムによ
り、とりわけ、この場合では患者の頭部１の周りでのＸ
線源２の回転により得られる２次元収集画像の組ＩＡ１
〜ＩＡｍを取得することが可能となる。
【００２１】血管造影では標準的なことであるが、収集
した各画像ＩＡｉは、例えば、脈管構造への造影剤注入
の前後に同じ入射角で撮影した２枚のＸ線像に対する標
準技法の対数サブトラクションにより取得されるような
サブトラクション画像であり、これらの３次元画像は、
記載するように、先ず静止像で次いでダイナミック像で
再構成させることが望ましい。
【００２２】より厳密には、造影剤注入の前と後で、患
者の頭部の周りでの回転による２系統を収集する。これ
を実行するためには、放射線撮影装置により、例えば、
概ね１０秒で１９３．５度の回転を２度させて、４．５
度ごとに４４枚の２次元画像を収集することができる。
画像ＩＡｉの収集は、静止と呼ぶ、すなわち、造影剤伝
播に関係しない脈管構造の画像の３次元再構成を目的と
している。これらの条件では、造影剤が検査対象血管内
でほとんど静止している時点でその収集を行うことが好
ましい。
【００２３】各々の収集した画像ＩＡｉは、「投影面Ｐ
Ｐｉ」と呼ぶ平面内でＸ線管の反対側に配置した、例え
ば放射線医学で使用されるルミネセンス増幅器タイプの
ような２次元放射線検出器から取得される。患者の頭部
の周りでの回転により検出器の異なる角度位置で様々な
投影面が得られる。
【００２４】投影面ＰＰｉに対する法線ＸＯｉは、収集
画像ＩＡｉの光軸を規定している。検出器は、とりわ
け、本発明で使用する画像再構成アルゴリズムを付属の
プログラム・メモリ内にソフトウェアとして組み込んで
いるマイクロプロセッサに接続させたサンプリング手段
と、一般に、別の画像再構成手段など、本方法の使用を
可能にしている機能手段のすべてと、を含んだ処理手段
３と接続している。
【００２５】２次元検出器を伴うＸ線源で構成したＸ線
イメージング・システムの場合では、収集画像の作成に
到るための幾何学演算は、解析する対象を、３次元空間
で展開し、その検出面に対応する投影面の空間にあたる
２次元空間上に円錐投影することである。様々な円錐投
影を完全に記述する幾何学的パラメータが既知である必
要がある。これらのパラメータに直接アクセスするこ
と、すなわち、例えば、収集システム上でＸ線源と検出
器の間の距離を直接計測することなどによりアクセスす
ることは、多くの場合不可能であり、また不正確であ
る。
【００２６】画像作成に到るための幾何学的パラメータ
に関する間接的な把握をもたらす演算を、イメージング
・システムの「較正(calibration) 」と呼ぶ。
【００２７】この原理は、標準的かつ周知であり、３次
元空間における周知の幾何学的ファントムを使用してそ
の２次元投影を収集することに基づいている。より厳密
には、この較正には以下のステージ、すなわち（１）既
知の対象物として、当該対象物に特有のマークを基準と
して計測した座標によってその空間的位置を知ることが
できる多数の特性点を提示する較正ファントムを使用す
るステージ、（２）較正の対象である観測点（または、
入射角）の幾何学条件下で当該ファントムの画像を収集
するステージ、及び（３）この画像において特性点の投
影を識別するステージ、という各ステージが必要であ
る。このためには、収集した投影画像において、対象物
の各特性点を結合させてトレース（軌跡）をつくる、す
なわち、（４）投影を記述している方程式体系を逆転さ
せ、さらに（５）最終的に、所与の観測点に対する投影
パラメータの組を取得する。
【００２８】よく使用される幾何学的較正ファントムの
形状は、その８つの隅にＸ線に対して不透明な金属球を
配置した立方体の形状である。較正は当業者には周知の
操作であるため、さらに詳細には説明しない。
【００２９】Ｘ線によるイメージング・システムの自動
的幾何学較正に関する周知の一方法としては、フランス
特許出願第９３００８０４号に記載された方法などを挙
げることができる。簡単に言えば、こうした自動的較正
のためにはファントムを使用するが、このファントムで
は、イメージング・システムの回転軸（特に、ファント
ムの軸）に沿って計測される球の高さが、例えばその順
列内の球の直列数に伴って単調に増加（または、減少）
するように、球を順列内の方々に分散させている。
【００３０】イメージング・システムの較正により、被
検体の周りのカメラの推定平均回転軸Ａｘの決定、並び
に線源２の位置や収集した様々な画像の光軸の幾何学的
特性の決定が可能となる。較正によりさらに、被検体１
を取り囲みかつ基本的ボリューム要素Ｖｉ（すなわち
「ボクセル」）に分解できる仮想ボリュームＶＶ（様々
な投影円錐の交差箇所）の規定が可能となる。このボリ
ュームＶＶ（したがって、各ボクセルＶｉ）は、以下に
おいて較正基準系と呼ぶ基準系において空間的にマーク
されている。この基準系の軸のうちの１つ（この場合で
は軸Ｚ）は推定回転軸Ａｘに合併している。ここで、収
集画像ＩＡｉを投影させる投影面ＰＰｉは一般に、軸Ｚ
に平行であることに留意されたい。
【００３１】較正によりさらに、収集した各画像ＩＡｉ
に対して投影マトリックスＰｉを規定することができ、
これにより、各ボクセルＶｉごとに、対応する収集画像
ＩＡｉにおける投影（ピクセル）の座標の決定が可能と
なる。
【００３２】ここで、得られた第１の組の投影静止２次
元画像ＩＡｉから、いわゆる「静止」３次元画像を再構
成するフェーズを、特に図４を参照しながら説明するこ
とにする。
【００３３】再形成した静止３次元画像ＩＦを、典型的
には３回の反復の後に出力するために、代数的画像再構
成アルゴリズムが画像ＩＡｉ（較正３０及び収集３１の
後に取得される）に直接適用される。
【００３４】ここで使用する画像再構成アルゴリズムの
基本原理は、ボリュームのボクセルを所定の初期値（例
えば、ゼロ値）に初期化すること、次いで多数回（典型
的には、３回）反復すること、からなることを想起され
たい。各反復では、収集した各画像ごとに、仮想画像を
取得するために収集した各画像の面にボクセルを投影す
ること、投影ボリューム（仮想画像）と対応する収集画
像との間の差を決定すること、並びにボリューム内に当
該の差を逆投影すること、が必要である。
【００３５】一般に、３回の反復の後、座標（ｉ，ｊ，
ｋ）の各ボクセルに対して、Ｘ線撮影した血管内に注入
された造影剤の密度を表す推定値ｃi,j,k （ここでは静
止推定値を用いて表記される）が得られ、これにより、
これらのＸ線撮影した血管（画像ＩＦ）の地図を３次元
で視覚化することが可能となる。
【００３６】もちろん、本発明の範囲内では実際上、こ
うして得られる静止３次元画像ＩＦがここでは再構成し
たダイナミック３次元画像を得るための中間的ステージ
（これについては、以下でさらに詳細に記載する）のみ
から構成されるため、静止３次元画像ＩＦを必ずしもイ
メージング・システムの表示スクリーン上で視覚化させ
る必要はない。
【００３７】視覚化を希望するならば、当業者は、画像
再構成の代数的アルゴリズムの使用に関するより詳細に
ついて、上述のフランス特許出願第９８０７３７１号を
参照することができる。
【００３８】第１の組の画像ＩＡｉの収集の外に、第２
の組の静止双方向投影画像ＩＢ１〜ＩＢｎ（図３）の収
集を実施する。この一連の画像ＩＢｉは、同じ観測点に
対応した同じカメラ位置ＰＤＶ１に対して収集される。
【００３９】したがって、ｎ回の連続する造影剤伝播時
刻に対応したｎ個の画像ＩＢｔが収集される。
【００４０】図５に示すように、収集ステージ５０の間
に収集したこれらの画像ＩＢｉから、並びに再形成した
静止３次元画像ＩＦから、ステージ５１において脈管構
造の３次元ダイナミック画像を再形成している。
【００４１】ここで、較正の後に収集３１（図４）と収
集５０（図５）を実行する順序は、本発明による方法を
使用するためには重要でないことに留意されたい。
【００４２】再形成したダイナミック３次元画像ＩＦＤ
は、実際には、ｎ個の基本ダイナミック３次元画像から
構成されており、その各画像に対して脈管構造内での造
影剤の伝播状態を所与の伝播時刻ｔにおいて視覚化する
ことができる。ｎ個の基本３次元画像を連続的に視覚化
することにより、脈管構造内での造影剤伝播に対するダ
イナミックな視覚化が可能となる。
【００４３】再形成した静止３次元画像に関して上述し
たのと同様に、再形成したダイナミック３次元画像の視
覚化は、各座標（ｉ，ｊ，ｋ）値に対する伝播時刻ｔに
おける造影剤密度の推定値により得られる。伝播時刻ｔ
において妥当な推定値をｃi, j,k （ｔ）で表す。
【００４４】この推定値は次の式 ｃi,j,k（ｔ）＝［ｃi,j,k ・ｐ（ｔ）］／Ｅ ［１］ で規定される。
【００４５】上式において、ｃi,j,k は、再形成した静
止３次元画像ＩＦの同じボクセルに対する造影剤の静止
推定値を意味し、ｐ（ｔ）は、伝播時刻ｔにおいて収集
し、かつその上に座標ｉ，ｊ，ｋのボクセルを投影させ
る画像ＩＢｔのピクセル強度を意味し、Ｅは、当該ピク
セルｐに対して、当該ピクセルｐに関連する観察線に沿
って位置している仮想ボリュームＶＶの各ボクセルの静
止密度推定値の積算値を意味している。
【００４６】Ｅは次式 Ｅ＝Σｃi,j,k λi,j,k ｉ，ｊ，ｋ ［２］ により規定される。式［２］において、λi,j,k の値
は、座標ｉ，ｊ，ｋのボクセルが観察線上に位置してい
ない場合には０に等しく、また観察線上に位置している
場合には１に等しい。
【００４７】ここで、所定の観察線に沿って撮影された
一連の静止画像に関連する投影マトリックスが既知であ
るため、ピクセルｐの決定が可能となることに留意され
たい。同様に、式［２］の値λもこの投影マトリックス
の値から既知となる。
【００４８】ここで、時刻ｔにおいて座標ｉ，ｊ，ｋの
ボクセルにまだ造影剤が到達していない場合、ｃi,j,k
（ｔ）の値は小さくなるはずであることに留意された
い。値ｐ（ｔ）が小さい場合、すなわち、投影したピク
セルの密度が小さい場合は、したがって式［１］より時
刻ｔにおけるボクセルの密度も小さくなることが十分に
分かる。
【００４９】一方、値ｐ（ｔ）が大きい場合には、当該
の観察線が、当該の座標ｉ，ｊ，ｋのボクセルによる密
度に食い込む(parasite)ことになるような別の高輝度(l
uminous)ボクセルと交差していることを意味する。した
がって、式［１］によりそのボクセルをそのままの場合
と比べてさらに利用しやすく変更している。
【００５０】座標ｉ，ｊ，ｋのボクセルに造影剤が到達
すると、上述した場合のように、当該のボクセルが造影
剤のまだ到達していないボクセルとが重なった場合に、
式［１］ではアーチファクトが生成されることがある。
【００５１】しかし、上述したアーチファクト生成が及
ぶ範囲がボリュームの全ボクセルの概ね１０％だけであ
ると考えると、上述の式［１］により、再形成したダイ
ナミック画像に対する極めて良好な推定が可能となる。
【００５２】しかしながら、こうした曖昧さを解決し、
結果として画像のダイナミック再構成を向上させるため
には、図６に示すように、第１の系列の画像ＩＢｉに対
する観測点と異なる観測点からここに投影される別系列
のｎ個の静止２次元画像の収集を行うことができる。画
像ＩＢｉ及びＩＣｉはそれぞれ同時に収集される。
【００５３】この場合では、図５に点線で示すように、
ステージ５２において収集した画像ＩＣｉからも、ステ
ージ５１におけるダイナミック３次元画像の再形成を実
行する。
【００５４】次いで、画像ＩＢｉを用いて、時刻ｔにお
けるボクセルｉ，ｊ，ｋに対するｃ i,j,k （ｔ）の第１
の値が計算され、さらに今度は、式［１］を画像ＩＣｉ
と共に用いることによってｃi,j,k （ｔ）の第２の値が
計算される。
【００５５】そしてここで、座標ｉ，ｊ，ｋのボクセル
に関連するｃi,j,k （ｔ）の最終値は、上述した２つの
値のうちの小さい方に等しくなる。
【００５６】当業者であれば、特許請求の範囲に示した
本発明の範囲及び領域を逸脱することなく、構造及び／
またはステップ及び／または機能に関して様々な修正を
行うことができる。
【図面の簡単な説明】
【図１】被検体の周りの２次元投影画像の組の模式図で
ある。
【図２】２次元投影画像の収集をより詳細に表した略図
である。
【図３】同じ観測点から収集される静止投影画像の組の
模式図である。
【図４】本方法の実施の一形態の流れ図である。
【図５】本方法の実施の一形態の流れ図である。
【図６】本発明の別の実施形態の模式図である。
【符号の説明】
１ 被検体 ２ Ｘ線源 ３ 処理手段 Ａｘ カメラの推定回転軸 ＩＡｉ 第１組の静止２次元画像 ＩＢｉ 第２組の静止２次元画像 ＩＣｉ 第３組の静止２次元画像 ＰＤＶ１ 第１のイメージング位置 Ｖｉ ボクセル ＶＶ 仮想ボリューム
───────────────────────────────────────────────────── フロントページの続き (72)発明者 ドロシー・フレイマン フランス、67000・ストラスブール、リ ュ・ド・ロンドル、９番 (72)発明者 ローラン・ローナイ フランス、78470・エスティー・ルミー・ レ・シェヴルユゼ、アンパス・ド・サル ジ、11番 (72)発明者 エルワン・ケリエン フランス、78000・ヴェルサイユ、リュ・ ドュ・マレシャル・ガリエニ、14番 Ｆターム(参考） 4C093 AA08 AA24 CA23 DA02 EB17 FE01 FE14 FF34 FF42 GA01 5B057 AA07 AA09 BA03 CA08 CA13 CA16 CB08 CB13 CB16 CE08 5B080 AA17 GA00

Claims (7)

    【特許請求の範囲】
  1. 【請求項１】 造影剤に覆われた被検体のダイナミック
    ３次元画像を再構成する方法であって、 （ａ）被検体の周りの異なるイメージング位置に対して
    それぞれ取得される、被検体の２次元ディジタル化投影
    画像の第１の組（ＩＡｉ）から被検体の静止３次元画像
    （ＩＦ）を再構成する（３３）ステップと、 （ｂ）同じ第１のイメージング位置（ＰＤＶ１）に対し
    てｎ回の連続する造影剤伝播時刻においてそれぞれ取得
    される、ｎ個の静止２次元投影画像からなる少なくとも
    第２の組（ＩＢｉ）の収集フェーズ（５０）のステップ
    と、 （ｃ）第２の組の静止２次元画像（ＩＢｉ）の各々と再
    構成した静止３次元画像（ＩＦ）とから、被検体のダイ
    ナミック３次元画像（ＩＦＤ）を再構成する（５１）ス
    テップと、を含む方法。
  2. 【請求項２】被検体を取り囲みかつボクセル（Ｖｉ）に
    分解される仮想ボリューム（ＶＶ）を生成させること、 次いで再構成した静止３次元画像（ＩＦ）は、その各値
    が仮想ボリュームのボクセルにそれぞれ関連しており、
    被検体に注入された造影剤の密度に関する対応するボク
    セルを表している静止推定値（ｃi,j,k）から構成され
    ていること、 前記ダイナミック３次元画像（ＩＦＤ）は実際にｎ回の
    造影剤伝播時刻に対応するｎ個の基本３次元画像から構
    成されていること、並びに、 （１）被検体に注入された造影剤の密度の静止推定値
    （ｃi,j,k）から、 （２）各ボクセル（ｉ，ｊ，ｋ）ごとに、該ボクセル
    （ｉ，ｊ，ｋ）をその上に投影させている現在の伝播時
    刻に対応する静止２次元画像のピクセル強度（ｐ
    （ｔ））から、 （３）前記ピクセルの各々に対して、当該ピクセルに関
    連する観察線に沿って位置する仮想ボリュームのボクセ
    ルの密度推定値の積算値（Ｅ）から、現在の伝播時刻に
    対応する現在の基本３次元画像を再構成すること、を特
    徴とする装置較正（３０）のステップを含む、請求項１
    に記載の方法。
  3. 【請求項３】 第１の位置とは別の同じ第２のカメラ位
    置に対して同じ連続する造影剤伝播時刻においてそれぞ
    れ取得された静止２次元投影画像からなる第３の組（Ｉ
    Ｃｉ）を収集する（５２）ステップを含むと共に、被検
    体のダイナミック３次元画像を静止３次元画像から、第
    ２の組の各画像から、並びに第３の組の各画像から再構
    成することを特徴とする、請求項１または２に記載の方
    法。
  4. 【請求項４】 造影剤に覆われた被検体のダイナミック
    ３次元画像を再構成させるためのデバイスであって、 （１）被検体の周りの異なるイメージング位置に対して
    それぞれ取得される、被検体の２次元ディジタル化投影
    画像の第１の組（ＩＡｉ）から、被検体の静止３次元画
    像（ＩＦ）を再構成するための第１の手段（３）と、 （２）ｎ回の連続する造影剤伝播時刻に対応して、同じ
    第１のイメージング位置からｎ個の静止２次元投影画像
    からなる少なくとも第２の組（ＩＢｉ）を収集するため
    の手段と、 （３）第２の組の各静止２次元画像から、及び再構成し
    た静止３次元画像から被検体のダイナミック３次元画像
    を再構成するための第２の手段（３）と、を備えるデバ
    イス。
  5. 【請求項５】 前記請求項１から３までのいずれか一項
    に記載の方法を使用するために、造影剤に覆われた被検
    体のダイナミック３次元画像を再構成させるためのデバ
    イス。
  6. 【請求項６】 プロセッサで実行させる際に前記請求項
    １から３までのいずれか一項の規定に従った方法を利用
    しているプログラムコード手段を備えるコンピュータ・
    プログラム。
  7. 【請求項７】 プロセッサにより読み取り可能であり、
    かつプロセッサで実行させる際に前記請求項１から３ま
    でのいずれか一項の規定に従った方法を利用できるプロ
    グラムコード手段を含んでいるサポート。
JP2001232982A 2000-08-02 2001-08-01 造影剤に覆われた被検体のダイナミック３次元画像を再構成するための方法及びデバイス Expired - Lifetime JP4745551B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR0010211A FR2812741B1 (fr) 2000-08-02 2000-08-02 Procede et dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste
FR0010211 2000-08-02

Publications (2)

Publication Number Publication Date
JP2002177254A true JP2002177254A (ja) 2002-06-25
JP4745551B2 JP4745551B2 (ja) 2011-08-10

Family

ID=8853236

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001232982A Expired - Lifetime JP4745551B2 (ja) 2000-08-02 2001-08-01 造影剤に覆われた被検体のダイナミック３次元画像を再構成するための方法及びデバイス

Country Status (4)

Country Link
US (1) US7113631B2 (ja)
JP (1) JP4745551B2 (ja)
DE (1) DE10137521A1 (ja)
FR (1) FR2812741B1 (ja)

Families Citing this family (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7583992B2 (en) * 2000-03-30 2009-09-01 Wisconsin Alumni Research Foundation Magnetic resonance angiography with automated vessel segmentation
DE10215808B4 (de) * 2002-04-10 2005-02-24 Siemens Ag Verfahren zur Registrierung für navigationsgeführte Eingriffe
EP1709617A2 (en) * 2003-12-30 2006-10-11 Trustees Of The Stevens Institute Of Technology Three-dimensional imaging system using optical pulses, non-linear optical mixers and holographic calibration
US7693318B1 (en) 2004-01-12 2010-04-06 Pme Ip Australia Pty Ltd Method and apparatus for reconstruction of 3D image volumes from projection images
CN1664861A (zh) * 2004-09-22 2005-09-07 冯彦 一种动态标识的制作及显示方法
US8189002B1 (en) 2004-10-29 2012-05-29 PME IP Australia Pty, Ltd. Method and apparatus for visualizing three-dimensional and higher-dimensional image data sets
US7778392B1 (en) 2004-11-02 2010-08-17 Pme Ip Australia Pty Ltd Method of reconstructing computed tomography (CT) volumes suitable for execution on commodity central processing units (CPUs) and graphics processors, and apparatus operating in accord with those methods (rotational X-ray on GPUs)
US7609884B1 (en) 2004-12-23 2009-10-27 Pme Ip Australia Pty Ltd Mutual information based registration of 3D-image volumes on GPU using novel accelerated methods of histogram computation
CN101115989B (zh) * 2005-02-04 2012-09-05 皇家飞利浦电子股份有限公司 用于确定血管几何形状和流特征的系统
US7623732B1 (en) 2005-04-26 2009-11-24 Mercury Computer Systems, Inc. Method and apparatus for digital image filtering with discrete filter kernels using graphics hardware
US8411914B1 (en) * 2006-11-28 2013-04-02 The Charles Stark Draper Laboratory, Inc. Systems and methods for spatio-temporal analysis
DE602008003941D1 (de) * 2007-04-13 2011-01-27 Brainlab Ag Verfahren und Vorrichtung für die Planung lokaler Verabreichung eines Arzneimittels
US7853061B2 (en) * 2007-04-26 2010-12-14 General Electric Company System and method to improve visibility of an object in an imaged subject
DE102007024451B4 (de) * 2007-05-25 2010-07-08 Siemens Ag Ermittlungsverfahren für örtlich dreidimensional aufgelöste Anwesenheitsverteilungen einer Substanz in einem Gefäßsystem und hiermit korrespondierende Einrichtungen
US8019151B2 (en) * 2007-06-11 2011-09-13 Visualization Sciences Group, Inc. Methods and apparatus for image compression and decompression using graphics processing unit (GPU)
US8392529B2 (en) 2007-08-27 2013-03-05 Pme Ip Australia Pty Ltd Fast file server methods and systems
US9904969B1 (en) 2007-11-23 2018-02-27 PME IP Pty Ltd Multi-user multi-GPU render server apparatus and methods
WO2009067675A1 (en) 2007-11-23 2009-05-28 Mercury Computer Systems, Inc. Client-server visualization system with hybrid data processing
US8319781B2 (en) 2007-11-23 2012-11-27 Pme Ip Australia Pty Ltd Multi-user multi-GPU render server apparatus and methods
US8548215B2 (en) 2007-11-23 2013-10-01 Pme Ip Australia Pty Ltd Automatic image segmentation of a volume by comparing and correlating slice histograms with an anatomic atlas of average histograms
US10311541B2 (en) 2007-11-23 2019-06-04 PME IP Pty Ltd Multi-user multi-GPU render server apparatus and methods
US8693745B2 (en) * 2009-05-04 2014-04-08 Duke University Methods and computer program products for quantitative three-dimensional image correction and clinical parameter computation in optical coherence tomography
US8538106B2 (en) * 2009-10-20 2013-09-17 Siemens Aktiengesellschaft Three-dimensional esophageal reconstruction
US8787521B2 (en) * 2009-12-23 2014-07-22 General Electric Company System and method of iterative image reconstruction for computed tomography
ES2392229B1 (es) 2010-08-27 2013-10-16 Telefónica, S.A. Método de generación de un modelo de un objeto plano a partir de vistas del objeto.
EP2737852B1 (de) * 2012-11-30 2015-08-19 GE Sensing & Inspection Technologies GmbH Verfahren zum Erfassen geometrischer Abbildungseigenschaften eines Flachbilddetektors, entsprechend eingerichtete Röntgenprüfanlage und Kalibrierkörper
US10540803B2 (en) 2013-03-15 2020-01-21 PME IP Pty Ltd Method and system for rule-based display of sets of images
US8976190B1 (en) 2013-03-15 2015-03-10 Pme Ip Australia Pty Ltd Method and system for rule based display of sets of images
US10070839B2 (en) 2013-03-15 2018-09-11 PME IP Pty Ltd Apparatus and system for rule based visualization of digital breast tomosynthesis and other volumetric images
US11244495B2 (en) 2013-03-15 2022-02-08 PME IP Pty Ltd Method and system for rule based display of sets of images using image content derived parameters
US9509802B1 (en) 2013-03-15 2016-11-29 PME IP Pty Ltd Method and system FPOR transferring data to improve responsiveness when sending large data sets
US11183292B2 (en) 2013-03-15 2021-11-23 PME IP Pty Ltd Method and system for rule-based anonymized display and data export
US11599672B2 (en) 2015-07-31 2023-03-07 PME IP Pty Ltd Method and apparatus for anonymized display and data export
US9984478B2 (en) 2015-07-28 2018-05-29 PME IP Pty Ltd Apparatus and method for visualizing digital breast tomosynthesis and other volumetric images
US10909679B2 (en) 2017-09-24 2021-02-02 PME IP Pty Ltd Method and system for rule based display of sets of images using image content derived parameters

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0678916A (ja) * 1992-04-02 1994-03-22 Toshiba Corp Ｃｔ装置及びｃｔ装置を用いた血液量分布測定方法
JPH06205771A (ja) * 1992-09-09 1994-07-26 Picker Internatl Inc 高速掃引を用いた四次元螺旋状体積撮像方法及び装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2644590B1 (fr) 1989-03-20 1994-08-19 General Electric Cgr Sa Procede d'acquisition de donnees radiologiques et de reconstruction de structures correspondant a ce corps
FR2656129B1 (fr) 1989-12-20 1992-03-13 Gen Electric Cgr Procede de reconstruction multi-echelle de l'image de la structure d'un corps.
US5287273A (en) * 1990-03-15 1994-02-15 Mount Sinai School Of Medicine Functional organ images
FR2700909B1 (fr) 1993-01-27 1995-03-17 Gen Electric Cgr Dispositif et procédé automatique de calibration géométrique d'un système d'imagerie par rayons X.
US5827187A (en) * 1996-04-23 1998-10-27 Mayo Foundation For Medical Education And Research Dynamic MR digital subtraction angiography with complex subtraction
US6073042A (en) * 1997-09-25 2000-06-06 Siemens Medical Systems, Inc. Display of three-dimensional MRA images in which arteries can be distinguished from veins
US6230040B1 (en) * 1997-11-21 2001-05-08 Cornell Research Foundation, Inc. Method for performing magnetic resonance angiography with dynamic k-space sampling
FR2779853B1 (fr) * 1998-06-11 2000-08-11 Ge Medical Syst Sa Procede de reconstruction d'une image tridimensionnelle d'un objet, en particulier une image tridimensionnelle angiographique
US6195579B1 (en) * 1998-12-18 2001-02-27 Wisconsin Alumni Research Foundation Contrast detection and guided reconstruction in contrast-enhanced magnetic resonance angiography
CN1321616C (zh) * 2001-08-10 2007-06-20 皇家飞利浦电子股份有限公司 用于根据投影图像重建三维数据集的x－射线检查装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0678916A (ja) * 1992-04-02 1994-03-22 Toshiba Corp Ｃｔ装置及びｃｔ装置を用いた血液量分布測定方法
JPH06205771A (ja) * 1992-09-09 1994-07-26 Picker Internatl Inc 高速掃引を用いた四次元螺旋状体積撮像方法及び装置

Also Published As

Publication number Publication date
US7113631B2 (en) 2006-09-26
FR2812741B1 (fr) 2003-01-17
US20020123680A1 (en) 2002-09-05
JP4745551B2 (ja) 2011-08-10
DE10137521A1 (de) 2002-02-14
FR2812741A1 (fr) 2002-02-08

Similar Documents

Publication Publication Date Title
JP4745551B2 (ja) 造影剤に覆われた被検体のダイナミック３次元画像を再構成するための方法及びデバイス
US9754390B2 (en) Reconstruction of time-varying data
Saba et al. Maximizing quantitative accuracy of lung airway lumen and wall measures obtained from X-ray CT imaging
JP5737725B2 (ja) 局所化および高度に限定された画像再構成法
JP3761094B2 (ja) 対象物の３次元画像を再構成する方法
JP5268365B2 (ja) 脈管幾何形状及び流れ特性を決定するシステム
US9449381B2 (en) Methods, systems, and media for generating and analyzing medical images having elongated structures
Manhart et al. Dynamic iterative reconstruction for interventional 4-D C-arm CT perfusion imaging
JP5566370B2 (ja) 医用画像処理装置及び方法
JP2006255412A (ja) 腫瘍量を監視する方法及びシステム
JP2022027757A (ja) ４ｄコンピュータ断層撮影から呼吸性血液量変化を決定するシステム及び方法
US7355605B2 (en) Method and system for automatic orientation of local visualization techniques for vessel structures
EP4026102B1 (en) Confidence map for neural network based limited angle artifact reduction in cone beam ct
US11317875B2 (en) Reconstruction of flow data
US11730384B2 (en) Fluid analysis apparatus, method for operating fluid analysis apparatus, and fluid analysis program
EP3607527B1 (en) Quantitative evaluation of time-varying data
Liang et al. GUEST EDITORIAL:: MEDICAL IMAGING INFORMATICS—AN INFORMATION PROCESSING FROM IMAGE FORMATION TO VISUALIZATION.
Manhart et al. Guided noise reduction with streak removal for high speed flat detector CT perfusion
Deuerling-Zheng et al. Dynamic Iterative Reconstruction for Interventional 4-D C-Arm CT Perfusion Imaging

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20080730

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20101116

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20101116

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20101227

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110208

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110404

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: 20110419

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110512

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140520

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 4745551

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term