(57)【要約】複数の高分解能画像スライスに対応するデータを備えている高分解能容量画像データ（７４）を取得（７０）し、かつ少なくとも一つの血管識別画像スライス（７６）に対応するデータを取得（７２）することを有する磁気共鳴血管造影方法であって、この取得されたデータは、動脈及び静脈の一つに関して選択的に強調されたコントラストを有している。高分解能容量画像表示（８０）は、取得された高分解能容量画像データ（７４）から再構築される。少なくとも一つの血管識別スライス画像表示（８２）は、少なくとも一つの血管識別画像スライス（７６）に対応する取得されたデータから再構築される。動脈開始点（８６）及び静脈開始点（８８）の少なくとも一つは、血管識別スライス画像表示（８２）に基づいて特定される。好ましくは、当該方法は、前記の少なくとも一つの開始点（１２６）を囲む開始表面を特定（１２４）し、かつこの表面が血管系の境界に本質的に一致するまで、開始表面を反復的に成長（１２８）させることをさらに含んでいる。
2004535869.tif 000002
【特許請求の範囲】
【請求項１】
高分解能画像容量に対応するデータを有する高分解能容量画像データを取得し；
少なくとも一つの血管識別画像スライスに対応するデータを取得し、該取得データは、動脈及び静脈の一つに関して選択的に強調されたコントラストを有し；
前記の取得された高分解能容量画像データから高分解能容量画像表示を再構築し；
少なくとも一つの血管識別画像スライスに対応する前記取得データから、少なくとも一つの血管識別スライス画像表示を再構築し；かつ
前記血管識別スライス画像表示に基づいて、前記高分解能画像において、少なくとも一つの、動脈開始点及び静脈開始点を特定する；
ことを有する磁気共鳴血管造影方法。
【請求項２】
前記の少なくとも一つの血管識別画像スライスに対応するデータを取得するステップは、前記の高分解能容量画像データの取得の間に分散されることを特徴とする請求項１に記載の磁気共鳴血管造影方法。
【請求項３】
前記の少なくとも一つの血管識別画像スライスに対応するデータを取得するステップは：
動脈及び静脈の一つに関して選択的に強調されたコントラストを供するように選択的前飽和を用いてデータを取得する；か
動脈及び静脈の一つに関して選択的に強調されたコントラストを供するように、選択的前飽和を用いてタイム・オブ・フライトＭＲＡデータを取得する；
かのいずれか一つを有することを特徴とする請求項１又は２に記載の磁気共鳴血管造影方法。
【請求項４】
前記の少なくとも一つの血管識別画像スライスに対応するデータを取得するステップは：
単位スライス当たり約１秒以下に対応するデータ取得率にて前記の少なくとも一つの血管識別画像データ、を取得することを有する；
ことを特徴とする請求項１乃至３のいずれか一項に記載の磁気共鳴血管造影方法。
【請求項５】
前記の少なくとも一つの血管識別画像スライスに対応するデータを取得するステップは：
動脈及び静脈の一つのルートに対応している領域のデータを取得する；
ことを有することを特徴とする請求項１乃至４のいずれか一項に記載の磁気共鳴血管造影方法。
【請求項６】
前記の少なくとも一つの血管識別画像スライスに対応するデータを取得するステップが含むのは：
少なくとも一つの血管識別スライスに対応するデータを取得することであり、該取得された血管識別画像スライスは、複数の前記高分解能画像のスライスの一つに空間的に対応している；
ことを特徴とする請求項１乃至５のいずれか一項に記載の磁気共鳴血管造影方法。
【請求項７】
前記の識別画像スライスに関するデータの取得及び前記の対応する高分解能画像スライスに関するデータの取得は分散されていることを特徴とする請求項６に記載の磁気共鳴血管造影方法。
【請求項８】
前記の血管識別スライス画像表示に基づいた動脈開始点及び静脈開始点の少なくとも一つを特定するステップは：
前記血管識別スライス画像表示に基づいて少なくとも一つの動脈開始点を特定し；かつ
前記血管識別スライス画像表示及び前記の少なくとも一つの動脈開始点に関する前記血管識別スライス画像表示に空間的に対応する前記高分解能画像スライスの比較に基づいて少なくとも一つの静脈開始点を特定する；
ことを有することを特徴とする請求項６又は７に記載の磁気共鳴血管造影方法。
【請求項９】
前記の血管識別スライス画像表示に基づいた動脈開始点及び静脈開始点の少なくとも一つを特定するステップは：
前記血管識別スライス画像表示に基づいて、少なくとも一つの静脈開始点を特定し；かつ
前記血管識別スライス画像表示及び前記の少なくとも一つの静脈開始点に関する前記血管識別スライス画像表示に空間的に対応する前記高分解能画像スライスを比較することに基づいて少なくとも一つの動脈開始点を特定する；
ことを有することを特徴とする請求項６又は７に記載の磁気共鳴血管造影方法。
【請求項１０】
前記の少なくとも一つの開始点における血管下位組織の開始を三次元で追跡する；
ことをさらに有する請求項１乃至９のいずれか一項に記載の磁気共鳴血管造影方法。
【請求項１１】
前記の追跡ステップは：
前記開始点を囲む開始表面を規定し；かつ
前記表面が本質的に前記血管系の境界に一致するまで前記開始表面を反復的に成長させる；
ことを有する請求項１０に記載の磁気共鳴血管造影方法。
【請求項１２】
前記の開始点を囲む開始表面を規定するステップは：
複数の最大値投影（ＭＩＰ）を演算し；
各前記ＭＩＰにおいて、前記開始点に対応する血管の見かけの径を演算し；
最小見かけ径を選択し；かつ
前記開始点に中心化された、前記の選択された最小見かけ径の約１／２の半径の球体表面として前記開始表面を規定する；
ことを有することを特徴とする請求項１１に記載の磁気共鳴血管造影方法。
【請求項１３】
前記の表面が本質的に前記血管系の境界に一致するまで前記開始表面を反復的に成長させるステップは：
成長方向を算出し；
局部的なグレースケール情報に基づいて成長速度を積分し；
前記成長速度に対応する距離だけ前記成長方向へと前記表面を移動し；かつ
前記表面が前記血管系の前記表面に本質的に一致するまで、前記算出、前記積分、及び前記移動を反復する；
ことを有することを特徴とする請求項１１に記載の磁気共鳴血管造影方法。
【請求項１４】
前記の成長方向を算出するステップは：
少なくとも一つのヘシアン型マトリックスを算出し；
前記の少なくとも一つのヘシアン型マトリックスの固有値及び固有ベクトルを算出し；かつ
前記最小固有値に対応する方向として前記成長方向を算出する；
ことを有することを特徴とする請求項１３に記載の磁気共鳴血管造影方法。
【請求項１５】
前記反復成長の間、第１ブランチ及び少なくとも一つの追加ブランチを有する血管分岐点を対抗させ；
前記の少なくとも一つの追加ブランチに対応する位置をタグ付けし；
前記第１ブランチに沿って前記反復成長を継続し；
前記のタグ付けされた位置に戻り、かつ、前記の少なくとも一つの追加ブランチに沿って前記反復成長を繰り返し；かつ
全てのタグ付けされた位置が訪問されるまで前記戻りステップを繰り返す；
ことを有することを特徴とする請求項１１乃至１４のいずれか一項に記載の磁気共鳴血管造影方法。
【請求項１６】
エッジプレサービング三次元グレースケールスムージング法を用いて、前記高分解能容量画像表示をスムージングする；
ことをさらに有する請求項１乃至１５のいずれか一項に記載の磁気共鳴血管造影方法。
【請求項１７】
前記の少なくとも一つの血管識別画像スライスに対応するデータを取得するステップが含むのは：
複数の血管識別スライスに対応するデータを取得し、該血管識別スライスは複数の選択位置にて分散されている；
ことを特徴とする請求項１乃至１６のいずれか一項に記載の磁気共鳴血管造影方法。
【請求項１８】
（ｉ）高分解能画像容量に対応するデータ、及び
（ｉｉ）少なくとも一つの血管識別画像スライスに対応するデータ、
を有する、画像化データを取得するための手段であって、前記の少なくとも一つの血管識別画像スライスに対応する取得された血管識別データは、動脈及び静脈の一つに関して選択的に強調されたコントラストを有しており；
（ｉ）前記の取得された高分解能容量画像データから高分解能容量画像表示、及び
（ｉｉ）前記の少なくとも一つの血管識別画像スライスに対応した取得されたデータから少なくとも一つの血管識別スライス画像表示、
を再構築するための手段；並びに
前記血管識別スライス画像表示に基づいた前記高分解能画像において、動脈開始点及び静脈開始点の少なくとも一つを特定するための手段；
を有する磁気共鳴血管造影装置。
【請求項１９】
前記の血管識別スライス画像表示に基づいた動脈開始点及び静脈開始点の少なくとも一つを特定するための手段は：
前記血管識別スライス画像表示に基づいて、動脈及び静脈の一つに対応する第１血管型の少なくとも一つの開始点を特定するための手段；及び
前記血管識別スライス画像表示と前記高分解能画像の対応するスライスとの比較に基づいて、第２血管型に対応する少なくとも一つの開始点を特定するための手段；
を有することを特徴とする請求項１８に記載の磁気共鳴血管造影装置。
【請求項２０】
前記の画像データを取得するための手段は、前記高分解能画像容量内の異なる位置及び方向の少なくとも一つに配置されている複数の血管識別スライスに対応するデータを取得することを特徴とする請求項１８又は１９に記載の磁気共鳴血管造影装置。
【請求項２１】
関係する患者に磁気コントラスト剤を投与するための手段であって、前記の少なくとも一つの血管識別画像スライスに対応するデータは、コントラストが強調された磁気共鳴血管造影データを有することを特徴とする、手段；
をさらに有していることを特徴とする請求項１８乃至２０のいずれか一項に記載の磁気共鳴血管造影装置。
【請求項２２】
前記の画像化データを取得するための手段は、
動脈及び静脈の一つに関して選択的に強調されたコントラストを供するように選択的前飽和を用いて取得されたデータ；か
動脈及び静脈の一つに関して選択的に強調されたコントラストを供するように選択的前飽和を用いて取得されたタイム・オブ・フライトＭＲＡデータ；
かの一つを有する、少なくとも一つの血管識別画像スライスに対応するデータ；
を取得することを特徴とする請求項１８乃至２１のいずれか一項に記載の磁気共鳴血管造影装置。
【請求項２３】
前記の画像化データを取得するための手段は、前記の高分解能容量画像データの取得の間の前記の少なくとも一つの血管識別画像スライスに対応するデータの取得を分散することを特徴とする請求項１８乃至２２のいずれか一項に記載の磁気共鳴血管造影装置。
【請求項２４】
前記の特定された開始点における血管下位組織の開始を三次元で追跡するための手段；
をさらに有する請求項１８乃至２３のいずれか一項に記載の磁気共鳴血管造影装置。
【請求項２５】
前記の追跡手段は：
前記開始点を囲む開始表面を規定するための手段；及び
前記血管下位組織の境界に実質的に共通して延びるように、前記開始表面を反復的に成長させるための手段；
を有していることを特徴とする請求項２４に記載の磁気共鳴血管造影装置。
【請求項２６】
前記の画像化データを取得するための手段は、磁気共鳴画像化（ＭＲＩ）スキャナーであることを特徴とする請求項１８乃至２５のいずれか一項に記載の磁気共鳴血管造影装置。
【発明の詳細な説明】
【０００１】
（本発明の背景）
本発明は、画像化及び磁気共鳴技術に関する。特に、磁気共鳴血管造影に関し、それに関する特定の参照文献とともに述べる。しかしながら、本発明は、その他の管状構造又はネットワークの画像化への適用を見出すことができ、その中では、同様の管状構造及びネットワークが有意に分化されている。
【０００２】
血管造影法は、血管及び血管系の画像化に関する。血管造影法は、外科的計画及び処置の向上、診断の向上及び慢性血管疾患の非侵襲なモニタリングの簡便性を可能とし、かつ、動脈瘤や血栓などの潜在的な致死的状況の早期の傾向を供することができる。
【０００３】
血管造影法は、複数の異なる画像化様相を用いて行われ、これに含まれるのは、二方向性Ｘ線／ＤＳＡ、磁気共鳴（ＭＲ）、コンピューター断層撮影（ＣＴ）、超音波、及びこれら技術の種々の組み合わせである。磁気共鳴血管造影法（ＭＲＡ）は、コントラストが強調したモードにて行われてもよく、ここでは、ガドリニウムヂチレントリアミンペンタアセテート（Ｇａｄｏｌｉｎｉｕｍ−ｄｉｔｈｙｌｅｎｅ−Ｔｒｉａｍｉｎｅ−Ｐｅｎｔａ−Ａｃｅｔａｔｅ）を患者に投与し、患者のＭＲコントラストを向上させてもよいし、コントラストを強調しないモードにて行ってもよい。血管のコントラストは、典型的には、タイム・オブ・フライト（ＴＯＦ；ｔｉｍｅ  ｏｆ  ｆｌｉｇｈｔ）、ブラックブラッド（ｂｌａｃｋ−ｂｌｏｏｄ）、位相コントラスト（ｐｈａｓｅ  ｃｏｎｔｒａｓｔ）、Ｔ２、Ｔ２＊画像化などのＭＲ画像化技術を用いて血流を画像化することにより得られる。
【０００４】
ＴＯＦ法は、ＭＲＡにて広く行き渡っている。ＴＯＦ画像化手順は、典型的に、９０°のＲＦパルスを用いて第１組織スライスに磁気共鳴を励起するステップ、及び、第２スライスの近傍に１８０°位相の再集束ＲＦパルスを適用するステップを含む。血流速度及び９０°と１８０°とのＲＦパルス間の方向に関連した時間のずれが存在する。この時間のずれの間、第１スライスから第２スライスへと流れた血流は、９０°の励起パルス及び１８０°の再集束パルスの両方を受け、このＴＯＦ技術にて選択的に画像化されるスピンエコーを発生する。ＴＯＦと同様に多くのその他のＭＲＡ手法は、グレースケールの三次元画像を発生し、その中で、血管（又は血管内の血液）は、周囲の組織よりも、より明るく（白血球血管造影技術）又はより暗く（ブラックブラッド血管造影技術）見える。
【０００５】
未処理のグレースケールのＭＲＡ画像の解析及び解釈は、種々の因子により複雑化されている。人体における血管ネットワークが高度に複雑であるがゆえ、複雑になっている。特定の画像は典型的に、曲がりくねったり閉塞した血管や、可変な形状、Ｓ型のターン、非常に高度の血管密度を有する領域、種々の血管径の範囲、又はトラッキングを複雑にしている血管径におけるギャップなどを含んでいる。ＭＲＡデータの取得は、限定された動的範囲、部分的な容量の平均化及びこれに類似するものに起因して、追加的な複雑性を導入している。
【０００６】
ＭＲＡデータ解釈における特定の困難な問題は、血管のタイプの区別であり、つまり、画像化における動脈や静脈のことである。医療関係者はしばしば二つの方法を望む、つまり、動脈のみを表示する動脈図と静脈のみを表示する静脈図とである。しかしながら、多くのＭＲＡ技術は、動脈と静脈血管系とで感知可能な区別を供しない。先行技術において、血管型の識別は、典型的に、二つの技術を用いて達成されている：選択的なデータ取得並びに動脈及び静脈を区別するための取得後の画像処理である。
【０００７】
選択的なデータ取得は、通常、ＴＯＦ手法に適用されており、そこでは、静脈又は動脈のどちらかに由来する信号を抑制するため、画像化容量の上方又は下方のいずれかに前飽和スラブが配置されている。その他の選択的な取得法は、ＥＣＧゲートを使用し、かつ発生した画像データの位相の調光及び規模に基づいて、動脈図及び静脈図を区別するために、心臓の拍出サイクルの収縮及び拡張の間、血流速度をコード化する位相コントラスト法を使用している。ＴＯＦ及び位相コントラスト法が想定するのは、画像容量における血流方向が、限定的に静脈流と動脈流とを区別することである。一時的な区別もまた、コントラスト強調ＭＲＡにおいて提案されてきた。この手法は、動脈及び動脈血管系における循環時間の差異を利用している。画像化された容量の全領域に関してこの時間の差異を測定し感知することにより、動脈図及び静脈図を分離して取得することができる。
【０００８】
選択的なデータ取得法は、複数の欠点を有している。動脈図及び静脈図を分離して取得し、かつ、この取得時間が効果的に倍化するゆえ、走査時間が増加される。頭部やその他の身体の末端部などの場合などのように、画像化容量が、循環の分岐点を含む場合、この方法は、典型的に貧弱に行う。これらの分岐点の近傍では、静脈は、動脈図において、不適切に画像化され、及び／又は例えば、流入が抑制されるなど、静脈図は、貧弱なコントラストを有する。コントラストが強調したＭＲＡにおいて動脈と静脈とを区別することは、主として、第１と第２位相との間の時間のずれを注意深く選択することにより達成され、この画像の質は、このタイミングに非常に感受的になる。
【０００９】
その他の欠点は、多くの選択的な取得法に使用される前飽和パルスが、典型的に特異吸収率（Ｓｐｅｃｉｆｉｃ  Ａｂｓｏｒｐｔｉｏｎ  Ｒａｔｉｏ；ＳＡＲ）を増加することである。前飽和技術において上昇したＳＡＲは、特に、高い磁場強度において画像化するには問題となる。前飽和は、使用されるより高いサンプリング帯域を押しやる効果的なサンプルリングウィンドーを減少させる。前飽和法は、また、血流パターンの複雑性に起因して、ブラックブラッドをベースにしたＭＲＡにおいて非効果的である。
【００１０】
取得後の処理法には、使用者による手動の識別を含む。手動識別は、動脈系の初期動脈（ｐｒｉｍａｒｙ  ａｒｔｅｒｙ）又はルート動脈（ｒｏｏｔ  ａｒｔｅｒｙ）を好ましく横断する単一のスライスを使用してもよい。代替的に、スライス毎に使用者により解析される。手動法は、画像解析工程の遅延、この工程に対する人為的なエラーの要素の導入、及び正確な手動による血管系識別を行うための人体の画像化された領域に関する十分な解剖学的知識を備えた実務者が必要であるなどの欠点を有する。
【００１１】
動脈及び静脈の取得後の区別に関する自動化された方法は、典型的に特定のＭＲＡ画像化法にのみ適用するという制限された欠点を含み、例えば、位相コントラストＭＲＡやコントラスト強調ＭＲＡのみである。特に、これらの取得後処理法は、ブラックブラッドベースのＭＲＡ技術に互換性がない。ブラックブラッド技術は、白血球ベースのＭＲＡに対して利点を有し、この中で、ブラックブラッド技術は、典型的にＴＯＦ  ＭＲＡに対してより正確な血管形態を典型的に供する。
【００１２】
本発明は、上述の制限その他を克服する向上されたＭＲＡシステム及び方法を意図している。
【００１３】
（本発明の概略）
本発明の一面によると、磁気共鳴血管造影法を供している。高分解能容量画像データが取得され、高分解能画像容量に対応しているデータを備えている。少なくとも一つの血管識別画像スライスに対応するデータを取得し、この取得された血管識別データは、動脈及び静脈の一つに関して選択的に強調されたコントラストを有している。高分解能容量画像表示は、この取得された高分解能容量画像データにより再構築される。少なくとも一つの血管識別スライス画像表示は、少なくとも一つの血管識別画像スライスに対応する取得されたデータから再構築される。動脈開始点及び静脈開始点の少なくとも一つは、この血管識別スライス画像表示をベースにした高分解能画像にて特定される。
【００１４】
本発明の別の面によると、磁気共鳴血管造影装置を開示している。画像データ取得に関する手段を供しており、これには：（ｉ）高分解能画像容量に対応した高分解能容量画像データ；及び（ｉｉ）少なくとも一つの血管識別画像スライスに対応したデータ；を有しており、少なくとも一つの血管識別画像スライスに対応した、この取得された血管識別データは、動脈及び静脈の一つに関して選択的に強調されたコントラストを有している。（ｉ）取得された高分解能容量画像データに由来する高分解能容量画像表示；及び（ｉｉ）少なくとも一つの血管識別画像スライスに対応した画像データに由来する少なくとも一つの血管識別スライス画像表示；を再構築するための手段が供されている。血管識別スライス画像表示をベースにした、高分解能画像における動脈開始点及び静脈開始点の少なくとも一つを特定するための手段が供されている。
【００１５】
本発明の一つの利点は、二つの選択的な画像よりもむしろ、単一画像において、静脈図及び動脈図の両方を同時に取得することにより走査時間を向上している点である。
【００１６】
本発明のその他の利点は、本質的に、いかなる高分解能ＭＲＡ画像化技術にも適合することである。
【００１７】
本発明のその他の利点は、前飽和ＲＦパルスを阻止する点である。結果として、特異吸収率（ＳＡＲ）レベルが減少される。このことは特に、高磁場における画像化の場合に有利である。
【００１８】
本発明のその他の利点は、三次元画像スペースにおいて血管追跡を行うことによる血管追跡及び区別工程の間、三次元血管接続性なる利点を取得する点である。
【００１９】
本発明のその他の利点は、より高速により正確に血管区別及び分離を達しすべく、局部的なグレースケール情報に、方向性成長法（ｄｉｒｅｃｔｉｏｎａｌ  ｇｒｏｗｔｈ  ｍｅｔｈｏｄ）を組み合わせる点である。
【００２０】
本発明のさらにその他の利点は、静脈及び動脈が一貫した血流方向を欠除する部位において、動脈／静脈の効果的な分離性を向上する点である。
【００２１】
本発明のよりさらにその他の利点は、コントラスト強調ＭＲＡに関して、より簡便に血管分離性を有する点である。この区別は、先行技術に対して、一時的な情報及び画像化パラメータに対してより感受性が低い。動きに関連した、あるいはコントラストに関連したアーチファクトを変化する時間は、特定され、及び／又は修正される。
【００２２】
本発明のさらなる利点及び有益性は、好適実施例に関する以下の詳細な記載を読むことにより、当業者に明らかになるであろう。
【００２３】
（図面の簡単な説明）
本発明は、種々の構成要素及び構成要素の配置並びに種々のステップ及びステップの配置をとってもよい。図面は、好適実施例を示す目的のためのみであり、本発明を限定するものとして解釈されるべきではない。
【００２４】
（好適実施例に関する詳細な記述）
図１を参照すると、磁気共鳴画像化（ＭＲＩ）スキャナー１０は、検査領域１４を介するｚ軸に沿った、実質的に単一で、一時的に一定な主磁場Ｂ０を発生する超伝導性又は抵抗性磁石１２を有している。図１には、ボア型磁石が示されているが、本発明は、オープン型マグネットシステム、垂直型磁場システム及びその他のタイプのＭＲＩスキャナーに同等に適用可能である。磁石１２は、主磁場コントロール１６により制御される。画像化は、画像化される対象に対して、磁気共鳴（ＭＲ）工程を実行することにより行われ、例えば、磁気共鳴血管造影（ＭＲＡ）診察の間、患者を少なくとも部分的に検査領域１４内部に配置し、典型的に対象領域をアイソセンター（ｉｓｏｃｅｎｔｅｒ）に配する。
【００２５】
磁気共鳴工程は、シリーズのＲＦ及び磁場グラジエントパルスを必要とし、磁気スピンを反転又は励起し、磁気共鳴を誘導し、磁気共鳴を再集束し、磁気共鳴をマニピュレートし、部分的その他、磁気共鳴をコード化し、スピンを飽和し、並びに同様に操作すべく対象物に適用される。さらに特に、グラジエントパルス増幅器２０は、検査領域１４のｘ−、ｙ−及びｚ−軸に沿って磁場グラジエントを発生すべく、電流パルスをフォールボディーグラジエントコイルアッセンブリ２２に適用する。
【００２６】
好ましくはデジタル型のＲＦトランスミッター２４は、ＲＦパルスを検査領域に送信すべく、ＲＦパルス又はパルスパケットをフォールボディーＲＦコイル２６に適用する。典型的なＲＦパルスは、短時間に急速で連続的なパルスセグメントのパケットで構成されており、互いとともに及び種々の適用されたグラジエントが選択された磁気共鳴マニピュレーションを達成している。このＲＦパルスは、検査領域１４の選択された部分において、飽和し、共鳴を励起し、磁気化を反転し、共鳴を再集束し、又は共鳴をマニピュレートすべく使用される。
【００２７】
全身への適用に関して、選択されたマニピュレーションの結果として発生された得られる共鳴信号は、また、典型的に、フォールボディーＲＦコイル２６により取り出される。代替的に、対象の限定された領域においてＲＦパルスを発生することに関して、局部的なＲＦコイルは、一般に、選択された領域に連続的に配置される。例えば、本技術分野にて公知なように、挿入可能ヘッドコイル２８は、ボアのアイソセンターにおいて、選択された脳領域の周囲に挿入される。その他の表面コイル又はその他の特定なＲＦコイルもまた、任意で使用される。例えば、このＲＦシステムは、任意で位相化アレイレシーブコイル（図示せず）を有しており、これにより、本技術分野公知の部分的平行画像化（Ｐａｒｔｉａｌｌｙ  ｐａｒａｌｌｅｌ  ｉｍａｇｉｎｇ；ＰＰＩ）が使用される。好ましくは、フォールボディーＲＦコイル２６は、共鳴を誘導し、局部的なＲＦコイル又はコイルアレイは、選択領域から放射する磁気共鳴信号を受信する。その他の実施例では、この局部的なＲＦコイルは、得られる磁気共鳴信号を励起しかつ受信する。
【００２８】
ＲＦコイルの形状及びこの適用に関係なく、ＲＦコイルにより取得される、得られるＲＦ磁気共鳴信号は、ＲＦレシーバー３２により受信され、復調される。好ましくは、シーケンスコントロールプロセッサー３４は、一体化されたＭＲＩパルスシーケンスを発生し、磁気共鳴（ＭＲ）信号及び任意的なエコーを発生し、得られるＭＲ反応を空間的にコード化すべく適切なコード化グラジエントを供し、かつＭＲピックアップ及びレシーブオペレーションを調整する波型を読み出すべく、少なくともグラジエントパルス増幅器２０、ＲＦトランスミッター２４及びＲＦレシーバーを制御する。
【００２９】
ＭＲＩシーケンスは、典型的に、ＲＦコイル２６、２８により発生された選択されたＲＦパルスに沿って、ｋ−スペースに位置する磁気共鳴エコーをもたらす、グラジエント増幅器２０により発生される磁場グラジエントパルス及び／又は掃引の複合体シリーズを含む。
【００３０】
磁気共鳴血管造影の実験において、患者３６は、得られる画像において血管系を特に強調する画像化条件を用いてＭＲＩシステム１０により画像化される。例示的な図１において、患者３６の頸動脈部分が画像化されている。任意で、この患者は、例えば、コントラスト強調ＭＲＡを実施する、血管のコントラストを向上すべく、例えば、ガドリニウムジチレンペンタアセテートのボーラス注射などの磁気共鳴コントラスト剤３８を受ける。タイム・オブ・フライト（ＴＯＦ）シーケンスなどのＭＲＡシーケンス、ブラックブラッド血管造影シーケンス又はこれに類似する方法は、シーケンスコントロールプロセッサーにより適用される。好適なＴＯＦ実施例に関するシーケンスパラメーターには：ＴＥ＝６．７ｍｓ、ＴＲ＝２７ｍｓ、ＦＯＶ＝２０ｃｍ、ＰＥ＝１９２、ＮＳＡ＝１、フリップ角＝３５、ギャップ＝０ｍｍ、マトリックスサイズ＝３８４×５１２を含む。もちろん、当業者は、特定の画像化状況に関して、画像パラメーターコントロール及びウェイティング（ｗｅｉｇｈｔｉｎｇ）のタイプを用いて、よりより画質を発生すべく、このパルスシーケンスプログラムを調節してもよい。
【００３１】
例示的なＭＲＡシーケンスを上記したが、理解すべきことは、この開示された発明は、種々のＭＲＡ画像化手法に適合可能であり、これらには、限定されず、非コントラスト強調ＭＲＡ、コントラスト強調ＭＲＡ、ＴＯＦ、白血球ＭＲＡ、ブラックブラッドＭＲＡ、Ｔ２−ウェイティッドＭＲＡなどを含む。好ましくは、このＭＲＡ画像化法は、実質的に非選択的な様式にて、動脈及び静脈の両者を画像化する。
【００３２】
好適実施例において、血管識別に使用されるより多くのスライスの追加的なデータの取得が行われる。追加的なより多くのスライスに関してコード化された共鳴信号は、好ましくは、高分解能ＭＲＡ容量画像スライスに関する共鳴信号に分散される。好ましくは、この分散された血管識別画像スライスは、低分解能かつ高い走査速度でおこなわれており、したがって、患者全体に対する画像化時間に対する効果を無視できる。好適実施例において、より多くの血管識別スライスに関する追加的なデータ取得は、一つの血管型のみを選択的に画像化するシーケンスを用いて行われ、例えば、動脈のみを選択的に画像化したり、静脈のみを選択的に画像化したりする。
【００３３】
得られる磁気共鳴データは、３９において仕分けられ、高分解能３−Ｄデータを保存する高分解能ｋ−スペースメモリ４０にて保存され、かつ、少なくとも一つの血管識別スライスに対応するデータを保存する血管識別ｋ−スペースメモリ４２にて保存される。ｋ−スペースデータ４０、４２は、典型的にはフーリエ逆変換プロセッサー４４や本技術分野公知の再構築プロセッサーである再構築プロセッサー４４にて処理される。好適実施例において、この再構築プロセッサー４４は、画像メモリ４６に保存されている高分解能際構築画像を生成する。また、再構築プロセッサーは、血管識別画像メモリ４８に保存されている、血管識別画像に対応するより多くの分散されたスライスに関するより多くの血管識別画像を生成する。
【００３４】
この高分解能ＭＲＡ画像は、好ましくは、患者３６の他の体部組織に相対して、血管系に関して良好なコントラストを有する、患者の検討領域の三次元グレースケール画像表示である。しかしながら、この動脈及び静脈データは、実質的に互いに区別されず、高分解能ＭＲＡ容量画像表示４６は、これに含まれる血管情報に対して処理されない。より多くの血管識別スライスは、好ましくは、より低い分解能を有し、かつ、その他の血管型に対して動脈及び静脈の両者を選択的に画像化する。
【００３５】
取得後プロセッサー５０は、画像化された血管系についての追加的な情報を抽出及び／又は強調すべく、血管識別スライス画像４８と組み合わせて高分解能３−ＤＭＲＡ容量画像表示４６を処理する。この後処理は、好ましくは、静脈及び動脈系の自動的な分離、血管系及び血管ネットワークに関する情報の抽出並びにこれに類似するものを含む。動脈及び静脈の自動的分離は、有利である。なぜなら、一つ以上の分散された血管識別スライスは、キーを供し、これにより、高分解能画像に関する実質的に非区別な動脈及び静脈が分離されるためである。
【００３６】
得られる動脈及び静脈情報は、好ましくは、３−Ｄ動脈図メモリ５２及び３−Ｄ静脈図メモリ５４にそれぞれ保存される。この血管情報は、適切なユーザーインターフェース５６上に図示的に表示されてもよく、かつ、もちろん、電気的に保存されてもよく、あるいは、磁気記憶媒体上に保存されてもよく、紙上に印刷されてもよい（図示せず）。
【００３７】
続いて図１、そして、さらに図２を参照すると、分離した動脈図及び静脈図画像表示を生成するための工程に関する好適実施例が述べられている。高分解能容量、つまり、３−ＤのＭＲＡ画像は、実質的に非優先的な様式にて、動脈及び静脈の両者を画像化する公知の２−Ｄ又は３−Ｄ画像化法を用いて、ステップ７０にて取得される。例示的な実施例において、患者３６の頸動脈領域に関するブラックブラッドＭＲＡ画像化が使用されている。しかしながら、強調されるのは、本発明はこれに限定されるよりもむしろ、公知の二次元又は三次元ＭＲＡ画像化技術及び実質的にその他の種々の身体部分の画像化に適用可能であるということである。コントラスト剤３８は、血管のコントラストを向上すべく、任意で使用される。磁気共鳴画像化は、好ましくは、適切な画像再構築工程を用いて、三次元画像を形成すべく組み合わせ可能な、複数の実質的に空間的に近接したスライスを取得する通常の手法の後に続く。
【００３８】
加えて、少なくとも一つの血管識別スライスに対応するデータは、ステップ７２にて取得される。好ましくは、血管識別スライスに関するこのデータ取得は、高分解能３−Ｄ画像スライスに関するデータ取得に分散される。血管識別スライス画像化は、好ましくは、動脈又は静脈のどちらか一方のみを選択的に画像化する画像化法を用いる。取得７２はまた、好ましくは、高速であり、低分解能画像であってもよい。なぜなら、動脈ルート（ａｒｔｅｒｙ  ｒｏｏｔ）又は静脈ルート（ｖｅｎｏｕｓ  ｒｏｏｔ）などの主要血管を特定するためだけに用いられるからである。
【００３９】
好ましくは、この血管識別スライスは、標的とされた血管系のルート（ｒｏｏｔ  ｏｆ  ｔｈｅ  ｔａｒｇｅｔｅｄ  ｖａｓｃｕｌａｒ  ｓｙｓｔｅｍ）（動脈又は静脈）が位置している領域にて取得される。したがって、例示的な頸動脈領域に関して、取得７２は、好ましくは、毛動脈及び静脈のルートが配置されている頚部領域にて行われる。しかしながら、血管は、効果的に、例えば、狭窄点など、身体の種々の位置にて終結する可能性があるので、血管識別スライスは、スライスの選択方向に沿った複数の選択点に有利に分散されており、これは、血管識別スライス取得７２において、狭窄点などの分離された種々の血管系が、開始点に対応する少なくとも一つを有することを確実にするためである。
【００４０】
少なくとも一つの血管識別スライスに関するデータの取得７２は、好ましくは、動脈を覆う静脈を選択的に画像化するか、あるいは、動脈を好ましく覆う静脈を選択的に画像化する種々のＭＲＡ法によってもよい。この取得７２は、有利に、高速で低分解能の画像化法を用いる。好適実施例において、血管識別取得７２は、動脈又は静脈のどちらかを選択する選択的な前飽和を伴った２−Ｄタイム・オブ・フライト（ＴＯＦ）ＭＲＡ法を用いる。この手法を用いて取得されたスライスに関する画像化時間は、典型的には、数秒のオーダーであり、これは、患者全体の画像化時間と比較して無視可能なものである。しかしながら、再度強調されるのは、位相コントラスト、又はマグニチュードコントラスト法（ｍａｇｎｉｔｕｄｅ  ｃｏｎｔｒａｓｔ  ｍｅｔｈｏｄ）などを血管識別取得７２に関して用いてもよい。
【００４１】
高分解能３−Ｄ画像データの取得７０及び血管識別スライスデータの取得７２は、メモリ４０に保存された画像ｋ−スペースデータ７４及びメモリ４２に保存された血管識別ｋ−スペースデータ７６を、それぞれ生成する。画像再構築ステップ７８において、このｋ−スペースデータセット７４、７６は、高分解能画像メモリ４６に保存される高分解能３−Ｄ画像表示８０を発生すべく、かつ、血管識別画像メモリ４８に保存される少なくとも一つの血管識別スライス画像表示８２を発生すべく、再構築プロセッサー４４により再構築される。少なくとも一つの血管識別スライス８２及び高分解能３−Ｄ画像表示８０を用いて、ポストプロセッサー５０は、一つ以上の動脈ルート（ａｒｔｅｒｙ  ｒｏｏｔ）８６及び一つ以上の静脈ルート（ｖｅｉｎ  ｒｏｏｔ）８８を特定する。開始点８６、８８を用いて、血管系は、動脈図メモリ５２に保存される動脈図９４を発生すべく、かつ、静脈図メモリ５４に保存される動脈図９６を発生すべく、ステップ９２において分割され、かつ分離される。
【００４２】
続いて図２、そして、さらに図３を参照すると、動脈及び静脈を特定するステップ８４に関する好適な方法が述べられている。血管識別スライス１０２は、ステップ１００において、一つ以上の血管識別スライス８２にから選択される。血管識別スライス１０２は、動脈又は静脈の一つのみを選択的に有利に画像化する。図示された例示的な図３の実施例において、血管識別スライス１０２は、動脈を選択的に画像化するものとして示されている。血管識別スライス８２は、好ましくは、３−Ｄ高分解能画像８０のスライスのうちに分散されるので、血管識別スライス１０２に近接的かつ空間的に対応する高分解能画像８０から選択された高分解能スライス１０４が有意に存在する。例えば、頸動脈領域のＭＲＡが増加の場合、ｚ軸方向に垂直な軸方向に配置されたスライスが典型的に取得される。したがって、選択された血管識別スライス１０２に関して、血管識別スライス１０２としてｚ軸に調整されたのと同様又は近似的に同様なスライス１０４に関して、高分解能３−Ｄ画像を備えているスライスにより製造される。理解されるであろうことは、高分解能画像８０及びしたがって、スライス１０４に対応する高分解能は、図１及び２にて参照した前述のように、実質的に非分離的形態において、動脈及び静脈の両者を好ましく画像化する、ということである。
【００４３】
図３を続けて参照すると、選択された血管識別スライス１０２は、好ましくは、ここに含有される動脈を特定すべく、種々の適切な画像処理技術により解析される。血管識別スライス１０２は、好ましくは、一つの血管型のみを含んでいるので、この解析は、大部分、血管型区分が含まれないものとして、単純化される。この様式において、動脈ルート（ａｒｔｅｒｙ  ｒｏｏｔ）又は動脈マスク（ａｒｔｅｒｉｅｓ  ｍａｓｋ）８６が得られる。もちろん、静脈選択的な血管識別スライスの場合においても、静脈はステップ１０６に類似するステップにて、好ましく特定されるであろう。
【００４４】
ステップ１０８では、動脈及び静脈の両方を含む選択された高分解能３−Ｄスライス１０４は、ここに含まれる静脈を特定すべく、例えば、高分解能スライス１０４から動脈マスク８６を効果的に減じることによるなど、特定された動脈８６に比較して解析される。この様式において、静脈ルート（ｖｅｉｎ  ｒｏｏｔ）又は静脈マスク８８が得られる。
【００４５】
血管識別スライス取得７２は、高分解能画像の取得７０に分散されているので、高分解能スライス１０４を伴った血管識別スライス１０２の空間的な記録は、典型的に、スライス１０２、１０４の取得の間、有意に患者が動く可能性が低いこととして、進められる。画像の記録性をさらに向上すべく、好適実施例において、血管識別スライス１０２は、高分解能スライス１０４として、同様の（機械的許容性の範囲内の）平面に対応しており、つまり、例示的な頸静脈の画像化に関して、スライス１０２、１０４は、同様の軸平面に採取され、例えば、連続した交換可能なデータ取得を伴った同様のｚ軸座標値である。
【００４６】
ステップ１００、１０６及び１０８は、少なくとも一つの血管識別スライス８２のそれぞれに関して、好ましく繰り返される。この様式において、複数の動脈開始点８６及び静脈開始点８８が得られる。認識されるであろうことは、スライス８２の分散は、好ましくは、例えば、例示的な頸動脈ＭＲＡセッションに関して、より下方の頚部領域においてなど、血管ルート（ｖａｓｃｕｌａｒ  ｒｏｏｔ）の領域に向かって重点が置かれている。しかしながら、狭窄の可能性及びその他の血管追跡工程９２に対する困難性に起因して、前記の少なくとも一つの血管識別スライス８２は、好ましくは、不連続な脈管構造の追跡を可能にすべく、血管ルート（ｖａｓｃｕｌａｒ  ｒｏｏｔ）からはなれた追加的な開始点を供すべく配置されている複数のスライスを含んでいる。
【００４７】
血管型識別ステップ８４の代替的に意図された実施例において、複数の血管識別スライスが取得され、その中では、いくつかのスライスが静脈選択的である一方、複数のスライスが動脈選択的である。この代替的な実施例において、スライスの選択及び解析ステップ１００、１０６は、順番に、各血管識別スライスに関して繰り返し行われる。したがって、動脈及び静脈の両者に関する開始点及びルート（ｒｏｏｔ）は、３−Ｄ画像８０由来の対応する高分解能スライスの解析を含むステップ１０８に依存することなく得られる。
【００４８】
ここで、図２に戻って参照すると、動脈又は静脈に関する開始位置８６、８８がステップ８４においていったん特定されると、３−Ｄ動脈又は静脈ツリー（ａｒｔｅｒｉａｌ  ｏｒ  ｖｅｎｏｕｓ  ｔｒｅｅ）は、種々の適用可能な処理後技術を用いて、ステップ９２において血管を区分しかつ分離することにより推定（ｅｓｔｉｍａｔｅｄ）される。容量取得７０は、三次元の容量画像データ８０を生成するので、追跡は、好ましくは、直接的な３−Ｄ区分を構築するための３−Ｄ容量８０において得られる血管ツリー（ｖａｓｃｕｌａｒ  ｔｒｅｅ）についての３−Ｄ的な接続性情報に関する利点を取得する。
【００４９】
ここで、図４を参照すると、血管区分及び分離ステップ９２に関する好適実施例が述べられており、ここでは、開始位置８６、８８に対応するシード（ｓｅｅｄ）は、血管系を推定すべく、エッジプレサーブド容量において三次元的に成長される。最小限の抵抗性の通路に沿って流動体が流れている流動流に類似性において、この成長工程は、血管領域に対応する最小限の抵抗の通路の後に続く。好ましくは、エッジプレサーブド容量１２２は、エッジプレサーブド３−Ｄ形体的スムージング（ｓｍｏｏｔｈｉｎｇ）及びアイソトロピック容量生成ステップ１２０を用いて、高分解能３−Ｄ血管造影画像８０から生成される。３−Ｄにこの端部を保存することは、この３−Ｄの数学的な形態学を用いて生じ、したがって、この３−Ｄ血管接続情報をこのスムージング工程に有意に取り込む。
【００５０】
初期的な開始点８６、８８は、３−Ｄエッジプレサーブド容量１２２内部に初期的な血管表面位置を空間的に同定すべくステップ１２４において用いられる。図４は、動脈の例示的な追跡を示しており、したがって、この動脈開始点８６が適用される。もちろん、静脈系に関し、静脈開始点８８は、初期的な表面配置ステップ１２４に変わって適用されるだろう。一つの実施例において、この追跡工程９２は、動脈図９４及び静脈図９６を生成すべく連続的に繰り返される。その他の実施例では、動脈図及び静脈図は、実質的に平行して動脈及び静脈に追跡工程９２を適用することにより実質的に同時に演算される。
【００５１】
継続して図４、そして、さらに図５を参照すると、反復的な血管成長工程１２８を用いて、分割された血管ツリー１３０は、反復的に初期的な血管表面が成長することにより推定される。単純な成長技術において、各シードボクセル（ｓｅｅｄ  ｖｏｘｅｌ）の周囲の３Ｄ画素は、前もって選択された許容性の範囲内のグレースケールにて、これら３Ｄ画素に類似するものを同定すべく検索される。この工程は、反復的な成長を達成すべくグレースケールに類似の同定された３Ｄ画素のそれぞれの周囲の３Ｄ画素に関して繰り返される。典型的に状況において、例えば、動脈血管系ツリー１３０に関する図示された例示的な成長において、複数の血管下位組織（ｓｕｂ−ｓｙｓｔｅｍ）は、典型的に、同時にあるいは連続的のいずれかにて成長される。図５は、二つの例示的な血管下位組織（ｓｕｂ−ｓｙｓｔｅｍ）１３０Ａ、１３０Ｂの一部分を示しており、例えば、図３のステップ１０６に従った、動脈選択的な血管識別スライス１０２から抽出された情報に基づいて配置されている初期的な動脈表面１４０Ａ、１４０Ｂに由来して成長している。
【００５２】
典型的に、血管系は数学的に成長されるにつれ、例えば、図５に示されている例示的な分岐位置１４２Ａ、１４２Ｂなど、分岐位置に突き当たる。この分岐位置及びこれらの三次元的な位置を示しており、例えば、分岐のそれぞれのブランチ１４４Ａ、１４４Ｂの下位地位置における新規な初期表面１４０Ｃ、１４０Ｄや、この分岐のブランチ１４６Ａ、１４６Ｂに沿った連続的な成長などにより示している。示された３−Ｄ位置は、再度訪問され、かつ、この成長工程は、開始表面１４０Ｃ、１４０Ｄにおいて、再度開始される。これは、成長が、前出の追跡されていないブランチ１４４Ａ、１４４Ｂに沿って継続しているためである。
【００５３】
血管ブランチのタグ付け、一つのブランチに沿った継続及び続いてタグ付けされたブランチへの戻り、並びに成長工程１２４、１２６の再開工程は、全てのタグ付けブランチが訪問され、この血管ツリーが完全に追跡されるまで、続く。この成長工程は、シード１４０の公算において行われる。これは、このシステムが高速であり、粗野であるからである。この公算は、使用者に制御され、かつ、端部成長境界の両方の側面上にて３から５ボクセル層に変化してもよい。もちろん、また理解されるであろうことは、ステップ９２に実質的に類似の工程は、好ましく、静脈ツリーの追跡のために使用され、多重的な開始表面、分岐位置、及びこれに類似のものなどの状況は、同様に静脈を追跡する際、遭遇する。
【００５４】
ここで、戻って図４を参照すると、エッジプレサーブド容量１２０を生成する工程に関する好適実施例を述べている。典型的に血管造影的３−Ｄ容量８０は、等方性ボクセルを備えた容量を形成すべく、再度分割される。取得した血管造影データ８０は、通常、体部の常套的な直行軸、正中矢状方向及び前頭面に沿って配置される。頚部血管系のＭＲＡが増加に関する例示的な場合では、画像スライスは、典型的には、軸方向に配向されたスライスであり、位相コード化及び読み出し方向は、正中矢状方向及び前頭面方向の一つに沿って横たわる軸平面である。したがって、各ボクセルは、軸、正中矢状、及び前頭面平面方向の側に斜方晶系であるが、取得したボクセルは、典型的には立方体ではない。したがって、このボクセルは、立方体とすべく、等方性データを生成すべく、再分割される。当業者に認識されるであろうことは、同様の再分割は、体部に関して、軸、正中矢状方向及び前頭面に一致しない平面に沿って取得されたＭＲＡデータにも適用可能であるということである。
【００５５】
続けて図４を参照すると、好適実施例において、３−Ｄ端部保存システム１２０は、三次元において、等方性容量の覆う数学的で形態的な手段により行われる。このことは、３−Ｄスペースにおいて、膨張、侵食、開口、及び閉口などの基礎的な数学的形態学的手段を組み合わせて行われてもよい。
【００５６】
画像信号としてＩ、及び、構築要素としてＳを定義することにより、グレースケール膨張が好ましく二次元において定義される：
【００５７】
【数１】
2004535869.tif 000003
グレースケール膨張は、画像において形体的構造を膨張するという特性を有している。
【００５８】
このグレースケール侵食は、好ましくは、以下の二つの式にて定義される：
【００５９】
【数２】
2004535869.tif 000004
グレースケール侵食は、画像において形態学的構造が収縮するという特定を有している。
【００６０】
グレースケールオープニング（ｇｒａｙ  ｓｃａｌｅ  ｏｐｅｎｉｎｇ）及びグレースケールクロージング（ｇｒａｙ  ｓｃａｌｅ  ｃｌｏｓｉｎｇ）は、それぞれ式（１）及び式（２）で与えた膨張及び侵食なるグレースケール形態的制御を組み合わせて実行されてもよい。グレースケールクロージング制御は、好ましくは（以下の式（３））のように記述され：
【００６１】
【数３】
2004535869.tif 000005
かつ、グレースケールオープニング制御は好ましくは（以下の式（４）のように）記述される。
【００６２】
【数４】
2004535869.tif 000006
図６を参照すると、２−Ｄにおけるグレースケールの数学的形態的な基礎的な上述の４つの方程式は、エッジプレサーブドスムージングされた３−Ｄ容量１２２を生成するステップ１２０を実行すべく三次元に拡張されてもよい。これは、当業者公知である、式（１）から式（４）の二次元的表現を三次元的表現へと拡張することにより行われる。例示的な３−Ｄ構造要素１５２に沿った形態学的スムージングの前の、等方性３−Ｄ容量１５０の例示的な部分が図６に示されている。上述の技術の有意な利点は、ノイズは本質的に除去される一方、非常に小型の血管が典型的に保存されるという点である。この制御は３−Ｄであるので、血管のトポロジーは、このノイズの除去の間保存され、かつ、３−Ｄ血管接続性は、３−Ｄにおいてエッジプレサーブドスムージングを実行することにより有意に取り込まれる。
【００６３】
続けて図４、そして、さらに図７を参照すると、初期的な血管表面を生成し、かつ３−Ｄオイラー型表現１２６を演算するステップ１２４に関する好適実施例が述べられている。この初期表面は、スムージングされた三次元容量１２２において、通学的な血管成長に関するシードとして機能している。好ましくは、円型表面は、動脈ルート（ａｒｔｅｒｙ  ｒｏｏｔ）８６の３−Ｄ位置に挿入される。（もちろん、静脈の追跡に関し、静脈ルート（ｖｅｉｎ  ｒｏｏｔ）８８は、代替的に使用されるであろう）。複数の動脈ルート（ａｒｔｅｒｙ  ｒｏｏｔ）８６に関する最大値投影（ＭＩＰ）は、ステップ１６０において構築される。このＭＩＰ図１６２は、好ましくは、複数の正投影に沿って得られ、各ＭＩＰ図１６２において特定された動脈ルート（ａｒｔｅｒｙ  ｒｏｏｔ）の位置を伴っている。動脈ルート（ａｒｔｅｒｙ  ｒｏｏｔ）位置座標（ｘ、ｙ、ｚ）の中心に位置している３−Ｄ円型の半徑は、以下のように、このＭＩＰ図から演算される。
【００６４】
動脈ルート（ａｒｔｅｒｙ  ｒｏｏｔ）８６の見かけの径は、複数の見かけの血管径１６６を発生すべく、ステップ１６４において各ＭＩＰ図にて算出される。最小見かけ半径ｄｍｉｎは、１６０を選択され、かつ、円型表面１７は、ｒ＝ｄｍｉｎ／２なる半徑を有するルート動脈（ａｒｔｅｒｙ  ｒｏｏｔ）座標８６に形成した円型として定義されている。もし、複数のＭＩＰ図又は座標（ｘ、ｙ、ｚ）に配置された動脈ルート（ａｒｔｅｒｙ  ｒｏｏｔ）８６の投射１６２がｎである場合、円型１７０の半径ｒは、好ましくは、以下のとおり演算される：
【００６５】
【数５】
2004535869.tif 000007
ここで、ｄ１、ｄ２、・・・ｄｎは、ｎ番目のＭＩＰ図１６２における血管ルート（ｖａｓｃｕｌａｒ  ｒｏｏｔ）８６の見かけの半径１６６である。初期円型表面１７０の算出において、最小見かけ半径ｄｍｉｎの適用とともに複数のＭＩＰ図１６２を使用することは、選択された円型表面１７９が不適切に血管領域の外側に横たわる領域を含む可能性を大いに減弱させる。以下に述べる数学的な血管成長が初期表面において開始するので、過大推定された初期表面は、典型的に初期表面の位置において、誤って大きい血管径を典型的にもたらすであろう。当業者は、３ｎのＭＩＰ図１６２を生成すべく、正中矢状、前頭、及び軸図データセットに対応するＭＩＰ投影１６２を使用してもよい。ｎに関する範囲は、６である。また、さらに、初期半径ｒの過大推定の可能性を減ずべく、追加的なＭＩＰ図を演算すべく、複数の回転を使用してもよい。
【００６６】
続けて図７を参照すると、初期球体１７０は、この球体１７０の示された距離トランスフォーム１７２を演算することによりオイラー的表現１２６に変換される。この示された距離トランスフォーム１７２は、工程であり、ここで、ボクセルは、球体１７０の端部から、＋ｖｅ又は−ｖｅなる距離値にて割り当てられている。例えば、球体１７０の第１内部層は、−ｖｅの一つの値を与えられている。同様の様式にて、この公算における全てのボクセルは、±ｄなる示された距離値を割り当てられており、ここで、ｄは、球体端部からの距離である。このトランスフォーム１７２の結果は、典型的に埋め込まれた球体１２６のオイラー的表現と呼ばれる。
【００６７】
ここで図８を参照すると、３−Ｄ成長工程１２８による反復的な表面推定に関する好適実施例が述べられている。このステップ１２８は、実質的にボクセルの接続性及び３−Ｄグレースケール容量分布に三次元的に基づいた初期表面１７０の成長を有している。この反復的な表面工程は、二つの主要な構成要素を有している：一つは、方向性構成要素１８０であり、もう一つは、表面が成長する位置における速度に対応している積分速度構成要素１８２である。方向性構成要素１８０は、演算されたヘシアン型マトリックス（Ｈｅｓｓｉａｎ  ｍａｔｒｉｃｅｓ）の固有ベクトルを用いて演算される。この三次元のヘシアン型マトリックスの最小固有値は、３−Ｄ成長工程演算の方向性構成要素１８０を与える。好適実施例において、この方向性構成要素１８０に関する算出におけるノイズの妨害は、方向性演算１８４の前に、ガウス型オペレーター（図示せず）のより高いオーダーの微分（ｈｉｇｈｅｒ  ｏｒｄｅｒ  ｄｅｒｉｖａｔｉｖｅ）による幾何学的容量１２２を畳み込む（ｃｏｎｖｏｌｖｉｎｇ）ことにより限定される。
【００６８】
成長及び積分された速度１８２は、表面又は成長容量の周囲に存在する領域に関する端部、領域及びグラジエント情報を推定することにより演算される。端部、領域及びグラジエント速度構成要素は、成長工程に関するレベルセットフレームワークに積分される。３−Ｄ表面１８８の反復的な伝播は、３−Ｄエッジプレサーブド容量１２２における血管径を数学的に成長すべく、演算され積分された速度１８２に従った方向に生じる。この様式において、方向性成長は、この成長速度を制御するグレースケールの局部情報と組み合わされる。従って、この方向性構成要素は、埋め込まれた球体１２６の初期オイラー的表現に由来する血管系を急速かつ正確に数学的に成長すべく、エッジプレサーブドスムーズド容量１２２に含有される統計的な端部情報及びグラジエント情報に融合される。
【００６９】
３−Ｄ血管が成長するので、例えば、図５に示した分岐点１４２Ａ、１４２Ｂなどの分岐点は、典型的には、遭遇される。図５を参照して上述したように、この分岐点、例えば１４４Ａ、１４４Ｂなど、の一つのブランチの３−Ｄ位置は、タグ付けされ、成長は、その他のブランチ、例えば１４６Ａ、１４６Ｂなどへと下方的に継続する。このタグ付けされた位置は、再度訪問され、この成長工程１２４、１２８は、順番にタグ付けされた点のそれぞれにおいて再開される。
【００７０】
成長工程は、表面の公算に限定されているので、成長工程１２８は、速く、粗野である。この公算は、好ましくは、ユーザーに制御され、かつ、典型的には、端部成長境界のそれぞれの側において３から５ボクセルにて、変化する。
【図面の簡単な説明】
【００７１】
【図１】本発明の好適実施例による、例示的な磁気共鳴血管造影装置を示している。
【図２】本発明の好適実施例による、血管造影画像再構築システムに関する概要を示している。
【図３】本発明の好適実施例による、動脈及び静脈の区別方法を示している。
【図４】本発明の好適実施例による、血管区別及び分離のための方法を示している。
【図５】血管識別スライスにおける多重的な動脈開始点と同様に二つの分岐点を有する例示的な血管系を示している。
【図６】グレースケールスムージング工程に関する例示的な幾何学的形状を示している。
【図７】本発明の好適実施例による、血管系開始点に対応する開始円型表面に関する初期的なオイラー的表現を得るための方法を示している。
【図８】本発明の好適実施例による、エッジプレサーブド容量（ｅｄｇｅ−ｐｒｅｓｅｒｖｅｄ  ｖｏｌｕｍｅ）における血管系を反復的に成長させるための方法を示している。
【図１】
2004535869.tif 000008
【図２】
2004535869.tif 000009
【図３】
2004535869.tif 000010
【図４】
2004535869.tif 000011
【図７】
2004535869.tif 000012
【図８】
2004535869.tif 000013
【国際公開パンフレット】
2004535869.tif 000014
2004535869.tif 000015
2004535869.tif 000016
2004535869.tif 000017
2004535869.tif 000018
2004535869.tif 000019
2004535869.tif 000020
2004535869.tif 000021
2004535869.tif 000022
2004535869.tif 000023
2004535869.tif 000024
2004535869.tif 000025
2004535869.tif 000026
2004535869.tif 000027
2004535869.tif 000028
2004535869.tif 000029
2004535869.tif 000030
2004535869.tif 000031
2004535869.tif 000032
2004535869.tif 000033
2004535869.tif 000034
2004535869.tif 000035
2004535869.tif 000036
2004535869.tif 000037
2004535869.tif 000038
2004535869.tif 000039
2004535869.tif 000040
2004535869.tif 000041
2004535869.tif 000042
2004535869.tif 000043
2004535869.tif 000044
2004535869.tif 000045
2004535869.tif 000046
2004535869.tif 000047
2004535869.tif 000048
2004535869.tif 000049
2004535869.tif 000050
2004535869.tif 000051
2004535869.tif 000052
2004535869.tif 000053
2004535869.tif 000054
【国際調査報告】
2004535869.tif 000055
2004535869.tif 000056
2004535869.tif 000057
2004535869.tif 000058