(57)【要約】【課題】光源の波数変化特性を考慮したリサンプリング処理を実現する。【解決手段】光源からの光の波長掃引を行って対象の物体に照射し、該物体からの反射光を、干渉計により干渉波形信号として取得する光干渉断層撮影法を適用したデータ処理装置におけるリサンプリング方法である。前記干渉波形信号と、前記光源の動作開始時からの経過時間とを取得する。前記光源からの光の波数−時間関数を確定するために、取得した前記経過時間に基づいて、前記光源の動作開始時からの経過時間に応じた波数−時間関数が格納されたデータテーブルを参照する。前記データテーブルを参照して得られた波数−時間関数の波数が一定間隔となる、前記干渉波形信号のリサンプリング点を求める。前記干渉波形信号を、前記リサンプリング点でリサンプリングする。【選択図】図９
2014228397.tif 000002
【特許請求の範囲】
【請求項１】
光源からの光の波長掃引を行って対象の物体に照射し、該物体からの反射光を、干渉計により干渉波形信号として取得する光干渉断層撮影法を適用したデータ処理装置におけるリサンプリング方法であって、
前記データ処理装置の取得手段が、前記干渉波形信号と、前記光源の動作開始時からの経過時間とを取得する第１ステップと、
前記データ処理装置の参照手段が、前記光源からの光の波数−時間関数を確定する第２ステップであって、前記第１ステップで取得した前記経過時間に基づいて、前記光源の動作開始時からの経過時間に応じた波数−時間関数が格納されたデータテーブルを参照する第２ステップと、
前記データ処理装置の算出手段が、前記データテーブルを参照して得られた波数−時間関数の波数が一定間隔となる、前記干渉波形信号のリサンプリング点を求める第３ステップと、
前記データ処理装置のリサンプリング手段が、前記干渉波形信号を、前記リサンプリング点でリサンプリングする第４ステップと
を備えたことを特徴とするリサンプリング方法。
【請求項２】
前記光源は、半導体光増幅器、ＫＴＮ（ＫＴａ1-ｘＮｂｘＯ3（０＜ｘ＜１））結晶からなるＫＴＮ光偏向器、回折格子およびミラーからなるLittman型外部共振器で構成されていることを特徴とする請求項１に記載のリサンプリング方法。
【請求項３】
光源からの光の波長掃引を行って対象の物体に照射し、該物体からの反射光を、干渉計により干渉波形信号として取得する光干渉断層撮影法を適用したデータ処理装置であって、
前記干渉波形信号と、前記光源の動作開始時からの経過時間とを取得する取得手段と、
前記光源からの光の波数−時間関数を確定する参照手段であって、前記取得手段で取得した前記経過時間に基づいて、前記光源の動作開始時からの経過時間に応じた波数−時間関数が格納されたデータテーブルを参照する参照手段と、
前記データテーブルを参照して得られた波数−時間関数の波数が一定間隔となる、前記干渉波形信号のリサンプリング点を求める算出手段と、
前記干渉波形信号を、前記リサンプリング点でリサンプリングするリサンプリング手段と
を備えたことを特徴とするデータ処理装置。
【請求項４】
前記光源は、半導体光増幅器、ＫＴＮ（ＫＴａ1-ｘＮｂｘＯ3（０＜ｘ＜１））結晶からなるＫＴＮ光偏向器、回折格子およびミラーからなるLittman型外部共振器で構成されていることを特徴とする請求項３に記載のデータ処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、データ処理装置およびリサンプリング方法に関し、より詳細には、主として医療で使用されている、３次元断層画像を取得するための光干渉断層撮影法のデータ処理装置において、波数掃引（波数が時間に対して線形的に変化する掃引）ではない波長可変光源を用いた場合に得た干渉波形信号を、波数掃引した波長可変光源を用いた場合と同等の干渉波形信号に変換（リサンプリング）する方法に関する。
【背景技術】
【０００２】
光干渉断層撮影法（OCT：Optical Coherence Tomography）は、赤外線の干渉を用いて断層画像を取得する方法であり、眼底検査などの医療分野で用いられている。赤外線光源として広帯域光源を用いて、多くの波長を含んだ光を対象の物体に照射し、反射した光をＣＣＤカメラ等の受光素子で取得する。反射した光は、波長によって位相のずれ方が変わるため、各々の波長における干渉強度が異なるので、分光器を通して波長分解することにより、各波長ごとの干渉信号を取得することができる。この干渉信号をフーリエ変換することにより、物体の反射面の情報を得ることができる。この方法は、SD-OCT（Spectral Domain-OCT）と呼ばれている。また、全ての波長の光を同時に入射するのではなく、光源からの出射光の波長掃引を行って、波長ごとに干渉波形信号を時間的に取得することも行われている。この方法は、SS-OCT（Swept Source-OCT）と呼ばれている。SS-OCTは、分光器が不要で、ＣＣＤカメラを用いる必要がないので、装置を小型化でき、SD-OCTよりも高速にデータを取得することができる点で優れている。
【０００３】
SS-OCT画像解析システムは、光源と受光素子とを含む光学系と、光源の波長掃引を制御、受光素子で取得された干渉波形信号の格納、および格納された干渉波形信号の解析を行うためのデータ処理装置とから構成されている。また、SS-OCT画像解析システムの光源として、ＫＴＮ（ＫＴａ1-ｘＮｂｘＯ3（０＜ｘ＜１））結晶からなる光偏向器を利用したＫＴＮ光源が知られている。
【０００４】
図１に、従来のSS-OCT光源に用いられるＫＴＮ光源の構成を示す。ＫＴＮ光源は、ＳＯＡ１０４と、ＳＯＡ１０４からの出射光を平行光にするコリメータレンズ１０６と、ＫＴＮ結晶からなるＫＴＮ光偏向器１０７と、反射型の回折格子１１１と、ミラー１１２とを備え、ＳＯＡ１０４の反射面１０５とミラー１１２との間で外部共振器、いわゆるLittman型外部共振器を構成している（例えば、非特許文献１参照）。ＳＯＡ１０４の反射面１０５からの出力光、すなわちＫＴＮ光源からの出力光は、集光レンズ１０３を介して光ファイバ１０１に結合させる。アイソレータ１０２は、光ファイバ１０１に接続された機器からの戻り光が、ＫＴＮ光源に入射するのを防いでいる。
【０００５】
ＫＴＮ光偏向器１０７は、電極１０８，１０９に電圧を印加することにより、ＫＴＮ光偏向器１０７内部の電界によって形成された屈折率分布によって、内部を透過するＳＯＡ１０４からの出射光の向きを変える。電源１１０は、電極１０８，１０９に正弦波の電圧を出力する。回折格子１１１は、コリメータレンズ１０６の光軸に対して角度αの傾きを持って配置されている。ミラー１１２は、コリメータレンズ１０６の光軸に垂直な軸に対して角度φの傾きを持って配置されている。なお、図１においては、時計回り方向が角度のプラス方向である。
【０００６】
ＳＯＡ１０４からの出射光は、コリメータレンズ１０６によって平行光となり、ＫＴＮ光偏向器１０７に入射する。ＫＴＮ光偏向器１０７に入射した光は、進行方向が曲り（偏向し）、コリメータレンズ１０６の光軸に対して角度Ψで出射される。ＫＴＮ光偏向器１０７から出射した光は、回折格子１１１に入射角α’で入射し、光の波長λと入射角α’で決まる出射角βで出射される。光の波長λ、入射角α’、出射角βの関係は、回折格子方程式と呼ばれる式で決められており、以下の関係がある。
【０００７】
【数１】
2014228397.tif 000003
【０００８】
ただし、Λは回折格子の溝間隔、ｍは回折の次数である。
【０００９】
回折格子１１１から出射した光は、ミラー１１２で垂直に反射された光だけが、回折格子１１１、ＫＴＮ光偏向器１０７、コリメータレンズ１０６を介してＳＯＡ１０４に戻る。従って、ＫＴＮ光偏向器１０７から出射する光の偏向角（出射角）Ψと、回折格子１１１の配置の角度αと、ミラー１１２の配置の角度φとによって、レーザ発振する波長λが以下のように決まる。
【００１０】
【数２】
2014228397.tif 000004
【００１１】
つまり、ＫＴＮ光偏向器１０７から出射する光の偏向角（出射角）Ψによって、レーザ発振する波長λを制御することができる。
【００１２】
図２に、従来のSS-OCTにおける光学系の基本構成を示す。ＫＴＮ光源１２１は、時間ｔによって、出射する光の波数ｋ（ｔ）が変化する光源である。光カプラ１２２は、ＫＴＮ光源１２１からの出射光を２方向に分ける。一方の光ｌｒは、サーキュレータ１２３とレンズ１２４とを介して、参照ミラー１２５に入射される。光ｌｒの光軸に対して垂直に置かれている参照ミラー１２５によって反射された光は、サーキュレータ１２３と光カプラ１２７とを透過して、受光素子１２６に入力される。光カプラ１２２と参照ミラー１２５との間は距離ｌだけ離れている。
【００１３】
光カプラ１２２で分岐された他方の光ｌｏは、サーキュレータ１２８、駆動ミラー１２９およびレンズ１３０を介して、測定対象物１３１の表面（以下、基準面）に入射される。基準面は、光ｌｏの光軸と垂直な面であり、光カプラ１２２と基準面との間は距離ｌだけ離れている。測定対象物１３１の表面、測定対象物１３１内の反射面で反射された光は、サーキュレータ１２８と光カプラ１２７とを透過して、受光素子１２６に入力される。受光素子１２６は、参照ミラー１２５からの反射光（ｌｒ）と測定対象物１３１の表面または測定対象物１３１内の反射面からの反射光（ｌｏ）の両方を同時に受光するので、両者の干渉信号を観測することができる。
【００１４】
この光学系には、ＫＴＮ光源１２１の波長掃引を制御、受光素子１２６で取得された干渉波形信号の格納、および格納された干渉波形信号の解析を行うためのデータ処理装置１３２と、駆動ミラー１２９を制御して測定対象物１３１の表面を走査するためのドライバ回路１３３とが接続されている。
【００１５】
参照ミラー１２５で反射された光の受光素子１２６上での電界を、Er(t)e-jk(t)lとし、測定対象物１３１の表面または測定対象物１３１内の反射面で反射された光の受光素子１２６上での電界を、E0(t)e-jk(t)(l+z)とする。ここで、Ｅｒ（ｔ）はミ参照ミラー１２５で反射された光の振幅、Ｅ０（ｔ）は測定対象物１３１の表面または測定対象物１３１内の反射面で反射された光の振幅、ｔは時間、ｅはネーピア数、ｊは虚数単位、ｋ（ｔ）は光の波数、ｌは光カプラ１２２から参照ミラー１２５までの距離、および、光カプラ１２２から基準面までの距離、ｚは基準面から測定対象物１３１の表面または測定対象物１３１内の反射面までの距離を表す。受光素子上の光の強度は、
【００１６】
【数３】
2014228397.tif 000005
【００１７】
で表される。ただし、γ（ｚ）はコヒーレンス関数であり、光路長差ｚの場合の干渉の強度を表す。受光素子で観測できる干渉波形信号ｓｆは、
【００１８】
【数４】
2014228397.tif 000006
【００１９】
となる。ただし、ηは受光素子の量子効率である。
【００２０】
もし、波数ｋ（ｔ）＝ｂｔ＋ｃ（ただし、ｂとｃは定数）であれば、
【００２１】
【数５】
2014228397.tif 000007
【００２２】
となるので、干渉波形信号ｓｆは時間周波数ｆ＝ｂｚ／πで振動する。したがって、干渉波形信号ｓｆの時間変化する信号ｓｆ（ｔ）の周波数ｆから、測定対象物１３１の表面または測定対象物１３１内の反射面と基準面との距離ｚが分かる。たとえば、干渉波形信号ｓｆ（ｔ）の時間軸方向のフーリエ変換結果Ｓｆ（ｆ）のピークとなる周波数をｆ０とすると、基準面から測定対象物１３１の表面または測定対象物１３１内の反射面までの距離はｚ０＝πｆ０／ｂと算出できる。このことから、Ｓｆ（ｆ）の周波数ｆ軸を距離ｚ＝πｆ／ｂと変換することにより、深さ方向の情報を得ることができる。基準面に対して平行な方向に少しずつ移動しながらＳｆ（ｆ）を取得すると断層画像が得られる。このとき、フーリエ変換結果Ｓｆ（ｆ）のピークが先鋭化するほど、鮮明な断層画像が得られる。
【００２３】
しかし、通常使われているSS-OCT光源は、上記のような波数ｋ（ｔ）＝ｂｔ＋ｃ（ただし、ｂとｃは定数）となってはいないことから、干渉波形信号ｓｆ（ｆ）の時間軸方向のフーリエ変換結果Ｓｆ（ｆ）は明確なピークをもたず、従って、鮮明な断層画像を得ることができない。この問題を解決する方法の一つとして、干渉波形信号ｓｆ（ｆ）を時間軸上で伸縮させて、等価的に波数がｋ（ｔ）＝ｂｔ＋ｃに従う光源を使った時と同じ干渉波形信号となるように、干渉波形信号ｓｆ（ｆ）を変換（リサンプリング）する処理が行われている。
【００２４】
図３に、干渉波形信号をリサンプリング処理を行った場合の効果を示す。図３（ａ）は、波数ｋ（ｔ）＝ｂｔ＋ｃとなってはいない光源の周波数変化特性を示す。横軸は時間、縦軸は相対波数である。図３（ｂ）は、このような光源により取得した干渉波形信号ｓｆ（ｔｕ）であり、SS-OCTで得られた信号の点広がり関数（ＰＳＦ：point spread function）を図３（ｃ）に示す。
【００２５】
一方、図３（ａ）に示した既知の周波数変化特性により、測定された干渉波形信号の時間間隔をリサンプリング処理した場合の干渉波形信号ｓｆ’’（ｔｕ）を図３（ｄ）に示す。図３（ｅ）は、リサンプリング処理した場合のＰＳＦであり、図３（ｃ）と比較すると、ＰＳＦが小さく、精度の高い画像が得られていることがわかる。
【００２６】
図４に、従来のデータ処理装置におけるリサンプリング処理のための機能ブロック図を示す（例えば、非特許文献１参照）。図５を参照して、従来の干渉波形信号の変換（リサンプリング）方法を説明する。
【００２７】
干渉波形信号取得手段１４１は、基準面から距離ｚ０（既知）の位置に反射面を持つ測定対象物１３１に対して、図２に示したような干渉計と受光素子からなる光学系を使って、第１の干渉波形信号（時間−強度の信号）
【００２８】
【数６】
2014228397.tif 000008
【００２９】
を取得して出力する（Ｓ１６１）。ただし、ｔｕは時間であり離散値であり、ｋ（ｔｕ）は波数である。また、ｕはその離散値を識別する整数であり、第１の干渉波形信号がＮ個ある場合は、ｕは１〜Ｎの値を取る。
【００３０】
リサンプリング点算出手段１４８は、波数変化特性データ参照手段１５０により予め取得されている波数変化特性（波数−時間特性）を、ｎ次多項式近似により導出したｋ（ｔ）の逆関数ｔ（ｋ）の近似関数ｔ’（ｋ）を使って、波数が一定間隔となる波数値ｋν（νは波数値同士を区別する整数）を選び（ｋν−ｋν＋１＝Δｋ; Δｋは一定値）、
【００３１】
【数７】
2014228397.tif 000009
【００３２】
を干渉波形ｓｆ（ｔｕ）のリサンプリング点として出力する（Ｓ１６８）。
【００３３】
リサンプリング手段１４９は、干渉波形信号取得手段１４１から出力された第２の干渉波形信号ｓｆ’’（ｔｕ）を得て（Ｓ１６９）、ｓｆ’’（ｔｕ）をリサンプリング点算出手段１４８から得たリサンプリング点τνにおいてリサンプリングして、リサンプリング後の干渉波形ｓｆ’’’（ｔｕ）を得て、ｓｆ’’’（ｔｕ）を出力する（Ｓ１７０）。ただし、ｓｆ’’（ｔｕ）は離散的なデータであるため、サンプリング点τνの時刻にサンプリングしたデータがあるとは限らない（τν＝ｔｕとなるデータがあるとは限らない）。この場合は、補間が必要である。例えば、τνに近い２つのデータからの線形補間（Neugebauerの方程式による補間等）を使うことができる。
【００３４】
あるｚ＝ｚ０に対する第１の干渉波形信号（時間−強度の信号）ｓｆ（ｔ）の干渉波形が変化しなければ、上記ステップＳ１６８で求めた
【００３５】
【数８】
2014228397.tif 000010
【００３６】
の関係は変わらない。したがって、予め第１の干渉波形信号ｓｆ（ｔ）を測定し、算出した近似式の係数をもとに、新たに取得した第２の干渉波形信号ｓｆ’’（ｔｕ）をリサンプリングする。
【００３７】
図６に、波数変化特性を導出するための機能ブロック図を示す。図７を参照して、波数変化特性を導出する方法を説明する。図４に示した波数変化特性データ参照手段１５０による処理（Ｓ１７１）について、詳細に説明する。
【００３８】
離散フーリエ変換手段１４２は、干渉波形信号取得手段１４１から出力された第１の干渉波形信号ｓｆ（ｔｕ）を得て、ｓｆ（ｔｕ）を離散フーリエ変換して、
【００３９】
【数９】
2014228397.tif 000011
【００４０】
を計算して、Ｓｆ（ων）を出力する（Ｓ１６２）。ただし、
【００４１】
【数１０】
2014228397.tif 000012
【００４２】
は離散フーリエ変換処理を示し、ｆ＊ｇはｆとｇの畳み込み処理、ωνは角周波数、νは角周波数を識別するための整数、ｊは虚数単位（−１の平方根）を示す。k(tu)=k0t+kΩ(t)とすると、（Ｂ−５）式は、
【００４３】
【数１１】
2014228397.tif 000013
【００４４】
となる。
【００４５】
もし、ωνに２ｋ０ｚが含まれるようであれば、以下のようになる。
【００４６】
【数１２】
2014228397.tif 000014
【００４７】
マイナス波数領域除去手段１４３は、離散フーリエ変換手段１４２から出力された、離散フーリエ変換した干渉波形の信号（周波数−複素強度の信号）Ｓｆ（ων）を得て、Ｓｆ（ων）の中の第１項
【００４８】
【数１３】
2014228397.tif 000015
【００４９】
を０に置き換えた信号
【００５０】
【数１４】
2014228397.tif 000016
【００５１】
を算出して、Ｓｆ’（ων）を出力する（Ｓ１６３）。
【００５２】
上記の第１項を０にする第１の手段としては、Ｓｆ（ων）のマイナスの周波数成分を０とする方法がある。もし、第１項にプラスの周波数成分が無く、かつ、（Ｂ−８）式右辺にマイナスの成分が無い場合は、正確に第１項を０にすることができる。もし、そうでなければ、プラスの周波数は（Ｂ−８）式右辺と上記第１項とが混在し、また、マイナスの周波数でも同様であるので、マイナスの周波数を単純に０にすると、（Ｂ−８）式右辺の一部が削除され、上記第１項の一部が残存してしまい、歪みが生じる。ただし、このような混在がユーザの決めた所定の許容値以内であれば、Ｓｆ（ων）のマイナスの周波数を０にする方法は有用である。
【００５３】
上記の第１項を０にする第２の手段としては、Ｓｆ（ων）のマイナスの周波数領域の強度が０のウィンドウ関数を、Ｓｆ（ων）に掛けることが考えられる。例えば、非特許文献１にはflat-topウィンドウ関数を掛けることが記載されている。
【００５４】
離散フーリエ逆変換手段１４４は、マイナス波数領域除去手段１４３から出力された信号Ｓｆ’（ων）を得て、Ｓｆ’（ων）を離散フーリエ逆変換した信号
【００５５】
【数１５】
2014228397.tif 000017
【００５６】
を算出して、Ｓｆ’（ｔｕ）を出力する（Ｓ１６４）。
【００５７】
Ｅｒ（ｔｕ）Ｅ０（ｔｕ）の位相が２ｋ（ｔｕ）ｚに対して無視できるほど小さいのであれば、上記離散フーリエ逆変換後の信号（時間−複素強度）ｓｆ’（ｔｕ）の位相φｍ（ｔｕ）が波数ｋ（ｔｕ）に比例するので（φｍ（ｔｕ）＝２ｋ（ｔｕ）ｚ）、これを利用して、変換後の信号を使って、再度時間的にサンプリングする点（リサンプリングする点）を求める。そのために、下記の手順を行う。
【００５８】
位相算出手段１４５は、離散フーリエ逆変換手段１４４から出力された変換後の信号ｓｆ’（ｔｕ）を得て、ｓｆ’（ｔｕ）から位相φｍ（ｔｕ）を算出し、φｍ（ｔｕ）を出力する（Ｓ１６５）。
【００５９】
波数算出手段１４６は、位相算出手段２５から出力された位相φｍ（ｔｕ）を得て、φｍ（ｔｕ）を２ｚ０で割って波数ｋ（ｔｕ）を算出し、ｋ（ｔｕ）を出力する（Ｓ１６６）。
【００６０】
波数−時間関数ｎ次多項式近似手段１４７は、波数算出手段１４６から出力された波数ｋ（ｔｕ）を得て、ｋ（ｔｕ）からｋ（ｔ）の逆関数（波数を変数ｋとした時間の関数）ｔ（ｋ）を、ｎ次多項式
【００６１】
【数１６】
2014228397.tif 000018
【００６２】
に近似した近似式を導出する（Ｓ１６７）。例えば、非特許文献１では、３次の多項式
t'(k)=a0+a1k+a2k2+a3k3
で近似している。
【００６３】
なお、
【００６４】
【数１７】
2014228397.tif 000019
【００６５】
であることから、上記ステップＳ１６６の位相φｍ（ｔｕ）から波数ｋ（ｔｕ）を計算する手順は必須ではなく、上記ステップＳ１６７以降の波数ｋ（ｔｕ）を位相φｍ（ｔｕ）に置き換えても同等の結果となる。
【先行技術文献】
【非特許文献】
【００６６】
【非特許文献１】Shogo Yagi, Kazunori Naganuma, Tadayuki Imai, Yasuo Shibata, Shigeo Ishibashi, Yuzo Sasaki, Masahiro Sasaura, Kazuo Fujiura and Kazutoshi Kato, “A Mechanical-free 150-kHz Repetition Swept Light Source Incorporated a ＫＴＮ Electro-optic Deflector”, Proc. of SPIE Vol. 7889 78891J-1
【非特許文献２】Y. Yasuno, V. D. Madjarova, S. Makita, M. Akiba, A. Morosawa, C. Chong, T. Sakai, K. Chan, M. Itoh, and T. Yatagai, “Three-dimensional and high-speed swpt-source optical coherence tomography for in vivo investigation of human anterior eye segments,” Optics Express , Vol. 13, No. 26, pp. 10652-10664 (2005)
【発明の概要】
【発明が解決しようとする課題】
【００６７】
ＫＴＮ光源に用いられるＫＴＮ結晶は、ＡＣ電圧を印可することにより、結晶内に注入された電荷を制御し、屈折率分布を変えることにより結晶内を透過する光のビームを偏向させる。ＫＴＮ光源は、図１に示したように、ＫＴＮ光偏向器１０７から回折格子１１１への出射光の入射角を変化させることにより、発振する光の波数を変化させている。ＫＴＮ結晶への電荷注入は、結晶に一定のＤＣ電圧を印可することにより実施される。注入された電荷は時間とともに変化するので、これに伴って、透過する光のビームの偏向速度、偏向角などの偏向特性も経時変的に変化してしまう。波数掃引は、ＫＴＮ光偏向器による偏向により実現しているため、偏向特性が変化すると、波数変化特性も変化してしまう。つまり、ＫＴＮ光源は、時間経過に伴い、波数変化特性が変化（ドリフト）する。
【００６８】
上述したように、リサンプリングは、リファレンスである波数変化特性をもとに、干渉信号の時間間隔を補正する。逐次リサンプリングを行うためには、光源の波数変化が、時間経過したとしても、リファレンスとなる波数変化特性と同等の特性である必要がある。
【００６９】
図８に、ＫＴＮ光源の波数変化特性のドリフト現象を示す。ＫＴＮ光源は、上述した理由により、波数変化特性は時間とともにドリフトしてしまう。このため、リサンプリング処理の十分な効果を得ることができず、得られた干渉波形信号と、波数掃引した光源を用いた場合の干渉波形信号との誤差が大きくなる。この誤差のため、SS-OCTで得られたＰＳＦの広がりが大きくなり、得られた画像の精度が劣化するという問題があった。
【００７０】
本発明の目的は、ＫＴＮ光源の波数変化特性を考慮したリサンプリング処理を実現することにある。
【課題を解決するための手段】
【００７１】
このような目的を達成するために、一実施態様は、光源からの光の波長掃引を行って対象の物体に照射し、該物体からの反射光を、干渉計により干渉波形信号として取得する光干渉断層撮影法を適用したデータ処理装置におけるリサンプリング方法であって、前記データ処理装置の取得手段が、前記干渉波形信号と、前記光源の動作開始時からの経過時間とを取得する第１ステップと、前記データ処理装置の参照手段が、前記光源からの光の波数−時間関数を確定する第２ステップであって、前記第１ステップで取得した前記経過時間に基づいて、前記光源の動作開始時からの経過時間に応じた波数−時間関数が格納されたデータテーブルを参照する第２ステップと、前記データ処理装置の算出手段が、前記データテーブルを参照して得られた波数−時間関数の波数が一定間隔となる、前記干渉波形信号のリサンプリング点を求める第３ステップと、前記データ処理装置のリサンプリング手段が、前記干渉波形信号を、前記リサンプリング点でリサンプリングする第４ステップとを備えたことを特徴とする。
【発明の効果】
【００７２】
以上説明したように、光源の動作開始時からの経過時間に応じた波数−時間関数を正確に把握できるため、リサンプリングしても波数掃引した光源を用いた場合と同等な干渉波形信号となる。従って、SS-OCTで得られた信号のＰＳＦの広がりが小さくなり、SS-OCTで得られた画像が鮮明となる。
【図面の簡単な説明】
【００７３】
【図１】従来のSS-OCT光源に用いられるＫＴＮ光源の構成を示す図である。
【図２】従来のSS-OCTにおける光学系の基本構成を示す図である。
【図３】干渉波形信号のリサンプリング処理を説明するための図である。
【図４】従来の従来のデータ処理装置におけるリサンプリング処理のための機能ブロック図である。
【図５】従来のリサンプリング方法を示すフローチャートである。
【図６】波数変化特性を導出するための機能ブロック図である。
【図７】波数変化特性を導出する方法を示すフローチャートである。
【図８】ＫＴＮ光源の波数変化特性のドリフト現象を示す図である。
【図９】本発明の一実施形態にかかるリサンプリング処理のための機能ブロック図である。
【図１０】本発明の一実施形態にかかるリサンプリング方法を示すフローチャートである。
【図１１】ＫＴＮ光源の波数変化特性のデータテーブルの一例を示す図である。
【図１２】本発明の一実施形態にかかるリサンプリング処理の結果を示す図である。
【発明を実施するための形態】
【００７４】
以下、図面を参照しながら本発明の実施形態について詳細に説明する。ＫＴＮ光源の動作時における波数変化特性は経時変化するが、温度等の環境条件が同じで、かつ、光源として動作を開始した後の経過時間が同じであれば、波数変化特性に伴うドリフトも同じように再現する。そこで、ＫＴＮ光源の動作開始時からの経過時間に応じて波数変化特性を予め測定し、データテーブルとし格納しておく。SS-OCTによる画像解析を行う場合には、ＫＴＮ光源の動作開始時からの経過時間に応じてデータテーブルを参照し、参照された波数変化特性をもとにリサンプリング処理を行う。波数変化特性が経時変化するＫＴＮ光源であったとしても、リサンプリング処理の効果を十分に得ることができ、SS-OCTで得られた信号のＰＳＦの広がりが小さくなり、SS-OCTで得られた画像が鮮明となる。
【００７５】
図９に、本発明の一実施形態にかかるリサンプリング処理のための機能ブロック図を示す。図１０を参照して、本発明の一実施形態にかかる干渉波形信号の変換（リサンプリング）方法を説明する。
【００７６】
干渉波形信号取得手段２０１は、図２に示した光学系の受光素子１２６から得られた干渉波形信号（時間−強度の信号）を取得して、記憶手段に格納しておく（Ｓ２２１）。このとき、基準面から距離ｚ０（既知）の位置に反射面（例えば鏡面）を持つ測定対象物１３１に対して、その干渉波形信号（時間−強度の信号）は、
【００７７】
【数１８】
2014228397.tif 000020
【００７８】
となる。ただし、ｔｕは時間であり離散値であり、ｋ（ｔｕ）は波数である。また、ｕはその離散値を識別する整数であり、第１の干渉波形信号がＮ個ある場合は、ｕは１〜Ｎの値を取る。このとき、干渉波形信号取得手段２０１は、ＫＴＮ光源１２１の動作開始時からの経過時間（Ｔ）も取得しておく。
【００７９】
データテーブル参照手段２１０は、測定された経過時間（Ｔ）に基づいて、予め記憶手段に格納されている波数変化特性（波数−時間関数）を取得する（Ｓ２３１）。
【００８０】
リサンプリング点算出手段２０８は、データテーブル参照手段２１０から得られた波数変化特性（波数−時間特性）を、ｎ次多項式近似により導出したｋ（ｔ）の逆関数ｔ（ｋ）の近似関数ｔ’（ｋ）を使って、波数が一定間隔となる波数値ｋν（νは波数値同士を区別する整数）を選び（ｋν−ｋν＋１＝Δｋ; Δｋは一定値）、
【００８１】
【数１９】
2014228397.tif 000021
【００８２】
を干渉波形ｓｆ（ｔｕ）のリサンプリング点として出力する（Ｓ２２８）。
【００８３】
リサンプリング手段２０９は、干渉波形信号取得手段２０１から出力された第２の干渉波形信号ｓｆ’’（ｔｕ）を得て（Ｓ２２９）、ｓｆ’’（ｔｕ）をリサンプリング点算出手段１０８から得たリサンプリング点τνにおいてリサンプリングして、リサンプリング後の干渉波形ｓｆ’’’（ｔｕ）を得て、ｓｆ’’’（ｔｕ）を出力する（Ｓ２３０）。ただし、ｓｆ’’（ｔｕ）は離散的なデータであるため、サンプリング点τνの時刻にサンプリングしたデータがあるとは限らない（τν＝ｔｕとなるデータがあるとは限らない）。この場合は、補間が必要である。例えば、τνに近い２つのデータからの線形補間（Neugebauerの方程式による補間等）を使うことができる。
【００８４】
ステップＳ２２１においては、干渉波形信号取得手段２０１が、取得した干渉信号波形を記憶手段に格納したが、記憶手段には格納せずに、干渉波形信号をリサンプリング手段２０９に出力してもよい。この構成によれば、記憶手段が不要となるので、データ処理装置のコスト削減、SS-OCT画像解析システムの小型化を図ることができる。
【００８５】
SS-OCT光源から出力される光の波数−時間関数ｔ（ｋ，ａ０，．．．，ａｎ）は、ＫＴＮ光源の波数変化特性の経時変化を考慮した波数−時間関数であることから、リサンプリングした干渉波形信号と、波数掃引した光源を用いた干渉波形信号との誤差は、わずかである。この場合には、誤差が小さくなることから、SS-OCTで得られた信号のＰＳＦの広がりが小さくなり、SS-OCTで得られた画像が鮮明になる。
【００８６】
図１１に、ＫＴＮ光源の波数変化特性のデータテーブルの一例を示す。ＫＴＮ光源の動作開始時からの経過時間０秒、１０秒、および３０秒ごとの波数変化特性を示している。この波数変化特性を、所定の時間間隔ごとの相対波数としてデータテーブルを予め作成し、記憶手段に格納しておく。
【００８７】
図１２に、本発明の一実施形態にかかるリサンプリング処理の結果を示す。図１２（ａ）は、ＫＴＮ光源の波数変化特性のドリフトを考慮せずに、リサンプリング処理を行った場合のＰＳＦを示している。図１２（ｂ）は、本実施形態にかかるリサンプリング処理を行った場合のＰＳＦを示している。上述したように、従来の方法では、SS-OCTで得られた信号のＰＳＦの広がりが大きく、明確なピークが認められない。これを画像として見ると、２つに分裂した、ぼやけた画像となる。
【００８８】
一方、本実施形態にかかるリサンプリング処理の場合、ＰＳＦの広がりは小さくなっており、シャープなピークが見られる。これを画像として見ると、コントラストの高い画像となる。
【００８９】
なお、ＫＴＮ光源に用いる光偏向器として、ＫＴＮ結晶の一部の元素を置換したＫＬＴＮ結晶（Ｋ１−ｙＬｉｙＴａ１−ｘＮｂｘＯ３（０＜ｘ＜１、０＜ｙ＜１））を用いた光偏向器も、SS-OCTの光学系に適用することができる。
【符号の説明】
【００９０】
１０１  光ファイバ
１０２  アイソレータ
１０３  集光レンズ
１０４  ＳＯＡ（semiconductor optical amplifier）
１０５  反射面
１０６  コリメータレンズ
１０７  ＫＴＮ光偏向器
１０８，１０９  電極
１１０  電源
１１１  回折格子
１２１  ＫＴＮ光源
１２２，１２７  光カプラ
１２３，１２８  サーキュレータ
１２４，１３０  レンズ
１２５  参照ミラー
１２６  受光素子
１２９  駆動ミラー
１３１  測定対象物
１３２  データ処理装置
１３３  ドライバ回路
【図１】
2014228397.tif 000022
【図２】
2014228397.tif 000023
【図３】
2014228397.tif 000024
【図４】
2014228397.tif 000025
【図５】
2014228397.tif 000026
【図６】
2014228397.tif 000027
【図７】
2014228397.tif 000028
【図８】
2014228397.tif 000029
【図９】
2014228397.tif 000030
【図１０】
2014228397.tif 000031
【図１１】
2014228397.tif 000032
【図１２】
2014228397.tif 000033