JPH11258144A - 粒子画像解析装置 - Google Patents

Abstract

(57)【要約】 【課題】 粒子撮像画像から粒子の体積情報を算出す
る。 【解決手段】 同一粒子に対して直交する２つの方向か
ら粒子を撮像し、各画像中の粒子領域に適合する幾何学
形状を見つけ出し、その適合幾何学形状に基いて粒子ご
とに体積情報を算出し、体積情報の分布パターンを表示
する。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、細胞、パウダー、
粉体等の粒子を撮像し、解析し、粒子の大きさに関する
情報を提供する粒子画像解析装置に関する。

【０００２】

【従来の技術】ファインセラミックス粒子、顔料、化粧
品用パウダーやトナー粒子等各種粉体の品質を管理する
上で、その粒子の形状やサイズを分析、測定することは
重要なことである。測定装置として、粒子をフローセル
に流しその粒子像を撮像し解析する粒子画像解析装置
（イメージフローサイトメータ）が知られている。しか
しながら、従来、この種の装置では、撮像粒子の粒子形
状と投影面積の情報は得られるが、粒子の体積情報は得
られないという問題があった。

【０００３】その問題を解決するため、（ａ）特開平９
−１７８６４４号公報では、電気的検知帯法を組み合わ
せることにより粒子像の奥行寸法情報を得ようとしてい
る。（ｂ）特開平１−２９６１３６号公報では、直交す
る方向から粒子を撮像し各粒子形状が共に円形であれ
ば、撮像粒子が球形であると判断している。

【０００４】

【発明が解決しようとする課題】上記（ａ）では、粒子
ごとに画像情報と電気的検知情報との対応がとれていな
いため、推定される奥行寸法情報としてはかなり幅のあ
るものとなっていた。そこで粒子体積についてもっと確
かな情報を提供できる装置が求められていた。上記
（ｂ）は、粒子が球形か否かを判定するものであり、こ
れも満足のいくものではなかった。ところで、（ｃ）特
開昭６０−５２７４１号公報では、粒子画像にマッチン
グした仮想円の情報から粒子径を求めることが開示され
ているが、粒子が円形ではないものについては無力であ
り、また粒子体積を算出することは開示されていない。
本発明は、より確かな粒子体積情報を提供することがで
きる粒子画像解析装置を提供することを課題とする。

【０００５】

【課題を解決するための手段】本発明は、粒子撮像画像
から粒子の体積に関する情報を得るために、同一粒子を
交差する２つの方向から撮像し、各画像中の粒子領域に
適合する幾何学形状を求め、その適合幾何学形状から得
た情報に基いて粒子ごとに体積情報を算出するように構
成した。

【０００６】本発明の粒子画像解析装置は、粒子含有試
料液をシース液で囲んで流すフローセルと、その流れ粒
子を撮像する粒子撮像系と、撮像された粒子画像を処理
し解析する解析手段と、解析結果等を表示する表示手段
とを備えている。粒子撮像系は２系列の粒子撮像系を含
んでおり、第１の粒子撮像系は粒子をある面側から見た
第１の画像を得、第２の粒子撮像系は同一粒子を上記あ
る面と交差する面側から見た第２の画像を得る。解析手
段は、撮像粒子ごとに、第１、第２の画像中の各粒子領
域の形状に適合する幾何学形状を決定し、それら適合幾
何学形状の寸法に基いて撮像粒子の体積情報パラメータ
の値を算出し、さらに、粒子ごとに算出された上記体積
情報パラメータの値を集計しその分布パターンを生成す
る。表示装置は、上記解析結果を表示する。

【０００７】

【発明の実施の形態】本発明が分析対象とする粒子とし
ては、ファインセラミックス、顔料、化粧品パウダーや
トナーのような無機物の粉体、エマルジョン、食品添加
物のような有機物の粉体、血球や細胞のような生体粒子
などのあらゆる粒子がその対象となる。粒子は染色処理
が施されていてよい。

【０００８】粒子撮像をより確実にかつ効率的に行うた
めに、フローセルの粒子撮像領域の上流側に粒子の到来
を検出する粒子検出系を設け、粒子検出に基づき粒子撮
像系で粒子撮像するよう構成するのが好ましい。粒子検
出系としては光学的、電気的にいろいろな実現形態が考
えられるが、一例として粒子の流れ領域に連続光を光照
射する光源と光照射領域に到達した粒子からの散乱光や
蛍光などの光を受光する光検出器との組合せを用いて実
現することができる。

【０００９】粒子検出系で生成された粒子検出信号に基
づき、第１、第２の粒子撮像系の光源をそれぞれ発光す
るよう制御することにより、粒子の各静止画像を確実に
かつ効率よく得ることができる。第２の粒子撮像系は第
１の粒子撮像系のわずか下流側に配置することが配置の
容易性の点で好ましく、また、パルス光発光のタイミン
グをずらすことができることから相互光干渉の低減の点
で好ましい。

【００１０】第１、第２の粒子画像は同一粒子を異なる
面側から撮像した画像であり、その方向は９０度に近い
程好ましい。

【００１１】撮像画像を解析する際に使用する幾何学形
状としては、長方形（正方形を含む）、楕円（円を含
む）などの比較的単純な形状である方が演算処理をする
上で好ましい。

【００１２】粒子の体積は、決定された２つの幾何学形
状を投影形状とする幾何学立体の体積を算出することに
よって得られる。例えば、幾何学形状をいずれも長方形
する場合に算出されるのは直方体の体積であり、幾何学
形状をいずれも楕円とする場合は楕円球体の体積であ
り、幾何学形状が長方形と楕円の場合は楕円柱体の体積
である。幾何学形状としては、対象とする粒子の形状に
よく適合すると予想されるものを選べばよいが、予期で
きない場合や不定形の場合には楕円とするのが無難であ
る。また、複数の幾何学形状の中から最もよく粒子画像
に適合するものを選択するようにすれば、体積情報をよ
り精度良く算出することができる。

【００１３】幾何学形状を粒子領域の形状に適合させる
ことについて説明する。適合幾何学形状は、各撮像画像
の各粒子領域について具体的に得られた長さや面積など
の実測値を用いて決定することができる。

【００１４】図５、６を参照して具体的に説明する。適
合幾何学形状が楕円の場合、粒子領域の粒子の流れ方向
の長さを一方の径長Ｌとし、第１の画像の粒子領域の面
積Ｓ１と同じ面積の楕円を第１の適合楕円とすることが
できる。同様に、粒子領域の粒子の流れ方向の長さを一
方の径長Ｌとし、第２の画像の粒子領域の面積Ｓ２と同
じ面積の楕円を第２の適合楕円とすることができる。第
１、第２の適合楕円の径長から楕円球体の３つの径長
Ｌ，Ｘ１，Ｘ２が判明すれば、楕円球体の体積ＶをＶ＝
π・Ｌ・Ｘ１・Ｘ２／６で算出することができる。

【００１５】長方形の場合には、同様に、粒子領域の粒
子の流れ方向の長さを一方の辺長Ｌとし、第１、第２の
粒子領域の面積Ｓ１，Ｓ２とそれぞれ同じ面積の長方形
を第１、第２の適合長方形とすることができる。第１、
第２の適合長方形の辺長から直方体の３つの辺長Ｌ，Ｙ
１，Ｙ２が判明すれば、直方体の体積ＶをＶ＝Ｌ・Ｙ１
・Ｙ２で算出することができる。

【００１６】撮像粒子の体積情報パラメータとしては、
上記算出された立体の体積Ｖそのものでもよいし、体積
の値Ｖを変換することによって得られるパラメータｆ
（Ｖ）であってもよい。例えば、体積球相当径ｄｖで表
わしてよい。これは対象粒子の体積Ｖと同体積の球の直
径であり、（６・Ｖ／π）の３乗根で表わされる。

【００１７】解析手段は、このようにして撮像粒子の体
積情報パラメータの値を算出し、集計しその体積情報パ
ラメータの分布パターンを生成する。表示手段はその分
布パターンを表示する。分布パターンとはそのパラメー
タの値の分布の状態がわかるように表現したものであっ
て、横軸をパラメータ、縦軸をそのパラメータの各値に
おける頻度として表わしたヒストグラム（１パラメータ
分布図）を用いることができる。

【００１８】さらに、解析手段で、粒子ごとに上記パラ
メータを第１パラメータとしその第１パラメータ以外の
第２パラメータの値を算出してもよい。そのとき、表示
手段は、第１パラメータ、第２パラメータの各分布パタ
ーンを表示してよい。あるいは、第１パラメータと第２
パラメータとで規定される２次元パラメータ空間の各点
における頻度を表現したスキャッタグラム（２パラメー
タ分布図）を分布パターンとして表示してもよい。

【００１９】第２パラメータは、粒子の何らかの特性を
表わすものであって、目的や必要に応じて適宜定義し算
出すればよい。一例として、粒子がどれだけ球に近いか
を表わす球形度を反映するパラメータが考えられる。球
形度は各幾何学形状の径長比あるいは円形度をそれぞれ
乗じたものとして算出することができる。あるいは第２
パラメータは粒子の投影面積を反映するパラメータであ
ってよい。そのパラメータとしては粒子領域の面積Ｓそ
のものでもよいし、面積Ｓを変換することによって得ら
れるパラメータｇ（Ｓ）であってもよい。例えば、面積
円相当径ｄｓで表わしてよい。これは対象粒子の投影面
積と同面積の円の直径であり、（４・Ｓ／π）の２乗根
で表わされる。

【００２０】

【実施例】第１図は本発明の一実施例のブロック図であ
る。３は粒子含有試料液をシース液で囲んで流すフロー
セルである。そのフローセル３に、粒子をある方向から
見た第１の画像を得るための第１の粒子撮像系３１と、
同じ粒子を上記ある方向とは実質的に９０度異なる方向
から見た第２の画像を得るための第２の粒子撮像系３２
が配置されている。本実施例においてはさらに粒子検出
系３４も配置されている。

【００２１】図２はフローセル３と関係づけられる光学
系部分、すなわち粒子検出系３４、第１、第２の粒子撮
像系３１、３２部分の平面図である（粒子は紙面表裏方
向に流れる）。第３図はフローセル付近の概略斜視図で
ある（ノズル４０から吐出された粒子４２は、シース液
に包まれながら紙面下から上へ流れる）。

【００２２】粒子検出系３４は、連続発光光源（波長５
３２ｎｍのグリーン半導体レーザ）１、粒子から前方方
向に発せられた光（前方散乱光）を検出するための第１
の光検出器（フォトダイオード）７及び粒子から側方方
向に発せられた光（側方蛍光）を検出するための第２の
光検出器（光電子増倍管）１１とを含む。

【００２３】第１の粒子撮像系３１は、第１のパルス光
源（波長７８０ｎｍの近赤外半導体レーザ）１８と第１
のビデオカメラ２３とを含む。第２の粒子撮像系３２
は、第２のパルス光源（波長７８０ｎｍの近赤外半導体
レーザ）２と第２のビデオカメラ２３とを含む。各光源
１、１８、２からの光２７、２８、２９はそれぞれ粒子
の流れ方向に対して直交する方向から流れ粒子に照射さ
れる。光源１、１８からの光２７、２８の方向と光源２
からの光２９の方向とは９０゜異なる。

【００２４】本実施例では図３に示すように、第１のパ
ルス光源１８からの光Ｌ２は連続発光光源１からの光ｗ
７の照射位置より例えば０．０２ｍｍ程度下流側に位置
するよう、また、第２のパルス光源２からの光２９はさ
らに例えば０．０２ｍｍ程度下流側に位置するように各
光学要素が配置される。

【００２５】連続発光光源１は、コリメートレンズ２２
で平行光に変換された後、ダイクロイックミラー２１で
反射され、コンデンサレンズ４によって細く絞られフロ
ーセル３のセンターを流れる試料流に向けて照射され
る。この光照射領域に該当する位置Ａ１に粒子が流れて
くると、その粒子からの前方散乱光、側方蛍光がそれぞ
れ集光レンズ５，８によって集光され、それぞれコンデ
ンサレンズ６，１０を経て第１光検出器７、第２光検出
器１１にて受光される。

【００２６】光検出器７、１１にて粒子個々に検出され
た信号は信号処理部３６にて特徴パラメータが抽出され
る（特徴パラメータは信号のピーク値やパルス幅などで
あるが、ここでは後者とした）。信号処理部３６で得ら
れた特徴パラメータは撮像制御部３７へ送られる。特徴
パラメータが条件に合致するものであった場合、撮像ト
リガ信号を発生させ第１、第２の粒子撮像系を駆動し、
位置Ａ１で検出された粒子が位置Ｂ１、Ｂ２に到来した
際にそれぞれ粒子撮像する。すなわち、予め撮像対象と
して指定された粒子だけを選別して撮像できるように制
御する。そのため、ある粒子の特徴パラメータが、撮像
対象としている種類の粒子の特徴パラメータに合致して
いるかどうかを実時間で比較し、その結果、撮像対象と
している粒子であると判定されれば、粒子撮像のための
撮像トリガ信号を生成しパルス光源１８、２に供給す
る。トリガ信号Ｔｓ１は位置Ａ１で検出された粒子が位
置Ｂ１に到達したときにちょうど第１パルス光がパルス
発光するよう発せられる。トリガ信号Ｔｓ２は位置Ａ１
で検出された粒子が位置Ｂ２に到達したときにちょうど
第２パルス光がパルス発光するよう発せられる。

【００２７】パルス光源１８、２は、それぞれ撮像制御
部３７からのトリガ信号Ｔｓ１、Ｔｓ２によって、一瞬
だけ（数十ナノ秒程度）発光するタイプの光源であり、
試料流の流速が数ｍ／秒と高速であっても、位置Ｂ１，
Ｂ２の粒子をブレ無く撮像することができる。また、位
置Ｂ１、Ｂ２間の距離は非常に短いのでその間に粒子が
回転して方向を変えるということもない。

【００２８】第２のビデオカメラ１７のビデオ同期信号
を第１のビデオカメラ２３のビデオ同期信号に対し粒子
が位置Ｂ１からＢ２に移動するのに要する時間（例えば
４μｓｅｃ）だけ位相をずらせておく。このことによ
り、第１、第２のビデオカメラで同一粒子をそれぞれ同
一のビデオタイミングで撮像することができ、同一粒子
であるとの対応づけが容易になる。

【００２９】パルス光源１８、２のパルス光は、それぞ
れ光ファイバー１８ａ，２ａ、コリメートレンズ１９、
１２、コールドミラー２０、１３、コンデンサレンズ
４、１４によって細く絞られて試料流に照射される。光
ファイバーを使用することにより、パルスレーザ光のコ
ヒーレンス性が落ち、回折縞の少ない粒子画像を撮像す
ることができる。

【００３０】試料流を通過したパルス光２８，２９は、
集光レンズ５、８、ダイクロイックミラー２６、９、バ
ンドパスフィルター２５、１６、投影レンズ２４，１５
によってそれぞれ第１ビデオカメラ２３、第２ビデオカ
メラ１７の受光面に結像され、同一粒子の透過光像が撮
像される。バンドパスフィルター２５、１６は、例え
ば、連続発光光源１から発せられる可視光領域の光を遮
り、各パルス光源からの光のみを検出する目的で使用さ
れる。ビデオカメラ２３、１７からのビデオ信号は、解
析手段としての機能を果たす画像処理部３８に送られ、
ディジタル画像としてそれぞれ記憶され、画像処理され
る。

【００３１】図４、５、６、７を参照しながら画像処理
について説明する。 〔１〕．第１、第２の画像の獲得 粒子４２を粒子の流れ方向と直交するある面側から撮像
し粒子４２が写っている第１の画像を得、同一粒子を上
記ある面側と直交する面側から撮像して粒子４２が写っ
ている第２の画像を得、それぞれ画像処理部３８の画像
メモリに格納する。 〔２〕．粒子領域の抽出 １セットの粒子画像、すなわち同一粒子に対して得られ
た第１、第２の画像中の各粒子領域のエッジを抽出し二
値化し、それぞれ第１の画像中の粒子領域４４、第２の
画像中の粒子領域４６を抽出する。 〔３〕．粒子領域における各種パラメータの算出 抽出された第１、第２の粒子領域４４、４６の各画素数
を実測し各領域についてそれぞれ円形度Ｅ１、Ｅ２、面
積Ｓ１、Ｓ２、長さ（流れ方向の径）Ｌ１、Ｌ２、幅
（流れと直交する方向の径）Ｗ１、Ｗ２を算出する。 〔４〕．幾何学形状の決定 適合幾何学形状を楕円とする場合は次の通りである。長
さＬ（＝（Ｌ１＋Ｌ２）／２）で幅Ｘ１（＝４・Ｓ１／
π・Ｌ）の楕円４５を第１の粒子領域４４に対する適合
幾何学形状とする。長さＬ（＝（Ｌ１＋Ｌ２）／２）で
幅Ｘ２（＝４・Ｓ２／π・ａ）の楕円４７を第２の粒子
領域４６に対する適合幾何学形状とする。 〔５〕．粒子体積情報パラメータ（体積球相当径）の算
出 第１、第２の適合幾何学形状の寸法Ｌ，Ｘ１，Ｘ２を用
いて粒子の体積Ｖ（＝π・Ｌ・Ｘ１・Ｘ２／６）および
体積球相当径ｄｖ（（６・Ｖ／π）の３乗根）算出す
る。 〔６〕．他のパラメータ（球形度）の算出 上記算出された第１の粒子領域４４の円形度Ｅ１と第２
の粒子領域４６の円形度Ｅ２とを乗じて球形度ｋとする
（ｋ＝Ｅ１・Ｅ２）。 〔７〕．パラメータの集計 他の粒子について、上記処理を繰り返し、上記粒子体積
情報パラメータおよび他のパラメータを記憶し、分布パ
ターンを生成する。 〔８〕．解析結果の表示 生成された分布パターンをＣＲＴに表示する。図７〜１
０はその一例である。

【００３２】図７、８は球形に近いポリスチレン粒子を
測定したときのスキャッタグラム（２パラメータ分布
図）、ヒストグラム（１パラメータ分布図）である。体
積球相当径ｄｖが揃っており、球径度ｋも１．０に近い
ことがわかる。

【００３３】図９、１０はトナー粒子をを測定したとき
のスキャッタグラム、ヒストグラムである。体積球相当
径ｄｖがばらつき、球径度ｋも１．０からばらついてい
ることがわかる。

【００３４】このように体積情報パラメータ以外に第２
のパラメータも算出し分布図を作成すればより詳細な粒
子特性を知ることができる。

【００３５】ところで、上記〔４〕幾何学形状の決定工
程において、体積情報をより精度高く算出するために、
楕円と長方形の中からより粒子領域に適合する形状を選
択するようにすることができる。その場合には、下記の
ように各幾何学形状に対して何らかの評価値を算出する
ようにしその値に基いて定ｋ号幾何学形状として楕円か
長方形かを選ぶようにする。 〔４−１〕．評価値の算出 適合幾何学形状の評価値Ｑを、その適合幾何学形状の幅
と実際の粒子領域の幅との偏差とする。第１の粒子領域
において楕円を適合幾何学形状をした場合の評価値Ｑ１
をその楕円４５の幅Ｘ１（＝４・Ｓ１／π・Ｌ）と実際
の第１の粒子領域４４の幅Ｗ１との偏差とする。すなわ
ちＱ１＝｜Ｘ１−Ｗ１｜。第１の粒子領域において長方
形を適合幾何学形状をした場合の評価値Ｒ１をその長方
形４８の幅Ｙ１ （＝Ｓ１／π・Ｌ）と実際の第１の粒
子領域４４の幅Ｗ１との偏差とする。すなわちＲ１＝｜
Ｙ１−Ｗ１｜。第２の粒子領域において楕円を適合幾何
学形状をした場合の評価値Ｑ２をその楕円４７の幅Ｘ２
（＝４・Ｓ２／π・Ｌ）と実際の第１の粒子領域４６の
幅Ｗ２との偏差とする。すなわちＱ２＝｜Ｘ２−Ｗ２
｜。第２の粒子領域において長方形を適合幾何学形状を
した場合の評価値Ｒ２をその長方形４９の幅Ｙ２（＝Ｓ
２／π・Ｌ）と実際の第１の粒子領域４６の幅Ｗ２との
偏差とする。すなわちＲ２＝｜Ｙ２−Ｗ２｜。 〔４−２〕．適合幾何学形状の選択 評価値Ｑ１とＲ１、Ｑ２とＲ２をそれぞれ比べ評価値が
小さい方を適合幾何学形状とする。すなわち、第１の粒
子領域４４に対しＱ１＜Ｒ１なら適合幾何学形状として
楕円、Ｑ１≧Ｒ１なら長方形とし、第２の粒子領域４６
に対しＱ２＜Ｒ２なら楕円、Ｑ２≧Ｒ２なら長方形とす
る。

【００３６】適合幾何学形状が楕円と楕円の場合は楕円
球体の体積情報を算出し、長方形と長方形の場合は直方
体の体積情報を算出し、楕円と長方形の場合には楕円柱
体の体積情報を算出する。

【００３７】

【発明の効果】本発明は、同一粒子に対して異なる２つ
の方向から粒子を撮像し、各画像中の粒子領域に適合す
る幾何学形状を求めているので、その適合幾何学形状に
基いて粒子ごとに体積情報を算出することができる。

【００３８】粒子画像に手寄贈させる幾何学形状とし
て、長方形（正方形を含む）、楕円（円を含む）などの
比較的単純な形状とする場合には演算処理が容易とな
る。体積情報パラメータ以外に第２のパラメータも算出
し分布図を作成することによりより詳細な粒子特性を知
ることができる。

【図面の簡単な説明】

【図１】本発明の粒子画像解析装置の一実施例のブロッ
ク図である。

【図２】本発明の実施例の光学系部分の平面図である。

【図３】本発明の実施例のフローセル部分の概略斜視図
である。

【図４】粒子モデルの図である。

【図５】粒子画像中の第１の粒子領域に幾何学形状を適
合させることを説明するための図である。

【図６】粒子画像中の第２の粒子領域に幾何学形状を適
合させることを説明するための図である。

【図７】ポリエステル粒子の測定結果の例である（スキ
ャッタグラム）。

【図８】ポリエステル粒子の測定結果の例である（ヒス
トグラム）。

【図９】トナー粒子の測定結果の例である（スキャッタ
グラム）。

【図１０】トナー粒子の測定結果の例である（ヒストグ
ラム）。

【符号の説明】

１ 連続発光光源 ２ 第２パルス光源 ３ フローセル ７ 第１光検出器 １１ 第２光検出器 １７ 第２ビデオカメラ １８ 第１パルス光源 ２３ 第１ビデオカメラ ３１ 第１の粒子撮像系 ３２ 第２の粒子撮像系 ３４ 粒子検出系 ３６ 信号処理部 ３７ 撮像制御部 ３８ 画像処理部 ３９ 表示部 ４２ 粒子 ４４ 第１の粒子領域 ４６ 第２の粒子領域 ４５、４７、４８、４９ 適合幾何学形状

Claims (5)

【特許請求の範囲】
1. 【請求項１】 粒子含有試料液をシース液で囲んで流す
   フローセルと、その流れ粒子を順次撮像する粒子撮像系
   と、撮像された粒子画像を処理し解析する解析手段と、
   解析結果等を表示する表示手段と、を備える粒子画像解
   析装置において、 粒子撮像系は、粒子をある面側から見た第１の画像を得
   る第１の粒子撮像系と、上記同一粒子を上記ある面と交
   差する面側から見た第２の画像を得る第２の粒子撮像系
   とを含むものであって、 解析手段は、撮像粒子ごとに、第１、第２の画像中の各
   粒子領域の形状に適合する幾何学形状を決定し、それら
   適合幾何学形状の寸法情報に基いて粒子の体積情報パラ
   メータの値を算出し、さらに、粒子ごとに算出された上
   記体積情報パラメータの値を集計しその分布パターンを
   生成するものであって、 表示手段は、上記解析結果を表示するものであることを
   特徴とする粒子画像解析装置。
2. 【請求項２】 各適合幾何学形状は楕円および長方形の
   いずれかから選ばれたものであり、体積情報パラメータ
   として直方体、楕円球体および円柱体のいずれかの体積
   情報を算出することを特徴とする請求項１記載の粒子画
   像解析装置。
3. 【請求項３】 体積情報パラメータが体積球相当径であ
   ることを特徴とする請求項１記載の粒子画像解析装置。
4. 【請求項４】 解析手段は、さらに、撮像粒子ごとに粒
   子の体積情報パラメータ以外の第２のパラメータの値も
   算出し、表示手段はそれら粒子の体積情報パラメータを
   含む複数のパラメータで規定される複数次元パラメータ
   空間における分布パターンを解析結果として表示するこ
   とを特徴とする請求項１記載の粒子画像解析装置。
5. 【請求項５】 第２のパラメータは粒子の球形度情報パ
   ラメータであることを特徴とする請求項４記載の粒子画
   像解析装置。