JP2018048954A - 撮像装置および撮像方法 - Google Patents

撮像装置および撮像方法 Download PDF

Info

Publication number
JP2018048954A
JP2018048954A JP2016185552A JP2016185552A JP2018048954A JP 2018048954 A JP2018048954 A JP 2018048954A JP 2016185552 A JP2016185552 A JP 2016185552A JP 2016185552 A JP2016185552 A JP 2016185552A JP 2018048954 A JP2018048954 A JP 2018048954A
Authority
JP
Japan
Prior art keywords
light
imaged
imaging
container
depth
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.)
Granted
Application number
JP2016185552A
Other languages
English (en)
Other versions
JP6728007B2 (ja
Inventor
佳祐 加藤
Keisuke Kato
佳祐 加藤
石川 直樹
Naoki Ishikawa
直樹 石川
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.)
Screen Holdings Co Ltd
Original Assignee
Screen Holdings Co Ltd
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 Screen Holdings Co Ltd filed Critical Screen Holdings Co Ltd
Priority to JP2016185552A priority Critical patent/JP6728007B2/ja
Priority to EP17191736.2A priority patent/EP3299761B1/en
Priority to US15/711,684 priority patent/US10317188B2/en
Publication of JP2018048954A publication Critical patent/JP2018048954A/ja
Application granted granted Critical
Publication of JP6728007B2 publication Critical patent/JP6728007B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02055Reduction or prevention of errors; Testing; Calibration
    • G01B9/02062Active error reduction, i.e. varying with time
    • G01B9/02067Active error reduction, i.e. varying with time by electronic control systems, i.e. using feedback acting on optics or light
    • G01B9/02069Synchronization of light source or manipulator and detector
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02055Reduction or prevention of errors; Testing; Calibration
    • G01B9/02062Active error reduction, i.e. varying with time
    • G01B9/02063Active error reduction, i.e. varying with time by particular alignment of focus position, e.g. dynamic focussing in optical coherence tomography
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02055Reduction or prevention of errors; Testing; Calibration
    • G01B9/02062Active error reduction, i.e. varying with time
    • G01B9/02064Active error reduction, i.e. varying with time by particular adjustment of coherence gate, i.e. adjusting position of zero path difference in low coherence interferometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02083Interferometers characterised by particular signal processing and presentation
    • G01B9/02087Combining two or more images of the same region
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/0209Low-coherence interferometers
    • G01B9/02091Tomographic interferometers, e.g. based on optical coherence

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Signal Processing (AREA)
  • Automation & Control Theory (AREA)
  • Optics & Photonics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

【課題】簡単な構成で、容器壁面での反射に起因する画像ノイズの影響のない画像品質の良好な断層画像を得る。
【解決手段】光透過性の壁部を有する容器内に収容された被撮像物を断層撮像するFD−OCT撮像装置において、深さ方向の反射光強度分布における物体光学系の焦点深さに対応する位置が、容器壁面のうち被撮像物側の第1表面に共役な位置と、これと反対側の第2表面に共役な位置との間となるように、焦点深さの設定に連動させて参照光の光路長を変更設定する。
【選択図】図6

Description

この発明は、被撮像物からの反射光と参照光との干渉光成分を検出して撮像を行う技術に関し、特に光透過性を有する容器内の被撮像物を撮像する技術に関するものである。
医学や生化学の技術分野では、容器中で培養された細胞や微生物を観察することが行われる。観察対象となる細胞等に影響を与えることなく観察を行う方法として、顕微鏡等を用いて細胞等を撮像する技術が提案されている。このような技術の1つとして、光コヒーレンストモグラフィ技術を利用したものがある。この技術は、光源から出射される低コヒーレンス光を照明光として被撮像物に入射させ、被撮像物からの反射光(信号光)と光路長が既知である参照光との干渉光を検出することで、被撮像物からの反射光の深さ方向における強度分布を求めて断層画像化するものである。
このうち、広帯域成分を含む信号光と参照光とを干渉させ、干渉光のスペクトルをフーリエ変換することで深さ方向の反射光強度分布を得る技術は、フーリエドメイン光干渉断層撮像技術(Fourier Domain Optical Coherence Tomography;FD−OCT)と呼ばれる。FD−OCT技術では、干渉光のスペクトルにおける波長が被撮像物の深さ方向に対応する情報となるため、光学系の被写界深度内の被撮像物の断層情報を一括して取得することができ、高速撮像が可能である。
FD−OCT技術においては、画像に複素共役ノイズと呼ばれるノイズが被撮像物の像に重畳されることが知られている。これは、フーリエ変換によって干渉光のスペクトルを深さ方向の反射光強度分布に変換するという原理上、不可避的に発生するものである。このようなノイズの影響を排除するための技術としては、例えば特許文献1に記載の技術がある。この技術は、複素共役ノイズをキャンセルするための従来技術である位相シフト法が有する複数回の撮像を要するという問題を解決するため、光路長の異なる2つの参照光をチョッパで切り替えながら信号光と干渉させるようにしたものである。
特開2010−164574号公報
容器中に担持された細胞等を被撮像物とする場合、光透過性を有する容器の壁部(例えば底部)を介して撮像が行われる場合がある。このような場合、容器壁面が強い反射面として作用するため、断層画像には容器壁面からの反射光に起因する複素共役ノイズが重畳されることになる。参照基準面の設定によっては、必要とされる被撮像物の像と容器壁面の複素共役像とが重なってしまうことがあり、被撮像物の正確な断層情報が得られないという問題が生じ得る。
上記従来技術のノイズ除去方法は、このようなケースにも一定の効果を奏し得ると考えられる。しかしながら、容器壁部からの強い反射がある場合を想定したものではないため、実効性は不明である。また、複数の参照系を有する特殊な装置構成が必要となり、装置が複雑化、大型化して高コストになるという問題がある。このことから、このような複雑な装置構成と必要とせずに、容器壁面からの反射光に起因する複素共役ノイズの影響を抑制することのできる技術の確立が望まれる。
この発明は上記課題に鑑みなされたものであり、被撮像物からの反射光と参照光との干渉を利用して容器内の被撮像物を撮像する技術において、簡単な構成で、容器壁面での反射に起因する画像ノイズの影響のない画像品質の良好な断層画像を得ることのできる技術を提供することを目的とする。
この発明の一の態様は、光透過性の壁部を有する容器内に収容された被撮像物を断層撮像する撮像装置において、上記目的を達成するため、光源から出射される広帯域の低コヒーレンス光が分岐された一の分岐光を照明光として前記壁部を介して前記被撮像物に入射させ、前記壁部を介して出射される前記被撮像物の反射光を物体光学系により集光した信号光と、他の一の分岐光から生成された参照光とが干渉して生じる干渉光を検出し、検出された前記干渉光に応じた干渉信号を出力する検出手段と、前記干渉信号に基づき、前記干渉光のスペクトルをフーリエ変換して前記照明光の入射方向における前記被撮像物の反射光強度分布を求め、該反射光強度分布から断層画像を作成する信号処理手段と、前記被撮像物に対する前記物体光学系の焦点深さを変更設定する制御手段とを備え、前記制御手段は、深さ方向の位置ごとの反射光強度を表す前記反射光強度分布において、前記焦点深さに対応する位置が、前記壁部の表面のうち前記被撮像物側の第1表面に共役な位置と、前記被撮像物とは反対側の第2表面に共役な位置との間となるように、前記焦点深さの設定に連動させて前記参照光の光路長を変更設定する。
また、この発明の他の態様は、光透過性の壁部を有する容器内に収容された被撮像物を断層撮像する撮像方法において、上記目的を達成するため、光源から出射される広帯域の低コヒーレンス光が分岐された一の分岐光を照明光として前記壁部を介して前記被撮像物に入射させ、前記壁部を介して出射される前記被撮像物の反射光を物体光学系により集光した信号光と、他の一の分岐光から生成された参照光とが干渉して生じる干渉光を検出し、検出された前記干渉光に応じた干渉信号を出力する工程と、前記干渉信号に基づき、前記干渉光のスペクトルをフーリエ変換して前記照明光の入射方向における前記被撮像物の反射光強度分布を求め、該反射光強度分布から断層画像を作成する工程とを備え、前記被撮像物に対する前記物体光学系の焦点深さを変更設定可能であり、深さ方向の位置ごとの反射光強度を表す前記反射光強度分布において、前記焦点深さに対応する位置が、前記壁部の表面のうち前記被撮像物側の第1表面に共役な位置と、前記被撮像物とは反対側の第2表面に共役な位置との間となるように、前記焦点深さの設定に連動させて前記参照光の光路長を変更設定する。
この発明における撮像原理は、広帯域成分を含む信号光と参照光とを干渉させ、干渉光のスペクトルをフーリエ変換することで深さ方向の反射光強度分布を得る、FD−OCT技術によるものである。容器壁部を介したFD−OCT撮像では、壁部の表面のうち被撮像物に近い側の第1表面、第1表面を挟んで被撮像物とは反対側の第2表面それぞれの複素共役像が被撮像物の像に重なってノイズとなり得る。
容器壁面の複素共役像は、容器壁面の像が現れる位置に対して共役な、すなわち当該位置に対し参照光の光路長で規定される参照基準面に関して対称な位置に現れるから、その出現位置は参照基準面の設定に依存する。このため、参照基準面を被撮像物から十分に離すことで容器壁面の複素共役像を被撮像物の像から分離するという対応も考えられるが、画像品質の点では参照基準面を被撮像物の近傍に設定することが望ましい。より具体的には、被撮像物からの反射光を集光する物体光学系の焦点面と参照基準面とができるだけ近いことが好ましい。
そこで本発明では、物体光学系の焦点深さに応じて参照基準面の位置が設定される。具体的には、物体光学系の光軸方向に沿った深さ方向の反射光強度分布において、第1表面および第2表面それぞれに共役な位置が物体光学系の焦点深さに対応する位置を挟む位置関係となるように、参照基準面を規定する参照光の光路長が設定される。断層画像に第1表面の複素共役像と第2表面の複素共役像とが現れたとしても、これらで挟まれる領域には容器壁面に起因するさらなる複素共役像は現れない。したがって、容器壁面の複素共役像の間に物体光学系の焦点深さに相当する位置が来るように、言い換えれば物体光学系の焦点面を挟む位置に容器壁面の複素共役像が出現するように参照光の光路長を設定すれば、少なくとも焦点面の近傍では容器壁面に起因する複素共役ノイズの影響を受けずに被撮像物の断層情報を取得することが可能である。
使用される容器の形状から壁部の厚さを予め把握しておくことが可能であり、撮像時において第1表面および第2表面が占める位置も推定可能である。したがって、第1表面および第2表面の複素共役像が上記の位置に出現するような参照光の光路長を、焦点深さの設定に連動させて設定することが可能である。参照基準面を焦点面から大きく離して設定する必要がないので、画像品質も良好なものとすることができる。
上記のように、本発明によれば、物体光学系の焦点深さに相当する位置の近傍において容器壁面に起因する複素共役ノイズの影響を排除することができる。物体光学系の焦点深さの設定に連動させて参照光の光路長を設定する機能があれば上記効果が得られるので、簡単な装置構成で複素共役ノイズの影響を排除した断層画像を取得することができる。
本発明にかかる画像処理装置の構成例を示す原理図である。 この画像処理装置における撮像原理を説明する図である。 OCT装置の具体的構成例を示す図である。 物体光学系の焦点深さと参照基準面との位置関係を模式的に示す図である。 参照基準面の位置と反射光強度分布との関係を示す図である。 この撮像装置における撮像動作の原理を説明する図である。 焦点深さに応じた参照基準面の設定原理を示す図である。 複素共役信号のピークが一定の広がりを有する場合の考え方を示す図である。 この撮像装置における撮像動作を示すフローチャートである。 焦点深さが容器底部厚さの半分未満であるときの反射光ピーク位置を示す図である。
図1は本発明にかかる画像処理装置の構成例を示す原理図である。この画像処理装置1は、培地中で培養されたスフェロイド(細胞集塊)を被撮像物として断層撮像し、得られた断層画像を画像処理して、スフェロイドの立体像を作成する。なお、ここでは培地中のスフェロイドを被撮像物とした例を説明するが、被撮像物はこれに限定されない。以下の各図における方向を統一的に示すために、図1に示すようにXYZ直交座標軸を設定する。ここでXY平面が水平面を表す。また、Z軸が鉛直軸を表し、より詳しくは(−Z)方向が鉛直下向き方向を表している。
画像処理装置1は保持部10を備えている。保持部10は、ガラス製または樹脂製の透明で均質な平底を有する浅皿状のディッシュと呼ばれる容器11を、その開口面を上向きにして略水平姿勢に保持する。容器11には予め適宜の培地Mが所定量注入されており、培地中では容器11の底部111にスフェロイドSpが培養されている。図1では1つのスフェロイドSpのみが記載されているが、1つの容器11内で複数のスフェロイドSpが培養されていてもよい。
保持部10により保持された容器11の下方に、撮像ユニット20が配置される。撮像ユニット20には、被撮像物の断層画像を非接触、非破壊(非侵襲)で撮像することが可能な光干渉断層撮像(Optical Coherence Tomography;OCT)装置が用いられる。詳しくは後述するが、OCT装置である撮像ユニット20は、被撮像物への照明光を発生する光源21と、ビームスプリッタ22と、物体光学系23と、参照ミラー24と、分光器25と、光検出器26とを備えている。
また、画像処理装置1はさらに、装置の動作を制御する制御ユニット30と、撮像ユニット20の可動機構を制御する駆動制御部40とを備えている。制御ユニット30は、CPU(Central Processing Unit)31、A/Dコンバータ32、信号処理部33、3D復元部34、インターフェース(IF)部35、画像メモリ36およびメモリ37を備えている。
CPU31は、所定の制御プログラムを実行することで装置全体の動作を司り、CPU31が実行する制御プログラムや処理中に生成したデータはメモリ37に保存される。A/Dコンバータ32は、撮像ユニット20の光検出器26から受光光量に応じて出力される信号をデジタルデータに変換する。信号処理部33は、A/Dコンバータ32から出力されるデジタルデータに基づき後述する信号処理を行って、被撮像物の断層画像を作成する。3D復元部34は、撮像された複数の断層画像の画像データに基づいて、撮像された細胞集塊の立体像(3D像)を作成する機能を有する。信号処理部33により作成された断層画像の画像データおよび3D復元部34により作成された立体像の画像データは、画像メモリ36により適宜記憶保存される。
インターフェース部35は画像処理装置1と外部との通信を担う。具体的には、インターフェース部35は、外部機器と通信を行うための通信機能と、ユーザからの操作入力を受け付け、また各種の情報をユーザに報知するためのユーザインターフェース機能とを有する。この目的のために、インターフェース部35には、装置の機能選択や動作条件設定などに関する操作入力を受け付け可能な例えばキーボード、マウス、タッチパネルなどの入力デバイス351と、信号処理部33により作成された断層画像や3D復元部34により作成された立体像など各種の処理結果を表示する例えば液晶ディスプレイからなる表示部352とが接続されている。
また、CPU31は駆動制御部40に制御指令を与え、これに応じて駆動制御部40は撮像ユニット20の可動機構に所定の動作を行わせる。次に説明するように、駆動制御部40により実行される撮像ユニット20の走査移動と、光検出器26による受光光量の検出との組み合わせにより、被撮像物であるスフェロイド(細胞集塊)の断層画像が取得される。
図2はこの画像処理装置における撮像原理を説明する図である。より具体的には、図2(a)は撮像ユニット20における光路を示す図であり、図2(b)はスフェロイドの断層撮像の様子を模式的に示す図である。前記したように、撮像ユニット20は光干渉断層撮像(OCT)装置として機能するものである。
撮像ユニット20では、例えば発光ダイオードまたはスーパールミネッセントダイオード(SLD)などの発光素子を有する光源21から、広帯域の波長成分を含む低コヒーレンス光ビームL1が出射される。光ビームL1はビームスプリッタ22に入射して分岐し、破線矢印で示すように一部の光L2が容器11に向かい、一点鎖線矢印で示すように一部の光L3が参照ミラー24に向かう。
容器11に向かった光L2は、物体光学系23を経て容器11に入射する。より具体的には、ビームスプリッタ22から出射される光L2は、物体光学系23を介して容器底部111に入射する。物体光学系23は、ビームスプリッタ22から容器11に向かう光L2を容器11内の被撮像物(この場合にはスフェロイドSp)に収束させる機能と、被撮像物から出射される反射光を集光してビームスプリッタ22に向かわせる機能とを有する。図では物体光学系23は単一の対物レンズにより代表的に表されているが、複数の光学素子が組み合わされたものであってもよい。
物体光学系23は、駆動制御部40に設けられた焦点調整機構41により、Z方向に移動可能に支持されている。これにより、被撮像物に対する物体光学系23の焦点位置がZ方向に変更可能となっている。以下、深さ方向(Z方向)における物体光学系23の焦点位置を「焦点深さ」と称する。物体光学系23の光軸は鉛直方向と平行であり、したがって平面状の容器底部111に垂直である。また、物体光学系23への照明光の入射方向は光軸と平行であり、その光中心が光軸と一致するように、物体光学系23の配置が定められている。
スフェロイドSpが光L2に対する透過性を有するものでなければ、容器底部111を介して入射した光L2はスフェロイドSpの表面で反射される。一方、スフェロイドSpが光L2に対してある程度の透過性を有するものである場合、光L2はスフェロイドSp内まで進入してその内部の構造物により反射される。光L2として例えば近赤外線を用いることで、入射光をスフェロイドSp内部まで到達させることが可能である。スフェロイドSpからの反射光は散乱光として種々の方向に放射される。そのうち物体光学系23の集光範囲内に放射された光L4が、物体光学系23で集光されてビームスプリッタ22へ送られる。
参照ミラー24は、駆動制御部40に設けられたミラー駆動機構42により、その反射面を光L3の入射方向に対し垂直姿勢に、しかも、該入射方向に沿った方向(図ではY方向)に移動可能に支持されている。参照ミラー24に入射した光L3は反射面で反射されて、入射光路を逆向きに辿るように進む光L5としてビームスプリッタ22に向かう。この光L5が参照光となる。ミラー駆動機構42により参照ミラー24の位置が変更されることにより、参照光の光路長が変化する。
スフェロイドSpの表面もしくは内部の反射面で反射された反射光L4と、参照ミラー24で反射された参照光L5とは、ビームスプリッタ22を介して光検出器26に入射する。このとき、反射光L4と参照光L5との間で位相差に起因する干渉が生じるが、干渉光の分光スペクトルは反射面の深さにより異なる。つまり、干渉光の分光スペクトルは被撮像物の深さ方向の情報を有している。したがって、干渉光を波長ごとに分光して光量を検出し、検出された干渉信号をフーリエ変換することにより、被撮像物の深さ方向における反射光強度分布を求めることができる。このような原理に基づくOCT撮像技術は、フーリエドメイン(Fourier Domain)OCT(FD−OCT)と称される。
この実施形態の撮像ユニット20は、ビームスプリッタ22から光検出器26に至る干渉光の光路上に分光器25が設けられている。分光器25としては、例えばプリズムを利用したもの、回折格子を利用したもの等を用いることができる。干渉光は分光器25により波長成分ごとに分光されて光検出器26に受光される。
光検出器26が検出した干渉光に応じて光検出器26から出力される干渉信号をフーリエ変換することで、スフェロイドSpのうち、光ビームL2の入射位置における深さ方向、つまりZ方向の反射光強度分布が求められる。容器11に入射する光ビームL2をX方向に走査することで、XZ平面と平行な平面における反射光強度分布が求められ、その結果から当該平面を断面とするスフェロイドSpの断層画像を作成することができる。以下、本明細書では、X方向へのビーム走査によってXZ平面と平行な断面における1つの断層画像Itを取得する一連の動作を、1回の撮像と称することとする。
また、Y方向におけるビーム入射位置を多段階に変更しながら、その都度断層画像の撮像を行うことで、図2(b)に示すように、スフェロイドSpをXZ平面と平行な断面で断層撮像した多数の断層画像Itを得ることができる。Y方向の走査ピッチを小さくすれば、スフェロイドSpの立体構造を把握するのに十分な分解能の画像データを得ることができる。X方向およびY方向へのビーム走査は、例えば図示しないガルバノミラー等の光路を変化させる光学部品を用いてビーム入射位置をXY方向に変化させる方法、スフェロイドSpを担持する容器11と撮像ユニット20とのいずれかをXY方向に移動させてこれらの相対位置を変化させる方法などにより実現可能である。
なお、上記の原理説明では、撮像ユニット20において光源21からの光を照明光と参照光とに分岐させる分波機能、および信号光と参照光とを合成して干渉光を生じさせる機能がビームスプリッタ22により実現されている。一方、近年では、OCT装置においてこのような分波・合波機能を担うものとして、以下に例示するような光ファイバカプラが用いられる場合がある。
図3はOCT装置の具体的構成例を示す図である。なお、理解を容易にするために、以下の説明では、上記した原理図の構成と同一のまたは相当する構成に同一符号を付すものとする。その構造および機能は、特に説明のない限り上記原理図のものと基本的に同じであり、詳しい説明は省略する。また、光ファイバカプラによる干渉光を検出するOCT撮像原理も基本的に上記と同じであるので、詳しい説明を省略する。
図3(a)に示す構成例では、撮像ユニット20aは、ビームスプリッタ22に代わる分波・合波器として光ファイバカプラ220を備えている。光ファイバカプラ220を構成する光ファイバの1つ221は光源21に接続されており、光源21から出射される低コヒーレンス光は、光ファイバカプラ220により2つの光ファイバ222,223への光に分岐される。光ファイバ222は物体系光路を構成する。より具体的には、光ファイバ222の端部から出射される光はコリメータレンズ223を介して物体光学系23に入射する。被撮像物からの反射光(信号光)は物体光学系23、コリメータレンズ223を介して光ファイバ222に入射する。
他の光ファイバ224は参照系光路を構成する。より具体的には、光ファイバ224の端部から出射される光はコリメータレンズ225を介して参照ミラー24に入射する。参照ミラー24からの反射光(参照光)はコリメータレンズ225を介して光ファイバ224に入射する。光ファイバ222を伝搬する信号光と光ファイバ224を伝搬する参照光とが光ファイバカプラ220において干渉し、干渉光が光ファイバ226および分光器25を介して光検出器26に入射する。光検出器26により受光された干渉光から被撮像物における反射光の強度分布が求められることは上記原理通りである。
図3(b)に示す例でも、撮像ユニット20bに光ファイバカプラ220が設けられる。ただし光ファイバ224は使用されず、光ファイバ222から出射される光の光路に対してコリメータレンズ223およびビームスプリッタ227が設けられる。そして、前述の原理通り、ビームスプリッタ227により分岐される2つの光路にそれぞれ物体光学系23、参照ミラー24が配置される。このような構成ではビームスプリッタ227により信号光と参照光とが合成され、それにより生じた干渉光が光ファイバ222,226を通って光検出器26へ導かれる。
これらの例では、図1の原理図では空間中を進行する各光の光路の一部が光ファイバに置き換えられているが動作原理は同じである。これらの例においても、焦点調整機構41が物体光学系23を容器11に対し接近・離間方向に移動させることにより、被撮像物に対する物体光学系23の焦点深さを調整することが可能である。また、ミラー駆動機構42が参照ミラー24を光の入射方向に沿って移動させることにより、参照光の光路長を変更可能である。
以下、この画像処理装置1による撮像動作について説明する。撮像ユニットの構成は、上記したビームスプリッタを用いるもの、光ファイバカプラを用いるもののいずれであっても同じ撮像動作が可能である。
図4は物体光学系の焦点深さと参照基準面との位置関係を模式的に示す図である。また、図5は参照基準面の位置と反射光強度分布との関係を示す図である。OCT撮像装置では、その原理上、信号光の光路長が参照光の光路長と同じになる位置が、画像における深さ方向の基準位置となる。
以下の説明では、図4に示すように、物体光学系23を介して照明光L2および信号光L4が伝搬する物体系光路において、参照系光路における参照ミラー24の反射面と対応する位置を参照基準面Srと称する。また、容器11の底部111のうち被撮像物であるスフェロイドSpに近い側の表面(つまり培地と接触する内底面)を上部底面Sa、容器11の底部111のうち上部底面Saとは反対側の外底面を下部底面Sbと称する。さらに、物体光学系23の焦点面、つまり物体光学系23の物体側焦点FPを含み物体光学系23の光軸AXに垂直な平面を符号Sfにより表す。
また、後の説明のために、上部底面Saと下部底面Sbとの距離、つまり容器底部111の厚さを符号Tにより表す。また、上部底面Saと焦点面Sfとの距離、つまり容器11の内底面から焦点FPまでの距離を符号Dにより表す。この距離Dは、容器11の内底面(上部底面Sa)を起点としたときの物体光学系23の焦点深さであるとも言える。
物体系光路において参照基準面Srに仮想的な反射面があるとき、当該反射面で反射される光の(光源21から光検出器26までの)光路長は、参照ミラー24の反射面で反射される光の(光源21から光検出器26までの)光路長と等しくなる。被撮像物近傍にある各反射面の深さ方向位置は、参照基準面SrからのZ方向距離により表される。
被撮像物(スフェロイドSp)が焦点面Sfに反射面を有するとき、フーリエ変換後の反射光強度分布においては、当該反射面からの反射光強度に対応する大きさの信号が、参照基準面Srから当該反射面(すなわち焦点面Sf)までの距離に対応する深さの位置に現れる。実際の被撮像物においては、種々の深さにある反射面からの反射光に対応する信号が各位置に現れ、反射光強度分布はそれらの信号が重畳されたものなる。ただし、ここでは理解を容易にするために、焦点面Sfにある反射面からの信号のみを考えることとする。
容器底部111の平坦な表面Sa、Sbも強い反射面であり、反射光強度分布においては、これらに対応する信号が参照基準面Srからの距離に応じた位置にそれぞれ現れる。例えば図5(a)に示すように、反射光強度分布においては焦点面Sfに対応する信号P1f、上部底面Saに対応する信号P1a、下部底面Sbに対応する信号P1bが、それぞれ参照基準面Srからの距離に応じた位置に現れる。これらの信号は、実空間において各反射面が離隔していれば反射光強度分布においても異なる位置に現れるため、画像処理装置1の分解能が十分に高ければ相互に分離が可能である。
一方、FD−OCT撮像技術において得られる反射光強度分布では、参照基準面Srに関して各反射面と対称な位置に、各反射面からの信号と複素共役な信号が原理的に現れる。すなわち、焦点面Sfに対応する信号P1fと共役な位置にその複素共役信号P2fが、上部底面Saに対応する信号P1aと共役な位置にその複素共役信号P2aが、下部底面Sbに対応する信号P1bと共役な位置にその複素共役信号P2bが、それぞれ現れる。以下では、実空間に存在する反射面からの信号を、複素共役信号と区別するために「実信号」ということがある。
図5(a)に示すように、参照基準面Srが容器11の下部底面Sbよりも下方にあるとき、複素共役信号P2f,P2a,P2bは計算上の仮想空間である(−Z)側に現れる。そのため、(+Z)側で検出される実空間における反射面からの反射光強度には影響を及ぼさない。しかしながら、図5(b)および図5(c)に示すように、参照基準面Srの設定によっては(−Z)側に現れる信号の複素共役信号が(+Z)側の共役位置に現れることがある。
例えば図5(b)に示す例では、焦点面Sfと上部底面Saとの間に参照基準面Srがあり、上部底面Saおよび下部底面Sbに対応する複素共役信号P2a,P2bが(+Z)側に現れている。この場合でも、焦点面Sfに対応する実信号P1fは他の複素共役信号からは離隔しているので検出可能である。しかしながら、参照基準面Sfの設定が僅かに異なる図5(c)に示す例では、焦点面Sfに対応する実信号P1fと上部底面Saに対応する複素共役信号P2aとが(+Z)領域で重なってしまっており、焦点面Sfに対応する信号P1fを検出することが不可能となっている。
このように、原理的には深さ方向の広い範囲を撮像可能なFD−OCT技術であっても、強い反射面として作用する容器底面を介した撮像においては、底面に対応する複素共役信号が被撮像物からの実信号と重なることで、実信号の検出が困難になる場合があり得る。特に底面が強い反射面として作用するとき、単なる画像ノイズのレベルを超えて、実信号自体を遮蔽してしまうことが起こり得る。したがって、容器底面に起因する複素共役ノイズの影響は、被撮像物自身に起因する複素共役ノイズよりも遥かに大きなものとなることがある。
その一方で、形状やサイズ、内部構造等が不特定であるスフェロイドSp等の被撮像物とは異なり、被撮像物を担持する容器11については、その底面の形状や寸法を予め把握することが可能である。すなわち、容器底部111において強い反射面となるのは上部底面Saと下部底面Sbとであり、両者の間に強い反射面は存在しない。したがって、現れる複素共役信号は上部底面Saおよび下部底面Sbに起因するもののみであり、それらの間には他の複素共役信号は現れない。
また、上部底面Saおよび下部底面Sbが平滑な面であれば、これに起因する複素共役信号P2a,P2bは鋭いピークとなって現れその広がりは小さい。そして、参照基準面Srがある位置に設定されるとき、上部底面Saおよび下部底面Sbに起因する複素共役信号のピークが反射光強度分布においてどの位置に現れるかは事前に推定することが可能である。さらに、それらのピークの位置は参照基準面Srの設定に応じて変化するが、両ピーク間の距離は容器底部111の厚さTにより決まり参照基準面Srの設定に対し不変である。
これらのことから、参照基準面Srを適宜に設定し複素共役信号の出現位置を管理することで、少なくとも撮像対象とする一部の領域については容器底面の影響が現れないようにすることができると期待される。具体的には、次に説明するように、物体光学系23の焦点面Sfおよびその近傍からできるだけ離れた位置に、容器底面に起因する複素共役信号が現れるようにすれば、少なくとも焦点深さおよびその近傍範囲で複素共役ノイズの影響を受けない反射光強度分布を取得することができる。
図6はこの撮像装置における撮像動作の原理を説明する図である。より具体的には、図6(a)はこの画像処理装置1における焦点位置と底面に起因する複素共役信号との位置関係を示す図であり、図6(b)は一の断層画像を複数の部分画像から構成する方法を模式的に示す図である。図6(a)に示すように、深さ方向(Z方向)の反射光強度分布においては、物体光学系23の焦点深さZfの位置に、焦点面Sfからの実信号P1fが現れる。容器11の上部底面Saに対応する複素共役信号P2aのピークおよび下部底面Sbに対応する複素共役信号P2bのピークは、参照基準面Srの設定に応じてZ方向に位置を変えるがピーク間距離は不変である。
図6(a)に示すように、複素共役信号の2つのピークが焦点深さZfに対応する位置を挟むように出現するとき、これらのピークに挟まれる領域Reの内部においては、容器底面に起因する複素共役ノイズの影響を受けることがない。この領域Reは焦点深さZfに対応する位置を含んでおり、深さ方向においては最も良好な分解能が得られる領域でもある。したがって、この領域Re内の反射光強度分布から得られる断層画像では、容器底面に起因する複素共役ノイズの影響がなく、しかも合焦状態で検出されているため画像品質も良好なものとなる。この意味において、この領域Reをここでは「有効領域」と称し、そのZ方向の長さを「有効高さ」と称し符号Zeにより表すこととする。
FD−OCT撮像では原理的には深さ方向に広い範囲を一括して撮像することができる。しかしながら、図6(a)から明らかなように、複素共役ノイズの影響を受けない断層画像を得ることができるのは、そのうち複素共役信号の2つのピークP2a,P2bに挟まれた有効領域Re内のみである。
撮像対象とされる範囲のZ方向高さが有効高さZeよりも大きい場合には、図6(b)に示すように、Z方向に撮像位置を異ならせて撮像された複数の部分画像IpをZ方向に並べて合成することにより、撮像対象範囲の全体をカバーする断層画像Itを作成することが可能である。
被撮像物がスフェロイドSpであるとき、その高さは概ね数μmから数百μm程度である。一方、容器底部111の厚さTは、数百μmないし数mmのものが一般的である。したがって、1回の撮像で得られる断層画像に被撮像物の全体が収まる場合とそうでない場合とがあり得る。1回の撮像で得られる断層画像に被撮像物の全体が収まらない場合には、上記したようにZ方向に撮像位置を異ならせて複数回撮像を行えばよい。
このとき、各部分画像Ipにおいて画像品質を良好なものとするために、各部分画像Ipの撮像ごとに、物体光学系23の焦点深さZfを撮像範囲に合わせて変更することが望ましい。ここで、焦点深さZfの変更に連動させて、参照基準面Srの設定も変更されることがさらに望ましい。容器底面に起因する複素共役ノイズの影響をなくすためには焦点深さZfの変更に伴って有効領域ReもZ方向に位置をシフトさせることが望ましく、そのために複素共役信号のピーク出現位置を焦点位置に連動させる必要があるからである。
言い換えれば、撮像時の焦点深さZfに連動させて参照基準面Srの位置を変更することで、容器底面に起因する複素共役ノイズの影響がなく、しかも合焦状態で撮像された部分画像Ipを取得することができる。そして、焦点深さZfを異ならせた複数の部分画像Ipを合成することで、画像全体において容器底面に起因する複素共役ノイズの影響がなく、しかも各深さに合焦した優れた画像品質の断層画像Itを得ることができる。
次に、上記条件に合致する参照基準面の設定方法について説明する。上記したように、この画像処理装置1の撮像動作においては、物体光学系23の焦点深さの設定に応じて参照基準面が設定される。参照基準面の設定は、参照光の光路長を設定することによってなされる。
図7は焦点深さに応じた参照基準面の設定原理を示す図である。図7左端に示すように、容器11の上部底面Saと下部底面Sbとの距離が容器底部の厚さTとして表され、物体光学系23の焦点FPを含む焦点面Sfと容器11の上部底面Saとの距離が、焦点深さDとして表される。
図7(a)に示されるように、フーリエ変換後の反射光強度分布においては、焦点面Sf、上部底面Saおよび下部底面Sbに対応する実信号P1f,P1a,P1bが、各面の深さ方向位置に応じた位置に現れる。実際にフーリエ変換を行って求められる反射光強度分布には、各信号の複素共役信号が重畳される。ここでは、焦点面Sfに対応する実信号Sfが容器底面に対応する複素共役信号のピーク間に挟まれるための条件の導出を試みる。
焦点面Sfに対応する実信号P1fが容器底面に起因する複素共役信号P2a,P2bに挟まれる条件を求めるために、図7(b)に示すように焦点面Sfに対応する実信号P1fが上部底面Saに起因する複素共役信号P2aと重なる条件と、図7(c)に示すように実信号P1fが下部底面Sbに起因する複素共役信号P2bと重なる条件とを考える。
まず、図7(b)に示すように、焦点面Sfに対応する実信号P1fと上部底面Saに起因する複素共役信号P2aとが重なる条件を考える。ここで、基準参照面SrのZ方向位置を符号Zrで表し、容器11の下部底面Sbから基準参照面Srまでの距離を符号Rにより表すとする。そうすると、図7(b)に示す位置関係から、次式:
R=D/2+T=(D+2T)/2 … (1)
の関係が得られる。値Dは撮像時の焦点深さが決まれば一意に決まる。値Tは容器11が決まれば一意に決まる。
一方、図7(c)に示すように、焦点面Sfに対応する実信号P1fと下部底面Sbに起因する複素共役信号P2bとが重なる条件を考える。このとき、図7(c)に示す位置関係から、次式:
R=(D+T)/2 … (2)
の関係が得られる。
焦点面Sfに対応する実信号Sfが容器底面に対応する複素共役信号のピークP2a,P2bの間に位置するための条件は、値Rが上記式(1)で表される値と式(2)で表される値との間にあることである。したがって、このための条件は次式:
(D+T)/2<R<(D+2T)/2 … (3)
により表される。したがって、値Rが上記式(3)の関係を満たすように、焦点深さに応じて参照基準面Srの位置を設定すればよいことがわかる。
特別な例として、図7(d)に示すように、実信号P1fが2つの複素共役信号P2a,P2bのちょうど中間位置に現れる条件を考えると、図7(d)に示す位置関係から、次式:
R=D/2+3T/4 … (4)
となるように、参照基準面Srの位置が設定されればよい。このような条件が満たされるとき、有効領域Reは焦点深さZfを中心として(+Z)方向および(−Z)方向に同程度の広がりを有することになる。物体光学系23が焦点位置FPを中心として(+Z)方向および(−Z)方向に同程度の合焦範囲(被写界深度内の範囲)を有する場合、上記式(4)の条件を満たすように参照基準面Srの位置を設定することで、物体光学系23の合焦範囲を最も有効に活用して撮像を行うことが可能となる。
なお、上記の考察は複素共役信号のピークP2a,P2bの広がりが無視できる場合に成り立つものである。これらのピークP2a,P2bの広がりが大きいとき、式(3)で表される値Rの範囲のうち上限および下限に近い領域では、これらのピークP2a,P2bが実信号P1fと不可分に重なってしまうことがあり得る。ピークの広がりを考慮した場合の条件は次のようにして求めることができる。
図8は複素共役信号のピークが一定の広がりを有する場合の考え方を示す図である。図8(a)に示すように、複素共役信号のピークP2a,P2bがそれぞれ2Δ程度の広がりを有するものとする。値Δは例えばピークの半値半幅によって特定することができる。図8(a)は図7(b)のケースにピーク幅の影響を加味した場合に相当する。すなわち、図8(a)は、実信号P1fが複素共役信号のピークP2aと完全に一致せず、ピークP2aの半値半幅Δ程度だけもう1つの複素共役信号のピークP1a側にシフトした位置に現れる条件を示している。このようなケースでは、実信号P1fが複素共役信号のピークP2aによって遮蔽されてしまうことが回避される。
このとき、図8(a)に示す位置関係から次式:
R=(D−Δ)/2+T=(D−Δ+2T)/2 … (1a)
が成り立つ。
図示を省略するが、図7(c)のケースについても同様にピーク幅を加味すると、次式:
R=(D+Δ+T)/2 … (2a)
の関係が成り立つ。これらの式(1a)、(2a)から、複素共役信号のピーク幅を考慮したときでも実信号P1fが遮蔽されないための条件は、次式:
(D+Δ+T)/2<R<(D−Δ+2T)/2 … (3a)
となる。値Rが取り得る範囲は式(3)で示される範囲より狭くなるが、複素共役信号のピークの広がりによって実信号P1fが遮蔽されてしまうことはより確実に防止される。
さらに実信号P1fと複素共役信号P2a,P2bとをより確実に分離可能とするために、図8(b)に示されるように、実信号P1fが複素共役信号のピークP2aの半値半幅Δの2倍程度シフトした位置に現れるようにしてもよい。このための条件は次式:
R=(D−2Δ)/2+T=(D−2Δ+2T)/2=D/2−Δ+T
… (1b)
で表される。同様に、図7(c)に対応するケースは次式:
R=(D+2Δ+T)/2 … (2b)
により表される。これらから、値Rが取り得る好ましい範囲は、次式:
(D+2Δ+T)/2<R<(D−2Δ+2T)/2 … (3b)
により表される。
被撮像物から得られる実信号と容器底面に起因する複素共役信号とをどの程度分離させる必要があるかに応じて、上記条件式(3)、(3a)、(3b)を使い分けることができる。具体的には、複素共役信号のピークの幅がほぼ無視できる場合には式(3)を適用することができる。一方、複素共役信号のピークが比較的大きな広がりを有している場合や、実信号のレベルが複素共役信号に対して小さい場合等には、複素共役ノイズの影響をより確実に排除するために、式(3b)が適用されることが好ましい。それらの中間的な場合には式(3a)が適用されればよい。
容器底面に起因する複素共役信号の広がりの程度は、容器底面の表面状態に依存する。すなわち、容器11の上部底面Saおよび下部底面Sbが高精度の平滑面である場合、複素共役信号のピークは鋭くなる。上部底面Saおよび下部底面Sbの表面粗さが大きくなるほどピーク幅が広くなる。したがって、複素共役信号の広がりは容器底面の状態から見積もることが可能である。
なお、上記の条件式(3)、(3a)、(3b)は焦点深さDおよび容器底部厚さTに対応する参照基準面の位置を表す値Rの好ましい範囲を示す式である。しかしながら、容器底部厚さTおよび焦点深さDに応じて参照基準面の位置を導出するに当たって、値Rにこのような幅を持たせることを要するものではない。すなわち、値Rは値D,Tを変数とする関数、すなわち次式:
R=F(D,T) … (5)
を用いて一意に定められればよい。要するに、式(5)に示される関数F(D,T)が返す値Rが、任意の値D,Tに対して条件式(3)、(3a)、(3b)のいずれかで表される範囲に入るように関数F(D,T)が定められていればよい。
この実施形態では、物体光学系23の焦点位置近傍の合焦範囲からの信号を有効に活用するため、上記の式(4)を用いて参照基準面Srの位置が定められるものとする。すなわち、容器底部の厚さTと、容器の上部底面Saを基準としたときの物体光学系23の焦点深さDとが与えられたとき、容器の下部底面Sbを基準としたときの参照基準面の位置Rが式(4)で示される値となるように、参照光の光路長を規定する参照ミラー24の位置が設定される。
図9はこの撮像装置における撮像動作を示すフローチャートである。被撮像物となるスフェロイドSpを担持する容器11が保持部10にセットされると(ステップS101)、CPU31は当該容器11の底部厚さTに関する情報を取得する(ステップS102)。この情報については、ユーザが入力デバイス351を介して入力する態様であってもよく、また予めメモリ37に登録されている容器に関するデータベースから、使用される容器11に応じた情報が読み出される態様であってもよい。このようにしてCPU31は容器11の底部厚さTに関する情報を取得する。
被撮像物の撮像は、焦点深さを多段階に変更設定しながら複数回実行される。すなわち、最初に物体光学系23の焦点深さDが、焦点調整機構41により予め定められた初期値に設定される(ステップS103)。次に、CPU31は、設定された焦点深さDが、容器底部厚さTの半分以上であるか否かを判定する(ステップS104)。
焦点深さDが容器底部厚さTの半分以上であるときには(ステップS104においてYES)、CPU31は上記した式(4)を用いて参照ミラー24の位置を導出する。これに応じて、ミラー駆動部42が参照ミラー24を求められた位置に位置決めする(ステップS105)。一方、焦点深さDが容器底部厚さTの半分より小さいときには(ステップS104においてNO)、CPU31は特殊なケースとして参照基準面の位置に対応する値Rが容器11の底部厚さTと一致するようにミラー駆動部42を制御し、この条件を満たすように参照ミラー24の位置が設定される(ステップS106)。このようにする理由については後述する。
こうして焦点深さおよび参照基準面が設定された状態で、光源21から広帯域の低コヒーレンス光が出射され、被撮像物から出射される信号光と参照光との干渉光が検出される。干渉光は分光器25により分光されて波長成分ごとの強度が光検出器26により検出され、光検出器26の出力信号がスペクトル情報としてCPU31に与えられる。被撮像物への光照射および干渉光の検出は、光入射位置をX方向へ走査しながら複数回実行され、その都度反射光の分光スペクトル情報が取得される(ステップS107)。
信号処理部33は、CPU31から与えられた分光スペクトル情報をフーリエ変換することで、焦点深さ付近に存在する反射面からの反射光強度分布を求める(ステップS108)。前記したように、被撮像物からの反射光強度を有効に求めることができるのは、容器底面に対応する複素共役信号のピークの間の有効領域Re内である。有効領域Reの外側の領域については反射光強度分布の算出対象から除外されてもよく、また反射光強度分布の算出時点では処理対象に含まれた上で後の画像合成時に削除されてもよい。
焦点深さの設定値が所定の最終値に到達するまで(ステップS109)、焦点深さを1ステップずつ段階的に変更しながら(ステップS110)、撮像が繰り返される。焦点深さの変更ステップ幅については、予め定められていてもよく、また容器11の底部厚さTに応じてステップ幅が設定されてもよい。画像全体で合焦した断層画像を得るためには、変更ステップ幅は有効領域Reの高さZe以下であることが望ましい。
なお、例えば物体光学系23が高倍率である、開口数が大きいなどの理由により、物体光学系23の被写界深度が、複素共役信号の間隔で決まる有効高さZeよりも小さい場合があり得る。この場合、有効領域Reの端部が物体光学系23の合焦範囲から部分的に外れることを許容した上で上記の有効領域Reを維持してもよく、また有効領域Reを物体光学系23の被写界深度で決まる深さの範囲に限定するようにしてもよい。
撮像が終了すると(ステップS109においてYES)、信号処理部33は、各焦点深さで取得された部分画像Ipを合成することで、XZ平面に平行な1つの断面における被撮像物の断層画像Itを作成する(ステップS111)。必要に応じて、Y方向に位置を変えながら上記の撮像動作を繰り返せば、Y方向に位置の異なる複数の断層画像Itが得られる。3D復元部34は、これらの断層画像Itから被撮像物の立体像を作成することができる。
次に、設定された焦点深さDが容器底部厚さの半分(T/2)よりも小さい場合に(ステップS104においてNO)、R=Tとして参照ミラー24の位置を設定する理由について説明する。これまで述べてきたように、式(4)に基づいて参照ミラー24の位置を設定する理由は、物体光学系23の合焦範囲を有効に利用するために、容器底面に起因する2つの複素共役信号の中央に焦点位置が来るようにするためである。このとき、複素共役信号のピーク幅を無視すれば、有効領域Reの範囲は焦点深さを中心として(+Z)方向および(−Z)方向にそれぞれ(T/2)ずつとなる(図7(d)参照)。
図10は焦点深さが容器底部厚さの半分未満であるときの反射光ピーク位置を示す図である。図10(a)に示すように、焦点面Sfが容器11の上部底面Saに近く、焦点深さDが容器底部厚さTの半分より小さい場合を考える。このとき、式(4)に従って2つの複素共役信号P2a,P2bの中央に実信号P1aが来るように参照基準面を設定すると、図10(b)に示すように、2つの複素共役信号P2a,P2bが、実信号P1aの位置からそれぞれ(T/2)ずつ離れた位置に出現する。
しかしながら、実信号P1aからDだけ離れた位置に、容器11の上部底面Saに対応する実信号P1aが現れている。D<(T/2)であるから、このときの有効領域Reは、容器11の上部底面Saに対応する実信号P1aと、下部底面Sbに対応する複素共役信号P2bとに挟まれた領域に限定されてしまう。つまり、2つの複素共役信号P2a,P2bの中央に実信号P1aが来るように参照ミラー24の位置を設定する意義がなく、むしろ有効領域Reを狭くしてしまう。
そこで、図10(c)に示すように、容器11の上部底面Saに対応する実信号P1aとこれに共役な複素共役信号P2aとが重なるように、つまりR=Tとなるように参照ミラー24の位置が設定される。この場合の有効領域Reは、容器11の上部底面Saに対応する実信号P1aまたはこれに共役で同じ位置に現れる複素共役信号P2aと、もう1つの複素共役信号P2bとで挟まれた領域となり、図10(b)に示す場合よりも有効領域Reを広くすることができる。
これを可能とするために、図9の撮像動作においては、設定された焦点深さDについてD<(T/2)の関係が成立するとき、式(4)によらずR=Tとなるように参照ミラー24の位置が設定される。このことは、参照基準面Srを容器11の上部底面Saに一致させることと等価である。
これにより、物体光学系23の合焦範囲を有効に利用して撮像を行うことができる。複数の部分画像Ipを合成して断層画像Itを作成する構成においては、各撮像における物体光学系23の合焦範囲を最大限に利用することで必要な部分画像Ipの数を削減し、画像の作成に要する時間を短縮することが可能になる。D<(T/2)の関係がそれぞれ成立する複数の焦点深さで撮像が行われる際にも、参照ミラー24は焦点深さの設定によらずR=Tとなる位置に固定される。参照ミラー24が移動する工程を省くことができるので、処理時間も短縮することが可能である。
以上のように、この実施形態の画像処理装置1は、フーリエドメイン(FD)OCT撮像原理を利用した撮像装置であるが、物体光学系23の焦点深さを変更するための機能を有している。そして、容器11底部を介して撮像が行われる際、焦点深さの設定に応じて、参照基準面を規定する参照光路の光路長が参照ミラー24の位置調整により設定される。参照基準面の位置を導出する計算式は、焦点深さDおよび容器11底部の厚さTの関数として与えられる。
このような構成によれば、深さ方向の反射光強度分布において、容器11の上部底面Saと下部底面Sbとに対応する複素共役信号が焦点深さに対応する位置を挟んで出現するようにすることができる。2つの複素共役信号の間隔は容器11底部の厚さTにより決まる一定値であるので、結果として、焦点位置およびその近傍の一定幅の範囲において複素共役ノイズの影響を排除した撮像を行うことが可能となる。
必要に応じて、深さ方向に焦点位置を異ならせた複数回の撮像で得られた結果を合成することで、深さ方向において2つの複素共役信号の間隔よりも広い範囲を含む断層画像を作成することが可能である。このとき、焦点位置の変更に伴って参照ミラー24の位置を変更することにより、焦点位置の近傍では常に複素共役ノイズの影響を排除することができる。しかも合焦状態で撮像された比較的狭い範囲の部分画像を合成してより広範囲の断層画像を作成することで、断層画像全体の画像品質を良好なものとすることができる。
以上説明したように、この実施形態の画像処理装置1は、スフェロイドSp等を「被撮像物」とする本発明の「撮像装置」に相当するものである。そして、撮像ユニット20,20a,20bがそれぞれ本発明の「検出手段」として機能し、制御ユニット30のうちCPU31が本発明の「制御手段」として、また信号処理部33が本発明の「信号処理手段」としてそれぞれ機能している。
また、上記実施形態においては、容器11の底部111が本発明の「壁部」に相当し、その上部底面Saおよび下部底面Sbがそれぞれ本発明の「第1表面」および「第2表面」に相当している。
なお、本発明は上記した実施形態に限定されるものではなく、その趣旨を逸脱しない限りにおいて上述したもの以外に種々の変更を行うことが可能である。例えば、上記実施形態における参照基準面Srの位置の導出処理では、物体光学系23の焦点深さが容器11の上部底面Saから起算した値Dによって、また参照基準面Srの位置が容器11の下部底面Sbから起算した値Rによって表されている。しかしながら、どの位置を基準として各位置を表すかは任意であり、それに応じて各式を適宜改変することにより、上記技術思想と実質的に等価な処理を実現することが可能である。
また例えば、上記実施形態の被撮像物は浅皿型のディッシュと呼ばれる容器11に担持されたスフェロイドSpであるが、被撮像物およびそれを担持する容器の種類はこれに限定されない。例えば、試料を担持可能な複数のウェルが1つのプレートに設けられたウェルプレートで培養された細胞等が被撮像物とされてもよい。
また例えば、上記実施形態の画像処理装置1は、容器11の底部111を介して容器内の被撮像物を撮像する。しかしながら、撮像の方向はこれに限定されず、例えば被撮像物を収容する容器の側壁面を介して撮像を行う場合にも、本発明を好適に適用することが可能である。
また例えば、上記実施形態では、焦点深さDが容器底部厚さTの半分よりも大きいか小さいかによって値Rの求め方が変更される。しかしながら、式(4)による値Rの求め方は、焦点深さDが容器底部厚さTの半分よりも小さい場合にも有効である。例えば有効領域Reの範囲が狭くなることが許容される状況であれば、焦点深さDの可変範囲の全体において一律に式(4)が適用されてもよい。
また例えば、上記実施形態では容器底面に起因する2つの複素共役信号のピークに挟まれる領域の中央に焦点位置が来るように式(4)が採用されているが、焦点位置が2つのピークのちょうど中間である必要は必ずしもない。例えば、2つのピークの間隔よりも物体光学系の被写界深度の方が大きければ、焦点位置が両ピークの中央からずれていても両ピーク間の領域全体を合焦範囲に含めることが可能であり、得られる画像の品質も上記実施形態と変わりないものとなる。式(4)に適宜のオフセットを加えた式、あるいは、図7(d)の値「T/2」を適宜に改変して式(4)と同様にして求めた式を用いて、このような焦点位置と複素共役信号のピークとの位置関係を実現することが可能である。
また例えば、上記実施形態では、焦点深さを多段階に設定して撮像された部分画像Ipを合成することで断層画像Itが得られる。しかしながら、例えば被撮像物の高さに対して容器底部厚さおよび物体光学系の被写界深度が十分に大きければ、1回の撮像で被撮像物の全体像を得ることが可能である。この場合にも、焦点位置の設定に応じて、上記した条件を満たすように参照光の光路長が設定されることにより、撮像結果に容器底面に起因する複素共役ノイズの影響が現れるのを効果的に防止することができる。
また、本発明において物体光学系の分解能は特に限定されないが、高倍率または高分解能での撮像を必要とする場合において本発明は特に有効である。その理由は以下の通りである。高倍率、高分解能を求めると物体光学系の被写界深度が浅くなるため、本来的には深さ方向に広く撮像可能なFD−OCT装置であっても、1回の撮像で得られる画像の深さ方向の範囲は限定される。したがって、本発明の技術思想とは直接関係のない理由で複数回の撮像が必要となる場合がある。このような場合に本発明を適用することで、1回の撮像で得られる画像は複素共役ノイズの影響がなく、また全体を合焦状態とすることができるので、良好な画像品質で撮像を行うことが可能となる。
また、上記実施形態の制御ユニット30としては、パーソナルコンピュータやワークステーション等の一般的な構成の汎用処理装置を用いることも可能である。すなわち、撮像ユニット20、駆動制御部40およびこれらの動作させるための最小限の制御機能を有する撮像装置と、上記処理内容を記述した制御プログラムを実行することで制御ユニット30として機能するパーソナルコンピュータ等との組み合わせにより、画像処理装置1が構成されてもよい。
以上、具体的な実施形態を例示して説明してきたように、この発明においては、一の撮像において得られる反射光強度分布から求められる断層画像のうち第1表面の共役像と第2表面の共役像との間の領域が、有効な画像領域とされてもよい。このような構成によれば、第1表面および第2表面の共役像に起因するノイズを含まない画像領域を有効なものとして抽出することができる。
また、壁部の厚さをT、第1表面から物体光学系の焦点までの距離をD、第2表面から、照明光の光路に垂直な平面であって当該平面までの照明光の光路長が参照光の光路長と等しい参照基準面までの距離をRとするとき、
(D+T)/2<R<(D+2T)/2
の関係が成立するように、参照光の光路長が設定されてもよい。値Rがこのような条件を満たすように参照基準面が設定されるとき、第1表面および第2表面の共役像の出現位置を、確実に物体光学系の焦点を挟む位置とすることができ、本発明の目的が達成される。
この場合の値Rは、予め定められた値DおよびTを変数とする関数F(D,T)に、値Dおよび値Tを代入することで求められるものであってもよい。このような関数を予め準備しておくことで、値D,Tが決まれば直ちに参照光の光路長を決定することが可能になる。
より具体的には、関数F(D,T)は次式:
F(D,T)=D/2+3T/4
で表されるものであってもよい。前述した通り、このような関数式に基づき参照光の光路長を設定すれば、第1表面および第2表面の共役像のちょうど中央に当たる位置が物体光学系の焦点位置となる。これにより焦点位置前後における物体光学系の合焦範囲の検出結果を有効に利用して、画像品質の良好な断層画像を作成することが可能になる。
また、同一の被撮像物に対し焦点深さを互いに異ならせて複数回干渉光の検出が行われ、焦点深さを設定する度に参照光の光路長が設定され、得られた複数の検出結果から各々の焦点近傍における反射光強度分布が求められる構成であってもよい。本発明の原理によれば、一の撮像で複素共役ノイズを排除した反射光強度分布が得られる深さの範囲は容器壁部の厚さにより制限される。焦点深さを互いに異ならせた複数回の検出結果を重ね合わせれば、容器壁部の厚さを超える範囲を含む、画像品質の良好な断層画像を得ることが可能となる。
また、この発明にかかる撮像装置は、検出手段が参照光の光路に配置されて参照光の光路長を規定する参照ミラーを備え、制御手段が、参照ミラーの位置を変化させて参照光の光路長を調整するミラー駆動機構と、物体光学系を駆動して焦点深さを調整する焦点調整機構とを備えるものであってもよい。このような構成は、物体光学系の焦点深さに応じて参照ミラーの位置を変化させることで参照光の光路長を変更設定することができるものであり、上記発明を実施するのに好適なものである。
また、この発明にかかる撮像方法は、参照光の光路長を設定する前に、壁部の厚さに関する情報を取得する工程を備えるものであってもよい。このような構成によれば、焦点深さの設定値と壁部の厚さに関する情報とに基づき、参照光の光路長を上記条件に合致するように調整することが容易となる。
この発明は、FD−OCT撮像技術全般に適用することができる。特に、ディッシュ等の容器中で培養された細胞や細胞集塊を撮像する医学・生化学・創薬の分野において好適に適用することができる。
1 画像処理装置(撮像装置)
11 容器
20,20a,20b 撮像ユニット(検出手段)
22 ビームスプリッタ
23 物体光学系
24 参照ミラー
30 制御ユニット
31 CPU(制御手段)
33 信号処理部(信号処理手段)
40 駆動制御部
41 焦点調整機構
42 ミラー駆動機構
111 (容器11の)底部(壁部)
220 光ファイバカプラ
Sa (容器11の)上部底面(第1表面)
Sb (容器11の)下部底面(第2表面)
Sf 焦点面
Sp スフェロイド(被撮像物)
Sr 基準参照面

Claims (11)

  1. 光透過性の壁部を有する容器内に収容された被撮像物を断層撮像する撮像装置において、
    光源から出射される広帯域の低コヒーレンス光が分岐された一の分岐光を照明光として前記壁部を介して前記被撮像物に入射させ、前記壁部を介して出射される前記被撮像物の反射光を物体光学系により集光した信号光と、他の一の分岐光から生成された参照光とが干渉して生じる干渉光を検出し、検出された前記干渉光に応じた干渉信号を出力する検出手段と、
    前記干渉信号に基づき、前記干渉光のスペクトルをフーリエ変換して前記照明光の入射方向における前記被撮像物の反射光強度分布を求め、該反射光強度分布から断層画像を作成する信号処理手段と、
    前記被撮像物に対する前記物体光学系の焦点深さを変更設定する制御手段と
    を備え、
    前記制御手段は、深さ方向の位置ごとの反射光強度を表す前記反射光強度分布において、前記焦点深さに対応する位置が、前記壁部の表面のうち前記被撮像物側の第1表面に共役な位置と、前記被撮像物とは反対側の第2表面に共役な位置との間となるように、前記焦点深さの設定に連動させて前記参照光の光路長を変更設定する撮像装置。
  2. 前記信号処理手段は、前記反射光強度分布から求められる前記断層画像のうち前記第1表面の共役像と前記第2表面の共役像との間の領域を有効な画像領域とする請求項1に記載の撮像装置。
  3. 前記壁部の厚さをT、
    前記第1表面から前記物体光学系の焦点までの距離をD、
    前記第2表面から、前記照明光の光路に垂直な平面であって当該平面までの前記照明光の光路長が前記参照光の光路長と等しい参照基準面までの距離をRとするとき、
    (D+T)/2<R<(D+2T)/2
    の関係が成立するように前記参照光の光路長が設定される請求項1または2に記載の撮像装置。
  4. 前記Rは、予め定められた前記Dおよび前記Tを変数とする関数F(D,T)に、値Dおよび値Tを代入することで求められる請求項3に記載の撮像装置。
  5. 前記関数F(D,T)は次式:
    F(D,T)=D/2+3T/4
    で表される請求項4に記載の撮像装置。
  6. 前記信号処理手段は、同一の前記被撮像物に対し前記焦点深さを互いに異ならせて得られた複数の検出結果から各々の焦点近傍における前記反射光強度分布を求める請求項1ないし5のいずれかに記載の撮像装置。
  7. 前記検出手段は、前記参照光の光路に配置されて前記参照光の光路長を規定する参照ミラーを備え、
    前記制御手段は、前記参照ミラーの位置を変化させて前記参照光の光路長を調整するミラー駆動機構と、前記物体光学系を駆動して前記焦点深さを調整する焦点調整機構とを備える請求項1ないし6のいずれかに記載の撮像装置。
  8. 光透過性の壁部を有する容器内に収容された被撮像物を断層撮像する撮像方法において、
    光源から出射される広帯域の低コヒーレンス光が分岐された一の分岐光を照明光として前記壁部を介して前記被撮像物に入射させ、前記壁部を介して出射される前記被撮像物の反射光を物体光学系により集光した信号光と、他の一の分岐光から生成された参照光とが干渉して生じる干渉光を検出し、検出された前記干渉光に応じた干渉信号を出力する工程と、
    前記干渉信号に基づき、前記干渉光のスペクトルをフーリエ変換して前記照明光の入射方向における前記被撮像物の反射光強度分布を求め、該反射光強度分布から断層画像を作成する工程と
    を備え、
    前記被撮像物に対する前記物体光学系の焦点深さを変更設定可能であり、
    深さ方向の位置ごとの反射光強度を表す前記反射光強度分布において、前記焦点深さに対応する位置が、前記壁部の表面のうち前記被撮像物側の第1表面に共役な位置と、前記被撮像物とは反対側の第2表面に共役な位置との間となるように、前記焦点深さの設定に連動させて前記参照光の光路長を変更設定する撮像方法。
  9. 前記反射光強度分布から求められる前記断層画像のうち前記第1表面の共役像と前記第2表面の共役像との間の領域を有効な画像領域とする請求項8に記載の撮像方法。
  10. 同一の前記被撮像物に対し前記焦点深さを多段階に互いに異ならせて複数回前記干渉光の検出を行い、前記焦点深さを設定する度に前記参照光の光路長を設定する請求項9に記載の撮像方法。
  11. 前記参照光の光路長を設定する前に、前記壁部の厚さに関する情報を取得する工程を備える請求項8ないし10のいずれかに記載の撮像方法。
JP2016185552A 2016-09-23 2016-09-23 撮像装置および撮像方法 Active JP6728007B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2016185552A JP6728007B2 (ja) 2016-09-23 2016-09-23 撮像装置および撮像方法
EP17191736.2A EP3299761B1 (en) 2016-09-23 2017-09-19 Imaging apparatus and imaging method
US15/711,684 US10317188B2 (en) 2016-09-23 2017-09-21 Imaging apparatus and imaging method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016185552A JP6728007B2 (ja) 2016-09-23 2016-09-23 撮像装置および撮像方法

Publications (2)

Publication Number Publication Date
JP2018048954A true JP2018048954A (ja) 2018-03-29
JP6728007B2 JP6728007B2 (ja) 2020-07-22

Family

ID=59930184

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016185552A Active JP6728007B2 (ja) 2016-09-23 2016-09-23 撮像装置および撮像方法

Country Status (3)

Country Link
US (1) US10317188B2 (ja)
EP (1) EP3299761B1 (ja)
JP (1) JP6728007B2 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10401793B2 (en) * 2010-06-17 2019-09-03 Purdue Research Foundation Digital holographic method of measuring cellular activity and measuring apparatus with improved stability
JP6702851B2 (ja) * 2016-12-26 2020-06-03 株式会社Screenホールディングス 撮像装置および撮像方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008541096A (ja) * 2005-05-13 2008-11-20 ザ ジェネラル ホスピタル コーポレイション 化学的試料および生体試料の高感度検出用スペクトル領域光コヒーレンス反射計測を実行可能な装置、システム、および方法
JP2010169503A (ja) * 2009-01-22 2010-08-05 Canon Inc 光断層画像撮像装置
JP2011024842A (ja) * 2009-07-28 2011-02-10 Canon Inc 光断層画像撮像装置
JP2012513260A (ja) * 2008-12-23 2012-06-14 カール ツァイス メディテック アクチエンゲゼルシャフト 掃引信号源光学コヒーレンスドメイン反射率測定用の装置
JP2012202774A (ja) * 2011-03-24 2012-10-22 Nikon Corp 観察装置及び観察方法
JP2014076403A (ja) * 2014-02-06 2014-05-01 Canon Inc 撮像装置、制御装置及び制御方法
US20140128731A1 (en) * 2012-11-02 2014-05-08 Optimedica Corporation Optical surface identification for laser surgery
JP2015072152A (ja) * 2013-10-02 2015-04-16 株式会社日立エルジーデータストレージ 光計測装置
JP2015175678A (ja) * 2014-03-14 2015-10-05 株式会社日立エルジーデータストレージ 光断層観察装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7554669B2 (en) * 2005-11-15 2009-06-30 Bioptigen, Inc. Spectral domain phase microscopy (SDPM) dual mode imaging systems and related methods and computer program products
JP4852651B2 (ja) 2010-02-24 2012-01-11 国立大学法人 筑波大学 多重化スペクトル干渉光コヒーレンストモグラフィー
US8564788B2 (en) * 2011-07-26 2013-10-22 Kabushiki Kaisha Topcon Optical imaging method and optical imaging apparatus
US20150026050A1 (en) * 2013-07-17 2015-01-22 Jvl Ventures, Llc Systems, methods, and computer program products for reporting contactless transaction data
JP6231958B2 (ja) * 2014-08-20 2017-11-15 株式会社日立エルジーデータストレージ 光画像計測装置

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008541096A (ja) * 2005-05-13 2008-11-20 ザ ジェネラル ホスピタル コーポレイション 化学的試料および生体試料の高感度検出用スペクトル領域光コヒーレンス反射計測を実行可能な装置、システム、および方法
JP2012513260A (ja) * 2008-12-23 2012-06-14 カール ツァイス メディテック アクチエンゲゼルシャフト 掃引信号源光学コヒーレンスドメイン反射率測定用の装置
JP2010169503A (ja) * 2009-01-22 2010-08-05 Canon Inc 光断層画像撮像装置
JP2011024842A (ja) * 2009-07-28 2011-02-10 Canon Inc 光断層画像撮像装置
JP2012202774A (ja) * 2011-03-24 2012-10-22 Nikon Corp 観察装置及び観察方法
US20140128731A1 (en) * 2012-11-02 2014-05-08 Optimedica Corporation Optical surface identification for laser surgery
JP2015072152A (ja) * 2013-10-02 2015-04-16 株式会社日立エルジーデータストレージ 光計測装置
JP2014076403A (ja) * 2014-02-06 2014-05-01 Canon Inc 撮像装置、制御装置及び制御方法
JP2015175678A (ja) * 2014-03-14 2015-10-05 株式会社日立エルジーデータストレージ 光断層観察装置

Also Published As

Publication number Publication date
EP3299761A1 (en) 2018-03-28
EP3299761B1 (en) 2019-06-12
JP6728007B2 (ja) 2020-07-22
US20180087892A1 (en) 2018-03-29
US10317188B2 (en) 2019-06-11

Similar Documents

Publication Publication Date Title
CA2934556C (en) Optical tomography apparatus and method
EP2798321B1 (en) Spectroscopic instrument and process for spectral analysis
US10113856B2 (en) Line-field imaging systems and methods incorporating planar waveguides
JP5975522B2 (ja) 動的焦点移動型光干渉断層顕微鏡
JP6762810B2 (ja) 画像処理装置および画像処理方法
CN107567305B (zh) 摄像设备
JP7113637B2 (ja) 光干渉断層撮像装置および光干渉断層撮像方法
JP5527478B2 (ja) 光干渉断層観察装置、画像間の相対位置決定方法および画像間の相対位置決定プログラム
JP6728007B2 (ja) 撮像装置および撮像方法
JP6702851B2 (ja) 撮像装置および撮像方法
JP5990224B2 (ja) 撮像装置及び撮像方法
WO2018150744A1 (ja) 撮像方法および撮像装置
EP3162282A1 (en) Tomography device
JP2018132466A (ja) 撮像方法および撮像装置
EP2565625A1 (en) Optical measurement system and method for operating an optical measurement system
JP2021173723A (ja) 撮像装置および撮像方法
EP3999897A1 (en) High sensitivity optical devices using pupil splittig means
JP2021004749A (ja) 補正方法、補正装置および撮像装置
JP2021524330A (ja) Octシステム及びoct方法

Legal Events

Date Code Title Description
RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20170725

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190624

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200311

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: 20200609

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200701

R150 Certificate of patent or registration of utility model

Ref document number: 6728007

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250