(57)【要約】【課題】  個別周波数制御走査レーザーの各周波数を高精度で設定する方法と装置。【解決手段】  個別周波数制御走査レーザーを走査させたままで制御電流を変更できる機構を持っている走査制御回路を用いて、個別周波数制御走査レーザーを連続走査し、光干渉計の干渉信号を干渉信号検出器により計測することによって、前記個別周波数制御走査レーザーの光周波数を高感度で測定し、走査した状態で走査制御回路により個別周波数制御走査レーザーの各周波数を高精度で設定できる。【選択図】図１
2015001452.tif 000002
【特許請求の範囲】
【請求項１】
走査制御回路によって走査が制御されている個別周波数制御走査レーザーの出力光が入力される光干渉計と、前記光干渉計が出力する干渉信号を検出する干渉信号検出器とからなるシステムにおいて、前記個別周波数制御走査レーザーを走査させた状態で前記干渉信号検出器が検出する信号が所望の形の干渉信号になるように前記走査制御回路を用いて前記個別周波数制御走査レーザーの制御電流の設定を可能にする方法。
【請求項２】
走査制御回路によって走査が制御されている個別周波数制御走査レーザーの出力光が入力される光干渉計と、前記光干渉計が出力する干渉信号を検出する干渉信号検出器とから構成されていて、前記個別周波数制御走査レーザーを走査させた状態で前記干渉信号検出器が検出する信号が所望の形の干渉信号になるように前記走査制御回路を用いて前記個別周波数制御走査レーザーの制御電流の設定を可能にする装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、周波数走査レーザーの各走査周波数を精密に設定する装置および方法に関するものである。
【背景技術】
【０００２】
オプティカル・コヒーレンス・トモグラフィー（英語はOptical coherence tomographyで、以降、ＯＣＴと略称する）は、赤外光線を被検体に照射し、後方反射または後方散乱された（以降、後方反射に後方散乱も含めるものとする）信号光を、参照光と干渉させて、照射光線に沿って、後方反射光強度の被検体内での奥行方向の距離依存性を測定し、さらに、照射光線を横方向に走査することにより、被検体の２次元や３次元断層画像を撮像する方法である。ＯＣＴは、プローブを被検体に接触させる必要がないので非接触であり、赤外線はX線とは異なり被験者の被爆が無く、規格内の照射強度で使用すれば被検体に影響を残さず無侵襲である。X線ＣＴ（コンピューテッド・トモグラフィー）、ＭＲＩ（磁気共鳴映像法）、超音波に比べてより高分解能で高速の測定ができる。ＯＣＴは、眼科、冠状動脈内壁、皮膚などの診断装置として既に臨床応用されており、歯科や内臓の癌などの診断装置や、工業生産品の検査装置等の実用化も進められていて、産業上の幅広い応用が期待される断層画像撮像方法である。
【０００３】
ＯＣＴの方法には、時間領域（英語はTime domainで、以降、ＴＤと略称する）ＯＣＴ、スペクトル領域（英語はSpectral domainで、以降、ＳＤと略称する）ＯＣＴ、周波数走査（英語はSwept  sourceで、以降、ＳＳと略称する）ＯＣＴの３つの方法がある。３つのＯＣＴの方法を比べると、同じ計測時間であれば、ＳＤ−ＯＣＴやＳＳ−ＯＣＴの方がＴＤ−ＯＣＴに比べて感度が高い。また、ＳＤ−ＯＣＴとＳＳ−ＯＣＴを比べると、ＳＳ−ＯＣＴの方が計測可能距離を長くできる。したがって、ＳＳ−ＯＣＴが多くの応用分野で用いられるようになってきている。
【０００４】
ＳＳ−ＯＣＴの構成図の一例を図３（Ａ）に示す。周波数走査光源５０１は、出力光を光ファイバ５０３に出力し、サンプル開始トリガー電気信号５４１とサンプリングクロック電気信号５４２をコンピュータ５１８に接続されたＡ／Ｄ変換器５１７に出力する。周波数走査光源５０１には、通常、周波数走査レーザーが用いられる。
【０００５】
周波数走査光源５０１からの出力光はＯＣＴ干渉計５０２に導かれる。図３（A）にはマッハ・ツェンダー干渉計を用いたＯＣＴ干渉計が例示されているが、マイケルソン干渉計を用いる場合も多い。ＯＣＴ干渉計５０２に導かれた周波数走査光源からの出力光は、分波カプラ５０５に導かれ、分割され、参照光路５３１と試料光路５３０に導かれる。両光路とも、分波カプラ５０５から始まり、合波カプラ５１１で終了する。参照光路の光はサーキュレータ５１２に入力ポート１から入力し、入出力ポート２から出力してコリメータ５１３とレンズ５１４を用いて参照ミラー５１５に照射され、参照ミラー５１５により反射された光は照射光路を逆に進み、サーキュレータ５１２の入出力ポート２に入力し出力ポート３に出力される。サーキュレータ５１２の出力ポート３から出力された光は、合波カプラ５１１に導かれる。
【０００６】
分波カプラ５０５から試料光路に導かれた光は、サーキュレータ５０６の入力ポート１に導かれ、サーキュレータ５０６の入出力ポート２から出力し、コリメータ５０７、ガルバノミラー５０８、レンズ５０９を介して試料５１０に照射される。試料５１０から後方反射または後方散乱された光（以降、後方反射に両方を含める）は、照射光路を逆に進み、サーキュレータ５０６の入出力ポート２に入力し、出力ポート３から出力して合波カプラ５１１に導かれる。ガルバノミラー５０８によって、試料への光の照射位置を移動し、２次元や３次元の断層画像を構成するため必要とされるデータを得る。
【０００７】
合波カプラ５１１に導かれた参照光路からの光と試料光路からの光は干渉し、干渉信号を出力する。出力された干渉信号は差動光検出器５１６で検出され電気信号に変換される。差動光検出器５１６で検出された電気信号はＡ／Ｄ変換器５１７でデジタル信号に変換され、コンピュータ５１８に収録される。コンピュータは、取り込まれた信号からOCT断層画像５１９を構成し、表示器５２０に表示する。
【０００８】
測定された干渉信号からOCT画像を構成する方法について説明する。周波数走査光源５０１は、図３（Ｂ）に示すように周波数走査を繰り返し、各周波数走査開始時にサンプル開始トリガー信号を出力し、等周波数間隔δfで周波数f1からfNまで走査される。走査する各周波数に同期してサンプリングリングができるように、サンプリングクロック信号を出力する。このような走査を連続して繰り返す。図３（A)において、コンピュータ５１８に接続されたＡ／Ｄ変換器５１７は、サンプル開始トリガー信号に同期して一連のサンプリングを開始し、サンプリングクロック信号に同期して各周波数における干渉信号をサンプリングする。1つのサンプル開始トリガー信号に続く一連の周波数f1からfNまでのサンプリングで、１つの奥行方向を走査する信号が得られる。この一つの奥行方向の走査をA走査と呼ぶ。
【０００９】
図３（Ａ）において、差動光検出器５１６で信号を検出すると、直流成分は差し引かれ、i番目の周波数fiおいて、次の式で与えられる干渉信号Siが検出される。
【数１】
2015001452.tif 000003
ここで、Cは定数、r(z)は照射光線に沿った位置zにおける試料の後方反射率、cは光速である。試料が、位置z0にある反射率r0の鏡だけの場合は、r(z)はδ関数を用いてr(z)=r0δ（z−z0）と表せる。この場合、数１式をzについて積分すると、次式のように単一の正弦波が得られる。
【数２】
2015001452.tif 000004
jは純虚数であり、j2=−1である。数２式に、次式のように離散的フーリエ変換を施す。
【数３】
2015001452.tif 000005
このフーリエ変換Di(Z)において、右辺括弧内の第一項の指数関数はZ=z0のとき１となり、加算するとNCr0/2になる。このとき、右辺第二項は、fiの変化ととともに激しく変わるプラスとマイナスの値をとり、加算したものは小さな値になる。従って、Di(Z)はZ=z0のとき実質的に値NCr0/2をとる。同様に、Di(Z)はZ＝−z0でも値NCr0/2になる。その他のZの値では、加算される指数関数がプラスの値とマイナスの値の間で激しく変化するので、加算した値はNCr0/2に比べて非常に小さな値になる。このように、Di(Z)はZ＝z0とZ＝−z0でピークを持ち、２つピークから離れると急激に小さな値になる。この結果は、観測された数２式の干渉信号に、位置Zを変数として数３式のフーリエ変換を繰り返して行い、その結果をZの関数としてプロットすれば、位置Z＝z0に反射率r0の反射面があるということが分かることを示している。（通常、Z>0の位置に試料を置く。）複数の反射率r0、 r1、 ・・・rMの反射面が、それぞれ試料内の位置z0、 z1、 ・・・zMにあれば、干渉信号SiはCrmcos(4πfizm/c)の形の式の和になり、フーリエ変換  Di(Z)は、Z＝z0、 z1、・・・zMでそれぞれ反射率r0、 r1、 ・・・rMに比例した振幅のピークを持つ。つまり、図３（Ａ）に示すようなＳＳ−ＯＣＴ装置で観測される干渉信号にフーリエ変換を行えば、試料内のどの位置にどのような反射率の面があるかが分かることになる。生体のように連続的に反射面があれば、試料の奥行方向の位置の関数として、反射率がZの関数として連続的に求められることになる。
【００１０】
通常は、数３式の絶対値の２乗をとったスペクトル強度をOCT画像表示に用いる。奥行方向Zの関数としてこのスペクトル強度を表したものをA走査信号として用いる。図３（Ａ）のガルバノミラー５０８を用い、試料上で光を横に走査してA走査を繰り返す走査をB走査と呼ぶ。B走査により得られたA走査信号を集め、２次元上に試料の奥行方向の位置と横走査方向の位置の関数としてスペクトル強度を明るさの濃淡に対応させて表した図が２次元ＯＣＴ断層画像である。さらに、ガルバノミラー５０８を用いてB走査の位置を移動させながら取得した２次元断層画像を集めて画像処理を行えば、３次元ＯＣＴ断層画像が構成できる。
【００１１】
ＳＳ−ＯＣＴの性能は光源に用いる周波数走査レーザーの性能に大きく依存するため、より高性能のＳＳ−ＯＣＴを実現するための周波数走査レーザーの開発が進められてきた。ＯＣＴの性能には、高速で高分解能で深い計測が求められる。超周期回折格子分布反射型レーザー（英語はSuper-structure grating distributed Bragg reflector laser、以降、ＳＳＧ−ＤＢＲレーザーとよぶ略称する）がＳＳ−ＯＣＴの周波数走査レーザー光源に利用できるという発明がなされた（特許文献１、特許文献２）。このレーザーは、数ナノ秒で周波数をスイッチでき、高速走査が可能である。また、複数のレーザーを組み合わせて走査することによって周波数走査幅を広くでき、高分解能が実現できる。更に、瞬間的なコヒーレンス長は数mにでき、深いＯＣＴ計測が可能である。このレーザーは、走査する各周波数を、それぞれ個別に制御電流を用いて設定でき、周波数を図３（Ｂ）に示す通りに、時間の関数としてステップ状（階段状）に周波数等間隔で走査できるという他の周波数走査レーザーにはない長所を持っている。等間隔の時間で正確に周波数が走査できるので、図３（Ａ）において、サンプリングクロックを用いずに、A/D変換器の内部クロックでも各データをサンプリングすることもできる
【００１２】
ＳＳＧ−ＤＢＲレーザーの概略図を図４（Ａ）に示す（非特許文献１）。このレーザーは、前方ＳＳＧ−ＤＢＲ領域２０３と後方ＳＳＧ−ＤＢＲ領域２００の２つの反射体を向かい合わせたファブリ・ペロー干渉計のキャビティ構造をしている。ＳＳＧ−ＤＢＲは、注入する制御電流に依存して、複数の光周波数において選択的に高い反射率を示す。前方と後方のＳＳＧ−ＤＢＲが反射する光の周波数が一致した周波数でレーザーは発振する。レーザーの出力周波数は、前方ＳＳＧ−ＤＢＲ領域２０３への制御電流（以降、前方ＳＳＧ−ＤＢＲ電流と呼ぶ）と後方ＳＳＧ−ＤＢＲ領域への制御電流（以降、後方ＳＳＧ−ＤＢＲ電流と呼ぶ）の２つの制御電流で粗調整し、位相領域２０１への制御電流（以降、位相電流と呼ぶ）で微調整する。ゲイン領域２０２に流す制御電流（以降、ゲイン電流と呼ぶ）を増加すると、レーザー出力が増加する。半導体光増幅器（英語はSemiconductor optical amplifier で、以降、略してＳＯＡと呼ぶ）領域２０４によって、外部に取り出す光の強度を増加させるとともに、ＳＯＡ領域の制御電流（以降、ＳＯＡ電流と呼ぶ）によって、周波数によらずに出力強度を一定にすることができる。
【００１３】
ＳＳＧ−ＤＢＲレーザーと類似のレーザーとして、図４（Ｂ）〜（Ｅ）に示す４種類が報告されている。図４（Ｂ）は、Sampled  Grating (ＳＧ) ＤＢＲレーザーで、ＳＧ−ＤＢＲがＳＳＧ−ＤＢＲに対応し、ＳＳＧ−ＤＢＲレーザーと類似の動作をする（非特許文献２）。図４（Ｃ）は、Digital Supermode (ＤＳ) ＤＢＲレーザーで、前方ＤＢＲ領域２２３が複数の部分に分けられ、独立した複数の制御電流で制御されている（非特許文献３）。ＤＳ−ＤＢＲレーザーの後方ＤＢＲ領域２２０と前方ＤＢＲ領域２２３が、それぞれ、ＳＳＧ−ＤＢＲレーザーの後方ＳＳＧ−ＤＢＲ領域２００と前方ＳＳＧ−ＤＢＲ領域２０３に対応し、ＤＳ−ＤＢＲレーザーの動作原理は、ＳＳＧ−ＤＢＲレーザーと類似である。図４（Ｄ）は、英語名はVertical Grating Assisted Co-directional Coupler with Sampled  Grating Reflector Laserで、略してＧＣＳＲレーザーと呼ばれる（非特許文献４）。構造は、前方ミラーにＤＢＲを用いずに、回折格子カプラ２３２を用い、後方ミラーにStructured（Ｓ）ＤＢＲを用いている。回折格子カプラ２３２は反射器ではなくて、周波数を選択的に透過する役割をする。ゲイン領域２３３の出力部にレーザーの前方ミラーの処理が施されている。しかし、制御方法は、ＧＣＳＲレーザーの回折格子カプラ２３２の制御電流と後方Ｓ−ＤＢＲの制御電流を、それぞれ、ＳＳＧ−ＤＢＲレーザーの前方ＳＳＧ−ＤＢＲ電流と後方ＳＳＧ−ＤＢＲ電流に対応させれば、ＳＳＧ−ＤＢＲレーザーと類似である。図４（Ｅ）は、Modulated Grating Y-branch (ＭＧＹ) レーザーの概略図である（非特許文献５）。このレーザーは、２つの特性の異なるＤＢＲ−Ａミラー２４０とＤＢＲ−Ｂミラー２４１をマルチモード干渉計カプラ２４２で結合した構造をしている。ゲイン領域２４４に前方ミラーの処理が施されている。MGYレーザーの制御方法は、ＤＢＲ−Ａの制御電流とＤＢＲ−Ｂの制御電流を、それぞれ、ＳＳＧ−ＤＢＲレーザーの前方ＳＳＧ−ＤＢＲ電流と後方ＳＳＧ−ＤＢＲ電流に対応させれば、ＳＳＧ−ＤＢＲレーザーと類似である。
【００１４】
ＳＳＧ−ＤＢＲレーザーにおいて、後方ＳＳＧ−ＤＢＲと前方ＳＳＧ−ＤＢＲが反射する光の周波数の電流依存性は異なっており、ノギスにおいて目盛間隔の異なる主尺と副尺（バーニア：Vernier）の一致した位置で正確に長い距離の測定を可能にするように、周波数の電流依存性の異なる２つのＳＳＧ−ＤＢＲの組み合わせで、幅広い周波数範囲の発振を可能にしている。これはバーニア効果と呼ばれる。図４に示す５つのレーザーは、どれも2つ以上の制御電流を用い、バーニア効果を利用してレーザーの発振周波数を粗調整で広帯域に走査し、位相電流を用いて発振周波数の微調整を行うことによって、正確な広帯域の周波数走査を可能にしているという共通点がある。ＳＳ−ＯＣＴの観点からすると、これらのレーザーは、ＳＳ−ＯＣＴに用いられている他の周波数走査レーザーと異なり、制御電流を用いて走査する各周波数を個別にそれぞれ制御するという点に特徴があるので、この発明では、図４に示す５つのレーザーを「個別周波数制御走査レーザー」と総称することにする。現状報告されているこの種のレーザーは図４に示す５種であるが、今後、これら以外のレーザーも同じ技術的思想に基づいて開発される可能性があるが、それらも個別周波数制御走査レーザーに含まれるものとする。個別周波数制御走査レーザーは、複数の制御電流を用いて、走査する周波数を個別にすべて制御しなければならないという複雑さがあるが、精密な周波数制御方法を開発することによって、個別の周波数をすべて正確に設定できるという利点を有する。
【００１５】
個別周波数制御走査レーザーのＯＣＴ撮像への有効性は、ＳＳＧ−ＤＢＲレーザーを用いて実証された（非特許文献６、非特許文献７）。24mmの深さのＯＣＴ計測が可能であることも実証された（非特許文献８）。しかし、従来技術では、ＳＳＧ−ＤＢＲレーザーを用いたＳＳ−ＯＣＴ撮像の画像の濃淡のダイナミックレンジが実際には悪く、高画質の撮像ができていないという問題があった。
【００１６】
ＳＳ−ＯＣＴ画像の濃淡のダイナミックレンジが悪くなる原因について説明する。一枚のミラーを試料の位置に置いたときの離散的フーリエ変換のスペクトル（数３式）の絶対値の二乗をとったスペクトル強度は、ポイント・スプレッド・ファンクション（英語はPoint Spread Functionで、以降ＰＳＦと略称する）と呼ばれる。ＳＳ−ＯＣＴの画像の濃淡のダイナミックレンジは、ＰＳＦのダイナミックレンジで評価できる。ＰＳＦのダイナミックレンジは、ＰＳＦのピークの強度とノイズフロアのレベルとの比をdBで表したもので定義される。数３式に関係する個別周波数制御走査レーザーの因子は、各出力における強度と周波数である。これら因子の理想値からのずれが、ＰＳＦのダイナミックレンジを劣化させる。これらの影響を、数値計算でシミュレーションする。ＳＳ−ＯＣＴでは、ＰＳＦのダイナミックレンジを向上させるために、通常、数３式に窓関数を用いて離散的フーリエ変換が行われる。代表的な窓関数は以下に示すハニング窓である。
【数４】
2015001452.tif 000006
ＳＳ−ＯＣＴの数値処理では、他の窓関数も用いられるが、この発明の有効性は窓関数には依存しない。数値計算シミュレーションでは、数４式の窓関数を用い、以下に示す離散的フーリエ変換を行った。
【数５】
2015001452.tif 000007
個別周波数制御走査レーザーの出力強度は、出力周波数に依存しないようにＳＯＡ電流を変化させて一定化できるが、制御回路の性能の限界により、0.1％〜0.5％程度の変動は避けられない。レーザー強度の変動の影響は、数１式の係数Cが出力番号iに依存し、Ci=Ca(1+δci)と変動するものとして表せる。ここで、Caは平均値、δciは変動の割合である。周波数の変動はないものとして、δciの標準偏差が、0.001、0.002、0.003、0.004、0.005のときのＰＳＦを、数５式を用いて数値計算し、ダイナミックレンジを求めると、それぞれ、87dB、80dB、78dB、75dB、73dBになる。このダイナミックレンジは深さZに依存しない。現状の技術によれば、強度の変動δciの標準偏差を0.005（0.5％）以下に抑えるのは可能であり、周波数の変動がなければ、70dB以上のダイナミックレンジの実現が可能になるはずである。
【００１７】
個別周波数制御走査レーザー用いた場合、出力する各周波数に誤差があると、どのようにＰＳＦのダイナミックレンジが劣化するかを、数値計算シミュレーションで評価する。出力周波数に誤差が無い時のi番目の周波数fiを次式で表す。
【数６】
2015001452.tif 000008
この式は、誤差がない理想的な場合は、初期周波数f1=f0+δfから最終の周波数fN=f0+Nδfまで、等周波数間隔δfで走査されることを示している。しかし、実際に走査されるi番目の周波数Fiには誤差Δiがあり、正確な周波数fiからずれている。この場合、実際に走査される周波数Fiは次式で表される。
【数７】
2015001452.tif 000009
図５に、i=1から５までの場合について、Fi、fi、Δiの関係を例示する。tsは、隣接する周波数走査の時間間隔である。干渉信号の数２式は、次式で表される。
【数８】
2015001452.tif 000010
数８式において、4πΔiz0/cは余弦関数の位相の一部であるが、これがランダムに変動するので、余弦関数は乱れたものになる。深さが深く、z0が大きいほど、位相の変動が大きくなり、余弦関数の乱れも大きくなる。余弦関数の乱れが大きくなると、ＰＳＦのダイナミックレンジは小さくなる。
【００１８】
具体例として、周波数走査の開始周波数f1＝182.50 THz、終了周波数fN＝192.20
THz、周波数間隔δf＝6.25 GHzの場合について、ＰＳＦを数値計算する。走査する周波数の数はN＝1552とする。この例の周波数走査の中心波長は１．６μmで、光通信の波長帯に含まれていて、光計測の検査機器がそろっているとともに、ＯＣＴ計測に適した波長帯の一つである。数５式の離散的フーリエ変換の指数関数の周期性から、計測可能距離Δzは次式で与えられる。
【数９】
2015001452.tif 000011
周波数間隔をδf = 6.25GHzとすると、計測可能距離はΔz = 12mmになる。
【００１９】
周波数の誤差Δiがガウス分布でランダムに変動するようにし、その標準偏差σfの値を変えて、３つのミラーの位置（z0=0.5mm、z0=6mm、z0=11.5mm）について計算したＰＳＦを図６に示す。標準偏差σfが大きいほど、周波数の設定精度が悪いことを意味する。周波数変動の影響を見るために、強度変動はないものとした。各グラフの横軸に示す測定位置の両端に対応する位置z0=0mmとz0=12mmでは、ピークが半分欠けるので、代表値として用いるのを避けた。各グラフの縦軸のスペクトル強度はdB（デシベル）で表し、z0=0mmでのスペクトル強度を0dBとして規格化した。図６（Ｄ）に例示するように、ＰＳＦはピーク強度とノイズフロア（図６中では記号ＮＦで示す）で特徴づけられる。図６（Ｄ）の場合、ピークは0dBでノイズフロアは−100dBである。それらの差の100dBが、ＰＳＦのダイナミックレンジ（図６中では記号ＤＲで示す）を表す。図６に示すように、ＰＳＦのダイナミックレンジは、個別周波数制御走査レーザーの設定周波数の精度に大きく依存する。ＯＣＴ計測の場合、ダイナミックレンジは、どの深さにおいても40dB以上が望ましく、画質はダイナミックレンジが大きくなるとともに向上するので、さらにより大きなダイナミックレンジが望ましい。周波数の誤差の標準偏差σfが０の場合は、図６（Ａ）〜（Ｃ）に示すように、奥行方向の位置z0が0.5から11.5mmの範囲のどの位置においても、ダイナミックレンジは160dB以上である。周波数の誤差の標準偏差がσf =0.01GHzとσf =0.05GHzの場合、それぞれ図６（Ｄ）〜（Ｆ）と図６（Ｇ）〜（Ｉ）に示すように、奥行方向の位置z0が0.5から11.5mmの範囲のどの位置においても、ダイナミックレンジは60dB以上である。これに対し、周波数の誤差の標準偏差がσf =0.5GHzまで劣化すると、図６（Ｊ）〜（Ｌ）に示すようにノイズフロアの増加が大きく、ダイナミックレンジも、z0=6mm以上では50dBよりも悪くなり、画質の劣化が起きる。周波数の誤差の標準偏差がさらに劣化し、σf =5GHzになると、図６（Ｍ）〜（Ｏ）に示すように、z0=0.5mm以上では、ダインミックレンジは50dB以下であり、z0=6mmでは21dB、z0=11.5mmでは4dBしかなく、OCT画像の劣化は著しく、z0の大きい深い位置では、OCT撮像そのものが困難になる。以上のように、光源の周波数を走査するＳＳ−ＯＣＴにおいては、設定する周波数の精度が極めて重要である。
【００２０】
図７に、３つの深さ位置における、ダイナミックレンジ（ＤＲ）と周波数の誤差の標準偏差σfの関係をまとめて示す。図７の中の曲線（Ａ）、（Ｂ）、（Ｃ）は、それぞれ、深さ位置z0=3mm、z0=6mm、z0=12mmにおける関係である。必要なダイナミックレンジを60dBとすると、ダイナミックレンジが60dBの値で引いた点線と各曲線との交点から実現しなければならない周波数の誤差の標準偏差が求まる。深さ位置3mm、6mm、12mmでDR＞60dBを実現するためには、周波数の誤差の標準偏差σfは、それぞれσf＜0.24GHz、σf＜0.12GHz、σf＜0.06GHzでなければならない。必要なダイナミックレンジを50dB以上とすると、必要な周波数の誤差の標準偏差σfは、深さ位置3mm、6mm、12mmに対し、σf＜0.78GHz、σf＜0.38GHz、σf＜0.18GHzでなければならない。必要なダイナミックレンジを40dB以上とすると、必要な周波数の誤差の標準偏差σfは、深さ位置3mm、6mm、12mmに対し、σf＜2.5GHz、σf＜1.2GHz、σf＜0.59GHzでなければならない。
【００２１】
以上のように、例えばＰＳＦのダイナミックレンジが少なくとも50dB以上の、濃淡の画質の比較的に良いOCT画像を、深さ位置12mmまで撮像するためには、周波数の誤差の標準偏差が0.2GHz以下になるように周波数調整がなされなければならないが、個別周波数制御走査レーザーの従来技術では、このような精度の周波数調整はできていなかった。従来技術によって実現された周波数の誤差の標準偏差についてＳＳＧ−ＤＢＲを例として説明する。発振周波数は、前後の２つのＳＳＧ−ＤＢＲ電流で粗調整され、位相電流で連続的に微調整される。図４（A)に示すＳＳＧ−ＤＢＲレーザーで、まず位相電流を零にし、ゲイン電流とSOA電流を一定値に設定し、前方ＳＳＧ−ＤＢＲ御電流と後方ＳＳＧ−ＤＢＲ電流を変化させ、各値でのレーザーの発振周波数を光スペクトルメーターで測定する。横軸に前方ＳＳＧ−ＤＢＲ電流の値、縦軸に後方ＳＳＧ−ＤＢＲ電流の値を取って、レーザー発振周波数の等高線を描くと、図８に細い実線で示す等高線マップが得られる。この等高線マップの数値の詳細は、素子によって異なるが、ここで議論する一般的な振る舞いは、素子に依存しない。図８において、斜め左下から右上に向かって発振周波数は増加する。太い実線は、レーザーの発振モードが異なる領域の境界を示す。各領域の白黒の濃淡は、濃い領域ほど低い周波数に対応している。太い実線の境界で、周波数の等高線は不連続になっているので、破線で示すような、同一モード領域内の線に沿って前方ＳＳＧ−ＤＢＲ電流および後方ＳＳＧ−ＤＢＲ電流を対で変化させて、周波数を走査する。領域内での破線の位置は、位相電流を調整してスペクトル純度が高くなるように設定する。レーザーの発振は、主スペクトルに対してサイドモードが弱い方がスペクトルの純度が高い。このスペクトル純度は、サイドモード・サプレッション・レシオ（英語はSide Mode Suppression Ratioで、以降、ＳＭＳＲと略称する）で判定される。ＳＭＳＲが極大になるところを、位相電流をわずかずつ変化させて求める。図８の破線はこの手順で決められた。この破線に沿って３つの電流を走査すれば、一つのモード領域内では、モードの飛びが無い走査が得られる。幅広い走査をするためには、１から８までの番号を付した破線に沿って制御電流値を順次走査する。このようにして、ＳＳＧ−ＤＢＲレーザーを走査するための制御電流は決定されてきた。これらの手順は、例えば非特許文献９や非特許文献１０で説明されている。
【００２２】
上述の手順で決められた制御電流で個別周波数制御走査レーザーを走査した場合、濃淡のダイナミックレンジが大きいOCT画像を撮影する為に必要な大きさのＰＳＦのダイナミックレンジは得られなかった（例えば非特許文献１０のFigure９）。この主たる原因は、図６や図７の数値シミュレーションで示された結果から必要とされる周波数の精度が実現できていなかったためである。図８の特性は、レーザーを一定周波数で連続発振し、熱平衡状態にあるレーザーについて測定したものである。ＤＢＲの光の反射特性は、温度に依存することが知られている。このため、ＳＳＧ−ＤＢＲレーザー全体は、一箇所に取り付けられた温度計で温度を測定し、その温度が一定になるように制御されている。しかし、レーザーを走査すると、周波数毎に各領域に注入する制御電流値が変化し、その領域に生じるジュール熱が変化し、その領域の温度が局所的に変化する。この局所的な温度変化が、レーザー全体に伝わって熱平衡値に達するには数ミリ秒かかることが知られている。ＳＳ−ＯＣＴの場合、A走査の繰り返し周波数はkHz程度以上必要で、１５００の周波数を走査するためには、一周波数あたり約１００ナノ秒程度で走査しなければならない。この時間は、レーザー全体の熱平衡への緩和時間である数ミリ秒よりはるかに短い。従って、制御電流値が変化した領域の温度は、局所的には制御電流の増減に従ってレーザー全体の熱平衡値とは異なった温度になったまま、レーザーの周波数は次々と走査される。つまり、レーザーは熱平衡状態で計測された図８に基づいた制御電流値から期待される周波数とは異なる周波数で発振することになる。
【００２３】
この問題を解決するために、藤原等は、図４（Ａ）に示すＳＳＧ−ＤＢＲレーザーにおいて、周波数走査に伴って変化する制御電流注入領域毎に熱補償ヒーターを導入した（特許文献３、非特許文献１１）。各制御電流を注入する領域において、制御電流が発生するジュール熱と熱補償ヒーターが発生するジュール熱の和が一定になるように計算値を求め、制御電流の変化に合わせて同時に熱補償ヒーターの発生する熱量を調整すれば、各領域への制御電流の変化によるレーザーの温度変化を補償して抑えることができると考えた。この改良の結果、１周波数当たり５００nsの速度で走査したときの、周波数の誤差の標準偏差σfを５GHz程度に抑えることが出来た。それ以前の周波数設定比べて大幅に改良されたが、標準偏差で５GHzばらつくと、図７示すように、計測深さ3mmにおいても、PSFのダイナミックレンジは40dBよりも小さく、悪い画質のOCT画像しか得られないという問題点が残った。熱補償ヒーターを用いた方法により、個別周波数制御走査レーザーの各周波数の設定精度は大幅に改良されたけれども、画像の濃淡のダイナミックレンジが大きい高画質のＯＣＴ撮像はできなかった。
【先行技術文献】
【特許文献】
【００２４】
【特許文献１】特許第3796550号
【００２５】
【特許文献２】US 7,564,565 B2
【００２６】
【特許文献３】WO 2008/108475 A1
【非特許文献】
【００２７】
【非特許文献１】Y. Tohmori, Y. Yoshikuni, H. Ishii, F. Kano, T. Tamamura, Y. Kondo,and M.Yamamoto, 「Broad-range  wavelength-tunable  superstructure  grating  (SSG)  DBR  lasers」IEEE出版社   Journal  of  Quantum Electronics誌、 １９９３年６月発行、ｐ．１８１７〜１８２３
【００２８】
【非特許文献２】V. Jayaraman, Z. M. Chuang, and L. A. Coldren著、「Theory, design, and performance  of  extended  tuning  range  semiconductor  lasers  with  sampled  gratings」IEEE出版社  Journal  of  Quantum  Electronics誌、１９９３年６月発行、ｐ．１８２４〜１８３４
【００２９】
【非特許文献３】A. J. Ward, D. J. Robbins, G. Busico, E. Barton, L. Ponnampalam, J.P. Duck, N. D. Whitbread, P. J. Williams, D. C. J. Reid, A. C. Carter, and M.J. Wale著、「Widely  tunable  DS-DBR  laser  with  monolithically  integrated  SOA:  design  and  performance」IEEE出版社  Journal  of  Selected  Topics  in  Quantum  Electronics誌、２００５年１〜２月号発行、ｐ．１４９〜１５６
【００３０】
【非特許文献４】M. Oberg, S. Nilsson, K. Streubel, J. Wallin, L. Backbom, and T.Klinga著、「74 nm  wavelength  tuning  range  of  an  InGaAsP/InP  vertical  grating  assisted  codirectional  coupler  laser  with  rear  sampled  grating  reflector」IEEE出版社  PhotonicsTechnology Letters誌、１９９３年７月発行、ｐ．７３５〜７３７
【００３１】
【非特許文献５】R. Laroy, G. Morthier, T. Mullane, M. Todd, and R. Baets著、「Stabilisation  and  control  of  widely  tunable  MG-Y  lasers  with  integrated  photodetectors」IET 出版社  Optoelectronics誌、２００７年２月発行、ｐ．３５〜３８
【００３２】
【非特許文献６】T. Amano, H. Hiro-oka, D. H. Choi, H. Furukawa, F. Kano, M. Takeda,M. Nakanishi, K. Shimizu and K. Ohbayashi著「Optical  frequency-domain  reflectometry  with  a  rapid  wavelength-scanning  superstructure-grating  distributed  Bragg  reflector  laser」OSA出版社  AppliedOptics誌、２００５年２月１０日発行、ｐ．８０８〜８１６
【００３３】
【非特許文献７】H. Kakuma, K. Ohbayashi and Y. Arakawa著、「Optical  imaging  of  hard  and  soft  dental  tissues  using  discretely  swept  optical  frequency  domain  reflectrometry  optical  coherence  tomography  at  wavelength  from  1560  to  1600 nm」SPIE出版社  Journal  of  Biomedical  Optics誌、２００８年１〜２月号、ｐ．０１４０１２−１〜０１４０１２−５
【００３４】
【非特許文献８】H. Kakuma, D. Choi, H. Furukawa, H. Hiro-Oka, and K. Ohbayashi著、「24 mm  depth  range  discretely  swept  optical  frequency  domain  imaging  in  dentistry」、SPIE出版社  Proceedings  of  SPIE  7162誌、２００９年２月１８日発行、ｐ．７１６２０８−１〜７１６２８−８
【００３５】
【非特許文献９】G. Sarlet, G. Morthier, and R. Baets著、「Control  of  Widely  Tunable  SSG-DBR  Lasers  for  Dense  Wavelength  Division  Multiplexing」、OSA出版社  Journal  of  Lightwave  Technology誌、２０００年８月発行、ｐ．１１２８〜１１３８
【００３６】
【非特許文献１０】S. O'Connor, M. A. Bernacil, A. DeKelaita, B. Maher, D. Derickson著「100 kHz  axial  scan  rate  swept-wavelength  OCT  using  sampled  grating  distributed  Bragg  reflector  lasers」、SPIE出版社  Proceedings  of  SPIE  7168誌、２００９年２月２０日発行、ｐ．７１６８２５−１〜７１６８２５−８
【００３７】
【非特許文献１１】N. Fujiwara, H. Ishii, H. Okamoto, Y. Kawaguchi,Y. Kondo, and H. Oohashi著、「Suppression  of  Thermal  Wavelength  Drift  in  Super-Structure  Grating  Distributed  Bragg  Reflector  (SSG-DBR)  Laser  with  Thermal  Drift  Compensator」、IEEE出版社  Journal  of  Selected  Topics  in  Quantum  Electronics誌、２００７年９〜１０月号発行、ｐ．１１６４〜１１６９
【発明の概要】
【発明が解決しようとする課題】
【００３８】
個別周波数制御走査レーザーの各周波数を高精度で設定するという課題を本発明は解決する。
【課題を解決するための手段】
【００３９】
個別周波数制御走査レーザーを計測に使用時と同じ条件で走査させた状態で各制御電流を変更することを可能にする走査制御回路と、個別周波数制御走査レーザーからの出力光を導き干渉信号を発生する光干渉計と、前記干渉信号を測定する干渉信号検出器を用いて前記干渉信号を測定し、個別周波数制御走査レーザーを走査させた状態で干渉信号が所望の形になるように各制御電流を設定する装置と方法。
【発明の効果】
【００４０】
図９は、当該発明の図１に示す周波数調整装置を用いて周波数調整した個別周波数制御走査レーザーを光源として用い、図３（Ａ）に示すＳＳ−ＯＣＴ装置で撮像した人の眼の前眼部のＯＣＴ画像である。周波数調整が不十分で、ＰＳＦのダイナミックレンジが約４０ｄＢ以下の場合に撮像したところ、図９（Ａ）に示すＯＣＴ画像が得られた。前眼部は全体として明瞭に撮像できているが、虹彩や強膜の強い後方反射光に伴うノイズフロアが画面の上下に走り、濃淡のダイナミックレンジを大きく出来ず、画質が比較的悪い画像になっている。これに対し、当該発明の方法で周波数調整の精度を上げてＰＳＦのダイナミックレンジを５５ｄＢ以上にしたところ、図９（Ｂ）に示す画像が撮像された。この画像では、虹彩や強膜の強い後方反射光に伴うノイズフロアは図９（Ａ）の画像と比べて格段に弱くなり、濃淡のダイナミックレンジが大きく取れ、画像が鮮明になった。このように、周波数調整の精度を高めＰＳＦのダイナミックレンジを大きくすることにより、濃淡のダイナミックレンジが大きく画質の良いＯＣＴ撮像画像を可能にできることが当該発明のＳＳ−ＯＣＴ計測における効果である。実現された５５ｄＢのダイナミックレンジは当該発明の原理的な限界ではなく、現状では図１における走査制御回路１０２が発生する電気的ノイズで決められていて、電気回路のノイズ低減によってＰＳＦの更なる向上ができ、発明の効果を高めることができる。（なお、図９（Ｂ）画像の左右の中心付近にみられる正面反射によるノイズフロアは、角膜の表面に垂直にプローブ光が入射し、正面反射した光の強度が強すぎて、差動光検出器の測定可能範囲を信号が超えてしまい、そのフーリエ変換のノイズフロアが強く観測されていることによるものである。これを避ける方法は、プローブ光の位置を正面反射が起こらないようにずらせばよく、当該発明で解決するノイズフロアの課題とは別である。）
【００４１】
当該発明は、ＳＳ−ＯＣＴへの応用を主目的としてなされたが、高速で高分解能の分光計測などにも応用できる効果がある。当該発明によって調整された個別周波数制御走査レーザーを、例えば光吸収などの分光計測に用いれば、分光器を用いずに、計測時間数十マイクロ秒以下の高速で、波数分解能数十メガヘルツ以下の高分解能で、分光計測を行うことができる。従って、当該発明の効果は、ＳＳ−ＯＣＴに限定されず、周波数設定が高精度なされた個別周波数制御走査レーザーを用いた、光吸収や光散乱計測にも及ぶ。
【図面の簡単な説明】
【００４２】
【図１】個別周波数制御走査レーザーの周波数調整装置
【図２】（Ａ）個別周波数制御走査レーザーの周波数調整装置、（Ｂ）干渉信号の例
【図３】（Ａ）周波数走査ＯＣＴ装置の構成図、（Ｂ）周波数走査ＯＣＴ装置の光源の周波数と時間の関係
【図４】個別周波数制御走査レーザーの例：（Ａ）ＳＳＧ−ＤＢＲレーザー、（Ｂ）ＳＧ−ＤＢＲレーザー、（Ｃ）ＤＳ−ＤＢＲレーザー、（Ｄ）ＧＣＳＲレーザー、  （Ｅ）ＭＧＹレーザー
【図５】個別周波数制御走査レーザーの周波数の誤差と時間変化の例示
【図６】周波数の誤差の様々な標準偏差(σf )値における、個別周波数制御走査レーザーのＰＳＦ
【図７】ＰＳＦのダイナミックレンジ（ＤＲ）と、個別周波数制御走査レーザーの周波数の誤差の標準偏差(σf )との関係
【図８】個別周波数制御走査レーザーの制御電流と発振周波数のマッピング
【図９】人の前眼部のＳＳ−ＯＣＴ画像：（Ａ）ＰＳＦのダイナミックレンジが約４０ｄＢ以下の場合、（Ｂ）ＰＳＦのダイナミックレンジが５５ｄＢ以上の場合
【図１０】光干渉計と干渉信号の例
【図１１】所望の形の干渉信号の例
【図１２】（Ａ）個別周波数制御走査レーザーの出力周波数の間隔がサンプリン周波数間隔よりも短い場合の所望の形の干渉信号、（Ｂ）個別周波数制御走査レーザーの出力が、実質的にはアナログ信号とみなせる場合の所望の形の干渉信号
【図１３】（Ａ）光スペクトル計測装置、（Ｂ）光スペクトルの例：ＳＭＳＲが大きい場合、（Ｃ）光スペクトルの例：ＳＭＳＲが小さい場合
【図１４】個別周波数制御走査レーザーにおける周波数調整手順の例
【発明を実施するための形態】
【００４３】
［発明の原理］
個別周波数制御走査レーザーをＳＳ−ＯＣＴ計測に用いる場合と同じ条件で連続走査し、出力光を光干渉計に接続して干渉信号を観測したところ、走査開始後数分で干渉信号は一定になり安定状態を続けることを見出した。この状態で制御電流を変えて周波数調整を行うと、調整された状態を保持したまま、安定に走査を続けることも見出した。上述したように、この条件では各個別の周波数において静的な熱平衡は達成されていない。ＳＳＧ−ＤＢＲ領域での発熱の変化がレーザー全体に達し熱平衡になる数ミリ秒に比べると、走査中での各周波数状態での滞在時間は数百ナノ秒以下で、熱平衡に向かって変化する前の断熱的状態の間に、次々と周波数が走査される状態である。しかし、各発振周波数は断熱状態で次の周波数に走査されるが、レーザー全体としては、同じ走査状態が繰り返されているので、熱平衡に達していて安定である。従って、この条件下で繰り返される各周波数は再現性良く安定に動作する。この状態で制御電流を調整して所望の周波数精度が達成できれば、ＳＳ−ＯＣＴ計測時と同じ走査条件で高精度の周波数調整が実現できることになり、高画質のＯＣＴ画像が取得できる。走査したままの個別周波数制御走査レーザーを光干渉計に接続し、干渉信号検出器で干渉信号を計測し、個別周波数制御走査レーザーの各周波数が正確に設定された場合に期待される所望の形の干渉信号が得られるように制御電流を調整すれば、正確な周波数調整が実現できる。ＳＳ−ＯＣＴ計測に用いるときと同じ走査条件の状態で正確な周波数調整ができるので、このように調整された個別周波数制御走査レーザーをＳＳ−ＯＣＴに用いれば、濃淡のダイナミックレンジが大きい高画質のＯＣＴ画像が撮像できる。これが当該発明の基本的な原理である。
【実施例１】
【００４４】
図１に第一の実施の形態を示す。個別周波数制御走査レーザー１０１は、走査制御回路１０２で制御されている。走査制御回路１０２は、個別周波数制御走査レーザー１０１を繰り返し走査させた状態で制御電流値を変更し、その変更された制御電流値に従って個別周波数制御走査レーザー１０１を継続して繰り返し走査させることができる機構を持っている。個別周波数制御走査レーザー１０１の出力光は光干渉計１０５に接続されている。光干渉計１０５の出力は、個別周波数制御走査レーザー１０１の出力する光の周波数の変化に従って時間変化し、光の干渉信号を出力する。光の干渉信号は、干渉信号検出器１０６により電気信号に変換される。干渉信号検出器１０６には、干渉信号を表示する機能を持たせることができ、その表示された干渉信号と所望の形の干渉信号との差を人が測定し、差が無くなるように手動で走査制御回路１０２を設定することもできる。あるいは、干渉信号検出器１０６に所望の形の干渉信号を参照値として入力しておいて、それと、光干渉計１０５から取得された干渉信号との差を自動で検出し、差が無くなるように計算された制御電流値を、信号経路１０３を介し走査制御回路１０２に転送するようにすることもできる。
【００４５】
光干渉計１０５が出力する「所望の形の干渉信号」について説明する。「所望の形の干渉信号」とは、数７式で表される個別周波数制御走査レーザーが実際に走査する各周波数Fiの誤差Δiが十分小さい条件で精度よく設定できたため、数６式に従って正確に周波数設定された個別周波数制御走査レーザーが走査されたときに期待できる干渉信号の形に近似的に一致する形の干渉信号を意味する。この時、誤差Δiが小さければ小さいほど、より望ましい「所望の形の干渉信号」になる。
【００４６】
図１０（Ａ）にファブリ・ペロー干渉計を示す。この干渉計は、間隔d離れて平行に置かれたミラーA８０１とミラーB８０２で構成されている。この平行ミラーに強度一定の入力光を入射すると、光の干渉によって、光の透過率は図１０（Ｄ）に示すように、周波数の関数として周期的に鋭い透過率のピーク示す。透過率のピークの周波数間隔は、fFP=c/2dで与えられ、周波数によらず一定である。個別周波数制御走査レーザーの周波数間隔がδf=fFPになるように間隔dを設定すれば、レーザーの各設定周波数fi=f0+iδfで干渉信号光は極大値を示すようにできる。δfとfFPの関係はδf=mfFP（m：１以上の整数）であれば、干渉信号を極大値に合わせることによって、所望の形の干渉信号が得られる。
【００４７】
図１０（Ｂ）はマイケルソン干渉計を示す。入力光をハーフミラー８０５で反射光と透過光に分割する。反射光はミラーC８０３で反射されて、光路C８０６を往復する。透過光はミラーD８０４で反射されて、光路D８０７を往復する。それぞれの光路を往復した光は、ハーフミラーで干渉し、光干渉信号を生成する。光干渉信号の変動部分は、入力光の周波数fに依存して次式で与えられる正弦波形に比例した変化を示す。
【数１０】
2015001452.tif 000012
この関係を図１０（Ｅ）に示す。Δlは、光路C８０６と光路D８０７の往復の光路長の差である。数１０式は、光周波数fMZ=c/Δlの周期で変動する。図１０（Ｃ）は光ファイバで接続されたマッハ・ツェンダー干渉計を示す。入力光は分波器８１０で光路E８１３と光路F８１４に分波される。光路F８１４との光路差がΔlになるように、光路E８１３には光ディレイ８１１が挿入されている。光路E８１３と光路F８１４は、分波器８１０から始まり合波器８１２で終わる。これら２つの光路を通った光は合波器８１２に導かれ干渉する。合波器８１２からは、位相がπ異なる２つの干渉信号光＋と干渉信号光−が出力される。これら干渉信号の変動部分は、数１０式に比例し、位相が互いにπだけ異なる。干渉信号光＋と干渉信号光−を差動検出すると、数１０式に比例する信号を検出できる。
【００４８】
干渉信号が数１０式で与えられるときの所望の形の干渉信号の例を図１１に示す。図１１（Ａ）は、光周波数が連続的に変化した場合、数１０式で与えられる信号は正弦波形で変化し、その周期がfMZであることを示している。走査させた個別周波数制御走査レーザーからの出力光を図１０（Ｂ）や図１０（Ｃ）に示す干渉計の入力光として用いると、出力される光の周波数は離散的であるので、図１１（Ａ）に破線で示すような連続的変化は示さない。図１１（Ｂ）の太い実線は、個別周波数制御走査レーザーの出力する周波数間隔δfが周期fMZに等しく（δf = fMZ）、かつ周波数f0がf0Δl=(m−1/4)c（mは正の整数）になるようにレーザーの各周波数を調整した場合の干渉信号である。この実線は、値が０なので、横軸の時間軸に沿っている。これが、この場合の所望の形の干渉信号である。正弦波が零レベルになる光周波数は、正弦波の変化が激しいところなので、零レベルからのずれの検出の感度が良い。周波数f0は、f0Δl=(m+1/4)c（mは正の整数）に設定してもよい。
【００４９】
図１１（Ｃ）は、δf = fMZであるが、周波数f0がf0Δl=mc（mは正の整数）になるようにレーザーの各周波数を調整した場合の所望の形の干渉信号で、干渉信号が極大値をとる周波数にレーザーの各出力周波数を合わせた場合である。同様に、周波数f0をf0Δl=(m+1/2)c（mは正の整数）に調整し、所望の形の干渉信号が極大値の替りに極小値をとるように調整してもよい。極大や極小を用いた調整は、零レベルに合わせて調整する場合よりも、設定周波数の検出感度が悪い。
【００５０】
図１１（Ｄ）は、fMZ=2δfで、かつ周波数f0がf0Δl=(m−1/4)c（mは正の整数）になるように調整した場合の所望の形の干渉信号で、レーザーのどの周波数でも干渉信号が零になる。図１１（Ｅ）もfMZ=2δfで、かつ周波数f0がf0Δl=mc（mは正の整数）になるようにレーザーの周波数を調整した場合の所望の形の干渉信号で、干渉信号は交互に極大と極小をとる。所望の形の干渉信号は、矩形波になる。図１１（Ｆ）に示す所望の形の干渉信号は、fMZ=4δfの場合で、順次、極大→零→極小→零→極大と繰り返す。
【００５１】
個別周波数制御走査レーザーの周波数間隔δfと、干渉計の周波数周期fMZが整数の比である必要はない。期待される所望の形の干渉信号は、周波数周期fMZが決まっていれば、設定したい周波数間隔δfに対し計算で求めることができる。その計算値と観測される干渉信号が一致するように、レーザーの各周波数を調整すればよい。調整した結果の所望の形の干渉信号を太線で示すと、図１１（Ｇ）に示すような図が得られる。
【００５２】
個別周波数制御走査レーザーの全出力周波数数をNとして、fMZ＞２Nδfに設定した周波数周期fMZの光干渉計を用いると、走査周波数の全領域で単調増加または単調減少する一価関数の干渉信号を所望の形の干渉信号にすることができ、自動調整が容易になる。また、図１１に示す干渉信号は、δfがfMZよりも小さい場合の例が示してあるが、mを２以上の正の整数として、fMZ=δf/mの条件を満たす光干渉計も利用できる。この場合、干渉信号が零レベルを通過するときの微分変化が、fMZ=δfのときより大きいので、零レベルの検出感度がより良い周波数調整が可能になる。この時の検出感度を計算する。数１０式は、cos(2πf/fMZ)と書ける。周波数fがΔf変動した時の信号の変化量は、この式の微分にΔfを乗じたものになり、fMZ=δf/mの条件を代入すると、−sin(2πf/fMZ)(2π/fMZ)Δf=−(2πm/δf)Δfが得られる。（ここで、周波数が調整できた時の条件sin(2πf/ fMZ)=1を用いた。）通常の電子回路はフルスケールの０．１％程度の精度で電圧を測定することは可能であり、干渉信号の振幅を１に規格化すれば10-4の信号の検出が、干渉信号検出器を用いて可能である。(2πm/δf)Δf =10-4とおいてΔfを求めると、周波数測定の感度の程度として、次式が得られる。
【数１１】
2015001452.tif 000013
例として、数値シミュレーションで用いた周波数の周期δf = 6.25GHzの場合、m=1でΔf=99kHz、m=2でΔf=50kHz、m=3でΔf=33kHzになる。このように、光干渉計を用いた周波数の検出感度は極めて高く、図６（Ｄ）、（Ｅ）、（Ｆ）で例示した標準偏差σf=0.01GHz=10MHzの精度の周波数調整の実現は、原理的に可能である。
【００５３】
これまでは、個別周波数制御走査レーザーが出力する周波数間隔とデータをサンプリングする周波数間隔δfは同じ場合について説明した。制御電流をステップ状（階段状）に変化（スイッチ）すると、電気回路やレーザー発振の応答の速さが有限であるために、レーザー発振の周波数は完全にはステップ状に変化せず、遷移現象を伴って変化する。この遷移の時間内では周波数が誤差を持つことになる。遷移現象の振幅は、制御電流のステップ状の変化が小さいほど小さい。そこで、複数回に分けてスイッチすれば、各スイッチでの遷移現象の振幅は小さくできる。図１２（Ａ）は、３回に分けてスイッチし、３つステップ毎に一つのデータをサンプリングする例を示す。太い実線がサンプルする周波数のデータであり、点線で示したデータは、サンプルしないデータである。図１に示す方法を用いた周波数調整では、図１２（Ａ）において点線で表したデータも含めて行い、所望の形の干渉信号は実線のデータと点線のデータを含めたものである。制御電流のスイッチ間隔を小さくすればするほど、遷移現象の振幅は小さくなる。スイッチ間隔の時間を、電気回路やレーザー発振の応答時間と同等またはそれ以下にすると、レーザーの出力は実質的には図１２（Ｂ）に実線で示すようにアナログ出力になる。電気回路やレーザー発振の応答の時間だけ、点線で示す制御電流の変化よりも遅れて出力される。このような場合でも、個別周波数制御走査レーザーの制御電流は各ステップ変化毎に設定しなければならず、所望の形の干渉信号が図１２（Ｂ）に示す実線の様に実質的にはアナログ信号になるという条件で、図１に示した周波数調整法が適用できる。
【００５４】
図１に示す個別周波数制御走査レーザー１０１を制御する「走査制御回路」１０２について説明する。個別周波数制御走査レーザー１０１は、必要な制御電流が供給されると、それに従って発振する。走査制御回路１０２が一定の制御電流を定常的に供給すれば、個別周波数制御走査レーザー１０１は定常的に連続発振（ＣＷ発振）する。走査制御回路１０２が供給する制御電流の値を変えると、個別周波数制御走査レーザー１０１はそれに従って発振状態を変える。短時間の間隔で次々と供給電流を変化させれば、個別周波数制御走査レーザー１０１は走査されることになる。このような走査は、通常はあらかじめ走査制御回路１０２内のメモリに書き込まれた制御電流のテーブル（一連のデータ）に従って行われ、走査の仕方を変える場合は、走査を止めて制御電流のテーブルを書き換えてから走査を開始するという手順が、従来技術では行われていた。当該発明では、走査制御回路１０２は、個別周波数制御走査レーザー１０１を走査させた状態で制御電流を変更することができる機構を持っているものとするが、その機構は通常よく知られている電子回路技術を用いて容易に実現できる。方法の一例について説明する。走査制御回路は、メモリの中に全ての制御電流の値を記憶していて、順次メモリのアドレスを走査の速度で変化させ、各アドレスで出力されるデジタルの制御電流の値をデジタル・アナログ変換して個別周波数制御走査レーザーの各制御電流を供給する。このとき、例えば制御電流値のメモリとして、走査メモリと入力メモリの２つを持ち、走査メモリは走査する目的に用い、入力メモリは制御電流変更値を入力する目的に用い、図３（Ｂ）において、最終周波数fNまで走査メモリに従って走査した直後に最初の周波数f1を出力する前に、走査中に変更された入力メモリの内容を走査メモリに短時間で転送して走査を続ければ、走査を乱すことなく、制御電流の値を変更することができる。また別の方法としては、走査メモリの読み出しと書き込みの機能を持たせ、制御電流変更値１０３を直接走査メモリに必要な時点で短時間内に書き込み、走査を続けるという機構も可能である。この場合、書き込みの瞬間だけ干渉信号が乱れるが、周波数調整時であれば支障はない。
【００５５】
図４に示す個別周波数制御走査レーザーの出力周波数は、３つの制御電流によって決定される。図１に示す周波数調整方法において、これら３つの制御電流を同時に変化させて調整すると手順が複雑になる。変化させるパラメータ（制御電流値）は一つであることが望ましい。位相領域に注入する制御電流だけを単一パラメータとして図１の周波数調整方法を可能にする手順について説明する。個別周波数制御走査レーザーが定常状態で連続発振する状態でのレーザー出力の光スペクトル測定を、制御電流を逐次変えながら繰り返し、制御電流の関数として一連のスペクトルデータを計測する装置を図１３（Ａ）に示す。走査制御回路１０２にコンピュータ３００から各制御電流設定命令を送り、個別周波数制御走査レーザー１０１を、各制御電流が一定の条件で連続発振（CW発振）させる。この出力光のスペクトルを、コンピュータからの光スペクトル取得命令に従って、光周波数の関数として光スぺクトラムアナライザ３０２で計測する。計測された光スペクトルデータをコンピュータに転送し、記憶する。この操作を、プログラムによって自動的に制御電流を変えながら繰り返し、一セットのデータをコンピュータのメモリ内に取得する。
【００５６】
具体例として、図４（Ａ）に示すＳＳＧ−ＤＢＲレーザーの場合について説明する。ゲイン電流とＳＯＡ電流は100mAの一定値に設定した。位相電流も、走査時の平均値に近い一定値6mAに設定した。後方SSG-DBR電流と前方SSG-DBR電流を、0.2mAの間隔で、0.2mAから50mAまで網羅的に変化させて、各制御電流での連続発振のスペクトルを計測した。計測されたスペクトルの例を図１３（Ｂ）と図１３（Ｃ）に示す。スペクトルが単一モードであるためには、最強のスペクトルの強度が、他のサイドモードよりも圧倒的な強い必要がある。最強のスペクトル強度と次に強いスペクトル強度との比であるＳＭＳＲが大きいほど、スペクトルの純度が高い。図１３（Ｂ）のＳＭＳＲは49.6dB 、図１３（Ｃ）のＳＭＳＲは28.7dBである。図１３（Ｃ）は、スペクトル純度が悪い場合である。レーザー発振の純度が比較的悪いＳＭＳＲが35dB以下のデータを除外し、スペクトルのピークの光周波数を横軸にとり、縦軸に前方ＳＳＧ−ＤＢＲ電流をプロットした図を図１４（Ａ）に示す。丸印で示したデータ点は、７本の特性曲線にそって狭い範囲に分布する。これらの特性曲線に沿ったデータ分布は狭く、最少二乗法によって各特性曲線の最適値を決定することを可能にする。７つのそれぞれの曲線を、３次関数を用いて最少二乗法で最適化した結果を実線で示す。後方ＳＳＧ−ＤＢＲについても同様に特性曲線の最適値が得られる。従来技術として図８に示した前方ＳＳＧ−ＤＢＲ電流と後方ＳＳＧ−ＤＢＲ電流の２次元上に発振周波数をマッピングする方法に比べると、この当該発明者が見出した図１４（Ａ）に示したプロットから最適特性曲線を得る手順に従えば、遥かに簡単なデータ処理で、各光周波数を発振するための前方ＳＳＧ−ＤＢＲ電流と後方ＳＳＧ−ＤＢＲ電流の最適値が得られる。
【００５７】
上記のデータ処理で得られた最適の特性曲線を、前方SSG-DBR電流については細い実線で、後方ＳＳＧ−ＤＢＲ電流については細い破線で、図１４（Ｂ）に示す。所望の光周波数の値で縦軸に平行に引いた直線と各最適特性曲線との交点から、所望の光周波数を出力するための前方ＳＳＧ−ＤＢＲ電流と後方ＳＳＧ−ＤＢＲ電流の対の値が得られる。一つの前方（または後方）ＳＳＧ−ＤＢＲ電流に対して２つの後方（または前方）ＳＳＧ−ＤＢＲ電流の値が可能の場合があるが、どちらを選択してもよい。電流の変化量が同じなら、前方ＳＳＧ−ＤＢＲ電流（実線）の方が、後方ＳＳＧ−ＤＢＲ電流（点線）よりも幅広い周波数領域をカバーするようになっている、つまり制御電流で設定する周波数の目盛の間隔が前方と後方のＳＳＧ−ＤＢＲで異なっている。これが前述のバーニア効果である。このバーニア効果により、前方ＳＳＧ−ＤＢＲ電流と後方ＳＳＧ−ＤＢＲ電流の一つの対により、幅広い光周波数の領域で一義的に出力周波数が決められるようになっている。図１４（Ｂ）において、太線は周波数範囲187.2815THzから192.2041THzまで6.25GHzの周波数間隔で走査するために最適化曲線にそって選択した各ＳＳＧ−ＤＢＲ電流の一例である。実線が前方ＳＳＧ−ＤＢＲ電流、破線が後方ＳＳＧ−ＤＢＲ電流である。このように各ＳＳＧ−ＤＢＲ電流を設定すると、記号ｓで示した光周波数の値で、片方または両方のＳＳＧ−ＤＢＲ電流を大きな値から小さな値に急激にスイッチしなければならない。この時、電流の不連続変化が大きいことと、レーザーの発振モードが飛ぶために、レーザーの発振が不安定になり、干渉信号の乱れが発生する。この乱れは従来技術では除去不可能とされていた。従来技術は非特許文献１０のFig.7や非特許文献１１のFig.７とFig.８で、この乱れ（ノイズ）が議論されている。この問題を解決するために、当該発明者は、電流をスイッチした後の電流設定値のままで繰り返し走査を続けると、やがて出力周波数は安定することを確認した。安定した後に次の周波数への走査を開始し、全体のデータをサンプリングしたあとで安定のために余分に挿入したデータをデータ解析からは取り除けば、急激な電流値のスイッチによるノイズが含まれたデータを取り除けることを見出した。
【００５８】
図１４（Ｃ）は、図１４（Ｂ）に示す太線で示した制御電流値に従ってレーザーを走査し、fMZ=351GHzのマッハ・ツェンダー干渉計を通して測定した信号である。記号ｓで記した制御電流の飛びの後は、同じ制御電流値を５回繰り返して干渉信号が安定するまで待つようにした。このため、横軸は光周波数ではなく、出力番号である。一見すると、図１４（Ｃ）の干渉信号には沢山の飛びがあり、所望の正弦波形の干渉信号から大きくずれているように見えるが、位相電流の調整によって、所望の正弦波形に容易に合わせることができる。図１４（Ｄ）の細線は、図１４（Ｃ）において破線で囲まれたデータの部分を、横軸を拡大して記したものである。太線は、fMZ=351GHzのマッハ・ツェンダー干渉計から期待させる所望の形の干渉波形である。観測された細線が周波数の低い方に太線からずれている場合は位相電流を一致するまで増加し、反対に周波数の高い方にずれている場合は位相電流を一致するまで減少することによって、所望の正弦波形に容易に合わせることができる。位相電流の初期値は、どの出力番号でも6mA に設定したが、調整の結果は図１４（Ｅ）に示すように6mAの上下を規則的に変化する値になる。このように、図１で示した方法による周波数の調整は、単一パラメータである位相電流を調整することのみで達成できる。
【００５９】
位相電流を調整した後の干渉信号を図１４（Ｆ）に示す。記号ｓで示す制御電流の急激な飛びの後に、乱れが見られる。これらの出力番号値では、制御電流の急激な飛びによって、所望の干渉信号の値に制御できない。上述のように、所望の干渉信号の値が得られるまで同じ制御電流値を繰り返した。この図の場合は、飛びの後に５回同じ制御電流値を繰り返した。従って、各制御電流の飛びの後の４つのデータは、最終データとしては不要で、解析から取除いた。その結果、図１４（Ｇ）に示すようにノイズのない所望の形の干渉信号が得られた。横軸は、余分なデータを取り除いた後の選択されたデータの番号である。
【００６０】
なお、各制御電流を増減すると、その領域の光吸収の増減が起き、レーザーの出力強度が変わる。出力強度が変動すると干渉信号の振幅が変わるので、振幅で合わせる操作に誤差が生じる。従って、図１４（Ｃ）の測定の前には、レーザーの出力がどの出力番号でも一定なるようにＳＯＡ電流を調整する必要がある。また、位相電流のわずかな調整でも、レーザーの出力強度はわずかに変わるので、図１４（Ｇ）に示す最終データ取得前にも、レーザー強度一定化の調整をすべきである。レーザー強度一定化の調整は、容易に自動化で行える。
【００６１】
図１３と図１４に示した手順に従って、数多くのＳＳＧ−ＤＢＲレーザーを図１に示す調整装置を用いて調整したところ、わずかな例ではあるが、前方ＳＳＧ−ＤＢＲ電流や後方ＳＳＧ−ＤＢＲ電流を、図１４（Ａ）と図１４（Ｂ）に示した手順で求めた値からわずか変化させた方がより安定な出力が得られるという場合もあった。位相電流のみを用いた単一パラメータを用いた調整が基本ではあるが、当該発明はそれに限定されるものではなく、前方ＳＳＧ−ＤＢＲ電流および後方ＳＳＧ−ＤＢＲ電流および位相電流のうちの、２つまたは３つを用いて調整する場合も含む。
【００６２】
図１４（Ｃ）から（Ｇ）で説明した例では、分かり易いように、所望の正弦波の波形が図１１（Ｇ）に示すように正弦波である場合を示した。周波数調整のための位相電流の変化は、図１４（Ｅ）に示すように規則的であるので、所望の干渉信号の変化が正弦波形の変化が顕在化する図１１（Ｇ）の場合だけでなく、図１１（Ｂ）から（Ｆ）に示すような場合でも、図１４で説明した周波数調整手順は可能である。また、図１において、光干渉計１０５を順次変えて、図１１に示した複数の種類の所望の干渉信号が観測できるようにして、調整を繰り返し、周波数調整の精度を向上してもよい。
【実施例２】
【００６３】
正弦波が零になるところでは、正弦波の変数による微分は最大になり、変数の変化に対してもっとも敏感である。従って、図１１に示す干渉信号の中では、周波数の所望の設定値からのずれの検出感度が最も高いのは図１１（Ｂ）の場合である。しかし、この場合、全ての周波数で零に合わせるので、周波数間隔δfの整数倍異なった周波数に誤って調整しても、信号の形としての乱れは観測されない。このとき、図１１（Ｇ）に示すようなゆっくり変化する信号では乱れが生じるので、その乱れから調整の誤りが判断できる。このように、異なる周波数間隔の光干渉計を用いれば、周波数調整の誤りを複数の光干渉計の信号の比較から検出でき、周波数調整の誤りをなくすことができる。
【００６４】
図２（Ａ）は、周波数間隔が互いに異なる３つの光干渉計を同時に用いる周波数調整器の例である。個別周波数制御走査レーザー１０１の出力光を、光分波器１１０で４分割し、互いに周波数間隔の異なる３つの光干渉計と強度測定のための強度信号検出器１１４に導く。光干渉計Ａ１１１、光干渉計Ｂ１１２、光干渉計Ｃ１１３の出力信号は、それぞれ、光信号検出器Ａ１１５、光信号検出器Ｂ１１６、光信号検出器Ｃ１１７で検出される。各検出器で検出された信号は、そのままオシロスコープなどで表示してもよいし、コンピュータにデータを取り込みながら、コンピュータのディスプレイに表示してもよい。
【００６５】
図１４に示す調整手順で説明したように、周波数調整の前に、レーザーの出力を出力番号によらずに一定することが望ましい。走査制御回路１０２は、入力された制御電流値に従って個別周波数制御走査レーザーを走査させる機能を持っていて、ＳＯＡ電流についても適用できる。強度信号検出器１１４で検出する強度が一定になるように、信号経路１０３を介して、走査制御回路１０２でＳＯＡ電流を設定し、個別周波数制御走査レーザー１０１の出力強度を一定にする。
【００６６】
個別周波数制御走査レーザー１０１の設定したい周波数間隔をδfとする。光干渉計Ａ１１１の周波数周期fMZAをδfに等しくすると、調整されたとき所望の形の干渉信号は図２（Ｂ）の干渉信号Ａようにいつも零の信号である。光干渉計Ｂ１１２の周波数周期fMZBを、fMZB＝４δfになるようにすると、図１１（Ｆ）に示すような信号が所望の形の干渉信号になり、一部の不完全部分を除いて、図２（Ｂ）に示す干渉信号Ｂになる。光干渉計Ｃ１１３の周波数周期fMZCを、fMZC＝６０δfになるようにすると、ゆっくり変動する正弦波形が所望の形の干渉信号になり、一部の不完全部分を除いて、図２（Ｂ）に示す干渉信号Ｃになる。図２（Ｂ）では、４番目の出力番号のデータが誤って３番目の出力番号の周波数と同じに設定された場合が示されており、干渉信号Ａでは誤りは分からず、干渉信号Ｃでも信号の変化が小さいところで誤りが不明瞭であるが、干渉信号Ｂを見れば、−１となるべきデータが０になっており、誤りがすぐに分かる。出力番号１６のデータは、誤って出力番号１２と同じ周波数に設定されている。この場合、干渉信号Ａにも干渉信号Ｂにも異常は見られないが、干渉信号Ｃでは、はっきりと異常が検出できる。このように、一つの光干渉計だけでは検出できない設定の誤りも、複数の周波数周期の異なる光干渉計を同時または順次に用いることにより検出できる。図２に示す光干渉計の数や組み合わせ一例であって、他にも様々な組み合わせが可能である。
【００６７】
実施例２で説明した同時に複数の光干渉計を用いる方法は、実施例１の方法の一つの例であり、図１に示す方法には、複数の光干渉計を用いる場合も含まれる。
【符号の説明】
【００６８】
１０１  個別周波数制御走査レーザー
１０２  走査制御回路
１０３  制御電流値の信号経路
１０５  光干渉計
１０６  干渉信号検出器
１１０  光分波器
１１１  光干渉計Ａ
１１２  光干渉計Ｂ
１１３  光干渉計Ｃ
１１４  強度信号検出器
１１５  干渉信号検出器Ａ
１１６  干渉信号検出器Ｂ
１１７  干渉信号検出器Ｃ
１１８  コンピュータ
２００  後方ＳＳＧ−ＤＢＲ領域
２０１  位相領域
２０２  ゲイン領域
２０３  前方ＳＳＧ−ＤＢＲ領域
２０４  ＳＯＡ領域
２１０  後方ＳＧ−ＤＢＲ領域
２１１  位相領域
２１２  ゲイン領域
２１３  前方ＳＧ−ＤＢＲ領域
２１４  ＳＯＡ領域
２２０  後方ＤＢＲ領域
２２１  位相領域
２２２  ゲイン領域
２２３  前方ＤＢＲ領域
２２４  ＳＯＡ領域
２３０  後方Ｓ−ＤＢＲ領域
２３１  位相領域
２３２  回折格子カプラ領域
２３３  ゲイン領域
２３４  ＳＯＡ領域
２４０  ＤＢＲ−Ａ領域
２４１  ＤＢＲ−Ｂ領域
２４２  マルチモード干渉計カプラ領域
２４３  位相領域
２４４  ゲイン領域
２４５  ＳＯＡ領域
３００  コンピュータ
３０２  光スペクトラムアナライザ
５０１  周波数走査光源
５０２  ＯＣＴ干渉計
５０３  光ファイバ
５０５  分波カプラ
５０６  サーキュレータ
５０７  コリメータ
５０８  ガルバノミラー
５０９  レンズ
５１０  試料
５１１  合波カプラ
５１２  サーキュレータ
５１３  コリメータ
５１４  レンズ
５１５  参照ミラー
５１６  差動光検出器
５１７  A/D変換器
５１８  コンピュータ
５１９  OCT断層画像
５２０  表示器
５３０  試料光路
５３１  参照光路
５４１  サンプル開始トリガー
５４２  サンプリングクロック
８０１  ミラーＡ
８０２  ミラーＢ
８０３  ミラーＣ
８０４  ミラーＤ
８０５  ハーフミラー
８０６  光路Ｃ
８０７  光路Ｄ
８１０  分波器
８１１  光ディレイ
８１２  合波器
８１３  光路Ｅ
８１４  光路Ｆ
【図１】
2015001452.tif 000014
【図２】
2015001452.tif 000015
【図３】
2015001452.tif 000016
【図４】
2015001452.tif 000017
【図５】
2015001452.tif 000018
【図６】
2015001452.tif 000019
【図７】
2015001452.tif 000020
【図８】
2015001452.tif 000021
【図９】
2015001452.tif 000022
【図１０】
2015001452.tif 000023
【図１１】
2015001452.tif 000024
【図１２】
2015001452.tif 000025
【図１３】
2015001452.tif 000026
【図１４】
2015001452.tif 000027