JP7427291B2 - イメージングフローサイトメーター - Google Patents

Description

本発明は、イメージングフローサイトメーターに関する。
本願は、２０２０年９月１０日に、日本に出願された特願２０２０－１５２３３１号に基づき優先権を主張し、その内容をここに援用する。

従来、流路を流れる観測対象物に照明光を照射して発せられる光を測定することによって観測対象物の特徴を評価するフローサイトメトリー法や、このフローサイトメトリー法を用いたフローサイトメーターが知られている（例えば、特許文献１）。また、観測対象物となる細胞・細菌等の微粒子を画像によって評価する蛍光顕微鏡やイメージングサイトメーターが知られている。
加えて、フローサイトメーターの一種で、従来のフローサイトメーターと同等のスループットで微粒子の形態情報を高速に撮像するイメージングフローサイトメーターが知られている（例えば、特許文献２）。イメージングフローサイトメーターでは、流路中に照明光の当たる計測野が設けられ、細胞が計測野を通り抜ける。つまりイメージングフローサイトメーターでは、スキャンされる際に、二次元または三次元画像の撮像、または画像に相当する細胞空間情報の計測を行う。
また、観測対象物となる細胞に帯状の励起光を照射して、この励起光が照射される照射面に対して異なる角度の面の画像や、観察対象物の３次元画像を得る多面光学顕微鏡が知られている（例えば、特許文献３及び特許文献４）。特許文献４に記載の多面光学顕微鏡では、Ａｘｉａｌ ｐｌａｎｅ ｏｐｔｉｃａｌ ｍｉｃｒｏｓｃｏｐｙ（ＡＰＯＭ）の技術が用いられている（例えば、非特許文献１）。ＡＰＯＭとは、走査を行うことなく観測対象物の対物レンズの光軸に対して平行な断面を直接撮像する多面光学顕微鏡の技術である。

特許第５５３４２１４号公報 米国特許第６２４９３４１号明細書 米国特許出願公開第２０１５／０１９２７６７号明細書 米国特許第８５８２２０３号明細書

Ｔｏｎｇｃａｎｇ Ｌｉ、Ｓａｄａｏ Ｏｔａ、Ｊｅｏｎｇｍｉｎ Ｋｉｍ、Ｚｉ Ｊｉｎｇ Ｗｏｎｇ、Ｙｕａｎ Ｗａｎｇ 、Ｘｉａｏｂｏ Ｙｉｎ、Ｘｉａｎｇ Ｚｈａｎｇ、「Ａｘｉａｌ ｐｌａｎｅ ｏｐｔｉｃａｌ ｍｉｃｒｏｓｃｏｐｙ」、「Ｓｃｉｅｎｔｉｆｉｃ Ｒｅｐｏｒｔｓ」、Ｓｐｒｉｎｇｅｒ Ｎａｔｕｒｅ、２０１４年１２月１日、４巻、７２５３号 Ｄａｎｕｔａ Ｍ． Ｂｕｋｏｗｓｋａ、Ｌａｄｉｓｌａｖ Ｄｅｒｚｓｉ、Ｓｚｙｍｏｎ Ｔａｍｂｏｒｓｋｉ、Ｍａｃｉｅｊ Ｓｚｋｕｌｍｏｗｓｋｉ、Ｐｉｏｔｒ Ｇａｒｓｔｅｃｋｉ、Ｍａｃｉｅｊ Ｗｏｊｔｋｏｗｓｋｉ、「Ａｓｓｅｓｓｍｅｎｔ ｏｆ ｔｈｅ ｆｌｏｗ ｖｅｌｏｃｉｔｙ ｏｆ ｂｌｏｏｄ ｃｅｌｌｓ ｉｎ ａ ｍｉｃｒｏｆｌｕｉｄｉｃ ｄｅｖｉｃｅ ｕｓｉｎｇ ｊｏｉｎｔ ｓｐｅｃｔｒａｌ ａｎｄ ｔｉｍｅ ｄｏｍａｉｎ ｏｐｔｉｃａｌ ｃｏｈｅｒｅｎｃｅ ｔｏｍｏｇｒａｐｈｙ」、「Ｏｐｔｉｃｓ Ｅｘｐｒｅｓｓ」、ＯＳＡ Ｐｕｂｌｉｓｈｉｎｇ、２０１３年１０月１日、２１巻、２０号、ｐ．２４０２５－２４０３８

上述したようなイメージングサイトメーターでは、流路を流れる流体の流速は、流路の幅方向、あるいは高さ方向についてばらつきが生じ得る。また、何もしなければ、観測対象物の位置は流路の幅方向、あるいは高さ方向についてばらついてしまう。その結果、流路に生じた流速のばらつきによって、流路を流れる観測対象物の速度（進行方向と速さ）のばらつき、あるいは回転が起こってしまう。その結果、撮像位置においても測対象物の流線方向に垂直な方向についての位置あるいは速度が安定しない、あるいは対象物が回転してしまうことがあった。そのため、観測対象物を流線方向に走査する場合に撮像条件が安定しない場合があった。撮像条件には、撮像過程における観測対象物の流線方向に垂直な方向における位置、回転、及び速度（進行方向と速さ）が含まれる。

本発明は上記の点に鑑みてなされたものであり、流路を流れる観測対象物を流線方向に走査する場合に撮像条件を安定させることができるイメージングフローサイトメーターを提供する。

本発明は上記の課題を解決するためになされたものであり、本発明の一態様は、観測対象物が流される流路であって、幅方向の長さが高さ方向の長さよりも長い流路と、前記流路に対して音波を定在波にして照射する音響素子と、前記流路に対して照明光を照射する光源と、前記照明光が照射される位置を通過した前記観測対象物を計測または撮像することにより、前記観測対象物が前記流路を流れる方向である流線方向についての前記観測対象物の断面に含まれる線を少なくとも撮像する撮像部と、前記撮像部が前記線を時系列に撮像して得られる複数の撮像画像に基づいて、前記観測対象物が前記流線方向に走査された画像を生成する走査画像生成部と、を備え、前記流路を流れる前記観測対象物の流速は前記撮像部に応じて決まっており、前記音響素子と前記撮像部との前記流路の長さ方向についての距離は、前記音響素子による前記音波の音響効果による複数の前記観測対象物の整列が持続する程度に短く、前記高さ方向について前記定在波の節の数は所定の個数であり、前記撮像部は、前記撮像部に備えられる撮像素子の画素列のうち前記高さ方向の前記定在波の節の位置に対応する画素列のみを用いて撮像を行うイメージングフローサイトメーターである。

また、本発明の一態様は、上記のイメージングフローサイトメーターにおいて、前記高さ方向について前記定在波の節の数は１個であり、前記幅方向について前記定在波の節の数は０個である。

また、本発明の一態様は、上記のイメージングフローサイトメーターにおいて、前記高さ方向について前記定在波の節の数は１個であり、前記幅方向について前記定在波の節の数は１個以上である。

また、本発明の一態様は、上記のイメージングフローサイトメーターにおいて、前記流路には、前記幅方向に整列された複数の前記観測対象物が流される。

本発明によれば、流路を流れる観測対象物を流線方向に走査する場合に撮像条件を安定させることができる。

本発明の第１の実施形態に係る細胞測定システムの構成の一例を示す図である。 本発明の第１の実施形態に係る細胞測定システムの機能構成の一例を示す図である。 本発明の第１の実施形態に係る流路の構成の一例を示す図である。 本発明の第１の実施形態に係る流路における幅方向の流速分布の一例を示す図である。 本発明の第１の実施形態に係る流路における高さ方向の流速分布の一例を示す図である。 本発明の第１の実施形態に係る音波によって流路を流れる複数の細胞が高さ方向について整列される様子の一例を示す図である。 本発明の第１の実施形態に係る音響素子が音波を流路に照射していない場合に複数の細胞が流路を流れる様子の一例を示す図である。 本発明の第１の実施形態に係る音波の定在波の節の数が幅方向について１以上である場合に複数の細胞が流路を流れる様子の一例を示す図である。 本発明の第１の実施形態に係る撮像素子が撮像する細胞の断面の一例を示す図である。 本発明の第１の実施形態に係る撮像素子において像が結像される領域の画素列の一例を示す図である。 本発明の第１の実施形態に係る３次元画像生成部が断面画像を結合する順番の一例を示す図である。 本発明の第１の実施形態に係るイメージングフローサイトメーターの動作の一例を示す流れ図である。 本発明の第１の実施形態に係る流路の高さ方向から複数の細胞が撮像された明視野画像の一例を示す図である。 本発明の第１の実施形態に係る高さ方向から所定の角度だけ傾いた面において複数の細胞が撮像された蛍光画像の一例を示す図である。 本発明の第１の実施形態に係る高さ方向から所定の角度だけ傾いた面が異なる時刻において撮像された断面画像の一例を示す図である。 本発明の第１の実施形態に係る３次元画像の一例を示す図である。 本発明の第１の実施形態に係るｘ方向とｙ方向についての断面の画像の一例を示す図である。 本発明の第１の実施形態に係る断面画像の各領域の拡大図の一例を示す図である。 本発明の第２の実施形態に係る細胞測定システムの機能構成の一例を示す図である。 本発明の第２の実施形態に係る２次元画像生成部が線画像を結合する順番の一例を示す図である。 本発明の実施例に係る細胞の３次元画像を２次元平面（ｘｙ平面）に射影した画像の一例を示す図である。 本発明の実施例に係る図２１に示した画像の領域を拡大した画像の一例を示す図である。 本発明の実施例に係るＤＡＰＩによって染色された細胞の細胞分裂の中期の３次元画像を２次元平面に射影した画像の一例を示す図である。 本発明の実施例に係るＤＡＰＩによって染色された細胞の細胞分裂の後期の３次元画像を２次元平面に射影した画像の一例を示す図である。 本発明の実施例に係る３次元画像のうちＤＡＰＩによって染色された領域の一例を示す図である。 本発明の実施例に係るｘｙ平面、ｙｚ平面、及びｘｚ平面それぞれに射影した場合の２次元画像の一例を示す図である。 本発明の実施例に係る座標軸を主成分分析によって得られた座標軸へと回転させた場合の２次元画像の一例を示す図である。 本発明の実施例に係る広い範囲の３次元画像の一例を示す図である。 本発明の実施例に係るＭＰＭ－２及びＤＡＰＩの強度を解析した結果を示す散布図の一例を示す図である。 本発明の実施例に係る２つのＤＡＰＩ染色された領域の３次元画像がｘｙ平面に射影された２次元画像の一例を示す図である。 本発明の実施例に係る１つのＤＡＰＩ染色された領域の３次元画像が平面に射影された２次元画像の一例を示す図である。 本発明の実施例に係る３次元イメージングによる画像と２次元イメージングに相当する画像とでＤＡＰＩ領域のアスペクト比を比較した結果の一例を示す図である。

（第１の実施形態）
以下、図面を参照しながら本発明の実施形態について詳しく説明する。図１は、本実施形態に係る細胞測定システム１の構成の一例を示す図である。

細胞測定システム１は、イメージングフローサイトメーター２０と、解析装置３０と、表示部１０とを備える。イメージングフローサイトメーター２０は、観測対象物が流される流路を備える。なお、本実施形態では、観測対象物が細胞である場合について説明するが、観測対象物は細胞に限られない。観測対象物は、細菌等の微粒子であってもよい。
イメージングフローサイトメーター２０が観測対象物を観察するために用いる顕微鏡には、蛍光顕微鏡、あるいは蛍光染色を必要としない各種の顕微鏡が用いられる。イメージングフローサイトメーター２０に蛍光顕微鏡が用いられる場合、観測対象物は蛍光染色される。イメージングフローサイトメーター２０に蛍光染色を必要としない各種の顕微鏡が用いられる場合、観測対象物は蛍光染色されない。蛍光染色を必要としない各種の顕微鏡とは、例えば、明視野顕微鏡、暗視野顕微鏡、位相差顕微鏡、位相顕微鏡、光散乱顕微鏡、ラマン散乱光顕微鏡、吸収位相顕微鏡などである。蛍光染色を必要としない顕微鏡を用いるイメージングフローサイトメーターでは、蛍光染色されていない観測対象物の画像を撮像するラベルフリーイメージングの技術が用いられ、細胞、細菌等を非侵襲的に観察することができる。イメージングフローサイトメーター２０が観測対象物を観察するイメージングは、上述したものに限られず、可視光、可視光以外の電磁波、熱、あるいは振動波などによる各種のイメージングが用いられてよい。

イメージングフローサイトメーター２０は、この流路を流される細胞の３次元画像を生成する。解析装置３０は、イメージングフローサイトメーター２０が生成した３次元画像を解析する。解析装置３０は、３次元画像を解析することによって、例えば、当該３次元画像として再構築された細胞の分類、あるいは当該細胞についての各種の情報の読み取りを実行する。解析装置３０は、例えば、パーソナルコンピュータ（ＰＣ：Ｐｅｒｓｏｎａｌ Ｃｏｍｐｕｔｅｒ）である。解析装置３０は、非破壊の細胞解析に好適に用いられる。
表示部１０は、イメージングフローサイトメーター２０が生成した３次元画像、あるいは解析装置３０による解析結果を表示する。表示部１０は、例えば、液晶ディスプレイを備えており、さまざまな画像を表示する。この表示部１０に表示される画像には、イメージングフローサイトメーター２０が生成する細胞の３次元画像が含まれている。

［イメージングフローサイトメーター２０の機能構成］
次に、図２を参照して、イメージングフローサイトメーター２０の機能構成について説明する。
図２は、本実施形態の細胞測定システム１の機能構成の一例を示す図である。図２には、３次元直交座標系として、ｘｙｚ座標系を示す。本実施形態において、ｘ軸方向は、流路２１の幅方向である。また、ｙ軸方向が流路２１における細胞ＣＬが流される方向である。細胞ＣＬは、ｙ軸方向の＋ｙ方向に流される。ｚ軸方向は、流路２１と直交する方向であって、流路２１の高さ方向である。

イメージングフローサイトメーター２０は、流路２１と、音響素子２２と、撮像部２３と、記憶部２５と、３次元画像生成部２６とを備える。

ここで図３から図８を参照し、流路２１、及び音響素子２２について説明する。図３は、本実施形態に係る流路２１の構成の一例を示す図である。流路２１は、観測対象物が流される。観測対象物は、流路２１を流れる流体と共に流れされる。本実施形態では、流路２１には、流路２１の幅方向に整列された複数の観測対象物が流される。

複数の観測対象物が幅方向に整列されるとは、高さ方向の位置が複数の観測対象物相互間において略等しくされて、複数の観測対象物が幅方向に１列に配置されることである。隣り合う観測対象物同士の間隔は任意であってよい。複数の観測対象物が幅方向に整列されている場合、流路２１の長さ方向についての位置は、複数の観測対象物相互間において揃っていることが好ましい。複数の観測対象物が幅方向に整列されている場合、流路２１の長さ方向についての位置は、複数の観測対象物について、少なくとも、ある観測対象物が他の観測対象物と長さ方向について重なり合う部分を有していることが好ましい。図３に示す例では、流路２１には、観測対象物として複数の細胞ＣＬｓが、流路２１の幅方向に整列されて流されている。
なお、流路２１には、１つの観測対象物が流されてもよい。

流路２１では、幅方向長さＦＷは高さ方向長さＦＨよりも長い。ここで、流路の幅方向の長さの高さ方向の長さに対する比率が非常に大きい場合、流路内における幅方向の流速分布は、幅方向の位置によらず一定となる（例えば、非特許文献２を参照）。一方、この場合、高さ方向の流速分布は、高さ方向の位置によって大きく変化する。

［流速分布について］
図４に、流路２１における幅方向の流速分布の一例を示す。流路２１では、幅方向長さＦＷは高さ方向長さＦＨよりも長いために、幅方向の流速分布は、流路２１の壁の付近を除いて、幅方向（ｘ軸方向）の位置によらずほぼ一定となる。

図５に、流路２１における高さ方向の流速分布の一例を示す。流路２１では、高さ方向の流速分布は、高さ方向（ｚ軸方向）の位置によって大きく変化する。ここで高さ方向の流速は、流路２１の壁から遠いほど大きくなる。一方、高さ方向の流速は、流路２１の壁に近いほど小さくなる。したがって、高さ方向の流速は、高さ方向（ｚ軸方向）の位置について中央の位置で最大となる。

このように流路２１では、高さ方向について流速のばらつきが生じ得る。ここで高さ方向の流速分布の結果としてもあり得るが、他の色々な原因（例えば、重力）によっても、流路２１を流れる複数の細胞ＣＬｓの高さ方向についての位置は元から（最初から）ばらついている。その結果、高さ方向の流速分布によって、流路２１を流れる複数の細胞ＣＬｓはそれぞれ、速度（進行方向と速さ）が変化し得る、あるいは回転し得る。それらのため、複数の細胞ＣＬｓそれぞれについて流線方向に垂直な方向における位置、あるいは向きが、複数の細胞ＣＬｓを撮像する過程において撮像部２３に対して変化してしまう。
また、複数の細胞ＣＬｓそれぞれについて流線方向に垂直な方向における位置がばらついた場合、複数の細胞ＣＬｓのうち撮像部２３に備えられる撮像素子２８によって撮像される位置以外の位置を通過してしまうものがある場合がある。その場合、１フレームあたりの細胞の数が少なくなってしまう。

イメージングフローサイトメーター２０では、音響素子２２による音波ＡＷを用いて観測対象物の流路２１の高さ方向の位置を揃える。音響素子２２は、一例として、トランスデューサを含んで構成される。

図６に、本実施形態に係る音波ＡＷによって流路２１を流れる複数の細胞ＣＬｓが高さ方向について整列される様子の一例を示す。音響素子２２は、流路２１に対して音波ＡＷを定在波にして照射する。音波ＡＷの波長は、流路２１の高さに応じて決まる。
複数の細胞ＣＬｓは、音響効果によって音波の定在波の節に集まるように整列される。音響効果によって音波の定在波の節に観測対象物を集中させる方法を音響フォーカシングという。音響素子２２から照射される音波ＡＷは、例えば、超音波である。なお、音響素子２２から照射される音波ＡＷは、超音波より低い周波数の音波であってもよい。

本実施形態では、音波ＡＷの定在波の節の数は、一例として、高さ方向について１個である。音波ＡＷの定在波の節の位置は、流路２１の高さ方向について高さの中央の位置である。音波ＡＷの波長は、流路２１の高さに等しい。
複数の細胞ＣＬｓは、音波ＡＷの定在波の節の数が高さ方向について１個であることに応じて、高さ方向について高さの中央の位置に１列に整列されている。図６では、一例として、複数の細胞ＣＬｓのうち細胞ＣＬ１の付近の位置に生じている音波ＡＷの振幅のパターンを示している。音波ＡＷの定在波は、流路２１の幅方向の他の位置についても、同様にして流路２１の高さの中央の位置に節をもつ。
一方、音波ＡＷの定在波の節の数は、一例として、幅方向について０個である。換言すれば、音波ＡＷは幅方向については定在波を形成しない。

比較のために、図７に、音響素子２２が音波ＡＷを流路２１に照射していない場合に複数の細胞ＣＬｓが流路２１を流れる様子の一例を示す。図７では、流路２１の長さ方向（ｙ軸方向）についてのある時刻ｔ１におけるある断面が示されている。複数の細胞ＣＬｓ－１は、複数の細胞ＣＬｓのうち流路２１の幅方向（ｘ軸方向）について中央付近を流れる速度の速い細胞である。細胞ＣＬｓ－２、ＣＬｓ－３は、複数の細胞ＣＬｓのうち流路２１の幅方向（ｘ軸方向）について中央から離れた位置を流れる速度の遅い細胞である。速度の遅い細胞ＣＬｓ－２、ＣＬｓ－３は、時刻ｔ１において、速度の速い細胞ＣＬｓ－１よりも上流を流れており、図７に示す断面の位置には到達していない。
また、複数の細胞ＣＬｓでは、流路２１の高さ方向（ｚ軸方向）についての位置が、高さ方向の流速分布のために揃っていない。

また、上述したように、本実施形態では、音波ＡＷの定在波の節の数は、幅方向について０個である場合の一例について説明するが、これに限られない。音波ＡＷの定在波の節の数が幅方向について１個以上であってもよい。図８に、音波ＡＷの定在波の節の数が幅方向について１以上である場合に複数の細胞ＣＬｓが流路２１を流れる様子の一例を示す。図８では、一例として、音波ＡＷの定在波の節の数は幅方向について５個である。複数の細胞ＣＬｓは、音響効果によって、幅方向について、定在波の幅方向の節の位置に整列する。

幅方向について定在波の節がある場合、流路２１を流れる複数の細胞ＣＬｓは、幅方向について当該節の位置を流れる。そのため、幅方向に並べられる複数の細胞ＣＬｓの数は、幅方向についての定在波の節の数と等しくなる。つまり、幅方向について定在波の節は、幅方向に並べられる複数の細胞ＣＬｓの数を制限してしまう。

なお、幅方向について定在波の節は、幅方向に並べられる複数の細胞ＣＬｓの数を制限する一方で、幅方向に並べられた複数の細胞ＣＬｓの隣り合う細胞同士の距離を所定の距離以上とすることができる。そのため、幅方向について定在波の節は、複数の細胞ＣＬｓ同士が相互作用してしまうことを抑制できる。

一方、本実施形態では、音波ＡＷの定在波の節の数は幅方向について０個であるため、幅方向に並べられる複数の細胞ＣＬｓの数は、当該節の数によっては制限されない。本実施形態では、幅方向に並べられる複数の細胞ＣＬｓの数は、それら複数の細胞ＣＬｓの大きさと流路２１の幅とによって制限される。

図２に戻り、音響素子２２が備えられる位置について説明する。
音響素子２２は、流路２１の長さ方向について撮像部２３が備えられる位置よりも上流側（－ｙ軸方向）であって、当該位置から所定の距離だけ離れた位置に備えられる。したがって、音響素子２２が照射する音波ＡＷによる定在波は、撮像部２３が備えられる位置よりも上流側（－ｙ軸方向）の位置において生じる。

音波ＡＷによって高さ方向について整列された複数の細胞ＣＬｓは、音波ＡＷが照射される位置を通過すると、下流側においてやがては流速分布の影響を受ける。つまり、一度整列された複数の細胞ＣＬｓには、音波ＡＷが照射される位置の下流側において、位置のばらつき、あるいは回転が生じ得る。イメージングフローサイトメーター２０では、撮像部２３は、一度整列された複数の細胞ＣＬｓに位置のばらつき、あるいは回転が生じる位置よりも上流側の位置に配置される。したがって、撮像部２３は、複数の細胞ＣＬｓが高さ方向について整列された状態を撮像する。
このように、音響素子２２と撮像部２３との流路２１の長さ方向についての距離は、音響素子２２による音波ＡＷの音響効果による複数の観測対象物の整列が持続する程度に短い。

なお、音響素子２２は、流路２１の長さ方向について撮像部２３が備えられる位置よりも上流側（－ｙ軸方向）であって、当該位置から所定の距離だけ離れた位置に備えられるため、音響素子２２の外観が撮像部２３が撮像する画像に映り込んでしまうことはない。

撮像部２３は、光源２７と、対物レンズＯＧと、撮像素子２８と、第１光学素子Ｌ１と、第２光学素子Ｌ２と、第３光学素子Ｌ３と、制御部２００とを備える。

光源２７は、流路２１に対して帯状の照明光ＬＳを照射する。具体的には、光源２７は、流路２１に対して帯状のコヒーレント光を照射する。帯状の照明光ＬＳとは、コヒーレント光を絞ることにより、帯状にした照明光である。帯状にした照明光とは、ｘ軸方向に試料の径よりも長い幅をもち、ｙ軸方向に試料の径よりも短い厚みに調整されたコヒーレント光である。本実施形態では、光源２７は、ダイクロイックミラーＭ１と、対物レンズＯＧとを介して、流路２１に対して帯状の照明光ＬＳを照射する。なお、照明光ＬＳは、コヒーレント光であっても、インコヒーレント光であってもよい。本実施形態では、光源２７が照射する照明光ＬＳは、一例として、コヒーレント光である。
光源２７は、ダイクロイックミラーＭ１と、対物レンズＯＧとを介して、流路２１に対して照明光ＬＳを照射する。

以下の説明では、流路２１の照明光ＬＳが照射される位置のことを、照射位置とも記載する。この照射位置を通過した細胞ＣＬは、照明光ＬＳが照射されると光Ｌを発する。例えば、細胞ＣＬが蛍光染色されている場合、照明光ＬＳとして励起光が用いられ、細胞ＣＬは、励起光である照明光ＬＳによって蛍光分子が励起されることにより発光する。この発光による光は蛍光であり、光Ｌとはこの蛍光である。

撮像部２３は、照明光ＬＳが照射される位置を通過した細胞ＣＬを撮像する。撮像部２３は、一例として、照明光ＬＳが照射される位置を通過した細胞ＣＬからの光Ｌを撮像する。これにより、撮像部２３は、細胞ＣＬのある断面を撮像する。ここで、ある断面とは、照明光ＬＳによって細胞ＣＬが光Ｌを発した面である。言い換えると、ある断面とは、細胞ＣＬが流路２１を流れる方向と交差する面である。なお、この細胞のある断面を撮像する光学系の構成は、例えば、米国特許出願公開第２０１５／０１９２７６７号明細書及び米国特許第８５８２２０３号明細書などに開示される構成である。またなお、撮像部２３は、細胞ＣＬの空間情報を計測してもよい。またなお、イメージングフローサイトメーター２０が観測対象物を観察するイメージングの種類に応じて、撮像部２３は、可視光、可視光以外の電磁波、熱、あるいは振動波などを計測または撮像に用いてよい。

対物レンズＯＧは、細胞ＣＬからの光Ｌを集光する。対物レンズＯＧとは、照射位置に焦点があう位置に配置される対物レンズである。以下の説明では、照射位置に焦点が合う位置のことを焦点位置ＦＰとも記載する。なお、対物レンズＯＧは、ドライ対物レンズあっても、液浸対物レンズであってもよい。液浸対物レンズとは、油浸レンズ、あるいは水浸レンズなどである。

対物レンズＯＧによって集光された光は、第１光学素子Ｌ１と第２光学素子Ｌ２と第３光学素子Ｌ３とを介して、撮像素子２８において結像される。この一例では、対物レンズＯＧによって集光された光は、第１光学素子Ｌ１を介して像ＶＲＩとして結像される。この像ＶＲＩは、ｘ軸方向とｚ０軸方向とを含む面についての細胞ＣＬの断面の像である。ここでｚ０軸とは、対物レンズＯＧによって集光された帯状の照明光ＬＳが細胞ＣＬに入射する方向である。帯状の照明光ＬＳが細胞ＣＬに入射する方向は、ｚ軸方向に対して所定の角度だけ傾いている。つまり、ｚ０軸はｚ軸に対して所定の角度だけ傾いている。撮像素子２８は、像ＶＲＩを、第２光学素子Ｌ２及び第３光学素子Ｌ３を介して撮像する。

［断面について］
ここで、図９を参照して、撮像素子２８が撮像する細胞ＣＬの断面の一例について説明する。
図９は、撮像素子２８が撮像する細胞ＣＬの断面を示す図である。図９（ａ）は、照明光ＬＳと、細胞ＣＬとの位置の関係を示す図である。光源２７から照射される照明光ＬＳは、ｚ０軸方向に照射される。
図９（ｂ）は、図９（ａ）から見た細胞ＣＬを、ｙ軸方向から見た図である。照明光ＬＳは、ｚ０軸方向とｘ軸方向とに帯状に広がる。細胞ＣＬは流路２１を流されることにより、この帯状に広がった照明光ＬＳを通過する。ここで、図９（ａ）に示す照明光ＬＳの幅Ｗ１は、図９（ｂ）に示す照明光ＬＳの幅Ｗ２よりも狭い。具体的には、幅Ｗ１は、２～３μｍである。
図９（ｃ）は、撮像素子２８と細胞ＣＬの断面の像である像ＶＲＩとの位置関係示す図である。撮像素子２８と像ＶＲＩとは平行である。撮像素子２８は、像ＶＲＩに垂直な方向（つまり、図２に示した第２光学素子Ｌ２、及び第３光学素子Ｌ３の光軸の方向）から像ＶＲＩを撮像する。

図２に戻り、イメージングフローサイトメーター２０の構成の説明を続ける。
撮像素子２８は、細胞ＣＬの断面を撮像する。ここで、撮像素子２８は、一例として、縦方向及び横方向に画素が配列されたカメラである。撮像素子２８は、縦方向及び横方向に配列された画素のうち、像が結像される領域の画素列の光の強さのみを取得することができる。例えば、撮像素子２８は、縦方向にＮ行、横方向にＮ列の画素が配列されている場合に、縦方向について中央のＭ行（ＭはＮより小さい数）分の画素列の光の強さのみを取得する。つまり、撮像素子２８は、Ｍ×Ｎ個の画素の光の強さのみを取得する。これにより、撮像素子２８は、像が結像されない領域との両方の画素の光の強さを取得する撮像素子と比較して、像が結像されない領域の画素列の光の強さを取得する時間を短縮することができる。また、撮像素子２８から取得した光の強さに基づいて、撮像画像を生成する画像処理の時間を短縮することができる。これにより、撮像素子２８は、高速に像を取得することができる。
なお、撮像素子２８は、像が結像される領域と、像が結像されない領域との両方の画素の光の強さを取得するカメラであってもよい。

図１０に、撮像素子２８において像が結像される領域の画素列の一例を示す。複数の像ＶＲＩそれぞれは、複数の細胞ＣＬｓそれぞれの像である。上述したように、複数の細胞ＣＬｓは、定在波である音波ＡＷが流路２１に照射されることによって、流路２１の高さ方向について整列されている。そのため、撮像素子２８は、複数の細胞ＣＬｓが幅方向について１列に整列された状態で複数の細胞ＣＬｓを撮像する。イメージングフローサイトメーター２０では、フレーム当たりの観測対象物の数を、観測対象物が高さ方向について１列に整列されていない場合に比べて多くできる。

本実施形態では、撮像素子２８は、ｓＣＭＯＳ（Ｓｃｉｅｎｔｉｆｉｃ ＣＭＯＳ；科学計測用ＣＭＯＳ）などから構成される撮像素子である。ｓＣＭＯＳは、従来のＣＣＤや、ＣＭＯＳによって構成される撮像素子よりも、高速かつ画質よく撮像することができる。
撮像素子２８は、撮像した撮像画像を、画像取得部２４に対して供給する。

制御部２００は、例えばＣＰＵや、ＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＦＰＧＡ（ｆｉｅｌｄ－ｐｒｏｇｒａｍｍａｂｌｅ ｇａｔｅ ａｒｒａｙ）などを備えており、種々の演算や情報の授受を行う。制御部２００は、画像取得部２４を、その機能部として備える。
画像取得部２４は、撮像素子２８から撮像画像を取得する。画像取得部２４は、撮像素子２８から取得した撮像画像を、断面画像ＰＩＣとして、記憶部２５に記憶させる。記憶部２５には、断面画像ＰＩＣが撮像された順に記憶される。

３次元画像生成部２６は、前記撮像部が細胞ＣＬの断面ＣＳを時系列に撮像して得られる複数の断面画像ＰＩＣに基づいて、細胞ＣＬの３次元画像を生成する。細胞ＣＬの３次元画像は、細胞ＣＬが流路２１を流れる方向である流線方向に細胞ＣＬが走査された画像である。３次元画像生成部２６は、記憶部２５から複数の断面画像ＰＩＣを取得する。３次元画像生成部２６は、記憶部２５から取得した複数の断面画像ＰＩＣに基づいて、細胞ＣＬの３次元画像を生成する。具体的には、３次元画像生成部２６は、断面画像ＰＩＣが撮像された順に－ｙ軸方向に向かって結合する。この断面画像ＰＩＣは、上述した面ＬＰを撮像した画像である。

ここで、図１１を参照して、３次元画像生成部２６が断面画像ＰＩＣを結合する順番について詳細に説明する。
図１１は、３次元画像生成部２６が断面画像ＰＩＣを結合する順番を示す図である。
図１１（ａ）は、細胞ＣＬの断面と断面画像と対応関係の一例を示す図である。
細胞ＣＬは、＋ｙ軸方向に向かって移動する。撮像部２３は、細胞ＣＬの断面ＣＳ１から断面ＣＳｎまでのｎ枚の断面画像ＰＩＣを、この順に時系列に撮像する。断面ＣＳｎのｎとは、１以上の整数である。
記憶部２５には、断面ＣＳ１の断面画像ＰＩＣ１から、断面ＣＳｎの断面画像ＰＩＣｎまでが記憶される。

図１１（ｂ）は、３次元画像生成部２６が断面画像ＰＩＣ１から断面画像ＰＩＣｎまでを結合する順番を示す図である。３次元画像生成部２６は、断面画像ＰＩＣが撮像された順に－ｙ軸方向に向かって結合する。これは、細胞ＣＬが＋ｙ軸方向に移動するためである。なお、結合の順番は、細胞ＣＬの移動方向に応じて変更する。
図１１（ｂ）には、３次元直交座標系として、ＸＹＺ座標系を示す。このＸＹＺ座標系は、３次元画像生成部２６が生成する３次元画像の座標系である。３次元画像生成部２６は、このＸＹＺ座標系と、ｘｙｚ座標系とを対応付けた状態で、断面画像ＰＩＣ１から断面画像ＰＩＣｎを結合する。具体的には、３次元画像生成部２６は、Ｘ方向と、ｘ方向とを同じ方向にして結合する。３次元画像生成部２６は、Ｙ軸方向とｙ軸方向とを同じ方向にして結合する。３次元画像生成部２６は、Ｚ軸方向とｚ軸方向とを同じ方向にして結合する。また、３次元画像生成部２６は、断面画像ＰＩＣ１から断面画像ＰＩＣｎを順に、－Ｙ軸方向に結合することにより、細胞ＣＬの３次元画像を生成する。

［イメージングフローサイトメーター２０の動作の概要］
次に、図１２を参照して、イメージングフローサイトメーター２０の動作手順の概要について説明する。
図１２は、イメージングフローサイトメーター２０の動作の一例を示す流れ図である。なお、ここに示す動作手順は、一例であって、動作手順の省略や動作手順の追加が行われてもよい。

（ステップＳ１０）音響素子２２は、流路２１に対して音波ＡＷを照射する。
（ステップＳ２０）撮像素子２８は、常に焦点位置ＦＰの像を撮像する。画像取得部２４は、撮像素子２８から信号を取得する。信号とは、上述した断面画像ＰＩＣを示す信号である。
（ステップＳ３０）画像取得部２４は、撮像素子２８から取得した信号に変化があるか否か判定する。具体的には、画像取得部２４には、流路を細胞ＣＬが流されていない状態の信号である検出前信号が記憶される。画像取得部２４は、検出前信号と、撮像素子２８から取得した信号とを比較する。具体的には、画像取得部２４は、検出前信号と撮像素子２８から取得した信号とが、所定の差が無い場合には、信号に変化が無いと判定する。画像取得部２４は、検出前信号と撮像素子２８から取得した信号とが、所定の差がある場合には、信号に変化があると判定する。

（ステップＳ３０；ＮＯ）画像取得部２４は、信号に変化が無いと判定する場合には、ステップＳ１０から処理を繰り返す。
（ステップＳ３０；ＹＥＳ）画像取得部２４は、信号に変化があると判定する場合には、撮像素子２８から取得した信号を、断面画像ＰＩＣとして、記憶部２５に記憶させる（ステップＳ４０）。

（ステップＳ５０）画像取得部２４は、撮像素子２８から信号を取得する。
（ステップＳ６０）画像取得部２４は、撮像素子２８から取得した信号に基づいて、細胞ＣＬの断面を撮像し終えたか判定する。具体的には、画像取得部２４は、検出前信号と撮像素子２８から取得した信号とが、所定の差がある場合には、細胞の断面を撮像し終えていないと判定する。画像取得部２４は、検出前信号と撮像素子２８から取得した信号とが、所定の差が無い場合には、細胞の断面を撮像し終えたと判定する。

（ステップＳ６０；ＮＯ）画像取得部２４は、細胞の断面を撮像し終えていないと判定する場合には、ステップＳ４０から処理を繰り返す。
（ステップＳ６０；ＹＥＳ）画像取得部２４は、細胞の断面を撮像し終えたと判定する場合には、ステップＳ７０の処理を実行する。

（ステップＳ７０）３次元画像生成部２６は、記憶部２５から複数の断面画像ＰＩＣを取得する。３次元画像生成部２６は、記憶部２５から取得した複数の断面画像ＰＩＣを、結合することにより、３次元画像を生成する。
（ステップＳ８０）３次元画像生成部２６は、生成した３次元画像を表示部１０に表示させる。イメージングフローサイトメーター２０は、処理を終了する。

ここで図１３及び図１４を参照し、実施形態に係るイメージングフローサイトメーター２０による音響効果について説明する。
図１３は、本実施形態に係る流路２１の高さ方向（ｚ軸方向）から複数の細胞が撮像された明視野画像の一例を示す図である。図１３（ａ）は、音波ＡＷが照射されていない場合の明視野画像であり、図１３（ｂ）は、音波ＡＷが照射された場合の明視野画像である。明視野画像において、複数の細胞は、高さ方向（ｚ軸方向）の位置に応じたコントラスト（明るさ）によって撮像されている。図１３に示す撮像結果では、流路２１として高さが２００マイクロメートル、幅が２ミリメートルのガラス毛細管が用いられている。流路２１には、蛍光染色された複数の細胞を１ミリリットル当たり１×１０７個含む懸濁液が毎分１０マイクロリットルの流速で流される。

図１３（ａ）に示すように、音波ＡＷが照射されていない場合、複数の細胞はコントラストがばらついており、流路２１の高さ方向（ｚ軸方向）の位置、及び速度がばらついている。また、音波ＡＷが照射されていない場合、複数の細胞には撮像過程において回転してしまったことが確かめられている。一方、図１３（ｂ）に示すように、音波ＡＷが照射された場合、複数の細胞はコントラストがばらついておらず、流路２１の高さ方向（ｚ軸方向）の位置及び速度がばらついておらず、かつ回転もみられなかった。

図１４は、本実施形態に係る高さ方向から所定の角度だけ傾いた面において複数の細胞が撮像された蛍光画像の一例を示す図である。高さ方向から所定の角度だけ傾いた面とは、図９において説明した面ＬＰである。図１４（ａ）は、音波ＡＷが照射されていない場合の蛍光画像であり、図１４（ｂ）は、音波ＡＷが照射された場合の蛍光画像である。図１４に示す撮像結果は、撮像素子２８の画素のうち像が結像される領域に対応する８００×２５６０個の画素を用いて毎秒２０フレームの撮像が行われたものである。
図１４（ａ）に示すように、音波ＡＷが照射されていない場合、複数の細胞は重力のために流路２１の高さ方向について中央よりも底に近い領域を流れている。一方、図１４（ｂ）に示すように、音波ＡＷが照射された場合、複数の細胞は流路２１の高さ方向について中央を流れている。

ここで図１５から図１８を参照し、イメージングフローサイトメーター２０による細胞の３次元画像の撮像結果（イメージング）について説明する。
図１５は、本実施形態に係る高さ方向から所定の角度だけ傾いた面（面ＬＰ）が異なる時刻において撮像された断面画像である。図１５に示す断面画像は、毎秒７２０フレームで撮像が行われている。この撮像では、流路２１の高さが２００マイクロメートルに対して、撮像素子２８の画素のうち縦方向について３６マイクロメートルの高さに対応する領域の撮像に用いられる画素のみが用いられている。
図１６は、図１５に示す断面画像から再構成された３次元画像の一例を示す図である。図１６は、毎秒７２０フレームで撮像された断面の画像について、６００フレームの画像から再構成された３次元画像である。図１６に示す３次元画像では、１９２１個の細胞が８３３ミリ秒をかけて撮像されている。スループットとして、毎秒２３０５個の細胞が撮像されたことに対応する。
図１７は、ｘ方向とｙ方向についての断面の画像である。図１８は、本実施形態に係る図１７に示した断面画像の各領域（領域Ｒ１、領域Ｒ２、領域Ｒ３、領域Ｒ４）の拡大図である。

（第２の実施形態）
次に、図１９を参照して、イメージングフローサイトメーターの第２の実施形態について説明する。なお、上述した第１の実施形態と同一の構成及び動作については、同一の符号を付してその説明を省略する。
図１９は、第２の実施形態の細胞測定システム１ａの機能構成の一例を示す図である。細胞測定システム１ａは、２次元画像を生成する点において、上述した第１の実施形態と異なる。

イメージングフローサイトメーター２０ａは、流路２１と、音響素子２２と、撮像部２３ａと、記憶部２５と、２次元画像生成部２６ａとを備える。
撮像部２３ａは、光源２７ｂと、対物レンズＯＧと、撮像素子２８と、第１光学素子Ｌ１と、第２光学素子Ｌ２と、第３光学素子Ｌ３と、制御部２００とを備える。

撮像部２３ａは、ラインスキャンカメラである。ラインスキャンカメラは、撮像素子２８の縦方向又は横方向に並んだ１列の画素列について、検出された光の強さを取得する。撮像部２３ａは、観測対象物の断面に含まれる線を時系列に撮像する。この線は、断面をある方向について線の集まりとみなした場合に、当該集まりに含まれるいずれかの線である。撮像素子２８としては、例えば、リニアアレイ型の光電子増倍管（ＰＭＴ：Ｐｈｏｔｏｍｕｌｔｉｐｌｉｅｒ Ｔｕｂｅ）などの光検出素子、あるいはｓＣＭＯＳカメラの一部の画素が用いられる。
２次元画像生成部２６ａは、撮像部２３ａが時系列に撮像した複数の線の画像に基づいて、観測対象物が流線方向に走査された２次元画像が生成する。

ここで、図２０を参照して、２次元画像生成部２６ａが線画像ＰＬＣを結合する順番について詳細に説明する。
図２０は、２次元画像生成部２６ａが線画像ＰＬＣを結合する順番を示す図である。
図２０（ａ）は、細胞ＣＬの断面に含まれる線と線画像と対応関係の一例を示す図である。
細胞ＣＬは、＋ｙ軸方向に向かって移動する。撮像部２３は、細胞ＣＬの断面ＣＳ１から断面ＣＳｎまでにそれぞれ含まれる線についてｎ枚の線画像ＰＬＣを、この順に時系列に撮像する。断面ＣＳｎのｎとは、１以上の整数である。
記憶部２５には、断面ＣＳ１に含まれる線の線画像ＰＬＣ１から、断面ＣＳｎに含まれる線の線画像ＰＬＣｎまでが記憶される。

図２０（ｂ）は、２次元画像生成部２６ａが線画像ＰＬＣ１から線画像ＰＬＣｎまでを結合する順番を示す図である。２次元画像生成部２６ａは、線画像ＰＬＣが撮像された順に－ｙ軸方向に向かって結合する。
図２０（ｂ）には、３次元直交座標系として、ＸＹＺ座標系を示す。このＸＹＺ座標系に含まれるＸ軸とＹ軸とからなる２次元座標系であるＸＹ座標系は、２次元画像生成部２６ａが生成する２次元画像の座標系である。２次元画像生成部２６ａは、このＸＹ座標系と、ｘｙｚ座標系のうちｘ軸とｙ軸とを対応付けた状態で、線画像ＰＬＣ１から線画像ＰＬＣｎを結合する。具体的には、２次元画像生成部２６ａは、Ｘ方向と、ｘ方向とを同じ方向にして結合する。２次元画像生成部２６ａは、Ｙ軸方向とｙ軸方向とを同じ方向にして結合する。また、２次元画像生成部２６ａは、線画像ＰＬＣ１から線画像ＰＬＣｎを順に、－Ｙ軸方向に結合することにより、細胞ＣＬの２次元画像を生成する。

［実施形態のまとめ］
以上に説明したように、上述した実施形態に係るイメージングフローサイトメーター２０、２０ａは、流路２１と、音響素子２２と、光源２７と、撮像部２３、２３ａ、２３ｂと、走査画像生成部（３次元画像生成部２６、２次元画像生成部２６ａ）とを備える。
流路２１は、観測対象物（上述した実施形態において、細胞ＣＬ）が流される流路であって、幅方向の長さ（上述した実施形態において、幅方向長さＦＷ）が高さ方向の長さ（上述した実施形態において、高さ方向長さＦＨ）よりも長い。
音響素子２２は、流路２１に対して音波ＡＷを定在波にして照射する。
光源２７、２７ｂは、流路２１に対して照明光ＬＳ、照明光ＬＳｂを照射する。
撮像部２３、２３ａ、２３ｂは、照明光ＬＳ、照明光ＬＳｂが照射される位置を通過した観測対象物（上述した実施形態において、細胞ＣＬ）を計測または撮像することにより、観測対象物（上述した実施形態において、細胞ＣＬ）が流路２１を流れる方向である流線方向についての観測対象物（上述した実施形態において、細胞ＣＬ）の断面ＣＳに含まれる線を少なくとも撮像する。
走査画像生成部（３次元画像生成部２６、２次元画像生成部２６ａ）は、撮像部２３、２３ａ、２３ｂが断面ＣＳに含まれる線を時系列に撮像して得られる複数の撮像画像（上述した実施形態において、断面画像ＰＩＣ、線画像ＰＬＣ）に基づいて、観測対象物（上述した実施形態において、細胞ＣＬ）が流線方向に走査された画像（上述した実施形態において、細胞ＣＬの３次元画像、細胞ＣＬの２次元画像）を生成する。流線方向は、観測対象物（上述した実施形態において、細胞ＣＬ）が流路２１を流れる方向である。

この構成により、上述した実施形態に係るイメージングフローサイトメーター２０、２０ａは、観測対象物が撮像される過程において、観測対象物の速度（進行方向と速さ）がばらついてしまうこと、観測対象物が回転してしまうことを音響効果によって抑制できるため、撮像条件を安定させることができる。撮像条件には、撮像過程における観測対象物の流線方向に垂直な方向における位置、回転、及び速度が含まれる。実施形態に係るイメージングフローサイトメーター２０、２０ａでは、観測対象物を流線方向に走査して走査画像を生成するため、観測対象物が撮像される過程において、観測対象物の流線方向に垂直な方向における位置、回転（観測対象物の向き）、あるいは速度が変化してしまうと、観測対象物を流線方向に走査できなくなってしまう。

また、実施形態に係るイメージングフローサイトメーター２０、２０ａでは、音響効果によって観測対象物が流路２１の高さ方向ついて撮像される位置以外の位置を通過ししてしまうことを抑制できるため、撮像部２３に備えられる撮像素子２８の画素列のうち流路２１の高さ方向の特定の位置に対応する画素列のみを用いて撮像を行うことができる。撮像部２３に備えられる撮像素子２８の画素列のうち流路２１の高さ方向の特定の位置に対応する画素列のみを用いて撮像を行う場合、観測対象物が流路２１の高さ方向ついて撮像される位置以外の位置を通過ししてしまうと、観測対象物の像自体が得られなくなってしまう。

また、上述したように撮像部２３は、カメラであり、撮像速度限界、つまり単位時間当たりに撮像可能なフレーム数に限界が存在する。そのため、イメージングフローサイトメーター２０、２０ａでは、流路２１を流れる流体の流速を低く抑えなければ、撮像対象に対して流れ方向の撮像回数が減り、解像度が落ちることになる。つまり、イメージングフローサイトメーター２０、２０ａでは、撮像部２３に応じて所望の流速は決まっている。
従来、観測対象物が流路２１の高さ方向についての位置を安定させるために、流体力学的フォーカシングが知られている。複数の流路を備えて、流体力学的フォーカシングを行う場合、流路２１を流れる流体の流速が、流体力学的フォーカシングを行う前に比べて速くなってしまい、流速が安定しない。また、撮像部に応じた流速が所望の流速よりも高くなりすぎるために撮像される像の流れ方向の解像度が低下してしまう。イメージングフローサイトメーター２０、２０ａでは、流速が速くなってしまう流体力学的フォーカシングとは異なり、音響効果によって流速を所望の流速から変化させることなく撮像が可能である。

イメージングフローサイトメーター２０、２０ａでは、流路を流れる観測対象物を流線方向に走査する場合に上述した撮像条件を安定させることができ、安定してかつ高速に走査画像を生成するイメージングが可能となる。

また、上述した実施形態に係るイメージングフローサイトメーター２０、２０ａでは、流路２１の高さ方向について音波ＡＷの定在波の節の数は１個であり、流路２１の幅方向について音波ＡＷの定在波の節の数は０個である。
この構成により、上述した実施形態に係るイメージングフローサイトメーター２０、２０ａでは、流路２１の幅方向について音波ＡＷの定在波の節の数は０個であるため、流路２１の幅方向に並べられる複数の細胞ＣＬｓの数は、当該節の数によっては制限されない。
なお、高さ方向について音波ＡＷの定在波の節の数が１個であり、流路２１の幅方向について音波ＡＷの定在波の節の数が１個以上とされてもよい。
なお、流路２１の高さ方向について音波ＡＷの定在波の節の数は、２個以上とされてもよい。

また、上述した実施形態に係るイメージングフローサイトメーター２０、２０ａでは、流路２１には、流路２１の幅方向に整列された複数の観測対象物（上述した実施形態において、複数の細胞ＣＬｓ）が流される。
この構成により、上述した実施形態に係るイメージングフローサイトメーター２０、２０ａは、流路２１に１個の観測対象物が流される場合に比べてフレーム当たりの観測対象物の数を多くできるため、単位時間当たりに、流路２１に１個の観測対象物が流される場合に比べて多い数の観測対象物の走査画像を生成できる。実施形態に係るイメージングフローサイトメーター２０、２０ａでは、スループット（単位時間当たりに計測可能な細胞の個数）と信号対ノイズ比を両立させることができる。スループットと信号対ノイズ比を両立させるとは、同じスループットで計測した際に、各瞬間で同時にたくさんの細胞を計測することによって、細胞当たりに掛けられる計測時間を増加させ、信号対ノイズ比を増加させることである。
また、上述した実施形態に係るイメージングフローサイトメーター２０、２０ａでは、複数の観測対象物は流路２１の幅方向に１列に整列されているため、撮像部２３に備えられる撮像素子２８として像が結像される領域の画素列の光の強さのみを取得する撮像素子を採用することが好適である。そのような撮像素子では、流路２１の高さ方向について、複数の観測対象物が整列されている位置以外の像が結像されない領域の画素列の光の強さは取得されないため、撮像画像を生成する画像処理の時間を短縮することができる。

なお、上述した実施形態に係るイメージングフローサイトメーター２０、２０ａにおいては、帯状の照明光ＬＳが用いられる場合の一例について説明したが、これに限られない。流路を流れる観測対象物の断面を撮像して３次元情報、または２次元情報を撮像する他のイメージング技術が用いられてもよい。照明光は帯状でなくてもよいし、照明光は、流路の高さ方向に対して傾いていなくてもよい。

なお、上述した実施形態に係るイメージングフローサイトメーター２０、２０ａにおいて撮像される観測対象物が流線方向に走査された３次元画像、または２次元画像に基づいて、イメージングフローサイトメーター２０、２０ａは観測対象物のソーティングをしてもよい。ソーティングとは、流路２１を流される観測対象物のうちから、所定の観測対象物を分取することである。この所定の観測対象物は、ユーザーによって予め選択されていればよい。

ここでソーティングについて説明する。例えば、観測対象物として所定の細胞と、ゴミや他の細胞などの所定の細胞とは異なる物とが、流路２１に流されることがある。観測対象物のうちから所定の細胞を選択して取り出すことが、ソーティングである。

つまり、イメージングフローサイトメーター２０、２０ａは、ユーザーによって予め選択された細胞の形態を示す情報と、流線方向に走査された３次元画像、または２次元画像に含まれる細胞の形態を示す情報とを比較することにより、ソーティング対象の細胞か否かを判定し、分取する。
イメージングフローサイトメーター２０、２０ａでは、流路２１を流れる流体の流速は、ソーティング対象の判定、及び分取に要する時間に対して、十分に低速である。

（実施例）
以下では、上述した第１の実施形態に係るイメージングフローサイトメーター２０を用いて細胞の３次元画像の生成（イメージング）を実施した例について説明する。
本実施例では、観測対象物として、有糸分裂期のヒトの白血病細胞（Ｋ５６２細胞）である細胞Ｃ１を解析した。細胞Ｃ１は、抗有糸分裂性のタンパク質であるＭＰＭ－２と、ＤＡＰＩとによって染色されている。当該染色された細胞Ｃ１を１ミリリットル当たり約１×１０^７個の濃度で含む懸濁液を用いた。当該懸濁液を、毎分１０マイクロリットルの流量でイメージングフローサイトメーター２０に備えられる流路２１に流した。

図２１は、上記の実験条件でイメージングフローサイトメーター２０によって８５７ミリ秒間の観察の結果、生成された細胞Ｃ１の３次元画像を２次元平面（ｘｙ平面）に射影した画像である。図２１に示す画像は、３次元画像に含まれるｚ方向の位置が異なる複数の断面の画像が、重ねて表示された画像である。
図２１に示す３次元画像には、１２４８個の細胞Ｃ１が３次元画像として再構築されている。図２１に示す３次元画像から、この母集団において６０個の細胞Ｃ１がＭＰＭ－２陽性であると判定された。ＤＡＰＩによって染色された細胞Ｃ１の３次元画像の解析によって、３４個の細胞Ｃ１が細胞分裂の前期、１０個の細胞Ｃ１が中期、２個の細胞Ｃ１が後期であり、２組の細胞Ｃ１が終期であること判定された。

図２１において、領域Ｒ１１、領域Ｒ１２、領域Ｒ１３、及び領域Ｒ１４にそれぞれ含まれる画像は前期、中期、後期、及び終期それぞれの細胞Ｃ１の画像である。領域Ｒ１１、領域Ｒ１２、及び領域Ｒ１４にそれぞれ含まれる細胞の画像は、ＤＡＰＩによって染色された細胞Ｃ１の画像である。領域Ｒ１３に含まれる細胞の画像は、ＭＰＭ－２によって染色された細胞Ｃ１の画像である。図２２に、領域Ｒ１１、領域Ｒ１２、領域Ｒ１３、及び領域Ｒ１４それぞれを拡大した画像を示す。

図２３及び図２４に、ＤＡＰＩによって染色された細胞Ｃ１の３次元画像を２次元平面に射影した画像を、細胞分裂の中期及び後期のそれぞれについて示す。図２３（Ａ）、（Ｂ）、及び（Ｃ）はそれぞれ、ｘｙ平面、ｙｚ平面、及びｘｚ平面それぞれに射影された中期の細胞Ｃ１の画像である。図２４（Ａ）、（Ｂ）、及び（Ｃ）はそれぞれ、ｘｙ平面、ｙｚ平面、及びｘｚ平面それぞれに射影された後期の細胞Ｃ１の画像である。ここで図２４（Ｂ）、及び（Ｃ）から、後期の細胞Ｃ１のｙｚ平面、及びｘｚ平面それぞれに射影された画像について、染色体の分裂を確認することができる。一方、図２４（Ａ）に示した後期の細胞Ｃ１のｘｙ平面に射影された画像では、分裂した染色体同士が重なって撮像されており、当該画像から染色体の分裂を確認することは難しい。

図２４（Ａ）に示したｘｙ平面に射影された画像は、従来の２次元イメージングによって得られる２次元画像に相当する画像である。従来の２次元イメージングでは中期板や染色体の分裂などの細胞周期まで判別可能な画像として生成することは、細胞の向きによっては難しかった。一方、イメージングフローサイトメーター２０による３次元イメージングによれば、図２４（Ｂ）及び（Ｃ）に示すように、染色体のｚ方向の位置が識別することができた。このように、イメージングフローサイトメーター２０では、ｚ方向も含めた３次元イメージングが可能であるため、中期板や染色体の分裂などの細胞周期まで判別可能な３次元画像を生成することができる。

ここで細胞Ｃ１のうちいくつかの細胞では、ｘｙ平面、ｙｚ平面、及びｘｚ平面それぞれの平面の画像では、細胞分裂の終期を判定するために染色体の配置を識別するために十分ではない。図２５から図２７を参照し、３次元画像を射影する方向を変更することによって染色体の配置を識別できるようになる場合があることについて説明する。

図２５は、イメージングフローサイトメーター２０によって生成された３次元画像のうちＤＡＰＩによって染色された領域の一例を示す。当該細胞Ｃ１は、細胞分裂の中期の細胞である。図２５には、ｘ軸、ｙ軸、及びｚ軸が示されている。これらのｘ軸、ｙ軸、及びｚ軸それぞれが示す方向は、例えば、図１に示したｘｙｚ座標系の座標軸それぞれが示す方向に対応する。図２５には、ｘ軸、ｙ軸、及びｚ軸とともに、λ_１軸、λ_２軸、λ_３軸が示されている。λ_１軸、λ_２軸、λ_３軸は、主成分分析によって得られる固有ベクトルに基づいてｘ軸、ｙ軸、及びｚ軸が回転されて得られた座標軸である。なお、λ_１軸、λ_２軸、λ_３軸それぞれの方向を示す単位ベクトルはそれぞれ、元のｘｙｚ座標系において［－０．５４，０．８２，０．２０］、［０．７１０．３１０．６３］、及び［０．４５，０．４８，－０．７５］である。

図２５に示す領域を、ｘｙ平面、ｙｚ平面、及びｘｚ平面それぞれに射影した場合、図２６に示す２次元画像が得られる。図２６（Ａ）、（Ｂ）、及び（Ｃ）に示すように、ｘｙ平面、ｙｚ平面、及びｘｚ平面の２次元画像では、染色体の配置を識別することは難しい。
図２７は、図２５に示す領域を、軸をｘ軸、ｙ軸、ｚ軸からλ_１軸、λ_２軸、λ_３軸へと回転させて、平面に射影した場合の２次元画像である。図２７（Ａ）、（Ｂ）、及び（Ｃ）はそれぞれ、λ_１λ_２平面、λ_１λ_３平面、λ_２λ_３平面それぞれに射影された２次元画像を示す。図２７（Ａ）では、染色体の配置が識別可能となっている。

イメージングフローサイトメーター２０では、３次元画像を生成しているため、観測対象物の構造を識別するために３次元画像を２次元画像に射影する方向を選択することができる。射影する方向を選択することは、イメージングフローサイトメーター２０が３次元画像を生成しているために可能である。２次元画像や、予め決められた複数の方向それぞれについての２次元画像を生成する場合には、画像を生成した後に解析する段階において、射影する方向を選択することはできない。

図２８に、上述した図２１に示した画像で撮像された範囲よりも広い範囲について３次元画像を解析した結果を示す。図２８では、流路２１を流れる１０^５個のオーダーの個数の細胞Ｃ１の３次元画像がイメージングフローサイトメーター２０によって生成されている。当該細胞Ｃ１は、ＭＰＭ－２またはＤＡＰＩによって染色されている。得られた３次元画像では、２８６秒の合計撮像時間で４０８９３７個の細胞Ｃ１の画像が検出された。図２８に示す画像では、ＭＰＭ－２及びＤＡＰＩによって染色された細胞Ｃ１の全体の１パーセントの３次元画像を２次元平面（ｘｙ平面）に射影した画像である。ｘｙ平面に射影された画像である。図２８に示す画像は、３次元画像に含まれるｚ方向の位置が異なる複数の断面の画像が、重ねて表示された画像である。スケールバーは、１００マイクロメートルの長さを示す。

まず、細胞Ｃ１それぞれについてのＭＰＭ－２及びＤＡＰＩの強度を解析した。図２９は、当該解析の結果を示す散布図である。当該解析によって、ＭＰＭ－２陽性の細胞Ｃ１が選別された。領域Ｒ２０は、選別されたＭＰＭ－２陽性の細胞Ｃ１を示す。
次に、領域Ｒ２０において、細胞Ｃ１内においてＤＡＰＩ染色された領域の数が計測された。当該計測の結果、２つのＤＡＰＩ染色された領域が計測された。図３０は、２つのＤＡＰＩ染色された領域の３次元画像がｘｙ平面に射影された２次元画像である。図３０に示す画像から細胞分裂の後期または終期であることが判別された。

１つのＤＡＰＩ染色された領域しかもたない細胞Ｃ１については、当該領域のアスペクト比を３次元画像から算出した。図３１（Ａ）に、１つのＤＡＰＩ染色された領域の３次元画像が平面に射影された２次元画像を示す。図３１（Ａ）では、ＤＡＰＩ染色された領域の長軸と短軸とのそれぞれの向きが垂直方向と水平方向とそれぞれ一致するように、３次元画像が射影される平面が選択されている。
比較のために図３１（Ｂ）に、１つのＤＡＰＩ染色された領域の３次元画像がｘｙ平面に射影された２次元画像を示す。図３１（Ｂ）に示す２次元画像は、２次元のイメージングを行うイメージングフローサイトメトリーによって生成される２次元画像を模したものである。つまり図３１（Ｂ）に示す２次元画像は、２次元イメージングに相当する画像である。
図３１（Ｂ）に示すＤＡＰＩ染色された領域では、図３１（Ａ）に示すＤＡＰＩ染色された領域に比べてアスペクト比が大きくなっている。

図３１（Ａ）に示すような３次元イメージングによる画像から算出されたアスペクト比と、図３１（Ｂ）に示すような２次元イメージングに相当する画像から算出されたアスペクト比とを比較した。図３２に比較結果を示す。図３２では、複数の細胞についてＤＡＰＩ領域のアスペクト比の分布を示すヒストグラムが示されている。ヒストグラムＨ１は、３次元イメージングによる画像から得られたアスペクト比を示す。ヒストグラムＨ２は、２次元イメージングに相当する画像から得られたアスペクト比を示す。図３２から、２次元イメージングに相当する画像では、３次元イメージングによる画像に比べてアスペクト比が過大に算出される傾向があることがわかる。当該傾向は、２次元イメージングでは、３次元イメージングに比べて次元が少ないことが原因であると考えられる。したがって、細胞の構造のプロファイルを正確に得るためには、高いスループットでの３次元イメージングが有用であると考えられる。イメージングフローサイトメーター２０は、そのような３次元イメージングを可能とする。

以上、本発明の実施形態を、図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、本発明の趣旨を逸脱しない範囲で適宜変更を加えることができる。

なお、上述の各装置は内部にコンピュータを有している。そして、上述した装置の各処理の過程は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータが読み出して実行することによって、上記処理が行われる。ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ－ＲＯＭ、ＤＶＤ－ＲＯＭ、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしてもよい。

また、上記プログラムは、前述した機能の一部を実現するためのものであってもよい。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であってもよい。

２０、２０ａ…イメージングフローサイトメーター、２１…流路、２２…音響素子、２７、２７ｂ…光源、２３…撮像部、２６…３次元画像生成部、２６ａ…２次元画像生成部、ＣＬ…細胞、ＬＳ、ＬＳｂ…照明光、ＡＷ…音波、ＰＩＣ…断面画像、ＰＬＣ…線画像

Claims (4)

1. 観測対象物が流される流路であって、幅方向の長さが高さ方向の長さよりも長い流路と、
   前記流路に対して音波を定在波にして照射する音響素子と、
   前記流路に対して照明光を照射する光源と、
   前記照明光が照射される位置を通過した前記観測対象物を計測または撮像することにより、前記観測対象物が前記流路を流れる方向である流線方向についての前記観測対象物の断面に含まれる線を少なくとも撮像する撮像部と、
   前記撮像部が前記線を時系列に撮像して得られる複数の撮像画像に基づいて、前記観測対象物が前記流線方向に走査された画像を生成する走査画像生成部と、
   を備え、
   前記流路を流れる前記観測対象物の流速は前記撮像部に応じて決まっており、
   前記音響素子と前記撮像部との前記流路の長さ方向についての距離は、前記音響素子による前記音波の音響効果による複数の前記観測対象物の整列が持続する程度に短く、
   前記高さ方向について前記定在波の節の数は所定の個数であり、
   前記撮像部は、前記撮像部に備えられる撮像素子の画素列のうち前記高さ方向の前記定在波の節の位置に対応する画素列のみを用いて撮像を行う
   イメージングフローサイトメーター。
2. 前記高さ方向について前記定在波の節の数は１個であり、前記幅方向について前記定在波の節の数は０個である
   請求項１に記載のイメージングフローサイトメーター。
3. 前記高さ方向について前記定在波の節の数は１個であり、前記幅方向について前記定在波の節の数は１個以上である
   請求項１に記載のイメージングフローサイトメーター。
4. 前記流路には、前記幅方向に整列された複数の前記観測対象物が流される
   請求項１から請求項３のいずれか一項に記載のイメージングフローサイトメーター。

Images