JP5616626B2 - Optical tomographic imaging apparatus and operating method thereof - Google Patents

Optical tomographic imaging apparatus and operating method thereof Download PDF

Info

Publication number
JP5616626B2
JP5616626B2 JP2009297876A JP2009297876A JP5616626B2 JP 5616626 B2 JP5616626 B2 JP 5616626B2 JP 2009297876 A JP2009297876 A JP 2009297876A JP 2009297876 A JP2009297876 A JP 2009297876A JP 5616626 B2 JP5616626 B2 JP 5616626B2
Authority
JP
Japan
Prior art keywords
light
optical
tomographic image
image data
optical tomographic
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
JP2009297876A
Other languages
Japanese (ja)
Other versions
JP2011137717A (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.)
Topcon Corp
Original Assignee
Topcon 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 Topcon Corp filed Critical Topcon Corp
Priority to JP2009297876A priority Critical patent/JP5616626B2/en
Publication of JP2011137717A publication Critical patent/JP2011137717A/en
Application granted granted Critical
Publication of JP5616626B2 publication Critical patent/JP5616626B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は光断層画像化装置及びその作動方法に係り、特にOCT(Optical Coherence Tomography)計測により光断層画像を生成する光断層画像化装置及びその作動方法に関する。 The present invention relates to an optical tomographic imaging apparatus and an operating method thereof, and more particularly, to an optical tomographic imaging apparatus that generates an optical tomographic image by OCT (Optical Coherence Tomography) measurement and an operating method thereof.

従来、生体組織の光断層画像を取得する際に、OCT計測を利用した光断層画像取得装置を用いることが提案されている。たとえば眼底や前眼部、皮膚の断層画像を取得する場合の他に、OCTプローブ(光プローブ)を用いる動脈血管壁の観察、内視鏡の鉗子チャンネルからOCTプローブを挿入する消化器管の観察など、様々な部位に応用されている。この光断層画像取得装置では、光源から射出された低コヒーレント光を測定光と参照光とに分割した後、該測定光が測定対象に照射されたときの測定対象からの反射光、もしくは後方散乱光と参照光とを合波し、該反射光と参照光との干渉光の強度に基づいて光断層画像を取得する。以下、測定対象からの反射光、後方散乱光をまとめて反射光と標記する。   Conventionally, it has been proposed to use an optical tomographic image acquisition apparatus using OCT measurement when acquiring an optical tomographic image of a living tissue. For example, in addition to acquiring tomographic images of the fundus, anterior segment, and skin, observation of an arterial blood vessel wall using an OCT probe (optical probe), observation of a digestive tract in which an OCT probe is inserted from a forceps channel of an endoscope It is applied to various parts. In this optical tomographic image acquisition apparatus, after the low-coherent light emitted from the light source is divided into measurement light and reference light, reflected light from the measurement object when the measurement light is irradiated to the measurement object, or backscattering The light and the reference light are combined, and an optical tomographic image is acquired based on the intensity of the interference light between the reflected light and the reference light. Hereinafter, the reflected light and the backscattered light from the measurement object are collectively referred to as reflected light.

上記のOCT計測には、大きくわけてTD−OCT(Time domain OCT)計測とFD−OCT(Fourier Domain OCT)計測の2種類がある。TD−OCT(Time domain OCT)計測は、参照光の光路長を変更しながら干渉光強度を測定することにより、測定対象の深さ方向の位置(以下、深さ位置という)に対応した反射光強度分布を取得する方法である。   The OCT measurement is roughly divided into two types: TD-OCT (Time domain OCT) measurement and FD-OCT (Fourier Domain OCT) measurement. In TD-OCT (Time domain OCT) measurement, reflected light corresponding to a position in the depth direction of a measurement target (hereinafter referred to as a depth position) is measured by measuring the intensity of interference light while changing the optical path length of the reference light. This is a method for obtaining an intensity distribution.

一方、FD−OCT計測は、参照光の光路長は変えることなく、光のスペクトル成分毎に干渉光強度を測定し、ここで得られたスペクトル干渉強度信号を計算機にてフーリエ変換に代表される周波数解析を行うことで、深さ位置に対応した反射光強度分布を取得する方法である。TD―OCTに存在する機械的な走査が不要となることで、高速な測定が可能となる手法として、近年注目されている。FD−OCT計測を行う装置構成で代表的なものとしては、SD−OCT(Spectral Domain OCT)装置とSS−OCT(Swept source OCT)の2種類が挙げられる。   On the other hand, in the FD-OCT measurement, the interference light intensity is measured for each spectral component of the light without changing the optical path length of the reference light, and the spectral interference intensity signal obtained here is represented by Fourier transform using a computer. This is a method of acquiring a reflected light intensity distribution corresponding to a depth position by performing frequency analysis. In recent years, it has attracted attention as a technique that enables high-speed measurement by eliminating the need for mechanical scanning existing in TD-OCT. As a typical apparatus configuration for performing FD-OCT measurement, there are two types, that is, an SD-OCT (Spectral Domain OCT) apparatus and an SS-OCT (Swept source OCT).

OCT計測技術は、特に細胞レベルでの測定による病変部位の特定を非侵襲で行う光バイオプシ(生検)への期待は高い。しかし、現状のOCTの分解能は約10μmであり、10〜20μm大の細胞レベルの観察には不十分であり、さらなる高分解能化が望まれている。   The OCT measurement technique has high expectations for optical biopsies (biopsy) that perform non-invasive specification of a lesion site by measurement at the cellular level. However, the resolution of the current OCT is about 10 μm, which is insufficient for observation of a cell level as large as 10 to 20 μm, and further higher resolution is desired.

OCTの軸分解能Δzは、以下の式(1)により使用する光源のスペクトル幅と中心波長により決定される。   The axial resolution Δz of OCT is determined by the spectral width and center wavelength of the light source used according to the following equation (1).

Δz=(2・ln2/π)・(λ2/Δλ) ……(1)
Δλ:スペクトル幅、 λ:中心波長
OCTの分解能を上げるために光源のスペクトル幅を広げることが必要であり、Ti:sapphire光源などの研究開発が進められている。しかしながら、現時点ではこれらの広帯域光源は高価であり、量産性も有していないという欠点がある。
Δz = (2 · ln2 / π) · (λ 2 / Δλ) (1)
Δλ: Spectral width, λ: Center wavelength It is necessary to widen the spectral width of the light source in order to increase the resolution of OCT, and research and development of Ti: sapphire light source and the like are underway. However, at present, these broadband light sources are expensive and do not have mass productivity.

この問題を解決するために、近年、複数の波長帯域を有するSLDを用いて、擬似的に広帯域光源を作り出すことによる高分解能化が検討されている。この擬似的な広帯域化での課題は、生体に一般的に用いられる中心波長1.0μmおよび1.3μmのSLD光源では全スペクトル帯域をカバーできないため、スペクトル帯域が分離してしまうことにある。スペクトル帯の分離は干渉信号を不連続に分断することとなり、OCT画像にアーティファクトを生成する原因となってしまう。特に、TD方式においては取得信号がOCT画像そのものとなってしまうため、これらのアーティファクト防ぐことは困難である。   In order to solve this problem, in recent years, high resolution has been studied by creating a pseudo broadband light source using an SLD having a plurality of wavelength bands. The problem with this pseudo-broadbanding is that the spectral bands are separated because the SLD light sources having center wavelengths of 1.0 μm and 1.3 μm that are generally used for living bodies cannot cover the entire spectral band. Spectral band separation discontinuously divides the interference signal, causing artifacts in the OCT image. In particular, in the TD system, the acquired signal becomes an OCT image itself, and it is difficult to prevent these artifacts.

そこで、FD方式を用いた合波OCT技術として、例えば、複数光源を用いた複数のOCTシステムを組み合わせる技術が開示されている(特許文献1)。   Thus, as a combined OCT technique using the FD method, for example, a technique of combining a plurality of OCT systems using a plurality of light sources is disclosed (Patent Document 1).

ところが、FD方式においては、検出光と参照光の光路差が干渉信号の周波数として表されるため、同一物体を測定したとしても、単純に2つ以上のシステムを連結させる構成の従来の合波OCT技術では、互いのシステムにおける参照面までの距離により、干渉信号の周波数が異なってしまう問題があり、高精度での参照面までの距離の位置あわせが必要となる。すなわち、OCT画像自体が10μm以下の分解能を有するため、その分解能以下での参照面までの距離の調整が必要となる。   However, in the FD method, since the optical path difference between the detection light and the reference light is expressed as the frequency of the interference signal, even if the same object is measured, the conventional multiplexing of a configuration in which two or more systems are simply connected. In the OCT technology, there is a problem that the frequency of the interference signal varies depending on the distance to the reference plane in each system, and the positioning of the distance to the reference plane with high accuracy is required. That is, since the OCT image itself has a resolution of 10 μm or less, it is necessary to adjust the distance to the reference surface with the resolution or less.

このような問題(高精度での参照面までの距離の位置あわせ)に対応するため、例えば低反射体を用いたキャリブレーション方法として、プローブを交換した場合の長さのばらつきを補正する方法が示されている(特許文献2)。   In order to cope with such a problem (positioning of the distance to the reference surface with high accuracy), for example, as a calibration method using a low reflector, there is a method of correcting variation in length when the probe is replaced. (Patent Document 2).

この特許文献2では、低反射体の信号からプローブの長さ位置によるオフセット量を読取り、参照面の位置の移動に使用しており、低反射体からの信号により、異なるプローブ長さ場合の参照面自動位置移動に適応できる。   In this patent document 2, the offset amount according to the length position of the probe is read from the signal of the low reflector and used to move the position of the reference surface. The reference when the probe length is different depending on the signal from the low reflector. It can be applied to automatic surface position movement.

特開2008−128708号公報JP 2008-128708 A 特開2006−215006号公報JP 2006-215006 A

しかしながら、複数光源を用いた複数のOCTシステムにおいて、上記特許文献2のキャリブレーション方法を適用する場合、実際に平均値などを用いて中間スペクトルの推定処理を行うと、2つの干渉信号間の位相の不一致により、アーティファクトの原因となる。   However, in the case of applying the calibration method of Patent Document 2 in a plurality of OCT systems using a plurality of light sources, if an intermediate spectrum estimation process is actually performed using an average value or the like, the phase between two interference signals This discrepancy causes artifacts.

すなわち、複数光源を用いた複数のOCTシステムでは、異なる中心波長を有する光源を使用するため、物理的な距離ではなく、光学的な距離を合わせる必要があり、ファイバ自体の曲げや環境などによる屈折率変化に対してもその分解能の低下を招く原因となってしまうため、リアルタイムでのキャリブレーションが必要なる。   That is, since a plurality of OCT systems using a plurality of light sources use light sources having different center wavelengths, it is necessary to match the optical distance, not the physical distance, and the refraction due to the bending of the fiber itself or the environment. Even if the rate changes, it causes a reduction in resolution, so calibration in real time is necessary.

本発明は、このような事情に鑑みてなされたもので、高精度な参照面の位置あわせが不要であり、複数の光源から得られる干渉信号を用いて分解能の低下がない高精細OCT画像を得ることのできる光断層画像化装置及びその作動方法を提供することを目的とする。 The present invention has been made in view of such circumstances, and does not require high-precision alignment of a reference surface, and uses a interference signal obtained from a plurality of light sources to produce a high-definition OCT image that does not have a reduction in resolution. It is an object of the present invention to provide an optical tomographic imaging apparatus and an operating method thereof.

前記目的を達成するために、本発明の第1態様に係る光断層画像化装置は、互いに波長帯域が異なり、それぞれ所定の前記波長帯域内で波長を掃引して複数の光束を射出する光源ユニットと、前記光源ユニットから射出された前記複数の各光束をそれぞれ測定光と参照光とに分割する光分割手段と、前記光分割手段により分割された複数の前記測定光を測定対象に照射する手段であって、前記測定光を導く光ファイバと、前記光ファイバの先端に固定され、前記光ファイバの先端から出射した前記測定光を前記測定対象に向けて偏向させる先端光学系とを有する測定光照射手段と、前記先端光学系と前記測定対象との間に配置され、前記測定光の一部を反射し他の一部を透過する基準反射手段と、前記測定光が測定対象に照射されたときの前記測定対象からの測定反射光及び前記基準反射手段からの基準反射光と、前記参照光とのそれぞれの干渉光信号を前記光束毎に検出する複数の干渉光検出手段と、前記干渉光検出手段にて検出された前記光束毎の前記干渉光信号をフーリエ変換して得られる光断層画像データのピーク位置を検出するピーク位置抽出手段と、前記複数の光束のうち一の光束についての光断層画像データのピーク位置と、他の光束についての光断層画像データのピーク位置との差をピーク情報として算出するピーク情報算出手段と、前記ピーク情報に基づいて前記他の光束についての光断層画像データを補正する光断層画像データ補正手段と、前記一の光束についての光断層画像データ、及び前記光断層画像データ補正手段による補正後の光断層画像データを逆フーリエ変換してそれぞれの干渉光信号を前記光束毎に再構築する干渉光再構築手段と、を備える。 In order to achieve the above object, the optical tomographic imaging apparatus according to the first aspect of the present invention is a light source unit that emits a plurality of light fluxes, each having a different wavelength band and sweeping a wavelength within the predetermined wavelength band. And a light dividing means for dividing the plurality of light beams emitted from the light source unit into measurement light and reference light, respectively, and means for irradiating the measurement object with the plurality of measurement lights divided by the light dividing means An optical fiber that guides the measurement light, and a measurement light that is fixed to the distal end of the optical fiber and that deflects the measurement light emitted from the distal end of the optical fiber toward the measurement object. irradiating means, wherein disposed between the distal optical system and the measurement object, and the reference reflecting means for over-permeability of another part reflects part of the measurement light, the measurement light is irradiated to the measurement object Before A reference light reflected from the measurement reflected light and the reference reflecting means from the measurement object, the respective interference light signal with the reference light and a plurality of interference light detecting means for detecting for each of the light beam, the interference light detecting means Peak position extraction means for detecting a peak position of optical tomographic image data obtained by Fourier transforming the interference light signal for each of the light beams detected in this manner, and optical tomographic image data for one of the plurality of light beams The peak information calculation means for calculating the difference between the peak position of the light beam and the peak position of the optical tomographic image data for the other light flux as peak information, and the optical tomographic image data for the other light flux is corrected based on the peak information The optical tomographic image data correcting means, the optical tomographic image data for the one light beam, and the optical tomographic image data corrected by the optical tomographic image data correcting means are reversed. It comprises a interference light reconstructing means for reconstructing the respective interference light signal for each of the light beams, a and Rie conversion.

本発明の第1態様によれば、前記複数の干渉光検出手段が前記測定光が測定対象に照射されたときの前記測定対象からの測定反射光及び前記基準反射手段からの基準反射光と、前記参照光とのそれぞれの干渉光信号を前記光束毎に検出することで、高精度な参照面の位置あわせが不要であり、複数の光源から得られる干渉信号を用いて分解能の低下がない高精細OCT画像を得ることを可能とする。 According to the first aspect of the present invention, when the plurality of interference light detection means is irradiated with the measurement light, the measurement reflected light from the measurement object and the reference reflected light from the reference reflection means, By detecting each interference light signal with the reference light for each light beam, it is not necessary to align the reference surface with high accuracy, and there is no reduction in resolution using interference signals obtained from a plurality of light sources. It is possible to obtain a fine OCT image.

本発明の第態様に係る光断層画像化装置は、第態様において、前記ピーク位置抽出手段は、前記測定光照射手段による前記測定光の照射方向に関して、前記光断層画像データが所定閾値を越える位置であって前記先端光学系に最も近い位置を前記ピーク位置として検出する。 The optical tomographic imaging apparatus according to a second aspect of the present invention is the optical tomographic imaging apparatus according to the first aspect, wherein the peak position extraction means has a predetermined threshold value for the optical tomographic image data with respect to the measurement light irradiation direction by the measurement light irradiation means. The position that is beyond and closest to the tip optical system is detected as the peak position.

本発明の第態様に係る光断層画像化装置は、第態様又は第態様において、前記干渉光再構築手段によって前記光束毎にそれぞれ再構築された前記干渉光信号を合成した合成干渉光信号を生成する合成干渉光信号生成手段と、前記合成干渉光信号より前記測定対象の光断層像情報を生成する光断層像情報生成手段と、をさらに備える。 The optical tomographic imaging apparatus according to a third aspect of the present invention is the combined interference light obtained by synthesizing the interference light signal reconstructed for each of the light beams by the interference light reconstruction means in the first aspect or the second aspect. And a combined interference light signal generating means for generating a signal, and an optical tomographic image information generating means for generating optical tomographic image information of the measurement object from the combined interference light signal.

本発明の第態様に係る光断層画像化装置は、第1態様ないし第態様のいずれかにおいて、前記複数の光束の波長掃引周期は、同期されている。 The optical tomographic imaging apparatus according to the fourth aspect of the present invention is the optical tomography apparatus according to any one of the first to third aspects, wherein the wavelength sweep periods of the plurality of light beams are synchronized.

本発明の第態様に係る光断層画像化装置の作動方法は、互いに波長帯域が異なり、それぞれ所定の前記波長帯域内で波長を掃引して複数の光束を射出する光源ユニットと、前記光源ユニットから射出された前記複数の各光束をそれぞれ測定光と参照光とに分割する光分割手段と、前記光分割手段により分割された複数の前記測定光を測定対象に照射する手段であって、前記測定光を導く光ファイバと、前記光ファイバの先端に固定され、前記光ファイバの先端から出射した前記測定光を前記測定対象に向けて偏向させる先端光学系とを有する測定光照射手段と、前記先端光学系と前記測定対象との間に配置され、前記測定光の一部を反射し他の一部を透過する基準反射手段と、複数の干渉光検出手段と、ピーク位置抽出手段と、ピーク情報算出手段と、光断層画像データ補正手段と、干渉光再構築手段と、を備える光断層画像化装置の作動方法であって、前記複数の干渉光検出手段が、前記測定光が測定対象に照射されたときの前記測定対象からの測定反射光及び前記基準反射手段からの基準反射光と、前記参照光とのそれぞれの干渉光信号を前記光束毎に検出前記ピーク位置抽出手段が、前記光束毎の前記干渉光信号をフーリエ変換して得られる光断層画像データのピーク位置を検出前記ピーク情報算出手段が、前記複数の光束のうち一の光束についての光断層画像データのピーク位置と、他の光束についての光断層画像データのピーク位置との差をピーク情報として算出前記光断層画像データ補正手段が、前記ピーク情報に基づいて前記他の光束についての光断層画像データを補正前記干渉光再構築手段が、前記一の光束についての光断層画像データ、及び前記補正後の光断層画像データを逆フーリエ変換してそれぞれの干渉光信号を前記光束毎に再構築する An optical tomographic imaging device operating method according to a fifth aspect of the present invention includes: a light source unit that emits a plurality of light fluxes each having a different wavelength band and sweeping a wavelength within a predetermined wavelength band; and the light source unit A light dividing means for dividing each of the plurality of light beams emitted from the light into measurement light and reference light, and means for irradiating the measurement object with the plurality of measurement lights divided by the light dividing means, Measuring light irradiation means comprising: an optical fiber that guides measurement light; and a tip optical system that is fixed to the tip of the optical fiber and deflects the measurement light emitted from the tip of the optical fiber toward the measurement object; and A reference reflecting means that is disposed between the tip optical system and the measurement object, reflects a part of the measurement light and transmits the other part, a plurality of interference light detection means, a peak position extraction means, a peak information Means leaving, an optical tomographic image data correcting means, a method of operating an optical tomography system including the interference light reconstructing means, wherein the plurality of interference light detecting means, the irradiation the measurement light to the measurement target When the measured reflected light from the measurement object and the reference reflected light from the reference reflecting means, and the interference light signal of the reference light is detected for each of the luminous flux, the peak position extracting means, A peak position of optical tomographic image data obtained by Fourier transforming the interference light signal for each light beam is detected, and the peak information calculation means is a peak position of optical tomographic image data for one light beam among the plurality of light beams. When the difference between the peak position of the optical tomographic image data for the other light beam is calculated as the peak information, the optical tomographic image data correcting means, optical disconnection of the other light beam on the basis of the peak information Correcting the image data, the interference light reconstructing means, optical tomographic image data, and the light flux of each of the coherent light signal an optical tomographic image data by performing inverse Fourier transform after pre Kiho positive for the one light beam Rebuild every time .

以上説明したように、本発明によれば、複数の光源から得られる干渉信号を用いてアーティファクトの発生がない高精細なOCT画像を得ることができるという効果がある。   As described above, according to the present invention, there is an effect that a high-definition OCT image free from artifacts can be obtained using interference signals obtained from a plurality of light sources.

光断層画像化装置の好ましい実施形態を示す模式図Schematic diagram showing a preferred embodiment of an optical tomographic imaging apparatus 図1の光プローブの要部の構成を示す図The figure which shows the structure of the principal part of the optical probe of FIG. 図1の演算部の構成を示すブロック図The block diagram which shows the structure of the calculating part of FIG. 2つのシステムの異なる参照面による干渉信号を示す図Diagram showing interference signals due to different reference planes of the two systems 図4の2つのシステムの干渉信号のキャリブレーションを説明する図The figure explaining the calibration of the interference signal of two systems of FIG. 図3の演算部の処理の流れを説明するフローチャートでFIG. 4 is a flowchart for explaining a processing flow of the arithmetic unit in FIG. 3. 図6の処理における反射光L3aにおけるピーク位置検出を説明する図The figure explaining the peak position detection in the reflected light L3a in the process of FIG. 図6の処理における反射光L3bにおけるピーク位置検出を説明する図The figure explaining the peak position detection in the reflected light L3b in the process of FIG. 図1の光断層画像化装置のOCTプローブが適用可能な内視鏡装置と併用した画像診断装置を示す図The figure which shows the diagnostic imaging apparatus used together with the endoscope apparatus which can apply the OCT probe of the optical tomographic imaging apparatus of FIG.

以下、添付図面を参照して、本発明に係る光断層画像化装置及びその作動方法について詳細に説明する。 Hereinafter, an optical tomographic imaging apparatus and an operation method thereof according to the present invention will be described in detail with reference to the accompanying drawings.

図1は光断層画像化装置の好ましい実施形態を示す模式図である。また、図2は図1の光プローブの要部の構成を示す図である。   FIG. 1 is a schematic diagram showing a preferred embodiment of an optical tomographic imaging apparatus. FIG. 2 is a diagram showing a configuration of a main part of the optical probe of FIG.

本実施形態の光断層画像化装置1は、例えば体腔内の生体組織や細胞等の測定対象の断層画像をSS−OCT(Swept source OCT)計測により取得するものであって、図1に示すように、光源ユニット2、光分割手段としての光分割部3A、3B、合波部6A、6B、干渉光検出手段としての干渉光検出部8A、8B、ピーク情報検出手段、差分算出手段及び干渉情報補正手段としての演算部20、表示装置11を有して構成される。なお、演算部20は、例えばパーソナルコンピュータ等により構成される。   The optical tomographic imaging apparatus 1 according to the present embodiment acquires a tomographic image of a measurement target such as a living tissue or a cell in a body cavity by SS-OCT (Swept source OCT) measurement, as shown in FIG. In addition, the light source unit 2, the light splitting units 3A and 3B as the light splitting means, the multiplexing units 6A and 6B, the interference light detection units 8A and 8B as the interference light detection means, the peak information detection means, the difference calculation means, and the interference information It has the calculating part 20 as a correction means, and the display apparatus 11. In addition, the calculating part 20 is comprised by a personal computer etc., for example.

光源ユニット2は、互いに分離した波長帯域Δλa、Δλbを有する複数の光束La、Lbを射出するものである。具体的には、光源ユニット2は波長帯域Δλa内において波長を一定の周期で掃引させながらレーザ光を射出する第1光源2Aと、波長帯域Δλb内において波長を一定の周期で掃引させながらレーザ光を射出する第2光源2Bを有している。したがって、各光源2A、2Bは波長を1周期掃引させたときにそれぞれ波長帯域Δλa、Δλbからなる光束La、Lbを射出したことになる。   The light source unit 2 emits a plurality of light beams La and Lb having wavelength bands Δλa and Δλb separated from each other. Specifically, the light source unit 2 includes a first light source 2A that emits laser light while sweeping the wavelength in the wavelength band Δλa, and a laser beam that sweeps the wavelength in the wavelength band Δλb. Has a second light source 2B. Accordingly, the light sources 2A and 2B emit light beams La and Lb having wavelength bands Δλa and Δλb, respectively, when the wavelength is swept for one period.

光分割部3A、3Bは、例えば、分岐比90:10の2×2の光カプラから構成されている。光分割部3Aは、光束Laを測定光L1aと参照光L2aとに分割し、光分割部3Bは、光束Lbを測定光L1bと参照光L2bとに分割する。このとき、光分割部3A、3Bは、測定光:参照光=90:10の割合で分割する。ここで、測定光L1a、L1bは、サーキュレータ4A、4B及び合分波部12を介して光プローブ600に導波される。   The light splitting units 3A and 3B are composed of, for example, a 2 × 2 optical coupler having a branching ratio of 90:10. The light splitting unit 3A splits the light beam La into measurement light L1a and reference light L2a, and the light splitting unit 3B splits the light beam Lb into measurement light L1b and reference light L2b. At this time, the light dividing units 3A and 3B divide at a ratio of measurement light: reference light = 90: 10. Here, the measurement lights L1a and L1b are guided to the optical probe 600 via the circulators 4A and 4B and the multiplexing / demultiplexing unit 12.

光プローブ600は、光ロータリコネクタ18を介して入射された測定光L1a、L1bを測定対象Sまで導波し、測定対象Sの同一部位に同時に照射する。また、光プローブ600は、測定光L1a、L1bが測定対象Sに照射されたときの測定対象Sからの反射光L3a、L3bを導波する。   The optical probe 600 guides the measurement lights L1a and L1b incident through the optical rotary connector 18 to the measurement target S and irradiates the same part of the measurement target S simultaneously. Further, the optical probe 600 guides the reflected lights L3a and L3b from the measurement object S when the measurement light L1a and L1b are irradiated on the measurement object S.

光プローブ600は、図2に示すように、モータ(不図示)により光ロータリコネクタ18に連結され、透明なシース601に覆われた先端にボールレンズ680を設けたファイバ部FB1が回転する構成となっている。それにより光プローブ600は、シース601を介して測定対象S上において円周状に光束をラジアル走査(図2の矢印R)し、2次元断層画像が計測可能となっている。さらに、光プローブ600は、モータ(不図示)により光プローブ600の先端が光路の走査円が形成する平面に対して垂直な方向に進退走査(図2の矢印D1及びD2)することにより、3次元断層画像の計測も可能となっている。   As shown in FIG. 2, the optical probe 600 is connected to the optical rotary connector 18 by a motor (not shown), and has a configuration in which a fiber part FB1 having a ball lens 680 provided at the tip covered with a transparent sheath 601 rotates. It has become. As a result, the optical probe 600 can scan the light beam radially on the measurement target S via the sheath 601 (arrow R in FIG. 2) and measure a two-dimensional tomographic image. Furthermore, the optical probe 600 is moved forward and backward (arrows D1 and D2 in FIG. 2) in a direction perpendicular to the plane formed by the scanning circle of the optical path by a motor (not shown). Measurement of dimensional tomographic images is also possible.

また、光プローブ600は、ボールレンズ680からの光束の一部を反射する反射部材700がシース601の外周面あるいは内周面の、全面または一部に設けられている。   Further, in the optical probe 600, a reflection member 700 that reflects a part of the light beam from the ball lens 680 is provided on the entire outer surface or the inner peripheral surface of the sheath 601.

なお、反射部材700は、測定光の波長に対して所定の反射率を有するガラスなどの低反射部材により構成される。   The reflecting member 700 is made of a low reflecting member such as glass having a predetermined reflectance with respect to the wavelength of the measurement light.

図1に戻り、光断層画像化装置1は、光分割部3Aと光プローブ600の間の光路、光分割部3Bと光プローブ600の間の光路には合分波部12を設けている。合分波部12は、設定されたカットオフ波長に応じて光を合分波する機能を有し、たとえばWDM(Wavelength Division Multiplexing:波長分割多重)カプラにより構成される。   Returning to FIG. 1, the optical tomographic imaging apparatus 1 includes the multiplexing / demultiplexing unit 12 in the optical path between the optical dividing unit 3 </ b> A and the optical probe 600 and in the optical path between the optical dividing unit 3 </ b> B and the optical probe 600. The multiplexing / demultiplexing unit 12 has a function of multiplexing and demultiplexing light according to a set cutoff wavelength, and is configured by, for example, a WDM (Wavelength Division Multiplexing) coupler.

合分波部12は、光分割部3A、3B側からそれぞれ入射された測定光L1a、L1bを合波して光プローブ600側に射出し、光プローブ600側から入射された反射光L3a、L3bの光を分波してそれぞれ合波部6A、6B側へ射出する。   The multiplexing / demultiplexing unit 12 multiplexes the measurement lights L1a and L1b incident from the light splitting units 3A and 3B, respectively, and emits them to the optical probe 600 side, and the reflected lights L3a and L3b incident from the optical probe 600 side. Are demultiplexed and emitted to the multiplexing sections 6A and 6B, respectively.

反射光L3aは合波部6Aにおいて参照光L2aと合波され、反射光L3bは合波部6Bにおいて参照光L2bと合波される。   The reflected light L3a is combined with the reference light L2a in the combining unit 6A, and the reflected light L3b is combined with the reference light L2b in the combining unit 6B.

なお、光分割部3Aから合波部6Aまでの参照光L2aのサーキュレータ5Aを介した光路には例えば反射ミラー(不図示)等の参照面を有する反射型の光路長調整部7Aが設けられ、光分割部3Bから合波部6Bまでの参照光L2bのサーキュレータ5Bを介した光路には例えば反射ミラー(不図示)等の参照面を有する反射型の光路長調整部7Bが設けられている。光路長調整部7A、7Bは、断層画像の取得を開始する位置を調整するために、それぞれ参照光L2a、L2bの光路長を変更するものである。   A reflection type optical path length adjustment unit 7A having a reference surface such as a reflection mirror (not shown) is provided in the optical path of the reference light L2a from the light dividing unit 3A to the multiplexing unit 6A via the circulator 5A. A reflection type optical path length adjustment unit 7B having a reference surface such as a reflection mirror (not shown) is provided in the optical path of the reference light L2b from the light splitting unit 3B to the multiplexing unit 6B via the circulator 5B. The optical path length adjusting units 7A and 7B change the optical path lengths of the reference lights L2a and L2b, respectively, in order to adjust the position where the tomographic image acquisition is started.

合波部6A、6Bは、例えば、分岐比50:50の2×2の光ファイバカプラから構成されている。合波部6Aは反射光L3aと参照光L2aを合波して、このとき生じた干渉光L4aを干渉光検出部8Aへ射出する。合波部6Bは反射光L3bと参照光L2bを合波して、このとき生じた干渉光L4bを干渉光検出部8Bへ射出する。   The multiplexing units 6A and 6B are composed of, for example, a 2 × 2 optical fiber coupler with a branching ratio of 50:50. The multiplexing unit 6A combines the reflected light L3a and the reference light L2a, and emits the interference light L4a generated at this time to the interference light detection unit 8A. The multiplexing unit 6B combines the reflected light L3b and the reference light L2b, and emits the interference light L4b generated at this time to the interference light detection unit 8B.

なお、本光断層画像化装置1においては、合波部6A、6Bはそれぞれ干渉光L4a、L4bを二分して干渉光検出部8A、8Bへ射出し、干渉光検出部8A、8Bは二分された干渉光L4a、L4bをそれぞれ2つの光検出素子を用いてバランス検波する。このバランス検波機構により、光強度ゆらぎの影響を抑え、より鮮明な画像を得ることができる。   In the optical tomographic imaging apparatus 1, the multiplexing units 6A and 6B bisect the interference lights L4a and L4b, respectively, and emit them to the interference light detection units 8A and 8B, and the interference light detection units 8A and 8B are bisected. The detected interference lights L4a and L4b are each subjected to balance detection using two photodetectors. By this balance detection mechanism, the influence of light intensity fluctuation can be suppressed and a clearer image can be obtained.

干渉光検出部8A、8Bは各干渉光L4a、L4bをそれぞれ光電変換し、各光束La、Lbの波長帯域Δλa、Δλbごとの複数の干渉信号ISa、ISbとして検出する機能を有している。ここで、干渉光検出部8A、8Bは、光源2A、2Bからの波長掃引のトリガ信号9A、9Bにより同期をとることで、対応する光束を認識するようになっている。このとき、干渉光検出部8A、8Bにおいて、各光束毎La、Lbの干渉信号ISa、ISbが観測されることになる。干渉信号ISa、ISbは、演算部20に出力される。   The interference light detectors 8A and 8B have a function of photoelectrically converting the interference lights L4a and L4b, respectively, and detecting them as a plurality of interference signals ISa and ISb for the wavelength bands Δλa and Δλb of the light beams La and Lb. Here, the interference light detection units 8A and 8B recognize the corresponding light flux by synchronizing with the wavelength sweep trigger signals 9A and 9B from the light sources 2A and 2B. At this time, the interference light detection units 8A and 8B observe the interference signals ISa and ISb of each light beam La and Lb. The interference signals ISa and ISb are output to the calculation unit 20.

そして、演算部20にて干渉信号ISa、ISbが信号処理され、断層画像が表示装置11に表示される。   The arithmetic unit 20 processes the interference signals ISa and ISb, and displays a tomographic image on the display device 11.

図3は図1の演算部の構成を示すブロック図である。光断層画像化装置1の演算部20は、図3に示すように、第1干渉波データ格納部21、第2干渉波データ格納部22、第1光断層画像データ生成部23、第2光断層画像データ生成部24、ピーク情報検出手段及び差分算出手段としてのピーク位置抽出部25、干渉情報補正手段としての第2光断層画像データ補正部26、第1干渉波データ再構築部27、第2干渉波データ再構築部28、合成干渉光信号生成手段としての合成干渉波生成部29及光断層像情報生成手段としてのFFT部30を備えて構成される。   FIG. 3 is a block diagram showing the configuration of the calculation unit of FIG. As shown in FIG. 3, the calculation unit 20 of the optical tomographic imaging apparatus 1 includes a first interference wave data storage unit 21, a second interference wave data storage unit 22, a first optical tomographic image data generation unit 23, and a second light. A tomographic image data generation unit 24, a peak position extraction unit 25 as a peak information detection unit and a difference calculation unit, a second optical tomographic image data correction unit 26 as an interference information correction unit, a first interference wave data reconstruction unit 27, a first 2 includes an interference wave data reconstruction unit 28, a combined interference wave generation unit 29 as a combined interference light signal generation unit, and an FFT unit 30 as an optical tomographic image information generation unit.

第1干渉波データ格納部21は、光束La側の干渉信号ISaのデジタルデータを格納するものであり、第2干渉波データ格納部22は、光束Lb側の干渉信号ISbのデジタルデータを格納するものである。   The first interference wave data storage unit 21 stores digital data of the interference signal ISa on the light beam La side, and the second interference wave data storage unit 22 stores digital data of the interference signal ISb on the light beam Lb side. Is.

第1光断層画像データ生成部23は、第1干渉波データ格納部21に格納された干渉信号ISaのデジタルデータをFFT(高速フーリエ変換)して、第1光断層画像データを生成するものである。   The first optical tomographic image data generation unit 23 generates first optical tomographic image data by performing FFT (Fast Fourier Transform) on the digital data of the interference signal ISa stored in the first interference wave data storage unit 21. is there.

第2光断層画像データ生成部24は、第2干渉波データ格納部22に格納された干渉信号ISbのデジタルデータをFFTして、第2光断層画像データを生成するものである。   The second optical tomographic image data generation unit 24 performs FFT on the digital data of the interference signal ISb stored in the second interference wave data storage unit 22 to generate second optical tomographic image data.

ピーク位置抽出部25は、第1光断層画像データ生成部23にて得られた第1光断層画像データ及び第2光断層画像データ生成部24にて得られた第2光断層画像データのそれぞれのデータより、反射部材700を示すピークの光断層画像上のピーク位置を検出し、第1光断層画像データのピーク位置Aと第2光断層画像データBのピーク位置との差(A−B)をピーク情報として算出するものである。   The peak position extraction unit 25 is configured to receive the first optical tomographic image data obtained by the first optical tomographic image data generation unit 23 and the second optical tomographic image data obtained by the second optical tomographic image data generation unit 24, respectively. From the data, the peak position on the optical tomographic image of the peak indicating the reflecting member 700 is detected, and the difference (AB) between the peak position A of the first optical tomographic image data and the peak position of the second optical tomographic image data B is detected. ) As peak information.

第2光断層画像データ補正部26は、ピーク情報(A−B)に基づいて第2光断層画像データ生成部24が生成した第2光断層画像データを補正するものである。   The second optical tomographic image data correction unit 26 corrects the second optical tomographic image data generated by the second optical tomographic image data generation unit 24 based on the peak information (AB).

なお、ピーク位置抽出部25及び第2光断層画像データ補正部26についての詳細は後述する。   Details of the peak position extraction unit 25 and the second optical tomographic image data correction unit 26 will be described later.

第1干渉波データ再構築部27は、第光断層画像データをIFFT(逆高速フーリエ変換)して、干渉信号ISaを再構築するものである。以下、再構築した干渉信号をISa‘と記す。 The first interference wave data reconstruction unit 27 reconstructs the interference signal ISa by performing IFFT (inverse fast Fourier transform) on the first optical tomographic image data. Hereinafter, the reconstructed interference signal is denoted as ISa ′.

第2干渉波データ再構築部28は、第2光断層画像データ補正部26にて補正した第2光断層画像データをIFFTして、干渉信号ISbを再構築するものである。以下、再構築した干渉信号をISb‘と記す。   The second interference wave data reconstruction unit 28 performs IFFT on the second optical tomographic image data corrected by the second optical tomographic image data correction unit 26 to reconstruct the interference signal ISb. Hereinafter, the reconstructed interference signal is referred to as ISb ′.

合成干渉波生成部29は、干渉信号ISa‘のデジタルデータと干渉信号ISb‘のデジタルデータとを波長帯域毎に合成して合成干渉波(合成干渉信号)を生成するものである。   The combined interference wave generation unit 29 combines the digital data of the interference signal ISa ′ and the digital data of the interference signal ISb ′ for each wavelength band to generate a combined interference wave (synthesized interference signal).

FFT部30は、合成干渉波(合成干渉信号)をFFTして測定対象Sの断層像情報(OCT画像)を生成し、表示装置11にOCT画像を表示させるものである。   The FFT unit 30 generates a tomographic image information (OCT image) of the measurement target S by performing FFT on the synthetic interference wave (synthetic interference signal), and causes the display device 11 to display the OCT image.

次に、このように構成された本実施形態の光断層画像化装置1における演算部20について説明する。   Next, the calculation unit 20 in the optical tomographic imaging apparatus 1 of the present embodiment configured as described above will be described.

図4は2つのシステムの異なる参照面による干渉信号を示す図であり、図5は図4の2つのシステムの干渉信号のキャリブレーションを説明する図である。   FIG. 4 is a diagram illustrating interference signals due to different reference surfaces of the two systems, and FIG. 5 is a diagram illustrating calibration of interference signals of the two systems in FIG.

干渉光検出部で得られる信号は、測定光及び反射光と、光路長調整部7の参照面にて反射された参照光との光学的距離差に依存した周波数をもつ、以下の式(1)に示すような干渉信号IDを生じる。 The signal obtained by the interference light detection unit has a frequency that depends on the optical distance difference between the measurement light and the reflected light and the reference light reflected by the reference surface of the optical path length adjustment unit 7 as follows: causing interference signal I D as shown in).

Figure 0005616626
Figure 0005616626

R:参照面位置、ZSn:測定対象位置、RR:参照面反射率、RSn:測定対象反射率
本実施形態のように、2つのシステムで同じ対象を測定した場合を考えたとき、それぞれのシステムから得られる干渉信号はそれぞれの参照面の位置に依存する。すなわち、干渉光検出部8A、8Bで得られる信号は、それぞれの測定光L1a,L1及び反射光L3a、L3bと、光路長調整部7A、7Bの参照面にて反射された参照光L2a,L2bとの光学的距離差に依存した周波数をもつ、以下の式(2)、(3)に示すような干渉信号ID1,ID2を生じる。
Z R : Reference surface position, Z Sn : Measurement object position, R R : Reference surface reflectance, R Sn : Measurement object reflectance When the same object is measured by two systems as in this embodiment The interference signal obtained from each system depends on the position of the respective reference plane. That is, the signals obtained by the interference light detection units 8A and 8B are the measurement light L1a and L1 and the reflected light L3a and L3b, and the reference light L2a and L2b reflected by the reference surfaces of the optical path length adjustment units 7A and 7B. Interference signals I D1 and I D2 as shown in the following formulas (2) and (3) having a frequency that depends on the difference in optical distance from the above are generated.

Figure 0005616626
Figure 0005616626

Figure 0005616626
Figure 0005616626

R1:光路長調整部7Aの参照面位置、ZR2:光路長調整部7Bの参照面位置、RR1:光路長調整部7Aの参照面反射率、RR2:光路長調整部7Bの参照面反射率
図4に示すように、結果的に参照面の位置の差ΔLにより、異なる周波数の干渉信号が得られることとなる。
Z R1 : Reference plane position of the optical path length adjustment unit 7A, Z R2 : Reference plane position of the optical path length adjustment unit 7B, R R1 : Reference plane reflectance of the optical path length adjustment unit 7A, R R2 : Reference of the optical path length adjustment unit 7B Surface Reflectance As shown in FIG. 4, as a result, interference signals having different frequencies are obtained due to the difference ΔL in the position of the reference surface.

つまり、これら2つの干渉信号ID1,ID2から高分解画像を作成しようとしたとき、同じもの見ているにも関わらず違う位置に信号がくるため、違う周波数の波と表現されることとなる。 In other words, when trying to create a high resolution image from these two interference signals I D1 and I D2 , the signals come in different positions even though they are looking at the same thing, so they are expressed as waves of different frequencies. Become.

本実施形態では、測定対象を光プローブ600先端に設けられた反射部材700とした際の、反射部材700の位置に対して常に上述した干渉信号ID1,ID2が取得できる。 In this embodiment, when the measurement target is the reflection member 700 provided at the tip of the optical probe 600, the above-described interference signals I D1 and I D2 can always be obtained with respect to the position of the reflection member 700.

この干渉信号ID1,ID2をフーリエ変換することでZR1−ZSnおよびZR2−ZSnにピークを持つ光断層画像を得ることができる。 An optical tomographic image having peaks at Z R1 -Z Sn and Z R2 -Z Sn can be obtained by Fourier transforming the interference signals I D1 and I D2 .

それぞれのピーク位置の差は
R1−ZSn−(ZR2−ZSn)=ZR1−ZR2 …(4)
となり、システム間(光路長調整部7A、7B)の参照面の位置の差であることがわかる。
The difference between the peak positions is Z R1 −Z Sn − (Z R2 −Z Sn ) = Z R1 −Z R2 (4)
Thus, it can be seen that this is the difference in the position of the reference plane between the systems (optical path length adjusting units 7A and 7B).

すなわち、2つのシステムをシステムA,Bとしたときに、図5に示すように、反射部材700の干渉信号ID1,ID2をキャリブレーション信号(A),(B)とし、キャリブレーション信号(A),(B)のピーク位置を検出し、その位置の差だけOCT信号(B)をシフトさせOCT信号(A)に一致させることで、システム間(光路長調整部7A、7B)の参照面の位置があった干渉信号を得ることができる。 That is, when the two systems are systems A and B, as shown in FIG. 5, the interference signals I D1 and I D2 of the reflecting member 700 are used as calibration signals (A) and (B), and the calibration signal ( By detecting the peak positions of A) and (B), and shifting the OCT signal (B) by the difference between the positions to match the OCT signal (A), reference between systems (optical path length adjustment units 7A and 7B) is made. An interference signal having a surface position can be obtained.

なお、反射部材700のピーク位置の検出には、反射部材700を測定対象Sの手前に配置することで、常に0の測定光の出射方向位置(図2参照)に最も近いピーク位置を閾値検出することが望ましい。   For the detection of the peak position of the reflecting member 700, the reflecting member 700 is arranged in front of the measuring object S, so that the peak position that is always closest to the position of the zero measurement light in the emission direction (see FIG. 2) is detected as a threshold. It is desirable to do.

本実施形態の演算部20は、上記のように、常に同じ位置にいる信号の検出することにより、リアルタイムに参照面間の差のキャリブレーションを行う。   As described above, the calculation unit 20 of the present embodiment calibrates the difference between the reference planes in real time by detecting a signal that is always at the same position.

続いて、本実施形態における演算部20の具体的な処理の流れを説明する。図6は図3の演算部の処理の流れを説明するフローチャートである。図7及び図8は図6の処理における反射光L3a,L3bにおけるピーク位置検出を説明する図である。   Next, a specific processing flow of the calculation unit 20 in the present embodiment will be described. FIG. 6 is a flowchart for explaining the flow of processing of the arithmetic unit in FIG. 7 and 8 are diagrams for explaining peak position detection in the reflected lights L3a and L3b in the process of FIG.

図6に示すように、演算部20は、第1干渉波データ格納部21に干渉信号ISaのデジタルデータを格納し、第2干渉波データ格納部22に干渉信号ISbのデジタルデータを格納する(ステップS1)。   As shown in FIG. 6, the calculation unit 20 stores the digital data of the interference signal ISa in the first interference wave data storage unit 21 and stores the digital data of the interference signal ISb in the second interference wave data storage unit 22 ( Step S1).

そして、演算部20は、第1光断層画像データ生成部23にて、第1干渉波データ格納部21に格納された干渉信号ISaのデジタルデータをFFTして、第1光断層画像データを生成すると共に、第2光断層画像データ生成部24にて、第2干渉波データ格納部22に格納された干渉信号ISbのデジタルデータをFFTして、第2光断層画像データを生成する(ステップS2)。   Then, the computing unit 20 generates the first optical tomographic image data by performing FFT on the digital data of the interference signal ISa stored in the first interference wave data storage unit 21 in the first optical tomographic image data generating unit 23. At the same time, the second optical tomographic image data generation unit 24 performs FFT on the digital data of the interference signal ISb stored in the second interference wave data storage unit 22 to generate second optical tomographic image data (step S2). ).

続いて、演算部20は、ピーク位置抽出部25にて、図7及び図8に示すように、第1光断層画像データ生成部23にて得られた第1光断層画像データ及び第2光断層画像データ生成部24にて得られた第2光断層画像データのそれぞれのデータより、反射部材700を示すピークの光断層画像上のピーク位置を検出する(ステップS3)。   Subsequently, the calculation unit 20 causes the peak position extraction unit 25 to obtain the first optical tomographic image data and the second light obtained by the first optical tomographic image data generation unit 23, as shown in FIGS. From the data of the second optical tomographic image data obtained by the tomographic image data generation unit 24, the peak position on the optical tomographic image indicating the reflecting member 700 is detected (step S3).

具体的には、図7及び図8に示すように、ピーク位置抽出部25は、参照光L2a,L2bにおける測定光の出射方向位置(図2参照)に対して得られたそれぞれの光断層画像データを所定の閾値により、出射方向位置が0に最も近いピーク位置を反射部材700のピーク位置として検出する。   Specifically, as shown in FIGS. 7 and 8, the peak position extraction unit 25 obtains each optical tomographic image obtained with respect to the measurement light emission direction position (see FIG. 2) in the reference lights L <b> 2 a and L <b> 2 b. The peak position whose emission direction position is closest to 0 is detected as the peak position of the reflecting member 700 based on the predetermined threshold value.

なお、所定の閾値は、ノイズの影響を十分受けないレベルに設定する。また、反射部材700の信号レベルは測定対象S手前での反射レベルを見ているため、光源の光量変動などの要因以外では測定対象Sに依存することなく安定的に検出可能である。   The predetermined threshold is set to a level that is not sufficiently affected by noise. In addition, since the signal level of the reflection member 700 is seen from the reflection level before the measurement object S, it can be stably detected without depending on the measurement object S except for factors such as fluctuations in the light amount of the light source.

そして、演算部20は、ピーク位置抽出部25にて、第1光断層画像データのピーク位置Aと第2光断層画像データBのピーク位置との差(A−B)をピーク情報として算出する(ステップS4)。   Then, the calculation unit 20 calculates, in the peak position extraction unit 25, the difference (A−B) between the peak position A of the first optical tomographic image data and the peak position of the second optical tomographic image data B as peak information. (Step S4).

次に、演算部20は、第2光断層画像データ補正部26は、ピーク情報(A−B)に基づいて第2光断層画像データ生成部24が生成した第2光断層画像データの出射方向位置をピーク情報(A−B)だけシフトさせ、第1光断層画像データの出射方向位置にキャリブレーションして補正第2光断層画像データを生成する(ステップS5)。   Next, the calculation unit 20 outputs the second optical tomographic image data generated by the second optical tomographic image data generation unit 24 based on the peak information (AB). The position is shifted by the peak information (A-B) and calibrated to the position in the emission direction of the first optical tomographic image data to generate corrected second optical tomographic image data (step S5).

また、演算部20は、第1干渉波データ再構築部27にて第光断層画像データをIFFTして、干渉信号ISaを再構築し干渉信号をISa‘を生成すると共に、第2干渉波データ再構築部28にて第2光断層画像データ補正部26にて補正した第2光断層画像データをIFFTして、干渉信号ISbを再構築し干渉信号をISb‘を生成する。そして、演算部20は、合成干渉波生成部29にて、干渉信号ISa‘のデジタルデータと干渉信号ISb‘のデジタルデータとを波長帯域毎に合成して合成干渉波(合成干渉信号)を生成する(ステップS6)。 In addition, the calculation unit 20 performs IFFT on the first optical tomographic image data in the first interference wave data reconstruction unit 27, reconstructs the interference signal ISa, generates the interference signal ISa ′, and generates the second interference wave. The data reconstruction unit 28 performs IFFT on the second optical tomographic image data corrected by the second optical tomographic image data correction unit 26, reconstructs the interference signal ISb, and generates the interference signal ISb ′. The arithmetic unit 20 then combines the digital data of the interference signal ISa ′ and the digital data of the interference signal ISb ′ for each wavelength band at the combined interference wave generation unit 29 to generate a combined interference wave (synthesized interference signal). (Step S6).

次に、演算部20は、FFT部30において、合成干渉波生成部29にて生成された、合成干渉波を高速フーリエ変換(FFT)して測定対象Sの断層像情報(OCT画像)を生成し、表示装置11にOCT画像を表示させる(ステップS7)。   Next, in the FFT unit 30, the calculation unit 20 performs fast Fourier transform (FFT) on the combined interference wave generated by the combined interference wave generation unit 29 to generate tomographic image information (OCT image) of the measurement target S. Then, the OCT image is displayed on the display device 11 (step S7).

このように本実施形態では、(1)反射部材700を示すピークの光断層画像上のピーク位置を検出し、そして第1光断層画像データのピーク位置Aと第2光断層画像データBのピーク位置との差(A−B)をピーク情報として算出し、(2)その後ピーク情報(A−B)に基づいて第2光断層画像データ生成部24が生成した第2光断層画像データの出射方向位置をピーク情報(A−B)だけシフトさせ、第1光断層画像データの出射方向位置にキャリブレーションするので、高精度な参照面の位置あわせが不要であり、このようにして得られた合成干渉波により、複数の光源から得られる干渉信号を用いた分解能の低下がない高精細OCT画像を得ることが可能となる。   As described above, in this embodiment, (1) the peak position on the optical tomographic image indicating the reflecting member 700 is detected, and the peak position A of the first optical tomographic image data and the peak of the second optical tomographic image data B are detected. A difference (AB) from the position is calculated as peak information, and (2) emission of the second optical tomographic image data generated by the second optical tomographic image data generation unit 24 based on the peak information (AB) thereafter. Since the position of the direction is shifted by the peak information (AB) and calibrated to the position of the emission direction of the first optical tomographic image data, it is not necessary to align the reference plane with high accuracy, and thus obtained. With the combined interference wave, it is possible to obtain a high-definition OCT image without a reduction in resolution using interference signals obtained from a plurality of light sources.

また、光ファイバ、プローブなどは、曲げなどの応力により屈折率が変化してしまうため、常に反射部材700の信号でリアルタイムにキャリブレーションすることで、外乱や環境変化に強いシステムの構築が可能である。   In addition, since the refractive index of optical fibers, probes, etc. changes due to stress such as bending, it is possible to construct a system that is resistant to external disturbances and environmental changes by always calibrating in real time with the signal of the reflecting member 700. is there.

なお、合成干渉波を用いて、見かけの周波数帯域を広げることでOCT計測の分解能を上げても良いし、光断層画像の差分により分光光断層画像に使用してもよい。   Note that the resolution of the OCT measurement may be increased by expanding the apparent frequency band using the synthetic interference wave, or may be used for the spectroscopic tomographic image by the difference of the optical tomographic image.

なお、本実施形態は、内視鏡装置と併用した画像診断装置に適用することが可能である。詳細に説明すると、図9に示すように、本実施形態の光プローブ600と内視鏡装置と併用した画像診断装置50は、主として内視鏡100、内視鏡プロセッサ200、光源装置300、生体断層画像生成装置としてのOCTプロセッサ400、及び表示手段としてのモニタ装置である画像表示部500とから構成されている。尚、内視鏡プロセッサ200は、光源装置300を内蔵するように構成されていてもよい。   Note that this embodiment can be applied to an image diagnostic apparatus used in combination with an endoscope apparatus. More specifically, as shown in FIG. 9, the diagnostic imaging apparatus 50 used in combination with the optical probe 600 and the endoscope apparatus of the present embodiment mainly includes an endoscope 100, an endoscope processor 200, a light source device 300, a living body. An OCT processor 400 serving as a tomographic image generation apparatus and an image display unit 500 serving as a monitor device serving as display means are included. The endoscope processor 200 may be configured to incorporate the light source device 300.

内視鏡100は、手元操作部112と、この手元操作部112に連設される挿入部114とを備える。術者は手元操作部112を把持して操作し、挿入部114を被検者の体内に挿入することによって観察を行う。   The endoscope 100 includes a hand operation unit 112 and an insertion unit 114 that is connected to the hand operation unit 112. The surgeon grasps and operates the hand operation unit 112 and performs observation by inserting the insertion unit 114 into the body of the subject.

手元操作部112には、鉗子挿入部138が設けられており、この鉗子挿入部138が挿入部114内に設けられている鉗子チャンネル(不図示)を介して先端部144の鉗子口156に連通されている。画像診断装置10では、OCTプローブである光プローブ600を鉗子挿入部138から挿入することによって、光プローブ600を鉗子口156から導出する。光プローブ600は、鉗子挿入部138から挿入され、鉗子口156から導出される挿入部602と、術者が光プローブ600を操作するための操作部604、及びコネクタ401を介して光プロセッサ400と接続されるケーブル606から構成されている。   The hand operation unit 112 is provided with a forceps insertion portion 138, and this forceps insertion portion 138 communicates with a forceps port 156 of the distal end portion 144 via a forceps channel (not shown) provided in the insertion portion 114. Has been. In the diagnostic imaging apparatus 10, the optical probe 600 is led out from the forceps opening 156 by inserting the optical probe 600 that is an OCT probe from the forceps insertion portion 138. The optical probe 600 is inserted from the forceps insertion portion 138 and inserted from the forceps port 156, the operation portion 604 for the operator to operate the optical probe 600, and the optical processor 400 via the connector 401. It consists of a cable 606 to be connected.

内視鏡100の先端部144には、観察光学系150、照明光学系152、及びCCD(不図示)が配設されている。   At the distal end portion 144 of the endoscope 100, an observation optical system 150, an illumination optical system 152, and a CCD (not shown) are disposed.

観察光学系150は、被検体を図示しないCCDの受光面に結像させ、CCDは受光面上に結像された被検体像を各受光素子によって電気信号に変換する。この実施の形態のCCDは、3原色の赤(R)、緑(G)、青(B)のカラーフィルタが所定の配列(ベイヤー配列、ハニカム配列)で各画素ごとに配設されたカラーCCDである。   The observation optical system 150 forms an image of a subject on a light receiving surface (not shown) of the CCD, and the CCD converts the subject image formed on the light receiving surface into an electric signal by each light receiving element. The CCD of this embodiment is a color CCD in which three primary color red (R), green (G), and blue (B) color filters are arranged for each pixel in a predetermined arrangement (Bayer arrangement, honeycomb arrangement). It is.

光源装置300は、可視光を図示しないライトガイドに入射させる。ライトガイドの一端はLGコネクタ120を介して光源装置300に接続され、ライトガイドの他端は照明光学系152に対面している。光源装置300から発せられた光は、ライトガイドを経由して照明光学系152から出射され、観察光学系150の視野範囲を照明する。   The light source device 300 causes visible light to enter a light guide (not shown). One end of the light guide is connected to the light source device 300 via the LG connector 120, and the other end of the light guide faces the illumination optical system 152. The light emitted from the light source device 300 is emitted from the illumination optical system 152 via the light guide, and illuminates the visual field range of the observation optical system 150.

内視鏡プロセッサ200には、CCDから出力される画像信号が電気コネクタ110を介して入力される。このアナログの画像信号は、内視鏡プロセッサ200内においてデジタルの画像信号に変換され、画像表示部500の画面に表示するための必要な処理が施される。   An image signal output from the CCD is input to the endoscope processor 200 via the electrical connector 110. The analog image signal is converted into a digital image signal in the endoscope processor 200, and necessary processing for displaying on the screen of the image display unit 500 is performed.

このように、内視鏡100で得られた観察画像のデータが内視鏡プロセッサ200に出力され、内視鏡プロセッサ200に接続された画像表示部500に画像が表示される。   In this manner, observation image data obtained by the endoscope 100 is output to the endoscope processor 200, and an image is displayed on the image display unit 500 connected to the endoscope processor 200.

以上、本発明の光断層画像化装置及びその作動方法について詳細に説明したが、本発明は、以上の例には限定されず、本発明の要旨を逸脱しない範囲において、各種の改良や変形を行ってもよいのはもちろんである。
The optical tomographic imaging apparatus and the operation method thereof according to the present invention have been described in detail above. However, the present invention is not limited to the above examples, and various improvements and modifications can be made without departing from the gist of the present invention. Of course you can go.

1…光断層画像化装置、2…光源ユニット、2A、2B…光源、3A、3B…光分割部、6A、6B…合波部、8A、8B…干渉光検出部、11…表示装置、20…演算部、21…第1干渉波データ格納部、22…第2干渉波データ格納部、23…第1光断層画像データ生成部、24…第2光断層画像データ生成部、25…ピーク位置抽出部、26…第2光断層画像データ補正部、27…第1干渉波データ再構築部、28…第2干渉波データ再構築部、29…合成干渉波生成部、30…FFT部、600…光プローブ   DESCRIPTION OF SYMBOLS 1 ... Optical tomographic imaging apparatus, 2 ... Light source unit, 2A, 2B ... Light source, 3A, 3B ... Light splitting part, 6A, 6B ... Multiplexing part, 8A, 8B ... Interference light detection part, 11 ... Display apparatus, 20 ... Calculation unit, 21 ... first interference wave data storage unit, 22 ... second interference wave data storage unit, 23 ... first optical tomographic image data generation unit, 24 ... second optical tomographic image data generation unit, 25 ... peak position Extraction unit 26 ... second optical tomographic image data correction unit 27 ... first interference wave data reconstruction unit 28 ... second interference wave data reconstruction unit 29 ... combined interference wave generation unit 30 ... FFT unit 600 ... Optical probe

Claims (5)

互いに波長帯域が異なり、それぞれ所定の前記波長帯域内で波長を掃引して複数の光束を射出する光源ユニットと、
前記光源ユニットから射出された前記複数の各光束をそれぞれ測定光と参照光とに分割する光分割手段と、
前記光分割手段により分割された複数の前記測定光を測定対象に照射する手段であって、前記測定光を導く光ファイバと、前記光ファイバの先端に固定され、前記光ファイバの先端から出射した前記測定光を前記測定対象に向けて偏向させる先端光学系とを有する測定光照射手段と、
前記先端光学系と前記測定対象との間に配置され、前記測定光の一部を反射し他の一部を透過する基準反射手段と、
前記測定光が測定対象に照射されたときの前記測定対象からの測定反射光及び前記基準反射手段からの基準反射光と、前記参照光とのそれぞれの干渉光信号を前記光束毎に検出する複数の干渉光検出手段と、
前記干渉光検出手段にて検出された前記光束毎の前記干渉光信号をフーリエ変換して得られる光断層画像データのピーク位置を検出するピーク位置抽出手段と、
前記複数の光束のうち一の光束についての光断層画像データのピーク位置と、他の光束についての光断層画像データのピーク位置との差をピーク情報として算出するピーク情報算出手段と、
前記ピーク情報に基づいて前記他の光束についての光断層画像データを補正する光断層画像データ補正手段と、
前記一の光束についての光断層画像データ、及び前記光断層画像データ補正手段による補正後の光断層画像データを逆フーリエ変換してそれぞれの干渉光信号を前記光束毎に再構築する干渉光再構築手段と、
を備えたことを特徴とする光断層画像化装置。
A light source unit that emits a plurality of light fluxes, each having a different wavelength band and sweeping a wavelength within a predetermined wavelength band;
A light splitting means for splitting the plurality of light beams emitted from the light source unit into measurement light and reference light, respectively;
A means for irradiating a measurement object with a plurality of the measurement lights divided by the light dividing means, and an optical fiber for guiding the measurement light, fixed to the tip of the optical fiber, and emitted from the tip of the optical fiber Measuring light irradiation means having a tip optical system that deflects the measuring light toward the measurement object;
Reference reflecting means that is disposed between the tip optical system and the measurement object, reflects a part of the measurement light and transmits the other part,
A plurality of interference light signals for each of the light fluxes, which are measured reflected light from the measuring object and the reference reflected light from the reference reflecting means when the measuring light is applied to the measuring object, and the reference light. Interference light detecting means,
Peak position extraction means for detecting a peak position of optical tomographic image data obtained by Fourier transforming the interference light signal for each of the light beams detected by the interference light detection means;
Peak information calculating means for calculating, as peak information, the difference between the peak position of the optical tomographic image data for one light beam among the plurality of light beams and the peak position of the optical tomographic image data for the other light beam;
Optical tomographic image data correction means for correcting optical tomographic image data for the other light flux based on the peak information;
Interfering light reconstruction for reconstructing each interference light signal for each light flux by inverse Fourier transforming the optical tomographic image data for the one light flux and the optical tomographic image data corrected by the optical tomographic image data correction means Means,
Optical tomography system characterized in that example Bei a.
前記ピーク位置抽出手段は、前記測定光照射手段による前記測定光の照射方向に関して、前記光断層画像データが所定閾値を超える位置であって前記先端光学系に最も近い位置を前記ピーク位置として検出することを特徴とする請求項に記載の光断層画像化装置。 The peak position extraction means detects, as the peak position, a position where the optical tomographic image data exceeds a predetermined threshold and is closest to the tip optical system with respect to the measurement light irradiation direction by the measurement light irradiation means. The optical tomographic imaging apparatus according to claim 1 . 前記干渉光再構築手段によって前記光束毎にそれぞれ再構築された前記干渉光信号を合成した合成干渉光信号を生成する合成干渉光信号生成手段と、
前記合成干渉光信号より前記測定対象の光断層像情報を生成する光断層情報生成手段と、
をさらに備えたことを特徴とする請求項又はに記載の光断層画像化装置。
Combined interference light signal generating means for generating a combined interference light signal obtained by combining the interference light signals reconstructed for each of the light beams by the interference light reconstructing means;
Optical tomographic information generating means for generating optical tomographic image information of the measurement object from the combined interference light signal;
Optical tomographic imaging apparatus according to claim 1 or 2, further comprising a.
前記複数の光束の波長掃引周期は、同期されていることを特徴とした請求項1ないしのいずれか1つに記載の光断層画像化装置。 Wherein the plurality of wavelength sweep cycle of the light beam, the optical tomographic imaging apparatus according to any one of claims 1 to 3 characterized by being synchronized. 互いに波長帯域が異なり、それぞれ所定の前記波長帯域内で波長を掃引して複数の光束を射出する光源ユニットと、
前記光源ユニットから射出された前記複数の各光束をそれぞれ測定光と参照光とに分割する光分割手段と、
前記光分割手段により分割された複数の前記測定光を測定対象に照射する手段であって、前記測定光を導く光ファイバと、前記光ファイバの先端に固定され、前記光ファイバの先端から出射した前記測定光を前記測定対象に向けて偏向させる先端光学系とを有する測定光照射手段と、
前記先端光学系と前記測定対象との間に配置され、前記測定光の一部を反射し他の一部を透過する基準反射手段と、
複数の干渉光検出手段と、ピーク位置抽出手段と、ピーク情報算出手段と、光断層画像データ補正手段と、干渉光再構築手段と、を備える光断層画像化装置の作動方法であって、
前記複数の干渉光検出手段が、前記測定光が測定対象に照射されたときの前記測定対象からの測定反射光及び前記基準反射手段からの基準反射光と、前記参照光とのそれぞれの干渉光信号を前記光束毎に検出
前記ピーク位置抽出手段が、前記光束毎の前記干渉光信号をフーリエ変換して得られる光断層画像データのピーク位置を検出
前記ピーク情報算出手段が、前記複数の光束のうち一の光束についての光断層画像データのピーク位置と、他の光束についての光断層画像データのピーク位置との差をピーク情報として算出
前記光断層画像データ補正手段が、前記ピーク情報に基づいて前記他の光束についての光断層画像データを補正
前記干渉光再構築手段が、前記一の光束についての光断層画像データ、及び前記補正後の光断層画像データを逆フーリエ変換してそれぞれの干渉光信号を前記光束毎に再構築することを特徴とする光断層画像化装置の作動方法。
A light source unit that emits a plurality of light fluxes, each having a different wavelength band and sweeping a wavelength within a predetermined wavelength band;
A light splitting means for splitting the plurality of light beams emitted from the light source unit into measurement light and reference light, respectively;
A means for irradiating a measurement object with a plurality of the measurement lights divided by the light dividing means, and an optical fiber for guiding the measurement light, fixed to the tip of the optical fiber, and emitted from the tip of the optical fiber Measuring light irradiation means having a tip optical system that deflects the measuring light toward the measurement object;
Reference reflecting means that is disposed between the tip optical system and the measurement object, reflects a part of the measurement light and transmits the other part,
An operation method of an optical tomographic imaging apparatus comprising a plurality of interference light detection means, peak position extraction means, peak information calculation means, optical tomographic image data correction means, and interference light reconstruction means,
The plurality of interference light detection means are interference lights of the measurement reflected light from the measurement object, the reference reflected light from the reference reflection means, and the reference light when the measurement light is irradiated onto the measurement object. A signal is detected for each luminous flux,
The peak position extraction means detects a peak position of the optical tomographic image data obtained by the interference light signal for each of the light beams by Fourier transform,
The peak information calculation means calculates as the peak information a difference between the peak position of the optical tomographic image data for light and the peak position of the tomographic image data, the other light beam for one of the light beams of the plurality of light beams,
The optical tomographic image data correcting means, an optical tomographic image data for the other light beam is corrected on the basis of the peak information,
The interference light reconstruction means, to rebuild the respective interference light signal for each of the beam optical tomographic image data, and an optical tomographic image data before and after Kiho forward and inverse Fourier transform for the one light beam method of operating an optical tomography system according to claim and this.
JP2009297876A 2009-12-28 2009-12-28 Optical tomographic imaging apparatus and operating method thereof Active JP5616626B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009297876A JP5616626B2 (en) 2009-12-28 2009-12-28 Optical tomographic imaging apparatus and operating method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009297876A JP5616626B2 (en) 2009-12-28 2009-12-28 Optical tomographic imaging apparatus and operating method thereof

Publications (2)

Publication Number Publication Date
JP2011137717A JP2011137717A (en) 2011-07-14
JP5616626B2 true JP5616626B2 (en) 2014-10-29

Family

ID=44349279

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009297876A Active JP5616626B2 (en) 2009-12-28 2009-12-28 Optical tomographic imaging apparatus and operating method thereof

Country Status (1)

Country Link
JP (1) JP5616626B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2852315A4 (en) * 2012-05-21 2016-06-08 Gen Hospital Corp Apparatus, device and method for capsule microscopy
CA2880038C (en) * 2012-07-27 2021-06-01 Thorlabs, Inc. Agile imaging system
JP6281059B2 (en) * 2013-07-31 2018-02-21 アダマンド並木精密宝石株式会社 Optical imaging probe
CN115078264B (en) * 2022-08-22 2022-11-18 广东大湾区空天信息研究院 Polarization sensitive optical coherence tomography system and method

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006215006A (en) * 2005-02-07 2006-08-17 Fujinon Corp Optical tomographic imaging system
JP4677636B2 (en) * 2005-12-13 2011-04-27 日本電信電話株式会社 Optical coherence tomography apparatus and variable wavelength light generator used therefor
JP4864511B2 (en) * 2006-03-31 2012-02-01 富士フイルム株式会社 Electronic endoscope apparatus and program
JP2008089349A (en) * 2006-09-29 2008-04-17 Fujifilm Corp Optical tomographic imaging system
JP4869877B2 (en) * 2006-11-17 2012-02-08 富士フイルム株式会社 Optical tomographic imaging system

Also Published As

Publication number Publication date
JP2011137717A (en) 2011-07-14

Similar Documents

Publication Publication Date Title
JP4986296B2 (en) Optical tomographic imaging system
JP4642681B2 (en) Optical tomographic imaging system
JP5069585B2 (en) Optical tomographic imaging system using an optical probe
JP4963708B2 (en) Optical coherence tomography device
US7835010B2 (en) Optical coherence tomography system and optical coherence tomography method
US8599386B2 (en) Optical tomographic imaging apparatus, interference signal processing method, and endoscope apparatus
JP4895277B2 (en) Optical tomographic imaging system
US20240065552A1 (en) Intraoral oct with color texture
WO2005031320A1 (en) Variable-wavelength light generator and light interference tomograph
JP2008253492A (en) Tomographic image processing method, apparatus and program
JP2008128709A (en) Optical tomographic imaging apparatus
JP2008142443A (en) Optical tomographic apparatus
JP2007275193A (en) Optical probe and optical tomographic imaging equipment
JP2007085931A (en) Optical tomographic imaging system
JP2007101263A (en) Optical tomographic imaging device
JP5373485B2 (en) OCT apparatus and interference signal level control method thereof
JP2009072280A (en) Optical tomography imaging system, contact area detecting method, and image processing method using the same
JP5616626B2 (en) Optical tomographic imaging apparatus and operating method thereof
JP5037215B2 (en) Compensation table generation method, apparatus, program, and tomographic image processing apparatus using the same
JP2008275529A (en) Method, apparatus and program for processing tomographic image
JP2008253493A (en) Tomographic image processing method, apparatus and program
JP2012013520A (en) Optical tomographic imaging device and method for estimating pressing force of its optical probe
JP2008128707A (en) Tomographic image processing method, device and program, and optical tomographic imaging system using it
JP2022524547A (en) High Speed Dental Optical Coherence Tomography System
JP2012200283A (en) Optical coherence tomographic imaging apparatus

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20120514

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20130528

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130529

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130628

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130823

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140320

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140513

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140606

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140718

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140912

R150 Certificate of patent or registration of utility model

Ref document number: 5616626

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250