JP2009502397A - 3d−2d適応型形状モデル支援による運動補償再構成 - Google Patents

3d−2d適応型形状モデル支援による運動補償再構成 Download PDF

Info

Publication number
JP2009502397A
JP2009502397A JP2008524630A JP2008524630A JP2009502397A JP 2009502397 A JP2009502397 A JP 2009502397A JP 2008524630 A JP2008524630 A JP 2008524630A JP 2008524630 A JP2008524630 A JP 2008524630A JP 2009502397 A JP2009502397 A JP 2009502397A
Authority
JP
Japan
Prior art keywords
interest
shape model
image
motion
dimensional
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
JP2008524630A
Other languages
English (en)
Other versions
JP5244592B2 (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 JP2009502397A publication Critical patent/JP2009502397A/ja
Application granted granted Critical
Publication of JP5244592B2 publication Critical patent/JP5244592B2/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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/412Dynamic

Abstract

関心構造(60)に対応する3次元(3D)画像を生成あるいは再構成する方法は、関心構造(60)に対応する複数の画像投影を収集する段階、選定された3次元シード点(64)において形状モデル(66)を適用する段階、及び関心構造(60)を表すように形状モデル(66)を適応させ、それにより適応化された形状モデルを生じさせる段階を含んでいる。また、3次元(3D)画像の生成及び再構成のためのシステム(10)は、関心構造(60)に対応する画像投影データを提供するように構成された撮像システム(12)、及び撮像システム(12)と通信することが可能なコントローラ(50)を含んでいる。コントローラ(50)は、投影データを受信し、選定された3次元シード点(64)において形状モデル(66)を適用し;且つ関心構造(60)を表すように形状モデル(66)を適応させ、それにより適応化された形状モデルを生じさせるように構成されている。

Description

本発明は、3D及び4D画像再構成において運動を補償する方法に関し、具体的には、適応型形状モデルに基づく、X線透視法により生成された画像などの運動補償及び増強に関する。
例えば電気生理的診療のためのようなX線誘導心臓診療においては、診療の計画及び誘導のために、標的心室構造の3次元(3D)及び4次元(4D)再構成がしばしば利用されている。現在、これらのデータは必要に応じて、診療前に異なる画像診断方法を用いて収集されている。しかしながら、この手法を用いるとデータは最新のものでない場合がある。また、異なる画像診断方法を用いて収集されたデータは、誘導に実際に使用される画像情報に対して整合されなければならず、コスト、時間及び複雑さを増すことになる。
さらに、如何なる種類の撮像においても、その間の患者の運動は一貫性のないデータを生じさせ、ひいては、例えばブラーリング(ぼけ)やゴースト画像などのアーチファクトを生じさせる。故に、患者の運動は回避あるいは補償されなければならない。実際には、例えば患者の固定など運動を回避することは、一般に、困難であるか不可能である。故に、患者の運動を補償することの方が実際的である。運動補償手法の多くは、全てのデータが同一の運動状態に属するような一貫性を有する投影データを如何にして取得するかに焦点を当てており、投影データのこのような部分集合を再構成に用いている。このような部分集合を複数用いることにより、被検体の相異なる運動状態を再構成することができる。例えば、或る1つの手法は、被検体の運動及びX線減衰の時間進展を補償するために平行リビニング円錐ビーム逆投影法を用いている。運動の場(フィールド)がスライディングウィンドウ再構成のブロックマッチングによって見積もられ、考慮中のボクセルに関する一貫したデータが、全ての投影角に対して、同一方向からの時間的に隣接するデータからの線形回帰によって近似される。そのボクセルに関するフィルタリングされた投影データが運動ベクトル場に従って選択される。その他の方法は、投影の演算子を修正し運動補償された再構成を計算するために、予め計算された運動の場を用いて、画像再構成における運動の影響に対処している。
これまでの努力に拘わらず、依然として、3D/4Dデータセットを生成するための効率的且つコスト効率の良い方法が望まれる。診療にも使用される撮像システム上で3D/4Dデータセットを生成することが特に有利である。また、追加の実験時間及び画像の整合に関する追加の費用及び時間を回避するためには、この撮像をそれと同時に行われる診療のために実施することが有利である。
本発明は、3D/4Dデータセットを生成するための効率的且つコスト効率の良い方法を提供することを目的とする。
典型的な一実施形態において、関心構造に対応する画像投影データに基づいて3次元(3D)画像を生成あるいは再構成する方法が開示される。この方法は:画像投影からの3次元画像の生成又は再構成を容易にするのに十分な角度範囲で、関心構造に対応する複数の画像投影を収集する段階;3次元シード点(起点)を選定する段階;3次元シード点において形状モデルを適用する段階;及び関心構造を表すように形状モデルを適応させ、それにより適応化された形状モデルを生じさせる段階;を含んでいる。
典型的な実施形態によれば、他の必要に応じての一実施形態において、上述の方法は更に:画像投影に関連付けて関心構造の運動を指し示すデータを収集する段階;関心構造の運動を指し示すデータに基づいて、関心構造を表すように形状モデルを適応させ、それにより別の適応化された形状モデルを生じさせる段階;及びこの別の適応化された形状モデルに基づいて関心構造の3次元画像を生成する段階;を含んでいてもよい。
さらに、他の典型的な一実施形態において、3次元(3D)画像の生成及び再構成のためのシステムが開示される。このシステムは:画像投影からの3次元画像の生成又は再構成を容易にするのに十分な角度範囲で、関心構造に対応する画像投影データを提供するように構成された撮像システム;及び撮像システムと通信することが可能なコントローラを含んでいる。このコントローラは:画像投影データを受信し;3次元シード点を選定し;3次元シード点において形状モデルを適用し;且つ関心構造を表すように形状モデルを適応させ、それにより適応化された形状モデルを生じさせる;ように構成されている。
また、更なる他の典型的な一実施形態において、3次元(3D)画像の生成又は再構成のためのシステムが開示される。このシステムは:画像投影からの3次元画像の生成又は再構成を容易にするのに十分な角度範囲で、関心構造に対応する複数の画像投影を収集する手段;3次元シード点を選定する手段;3次元シード点において形状モデルを適用する手段;及び関心構造を表すように形状モデルを適応させ、それにより適応化された形状モデルを生じさせる手段;を含んでいる。
さらに、更なる他の典型的な一実施形態において、3次元(3D)画像を生成あるいは再構成する上述の方法をコンピュータに実行させる命令を含んだ機械読み取り可能なコンピュータプログラムコードで符号化された記憶媒体が開示される。
また、更なる他の典型的な一実施形態において、3次元(3D)画像を生成あるいは再構成する上述の方法をコンピュータに実行させる命令を有するコンピュータプログラムが開示される。
開示された方法に関する更なる特徴、機能及び利点は、以下の詳細な説明により、特にそれが添付図面と併せて精査されることにより明らかになる。
開示された実施形態を当業者が製造・使用することを支援するため、添付の図面を参照する。図面において、似通った参照は同様の符号を付されている。
ここで説明されるように、この開示は、特に電気生理的(electro physiology;EP)診療のための、とりわけ心室構造の、3次元(3D)回転X線撮像を可能且つ容易にするものである。さらに、この開示は、少ない数の投影からの形状モデルに基づく再構成を可能且つ容易にし、低ドーズ量の4D(例えば、心臓の位相を有する3D)X線再構成を実現するものである。
本発明は、3D/4D撮像の様々な種類の用途に使用され得る。本発明の好適な実施形態は、例示として、ここでは電気生理的診療に利用されるX線撮像に適用され得るものとして説明される。好適な実施形態が例として且つX線撮像及び診療を参照して示され且つ説明されるが、当業者に認識されるように、本発明はX線撮像又は診療のみに限定されるものではなく、複数の撮像システム及び用途に適用され得るものである。さらに、認識されるように、ここで開示される用途は診療のみに限定されるものではなく、実際には、一般に3D/4D撮像が望まれる如何なる用途にも適用され得るものである。
更に認識されるように、典型的な一実施形態を説明するために特定のセンサ及び用語が列挙されるが、このようなセンサは単に例示のために説明されており、限定的なものではない。ここでの開示を熟考する者には数多くの変形、代替及び均等物が明らかになるであろう。
典型的な一実施形態において、関心ある心室構造の3D回転X線データは、患者の心電図(ECG)の測定と同時に収集される。標的構造の起点(シード点)が選定され、このシードの周りに、患者の位置に適応された向きと好ましくは標的血管構造を十分に表すことが知られた形状とを有する適応型形状モデルが設けられる。投影データに表された実際の患者データに従って、形状モデルは複数の心臓位相に適応される。得られた4D心室モデルは診療誘導及び心室パラメータの評価に直接的に利用されてもよい。他の例では、形状表面についてのモデル化された3D運動を使用して局所的な運動ベクトル場を生成することが可能であり、この局所的な運動ベクトル場が、4Dデータセットの局所的運動補償された再構成をもたらす回転X線データの補償を実現するために使用されてもよい。
図1を参照するに、本発明の典型的な一実施形態に従ったシステムが示されている。システム10は、C字アーム14を具備したX線装置12を含んでおり、C字アーム14は第1の端部に配置されたX線管16と、その他方の端部に配置された例えば画像増強管であるX線検出器18とを有している。このX線装置12は、テーブル22上に配置された患者20の、相異なるX線位置からのX線投影画像を形成するのに適したものであり、この目的のため、C字アーム14の位置は様々な向きに変化させられることができ、C字アーム14はまた空間内の3つの軸、すなわち、図示されたX軸及びZ軸と図示されていないY軸との周りで回転可能なように構築されている。C字アーム14は、支持装置24、ピボット26及びレイルシステム30内で水平方向に移動可能なスライド28によって天井に取り付けられ得る。相異なるX線位置からの投影の収集のためのこれらの運動の制御、及びデータ収集の制御は、制御ユニット50によって行われる。
以下に限られないがプローブ、手術針、カテーテル、誘導ワイヤ、これらに類するもの、及びこれらの少なくとも1つを含む組み合わせ、を含む医療機器32が、例えば血管造影手順、生検、又は介入療法などの間に患者20に導入される。典型的な一実施形態に従って、患者20の検査領域の3次元画像データセットに対する医療機器32の位置が、位置測定システム(図示せず)を用いて収集・測定され、且つ/或いはここで説明されるように再構成された3D/4D画像上に重ね合わされる。
また、心電図(ECG)測定システム46がシステム10の一部としてX線装置12に備えられている。典型的な一実施形態において、ECG測定システム46は制御ユニット50にインターフェース接続されている。患者20のECGは、好ましくは、心臓位相の決定を容易にするようにX線データの収集中に測定且つ記録される。典型的な一実施形態において、心臓位相の情報はX線投影データを区分けし区別するために使用される。認識されるように、典型的な一実施形態はここでは心臓位相を解明するためのECG測定を参照して説明されているが、その他の手法も可能である。例えば、心臓位相及び/又は投影データの区分けは、X線データのみ、その他のパラメータ、又は更なる検知データに基づいて達成されてもよい。
制御ユニット50はX線装置12を制御し、画像捕捉を容易にし、画像再構成を促進するための機能及び処理を提供する。制御ユニット50は、算術ユニット52にて処理されるように収集データ(以下に限られないが、X線画像及び位置データなどを含む)を受信する。算術ユニット52はまた、制御ユニット50にインターフェース接続されており、それによって制御される。診療中に医師を支援するために様々な画像がモニター54上に表示され得る。
所定の機能、所望の処理、及びそれらのための計算(例えば、X線制御及び画像再構成など)を実行するため、制御ユニット50、算術ユニット52、モニター54及び再構成ユニット56等は、以下に限られないが、プロセッサ、コンピュータ、メモリ、記憶装置、レジスタ、タイミング、中断、通信インターフェース、入/出力信号インターフェース、これらに類するもの、及びこれらの少なくとも1つを有する組み合わせを含み得る。例えば、制御ユニット50、算術ユニット52、モニター54及び再構成ユニット56等は、X線投影の生成及びそれからの3D/4D画像の再構成を容易にするために必要な、X線信号の正確な標本化、変換、収集又は発生を可能にする信号インターフェースを含んでいてもよい。制御ユニット50、算術ユニット52、モニター54及び再構成ユニット56等の更なる特徴について、ここで十分に説明する。
図示されたX線装置12は、診療に先立って、且つ/或いは診療と同時に行われる典型的な一実施形態の例において、相異なるX線位置からの一連のX線投影画像を形成することに適したものである。X線投影画像から、3次元画像データセット、3次元再構成画像、及び望まれる場合にはそれらからのX線断層画像が生成され得る。収集された投影は、本発明の典型的な一実施形態に従って、この方法に適合した算術ユニット52に与えられ、その後、後述の運動補償に基づいて投影からそれぞれの再構成画像を形成する再構成ユニット56に与えられる。得られた3D画像はモニター54上に表示されることができる。最終的に、3次元画像データセット、3次元再構成画像、X線投影画像及びこれらに類するものは記憶ユニット58に保存・格納され得る。
続いて図2及び3を参照するに、図2は、開示された方法の典型的な一実施形態を例示するブロック図100を示している。図3は、心臓の図解に適用された典型的な実施形態の一例を示している。最初に、ブロック102にて示されるように、関心構造60(例えば、これに限られないが左心室を含む心室構造)の3D回転X線データが、画像投影からの3D画像の生成又は再構成を容易にするのに十分な角度範囲を有する軌道に沿って収集される。一実施形態において、少なくとも180°に扇角を足し合わせた値を有する撮像対象角度が用いられる。3次元X線データの収集を容易にするため、関心構造60内に含まれる血液ボリュームのコントラストが、参照符号62として大まかに示された造影剤によって強調される。造影剤62は、静脈内に投与されてもよいが、好ましくは、関心構造60が完全な回転収集に従って充たされるようにカテーテルを介して関心構造60に直接的に供給される。回転X線データの収集と並行して、患者20のECGが測定される。
ブロック104にて、図示されるように、関心構造60に対応する3D空間内の3Dシード点(例えば、モデルの開始点)64が選定される。典型的な一実施形態において、投影の最小強度の中心が用いられる。しかしながら、その他のシード点64及びその選定方法も可能である。
続いて図4をも参照するに、測定された投影の各々の上で3D空間内の3Dシード点64を選定するため、円錐ビーム形状に沿って2次元の線積分の組が測定される。例えば、3D空間内の造影剤で充たされた心室などの関心領域60に適用されるとき、この構造を貫く対応する線積分は、適切な較正後、検出器上に高い値で現れる。認識されるように、元々は強度が測定されるが、主ビームの強度が分かると、ランベルト・ベール(Lambert Beer)の法則を逆に用いることによって吸収係数を介した線積分が計算され得る。
m1と表記される投影内の例えば心室である関心構造60のおおよその中心は、その投影の最大の線積分を採ることによって、あるいは雑音を除去するために低域通過フィルタを用いて投影を畳み込んだ後に最大値を採ることによって、各投影に対して決定される。他の例では、投影面内で或る一定の形状を探索し、この形状内で、あるいは別の手法によって、線積分の重心を計算するセグメント化手法が適用されることができる。
各投影面内で関心構造60のおおよその中心(m1)が決定されると、同一の心臓位相に属する少なくとも3つ、好ましくはそれより多くの投影が、例えば10%RR間隔に属する投影群といった投影セットから選択される。この例ではθi、θj及びθkと表記されたゲーティングされた投影セットに含まれる、選択された“角度”又はこの投影に関する角度群に対応する1つの投影、例えばθi、が選択され、関心構造60の投影m1iの中心から放射線源までの放射線Siが採られる。同一位相に対応するその他全ての投影θj、θk、・・・から、対応する放射線Sj、Sk、・・・が生成される。例えばSj、Skであるその他の放射線から放射線Siへの最短距離が3D空間内で計算される。θiに関するこの放射線Si上の一組の点群d(i,k)、d(i,j)が得られ、加重和M1(i)が以下の等式:
M1(i)=(Σnd(i,n)w(i,n))/Σnw(i,n) (1)
w(i,n)=sin(θn−θi) (2)
に従って計算される。
その後、ゲーティングされた投影セットの例えばθjからの放射線Sjである第2の放射線に関してこの手順が繰り返される。θjに関するこの放射線Sj上の一組の点群d(j,k)、d(j,i)が得られ、加重和M1(j)が計算され、さらに、選択された全ての投影に関して続けられる。同一の心臓位相に属する投影群、故に対応する放射線群の各々に対し、各放射線上の3D点(M1(i)、M1(j)等)が得られる。最終的に、3Dシード点60が、やはり加重平均又は加重和の手法によってもたらされるこれらの点の“質量”中心として得られる。同一の手順は、心臓ゲーティングを用いるか否かに拘わらず、全ての投影に適用されることが可能である。しかしながら、心臓の運動のため、ゲーティングを用いないと結果が乱されることがある。有利には、上述の方法は心臓位相ごとに単一且つ固有のシード点60を生じさせる。
図2のブロック106を参照するに、3Dシード点の選定後、このシード点64の周りに適応型形状モデル66が適用される。好ましくは、必ずしもこうである必要はないが、適応型形状モデル66は、患者の位置に適応された向きと標的の関心構造60を十分に表すことが知られた形状とを有するように配置される。例えば、心臓診療のための血管への適用において、形状モデル66は好ましくは、例えば左心室などの撮像される血管構造と同様になるように位置付けられ、且つ成形される。
ブロック108に移り、典型的な一実施形態において、様々な投影データセット上への関心構造60の形状の順方向評価投影を生成するための適応処理の一部として、関心構造60に関する投影配置の正確な知識が用いられる。例えば、形状モデル66は、投影データに表される実際の患者データ、すなわち、関心構造60(例えば、心室)の厚さ及び吸収を3次元で表す該構造の境界及び線積分の値に基づいて、単一の心臓位相に適応される。
ここで図5A−5Cをも参照するに、典型的な一実施形態において、3D形状モデル66の適応化は、モデル66を複数の投影のうちの選択された数の投影に別々に適応させることによって達成される。必要に応じて、同様の運動の特徴を示す画像投影データの部分集合を取り扱うために、同一の心臓位相に属する全ての投影に同時に適応させることが用いられてもよい。3D適応型形状モデル66が、多数の接続線を有する形状の表面上に分布された多数の点を含んでいるとすると、適応化は図5A−5Cに示されるように系統立てられる。
最初に、3D適応型形状モデル66の表面点群が考慮中の投影面に順投影される。図5Aは例示的な形状モデルと関心ある投影上への順投影とを示している。検出面上でその投影内の点群の境界を定める3D表面点群が、図5Bに示されるように特定される。検出面内で隣接し合う境界点の間の接続が境界線を表している。3D境界線すなわち投影された3D境界線に対して垂直な方向にある境界点群の各々に対し、投影データ内の構造の境界線を表し得るエッジが探索される。境界点群はそれに従って図5Cに示されるように修正される。認識されるように、複数のゲーティングされた投影が利用可能である例においては、この適応化は好ましくは、同一の心臓位相に属する投影群の各々にて実行される。境界点群の新たな3D位置を決定した後、適応型形状モデルのその他の点群が、3D形状の所与の内部エネルギーの項に従って修正される。
形状モデル66を適応させる更なる情報として、考慮中の投影の方向に3D空間で形状モデル66を貫く線積分が考慮に入れられ得る。例えば、適応型3D形状モデル66を貫く線積分が計算されることができ、特定の投影における線積分の対応する2次元分布は、2D境界の修正に基づく最適3D形状の適応化を決定するように測定値を用いて補正されることができる。
この形状の適応化は、双方の測定(エッジ検出及び線積分分布)に基づく単一の適応化段階にて実行されることも可能であるし、反復的に実行されることも可能である。
認識されるように、形状モデル66の適応化を更に制約する、あるいは容易にする様々な既知の、あるいは推察可能な情報が有利に活用される。例えば、典型的な一実施形態において、適応化を容易にするため、既知である患者の向きは関心構造60の“もっともらしい”向きに対応する或る一定の論理的仮定又は“知識に裏付けられた推測”を可能にする。同様に、形状モデル66の適応化を更に容易にするために、個々の患者の関心構造に関する既知の情報が用いられてもよい。さらに、認識されるように、その他の心臓位相に関する同時又はその後の適応化は、場合により、隣接する位相における形状モデル又は関心構造60の形状についての知識を用いて行われてもよく、それにより形状モデル表面の連続的な運動をもたらすように形状変化を制限し得る。また、認識されるように、形状モデル66の一層正確な配置及び初期形状は、関心構造60の実際の形状とモデル化された形状との間の差を最小化すること、及び満足な結果を達成するために適応化に必要な繰り返しを削減することによって適応化を改善する。
最終的に、得られた4D心室モデルは、診療誘導を容易にするために直接的に使用されてもよく、また、当初のX線投影に基づいて適応された形状モデル66から生成された4D画像を提供することによって、心室パラメータの評価に用いられてもよい。この手法は図2にブロック110として示されている。他の例では、形状モデル表面の3D運動は、回転X線データの運動補償再構成をもたらす再構成処理において適用可能な局所的な運動ベクトル場を生成するために用いられることができる。それにより、全ての利用可能な投影が再構成処理中に、ある一定の参照状態に対して運動補償される。この手法は図2のブロック112にて示されている。
なお、開示された方法の動的でない部分は、静止構造の3Dモデルを生成するために使用されてもよい。また、認識されるように、ここで開示された技術は、回転データ収集中に周期運動によって形状が変化あるいは移動させられる如何なる用途にも容易に適用され得るものである。
開示された発明の更なる一実施形態は、必要とされる心室情報をモデル化によって生成するものである。このとき、同一の心臓位相内に相異なる投影方向から得られた複数の投影内で心室の輪郭が定められる。輪郭を付された心室構造は心室の3D形状の計算に使用される。この技術はまた、機能的情報を提供する4Dモデル化に拡張されることが可能である。
要するに、開示された発明は有利なことに、特に電気生理的(electro physiology;EP)診療のための、とりわけ心室構造の、3次元(3D)回転X線撮像を可能且つ容易にするものである。さらに、この開示は、少ない数の投影からの形状モデルに基づく再構成を可能且つ容易にし、低ドーズ量の4D(例えば、心臓の位相を有する3D)X線再構成を実現するものである。開示されたシステム及び方法は、電気生理的診療における誘導及びナビゲーションのために3D/4D再構成に頼る操作者、特に医師、に有意な利点を提供する。実際、開示されたシステム及び方法は、特に心サイクル内で生じる運動の補償を取り扱う3D/4D画像データのモデル化及び/又は再構成を提供する。開示されたシステム及び方法の更なる利点は、このモデル化及び再構成は、削減されたX線投影セットに基づいて実行されることができ、患者へのドーズ量の低減をもたらすことである。
明らかであるように、ここでは具体的に線積分、フィルタ、最大値を採ること、及び和を取ることとして参照された数学関数の実現には、技術的に数多くの数値計算方法が存在する。考え得る実現法が数多く存在しており、典型的な実施形態を例示するために用いられた具体的な実現法は限定的なものと見なされるべきではない。
以上の数多くの実施形態にて説明されたシステム及び方法は、特に心サイクル内で生じる運動の補償を取り扱う3D/4D画像データのモデル化及び/又は再構成のためのシステム及び方法を提供する。また、開示された発明は、コンピュータにより実現される処理、及びそのような処理を実行する装置の形態で具現化されてもよい。本発明はまた、例えばフロッピー(登録商標)ディスク、CD−ROM、ハードドライブ、またはその他の何らかのコンピュータ読み取り可能記憶媒体にて具現化された命令を含み、コンピュータにロードされ且つそれによって実行されるときに該コンピュータが本発明を実施する装置となるようなコンピュータプログラムコードの形態で具現化されることも可能である。本発明はまた、例えば、記憶媒体に格納され、コンピュータにロードされ且つ/或いはそれによって実行され、あるいは例えば電気配線若しくはケーブル上で、光ファイバを介して、あるいは電磁放射線を介してなど、何らかの伝送媒体上で変調された或いは変調されていない搬送波で伝送されるデータ信号として、コンピュータにロードされ且つそれによって実行されるときに該コンピュータが本発明を実施する装置となるようなコンピュータプログラムコードの形態で具現化されることも可能である。汎用マイクロプロセッサ上に実装されるとき、コンピュータプログラムコードのセグメント群は、特有の論理回路群を作り出すようにマイクロプロセッサを構成する。
認識されるように、同様の品目を表記するために“第1の”、“第2の”又はその他の同様の用語を用いていることは、特に断らない限り、何らかの特定の順序を指定したり伝えたりするものではない。同様に、“或る”(“a”若しくは“an”)又はその他の同様な用語を用いていることは、特に断らない限り、“1つ以上”を意味するものである。
本発明をその典型的な実施形態を参照して説明してきたが、当業者に理解されるように、この開示はそのような典型的な実施形態に限定されるものではなく、それらに含まれる要素は本発明の範囲を逸脱することなく様々に変更され、また均等物で代用され得るものである。また、具体的な状況又は材料を本発明による教示に適応させるために、本発明の基本精神又は範囲を逸脱することなく多様な変更、改善及び/又は変形が為され得る。故に、本発明は、それを実施するために意図されるベストモードとして開示された特定の実施形態に限定されるものではなく、添付の特許請求の範囲に入る全ての実施形態を含むものである。
本発明の典型的な一実施形態に従ったX線撮像システムを示す図である。 開示された方法の一例を示すブロック図である。 心臓の図解に適用される典型的な一実施形態に係る一例を示す図である。 本発明の典型的な一実施形態に従ったシード点を決定する方法を例示する図である。 例示的な形状モデルと関心ある投影上への順投影とを示す図である。 典型的な一実施形態に従った境界決定を示す図である。 典型的な一実施形態に従った境界点の修正を示す図である。

Claims (29)

  1. 関心構造に対応する画像投影データに基づいて3次元画像を生成あるいは再構成する方法であって:
    関心構造に対応する複数の画像投影を収集する収集段階であり、該画像投影からの3次元画像の生成又は再構成を容易にするのに十分な角度範囲で収集する収集段階;
    3次元シード点を選定する選定段階;
    前記3次元シード点において形状モデルを適用する適用段階;及び
    前記関心構造を表すように前記形状モデルを適応させる適応化段階であり、それにより適応化された形状モデルを生じさせる適応化段階;
    を有する方法。
  2. 前記適応化された形状モデルに基づいて前記関心構造の前記3次元画像を生成する生成段階;
    を更に有する請求項1に記載の方法。
  3. 前記生成段階は、前記関心構造の選択された運動位相に基づく、請求項2に記載の方法。
  4. 前記運動位相は心臓位相に相当する、請求項3に記載の方法。
  5. 前記運動位相を指し示すデータは心電図である、請求項4に記載の方法。
  6. 前記画像投影に関連付けて前記関心構造の運動を指し示すデータを収集する段階;
    前記関心構造の運動を指し示す前記データに基づいて、前記関心構造を表すように前記形状モデルを適応させる段階であり、それにより別の適応化された形状モデルを生じさせる段階;及び
    前記別の適応化された形状モデルに基づいて前記関心構造の前記3次元画像を生成する段階;
    を更に有する請求項1に記載の方法。
  7. 前記適応化された形状モデル及び前記関心構造の運動に基づいて運動ベクトル場を生成する段階;及び
    再構成を容易にするために前記運動ベクトル場に基づいて前記複数の画像投影を補償する段階;
    を更に有する請求項1に記載の方法。
  8. 前記運動ベクトル場を生成する段階は、前記関心構造の選択された運動位相に基づく、請求項7に記載の方法。
  9. 運動を指し示すデータを収集する段階が、前記複数の画像投影を収集する前記収集段階と実質的に同時に行われる、請求項1に記載の方法。
  10. 前記関心構造は冠状血管又は心室である、請求項1に記載の方法。
  11. 前記3次元シード点は、前記複数の画像投影のうちの選択された組の最小強度の重心に対応する、請求項1に記載の方法。
  12. 前記適応化段階は前記複数の画像投影の部分集合に基づく、請求項1に記載の方法。
  13. 前記角度範囲は、少なくとも約180°に扇角を足し合わせた角度に相当する、請求項1に記載の方法。
  14. 3次元画像の生成及び再構成のためのシステムであって:
    関心構造に対応する画像投影データを提供するように構成された撮像システムであり、該画像投影からの前記3次元画像の生成又は再構成を容易にするのに十分な角度範囲で画像投影データを提供する撮像システム;及び
    前記撮像システムと通信することが可能なコントローラであり:
    前記画像投影データを受信し;
    3次元シード点を選定し;
    前記3次元シード点において形状モデルを適用し;且つ
    前記関心構造を表すように前記形状モデルを適応させ、それにより適応化された形状モデルを生じさせる;
    ように構成されたコントローラ;
    を有するシステム。
  15. 前記コントローラは更に、前記適応化された形状モデルに基づいて前記関心構造の前記3次元画像を生成するように構成されている、請求項14に記載のシステム。
  16. 前記3次元画像は更に、前記関心構造の選択された運動位相に基づく、請求項15に記載のシステム。
  17. 前記運動位相は心臓位相に相当する、請求項16に記載のシステム。
  18. 前記運動位相を指し示すデータは心電図である、請求項17に記載のシステム。
  19. 前記適応化された形状モデル及び前記関心構造の運動の少なくとも一方に基づいて運動ベクトル場を生成し;且つ
    再構成を容易にするために前記運動ベクトル場に基づいて前記複数の画像投影を補償する;
    ように構成されている、請求項14に記載のシステム。
  20. 前記運動ベクトル場は、前記関心構造の選択された運動位相に基づく、請求項19に記載のシステム。
  21. 前記コントローラ又は前記撮像システムに動作的に接続され、前記画像投影データに関連付けて前記関心構造の運動を指し示すデータを提供するように構成された測定システム;
    を更に含み;
    前記コントローラは更に、前記関心構造の運動を指し示す前記データに基づいて、前記関心構造を表すように前記形状モデルを適応させ、それにより別の適応化された形状モデルを生じさせるように構成されており;且つ
    前記関心構造の前記3次元画像を生成することは、前記別の適応化された形状モデルに基づく;
    請求項14に記載のシステム。
  22. 運動を指し示す前記データは、前記複数の画像投影と実質的に同時に収集される、請求項21に記載のシステム。
  23. 前記関心構造は冠状血管又は心室である、請求項14に記載のシステム。
  24. 前記3次元シード点は、前記複数の画像投影のうちの選択された組の最小強度の重心に対応する、請求項14に記載のシステム。
  25. 前記撮像システムはX線撮像システムである、請求項14に記載のシステム。
  26. 前記角度範囲は、少なくとも約180°に扇角を足し合わせた角度に相当する、請求項14に記載のシステム。
  27. 3次元画像の生成又は再構成のためのシステムであって:
    関心構造に対応する複数の画像投影を収集する収集手段であり、該画像投影からの前記3次元画像の生成又は再構成を容易にするのに十分な角度範囲で収集する収集手段;
    3次元シード点を選定する選定手段;
    前記3次元シード点において形状モデルを適用する適用手段;及び
    前記関心構造を表すように前記形状モデルを適応させる適応化手段であり、それにより適応化された形状モデルを生じさせる適応化手段;
    を有するシステム。
  28. 3次元画像を生成あるいは再構成する方法をコンピュータに実行させる命令を含んだ機械読み取り可能なコンピュータプログラムコードで符号化された記憶媒体であって、該方法は:
    関心構造に対応する複数の画像投影を収集する収集段階であり、該画像投影からの3次元画像の生成又は再構成を容易にするのに十分な角度範囲で収集する収集段階;
    3次元シード点を選定する選定段階;
    前記3次元シード点において形状モデルを適用する適用段階;及び
    前記関心構造を表すように前記形状モデルを適応させる適応化段階であり、それにより適応化された形状モデルを生じさせる適応化段階;
    を有する、記憶媒体。
  29. 3次元画像を生成あるいは再構成する方法をコンピュータに実行させる命令を有するコンピュータプログラムであって、該方法は:
    関心構造に対応する複数の画像投影を収集する収集段階であり、該画像投影からの3次元画像の生成又は再構成を容易にするのに十分な角度範囲で収集する収集段階;
    3次元シード点を選定する選定段階;
    前記3次元シード点において形状モデルを適用する適用段階;及び
    前記関心構造を表すように前記形状モデルを適応させる適応化段階であり、それにより適応化された形状モデルを生じさせる適応化段階;
    を有する、コンピュータプログラム。
JP2008524630A 2005-08-04 2006-07-12 3d−2d適応型形状モデル支援による運動補償再構成 Expired - Fee Related JP5244592B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US70561605P 2005-08-04 2005-08-04
US60/705,616 2005-08-04
PCT/IB2006/052375 WO2007015181A1 (en) 2005-08-04 2006-07-12 3d-2d adaptive shape model supported motion compensated reconstruction

Publications (2)

Publication Number Publication Date
JP2009502397A true JP2009502397A (ja) 2009-01-29
JP5244592B2 JP5244592B2 (ja) 2013-07-24

Family

ID=37461481

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008524630A Expired - Fee Related JP5244592B2 (ja) 2005-08-04 2006-07-12 3d−2d適応型形状モデル支援による運動補償再構成

Country Status (8)

Country Link
US (1) US8446408B2 (ja)
EP (1) EP1913558B1 (ja)
JP (1) JP5244592B2 (ja)
KR (1) KR20080031358A (ja)
CN (1) CN101238488B (ja)
AT (1) ATE551677T1 (ja)
CA (1) CA2617382A1 (ja)
WO (1) WO2007015181A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011527584A (ja) * 2008-05-09 2011-11-04 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 運動物体の画像を生成するための装置

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102007016902B4 (de) * 2007-04-10 2014-05-28 Biosense Webster, Inc. Verfahren und Vorrichtung zur Gewinnung eines Volumendatensatzes von einem bewegten Gewebe oder Organ eines Patienten
JP5296057B2 (ja) * 2007-05-08 2013-09-25 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 低線量カルシウムスコアリングスキャンを使用する冠動脈選択的カルシウム割り当て
EP2254655B1 (en) 2008-01-31 2014-09-03 Medtronic, Inc. Automated programming of electrical stimulation electrodes using post-implant imaging
US8995731B2 (en) * 2008-11-26 2015-03-31 Medtronic, Inc. Image-based characterization of implanted medical leads
WO2010108146A2 (en) 2009-03-20 2010-09-23 Orthoscan Incorporated Moveable imaging apparatus
US20110093051A1 (en) * 2009-10-21 2011-04-21 Medtronic, Inc. Assignment and manipulation of implantable leads in different anatomical regions with image background
CZ2010225A3 (cs) * 2010-03-26 2011-07-20 Hanzelka@Tomáš Zpusob zvýšení presnosti rekonstrukce 3D rentgenového obrazu
CZ302594B6 (cs) * 2010-04-09 2011-07-27 Univerzita Tomáše Bati ve Zlíne Zpusob replikace povrchových struktur
US9814885B2 (en) 2010-04-27 2017-11-14 Medtronic, Inc. Stimulation electrode selection
RU2595757C2 (ru) * 2010-12-02 2016-08-27 Конинклейке Филипс Электроникс Н.В. Устройство совмещения изображений
US9125611B2 (en) 2010-12-13 2015-09-08 Orthoscan, Inc. Mobile fluoroscopic imaging system
US8406890B2 (en) 2011-04-14 2013-03-26 Medtronic, Inc. Implantable medical devices storing graphics processing data
CN103533892B (zh) * 2011-05-12 2016-08-17 皇家飞利浦有限公司 运动补偿成像
US10737116B2 (en) * 2013-08-07 2020-08-11 Koninklijke Philips N.V. Therapy planning
US9959631B2 (en) 2014-02-21 2018-05-01 Samsung Electronics Co., Ltd. Tomography apparatus and method for reconstructing tomography image thereof
WO2017202712A1 (en) * 2016-05-24 2017-11-30 Koninklijke Philips N.V. Depth-enhanced tomosynthesis reconstruction
US11367193B2 (en) * 2019-12-31 2022-06-21 Uih America, Inc. Systems and methods for medical imaging

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004033749A (ja) * 2002-04-15 2004-02-05 General Electric Co <Ge> Pet腫瘍画像に関する半自動セグメント分割アルゴリズム
JP2006527619A (ja) * 2003-06-16 2006-12-07 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 時系列イメージにおけるイメージセグメンテーション

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1149388C (zh) * 2001-02-23 2004-05-12 清华大学 一种基于移相法的数字投影三维轮廓重构方法
DE10130645A1 (de) * 2001-06-27 2003-01-16 Claas Selbstfahr Erntemasch Stützradanordnung für eine landwirtschaftliche Arbeitsmaschine
KR20040076854A (ko) * 2001-10-24 2004-09-03 뉴어오케이 엘엘씨 3차원 영상들의 투사
US6904118B2 (en) * 2002-07-23 2005-06-07 General Electric Company Method and apparatus for generating a density map using dual-energy CT
US7289841B2 (en) * 2002-10-25 2007-10-30 Koninklijke Philips Electronics N.V. Method and apparatus for volumetric cardiac computed tomography imaging
EP1593093B1 (en) * 2003-01-31 2006-08-16 Philips Intellectual Property & Standards GmbH Method for the reconstruction of three-dimensional objects
WO2005020155A1 (en) * 2003-08-21 2005-03-03 Philips Intellectual Property & Standards Gmbh Device and method for generating a three-dimensional vascular model

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004033749A (ja) * 2002-04-15 2004-02-05 General Electric Co <Ge> Pet腫瘍画像に関する半自動セグメント分割アルゴリズム
JP2006527619A (ja) * 2003-06-16 2006-12-07 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 時系列イメージにおけるイメージセグメンテーション

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JPN6012003655; Vidya Elangovan and Ross T. Whitaker: 'From Sinograms to Surfaces: A Direct Approach To The Segmentation Of Tomographic Data' MICCAI '01 Proceedings of the 4th International Conference on Medical Image Computing and Computer-A , 200110, p.213-223, Springer-Verlag Lo&#xFF4E *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011527584A (ja) * 2008-05-09 2011-11-04 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 運動物体の画像を生成するための装置

Also Published As

Publication number Publication date
ATE551677T1 (de) 2012-04-15
WO2007015181A1 (en) 2007-02-08
US8446408B2 (en) 2013-05-21
CN101238488A (zh) 2008-08-06
EP1913558A1 (en) 2008-04-23
JP5244592B2 (ja) 2013-07-24
CA2617382A1 (en) 2007-02-08
EP1913558B1 (en) 2012-03-28
US20080218510A1 (en) 2008-09-11
KR20080031358A (ko) 2008-04-08
CN101238488B (zh) 2010-08-04

Similar Documents

Publication Publication Date Title
JP5244592B2 (ja) 3d−2d適応型形状モデル支援による運動補償再構成
US7180976B2 (en) Rotational angiography based hybrid 3-D reconstruction of coronary arterial structure
US8184883B2 (en) Motion compensated CT reconstruction of high contrast objects
EP2467832B1 (en) System and method for four dimensional angiography and fluoroscopy
EP1606770B1 (en) Motion-corrected three-dimensional volume imaging method
US9633454B2 (en) Fluid-dynamic analysis of a vascular tree using angiography
US20100189337A1 (en) Method for acquiring 3-dimensional images of coronary vessels, particularly of coronary veins
JP2011507584A (ja) 身体及び体輪郭の3次元再構成
JP2006516440A (ja) 三次元物体の再構築法
US20230097849A1 (en) Creation method of trained model, image generation method, and image processing device
US8855391B2 (en) Operating method for an imaging system for the time-resolved mapping of an iteratively moving examination object
KR20170105876A (ko) 단층 촬영 장치 및 그에 따른 단층 영상 재구성 방법
CN107886554B (zh) 流数据的重构
US20090238412A1 (en) Local motion compensated reconstruction of stenosis
US7116808B2 (en) Method for producing an image sequence from volume datasets
Fahrig et al. 14 C-arm CT in the interventional suite: Current status and future directions
Lehmann et al. Angle‐independent measure of motion for image‐based gating in 3D coronary angiography
US10463334B2 (en) System and method for non-invasive, quantitative measurements of blood flow parameters in vascular networks
Zhu et al. Cone-Beam CT for Vascular Imaging
Giordano Perfusion imaging in the peripheral vasculature using interventional C-arm systems

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20090709

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120131

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20120425

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20120507

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120731

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120821

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20121119

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20121127

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130220

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130408

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

Free format text: PAYMENT UNTIL: 20160412

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20160412

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees