JP4982697B2 - 粒子計測装置 - Google Patents
粒子計測装置 Download PDFInfo
- Publication number
- JP4982697B2 JP4982697B2 JP2007104931A JP2007104931A JP4982697B2 JP 4982697 B2 JP4982697 B2 JP 4982697B2 JP 2007104931 A JP2007104931 A JP 2007104931A JP 2007104931 A JP2007104931 A JP 2007104931A JP 4982697 B2 JP4982697 B2 JP 4982697B2
- Authority
- JP
- Japan
- Prior art keywords
- particle
- particles
- unit
- measuring
- image
- 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
Links
- 239000002245 particle Substances 0.000 title claims description 304
- 239000012530 fluid Substances 0.000 claims description 36
- 238000005259 measurement Methods 0.000 claims description 30
- 238000003384 imaging method Methods 0.000 claims description 17
- 230000006399 behavior Effects 0.000 claims description 13
- 239000007787 solid Substances 0.000 claims description 11
- 239000007788 liquid Substances 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 21
- 238000000034 method Methods 0.000 description 6
- 210000001772 blood platelet Anatomy 0.000 description 5
- 239000004793 Polystyrene Substances 0.000 description 4
- 229920002223 polystyrene Polymers 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 210000004369 blood Anatomy 0.000 description 3
- 239000008280 blood Substances 0.000 description 3
- 210000000601 blood cell Anatomy 0.000 description 3
- 238000004220 aggregation Methods 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 210000003743 erythrocyte Anatomy 0.000 description 2
- 238000001033 granulometry Methods 0.000 description 2
- 230000004931 aggregating effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 229910052736 halogen Inorganic materials 0.000 description 1
- 150000002367 halogens Chemical class 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 210000000265 leukocyte Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000002381 plasma Anatomy 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 229910052724 xenon Inorganic materials 0.000 description 1
- FHNFHKCVQCLJFQ-UHFFFAOYSA-N xenon atom Chemical compound [Xe] FHNFHKCVQCLJFQ-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Description
粒子が混入された流体を層流状態で流下させる微小流路(2a)と、
該流体の流下方向(x)に略直交する方向から、前記微小流路(2a)を流下される粒子を撮影する撮影手段(6)と、
該撮影手段(6)が撮影した画像データ(GD)に基づき、流下される粒子の速度(u)を計測する粒子速度計測手段(70)と、
前記撮影手段(6)が撮影した画像データ(GD)に基づき、所定時間内に流下される粒子の数(n)を計測する粒子数計測手段(71)と、
前記撮影手段(6)が撮影した画像データ(GD)に基づき、流下される粒子の径を計測する粒径計測手段(72)と、
各計測手段(70,71,72)が計測を行なったときの時間のデータに基づき、各計測手段(70,71,72)が計測した計測値を時間に関連付けて管理するデータ関連付け手段(73)と、を備え、
前記粒子数計測手段(71)は、
前記撮影手段(6)が撮影した2枚の静止画像から差分画像を形成する差分画像形成部(711)と、
どの粒子がどの位置に移動したかの対応付けを行なう対応付け部(712)と、
該対応付け部(712)からのデータに基づき粒子の数(n)を計数する粒子数計数部(713)と、
前記粒子速度計測手段(70)で求めた各粒子の速度(u)から、次の静止画像における粒子の位置を予測する粒子位置予測部(714)と、
該粒子位置予測部(714)のデータに基づき、次の静止画像における粒子の挙動であるオーバーラップやフレーム・インやフレーム・アウトを予測する粒子挙動予測部(715)と、
該粒子挙動予測部(715)からのデータに基づき前記粒子数計数部(713)の計測データを補正する粒子数補正部(716)と、
を備え、
前記フレーム・インとは、粒子が1枚目の静止画像に写っておらず、次の静止画像で現れることをいい、
前記フレーム・アウトとは、粒子が1枚目の静止画像に写っていても、次の静止画像で写らなくなることをいい、
前記オーバーラップとは、粒子が、他の粒子に追いついて追い越すことをいうことを特徴とする。
・ 所定時間内に流下される粒子の数と、
・ 粒子の速度と、
・ 粒子の径
とを計測するものである。
・ 液体、
・ 気体、
・ 液体−気体の混相流体、
・ 液体−固体の混相流体、
・ 気体−固体の混相流体、又は
・ 液体−気体−固体の混相流体
を言う。本発明は、例えば、血液中の血小板の粒径を計測するために用いることができる。その血液中には血漿や血小板や赤血球や白血球が存在するが、その場合には、血漿が前記液体に該当し、血小板が前記粒子に該当し、赤血球が前記固体に該当することとなり、血液自体は液体−気体−固体の混相流体に該当することとなる。
(1)
前記撮影手段6が撮影した画像データGDに基づき、流下される粒子の速度uを計測する粒子速度計測手段70と、
(2)
該撮影手段6が撮影した画像データGDに基づき、所定時間内に流下される粒子の数nを計測する粒子数計測手段71と、
(3)
前記撮影手段6が撮影した画像データGDに基づき、流下される粒子の径dp3を計測する粒径計測手段72と、
(4)
各計測手段70,71,72が計測を行なったときの時間のデータに基づき、各計測手段70,71,72が計測した計測値u,n,dp3を時間に関連付けて管理するデータ関連付け手段73と、
が配置されている。以下、これらの手段70,71,72,73について詳述する。
微小流路2aでは流体は上述のように層流状態に流れるので、各粒子は微小流路2aの壁面との距離がほぼ一定のままで流下され(つまり、壁面に近づいたり遠ざかったりすることはほとんど無く)、一つ一つの粒子の速度自体はほとんど変化しない。したがって、前記粒子速度計測手段70は、
・ 前記撮影手段6が撮影した複数枚の静止画像を画像処理して各粒子の移動距離を算出する移動距離算出部701と、
・ 該算出した移動距離を撮影間隔Δtで割ることによりその速度を算出する粒子速度算出部702と、
により構成すると良い。前記移動距離算出部701は、例えば、後述する差分画像形成部711と同様に差分画像を形成し、公知の画像処理を施すことによって移動距離を算出すると良い。
この粒子数計測手段71は、
・ 前記撮影手段6が撮影した2枚の静止画像(画像データ)から差分画像を形成する差分画像形成部711と、
・ どの粒子がどの位置に移動したかの対応付けを行なう対応付け部712と、
・ 該対応付け部712からのデータに基づき、粒子の数を計数する粒子数計数部713と、
により構成すると良い。ここで、図3(a)
(b) は、撮影手段6が撮影した2枚の静止画像の一例を示すものである。図3(a) に示す4つの粒子は、所定時間後には同図(b) に示す位置まで移動している。また、図3(c)
は、それら2枚の静止画像の差分画像であって、適当なオフセットを加えたものである。その差分画像に表示される白い点は、後の画像で撮影された粒子であり、黒い点は、前の画像で撮影された粒子である。いま、バックグラウンドの灰色よりも黒い色と、該灰色よりも白い色の2つの閾値を使って画像処理を行ない、黒い点及び白い点を抽出する。粒子を抽出後、前記対応付け部712によってPIV(Particle Image Velocimetry:粒子画像流速測定法)処理を行なって粒子の対応付けを行ない、前記粒子数計数部713により粒子数の計数を行なう。
上述したように流体は層流状態で流下されて一つ一つの粒子の速度自体はほとんど変化しない。しかし、層流状態で流れる流体の流速そのものは、流路壁面からの距離によって異なり、粒子の速度も、流路壁面からの距離によって異なる。図5は、粒子の速度と第3軸方向位置との関係の一例を示す図であり、第3軸方向位置が0及び50μmは流路壁面で、第3軸方向位置が25μmでは流路の真ん中ということになる。つまり、粒子の速度は壁面から離れるほど速くなり、上に凸の放物線のような曲線を描くこととなる。いま、図6に示す微小流路2aにおいて仮想面Bを考えると、この仮想面Bに沿って流下される粒子について図5の関係が成立し(当然ながら、該仮想面Bに直交するxy面についても図5の関係は成立する)、速度分布曲線は符号Cで示すようになり、壁面に近いほど粒子速度が遅くなり、壁面から離れるほど粒子速度は速くなる。ここで、
・ 流路上壁面に近い位置を流れる粒子2
・ 該上壁面から離れた位置を流れる粒子1
の2つの粒子について考察すると、上側の粒子2の速度は下側の粒子1の速度に比べると遅くなることが分る。
に示す粒子1は、粒子2よりも速度が速いので、次の撮影画像(同図(b) 参照)では粒子2に追いつき、さらに次の撮影画像(同図(c) 参照)では該粒子2を追い超すこととなる。
・ 前記粒子速度計測手段70で求めた各粒子の速度uから、次の静止画像における粒子の位置を予測する粒子位置予測部714と、
・ 該粒子位置予測部714のデータに基づき粒子の挙動(粒子同士の重なりやオーバーラップやフレーム・インやフレーム・アウト等)を予測する粒子挙動予測部715と、
・ 該粒子挙動予測部715からのデータに基づき前記粒子数計数部713の計測データを補正する粒子数補正部716と、
を設けておくと良い。
・ 粒子についての第2軸方向に関する位置yを前記画像データGDから測定する第2軸方向位置測定部723と、
・ 前記流量測定部4の測定結果Q、前記粒子速度計測手段70の算出結果u、前記第2軸方向位置測定部723による測定結果y、及び前記微小流路2aを流れる流体についてのナビエ・ストークス(Navier−Stokes)の方程式から粒子の第3軸方向位置zを算出する第3軸方向位置算出部724と、
・ 該第3軸方向位置算出部724の算出結果zに基づき、前記撮影手段6の焦点面60aを基準とした粒子の第3軸方向についてのズレ量Δzを算出するズレ量算出部725と、
・ 該ズレ量算出部725の算出結果Δzと前記第1粒径算出手段72Aの算出結果de’とから、ズレ量Δzに伴う補正を行なった第2粒径dp1を算出する第2粒径算出部726と、
を備えている。
)との関係は図5に示すようになり、粒子速度計測手段70の演算結果uに基づき第3軸方向位置zが算出されることとなる。
・ 径が2μmのポリスチレン粒子を純水に混入したもの
・ 径が5μmのポリスチレン粒子を純水に混入したもの
・ 径が10μmのポリスチレン粒子を純水に混入したもの
・ 径が20μmのポリスチレン粒子を純水に混入したもの
を用意し、前記粒径計測手段72により各試料の粒径を計測した。その結果は、図14に示すようになり、本装置の計測精度が高いことが分った。
2a 微小流路
6 撮影手段
70 粒子速度計測手段
71 粒子数計測手段
72 粒径計測手段
72A 第1粒径算出手段
72B 粒径補正部
73 データ関連付け手段
701 移動距離算出部
702 粒子速度算出部
711 差分画像形成部
712 対応付け部
713 粒子数計数部
714 粒子位置予測部
715 粒子挙動予測部
716 粒子数補正部
GD 画像データ
Claims (5)
- 画像データに基づいて粒子を計測する粒子計測装置において、
粒子が混入された流体を層流状態で流下させる微小流路と、
該流体の流下方向に略直交する方向から、前記微小流路を流下される粒子を撮影する撮影手段と、
該撮影手段が撮影した画像データに基づき、流下される粒子の速度を計測する粒子速度計測手段と、
前記撮影手段が撮影した画像データに基づき、所定時間内に流下される粒子の数を計測する粒子数計測手段と、
前記撮影手段が撮影した画像データに基づき、流下される粒子の径を計測する粒径計測手段と、
各計測手段が計測を行なったときの時間のデータに基づき、各計測手段が計測した計測値を時間に関連付けて管理するデータ関連付け手段と、
を備え、
前記粒子数計測手段は、
前記撮影手段が撮影した2枚の静止画像から差分画像を形成する差分画像形成部と、
どの粒子がどの位置に移動したかの対応付けを行なう対応付け部と、
該対応付け部からのデータに基づき粒子の数を計数する粒子数計数部と、
前記粒子速度計測手段で求めた各粒子の速度から、次の静止画像における粒子の位置を予測する粒子位置予測部と、
該粒子位置予測部のデータに基づき、次の静止画像における粒子の挙動であるオーバーラップやフレーム・インやフレーム・アウトを予測する粒子挙動予測部と、
該粒子挙動予測部からのデータに基づき前記粒子数計数部の計測データを補正する粒子数補正部と、
を備え、
前記フレーム・インとは、粒子が1枚目の静止画像に写っておらず、次の静止画像で現れることをいい、
前記フレーム・アウトとは、粒子が1枚目の静止画像に写っていても、次の静止画像で写らなくなることをいい、
前記オーバーラップとは、粒子が、他の粒子に追いついて追い越すことをいう、
ことを特徴とする粒子計測装置。 - 前記粒子速度計測手段は、
前記撮影手段が撮影した複数枚の静止画像を画像処理して各粒子の移動距離を算出する移動距離算出部と、
該算出した移動距離を撮影間隔で割ることによりその速度を算出する粒子速度算出部と、
からなることを特徴とする請求項1に記載の粒子計測装置。 - 前記粒径計測手段は、前記撮影手段が撮影した画像データに基づき、画像上での粒径を算出する第1粒径算出手段を備えた、
ことを特徴とする請求項1又は2に記載の粒子計測装置。 - 前記粒径計測手段は、前記撮影手段の焦点面を基準とした粒子のズレ量に基づき前記第1粒径算出手段の算出結果を補正する粒径補正部を備えた、
ことを特徴とする請求項1乃至3のいずれか1項に記載の粒子計測装置。 - 前記流体は、液体、気体、液体−気体の混相流体、液体−固体の混相流体、気体−固体の混相流体、又は液体−気体−固体の混相流体である、
ことを特徴とする請求項1乃至4のいずれか1項に記載の粒子計測装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007104931A JP4982697B2 (ja) | 2007-04-12 | 2007-04-12 | 粒子計測装置 |
US12/595,348 US8218132B2 (en) | 2007-04-12 | 2008-04-08 | Particle measuring device and particle size measure device |
PCT/JP2008/056939 WO2008132995A1 (ja) | 2007-04-12 | 2008-04-08 | 粒子計測装置、及び粒径計測装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007104931A JP4982697B2 (ja) | 2007-04-12 | 2007-04-12 | 粒子計測装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2008261740A JP2008261740A (ja) | 2008-10-30 |
JP4982697B2 true JP4982697B2 (ja) | 2012-07-25 |
Family
ID=39984325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007104931A Expired - Fee Related JP4982697B2 (ja) | 2007-04-12 | 2007-04-12 | 粒子計測装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4982697B2 (ja) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016063912A1 (ja) * | 2014-10-24 | 2016-04-28 | 国立大学法人東京大学 | 粒子検出方法、粒子検出装置および粒子検出システム |
WO2016171198A1 (ja) * | 2015-04-21 | 2016-10-27 | 国立大学法人東京大学 | 微粒子検出システム及び微粒子検出プログラム |
JP2019100988A (ja) * | 2017-12-08 | 2019-06-24 | 株式会社島津製作所 | 粒子画像解析装置及び粒子画像解析方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3096099B2 (ja) * | 1991-08-13 | 2000-10-10 | 科学技術振興事業団 | 粒子測定装置 |
JP3260469B2 (ja) * | 1992-04-01 | 2002-02-25 | シスメックス株式会社 | 粒子分析装置 |
JPH0989753A (ja) * | 1995-09-25 | 1997-04-04 | Hitachi Ltd | 粒子分析装置 |
JP2007010524A (ja) * | 2005-06-30 | 2007-01-18 | Univ Of Tokyo | 微小流路内速度分布計測装置及び方法 |
-
2007
- 2007-04-12 JP JP2007104931A patent/JP4982697B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2008261740A (ja) | 2008-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8218132B2 (en) | Particle measuring device and particle size measure device | |
JP5146462B2 (ja) | 血液流動性計測システム | |
TWI674401B (zh) | 粒子測量裝置和粒子測量方法 | |
CN204116189U (zh) | 智能光微流成像及检测系统 | |
CN105765440A (zh) | 用于显微的样品处理改进 | |
Gu et al. | Simultaneous vision-based shape and motion analysis of cells fast-flowing in a microchannel | |
JP6497039B2 (ja) | 粒子画像流速測定方法および粒子画像流速測定装置 | |
JP2005331506A (ja) | マイクロチップ、及び該マイクロチップを備えた蛍光粒子計数装置 | |
CN108169510B (zh) | 基于单光场相机的微尺度流动三维速度场测量装置和方法 | |
Dehaeck et al. | Laser marked shadowgraphy: a novel optical planar technique for the study of microbubbles and droplets | |
JP4982697B2 (ja) | 粒子計測装置 | |
WO2021115239A1 (zh) | 藻类分析方法 | |
JP2004069431A (ja) | 液中粒子の画像解析装置 | |
CN103675333A (zh) | 一种实时测量微流体速度场的装置及方法 | |
JP6726687B2 (ja) | 粒子分析装置及び粒子分析方法 | |
US10168196B2 (en) | Immersion-type online multiphase measuring instrument and method | |
JP4863083B2 (ja) | 粒径計測装置 | |
CN111458311B (zh) | 固着液滴的内部三维瞬态速度场获取方法 | |
TW202028716A (zh) | 粒子測定裝置、校正方法以及測定裝置 | |
CN105806780B (zh) | 基于多视觉的小通道气液两相流相含率测量系统的相含率测量方法 | |
US20220364970A1 (en) | Particle measuring device and particle measuring method | |
CN110261642A (zh) | 适用于气液界面的立体式粒子图像测速方法 | |
Paz et al. | On the application of image processing methods for bubble recognition to the study of subcooled flow boiling of water in rectangular channels | |
He et al. | Full-field convection flow visualization in pendant droplets by tilt-angle imaging | |
CN109870401A (zh) | 一种流式细胞仪及超分辨率细胞图像的采集方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20100305 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110927 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20111116 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20120306 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20120321 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20120402 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20120402 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150511 Year of fee payment: 3 |
|
LAPS | Cancellation because of no payment of annual fees |