(57)【要約】【課題】デジタルホログラフィ３次元撮像においてノイズを低減する。【解決手段】撮像装置は、光源１０１からの光を光路分割素子１０２により物体光と参照光とに分割し、物体を経た物体光と参照光路１１０を通った参照光とを光路合成素子１０８を介して干渉させることで撮像素子１０９上に干渉縞を形成する光学系を有し、該干渉縞の撮像を物体光の物体に対する入射角を変化させて複数回行い、該複数回の撮像により得られた複数の干渉縞から入射角ごとの透過波面を求める波面再生処理と、入射角ごとの透過波面から３次元の屈折率分布を算出する３次元再構成処理とを行う。光源から光路分割素子までの光路、参照光路および光路合成素子から撮像素子までの光路のうちいずれかに、光の位相分布を変化させることが可能な位相分布変更手段２０１を配置する。少なくとも１回の撮像において位相分布変更手段を用いて位相分布を変化させる。【選択図】図２
2015118290.tif 000002
【特許請求の範囲】
【請求項１】
光源からの光を光路分割素子により物体光と参照光とに分割し、物体を経た前記物体光と参照光路を通った前記参照光とを光路合成素子を介して干渉させることで撮像素子上に干渉縞を形成する光学系を有し、該干渉縞の撮像を前記物体光の前記物体に対する入射角を変化させて複数回行い、
該複数回の撮像により得られた複数の前記干渉縞から前記入射角ごとの透過波面を求める波面再生処理と、前記入射角ごとの前記透過波面から３次元の屈折率分布を算出する３次元再構成処理とを行う撮像装置であって、
前記光源から前記光路分割素子までの光路、前記参照光路、および前記光路合成素子から前記撮像素子までの光路のうちいずれかに、光の位相分布を変化させることが可能な位相分布変更手段を配置し、 
前記複数回の撮像のうち少なくとも１回において前記位相分布変更手段を用いて前記位相分布を変化させることを特徴とする撮像装置。
【請求項２】
前記位相分布変更手段は、光を拡散させる光拡散素子であることを特徴とする請求項１に記載の撮像装置。 
【請求項３】
前記位相分布変更手段は、前記参照光路の光路または光路長を変化させる手段であることを特徴する請求項１に記載の撮像装置。
【請求項４】
前記複数回の撮像における撮像ごとに前記位相分布変更手段を用いて前記位相分布を変化させることを特徴とする請求項１から３のいずれか一項に記載の撮像装置。
【請求項５】
光源からの光を光路分割素子により物体光と参照光とに分割し、物体を経た前記物体光と参照光路を通った前記参照光とを光路合成素子を介して干渉させることで撮像素子上に干渉縞を形成する光学系を有し、該干渉縞の撮像を前記物体光の前記物体に対する入射角を変化させて複数回行い、
該複数回の撮像により得られた複数の前記干渉縞から前記入射角ごとの透過波面を求める波面再生処理と、前記入射角ごとの前記透過波面から３次元の屈折率分布を算出する３次元再構成処理とを行う撮像装置であって、
前記光学系の収差を変化させることが可能な収差変更手段を有し、
前記複数回の撮像のうち少なくとも１回において前記収差変更手段を用いて前記収差を変化させ、 
前記波面再生処理または前記３次元再構成処理において前記収差を補正する演算を行うことを特徴とする撮像装置。 
【請求項６】
前記収差変更手段は、前記物体の位置、前記撮像素子の位置または前記光学系の位置を変化させることが可能な駆動装置であることを特徴とする請求項５に記載の撮像装置。
【請求項７】
前記収差変更手段は、前記光学系に配置される空間変調器であることを特徴とする請求項５に記載の撮像装置。 
【請求項８】
前記収差変更手段は、前記光学系に配置される瞳フィルタであることを特徴とする請求項５に記載の撮像装置。 
【請求項９】
前記複数回の撮像における撮像ごとに前記収差変更手段を用いて前記位相分布を変化させることを特徴とする請求項５から８のいずれか一項に記載の撮像装置。
【請求項１０】
光源からの光を光路分割素子により物体光と参照光とに分割し、物体を経た前記物体光と参照光路を通った前記参照光とを光路合成素子を介して干渉させることで撮像素子上に干渉縞を形成する光学系を有し、該干渉縞の撮像を前記物体光の前記物体に対する入射角を変化させて複数回行い、
該複数回の撮像により得られた複数の前記干渉縞から前記入射角ごとの透過波面を求める波面再生処理と、前記入射角ごとの前記透過波面から３次元の屈折率分布を算出する３次元再構成処理とを行う撮像装置であって、
前記波面再生処理または前記３次元再構成処理において、前記光学系が有する既知の収差を補正する演算を行うことを特徴とする撮像装置。
【請求項１１】
光源からの光を光路分割素子により物体光と参照光とに分割し、物体を経た前記物体光と参照光路を通った前記参照光とを光路合成素子を介して干渉させることで撮像素子上に干渉縞を形成する光学系を用いて、該干渉縞の撮像を前記物体光の前記物体に対する入射角を変化させて複数回行い、
該複数回の撮像により得られた複数の前記干渉縞から前記入射角ごとの透過波面を求める波面再生処理と、前記入射角ごとの前記透過波面から３次元の屈折率分布を算出する３次元再構成処理とを行う撮像方法であって、
前記複数回の撮像のうち少なくとも１回において、前記光源から前記光路分割素子までの光路、前記参照光路、および前記光路合成素子から前記撮像素子までの光路のうちいずれかにおける光の位相分布を変化させることを特徴とする撮像方法。 
【請求項１２】
光源からの光を光路分割素子により物体光と参照光とに分割し、物体を経た前記物体光と参照光路を通った前記参照光とを光路合成素子を介して干渉させることで撮像素子上に干渉縞を形成する光学系を用いて、該干渉縞の撮像を前記物体光の前記物体に対する入射角を変化させて複数回行い、
該複数回の撮像により得られた複数の前記干渉縞から前記入射角ごとの透過波面を求める波面再生処理と、前記入射角ごとの前記透過波面から３次元の屈折率分布を算出する３次元再構成処理とを行う撮像方法であって、
前記複数回の撮像のうち少なくとも１回において前記光学系の収差を変化させ、 
前記波面再生処理または前記３次元再構成処理において前記収差を補正する演算を行うことを特徴とする撮像方法。 
【請求項１３】
光源からの光を光路分割素子により物体光と参照光とに分割し、物体を経た前記物体光と参照光路を通った前記参照光とを光路合成素子を介して干渉させることで撮像素子上に干渉縞を形成する光学系を用いて、該干渉縞の撮像を前記物体光の前記物体に対する入射角を変化させて複数回行い、
該複数回の撮像により得られた複数の前記干渉縞から前記入射角ごとの透過波面を求める波面再生処理と、前記入射角ごとの前記透過波面から３次元の屈折率分布を算出する３次元再構成処理とを行う撮像方法であって、
前記波面再生処理または前記３次元再構成処理において、前記光学系が有する既知の収差を補正する演算を行うことを特徴とする撮像方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、デジタルホログラフィを用いた顕微鏡装置等の３次元撮像装置に関する。
【背景技術】
【０００２】
医学や生物学等の分野において生体試料を染色することなく顕微鏡装置により観察する方法として、非特許文献１には、デジタルホログラフィを利用する方法が開示されている。デジタルホログラフィは、試料を透過した透過光と参照光との干渉により生じた干渉縞を撮像素子でデジタル信号として取得し、該デジタル信号に対する計算機による後処理によって透過光の波面を算出するものである。この方法を用いることで試料による波面の位相変化を定量的に求めることができるため、染色を行うことなく試料を高いコントラストで観察することができる。
【０００３】
また、一般に生体試料は複雑な立体構造を持っており、医師等の観察者はその立体構造に基づいて病気の診断を行う。特に癌等では、細胞の立体的な配置が不規則になるため、迅速な診断を行うためにはこの立体構造を素早く認識する必要がある。このためには、試料の３次元的な情報を取得し、それを観察者に明確に提供することが望まれる。これに対し、非特許文献２には、試料を異なる光入射角で照明して撮像素子により取得した複数の干渉縞の画像から、再構成処理によって試料の３次元的な屈折率分布を取得する３次元撮像方法が開示されている。
【０００４】
干渉縞を生成するためには通常はレーザ等の干渉性の高い光を用いるが、干渉性が高い光では、該光が通過する光学素子が有する欠陥やその光学素子の表面粗さ等によって光にランダムな位相分布が生じ、これにより本来一様な照度分布が不均一な分布となる。これは、スペックルまたはスペックルノイズと呼ばれ、像の劣化を招く。このようなスペックルノイズを低減する方法として、非特許文献３には、拡散板を移動させることでスペックルを平均化する方法が開示されている。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】C. D. Depeursinge, E. Cuche, P. Marquet, T. Colomb, P. Dahlgren, A. M. Marian, F. Montfort and P. J. Magistretti 共著, ”Digital Holography Applied to Microscopy” Proceedings of SPIE, Vol.4659, Practical Holography XVI and Holographic Materials VIII, pp.30-34, June 3, 2002, USA
【非特許文献２】C. Fang-Yen, W. Choi, Y. Sung, C.  J. Holbrow, R. R. Dasari,  M. S. Feld, ”Video-rate tomographic phase microscopy” Journal of Biomedical Optics, Vol.16(1),  pp.011005-1 - 011005-5, 14 January, 2011, USA
【非特許文献３】J. W. Goodman, ”Speckle Phenomena In Optics” Roberts and Company Publishers, Chapter 5, USA
【発明の概要】
【発明が解決しようとする課題】
【０００６】
非特許文献２にて開示された３次元撮像方法について、本願発明の発明者が検討したところ、撮像素子と共役な面（焦点面）にノイズが集中的に発生し、これが取得される屈折率分布に影響するという課題を新たに発見した。非特許文献２においても、これに類似する課題として、物体光の入射角を変えた際に固定パターンのノイズが出現することが指摘されているが、焦点面に集中的に発生するノイズに関する記載はない。焦点面には観察者が求める試料の情報が最も多く含まれており、焦点面で画像が劣化することは好ましくない。また、この３次元撮像方法では複数枚の干渉縞の画像を取得する必要があるため、本質的にデータ取得に時間を要する。
【０００７】
さらに、非特許文献３にて開示されたスペックルを単純に平均化する方法では、撮像時間の増大を招くため、ノイズ低減手法としては好ましくない。
【０００８】
本発明は、デジタルホログラフィを用いた３次元撮像を行う際に焦点面近傍に発生するノイズを低減できるようにした３次元撮像装置および３次元撮像方法を提供する。
【課題を解決するための手段】
【０００９】
本発明の一側面としての撮像装置は、光源からの光を光路分割素子により物体光と参照光とに分割し、物体を経た物体光と参照光路を通った参照光とを光路合成素子を介して干渉させることで撮像素子上に干渉縞を形成する光学系を有し、該干渉縞の撮像を物体光の物体に対する入射角を変化させて複数回行い、該複数回の撮像により得られた複数の干渉縞から入射角ごとの透過波面を求める波面再生処理と、入射角ごとの透過波面から３次元の屈折率分布を算出する３次元再構成処理とを行う。該撮像装置は、光源から光路分割素子までの光路、参照光路および光路合成素子から撮像素子までの光路のうちいずれかに、光の位相分布を変化させることが可能な位相分布変更手段を配置する。そして、複数回の撮像のうち少なくとも１回において位相分布変更手段を用いて位相分布を変化させることを特徴とする。
【００１０】
本発明の他の一側面としての撮像装置は、光源からの光を光路分割素子により物体光と参照光とに分割し、物体を経た物体光と参照光路を通った参照光とを光路合成素子を介して干渉させることで撮像素子上に干渉縞を形成する光学系を有し、該干渉縞の撮像を物体光の物体に対する入射角を変化させて複数回行い、該複数回の撮像により得られた複数の干渉縞から入射角ごとの透過波面を求める波面再生処理と、入射角ごとの透過波面から３次元の屈折率分布を算出する３次元再構成処理とを行う。該撮像装置は、光学系の収差を変化させることが可能な収差変更手段を有し、複数回の撮像のうち少なくとも１回において収差変更手段を用いて収差を変化させ、 波面再生処理または３次元再構成処理において収差を補正する演算を行うことを特徴とする。
【００１１】
本発明の他の一側面としての撮像装置は、光源からの光を光路分割素子により物体光と参照光とに分割し、物体を経た物体光と参照光路を通った参照光とを光路合成素子を介して干渉させることで撮像素子上に干渉縞を形成する光学系を有し、該干渉縞の撮像を物体光の物体に対する入射角を変化させて複数回行い、該複数回の撮像により得られた複数の干渉縞から入射角ごとの透過波面を求める波面再生処理と、入射角ごとの透過波面から３次元の屈折率分布を算出する３次元再構成処理とを行う。該撮像装置は、波面再生処理または３次元再構成処理において、光学系が有する既知の収差を補正する演算を行うことを特徴とする。
【００１２】
さらに、本発明の他の一側面としての撮像方法は、光源からの光を光路分割素子により物体光と参照光とに分割し、物体を経た物体光と参照光路を通った参照光とを光路合成素子を介して干渉させることで撮像素子上に干渉縞を形成する光学系を用いて、該干渉縞の撮像を物体光の物体に対する入射角を変化させて複数回行い、該複数回の撮像により得られた複数の干渉縞から入射角ごとの透過波面を求める波面再生処理と、入射角ごとの透過波面から３次元の屈折率分布を算出する３次元再構成処理とを行う。該撮像方法は、複数回の撮像のうち少なくとも１回において、光源から前記光路分割素子までの光路、参照光路、および光路合成素子から撮像素子までの光路のうちいずれかにおける光の位相分布を変化させることを特徴とする。 
【００１３】
本発明の他の一側面としての撮像方法は、光源からの光を光路分割素子により物体光と参照光とに分割し、物体を経た物体光と参照光路を通った参照光とを光路合成素子を介して干渉させることで撮像素子上に干渉縞を形成する光学系を用いて、該干渉縞の撮像を物体光の物体に対する入射角を変化させて複数回行い、該複数回の撮像により得られた複数の干渉縞から入射角ごとの透過波面を求める波面再生処理と、入射角ごとの透過波面から３次元の屈折率分布を算出する３次元再構成処理とを行う。該撮像方法は、複数回の撮像のうち少なくとも１回において光学系の収差を変化させ、波面再生処理または３次元再構成処理において収差を補正する演算を行うことを特徴とする。
【００１４】
本発明の他の一側面としての撮像方法は、光源からの光を光路分割素子により物体光と参照光とに分割し、物体を経た物体光と参照光路を通った参照光とを光路合成素子を介して干渉させることで撮像素子上に干渉縞を形成する光学系を用いて、該干渉縞の撮像を物体光の物体に対する入射角を変化させて複数回行い、該複数回の撮像により得られた複数の干渉縞から入射角ごとの透過波面を求める波面再生処理と、入射角ごとの透過波面から３次元の屈折率分布を算出する３次元再構成処理とを行う。該撮像方法は、波面再生処理または３次元再構成処理において、光学系が有する既知の収差を補正する演算を行うことを特徴とする。
【発明の効果】
【００１５】
本発明によれば、デジタルホログラフィを用いた３次元撮像を行う際に焦点面近傍に出現するノイズを良好に低減することができ、高精度の３次元屈折率分布を得ることができる。
【図面の簡単な説明】
【００１６】
【図１】本発明の実施例１〜４に共通するデジタルホログラフィ３次元撮像装置の構成を示す図。
【図２】実施例１の撮像装置の構成を示す図。
【図３】実施例１における入射角の定義を示す図。
【図４】参照光がスペックルを有する場合の３次元屈折率分布を示す図。
【図５】実施例１で取得される３次元屈折率分布を示す図。
【図６】実施例２の撮像装置の構成を示す図。
【図７】実施例２で取得される３次元屈折率分布を示す図。
【図８】実施例３の撮像装置の構成を示す図。
【図９】実施例３で取得される３次元屈折率分布を示す図。
【図１０】実施例４の撮像装置の構成を示す図。
【図１１】実施例４で取得される３次元屈折率分布を示す図。
【発明を実施するための形態】
【００１７】
以下、本発明の実施例について図面を参照しながら説明する。
【００１８】
まず、後述する具体的な実施例１〜４において共通する事項について説明する。各実施例では、算出される透過波面に含まれる参照光の位相を物理的または数値的に物体（試料）への入射角ごとに変化させることで、焦点面近傍にて発生するノイズを低減する。図１には、各実施例にて説明するノイズ低減方法が適用されるデジタルホログラフィ３次元撮像装置の光学系（以下、撮像光学系という）の基本的な構成を示している。本撮像装置は、参照光と物体光が異なる角度で撮像素子に入射する、いわゆるｏｆｆ−ａｘｉｓ型の撮像光学系を採用する。
【００１９】
レーザ光源１０１から発せれられたレーザ光は、ビームスプリッタ（光路分割素子）１０２によって参照光１２２と物体光１２１とに分割される。物体光１２１は、ガルバノミラー１０３によって反射され、コンデンサーレンズ１０４を通して物体としての試料１０５に照射される。ガルバノミラー１０３は、試料１０５と共役な面に配置されており、ガルバノミラー１０３の角度を変化させることで物体光１２１の試料１０５に対する入射角が変化する。試料１０５を経た（透過した）物体光１２１は、対物レンズ１０６とリレーレンズ１０７を透過した後、ビームスプリッタ（光路合成素子）１０８を介して撮像素子１０９に到達する。
【００２０】
一方、参照光１２２は、参照光路１１０を通り、ビームスプリッタ１０８を介して、試料１０５を透過した後の物体光１２１と同じ光路に導かれて撮像素子１０９に到達する。物体光１２１と参照光１２２は撮像素子１０９の撮像面上（撮像素子上）で干渉して干渉縞を形成する。この干渉縞は、撮像素子１０９による撮像（光電変換）によってデジタル信号に変換され、不図示の計算機（コンピュータ）に保存される。保存された干渉縞のデータを、以下の説明では、ホログラムデータという。各実施例では、物体光１２１の試料１０５に対する入射角（以下、単に物体光の入射角という）を変化させながら撮像を複数回行い、複数のホログラムデータ（画像）を取得する。
【００２１】
そして、計算機は、該複数のホログラムデータから透過波面を求め、その後３次元的な屈折率の分布を算出する。透過波面の算出方法は複数存在するが、各実施例では、高速フーリエ変換（ＦＦＴ）を用いた方法を採用する。具体的には、まず計算機は、ホログラムデータに対してＦＦＴを行ってホログラムデータの空間スペクトルを得る。次に、計算機は、得られた空間スペクトルから実像の空間スペクトルを抽出することで試料１０５を透過した波面の空間スペクトルを得る。さらに、計算機は、その空間スペクトルのピークが原点になるように空間スペクトル全体をシフトさせることにより、試料１０５への斜入射の影響を除去する。このとき、空間スペクトルに適当な位相分布を与えることで撮像光学系が有する収差や試料１０５の撮像位置を補正することも可能である。計算機は、シフトされた空間スペクトルに逆ＦＦＴを行うことで透過波面を算出する。この過程を、以下の説明では、波面再生処理（または単に波面再生）と称する。計算機は、物体光の入射角を変化させて取得したホログラムデータのすべてに波面再生を行って、すべての入射角に対する透過波面を算出する。
【００２２】
複数の透過波面から３次元的な屈折率分布を算出する方法も種々存在するが、各実施例ではＲｙｔｏｖ近似に基づく逆伝搬法を用いる。この逆伝搬法では、まず計算機は、透過波面の対数、すなわち位相分布を求め、その空間スペクトルを算出する。そして、計算機は、得られた空間スペクトルに適切な係数を乗じた後、空間周波数空間で半径がｎｅ／λであり、中心位置が（−ｆｘ０，−ｆｙ０，−ｆｚ０）であるＥｗａｌｄ球の球面上にスペクトルを貼り付ける。ｎｅは媒質の屈折率、λはレーザ光の波長、ｆｘ０，ｆｙ０，ｆｚ０は試料に入射する物体光の空間周波数である。
【００２３】
計算機は、この工程をすべての透過波面に対して行い、得られたすべてのＥｗａｌｄ球の和を算出することで屈折率分布の３次元的な空間スペクトルを得る。さらに、計算機は、この３次元空間スペクトルを逆フーリエ変換することで３次元的な屈折率分布を算出する。このとき、逆フィルタ等を用いることで、屈折率分布を鮮鋭化することも可能である。このようにして複数の透過波面から３次元の屈折率分布を算出することを、以下の説明において３次元再構成処理（または単に３次元再生）と称する。
【００２４】
ホログラムデータを取得するためには、通常レーザ光等、干渉性の高い光を用いる。干渉性が高い光を用いると、この光が入射する光学素子が有する欠陥や該光学素子の表面粗さ等によって該光にランダムな位相分布が生じる。このランダムな位相分布によって、本来一様な照度分布が不均一な分布となる。前述したように、この不均一な照度分布をスペックルまたはスペックルノイズという。スペックルは、取得するホログラムデータの劣化を招き、さらには最終的に得られる３次元の屈折率分布を劣化させる。
【００２５】
スペックルの分布は、光の伝搬や介在する光学素子の変化等に伴う光の位相分布の変動によって変化する。参照光および物体光はともにスペックルを有するが、物体光の入射角が変わるごとに、該物体光が通過する光路長や該物体光の光学素子への入射位置が変化するため、スペックルの分布も物体光の入射角に応じて変化する。しかし、参照光は、撮像中において常に一定の光路を通るため、スペックルの分布は変化しない。その結果、ホログラムデータは、物体光の入射角にかかわらず常に一定のスペックルを有することになる。このことは、波面再生の結果として得られる透過波面に、物体光の入射角に依存しない一定のスペックルが付加されることを意味する。
【００２６】
さらに、透過波面が物体光の入射角に依存しない一定のスペックルを有すると、３次元再構成の結果得られる３次元の屈折率分布は、撮像素子と共役な面上に集中的にノイズを有することになる。この理由は以下の通りである。
【００２７】
先に述べた３次元再構成は、大きく分けると、位相分布の空間スペクトルをＥｗａｌｄ球へ貼り付ける第１のステップと、各入射角に対して算出したＥｗａｌｄ球を足し合わせる第２のステップとに分けられる。このうち第１のステップは、実空間で考えると、位相分布を入射方向に沿って伝搬させる過程と考えることができる。撮像素子と共役な面では、この演算によってデータが変化しないため、透過波面が有する一定のスペックルも物体光の入射角にかかわらず必ず同じ分布のままとなる。
【００２８】
しかし、撮像素子と共役な面以外の領域では、位相分布を入射方向に沿って伝搬させることにより、スペックルは物体光の入射角に依存して異なる分布を有することになる。また、第２のステップは、伝搬したそれぞれの位相分布を積算する過程と考えることができる。これにより、撮像素子と共役な面以外の領域では、異なる分布のスペックルを積算することになり、その結果としてスペックルが低減される。しかし、撮像素子と共役な面ではスペックルは物体光の入射角にかかわらず同じ分布であるため、積算してもスペックルが低減されることはない。このようにして、焦点面近傍にのみ集中的にスペックルが発生するという現象が発生する。以下、このスペックルノイズを、焦点ノイズと称する。
【００２９】
焦点ノイズは、参照光のスペックルが物体光の入射角によらず常に一定であることに起因している。したがって、焦点ノイズを低減するためには、物体光の入射角を変化させる際に参照光の位相に変化を与えてスペックルの分布を変化させればよい。これを実現するためには、例えば参照光路１１０に入射した光を拡散させる光拡散板等の光拡散素子を配置し、物体光の入射角の変化に合わせてこの光拡散素子を移動させればよい。
【００３０】
また、参照光路１１０に、ガラス製の平行平板等、光路長を変更可能な素子を参照光路に対して挿入または抜き取ることで参照光路長を変化させてもよい。他にも、例えば参照光路１１０に配置したミラーの角度や位置を変化させることで、参照光路の光路を変化させることでも実現できる。
【００３１】
ただし、これらはスペックルの分布を変化させるための手段の具体例であり、本質的には、参照光の光路に光の位相分布を変化させることが可能な位相分布変更手段を配置し、これを用いて物体光の入射角の変化に伴い参照光のスペックル分布を変化させればよい。
【００３２】
さらに、別の方法によってスペックルの分布を変化させてもよい。先に説明したように、各実施例では、波面再生処理や３次元再構成処理等の演算処理を行う。このため、後の数値演算時に参照光のスペックル分布を変化させてもよい。これを実現する方法としては、物体光の入射角の変化に合わせて試料を移動させ、位置を数値的に補正する方法が挙げられる。例えば、物体光の入射角ごとに試料を光軸方向の異なる位置に移動させる。このとき、試料の像は光軸方向に移動するが、参照光の光路に変化はないため、参照光のスペックル分布は変化しない。
【００３３】
デジタルホログラフィでは、後の演算によって透過波面を任意の撮像位置での波面に変換することができるため、取得した透過波面から同じ撮像位置で透過波面を算出することができる。この演算は、取得された透過波面、すなわち参照光のスペックルを含む波面に対して行われるため、参照光のスペックル分布も同時に変化する。これは、入射角ごとに参照光の位相分布を変化させたことと等価であるため、先に述べたように、この後に３次元再構成を行えば、焦点ノイズが低減される。
【００３４】
試料を移動させることは、撮像光学系に波面収差を与えることと等価である。例えば、横方向への移動はチルト収差を与え、光軸方向への移動はデフォーカス収差を与えることに相当する。このことから、上述した方法は、物体光の入射角ごとに異なる収差を与え、波面再生時に数値的に収差補正を行うことと等価と考えることもできる。したがって、撮像光学系に、その収差を変化させることが可能な収差変更手段を導入し、物体光の入射角ごとに該光学素子によって収差を変化させ、波面再生時または３次元再構成時に収差を補正する演算を行うことでも焦点ノイズを低減することができる。
【００３５】
なお、撮像光学系に収差を与えれば、物体光の入射角ごとに収差を変更する必要はない。これは、物体光の入射角を変化させると、透過波面の回折光は対物レンズ１０６の異なる位置を通ることになるため、物体光の入射角によって透過波面が受ける収差が異なるからである。したがって、例えば、撮像光学系の瞳面で回転方向に値が変化するような収差、例えば非点収差を与え、波面再生時にこれを補正すれば、撮像ごとに光学素子を移動させることなく、焦点ノイズを低減することができる。
【００３６】
各実施例では、物体光の入射角を変化させるとともに参照光の位相を物理的または数値的に変化させることを本質とする具体的手段を採用することで、焦点ノイズを低減し、高精度の３次元屈折率分布を得る。
【００３７】
また、先に述べた非特許文献３にて開示された方法では、１つの入射角に対して平均化を行うことになるため、露光時間の増大や撮像回数の増加を招くが、各実施例ではこれらの問題は生じない。
【実施例１】
【００３８】
本発明の実施例１であるノイズ低減方法（撮像方法）について説明する。図２には、図１に示した撮像装置に本実施例のノイズ低減方法を適用したときの撮像光学系の構成を示す。図２において、図１中に示した構成要素については、図１と同符号を付している。本実施例では、参照光路１１０中に拡散板２０１を配置し、ガルバノミラー１０３を回動させるとともに、位相分布変更手段としての拡散板２０１を参照光路１１０の光軸に直交する方向に移動させる。これにより、物体光１２１の試料１０５に対する入射角が変化するとともに、参照光１２２の位相が変化し、スペックル分布が変化する。拡散板２０１を移動方向は参照光路１１０の光軸に直交する方向に限らず光軸に沿った方向でもよいし、また拡散版２０１を光軸に沿って回転させてもよい。
【００３９】
本実施例のノイズ低減効果をシミュレーションによって示す。まず、座標の定義として、図３に示すように、撮像光学系の光軸をｚ軸とし、光軸（ｚ軸）に直交し、かつ互いに直交する２つの軸をｘ軸およびｙ軸とする。また、物体光の試料に対する入射角の定義を、図３を用いて説明する。図３中の３０１は物体光が試料に対して入射する方向を示している。物体光の進行方向とｚ軸とがなす角度をθとし、ｘｙ平面に射影された該進行方向とｘ軸とがなす角度をφとし、θとφをともに入射角という。
【００４０】
シミュレーションでは、光源を波長λ＝５４３ｎｍのレーザとし、撮像光学系のＮＡを１．４とし、撮像光学系の結像倍率を１とした。また、物体光の入射角θを５０°で固定し、φを０°から３６０°まで３°刻みに変化させて１２０方向から照明を行ったときのホログラムデータを算出した。試料は、直径が１０μｍで屈折率が１．５９７のポリスチレンビーズを、屈折率が１．５７５のオイルに浸したものである。波面再生時に実像スペクトルを抽出する領域は、入射角によって決まる空間周波数を中心とした円内とし、その半径（数値ＮＡ）は０．７とした。
【００４１】
参照光が持つランダムな位相分布を再現するため、参照光の位相に平均値ｍが０ｒａｄで標準偏差σが０．３ｒａｄのガウシアンノイズを付加した。ガウシアンノイズとは、各点でのノイズ値ｉの出現頻度Ｐ（ｉ）が、標準偏差σで平均値ｍであるガウス分布に従うガウス乱数である。
【００４２】
まず、参照光の位相分布が撮像ごとにすべて同じである場合のシミュレーション結果として算出された３次元の屈折率分布を図４に示す。図４（ａ）はビーズの中心を通るｘｙ断面図、図４（ｂ）はビーズの下方２μｍ（ｚ＝２μｍ）を通るｘｙ断面、図４（ｃ）はビーズの中心を通るｘｚ断面図である。図４（ａ）では画像にノイズが見られるが、図４（ｂ）ではこのノイズが見られない。さらに図４（ｃ）から、焦点面近傍に集中的にノイズが生じていることが分かる。
【００４３】
本実施例の効果を示すため物体光の入射角を変化させるとともに、参照光の位相に付加するノイズを変化させた。付加したノイズは平均値ｍが０ｒａｄで標準偏差σが０．３ｒａｄのガウシアンノイズである。
【００４４】
算出された３次元の屈折率分布を図５に示す。図５（ａ）はビーズの中心を通るｘｙ断面図、図５（ｂ）はビーズの下方２μｍ（ｚ＝２μｍ）を通るｘｙ断面、図５（ｃ）はビーズの中心を通るｘｚ断面図である。図５（ａ）を図４（ａ）と比較すると、図４（ａ）で生じていたノイズが図５（ａ）では低減されている。さらに、図５（ｂ）では図４（ｂ）と同程度のノイズしか現れていない。また、図５（ｃ）から、焦点面近傍に集中的にノイズが生じる現象は起こっていないことが分かる。本シミュレーション結果から、拡散板２０１によって参照光１２２のスペックル分布を物体光１２１の入射角の変化とともに変化させることで、焦点ノイズが低減されることが分かる。
【００４５】
本実施例では拡散板２０１を用いて参照光１２２の位相を変化させる方法について説明したが、他の方法を用いてもよい。例えば、ガラス製の平行平板を参照光の光路に対して挿抜することで光路長を変化させてもよく、また参照光の光路に配置したミラーの角度や位置を変化させることで光路を変化させてもよい。
【００４６】
さらに、本実施例では拡散板２０１を参照光路１１０中に配置する方法について説明したが、拡散板を他の位置に配置してもよい。例えば、光源１０１からビームスプリッタ１０２までの間に配置したり、ビームスプリッタ１０８から撮像素子１０９の間に配置したりしてもよい。なお、拡散板は、撮像素子の瞳位置またはその共役面に配置することが望ましい。
【実施例２】
【００４７】
次に、本発明の実施例２であるノイズ低減手法（撮像方法）について説明する。図６には、図１に示した撮像装置に本実施例のノイズ低減方法を適用したときの撮像光学系の構成を示す。図６おいて、図１中に示した構成要素については、図１と同符号を付している。本実施例では、試料１０５を配置するステージを移動させる駆動装置６０１を有する。ガルバノミラー１０３を回動させるとともに駆動装置６０１によってステージを物体光１２１の光路の光軸方向（ｚ方向）に移動させる、つまりは試料１０５を移動させることで、物体光１２１の入射角を変化させるとともに試料１０５の撮像位置が変化する。この変化量をホログラムデータとともに保存しておき、波面再生時に数値的に補正する。
【００４８】
本実施例のノイズ低減効果をシミュレーションによって示す。試料の撮像位置と参照光に与えるノイズ以外の計算条件は実施例１と同じである。本実施例では、試料を撮像するｚ方向での位置を−１２μｍから１２μｍまで、物体光の入射角を変化させるとともに０．２μｍ刻みに変化させた。波面再生時にこのｚ方向位置に対応するデフォーカス収差を除くことで数値的にフォーカス位置を補正した。参照光の位相には入射角に依存しないノイズを与えた。付加したノイズは、平均値ｍが０ｒａｄで標準偏差σが０．３ｒａｄのガウシアンノイズである。
【００４９】
シミュレーション結果として算出された３次元の屈折率分布を図７に示す。図７（ａ）はビーズの中心を通るｘｙ断面図、図７（ｂ）はビーズの下方２μｍ（ｚ＝２μｍ）を通るｘｙ断面、図７（ｃ）はビーズの中心を通るｘｚ断面図である。図７（ａ）を図４（ａ）と比較すると、図４（ａ）で生じていたノイズが図７（ａ）では低減されている。さらに図７（ｂ）では、図４（ｂ）と同程度のノイズしか現れていない。さらに、図７（ｃ）から、焦点面近傍に集中的にノイズが生じる現象が起こっていないことが分かる。本シミュレーション結果から、物体光１２１の入射角の変化とともに試料１０５の位置を変化させ、後に数値的にこれを補正することで、焦点ノイズが低減されることが分かる。
【００５０】
本実施例では、試料１０５が配置されるステージを移動させて試料１０５の位置を変化させる駆動装置６０１を有する場合について説明したが、これに代えて、撮像素子１０９の位置を変化させたり撮像光学系の位置を変化させたりする駆動装置を設けてもよい。
【００５１】
また、本実施例では試料１０５をｚ方向にのみ移動させる場合について説明したが、ｚ方向ではなくｘｙ方向に移動させてもよいし、ｚ方向およびｘｙ方向に移動させてもよい。さらに、試料を傾けたり回転させたりしてもよい。
【実施例３】
【００５２】
次に、本発明の実施例３であるノイズ低減手法（撮像方法）について説明する。図８には、図１に示した撮像装置に本実施例のノイズ低減方法を適用したときの撮像光学系の構成を示す。図８おいて、図１中に示した構成要素については、図１と同符号を付している。本実施例では、対物レンズ１０６の瞳位置に収差変更手段としての空間変調器８０１を配置している。ガルバノミラー１０３を回動させるとともに空間変調器８０１の空間変調量を変化させることで、物体光の入射角を変化させるとともに撮像光学系の収差を変化させる。この収差をホログラムデータとともに保存しておき、波面再生時に数値的に補正する。
【００５３】
本実施例によるノイズ低減効果をシミュレーションによって示す。撮像光学系に与える収差と参照光に与えるノイズ以外の計算条件は実施例１と同じである。本実施例では、物体光の入射角を変化させるとともに、撮像光学系に与える非点収差の量を変化させた。非点収差は、瞳座標を（ρｘ，ρｙ）としたとき、以下の式（１）で与えられる。
【００５４】
【数１】
2015118290.tif 000003
【００５５】
ただし、Ｃは定数であり、この値を−４πから４πまでπ／１５刻みに変化させた。さらに、波面再生時に数値的にこの収差を除去し、また参照光の位相には入射角に依存しないノイズを与えた。付加したノイズは、平均値ｍが０ｒａｄで標準偏差σが０．３ｒａｄのガウシアンノイズである。
【００５６】
シミュレーション結果として算出された３次元の屈折率分布を図９に示す。図９（ａ）はビーズの中心を通るｘｙ断面図、図９（ｂ）はビーズの下方２μｍ（ｚ＝２μｍ）を通るｘｙ断面、図９（ｃ）はビーズの中心を通るｘｚ断面図である。図９（ａ）を図４（ａ）と比較すると、図４（ａ）で生じていたノイズが図９（ａ）では低減されている。さらに、図９（ｂ）では、図４（ｂ）と同程度のノイズしか現れていない。また、図９（ｃ）から、焦点面近傍に集中的にノイズが生じる現象は起こっていないことが分かる。本シミュレーション結果から、物体光の入射角の変化とともに収差を変化させ、後に数値的にこれを補正することで、焦点ノイズが低減されることが分かる。
【００５７】
本実施例では撮像光学系に非点収差を与えた場合について説明したが、球面収差やコマ収差等の他の収差を与えてもよい。
【００５８】
また、本実施例では撮像光学系に収差を与えるために空間変調器を用いたが、他の収差変更手段により収差を与えてもよい。例えば、撮像光学系に含まれるレンズの配置を変化させたり撮像ごとに瞳フィルタを入れ替えたりして収差を変化させてもよい。
【実施例４】
【００５９】
次に、本発明の実施例４であるノイズ低減手法（撮像方法）について説明する。図１０には、図１に示した撮像装置に本実施例のノイズ低減方法を適用したときの撮像光学系の構成を示す。図１０おいて、図１中に示した構成要素については、図１と同符号を付している。本実施例では、対物レンズ１０６の瞳位置に瞳フィルタ１００１を配置している。そして、この瞳フィルタが与える収差を予めコンピュータに保存しておき、波面再生時にこの収差を数値的に補正する。つまり、撮像光学系が有する既知の収差を補正する演算を行う。
【００６０】
本実施例によるノイズ低減効果をシミュレーションによって示す。撮像光学系に与える収差と参照光に与えるノイズ以外の計算条件は実施例１と同じである。撮像光学系に実施例３中の式（１）で表わされる非点収差を与え、波面再生時にこの収差を数値的に除去した。式（１）中のＣは４πとした。また、参照光の位相には入射角に依存しないノイズを与えた。付加したノイズは、平均値ｍが０ｒａｄで標準偏差σが０．３ｒａｄのガウシアンノイズである。
【００６１】
シミュレーション結果として算出された算出された３次元の屈折率分布を図１１に示す。図１１（ａ）はビーズの中心を通るｘｙ断面図、図１１（ｂ）はビーズの下方２μｍ（ｚ＝２μｍ）を通るｘｙ断面、図１１（ｃ）はビーズの中心を通るｘｚ断面図である。図１１（ａ）を図４（ａ）と比較すると、図４（ａ）で生じていたノイズが図１１（ａ）では低減されている。さらに、図１１（ｂ）では、図４（ｂ）と同程度のノイズしか現れていない。また、図１１（ｃ）から、焦点面近傍に集中的にノイズが生じる現象は起こっていないことが分かる。本シミュレーション結果から、撮像光学系に収差を与え、後に数値的にこれを補正することで、焦点ノイズが低減されることが分かる。
【００６２】
本実施例では瞳フィルタが非点収差を持つ場合について説明したが、瞳フィルタは非点収差以外の収差、例えば、球面収差やコマ収差を有してもよい。
【００６３】
また、本実施例では瞳フィルタを用いて収差を撮像光学系に与える場合について説明したが、撮像光学系に含まれるレンズを移動させて収差を与えてもよいし、対物レンズに予め収差を与えておいてもよい。また、空間変調器を用いてもよい。
【００６４】
さらに、上記各実施例ではｏｆｆ−ａｘｉｓ型の撮像光学系を備えた撮像装置について説明したが、参照光の光路長差が異なる複数の干渉縞のデータから波面を生成する位相シフト型のホログラムを用いてもよい。
【００６５】
また、上記各実施例では物体光の入射角が異なる複数回の撮像において入射角ごと（つまりは撮像ごと）に参照光の位相分布を変化させる場合について説明した。しかし、複数回の撮像のうち少なくとも１回は他の回の撮像に対して参照光の位相分布を変化させる（例えば、実施例１にて説明した拡散板を移動させたり、実施例２で説明した平行平板を挿抜したりする）ことで、ノイズの低減効果を得ることができる。
【００６６】
また、上記各実施例では試料の位置または収差を波面再生時に補正する演算を行う場合について説明したが、３次元再構成時にこれらを補正する演算を行ってもよい。
【００６７】
さらに、上記各実施例では試料の像に対して干渉縞を生成するイメージホログラムを用いる場合について説明したが、試料の透過波面の回折光の分布を取得するフーリエホログラム等を用いてもよい。
【００６８】
以上説明した各実施例は代表的な例にすぎず、本発明の実施に際しては、各実施例に対して種々の変形や変更が可能である。
【産業上の利用可能性】
【００６９】
デジタルホログラフィを用いた３次元顕微鏡装置等の撮像装置を提供することができる。
【符号の説明】
【００７０】
１０５  試料（物体）
１０９  撮像素子
１１０  参照光路
１２１  物体光
１２２  参照光
２０１  拡散板
【図１】
2015118290.tif 000004
【図２】
2015118290.tif 000005
【図３】
2015118290.tif 000006
【図４】
2015118290.tif 000007
【図５】
2015118290.tif 000008
【図６】
2015118290.tif 000009
【図７】
2015118290.tif 000010
【図８】
2015118290.tif 000011
【図９】
2015118290.tif 000012
【図１０】
2015118290.tif 000013
【図１１】
2015118290.tif 000014