JP6628103B2 - デジタルホログラフィ記録装置、デジタルホログラフィ再生装置、デジタルホログラフィ記録方法、およびデジタルホログラフィ再生方法 - Google Patents

Description

本発明はデジタルホログラフィ装置およびデジタルホログラフィ方法に関する。

以後の文章中で位相の単位はラジアンで表す。光の干渉を利用した干渉計測技術、特にデジタルホログラフィは、非接触かつ非破壊で、物体の３次元情報を得ることができるため、近年、注目を集めている測定法の一つとなっている。

デジタルホログラフィは、３次元物体への光照射によって得られる干渉縞から、コンピュータを用いて３次元物体の像を再生する技術である。一般的には例えば、３次元物体への光照射によって得られる物体光と、該物体光に対して可干渉（コヒーレント）である参照光とが作る干渉縞を、ＣＣＤ（charge coupled device）またはＣＭＯＳ（complementary metal oxide semiconductor）等の撮像素子を用いて記録する。記録された干渉縞に基づいて、コンピュータで３次元物体の像を再生する。

特許文献１には、デジタルホログラフィにより入手した干渉計信号と、試料からの蛍光信号を順次別々に測定して記録する測定および記録ステップと、干渉計信号および蛍光信号を、３次元画像を再構成するために結合するステップとを行う技術が記載されている。

特許文献２には、物体からの蛍光を光学的に検出する技術が記載されている。

特許文献３には、試料から発せられる蛍光の偏光成分を計測する技術が記載されている。

特許文献４〜７には、蛍光検出を利用する技術が記載されている。

日本国公表特許公報「特表２００４−５３８４５１号公報（２００４年１２月２４日公表）」 日本国公開特許公報「特開２００５−０３７３８８号公報（２００５年２月１０日公開）」 日本国公開特許公報「特開２００８−２２５０１３号公報（２００８年９月２５日公開）」 日本国公開特許公報「特開２００２−０５５０５０号公報（２００２年２月２０日公開）」 日本国公開特許公報「特開２００２−０７２８３５号公報（２００２年３月１２日公開）」 日本国公開特許公報「特開２０１１−１８５７５７号公報（２０１１年９月２２日公開）」 日本国公表特許公報「特表２０１３−５１１３４１号公報（２０１３年４月４日公開）」

しかしながら、従来の技術では、物体の３次元の形状情報を含む物体光によるホログラムと、物体からの蛍光とは、個別に撮像されていた。例えば、第１時刻に、レーザ光源を用いて物体から物体光を生じさせ、物体光と参照光とが干渉することにより形成されるホログラムを撮像する。次に、第２時刻に、レーザ光源をオフし、励起光源を用いて物体から蛍光を生じさせ、蛍光のみを撮像する。このように、従来の技術では、個別の光源を用いて、ホログラムと蛍光とを個別に撮像する。そのため、カメラのフレームレートでホログラムと蛍光との両方の動画像を撮像することはできない。

本発明の一態様では、ホログラムと蛍光との両方を分離再生可能な状態で同時に撮像することができる。

本発明に係るデジタルホログラフィ記録装置は、物体光が生じるように物体に物体照明光を照射する光源と、参照光と前記物体光とが干渉することにより形成されるホログラムと、蛍光の像とを撮像する撮像装置とを備え、前記物体照明光はまた、前記物体に含まれる蛍光体を励起する。

本発明に係るデジタルホログラフィ記録装置は、物体に物体照明光を照射する光源と、参照光と前記物体からの物体光とが干渉することにより形成されるホログラムと、前記物体に含まれる蛍光体が発する蛍光の像とが重畳した画像を撮像する撮像装置とを備える。

本発明に係るデジタルホログラフィ再生装置は、参照光と物体光とが干渉することにより形成されるホログラムと、インコヒーレント光の像とが重畳した画像から、空間的位相シフト法を用いて前記物体光の複素振幅を求め、前記物体光の複素振幅から前記ホログラムの強度を求め、前記画像から前記ホログラムを除去することにより、前記インコヒーレント光の像を得る。

本発明に係るデジタルホログラフィ再生装置は、参照光と物体光とが干渉することにより形成されるホログラムと、インコヒーレント光の像とが重畳した画像をフーリエ変換し、フーリエ変換された前記画像から、前記物体光の空間スペクトルを抽出し、抽出された前記物体光の空間スペクトルを逆フーリエ変換することで、前記物体光の複素振幅を求め、前記物体光の複素振幅から０次回折光成分を求め、フーリエ変換された前記画像から０次回折光成分を除去し、０次回折光成分が除去されたフーリエ変換された前記画像から、前記インコヒーレント光の像の空間スペクトルを抽出する。

本発明に係るデジタルホログラフィ再生装置は、参照光と物体光とが干渉することにより形成されるホログラムと、インコヒーレント光の像とが重畳した画像をフーリエ変換し、フーリエ変換された前記画像から、前記物体光の空間スペクトルを抽出し、抽出された前記物体光の空間スペクトルを逆フーリエ変換することで、前記物体光の複素振幅を求め、前記物体光の複素振幅から前記ホログラムの強度を求め、前記画像から前記ホログラムを除去することにより、前記インコヒーレント光の像を得る。

本発明に係るデジタルホログラフィ記録方法は、光源から出射された物体照明光を物体に照射することにより物体光を生じさせるとともに、前記物体照明光によって前記物体に含まれる蛍光体を励起するステップと、参照光と前記物体光とが干渉することにより形成されるホログラムと、前記蛍光体が発する蛍光の像とを撮像するステップとを含む。

本発明に係るデジタルホログラフィ記録方法は、物体照明光を物体に照射するステップと、参照光と前記物体からの物体光とが干渉することにより形成されるホログラムと、前記物体に含まれる蛍光体が発する蛍光の像とが重畳した画像を撮像するステップとを含む。

本発明に係るデジタルホログラフィ再生方法は、参照光と物体光とが干渉することにより形成されるホログラムと、インコヒーレント光の像とが重畳した画像から、空間的位相シフト法を用いて前記物体光の複素振幅を求めるステップと、前記物体光の複素振幅から前記ホログラムの強度を求めるステップと、前記画像から前記ホログラムを除去することにより、前記インコヒーレント光の像を得るステップとを含む。

本発明に係るデジタルホログラフィ再生方法は、参照光と物体光とが干渉することにより形成されるホログラムと、インコヒーレント光の像とが重畳した画像をフーリエ変換するステップと、フーリエ変換された前記画像から、前記物体光の空間スペクトルを抽出するステップと、抽出された前記物体光の空間スペクトルを逆フーリエ変換することで、前記物体光の複素振幅を求めるステップと、前記物体光の複素振幅から０次回折光成分を求め、フーリエ変換された前記画像から０次回折光成分を除去するステップと、０次回折光成分が除去されたフーリエ変換された前記画像から、前記インコヒーレント光の像の空間スペクトルを抽出するステップとを含む。

本発明に係るデジタルホログラフィ再生方法は、参照光と物体光とが干渉することにより形成されるホログラムと、インコヒーレント光の像とが重畳した画像をフーリエ変換するステップと、フーリエ変換された前記画像から、前記物体光の空間スペクトルを抽出するステップと、抽出された前記物体光の空間スペクトルを逆フーリエ変換することで、前記物体光の複素振幅を求めるステップと、前記物体光の複素振幅から前記ホログラムの強度を求め、前記画像から前記ホログラムを除去することにより、前記インコヒーレント光の像を得るステップとを含む。

本発明の一態様によれば、ホログラムと蛍光の像との両方を分離再生可能な状態で同時に撮像することができる。

本発明の一態様によれば、ホログラムと蛍光の像とが重畳した画像から、物体光および蛍光の像を分離して再生することができる。

本発明の一実施形態に係るデジタルホログラフィ装置の構成を示す模式図である。 本発明の一実施形態に係る再生装置が行う再生処理のフローを示す図である。 記録された画像を２次元フーリエ変換したものを概略的に示す図である。 本発明の他の実施形態に係るデジタルホログラフィ装置の構成を示す模式図である。 記録された画像を２次元フーリエ変換したものを概略的に示す図である。 本発明のさらに他の実施形態に係るデジタルホログラフィ装置の構成を示す模式図である。 本発明の一実施形態に基づくシミュレーションに用いる被写体を示す。 通常のカメラによって認識される上記被写体の画像を示す。 （ａ）は上記シミュレーションにおける記録画像を示し、（ｂ）はその一部を拡大したものを示す。 （ａ）はフーリエ変換画像を示し、（ｂ）は再生された反射光画像を示し、（ｃ）は再生された位相分布を示す。 記録画像から０次回折光成分を減算する様子を示す。 （ａ）は０次回折光成分が除去されたフーリエ変換画像を示し、（ｂ）は再生された蛍光画像を示す。 （ａ）は本発明のさらに他の実施形態に基づくシミュレーションに用いる被写体の反射光のＲチャネルの画像を示し、（ｂ）は被写体の反射光のＧチャネルの画像を示し、（ｃ）は被写体の反射光のＢチャネルの画像を示し、（ｄ）はＲＧＢが色合成された被写体の反射光画像を示す。 （ａ）は被写体の蛍光のＲチャネルの画像を示し、（ｂ）は被写体の蛍光のＧチャネルの画像を示し、（ｃ）は被写体の蛍光のＢチャネルの画像を示し、（ｄ）はＲＧＢが色合成された被写体の蛍光画像を示す。 上記被写体の高さ分布を示す。 通常のカメラによって認識される上記被写体の画像を示す。 （ａ）（ｂ）（ｃ）は、それぞれＲチャネルの記録画像より得た、反射光画像（再生像）、蛍光画像（再生像）、位相分布である。 （ａ）（ｂ）（ｃ）は、それぞれＧチャネルの記録画像より得た、反射光画像（再生像）、蛍光画像（再生像）、位相分布である。 （ａ）（ｂ）（ｃ）は、それぞれＢチャネルの記録画像より得た、反射光画像（再生像）、蛍光画像（再生像）、位相分布である。 （ａ）は、撮像装置の画素と参照光との関係を説明する模式図であり、（ｂ）は撮像装置の撮像面の一部における参照光の位相シフト量を説明する模式図である。 本発明のさらに他の実施形態に係る再生装置が行う再生処理のフローを示す図である。 上記再生処理のフローの概要を示す図である。 拡大した記録画像の一部における参照光の位相シフト量を示す図である。 （ａ）は、物体光の複素振幅を用いて再生した被写体の反射光画像（再生像）であり、（ｂ）は、物体光の複素振幅を用いて計算した位相分布である。 再生された蛍光画像を示す。

本発明の実施形態について図に基づいて以下に説明する。説明の便宜上、各項目において、上述の項目に示した部材と同一の機能を有する部材については、同一の符号を付し、適宜その説明を省略することがある。

〔実施形態１〕
本実施形態は、物体光の情報と蛍光の情報とを単一露光で同時に記録するデジタルホログラフィ装置に関する。そして、デジタルホログラフィ装置は、記録したホログラムを用いて、物体光の情報と蛍光の情報とを互いに分離して再生することができる。

（デジタルホログラフィ装置１の構成）
図１は、本実施形態のデジタルホログラフィ装置１の構成を示す模式図である。デジタルホログラフィ装置１は、軸外し型（off-axis型）のデジタルホログラフィ装置である。デジタルホログラフィ装置１は、記録装置１０（デジタルホログラフィ記録装置）と再生装置１１（デジタルホログラフィ再生装置）とを備える。記録装置１０は、撮像装置１２、波長λ１のレーザ光源ＬＳ１（第１光源）、および光学系を備える。再生装置１１は、コンピュータ等の計算機によって構成することができる。なお、レーザ光としては、可視光のみならず、不可視光（赤外線、紫外線、Ｘ線等）も利用することができる。

光学系は、ミラー等の複数の光学素子を備え、波長λ１のレーザ光（コヒーレント光）を被写体１３（物体）および撮像装置１２へ導く。具体的には、光学系は、複数の光学素子として、ビームスプリッタＢＳ１〜ＢＳ２、およびミラーＭ１〜Ｍ３、ビームエキスパンダＢＥ１〜ＢＥ２、および結像素子１５（結像光学素子）を備える。ビームエキスパンダＢＥ１〜ＢＥ２はそれぞれ、対物レンズＢＥａ、ピンホールＢＥｂ、およびコリメータレンズＢＥｃを備える。ビームスプリッタＢＳ１〜ＢＳ２は、ハーフミラーからなる。結像素子１５は、レンズを含む。ただし、これに限らず、結像させるための任意の光学素子を結像素子１５として用いることができる。

撮像装置１２は、撮像するための複数の画素がｘ方向およびｙ方向に配列した撮像面を有し、撮像面に到達した光の強度を記録する。ｘ方向はｙ方向に対して垂直である。ｚ方向はｘ方向およびｙ方向に対して垂直である。撮像装置１２は、ＣＣＤまたはＣＭＯＳ等の撮像素子を有する。撮像装置１２は、撮像面上に形成された干渉縞を記録する。この干渉縞は物体光の情報を有するホログラムである。なお、１つの画素は有限の受光領域を有する。そのため、１つの画素によって検出される光強度は、受光領域における干渉縞の光強度の積分である。撮像装置１２はカラーフィルタを備えないため、撮像装置１２の１つの画素は同時に複数の波長の光を受光することができる。すなわち撮像装置１２はモノクロ撮像装置である。撮像装置１２は、撮像した画像を示す画像データを再生装置１１に出力する。再生装置１１の詳細については後述する。

（物体光、参照光および蛍光）
レーザ光源ＬＳ１から出射された波長λ１のレーザ光は、ビームスプリッタＢＳ１によって参照光と物体照明光とに分けられる。

ビームスプリッタＢＳ１によって分けられた波長λ１の物体照明光は、ミラーＭ１を経由して、被写体１３に照射される。物体照明光が被写体１３によって散乱（反射または透過）または回折されたものが物体光である。ここでは、被写体１３による反射光を物体光として利用しているが、被写体１３の背面から物体照明光を照射し、被写体１３を透過した光を物体光として利用してもよい。物体照明光を被写体１３に照射することにより生じた物体光は、結像素子１５およびビームスプリッタＢＳ２を通過し、撮像装置１２の撮像面に入射する。

ここで、被写体１３は、蛍光体１４を含む。波長λ１の物体照明光はまた、蛍光体１４を励起する役割も果たす。波長λ１の物体照明光によって励起された蛍光体１４は、蛍光分子に応じた所定の波長の蛍光を発する。波長λ１は、蛍光の波長より短い。ただし、２光子励起または多光子励起を利用する場合、波長λ１は、蛍光の波長以上となり得る。蛍光体１４から発せられた蛍光は、物体光と同様に、結像素子１５およびビームスプリッタＢＳ２を通過し、撮像装置１２の撮像面に入射する。

ビームスプリッタＢＳ１によって分けられた波長λ１の参照光は、ミラーＭ２、ビームエキスパンダＢＥ２、ミラーＭ３を経由し、ビームスプリッタＢＳ２で反射され、撮像装置１２の撮像面に入射する。物体光および参照光が撮像面に入射する角度は、光学系に含まれる光学素子（ビームスプリッタ、ミラー等）によって調整できる。

ここでは、被写体１３の中心からの物体光は、撮像装置１２の撮像面の中心に対して、垂直に入射すると仮定している。言い換えれば、被写体１３の中心は、撮像面の中心を通り撮像面に垂直な線上に位置する。これに対し、参照光は、撮像面に斜めに入射する。すなわち、撮像面に入射する物体光の光軸と参照光の光軸との間には角度差がある。

物体光および蛍光は、結像素子１５によって、撮像装置１２の撮像面上に像を結ぶように結像される。物体光は参照光と干渉することにより、撮像面上にホログラム（結像型ホログラム）を形成する。一方、蛍光は、インコヒーレント光と呼ばれるコヒーレント性が低い光である。また、蛍光は、物体光および参照光とは波長が異なる。そのため、蛍光は、物体光および参照光とは干渉しない。撮像装置１２は、ホログラムと結像された蛍光の像とが重畳した画像を撮像することになる。なお、レーザ光の波長の光を減衰させ、蛍光の波長の光を通過させる光学フィルタ等を撮像面の手前に設けてもよい。蛍光体１４から生じる蛍光の強度が物体光の強度に比べて非常に小さい場合、蛍光の像がホログラムに埋もれないように、物体光または参照光または物体光と参照光との両方の強度を光学フィルタによって小さくしてもよい。

（再生装置１１の構成）
図２は、本実施形態の再生装置１１が行う再生処理のフローを示す図である。再生装置１１は、撮像装置１２によって記録された画像（ホログラム＋蛍光）を用いて、被写体１３の像の再生および蛍光の像の再生を行う。ただし、記録された画像は、物体光の情報（すなわち被写体１３の３次元形状情報）を含むホログラムと、蛍光とが重畳した画像であるので、そのままでは公知のデジタルホログラフィ技術で再生することができない。

ここで、記録された画像の各画素の値は、ホログラムの光強度Ｉｈおよび蛍光の光強度Ｉｆの和である。ＩｈおよびＩｆは、画素毎に異なる。蛍光は撮像面上に結像しているので、蛍光の光強度Ｉｆは、撮像装置１２から見た蛍光の分布を表す。ホログラムの光強度Ｉｈは、以下のようになる。
Ｉｈ＝｜Ｏ｜^２＋｜Ｒ｜^２＋２｜Ｏ｜｜Ｒ｜cosφ （１）
ここで、｜Ｏ｜は物体光の振幅、｜Ｒ｜は参照光の振幅、｜Ｏ｜^２は物体光の強度、｜Ｒ｜^２は参照光の強度、φは参照光に対する物体光の位相差である。｜Ｏ｜、｜Ｒ｜、およびφは、画素毎に異なり得る。なお、参照光はレーザ光を引き伸ばした平面波であるので、参照光の強度は一様な値または所定の分布（例えばガウス分布）の値であると仮定（推定）してもよい。なお、より精度よく像を再生するために、事前または事後に参照光のみの強度分布を測定してもよい。参照光の強度分布を測定することは再生の精度を上げるためのキャリブレーションに相当する。それゆえ｜Ｒ｜^２は既知の値として扱うことができる。

具体的には、再生装置１１は、記録された画像を２次元フーリエ変換する（Ｓ１）。ホログラムは、参照光と物体光とが干渉することにより形成される、様々な間隔の干渉縞を含む。このホログラムを含む画像をフーリエ変換すると、縞の間隔に応じて、物体光波（１次回折光）、０次回折光、および共役像（−１次回折光）のそれぞれのスペクトルが現れる。以後、２次元フーリエ変換した後の平面を空間周波数面と呼び、スペクトルを空間スペクトルと呼ぶ。

図３は、記録された画像を２次元フーリエ変換したものを概略的に示す図である。図３において、横軸はｘ軸の空間周波数ｆｘを示し、縦軸はｙ軸の空間周波数ｆｙを示す。ホログラムと蛍光の像とを含む画像の空間周波数面において、波長λ１の０次回折光の空間スペクトル、波長λ１の物体光波（１次回折光）の空間スペクトル、波長λ１の共役像（−１次回折光）の空間スペクトル、ならびに、蛍光の空間スペクトルが存在する。なお、共役像の空間スペクトルは物体光波の空間スペクトルに対して原点対称の位置に現れる。なお、ここでは各空間スペクトルは円形である。物体光の経路に同様の形の光透過フィルタ（開口または絞り）を配置したとき等に、このような空間スペクトルが得られる。撮像装置において記録可能な空間周波数帯域は、空間周波数面において原点を中心とした幅１／ｄの範囲である。ｄは画素間隔である。

０次回折光の空間スペクトルは原点を中心に分布する。蛍光の空間スペクトルもまた、その多くが原点を中心に分布する。それゆえ、０次回折光の空間スペクトルと蛍光の空間スペクトルとは、空間周波数面において、少なくとも一部が互いに重なって分布する。なお、記録装置１０において、被写体１３と撮像装置１２との間（例えば結像素子１５付近）に開口または絞りを配置することで、蛍光の空間スペクトルの高周波成分をカットすることができる。これにより、蛍光の空間スペクトルが主に低周波領域（原点付近）に位置するようにすることができる。

一方、λ１の物体光波の空間スペクトルは、０次回折光の空間スペクトルよりも高周波側に現れる。記録装置１０において、撮像面に入射する物体光の光軸と参照光の光軸との角度差が大きくなると、干渉縞の縞の間隔が短くなるため、物体光波の空間スペクトルが高周波側にシフトする。物体光の光軸と参照光の光軸との間にある程度の角度差を設けることにより、空間周波数面において物体光波の空間スペクトルと０次回折光の空間スペクトルとが重ならないようにすることができる。

再生装置１１は、フーリエ変換された画像から、物体光波の空間スペクトルを抽出する（Ｓ２）。例えば、再生装置１１は、フーリエ変換された画像から、あらかじめ定められた範囲を物体光波の空間スペクトルとして抽出する。物体光波の空間スペクトルが分布する範囲は、物体光の波長、物体光および参照光の入射角、および撮像装置１２から見た被写体１３の視野角（開口または絞りで制限することができる）に依存する。そのため、物体光波の空間スペクトルが分布する範囲は、あらかじめ知ることができる。

再生装置１１は、抽出した物体光波の空間スペクトルを逆フーリエ変換する。また、参照光の傾け角度に応じた位相変調量を補正する。これらにより、再生装置１１は、各画素の物体光の複素振幅（｜Ｏ｜、およびφから参照光による変調量が補正された位相）を求める（Ｓ３）。物体光の複素振幅が得られれば、再生装置１１は、物体光の複素振幅から被写体１３の像を再生することができる。

なお、上記逆フーリエ変換された結果は、記録された画像（ホログラム＋蛍光）から０次回折光、共役像、および蛍光を除いたものを表す。画素毎に参照光の位相は異なる（参照光が平面波であれば参照光の位相は周期的に変化する）。ここで、位相変調量を補正するとは、上記逆フーリエ変換された結果から得られる物体光のφに含まれる（参照光の位相に応じた）位相変調量を、参照光の位相分布に応じて補正（キャンセル）することを意味する。

再生装置１１は、物体光の複素振幅から、各画素の｜Ｏ｜^２を求める（Ｓ４）。参照光の強度｜Ｒ｜^２は上述のように既知である。再生装置１１は、物体光の複素振幅から０次回折光成分（｜Ｏ｜^２＋｜Ｒ｜^２）を求める。式（１）の右辺は全て既知となるため、再生装置１１は、ホログラムの光強度Ｉｈを特定することができる。

再生装置１１は、記録された画像の画素値から０次回折光成分（｜Ｏ｜^２＋｜Ｒ｜^２）を減算する。再生装置１１は、０次回折光成分が除去された画像をフーリエ変換する（Ｓ５）。

空間周波数面において、物体光波の空間スペクトルと蛍光の空間スペクトルとは分離されている。空間周波数面において蛍光と重畳する０次回折光成分は除去してある。そのため、再生装置１１は、０次回折光成分が除去されたフーリエ変換された前記画像から原点を中心とした所定の範囲を抽出することにより、蛍光の空間スペクトルを抽出することができる。再生装置１１は、抽出した蛍光の空間スペクトルを逆フーリエ変換することにより、蛍光の強度分布（蛍光の像）を求める（Ｓ６）。なお、上述のように、記録された画像から０次回折光成分を減算してもよいし、フーリエ変換された画像からフーリエ変換された０次回折光成分を減算してもよい。０次回折光成分の除去は、実空間およびフーリエ空間のいずれかで行うことができる。

再生装置１１は、このようにして、ホログラムと蛍光の像とが重畳した画像から、蛍光の像を分離して再生することができる。また、再生装置１１は、物体光の複素振幅から被写体１３の像を再生し、被写体１３の３次元の形状情報を得ることができる。なお、抽出した物体光波の空間スペクトルから像を再生する方法としては公知の技術を利用することができる。

本実施形態では、記録装置１０は、単一のレーザ光源ＬＳ１から出射される物体照明光によって、物体光の生成と蛍光体１４の励起とを行う。蛍光体１４を励起するためだけの光源を別途設ける必要がないので、記録装置１０の小型化を実現することができる。

また、記録装置１０は、参照光の光軸を物体光の光軸に対して傾けることにより、物体光と蛍光とを分離再生可能なように、ホログラムと蛍光の像とを同時に撮像する。すなわち、記録装置１０は、空間周波数面において、物体光波の空間スペクトルと０次回折光および蛍光の空間スペクトルとが分離されるように撮像する。このようにして記録されたホログラムから、再生装置１１は、上述のように物体光と蛍光とを分離・再生する。デジタルホログラフィ装置１は、ホログラムと蛍光とを同時に１つの画像として撮像するため、撮像装置１２のフレームレートで動画像を記録することが可能となる。また、記録装置１０の光学系を、反射型の光学系の代わりに、透過型の光学系にしてもよく、反射光の代わりに透過光、回折光、または散乱光を含む光を記録することもできる。そのため、デジタルホログラフィ装置１は、例えば、細胞内の標的の動きの解析、蛍石または塗料等の透明材料の高速特性評価、および、魚類、植物または果物等の粗面体または散乱体の高速特性評価等に利用することができる。

なお、撮像装置１２として、カラー撮像装置を用いることもできる。カラー撮像装置は、ＲＧＢ（赤緑青）の各色の画像を記録する。この場合、物体光と蛍光との波長が大きく離れていなければ、１つの色の画像に、ホログラムと蛍光の像とが重畳して含まれ得る。また、光源の波長は複数あるため、光源の波長と同じ波長の蛍光が発生し得る。それゆえ、例えば、再生装置１１は、ＲＧＢの各色の画像において、上述の物体光および蛍光の分離・再生を行う。例えば、蛍光について、ＲＧＢに対応した３つの再生画像が得られる。これにより、被写体１３に含まれる蛍光体１４が既知でない場合でも、蛍光の色を特定することができる。そのため、カラー撮像装置を用いたデジタルホログラフィ装置は、例えば、物質の反射光と蛍光の波長スペクトル同時解析に利用することができる。なお、例えば物体照明光が赤外線の場合、蛍光は物体照明光の波長より長い波長の不可視光である。この場合、例えば赤外線撮像装置を用いることができる。また例えば、カラー撮像装置のように、撮像装置の画素毎に異なる波長選択フィルタ（不可視光領域の波長を選択的に透過するフィルタ）を設け、不可視光である物体光、参照光、および蛍光を撮像してもよい。

このように、記録装置１０は、コヒーレント光によって形成されるホログラムと、被写体が発光するインコヒーレント光によって形成される被写体の像とが重畳した画像を記録（撮像）する。再生装置１１は、上記画像からインコヒーレント光によって形成される被写体の像を分離して再生する。ここでは、被写体が発光するインコヒーレント光が蛍光である場合を例に挙げて説明したが、もちろん、被写体が発光するインコヒーレント光は、蛍光に限らず、任意のインコヒーレント光（燐光、放電による発光等）であってもよい。このことは他の実施形態でも同様である。

なお、図２のＳ３のフローの後、再生装置１１は、求めた物体光の複素振幅からホログラムの光強度Ｉｈを求めてもよい。再生装置１１は、記録画像からホログラムの光強度Ｉｈを減算することにより、蛍光の像のみを抽出することができる。この場合、再生装置１１は、２回目のフーリエ変換（Ｓ５）・逆フーリエ変換（Ｓ６）の処理を省略することができる。フーリエ変換の処理は計算処理の負荷が大きい。そのため、上記のように処理することで、再生装置１１は、蛍光の像の再生処理を高速に行うことができる。

ただし、ホログラムの光強度Ｉｈを正確に算出するためには、各画素における参照光の位相を正確に特定し、再生装置１１に入力する必要がある。（１）光学系において物体光の光軸に対する参照光の光軸の角度を測定する、（２）記録画像のフーリエ変換画像において、原点から物体光波の空間スペクトルの中心までの距離および波長から参照光の光軸の上記角度を算出する、または、（３）参照光の光軸の上記角度を推定（仮定）する、ことによってユーザは参照光の位相を特定することができる。（１）においては、計測（ホログラムおよび蛍光の記録）の前または後で、キャリブレーション処理を行い角度差の情報を求めておくことも考えられる。キャリブレーション処理では、光軸の角度を測定するためにミラーまたは基準となる物体を物体光路に設置する、または光軸に沿った方向に伝搬する平面波または球面波を物体光として設定し、あらかじめホログラムを記録し、記録したホログラムをフーリエ変換し、物体光の空間スペクトルを抽出して位相情報を得る。平面波など複雑でない位相分布を持つ光であれば、参照光による位相変調量の情報抽出が容易であり、より高い精度で角度差の測定（算出）ができる。測定または推定が正確でない場合、算出されるホログラムの光強度の精度が低くなる。それゆえ、得られる蛍光の像の精度が低くなる。一方、０次回折光の強度は、参照光の位相に関係なく求めることができる。そのため、図２のフローのように処理を行うことにより、得られる蛍光の像の精度を高くすることができる。また、本実施形態の再生方法は、後述する実施形態４の再生方法よりも、蛍光の像の精度を高くすることができる。

また、記録装置１０は、物体光、参照光および蛍光が発生していない（遮光された）状態で撮像を行い、キャリブレーションのための画像として、再生装置１１に出力してもよい。この画像は、ノイズとなる迷光を撮像したものになる。再生装置１１は、記録画像の画素値（強度）から迷光の画像の画素値を減算したものを、較正された記録画像として用いてもよい。これによれば、迷光等のノイズの影響を排除し、より再生の精度を高くすることができる。

〔実施形態２〕
本実施形態では、複数の波長のレーザ光源を用いる場合について説明する。

（物体光および蛍光の記録）
図４は、本実施形態のデジタルホログラフィ装置２の構成を示す模式図である。デジタルホログラフィ装置２は、軸外し型（off-axis型）のデジタルホログラフィ装置である。デジタルホログラフィ装置２は、記録装置２０（デジタルホログラフィ記録装置）と再生装置１１とを備える。記録装置２０は、撮像装置１２、波長λ１のレーザ光源ＬＳ１、波長λ２のレーザ光源ＬＳ２（第１光源）、波長λ３のレーザ光源ＬＳ３（第２光源）、および光学系を備える。

光学系は、ミラー等の複数の光学素子を備え、波長λ１〜λ３のレーザ光（コヒーレント光）を被写体１３および撮像装置１２へ導く。具体的には、光学系は、複数の光学素子として、ビームスプリッタＢＳ１〜ＢＳ２、およびミラーＭ１〜Ｍ９、ビームエキスパンダＢＥ１〜ＢＥ４、および結像素子１５を備える。ミラーＭ５〜Ｍ９は、ダイクロイックミラーまたは偏光ビームスプリッタからなってもよい。

波長λ１〜λ３のレーザ光は、ミラーＭ４、Ｍ５およびビームスプリッタＢＳ１によって同軸に揃えられる。波長λ１〜λ３のレーザ光は、ビームスプリッタＢＳ１によって物体照明光と参照光とに分割される。

波長λ１〜λ３の物体照明光は被写体１３に照射される。被写体１３からの波長λ１〜λ３の物体光は、結像素子１５およびビームスプリッタＢＳ２を通過して、撮像装置１２の撮像面に入射する。

波長λ２、λ３の物体照明光はまた、蛍光体１４を励起する役割（励起光の役割）も果たす。加えて、蛍光体１４は、波長λ１の物体照明光によっても励起されてもよい。蛍光体１４は複数の物体照明光によって励起されるので、発する蛍光の強度が増強される。蛍光体１４から発せられた蛍光は、物体光と同様に、結像素子１５およびビームスプリッタＢＳ２を通過し、撮像装置１２の撮像面に入射する。

波長λ１の参照光は、ミラーＭ６、Ｍ７を通過し、ミラーＭ２、ビームエキスパンダＢＥ２およびミラーＭ３を経由し、ミラーＭ９、Ｍ８を通過し、ビームスプリッタＢＳ２で反射され、撮像装置１２の撮像面に入射する。

波長λ２の参照光は、ミラーＭ６を通過し、ミラーＭ７で反射され、ビームエキスパンダＢＥ３を経由し、ミラーＭ９で反射され、ミラーＭ８を通過し、ビームスプリッタＢＳ２で反射され、撮像装置１２の撮像面に入射する。

波長λ３の参照光は、ミラーＭ６で反射され、ビームエキスパンダＢＥ４を経由し、ミラーＭ８で反射され、ビームスプリッタＢＳ２で反射され、撮像装置１２の撮像面に入射する。

記録装置２０は、波長毎に参照光の経路を分けることで、波長毎に異なる入射角（ｘ−ｚ平面における入射角およびｙ−ｚ平面における入射角）で参照光を撮像面に入射させる。

撮像装置１２はモノクロ撮像装置である。撮像装置１２は、複数の波長に対応する複数のホログラムと結像された蛍光の像とが重畳した画像を撮像することになる。

（物体光および蛍光の再生）
図５は、記録された画像を２次元フーリエ変換したものを概略的に示す図である。図５において、横軸はｘ軸の空間周波数ｆｘを示し、縦軸はｙ軸の空間周波数ｆｙを示す。

再生装置１１は、記録された画像を２次元フーリエ変換することにより、フーリエ変換された画像（図５）を得る。ホログラムと蛍光の像とを含む画像の空間周波数面において、各波長λ１〜λ３の０次回折光の空間スペクトル、各波長λ１〜λ３の物体光波（１次回折光）の空間スペクトル、各波長λ１〜λ３の共役像（−１次回折光）の空間スペクトル、ならびに、蛍光の空間スペクトルが存在する。

各波長λ１〜λ３の０次回折光の空間スペクトルは原点を中心に分布し、互いに重なっている。蛍光の空間スペクトルもまた、その多くが原点を中心に分布する。それゆえ、０次回折光の空間スペクトルと蛍光の空間スペクトルとは、空間周波数面において、少なくとも一部が互いに重なって分布する。

一方、λ１〜λ３の物体光波の空間スペクトルは、各波長の参照光の入射角（物体光の光軸と参照光の光軸との間の角度）に応じて、互いに異なる位置に現れる。例えば撮像面に入射する参照光がｘ軸方向に傾いていれば、ｆｘの高周波側に物体光波の空間スペクトルが現れる。同様に撮像面に入射する参照光がｙ軸方向に傾いていれば、ｆｙの高周波側に物体光波の空間スペクトルが現れる。各波長の物体光波の空間スペクトルが分離可能なように（重ならないように）、各波長の入射角が異なるように記録装置２０は構成される。

再生装置１１は、フーリエ変換された画像から、波長毎に物体光波の空間スペクトルを抽出する。それゆえ、再生装置１１は、波長毎に、抽出した物体光波の空間スペクトルから各画素の物体光の複素振幅を求めることができる。

再生装置１１は、各波長の０次回折光成分を求め、フーリエ変換された画像から各波長の０次回折光成分を除去することにより、蛍光の空間スペクトルを抽出することができる。このようにして、複数の波長のホログラムと蛍光の像とが重畳している場合でも、再生装置１１は、各波長の物体光および蛍光を分離し、再生することができる。デジタルホログラフィ装置２は、被写体１３の３次元形状のカラーイメージングを行うと共に、蛍光の画像（モノクロ）を得ることができる。また、各波長の参照光の入射角が同じ場合でも、波長差に応じて空間周波数が異なることを利用して各波長の物体光波の空間スペクトルを分離しても良い。また、物体光と蛍光、または各波長の物体光同士の空間スペクトルが一部重畳する場合、重畳する部分を除いて空間スペクトルを抽出し像を再生しても良い。

〔実施形態３〕
本実施形態では、複数の波長のレーザ光源およびカラー撮像装置を用いる場合について説明する。

（物体光および蛍光の記録）
図６は、本実施形態のデジタルホログラフィ装置３の構成を示す模式図である。デジタルホログラフィ装置３は、軸外し型（off-axis型）のデジタルホログラフィ装置である。デジタルホログラフィ装置３は、記録装置３０（デジタルホログラフィ記録装置）と再生装置１１とを備える。記録装置３０は、撮像装置３２、波長λ１のレーザ光源ＬＳ１、波長λ２のレーザ光源ＬＳ２、波長λ３のレーザ光源ＬＳ３、および光学系を備える。

光学系は、ミラー等の複数の光学素子を備え、波長λ１〜λ３のレーザ光（コヒーレント光）を被写体１３および撮像装置３２へ導く。具体的には、光学系は、複数の光学素子として、ビームスプリッタＢＳ１〜ＢＳ２、およびミラーＭ１〜Ｍ５、ビームエキスパンダＢＥ１〜ＢＥ２、および結像素子１５を備える。

物体照明光、物体光および蛍光の経路は、実施形態２と同様である。本実施形態では、実施形態２とは異なり、波長λ１〜λ３の参照光は、同じ経路を通り、互いに同軸で撮像面に入射する。ただし、物体光の光軸と参照光の光軸との間には角度差がある。

本実施形態では、実施形態２とは異なり、撮像装置３２はカラー撮像装置である。撮像装置３２は、Ｒ（赤）のサブ画素、Ｇ（緑）のサブ画素、およびＢ（青）のサブ画素毎に撮像を行い、カラー画像を生成する。カラー画像には、Ｒチャネルの情報、Ｇチャネルの情報、Ｂチャネルの情報が含まれており、カラー画像から各色の画像を分離することができる。ここで、波長λ１の光はＲチャネルにのみ記録され、波長λ２の光はＧチャネルにのみ記録され、波長λ３の光はＢチャネルにのみ記録されるとする。蛍光はその波長に応じて１つまたは複数のチャネルに記録される。

例えば、Ｒチャネルの画像は、波長λ１の波長のホログラムと蛍光の像とが重畳した画像である。例えば、Ｇチャネルの画像は、波長λ２の波長のホログラムと蛍光の像とが重畳した画像である。例えば、Ｂチャネルの画像は、波長λ３の波長のホログラムのみを含む画像である。

再生装置１１は、Ｒチャネルの画像から、上記実施形態と同様に、波長λ１の物体光と蛍光とを分離し、再生することができる。Ｇチャネルの画像についても同様である。再生装置１１は、Ｂチャネルの画像から、波長λ３の物体光を再生することができ、結果的に、蛍光の画像としては真っ暗な画像を得る。

再生装置１１は、物体光または蛍光について、ＲＧＢの再生像を合成することで、カラーの再生像を得ることができる。デジタルホログラフィ装置３は、被写体１３の３次元形状のカラーイメージングを行うと共に、蛍光のカラー画像を得ることができる。もちろん、ＲチャネルおよびＧチャネルに加えて、Ｂチャネルの画像が、波長λ３の波長のホログラムと蛍光の像とが重畳した画像である場合も、同様にして蛍光のカラー画像を得ることができる。

〔実施形態４〕
本実施形態では、フーリエ変換を用いない再生方法について説明する。本実施形態の再生装置１１は、参照光の位相が空間的に異なる（撮像面の画素毎に参照光の位相が異なる）場合に、空間的位相シフト法を用いて物体光の複素振幅を求める。なお、記録装置は、実施形態１と同様の構成とすることができる。以下に具体的に説明する。

（記録装置）
図２０の（ａ）は、撮像装置の画素と参照光との関係を説明する模式図であり、図２０の（ｂ）は撮像装置の撮像面（または記録画像）の一部における参照光の位相シフト量を説明する模式図である。図２０の（ｂ）の１つの矩形は１つの画素に対応する。

撮像装置の撮像面には、複数の画素１２ａが配列している。例えば、参照光がｚ軸に対してｘ軸方向に傾いている場合、ｘ軸方向において、画素１２ａ毎に参照光の位相が異なる。例えば、ｘ−ｚ平面における入射角θが４ｄ・sinθ＝λを満たす場合、ある画素１２ａにおける参照光の位相を基準（０）とすると、ｘ軸方向に沿って、参照光の位相が０、π／２、π、３π／２である画素が並ぶ。ｄは画素ピッチである。この基準からの位相のずれを位相シフト量と呼ぶ。記録画像のホログラムにおける参照光の位相シフト量も同様に、図２０の（ｂ）に示すようになる。

なお、参照光の光軸をｘ軸およびｙ軸の両方向に傾けると、参照光の位相は斜めにずれていく。位相シフト量の差はπ／２に限らない。また、物体光と参照光の光軸が揃ったin-line型の光学系において、空間光変調素子アレイ、ガラスまたは波長板等の位相をシフトさせる素子のアレイ、タルボット効果を起こす回折素子、あるいは、偏光子アレイ等を利用して図２０の（ｃ）に示すような参照光の位相分布の記録画像を得てもよい。また、参照光の代わりに、物体光の光軸を傾けるか、物体光路に位相シフト素子アレイを設けるかによって、物体光に位相ずれを生じさせてもよい。

（再生装置）
図２１は、本実施形態の再生装置１１が行う再生処理のフローを示す図である。図２２は、再生処理のフローの概要を示す図である。図２２における矩形は一部を拡大して示される画像の画素を表し、矩形の中の数値は位相シフト量を示す。ここでは、ｘ軸方向（水平方向）の隣接画素において、位相シフト量の差が２π／３である場合を例にして説明する。記録画像には、物体光および位相がシフトした参照光が干渉したホログラムと、蛍光の像とが重畳して記録されている。

再生装置１１は、記録画像から、位相シフト量が同じ画素を抽出し、位相シフト量毎に、複数の画像を得る。再生装置１１は、複数の画像について、画素の間の画素値を補間し、補間された複数の画像を得る（Ｓ１１）。

再生装置１１は、空間的位相シフト法を用いて、互いに位相シフト量が異なる補間された複数の画像から物体光の複素振幅の分布を得る（Ｓ１２）。この際、ホログラムのみから物体光の複素振幅を求める空間的位相シフト法の公知の式を利用することができる。再生装置１１は、物体光の複素振幅分布から、被写体の再生像を得ることができる。

再生装置１１は、物体光の複素振幅分布と既知の参照光の強度分布と参照光の位相シフト量とから、ホログラムの光強度Ｉｈの分布を求める（Ｓ１３）。再生装置１１は、記録画像からホログラムの光強度分布を減算（ホログラムを除去）することにより、蛍光の像を得る（Ｓ１４）。このようにして、再生装置１１は、ホログラムと蛍光の像とが重畳した記録画像から、物体光と蛍光の像とを再生することができる。

フーリエ変換の処理は計算処理の負荷が大きい。本実施形態によれば、再生処理においてフーリエ変換を必要としないので、高速に物体光および蛍光の像の分離・再生を行うことができる。そのため、本実施形態の再生装置１１は、物体光および蛍光の像のリアルタイムでの動画像再生を可能にする。

（変形例）
上記では、画素を抽出・補間する場合を例に挙げて説明した。さらなる再生処理負荷の低減を目的として、補間処理を行わずに物体光と蛍光の像とを分離・再生する例について説明する。

図２３は、拡大した記録画像の一部における参照光の位相シフト量を示す図である。参照光の光軸をｘ軸およびｙ軸の両方向に傾けると、参照光の位相シフト量の分布は、図２３に示すような分布になる。

記録画像の画素値Ｉ（ｘ、ｙ）は、蛍光の像の光強度Ｉｆ（ｘ、ｙ）およびホログラムの光強度Ｉｈ（ｘ、ｙ）の和である。ｘ、ｙは画素の座標を示す。
Ｉ(x,y)＝Ｉｆ(x,y)＋Ｉｈ(x,y) （２）
Ｉｈ(x,y)＝｜Ｕｏ(x,y)｜^２＋｜Ｕｒ(x,y)｜^２
＋２Ａｏ(x,y)Ａｒ(x,y)cos｛φｏ(x,y)−φｒ(x,y)｝ （３）
ここで、物体光の複素振幅をＵｏ(x,y)＝Ａｏ(x,y)exp｛ｊφｏ(x,y)｝とし、参照光の複素振幅をＵｒ(x,y)＝Ａｒ(x,y)exp｛ｊφｒ(x,y)｝とした。Ａｏ、Ａｒは振幅を表し、φｏ、φｒは位相を表す。ｊは虚数単位である。Ｍをｘ軸方向の画素数、Ｎをｙ軸方向の画素数とし、０≦ｘ≦Ｍ−１、０≦ｙ≦Ｎ−１とする（ｘ、ｙは画素番地（整数））。

再生装置１１は、対象画素（ｘ、ｙ）の近傍にある画素の値を利用して、空間的位相シフト法によって物体光の複素振幅を算出する。具体的には以下の式にて対象画素の物体光を算出する。
・（ｘ＋ｙ）÷４の余りが０のとき、
Ｕｏ(x,y)＝
[{Ｉ(x,y)−Ｉ(x+1,y+1)}＋ｊ{Ｉ(x+1,y)＋Ｉ(x,y+1)−Ｉ(x,y)−Ｉ(x+1,y+1)}]
／{４Ａｒ(x,y)} （４）
・（ｘ＋ｙ）÷４の余りが１のとき、
Ｕｏ(x,y)＝
[−{Ｉ(x+1,y)＋Ｉ(x,y+1)−Ｉ(x,y)−Ｉ(x+1,y+1)}＋ｊ{Ｉ(x,y)−Ｉ(x+1,y+1)}]
／{４Ａｒ(x,y)} （５）
・（ｘ＋ｙ）÷４の余りが２のとき、
Ｕｏ(x,y)＝
[−{Ｉ(x,y)−Ｉ(x+1,y+1)}−ｊ{Ｉ(x+1,y)＋Ｉ(x,y+1)−Ｉ(x,y)−Ｉ(x+1,y+1)}]
／{４Ａｒ(x,y)} （６）
・（ｘ＋ｙ）÷４の余りが３のとき、
Ｕｏ(x,y)＝
[{Ｉ(x+1,y)＋Ｉ(x,y+1)−Ｉ(x,y)−Ｉ(x+1,y+1)}−ｊ{Ｉ(x,y)−Ｉ(x+1,y+1)}]
／{４Ａｒ(x,y)} （７）
となる。これによって、再生装置１１は、物体光の複素振幅を求めることができる。再生装置１１は、物体光の複素振幅分布と既知の参照光の強度分布と参照光の位相シフト量（参照光の位相）とから、ホログラムの光強度Ｉｈの分布を求める。式（２）より、再生装置１１は、記録画像の画素値Ｉからホログラムの光強度Ｉｈを減算することにより、蛍光の像の光強度Ｉｆを求める。また、参照光の強度分布にむらがある場合、記録画像から参照光の強度分布を減算した後に空間的位相シフト法の計算をしてもよい。これにより、参照光の強度のむらによる出力される像の画質劣化の問題を緩和することができる。

この変形例によれば、フーリエ変換および画素の補間処理を行わずに済むので、より再生処理の負荷を低減することができ、より高速に再生を行うことができる。

〔シミュレーション１〕
本発明の実施形態１に基づくホログラムの記録および再生のシミュレーション結果について説明する。ここでは、図１に示すデジタルホログラフィ装置１を用いると仮定する。

図７は、シミュレーションに用いる被写体を示す。図７の（ａ）は、被写体の反射光画像（波長λ１の物体光の画像）を示し、図７の（ｂ）は、被写体の蛍光の画像を示し、図７の（ｃ）は、被写体の高さ分布を示す。例えば、「Ｃ」の文字の箇所からは、波長λ１の物体光はほとんど反射されないが、蛍光は強く出射される。通常のカメラまたは人の目で被写体を観察すると、図８に示す画像が認識され、反射光（物体光）と蛍光とを分離して識別することができない。なお、図７の（ｃ）において、被写体の高さ（ｚ方向における撮像装置側への高さ）は、明暗で表されている。明るい箇所は暗い箇所よりも高い。

シミュレーションの条件は以下の通りである。光源の波長はλ＝５３２ｎｍである。撮像素子の画素数は２０４８×２０４８（横×縦）、画素ピッチは横方向２．２μｍ、縦方向２．２μｍである。ｏｆｆ−ａｘｉｓ型のホログラム記録配置を仮定し、参照光の入射角度はｚ軸からｘ軸に（ｘ−ｚ平面で）、ｙ軸に（ｙ−ｚ平面で）それぞれ６．０５×１０^−２ｒａｄである。物体照明光と参照光の強度比は１：１である。参照光強度分布画像は、計測の事前または事後に記録する。反射光および蛍光の強度は０〜２５５（階調数は２５６）である。物体の最大高さは２００ｎｍである。単一の光源により反射光と蛍光を得るため、本シミュレーションでは以下の条件を設ける。反射光と蛍光のエネルギー総和は物体照明光のエネルギーよりも少ないものとする。単一光源より射出された光は励起光の役割と反射光およびそのホログラムを生成する役割とを兼ねる。また、物体光と蛍光とを分離するための、物体光または蛍光を遮光するフィルタを設けない。そのため、フィルタによる蛍光の光量損失がなく明るい画像記録が期待できる。

上記の条件の下、計算機によって、図９の（ａ）に示す記録画像を記録し、再生装置１１により再生像を計算するシミュレーションを行った。なお、ホログラムの記録も計算機によるシミュレーションによって行っている。記録画像は、ホログラムと蛍光の像とが重畳した画像である。図９の（ｂ）に記録画像の一部を拡大したものを示す。このように、記録画像には、干渉縞（ホログラム）が記録されている。

図９の（ａ）の記録画像をフーリエ変換すると、図１０の（ａ）に示す画像（フーリエ変換画像）が得られる。反射光の成分は、ホログラムとして記録されているため、蛍光画像の成分とは分離して取得できる。フーリエ変換画像から反射光成分の空間スペクトル（右下の丸囲みの範囲）を抽出したものを逆フーリエ変換する。逆フーリエ変換された画像から、図１０の（ｂ）に示す反射光画像（再生像）と、図１０の（ｃ）に示す位相分布とが求められる。位相分布は、被写体の高さを用いた光の位相で表したものである。位相分布には３次元形状の情報が含まれるため、図１０の（ｃ）より被写体の高さ情報が分かる。

次に、図１１に示すように、記録画像（図９の（ａ））から、０次回折光成分を減算し、フーリエ変換すると、図１２の（ａ）の画像を得る。このとき、０次回折光成分を減算せずに得られたフーリエ変換画像（図１０の（ａ））と異なり、画像中心部にある空間スペクトルの分布が変わる。これは、０次回折光成分が除去されたことにより、蛍光画像の空間スペクトルのみが中心部に残っていることを示す。すなわち、図１２の（ａ）のフーリエ変換画像から、蛍光画像の空間スペクトル（中央の丸囲みの範囲）のみを抽出し、抽出した空間スペクトルを逆フーリエ変換することにより、図１２の（ｂ）の蛍光画像（再生像）を得た。

定量的評価として、元の画像（図７）に対する再生像の平均二乗誤差、相互相関係数、および信号対雑音比を計算した。反射光画像（再生像）の平均二乗誤差は８．６×１０^−２、相互相関係数は１．００、信号対雑音比は５６ｄＢであった。蛍光画像（再生像）の平均二乗誤差は８．５×１０^−２、相互相関係数は１．００、信号対雑音比は４８ｄＢであった。高さ分布（位相分布）の平均二乗誤差は２．７×１０^−２ｎｍ、相互相関係数は１．００、信号対雑音比は７９ｄＢであった。以上より、動的な物体の瞬間の反射光画像、蛍光画像、３次元形状の情報を単一光源の記録装置により高精度に得られる。

〔シミュレーション２〕
本発明の実施形態３に基づくホログラムの記録および再生のシミュレーション結果について説明する。ここでは、図６に示すデジタルホログラフィ装置３を使用すると仮定する。

図１３は、シミュレーションに用いる被写体の反射光画像を示す。図１３の（ａ）は、被写体の反射光のＲ（赤）チャネルの画像を示し、図１３の（ｂ）は、被写体の反射光のＧ（緑）チャネルの画像を示し、図１３の（ｃ）は、被写体の反射光のＢ（青）チャネルの画像を示し、図１３の（ｄ）は、ＲＧＢが色合成された被写体の反射光画像を示す。図１３の（ｄ）では各色が単色（グレー）で表示されているため、色合成された画像のＲＢの成分が薄く見えるが、図１３の（ｄ）は、各色の画像が合成されたものである。なお、他の色合成された画像でも同様である。

図１４は、シミュレーションに用いる被写体の蛍光の画像を示す。図１４の（ａ）は、被写体の蛍光のＲチャネルの画像を示し、図１４の（ｂ）は、被写体の蛍光のＧチャネルの画像を示し、図１４の（ｃ）は、被写体の蛍光のＢチャネルの画像を示し、図１４の（ｄ）は、ＲＧＢが色合成された被写体の蛍光画像を示す。

図１５は、被写体の高さ分布を示す。図１５において被写体の高さは、明暗で表されている。通常のカメラや人の目で被写体を観察すると、図１６に示す画像が認識され、反射光（物体光）と蛍光とを分離して識別することができない。

シミュレーションの条件は以下の通りである。光源の波長はλ１＝６４０ｎｍ、λ２＝５３２ｎｍ、λ３＝４７３ｎｍである。Ｂａｙｅｒ型のカラーフィルタアレイが付されたカラーイメージセンサを撮像装置３２として用いる。カラーフィルタアレイの各チャネルはλ１、λ２、λ３の３波長の光のうち２つの波長の光を遮光する。撮像素子の画素数は２０４８×２０４８（横×縦）、画素ピッチは横方向１．４μｍ、縦方向１．４μｍである。ｏｆｆ−ａｘｉｓ型のホログラム記録配置を仮定し、参照光の入射角度はｘ軸、ｙ軸それぞれ６．０５×１０^−２ｒａｄである。物体照明光と参照光の強度比は１：１である。各波長の参照光の強度分布画像は、計測の事前または事後に各波長毎に独立して記録する。反射光および蛍光の強度は０から２５５（階調数２５６）である。被写体の最大高さは２４０ｎｍである。各光源による光により反射光を得られ、３つの光源のうち２つの光源の光により蛍光体が励起されると仮定する。本シミュレーションでは各光源において、反射光と蛍光のエネルギー総和は物体照明光のエネルギーよりも少ないものとしている。また、反射光と蛍光で等しい波長を有する場合を想定している。例えば、λ３により発生した蛍光と、λ１による反射光は、等しい波長成分を含む条件を想定している。すなわち、反射光と蛍光を波長フィルタにより分離できない条件を想定している。なお、波長λ３の物体照明光により、シアンの蛍光（文字「Ａ」）、黄緑色の蛍光（文字「０」）、黄色の蛍光（文字「Ｂ」）、および赤色の蛍光（文字「１」「２」）が生じるとする。波長λ２の物体照明光により、黄緑色の蛍光、黄色の蛍光、および赤色の蛍光が生じるとする。波長λ１の物体照明光は、蛍光体を励起しないとする。

上記の条件の下、計算機によって、ＲＧＢの記録画像をカラーイメージセンサの各チャネルで記録し、再生像を計算するシミュレーションを行った。なお、ホログラムの記録も計算機によるシミュレーションによって行っている。各記録画像は、ホログラムと蛍光の像とが重畳した画像である。

カラーイメージセンサのＲＧＢ各チャネルで記録された各色の記録画像から、色毎に、シミュレーション１と同様の手続きにより反射光画像、蛍光画像、高さ分布（位相分布）を再生する。Ｒチャネルの記録画像より図１７に示す画像を得た。Ｇチャネルの記録画像より図１８に示す画像を得た。Ｂチャネルの記録画像より図１９に示す画像を得た。

図１７の（ａ）は、Ｒチャネルの記録画像より得た反射光画像（再生像）であり、図１７の（ｂ）は、Ｒチャネルの記録画像より得た蛍光画像（再生像）であり、図１７の（ｃ）は、Ｒチャネルの記録画像より得た位相分布である。図１８の（ａ）（ｂ）（ｃ）は、それぞれＧチャネルの記録画像より得た、反射光画像（再生像）、蛍光画像（再生像）、位相分布である。図１９の（ａ）（ｂ）（ｃ）は、それぞれＢチャネルの記録画像より得た、反射光画像（再生像）、蛍光画像（再生像）、位相分布である。図１８の（ｃ）および図１９の（ｃ）の位相分布では、位相の折り畳み（位相が２πを超える位置で明から暗に変化する）が生じているため、中央が暗く表示されている。これらの位相分布からでも、位相接続を行うことで３次元形状の情報を得ることができる。

各チャネルにおいて、反射光と蛍光とを分離・再生し、高さ分布（３次元形状）の情報も得られることが分かる。定量的評価として、元の画像（図１３〜１５）に対する各チャネルにおける再生像の平均二乗誤差、相互相関係数、および信号対雑音比を計算した。Ｒチャネルの反射光画像の平均二乗誤差は５．８×１０^−１、相互相関係数は１．００、信号対雑音比は３０ｄＢであった。Ｒチャネルの蛍光画像の平均二乗誤差は９．６×１０^−１、相互相関係数は１．００、信号対雑音比は２９ｄＢであった。Ｇチャネルの反射光画像の平均二乗誤差は６．０×１０^−１、相互相関係数は１．００、信号対雑音比は３１ｄＢであった。Ｇチャネルの蛍光画像の平均二乗誤差は１．０、相互相関係数は１．００、信号対雑音比は３０ｄＢであった。Ｂチャネルの反射光画像の平均二乗誤差は６．９×１０^−１、相互相関係数は１．００、信号対雑音比は３０ｄＢであった。Ｂチャネルの蛍光画像の平均二乗誤差は５．９×１０^−１、相互相関係数は１．００、信号対雑音比は２９ｄＢであった。高さ分布（位相分布）の平均二乗誤差は３．４×１０^−１ｎｍ、相互相関係数は１．００、信号対雑音比は５８ｄＢであった。以上より、動的な物体の瞬間のカラー反射光画像、カラー蛍光画像、３次元形状の情報を得られる。

〔シミュレーション３〕
本発明の実施形態４に基づくホログラムの記録および再生のシミュレーション結果について説明する。ここでは、図１に示すデジタルホログラフィ装置１を使用すると仮定する。

シミュレーションに用いる被写体（反射光画像、蛍光画像、および高さ分布）は、シミュレーション１（図７）のものと同じである。また、記録におけるシミュレーションの条件もシミュレーション１と同じである。すなわち、本シミュレーションに用いる記録画像はシミュレーション１（図９）のものと同じである。

上記の条件の下、計算機によって、図９の（ａ）に示す記録画像を記録し、再生装置１１により再生像を計算するシミュレーションを行った。なお、ホログラムの記録も計算機によるシミュレーションによって行っている。記録画像は、ホログラムと蛍光の像とが重畳した画像である。

再生装置１１により、空間的位相シフト法（上式（４）〜（７））を用いて、記録画像から物体光の複素振幅（の分布）を求めた。図２４の（ａ）は、物体光の複素振幅を用いて再生した被写体の反射光画像（再生像）であり、図２４の（ｂ）は、物体光の複素振幅を用いて計算した位相分布である。反射光画像および位相分布の成分（物体光の情報）はホログラムとして記録されており、空間的位相シフト法の信号処理を施すと蛍光画像から分離してこれらを得ることができる。位相分布には３次元形状の情報が含まれるため、図２４の（ｂ）より被写体の高さ情報が分かる。

再生装置１１により、物体光の複素振幅、参照光の強度分布および参照光の位相から、ホログラムの画像（ホログラムの光強度）を求めた。そして、記録画像からホログラムの画像を減算することにより、図２５に示す蛍光画像（再生）を得た。これは、ホログラムと蛍光の像とが重畳した記録画像から、ホログラムの成分を除去する処理に対応する。このように、再生装置１１によって、ホログラムと蛍光の像とが重畳した記録画像から、反射光画像（物体光の再生像）と蛍光画像とを分離して再生することができた。

定量的評価として、元の画像（図７）に対する再生像の平均二乗誤差、相互相関係数、および信号対雑音比を計算した。反射光画像の平均二乗誤差は５．８×１０^−１、相互相関係数は０．９９９、信号対雑音比は３９ｄＢであった。蛍光画像の平均二乗誤差は７．５×１０^−１、相互相関係数は０．９９９、信号対雑音比は２９ｄＢであった。高さ分布（位相分布）の平均二乗誤差は３．４×１０^−１ｎｍ、相互相関係数は１．００、信号対雑音比は５７ｄＢであった。以上より、動的な物体の瞬間の反射光画像、蛍光画像、３次元形状の情報を単一光源の記録装置により高精度に得られる。

〔ソフトウェアによる実現例〕
再生装置１１の処理は、集積回路（ＩＣチップ）等に形成された論理回路（ハードウェア）によって実現してもよいし、ＣＰＵ（Central Processing Unit）を用いてソフトウェアによって実現してもよい。

後者の場合、再生装置１１は、各機能を実現するソフトウェアであるプログラムの命令を実行するＣＰＵ、上記プログラムおよび各種データがコンピュータ（またはＣＰＵ）で読み取り可能に記録されたＲＯＭ（Read Only Memory）または記憶装置（これらを「記録媒体」と称する）、上記プログラムを展開するＲＡＭ（Random Access Memory）などを備えている。そして、コンピュータ（またはＣＰＵ）が上記プログラムを上記記録媒体から読み取って実行することにより、本発明の目的が達成される。上記記録媒体としては、「一時的でない有形の媒体」、例えば、テープ、ディスク、カード、半導体メモリ、プログラマブルな論理回路などを用いることができる。また、上記プログラムは、該プログラムを伝送可能な任意の伝送媒体（通信ネットワークや放送波等）を介して上記コンピュータに供給されてもよい。なお、本発明は、上記プログラムが電子的な伝送によって具現化された、搬送波に埋め込まれたデータ信号の形態でも実現され得る。

〔まとめ〕
本発明の一態様に係るデジタルホログラフィ記録装置は、物体光が生じるように物体に物体照明光を照射する光源と、参照光と前記物体光とが干渉することにより形成されるホログラムと、蛍光の像とを撮像する撮像装置とを備え、前記物体照明光はまた、前記物体に含まれる蛍光体を励起する。

上記の構成によれば、ホログラムと蛍光の像との両方を分離再生可能な状態で同時に撮像することができる。また、１つの光源から出射される物体照明光が、物体光を生じさせる役割と、蛍光体を励起する励起光としての役割とを兼ねる。そのため、上記デジタルホログラフィ記録装置は小型化に有利である。

前記光源は、単一波長の前記物体照明光を出射する構成であってもよい。

前記撮像装置は、前記物体光と参照光とが干渉することにより形成されるホログラムと、前記蛍光体が発する蛍光の像とが重畳した画像を撮像する構成であってもよい。

前記蛍光を前記撮像装置の撮像面に結像させる結像光学素子を備える構成であってもよい。

上記の構成によれば、記録された画像から結像した蛍光の像を再生することができる。

前記撮像装置は、カラーフィルタを用いないモノクロ撮像装置であってもよい。

上記の構成によれば、カラーフィルタを用いずとも、物体光と蛍光とを分離・再生することができる。また、撮像装置は、カラーフィルタによる蛍光の光量損失がなく、明るい画像を記録することができる。

前記光源は、前記物体に第１波長の物体照明光を照射する第１光源であり、前記第１波長とは異なる第２波長の物体光が生じるように前記物体に前記第２波長の物体照明光を照射する第２光源を備え、前記撮像装置は、前記第１波長のホログラムと、前記第２波長のホログラムと、前記蛍光の像とを撮像し、前記第１波長の物体照明光および前記第２波長の物体照明光の両方は、前記蛍光体を励起する構成であってもよい。

上記の構成によれば、２つの波長の物体光による複数のホログラムを得ることができ、物体の分光計測を行うことができる。また、２つの波長の物体照明光によって蛍光体を励起することができるため、蛍光体が発する蛍光の強度を増強することができる。

前記撮像装置に入射する前記物体光の光軸に対して、前記撮像装置に入射する前記参照光の光軸は傾いている構成であってもよい。

上記の構成によれば、記録画像の空間周波数面において、精度よく物体光波の空間スペクトルを抽出することができる。そのため、物体光の再生の精度を高くすることができる。

本発明の一態様に係るデジタルホログラフィ記録装置は、物体に物体照明光を照射する光源と、参照光と前記物体からの物体光とが干渉することにより形成されるホログラムと、前記物体に含まれる蛍光体が発する蛍光の像とが重畳した画像を撮像する撮像装置とを備える。

上記の構成によれば、ホログラムと蛍光の像との両方を分離再生可能な状態で同時に撮像することができる。

本発明の一態様に係るデジタルホログラフィ再生装置は、参照光と物体光とが干渉することにより形成されるホログラムと、蛍光の像とが重畳した画像から、空間的位相シフト法を用いて前記物体光の複素振幅を求め、前記物体光の複素振幅から前記ホログラムの強度（光強度分布）を求め、前記画像から前記ホログラムを除去する（ホログラムのの光強度分布を減算する）ことにより、前記蛍光の像（を示す画像）を得る。

上記の構成によれば、ホログラムと蛍光の像とが重畳した画像から、高速に物体光の情報および蛍光の情報を分離することができる。

本発明の一態様に係るデジタルホログラフィ再生装置は、参照光と物体光とが干渉することにより形成されるホログラムと、蛍光の像とが重畳した画像をフーリエ変換し、フーリエ変換された前記画像から、前記物体光の空間スペクトルを抽出し、抽出された前記物体光の空間スペクトルを逆フーリエ変換することで、前記物体光の複素振幅を求め、前記物体光の複素振幅から０次回折光成分を求め、フーリエ変換された前記画像から０次回折光成分を除去し、０次回折光成分が除去されたフーリエ変換された前記画像から、前記蛍光の像の空間スペクトルを抽出する。

上記の構成によれば、ホログラムと蛍光の像とが重畳した画像から、精度よく物体光の情報および蛍光の情報を分離することができる。

前記デジタルホログラフィ再生装置は、抽出された前記蛍光の像の空間スペクトルを逆フーリエ変換することで、前記蛍光の像を得る構成であってもよい。

上記の構成によれば、ホログラムと蛍光の像とが重畳した画像から、物体光および蛍光の像を分離して再生することができる。

本発明の一態様に係るデジタルホログラフィ再生装置は、参照光と物体光とが干渉することにより形成されるホログラムと、蛍光の像とが重畳した画像をフーリエ変換し、フーリエ変換された前記画像から、前記物体光の空間スペクトルを抽出し、抽出された前記物体光の空間スペクトルを逆フーリエ変換することで、前記物体光の複素振幅を求め、前記物体光の複素振幅から前記ホログラムの強度を求め、前記画像から前記ホログラムを除去することにより、前記蛍光の像を得る。

本発明の一態様に係るデジタルホログラフィ記録方法は、光源から出射された物体照明光を物体に照射することにより物体光を生じさせるとともに、前記物体照明光によって前記物体に含まれる蛍光体を励起するステップと、参照光と前記物体光とが干渉することにより形成されるホログラムと、前記蛍光体が発する蛍光の像とを撮像するステップとを含む。

本発明の一態様に係るデジタルホログラフィ記録方法は、物体照明光を物体に照射するステップと、参照光と前記物体からの物体光とが干渉することにより形成されるホログラムと、前記物体に含まれる蛍光体が発する蛍光の像とが重畳した画像を撮像するステップとを含む。

本発明の一態様に係るデジタルホログラフィ再生方法は、参照光と物体光とが干渉することにより形成されるホログラムと、蛍光の像とが重畳した画像から、空間的位相シフト法を用いて前記物体光の複素振幅を求めるステップと、前記物体光の複素振幅から前記ホログラムの強度を求めるステップと、前記画像から前記ホログラムを除去することにより、前記蛍光の像を得るステップとを含む。

本発明の一態様に係るデジタルホログラフィ再生方法は、参照光と物体光とが干渉することにより形成されるホログラムと、蛍光の像とが重畳した画像をフーリエ変換するステップと、フーリエ変換された前記画像から、前記物体光の空間スペクトルを抽出するステップと、抽出された前記物体光の空間スペクトルを逆フーリエ変換することで、前記物体光の複素振幅を求めるステップと、前記物体光の複素振幅から０次回折光成分を求め、フーリエ変換された前記画像から０次回折光成分を除去するステップと、０次回折光成分が除去されたフーリエ変換された前記画像から、前記蛍光の像の空間スペクトルを抽出するステップとを含む。

本発明の一態様に係るデジタルホログラフィ再生方法は、参照光と物体光とが干渉することにより形成されるホログラムと、蛍光の像とが重畳した画像をフーリエ変換するステップと、フーリエ変換された前記画像から、前記物体光の空間スペクトルを抽出するステップと、抽出された前記物体光の空間スペクトルを逆フーリエ変換することで、前記物体光の複素振幅を求めるステップと、前記物体光の複素振幅から前記ホログラムの強度を求め、前記画像から前記ホログラムを除去することにより、前記蛍光の像を得るステップとを含む。

本発明は上述した各実施形態に限定されるものではなく、請求項に示した範囲で種々の変更が可能であり、異なる実施形態にそれぞれ開示された技術的手段を適宜組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。

本発明は、デジタルホログラフィ装置、ハイスピード３Ｄカメラ、および蛍光３次元顕微鏡等に利用することができる。

１〜３ デジタルホログラフィ装置
１０、２０、３０ 記録装置（デジタルホログラフィ記録装置）
１１ 再生装置（デジタルホログラフィ再生装置）
１２、３２ 撮像装置
１３ 被写体（物体）
１４ 蛍光体
１５ 結像素子（結像光学素子）
ＬＳ１〜ＬＳ３ レーザ光源（第１光源、第２光源）

Claims (17)

1. 物体光が生じるように物体に物体照明光を照射する光源と、
   参照光と前記物体光とが干渉することにより形成され、画素に応じて変化する前記参照光および前記物体光についての位相シフト量が導入されたホログラムと、蛍光の像とを同時に撮像する撮像装置とを備え、
   前記物体照明光はまた、前記物体に含まれる蛍光体を励起することを特徴とするデジタルホログラフィ記録装置。
2. 前記光源は、単一波長の前記物体照明光を出射することを特徴とする請求項１に記載のデジタルホログラフィ記録装置。
3. 前記撮像装置は、前記物体光と参照光とが干渉することにより形成されるホログラムと、前記蛍光体が発する蛍光の像とが重畳した画像を撮像することを特徴とする請求項１または２に記載のデジタルホログラフィ記録装置。
4. 前記蛍光を前記撮像装置の撮像面に結像させる結像光学素子を備えることを特徴とする請求項３に記載のデジタルホログラフィ記録装置。
5. 前記撮像装置は、カラーフィルタを用いないモノクロ撮像装置であることを特徴とする請求項１から４のいずれか一項に記載のデジタルホログラフィ記録装置。
6. 前記光源は、前記物体に第１波長の物体照明光を照射する第１光源であり、
   前記第１波長とは異なる第２波長の物体光が生じるように前記物体に前記第２波長の物体照明光を照射する第２光源を備え、
   前記撮像装置は、前記第１波長のホログラムと、前記第２波長のホログラムと、前記蛍光の像とを撮像し、
   前記第１波長の物体照明光および前記第２波長の物体照明光の両方は、前記蛍光体を励起することを特徴とする請求項１から５のいずれか一項に記載のデジタルホログラフィ記録装置。
7. 前記撮像装置に入射する前記物体光の光軸に対して、前記撮像装置に入射する前記参照光の光軸は傾いていることを特徴とする請求項１から６のいずれか一項に記載のデジタルホログラフィ記録装置。
8. 物体に物体照明光を照射する光源と、
   参照光と前記物体からの物体光とが干渉することにより形成され、画素に応じて変化する前記参照光および前記物体光についての位相シフト量が導入されたホログラムと、前記物体に含まれる蛍光体が発する蛍光の像とを同時に撮像することにより、前記ホログラムと前記蛍光の像とが重畳した画像を撮像する撮像装置とを備えることを特徴とするデジタルホログラフィ記録装置。
9. 参照光と物体光とが干渉することにより形成されるホログラムと、インコヒーレント光の像とが重畳した画像から、空間的位相シフト法を用いて前記物体光の複素振幅を求め、
   前記物体光の複素振幅から前記ホログラムの強度を求め、
   前記画像から前記ホログラムを除去することにより、前記インコヒーレント光の像を得ることを特徴とするデジタルホログラフィ再生装置。
10. 参照光と物体光とが干渉することにより形成されるホログラムと、インコヒーレント光の像とが重畳した画像をフーリエ変換し、
    フーリエ変換された前記画像から、前記物体光の空間スペクトルを抽出し、
    抽出された前記物体光の空間スペクトルを逆フーリエ変換することで、前記物体光の複素振幅を求め、
    前記物体光の複素振幅から０次回折光成分を求め、フーリエ変換された前記画像から０次回折光成分を除去し、
    ０次回折光成分が除去されたフーリエ変換された前記画像から、前記インコヒーレント光の像の空間スペクトルを抽出することを特徴とするデジタルホログラフィ再生装置。
11. 抽出された前記インコヒーレント光の像の空間スペクトルを逆フーリエ変換することで、前記インコヒーレント光の像を得ることを特徴とする請求項１０に記載のデジタルホログラフィ再生装置。
12. 参照光と物体光とが干渉することにより形成されるホログラムと、インコヒーレント光の像とが重畳した画像をフーリエ変換し、
    フーリエ変換された前記画像から、前記物体光の空間スペクトルを抽出し、
    抽出された前記物体光の空間スペクトルを逆フーリエ変換することで、前記物体光の複素振幅を求め、
    前記物体光の複素振幅から前記ホログラムの強度を求め、前記画像から前記ホログラムを除去することにより、前記インコヒーレント光の像を得ることを特徴とするデジタルホログラフィ再生装置。
13. 光源から出射された物体照明光を物体に照射することにより物体光を生じさせるとともに、前記物体照明光によって前記物体に含まれる蛍光体を励起するステップと、
    参照光と前記物体光とが干渉することにより形成され、画素に応じて変化する前記参照光および前記物体光についての位相シフト量が導入されたホログラムと、前記蛍光体が発する蛍光の像とを同時に撮像するステップとを含むことを特徴とするデジタルホログラフィ記録方法。
14. 物体照明光を物体に照射するステップと、
    参照光と前記物体からの物体光とが干渉することにより形成され、画素に応じて変化する前記参照光および前記物体光についての位相シフト量が導入されたホログラムと、前記物体に含まれる蛍光体が発する蛍光の像とを同時に撮像することにより、前記ホログラムと前記蛍光の像とが重畳した画像を撮像するステップとを含むことを特徴とするデジタルホログラフィ記録方法。
15. 参照光と物体光とが干渉することにより形成されるホログラムと、インコヒーレント光の像とが重畳した画像から、空間的位相シフト法を用いて前記物体光の複素振幅を求めるステップと、
    前記物体光の複素振幅から前記ホログラムの強度を求めるステップと、
    前記画像から前記ホログラムを除去することにより、前記インコヒーレント光の像を得るステップとを含むことを特徴とするデジタルホログラフィ再生方法。
16. 参照光と物体光とが干渉することにより形成されるホログラムと、インコヒーレント光の像とが重畳した画像をフーリエ変換するステップと、
    フーリエ変換された前記画像から、前記物体光の空間スペクトルを抽出するステップと、
    抽出された前記物体光の空間スペクトルを逆フーリエ変換することで、前記物体光の複素振幅を求めるステップと、
    前記物体光の複素振幅から０次回折光成分を求め、フーリエ変換された前記画像から０次回折光成分を除去するステップと、
    ０次回折光成分が除去されたフーリエ変換された前記画像から、前記インコヒーレント光の像の空間スペクトルを抽出するステップとを含むことを特徴とするデジタルホログラフィ再生方法。
17. 参照光と物体光とが干渉することにより形成されるホログラムと、インコヒーレント光の像とが重畳した画像をフーリエ変換するステップと、
    フーリエ変換された前記画像から、前記物体光の空間スペクトルを抽出するステップと、
    抽出された前記物体光の空間スペクトルを逆フーリエ変換することで、前記物体光の複素振幅を求めるステップと、
    前記物体光の複素振幅から前記ホログラムの強度を求め、前記画像から前記ホログラムを除去することにより、前記インコヒーレント光の像を得るステップとを含むことを特徴とするデジタルホログラフィ再生方法。

Images