(57)【要約】【課題】本発明の課題は、シミュレーション結果の可視化処理に必要十分なデータを抽出することを目的とする。【解決手段】上記課題は、現象部分に係る線分の近傍の断面毎に、断面形状を示す断面形状情報と、該断面のテクスチャを表わすテクスチャ情報とを記憶した記憶部から、視点方向に対応する断面形状情報とテクスチャ情報とを選択する選択部と、選択した前記断面形状情報と前記テクスチャ情報とを用いて、前記線分を３次元画像で可視化する可視化データを生成し、該可視化データにより該線分を表示部に表示させる表示部と、を有することを特徴とする情報処理装置により達成される。【選択図】図５
2018005671.tif 000002
【特許請求の範囲】
【請求項１】
現象部分に係る線分の近傍の断面毎に、断面形状を示す断面形状情報と、該断面のテクスチャを表わすテクスチャ情報とを記憶した記憶部から、視点方向に対応する断面形状情報とテクスチャ情報とを選択する選択部と、
選択した前記断面形状情報と前記テクスチャ情報とを用いて、前記線分を３次元画像で可視化する可視化データを生成し、該可視化データにより該線分を表示部に表示させる表示部と、
を有することを特徴とする情報処理装置。
【請求項２】
前記線分の中心と、該線分の一端を含む第１平面と、該線分の他端を含む第２平面とを、該第１平面と該第２平面とが交わり、かつ、該中心を含む直線で結合して、前記断面形状情報を作成する断面化処理部を有することを特徴とする請求項１記載の情報処理装置。
【請求項３】
前記近傍の断面を包含する３次元の矩形領域に対して、垂直な各方向で該矩形領域を等分する切断面を設定し、該切断面の頂点を用いて曲面を形成することで、前記断面形状情報を作成する断面化処理部を有することを特徴とする請求項１記載の情報処理装置。
【請求項４】
現象部分に係る線分の近傍の断面毎に、断面形状を示す断面形状情報と、該断面のテクスチャを表わすテクスチャ情報とを記憶した記憶部から、視点方向に対応する断面形状情報とテクスチャ情報とを選択し、
選択した前記断面形状情報と前記テクスチャ情報とを用いて、前記線分を３次元画像で可視化する可視化データを生成し、該可視化データにより該線分を表示部に表示させる、
処理をコンピュータが行うシミュレータ結果表示方法。
【請求項５】
現象部分に係る線分の近傍の断面毎に、断面形状を示す断面形状情報と、該断面のテクスチャを表わすテクスチャ情報とを記憶した記憶部から、視点方向に対応する断面形状情報とテクスチャ情報とを選択し、
選択した前記断面形状情報と前記テクスチャ情報とを用いて、前記線分を３次元画像で可視化する可視化データを生成し、該可視化データにより該線分を表示部に表示させる、
処理をコンピュータに行わせるシミュレータ結果表示プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、情報処理装置、シミュレータ結果表示方法、及びシミュレータ結果表示プログラムに関する。
【背景技術】
【０００２】
近年、計算機上で実施する数値シミュレーションにより、現象の再現や予測が行われている。数値シミュレーション後、数値データを直接確認するのでは得られない可視化処理が行われる。可視化にはコンピュータグラフィックスの技術を利用した３次元可視化、動画化による時間変化の可視化が一般的に行われている。
【０００３】
一方、スーパーコンピュータを用いた大規模計算処理によるシミュレーションの結果は、ペタバイトオーダのデータサイズを超える超巨大データとなって出力される。これらのデータに対し、可視化方法の検討が数多く行われている。
【０００４】
大規模データからスカラー場を作成し、直接ボリュームレンダリングという手法で可視化する技術、ある瞬間における各点の流体分子の速度ベクトルに沿ってできる流線の節点間ごとに四辺形を生成し、その四辺形上のポリゴンにテクスチャマッピングする技術等が知られている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】再公表ＷＯ２０１３／０３８５４８号公報
【特許文献２】特開平１１−１０４０７２号公報
【特許文献３】特開平１１−０６６２７７号公報
【非特許文献】
【０００６】
【非特許文献１】松岡、荒木、「大規模シミュレーションデータ可視化における研究の動向」、JAMSTEC  Rep．Res．Dev．、Volume  13、September  2011、35_63
【非特許文献２】“名古屋大学大型計算機センターライブラリープログラム使用法説明書  DAFM：3次元非構造格子自動生成プログラム”、［online］、［平成２８年４月２７日検索］、インターネット＜URL：http：//iij．sinsakijun．com/information/newtech．html＞http：//www2．itc．nagoya‐u．ac．jp/sys_riyou/lib/dafm．pdf
【非特許文献３】Motoi  Kinishi  et  al．，“Analysis  of  Centerline  Extraction  in  Three‐Dimensional  Scale  Space‐Extracting  Centerline  of  Vessels  in  Hepatic  Artery”，2009  Fifth  International  Conference  on  Natural  Computation，Volume  6，14‐16  Aug．2009
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、俯瞰的な可視化を実現する上で、データサイズの削減は大きな課題である。そのため、一般的には、格子情報の間引きが行われるが、大規模シミュレーションの解析精度のまま表示することができないといった問題がある。
【０００８】
したがって、１つの側面では、本発明は、シミュレーション結果の可視化処理に必要十分なデータを抽出することを目的とする。
【課題を解決するための手段】
【０００９】
一態様によれば、現象部分に係る線分の近傍の断面毎に、断面形状を示す断面形状情報と、該断面のテクスチャを表わすテクスチャ情報とを記憶した記憶部から、視点方向に対応する断面形状情報とテクスチャ情報とを選択する選択部と、選択した前記断面形状情報と前記テクスチャ情報とを用いて、前記線分を３次元画像で可視化する可視化データを生成し、該可視化データにより該線分を表示部に表示させる表示部と、を有することを特徴とする情報処理装置が提供される。
【００１０】
また、上記課題を解決するための手段として、シミュレータ結果表示方法、及び、シミュレータ結果表示プログラムとすることもできる。
【発明の効果】
【００１１】
シミュレーション結果の可視化処理に必要十分なデータを抽出することができる。
【図面の簡単な説明】
【００１２】
【図１】特徴抽出から可視化までの処理例を説明するための図である。
【図２】３次元可視化の例を示す図である。
【図３】カメラ位置に応じた平面選択による現象領域の表示変化例を示す図である。
【図４】情報処理装置のハードウェア構成を示す図である。
【図５】本実施例における情報処理装置の機能構成例を示す図である。
【図６】計算空間から抽出された領域例を示す図である。
【図７】本実施例における断面化方法の概要を説明するための図である。
【図８】本実施例における３次元可視化処理を説明するためのフローチャート図である。
【図９】構造格子及び非構造格子の例を示す図である。
【図１０】断面化方法１及び断面化方法２による断面化例を示す図である。
【図１１】断面化方法１の処理例を説明するための図である。
【図１２】断面化方法２の処理例を説明するための図である。
【図１３】処理結果画像の例を示す図である。
【図１４】図１３の処理結果画像に対して一部断面の設定を行った例を示す図である。
【発明を実施するための形態】
【００１３】
以下、本発明の実施の形態を図面に基づいて説明する。図１は、特徴抽出から可視化までの処理例を説明するための図である。図１では、数値シミュレーションにより得られた大規模データから作成されたスカラー場のデータ２を用いて可視化する処理の一例を説明する。
【００１４】
先ず、スカラー場データ２から現象を表わす現象領域１の中心線３を取得して、中心線３上に節点Ｉ０、Ｉ１、・・・を与え、各節点Ｉｋ（ｋ＝０、１、・・・）におけるベクトルｖ×（Ｉ０−Ｉ１）を求める（１ａ）。節点Ｉｋ毎に求めたベクトルを用いて捻じれ情報１ｂを生成し、生成した捻じれ情報１ｂを用いて、曲面１ｃを可視化する。このような３次元可視化の例を図２に示す。
【００１５】
図２は、３次元可視化の例を示す図である。図２では、シミュレーションにより得られた結果を例として示している。図２（Ａ）に、キャビティ流れ方向を示している。図２（Ｂ）では、数値シミュレーションにより得られたスカラー場データ２から抽出した現象領域１の例を示している。現象領域１は、例えば、図１の曲面１ｃのように表される。図２（Ｃ）では、現象領域１のカメラ位置に応じた３次元可視化の例を示している
図３は、カメラ位置に応じた平面選択による現象領域の表示変化例を示す図である。カメラ位置と、各節点Ｉｋのベクトルとに基づいて、中心線３上の隣接する２つの節点Ｉｋを含む平面を形成し、隣接する側の辺を平面同士で共有することにより形成される。
【００１６】
図３（Ａ）は、カメラ位置ＰＡに応じて選択された平面で描画された現象領域１の例を示している。図３（Ａ）で示す注視点の移動を表わすカメラ移動方向２ｄに従って、図３（Ｂ）〜図３（Ｅ）に、変化する現象領域１の表示例を示す。
【００１７】
図３（Ｂ）は、カメラ位置ＰＢに応じて選択された平面で描画された現象領域１の例を示している。図３（Ｃ）は、カメラ位置ＰＣに応じて選択された平面で描画された現象領域１の例を示している。図３（Ｄ）は、カメラ位置ＰＤに応じて選択された平面で描画された現象領域１の例を示している。図３（Ｅ）は、カメラ位置ＰＥに応じて選択された平面で描画された現象領域１の例を示している。
【００１８】
可視化する現象として、形状の存在しない圧力、気流等が考えられるが、本実施例は、これら現象に限定されず、種々の現象領域１を３次元に可視化可能である。
【００１９】
本実施例における３次元可視化は、図４に示すような情報処理装置１００によって行われる。図４は、情報処理装置のハードウェア構成を示す図である。図４において、情報処理装置１００は、コンピュータによって制御される情報処理装置であって、ＣＰＵ（Central  Processing  Unit）１１と、主記憶装置１２と、補助記憶装置１３と、入力装置１４と、表示装置１５と、通信Ｉ／Ｆ（インターフェース）１７と、ドライブ装置１８とを有し、バスＢに接続される。
【００２０】
ＣＰＵ１１は、主記憶装置１２に格納されたプログラムに従って情報処理装置１００を制御するプロセッサに相当する。主記憶装置１２には、ＲＡＭ（Random  Access  Memory）、ＲＯＭ（Read  Only  Memory）等が用いられ、ＣＰＵ１１にて実行されるプログラム、ＣＰＵ１１での処理に必要なデータ、ＣＰＵ１１での処理にて得られたデータ等を記憶又は一時保存する。
【００２１】
補助記憶装置１３には、ＨＤＤ（Hard  Disk  Drive）等が用いられ、各種処理を実行するためのプログラム等のデータを格納する。補助記憶装置１３に格納されているプログラムの一部が主記憶装置１２にロードされ、ＣＰＵ１１に実行されることによって、各種処理が実現される。
【００２２】
入力装置１４は、マウス、キーボード等を有し、ユーザが情報処理装置１００による処理に必要な各種情報を入力するために用いられる。表示装置１５は、ＣＰＵ１１の制御のもとに必要な各種情報を表示する。入力装置１４と表示装置１５とは、一体化したタッチパネル等によるユーザインタフェースであってもよい。通信Ｉ／Ｆ１７は、有線又は無線などのネットワークを通じて通信を行う。通信Ｉ／Ｆ１７による通信は無線又は有線に限定されるものではない。  
情報処理装置１００によって行われる処理を実現するプログラムは、例えば、ＣＤ−ＲＯＭ（Compact  Disc  Read‐Only  Memory）等の記憶媒体１９によって情報処理装置１００に提供される。
【００２３】
ドライブ装置１８は、ドライブ装置１８にセットされた記憶媒体１９（例えば、ＣＤ−ＲＯＭ等）と情報処理装置１００とのインターフェースを行う。
【００２４】
また、記憶媒体１９に、後述される本実施の形態に係る種々の処理を実現するプログラムを格納し、この記憶媒体１９に格納されたプログラムは、ドライブ装置１８を介して情報処理装置１００にインストールされる。インストールされたプログラムは、情報処理装置１００により実行可能となる。
【００２５】
尚、プログラムを格納する記憶媒体１９はＣＤ−ＲＯＭに限定されず、コンピュータが読み取り可能な、構造（structure）を有する１つ以上の非一時的（non‐transitory）な、有形（tangible）な媒体であればよい。コンピュータ読取可能な記憶媒体として、ＣＤ−ＲＯＭの他に、ＤＶＤディスク、ＵＳＢメモリ等の可搬型記録媒体、フラッシュメモリ等の半導体メモリであっても良い。
【００２６】
図５は、本実施例における情報処理装置の機能構成例を示す図である。図５において、情報処理装置１００は、空間内挿部４１と、現象領域抽出部４２と、中心線抽出部４３と、断面化処理部４４と、テクスチャ生成部４５と、レンダリング部４６とを有する。空間内挿部４１と、現象領域抽出部４２と、中心線抽出部４３と、断面化処理部４４と、テクスチャ生成部４５と、レンダリング部４６とは、情報処理装置１００にインストールされたプログラムが、情報処理装置１００のＣＰＵ１１に実行させる処理により実現される。
【００２７】
また、記憶部１３０には、ベクトル場データ５１、スカラー場のλ２分布データ５２、断面形状情報５３、テクスチャ情報５４、可視化データ５５、カメラ位置情報５９等が記憶される。
【００２８】
空間内挿部４１は、ベクトル場データ５１に基づく計算空間７（図６）全体を矩形領域により包含し、構造格子状に内挿処理を行う。ベクトル場データ５１は、構造格子又は非構造格子のいずれのデータであってもよく、また、混在していてもよい。計算空間７（図６）は、ユークリッド空間とする。
【００２９】
ベクトル場データ５１が構造格子のデータである場合、内挿処理は行われない。ベクトル場データ５１のうち、非構造格子のデータに対して、空間内挿部４１による内挿処理が行われる。
【００３０】
現象領域抽出部４２は、Jing等のλ２法を用いて、現象領域（例えば、渦領域又は血管領域等）１を抽出する。λ２法により得られたスカラー場によるλ２分布データ５２を用いて、現象領域１の境界４（図６）が特定される。現象領域１の抽出には、渦領域と渦領域でない領域とを分けること、血管の領域と血管でない領域とを分けることを含む。スカラー場によるλ２分布データ５２が記憶部１３０に出力される。
【００３１】
中心線抽出部４３は、非特許文献３に基づく既存技術により、現象領域抽出部４２によって抽出された領域の中心線３を抽出する。
【００３２】
断面化処理部４４は、中心線３に基づいて抽出した領域に対する断面化処理を行う。断面形状情報５３が記憶部１３０に記憶される。断面形状情報５３は、複数の断面の形状に係る情報を含む。
【００３３】
本実施例における断面化処理部４４では、ユーザの選択に応じて、下記の断面化方法１又は断面化方法２を行う。
【００３４】
・断面化方法１
中心線３の端点毎に、その端点と中心線３の領域中心５（図６）とを含む複数の平面を生成して、断面化処理を行う。
【００３５】
・断面化方法２
抽出された領域の最大最小の空間情報を設定し、空間情報から得られる中点情報を用いた曲面を生成して、断面化処理を行う。空間情報は、立方体又は直方体等の３次元の矩形領域を表わす。以下、単に「矩形領域」という。
【００３６】
テクスチャ生成部４５は、断面形状情報５３を用いて、各断面形状に応じたテクスチャを生成する。テクスチャ情報５４が記憶部１３０に出力される。テクスチャ情報５４は、三次元空間の所定位置に配置された光源に対する帯状面の反射光を表わすテクスチャ等の情報である。
【００３７】
レンダリング部４６は、記憶部１３０から、カメラ位置情報５９に応じた断面形状情報５３とテクスチャ情報５４とを選択する選択部と、可視化データ５５を記憶部１３０に生成し、表示装置１５に表示させる表示部とを有する。
【００３８】
図６は、計算空間から抽出された領域例を示す図である。図６において、内挿処理後の計算空間７において、抽出された現象領域１の例を示している。抽出された各現象領域１から中心線３が抽出され、更に、抽出された中心線３の中心を特定することで領域中心５が定められる。図６では、抽出された現象領域１は、境界４で示される。以下の説明においても同様である。
【００３９】
抽出された領域毎に、断面化方法１又は断面化方法２が行われる。断面化方法１及び断面化方法２について説明する。
【００４０】
図７は、本実施例における断面化方法の概要を説明するための図である。図７（Ａ）は、中心線３に基づく曲率判定により平面指定を行う断面化方法１を示している。断面化方法１では、領域から抽出された中心線３の中心を求めて、領域中心５を設定する。
【００４１】
領域中心５を含む一辺を共有する、中心線３の一方の端点と領域中心５とを含む平面と、中心線３の他方の端点と領域中心５とを含む平面の２つの平面のペアを設定することで断面化処理が行われる。異なる角度θで、領域中心５を含む一辺を共有する複数の平面のペアが作成されればよい。この例ではθ＝９０°の例を示しているが、任意のθを設定可能である。
【００４２】
図７（Ｂ）は、抽出された領域の最大最小の空間情報を設定し、空間情報から得られる中点情報を用いた曲面を生成する断面化方法２を示している。断面化方法１と同様に、現象領域１の領域中心５を設定する。また、境界４に基づいて予測した現象領域１の３次元形状に対して最大最小の矩形領域６を決定する。
【００４３】
矩形領域６の各辺の中心点を用いて曲線８を定め、定めた複数の曲線８により曲面を生成する。
【００４４】
図８は、本実施例における３次元の可視化処理を説明するためのフローチャート図である。図８において、空間内挿部４１は、計算空間７における現象領域１に係るベクトル場データ５１を用いて、空間内挿処理を行う（ステップＳ７１）。非構造格子に対して空間内処理が実施される。
【００４５】
現象領域抽出部４２は、計算空間７のλ２値分布を算出してλ２分布データ５２を生成し、現象領域１を抽出する（ステップＳ７２）。現象領域１の抽出によって、少なくとも境界４が抽出される。そして、中心線抽出部４３は、計算空間７からλ２分布データ５２に基づいて、境界４と中心線３とを抽出する（ステップＳ７３）。
【００４６】
その後、断面化処理部４４は、ユーザが指定した断面化方法１又は２により、中心線抽出部４３が抽出した境界４と中心線３に基づいて断面化処理を行う（ステップＳ７４）。断面形状情報５３が出力される。
【００４７】
そして、テクスチャ生成部４５は、断面形状情報５３に基づいて、境界４で特定される各領域に対して断面形状毎にテクスチャを生成し、生成したテクスチャを表わすテクスチャ情報５４を生成する（ステップＳ７５）。
【００４８】
レンダリング部４６は、断面形状情報５３及びテクスチャ情報５４を用いて、カメラ位置情報５９に基づいて、可視化データ５５を作成する（ステップＳ７６）。可視化データ５５により、表示装置１５に３次元画像又は動画が表示される。
【００４９】
次に、空間内挿部４１による空間内挿処理について説明する。空間内挿部４１は、構造格子、非構造格子、又はそれらが混在したベクトル場データ５１を入力する。
【００５０】
図９は、構造格子及び非構造格子の例を示す図である。図９（Ａ）は、構造格子の例を示す。構造格子は、六面体で表される。図９（Ｂ）は、非構造格子の例を示す。非構造格子は、四面体（三角錐）又は五面体（四角錐又は三角錐）で表される。
【００５１】
空間内挿部４１は、計算空間７の全体を矩形領域により包含し、構造格子状に内挿処理を行う。内挿処理は、以下の手順で行えばよい。
【００５２】
１）非構造格子において、辺の長さを基準に自動的にヒストグラム化する。つまり、
（min，max）＝（1．0e‐3，1．0e+2）であれば、対数で１０分割する。
【００５３】
２）最初のヒストグラム範囲のカウント値（例えばａ）と、最大のカウント値（例えばｂ）を数１に導入し、代表値Ｌを得る。この値に基づき、空間を均等に分割する。
【００５４】
【数１】
2018005671.tif 000003
【００５５】
入力されるデータ構造は構造格子及び非構造格子の両方であっても良い。
【００５６】
非構造格子のデータの内挿処理では、非特許文献２の四面体要素を用いて三次元構造格子を生成する手法を用い、ユーザが指定する分割数に内挿する。現象領域１の抽出を、内挿精度にあまり依存せず行える。
【００５７】
現象領域抽出部４２は、以下の手順で現象領域１とそうでない領域とを分ける。Jing等によるλ２法を用いた例で説明するが、現象領域１の抽出は、λ２法に限定されない。λ２法を用いた場合、出力される値はλ２値（スカラー値）となる。
【００５８】
１）構造格子の節点上に速度ベクトル（u，v，w）を配置する。非構造格子の場合は構造格子上の節点に内挿する。
【００５９】
２）（ui，vi，wi）から節点ｉ上の勾配を計算し、
【００６０】
【数２】
2018005671.tif 000004
行列Ｊを作成する。
【００６１】
３）そして、
Ａ＝Ｓ２＋Ｑ２
を計算する。ここで、Ｓは、
【００６２】
【数３】
2018005671.tif 000005
また、Ｑは、
【００６３】
【数４】
2018005671.tif 000006
である。
【００６４】
４）行列Ａの固有値を計算し、
５）λ２を選択する。つまり、λ１、λ２、及びλ３から２番目の値を選択する。
【００６５】
中心線抽出処理部４３による中心線抽出処理は、既存技術で行う。一例として、非特許文献３の手法を用いれば良い。中心線３が抽出される。
【００６６】
次に、本実施例に係る断面化処理とテクスチャ生成処理について説明する。断面化処理部４４による断面化処理は、上述したように、断面化方法１又は断面化方法２により行う。
【００６７】
図１０は、断面化方法１及び断面化方法２による断面化例を示す図である。図１０（Ａ）では、現象領域１の領域中心５を含む一辺を共有する２つの平面９ａによって１つの断面を作成する。領域中心５を含む一辺を共有する２つの平面９ａのペアによって一断面が表される。
【００６８】
図１０（Ｂ）では、最大最小の矩形領域６の中点情報を用いて、生成した２つの曲面９ｂによって１つの断面を作成する。矩形領域６の長さ方向に平行な４面の各々において、幅方向で対向する２辺の各中点を含む直線を共有する辺とした、上下２つの曲面９ｂによって、１つの断面が表わされる。
【００６９】
図１１は、断面化方法１の処理例を説明するための図である。図１１（ａ）は、現象領域１の領域中心５を設定した状態を示している。
【００７０】
図１１（ｂ）において、中心線３の両端５ａ及び５ｂの３次元座標（xa，ya，za）及び（xb，yb，zb）を取得する。領域中心５と端５ａとを含む四辺形９ａ１１と、領域中心５と端５ｂとを含む四辺形９ａ１２とを定める。四辺形９ａ１１と四辺形９ａ１２とによって、境界４を含む必要最小限の領域を設定する。
【００７１】
四辺形９ａ１１と四辺形９ａ１２とに基づいて、領域中心５を含む直線を一辺として共有する２つの平面９ａのペアを作成する。角度θを変えて複数の平面９ａのペアを作成することで、ｎ断面を作成する。
【００７２】
単純には、中心線３の各軸方向の最大値及び最小値を求め、最も距離の長い軸を選択し、その距離の半分を求めて中心位置とする。その中心位置と端５ａを共有する平面（四辺形９ａ１１）、或いは、端５ｂと領域中心５（四辺形９ａ１２）を共有する平面を、２点を結ぶ軸で角度θ毎に回転させることで、ｎ断面を作成する。
【００７３】
図１１（ｃ）において、領域中心５と端５ａを含む四辺形９ａ１１と、端５ｂと領域中心５とを含む四辺形９ａ１２とを生成する。
【００７４】
図１１（ｄ）において、領域中心５を通り、四辺形９ａ１１を含む平面と、四辺形９ａ１２を含む平面とが交わる直線５γ１を求めて、２つの平面を結合したペアの平面９ａを作成する。
【００７５】
図１１（ｅ）では、領域中心５を軸として、角度θ回転させ、図１１（ｃ）と同様に、領域中心５と端５ａとを含む四辺形９ａ２１及び端５ｂと領域中心５とを含む四辺形９ａ２２を生成する。
【００７６】
図１１（ｆ）において、領域中心５を通り、四辺形９ａ２１を含む平面と、四辺形９ａ２２を含む平面とが交わる直線５γ２を求めて、２つの平面を結合したペアの平面９ａを作成する。
【００７７】
更に、図１１（ｇ）では、直線５γ１及び直線５γ２を含む平面において、直線５γ２に対して更に角度θ回転させ、領域中心５と端５ａとを含む四辺形９ａ３１及び端５ｂと領域中心５とを含む四辺形９ａ３２を生成する。
【００７８】
図１１（ｈ）において、領域中心５を通り、四辺形９ａ３１を含む平面と、四辺形９ａ３２を含む平面とが交わる直線５γ３を求めて、２つの平面を結合したペアの平面９ａを作成する。
【００７９】
よって、角度θ＝６０°とした場合、図１１（ｉ）に示すように、３つのペアの平面９ａが定められる。つまり、６方向の断面を定めたことになる。
【００８０】
図１２は、断面化方法２の処理例を説明するための図である。図１２（ａ）において、現象領域１の空間的な最大最小の矩形領域６を設定した状態を示している。この矩形領域６の情報を使って面構成を行う。
【００８１】
図１２（ｂ）において、先ず、不定形の現象領域１に対して、ｘ、ｙ、及びｚ軸方向の最小値及び最大値をそれぞれ探索し、矩形領域６の角に相当する点ｅ１（xe1，ye1，ze1）及び点ｅ２（xe2，ye2，ze2）を設定する。
【００８２】
点ｅ１及び点ｅ２間の距離で、ｘ、ｙ、及びｚの値が最大の座標軸を選択する。ｙが最も長い場合、（ye2‐ye1）/２の位置で、矩形領域６を１／２にする平面Ｄ１を設定する。
【００８３】
図１２（ｃ）において、また、ｘ軸に対しても中間の位置（xe2‐xe1）/２の位置で、矩形領域６を１／２にする平面Ｄ２を設定する。平面Ｄ２が矩形領域６の辺と交わる点ｐ１（xp1，yp1，zp1）と、平面Ｄ１が矩形領域６の辺とが交わる点ｑ１（xq1，yq1，zq1）とを繋ぐ二次曲線を設定する。また、平面Ｄ２が矩形領域６の辺と交わる点ｐ２（xp2，yp2，zp2）と、平面Ｄ１が矩形領域６の辺とが交わる点ｑ２（xq2，yq2，zq2）とを繋ぐ二次曲線を設定する。点ｐ１と点ｑ１とを繋ぐ二次曲線と、点ｐ２と点ｑ２とを繋ぐ二次曲線と、点ｐ１と点ｐ２とを繋ぐ直線と、点ｑ１と点ｑ２とを繋ぐ直線とにより曲面９ｂが作成される。
【００８４】
図１２（ｄ）において、同様に、点ｐ１と点ｑ３とを繋ぐ二次曲線を設定し、点ｐ２と点ｑ４とを繋ぐ二次曲線を設定する。点ｐ１と点ｑ３とを繋ぐ二次曲線と、点ｐ２と点ｑ４とを繋ぐ二次曲線と、点ｐ１と点ｐ２とを繋ぐ直線と、点ｑ３と点ｑ４とを繋ぐ直線とにより曲面９ｂが作成される。図１２（ｃ）で作成した曲面９ｂの点ｐ１と点ｐ２とを繋ぐ直線は、共有される。
【００８５】
また、点ｐ３と点ｑ１とを繋ぐ二次曲線を設定し、点ｐ４と点ｑ２とを繋ぐ二次曲線を設定する。点ｐ３と点ｑ１とを繋ぐ二次曲線と、点ｐ４と点ｑ２とを繋ぐ二次曲線と、点ｐ３と点ｐ４とを繋ぐ直線と、点ｐ３と点ｐ４とを繋ぐ直線とにより他の曲面９ｂが作成される。
【００８６】
更に、点ｐ３と点ｑ３とを繋ぐ二次曲線を設定し、点ｐ４と点ｑ４とを繋ぐ二次曲線を設定する。曲面９ｂが作成される。
【００８７】
このように作成された４つの曲面９ｂは、隣接する曲面９ｂとの間で中点情報で得られる直線を共有する。
【００８８】
図１２（ｅ）では、ｚ軸に対して、中間の位置（ze2‐Ze1）/２の位置で、矩形領域６を１／２にする平面Ｄ３を設定する。平面Ｄ３が矩形領域６の辺と交わる点ｒ１と、平面Ｄ１が矩形領域６の辺とが交わる点ｑ３とを繋ぐ二次曲線を設定する。そして、平面Ｄ３が矩形領域６の辺と交わる点ｒ２と、平面Ｄ１が矩形領域６の辺とが交わる点ｑ１とを繋ぐ二次曲線を設定する。曲面９ｂが作成される。
【００８９】
同様に、点ｒ１と点ｑ４とを繋ぐ二次曲線を設定し、点ｒ２と点ｑ３とを繋ぐ二次曲線を設定して曲面９ｂを作成する。
【００９０】
また、点ｒ３と点ｑ３とを繋ぐ二次曲線を設定し、点ｒ４と点ｑ１とを繋ぐ二次曲線を設定して曲面９ｂを作成する。更に、点ｒ３と点ｑ４とを繋ぐ二次曲線を設定し、点ｒ４と点ｑ２とを繋ぐ二次曲線を設定して曲面９ｂを作成する。断面化方法２では、断面を２つの曲面９ｂのペアで表す。
【００９１】
本実施例による可視化処理での表示例を図１３及び図１４で示す。図１３は、処理結果画像の例を示す図である。図１４は、図１３の処理結果画像に対して一部断面の設定を行った例を示す図である。
【００９２】
図１４では、断面化方法１での断面設定の例を示している。上述したような本実施例による可視化処理では、データ量を削減可能であることが分かっている。
【００９３】
入力は、
・シミュレーション結果
節点の座標、要素情報、要素上の物理値
・カラーマップデータ
・ユーザ指定のパラメータ
断面化方法１又は断面化方法２
である。
【００９４】
本実施例において、可視化処理に必要となるデータ量は、
元の３次元ポリゴン数をｎ、各断面での使用ポリゴンが２、１オブジェクト当たりの断面数（平均）が５、オブジェクト数がｍとすると、
ｎ＞＞ｍ×２×５  ＋  テクスチャデータ量
で表される。本実施例の可視化処理では、現象領域１を含む最小の領域に相当する平面９ａ又は曲面９ｂの領域に相当するデータを抽出することで、現象の三次元の可視化処理に利用するデータ量を少なくすることができる。本実施例の可視化処理を用いない既存の技術と比べて最大１／１００までデータ量を削減可能である。
【００９５】
従って、本実施例では、構造・非構造格子の３次元大規模データを用いたシミュレーション結果の可視化処理において、データリダクションを可能とする。
【００９６】
本発明は、具体的に開示された実施例に限定されるものではなく、特許請求の範囲から逸脱することなく、主々の変形や変更が可能である。
【００９７】
以上の実施例を含む実施形態に関し、更に以下の付記を開示する。
（付記１）
現象部分に係る線分の近傍の断面毎に、断面形状を示す断面形状情報と、該断面のテクスチャを表わすテクスチャ情報とを記憶した記憶部から、視点方向に対応する断面形状情報とテクスチャ情報とを選択する選択部と、
選択した前記断面形状情報と前記テクスチャ情報とを用いて、前記線分を３次元画像で可視化する可視化データを生成し、該可視化データにより該線分を表示部に表示させる表示部と、
を有することを特徴とする情報処理装置。
（付記２）
前記線分の中心と、該線分の一端を含む第１平面と、該線分の他端を含む第２平面とを、該第１平面と該第２平面とが交わり、かつ、該中心を含む直線で結合して、前記断面形状情報を作成する断面化処理部を有することを特徴とする付記１記載の情報処理装置。
（付記３）
前記断面化処理部は、
前記第１平面又は前記第２平面において、前記中心を通る直線を軸として所定角度毎に回転させて、角度毎の前記断面形状情報を作成する
ことを特徴とする付記２記載の情報処理装置。
（付記４）
前記近傍の断面を包含する３次元の矩形領域に対して、垂直な各方向で該矩形領域を等分する切断面を設定し、該切断面の頂点を用いて曲面を形成することで、前記断面形状情報を作成する断面化処理部を有することを特徴とする付記１記載の情報処理装置。
（付記５）
前記３次元の矩形領域の各面において、前記切断面により等分された領域毎に前記曲面が作成され、
前記面を等分する直線は、該曲面の一部として形成され、もう一方の領域に作成された曲面と共有されることを特徴とする付記４記載の情報処理装置。
（付記６）
現象部分に係る線分の近傍の断面毎に、断面形状を示す断面形状情報と、該断面のテクスチャを表わすテクスチャ情報とを記憶した記憶部から、視点方向に対応する断面形状情報とテクスチャ情報とを選択し、
選択した前記断面形状情報と前記テクスチャ情報とを用いて、前記線分を３次元画像で可視化する可視化データを生成し、該可視化データにより該線分を表示部に表示させる、
処理をコンピュータが行うシミュレータ結果表示方法。
（付記７）
現象部分に係る線分の近傍の断面毎に、断面形状を示す断面形状情報と、該断面のテクスチャを表わすテクスチャ情報とを記憶した記憶部から、視点方向に対応する断面形状情報とテクスチャ情報とを選択し、
選択した前記断面形状情報と前記テクスチャ情報とを用いて、前記線分を３次元画像で可視化する可視化データを生成し、該可視化データにより該線分を表示部に表示させる、
処理をコンピュータに行わせるシミュレータ結果表示プログラム。
【符号の説明】
【００９８】
１            現象領域
３            中心線
４            境界
５            領域中心
６            矩形領域
７            計算空間
８            曲線
９ａ          平面
９ｂ          曲面
４１          空間内挿部
４２          現象領域抽出部
４３          中心線抽出部
４４          断面化処理部
４５          テクスチャ生成部
４６          レンダリング部
５１          ベクトル場データ
５２          λ２分布データ
５３          断面形状情報
５４          テクスチャ情報
５５          可視化データ
５９          カメラ位置情報
１００        情報処理装置
【図１】
2018005671.tif 000007
【図２】
2018005671.tif 000008
【図３】
2018005671.tif 000009
【図４】
2018005671.tif 000010
【図５】
2018005671.tif 000011
【図６】
2018005671.tif 000012
【図７】
2018005671.tif 000013
【図８】
2018005671.tif 000014
【図９】
2018005671.tif 000015
【図１０】
2018005671.tif 000016
【図１１】
2018005671.tif 000017
【図１２】
2018005671.tif 000018
【図１３】
2018005671.tif 000019
【図１４】
2018005671.tif 000020