JP3130628B2 - 粒子判定装置 - Google Patents

粒子判定装置

Info

Publication number
JP3130628B2
JP3130628B2 JP04042319A JP4231992A JP3130628B2 JP 3130628 B2 JP3130628 B2 JP 3130628B2 JP 04042319 A JP04042319 A JP 04042319A JP 4231992 A JP4231992 A JP 4231992A JP 3130628 B2 JP3130628 B2 JP 3130628B2
Authority
JP
Japan
Prior art keywords
particle
target particle
value
determination
characteristic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP04042319A
Other languages
English (en)
Other versions
JPH05209821A (ja
Inventor
時弘 小坂
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sysmex Corp
Original Assignee
Sysmex Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sysmex Corp filed Critical Sysmex Corp
Priority to JP04042319A priority Critical patent/JP3130628B2/ja
Priority to EP93300046A priority patent/EP0553951A1/en
Priority to US08/011,812 priority patent/US5469375A/en
Publication of JPH05209821A publication Critical patent/JPH05209821A/ja
Application granted granted Critical
Publication of JP3130628B2 publication Critical patent/JP3130628B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/1012—Calibrating particle analysers; References therefor
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00—Scenes; Scene-specific elements
    • G06V20/60—Type of objects
    • G06V20/69—Microscopic objects, e.g. biological cells or cellular parts
    • G06V20/698—Matching; Classification
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1429—Signal processing
    • G01N15/1433—Signal processing using image recognition
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1456—Optical investigation techniques, e.g. flow cytometry without spatial resolution of the texture or inner structure of the particle, e.g. processing of pulse signals
    • G01N15/1459—Optical investigation techniques, e.g. flow cytometry without spatial resolution of the texture or inner structure of the particle, e.g. processing of pulse signals the analysis being performed on a sample stream
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/1012—Calibrating particle analysers; References therefor
    • G01N2015/1014—Constitution of reference particles

Landscapes

  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • Dispersion Chemistry (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Description

【発明の詳細な説明】
【０００１】
【産業上の利用分野】本発明は、検出した粒子が、捜そ
うとしている粒子であるか否か判定する装置に関し、特
に粒子を検出して得た検出信号に基づいて特徴パラメー
タを求め、これに基づいて判定するものに関する。
【０００２】
【従来の技術】従来、様々な種類の粒子が混在している
検体中から、特定の種類の粒子を抽出し、その粒子の二
次元画像を撮像したり、その粒子自体を捕捉収拾したり
して、粒子を分析する粒子分析装置が使用されている。
このような粒子分析装置には、特定の種類の粒子の抽出
用に、粒子判定装置が設けられている。粒子判定装置で
は、粒子を検出領域に流すことによって、個々の粒子か
ら得られた信号を各種信号処理し、各粒子ごとに様々な
特徴パラメータを得ている。
【０００３】これら特徴パラメータとしては、例えば、
散乱光強度、蛍光強度、吸光量、粒子径、粒子面積、Ｎ
／Ｃ比（粒子が血球等のように核と細胞質とからなるも
のの場合、その核と細胞質との比）、円形度等がある。
これら特徴パラメータ値は、各粒子の種類ごとにそれぞ
れ固有の値を持つので、これら特徴パラメータの違いに
より、粒子判定装置は、検出した粒子が目的とする種類
の粒子か否かをリアルタイムに判定している。
【０００４】このような粒子判定装置としては、例えば
図７に示すようなものがあった。この装置では、２つの
特徴パラメータＸ１、Ｘ２を用い、特徴パラメータＸ
１、Ｘ２に対し、上限値Ｘ１ｍａｘ、Ｘ２ｍａｘを上限
値設定レジスタ１０、１２に、下限値Ｘ１ｍｉｎ、Ｘ２
ｍｉｎを、下限値設定レジスタ１４、１６にそれぞれ設
定し、特徴パラメータＸ１が上限値Ｘ１ｍａｘ及び下限
値Ｘ１ｍｉｎによって規定される範囲内にあるとき、そ
れぞれＨレベルの信号を発生するように 比較器１８、
２０を構成し、特徴パラメータＸ２が上限値Ｘ２ｍａｘ
及び下限値Ｘ２ｍｉｎによって規定される範囲内にある
とき、それぞれＨレベルの信号を発生するように、比較
器２２、２４を構成してある。
【０００５】さらに、比較器１８、２０の出力をＡＮＤ
回路２６に供給し、比較器２２、２４の出力をＡＮＤ回
路２８に供給し、これらＡＮＤ回路２６、２８の出力を
ＡＮＤ回路３０に供給してある。従って、特徴パラメー
タＸ１が上限値Ｘ１ｍａｘ及び下限値Ｘ１ｍｉｎによっ
て規定される範囲内にあり、かつ特徴パラメータＸ２が
上限値Ｘ２ｍａｘ及び下限値Ｘ２ｍｉｎによって規定さ
れる範囲内にあるとき、ＡＮＤ回路３０がＨレベルの信
号を発生し、検出した粒子が目的粒子であると判定す
る。
【０００６】
【発明が解決しようとする課題】しかし、上記のような
粒子判定装置では、特徴パラメータＸ１、Ｘ２の内いず
れかが上下限値によって規定された範囲からはずれる
と、目的粒子とは判定されない。図９では、説明を簡略
化するために、特徴パラメータを２種類としたが、実際
にはより多くの特徴パラメータを使用する。従って、わ
ずか１つの特徴パラメータが上下限値から外れただけ
で、目的粒子でないと判定すると、実際には目的粒子で
あるにもかかわらず、目的粒子でないと誤判定する可能
性が高くなる。これを防止しようとすると、上下限値に
よって規定される範囲を広げる必要があるが、そうする
と逆に目的外の粒子を目的粒子と誤判断する可能性が高
くなる。また、このような粒子判定装置では、判定対象
の粒子の目的粒子に対する確からしさを数値で表すこと
ができなかった。
【０００７】
【課題を解決するための手段】本発明は、リアルタイム
で粒子が目的粒子であるか否か判定することができる上
に、判定対象の粒子の目的粒子に対する確からしさを数
値で表すことができ、且つその判定精度を向上させた粒
子判定装置を提供することを目的とするものである。
【０００８】そのため、第１の発明は、個々の粒子をそ
れぞれ検出した信号に基づいて、各粒子ごとにその特徴
パラメータを求め、その違いに基づいて目的粒子である
か否かの判定を行う装置に、データ変換部を備えたもの
で、このデータ変換部は、求めた特徴パラメータを、目
的粒子の特徴パラメータの分布特性またはその近似値に
基づいてデータ変換することにより、上記求めたパラメ
ータに対応する粒子の上記目的粒子に対する確からしさ
を示す値を得るものである。このデータ変換部は、目的
粒子の特徴パラメータの平均値が設定されている設定手
段と、目的粒子の特徴パラメータの分布幅が設定される
設定手段と、求めた特徴パラメータと上記平均値との差
分値を得る差分手段と、上記差分値と上記分布幅とから
上記差分値を正規化した偏差を得る正規化手段と、上記
正規化した偏差を入力して変換する関数手段とを、具備
している。
【０００９】また、第２の発明は、第１の発明における
データ変換部を、各特徴パラメータごとに設け、これら
データ変換部からの出力を演算、比較することにより上
記目的粒子であるか否か判定する判定部を設けたもので
ある。各関数手段は、それらの出力にそれぞれ所定の重
み係数に基づいて重み付けを行うことができる。また、
この重み付けは、各関数手段に複数の重み係数のうち１
つを選択する設定手段を、それぞれ設けるようにするこ
ともできる。さらに、各関数手段は、それぞれ異なる関
数のうち選択されたものに基づいて、変換するようにす
ることもできる。さらに、これらデータ変換部と判定部
とを複数設け、これら判定部の出力を論理演算するよう
にすることもできる。
【００１０】さらに、判定部は、各変換部からの出力を
加算して判定評価値を求める演算手段と、判定基準値が
設定される設定部と、判定評価値と判定基準値とを比較
する比較手段とを、備えたものとすることができる。
【００１１】
【作用】第１の発明によれば、個々の粒子ごとに得られ
た特徴パラメータが、変換部に供給され、ここで逐次デ
ータ変換され、リアルタイムで目的粒子に対する確から
しさを示す値が求められる。また、第２の発明によれ
ば、個々の粒子ごとに得られた複数の特徴パラメータ
が、各パラメータに対応して設けられた各変換部におい
て、目的粒子に対する確からしさを示す値に変換され
る。これら各特徴パラメータごとに求められた目的粒子
に対する確からしさを表す値が、判定部において演算、
比較され、目的粒子であるか否か判定される。
【００１２】
【実施例】図１に、本発明による粒子判定装置の１実施
例の概略ブロック図を示す。この粒子判定装置は、図示
しない検出装置、例えばイメージングフローサイトメー
タによって１個の粒子を検出するごとに発生される特徴
パラメータＸ１、Ｘ２・・・の各チャンネルに対応した
変換部３２−１、３２−２・・・を有し、これら変換部
３２−１、３２−２・・・からの出力信号は、判定部３
４に供給される。
【００１３】変換部３２−１、３２−２・・・は、対応
する特徴パラメータＸ１、Ｘ２・・・を入力し、これを
データ変換して、被検粒子の目的粒子に対する確からし
さを表す値Ｙ１、Ｙ２・・・を求める。値Ｙ１、Ｙ２・
・・は被検粒子がどれだけ目的粒子に近いかを表すもの
である。これらデータ変換は、目的粒子の特徴パラメー
タＸ１、Ｘ２・・・の分布特性（目的粒子であれば、ど
のような値を取るかを示す分布特性、またはその近似特
性）に基づく関数Ｆ１、Ｆ２・・・により、Ｙ１＝Ｆ１
（Ｘ１）、Ｙ２＝Ｆ２（Ｘ２）・・・によって行われ
る。また、これらの変換値Ｙ１、Ｙ２・・・は、更に或
る重み付け係数ａ１、ａ２・・・により重み付けされ、
Ｙａ１（＝ａ１＊Ｙ１）、Ｙａ２（＝ａ２＊Ｙ２）・・
・として出力される。係数ａ１、ａ２・・・の値は、そ
れぞれ０乃至１の定数であり、後述する判定ルールに基
づいて決定される。
【００１４】判定部３４は、上記Ｙａ１、Ｙａ２・・・
を入力し、演算、比較等の処理を行って、その被検粒子
が目的粒子か否かを示す判定結果Ｚ１を出力する。これ
ら変換部３２−１、３２−２・・・、判定部３４での処
理は、リアルタイムでなされる。実際には、各特徴パラ
メータＸ１、Ｘ２・・・の値が揃ってから目的粒子か否
かの判定結果がでるまでに要する時間は、数百ｎｓｅｃ
以内である。
【００１５】図２に、データ変換部３２−１の詳細なブ
ロック図を示す。データ変換部３２−１では、目的粒子
の特徴パラメータＸ１の平均値Ｕ１を平均値レジスタ３
６に、分布幅Ｗ１を分布幅レジスタ３８に設定してあ
る。分布幅Ｗ１は、例えば目的粒子の特徴パラメータＸ
１の標準偏差ＳＤの３倍の値３ＳＤに設定してある。
【００１６】検出された特徴パラメータＸ１と、平均値
Ｕ１とは差分手段４０に入力され、両者の差の絶対値｜
Ｘ１−Ｕ１｜が、差分手段４０より出力される。この差
分手段４０の出力｜Ｘ１−Ｕ１｜と、分布幅レジスタ３
８の分布幅Ｗ１は、正規化手段４２に入力され、正規化
手段４２は、正規化された偏差ΔＸ１を、 ΔＸ１＝ｋ１（｜Ｘ１−Ｕ１｜／Ｗ１） の演算を行って出力する。なお、ｋ１は係数である。
【００１７】正規化偏差ΔＸ１は、関数手段４４に入力
される。関数手段４４は、特徴パラメータＸ１におい
て、その被検粒子が目的粒子にどの程度近いかを表す関
数ｆ１を有し、この関数ｆ１によって目的粒子である確
率を表す値Ｙ１〔＝ｆ１（ΔＸ１）〕を出力する。この
関数ｆ１は、目的粒子がどのような特徴パラメータ値Ｘ
１を取りうるのかという分布情報（ヒストグラム）等を
元に設定されている。
【００１８】このような関数ｆ１の例を図３（ａ）乃至
（ｄ）に示す。それぞれ横軸は正規化偏差ΔＸ１を表
し、縦軸は目的粒子である確率を表しており、同図
（ａ）乃至（ｃ）はいずれも減少関数であり、入力され
る正規化偏差ΔＸ１が大きくなるに従って、被検粒子が
目的粒子である確率が低くなることを示している。なお
同図（ｄ）に示すような増加関数を用いれば、関数手段
４４は目的粒子でない確率を表す値を出力する。
【００１９】目的粒子の特徴パラメータの頻度分布曲線
は、一般には正規分布曲線となるので、関数手段４４で
使用する関数として図３（ｃ）に示す正規分布関数が適
当であるが、これと近似している同図（ｂ）に示す台形
関数でもよい。この台形関数を用いる場合、Ｙ１の値を
正規化関数ΔＸのどの値から減少させるか、またその減
少の勾配をどのように定めるかは実際の評価結果によっ
て定めることができる。
【００２０】同図（ａ）に示す直線関数も使用可能であ
るが、これを使用した場合には、正規化偏差ΔＸ１と目
的粒子である確率との間に線形関係があるとするもので
あり、多変変量解析における重回帰分析法に相当する。
【００２１】図３（ａ）乃至（ｃ）に示した関数は、い
ずれも関数手段４４に設定されており、関数選択手段４
６からの関数選択信号によって選択されたものが関数ｆ
１として使用される。このように関数ｆ１を選択できる
ようにしたのは、汎用性を向上させるためと、これら関
数のうち最も判定精度が良い関数を選択するためであ
る。
【００２２】このように関数手段４４では、正規化偏差
ΔＸ１を基に目的粒子である確率を表す値Ｙ１を求めら
れるが、さらに、関数手段４４は、Ｙ１に或る重み係数
ａ１を乗算した値Ｙａ１（＝ａ１＊Ｙ１）を出力できる
ように構成されている。この重み係数を乗算するのは、
判定部３４での判定ルールに基づく。従って、判定部３
４での判定ルールを変更すると、それに応じて重み係数
ａ１も変更する必要があるので、関数手段４４には複数
の重み係数が設定されており、重み係数設定手段４８か
らの選択信号に応じて、いずれかの重み係数をａ１とし
て使用することができるようになっている。
【００２３】差分手段４０、正規化手段４２、関数手段
４４としては、ルックアップテーブルを用いるのが望ま
しい。ルックアップテーブルは予め計算された数表テー
ブルをメモリに書き込んだもので、リアルタイムの演算
が要求されるときに有効である。
【００２４】上記のようにデータ変換部３２−１は構成
されているが、他のデータ変換部３２−２・・・も同様
な構成である。但し、使用する平均値、分布幅、関数
は、入力される特徴パラメータＸ２・・・に応じたもの
に変更されている。
【００２５】図４に判定部３４のブロック図を示す。各
変換部３２−１、３２−１・・・から各特徴パラメータ
Ｘ１、Ｘ２・・・を変換して得た、目的粒子である確率
を表す値Ｙ１、Ｙ２・・・に、所定の重み係数ａ１、ａ
２・・・を乗算したＹａ１、Ｙａ２・・・は、加算器５
０で加算され、判定評価値Ｔ（＝Ｙａ１＋Ｙａ２・・
・）として出力され、判定基準値設定手段５２に設定さ
れている判定基準値Ｓと比較器５４で比較される。判定
基準値Ｓとしては、例えばＹａ１、Ｙａ２・・・それぞ
れの予測最大値を加算した値に、適当な係数（０乃至
１）を乗算した値を用いている。
【００２６】比較器５４は、判定評価値Ｔが判定基準値
Ｓよりも大きいとき、Ｈレベルの出力Ｖを発生するよう
に構成されている。従って、比較器５４の出力ＶがＨレ
ベルのとき、被検粒子が目的粒子と判定でき、出力Ｖが
Ｌレベルのとき、被検粒子が目的粒子でないと判定でき
る。
【００２７】なお、判定には被検粒子が目的粒子である
と判定したい場合と、被検粒子が目的粒子でないと判定
したい場合とがあり、いずれの判定を行う場合でも、そ
の目的とする判定結果が得られたときには、Ｈレベルま
たはＬレベルいずれか一方の出力が得られるようにした
い場合がある。そのため、この実施例では、比較器５４
の出力Ｖを排他的論理和回路５６の一方の入力に供給
し、他方の入力にはスイッチ５８のオン・オフの切換に
よってＨレベルまたはＬレベルの信号が供給されるよう
に構成してある。
【００２８】即ち、スイッチ５８をオンとして、排他的
論理和回路５６の他方の入力にＬレベルの信号を供給す
ると、比較器５４の出力Ｖが反転されるので、被検粒子
が目的粒子であると判定された場合には、排他的論理和
回路５６からＬレベルの出力Ｚ１が得られる。また、ス
イッチ５８をオフとして、排他的論理和回路５６の他方
の入力にＨレベルの信号を供給しておくと、比較器５４
の出力Ｖがそのまま排他的論理和回路５６から出力され
るので、被検粒子が目的粒子でないと判定されたとき
に、排他的論理和回路５６からＬレベルの出力Ｚ１が得
られる。
【００２９】このように、この粒子判定装置によって被
検粒子が目的粒子であるか否かの判定が行われるが、そ
の判定に基礎となっているルールとしては、例えば次の
ようなものがある。これらの例は、特徴パラメータとし
てＸ１乃至Ｘ３の３つを使用した場合のものである。
【００３０】（例１） ルール：条件１ 特徴パラメータＸ１が目的粒子の特徴
パラメータＸ１の平均値に近く、 条件２ 特徴パラメータＸ２が目的粒子の特徴パラメー
タＸ２の平均値に近く、 条件３ 特徴パラメータＸ３が目的粒子の特徴パラメー
タＸ３の平均値に近い場合、 その粒子は目的粒子Ａである。被検粒子がＡであるなら
ば、判定出力（排他的論理和回路５６の出力）Ｚ１をＬ
レベルにする。但し、条件１乃至３の重要度は同じとす
る。 判定条件：条件１、２、３に示すように本例では特徴パ
ラメータＸ１、Ｘ２、Ｘ３を用いて被検粒子が目的の粒
子Ａであるか否かの判定を行うものであるので、予め、
粒子Ａの特徴パラメータＸ１、Ｘ２、Ｘ３に関する設定
値が設定される。具体的には変換部３２−１、３２−
２、３２−３における平均値、分布幅、関数、重み係数
が設定される。なお、条件１ないし３の重要度を同じに
するため、重み係数ａ１乃至ａ３には同一の値を用い
る。また、被検粒子がＡである場合にＺ１をＬレベルに
するため、スイッチ５６をオンとしておく。このような
状態において、粒子を測定して粒子ごとに特徴パラメー
タＸ１、Ｘ２、Ｘ３を得、変換部３２−１、３２−２、
３２−３により値Ｙａ１、Ｙａ２、Ｙａ３を求め、さら
に、 判定評価値Ｔ＝Ｙａ１＋Ｙａ２＋Ｙａ３ を求める。そして、この判定評価値Ｔについて 判定評価値Ｔ≧判定基準値Ｓ の判定を行うことにより被検粒子が目的粒子Ａの場合に
判定出力Ｚ１をＬレベルとすることができる。なお、被
検粒子がＡでない場合に、Ｚ１をＬレベルにする場合、
スイッチ５６をオフとする。
【００３１】（例２） ルール：条件１ 特徴パラメータＸ１が目的粒子の特徴
パラメータＸ１の平均値に近く、 条件２ 特徴パラメータＸ２が目的粒子の特徴パラメー
タＸ２の平均値に近く、 条件３ 特徴パラメータＸ３が目的粒子の特徴パラメー
タＸ３の平均値に近い場合、 その粒子は目的粒子Ａである。被検粒子がＡであるなら
ば、判定出力（排他的論理和回路５６の出力）Ｚ１をＬ
レベルにする。但し、条件１乃至３の重要度は、この順
に高、中、低とする。 判定条件：（例１）と同様に、変換部３２−１、３２−
２、３２−３の各手段の設定をする。ただし、条件１乃
至３の重要度を高、中、低とするため、重み係数ａ１乃
至ａ３は ａ１＞ａ２＞ａ３ とする。そして、被検粒子がＡである場合にＺ１をＬレ
ベルにするため、スイッチ５６をオンとする。これによ
って 判定評価値Ｔ（＝Ｙａ１＋Ｙａ２＋Ｙａ３）≧判定基準
値Ｓ で判定出力Ｚ１をＬレベルとすることができる。なお、
被検粒子がＡでない場合に、Ｚ１をＬレベルにするに
は、スイッチ５６をオフとする。
【００３２】このように変換部３２−１、３２−２・・
・の平均値、分布幅、関数によって求められた目的粒子
である確率を表す値Ｙ１、Ｙ２・・・に、それぞれ重み
係数ａ１、ａ２・・・を乗算して、重み付けすることに
よって、粒子判定における条件の重要度を加味して、精
度の良い判定を行うことができる。なお、重み係数、例
えばａ１を０とした場合、この重み係数ａ１に対応する
特徴パラメータＸ１を考慮しない、即ち無視することに
なる。
【００３３】上記の例１、例２は、いずれも判定ルール
が１つの場合であるが、或る被検粒子が目的粒子である
か判定するのに複数のルールを使用することがある。例
えば、被検粒子が目的粒子であると判定できる場合とし
て、（１）被検粒子の特徴パラメータＸ１、Ｘ２の値
が、それぞれ目的粒子の特徴パラメータＸ１、Ｘ２の平
均値に近く、または（２）被検粒子の特徴パラメータＸ
３、Ｘ４の値が、それぞれ目的粒子の特徴パラメータＸ
３、Ｘ４の平均値に近い場合がある。なお、各特徴パラ
メータの重要度はいずれも同じとする。このような場合
には、図５に示すように、（１）に対して１つの判定ブ
ロック６０−１を、（２）に対して１つの判定ブロック
６０−２・・・を、それぞれ設けなければならない。
【００３４】無論、判定ブロック６０−１、６０−２・
・・は、いずれも図１に示したように、変換部と、判定
部とからなる。これら変換部は、図２に示したものと同
様に構成され、判定部は、図４に示したものと同様に構
成されている。
【００３５】各判定ブロック６０−１、６０−２・・・
からの判定出力（各判定ブロック内の判定部からの出
力）Ｚ１、Ｚ２・・・は、判定条件を満たすときＬレベ
ルになるようにワイヤードＯＲライン６２によってワイ
ヤードＯＲ接続されており、このワイヤードＯＲ出力Ｗ
Ｚは、排他的論理和回路６４の一方の入力に供給されて
いる。この排他的論理和回路６４の他方の入力には、ス
イッチ６６のオンオフに応じてＨレベルまたはＬレベル
の出力が供給されている。
【００３６】従って、スイッチ６６がオンのときには、
ワイヤードＯＲ出力ＷＺが反転され、スイッチ６６がオ
フのときには、ワイヤード出力ＷＺがそのまま出力され
る。スイッチ６６は、通常にはオンとされ、判定ブロッ
ク６０−１、６０−２のいずれかがＬレベルの判定出力
を発生したとき（図４の判定部３４のスイッチ５８に相
当する、判定ブロック６０−１、６０−２の判定部のス
イッチ５８−１、５８−２をオンとしている場合）、即
ち２つの判定ルールのうちいずれか一方を満たしたと
き、Ｈレベルの最終判定結果Ｚを出力する。
【００３７】もし、上記の目的粒子以外の粒子を判定す
る場合には、スイッチ６６をオフとすればよい。これに
よって、いずれの判定ブロック６０−１、６０−２もＨ
レベルの判定出力が発生したとき、排他的論理和回路６
４はＨレベルの最終判定結果を出力する。
【００３８】このように判定ルールごとに判定ブロック
６０−１、６０、２を設け、これら判定ブロックの出力
をワイヤードＯＲによって論理演算して、排他的論理和
回路６４に入力し、かつスイッチ６６のオン、オフに応
じて最終判定結果Ｚを正転／逆転させているので、リア
ルタイムで様々な判定ルールに対応することが可能とな
る。また、上記のように論理演算をワイヤードＯＲによ
って行っているので、判定ブロックの追加、削除が容易
に行える。
【００３９】図６は、粒子画像の撮像機能を備えたフロ
ーサイトメータであるイメージングフローサイトメータ
に、本発明による粒子判定装置を使用した場合のブロッ
ク図であり、この場合には、目的粒子のみを選択して、
それの静止画像を撮像することが可能となる。
【００４０】イメージングフローサイトメータ本体７０
は、試料の偏平なシースフローを形成するためのフロー
セル機構、粒子個々の特徴を検出するための光源、光検
出器、粒子を撮像するためのストロボ光源、ビデオカメ
ラ等を備えている。
【００４１】フローセル中を流れる個々の粒子に対して
得られた検出信号は、信号処理装置７２に供給され、信
号処理装置７２は各粒子ごとにそれぞれ特徴パラメータ
Ｘ１、Ｘ２・・・をリアルタイムで求める。これら特徴
パラメータＸ１、Ｘ２・・・は、判定ブロック６００−
１、６００−２、６００−３・・・に供給される。これ
ら判定ブロック６００−１、６００−２、６００−３・
・・は、図７に示した判定ブロック６０−１、６０−２
・・・にそれぞれ対応するものである。これら判定ブロ
ック６００−１、６００−２、６００−３・・・は、特
徴パラメータが供給されてから数百ｎｓｅｃ以内に判定
出力を出力する。
【００４２】これら判定出力は、符号６２１で示すよう
にワイヤードＯＲされて、信号処理装置７２に戻され
る。信号処理装置７２内には、図５の排他的論理和回路
６４、スイッチ６６に対応する回路が設けられており、
これらによって最終判定される。この最終判定の結果、
目的粒子であると判定されると、被検粒子を撮像するた
めのストロボトリガ信号をイメージングフローサイトメ
ータ本体７０に供給し、この本体７０では、ストロボ光
源が発光し、被検粒子の撮像がビデオカメラによって行
われる。
【００４３】信号処理装置７２で求めた各特徴パラメー
タは、データ解析装置７４にも供給され、このデータ解
析装置７４によって各特徴パラメータごとにヒストグラ
ムや多次元スキャッタグラム等の解析が行われる。この
解析結果を用いて、粒子判定ルールを決めたり、各特徴
パラメータの平均値、分布幅、重み、判定基準値等の判
定条件を決めて、それらの値を判定ブロックに設定する
こともできる。また、関数手段４４への非線形関数テー
ブルのロードや判定結果の反転等のコントロールを行う
こともできる。
【００４４】なお、本発明の粒子判定装置をセルソータ
に適用した場合には、目的粒子のみを採集することがで
きる。
【００４５】上記の実施例では、複数の特徴パラメータ
を用いて被検粒子が目的粒子であるか否かを判定したの
で、各特徴パラメータごとに変換部を設けたが、１つの
特徴パラメータのみによって被検粒子が目的粒子である
か否かを判定する場合には、変換部は１台だけでよく、
その場合、重み付けは不要となる。
【００４６】
【発明の効果】請求項１、２記載の発明によれば、被検
粒子の目的粒子に対する確からしさを数値をリアルタイ
ムによって得ることができ、大容量のメモリを用いず
に、比較的小容量のメモリだけでデータ変換を行うこと
ができ、効率的で経済的である。また、各手段で処理す
るデータも単純化できる。しかも、関数手段は、現実の
分布特性に基づいた非線形関数であるので、粒子の確か
らしさを示す値をより高精度に求めることができる。
【００４７】さらに、請求項３記載の発明によれば、各
変換部からの出力に重み付けを行っているので、各特徴
パラメータの重要度を加味した判定が行え、高精度の判
定を行える。請求項４，５記載の発明では、重み係数や
関数を選択あるいは設定することができるので、汎用性
を増すことができる。請求項６記載の発明によれば、粒
子ごとに得た複数の特徴パラメータからリアルタイムで
被検粒子が目的粒子であるか否かを判定できる。請求項
７記載の発明によれば、リアルタイムで様々な判定ルー
ルに対応することができる。
【図面の簡単な説明】
【図１】本発明による粒子判定装置の１実施例のブロッ
ク図である。
【図２】同実施例のデータ変換部の詳細なブロック図で
ある。
【図３】同実施例のデータ変換部で使用する関数を示す
図である。
【図４】同実施例の判定部の詳細なブロック図である。
【図５】同実施例の粒子判定装置を複数の判定ルールに
基づくものに変更した場合のブロック図である。
【図６】同次子異例の粒子判定装置をイメージングフロ
ーサイトメータと共に使用する場合のブロック図であ
る。
【図７】従来の粒子判定装置のブロック図である。
【符号の説明】
３２−１、３２−２ データ変換部 ３４ 判定部 ３６ 平均値レジスタ ３８ 分布幅レジスタ ４０ 差分手段 ４２ 正規化手段 ４４ 関数手段 ４６ 関数設定手段 ４８ 重み係数設定手段

Claims (7)

    (57)【特許請求の範囲】
  1. 【請求項１】 個々の粒子をそれぞれ検出した信号に基
    づいて上記各粒子ごとにその特徴パラメータを求め、そ
    の違いに基づいて目的粒子であるか否かの判定を行う装
    置に、上記求めた特徴パラメータを、目的粒子の特徴パ
    ラメータの分布特性またはこれの近似値に基づいてデー
    タ変換することにより、上記求めたパラメータに対応す
    る粒子の上記目的粒子に対する確からしさを示す値を得
    るデータ変換部を、備え、上記データ変換部が、上記目
    的粒子の特徴パラメータの平均値が設定されている設定
    手段と、上記目的粒子の特徴パラメータの分布幅が設定
    される設定手段と、上記求めた特徴パラメータと上記平
    均値との差分値を得る差分手段と、上記差分値と上記分
    布幅とから上記差分値を正規化した偏差を得る正規化手
    段と、上記正規化した偏差を入力して変換する非線形関
    数手段とを、具備する粒子判定装置。
  2. 【請求項２】 個々の粒子をそれぞれ検出した信号に基
    づいて上記各粒子ごとにそれぞれ異なる複数の特徴パラ
    メータを求め、これら各特徴パラメータの違いに基づい
    て目的粒子であるか否かの判定を行う装置において、 上記求めた各特徴パラメータを、上記目的粒子の各特徴
    パラメータ分布特性またはこれの近似値に基づいてデー
    タ変換することにより、上記求めた各特徴パラメータご
    とに、これら求めた特徴パラメータに対応する粒子の上
    記目的粒子に対する確からしさを示す値をそれぞれ求め
    るデータ変換部と、 これらデータ変換部からの出力を演算、比較することに
    より上記目的粒子であるか否か判定する判定部と、 を具備し、上記各データ変換部が、上記目的粒子の特徴
    パラメータの平均値が設定される設定手段と、上記求め
    た特徴パラメータ値と上記平均値の差分値を求める差分
    手段と、上記目的粒子の特徴パラメータの分布幅が設定
    される設定手段と、上記差分手段からの差分値と上記特
    徴パラメータの分布幅とから上記差分値を正規化した偏
    差を求める正規化手段と、上記正規化した偏差を入力し
    て非線形変換する関数手段とを、それぞれ具備すること
    を特徴とする粒子判定装置。
  3. 【請求項３】 請求項２記載の粒子判定装置において、
    上記各関数手段が、それらの出力にそれぞれ所定の重み
    係数に基づいて重み付けを行うことを特徴とする粒子判
    定装置。
  4. 【請求項４】 請求項３記載の粒子判定装置において、
    上記各関数手段は複数の重み係数を有し、これら重み係
    数のうちいずれかを選択する設定手段を設けたことを特
    徴とする粒子判定装置。
  5. 【請求項５】 請求項２記載の粒子判定装置において、
    上記各関数手段は、それぞれ異なる関数のうち選択され
    たものに基づいて変換することを特徴とする粒子判定装
    置。
  6. 【請求項６】 個々の粒子をそれぞれ検出した信号に基
    づいて上記各粒子ごとにそれぞれ異なる複数の特徴パラ
    メータを求め、これら各特徴パラメータの違いに基づい
    て目的粒子であるか否かの判定を行う装置において、 上記求めた各特徴パラメータを、上記目的粒子の各特徴
    パラメータ分布特性またはこれの近似値に基づいてデー
    タ変換することにより、上記求めた各特徴パラメータご
    とに、これら求めた特徴パラメータに対応する粒子の上
    記目的粒子に対する確からしさを示す値をそれぞれ求め
    るデータ変換部と、 これらデータ変換部からの出力を演算、比較することに
    より上記目的粒子であるか否か判定する判定部と、 を具備し、上記判定部が、上記各変換部からの出力を加
    算して判定評価値を求める演算手段と、判定基準値が設
    定される設定部と、上記判定評価値と上記判定基準値と
    を比較する比較手段とを、備えていることを特徴とする
    粒子判定装置。
  7. 【請求項７】 個々の粒子をそれぞれ検出した信号に基
    づいて上記各粒子ごとにそれぞれ異なる複数の特徴パラ
    メータを求め、これら各特徴パラメータの違いに基づい
    て目的粒子であるか否かの判定を行う装置において、 上記求めた各特徴パラメータを、上記目的粒子の各特徴
    パラメータ分布特性またはこれの近似値に基づいてデー
    タ変換することにより、上記求めた各特徴パラメータご
    とに、これら求めた特徴パラメータに対応する粒子の上
    記目的粒子に対する確からしさを示す値をそれぞれ求め
    るデータ変換部と、 これらデータ変換部からの出力を演算、比較することに
    より上記目的粒子であるか否か判定する判定部と、 を具備する粒子判定装置を、複数のルールに対応して複
    数台設け、これら粒子判定装置の出力を論理演算するこ
    とを特徴とする粒子判定装置。
JP04042319A 1992-01-30 1992-01-30 粒子判定装置 Expired - Fee Related JP3130628B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP04042319A JP3130628B2 (ja) 1992-01-30 1992-01-30 粒子判定装置
EP93300046A EP0553951A1 (en) 1992-01-30 1993-01-06 Particle judging device
US08/011,812 US5469375A (en) 1992-01-30 1993-02-01 Device for identifying the type of particle detected by a particle detecting device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP04042319A JP3130628B2 (ja) 1992-01-30 1992-01-30 粒子判定装置

Publications (2)

Publication Number Publication Date
JPH05209821A JPH05209821A (ja) 1993-08-20
JP3130628B2 true JP3130628B2 (ja) 2001-01-31

Family

ID=12632700

Family Applications (1)

Application Number Title Priority Date Filing Date
JP04042319A Expired - Fee Related JP3130628B2 (ja) 1992-01-30 1992-01-30 粒子判定装置

Country Status (3)

Country Link
US (1) US5469375A (ja)
EP (1) EP0553951A1 (ja)
JP (1) JP3130628B2 (ja)

Families Citing this family (42)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3111706B2 (ja) * 1992-02-18 2000-11-27 株式会社日立製作所 粒子分析装置及び粒子分析方法
JP2826449B2 (ja) * 1993-09-17 1998-11-18 株式会社日立製作所 フロー式粒子画像解析方法およびフロー式粒子画像解析装置
EP0720012B1 (en) * 1994-12-26 2004-09-08 Sysmex Corporation Flow cytometer
JP3504030B2 (ja) * 1995-07-04 2004-03-08 シスメックス株式会社 粒子判定基準の決定方法およびその装置並びにその判定基準を用いた粒子分析装置
JP3504029B2 (ja) * 1995-07-04 2004-03-08 シスメックス株式会社 粒子分析装置
GB9606423D0 (en) * 1996-03-27 1996-06-05 Univ Hertfordshire An instrument for the real-time classification of particle shape within clouds and aerosols
WO1998034094A1 (en) 1997-01-31 1998-08-06 The Horticulture & Food Research Institute Of New Zealand Ltd. Optical apparatus
US6149867A (en) 1997-12-31 2000-11-21 Xy, Inc. Sheath fluids and collection systems for sex-specific cytometer sorting of sperm
US6887711B1 (en) * 1999-01-22 2005-05-03 Cytospec Method of characterization of biological entities
JP3817091B2 (ja) * 1999-07-02 2006-08-30 テルモ株式会社 細胞数測定装置
US7208265B1 (en) 1999-11-24 2007-04-24 Xy, Inc. Method of cryopreserving selected sperm cells
GB2368390B (en) * 2000-03-16 2004-01-21 Univ Loughborough Analysing data
AU2002237689B2 (en) 2000-11-29 2008-01-10 Xy, Llc. System to separate frozen-thawed spermatozoa into X-chromosome bearing and Y-chromosome bearing populations
US7713687B2 (en) 2000-11-29 2010-05-11 Xy, Inc. System to separate frozen-thawed spermatozoa into x-chromosome bearing and y-chromosome bearing populations
US8486618B2 (en) 2002-08-01 2013-07-16 Xy, Llc Heterogeneous inseminate system
EP2275533B9 (en) 2002-08-01 2016-10-19 Xy, Llc Method of assessing sperm cells
WO2004017041A2 (en) 2002-08-15 2004-02-26 Xy, Inc. High resolution flow cytometer
US7169548B2 (en) 2002-09-13 2007-01-30 Xy, Inc. Sperm cell processing and preservation systems
ES2561816T3 (es) 2003-03-28 2016-03-01 Inguran, Llc Aparatos, métodos y procesos para clasificar partículas y para proporcionar esperma animal clasificado por sexo
CA2566749C (en) 2003-05-15 2017-02-21 Xy, Inc. Efficient haploid cell sorting for flow cytometer systems
MXPA06011344A (es) 2004-03-29 2006-12-15 Monsanto Technology Llc Suspensiones de espermatozoides para separar en poblaciones enriquecidas con cromosomas x o y.
BRPI0513685A (pt) 2004-07-22 2008-05-13 Monsanto Technology Llc processo para enriquecimento de uma população de células de esperma
JP4649231B2 (ja) * 2005-02-28 2011-03-09 株式会社カネカ フローサイトメータ、細胞の解析方法、細胞解析プログラム、蛍光検出器の感度設定方法および陽性率判定法における基準ゲート設定方法
US7996188B2 (en) * 2005-08-22 2011-08-09 Accuri Cytometers, Inc. User interface for a flow cytometer system
US8303894B2 (en) 2005-10-13 2012-11-06 Accuri Cytometers, Inc. Detection and fluidic system of a flow cytometer
US8017402B2 (en) 2006-03-08 2011-09-13 Accuri Cytometers, Inc. Fluidic system for a flow cytometer
US7299135B2 (en) 2005-11-10 2007-11-20 Idexx Laboratories, Inc. Methods for identifying discrete populations (e.g., clusters) of data within a flow cytometer multi-dimensional data set
US7780916B2 (en) * 2006-03-08 2010-08-24 Accuri Cytometers, Inc. Flow cytometer system with unclogging feature
US8283177B2 (en) * 2006-03-08 2012-10-09 Accuri Cytometers, Inc. Fluidic system with washing capabilities for a flow cytometer
US8077310B2 (en) * 2006-08-30 2011-12-13 Accuri Cytometers, Inc. System and method of capturing multiple source excitations from a single location on a flow channel
US8715573B2 (en) * 2006-10-13 2014-05-06 Accuri Cytometers, Inc. Fluidic system for a flow cytometer with temporal processing
US8445286B2 (en) 2006-11-07 2013-05-21 Accuri Cytometers, Inc. Flow cell for a flow cytometer system
US7739060B2 (en) * 2006-12-22 2010-06-15 Accuri Cytometers, Inc. Detection system and user interface for a flow cytometer system
JP4389991B2 (ja) * 2007-10-26 2009-12-24 ソニー株式会社 微小粒子の光学的測定方法及び光学的測定装置
US8432541B2 (en) * 2007-12-17 2013-04-30 Accuri Cytometers, Inc. Optical system for a flow cytometer with an interrogation zone
US8004674B2 (en) * 2009-06-02 2011-08-23 Accuri Cytometers, Inc. Data collection system and method for a flow cytometer
US8507279B2 (en) * 2009-06-02 2013-08-13 Accuri Cytometers, Inc. System and method of verification of a prepared sample for a flow cytometer
US8779387B2 (en) * 2010-02-23 2014-07-15 Accuri Cytometers, Inc. Method and system for detecting fluorochromes in a flow cytometer
US9551600B2 (en) 2010-06-14 2017-01-24 Accuri Cytometers, Inc. System and method for creating a flow cytometer network
ES2897531T3 (es) 2010-10-25 2022-03-01 Accuri Cytometers Inc Sistemas e interfaz de usuario para la recopilación de un conjunto de datos en un citómetro de flujo
CN110462372B (zh) 2017-05-25 2022-06-14 佛罗乔有限责任公司 大型多参数数据集的可视化、比较分析和自动差异检测
EP3867625A4 (en) * 2018-10-17 2022-08-24 Becton, Dickinson and Company ADAPTIVE SORTING FOR PARTICLE ANALYZERS

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3826364A (en) * 1972-05-22 1974-07-30 Univ Leland Stanford Junior Particle sorting method and apparatus
US3870872A (en) * 1974-02-01 1975-03-11 Abbott Lab Probability analog function computer
US4008388A (en) * 1974-05-16 1977-02-15 Universal Monitor Corporation Mass spectrometric system for rapid, automatic and specific identification and quantitation of compounds
US4284412A (en) * 1979-07-13 1981-08-18 Ortho Diagnostics, Inc. Method and apparatus for automated identification and enumeration of specified blood cell subclasses
US4741035A (en) * 1983-06-01 1988-04-26 Head Stereo Gmbh Wide band, low noise artificial head for transmission of aural phenomena
US4661913A (en) * 1984-09-11 1987-04-28 Becton, Dickinson And Company Apparatus and method for the detection and classification of articles using flow cytometry techniques
CA1233250A (en) * 1985-09-23 1988-02-23 Gerald B. Price Method and apparatus for multiple category threshold and window events counting in flow cytometric studies
US5041733A (en) * 1987-03-20 1991-08-20 Agency Of Industrial Science & Technology Method and apparatus for identifying chromosomes or cells
US5124932A (en) * 1988-03-10 1992-06-23 Indiana University Foundation Method for analyzing asymmetric clusters in spectral analysis
JP2667867B2 (ja) * 1988-03-30 1997-10-27 東亜医用電子株式会社 粒子解析装置
JP2674705B2 (ja) * 1988-06-10 1997-11-12 東亜医用電子株式会社 一次元分布分画方法
US5040112A (en) * 1988-12-07 1991-08-13 Serono-Baker Diagnostics, Inc. Method of separating the three major types of blood cells from a white blood cell histogram

Also Published As

Publication number Publication date
EP0553951A1 (en) 1993-08-04
US5469375A (en) 1995-11-21
JPH05209821A (ja) 1993-08-20

Similar Documents

Publication Publication Date Title
JP3130628B2 (ja) 粒子判定装置
US7668388B2 (en) System and method for single image focus assessment
US5911002A (en) Pattern recognition system
JP6391083B2 (ja) 自動欠陥分類における分類器の準備及び保守
CA1223355A (en) Method and apparatus for automatic quantitative measurement of textures by image analysis
JP3581149B2 (ja) 境界ピクセル・パラメータの規則正しいシーケンスを用いた物体を識別するための方法及び装置
US6941009B2 (en) Method for evaluating pattern defects on a water surface
US20040234120A1 (en) Defect classification method
JP4961161B2 (ja) 検査装置
US20100272351A1 (en) Information processing apparatus and method for detecting object in image data
JP2006098155A (ja) 検査方法及びその装置
CN109063574B (zh) 一种基于深度神经网络检测的包络框的预测方法、***及设备
CN115311484B (zh) 一种焊缝数据的优化聚类分割方法
KR102629852B1 (ko) 다중-모드 검사를 위한 통계적 학습-기반 모드 선택
US20140185918A1 (en) Charged particle beam apparatus that performs image classification assistance
JP7393313B2 (ja) 欠陥分類装置、欠陥分類方法及びプログラム
CN117635604B (zh) 一种建筑工程钢结构焊接质量检测方法
US7136809B2 (en) Method for performing an empirical test for the presence of bi-modal data
JP2001188906A (ja) 画像自動分類方法及び画像自動分類装置
JP2006292615A (ja) 外観検査装置、外観検査方法、コンピュータを外観検査装置として機能させるためのプログラムおよび記録媒体
JP3002325B2 (ja) 表面検査装置
CN114758470B (zh) 一种基于消防工程的火灾预警方法及***
US20120044262A1 (en) Surface observation apparatus and surface observation method
JP7390851B2 (ja) 欠陥分類装置、欠陥分類プログラム
JP2001174417A (ja) 欠陥検出装置及び方法

Legal Events

Date Code Title Description
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20001107

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees