(57)【要約】      （修正有）【課題】多孔膜上で細胞を非侵襲的に観察する定量位相画像生成装置を提供する。【解決手段】定量位相画像生成装置１は、顕微鏡１００と、情報処理部４０とを備える。顕微鏡１００は、対象物Ｓに照明光を照射する透過照明光学系１０と、ステージ８と、結像光学系７と、落射照明光学系１１０と、検出部９とを備え、結像光学系７は、対物光学系２０と、フィルターキューブ１２０と、リレー光学系３０とを備え、透過照明光学系１０は、光源１１と、レンズ１２と、バンドパスフィルター１３と、視野絞り１４と、レンズ１５と、開口絞り１６と、集光レンズ１７とを備え、対物光学系２０は、複数の対物レンズ２１を備え、リレー光学系３０は、結像レンズ３１と、ビームスプリッター３２と、ミラー３３ａ，３３ｂ，３３ｃと、レンズ３４ａ，３４ｂ，３４ｃと、接眼レンズ３５とを備える。【選択図】図１
2020188717.tif 000002
【特許請求の範囲】
【請求項１】
多孔膜上に配置された細胞の数、形態又は形状を測定する方法であって、
非干渉方式の定量位相イメージングにより前記細胞の定量位相画像を生成する工程と、
前記定量位相画像に基づいて前記細胞の数、形態又は形状を測定する工程と、を含む、方法。
【請求項２】
前記多孔膜は細胞培養等に通常用いられるポアサイズ、ポア密度、膜厚で構成される、請求項１に記載の方法。
【請求項３】
前記細胞の数、形態又は形状を測定する工程は、前記細胞の、面積、体積、真円度若しくは縦横比、又は、前記細胞の間隙の面積を測定する工程を含む、請求項１又は２に記載の方法。
【請求項４】
前記細胞の数、形態又は形状を測定する工程が、前記定量位相画像の輝度情報に基づいて行われる、請求項１〜３のいずれか一項に記載の方法。
【請求項５】
前記多孔膜が、細胞培養容器内に配置した細胞培養用インサートの底面に配置されている、請求項１〜４のいずれか一項に記載の方法。
【請求項６】
前記多孔膜が、Ｏｒｇａｎ−ｏｎ−ａ−Ｃｈｉｐ内に配置されている、請求項１〜５のいずれか一項に記載の方法。
【請求項７】
前記定量位相画像を生成する工程において、異なる時間における複数の定量位相画像を生成し、
前記細胞の数、形態又は形状を測定する工程において、前記細胞の数、形態又は形状の変化を測定する、請求項１〜６のいずれか一項に記載の方法。
【請求項８】
前記定量位相画像を生成する工程が、
対物レンズの光軸方向に前記細胞を配置する工程と、
前記細胞に照明光を照射する工程と、
前記光軸に沿って互いに間隔Δｚ離れた複数の位置のそれぞれに前記対物レンズの焦点を配置して前記細胞からの光を検出する工程と、
検出した前記光に基づいて、前記複数の位置のそれぞれに対応する光強度分布データを生成する工程と、
前記光強度分布データに基づいて前記定量位相画像を生成する工程と、を含み、
前記間隔Δｚは、前記対物レンズの開口数ＮＡ、前記照明光の波長λ、及び、前記対物レンズと前記細胞との間の屈折率ｎの少なくとも１つに基づいて設定される、請求項１〜７のいずれか一項に記載の方法。
【請求項９】
前記間隔Δｚは、更にパラメータｋに基づいて設定され、前記ｋは、生成する前記定量位相画像において位相を復元する空間周波数を示す、請求項８に記載の方法。
【請求項１０】
前記Δｚ、前記ＮＡ、前記λ、前記ｎ、前記ｋは下記式（１００）を満たし、前記ｋの値は１〜２５である、請求項９に記載の方法。
【数１】
2020188717.tif 000009
【請求項１１】
多孔膜上に配置された細胞の数、形態又は形状を測定するための定量位相画像を生成する定量位相画像生成装置であって、
顕微鏡と、装置制御部と、検出器と、画像生成部とを備え、
前記顕微鏡は、多孔膜上に配置された細胞が配置されるステージと、前記細胞に照明光を照射する照明光学系と、対物レンズとを有し、
前記装置制御部は、前記対物レンズの光軸に沿って互いに間隔Δｚ離れた複数の位置のそれぞれに前記対物レンズの焦点が配置されるように前記焦点の位置を調整するものであり、
前記検出器は、前記照明光学系から照射された照明光により、前記細胞から生じる光を検出するものであり、
前記画像生成部は、前記検出器が検出した前記光に基づいて、前記複数の位置に対応する光強度分布データを生成し、前記光強度分布データに基づいて、定量位相画像を生成するものであり、
前記装置制御部は、前記対物レンズの開口数ＮＡ、前記照明光の波長λ、及び、前記対物レンズと前記細胞との間の屈折率ｎの少なくとも１つに基づいて前記間隔Δｚを設定するものである、定量位相画像生成装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、細胞の数、形態又は形状を測定する方法及び装置に関する。
【背景技術】
【０００２】
定量位相イメージングは、測定試料の物理的／幾何学的な厚みの変化と屈折率の変化を積とする位相を測定することができる観察技術である。定量位相イメージングでは、試料を保持する容器、試料を含む液体の液面等により生じる光の干渉により、定量位相画像にノイズが生じる場合がある。特許文献１には、このようなノイズを減少させることを目的とした、デジタルホログラフィック顕微鏡用観察容器が記載されている。
【０００３】
ところで、Ｏｒｇａｎ−ｏｎ−ａ−Ｃｈｉｐ（オーガンチップ）とは、複数のヒト臓器細胞を一つのチップ上に集約したデバイスであり、薬剤スクリーニングや、ヒトに薬物を投与したときの薬効、薬物動態、安全性などの評価等に利用される。オーガンチップの一形態としては、流路の中央に多孔膜を配置して多孔膜上で細胞を培養し、その上下に異なる培地や試薬を流す構成を有するものがある。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特表２０１１−５２５２５２号公報
【発明の概要】
【０００５】
一実施形態に係る方法は、多孔膜上に配置された細胞の数、形態又は形状を測定する方法であって、非干渉方式の定量位相イメージングにより前記細胞の定量位相画像を生成する工程と、前記定量位相画像に基づいて前記細胞の数、形態又は形状を測定する工程とを含む方法である。
【０００６】
一実施形態に係る定量位相画像生成装置は、多孔膜上に配置された細胞の数、形態又は形状を測定するための定量位相画像を生成するものであり、顕微鏡と、装置制御部と、検出器と、画像生成部とを備える。前記顕微鏡は、多孔膜上に配置された細胞が配置されるステージと、前記細胞に照明光を照射する照明光学系と、対物レンズとを有する。前記装置制御部は、前記対物レンズの光軸に沿って互いに間隔Δｚ離れた複数の位置のそれぞれに前記対物レンズの焦点が配置されるように前記焦点の位置を調整するものである。前記検出器は、前記照明光学系から照射された照明光により、前記細胞から生じる光を検出するものである。前記画像生成部は、前記検出器が検出した前記光に基づいて、前記複数の位置に対応する光強度分布データを生成し、前記光強度分布データに基づいて、定量位相画像を生成するものである。前記装置制御部は、前記対物レンズの開口数ＮＡ、前記照明光の波長λ、及び、前記対物レンズと前記細胞との間の屈折率ｎの少なくとも１つに基づいて前記間隔Δｚを設定するものである。
【図面の簡単な説明】
【０００７】
【図１】定量位相画像生成装置の構成を示す概念図である。
【図２】定量位相画像生成装置の情報処理部の構成を示す概念図である。
【図３】対象物による照明光の位相の変化により、計測される光の強度が変化することを示す概念図である。
【図４】計測面を説明する概念図である。
【図５】焦点間隔を説明する概念図である。
【図６】細胞培養容器及び細胞培養用インサートを用いた細胞培養を説明する模式断面図である。
【図７】オーガンチップの一形態の構成を説明する模式断面図である。
【図８】（ａ）及び（ｂ）は、実験例１で撮影した、通常の細胞培養プレート上に配置された細胞の顕微鏡写真である。（ａ）は位相差画像であり、（ｂ）は定量位相画像である。（ｃ）は、（ａ）及び（ｂ）の各画像に示す矢印に沿って、横軸に距離（μｍ）、縦軸に輝度（相対値）をそれぞれプロットしたグラフである。
【図９】（ａ）及び（ｂ）は、実験例１で撮影した、多孔膜上に配置された細胞の顕微鏡写真である。（ａ）は位相差画像であり、（ｂ）は定量位相画像である。（ｃ）は、（ａ）及び（ｂ）の各画像に示す矢印に沿って、横軸に距離（μｍ）、縦軸に輝度（相対値）をそれぞれプロットしたグラフである。
【図１０】実験例２における計測面を説明する概念図である。
【発明を実施するための形態】
【０００８】
以下、場合により図面を参照しつつ、本発明の実施形態について詳細に説明する。なお、図面中、同一又は相当部分には同一又は対応する符号を付し、重複する説明は省略する。なお、各図における寸法比は、説明のため誇張している部分があり、必ずしも実際の寸法比とは一致しない。
【０００９】
［多孔膜上に配置された細胞の数、形態又は形状を測定する方法］
一実施形態において、本発明は、多孔膜上に配置された細胞の数、形態又は形状を測定する方法であって、非干渉方式の定量位相イメージングにより前記細胞の定量位相画像を生成する工程と、前記定量位相画像に基づいて前記細胞の数、形態又は形状を測定する工程と、を含む、方法を提供する。
【００１０】
多孔膜上で細胞を非侵襲的に観察する需要がある。例えば、マルチウェルプレート等の細胞培養容器のウェルの内部に細胞培養用インサート（カルチャーインサート）を配置し、インサートの内部で細胞を培養することが広く行われている。また、インサートの外部（細胞培養容器のウェル）でインサートの内部とは異なる細胞を培養する場合もある。細胞培養用インサートの底部には多孔膜が配置されており、インサートの内部と外部（細胞培養容器のウェルの内部）との間で物質交換が可能となっている。
【００１１】
図６は、一例として、細胞培養容器及び細胞培養用インサートを用いた細胞培養を説明する模式断面図である。図６に示すように、細胞培養容器６００は、基板６０１及び基板６０１の一方面上に配置された側壁６０２から構成される空間（ウェル）を有している。基板６０１と側壁６０２は密着しており、ウェル内に培地等の液体を保持することができる。そして、ウェルの内部には、細胞培養用インサート６１０が配置されている。
【００１２】
細胞培養用インサート６１０の底部には、多孔膜６１１が配置されている。細胞培養用インサート６１０の内部には培地６２０が保持されている。また、細胞培養容器６００のウェルの内部にも培地６２１が保持されている。
【００１３】
また、多孔膜６１１の上には細胞Ｃが配置されている。細胞Ｃは接着細胞であってもよいし、浮遊細胞であってもよい。浮遊細胞であっても、培養容器を静置すれば多孔膜６１１の上に沈み、観察することができる。
【００１４】
図７は、一例として、オーガンチップの構成を説明する模式断面図である。図７に示すように、オーガンチップ７００は、一方面がオーガンチップ７００の天井面７１１となる第１の基板７１０と、一方面がオーガンチップ７００の底面７２１となる第２の基板７２０と、第１の基板７１０と第２の基板７２０との間に、第１の基板７１０及び第２の基板７２０に略平行に配置された多孔膜６１１と、第１の基板７１０、第２の基板７２０及び多孔膜６１１とを支持する支持体（図示せず）とを有する。第１の基板７１０と多孔膜６１１との間の空間は第１の流路７４０を形成し、第２の基板７２０と多孔膜６１１との間の空間は第２の流路７５０を形成する。多孔膜６１１のうち、第１の基板７１０に対向する面上では、細胞Ｃが培養される。また、第１の流路７４０には培地６２０が流される。また、第２の流路７５０には培地６２１が流される。培地６２０又は培地６２１には試薬等が添加されてもよい。
【００１５】
図６や図７に例示されるような多孔膜６１１上の細胞Ｃを観察する場合、従来の光学顕微鏡を用いて、細胞Ｃの明視野画像、位相差画像、微分干渉画像等を撮影しても、細胞Ｃの形状を鮮明に観察することが困難である。
【００１６】
これは、観察対象である細胞Ｃ以外に、多孔膜６１１、培地６２０及び６２１、細胞培養容器６１０の底面を構成する基板６０１、オーガンチップを構成する第１の基板７１０、第２の基板７２０等が、光学顕微鏡の光路に入ってしまうためである。
【００１７】
また、定量位相イメージングには光干渉を利用する方式（干渉方式）と、非干渉方式が存在する。図６や図７に例示されるような多孔膜６１１上の細胞Ｃを観察する場合に、干渉方式の定量位相イメージングを行っても、やはり細胞Ｃの形状を鮮明に観察することが困難である。また、特許文献１に記載されているような容器を用いても、干渉方式の定量位相イメージングを行えば、やはり細胞Ｃの形状を鮮明に観察することが困難である。
【００１８】
これは、上述したように、多孔膜６１１、培地６２０及び６２１、基板６０１多孔膜６１１、培地６２０及び６２１、基板６０１、基板７１０、基板７２０等が顕微鏡の光路に入って迷光となり、定量位相画像にノイズが生じるためである。特に、多孔膜６１１による干渉により、細胞Ｃの画像が不鮮明になる傾向にある。
【００１９】
これに対し、実施例において後述するように、発明者らは、非干渉方式の定量位相イメージングによれば、多孔膜上に配置された細胞の形状を正確に測定することができることを見出し、本発明を完成させた。
【００２０】
本明細書において細胞の形状を正確に測定するとは、細胞の形状を明瞭に観察するといいかえることもできる。ここで、「正確に」とは、例えば、明視野画像、位相差画像、微分干渉画像等と比較して、実際の細胞の形態や形状の測定精度が高いことを意味する。また、細胞の形態とは大きさや厚み等を意味する。また、細胞の形状とは、細胞の外面形状、細胞内部の細胞小器官の形状等を意味する。本実施形態の方法によれば、特に細胞の外面形状を明瞭に観察することができる。このため、細胞の数を正確に測定することもできる。ここで、細胞の形態や形状を正確に測定するとは、例えば、視野内に存在する細胞の高さや幅等の大きさ、細胞の位置等を正確にマッピングすることであってよい。本実施形態の方法によれば、多孔膜上で細胞を非侵襲的に観察することができる。
【００２１】
多孔膜としては、細胞培養等に通常用いられるものであれば特に限定されず、例えば、平均ポアサイズ０．１〜１０μｍの多孔膜が挙げられる。ポア密度は１０５〜１０９個／ｃｍ２程度であってよい。また、多孔膜の厚さは、０．０１〜１．０ｍｍ程度であってよい。また、多孔膜の材質としては、例えば、ポリカーボネート、ポリエステル、ポリエチレンテレフタレート、ポリテトラフルオロエチレン等が挙げられる。
【００２２】
多孔膜は、細胞培養容器内に配置した細胞培養用インサートの底面に配置されているものであってもよい。また、細胞培養容器はオーガンチップを構成していてもよい。すなわち、多孔膜は、オーガンチップ内に配置されているものであってもよい。
【００２３】
本実施形態の方法は、光干渉の影響を実質的に受けないことから、細胞培養容器ごと、あるいは、オーガンチップごと顕微鏡のステージに載せて細胞を観察することができる。また、細胞培養容器の材質は、通常の細胞培養等に通常用いられる光透過性の高いものであれば特に限定されず、例えば、ポリスチレン、ポリプロピレン、ガラス、ポリエステル、ポリジメチルポリシロキサン（ＰＤＭＳ）等が挙げられる。
【００２４】
また、細胞培養容器の表面には、細胞の接着性を高めるための表面処理が施されていてもよい。このような表面処理としては、プラズマ処理等の物理的な浸水化処理；ゼラチン、フィブロネクチン、コラーゲン、細胞外マトリックス等の細胞接着物質の塗布等が挙げられる。
【００２５】
細胞培養容器の形状も特に限定されず、シャーレ；６ウェル、１２ウェル、２４ウェル、４８ウェル、９６ウェル等のウェルプレート；流体デバイス等が挙げられる。
【００２６】
上述したように、本実施形態の方法は、非干渉方式の定量位相イメージングにより細胞の定量位相画像を生成する工程と、前記定量位相画像に基づいて前記細胞の数、形態又は形状を測定する工程とを含む。
【００２７】
定量位相画像を生成する工程は、対物レンズの光軸方向に細胞を配置する工程と、前記細胞に照明光を照射する工程と、前記光軸に沿って互いに間隔Δｚ離れた複数の位置のそれぞれに前記対物レンズの焦点を配置して前記細胞からの光を検出する工程と、検出した前記光に基づいて、前記複数の位置のそれぞれに対応する光強度分布データを生成する工程と、前記光強度分布データに基づいて前記定量位相画像を生成する工程と、を含み、前記間隔Δｚは、前記対物レンズの開口数ＮＡ、前記照明光の波長λ、及び、前記対物レンズと前記細胞との間の屈折率ｎの少なくとも１つに基づいて設定されてもよい。
【００２８】
また、前記間隔Δｚは、更にパラメータｋに基づいて設定され、前記ｋは、生成する前記定量位相画像において位相を復元する空間周波数を示す値であってもよい。
【００２９】
一態様において、前記間隔Δｚは、下記式（１００）に基づいて設定されてもよい。以下、各工程について説明する。
【００３０】
【数１】
2020188717.tif 000003
［式（１００）中、ＮＡは前記対物レンズの開口数を表し、λは前記照明光の波長を表し、ｎは前記対物レンズと前記細胞との間の屈折率を表し、ｋは前記定量位相画像の生成において位相を復元する空間周波数を示すパラメータを表し、前記パラメータｋの値は１〜２５である。］
【００３１】
定量位相画像とは、対象物における厚さの変化と屈折率の変化を積とする位相を表す（画像化した）画像である。既存の位相差画像や微分干渉画像は、位相に振幅が混ざる干渉画像であるため、対象物の位相を完全に定量化することができないが（振幅が無視でき、かつ位相が小さい条件のみ定量化することができる。）、定量位相画像は、対象物の位相が２倍、３倍、４倍等の様々な値に変化しても、定量的な位相を表すことができるため、対象物の位相を定量化した画像ということができる。
【００３２】
定量位相イメージングは、対象試料を染色する必要がなく、本来無色透明で見えにくい試料に明暗を付けて可視化できる利点に加え、生体細胞構造を反映した高さ情報を測定することができる。また、細胞の屈折率を固定値とすれば、細胞厚が低いものと高いものを比較でき、三次元構造を俯瞰することができる。このため、定量位相画像を培養細胞の品質管理等に応用することもできる。
【００３３】
図１は、本実施形態の方法に用いることができる定量位相画像生成装置の構成を示す概念図である。定量位相画像生成装置１は、顕微鏡、すなわち定量位相顕微鏡であり、顕微鏡（顕微鏡本体部）１００と、情報処理部４０とを備える。
【００３４】
顕微鏡１００は、対象物（細胞）Ｓに照明光を照射する透過照明光学系１０と、ステージ８と、結像光学系７と、落射照明光学系１１０と、検出部９とを備える。結像光学系７は、対物光学系２０と、フィルターキューブ１２０と、リレー光学系３０とを備える。なお、情報処理部４０の機能は顕微鏡１００から物理的に離れたＰＣやサーバ等の電子計算機等に配置してもよい。
【００３５】
透過照明光学系１０は、光源１１と、レンズ１２と、バンドパスフィルター１３と、視野絞り１４と、レンズ１５と、開口絞り１６と、集光レンズ１７とを備える。対物光学系２０は、複数の対物レンズ２１を備える。対物レンズ２１としては、対物レンズ２１ａ，２１ｂ，２１ｃの複数の対物レンズがターレット等により保持されており、いずれかを適宜選択して使用する。リレー光学系３０は、結像レンズ３１と、ビームスプリッター３２と、ミラー３３ａ，３３ｂ，３３ｃと、レンズ３４ａ，３４ｂ，３４ｃと、接眼レンズ３５とを備える。接眼レンズ３５はユーザの眼Ｅを近づけて覗き込むことができるように構成されている。
【００３６】
なお、交換可能な対物レンズ２１の数は特に限定されない。また、顕微鏡１００に含まれる光学系の態様は、対象物Ｓの所望の面の画像を取得することができれば特に限定されない。
【００３７】
図１において、顕微鏡１００の光学系の光軸を一点鎖線Ｌ１で示し、光源１１からの照明光（光束）を二点鎖線Ｌ２で模式的に示した。対物レンズ２１のうち、対象物Ｓの観察に使用されている対物レンズを対物レンズ２１ａとすると、後に詳述する対物レンズ２１ａの光軸Ｌｏは、透過照明光学系１０の光軸Ｌ１と一致しているため、これを図１においてＬ１（Ｌｏ）と示した。検出部９で検出された光の検出信号は、実線矢印Ａ１で示されるように、情報処理部４０に出力される。また、情報処理部４０は顕微鏡１００を制御するが、その制御を実線矢印Ａ２で模式的に示した。一点鎖線Ｌ１，Ｌ２、実線矢印Ａ１，Ａ２は以下の実施形態でも同様のものを示す。
【００３８】
光源１１は、ハロゲンランプ等の非コヒーレント光源を含み、情報処理部４０による制御のもと、対象物Ｓに照射する照明光Ｌ２として非コヒーレント光を出射する。本実施形態の透過照明光学系１０では、光源１１から出射された非コヒーレント光が、後述する開口絞り１６を含む透過照明光学系１０の各光学部材により、位相復元が可能な光軸に略垂直な波面を有する光とされて対象物Ｓに照射される構成となっている。
【００３９】
「位相復元」とは、後に詳述するが、計測した対象物Ｓからの光の強度から、強度輸送方程式を用いて対象物Ｓの構造を反映した位相値を算出することである。照明光Ｌ２は、対象物Ｓの構造を反映した位相の計測が可能であればどのような波長の光でも構わないが、可視光で見えづらい対象物を位相計測により可視化できる観点では、そのまま可視光を使用することが簡素で好ましい。
【００４０】
なお、対象物Ｓに照射する照明光Ｌ２は、可視光でなくてもよく、紫外光や赤外光であってもよい。さらに、照明光Ｌ２の波面は透過照明光学系１０の光軸に略垂直でなくとも、波面の形状が既知であればよく、例えば、照明光Ｌ２の波面は略球面であってもよい。また、光源１１は、パルスレーザや連続発振（ＣＷ）レーザ等のコヒーレント光源を含み、照明光Ｌ２としてコヒーレント光を出射してもよい。
【００４１】
光源１１から出射した照明光Ｌ２は、レンズ１２に入射する。レンズ１２に入射した照明光Ｌ２は、レンズ１２により屈折され略平行な光となってレンズ１２を出射し、バンドパスフィルター１３に入射する。バンドパスフィルター１３に入射した照明光Ｌ２は、所望の波長範囲の波長成分の光のみがバンドパスフィルター１３により透過されて、視野絞り１４に入射する。ここで、バンドパスフィルター１３により透過される光の波長範囲は、軸上色収差によって発生する、位相を復元する際の、位相の実測値と実際の位相との誤差（後述の位相復元誤差）が大きくなり過ぎないように適宜設定する。バンドパスフィルター１３は透過照明光学系１０の光路に挿脱可能に構成されており、情報処理部４０の制御のもと、適宜光路外の位置Ｐ１に退避可能に構成されている。バンドパスフィルター１３を透過照明光学系１０の光路に挿脱可能するバンドパスフィルター駆動部は、簡略化のため、図示を省略している。なお、バンドパスフィルター駆動部は必ずしも必要な構成でなく、例えば、ユーザが手動で、バンドパスフィルター１３を透過照明光学系１０の光路に挿脱してもよい。
【００４２】
なお、光源１１としてレーザ光等の波長範囲が狭い光源を用いる際は、バンドパスフィルター１３を光路に配置する必要は無い。また、バンドパスフィルター１３を用いずに、対象物Ｓに対して結像側に配置されたフィルターキューブ１２０のフィルターにより、検出部９で検出する波長範囲を制限してもよい。
【００４３】
視野絞り１４に入射した照明光Ｌ２は、その光束径が調節されて視野絞り１４を出射し、レンズ１５に入射する。レンズ１５に入射した照明光Ｌ２は、レンズ１５により収束されてレンズ１５を出射し、開口絞り１６に入射する。開口絞り１６に入射した照明光Ｌ２は、その波面が球面状になるように変換され、開口絞り１６を出射し、集光レンズ１７に入射する。集光レンズ１７に入射した照明光Ｌ２は、集光レンズ１７により屈折され、対象物Ｓに照射される際に、透過照明光学系１０の光軸に略垂直な波面を有する光となり、対象物Ｓに照射される。
【００４４】
対象物Ｓは、ステージ８上に配置される。対象物Ｓは、特に限定されないが、照明光が対象物Ｓに照射され透過する際に、振幅の変化の割合は比較的小さく、一方、位相の変化の割合は大きい場合に、本実施形態に係る位相計測の効果がより顕著に得られるので好ましい。このような物体を位相物体と呼ぶ。すなわち、位相物体は、透過光においては光の強度を変えず、位相のみを変える無色透明な物体である。培養細胞等の細胞は位相物体である。
【００４５】
ステージ８は、対物レンズ２１ａの光軸、及び当該光軸に垂直な軸に沿って移動可能に構成されている。ここで、対物レンズ２１ａの光軸とは、対物レンズ２１ａの対象物Ｓ側のレンズ面の光学中心と対物レンズ２１ａの焦点を通る直線で示される光軸Ｌｏ（図５参照）を指す。ステージ８は、モータ等の移動装置による電動駆動により、対象物Ｓを撮像する際に、対物レンズ２１ａの光軸Ｌｏに沿った移動が可能に構成されている。図５では、対物レンズ２１ａの光軸Ｌｏが、透過照明光学系１０の光軸Ｌ１と一致している点を、Ｌｏ（Ｌ１）と示した。
【００４６】
以下の実施形態では、図１に座標軸９００として示すように、顕微鏡１００の光軸Ｌ１、すなわち対物レンズ２１ａの光軸Ｌｏに平行にｚ軸をとり、ｚ軸に垂直であって紙面に平行にｘ軸を取り、ｘ軸及びｚ軸に垂直にｙ軸をとる。情報処理部４０は、図示しないステージ駆動部を駆動して、対象物Ｓにおける観察対象部分またはその近傍に、透過照明光学系１０および対物レンズ２１ａの光軸Ｌｏ（Ｌ１）が一致するようにステージ８を移動する。
【００４７】
対物光学系２０は、対物レンズ２１を含んで構成される。本実施形態では、対物レンズ２１としては、複数の対物レンズ２１ａ，２１ｂ，２１ｃを備える。後述するように、定量位相画像生成装置１では、開口数ＮＡ等の異なる対物レンズ２１ａ，２１ｂ，２１ｃを交換して対象物Ｓを撮像して、それぞれの対物レンズ２１により得たデータに基づいて位相計測を行う際、位相復元パラメータを揃えることによって、得られる位相値のばらつきを抑制するように構成されている。
【００４８】
対物レンズ２１を変更した際に計測される位相値の変化が大きいと、異なる対物レンズ２１を用いて取得したデータの間の比較が難しくなる。また、異なる対物レンズ２１を用いて取得したデータのうち、どの対物レンズ２１を用いて取得したデータが位相を復元する際の精度が高いのか判断が難しい等の理由により、一連の測定に係る作業に負担が生じる。
【００４９】
落射照明光学系１１０は、水銀ランプ等の光源１１１や、フィルターキューブ１２０等を含んで構成され、蛍光観察等のための励起光を出射する。光源１１１から出射した光は、フィルターキューブ１２０に入射する。
【００５０】
フィルターキューブ１２０は、励起フィルター１２１と、ダイクロイックミラー１２２と、吸収フィルター１２３とを含んで構成され、光源１１１から入射した光のうち、所定波長の光を励起光として対象物Ｓへと反射する一方、対象物Ｓから発せられる蛍光をリレー光学系３０へと透過させる。フィルターキューブ１２０は、落射光学系１１０の光路に挿脱可能に構成されており、情報処理部４０の制御のもと、光路外の位置Ｐ２に適宜退避できるように構成されている。フィルターキューブ１２０を落射光学系１１０の光路に挿脱するフィルターキューブ駆動部は、簡略化のため、図示を省略している。なお、フィルターキューブ駆動部は必ずしも必要な構成でなく、例えば、ユーザが手動で、フィルターキューブ１２０を落射光学系１１０の光路に挿脱してもよい。
【００５１】
励起フィルター１２１は、光源１１１から入射した光のうち、所定の波長範囲の光を励起光として透過する。励起フィルター１２１の透過波長領域は、ダイクロイックミラー１２２の反射波長領域に一致する。ダイクロイックミラー１２２は、励起フィルター１２１を透過した励起光を対象物Ｓに向けて反射し、対物レンズ２１ａから入射した対象物Ｓからの蛍光を吸収フィルター１２３へと透過する。吸収フィルター１２３は、ダイクロイックミラー１２２から入射した光のうち、対象物Ｓや光学系からの不要な散乱光や迷光を吸収し、必要な光、すなわち、対象物Ｓからの蛍光のみをリレー光学系３０の結像レンズ３１に出射する。
【００５２】
リレー光学系３０の結像レンズ３１は、フィルターキューブ１２０から入射した光を検出部９に結像するように屈折させてビームスプリッター３２に出射する。ビームスプリッター３２は、フィルターキューブ１２０から入射した光の一部を検出部９へと反射し、残りは透過させてミラー３３ａに出射する。ミラー３３ａで反射した光は、レンズ３４ａ、ミラー３３ｂ、レンズ３４ｂ、レンズ３４ｃ及びミラー３３ｃの順で、ミラーでの反射又はレンズでの屈折を経て接眼レンズ３５に入射する。接眼レンズ３５に入射した光は接眼レンズ３５により屈折されてユーザの眼Ｅに入射して知覚される。
【００５３】
検出部９は、ＣＣＤやＣＭＯＳ等の撮像素子の検出器を含んで構成され、リレー光学系３０のビームスプリッター３２により反射された光を検出する。撮像素子はエリアセンサであってもよいし、リニアセンサであってもよいし、ポイントセンサであってもよい。検出した光に対応する検出信号は不図示のＡ／Ｄ変換器等により適宜Ａ／Ｄ変換されて情報処理部４０に出力される。いいかえると、検出部９は、対象物Ｓの像を撮像する。
【００５４】
なお、定量位相画像生成装置１では、対象物Ｓの構造を反映した位相の計測に際し、検出部９が、透過照明光学系１０からの照明光Ｌ２が照射された対象物Ｓの像を撮像する構成にしたが、落射照明光学系１１０からの照明光が照射された対象物Ｓの像を撮像する構成にしてもよい。この場合、落射照明光学系１１０からの照明光は、対物レンズ２１ａ（対物光学系２０）を介して対象物Ｓに照射され、対象物Ｓで反射及び散乱した光に基づく対象物Ｓの像を検出部９で撮像する。
【００５５】
図２は、情報処理部４０の構成を示した図である。情報処理部４０は、入力部４１と、表示部４２と、通信部４３と、記憶部４４と、制御部５０とを備える。制御部５０は、装置制御部５１と、解析部５２とを備える。これらはバス等により接続されているが、簡略化のため省略している。装置制御部５１は、最適条件計算部５１１を備える。解析部５２は、位相復元処理部５２１と、画像構築部５２２と、画像解析部５２３とを備える。
【００５６】
入力部４１は、キーボード、マウス及び／又はタッチパネル等の入力装置により構成され、顕微鏡１００による対象物Ｓの撮像や解析部５２による当該撮像により得られたデータの解析に必要な情報等を含む入力データを受け付ける。入力部４１は、受け付けた入力データを、後述の記憶部４４に適宜記憶させる。なお、入力データは、後述の通信部４３を介して取得してもよい。
【００５７】
表示部４２は、液晶モニタ等の表示装置により構成され、顕微鏡１００による撮像の条件、当該撮像により得られたデータに基づいて解析部５２が生成した解析結果及び定量位相画像等を表示する。
【００５８】
通信部４３は、インターネット等の通信網を利用して通信を行う通信装置により構成され、顕微鏡１００による撮像の条件、当該撮像により得られたデータに基づいて解析部５２が生成した解析結果及び定量位相画像等を送信したり、適宜必要なデータを送受信したりする。
【００５９】
記憶部４４は、不揮発性メモリ等の記憶装置により構成され、制御部５０に処理を行わせるプログラム、顕微鏡１００による撮像に必要なデータ、当該撮像により得られたデータ及び当該データに基づいて解析部５２が生成した解析結果及び定量位相画像等を記憶する。
【００６０】
制御部５０は、ＣＰＵ等のマイクロプロセッサ等を含むプロセッサにより構成され、定量位相画像生成装置１を制御する動作の主体として機能する。すなわち、記憶部４４等に搭載されているプログラムを実行することにより、顕微鏡１００による対象物Ｓの撮像を行う装置制御処理、当該撮像により得られたデータの位相復元処理等の解析処理、及び出力処理等の各種処理を行う。
【００６１】
制御部５０の装置制御部５１は、ユーザにより入力された入力データに基づいて顕微鏡１００による対象物Ｓの撮像に必要なパラメータを適宜算出するとともに、取得した当該パラメータ及び入力データ等に基づいて顕微鏡１００の各部の動作を制御する。
【００６２】
最適条件計算部５１１は、顕微鏡１００が対象物Ｓにおける複数の位置のそれぞれに対物レンズ２１ａの焦点を順次配置して検出部９により対象物Ｓからの光を検出する際に設定するための、当該複数の位置の間隔Δｚを算出する。間隔Δｚは透過照明光学系１０の光軸方向（対物レンズの光軸方向）に設定される。この間隔Δｚを焦点間隔と呼ぶ。焦点間隔Δｚは、精度よく位相計測を行うために重要なパラメータである。以下、この点に関し、強度輸送方程式を利用した位相計測の方法に基づいて説明する。
【００６３】
図３は、位相物体を透過した光の強度分布が変化する点を説明する概念図である。対象物Ｓは、細胞等の位相物体であり、対象物Ｓには、波面Ｗ１，Ｗ２，Ｗ３により模式的に示されているように、互いに平行で、透過照明光学系１０／対物レンズの光軸Ｌｏ（Ｌ１）に略垂直な波面を有する平面波が照射されている。平面波の進行方向は矢印Ａ３で示した。対象物Ｓを透過した光は、振幅は大きく変化しないものの、位相が変化することにより、等位相面、すなわち波面が変化する（波面Ｗ４）。
【００６４】
ホイヘンス＝フレネルの原理で、曲線状の波面を構成する光は、実線矢印Ａ４で示されるように計測面ｉ１に到達する。ここで、後に詳述するように、一例として対物レンズ２１ａの焦点が配置される対象物Ｓにおける特定の面（いいかえると、対物レンズ２１ａの合焦面）を計測面ｉ１と称する。計測面ｉ１は、対物レンズ２１ａの光軸に実質的に垂直な面である。
【００６５】
図３では、計測面ｉ１は平面波の波面Ｗ１〜Ｗ３と実質的に平行である。計測面ｉ１上では、上述のように対象物Ｓによる波面の変化に伴い所定の光の強度分布（以下、「光強度分布」という場合がある。）が形成される。対象物Ｓによる透過光の位相の変化と光強度分布との関係に基づいて、透過光の強度を解析することにより、位相の変化を計測する方法が提案されている。本実施形態では、以下に説明する強度輸送方程式を用いた方法を利用する。
【００６６】
伝播する波動における、強度Ｉと位相φの関係は、強度輸送方程式により記述される。強度輸送方程式は、以下の式（１）及び（２）による連立方程式である。詳細は、Paganin D. and Nugent K.A., Noninterferometric Phase Imaging with Partially Coherent Light, Physical Review Letters, Volume 80, 2586-2589, 1998.を参照されたい。なお、この文献は、特段の定めのない限り、参照としてここに組み込まれる。
【００６７】
【数２】
2020188717.tif 000004
【００６８】
【数３】
2020188717.tif 000005
【００６９】
ここで、∇の添字であるｘｙは、ＸＹ平面（光の伝搬方向に垂直な平面、すなわち本実施形態ではｚ軸に垂直な平面。ここで、ｚ軸は透過照明光学系１０／対物レンズの光軸Ｌｏ（Ｌ１）と一致する。）を示す。すなわち、式（２）の左辺に現れる∇ｘｙφは、ＸＹ平面における位相φの勾配（ｇｒａｄｉｅｎｔ）を示す。式（１）の左辺と式（２）の右辺に現れるΦは、式（１）をポアソン方程式の形にして計算しやすいように導入された関数であり、検出したＸＹ平面における強度Ｉのｚに関する微分係数の分布を得て、ポアソン方程式（１）を解くことにより導出できる。式（１）を解いて得た関数Φと光強度分布から、式（２）を解いて位相分布φを算出することができる。
【００７０】
図４は、位相分布φの算出のために行う対象物Ｓからの光の検出方法を説明する図である。まず、対象物Ｓの計測面ｉ１における光強度分布及び計測面ｉ１におけるｚに関する微分係数分布を得る必要がある。このため、装置制御部５１は、図４で示すように、対物レンズ２１ａの焦点が所望の計測面ｉ１に含まれるように設定して、結像光学系７により検出部９の検出面に対象物９の像を結像させて対象物Ｓからの光を検出部９に検出させる。検出部９の各画素からの検出信号は、制御部５０（図２）に入力され、位相復元処理部５２１により各画素の位置と、各画素で検出された光強度とが対応付けられた光強度分布データが生成される。
【００７１】
光強度分布データは、ある座標ｚの値に対応する計測面上での座標（ｘ、ｙ）の位置に対応する、検出部９における各画素で検出された光の強度分布を表すデータである。光強度分布データは、当該座標ｘ、ｙ、ｚの値に対応する光強度を示すデータであり、ルックアップテーブルの形式で構築されている。例えば、あるｚの値に対応する計測面上での光強度分布を色又は階調で区別して二次元にマッピングすることにより、当該ｚの位置における光の強度分布を示す画像（以下、「光強度画像」という場合がある。）を作成することができる。なお、所定の座標ｘ、ｙ、ｚの値に対応する光強度の値を取り出すことができれば、光強度分布データのデータ構造は特に限定されず、他の既存のデータ構造であってもよい。
【００７２】
計測面ｉ１の位置（いいかえると、対物レンズ２１ａの焦点位置）は、検出部９により検出される対象物Ｓからの光の光強度のコントラストに基づいた位置に設定されることが好ましい。装置制御部５１は、事前に（本実施形態に係る位相計測の前に）取得した対象物Ｓの三次元の光強度分布データから算出したパラメータ（以下に示す分散ｖｚ等の対象物Ｓからの光の光強度のコントラストを示すパラメータ）に基づいて、計測面ｉ１の位置を設定することができる。上記三次元の光強度分布データを事前に撮像により取得する場合、この光強度分布データの取得のための撮像の際の対象物Ｓにおける焦点の位置は特に限定されず、適宜設定すればよい。
【００７３】
装置制御部５１は、事前に取得した三次元の光強度分布データについて、各ｚの値に対応した、ｘ方向及びｙ方向についての二次元の光強度分布をＩｚ（ｘ、ｙ）としたとき、各ｚの値に対応する光強度分布の分散値ｖｚを以下の式（１１）により算出する。
【００７４】
【数４】
2020188717.tif 000006
【００７５】
ここで、Ｎｘは事前に取得した光強度分布データのｘ方向のピクセル数、Ｎｙは事前に取得した光強度分布データのｙ方向のピクセル数、Ｉｚの上にバーが示されたものは各ｚの値に対応する、当該光強度分布データ上の全ての（ｘ，ｙ）についてのＩｚ（ｘ，ｙ）の算術平均等の平均値である。装置制御部５１は、算出した各ｚの値に対応する分散ｖｚに基づいて計測面ｉ１の位置、すなわち、対物レンズ２１ａの焦点位置を設定する。例えば、装置制御部５１は、算出した分散ｖｚのうち、極小となる分散ｖｚに対応するｚ方向の位置に計測面ｉ１の位置を設定する。いいかえると、装置制御部５１は、計測面ｉ１の位置をコントラストが低い光強度分布データに対応するｚの値に設置する。
【００７６】
市販のウェルプレートの場合、細胞（位相物体）内をｚ軸方向に異なる複数のｚ位置に順次移動させて、前記複数のｚ位置の各々で光強度分布データを取得して、極小となる分散ｖｚに対応するｚ方向の位置を算出して、計測面ｉ１の位置を設定できる。
【００７７】
トランズウェルの場合、多孔膜６１１の上にいる細胞に対して、計測面ｉ１の位置を設定する必要がある。光強度分布データから、空間周波数が低い背景成分を除去して、ポア形状の円形度でフィルタリングした画像を生成して、細胞がしない背景領域に存在するポアに注目して、極小となる分散ｖｚに対応するｚ方向の位置を算出した後、多孔膜の厚さ平均値分をオフセットさせたｚ位置を、計測面ｉ１の位置に設定することが好ましい。
【００７８】
多孔膜６１１の厚さはｘｙ平面に一定ではない（平坦でない）ために、ポアサイズとポア密度に依存して、局所的な厚さ変化が大きくなる傾向が確認されて、撮影視野内でも計測面ｉ１の位置が一つに決定されない場合がある。
【００７９】
局所的な膜厚ばらつきが大きい場合には、撮影視野サイズ内の画素群を一度、四隅部及び中心部の画素領域に細かく分割して、分割した各々の画素領域のポアを抽出する。各々画素領域で算出した、極小となる分散ｖｚに対応するｚ位置データ群の平均値を採用して、多孔膜の厚さ分をオフセットさせたｚ位置を、計測面ｉ１の位置に設定できる。
【００８０】
装置制御部５１は、最適条件計算部５１１が算出した焦点間隔Δｚに基づいて、対物レンズ２１ａの焦点位置が、計測面ｉ１上の位置からｚ軸に沿って−Δｚ及び＋Δｚの距離ずれたそれぞれの位置となるように計測面ｉ２及びｉ３上に当該焦点位置を順次設定する。焦点間隔Δｚの具体的な値の決定方法は後述する。装置制御部５１は、これらの計測面ｉ２及びｉ３にそれぞれ対物レンズ２１ａの焦点を配置した際に、結像光学系７により検出部９の検出面に対象物Ｓの像を結像させて対象物Ｓからの光を検出部９に検出させる。
【００８１】
計測面ｉ２及びｉ３に対物レンズ２１ａの焦点を配置したそれぞれの場合に対応する検出部９の各画素からの検出信号は、制御部５０（図２）に入力され、位相復元処理部５２１により各画素の位置と検出信号に基づいた光強度とが対応付けられた光強度分布データがそれぞれ生成される。なお、計測面ｉ１〜ｉ３上に配置される、それぞれの対物レンズ２１ａの焦点の位置は、それぞれ計測面ｉ１〜ｉ３上に存在すれば、ＸＹ平面上の位置は特に限定されない。
【００８２】
生成された計測面ｉ２と計測面ｉ３のそれぞれの計測面に対応するそれぞれの光強度分布データは、計測面ｉ１における光強度のｚに関する微分係数を算出するのに用いられる。位相復元処理部５２１は、計測面ｉ２上の点及び計測面ｉ３上の点であって、ＸＹ平面上での座標が同じ位置等の互いに対応する位置にある２点の強度の値の差分を、計測面ｉ２及びｉ３の間の距離である２Δｚで割ることにより、計測面ｉ１における光強度のｚに関する微分係数に相当するｄＩ／ｄｚ＝（Ｉ３−Ｉ２）／２Δｚを算出する。
【００８３】
位相復元処理部５２１は、算出された光強度のｚに関する微分係数分布に対応するデータ（以下、「微分係数分布データ」という場合がある。）を適宜記憶部４４等に記憶させる。微分係数分布データは、座標ｘ、ｙ、ｚの値に対応する光強度のｚに関する微分係数の分布を示すデータであり、ルックアップテーブルの形式で構築されている。なお、所定の座標ｘ、ｙ、ｚの値に対応するｚに関する微分係数の値を取り出すことができれば、微分係数分布データのデータ構造は特に限定されず、他の既存のデータ構造であってもよい。
【００８４】
図５は、対象物Ｓを撮像する際の対物レンズ２１ａの焦点位置Ｆを計測面ｉ１としたときの配置例を示す概念図である。図５では、ステージ８上に対象物Ｓが載置されている。図５では、対物レンズ２１ａの光軸Ｌｏに平行にｚ軸が設定され、対象物Ｓから対物レンズ２１ａへ向かう向きをｚ軸の+の向きとしている。また、焦点位置Ｆから対物レンズ２１ａに入射する光の光束を二点鎖線Ｌ２００で模式的に示した。図５では、対物レンズ２１ａの焦点位置Ｆをｚ＝ｚ０とし、ｚ＝ｚ０の点を含んでｚ軸に垂直な面を計測面ｉ１として示している。
【００８５】
装置制御部５１は、モータ等の移動装置を介して電動駆動することによりステージ８を移動させ、計測面ｉ１，ｉ２，ｉ３上にそれぞれ対物レンズ２１ａの焦点位置を設定する。装置制御部５１は、照明光学系１０及び対物レンズ２１ａの光軸Ｌ１、Ｌｏに沿って互いに焦点間隔Δｚ離れた複数の位置のそれぞれに対物レンズ２１ａの焦点を順次配置する。ここで、「対物レンズ２１ａの光軸Ｌｏに沿って互いに焦点間隔Δｚ離れた複数の位置」とは、当該複数の位置を対物レンズ２１ａの光軸Ｌｏに射影した際に、互いにΔｚの距離離れる複数の位置を指す。また、対物レンズ２１ａは結像光学系７においてステージ８側にあるため、対物レンズ２１ａの焦点位置Ｆは対物光学系２０又は結像光学系７の焦点位置といいかえることもできる。
【００８６】
装置制御部５１の最適条件計算部５１１は、本実施形態における顕微鏡１００の設定情報と、後に詳述するパラメータｋとを用い、焦点間隔Δｚを以下の式（１００）により算出する。ここで、顕微鏡１００の設定情報とは、定量位相画像の生成のために顕微鏡１００に設定される情報であり、例えば、対物レンズ２１ａの開口数、照明光Ｌ２の波長、対物レンズ２１ａと対象物Ｓとの間の屈折率の情報である。
【００８７】
【数５】
2020188717.tif 000007
【００８８】
ここで、λは照明光Ｌ２の波長であり、入力部４１へのユーザの入力等に基づいて装置制御部５１が設定する。ｎは対物レンズ２１ａと対象物Ｓとの間の屈折率であり、入力部４１へのユーザの入力等に基づいて装置制御部５１が設定する。装置制御部５１は、対物レンズ２１ａが乾燥対物レンズの場合、対物レンズ２１ａと対象物Ｓとの間の雰囲気は空気であるため、空気の屈折率として例えばｎ＝１．００に設定する。一方、対物レンズ２１ａが液浸対物レンズの場合、装置制御部５１は、対物レンズ２１ａと対象物Ｓとの間に満たされる浸液の屈折率を屈折率ｎとして設定する。ＮＡは対物レンズ２１ａの開口数であり、入力部４１へのユーザの入力等に基づいて装置制御部５１が設定する。
【００８９】
そして、最適条件計算部５１１は、顕微鏡１００の設定情報に基づいて、焦点間隔Δｚを算出する。具体的には、最適条件計算部５１１は、式（１００）に示されたように、対物レンズ２１ａの開口数ＮＡ、照明光Ｌ２の波長λ、及び対物レンズ２１ａと対象物Ｓとの間の屈折率ｎの情報に基づいて焦点間隔Δｚを算出する。
【００９０】
また、パラメータｋは１以上の値であり、式（１００）において、ｋ＝１の場合には、対物レンズ２１ａの性能を最大限に生かす遮断空間周波数まで、強度輸送方程式により位相復元を行うことができる。ｋの値が１より大きくなるにしたがって、強度輸送方程式により位相復元を行うことができる空間周波数は低くなる。以下、パラメータｋを位相復元パラメータという場合がある。
【００９１】
なお、対象物Ｓの構造を反映した位相値と測定した位相値との差を位相復元誤差という。高い空間周波数で位相を復元すると、基本的に位相復元誤差が小さくなる。位相復元誤差が小さい場合には、実際の位相値により近い位相を復元することになる（つまり、より高精度に位相を復元することになる。）ので、位相を復元する際の精度（以下、「位相復元精度」という場合がある。）が高いということができる。
【００９２】
最適条件計算部５１１は、位相復元パラメータｋを、２５以下に設定することが好ましく、１６以下に設定することがより好ましい。位相復元パラメータｋが大きくなると、上記の通り、位相復元可能な空間周波数が低くなる。
【００９３】
ｋ＝２５の場合における焦点間隔Δｚを採用した場合では、対物レンズ２１ａの遮断空間周波数の５分の１の空間周波数まで復元することに相当する性能となる。また、ｋ＝１６の場合における焦点間隔Δｚを採用した場合では、対物レンズ２１ａの遮断空間周波数の４分の１の空間周波数まで復元することに相当する性能となる。可視光での数μｍ〜数百μｍ程度の大きさの細胞等の対象物Ｓ（位相物体）の観察では、位相復元パラメータｋは２５以下の比較的大きい値に設定すると細胞全体を検出する上で好ましく、パラメータｋが１６以下だと、核等の細胞内の比較的大きい構造物を検出しやすくなるため更に好ましい。
【００９４】
なお、上記の例においては、対物レンズ２１ａの遮断空間周波数の５分の１や４分の１の空間周波数まで復元することに相当するようにパラメータｋの値を決定したが、パラメータｋの値はこれらの値に限られることはなく、対象物Ｓの構造を反映した位相の計測に必要な分解能に応じて設定すればよい。
【００９５】
最適条件計算部５１１が、照明光Ｌ２の波長λが３４０ｎｍ、４００ｎｍ、５５０ｎｍ、７００ｎｍ及び１３００ｎｍにおける焦点間隔Δｚを、位相復元パラメータｋを１〜１６の範囲として式（１００）に基づいて設定する場合の数値範囲の例を以下の表１に示す。
【００９６】
【表１】
2020188717.tif 000008
【００９７】
最適条件計算部５１１は、特に可視光を照明光とした位相の計測に際し、各種の対物レンズ２１に対し、位相復元パラメータｋを１〜１６の範囲として焦点間隔Δｚを算出した表１の結果に基づき、焦点間隔Δｚを０．０７μｍ以上９０．３μｍ以下に設定することが好ましい。焦点間隔Δｚを０．０７μｍ以上９０．３μｍ以下に設定すると、開口数ＮＡが０．３〜１．４０等の範囲におけるいずれかの対物レンズ２１を用いて、遮断空間周波数の４分の１から遮断空間周波数の全体までを精度よく位相復元することができる。つまり、多くの種類の対物レンズ２１について、対物レンズ２１の遮断空間周波数の大部分を精度よく（いいかえると、対物レンズ２１の性能を活かして）位相復元することができる。
【００９８】
なお、対物レンズ２１の開口数ＮＡを０．３〜１．４０の範囲としたが、この範囲に限られず、任意の開口数ＮＡを有する対物レンズ２１を用いてもよい。また、照明光Ｌ２の波長を、可視光を中心に３４０ｎｍ、４００ｎｍ、５５０ｎｍ、７００ｎｍ及び１３００ｎｍとしたが、この範囲に限られず、可視域の他の波長を含んでもよいし、紫外域や赤外域の波長を含んでもよい。
【００９９】
細胞の定量位相画像を生成した後、当該定量位相画像に基づいて細胞の数、形態又は形状を測定する。例えば、定量位相画像に基づいて、視野内に存在する細胞の、面積、体積、真円度若しくは真球度若しくは縦横比、高さや幅等の大きさ、細胞の位置等を算出することができる。また、多孔膜上に配置された細胞が複数個存在する場合には、細胞同士の間の間隙の面積を測定することもできる。
【０１００】
実施例において後述するように、細胞の数、形態又は形状を測定する工程は、定量位相画像の輝度情報に基づいて行われてもよい。
【０１０１】
また、定量位相画像を生成する工程において、異なる時間における複数の定量位相画像を生成し、細胞の数、形態又は形状を測定する工程において、異なる時間における複数の定量位相画像に基づいて、細胞の数、形態又は形状の変化（経時変化）を測定することもできる。
【０１０２】
［定量位相画像生成装置］
一実施形態において、本発明は、多孔膜上に配置された細胞の数、形態又は形状を測定するための定量位相画像を生成する定量位相画像生成装置であって、顕微鏡と、装置制御部と、検出器と、画像生成部とを備え、前記顕微鏡は、多孔膜上に配置された細胞が配置されるステージと、前記細胞に照明光を照射する照明光学系と、対物レンズとを有し、前記装置制御部は、前記対物レンズの光軸に沿って互いに間隔Δｚ離れた複数の位置のそれぞれに前記対物レンズの焦点が配置されるように前記焦点の位置を調整するものであり、前記検出器は、前記照明光学系から照射された照明光により、前記細胞から生じる光を検出するものであり、前記画像生成部は、前記検出器が検出した前記光に基づいて、前記複数の位置に対応する光強度分布データを生成し、前記光強度分布データに基づいて、定量位相画像を生成するものであり、前記装置制御部は、前記対物レンズの開口数ＮＡ、前記照明光の波長λ、及び、前記対物レンズと前記細胞との間の屈折率ｎの少なくとも１つに基づいて前記間隔Δｚを設定するものである、定量位相画像生成装置を提供する。
【０１０３】
本実施形態に係る定量位相画像生成装置は、図１等を参照して上述した定量位相画像生成装置と同様の構成を有しており、顕微鏡１００と、装置制御部５１と、検出部（検出器）９と、解析部（画像生成部）５２とを備える。検出器９は、照明光学系（透過照明光学系１０又は落射照明光学系１１０）から照射された照明光により、ステージ上の対象物（細胞、多孔膜等）から生じる光を検出する。
【０１０４】
顕微鏡１００は、多孔膜上に配置された細胞が配置されるステージ８と、前記細胞に照明光を照射する照明光学系（透過照明光学系１０又は落射照明光学系１１０）と、対物レンズ２１を有する。
【０１０５】
本実施形態の定量位相画像生成装置のステージ８には、多孔膜上に配置された細胞が配置される。多孔膜上に配置された細胞としては、細胞培養容器、オーガンチップ等が挙げられる。
【０１０６】
本実施形態の定量位相画像生成装置によれば、従来観察が困難であった、多孔膜上の細胞の数、形態又は形状を明瞭に測定することができる。また、多孔膜上の細胞を非侵襲的に観察することができる。
【実施例】
【０１０７】
次に実施例を示して本実施形態を説明するが、本発明は以下の実施例に限定されるものではない。
【０１０８】
［実験例１］
２４ウェルプレート（コーニング社）の各ウェルに０．６ｍＬの細胞培養液を添加した。続いて、各ウェル内に細胞培養用インサート（コーニング社）を設置した。細胞培養用インサートの底面には、ポリエチレンテレフタレート製の多孔膜（ポアサイズ０．４μｍ、ポア密度４×１０４個／ｍｍ２）が配置されていた。続いて、各細胞培養用インサートに、細胞培養液０．１ｍＬに懸濁した腎臓近位尿細管上皮細胞（ＲＰＴＥＣ細胞）を１×１０４個／ｃｍ２の細胞密度で播種した。
【０１０９】
また、比較のために、別の２４ウェルプレート（コーニング社）の各ウェルに、細胞培養用インサートを使用せずに、ＲＰＴＥＣ細胞を１×１０４個／ｃｍ２の細胞密度で播種した。
【０１１０】
続いて、培養した各細胞を顕微鏡で観察した。観察は２４ウェルプレートごと顕微鏡のステージに載せて行い、同一視野において、明視野画像、定量位相画像及び位相差画像を撮影した。定量位相画像の撮像条件は次の通りであった。焦点位置Ｆは細胞とし、Δｚ＝５μｍ、ＮＡ＝０．３、λ＝５５０ｎｍ、ｎ＝１．０、ｋ＝２．１。
【０１１１】
図８（ａ）〜（ｃ）は、細胞培養用インサートを使用していない、通常の細胞培養用プレートで培養した細胞の代表的な写真である。図８（ａ）に示す画像Ａは、位相差画像である。図８（ｂ）に示す画像Ｂは、定量位相画像である。また、図８（ｃ）は、図８（ａ）及び（ｂ）に示す矢印に沿って、横軸に距離（相対値）、縦軸に輝度（相対値）をそれぞれプロットしたグラフ（ラインプロファイル）である。図８（ｃ）中、「画像Ａ／位相差」は図８（ａ）に示す位相差画像のラインプロファイルであることを示し、「画像Ｂ／定量位相」は図８（ｂ）に示す定量位相画像のラインプロファイルであることを示す。
【０１１２】
その結果、多孔膜上でない、通常の細胞培養用プレート上の細胞は、位相差画像でも定量位相画像でも輝度情報に基づくセグメンテーションが容易であり、細胞の数、形態又は形状を明瞭に観察することができることが明らかとなった。図８（ａ）〜（ｃ）の例では、細胞ａ〜細胞ｅの５個の細胞が観察された。また、位相差画像と異なり、定量位相画像では、輝度情報が細胞の高さ情報を示すことが明らかとなった。
【０１１３】
また、図９（ａ）〜（ｃ）は、多孔膜上で培養した細胞の代表的な写真である。図９（ａ）に示す画像Ｃは、位相差画像である。図９（ｂ）に示す画像Ｄは、定量位相画像である。また、図９（ｃ）は、図９（ａ）及び（ｂ）に示す矢印に沿って、横軸に距離（相対値）、縦軸に輝度（相対値）をそれぞれプロットしたグラフ（ラインプロファイル）である。図９（ｃ）中、「画像Ｃ／位相差」は図９（ａ）に示す位相差画像のラインプロファイルであることを示し、「画像Ｄ／定量位相」は図９（ｂ）に示す定量位相画像のラインプロファイルであることを示す。
【０１１４】
その結果、多孔膜上の細胞の位相差画像では、輝度情報に基づくセグメンテーションが不可能であることが明らかとなった。これに対し、定量位相画像では、多孔膜上の細胞であっても輝度情報に基づくセグメンテーションが容易であり、細胞の数、形態又は形状を明瞭に観察することができることが明らかとなった。図９（ａ）〜（ｃ）の例では、細胞ｓ、細胞ｔ、細胞ｕの３個の細胞が観察された。
【０１１５】
これは、位相差画像では、多孔膜の周囲の輝度が高く、セグメンテーション上のノイズになるのに対して、定量位相画像では、多孔膜は輝度値の低い領域となるため、輝度情報への影響が少ないためであると考えられた。図９（ｂ）における黒い点は多孔膜に対応していると考えられた。また、多孔膜上の細胞においても、定量位相画像では、輝度情報を細胞の高さ情報とすることができると考えられた。
【０１１６】
［実施例２］
（多孔膜領域の定量位相画像の検討）
図１０に示すように、細胞（位相物体）内に計測面ｉ１を設定して、複数のｚ位置の光強度分布データから復元した定量位相画像と、多孔膜６１１のポアＰ内に計測面ｉ１００を設定して、複数のｚ位置の光強度分布データから復元した定量位相画像とを演算処理し、細胞が存在しない背景領域に計測されるポアを除去した定量位相画像を構築した。
【０１１７】
細胞（位相物体）内への計測面ｉ１の設定は、上述したように、複数のｚ位置の各々で光強度分布データを取得した後、ｚ位置に関する分散ｖｚデータに基づいて決定した。撮影条件は、Δｚ＝５μｍ、ＮＡ＝０．３、λ＝５５０ｎｍ、ｎ＝１．０、ｋ＝２．１とした。計測面ｉ１を決定した後、計測面ｉ１の上下の異なるｚ位置の光強度分布データ複数枚を再度撮像する必要はなく、分散ｖｚデータを生成するために撮像した複数枚のｚ位置の光強度分布データの中から、定量位相画像を構築する用途に必要な枚数だけを選択して、転用することができた。細胞内に設定した計測面ｉ１のｚ位置の上下の異なるｚ位置の光強度分布データ複数枚から、定量位相画像を生成した。
【０１１８】
続いて、多孔膜６１１のポアＰ内に計測面ｉ１００を設定した定量位相画像を準備した。ポアＰ内の計測面ｉ１００は、上述したようにｚ位置に関する分散ｖｚデータに基づいて決定した。計測面ｉ１００を決定した後、計測面ｉ１００の上下の異なるｚ位置の光強度分布データ複数枚を、再度撮像する必要はなく、分散ｖｚデータを生成するために撮像した複数枚のｚ位置の光強度分布データの中から、定量位相画像を構築する用途に必要な枚数だけを選択して、転用することができた。細胞内に設定した計測面ｉ１００のｚ位置の上下の異なるｚ位置の光強度分布データ複数枚から、定量位相画像を生成した。
【０１１９】
細胞領域に注目すると、計測面ｉ１を含む光強度分布データから復元した定量位相画像は、多孔膜６１１のポアＰ内に計測面ｉ１００を含む光強度分布データから復元した定量位相画像と比較して、細胞構造の細かい凹凸（高周波成分）がより精度良く復元されていた。
【０１２０】
一方、細胞が存在しない背景領域に計測されたポア部に注目すると、ポアＰ内に計測面ｉ１００を含む光強度分布データから復元した定量位相画像の方が、ポアＰの高周波成分がより良好に復元されて、ポア形状がより明瞭であった。
【０１２１】
また、計測面ｉ１を含む光強度分布データから復元した定量位相画像は、ポア形状が明瞭でないため、円形度等の形態に基づくフィルタリングを使用して、細胞ＣとポアＰを分離するセグメンテーションが困難であった。
【０１２２】
また、多孔膜６１１のポアＰ内に計測面ｉ１００を含む光強度分布データから復元した定量位相画像は、背景領域のポアＰがより明瞭であるため、高周波成分を抽出する周波数フィルタリングと、円形度等の形態に基づくフィルタリングを実施することで、ポア部だけの定量位相画像を生成することができた。また、計測面ｉ１を含む光強度分布データから復元した定量位相画像は、ポア部だけの定量位相画像を減算することにより、ポアＰを除去した定量位相画像を生成することができた。
【符号の説明】
【０１２３】
７…結像光学系、８…ステージ、９…検出部（検出器）、１０…透過照明光学系、１１…光源、１２，１５…レンズ、１３…バンドパスフィルター、１４…視野絞り、１６…開口絞り、１７…集光レンズ、２０…対物光学系、２１，２１ａ，２１ｂ，２１ｃ…対物レンズ、３０…リレー光学系、３１…結像レンズ、３２…ビームスプリッター、３３ａ，３３ｂ，３３ｃ…ミラー、３４ａ，３４ｂ，３４ｃ…レンズ、３５…接眼レンズ、４０…情報処理部、４１…入力部、４２…表示部、４３…通信部、４４…記憶部、５０…制御部、５１…装置制御部、５２…解析部（画像生成部）、１００…顕微鏡１００、１１０…落射照明光学系、１１１…水銀ランプ（光源）、１２０…フィルターキューブ、１２１…励起フィルター、１２２…ダイクロイックミラー、１２３…吸収フィルター、５１１…最適条件計算部、５２１…位相復元処理部、５２２…画像構築部、５２３…画像解析部、６００…細胞培養容器、６０１，７１０，７２０…基板、６０２…側壁、６１０…細胞培養用インサート、６１１…多孔膜、６２０，６２１…培地、７００…オーガンチップ、７１１，７２１…面、７４０，７５０…流路、９００…座標軸、Ｌ１（Ｌｏ）…光軸、Ｌ２…照明光、Ｐ１，Ｐ２…位置、Ｐ…ポア、Ｅ…眼、ｉ１，ｉ２，ｉ３，ｉ１００，ｉ１０１，ｉ１０２…計測面、Ｌ２００…光束、Ｗ１，Ｗ２，Ｗ…３波面、Ｃ…細胞、Ｆ…焦点位置、Ｓ…対象物。
【図１】
2020188717.tif 000010
【図２】
2020188717.tif 000011
【図３】
2020188717.tif 000012
【図４】
2020188717.tif 000013
【図５】
2020188717.tif 000014
【図６】
2020188717.tif 000015
【図７】
2020188717.tif 000016
【図８】
2020188717.tif 000017
【図９】
2020188717.tif 000018
【図１０】
2020188717.tif 000019