The present invention relates to a wavelength swept light measurement technique for measuring a profile of a sweep frequency width and the like for a wavelength swept light source.
Depth distance measurement methods including a frequency modulated continuous wave (FMCW) radar method, a swept source optical coherence tomography (SS-OCT) method, and an optical frequency domain reflectometry (OFDR) method use a wavelength swept light source, an interferometer, and a signal processing apparatus. A resolution of the depth distance in these methods is inversely proportional to a sweep frequency width of the wavelength swept light sources. Thus, in configuring a depth distance measurement apparatus using such methods, it is important to measure the sweep frequency width and a swept wavelength width to ensure the performance of the wavelength swept light source.
In measuring the sweep frequency width, an optical spectrum analyzer is typically employed to acquire a sweeping spectrum and obtain the sweep frequency width. NPLs 1 and 2 describe a dispersion spectroscopy and an interference spectroscopy of an optical spectrum analyzer.
In the dispersion spectroscopy, a dispersion element such as a prism or a diffraction grating, a slit, a photodetector, and a signal processing apparatus are used. When light to be measured enters the dispersion element, the light is divided (dispersed) according to an emission direction from the dispersion element and a location of the light depending on a wavelength of the light, and light of a certain wavelength is separated from the light dispersed by the slit, to obtain the intensity of the separated light of the certain wavelength by the photodetector. At this time, when the incident angle of the light to be measured on the dispersion element is changed by, for example, rotating the dispersion element, a wavelength passing the slit changes, and a wavelength-to-light intensity property (wavelength spectrum) is obtained by the signal processing apparatus, from a relationship between the incident angle and the wavelength passing the slit.
The interference spectroscopy includes two-beam interferometry using a Michelson interferometer, a Mach-Zehnder interferometer, or the like, and multi-beam interferometry using a Fabry-Perot interferometer.
In the two-beam interferometry, a two-beam interferometer, a photodetector, and a signal processing apparatus are used. In the interferometer, light to be measured is demultiplexed, and after each of the demultiplexed light beams passes through each of two optical paths, the demultiplexed light beams are multiplexed. A function of a multiplexed light intensity with optical path lengths of the two optical paths as variables is a sum of an autocorrelation function of the light to be measured and the intensities of each of the light beams. The multiplexed light intensity is acquired by the photodetector and subsequently Fourier-transformed by the signal processing apparatus, to obtain a power spectrum of the light to be measured.
In the multi-beam interferometry, a scanning Fabry-Perot interferometer is used. The scanning Fabry-Perot includes two high-reflective mirrors faced each other and being capable of continuously changing an interval between the high-reflective mirrors, a photodetector, and a signal processing apparatus are used.
Transmission spectral characteristics of the scanning Fabry-Perot interferometer have comb-teeth-shaped transmission characteristics in which light is passed at equal frequency intervals (with a free spectral range (FSR)). Each of the comb teeth is referred to as a vertical mode, and the frequency of each vertical mode is denoted by nFSR (n being an integer of 0 or greater). This n is referred to as a vertical mode order. The FSR varies in inverse proportion to the interval between the two high-reflective mirrors of the Fabry-Perot interferometer, and thus, this property is utilized by the optical spectrum analyzer.
When the frequency band of the light to be measured is smaller than the FSR, light passing only one order is passed from the Fabry-Perot interferometer. When the interval between the two high-reflective mirrors of the Fabry-Perot interferometer varies, the frequency of the order varies in inverse proportion to the variation. After the photodetector acquires the light intensity of the light passing the Fabry-Perot interferometer, the signal processing apparatus associates the light intensity and the frequency of the order to obtain a wavelength-to-light intensity property (wavelength spectrum). When only one Fabry-Perot interferometer is used in obtaining a spectrum having a wider band than in FSR, light passing a plurality of orders of vertical modes enters the photodetector. In order to avoid this, for example, it is necessary to use a device in which a plurality of Fabry-Perot interferometers are connected in series to limit the number of orders to one order.
NPL 1: Akihisa Mikami, “Multimedia-Its Impact and Technology. Technologies and Applications of Speech Recognition and Speech Synthesis” Journal of the Society of Instrument and Control Engineers, Vol. 35, No. 1, pp. 29-32, published Jan. 10, 1996
NPL 2: Takao Tanimoto, “Various Measurement Technology in Laser Experiment—Principle of Optical Spectrum Analyzer and Notes in Measurement—”, The Review of Laser Engineering, Vol. 39, No. 5, pp. 354-361, published May 15, 2011.
The measurement apparatus described above acquires a temporally static optical spectrum, and is not suitable for a wavelength swept light source and the like that an optical spectrum varies over time, so that it is inherently difficult to accurately acquire a sweeping spectrum by the measurement apparatus. As a result, there is a problem in that the measurement apparatus cannot accurately measure a sweep frequency width of the wavelength swept light source.
Furthermore, when the sweeping speed of the wavelength swept light source is sufficiently fast with respect to a frequency sweeping speed of the optical spectrum analyzer in the dispersion spectroscopy, the swept wavelength width can be measured with a certain degree of accuracy. However, as the sweeping speed of the wavelength swept light source approaches the wavelength measurement speed of the optical spectrum analyzer, an unevenness occurs in the observed spectrum, and thus, there occurs a problem that it is not possible to accurately measure the sweep frequency width.
In order to clearly explain this problem, results of measuring a spectrum (wavelength spectrum) of a wavelength swept light source by using an optical spectrum analyzer employing the dispersion spectroscopy are illustrated in
As illustrated in
The wavelength swept light source employed in the spectral measurement has a configuration in which a current obtained by adding a direct current and an alternating current is applied to a distributed feedback (DFB) laser. The alternating waveform has a sine waveform, and a frequency of the sine waveform is a sweep frequency of the light source. During the measurement, the sweep frequency of the optical spectrum analyzer being used is about 1 Hz.
As illustrated in
Note that the reason why the swept wavelength (sweep frequency) width is narrower at 100 kHz than at 10 Hz is due to the properties of a DFB laser. Furthermore, as seen in the results at 100 kHz, the light intensity level at both ends of the swept wavelength rises like a horn because the alternating waveform of the current applied to the DFB laser has a sine waveform, and thus, the time transition of the wavelength becomes slow toward the end of the swept wavelength, so that in the photodetector receiving the light inside the spectrum analyzer, the irradiation time per unit wavelength increases.
Embodiments of the present invention have been contrived to solve the above problem and an object thereof is to provide a wavelength swept light measurement technique capable of accurately measuring a profile such as a sweep frequency width for wavelength swept light, without being affected by a sweeping speed of a wavelength swept light source.
In order to achieve this object, a wavelength swept light measurement system according to embodiments of the present invention includes a photoelectric conversion apparatus including an interferometer configured to cause interference in wavelength swept light output from a wavelength swept light source, the photoelectric conversion apparatus being configured to convert the interfered wavelength swept light by photoelectric conversion, and a signal processing apparatus configured to calculate, in chronological order, relative frequencies indicating frequencies relative to interference signals obtained in the photoelectric conversion, and measure a difference between a maximum value and a minimum value of the relative frequencies as a sweep frequency width of the wavelength swept light.
Furthermore, another wavelength swept light measurement system according to embodiments of the present invention includes a photoelectric conversion apparatus including an interferometer configured to cause interference in wavelength swept light output from a wavelength swept light source, the photoelectric conversion apparatus being configured to convert, by photoelectric conversion, the interfered wavelength swept light and specific wavelength light detected by a narrow band wavelength filter from the wavelength swept light, and a signal processing apparatus configured to calculate, in chronological order, relative frequencies indicating frequencies relative to interference signals obtained by the photoelectric conversion of the interference light, calculate predicted frequencies indicating absolute frequencies for the relative frequencies, by referring to a detection timing of the specific wavelength light obtained by the photoelectric conversion of the specific wavelength light, and measure a difference between a maximum value and a minimum value of predicted wavelengths corresponding to the predicted frequencies, as a swept wavelength width of the wavelength swept light.
According to embodiments of the present invention, it is possible to measure a profile for wavelength swept light with extremely high precision and accuracy, without being affected by a sweeping speed of a wavelength swept light source.
Embodiments of the present invention will be described below with reference to the drawings.
First, a wavelength swept light measurement system 100 according to a first embodiment of the present invention will be described with reference to
The wavelength swept light measurement system 100 according to embodiments of the present invention is a system for measuring a profile such as a swept wavelength width Δλ for a wavelength swept light source X, on the basis of wavelength swept light Lx, to be measured, output from the wavelength swept light source X.
The wavelength swept light source X is a laser light source capable of sweeping a wavelength of oscillating light at high speed and in a wide range. In the wavelength swept light source X, a time-optical frequency pattern is the same for every sweep operation, and a trigger electrical signal Trg is output in synchronization with the sweep operation.
As illustrated in
The photoelectric conversion apparatus 10 includes an interferometer ii that generates interference light iL by causing an interference in the wavelength swept light Lx, to be measured, output from the wavelength swept light source X, and a balanced photodetector 12 that outputs an interference electrical signal iE(t) by photoelectrically converting the obtained interference light iL.
The interferometer 11 illustrated in
The balanced photodetector 12 is a typical differential amplification-type photodetector and differentially amplifies and photoelectrically converts the two light beams of the interference light iL branched by the C1, to output the interference electrical signal iE(t).
The interferometer 11 is not limited to the Mach-Zehnder type interferometer, and may be a Michelson-type or a Fabry-Perot interferometer. In the case of the Fabry-Perot type interferometer, a normal photodetector may be used instead of the balanced photodetector 12. Furthermore, in
The A/D converter (ADC) 30 is a typical A/D converter having a trigger function. The ADC 30 A/D-converts a trigger electrical signal trE(t) that is the trigger electrical signal Trg of the wavelength swept light source X and input from the photoelectric conversion apparatus 10, at each time t, and outputs, in chronological order, trigger signals tr(t) composed of digital data. Furthermore, the ADC 30 A/D-converts the interference electrical signal iE(t) input from the photoelectric conversion apparatus 10, at each time t, and outputs, in chronological order, interference signals i(t) composed of digital data.
At this time, the ADC 30 includes therein a memory 31 and saves, in the memory 31 in chronological order, the interference signals i(t) during an effective period Tmem having a duration of one or more sweeps from a sweep start specified by the trigger electrical signal trE(t), and in response to a request from the signal processing apparatus 20, reads and outputs the interference signals i(t) for a specified period of the effective period Tmem from the memory 31. The memory 31 may be an external memory externally connected to the ADC 30.
The signal processing apparatus 20 includes a microprocessor such as a digital service unit (DSU) and a CPU. The signal processing apparatus 20 realizes, when the microprocessor and a program stored in a storage apparatus 32 are jointly operated, various types of signal processes for calculating a profile for the wavelength swept light Lx of the wavelength swept light source X, that is, in this case, a sweep frequency width Δf.
As illustrated in
The target extraction unit 21 of the signal processing apparatus 20 illustrated in
Specifically, first, the target extraction unit 21 detects an output timing at which the trigger electrical signal Trg is output from the wavelength swept light source X, on the basis of an intensity of the trigger signal tr(t) from the ADC 30. Subsequently, the target extraction unit 21 identifies, each time the output timing is detected, the sweep interval Taq by referring to a trigger time Ttrg indicating the output timing, and extracts, from the ADC 30, the target interference signal it(t) corresponding to the sweep interval Taq.
Thus, it is only required in the ADC 30 that the interference signal i(t) is sequentially updated and held in the memory 31 of the ADC 30 for the predetermined effective period Tmem, for example. Consequently, when the target extraction unit 21 reads the interference signal i(t) from the memory 31 of the ADC 30 during an interval from the maximum frequency to the minimum frequency of the wavelength swept light Lx, that is, during the sweep interval Taq, the target extraction unit 21 can extract the target interference signal it(t).
A signal extraction operation in the target extraction unit 21 will be described with reference to
In
The effective period Tmem has a duration longer than a wavenumber period of the wavelength swept light Lx, and thus, even when the trigger time Ttrg shifts before and behind to some extent, it is possible to stably extract the target interference signal it(t) during the sweep interval Taq from the interference signal i(t).
The relative frequency calculation unit 22 of the signal processing apparatus 20 illustrated in
On the basis of the target interference signal it(t) for the sweep interval Taq extracted by the target extraction unit 21, the relative frequency calculation unit 22 is configured to calculate a difference at each time t between a reference frequency f(t2) at a reference time and a frequency f(t) of the target interference signal it(t) at each time, where the reference time is an initial time of the sweep interval Taq, that is, the time t2 at which the maximum frequency is output in
A relative frequency calculation operation by the relative frequency calculation unit 22 will be described with reference to
In
Math. 1
f
r(t)=f(t)−f(0)=f(t)−f0 (1)
where f0 is actually a value of the frequency f(0) of light output from a light source at a time t=0.
The relative frequency calculation unit 22 determines the relative frequency fr(t) from the phase of the target interference signal it(t). This is on the basis of the following principle.
According to reference literature (Yoshiaki Yasuno, Violeta Dimitrova Madjarova, Shuichi Makita, Masahiro Akiba, Atsushi Morosawa, Changho Chong, Toni Sakai, Kin-Pui Chan, Masahide Itoh, and Toyohiko Yatagai, “Three-dimensional and high-speed swept-source optical coherence tomography for in vivo investigation of human anterior eyesegments”, OPTICS EXPRESS, Vol. 13, No. 26, pp. 10652-10664, 2005), the interference signal i(t) is expressed by Equation (2) below.
where z is a difference of the optical path lengths in the interferometer.
In Equation (2), a(t) is an amplitude and is, according to the reference literature, proportional to the photon efficiency of the photodetector, the coherence function of the light of the light source, and the magnitude of the electric field of the two optical paths in the interferometer. Furthermore, c is a speed of light. When the amplitude a(t) is constant, a phase θ(t) of the interference signal i(t) is expressed by Equation (3) below. Consequently, if the phase θ(t) is determined, a frequency f(x) can be obtained by Equation (4) below.
Incidentally, the interference signal i(t) may be deformed as in Equation (5) below, and when the amplitude a(t) is substantially a constant value a0, the interference signal i(t) is given by Equation (6) below.
where j is in units of imaginary numbers.
The interference signal i(t) can be expressed by Equation (7) below. If the exponential part of the Napiar number e in Equation (7) is compared with that in Equation (6), the instantaneous frequency fi(t) of the interference signal i(t) is expressed by Equation (8) below.
where fi(t) is an instantaneous frequency of the interference signal i(t) at a time t.
Thus, it can be understood that the interference signal i(t) is a signal having a frequency of ±(z/c) df/dt|t. A signal i′(t) that is a signal obtained by deleting a negative frequency component from the interference signal i(t)is expressed by Equation (9) below.
As in Equation (9), the signal i′(t) is a complex number, and thus, the phase θ(t)=(4πz/c) f(t) is determined as an argument of the signal i′(t). The phase θ(t) is expressed by Equation (10) below.
Math. 10
θ(t)=arg(i′(t)) (10)
where arg(.) is a function for determining the argument from a complex number.
It should be noted that an argument obtained from the actual signal i′(t) can only be obtained from a principal value (value in a range from 0 to 2π). That is, the argument obtained from the actual signal i′(t) is not the original value, but a value wrapped in the range from 0 to 2π is obtained. Here, an original argument θ(t) is expressed by Equation (11) below.
Math. 11
θ(t)=θr(t)+θ0=unwrap(Arg(i′(t)))+θ0 (11)
where θ(0)=θ0 is the original argument at time t=0, Arg(.) is a function for obtaining the principal value of the argument from a complex number, and unwrap(.) is a function for unwrapping a value.
Note that, an argument not including θ0 is referred to as a relative phase θr(t) so as to distinguish from the original phase.
Consequently, from Equations (4) and (11), the frequency f(x) of the wavelength swept light Lx is expressed by Equation (12) below.
In comparison between Equation (1) and Equation (12), it can be understood that the relative frequency fr(t) and the frequency f0 (=f(0)) of the light output from the light source at time t=0 correspond to each other as in Equations (13) below.
Thus, it can be understood that the principal value of the argument of the signal i′(t) obtained by removing the negative frequency component from the interference signal i(t) is determined, and the argument is multiplied by c/4πz to obtain the relative frequency fr(t) of the wavelength swept light Lx.
An exemplary configuration of the relative frequency calculation unit 22 in the signal processing apparatus 20 in
In the relative frequency calculation unit 22, the negative frequency component deletion unit 41 deletes a negative frequency component from the input interference signal i(t) and outputs the signal i′(t). The argument calculation unit 42 determines and outputs a relative argument θr(t) of the signal i′(t) output from the negative frequency component deletion unit 41. The argument-to-frequency conversion unit 43 determines and outputs the relative frequency fr(t) of the wavelength swept light Lx, on the basis of the relative argument θr(t) output from the argument calculation unit 42.
An exemplary configuration of the negative frequency component deletion unit 41 in the relative frequency calculation unit 22 in
In the negative frequency component deletion unit 41, the Fourier transform unit 45 transforms the input interference signal i(t) by Fourier transform (frequency transformation), and outputs a signal i(fi) resulting from the Fourier transform. The negative frequency component replacement unit 46 replaces a negative frequency component of the signal i(fi) output from the Fourier transform unit 45 with zero, and outputs a signal i′(fi) resulting from the replacement. The inverse Fourier transform unit 47 transforms the signal i′(fi) output from the negative frequency component replacement unit 46 by inverse Fourier transform, and outputs a signal i′(t) resulting from the inverse Fourier transform.
The sweep frequency width measurement unit 23 of the signal processing apparatus 20 illustrated in
As described above, in the present embodiment, the photoelectric conversion apparatus 10 photoelectrically converts the interference light iL obtained as a result of interference of the wavelength swept light Lx output from the wavelength swept light source X in the interferometer 11, and outputs the photoelectrically converted interference light iL. The signal processing apparatus 20 calculates, in chronological order, relative frequencies fr(t) indicating frequencies relative to the interference signals i(t) obtained by the photoelectric conversion of the interference light iL, and measures a difference between a maximum value and a minimum value of the relative frequencies fr(t), as the sweep frequency width Δf of the wavelength swept light Lx.
Consequently, unlike a case where the sweep frequency width of a wavelength swept light source is measured on the basis of an optical spectrum acquired by dispersion spectroscopy, it is possible to measure the sweep frequency width Δf of the wavelength swept light Lx with extremely high precision and accuracy, regardless of a sweeping speed of the wavelength swept light source X.
Furthermore, in the present embodiment, in the signal processing apparatus 20, the target extraction unit 21 may extract, from the interference signals i(t), a target interference signal it(t) corresponding to the sweep interval Taq in which the wavelength swept light Lx is swept from a maximum frequency to a minimum frequency, the relative frequency calculation unit 22 may calculate a difference between the frequency of the extracted target interference signal it(t) at a reference time and the frequency of the target interference signal it(t) at each time, as the relative frequencies ft(t), and the sweep frequency width measurement unit 23 may measure a difference between a maximum value and a minimum value of the obtained relative frequencies fr(t), as the sweep frequency width Δf of the wavelength swept light Lx. Consequently, the sweep frequency width Δf of the wavelength swept light Lx can be measured with a relatively simple process.
Furthermore, there is further provided the ADC 30 that A/D-converts the interference electrical signal iE(t) output from the photoelectric conversion apparatus 10 into an interference signal i(t) composed of digital data in the present embodiment, and the ADC 30 is required to A/D-convert the interference electrical signal iE(t) at each time t during the effective period Tmem having a duration of one or more sweeps from a start of sweeping of the wavelength swept light Lx, save the obtained interference signal i(t) into the memory 31 in chronological order, and, in response to a request from the signal processing apparatus 20, read the interference signal i(t) for a specified period of the effective period Tmem from the memory 31, to output the read interference signal i(t). Consequently, the interference signals i(t) required for measuring the sweep frequency width Δf can be easily provided to the signal processing apparatus 20.
Next, a wavelength swept light measurement system 101 according to a second embodiment of the present invention will be described with reference to
As illustrated in
Note that the sweep frequency width measurement unit 23 of the signal processing apparatus 20 is similar to that in
The signal processing apparatus 20 includes a microprocessor such as a digital service unit (DSU) and a CPU, and realizes, when the microprocessor and a program stored in the storage apparatus 32 are jointly operated, various types of signal processes for calculating a profile for the wavelength swept light Lx of the wavelength swept light source X, that is, in this case, a sweep frequency width Δf.
As illustrated in
The relative frequency calculation unit 22A of the signal processing apparatus 20 illustrated in
The target extraction unit 21A of the signal processing apparatus 20 illustrated in
As illustrated in
The sweep frequency width measurement unit 23 of the signal processing apparatus 20 illustrated in
As described above, in the present embodiment, in the signal processing apparatus 20, the relative frequency calculation unit 22A calculates a difference between a frequency of the interference signals i(t) at a reference time and a frequency of the interference signals i(t) at each time, as the relative frequencies fr(t), the target extraction unit 21A extracts, from the relative frequencies fr(t), target relative frequencies frt(t) corresponding to the sweep interval Taq in which the wavelength swept light Lx is swept from a maximum frequency to a minimum frequency, and the sweep frequency width measurement unit 23 measures a difference between a maximum value and a minimum value of the obtained target relative frequencies frt(t), as the sweep frequency width Δf of the wavelength swept light Lx. Thus, unlike the first embodiment, the sweep interval Taq can be automatically identified.
Next, a wavelength swept light measurement system 102 according to a third embodiment of the present invention will be described with reference to
As illustrated in
Note that the relative frequency calculation unit 22A of the signal processing apparatus 20 is similar to that in
The signal processing apparatus 20 includes a microprocessor such as a digital service unit (DSU) and a CPU, and realizes, when the microprocessor and a program stored in the storage apparatus 32 are jointly operated, various types of signal processing units for calculating a profile for the wavelength swept light Lx of the wavelength swept light source X, that is, in this case, a sweep frequency width Δf.
As illustrated in
The relative frequency calculation unit 22A of the signal processing apparatus 20 illustrated in
The sweep frequency width measurement unit 23A is configured to first calculate the time differential value dfr(t)/dt of the relative frequency change curve obtained from the relative frequency fr(t), and to measure and output a difference between a maximum value and a minimum value of the relative frequency fr(t) identified from the time differential values dfr(t)/dt, as the sweep frequency width Δf.
As illustrated in
As described above, in the present embodiment, in the signal processing apparatus 20, the relative frequency calculation unit 22 calculates a difference between the frequency of the interference signals i(t) at the reference time and the frequency of the interference signals i(t) at each time, as the relative frequencies fr(t), and the sweep frequency width measurement unit 23A calculates time differential values dfr(t)/dt of the relative frequency change curve obtained from the relative frequencies fr(t), and measures a difference between a maximum value and a minimum value of the relative frequencies fr(t) identified from the time differential values dfr(t)/dt, as the sweep frequency width Δf. Consequently, it is possible to omit a process of obtaining the maximum value and the minimum value from the relative frequency fr(t), and thus, unlike the second embodiment, a processing amount can be reduced in the signal processing apparatus 20.
Next, a wavelength swept light measurement system 103 according to a fourth embodiment of the present invention will be described with reference to
As illustrated in
In addition to the interferometer 11 and the balanced photodetector 12 as with those illustrated in
In the first embodiment, the wavelength swept light Lx is directly input to the interferometer 11, but in the present embodiment, the wavelength swept light Lx is branched into two light beams by using a coupler C3, and the two light beams are simultaneously input to each of the interferometer 11 and the narrow band wavelength filter 13.
The specific wavelength light λbL detected by the narrow band wavelength filter 13 is photoelectrically converted by the photodetector 14 into a specific wavelength electrical signal λbE(t) and then input to the ADC 30.
The narrow band wavelength filter 13 passes light centered on the previously set wavelength λbase among the wavelength swept light Lx. Note that the wavelength λbase is within a swept wavelength band of the wavelength swept light source X. Consequently, at a time when a center of the wavelength of the wavelength swept light Lx is the wavelength λbase, the specific wavelength light λbL, the specific wavelength electrical signal λbE(t), and a specific wavelength signal λb(t) each have a maximum value.
AS in
The signal processing apparatus 20 includes a microprocessor such as a digital service unit (DSU) and a CPU, and realizes, when the microprocessor and a program stored in the storage apparatus 32 are jointly operated, various types of signal processes for calculating a profile for the wavelength swept light Lx of the wavelength swept light source X, that is, in this case, a sweep frequency width Δf.
As illustrated in
The predicted frequency calculation unit 24 of the signal processing apparatus 20 illustrated in
As illustrated in
A relative frequency fr(tbase) (=fr, base) corresponds to a frequency fr(tbase), and thus, fr(0)(=f0) is expressed as in Equation (15) below. Furthermore, the center frequency of the wavelength swept light Lx is expressed as in Equation (16) below.
According to this scheme, the predicted frequency calculation unit 24 calculates a predicted (center) frequency, that is, the predicted frequency fp(t) of the wavelength swept light Lx in accordance with Equation (17) below.
The predicted wavelength calculation unit 25 of the signal processing apparatus 20 illustrated in
The predicted wavelength calculation unit 25 converts the predicted frequency fp(t) into the predicted wavelength λp(t) in accordance with Equation (18) below.
The swept wavelength width measurement unit 26 of the signal processing apparatus 20 illustrated in
As described above, in the present embodiment, the photoelectric conversion apparatus 10 photoelectrically converts and outputs each of the interference light iL obtained as a result of interference of the wavelength swept light Lx output from the wavelength swept light source X in the interferometer 11, and a specific wavelength light λbL detected by the narrow band wavelength filter 13 from the wavelength swept light Lx, and the signal processing apparatus 20 calculates, in chronological order, the relative frequencies fr(t) for interference signals i(t) obtained by photoelectric conversion of the interference light iL, calculates the predicted frequencies fp(t) for the relative frequencies fr(t), on the basis of the detection timing of the specific wavelength light λbL obtained by photoelectric conversion of the specific wavelength light λbL, and measures a difference between the maximum value and the minimum value of the predicted wavelengths λp(t) corresponding to the obtained predicted frequencies fp(t), as the swept wavelength width Δλ of the wavelength swept light Lx.
Consequently, it is possible to measure the swept wavelength width Δλ of the wavelength swept light Lx with extremely high precision and accuracy, regardless of a sweeping speed of the wavelength swept light source X.
Furthermore, in the present embodiment, in the signal processing apparatus 20, the predicted frequency calculation unit 24 may calculate the predicted frequency fp(t), at each time t, of the relative frequency fr(t) calculated by the relative frequency calculation unit 22, on the basis of the detection timing of the specific wavelength light λbL, the predicted wavelength calculation unit 25 may calculate the predicted wavelength λp(t) of the predicted frequency fp(t) at each time t, and the swept wavelength width measurement unit 26 may measure a difference between a maximum value and a minimum value of the predicted wavelength λp(t) as the swept wavelength width Δλ of the wavelength swept light Lx. Consequently, the swept wavelength width Δλ of the wavelength swept light Lx can be measured by a relatively simple process.
Furthermore, in the present embodiment, in the signal processing apparatus 20, the predicted frequency calculation unit 24 may calculate the predicted frequency fp(t), at each time t, of the relative frequency fr(t) calculated by the relative frequency calculation unit 22, on the basis of the detection timing of the specific wavelength light λbL, and the swept wavelength width measurement unit 26A may calculate only two predicted wavelengths λp(t) corresponding to a maximum value and a maximum value of the predicted frequencies fp(t), and measure a difference between the two predicted wavelengths λp(t) as the swept wavelength width Δλ of the wavelength swept light Lx. Consequently, it is possible to reduce a processing load for calculating the predicted wavelength λp(t) in the signal processing apparatus 20.
Furthermore, in the present embodiment, in the signal processing apparatus 20, the predicted frequency calculation unit 24 may calculate only two predicted frequencies fp(t) corresponding to a maximum value and a maximum value of the relative frequencies fr(t) calculated by the relative frequency calculation unit 22, on the basis of the detection timing of the specific wavelength light λbL, and the swept wavelength width measurement unit 26B may calculate only two predicted wavelengths λp(t) corresponding to a maximum value and a maximum value of the predicted frequencies fp(t), and measure a difference between the two predicted wavelengths λp(t), as the swept wavelength width Δλ of the wavelength swept light Lx. Consequently, it is possible to reduce a processing load for calculating the predicted frequency fp(t) and the predicted wavelength λp(t) in the signal processing apparatus 20.
The present embodiment is described on the assumption that the specific wavelength light λbL, the specific wavelength electrical signal λbE(t), and the specific wavelength signal λb(t) have a peak at a time when the center of the wavelength of the wavelength swept light Lx output from the wavelength swept light source X to be measured is λbase. However, in the ADC 30, the peak of the specific wavelength electrical signal λbE(t) may not always be sampled. That is, the peak time tbase of the specific wavelength electrical signal λbE(t) may not be synchronized with a sampling time of the ADC 30, but may be between two adjacent sampling times.
Consequently, for example, the predicted frequency calculation unit 24 illustrated in
Next, a wavelength swept light measurement system 104 according to a fifth embodiment of the present invention will be described with reference to
As illustrated in
In addition to the interferometer 11, the balanced photodetector 12, the narrow band wavelength filter 13, and the photodetector 14 as with the configuration illustrated in
The wavelength swept light Lx is branched into two light beams by a coupler C4, and one of the branched light beams is further branched by the coupler C3, as with the configuration illustrated in
As in
The signal processing apparatus 20 includes a microprocessor such as a digital service unit (DSU) and a CPU, and realizes, when the microprocessor and a program stored in the storage apparatus 32 are jointly operated, various types of signal processes for calculating a profile for the wavelength swept light Lx of the wavelength swept light source X, that is, in this case, a sweep frequency width Δf.
In each of the embodiments described above, for ease of understanding, a case is described where the interference signal i(t) is constant over time, however, the interference signal i(t) is actually not constant over time, but varies with respect to time, as illustrated in
In depth distance measurement methods such as an FMCW radar method, an SS-OCT method, and an OFDR method, a resolution of the depth distance is inversely proportional to the sweep frequency width Δf of the wavelength swept light Lx, and substantially depends on a shape of the optical spectrum of the wavelength swept light Lx. The following description states that the width of the optical spectrum of the wavelength swept light Lx has an important influence on the performance of depth distance measurement methods such as the FMCW radar method, the SS-OCT method, and the OFDR method.
According to the reference literature mentioned above, when a coherence length does not differ for each frequency of the light, an envelope of the interference signal i(t) coincides with a shape of the light intensity p(t).
In the depth distance measurement methods such as the FMCW radar method, the SS-OCT method, and the OFDR method, the rescaled interference signal i′(t) is Fourier-transformed to obtain depth information.
The depth information has one peak with respect to one reflection surface. A signal for the one reflection surface is called a point spread function (PSF). The interference signal i′(t) generated by one reflection surface is an AM-modulated signal obtained by resealing the light intensity p(t) to have a sine wave having a frequency proportional to a magnitude of a position z on the reflection surface. The PSF value is a result of the Fourier transform of the interference signal i′(t), and thus, the shape of the point spread function PSF coincides with the shape of the frequency obtained by Fourier transform of the interference signal i′(t).
A distance L on the horizontal axis of the graph of the point spread function PSF illustrated in
where c is a speed of light.
A time axis t′ of the resealed interference signal i′(t) is proportional to the frequency of the light, and thus, the envelope of the interference signal i(t) coincides with the shape of the optical spectrum of the wavelength swept light Lx. That is, it can be understood that a result of the Fourier transform of a spectral shape of the light coincides with the shape of the point spread function PSF. Thus, as the optical spectrum is wider, a width of the point spread function PSF is narrower, and thus, the resolution of the depth information increases.
The sweep frequency width Δf of the wavelength swept light source X and the resolution of the depth information (distance) will be also described below. Discrete Fourier transform is usually used for the Fourier transform of the interference signal i(t) in determining the depth information, and frequencies of the discrete Fourier transform are expressed by an integer, and an interval (that is, “11”) between the frequencies is 1/Tsw at an actual frequency. Thus, the depth information is expressed in increments of a frequency 1/Tsw. n appearing in the equation z=n*Δz=n*c/(2*Δf) expressed above represents a frequency of the discrete Fourier transform, and the distance z is expressed in increments of c/(2*Δf). This indicates that, as the sweep frequency width Δf of the wavelength swept light source X is wider, the resolution of the depth information (distance) increases in inverse proportion to the sweep frequency width Δf.
As illustrated in
The frequency-to-light intensity calculation unit 27 is configured to calculate a frequency-to-light intensity spf(f), that is, a frequency spectrum, indicating the light intensity of the wavelength swept light Lx at each frequency of the wavelength swept light Lx, on the basis of the predicted frequency fp(t) from the predicted frequency calculation unit 24 and the light intensity p(t) from the ADC 30. In the frequency-to-light intensity calculation unit 27, the frequency-to-light intensity spf(f) may be output from the signal processing apparatus 20, in accordance with a request from a user.
The frequency-to-light intensity spf(f) is a value representing the light intensity for the wavelength swept light Lx at each frequency, and thus, may be regarded as the frequency spectrum of the wavelength swept light Lx. In a method for obtaining the frequency-to-light intensity spf(f), it may be possible to combine the predicted frequency fp(t) and the light intensity p(t) at the same time t to use the light intensity p(t) for the predicted frequency fp(t) as the frequency-to-light intensity spf(f), for example. In another method, it may be possible to determine an inverse function predicted frequency t(f) of the predicted frequency fp(t) and use p(t(f)) as the frequency-to-light intensity spf(f). In the both methods, the frequency-to-light intensity spf(f) is discrete data obtained from sampling by the ADC 30, and thus, may be interpolated, as necessary.
The sweep frequency width measurement unit 23B is configured to measure and output the sweep frequency width Δf of the wavelength swept light Lx from the frequency-to-light intensity spf(f). For example, it may be possible to determine the sweep frequency width Δf as the full width at half maximum or as a width that is 1/e2 of the maximum value of the frequency-to-light intensity spf(f).
As described above, in the present embodiment, the photoelectric conversion apparatus 10 photoelectrically converts each of the wavelength swept light Lx and the specific wavelength light λbL detected by the narrow band wavelength filter 13 from the wavelength swept light Lx, and the signal processing apparatus 20 calculates the frequency-to-light intensity spf(f) indicating the light intensity of the wavelength swept light Lx for each frequency of the wavelength swept light Lx, on the basis of the predicted frequency fp(t) obtained in much the same way as the above description and the detection timing of the specific wavelength light λbL detected by the photoelectric conversion apparatus 10, and measures the sweep frequency width Δf of the wavelength swept light Lx, on the basis of the frequency-to-light intensity spf(f).
Consequently, it is possible to measure the sweep frequency width Δf of the wavelength swept light source X for the resolution of depth distance measurement methods such as the FMCW radar method, the SS-OCT method, and the OFDR method. Furthermore, when the frequency-to-light intensity spf(f) is output to the outside of the signal processing apparatus 20, it is also possible to obtain a frequency spectrum for the wavelength swept light Lx of the wavelength swept light source X.
Next, a wavelength swept light measurement system 105 according to a sixth embodiment of the present invention will be described with reference to
As illustrated in
As illustrated in
The specific wavelength electrical signal λbE(t) from the photoelectric conversion apparatus 10 is omitted from the configuration illustrated in
The signal processing apparatus 20 includes a microprocessor such as a digital service unit (DSU) and a CPU, and realizes, when the microprocessor and a program stored in the storage apparatus 32 are jointly operated, various types of signal processes for calculating a profile for the wavelength swept light Lx of the wavelength swept light source X, that is, in this case, a sweep frequency width Δf.
As illustrated in
The frequency-to-light intensity calculation unit 27A is configured to calculate a frequency-to-light intensity spf(fr), that is, a frequency spectrum, indicating the light intensity of the wavelength swept light Lx for each relative frequency fr(t) of the wavelength swept light Lx, on the basis of the relative frequency fr(t) from the relative frequency calculation unit 22 and the light intensity p(t) from the ADC 30. In the frequency-to-light intensity calculation unit 27A, the frequency-to-light intensity spf(fr) may be output from the signal processing apparatus 20, in accordance with a request from a user.
As described above, in the present embodiment, the signal processing apparatus 20 calculates, on the basis of the relative frequencies fr(t) obtained as in
Consequently, the sweep frequency width Δf of the wavelength swept light source X for the resolution of depth distance measurement methods such as the FMCW radar method, the SS-OCT method, and the OFDR method, can be measured by a simpler configuration than in
Next, a wavelength swept light measurement system 106 according to a seventh embodiment of the present invention will be described with reference to
As illustrated in
Note that, in the present embodiment, compared to the fifth embodiment in
Note that the photoelectric conversion apparatus 10, the ADC 30, and the storage apparatus 32 in the wavelength swept light measurement system 106 according to the present embodiment are similar to those in the fifth embodiment, and thus, detailed description thereof will be omitted here.
The signal processing apparatus 20 includes a microprocessor such as a digital service unit (DSU) and a CPU, and realizes, when the microprocessor and a program stored in the storage apparatus 32 are jointly operated, various types of signal processes for calculating a profile for the wavelength swept light Lx of the wavelength swept light source X, that is, in this case, a swept wavelength width Δλ.
As illustrated in
The wavelength-to-light intensity calculation unit 28 is configured to calculate the wavelength-to-light intensity spec(λ), that is, a wavelength spectrum, indicating the light intensity of the wavelength swept light Lx for individual wavelengths, on the basis of the predicted wavelength λp(t) from the predicted wavelength calculation unit 25 and the light intensity p(t) from the ADC 30. In the wavelength-to-light intensity calculation unit 28, the wavelength-to-light intensity spec(λ) may be output from the signal processing apparatus 20, in accordance with a request from a user.
The swept wavelength width measurement unit 26A is configured to measure a difference between a maximum value and a minimum value of the wavelength-to-light intensity spec(λ) output from the wavelength-to-light intensity calculation unit 28, as the swept wavelength width Δλ.
As described above, in the present embodiment, the signal processing apparatus 20 calculates the wavelength-to-light intensity spec(λ) indicating the light intensity of the wavelength swept light Lx for the predicted wavelength λp(t), on the basis of the predicted wavelength λp(t) obtained as in
Consequently, it is possible to measure the swept wavelength width Δλ of the wavelength swept light Lx with extremely high precision and accuracy, regardless of a sweeping speed of the wavelength swept light source X. Furthermore, when the wavelength-to-light intensity spec(λ) is output to the outside of the signal processing apparatus 20, it is also possible to obtain a wavelength spectrum for the wavelength swept light Lx of the wavelength swept light source X.
Next, a wavelength swept light measurement system 107 according to an eighth embodiment of the present invention will be described with reference to
As illustrated in
The signal processing apparatus 20 includes a microprocessor such as a digital service unit (DSU) and a CPU, and realizes, when the microprocessor and a program stored in the storage apparatus 32 are jointly operated, various types of signal processes for calculating a profile for the wavelength swept light Lx of the wavelength swept light source X, that is, in this case, a sweep frequency width Δf.
As illustrated in
The relative frequency calculation unit 22B is configured to calculate the relative frequency fr(t) and the light intensity pp(t) of the wavelength swept light Lx, on the basis of the target interference signal it(t) output from the target extraction unit 21.
Among these processing units, the negative frequency component deletion unit 41, the argument calculation unit 42, and the argument-to-frequency conversion unit 43 are similar to those in
The light intensity calculation unit 44 calculates an intensity of a signal i′(t) output from the negative frequency component deletion unit 41 and outputs a result of the calculation as the light intensity pp(t) of the wavelength swept light Lx. The signal i′(t) is a complex number, and the light intensity pp(t) can be calculated, for example, by Equation (21) below.
Math. 21
p
p(t)=√
where R(x) and I(x) represent a real part and an imaginary part of x, respectively.
The frequency-to-light intensity calculation unit 27B is configured to calculate the frequency-to-light intensity spf(fr), that is, a frequency spectrum, indicating the light intensity of the wavelength swept light Lx for the relative frequency fr(t), on the basis of the predicted frequency fp(t) from the predicted frequency calculation unit 24 and the light intensity pp(t) output from the relative frequency calculation unit 22B.
As described above, in the present embodiment, the signal processing apparatus 20 calculates, on the basis of the detection timing of the specific wavelength light λbL detected by the photoelectric conversion apparatus 10 and the relative frequencies fr(t), the predicted frequency fp(t) indicating an absolute frequency for these relative frequencies fr(t), calculates the light intensity pp(t) of the wavelength swept light Lx from the target interference signal it(t), calculates the frequency-to-light intensity spf(f) indicating the light intensity of the wavelength swept light Lx for the predicted frequency fp(t), on the basis of the light intensity pp(t) and the predicted frequency fp(t), and measures the sweep frequency width Δf, on the basis of the frequency-to-light intensity spf(f).
Consequently, compared to the fifth embodiment of
Next, a wavelength swept light measurement system 108 according to a ninth embodiment of the present invention will be described with gathering
In the eighth embodiment of
The signal processing apparatus 20 includes a microprocessor such as a digital service unit (DSU) and a CPU, and realizes, when the microprocessor and a program stored in the storage apparatus 32 are jointly operated, various types of signal processes for calculating a profile for the wavelength swept light Lx of the wavelength swept light source X, that is, in this case, a sweep frequency width Δf.
As illustrated in
The frequency-to-light intensity calculation unit 27C is configured to calculate the frequency-to-light intensity spf(fr), that is, a frequency spectrum, indicating the light intensity of the wavelength swept light Lx for the relative frequency fr(t), on the basis of the relative frequency fr(t) and the light intensity pp(t) of the wavelength swept light Lx output from the relative frequency calculation unit 22B.
As described above, in the present embodiment, the signal processing apparatus 20 calculates the light intensity pp(t) of the wavelength swept light Lx from the target interference signal it(t), calculates the frequency-to-light intensity spf(fr) indicating the light intensity of the wavelength swept light Lx for the relative frequency fr(t), on the basis of the light intensity pp(t) and the relative frequency fr(t), and measures the sweep frequency width Δf on the basis of the frequency-to-light intensity spf(fr).
Consequently, compared to the sixth embodiment of
Next, a wavelength swept light measurement system 109 according to a tenth embodiment of the present invention will be described with reference to
As illustrated in
The signal processing apparatus 20 includes a microprocessor such as a digital service unit (DSU) and a CPU, and realizes, when the microprocessor and a program stored in the storage apparatus 32 are jointly operated, various types of signal processes for calculating a profile for the wavelength swept light Lx of the wavelength swept light source X, that is, in this case, a swept wavelength width Δλ.
As illustrated in
The wavelength-to-light intensity calculation unit 28A is configured to calculate the wavelength-to-light intensity spec(λ), that is, a wavelength spectrum, indicating the light intensity of the wavelength swept light Lx for individual wavelengths, on the basis of the predicted wavelength λp(t) from the predicted wavelength calculation unit 25 and the light intensity pp(t) from the relative frequency calculation unit 22B. In the wavelength-to-light intensity calculation unit 28A, the wavelength-to-light intensity spec(λ) may be output from the signal processing apparatus 20, in accordance with a request from a user.
As described above, in the present embodiment, the signal processing apparatus 20 calculates the wavelength-to-light intensity spec(λ) indicating the light intensity of the wavelength swept light Lx for the predicted wavelength λp(t), on the basis of the predicted wavelength λp(t) obtained as in
Consequently, it is possible to measure the swept wavelength width Δλ of the wavelength swept light Lx with extremely high precision and accuracy, regardless of a sweeping speed of the wavelength swept light source X. Furthermore, when the wavelength-to-light intensity spec(λ) is output to the outside of the signal processing apparatus 20, it is also possible to obtain a wavelength spectrum for the wavelength swept light Lx of the wavelength swept light source X.
The present invention has been described above with reference to the embodiments, but the present invention is not limited to the above-described embodiments. Various changes understood by a person skilled in the art within the scope of the present invention can be made to the configurations and details of the present invention. Furthermore, the embodiments can be freely combined within a range where no inconsistency occurs.
100, 101, 102, 103, 104, 105, 106, 107, 108, 109 . . . Wavelength swept light measurement system
10 . . . Photoelectric conversion apparatus
11 . . . Interferometer
12 . . . Balanced photodetector
13 . . . Narrow band wavelength filter
14, 15 . . . Photodetector
C1, C2, C3, C4 . . . Coupler
20 . . . Signal processing apparatus
21, 21A . . . Target extraction unit
22, 22A, 22B . . . Relative frequency calculation unit
23, 23A, 23B . . . Sweep frequency width measurement unit
24 . . . Predicted frequency calculation unit
25 . . . Predicted wavelength calculation unit
26, 26A, 26B . . . Swept wavelength width measurement unit
27, 27A, 27B, 27C . . . Frequency-to-light intensity calculation unit
28, 28A . . . Wavelength-to-light intensity calculation unit
30 . . . A/D converter (ADC)
31 . . . Memory
32 . . . Storage apparatus
41 . . . Negative frequency component deletion unit
42 . . . Argument calculation unit
43 . . . Argument-to-frequency conversion unit
44 . . . Light intensity calculation unit
45 . . . Fourier transform unit
46 . . . Negative frequency component replacement unit
47 . . . Inverse Fourier transform unit
X . . . Wavelength swept light source
Lx . . . Wavelength swept light
Trg, trE(t) . . . Trigger electrical signal
iL . . . Interference light
iE(t) . . . Interference electrical signal
λbL . . . Specific wavelength light
λbE(t) . . . Specific wavelength electrical signal
pE(t) . . . Light intensity electrical signal
tr(t) . . . Trigger signal
i(t), i′(t) . . . Interference signal
it(t) . . . Target interference signal
λb(t) . . . Specific wavelength signal
p(t), pp(t) . . . Light intensity
fr(t) . . . Relative frequency
frt(t) . . . Target relative frequency
fp(t) . . . Predicted frequency
λp(t) . . . Predicted wavelength
spf(f), spf(fr) . . . Frequency-to-light intensity
spec(λ) . . . Wavelength-to-light intensity
Δf . . . Sweep frequency width
Δλ . . . Swept Wavelength width
Tmem . . . Effective period
Taq . . . Sweep interval
Ttrg . . . Trigger time
Tpre . . . Pre-time
Tpos . . . Post-time.
This application is a national phase entry of PCT Application No. PCT/JP2019/028300, filed on Jul. 18, 2019, which application is hereby incorporated herein by reference.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2019/028300 | 7/18/2019 | WO | 00 |