ホーム > 特許ランキング > 浜松ホトニクス株式会社
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-09-11
(45)【発行日】2024-09-20
(54)【発明の名称】観察装置および観察方法
(51)【国際特許分類】
G02B 21/36 20060101AFI20240912BHJP
G01N 21/17 20060101ALI20240912BHJP
G01N 21/45 20060101ALI20240912BHJP
G02B 21/06 20060101ALI20240912BHJP
【ＦＩ】
G02B21/36
G01N21/17 610
G01N21/45 A
G02B21/06
【請求項の数】
18
(21)【出願番号】P 2023023358
(22)【出願日】2023-02-17
(65)【公開番号】P2024117323
(43)【公開日】2024-08-29
【審査請求日】2024-06-13
【早期審査対象出願】
(73)【特許権者】
【識別番号】000236436
【氏名又は名称】浜松ホトニクス株式会社
(74)【代理人】
【識別番号】100088155
【弁理士】
【氏名又は名称】長谷川 芳樹
(74)【代理人】
【識別番号】100113435
【弁理士】
【氏名又は名称】黒木 義樹
(74)【代理人】
【識別番号】100140442
【弁理士】
【氏名又は名称】柴山 健一
(74)【代理人】
【識別番号】100110582
【弁理士】
【氏名又は名称】柴田 昌聰
(72)【発明者】
【氏名】安彦 修
(72)【発明者】
【氏名】竹内 康造
【審査官】瀬戸 息吹
(56)【参考文献】
【文献】国際公開第２０２２／０５４３０５（ＷＯ，Ａ１）
【文献】特開２０１９－１４４１２１（ＪＰ，Ａ）
【文献】国際公開第２０２１／２６１１５０（ＷＯ，Ａ１）
【文献】特開２０１９－１２４６７９（ＪＰ，Ａ）
【文献】特開２０１５－１５２６６０（ＪＰ，Ａ）
【文献】国際公開第２０２０／０１３３２５（ＷＯ，Ａ１）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｇ０２Ｂ １９／００ ー ２１／３６
Ｇ０１Ｎ ２１／００ ー ２１／０１
Ｇ０１Ｎ ２１／１７ ー ２１／６１
(57)【特許請求の範囲】
【請求項1】
複数の光照射方向それぞれに沿って光が照射された観察対象物の干渉強度画像を取得する干渉強度画像取得部と、複数の光照射方向それぞれについて、前記干渉強度画像に基づいて、平面波の光が第１端面および第２端面のうちの何れか一方に入射して他方まで到達したときの光の第１複素振幅画像を生成する第１複素振幅画像生成部と、
複数の光照射方向それぞれの前記第１複素振幅画像に基づいて、前記第１端面から前記第１端面と前記第２端面との間の第１中間面までの区間の第１屈折率分布画像を生成する屈折率分布画像生成部と、
前記第１屈折率分布画像に基づいて、複数の光照射方向それぞれについて、平面波の光が前記第１端面に入射して前記第１中間面まで到達したときの光の第２複素振幅画像を生成する第２複素振幅画像生成部と、
複数の光照射方向それぞれの前記第１複素振幅画像に基づいてトランスミッション行列Ｔを生成するとともに、複数の光照射方向それぞれの前記第２複素振幅画像に基づいてトランスミッション行列Ｔ１を生成するトランスミッション行列生成部と、
前記トランスミッション行列Ｔと前記トランスミッション行列Ｔ１の逆行列との積に基づいて、複数の光照射方向それぞれについて、平面波の光が前記第１中間面に入射して前記第２端面に到達したときの光の第３複素振幅画像を生成する第３複素振幅画像生成部と、
を備える観察装置。
【請求項2】
前記干渉強度画像取得部は、複数の光照射方向それぞれに沿って前記観察対象物に照射されて前記観察対象物を経た光と参照光との干渉による干渉強度画像を撮像した撮像部から、前記複数の光照射方向それぞれの前記干渉強度画像を取得する、請求項１記載の観察装置。
【請求項3】
前記干渉強度画像取得部は、複数の光照射方向それぞれに沿って前記観察対象物に照射されて前記観察対象物を経た光と、一定の照射方向に沿って前記観察対象物に照射されて前記観察対象物を経た光と、の干渉による干渉強度画像を取得する、請求項１に記載の観察装置。
【請求項4】
前記屈折率分布画像生成部は、複数の光照射方向それぞれの前記第３複素振幅画像に基づいて、前記第１中間面から前記第２端面での区間の屈折率分布画像を生成する、請求項１～３いずれか一項に記載の観察装置。
【請求項5】
前記屈折率分布画像生成部は、複数の光照射方向それぞれの前記第３複素振幅画像に基づいて、前記第１中間面から前記第１中間面と前記第２端面との間の第２中間面までの区間の第２屈折率分布画像を生成し、前記第２複素振幅画像生成部は、前記第２屈折率分布画像に基づいて、複数の光照射方向それぞれについて、平面波の光が前記第１中間面に入射して前記第２中間面まで到達したときの光の第４複素振幅画像を生成し、
前記トランスミッション行列生成部は、複数の光照射方向それぞれの前記第３複素振幅画像に基づいてトランスミッション行列Ｔ３２を生成するとともに、複数の光照射方向それぞれの前記第４複素振幅画像に基づいてトランスミッション行列Ｔ２を生成し、
前記第３複素振幅画像生成部は、前記トランスミッション行列Ｔ３２と前記トランスミッション行列Ｔ２の逆行列との積に基づいて、複数の光照射方向それぞれについて、平面波の光が前記第２中間面に入射して前記第２端面に到達したときの光の第５複素振幅画像を生成する、
請求項１～３いずれか一項に記載の観察装置。
【請求項6】
前記屈折率分布画像生成部は、複数の光照射方向それぞれの前記第３複素振幅画像に基づいて、前記第１中間面から前記第１中間面と前記第２端面との間の第２中間面までの区間の第２屈折率分布画像を生成し、前記第２複素振幅画像生成部は、前記第２屈折率分布画像に基づいて、複数の光照射方向それぞれについて、前記第２複素振幅画像で表される光が前記第１中間面に入射して前記第２中間面まで到達したときの光の第６複素振幅画像を生成し、
前記トランスミッション行列生成部は、複数の光照射方向それぞれの前記第６複素振幅画像に基づいてトランスミッション行列Ｔ２１を生成し、
前記第３複素振幅画像生成部は、前記トランスミッション行列Ｔと前記トランスミッション行列Ｔ２１の逆行列との積に基づいて、複数の光照射方向それぞれについて、平面波の光が前記第２中間面に入射して前記第２端面に到達したときの光の第７複素振幅画像を生成する、
請求項１～３いずれか一項に記載の観察装置。
【請求項7】
前記干渉強度画像取得部は、第１期間および第２期間それぞれにおいて複数の光照射方向それぞれの前記干渉強度画像を取得し、前記第１複素振幅画像生成部は、前記第１期間および前記第２期間それぞれについて、複数の光照射方向それぞれの前記干渉強度画像に基づいて、複数の光照射方向それぞれの前記第１複素振幅画像を生成し、
前記屈折率分布画像生成部は、前記第１期間の複数の光照射方向それぞれの前記第１複素振幅画像に基づいて、前記第１屈折率分布画像を生成し、
前記第２複素振幅画像生成部は、前記第１屈折率分布画像に基づいて、複数の光照射方向それぞれの前記第２複素振幅画像を生成し、
前記トランスミッション行列生成部は、前記第２期間の複数の光照射方向それぞれの前記第１複素振幅画像に基づいてトランスミッション行列Ｔを生成するとともに、前記第１期間の複数の光照射方向それぞれの前記第２複素振幅画像に基づいてトランスミッション行列Ｔ１を生成し、
前記第３複素振幅画像生成部は、前記第２期間のトランスミッション行列Ｔと前記第１期間のトランスミッション行列Ｔ１の逆行列との積に基づいて、複数の光照射方向それぞれの前記第３複素振幅画像を生成する、
請求項１～３いずれか一項に記載の観察装置。
【請求項8】
前記干渉強度画像取得部は、波数空間において光照射方向を表す波数ベクトルの位置が離散的かつ周期的に分布する前記複数の光照射方向それぞれに沿って光が観察対象物に照射されたときに前記干渉強度画像を撮像した撮像部から、前記複数の光照射方向それぞれの前記干渉強度画像を取得し、前記第３複素振幅画像生成部は、前記複数の光照射方向それぞれの波数ベクトルの前記波数空間における位置の周期的分布に基づいて区分される複数の領域それぞれにおいて前記トランスミッション行列Ｔと前記トランスミッション行列Ｔ１の逆行列との積を求め、これら複数の領域それぞれの前記積に基づいて複数の光照射方向それぞれの前記第３複素振幅画像を生成する、
請求項２に記載の観察装置。
【請求項9】
複数の光照射方向それぞれに沿って光が照射された観察対象物の干渉強度画像を取得する干渉強度画像取得ステップと、複数の光照射方向それぞれについて、前記干渉強度画像に基づいて、平面波の光が第１端面および第２端面のうちの何れか一方に入射して他方まで到達したときの光の第１複素振幅画像を生成する第１複素振幅画像生成ステップと、
複数の光照射方向それぞれの前記第１複素振幅画像に基づいて、前記第１端面から前記第１端面と前記第２端面との間の第１中間面までの区間の第１屈折率分布画像を生成する屈折率分布画像生成ステップと、
前記第１屈折率分布画像に基づいて、複数の光照射方向それぞれについて、平面波の光が前記第１端面に入射して前記第１中間面まで到達したときの光の第２複素振幅画像を生成する第２複素振幅画像生成ステップと、
複数の光照射方向それぞれの前記第１複素振幅画像に基づいてトランスミッション行列Ｔを生成するとともに、複数の光照射方向それぞれの前記第２複素振幅画像に基づいてトランスミッション行列Ｔ１を生成するトランスミッション行列生成ステップと、
前記トランスミッション行列Ｔと前記トランスミッション行列Ｔ１の逆行列との積に基づいて、複数の光照射方向それぞれについて、平面波の光が前記第１中間面に入射して前記第２端面に到達したときの光の第３複素振幅画像を生成する第３複素振幅画像生成ステップと、
を備える観察方法。
【請求項10】
前記干渉強度画像取得ステップにおいて、複数の光照射方向それぞれに沿って観察対象物に照射されて前記観察対象物を経た光と参照光との干渉による干渉強度画像を撮像した撮像部から、前記複数の光照射方向それぞれの前記干渉強度画像を取得する、請求項９記載の観察方法。
【請求項11】
前記干渉強度画像取得ステップにおいて、複数の光照射方向それぞれに沿って前記観察対象物に照射されて前記観察対象物を経た光と、一定の照射方向に沿って前記観察対象物に照射されて前記観察対象物を経た光と、の干渉による干渉強度画像を取得する、請求項９に記載の観察方法。
【請求項12】
前記屈折率分布画像生成ステップにおいて、複数の光照射方向それぞれの前記第３複素振幅画像に基づいて、前記第１中間面から前記第２端面での区間の屈折率分布画像を生成する、請求項９に記載の観察方法。
【請求項13】
前記屈折率分布画像生成ステップにおいて、複数の光照射方向それぞれの前記第３複素振幅画像に基づいて、前記第１中間面から前記第１中間面と前記第２端面との間の第２中間面までの区間の第２屈折率分布画像を生成し、前記第２複素振幅画像生成ステップにおいて、前記第２屈折率分布画像に基づいて、複数の光照射方向それぞれについて、平面波の光が前記第１中間面に入射して前記第２中間面まで到達したときの光の第４複素振幅画像を生成し、
前記トランスミッション行列生成ステップにおいて、複数の光照射方向それぞれの前記第３複素振幅画像に基づいてトランスミッション行列Ｔ３２を生成するとともに、複数の光照射方向それぞれの前記第４複素振幅画像に基づいてトランスミッション行列Ｔ２を生成し、
前記第３複素振幅画像生成ステップにおいて、前記トランスミッション行列Ｔ３２と前記トランスミッション行列Ｔ２の逆行列との積に基づいて、複数の光照射方向それぞれについて、平面波の光が前記第２中間面に入射して前記第２端面に到達したときの光の第５複素振幅画像を生成する、
請求項９に記載の観察方法。
【請求項14】
前記屈折率分布画像生成ステップにおいて、複数の光照射方向それぞれの前記第３複素振幅画像に基づいて、前記第１中間面から前記第１中間面と前記第２端面との間の第２中間面までの区間の第２屈折率分布画像を生成し、前記第２複素振幅画像生成ステップにおいて、前記第２屈折率分布画像に基づいて、複数の光照射方向それぞれについて、前記第２複素振幅画像で表される光が前記第１中間面に入射して前記第２中間面まで到達したときの光の第６複素振幅画像を生成し、
前記トランスミッション行列生成ステップにおいて、複数の光照射方向それぞれの前記第６複素振幅画像に基づいてトランスミッション行列Ｔ２１を生成し、
前記第３複素振幅画像生成ステップにおいて、前記トランスミッション行列Ｔと前記トランスミッション行列Ｔ２１の逆行列との積に基づいて、複数の光照射方向それぞれについて、平面波の光が前記第２中間面に入射して前記第２端面に到達したときの光の第７複素振幅画像を生成する、
請求項９に記載の観察方法。
【請求項15】
前記干渉強度画像取得ステップにおいて、第１期間および第２期間それぞれにおいて複数の光照射方向それぞれの前記干渉強度画像を取得し、前記第１複素振幅画像生成ステップにおいて、前記第１期間および前記第２期間それぞれについて、複数の光照射方向それぞれの前記干渉強度画像に基づいて、複数の光照射方向それぞれの前記第１複素振幅画像を生成し、
前記屈折率分布画像生成ステップにおいて、前記第１期間の複数の光照射方向それぞれの前記第１複素振幅画像に基づいて、前記第１屈折率分布画像を生成し、
前記第２複素振幅画像生成ステップにおいて、前記第１屈折率分布画像に基づいて、複数の光照射方向それぞれの前記第２複素振幅画像を生成し、
前記トランスミッション行列生成ステップにおいて、前記第２期間の複数の光照射方向それぞれの前記第１複素振幅画像に基づいてトランスミッション行列Ｔを生成するとともに、前記第１期間の複数の光照射方向それぞれの前記第２複素振幅画像に基づいてトランスミッション行列Ｔ１を生成し、
前記第３複素振幅画像生成ステップにおいて、前記第２期間のトランスミッション行列Ｔと前記第１期間のトランスミッション行列Ｔ１の逆行列との積に基づいて、複数の光照射方向それぞれの前記第３複素振幅画像を生成する、
請求項９に記載の観察方法。
【請求項16】
前記干渉強度画像取得ステップにおいて、波数空間において光照射方向を表す波数ベクトルの位置が離散的かつ周期的に分布する前記複数の光照射方向それぞれに沿って光が観察対象物に照射されたときに前記干渉強度画像を撮像した撮像部から、前記複数の光照射方向それぞれの前記干渉強度画像を取得し、前記第３複素振幅画像生成ステップにおいて、前記複数の光照射方向それぞれの波数ベクトルの前記波数空間における位置の周期的分布に基づいて区分される複数の領域それぞれにおいて前記トランスミッション行列Ｔと前記トランスミッション行列Ｔ１の逆行列との積を求め、これら複数の領域それぞれの前記積に基づいて複数の光照射方向それぞれの前記第３複素振幅画像を生成する、
請求項１０に記載の観察方法。
【請求項17】
請求項９～１６の何れか１項に記載の観察方法の各ステップをコンピュータに実行させるためのプログラム。
【請求項18】
請求項１７に記載のプログラムを記録したコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、観察装置および観察方法に関するものである。
【背景技術】
【0002】
近年、スフェロイドやオルガノイドと呼ばれる３次元の細胞組織を作製する技術が進歩している。また、これらの３次元細胞組織を創薬や再生医療などに応用する研究が進んでいる。これらの３次元細胞組織は光学的に透明な多重散乱体である。このような光学的に透明な散乱体をイメージングする技術として、これまでに多種多様な手法が提案されている。そのうち蛍光プローブを用いるイメージング技術としては、共焦点顕微鏡、多光子顕微鏡、ライトシート顕微鏡が挙げられる。一方、蛍光プローブを用いない非染色・非侵襲のイメージング技術としては、光コヒーレンス・トモグラフィ（Optical Coherence Tomography、ＯＣＴ）などが知られている。
【0003】
スフェロイドやオルガノイドなどのような観察対象物については非染色・非侵襲のイメージングが望まれる場合が多いものの、これらの観察対象物のイメージングにＯＣＴが適用されたという報告例は多くない。その理由としては、ＯＣＴによるイメージングの分解能が低いこと、および、ＯＣＴによるイメージングにより得られた信号の解釈が難しいこと、が考えられる。したがって、現時点では、ゴールドスタンダードとなりうる非染色の３次元細胞組織のイメージング技術は確立されていないと言ってよい。
【0004】
観察対象物の光路長を非染色・非侵襲でイメージングすることができる技術として、定量位相イメージング（Quantitative Phase Imaging、ＱＰＩ）も知られている。ＱＰＩは、観察対象物（例えば細胞）の光路長という物理的な情報を取得することができることから、生物分野で応用が進んでいる。ＱＰＩにより取得した画像を用いて、微分干渉画像や位相差顕微鏡画像などの他の種類の画像を生成することができる。ＱＰＩは、情報量が比較的多い画像を取得することができる技術であり、従来の明視野画像を用いた解析よりハイコンテントな解析にも適用することができると期待されている。また、近年の機械学習による画像認識精度の向上により非染色のイメージング技術を使ったハイコンテントな解析が盛んに研究されており、今後、多重散乱体の非染色イメージングは重要な役割を担うことが期待される。しかし、ＱＰＩは、取得される画像があくまで光路長の２次元への投影であるので、真の３次元の構造を把握できない。
【0005】
また、観察対象物の光路長を非染色・非侵襲でイメージングすることができる技術として、特許文献１に記載されている光回折トモグラフィ（Optical Diffraction Tomography、ＯＤＴ）も知られている。ＯＤＴは、ＱＰＩを３次元イメージング可能な技術に発展させたものであり、観察対象物の３次元屈折率トモグラフィを実現することができる。ＯＤＴを用いて細胞観察を行うことにより、細胞核やミトコンドリアなどの細胞小器官の同定が可能になり、また、３次元的な形態変化の追跡が可能になって、ＱＰＩより更にハイコンテントな解析ができることが期待されている。
【先行技術文献】
【特許文献】
【0006】
【文献】特開２０１７－２１９８２６号公報
【非特許文献】
【0007】
【文献】Pritam Pai, et al, "Opticaltransmission matrix measurement sampled on a dense hexagonal lattice," OSAContinuum, Vol.3, No.3, pp.637-648 (2020).
【文献】Duygu Akbulut, et al,"Optical transmission matrix as a probe of the photonic strength,"PHYSICAL REVIEW A 94, 043817 (2016).
【文献】Elbert G. van Putten, at al."The information age in optics: Measuring the transmission matrix", ［online］，［令和４年１２月２１日検索］，インターネット＜URL：https://physics.aps.org/articles/v3/22?referer=apshome＞
【発明の概要】
【発明が解決しようとする課題】
【0008】
しかしながら、従来のＯＤＴは、数個からなる細胞の観察に適用され得るものの、上記のような３次元細胞組織などの多重散乱体の観察には適用が困難である。何故なら、従来のＯＤＴでは、観察対象物で生じる多重散乱光が多い場合に、取得される画像に多重散乱光の影響が大きく現れるからである。
【0009】
光の散乱とは、光が対象物と相互作用することによって光の進行方向が変えられる現象をいう。特に対象物の屈折率の空間的な不均一さが増大すると、光は対象物を通過する間に対象物と多数回相互作用するようになる。このように対象物と多数回相互作用した光は多重散乱光と呼ばれる。これに対して、対象物と一回のみ相互作用した光は単一散乱光と呼ばれる。多重散乱光は、スペックルの増大および単一散乱-多重散乱比（Single-scattering to Multi-scatteringRatio、ＳＭＲ）の悪化の原因となり、測定の障壁となることが知られている。
【0010】
スペックルは、光が時間的かつ空間的にコヒーレントである場合に、多重散乱光の干渉によって空間的に強度または位相の大きな変化が引き起こされることにより生じる。スペックル発生を抑制するには、時間的または空間的にインコヒーレントである光を出力する光源を用いればよい。例えば、位相差顕微鏡などの通常の明視野顕微鏡は、ハロゲンランプや発光ダイオードなどの空間的かつ時間的にインコヒーレントな光源を用いることで、スペックルのない画像を取得している。
【0011】
ＳＭＲの悪化は、単一散乱光よりも多重散乱光が支配的になって、単一散乱光が多重散乱光の中に埋もれてしまうことにより生じる。観察対象物が大きくなって観察深度が深いほど、単一散乱光の成分は指数関数的に減少する一方で、これと対照的に多重散乱光の成分が増大する。単一散乱光は、その散乱方向が対象物の構造と直接的な対応関係を持っていることから、対象物の構造の測定に用いやすい。一方で、多重散乱光は、対象物の構造との関係が複雑であり、対象物の構造の情報を抽出することが難しい。それ故、単一散乱光を利用したイメージング技術では、多重散乱光の中に単一散乱光が埋もれると（すなわち、ＳＭＲが悪化すると）測定が失敗することが知られている。
【0012】
ＳＭＲ悪化の抑制は、単一散乱光および多重散乱光のうち単一散乱光を選択的に検出するゲーティングと呼ばれる技術により可能である。ゲーティングにより多重散乱光が抑制されるので、ＳＭＲ悪化の抑制と同時にスペックルの抑制も可能である。ゲーティングは、空間、時間および偏光などの自由度を用いて実現される。共焦点顕微鏡は、空間的ゲーティングの一例である。ＯＣＴは、時間的および空間的なゲーティングの一例である。
【0013】
従来のＯＤＴは、多重散乱光の影響を除去していないことから、観察対象物で生じる多重散乱光が多い場合に、取得される画像においてスペックルが増大し、また、ＳＭＲが悪化する。それ故、従来のＯＤＴは、多重散乱光の発生が少ない数個からなる細胞の観察に適用され得るものの、多重散乱光の発生が多い３次元細胞組織などの多重散乱体の観察には適用が困難である。
【0014】
本発明は、上記問題点を解消する為になされたものであり、観察対象物が多重散乱体である場合であっても、多重散乱光の影響を低減して観察対象物を観察することができる観察装置および観察方法を提供することを目的とする。
【課題を解決するための手段】
【0015】
本発明の観察装置の第１態様は、(1) 複数の光照射方向それぞれに沿って光が照射された観察対象物の干渉強度画像を取得する干渉強度画像取得部と、(2) 複数の光照射方向それぞれについて、干渉強度画像に基づいて、平面波の光が第１端面および第２端面のうちの何れか一方に入射して他方まで到達したときの光の第１複素振幅画像を生成する第１複素振幅画像生成部と、(3) 複数の光照射方向それぞれの第１複素振幅画像に基づいて、第１端面から第１端面と第２端面との間の第１中間面までの区間の第１屈折率分布画像を生成する屈折率分布画像生成部と、(4) 第１屈折率分布画像に基づいて、複数の光照射方向それぞれについて、平面波の光が第１端面に入射して第１中間面まで到達したときの光の第２複素振幅画像を生成する第２複素振幅画像生成部と、(5) 複数の光照射方向それぞれの第１複素振幅画像に基づいてトランスミッション行列Ｔを生成するとともに、複数の光照射方向それぞれの第２複素振幅画像に基づいてトランスミッション行列Ｔ１を生成するトランスミッション行列生成部と、(6) トランスミッション行列Ｔとトランスミッション行列Ｔ１の逆行列との積に基づいて、複数の光照射方向それぞれについて、平面波の光が第１中間面に入射して第２端面に到達したときの光の第３複素振幅画像を生成する第３複素振幅画像生成部と、を備える。
【0016】
本発明の観察装置は次のような態様としてもよい。
第２態様では、第１態様に加えて、干渉強度画像取得部は、複数の光照射方向それぞれに沿って光が照射された観察対象物の干渉強度画像を取得する。
第２態様では、第１態様に加えて、干渉強度画像取得部は、複数の光照射方向それぞれに沿って光が照射された観察対象物の干渉強度画像を取得する。
【0017】
第３態様では、第１態様に加えて、干渉強度画像取得部は、複数の光照射方向それぞれに沿って観察対象物に照射されて観察対象物を経た光と、一定の照射方向に沿って観察対象物に照射されて観察対象物を経た光と、の干渉による干渉強度画像を取得する。
【0018】
第４態様では、第１～第３の態様の何れかに加えて、屈折率分布画像生成部は、複数の光照射方向それぞれの第３複素振幅画像に基づいて、第１中間面から第２端面での区間の屈折率分布画像を生成する。
【0019】
第５態様では、第１～第４の態様の何れかに加えて、次のような構成とする。屈折率分布画像生成部は、複数の光照射方向それぞれの第３複素振幅画像に基づいて、第１中間面から第１中間面と第２端面との間の第２中間面までの区間の第２屈折率分布画像を生成する。第２複素振幅画像生成部は、第２屈折率分布画像に基づいて、複数の光照射方向それぞれについて、平面波の光が第１中間面に入射して第２中間面まで到達したときの光の第４複素振幅画像を生成する。トランスミッション行列生成部は、複数の光照射方向それぞれの第３複素振幅画像に基づいてトランスミッション行列Ｔ３２を生成するとともに、複数の光照射方向それぞれの第４複素振幅画像に基づいてトランスミッション行列Ｔ２を生成する。第３複素振幅画像生成部は、トランスミッション行列Ｔ３２とトランスミッション行列Ｔ２の逆行列との積に基づいて、複数の光照射方向それぞれについて、平面波の光が第２中間面に入射して第２端面に到達したときの光の第５複素振幅画像を生成する。
【0020】
第６態様では、第１～第４の態様の何れかに加えて、次のような構成とする。屈折率分布画像生成部は、複数の光照射方向それぞれの第３複素振幅画像に基づいて、第１中間面から第１中間面と第２端面との間の第２中間面までの区間の第２屈折率分布画像を生成する。第２複素振幅画像生成部は、第２屈折率分布画像に基づいて、複数の光照射方向それぞれについて、第２複素振幅画像で表される光が第１中間面に入射して第２中間面まで到達したときの光の第６複素振幅画像を生成する。トランスミッション行列生成部は、複数の光照射方向それぞれの第６複素振幅画像に基づいてトランスミッション行列Ｔ２１を生成する。第３複素振幅画像生成部は、トランスミッション行列Ｔとトランスミッション行列Ｔ２１の逆行列との積に基づいて、複数の光照射方向それぞれについて、平面波の光が第２中間面に入射して第２端面に到達したときの光の第７複素振幅画像を生成する。
【0021】
第７態様では、第１～第６の態様の何れかに加えて、次のような構成とする。干渉強度画像取得部は、第１期間および第２期間それぞれにおいて複数の光照射方向それぞれの干渉強度画像を取得する。第１複素振幅画像生成部は、第１期間および第２期間それぞれについて、複数の光照射方向それぞれの干渉強度画像に基づいて、複数の光照射方向それぞれの第１複素振幅画像を生成する。屈折率分布画像生成部は、第１期間の複数の光照射方向それぞれの第１複素振幅画像に基づいて、第１屈折率分布画像を生成する。第２複素振幅画像生成部は、第１屈折率分布画像に基づいて、複数の光照射方向それぞれの第２複素振幅画像を生成する。トランスミッション行列生成部は、第２期間の複数の光照射方向それぞれの第１複素振幅画像に基づいてトランスミッション行列Ｔを生成するとともに、第１期間の複数の光照射方向それぞれの第２複素振幅画像に基づいてトランスミッション行列Ｔ１を生成する。第３複素振幅画像生成部は、第２期間のトランスミッション行列Ｔと第１期間のトランスミッション行列Ｔ１の逆行列との積に基づいて、複数の光照射方向それぞれの第３複素振幅画像を生成する。
【0022】
第８態様では、第２態様に加えて、次のような構成とする。干渉強度画像取得部は、波数空間において光照射方向を表す波数ベクトルの位置が離散的かつ周期的に分布する複数の光照射方向それぞれに沿って光が観察対象物に照射されたときに干渉強度画像を撮像した撮像部から、複数の光照射方向それぞれの干渉強度画像を取得する。そして、第３複素振幅画像生成部は、複数の光照射方向それぞれの波数ベクトルの波数空間における位置の周期的分布に基づいて区分される複数の領域それぞれにおいてトランスミッション行列Ｔとトランスミッション行列Ｔ１の逆行列との積を求め、これら複数の領域それぞれの積に基づいて複数の光照射方向それぞれの第３複素振幅画像を生成する。
【0023】
本発明の観察方法の第１態様は、(1) 複数の光照射方向それぞれに沿って光が照射された観察対象物の干渉強度画像を取得する干渉強度画像取得ステップと、(2) 複数の光照射方向それぞれについて、干渉強度画像に基づいて、平面波の光が第１端面および第２端面のうちの何れか一方に入射して他方まで到達したときの光の第１複素振幅画像を生成する第１複素振幅画像生成ステップと、(3) 複数の光照射方向それぞれの第１複素振幅画像に基づいて、第１端面から第１端面と第２端面との間の第１中間面までの区間の第１屈折率分布画像を生成する屈折率分布画像生成ステップと、(4) 第１屈折率分布画像に基づいて、複数の光照射方向それぞれについて、平面波の光が第１端面に入射して第１中間面まで到達したときの光の第２複素振幅画像を生成する第２複素振幅画像生成ステップと、(5) 複数の光照射方向それぞれの第１複素振幅画像に基づいてトランスミッション行列Ｔを生成するとともに、複数の光照射方向それぞれの第２複素振幅画像に基づいてトランスミッション行列Ｔ１を生成するトランスミッション行列生成ステップと、(6) トランスミッション行列Ｔとトランスミッション行列Ｔ１の逆行列との積に基づいて、複数の光照射方向それぞれについて、平面波の光が第１中間面に入射して第２端面に到達したときの光の第３複素振幅画像を生成する第３複素振幅画像生成ステップと、を備える。
【0024】
本発明の観察方法は次のような態様としてもよい。
第２態様では、第１態様に加えて、干渉強度画像取得部は、複数の光照射方向それぞれに沿って光が照射された観察対象物の干渉強度画像を取得する。
第２態様では、第１態様に加えて、干渉強度画像取得部は、複数の光照射方向それぞれに沿って光が照射された観察対象物の干渉強度画像を取得する。
【0025】
第３態様では、第１態様に加えて、干渉強度画像取得ステップにおいて、複数の光照射方向それぞれに沿って観察対象物に照射されて観察対象物を経た光と、一定の照射方向に沿って観察対象物に照射されて観察対象物を経た光と、の干渉による干渉強度画像を取得する。
【0026】
第４態様では、第１～第３の態様の何れかに加えて、屈折率分布画像生成ステップにおいて、複数の光照射方向それぞれの第３複素振幅画像に基づいて、第１中間面から第２端面での区間の屈折率分布画像を生成する。
【0027】
第５態様では、第１～第４の態様の何れかに加えて、次のような構成とする。屈折率分布画像生成ステップにおいて、複数の光照射方向それぞれの第３複素振幅画像に基づいて、第１中間面から第１中間面と第２端面との間の第２中間面までの区間の第２屈折率分布画像を生成する。第２複素振幅画像生成ステップにおいて、第２屈折率分布画像に基づいて、複数の光照射方向それぞれについて、平面波の光が第１中間面に入射して第２中間面まで到達したときの光の第４複素振幅画像を生成する。トランスミッション行列生成ステップにおいて、複数の光照射方向それぞれの第３複素振幅画像に基づいてトランスミッション行列Ｔ３２を生成するとともに、複数の光照射方向それぞれの第４複素振幅画像に基づいてトランスミッション行列Ｔ２を生成する。第３複素振幅画像生成ステップにおいて、トランスミッション行列Ｔ３２とトランスミッション行列Ｔ２の逆行列との積に基づいて、複数の光照射方向それぞれについて、平面波の光が第２中間面に入射して第２端面に到達したときの光の第５複素振幅画像を生成する。
【0028】
第６態様では、第１～第４の態様の何れかに加えて、次のような構成とする。屈折率分布画像生成ステップにおいて、複数の光照射方向それぞれの第３複素振幅画像に基づいて、第１中間面から第１中間面と第２端面との間の第２中間面までの区間の第２屈折率分布画像を生成する。第２複素振幅画像生成ステップにおいて、第２屈折率分布画像に基づいて、複数の光照射方向それぞれについて、第２複素振幅画像で表される光が第１中間面に入射して第２中間面まで到達したときの光の第６複素振幅画像を生成する。トランスミッション行列生成ステップにおいて、複数の光照射方向それぞれの第６複素振幅画像に基づいてトランスミッション行列Ｔ２１を生成する。第３複素振幅画像生成ステップにおいて、トランスミッション行列Ｔとトランスミッション行列Ｔ２１の逆行列との積に基づいて、複数の光照射方向それぞれについて、平面波の光が第２中間面に入射して第２端面に到達したときの光の第７複素振幅画像を生成する。
【0029】
第７態様では、第１～第６の態様の何れかに加えて、次のような構成とする。干渉強度画像取得ステップにおいて、第１期間および第２期間それぞれにおいて複数の光照射方向それぞれの干渉強度画像を取得する。第１複素振幅画像生成ステップにおいて、第１期間および第２期間それぞれについて、複数の光照射方向それぞれの干渉強度画像に基づいて、複数の光照射方向それぞれの第１複素振幅画像を生成する。屈折率分布画像生成ステップにおいて、第１期間の複数の光照射方向それぞれの第１複素振幅画像に基づいて、第１屈折率分布画像を生成する。第２複素振幅画像生成ステップにおいて、第１屈折率分布画像に基づいて、複数の光照射方向それぞれの第２複素振幅画像を生成する。トランスミッション行列生成ステップにおいて、第２期間の複数の光照射方向それぞれの第１複素振幅画像に基づいてトランスミッション行列Ｔを生成するとともに、第１期間の複数の光照射方向それぞれの第２複素振幅画像に基づいてトランスミッション行列Ｔ１を生成する。第３複素振幅画像生成ステップにおいて、第２期間のトランスミッション行列Ｔと第１期間のトランスミッション行列Ｔ１の逆行列との積に基づいて、複数の光照射方向それぞれの第３複素振幅画像を生成する。
【0030】
第８態様では、第２態様に加えて、次のような構成とする。干渉強度画像取得ステップにおいて、波数空間において光照射方向を表す波数ベクトルの位置が離散的かつ周期的に分布する複数の光照射方向それぞれに沿って光が観察対象物に照射されたときに干渉強度画像を撮像した撮像部から、複数の光照射方向それぞれの干渉強度画像を取得する。そして、第３複素振幅画像生成ステップにおいて、複数の光照射方向それぞれの波数ベクトルの波数空間における位置の周期的分布に基づいて区分される複数の領域それぞれにおいてトランスミッション行列Ｔとトランスミッション行列Ｔ１の逆行列との積を求め、これら複数の領域それぞれの積に基づいて複数の光照射方向それぞれの第３複素振幅画像を生成する。
【0031】
本発明のプログラムは、上記の本発明の観察方法の各ステップをコンピュータに実行させるためのものである。本発明の記録媒体は、上記の本発明のプログラムを記録したコンピュータ読み取り可能なものである。
【発明の効果】
【0032】
本発明によれば、観察対象物が多重散乱体である場合であっても、多重散乱光の影響を低減して観察対象物を観察することができる。
【図面の簡単な説明】
【0033】
【図17】図１７は、光照射方向の数が複素振幅画像の画素数と等しい場合の（ａ）ｋｘｋｙ波数空間での光照射方向の分布および（ｂ）或る１点ｒ０に光を集光照射したときのトランスミッション行列Ｔrr（ｒout；ｒin＝ｒ０）のｘｙ空間での分布を示す図である。
【発明を実施するための形態】
【0034】
以下、添付図面を参照して、本発明を実施するための形態を詳細に説明する。なお、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。本発明は、これらの例示に限定されるものではなく、特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
【0035】
図１は、観察装置１Ａの構成を示す図である。この観察装置１Ａは、光源１１、レンズ１２、レンズ２１、ミラー２２、レンズ２３、コンデンサレンズ２４、対物レンズ２５、ビームスプリッタ４１、レンズ４２、撮像部４３および解析部５０などを備える。
【0036】
光源１１は、空間的・時間的にコヒーレントな光を出力するものであり、好適にはレーザ光源である。レンズ１２は、光源１１と光学的に接続されており、光源１１から出力された光を光ファイバ１４の光入射端１３に集光して、その光を光入射端１３に入射させる。光ファイバ１４は、レンズ１２により光入射端１３に入射された光をファイバカプラ１５へ導光する。ファイバカプラ１５は、光ファイバ１４と光ファイバ１６，１７との間で光を結合するものであり、光ファイバ１４により導光されて到達した光を２分岐して、一方の分岐光を光ファイバ１６により導光させ、他方の分岐光を光ファイバ１７により導光させる。光ファイバ１６により導光された光は光出射端１８から発散光として出射される。光ファイバ１７により導光された光は光出射端１９から発散光として出射される。
【0037】
レンズ２１は、光出射端１８と光学的に接続されており、光出射端１８から発散光として出力された光をコリメートする。ミラー２２は、レンズ２１と光学的に接続されており、レンズ２１から到達した光をレンズ２３へ反射させる。ミラー２２の反射面の方位は可変である。レンズ２３は、ミラー２２と光学的に接続されている。コンデンサレンズ２４は、レンズ２３と光学的に接続されている。レンズ２３およびコンデンサレンズ２４は、好適には４ｆ光学系を構成している。レンズ２３およびコンデンサレンズ２４は、ミラー２２の反射面の方位に応じた光照射方向から観察対象物Ｓに対して光を照射する。対物レンズ２５は、コンデンサレンズ２４と光学的に接続されている。対物レンズ２５とコンデンサレンズ２４との間に観察対象物Ｓが配置される。対物レンズ２５は、コンデンサレンズ２４から出力されて観察対象物Ｓを経た光（物体光）を入力し、その光をビームスプリッタ４１へ出力する。
【0038】
ビームスプリッタ４１は、対物レンズ２５と光学的に接続され、また、光出射端１９とも光学的に接続されている。ビームスプリッタ４１は、対物レンズ２５から出力されて到達した光（物体光）と、光出射端１９から出力されて到達した光（参照光）とを合波して、両光をレンズ４２へ出力する。レンズ４２は、ビームスプリッタ４１と光学的に接続されており、ビームスプリッタ４１から到達した物体光および参照光それぞれをコリメートして撮像部４３へ出力する。撮像部４３は、レンズ４２と光学的に接続されており、レンズ４２から到達した物体光と参照光との干渉による干渉縞像（干渉強度画像）を撮像する。撮像部４３の撮像面への物体光の入射方向に対して参照光の入射方向は傾斜している。ビームスプリッタ４１により物体光と参照光とが合波される位置は、結像レンズより後段であってもよいが、収差の影響を考慮すると、図に示されるように対物レンズ２５とレンズ４２との間であるのが望ましい。
【0039】
解析部５０は、撮像部４３と電気的に接続されており、撮像部４３により撮像された干渉強度画像を入力して、この干渉強度画像に基づいて所要の処理をする。解析部５０の詳細については後述する。
【0040】
【0041】
レンズ３１は、光出射端１９と光学的に接続されており、光出射端１９から発散光として出力された光（参照光）をコリメートする。ミラー３２は、レンズ３１と光学的に接続されており、レンズ３１から到達した光をレンズ３４へ反射させる。レンズ３４は、ミラー３２と光学的に接続されており、ミラー３２から到達した光をビームスプリッタ４１へ出力する。レンズ３４から出力された光は、ビームスプリッタ４１の手前で一旦集光された後、発散光としてビームスプリッタ４１に入力される。ビームスプリッタ４１は、対物レンズ２５から出力されて到達した光（物体光）と、レンズ３４から出力されて到達した光（参照光）とを合波して、両光を同軸にしてレンズ４２へ出力する。撮像部４３は、レンズ４２から到達した物体光と参照光との干渉による干渉縞像（干渉強度画像）を撮像する。撮像部４３の撮像面への物体光の入射方向に対して参照光の入射方向は平行である。
【0042】
駆動部３３は、ミラー３２の反射面に垂直な方向にミラー３２を移動させる。駆動部３３は例えばピエゾアクチュエータである。このミラー３２の移動により、ファイバカプラ１５における光分岐からビームスプリッタ４１における合波に至るまでの物体光および参照光それぞれの光路長の差（位相差）を変化させる。この光路長差が異なると、撮像部４３により撮像される干渉強度画像も異なる。
【0043】
観察装置は、図１および図２の構成例に限定されるものではなく、様々な変形が可能である。観察装置１Ａ（図１）および観察装置１Ｂ（図２）の構成では観察対象物Ｓを透過した光を物体光としたが、以下に説明する観察装置１Ｃ（図３）の構成のように観察対象物Ｓで反射された光を物体光としてもよい。
【0044】
図３は、観察装置１Ｃの構成を示す図である。観察装置１Ｃは、光源１１、レンズ１２、レンズ２１、ミラー２２、レンズ２３、対物レンズ２５、ビームスプリッタ４１、レンズ４２、撮像部４３および解析部５０などを備える。以下では、観察装置１Ａ（図１）と相違する点について主に説明する。
【0045】
レンズ２１は、光ファイバ１６の光出射端１８と光学的に接続されており、光出射端１８から発散光として出力された光をコリメートする。ミラー２２は、レンズ２１と光学的に接続されており、レンズ２１から到達した光をレンズ２３へ反射させる。ミラー２２の反射面の方位は可変である。レンズ２３は、ミラー２２と光学的に接続されている。対物レンズ２５は、レンズ２３と光学的に接続されている。レンズ２３と対物レンズ２５との間にビームスプリッタ４１が配置されている。レンズ２３および対物レンズ２５は、好適には４ｆ光学系を構成している。レンズ２３および対物レンズ２５は、ミラー２２の反射面の方位に応じた光照射方向から観察対象物Ｓに対して光を照射する。対物レンズ２５は、観察対象物Ｓで反射された光（物体光）を入力し、その光をビームスプリッタ４１へ出力する。
【0046】
ビームスプリッタ４１は、対物レンズ２５と光学的に接続され、また、光ファイバ１７の光出射端１９とも光学的に接続されている。ビームスプリッタ４１は、対物レンズ２５から出力されて到達した光（物体光）と、光出射端１９から出力されて到達した光（参照光）とを合波して、両光をレンズ４２へ出力する。レンズ４２は、ビームスプリッタ４１と光学的に接続されており、ビームスプリッタ４１から到達した物体光および参照光それぞれをコリメートして撮像部４３へ出力する。撮像部４３は、レンズ４２と光学的に接続されており、レンズ４２から到達した物体光と参照光との干渉による干渉縞像（干渉強度画像）を撮像する。撮像部４３の撮像面への物体光の入射方向に対して参照光の入射方向は傾斜している。ビームスプリッタ４１により物体光と参照光とが合波される位置は、結像レンズより後段であってもよいが、収差の影響を考慮すると、図に示されるように対物レンズ２５とレンズ４２との間であるのが望ましい。
【0047】
観察装置１Ｃ（図３）の構成において、観察装置１Ｂ（図２）と同様に参照光の光路長を変化させる機構（図２中のレンズ３１、ミラー３２、駆動部３３およびレンズ３４）を設けて、ファイバカプラ１５における光分岐からビームスプリッタ４１における合波に至るまでの物体光および参照光それぞれの光路長の差（位相差）を変化させてもよい。この場合、撮像部４３の撮像面への物体光の入射方向に対して参照光の入射方向は平行であってよい。
【0048】
図４は、観察装置１Ａ～１Ｃの解析部５０の構成を示す図である。解析部５０は、干渉強度画像取得部５１、第１複素振幅画像生成部５２、屈折率分布画像生成部５３、第２複素振幅画像生成部５４、トランスミッション行列生成部５５、第３複素振幅画像生成部５６、表示部５７および記憶部５８を備える。解析部５０はコンピュータであってよい。解析部５０は、ＣＰＵ、ＧＰＵ、ＤＳＰまたはＦＰＧＡ等の処理装置、および、ハードディスクドライブ、フラッシュメモリ、ＲＡＭまたはＲＯＭ等の記憶装置を含む。
【0049】
干渉強度画像取得部５１は、ミラー２２の反射面の方位を変化させることにより、観察対象物Ｓに対して複数の光照射方向それぞれに沿って光を照射させる。また、干渉強度画像取得部５１は、複数の光照射方向それぞれについて基準位置における干渉強度画像を撮像部４３から取得する。干渉強度画像取得部５１は、ミラー２２の反射面の方位を変化させる為の制御信号を出力する出力ポートを有し、また、撮像部４３から干渉強度画像を入力する入力ポートを有する。対物レンズ２５を光軸方向に移動させる必要はない。基準位置は、撮像部４３の撮像面に対して共役関係にある像面位置である。
【0050】
第１複素振幅画像生成部５２、屈折率分布画像生成部５３、第２複素振幅画像生成部５４、トランスミッション行列生成部５５および第３複素振幅画像生成部５６は、干渉強度画像に基づいて処理を行う。表示部５７は、処理すべき画像、処理途中の画像および処理後の画像などを表示するものであり、例えば液晶ディスプレイを含む。記憶部５８は、各種の画像のデータを記憶する。第１複素振幅画像生成部５２、屈折率分布画像生成部５３、第２複素振幅画像生成部５４、トランスミッション行列生成部５５、第３複素振幅画像生成部５６、表示部５７および記憶部５８は、クラウドコンピューティングによって構成されてもよい。
【0051】
記憶部５８は、干渉強度画像取得部５１、第１複素振幅画像生成部５２、屈折率分布画像生成部５３、第２複素振幅画像生成部５４、トランスミッション行列生成部５５および第３複素振幅画像生成部５６に各処理を実行させるためのプログラムをも記憶する。このプログラムは、観察装置の製造時または出荷時に記憶部５８に記憶されていてもよいし、出荷後に通信回線を経由して取得されたものが記憶部５８に記憶されてもよいし、コンピュータ読み取り可能な記録媒体２に記録されていたものが記憶部５８に記憶されてもよい。記録媒体２は、フレキシブルディスク、ＣＤ-ＲＯＭ、ＤＶＤ-ＲＯＭ、ＢＤ-ＲＯＭ、ＵＳＢメモリなど任意である。
【0052】
第１複素振幅画像生成部５２、第２複素振幅画像生成部５４および第３複素振幅画像生成部５６は、何れも光の複素振幅画像を生成するものであるが、複素振幅画像の生成の為の処理の内容の点で相違する。干渉強度画像取得部５１、第１複素振幅画像生成部５２、屈折率分布画像生成部５３、第２複素振幅画像生成部５４、トランスミッション行列生成部５５、第３複素振幅画像生成部５６それぞれの処理の詳細については後述する。
【0053】
図５は、観察装置１Ａ～１Ｃを用いて行う観察方法の概要を説明する図である。第１複素振幅画像生成部５２は、干渉強度画像取得部５１により取得された複数の光照射方向それぞれの干渉強度画像に基づいて第１複素振幅画像を生成し、さらに、トランスミッション行列生成部５５は、この第１複素振幅画像に基づいてトランスミッション行例Ｔを生成する。レンズ２５の光軸に垂直な複数の面（図中の例では第１端面、第１中間面、第２中間面、第２端面）により観察対象物を含む領域を複数（図中の例では３個）のブロックに区分して、第１端面と第１中間面との間の第１ブロックについてのトランスミッション行例をＴ１とし、第１中間面と第２中間面との間の第２ブロックについてのトランスミッション行例をＴ２とし、第２中間面と第２端面との間の第３ブロックについてのトランスミッション行例をＴ３とする。全体のトランスミッション行列Ｔは、トランスミッション行例Ｔ１，Ｔ２，Ｔ３の積（Ｔ３Ｔ２Ｔ１）で近似的に表される。このことを利用して、第１～第３のブロックのうち、先ず、第１ブロックの屈折率分布画像を生成して、これを用いて第１中間面における複素振幅画像を生成し、次に、第２ブロックの屈折率分布画像を生成して、これを用いて第２中間面における複素振幅画像を生成する、というような処理を繰り返す。このようにすることで、深達度の高い屈折率分布画像を生成する。
【0054】
図６は、観察方法のフローチャートである。このフローチャートに示される観察方法は、観察装置１Ａ（図１）、観察装置１Ｂ（図２）および観察装置１Ｃ（図３）の何れを用いた場合においても可能なものである。図７は、観察方法における各画像の生成手順を説明する図である。この観察方法は、屈折率分布画像を生成するものであって、干渉強度画像取得ステップＳ５１、第１複素振幅画像生成ステップＳ５２、第１屈折率分布画像生成ステップＳ５３、第２複素振幅画像生成ステップＳ５４、トランスミッション行列生成ステップＳ５５および第３複素振幅画像生成ステップＳ５６を備える。
【0055】
干渉強度画像取得ステップＳ５１の処理は干渉強度画像取得部５１により行われる。第１複素振幅画像生成ステップＳ５２の処理は第１複素振幅画像生成部５２により行われる。第１屈折率分布画像生成ステップＳ５３の処理は屈折率分布画像生成部５３により行われる。第２複素振幅画像生成ステップＳ５４の処理は第２複素振幅画像生成部５４により行われる。トランスミッション行列生成ステップＳ５５の処理はトランスミッション行列生成部５５により行われる。第３複素振幅画像生成ステップＳ５６の処理は第３複素振幅画像生成部５６により行われる。
【0056】
干渉強度画像取得ステップＳ５１において、干渉強度画像取得部５１は、ミラー２２の反射面の方位を変化させることにより、観察対象物Ｓに対して複数の光照射方向それぞれに沿って光を照射させる。そして、干渉強度画像取得部５１は、複数の光照射方向それぞれについて基準位置における干渉強度画像を撮像部４３から取得する。
【0057】
図１，図２および図３それぞれにおいて説明の便宜のためにｘｙｚ直交座標系が示されている。ｚ軸は対物レンズ２５の光軸に対し平行である。基準位置は、撮像部４３の撮像面に対して共役関係にある像面位置である。この位置をｚ＝０とする。観察対象物Ｓへの光照射方向は、その照射光の波数ベクトル（ｋｘ，ｋｙ，ｋｚ）のうちのｋｘ及びｋｙにより表すことができる。
【0058】
図８は、干渉強度画像取得ステップＳ５１における観察対象物Ｓへの光照射方向の走査の例を示す図である。この図では、横軸をｋｘとし縦軸をｋｙとしたｋｘｋｙ平面において各々の丸印の位置が光照射方向を表している。光照射方向の走査は、ｋｘｋｙ平面において、（ａ）矩形格子状に離散的かつ周期的に配置されるようなものであってもよいし、（ｂ）ハニカム状に離散的かつ周期的に配置されるようなものであってもよいし、（ｃ）六方格子状に離散的かつ周期的に配置されるようなものであってもよいし、（ｄ）同心の複数の円それぞれの周上に離散的に配置されるようなものであってもよいし、（ｅ）螺旋状に離散的に配置されるようなものであってもよい。何れの場合にも、図１および図２の構成におけるコンデンサレンズ２４または図３の構成における対物レンズ２５のＮＡの許す限りで光照射方向の走査が可能である。ラスタスキャンおよびランダムスキャンの何れであってもよい。ラスタスキャンの場合には、戻りスキャンが有ってもよいし無くてもよい。
【0059】
第１複素振幅画像生成ステップＳ５２において、第１複素振幅画像生成部５２は、複数の光照射方向それぞれについて、干渉強度画像取得部５１により取得された干渉強度画像に基づいて第１複素振幅画像を生成する。観察装置１Ａ（図１）および観察装置１Ｃ（図３）の場合には、第１複素振幅画像生成部５２は、フーリエ縞解析法により、１枚の干渉強度画像に基づいて第１複素振幅画像を生成することができる。観察装置１Ｂ（図２）の場合には、第１複素振幅画像生成部５２は、位相シフト法により、物体光と参照光との間の光路長差（位相差）が互いに異なる３枚以上の干渉強度画像に基づいて第１複素振幅画像を生成することができる。この第１複素振幅画像は、平面波の光が第１端面および第２端面のうちの何れか一方に入射して他方まで到達したときの光の複素振幅画像である。
【0060】
第１屈折率分布画像生成ステップＳ５３において、屈折率分布画像生成部５３は、第１複素振幅画像生成部５２により生成された複数の光照射方向それぞれの第１複素振幅画像に基づいて、第１端面から第１中間面までの区間の第１屈折率分布画像を生成する。屈折率分布画像を生成する処理の詳細については後述する。
【0061】
第２複素振幅画像生成ステップＳ５４において、第２複素振幅画像生成部５４は、屈折率分布画像生成部５３により生成された第１屈折率分布画像に基づいて、複数の光照射方向それぞれについて、平面波の光が第１端面に入射して第１中間面まで到達したときの光の第２複素振幅画像を生成する。第２複素振幅画像を生成する処理の詳細についても後述する。
【0062】
トランスミッション行列生成ステップＳ５５において、トランスミッション行列生成部５５は、第１複素振幅画像生成部５２により生成された複数の光照射方向それぞれの第１複素振幅画像に基づいてトランスミッション行列Ｔを生成する。また、トランスミッション行列生成部５５は、第２複素振幅画像生成部５４により生成された複数の光照射方向それぞれの第２複素振幅画像に基づいてトランスミッション行列Ｔ１を生成する。トランスミッション行列の詳細についても後述する。
【0063】
第３複素振幅画像生成ステップＳ５６において、第３複素振幅画像生成部５６は、トランスミッション行列生成部５５により生成されたトランスミッション行列Ｔとトランスミッション行列Ｔ１の逆行列との積（ＴＴ１
-1＝Ｔ３Ｔ２）に基づいて、複数の光照射方向それぞれについて、平面波の光が第１中間面に入射して第２端面に到達したときの光の第３複素振幅画像を生成する。
【0064】
第３複素振幅画像生成ステップＳ５６の後、屈折率分布画像生成部５３は、第３複素振幅画像生成部５６により生成された複数の光照射方向それぞれの第３複素振幅画像に基づいて、第１中間面から第２端面までの区間の屈折率分布画像を生成してもよい。或いは、図９～図１２を用いて説明するとおり、屈折率分布画像生成部５３は、第１中間面から第２中間面までの区間の第２屈折率分布画像を生成し続いて第２中間面から第２端面までの区間の第３屈折率分布画像を生成してもよい。このように、一方の端面の側のブロックから順に屈折率分布画像を生成していくことで、深達度の高い屈折率分布画像を生成することができる。
【0065】
【0066】
第２屈折率分布画像生成ステップＳ６３において、屈折率分布画像生成部５３は、第３複素振幅画像生成部５６により生成された複数の光照射方向それぞれの第３複素振幅画像に基づいて、第１中間面から第２中間面までの区間の第２屈折率分布画像を生成する。
【0067】
第４複素振幅画像生成ステップＳ６４において、第２複素振幅画像生成部５４は、第２屈折率分布画像に基づいて、複数の光照射方向それぞれについて、平面波の光が第１中間面に入射して第２中間面まで到達したときの光の第４複素振幅画像を生成する。
【0068】
トランスミッション行列生成ステップＳ６５において、トランスミッション行列生成部５５は、複数の光照射方向それぞれの第３複素振幅画像に基づいてトランスミッション行列Ｔ３２（＝Ｔ３Ｔ２）を生成し、また、複数の光照射方向それぞれの第４複素振幅画像に基づいてトランスミッション行列Ｔ２を生成する。
【0069】
第５複素振幅画像生成ステップＳ６６において、第３複素振幅画像生成部５６は、トランスミッション行列Ｔ３２とトランスミッション行列Ｔ２の逆行列との積（Ｔ３２Ｔ２
-1＝Ｔ３）に基づいて、複数の光照射方向それぞれについて、平面波の光が第２中間面に入射して第２端面に到達したときの光の第５複素振幅画像を生成する。
【0070】
第３屈折率分布画像生成ステップＳ６７において、屈折率分布画像生成部５３は、第５複素振幅画像生成部６６により生成された複数の光照射方向それぞれの第５複素振幅画像に基づいて、第２中間面から第２端面までの区間の第３屈折率分布画像を生成する。
【0071】
【0072】
第２屈折率分布画像生成ステップＳ７３において、屈折率分布画像生成部５３は、第３複素振幅画像生成部５６により生成された複数の光照射方向それぞれの第３複素振幅画像に基づいて、第１中間面から第２中間面までの区間の第２屈折率分布画像を生成する。
【0073】
第６複素振幅画像生成ステップＳ７４において、第２複素振幅画像生成部５４は、第２屈折率分布画像に基づいて、複数の光照射方向それぞれについて、第２複素振幅画像生成部５４により生成された第２複素振幅画像で表される光が第１中間面に入射して第２中間面まで到達したときの光の第６複素振幅画像を生成する。
【0074】
トランスミッション行列生成ステップＳ７５において、トランスミッション行列生成部５５は、複数の光照射方向それぞれの第６複素振幅画像に基づいてトランスミッション行列Ｔ２１（＝Ｔ２Ｔ１）を生成する。
【0075】
第７複素振幅画像生成ステップＳ７６において、第３複素振幅画像生成部５６は、トランスミッション行列Ｔとトランスミッション行列Ｔ２１の逆行列との積（ＴＴ２１
-1＝Ｔ３）に基づいて、複数の光照射方向それぞれについて、平面波の光が第２中間面に入射して第２端面に到達したときの光の第７複素振幅画像を生成する。
【0076】
第３屈折率分布画像生成ステップＳ７７において、屈折率分布画像生成部５３は、第７複素振幅画像生成部７６により生成された複数の光照射方向それぞれの第７複素振幅画像に基づいて、第２中間面から第２端面までの区間の第３屈折率分布画像を生成する。
【0077】
次に、トランスミッション行列の詳細について説明する。図１３は、撮像部により干渉強度画像を撮像するときの入力光Ｕin（ｋin）および出力光ｕout（ｒout）を示す図である。Ｕin（ｋin）は、観察対象物へ平面波として照射される光の波数ｋinの複素振幅を表す。ｕout（ｒout）は、観察対象物から出力される光の位置ｒoutの複素振幅を表す。
【0078】
Ｕin（ｋin）とｕout（ｒout）との間の関係は、下記（１）式で表される。列ベクトルＵinの第ｎ要素Ｕin（ｋin
n）は、波数ｋin
nの平面波の複素振幅を表す。列ベクトルｕoutの第ｎ要素ｕout（ｒout
n）は、位置ｒout
nで観測される光の複素振幅を表す。この式にあるＮ行Ｎ列の行列Ｔrk（ｒout；ｋin）は、Ｕin（ｋin）とｕout（ｒout）との間の線形な関係を表すものであって、トランスミッション行列と呼ばれる。このようなトランスミッション行列により、観察対象物における光の散乱過程を表すことができる。行列Ｔrk（ｒout；ｋin）の第ｎ１行第ｎ２列の要素Ｔn1,n2は、波数ｋin
n2で振幅１の平面波が入力されたときに位置ｒout
n1で観測される光の複素振幅を表す。
【0079】
【数1】
【0080】
複数の光照射方向それぞれについて撮像部により干渉強度画像を撮像するときの第ｎの光照射方向の入力光のベクトルＵin
n（ｋin）は、下記（２）式で表され、第ｎ要素の値のみが１であって、他の要素の値が０である。この入力光Ｕin
n（ｋin）に対して、出力光ｕout
n（ｒout）は、下記（３）式で表される。この（３）式は、第ｎの光照射方向の際に得られた複素振幅に対応する。
【0081】
【数2】
【0082】
【数3】
【0083】
この（２）式および上記（１）式から、下記（４）式が得られる。そして、複数の光照射方向それぞれについて同様に求めると、下記（５）式が得られる。このようにして、トランスミッション行列Ｔrk（ｒout；ｋin）を求めることができる。
【0084】
【数4】
【0085】
【数5】
【0086】
このトランスミッション行列Ｔrk（ｒout；ｋin）は、入力が波数基底であって、出力が位置基底のものである。下記（６）式により、このトランスミッション行列Ｔrk（ｒout；ｋin）を入力に関してフーリエ変換することで、入力および出力の双方が位置基底であるトランスミッション行列Ｔrr（ｒout；ｒin）を生成することができる。
【0087】
【数6】
【0088】
このようにして生成されたトランスミッション行列Ｔrk（ｒout；ｋin），Ｔrr（ｒout；ｒin）は、次のような光学的意味を有するものである（非特許文献１～３）。入力が波数基底であって出力が位置基底であるトランスミッション行列Ｔrk（ｒout；ｋin）は、観察対象物が存在しない場合（入力）の第１面上の各波数ｋinの光の波面に対して、観察対象物が存在する場合（出力）の第２面上の各位置ｒoutの光の波面に対する透過率を与える行列である。入力および出力の双方が位置基底であるトランスミッション行列Ｔrr（ｒout；ｒin）は、観察対象物が存在しない場合（入力）の第１面上の各位置ｒinの光の波面に対して、観察対象物が存在する場合（出力）の第２面上の各位置ｒoutの光の波面に対する透過率を与える行列である。
【0089】
第３複素振幅画像生成部５６は、或るトタンスミッション行例と他の或るトタンスミッション行例の逆行列との積を計算する際に、両者ともに入力が波数基底であって出力が位置基底であるトランスミッション行列Ｔrkとしておくか、或いは、両者ともに入力および出力の双方が位置基底であるトランスミッション行列Ｔrrとしておく。また、例えば、トランスミッション行列Ｔとトランスミッション行列Ｔ１の逆行列との積を計算する際に、光の自由伝搬の数値計算に基づいてトランスミッション行列Ｔ１を第１中間位置におけるトランスミッション行列としておく。他のトランスミッション行列と他のトランスミッション行列の逆行列との積の際も同様である。
【0090】
複素微分干渉画像Ｗ（ｒout）は、トランスミッション行列Ｔrrを用いて、下記（７）式または（８）式により生成され得る。Ｔrr（ｒout；ｒin）に対してｒoutおよびｒinの双方をδｒだけシアーしたＴrr（ｒout－δｒ；ｒin－δｒ）の複素共役をＴrr
*（ｒout－δｒ；ｒin－δｒ）とする。
【0091】
【数7】
【0092】
【数8】
【0093】
（７）式は、第２面上の各位置ｒoutについて、第１面においてＴrr（ｒout；ｒin）とＴrr
*（ｒout－δｒ；ｒin－δｒ）との要素ごとの積の総和を求めることで、複素微分干渉画像Ｗ（ｒout）を生成する。この総和計算は、出力側空間分解の複素微分干渉画像Ｗ（ｒout）を生成するものである。すなわち、入力に関しては空間分解能を有しない一方で、出力に関しては空間分解能を有するとして、入力の基底について総和を求める。これは、観察対象物の全体に同時に光を照射する（すなわち、入力に関しては空間分解能を有しない）こととし、観察対象物の二次元像を取得する（すなわち、出力に関しては空間分解能を有する）系と相似している。この系では、取得された画像では、観察対象物と受光部との間に存在する散乱体の影響が大きく現れる。
【0094】
（８）式は、第１面上の各位置ｒinについて、第２面においてＴrr（ｒout；ｒin）とＴrr
*（ｒout－δｒ；ｒin－δｒ）との要素ごとの積の総和を求めることで、複素微分干渉画像Ｗ（ｒin）を生成する。この総和計算は、入力側空間分解の複素微分干渉画像Ｗ（ｒin）を生成するものである。すなわち、入力に関しては空間分解能を有する一方で、出力に関しては空間分解能を有しないとして、出力の基底について総和を求める。これは、観察対象物上の各位置に光を集光照射するとともに当該集光点を走査する（すなわち、入力に関しては空間分解能を有する）こととし、各集光点について観察対象物を経た光の全光量を測定する（すなわち、出力に関しては空間分解能を有しない）系と相似している。この系では、観察対象物と受光部との間に存在する散乱体の影響が抑制された観察対象物の画像を取得することができる。
【0095】
δｒのｘ成分δｘおよびｙ成分δｙのうち少なくとも一方は非０である。δｘ≠０，δｙ＝０であれば、ｘ方向をシアー方向とする複素微分干渉画像が得られる。δｘ＝０，δｙ≠０であれば、ｙ方向をシアー方向とする複素微分干渉画像が得られる。δｘ≠０，δｙ≠０であれば、δｘとδｙとの比に応じた方向をシアー方向とする複素微分干渉画像が得られる。
【0096】
出力側空間分解の複素微分干渉画像Ｗｎ（ｒout）に基づいて、下記（９）式により、位相微分画像を生成することができる。或いは、入力側空間分解の複素微分干渉画像Ｗｎ（ｒin）に基づいて、下記（１０）式により、位相微分画像を生成することができる。複素微分干渉画像Ｗｎ（ｒout），Ｗｎ（ｒin）の位相が位相微分画像に相当する。
【0097】
【数9】
【0098】
【数10】
【0099】
次に、屈折率分布画像生成部５３による屈折率分布画像の生成処理の詳細について説明する。図１４は、屈折率分布画像生成処理のフローチャートである。
【0100】
ステップＳ８１において、複数の光照射方向それぞれについて、位置ｚに光波面を自由伝搬させて、位置ｚにおける複素振幅画像ｕを求める。例えば、ｚ＝０の位置からｚ＝ｄの位置への光波面の自由伝搬を想定する場合、ｚ＝０の位置の複素振幅画像をｕ（ｘ，ｙ，０）とし、このｕ（ｘ，ｙ，０）の２次元フーリエ変換をＵ（ｋｘ，ｋｙ，０）とし、ｚ＝ｄの位置の複素振幅画像をｕ（ｘ，ｙ，ｄ）とし、このｕ（ｘ，ｙ，ｄ）の２次元フーリエ変換をＵ（ｋｘ，ｋｙ，ｄ）とすると、自由伝搬は下記（１１）式，（１２）式を含む処理により計算で求めることができる。ｉは虚数単位であり、ｋ０は観察対象物中における光の波数である。
【0101】
【数11】
【0102】
【数12】
【0103】
ステップＳ８２において、下記（１３）式により、複素微分干渉画像Ｗ（ｒ）を生成する。δｒはシアーを表す。ステップＳ８３において、この複素微分干渉画像Ｗ（ｒ）の位相で表される位相微分画像を生成する。
【0104】
【数13】
【0105】
ステップＳ８４において、ｚにδｚを加算したものを新たなｚとし、ステップＳ８５において、その新たなｚが最終値ｚendに達したか否かを判定する。新たなｚが最終値ｚendに達していないとステップＳ８５において判定された場合には、ステップＳ８１に戻る。新たなｚが最終値ｚendに達したとステップＳ８５において判定された場合には、ステップＳ８６へ進む。
【0106】
ステップＳ８１からステップＳ８５までの各処理は、ｚの初期値から最終値ｚendまでδｚ刻みで繰り返し行われる。この繰り返し処理により、ｚの初期値から最終値ｚendまでのδｚ刻みの各ｚ位置における位相微分画像、すなわち、３次元位相微分画像が生成される。ステップＳ８５から進んだステップＳ８６において、３次元位相微分画像に基づいて、デコンボリューションにより屈折率分布画像を生成する。
【0107】
次に、第２複素振幅画像生成部５４による複素振幅画像の生成処理の詳細について説明する。第１複素振幅画像生成部５２は、実測データである干渉強度画像に基づいて複素振幅画像を生成する。第３複素振幅画像生成部５６は、トランスミッション行列が複素振幅画像を表現するものであることを利用して、トランスミッション行列の演算に基づいて複素振幅画像を生成する。これに対して、第２複素振幅画像生成部５４は、次のようにして、屈折率分布画像生成部５３により生成された屈折率分布画像を用いて、この屈折率分布画像で表される媒質中における光の伝搬を数値計算で求めることにより複素振幅画像を生成する。
【0108】
各ブロックの屈折率分布を考慮した数値計算により、第１位置（ｚ＝ｚｊ－１）の光波面を該ブロックの内部を伝搬させることで、第２位置（ｚ＝ｚｊ）の光波面を求める。すなわち、複数の光照射方向それぞれについて、第ｊブロックの第１位置（ｚ＝ｚｊ－１）の複素振幅画像および第ｊブロックの屈折率分布に基づいて、第ｊブロックの第２位置（ｚ＝ｚｊ）の複素振幅画像を生成する。この処理に際して、媒質の屈折率分布を考慮して光波面の伝搬を数値計算する手法が用いられる。このような不均一媒質伝搬の数値計算手法として、ＢＰＭ（Beam Propagation Method）、Split-Step Non-Paraxial（ＳＳＮＰ） MethodおよびMulti-layer Bornmultiple-scattering model等が知られている。以下では、ＢＰＭを用いた処理について説明する。
【0109】
図１５は、ＢＰＭの処理内容を説明する図である。この図は、或る一つのブロックを示している。この図に示されるように、光伝搬経路（ｚ方向）に沿った撮像部からの距離に基づいてブロックをＭ個（この図では７個）のスライス（第１～第Ｍのスライス）に区分する。各スライスの厚みは波長程度である。
【0110】
各スライスの厚みは一定であってもよい。ここでは、各スライスの厚みを一定値のΔｚとする。第１～第Ｍのスライスのうちの第ｍスライスは、位置（ｚｊ－１＋(ｍ－１)Δｚ）から位置（ｚｊ－１＋ｍΔｚ）までである。第１位置（ｚ＝ｚｊ－１）から第２位置（ｚ＝ｚｊ）へ向かって順に、第１～第Ｍのスライスにおいて順次に、屈折率分布に応じた位相変化を与えて光波面をΔｚだけ伝搬させる。
【0111】
位置ｚにある厚みΔｚのスライスを伝搬する際に光波面に与えられる位相変化ｏ（ｘ，ｙ，ｚ）は、下記（１４）式で表される。この（１４）式中のｋｖは真空中での光の波数である。δｎ（ｘ，ｙ，ｚ）は、位置ｚにおける観察対象物の屈折率分布ｎ（ｘ，ｙ，ｚ）と背景（媒質）の屈折率ｎｂとの差分であり、下記（１５）式で表される。また、cosθは下記（１６）式で表される。
【0112】
【数14】
【0113】
【数15】
【0114】
【数16】
【0115】
第ｍスライスの位置（ｚ＝ｚｊ－１＋(ｍ－１)Δｚ）における光の複素振幅をｕ（ｘ，ｙ，ｚ）とすると、第ｍスライスの内部を光が伝搬した後の位置（ｚ＋Δｚ）における光の複素振幅ｕ（ｘ，ｙ，ｚ＋Δｚ）は、下記（１７）式で表される。この（１７）式中のＰ（ｋｘ，ｋｙ；Δｚ）は、下記（１８）式で表される。（１７）式は、光の複素振幅ｕ（ｘ，ｙ，ｚ）と位相変化ｏ（ｘ，ｙ，ｚ）との積をフーリエ変換し、このフーリエ変換の結果とＰ（ｋｘ，ｋｙ；Δｚ）との積を逆フーリエ変換することで、厚みΔｚのスライスを伝搬した後の位置（ｚ＋Δｚ）における光の複素振幅ｕ（ｘ，ｙ，ｚ＋Δｚ）を求めることを表している。ＰΔｚは、Δｚの光伝搬の演算を行う関数である。
【0116】
【数17】
【0117】
【数18】
【0118】
ブロックの各スライスにおける光波面の伝搬は、下記（１９）式～（２１）式で表される。すなわち、ブロックの第１位置（ｚ＝ｚｊ－１）における光の複素振幅をｕ（ｘ，ｙ，ｚｊ－１）とすると、ブロックの第１スライスを伝搬した後の光の複素振幅ｕ（ｘ，ｙ，ｚｊ－１＋Δｚ）は、下記（１９）式で表される。ブロックの第(ｍ－１)スライスを伝搬した後の光の複素振幅をｕ（ｘ，ｙ，ｚｊ－１＋(ｍ－１)Δｚ）とすると、ブロックの第ｍスライスを伝搬した後の光の複素振幅ｕ（ｘ，ｙ，ｚｊ－１＋ｍΔｚ）は、下記（２０）式で表される。ブロックの第(Ｍ－１)スライスを伝搬した後の光の複素振幅をｕ（ｘ，ｙ，ｚｊ－１＋(Ｍ－１)Δｚ）とすると、ブロックの第Ｍスライスを伝搬した後の第２位置（ｚ＝ｚｊ）における光の複素振幅ｕ（ｘ，ｙ，ｚｊ）は、下記（２１）式で表される。
【0119】
【数19】
【0120】
【数20】
【0121】
【数21】
【0122】
このようにして、ブロックの屈折率分布を考慮した数値計算により、第１位置（ｚ＝ｚｊ－１）の光波面をブロック内部をスライス毎に順次に伝搬させることで、第２位置（ｚ＝ｚｊ）の光波面を求めることができる。
【0123】
図１６は、第２複素振幅画像生成ステップＳ５４のフローチャートである。ステップＳ９１において、位置ｚを、ブロックの第１位置（ｚ＝ｚｊ－１）に初期化する。ステップＳ９２において、位置ｚにおける光の複素振幅ｕ（ｘ，ｙ，ｚ）と位相変化ｏ（ｘ，ｙ，ｚ）との相互作用を求める。ステップＳ９３において、その相互作用後の光の波面を距離Δｚだけ伝搬させて、位置ｚ＋Δｚにおける光の複素振幅ｕ（ｘ，ｙ，ｚ＋Δｚ）を求める。ステップＳ９４において、Δｚを加算したｚを新たなｚとする。ステップＳ９５において、位置ｚがブロックの第２位置（ｚ＝ｚｊ）に未だ到達していないと判断されれば、ステップＳ９２に戻って、ステップＳ９２～Ｓ９４を繰り返す。ステップＳ９５において、位置ｚがブロックの第２位置（ｚ＝ｚｊ）に到達したと判断されれば、第２複素振幅画像生成ステップＳ５４の処理を終了する。終了時に取得された光の複素振幅がブロックの第２位置（ｚ＝ｚｊ）における複素振幅となる。
【0124】
ところで、干渉強度画像取得ステップＳ５１の際の観察対象物への光照射方向の数と、複素振幅画像の画素数とは、互いに等しいことが理想的である。しかし、現実的には、光照射方向の数は複素振幅画像の画素数より少ない（アンダーサンプリング）。例えば、複素振幅画像の画素数が１０２４×１０２４であるとすると、これと同数の光照射方向を実現するのは、可能ではあっても、容易でない。或いは、撮像部により得られた画像のうち一部範囲の画像（光照射方向の数と同数の画素）のみを以降の処理に用いることも考えられるが、解像度の低下につながるので、好ましくない。
【0125】
入力および出力の双方が位置基底であるトランスミッション行列Ｔrr（ｒout；ｒin）を求める上記（６）式の積分計算は、実際の数値計算の際には下記（２２）式のような総和計算となる。光照射方向の数が複素振幅画像の画素数より少ない場合、観察対象物への光照射方向の走査が図８の（ａ）～（ｃ）に示されるようにｋｘｋｙ平面において離散的かつ周期的であると、理想的なトランスミッション行列Ｔrr（ｒout；ｒin）ではなく、これに替えて下記（２３）式で表されるようなトランスミッション行列Ｔ’rr（ｒout；ｒin）が得られる。（２３）式中のｒｐは、波数空間における複数の光照射方向の位置の周期的分布に応じたものであり、その周期の逆数に相当する。（２３）式のＴ’rr（ｒout；ｒin）は、理想的なＴrr（ｒout；ｒin）をｒｐだけ繰り返してシフトさせて足し合わせたものに相当する。
【0126】
【数22】
【0127】
【数23】
【0128】
図１７は、光照射方向の数が複素振幅画像の画素数と等しい場合の（ａ）ｋｘｋｙ波数空間での光照射方向の分布および（ｂ）或る１点ｒ０に光を集光照射したときのトランスミッション行列Ｔrr（ｒout；ｒin＝ｒ０）のｘｙ空間での分布を示す図である。図１８は、光照射方向の数が複素振幅画像の画素数より少ない場合の（ａ）ｋｘｋｙ波数空間での光照射方向の分布および（ｂ）或る１点ｒ０に光を集光照射したときのトランスミッション行列Ｔ’rr（ｒout；ｒin＝ｒ０）のｘｙ空間での分布を示す図である。図１８中のｒpx，ｒpyは、波数空間における複数の光照射方向の位置の周期的分布に応じたものであり、その周期の逆数に相当する。この図に示されるように、光照射方向の数が複素振幅画像の画素数より少ない場合、多点照射の場合の出力となっている。このような（２３）式のトランスミッション行列Ｔ’rr（ｒout；ｒin）を用いて以降の各ステップの処理を行うと、空間の複数の点の結果が足し合わされてしまうことになって、複素微分干渉画像、位相微分画像および屈折率分布画像を正しく求めることができない。
【0129】
そこで、このようなアンダーサンプリングの場合には、干渉強度画像取得ステップＳ５１において、干渉強度画像取得部５１は、図８の（ａ）～（ｃ）に示されるように波数空間において光照射方向を表す波数ベクトルの位置が離散的かつ周期的に分布する複数の光照射方向それぞれに沿って光が観察対象物に照射されたときに干渉強度画像を撮像した撮像部から、複数の光照射方向それぞれの干渉強度画像を取得する。そして、第３複素振幅画像生成ステップＳ５６において、第３複素振幅画像生成部５６は、複数の光照射方向それぞれの波数ベクトルの波数空間における位置の周期的分布に基づいて区分される複数の領域それぞれにおいてトランスミッション行列Ｔとトランスミッション行列Ｔ１の逆行列との積を求め、これら複数の領域それぞれの積に基づいて複数の光照射方向それぞれの第３複素振幅画像を生成する。
【0130】
すなわち、アンダーサンプリングの場合、第３複素振幅画像生成ステップＳ５６において、図１９に示されるように、トランスミッション行列Ｔ’rr（ｒout；ｒin＝ｒ０）のｘｙ空間を、各照射点を含む領域Ｄ（ｒin－ｒｐ）に分割する。そして、これら複数の領域Ｄ（ｒin－ｒｐ）それぞれでトランスミッション行列Ｔとトランスミッション行列Ｔ１の逆行列との積を求めた後に、複数の領域Ｄ（ｒin－ｒｐ）それぞれで求めた積を貼り合わせることで、全体のトランスミッション行列Ｔとトランスミッション行列Ｔ１の逆行列との積を求める（下記（２４）式）。
【0131】
【数24】
【0132】
本実施形態の観察装置および観察方法は、対物レンズの光軸に垂直な複数の面で観察対象物をブロック分割して、照射側または撮像側から順次に各ブロックの屈折率分布画像を生成することができる。これまでの説明では、照射側から順次に各ブロックの屈折率分布画像を生成する態様について説明してきた。しかし、図２０を用いて説明するとおり、撮像側から順次に各ブロックの屈折率分布画像を生成する態様とすることもできる。また、或る第１期間の測定データに基づいて生成された屈折率分布画像を用いて、他の第２期間の測定データに基づいて第３複素振幅画像を生成する態様も可能である。
【0133】
図２０は、観察方法の他の態様における各画像の生成手順を説明する図である。互いに異なる第１期間および第２期間それぞれにおいて、干渉強度画像取得部５１は、撮像部から複数の光照射方向それぞれの干渉強度画像を取得する。第１複素振幅画像生成部５２は、干渉強度画像取得部５１により取得された複数の光照射方向それぞれの干渉強度画像に基づいて、複数の光照射方向それぞれの第１複素振幅画像を生成する。
【0134】
屈折率分布画像生成部５３は、第１複素振幅画像生成部５２により生成された第１期間の複数の光照射方向それぞれの第１複素振幅画像に基づいて、第１端面から第１中間面までの区間の第１屈折率分布画像を生成する。第２複素振幅画像生成部５４は、屈折率分布画像生成部５３により生成された第１屈折率分布画像に基づいて、複数の光照射方向それぞれについて、平面波の光が第１端面に入射して第１中間面まで到達したときの光の第２複素振幅画像を生成する。
【0135】
トランスミッション行列生成部５５は、第２期間の複数の光照射方向それぞれの第１複素振幅画像に基づいてトランスミッション行列Ｔを生成する。また、トランスミッション行列生成部５５は、第１期間の複数の光照射方向それぞれの第２複素振幅画像に基づいてトランスミッション行列Ｔ１を生成する。第３複素振幅画像生成部５６は、第２期間のトランスミッション行列Ｔと第１期間のトランスミッション行列Ｔ１の逆行列との積に基づいて、複数の光照射方向それぞれについて、平面波の光が平面波の光が第１中間面に入射して第２端面に到達したときの光の第３複素振幅画像を生成する。屈折率分布画像生成部５３は、第３複素振幅画像生成部５６により生成された複数の光照射方向それぞれの第３複素振幅画像に基づいて、第１中間面から第２中間面までの区間の第２屈折率分布画像を生成する。
【0136】
この観察方法の態様は、例えば、観察対象物の中心領域等を速いタイムラプスイメージングをする際に、全体としては屈折率分布が大きくは変化しないと仮定できるときに有効である。第１期間の測定データに基づいて第１端面と第１中間面との間の区間の屈折率分布画像を生成してトランスミッション行列Ｔ１を予め求めておけば、第２期間の測定データに基づいて生成された第１複素振幅画像およびトランスミッション行列Ｔ１に基づいて第３振幅画像を高速に生成することができる。
【0137】
【0138】
シミュレーションＡでは、図２１に示されるように、５種類の位相画像を互いに間隔をあけて並列配置したものを観察対象物として用いてシミュレーションを行った。図２１は、シミュレーションＡ時の配置を模式的に説明する図である。第１複素振幅画像生成部５２により生成される第１複素振幅画像は、５種類の位相画像の全てを通過した後の光の複素振幅画像である。第２複素振幅画像生成部５４により生成される第２複素振幅画像は、５種類の位相画像のうち光照射側の３枚の位相画像を通過した位置での光の複素振幅画像である。第３複素振幅画像生成部５６により生成される第３複素振幅画像は、上記の第１複素振幅画像および第２複素振幅画像に基づくトランスミッション行列の演算により生成されたものであって、５種類の位相画像のうち撮像側の２枚の位相画像を通過した光の複素振幅画像に対応する。
【0139】
図２２は、シミュレーションＡの結果を示す図である。図２２において、最上段は、厳密解の位相微分画像を示し、第２段は、第１複素振幅画像に基づいて生成された入力側空間分解の位相微分画像を示し、第３段は、第２複素振幅画像に基づいて生成された入力側空間分解の位相微分画像を示し、最下段は、第３複素振幅画像に基づいて生成された入力側空間分解の位相微分画像を示す。第１複素振幅画像または第２複素振幅画像に基づいて生成された入力側空間分解の位相微分画像は、撮像部から遠いほど鮮明な位相微分画像が得られた。一方で、第３複素振幅画像に基づいて生成された入力側空間分解の位相微分画像は、撮像部に近いほど鮮明な位相微分画像が得られた。シミュレーションＡでは、第１複素振幅画像、第２複素振幅画像および第３複素振幅画像の何れかに基づいて位相微分画像を生成することで、深達度の高い観察対象物の観察が可能であることが示された。
【0140】
シミュレーションＢでは、図２３に示されるように、細胞塊を模擬したものを観察対象物として用いてシミュレーションを行った。図２３は、シミュレーションＢ時の配置を模式的に説明する図である。図２４および図２５それぞれは、シミュレーションＢの結果を示す図である。
【0141】
図２４において、上段は、第１複素振幅画像のみに基づいて生成された入力側空間分解の位相微分画像を示し、下段は、本実施形態の方法により生成された入力側空間分解の位相微分画像を示す。また、図２４において、左端はｚ軸に平行な断面における位相微分画像を示し、その他はｚ方向の三つの位置それぞれにおけるｘｙ断面における位相微分画像を示す。
【0142】
図２５において、上段は厳密解の屈折率分布画像を示し、中段は、第１複素振幅画像のみに基づいて生成された入力側空間分解の位相微分画像から生成された屈折率分布画像を示し、下段は、本実施形態の方法により生成された入力側空間分解の位相微分画像から生成された屈折率分布画像を示す。また、図２５において、左端はｙ方向への最大値投影画像を示し、その他はｚ方向の三つの位置それぞれにおけるｘｙ断面における屈折率分布画像を示す。
【0143】
このシミュレーションＢでも、第１複素振幅画像、第２複素振幅画像および第３複素振幅画像の何れかに基づいて位相微分画像を生成することで、深達度の高い観察対象物の観察が可能であることが示された。
【0144】
以上のとおり、本実施形態によれば、観察対象物が多重散乱体である場合であっても、良好な深達度を有することができ、多重散乱光の影響を低減して観察対象物を観察することができる。
【0145】
これまで説明してきた観察装置および観察方法は、光源から出力されて２分岐された第１分岐光および第２分岐光のうち、一方の第１分岐光については観察対象物を経ることなく参照光として撮像部により受光し、他方の第２分岐光については観察対象物を経て物体光として撮像部により受光して、第２分岐光（物体光）の複数の光照射方向それぞれについて、第１分岐光（参照光）と第２分岐光（物体光）との干渉による干渉強度画像を撮像し、この干渉強度画像に基づいて複素振幅画像を生成するものであった。
【0146】
これに限られず、以下に図２６および図２７を用いて説明するように、光源から出力されて２分岐された第１分岐光および第２分岐光の双方について観察対象物を経た後に撮像部により受光して、第２分岐光の複数の光照射方向それぞれについて、第１分岐光と第２分岐光との干渉による干渉強度画像を撮像し、この干渉強度画像に基づいて複素振幅画像を生成する構成としてもよい。
【0147】
図２６は、観察装置１００の構成を示す図である。この観察装置１００は、光源１１０、照射部１３１および撮像部１５０等を備える。光源１１０は、空間的にコヒーレントな光を出力する。光源１１０から出力される光は、時間的にコヒーレントであってもよいし、時間的にコヒーレントでなくてもよい。光源１１０は、レーザ光源であってもよいし、例えばＳＬＤ（Super Luminescent Diode）、ＳＣ（Super Continuum）光源および光周波数コム光源などの光源であってよい。また、ＬＥＤ（Light Emitting Diode）や水銀ランプなどから出力された空間インコヒーレントな光を、ピンホール等を通すことで空間的なコヒーレンスを高めてもよい。
【0148】
レンズ１２１は、光源１１０と光学的に接続されており、光源１１０から出力された光を光ファイバ１２３の光入射端１２２に集光して、その光を光入射端１２２に入射させる。光ファイバ１２３は、光入射端１２２に入射された光を光出射端１２４まで導光する。光ファイバ１２３により導光された光は光出射端１２４から発散光として出射される。レンズ１２５は、光出射端１２４と光学的に接続されており、光出射端１２４から発散光として出力された光を入力してコリメートし、そのコリメートした光を照射部１３１へ出力する。
【0149】
照射部１３１は、光源１１０から出力されてレンズ１２１，光ファイバ１２３およびレンズ１２５を経た光を入力し、その入力した光を２分岐して第１分岐光および第２分岐光とする。照射部１３１は、これら第１分岐光および第２分岐光を互いに重ね合わせて観察対象物Ｓに照射する。照射部１３１は、観察対象物Ｓに対して第１分岐光を一定の光照射方向に沿って照射し、観察対象物Ｓに対して第２分岐光を複数の光照射方向それぞれに沿って照射する。
【0150】
照射部１３１は、ビームスプリッタ３１１、位相変調型の空間光変調器３１３、偏光子３１４、２分の１波長板３１５、偏光子３１６、レンズ３１８および対物レンズ３１９を含む。
【0151】
ビームスプリッタ３１１は、レンズ１２５との間に設けられた偏光子３１４および２分の１波長板３１５を経て到達した光を空間光変調器３１３へ反射させる。また、ビームスプリッタ３１１は、空間光変調器３１３から到達した光を入力して、この光を偏光子３１６へ出力する。
【0152】
空間光変調器３１３は、変調面に入射される互いに直交する第１方位および第２方位それぞれの直線偏光の光のうち、第１方位の直線偏光の光に対して位相変調することなく、第２方位の直線偏光の光に対して選択的に位相変調する。偏光子３１４および２分の１波長板３１５は、ビームスプリッタ３１１から空間光変調器３１３の変調面へ入射される光が第１方位および第２方位それぞれの直線偏光成分を互いに同程度に含むように、光の偏光状態を設定する。
【0153】
偏光子３１６は、空間光変調器３１３からビームスプリッタ３１１を経て到達した光を入力し、その光に含まれる第１方位および第２方位それぞれの直線偏光の光を干渉可能とする。レンズ３１８および対物レンズ３１９は、偏光子３１６から出力された第１分岐光および第２分岐光それぞれを平面波として観察対象物Ｓに照射する。
【0154】
このような構成を有する照射部１３１は、空間光変調器３１３により位相変調されなかった第１方位の直線偏光の光を第１分岐光として、この第１分岐光を観察対象物Ｓに対して一定の光照射方向に沿って照射することができる。照射部１３１は、空間光変調器３１３により位相変調された第２方位の直線偏光の光を第２分岐光として、この第２分岐光を観察対象物Ｓに対して複数の光照射方向それぞれに沿って照射することができる。観察対象物Ｓへの第２分岐光の光照射方向は、空間光変調器３１３の変調面における位相変調パターンの向き及び間隔により設定することができる。また、第１分岐光と第２分岐光との間の位相差は、空間光変調器３１３の変調面における位相変調パターンのシフトにより設定することができる。
【0155】
なお、空間光変調器３１３の位置により位相差を調整することもできる。しかし、空間光変調器３１３の変調面における位相変調パターンのシフトにより位相差を設定する場合の方が、構成要素の機械的な移動がない点で好ましい。
【0156】
対物レンズ１４１は、照射部１３１により観察対象物Ｓに照射されて観察対象物Ｓを経た光（第１分岐光、第２分岐光）を入力し、その光をミラー１４２へ出力する。レンズ１４３は、対物レンズ１４１から出力されてミラー１４２により反射された光を入力し、その光を撮像部１５０の撮像面に入射させる。
【0157】
撮像部１５０は、レンズ１４３から撮像面に到達した第１分岐光および第２分岐光の双方を受光して、第１分岐光と第２分岐光との干渉による干渉強度画像を撮像する。撮像部１５０は、第２分岐光の複数の光照射方向それぞれについて、第１分岐光と第２分岐光との間の位相差が複数の位相差それぞれに設定されたときの干渉強度画像を撮像する。この撮像部１５０により撮像された干渉強度画像に基づいて所要の処理をすることで、複素振幅画像などを生成することができる。
【0158】
この観察装置１００（図２６）を用いた場合の複素振幅画像生成部の処理の内容は次のとおりである。図２７は、観察装置１００において、観察対象物Ｓへの第１分岐光および第２分岐光の入射、ならびに、観察対象物Ｓを経た後の撮像部１５０への第１分岐光および第２分岐光の入射、を模式的に示す図である。照射部１３１は、観察対象物Ｓに対して第１分岐光および第２分岐光を重ねて照射する。このとき、観察対象物Ｓに対する第１分岐光の光照射方向を一定とし、観察対象物Ｓに対する第２分岐光の光照射方向を複数の光照射方向それぞれとし、第１分岐光と第２分岐光との間の位相差φを各値に設定する。
【0159】
観察対象物Ｓへ入射する第１分岐光の波面をｕ０,in（ｒ）と表す。第２分岐光の複数（Ｎ個）の光照射方向のうち第ｎ光照射方向（ｎ＝１～Ｎ）に沿って観察対象物Ｓへ入射する第２分岐光の波面をｕｎ,in（ｒ）exp（ｉφ）と表す。ｒは位置を表す変数である。φは、第１分岐光と第２分岐光との間の位相差である。撮像部１５０の撮像面または焦点面（撮像面に対して光学的に共役な面）における第１分岐光の波面をｕ０（ｒ）と表し、第２分岐光の波面をｕｎ（ｒ）exp（ｉφ）と表す。
【0160】
撮像部１５０による撮像により取得される干渉強度画像Ｉｎ（ｒ，φ）は、ｕ０（ｒ）とｕｎ（ｒ）exp（ｉφ）との和の絶対値の二乗で表される。干渉強度画像Ｉｎ（ｒ，φ）は、第１分岐光と第２分岐光との間の位相差をφとして、観察対象物に対して第１分岐光を一定の光照射方向に沿って観察対象物Ｓへ入射させ、観察対象物に対して第２分岐光を第ｎ光照射方向に沿って入射させたときに、撮像部１５０による撮像により取得される干渉強度画像である。焦点面（撮像面に対して光学的に共役な面）は、観察対象物Ｓにあってもよいし、観察対象物Ｓより撮像部１５０側にあってもよいし、観察対象物Ｓより照射部１３１側にあってもよい。
【0161】
第２分岐光の複数の光照射方向それぞれについて、複数の位相差φそれぞれに設定されたときに撮像部１５０により取得された干渉強度画像に基づいて、位相シフト法により干渉項Ｃｎ（ｒ）＝ｕ０
＊（ｒ）・ｕｎ（ｒ）を求める。干渉項ｕ０（ｒ）・ｕｎ
＊（ｒ）を求めてもよい。この干渉項Ｃｎ（ｒ）を第２分岐光の複数の光照射方向それぞれについて（すなわち、各ｎ（＝１～Ｎ））について求める。
【0162】
第２分岐光の複数の光照射方向それぞれについて求めた干渉項Ｃｎ（ｒ）に基づいて、第１分岐光の複素振幅画像を生成する。第１分岐光の複素振幅ｕ０（ｒ）の位相φ０（ｒ）は、第１分岐光と第２分岐光との間の位相傾き（光入射方向の差異）を補正した後、その補正後の干渉項のコヒーレント和Ｃsum（ｒ）を求め、このコヒーレント和Ｃsum（ｒ）の位相で近似的に表すことができる。第１分岐光の複素振幅ｕ０（ｒ）の振幅Ａ０（ｒ）は、第２分岐光を観察対象物Ｓに照射せずに、第１分岐光のみを観察対象物Ｓに照射したときに、撮像部１５０により撮像された強度画像｜ｕ０（ｒ）｜２から求めることができる。または、第１分岐光の複素振幅ｕ０（ｒ）の振幅Ａ０（ｒ）は、干渉項Ｃｎ（ｒ）の強度和Ｉsum（ｒ）の平方根で近似的に表すことができる。
【0163】
以上のようにして求めた第１分岐光の複素振幅ｕ０（ｒ）の位相φ０（ｒ）および振幅Ａ０（ｒ）に基づいて、第１分岐光の複素振幅画像ｕ０（ｒ）を生成することができる。そして、第１分岐光の複素振幅画像ｕ０（ｒ）および干渉項Ｃｎ（ｒ）に基づいて、複数の光照射方向それぞれの第２分岐光の複素振幅画像ｕｎ（ｒ）を生成することができる。以降の処理は、既に説明したものと同様である。
【符号の説明】
【0164】
１Ａ～１Ｃ…観察装置、２…記録媒体、１１…光源、１２…レンズ、１３…光入射端、１４…光ファイバ、１５…ファイバカプラ、１６，１７…光ファイバ、１８，１９…光出射端、２１…レンズ、２２…ミラー、２３…レンズ、２４…コンデンサレンズ、２５…対物レンズ、３１…レンズ、３２…ミラー、３３…駆動部、３４…レンズ、４１…ビームスプリッタ、４２…レンズ、４３…撮像部、５０…解析部、５１…干渉強度画像取得部、５２…第１複素振幅画像生成部、５３…屈折率分布画像生成部、５４…第２複素振幅画像生成部、５５…トランスミッション行列生成部、５６…第３複素振幅画像生成部、５７…表示部、５８…記憶部、１００…観察装置、１１０…光源、１３１…照射部、１５０…撮像部。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
【図26】
【図27】
