(57)【要約】【課題】医用画像データにより構築される空間内において、右心耳および左心耳を含まない心臓の存在領域を自動抽出可能な心臓存在領域抽出装置およびプログラムを得る。【解決手段】医用画像データに基づき冠状動脈の芯線Ｌを求め、その位置情報に基づき、芯線Ｌに包絡される閉曲面を、複数の放射基底関数の線形和により任意の関数を近似するＲＢＳ手法により最適化し、それをマスク用閉曲面Ｗとして設定する。このマスク用閉曲面Ｗの内部領域を心臓存在領域として抽出し、抽出された心臓存在領域に対してボリュームレンダリング法を適用し、マスク用閉曲面の内部領域に対応した心臓全体のＶＲ画像を形成する。【選択図】図７
2012245106.tif 000002
【特許請求の範囲】
【請求項１】
医用画像データに基づき構築される空間内において、心臓が存在する領域を抽出する心臓存在領域抽出装置であって、
前記医用画像データに基づき求められた、冠状血管の走行ラインの位置情報に基づき、前記走行ラインに包絡される閉曲面をマスク用閉曲面として設定するマスク用閉曲面設定手段と、
前記マスク用閉曲面の内部領域を前記心臓が存在する領域として抽出する領域抽出手段と、を備えたことを特徴とする心臓存在領域抽出装置。
【請求項２】
前記閉曲面設定手段は、前記走行ラインに沿って複数の制御点を設定する制御点設定手段と、複数の放射基底関数の線形和により任意の関数を近似する手法を用いて、前記複数の制御点を通るという条件を満足する閉曲面を最適化する閉曲面最適化手段とを有してなり、前記閉曲面最適化手段により最適化された閉曲面を前記マスク用閉曲面として設定するように構成されている、ことを特徴とする請求項１に記載の心臓存在領域抽出装置。
【請求項３】
ボリュームレンダリング法を用いて前記マスク用閉曲面の内部領域に対応したボリュームレンダリング画像を形成する画像形成手段を備えている、ことを特徴とする請求項１または２に記載の心臓存在領域抽出装置。
【請求項４】
医用画像データに基づき構築される空間内において、心臓が存在する領域を抽出する処理を、コンピュータにおいて実行せしめる心臓存在領域抽出プログラムであって、
前記医用画像データに基づき求められた、冠状血管の走行ラインの位置情報に基づき、前記走行ラインに包絡される閉曲面をマスク用閉曲面として設定するマスク用閉曲面設定ステップと、
前記マスク用閉曲面の内部領域を前記心臓が存在する領域として抽出する領域抽出ステップと、を前記コンピュータにおいて実行せしめることを特徴とする心臓存在領域抽出プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像診断を支援する技術に関するものであり、特に、Ｘ線ＣＴ（Computed Tomography）、ＭＲＩ（Magnetic Resonance Imaging）等の画像診断装置により得られた医用画像データに基づき構築される空間内において、心臓が存在する領域を抽出する心臓存在領域抽出装置およびプログラムに関する。
【背景技術】
【０００２】
Ｘ線ＣＴやＭＲＩ等の画像診断装置によって、高精度な３次元画像データ（ボリュームデータ）を短時間に得ることが可能となったことに伴い、得られた３次元画像データに基づき構築される空間内において臓器や腫瘍等の生体部位の領域を抽出し、これらを見やすく可視化して対患者や学術用としての説明に使用したり、手術計画を立てたりするなど、様々な目的で使用するという需要が高まっている。
【０００３】
従来、輪郭や病変部が比較的に抽出しやすい肺野、脳といった領域では、コンピュータ等の計算機を用いて生体領域を自動抽出するための研究が盛んに行われていたが、近年では、これまで自動抽出が困難とされていた心臓領域においても、心臓全体や冠動脈といった所定領域を自動抽出するための研究が進んでいる。
【０００４】
冠動脈の自動抽出に関しては、３次元画像データの各ボクセルに対応した画像信号値（ボクセル値）に基づき、冠動脈の芯線（中心線）の位置を特定した後、この芯線に対し垂直な平面内において冠動脈の輪郭を識別するという手法（例えば、下記特許文献１、下記非特許文献１等を参照）が開発されており、ＣＴ装置により得られた医用画像（以下「ＣＴ画像」と称する）のみならず、ＭＲＩ装置により得られた医用画像（以下「ＭＲＩ画像」と称する）においても、高精度な抽出結果が得られるようになっている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１１−２５００５号公報
【非特許文献】
【０００６】
【非特許文献１】Stephen R. Aylward, Member, IEEE and Elizabeth Bullitt, “Initialization, Noise, Singularities, and Scale in Height Ridge Traversal for Tubular Object Centerline Extraction” IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 21, NO.2, FEBRUARY 2002
【発明の概要】
【発明が解決しようとする課題】
【０００７】
一方、心臓全体を自動抽出する手法としては、造影剤を注入して撮影された胸部のＣＴ画像データを利用するものが知られている。この手法では、ＣＴ画像データにより構築された空間内において心臓の存在領域を大まかに指定し、その存在領域外に位置する骨や肺野といった不要な組織を削除するとともに、ＣＴ画像データの各データ値（ＣＴ値）等に基づき、造影剤が注入された冠動脈部分とそれに隣接する軟組織部分（心筋や右心系）を選択し、これを心臓として抽出するようにしている。
【０００８】
この手法によれば、冠動脈の一部に覆い被さるように存在する右心耳や左心耳も心臓の一部として抽出される（後述の図９（Ｂ）参照）。右心耳および左心耳を心臓の一部とすることは解剖学的には正しいことであるが、冠動脈の走行状態を観察したい場合においては、右心耳や左心耳の存在が冠動脈の走行確認の妨げとなるため、これらを画像中から削除したいという要望がある。
【０００９】
しかしながら、右心耳および左心耳は、冠動脈と密着するように存在しており、また、心筋等の他の心臓組織とＣＴ値が類似しているため、右心耳および左心耳を含まない状態の心臓領域を自動的に抽出することは困難であった。このため、従来は、解剖学的な専門知識を有する医師等が、目視に基づく手動操作により、右心耳および左心耳を画像中から削除することも行われているが、このような人為的操作が必要とされることにより、画像生成に多大な時間を要するとともに、得られる冠動脈の走行情報等に悪影響を及ぼす虞が生じる。
【００１０】
また、上述した心臓の自動抽出手法は、撮影機器や撮影条件に影響を受けることなく所定の画像信号値（ＣＴ値）が得られるＣＴ画像データに適用し得るように開発されたものであるため、撮影機器や撮影条件によって画像信号値に大きな変動が生じるＭＲＩ画像データに適用した場合には、安定した抽出結果を得ることは困難であった（後述の図１０（Ｂ）参照）。
【００１１】
本発明は、このような事情に鑑みなされたものであり、ＣＴ装置またはＭＲＩ装置により得られた医用画像データに基づき構築される空間内において、右心耳および左心耳を含まない心臓の存在領域を自動的に抽出することが可能な心臓存在領域抽出装置およびプログラムを提供することを目的とする。
【課題を解決するための手段】
【００１２】
上記目的を達成するため、本発明に係る心臓存在領域抽出装置およびプログラムは、以下のように構成されている。
【００１３】
すなわち、本発明に係る心臓存在領域抽出装置は、
医用画像データ（例えば、ＣＴ画像データやＭＲＩ画像データ）に基づき構築される空間内において、心臓が存在する領域を抽出する心臓存在領域抽出装置であって、
前記医用画像データに基づき求められた、冠状血管の走行ラインの位置情報に基づき、前記走行ラインに包絡される閉曲面をマスク用閉曲面として設定するマスク用閉曲面設定手段と、
前記マスク用閉曲面の内部領域を前記心臓が存在する領域として抽出する領域抽出手段と、を備えたことを特徴とする。
【００１４】
本発明に係る心臓存在領域抽出装置において、
前記閉曲面設定手段は、前記走行ラインに沿って複数の制御点を設定する制御点設定手段と、複数の放射基底関数の線形和により任意の関数を近似する手法を用いて、前記複数の制御点を通るという条件を満足する閉曲面を最適化する閉曲面最適化手段とを有してなり、前記閉曲面最適化手段により最適化された閉曲面を前記マスク用閉曲面として設定するように構成されている、とすることができる。
【００１５】
また、ボリュームレンダリング法を用いて前記マスク用閉曲面の内部領域に対応したボリュームレンダリング画像を形成する画像形成手段を備えるようにしてもよい。
【００１６】
また、本発明に係る心臓存在領域抽出プログラムは、
医用画像データに基づき構築される空間内において、心臓が存在する領域を抽出する処理を、コンピュータにおいて実行せしめる心臓存在領域抽出プログラムであって、
前記医用画像データに基づき求められた、冠状血管の走行ラインの位置情報に基づき、前記走行ラインに包絡される閉曲面をマスク用閉曲面として設定するマスク用閉曲面設定ステップと、
前記マスク用閉曲面の内部領域を前記心臓が存在する領域として抽出する領域抽出ステップと、を前記コンピュータにおいて実行せしめることを特徴とする。
【００１７】
上記「冠状血管」とは、心臓の外表面に沿って走行する血管を意味しており、冠動脈の他に、冠静脈をも含む概念である。
【００１８】
上記「冠状血管の走行ライン」とは、冠動脈等の冠状血管の３次元的な走行状態を示す各ラインを意味しており、冠状血管の芯線（中心線）の他に、冠状血管の表面に沿ったライン等をも含む概念である。
【００１９】
上記「走行ラインに包絡される閉曲面」とは、各走行ラインの一部と接しながら、各走行ラインの全体により内包される閉曲面を意味しているが、閉曲面の内部に走行ラインの一部が含まれる場合を排除するものではない。
【発明の効果】
【００２０】
本発明に係る心臓存在領域抽出装置およびプログラムは、医用画像データに基づき求められた、冠状血管の走行ラインの位置情報に基づき、走行ラインに包絡される閉曲面をマスク用閉曲面として設定し、このマスク用閉曲面の内部領域を心臓が存在する領域として抽出するものであり、これにより、以下のような作用効果を奏する。
【００２１】
すなわち、冠動脈等の冠状血管は心臓の外面に沿うように走行しているため、冠状血管の走行ラインに包絡される閉曲面は、心臓の外形に近似した形状を有するものとなり、その内部領域には、心臓の主要部（心筋部分）が含まれることとなる。一方、右心耳および左心耳は解剖学的に冠動脈の一部を覆う位置、すなわち冠動脈の走行ラインの外側に位置するため、閉曲面の内部領域には含まれないこととなる。
【００２２】
また、ＣＴ画像データであってもＭＲＩ画像データであっても、冠状血管の走行ラインの位置情報を自動的に計算することは可能とされており、さらに、走行ラインに包絡される閉曲面を求めることについても、画像処理技術において従来用いられてきた３次元曲面の生成手法を適用することにより、自動的に計算することが可能である。しかも、この閉曲面を求めるための演算を実行するには、走行ラインの位置情報が必要とされるだけで、医用画像データの画像信号値（ボクセル値）の分布情報等は必要とされない。このため、ＣＴ画像データであってもＭＲＩ画像データであっても、閉曲面を自動的に求めることが可能である。
【００２３】
したがって、冠状血管の走行ラインに包絡される閉曲面をマスク用閉曲面として設定し、このマスク用閉曲面の内部領域を心臓が存在する領域として抽出するという特徴構成を有する本発明によれば、ＣＴ画像データまたはＭＲＩ画像データに基づき構築される空間内において、右心耳および左心耳を含まない心臓の存在領域を自動的に抽出することが可能となる。
【００２４】
また、本発明によれば、右心耳および左心耳を画像中から削除するために従来行われていた人為的操作が不要となるので、効率的かつ安定的に、適正な情報（冠動脈の走行情報等）を担持した心臓抽出画像を生成することが可能となる。
【図面の簡単な説明】
【００２５】
【図１】本発明の一実施形態に係る心臓存在領域抽出装置の概略構成を示すブロック図である。
【図２】図１に示す演算処理装置の機能的構成を示すブロック図である。
【図３】図１に示す演算処理装置が実行する処理の流れを示すフローチャートである。
【図４】マスク用閉曲面設定処理の流れを示すフローチャートである。
【図５】求められた冠動脈の走行ラインを示す図である。
【図６】冠動脈の走行ライン上に設定される複数の制御点の設定態様を示す図である。
【図７】冠動脈の走行ラインに包絡される閉曲面を示す図である。
【図８】冠動脈抽出画像を示す図である。
【図９】ＣＴ画像データに基づく心臓抽出画像の一例を示す図（（Ａ）は本発明手法によるものであり、（Ｂ）は従来手法によるもの）である。
【図１０】ＭＲＩ画像データに基づく心臓抽出画像の一例を示す図（（Ａ）は本発明手法によるものであり、（Ｂ）は従来手法によるもの）である。
【発明を実施するための形態】
【００２６】
以下、本発明の実施形態について上述の図面を参照しつつ詳細に説明する。まず、図１および図２を参照しながら本発明の一実施形態に係る心臓存在領域抽出装置について説明する。
【００２７】
図１に示す心臓存在領域抽出装置１は、医用画像データに基づき構築される空間（以下「画像データ空間」と称する）内において心臓が存在する領域を抽出し、この領域内に位置する心臓の全体画像を形成するものであり、コンピュータ等からなる演算処理装置１０と、マウスやキーボード等からなる入力装置２０と、画像表示装置等からなる出力装置３０とを備えてなる。
【００２８】
上記演算処理装置１０は、図２に示すように、医用画像データ記憶手段１１、走行ライン特定手段１２、閉曲面設定手段１３、領域抽出手段１４および画像形成手段１５を備えており、閉曲面設定手段１３は制御点設定手段１３１および閉曲面最適化手段１３２を有している。これらの手段は、各種の演算処理を行うＣＰＵ、ハードディスクやＲＯＭ等の記憶装置、該記憶装置に格納された制御プログラム（本発明の一実施形態に係る心臓存在領域抽出プログラムを含む）およびＲＡＭ等の一時記憶装置等により構成されるものを、概念的に示したものである。
【００２９】
上記医用画像データ記憶手段１１は、Ｘ線ＣＴ装置またはＭＲＩ装置により得られた、被検体胸部の３次元医用画像データ（ＣＴ画像データまたはＭＲＩ画像データ）と、上述の走行ライン特定手段１２、閉曲面設定手段１３、領域抽出手段１４および画像形成手段１５により得られた各種データと、を記憶するものである。
【００３０】
上記走行ライン特定手段１２は、上記画像データ空間内において、心臓の冠状血管（本実施形態では冠動脈）の全体的な走行ラインを特定し、その３次元的な位置情報を求めるものである。
【００３１】
上記閉曲面設定手段１３は、上記走行ライン特定手段１２により求められた、冠動脈の走行ラインの位置情報に基づき、この走行ラインに包絡される閉曲面を求め、これをマスク用閉曲面として設定するものである。具体的には、上記制御点設定手段１３１により、上記走行ラインに沿った複数の制御点を設定し、この複数の制御点を通るという条件を満足する閉曲面を、上記閉曲面最適化手段１３２により最適化し、この最適化された閉曲面を上記マスク用閉曲面として設定するように構成されている。
【００３２】
上記領域抽出手段１４は、上記閉曲面設定手段１３により設定された上記マスク用閉曲面の内部領域を、心臓が存在する領域として抽出するものである。
【００３３】
上記画像形成手段１５は、ボリュームレンダリング法を用いて、上記マスク用閉曲面の内部領域に対応した心臓全体のボリュームレンダリング画像（以下「ＶＲ画像」と称する）を形成するものである。
【００３４】
次に、本実施形態の心臓存在領域抽出装置１による心臓存在領域の抽出手順および心臓全体のＶＲ画像の形成手順について、主に図３〜図１０を参照しながら説明する。なお、次述する心臓存在領域の抽出手順は、本発明の一実施形態に係る心臓存在領域抽出プログラムに従って実行されるものである。
【００３５】
〈１〉被検体胸部の医用画像データ（ＣＴ画像データまたはＭＲＩ画像データ）を取得する（医用画像データ取得ステップ；図３のステップＳ１参照）。本実施形態では、必要な医用画像データが、予め上記医用画像データ記憶手段１１（図２参照）に記憶されていることを想定しているが、ＣＤやＤＶＤ等の情報記憶媒体に記録された医用画像データを必要に応じて読み出すようにしてもよい。
【００３６】
〈２〉取得された医用画像データに基づき、該医用画像データにより構築される画像データ空間内において、心臓冠動脈の全体的な走行ラインを特定し、その３次元的な位置情報を求める（走行ライン特定ステップ；図３のステップＳ２、図５参照）。本実施形態では、前掲の非特許文献１に開示された手法により、冠動脈の芯線（中心線）を走行ラインとして特定する。
【００３７】
具体的には、まず、医用画像データの各データ値（ボクセル値）に基づき、画像データ空間内の局所領域において冠動脈の芯線上に位置すると推定される点（以下「経由点」と称する）を設定する。次に、経由点を含む局所領域におけるヘッセ（Hesse）行列の固有値を算出し、算出された固有値に基づき冠動脈の芯線の主軸方向（芯線の延びる方向）を求め、その主軸方向に新たな経由点を設定する。この手順を繰り返すことにより多数の経由点を順次設定し、設定された経由点同士を滑らかに連結することにより、冠動脈の全体的な芯線Ｌ（図５参照）を特定する。
【００３８】
なお、解剖学的に冠動脈は、大動脈起始部の膨大部から分岐する右冠動脈（ＲＣＡ）と左冠動脈（ＬＣＡ）とに大別される。右冠動脈は、右房室間溝を通り心臓の右辺縁を下って後室間溝で後下行枝と房室枝とに枝分かれし、それぞれの枝がさらに細かく枝分かれしている。左冠動脈はその根元近傍部から心臓の前側に下る左前下行枝（ＬＡＤ）と背中側に回り込む左回旋枝（ＬＣＸ）とに枝分かれし、それぞれの枝がさらに細かく枝分かれしている。上記芯線Ｌの特定は、冠動脈の各枝の先端部に至るまで詳細に行うことが好ましいが、主要な枝のみについて行うなど、後述の閉曲面の特定に支障がない範囲で適宜省略してもよい。
【００３９】
〈３〉求められた冠動脈の芯線Ｌの位置情報に基づき、この芯線Ｌに包絡される閉曲面を求め、これをマスク用閉曲面として設定するマスク用閉曲面設定処理を実行する（マスク用閉曲面設定ステップ；図３のステップＳ３参照）。
【００４０】
具体的には、まず、芯線Ｌ上に複数の制御点Ｐを設定する（制御点設定ステップ；図４のステップＳ３１、図６参照）。制御点Ｐの設定は、芯線Ｌに沿って任意の間隔で適宜行うことが可能であるが、本実施形態では、芯線Ｌの各部の曲率に応じて、曲率が大きい部位では間隔を狭く曲率が小さい部位では間隔を広げ、平均的には実寸換算で数ミリメートル置きに制御点Ｐを自動設定する。なお、図６では、概略的に制御点Ｐを「・」（一部のみに記号「Ｐ」を添付）で示しているが、実際には、図示した状態よりも高密度に制御点Ｐが設定される。
【００４１】
次に、複数の放射基底関数（Radial Basis Functions）の線形和により任意の関数を近似する手法（以下「ＲＢＦ手法」と称する）を用いて、上述の複数の制御点Ｐを通るという条件を満足する閉曲面を最適化し、この最適化された閉曲面をマスク用閉曲面Ｗとして設定する（閉曲面最適化ステップ；図４のステップＳ３２、図７参照）。なお、図７では、マスク用閉曲面Ｗの外形線を一点鎖線により簡略化して示している。
【００４２】
ＲＢＦ手法は、３次元空間内の任意の点ｘ（ｘはベクトル（３次元点））の関数ｓ(ｘ)を用いて、ｓ(ｘ)＝０の陰関数表現により、任意の曲面を近似表現する手法であり、自由度の高い様々な曲面表現が可能であることや、曲面の局所平滑性を容易に調整することができるといった利点を有している。
【００４３】
ＲＢＦ手法に用いられる関数ｓ(ｘ)は、下式（１）によって表される。
【００４４】
【数１】
2012245106.tif 000003
【００４５】
放射基底関数φ（ｒ）としては、thin-plate関数、ガウス関数、multi-quadratic関数、bi-harmonic関数、tri-harmonic関数等を用いることが可能であるが、本実施形態では、放射基底関数φ（ｒ）としてbi-harmonic関数を用いている。この場合、放射基底関数φ（ｒ）は、下式（２）で表され、上記多項式関数ｐ(ｘ)は下式（３）で表される。
【００４６】
【数２】
2012245106.tif 000004
【００４７】
本実施形態において、上述の複数の制御点Ｐを通るという条件を満足するマスク用閉曲面Ｗ（図７参照）を求めるには、上式（１）における各点ｘの座標値（ｘ，ｙ，ｚ）として各制御点Ｐの座標値を代入し、下式（４）および（５）により定義されるベクトルΛおよびｅを求める演算処理を行えばよい。
【００４８】
【数３】
2012245106.tif 000005
【００４９】
この演算処理は、下式（６）の追加条件の下で、下式（７）の線形方程式を解くことに等しいことが知られている（下記参考文献１参照）。
【００５０】
【数４】
2012245106.tif 000006
【００５１】
参考文献１：Yoshitaka Masutani Image Computing and Analysis Laboratory, Department of Radiology, University of Tokyo(UT-RAD/ICAL) “RBF-based Representation of Volumetric Data: Application in Visualization and Segmentation”, MICCAI 2002, LNCS 2489, pp.300-307
【００５２】
〈４〉上述のＲＢＦ手法による演算処理により求められたマスク用閉曲面Ｗの内部領域を心臓が存在する領域（以下「心臓存在領域」と称する）として抽出する（領域抽出ステップ；図３のステップＳ４参照）。
【００５３】
〈５〉抽出された心臓存在領域に対してボリュームレンダリング法を適用し、上記マスク用閉曲面Ｗの内部領域に対応した心臓全体のＶＲ画像を形成する（ＶＲ画像形成ステップ；図３のステップＳ５参照）。なお、ボリュームレンダリング法とは、３次元医用画像データのボクセル毎に、画像信号値に対応させて色や不透明度を付与して画像化する手法である。
【００５４】
本実施形態では、冠動脈の芯線Ｌを走行ラインとして用いているため、上記マスク用閉曲面Ｗは、冠動脈を径方向に２分する位置を通過するものとなる。このため、上述のＶＲ画像形成ステップにおいて形成されるＶＲ画像（以下「心臓全体画像」と称する）では、芯線Ｌに沿って切断された状態の冠動脈が描出されることとなる。冠動脈の表面部分が描出された画像を形成したい場合には、上述の走行ライン特定ステップにおいて芯線Ｌを抽出する際に、前掲の非特許文献１に開示された手法に基づき、冠動脈のみが抽出されたＶＲ画像（以下「冠動脈抽出画像」と称する）を形成しておき、この冠動脈抽出画像（図８参照）を、上記心臓全体画像と合成すればよい。あるいは、芯線Ｌに沿って上記マスク用閉曲面Ｗの一部を拡張し、上記心臓存在領域内に冠動脈の全体が含まれるように設定することにより、冠動脈の表面部分が描出された画像を形成するようにしてもよい。
【００５５】
以上説明した本発明手法によれば、ＣＴ画像データまたはＭＲＩ画像データに基づき構築される画像データ空間内において、右心耳および左心耳を含まない心臓の存在領域を自動的に抽出することができ、冠動脈の走行状態を観察するのに適した心臓全体の抽出画像を形成することが可能である。
【００５６】
具体的には、ＣＴ画像データを用いた場合、従来手法では、図９（Ｂ）に示すように、右心耳３２や左心耳３３、肺動脈３４といった、冠動脈の走行状態を観察する際には邪魔となる組織（以下「観察時不要組織」と称する）が描出されているのに対し、本発明手法によれば、図９（Ａ）に示すように、これらの観察時不要組織が描出されていない心臓全体の抽出画像を形成することができる。
【００５７】
なお、図９（Ａ）は、マスク用閉曲面Ｗの内部領域に対応した心臓全体画像と、前掲の非特許文献１に開示された手法に基づき形成された冠動脈抽出画像とを合成することにより形成されたものである。図９（Ａ）では、大動脈３１の一部が描出されているが、これは、冠動脈抽出画像を形成する際に、大動脈の一部が抽出されていることによる。このことは、次述する図１０（Ａ）についても同様である。
【００５８】
また、ＭＲＩ画像データを用いた場合、従来手法では、図１０（Ｂ）に示すように、心臓の周囲に存在する観察時不要組織が多数描出されてしまう可能性が高いのに対し、本発明手法によれば、図１０（Ａ）に示すように、これらの観察時不要組織（右心耳、左心耳および肺動脈を含む）が描出されていない心臓全体の抽出画像を形成することができる。
【００５９】
以上、本発明の一実施形態について説明したが、本発明は上記実施形態に限定されるものではなく、種々に態様を変更することが可能である。
【００６０】
例えば上記実施形態では、冠動脈の芯線（走行ライン）に包絡されるマスク用閉曲面を求める手法として、上述のＲＢＦ手法を用いているが、スプライン関数を用いて曲面を近似する手法など、他の３次元曲面生成手法を用いてマスク用閉曲面を求めるようにしてもよい。
【００６１】
また、上記実施形態では、冠動脈の走行ラインとしての芯線を求めるために、前掲の非特許文献１に開示された手法を用いているが、他の手法を用いて、冠動脈の走行ラインを求めるようにしてもよい。
【符号の説明】
【００６２】
１  心臓存在領域抽出装置
１０  演算処理装置
２０  入力装置
３０  出力装置
１１  医用画像データ記憶手段
１２  走行ライン特定手段
１３  閉曲面設定手段
１４  領域抽出手段
１５  画像形成手段
３１  大動脈
３２  右心耳
３３  左心耳
３４  肺動脈
１３１  制御点設定手段
１３２  閉曲面最適化手段
Ｌ  芯線（走行ライン）
Ｐ  制御点
Ｗ  マスク用閉曲面
【図１】
2012245106.tif 000007
【図２】
2012245106.tif 000008
【図３】
2012245106.tif 000009
【図４】
2012245106.tif 000010
【図５】
2012245106.tif 000011
【図６】
2012245106.tif 000012
【図７】
2012245106.tif 000013
【図８】
2012245106.tif 000014
【図９】
2012245106.tif 000015
【図１０】
2012245106.tif 000016