JP5892548B2 - 描画プログラム、描画装置、および描画方法 - Google Patents

Description

本発明は、描画プログラム、描画装置、および描画方法に関する。

人間や動物の心臓は血液を拍出する臓器であり、心臓の機能を再現する数値解析が従来から行われている。心臓の数値解析では、心筋の任意領域が時刻変化する心筋梗塞のシミュレーションが行われている。このような心筋の任意領域の変化を観察する場合、心筋の任意領域の抽出技術が利用可能である（たとえば、下記非特許文献１を参照。）。

Ｋｉｔｗａｒｅ，Ｉｎｃ．著 ｖｔｋＭａｒｃｈｉｎｇＣｏｎｔｏｕｒＦｉｌｔｅｒ Ｃｌａｓｓ Ｒｅｆｅｒｅｎｃｅ．ｐｄｆ ２０１０年９月２７日発行 ［ＵＲＬ］ｗｗｗ．ｖｔｋ．ｏｒｇ／ｄｏｃ／ｒｅｌｅａｓｅ／５．６／ｈｔｍｌ／ａ０１１０３．ｈｔｍｌ ［検索日］２０１２年３月１２日

しかしながら、上述した従来技術では、任意領域内の要素が持つ仕事率や仕事力などの物理値の上限および下限の許容範囲を心筋梗塞の領域として、ユーザが指定することになる。心筋梗塞の領域を特定する物理値の許容範囲は、患者や患者の状態により異なるため、心筋梗塞の許容範囲を指定するのは困難である。このような問題は、心臓に限らず他の臓器についても同様であり、心筋梗塞以外の他の疾患についても同様である。

本発明は、臓器における疾患部位の特定の容易化を図ることを目的とする。

本発明の一側面によれば、シミュレーション対象の位置毎に物理値をそれぞれ有する要素の集合であるシミュレーションモデルにおいて指定された領域の描画対象である描画対象要素群が有する物理値の範囲を検出し、前記描画対象要素群に含まれる各描画対象要素に隣接する要素の集合である隣接要素群を前記シミュレーションモデルから特定し、特定された前記隣接要素群に含まれる各隣接要素の物理値が前記範囲に含まれるかを判断し、物理値が前記範囲に含まれると判断された隣接要素を前記描画対象要素群に追加させ、隣接要素が追加された描画対象要素群に応じて、前記シミュレーションモデルを描画する描画プログラム、描画装置、および描画方法が提案される。

本発明の一側面によれば、臓器における疾患部位の特定の容易化を図ることができるという効果を奏する。

図１は、実施の形態にかかる描画装置による心臓モデルの描画例を示す説明図である。 図２は、実施の形態にかかる描画装置のハードウェア構成例を示すブロック図である。 図３は、実施の形態にかかる描画装置の機能的構成例を示すブロック図である。 図４は、描画対象要素群のデータ構造例を示す説明図である。 図５は、図４に示したテトラ要素ｔｉを示す説明図である。 図６は、ヒストグラムの一例を示す説明図である。 図７は、隣接要素群のデータ構造例を示す説明図である。 図８は、追加要素群のデータ構造例を示す説明図である。 図９は、実施の形態にかかる描画装置による描画処理手順例を示すフローチャートである。 図１０は、図９の拡張処理（ステップＳ９０７）の詳細な処理手順例を示すフローチャートである。

以下に添付図面を参照して、本発明にかかる描画プログラム、描画装置、および描画方法の実施の形態を詳細に説明する。なお、本実施の形態では、臓器の一例として人間の心臓の３次元モデルを用いて説明するが、心臓以外の人間や動物の他の臓器の３次元モデルでも実現可能である。

＜心臓モデルの描画例＞
図１は、実施の形態にかかる描画装置による心臓モデルの描画例を示す説明図である。図１において、（Ａ）は心臓の３次元モデルを示す。心臓の３次元モデルを、以下、「心臓モデル」と称す。心臓モデル１００は、たとえば、テトラ要素と呼ばれる四面体の要素（非構造格子データ）の集合であり、各テトラ要素は、当該テトラ要素の配置位置に応じた物理値を持つ。物理値は、テトラ要素に対応する心筋の挙動を示す値であり、たとえば、仕事力［ＫＰａ］や、仕事［Ｊ／ｍｌ］、仕事率［Ｊ／ｓ・ｍｌ］が物理値に採用される。

描画装置は、ユーザによるマウス等のポインティングデバイスの操作により、心臓モデル１００においてある領域を指定する。指定された領域を「指定領域」と称す。具体的には、たとえば、描画装置は、ユーザの操作により指定された矩形領域を心臓モデル１００の表面に写像することにより、当該矩形領域に対応する心臓モデル１００の表面を指定領域Ａに設定する。

（Ｂ）は、（Ａ）で指定された指定領域Ａの拡大図である。指定領域Ａはテトラ要素群を含む。指定領域Ａ内のテトラ要素群を描画対象要素群Ｕと称す。（Ｂ）において、三角形の図形はテトラ要素を示す。テトラ要素は四面体であるが、ここでは便宜的に三角形とする。描画装置は、指定領域Ａを、太矢印で示したテトラ要素の隣接方向に拡大していく。

この場合、（Ｂ）の初期の描画対象要素群Ｕにおける物理値の上限と下限との間の許容範囲に物理値をもつ隣接要素が、描画対象要素群Ｕに含まれる。また、当該許容範囲外の隣接要素であっても、後述する一定条件下で描画対象要素群Ｕに含まれる。（Ｃ）は、拡大後の描画対象要素群Ｕである。拡張後の描画対象要素群Ｕが、観測したい心筋梗塞等の病態を表示した領域となる。

（Ｄ）は、（Ｃ）の拡張後の描画対象要素群Ｕを特徴づけて心臓モデル１００を描画した描画例である。このように、（Ａ）において、指定領域Ａを指定するだけで、患者や患者の状態により変化する心筋梗塞等の病態の領域を容易に特定することができる。

＜描画装置のハードウェア構成例＞
図２は、実施の形態にかかる描画装置のハードウェア構成例を示すブロック図である。図２において、描画装置は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）２０１と、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）２０２と、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）２０３と、磁気ディスクドライブ２０４と、磁気ディスク２０５と、光ディスクドライブ２０６と、光ディスク２０７と、ディスプレイ２０８と、Ｉ／Ｆ（Ｉｎｔｅｒｆａｃｅ）２０９と、キーボード２１０と、マウス２１１と、スキャナ２１２と、プリンタ２１３と、を備えている。また、各構成部はバス２００によってそれぞれ接続されている。

ここで、ＣＰＵ２０１は、描画装置の全体の制御を司る。ＲＯＭ２０２は、ブートプログラムなどのプログラムを記憶している。ＲＡＭ２０３は、ＣＰＵ２０１のワークエリアとして使用される。磁気ディスクドライブ２０４は、ＣＰＵ２０１の制御にしたがって磁気ディスク２０５に対するデータのリード／ライトを制御する。磁気ディスク２０５は、磁気ディスクドライブ２０４の制御で書き込まれたデータを記憶する。

光ディスクドライブ２０６は、ＣＰＵ２０１の制御にしたがって光ディスク２０７に対するデータのリード／ライトを制御する。光ディスク２０７は、光ディスクドライブ２０６を介して書き込まれたデータを記憶したり、光ディスクドライブ２０６を介して記憶されたデータをコンピュータに読み取らせたりする。

ディスプレイ２０８は、カーソル、アイコンあるいはツールボックスをはじめ、文書、画像、機能情報などのデータを表示する。このディスプレイ２０８は、たとえば、液晶ディスプレイ、プラズマディスプレイなどを採用することができる。

インターフェース（以下、「Ｉ／Ｆ」と略する。）２０９は、通信回線を通じてＬＡＮ（Ｌｏｃａｌ Ａｒｅａ Ｎｅｔｗｏｒｋ）、ＷＡＮ（Ｗｉｄｅ Ａｒｅａ Ｎｅｔｗｏｒｋ）、インターネットなどのネットワーク２１４に接続され、このネットワーク２１４を介して他の装置に接続される。そして、Ｉ／Ｆ２０９は、ネットワーク２１４と内部のインターフェースを司り、外部装置からのデータの入出力を制御する。Ｉ／Ｆ２０９には、たとえばモデムやＬＡＮアダプタなどを採用することができる。

キーボード２１０は、文字、数字、各種指示などの入力のためのキーを備え、データの入力をおこなう。また、タッチパネル式の入力パッドやテンキーなどであってもよい。マウス２１１は、カーソルの移動や許容範囲選択、あるいはウィンドウの移動やサイズの変更などをおこなう。ポインティングデバイスとして同様に機能を備えるものであれば、トラックボールやジョイスティックなどであってもよい。

スキャナ２１２は、画像を光学的に読み取り、描画装置内に画像データを取り込む。なお、スキャナ２１２は、ＯＣＲ（Ｏｐｔｉｃａｌ Ｃｈａｒａｃｔｅｒ Ｒｅａｄｅｒ）機能を持たせてもよい。また、プリンタ２１３は、画像データや文書データを印刷する。プリンタ２１３には、たとえば、レーザプリンタやインクジェットプリンタを採用することができる。なお、光ディスクドライブ２０６、光ディスク２０７、ディスプレイ２０８、キーボード２１０、マウス２１１、スキャナ２１２、およびプリンタ２１３の少なくともいずれか１つは、なくてもよい。

＜描画装置の機能的構成例＞
図３は、実施の形態にかかる描画装置の機能的構成例を示すブロック図である。描画装置３００は、検出部３０１と、特定部３０２と、判断部３０３と、追加部３０４と、描画部３０５と、を有する。検出部３０１、特定部３０２、判断部３０３、追加部３０４、および描画部３０５は、具体的には、たとえば、図２に示したＲＯＭ２０２、ＲＡＭ２０３、磁気ディスク２０５、光ディスク２０７などの記憶装置に記憶されたプログラムをＣＰＵ２０１に実行させることにより、または、Ｉ／Ｆ２０９により、それらの機能を実現する。

検出部３０１は、シミュレーション対象の位置毎に物理値をそれぞれ有する要素の集合であるシミュレーションモデルにおいて指定された領域の描画対象である描画対象要素群Ｕが有する物理値の範囲を検出する。具体的には、検出部３０１は、各々臓器での位置に応じた物理値を有する要素の集合である臓器モデルにおける指定領域Ａ内の描画対象要素群Ｕにおける物理値の許容範囲を検出する。臓器モデルにおける指定領域Ａ内の描画対象要素群Ｕのデータ構造は、たとえば、心臓モデル１００を構成するテトラ要素群のデータ構造の一部である。

図４は、描画対象要素群Ｕのデータ構造例を示す説明図である。描画対象要素群Ｕのデータ構造４００は、ＩＤ項目と、第１頂点項目〜第４頂点項目と、重心項目と、物理値項目と、を有し、テトラ要素ごとに、各項目の値を有するデータ構造である。

ＩＤ項目には、指定領域Ａ内のテトラ要素のＩＤ：ｉ（１≦ｉ≦ｎ）が格納される。ＩＤとは、テトラ要素を一意に識別する識別情報である。ＩＤ：ｉのテトラ要素を、テトラ要素ｔｉとする。第１頂点項目〜第４頂点項目には、テトラ要素ｔｉの第１頂点ｖ１（ｔｉ）〜第４頂点ｖ４（ｔｉ）の頂点座標値が格納される。便宜上、座標値についても、ｖ１（ｔｉ）〜ｖ４（ｔｉ）と表記する。第１頂点ｖ１（ｔｉ）〜第４頂点ｖ４（ｔｉ）とは、四面体であるテトラ要素の各頂点である。

重心項目には、テトラ要素ｔｉの重心ｇ（ｔｉ）となる重心座標値が格納される。便宜上、座標値についても、ｇ（ｔｉ）と表記する。物理値項目には、テトラ要素ｔｉに割り当てられた物理値ｅ（ｔｉ）が格納される。上述したように、物理値は、テトラ要素に対応する心筋の挙動を示す値であり、たとえば、仕事力［ＫＰａ］や、仕事［Ｊ／ｍｌ］、仕事率［Ｊ／ｓ・ｍｌ］が採用される。

図５は、図４に示したテトラ要素ｔｉを示す説明図である。テトラ要素ｔｉの形状は第１頂点ｖ１（ｔｉ）〜第４頂点ｖ４（ｔｉ）を有する四面体である。また、重心ｇ（ｔｉ）には、物理値ｅ（ｔｉ）が割り当てられる。

図３に戻る。検出部３０１は、具体的には、図１の（Ａ）に示したように、指定領域Ａ内の描画対象要素群Ｕの各テトラ要素が持つ物理値を抽出する。たとえば、検出部３０１は、描画対象要素群Ｕの各テトラ要素が持つ物理値のヒストグラムを作成し、物理値の上限と下限を検出する。

図６は、ヒストグラムの一例を示す説明図である。検出部３０１は、たとえば、ヒストグラム６００において、上限ＨＬから下限ＬＬまでの許容範囲６０１を検出する。ここで、上限ＨＬとは、物理値の最大値であり、下限ＬＬとは、物理値の最小値である。また、上限ＨＬについては、要素数が所定数以上となる物理値の最大値としてもよい。同様に、下限ＬＬについても、要素数が所定数以上となる物理値の最小値としてもよい。また、あらかじめ設定された許容値により許容範囲６０１を拡大してもよい。たとえば、許容範囲６０１が［ＬＬ，ＨＬ］の場合、許容値α（＞０）を考慮すると、許容範囲６０１が［ＬＬ−α，ＨＬ＋α］となる。

図３に戻る。特定部３０２は、描画対象要素群Ｕに含まれる各描画対象要素に隣接する要素の集合である隣接要素群を前記シミュレーションモデルから特定する。具体的には、特定部３０２は、描画対象要素群Ｕの各要素に隣接する要素の集合である隣接要素群を臓器モデルから特定する。たとえば、特定部３０２は、隣接要素群を心臓モデル１００から特定する。隣接要素とは、隣接元のテトラ要素に隣接するテトラ要素であり、たとえば、隣接元のテトラ要素の面、辺、頂点のいずれかに接触するテトラ要素である。特定された隣接要素群は、図７に示すデータ構造で、記憶装置に保持される。

図７は、隣接要素群のデータ構造例を示す説明図である。隣接要素群のデータ構造７００は、描画対象要素群Ｕのデータ構造４００を構成する項目群に、選択項目が追加されたデータ構造である。描画対象要素群Ｕのデータ構造４００を構成する項目群と同一項目については説明を省略する。

選択項目には、その隣接要素が選択済みか未選択かを示す選択フラグが格納される。ここでは、例として、「１」が選択済み、「０」が未選択とする。図７では、すべて「０」であるため、まだ１つの隣接要素も選択されていない初期状態を示している。

図３に戻る。判断部３０３は、特定された隣接要素群に含まれる各隣接要素の物理値が前記範囲に含まれるかを判断する。具体的には、判断部３０３は、特定部３０２によって特定された隣接要素群の各要素の物理値が許容範囲６０１内であるか否かを判断する。具体的には、たとえば、判断部３０３は、図７の隣接要素群のデータ構造７００から１レコードずつ隣接要素を選択する。選択された隣接要素の選択フラグは、「０」から「１」に遷移する。そして、判断部３０３は、選択された隣接要素の物理値が、検出部３０１で検出された許容範囲６０１内の値か否かを判断する。

追加部３０４は、物理値が前記範囲に含まれると判断された隣接要素を描画対象要素群Ｕに追加する。具体的には、追加部３０４は、判断部３０３によって許容範囲６０１内と判断された要素を描画対象要素群Ｕに追加する。たとえば、追加部３０４は、物理値が許容範囲６０１内と判断された隣接要素を、追加要素として、追加要素群のデータ構造に追加する。

図８は、追加要素群のデータ構造例を示す説明図である。追加要素群のデータ構造８００は、描画対象要素群Ｕのデータ構造４００を構成する項目群と同一項目であるため、説明を省略する。追加要素ｗｋは、隣接要素群のうち物理値が許容範囲６０１内であるテトラ要素である。

図３に戻る。追加部３０４は、隣接要素群のうち許容範囲６０１外と判断された要素群の要素が、描画対象要素群Ｕおよび隣接要素群の中の複数の要素に接する場合、描画対象要素群Ｕに追加する。具体的には、たとえば、追加部３０４は、物理値が許容範囲６０１外の隣接要素であっても、追加後の描画対象要素群Ｕの中の複数のテトラ要素に囲まれている場合は、物理値が許容範囲６０１内にあるとみなして、追加要素として、追加要素群のデータ構造８００に追加する。

心臓モデル１００を生成するミュレーションでは、テトラ要素の細分化により、物理値が許容範囲６０１外のテトラ要素が、物理値が許容範囲６０１内のテトラ要素群の中に離散的に存在する場合がある。一方、ユーザが観測したい心筋梗塞などの疾患部位は連続的である。したがって、描画対象要素群Ｕの連続性を担保するため、物理値が許容範囲６０１外の隣接要素であっても、描画対象要素群Ｕ（追加されたものも含む）の複数のテトラ要素に接する場合は、追加要素とする。

ここで、複数のテトラ要素との接触については、たとえば、面どうしが接触する場合、辺どうしが接触する場合、頂点どうしが接触する場合がある。どの接触をもって複数のテトラ要素との接触と規定するかは、設定により決まる。なお、面、辺、頂点のうち、面が最も接触部位が多く、頂点が最も接触部位が少ない。したがって、複数のテトラ要素との接触については、いずれのテトラ要素とも面のみ接触する場合に、複数のテトラ要素との接触としてもよい。また、いずれのテトラ要素とも面または辺で接触する場合に、複数のテトラ要素との接触としてもよい。また、面、辺または頂点のいずれかで接触する場合に、複数のテトラ要素との接触としてもよい。

また、追加部３０４は、許容範囲６０１外と判断された要素群の要素が、複数の要素のうち少なくともいずれか一つの要素から所定距離以内に存在する場合、描画対象要素群Ｕに追加する。具体的には、追加部３０４は、許容範囲６０１外と判断された要素群のテトラ要素と、接触する複数の描画対象要素のうち少なくともいずれか一つのテトラ要素との重心間距離が所定距離以内である場合、許容範囲６０１外と判断された要素群のテトラ要素を描画対象要素群Ｕに追加する。換言すれば追加部３０４は、複数のテトラ要素から所定距離離れている場合は、追加要素にはしない。このように、所定距離による制限を規定しておくことにより、描画対象要素群Ｕの過剰な拡大を抑制することができ、指定領域Ａの自動特定の精度向上を図ることができる。

なお、ここでは、許容範囲６０１外と判断された要素群の要素が、複数の要素のうち少なくともいずれか一つの要素から所定距離以内に存在する場合を条件とするため、複数の要素のいずれに対しても所定距離以内に存在する場合に制限してもよい。

描画部３０５は、隣接要素が追加された描画対象要素群Ｕに応じて、シミュレーションモデルを描画する。具体的には、描画部３０５は、追加部３０４によって追加された追加後の描画対象要素群Ｕに基づいて、臓器モデルを描画する。たとえば、描画部３０５は、描画対象要素群Ｕを描画する。これにより、ユーザは、疾患部位がどのような３次元形状であるかを確認することができる。また、描画部３０５は、心臓モデル１００を描画し、心臓モデル１００のうち描画対象要素群Ｕを強調表示するように描画してもよい。これにより、ユーザは、心臓モデル１００のうちどの程度の領域が疾患部位であるかを確認することができる。また、描画部３０５は、心臓モデル１００のうち描画対象要素群Ｕ以外の残余の要素群を描画してもよい。これにより、ユーザは、疾患部位を除いた心臓モデル１００や、疾患部位の境界を確認することができる。

また、描画部３０５は、追加要素群がなくなった場合、追加後の描画対象要素群Ｕに基づいて、臓器モデルを描画することとしてもよい。これにより、指定領域Ａの自動拡大を停止することができる。一方、追加要素群が存在する場合は、描画部３０５は、まだ描画しない。この場合、上述した特定部３０２は、追加分のテトラ要素群の各要素に隣接する隣接要素群を特定することになる。具体的には、特定部３０２は、あらたに特定された隣接要素群を隣接要素群のデータ構造に追加する。この場合、選択フラグは「０」に設定される。したがって、判断部３０３および追加部３０４は、新たに特定された隣接要素群について判断処理および追加処理を実行することになる。したがって、指定領域Ａの自動拡大をおこなうことができる。

＜描画処理手順＞
図９は、実施の形態にかかる描画装置３００による描画処理手順例を示すフローチャートである。まず、描画装置３００は、心臓モデル１００に対してユーザ操作により領域の指定を受け付け（ステップＳ９０１）、指定領域Ａ上の描画対象要素群Ｕを検出する（ステップＳ９０２）。つぎに、描画装置３００は、描画対象要素群Ｕの物理値についてヒストグラム６００を生成し（ステップＳ９０３）、物理値の許容範囲６０１を検出する（ステップＳ９０４）。

そして、描画装置３００は、描画対象要素群Ｕの各要素の隣接要素群を特定し（ステップＳ９０５）、隣接要素群に未選択要素があるか否かを判断する（ステップＳ９０６）。未選択要素がある場合（ステップＳ９０６：Ｙｅｓ）、拡張処理を実行して（ステップＳ９０７）、ステップＳ９０６に戻る。拡張処理（ステップＳ９０７）の詳細については図１０で後述する。これにより、未選択要素ごとに拡張処理（ステップＳ９０７）が実行されることになる。

また、ステップＳ９０６において、未選択要素がない場合（ステップＳ９０６：Ｎｏ）、描画装置３００は、新規追加分の追加要素群があるか否かを判断する（ステップＳ９０８）。新規追加分の追加要素群がある場合（ステップＳ９０８：Ｙｅｓ）、描画装置３００は、新規追加分の追加要素群の各要素に隣接する隣接要素群を特定し（ステップＳ９０９）、ステップＳ９０６に戻る。これにより、描画装置３００は、新規追加分の追加要素群が得られなくなるまで、隣接要素群が特定されるため、描画対象要素群Ｕを拡張することができる。

また、ステップＳ９０８において、新規追加分の追加要素群がない場合（ステップＳ９０８：Ｎｏ）、描画装置３００は、描画部３０５により描画処理を実行する（ステップＳ９１０）。これにより、図１の（Ｄ）に示したように心臓モデル１００の疾患部位を表示することができる。

図１０は、図９の拡張処理（ステップＳ９０７）の詳細な処理手順例を示すフローチャートである。描画装置３００は、隣接要素群から未選択要素を１つ選択する（ステップＳ１００１）。すなわち、描画装置３００は、選択フラグが「０」のテトラ要素を選択する。選択された隣接要素を「選択隣接要素」と称す。この選択により選択フラグが「０」から「１」に遷移する。描画装置３００は、選択隣接要素の物理値が許容範囲６０１内か否かを判断する（ステップＳ１００２）。許容範囲６０１内である場合（ステップＳ１００２：Ｙｅｓ）、描画装置３００は、選択隣接要素を、追加要素群のデータ構造に追加し（ステップＳ１００３）、ステップＳ９０６に移行する。これにより、選択隣接要素が描画対象要素群Ｕに加わるため、描画対象要素群Ｕを拡張することができる。

一方、ステップＳ１００２において、選択隣接要素の物理値が許容範囲６０１外である場合（ステップＳ１００２：Ｎｏ）、描画装置３００は、選択隣接要素が複数の描画対象要素に接触するか否かを判断する（ステップＳ１００４）。接触する場合（ステップＳ１００４：Ｙｅｓ）、描画装置３００は、選択隣接要素と接触先の複数の指定領域Ａの少なくともいずれか一方のテトラ要素との距離が所定距離ｄ以内であるか否かを判断する（ステップＳ１００５）。

所定距離ｄ以内である場合（ステップＳ１００５：Ｙｅｓ）、ステップＳ１００３に移行して、描画装置３００は、選択隣接要素を追加要素群のデータ構造に追加する（ステップＳ１００３）。これにより、物理値が許容範囲６０１外であっても疾患部位の連続性を担保することができる。

また、ステップＳ１００５において、所定距離ｄ以内でない場合（ステップＳ１００５：Ｎｏ）、描画装置３００は、選択隣接要素を追加要素群のデータ構造に追加せず、ステップＳ９０６に移行する。所定距離ｄ以内という条件があるため、描画対象要素群Ｕの過剰な拡張を抑制することができる。

また、ステップＳ１００４において、選択隣接要素が複数の描画対象要素に接触しない場合（ステップＳ１００４：Ｎｏ）、描画装置３００は、選択隣接要素を追加要素群のデータ構造に追加せず、ステップＳ９０６に移行する。これにより、拡張処理（ステップＳ９０７）を終了する。

このように、本実施の形態によれば、疾患部位の特定の容易化を図ることができるという効果を奏する。具体的には、ユーザは患者や患者の状態ごとに疾患部位を特定する必要がなく、観測したい箇所を指定領域Ａとして指定するだけで、疾患部位にふさわしい描画対象要素を自動的に拡張することができる。

また、選択隣接要素の物理値が許容範囲６０１外である場合であっても、所定の条件下で許容範囲６０１内のテトラ要素とみなすことができるため、疾患部位の連続性を担保することができる。また、この場合、描画対象要素と所定距離以内でない場合は許容範囲６０１内のテトラ要素とみなさないため、描画装置３００は、描画対象要素群Ｕの過剰な拡張を抑制することができる。したがって、描画装置３００は、指定領域Ａの特定精度の向上を図ることができる。

また、描画対象要素への追加分のテトラ要素がなくなれば、拡張処理が終了する。このため、ユーザは最初に指定領域Ａを指定するだけで、どのくらいの領域までが疾患部位であるか、といった判断をおこなうことなく、最終的に拡張された描画対象要素群Ｕを疾患部位として得ることができる。したがって、疾患部位特定の利便性の向上を図ることができる。

１００ 心臓モデル
３００ 描画装置
３０１ 検出部
３０２ 特定部
３０３ 判断部
３０４ 追加部
３０５ 描画部
４００，７００，８００ データ構造
６００ ヒストグラム
６０１ 許容範囲
Ｕ 描画対象要素群

Claims (8)

1. シミュレーションの結果を描画する描画プログラムにおいて、
   コンピュータに、
   シミュレーション対象の位置毎に物理値をそれぞれ有する要素の集合であるシミュレーションモデルにおいて指定された領域の描画対象である描画対象要素群が有する物理値の範囲を検出させ、
   前記描画対象要素群に含まれる各描画対象要素に隣接する要素の集合である隣接要素群を前記シミュレーションモデルから特定させ、
   特定された前記隣接要素群に含まれる各隣接要素の物理値が前記範囲に含まれるかを判断させ、
   物理値が前記範囲に含まれると判断された隣接要素を前記描画対象要素群に追加させ、
   隣接要素が追加された描画対象要素群に応じて、前記シミュレーションモデルを描画させることを特徴とする描画プログラム。
2. 前記隣接要素群の特定では、
   前記描画対象要素群への追加分の要素群の各要素に隣接する要素の集合を前記隣接要素群としてあらたに前記シミュレーションモデルから特定することを特徴とする請求項１に記載の描画プログラム。
3. 前記隣接要素の前記描画対象要素群への追加では、
   前記隣接要素群のうち前記範囲外と判断された要素群の要素が、前記描画対象要素群の中の複数の要素に接する場合、前記描画対象要素群に追加することを特徴とする請求項２に記載の描画プログラム。
4. 前記隣接要素の前記描画対象要素群への追加では、
   前記範囲外と判断された要素群の要素が、前記複数の要素のうち少なくともいずれか一つの要素から所定距離以内に存在する場合、前記描画対象要素群に追加することを特徴とする請求項３に記載の描画プログラム。
5. 前記シミュレーションモデルの描画では、
   前記追加分の要素群がなくなった場合、前記追加後の描画対象要素群に基づいて、臓器モデルを描画することを特徴とする請求項１〜４のいずれか一つに記載の描画プログラム。
6. 前記シミュレーションモデルの描画では、
   前記追加後の描画対象要素群と前記シミュレーションモデルのうち前記追加後の描画対象要素群を除く残余の要素群とのうち、少なくともいずれか一方の要素群を描画することを特徴とする請求項１〜５のいずれか一つに記載の描画プログラム。
7. シミュレーションの結果を描画する描画方法において、
   コンピュータが、
   シミュレーション対象の位置毎に物理値をそれぞれ有する要素の集合であるシミュレーションモデルにおいて指定された領域の描画対象である描画対象要素群が有する物理値の範囲を検出し、
   前記描画対象要素群に含まれる各描画対象要素に隣接する要素の集合である隣接要素群を前記シミュレーションモデルから特定し、
   特定された前記隣接要素群に含まれる各隣接要素の物理値が前記範囲に含まれるかを判断し、
   物理値が前記範囲に含まれると判断された隣接要素を前記描画対象要素群に追加し、
   隣接要素が追加された描画対象要素群に応じて、前記シミュレーションモデルを描画することを特徴とする描画方法。
8. シミュレーションの結果を描画する描画装置において、
   シミュレーション対象の位置毎に物理値をそれぞれ有する要素の集合であるシミュレーションモデルにおいて指定された領域の描画対象である描画対象要素群が有する物理値の範囲を検出する検出部と、
   前記描画対象要素群に含まれる各描画対象要素に隣接する要素の集合である隣接要素群を前記シミュレーションモデルから特定する特定部と、
   特定された前記隣接要素群に含まれる各隣接要素の物理値が前記範囲に含まれるかを判断する判断部と、
   物理値が前記範囲に含まれると判断された隣接要素を前記描画対象要素群に追加する追加部と、
   隣接要素が追加された描画対象要素群に応じて、前記シミュレーションモデルを描画する描画部と、
   を有することを特徴とする描画装置。

Images