(57)【要約】  反射干渉計測用光源１０６と、定量位相計測用光源１０５と、反射干渉計測用光源１０６から放射され細胞１０１から反射される光を撮像することにより、反射干渉画像を生成する反射干渉検出用カメラ１１０と、定量位相計測用光源１０５から放射され細胞１０１を透過する光を撮像することにより、定量位相画像を生成する定量位相検出用カメラ１１２と、反射干渉画像と定量位相画像との間で、空間的な位置を一致させることにより、両画像の位置合わせを行う画像位置合わせ部２０１と、定量位相画像との間で位置合わせが行われた反射干渉画像より、第１パラメータを抽出する第１抽出部２０４と、反射干渉画像との間で位置合わせが行われた定量位相画像より、第２パラメータを抽出する第２抽出部２０５と、を備える。
2011132587.tif 000002
【特許請求の範囲】
【請求項１】
反射干渉計測用光源と、
定量位相計測用光源と、
前記反射干渉計測用光源から放射され細胞から反射される光を撮像することにより、反射干渉画像を生成する反射干渉撮像手段と、
前記定量位相計測用光源から放射され前記細胞を透過する光を撮像することにより、定量位相画像を生成する定量位相撮像手段と、
前記反射干渉画像と前記定量位相画像との間で、空間的な位置を一致させることにより、両画像の位置合わせを行う画像位置合わせ手段と、
前記画像位置合わせ手段により前記定量位相画像との間で前記位置合わせが行われた前記反射干渉画像より、第１パラメータを抽出する第１抽出手段と、
前記画像位置合わせ手段により前記反射干渉画像との間で前記位置合わせが行われた前記定量位相画像より、第２パラメータを抽出する第２抽出手段と、
を備えることを特徴とする細胞観察装置。
【請求項２】
前記定量位相画像に基づき、前記細胞の輪郭を抽出する輪郭抽出手段と、
前記輪郭抽出手段が抽出した前記輪郭を前記反射干渉画像に適用することにより、輪郭適用後の反射干渉画像を生成する輪郭適用手段と、
前記輪郭適用後の反射干渉画像より、第３パラメータを抽出する第３抽出手段と、
を更に備えることを特徴とする請求項１に記載の細胞観察装置。
【請求項３】
前記第１パラメータは、前記細胞が置かれている基板と、前記細胞との間の接着状態に基づく情報である、
ことを特徴とする請求項１または２に記載の細胞観察装置。
【請求項４】
前記第２パラメータは、前記細胞の光学的厚さ、面積、体積、または前記細胞内での屈折率の変化に基づく情報である、
ことを特徴とする請求項１〜３の何れか１項に記載の細胞観察装置。
【請求項５】
前記第３パラメータは、前記細胞の前記輪郭の中での、前記基板への前記接着状態に基づく情報である、
ことを特徴とする請求項２〜４の何れか１項に記載の細胞観察装置。
【請求項６】
予め既知の種類または状態の前記細胞に対して抽出したパラメータをリファレンスデータとして記憶するリファレンス記憶手段と、
前記リファレンスデータに基づき、未知の細胞の種類または状態を判別する解析手段と、
を更に備えることを特徴とする請求項１〜５の何れか１項に記載の細胞観察装置。
【請求項７】
未知の細胞に対して抽出したパラメータの中から所定のパラメータを選択し、当該選択された所定のパラメータを用いて、前記未知の細胞の種類または状態を判別する解析手段を更に備える、
ことを特徴とする請求項１〜５の何れか１項に記載の細胞観察装置。
【請求項８】
前記第１抽出手段または前記第２抽出手段が３つ以上のパラメータを抽出した場合に、前記３つ以上のパラメータに対して主成分分析を行うことにより、未知の細胞の種類または状態を判別する解析手段を更に備える、
ことを特徴とする請求項１〜５の何れか１項に記載の細胞観察装置。
【請求項９】
前記定量位相計測用光源として、低コヒーレント光を使用する、
ことを特徴とする請求項１〜８の何れか１項に記載の細胞観察装置。
【請求項１０】
前記反射干渉計測用光源として、低コヒーレント光を使用する、
ことを特徴とする請求項１〜９の何れか１項に記載の細胞観察装置。
【請求項１１】
前記反射干渉計測用光源から放射され前記細胞から反射される光が集光される対物レンズと、
前記対物レンズの前記反射干渉計測用光源側の開口絞りと共役な位置に、リング状のスリットと、
を更に備えることを特徴とする請求項１〜１０の何れか１項に記載の細胞観察装置。
【請求項１２】
前記細胞を収容する容器を更に備え、
前記容器の前記細胞の接着面の反対側に反射防止コートが施されている、
ことを特徴とする請求項１〜１１の何れか１項に記載の細胞観察装置。
【請求項１３】
前記反射干渉計測用光源と前記反射干渉撮像手段との間および前記定量位相計測用光源と前記反射干渉撮像手段との間に位置し、波長により反射対透過の比率を異ならせることが可能なミラーを更に備える、
ことを特徴とする請求項１〜１２の何れか１項に記載の細胞観察装置。
【請求項１４】
反射干渉撮像手段が、反射干渉計測用光源から放射され細胞から反射される光を撮像することにより、反射干渉画像を生成する反射干渉撮像ステップと、
定量位相撮像手段が、定量位相計測用光源から放射され前記細胞を透過する光を撮像することにより、定量位相画像を生成する定量位相撮像ステップと、
画像位置合わせ手段が、前記反射干渉画像と前記定量位相画像との間で、空間的な位置を一致させることにより、両画像の位置合わせを行う画像位置合わせステップと、
第１抽出手段が、前記画像位置合わせ手段により前記定量位相画像との間で前記位置合わせが行われた前記反射干渉画像より、第１パラメータを抽出する第１抽出ステップと、
第２抽出手段が、前記画像位置合わせ手段により前記反射干渉画像との間で前記位置合わせが行われた前記定量位相画像より、第２パラメータを抽出する第２抽出ステップと、
を備えることを特徴とする細胞観察方法。
【請求項１５】
輪郭抽出手段が、前記定量位相画像に基づき、前記細胞の輪郭を抽出する輪郭抽出ステップと、
輪郭適用手段が、前記輪郭抽出手段が抽出した前記輪郭を前記反射干渉画像に適用することにより、輪郭適用後の反射干渉画像を生成する輪郭適用ステップと、
第３抽出手段が、前記輪郭適用後の反射干渉画像より、第３パラメータを抽出する第３抽出ステップと、
を更に備えることを特徴とする請求項１４に記載の細胞観察方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、細胞観察装置および細胞観察方法に関するものである。
【背景技術】
【０００２】
従来、細胞に例えば蛍光ラベルを付するなど、細胞の状態を侵襲的な方法により判別・評価することが行われていた。しかし、このような方法では、細胞に毒性のある色素を用いたり、細胞を固定液を用いて固定して処理したりするために、処理後に細胞を再利用することはできず、細胞を自然な状態で評価したり、評価後に培養を継続したり、移植に用いたりすることは難しい。したがって、非侵襲的な方法で細胞の状態を判別・評価することが好ましいといえる。特許文献１、および非特許文献１〜３は、非侵襲的な方法で細胞を判別・評価することについて開示している。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】国際公開ＷＯ２００５／０２９４１３パンフレット
【非特許文献】
【０００４】
【非特許文献１】Y.Usson,A.Guignandon,N. Laroche,M-H. Lafage-Proust,L.Vico,“Quantitation of Cell-Matrix Adhesion Using Confocal Image Analysis of FocalContact Associated Proteins and Interference Reflectioｎｍicroscopy.”, Cytometry,２８,２９８-３０４,（１９９７）
【非特許文献２】C.K.Choi, C.H.Margraves,A.E.English, K.D.Kihm,  “Multicontrast microscopy technique to dynamicallyfingerprint  live-cell focal contacts during exposure and replacement ofcytotxic medium.”, J. Biomedical Optics, １３（５）, （２００８）
【非特許文献３】I.Weber, R.Albrecht, “Imageprocessing for combined bright-field and reflection interference contrast videomicroscopy.”, Computer Methods and Programs in Biomedicine, ５３, １１３-１１８, （１９９７）
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１は、非侵襲的な方法で細胞の状態を判別・評価する方法の一つとして、定量位相顕微鏡を用いる方法について開示している。この方法では、細胞の存否による光路長の差を用いて、細胞の状態を非侵襲的に判別・評価できる。しかし、この方法によって得られる情報は、細胞の光学的厚さ、面積、体積、および細胞内の屈折率の変化に基づく情報のみに限られており、細胞の状態を適切に判別・評価する上で必要な情報量として決して十分とはいえない。【０００６】
非特許文献１〜３は、非侵襲的な方法の他の一例として、反射干渉顕微鏡を用いる方法について開示している。この方法では、個々の細胞の基板への接着状態に基づき、基板と培養液との間の界面からの反射光と、培養液と細胞膜との間の界面からの反射光との干渉により明暗のコントラストが現れる現象を利用して、細胞の状態を非侵襲的に判別・評価できる。しかし、この方法によって得られる情報は、細胞の基板への接着状態に基づく情報のみに限られており、細胞の状態を適切に判別・評価する上で必要な情報量として決して十分とはいえない。【０００７】
このように、定量位相顕微鏡、反射干渉顕微鏡の何れか一方を単独で用いる方法では、細胞の光学的厚さ、面積、体積、および細胞内の屈折率の変化に基づく情報だけ、または細胞の基板への接着状態に基づく情報だけが得られる。それら単独の情報だけでは、細胞の状態を適切に判別・評価する上で必要な情報量として決して十分とはいえない。また、定量位相顕微鏡および反射干渉顕微鏡を一緒に用いることを開示または示唆するような文献は発見されておらず、両顕微鏡は機能が異なる別個の顕微鏡であることから、両顕微鏡を一緒に用いるには何らかの工夫が必要となる。【０００８】
そこで、本発明は上記に鑑みてなされたもので、細胞の状態を適切に判別・評価するためのより多くの情報を得ることを可能とする細胞観察装置および細胞観察方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記課題を解決するために、本発明の細胞観察装置は、反射干渉計測用光源と、定量位相計測用光源と、前記反射干渉計測用光源から放射され細胞から反射される光を撮像することにより、反射干渉画像を生成する反射干渉撮像手段と、前記定量位相計測用光源から放射され前記細胞を透過する光を撮像することにより、定量位相画像を生成する定量位相撮像手段と、前記反射干渉画像と前記定量位相画像との間で、空間的な位置を一致させることにより、両画像の位置合わせを行う画像位置合わせ手段と、前記画像位置合わせ手段により前記定量位相画像との間で前記位置合わせが行われた前記反射干渉画像より、第１パラメータを抽出する第１抽出手段と、前記画像位置合わせ手段により前記反射干渉画像との間で前記位置合わせが行われた前記定量位相画像より、第２パラメータを抽出する第２抽出手段と、を備えることを特徴とする。【００１０】
また、本発明の細胞観察方法は、反射干渉撮像手段が、反射干渉計測用光源から放射され細胞から反射される光を撮像することにより、反射干渉画像を生成する反射干渉撮像ステップと、定量位相撮像手段が、定量位相計測用光源から放射され前記細胞を透過する光を撮像することにより、定量位相画像を生成する定量位相撮像ステップと、画像位置合わせ手段が、前記反射干渉画像と前記定量位相画像との間で、空間的な位置を一致させることにより、両画像の位置合わせを行う画像位置合わせステップと、第１抽出手段が、前記画像位置合わせ手段により前記定量位相画像との間で前記位置合わせが行われた前記反射干渉画像より、第１パラメータを抽出する第１抽出ステップと、第２抽出手段が、前記画像位置合わせ手段により前記反射干渉画像との間で前記位置合わせが行われた前記定量位相画像より、第２パラメータを抽出する第２抽出ステップと、を備えることを特徴とする。【００１１】
このような本発明の細胞観察装置および細胞観察方法によれば、反射干渉計測用光源、反射干渉撮像手段、第１抽出手段（以下、「反射干渉計測ユニット」という。）を備えることにより、細胞からの反射光をもとに、第１パラメータを得る。また、定量位相計測用光源、定量位相撮像手段、第２抽出手段（以下、「定量位相計測ユニット」という。）を備えることにより、細胞からの透過光をもとに、第２パラメータを得る。このように、本発明の細胞観察装置および細胞観察方法では、反射干渉計測ユニットおよび定量位相計測ユニットの両方を備えることにより、第１パラメータおよび第２パラメータの両方を得ることができ、これにより使用者は細胞の状態を適切に判別・評価するためのより多くの情報を得ることとなる。更に、第１パラメータは定量位相画像との間で位置合わせが行われた反射干渉画像より抽出されるものであり、第２パラメータは反射干渉画像との間で位置合わせが行われた定量位相画像より抽出されるものであるため、両パラメータ間には細胞のどの部分に関するパラメータであるかについての整合性が取れているといえる。つまり、本発明の画像位置合わせ手段による両画像の位置合わせにより、両パラメータは細胞の同一部分の状態を判別・評価するためのパラメータとして活用可能となる。【００１２】
また、本発明において、前記定量位相画像に基づき、前記細胞の輪郭を抽出する輪郭抽出手段と、前記輪郭抽出手段が抽出した前記輪郭を前記反射干渉画像に適用することにより、輪郭適用後の反射干渉画像を生成する輪郭適用手段と、前記輪郭適用後の反射干渉画像より、第３パラメータを抽出する第３抽出手段と、を更に備えても良い。【００１３】
また、本発明において、輪郭抽出手段が、前記定量位相画像に基づき、前記細胞の輪郭を抽出する輪郭抽出ステップと、輪郭適用手段が、前記輪郭抽出手段が抽出した前記輪郭を前記反射干渉画像に適用することにより、輪郭適用後の反射干渉画像を生成する輪郭適用ステップと、第３抽出手段が、前記輪郭適用後の反射干渉画像より、第３パラメータを抽出する第３抽出ステップと、を更に備えても良い。【００１４】
この発明によれば、第１パラメータおよび第２パラメータに加え、第３パラメータが更に得られる。この第３パラメータは、定量位相画像と反射干渉画像との間で細胞の輪郭を合わせた上で得られるパラメータであるため、第１パラメータおよび第２パラメータとは異なる性質のパラメータであり、使用者はこの第３パラメータを更に得ることにより、細胞の状態を適切に判別・評価するためのより多くの情報を得ることとなる。【００１５】
また、定量位相画像から細胞の輪郭を抽出し、当該輪郭を反射干渉画像に適用することにより、輪郭適用後の反射干渉画像を生成する処理の中で使用者の介入はない。一方、従来の反射干渉顕微鏡を用いる方法（以下、「反射干渉法」という。）では、個々の細胞の輪郭を予め認識しておく必要があるため、反射干渉法は輪郭認識のための他の方法との組み合わせで用いられることが多い。上記の非特許文献１は、蛍光法と反射干渉法の組み合わせについて開示しているが、細胞の輪郭を自動的に抽出することが困難であったため、操作者が手描きにより細胞の輪郭を定めている。また、非特許文献２は、透過照明法と反射干渉法の組み合わせについて開示しているが、細胞の輪郭を人の手に寄らず自動的に定めることについての言及はない。更に、非特許文献３は、明視野法と反射干渉法の組み合わせについて開示しており、細胞輪郭抽出の自動化について一応の言及はしているものの、例えば、使用者が明視野画像を見ながら輪郭抽出のための最適な閾値を手動で設定する必要があることが記載されているなど、使用者のある程度の介入の必要性を認め、細胞の輪郭を人の手に寄らず自動的に定めるには至っていない。これに対して、本発明では、定量位相画像からの細胞輪郭抽出、当該輪郭の反射干渉画像への適用、輪郭適用後の反射干渉画像の生成までの処理が、輪郭抽出手段および輪郭適用手段により人の手に寄らず行われる。これらの処理が人の手に寄らず行われることにより作業効率は著しく向上し、それに伴う作業時間の顕著な短縮が達成される。【００１６】
また、本発明において、前記第１パラメータは、前記細胞が置かれている基板と、前記細胞との間の接着状態に基づく情報であっても良い。【００１７】
この発明によれば、細胞の基板への接着状態に基づく情報として、細胞の基板への接着面積、撮像範囲全体に対した接着面積の割合、細胞の基板への接着模様（パターン・２次元分布)等の情報が得られ、使用者はこれらの情報を用いて細胞の状態を適切に判別・評価することができる。【００１８】
また、本発明において、前記第２パラメータは、前記細胞の光学的厚さ、面積、体積、または前記細胞内での屈折率の変化に基づく情報であっても良い。【００１９】
この発明によれば、細胞の光学的厚さ、面積、体積、または細胞内での屈折率の変化に基づく情報が得られ、使用者はこれらの情報を用いて細胞の状態を適切に判別・評価することができる。【００２０】
また、本発明において、前記第３パラメータは、前記細胞の前記輪郭の中での、前記基板への前記接着状態に基づく情報であっても良い。【００２１】
この発明によれば、細胞の輪郭の中での、つまり細胞が占める空間の範囲内での基板への接着状態に基づく情報として、細胞の全体面積に対した接着面積の割合等の情報が得られ、使用者はこのような情報を用いて細胞の状態を適切に判別・評価することができる。【００２２】
また、本発明において、予め既知の種類または状態の前記細胞に対して抽出したパラメータをリファレンスデータとして記憶するリファレンス記憶手段と、前記リファレンスデータに基づき、未知の細胞の種類または状態を判別する解析手段と、を更に備えても良い。【００２３】
この発明によれば、リファレンスデータに基づき、未知の被検体の細胞の種類または状態を判別することができる。【００２４】
また、本発明において、未知の細胞に対して抽出したパラメータの中から所定のパラメータを選択し、当該選択された所定のパラメータを用いて、前記未知の細胞の種類または状態を判別する解析手段を更に備えても良い。【００２５】
この発明によれば、解析手段は、リファレンスデータによらずに、処理を行うことができる。この場合には、リファレンス記憶手段を設けなくて済むので、装置構成が簡潔になる。なお、解析手段は、特異な値を示すパラメータを所定のパラメータとして自動的に選択するように構成されていても良く、ユーザからの入力に基づき所定のパラメータを選択するように構成されていても良い。【００２６】
また、本発明において、前記第１抽出手段または前記第２抽出手段が３つ以上のパラメータを抽出した場合に、前記３つ以上のパラメータに対して主成分分析を行うことにより、未知の細胞の種類または状態を判別する解析手段を更に備えても良い。【００２７】
この発明によれば、抽出した多数のパラメータに対して主成分分析を行うことにより、当該多数のパラメータが細胞判別に適切に影響を及ぼうようにすることができる。主成分分析によれば、抽出した多数のパラメータのうち２つを選択する代わりに、当該多数のパラメータの全部または３つ以上の一部が適切に反映された２つの主成分により細胞判別を行うからである。【００２８】
また、本発明において、前記定量位相計測用光源として、低コヒーレント光を使用しても良い。【００２９】
この発明によれば、定量位相計測において、波長幅が広く干渉性の低い低コヒーレント光を使用する。これにより、光学系に由来する干渉ノイズを低減でき、より安定した計測をすることができる。【００３０】
また、本発明において、前記反射干渉計測用光源として、低コヒーレント光を使用しても良い。【００３１】
この発明によれば、反射干渉計測において、波長幅の広く干渉性の低い低コヒーレント光を使用している。広い波長幅の照明光を用いることにより、干渉が起きる距離を狭くすることができ、細胞の基板への接着面に限定した反射干渉画像を取り出すことが可能となる。【００３２】
また、本発明において、前記反射干渉計測用光源から放射され前記細胞から反射される光が集光される対物レンズと、前記対物レンズの前記反射干渉計測用光源側の開口絞りと共役な位置に、リング状のスリットと、を更に備えても良い。【００３３】
この発明によれば、リング状に開口したスリットを反射干渉計測用光源からの照明光が通り、照明光は対物レンズの中心は通らず周辺を通って照明され、角度のついた高いＮＡの光のみを使って細胞が照明されるため、細胞上部の溶液からの反射光の影響を低減することができる。また、リング状のスリットを用いることにより、対物レンズ内部の反射による背景光を低減することもできる。【００３４】
また、本発明において、前記細胞を収容する容器を更に備え、前記容器の前記細胞の接着面の反対側に反射防止コートが施されていても良い。【００３５】
この発明によれば、ドライ系の対物レンズを用いる場合であっても、高いコントラストで細胞接着面の反射干渉画像を得ることができる。【００３６】
また、本発明において、前記容器の上部に、前記細胞を浸す培養液に接触するような高さに観察窓を備えても良い。【００３７】
この発明によれば、観察窓を用いて培養液の高さを一定に保つことにより、定量位相計測を安定な条件下で行うことができる。【００３８】
また、本発明において、前記観察窓において、前記培養液に接触する面の反対側に反射防止コートが施されていても良い。【００３９】
この発明によれば、容器上面からの反射光の影響を抑えて、高いコントラストで細胞接着面の反射干渉画像を得ることができる。【００４０】
また、本発明においては、前記反射干渉計測用光源と前記反射干渉撮像手段との間および前記定量位相計測用光源と前記反射干渉撮像手段との間に位置し、波長により反射対透過の比率を異ならせることが可能なミラーを更に備えても良い。【００４１】
この発明によれば、波長により反射対透過の比率を異ならせることが可能なミラーを備えることにより、例えば定量位相計測用光源の波長と反射干渉計測用光源の波長とが異なる場合に、反射干渉計測用光源の波長に対してはハーフミラーとして機能するような反射対透過の比率と設定し、定量位相計測用光源の波長に合わせては透過率を高めるような構成をとることができる。このような構成により、反射干渉計測においてはハーフミラーとして機能しながらも、定量位相計測においては観察光のロスを低減させることができる。
【発明の効果】
【００４２】
本発明によれば、細胞の状態を適切に判別・評価するためのより多くの情報を得ることを可能な細胞観察装置および細胞観察方法を提供することができる。
【図面の簡単な説明】
【００４３】
【図１】細胞観察装置１の全体構成を示す概要図である。
【図２】処理部２０のハードウェア構成図である。
【図３】処理部２０の機能および動作を示すフローチャートである。
【図４】画像位置合わせに用いる標本３０を示す図である。
【図５】輪郭抽出処理および輪郭適用処理の一例を示す図である。
【図６】細胞から得られたパラメータの一例を示す図である。
【図７】抽出したパラメータを用いて作成した２成分散布図の一例である。
【図８】リファレンスデータを用いて細胞判別を行った一例を示す図である。
【図９】第３実施例による効果を示すための図である。
【図１０】第４実施例にかかる細胞観察装置１Ａの構成概要図である。
【図１１】第４実施例による効果を示すための図である。
【図１２】第５実施例における容器１０２Ａを示す図である。
【図１３】第５実施例による効果を示すための図である。
【図１４】第５実施例による効果を示すための図である。
【図１５】第５実施例による効果を示すための図である。
【図１６】第５実施例の変形例を示す図である。
【図１７】第６実施例における画像取得部１０Ａを示す図である。
【図１８】第７実施例における画像取得部１０Ｂを示す図である。
【図１９】第８実施例における画像取得部１０Ｃを示す図である。
【図２０】第９実施例における画像取得部１０Ｄを示す図である。
【図２１】第１０実施例における画像取得部１０Ｅを示す図である。
【図２２】第１１実施例における画像取得部１０Ｆを示す図である。
【図２３】第１２実施例における画像取得部１０Ｇを示す図である。
【図２４】第１３実施例における画像取得部１０Ｈを示す図である。
【図２５】第１４実施例における画像取得部１０Ｉを示す図である。
【図２６】第１５実施例においてユーザに表示される画面の一例を示す図である。
【図２７】第１６実施例において、リファレンスとして用いた各細胞の反射干渉画像および定量位相画像を示す図である。
【図２８】第１６実施例において、抽出した７つのパラメータの実際の計測値を示す図である。
【図２９】第１６実施例において、パラメータの情報量として母集団に対した偏差平方和を求めた結果を示す図である。
【図３０】第１６実施例において、各パラメータに対して標準化を行った結果を示す図である。
【図３１】第１６実施例において、第１主成分ｆの係数を求めた結果を示す図である。
【図３２】第１６実施例において、第２主成分ｇの係数を求めた結果を示す図である。
【図３３】第１６実施例において、各細胞について第１主成分ｆおよび第２主成分ｇを求めた結果を示す図である。
【図３４】第１６実施例において、第１主成分ｆおよび第２主成分ｇの値を用いて、２成分散布図に各細胞をプロットした結果を示す図である。
【発明を実施するための形態】
【００４４】
以下、添付図面を参照して本発明にかかる細胞観察装置および細胞観察方法の好適な実施形態を詳細に説明する。まず、細胞観察装置および細胞観察方法の概要について全体的に説明し、続いて、第１実施例〜第１６実施例の欄を用いてより詳細に説明する。なお、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。
【００４５】
［細胞観察装置１の全体構成］
まず、本発明の実施形態に係る細胞観察装置１の全体構成について、図１を参照しながら説明する。図１は、細胞観察装置１の全体構成を示す概要図である。図１に示すように、細胞観察装置１は、画像取得部１０および処理部２０から構成されている。
【００４６】
画像取得部１０は、試料となる細胞１０１を収容する容器１０２、細胞１０１の培養環境が維持された培養空間１０３、対物レンズ１０４、定量位相計測用光源１０５、照明光絞り部１０５Ａ、反射干渉計測用光源１０６、ハーフミラー１０７、ダイクロイックミラー１０８、全反射ミラー１０９、反射干渉検出用カメラ１１０（特許請求の範囲の「反射干渉撮像手段」に相当）、回折型干渉光学系１１１、および定量位相検出用カメラ１１２（特許請求の範囲の「定量位相撮像手段」に相当）を備える。ハーフミラー１０７は反射干渉計測用光源１０６からの光を細胞１０１へ導くための反射干渉入射用光学系であり、また、ダイクロイックミラー１０８は、反射干渉検出用カメラ１１０へ細胞１０１からの光を導くための反射干渉計測用光学系である。また、対物レンズ１０４およびハーフミラー１０７、ダイクロイックミラー１０８は、細胞１０１の同一範囲からの光を反射干渉検出用カメラ１１０および定量位相検出用カメラ１１２へ導くための共通光学系である。つまり、本実施形態では、反射干渉入射用光学系と反射干渉計測用光学系が共通光学系となっている。また、全反射ミラー１０９および回折型干渉光学系１１１は、定量位相検出用カメラ１１２へ光を導くための定量位相光学系である。反射干渉検出用カメラ１１０は、反射干渉計測用光源１０６から放射され、細胞１０１から反射される光を撮像することにより、反射干渉画像を生成するものである（特許請求の範囲の「反射干渉撮像ステップ」に相当）。定量位相検出用カメラ１１２は、定量位相計測用光源１０５から放射され照射光が、ピンホールやアパーチャーなどの照明光絞り部１０５Ｂを通過することで点光源に近い照射光となり、細胞１０１を透過する光を撮像することにより、定量位相画像を生成するものである（特許請求の範囲の「定量位相撮像ステップ」に相当）。なお、定量位相計測用光源１０５としては、ハロゲンランプやキセノンランプなどの広い波長範囲にわたって放射感度を有する光源を用いることができる。また、ＬＥＤ（発光ダイオード）や半導体レーザー（レーザーダイオード）、ＳＬＤ（スーパールミネッセントダイオード）などの光源を用いてもよい。レーザーやＳＬＤの場合には、光源サイズが十分小さいため、照明光絞り部１０５Ａは不要になる。
【００４７】
処理部２０は、画像位置合わせ部２０１（特許請求の範囲の「画像位置合わせ手段」に相当）、第１抽出部２０４（特許請求の範囲の「第１抽出手段」に相当）、第２抽出部２０５（特許請求の範囲の「第２抽出手段」に相当）からなる第１処理部、輪郭抽出部２０２（特許請求の範囲の「輪郭抽出手段」に相当）、輪郭適用部２０３（特許請求の範囲の「輪郭適用手段」に相当）、第３抽出部２０６（特許請求の範囲の「第３抽出手段」に相当）からなる第２処理部、解析部２０７（特許請求の範囲の「解析手段」に相当）、およびリファレンス記憶部２０８（特許請求の範囲の「リファレンス記憶手段」に相当）からなる第３処理部を備える。なお、第１処理部および第２処理部、第３処理部は、それぞれ別の演算装置であってもよいし、同一の演算装置であってもよい。
【００４８】
画像位置合わせ部２０１は、反射干渉画像と定量位相画像との間で、空間的な位置を一致させることにより、両画像の位置合わせを行うものである（特許請求の範囲の「画像位置合わせステップ」に相当）。輪郭抽出部２０２は、定量位相画像に基づき、細胞１０１の範囲である輪郭を抽出するものである（特許請求の範囲の「輪郭抽出ステップ」に相当）。輪郭適用部２０３は、輪郭抽出手段が抽出した輪郭を反射干渉画像に適用することにより、輪郭適用後の反射干渉画像を生成するものである（特許請求の範囲の「輪郭適用ステップ」に相当）。
【００４９】
第１抽出部２０４は、画像位置合わせ手段により定量位相画像との間で位置合わせが行われた反射干渉画像より、第１パラメータを抽出するものである（特許請求の範囲の「第１抽出ステップ」に相当）。第１パラメータは、細胞１０１が置かれている基板と、細胞１０１との間の接着状態に基づく情報である。細胞１０１の基板への接着状態に基づく情報としては、例えば、細胞１０１の基板への接着面積、撮像範囲全体に対した接着面積の割合、細胞の基板への接着模様（パターン・２次元分布）等の情報が挙げられる。
【００５０】
第２抽出部２０５は、画像位置合わせ手段により反射干渉画像との間で位置合わせが行われた定量位相画像より、第２パラメータを抽出するものである（特許請求の範囲の「第２抽出ステップ」に相当）。第２パラメータは、細胞１０１の光学的厚さ、面積、体積、または細胞１０１内での屈折率の変化に基づく（または表す）情報である。
【００５１】
第３抽出部２０６は、輪郭適用後の反射干渉画像より、第３パラメータを抽出するものである（特許請求の範囲の「第３抽出ステップ」に相当）。第３パラメータは、細胞１０１の輪郭の中での、基板への接着状態に基づく情報である。細胞１０１が占める空間の範囲内での基板への接着状態に基づく情報としては、例えば、細胞１０１の全体面積に対した接着面積の割合等の情報が挙げられる。
【００５２】
リファレンス記憶部２０８は、予め既知の種類または状態の細胞に対して抽出したパラメータ（上記の第１、第２、第３パラメータ）をリファレンスデータとして記憶するものである。解析部２０７は、第１抽出部２０４、第２抽出部２０５、第３抽出部２０６で抽出された第１、第２、第３パラメータを基にして、未知の被検体の細胞の種類または状態を判別するものである。具体的には、第１、第２、第３パラメータのうち、所定のパラメータを選択し、選択された所定のパラメータを用いて、未知の被検体の細胞の種類または状態を判別する。所定のパラメータは、解析部２０７により（または図示しない別途の処理手段により）自動的に選択された特異な値を示すパラメータであってもよいし、ユーザーが選択したパラメータであってもよい。解析部２０７は、このように選択された所定のパラメータを用いて、未知の被検体の細胞の種類または状態を判別する。また、解析部２０７が、リファレンス記憶部２０８に記憶されたリファレンスデータに基づき、所定のパラメータを選択し、未知の被検体の細胞の種類または状態を判別しても良い。更に、第１抽出部２０４、第２抽出部２０５、および第３抽出部２０６のそれぞれ、または組み合わせにより、３つ以上のパラメータが抽出された場合には、解析部２０７が当該３つ以上のパラメータに対して主成分分析を行うことにより、未知の細胞の種類または状態を判別しても良い。なお、以下の記載においては、解析部２０７がリファレンスデータに基づき所定のパラメータを選択し、未知の被検体の細胞の種類または状態を判別する態様を［第１実施例］の欄で記載する。一方で、解析部２０７が、リファレンスデータに基づかないで、特異な値またはユーザ指示により所定のパラメータを選択し、未知の被検体の細胞の種類または状態を判別する態様を［第１５実施例］の欄で記載する。更に、解析部２０７が主成分分析により未知の細胞の種類または状態を判別する態様を［第１６実施例］の欄で記載する。
【００５３】
図２は、以上のような機能的構成要素を備える処理部２０のハードウェア構成図である。図２に示すように、処理部２０は、物理的には、ＣＰＵ２１、ＲＯＭ２２及びＲＡＭ２３等の主記憶装置、キーボード及びマウス等の入力デバイス２４、ディスプレイ等の出力デバイス２５、画像取得部１０との間でデータの送受信を行うためのネットワークカード等の通信モジュール２６、ハードディスク等の補助記憶装置２７などを含む通常のコンピュータシステムとして構成される。処理部２０の各機能は、ＣＰＵ２１、ＲＯＭ２２、ＲＡＭ２３等のハードウェア上に所定のコンピュータソフトウェアを読み込ませることにより、ＣＰＵ２１の制御の元で入力デバイス２４、出力デバイス２５、通信モジュール２６を動作させると共に、主記憶装置２２，２３や補助記憶装置２７におけるデータの読み出し及び書き込みを行うことで実現される。
【００５４】
［第１実施例］
引き続き、本発明の第１実施例について図１を再び参照しながら詳細に説明する。図１は、第１実施例にかかる細胞観察装置１の構成概要図である。
【００５５】
（画像取得部１０の説明）
図１に示されるように、測定対象となる細胞１０１の入った容器１０２を細胞１０１の培養環境が維持された培養空間１０３に静置する。培養環境は、細胞１０１の生育、状態維持に適した温度、湿度、炭酸ガス濃度などが整った環境である。
【００５６】
定量位相計測について述べる。細胞１０１の入った容器１０２の上部に設置した定量位相計測用光源１０５から放射された照明光は、細胞１０１の入った容器１０２を透過し、対物レンズ１０４にて集光される。そして、ハーフミラー１０７を介し、定量位相像と反射干渉像を波長で分離するためのダイクロイックミラー１０８を透過して、更に全反射ミラー１０９を介して、位相計測のための回折型干渉光学系１１１にて物体光と参照光の干渉像を形成し、定量位相検出用カメラ１１２にて干渉縞画像が撮像される。
【００５７】
反射干渉計測について述べる。反射干渉計測用光源１０６から放射された照明光は、ハーフミラー１０７にて反射して、対物レンズ１０４を通って、測定対象となる細胞１０１を入れた容器１０２の底面側から入射する。ハーフミラー１０７に限らず、照明光強度が十分ある場合には５：９５（反射：透過）や２０：８０（反射：透過）にするなど、反射比率を低減させたビームスプリッタを用いてもよい。また、波長によって反射率・透過率が異なるダイクロイックミラーを用いてもよい。容器１０２の底面の細胞１０１の接着面から得られる反射光は細胞１０１の接着距離に応じて干渉を生じ、得られた反射干渉光は再び対物レンズ１０４にて集光され、ハーフミラー１０７を介し、定量位相像と反射干渉像を波長で分離するためのダイクロイックミラー１０８にて反射干渉像のみを反射して、反射干渉検出用カメラ１１０にて撮像される。なお、ダイクロイックミラー１０８の代わりにハーフミラーなどのビームスプリッタを用いる場合は、ビームスプリッタと反射干渉検出用カメラ１１０との間に反射干渉検出に用いる波長を選択するためのフィルターを配置すればよい。容器１０２の底面の細胞１０１の接着面から反射した光は、細胞１０１の接着距離に応じて干渉した光の振幅が異なり、明暗のコントラストとして撮影される。なお、定量位相像および反射干渉像は、共通の対物レンズ１０４を介して取得されるので、細胞１０１の撮像範囲は、定量位相計測および反射干渉計測においてほぼ同一範囲となる。
【００５８】
定量位相画像を取得するカメラ１１２と反射干渉画像を取得するカメラ１１０は同じ性能、同じ画素分解能のカメラに限定する必要はない。それぞれのカメラに入射する光量や入射波長は互いに異なるので、それぞれに適した性能、空間分解能のカメラを用いてもよい。例えば定量位相検出用カメラ１１２には８３０ｎｍに感度が高く、画素サイズが大きめの感度優先のカメラを用い、反射干渉検出用カメラ１１０には可視域に感度が高く、画素サイズが小さめの空間分解能優先のカメラを利用するなどしてよい。このような場合に二つのカメラの空間座標を一致させる処理が必要となるが、これについては後述する。
【００５９】
定量位相計測と反射干渉計測を波長で区別できるように、反射干渉計測用光源１０６と定量位相計測用光源１０５のそれぞれの照明波長には、異なる波長域の光を用いても良い。ダイクロイックミラー１０８にて特定の波長で定量位相像と反射干渉像を分離選択することができる。したがって、定量位相計測と反射干渉計測で、クロストークなく同時刻の像を得ることができる。
【００６０】
できるだけ多くの細胞１０１を見るために、観察する位置を移動させる機構を有しても良い。細胞１０１への影響を少なくするため、および定量位相計測において液面の揺れを抑えるために、細胞１０１を入れた容器は静止したままで、照明光学系と観察光学系を一体とした画像取得部１０の本体をＸＹ平面上で移動させることによって観察位置を変える方法が望ましい。これとあわせて、観察しているＸＹの空間平面座標を画像に記録することが望ましい。
【００６１】
（処理部２０の説明）
図３のフローチャートを更に参照しながら、処理部２０の機能および動作について説明する。
【００６２】
最初に、定量位相検出用カメラ１１２にて参照光と細胞１０１を透過した物体光の干渉縞画像が得られる（ステップＳ１０１、特許請求の範囲の「定量位相撮像ステップ」に相当）。干渉縞画像から公知の演算方法を用いて定量的な位相画像を形成する。定量位相画像を得るためには、細胞１０１のない背景領域のオフセット補正および背景領域の視野におけるシェーディング補正を施し、背景部分は空間的に均一で、背景部の位相の値を０に補正し、細胞１０１の位相（光路長）の２次元マップを得る。
【００６３】
一方で、ステップＳ１０１とともに、反射干渉検出用カメラ１１０にて細胞１０１の接着面の反射干渉画像が得られる（ステップＳ１０２、特許請求の範囲の「反射干渉撮像ステップ」に相当）。容器１０２の底面に接着している細胞１０１の容器１０２の底面からの距離に応じて干渉光の振幅が異なり、反射干渉画像は明暗のコントラストとして撮影される。反射干渉画像の視野内の反射光のシェーディングを補正する。あわせて細胞１０１のない背景の値が時間的に変動しないように、時間ごとに背景部のオフセット補正を行う。これらの画像演算補正によって、空間的、時間的に変動の少ない定量位相画像および反射干渉画像を得ることができる。
【００６４】
次に、定量位相画像および反射干渉画像の空間的位置補正が行われる（ステップＳ１０３、特許請求の範囲の「画像位置合わせステップ」に相当）。図４は、画像位置合わせに用いる標本３０を示す。位置合わせの標本３０として、ガラス基板に格子３１が刻まれたマイクロスケールやドットパターンなどを用いる。マイクロスケールやドットパターンには２つのカメラで撮影したときに同じ位置が分かるように、少なくとも１か所以上に目印が付いていると操作が簡便である。図４の例では、三角と丸の二つの目印が付いている。
【００６５】
ここで、反射干渉側のカメラ１１０の像を定量位相の像に合わせて位置をずらす手順について詳細に説明する。まず、細胞１０１の標本の代わりにマイクロスケールまたはドットパターンを置き、反射干渉および定量位相のそれぞれの撮影用カメラで撮影する。次に、マイクロスケールまたはドットパターンを撮影して得られた反射干渉画像と定量位相画像をそれぞれ異なる疑似カラー（例えば反射干渉画像は緑、定量位相画像は赤など）で重ね合わせる。場合によっては画像の輝度を反転すると見やすくなる。続いて、反射干渉画像の格子像またはドット像を定量位相側の格子像またはドット像に合致させるために、重ね合わせた画像を見ながら、反射干渉画像に対して（１）拡大縮小倍率、（２）水平方向の移動量（画素数）、（３）垂直方向の移動量（画素数）、（４）回転角度、（５）回転方向、（６）左右反転などを微調整し、調整量を決定する。調整量は、重ね合わせて表示した格子像またはドット像が重なり、エッジが出なくなるような位置に決定されればよく、例えば、それぞれの格子像またはドット像が重なり合って緑色と赤色が黄色になるような位置で決定される。
【００６６】
撮影した２つのカメラの像に空間歪があり、両者の空間歪のパターンが大きく異なる場合は、上述の（１）から（６）のデータだけで一律にすべての画素の位置を合致させることは難しい場合がある。その場合は視野内のすべての画素における移動量が異なるため、各々の画素に対して移動量を定める必要がある。具体的には、格子像の各交点またはドット像の各ドットの中心座標が２つのカメラで一致するように、各点のずらし量を決定する。格子またはドットの存在する座標以外の空間座標は補間することによってずらし量を決定する。このようにして空間的に全画素（全座標）について１画素単位でずらし量を記憶したテーブルを作成して、位置合わせデータとする。得られた位置合わせデータをファイルで保存する。新規に測定を行い画像を取得したときに、ファイルで保存した位置合わせデータを用いて、画像間の位置合わせを行い、位置を補正した画像を出力する。また、新規に測定を行い取得した画像データをファイルに保存するときには、位置合わせデータを画像データと一緒にファイルに保存することが好ましい。こうすることによって、画像ファイルを呼び出した時にも、画像を新規に取得した時の最も適切であった位置合わせデータを援用して、当該呼び出した画像を補正することができる。
【００６７】
なお、以上のような位置合わせは定量位相や反射干渉の計測のたびに行う必要はなく、同じ光学系を利用するのであれば、使用環境や振動の影響による位置ずれを考慮して、例えば１カ月に１回程度の頻度で行えば良い。また、観察側の光学系に観察用のフィルタを選択する部品（例えば電動フィルタホイールなど）が装着されている場合は、位置合わせは使用するフィルタごとに行うことが望ましい。なぜなら、使用するフィルタの傾きや平行度によって像の位置ずれ量や位置ずれの方向が異なるからである。
【００６８】
なお、繰り返しとなるため説明は省略するが、以上とは逆に、定量位相側のカメラ１１２の像を反射干渉の像に合わせて位置をずらす場合にも、上述の手順を適宜適用することができる。また、以上の手順は、自動化された画像処理によって、人の手に寄らずに、行うことができる。
【００６９】
図３に戻り、ステップＳ１０３の位置合わせの後に、細胞１０１の輪郭領域を抽出する処理（以下、「セグメンテーション」ともいう。）が行われる（ステップＳ１０４およびＳ１０５、特許請求の範囲の「輪郭抽出ステップ」および「輪郭適用手段」に相当）。
【００７０】
まず、図５に示すように、同時刻に撮像され、位置合わせが行われた後の定量位相画像および反射干渉画像のうち、定量位相画像から個々の細胞１０１の輪郭となる領域を画像処理により検出する（ステップＳ１０４、図５の画像Ａ、特許請求の範囲の「輪郭抽出ステップ」に相当）。つまり、定量位相画像では、細胞１０１のない背景となる溶液を通る光に対して細胞１０１を通る光は、細胞１０１の屈折率が溶液の屈折率よりも大きいために、光路長が長くなる。そのため、細胞１０１の存在する領域の画素の位相の値は背景よりも大きな値となる。そこで、適正な閾値または空間フィルタ処理を用いることにより、人の手に寄らず自動的に、背景と細胞１０１を分離することができる。そして、一つ一つの細胞１０１に対応した輪郭を決定し、それぞれの細胞１０１の占める領域に対応する画素座標の領域を決定することができる。
【００７１】
次に、ステップＳ１０４で得られた個々の細胞１０１の輪郭領域の画素座標を空間座標の一致した反射干渉画像に適応させ、つまりステップＳ１０４で得られたセグメンテーション領域を反射干渉画像にコピーすることにより、定量位相画像で決定された個々の細胞１０１の輪郭領域を反射干渉画像に当てはめる（ステップＳ１０５、図５の画像Ｂ、特許請求の範囲の「輪郭適用ステップ」に相当）。このようにすることで、図５に示すように、定量位相画像と反射干渉画像の両方の画像Ａ，Ｂに対して同じ輪郭領域を決定することができる。
【００７２】
次に、ステップＳ１０４およびＳ１０５のセグメンテーションの後に、パラメータを抽出する処理が行われる（ステップＳ１０６、特許請求の範囲の「第１抽出ステップ」、「第２抽出ステップ」および「第３抽出ステップ」に相当）。
【００７３】
すなわち、ステップＳ１０４で得られた一つ一つの細胞１０１の領域において、第２抽出部２０５が、定量位相画像から（１）細胞１０１の光学的厚さ（光路長）の平均値、（２）光学的厚さの標準偏差、（３）細胞１０１の面積、（４）細胞１０１の光学的体積（光学的厚さの合計値）を求める（ステップＳ１０６、特許請求の範囲の「第２抽出ステップ」に相当）。更に、ステップＳ１０３で位置合わせした反射干渉画像、およびステップＳ１０５でセグメンテーション領域をコピーした反射干渉画像から、第１抽出部２０４および第３抽出部２０６が、それぞれ、（５）明暗の平均輝度、（６）強い接着を示す暗領域の面積、（７）弱い接着を示す明領域の面積、（８）撮像範囲全体および細胞１０１の面積に占める暗領域の割合、（９）撮像範囲全体および細胞１０１の面積に占める明領域の割合を求める（ステップＳ１０６、特許請求の範囲の「第１抽出ステップ」および「第３抽出ステップ」に相当）。
【００７４】
また、接着部位の２次元分布パターンを求める方法として、反射干渉画像の明暗の２次元濃淡分布パターンから求めたテクスチャ解析の定量解析を行う。テクスチャ解析には、画像における濃度ｉの画素がθ方向に距離ｄ画素離れた画素の濃度がｊである確率Ｐ（ｉ，ｊ）（ｉ，ｊ＝０，１，２，３，．．．ｎ−１）を要素とする同時生起行列から求めたテクスチャの特徴を得る同時生起行列法、画像における濃度ヒストグラムＰ（ｉ）からテクスチャの特徴を求める濃度ヒストグラム法等を適用する。この解析処理は解析部２０７により行われても良い。
【００７５】
同時生起行列法から得られるテクスチャの特徴パラメータとしては、（１０）エネルギー、（１１）エントロピー、（１２）相関、（１３）局所一様性、（１４）慣性等が挙げられる。この中で、エネルギーは同時生起行列の確率分布が特定の値に集中しているかを示す指標となり、エントロピーは同時生起行列の確率分布が広い範囲の値に分布しているかを示す指標となる。一方、濃度ヒストグラム法から得られるテクスチャの特徴パラメータとしては、（１５）平均、（１６）分散、（１７）歪度、（１８）尖度等が挙げられる。この中で、歪度は濃度ヒストグラムの形状が対称な形からどれだけずれているかを示す指標となり、尖度は濃度ヒストグラムの分布がどれだけ平均値の周りに集中しているのかを示す指標となる。
【００７６】
このうち、同時生起行列法の（１２）相関や濃度ヒストグラム法の（１５）平均、（１６）分散のパラメータは画像取得時の輝度の明るさの影響を受けやすいので、評価する場合には同じ露光時間や同じ照明輝度で取得したものに限るなどの条件を把握していた方が良い。この場合、分散のかわりに（１９）「√（分散）／平均輝度」を用いることにより、画像取得時の明るさの条件に影響を受けないパラメータとして利用することができる。また、同時生起行列法では、θ方向、距離ｄについて複数の値を得ることができるが、あらかじめどの方向、どの距離のパラメータが細胞間で差が出るかを予備実験で検討および決定しておくと、膨大なパラメータの数を最適化することができる。今回はθ方向＝０度、距離ｄ＝８画素に固定して結果を得ることにする。このようにして一つ一つの細胞１０１について（１）から（１９）を計測結果として数値として記憶する。
【００７７】
図６に個々の細胞から得られたパラメータの一例を示す。またこれらのパラメータのうちの２つの成分を用いて散布図を描くことができる。図７は骨髄幹細胞（stem cells）と分化誘導した破骨細胞（osteoclasts）の２種類の細胞について、識別を行うために特徴的なパラメータを用いてプロットしたものである。図７（Ａ）は、横軸に反射干渉の明暗の平均輝度（５）をとり、縦軸に定量位相の細胞の光学的厚さの平均値（１）をとって、それぞれの細胞をプロットした２成分散布図である。一方、図７（Ｂ）は横軸に反射干渉のテクスチャ解析の濃度ヒストグラム法の尖度（１６）をとり、縦軸に反射干渉のテクスチャ解析の同時生生起行列法の慣性（１３）をとってプロットしたものである。何れの２成分散布図においても、骨髄幹細胞と破骨細胞はそれぞれ集団を形成しており、種類の異なる細胞を識別する有効なパラメータとなっている。
【００７８】
図３に戻り、ステップＳ１０６にて個々の細胞１０１から求めた計測情報（パラメータ）をもとに、未知の被検体の細胞がどの種類の細胞集団に類似しているのか、またはどんな状態の細胞集団に類似しているのかを判別する（ステップＳ１０７）。そのためには、あらかじめ既知の状態の細胞集団の計測情報をリファレンスデータとして格納しておき、被検体細胞とリファレンスデータとを比較することにより被検体細胞の種類や状態を判別する。細胞観察装置１には、図１に示すように、リファレンスデータを格納するためのリファレンス記憶部２０８が設けられている。リファレンス記憶部２０８には、複数種類の細胞または複数種類の状態の細胞について記憶することができ、初期値として代表的な細胞群についてリファレンスを持つこともできるし、作業者にとって必要な細胞群について新たに計測を行い追加することもできる。
【００７９】
リファレンスデータを取得する手順としては、最初に、あらかじめ既知の種類の細胞集団または既知の状態の細胞集団について複数個の細胞を用意する。これらは細胞間のばらつきが少なく、できるだけ均質なものを用意する。次に、これらの細胞に対して上述したステップＳ１０１〜Ｓ１０６の処理を画像取得部１０および処理部２０の協同により行い、求めた計測結果をリファレンスデータとしてリファレンス記憶部２０８に記憶する。
【００８０】
リファレンスデータを用いた細胞判別の一例として、被検体となる未知の細胞が二つの細胞集団のうちどちらに属するかを判別する場合を考える。なお、以下の判別処理は、解析部２０７によって行われる。最初に、リファレンスとして格納した２種類の細胞集団について、個々の細胞の（１）から（１９）の計測結果のうち、任意の２成分を選んでそれぞれを縦軸、横軸にとった２成分散布図を作成する。作成した２成分散布図の中で二つの細胞集団を単純な区画線で簡便に識別できるような２成分散布図を、被検体の細胞を識別する際の判断基準として利用する。例えば図８（Ａ）に示す一例のリファレンスデータでは、細胞Ａと細胞Ｂを識別する手段として、横軸に反射干渉のテクスチャ解析の濃度ヒストグラム法の尖度（１６）をとり、縦軸に反射干渉のテクスチャ解析の同時生生起行列法の慣性（１３）をとってプロットした。この２成分の散布図を用いると、細胞Ａと細胞Ｂを単純な区画線Ｌにて識別することができる。この区画線Ｌを用いた細胞Ａの識別の性能としては、感度１００％、特異度１００％で識別が可能である。
【００８１】
次に、被検体となる未知の細胞について計測を行い、未知の細胞一つ一つについて上述の（１）から（１９）の計測結果のパラメータを得る。被検体の細胞が細胞Ａか細胞Ｂであるかを識別したい場合、細胞Ａと細胞Ｂを識別するためのリファレンスデータの２成分散布図を呼び出し、これを参照する。図８（Ｂ）および（Ｃ）のように、リファレンスデータの２成分散布図に被検体の細胞をプロットし、区画線Ｌを用いて細胞Ａおよび細胞Ｂに相当する領域を明示する。すると、被検体の細胞は区画線Ｌにて区切られた細胞Ａの領域または細胞Ｂの領域に存在し、どちらの細胞集団に属しているかを判別することができる。図８（Ｃ）においては、１３個の被検体の細胞のうち、７個の細胞は細胞Ａに分類され、残る６個の細胞は細胞Ｂに分類されるという結果が得られる。もちろん、図８（Ｄ）および（Ｅ）に示すように、２成分散布図にプロットされた一つ一つのデータが、定量位相画像および反射干渉画像においてどの細胞であるかを関連付けて明示することもできる。
【００８２】
なお、区画線は２成分散布図において目視により縦軸横軸に水平な矩形領域を設けてもよいし、統計解析の判別分析の手法を用いて２集団を弁別する線形判別式を用いてもよい。線形判別式にて細胞を識別する場合には、個々の細胞について、得られた線形判別式をもとに判別得点を求めると、２群における判別得点の一方は正の値、他方は負の値をとるので、簡単に細胞を識別することができる。
【００８３】
（第１実施例の作用及び効果）
続いて、以上まで説明した第１実施例にかかる細胞観察装置１の作用及び効果について説明する。本実施例の細胞観察装置１によれば、反射干渉計測用光源１０６、反射干渉検出用カメラ１１０、および第１抽出部２０４の反射干渉計測ユニットを備えることにより、細胞からの反射光をもとに、第１パラメータを得る。また、定量位相計測用光源１０５、定量位相検出用カメラ１１２、および第２抽出部２０５の定量位相計測ユニットを備えることにより、細胞からの透過光をもとに、第２パラメータを得る。このように、本実施例の細胞観察装置１では、反射干渉計測ユニットおよび定量位相計測ユニットの両方を備えることにより、第１パラメータおよび第２パラメータの両方を得ることができ、これにより使用者は細胞の状態を適切に判別・評価するためのより多くの情報を得ることとなる。更に、第１パラメータは定量位相画像との間で位置合わせが行われた反射干渉画像より抽出されるものであり、第２パラメータは反射干渉画像との間で位置合わせが行われた定量位相画像より抽出されるものであるため、両パラメータ間には細胞のどの部分に関するパラメータであるかについての整合性が取れているといえる。つまり、本実施例の画像位置合わせ部２０１による両画像の位置合わせにより、両パラメータは細胞の同一部分の状態を判別・評価するためのパラメータとして活用可能となる。
【００８４】
また、本実施例によれば、第１パラメータおよび第２パラメータに加え、第３パラメータが更に得られる。この第３パラメータは、定量位相画像と反射干渉画像との間で細胞の輪郭を合わせた上で得られるパラメータであるため、第１パラメータおよび第２パラメータとは異なる性質のパラメータであり、使用者はこの第３パラメータを更に得ることにより、細胞の状態を適切に判別・評価するためのより多くの情報を得ることとなる。
【００８５】
また、定量位相画像から細胞の輪郭を抽出し、当該輪郭を反射干渉画像に適用することにより、輪郭適用後の反射干渉画像を生成する処理の中で使用者の介入はない。一方、従来の反射干渉顕微鏡を用いる方法（反射干渉法）では、個々の細胞の輪郭を予め認識しておく必要があるため、反射干渉法は輪郭認識のための他の方法との組み合わせで用いられることが多い。上述したように、上記の非特許文献１は、蛍光法と反射干渉法の組み合わせについて開示しているが、細胞の輪郭を自動的に抽出することが困難であったため、操作者が手描きにより細胞の輪郭を定めている。また、非特許文献２は、透過照明法と反射干渉法の組み合わせについて開示しているが、細胞の輪郭を人の手に寄らず自動的に定めることについての言及はない。更に、非特許文献３は、明視野法と反射干渉法の組み合わせについて開示しており、細胞輪郭抽出の自動化について一応の言及はしているものの、例えば、使用者が明視野画像を見ながら輪郭抽出のための最適な閾値を手動で設定する必要があることが記載されているなど、使用者のある程度の介入の必要性を認め、細胞の輪郭を人の手に寄らず自動的に定めるには至っていない。これに対して、本実施例では、定量位相画像からの細胞輪郭抽出、当該輪郭の反射干渉画像への適用、輪郭適用後の反射干渉画像の生成までの処理が、輪郭抽出部２０２および輪郭適用部２０３により人の手に寄らず行われる。これらの処理が人の手に寄らず行われることにより作業効率は著しく向上し、それに伴う作業時間の顕著な短縮が達成される。
【００８６】
輪郭抽出部２０２および輪郭適用部２０３によるセグメンテーションの自動化により、どれぐらい作業時間が短縮できるかについて、非特許文献３の場合と比べる実験を実施した。視野中に１００個ほど存在する細胞を明視野（非特許文献３の場合）と定量位相（本実施例の場合）にてそれぞれ撮影し、セグメンテーションの処理にかかる時間を計測した。明視野の場合、明るさの閾値を用いたセグメンテーションが困難であったため、マウスにて細胞一つ一つの輪郭を通常の作業能力を有する作業者の手によりトレースしてセグメンテーションを行った。この場合、１００個の細胞をセグメンテーションするのに１５０秒を要した。一方、定量位相にて得られた細胞の位相画像においては、輪郭抽出部２０２および輪郭適用部２０３における自動閾値と画像処理によってセグメンテーションを行った。この場合には、３秒の演算処理時間で１０８個の細胞をセグメンテーションすることができ、１５０秒から３秒へと大幅な時間短縮が可能であった。以上の実験結果から分かるように、本実施例の輪郭抽出部２０２および輪郭適用部２０３を備え、人の手に寄らずセグメンテーションを行うことにより、作業効率は著しく向上し、それに伴う作業時間の顕著な短縮が達成される。
【００８７】
また、本実施例によれば、細胞の基板への接着状態に基づく情報として、細胞の基板への接着面積、撮像範囲全体に対した接着面積の割合、細胞の基板への接着模様（パターン・２次元分布）等の情報が得られ、使用者はこれらの情報を用いて細胞の状態を適切に判別・評価することができる。
【００８８】
また、本実施例によれば、細胞の光学的厚さ、面積、体積、または細胞内での屈折率の変化に基づく情報が得られ、使用者はこれらの情報を用いて細胞の状態を適切に判別・評価することができる。
【００８９】
また、本実施例によれば、細胞の輪郭の中での、つまり細胞が占める空間の範囲内での基板への接着状態に基づく情報として、細胞の全体面積に対した接着面積の割合等の情報が得られ、使用者はこのような情報を用いて細胞の状態を適切に判別・評価することができる。
【００９０】
また、本実施例によれば、リファレンスデータに基づき、未知の被検体の細胞の種類または状態を判別することができる。
【００９１】
［第２実施例］
引き続き、本発明の第２実施例について説明する。第２実施例では、図１に示した第１実施例の構成要素を全て含み、定量位相計測用光源１０５に更なる特徴がある。
【００９２】
従来、位相変化を求めるために干渉性の高いレーザ光が光源としてよく使用されていた。しかしながら、レーザ光を使用する場合には、対物レンズなどの光学系に由来する過干渉による背景ノイズや、細胞内の高屈折率の顆粒などからの散乱光ノイズが、計測に影響を与えることがしばしばある。本実施例の細胞観察装置１の定量位相計測においては、照明用光源として、波長幅が広く干渉性の低い低コヒーレント光を使用することに特徴がある。この特徴により、過干渉によるノイズを低減でき、より安定した計測をすることができる。一例では、中心波長８３０ｎｍ、波長幅２０ｎｍ程度のＳＬＤ（スーパールミネセンスダイオード）を使用することにより、過干渉によるノイズを低減でき、より安定した計測に成功した。
【００９３】
［第３実施例］
引き続き、本発明の第３実施例について説明する。第３実施例では、図１に示した第１実施例の構成要素を全て含み、反射干渉計測用光源１０６に更なる特徴がある。
【００９４】
従来、コントラストを得るためにある程度波長幅を制限したバンドパスフィルタを用いた照明光が利用されていた。しかしながら、バンドパスされた照明光では、干渉性が高く、細胞の接着面とは関係のない細胞の上部の細胞膜と培養液との界面からの反射による干渉縞が映り込むことがしばしばある。本実施例の細胞観察装置１の反射干渉計測においては、波長幅の広く干渉性の低い低コヒーレント光を使用している。広い波長幅の照明光を用いることにより、干渉が起きる距離を狭くすることができ、細胞の基板への接着面に限定した反射干渉画像を取り出すことが可能となる。
【００９５】
図９に４２０ｎｍ〜７５０ｎｍの波長幅の広い照明光を用いた場合と４８０ｎｍ、５３０ｎｍを中心波長として波長幅３０ｎｍ程度に狭い波長範囲でバンドパスされた照明光を用いた場合の反射干渉画像の違いを示す。図９（Ａ）と図９（Ｂ）のように狭い波長範囲でバンドパスされた照明光では、細胞の接着面ではない細胞の上部の細胞膜と培養液との界面からの反射が干渉縞のように重なって観察されるが、図９（Ｃ）のように波長幅の広い照明光を用いることにより、細胞の上部からの反射光による干渉は観察されず、細胞の接着面により限定した情報のみを取り出すことが可能である。
【００９６】
本実施例では、広い放射波長を有するハロゲンランプを光源として用い、出力光に４２０ｎｍから８００ｎｍ程度の可視から近赤外の波長域において任意の広い波長幅のバンドパスフィルタを用いている。中心波長が５００ｎｍから１０００ｎｍで、半値全幅が１００ｎｍ以上のバンドパスされた光を用いると、細胞の接着面に限定した反射干渉画像を取り出すことができる。このような照明法の工夫は細胞からの定量的なパラメータを取得する上で非常に有効である。
【００９７】
［第４実施例］
引き続き、本発明の第４実施例について説明する。第４実施例では、図１に示した第１実施例の構成要素を全て含み、反射干渉画像取得のための照明方法に更なる特徴がある。図１０は、第４実施例にかかる細胞観察装置１Ａの構成概要図である。図１０に示すように、細胞観察装置１Ａには、リングスリット１１３が更に備えられている。
【００９８】
一般に、統計的に有意なデータを取得したり、多くの細胞の中にまれに混在する異種細胞を発見するためには、できるだけ多くの細胞について計測を行う必要がある。そのため、できるだけ低倍率の対物レンズを用いて一度に広い視野を計測することが望まれる。しかしながら、低倍率の対物レンズはＮＡが低く、ＮＡの低い照明光は標本に垂直に照射される成分が多く、細胞上部に存在する溶液と空気の界面で反射し、反射光によって細胞が照明されるという現象を引き起こす。これによって接着面とは無関係な細胞の形態像が観察側に含まれるようになる。また、対物レンズの中心部分を通過した照明光は、対物レンズ内部で反射され、反射光が観察側に多く含まれて高い背景光となり、接着面の反射干渉画像のコントラストを低下させる原因となる。
【００９９】
本実施例にかかる細胞観察装置１Ａの反射干渉計測において、図１０に示すように、低いＮＡの対物レンズ１０４Ａを使う場合は、対物レンズ１０４Ａの反射干渉計測用光源１０６側の開口絞りと共役な位置にリング状のスリット１１３を設ける。リング状に開口したスリット１１３を反射干渉計測用光源１０６からの照明光が通り、照明光は対物レンズ１０４Ａの中心は通らず周辺を通って照明され、角度のついた高いＮＡの光のみを使って細胞１０１が照明されるため、細胞１０１上部の溶液からの反射光の影響を低減することができる。また、リング状のスリット１１３を用いた照明は、低ＮＡの対物レンズ１０４Ａのみならず、広く対物レンズ１０４Ａ内部の反射による背景光を低減することができる。リング状に開口したスリット１１３は使用する対物レンズ１０４Ａの瞳径に適合して対物レンズ１０４Ａごとに変えられるように、使用する対物レンズ１０４Ａに適合したリングスリット１１３が円盤上に複数個取り付けられて、必要に応じて、円盤を回転させることによって、使用者がリングスリット１１３を選択できるように構成されている。
【０１００】
図１１にリングスリット１１３を用いてない場合と、用いた場合の反射干渉画像を示す。図１１（Ａ）はリングスリット１１３を用いてない場合の反射干渉画像を示す。細胞の上部の溶液と空気との界面からの反射光によって細胞が照明され、矢印で示すように、細胞の輪郭像が細胞の接着画像に重なって観察されてしまう。一方、図１１（Ｂ）のように照明光学系にリングスリット１１３を設けた場合には、細胞の上部の溶液と空気との界面に直進する照明成分がカットされるために、細胞が上部からの反射光によって照明されることがなくなり、細胞の接着面の情報をコントラスト良く観察することができる。このような照明法の工夫は細胞からの定量的なパラメータを取得する上で非常に有効である。
【０１０１】
［第５実施例］
引き続き、本発明の第５実施例について説明する。第５実施例では、図１に示した第１実施例の構成要素を全て含み、細胞１０１の入った容器に更なる特徴がある。
【０１０２】
図１２は、第５実施例における容器１０２Ａを示す図である。容器１０２Ａの有する特徴の一つは、細胞１０１が入れられた容器１０２Ａの細胞１０１の接着面の反対側（容器１０２Ａの観察領域の対物レンズ側）に反射防止コート１１４が施されたことである（以下、「特徴１」ともいう。）。容器１０２Ａの有する他の特徴の一つは、容器１０２Ａの上部に、培養液１１５の高さを一定にできるように、培養液１１５に接触するような高さに平行にガラスなどの観察窓１１６を備えたことである（以下、「特徴２」ともいう。）。容器１０２Ａの有する更に他の特徴の一つは、観察窓１１６において、容器１０２Ａの外側（培養液１１５に接する面の反対側の空気側）に反射防止コート１１７が施されたことである（以下、「特徴３」ともいう。）。
【０１０３】
特徴１は、反射干渉画像を取得する上で著しい効果を発揮する。できるだけ多くの細胞１０１を計測する目的では、低倍率対物レンズを利用することが良い。しかし、低倍率の対物レンズはＮＡが低く、通常油浸や水浸ではなく、ドライ系の対物レンズが一般的である。しかしながら、反射干渉画像において、ドライ系の対物レンズを用いると、対物レンズから出射した照明光は細胞１０１の入っている容器の底面で大きく反射する。これは空気と容器の底面のガラスなどの屈折率差が大きいためである。そのため、大幅に背景光が増加して細胞１０１の接着面の反射干渉画像を見ることはほとんど困難となる。従来、反射干渉観察に、油浸や水浸の対物レンズが用いられてきた理由はここにある。
【０１０４】
そこで、第５実施例では、細胞１０１が入れられた容器１０２Ａの細胞１０１の接着面の反対側（容器１０２Ａの観察領域の対物レンズ側）に反射防止コート１１４を施すことによって、ドライ系の対物レンズでの反射干渉画像の取得を可能にした。図１３に示すように、反射防止コート１１４をしていない状態で対物レンズから照射された照明光は空気とガラスとの界面で４％程度の反射があるが、ガラス底面に照明光の波長範囲（４２０ｎｍ〜７２０ｎｍ）において反射率Ｒを０．５％程度に抑えた反射防止コート処理を行うことによって、背景光を１／８以下に低減することができる。このため、第５実施例によれば、ドライ系の対物レンズを用いる場合であっても、高いコントラストで細胞接着面の反射干渉画像を得ることができる。
【０１０５】
図１４に容器１０２Ａの底面に反射防止コート１１４を施した場合のドライ対物レンズを用いた反射干渉画像を示す。図１４（Ａ）のように反射防止コート１１４を施さないと、容器１０２Ａの底面からの反射が大きいために、反射干渉画像のコントラストは極めて低下してしまう。これに対して、図１４（Ｂ）のように底面に反射防止コート１１４を施した容器１０２Ａを用いることによって、ドライ系の対物レンズでも良好なコントラストを得ることができる。
【０１０６】
特徴２は、定量位相画像を取得する上で極めて有効である。定量位相画像を取得する上で、溶液の揺れによって生じる液面の高さの変化は光路長の変化となり、計測上のノイズを生じる。そのため、液面の高さが変わらないようにして計測を行うことが望まれる。しかしながら、継続的に観察を行う場合、途中で薬液を投与してその後の変化を計測する場合があり、薬液投与によって液面の揺れが生じることは避けられない。そこで、第５実施例では、細胞１０１の入った容器１０２Ａの上部（または蓋）に特徴がある。
【０１０７】
すなわち、容器１０２Ａの上部（または蓋）には、培養液１１５の高さを一定にできるように、培養液１１５に接触するような高さに平行にガラスなどの観察窓１１６が設けられている。観察窓１１６の存在により、細胞１０１の入った容器１０２Ａの観察領域における溶液の高さは、細胞１０１の接着した観察面である容器１０２Ａの底面と容器１０２Ａの上部の観察窓１１６との間で一定となる。観察窓１１６は容器１０２Ａの底面の観察領域に対応した範囲に設けられ、周囲は開口部を設けて薬液の投与を可能にしている。なお、定量位相画像を取得するにあたり、観察窓１１６は標本の底面に対して並行であることが望ましい。容器１０２Ａの上部に培養液１１５と接して培養液１１５の高さを一定に保つ観察窓１１６を設けた場合には、定量位相画像は時間的に安定であり、細胞１０１の光路長の測定結果を常に安定な条件下で得ることが可能となる。なお、容器１０２Ａ上部の観察窓１１６は容器１０２Ａの蓋を兼ねており、蓋を閉じた状態で無菌状態を維持することが可能である。また、薬液を投与する場合を考慮して蓋部分には薬液投与用の開口部を設けてもよい。
【０１０８】
特徴３は、反射干渉画像を取得するにあたり、観察窓１１６の容器１０２Ａの外側（培養液１１５に接する面の反対側の空気側）に、反射防止コート１１７が施されたことである。反射防止コート１１７を施さなかった場合には、細胞１０１上部の溶液に接する観察窓１１６のガラスと外側の空気との屈折率の差が大きいために、反射干渉画像を得る目的で対物レンズから照射された照明光が観察窓１１６において４％程度反射される。その反射光により細胞１０１が照明されて反射干渉画像に細胞１０１の輪郭像が重ねて投影される。このことは反射干渉画像から接着面に限定した情報を取り出す上で邪魔な情報となり、クリアな反射干渉画像を得ることの妨げとなる。そこで、第５実施例では、観察窓１１６の容器１０２Ａの外側に照明光の波長範囲において反射率を０.５％程度に抑えるような反射防止コート処理を施すことによって、反射率を１/８程度に低減することができる。このため、上面からの反射光の影響を抑えて、高いコントラストで細胞接着面の反射干渉画像を得ることができる。
【０１０９】
図１５に容器１０２Ａ上部の観察窓１１６に反射防止コート１１７を施した場合の反射干渉画像の改善を示す。図１５（Ａ）のように無処理の場合では、容器１０２Ａ上部の観察窓１１６と空気との界面からの反射光によって細胞１０１が上部から照明されて細胞１０１の輪郭像が映り込んでしまう。一方で、図１５（Ｂ）のように上部の観察窓１１６の空気側に反射防止コート１１７を施すことによって、上部からの反射は低減され、細胞１０１の接着面の反射干渉画像をコントラスト良く観察することができる。
【０１１０】
なお、細胞１０１の入った容器１０２Ａに培養液１１５の乱れが起きないように、例えば測定中に試薬などを投与しない状態で静置観察を行うのであれば、上述の観察窓１１６を設けることは必須ではなく、図１６に示すように、容器の底面に反射防止コート１１４を施した特徴１のみを有する容器１０２Ｂとしても良い。
【０１１１】
［第６実施例］
引き続き、本発明の第６実施例について説明する。第６実施例では、図１に示した第１実施例の構成要素のうち、画像取得部１０に相違点がある。つまり、定量位相計測および反射干渉計測を同時に行うための光学系の構成に特徴がある。
【０１１２】
図１７は、第６実施例における画像取得部１０Ａを示す。第１実施例の画像取得部１０と比べると、全反射ミラー１０９を用いずに構成されている。つまり、定量位相光学系は回折型干渉光学系１１１のみとなる。このような構成により、簡素な構成が実現でき、装置化しやすいというメリットがある。
【０１１３】
定量位相計測において、細胞１０１の入った容器１０２の上部に設置した定量位相計測用光源１０５から放射された照明光は、細胞１０１の入った容器１０２を透過し、対物レンズ１０４にて集光される。そして、ハーフミラー１０７を介し、定量位相像と反射干渉像を波長で分離するためのダイクロイックミラー１０８を透過して、位相計測のための回折型干渉光学系１１１にて物体光と参照光の干渉像を形成し、定量位相検出用カメラ１１２にて干渉縞画像が撮像される。
【０１１４】
反射干渉計測において、反射干渉計測用光源１０６から放射された照明光は、ハーフミラー１０７にて反射して、対物レンズ１０４を通って、測定対象となる細胞１０１を入れた容器１０２の底面側から入射する。容器１０２の底面の細胞１０１の接着面から得られる反射光は細胞１０１の接着距離に応じて干渉を生じ、得られた反射干渉光は再び対物レンズ１０４にて集光され、ハーフミラー１０７を介し、定量位相像と反射干渉像を波長で分離するためのダイクロイックミラー１０８にて反射干渉像のみを反射して、反射干渉検出用カメラ１１０にて撮像される。
【０１１５】
［第７実施例］
引き続き、本発明の第７実施例について説明する。第７実施例では、図１に示した第１実施例の構成要素のうち、画像取得部１０に相違点がある。つまり、定量位相計測および反射干渉計測を同時に行うための光学系の構成に特徴がある。
【０１１６】
図１８は、第７実施例における画像取得部１０Ｂを示す。第１実施例の画像取得部１０と比べると、全反射ミラー１０９を用いずに構成されている。従って、簡素な構成が実現でき、装置化しやすいというメリットがある。また、反射干渉入射用光学系であるハーフミラー１０７が共通光学系となっていないため、定量位相像はハーフミラー１０７を通過せずに撮像され、定量位相側の光量ロスを低減できるというメリットがある。
【０１１７】
定量位相計測において、細胞１０１の入った容器１０２の上部に設置した定量位相計測用光源１０５から放射された照明光は、細胞１０１の入った容器１０２を透過し、対物レンズ１０４にて集光される。そして、ハーフミラー１０７を介することなく、定量位相像と反射干渉像を波長で分離するためのダイクロイックミラー１０８を透過して、位相計測のための回折型干渉光学系１１１にて物体光と参照光の干渉像を形成し、定量位相検出用カメラ１１２にて干渉縞画像が撮像される。
【０１１８】
反射干渉計測において、反射干渉計測用光源１０６から放射された照明光は、ハーフミラー１０７およびダイクロイックミラー１０８にて反射して、対物レンズ１０４を通って、測定対象となる細胞１０１を入れた容器１０２の底面側から入射する。容器１０２の底面の細胞１０１の接着面から得られる反射光は細胞１０１の接着距離に応じて干渉を生じ、得られた反射干渉光は対物レンズ１０４にて集光され、ダイクロイックミラー１０８およびハーフミラー１０７を再び介し、反射干渉検出用カメラ１１０にて撮像される。
【０１１９】
［第８実施例］
引き続き、本発明の第８実施例について説明する。第８実施例では、図１に示した第１実施例の構成要素のうち、画像取得部１０に相違点がある。つまり、定量位相計測および反射干渉計測を同時に行うための光学系の構成に特徴がある。
【０１２０】
図１９は、第８実施例における画像取得部１０Ｃを示す。第１実施例の画像取得部１０と比べると、全反射ミラー１０９を用いずに構成されている。つまり、定量位相光学系は回折型干渉光学系１１１のみとなり、簡素な構成が実現でき、装置化しやすいというメリットがある。また、反射干渉入射用光学系であるハーフミラー１０７が共通光学系となっていないため、定量位相像はハーフミラー１０７を通過せずに撮像されるので、定量位相側の光量ロスを低減できるというメリットがある。
【０１２１】
定量位相計測において、細胞１０１の入った容器１０２の上部に設置した定量位相計測用光源１０５から放射された照明光は、細胞１０１の入った容器１０２を透過し、対物レンズ１０４にて集光される。そして、ハーフミラー１０７を介することなく、定量位相像と反射干渉像を波長で分離するためのダイクロイックミラー１０８から反射して、位相計測のための回折型干渉光学系１１１にて物体光と参照光の干渉像を形成し、定量位相検出用カメラ１１２にて干渉縞画像が撮像される。
【０１２２】
反射干渉計測において、反射干渉計測用光源１０６から放射された照明光は、ハーフミラー１０７およびダイクロイックミラー１０８にて反射して、対物レンズ１０４を通って、測定対象となる細胞１０１を入れた容器１０２の底面側から入射する。容器１０２の底面の細胞１０１の接着面から得られる反射光は細胞１０１の接着距離に応じて干渉を生じ、得られた反射干渉光は対物レンズ１０４にて集光され、ダイクロイックミラー１０８およびハーフミラー１０７を再び介し、反射干渉検出用カメラ１１０にて撮像される。
［第９実施例］
引き続き、本発明の第９実施例について説明する。第９実施例では、図１に示した第１実施例の構成要素のうち、画像取得部１０に相違点がある。つまり、定量位相計測および反射干渉計測を同時に行うための光学系の構成に特徴がある。
【０１２３】
図２０は、第９実施例における画像取得部１０Ｄを示す。第１実施例の画像取得部１０と比べると、全反射ミラー１０９およびダイクロイックミラー１０８を用いずに構成されており、定量位相計測用光源１０５および反射干渉計測用光源１０６を一つの光源１０５Ｄに共通化している。従って、簡素な構成が実現でき、装置化しやすいというメリットがある。ただし、二つ目の対物レンズ１０４Ｄが必要である。
【０１２４】
定量位相計測において、共通の光源１０５Ｄから放射された照明光は、ハーフミラー１０７にて反射して、対物レンズ１０４を通って、細胞１０１の入った容器１０２を透過し、第２対物レンズ１０４Ｄにて集光される。そして、回折型干渉光学系１１１にて物体光と参照光の干渉像を形成し、定量位相検出用カメラ１１２にて干渉縞画像が撮像される。
【０１２５】
反射干渉計測において、共通の光源１０５Ｄから放射された照明光は、ハーフミラー１０７にて反射して、対物レンズ１０４を通って、測定対象となる細胞１０１を入れた容器１０２の底面側から入射する。容器１０２の底面の細胞１０１の接着面から得られる反射光は細胞１０１の接着距離に応じて干渉を生じ、得られた反射干渉光は対物レンズ１０４にて集光され、ハーフミラー１０７を再び介し、反射干渉検出用カメラ１１０にて撮像される。
【０１２６】
［第１０実施例］
引き続き、本発明の第１０実施例について説明する。第１０実施例では、図１に示した第１実施例の構成要素のうち、画像取得部１０に相違点がある。つまり、定量位相計測および反射干渉計測を同時に行うための光学系の構成に特徴がある。
【０１２７】
図２１は、第１０実施例における画像取得部１０Ｅを示す。第１実施例の画像取得部１０と比べると、全反射ミラー１０９、ハーフミラー１０７、およびダイクロイックミラー１０８を用いずに構成され、定量位相計測用の光と反射干渉計測用の光とを波長ではなく角度で分けている。また、二つ目の対物レンズ１０４Ｅが必要である。
【０１２８】
定量位相計測において、細胞１０１の入った容器１０２の上部に設置した定量位相計測用光源１０５から放射された照明光は、細胞１０１の入った容器１０２を透過し、対物レンズ１０４にて集光される。そして、ハーフミラー１０７やダイクロイックミラー１０８を介することなく、回折型干渉光学系１１１にて物体光と参照光の干渉像を形成し、定量位相検出用カメラ１１２にて干渉縞画像が撮像される。
【０１２９】
反射干渉計測において、反射干渉計測用光源１０６から放射された照明光は、測定対象となる細胞１０１を入れた容器１０２の底面側から斜めに入射する。容器１０２の底面の細胞１０１の接着面から得られる反射光は細胞１０１の接着距離に応じて干渉を生じ、得られた反射干渉光は再び斜めで第２対物レンズ１０４Ｅにて集光され、ハーフミラー１０７やダイクロイックミラー１０８を介することなく、反射干渉検出用カメラ１１０にて撮像される。反射干渉計測用光源１０６と、第２対物レンズ１０４Ｅおよび反射干渉検出用カメラ１１０との間の角度調整が必要であるが、全反射ミラー１０９、ハーフミラー１０７、およびダイクロイックミラー１０８を用いずに構成できるため、簡素な構成が実現できるというメリットがある。
【０１３０】
［第１１実施例］
引き続き、本発明の第１１実施例について説明する。第１１実施例では、図１に示した第１実施例の構成要素のうち、画像取得部１０に相違点がある。つまり、定量位相計測および反射干渉計測を同時に行うための光学系の構成に特徴がある。
【０１３１】
図２２は、第１１実施例における画像取得部１０Ｆを示す。第１実施例の画像取得部１０と比べると、全反射ミラー１０９、ハーフミラー１０７、およびダイクロイックミラー１０８を用いずに構成され、定量位相計測用の光と反射干渉計測用の光とを波長ではなく角度で分けている。また、二つ目の対物レンズ１０４Ｆが必要である。
【０１３２】
定量位相計測において、定量位相計測用光源１０５から放射された照明光は、測定対象となる細胞１０１を入れた容器１０２の底面側から入射する。照明光は細胞１０１の入った容器１０２を透過し、対物レンズ１０４にて集光される。そして、全反射ミラー１０９、ハーフミラー１０７、またはダイクロイックミラー１０８を介することなく、回折型干渉光学系１１１にて物体光と参照光の干渉像を形成し、定量位相検出用カメラ１１２にて干渉縞画像が撮像される。
【０１３３】
反射干渉計測において、反射干渉計測用光源１０６から放射された照明光は、測定対象となる細胞１０１を入れた容器１０２の底面側から斜めに入射する。容器１０２の底面の細胞１０１の接着面から得られる反射光は細胞１０１の接着距離に応じて干渉を生じ、得られた反射干渉光は再び斜めで第２対物レンズ１０４Ｆにて集光され、全反射ミラー１０９、ハーフミラー１０７、またはダイクロイックミラー１０８を介することなく、反射干渉検出用カメラ１１０にて撮像される。反射干渉計測用光源１０６と、第２対物レンズ１０４Ｆおよび反射干渉検出用カメラ１１０との間の角度調整が必要であるが、全反射ミラー１０９、ハーフミラー１０７、およびダイクロイックミラー１０８を用いずに構成できるため、簡素な構成が実現できるというメリットがある。
【０１３４】
［第１２実施例］
引き続き、本発明の第１２実施例について説明する。第１２実施例では、図１に示した第１実施例の構成要素のうち、画像取得部１０に相違点がある。つまり、定量位相計測および反射干渉計測を同時に行うための光学系の構成に特徴があり、具体的には、定量位相の計測方法を２光束方式に変えたものである。
【０１３５】
図２３は、第１２実施例における画像取得部１０Ｇを示す。第６実施例の画像取得部１０Ａと比べると、回折型干渉光学系１１１を用いずに、定量位相計測用光源１０５からの照明光および参照光を光路１１８Ｇおよび光路１１９Ｇを用いて定量位相検出用カメラ１１２に到達させている点が相違する。回折型干渉光学系１１１を備えなくても装置構成ができるというメリットがある。定量位相の計測方法を回折型干渉光学系１１１を用いる代わりに、２光束方式を採用し、参照光を別の光路で作成し、最終的に物体光と参照光を干渉させる構成である。
【０１３６】
定量位相計測において、定量位相計測用光源１０５からの照明光は、照明光のための光路１１８Ｇを通って、容器１０２の上部に設置したレンズ１２０Ｇを介して、容器１０２内の細胞１０１を透過し、対物レンズ１０４にて集光される。そして、ハーフミラー１０７を介し、ダイクロイックミラー１０８を透過して、ハーフミラー１２１Ｇおよびレンズ１２２Ｇを更に介し、定量位相検出用カメラ１１２に到達する。一方、定量位相計測用光源１０５からの参照光は、参照光のための光路１１９Ｇを通って、レンズ１２３Ｇ、ハーフミラー１２１Ｇ、およびレンズ１２２Ｇを順に介し、定量位相検出用カメラ１１２に到達する。
【０１３７】
反射干渉計測において、反射干渉計測用光源１０６から放射された照明光は、ハーフミラー１０７にて反射して、対物レンズ１０４を通って、測定対象となる細胞１０１を入れた容器１０２の底面側から入射する。容器１０２の底面の細胞１０１の接着面から得られる反射光は細胞１０１の接着距離に応じて干渉を生じ、得られた反射干渉光は再び対物レンズ１０４にて集光され、ハーフミラー１０７を介し、定量位相像と反射干渉像を波長で分離するためのダイクロイックミラー１０８にて反射干渉画像のみを反射して、反射干渉検出用カメラ１１０にて撮像される。
【０１３８】
［第１３実施例］
引き続き、本発明の第１３実施例について説明する。第１３実施例では、図１に示した第１実施例の構成要素のうち、画像取得部１０に相違点がある。つまり、定量位相計測および反射干渉計測を同時に行うための光学系の構成に特徴がある。
【０１３９】
図２４は、第１３実施例における画像取得部１０Ｈを示す。第６実施例の画像取得部１０Ａと比べると、照明波長切替装置１２４Ｈおよび観察波長切替装置１２５Ｈが更に備えられている点が相違する。照明波長切替装置１２４Ｈは、反射干渉計測用光源１０６からの照明光の波長を切り替える装置である。観察波長切替装置１２５Ｈは、ダイクロイックミラー１０８からの観察光の波長を切り替える装置である。これらの装置を備えることにより、反射干渉検出用カメラ１１０にて所望の蛍光画像等を取得することができる。なお、図２４においては、回転により波長の切り替えを行う照明波長切替装置１２４Ｈおよび観察波長切替装置１２５Ｈを一例として示している。
【０１４０】
定量位相計測において、細胞１０１の入った容器１０２の上部に設置した定量位相計測用光源１０５から放射された照明光は、細胞１０１の入った容器１０２を透過し、対物レンズ１０４にて集光される。そして、ハーフミラー１０７を介し、定量位相像と反射干渉像を波長で分離するためのダイクロイックミラー１０８を透過して、位相計測のための回折型干渉光学系１１１にて物体光と参照光の干渉像を形成し、定量位相検出用カメラ１１２にて干渉縞画像が撮像される。
【０１４１】
反射干渉計測において、反射干渉計測用光源１０６から放射された照明光は、照明波長切替装置１２４Ｈにて所望の波長に切り替えられ、ハーフミラー１０７にて反射して、対物レンズ１０４を通って、測定対象となる細胞１０１を入れた容器１０２の底面側から入射する。容器１０２の底面の細胞１０１の接着面から得られる反射光は細胞１０１の接着距離に応じて干渉を生じ、得られた反射干渉光は対物レンズ１０４にて集光され、ハーフミラー１０７を再び介し、ダイクロイックミラー１０８にて反射干渉画像のみを反射して、当該反射光は観察波長切替装置１２５Ｈにて所望の波長に切り替えられ、反射干渉検出用カメラ１１０にて撮像される。
【０１４２】
［第１４実施例］
引き続き、本発明の第１４実施例について説明する。第１４実施例では、図１に示した第１実施例の構成要素のうち、ハーフミラー１０７に相違点がある。つまり、第１４実施例は、第１実施例の構成要素を全て含み、ハーフミラー１０７の特性をより向上させたものである。
【０１４３】
通常のハーフミラーは全波長域において同等の分離比率を有するものが一般的である。このようなハーフミラーを本実施形態におけるハーフミラー１０７として用いた場合では、反射干渉の目的で光路中に設置したハーフミラー１０７を定量位相の像が通過する際に、定量位相画像の透過率はハーフミラー１０７によって制限されてしまう。そこで、第１４実施例では、図２５の画像取得部１０Ｉに示すように、反射干渉計測用光源１０６とカメラ１１０との間および定量位相計測用光源１０５とカメラ１１０との間に位置し、波長により反射対透過の比率を異ならせることが可能なハーフミラー１０７Ａを用いる。ここで、第１４実施例においては、反射干渉と定量位相の照明光の波長が重ならないように、定量位相計測用光源１０５および反射干渉計測用光源１０６はそれぞれ別の波長の照明光を用いることが前提である。ハーフミラー１０７Ａの波長特性とは、つまり、反射干渉の観察波長域においてはハーフミラー（図２５では透過：反射＝８０：２０）として働き、一方、定量位相の観察波長域においては最大限光を透過する波長特性である。こうすることにより、定量位相の観察光の光量のロスが少なくて済む。
【０１４４】
例えば、反射干渉では可視光（４２０ｎｍから７５０ｎｍ）を用い、定量位相では中心波長８３０ｎｍで半値全幅２０ｎｍの光を用いる場合には、ハーフミラー１０７Ａの波長特性を可視光（４２０ｎｍから７５０ｎｍ）では反射：透過の比率を２０：８０に設定する一方で、８００ｎｍ以上では反射：透過の比率を５：９５に設定する。こうすることで、反射干渉では観察光の８０％を得ることができるとともに、定量位相については、ハーフミラー１０７Ａの反射干渉時の透過率に左右されず、９５％の透過率で信号を得ることが可能となる。
【０１４５】
［第１５実施例］
引き続き、本発明の第１５実施例について説明する。第１５実施例では、図１に示した第１実施例の構成と同様の構成を有するが、解析部２０７がリファレンス記憶部２０８によらずに処理を行うため、図示まではしないが、図１におけるリファレンス記憶部２０８を除いた構成としても良い。
【０１４６】
第１５実施例において、解析部２０７は、リファレンスデータに基づかないで、未知の被検体の細胞の種類または状態を判別する。つまり、解析部２０７は、未知の細胞に対して抽出したパラメータの中から所定のパラメータを選択し、当該選択された所定のパラメータを用いて、未知の細胞の種類または状態を判別する。ここで、解析部２０７は、特異な値を示すパラメータを所定のパラメータとして自動的に選択するように構成されていても良い。また、当該パラメータ選択処理を行うために、別途の処理手段が設けられていても良い。一方で、解析部２０７が、ユーザからの入力に基づき、所定のパラメータを選択するように構成されていても良い。
【０１４７】
図２６は、第１５実施例において、解析部２０７がユーザからの入力に基づき所定のパラメータを選択するように構成された場合に、ユーザに表示される画面４０を示す。画面４０Ａは、ユーザにパラメータの種類を示し、ユーザの選択を受けるための画面である。画面４０Ａの例において、ユーザは、慣性と尖度を所定のパラメータとして選択している。なお、選択ツール４０Ａ１等を用いて、ユーザが画面４０Ａに表示されていない他のパラメータを選択できるようにしても良い。画面４０Ｂは、ユーザが選択したパラメータを反映してグラフに表示された例を示す。ユーザは画面４０Ｂを確認してから、画面４０Ａに戻りパラメータを選択し直しても良く、そのまま画面４０Ｃに進んでも良い。画面４０Ｃは、細胞の画像を表示する画面である。
【０１４８】
このような第１５実施例によれば、解析部２０７は、リファレンスデータによらずに、処理を行うことができる。したがって、リファレンス記憶部２０７を設けなくて済むので、装置構成が簡潔になる。このような実施態様は、例えば、試料中の細胞が骨髄幹細胞と破骨細胞とに分かれていることが事前に分かっているときに、細胞それぞれを分別する場合に好適である。
【０１４９】
［第１６実施例］
引き続き、本発明の第１６実施例について説明する。第１６実施例では、図１に示した第１実施例の構成要素のうち、主に、解析部２０７の動作に相違点がある。第１６実施例では、第１抽出部２０４、第２抽出部２０５、および第３抽出部２０６のそれぞれ、または組み合わせにより、３つ以上のパラメータが抽出された場合を前提としており、この場合に解析部２０７は、当該３つ以上のパラメータに対して主成分分析を行って、未知の細胞の種類または状態を判別する。
【０１５０】
第１６実施例における動作を簡単にまとめると、次のようである。最初に、第１抽出部２０４、第２抽出部２０５、および第３抽出部２０６のそれぞれ、または組み合わせが、３つ以上のパラメータを抽出して、解析部２０７に出力する（以下「手順１」）。次に、解析部２０７が、当該３つ以上のパラメータに対して、標準化を行う（以下「手順２」）。次に、解析部２０７が、当該標準化データに基づき、第１主成分および第２主成分を求める（以下「手順３」）。次に、解析部２０７が、第１主成分および第２主成分の２成分散布図に細胞をプロットし、細胞判別を行う（以下「手順４」）。以下、実際の実験例に基づき、各手順の詳細について説明する。
【０１５１】
（手順１：パラメータ抽出）
本実験においては、それぞれ１０個ずつの４種類の細胞をリファレンスとして用い、３つの未知の細胞に対して識別を実施した。なお、３つの未知の細胞は４種の細胞の何れかに含まれるものとして実験を行い、実験の目的は３つの未知の細胞が４種の細胞の何れに含まれるものかを識別することである。以下、リファレンスとして用いた各細胞の名称と、表記名と、個数（ｎ）とを順に示す。
ラット膵β細胞株  ＩＮＳ−１  （ｎ＝１０）
マウス膵β細胞株  ＭＩＮ−６  （ｎ＝１０）
ヒト膵癌細胞株  ＭＩＡ ＰａＣａ−２  （ｎ＝１０）
ヒト子宮頸癌細胞株  ＨｅＬａ  （ｎ＝１０）
【０１５２】
図２７は、リファレンスとして用いた各細胞の反射干渉画像および定量位相画像を示す。図２７（Ａ）および（Ｂ）はラット膵β細胞株（ＩＮＳ−１）の反射干渉画像および定量位相画像を順に示す。図２７（Ｃ）および（Ｄ）はマウス膵β細胞株（ＭＩＮ−６）の反射干渉画像および定量位相画像を順に示す。図２７（Ｅ）および（Ｆ）はヒト膵癌細胞株（ＭＩＡＰａＣａ−２）の反射干渉画像および定量位相画像を順に示す。図２７（Ｇ）および（Ｈ）はヒト子宮頸癌細胞株（ＨｅＬａ）の反射干渉画像および定量位相画像を順に示す。
【０１５３】
本実験では、第１抽出部２０４、第２抽出部２０５、および第３抽出部２０６のそれぞれ、または組み合わせが、図２７に示した反射干渉画像および定量位相画像から以下に示す７つのパラメータを抽出した。下記の（１）と（２）が定量位相画像より得られたパラメータであり、（３）〜（７）が反射干渉画像より得られたパラメータである。図２８は、７つのパラメータの実際の計測値を示す。
（１）細胞の面積
（２）細胞の厚さ
（３）接着面のテクスチャ：同時生起行列／局所一様性
（４）接着面のテクスチャ：同時生起行列／慣性
（５）接着面のコントラスト：濃度ヒストグラム／歪度
（６）接着面のコントラスト：濃度ヒストグラム／尖度
（７）接着面のコントラスト：濃度ヒストグラム／√（分散）／平均
【０１５４】
（手順２：標準化）
上記手順１で抽出した複数のパラメータから、要約した特徴を示す成分（主成分）を求めるには、各パラメータを含む１次元の式を作成する必要がある。しかし、図２８に示す計測値では、パラメータによって単位が異なるために、取りうる値の大きさがまちまちである。これらの値をそのまま取り扱うと、値の大きいパラメータに情報が偏り、各パラメータから均等に情報を取り出すことが難しい。そこで、各パラメータから均等に情報を取り出すために、各パラメータの情報量を求めた。情報量として、母集団に対した偏差平方和を図２８に示す計測値から求めてみたところ、図２９に示すような計算結果が出た。しかし、図２９に示されるように、各パラメータ間の情報量に大きな差があった。各パラメータの情報量を等しくすることが望ましいので、本実験では、各パラメータにつき、以下の数式（１）に基づき、データの標準化を行った。
Ｘ’＝（Ｘｉ―Ｘ）／Ｘｓｄ…（１）
なお、数式（１）において、Ｘ’は標準化後のデータであり、Ｘｉは図２８で示した各計測値であり、ＸはＸｉにおける平均値であり、ＸｓｄはＸｉにおける標準偏差である。
【０１５５】
図３０は標準化後の各パラメータの値を示す。図３０に示されるように、標準化の結果、各パラメータの情報量はパラメータの種類によらず等しくなっている。このような標準化後のデータを用いることにより、平等にデータを扱うことができるようになる。
【０１５６】
（手順３：主成分の計算）
次に、解析部２０７が、図３０の標準化後のデータに基づき、第１主成分ｆおよび第２主成分ｇを求める。まず、標準化したデータをそれぞれ、面積＝Ｘ１、光学厚さ＝Ｘ２、…、√（分散）／平均＝Ｘ７とし、それらから統合した特徴を表す２つの主成分であるｆとｇを以下の数式（２）および（３）に基づき計算する。なお、数式（２）および（３）において、ａ１からａ７、ｂ１からｂ７は係数である。
ｆ＝ａ１＊Ｘ１＋ａ２＊Ｘ２＋ … ＋ａ７＊Ｘ７…（２）
ｇ＝ｂ１＊Ｘ１＋ｂ２＊Ｘ２＋ … ＋ｂ７＊Ｘ７…（３）
【０１５７】
このとき、係数には以下の条件を設定する。
ａ１２＋ａ２２＋ … ＋ａ７２＝１…（４）
ｂ１２＋ｂ２２＋ … ＋ｂ７２＝１…（５）
数式（４）および（５）で示す条件を設定する理由は、後の手順で主成分ｆ、ｇの分散を最大にするような係数ａ１からａ７、ｂ１からｂ７を求めるので、数式（４）および（５）の条件にて係数を制限しておかないと、係数が制限なく大きくなってしまうからである。「ａ１２＋ａ２２＋・・・ ＋ａ７２」および「ｂ１２＋ｂ２２＋… ＋ｂ７２」は係数のベクトルの大きさに相当するものであり、ベクトルの大きさを１とすることで、元データの情報量を変えずに、新しい判断軸としての第１主成分ｆおよび第２主成分ｇを得ることができる。
【０１５８】
また、ｆとｇを相関のない互いに直交する直線とするために、以下の条件を更に設定する。
ａ１＊ｂ１＋ａ２＊ｂ２＋ … ＋ａ７＊ｂ７＝０…（６）
【０１５９】
ここで、ｆについて元になるデータの情報量をできるだけ多く反映した１次式となるようにするには、ｆで与えられる成分の値のばらつきが最も大きくなるようにすることが考えられる。ばらつきは分散と考えることができるので、ｆの分散が最大になるようにする。標本数４３の細胞それぞれのｆの値をｆ１、ｆ２、ｆ３、…、ｆ４３とすると、ｆの分散Ｖｆは以下の数式（７）で与えられ、分散Ｖｆが最大になる条件を満たすｆを計算する。
【数１】
2011132587.tif 000003
なお、上記の数式（７）において、
【数２】
2011132587.tif 000004
はｆの平均値である。
【０１６０】
以上の方法により、第１主成分ｆについて、上記数式（４）の条件を満たし、且つ上記数式（７）の条件を満たす係数ａ１、ａ２、…、ａ７を計算する。図３１は、以上で説明した計算方法に基づき計算された主成分ｆの係数を示す。
【０１６１】
次に、第２主成分ｇについて、上記数式（５）の条件を満たし、且つ上記数式（７）と同趣旨の条件（つまり分散が最大になる条件）を満たし、更に、上記数式（６）の条件を満たすような係数ｂ１、ｂ２、…、ｂ７を計算する。図３２は、以上で説明した計算方法に基づき計算された第２主成分ｇの係数を示す。
【０１６２】
次に、各細胞についてｆおよびｇの値を求め、それぞれ第１主成分および第２主成分とする。図３３は、各細胞について第１主成分ｆおよび第２主成分ｇを求めた結果を示す。図３３（Ａ）はラット膵β細胞株（ＩＮＳ−１）およびマウス膵β細胞株（ＭＩＮ−６）の第１主成分ｆおよび第２主成分ｇを示しており、図３３（Ｂ）はヒト膵癌細胞株（ＭＩＡＰａＣａ−２）およびヒト子宮頸癌細胞株（ＨｅＬａ）の第１主成分ｆおよび第２主成分ｇを示している。
【０１６３】
（手順４：プロットおよび細胞判別）
解析部２０７が、図３３に示した第１主成分ｆおよび第２主成分ｇの値を用いて、２成分散布図に各細胞をプロットし、細胞判別を行う。図３４はプロットした結果を示しており、第１主成分ｆを横軸に、第２主成分ｇを縦軸にして、各細胞をプロットした例である。図３４のプロット結果で示されるように、第１主成分ｆと第２主成分ｇの散布図（主成分散布図）により、４種類の細胞集団はその特徴量に応じて容易に識別することのできる集団を形成した。これは主成分の分散が最も大きくなるように求めることにより（上記数式（７））、各細胞の特徴が明示できたことを示すものである。その結果、主成分散布図において未知の細胞（○で表示したｕｎｋｎｏｗｎ）が△で表示したＭＩＡ−ＰａＣａ細胞の集団に属すると判別することができた。実際に未知の細胞はＭＩＡ−ＰａＣａであったため、本実験による細胞判別は正しかったことが確認できた。
【０１６４】
以上で説明した実施例１６は、抽出した多数のパラメータの全部また３つ以上であって抽出パラメータの一部を細胞判別に有効活用したい場合に特に好適である。例えば、判別したい細胞の種類が多くて多数のパラメータを用いて高い精度で判別すべき場合や、抽出した多数のパラメータから２つのパラメータを選択することが容易ではない場合等に、特に好適である。実施例１６では、抽出した多数のパラメータに対して主成分分析を行うことにより、当該多数のパラメータが細胞判別に適切に影響を及ぼうようにすることができる。主成分分析によれば、抽出した多数のパラメータのうち２つを選択する代わりに、当該多数のパラメータの全部または３つ以上の一部が適切に反映された２つの主成分により細胞判別を行うからである。
【０１６５】
なお、以上の説明では、パラメータの数が７つである場合を例に挙げて説明したが、これに限らず、３つ以上のパラメータであれば本実施例を適用することができる。重複するため説明を省略するが、実際１０個のパラメータを用いて上記実験と同様の実験を行った結果、７つのパラメータを用いた場合と同様の結果が得られた。
【符号の説明】
【０１６６】
１，１Ａ…細胞観察装置、１０，１０Ａ，１０Ｂ，１０Ｃ，１０Ｄ，１０Ｅ，１０Ｆ，１０Ｇ，１０Ｈ…画像取得部、２０…処理部、３０…標本、３１…格子、４０，４０Ａ，４０Ｂ，４０Ｃ…画面、１０１…細胞、１０２，１０２Ａ，１０２Ｂ…容器、１０３…培養空間、１０４，１０４Ａ，１０４Ｄ，１０４Ｅ，１０４Ｆ…対物レンズ、１０５…定量位相計測用光源、１０５Ａ…照明光絞り部、１０５Ｄ…共通の光源、１０６…反射干渉計測用光源、１０７，１０７Ａ…ハーフミラー、１０８…ダイクロイックミラー、１０９…全反射ミラー、１１０…反射干渉検出用カメラ、１１１…回折型干渉光学系、１１２…定量位相検出用カメラ、１１３…リングスリット、１１４…反射防止コート、１１５…培養液、１１６…観察窓、１１７…反射防止コート、１１８Ｇ…定量位相計測における照明光のための光路、１１９Ｇ…定量位相計測における参照光のための光路、１２０Ｇ，１２２Ｇ，１２３Ｇ…レンズ、１２１Ｇ…ハーフミラー、１２４Ｈ…照明波長切替装置、１２５Ｈ…観察波長切替装置、２０１…画像位置合わせ部、２０２…輪郭抽出部、２０３…輪郭適用部、２０４…第１抽出部、２０５…第２抽出部、２０６…第３抽出部、２０７…解析部、２０８…リファレンス記憶部。
【産業上の利用可能性】
【０１６７】
本発明は、細胞の状態を適切に判別・評価するためのより多くの情報を得ることができる細胞観察装置および細胞観察方法を提供する。
【図１】
2011132587.tif 000005
【図２】
2011132587.tif 000006
【図３】
2011132587.tif 000007
【図４】
2011132587.tif 000008
【図６】
2011132587.tif 000009
【図１０】
2011132587.tif 000010
【図１２】
2011132587.tif 000011
【図１３】
2011132587.tif 000012
【図１６】
2011132587.tif 000013
【図１７】
2011132587.tif 000014
【図１８】
2011132587.tif 000015
【図１９】
2011132587.tif 000016
【図２０】
2011132587.tif 000017
【図２１】
2011132587.tif 000018
【図２２】
2011132587.tif 000019
【図２３】
2011132587.tif 000020
【図２４】
2011132587.tif 000021
【図２５】
2011132587.tif 000022
【図２８】
2011132587.tif 000023
【図２９】
2011132587.tif 000024
【図３０】
2011132587.tif 000025
【図３１】
2011132587.tif 000026
【図３２】
2011132587.tif 000027
【図３３】
2011132587.tif 000028
【図３４】
2011132587.tif 000029
【図５】
2011132587.tif 000030
【図７】
2011132587.tif 000031
【図８】
2011132587.tif 000032
【図９】
2011132587.tif 000033
【図１１】
2011132587.tif 000034
【図１４】
2011132587.tif 000035
【図１５】
2011132587.tif 000036
【図２６】
2011132587.tif 000037
【図２７】
2011132587.tif 000038
【国際調査報告】
2011132587.tif 000039
2011132587.tif 000040
2011132587.tif 000041
2011132587.tif 000042