JP5537026B2 - マルチサーフェスモデリング - Google Patents

マルチサーフェスモデリング Download PDF

Info

Publication number
JP5537026B2
JP5537026B2 JP2008508386A JP2008508386A JP5537026B2 JP 5537026 B2 JP5537026 B2 JP 5537026B2 JP 2008508386 A JP2008508386 A JP 2008508386A JP 2008508386 A JP2008508386 A JP 2008508386A JP 5537026 B2 JP5537026 B2 JP 5537026B2
Authority
JP
Japan
Prior art keywords
surface mesh
mesh
interest
constructing
model
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.)
Active
Application number
JP2008508386A
Other languages
English (en)
Other versions
JP2008538970A (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 JP2008538970A publication Critical patent/JP2008538970A/ja
Application granted granted Critical
Publication of JP5537026B2 publication Critical patent/JP5537026B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation

Description

発明の詳細な説明
本発明はサーフェスモデリング(surface modeling)の分野に関する。具体的には、本発明は、関心対象の多次元データセットからのマルチサーフェスモデルの構成方法と、画像処理装置と、検査装置と、コンピュータ読み取り可能媒体と、プログラム要素とに関する。
臓器表面のモデリングは、コンピュータ支援教育で使用される方法であるが、画像ベースの医療診断や臨床的介入でも使用される方法である。この場合、モデルは、患者を診察する前に、またはそれに加えて、解剖学的情報を提供し得るものである。変形可能なモデルを提供してそれを画像に適合させることにより、臓器の検出、セグメンテーション、及びラベリングにおいて、モデルを使用する。最初のモデルはサンプル画像から事前に構成される。モデル条件(例えば所定トポロジ(topology)を有する所定数の表面三角形)を決定し、モデルジオメトリを各トレーニング画像に適応させる。このように、モデル条件(すなわち、頂点または三角形の位置)により異なる画像中の解剖学的エンティティ間の対応を技術的に確立する。
かかるモデリングには、2次元の(ほとんど閉じた)多様体を使用してもよい。また、複数のかかるモデルを使用してもよい。この場合モデル間の距離もモデル化してもよい。
しかし、心臓領域では、キャビティ(cavities)(心房、心室、及びその入出流)または筋肉の一方だけをモデル化しても十分ではない。両方とも重要であり、両方の形状が相互に依存しているからである。
サーフェスモデリング(surface modeling)を改良することが望ましい。
本発明によると、関心対象の多次元データセットからマルチサーフェスモデルを構成する方法が提供される。該方法は、前記関心対象の第1の部分に対応する適応的第1のサーフェスメッシュを構成する段階と、前記関心対象の第2の部分に対応する適応的第2のサーフェスメッシュを構成する段階と、前記第1のサーフェスメッシュと前記第2のサーフェスメッシュとを結合して、前記関心対象のマルチサーフェスモデルを求める段階とを有する。
このように、本発明によると、前記マルチサーフェスモデルの構成を自動化することが
できる。単連結(single)の臓器部分のベーシック(basic)サーフェスメッシュを構成
してもよい。これらのメッシュを結合して、それにより関心対象のマルチサーフェスモデ
ルを形成することができる。このマルチサーフェスモデルは2次元でない多様体メッシュであってもよい。
本発明の実施形態は従属項に開示されている。
本発明の他の一実施形態によると、前記第1のサーフェスメッシュは第1の2次元多様体であり、前記第2のサーフェスメッシュは第2の2次元多様体であり、前記第1のサーフェスメッシュを構成する段階は、前記第1のサーフェスメッシュを前記関心対象の前記第1の部分に適合させる段階を有し、前記第2のサーフェスメッシュを構成する段階は、前記第2のサーフェスメッシュを前記関心対象の前記第1の部分に適合させる段階を有する。
こうすることにより、心室、心房、付属血管等の異なる解剖学的構成要素(anatomical entities)のモデリングを自動化できる。
本発明の他の一実施形態によると、前記第1のサーフェスメッシュと前記第2のサーフェスメッシュを結合する段階は、第1の演算、第2の演算、及び第3の演算を含み、前記第1、第2、及び第3の演算は2項演算(binary)である。さらに、前記第1の演算は和集合演算であり、前記第2の演算は共通部分演算であり、前記第3の演算は差集合演算である。
第1のサーフェスメッシュと前記第2のサーフェスメッシュの結合の際に3つの演算を適用することにより、マルチサーフェスモデルを効率的かつ正確に構成することができる。
本発明の他の一実施形態によると、前記第1のサーフェスメッシュは第1のメッシュプリミティブを有し、前記第2のサーフェスメッシュは第2のメッシュプリミティブを有する。さらに、前記第1のメッシュプリミティブは前記第2のメッシュプリミティブに対応する。前記第1のサーフェスメッシュと前記第2のサーフェスメッシュを結合する段階は第4の演算を含み、その第4の演算は前記第1のメッシュプリミティブと前記第2のメッシュプリミティブのマージである。
マージ又はジョイン(merging or join)演算を行うことにより、両方のメッシュにある対応する面(faces)が一体化できる。
本発明の他の一実施形態によると、該方法は、前記第1のサーフェスメッシュの第3のメッシュプリミティブに対応するラベルを生成する段階をさらに有し、前記ラベルは前記第3のメッシュプリミティブに関する情報を含む。
再構成段階において、マルチサーフェスメッシュの各メッシュプリミティブには、それがどのボリュームを分けるか、またはどのサーフェスに属するかを決定する対応ラベルが設けられる。
本発明の他の実施形態によると、本方法は、前記関心対象の結果として得られるマルチサーフェスモデルの精度が所定値より低いとき、前記第1のサーフェスメッシュと前記第2のサーフェスメッシュの少なくとも1つを細かくする段階と、前記関心対象について得られるマルチサーフェスモデルの精度が所定値より高いとき、前記第1のサーフェスメッシュと前記第2のサーフェスメッシュの少なくとも1つを粗くする段階とのうち少なくとも一方をさらに有する。
それゆえ、共通部分領域であってエッジ三角形またはメッシュプリミティブが何回か分割されているところでは、局所的にオーバーサンプリングとなっているので、これを補正することができる。
本発明の他の一実施形態によると、前記第1のサーフェスメッシュを構成する段階と、前記第2のサーフェスメッシュを構成する段階と、前記第1のサーフェスメッシュと前記第2のサーフェスメッシュを結合する段階はインターラクティブにユーザに誘導されて行われる。
こうすることにより、例えば非常にフレキシブルになる。例えば、ユーザは関心のある部分をインターラクティブに選択したり、ズームしたりできる。様々なユーザによるインターラクション(interactions)の間に、サーフェス表示やトレーニング画像(training image)へのサーフェスメッシュの適応が行われてもよい。
本発明の他の一実施形態によると、前記第1のサーフェスメッシュを構成する段階と、前記第2のサーフェスメッシュを構成する段階と、前記第1のサーフェスメッシュと前記第2のサーフェスメッシュを一体化する段階はバッチモードで行われ、前記第1の演算と第2の演算と第3の演算と第4の演算よりなる群から選択された少なくとも1つの演算のシーケンスを行う。
このように、全体的な構成はバッチモードで実行してもよい。この場合、必要な演算のシーケンスはプログラミング言語等によりコード化される。こうすることにより自動化のレベルを高くすることができる。
本発明の他の一実施形態によると、関心対象の多次元データセットからマルチサーフェスモデルを構成する方法を提供する画像処理装置が提供される。画像処理装置は、多次元データセットを記憶するメモリと、計算ユニットとを有する。計算ユニットは、本発明の一実施形態による方法を実行するように構成されている。
本発明の他の一実施形態によると、関心対象の多次元データセットからマルチサーフェスモデルを構成する検査装置が提供される。該装置は、上記の方法ステップを実行するように構成された計算ユニットを有する。
本発明の他の一実施形態によると、検査装置は、超音波画像化システム、CT(コンピュータ断層撮影)画像化システム、CSCT(コヒーレントスキャタコンピュータトモグラフィ)画像化システム、PET(陽電子放出トモグラフィ)画像化システム、SPECT(単一電子放出コンピュータトモグラフィ)画像化システム、超音波、及びMR(磁気共鳴)画像化システムよりなる群から選択することができる。それゆえ、複数の異なる検査(例えば、悪性と良性の病変)をできる検査ツールを提供することができる。
本発明の他の一実施形態によると、本検査装置は、さらに、前記関心対象に電磁放射を発するように構成された電磁放射源と、前記電磁放射源と検出器との間に配置されたコリメータとを有し、前記コリメータは前記電磁放射源により発せられた電磁ビームを平行化してファンビームまたはコーンビームを形成するように構成される。
さらに、本発明の他の一実施形態によると、前記検査装置は、手荷物検査装置、医療アプリケーション装置、材料テスト装置、または材料科学分析装置に適用され得る。本発明の応用分野は手荷物検査である。本発明の機能性により、手荷物の内容の安全かつ信頼性の高い分析が可能となり、不審物の検出やかかる手荷物の中の材料の種類の決定も可能となる。
本発明の他の一実施形態によると、コンピュータ読み取り可能媒体が提供される。前記媒体には検査装置により関心対象の多次元データセットからマルチサーフェスモデルを構成するコンピュータプログラムが格納され、プロセッサにより実行されたとき、上記の方法ステップを実行するように構成される。
また、本発明は、関心対象の多次元データセットからマルチサーフェスモデルを構成するプログラム要素にも関する。このプログラム要素は、コンピュータ読み取り可能媒体に格納することもできる。プログラム要素は、次のステップを実行するように構成され得る。前記関心対象の第1の部分に対応する適応的第1のサーフェスメッシュを構成する段階と、前記関心対象の第2の部分に対応する適応的第2のサーフェスメッシュを構成する段階と、前記第1のサーフェスメッシュと前記第2のサーフェスメッシュとを結合して、前記関心対象のマルチサーフェスモデルを求める段階。
本発明の一実施形態によるプログラム要素は、好ましくはデータプロセッサの作業記憶領域にロードされる。このように、データプロセッサは本発明の方法の実施形態例を実行するように構成され得る。コンピュータプログラムは、C++等のどの好適なプログラミング言語で記述してもよく、CD−ROM等のコンピュータ読み取り可能媒体に格納してもよい。また、コンピュータプログラムは、ワールドワイドウェブ等のネットワークから入手することもできる。本コンピュータプログラムをネットワークから画像処理ユニットまたはプロセッサ、またはその他の好適なコンピュータにダウンロードしてもよい。
本発明の一実施形態の要点は、単連結かつベーシック(basic)な2次元多様体のサーフェスメッシュを結合して、マルチサーフェスモデルを形成することである。このマルチサーフェスモデルは2次元でない多様体メッシュでなくなることもある。結合(combination)には、両方の曲面集合を一体化して両方のメッシュにある対応する曲面を一体化するマージ又はジョイン(merging or join)演算も含む。
本発明の上記その他の態様を、以下に説明する実施形態を参照して明らかにし、詳しく説明する。
添付図面を参照して、本発明の実施形態を以下に説明する。
図面中の図は概略である。別の図面においても、同様または同一の要素には同じ参照数字を付した。
図1は、本発明の一実施形態によるCTスキャナシステムの一実施形態例を示す図である。この実施形態例を参照して、本発明を医療画像化に応用する場合について説明する。しかし、留意すべきことは、本発明はこのアプリケーションに限定はされず、手荷物検査やその他の産業用途(例えば、材料検査等)にも適用できることである。
図1に示したコンピュータトモグラフィ装置100はコーンビームCTスキャナである。しかし、本発明はファンビームジオメトリ(geometry)でも実行できる。図1に示したCTスキャナはガントリ101を有する。このガントリは回転軸102の回りに回転可能である。ガントリ101はモータ103により駆動される。参照数字104は、X線光源等の放射源を示し、本発明の一態様では、多色放射光(polychromatic radiation)を放射する。
参照数字105はアパチャー系(aperture system)を示す。このアパチャー系は、放射光源から放射された放射光ビームをコーン形状の放射光ビーム106に形成する。コーンビーム106の方向は、ガントリ101の中心すなわちCTスキャナの検査領域に配置された関心対象107を透過して、ディテクタ108に当たるようになっている。図1から分かるように、ディテクタ108は、ガントリ101上に放射光源104に対向して配置され、ディテクタ108の表面はコーンビーム106で覆われる。ディテクタ108は、図1に示したように、複数のディテクタ要素123を有する。各ディテクタ要素は、関心対象107を透過したX線または個々の光子をエネルギー分解して検出できる。
関心対象107をスキャンする際、放射光源104、アパチャー系105、及びディテクタ108は、矢印116で示される方向にガントリ101と共に回転する。放射光源104、アパチャー系105、及びディテクタ108を有するガントリ101を回転させるため、モータ103がモータ制御部117に接続されている。このモータ制御部は計算・決定ユニット118に接続されている。
図1では、関心対象107はコンベヤベルト119上に配置されている患者または手荷物である。関心対象107のスキャンの際、ガントリ101は手荷物107の回りを回転するが、コンベヤベルト119は関心対象107をガントリ101の回転軸102と平行な方向に並進移動する。これにより、関心対象107はヘリカルスキャンパスに沿ってスキャンされる。コンベヤベルト119は、スキャン中に停止して、単一スライスを測定することもある。コンベヤベルト119を設ける替わりに、例えば、医療への応用では、関心対象107が患者であり、可動テーブルを使用する。しかし、留意すべきことは、説明する全ての場合において、円状スキャンを実行することが可能である。この場合、回転軸102に平行な方向での並進移動はなく、回転軸102の回りをガントリ101が回転するだけである。
さらに、強調しておくが、図1に示したコーンビーム(cone-beam)構成に替えて、本発明はファンビーム(fan-beam)構成により実施することもできる。主ファンビームを発生するために、アパチャー系105をスリットコリメータとして構成し得る。
ディテクタ108は計算ユニット部118に接続されている。計算ユニット118は、検出結果、すなわちディテクタ108のディテクタ要素123からのリードアウト(read-outs)を受け取り、そのリードアウトに基づいてスキャン結果を決定する。さらに、計算ユニット118は、モータ制御部117と通信する。これは、モータ103、120によるガントリ101の動きをコンベヤベルト119と調和させるためである。
本発明の一実施形態では、計算ユニット118は、ディテクタ108のリードアウトから統計的方法を用いて画像を構成するように構成され得る。計算ユニット118により生成された再構成画像は、インターフェイス122を介してディスプレイ(図1には示さず)に出力される。
計算ユニット118は、データプロセッサにより実施され、ディテクタ108のディテクタ要素123からのリードアウト(read-outs)を処理する。
さらに、図1から分かるように、計算ユニット118はスピーカ121に接続され、例えば、手荷物107中に不審物を検出した場合、自動的に警報を発する。
コンピュータトモグラフィ装置100は、関心対象107の検査をするものであり、ディテクタ108を含んでいる。このディテクタはマトリックス状に配置された複数の検出要素123を有し、各検出要素123はX線を閾値ベースで検出し、よってエネルギー分解して検出するように構成されている。さらに、コンピュータトモグラフィ装置100は決定ユニットすなわち再構成ユニット118を有する。この決定ユニットすなわち再構成ユニットは、関心対象107の多次元データセットからマルチサーフェスモデルを構成することができる。
コンピュータトモグラフィ装置100は、関心対象107にX線を放射するように構成されたX線源104を有する。コリメータ105が電磁放射源104と検出要素123との間に設けられ、電磁放射源104から放射された電磁放射ビームを平行にして(collimate)、コーンビームを形成する。あるいは、図1には図示していないが、コリメータ105の代わりにスリットコリメータを使用してファンビームを作ることもできる。検出要素123は、マルチスライスのディテクタ配列108を形成する、コンピュータトモグラフィ装置100は医療用画像化装置として構成することもできる。
図2は心臓CT画像であり、心臓の画像スライスを示す図である。図2から分かるように、動脈201、右心房202、左心室(left blood pool)203、大動脈204の断面と、心筋205の断面とがはっきりと見える。それゆえ、モデルのトポロジ(topology)は多様な種類の近傍(manifold kinds of neighbourhoods)となる。留意すべきことは、画像はマルチサーフェスモデルをまだ見ぬ画像(unseen image)に適応し、その人体構造上のコンポーネントをセグメント化し、ラベル付けしたものであることである。
図3は、本発明による方法の一実施形態によるモデル構成プロセスの、異なる6つの段階を示す図である。段階(a)はベーシックメッシュ301を示し、このベーシックメッシュ(basic mesh)を関心対象(すなわち、心臓の左心房)の所定部分に適用可能であり、適応されたものであってもよい。段階(b)乃至(d)は中間メッシュを示す。この中間メッシュは本発明の一実施形態により適応または結合されたものである。段階(b)は左血管内貯留(left blood pool)302を示し、段階(c)は左筋303を示し、段階(d)は右血管内貯留(right blood pool)304を示している。段階(e)は心臓の最終的なマルチサーフェスモデルを示す。このモデルは、左血管内貯留302、左筋303、及び右血管内貯留304を有する。段階(f)は左心室305の内部のみを示す。この左心室は、本発明の他の実施形態では、段階(e)のマルチサーフェスモデルからサブモデルとして構成することもできる。
図4は、本発明の一実施形態による関心対象(すなわち心臓)の多次元データセットからマルチサーフェスモデルを構成する方法の一実施形態を示すフローチャートである。モデリングフレームワークは、マルチサーフェスモデルの構成を自動化することもできる。単連結の臓器部分のベーシック(basic)な2次元多様体であるベーシックサーフェスメッシュを構成してもよい。(図3の段階(a)に示したように)、球や管などの簡単な幾何学形状から初めて、サーフェスメッシュをトレーニング画像(training image)の目標となる臓器の対応部分に適応させることができる。
フレームワークの中心部分は、これらのベーシックメッシュを結合して、ベーシックメッシュからマルチサーフェスモデルを構成する。このマルチサーフェスモデルは2次元でない多様体であってもよい。関心対象のマルチサーフェスモデルは異なるボリューム(volumnes)を含んでいてもよく、これらのボリュームは結合的(共通面を共有する)であっても、非結合的であっても、重なっていてもよい。
再構成段階において、マルチサーフェスメッシュの各メッシュプリミティブ(mesh primitive)にはラベルが付される。メッシュプリミティブは例えば三角形である。ラベル付けは、各メッシュプリミティブがどのボリュームに分かれるか、または各メッシュプリミティブがどのサーフェスに属するかを決めるものである。これは、モデルの可視化のみでなく、画像への適応にとっても重要な特徴(feature)である。画像特徴が異なれば異なる臓器部分であると見なされる。
本発明の一実施形態では、本フレームワークは、血液で満たされているすべてのキャビティと心筋を含むマルチサーフェス三角形モデルを構成し、このモデルを未だ見ぬ画像(unseen images)に適応させる。
マルチサーフェスメッシュ構成段階では、4つの基本演算子を使用する。これらは2項(binary)演算子であり、すなわち2つのメッシュを1つに結合するものである。それらの演算子のうち3つ(和集合(union)演算子、積集合(intersection)演算子、差集合(difference)演算子)は、2つのメッシュに含まれるボリュームを考慮し、それに集合演算を適用し、その集合演算のサーフェスを与える。本発明の一態様によると、これらの演算子は、変形可能サーフェスモデルの体積画像(volumetric images)の場合に使用できる。
オペランドも演算結果もともに2次元多様体である。4番目の演算子(ジョイン又はマージ(join or merge))の結果は2次元でない多様体である。4番目の演算子は、三角形の集合や(三角形ではなく、例えば多角形の)メッシュプリミティブを一体化するが、両方のメッシュに含まれる対応する三角形を一体化してもよい。これらの演算子は、ベーシックサーフェスメッシュ(basic surface meshes)に正しく適用されると、最終的なマルチサーフェスモデルを構成するために十分である。また、結果のモデルの精度が所定の空間的精度より高いか低いとき、メッシュを細かくしたり粗くしたりしてもよい。このことは共通部分領域において特に重要である。この場合、エッジ(edge)の三角形やメッシュプリミティブは通常何回か分割されており、そのために局所的にオーバーサンプリングになっており、補正が必要である。
留意すべきこととして、本発明の一態様では、第1のサーフェスメッシュまたは第2のサーフェスメッシュは、(少なくとも1つの基本演算がすでに適用された)複数のベーシックサーフェスメッシュを含んでもよい。それゆえ、関心対象107のマルチサーフェスモデルを反復的または段階的に構成することができる。
本発明の一実施形態では、ユーザがインターラクティブに誘導してすぐにサーフェス表示とトレーニング画像への適応をさせて構成することができる。構成プロセスはバッチモードで実行してもよい。この場合、必要な演算のシーケンスはプログラミング言語のようにコード化される。
本発明の一態様では、モデルは心臓の左側と右側の血液プール(blood pool)をカバーする。心室の血液プールは心房の血液プールと区別してもよい。付属している血管(すなわち、大動脈、肺動脈、大静脈、主要肺静脈)もすべてモデル化してもよい。さらに、左心筋もモデルに表すことができる。
隣接する体積エンティティを含めるには、2次元多様体を越えるサーフェスモデリングスキームが必要である。本発明の一態様では、メッシュプリミティブは三角形の面である。4つ以上の近傍(neighbours)を有し、どの面でも複数のサーフェス(surfaces)がエッジ(edge)を共有する面(faces)がある。マルチスケールまたはマルチ解像度のアプローチを実現するため、または正確性と複雑性の理想的なトレードオフを見つけるため、サーフェス離散化をマルチ解像度で表示してもよい。
最初のステップS1において、単連結のベーシックな形状が構成される。このベーシックな形状は、例えば、心房となる球面、付属血管となる管、心室となる開いた楕円等である。それぞれの適応されるベーシック形状はサーフェスメッシュであり、このサーフェスメッシュは複数のメッシュプリミティブ(三角形等)を有し、2次元多様体である。それぞれのベーシック形状は人体構造上の構成要素(anatomical entity)をモデル化するものであってもよい。
第2のステップS2において、ベーシック形状(basic shapes)がトレーニング画像(training image)に配置され、対応する構成要素(entities)に適応される。さらに、所定レベルの精度を得るため、再サンプリング(re-sampling)をしてもよい。この再サンプリングには、サーフェスメッシュが所定の精度値より低い場合にサーフェスメッシュを細かくすることが含まれる。さらに、再サンプリングには、サーフェスメッシュの精度が所定の値より高い場合にサーフェスメッシュを粗くすることも含まれる。
第3のステップS3において、単連結ベーシック2次元多様体(single basic two-dimensional manifolds)が結合されて、マルチサーフェスモデルが形成される。結合ステップS3は、サーフェスメッシュに4つの基本演算を連続的に適用して行われる。最初はベーシックメッシュからスタートする。2つのメッシュの体積を考慮する体積集合演算(volumetric set operations)において、共通部分(∪)演算または差(\)演算を施し、その結果、サーフェスメッシュを作る。これらの演算はB×B→Bと定義される。ここで、Bは2次元多様体メッシュである。これらの演算におけるさらに別の制約として、両方のメッシュ間の共通部分を閉じた多角形としてもよい。これには、近傍(例えば、心房を有する心室)と完全に重なる開いたベーシックメッシュが必要である。
ジョイン又はマージ(join or merge)演算(◇)は、B×B→Mと定義してもよい。ここで、Mは2次元でない多様体メッシュである。ジョイン演算子は、両方の面の集合を一体化するだけでもよいし、両方のメッシュにある対応する面を一体化してもよい。
さらに、単項演算子c1(M→M)は、パラメータ1により与えられるより小さなすべてのエッジ(edges)を削除し、三角形のラベルを保存する。単項演算子は、体積集合演算により発生する補助三角形を置換するために使用される。左血液プールplBは次式で構成される。
Figure 0005537026
完全なマルチサーフェスモデルを構成するために、第1の中間メッシュa1とa2を次式で構成することができる。
Figure 0005537026
ここで、a1はちょうど左心筋を取り囲んでいる。完全な左心臓モデルは次式により構成される。
Figure 0005537026
右血液プールは次式で構成される。
Figure 0005537026
左右の部分を融合して次式としてもよい。
Figure 0005537026
図3を参照して、ベーシックメッシュ(a)、中間メッシュ(b−d)、最終メッシュ(e)について、
外1
Figure 0005537026
を同一の視点から表示してもよい。エッジのサイズは2.5mmから5mmの間に設定される。また、最終メッシュを構成するのに必要なもの以外のサブメッシュを構成してもよい。例えば、左心房を除いた左心室の血液プールは、次式で構成される。
Figure 0005537026
結果として得られる関心対象のマルチサーフェスモデルにおいて、各面には、それが属する解剖学的構造を示すラベルが割り当てられる。この情報は、面(face)が最初のどのベーシック形状であったかを記憶することにより求めることができる。ベーシックな形状は、左心房al、左心室心内膜(内部vi)、左心室心外膜(外部vo)、大動脈a、上大静脈vs、下大静脈vl、右心房ar、右心室vr、肺動脈ap(右ブランチのみ)、左心房に排出される肺静脈(v1、v2、v3、v4)である。
図5は、本発明による方法の一実施形態例を実行する、本発明による画像処理装置の一実施形態例を示す図である。図5に示した画像処理装置400は、中央処理装置(CPU)または画像プロセッサ401を含む。この中央処理装置または画像プロセッサは、患者や手荷物等の関心対象を示す画像を格納するメモリ402に接続されている。データプロセッサ401は、診断装置(CT装置等)の複数の入出力ネットワークに接続されている。データプロセッサ401は、さらに、データプロセッサ401で計算、または適応された情報または画像を表示する表示装置403(コンピュータモニタ等)に接続されている。オペレータまたはユーザは、キーボード404及び/またはその他の出力装置(図5には図示せず)を介してデータプロセッサ401とインターラクト(interact)する。さらに、バスシステム405を介して、画像処理及び制御プロセッサ401を例えば動きモニターに接続することも可能である。この動きモニターは関心対象の動きを監視する。例えば、患者の肺を画像化する場合、動きセンサーは例えば呼気センサーである。心臓を画像化する場合、動きセンサーは、例えば心臓ゲート法(method of cardiac gating)により心臓の動きを補償する基礎として、心電図を形成する心電計であってもよい。
本発明による関心対象の検査により、マルチサーフェスモデル、関心対象の多次元データセットを構成でき、複数メッシュにある対応するメッシュプリミティブ面をマージすることができる。関心対象の画像は、CTその他のスキャナシステムにより取得することができる。例えば、MRI(磁気共鳴画像化)スキャナシステム、PET(陽電子放射トモグラフィ)スキャナシステム、SPECT(単光子放出コンピュータ断層撮影)スキャナシステム、または超音波画像化システム等である。
本発明の実施形態は、CTスキャナコンソール、画像化ワークステーション、診断医療ワークステーション、またはPACS(画像アーカイブ・通信システム)ワークステーションに対するソフトウェアオプションとして販売することができる。アプリケーションは、例えば、心臓・血管スクリーニングや心臓電気生理学等である。
留意すべきことは、「有する」という用語は他の要素やステップを排除するものではなく、「1つの」という用語は複数の場合を排除するものではないことである。また、請求項に記載されたいくつかの手段またはユニットの機能を単一のプロセッサまたはシステムが満たすことができることである。また、異なる実施形態に関して説明した要素を組み合わせてもよい。
請求項中の参照符号は、その請求項の範囲を限定するものと解釈してはならないことにも留意すべきである。
本発明によるCTスキャナの一実施形態を簡略化して示す概略図である。 患者の心臓の心臓CT画像のスライスを示す図である。 本発明による方法の一実施形態によるモデル構成プロセスの異なる6つの段階を示す図である。 本発明による関心対象の多次元データセットからマルチサーフェスモデルを構成する方法の一実施形態を示すフローチャートである。 本発明による方法の一実施形態例を実行する、本発明による画像処理装置の一実施形態を示す図である。

Claims (15)

  1. 関心対象の多次元データセットからマルチサーフェスモデルを構成する方法であって、
    前記関心対象の第1の部分に対応する適応的な第1のサーフェスメッシュを構成する段階と、
    前記関心対象の第2の部分に対応する適応的な第2のサーフェスメッシュを構成する段階と、
    前記第1のサーフェスメッシュと前記第2のサーフェスメッシュとを、対応する面を一体化することにより結合して、前記関心対象のマルチサーフェスモデルを求める段階とを有する方法。
  2. 前記第1のサーフェスメッシュは第1の2次元多様体であり、前記第2のサーフェスメッシュは第2の2次元多様体であり、
    前記第1のサーフェスメッシュを構成する段階は、前記第1のサーフェスメッシュを前記関心対象の前記第1の部分に適合させる段階を有し、
    前記第2のサーフェスメッシュを構成する段階は、前記第2のサーフェスメッシュを前記関心対象の前記第2の部分に適合させる段階を有する、請求項1に記載の方法。
  3. 前記第1のサーフェスメッシュと前記第2のサーフェスメッシュを結合する段階は、第1の演算、第2の演算、及び第3の演算を含み、
    前記第1、第2、及び第3の演算は2項演算であり、
    前記第1の演算は和集合演算であり、
    前記第2の演算は共通部分演算であり、
    前記第3の演算は差集合演算である、請求項1に記載の方法。
  4. 前記第1のサーフェスメッシュは第1のメッシュプリミティブを有し、
    前記第2のサーフェスメッシュは第2のメッシュプリミティブを有し、
    前記第1のメッシュプリミティブは前記第2のメッシュプリミティブに対応し、
    前記第1のサーフェスメッシュと前記第2のサーフェスメッシュ、対応する面を一体化することにより結合する段階は、第4の演算を含み、
    前記第4の演算は前記第1のメッシュプリミティブと前記第2のメッシュプリミティブのマージである請求項1に記載の方法。
  5. 前記第1のサーフェスメッシュの第3のメッシュプリミティブに対応するラベルを生成する段階をさらに有し、
    前記ラベルは前記第3のメッシュプリミティブに関する情報を含む、請求項1に記載の方法。
  6. 前記関心対象について得られるマルチサーフェスモデルの精度が所定値より低いとき、前記第1のサーフェスメッシュと前記第2のサーフェスメッシュの少なくとも1つを細かくする段階と、
    前記関心対象について得られるマルチサーフェスモデルの精度が所定値より高いとき、前記第1のサーフェスメッシュと前記第2のサーフェスメッシュの少なくとも1つを粗くする段階とのうち少なくとも一方をさらに有する、請求項1に記載の方法。
  7. 前記第1のサーフェスメッシュを構成する段階と、前記第2のサーフェスメッシュを構成する段階と、前記第1のサーフェスメッシュと前記第2のサーフェスメッシュを結合する段階はインターラクティブにユーザに誘導されて行われる、請求項1に記載の方法。
  8. 前記第1のサーフェスメッシュを構成する段階と、前記第2のサーフェスメッシュを構成する段階と、前記第1のサーフェスメッシュと前記第2のサーフェスメッシュを一体化する段階はバッチモードで行われる、
    前記第1の演算と第2の演算と第3の演算と第4の演算よりなる群から選択された少なくとも1つの演算のシーケンスを行う、請求項3に記載の方法。
  9. 関心対象の多次元データセットからマルチサーフェスモデルを構成する画像処理装置であって、
    多次元データセットを記憶するメモリと、
    前記関心対象の第1の部分に対応する適応的な第1のサーフェスメッシュを構成し、
    前記関心対象の第2の部分に対応する適応的な第2のサーフェスメッシュを構成し、
    前記第1のサーフェスメッシュと前記第2のサーフェスメッシュとを、対応する面を一体化することにより結合して、前記関心対象のマルチサーフェスモデルを求めるように構成された計算ユニットとを有する装置。
  10. 関心対象の多次元データセットからマルチサーフェスモデルを構成する検査装置であって、
    前記関心対象の第1の部分に対応する適応的な第1のサーフェスメッシュを構成し、
    前記関心対象の第2の部分に対応する適応的な第2のサーフェスメッシュを構成し、
    前記第1のサーフェスメッシュと前記第2のサーフェスメッシュとを、対応する面を一体化することにより結合して、前記関心対象のマルチサーフェスモデルを求めるように構成された計算ユニットとを有する装置。
  11. 前記検査装置は、超音波画像化システム、CT画像化システム、CSCT画像化システム、PET画像化システム、SPECT画像化システム、及びMR画像化システムよりなる群から選択される、請求項10に記載の検査装置。
  12. 前記関心対象に電磁放射を発するように構成された電磁放射源と、
    前記電磁放射源と検出器との間に配置されたコリメータとを有し、
    前記コリメータは前記電磁放射源により発せられた電磁ビームを平行化してファンビームまたはコーンビームを形成するように構成された、請求項10に記載の検査装置。
  13. 手荷物検査装置、医療アプリケーション装置、材料テスト装置、及び材料科学分析装置よりなる群のうちの1つとして構成されている、請求項10に記載の検査装置。
  14. コンピュータ読み取り可能媒体であって、前記媒体には検査装置により関心対象の多次元データセットからマルチサーフェスモデルを構成するコンピュータプログラムが格納され、プロセッサにより実行されたとき、
    前記関心対象の第1の部分に対応する適応的な第1のサーフェスメッシュを構成する段階と、
    前記関心対象の第2の部分に対応する適応的な第2のサーフェスメッシュを構成する段階と、
    前記第1のサーフェスメッシュと前記第2のサーフェスメッシュとを、対応する面を一体化することにより結合して、前記関心対象のマルチサーフェスモデルを求める段階とを実行するように構成されたコンピュータ読み取り可能媒体。
  15. 関心対象の多次元データセットからマルチサーフェスモデルを構成するプログラム要素であって、プロセッサにより実行されたとき、
    前記関心対象の第1の部分に対応する適応的な第1のサーフェスメッシュを構成する段階と、
    前記関心対象の第2の部分に対応する適応的な第2のサーフェスメッシュを構成する段階と、
    前記第1のサーフェスメッシュと前記第2のサーフェスメッシュとを、対応する面を一体化することにより結合して、前記関心対象のマルチサーフェスモデルを求める段階とを実行するように構成されたプログラム要素。
JP2008508386A 2005-04-29 2006-04-25 マルチサーフェスモデリング Active JP5537026B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP05103591.3 2005-04-29
EP05103591 2005-04-29
PCT/IB2006/051284 WO2006117719A2 (en) 2005-04-29 2006-04-25 Multi-surface modelling

Publications (2)

Publication Number Publication Date
JP2008538970A JP2008538970A (ja) 2008-11-13
JP5537026B2 true JP5537026B2 (ja) 2014-07-02

Family

ID=37308374

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008508386A Active JP5537026B2 (ja) 2005-04-29 2006-04-25 マルチサーフェスモデリング

Country Status (5)

Country Link
US (1) US20090267941A1 (ja)
EP (1) EP1877984B1 (ja)
JP (1) JP5537026B2 (ja)
CN (1) CN101167105B (ja)
WO (1) WO2006117719A2 (ja)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9405996B2 (en) * 2008-09-18 2016-08-02 Siemens Aktiengesellschaft Method and system for generating a personalized anatomical heart model
JP6001280B2 (ja) * 2012-03-09 2016-10-05 学校法人金沢医科大学 心臓矯正ネットの製造方法
TWI472792B (zh) * 2012-04-27 2015-02-11 Iner Aec Executive Yuan 輻射偵測信號處理方法與系統
JP5700712B2 (ja) * 2013-03-26 2015-04-15 国立大学法人愛媛大学 3次元ダイナミック心筋核医学画像データのパトラック解析
CN105814605B (zh) * 2013-12-10 2019-09-24 皇家飞利浦有限公司 对解剖结构的基于模型的分割
JP6619431B2 (ja) * 2014-11-28 2019-12-11 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 解剖学的構造のモデルベースセグメンテーション
US11478212B2 (en) * 2017-02-16 2022-10-25 Siemens Healthcare Gmbh Method for controlling scanner by estimating patient internal anatomical structures from surface data using body-surface and organ-surface latent variables
US10388015B2 (en) 2017-09-06 2019-08-20 International Business Machines Corporation Automated septal defect detection in cardiac computed tomography images
CN110327033B (zh) * 2019-04-04 2022-05-03 浙江工业大学 一种基于深度神经网络的心肌梗死心电图的筛查方法
US11831931B2 (en) * 2021-04-14 2023-11-28 Microsoft Technology Licensing, Llc Systems and methods for generating high-resolution video or animated surface meshes from low-resolution images
US11849220B2 (en) 2021-04-14 2023-12-19 Microsoft Technology Licensing, Llc Systems and methods for generating depth information from low-resolution images

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6184897B1 (en) * 1997-01-15 2001-02-06 International Business Machines Corporation Compressed representation of changing meshes and method to decompress
US6262737B1 (en) * 1998-01-30 2001-07-17 University Of Southern California 3D mesh compression and coding
US6879324B1 (en) * 1998-07-14 2005-04-12 Microsoft Corporation Regional progressive meshes
JP3844174B2 (ja) * 1998-10-26 2006-11-08 マツダ株式会社 Cadデータ処理装置及び記憶媒体
EP1264281A4 (en) * 2000-02-25 2007-07-11 Univ New York State Res Found ARRANGEMENT AND METHOD FOR PROCESSING AND PLAYING A VOLUME
EP1221674A3 (en) 2001-01-05 2003-09-24 Interuniversitair Microelektronica Centrum Vzw System and method to obtain surface structures of multidimensional objects, and to represent those surface structures for animation, transmission and display
JP2002328958A (ja) * 2001-04-27 2002-11-15 Fuji Heavy Ind Ltd メッシュ生成方法、プログラム、記録媒体およびメッシュ生成システム
US7324116B2 (en) * 2002-06-20 2008-01-29 Microsoft Corporation Systems and methods for providing controllable texture sampling
JP4575157B2 (ja) * 2002-07-19 2010-11-04 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ メッシュ適応による多数の対象又は複合的な対象の同時セグメンテーション

Also Published As

Publication number Publication date
EP1877984B1 (en) 2014-05-07
EP1877984A2 (en) 2008-01-16
WO2006117719A3 (en) 2007-03-15
JP2008538970A (ja) 2008-11-13
CN101167105B (zh) 2013-03-27
US20090267941A1 (en) 2009-10-29
WO2006117719A2 (en) 2006-11-09
CN101167105A (zh) 2008-04-23

Similar Documents

Publication Publication Date Title
JP5537026B2 (ja) マルチサーフェスモデリング
US11501485B2 (en) System and method for image-based object modeling using multiple image acquisitions or reconstructions
US9576391B2 (en) Tomography apparatus and method of reconstructing a tomography image by the tomography apparatus
US9514530B2 (en) Systems and methods for image-based object modeling using multiple image acquisitions or reconstructions
JP7139479B2 (ja) 画像処理装置、医用画像診断装置、及びプログラム
US10937209B2 (en) Tomography imaging apparatus and method of reconstructing tomography image
JP5345934B2 (ja) 観察のための三次元レンダリングからのデータ集合の選択
EP3107457B1 (en) Tomography apparatus and method of reconstructing a tomography image by the tomography apparatus
JP5919296B2 (ja) イメージレジストレーション装置
US10275946B2 (en) Visualization of imaging uncertainty
CN108876794A (zh) 体积图像数据中的动脉瘤与载瘤血管的隔离
KR101775556B1 (ko) 단층 촬영 장치 및 그에 따른 단층 영상 처리 방법
CN111311551B (zh) 基于机器学习的生理运动测量系统和方法
JP2011506032A (ja) 一貫性指標に基づく画像レジストレーション
EP2220618B1 (en) Apparatus for determining a parameter of a moving object
Medina et al. Left ventricle myocardium segmentation in multi-slice computerized tomography
JP7318058B2 (ja) 画像処理装置
Salgo Three-dimensional echocardiography
US11704795B2 (en) Quality-driven image processing
Gallenda Reconstruction of right ventricle morphology and displacements by merging time resolved MRI series
Abboud Segmenting the Right Ventricle Cavity from 4D Echocardiography Images for Stroke Volume Measurement

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20090422

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20111206

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20120306

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20120313

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120406

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120724

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20121211

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130409

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130820

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140425

R150 Certificate of patent or registration of utility model

Ref document number: 5537026

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250