(57)【要約】  粒子レンダリングの処理負荷の増大を抑制する。粒子レンダリング処理装置は、格子の節点に数値データが定義された格子データを有して構成される対象データ２１を記憶する記憶部２０と、対象データ２１に対する粒子レンダリング処理１００を行う処理部１０と、を備える。粒子レンダリング処理１００は、節点ｉを共有する１又は複数の格子Ｅの体積総和Ｖｉに基づいて、節点ｉに一つの粒子Ｐｉが配置される場合における節点ｉの粒子密度ρｉを推定する処理１０１と、節点ｉに定義された数値データＳｉから、節点ｉにおける不透明度αｉを決定する処理１０３と、節点ｉの粒子密度ρｉと不透明度αｉとに基づいて決定された大きさを有する一つの粒子Ｐｉを節点ｉにおいて生成する処理１０４と、生成された粒子が画像面上に投影された画像データ４０１を生成する処理１０５と、を含む。
2016125754.tif 000002
【特許請求の範囲】
【請求項１】
粒子レンダリング処理装置であって、
格子の節点に数値データが定義された格子データを有して構成される対象データを記憶する記憶部と、
前記対象データに対する粒子レンダリング処理を行う処理部と、
を備え、
前記粒子レンダリング処理は、
前記節点を共有する１又は複数の格子の体積総和に基づいて、前記節点に一つの粒子が配置される場合における前記節点の粒子密度を推定する第１処理と、
前記節点に定義された数値データから、前記節点における不透明度を決定する第２処理と、
前記節点の前記粒子密度と前記不透明度とに基づいて決定された大きさを有する一つの粒子を前記節点において生成する第３処理と、
生成された粒子が画像面上に投影された画像データを生成する第４処理と、
を含む粒子レンダリング処理装置。
【請求項２】
前記処理部は、前記第１処理、前記第２処理、前記第３処理、及び前記第４処理を実行した後に、前記第２処理、前記第３処理、及び前記第４処理を繰り返し実行できるよう構成されている
請求項１記載の粒子レンダリング処理装置。
【請求項３】
前記不透明度は、前記第２処理の繰り返しの度に更新される
請求項２記載の粒子レンダリング処理装置。
【請求項４】
前記粒子レンダリング処理は、
前記対象データに対して、第１粒子レンダリング処理及び第２粒子レンダリング処理のうちのいずれの処理を実行するかを決定する処理選択を更に含み、
前記第１レンダリング処理は、前記第１処理、前記第２処理、前記第３処理、前記第４処理を含み、
前記第２レンダリング処理は、
前記節点に定義された数値データから、前記節点における不透明度を決定する処理と、
前記不透明度とユーザ指定の粒子半径とに基づいて、前記格子の粒子密度を推定する処理と、
前記格子の前記粒子密度に従った数の粒子を生成し、生成された粒子を前記格子内に配置する処理と、
生成された粒子が画像面上に投影された画像データを生成する処理と、
を含む
請求項１〜３のいずれか１項に記載の粒子レンダリング処理装置。
【請求項５】
前記処理選択は、ユーザの選択入力に基づいて決定される
請求項４記載の粒子レンダリング処理装置。
【請求項６】
前記第１処理において、前記粒子密度は、前記節点を共有する１又は複数の格子の体積総和の逆数である
請求項１〜５のいずれか１項に記載の粒子レンダリング処理装置。
【請求項７】
前記第３処理において、前記粒子の大きさは、以下の式に基づいて算出される粒子半径として決定される
請求項１〜６のいずれか１項に記載の粒子レンダリング処理装置。
【数１】
2016125754.tif 000005
ここで、
ｉは、前記節点の識別子
ｒｉは、前記節点ｉの粒子半径
αｉは、前記節点ｉの不透明度
ρｉは、前記節点ｉの粒子密度
Δｔは、レイキャスティングで使用されるレイセグメント長さ
【請求項８】
粒子レンダリング方法であって、
格子の節点に数値データが定義された格子データを有して構成される対象データに対する粒子レンダリング処理を行うことを含み、
前記粒子レンダリング処理は、
前記節点を共有する１又は複数の格子の体積総和に基づいて、前記節点に一つの粒子が配置される場合における前記節点の粒子密度を推定する第１処理と、
前記節点に定義された数値データから、前記節点における不透明度を決定する第２処理と、
前記節点の前記粒子密度と前記不透明度とに基づいて決定された大きさを有する一つの粒子を前記節点において生成する第３処理と、
生成された粒子が画像面上に投影された画像データを生成する第４処理と、
を含む粒子レンダリング方法。
【請求項９】
前記第１処理、前記第２処理、前記第３処理、及び前記第４処理を実行した後に、前記第２処理、前記第３処理、及び前記第４処理を繰り返し実行する
請求項８記載の粒子レンダリング方法。
【請求項１０】
前記不透明度は、前記第２処理の繰り返しの度に更新される
請求項９記載の粒子レンダリング方法。
【請求項１１】
前記対象データに対して、第１粒子レンダリング処理及び第２粒子レンダリング処理のうちのいずれの処理を実行するかを決定することを更に含み、
前記第１粒子レンダリング処理は、前記第１処理、前記第２処理、前記第３処理、前記第４処理を含み、
前記第２レンダリング処理は、
前記節点に定義された数値データから、前記節点における不透明度を決定する処理と、
前記不透明度とユーザ指定の粒子半径とに基づいて、前記格子の粒子密度を推定する処理と、
前記格子の前記粒子密度に従った数の粒子を生成し、生成された粒子を前記格子内に配置する処理と、
生成された粒子が画像面上に投影された画像データを生成する処理と、
を含む
請求項８〜１０のいずれか１項に記載の粒子レンダリング方法。
【請求項１２】
前記第１粒子レンダリング処理及び前記第２粒子レンダリング処理のうちのいずれの処理を実行するかを決定することは、ユーザの選択入力に基づく
請求項１１記載の粒子レンダリング方法。
【請求項１３】
前記第１処理において、前記粒子密度は、前記節点を共有する１又は複数の格子の体積総和の逆数である
請求項８〜１２のいずれか１項に記載の粒子レンダリング方法。
【請求項１４】
前記第３処理において、前記粒子の大きさは、以下の式に基づいて算出される粒子半径として決定される
請求項８〜１３のいずれか１項に記載の粒子レンダリング処理装置。
【数２】
2016125754.tif 000006
ここで、
ｉは、前記節点の識別子
ρｉは、前記節点ｉの粒子密度
ｒｉは、前記節点ｉの粒子半径
αｉは、前記節点ｉの不透明度
Δｔは、レイキャスティングで使用されるレイセグメント長さ
【請求項１５】
格子の節点に数値データが定義された複数の格子データを有して構成される対象データに対する粒子レンダリング処理を、コンピュータに行わせるコンピュータプログラムであって、
前記粒子レンダリング処理は、
前記節点を共有する１又は複数の格子の体積総和に基づいて、前記節点に一つの粒子が配置される場合における前記節点の粒子密度を推定する第１処理と、
前記節点に定義された数値データから、前記節点における不透明度を決定する第２処理と、
前記節点の前記粒子密度と前記不透明度とに基づいて決定された大きさを有する一つの粒子を前記節点において生成する第３処理と、
生成された粒子が画像面上に投影された画像データを生成する第４処理と、
を含むコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、粒子レンダリングに関するものである。
【背景技術】
【０００２】
大規模なボリュームデータの可視化技術として粒子レンダリングが知られている（非特許文献１参照）。粒子レンダリングは、与えられたボリュームデータを、不透明な発光粒子で構成しレンダリングを行う技術である。【０００３】
粒子レンダリングの処理は、粒子生成及び粒子投影という２つの処理を含む。粒子生成の前に、ボリュームデータを構成する格子ごとに粒子密度の推定が必要である。粒子密度とは、単位体積あたりの粒子の数である。格子ごとの粒子密度は、格子に定義された数値データから、不透明度及び粒子半径などを用いて求められる。非特許文献１において、不透明度は、格子に定義された数値データから、ユーザ指定の伝達関数を用いて算出される。【０００４】
粒子生成の処理では、粒子密度に従った数の粒子が生成され、生成された粒子が格子内に配置される。粒子投影の処理では、生成された粒子が画像面に投影される。粒子は画像面上に投影され、画像の各ピクセルで画素値が計算される。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】小山田耕二・坂本尚久共著、”粒子ボリュームレンダリング  理論とプログラミング”コロナ社、２０１０年
【発明の概要】
【発明が解決しようとする課題】
【０００６】
従来の粒子生成は、粒子密度に従った数の粒子を格子内へ配置することを必要とするため、粒子生成手法によっては、処理負荷が非常に大きくなることがある。可視化では、数値データと不透明度の関係を対話的に変更しながらデータの分布や特徴を理解するので、変更のたびに粒子生成処理が必要となる。したがって、非常に大規模なデータを対象とする場合には、処理に要する時間が長くなることがある。【０００７】
このため、粒子レンダリングの処理負荷の増大を抑制することが望まれる。
【課題を解決するための手段】
【０００８】
本発明のある態様は、粒子レンダリング処理装置であって、格子の節点に数値データが定義された格子データを有して構成される対象データを記憶する記憶部と、前記対象データに対する粒子レンダリング処理を行う処理部と、を備え、前記粒子レンダリング処理は、前記節点を共有する１又は複数の格子の体積総和に基づいて、前記節点に一つの粒子が配置される場合における前記節点の粒子密度を推定する第１処理と、前記節点に定義された数値データから、前記節点における不透明度を決定する第２処理と、前記節点の前記粒子密度と前記不透明度とに基づいて決定された大きさを有する一つの粒子を前記節点において生成する第３処理と、生成された粒子が画像面上に投影された画像データを生成する第４処理と、を含む。本発明の他の態様は、粒子レンダリング処理を行うことを含む粒子レンダリング方法である。本発明の他の態様は、コンピュータに、粒子レンダリング処理の一部又は全部を実行させるコンピュータプログラムである。コンピュータプログラムは、コンピュータ読み取り可能な記録媒体に記録される。
【発明の効果】
【０００９】
本発明によれば、粒子レンダリングの処理負荷の増大を抑制することができる。
【図面の簡単な説明】
【００１０】
【図１】粒子レンダリング処理装置のブロック図である。
【図２】粒子レンダリング処理の説明図である。
【図３】粒子レンダリングの処理時間を示すグラフである。
【発明を実施するための形態】
【００１１】
以下、本発明の好ましい実施形態について添付図面を参照しながら説明する。
［１．粒子レンダリング処理装置］
図１に示すように、実施形態に係る粒子レンダリング処理装置１は、処理部１０と、記憶部２０と、を備えている。粒子レンダリング処理装置１は、処理部１０及び記憶部２０を有するコンピュータによって構成されている。処理部１０は、記憶部２０に記憶されたコンピュータプログラムを実行することにより、コンピュータに粒子レンダリング処理を行わせ、コンピュータを粒子レンダリング処理装置１として機能させる。処理部１０は、ＣＰＵ（Central Processing Unit）を有して構成されている。処理部１０は、ＣＰＵのほか、画像処理用のＧＰＵ（Graphics Processing Unit）を有していても良い。処理の高速化のため、粒子レンダリング処理装置１が行う処理のうち、粒子投影処理１０５，２０５は、ＣＰＵではなく、ＧＰＵが行うのが好ましい。
【００１２】
処理部１０が行う粒子レンダリング処理は、記憶部２０に記憶された対象データ２１から、可視化した画像データ４０１，４０２を生成する処理である。生成された画像データ４０１，４０２は、記憶部２０に格納される。粒子レンダリング装置１は、画像を表示するためのモニタを有し、生成された画像データ４０１，４０２をモニタに表示する。
【００１３】
粒子レンダリング処理の対象となる対象データ２１は、例えば、ボリュームデータである。ボリュームデータは、複数の格子データを有して構成される。格子データは、格子の節点（頂点）に数値データが定義されてなるデータである。格子は、規則格子、構造格子、不規則格子のいずれであってもよい。格子の形状は、例えば、四面体、六面体であり、その他の３次元形状であってもよい。数値データは、例えば、スカラデータである。スカラデータは、例えば、温度、圧力、密度、エネルギーを示す値である。
【００１４】
ボリュームデータ２１の例を図２中に示した。ボリュームデータ２１が定義される３次元空間は、格子Ｅと呼ばれる小さな３次元立体に分割される。格子Ｅは複数の節点を有する。図２の格子Ｅは、四面体である。ここでは、ｉを格子Ｅの節点（頂点）の識別子とする。図２の格子Ｅは、節点識別子ｉ＝０〜３で示される４つの節点を有する。各節点ｉ＝０〜３は、スカラデータＳｉなどの数値データが定義されている。つまり、ボリュームデータ２１は、複数の節点ｉ＝０〜３それぞれの位置（座標）と、各節点ｉ＝０〜３に対応付けられたスカラデータＳｉと、を有して構成されている。図２のボリュームデータ２１から明らかなように、ある格子Ｅの節点は、他の格子の節点でもある。つまり、節点は、複数の格子によって共有されることがある。
【００１５】
対象データ２１は、３次元データであるボリュームデータに限られず、２次元データ（面データ）又は１次元データ（線データ）であってもよい。対象データ２１が２次元データである場合、２次元データを構成する格子は、面状の２次元格子となり、厳密には格子厚みが０となる。この場合、２次元格子を、生成される粒子の直径と等しい厚みを有する板状３次元格子とみなして、粒子生成が行われる。同様に、対象データ２１が１次元データある場合、線状の１次元格子を、生成される粒子の直径と等しい直径を有する棒状の３次元格子とみなして、粒子生成が行われる。
【００１６】
処理部１０が行う粒子レンダリング処理は、第１粒子レンダリング処理１００と第２粒子レンダリング処理２００とを含む。後述のように、第１粒子レンダリング処理１００は、節点Ｖの位置だけに粒子Ｐを生成する。このため、第１粒子レンダリング処理１００では、粒子生成を繰り返した場合の処理負荷の増大が抑えられている。また、生成された粒子Ｐのデータを保持するためのメモリ領域の増加も抑えられている。
【００１７】
一方、第２粒子レンダリング処理２００は、従来の粒子レンダリング処理と同様に、粒子密度に従った数の粒子を生成し、生成した粒子Ｐを格子Ｅ内に配置する。このため、第２粒子レンダリング処理２００では、粒子生成を繰り返した場合に処理負荷が増大しやすい。また、生成された粒子のデータを保持するために、多くのメモリ領域を必要とする。ただし、第２粒子レンダリング処理２００では、格子Ｅ内に粒子密度に応じた粒子Ｐの分布が得られるため、生成される画像データ４０２の画質は、第１粒子レンダリング処理１００よりも良好になることが期待される。
【００１８】
処理部１０が行う粒子レンダリング処理では、第１粒子レンダリング処理１００及び第２粒子レンダリング処理１００のうちの、いずれの処理を実行するかを決定するための処理選択３００が行われる。処理選択３００では、対象データ２１に対して、第１粒子レンダリング処理１００を実行するのか、第２粒子レンダリング処理１００を実行するのかが、ユーザの選択入力に基づいて決定される。
【００１９】
処理負荷の小さい第１粒子レンダリング処理１００は、プレビュー用の画像データ４０１を生成するためのプレビュー処理として用いることができる。ユーザは、対象データ２１を可視化した画像データ４０１をプレビュー画像として得たい場合には、第１粒子レンダリング処理（プレビュー処理）を選択することで、比較的短時間で、効率よくプレビュー画像４０１を得ることができる。プレビュー画像４０１の確認後に、第２粒子レンダリング処理２００で画質の良い画像データ４０１を得ることで、ユーザは可視化作業を効率よく行うことができる。
なお、粒子レンダリング処理装置１は、第２粒子レンダリング処理２００を実行する機能を有していなくてもよい。
【００２０】
［２．第１粒子レンダリング処理と第２粒子レンダリング処理］
以下、第１粒子レンダリング処理１００と第２粒子レンダリング処理２００について説明する。以下では、一つの格子に関して説明するが、対象データ２１を構成する他の格子についても同様の処理が行われる。
【００２１】
［２．１  第１粒子レンダリング処理］
第１粒子レンダリング処理１００は、従来の粒子レンダリング処理のように粒子密度に従った数の粒子を格子内に配置するのではなく、節点ｉの位置に粒子Ｐｉを生成する。
【００２２】
図１に示すように、第１粒子レンダリング処理１００では、第１処理として粒子密度推定処理１０１が行われる。粒子密度推定処理１０１は、節点ｉに一つの粒子Ｐｉが配置される場合における前記節点ｉの粒子密度ρｉを推定する処理である。第１粒子レンダリング処理１００で推定される粒子密度ρｉは、各節点ｉ＝１〜３を共有する１又は複数の格子Ｅの体積総和Ｖｉに基づいて推定される。第１粒子レンダリング処理１００では、粒子Ｐｉは各節点ｉに一つずつ配置されるため、粒子密度ρｉは節点ｉの周囲の体積、すなわち当該節点ｉを共有する１又は複数の格子Ｅの体積総和Ｖｉ、に基づいて計算される。より具体的には、当該節点ｉを共有する１又は複数の格子Ｅの体積総和Ｖｉの逆数が、当該節点ｉに一つの粒子Ｐｉが配置される場合における当該節点ｉの粒子密度ρｉとなる。例えば、ある節点ｉを共有する格子Ｅが４つ存在する場合、その４つの格子Ｅの体積総和Ｖｉの逆数である１／Ｖｉが、当該節点ｉの粒子密度ρｉとなる。
【００２３】
第１粒子レンダリング処理１００では、さらに、第２処理として不透明度更新処理（不透明度決定処理）１０３が行われる。不透明度更新処理１０３は、各節点ｉ＝０〜３の不透明度αｉを決定（更新）する処理である。不透明度αｉは、各節点ｉに定義された数値データ（スカラデータＳｉ）から、ユーザ指定の伝達関数を用いて算出される。不透明度の算出は、従来の粒子レンダリング処理と同様である。
【００２４】
第１粒子レンダリング処理１００では、さらに、第３処理として粒子半径計算処理１０４が行われる。粒子半径計算処理１０４は、節点ｉに配置される粒子の大きさを決定し、決定された大きさの粒子を生成する処理である。粒子の大きさは、例えば、粒子Ｐｉの半径ｒｉによって決定される。ある節点ｉの粒子半径ｒｉは、当該節点ｉの粒子密度ρｉと不透明度αｉに基づいて決定される。より具体的には、粒子半径ｒｉは以下の式（１）に基づいて算出される。
【数１】
2016125754.tif 000003
【００２５】
上記式（１）は、特許文献１に記載の従来の粒子密度推定式を変形して得たものである。従来の粒子密度推定式は、以下の式（２）のとおりである。
【数２】
2016125754.tif 000004
【００２６】
式（２）に示す従来の粒子密度推定式では、不透明度αと粒子半径ｒとが決まると、粒子密度ρを求めることができる。なお、Δｔはレイキャスティングで使用されるレイセグメント長さである（非特許文献１参照）。また、従来の式（２）では、粒子半径ｒは、ユーザ指定の値であり、不透明度αは、節点に定義された数値データから、ユーザ指定の伝達関数を用いて算出される値である。
【００２７】
従来の式（２）では、不透明度αと粒子半径ｒとから粒子密度ρを求めていたのに対して、式（１）では、節点ｉの粒子半径ｒｉを、粒子密度推定処理１０１で求めた粒子密度ρｉと、不透明度更新処理１０３で求めた不透明度αｉと、から求める。
図２に示すように、求めた粒子半径ｒｉを有する粒子Ｐｉが、節点ｉの位置に生成される。
【００２８】
第１粒子レンダリング処理１００では、粒子半径計算処理１０４に続いて、第４処理として粒子投影処理１０５が行われる。粒子投影処理１０５は、生成された粒子が画像面に投影された画像データを生成する処理である。粒子投影処理１０５の詳細は、非特許文献１に記載されおり、本実施形態の粒子投影処理１０５も、非特許文献１の記載に従って行われる。
なお、処理部１０は、不透明度更新処理１０３及び粒子半径計算処理１０４を含む処理１０２と、粒子投影処理１０５と、を複数回繰り返し実行することができる。処理１０２を繰り返し行う際には、繰り返しの度に、不透明度αｉを変化（更新）させて、生成される画像データを変化させることができる。
【００２９】
従来の粒子レンダリング処理では、粒子密度ρは、式（２）に示すように、不透明度αに応じて決まるため、不透明度αを更新すると、粒子密度ρを推定し直して、粒子密度ρに応じた数の粒子を生成する処理を再度行う必要が生じる。したがって、不透明度を更新して処理を繰り返すと、その分、処理負荷が増大する。
【００３０】
これに対し、第１粒子レンダリング処理１００では、粒子Ｐｉは、節点ｉだけに配置されるため、単位体積あたりの粒子数を示す粒子密度ρｉは、不透明度αｉとは無関係に決まり、不透明度αｉが更新されても、変化しない。また、第１粒子レンダリング処理１００では、粒子Ｐｉは、節点ｉだけに配置されるため、粒子密度に応じた数の粒子を生成する処理も不要である。したがって、処理を繰り返しても、処理負荷の増大が抑制される。
【００３１】
なお、第１レンダリング処理では、処理１０２の複数回の繰り返しで得られた複数の画素値を加算していき、最後に繰り返し数で割り算して得られたアンサンブル平均を最終画素値とすることができる（非特許文献１参照）。
【００３２】
［２．２  第２粒子レンダリング処理］
第２粒子レンダリング処理は、従来の粒子レンダリング処理と同様に、粒子密度に従った数の粒子を格子内に配置する。
【００３３】
図１に示すように、第２粒子レンダリング処理２００では、不透明度更新処理（不透明度決定処理）２０２と、粒子密度推定処理２０３と、粒子生成処理２０４と、粒子投影処理２０５と、が行われる。これらの処理は、従来の粒子レンダリング処理と同様である。
つまり、不透明度更新処理２０２は、各節点ｉ＝０〜３の不透明度αを決定（更新）する処理である。不透明度αｉは、各節点ｉに定義された数値データ（スカラデータＳｉ）から、ユーザ指定の伝達関数を用いて算出される。
【００３４】
粒子密度推定処理２０３では、格子Ｅごとの粒子密度ρが推定される（非特許文献１参照）。粒子密度ρの推定は、前述の式（２）を用いて、不透明度αとユーザ指定の粒子半径ｒとに基づいて行われる。
粒子生成処理２０４では、粒子密度ρに従った数の粒子Ｐが生成され、生成された粒子Ｐが格子Ｅ内に配置される。粒子生成処理２０４において採用される粒子生成手法は、例えば、一様サンプリング法、メトロポリスサンプリング法である（非特許文献１参照）。一様サンプリング法は、格子Ｅ内の粒子密度を一定と仮定して粒子を生成する手法である。メトロポリスサンプリング法は、格子Ｅ内部における粒子密度変化を考慮に入れた粒子生成手法である。図２では、一様サンプリング法及びメトロポリスサンプリング法それぞれで生成される粒子Ｐの例を示した。
【００３５】
第２粒子レンダリング処理２００では、図２に示すように、粒子Ｐが粒子密度ρに従って格子Ｅ内に多数配置されるため、画像データ４０２の画質が良好になることが期待される。
【００３６】
第２粒子レンダリング処理２００では、粒子生成処理２０４に続いて、粒子投影処理２０５が行われる。粒子投影処理２０５は、生成された粒子が画像面に投影された画像データを生成する処理である（非特許文献１参照）。
なお、処理部１０は、不透明度更新処理２０２、粒子密度推定処理２０３及び粒子生成処理２０４からなる処理２０１と、粒子投影処理２０５と、を複数回繰り返し実行することができる。また、処理部１０は、複数回の繰り返しで得られた複数の画素値を加算していき、最後に繰り返し数で割り算して得られたアンサンブル平均を最終画素値とすることができる（非特許文献１参照）。
【００３７】
［３．処理時間の比較］
図３は、第１粒子レンダリング処理１００と第２粒子レンダリング処理２００それぞれの処理時間を示している。図３において、「初期描画」は、第１粒子レンダリング処理１００については、粒子密度推定処理１０１から粒子投影処理１０５までの一連の処理を１回行った場合の処理時間を示し、第２粒子レンダリング処理２００については、不透明度更新処理２０２から粒子投影処理２０５までの処理を１回行った場合の処理時間を示している。また、図３において、「第１回繰り返し」は、第１粒子レンダリング処理１００については、「初期描画」の後に、不透明度更新１０３から粒子投影処理１０５までの一連の処理を再度行った場合の処理時間を示し、第２粒子レンダリング処理１００については、「初期描画」の後に不透明度更新処理２０２から粒子投影処理２０５までの処理を再度行った場合の処理時間を示している。さらに、図３において、「第２回繰り返し」は、第１粒子レンダリング処理１００については、「第１回繰り返し」の後に、不透明度更新１０３から粒子投影処理１０５までの一連の処理を再度行った場合の処理時間を示し、第２粒子レンダリング処理１００については、「第１回繰り返し」の後に不透明度更新処理２０２から粒子投影処理２０５までの処理を再度行った場合の処理時間を示している。
【００３８】
図３から明らかなように、第２粒子レンダリング処理２００では、繰り返しの度に比較的長い処理時間を要するのに対し、第１粒子レンダリング処理１００では、初期描画後の処理の繰り返しでは、処理時間を短縮できており、処理負荷の増大を抑制することができる。
【符号の説明】
【００３９】
１  粒子レンダリング処理装置
１０  処理部
２０  記憶部
２１  対象データ
１０１  粒子密度推定処理
１０３  不透明度更新処理
１０４  粒子半径計算処理
１０５  粒子投影処理
４０１  画像データ
【図１】
2016125754.tif 000007
【図２】
2016125754.tif 000008
【図３】
2016125754.tif 000009
【国際調査報告】
2016125754.tif 000010
2016125754.tif 000011