(57)【要約】【課題】物理厚さデータを推定することのできる物理厚さ推定プログラム、及びこれを用いた物理厚さデータ推定システム、並びに、これらを用いる屈折率データ推定プログラム、並びに、屈折率データ推定システムを提供する。【解決手段】本発明に係る屈折率データ推定プログラムは、コンピュータに、ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップ、ホログラム画像データに基づき測定対象物の位相差データを求めるステップ、前記物理厚さデータ及び前記位相差データに基づき前記測定対象物の屈折率データを求めるステップ、を実行させることを特徴とする。【選択図】  図１
2018128400.tif 000002
【特許請求の範囲】
【請求項１】
コンピュータに、
ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップ、
ホログラム画像データに基づき測定対象物の位相差データを求めるステップ、
前記物理厚さデータ及び前記位相差データに基づき前記測定対象物の屈折率データを求めるステップ、を実行させるための、屈折率データ推定プログラム。
【請求項２】
前記屈折率データに基づき、屈折率マップデータを作成する請求項１記載の屈折率データ推定プログラム。
【請求項３】
前記ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップは、
前記ホログラム画像データに基づき焦点位置の異なる複数の再生画像データを作成し、
前記複数の再生画像データ各々における焦点位置抽出画像データを作成し、
複数の前記焦点位置抽出画像データに基づき前記物理厚さを求める請求項１記載の屈折率データ推定プログラム。
【請求項４】
測定対象物を配置する測定対象物配置部と、
前記測定対象物配置部に配置された測定対象物を撮影し、ホログラム画像データを取得するためのホログラム画像データ取得装置と、
前記画像データ取得装置に接続され、前記画像データの処理を行う情報処理装置と、を備える屈折率推定システムであって、
前記情報処理装置は、ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップ、ホログラム画像データに基づき測定対象物の位相差データを求めるステップ、前記物理厚さデータ及び前記位相差データに基づき前記測定対象物の屈折率データを求めるステップ、を実行する、屈折率データ推定システム。
【請求項５】
コンピュータに、
ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップ、を実行させるための物理厚さデータ推定プログラム。
【請求項６】
測定対象物を配置する測定対象物配置部と、
前記測定対象物配置部に配置された測定対象物を撮影し、ホログラム画像データを取得するためのホログラム画像データ取得装置と、
前記画像データ取得装置に接続され、前記画像データの処理を行う情報処理装置と、を備える物理厚さデータ推定システムであって、
前記情報処理装置は、ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップ、を実行する、物理厚さデータ推定システム。  
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、物理厚さ推定プログラム及び屈折率データ推定プログラム、並びに、物理厚さデータ推定システム及び屈折率データ推定システムに関する。
【背景技術】
【０００２】
現在、顕微鏡下で、様々な物質の屈折率分布を計測することが多くの産業において求められている。例えば、細胞評価の分野においては、細胞の屈折率分布を得ることで、細胞を構成する成分量の定量化を行うことができる。すなわち、従来の技術によると、試料の光学的厚さを計測することは可能であった。例えば下記特許文献１には、光学的な厚さを得ることのできる細胞観察装置が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】ＷＯ２０１１／１３２５８７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、上記特許文献１に記載のように、従来の装置では光学的厚さを図ることができるとしても、屈折率の値の空間的な分布を得ることはできず、物理的な厚さを得ることは極めて困難であった。
【０００５】
そこで、本発明は、上記課題に鑑み、物理厚さデータを推定することのできる物理厚さ推定プログラム、及びこれを用いた物理厚さデータ推定システム、並びに、これらを用いる屈折率データ推定プログラム、並びに、屈折率データ推定システムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
上記課題を解決する本発明の一観点に係る屈折率データ推定プログラムは、コンピュータに、ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップ、ホログラム画像データに基づき測定対象物の位相差データを求めるステップ、物理厚さデータ及び前記位相差データに基づき前記測定対象物の屈折率データを求めるステップ、を実行するためのものである。
【０００７】
また、本発明の他の一観点に係る屈折率データ推定システムは、測定対象物を配置する測定対象物配置部と、測定対象物配置部に配置された測定対象物を撮影し、ホログラム画像データを取得するためのホログラム画像データ取得装置と、画像データ取得装置に接続され、画像データの処理を行う情報処理装置と、を備える屈折率推定システムであって、情報処理装置は、ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップ、ホログラム画像データに基づき測定対象物の位相差データを求めるステップ、物理厚さデータ及び位相差データに基づき測定対象物の屈折率データを求めるステップ、を実行するものである。
【０００８】
また、本発明の他の一観点に係る物理厚さデータ推定プログラムは、コンピュータに、ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップ、を実行させるためのものである。
【０００９】
また、本発明の他の一観点に係る物理厚さ推定システムは、測定対象物を配置する測定対象物配置部と、測定対象物配置部に配置された測定対象物を撮影し、ホログラム画像データを取得するためのホログラム画像データ取得装置と、画像データ取得装置に接続され、前記画像データの処理を行う情報処理装置と、を備え、情報処理装置は、ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップ、を実行するものである。
【発明の効果】
【００１０】
以上、本発明により、物理厚さデータを推定することのできる物理厚さ推定プログラム、及びこれを用いた物理厚さデータ推定システム、並びに、これらを用いる屈折率データ推定プログラム、並びに、屈折率データ推定システムを提供することができる。
【図面の簡単な説明】
【００１１】
【図１】実施形態に係る物理厚さデータ推定システムの概略を示す図である。
【図２】実施形態に係る測定対象物配置部のイメージを示す図である。
【図３】実施形態において焦点高さの異なる複数の再生画像データを作成する場合のイメージ図である。
【図４】実施形態において各再生画像データに対し輪郭抽出処理を行う場合のイメージ図である。
【図５】実施形態における位相差データのイメージを示す図である。
【図６】実施形態における応用例に係るミドリムシの体内における脂質の分布のイメージ図である。
【図７】実施例において作成した再生画像を示す図である。
【図８】実施例において作成した焦点位置抽出画像を示す図である。
【図９】実施例において作成した焦点位置抽出画像で塗りつぶした場合を示す図である。
【図１０】実施例において作成した焦点位置抽出画像を合成した場合の図である。
【図１１】実施例において作成した屈折率差のマップ画像を示す図である。
【発明を実施するための最良の形態】
【００１２】
以下、本発明の実施形態について図面を用いて詳細に説明する。ただし、本発明は多くの異なる形態による実施が可能であり、以下に示す実施形態、実施例における具体的な例示にのみ限定されるわけではない。
【００１３】
（物理厚さデータ推定システム）
図１は、本実施例に係る物理厚さデータ推定システム（以下「本システム」という。）Ｓ１の概略を示す図である。本図で示すように、本システムＳ１は、測定対象物Ｍを配置する測定対象物配置部１と、対象物配置部１に配置された測定対象物Ｍを撮影しホログラム画像データを取得するためのホログラム画像データ取得装置２と、ホログラム画像データ取得装置２に接続されホログラム画像データの処理を行う情報処理装置３と、を備える物理厚さデータ推定システムである。また、本システムＳ１において、情報処理装置３は、（Ｐ１）ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップ、を実行することができる。
【００１４】
本システムＳ１において、測定対象物配置部１は、文字通り測定対象物Ｍを保持することができるものである。ここで測定対象物Ｍは、特に限定されるわけではないが、光を透過することができるものであることが好ましい。具体的には、植物や動物等の生物細胞、ミドリムシ、ゾウリムシ及びミジンコ等の動物プランクトン、及び、ミカヅキモ、ハネケイソウ、ボルボックス等の植物性プランクトンといったプランクトン等を例示することができる。また、生物以外であっても、微小なレンズ等の透過物体を含む工業製品等に対しても応用が可能である。なお、光透過性の低い物体の場合や、光を反射する物体である場合、反射顕微鏡用のデジタルホログラムを用いることで、測定対象物とすることは可能である。
【００１５】
また、本システムＳ１における測定対象物配置部１の具体的な構造としては、測定時において安定的に測定対象物を配置することができる限りにおいて限定されるわけではないが、例えば平滑な表面を備えた基板、光透過性を備え測定対象物を内部に保持することができる中空の管、容器等を用いることができる。また、測定対象物配置部１が中空の管である場合、この中空の部分に測定対象物１を流す構成としてもよい。測定対象物１を流すことのできる構成とすれば、流れに合わせて連続的に撮影し、大量の測定対象物を高速に処理することが可能となるといった利点がある。
【００１６】
また、本システムＳ１における測定対象物配置部１において、測定対象物Ｍは、既知の屈折率からなる流体により均一に囲まれた状態となっていることが好ましい。この理由については後述の記載から明らかとなるが、このようにすることで、位相差を求める精度を向上させることができる。この場合において、流体としては、測定対象物Ｍを安定的に保持し、測定対象物Ｍの状態に対し影響を与えないものであることが好ましく、例えばプランクトン等の生体である場合、純水や生理食塩水であることが好ましい例であり、空気中で安定的に配置することができるのであれば空気、窒素等の非酸化性のガス等であることが好ましい。この場合のイメージを図２に示しておく。なお本図では、スライドガラスに測定対象物としてミドリムシ等の微生物を配置した場合のイメージを示している。
【００１７】
また、本システムＳ１において、ホログラム画像データ取得装置２は、上記の通り、対象物配置部１に配置された測定対象物Ｍを撮影し、ホログラム画像データを取得するために用いられるものである。またここで「ホログラム画像データ」とは、測定対象物の情報を干渉縞として記録した画像データであり、具体的には、位置データと、この位置データそれぞれにおける位相データ及び振幅データの組を複数備えたものである。ホログラム画像データは、複数の位置データを備えているが、イメージとして、ホログラム画像は複数の画素の集合であると考えることがデータ処理上非常に効率的であるため、位置データは、複数の画素の位置データとして考えることが好ましい。本システムＳ１では、ホログラム画像データに基づき、回折計算を施すことによって、三次元的な形状を定量的に求めることが可能となる。
【００１８】
また、ホログラム画像データ取得装置２の具体的な構造としては、限定されるわけではないが、例えば、レーザー光を発する光源と、このレーザー光を参照光と測定対象物照射用の光にそれぞれ照射するよう分割するビームスプリッタと、参照光と測定対象物に照射され反射した光を干渉させるための光学部材と、この干渉後の光を画像として記録する撮像装置と、を備えたものを例示することができる。
【００１９】
また、本システムＳ１では、微小な物の厚みをより正確に測定する観点から、上記構成にさらに拡大レンズ（対物レンズ）を備えた、いわゆるホログラフィック顕微鏡を用いることが好ましい。
【００２０】
また、本システムＳ１において、情報処理装置３は、上記の通り、ホログラム画像データ取得装置２に接続され、ホログラム画像データの処理を行うことができる装置である。
【００２１】
また、本システムＳ１の情報処理装置３としては、入力されたデータに基づき所望の計算処理を施し、所望のデータを得るために用いられる装置であって、本システムＳ１の機能を達成するための特定の処理及びデータの格納を行うことができる集積回路を組み合わせて基板上に配置した半導体基板を備えた特別の装置であってもよいが、いわゆるコンピュータであることは、データ処理の最適化のための改変、様々なデータ処理機能の付加を容易に行うことができる観点から好ましい。
【００２２】
本システムＳ１の情報処理装置３がコンピュータである場合、限定されるわけではないが、例えば中央演算装置（ＣＰＵ）、ハードディスクやメモリ等の記録媒体、キーボードやマウス等の入力装置、液晶ディスプレイ等の表示装置、これらを電気的に接続するために用いられるバスライン、を備えていることが好ましい。そして、上記の記録媒体に、所定の処理を実行させるためのプログラムを格納し、使用者の要求に応じてこのプログラムを実行することで所望の処理を実行することができ、必要に応じ、様々なデータ処理機能等を付加させやすくなる。
【００２３】
本システムＳ１の情報処理装置３は、上記の通り、（Ｐ１）ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップ、を実行することができる。
【００２４】
またこのステップは、限定されるわけではないが、例えば、（Ｐ１−１）ホログラム画像データに基づき焦点高さの異なる複数の再生画像データを作成し、（Ｐ１−２）複数の再生画像データ各々における焦点位置抽出画像データを作成し、（Ｐ１−３）複数の焦点位置抽出画像データに基づき物理厚さを求める、という細かなステップにより構成されていることが好ましい。
【００２５】
ここでより具体的に説明すると、まず、（Ｐ１−１）ホログラム画像データに基づき焦点高さの異なる複数の再生画像データを作成する。この場合のイメージ図を図３に示しておく。
【００２６】
本システムＳ１においてホログラム画像データは、各画素の位置において位相データ及び振幅データを備えているため、焦点の高さを変更しながら二次元の再生画像データとしてそれぞれ再生することができる。ここで「焦点の高さを変更しながら再生する」とは、特定の方向に沿った異なる位置（具体的には撮影方向に沿った異なる深さ）に焦点を合わせた二次元の再生画像を複数作成することをいう。ここで、特定の方向としては特に限定されるわけではないが、画像の深さ方向であることが好ましい。またこの場合において、焦点の高さの間隔を「焦点間隔」という。
【００２７】
また、本システムＳ１では、（Ｐ１−２）複数の再生画像データ各々における焦点位置抽出画像データを作成する。上記の処理において作成された再生画像データには、測定対象物において焦点が合っている部分と焦点が合っていない部分が存在する。ここで焦点が合っている位置は、測定対象物の輪郭の部分であり、他の部分と明確に区別することができる。そのため、本処理では、焦点が合っている位置を抽出することで、その焦点位置における輪郭抽出処理が可能となる。
【００２８】
より具体的に本処理は、各再生画像データに対し、輪郭抽出処理を行い、焦点位置抽出画像データを作成する。このイメージを図４に示しておく。
【００２９】
またこの処理において、輪郭抽出処理を行った後、この輪郭内の領域に所定の値を付しておく処理を行うことが好ましい。より具体的には、この輪郭内の領域を塗りつぶしておくことが好ましい。このようにすることで、後に各再生画像データを合成したときに、その値により物理的厚さを判定することができるようになる。
【００３０】
また、本システムＳ１では、（Ｐ１−３）複数の前記焦点位置抽出画像データに基づき物理厚さデータを求める処理を行う。具体的には、上記処理により得た各再生画像データを合成する、更に具体的には、上記作成した各再生画像データの各画素データにおける値を合計し、その合計値に焦点間隔に応じた値を乗じて、各画素における物理的厚みデータを得ることができる。
【００３１】
以上、本システムＳ１により、物理厚さデータを求めることのできる物理厚さ推定システム及びこれに用いられる物理厚さ推定プログラムを実現することができる。
【００３２】
（屈折率推定システム）
次に、上記物理厚さ推定システムを用いて、屈折率推定を行う屈折率推定システム（以下「本システム」という。）Ｓ２について説明する。本システムＳ２の概略は、基本的には上記図１と同様の構成となっているが、情報処理装置３の内部処理が異なる。この点について説明する。
【００３３】
本システムＳ２における情報処理装置３は、（Ｐ１）ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップ、（Ｐ２）ホログラム画像データに基づき測定対象物の位相差データを求めるステップ、（Ｐ３）物理厚さデータ及び位相差データに基づき前記測定対象物の屈折率データを求めるステップ、を実行する。なお本システムＳ２における（Ｐ１）ホログラム画像データに基づき測定対象物の物理厚さデータを求めるステップは、上記システムＳ１と同様であるため、その説明は省略する。
【００３４】
本システムＳ２によると、後述の説明から明らかとなるが、測定対象物の屈折率を推定することが可能となる。これは測定対象物全体の屈折率を求めることもでき、また、測定対象物の一部領域における屈折率を求めることができる。もちろん、本システムＳ２によると、屈折率そのものだけでなく、測定対象物の周囲に存在する流体との屈折率差を求めておくこともできる。その意味において「屈折率」には「屈折率差」を含むものとする。
【００３５】
本システムＳ２における情報処理装置３は、上記の通り（Ｐ２）ホログラム画像データに基づき測定対象物の位相差データを求めるステップを実行する。
【００３６】
ここで「位相差データ」とは、測定対象物の位置における位相と、測定対象物を透過していない部分の光の位相の差に関するデータをいう。この場合における位相差のイメージを図５に示しておく。
【００３７】
本図で示すように、より具体的に説明すると、測定対象物が存在する位置における位相をφｍ、測定対象物が存在していない位置の位相をφｒｅｆとすると、この位相差は、下記式であらわされる。
【数１】
2018128400.tif 000003
【００３８】
ここで、測定対象物の存在していない位置の位相については、限定されるわけではないが、例えば測定対象物の存在していない特定の領域を抽出し、その箇所における位相の値の平均値を取ることが精度向上の観点から好ましい。
【００３９】
また本システムＳ２における情報処理装置３は、上記の通り、（Ｐ３）物理厚さデータ及び前記位相差データに基づき前記測定対象物の屈折率データを求めるステップを実行する。
【００４０】
ここで、測定対象物の屈折率データを求める方法について上記図５を改めて用いて説明する。
【００４１】
まず、測定対象物が存在する位置及び存在しない位置にそれぞれ波長λの光が侵入する。すると、存在しない位置では、流体のみの位相変化を受けることにある。この結果、測定対象物が存在しない領域の位相はφｒｅｆとなる。一方、存在する領域では、流体と測定対象物を通り抜け、位相はφｍとなる。この位相の差が上記位相差となる。なおこの位相差は、純粋に測定対象物の物理的高さｄに依存した値となる。
【００４２】
一方、流体の屈折率は予め求めることができるため、ｎｒｅｆは既知の値にできる。測定対象物の屈折率をｎｍとすると、光路差は、下記の通りとなる。
【数２】
2018128400.tif 000004
【００４３】
また光路差は、上記位相差を用いて表すことも可能であるため、下記の指揮を満たすことになる。
【数３】
2018128400.tif 000005
【００４４】
そして、上記式において、未知数は測定対象物の屈折率ｎｒｅｆのみとなるため、式を整理すると下記のようになる。
【数４】
2018128400.tif 000006
【００４５】
以上の通り、測定対象物の物理厚さｄを取得することによって、測定対象物の屈折率を求めることができる。しかも、この屈折率は測定対象物の各位置において可能であるため、測定対象内の屈折率分布を求めることが可能である。この場合、屈折率分布を表すマップ画像データとして得ることが可能である。
【００４６】
もちろん、ｎｍ−ｎｒｅｆをそのまま屈折率差として求めておくことも可能である。
【００４７】
以上、本実施形態により、物理厚さデータを推定することのできる物理厚さ推定プログラム、及びこれを用いた物理厚さデータ推定システム、並びに、これらを用いる屈折率データ推定プログラム、並びに、屈折率データ推定システムを提供することができる。
【００４８】
ところで、上記システムを用いると、測定対象物内の屈折率及びその分布を求めることが可能となるため、様々な判断を行うことができ、例えば、測定対象物をミドリムシとし、このミドリムシの屈折率及びその分布を求めることで、個体間の脂質含有量を判断することができるようになる。これを用いた応用について以下具体的に説明する。
【００４９】
ミドリムシは体内に脂質を有しており、この脂質は高エネルギーであるため、バイオ燃料等に応用が期待されている。そして、体内に脂質を多く含有するミドリムシを選択的に抽出し、繁殖させ、更に高い資質を含むミドリムシを選択するといった工程を繰り返すことで、非常に高い脂質含有量のミドリムシを得ることができるようになる。
【００５０】
そして、この脂質はミドリムシ本体とは異なる組織成分を有しているため屈折率が異なる。より具体的には脂質の屈折率はそれ以外の組織成分よりも高いと考えられているため、ミドリムシ全体の屈折率が高いほど、脂質の割合が大きいと判断できる。この場合のイメージを図６に示しておく。
【００５１】
以上、本システムは、様々な応用が可能である。
【実施例】
【００５２】
ここで、上記システムに関し、二匹のミドリムシに対し、上記に従い実際の測定を行いその効果を確認した。測定は、スライドガラス上にミドリムシと水を配置し、ホログラフィック顕微鏡によってホログラム画像を取得し、各処理を施すことで行った。
【００５３】
まず、ホログラム画像データから、焦点高さの異なる複数の再生画像データを作成した。この作製した結果のうちの数枚を図７に示す。
【００５４】
また、この再生画像から焦点位置抽出画像データを作成した。具体的には、輪郭抽出を行った。この結果の画像を図８に示す。なお、この抽出した輪郭内について塗りつぶした画像を図９に示しておく。なお、この焦点位置抽出画像の焦点間隔は０．２９μｍである。
【００５５】
そして、上記塗りつぶした画像を合成し、高さを明暗で表した。この結果を図１０に示しておく。この結果、厚さとしては８．７μｍであった。
【００５６】
次に、屈折率差について算出した。具体的には、位相画像において、水部分の位相の平均を１．２７ｒａｄと求め、水の屈折率を１．３３、波長を０．５３２μｍとし、物理厚さをミドリムシのそれぞれの位置において屈折率差の分布を求めた。この屈折率差のマップ画像を図１１に示しておく。本図で示すミドリムシでは、屈折率差が濃淡により表現され、白く明るい部分が屈折率差の大きい部分すなわち脂質の割合が大きい部分と判定できた。この結果、ミドリムシ全体の屈折率は１．３５４４であると求めることができた。なおこれを、複数のミドリムシで行ったところ、様々な形状のミドリムシにそれぞれ異なる屈折率分布のものが存在していることが確認でき、より多くの脂質を含むミドリムシを選別できる可能性があることを確認した。
【００５７】
以上、本発明の効果について確認することができた。
【産業上の利用可能性】
【００５８】
本発明は、物理厚さ推定プログラム及び屈折率データ推定プログラム、並びに、物理厚さデータ推定システム及び屈折率データ推定システムとして産業上の利用可能性がある。
【図１】
2018128400.tif 000007
【図２】
2018128400.tif 000008
【図３】
2018128400.tif 000009
【図４】
2018128400.tif 000010
【図５】
2018128400.tif 000011
【図６】
2018128400.tif 000012
【図７】
2018128400.tif 000013
【図８】
2018128400.tif 000014
【図９】
2018128400.tif 000015
【図１０】
2018128400.tif 000016
【図１１】
2018128400.tif 000017