The invention starts with a method for signal processing in optical coherence tomography by means of a tunable light source (swept source), wherein the method comprises the steps:
tuning the light source and sensing a signal intensity of the light source in linear dependence on the respective wave number (k) of the tunable light source and producing a signal intensity distribution in dependence on k;
applying a window function to the sensed signal intensity distribution and producing a weighted signal intensity distribution; and
applying a Fast Fourier Transform (FFT) to the weighted signal intensity distribution.
This section provides background information related to the present disclosure which is not necessarily prior art.
The window function leads to portions of the signal not being used, since they are “windowed out,” i.e., filtered out. To increase the usable spectrum range and to amplify the signal intensity, it is desirable to minimize the lost proportion by applying the window function.
This section provides a general summary of the disclosure, and is not a comprehensive disclosure of its full scope or all of its features.
To achieve this, it is proposed to limit the tuned frequency spectrum to a passband of the window function in the tuning of the light source.
The window function applied can be a von Hann window. In principle, a light source of which the spectrum shape corresponds to the curve of the window function would be optimal. In order to obtain at the same time, the best possible side wall suppression, a Gaussian light distribution in the k-space can be used as a suitable compromise. Thus, the tunable light source can have a Gaussian light intensity distribution in the k-space. In particular, the spectrum shape of the light source in the k-space can correspond to the shape of the curve of the window function or at least approximate said curve.
In the same way, for a determination of the spectrum shape of the light source in the k-space, the temporal curve of the signal intensity and the sweep speed curve can be measured. The temporal curve of the signal intensity can be recorded, for example, with a photodiode.
The sweep speed curve can be determined, for example, with the aid of a Mach-Zehnder interferometer (MZI), in that two fiber Bragg gratings of known and differing reflection wavelength are arranged in the beam propagation path between interferometer and detector, and a modulation frequency of the tunable light source, which is proportional to the sweep speed, is sensed and plotted versus time, so that the modulation frequency curve has two dips in each case at a certain wavelength or wave number, and wherein, the sweep speed curve can be determined from the obtained MZI signal curve by Short Time Fourier Transform (STFT).
The envelope of the modulation frequency curve can be determined, and its deviation from a reference shape from the window function, in particular from the Gaussian shape, can be sensed. Moreover, for an adaptation of the envelope of the modulation frequency curve to the reference shape of the window function, the current of a Semiconductor Optical Amplifier (SOA) of the light source can be modulated or the sweep speed curve can be varied.
The data sensed for the determination of the envelope of the modulation frequency curve can be used via a feedback loop for the fine tuning of the swept source, in particular of a VCSEL light source.
Further areas of applicability will become apparent from the description provided herein. The description and specific examples in this summary are intended for purposes of illustration only and are not intended to limit the scope of the present disclosure.
The drawings described herein are for illustrative purposes only of selected embodiments and not all possible implementations, and are not intended to limit the scope of the present disclosure.
Additional details of the invention are explained in reference to the following figures. In the drawings:
Example embodiments will now be described more fully with reference to the accompanying drawings.
The spectrum shape of the tunable light source in the k-space is not directly measurable. However, the measurement arrangements reproduced diagrammatically in
In a first step, by exposure of the signal to a photodiode 1 or to other light-sensitive measurement cells, the envelope of the light intensity of the tunable light source versus time can be sensed.
The measurement of the sweep speed curve occurs with the aid of a Mach-Zehnder interferometer (MZI) 3, in that the modulation frequency, which is proportional to the sweep speed, is analyzed and plotted versus time. In addition, two sites in the sweep are marked, in that two fiber Bragg gratings 5 of known and differing reflection wavelength are arranged between interferometer 4 and detector 3. By means of said gratings, signal dips in the MZI raw signal are produced, which correspond to a corresponding wavelength or wave number.
From the MZI signal, with the aid of a STFT, a speed curve can be produced. From the data obtained as described above, a conclusion can be drawn as to the power distribution in the k-space. This makes it possible to determine the deviation of the spectrum shape of the tunable light source from the optimal Gaussian shape.
For the optimization of the spectrum shape of the tunable light source, for example, the current of a semiconductor optical amplifier (SOA) of the light source can then be modulated or the sweep speed curve can be varied. The two ways are in principle equally suitable, although a speed of the sweep that is as constant as possible is considered to be advantageous in order to be able to save detector bandwidth, which in turn improves the measurement sensitivity.
The dated obtained by means of the above-described method can be used via feedback loop for automatic fine tuning of the VCSEL light source or any other suitable swept source.
The features of the invention disclosed in the above description, in the drawings and in the claims can be essential both individually and also in any combination for the implementation of the invention.
Number | Date | Country | Kind |
---|---|---|---|
10 2016 103 605.6 | Mar 2016 | DE | national |
This application is a 371 U.S. National Stage of International Application No. PCT/DE2017/100026, filed on Jan. 19, 2017, which claims priority to German Application 10 2016 103 605.6, filed Mar. 1, 2016. The entire disclosures of the above applications are incorporated herein by reference.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/DE2017/100026 | 1/19/2017 | WO | 00 |