"微小粒子の測定データを保存するデータ格納部と、前記測定データから独立した3種の変数を選択する入力部と、前記3種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す3D立体画像を作成するデータ処理部と、前記3D立体画像を表示する表示部と、を有する3Dデータ解析装置。前記入力部からの入力信号に基づいて、前記表示部において前記3D立体画像を回転、拡大および縮小させて表示する請求項1記載の3Dデータ解析装置。前記入力部からの入力信号に基づいて、前記表示部において前記3D立体画像を回転させ、任意に選択された座標軸方向からの該3D立体画像の立体像観察を可能とする請求項2記載の3Dデータ解析装置。前記図形は、所定形状のポリゴンの組み合わせからなる多面体である請求項3記載の3Dデータ解析装置。前記3D立体画像において、該3D立体画像の立体像観察時に手前に観察される前記図形をより濃く、奥に観察される前記図形をより淡く表示する請求項4記載の3Dデータ解析装置。前記3D立体画像において、前記座標空間を構成する立体形状の各辺に前記座標軸を表示する請求項5記載の3Dデータ解析装置。前記3D立体画像において、前記座標軸のうち、該3D立体画像の立体像観察時に手前に観察される部分をより太く、奥に観察される部分をより細く表示する請求項6記載の3Dデータ解析装置。前記座標軸の少なくとも一以上が、バイエクスポネンシャル軸である請求項7記載の3Dデータ解析装置。前記3D立体画像を前記図形が揺動する動画像により表示し、該動画像の立体像観察時に奥に観察される前記図形に比して手前に観察される図形をより大きく揺動させて表示する請求項8記載の3Dデータ解析装置。前記3D立体画像を前記図形が点滅する動画像により表示し、該動画像の立体像観察時に奥に観察される前記図形に比して手前に観察される図形をより高頻度に点滅させて表示する請求項9記載の3Dデータ解析装置。前記表示部がディスプレイであり、該ディスプレイに表示される前記3D立体画像を立体視するためのメガネを備える請求項10記載の3Dデータ解析装置。請求項1~11に記載のいずれか一項に記載の3Dデータ解析装置と微小粒子測定装置とを連設してなる微小粒子解析システム。微小粒子の測定データから独立した3種の変数を選択する手順と、前記3種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す3D立体画像を作成する手順と、前記3D立体画像を表示させる手順と、を含む3Dデータ解析方法。微小粒子の測定データから選択される独立した3種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す3D立体画像を作成するステップと、前記3D立体画像を表示するステップと、をコンピュータに実行させる3Dデータ解析プログラム上述のように、従来、フローサイトメータに用いられているデータ解析装置では、ゲーティングの際、ユーザは多数のヒストグラムあるいはサイトグラムを参照したり、2つのサイトグラムを組み合わせた立体的な分布図三次元分布図を想起したりしながら、分布図上における解析対象細胞小集団の位置を特定する必要があった。そこで、本発明は、多数のヒストグラムあるいはサイトグラムを参照したり、三次元分布図を想起したりすることなく、解析対象とする微小粒子ならびに微小粒子小集団を分布図上で容易かつ直感的に特定できるデータ解析装置を提供することを主な目的とする。上記課題解決のため、本発明は、微小粒子の測定データを保存するデータ格納部と、前記測定データから独立した3種の変数を選択する入力部と、前記3種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す3D立体画像を作成するデータ処理部と、前記3D立体画像を表示する表示部と、を有する3Dデータ解析装置を提供する。この3Dデータ解析装置によれば、任意に選択した3種のパラメータを座標軸とした三次元分布図を立体視しながら測定データの解析を行うことができる。また、本発明は、微小粒子の測定データから独立した3種の変数を選択する手順と、前記3種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す3D立体画像を作成する手順と、前記3D立体画像を表示させる手順と、を含む3Dデータ解析方法と、微小粒子の測定データから選択される3種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す3D立体画像を作成するステップと、前記3D立体画像を表示するステップと、をコンピュータに実行させる3Dデータ解析プログラムをも提供する。本発明において、微小粒子には、細胞、微生物およびリポソーム等の生体関連粒子、ラテックス粒子、ゲル粒子、工業用粒子等の合成粒子などの微小粒子が広く含まれるものとする。細胞には、動物細胞血球系細胞などおよび植物細胞が含まれる。微生物には、大腸菌等の細菌類、タバコモザイクウイルス等のウイルス類、イースト菌等の菌類などが含まれる。生体関連粒子には、各種細胞を構成する染色体、リポソーム、ミトコンドリア、オルガネラなども含まれる。さらに、生体関連粒子には、核酸やタンパク質、これらの複合体などの生体関連高分子も包含され得るものとする。工業用粒子は、有機もしくは無機高分子材料、金属などであってもよい。有機高分子材料には、ポリスチレン、スチレンジビニルベンゼンおよびポリメチルメタクリレートなどが含まれる。無機高分子材料には、ガラス、シリカおよび磁性体材料などが含まれる。金属には、金コロイドおよびアルミなどが含まれる。これら微小粒子の形状は、一般には球形であるが、非球形であってもよい。また、微小粒子の大きさや質量なども特に限定されない。本発明により、多数のヒストグラムあるいはサイトグラムを参照したり、三次元分布図を想起したりすることなく、解析対象とする微小粒子ならびに微小粒子小集団を分布図上で容易かつ直感的に特定できるデータ解析装置が提供される。以下、本発明を実施するための好適な形態について図面を参照しながら説明する。なお、以下に説明する実施形態は、本発明の代表的な実施形態の一例を示したものであり、これにより本発明の範囲が狭く解釈されることはない。なお、説明は以下の順序で行う。1.3Dデータ解析装置の構成2.3D立体画像の表示3.3D立体画像の特徴3-1図形の形状3-2図形のシェード処理3-3座標軸3-4動画像4.3Dデータ解析プログラム1.3Dデータ解析装置の構成図1に、本発明に係る3Dデータ解析装置の装置構成を示す。ここでは、3Dデータ解析装置を微小粒子測定装置に連設して微小粒子解析システムとして構成した実施形態を示す。また、図2に、この微小粒子解析システムの機能的構成を示す。以下では、微小粒子測定装置としてフローサイトメータを用いる場合を例に説明する。図中符号1で示す3Dデータ解析装置は、通信ケーブル4によりフローサイトメータ2と接続されて微小粒子解析システム3を構成している。3Dデータ解析装置1は、中央処理装置CPU10、メモリ20、ハードディスク30、ユーザインターフェイスなどを含んでいる。ハードディスク30には、3Dデータ解析プログラム31、微小粒子の測定データ32およびオペレーティングシステムOS33などが格納保持されている。ユーザインターフェイスには、ユーザからの情報の入力を受け付けるマウス41およびキーボード42と、ユーザに対して情報を出力するディスプレイ43およびプリンタ44などが含まれる。なお、マウス41およびキーボード42に替えて、あるいはこれらとともに、スティックコントローラおよびペンタブレットなどの入力デバイスを設けてもよい。データ格納部130ハードディスク30は、フローサイトメータ2から出力される微小粒子細胞の測定データ32を保存する。フローサイトメータ2の入出力インターフェイス250から出力される測定データは、通信ケーブル4を介して3Dデータ解析装置1の入出力インターフェイス150に入力され、データ格納部30ハードディスク30に格納される。測定データ32は、データ処理部120において処理される。データ処理部120は、入力部141マウス41あるいはキーボード42などからのユーザの入力を受けて処理を開始する。すなわち、データ処理部120は、ユーザによって測定データ32から独立した3種の変数パラメータが選択され入力されると、選択されたパラメータを座標軸として微小粒子の特性分布を表す三次元分布図を作成する。三次元分布図は、選択されたパラメータを座標軸とした座標空間内に微小粒子をプロットすることにより作成される。微小粒子のプロットは、選択されたパラメータの測定値から各微小粒子の座標空間内の位置と図形を計算し、算出された位置に算出された図形を描画することによって行われる。ここで、独立したパラメータとは、微小粒子の前方散乱光FS、側方散乱SS、蛍光FLおよびインピーダンスなどから選択される互いに異なるパラメータを意味する。蛍光FLは、微小粒子に標識された蛍光色素の波長域毎に異なるパラメータとして扱うことができ、FL1,FL2~FLnnは3以上の整数などと表記される。独立した3種のパラメータとしては、例えば、前方散乱光FS、側方散乱SSおよび蛍光FL1の組み合わせや、前方散乱光FS、側方散乱SSおよびインピーダンスの組み合わせが例示される。この他、独立した3種のパラメータは、測定データから任意に選択して組み合わせることができる。データ処理部120により作成された三次元分布図は、表示部142ディスプレイ43に3D立体画像として表示される。表示部142に表示される3D立体画像は、一または二以上とできる。二以上の3D立体画像を表示する場合、同一の三次元分布図について複数の異なる方向から観察した3D立体画像を表示しても、選択された3種のパラメータのうち少なくとも一つを異にする複数の三次元分布図の3D立体画像を表示してもよい。この3D立体画像は、次に詳しく説明する二眼式ステレオ立体画像とされる。さらに、測定データ32に複数の異なる時刻タイムポイントにおける測定値が含まれる場合には、表示部142は、複数のタイムポイントにおける微小粒子の特性分布を表す三次元分布図を3D立体画像により表示してもよい。複数のタイムポイントにおける測定値を含む測定データとしては、例えば、細胞表面分子複合体の会合若しくは解離を蛍光共鳴エネルギー移動FRET法を用いて経時的に測定したデータ、細胞膜の変化を細胞膜の電荷に対応して蛍光波長が変化する蛍光色素を用いて経時的に測定したデータ、又は細胞表面分子の発現強度を細胞内カルシウムの流入レスポンスと相関させて測定したデータなどがある。複数のタイムポイントにおける三次元分布図の3D立体画像は、並べて一度に表示させてもよく、あるいは切り換えながら一つずつ表示させてもよい。3D立体画像の切り換え表示を行う場合、切り換えは自動的に行われてもよく、ユーザの入力信号に基づいて行われてもよい。複数のタイムポイントにおける三次元分布図の3D立体画像を表示させることで、ユーザが微小粒子の特性分布の経時的変化を確認しながらデータ解析を行うことができ、3種のパラメータ座標軸に時間時間軸を加えたより多元的な解析が可能となる。表示部142への3D立体画像の表示は、入力部141マウス41あるいはキーボード42などからのユーザの入力信号に基づき任意に回転、拡大あるいは縮小させて行ってもよい。また、入力部141からの入力信号に基づいて三次元分布図の座標空間内にゲーティングのための分画領域が設定される場合には、3D立体画像は、その3D立体画像中に表示された分画領域を示す立体形状とともに回転、拡大あるいは縮小される。フローサイトメータ2は、従来公知の装置と同様の構成あるいはこれを適宜改変した構成とでき、具体的には、制御部210、フロー系220、検出系230、入出力インターフェイス250等から構成される。フロー系220は、フローセルやマイクロチップに形成された流路内において、微小粒子を含むサンプル液層流をシース液層流の中心に流し、層流中に微小粒子を一列に配列させる。検出系230は、流路を通流する微小粒子の特性を示すパラメータ値を取得する。具体的には、光学検出部231は、通流する微小粒子に光を照射し、微小粒子から生じる散乱光や蛍光などを検出し、その強度を取得する。光学検出部231は、レーザー光源、レンズ、ミラー、フィルター、CCDおよびCMOS素子等のエリア撮像素子あるいはPMTphotomultipliertubeなどを含んでなる。また、電気検出部232は、通流する微小粒子に対向して配された電極を含んでなり、微小粒子のインピーダンス、容量値キャパシタンスおよびインダクタンスなどを取得する。フローサイトメータ2は、解析の結果所望の特性を有すると判定された微小粒子を分取するための分取系240を備えていてもよい。分取系240には、例えば、微小粒子を含む液滴をフローセル外の空間に吐出し、液滴の移動方向を制御することにより、所望の微小粒子のみを容器に回収する方式を採用できる。検出系230において検出された散乱光および蛍光などの強度の測定値、あるいはインピーダンス、容量値キャパシタンスおよびインダクタンスなどの測定値は、電気信号に変換され、測定データとして入出力インターフェイス250から出力される。2.3D立体画像の表示図3に、本発明に係る3Dデータ解析装置により表示される三次元分布図を模式的に示す。この三次元分布図は、表示部142に3D立体画像により表示され、ユーザによって立体的に視認され得るものである。三次元分布図5は、ユーザにより選択された3種のパラメータを座標軸とした座標空間6内に微小粒子の特性分布を示す。三次元分布図5では、選択されたパラメータの測定値から算出される位置に各微小粒子に対応する図形7が描画される。図では、3種のパラメータを、前方散乱光FSLin:X軸、側方散乱SSLin:Y軸および第一の蛍光FL1Lin:Z軸の組み合わせとした場合を例示した。各座標軸にとるパラメータは、任意に選択される組み合わせとできる。例えば、X軸に第一の蛍光FL1、Y軸に第二の蛍光FL2、Z軸にインピーダンスをとることもできる。三次元分布図5の3D立体表示は、二眼式ステレオ立体画像により行われる。図4に、本発明に係る3Dデータ解析装置により表示される二眼式ステレオ立体画像を模式的に示す。データ処理部120は、ユーザによってパラメータが選択されると三次元分布図5を作成して、該分布図を左眼で見たときの画像左眼用画像5Lと、右眼で見たときの画像右眼用画像5Rを作成する。表示部142ディスプレイ43は、左眼用画像5Lおよび右眼用画像5Rを同時に表示し、左眼用画像5Lを左眼だけに、右眼用画像5Rを右眼だけに分離提示する。分離提示は、例えばメガネ方式の一つである時分割方式では、左眼用画像5Lと右眼用画像5Rとをわずかな時間差で交互に表示させ、これにシャッタメガネ8を同期させることによって行うことができる。この他、分離提示は、アナグリフ方式および偏光フィルタ方式などの他のメガネ方式、パララックスバリア方式およびレンティキュラ方式などの裸眼方式、ステレオスコープ方式およびヘッドマウント方式などのビューア方式を採用してもよい。ディスプレイ43は、左眼用画像5Lおよび右眼用画像5Rを分離提示することによって、三次元空間で三次元分布図を見ている際に目に映る映像を再現し、ユーザに分布図を立体視させる。このように、3Dデータ解析装置1では、ユーザは、任意に選択した3種のパラメータを座標軸とした三次元分布図を立体視しながら測定データの解析を行うことができる。このため、3Dデータ解析装置1では、分布図上において解析対象とする微小粒子ならびに微小粒子小集団を容易かつ直感的に特定することができ、従来のように多数のヒストグラムあるいはサイトグラムを参照したり、三次元分布図を想起したりする必要がない。また、座標軸に用いるパラメータを任意に組み合わせて三次元分布図を表示させることで、一つのグラフで微小粒子の3つの特性に関する情報を得ることができる。さらに、同一の三次元分布図について複数の異なる方向から観察した3D立体画像を表示させたり、選択された3種のパラメータのうち少なくとも一つを異にする複数の三次元分布図の3D立体画像を表示させたりすることで、より多くの情報を得ることができる。このため、3Dデータ解析装置1では、従来のヒストグラムあるいはサイトグラムによる表示に比べて、参照すべきグラフ数を減らして効率的な解析を行うことができる。3.3D立体画像の特徴以下、本発明に係る3Dデータ解析装置により表示される3D立体画像の特徴について順に説明する。3-1図形の形状図3中符号7で示した微小粒子に対応する図形は、所定形状のポリゴンの組み合わせからなる多面体として計算され、3D立体画像中に表示される。上記のように、データ処理部120は、ユーザが選択したパラメータの測定値に基づいて各微小粒子の座標空間内の位置と図形7を計算し、三次元分布図を作成する。この際、図形7を所定形状のポリゴンの組み合わせからなる多面体として計算することで、データ処理部120における計算負荷を軽減できる。また、図形7を所定形状のポリゴンの組み合わせからなる多面体として3D立体画像中に表示することで、立体視した際の画像の立体感を増強することができる。所定形状のポリゴンの組み合わせからなる多面体としては、例えば、図5のAに示すような、三角形のポリゴンが6つ組み合わされてなる6面体や、図5のBに示す該ポリゴンが8つ組み合わされてなる8面体を採用できる。図形7の形状は、所定形状のポリゴンの組み合わせからなる多面体である限りにおいて特に限定されないが、計算負荷の軽減と立体感の観点から6面体あるいは8面体が好ましい。3-2図形のシェード処理3D立体画像において、図形7は、立体視した際に手前に観察される図形ほど濃く、奥に観察される図形ほど淡く表示される。このように図形7の濃淡を変化させる処理を、以下シェード処理と称するものとする。シェード処理された図形7の立体観察像以下単に立体像というの概念図を図6に示す。図中矢印の方向に従って、手前に観察される図形7ほど濃く、奥に観察される図形7ほど淡くなっている。このように、図形7のシェード処理を行うことによって、3D立体画像の立体像に奥行きを与えて立体感を向上できる。図7を参照して、シェード処理の処理方法を説明する。ディスプレイ43には左眼用画像と右眼用画像が同時に表示され、立体視した際にディスプレイ43の画面の位置に観察される図形70の左眼用画像と右眼用画像は重なった状態で表示される図7のB参照。ディスプレイ43に表示される左眼用画像が右眼用画像よりも右側にある場合図7のA参照、図形は、ディスプレイ43の画面の位置よりも手前に立体視される。画面位置から飛び出して観察される図形の立体像を図中符号71で示し、ディスプレイ43に表示される図形71の左眼用画像を符号71Lで、右眼用画像を符号71Rで示す。一方、ディスプレイ43に表示される左眼用画像が右眼用画像よりも左側にある場合図7のC参照、図形は、ディスプレイ43の画面の位置よりも奥に立体視される。画面位置から飛び出して観察される図形の立体像を図中符号72で示し、ディスプレイ43に表示される図形71の左眼用画像を符号72Lで、右眼用画像を符号72Rで示す。シェード処理では、手前に観察される図形71の左眼用画像71Lおよび右眼用画像71Rをより濃く、奥に観察される図形72の左眼用画像72Lおよび右眼用画像72Rをより淡く表示するようにする。3-3座標軸3D立体画像において、座標軸は、立体視した際に手前に観察される部分ほど太く、奥に観察される分ほど細く表示される。太さを変化させた座標軸の立体像の概念図を図8に示す。このように、座標軸の太さを変化させることによって、3D立体画像の立体像に奥行きを与えて立体感を向上できる。また、図8に示すように、座標軸の目盛間隔を立体視した際に手前に観察される部分ほど広く、奥に観察される分ほど狭く表示することによって、立体像にさらに奥行きを与えられる。また、座標軸名図ではSSLinおよび目盛数値図では2004006008001000の文字を手前ほど大きく、奥ほど小さく表示したりすることによっても同様の効果が得られる。なお、座標軸の太さ、目盛間隔および文字の大きさを変化させる処理は、上述のシェード処理を応用して行うことが可能である。座標軸は、線形軸リニア軸と対数軸とが組み合わされた特性を有するバイエクスポネンシャル軸としてもよい。このバイエクスポネンシャル軸の詳細は、例えば、ANewLogicleDisplayMethodAvoidsDeceptiveEffectsofLogarithmicScalingforLowSignalsandCompensatedDataCytometryPartA69A5415512006に記載されている。バイエクスポネンシャル軸では、座標軸として選択されたパラメータの測定値が所定値より小さいデータに対しては、線形関数を主たる関数要素とする関数を適用して微小粒子に対応する図形7の位置が算出される。また、測定値が所定値より大きいデータに対しては、対数関数を主たる関数要素とする関数を適用して図形7の位置を算出する。より簡便には、バイエクスポネンシャル軸は、所定値よりも大きい領域を対数軸、所定値よりも小さい領域を線形軸としたものとできる。三次元分布図の座標軸をバイエクスポネンシャル軸とすることで、対数軸の特性を活かした広いダイナミックレンジの表示が可能となり、同時に線形軸の特性により負数の表示も可能となる。なお、三次元分布図の座標軸は、少なくとも一つがバイエクスポネンシャル軸とされればよいものとする。3-4動画像既に説明したように、表示部142ディスプレイ43への3D立体画像の表示は、入力部141マウス41あるいはキーボード42などからのユーザの入力信号に基づき任意に回転、拡大あるいは縮小させて行ってもよい。3D立体画像を回転する際には、図3に示したように、座標空間6を構成する立体形状図では立方体の各辺に座標軸が表示されていることが好ましい。これらの座標軸により座標空間6の立体形状が明確になるため、3D立体画像を回転させた際の三次元分布図の向きの変化を認識し易くなる。ディスプレイ43に表示される3D立体画像は、ユーザからの入力により随意に回転される他、常時一定方向あるいは不特定方向にゆっくりと回転していてもよい。3D立体画像を常時回転する動画像により表示することで、静止画像による表示に比べて立体感を高められる。さらに、ディスプレイ43に表示される3D立体画像は、ユーザによる回転操作中及び自動回転中の任意のタイミングで、ユーザの入力信号に基づいて、ユーザにより選択された座標軸方向からの立体観察像を提供する向きにまで自動的に回転される。図9に、各座標軸方向からの三次元分布図の立体観察像を示す。図9のAはZ軸方向からの観察像、図9のBはX軸方向からの観察像、図9のCはY軸方向からの観察像を示す。各座標軸方向からの視点の切り換えは、例えばキーボード42からのZキーの入力によって画像がZ軸方向視点に回転し、Xキーの入力によって画像がZ軸方向視点からX軸方向視点へ回転するようにすればよい。また、各座標軸方向からの視点の切り換えは、例えばディスプレイ43に表示されたアイコンをマウス41によってクリックすることによって行われるようにしてもよい。このように、簡単な入力によって各座標軸方向からの視点を切り換えて3D立体画像を観察できるようにすることで、三次元分布図における微小粒子の特性分布をユーザが把握し易くなる。なお、ディスプレイ43において3D立体画像を常時回転させて表示する場合、3D立体画像は、三次元分布図の上下方向を保つように回転されることが好ましい。すなわち、3D立体画像は、三次元分布図のXY平面、YZ平面及びZX平面から選択されるいずれかの面を常に分布図下方に向けた状態で回転されることが好ましい。具体的には、例えば図9のAに示す3D立体画像を常時回転させる場合、ZX平面が三次元分布図の底面に常に位置するように画像を回転させる。この際、三次元分布図の回転軸を傾けたり、傾き角度を変化させたりしながら、画像を回転させてもよい。このように、3D立体画像の回転方向に一定の制限を加えたほうが、ユーザが三次元分布図に対する自身の視点方向を知覚し易く、三次元分布図の向きを把握できなくなるのを防止できる。ディスプレイ43に表示される3D立体画像は、微小粒子に対応する図形が揺動する動画像により表示されてもよい。このとき、立体視した際に奥に観察される図形に比して手前に観察される図形を大きく揺動させて表示させる。揺動動作が付与された図形の立体像の概念図を図10に示す。図形71,72は、図中矢印に示すように左右に揺動し、左右への揺動幅は手前に観察される図形71で大きく、奥に観察される図形72で小さくなっている。このように、立体視した際に奥に観察される図形に比して手前に観察される図形をより大きく揺動させて表示することで、3D立体画像の立体像に奥行きを与えて立体感を向上できる。また、3D立体画像を動画像により表示する場合、微小粒子に対応する図形が点滅するようにしてもよい。このとき、立体視した際に奥に観察される図形に比して手前に観察される図形をより高頻度に点滅させて表示させることで、3D立体画像の立体感をさらに向上できる。さらに、測定データ32に複数のタイムポイントにおける測定値が含まれる場合には、各タイムポイントにおける三次元分布図の3D立体画像を動画像により表示させることもできる。これにより、例えば上述の細胞表面分子複合体の会合又は解離を測定する例においては、細胞表面分子複合体の会合等の経時的変化を動画像によって確認することができる。以上のように、本発明に係る3Dデータ解析装置は、表示される3D立体画像の立体感を向上させるための工夫がなされている。従って、点微小粒子に対応する図形と線座標軸のみからなる三次元分布図であっても、ユーザはその立体像を良好に視認しながら測定データの解析を行うことができ、分布図上において解析対象とする微小粒子ならびに微小粒子小集団を容易かつ直感的に特定できる。4.3Dデータ解析プログラム本発明に係る3Dデータ解析プログラムは、微小粒子の測定データから選択される3種の変数を座標軸とする座標空間内における位置と図形を計算し、微小粒子の特性分布を表す3D立体画像を作成するステップと、前記3D立体画像を表示するステップと、をコンピュータに実行させる。図1および図2を再度参照して上述した実施形態に基づき説明すると、3Dデータ解析プログラムは、ハードディスク30に格納保持される図中符号31参照。3Dデータ解析プログラムは、CPU10およびオペレーティングシステムOS33の制御の下でメモリ20に読み込まれて、データ処理部120における三次元分布図の3D立体画像の作成処理と、表示部142への3D立体画像の表示処理を実行する。3Dデータ解析プログラムは、コンピュータ読み取り可能な記録媒体に記録されたものとできる。記録媒体としては、コンピュータ読み取り可能な記録媒体であれば特に制限はないが、例えば、フレキシブルディスクやCD-ROM等の円盤形記録媒体が用いられる。また、磁気テープ等のテープ型記録媒体を用いてもよい。"