拡散場における複数のパーティクルを可視化する可視化処理プログラムにおいて、コンピュータに、前記複数のパーティクルのうち、パーティクルの位置及び速度のデータに基づき、パーティクル間の距離が所定距離以下且つ前記パーティクル間の速度ベクトルの差の絶対値が所定値以下であるパーティクルが同じグループになるように、前記複数のパーティクルのグループ分けを行わせ、前記グループ分けされたグループの各々について、グループに属するパーティクルを含むグループ表示データを生成させ、前記グループに属するパーティクルを表示するパーティクル表示データに前記グループ表示データを付加した表示データを生成させることを特徴とする可視化処理プログラム。前記可視化処理プログラムはさらに、前記コンピュータに、前記グループ分けされたグループの各々について、グループに属するパーティクルの位置のデータを用いて、第1の軸方向におけるグループの位置を特定させ、前記特定された第1の軸方向におけるグループの位置に応じた色に、グループ表示の色を設定させることを特徴とする請求項1記載の可視化処理プログラム。前記可視化処理プログラムはさらに、前記コンピュータに、前記グループ分けされたグループの各々について、グループに属するパーティクルのうち、グループに属するパーティクルの位置の平均に最も近いパーティクルを特定させ、前記特定されたパーティクルの位置から所定時間後の前記特定されたパーティクルの位置までの流線のデータを生成させ、前記生成された流線のデータを前記グループ表示データに付加させることを特徴とする請求項1又は2記載の可視化処理プログラム。前記可視化処理プログラムはさらに、前記コンピュータに、前記グループ分けされたグループの各々について、グループに属するパーティクルの間を結ぶ線分を、ループを形成しないように、短い順にN-1Nは当該グループに属するパーティクルの数本接続させることにより、第1の構造のデータを生成させ、前記第1の構造に含まれる線分の太さを太くさせ、前記第1の構造に含まれる凹み部分を平滑化させることにより、前記グループ表示データを生成させることを特徴とする請求項1乃至3のいずれか1つ記載の可視化処理プログラム。前記可視化処理プログラムはさらに、前記コンピュータに、前記グループ分けされたグループの各々について、グループに属するパーティクルの位置のデータを用いて、最小二乗法に基づいた直線を特定させ、前記特定された直線とグループに属するパーティクルとの位置関係に基づき、グループに属するパーティクルを第1のグループと第2のグループとにグループ分けさせ、前記第1のグループに属する各パーティクルから前記直線に引いた垂線が前記直線と交わる点の位置を特定させ、各々から引かれた垂線が前記直線と交わる点の位置が隣り合っているパーティクルの間を線分で接続させることにより、第2の構造のデータを生成させ、前記第2のグループに属する各パーティクルから前記直線に引いた垂線が前記直線と交わる点の位置を特定させ、各々から引かれた垂線が前記直線と交わる点の位置が隣り合っているパーティクルの間を線分で接続させることにより、第3の構造のデータを生成させ、前記第2の構造における端と第3の構造における端との間を線分で接続させることにより、第4の構造のデータを生成させ、前記第4の構造の凹み部分を平滑化させることにより、前記グループ表示のデータを生成させることを特徴とする請求項1乃至3のいずれか1つ記載の可視化処理プログラム。拡散場における複数のパーティクルを可視化する可視化処理方法において、コンピュータが、前記複数のパーティクルのうち、各パーティクルの位置及び速度のデータに基づき、パーティクル間の距離が所定距離以下且つ前記パーティクル間の速度ベクトルの差の絶対値が所定値以下であるパーティクルが同じグループになるように、前記複数のパーティクルをグループ分けし、前記グループ分けされたグループの各々について、グループに属するパーティクルを含むグループ表示データを生成し、前記グループに属するパーティクルを表示するパーティクル表示データに前記グループ表示データを付加した表示データを生成することを特徴とする可視化処理方法。拡散場における複数のパーティクルの位置及び速度のデータに基づき、パーティクル間の距離が所定距離以下且つ前記パーティクル間の速度ベクトルの差の絶対値が所定値以下であるパーティクルが同じグループになるように、前記複数のパーティクルをグループ分けするグルーピング部と、前記グループ分けされたグループの各々について、グループに属するパーティクルを含むグループ表示データを生成するグループ表示生成部と、前記グループに属するパーティクルを表示するパーティクル表示データに前記グループ表示データを付加した表示データを生成する表示データ生成部と、を有することを特徴とする可視化処理装置。本技術の目的は、一側面においては、拡散場におけるパーティクルの動きを容易に把握できるようにするための技術を提供することである。本技術の一側面に係る可視化処理方法は、コンピュータが、A複数のパーティクルのうち、各パーティクルの位置及び速度のデータに基づき、パーティクル間の距離が所定距離以下且つパーティクル間の速度ベクトルの差の絶対値が所定値以下であるパーティクルが同じグループになるように、複数のパーティクルをグループ分けし、Bグループ分けされたグループの各々について、グループに属するパーティクルを含むグループ表示データを生成し、Cグループに属するパーティクルを表示するパーティクル表示データにグループ表示データを付加した表示データを生成することを特徴とする。拡散場におけるパーティクルの動きを容易に把握できるようになる。[実施の形態1]図3に、第1の実施の形態に係る情報処理装置1の機能ブロック図を示す。情報処理装置1は、入力部101と、描画オブジェクト格納部102と、パーティクルフロー生成部103と、第1描画データ格納部104と、シミュレーション結果格納部105と、速度算出部106と、速度データ格納部107と、群表示生成部108と、第2描画データ格納部115と、描画処理部116と、第3描画データ格納部117と、グラフィックス処理部118と、表示部119とを含む。また、群表示生成部108は、群化処理部109と、オブジェクト生成部110と、流線生成部111と、色付け処理部112と、群化テーブル格納部113と、パーティクルテーブル格納部114とを含む。入力部101は、情報処理装置1のユーザから可視化の指示を受け付けると、パーティクルフロー生成部103に可視化処理の要求を行う。速度算出部106は、シミュレーション結果格納部105に格納されているシミュレーションデータを用いて速度場データ例えば、3次元におけるパーティクルの速度のデータを計算し、計算した速度場データを速度データ格納部107に格納する。パーティクルフロー生成部103は、シミュレーション結果格納部105に格納されているシミュレーションデータ、描画オブジェクト格納部102に格納されている描画オブジェクトデータ及び速度データ格納部107に格納されている速度場データを用いて、パーティクルフローによって流れを表示するための描画データを生成し、第1描画データ格納部104に格納する。群化処理部109は、速度データ格納部107に格納されている速度場データを用いてパーティクルの群化処理を行い、群化処理の結果である群化データを群化テーブル格納部113に格納する。オブジェクト生成部110は、群化テーブル格納部113に格納されている群化データ及びパーティクルテーブル格納部114に格納されているパーティクルテーブルのデータを用いて群表示のデータを生成する処理を行い、生成した群表示のデータを第2描画データ格納部115に格納する。流線生成部111は、オブジェクト生成部110によって生成された群表示の流線を生成する処理を行う。色付け処理部112は、オブジェクト生成部110によって生成された群表示に色付けする処理を行う。描画処理部116は、第1描画データ格納部104に格納されている描画データ及び第2描画データ格納部115に格納されている群表示の描画データを用いて処理を行い、処理結果である表示データを第3描画データ格納部117に格納する。グラフィックス処理部118は、第3描画データ格納部117に格納されているデータを表示するための周知の処理を行う。表示部119は、グラフィックス処理部118の処理結果を表示する。次に、図4乃至図34を用いて、図3に示した情報処理装置1の動作について説明する。なお、シミュレーション結果格納部105には、拡散方程式等を用いて計算された、拡散場におけるパーティクルの位置に関するシミュレーション結果が格納されているとする。本実施の形態では、時刻t=0から時刻t=NNは予め定められた自然数までについてのシミュレーション結果が格納されているとする。また、計算はシミュレータ図示しないによって行われ、その結果がシミュレーション結果格納部105に格納されているとする。なお、以下ではt及びgは0以上の整数、m及びnは自然数であるとする。まず、速度算出部106は、例えばユーザからの指示に応じて、シミュレーション結果格納部105に格納されているデータを用いて速度場を計算し、計算結果である速度場データを速度データ格納部107に格納する図4:ステップS1。図5に、速度データ格納部107に格納されている速度場データの一例を示す。図5の例では、パーティクル番号と、パーティクルのx軸方向の位置と、パーティクルのy軸方向の位置と、パーティクルのz軸方向の位置と、パーティクルのx軸方向の速度と、パーティクルのy軸方向の速度と、パーティクルのz軸方向の速度とが格納されている。このような、各パーティクルについての位置及び速度を含む速度場のデータが、各時刻t=0乃至Nについて格納されている。そして、群表示生成部108における群化処理部109は、時刻tの初期値を0に設定しステップS3、また群番号gの初期値を0に設定するステップS5。群化処理部109は、パーティクル番号nを1に設定するステップS7。そして、m≧2且つm>nを満たす最小の自然数を特定し、特定された自然数をパーティクル番号mに設定するステップS9。そして、群化処理部109は、速度データ格納部107に格納されている位置のデータを用いて、時刻tにおける2つのパーティクルここでは、パーティクルmとパーティクルnの位置の差すなわち距離Δpを計算し、メインメモリ等の記憶装置図示しないに格納するステップS11。そして、群化処理部109は、ΔpがΔp<pgを満たすか判断するステップS13。ここでpgは、予め設定された、距離についての基準値である。Δp<pgを満たさないと判断された場合ステップS13:Noルート、パーティクルmはパーティクルnと同じ群には属さないので、処理は端子Aを介して図6のステップS21に移行する。一方、Δp<pgを満たすと判断された場合ステップS13:Yesルート、群化処理部109は、速度データ格納部107に格納されている速度のデータを用いて、時刻tにおけるパーティクルmの速度ベクトルとパーティクルnの速度ベクトルとの差の速度ベクトルの絶対値Δvを計算し、メインメモリ等の記憶装置図示しないに格納するステップS15。そして、群化処理部109は、ΔvがΔv<vgを満たすか判断するステップS17。ここでvgは、予め設定された、速度の大きさについての基準値である。なお、ステップS15においては2つの速度ベクトルの差分の速度ベクトルを求め、求められた速度ベクトルの絶対値を計算している。そのため、2つの速度ベクトルの向きが大きく異なる場合には、2つの速度ベクトルの向きが近い場合と比べると、差分の速度ベクトルの絶対値が大きくなる。そのため、2つの速度ベクトルの向きが大きく異なる場合には、2つのパーティクルは異なる群に属すると判断される。Δv<vgを満たさないと判断された場合ステップS17:Noルート、パーティクルmはパーティクルnと同じ群には属さないので、処理は端子Aを介して図6のステップS21に移行する。一方、Δv<vgを満たすと判断された場合ステップS17:Yesルート、パーティクルmはパーティクルnと同じ群に属するので、群化処理部109は、パーティクルテーブル格納部114に格納されているパーティクルテーブルにおいて、パーティクルm及びnの群番号をgに設定するステップS19。処理は端子Aを介して図6のステップS21に移行する。図7に、パーティクルテーブル格納部114に格納されているパーティクルテーブルの一例を示す。図7の例では、パーティクル番号と、群番号と、流線フラグとが格納されている。流線フラグについては、後で説明する。図6の説明に移行し、群化処理部109は、未処理のmが有るか判断するステップS21。未処理のパーティクルmが有ると判断された場合ステップS21:Yesルート、次のmについて処理するため、群化処理部109は、m=m1と設定するステップS23。そして処理は端子Bを介してステップS11に戻る。一方、未処理のmが無いと判断された場合ステップS21:Noルート、群化処理部109は、群を表す群番号gをg=g1と設定するステップS25。そして、未処理のnが有るか判断するステップS27。未処理のnが有ると判断された場合ステップS27:Yesルート、次のnについて処理するため、群化処理部109は、n=n1と設定するステップS29。そして処理は端子Cを介してステップS9に戻る。一方、未処理のnが無いと判断された場合ステップS27:Noルート、オブジェクト生成部110、流線生成部111及び色付け処理部112は、時刻tについて群表示生成処理を実施するステップS31。群表示生成処理については、図8乃至図30を用いて説明する。まず、オブジェクト生成部110は、群番号gを0に設定し図8:ステップS41、またパーティクル番号nを1に設定するステップS43。そして、m≧2且つm>nを満たす最小の自然数を特定し、特定された自然数をパーティクル番号mに設定するステップS45。そして、オブジェクト生成部110は、パーティクルテーブルをパーティクル番号m及びnで検索し、パーティクルm及びパーティクルnの群番号がいずれもgであるか判断するステップS47。すなわち、パーティクルm及びパーティクルnが同じ群に属するかを判断する。パーティクルm及びパーティクルnの群番号のうちいずれかがgではないと判断された場合ステップS47:Noルート、ステップS51の処理に移行する。一方、パーティクルm及びパーティクルnの群番号がいずれもgであると判断された場合ステップS47:Yesルート、オブジェクト生成部110は、速度データ格納部107に格納されている位置のデータを用いて、時刻tにおけるパーティクルmとパーティクルnとを結ぶ線分の長さ、つまりパーティクル間の距離を算出し、群化テーブル格納部113に格納するステップS49。図9に、群化テーブル格納部113に格納されている群化テーブルの一例を示す。図9の例では、パーティクルを結ぶ線分の識別子と、線分の両端にあるパーティクルの番号と、線分の長さと、通過フラグとが格納されている。通過フラグについては後で説明する。なお、このような群化テーブルは、各群について生成される。そして、オブジェクト生成部110は、未処理のmが有るか判断するステップS51。未処理のmが有ると判断された場合ステップS51:Yesルート、次のmについて処理するため、群化処理部109は、m=m1と設定するステップS53。そしてステップS47の処理に戻る。一方、未処理のmが無いと判断された場合ステップS51:Noルート、オブジェクト生成部110は、未処理のnが有るか判断するステップS55。未処理のnが有ると判断された場合ステップS55:Yesルート、次のnについて処理するため、オブジェクト生成部110は、n=n1と設定するステップS57。そしてステップS45の処理に戻る。一方、未処理のnが無いと判断された場合ステップS55:Noルート、処理は端子Fを介して図10のステップS59に移行する。図10の説明に移行し、オブジェクト生成部110は、群番号gの群について、群化テーブルのデータを線分の長さでソートするステップS59。例えば図9に示した群化テーブルのデータをソートすると、図11に示すような状態の群化テーブルとなる。そして、オブジェクト生成部110は、スケルトン構造に利用するか判定する処理が未だ行われていない線分のうち、長さが最も短い線分を特定するステップS61。また、オブジェクト生成部110は、特定された線分を、生成途中のスケルトン構造スケルトン構造については、後で詳細に説明するに追加することによって当該構造内にループが形成されるか判断するステップS63。線分を追加することによってループが形成されるか否かの判断については甲jつする。ループが形成されると判断された場合ステップS63:Yesルート、オブジェクト生成部110は、処理対象の線分のデータを群化テーブルから削除するステップS64。そして、次の線分について処理するため、ステップS61の処理に戻る。一方、ループが形成されないと判断された場合ステップS63:Noルート、オブジェクト生成部110は、特定された線分の描画データを生成し、記憶装置図示しないに格納するステップS65。そして、オブジェクト生成部110は、描画データを生成した線分の本数がPnー1に達したか判断するステップS67。ここで、Pnは群番号gの群に属するパーティクルの数とする。描画データを生成した線分の本数がPnー1に達していないと判断された場合ステップS67:Noルート、次の線分について処理するため、ステップS61の処理に戻る。ここで、ステップS61乃至S67において行われる処理について、具体例を用いて説明する。ステップS61乃至S67の処理は、パーティクル間に線分を引いて、図12に示すような構造以下、スケルトン構造と呼ぶのデータを生成する処理である。なお、図12及びその他の図面においては、説明の便宜上、スケルトン構造にパーティクルが含まれているが、本実施の形態におけるスケルトン構造には本来、線分は含まれるがパーティクルは含まれないものとする。図13に示すような、7つのパーティクルを含む群についてスケルトン構造を生成することを考える。ステップS61乃至S67の処理によれば、長さが短い線分から順に線分の描画データが生成されるので、図13に示した群に対して処理を行うと、途中の段階で図14に示すような構造が生成される。そして、群化テーブルに図11に示すようなデータが格納されている場合には、次に処理対象となる線分は、パーティクル1とパーティクル4とを結ぶ線分L3である。ステップS63においては、生成途中のスケルトン構造への線分の追加によってループが形成されるかを判断するが、ここではパーティクル1、2及び4を相互に接続する線分によりループが形成されると判断される。この判断は、具体的には以下の方法によって行われる。まず、追加しようとする線分L3の両端のパーティクルパーティクル1及びパーティクル4のうちいずれかのパーティクルを選択する。ここでは、パーティクル1を選択するものとする。そして、パーティクル1から始まる探索経路を考える。具体的には、群化テーブルのデータを用いて、パーティクル1を根とするn分木構造を生成する。はじめに、既に描画データが生成された、パーティクル1とパーティクル2とを結ぶ線分L1と、今回追加する線分L3とによって図15に示すような木構造を生成する。この段階において、群化テーブルには、例えば図16に示すようなデータが格納される。図16の例では、線分L1及び線分L3に通過フラグ1が対応付けられている。なお、通過フラグは、n分木構造に含まれる線分を特定するためのデータである。次に、n分木構造の末端であるパーティクル2及び4について、既に描画データが生成された線分であるL8及びL7を特定し、n分木構造に追加する。線分L8は、図14においてパーティクル2とパーティクル4とを結ぶ線分であり、線分L7は、図14においてパーティクル2とパーティクル3とを結ぶ線分である。この処理によって、図17に示すような構造が生成される。この段階において、群化テーブルには、例えば図18に示すようなデータが格納される。図18の例では、図17のn分木構造に含まれる線分に、通過フラグ1が対応付けられている。ここでは、新たに線分L7及び線分L8に通過フラグ1が対応付けられている。ここで、図17に示したn分木構造には、パーティクル2が2つ含まれている。これは、スケルトン構造にループが形成されたことを示している。従って、この段階において、スケルトン構造に線分L3を追加することはできないと判断され、線分L3のデータは群化テーブルから削除される。なお、ステップS63の処理が終了すると、群化テーブルにおける通過フラグの列は初期化され、図19に示すような状態になる。図19は、線分L3のデータが削除された状態の群化テーブルを示している。そして、最終的には、図20のようなスケルトン構造が生成される。以上のように、パーティクル間をつなぐ線分を、ループが形成されないように、短いものから順につないでいくことにより、スケルトン構造を生成する。図10の説明に戻り、描画データを生成した線分の本数がPnー1に達したと判断された場合ステップS67:Yesルート、オブジェクト生成部110は、未処理の群が有るか判断するステップS69。未処理の群が有ると判断された場合ステップS69:Yesルート、次の群について処理するため、オブジェクト生成部110は、群番号gをg=g1と設定するステップS71。そして処理は端子Eを介してステップS43に戻る。一方、未処理の群が無いと判断された場合ステップS69:Noルート、処理は端子Gを介して図21のステップS73に移行する。図21の説明に移行し、オブジェクト生成部110は、時刻tについて生成されたスケルトン構造に対する肉付け処理を実施し、処理結果を記憶装置図示しないに格納するステップS73。ステップS73においては、モルフォロジmorphology処理例えば、ダイレーションdilation処理を実施する。このモルフォロジ処理はよく知られた処理であるので、ここでは説明を省略する。例えば図12に示したスケルトン構造に対してステップS73の処理を実施すると、図22に示すような構造が生成される。そして、オブジェクト生成部110は、ステップS73の処理によって生成された構造に対して凹み部分の補正及び角を丸める処理を行うことにより、群表示のデータを生成し、記憶装置図示しないに格納するステップS75。凹み部分の補正とは、肉付け処理された部分における窪み部分を、さらに肉付け処理する等して平滑化する処理である。角を丸める処理とは、例えば線分の端点であった部分等、尖った部分の周りに肉付け処理を行うなどして、尖った部分を無くす処理である。ステップS75の処理を行うと、群表示は楕円体又は球体に近い形状になる。例えば図22に示した構造に対してステップS75の処理を実施すると、図23に示すような群表示が生成される。色付け処理部112は、速度データ格納部107に格納されている位置のデータを用いて、時刻tにおける各群のz軸方向の位置を特定するステップS76。ここでは、群に属するパーティクルのz軸方向の座標値の平均を、群のz軸方向の位置の座標値であるとする。そして、色付け処理部112は、z軸方向における群の位置に応じて、群表示に対して色付け処理を実施し、処理結果を記憶装置図示しないに格納するステップS77。例えば図24に示すように、z軸方向の群の位置とカラースケール上の色とを対応付けることにより、各群の色を決定する。図24の例では、群表示241の色は赤であり、群表示242の色は黄緑色であり、群表示243の色は青である。図25は、図24に示した座標空間を、図24の上部に示した、z座標が正の位置から原点に向かう視線方向から観察したものである。このように色付け処理を行うことによって、流れの三次元的な構造をより容易に把握できるようになる。そして、流線生成部111は、流線描画処理を実施するステップS79。流線描画処理については後で説明する。そして元の処理に戻る。以上のようにして、同一の群に含まれるパーティクルを覆う群表示を生成する。このような表示を用いれば、単純にパーティクルのみを表示するよりも、流れの状況を容易に把握できるようになる。図26に、ここまで説明した処理によって生成される群表示の一例を示す。図26の例では、同一の群に属するパーティクルを覆う群表示が単一色で示されている。なお、図を見やすくするため、生成された群表示のうち一部の群表示のみを示している。また、視線の方向はz軸方向であるとする。また、図26に示した群表示に、さらにz軸方向の位置に応じて色付け処理を行ったものを図27に示す。このような色付け処理を行うことにより、例えばz軸方向に一様な流れになっている場合においても、流れの三次元的な構造を適切に把握することができるようになる。次に、図28を用いて、流線描画処理ステップS79について説明する。まず、流線生成部111は、パーティクルテーブルから未処理の群以下、処理対象の群と呼ぶを1つ特定する図28:ステップS81。また、流線生成部111は、処理対象の群に含まれるパーティクルの平均位置Paを、速度データ格納部107に格納されている位置のデータを用いて算出し、算出結果をメインメモリ等の記憶装置図示しないに格納するステップS83。そして、流線生成部111は、処理対象の群の中から未処理のパーティクルを1つ特定するステップS85。また、ステップS85において特定されたパーティクルと、ステップS83において算出されたPaとの距離を算出し、メインメモリ等の記憶装置図示しないに格納するステップS87。そして、流線生成部111は、未処理のパーティクルが有るか判断するステップS89。未処理のパーティクルが有ると判断された場合ステップS89:Yesルート、次のパーティクルについて処理するため、ステップS85の処理に戻る。一方、未処理のパーティクルが無いと判断された場合ステップS89:Noルート、流線生成部111は、処理対象の群に属するパーティクルのうち、ステップS87において算出された距離が最短であるパーティクルについて、パーティクルテーブルにおける流線フラグの列に1を設定するステップS91。そして、流線生成部111は、未処理の群が有るか判断するステップS93。未処理の群が有ると判断された場合ステップS93:Yesルート、次の群について処理するため、ステップS81の処理に戻る。一方、未処理の群が無いと判断された場合ステップS93:Noルート、流線生成部111は、時刻tにおける各パーティクルについて、速度データ格納部107に格納されているデータを用いて時刻tからΔt秒後の位置を特定するステップS95。そして、流線生成部111は、時刻tにおける各パーティクルについて、時刻tにおける位置からΔt秒後の位置までの流線のデータを生成して群表示のデータに付加し、第2描画データ格納部115に格納するステップS97。なお、パーティクルテーブルにおいて流線フラグが1に設定されている代表パーティクルについては、流線の色を他のパーティクルの流線の色とは異なる色にする。そして元の処理に戻る。図29に、群表示に流線を付加した場合の表示例を示す。図29の例では、群に属するパーティクルの各々に流線が付されており、特に主要な流れの方向に近い方向に動く代表パーティクルには他の流線とは異なる色の流線が付されている。これによって、群の大まかな移動方向を直感的に把握しやすくなる。なお、説明をわかりやすくするためパーティクルを表示しているが、群表示生成処理によって生成されるデータには、パーティクルを表示するためのデータは含まれない。また、上では時刻tにおける各パーティクルについて生成した流線を、時刻tにおける各パーティクルに付加する例を示したが、時刻tΔtにおける各パーティクルに付加するようにしてもよい。この場合の表示例を図30に示す。図30のような表示を、図29に示した表示の後に示すことで、時間の経過に伴う群の位置の推移をより容易に把握できるようになる。以上のような処理を実施することで、群表示及び流線のデータが生成され、第2描画データ格納部115に格納される。一方で、パーティクルフロー生成部103は、シミュレーション結果格納部105、描画オブジェクト格納部102及び速度データ格納部107に格納されているデータを用いて、パーティクルフローによって流れを表示するためのデータを生成する。そして、パーティクルフロー生成部103は、生成したデータを第1描画データ格納部104に格納する。例えば、まず図1及び図2に示すような表示を行うためのデータを生成する。パーティクルフローによって流れを表示するためのデータを生成する方法はよく知られているので、詳細な説明は省略する。そして、z軸方向の位置に応じてパーティクル及び流線の色付けを行うことにより、図31及び図32に示すような表示を行うためのデータを生成する。なお、色付けのルールは、ステップS77の処理の説明で述べたものと同様である。このようにして生成されたデータは、描画処理部116の処理に利用される。図6の説明に戻り、群表示生成処理の後、群化処理部109は、未処理のtが有るか判断するステップS33。未処理のtが有ると判断された場合ステップS33:Yesルート、次の時刻について処理するため、群化処理部109は、t=t1と設定するステップS35。そして処理は端子Dを介してステップS7の処理に戻る。一方、未処理のtが無いと判断された場合ステップS33:Noルート、描画処理部116は、パーティクルフローによって流れを表示するためのデータを第1描画データ104から読み出し、また群表示のデータを第2描画データ格納部115から読み出す。そして、描画処理部116は、パーティクル及び流線と群表示とを重ね合わせて表示するための処理を行い、処理結果である表示データを第3描画データ格納部117に格納するステップS37。ここまでの処理によって、例えば図33及び図34に示すような表示を行うための表示データが生成される。図33の例では、同一の群に属するパーティクルに群表示が付されている。そして、図34の例では、図33の例における群表示にさらに色付け処理が行われている。なお、図33及び図34においては、図を見やすくするため、群表示生成部108によって生成された流線のデータは使用していない。そして、グラフィックス処理部118は、第3描画データ格納部117に格納されている表示データを表示部119に表示させる処理を実施するステップS39。そして処理を終了する。ステップS39においては、例えば以下のような処理を行う。すなわち、図33及び図34に示すような表示を行うための表示データは各時刻tについて生成されているので、これらを連続的に表示する。また、第3描画データ格納部117に格納されている表示データのうち、パーティクルフローによって流れを表示するためのデータについては表示しないようにしてもよい。例えば、図26及び図27に示すような群表示に、ステップS79又はS127において生成された流線のデータを付加し、連続的に表示するようにしてもよい。以上のような処理を実施すれば、時間Δt前に同じ位置にあったと推定されるパーティクルが同じ群に属するようになるので、拡散場におけるパーティクルの動きを容易に把握できるようになる。[実施の形態2]次に、第2の実施の形態について説明する。第2の実施の形態においては、第1の実施の形態とは異なる群表示生成処理によって群表示を生成する。以下、第2の実施の形態における群表示生成処理について説明する。まず、群表示生成部108におけるオブジェクト生成部110は、群番号gの初期値を0に設定する図35:ステップS101。そして、オブジェクト生成部110は、群番号gに属するパーティクルの位置のデータを速度データ格納部107から読み出し、各パーティクルからの距離の自乗の総和が最小となるような直線を最小自乗法によって決定し、決定された直線のデータを記憶装置図示しないに格納するステップS103。ステップS103の処理によって、図36に示すような直線のデータが生成される。そして、オブジェクト生成部110は、群番号gに属するパーティクルのうち未処理のパーティクル以下、処理対象のパーティクルと呼ぶを1つ特定するステップS105。また、オブジェクト生成部110は、処理対象のパーティクルから直線へ垂線を引いた場合における垂線と直線との交点の位置を特定し、群化テーブルに格納するステップS107。また、オブジェクト生成部110は、直線上に設定された点以下、端点と呼ぶから、垂線と直線との交点までの距離を算出し、群化テーブルに格納するステップS109。端点は、例えば図37における370のような位置に設定する。そして、オブジェクト生成部110は、処理対象のパーティクルの位置のデータを用いて、処理対象のパーティクルが直線の右又は左のいずれに存在するかを判断し、判断結果を群化テーブルに格納するステップS111。なお、本実施の形態においては、z軸におけるz座標が正の位置から原点に向かう方向から座標空間を観察した場合を想定しており、観察者にとっての右を右とし、観察者にとっての左を左とする。例えば図37の場合、パーティクル4、5及び1は直線の左であると判定し、パーティクル3及び2は直線の右であると判断する。図38に、第2の実施の形態における群化テーブルに格納されるデータの一例を示す。図38の例では、パーティクル番号と、垂線と直線との交点の位置と、端点からの距離と、左右の判定結果とが格納されるようになっている。左右の判定結果が0である場合にはパーティクルは直線の左に位置し、左右の判定結果が1である場合にはパーティクルは直線の右に位置しているとする。なお、このようなテーブルは各群について生成される。そして、オブジェクト生成部110は、未処理のパーティクルが有るか判断するステップS113。未処理のパーティクルが有ると判断された場合ステップS113:Yesルート、次のパーティクルについて処理するため、ステップS105の処理に戻る。一方、未処理のパーティクルが無いと判断された場合ステップS113:Noルート、処理は端子Iを介して図39のステップS115に移行する。図39の説明に移行し、オブジェクト生成部110は、群化テーブルのデータを、距離及び左右判定に基づきソートするステップS115。例えば、図38に示したデータをソートすると、図40に示すようなデータになる。すなわち、第1行目乃至第3行目のパーティクルは直線の左に位置するパーティクルであり、パーティクル4、パーティクル5、パーティクル1の順に端点に近い位置にある。また、第4行目及び第5行目のパーティクルは直線の右に位置するパーティクルであり、パーティクル3の方がパーティクル2よりも端点に近い位置にある。そして、オブジェクト生成部110は、パーティクル間を結ぶ線分の描画データを生成し、記憶装置図示しないに格納するステップS117。例えば図40の例では、まず直線の左側に位置するパーティクルについて、垂線と直線との交点が隣り合っているパーティクル4とパーティクル5との間の線分の描画データを生成し、また垂線と直線との交点が隣り合っているパーティクル5とパーティクル1との間の線分の描画データを生成する。また、直線の右側に位置するパーティクルについて、垂線と直線との交点が隣り合っているパーティクル3とパーティクル2との間の線分の描画データを生成する。さらに、直線の左側に生成された構造の一方の端であるパーティクル4と、直線の右側に生成された構造の端に位置するパーティクルのうちパーティクル4に近い方のパーティクルであるパーティクル3とを結ぶ線分の描画データを生成する。また、直線の左側に生成された構造のもう一方の端であるパーティクル1と、直線の右側に生成された構造のもう一方の端であるパーティクル2とを結ぶ線分の描画データを生成する。以上のようにして生成された構造に色を付すと、図41に示すような構造が生成される。そして、オブジェクト生成部110は、ステップS117の処理によって生成された構造に対して凹み部分の補正及び角を丸める処理を行い、群表示のデータを生成し、記憶装置図示しないに格納するステップS119。これらの処理は、第1の実施の形態において説明したものと同様である。ステップS119までの処理によって、図42に示すようなデータが生成される。また、色付け処理部112は、速度データ格納部107に格納されている位置のデータを用いて、群番号gの群のz軸方向の位置を特定するステップS120。ここでは、群に属するパーティクルの位置の平均を、群の位置であるとする。そして、色付け処理部112は、z軸方向における群の位置に応じて、群表示に対して色付け処理を実施し、処理結果を記憶装置図示しないに格納するステップS121。色付け処理については、第1の実施の形態において説明したものと同様である。そして、オブジェクト生成部110は、未処理の群が有るか判断するステップS123。未処理の群が有ると判断された場合ステップS123:Yesルート、次の群について処理するため、オブジェクト生成部110は、群番号gをg=g1と設定するステップS125。そして処理は端子Hを介してステップS103に戻る。一方、未処理の群が無いと判断された場合ステップS123:Noルート、流線生成部111は、流線描画処理を実施するステップS127。流線描画処理については、第1の実施の形態において説明したとおりである。そして元の処理に戻る。以上のような処理を実施すれば、適切な形状の群表示を生成することができるようになる。以上本技術の一実施の形態を説明したが、本技術はこれに限定されるものではない。例えば、上で説明した情報処理装置1の機能ブロック図は必ずしも実際のプログラムモジュール構成に対応するものではない。また、上で説明した各テーブルの構成は一例であって、必ずしも上記のような構成でなければならないわけではない。さらに、処理フローにおいても、処理結果が変わらなければ処理の順番を入れ替えることも可能である。さらに、並列に実行させるようにしても良い。なお、上で述べた例では、z軸方向の位置に応じて群表示の色を設定しているが、色を設定する方向はz軸方向に限られるわけではない。座標空間における流れの特定に応じてx軸、y軸又は任意の軸の方向を決定し、決定した軸の方向における位置に応じて色を変更するようにしてもよい。また、上で述べた例では、凹み部分を平滑化する処理及び角を丸める処理を実施しているが、処理時間を短縮する場合には、これらの処理を省略してもよい。また、第2の実施の形態では最小二乗法を利用する例を示したが、主成分分析の手法を利用するようにしてもよい。また、上で述べた例では一台のコンピュータによって処理を実行しているが、複数台のコンピュータによって処理を実行するようにしてもよい。なお、上で述べた情報処理装置1は、コンピュータ装置であって、図43に示すように、メモリ2501とCPU2503とハードディスクドライブHDD2505と表示装置2509に接続される表示制御部2507とリムーバブルディスク2511用のドライブ装置2513と入力装置2515とネットワークに接続するための通信制御部2517とがバス2519で接続されている。オペレーティングシステムOS:OperatingSystem及び本実施例における処理を実施するためのアプリケーションプログラムは、HDD2505に格納されており、CPU2503により実行される際にはHDD2505からメモリ2501に読み出される。CPU2503は、アプリケーションプログラムの処理内容に応じて表示制御部2507、通信制御部2517、ドライブ装置2513を制御して、所定の動作を行わせる。また、処理途中のデータについては、主としてメモリ2501に格納されるが、HDD2505に格納されるようにしてもよい。本技術の実施例では、上で述べた処理を実施するためのアプリケーションプログラムはコンピュータ読み取り可能なリムーバブルディスク2511に格納されて頒布され、ドライブ装置2513からHDD2505にインストールされる。インターネットなどのネットワーク及び通信制御部2517を経由して、HDD2505にインストールされる場合もある。このようなコンピュータ装置は、上で述べたCPU2503、メモリ2501などのハードウエアとOS及びアプリケーションプログラムなどのプログラムとが有機的に協働することにより、上で述べたような各種機能を実現する。なお、図3に示した各処理部は、図43におけるCPU2503及びプログラムの組み合わせ、すなわち、CPU2503がプログラムを実行することにより実現してもよい。また、図3に示した各データ格納部は、図43におけるメモリ2501やHDD2505等として実現してもよい。以上述べた本技術の実施の形態をまとめると以下のようになる。本実施の形態に係る可視化処理は、A拡散場におけるパーティクルの位置及び速度のデータに基づき、パーティクル間の距離と速度ベクトルの差の絶対値とがいずれも所定の基準値以下であるパーティクルが同じグループになるように、パーティクルのグループ分けを行い、Bグループ分けにより得られたグループの各々について、当該グループに属するパーティクルを覆うグループ表示のデータを生成し、Cパーティクルを表示するためのデータにグループ表示のデータを付加した表示データを生成する処理をコンピュータが実行する。このようにすれば、微小時間前に同じ位置にあったと推定されるパーティクルが同じグループに属するようになるので、拡散場におけるパーティクルの動きを容易に把握できるようになる。また、コンピュータが、Dグループ分けにより得られたグループの各々について、当該グループに属するパーティクルの位置のデータを用いて、第1の軸方向における当該グループの位置を特定し、E特定された第1の軸方向におけるグループの位置に応じて、グループ表示の色を変更する処理をさらに実行するようにしてもよい。例えば、深さ方向に一様な流れになっている場合には、視線の方向によっては、流れの三次元的な構造を把握しにくい場合がある。このような場合に、深さ方向の位置に応じてグループ表示の色を変更すれば、流れの三次元的な構造がよりわかりやすくなる。また、コンピュータが、Fグループ分けにより得られたグループの各々について、当該グループに属するパーティクルのうち、当該グループに属するパーティクルの位置の平均に最も近いパーティクルを特定し、G特定されたパーティクルの位置から所定時間後の当該パーティクルの位置までの流線のデータを生成し、当該流線のデータをグループ表示のデータに付加する処理をさらに実行するようにしてもよい。このように、各グループについての流線を表示できるようにすれば、大まかな流れの状況を把握することも容易になる。また、上記グループ表示を生成する処理において、b1グループに属するパーティクルの間を結ぶ線分を、ループを形成しないように、短い順にN-1Nは当該グループに属するパーティクルの数本つなぐことにより、第1の構造のデータを生成し、b2少なくとも第1の構造に含まれる線分の太さを太くする処理と第1の構造に含まれる凹み部分を平滑化する処理とを実行することにより、グループ表示のデータを生成するようにしてもよい。このようにすれば、グループ表示は球体又は楕円体に近い自然な形状になる。また、上記グループ表示を生成する処理において、b3グループに属するパーティクルの位置のデータを用いて、最小自乗法によって直線を特定し、b4特定された直線とグループに属するパーティクルとの位置関係に基づき、グループに属するパーティクルを第1のグループと第2のグループとにグループ分けし、b5第1のグループに属する各パーティクルから直線に引いた垂線が直線と交わる点の位置を特定し、当該垂線が直線と交わる点の位置が隣り合っているという条件を満たすパーティクルの間を線分でつなぐことにより、第2の構造のデータを生成し、b5第2のグループに属する各パーティクルから直線に引いた垂線が直線と交わる点の位置を特定し、当該垂線が直線と交わる点の位置が隣り合っているという条件を満たすパーティクルの間を線分でつなぐことにより、第3の構造のデータを生成し、b6第2の構造における端と第3の構造における端との間を線分でつなぐことにより、第4の構造のデータを生成し、b7少なくとも第4の構造の凹み部分を平滑化する処理を実行することにより、グループ表示のデータを生成するようにしてもよい。このような方法であっても、適切な形状のグループ表示になる。なお、上記方法による処理をコンピュータに行わせるためのプログラムを作成することができ、当該プログラムは、例えばフレキシブルディスク、CD-ROM、光磁気ディスク、半導体メモリ、ハードディスク等のコンピュータ読み取り可能な記憶媒体又は記憶装置に格納される。尚、中間的な処理結果はメインメモリ等の記憶装置図示しないに一時保管される。以上の実施例を含む実施形態に関し、さらに以下の付記を開示する。付記1拡散場における複数のパーティクルを可視化する可視化処理プログラムにおいて、コンピュータに、前記複数のパーティクルのうち、パーティクルの位置及び速度のデータに基づき、パーティクル間の距離が所定距離以下且つ前記パーティクル間の速度ベクトルの差の絶対値が所定値以下であるパーティクルが同じグループになるように、前記複数のパーティクルのグループ分けを行わせ、前記グループ分けされたグループの各々について、グループに属するパーティクルを含むグループ表示データを生成させ、前記グループに属するパーティクルを表示するパーティクル表示データに前記グループ表示データを付加した表示データを生成させることを特徴とする可視化処理プログラム。付記2前記可視化処理プログラムはさらに、前記コンピュータに、前記グループ分けされたグループの各々について、グループに属するパーティクルの位置のデータを用いて、第1の軸方向におけるグループの位置を特定させ、前記特定された第1の軸方向におけるグループの位置に応じた色に、グループ表示の色を設定させることを特徴とする付記1記載の可視化処理プログラム。付記3前記可視化処理プログラムはさらに、前記コンピュータに、前記グループ分けされたグループの各々について、グループに属するパーティクルのうち、グループに属するパーティクルの位置の平均に最も近いパーティクルを特定させ、前記特定されたパーティクルの位置から所定時間後の前記特定されたパーティクルの位置までの流線のデータを生成させ、前記生成された流線のデータを前記グループ表示データに付加させることを特徴とする付記1又は2記載の可視化処理プログラム。付記4前記可視化処理プログラムはさらに、前記コンピュータに、前記グループ分けされたグループの各々について、グループに属するパーティクルの間を結ぶ線分を、ループを形成しないように、短い順にN-1Nは当該グループに属するパーティクルの数本接続させることにより、第1の構造のデータを生成させ、前記第1の構造に含まれる線分の太さを太くさせ、前記第1の構造に含まれる凹み部分を平滑化させることにより、前記グループ表示データを生成させることを特徴とする付記1乃至3のいずれか1つ記載の可視化処理プログラム。付記5前記可視化処理プログラムはさらに、前記コンピュータに、前記グループ分けされたグループの各々について、グループに属するパーティクルの位置のデータを用いて、最小二乗法に基づいた直線を特定させ、前記特定された直線とグループに属するパーティクルとの位置関係に基づき、グループに属するパーティクルを第1のグループと第2のグループとにグループ分けさせ、前記第1のグループに属する各パーティクルから前記直線に引いた垂線が前記直線と交わる点の位置を特定させ、各々から引かれた垂線が前記直線と交わる点の位置が隣り合っているパーティクルの間を線分で接続させることにより、第2の構造のデータを生成させ、前記第2のグループに属する各パーティクルから前記直線に引いた垂線が前記直線と交わる点の位置を特定させ、各々から引かれた垂線が前記直線と交わる点の位置が隣り合っているパーティクルの間を線分で接続させることにより、第3の構造のデータを生成させ、前記第2の構造における端と第3の構造における端との間を線分で接続させることにより、第4の構造のデータを生成させ、前記第4の構造の凹み部分を平滑化させることにより、前記グループ表示のデータを生成させることを特徴とする請求項1乃至3のいずれか1つ記載の可視化処理プログラム。付記6拡散場における複数のパーティクルを可視化する可視化処理方法において、コンピュータが、前記複数のパーティクルのうち、各パーティクルの位置及び速度のデータに基づき、パーティクル間の距離が所定距離以下且つ前記パーティクル間の速度ベクトルの差の絶対値が所定値以下であるパーティクルが同じグループになるように、前記複数のパーティクルをグループ分けし、前記グループ分けされたグループの各々について、グループに属するパーティクルを含むグループ表示データを生成し、前記グループに属するパーティクルを表示するパーティクル表示データに前記グループ表示データを付加した表示データを生成することを特徴とする可視化処理方法。付記7拡散場における複数のパーティクルの位置及び速度のデータに基づき、パーティクル間の距離が所定距離以下且つ前記パーティクル間の速度ベクトルの差の絶対値が所定値以下であるパーティクルが同じグループになるように、前記複数のパーティクルをグループ分けするグルーピング部と、前記グループ分けされたグループの各々について、グループに属するパーティクルを含むグループ表示データを生成するグループ表示生成部と、前記グループに属するパーティクルを表示するパーティクル表示データに前記グループ表示データを付加した表示データを生成する表示データ生成部と、を有することを特徴とする可視化処理装置。