JP6739391B2 - 擬似スペックルパターン生成装置、擬似スペックルパターン生成方法、観察装置および観察方法 - Google Patents

擬似スペックルパターン生成装置、擬似スペックルパターン生成方法、観察装置および観察方法 Download PDF

Info

Publication number
JP6739391B2
JP6739391B2 JP2017077474A JP2017077474A JP6739391B2 JP 6739391 B2 JP6739391 B2 JP 6739391B2 JP 2017077474 A JP2017077474 A JP 2017077474A JP 2017077474 A JP2017077474 A JP 2017077474A JP 6739391 B2 JP6739391 B2 JP 6739391B2
Authority
JP
Japan
Prior art keywords
pattern
light
intensity modulation
pseudo
intensity
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2017077474A
Other languages
English (en)
Other versions
JP2018180177A (ja
Inventor
寛人 酒井
寛人 酒井
太郎 安藤
太郎 安藤
豊田 晴義
晴義 豊田
良幸 大竹
良幸 大竹
優 瀧口
優 瀧口
知子 兵土
知子 兵土
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hamamatsu Photonics KK
Original Assignee
Hamamatsu Photonics KK
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hamamatsu Photonics KK filed Critical Hamamatsu Photonics KK
Priority to JP2017077474A priority Critical patent/JP6739391B2/ja
Priority to PCT/JP2018/012130 priority patent/WO2018190111A1/ja
Priority to DE112018001943.8T priority patent/DE112018001943T5/de
Priority to US16/499,924 priority patent/US11300800B2/en
Priority to CN201880023619.0A priority patent/CN110494794B/zh
Publication of JP2018180177A publication Critical patent/JP2018180177A/ja
Application granted granted Critical
Publication of JP6739391B2 publication Critical patent/JP6739391B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/48Laser speckle optics
    • GPHYSICS
    • G02OPTICS
    • G02FOPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
    • G02F1/00Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
    • G02F1/01Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics for the control of the intensity, phase, polarisation or colour 
    • G02F1/13Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics for the control of the intensity, phase, polarisation or colour  based on liquid crystals, e.g. single liquid crystal display cells
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B21/00Microscopes
    • G02B21/06Means for illuminating specimens
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B26/00Optical devices or arrangements for the control of light using movable or deformable optical elements
    • G02B26/02Optical devices or arrangements for the control of light using movable or deformable optical elements for controlling the intensity of light

Description

本発明は、擬似スペックルパターン生成装置、この擬似スペックルパターン生成装置を備える観察装置、擬似スペックルパターン生成方法、および、この擬似スペックルパターン生成方法を含む観察方法に関するものである。
レーザ光などのコヒーレントな光が拡散板等のディフューザ(散乱媒体)により散乱されると、その散乱光の干渉によりスペックルパターンが生成される。スペックルパターンにおける光強度の空間分布が乱数パターンに近い性質を有していることから、スペックルパターンは、干渉計測技術、超解像技術および光計測技術などにおいて利用され、また、最近では原子や粒子コロイドのトラップ等の光操作技術においても利用されている。
スペックルパターンは、空間構造および光強度統計分布により特徴付けられる。空間構造は、スペックルパターンの自己相関関数の空間的形状により表され、光学系の点像分布を示す。光強度統計分布は、スペックルパターンにおける光強度のヒストグラムにより表され、指数分布に従う。なお、理想的な乱数を2次元配列した2次元擬似乱数パターンの空間構造(自己相関関数)はデルタ関数となる。
非特許文献1に記載された技術は、ディフューザおよびレンズを備える構成において、レンズの後焦点面と異なる面にスペックルパターンを生成する。非特許文献2に記載された技術は、ディフューザおよびレンズを備える構成において、ディフューザの位置の変更又はディフューザ自体の交換によって様々なスペックルパターンを生成する。
非特許文献3に記載された技術は、外部から与えられる信号により設定可能な2次元の変調分布を有する空間光変調器をディフューザに替えて用い、この空間光変調器およびレンズを備える構成によりスペックルパターンを生成する。空間光変調器には、2次元擬似乱数パターンに応じた変調分布が設定される。
Vladlen G. Shvedov et.al,"Selective trapping of multiple particles by volume speckle field," OPTICSEXPRESS, Vol.18, No.3, pp.3137-3142 (2010). E. Mudry et.al, "Structuredillumination microscopy using unknown speckle patterns," NATURE PHOTONICS, VOL6, pp.312-315 (2012). Yaron Bromberg and Hui Cao,"Generating Non-Rayleigh Speckles with Tailored Intensity Statistics," PHYSICALREVIEW LETTERS 112, 213904 (2014).
非特許文献1に記載された技術は、1つのディフューザを用いるのみであるから、生成されるスペックルパターンの空間構造および光強度統計分布の設定の自由度が無い。
非特許文献2に記載された技術は、ディフューザを交換することによってスペックルパターンを変更することができる。しかし、この技術においては、多様なスペックルパターンを生成するためには多くのディフューザを用意しておくことが必要であり、生成されるスペックルパターンの空間構造および光強度統計分布の設定の自由度が低い。また、ディフューザ位置の変更によりスペックルパターンを設定する場合、同一のパターンを再現することは著しく困難である。
非特許文献3に記載された技術は、2次元擬似乱数パターンに応じて空間光変調器に呈示する変調分布を変更することによってスペックルパターンを変更することができる。しかし、この技術においても、スペックルパターンのコントラストが調整されるのみであり、生成されるスペックルパターンの空間構造および光強度統計分布の設定の自由度が低い。
本発明は、上記問題点を解消する為になされたものであり、生成されるスペックルパターンの空間構造または光強度統計分布の設定の自由度が高く、生成されるパターンが再現性を有する擬似スペックルパターン生成装置および擬似スペックルパターン生成方法を提供することを目的とする。また、本発明は、このような擬似スペックルパターン生成装置を備える観察装置、および、擬似スペックルパターン生成方法を含む観察方法を提供することを目的とする。
なお、本発明で生成されるパターンは、ディフューザを用いて生成される通常のスペックルパターンと異なる空間構造または光強度統計分布を有することも可能であるので、「擬似スペックルパターン」と呼ぶことにする。
本発明の擬似スペックルパターン生成装置は、(1) 光を出力する光源と、(2) 擬似乱数パターンに基づく第1強度変調分布を有し、光源から出力された光を第1強度変調分布に応じて空間的に強度変調して、その変調後の光を出力する第1空間光変調器と、(3) 第1空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを生成する第1光学系と、(4) 第1光学系によりフーリエ変換パターンが生成される面に設けられ、フィルタ関数に基づく第2強度変調分布を有し、第1光学系を経て到達した光を第2強度変調分布に応じて空間的に強度変調して、その変調後の光を出力する第2空間光変調器と、(5) 第2空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを擬似スペックルパターンとして光パターン生成面に生成する第2光学系と、(6) 第1空間光変調器の第1強度変調分布および第2空間光変調器の第2強度変調分布の双方または何れか一方を設定する制御部と、を備える。
本発明の擬似スペックルパターン生成装置の一側面では、制御部は、生成すべき擬似スペックルパターンの光強度統計分布に応じた統計分布を有する擬似乱数パターンに基づいて第1強度変調分布を第1空間光変調器に設定するのが好適である。また、制御部は、生成すべき擬似スペックルパターンの自己相関関数に応じた相関関数のフーリエ変換の平方根であるフィルタ関数に基づいて第2強度変調分布を第2空間光変調器に設定するのが好適である。
本発明の観察装置は、(1) 上記の本発明の擬似スペックルパターン生成装置と、(2) 擬似スペックルパターン生成装置により擬似スペックルパターンが生成される光パターン生成面に照射される観察用の光を出力する観察用光源と、(3) 光パターン生成面への観察用の光の照射に応じて生じた光を受光して撮像するカメラと、を備える。
本発明の擬似スペックルパターン生成方法は、(1) 擬似乱数パターンに基づく第1強度変調分布を有する第1空間光変調器により、光源から出力された光を第1強度変調分布に応じて空間的に強度変調して、その変調後の光を出力し、(2) 第1光学系により、第1空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを生成し、(3) 第1光学系によりフーリエ変換パターンが生成される面に設けられフィルタ関数に基づく第2強度変調分布を有する第2空間光変調器により、第1光学系を経て到達した光を第2強度変調分布に応じて空間的に強度変調して、その変調後の光を出力し、(4) 第2光学系により、第2空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを擬似スペックルパターンとして光パターン生成面に生成し、(5) 第1空間光変調器の第1強度変調分布および第2空間光変調器の第2強度変調分布の双方または何れか一方を設定する。
本発明の擬似スペックルパターン生成方法の一側面では、生成すべき擬似スペックルパターンの光強度統計分布に応じた統計分布を有する擬似乱数パターンに基づいて第1強度変調分布を第1空間光変調器に設定するのが好適である。また、生成すべき擬似スペックルパターンの自己相関関数に応じた相関関数のフーリエ変換の平方根であるフィルタ関数に基づいて第2強度変調分布を第2空間光変調器に設定するのが好適である。
本発明の観察方法は、上記の本発明の擬似スペックルパターン生成方法により擬似スペックルパターンが生成される光パターン生成面に、観察用光源から出力される観察用の光を照射し、光パターン生成面への観察用の光の照射に応じて生じた光をカメラにより受光して撮像する。
本発明によれば、生成される擬似スペックルパターンの空間構造または光強度統計分布の設定の自由度が高く、生成されるパターンが再現性を有する。
図1は、擬似スペックルパターンを計算により求める場合の計算手順のフローチャートである。 図2は、相関関数cおよびフィルタ関数Fの例を纏めた表である。 図3(a)は、統計分布が正規分布に従う2次元擬似乱数パターンa(x,y)の一例を示す図である。図3(b)は、図3(a)の2次元擬似乱数パターンa(x,y)および相関関数c(r)=exp(-r)を用いた場合に得られる擬似スペックルパターンb(x,y)を示す図である。 図4(a)は、図3(a)の2次元擬似乱数パターンa(x,y)および相関関数c(r)=exp(-r/9)を用いた場合に得られる擬似スペックルパターンb(x,y)を示す図である。図4(b)は、図4(a)の擬似スペックルパターンb(x,y)の光強度統計分布(輝度ヒストグラム)を示す図である。 図5(a)は、ディフューザを用いて生成されたスペックルパターンの一例を示す図である。図5(b)は、図5(a)のスペックルパターンの光強度統計分布(輝度ヒストグラム)を示す図である。 図6(a)は、図3(a)の2次元擬似乱数パターンa(x,y)および相関関数c(r)=exp(-r/3)を用いた場合の擬似スペックルパターンb(x,y)を示す図である。図6(b)は、図6(a)の擬似スペックルパターンb(x,y)の自己相関画像を示す図である。 図7は、図6(b)の自己相関画像における強度プロファイルを示す図である。 図8は、本実施形態の擬似スペックルパターン生成装置1Aの構成を示す図である。 図9は、他の本実施形態の擬似スペックルパターン生成装置1Bの構成を示す図である。 図10は、本実施形態の観察装置2Aの構成を示す図である。 図11は、他の実施形態の観察装置2Bの構成を示す図である。
以下、添付図面を参照して、本発明を実施するための形態を詳細に説明する。なお、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。本発明は、これらの例示に限定されるものではなく、特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
先ず、本実施形態において生成される擬似スペックルパターンを説明する。図1は、擬似スペックルパターンを計算により求める場合の計算手順のフローチャートである。なお、以下に登場する関数のうち、英小文字で表した関数は実空間における関数であり、英大文字で表した関数はフーリエ空間における関数である。また、(x,y)は実空間中の直交座標系で表した位置を示し、(u,v)はフーリエ空間中の直交座標系で表した位置を示す。
ステップS11では、2次元擬似乱数パターンa(x,y)を生成する。具体的には、任意の生成方法により所定のシード数および統計分布に従う1次元の擬似乱数列を生成し、これらの擬似乱数を2次元に配列することで2次元擬似乱数パターンa(x,y)を生成する。1次元擬似乱数列の生成の際のシード数および統計分布ならびに2次元配列の際の配列規則が同じであれば、2次元擬似乱数パターンa(x,y)の生成は再現性を有する。ステップS12では、a(x,y)をフーリエ変換してA(u,v)を得る((1)式)。
ステップS21では、相関関数c(x,y)を準備する。ステップS22では、c(x,y)をフーリエ変換してC(u,v)を得る((2)式)。ステップS23では、C(u,v)の平方根を計算してフィルタ関数F(u,v)を得る((3)式)。
ステップS31では、A(u,v)とF(u,v)とを乗算してB(u,v)を得る((4)式)。そして、ステップS32では、B(u,v)を逆フーリエ変換して擬似スペックルパターンb(x,y)を得る((5)式)。なお、フーリエ変換と逆フーリエ変換とは、数値計算をする場合には同様の計算処理で実現可能であり、また、光学系により実現する場合には同様の光学系により実現可能であるから、両者を区別する必要はない。
このようにして生成された擬似スペックルパターンb(x,y)は、2次元擬似乱数パターンa(x,y)の統計分布に応じた光強度統計分布と、フィルタ関数F(u,v)に対応する相関関数c(x,y)に応じた空間構造(自己相関関数)と、を有する。なお、ステップS23においてC(u,v)の平方根を計算してフィルタ関数F(u,v)を求めるのは、擬似スペックルパターンb(x,y)の自己相関関数と相関関数c(x,y)とを互いに一致させるためである。
図2は、相関関数cおよびフィルタ関数Fの例を纏めた表である。この表に示される相関関数cは、実空間における原点(0,0)からの距離rのみを変数とする関数c(r)である。対応するフィルタ関数Fも、フーリエ空間における原点(0,0)からの距離kのみを変数とする関数F(k)である。各式において、r,kは非負の実数であり、η,ξは正の実数である。
相関関数c(x,y)としては、フーリエ変換の結果であるC(u,v)が非負の実数の関数となるような関数が選ばれ、また、距離rが大きくなるに従い関数値が漸減する関数であるのが好ましい。相関関数cは、距離rのみを変数とする関数であるのが好ましい。相関関数cは、xおよびyを変数とする関数(例えばexp(-|x+y|/ξ))であってもよい。相関関数cは、原点(0,0)における関数値が1であってもよいが、その関数値が任意であってよい。相関関数cは、数式で表現することができないものであってもよい。
図3(a)は、統計分布が正規分布に従う2次元擬似乱数パターンa(x,y)の一例を示す図である。図3(b)は、図3(a)の2次元擬似乱数パターンa(x,y)および相関関数c(r)=exp(-r)を用いた場合に得られる擬似スペックルパターンb(x,y)を示す図である。図4(a)は、図3(a)の2次元擬似乱数パターンa(x,y)および相関関数c(r)=exp(-r/9)を用いた場合に得られる擬似スペックルパターンb(x,y)を示す図である。図4(b)は、図4(a)の擬似スペックルパターンb(x,y)の光強度統計分布(輝度ヒストグラム)を示す図である。
図5(a)は、ディフューザを用いて生成されたスペックルパターンの一例を示す図である。図5(b)は、図5(a)のスペックルパターンの光強度統計分布(輝度ヒストグラム)を示す図である。ディフューザを用いて生成されたスペックルパターンの光強度統計分布(図5(b))は指数関数分布で近似することができるのに対して、擬似スペックルパターンb(x,y)の光強度統計分布(図4(b))は、正規分布で近似することができ、2次元擬似乱数パターンa(x,y)の統計分布と略一致している。なお、図4(b)および図5(b)の光強度統計分布(輝度ヒストグラム)の横軸は、256階調の光強度を表している。
図6(a)は、図3(a)の2次元擬似乱数パターンa(x,y)および相関関数c(r)=exp(-r/3)を用いた場合の擬似スペックルパターンb(x,y)を示す図である。図6(b)は、図6(a)の擬似スペックルパターンb(x,y)の自己相関画像を示す図である。図7は、図6(b)の自己相関画像における強度プロファイルを示す図である。x方向およびy方向それぞれの強度プロファイルは、互いに略一致しており、また、相関関数c(r)=exp(-r/3)とも略一致している。すなわち、擬似スペックルパターンb(x,y)の自己相関関数は、相関関数cと略一致している。
本実施形態の擬似スペックルパターン生成装置および擬似スペックルパターン生成方法では、2次元擬似乱数パターンa(x,y)に基づく第1強度変調分布を有する第1空間光変調器と、フィルタ関数F(u,v)に基づく第2強度変調分布を有する第2空間光変調器とを用いて、擬似スペックルパターンb(x,y)を光学的に生成する。
図8は、本実施形態の擬似スペックルパターン生成装置1Aの構成を示す図である。擬似スペックルパターン生成装置1Aは、制御部10、光源11、ビームエキスパンダ12、第1空間光変調器13、レンズ14、第2空間光変調器15およびレンズ16を備え、光パターン生成面Pに擬似スペックルパターンb(x,y)を生成する。第1空間光変調器13から光パターン生成面Pに至るまでの光学系は4f結像光学系を構成している。
光源11は、光を出力する。光源11として、例えば、レーザ光源、ランプ光源、SLD(SuperluminescentDiode)光源等が用いられる。ビームエキスパンダ12は、光源11と光学的に接続され、光源11から出力される光のビーム径を拡大して出力する。このとき、ビームエキスパンダ12から出力される光は、ビーム断面において強度が均一化されるのが好ましい。
第1空間光変調器13は、強度変調型のものであって、ステップS11で生成される2次元擬似乱数パターンa(x,y)に基づく第1強度変調分布を有する。第1空間光変調器13として、例えば、変調媒体が液晶である変調器や、デジタルミラーデバイス(DMD)が用いられる。第1空間光変調器13は、ビームエキスパンダ12と光学的に接続され、光源11から出力されてビームエキスパンダ12によりビーム径を拡大された光を入力し、その入力した光を第1強度変調分布に応じて空間的に強度変調して、その変調後の光を出力する。
レンズ14は、第1空間光変調器13と光学的に接続され、第1空間光変調器13から出力された光のパターンを光学的にフーリエ変換する第1光学系を構成している。このレンズ14により生成されるフーリエ変換パターンは、ステップS12で得られるA(u,v) に応じたものである。
第2空間光変調器15は、強度変調型のものであって、ステップS23で生成されるフィルタ関数F(u,v)に基づく第2強度変調分布を有する。第2空間光変調器15として、例えば、変調媒体が液晶である変調器や、デジタルミラーデバイス(DMD)が用いられる。第2空間光変調器15は、レンズ(第1光学系)14によりフーリエ変換パターンが生成される面に設けられ、レンズ14を経て到達した光を第2強度変調分布に応じて空間的に強度変調して、その変調後の光を出力する。第2空間光変調器15から出力される光のパターンは、ステップS31で得られるB(u,v)に応じたものである。
レンズ16は、第2空間光変調器15と光学的に接続され、第2空間光変調器15から出力された光のパターンを光学的にフーリエ変換する第2光学系を構成している。このレンズ16により光パターン生成面Pにおいて生成されるフーリエ変換パターンは、ステップS32で生成される擬似スペックルパターンb(x,y)となる。
このようにして生成される擬似スペックルパターンb(x,y)は、2次元擬似乱数パターンa(x,y)の統計分布(すなわち、第1空間光変調器13の第1強度変調分布)に応じた光強度統計分布を有するとともに、フィルタ関数F(u,v)(すなわち、第2空間光変調器15の第2強度変調分布)に対応する相関関数c(x,y)に応じた空間構造を有する。
制御部10は、第1空間光変調器13および第2空間光変調器15の双方または何れか一方と電気的に接続され、第1空間光変調器13の第1強度変調分布および第2空間光変調器15の第2強度変調分布の双方または何れか一方を設定する。具体的には、制御部10は、生成すべき擬似スペックルパターンb(x,y)の光強度統計分布に応じた統計分布を有する2次元擬似乱数パターンa(x,y)に基づいて第1強度変調分布を第1空間光変調器13に設定する。或いは、制御部10は、生成すべき擬似スペックルパターンb(x,y)の自己相関関数に応じた相関関数c(x,y)に対応するフィルタ関数F(u,v)に基づいて第2強度変調分布を第2空間光変調器15に設定する。
制御部10は、例えばコンピュータにより構成され、第1空間光変調器13および第2空間光変調器15の双方または何れか一方と電気的に接続されて通信を行う通信部を備える他、演算部、記憶部、入力部および表示部を備える。演算部は、CPUまたはDSP等を含み、2次元擬似乱数パターンa(x,y)や相関関数c(x,y)の生成、フーリエ変換、平方根計算、乗算、自己相関関数の演算、強度統計分布の算出等を行う。記憶部は、例えばハードディスクまたはメモリなどを含み、2次元擬似乱数パターンa(x,y)や相関関数c(x,y)の生成条件や各パターンの演算結果等を記憶する。入力部は、例えばキーボードまたはポインティングデバイスなどを含み、上記生成条件等の入力を受け付ける。表示部は、例えば液晶ディスプレイを含み、a(x,y)、A(u,v)、c(x,y)、C(u,v)、F(u,v)、B(u,v)およびb(x,y)等の各パターンを表示する。
本実施形態では、2つの空間光変調器(第1空間光変調器13、第2空間光変調器15)を用い、第1空間光変調器13の第1強度変調分布および第2空間光変調器15の第2強度変調分布の双方または何れか一方を設定することができるので、生成される擬似スペックルパターンb(x,y)の空間構造または光強度統計分布の設定の自由度が高い。生成すべき擬似スペックルパターンb(x,y)の光強度統計分布に応じた統計分布を有する2次元擬似乱数パターンa(x,y)に基づいて第1強度変調分布を第1空間光変調器13に設定することにより、所望の光強度統計分布を有する擬似スペックルパターンb(x,y)を生成することができる。また、生成すべき擬似スペックルパターンb(x,y)の自己相関関数に応じた相関関数c(x,y)のフーリエ変換の平方根であるフィルタ関数F(u,v)に基づいて第2強度変調分布を第2空間光変調器15に設定することにより、所望の空間構造を有する擬似スペックルパターンb(x,y)を生成することができる。
本実施形態の擬似スペックルパターン生成装置または擬似スペックルパターン生成方法では、生成される擬似スペックルパターンb(x,y)の空間構造または光強度統計分布の設定の自由度が高い。加えて、擬似スペックルパターンb(x,y)の空間構造または光強度統計分布を、制御部10により、再現性よく設定することができ、迅速に設定を変更することができる。本実施形態の擬似スペックルパターン生成装置または擬似スペックルパターン生成方法は、計測対象や光操作対象に適した擬似スペックルパターンb(x,y)を生成することができ、様々な検証方法に有効な計測技術や光操作技術において好適に用いることができる。また、本実施形態の擬似スペックルパターン生成装置または擬似スペックルパターン生成方法は、擬似スペックルパターンを構造化照明として用いる場合に、平均輝度が高く明るい擬似スペックルパターンを生成することもできるので、照明回数を少なくして、高感度で短時間の照明が可能となる。
本実施形態により所望の空間構造および光強度統計分布を有する擬似スペックルパターンb(x,y)を生成することができることは、次のようにして確認することができる。すなわち、光パターン生成面Pに撮像面を有するカメラを用いて、光パターン生成面Pに生成された擬似スペックルパターンb(x,y)を撮像する。そして、制御部10により、撮像により得られた擬似スペックルパターンb(x,y)の自己相関関数を求めて、これを相関関数c(x,y)と比較する。また、制御部10により、撮像により得られた擬似スペックルパターンb(x,y)の光強度統計分布を求めて、これを2次元擬似乱数パターンa(x,y)の統計分布と比較する。
図9は、他の本実施形態の擬似スペックルパターン生成装置1Bの構成を示す図である。擬似スペックルパターン生成装置1Bは、制御部10、光源11、ビームエキスパンダ12、位相変調型の空間光変調器17、偏光板18,19、レンズ14、第2空間光変調器15およびレンズ16を備え、光パターン生成面Pに擬似スペックルパターンb(x,y)を生成する。図8に示された擬似スペックルパターン生成装置1Aの構成と比較すると、この図9に示される擬似スペックルパターン生成装置1Bは、強度変調型の第1空間光変調器13に替えて位相変調型の空間光変調器17および偏光板18,19を備える点で相違する。位相変調型の空間光変調器17は、例えば、変調媒体が液晶である変調器が用いられる。
位相変調型の空間光変調器17を挟んで偏光板18および偏光板19が設けられている。一般に変調媒体が液晶である空間光変調器17は特定の配向方向の液晶分子を有する。偏光板18,19は、この液晶分子の配向方向に対して相対的に45°の角度である偏光方向となるよう配置される。偏光板18,19それぞれの偏光方向は、互いに平行であってもよいし、互いに垂直であってもよい。このように位相変調型の空間光変調器17および偏光板18,19が配置されることで、これらは実質的に強度変調型の第1空間光変調器13Bとして機能する。
制御部10は、第1空間光変調器13B(これを構成する空間光変調器17)および第2空間光変調器15の双方または何れか一方と電気的に接続され、第1空間光変調器13Bの第1強度変調分布(空間光変調器17の位相変調分布)および第2空間光変調器15の第2強度変調分布の双方または何れか一方を設定する。
第2空間光変調器15も、同様に、位相変調型の空間光変調器を挟んで2つの偏光板を設ける構成とすることができる。
なお、擬似スペックルパターン生成装置1A,1Bにおいて、擬似スペックルパターンb(x,y)が生成される光パターン生成面Pは、レンズ(第2光学系)16の後焦点面であってもよいが、その後焦点面と異なる面とすることもできる。例えば、第2強度変調分布にフレネルレンズパターンを加えたパターンを第2空間光変調器15に設定することにより、レンズ(第2光学系)16の後焦点面と異なる位置にある光パターン生成面Pに擬似スペックルパターンb(x,y)を生成することができる。
次に、本実施形態の擬似スペックルパターン生成装置を備える観察装置の実施形態、および、本実施形態の擬似スペックルパターン生成方法を含む観察方法の実施形態について説明する。本実施形態の観察装置または観察方法は、上記の本実施形態の擬似スペックルパターン生成装置または擬似スペックルパターン生成方法により擬似スペックルパターンが生成される光パターン生成面Pに観察用の光を照射し、この光パターン生成面Pへの光の照射に応じて発生する光(透過光、反射光、散乱光、蛍光など)を受光して、擬似スペックルパターンが照射された光パターン生成面P上の微小体等を撮像する。
図10は、本実施形態の観察装置2Aの構成を示す図である。観察装置2Aは、制御部10、光源11、ビームエキスパンダ12、第1空間光変調器13、レンズ14、第2空間光変調器15およびレンズ16を備え、光パターン生成面Pに擬似スペックルパターンを生成する。これらは、本実施形態の擬似スペックルパターン生成装置の構成と同様である。レンズ16として対物レンズが用いられる。
観察装置2Aは、擬似スペックルパターン生成装置の構成に加えて、ダイクロイックミラー21、観察用光源22、カメラ23およびレンズ24を更に備える。ダイクロイックミラー21は、第2空間光変調器15とレンズ16との間の光路上に設けられ、光源11から出力される光を透過させ、観察用光源22から出力される光を反射させる。
観察用光源22は、光源11から出力される光の波長と異なる波長の観察用の光を出力する。観察用光源22は、光パターン生成面Pに対してレンズ16が設けられる側と反対の側に設けられ、光パターン生成面Pに観察用の光を照射する。観察用光源22による観察用の光の照射に応じて発生した光は、レンズ16を経て、ダイクロイックミラー21により反射され、レンズ24を経て、カメラ23の撮像面に到達する。レンズ16およびレンズ24は、光パターン生成面Pで発生した光をカメラ23の撮像面に結像する。カメラ23は、擬似スペックルパターンが照射された光パターン生成面P上の微小体等を撮像する。カメラ23は、例えばCCDカメラやCMOSカメラ等である。
図11は、他の実施形態の観察装置2Bの構成を示す図である。観察装置2Bも、制御部10、光源11、ビームエキスパンダ12、第1空間光変調器13、レンズ14、第2空間光変調器15およびレンズ16を備え、光パターン生成面Pに擬似スペックルパターンを生成する。これらは、本実施形態の擬似スペックルパターン生成装置の構成と同様である。レンズ16として対物レンズが用いられる。
観察装置2Bは、擬似スペックルパターン生成装置の構成に加えて、ダイクロイックミラー21、観察用光源22、カメラ23、レンズ25、対物レンズ26およびレンズ27を更に備える。ダイクロイックミラー21は、第2空間光変調器15とレンズ16との間の光路上に設けられ、光源11から出力される光を透過させ、観察用光源22から出力される光を反射させる。
観察用光源22は、光源11から出力される光の波長と異なる波長の観察用の光を出力する。観察用光源22から出力された光は、レンズ25を経て、ダイクロイックミラー21により反射され、レンズ16を経て、光パターン生成面Pに照射される。カメラ23、対物レンズ26およびレンズ27は、光パターン生成面Pに対してレンズ16が設けられる側と反対の側に設けられている。観察用光源22による観察用の光の照射に応じて発生した光は、対物レンズ26およびレンズ27を経て、カメラ23の撮像面に到達する。対物レンズ26およびレンズ27は、光パターン生成面Pで発生した光をカメラ23の撮像面に結像する。カメラ23は、擬似スペックルパターンが照射された光パターン生成面P上の微小体等を撮像する。
これらの観察装置2A,2Bは、擬似スペックルパターン生成装置を光トラップ技術および構造化照明技術をはじめとする光操作技術あるいは光イメージング技術に応用する場合に、顕微イメージング装置として好適に用いられる。本実施形態の観察装置または観察方法は、擬似スペックルパターン生成装置または擬似スペックルパターン生成方法により生成される擬似スペックルパターンを微小体のトラップやソーティングなどの光操作に用いることで、様々な態様の光操作が可能となって、その光操作の様子を観察することができる。
なお、上記の実施形態では擬似乱数パターンおよび擬似スペックルパターンは2次元であるとしたが、本発明では擬似乱数パターンおよび擬似スペックルパターンは1次元または3次元であってもよい。
1A,1B…擬似スペックルパターン生成装置、2A,2B…観察装置、10…制御部、11…光源、12…ビームエキスパンダ、13…第1空間光変調器、14…レンズ(第1光学系)、15…第2空間光変調器、16…レンズ(第2光学系)、17…空間光変調器、18,19…偏光板、21…ダイクロイックミラー、22…観察用光源、23…カメラ、24…レンズ、25…レンズ、26…対物レンズ、27…レンズ、P…光パターン生成面。

Claims (8)

  1. 光を出力する光源と、
    擬似乱数パターンに基づく第1強度変調分布を有し、前記光源から出力された光を前記第1強度変調分布に応じて空間的に強度変調して、その変調後の光を出力する第1空間光変調器と、
    前記第1空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを生成する第1光学系と、
    前記第1光学系により前記フーリエ変換パターンが生成される面に設けられ、フィルタ関数に基づく第2強度変調分布を有し、前記第1光学系を経て到達した光を前記第2強度変調分布に応じて空間的に強度変調して、その変調後の光を出力する第2空間光変調器と、
    前記第2空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを擬似スペックルパターンとして光パターン生成面に生成する第2光学系と、
    前記第1空間光変調器の前記第1強度変調分布および前記第2空間光変調器の前記第2強度変調分布の双方または何れか一方を設定する制御部と、
    を備え、
    前記制御部が、生成すべき擬似スペックルパターンの光強度統計分布に応じた統計分布を有する擬似乱数パターンに基づいて前記第1強度変調分布を前記第1空間光変調器に設定する、
    擬似スペックルパターン生成装置。
  2. 光を出力する光源と、
    擬似乱数パターンに基づく第1強度変調分布を有し、前記光源から出力された光を前記第1強度変調分布に応じて空間的に強度変調して、その変調後の光を出力する第1空間光変調器と、
    前記第1空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを生成する第1光学系と、
    前記第1光学系により前記フーリエ変換パターンが生成される面に設けられ、フィルタ関数に基づく第2強度変調分布を有し、前記第1光学系を経て到達した光を前記第2強度変調分布に応じて空間的に強度変調して、その変調後の光を出力する第2空間光変調器と、
    前記第2空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを擬似スペックルパターンとして光パターン生成面に生成する第2光学系と、
    前記第1空間光変調器の前記第1強度変調分布および前記第2空間光変調器の前記第2強度変調分布の双方または何れか一方を設定する制御部と、
    を備え、
    前記制御部が、生成すべき擬似スペックルパターンの自己相関関数に応じた相関関数のフーリエ変換の平方根であるフィルタ関数に基づいて前記第2強度変調分布を前記第2空間光変調器に設定する、
    擬似スペックルパターン生成装置。
  3. 光を出力する光源と、
    擬似乱数パターンに基づく第1強度変調分布を有し、前記光源から出力された光を前記第1強度変調分布に応じて空間的に強度変調して、その変調後の光を出力する第1空間光変調器と、
    前記第1空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを生成する第1光学系と、
    前記第1光学系により前記フーリエ変換パターンが生成される面に設けられ、フィルタ関数に基づく第2強度変調分布を有し、前記第1光学系を経て到達した光を前記第2強度変調分布に応じて空間的に強度変調して、その変調後の光を出力する第2空間光変調器と、
    前記第2空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを擬似スペックルパターンとして光パターン生成面に生成する第2光学系と、
    前記第1空間光変調器の前記第1強度変調分布および前記第2空間光変調器の前記第2強度変調分布の双方または何れか一方を設定する制御部と、
    を備え、
    前記制御部が、生成すべき擬似スペックルパターンの光強度統計分布に応じた統計分布を有する擬似乱数パターンに基づいて前記第1強度変調分布を前記第1空間光変調器に設定し、生成すべき擬似スペックルパターンの自己相関関数に応じた相関関数のフーリエ変換の平方根であるフィルタ関数に基づいて前記第2強度変調分布を前記第2空間光変調器に設定する、
    擬似スペックルパターン生成装置。
  4. 請求項1〜3の何れか1項に記載の擬似スペックルパターン生成装置と、
    前記擬似スペックルパターン生成装置により前記擬似スペックルパターンが生成される前記光パターン生成面に照射される観察用の光を出力する観察用光源と、
    前記光パターン生成面への前記観察用の光の照射に応じて生じた光を受光して撮像するカメラと、
    を備える観察装置。
  5. 擬似乱数パターンに基づく第1強度変調分布を有する第1空間光変調器により、光源から出力された光を前記第1強度変調分布に応じて空間的に強度変調して、その変調後の光を出力し、
    第1光学系により、前記第1空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを生成し、
    前記第1光学系により前記フーリエ変換パターンが生成される面に設けられフィルタ関数に基づく第2強度変調分布を有する第2空間光変調器により、前記第1光学系を経て到達した光を前記第2強度変調分布に応じて空間的に強度変調して、その変調後の光を出力し、
    第2光学系により、前記第2空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを擬似スペックルパターンとして光パターン生成面に生成し、
    生成すべき擬似スペックルパターンの光強度統計分布に応じた統計分布を有する擬似乱数パターンに基づいて前記第1空間光変調器の前記第1強度変調分布を設定する、
    擬似スペックルパターン生成方法。
  6. 擬似乱数パターンに基づく第1強度変調分布を有する第1空間光変調器により、光源から出力された光を前記第1強度変調分布に応じて空間的に強度変調して、その変調後の光を出力し、
    第1光学系により、前記第1空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを生成し、
    前記第1光学系により前記フーリエ変換パターンが生成される面に設けられフィルタ関数に基づく第2強度変調分布を有する第2空間光変調器により、前記第1光学系を経て到達した光を前記第2強度変調分布に応じて空間的に強度変調して、その変調後の光を出力し、
    第2光学系により、前記第2空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを擬似スペックルパターンとして光パターン生成面に生成し、
    生成すべき擬似スペックルパターンの自己相関関数に応じた相関関数のフーリエ変換の平方根であるフィルタ関数に基づいて前記第2空間光変調器の前記第2強度変調分布を設定する、
    擬似スペックルパターン生成方法。
  7. 擬似乱数パターンに基づく第1強度変調分布を有する第1空間光変調器により、光源から出力された光を前記第1強度変調分布に応じて空間的に強度変調して、その変調後の光を出力し、
    第1光学系により、前記第1空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを生成し、
    前記第1光学系により前記フーリエ変換パターンが生成される面に設けられフィルタ関数に基づく第2強度変調分布を有する第2空間光変調器により、前記第1光学系を経て到達した光を前記第2強度変調分布に応じて空間的に強度変調して、その変調後の光を出力し、
    第2光学系により、前記第2空間光変調器から出力された光のパターンを光学的にフーリエ変換して当該フーリエ変換パターンを擬似スペックルパターンとして光パターン生成面に生成し、
    生成すべき擬似スペックルパターンの光強度統計分布に応じた統計分布を有する擬似乱数パターンに基づいて前記第1空間光変調器の前記第1強度変調分布を設定し、生成すべき擬似スペックルパターンの自己相関関数に応じた相関関数のフーリエ変換の平方根であるフィルタ関数に基づいて前記第2空間光変調器の前記第2強度変調分布を設定する、
    擬似スペックルパターン生成方法。
  8. 請求項5〜7の何れか1項に記載の擬似スペックルパターン生成方法により前記擬似スペックルパターンが生成される前記光パターン生成面に、観察用光源から出力される観察用の光を照射し、
    前記光パターン生成面への前記観察用の光の照射に応じて生じた光をカメラにより受光して撮像する、
    観察方法。
JP2017077474A 2017-04-10 2017-04-10 擬似スペックルパターン生成装置、擬似スペックルパターン生成方法、観察装置および観察方法 Active JP6739391B2 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2017077474A JP6739391B2 (ja) 2017-04-10 2017-04-10 擬似スペックルパターン生成装置、擬似スペックルパターン生成方法、観察装置および観察方法
PCT/JP2018/012130 WO2018190111A1 (ja) 2017-04-10 2018-03-26 擬似スペックルパターン生成装置、擬似スペックルパターン生成方法、観察装置および観察方法
DE112018001943.8T DE112018001943T5 (de) 2017-04-10 2018-03-26 Pseudo-specklemuster-erzeugungsvorrichtung, pseudo-specklemuster-erzeugungsverfahren, beobachtungsvorrichtung und beobachtungsverfahren
US16/499,924 US11300800B2 (en) 2017-04-10 2018-03-26 Pseudo speckle pattern generation device, pseudo speckle pattern generation method, observation device, and observation method
CN201880023619.0A CN110494794B (zh) 2017-04-10 2018-03-26 伪散斑图案生成装置、伪散斑图案生成方法、观察装置和观察方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017077474A JP6739391B2 (ja) 2017-04-10 2017-04-10 擬似スペックルパターン生成装置、擬似スペックルパターン生成方法、観察装置および観察方法

Publications (2)

Publication Number Publication Date
JP2018180177A JP2018180177A (ja) 2018-11-15
JP6739391B2 true JP6739391B2 (ja) 2020-08-12

Family

ID=63792641

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017077474A Active JP6739391B2 (ja) 2017-04-10 2017-04-10 擬似スペックルパターン生成装置、擬似スペックルパターン生成方法、観察装置および観察方法

Country Status (5)

Country Link
US (1) US11300800B2 (ja)
JP (1) JP6739391B2 (ja)
CN (1) CN110494794B (ja)
DE (1) DE112018001943T5 (ja)
WO (1) WO2018190111A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL273248B2 (en) * 2019-03-19 2023-05-01 Continuse Biometrics Ltd System and method for optical monitoring based on symmetry conditions of the optical field
CN111045230B (zh) * 2019-12-30 2023-01-03 上海天马微电子有限公司 空间光调制器及其制作方法及全息3d显示设备
CN112804513B (zh) * 2021-01-05 2023-02-17 暨南大学 一种光场相机及成像方法
CN113281305B (zh) * 2021-05-17 2023-04-14 太原理工大学 一种基于散射介质实现超分辨显微成像方法及其装置

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3475947B2 (ja) * 1991-05-21 2003-12-10 セイコーエプソン株式会社 光学装置
JP3264544B2 (ja) * 1993-02-16 2002-03-11 浜松ホトニクス株式会社 光学的変位量測定装置
GB0223119D0 (en) * 2002-10-05 2002-11-13 Holographic Imaging Llc Reconfigurable spatial light modulators
KR100981401B1 (ko) * 2004-04-22 2010-09-10 고쿠리쓰다이가쿠호진 덴키쓰신다이가쿠 미소 변위 계측법 및 장치
JP5001286B2 (ja) * 2005-10-11 2012-08-15 プライム センス リミティド 対象物再構成方法およびシステム
GB2456170B (en) * 2008-01-07 2012-11-21 Light Blue Optics Ltd Holographic image display systems

Also Published As

Publication number Publication date
CN110494794B (zh) 2021-07-27
US11300800B2 (en) 2022-04-12
JP2018180177A (ja) 2018-11-15
WO2018190111A1 (ja) 2018-10-18
US20200117018A1 (en) 2020-04-16
CN110494794A (zh) 2019-11-22
DE112018001943T5 (de) 2020-04-09

Similar Documents

Publication Publication Date Title
JP6739391B2 (ja) 擬似スペックルパターン生成装置、擬似スペックルパターン生成方法、観察装置および観察方法
JP4820750B2 (ja) 複数の位相コントラストフィルタによる所望波面の生成
US9599805B2 (en) Optical imaging system using structured illumination
JP6377247B2 (ja) 物体を画像化する方法及び装置
JP2014032371A (ja) スペックルコントラスト発生器及びスペックルコントラスト評価装置
US20230101292A1 (en) Pseudo speckle pattern generation device, pseudo speckle pattern generation method, observation device, and observation method
JP2018180179A (ja) 量子シミュレータおよび量子シミュレーション方法
KR101492343B1 (ko) 광 산란 패턴 검출 장치 및 광 산란 패턴 검출 방법
US11867625B2 (en) System and method for imaging via scattering medium
US20210255488A1 (en) Wavefront Shaping With A 1D Spatial Light Modulator
JP2021103857A (ja) 撮像装置および撮像方法
Kirby et al. Assessment of incident intensity on laser speckle contrast imaging using a nematic liquid crystal spatial light modulator
Ostrovsky et al. Modulation of spatial coherence of optical field by means of liquid crystal light modulator
JP7412166B2 (ja) 撮像装置および撮像方法
CN113406004A (zh) 电磁波相位振幅生成装置、方法和非临时记录介质
Liu et al. Experimental implementation of wavefront sensorless real-time adaptive optics aberration correction control loop with a neural network
Cassina et al. Speckled-speckle field as a resource for imaging techniques
Abregana et al. Digital micromirror device as amplitude diffuser for multiple-plane phase retrieval
Bañas et al. GPC: Recent developments
CN115981018A (zh) 一种非瑞利散斑生成及演化的装置和方法
Voelz et al. Wave Optics Modeling and Laboratory Generation of “Exotic” Partial Coherent Beams
Herrera-Ramirez et al. Evaluating polarization diversity for speckle reduction in retinal image quality estimates
Stoykova et al. SLM-based Testing of Algorithms in Dynamic Speckle Metrology
Ji et al. Speckle Pattern Calculation and Scattering Imaging for General Surface Scattering Problem
Feldkhun Doppler encoded excitation patterning (deep) microscopy

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20191204

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200519

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200701

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200714

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200721

R150 Certificate of patent or registration of utility model

Ref document number: 6739391

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150