(57)【特許請求の範囲】
【請求項１】
モデルをベースとした画像分割の際に３Ｄ面のモデルを画像特性に自動的に適合させるためにデータ処理装置において使用される画像データ処理方法であって、
指定された一組の点で構成された３Ｄパスに適応する変形可能な管状のメッシュモデルを作る段階と、
前記３Ｄパスの曲率、前記３Ｄパスにおける点同士のサンプル距離及び所定の入力半径に基づいて、前記変形可能なメッシュモデルの半径を自動的に適合させる段階と
を有し、前記３Ｄパスに適応する変形可能な管状のメッシュモデルを作る段階が、
セグメントに対する管状の関心対象物の中心線に対応する前記３Ｄパスを算出し、前記３Ｄパス上のセグメントを規定する段階と、
前記３Ｄパスの長さに等しい長手方向軸に沿って規定された長さを有する最初は真っ直ぐで変形可能な円筒形のメッシュモデルを作る段階と、
最初は真っ直ぐな前記メッシュモデルを、前記３Ｄパスの個々のセグメントに関連した長さのセグメントに分割する段階と、
前記メッシュモデルの各セグメントについて、前記最初は真っ直ぐで変形可能な円筒形のメッシュモデルとして規定された前記メッシュモデルの最初の方向を、前記３Ｄパスにおいて対応するセグメントの方向に変換する剛体変換を算出し、及び前記セグメントに対応する前記メッシュモデルの頂点に前記剛体変換を適用する段階と
を有する、画像データ処理方法。
【請求項２】
連続するセグメントに関連する剛体変換をさらに算出する段階を有し、該剛体変換は、２つの連続するセグメントの間で融合される変換であり、前記メッシュモデルの屈曲部分の間の自己交差を制限する、請求項１記載の画像データ処理方法。
【請求項３】
連続するセグメントの間の剛体変換用の回転を算出する段階をさらに有し、
線形補間が、３Ｄ剛体変換融合に関する２つの回転の間で使用され、前記メッシュモデルの屈曲部分の間の自己交差を制限する、請求項２記載の画像データ処理方法。
【請求項４】
メッシュの捩れを最低限化するように、前記最初は真っ直ぐで変形可能な円筒形のメッシュモデルの方向である最初のメッシュモデル方向から、ターゲットである前記３Ｄパスのセグメントへの最低限の３Ｄ回転を算出する段階を有する、請求項１ないし３のうちいずれか一項に記載の画像データ処理方法。
【請求項５】
軸パラメータ及び回転角度パラメータを利用して、２つのセグメント間の回転を規定する段階を有し、当該画像データ処理方法は、現在のセグメントに対する新しい回転が、以前のセグメントに対して算出された回転と、前記以前のセグメント及び前記現在のセグメントの間の最低限の回転との複合として算出されるように、１つのセグメントから他のセグメントへ反復的にパラメータを算出する段階をさらに有する、請求項４記載の画像データ処理方法。
【請求項６】
医療用観察システムであって、
ユーザが当該医療用観察システムとやり取りできるように動作する入力装置と、
映像化装置と、
前記入力装置及び前記映像化装置に関するデータ処理手段と
を有し、前記データ処理手段は、請求項１ないし５のうちいずれか一項に記載の画像データ処理方法を実行するようプログラムされている、医療用観察システム。
【請求項７】
体の器官の三次元画像を取得する手段と、
請求項６による医療用観察システムと
を有する医療用検査装置。
【請求項８】
請求項１ないし５のうちいずれか一項に記載の画像データ処理方法を実行する一組の命令を有する、コンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医療用画像における物体の略管状の面に三次元（３Ｄ）の変形可能なメッシュモデルを自動的に適合させる画像処理方法に係る。本発明は、三次元画像で三次元の管状の物体の分割の方法に適用し、三次元の別個の変形可能なモデルを前述の三次元の管状の物体に適合させる操作をしばしば有する。本発明は更に、結腸又は動脈等の体器官である物体の分割は、臓器の病状を研究又は検出をするようにされた物体の分割に関して、医療用画像装置若しくはシステムに、又は、それらの装置若しくはシステムによって作られた医療用三次元画像を処理するプログラムに係る。本発明は、医療用画像検査方法の分野、プログラム、及び、装置又はシステムにおいて、特定な用途を見つける。
【背景技術】
【０００２】
３Ｄ物体のかたどりの技術は、既に、「processing  of  the  International  Conference  on  Computer  Vision  and  Pattern  Recognition  (CVPR'94)、１９９４年６月２０日−２４日、米国シアトル」会議での出版物「Simplex  Meshes:  a  General  Representation  for  3D  shape  Reconstruction」（非特許文献１）で、H.  Delingetteによって開示される。かかる文献では、三次元の物体を再生する物理的なアプローチが提起される。かかるアプローチは、「単一メッシュ（Ｓｉｍｐｌｅｘ  Ｍｅｓｈ）」の形状に基づく。メッシュの弾性的挙動は、各頂点（メッシュの結節）で抽出された単一角度を介して平均曲率を制御する、局所的安定化機能によってかたどられる。それらの機能は、視座が不変であり、内因性で、スケール依存である。単一メッシュは、不変の頂点持続性を有する。３Ｄ面を示すには、双単一メッシュ（ｔｗｏ−Ｓｉｍｐｌｅｘ  Ｍｅｓｈ）と呼ばれる単一メッシュが使用され、各頂点が３つの近接する頂点と接続される。単一メッシュの構造は、非特許文献１に関し図１が例示する通り、三角形の構造に対して二元的である。単一メッシュの輪郭は、閉多角形鎖として定義され、単一メッシュ上の近接する頂点から構成される。４つの独立した変換は、可能なメッシュ変換の全範囲を達成するよう定義される。それらは、面でのエッジ部の挿入又は削除にある。単一メッシュの説明は、また、平面形状で使用される角度を一般化した単一角度の定義と、頂点が他の２つの近接する頂点に対してどのように位置づけられるかを説明する測定パラメータの定義を有する。各頂点のダイナミックは、ニュートンの運動の法則によって与えられる。変形とは、形状が平滑になるようにさせる力、及び、メッシュを３Ｄの物体に近づけさせる力を示す。内力は、外的制約に対する物理的なモデルの応答を決定する。内力は、それらが内因性の視座が不変でスケール依存であるように表される。類似した種類の制約は、輪郭を維持する。従って、引用された文献は、所定の３Ｄ物体を示すよう単純なモデルを与える。力は、関心対象の３Ｄ物体上でモデルを再形成し調整するよう、適用されると定義される。
【０００３】
医療用画像においては、結腸又は動脈等の管状の器官を分割することが頻繁に求められる。離れた変形可能なモデルに基づく分割は、研究された器官の直径又は容量等の臨床的指標の抽出を可能にする。問題が起きるのは、双単一メッシュと称される種類のモデル、三角形メッシュ、又は他のいかなる種類のアクティブ輪郭のモデルに関わらず、離れた変形可能なモデルが、非常に多くの屈曲部を示すかかる管状の器官に適合しなければならないときである。かかる屈曲部は、小さい又は大きい曲率半径によって定義された異なる曲率を有し得る。
【非特許文献１】H.  Delingette、「Simplex  Meshes:  a  General  Representation  for  3D  shape  Reconstruction」、（米国）、１９９４年６月２０日−２４日
【発明の開示】
【発明が解決しようとする課題】
【０００４】
本発明は、変形可能なモデルを作る画像処理方法を提案することを目的とする。モデルは、双単一メッシュ又は三角形メッシュ又は他の種類のメッシュを備え、指定された一式の点で構成された以前に定義付けられた３Ｄパスを適合させる管状の構造を有する。かかる３Ｄパスは、関心対象の物体、即ち屈曲部を示す管状の器官の中心線に対応し得る。従って、前出の３Ｄパスは、全ての種類の曲率を示す。
【０００５】
本発明は更に、３Ｄの変形可能な管状モデルを、異なる屈曲部を示す関心対象の管状の物体の３Ｄ面へ自動的にマッピングする、画像処理方法を提案することを目的とする。関心対象の物体は、グレーレベルの３Ｄ画像で示され得る。かかる方法は、分割する関心対象の管状の物体の中心線に対応する３Ｄパスを算出する段階と、前出の３Ｄパスのセグメントを定義付ける段階と、３Ｄパスの長さと同一の長手方向軸に沿って定義された長さを備えた最初の真っ直ぐで変形可能な円筒形のメッシュモデル、又は他のいかなる種類のメッシュを作成する段階と、かかる最初のメッシュモデルを３Ｄパスの異なるセグメントに関連する長さのセグメントに分割する段階と、メッシュの各セグメントに関しメッシュの最初の方向を３Ｄパスの関連するセグメントの方向に変換する剛体変換を算出する段階と、この変換をセグメントに対応するメッシュの頂点に適用する段階と、を有する。
【課題を解決するための手段】
【０００６】
ここで、主要な問題は、３Ｄの２つの連続するセグメントの間の方向が素早く変化するとき、３Ｄパスの屈曲部によるメッシュモデルの自己交差の潜在的な出現にある。
【０００７】
本発明は、変形可能なメッシュモデルの自己交差を制限することができる画像処理方法を提案することを目的とする。自己交差を制限するにあたり、提案された方法は各セグメントに対して固有の変換を使用しないが、その代わりに、提案された方法は、続くセグメントに関連した変換を算出する段階を有し、変換は２つの連続するセグメントの間で融合される。この技術は、３Ｄ剛体変換融合と称される。回転がセグメントの方向に従うよう使用されるとき、線形補間が、回転融合の２つの回転の間で使用され得る。１つの３Ｄ剛体変換から他へのかかる線形補間は、四元数の形式主義を使用して解決され得る。
【０００８】
かかる３Ｄ剛体変換融合の技術の利点は、提案された方法に従ったメッシュモデルの作成が、他の種類のものよりもメッシュモデルの種類を与えることを必要としないことである。提案された方法は、双単一メッシュ、三角形メッシュ、又は他の種類の変形可能なモデル等のいかなる種類のメッシュモデルにも適用され得る。
【０００９】
１つのセグメントから他のセグメントへの３Ｄ剛体変換の線形補間は、自己交差を完全に防ぐのに常に十分ではない。明らかに、かかる自己交差もまた、３Ｄパスの局所的曲率と作られた離れている変形可能なモデルの所望される半径との間の関係に依存する。後者が曲率の局所的な半径より大きい場合は、自己交差が発生する。
【００１０】
本発明は、屈曲部を示す管状の離れたモデルの領域で可能性のある自己交差を防ぐことが出来る画像処理方法を提案することを目的とする。本発明の方法は、局所的曲率に従って、円筒形の変形可能なモデルの半径を調整する段階を有する。最初に、ターゲットである３Ｄパスは、単なる点のリスト異常のもではあり得ないため、局所的曲率が概算される。続いて、メッシュモデルの半径は、曲率の機能で限定される。しかしながら、モデルの１つのセグメントから他のセグメントへの急な半径の変化を避けるよう、線形融合、又は、１つの半径から他への双三次スプライン補間等の半径調整の技術が、メッシュモデルの半径に適用される。
【００１１】
半径調整のかかる技術の利点は、メッシュモデルの表面が、半径調整によって得られた限定された半径の部分において、いかなる断続も示さないことである。面は、所定の半径の一部分から限定された半径の他の部分まで平滑である。従って、メッシュモデルは変形されない。半径調整のかかる技術の他の利点は、関心対象の複雑な器官の改善された分割を多数の折畳みが映像化又はかかる器官の全部分を追うことを難しくし得る領域において、可能にすることである。映像化は改善される。
【００１２】
半径調整のかかる技術の他の利点は、メッシュモデルの半径が限定されるため、提案された方法に従ったメッシュモデルの作成が、屈曲部の自己交差部分への切断の適用を必要としないことである。切断という防がれた代替案は、当業者は自己交差する屈曲部の部分に適用することを知っている種類の操作である。かかる切断は、自己交差を防ぐ一解決方法であるが、曲率の領域での情報の抑圧及び不自然な突然の折畳み等の問題点を示し、それらを本発明が防ぐ利点を有する。
【００１３】
更に他の問題は、上述の変換の連続性制御の欠落であり得、モデルのメッシュ捩れとして見られ得る。
【００１４】
本発明は更に、メッシュ捩れを最低限化することが出来る画像処理方法を提案することを目的とする。メッシュ捩れは、剛体変換の２つの連続する回転の間の距離が最小限であるときに、最低限化される。提案された方法は、最初のメッシュ方向からターゲットセグメントへの最小限の３Ｄ回転を算出する段階を有する。この回転は、軸のパラメータ及び回転角度のパラメータに沿って定義される。提案された方法は、これらのパラメータを一セグメントから他へ反復して算出する段階を有するため、最新のセグメントの新しい回転は、以前のセグメントに見つけられた回転と、以前及び最新のセグメントからの最小限の回転との複合として算出される。
【００１５】
提案された画像処理方法の主要な特徴が、以下において説明される。本発明の方法は、望ましくは、管状の体器官の三次元医療画像の分割に適用され得る。本発明はまた、３Ｄ画像処理手段を有する医療用検査撮像装置に係る。医療用撮像装置は、撮像手段を有するいずれの３Ｄ医療用検査装置でもあり得る。本発明は、更には、方法を実行するプログラム又はプログラム・パッケージに係る。
本発明の一形態による画像データ処理方法は、
モデルをベースとした画像分割の際に、３Ｄ面のモデルを画像特性に自動的に適合させるためにデータ処理装置において使用される画像データ処理方法であって、
指定された一組の点で構成された３Ｄパスに適応する変形可能な管状のメッシュモデルを作る段階と、
前記３Ｄパスの曲率、前記３Ｄパスにおける点同士のサンプル距離及び所定の入力半径に基づいて、前記変形可能なメッシュモデルの半径を自動的に適合させる段階と
を有し、前記３Ｄパスに適応する変形可能な管状のメッシュモデルを作る段階が、
セグメントに対する管状の関心対象物の中心線に対応する前記３Ｄパスを算出し、前記３Ｄパス上のセグメントを規定する段階と、
前記３Ｄパスの長さに等しい長手方向軸に沿って規定された長さを有する最初は真っ直ぐで変形可能な円筒形のメッシュモデルを作る段階と、
最初は真っ直ぐな前記メッシュモデルを、前記３Ｄパスの個々のセグメントに関連した長さのセグメントに分割する段階と、
前記メッシュモデルの各セグメントについて、前記最初は真っ直ぐで変形可能な円筒形のメッシュモデルとして規定された前記メッシュモデルの最初の方向を、前記３Ｄパスにおいて対応するセグメントの方向に変換する剛体変換を算出し、及び前記セグメントに対応する前記メッシュモデルの頂点に前記剛体変換を適用する段階と
を有する、画像データ処理方法である。
【発明を実施するための最良の形態】
【００１６】
本発明は、以下の図面を参照して下文に説明される。
【００１７】
本発明は、グレーレベルで示される三次元（３Ｄ）デジタル画像等に適用される画像処理方法に係る。画像は、関心対象の物体とされる器官のノイズのある三次元面を示し得る。ユーザが関心対象の物体を例えば背景に対してよりよく見えるにするよう、かかる物体は分割される。分割された画像によって、ユーザは、器官の異常又は病気をよりよく研究又は検出できる。本例では、関心対象の物体は、結腸又は動脈又は気管支等の管状の器官である。三次元の画像が得られる方法は、本発明の一部ではない。分割の方法は、超音波システム若しくはＸ線装置によって、又は当業者に既知である他のシステムによって得られ得る器官の三次元デジタル画像に適用され得る。
【００１８】
三次元の管状である関心対象の物体を示す三次元画像の獲得後、前出の画像は分割される。画像分割方法は、アクティブ輪郭と称される３Ｄの変形可能なモデルの活用に基く。本発明によれば、３Ｄの変形可能なモデルを作成するいかなる技術も、限定無く使用され得る。分割の操作は、３Ｄの変形可能なモデルを、結腸又は動脈等の３Ｄの略管状である関心対象の物体にマッピングすることにある。本例では、関心対象の物体は、多数の屈曲部を有する複雑な管状形を示す。
【００１９】
アクティブ輪郭のフィールドでは、最初のメッシュモデルが与えられなければならない。メッシュモデルのいかなる任意形状からも開始することが常に可能であるとしても、分割されるべき器官の望ましい形に形が近いメッシュモデルから開始したほうがより着実でより早い。本発明によれば、単一メッシュと称される種類の最初のメッシュモデル、三角形メッシュ、又はメッシュモデルの他の種類を作成することが提案され、一組の指定された点が形成する３Ｄパスを適合させる管状構造を伴うようにされる。この方法は、長く真っ直ぐな円筒の作成に基づき、円筒は同様に３Ｄパスを適合させるよう屈曲又は変形される。
【００２０】
困難な点は、形の複雑性を鑑みると、かかる管状の体器官の複雑な形の面全体を最終的に正確にマッピングするよう、真っ直ぐな最初の管状の変形可能なモデルを適切に変形する作業にある。
【００２１】
１）  管状のメッシュモデルを作る段階：
図１によれば、血管又は結腸等の管状構造の分割は、最初に３ＤパスＰを作ることを有する。かかる３Ｄパスは、有利には、かかる管状構造に対応する器官の中心線であり得る。かかる３Ｄパスの作成は、当業者には既知であるパス・トラッキング技術を使用して実施され得、管状構造に関連する一組の指定された点を与えるようにされる。本発明によれば、３ＤパスＰは、管状のメッシュモデルに転用される。最初のメッシュモデルは、円筒形のメッシュモデルである。かかる最初の円筒形のメッシュモデルは、管状の器官構造のデータに向かって変形される。これに関して、従来技術の出版物（Ｄｅｌｉｎｇｅｔｔｅ）にある通りメッシュモデルを直接物体の面とともに初期化する代わりに、３ＤパスＰ等からメッシュモデルを初期化するよう機能性が必要とされる。
【００２２】
冒頭で述べた通り、管状構造を分割しようとする用途はいずれも、管状形を有する最初のメッシュモデルから利益を得ると考えられる。パス・トラッキングのツールは、既に当業者には既知であるため、ツールは、分割されるべき関心対象の管状の物体の中心線等の関心対象の物体に関連する点の線を決定するよう使用され得る。かかる３ＤパスＰ又は中心線は、メッシュモデルの初期化工程を信頼性高く案内することができ、より早いだけではなくより堅固な環状構造用の分割ツールをもたらす。
【００２３】
本発明によれば、解決すべき問題は、分割されるべき管状の器官を適合させる管状のメッシュモデルの作成である。問題である異なる入力は以下の通りである：
１）  ３ＤパスＰに沿ってある点のソートされたリスト。これらの点の規則性及び間隔に概算は求められないが、かかる制約が平滑なメッシュモデルを得る十分な助けとなり得ることが更に明示されるであろう。
【００２４】
２）  円筒形の半径ｒ。
【００２５】
３）  セルの分解、即ち「セルの再分割の数」。
【００２６】
自然出力はメッシュ構造である。
【００２７】
図１によれば、円筒の基本形を作成する機能性が提案される。この機能性は、予め定義された参照符号Ｏｘ，Ｏｙ，Ｏｚでできたｚ軸に沿って、最初の円筒形のメッシュモデルの円形断面上にある一組の点を作成する段階と、続いて単一メッシュ構造を作成するよう一連の点を全てつなげる段階と、最後に、最初はｚ軸上に整列された作成された円筒を望ましい３Ｄ位置まで動かすよう剛体変換を適用する段階とを有する。一般化された円筒と称される３Ｄの可撓の管を作成するにあたり、本発明の方法は、真っ直ぐの円筒から開始する段階を有し、円筒はｚ軸上に整列され、３Ｄターゲット・パスＰの全長と同一である長さｓを有する。続いて、かかる方法は、所定の３Ｄパスを適合させるよう、かかる円筒を弾性的に歪める段階を有する。
【００２８】
図１によれば、本方法はまず、
３Ｄパスは分割される関心対象の管状の物体の中心線である等の、関心対象の管状の物体に関連する３Ｄパスを算出する段階と、
かかる３Ｄパス上のセグメントを定義付ける段階と、
最初の真っ直ぐで変形可能な円筒形のあらゆる種類のメッシュモデルを、３Ｄパスの長さと同一の長手方向軸に沿って定義された長さを備えて作成する段階と、
かかる最初のメッシュモデルを３Ｄパスの異なるセグメントに関連する長さのセグメントに分割する段階と、
メッシュの各セグメントに関し、メッシュの最初の方向を３Ｄパスの関連するセグメントの方向に変換する剛体変換を算出する段階、及び、かかる変換をセグメントに対応するメッシュの頂点に適用する段階と、
を有する。
【００２９】
真っ直ぐの円筒の開始位置がｚ座標０になるよう選択された場合、各円筒の点のｚ座標は、３Ｄパスの弦長パラメータ化と一致して使用され得る。パスの各位置ｐ（ｓ）に関し、局所的タンジェントの方向ｕ（ｓ）が決定される。続いて、ｓと同等のｚ座標を有する円筒の点は、ｐ（ｓ）から出た平面へとゆがめられ、ｕ（ｓ）に対して直交する。これらの基本的な操作は、図１に例示される。
【００３０】
しかしながら、いくつかの不適切な点が起こり得る。
【００３１】
まず、例えば、ポリラインのパスの２つの連続するセグメントの間に大きな角度があるために３Ｄパスが滑らかではない場合は、ゆがめられた円筒は自己交差し得、望ましくない自己交差するメッシュをもたらす。このように、主要な問題は、２つの連続するセグメントの間の方向が急激に変化したときの自己交差の潜在的な出現である。
【００３２】
続いて、ｚ軸をｕ（ｓ）の向きに一致させる剛体変換の選択がされたとき、これは結果もたらされるメッシュの望ましくない捩れにつながり得る。このように、他の問題は、変換中にメッシュの捩れとして見られ得る持続性制御の欠如である。
【００３３】
２）自己交差を制限又は防止すること：
自己交差の第１の原因は、ポリライン・パスの２つの連続するセグメントの間の方向ｕ（ｓ）の急速な変化である。局所的タンジェントｕ（ｓ）とｚ軸との間の参照符号θで示される角度が座標で示される場合、図２に示す階段型の階段状カーブが得られる。階段は、水平方向のセグメントＡ，Ｂ，．．．であり、図２中、実線で示される。かかる角度θは、図１に示される通り、各円の剛体変換を算出するよう直接使用されるため、階段の効果の低減は自己交差を妨げるのに役立つ。
【００３４】
図２によれば、階段の効果を低減させるよう、捩れを行うための算出された回転は、３Ｄパスのセグメント「Ａ」に対応し、一定であり、ｚ＝０とｚ＝１との間に位置付けられた全ての円は、同一の３Ｄ回転で回転される。同様に、セグメント「Ｂ」に対応する、即ち範囲［１．．２］内のｚ座標を有する全ての円は、同一の剛体変換を伴って動く。３Ｄ回転角度だけではなく、回転角度と同様に回転軸も算出されなければならず、回転軸はｓ座標に沿って区分的に一定であることが留意されるべきである。
【００３５】
２ａ）線形融合を使用することによって自己交差を制限する：
自己交差は、固有の変換が各セグメントに適用されない場合に防止され得る。その代わりに、変換は、２つの回転の間の線形補間を使用して２つの連続するセグメントの間に融合される。１つの３Ｄ剛体変換から他の３Ｄ剛体変換へのかかる線形補間は、四元数の形式主義を使用して容易に解決され得る。従って、本発明によれば、自己交差を制限するよう、方法は、各セグメント間の３Ｄ変換を線形に融合する段階、即ち、セグメントＡ，Ｂ等で形成された階段状カーブを、図２中点線で示される連続するカーブによって取り替える段階を有する。この線形融合は、四元数（Ｑ１及びＱ２で示す）として３Ｄ回転を考慮することによって達成される、（１）に示される形の四元法算法を伴う：
（１）
【００３６】
【数１】
0004545140.tif 000002
式（１）中、θは回転アングルであり、（ｕｘ，ｕｙ，ｕｚ）は回転軸である。融合自体は、真っ直ぐに、２つの四元数Ｑ１，Ｑ２の正規化された線形の組合せである：
（２）
【００３７】
【数２】
0004545140.tif 000003
非線形融合は、調査され得る。しかし、例えば管状器官分割では、線形融合は十分安定していることが証明された。
【００３８】
図３Ａ及び図３Ｂは、８個の点のみ含有する３Ｄパス上の回転融合の効果を示すものであり、１つのセグメントから他のセグメントへ非常に大きな角度変更を有する。図３Ａ中、３Ｄ回転融合は無く、１ａ，２ａ，３ａ等の結節点において異なる円が交差し、生成された単一メッシュはいくつかの自己交差を含む。図３Ｂ中、回転の線形融合によって、異なる円が一方向から次の方向へと平滑に変形され、結果的に、点１ｂ，２ｂ，３ｂに見られる通り、より定まったメッシュがもたらされる。
【００３９】
２ｂ）自動半径調整によって自己交差を防止する。
【００４０】
１つのセグメントから他のセグメントへの３Ｄ剛性変換の線形融合は、自己交差を避けるように常に十分であるわけではない。明らかに、かかる自己交差もまた、３Ｄパスの局所的曲率と作られたメッシュの所望される半径との間の回転に依存する。後者が曲率の局所的半径より大きい場合、曲率の半径が曲率に反比例することは判明しているため、曲率が高いときには後者は小さく、続いて自己交差が発生する。このように、ｓ座標に沿った剛体変換の平滑な進化が、上述された線形融合の操作によって保証される場合、なんらかの自己交差がまだ発生し得る。最初の円筒のｒで示される半径と、２つの連続する個別の円の距離と、３Ｄパスのｃで示される曲率との間にある関係は、かかる自己交差の作成に影響を及ぼす。非常に曲がったパスで大きな半径ｒを備えた円筒をゆがめようとすると、なんらかの深刻な問題が起きるであろう。それ故に、非常にカーブしたゾーンで円筒の直径を局所的に自動的に低減することが望ましい。本発明によれば、メッシュの半径は、曲率及び点のサンプル距離、及び所望される入力半径に基づき、自動的に適合される。３Ｄパスの曲率はいずれの点でも算出され得るという条件で、３Ｄ回転と組み合わされた縮小要因が計算され、相似変換につながる。本発明は器官に関わるため、与えられたポリラインは、単純な概算を使用するのに十分に平滑であることが予測される。３Ｄパスの曲率座標ｓ＝ｔで曲率ｃ（ｔ）が以下のと通り算出される。
【００４１】
（３）
【００４２】
【数３】
0004545140.tif 000004
式（３）中、パスの一次及び二次導関数は、（４）及び（５）の公式で推定される。
【００４３】
（４）
【００４４】
【数４】
0004545140.tif 000005
（５）
【００４５】
【数５】
0004545140.tif 000006
これは、各セグメントの長さがゼロに近く、それらが非常に均一であるときにのみ、有効である。かかる小さな規則性の問題を解決し、推定される曲率の大きな変化を避けるため、曲率値のガウス・フィルタが適用され得る。最後に、３Ｄパスの各点の算出された曲率値は、円筒の異なる円に対する縮小要因を適用するよう使用される。ｋで示されるかかる縮小要因は、最初の円筒の半径ｒ及び１／ｃと等しい３Ｄパスの曲率の推定半径に依存する：
（６）
【００４６】
【数６】
0004545140.tif 000007
状況に応じて、上述の回転融合と同様に、縮小要因はもまた融合され得る。半径調整の技術は、
（７）
【００４７】
【数７】
0004545140.tif 000008
に基づいて線形融合を使用し得、
（８）
【００４８】
【数８】
0004545140.tif 000009
であるか、又は、かかる半径調整の技術は、双三次スプライン補間を使用し得、ｕはｖによって交換されるため、
（９）
【００４９】
【数９】
0004545140.tif 000010
となる。
【００５０】
図５は、算出された曲率嬢の双三次保管の一例を図示する。曲線を成す座標ｓは、セグメントのインデックスに関連する。太線は、曲率１／ｃの半径のいくつかのサンプル値を表し、参照符号ｃは、ｓ＝０乃至ｓ＝５の位置で算出された曲率である。細い線で示されるカーブは、これらのまばらなデータの双三次補間を図示する。直線のものと比較してここで双三次補間を使用する利点は、機能が、連続する導関数（結節点でのタンジェントは水平である）を有することであるため、円筒の半径の変化はよりゆっくりと行なわれる。
【００５１】
これらの予防措置後もまだ自己交差がある場合は、従来技術に説明された変換に関連して冒頭部分で説明した通り、アクティブな輪郭アルゴリズムの内力を伴った自動的なメッシュ修繕及び平滑化が適用され得る。
【００５２】
図６Ａ及び図６Ｂは、自動的な半径調整の利点を例示する。図６Ａ及び図６Ｂは、結腸でのパスからのメッシュ作成の一例を図示する。図６Ａ中、メッシュモデルは半径調整を使用することなく作成され、結腸のメッシュは底部１０ａで自己交差し、３Ｄパスは非常に蛇行する。また、他の領域の屈曲部によって隠されている結腸のいくつかの領域を映像化するのは非常に困難である。図６Ｂ中、メッシュモデルは半径調整アルゴリズムを使用して作成され、自己交差は１０ｂで示される通り非常に低減される。しかしながら、結腸の一般的な形は、制限された半径の領域で摂動されない。１１ａ等の他の部分では、図６Ａと図６Ｂとの間で半径は変化しない。制限された半径の領域では、結腸の異なる領域の映像化及び追跡は、非常に向上する。
【００５３】
３）最低限の捩れを課す：
メッシュの捩れは、２つの連続する回転、即ち剛体変換の距離が最低限であるときに最低限化される。３Ｄ回転は、ｚ軸である最初のメッシュ方向から目標であるセグメントｕ（ｓ）までの最低限の回転として算出される。この回転は、軸Ｒ＝ｚ＾ｕ（ｓ）及び角度θ＝（ｚ，ｕ（ｓ））を伴って定義される。本発明によれば、軸Ｒ及び角度θによって構成されるパラメータは、１つのセグメントから他のセグメントへ反復して算出される。従って、かかる方法は、以前のセグメントｕ（ｓ）に見つけられた回転と以前のセグメントｕ（ｓ）から最新のセグメントｕ（ｓ＋１）への最低限の回転との複合として、最新のセグメントｕ（ｓ＋１）の新しい回転を算出する段階を有する。実際、ｚ軸を３Ｄパスの局所的タンジェントであるｕ（ｓ）方向にマッピングする固有な３Ｄ剛体変換は無い。並進が必要な場合は、かかる併進を算出するのは容易である。しかし、真の問題は、
【００５４】
【数１０】
0004545140.tif 000011
ベクトル（０，０，１）を
【００５５】
【数１１】
0004545140.tif 000012
へと変換する３Ｄ回転を見つけることである。最も単純な回転は、
【００５６】
【数１２】
0004545140.tif 000013
まで低減される軸の周囲での回転であり、以下の式、
（１０）
【００５７】
【数１３】
0004545140.tif 000014
で表される。回転角度θは、以下の式、
（１１）
【００５８】
【数１４】
0004545140.tif 000015
で表される。単純性を要求され得るかかる回転は、以下の式、
（１２）
【００５９】
【数１５】
0004545140.tif 000016
で表され、
【００６０】
【数１６】
0004545140.tif 000017
ベクトルを
【００６１】
【数１７】
0004545140.tif 000018
方向へマッピングする。しかし、
【００６２】
【数１８】
0004545140.tif 000019
に追随されるｚ軸周囲の回転（ｚ軸には触れないまま）で構成された他の回転も同様にマッピングする。それにもかかわらず、１つの回転のみが、最低限の変位の条件を充たし、正確に
【００６３】
【数１９】
0004545140.tif 000020
である。問題は、
【００６４】
【数２０】
0004545140.tif 000021
と
【００６５】
【数２１】
0004545140.tif 000022
との間では、
【００６６】
【数２２】
0004545140.tif 000023
から
【００６７】
【数２３】
0004545140.tif 000024
への各段階の間よりも、最低限の根実施例が見つからないという事実に留意しなければならないことである。従って算出されるべき回転は、
【００６８】
【数２４】
0004545140.tif 000025
ほど容易ではないが、以下の式、
（１３）
【００６９】
【数２５】
0004545140.tif 000026
で表される付加的な複合に対応し、式（１３）中、
【００７０】
【数２６】
0004545140.tif 000027
は、
【００７１】
【数２７】
0004545140.tif 000028
から
【００７２】
【数２８】
0004545140.tif 000029
への最低限の回転のマッピングを指定する。
【００７３】
図７Ａ及び図７Ｂは、最低限の捩れを課することを例示する。図７Ａは、ｚ軸とｕ（ｓ）との間の最低限の回転のみを使用してメッシュを作成する一例を図示する。図７Ｂは、ｚ軸とｕ（ｓ）との間の最低限の回転、続いて、最低限の捩れに繋がる付加的な回転Ｒ（ｓ）を使用するメッシュの作成の一例を図示する。図７Ａ中、領域４ａ及び５ａ等で、セルが結節点でよじれているためにメッシュ上に捩れが表れているのが判る。その代わりに、図７Ｂ中、セルはメッシュ全体に亘って整列されている。特に、セルは図７Ａの４ａ及び５ａに対応する領域４ｂ及び５ｂでよく整列される。
【００７４】
蛇行性合成物質の例では、最低限の捩れアルゴリズムがよい結果を出していることも観察され得る。図８Ａは、最初の管状のメッシュモデルを示し、そこから蛇行整合性物質の例が作成されて図８Ｂに表示され、図８Ｃでは異なる角度から図示される。図８Ｄは、同一の例を３Ｄパスのよりよい分解度を伴って示す。図８Ｂ及び図８Ｃが３Ｄパスに５０のセグメントのみを使用したのに対し、図８Ｄ中は２００セグメントが３Ｄパスに使用された。図８Ａ乃至図８Ｄ中、直線は、最初の円筒１０上に白で示され、続く対応する線１１，１２，１３での最終的なメッシュ上での捩れを算定することができる。
【００７５】
４）  所望される最初の状態
上述された方法は、異なる種類の３Ｄパスと連動する。しかしながら、最善の結果は、鋭角が存在しないときに観察される。従って、当業者に既知であるなんらかの平滑化技術を使用して入力３Ｄパスを予め平滑化したほうがよい。更によりよい結果は、パスのセグメントの長さが一様であるときに得られる。
【００７６】
５）  医療用観察システム及び装置
上述された段階は、本発明の観察システムによって実行され得る。図９は、医療用検査装置に併合された、本発明に従った画像観察システムの一実施例の基本的な構成要素を図示する。医療用検査装置１５１は、患者が横たわるベッド１１０、又は、画像装置に対して患者の位置を特定する他の要素を有し得る。医療用画像装置１５１は、ＣＴスキャナ、又は、Ｘ線若しくは超音波装置等の他の医療用画像装置であり得る。装置１５１によって作られたデータは、方法の段階を実行する汎用コンピュータ等のデータ処理手段１５３に与えられる。データ処理手段１５３は、通例、モニタ１５４等の映像化装置、キーボード等の入力装置１５５又はマウス１５６、位置決め装置等に関連付けられ、ユーザによって操作され得、システムと情報をやりとりする。データ処理装置１５３は、本発明に従って医療用画像データを処理する方法を実施するようプログラムされる。特には、データ処理装置１５３は、本方法の段階を実行するよう算出手段と記憶手段とを有する。方法を実行するよう予めプログラムされた命令を有するコンピュータプログラムが、実行され得る。
【００７７】
図面及び前出の説明は、本発明を制限するよりも例示するものである。添付の請求項の範囲内に入る代替案は多数あることは明らかである。更に、本発明は、表示用の画像データを生成することに関して説明されているが、本発明は、ディスプレイ装置上の表示、及び印刷を制限的にではなく含んだ画像データの映像化の略いかなる形式も含む。請求構内のいかなる参照符号も、請求項を制限するものとして解釈されるべきではない。
【図面の簡単な説明】
【００７８】
【図１】指定された点の所定のパスに基き、メッシュ融合の段階をセグメント毎に例示する。
【図２】最初の管状モデルに関連のｚ軸に対して３Ｄパスへの局所的なタンジェントｕ（ｓ）の角度の偏差を例示する。
【図３Ａ】線形変換融合無しのメッシュの作成と、線形変換融合有りのメッシュの作成を夫々環状に例示する。
【図３Ｂ】線形変換融合無しのメッシュの作成と、線形変換融合有りのメッシュの作成を夫々環状に例示する。
【図４Ａ】線形変換融合無しのメッシュの作成と、線形変換融合有りのメッシュの作成を夫々単一メッシュで例示する。
【図４Ｂ】線形変換融合無しのメッシュの作成と、線形変換融合有りのメッシュの作成を夫々単一メッシュで例示する。
【図５】算出された曲率での双三次スプライン補間の例を図示する。
【図６Ａ】線形変換融合及び捩れ最低限化を伴うメッシュモデル作成の例を、結腸の中心線のパスから図示する。
【図６Ｂ】結腸の中心線のパスに続くメッシュモデルに関連し、半径調整の段階によって更に修正された図６Ａの例を図示する。
【図７Ａ】付加的な回転Ｒ（ｓ）無し及び有りでの、局所的なタンジェントｕ（ｓ）とｚ軸との間の最低限の回転を使用する、メッシュの作成を例示する。
【図７Ｂ】付加的な回転Ｒ（ｓ）無し及び有りでの、局所的なタンジェントｕ（ｓ）とｚ軸との間の最低限の回転を使用する、メッシュの作成を例示する。
【図８Ａ】最初の管状メッシュモデルを図示する。
【図８Ｂ】捩れ最低限化につながる、線形変換融合及び付加的な回転Ｒ（ｓ）を伴うメッシュ作成の合成の蛇行の例を図示する。
【図８Ｃ】図８Ｂと同一の例を、他の方向で図示する。
【図８Ｄ】図６Ｂと同一の例を、３Ｄパスをより高い分解度で図示する。
【図９】画像処理方法を実行する手段を有する、医療用観察システムを例示する。
【図１】
0004545140.tif 000030
【図２】
0004545140.tif 000031
【図３Ａ】
0004545140.tif 000032
【図３Ｂ】
0004545140.tif 000033
【図４Ａ】
0004545140.tif 000034
【図４Ｂ】
0004545140.tif 000035
【図５】
0004545140.tif 000036
【図６Ａ】
0004545140.tif 000037
【図６Ｂ】
0004545140.tif 000038
【図７Ａ】
0004545140.tif 000039
【図７Ｂ】
0004545140.tif 000040
【図８Ａ】
0004545140.tif 000041
【図８Ｂ】
0004545140.tif 000042
【図８Ｃ】
0004545140.tif 000043
【図８Ｄ】
0004545140.tif 000044
【図９】
0004545140.tif 000045