JP2008531108A - カテーテルの経路の予測のための方法 - Google Patents

カテーテルの経路の予測のための方法 Download PDF

Info

Publication number
JP2008531108A
JP2008531108A JP2007556700A JP2007556700A JP2008531108A JP 2008531108 A JP2008531108 A JP 2008531108A JP 2007556700 A JP2007556700 A JP 2007556700A JP 2007556700 A JP2007556700 A JP 2007556700A JP 2008531108 A JP2008531108 A JP 2008531108A
Authority
JP
Japan
Prior art keywords
tube
microcatheter
corridor
catheter
probe
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.)
Pending
Application number
JP2007556700A
Other languages
English (en)
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 JP2008531108A publication Critical patent/JP2008531108A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G09EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
    • G09BEDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
    • G09B23/00Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes
    • G09B23/28Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes for medicine
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass
    • GPHYSICS
    • G09EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
    • G09BEDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
    • G09B23/00Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes
    • G09B23/28Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes for medicine
    • G09B23/285Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes for medicine for injections, endoscopy, bronchoscopy, sigmoidscopy, insertion of contraceptive devices or enemas
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10Computer-aided planning, simulation or modelling of surgical operations
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M25/00Catheters; Hollow probes
    • A61M25/01Introducing, guiding, advancing, emplacing or holding catheters

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Medicinal Chemistry (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Educational Technology (AREA)
  • Educational Administration (AREA)
  • Pure & Applied Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Biomedical Technology (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Radiology & Medical Imaging (AREA)
  • Epidemiology (AREA)
  • Pulmonology (AREA)
  • Primary Health Care (AREA)
  • Pathology (AREA)
  • Media Introduction/Drainage Providing Device (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本発明は、血管系における開始位置と目的位置との間のカテーテルの経路の予測のための方法に関する。好適な実施形態にしたがって、マイクロカテーテルは、血管系を通るマイクロカテーテル中心線(MC)を追従するマイクロカテーテル管(MT)によりモデル化され、前記中心線は、直線部分と曲線部分との交互の連続を有する。曲線部分は、マイクロカテーテル管が血管壁と接する及び/又は血管系の側枝に入る位置に導入される。

Description

本発明は、血管系内の(マイクロ)カテーテルの経路の予測のための方法、カテーテルの製造のための方法、予測方法の実行のためのデータ処理ユニット及び予測方法を実行するプログラムを有する記録担体に関する。
患者の血管系内のカテーテルの介入の成功性は、カテーテルの介入が適切なモデル化手順を得て計画され且つ準備される場合に、改善される。前記手順は、典型的には、患者の血管系の三次元の幾何学的モデルを必要とし、例えば、三次元回転血管造影により取得されることが可能である。カテーテルの介入の典型的な例は動脈瘤の治療であり、モデル化三次元血管系における動脈瘤のボクセルの完全自動化のための方法については、文献(“Fully−automatic labelling of aneurysm voxels for volume estimation”,by J.Bruijns,Proc.Bildverarbeitung fuer die Medizin,pages 51−55,Erlargen,Germany,March 2003参照)に記載されている。動脈瘤のラベリングの後、適切な品質、例えば、直径及び弾性を有する化テータルの選択を含む治療計画が、医師によりなされる必要がある。
"Semi−automatic shape extraction from tube−like geometry",by J.Bruijns,Proc.VMV,pages347−355,Saarbruecken,Germany,November 2000 "Fully−automatic branch labelling of voxel vessel structures",by J.Bruijns,Proc.VMV,pages341−350,Stuttgart,Germany,November 2001 "Fully−automatic labelling of aneurysm voxels for volume estimation",by J.Bruijns,In Proc.BVM,pages51−55,Erlangen,Germany,March 2001 "3−d−rotational angiography:First clinical application with use of a standard philips c−arm system",by R.Kemkers,J.Op de Beek,H.Aerts,R.Koppe,E.Klotz,M.Grasse,and J.Moret,In Proc.CAR,Tokyo,June 1998 "3d rotational angiography:Clinical value in endovascular treatment",by J.Moret,R.Kemkers,J.Op de Beek,R.Koppe,and M.Grass,Medicamundi,42(3):8−14,November 1998 Philips Medical Systems Nederland,Integris 3d−ra.instructions for use.release2.2.Technical Report 9896 001 32943,Philips Medical Systems Nederland,Best, The Netherlands, January 2001 "Pondering on discrete smooth interpolation",by C.W.A.M.van Overveld,Computer−aided Design,27(5):377−384,November 1995
この状況に基づいて、本発明の目的は、カテーテルの介入についての支援された且つ改善された計画のための手段を提供することである。
この目的は、請求項1に記載の方法、請求項9に記載のデータ処理ユニット及び請求項11に記載の記録担体により達成される。
第1の特徴にしたがって、本発明は、モデル化血管系における所定の開始位置(例えば、カテーテルが体内に導入された切開)と所定の目的位置(例えば、動脈瘤)との間のカテーテルの経路の予測のための方法に関する。用語“カテーテル”は、原理的に、患者の血管系を通って進むことができる何れかの長方形の器具を有する。カテーテルの経路は、“経路管”と呼ばれる管状物体により表され、前記管は、開始位置から目的位置まで繋が関連“経路中心線”に沿って進む。この方法は、次のステップを有する。
a)開始位置から目的位置まで繋がる血管系を通る経路の決定及び最初の経路中心線の前記経路との同一化。血管系は、例えば、中心線を有する管状物体によりモデル化される場合、その経路は前記血管中心線を追従することが可能である。
b)この中心線と関連する経路管が血管系内にあるように、上記の最初の経路中心線の調節。好適には、結果として得られる経路管は、更なる(最適)基準を満たし、また、例えば、屈曲エネルギーを最小化する構成を有する。
上記のような種類の方法により、前記血管系に導入されるカテーテルの経路が血管系について個々に予測されることが可能であるため、血管系における介入の計画を改善する 及び自動化することが可能である。このことは、医師が介入の実行可能性に関して判定することを及び最適な方法が実行されることを支援する。更に、その方法は、特定のタスクのためにカテーテルの最適なタイプを選択し、準備することが可能であるように支援する。
提案された方法の第1のアプリケーションにおいては、経路管は、カテーテルが開始位置から目的位置まで血管系を通って進むことが可能であるコリドール(corridor)を表す“コリドール管”であることが可能である。介入は、残りの血液流のために空間を尚も残しつつ、コリドールがカテーテルを受け入れるのに十分大きい場合に、実行可能である。
提案された方法の第2のアプリケーションにおいては、経路管は、開始位置から目的位置まで血管系を通って進むマイクロカテーテルの形状を表す“マイクロカテーテル管”であることが可能である。表記“マイクロカテーテル”は、このアプリケーションが小さく細長いカテーテルについて、特に適切であることを示す。しかしながら、その用語は制限するように意図されているのではなく、患者の血管系を通って進められる何れかの長方形の器具を原理的に有するものである。
上記のアプリケーションの組合せにおいては、コリドール管が先ず、決定され、マイクロカテーテル管が、前記コリドール管内にあるように、次に決定される。マイクロカテーテル管の決定において、上記方法のステップa)において必要な血管系を通る経路は、好適には、コリドール管の中心線として規定される。
本発明の好適な実施形態においては、マイクロカテーテルの中心線は、交互の直線部分及び曲線部分を有する。直線部分については、関連管部分は血管系の内側において規定されたところの傍(即ち、血管系の壁から遠い距離のいたるところ)にあり、それとは対照的に、曲線部分については、関連管部分は血管壁に接触する(血管系の周りの組織に入り込むことなく)及び/又は血管系の分枝の方に曲がる。それらの名前で示されるように、直線部分は、好適には、(略)直線的である一方、曲線部分は曲がっている。一続きの直線及び曲線部分は、マイクロカテーテルが血管壁に接触する又は血管系の枝に入らなければならなくなるまで、まっすぐに伸びる、薄く、細長いマイクロカテーテルの記述について特に適切である。
上記の一続きの直線及び曲線部分は、例えば、直線位置における直線部分から始まり、特に繰り返す方式で決定される。この繰り返しの間、直線部分は、その場合、曲線部分の導入が血管系の内部にマイクロカテーテルを戻す又は分枝に入ることが必要になるまで、伸ばされる。
上記方法の好適な実施形態においては、各々の繰り返しステップは次のステップを有する。
aa)“カテーテルコーナー”の決定。前記カテーテルコーナーは、i)繰り返しステップにおいて現在考慮されている直線部分を囲む血管壁との交差部分として、又はii)マイクロカテーテルが追従する側枝の最も遠い血管壁と、前記部分の始まりから同じ距離にある現在の直線部分における点として規定される(選択肢(i)、(ii)のどちらかが近い)。カテーテルコーナーは、それ故、現在の直線部分の直線状経路が終了する必要がある血管系の点を示す。
bb)ステップaa)において決定されたカテーテルコーナーの方への“最初の移動ベクトル”だけカテーテルコーナーの方に近い現在の直線部分の点の移動。典型的には、この点は、境界条件(例えば、マイクロカテーテル管が血管系の内側に留まる必要性)に対してできるだけ近いカテーテルコーナーまで移動される。
cc)現在の直線部分の上記のように移動される点において現在の直線部分から後続の曲線部分への移行の導入。曲線部分の経路は、その場合、所定の境界条件にしたがって決定される必要がある。
上記の繰り返す方法の好適な状態においては、ステップcc)において導入される後続の曲線部分は連続していて、ステップbb)の初期移動ベクトルの方向に少しずつ移動され、移動する長さは、関連マイクロカテーテル管がその血管系の壁を貫くことなく血管系の壁に単に接するように単調に減少する。移動長の単調な減少はマイクロカテーテル管の局所的メアンダーを回避する。更に、後続の直線部分は、現在の移動長だけ移動されるマイクロカテーテル管が初めて血管壁への接触を失う点から始める。換言すれば、マイクロカテーテルの中心線は血管壁により必要とされるのと同程度に移動され、マイクロカテーテル管が再び、血管系の内部で自由に進むことができる、次の直線部分に変化する。
上記方法において用いられる管又は管状物体(例えば、マイクロカテーテル管又はコリドール管)は、好適には、一連のプローブで表され、各々のプローブは中心及び関連平面を有する球を有する。前記球の中心はモデル化管の中心線上にあり、関連平面は前記中心を有し、その管の中心線に対して垂直になっている。更に、プローブは、更なるパラメータ、例えば、その管の断面に対応する楕円形断面の半径により特徴付けられることが可能である。
本発明は更に、カテーテル、好適には、マイクロカテーテルの製造のための方法であって、次のステップを有する方法に関する。
a)意図された上記の種類の方法により意図された繰り返しの間のカテーテルの経路の予測。
b)予測経路にしたがった、カテーテルの準備、好適には、前モデル化。
この方法により、(マイクロ)カテーテルは、特定の介入及び特定の患者について個別にデザインされることが可能である。これは、実質的に介入を容易にし、異なるケースの治療を可能にし、複雑化のリスクを低減する。
本発明は、上記の種類の予測方法を実行するように適合されるデータ処理ユニットに更に関連する。データ処理ユニットは、中央演算処理装置、記憶装置、I/Oインターフェースのような通常のコンピュータ構成要素及び関連コンピュータプログラムを伴うものを有することが可能である。
最終的に、本発明は、記録担体、例えば、フロッピー(登録商標)ディスク、ハードディスク又はコンパクトディスク(CD)を有し、それらの記憶担体において、上記の種類の方法にしたがったカテーテルの経路の予測のためのコンピュータプログラムが記憶される。
本発明の上記の及び他の特徴は、上記の実施形態から明らかであり、そして以下、詳述する実施形態に関連して理解することができる。
以下、本発明について、添付図の助けにより例示として説明する。
以下の図及び本発明の好適な実施形態においては、本発明者による文献を引用する。
1.序論
三次元回転血管造影法(参照文献4、5を参照されたい)により得られる血管のボリューム表現は、組織(組織は、血管以外の全て)と血管ボクセルとの間の濃淡値(濃淡値はX線吸収量を表す)における明確な違いを有する。それ故、それらのボクセルの表現は、動脈瘤の診断、血管の局所的で全方位的な広がりについて非常に適する(図5.1を参照されたい)。参照文献3において、動脈瘤ボクセルの完全自動化ラベリングの方法について記載されている(図5.2を参照されたい)。
ここで好適に用いられるモデル化血管系は、次のような構成要素を有する(参照文献1、2及び3を参照されたい)。
1.周期的な三次元グリッドの各々の点について、この点が血管に属するかどうか、その点が“標準的な”血管点か動脈瘤における点である血管点の場合、及びその点がどの分岐点又は合流点(“分岐”)に属するかの“標準的な”血管点の場合の表示を有する三次元ボリュームモデル。
2.血管と非血管との間の境界を表す表面モデル。この表面モデルの各々の頂点は、位置ばかりでなく、頂点が動脈瘤の境界の一部であるかどうか又はその頂点が分岐点又は合流点に属するかを示すnormal及びlabelを有する必要がある。
3.合流点と分岐点との間の関係を表すグラフ。
動脈瘤がラベリングされた後に、次のステップは治療計画を作成するようになっている。医師は、先ず、動脈瘤の内部にカテーテルを移動させ、次に、動脈瘤内にカテーテルを通してコイル又はグルーを注入することにより、その動脈瘤を治療することが可能である。“コリドール管”によりカテーテルについて血管を通るコリドールをモデル化する(図5.5を参照されたい)。コリドール管の中心線はコリドール血管を通る中心線を表す。コリドール管の直径はコリドール血管の直径を表す。コリドール管は、正しい品質(例えば、直径、弾性)を有するカテーテルを選択するように用いられることが可能である。コリドール管の演算については、第2節で説明される。
動脈瘤がカテーテルにより満たされる前に、マイクロカテーテルは、動脈瘤中に血管を通して移動される。マイクロカテーテルは、血管に比べて非常に細長い物体である。それ故、コリドール管の中心線はマイクロカテーテルの中心曲線とは異なる。実際には、マイクロカテーテルは、血管壁に多かれ少なかれ追従し、血管が折れ曲がるときに交差する(図5.3を図5.4と比較する)。マイクロカテーテルが、動脈瘤中への容易な移動のために選択され、前モデル化されたために、コリドール管からマイクロカテーテルの形状を演算する方法が開発された。マイクロカテーテルの形状の演算については、第3節で説明される。
2.コリドール管
医師は、先ず、動脈瘤内部にカテーテルを移動させ、次に、動脈瘤中にカテーテルを通してコイル又はグルーを注入することにより動脈瘤を治療することが可能である。“コリドール管”によりカテーテルについて血管を通るコリドールがモデル化される。コリドール管の中心曲線はコリドール血管を通る中心曲線を表す。コリドール管の直径はコリドール血管の直径を表す。コリドール管は、正しい品質(例えば、直径、弾性)を有するカテーテルを選択するように用いられることが可能である。実際には、コリドール管の最も小さい直径はカテーテルについての上限を与える。最も小さいコリドール管と選択されたカテーテルとの間の最も小さい断面積の差は、残りの流量能力を表す。コリドール管により動脈瘤を満たすことはシミュレートされないことに留意されたい。
本発明のシステムにおいては、管状物体(省略して管という)は一連のプローブを有する(参照文献1を参照されたい)。プローブは、球、球の中心を通る平面及び複数の形状パラメータの組合せである。管が完全自動化された血管トレーシング(参照文献2を参照されたい)により形成される場合、各々のプローブの球の中心は血管の中心軸の近傍にあり、各々のプローブの平面は血管に対して略垂直であり、各々のプローブの形状パラメータは局所的断面を近似する楕円を含む。管の表面の近似表現としてプローブの楕円を用いることができる。
コリドール管は2つの部分、即ち、血管及び延長管を有する。血管は“通常の”血管部分を通るコリドールを表す。延長管は。血管の端部から動脈瘤へのコリドールを表す。
コリドール管の始点及び終点は2つのプローブを開始することにより生成される。先ず、ユーザは、動脈瘤に接続された“通常の”血管部分の表面の二次元画像における点を選択する。次に、本発明のシステムは、選択された表面の点を通るビュー光線に対して最近接の中心軸上の血管ボクセルの方に第1プローブを移動させる。その後、ユーザは、動脈瘤の表面において点を選択する。この第2の点を通るビュー光線は、動脈瘤の前と後との間の線部分を規定する。第2プローブは、この線部分の中心に対して最近接の血管ボクセルの方に移動される。開始点及び終了点が選択された後に、コリドール管が次のアルゴリズムにより生成される。
1.第1プローブに対して最近接の(即ち、“標準的な”血管部分に沿って最も短い経路を有する)動脈瘤のネックを求める。動脈瘤のネックは、動脈瘤と“標準的な”血管部分との間の接続であり、例えば、各々の“標準的な”血管ボクセルは少なくとも1つの動脈瘤ボクセルに接続された面である“標準的な”血管ボクセル(“ネックボクセル”という)の接続された集合によりモデル化されることが可能である。即ち、2つ又はそれ以上のネックボクセルの離散的に接続された集合が存在する場合、動脈瘤は2つ以上のネックを有することが可能である。
2.第1プローブからこのネックの中心まで完全自動化血管トレーシング(参照文献2を参照されたい)により血管を生成する。この血管は次のように規定される。
(a)血管の中心線(即ち、プローブの球の中心)は条件付き緩和によりスムージングされる(参照文献7を参照されたい)。
(b)各々の楕円はそれらの楕円と同じ面積を有する円と置き換えられる。
(c)それらの円の半径は、それらの半径のスムーズな(例えば、最小二乗)近似により置き換えられる。所定の近似関数の値によりプローブ数の関数として強く変化する可能性のある半径の集合を置き換えることは、よりスムーズに変化する半径の集合をもたらす。一次関数、三次関数、スプライン関数等がオリジナルのデータを近似するように用いられることが可能である。
3.ネック中心から第2プローブに延長管を生成する。延長管の中心線を生成するように二次ベジェ曲線を用いる。このベジェ曲線は、ネック中心、第2プローブの位置及び動脈瘤中心とネック中心との間の正規化方向により規定される。残りの自由度は、制御多角形の2つの辺は等しい長さを有するという任意の、しかし妥当な制約により除かれる。楕円(プローブの形状パラメータ)の半径は、血管の最終の楕円(即ち、円)の半径に等しい。
4.コリドール管は、血管及び延長管の連結である。このコリドール管はまた、血管について(a)乃至(c)で記載されているのと同様な方法で改善される。
コリドール管(即ち、その表面)の例について、図5.5に示している。
3.マイクロカテーテル
マイクロカテーテルの形状は管状物体により表される。マイクロカテーテルはカテーテルと同じ血管を通るコリドールに追従するために、そしてこのコリドールはコリドール管により表されるために、マイクロカテーテル管は、マイクロカテーテルの半径により置き換えられる全ての半径を有するコリドール管をコピーすることにより初期設定される。
マイクロカテーテルの(及びマイクロカテーテル管の)最終的な中心曲線は、交互の直性部分及び曲線部分を有する。マイクロカテーテルの剛性によりもたらされる直線部分は、マイクロカテーテルがもはや、血管壁により曲げられないところで始まる。曲線部分は、直線部分が血管壁と衝突するところ(図1の左側のピクチャ)で、又はマイクロカテーテルが側枝に追従するところ(図1の右側のピクチャ)で始まる。
マイクロカテーテル管のこの最終的な中心曲線は、次のように、反復アルゴリズムにおいてマイクロカテーテルのプローブに一連の移動ベクトルを適用することにより演算される。
1.マイクロカテーテル管の開始位置及び方向に対して次の直線部分の開始位置及び方向を設定する(通常は、動脈瘤から最も遠く離れた点)。
2.次の直線部分の開始位置及び方向に有効な初期移動及び更新を適用する。
3.新しい開始位置及び方向が求められる一方、
(a)次の直線部分は現在の直線部分になる。
(b)現在の直線部分の次の曲線部分(図1における矢印の点)への移行を決定するカテーテルコーナーを求める。
(c)このカテーテルコーナーに値足手マイクロカテーテルの中心曲線を調節する。
(d)次の直線部分の開始位置及び方向を求める。
4.マイクロカテーテル管の延長部分を調節する。延長部分は、“通常の”血管部分からネックを通して動脈瘤に続く部分である。
直線部分の後続の曲線部分(図1における矢印の点)への移行を決定するカテーテルコーナーを求めるための方法は、第3.1節において説明されている。カテーテルコーナーに対するマイクロカテーテル管の中心曲線の調節については、第3.2節において説明されている。次の直線部分の開始位置及び方向の演算の仕方については、第3.3節において説明されている。マイクロカテーテル管の開始部分及び延長部分の調節については、第3.4節において説明されている。第4節においては、我々の結果が与えられ、考慮のための一部の考慮が与えられている。
3.1 カテーテルコーナーの探索
直線部分の後続の曲線部分(図1における矢印の点)への移行を決定するカテーテルコーナーは、3つのテストプローブ(冒頭部分で既に説明したように、プローブは、球、その球の中心を通る平面及び複数の球パラメータの組合せである)を用いて求められる。直線部分の開始位置及び方向が演算された(第3.3節において説明している)後は、第1テストプローブの位置はこの直線部分の開始位置である。第1テストプローブの垂線は、この直線部分の正規化された方向である。第1テストプローブはまた、一次テスト光線を規定する。一次テスト光線は、第1テストプローブの垂線方向における第1テストプローブの位置から開始する。
第2テストプローブの最初の位置は、血管表面との一次テスト光線の最近接の交差部分により与えられる。交差部分が求まらない(図1の右側のピクチャの場合のような)場合、第2テストプローブは、第1テストプローブと第2テストプローブとの間の距離が表面バウンディングボックスの最大対角線に等しいように、一次テスト光線において位置付けられる。この場合、第2テストプローブは、血管壁の表面モデルの何れの三角形の頂点のように第1テストプローブから常に遠く離れている。第2テストプローブの垂線は、第1テストプローブの逆の垂線に等しい。
コリドール管(それ故、マイクロカテーテル管)が、一次テスト光線が血管壁と交差する前に、側枝を追従する場合(図2を参照されたい)、第2テストプローブは、コリドール管からは、それ故、マイクロカテーテル管の将来の中心線からは、遠過ぎる。この場合、第2テストプローブの最初の位置がカテーテルコーナーとして用いられる場合、マイクロカテーテル管の中心曲線はメアンダーを得る。
第2テストプローブは、第2テストプローブの平面までのこのコリドールプローブの距離が次式のように十分に小さいように、コリドールプローブ(指数kで表されている)が存在する場合、コリドール管に対して十分に近接していて、次のように表され、
t,2T(p−pt,2)≦rx1.1 (1)
ここで、
−nt,2は第2テストプローブの平面の垂線であり、
−pt,2は第2テストプローブの球の中心の位置であり、
−pはコリドールプローブkの球の中心の位置であり、
−rは、コリドールプローブkの楕円の主半径であり、係数1.1は局所的な表面の副規則性を可能にするように用いられる。
指数ibeginを有するコリドールプローブを用いて検査を開始することができる。このコリドールプローブは現在の直線部分の開始に対応している。
勿論、コリドールプローブibeginとコリドールプローブkとの間でコリドール管は側枝に進むことが可能である。それ故、検査されているコリドールプローブと第1及び第2テストプローブにより規定される線との間の距離は、次式のように、十分に小さい必要がある。
Figure 2008531108
ここで、
−l12は、第1テストプローブと第2テストプローブとの間の線であり、
−d(P,l12)はPとl12との間の距離である。
式(2)の違反は、コリドール管の中心曲線が側枝において遠くに曲がることを意味する。式(2)が違反する前の式(1)の正確さは、コリドール管の表面が第2テストプローブに隣接する血管壁に近いことを示している。
第2テストプローブの最初の位置がコリドール管から遠くに離れ過ぎている(コリドール管が側枝に追従することを表す)場合、第1プローブに近接する一次テスト光線における点を必要とする。図2は、一次テスト光線と側枝の挿入された上部表面(即ち、第1テストプローブの面から最も遠く離れた側枝の表面)との間の交差部分を必要とすることを示している。この交差部分と第1テストプローブとの間の距離は、側枝の上部表面の開始(即ち、一次テスト光線に最近接の上部表面の点)と第1テストプローブの面との間の距離に等しい。側枝におけるコリドールプローブと側枝の上部表面との間の距離は、コリドールプローブの楕円の主半径に略等しいために、側枝の上部表面の開始は、一次テスト構成の方向における血管表面との交差部分についてコリドールプローブを調べることにより求められることが可能である。
各々のコリドールプローブについて二次テスト光線を生成することができる。この二次テスト光線は、第1テストプローブの垂線方向において調べられるコリドールプローブの位置から始まる(図2を参照されたい)。二次テスト光線の血管表面との最近接の交差部分は、次式のように、調べられるコリドールプローブの位置p及び第1テストプローブの垂線nt,1の関数として第3テストプローブの位置pt,3(即ち、上部表面の点)を与える。
t,3=pt,3(p,nt,1 ) (3)
次式のように、第3テストプローブの対応する面までの距離が十分に小さい(このコリドールプローブが上部表面と交差することを意味する)第1コリドールプローブは第3テストプローブの最終位置を与える(図2を参照されたい)。
t,3 (p−pt,3(p,nt,1 ))≦rx1.1 (4)
ここで、kはこのコリドールプローブの指数であり、nt,3は第3プローブの垂線(一貫した距離基準についての第2テストプローブの垂線に等しい)である。
上記のように、第2テストプローブ(それ故、カテーテルコーナー)の第2テストプローブの最終位置と第1テストプローブの位置との間の距離は、次式のように、第3テストプローブの最終位置(即ち、側枝の上部表面の開始)から第1テストプローブにより規定される面までの距離に等しい(図2を参照されたい)。
Figure 2008531108
式(4)を満足するコリドールプローブが存在することが可能である。それ故、側枝の上部表面の開始のための探索は、次の2つの条件のうちの一が満たされる場合に停止される。
1.調べられているコリドールプローブ(指数iで表されている)から第2テストプローブの最初の面までの距離は、次式のように、十分に小さい。
t,2 (p−pt,2)≦rx1.1 (6)
この場合、側枝におけるコリドール管の中心曲線は第2テストプローブの面に非常に接近し、第2テストプローブの面の他の側においてさえ、継続されることが可能である。
2.調べられているコリドールプローブから第1及び第2テストプローブにより規定される線までの距離は、次式のように、大き過ぎる。
d(p,l12)≦2xmax(r,j∈[1;Nprobes]) (7)
この場合、側枝におけるコリドール管の中心曲線は、第1プローブと第2プローブとの間の線から大きく外れて進む。
側枝の上部表面の開始のための探索が停止される場合。第2テストプローブの最終位置と第1テストプローブの位置との間の距離は、次式のように、第3テストプローブの位置から第1テストプローブにより規定される面までの距離を最小にするように設定される。
Figure 2008531108
ここで、iは調べられているコリドールプローブの指数である。
3.2 カテーテルコーナーに対する調節
カテーテルコーナーが求められた後、マイクロカテーテル管の残りの部分は、マイクロカテーテルプローブに有効に変化する移動ベクトルを適用することにより、このコーナーに対して調節される必要がある。この残りの部分は、現在の直線部分の開始に対応するマイクロカテーテルプローブ(即ち、式(1)及び式(2)において用いられるような指数ibeginを有するマイクロカテーテルプローブ)から始まる
マイクロカテーテル管の中心曲線の先行する部分(即ち、ibeginより小さい指数を有するマイクロカテーテルプローブの球の中心)は既に調節されていることに留意されたい。
残りの部分は2つの部分に副分割される。第1部分は、現在の直線部分からカテーテルコーナーまで続くマイクロカテーテル管の一部である。第2部分は、カテーテルコーナーからマイクロカテーテルの終端まで続くマイクロカテーテル管の部分である。(第1及び第2部分への副分割については、下で更に詳細に説明される)。
マイクロカテーテル管は、次の目的が達成されるように調節される。
1.マイクロカテーテルプローブの球の中心と血管壁との間の距離がマイクロカテーテル管の半径より略大きい。実際には、マイクロカテーテル管は、実際には、コリドール管の内側(即ち、血管を通るコリドールの内側)にある必要がある。
2.視覚的不連続性がない。視覚的不連続性を回避するように、移動ベクトルの方向は一定に保たれる。
3.マイクロカテーテルプローブの球の中心は、対応するコリドールプローブの球の中心にできるだけ近い。それ故、移動ベクトルの最大の大きさはできるだけ小さい必要がある。
4.第1部分の中心曲線は、現在の直線部分の開始とカテーテルコーナーとの間の直線部分にできるだけ近い。
5.第2部分の開始部分は、マイクロカテーテルを曲げる血管壁にできるだけ近い。
6.局所的メアンダーは存在しない。局所的メアンダーを回避するように、第2部分の移動ベクトルの大きさは、マイクロカテーテルプローブの指数の関数として、単調に減少する。
7.第1部分と第2部分との間の接続はスムーズである。
第1部分及び第2部分におけるマイクロカテーテル管の残りの部分の副分割の仕方について、先ず、説明する。次に、第2部分の制約された移動について説明する。その後、第1部分と第2部分との間のスムーズな移行を有する第1部分の調節について説明する。
第1及び第2部分における副分割
マイクロカテーテルの残りの部分の第1及び第2部分への副分割は目的3に基づいている。結局、移動ベクトルの最大の大きさは、カテーテルコーナーに対して最近接のマイクロカテーテルプローブ(以下、“コーナープローブ”という)がカテーテルコーナーの方への移動について選択される場合に、最小である。ヘアピン状に曲がっている場合、プローブが、カテーテルコーナーの方に直線に沿って組織により分離され、血管において遠く離れているように選択されることが可能である。それ故、次式、即ち
Figure 2008531108
及び
Figure 2008531108
及び
Figure 2008531108
であるように、指数iを有するコーナープローブが選択される。
最後の式は、テストするために用いられる最後のプローブの指数iendを決定する。
指数ibeginを有する第1プローブがこの式を既に違反している場合、第1プローブはコーナープローブとして用いられる。
カテーテルコーナーにできるだけ近い中心曲線を得るように、最初の移動ベクトルは、次式のように、カテーテルコーナーとコーナープローブとの間のベクトルに等しい必要がある。
initial=picorner−pt,2 (12)
第2部分の制約された移動
第2部分のマイクロカテーテルプローブは、マイクロカテーテル管が実際にコリドール管の内側に保たれるように、移動ベクトルの方向に移動される(3.2節の目的1)。移動されるべきマイクロカテーテルプローブの古い位置poldは、それが最初の位置、即ち、対応するコリドールプローブの球の中心に等しいためにか又は、前に制約された移動の結果であるために、既に許容されていることに留意されたい。
典型的な構成については、図3に示されている(水平線部分については下で説明される)。移動されるべき現在のマイクロカテーテルプローブの球の中心の古い位置poldは、コリドールプローブil,oldの面とil+1,oldの面との間に位置付けられる。現在の移動ベクトルvcurrent(即ち、最初の移動ベクトルvinitialか又は前のマイクロカテーテルプローブについて許容される新しい位置を与えられた移動ベクトル)を有する一時的な新しい位置、即ち
new=pold+vcurrent (13)
がコリドールプローブil,newの面とil+1,newの面との間に位置付けられる。
移動ベクトルの最大の許容された大きさ
コリドール管の表面を有する線部分pold→pnewの正確な位置の演算は非常に複雑であり、間違いを起こし易く、時間を要するタスクである。それ故、最小min(il,old,il,new)と最大max(il+1,old,il+1,new)との間のコリドールプローブの面におけるこの線部分の投影から移動ベクトルの最大の許容される大きさが近似される。
典型的な投影の結果については、図4に示されている。pp,oldは、検査されるコリドールプローブの面におけるpoldの投影であり、pp,newはpnewの投影である。円(下で説明される)の中心位置cはコリドールプローブの球の中心位置と同じである。その円の半径は、楕円の短半径rとマイクロカテーテルの半径rとの間の差分に等しい。
ところで、投影pp,oldがその円の外側に位置している(古い位置poldは、実際には、コリドール管の内側にあるために、投影がその円に非常に近くに位置している必要がある)場合、投影は、次式のように、円の内側に位置するようになるまで、円の中心の方に移動される。
Figure 2008531108
検査されるコリドールプローブに垂直な面及び円により規定される局所的な円筒により、マイクロカテーテル管の中心曲線が位置付けられる必要があるコリドール管の局所的表面が近似される(図3における水平線部分はそれらの局所的円筒の上部の境界を示している)。
この場合、poldと、この局所的円筒との交差部分と、の間の線部分pold→pnewの一部は、pp,oldとその円を有する交差部分との間の線部分pp,old→pp,newの一部fに等しい。それ故、この投影からの移動ベクトルの最大の許容される大きさは、次式のようになる。
Figure 2008531108
有する線部分pp,old→pp,newの円との交差部分は、その円の半径に対してそのベクトルの長さを等しくすることにより、次式のように、与えられる。
Figure 2008531108
図4に示す投影の結果については、式(16)が、1.0より小さいfについての1つの正の解を与える。投影pp,newがその円の内側に位置する場合、式(16)は、1.0より大きいfについての1つの正の解を与える。
古い位置及び新しい位置又はそれらの投影が一致する場合、移動ベクトルの最大の許容される大きさは、投影pp,newの位置に依存する。この投影が円の内側に位置している(線部分pold→pnewが、その場合、局所的円筒の内側に位置している)場合、現在の移動ベクトルvcurrentは、検査されるコリドールプローブについて許容可能である。これを示すように、移動ベクトルの最大の許容される大きさは、現在の移動ベクトルの大きさより大きい値に設定される。投影pp,newが円の外側に位置している場合、移動ベクトルの最大の許容される大きさは0に設定される。
安全性の理由のために、移動ベクトルの最終の最大の許容される大きさは、含まれるコリドールプローブの面における投影から演算される値の最小である。
2つの連続するコリドールプローブの面の間の距離は、2つのボクセル間の距離に略等しいことに留意されたい。それ故、局所的円筒による近似のための誤差は、ボクセル間の距離と同じ大きさである。
マイクロカテーテルプローブの位置の更新
移動ベクトル
Figure 2008531108
の最大の許容される大きさが、現在の移動ベクトルの大きさ
Figure 2008531108
より大きい又はそれに等しい場合、一時的な新しい位置pnewは最終の新しい位置になる。移動ベクトルの最大の許容される大きさが現在の移動ベクトルの大きさより小さい場合、現在の移動ベクトルは、次式のように調節され、
Figure 2008531108
最終の新しい位置は、次式のようになる。
new=pold+vcurrent (18)
第1部分の調節
コーナープローブ(第3.2節を参照されたい)がカテーテルコーナー(第3.1節を参照されたい)の方にできるだけ近接するように移動された(第3.2節を参照されたい)後に、第1部分の有効なマイクロカテーテルプローブは、次式のように、第1部分の最初のプローブと最後のプローブとの間で一時的に変わる大きさを有する同じ方向に移動され、
Figure 2008531108
ここで、vcornerはカテーテルコーナーの方にコーナープローブを移動させるように用いられる移動ベクトルである。
コーナープローブはできるだけカテーテルコーナーに近接しているために、そして第1部分の第1プローブとカテーテルコーナーとの間の直線部分は血管の内側に位置しているために、式(19)の移動は、第3.2節の目的1、2及び4と適合する。
第1及び第2部分のマイクロカテーテルプローブがカテーテルコーナーに対して調節された後、マイクロカテーテル管の全体の中心曲線は、第3.2節の目的7と適合するように制約された緩和アルゴリズム(参照文献7を参照されたい)によりスムーズにされる。緩和の間に用いられる制約は、反復kにおいて与えられるマイクロカテーテルプローブ
iの新しい位置pi,k+1は、次式、
Figure 2008531108
及び
Figure 2008531108
のように、コリドール管の内側に位置される必要があり、ここで、コリドールプローブl及びl+1は、マイクロカテーテルプローブの一時的な新しい位置pi,k+1がそれらのコリドールプローブ(p及びpl+1はがそれら2つのコリドールプローブの位置である)の面の間に位置している。任意に選択された係数0.9は、マイクロカテーテル管の小さい一部がコリドール管の外側に位置付けられることが可能であるために、マイクロカテーテル管の中心曲線の良好なスムージングを可能にする。
3.3 新しい直線部分
第3.2節に記載されている制約された移動は、移動ベクトル
Figure 2008531108
の最大の許容された大きさを第2部分の各々のマイクロカテーテルプローブについて演算する。現在の移動ベクトル
Figure 2008531108
の大きさが、移動ベクトルの最大の許容された大きさより大きい又はそれに等しい限り、マイクロカテーテル管は、コリドール管の表面により表されるような血管壁により曲げられる。現在の移動ベクトルの大きさが移動ベクトルの最大の許容された大きさより小さい第2部分の第1マイクロカテーテルプローブは、第1の制約されていないマイクロカテーテルプローブである。
次の直線部分の開始位置として第1の制約されていないマイクロカテーテルプローブに先行するマイクロカテーテルプローブの位置が用いられる。このマイクロカテーテルプローブ及びそれに先行するマイクロカテーテルプローブは、一般に、血管壁により曲げられるため、次の直線部分の方向としてそれらのマイクロカテーテルプローブの位置の間の正規化ベクトルが用いられる。実際には、この正規化ベクトルは、マイクロカテーテルが血管壁を出る前に、血管壁により誘導される最後の修正操舵を表す。
第2部分の全てのマイクロカテーテルプローブが制約を受ける場合、マイクロカテーテル管の演算は終了する。
3.4 マイクロカテーテル管の先端
第2節において説明しているように、コリドール管は、血管(即ち、ネック中心までの“標準的な”血管における部分)と延長管(即ち、ネック中心から動脈瘤までの部分)の連結である。前の節において説明されたアルゴリズムは、血管に対応するマイクロカテーテル管の一部のみに適用される。実際には、延長管に対応するマイクロカテーテル管の一部にこのアルゴリズムを適用することにより、動脈瘤の内側における選択された端部位置から動脈瘤の境界の方に離れるようにこの部分が移動される。実際には、延長部分は、マイクロカテーテル管整形アルゴリズムが適用される前に、最初のマイクロカテーテル管(即ち、マイクロカテーテルの半径により置き換えられる半径を有するコリドール管のコピー)から取り除かれる。
マイクロカテーテル管の血管部分の最後のプローブの位置(及び方向)はマイクロカテーテル管整形アルゴリズムにより変化されることが可能であるために、古い延長部分を有する新しいマイクロカテーテル管の容易な連結は、マイクロカテーテル管の中心曲線及び表面における視覚的不連続性をもたらす可能性がある。それ故、マイクロカテーテル管についての延長部分は、マイクロカテーテル管の血管部分の最後のプローブを用いて生成される。
例えば、三次元回転血管造影(参照文献6を参照されたい)により生成されるボリュームは、明確化のために、全体的な血管構造の副集合のみを主に有する。それ故、コリドール管は、一般に、イントロデューサシースから遠く離れた血管部分のどこかから開始する。それ故、最初のコリドールプローブの位置及び方向は、血管構造のその位置におけるマイクロカテーテルの実際の位置及び方向とは異なる可能性がある。実際には、マイクロカテーテルの表されていない先行する部分に適用される血管壁の曲がりは、その血管の中心軸に近いことに代えて、血管壁に近いマイクロカテーテルの表されている部分の開始位置をもたらす。
マイクロカテーテルの表されている部分の開始位置及び方向を改善するように、任意の初期移動ベクトルは、第1カテーテルコーナーが探索される前に、マイクロカテーテル管に適用されることが可能である(第3.2節において記載されている制約された移動アルゴリズムを用いて)。本発明で示すプログラム(上記のように、適切なアルゴリズムは任意の初期移動ベクトルを可能にする)は、次のような5つの所定の初期移動ベクトルを有し、
initial=rx(uxuaxis+vxvaxis) (22)
(u,v)∈[(0,0),(1,0),(0,1),(−1,0),(0,−1)] (23)
ここで、rは第1コリドールプローブの半径よりかなり大きい半径であり、uaxis及びvaxisは第1コリドールプローブの面における局所的座標系である。
4.結果及び結論
三次元Integris system(参照文献6を参照されたい)により取得される28個の医療用ボリュームデータ集合に対してコリドール管からのマイクロカテーテル管の演算のための方法が適用される。ボリュームの大きさは128x128x128である。18個の大動脈瘤は分岐に位置付けられ、10個は単一の血管部分に位置付けられている。
コリドール管の演算のための平均経過時間はSGI Octane(3000MHz MIPS R12000+MIPS R12010 FPU)において2.5秒である。マイクロカテーテル管の演算についての経過時間は、平均して、対応するコリドール管についての演算時間の20%である。
図5.4は中央曲線を示し、図5.6は、図5.3に示されている中央曲線及び図5.5に示されている表面を有するコリドール管からもたらされるマイクロカテーテル管の表面を示している。
本発明の方法を評価するように、次式のようにコリドール管の表面とマイクロカテーテル管の表面との間の相対的距離rdを各々のマイクロカテーテルプローブについて評価し、
Figure 2008531108
ここで、コリドールプローブl及びl+1は、マイクロカテーテルプローブの位置pがそれらのコリドールプローブの面間に位置するように選択されている。p及びpl+1はそれらの位置であり、r及びrl+1はそれら2つのコリドールプローブの楕円の短半径である。pi,l及びpi,l+1は、コリドールプローブの面に対するマイクロカテーテルプローブの位置pの投影である。
この相対的距離は、マイクロカテーテル管が部分的にコリドール管の外側にある場合に不であり、それらの表面が一致する場合に0であり、そしてマイクロカテーテル管が局所的に完全にコリドール管の内側にある場合に正である。この相対的な距離は、マイクロカテーテル管及びコリドールプローブの中心位置が一致する(即ち、マイクロカテーテル管の初期状態の)場合、1.0に等しい。
場合毎の及び4つの所定の初期移動ベクトルについての平均相対的距離が演算される(第3.4節における式(22)及び(23)を参照されたい)。それらの平均相対的距離の統計データが表1に与えられている。
Figure 2008531108
次のような結論が、上記結果、ピクチャ及びテスト中に収集された経験から導き出される。
1. コリドール管からマイクロカテーテルを演算するための方法は、視覚的に許容可能な結果を与える。医療的妥当性は、複数の診療所で開始されている(診療的評価については、後続の論文において報告されることになっている)。
2. その相対的距離(表1を参照されたい)は、本発明の方法の高効率性は、選択された初期移動に統計的に依存しない。
3. 血管が折れ曲がっているときに、反対側にコリドール管を交差させるマイクロカテーテル管の一部からみて、その相対的距離は、本発明の方法の効率が非常に良好であることを示している。しかしながら、“golden standard”は(未だに)利用可能でないため、視覚的検査による粗い検証のみが可能である。
4. マイクロカテーテル管(図5.6を参照されたい)は、動脈瘤内に容易に移動させるための実際のマイクロカテーテルの前モデル化及び選択のための開始点として用いられることができる。
コリドール管からのマイクロカテーテル管の演算のための本発明の方法は、オペレータ管の変動及びオペレータ内の変動を最小化することができるため、選択された及び前モデル化された実際のマイクロカテーテルは良好な品質を有する及び/又は患者の治療のために迅速に利用可能であることが期待できる。
最後に、本出願において、用語“を有する”は他の要素又はステップを排除するものではなく、単数表現は複数の存在を排除する者ではなく、そして単独の処理器又は他のユニットは複数の手段の機能を実行することが可能であることに注意する必要がある。本発明は、各々の及び全ての新規な特徴的な特性並びに各々の及び全ての特徴的な特性の組合せを網羅している。
血管系の部分(左側)及び血管系の分岐部分(右側)のそれぞれにおける中心線CC及びMCと共にコリドール管CT及びマイクロカテーテル管MTを模式的に示す図である。 血管系の側枝におけるカテーテルコーナーの決定であって、直線部分は前記カテーテルコーナーにおいて曲線部分に移行する、決定を示す図である。 曲線部分における位置poldから位置pnewへのマイクロカテーテルプローブの移動であって、ドットは、この領域における対応するコリドール管を表すプローブの中心を表す、移動を示す図である。 最大の移動ベクトルの長さの演算のためにコリドール管プローブの面に図3の移動ベクトルの投影を示す図である。 動脈瘤を有する血管系の異なる三次元表現であって、即ち、1.左上:濃淡値ボリューム2.右上:レベル化ボリューム、動脈瘤は黒色で印付けされている3.左中央:コリドール管の中心曲線4.右中央:マイクロカテーテル管の中心曲線5.左下:コリドール管の表面6.右下:マイクロカテーテル管の表面

Claims (12)

  1. モデル化血管系における開始位置と目的位置との間のカテーテルの経路の予測のための方法であって、前記経路は、関連経路中心線に沿って進む経路管により表される、方法であり:
    a)前記開始位置から前記目的位置まで前記血管系を通る経路を決定し、初期経路中心線を前記経路と一致させる段階;及び
    b)関連する前記経路管が前記血管系内にあるように、前記初期経路中心線を調節する段階;
    を有することを特徴とする方法。
  2. 請求項1に記載の方法であって、前記経路管は、カテーテルが前記開始位置から前記血管系を通って前記目的位置まで進むことが可能であるコリドールを表すコリドール管である、ことを特徴とする方法。
  3. 請求項1に記載の方法であって、前記経路管は、前記開始位置から前記血管系を通って前記目的位置まで進むマイクロカテーテルの予測形状を表すマイクロカテーテル管である、ことを特徴とする方法。
  4. 請求項2又は3に記載の方法であって、前記コリドール管が先ず、決定され、次に、マイクロカテーテル管が、コリドール管内にあるように決定される、ことを特徴とする方法。
  5. 請求項3に記載の方法であって、前記マイクロカテーテル中心線は、関連管の部分が前記血管系の内側にある直線部分と、前記関連管の部分が血管壁に接する及び/又は前記血管系の側枝に入る曲線部分とが交互に繋がっている連続を有する、ことを特徴とする方法。
  6. 請求項5に記載の方法であって、前記連続は交互に決定され、好適には、前記開始位置における直線部分から始まる、ことを特徴とする方法。
  7. 請求項6に記載の方法であって、各々の反復段階は:
    aa)(i)周りの血管壁との現在の直線部分の交差部分として、又は(ii)マイクロカテーテルが後続する前記側枝の最も遠い血管壁と前記部分の開始から同じ距離にある前記現在の直線部分における点として、カテーテルコーナーを決定する段階;
    bb)前記カテーテルコーナーの方に初期移動ベクトルだけ前記カテーテルコーナーの方に近い前記現在の直線部分の点を移動させる段階;及び
    cc)前記移動された点において前記現在の直線部分から後続の曲線部分への移行を導入する段階;
    を有する、ことを特徴とする方法。
  8. 請求項7に記載の方法であって、前記後続の曲線部分の前記マイクロカテーテル中心線は、前記関連管が前記血管系の壁に接するように単調に減少された移動長さを有する前記初期移動ベクトルの方向に部分毎に移動され、前記後続の直線部分は、前記血管壁に接触していないところで開始する、ことを特徴とする方法。
  9. 請求項1に記載の方法であって、前記血管系を通る管は、前記管の前記中心線上にある中心を有し、前記中心を有し且つ前記管の前記中心線に対して垂直に進む面を有する球を有するプローブによりモデル化される、ことを特徴とする方法。
  10. カテーテル、特に、マイクロカテーテルの製造のための方法であって:
    a)請求項1乃至9の何れの一項に記載の方法により前記カテーテルの経路を予測する段階;及び
    b)前記予測された経路にしたがって前記カテーテルを準備する段階;
    を有する方法。
  11. 請求項1乃至9の何れの一項に記載の方法を実行するように適合されたデータ処理ユニット。
  12. カテーテルの経路の予測のためのコンピュータプログラムを記憶する記録担体であって、前記プログラムは請求項1乃至9の何れの一項に記載の方法を実行するように適合されている、記録担体。
JP2007556700A 2005-02-23 2006-02-20 カテーテルの経路の予測のための方法 Pending JP2008531108A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP05101378 2005-02-23
PCT/IB2006/050547 WO2006090324A2 (en) 2005-02-23 2006-02-20 Method for the prediction of the course of a catheter

Publications (1)

Publication Number Publication Date
JP2008531108A true JP2008531108A (ja) 2008-08-14

Family

ID=36702642

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007556700A Pending JP2008531108A (ja) 2005-02-23 2006-02-20 カテーテルの経路の予測のための方法

Country Status (5)

Country Link
US (1) US20080160489A1 (ja)
EP (1) EP1856642A2 (ja)
JP (1) JP2008531108A (ja)
CN (1) CN101128829B (ja)
WO (1) WO2006090324A2 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015102551A (ja) * 2013-11-21 2015-06-04 バイオセンス・ウエブスター・(イスラエル)・リミテッドBiosense Webster (Israel), Ltd. インピーダンス測定を用いた挿入点から心臓までのカテーテルの追跡
JP5890055B1 (ja) * 2015-07-09 2016-03-22 株式会社アルム 血管画像処理装置、血管画像処理プログラム、および血管画像処理方法
JP2019030678A (ja) * 2013-03-15 2019-02-28 ライトラボ・イメージング・インコーポレーテッド 較正及び画像処理装置、方法、並びにシステム
JP2022512528A (ja) * 2018-10-16 2022-02-07 ユニオン・ストロング(ベイジン)・テクノロジー・カンパニー・リミテッド 血管内でのカテーテルの移動経路をシミュレートするための方法、装置、およびデバイス
JP7462051B2 (ja) 2020-01-09 2024-04-04 キヤノン ユーエスエイ,インコーポレイテッド 湾曲器具によるプランニング及び可視化の改良、並びにその湾曲器具

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2138095A1 (de) * 2008-06-25 2009-12-30 BrainLAB AG Verfahren zur Ermittlung der Position eines medizinischen Instruments in einem Körper
US9286671B2 (en) * 2011-01-20 2016-03-15 Medtronic Bakken Research Center B.V. Method for determining at least one applicable path of movement for an object in tissue
WO2018002250A1 (en) * 2016-06-30 2018-01-04 Koninklijke Philips N.V. Catheter type selection
JP2019524200A (ja) 2016-06-30 2019-09-05 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. カテーテルタイプの選択
CN109452972B (zh) * 2018-10-16 2022-03-01 复旦大学附属华山医院 一种导管塑形器形状模拟方法、装置及设备
CN109452971B (zh) * 2018-10-16 2021-09-17 复旦大学附属华山医院 一种导管在血管内的行进路线模拟方法、装置及设备
CN109512510B (zh) * 2018-10-16 2021-09-28 复旦大学附属华山医院 一种导管在血管内的行进路线模拟方法、装置及设备
NL2021849B1 (en) 2018-10-22 2020-05-13 Mat Nv System and method for catheter based intervention
CN110179548B (zh) * 2019-06-27 2021-01-08 浙江大学医学院附属第一医院 一种预测前端修剪picc导管置管长度的方法
CN114366296B (zh) * 2021-12-31 2023-05-30 杭州脉流科技有限公司 改进的微导管路径生成方法、塑形方法、设备和存储介质

Family Cites Families (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5694530A (en) * 1994-01-18 1997-12-02 Hitachi Medical Corporation Method of constructing three-dimensional image according to central projection method and apparatus for same
US5920319A (en) * 1994-10-27 1999-07-06 Wake Forest University Automatic analysis in virtual endoscopy
SG64340A1 (en) * 1996-02-27 1999-04-27 Inst Of Systems Science Nation Curved surgical instruments and methods of mapping a curved path for stereotactic surgery
US6167296A (en) * 1996-06-28 2000-12-26 The Board Of Trustees Of The Leland Stanford Junior University Method for volumetric image navigation
US5971767A (en) * 1996-09-16 1999-10-26 The Research Foundation Of State University Of New York System and method for performing a three-dimensional virtual examination
US6343936B1 (en) * 1996-09-16 2002-02-05 The Research Foundation Of State University Of New York System and method for performing a three-dimensional virtual examination, navigation and visualization
US6331116B1 (en) * 1996-09-16 2001-12-18 The Research Foundation Of State University Of New York System and method for performing a three-dimensional virtual segmentation and examination
US6016439A (en) * 1996-10-15 2000-01-18 Biosense, Inc. Method and apparatus for synthetic viewpoint imaging
US5891030A (en) * 1997-01-24 1999-04-06 Mayo Foundation For Medical Education And Research System for two dimensional and three dimensional imaging of tubular structures in the human body
US6369812B1 (en) * 1997-11-26 2002-04-09 Philips Medical Systems, (Cleveland), Inc. Inter-active viewing system for generating virtual endoscopy studies of medical diagnostic data with a continuous sequence of spherical panoramic views and viewing the studies over networks
US6928314B1 (en) * 1998-01-23 2005-08-09 Mayo Foundation For Medical Education And Research System for two-dimensional and three-dimensional imaging of tubular structures in the human body
DE69805209T2 (de) * 1998-02-23 2002-11-28 Algotec Systems Ltd System und methode zur automatischen wegplanung
US6606091B2 (en) * 2000-02-07 2003-08-12 Siemens Corporate Research, Inc. System for interactive 3D object extraction from slice-based medical images
EP1264281A4 (en) * 2000-02-25 2007-07-11 Univ New York State Res Found ARRANGEMENT AND METHOD FOR PROCESSING AND PLAYING A VOLUME
US6511418B2 (en) * 2000-03-30 2003-01-28 The Board Of Trustees Of The Leland Stanford Junior University Apparatus and method for calibrating and endoscope
WO2002086797A1 (en) * 2001-03-06 2002-10-31 The John Hopkins University School Of Medicine Simulation method for designing customized medical devices
US7190365B2 (en) * 2001-09-06 2007-03-13 Schlumberger Technology Corporation Method for navigating in a multi-scale three-dimensional scene
US20030152897A1 (en) * 2001-12-20 2003-08-14 Bernhard Geiger Automatic navigation for virtual endoscopy
SG165160A1 (en) * 2002-05-06 2010-10-28 Univ Johns Hopkins Simulation system for medical procedures
US7187790B2 (en) * 2002-12-18 2007-03-06 Ge Medical Systems Global Technology Company, Llc Data processing and feedback method and system
US7081088B2 (en) * 2003-01-30 2006-07-25 Siemens Corporate Research, Inc. Method and apparatus for automatic local path planning for virtual colonoscopy
ES2670344T3 (es) * 2003-09-05 2018-05-30 Infineum International Limited Composiciones de aditivo para combustible diésel estabilizado
US20050183325A1 (en) * 2004-02-24 2005-08-25 Sutkowski Andrew C. Conductivity improving additive for fuel oil compositions
US8457373B2 (en) * 2009-03-16 2013-06-04 Siemens Aktiengesellschaft System and method for robust 2D-3D image registration

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2019030678A (ja) * 2013-03-15 2019-02-28 ライトラボ・イメージング・インコーポレーテッド 較正及び画像処理装置、方法、並びにシステム
US11435233B2 (en) 2013-03-15 2022-09-06 Lightlab Imaging, Inc. Calibration and image processing devices, methods, and systems
JP2020199317A (ja) * 2013-03-15 2020-12-17 ライトラボ・イメージング・インコーポレーテッド 較正及び画像処理装置、方法、並びにシステム
JP2015102551A (ja) * 2013-11-21 2015-06-04 バイオセンス・ウエブスター・(イスラエル)・リミテッドBiosense Webster (Israel), Ltd. インピーダンス測定を用いた挿入点から心臓までのカテーテルの追跡
JP2017018234A (ja) * 2015-07-09 2017-01-26 株式会社アルム 血管画像処理装置、血管画像処理プログラム、および血管画像処理方法
CN107847207B (zh) * 2015-07-09 2018-12-18 Pentas股份有限公司 血管图像处理装置、存储介质和血管图像处理方法
CN107847207A (zh) * 2015-07-09 2018-03-27 Pentas股份有限公司 血管图像处理装置、血管图像处理程序和血管图像处理方法
US10373326B2 (en) 2015-07-09 2019-08-06 Pentas Inc. Blood vessel image processing apparatus, blood vessel image processing program, and blood vessel image processing method
WO2017006732A1 (ja) * 2015-07-09 2017-01-12 株式会社Pentas 血管画像処理装置、血管画像処理プログラム、および血管画像処理方法
JP5890055B1 (ja) * 2015-07-09 2016-03-22 株式会社アルム 血管画像処理装置、血管画像処理プログラム、および血管画像処理方法
JP2022512528A (ja) * 2018-10-16 2022-02-07 ユニオン・ストロング(ベイジン)・テクノロジー・カンパニー・リミテッド 血管内でのカテーテルの移動経路をシミュレートするための方法、装置、およびデバイス
JP7425039B2 (ja) 2018-10-16 2024-01-30 強聯智創(北京)科技有限公司 血管内でのカテーテルの移動経路をシミュレートするための方法、装置、およびデバイス
JP7462051B2 (ja) 2020-01-09 2024-04-04 キヤノン ユーエスエイ,インコーポレイテッド 湾曲器具によるプランニング及び可視化の改良、並びにその湾曲器具

Also Published As

Publication number Publication date
WO2006090324A2 (en) 2006-08-31
WO2006090324A3 (en) 2007-03-15
EP1856642A2 (en) 2007-11-21
US20080160489A1 (en) 2008-07-03
CN101128829A (zh) 2008-02-20
CN101128829B (zh) 2010-05-19

Similar Documents

Publication Publication Date Title
JP2008531108A (ja) カテーテルの経路の予測のための方法
JP7267743B2 (ja) 脈管内においてカテーテルのような導管ラインをルーティングするためのシステムおよび方法
US10580139B2 (en) Automated centerline extraction method for determining trajectory
US7371067B2 (en) Simulation method for designing customized medical devices
US7167180B1 (en) Automatic path planning system and method
US20190251755A1 (en) Interactive voxel manipulation in volumetric medical imaging for virtual motion, deformable tissue, and virtual radiological dissection
US9830427B2 (en) Method for intracranial aneurysm analysis and endovascular intervention planning
US8867801B2 (en) Method for determining properties of a vessel in a medical image
JP6797200B2 (ja) 血管構造内で血管内器具を誘導する助けとなるためのシステム及び同システムの作動方法
US7565000B2 (en) Method and apparatus for semi-automatic segmentation technique for low-contrast tubular shaped objects
US9443317B2 (en) Image display of a centerline of tubular structure
EP3025303B1 (en) Multi-modal segmentation of image data
US10535189B2 (en) Image display of a centerline of tubular structure
US11295537B1 (en) Method and apparatus for inserting virtual contrast into CT and MRI examinations
US20230210600A1 (en) Microcatheter path generation method, shaping method of mandrel, computer equipment, readable storage medium and program product
Bruijns et al. Fully automatic computation of the shape of a micro-catheter
Subramanian et al. A prototype virtual reality system for preoperative planning of neuro-endovascular interventions
Nguyen et al. Vascular flow rendering for interactive simulation of contrast and drugs injection

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20090218

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110628

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110928

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20120327