(57)【要約】【課題】マイクロ流体チップによって製造された液滴の粒径解析が簡単かつ迅速に行えるようにする。【解決手段】動画データの入力部１、テンプレート画像格納部２、動画データを前処理する前処理部３、動画データの各フレームに対する検出エリアとエッジ検出用検出線分を設定する検出エリア設定部４及びテンプレート画像を用いて検出エリア内でパターンマッチングを行う判定部５を備える。検出エリア内に新たな分散相が検出される度に分散相を計数する分散相計数部６と、検出エリア内に新たな分散相が検出される度に関係する１つのフレームに対しエッジ検出を行い、分散相長さを測定する分散相長さ計測部７と、分散相長さの測定値とマイクロ流路の断面積値に基づき分散相を球状とみなした場合の粒径を算出する粒径算出部８と、分散相の粒径の総和を分散相計数部の計数値で除算し、分散相の平均粒径を算出する平均粒径算出部９を備える。【選択図】図１
2018169224.tif 000002
【特許請求の範囲】
【請求項１】
マイクロ流路内を流れるエマルションが撮影された動画データに基づいて、前記エマルションの分散相のサイズを計測するシステムであって、
前記動画データ、および前記マイクロ流路の断面積値の入力を受ける入力部と、
予め登録された前記分散相の参照画像が格納された参照画像格納部と、
前記入力部に入力された前記動画データのキャリブレーションを行うとともに、前記動画データがカラー画像データである場合は、前記動画データをグレースケール変換してモノクロ画像データにする前処理部と、
前記前処理部で得られた前記動画データの各フレームに対して固定した矩形状の検出エリアを設定し、前記検出エリア内にはエッジ検出のための検出線分を設定する検出エリア設定部と、を備え、
前記検出エリアは前記マイクロ流路の直線状部分の一定長をカバーし、かつ前記検出エリアの長辺が前記マイクロ流路に平行となるように配置され、前記検出線分は前記マイクロ流路の中心軸に一致するように配置され、前記検出エリアの長辺が同一フレーム内の異なる２つの前記分散相が同時に前記検出エリア内に位置し得ないような長さに設定されており、さらに、
前記参照画像を用いて、前記前処理部で得られた前記動画データの各フレームに対しフレーム番号順に前記検出エリア内においてパターンマッチングを行い、前記検出エリア内に前記参照画像に一致する画像が存在する前記フレームを分散相検出フレームとする一方、前記検出エリア内に前記参照画像に一致する画像が存在しない前記フレームを分散相不検出フレームとする判定部と、
前記判定部において前記分散相不検出フレームの次に前記分散相検出フレームが出現した後別の前記分散相不検出フレームが出現するたびに、その都度、当該分散相不検出フレームおよび当該別の分散相不検出フレーム間に位置する前記分散相検出フレームのうちの１つに対し前記エッジ検出を行って、関係する前記分散相の前後端縁間の距離を測定する分散相長さ計測部と、
前記分散相長さ計測部による測定値と、前記入力部に入力された前記マイクロ流路の断面積値とに基づき、関係する前記分散相のサイズを算出して出力する分散相サイズ算出部と、を備えたものであることを特徴とするシステム。
【請求項２】
前記分散相サイズ算出部は、前記分散相長さ計測部による前記測定値と、前記入力部に入力された前記マイクロ流路の断面積値とに基づき、関係する前記分散相を球状とみなしたときの粒径を算出して出力する粒径算出部を有していることを特徴とする請求項１に記載のシステム。
【請求項３】
前記判定部において前記分散相不検出フレームの次に前記分散相検出フレームが出現するたびに、その都度、前記分散相の個数をカウントする分散相計数部と、
前記粒径算出部から出力された前記分散相の粒径の総和を、前記分散相計数部によるカウント値によって割り算することにより、前記分散相の平均粒径を算出して出力する平均粒径算出部と、を備えたものであることを特徴とする請求項２に記載のシステム。
【請求項４】
前記入力部が前記動画データのフレームレートの入力を受けるようになっており、
前記判定部において前記分散相不検出フレームの次に前記分散相検出フレームが出現した後、別の前記分散相不検出フレームの次に別の前記分散相検出フレームが出現するたびに、その都度、当該分散相不検出フレームおよび当該別の分散相検出フレーム間に位置するフレームの数を算出するフレーム数算出部と、
前記フレーム数算出部による前記フレームの数の算出値と、前記入力部に入力された前記フレームレートとに基づき分散相生成速度を算出して出力する分散相生成速度算出部と、を備えたものであることを特徴とする請求項１〜請求項３のいずれかに記載のシステム。
【請求項５】
前記判定部において前記分散相検出フレームが出現するたびに、その都度、前記検出エリア内における前記分散相の位置座標を検出し、前記位置座標を当該分散相検出フレームのフレーム番号とともに記録する分散相位置検出部と、
前記判定部において前記分散相不検出フレームの次に前記分散相検出フレームが出現した後、別の前記分散相不検出フレームが出現するたびに、その都度、当該分散相不検出フレームおよび当該別の分散相不検出フレーム間に位置する前記分散相検出フレームのうちの２つの前記分散相検出フレームに関する前記分散相の位置座標と、前記２つの分散相検出フレーム間に位置するフレームの数と、前記入力部に入力された前記フレームレートとに基づき、前記分散相の線速度を算出して出力する分散相線速度算出部と、を備えたものであることを特徴とする請求項４に記載のシステム。
【請求項６】
前記判定部において前記分散相不検出フレームの次に前記分散相検出フレームが出現した後、別の前記分散相不検出フレームの次に別の前記分散相検出フレームが出現するたびに、その都度、当該分散相不検出フレームおよび当該別の分散相検出フレーム間に位置する前記分散相検出フレームの数と前記分散相不検出フレームの数との比を分散相・連続相流量比として算出して出力する流量比算出部を、備えたものであることを特徴とする請求項１〜請求項５のいずれかに記載のシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、マイクロ流路内を流れるエマルションが撮影された動画データに基づいて、エマルションの分散相のサイズを計測するシステムに関するものである。
【背景技術】
【０００２】
バイオテクノロジーや薬品製造等の分野では、分散相と連続相が共に液体である分散系溶液（エマルション）中の分散相（液滴）が極微小反応容器やマイクロカプセルの形成材等として広く用いられており、液滴の製造には通常マイクロ流体チップが使用され、マイクロ流体チップによれば、液滴を高速かつ大量に調製することができる（例えば、特許文献１、２参照）。
【０００３】
ところで、製造された液滴の集合体については、粒径解析を行う必要があるが、従来技術においては、この粒径解析は顕微鏡を使用して専ら手作業でなされていた。
しかしながら、対象となる液滴集合体中の液滴数は膨大であるため、手作業での粒径解析には多大な時間と労力を要していた。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００７−２９９０９号公報
【特許文献２】特開２０１５−２１１９３１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
したがって、本発明の課題は、マイクロ流体チップによって製造された液滴の粒径解析を簡単にかつ高速で行えるようにすることにある。
【課題を解決するための手段】
【０００６】
上記課題を解決するため、本発明は、マイクロ流路内を流れるエマルションが撮影された動画データに基づいて、前記エマルションの分散相のサイズを計測するシステムであって、前記動画データ、および前記マイクロ流路の断面積値の入力を受ける入力部と、予め登録された前記分散相の参照画像が格納された参照画像格納部と、前記入力部に入力された前記動画データのキャリブレーションを行うとともに、前記動画データがカラー画像データである場合は、前記動画データをグレースケール変換してモノクロ画像データにする前処理部と、前記前処理部で得られた前記動画データの各フレームに対して固定した矩形状の検出エリアを設定し、前記検出エリア内にはエッジ検出のための検出線分を設定する検出エリア設定部と、を備え、前記検出エリアは前記マイクロ流路の直線状部分の一定長をカバーし、かつ前記検出エリアの長辺が前記マイクロ流路に平行となるように配置され、前記検出線分は前記マイクロ流路の中心軸に一致するように配置され、前記検出エリアの長辺が同一フレーム内の異なる２つの前記分散相が同時に前記検出エリア内に位置し得ないような長さに設定されており、さらに、前記参照画像を用いて、前記前処理部で得られた前記動画データの各フレームに対しフレーム番号順に前記検出エリア内においてパターンマッチングを行い、前記検出エリア内に前記参照画像に一致する画像が存在する前記フレームを分散相検出フレームとする一方、前記検出エリア内に前記参照画像に一致する画像が存在しない前記フレームを分散相不検出フレームとする判定部と、前記判定部において前記分散相不検出フレームの次に前記分散相検出フレームが出現した後別の前記分散相不検出フレームが出現するたびに、その都度、当該分散相不検出フレームおよび当該別の分散相不検出フレーム間に位置する前記分散相検出フレームのうちの１つに対し前記エッジ検出を行って、関係する前記分散相の前後端縁間の距離を測定する分散相長さ計測部と、前記分散相長さ計測部による測定値と、前記入力部に入力された前記マイクロ流路の断面積値とに基づき、関係する前記分散相のサイズを算出して出力する分散相サイズ算出部と、を備えたものであることを特徴とするシステムを構成したものである。
【０００７】
上記構成において、好ましくは、前記分散相サイズ算出部は、前記分散相長さ計測部による前記測定値と、前記入力部に入力された前記マイクロ流路の断面積値とに基づき、関係する前記分散相を球状とみなしたときの粒径を算出して出力する粒径算出部を有している。
さらに好ましくは、前記システムは、前記判定部において前記分散相不検出フレームの次に前記分散相検出フレームが出現するたびに、その都度、前記分散相の個数をカウントする分散相計数部と、前記粒径算出部から出力された前記分散相の粒径の総和を、前記分散相計数部によるカウント値によって割り算することにより、前記分散相の平均粒径を算出して出力する平均粒径算出部とを備えている。
【０００８】
また、上記構成において、前記入力部が前記動画データのフレームレートの入力を受けるようになっていて、前記システムが、さらに、前記判定部において前記分散相不検出フレームの次に前記分散相検出フレームが出現した後、別の前記分散相不検出フレームの次に別の前記分散相検出フレームが出現するたびに、その都度、当該分散相不検出フレームおよび当該別の分散相検出フレーム間に位置するフレームの数を算出するフレーム数算出部と、前記フレーム数算出部による前記フレームの数の算出値と、前記入力部に入力された前記フレームレートとに基づき分散相生成速度を算出して出力する分散相生成速度算出部と、を備えていてもよい。
【０００９】
また、上記構成において、前記システムが、さらに、前記判定部において前記分散相検出フレームが出現するたびに、その都度、前記検出エリア内における前記分散相の位置座標を検出し、前記位置座標を当該分散相検出フレームのフレーム番号とともに記録する分散相位置検出部と、前記判定部において前記分散相不検出フレームの次に前記分散相検出フレームが出現した後、別の前記分散相不検出フレームが出現するたびに、その都度、当該分散相不検出フレームおよび当該別の分散相不検出フレーム間に位置する前記分散相検出フレームのうちの２つの前記分散相検出フレームに関する前記分散相の位置座標と、前記２つの分散相検出フレーム間に位置するフレームの数と、前記入力部に入力された前記フレームレートとに基づき、前記分散相の線速度を算出して出力する分散相線速度算出部と、を備えていてもよい。
ここで、分散相の線速度は連続相の線速度に等しいので、「分散相の線速度」に関する記載は、「分散相」を「連続相」に置き換えることで「連続相の線速度」にも適用し得る。以下同様。
【００１０】
また、上記構成において、前記システムが、さらに、前記判定部において前記分散相不検出フレームの次に前記分散相検出フレームが出現した後、別の前記分散相不検出フレームの次に別の前記分散相検出フレームが出現するたびに、その都度、当該分散相不検出フレームおよび当該別の分散相検出フレーム間に位置する前記分散相検出フレームの数と前記分散相不検出フレームの数との比を分散相・連続相流量比として算出して出力する流量比算出部を、備えていてもよい。
【発明の効果】
【００１１】
本発明によれば、マイクロ流路内を流れるエマルションが撮影された動画データに基づいて、エマルションの分散相のサイズを自動的に計測する。
そのため、マイクロ流体チップを用いて液滴を製造する場合、エマルションがマイクロ流体チップのマイクロ流路内を流れる様子を動画撮影し、得られた動画データを記録しさえすれば、その動画データに基づいて分散相のサイズを自動的に高速で算出することによって、製造した液滴の粒径解析を簡単かつ迅速に行うことができる。
【図面の簡単な説明】
【００１２】
【図１】本発明の１実施例による粒径解析システムの概略構成を示すブロック図である。
【図２】（Ａ）は図１のシステムのテンプレート画像格納部に格納されたテンプレート画像を示す図であり、（Ｂ）は図１のシステムの検出エリア設定部によって１つのフレームの画像に設定された検出エリアを示す図である。
【図３】（Ａ）は動画データ中の分散相検出フレームおよび分散相不検出フレームの出現シーケンスの一例を示す図であり、（Ｂ）〜（Ｄ）は、（Ａ）のうちのｉ＋１番目、ｉ＋４番目およびｉ＋１０番目のフレームに対応するフレーム画像を示す図２Ｂに類似の図である。
【図４】図１のシステムから出力された平均粒径の時間変化を表すグラフである。
【図５】図１のシステムから出力された分散相生成速度および分散相数の時間変化を表すグラフである。
【発明を実施するための形態】
【００１３】
以下、添付図面を参照しつつ、本発明の構成を好ましい実施例に基づいて説明する。
図１は、本発明による粒径解析システムの概略構成を示すブロック図である。
本発明による粒径解析システムは、マイクロ流路内を流れるエマルションが撮影された動画データに基づいてエマルションの分散相のサイズを計測するものであり、図１に示すように、当該動画データ、動画データのフレームレートおよびマイクロ流路の断面積値の入力を受ける入力部１と、予め登録された分散相の参照画像（この実施例では「テンプレート画像」と呼ぶ。図２Ａ参照）が格納されたテンプレート画像格納部２が備えられる。
なお、図２Ａに示すように、この実施例では、マイクロ流路内を流れる分散相はプラグ状を有している。
【００１４】
本発明によれば、また、入力部１に入力された動画データのキャリブレーションを行うとともに、動画データがカラー画像データである場合は、動画データをグレースケール変換してモノクロ画像データにする前処理部３が備えられる。
前処理部３による動画データのキャリブレーションにより、動画データに基づいて測定した長さ、面積および中心位置等の測定値を指定した単位で算出できるようになる。
【００１５】
また、前処理部３で得られた動画データの各フレームに対して固定した矩形状の検出エリアを設定し、検出エリア内にはエッジ検出のための検出線分を設定する検出エリア設定部４が備えられる。
【００１６】
図２Ｂは、検出エリア設定部４によって１つのフレームの画像に設定された検出エリアを示した平面図であり、図２Ｂ中、１５はマイクロ流路であり、１６は（プラグ状の）分散相であり、１７は検出エリアであり、１８は検出線分である。
図２Ｂを参照して、検出エリア１７は、マイクロ流路１５の直線状部分の一定長をカバーし、かつ検出エリア１７の長辺がマイクロ流路１５に平行となるように、かつ検出線分１８がマイクロ流路１５の中心軸に一致するように配置される。さらに、検出エリア１７の長辺が同一フレーム内の異なる２つの分散相１６が同時に検出エリア１７内に位置し得ないような長さに設定される。
【００１７】
さらに、本発明によれば、テンプレート画像（図２Ａ参照）を用いて、前処理部３で得られた動画データの各フレームに対しフレーム番号順に検出エリア１７内においてパターンマッチング（テンプレートマッチング）を行い、検出エリア１７内にテンプレート画像に一致する画像が存在するフレームを分散相検出フレームとする一方、検出エリア１７内にテンプレート画像に一致する画像が存在しないフレームを分散相不検出フレームとする判定部５が備えられる。
【００１８】
さらに、判定部５において分散相不検出フレームの次に分散相検出フレームが出現するたびに、その都度、分散相１６の個数をカウントするとともに、カウント値を出力する分散相計数部６と、判定部５において分散相不検出フレームの次に分散相検出フレームが出現した後別の分散相不検出フレームが出現するたびに、その都度、当該分散相不検出フレームおよび当該別の分散相不検出フレーム間に位置する分散相検出フレームのうちの１つに対しエッジ検出を行うことにより、関係する分散相１６の前後端縁間の距離を測定する分散相長さ計測部７が備えられる。
【００１９】
以下、分散相計数部６および分散相長さ計測部７のこの動作を、図面を参照して具体的に説明する。
図３Ａは、動画データ中の分散相検出フレームおよび分散相不検出フレームの出現シーケンスの一例を示す図であり、図３Ａ中、各正方形枠はフレームを表し、各正方形枠の下側の記号はフレーム番号を表し、○印はフレームが分散相検出フレームであることを表し、×印はフレームが分散相不検出フレームであることを表している。
また、図３Ｂ〜図３Ｄは、図３Ａのうちの（ｉ＋１）番目、（ｉ＋４）番目および（ｉ＋１０）番目のフレームに対応するフレーム画像を示す図２Ｂに類似の図である。
図３Ｂでは、１つの分散相１６ａが検出エリア１７内において初めて検出されており、図３Ｃでは、当該分散相１６ａが検出エリア１７内において最後に検出され、図３Ｄでは、次の分散相１６ｂが検出エリア内において初めて検出されている。
【００２０】
図３においては、分散相不検出フレームであるフレーム（ｉ）の次に分散相検出フレームであるフレーム（ｉ＋１）が出現し、その後、分散相不検出フレームであるフレーム（ｉ＋９）の次に分散相検出フレームであるフレーム（ｉ＋１０）が出現する。
よって、分散相計数部６は、フレーム（ｉ＋１）が出現したとき分散相数１をカウントし、次に、フレーム（ｉ＋１０）が出現したき分散相数１をカウントする。
【００２１】
図３においては、分散相不検出フレームであるフレーム（ｉ）の次に分散相検出フレームであるフレーム（ｉ＋１）が出現した後、別の分散相不検出フレームであるフレーム（ｉ＋５）が出現する。
よって、分散相長さ計測部７は、フレーム（ｉ）とフレーム（ｉ＋５）の間に位置するフレーム（分散相検出フレーム）のうちの１つに対しエッジ検出を行い、分散相１６の前後端縁間の距離、すなわち分散相の長さを測定する。
【００２２】
また、本発明によれば、分散相長さ計測部７による測定値と、入力部１に入力されたマイクロ流路の断面積値とに基づき、分散相１６のサイズを算出して出力する分散相サイズ算出部が備えられる。
分散相サイズ算出部は、この実施例では、分散相長さ計測部７による測定値と、入力部１に入力されたマイクロ流路の断面積値とに基づき、関係する分散相１６を球状とみなしたときの粒径を算出して出力する粒径算出部８を有している。
【００２３】
さらに、本発明によれば、粒径算出部８から出力された分散相１６の粒径の総和を、分散相計数部６によるカウント値によって割り算することにより、分散相１６の平均粒径を算出して出力する平均粒径算出部９が備えられる。
【００２４】
粒径算出部８による分散相１６の粒径の算出および平均粒径算出部９による平均粒径の算出は、次のようにしてなされる。
本発明では、分散相１６は、マイクロ流路１５内では空間的制限を受けてプラグ状または球状またはそれらの中間の形状を有するが、マイクロ流体チップの液溜や別の容器に捕集された状態では、空間的制限を受けずに球状の液滴となる点に鑑み、エマルションの分散相のサイズの計測が、分散相の長さの測定値とマイクロ流路の断面積値とに基づいて分散相を球状とみなしたときの粒径を算出することによってなされる。
【００２５】
今、分散相１６が空間的制限を受けない状態で直径（粒径）Ｄの完全な球であり、分散相１６がマイクロ流路１５内で（空間的制限を受けた状態で）断面積Ｓの完全な直方体であるとし、分散相長さ計測部７による測定値をＬとすれば、
Ｓ×Ｌ＝（４／３）×π（Ｄ／２）３                  （１）
が成立する。（１）式をＤについて解けば
Ｄ＝（６×Ｓ×Ｌ／π）１／３                         （２）
となる。
【００２６】
次に、分散相１６の粒径の総和をＤｔｏｔａｌとし、分散相計数部６のカウント値をＮとすれば、平均粒径Ｄａｖは、
Ｄａｖ＝Ｄｔｏｔａｌ／Ｎ                                （３）
によって求められる。
そして、粒径算出部８は（２）式に基づいて粒径を算出し、平均粒径算出部９は（３）式に基づいて平均粒径を算出する。
【００２７】
図４には、平均粒径算出部９から出力された、平均粒径の時間変化を表すグラフの一例を示した。図４のグラフ中、縦軸は平均粒径（μｍ）を表し、横軸は時間（ｍｉｎ）を表わす。
【００２８】
本発明によれば、また、判定部５において分散相不検出フレームの次に分散相検出フレームが出現した後、別の前記分散相不検出フレームの次に別の前記分散相検出フレームが出現するたびに、その都度、当該分散相不検出フレームおよび当該別の分散相検出フレーム間に位置するフレームの数を算出するフレーム数算出部１０と、フレーム数算出部１０によるフレームの数の算出値と、入力部１に入力されたフレームレートとに基づき分散相生成速度を算出して出力する分散相生成速度算出部１１が備えられる。
【００２９】
フレーム数算出部１０および分散相生成速度算出部１１の動作を図３に基づいて説明すると、図３においては、分散相不検出フレームであるフレーム（ｉ）の次に分散相検出フレームであるフレーム（ｉ＋１）が出現した後、別の分散相不検出フレームであるフレーム（ｉ＋９）の次に別の分散相検出フレームであるフレーム（ｉ＋１０）が出現する。
【００３０】
よって、フレーム数算出部１０は、フレーム（ｉ）およびフレーム（ｉ＋１０）間に位置するフレーム数
（ｉ＋９）−ｉ＝９
を算出する。
そして、分散相生成速度算出部１１は、フレームレートをαとして
α／９
によって、分散相生成速度を算出する。
【００３１】
図５には、分散相生成速度算出部１１から出力された、分散相生成速度の時間変化を表すグラフ（点線で描いた）を示した。図５のグラフ中、縦軸（左側目盛）は分散相生成速度（Ｈｚ）を表し、横軸は時間（ｍｉｎ）を表している。
なお、図５には、分散相数の時間変化のグラフ（実線で描いた）も同時に示した（縦軸右側目盛参照）。
【００３２】
また、本発明によれば、判定部５において分散相検出フレームが出現するたびに、その都度、検出エリア１７内における分散相１６の位置座標を検出し、検出した位置座標を当該分散相検出フレームのフレーム番号とともに記録する分散相位置検出部１２と、判定部５において分散相不検出フレームの次に分散相検出フレームが出現した後、別の分散相不検出フレームが出現するたびに、その都度、当該分散相不検出フレームおよび当該別の分散相不検出フレーム間に位置する分散相検出フレームのうちの２つの分散相検出フレームに関する分散相１６の位置座標と、２つの分散相検出フレーム間に位置するフレームの数と、フレームレートとに基づき、分散相１６の線速度を算出して出力する分散相線速度算出部１３が備えられる。
【００３３】
この場合、分散相位置検出部１２は、分散相検出フレームが出現するたびに、検出エリア１７の左上隅角を原点として検出エリア１７内に設定したｘｙ座標系（ｘ軸が検出エリアの長辺に平行）に基づき、分散相１６の中心の位置座標を検出し、検出した位置座標を該当するフレーム番号とともに記録する。
【００３４】
分散相線速度算出部１３による分散相線速度の算出法を図３に基づいて説明すると、図３においては、分散相不検出フレームであるフレーム（ｉ）の次に分散相検出フレームであるフレーム（ｉ＋１）が出現した後、別の分散相不検出フレームであるフレーム（ｉ＋５）が出現する。
よって、分散相線速度算出部１３は、フレーム（ｉ）およびフレーム（ｉ＋５）間に位置するフレーム（分散相検出フレーム）のうちの２つのフレーム、例えば、フレーム（ｉ＋１）とフレーム（ｉ＋４）に関する分散相１６のｘ座標（ｘ座標）を分散相位置検出部１２から抽出する。
【００３５】
今、フレーム（ｉ＋１）に関する分散相１６のｘ座標をｘ１とし、フレーム（ｉ＋４）の分散相１６のｘ座標をｘ２とする。
次いで、フレーム（ｉ＋１）およびフレーム（ｉ＋４）間のフレーム数Ｆを、
Ｆ＝（ｉ＋４）−ｉ＝４
によって求め、さらに、フレームレートをαとして
Ｆ／α＝４／α
によって、フレーム（ｉ＋１）およびフレーム（ｉ＋４）間の経過時間を算出する。
そして、
（ｘ２−ｘ１）／（４／α）
によって、分散相線速度を算出する。
なお、分散相の線速度は連続相の線速度に等しいので、「分散相の線速度」に関する記載は、「分散相」を「連続相」に置き換えることで「連続相の線速度」にも適用し得る。
【００３６】
本発明によれば、さらに、判定部５において分散相不検出フレームの次に分散相検出フレームが出現した後、別の分散相不検出フレームの次に別の分散相検出フレームが出現するたびに、その都度、当該分散相不検出フレームおよび当該別の分散相検出フレーム間に位置する分散相検出フレームの数と分散相不検出フレームの数との比を、分散相・連続相流量比として算出して出力する流量比算出部１４が備えられる。
【００３７】
流量比算出部１４の動作を図３に基づいて説明すると、図３においては、分散相不検出フレームであるフレーム（ｉ）の次に分散相検出フレームであるフレーム（ｉ＋１）が出現した後、別の分散相不検出フレームであるフレーム（ｉ＋９）の次に別の分散相検出フレームであるフレーム（ｉ＋１０）が出現する。
【００３８】
よって、流量比算出部１４は、フレーム（ｉ）およびフレーム（ｉ＋１０）の間に位置する分散相検出フレームの数
（ｉ＋４）−ｉ＝４
と、分散相不検出フレームの数
（ｉ＋９）−ｉ−４＝５
との比４／５を算出して出力する。
【００３９】
本発明によれば、マイクロ流体チップを用いて液滴を製造する場合、エマルションがマイクロ流体チップのマイクロ流路内を流れる様子を動画撮影し、得られた動画データを記録しさえすれば、その動画データに基づいて分散相の個数、粒径および平均粒径を自動的に高速で算出することによって、製造した液滴の粒径解析が簡単かつ迅速に行える。
加えて、マイクロ流体チップによる液滴製造の間におけるマイクロ流路内を流れる分散相の生成速度、分散相の線速度、分散相および連続相の流量比も自動的に算出することができる。
【００４０】
以上、本発明の１実施例について説明したが、本発明の構成は上記実施例に限定されず、特許請求の範囲に記載した構成の範囲内で種々の変形例を案出できることは言うまでもない。
例えば、上記実施例では、システムは、分散相の個数、粒径および平均粒径に加えて、分散相生成速度、分散相線速度および分散相・連続相流量比を算出するようになっているが、分散相の平均粒径、分散相生成速度、分散相線速度および分散相・連続相流量比を算出するための構成については、必要に応じてシステムに備えられればよい。
また、例えば、システムが分散相のサイズのみを算出し出力するような構成とすることもできる。
【符号の説明】
【００４１】
１  入力部
２  テンプレート画像格納部
３  前処理部
４  検出エリア設定部
５  判定部
６  分散相計数部
７  分散相長さ計測部
８  粒径算出部
９  平均粒径算出部
１０  フレーム数算出部
１１  分散相生成速度算出部
１２  分散相位置検出部
１３  分散相線速度算出部
１４  流量比算出部
１５  マイクロ流路
１６、１６ａ、１６ｂ  分散相
１７  検出エリア
１８  検出線分
【図１】
2018169224.tif 000003
【図２】
2018169224.tif 000004
【図３】
2018169224.tif 000005
【図４】
2018169224.tif 000006
【図５】
2018169224.tif 000007