JP2009522586A - 撮像装置を用いて層流中の粒子を計数する方法とシステム - Google Patents
撮像装置を用いて層流中の粒子を計数する方法とシステム Download PDFInfo
- Publication number
- JP2009522586A JP2009522586A JP2008550360A JP2008550360A JP2009522586A JP 2009522586 A JP2009522586 A JP 2009522586A JP 2008550360 A JP2008550360 A JP 2008550360A JP 2008550360 A JP2008550360 A JP 2008550360A JP 2009522586 A JP2009522586 A JP 2009522586A
- Authority
- JP
- Japan
- Prior art keywords
- counting
- particles
- algorithm
- particle density
- flow
- 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.)
- Pending
Links
- 239000002245 particle Substances 0.000 title claims abstract description 243
- 238000003384 imaging method Methods 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims description 55
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 90
- 239000007788 liquid Substances 0.000 claims abstract description 37
- 238000005286 illumination Methods 0.000 claims abstract description 20
- 238000005070 sampling Methods 0.000 claims abstract description 9
- 239000012530 fluid Substances 0.000 claims description 17
- 230000004044 response Effects 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 4
- 238000012544 monitoring process Methods 0.000 claims description 4
- 230000002572 peristaltic effect Effects 0.000 claims description 4
- 239000003086 colorant Substances 0.000 claims description 3
- 210000004027 cell Anatomy 0.000 description 73
- 238000001514 detection method Methods 0.000 description 20
- 238000012545 processing Methods 0.000 description 18
- 241000894006 Bacteria Species 0.000 description 14
- 238000000684 flow cytometry Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 9
- 238000005259 measurement Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 239000007850 fluorescent dye Substances 0.000 description 5
- 230000001580 bacterial effect Effects 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 244000005700 microbiome Species 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 241000588724 Escherichia coli Species 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000000427 antigen Substances 0.000 description 2
- 102000036639 antigens Human genes 0.000 description 2
- 108091007433 antigens Proteins 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 201000011510 cancer Diseases 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000004163 cytometry Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000005375 photometry Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 241000224489 Amoeba Species 0.000 description 1
- 210000001744 T-lymphocyte Anatomy 0.000 description 1
- 210000001766 X chromosome Anatomy 0.000 description 1
- 210000002593 Y chromosome Anatomy 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000003975 animal breeding Methods 0.000 description 1
- 239000007864 aqueous solution Substances 0.000 description 1
- XKRFYHLGVUSROY-UHFFFAOYSA-N argon Substances [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 description 1
- 229910052786 argon Inorganic materials 0.000 description 1
- 210000000601 blood cell Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000001444 catalytic combustion detection Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 239000000975 dye Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 244000078673 foodborn pathogen Species 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 210000002443 helper t lymphocyte Anatomy 0.000 description 1
- 238000010166 immunofluorescence Methods 0.000 description 1
- 239000006194 liquid suspension Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 244000000010 microbial pathogen Species 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000000582 semen Anatomy 0.000 description 1
- 238000011896 sensitive detection Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007790 solid phase Substances 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 210000004881 tumor cell Anatomy 0.000 description 1
- 238000004148 unit process Methods 0.000 description 1
Images
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/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/14—Optical investigation techniques, e.g. flow cytometry
- G01N2015/1486—Counting the particles
Landscapes
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Physics & Mathematics (AREA)
- Dispersion Chemistry (AREA)
- Analytical Chemistry (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Apparatus Associated With Microorganisms And Enzymes (AREA)
Abstract
ビデオカメラ等を用いて流液(214)を撮像することによって、流動(層)液体中の蛍光粒子(210)の濃度を測定する発明が開示される。照明ビームは、対象粒子を照明するのに使用される。撮像オプティクス(408、412)は、流液ストリーム中の複数蛍光試料粒子の画像を形成するために焦点面を観測するように配置される。カメラは撮像オプティクスによって形成された画像を記録し、計数アルゴリズム(812)は粒子を計数する。粒子密度(828)の初期推定値、例えば、流量(824a)、露出時間(824)およびサンプリング間隔(824a)に応じて、システムの動作パラメータが調整される(824)。さらに、推定粒子密度によって、計数アルゴリズムが選択される。
Description
本発明は、半透明または透明な流動液体中の特定対象粒子を計数する装置ならびに方法に関する。特に、本発明は画像化された粒子を計数することに関し、流液中において推定された粒子密度に応じて、装置パラメータおよび信号処理が調整される。
体積の大きい流体中において、密度の低い選択された対象粒子、特に細胞を撮像しおよび分類することに関する出願は多くある。このような出願には、1.)バイオテロと生物戦争防御、2.)食物と水質保全、3.)癌細胞の臨床的検出、および4.)環境監視を含む。これまで開発された細胞撮像計数システムにおいては、1.)高コスト、2.)不十分な感度、3.)遅い、4.)大きいサイズ、5.)不十分なスペクトルおよび/もしくは空間的分解能、並びに/または、6.)労働集約的な準備ステップの問題を通常抱えている。
従来のフローサイトメトリを使用して、細胞を直接的に計数しても良い。フローサイトメトリは一般的に使用されている技術であり、細胞の化学的または物理的な性質を測定することができる。細胞は、流体中、大抵は空気または水中において浮遊されたまま、測定装置によって一列となって流動する。免疫蛍光フローサイトメトリにおいては、各細胞に蛍光抗体を付着させることによって、以下のように、細胞を特定することができる。
・蛍光性分子または蛍光光色素を用いて対象となる細胞に対して特定の抗体をラベルする。
・対象となる細胞と共に、ラベルされた抗体を溶液中で混ぜる。抗体が(抗原と称される)細胞の特定のサイトに付着する。
・細胞を(複数の)レーザ付近の溶液流中に一列で通過させ、レーザは蛍光色素を照明して、細胞は異なる波長で蛍光する。
・光電子増倍管またはフォトダイオードを用いて、マークされた細胞が検出装置の前を通る毎に蛍光放射散乱が検出される。
・マークされた細胞の数を計数することができる。対象となる細胞に特有の抗体を選択することができる。
・蛍光性分子または蛍光光色素を用いて対象となる細胞に対して特定の抗体をラベルする。
・対象となる細胞と共に、ラベルされた抗体を溶液中で混ぜる。抗体が(抗原と称される)細胞の特定のサイトに付着する。
・細胞を(複数の)レーザ付近の溶液流中に一列で通過させ、レーザは蛍光色素を照明して、細胞は異なる波長で蛍光する。
・光電子増倍管またはフォトダイオードを用いて、マークされた細胞が検出装置の前を通る毎に蛍光放射散乱が検出される。
・マークされた細胞の数を計数することができる。対象となる細胞に特有の抗体を選択することができる。
流動細胞計測は、現在、種々のアプリケーションに使用されており、HIV治療をモニタするためのヘルパーTリンパ球数を測定するアプリケーションや、癌治療の判定において腫瘍細胞DNA含有量を測定するアプリケーションや、動物育種のXとY染色体ベアリング精液を分離するアプリケーションに使用されている。
図1(先行技術)は、典型的なフローサイトメトリシステム(非特許文献1)を示している。サイトメトリを実行に移すことは、2つの同心円筒の流液を用いることを含む。内部流または中心流液は、サンプリングされるべき細胞を含んでいる。外側のストリームまたはシース流液を用いる目的は、中心流の直径を小さくすることである。中心およびシース流体が流液のテーパ状の領域に達するので、中心流の断面面積は小さくなる。小口径の中心流(約20ミクロン)によって、小さい直径のレーザ光線により照明された流液中の細胞が光学的に精度良く測定される。すなわち、細胞のすべてが、ビームのほとんど同一の部分を通り抜けて、均一に照明されるだろう。なぜ、細胞は小口径透明チューブを通過しないか?一般に、小さい直径の孔は機能しない。なぜならば、頻繁に目詰まりするからである。市販のフローサイトメータの全てが現在、シース/中心流配置を採用している。
フローサイトメータ内の蛍光ラベルのレーザ励起蛍光は、食物が媒介する病原体等の特定の微生物を、迅速、確実でかつ比較的明確に検知する唯一の強力な方法である。数点の研究論文がローサイトメータの方法およびアプリケーションについて説明している(例えば、A.L.GivanによるFlow Cytometry:First Principles、1992年、参照文献)。
歴史的には、フローサイトメータは非常に大きくて、高価な実験室ベースの器具である。フローサイトメータは、大電力を消費して、複雑なエレクトロニクスを採用している。フローサイトメータは、携帯型装置の分野内では通常考慮されていない。従来のフローサイトメータのサイズ(最小でもデスクトップサイズ)、所要電力量、および目詰まりのしやすさの故に、水の生物汚染のフィールドモニタリング等の多くのアプリケーションに利用することはできない。
特許文献1においては、ボイラー中の発光粒子を計数する撮像手段が開示されているものの、
・一般のフローサイトメータまたは特に噴水流フローサイトメータには適用しておらず、
・粒子を取り込むための流路を含んでおらず、
・対象粒子中の蛍光を励起するための照明装置を含んでおらず、
・ポンピング機構およびポンプ率を制御する手段を含んでおらず、
・粒子が1以上の回数で計数されないことを保証する方法を含んでおらず、
・システム操作パラメータを急いで最適化するのに使用されるべき測定計数をなす手段に備えておらず、
・空間的な錯誤に対する計数を修正することによって、正確に粒子を計数する方法を含んでおらず、
・計数評価基準として、粒子サイズを含んでおらず、
・粒子密度に基づいて複数の計数アルゴリズムを使用する手段を含んでいない。
・一般のフローサイトメータまたは特に噴水流フローサイトメータには適用しておらず、
・粒子を取り込むための流路を含んでおらず、
・対象粒子中の蛍光を励起するための照明装置を含んでおらず、
・ポンピング機構およびポンプ率を制御する手段を含んでおらず、
・粒子が1以上の回数で計数されないことを保証する方法を含んでおらず、
・システム操作パラメータを急いで最適化するのに使用されるべき測定計数をなす手段に備えておらず、
・空間的な錯誤に対する計数を修正することによって、正確に粒子を計数する方法を含んでおらず、
・計数評価基準として、粒子サイズを含んでおらず、
・粒子密度に基づいて複数の計数アルゴリズムを使用する手段を含んでいない。
特許文献4においては、透明媒体中の蛍光ラベルされたDNAを高スループットで分析する発明が開示されている。この発明は、流液中で撮像装置に近づく細胞を検出するデバイスである。流液は、非常に細長いビームを伴うレーザによって側面から焦点面に照明された透明チューブ内にある。
この発明は代替技術として上に示された欠点を有していないものの、流液媒体が透明でないアプリケーションには適していない。また、それは撮像技術ではなく、むしろ単一ポイントの測光的な検出および解析に適した技術である。
特許文献5および特許文献6には、細胞を計数する発明が開示されており、その計数は、比較的平坦な透明チャンバに細胞を流し、それらの蛍光を励起させ、カメラを用いて単一画像でチャンバを撮像することによってなされる。カメラの画像はデジタル化され、チャンバ内の細胞数が計数される。
特許文献3においては、特許文献5と同様に細胞を撮像するための発明が開示されており、試料がフローチャンバに流れて、撮像される。
特許文献2においては、所望の不確実性が得られるまで、一定の不確実性に定量化されるべき細胞を含む液体サンプルは、一度に測定される1個のサブサンプルである発明が開示されている。
図2(従来技術)は、本願の発明者による特許文献7に詳細に開示されている(そして本願明細書に参照として組み入れられた)噴水フローサイトメータと称されるシステムの概略を示している。手短に説明すれば、以下のとおりである。
−蛍光ラベル細胞(210)のサンプルは、CCD/CMOSカメラおよび前部オプティクス(208)に向かってチューブ(206)の上方に流れる。
−細胞は、レーザ(228)により焦点面内で照明される。
−細胞が、CCD/CMOSカメラ焦点面(234)を通り抜けるとき、透明窓を介するCCD/CMOSカメラ(218)およびレンズアセンブリ(212)、並びに蛍光発光の波長を分離するフィルタ(214)によって撮像される。
−細胞が浮遊している流体は窓のそばを通過し、排液管(230)から出る。
−蛍光ラベル細胞(210)のサンプルは、CCD/CMOSカメラおよび前部オプティクス(208)に向かってチューブ(206)の上方に流れる。
−細胞は、レーザ(228)により焦点面内で照明される。
−細胞が、CCD/CMOSカメラ焦点面(234)を通り抜けるとき、透明窓を介するCCD/CMOSカメラ(218)およびレンズアセンブリ(212)、並びに蛍光発光の波長を分離するフィルタ(214)によって撮像される。
−細胞が浮遊している流体は窓のそばを通過し、排液管(230)から出る。
通常、サンプルは垂直に上向きに流れる。(1実施例では、サンプルは水平な排液管(230)から流れ、ガラスカバースリップ(220)下の流路(222)から流れるだろう)。体積(volume)226は、我々の発明においては、空気または透明な液体で満たされる。これにより、溶媒を介在せず、介在する流体の光学的深度を0付近まで低減させずに、照明された体積(224)を見ることができる。照明された体積(224)は、薄い必要ではなく(すなわち、リボンとして形成されず)、またレーザを通して照明される必要はない。発光ダイオード等のインコヒーレント光を介した照明は、使用可能であるだけでなく、いくつかの点においては好ましい。照明された体積における照明の不均一性を増大させる際の干渉の効果を最小にするからである。
サンプル(206)の列は、エンドオンで観察されおよび照明される(これを実行する1つの方法としては、観察路と照明路の両方をフロー軸と一致させ、反対方向においては、エピ蛍光顕微鏡と同様の二色性ミラーを用いる方法がある)。撮像された体積は、カメラ前部オプティクスの被写界深度によって制御されて、通常、結果のデジタル画像を縮小する逆効果を防ぐために小さく維持されるであろう。理想的な状況は、流管に流れる粒子が少数ピクセルで撮像されるところにある。
利便性良く撮像される流管の断面は、サンプル流体の不透明度に依存しておらず、非常に大きく(3mm×3mm平方)することができるので、サンプルの流出量を高くすることができる。我々は、直径8mmの噴水流口を使用して、100ml/分より大きい流量で、LEDによって照明したアメーバを検知した。
チューブ壁付近の粒子の流れが遅いので、壁付近の粒子がより低速で移動し、その結果、複数の画像中に現れる。デバイスを較正するときには、この効果を考慮しなければならない。
選択的に微生物にタグ付けをする蛍光抗体の使用に組み合わされる場合、この発明により、2つの応用例を例示すれば、食物および水中の病原微生物の検出が可能となる。サイズ判定に基づく選別を可能にして、高解像度の画像により、蛍光タグ付けされた粒子のサイズを測定することができる。リアルタイムでの検出、計数およびサイズ判定を実行することができる。
しかしながら、粒子密度の大きいダイナミックレンジに渡って対象粒子密度を正確に測定するために、この技術を実用化するためには、画像処理手段および計数手段において、以下に示す点を改良しなければならない。
・ビデオ画像から、ゆっくり変化するものの比較的強度の高いバックグラウンドを低減させて、バックグラウンドと対象粒子とを識別できるようにする。
・ノイズを減らす。
・候補粒子のサイズおよび測光的強度を含み、感度ならびに客観的基準に基づいて、粒子を検出する。
・ソース錯誤の効果を減らす、すなわち、2つの粒子が1つの粒子として計数される1つのフレーム内に十分近接して共に出現する統計的確率を訂正する。
・1つ以上の画像に現れていたとしても、その粒子を1回だけ計数する。
・検出装置の焦点面を通り抜けるとき、対象粒子の各々の測光的強度およびそのサイズを正確に測定する。
・システムの操作パラメータを変えて、対象粒子密度に迅速に適応し、粒子密度を測定するシステムのダイナミックレンジを拡張する。これらの操作パラメータは、流速、露光時間、照明の明度、並びに画像の処理および計数に用いられるアルゴリズムを含んでいる。
・ビデオ画像から、ゆっくり変化するものの比較的強度の高いバックグラウンドを低減させて、バックグラウンドと対象粒子とを識別できるようにする。
・ノイズを減らす。
・候補粒子のサイズおよび測光的強度を含み、感度ならびに客観的基準に基づいて、粒子を検出する。
・ソース錯誤の効果を減らす、すなわち、2つの粒子が1つの粒子として計数される1つのフレーム内に十分近接して共に出現する統計的確率を訂正する。
・1つ以上の画像に現れていたとしても、その粒子を1回だけ計数する。
・検出装置の焦点面を通り抜けるとき、対象粒子の各々の測光的強度およびそのサイズを正確に測定する。
・システムの操作パラメータを変えて、対象粒子密度に迅速に適応し、粒子密度を測定するシステムのダイナミックレンジを拡張する。これらの操作パラメータは、流速、露光時間、照明の明度、並びに画像の処理および計数に用いられるアルゴリズムを含んでいる。
半透明な流動液体中の特定粒子を高スループット画像ベースで計数する改善された装置および方法の必要性が、当技術分野にある。流液中の算出された粒子密度に応じて、装置パラメータおよび信号処理が調整される必要性もある。
本発明の目的は、半透明または透明な流動液体中における特定粒子の高スループットイメージングベースの計数のための改良された装置と方法を提供することにある。流液中の推定された粒子密度に応じて、装置パラメータおよび信号処理が調整される。
この技術の1つの好適な実施例は、特許文献7の噴水流発明を利用するフローサイトメトリに関連している。しかし、この実施例に制限されず、流液の断面を抽出する焦点面を有する撮像装置によってモニタされるいかなる流液においても、対象粒子の計数に使用される。これは、従来のフローサイトメトリおよび特許文献4による発明を含んでいる。
本明細書において説明する発明は、新規な統合されたシステムであり、撮像装置と、照明器と、流液路と、バックグラウンドおよびノイズを低減させるイメージプロセッシング手段と、粒子を計数する計数手段と、計数値を修正して正確に粒子密度を測定する分析手段と、システムパラメータを調整して短い時間間隔で正確な測定を最適化する制御手段とを含む。
本発明のいくつかの実施例は、半透明または透明な流動液体中のサンプルの高スループット、高感度検出および識別をする装置ならびに方法を提供する。これは、比較的大きい断面図の軸流を与えることによって達成され、軸流中においては、流液中に浮遊する細胞または対象粒子は、撮像装置の焦点面を通して流れるときに、観測される。広いダイナミックレンジに亘る対象粒子密度の測定は、対象粒子密度の迅速な測定に基づいて、リアルタイムでシステムの動作パラメータを調整することによって、達成される。
本発明には、固相サイトメトリ(Lemarchand et al., Aquat. Microb. Ecol., 25: 301−309, 200)が可能であるという利点があり、はるかに低いコスト(<1万ドル)で個々の微生物を検出できるという利点も含まれる。本発明は、免疫蛍光色素を含む蛍光色素を使用した分類学的識別に基づいている。検出ステップは、以下のステップを含む。すなわち、
・対象細胞種に特有の抗体が、蛍光分子または蛍光色素によってラベル付けされる。
・ラベル付けされた抗体は、対象細胞種と共に溶液中で混ぜられる。抗体が細胞(抗原)上の特定位置に付着する。
・細胞の液体浮遊液はレーザまたはLEDに向かう流路内に送り込まれ、レーザまたはLEDは蛍光色素を照明し、これにより、細胞が長波長で蛍光を発する。対象細胞に対して、極めて特有の抗体を選択することができる。
・低コストCCD(電荷撮像素子)またはCMOS二次元検出装置が流路を撮像し、マークされたセルが検出装置の焦点面を通り抜ける間において、そのセルがレーザーダイオードによって照明される度に蛍光放射のフラッシュが撮像される。
・マークされた細胞の数がコンピュータを介して計数され、その結果は、測定される対象粒子の密度用にシステムのパラメータを最適化するのに用いられる。
・対象細胞種に特有の抗体が、蛍光分子または蛍光色素によってラベル付けされる。
・ラベル付けされた抗体は、対象細胞種と共に溶液中で混ぜられる。抗体が細胞(抗原)上の特定位置に付着する。
・細胞の液体浮遊液はレーザまたはLEDに向かう流路内に送り込まれ、レーザまたはLEDは蛍光色素を照明し、これにより、細胞が長波長で蛍光を発する。対象細胞に対して、極めて特有の抗体を選択することができる。
・低コストCCD(電荷撮像素子)またはCMOS二次元検出装置が流路を撮像し、マークされたセルが検出装置の焦点面を通り抜ける間において、そのセルがレーザーダイオードによって照明される度に蛍光放射のフラッシュが撮像される。
・マークされた細胞の数がコンピュータを介して計数され、その結果は、測定される対象粒子の密度用にシステムのパラメータを最適化するのに用いられる。
本発明は、流路を通り抜けるサンプル中の複数の蛍光粒子を撮像する装置および方法を含む。1つの好適な実施例では、流路は、フローストリーム中のサンプルに対して流れ方向を定義し、流液を分岐する観測面を有する。照明ビームは、流液中の画像焦点面において、対象粒子からの蛍光を励起するのに用いられる。撮像オプティクスは、フローストリーム中の複数の蛍光サンプル粒子の画像を形成するために焦点面を観測するように配置されている。カメラは、撮像オプティクスによって形成された画像を記録する。コンピュータは、カメラによって記録されたデジタル画像上において、粒子検出、計数、システム制御およびデータ解析を実行する。
撮像要素は、カラーフィルタと、オプティクスと、CCDまたはCMOSカメラ等の撮像要素とを含んでも良い。サンプルの流量を維持するポンプ要素は、シリンジ・ポンプ、ぜん動ポンプ、または他のコンピュータの制御可能なポンプから成っても良い。レーザまたはLEDにより照明される。
本発明は流液中の粒子を計数する方法を提供し、流液は、照明された流液の断面を撮像する焦点面を有する撮像装置を含む撮像システムによってモニタされる。その方法は、
(a)撮像装置を用いて1つ以上のフレームを撮像するステップと、
(b)既定の計数アルゴリズムを使用して、撮像されたフレームに基づいて粒子密度の予備的推定を行なうステップと、
(c)撮像されたフレームの密度が収まる2以上の密度カテゴリを決定するステップと、
(d)決定された密度カテゴリに応じて、
−計数アルゴリズム、
−撮像システムのパラメータ、
−流速
のうち少なくとも1つの要素を調整するステップと、
(e)ステップ(c)の調整をした後に、流液中の粒子を計数するステップとを含む。
(a)撮像装置を用いて1つ以上のフレームを撮像するステップと、
(b)既定の計数アルゴリズムを使用して、撮像されたフレームに基づいて粒子密度の予備的推定を行なうステップと、
(c)撮像されたフレームの密度が収まる2以上の密度カテゴリを決定するステップと、
(d)決定された密度カテゴリに応じて、
−計数アルゴリズム、
−撮像システムのパラメータ、
−流速
のうち少なくとも1つの要素を調整するステップと、
(e)ステップ(c)の調整をした後に、流液中の粒子を計数するステップとを含む。
上記の調整ステップは、露光時間または照明強度を含む撮像システムのパラメータを調整しても良い。また、上記の調整するステップは、異なる粒子密度カテゴリで動作するように設計された複数の計数アルゴリズムから選択することによって、計数アルゴリズムを調整しても良い。上述の方法は、粒子の強度分布を測定し、測定された強度分布に応じて照明強度、露光時間およびポンプスピードを調整する。
計数アルゴリズムが修正されるケースでは、計数アルゴリズムは少なくとも低位粒子密度アルゴリズムおよび高粒子密度アルゴリズムを含む。低位粒子密度アルゴリズムは、2つの連続した画像レームの差分を取るステップと、所定の閾値強度を超えたピクセルのグループをフラグ付け(flagging)ステップと、フラグ付けされたグループを粒子として計数するステップと、空間的な錯誤に対して計数を訂正するステップとを実行する。一般に、粒子計数ステップは、粒子を選別するためにフラグ付けされたグループのサイズを利用している。空間的な錯誤に対して計数を訂正するステップは、一般に、計数された粒子間の距離に基づいて、所定の評価基準に応じて計数を減らすことになる。
好適な実施例では、高粒子密度アルゴリズムは、所定の閾値強度を超えたピクセルのグループをフラグ付けするステップと、フラグ付けされたグループを粒子として計数するステップと、計数に基づく数式に応じて計数に空間的な錯誤補正因子を適用するステップとを含む。一般に、フラグ付けするステップに先立って、その処理はフレームに高域通過フィルタをかける。
好適な実施例では、低位粒子密度アルゴリズムは、低粒子密度アルゴリズムおよび中間粒子密度アルゴリズムにさらに分類される。この場合、中間粒子密度アルゴリズムは、計数された粒子が先のフレーム内で計数されているかを判定するステップと、それに応じて計数を減らすステップとを含む。また、中間粒子密度アルゴリズムは、連続したフレーム内に検出された運動粒子を1つの粒子として計数するために、連続した画像内のx、y座標を追跡するステップを含んでも良い。複数のフレーム内で検出された粒子の明度が最大測定明度となるようにとられる。
調整するステップは、以下の様に撮像システムの露光時間およびサンプリング時間を調整しても良い。
−低粒子密度アルゴリズムにおいては、調整ステップはポンプスピードを高い値に設定しおよび露光時間を長い時間に設定する。
−中間粒子密度アルゴリズムにおいては、調整ステップはポンプスピードを中位に設定しおよび露光時間を中位に設定する。
−高粒子密度アルゴリズムにおいては、調整ステップはポンプスピードを低い値に設定しおよび露光時間を短い時間に設定する。
−低粒子密度アルゴリズムにおいては、調整ステップはポンプスピードを高い値に設定しおよび露光時間を長い時間に設定する。
−中間粒子密度アルゴリズムにおいては、調整ステップはポンプスピードを中位に設定しおよび露光時間を中位に設定する。
−高粒子密度アルゴリズムにおいては、調整ステップはポンプスピードを低い値に設定しおよび露光時間を短い時間に設定する。
場合によっては、上述の流速が、流液の所定の体積内の粒子の全てが実質的に1つの画像内で検出されるように調節される。いくつかの実施例においては、フレーム撮像ステップは、2色のフレームを撮像し、および2色の各々の強度に対して対象粒子の各々を測定するステップをさらに含む。この方法は、撮像されたフレームの明度を既定の状態で測定するステップと、測定された明度に応じて照明または撮像システム露出時間を調整するステップとをさらに含む。都合がよければ、計数ステップは、計数アルゴリズムを決定する画像からのデータを用いて、撮像が完了した後に、一組の撮像されたフレームにおいて実行される。
本発明に係る流液中の粒子を計数する装置は、以下の(1)−(3)を含む。
(1)1つ以上の画像フレームを記録しおよび照明された流液の断面を撮像する焦点面を有する撮像装置を含み、流液をモニタする撮像システム。
(2)以下のモジュールを実装するプロセッサであって、そのモジュールは、以下のとおりある。
−プロセッサ内に記憶された既定の計数アルゴリズムを使用して、撮像システムからの撮像されたフレームに基づいて、粒子密度を推定する予備的推定モジュールと、
−推定された粒子密度に基づいて、撮像された前記フレームの粒子密度が収まる2つ以上の所定の密度カテゴリを判定する密度カテゴリ化モジュール
−判定された密度カテゴリに応じて、
−プロセッサ内に記憶された計数アルゴリズム、
−撮像システムのパラメータ、および/または、
−流速
の中から選択された1つ以上の要素を調整する制御信号を生成する制御システムモジュール
−要素が調整された後に、計数アルゴリズムに応じて流液中の粒子を計数する計数モジュール
(3)前記制御信号に応じて、選択された要素を調整する調整ユニット。
(1)1つ以上の画像フレームを記録しおよび照明された流液の断面を撮像する焦点面を有する撮像装置を含み、流液をモニタする撮像システム。
(2)以下のモジュールを実装するプロセッサであって、そのモジュールは、以下のとおりある。
−プロセッサ内に記憶された既定の計数アルゴリズムを使用して、撮像システムからの撮像されたフレームに基づいて、粒子密度を推定する予備的推定モジュールと、
−推定された粒子密度に基づいて、撮像された前記フレームの粒子密度が収まる2つ以上の所定の密度カテゴリを判定する密度カテゴリ化モジュール
−判定された密度カテゴリに応じて、
−プロセッサ内に記憶された計数アルゴリズム、
−撮像システムのパラメータ、および/または、
−流速
の中から選択された1つ以上の要素を調整する制御信号を生成する制御システムモジュール
−要素が調整された後に、計数アルゴリズムに応じて流液中の粒子を計数する計数モジュール
(3)前記制御信号に応じて、選択された要素を調整する調整ユニット。
上述のプロセッサは制御信号を生成して、撮像システムのパラメータを調整し、その場合、調整ユニットが制御信号に応じて撮像システムの露出時間または照明強度を変えても良い。または、プロセッサは、制御信号を生成して計数アルゴリズムを調整し、その場合、調整ユニットは制御信号に応じて計数モジュールにより用いられた計数アルゴリズムを変更する。
上述の撮像システムはレンズを含み、撮像装置はCCD画像要素、CMOS画像要素またはビデオカメラのうちいずれか1つを含んでも良い。流液はフローサイトメータ内で生じても良い。焦点面は前記流液の方向に対して実質的に垂直であり、実質的に平行であり、またはその他の構成であっても良い。流速が、流液の所定体積内の粒子の全てが実質的に1つの画像内で検出されるように調節されても良い。制御可能な流量で流速を決定するシリンジ・ポンプまたは蠕動ポンプのうちいずれか一方をさらに含んでも良い。
本発明は、流動液体中の粒子を計数する装置および方法を含む。流動液体中の粒子の濃度に応じて、装置パラメータおよび信号処理パラメータの両方が調整される。
図3Aおよび図3Bは、本発明に係る噴水流サイトメータから得られた一連の画像である。図3Cは、2つの画像の差分を示してある。一連の画像の差分を取ることが図9と図10とともに以下に説明するように利用される。したがって、図3Cは、(低粒子密度に対する)アルゴリズム808Aまたは(中間の粒子密度に対する)アルゴリズム808Bから生じ得る図を示している。
図4、図5および図6は、流液中の粒子濃度に応じて、粒子検出がどのように変わるかを示している。低濃度では、各々の粒子が計数され得る。高濃度では、検出粒子の低い数から実際の粒子の数を推定しなければならない。
図4は、実際の粒子密度と、条件1と称する低濃度、高流量速度条件下で測定された粒子密度との関係を示すプロットである。予測された細胞計数と測定された細胞計数との間には、強い線形相関があり、傾きは約1/2である。
条件1は流量が十分に高い条件であるので、カメラが1つの画像を表示するのにかかる時間、すなわち440ミリ秒よりもはるかに短時間内に、粒子が照明された焦点面内に留まる。例えば、細胞が<<440ミリ秒で焦点面を通過する場合、細胞が検出される見込みは約露光時間440ミリ秒である。焦点面内のセル残留時間は、検出には十分長いと想定される。これらの条件下では、個々の細胞は複数のフレームには検出されないだろう。この条件では、焦点面を通して流れる細胞の実際の数値に検出(p)確率を掛けることによって計数されている。いくつかのカメラ/条件に対して、サイクルタイムは露光時間に等しく、検出確率は1(unity)である(その場合は、連続した読み出しを可能にするCMOSイメージャを用いた場合である。)。確率が1であるとき、条件1および条件2は同一となる。
図5は、実際の粒子密度と、中濃度、低流量速度条件(条件2)下で測定された粒子密度との関係を示すプロット図である。
条件2下では、予測された細胞計数と測定された細胞計数とには、ほぼ線形的な関係があり、傾きは約1である。ソース錯誤が高濃度で効いてくるまで、条件2では、照明された焦点面における細胞残留時間がカメラのサイクルタイムより大であるかまたは等しく、画像面を通過するあらゆる細胞が、少なくとも一度撮像される。この条件では、複数フレーム(フレーム間において空間的に同時発生:coincidence)内で撮像された細胞に対して、訂正(correction)をしなければならない。さらに、個々の画像内の細胞濃度は十分低いと想定されるので、1つのフレーム内の空間的な同時発生、すなわち、ソース錯誤は顕著ではない。条件2の下では、(引き続くフレーム内のほぼ同一の座標で見られる1つのフレーム内にあるいかなる細胞をも、1つの検出として、計数して)フレーム間の空間的な同時発生を訂正することによって、焦点面を通して流れる粒子の実際の数値を測定することができる。
図6は、実際の粒子密度と、高濃度、低流量速度条件(条件3)において測定された粒子密度との関係を示すプロット図である。
条件3では、1つのフレーム内の細胞濃度が十分高いので、ソース錯誤が顕著となる。フレーム間の空間的な同時発生が非常に大きいので、それらは無視されている。これは、条件1下のように焦点面を通して流れる細胞の実際の数値をもはや測定していないことを意味している。しかし、実際の細胞数と計数された細胞数との間における経験的な関係が判定されている。この判定は、所定の対象粒子密度の静的な較正液体(calibration liquid)を連続的にスナップショットすることによってなされる。言い換えれば、高濃度な静的サンプルのスナップショットによって、焦点面内の何百もの細胞を検知することができる。異なった濃度のサンプルをスナップショットすることによって、計数と実際の濃度との間において経験的な関係を得ることができる。ソース錯誤はこの関係に影響するであろう。
条件3の下では、ソース錯誤を示すほど高くない濃度においてのみ、測定された計数と実際の計数との間には、線形的な相関がある。ソース錯誤が顕著になるにつれ、測定された計数は飽和し始める。
図7は、本発明に係る粒子計数装置の1つの実施例の構成要素を示すブロックダイアグラムである。図8は、粒子計数プロセスで実行される基本ステップを示すフロー図である。2つの図の間には密接な関係があるので、これらの2つの図を一緒に議論する。
図7は、本発明に係る粒子計数装置の実施例を示すブロックダイアグラムである。この実施例は、フローセル配列、2つの照明要素および2つの撮像要素を含んでいるが、他の多くの構成も可能である。図7に示す好適な実施例は、(本明細書に参考として共に援用されている)特許文献7および特許文献8の噴水流の発明を利用するフローサイトメトリに関連している。しかし、この実施例に制限されることはなく、いなかる流液中にある対象粒子の計数に使用することができる。流液の断面は流液の断面をサンプリングする焦点面を有する撮像装置によってモニタされる。これは、従来のフローサイトメトリおよび特許文献4の発明を含んでいても良い。
図7の実施例では、ポンプ460によって、流体はフローセル300内を流れる。流体は、光学要素452、454、442および444を介して、レーザ448、450によって照明される。流液から放射された光434は、光学要素444、442、および412を介してCCD408、408bに進み、ビームスプリッタ456によって分岐される。
CCD408およびCCD408bは装置の画像処理部に画像データを与え、画像処理部は、図8に示されているようにデータを処理する。処理部によって決定された粒子密度に応じて、処理部は、順番に、露光時間、サンプリング時間およびポンプスピード等の装置パラメータを調整する。
エシュリキア属大腸菌(Escherichia coli)は、以下のパラメータを使用して、本発明に係る噴水流システムによって検出される。大腸菌は、水溶液中で浮遊している間、SYBRゴールド染色液(Invitrogen, Eugene, Oregon)を使用して検出される。2mmの開口部は、22ミリワット電力のアルゴンイオンレーザによってフロー開口部において照明された。画像は、Omega XF3105フィルタを介したElectrim 1000L CCDカメラによって撮られた。
これらのパラメータにおける低密度とは、8.33x10−3ml/sの流量と100msの露光時間において、1mlあたり1万未満の細胞が存在する密度であった。これらのパラメータにおける中間密度とは、1.17x10−3ml/sの流量と100msの露光時間において、1mlあたり1万〜10万の細胞が存在する密度であった。これらのパラメータにおける高密度とは、5.83x10−4ml/sの流量、25〜100 msの露光時間において、1mlあたり10万より多い細胞が存在する密度であった。
当業者であれば、特定の構成および条件に応じて、特定の粒子密度、流量および露光時間は変動することを理解するであろう。本発明は、おおよその粒子密度を決定し、システムのシステムパラメータを変えて、その密度における粒子を最善に計数する。
図8は、粒子計数の全体的なプロセスを示すフローチャートである。このプロセスは、ステップ802で開始する。ステップ804において、CCD408(図7を参照)から得られた初期画像404が与えられる。初期信号処理ステップ806は、初期画像中の粒子密度を推定する。サブプロセス808は、推定された粒子密度に応じて変わるステップ810、812および814を含んでいる。図9−図11を参照されたい。
一般的に説明すると、ステップ810は、流速および露光時間等の装置パラメータを選択して、これらのパラメータをフロー−撮像機器に与える(図7を参照されたい)。ステップ812は、粒子密度環境に対して、正しい画像プロセッシング/粒子計数アルゴリズムを選択する(図9−11を参照されたい)。ステップ814は、選択されたアルゴリズムに応じて、CCD408からの画像404に対して全ての画像処理および粒子計数を遂行する。その結果826がステップ816に与えられ、ステップ816においては、ステップ806で推定された粒子密度が、ステップ814で検出されたより正確な粒子密度に対して十分に近接しているかどうかが判定される。近接している場合(818)、粒子データはステップ822において利用するために出力される。近接していない場合、フィードバックループ820において、新たな粒子密度計数がステップ810に与えられ、より正確な結果が求められる。
本発明の好適な実施例では、自動化されたバクテリア認識および計数コンピュータプログラムが、例えば、CCD/CMOSカメラ画像のバクテリアを計数する。
画像内で検出された粒子密度に応じて、プログラムによって適用されたアルゴリズムが選択される。低粒子密度が検出される場合、図9に示されたアルゴリズム808Aが使用される。中間の粒子密度が検出される場合、図10に示されたアルゴリズム808Bが使用される。高粒子密度が検出される場合、アルゴリズム808Cが使用される。低密度および中間の密度アルゴリズムによって、2つの連続した画像の差分を取りおよび負ピクセルをゼロに設定することによって、バックグラウンドが除去される。
噴水流焦点面の照明の不均一性が原因で、CCD画像の光度測定のバックグラウンドは不均一なため、2つの連続する画像の差分を取りおよび差分における負ピクセルをゼロに設定することによって、バックグラウンドは、プログラム(図3)内で排除される。対象粒子(例えば、バクテリア)およびノイズだけが、差分画像に残っている。対象粒子が1つの画像からその次の画像まで横方向に動くので、高濃度のバクテリア以外においては、1つのフレーム内の粒子は連続したフレーム内の粒子からは、通常差し引かれない。
例として、アルゴリズム808Aおよびアルゴリズム808Bの画像差分ステップは、表1中に示された3つのパラメータについてユーザに問い合わせを行う。(CCD検出器がコラムからコラムへの固定したパターンオフセットを示すように、)プログラムは、画像内の各コラムに対する中央値および標準偏差を測定する。次に、プログラムは、ノイズを超過したそれらの信号強度に基づいて、バクテリアと思われるピクセルを識別する。次に、プログラムは、バクテリアの候補を識別しおよび計数する。その識別および計数は、1つの連続したグループ内で検出された候補バクテリアピクセルの数と、そのグループの中心点の特定半径内の全積分光度測定強度とに基づいている。そのようなプログラムは多くのCCD画像でテストされており、目によって識別されたバクテリアの数はプログラムによって計数されたバクテリアの数とほぼ等しい。最終的に、プログラムは複数の連続する画像内で検出されたバクテリアの座標を比較して、複数の連続する画像内で検出された明度の高いスポットのどれが、(所定の距離内にそれらの座標があるので)同一のバクテリアであるのかを判定する。
画像の差分を取ることは、バクテリアの濃度が高い画像に対しては上手く機能しない。第1に、バクテリア密度が高い場合、複数の連続する画像内における同様のx、y座標で定義される粒子が同一のバクテリアを示しているのかどうか、または同一の座標を一致して共有する2つのバクテリアを示しているのかを、正確に判定することはできない。この効果は、ソース錯誤と呼ばれている。第2に、1つのフレームをその次のフレームから差し引くことによって、1つのバクテリアを連続する画像内における別のバクテリアから差し引いてしまう可能性がある。バクテリアの密度が高い場合にあっては、連続する画像の差分を取ることは実行されない。代わりに、アルゴリズム808Cは、座標のフレーム内の空間的な同時存在(coincidence)を無視して、これらに対して修正はしない。アルゴリズム808Cは、測定された噴水流計数と較正サンプルの計数との間における経験的な関係を構築することに依存している。アルゴリズム808Cは、バックグラウンドを差し引くのに連続する画像の差分を使用していない。代わりに、メジアンフィルタリングが、画像内の点放射源(バクテリア)から生ずる広範のバックグラウンドを取り除くのに使用される。言い換えれば、メジアン(median)フィルタ(1つの実施例では10×10ピクセル平方)を適用して、データから低空間周波数が取り除かれる。
このように、粒子を含む画像内の粒子密度に応じて、本発明は異なったアルゴリズムを粒子計数に使用している。理想的には、流速および画像露光時間等の装置パラメータも、粒子計数を最適化するように調整される。図8は、本発明の全体的な粒子計数アルゴリズムを例示しており、図9−図11は、異なった粒子密度条件における粒子計手法の変形例を例示している。
流動内の対象粒子の密度測定タスクの手順が、以下に与えられている。以下に与えられた画像処理、粒子検出および計数ステップと同時に、撮像(データ取得)手順が実行されるの。したがって、画像収集中に、すなわちリアルタイムで計数を表示することができる。好適な実施例は、特許文献7に記載されているような噴水流サイトメトリの使用を想定している。
ステップ806において、予備的な粒子計数が実行される。例として、このプロセスは以下のステップを含んでも良い。
1.低速で(または定常流で)送り込み、流液の「スナップショット」(1つの画像)を取る。
2.画像または画像中の列または行に対して平均値(または、中央値)および標準偏差を測定する。後者は、列から列へまたは行から行へと変動する画像内の固定されたパターンオフセットを取り扱うときに役立つ。
3.画像にメジアンフィルタをかける。
4.隣接ピクセルの最小数で構成された画像内の明るいスポットまたは縞の数を計数する。隣接ピクセルの各々は、平均値(または、中間値)を超えた標準偏差の選択数よりも大きい強度を有する。
5.明るいスポットを含む明るいピクセルのすべての強度を加算することによって、明るいスポット(または縞)の各々の概ねの強度を計算する。
6.画像中の対象粒子の密度(撮像された体積によって割られた計数値であり、体積=領域の深さ×照明されおよび撮像された焦点面の面積)を計算する。
7.対象粒子の密度および強度分布を使用して、画像処理、計数、分析並びに最適な流量、露光時間および照明強度の設定のために使用されるフローおよびアルゴリズムの適切な条件番号(number)を見つける。
8.粒子密度に応じて、操作パラメータを変える。
9.適切なアルゴリズム8A、B、またはCを使用して、最終的な測定をする。
図9は、低粒子密度条件下における図8のアルゴリズムで利用された特定のアルゴリズムを示すフローチャートである。図10は、中間的な粒子密度条件下における図8のアルゴリズムで利用された特定のアルゴリズムを示すフローチャートである。これらアルゴリズムの両方が、流液中の個別粒子の大部分またはすべてを検出するのに画像差分を利用している。
図9に図示したアルゴリズム808Aは、以下のステップを含んでいる。
ステップ902は、制御信号824aを介してポンプスピードを高い(High)に設定する。
ステップ904は、制御信号824bを介して露光時間およびサンプリング時間を長い(long)設定する。
ステップ906は、データストリーム404を介してCCD408からの画像を取得する。対象粒子密度の定量を所望の精度にするために、一連の画像が撮られる。
ステップ908は、前の画像を差し引くことによって各画像の差分をとり、バックグラウンドが差し引かれた一連の画像を生成する。負ピクセルがゼロに設定され、差し引かれた画像から対象粒子の影響が取り除かれる。
ステップ910は、「ホット」ピクセルにフラグ付け(flag)をする。ここで、「ホット」ピクセルとは、ある明度閾値を超えたピクセルと定義される。ステップ912は、ピクセルのグループを見つける。このピクセルのグループは、先に選択された最小数隣接ホットピクセル(以下、スポットまたは明るいスポットと称する)を含み、検出されたスポット群の重心におけるx、y画像座標を決定する。
ステップ916は、開口測光(aperture photometry)を使用することによって各スポットの明度を測定する。
ステップ918は、粒子を計数する。これは、予め選択された閾値より大きい強度を有する画像内の明るいスポット数である。明るいスポットの各々のおよその強度は、そのスポットを含むホットピクセルのすべての強度を加算することによって、計算される。開口測光を使用して、強度についてより正確な測定をすることができる。(すなわち、2以上のスポットを1つの粒子として計数することによって、所定半径内で検出された同一粒子から生ずる同じ粒子を二度計数することを訂正する)。一般に、粒子座標、粒子強度および粒子サイズが取得され、そして当初の画像は消去される。
ステップ922は、全画像が収集されたかどうかを判定する。収集された場合(926)、データ932はブロック822で出力される(図8)。粒子密度930がステップ928において計算され、およびブロック816に出力される。画像収集が完了していない場合(924)、処理はステップ906に戻る。
図10に示すアルゴリズム808Bは、以下のステップを含んでいる。
ステップ1002は、制御信号824aを介してポンプスピードを中間に設定する。ステップ1004は、露光時間およびサンプリング時間を制御信号824bを介して設定する。ポンプ460は、カメラサイクル時間に合致する率で動作される。粒子の各々は撮像された体積内に残留しているので、粒子の各々は、少なくとも1つのカメラ露光時間中に検出される。同時に、対象粒子の運動は、2つの連続した画像の差分において識別できる程度に十分に検出されるであろう。
ステップ1006は、データストリーム404を介してCCD408から一連の画像を取得する。ステップ1008は、前の画像を差し引くことによって、各画像の差分を取り、バックグランドが差し引かれた一連の画像を生成する。負ピクセルがゼロに設定され、差し引かれた画像から対象粒子の影響が取り除かれる。
ステップ1010はホットピクセルをフラグ付けし、ステップ1012は隣接ホットピクセルのグループを検出し、検出された全スポットの重心に関するx、y画像座標を決定する。
ステップ1016は、開口測光を使用することによって、各スポットの明度を測定する。ステップ1018は、強度および/またはサイズ評価基準を使用して、粒子を計数する。ステップ1020は、フレーム内における空間的な同時存在を訂正する。
ステップ1021は、前のフレーム内に検出された粒子に対する粒子数を減らす。粒子が複数のフレーム内に検出される場合、最大明度がその粒子の明度として記録される。
ステップ1022は、全画像が収集されたかどうか判定する。収集されている場合(1026)、データ1030はブロック822(図8)で出力される。粒子密度1032はステップ1028において計算され、ブロック816に出力される。画像収集が完了していない場合(1024)、処理はステップ1006へ戻る。
図11は、高粒子密度条件における図8のアルゴリズムで利用された特定のアルゴリズムを示すフローチャートである。このアルゴリズムは、低い密度および中間的な密度に対するアルゴリズムとは全く異なっている。なぜならば、画像について差分を取ることができないからである。処理は、以下の通りである。
ステップ1102は、制御信号824aを介してポンプスピードを低く設定する。ステップ1104は、制御信号824bを介して露光時間を速い(fast)に設定しおよびサンプリング時間(ポンプスピードおよびサンプリング時間に合致した露光時間)を短い(short)に設定する。
高濃度においては、低ポンプ率が使用されて、焦点面は少数の独立画像を用いてサンプリングされる(すなわち、被写界深度を横切って粒子を運ぶ流速に対して、露出間で十分な時間が与えられる)。
ステップ1106は、データストリーム404を介してCCD408から一連の画像を取得する。
ステップ1108は、画像にメジアンフィルタをかけ、画像中の低空的周波数バックグラウンドが取り除かれる。
ステップ1110はホットピクセルをフラグ付けし、ステップ1112は、隣接ホットピクセルのグループを検出しおよび検出されたスポットの重心に関してx、y画像座標を測定する。ステップ1118は、強度および/またはサイズ評価基準を使用して、粒子を計数する。ステップ1120は、フレーム間における空間的な同時存在を訂正する。
データ1132は、ブロック822(図8)での出力である。粒子密度1130はステップ1128で計算され、ブロック816に出力される。
粒子密度は、撮像された体積によって割られた計数として、1128において計算される。ここで、撮像された体積=被写界深度×照明されおよび撮像された焦点面の面積、である。ソース錯誤に対して、測定された粒子密度が訂正される(以下の数式1)。
10×10スーパーピクセルに分割されるような600×500の画像を見る場合を想定する。このとき、1つのこの画像内に3000スーパーピクセルが含まれることになる。同一のスーパーピクセル内で検出された2つの細胞について混同が生じる場合、すなわち、2つの細胞が1つの細胞として認識される程度に密接している場合を想定する。この画像が1500の細胞画像を含むとき、細胞数が1つだけ増えることは、この細胞がソース錯誤が原因で検出されないだろうという〜50%の可能性を意味する。細胞計数が増えると、更に、細胞濃度に対する非線形な応答が生じる。検出装置のサイズを600×500から1200×1000まで増やすと、ダイナミックレンジが4倍に増大するのは明らかである。
この状況は、放射能測定器にける不感時間統計と類似している。ここで、検知されたn'ピクセル(見かけのカウント)を計数した60×50のスーパーピクセルアレイ内を想定する。n'、「充満」、または「デッド」ピクセル(検出されたピクセル)がある。計数の真の数はnであるべきとする。デッドピクセルの割合はn'/3000である。同時存在(coincidence)の数または同一のスーパーピクセル内でなされた2つの検知のものでないカウント数は、およそn*n'/3000であるだろう。真の計数と見かけの計数との差は、 n-n'= n*n'/3000、またはn = n'/(1-n'/3000)によって与えられる。これは、まさに不感時間の式と同一である。より一般には、以下の通りとなる。
数式 1) n = n'/(1−n'/(α/ΔxΔy))
ここで、αは撮像された体積上で撮像された検出ピクセルの数であり、ΔxおよびΔyはそれぞれx-とy-方向におけるピクセルにおける「識別距離」(2つの粒子として識別するために画像内の2つの粒子の中心間の最小離間距離)である。画像における見かけ計数と真の計数との関係が、図6に示されている。測定される溶液中の真の細胞密度をより良く近似するのに、この関係を使用することができる。真の計数と見かけ計数との関係は非線形であるので、条件3の下での粒子計数は装置のダイナミックレンジを拡大させることに注意すべきである。
数式 1) n = n'/(1−n'/(α/ΔxΔy))
ここで、αは撮像された体積上で撮像された検出ピクセルの数であり、ΔxおよびΔyはそれぞれx-とy-方向におけるピクセルにおける「識別距離」(2つの粒子として識別するために画像内の2つの粒子の中心間の最小離間距離)である。画像における見かけ計数と真の計数との関係が、図6に示されている。測定される溶液中の真の細胞密度をより良く近似するのに、この関係を使用することができる。真の計数と見かけ計数との関係は非線形であるので、条件3の下での粒子計数は装置のダイナミックレンジを拡大させることに注意すべきである。
データ解析(条件番号)のアルゴリズムの選択は、数式 2)−数式 4)で以下に与えるように、流速、被写界深度、露光時間、孔のサイズ、検出形式および粒子密度に依存している。α/ΔxΔy(数式1)と比べて、撮像された体積中の粒子数が小である場合、フレーム間におけるソース錯誤は無視できる。撮像された体積中の粒子数は、粒子密度の撮像された体積倍に等しい。無視できるソース錯誤に対する評価基準(条件1または条件2の適切な使用)は、「静止した」流液に対しては数式2によって与えられる。すなわち、撮像されたボリュームを通過する粒子に対する時間は、露光時間よりはるかに短い。どのアルゴリズム(条件)を用いるべきかを判定するため、流液の「スナップショット」がnを推定するのに使用されるときに、数式2が用いられる。
数式 2) n =ρAd << (α/ΔxΔy) 静止流液中の無視可能なソース錯誤
ここで、nは1つのフレーム中の真の計数であり、ρは粒子密度であり、Aは孔の断面積であり、dは被写界深度である。数式2の不等式が満たされていない場合、条件3を適用でき、そうでなければ、条件1または条件2が適している。時間内にいかなる間隔で撮像された体積中の粒子数が、露光時間の間に撮像されたボリュームを通過する粒子数よりはるかに小である場合、数式2aの「動的流液」近似を適用できる。これは、一般に条件1および条件2の流液のケースである。さらに、n<<(α/ΔxΔy)である場合にのみ、n=n'であり、さもなければ、数式1を使用して訂正をしなければならない。
数式 2a) n =ρAυΔtexp (n<<α/ΔxΔy) 運動流液
ここで、υは流速であり、Δtexpはカメラサイクル時間である。カメラサイクル時間Δtは、露光の長さΔtexpと露光間の不感時間の長さΔtdtとを加えたものに等しい。条件1または条件2が測定に適切であるかどうかが、粒子検出に必要な最小露光時間によって判定される。焦点面(被写界深度/流れ速度)内の粒子の滞在時間がカメラサイクル時間以下である場合、条件1を適用できる。すなわち、
数式 3)Δtexp+Δtdt≦d/2υ 条件1
となる。
数式 2) n =ρAd << (α/ΔxΔy) 静止流液中の無視可能なソース錯誤
ここで、nは1つのフレーム中の真の計数であり、ρは粒子密度であり、Aは孔の断面積であり、dは被写界深度である。数式2の不等式が満たされていない場合、条件3を適用でき、そうでなければ、条件1または条件2が適している。時間内にいかなる間隔で撮像された体積中の粒子数が、露光時間の間に撮像されたボリュームを通過する粒子数よりはるかに小である場合、数式2aの「動的流液」近似を適用できる。これは、一般に条件1および条件2の流液のケースである。さらに、n<<(α/ΔxΔy)である場合にのみ、n=n'であり、さもなければ、数式1を使用して訂正をしなければならない。
数式 2a) n =ρAυΔtexp (n<<α/ΔxΔy) 運動流液
ここで、υは流速であり、Δtexpはカメラサイクル時間である。カメラサイクル時間Δtは、露光の長さΔtexpと露光間の不感時間の長さΔtdtとを加えたものに等しい。条件1または条件2が測定に適切であるかどうかが、粒子検出に必要な最小露光時間によって判定される。焦点面(被写界深度/流れ速度)内の粒子の滞在時間がカメラサイクル時間以下である場合、条件1を適用できる。すなわち、
数式 3)Δtexp+Δtdt≦d/2υ 条件1
となる。
数式3中の1/2の因子は、粒子が絶え間なく焦点面を通して移動しているという事実から生じる。検出において最悪の場合では、2つの連続したフレーム内に、1つの粒子が同じ期間、観察されるであろう。全粒子が露光時間の最小期間に露出されることを確実とするために、粒子がカメラサイクル時間内で移動する距離は、被写界深度の半分未満でなければならない。
1つの目的は、サンプリング時間を最小とすることである。露光時間を、対象粒子の検出のための最小許容時間Δtminに設定する場合、粒子検出においては、粒子が≧2(Δtmin+Δtdt)の撮像された体積中に残存しなければならないことが必要である。
この要件は、(条件1および2に対して)流液中の粒子検出を許容する最大流速度υmaxを順番に決定する。
数式 3a)υmax=d/2(Δtmin+Δtdt) 条件1および2に対する最大流速
一般に、粒子が検出される回数はkであり、ここで、
数式 4) k=d/υ(Δtmin+Δtdt) 1つの粒子の検知回数
である。k>1である場合、条件2を適用できる、そうでなければ、条件1が当てはまる。
この要件は、(条件1および2に対して)流液中の粒子検出を許容する最大流速度υmaxを順番に決定する。
数式 3a)υmax=d/2(Δtmin+Δtdt) 条件1および2に対する最大流速
一般に、粒子が検出される回数はkであり、ここで、
数式 4) k=d/υ(Δtmin+Δtdt) 1つの粒子の検知回数
である。k>1である場合、条件2を適用できる、そうでなければ、条件1が当てはまる。
粒子計数の当業者であれば、好適な実施例の図面および説明が本発明を示すのに有用であるが、多くの他の構成が本発明の精神の範囲内にあるものであることも理解できるだろう。本発明の核心は、画像データを保存し、既定のアルゴリズムに基づいて画像内の粒子密度を推定し、装置パラメータおよび/または計数アルゴリズムを調整して、より正確に計数するという概念にある。
Claims (34)
- 流液(314)中の粒子を計数する方法であって、前記流液は、前記流液の照明断面を撮像する焦点面を有する撮像装置を含む撮像システム(408、412)によってモニタされ、
前記方法は、
(a)前記撮像装置(408)によって1つ以上のフレームを撮像するステップ(404、408)と、
(b)既定の計数アルゴリズムを使用して、撮像されたフレームに基づいて粒子密度の予備的推定を行なうステップ(828)と、
(c)撮像された前記フレームの粒子密度が収まる2つ以上の所定の密度カテゴリを判定するステップ(810)と、
(d)判定された前記密度カテゴリに応じて、
−前記計数アルゴリズム(812)、
−前記撮像システム(824b)のパラメータ、および/または
−流速(824a)
の少なくとも1つを調整するステップと、
(e)前記ステップ(c)の調整をした後に、前記溶液中の粒子を計数するステップ(814)と
を備えることを特徴とする方法。 - 前記調整するステップは、露光時間または照明強度を含む前記撮像システムのパラメータを調整することを特徴とする請求項1に記載の方法。
- 前記調整するステップは、異なる粒子密度カテゴリで動作するように設計された複数の計数アルゴリズムから選択することによって、前記計数アルゴリズムを調整することを特徴とする請求項1に記載の方法。
- 前記計数アルゴリズムは、低位粒子密度アルゴリズム(808A、808B)および高粒子密度アルゴリズム(808C)を含むことを特徴とする請求項3に記載の方法。
- 前記低位粒子密度アルゴリズムは、
2つの連続した画像フレーム(908、1008)の差分を取るステップと、
所定の閾値強度(910、1010)を越えたピクセルのグループをフラグ付けするステップと、
フラグ付けされたグループを粒子(918、1018)として計数するステップと、
空間的な錯誤(920、1020)に対して計数を修正するステップと
を含むことを特徴とする請求項4に記載の方法。 - 粒子を計数する前記ステップは、フラグ付けされたグループのサイズを利用して粒子を識別するステップをさらに含むことを特徴とする請求項5に記載の方法。
- 空間的な錯誤(920、1020)に対して計数を修正する前記ステップは、計数された粒子間の距離に基づいて所定の評価基準に応じて前記計数を減らすことを特徴とする請求項5に記載の方法。
- 前記低位粒子密度アルゴリズムは、低粒子密度アルゴリズム(808A)および中間粒子密度アルゴリズム(808B)に分類されており、前記中間粒子密度アルゴリズムは、どの計数された粒子が先のフレーム(1020)内で計数されていたかを判定するステップと、それに応じて前記計数を判定するステップ(1021)とをさらに含むことを特徴とする請求項5に記載の方法。
- 前記調整するステップは、
前記低粒子密度アルゴリズムにおいては、前記調整するステップは前記ポンプスピードを高く設定しかつ前記露光時間を長く設定し、
前記中間粒子密度アルゴリズムにおいては、前記調整するステップは前記ポンプスピードを中位に設定しかつ前記露光時間を中位に設定し、
前記高粒子密度アルゴリズムにおいては、前記調整するステップは前記ポンプスピードを低く設定しかつ前記露光時間を速く設定して、
前記撮像システムの露光時間およびサンプリング時間を調整することを特徴とする請求項8に記載の方法。 - 前記中間粒子密度アルゴリズムは、連続したフレーム内で検出された運動粒子を1つの粒子(1021)として計数するために、連続した画像のx、y座標を追跡するステップを含むことを特徴とする請求項8に記載の方法。
- 前記高粒子密度アルゴリズムは、
所定の閾値強度(1110、1112)を越えたピクセルのグループをフラグ付けするステップと、
フラグ付けされた粒子(l118)を計数するステップと、
前記計数(1120、1128)に基づいた数式によって、空間的な錯誤訂正因子を前記計数に対して適用するステップとを含むことを特徴とする請求項4の方法。 - 前記フラグ付けするステップ(1108)に先立って、前記フレームに高域通過フィルタをかけるステップをさらに含むことを特徴とする請求項11に記載の方法。
- 前記撮像システムはレンズ(412)を含み、前記撮像装置はCCD画像要素、CMOS撮像要素またはビデオカメラのうちいずれか1つを含むことを特徴とする請求項1の方法。
- 前記撮像システムはレーザまたはLEDのうちいずれか1つを含む照明要素を含むことを特徴とする請求項13に記載の方法。
- 前記流液は、フローサイトメータ内で生ずることを特徴とする請求項1に記載の方法。
- 前記焦点面(234)は、前記流液の方向に対して実質的に垂直であることを特徴とする請求項15に記載の方法。
- 前記焦点面は、前記流液の方向に対して実質的に水平であることを特徴とする請求項15に記載の方法。
- 前記流速が、前記流液の所定の体積内の前記粒子の全てが実質的に1つの画像内で検出されるように調節される(824a)ことを特徴とする請求項1に記載の方法。
- 一連のフレーム群内において対象粒子を追跡するステップをさらに含むことを特徴とする請求項1に記載の方法。
- 前記フレームを撮像するステップは、2色(408、408b)のフレームを撮像し、および前記2色の各々の強度に対して対象粒子の各々を測定するステップをさらに備えることを特徴とする請求項1に記載の方法。
- 撮像されたフレームの明度を既定の状態で測定するステップと、
前記測定された明度に応じて照明または撮像システム露出時間(824b)を調整するステップと
をさらに備えることを特徴とする請求項1に記載の方法。 - 前記計数するステップは、前記計数アルゴリズムを決定する前記画像からのデータを用いて、撮像が完了した後に、一組の撮像されたフレームにおいて実行されることを特徴とする請求項1に記載の方法。
- 粒子の強度分布を測定するステップと、
前記強度分布に応じて照明強度、露出時間およびポンプ速度を調整するステップと
をさらに備えることを特徴とする請求項1に記載の方法。 - 前記流速は、制御可能な流量でシリンジ・ポンプ(460)または蠕動ポンプによって保たれることを特徴とする請求項1に記載の方法。
- 複数のフレーム内で検出された1つの粒子の明度が、最大測定明度となるように設定されることを特徴とする請求項10に記載の方法。
- 流液中の粒子を計数する装置であって、
1つ以上の画像フレームを記録しおよび照明された前記流液(314)の断面を撮像する焦点面を有する撮像装置を含み、前記流液をモニタする撮像システム(412、408)、
プロセッサ内に記憶された既定の計数アルゴリズムを使用して、前記撮像システムからの撮像されたフレーム(404)に基づいて、粒子密度を推定する予備的推定モジュール(806)、
推定された前記粒子密度に基づいて、撮像された前記フレームの前記粒子密度が収まる2つ以上の所定の密度カテゴリを判定する密度カテゴリ化モジュール(810)、
判定された前記密度カテゴリに応じて、前記プロセッサ内に記憶された計数アルゴリズム、撮像システムのパラメータおよび/または流速の中から選択された1つ以上の要素を調整する制御信号(824)を生成する制御システムモジュール(810)および
前記要素が調整された後に、計数アルゴリズムに応じて前記流液中の粒子を計数する計数モジュール(814)
を実装するプロセッサと、
前記制御信号に応じて、選択された要素を調整する調整ユニット(812、408、460)と
を備えたことを特徴とする装置。 - 前記プロセッサは制御信号(824b)を生成して、前記撮像システムのパラメータを調整し、前記調整ユニット(812)は前記制御信号に応じて前記撮像システムの前記露出時間または照明強度を変えることを特徴とする請求項26に記載の装置。
- 前記プロセッサは制御信号(824c)を生成して前記計数アルゴリズムを調整し、前記調整ユニット(812)は前記制御信号に応じて前記計数モジュールにより用いられた前記計数アルゴリズムを変えることを特徴とする請求項26に記載の装置。
- 前記撮像システムはレンズ(412)を含み、前記撮像装置はCCD画像要素、CMOS画像要素またはビデオカメラのうちいずれか1つを含むことを特徴とする請求項26の装置。
- 前記流液は、フローサイトメータ内で生ずることを特徴とする請求項26に記載の装置。
- 前記焦点面は、前記流液の方向に対して実質的に垂直であることを特徴とする請求項30に記載の装置。
- 前記焦点面は、前記流液の方向に対して実質的に平行であることを特徴とする請求項30に記載の装置。
- 前記流速が、前記流液の所定体積内の前記粒子の全てが実質的に1つの画像内で検出されるように調節される(824a)ことを特徴とする請求項26に記載の装置。
- 制御可能な流量で前記流速を決定するシリンジ・ポンプ(460)または蠕動ポンプのうちいずれか一方をさらに備えたことを特徴とする請求項26に記載の装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/328,033 US8798338B2 (en) | 2006-01-09 | 2006-01-09 | Method and system for counting particles in a laminar flow with an imaging device |
PCT/US2007/000470 WO2007081920A2 (en) | 2006-01-09 | 2007-01-08 | Method and system for counting particles in a laminar flow with an imaging device |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2009522586A true JP2009522586A (ja) | 2009-06-11 |
Family
ID=38098580
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008550360A Pending JP2009522586A (ja) | 2006-01-09 | 2007-01-08 | 撮像装置を用いて層流中の粒子を計数する方法とシステム |
Country Status (5)
Country | Link |
---|---|
US (1) | US8798338B2 (ja) |
EP (1) | EP1974198A2 (ja) |
JP (1) | JP2009522586A (ja) |
CA (1) | CA2634238A1 (ja) |
WO (1) | WO2007081920A2 (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013167582A (ja) * | 2012-02-16 | 2013-08-29 | Shimadzu Corp | フローサイトメータ |
JP2018009956A (ja) * | 2016-07-01 | 2018-01-18 | 株式会社リコー | 粒子計数装置及び粒子計数方法 |
JP2018017700A (ja) * | 2016-07-29 | 2018-02-01 | 株式会社リコー | 粒子計数装置及び粒子計数方法、並びに液滴形成装置、及び分注装置 |
WO2019187754A1 (ja) * | 2018-03-30 | 2019-10-03 | シスメックス株式会社 | フローサイトメーター及び粒子検出方法 |
KR102598764B1 (ko) * | 2023-06-26 | 2023-11-06 | 주식회사 큐리오시스 | 교정용 테이블을 이용하여 이미지의 밝기를 보정하는미세입자 계수 장치 및 방법 |
Families Citing this family (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7340957B2 (en) | 2004-07-29 | 2008-03-11 | Los Alamos National Security, Llc | Ultrasonic analyte concentration and application in flow cytometry |
US7835000B2 (en) * | 2006-11-03 | 2010-11-16 | Los Alamos National Security, Llc | System and method for measuring particles in a sample stream of a flow cytometer or the like |
JP2010509610A (ja) * | 2006-11-13 | 2010-03-25 | ユニバーシティ オブ ワイオミング | 多細胞生物の蛍光測定装置及びその蛍光測定方法 |
US8068673B2 (en) * | 2006-12-01 | 2011-11-29 | Beihang University | Rapid and high precision centroiding method and system for spots image |
US7837040B2 (en) | 2007-04-09 | 2010-11-23 | Los Alamos National Security, Llc | Acoustic concentration of particles in fluid flow |
US8083068B2 (en) | 2007-04-09 | 2011-12-27 | Los Alamos National Security, Llc | Apparatus for separating particles utilizing engineered acoustic contrast capture particles |
JP5324598B2 (ja) * | 2007-12-04 | 2013-10-23 | パーティクル・メージャーリング・システムズ・インコーポレーテッド | 非直角粒子検出システム及び方法 |
US8266950B2 (en) | 2007-12-19 | 2012-09-18 | Los Alamos National Security, LLP | Particle analysis in an acoustic cytometer |
WO2010104993A2 (en) | 2009-03-10 | 2010-09-16 | The Regents Of The University Of California | Fluidic flow cytometry devices and particle sensing based on signal-encoding |
US9645010B2 (en) * | 2009-03-10 | 2017-05-09 | The Regents Of The University Of California | Fluidic flow cytometry devices and methods |
EP2439511B1 (en) * | 2009-06-03 | 2015-08-19 | Hitachi High-Technologies Corporation | Flow type particle image analysis method and device |
US9500645B2 (en) | 2009-11-23 | 2016-11-22 | Cyvek, Inc. | Micro-tube particles for microfluidic assays and methods of manufacture |
US9700889B2 (en) * | 2009-11-23 | 2017-07-11 | Cyvek, Inc. | Methods and systems for manufacture of microarray assay systems, conducting microfluidic assays, and monitoring and scanning to obtain microfluidic assay results |
US10022696B2 (en) | 2009-11-23 | 2018-07-17 | Cyvek, Inc. | Microfluidic assay systems employing micro-particles and methods of manufacture |
US10065403B2 (en) | 2009-11-23 | 2018-09-04 | Cyvek, Inc. | Microfluidic assay assemblies and methods of manufacture |
KR101885936B1 (ko) | 2010-10-21 | 2018-09-10 | 더 리젠츠 오브 더 유니버시티 오브 캘리포니아 | 무선 전력이 공급되는 전자 회로를 포함하는 미세 유체 |
WO2012098287A1 (en) * | 2011-01-19 | 2012-07-26 | Teknologian Tutkimuskeskus Vtt | Method and system for determining particle size information |
CA2830533C (en) | 2011-03-22 | 2020-02-18 | Cyvek, Inc. | Microfluidic devices and methods of manufacture and use |
CN102609762A (zh) * | 2012-01-10 | 2012-07-25 | 林美祥 | 一种led芯片数量统计方法 |
TW201404878A (zh) * | 2012-07-27 | 2014-02-01 | Hsian-Chang Chen | 自動快速分析生物細胞的裝置和其相關的方法 |
EP2906928A4 (en) | 2012-10-15 | 2016-11-09 | Nanocellect Biomedical Inc | SYSTEMS, APPARATUS AND METHODS FOR SORTING PARTICLES |
EP3108219B1 (en) * | 2014-02-20 | 2020-04-01 | Malvern Panalytical Limited | Heterogeneous fluid sample characterization |
US10151677B2 (en) | 2014-07-08 | 2018-12-11 | Halliburton Energy Services, Inc. | Real-time optical flow imaging to determine particle size distribution |
US10801959B2 (en) * | 2015-06-18 | 2020-10-13 | Konica Minolta, Inc. | Objective biological substance analysis device, analysis system, analysis method, and recording medium storing computer readable analysis program |
US10228367B2 (en) | 2015-12-01 | 2019-03-12 | ProteinSimple | Segmented multi-use automated assay cartridge |
US10241024B1 (en) * | 2016-05-02 | 2019-03-26 | Rashid Mavliev | System and method for characterization of inclusions in liquid samples |
BE1023800B1 (fr) * | 2016-06-24 | 2017-07-26 | Occhio | Dispositif optique de mesure d'une charge en particules d'un echantillon |
US10809176B2 (en) | 2016-10-24 | 2020-10-20 | The Trustees Of The University Of Pennsylvania | Ultra-high throughput detection of fluorescent droplets using time domain encoded optofluidics |
WO2018227193A1 (en) * | 2017-06-09 | 2018-12-13 | miDiagnostics NV | Systems and methods for counting particles |
US11195299B2 (en) * | 2017-09-12 | 2021-12-07 | Nec Corporation | Image processing device, image processing method and storage medium |
EP3859007A4 (en) * | 2018-09-28 | 2021-11-24 | FUJIFILM Corporation | DETERMINATION PROCEDURE |
CN112752964A (zh) * | 2018-10-01 | 2021-05-04 | 惠普发展公司,有限责任合伙企业 | 批量颗粒分选 |
KR101970220B1 (ko) * | 2018-11-20 | 2019-04-18 | 주식회사 제덱스 | 여과재를 이용한 유체 중의 실시간 이물 계수 장치 |
US10876949B2 (en) * | 2019-04-26 | 2020-12-29 | Honeywell International Inc. | Flow device and associated method and system |
US11835432B2 (en) | 2020-10-26 | 2023-12-05 | Honeywell International Inc. | Fluid composition sensor device and method of using the same |
US20220243246A1 (en) * | 2021-02-01 | 2022-08-04 | Arizona Board Of Regents On Behalf Of Arizona State University | Rapid antimicrobial susceptibility testing by video-based object scattering intensity detection |
WO2024077048A1 (en) * | 2022-10-06 | 2024-04-11 | Regents Of The University Of Minnesota | Suspended particle detection and analysis |
CN116698709B (zh) * | 2023-06-09 | 2024-03-12 | 深圳市益希医疗器械有限公司 | 一种流式细胞仪的数据处理方法及流式细胞仪 |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4075462A (en) * | 1975-01-08 | 1978-02-21 | William Guy Rowe | Particle analyzer apparatus employing light-sensitive electronic detector array |
US4112301A (en) * | 1976-06-03 | 1978-09-05 | American Science And Engineering, Inc. | Moving particles suspended in a carrier fluid through a flow channel having an input end under gas pressure |
US4814868A (en) * | 1987-10-02 | 1989-03-21 | Quadtek, Inc. | Apparatus and method for imaging and counting moving particles |
US5030002A (en) * | 1989-08-11 | 1991-07-09 | Becton, Dickinson And Company | Method and apparatus for sorting particles with a moving catcher tube |
US5249238A (en) * | 1991-03-19 | 1993-09-28 | Komerath Narayanan M | Spatial cross-correlating velocimeter |
US5978435A (en) * | 1995-08-21 | 1999-11-02 | Foss Electric A/S | Method and an apparatus for determining the number of particles or cells in a liquid sample |
JP3127111B2 (ja) * | 1996-02-22 | 2001-01-22 | 株式会社日立製作所 | フロー式粒子画像解析方法および装置 |
EP1180676A3 (en) * | 1997-05-05 | 2003-05-21 | ChemoMetec A/S | A method and a system for determination of particles in a liquid sample |
US6731100B1 (en) * | 1997-05-05 | 2004-05-04 | Chemometec A/S | Method and a system for determination of somatic cells in milk |
DE19735066C1 (de) * | 1997-08-13 | 1999-01-28 | Hydac Filtertechnik Gmbh | Auswerteverfahren für einen Partikelzähler und Vorrichtung zum Durchführen des Verfahrens |
US6198839B1 (en) * | 1997-09-05 | 2001-03-06 | Tripath Imaging, Inc. | Dynamic control and decision making method and apparatus |
US6115119A (en) * | 1997-10-21 | 2000-09-05 | Bigelow Laboratory For Ocean Science | Device and method for studying particles in a fluid |
US6309886B1 (en) * | 1999-06-04 | 2001-10-30 | The Regents Of The University Of California | High throughput analysis of samples in flowing liquid |
AU2002364580A1 (en) * | 2001-12-18 | 2003-06-30 | University Of Wyoming | Apparatus and methods for high throughput analysis of particles in a translucent flowing liquid |
US7161665B2 (en) * | 2002-12-18 | 2007-01-09 | University Of Wyoming | High resolution imaging fountain flow cytometry |
AU2004254632B2 (en) * | 2003-07-02 | 2009-02-12 | Terumo Bct, Inc. | Monitoring and control system for blood processing |
GB0519323D0 (en) * | 2005-09-22 | 2005-11-02 | Bann John R | Scattering centre detector assembly and method |
KR101165162B1 (ko) * | 2008-04-25 | 2012-07-11 | 신닛뽄세이테쯔 카부시키카이샤 | 금속 재료 중 미립자의 입도 분포 측정 방법 |
CN102782480B (zh) * | 2010-03-01 | 2013-09-11 | 奥林巴斯株式会社 | 光学分析装置、光学分析方法和用于光学分析的计算机程序 |
-
2006
- 2006-01-09 US US11/328,033 patent/US8798338B2/en active Active
-
2007
- 2007-01-08 EP EP07716437A patent/EP1974198A2/en not_active Withdrawn
- 2007-01-08 JP JP2008550360A patent/JP2009522586A/ja active Pending
- 2007-01-08 CA CA002634238A patent/CA2634238A1/en not_active Abandoned
- 2007-01-08 WO PCT/US2007/000470 patent/WO2007081920A2/en active Application Filing
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013167582A (ja) * | 2012-02-16 | 2013-08-29 | Shimadzu Corp | フローサイトメータ |
JP2018009956A (ja) * | 2016-07-01 | 2018-01-18 | 株式会社リコー | 粒子計数装置及び粒子計数方法 |
JP2018017700A (ja) * | 2016-07-29 | 2018-02-01 | 株式会社リコー | 粒子計数装置及び粒子計数方法、並びに液滴形成装置、及び分注装置 |
WO2019187754A1 (ja) * | 2018-03-30 | 2019-10-03 | シスメックス株式会社 | フローサイトメーター及び粒子検出方法 |
JP2019178900A (ja) * | 2018-03-30 | 2019-10-17 | シスメックス株式会社 | フローサイトメーター及び粒子検出方法 |
JP7109229B2 (ja) | 2018-03-30 | 2022-07-29 | シスメックス株式会社 | フローサイトメーター及び粒子検出方法 |
KR102598764B1 (ko) * | 2023-06-26 | 2023-11-06 | 주식회사 큐리오시스 | 교정용 테이블을 이용하여 이미지의 밝기를 보정하는미세입자 계수 장치 및 방법 |
Also Published As
Publication number | Publication date |
---|---|
EP1974198A2 (en) | 2008-10-01 |
WO2007081920A3 (en) | 2007-09-20 |
CA2634238A1 (en) | 2007-07-19 |
US8798338B2 (en) | 2014-08-05 |
US20070159627A1 (en) | 2007-07-12 |
WO2007081920A2 (en) | 2007-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8798338B2 (en) | Method and system for counting particles in a laminar flow with an imaging device | |
US7307721B2 (en) | Particle imaging system with a varying flow rate | |
US7161665B2 (en) | High resolution imaging fountain flow cytometry | |
CA2395627C (en) | Imaging and analyzing parameters of small moving objects such as cells | |
CA2299325C (en) | A digital imaging system for assays in well plates, gels and blots | |
JP3290786B2 (ja) | 粒子分析装置 | |
US11125675B2 (en) | Fluid suspended particle classifier | |
Alvarez et al. | SP8 FALCON: a novel concept in fluorescence lifetime imaging enabling video-rate confocal FLIM | |
US11874228B2 (en) | Methods for identification of particles in a fluid sample | |
JP2002535614A (ja) | 改良したフローサイトメトリの装置および方法 | |
US20210010925A1 (en) | Flow cytometer and particle detection method | |
US7217937B2 (en) | Automatic identification of suspended particles | |
US20120002029A1 (en) | System and method for sweeping mirror enhanced imaging flow cytometry | |
US20110017915A1 (en) | Drift scanner for rare cell detection | |
JP3984132B2 (ja) | 蛍光分光分析装置 | |
US20090290780A1 (en) | Analysis method for chemical and/or biological samples | |
Balsam et al. | Cell streak imaging cytometry for rare cell detection | |
JP2008116422A (ja) | 微粒子検出装置 | |
WO2022080482A1 (ja) | 粒子検出装置、粒子検出システム、及び粒子検出方法 | |
US20230417671A1 (en) | Biological sample analysis system, information processing device, information processing method, and biological sample analysis method | |
CN114778419A (zh) | 一种高倍率光学放大成像流式细胞仪 | |
EP3922978A1 (en) | Optical device for detecting fluorescence emission | |
WO2023189819A1 (ja) | 粒子分取システム、及び粒子分取方法 | |
WO2023223851A1 (ja) | 生体試料分析システム、情報処理装置、情報処理方法及び生体試料分析方法 | |
WO2023276298A1 (ja) | 生体試料分析システム、情報処理装置、情報処理方法及び生体試料分析方法 |