JP6259370B2 - Optical coherence tomography device - Google Patents

Optical coherence tomography device Download PDF

Info

Publication number
JP6259370B2
JP6259370B2 JP2014151157A JP2014151157A JP6259370B2 JP 6259370 B2 JP6259370 B2 JP 6259370B2 JP 2014151157 A JP2014151157 A JP 2014151157A JP 2014151157 A JP2014151157 A JP 2014151157A JP 6259370 B2 JP6259370 B2 JP 6259370B2
Authority
JP
Japan
Prior art keywords
light
signal
time
light source
phase
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
JP2014151157A
Other languages
Japanese (ja)
Other versions
JP2016024176A (en
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.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
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 Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2014151157A priority Critical patent/JP6259370B2/en
Publication of JP2016024176A publication Critical patent/JP2016024176A/en
Application granted granted Critical
Publication of JP6259370B2 publication Critical patent/JP6259370B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Description

本発明は、光干渉断層装置に関し、詳細には、主に医療で使用されている3次元断層画像を取得するための光干渉断層装置に関する。   The present invention relates to an optical coherence tomography apparatus, and more particularly to an optical coherence tomography apparatus for acquiring a three-dimensional tomographic image mainly used in medicine.

時間に対して光の波数が連続的に変化(掃引)する光源を波長掃引光源と呼ぶが、その様な光源において、時間に対して光の波数が線形的に(1次関数的に)掃引する光源を波数リニアな波長掃引光源と呼ぶ。   A light source in which the wave number of light continuously changes (sweep) with respect to time is called a wavelength-swept light source. In such a light source, the wave number of light sweeps linearly (linearly) with respect to time. The light source to be called is called a wave number linear wavelength swept light source.

光干渉断層装置としては、波数リニアではない波長可変光源を用いた場合に得た干渉波形信号を、波数リニアとなる波長可変光源を用いた場合と同等の干渉波形信号に変換(リサンプリング)して3次元断層画像を取得するSwept Source Optical Coherence Tomography(SS−OCT)システムが知られている(非特許文献1)。非特許文献1には図1に示す構成のSS−OCTシステムが記載されている。図1に示すSS−OCTシステムは、櫛歯型のスペクトルとなるように波長掃引を行なう光源101と、カップラ102と、干渉計103と、フォトディテクタ104と、ピーク時刻取得部1106と、リサンプリング部107と、フーリエ変換部108とを備えている。さらに干渉計103は、バランスフォトディテクタ1031と、サーキュレータ1032と、カップラ1033と、ファイバコリメータ1038と、ミラー1035と、ファイバコリメータ1039と、光偏向器1037とを有している。図中の各構成は光ファイバ(二重線)または電気ケーブル(一本線)で接続されている。カップラ1033からファイバコリメータ1038を通ってミラー1035までの光路を参照アーム1034と呼び、また、カップラ1033からファイバコリメータ1039と光偏向器1037を通って観察対象Oの中の反射面Bまでの光路をサンプルアーム1036呼ぶこととする。参照アーム1034にあるファイバコリメータ1038からミラー1034の間と、サンプルアーム1036にあるファイバコリメータ1039から光偏向器1037を経由した反射面Bまでの光路は空間系となっており、空間中(通常は空気中)を光が伝播する。   As an optical coherence tomography device, the interference waveform signal obtained when using a wavelength variable light source that is not wave number linear is converted (resampled) into an interference waveform signal equivalent to that when a wavelength variable light source that is wave number linear is used. There is known a Swept Source Optical Coherence Tomography (SS-OCT) system that acquires a three-dimensional tomographic image (Non-Patent Document 1). Non-Patent Document 1 describes an SS-OCT system configured as shown in FIG. The SS-OCT system shown in FIG. 1 includes a light source 101 that performs wavelength sweeping so as to obtain a comb-like spectrum, a coupler 102, an interferometer 103, a photodetector 104, a peak time acquisition unit 1106, and a resampling unit. 107 and a Fourier transform unit 108. Further, the interferometer 103 includes a balance photodetector 1031, a circulator 1032, a coupler 1033, a fiber collimator 1038, a mirror 1035, a fiber collimator 1039, and an optical deflector 1037. Each component in the figure is connected by an optical fiber (double line) or an electric cable (single line). The optical path from the coupler 1033 through the fiber collimator 1038 to the mirror 1035 is referred to as a reference arm 1034, and the optical path from the coupler 1033 through the fiber collimator 1039 and the optical deflector 1037 to the reflecting surface B in the observation object O. The sample arm 1036 will be called. The optical path from the fiber collimator 1038 to the mirror 1034 in the reference arm 1034 and from the fiber collimator 1039 in the sample arm 1036 to the reflecting surface B via the optical deflector 1037 is a spatial system, and is in space (usually Light propagates in the air.

光源101は上述の通り櫛歯型のスペクトルとなるように波長掃引を行う光源であるが、これを実現するため、レーザ共振器内に、所定の自由スペクトル間隔(free spectral range:FSR)を持つフィルタ(櫛歯型スペクトルの透過率を持つフィルタ:櫛歯フィルタ)を入れている。このようなフィルタをレーザ共振器内に挿入したレーザは、フィルタの透過率に応じて出力光強度が変化する。たとえば、フィルタの透過率が極大/極小となる光の周波数となった場合は、出力光強度は極大/極小となる。このため、レーザ共振器の光の共振周波数を連続的に時間変化させた場合に、共振する光の周波数に応じて出力光強度が変わるが、この様子を分散分光方式(モノクロメータを使用する方式)の光スペクトルアナライザで観測すると、光源101のスペクトルは、FSRの間隔でピークを持つ櫛歯型となる。   The light source 101 is a light source that performs wavelength sweeping so as to have a comb-shaped spectrum as described above. To realize this, the laser resonator has a predetermined free spectral range (FSR) in the laser resonator. A filter (a filter having a comb-shaped spectrum transmittance: a comb-tooth filter) is included. In a laser in which such a filter is inserted in the laser resonator, the output light intensity changes according to the transmittance of the filter. For example, when the filter has a light frequency at which the transmittance is maximum / minimum, the output light intensity is maximum / minimum. For this reason, when the resonance frequency of the laser resonator light is continuously changed over time, the output light intensity changes according to the frequency of the resonating light. ), The spectrum of the light source 101 has a comb shape having peaks at the intervals of FSR.

SS−OCTでは、光源のコヒーレンス長が、SS−OCTで得られる断層画像で観測できる奥行(深達長)を決め、コヒーレンス長が長ければ深達長が長くなる。このため、コヒーレンス長の長い光源であることがSS−OCTにおいて重要である。   In SS-OCT, the coherence length of the light source determines the depth (depth depth) that can be observed in the tomographic image obtained by SS-OCT, and if the coherence length is long, the depth depth becomes long. For this reason, it is important in SS-OCT that the light source has a long coherence length.

上記に示す櫛歯型フィルタを共振器に挿入したレーザは、フィルタの透過率が極大、つまり、レーザ出力が極大となる波長の出力光のコヒーレンス長が長くなることが知られている。したがって、光源光強度の強い光、つまり、コヒーレンス長が長い光で得られた干渉計出力のみを抽出して信号処理すると、深達長の長い観測が可能となる。以下に説明する方法は上記抽出を行う方法であり、光源光101の出力光の強度変調のピーク値で干渉計103の出力をサンプリングする方法である。   It is known that the laser in which the comb filter shown above is inserted into a resonator has a maximum filter transmittance, that is, the coherence length of output light having a wavelength at which the laser output is maximized. Therefore, if only the output of the interferometer obtained by the light having a strong light source light intensity, that is, the light having a long coherence length is extracted and signal-processed, observation with a long depth can be performed. The method described below is a method of performing the above extraction, and is a method of sampling the output of the interferometer 103 with the peak value of intensity modulation of the output light of the light source light 101.

光源101からは、周期的に波長が掃引された光が出力される。その光強度は、レーザ共振器内に挿入された櫛歯型フィルタの透過率とレーザ共振器の共振周波数に応じて、時間的に変化している。その光強度の極大値(ピーク値)の間隔は、時間をt、光源出力光の波数をk(t)、hを定数としたときの、cos(2hk(t))の正のピーク間隔と同じである。hの値は観測対象において観測したいもっとも深い位置についての情報によって設計するが、これについては後述する。   The light source 101 outputs light whose wavelength is periodically swept. The light intensity changes with time according to the transmittance of the comb filter inserted in the laser resonator and the resonance frequency of the laser resonator. The interval between the maximum values (peak values) of the light intensity is the positive peak interval of cos (2hk (t)) where time is t, the wave number of the light source output light is k (t), and h is a constant. The same. The value of h is designed by information about the deepest position to be observed in the observation target, which will be described later.

光源101からの出力光はカップラ102で分岐する。カップラ102から出た一方の光を干渉計103に入れ、サーキュレータ1032を経由後、カップラ1033にて光を2つに分け、参照アーム1034とサンプルアーム1036に入れる。参照アーム1034ではミラー1035にて反射する。サンプルアーム1036では観測対象の物体Oの中の反射面にて反射される。参照アーム1034とサンプルアーム1036から返ってきた光は再度カップラ1033に入って合波されて干渉光(フリンジ光)となる。この干渉光(フリンジ光)は一方はサーキュレータ1032を経由し、もう一方は直接バランスフォトディテクタ1031に入れ光電変換し、それによって生じた電気信号であるフリンジ信号をリサンプリング部107に入れる。   Output light from the light source 101 is branched by the coupler 102. One light emitted from the coupler 102 enters the interferometer 103, passes through the circulator 1032, is divided into two by the coupler 1033, and enters the reference arm 1034 and the sample arm 1036. The reference arm 1034 is reflected by the mirror 1035. The sample arm 1036 is reflected by the reflection surface in the object O to be observed. The light returned from the reference arm 1034 and the sample arm 1036 enters the coupler 1033 again and is combined to become interference light (fringe light). One of the interference light (fringe light) passes through the circulator 1032, and the other is directly input to the balance photodetector 1031 for photoelectric conversion, and the fringe signal which is an electric signal generated thereby is input to the resampling unit 107.

参照アーム1034とサンプルアーム1036との光路長差zは、参照アーム中のミラーの配置位置である基準面Aと物体O内の反射面Bとの間の距離によって決まる光路長差である。そのフリンジ信号には、光路長差zに比例する周波数を持つ信号成分が含まれている。たとえば、時間をtとしてk(t)という関数で光源101の出力光の波数を表すと、フリンジ信号にはcos(2zk(t))の信号が含まれる。つまり、時刻tにはzk’(t)/πとなる周波数の信号が含まれる。ここで、k’(t)はk(t)を時間微分した関数であり、k’(t)=dk/dtである。物体Oの中の反射面Bは、一般に1つとは限らず、複数ある場合がある。たとえば、z1、…、zNのN個の反射面がある場合は、フリンジ信号はcos(2z1k(t))、…、cos(2zNk(t))のN個の成分を持つこととなる。 The optical path length difference z between the reference arm 1034 and the sample arm 1036 is an optical path length difference determined by the distance between the reference plane A that is the position of the mirror in the reference arm and the reflecting surface B in the object O. The fringe signal includes a signal component having a frequency proportional to the optical path length difference z. For example, when the wave number of the output light of the light source 101 is expressed by a function k (t) where time is t, the fringe signal includes a signal of cos (2zk (t)). That is, at time t, a signal having a frequency zk ′ (t) / π is included. Here, k ′ (t) is a function obtained by time-differentiating k (t), and k ′ (t) = dk / dt. In general, the reflection surface B in the object O is not limited to one and may be plural. For example, if there are N reflecting surfaces of z 1 ,..., Z N , the fringe signal has N components of cos (2z 1 k (t)), ..., cos (2z N k (t)). Will have.

カップラ102から出たもう一方の光はフォトディテクタ104が光電変換し、ピーク時刻取得部1106へ渡す。ピーク時刻取得部1106は、入力した電気信号のピークの時刻を測定してリサンプリング部107に出力する。   The other light emitted from the coupler 102 is photoelectrically converted by the photodetector 104 and passed to the peak time acquisition unit 1106. The peak time acquisition unit 1106 measures the peak time of the input electric signal and outputs it to the resampling unit 107.

図2は図1に示す構成のSS−OCTシステムの、フォトディテクタ104と干渉計103から出力された電気信号の処理手順である。   FIG. 2 shows a procedure for processing electrical signals output from the photodetector 104 and the interferometer 103 in the SS-OCT system having the configuration shown in FIG.

非特許文献1に記載のピーク時刻取得方法は、ピーク検出対象のデータ(フォトディテクタ104から出力された電気信号)をアップサンプリングした上でローパスフィルタをかけるものであり、高周波ノイズを除去してアップサンプリングのサンプリング間隔の時間精度でピーク時刻を取得する(S−1201)。   The peak time acquisition method described in Non-Patent Document 1 applies a low-pass filter after up-sampling peak detection target data (electric signal output from the photodetector 104), and removes high-frequency noise to up-sample. The peak time is acquired with the time accuracy of the sampling interval (S-1201).

リサンプリング部107は、ピーク時刻取得部106から入力したピーク時刻に従って、バランスフォトディテクタ1031から入力したフリンジ信号をサンプリングし、フーリエ変換部108に出力する。(S−1202) フーリエ変換部108は、リサンプリング部107から入力した、フリンジ信号のサンプリング信号をフーリエ変換して出力する。(S−1203)
リサンプリング部107はフリンジ信号に含まれるcos(2zk(t))の成分がcos(2z(Δk/ΔT)t)となるようにフリンジ信号を時間的に伸縮するものである。ここで、Δkは光源の掃引波数幅、ΔTはΔkだけ掃引する時間である。リサンプリング前のフリンジ信号にはzk’(t)/πで表される周波数を含み、k’(t)がtに対して線形でなければ(tの一次関数でなければ)時間によって周波数が変わる。リサンプリングはフリンジ信号をz(Δk/ΔT)t/πという時間変動しない周波数を持つ信号になるように、フリンジ信号を時間的に伸縮処理するものである。これは、光源101の出力光がピークとなる時刻がcos(2hk(t))の正のピークとなる時刻となる性質を使って、そのピーク時刻でフリンジ信号をサンプリングことにより得られる。このとき、リサンプリングによってフリンジ信号は、時間当たりΔk/ΔTとなるような、kの時間変動が時間に対して線形となるように変換される。このようにkが時間に対して線形になることを波数リニアと呼ぶ。
The resampling unit 107 samples the fringe signal input from the balance photodetector 1031 according to the peak time input from the peak time acquisition unit 106 and outputs the sampled signal to the Fourier transform unit 108. (S-1220) The Fourier transform unit 108 performs Fourier transform on the sampling signal of the fringe signal input from the resampling unit 107 and outputs the result. (S-1203)
The resampling unit 107 temporally expands and contracts the fringe signal so that the component of cos (2zk (t)) included in the fringe signal becomes cos (2z (Δk / ΔT) t). Here, Δk is a sweep wave number width of the light source, and ΔT is a time for sweeping by Δk. The fringe signal before resampling includes a frequency represented by zk ′ (t) / π, and if k ′ (t) is not linear with respect to t (if it is not a linear function of t), the frequency varies with time. change. Resampling is to temporally expand and contract the fringe signal so that the fringe signal becomes a signal having a frequency that does not vary with time, such as z (Δk / ΔT) t / π. This is obtained by sampling the fringe signal at the peak time using the property that the time when the output light of the light source 101 reaches the peak is the positive peak of cos (2hk (t)). At this time, the resampling converts the fringe signal so that the time variation of k is linear with respect to time, such that Δk / ΔT per time. The fact that k is linear with respect to time is called wave number linear.

光源101の出力光の強度変調の周期を決めるhは、観測対象において観測したいもっとも深い位置における、参照アーム1034とサンプルアーム1036の光路長差zmaxによって決める。リサンプリング部107は、光源101の出力光の強度のピーク時刻でフリンジ信号をサンプリングするが、zに関するフリンジ光の強度変調の時間tにおける周波数はzk’(t)/πであり、光源出力光の強度変調のそれはhk’(t)/πであるため、zが0からzmaxとなる全てのフリンジ信号cos(2zk(t))(ただし、0≦z≦zmax)を取得するために、サンプリング定理に基づいてサンプリング周波数はzmaxの周波数の2倍以上、つまり、2zmaxk’(t)/π≦hk’(t)/π、すなわち、2z≦hとなるようにhを決める。光源のhをどのように実現するかについては、たとえば、非特許文献1では、ファイバレーザ中にファイバファブリーペローフィルタ(fiber Fabry−Perot filter:FFPF)のFSRがc/(2h)となるように設計した例がある。ここでcは光速である。 H that determines the period of intensity modulation of the output light of the light source 101 is determined by the optical path length difference z max between the reference arm 1034 and the sample arm 1036 at the deepest position to be observed in the observation target. The resampling unit 107 samples the fringe signal at the peak time of the intensity of the output light from the light source 101. The frequency at the time t of the intensity modulation of the fringe light with respect to z is zk ′ (t) / π, and the light source output light In order to obtain all the fringe signals cos (2zk (t)) (where 0 ≦ z ≦ z max ) where z is from 0 to z max , it is hk ′ (t) / π. Based on the sampling theorem, the sampling frequency is determined to be at least twice the frequency of z max , that is, 2z max k ′ (t) / π ≦ hk ′ (t) / π, that is, 2z ≦ h. . Regarding how to realize h of the light source, for example, in Non-Patent Document 1, the FSR of a fiber Fabry-Perot filter (FFPF) in a fiber laser is c / (2h). There is an example designed. Here, c is the speed of light.

フーリエ変換部108から出力された信号は、SS−OCTの深さ方向のスキャン信号であり、この深さ信号のスキャンのことをAスキャンという。また、このスキャン信号のことをAスキャン信号と呼ぶこととする。Aスキャンとは光偏向器等で偏向する等、物理的な光の操作を意味するのではなく、上記の通り、フリンジ信号をフーリエ変換することにより奥行方向の情報を得ることを言う。   The signal output from the Fourier transform unit 108 is a scan signal in the depth direction of SS-OCT, and this depth signal scan is referred to as A scan. Also, this scan signal is called an A scan signal. The A-scan does not mean physical light operation such as deflecting with an optical deflector or the like, but means obtaining depth direction information by Fourier-transforming the fringe signal as described above.

Aスキャン信号はリサンプリングされたフリンジ信号のフーリエ変換であるため、周波数に対する強度を示す信号であるが、周波数fは深さ方向の距離zに換算でき、それらの関係はz=πf/k’と表される。ここで、k’はリサンプリング後のAスキャン信号においては、理想的にはk’=Δk/ΔTとなるので、リサンプリング後はz=πfΔT/Δkとなる。   Since the A scan signal is a Fourier transform of the resampled fringe signal, it is a signal indicating the intensity with respect to the frequency, but the frequency f can be converted to a distance z in the depth direction, and their relationship is z = πf / k ′. It is expressed. Here, k ′ is ideally k ′ = Δk / ΔT in the A-scan signal after resampling, and therefore z = πfΔT / Δk after resampling.

観測対象の物体O内の反射面は、Aスキャン信号では信号強度で表されるが、Aスキャン信号上では、反射面がある位置に相当する周波数の信号が強くなる。反射面が複数ある場合は、それら反射面に対応する、参照アームとサンプルアームの光路長差に応じた成分の信号がフリンジ信号に含まれる。たとえば、N個の反射点があり、それらに対応する光路長差がz1、…、zNであった場合、フリンジ信号にはcos(2z1k(t))、…、cos(2zNk(t))のN個の成分が含まれる。リサンプリングされたフリンジ信号にも同様にcos(2z1Δk/ΔTt)、…、cos(2zNΔk/ΔTt)のN個の成分が含まれ、それらの周波数はf1=z1Δk/(πΔT)、…、fN=zNΔk/(πΔT)となる。このリサンプリングされたフリンジ信号をフーリエ変換すると、観測対象の物体O内の反射強度に応じて、上記f1、…、fNに相当する周波数の位置の強度が大きい信号が現れる。これらの周波数をz=πfΔT/Δkで換算すると、奥行距離と反射の強度の対応を取ることができる。 The reflection surface in the object O to be observed is represented by the signal intensity in the A scan signal, but on the A scan signal, a signal having a frequency corresponding to the position where the reflection surface is located becomes strong. When there are a plurality of reflection surfaces, the fringe signal includes a component signal corresponding to the reflection surface and corresponding to the difference in optical path length between the reference arm and the sample arm. For example, there are N of the reflection point, the optical path length difference z 1 corresponding to them, ..., if it was z N, the fringe signal cos (2z 1 k (t) ), ..., cos (2z N N components of k (t)) are included. Similarly, the resampled fringe signal includes N components of cos (2z 1 Δk / ΔTt),..., Cos (2z N Δk / ΔTt), and their frequencies are f 1 = z 1 Δk / ( πΔT),..., f N = z N Δk / (πΔT). When the resampled fringe signal is Fourier-transformed, a signal having a high intensity at a frequency position corresponding to f 1 ,..., F N appears in accordance with the reflection intensity in the object O to be observed. When these frequencies are converted into z = πfΔT / Δk, the correspondence between the depth distance and the reflection intensity can be taken.

観測対象の物体O内の反射面が1つのみであり、その反射面がこの時の参照光とサンプル光の光路長差zに相当する箇所にある場合は、Aスキャン信号のzに相当する箇所、つまり、周波数f=z(Δk/ΔT)/πの信号が強くなり、それ以外の周波数の信号は0(デルタ関数)になることが望ましい。このような、1つの反射点を持つ場合のフリンジ信号から得られるAスキャン信号を点広がり関数(point spread function:PSF)と呼ぶ。しかし、ΔT内での光源光強度の変調や、光学的および電気的なノイズ、リサンプリング等の信号処理の不十分さ等が原因で、実際のPSFはデルタ関数とはならず、PSFの半値幅は大きくなり、また、SNR(PSFのピーク/ノイズフロアレベル)は低く(悪く)なる。   When there is only one reflecting surface in the object O to be observed and the reflecting surface is at a position corresponding to the optical path length difference z between the reference light and the sample light at this time, it corresponds to z of the A scan signal. It is desirable that the signal at the location, that is, the frequency f = z (Δk / ΔT) / π is strong, and the signal at the other frequency is 0 (delta function). Such an A-scan signal obtained from a fringe signal in the case of having one reflection point is called a point spread function (PSF). However, due to modulation of the light source light intensity within ΔT, optical and electrical noise, inadequate signal processing such as resampling, etc., the actual PSF does not become a delta function, and half of the PSF. The value range becomes large, and the SNR (PSF peak / noise floor level) becomes low (bad).

リサンプリングとは、波数リニアでないフリンジ信号を波数リニアとなるようにする処理であるが、サンプリングする時刻に誤差が生じると(正確なリサンプリングができないと)、正確な波数リニアなフリンジ信号を生成することはできなくなり、PSFの半値幅は大きくなり、また、PSFのSNRは悪化する。   Resampling is a process that makes a wavenumber linear that is not wavenumber linear, but if there is an error in the sampling time (if accurate resampling is not possible), an accurate wavenumber fringe signal is generated. The half width of the PSF becomes large, and the SNR of the PSF deteriorates.

Tsung-Han Tsaia, Chao Zhoua, Desmond Adlera, and James G. Fujimotoa, “Frequency Comb Swept Lasers for Optical Coherence Tomography,” Proc. of SPIE , Vol. 7554, pp. 7554E-1-10 (2010).Tsung-Han Tsaia, Chao Zhoua, Desmond Adlera, and James G. Fujimotoa, “Frequency Comb Swept Lasers for Optical Coherence Tomography,” Proc. Of SPIE, Vol. 7554, pp. 7554E-1-10 (2010).

ピーク時刻取得部1106において光源光の強度変調信号のピークが現れる時刻を取得する際のピーク時刻の検出精度はサンプリング周波数に依存する。したがってサンプリング周波数を高くすればピーク時刻の検出精度は向上するが、その分ピーク時刻検出処理に多大なリソースを要することになる。一方、低いサンプリング周波数でピーク時刻の検出を行うと、ピーク時刻検出精度が劣化してしまう。   The detection accuracy of the peak time when acquiring the time at which the peak of the intensity modulation signal of the light source light appears in the peak time acquisition unit 1106 depends on the sampling frequency. Therefore, if the sampling frequency is increased, the detection accuracy of the peak time is improved, but much resources are required for the peak time detection process. On the other hand, if peak time is detected at a low sampling frequency, the peak time detection accuracy deteriorates.

すなわち、ピーク時刻検出精度が低い場合は、その後の処理で正確なリサンプリングが難しくなるため、光源101の出力光のコヒーレンス長の長い光で干渉計103で干渉させて得たフリンジ信号のみを抽出することが難しくなるので、Aスキャン信号の深達長が長くならない、つまり、断層画像が深くまで見ることができないという問題がある。また、同じ原因により、フーリエ変換部108から出力されるAスキャン信号のPSFの半値幅が大きくなりやすいという問題と、SNRが悪くなりやすいという問題がある。   That is, when the peak time detection accuracy is low, it is difficult to perform accurate resampling in the subsequent processing. Therefore, only the fringe signal obtained by causing the interferometer 103 to interfere with light having a long coherence length of the output light of the light source 101 is extracted. Therefore, there is a problem that the depth of the A scan signal does not become long, that is, the tomographic image cannot be seen deeply. Further, due to the same cause, there is a problem that the half-value width of the PSF of the A scan signal output from the Fourier transform unit 108 is likely to be large, and a problem that the SNR is likely to be deteriorated.

本発明は上記従来の問題に鑑みてなされたものであり、本発明の課題は、簡素な処理で正確なリサンプリングを行うことができる光干渉断層装置を提供することにある。   The present invention has been made in view of the above-described conventional problems, and an object of the present invention is to provide an optical coherence tomography apparatus that can perform accurate resampling with a simple process.

上記の課題を解決するために、一実施形態に記載された発明は、櫛歯型のスペクトルとなるように波長掃引する光源と、前記光源から出力された光を2分岐する第1のカップラと、前記2分岐された一方の光を観測対象の物体に入射して観測対象の物体の構造に応じた干渉信号であるフリンジ信号を出力する干渉計と、前記2分岐された他方の光の強度信号の位相を算出する位相算出部と、前記位相算出部で算出された位相が2nπとなる時刻であるピーク時刻を取得する位相2nπ時刻取得部と、前記ピーク時刻におけるフリンジ信号をサンプリングするリサンプリング部とを備えることを特徴する光干渉断層装置である。   In order to solve the above-described problem, the invention described in one embodiment includes a light source that sweeps a wavelength so as to have a comb-shaped spectrum, and a first coupler that splits the light output from the light source into two. , An interferometer that makes one of the two branched lights incident on an object to be observed and outputs a fringe signal that is an interference signal according to the structure of the object to be observed; and the intensity of the other two branched lights A phase calculation unit that calculates the phase of the signal, a phase 2nπ time acquisition unit that acquires a peak time that is a time when the phase calculated by the phase calculation unit is 2nπ, and resampling that samples the fringe signal at the peak time And an optical coherence tomography apparatus.

従来のSS−OCTシステムの構成例を示す図である。It is a figure which shows the structural example of the conventional SS-OCT system. 図1のSS−OCTシステムのフォトディテクタから出力された信号の処理フローを示す図である。It is a figure which shows the processing flow of the signal output from the photodetector of the SS-OCT system of FIG. 第1の実施形態のSS−OCTシステムの構成例を示す図である。It is a figure which shows the structural example of the SS-OCT system of 1st Embodiment. 空間光学系で構成した光源の構成例を示す図である。It is a figure which shows the structural example of the light source comprised by the spatial optical system. 図3の装置のフォトディテクタから出力された信号の処理フローを示す図である。It is a figure which shows the processing flow of the signal output from the photodetector of the apparatus of FIG. 強度信号と位相の経時変化を示す図である。It is a figure which shows a time-dependent change of an intensity | strength signal and a phase.

以下、本発明の実施の形態について、詳細に説明する。本発明の光干渉断層装置は、櫛歯型のスペクトルとなるように波長掃引する光源と、前記光源から出力された光を2分岐する第1のカップラと、前記2分岐された一方の光を観測対象の物体に入射して観測対象の物体の構造に応じた干渉信号であるフリンジ信号を出力する干渉計と、前記2分岐された他方の光の強度信号の位相を算出する位相算出部と、前記位相算出部で算出された位相が2nπとなる時刻であるピーク時刻を取得する位相2nπ時刻取得部と、前記ピーク時刻におけるフリンジ信号をサンプリングするリサンプリング部とを備えている。   Hereinafter, embodiments of the present invention will be described in detail. The optical coherence tomography apparatus of the present invention includes a light source that sweeps a wavelength so as to obtain a comb-shaped spectrum, a first coupler that bifurcates light output from the light source, and one of the two branched lights. An interferometer that is incident on the object to be observed and outputs a fringe signal that is an interference signal corresponding to the structure of the object to be observed; a phase calculator that calculates the phase of the intensity signal of the other two branched light; A phase 2nπ time acquisition unit that acquires a peak time, which is a time when the phase calculated by the phase calculation unit becomes 2nπ, and a resampling unit that samples a fringe signal at the peak time.

上記光干渉断層装置によれば、光源出力光のピーク時刻を正確に得られるため、正確なリサンプリングを行うことができる。これにより、光源の出力光強度のピーク時刻に現れるコヒーレンス長の長い光の干渉信号を正確にサンプリングすることが可能となり、断層画像の深達長を長くすること、つまり、観察対象の物体を深くまで観測することが可能となる。   According to the optical coherence tomography apparatus, since the peak time of the light source output light can be obtained accurately, accurate resampling can be performed. This makes it possible to accurately sample the interference signal of light with a long coherence length that appears at the peak time of the output light intensity of the light source, and to increase the depth of the tomographic image, that is, to deepen the object to be observed. It becomes possible to observe up to.

また、上記の通り、正確なリサンプリングが可能となることにより、リサンプリング部から出力されるフリンジ信号は光源が正確に波数リニアな掃引した場合のフリンジ信号と同一となるため、フーリエ変換部から出力されるAスキャン信号のPSFの半値幅は細くなり、また、SNRが大きくなるため、鮮明度は高くなる。   In addition, as described above, since accurate resampling is possible, the fringe signal output from the resampling unit is the same as the fringe signal when the light source is accurately swept in a wavenumber linear form. Since the half width of the PSF of the output A scan signal becomes narrower and the SNR becomes larger, the sharpness becomes higher.

(第1の実施形態)
第1の実施形態の光干渉断層装置について説明する。図3は、第1の実施形態のSS−OCTシステムの構成例を示す図である。本実施形態のSS−OCTシステムは、図3に示すように、櫛歯型のスペクトルとなるように波長掃引を行なう光源101と、カップラ102と、干渉計103と、フォトディテクタ104と、位相算出部105と、位相2nπ時刻取得部106と、リサンプリング部107と、フーリエ変換部108とを備えている。さらに干渉計103は、バランスフォトディテクタ1031と、サーキュレータ1032と、カップラ1033と、ファイバコリメータ1038と、ミラー1035と、ファイバコリメータ1039と、光偏向器1037とを有している。カップラ1033からファイバコリメータ1038を通ってミラー1035までの光路が参照アーム1034であり、カップラ1033からファイバコリメータ1039と光偏向器1037を通って観察対象Oの中の反射面Bまでの光路がサンプルアーム1036である。図中の各構成は光ファイバ(二重線)または電気ケーブル(一本線)で接続されている。
(First embodiment)
The optical coherence tomography apparatus according to the first embodiment will be described. FIG. 3 is a diagram illustrating a configuration example of the SS-OCT system according to the first embodiment. As shown in FIG. 3, the SS-OCT system of the present embodiment includes a light source 101 that performs wavelength sweeping so as to obtain a comb-shaped spectrum, a coupler 102, an interferometer 103, a photodetector 104, and a phase calculation unit. 105, a phase 2nπ time acquisition unit 106, a resampling unit 107, and a Fourier transform unit 108. Further, the interferometer 103 includes a balance photodetector 1031, a circulator 1032, a coupler 1033, a fiber collimator 1038, a mirror 1035, a fiber collimator 1039, and an optical deflector 1037. The optical path from the coupler 1033 through the fiber collimator 1038 to the mirror 1035 is a reference arm 1034, and the optical path from the coupler 1033 through the fiber collimator 1039 and the optical deflector 1037 to the reflecting surface B in the observation object O is the sample arm. 1036. Each component in the figure is connected by an optical fiber (double line) or an electric cable (single line).

光源101は、時間的に波長(周波数)が推移する(掃引する)レーザであり、その発振スペクトルのピークは周波数に対して等間隔に離散的な値をとる。その発振スペクトルを1掃引分重ね合わせた(1掃引分を時間積分した)スペクトルは、周波数に対して等間隔にピークがある櫛歯型のスペクトルとなる。この光源101は、周波数が時間的に推移するので、レーザ発振は時間的に強くなったり弱くなったりを繰り返す。つまりこの光源101の出力光は強度変調する。この強度変調はcos(2hk(t))の基本成分とその高調波成分を持つ。ここで、cos(x)はxの余弦、hは定数、k(t)は光の波数であり、時間tにより変化する。このhは光源内のレーザ共振器中にある光学部品のパラメータであるが、測定対象の物体の最大深さに基づいて決定される定数である。   The light source 101 is a laser whose wavelength (frequency) changes (sweeps) in time, and the peak of its oscillation spectrum takes discrete values at equal intervals with respect to the frequency. A spectrum obtained by superimposing the oscillation spectrum for one sweep (a time integral of one sweep) is a comb-shaped spectrum having peaks at equal intervals with respect to the frequency. Since the frequency of the light source 101 changes with time, the laser oscillation repeatedly increases and decreases with time. That is, the intensity of the output light from the light source 101 is modulated. This intensity modulation has a fundamental component of cos (2hk (t)) and its harmonic component. Here, cos (x) is the cosine of x, h is a constant, k (t) is the wave number of light, and changes with time t. This h is a parameter of the optical component in the laser resonator in the light source, and is a constant determined based on the maximum depth of the object to be measured.

強度変調がcos(2hk(t))となる光源101は、その光源101で生じる櫛歯型のスペクトルの隣同士の櫛歯の周波数間隔(自由スペクトル間隔:FSR)がc/(2h)となる光源である。ここでcは光速である。その光源101の構成としては、ファイバレーザ中にファイバファブリーペローフィルタ(fiber Fabry−Perot filter:FFPF)を挿入したものや、空間光学系で構成したレーザ共振器内にエタロンを挿入したものがある。このときのFFPFやエタロンはそのFSRがc/(2h)となるものを使用することで、FSRがc/(2h)となる光源が実現する。   In the light source 101 whose intensity modulation is cos (2hk (t)), the frequency interval (free spectral interval: FSR) between adjacent comb teeth spectrum generated by the light source 101 is c / (2h). Light source. Here, c is the speed of light. As a configuration of the light source 101, there are a fiber laser in which a fiber Fabry-Perot filter (FFPF) is inserted, and an etalon inserted in a laser resonator constituted by a spatial optical system. The FFPF or etalon at this time uses a light source having an FSR of c / (2h), thereby realizing a light source having an FSR of c / (2h).

図4に空間光学系で構成した光源1の構成例を示す。図4に示す構成は空間光学系で構成されたLittman―Metcalf型共振器であり、半導体光増幅器(semiconductor optical amplifiler:SOA)304から発した光は、コリメータレンズ306でコリメートされた後、エタロン307を通して光偏向器308で偏向されて回折格子310に入射し、回折した光をミラー311に送る。ミラー311に垂直入射された光のみが、逆の経路をたどってSOA304に戻り、SOA反射面305で折り返され、最終的にはSOA反射面305とミラー311の間でレーザ共振する。SOA反射面305はハーフミラーとなっており、レーザ発振した光の一部がSOA反射面を通って結合レンズ303、アイソレータ302を介して光ファイバ301に至り、共振器外部に出力される。なお、エタロンの光路長はhである物を使う。たとえば、エタロンの屈折率をn、エタロンの厚みをdとすると、h=ndとなる。レーザ共振器の光の共振周波数は、光偏向器308から回折格子310への入射角度によって変化するので、光偏向器制御器309によって光偏向器308から出射する光の角度(偏向角)を時間的に連続的に変えることによって、光源101の出力光の周波数を連続的に変える(掃引する)。   FIG. 4 shows a configuration example of the light source 1 configured by a spatial optical system. The configuration shown in FIG. 4 is a Littman-Metcalf type resonator composed of a spatial optical system. Light emitted from a semiconductor optical amplifier (SOA) 304 is collimated by a collimator lens 306 and then etalon 307. Then, the light is deflected by the optical deflector 308 and incident on the diffraction grating 310, and the diffracted light is sent to the mirror 311. Only the light that is perpendicularly incident on the mirror 311 returns to the SOA 304 along the reverse path, is folded back at the SOA reflecting surface 305, and finally laser-resonates between the SOA reflecting surface 305 and the mirror 311. The SOA reflecting surface 305 is a half mirror, and part of the laser-oscillated light passes through the SOA reflecting surface, reaches the optical fiber 301 via the coupling lens 303 and the isolator 302, and is output to the outside of the resonator. An etalon having an optical path length of h is used. For example, when the refractive index of the etalon is n and the thickness of the etalon is d, h = nd. Since the resonance frequency of the light of the laser resonator varies depending on the incident angle from the optical deflector 308 to the diffraction grating 310, the angle (deflection angle) of the light emitted from the optical deflector 308 by the optical deflector controller 309 is changed over time. Therefore, the frequency of the output light from the light source 101 is continuously changed (swept).

カップラ102は、光源101から入力された光を2分岐し、一方を干渉計103に入力し、他方をフォトディテクタ104に入力する。   The coupler 102 splits the light input from the light source 101 into two parts, inputs one to the interferometer 103, and inputs the other to the photodetector 104.

カップラ102で分岐された一方の光が入力される干渉計103では、入力された光をサーキュレータ1032を介してカップラ1033に入れ、カップラ1033で2分岐した後、参照アーム1034とサンプルアーム1036に入力する。カップラ1033で分岐された一方の光は、参照アーム1034の基準位置Aに設けられたミラー1035で反射して再びカップラ1033に戻る。カップラ1033で分岐されたもう一方の光は、観測対象の物体O内の反射面B(反射面Bは物体内の一か所とは限らない)で反射されて再びカップラ1033に戻る。光源が1掃引するごとに光偏向器1037で幅方向WにAスキャンを行う位置を少しずつ移動させる。このように、物理的に光を観測対象にあてる位置(Aスキャンを行う位置)を連続的に変える操作をBスキャンという。   In the interferometer 103 to which one of the lights branched by the coupler 102 is input, the input light is input to the coupler 1033 via the circulator 1032 and branched into two by the coupler 1033 and then input to the reference arm 1034 and the sample arm 1036. To do. One light branched by the coupler 1033 is reflected by the mirror 1035 provided at the reference position A of the reference arm 1034 and returns to the coupler 1033 again. The other light branched by the coupler 1033 is reflected by the reflection surface B in the object O to be observed (the reflection surface B is not necessarily one place in the object) and returns to the coupler 1033 again. Each time the light source sweeps one time, the position where the A-scan is performed in the width direction W by the optical deflector 1037 is gradually moved. In this way, the operation of continuously changing the position where light is physically applied to the observation target (the position where the A scan is performed) is called B scan.

カップラ1033では参照アーム1034およびサンプルアーム1036からの2つの戻り光を合波して干渉させ、その干渉光(フリンジ光)を、お互いに位相が反転した状態で2分岐し、一方はサーキュレータ1032を介して、他方は直接バランスフォトディテクタ1031に入射する。バランスフォトディテクタ1031は、それぞれの入射光を電気信号に変換し、それらの差分信号をリサンプリング部107に出力する。この出力信号を干渉信号(フリンジ信号)と呼ぶこととする。   In the coupler 1033, the two return lights from the reference arm 1034 and the sample arm 1036 are combined to interfere with each other, and the interference light (fringe light) is branched into two with the phases inverted from each other. The other is directly incident on the balance photodetector 1031. The balance photodetector 1031 converts each incident light into an electrical signal and outputs the difference signal to the resampling unit 107. This output signal is called an interference signal (fringe signal).

干渉計103から出力されるフリンジ信号(電気信号)は、参照アームとサンプルアームとの光路長差zに応じて変化する信号である。具体的には、フリンジ信号はcos(2zk(t))の成分を持つ。また、光源101からの出力光はcos(2hk(t))で強度変調しているので、フリンジ信号はcos(2zk(t))およびcos(2hk(t))で表される2つの基本成分とその高調波成分を持つこととなる。観測対象の物体O内の反射面の基準面Aからの最深の位置zをzmaxとすると、光源101の強度変調の周波数を決定し、また、リサンプリングする時間間隔を決定するhは2zmax以上(h≧2zmax)となるようにする。このようにする理由は次に示すような原理による。zに対するフリンジ信号はcos(2zk(t))で変調されており、その時刻tにおける周波数fz(t)はzk’(t)/πとなる。ただし、k’(t)は時刻tにおけるk(t)の微分した関数でありk’(t)=dk/dtである。一方、光源101の光出力はcos(2hk(t))で変調されており、その時刻tにおける周波数fh(t)はhk’(t)/πとなる。光源101の光出力のピークの時刻でフリンジ信号をサンプリングするので、時刻tあたりではhk’(t)/πで表されるサンプリング周波数でフリンジ信号をサンプリングすることになる。観測対象の物体O内の反射点によるフリンジの最高周波数はzmaxとなる反射点からの光で生じる干渉光の周波数zmaxk’(t)/πであるので、この周波数の信号成分の情報をすべて取得するには、サンプリング定理から、2zmaxk’(t)/π≦hk’(t)/πである必要がある。したがって、hは2zmax以上(、h≧2zmax)となるようにする。 The fringe signal (electric signal) output from the interferometer 103 is a signal that changes according to the optical path length difference z between the reference arm and the sample arm. Specifically, the fringe signal has a component of cos (2zk (t)). Since the output light from the light source 101 is intensity-modulated by cos (2hk (t)), the fringe signal has two basic components represented by cos (2zk (t)) and cos (2hk (t)). And its harmonic components. If z max is the deepest position z from the reference plane A of the reflecting surface in the object O to be observed, the frequency of intensity modulation of the light source 101 is determined, and h for determining the time interval for resampling is 2z max More than (h ≧ 2z max ). The reason for this is based on the following principle. The fringe signal for z is modulated by cos (2zk (t)), and the frequency fz (t) at time t is zk ′ (t) / π. However, k ′ (t) is a differentiated function of k (t) at time t, and k ′ (t) = dk / dt. On the other hand, the light output of the light source 101 is modulated by cos (2hk (t)), and the frequency fh (t) at the time t is hk ′ (t) / π. Since the fringe signal is sampled at the peak time of the light output of the light source 101, the fringe signal is sampled at a sampling frequency represented by hk ′ (t) / π around time t. Since the highest frequency of the fringes due to reflection points in the object O to be observed is the z frequency of the interference light generated by the light from the max become reflection point z max k '(t) / π, the information of the signal component of the frequency In order to obtain all of the above, it is necessary from the sampling theorem that 2z max k ′ (t) / π ≦ hk ′ (t) / π. Therefore, h is set to 2z max or more (where h ≧ 2z max ).

カップラ102で分岐された他方の光が入力されるフォトディテクタ104では、入力された光を光強度に応じた光源強度信号(電気信号)に変換し、位相算出部105に送る。   In the photodetector 104 to which the other light branched by the coupler 102 is input, the input light is converted into a light source intensity signal (electrical signal) corresponding to the light intensity and sent to the phase calculation unit 105.

図5は、位相算出部105、位相2nπ時刻取得部106、リサンプリング部107によるフォトディテクタ104と干渉計103とから出力された電気信号に対する処理手順である。   FIG. 5 is a processing procedure for the electrical signals output from the photodetector 104 and the interferometer 103 by the phase calculation unit 105, the phase 2nπ time acquisition unit 106, and the resampling unit 107.

位相算出部105は、フォトディテクタ104から出力された光源光強度信号を入力し、信号の時刻ごとの位相を算出し、位相2nπ時刻取得部106に出力する(S−201)。   The phase calculation unit 105 receives the light source light intensity signal output from the photodetector 104, calculates the phase of each signal time, and outputs it to the phase 2nπ time acquisition unit 106 (S-201).

位相算出の一例を以下に示す。ただし、フォトディテクタ出力信号が光源強度に比例する場合は、位相算出部105の入力信号はcos(2hk(t))に比例するが、説明を簡単にするため、ここでは余弦関数の係数を1として、つまり、入力信号s(t)をcos(2hk(t))として位相算出手段とその原理を以下に説明する。   An example of phase calculation is shown below. However, when the photodetector output signal is proportional to the light source intensity, the input signal of the phase calculation unit 105 is proportional to cos (2hk (t)), but for the sake of simplicity of explanation, the coefficient of the cosine function is set to 1 here. That is, the phase calculation means and its principle will be described below with the input signal s (t) as cos (2hk (t)).

1)入力信号s(t)をフーリエ変換する。その結果S(ω)は以下のようになる。   1) Fourier transform the input signal s (t). As a result, S (ω) is as follows.

Figure 0006259370
Figure 0006259370

ここで、   here,

Figure 0006259370
Figure 0006259370

はフーリエ変換を示し、jは虚数単位を表す。上記式から、S(ω)は時刻tにおいて、±zk’(t)/πの2つの周波数(±2zk’(t)の2つの角周波数)を持つ信号であることが分かる。 Indicates a Fourier transform, and j indicates an imaginary unit. From the above equation, it can be seen that S (ω) is a signal having two frequencies ± zk ′ (t) / π (two angular frequencies ± 2zk ′ (t)) at time t.

2)S(ω)のマイナスの周波数成分を削除する。つまり、S(ω)の   2) Delete the negative frequency component of S (ω). That is, S (ω)

Figure 0006259370
Figure 0006259370

を削除する。マイナスの周波数成分が削除された信号をS’(ω)とすると、S’(ω)は以下のようになる。 Is deleted. Assuming that a signal from which the negative frequency component is deleted is S ′ (ω), S ′ (ω) is as follows.

Figure 0006259370
Figure 0006259370

3)S’(ω)を逆フーリエ変換する。その結果をs’(t)とするとs’(t)は以下のようになる。   3) Inverse Fourier transform S ′ (ω). If the result is s '(t), s' (t) is as follows.

Figure 0006259370
Figure 0006259370

4)s’(t)の位相を算出する。具体的な位相の算出方法としては、たとえば、s’(t)の実部と虚部から位相を算出する。算出した位相は、上式中の2zk(t)と一致する。   4) The phase of s ′ (t) is calculated. As a specific method for calculating the phase, for example, the phase is calculated from the real part and the imaginary part of s ′ (t). The calculated phase coincides with 2zk (t) in the above equation.

位相2nπ時刻取得部106は、位相算出部105から得た位相が2nπとなる(ただしnは整数)となる時刻をピーク時刻としてリサンプリング部107に出力する。位相が2nπとなる時刻について説明すると、光源強度信号はcos(2zk(t))に比例するため、2zk(t)=2nπとなる時刻tとは光源強度は正のピークとなる時刻のことである(S−202)。   The phase 2nπ time acquisition unit 106 outputs the time when the phase obtained from the phase calculation unit 105 becomes 2nπ (where n is an integer) to the resampling unit 107 as a peak time. The time when the phase becomes 2nπ will be described. Since the light source intensity signal is proportional to cos (2zk (t)), the time t when 2zk (t) = 2nπ is the time when the light source intensity reaches a positive peak. Yes (S-202).

図6は強度信号s(t)と位相2zk(t)を示したものである。ただし、説明を簡単化するために、位相を2πでラッピングして表記おり、図6では位相は−π〜πの範囲で示されている。したがって、図6では位相2nπは0で示されている。図6に示すように、信号の正のピーク位置において位相が0となっていることが分かる。したがって、位相が0、すなわち2nπとなる時刻を強度信号s(t)が正のピーク位置となるピーク時刻として検出することができる。   FIG. 6 shows the intensity signal s (t) and the phase 2zk (t). However, in order to simplify the description, the phase is shown by being wrapped by 2π, and in FIG. 6, the phase is shown in a range of −π to π. Therefore, the phase 2nπ is shown as 0 in FIG. As shown in FIG. 6, it can be seen that the phase is 0 at the positive peak position of the signal. Therefore, the time when the phase is 0, that is, 2nπ, can be detected as the peak time when the intensity signal s (t) is at the positive peak position.

図6では、強度信号s(t)に重畳したノイズのために、位相が0付近で行ったり来たりして、強度信号s(t)のピーク位置の付近で複数回位相が0となる現象がでているが、この場合は、ローパスフィルタ等で強度信号s(t)から高周波ノイズを除去することにより、ピーク位置付近で位相を1回のみ0にすることができる。   In FIG. 6, due to the noise superimposed on the intensity signal s (t), the phase goes back and forth near 0, and the phase becomes 0 several times near the peak position of the intensity signal s (t). However, in this case, the phase can be set to 0 only once near the peak position by removing high-frequency noise from the intensity signal s (t) with a low-pass filter or the like.

リサンプリング部107は、位相2nπ時刻取得部106で取得した時刻で、干渉計103から出力されたフリンジ信号をサンプリングする。(S−203)
フーリエ変換部108は、リサンプリング部107でサンプリングしたフリンジ信号をフーリエ変換する。このフーリエ変換後の信号がSS−OCTのAスキャン信号となる。(S−204)
The resampling unit 107 samples the fringe signal output from the interferometer 103 at the time acquired by the phase 2nπ time acquisition unit 106. (S-203)
The Fourier transform unit 108 performs a Fourier transform on the fringe signal sampled by the resampling unit 107. The signal after the Fourier transform becomes an A-scan signal of SS-OCT. (S-204)

以上のOCTシステムによれば、サンプリング周波数を高くすることなく、また、それによってピーク検出に多くの計算リソースを使用することなく、光源出力光のピーク時刻を正確に得られるため、簡素な処理で正確なリサンプリングを行うことができる。   According to the above OCT system, the peak time of the light source output light can be accurately obtained without increasing the sampling frequency and without using many calculation resources for peak detection. Accurate resampling can be performed.

101 光源
102 カップラ
103 干渉計
1031 バランスフォトディテクタ
1032 サーキュレータ
1033 カップラ
1034 参照アーム
1035 ミラー
1036 サンプルアーム
1037 光偏向器
1038 ファイバコリメータ
1039 ファイバコリメータ
104 フォトディテクタ
105 位相算出部
106 位相2nπ時刻取得部
107 リサンプリング部
108 フーリエ変換部
1106 ピーク時刻取得部
301 光ファイバ
302 アイソレータ
303 結合レンズ
304 SOA
305 SOA反射面
306 コリメータレンズ
307 エタロン
308 光偏向器
309 光偏向器制御器
310 回折格子
311 ミラー
DESCRIPTION OF SYMBOLS 101 Light source 102 Coupler 103 Interferometer 1031 Balanced photo detector 1032 Circulator 1033 Coupler 1034 Reference arm 1035 Mirror 1036 Sample arm 1037 Optical deflector 1038 Fiber collimator 1039 Fiber collimator 104 Photo detector 105 Phase calculation unit 106 Phase 2nπ time acquisition unit 107 Conversion unit 1106 Peak time acquisition unit 301 Optical fiber 302 Isolator 303 Coupled lens 304 SOA
305 SOA reflecting surface 306 Collimator lens 307 Etalon 308 Optical deflector 309 Optical deflector controller 310 Diffraction grating 311 Mirror

Claims (2)

櫛歯型のスペクトルとなるように波長掃引する光源と、
前記光源から出力された光を2分岐する第1のカップラと、
前記2分岐された一方の光を観測対象の物体に入射して観測対象の物体の構造に応じた干渉信号であるフリンジ信号を出力する干渉計と、
前記2分岐された他方の光の強度信号の時刻ごとの位相を算出する位相算出部と、
前記位相算出部で算出された位相が2nπとなる時刻であるピーク時刻を取得する位相2nπ時刻取得部と、
前記ピーク時刻におけるフリンジ信号をサンプリングするリサンプリング部とを備えることを特徴する光干渉断層装置。
A light source that sweeps the wavelength so that it has a comb-shaped spectrum;
A first coupler for bifurcating light output from the light source;
An interferometer that emits a fringe signal that is an interference signal in accordance with the structure of the object to be observed by making one of the two branched lights incident on the object to be observed;
A phase calculation unit that calculates a phase for each time of the intensity signal of the other branched light;
A phase 2nπ time acquisition unit that acquires a peak time that is a time when the phase calculated by the phase calculation unit becomes 2nπ;
An optical coherence tomography apparatus comprising: a resampling unit that samples a fringe signal at the peak time.
前記干渉計は、入力された光を2分岐し一方の光を参照アームに入力して基準位置Aに配置されたミラーで反射させ、他方の光をサンプルアームに入力して観測対象の物体の反射面Bで反射させた2つの光を合波して干渉した干渉光を出力する第2のカップラと、該干渉光をフリンジ信号に変換して出力するバランスフォトディテクタとを備えることを特徴とする請求項1記載の光干渉断層装置。   The interferometer splits the input light into two, inputs one light to the reference arm and reflects it by a mirror arranged at the reference position A, and inputs the other light to the sample arm to input the object to be observed. A second coupler that outputs interference light obtained by combining two lights reflected by the reflection surface B and outputs interference light, and a balance photodetector that converts the interference light into a fringe signal and outputs the fringe signal. The optical coherence tomography apparatus according to claim 1.
JP2014151157A 2014-07-24 2014-07-24 Optical coherence tomography device Active JP6259370B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2014151157A JP6259370B2 (en) 2014-07-24 2014-07-24 Optical coherence tomography device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2014151157A JP6259370B2 (en) 2014-07-24 2014-07-24 Optical coherence tomography device

Publications (2)

Publication Number Publication Date
JP2016024176A JP2016024176A (en) 2016-02-08
JP6259370B2 true JP6259370B2 (en) 2018-01-10

Family

ID=55271023

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014151157A Active JP6259370B2 (en) 2014-07-24 2014-07-24 Optical coherence tomography device

Country Status (1)

Country Link
JP (1) JP6259370B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6539235B2 (en) * 2016-08-26 2019-07-03 日本電信電話株式会社 Signal processing apparatus and rescaling method
CN109807471B (en) * 2019-02-01 2024-03-26 佛山科学技术学院 Laser marking device and method

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2171396B1 (en) * 2007-07-12 2020-05-13 Volcano Corporation Apparatus and methods for uniform frequency sample clocking
WO2011037980A2 (en) * 2009-09-22 2011-03-31 Bioptigen, Inc. Systems for extended depth frequency domain optical coherence tomography (fdoct) and related methods
JP5939866B2 (en) * 2012-04-05 2016-06-22 キヤノン株式会社 Optical coherence tomography imaging apparatus and imaging method

Also Published As

Publication number Publication date
JP2016024176A (en) 2016-02-08

Similar Documents

Publication Publication Date Title
JP6166697B2 (en) Optical coherence tomography device
US11098997B2 (en) Interferometric distance measurement based on compression of chirped interferogram from cross-chirped interference
US20040239943A1 (en) System and method for low coherence broadband quadrature interferometry
JP5984693B2 (en) Optical coherence tomography apparatus and optical coherence tomography method
JP2016504589A (en) Optical coherence tomography system reconfigurable between different imaging modes
WO2012103557A2 (en) Spectral phase analysis for precision ranging
US9441947B2 (en) N-wavelength interrogation system and method for multiple wavelength interferometers
US9494410B2 (en) Method for measuring characteristics of sample
JPWO2018070442A1 (en) Optical angle modulation measuring apparatus and measuring method
RU2654379C1 (en) Instant optical coherent tomography in temporary area
JP7222433B2 (en) rangefinder
JP5724133B2 (en) Structure measuring method and structure measuring apparatus
JP6259370B2 (en) Optical coherence tomography device
Amin et al. Improved displacement sensing by spectral processing of laser self-mixing interferometry signal phase
JP6539235B2 (en) Signal processing apparatus and rescaling method
JP7061364B2 (en) Distance measuring device and distance measuring method
JP6467279B2 (en) Signal processing device
WO2021245778A1 (en) Ranging device
JP6047059B2 (en) Data processing apparatus and resampling method
JP2021032661A (en) Interferometer
JP6088338B2 (en) Data processing apparatus and resampling method
JP6168652B2 (en) Data processing apparatus and resampling method
JP7420140B2 (en) Wavelength swept optical measurement system
JP2019537726A (en) Method and apparatus for optimizing optical performance of an interferometer
JP7069993B2 (en) Optical spectrum line width calculation method, device and program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160913

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170530

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170531

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170719

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20171208

R150 Certificate of patent or registration of utility model

Ref document number: 6259370

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150