WO2018105549A1 - 波長掃引光源、波長掃引光源のための駆動データ作成方法および光偏向器 - Google Patents

波長掃引光源、波長掃引光源のための駆動データ作成方法および光偏向器 Download PDF

Info

Publication number
WO2018105549A1
WO2018105549A1 PCT/JP2017/043446 JP2017043446W WO2018105549A1 WO 2018105549 A1 WO2018105549 A1 WO 2018105549A1 JP 2017043446 W JP2017043446 W JP 2017043446W WO 2018105549 A1 WO2018105549 A1 WO 2018105549A1
Authority
WO
WIPO (PCT)
Prior art keywords
wavelength
optical
light source
soa
gain
Prior art date
Application number
PCT/JP2017/043446
Other languages
English (en)
French (fr)
Inventor
上野 雅浩
明晨 陳
詔子 辰己
坂本 尊
雄三 佐々木
豊田 誠治
勇一 赤毛
山口 城治
匡 阪本
Original Assignee
日本電信電話株式会社
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 日本電信電話株式会社 filed Critical 日本電信電話株式会社
Priority to EP17878420.3A priority Critical patent/EP3553903B1/en
Priority to US16/466,608 priority patent/US11165219B2/en
Priority to EP23153071.8A priority patent/EP4191808A1/en
Priority to JP2018554986A priority patent/JP6774504B2/ja
Publication of WO2018105549A1 publication Critical patent/WO2018105549A1/ja
Priority to US17/380,309 priority patent/US11721948B2/en

Links

Images

Classifications

    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01SDEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
    • H01S3/00Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
    • H01S3/10Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating
    • H01S3/13Stabilisation of laser output parameters, e.g. frequency or amplitude
    • H01S3/131Stabilisation of laser output parameters, e.g. frequency or amplitude by controlling the active medium, e.g. by controlling the processes or apparatus for excitation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J9/00Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01SDEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
    • H01S3/00Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
    • H01S3/05Construction or shape of optical resonators; Accommodation of active medium therein; Shape of active medium
    • H01S3/08Construction or shape of optical resonators or components thereof
    • H01S3/08004Construction or shape of optical resonators or components thereof incorporating a dispersive element, e.g. a prism for wavelength selection
    • H01S3/08009Construction or shape of optical resonators or components thereof incorporating a dispersive element, e.g. a prism for wavelength selection using a diffraction grating
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01SDEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
    • H01S3/00Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
    • H01S3/10Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating
    • H01S3/10007Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating in optical amplifiers
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01SDEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
    • H01S3/00Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
    • H01S3/10Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating
    • H01S3/10069Memorized or pre-programmed characteristics, e.g. look-up table [LUT]
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01SDEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
    • H01S3/00Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
    • H01S3/10Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating
    • H01S3/106Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating by controlling devices placed within the cavity
    • H01S3/107Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating by controlling devices placed within the cavity using electro-optic devices, e.g. exhibiting Pockels or Kerr effect
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01SDEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
    • H01S3/00Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
    • H01S3/10Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating
    • H01S3/13Stabilisation of laser output parameters, e.g. frequency or amplitude
    • H01S3/136Stabilisation of laser output parameters, e.g. frequency or amplitude by controlling devices placed within the cavity
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01SDEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
    • H01S5/00Semiconductor lasers
    • H01S5/06Arrangements for controlling the laser output parameters, e.g. by operating on the active medium
    • H01S5/0617Arrangements for controlling the laser output parameters, e.g. by operating on the active medium using memorised or pre-programmed laser characteristics
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01SDEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
    • H01S5/00Semiconductor lasers
    • H01S5/10Construction or shape of the optical resonator, e.g. extended or external cavity, coupled cavities, bent-guide, varying width, thickness or composition of the active region
    • H01S5/14External cavity lasers
    • H01S5/141External cavity lasers using a wavelength selective device, e.g. a grating or etalon
    • H01S5/143Littman-Metcalf configuration, e.g. laser - grating - mirror
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02001Interferometers characterised by controlling or generating intrinsic radiation properties
    • G01B9/02002Interferometers characterised by controlling or generating intrinsic radiation properties using two or more frequencies
    • G01B9/02004Interferometers characterised by controlling or generating intrinsic radiation properties using two or more frequencies using frequency scans
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/0209Low-coherence interferometers
    • G01B9/02091Tomographic interferometers, e.g. based on optical coherence
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01SDEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
    • H01S5/00Semiconductor lasers
    • H01S5/02Structural details or components not essential to laser action
    • H01S5/022Mountings; Housings
    • H01S5/0225Out-coupling of light
    • H01S5/02251Out-coupling of light using optical fibres
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01SDEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
    • H01S5/00Semiconductor lasers
    • H01S5/06Arrangements for controlling the laser output parameters, e.g. by operating on the active medium
    • H01S5/062Arrangements for controlling the laser output parameters, e.g. by operating on the active medium by varying the potential of the electrodes
    • H01S5/06209Arrangements for controlling the laser output parameters, e.g. by operating on the active medium by varying the potential of the electrodes in single-section lasers
    • H01S5/06213Amplitude modulation
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01SDEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
    • H01S5/00Semiconductor lasers
    • H01S5/06Arrangements for controlling the laser output parameters, e.g. by operating on the active medium
    • H01S5/068Stabilisation of laser output parameters
    • H01S5/0683Stabilisation of laser output parameters by monitoring the optical output parameters
    • H01S5/0687Stabilising the frequency of the laser
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01SDEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
    • H01S5/00Semiconductor lasers
    • H01S5/10Construction or shape of the optical resonator, e.g. extended or external cavity, coupled cavities, bent-guide, varying width, thickness or composition of the active region
    • H01S5/14External cavity lasers
    • H01S5/141External cavity lasers using a wavelength selective device, e.g. a grating or etalon

Definitions

  • the present invention relates to a wavelength swept light source. More specifically, the present invention relates to a wavelength swept light source that can be used in an optical coherence tomography apparatus.
  • Wavelength swept light sources are widely used in optical equipment and electronic equipment using imaging technology that uses light. For example, it is most often used in consumer electronic devices such as cameras, printers, and facsimiles, and has spread to the medical field.
  • OCT optical coherence tomography
  • OCT has been attracting attention in order to non-invasively image a tomography inside a living body. It was.
  • OCT has also been used in the field of molecular (molecular) imaging for elucidating biological events at the molecular and cellular levels in vivo.
  • the wavelength swept light source is an important key element that supports the OCT technology, and greatly affects the performance of the OCT apparatus.
  • OCT apparatuses that acquire tomographic images using coherence interference include time domain OCT (TD-OCT: Time Domain Optical Coherence Tomography) and Fourier domain OCT (FD-OCT: Fourier Domain Optical Coherence Tomography).
  • the FD-OCT is further classified into a spectral domain OCT (SD-OCT: Spectral Domain Optical Coherence Tomography) and a wavelength sweep OCT (SS-OCT: Swept-Source Optical Coherence Tomography).
  • SD-OCT Spectral Domain Optical Coherence Tomography
  • SS-OCT Swept-Source Optical Coherence Tomography
  • the wavelength swept light source is required to operate according to the application field, such as a wave number linear light source in which the swept wave number changes linearly with time.
  • Patent Document 1 shows a configuration example of a wavelength swept light source using an optical fiber loop.
  • the wavelength swept light source includes a gain medium (12) and a circulator (13) in an optical fiber loop (11).
  • Light is extracted from the circulator (13) through the collimating lens (22) into the space, and the light reaches the mirror (25) through the bandpass filter (24).
  • the reflected light from the mirror is incident again on the collimating lens through the band-pass filter (24), whereby a reciprocating optical path of the resonator is formed.
  • the bandpass filter (24) for example, an optical interference type dielectric multilayer filter can be used, and the transmission wavelength can be changed by changing the incident angle of light to the filter.
  • the galvanometer (26) By changing the angle of the bandpass filter with respect to the optical axis at high speed can be realized.
  • the circular gain (gain) and the oscillation threshold vary depending on the wavelength.
  • this wavelength bandwidth is referred to as an unsaturated gain bandwidth.
  • the wavelength swept light source for SS-OCT be maintained in a state where the coherence length is long at all wavelengths being swept. In order to realize this, it is important that the non-saturation gain bandwidth (wavelength width capable of laser oscillation) is narrow over the entire wavelength being swept.
  • the operation state is such that the peak of the circular gain spectrum matches the oscillation spectrum with respect to all sweep wavelengths.
  • the amount of current injected into the gain medium since the amount of current injected into the gain medium is constant over time, it is difficult to realize such an operating state.
  • the amount of current injected into the gain medium determines the shape and peak value of the loop gain spectrum. For this reason, in the conventional wavelength swept light source in which the injection current amount is constant over time, the circular gain peak value does not always coincide with the oscillation threshold value over the entire swept wavelength. The reason can be further explained as follows.
  • the circular gain and the oscillation threshold are determined by the gain of the gain medium, various optical losses in the laser resonator, the filtering characteristics of the bandpass filter, the time taken for oscillation for each wavelength, and the like.
  • the time required for oscillation of each wavelength and the filtering characteristics of the bandpass filter are constant regardless of the wavelength, but the gain of the gain medium is the same as the sweep wavelength.
  • the value decreases as the distance from the center wavelength increases. Therefore, the oscillation threshold value near the center wavelength is small, and the oscillation threshold value increases as the distance from the center wavelength increases.
  • the amount of current injected into the gain medium is constant for any wavelength being swept, so that a swept wavelength width necessary for a wide band (for example, about 100 nm in SS-OCT) is secured. Therefore, it is necessary to determine the amount of current injected into the gain medium in accordance with the wavelength having a large oscillation threshold. In this way, when a constant current is injected into the gain medium, the wavelength width (non-saturation gain bandwidth) of the circular gain exceeding the threshold is large at a wavelength where the oscillation threshold is small.
  • the coherence length is long near the shortest wavelength and the longest wavelength. Yes.
  • the amount of injected current becomes excessive, and the degree to which the peak value of the circular gain exceeds the oscillation threshold becomes significant. For this reason, in the central part of the sweep wavelength band, the unsaturated gain bandwidth is widened, the oscillation line width is widened, and the coherence length is short.
  • the intensity of the interference signal increases as the coherence length of light from the wavelength swept light source increases.
  • the intensity of one pixel of the OCT image obtained by frequency conversion of the interference signal that is, the point spread function (PSF)
  • PSF point spread function
  • the energy (time integration value) of the interference signal is proportional to the energy (time integration value) of the interference signal. Therefore, in order to increase the intensity of the PSF, it is necessary to strengthen the entire interference signal in the entire wavelength band in which the wavelength is swept, and the coherence length must be long in the entire wavelength band to be swept.
  • the conventional wavelength swept light source in which a constant current is injected into the gain medium, the coherence length cannot be maintained in a long state over the entire wavelength band in which the wavelength is swept. As a result, there is a problem that the PSF intensity that directly affects the quality of the output image of the SS-OCT apparatus is lowered.
  • the present invention has been made in view of such a problem, and an object thereof is to provide a light source capable of maintaining a long coherence length over a wide sweep wavelength band. It is in. Also disclosed are a method for creating semiconductor optical amplifier (SOA) current waveform data (optical gain control data) for optimally operating the wavelength-swept light source of the present invention, and a method for determining the coherence length. An optical deflector suitable for implementing the wavelength swept light source of the present invention is also disclosed.
  • SOA semiconductor optical amplifier
  • a gain medium for amplifying the light, a means for changing the oscillation wavelength, and an optical gain of the gain medium are changed.
  • a laser resonator including means, a wavelength control means for sweeping the oscillation wavelength by controlling the means for changing the oscillation wavelength by a wavelength control signal, and light in the laser resonator over the entire wavelength band of the wavelength sweep.
  • a wavelength-swept light source comprising: an optical gain control means for controlling the means for changing the optical gain so that the gain is slightly higher than the optical gain in the laser oscillation threshold state. It is.
  • the optical gain control means includes a plurality of times obtained from the start of the sweep obtained from the wavelength control means, and light having a longest coherence length of the wavelength sweep light at each of the plurality of times. Based on the optical gain control data combined with the gain, the optical gain control signal may be configured to control the means for changing the optical gain.
  • the means for changing the oscillation wavelength includes a diffraction grating, a mirror on which the oscillation light from the diffraction grating is incident, and a mirror that vertically reflects the incident oscillation light, and the incidence of the oscillation light on the diffraction grating And an optical deflector that changes the angle.
  • the optical gain control means detects an output light level from the laser resonator, and based on the detected output light level, the output light level is a predetermined output level. It may be configured to generate an optical gain control signal that becomes an optical level value, and to control the means for changing the optical gain by the optical gain control signal.
  • the means for changing the oscillation wavelength may comprise a diffraction grating and an optical deflector for changing an incident angle of the oscillation light to the diffraction grating.
  • a gain medium for amplifying the light in another embodiment, a gain medium for amplifying the light, means for changing the oscillation wavelength, and an optical gain of the gain medium are changed.
  • a laser resonator including means, a wavelength control means for sweeping the oscillation wavelength by controlling the means for changing the oscillation wavelength by a wavelength control signal, and a plurality of times from the start of sweep obtained from the wavelength control means And means for changing the optical gain by an optical gain control signal based on optical gain control data combined with an optical gain having the longest coherence length of the wavelength-swept light at each of the plurality of times. This is a wavelength-swept light source.
  • the above-described optical deflector is an optical deflector in which at least two electrodes are formed on opposite surfaces of an electro-optic crystal, and a voltage is applied to the at least two electrodes. Then, a non-uniform refractive index distribution is generated in the electro-optic crystal due to the electro-optic effect, so that an optical deflector that bends the path of light orthogonal to the electric field formed by the voltage, and the incidence of the electro-optic crystal
  • An optical element disposed on at least one of the side and the emission side, and having a light collecting characteristic spatially distributed over a range in which the light is deflected.
  • the light deflector has a condensing characteristic that changes depending on a deflection angle, and the optical element compensates for a change in at least one of a beam waist position and a beam waist diameter on an emission side over a range in which the light is deflected.
  • the optical characteristics can be spatially distributed.
  • the coherence length of the output light during the wavelength sweep period can always be maintained long by the wavelength sweep light source, the drive data creation method, and the optical deflector according to the present invention.
  • FIG. 7 is a diagram illustrating another specific configuration example of the interferometer illustrated in FIG. 6.
  • FIG. 7 is a diagram illustrating another specific configuration example of the interferometer illustrated in FIG. 6.
  • FIG. 7 is a diagram illustrating another specific configuration example of the interferometer illustrated in FIG. 6. It is the figure which showed the relationship between SOA injection current, optical output power, and coherence length.
  • It is a figure which shows the structure of the wavelength sweep light source of the 2nd Embodiment of this invention.
  • It is a figure which shows the structure of the wavelength sweep light source device containing a control circuit. It is a figure explaining the relationship between SOA injection current and an optical output level.
  • FIG. 6 is a diagram showing a configuration of an optical deflector suitable for a wavelength swept light source of Example 4.
  • FIG. 10 is a diagram showing a configuration of an optical deflector suitable for a wavelength swept light source of Example 5. It is the figure which showed the applied voltage dependence of the focal distance of Example 5. It is the figure which showed the relationship between the free-form surface mirror focal distance of Example 5, and an applied voltage.
  • FIG. 10 is a diagram showing a configuration of an optical deflector suitable for a wavelength swept light source of Example 6.
  • FIG. 10 is a diagram showing a configuration of an optical deflector suitable for a wavelength swept light source of Example 7. It is the figure which showed the relationship between the focal distance of an optical element, and the beam waist position at the time of Gaussian beam incidence. It is the figure which showed the relationship between the beam waist position and beam diameter magnification at the time of changing a focal distance.
  • a gain medium for amplifying the light a means for changing the oscillation wavelength, and an optical gain of the gain medium are changed.
  • a laser resonator including means, a wavelength control means for sweeping the oscillation wavelength by controlling the means for changing the oscillation wavelength by a wavelength control signal, and light in the laser resonator over the entire wavelength band of the wavelength sweep.
  • the wavelength-swept light source includes optical gain control means for controlling the means for changing the optical gain so that the gain is slightly higher than the optical gain in the laser oscillation threshold state.
  • optical gain control means configuration examples (first embodiment and second embodiment) by different methods are disclosed.
  • One is to obtain an SOA current that maintains the coherence length to the maximum and operate the light source based on the SOA current waveform data.
  • the other is to detect the output light level from the laser resonator and control the optical gain (SOA current) based on the detected output light level so that the output light level becomes a predetermined output light level value. To do.
  • the wavelength swept light source of the first embodiment of the present invention can maintain the coherence length of the output light in a long state over the entire swept wavelength band.
  • the gain of the gain medium is changed with time, so that the coherence length is maintained at the maximum state. That is, the gain medium gain is maintained in the vicinity of the oscillation threshold at each wavelength over the entire sweep wavelength band, and the unsaturated gain bandwidth of the gain medium is narrowed.
  • the optical gain is controlled so that the gain (optical gain) in the laser resonator is slightly higher than the gain in the laser oscillation threshold state (oscillation threshold gain or oscillation gain).
  • the laser gain control means By controlling the laser gain control means so as to lower the loop gain at the approximate center of the swept wavelength band that had excessive gain in the conventional wavelength swept light source, the oscillation line width of the output oscillation light is always reduced. Keep it narrow.
  • the circular gain of the gain medium can be changed by the SOA current flowing in a semiconductor optical amplifier (SOA) injected into the gain medium in the laser resonator.
  • SOA semiconductor optical amplifier
  • the SOA current that maintains the maximum coherence length in the wavelength swept light source of the present invention is generated based on the SOA current waveform data.
  • the coherence length corresponding to each wavelength is obtained from the interference light of the light source output light while associating with the sweep time information in a state where the wavelength swept light source is swept for different SOA currents.
  • the SOA current that gives the maximum coherence length is selected from among the different SOA currents.
  • the coherence length is obtained based on a plurality of interference lights obtained from the wavelength sweep light, and the SOA current that maintains the coherence length at the maximum is obtained.
  • the SOA current for maintaining the coherence length at the maximum is supplied from a current source capable of outputting an arbitrary current waveform, and the wavelength swept light source is a SOA current waveform for realizing the SOA current waveform stored in the storage means (memory). Operates based on data.
  • the gain change of the gain medium of the laser resonator can be realized by a gain variable means other than the SOA drive current.
  • the present invention can also be implemented as a method for obtaining SOA current waveform data for driving a wavelength swept light source while maintaining a long coherence length of output oscillation light in the entire swept wavelength band. Also disclosed is a novel method of measuring the coherence length unique to the present invention adapted to the wavelength swept light source of the present invention for determining SOA current waveform data.
  • the wavelength-swept light source of the present invention a method for creating the drive data thereof, and a method for measuring the coherence length will be described together with each example.
  • Example 1 of the wavelength swept light source of the first embodiment Example 1 of the wavelength swept light source of the first embodiment.
  • FIG. 1 is a diagram showing a configuration example of a wavelength swept light source of the present invention.
  • the wavelength swept light source 100 of the present invention will be described as an example having a Littman type laser resonator 101.
  • the present invention can be applied to any type of laser resonator as long as the resonator has a gain medium having a wavelength-dependent circular gain in the resonator.
  • the light output from the SOA 102 is collimated into parallel light by the lens 105 and deflected through the optical deflector 106.
  • the optical deflector 106 deflects the optical path in the vertical direction of the drawing as shown in FIG. 1, for example, by a drive signal 115 from an optical deflector driver described later.
  • FIG. 1 two cases of the optical path 121 in a state of being deflected to the maximum on the upper side and the optical path 122 in a state of being deflected to the maximum on the lower side are illustrated.
  • the state of the optical path 121 and the state of the optical path 122 correspond to the longest wavelength state and the shortest wavelength state, respectively.
  • the light deflected by the optical deflector 106 is diffracted by the diffraction grating 107 and reaches the mirror 108. Only the light of the wavelength perpendicularly incident on the mirror 108 follows the optical path in the direction opposite to the optical path described above, and returns to the SOA 102.
  • the light returning from the mirror 108 to the SOA 102 is reflected by the half mirror 104 between the SOA 102 and the coupler 103.
  • a reciprocating optical path is formed between the mirror 108 and the half mirror 104, and laser oscillation occurs. A part of the oscillation light passes through the half mirror 104 from the SOA 102 and is output to the outside of the laser resonator 101 as output light through the coupler 103 and the optical fiber 120.
  • the laser oscillation wavelength is changed depending on the incident angle of the light to the diffraction grating 107. be able to.
  • the configuration between the SOA 102 and the mirror 108 forms, for example, a band-pass filter having a generally quadratic function shape, and determines the non-saturation gain bandwidth of the circular gain.
  • the optical deflector driver 113 controls the deflection direction (deflection angle) of the light emitted from the optical deflector 106 by the drive signal 115.
  • the drive signal 115 changes with time, the deflection direction also changes gradually.
  • the drive signal is periodically and repeatedly changed so that the wavelength of the oscillation light is swept in a direction that smoothly increases or decreases in one direction.
  • the wavelength sweep operation is realized by changing the oscillation wavelength continuously in time by the drive signal 115.
  • a synchronizing signal 114 synchronized with the wavelength sweep is output together with the drive signal. For example, when the wavelength sweep is performed from a short wavelength to a long wavelength, the synchronization signal 114 is output to the arbitrary waveform current source 111 at a timing at which the deflection angle corresponding to the shortest wavelength is obtained.
  • Arbitrary waveform current source 111 uses SOA current waveform data for generating SOA current 116 to be applied to SOA 102 in an associated memory (storage means) not shown in FIG. I have.
  • This memory may be one that reads current waveform data from another and stores it temporarily, or one that stores it semi-permanently. Details of the contents of the current waveform data will be described in the second and third embodiments, but information for generating the SOA current waveform data is recorded in the memory.
  • the memory may be provided in the arbitrary waveform current source 111 or may be provided anywhere in the wavelength sweep light source 100 outside the arbitrary waveform current source 111.
  • the arbitrary waveform current source 111 generates an SOA current waveform from current waveform data or information recorded in the memory in synchronization with the synchronization signal 114 from the optical deflector driver 113, and generates an SOA current according to the SOA current waveform. Flow to SOA 102.
  • a temperature controller for controlling the temperature of the SOA is normally connected to the SOA 102, but is omitted here. The same applies to the SOA described later.
  • the wavelength sweep light source further includes a control unit including a CPU for controlling the entire light source, including control of the wavelength sweep operation.
  • the control unit usually includes a memory, and an electronic control device that operates using data stored in or read from the memory usually has the control unit. Therefore, it should be noted that the controller is not shown in FIG.
  • the time-waveform itself may be recorded, or any form of data for generating the time-waveform may be recorded. That is, the time-waveform data itself, the time-waveform data and the coefficient combined therewith, the time-waveform data and the bias combined therewith, or the time-waveform data combined with the coefficient and the bias are output. May be.
  • the information recorded in the memory associated with the arbitrary waveform current source 111 can be data obtained by Fourier transforming the current waveform. Also in this case, the data obtained by inverse Fourier transform of the Fourier transformed data itself, the inverse Fourier transformed data and the coefficient combined therewith, the inverse Fourier transformed data and the bias combined therewith, or the inverse Fourier transformed data and the coefficient and bias are used. A combination may be output.
  • the optical deflector driver 113 deflects the optical deflector 106. While starting the deflection operation, the arbitrary waveform current source 111 reads out the current waveform data to be applied to the SOA 102 to the memory associated with the arbitrary waveform current source 111. The optical deflector driver 113 causes the optical deflector 106 to perform a deflection operation. The arbitrary waveform current source 111 applies the SOA current 116 corresponding to the current waveform data to the SOA 102 in synchronization with the synchronization signal 114 from the optical deflector driver 111 based on the current waveform data stored in the memory. The operation principle of the present invention will be described in more detail.
  • the fundamental narrowing operation of the oscillation line width is performed by the wavelength filter composed of the diffraction grating 107 and the mirror 108, as in the prior art.
  • the circular gain is dynamically adjusted in synchronization with the sweep, and the operation of the SOA 102 is performed so that the unsaturated gain bandwidth becomes as narrow as possible at any time during the wavelength sweep. To control.
  • the SOA current 116 is set so that the non-saturation gain bandwidth is as narrow as possible with the SOA 102 in a state where the circular gain is slightly larger than the oscillation threshold at an oscillation wavelength that changes every moment.
  • the coherence length of the oscillation light is maintained in a long state at any time during one wavelength sweep period.
  • the magnitude of the circular gain varies depending on the wavelength. Therefore, if the optical gain of the SOA is dynamically changed along with the wavelength during the wavelength sweep, the coherence length can be maintained long at any time during one wavelength sweep period.
  • the wavelength every moment during the wavelength sweep period corresponds to the time from the start of the sweep, that is, the sweep time one-to-one.
  • the wavelength swept light source of the present invention operates the SOA by using the current waveform data of the time-SOA current waveform obtained in the second embodiment.
  • FIG. 2 is a diagram for explaining the operation of the SOA current in the wavelength swept light source of the present invention.
  • the horizontal axis represents the time from the sweep start time, and the time is synchronized between the upper and lower graphs.
  • the upper graph in FIG. 2 shows the time variation of the oscillation wavelength during one wavelength sweep period.
  • the oscillation wavelength is approximately 1280 nm (approximately 4.9E6 rad / m in terms of wave number), and at the end of the sweep, the oscillation wavelength is approximately 1400 nm (approximately 4.5E6 rad / m in terms of wave number).
  • the bandwidth has reached 120 nm.
  • the SOA current 116 is changed according to the present invention in one wavelength sweep period. It can be seen that the SOA current 116 is changed according to the wavelength. Specifically, the SOA current in the central portion of the sweep wavelength band is reduced from at least one SOA current value at both ends of the sweep wavelength band (corresponding to the sweep start time and sweep end time). Thereby, the circular gain of the gain medium in the central portion of the sweep wavelength band is dynamically lowered, and the SOA 102 is controlled so that the unsaturated gain bandwidth becomes as narrow as possible at any time during the wavelength sweep. As is clear from the upper diagram of FIG. 2, there is a one-to-one relationship between the sweep time and the wavelength, and the instantaneous oscillation wavelength can be specified by the time (time) from the start of the sweep. I understand.
  • the wavelength swept light source 100 of FIG. 1 of this embodiment operates in a state where time-SOA current waveform data already exists for driving the SOA while dynamically changing the SOA current over time.
  • the method and apparatus for creating the SOA current waveform data will be described in more detail in the second and third embodiments.
  • the present invention can be any as long as it can include means for changing the circular gain.
  • the present invention can also be applied to a wavelength swept light source having such a type of laser resonator. That is, the present invention can also be applied to a wavelength swept light source equipped with a Littrow external resonator, a fiber laser, a vertical cavity surface emitting laser (VCSEL), a MEMS laser, or the like.
  • the laser resonator includes means for changing the wavelength (in FIG. 1, the optical deflector 106, the diffraction grating 107 and the mirror 108) and means for changing the optical gain (the SOA 102 in FIG. 1),
  • the present invention can be applied. Therefore, generalization is possible as follows.
  • the means for changing the wavelength receives the wavelength control signal (deflector drive signal 115) output from the wavelength control means corresponding to the optical deflector driver 113 in FIG. 1, and controls the laser oscillation wavelength. If the wavelength control unit controls the wavelength to change continuously with time, wavelength sweeping is possible. As already described, in the configuration example of the wavelength swept light source of the present embodiment in FIG. 1, the oscillation wavelength sweep operation occurs by continuously changing the light deflection angle of the optical deflector 106.
  • a KTN deflector using a potassium tantalate niobate (KTN) crystal as an electro-optic material can be used.
  • KTN potassium tantalate niobate
  • the means for changing the optical gain is controlled in response to the optical gain control signal (SOA current 116) output from the optical gain control means corresponding to the arbitrary waveform current source 111 in FIG. 1, and the optical gain is controlled. .
  • the optical gain control means changes the optical gain control signal with time
  • the optical gain also changes with time.
  • the optical gain of the SOA 102 circular gain of the gain medium
  • the SOA current waveform data in the wavelength swept light source 100 shown in FIG. 1 is data of time transition of the SOA current 116 during one sweep (data combining time and SOA current).
  • the SOA current waveform data is the time transition data (time and light) of the optical gain control signal during one wavelength sweep. It can be said to be an example of data combining gain control signals: optical gain control data.
  • the optical gain control data in which the optical gain control signal value is dynamically controlled along with the wavelength can be obtained by replacing the vertical axis of the lower graph with the optical gain control signal value.
  • the wavelength swept light source of the present invention is a wavelength swept light source that outputs a swept wavelength light whose oscillation wavelength continuously changes, a gain medium that amplifies the light, means for changing the oscillation wavelength (106, 107, 108), and A laser resonator including means (102) for changing the optical gain of the gain medium and a wavelength control means (113) for sweeping the oscillation wavelength by controlling the means for changing the oscillation wavelength by a wavelength control signal (115). ) And a plurality of times from the start of the sweep obtained from the wavelength control means, and an optical gain control data combining the optical gain with the longest coherence length of the wavelength swept light at each of the plurality of times.
  • the optical gain control signal (116) can be implemented as an optical gain control means (111) for controlling the means for changing the optical gain.
  • the means for changing the oscillation wavelength includes a diffraction grating (107), a mirror (108) on which the oscillation light from the diffraction grating is incident and vertically reflecting the incident oscillation light, and the diffraction grating And an optical deflector (106) for changing the angle of incidence of the oscillation light on.
  • the means for changing the optical gain of the gain medium is a semiconductor optical amplifier (SOA)
  • the optical gain control signal is an SOA current for driving the SOA
  • the optical gain control data is obtained at the plurality of times. It can be implemented as SOA current waveform data representing the selected SOA current for each.
  • the present invention provides optical gain control data.
  • the present invention can also be applied to a wavelength swept light source that operates based on the above. The same applies to the embodiments of the present invention described hereinafter.
  • the description has been made on the assumption that the SOA current waveform data for operating the wavelength swept light source 100 of the present invention already exists. However, in the next embodiment, this SOA current waveform data is generated. The apparatus and a series of procedures will be described.
  • FIG. 3 is a diagram showing a configuration of an apparatus for generating SOA current waveform data to be stored in a memory associated with the wavelength swept light source of the present invention.
  • the apparatus for generating the SOA current waveform data in FIG. 3 includes all the elements of the wavelength swept light source 101 of the first embodiment, and in addition, includes an element for acquiring the SOA current waveform data.
  • the additional elements in FIG. 3 represent functional blocks that indicate processing processes that can be realized by software arithmetic processing by a CPU or the like, except for the interferometer, as will be described later in the third embodiment. Therefore, it should be noted that the additional elements in FIG. 3 can be implemented as software arithmetic processing without any special hardware elements other than the processor, except for the interferometer.
  • the control unit 315 in FIG. 3 includes a processor such as a CPU and hardware elements such as an associated memory.
  • the current-time-coherence length correspondence data acquisition unit 313, the maximum coherence length correspondence current selection unit 314, and the like conveniently use elements of a series of procedures for acquiring SOA current waveform data, which will be described later together with the flowchart of FIG. Can be understood as shown. Therefore, the current-time-coherence length correspondence data acquisition unit 313 and the maximum coherence length correspondence current selection unit 314 are depicted as separate units for convenience, but a series of SOA currents to be described later under the control of the control unit 315. It should be understood that the procedure for acquiring waveform data (FIG. 4) is performed.
  • connection relation, inclusion relation, arrows, etc. between the control unit 315, the current-time-coherence length correspondence data acquisition means 313 and the maximum coherence length correspondence current selection means 314 are exchanges of control signals, control data, calculation result data, etc. Note that this is a rough conceptual representation.
  • the time-coherence length correspondence data acquisition means 312 included in the current-time-coherence length correspondence data acquisition means 313 includes a procedure for acquiring the coherence length of the oscillation output light, particularly as one explicit functional block. It is expressed as In the wavelength swept light source of the present invention, it is necessary to acquire the coherence length of the oscillation output light in order to determine the SOA current value at each wavelength when the SOA current is dynamically changed.
  • the wavelength swept light source of the present invention includes a novel aspect specific to the method and configuration for obtaining the coherence length as a premise for acquiring the SOA current waveform data, and this is disclosed as Example 3. Therefore, the configuration of the time-coherence length correspondence data acquisition unit 312 and the operation / procedure for obtaining the coherence length will be described in detail in the third embodiment.
  • the arbitrary waveform current source 111 operates as a current source that allows a constant current to flow through the SOA 102 without changing with time.
  • the current source 111 is used in the following description.
  • the SOA current value in the current source 111 can be set from the outside by SOA current set value data.
  • the current-time-coherence length correspondence data acquisition unit 313 is set, but this is an example, and the SOA current set value data may be set from the control unit 315.
  • the current-time-coherence length correspondence data acquisition unit 313 inputs SOA current value data to the current source 111, sets a predetermined SOA current value in the SOA 102, and sets time-coherence length correspondence data for each set SOA current value. get. That is, for each set SOA current value, a data set (time, coherence length) of “time” and “coherence length” is obtained.
  • time here refers to the time (time) from the sweep disclosure time in one wavelength sweep period.
  • the acquisition of the coherence length under the condition of the set specific SOA current value and the specific time is performed by the time-coherence length correspondence data acquisition means 312.
  • the time-coherence length correspondence data acquisition unit 312 simply represents a functional block that measures the coherence length under a specific condition.
  • the time-coherence length correspondence data acquisition unit 312 outputs time-coherence length correspondence data 316, that is, a data set (time, coherence length) 316.
  • the current-time-coherence length correspondence data acquisition unit 313 acquires a data set (time, coherence length) 316 for each set SOA current value, and finally “SOA current value”, “time”, “coherence length”. ”Is obtained (SOA current value, time, coherence length). Details of the above-described data set can be easily understood together with a flowchart of FIG. 4 and an example of acquired data of FIG.
  • the time-coherence length correspondence data acquisition unit 312 obtains the wavelength swept light from the laser resonator 101 in synchronization with the synchronization signal 318 from the optical deflector driver 113, and the wavelength sweep from the laser resonator 101. Measure the coherence length for each “time” of light. At this time, the time-coherence length correspondence data acquisition unit 312 generates different interference light under different interference conditions from the wavelength swept light in a state where the SOA 102 is set to one SOA current value, so that one SOA set current value is obtained. A coherence length is calculated for each “time” from a plurality of interference light data.
  • the SOA current waveform data for the wavelength swept light source of the present invention is obtained from “wavelength swept light” that is performing a wavelength sweep operation.
  • the wavelength sweep light source is not set to a specific oscillation wavelength and is operated statically to obtain the coherence length, but the wavelength sweep is based on the wavelength sweep light obtained under the actual wavelength sweep operation conditions.
  • the coherence length for each “time” of the period is calculated. As shown in the upper graph of FIG. 2, each “time” of the wavelength sweep period has a one-to-one correspondence with each wavelength in the sweep wavelength band. Therefore, the time-coherence length correspondence data acquisition unit 312 obtains the coherence length for each wavelength from the wavelength sweep light in the “actual wavelength sweep operation state” in which the oscillation wavelength dynamically changes every moment in the wavelength sweep light source. become.
  • the maximum coherence length correspondence current selection means 314 obtains the current-time-coherence length correspondence data 316 obtained by the current-time-coherence length correspondence data acquisition means 313, and for each “time” of the wavelength sweep period, The SOA current value that maximizes the coherence length is selected and determined. In other words, among three types of data sets (SOA current value, time, coherence length) of “SOA current value”, “time”, and “coherence length”, coherence is selected from data having the same “time”. The SOA current value having the maximum length is selected. Finally, the SOA current value data selected in each “time” (time-maximum coherence current corresponding data), that is, current waveform data 317 is output. Next, a detailed procedure for obtaining the SOA current waveform data will be described with a flowchart.
  • FIG. 4 is a flowchart showing a procedure for generating SOA current waveform data to be stored in a memory associated with the wavelength swept light source of the present invention.
  • the reference numerals assigned to the elements (boxes) in the flowchart of FIG. 4 correspond to the numbers of the following steps.
  • step 401 various parameters are set.
  • the initial value, final value, and current change step of the SOA current value data set in the current source 111 are set in the current-time-coherence length correspondence data acquisition unit 313 from the control unit 315.
  • What is set for the current source 111 is not the current but the setting data corresponding to the SOA current value to be set, but in the following description, it will be referred to as the SOA current value or the like for simplicity.
  • the SOA current initial value is set in step 402.
  • the current-time-coherence length correspondence data acquisition unit 313 sets the initial value of the SOA current value in the current source 111. Starting from the initial value, the SOA current value is sequentially changed in the current change step, and the coherence length for each time in one wavelength sweep is measured for each SOA current value in step 403 described later.
  • the SOA current value may be gradually increased as a small value first, or vice versa.
  • the procedure 400 obtains time-coherence length correspondence data. Details of the procedure for obtaining the coherence length will be described later as a third embodiment.
  • the time-coherence length correspondence data acquisition unit 312 acquires the wavelength swept light for one wavelength sweep from the laser resonator 101 in synchronization with the synchronization signal 318 from the optical deflector driver 113. After acquiring the wavelength swept light for one wavelength sweep, the interference condition is changed by changing the optical path length difference between two arms (sample arm and reference arm) of the interferometer described later in the third embodiment, and different interference is performed. A plurality of interference signals corresponding to the conditions are acquired.
  • a PSF for each optical path length difference is calculated for each time during wavelength sweeping. Thereafter, for each time, for each time, the difference between two optical path length differences 2z PSFHALF1 and 2z PSFHALF2 having a PSF peak value at which the PSF peak value (intensity) is 1 ⁇ 2 of the maximum value PSF MAX. Find the value.
  • of the difference value is defined as a coherence length Lc.
  • the signs of the optical path length differences 2z PSFHALF1 and 2z PSFHALF2 are reversed. The following describes the positive / negative optical path length difference and the measurement mechanism.
  • the optical path length difference is negative when the optical path length of the sample arm is short, and vice versa.
  • the length difference is positive.
  • the PSF is calculated for each optical path length difference. The closer the optical path length difference is to 0, the stronger the interference, so the PSF peak value increases. The longer the optical path length difference, the weaker the interference, so the PSF peak. The value becomes smaller. In the above-described measurement, both the positive and negative optical path length differences are measured. For this reason, the PSF peak value becomes maximum when the optical path length difference is near 0, and becomes smaller as the distance from the optical path length difference is 0. Since the coherence length is defined as a coherence distance until the interference intensity is halved from the maximum value, as described above, the coherence length is obtained as an optical path length difference at which the PSF peak value is halved from the maximum value.
  • the measurement may be performed only on the positive or negative side of the optical path length difference.
  • the PSF peak value (intensity) where the optical path length difference 2z is close to 0 is set as the PSF maximum value PSF MAX, and the optical path length difference 2z PSFHALF when the PSF peak value is PSF MAX / 2 is doubled.
  • this length 4z PSFHALF may be used as the coherence length at each time.
  • the PSF is calculated from the interference signal, and substantially the Fourier transform result of the envelope of the interference signal after rescaling becomes the shape of the PSF. This is because when the optical path length difference 2z is 0, the frequency corresponding to the carrier wave of the interference signal is 0, and it is difficult to stably output the envelope of the interference signal.
  • data corresponding to each time and coherence length during wavelength sweeping (time-coherence length correspondence data 316) is output.
  • the procedure 400 acquires current-time-coherence length correspondence data.
  • the current-time-coherence length correspondence data acquisition unit 313 associates the time-coherence length correspondence data 316 output in step 402 with the SOA current value currently set in the current source 111. If current-time-coherence length correspondence data already exists, the time-coherence length correspondence data newly associated with the SOA current in this step is added to the current-time-coherence length correspondence data. If current-time-coherence length correspondence data does not exist, the time-coherence length correspondence data newly associated with the SOA current in this step is set as current-time-coherence length correspondence data.
  • the current-time-coherence length correspondence data is a data set (SOA current value, time, coherence length) each including three types of data, “SOA current value”, “time”, and “coherence length”. ).
  • step 405 the procedure 400 determines whether the SOA current is the final value.
  • the current-time-coherence length correspondence data acquisition unit 313 determines whether the SOA current value currently set in the current source 111 is the final value. If yes, go to step 407, if no, go to step 406.
  • Procedure 400 changes the SOA current value in Step 406.
  • the current-time-coherence length correspondence data acquisition unit 313 changes the SOA current value set in the current source 111 from the current SOA current value to a new SOA current value according to the current change step. Then, the process returns to step 403, and the time-coherence length correspondence data 316 is acquired at the new SOA current value. That is, a data set (time, coherence length) 316 is obtained.
  • step 407 the procedure 400 acquires time-maximum coherence length-corresponding current data (current waveform data).
  • the maximum coherence length-corresponding current selection unit 314 selects the SOA current value that maximizes the coherence length for each wavelength sweep time with respect to the current-time-coherence length correspondence data acquired in step 404. Then, for each time, the SOA current value data (time-maximum coherence length current correspondence data 317), that is, SOA current waveform data is output.
  • FIG. 5 is a table showing the entire coherence length data acquired through the steps shown in the flowchart of FIG. 4 in a simplified manner.
  • the set SOA current value is shown in the vertical direction at the left end of the table. That is, the minimum SOA current value is 30 mA, the maximum SOA current value is 100 mA, and the current change step is 10 mA.
  • the wave number k for which the coherence length was obtained is shown in the horizontal direction at the top of the table shown in FIG.
  • the wave number has a one-to-one correspondence with the wavelength sweep time or wavelength.
  • the wave number is a center wave number k (0 ⁇ ⁇ t) to k (7 ⁇ ⁇ t) at equal time intervals. Seeking coherence length.
  • the coherence length is shown in each cell of Table 5 for the combination of SOA current value (mA) and wave number k of Table 5.
  • a coherence length is obtained for each wave number from k (0 ⁇ ⁇ t) to k (7 ⁇ ⁇ t).
  • eight coherence length data from 9.32 to 10.65 from the left to the right with respect to the center wave numbers k (0 ⁇ ⁇ t) to k (7 ⁇ ⁇ t). Is required.
  • These eight coherence length data correspond to the time-coherence length correspondence data 316 obtained in step 404 in the flowchart of FIG.
  • the wave number k has a one-to-one correspondence with the time (time) in the wavelength sweep period.
  • step 405 in the flowchart of FIG. 4 it is determined whether or not the SOA current value is the final value.
  • the SOA current value is changed to a new value.
  • the coherence is sequentially applied to different SOA current values. This corresponds to repeated acquisition of the length.
  • one SOA current value that maximizes the coherence length is selected for the eight coherence length data arranged in the vertical direction corresponding to one wave number in Table 5.
  • the maximum coherence length 11.11 corresponds to a hatched cell and is obtained when the SOA current is 60 mA. Therefore, 60 mA is selected as the SOA current value at which the coherence length is maximum when the wave number is k (0 ⁇ ⁇ t).
  • the SOA current value that maximizes the coherence length is selected for each of the center wave numbers k (0 ⁇ ⁇ t) to k (7 ⁇ ⁇ t).
  • the SOA current value of the hatched cell is selected.
  • the coherence length calculation becomes abnormal because the interference light is not stable under the conditions of the SOA current value and the center wave number, the oscillation itself disappears. Means.
  • the wavelength swept light source of the present invention further includes a memory associated with the optical gain control means (111) and storing the optical gain control data, and the optical gain control data includes means for changing the optical gain.
  • the means for changing the optical gain by the optical gain control signal having a constant value is controlled to a set optical gain value, and a sweep is performed for each of the set optical gain values having a plurality of different values.
  • a set optical gain value that gives a maximum coherence length is selected from the measured coherence lengths at the one time of the plurality of times, and the one time
  • Each said can also be implemented as being generated from the selected set light gain values.
  • the optical deflector 106 in FIG. 3 is used as a “means for changing the wavelength”
  • the optical deflector driver 113 is used as the wavelength control means
  • the deflector drive signal is used as the wavelength control signal
  • the SOA 102 is used as the “means for changing the optical gain”. Even if the SOA current is replaced with the optical gain control signal or the optical gain, the concept of the acquisition procedure of the SOA current waveform data described with reference to FIGS. 3 to 5 can be applied as it is.
  • the current source 111 in FIG. 3 operates so as to obtain a SOA current value and flow a SOA current having a constant value over time.
  • the SOA 102 is used as a means for changing the optical gain, the optical gain is constant over time. It is sufficient to provide means for controlling the optical gain changing means (constant optical gain control means) so that
  • the current-time-coherence length correspondence data acquisition unit 313 in FIG. 3 operates so as to pass an applied current value for setting the SOA current of the SOA 102 to the current source 111, but a constant optical gain control unit is used instead of the current source 111.
  • the optical gain value is input to the constant optical gain control means to change the optical gain of the optical gain changing means, and the time-coherence length correspondence data acquisition means 312 is used for each optical gain value.
  • time-coherence length correspondence data for each optical gain that is, optical gain-time-coherence length correspondence data is output.
  • the current-time-coherence length correspondence data acquisition unit 313 can be referred to as an optical gain-time-coherence length correspondence data acquisition unit.
  • the maximum coherence length-corresponding current selection means 3114 and the optical gain-time-coherence length correspondence data acquisition means 313 obtained from the optical gain-time-coherence length correspondence data acquisition means 313 is obtained for each time during one wavelength sweep.
  • the optical gain that maximizes the coherence length is selected. Then, the selected optical gain data for each time (time-maximum coherence length gain correspondence data), that is, optical gain control data is output.
  • time interval ⁇ t of data acquisition of the current-time-coherence length correspondence data acquired by the maximum coherence length-corresponding current selection means 314 is different from the time interval ⁇ t ′ of the current controlled by the current source 111, these time intervals make sure that the time is consistent.
  • the time interval of the time-maximum coherence length current correspondence data generated in the maximum coherence length correspondence current selection means 314 is ⁇ t, but the time interval can be recreated to ⁇ t ′ by linear or nonlinear interpolation.
  • the wavelength swept light source of the present invention shown in Example 1 uses the SOA current waveform data obtained by the procedure shown in Example 2. Further, the SOA current waveform data obtained in the second embodiment is obtained using the apparatus and procedure for obtaining the coherence length of the third embodiment described below.
  • the control unit 315 including a CPU or the like is originally provided in the wavelength sweep light source shown in FIG.
  • both of the optical gain-time-coherence length correspondence data acquisition means 313 and the maximum coherence length correspondence current selection means 314 can be realized by arithmetic processing by software using a CPU or the like. Therefore, the wavelength swept light source of the present invention can be implemented even with all the components shown in FIG. Since the procedure of the flowchart of FIG. 4 of the second embodiment is performed using the wavelength swept light, the SOA current waveform data is acquired even when the wavelength swept light source is being used for an actual light source. It is also possible to update the SOA current waveform data.
  • the wavelength sweeping light source is always updated by updating the SOA current waveform data when the wavelength sweeping light source is in an operating environment in accordance with a change in environmental conditions such as temperature when using the wavelength sweeping light source. It can also be operated in an optimal state.
  • FIG. 6 is a diagram showing a configuration of an apparatus for measuring a coherence length for obtaining SOA current waveform data used in the wavelength swept light source of the present invention.
  • a detailed configuration of the time-coherence length correspondence data acquisition unit 312 shown in FIG. 3 is shown in a dotted line region.
  • the laser resonator 101, current source 111, and optical deflector driver 113 in FIG. 6 are the same as those in FIG. 3 of the second embodiment.
  • the outline of the procedure for obtaining the coherence length is as follows.
  • the time-coherence length correspondence data acquisition unit 312 acquires the wavelength swept light 123 for one wavelength sweep from the laser resonator 101 in synchronization with the synchronization signal 318 from the optical deflector driver 113. After acquiring the wavelength swept light for one wavelength sweep, the interference condition is changed by changing the optical path length difference between the two arms (sample arm and reference arm) of the interferometer 501, and a plurality of corresponding interference conditions are changed. Acquire the interference signal. After rescaling the plurality of acquired interference lights 521, the time-optical path length difference-PSF correspondence data acquisition means 515 calculates the PSF for each optical path length difference at each time during wavelength sweeping.
  • the optical path length difference is negative when the optical path length of the sample arm is short, and vice versa.
  • the length difference is positive.
  • the PSF is calculated for each optical path length difference. The closer the optical path length difference is to 0, the stronger the interference, so the PSF peak value increases. The longer the optical path length difference, the weaker the interference, so the PSF peak. The value becomes smaller. In the above-described measurement, both the positive and negative optical path length differences are measured. For this reason, the PSF peak value becomes maximum when the optical path length difference is near 0, and decreases as the distance from the PSF peak value increases. Since the coherence length is defined as a coherent distance until the interference intensity is halved from the maximum value, as described above, the coherence length is obtained as an optical path length difference at which the PSF peak value is halved from the maximum value.
  • the optical path length difference 2z 0 is the center and the PSF peak value attenuation method is bilaterally symmetric with respect to either positive or negative, even if the optical path length difference is measured only on the positive side or the negative side, good.
  • the PSF peak value (intensity) where the optical path length difference 2z is close to 0 is set as the PSF maximum value PSF MAX, and the optical path length difference 2z PSFHALF when the PSF peak value is PSF MAX / 2 is doubled.
  • the length 4z PSFHALF may be used as the coherence length Lc at each time.
  • the reason why the PSF peak value where the optical path length difference 2z is close to 0 is set as the PSF maximum value PSF MAX because the PSF when the optical path length difference 2z is 0 cannot be calculated.
  • the PSF is calculated from the interference signal, and substantially the Fourier transform result of the envelope of the interference signal after rescaling becomes the shape of the PSF. This is because when the optical path length difference 2z is 0, the frequency corresponding to the carrier wave of the interference signal is 0, and it is difficult to stably output the envelope of the interference signal.
  • a plurality of PSF peak points (2zq, PSF peak_q ) having different optical path length differences (where q is a plurality of PSF peak points) instead of the PSF at the position where the optical path length difference 2z is 0 instead of the maximum PSF peak_MAX_PSF. It is also possible to perform function f (2z) fitting using (integer) and set the function value f (0) at which the optical path length difference 2z 0 to be the PSF maximum value PSF peak_MAX .
  • the function f (x) is preferably, for example, a Gaussian, a Laurentian, a function polynomial having x as an argument, or the like.
  • correspondence data (time-coherence length correspondence data 316) of each time during the wavelength sweep and the coherence length Lc is output.
  • the elements in the time-coherence length correspondence data acquisition unit 312 in FIG. 6 represent step elements of a procedure for obtaining the coherence length Lc as functional blocks. That is, elements other than the interferometer 501 represent functional blocks indicating processing processes that can be realized by software arithmetic processing by a CPU or the like. Therefore, it should be noted that, except for the interferometer 501, it can be realized as software arithmetic processing without any special hardware element other than the processor.
  • FIGS. 9A to 9D A specific configuration example of the interferometer 501 will be described last with reference to FIGS. 9A to 9D. Hereinafter, functions and operations of the respective processing blocks in FIG. 6 will be described in detail.
  • the wavelength swept light 123 from the laser resonator 101 is supplied to the interferometer 501.
  • the interferometer 501 emits the wavelength swept light 123 from the laser resonator 101 to the space by a fiber collimator (FC) 502.
  • FC fiber collimator
  • the emitted light is divided into two paths by a beam splitter (BS) 503 on the reference arm side and the sample arm side.
  • a mirror 504 and a mirror 505 are installed on the reference arm and the sample arm, respectively.
  • Light incident on each arm is reflected and returned to the BS 503, and is combined and interfered by the BS 503.
  • the interference light enters a photodetector (PD) 507 and is output from the PD 507 to the outside of the interferometer 501 as an electric signal (interference signal) 521.
  • PD photodetector
  • the mirror 505 on the sample arm side is installed on the mirror position moving means 506, and the distance between the BS 503 and the sample arm side mirror 505 can be controlled from the outside. That is, the difference z between the reference arm length and the sample arm length can be controlled from the outside of the interferometer 501.
  • the interference signal acquisition unit 511 receives a mirror movement completion signal 532 from a time-optical path length difference-PSF correspondence data acquisition unit 515 described later, and then synchronizes with the synchronization signal obtained from the optical deflector driver 113.
  • the interference signal 123 for the wavelength sweep is acquired and output.
  • the mirror movement completion signal 532 is a signal indicating that the mirror position movable means 506 in the interferometer 501 has finished one movement of the sample arm side mirror for changing the optical path length difference.
  • the rescaling means 512 rescals the interference signal 522 obtained from the interference signal acquisition means using the time-wave number correspondence data, and the rescaled interference signal 523, i (k) (k is the wave number). Is output.
  • the rescaling is described in detail in Non-Patent Document 1, but operates so as to resample the interference signal 522 (resampling) and expand and contract the time axis. As a result, the interference signal 522 is converted into a signal having a constant frequency with respect to time.
  • the extraction position calculation means 517 obtains the wave number corresponding to the time 531 obtained from the processing time generation means 518, that is, p ⁇ t (p is 0 to n ⁇ 1), that is, the extraction position 536 from the time-wave number correspondence data, and outputs this. To do. Specifically, the extraction position (wave number) k (p ⁇ t) is output.
  • the processing range extraction means 513 uses the extraction width ⁇ k and the extraction position k (p ⁇ t) obtained from the extraction position calculation means 518 to convert the resected interference signal 523 obtained from the rescaling means 512, i.e., i (k), to k. Extract and output within a range of ⁇ ⁇ k / 2 with (p ⁇ t) as the center.
  • extraction may be performed using a window function (for example, a square window, a Gaussian window, a Hanning window, a Hamming window, a Blackman window, or the like) having a center k (p ⁇ t) and a half width ⁇ ⁇ k / 2.
  • the window function / frequency conversion means 514 multiplies the interference signal 524 obtained from the processing range extraction means 513 with a window function, and uses a Fourier transform to obtain a power spectrum (or frequency distribution excluding phase, that is, the square root of the power spectrum). And only the positive side is output.
  • the output signal 525 is called a point spread function (PSF).
  • the window function need not be used when the processing range extracting means 513 uses a window whose intensity decreases as it goes away from the center, such as a Gaussian window, Hanning window, Hamming window, or Blackman window. Further, the presence / absence of the window and the specification / type may be determined according to the discontinuous point of the interference signal.
  • the time-optical path length difference-PSF correspondence data acquisition means 515 obtains the PSF data from the window function / frequency conversion means 514 and outputs the time p ⁇ t (p is 0 to n ⁇ 1), which is output from the processing time generation means 518.
  • the optical path length difference between the reference arm and the sample arm from FC502 to PD507 is 2z.
  • the mirror position signal 533 is sent to the mirror position moving means 506 incorporated in the interferometer 501 to change the position of the sample arm side mirror 505.
  • the mirror movement completion signal 532 is obtained as an interference signal.
  • m ⁇ z / ⁇ z
  • q integer of ⁇ m to m.
  • the time-optical path length difference-PSF correspondence data acquisition means 515 detects whether the movement of the position of the sample arm side mirror 505 has been completed.
  • the movement completion signal 532 is obtained from the mirror position movable means 506 and the movement is completed. You can judge. Further, the position of the mirror on the sample arm side of the mirror position movable means 506 is always observed, and if it matches the mirror position signal or falls within a preset error, it is determined that the movement is completed. Also good. Further, it may be determined that the movement of the mirror is completed after a sufficient time has passed for the movement of the mirror.
  • the time-optical path length difference-PSF correspondence data acquisition means 515 outputs the time-optical path length difference-PSF correspondence data 534 when p at time p ⁇ t obtained from the processing time generation means 518 reaches n ⁇ 1.
  • the output time-optical path length difference-PSF correspondence data 534 stores PSF data in which the time is in the range of 0 to (n ⁇ 1) ⁇ t and the optical path length difference is in the range of ⁇ m ⁇ z to m ⁇ z.
  • the time-coherence length correspondence data generation means 516 uses the time-optical path length difference-PSF correspondence data acquisition means 515 to obtain each time p ⁇ t (p is 0 to n-1) from the time-optical path length difference-PSF correspondence data 534.
  • the coherence length Lc is calculated, and p ⁇ t and Lc are associated and stored.
  • the time-coherence length correspondence data 316 is output.
  • the peak value of the PSF corresponding to each optical path length difference is obtained every time p ⁇ t, and first, the PSF peak value that is 1/2 of the maximum PSF MAX value PSF MAX is obtained.
  • 2 ⁇ z
  • of the difference value is set as a coherence length Lc.
  • the optical path length difference is negative when the optical path length of the sample arm is short, and vice versa.
  • the length difference is positive.
  • the PSF is calculated for each optical path length difference. The closer the optical path length difference is to 0, the stronger the interference, so the PSF peak value increases. The longer the optical path length difference, the weaker the interference, so the PSF peak. The value becomes smaller. In the above-described measurement, both the positive and negative optical path length differences are measured. For this reason, the PSF peak value becomes maximum when the optical path length difference is near 0, and becomes smaller as the distance from the optical path length difference is 0. Since the coherence length is defined as a coherent distance until the interference intensity is halved from the maximum value, as described above, the coherence length is obtained as an optical path length difference at which the PSF peak value is halved from the maximum value.
  • the peak value (intensity) of the PSF whose optical path length difference 2z is close to 0 is set as the PSF maximum value PSF MAX, and the optical path length difference 2z PSFHALF when the PSF peak value is PSF MAX / 2 is doubled.
  • this length 4z PSFHALF may be used as the coherence length Lc.
  • a plurality of PSF peak positions and peak values having different optical path length differences can be Gaussian fitted, and the full width at half maximum of Gaussian can be set as the coherence length Lc.
  • Gaussian Laurentian, polynomial of function, etc. are suitable.
  • the PSF peak value (intensity) where the optical path length difference 2z is close to 0 is set as the PSF maximum value PSF MAX, and the optical path length difference 2z PSFHALF when the PSF peak value is PSF MAX / 2 is doubled.
  • the length 4z PSFHALF may be used as the coherence length at each time. The reason why the PSF peak value where the optical path length difference 2z is close to 0 is used as the PSF maximum value PSF peak_MAX is that a PSF having an optical path length difference 2z of 0 cannot be calculated.
  • a plurality of PSF peak points (2zq, PSF peak_q ) having different optical path length differences (where q is a plurality of PSF peak points) instead of the PSF at the position where the optical path length difference 2z is 0 instead of the maximum PSF peak_MAX_PSF. It is also possible to perform function f (2z) fitting using (integer) and set the function value f (0) at which the optical path length difference 2z 0 to be the PSF maximum value PSF peak_MAX . Examples of the function f (x) are Gaussian, Laurentian, a polynomial of a function having x as an argument, and the like.
  • FIG. 7 is a flowchart showing a procedure for measuring a coherence length for obtaining SOA current waveform data used in the wavelength swept light source of the present invention.
  • the reference numerals assigned to the elements (boxes) in the flowchart of FIG. 7 correspond to the numbers of the following steps.
  • the procedure 600 for obtaining the coherence length Lc starts by setting various parameters in step 601.
  • step 601 the following parameters are set.
  • the mirror position range ⁇ ⁇ z and the mirror position step ⁇ z are set in the time-optical path length difference-PSF correspondence data acquisition means 515.
  • the processing time generating means 518 one sweep time T and time step ⁇ t are set.
  • the time-wave number correspondence data k (t) is set in the rescaling means 512 and the extraction position calculating means 517.
  • step 602 the procedure 600 sets an initial optical path length difference.
  • the time-optical path length difference-PSF correspondence data acquisition unit 515 outputs an initial mirror position signal 533 to the mirror moving unit 506.
  • the mirror position moving means 506 moves the sample arm side mirror 505 to a position (initial position) according to the mirror position signal 533.
  • the mirror position signal 533 is any signal representative of the mirror position z.
  • m ⁇ z / ⁇ z
  • q is an integer from ⁇ m to m. In the ascending order, the sample arm is moved to the position of ⁇ m ⁇ z, and in the descending order, it is moved to the position of m ⁇ z.
  • the time-optical path length difference-PSF correspondence data acquisition unit 515 outputs a mirror movement completion signal 532 to the interference signal acquisition unit 511 after moving the mirror.
  • step 603 the procedure 600 executes acquisition of an interference signal.
  • the interference signal acquisition unit 511 receives the mirror movement completion signal 532 from the time-optical path length difference-PSF correspondence data acquisition unit 515, the interference signal acquisition unit 511 receives 1 from the interferometer 501 in synchronization with the synchronization signal 318 obtained from the optical deflector driver 113. An interference signal 521 corresponding to the number of wavelength sweeps is acquired and output.
  • Procedure 600 performs rescaling in step 604.
  • the rescaling unit 512 rescals the interference signal 521 for one wavelength sweep obtained from the interference signal acquisition unit 511 using the time-wavenumber correspondence data.
  • Procedure 600 generates an initial processing time in step 605.
  • Procedure 600 calculates the extraction position in step 606. Based on the time-wave number correspondence data, the extraction position calculation means 517 obtains the wave number 536 corresponding to the time p ⁇ t (p is 0 to n ⁇ 1) obtained from the processing time generation means 518, that is, k (p ⁇ t). Is output.
  • the extraction position means the time for extracting the interference signal in the subsequent steps. Recall that in the wavelength sweep period of the wavelength swept light source, the time (time) from the start of the wavelength sweep has a one-to-one correspondence with the wavelength or wave number of the oscillation output light.
  • the procedure 600 extracts a processing range.
  • the processing range extraction unit 513 uses the extraction width ⁇ k and the extraction position 536 obtained from the extraction position calculation unit 517, that is, k (p ⁇ t), to use the resected interference signal 523 obtained from the rescaling unit 512, that is, i (k). Are extracted in a range of ⁇ ⁇ k / 2 around k (p ⁇ t), and the extracted interference signal 524 is output.
  • a window function having a center k (p ⁇ t) and a half width ⁇ ⁇ k / 2, for example, a rectangular window, a Gaussian window, a Hanning window, a Hamming window, a Blackman window, or the like may be used.
  • step 608 the procedure 600 performs frequency conversion by applying a window function to the extracted interference light.
  • the window function / frequency conversion means 514 multiplies the interference signal 524 obtained from the processing range extraction means 513 with a window function, and uses Fourier transform to obtain the power spectrum (or frequency distribution excluding the phase, that is, the square root of the power spectrum). ) And output only the positive frequency.
  • the output signal 525 is a point spread function (PSF).
  • PSF point spread function
  • the use of the window function in this step is not essential and may not be used. In particular, when a window other than the rectangular window is used in step 607, it is not necessary to use the window function in this step.
  • step 609 the procedure 600 obtains and accumulates time-optical path length difference-PSF correspondence data.
  • a memory associated with the wavelength swept light source can be used.
  • Procedure 600 performs time update in step 610.
  • the processing time generation means 518 updates the variable p at the time p ⁇ t. Specifically, p + 1 is updated to p. That is, (p + 1) ⁇ t is set as the updated time.
  • step 611 the procedure 600 determines the end time.
  • step 612 the procedure 600 updates the optical path length difference.
  • a value obtained by adding 1 to q is updated to become a new q.
  • a value obtained by subtracting 1 from q is updated to become a new q.
  • step 613 the procedure 600 determines whether the optical path length difference update is completed.
  • step 614 the procedure 600 generates time-coherence length correspondence data.
  • the time-coherence length correspondence data generation means 516 calculates the time p ⁇ t (p is 0 to n-1) from the time-optical path length difference-PSF correspondence data acquisition means 515 from the time-optical path length difference-PSF correspondence data.
  • a coherence length Lc is calculated, and p ⁇ t and Lc are associated and stored.
  • the time-coherence length correspondence data 316 is output.
  • step 614 is completed, the coherence length for each wave number is obtained for one SOA current. That is, referring to FIG. 5 again, for one SOA current value, the coherence length is obtained for each wave number from k (0 ⁇ ⁇ t) to k (7 ⁇ ⁇ t).
  • the optical gain control data is output from the laser resonator set to the set optical gain value, and is a wavelength swept light obtained from one wavelength sweep in the wavelength sweep period ( 123) based on a plurality of interference signals (521) corresponding to each of the interference conditions generated from an interferometer that causes interference under different interference conditions, the coherence length (316) determined from the plurality of interference signals ).
  • FIG. 8 is a diagram for explaining the operation of the SOA current obtained at equal time intervals in the wavelength swept light source of the present invention.
  • FIG. 8 shows the SOA current waveform obtained according to the procedure for obtaining the SOA current waveform data of Example 2 based on the coherence length Lc obtained in Example 3 described above.
  • the coherence length is calculated at equal wave number intervals, and the SOA current that maximizes the coherence length is determined.
  • the determined SOA current intervals are non-uniform when viewed on the time axis.
  • the coherence length may be calculated even at equal wavelength intervals, but in this case as well, the determined time interval of the SOA current is not uniform.
  • the SOA current waveform of FIG. 8 determines the calculation time so that the time interval ( ⁇ t) for calculating the coherence length is uniform, as described in the third embodiment, at equal time intervals. is doing. Therefore, it is possible to perform the Fourier transform on the zero-padded data and to interpolate the SOA current more simply and at high speed.
  • FIG. 2 and FIG. 8 only the number of wavenumbers k for calculating the coherent length is drawn, but actually, the coherence length is calculated for a maximum of 5000 to 10000 points depending on the wavelength sweep speed. It is preferable to do.
  • the SOA current is taken in a current variable step of 10 mA, but it goes without saying that it can be taken in finer steps.
  • FIGS. 9A to 9D are diagrams showing specific configuration examples of the interferometer in the time-coherence length correspondence data acquisition unit shown in FIG.
  • FIG. 6 shows a configuration example of the Michelson interferometer 501, it can be realized by various configurations of Mach-Zehnder type as shown in FIGS. 9A to 9D.
  • 9A shows a Mach-Zehnder type (differential output type, fiber interferometer type)
  • FIG. 9B shows a Mach-Zehnder type (differential output type)
  • FIG. 9C shows a Mach-Zehnder type (differential output type, fiber interferometer type)
  • FIG. 9D shows a Mach-Zehnder type.
  • FC fiber collimator
  • CPL coupler
  • CCL circulator
  • BPD balanced detector
  • Ms sample arm mirror
  • M R reference arm mirror
  • the coherence length Lc in the wavelength swept light source of the present invention in a laser resonator that uses a galvano scanner and the SOA current can be varied in the range of 0 to 600 mA, in the wavelength sweep range of the center wavelength of 1310 nm ⁇ 40 nm.
  • the coherence length Lc could be increased 1.7 times over the entire sweep wavelength band.
  • the coherence length of the wavelength swept output light during the wavelength sweep period can always be maintained long by the wavelength swept light source and the method of creating the drive data of the present invention.
  • the PSF intensity in the SS-OCT apparatus can be increased and the image detection sensitivity of the SS-OCT apparatus can be increased.
  • the optical gain in the laser resonator is slightly higher than the optical gain in the laser oscillation threshold state over the entire wavelength band of the wavelength sweep.
  • the optical gain changing means SOA
  • the coherence length may be maximized at an optical gain value that is significantly higher than that near the laser oscillation threshold. Even in such a case, by obtaining the SOA current waveform that gives the maximum coherence length as in the wavelength swept light source of the present embodiment, the state in which the coherence length of the output light is long over the entire sweep wavelength band. Can be maintained.
  • SOA current waveform data of the SOA used for the wavelength swept light source is obtained in advance, and the wavelength swept light source is operated based on the acquired data. Therefore, before the wavelength swept light source is normally operated, a process of acquiring this SOA current waveform data is necessary.
  • the procedure for obtaining the coherence length for acquiring the SOA current waveform data described in the third embodiment requires the interferometer 501 in the processing block diagrams shown in FIGS. 6 and 9A to 9D. As a device of this type, originally unnecessary components are required.
  • the wavelength sweep output during the wavelength sweep period is obtained by a simpler component and processing based on the output level of the light source without obtaining the coherence length at each sweep wavelength, as in the first embodiment.
  • An example is shown in which the coherence length of light is always kept long.
  • the SOA is characterized by a wide bandwidth and a short carrier lifetime, and is excellent in practicality that can be driven by an electric signal directly by current injection without requiring optical excitation.
  • FIG. 10 is a diagram showing the relationship between the injection current into the SOA, the optical output power, and the coherence length in the wavelength swept light source.
  • the horizontal axis represents the injection current into the SOA, and the vertical axis represents the coherence length and output light level at a specific wavelength.
  • the output light level increases as the injection current into the SOA increases, whereas the coherence length decreases conversely.
  • the coherence length rapidly increases in a state where the SOA injection current is close to the oscillation threshold value of 25 mA or less at a specific wavelength.
  • the wavelength swept light source of the present embodiment by monitoring the output light level in real time during the normal operation of the wavelength swept light source, it is possible to maintain a long coherence length over the entire swept wavelength band by a simpler control method.
  • the injection current to the SOA is controlled based on the level of the output light in conjunction with the operation of the means (wavelength selection element) for changing the oscillation wavelength in the wavelength swept light source. This drives the SOA near the oscillation threshold at each wavelength within the sweep wavelength band, thereby narrowing the instantaneous oscillation line width of the light source and maximizing the coherence length.
  • the wavelength swept light source according to the present embodiment immediately approaches a state where the coherence length is extended at least to a certain extent if a control step linked with wavelength sweep described later is repeated several times or more, and during the wavelength sweep operation that is a normal operation of the light source. There is no need to measure the coherence length.
  • the determination of the target value of the optical output power, which will be described later, can be performed only in the manufacturing adjustment stage of the wavelength swept light source device, and the interferometer required for obtaining the coherence length in the wavelength swept light source of the first embodiment is also available. It is unnecessary.
  • FIG. 11 is a diagram showing a configuration of a wavelength swept light source according to the second embodiment of the present invention.
  • the wavelength swept light source 700 includes, as components, an optical amplifying element (SOA) 701 that amplifies light, a means (wavelength selecting element) 702 that changes an oscillation wavelength for selectively filtering light in a specific wavelength band, An optical feedback element 705 such as a half mirror for branching the part back to the light source, a driver apparatus 709 for driving the means 702 for changing the oscillation wavelength, and a driver apparatus 708 for driving the optical amplification element (SOA) 701 are included.
  • the plurality of optical elements included in the means (wavelength selection element) 702 for changing the oscillation wavelength forms an optical reciprocating optical path as will be described later, and an optical resonator capable of repeatedly reciprocating light between these elements. It is configured as.
  • the wavelength selection element 702 includes an optical deflector 703 and a diffraction grating 704.
  • the diffraction grating has a Littman type arrangement, but in the wavelength swept light source of the present embodiment, the diffraction grating 704 is arranged under a Littrow condition in which the diffracted light is diffracted in the opposite direction to the incident light. ing.
  • an element based on another wavelength selection principle such as a Fabry-Perot filter may be used.
  • the optical deflector 703 As the optical deflector 703, a galvanometer mirror, a MEMS mirror, a KTN optical deflector, or the like can be used. A part of the other output light of the SOA 701 is reflected by the optical feedback element 105 and fed back to the resonator formed between the diffraction grating 704 and the optical feedback element 705. The light transmitted through the optical feedback element 705 is collected as it is as output light by the lens 706 and coupled to the optical fiber 707 to become output light 712 of the light source.
  • the resonator in the wavelength swept light source of FIG. 11 has a linear configuration, this is an example, and each resonator may be configured by a ring resonator in which optical elements are connected by optical fibers.
  • the operation of narrowing the line width of the output light in the resonator of the wavelength swept light source is explained as follows.
  • the light output from the SOA 701 has a wide optical spectrum.
  • the spectrum filtered by the wavelength selection element 702 and cut out by the filter shape of the wavelength selection element 702 is incident on the SOA 701 again.
  • the signal is amplified by the SOA 710 and output from the end face on the opposite side of the SOA 701.
  • a part of the light is returned again to the SOA 701 by the optical feedback element 705, and thus circulates in the optical resonator a plurality of times during a certain period of time.
  • the selected wavelength of the wavelength selection element 702 changes, and this certain time is determined by the time that the filter spectrum crosses a specific wavelength. As the light passes through the wavelength selection element 102 a plurality of times within the crossing time, the line width of the output light is narrowed by the cumulative effect of filtering.
  • the wavelength dependence of the laser oscillation threshold occurs due to the influence of the gain spectrum of the SOA 701 used and the wavelength dependence of the resonator loss. Therefore, the optimum SOA drive current value for maintaining a long coherence length differs for each wavelength.
  • the optical output power from the light source is monitored and controlled in order to set the SOA to an optimum SOA drive current value that differs for each wavelength.
  • the optical output power from the wavelength swept light source depends on the number of times the SOA 701 receives amplification.
  • the coherence length also depends on the number of passes through the wavelength selection element 702. For this reason, the optical output power and the coherence length are respectively associated with the number of oscillations of the oscillation light in the resonator. If the transmission characteristic of the wavelength selection element 702 has a Gaussian spectrum, the shape function F ( ⁇ ) is expressed by the following equation.
  • represents the maximum transmittance of the wavelength selection element 702
  • ⁇ C represents the center wavelength
  • ⁇ F represents the full width at half maximum of the spectrum.
  • the state in which the net gain G net (net gain) generated in a portion other than the wavelength selection element 702 when light makes one round of the optical resonator balances the loss generated in the wavelength selection element 702 is the state of the laser oscillation threshold. It becomes. Since the loss of the optical resonator is determined by the characteristics of each component, the net gain G net mainly varies with the gain of the SOA 701.
  • the loss L in the wavelength selection element 702 is expressed by the following expression in dB.
  • Equation (4) represents the half width of the wavelength band in which the net gain G net is larger than the loss L and exceeds the oscillation threshold. Therefore, the full width ⁇ eff of the wavelength region exceeding the oscillation threshold is obtained by doubling the half width of the equation (4) and combining the equation (5).
  • ⁇ eff in Expression (6) represents an effective wavelength selection width related to laser oscillation, and is also referred to as an effective wavelength selection width. It should be noted here that when the spectrum width ⁇ F of the wavelength selection element 702 can be regarded as invariant regardless of the wavelength (condition 1), the effective wavelength selection width ⁇ eff related to laser oscillation is the net gain G net, It is determined by the drive current of the SOA 701. This means that the effective wavelength selection width ⁇ eff can be kept constant by making the net gain G net constant at each wavelength. Further, assuming that the selection wavelength of the wavelength selection element 702 changes at a constant speed within the sweep wavelength band (condition 2), if the effective wavelength selection width ⁇ eff is constant, the number of times of light circulation in the resonator at each wavelength Is also constant.
  • the control voltage of the optical deflector 703 is preferably a waveform in which the selected wavelength changes at a constant speed within the sweep wavelength band.
  • the wavelength selection element 102 has a sine wave as a control signal (wavelength sweep signal).
  • the sweep speed varies depending on the wavelength, and even if the effective wavelength selection width ⁇ eff is constant, the number of times of light circulation in the resonator varies depending on the wavelength. Therefore, even if the optical output power is made constant, the coherence length is not always constant. However, the change speed of the selected wavelength of the wavelength selection element 702 near the center of the sweep wavelength band is substantially constant, and even when driven by a sine wave, the coherence length is almost constant around the center of the sweep wavelength band. Is obtained. Although the above discussion is expressed in the wavelength space, the effect of the present invention can be obtained in the same manner even if these are replaced with the frequency space.
  • the effective wavelength selection width ⁇ eff is further narrowed, and oscillation occurs at each wavelength as shown in FIG.
  • the coherence length can be maximized in the vicinity of the threshold state.
  • the coherence length seems to increase exponentially as it approaches the oscillation threshold state (dotted line), but actually the oscillation state becomes unstable as it approaches the oscillation threshold. For this reason, the coherence length is deteriorated by mode competition between the laser longitudinal modes. Considering this instability, the optimum values of the net gain G net and the SOA drive current value are slightly higher than the oscillation threshold (on the right side of the dotted line in FIG. 10).
  • the SOA drive current value is controlled to increase or decrease in time (synchronized with the corresponding wavelength) corresponding to the wavelength sweep which is the actual operation of the light source.
  • the wavelength swept light source of the present embodiment can be implemented as a wavelength swept light source device including this control circuit.
  • FIG. 12 shows a configuration of a wavelength sweep light source device including a wavelength sweep light source and a control circuit that controls the injection current of the SOA to keep the coherence length of the wavelength sweep output light always long.
  • the wavelength swept light source device 720 includes the wavelength swept light source 700 shown in FIG. 11, and further includes an optical branching element 722 that branches out a part of the output light from the wavelength swept light source 700, and the optical output power of one branched light 728.
  • a photo detector 723 for measuring the output signal
  • an A / D converter 724 for converting the output signal of the photo detector 723 into a digital signal
  • processing the digital signal to generate control waveform data to be a control signal to the SOA driver device 708
  • an arbitrary waveform generator 726 that outputs a waveform signal 711 generated by the calculator.
  • a part of the output light from the wavelength swept light source 700 is taken out by the optical branching element 722, one branch light is detected as the light source output light 727, and the other branch light 728 is detected by the photodetector 723 as the monitor signal light.
  • the monitoring signal light is converted into a voltage signal 729 and sampled by the A / D converter 724.
  • the sampled signal is input to the arithmetic unit 725.
  • the arithmetic unit 725 compares the input signal with a set value 730 described later, generates an error signal, and generates a control waveform 731 corresponding to the error signal.
  • Set value 730, and the threshold current for each wavelength below, the initial current, tolerance D, include such target optical output power value V t, may be read from memory or a storage device, not shown.
  • the control waveform 731 is input to an arbitrary waveform generator 726, and the arbitrary waveform generator 726 outputs the data waveform of the input waveform as a control signal 711 in the form of a voltage waveform in which data points are smoothly interpolated.
  • the current driver 708 in the wavelength swept light source 700 has a function of converting the voltage waveform input from the arbitrary waveform generator 726 into an SOA drive current waveform, and injects the SOA drive current corresponding to the input voltage waveform into the SOA 701. To do.
  • a driver device 709 that drives a means (for example, an optical deflector) 702 that changes the oscillation wavelength supplies a trigger signal 710 synchronized with the wavelength sweep to the arbitrary waveform generator 726.
  • the calculator 725 is described as including a D / A conversion function and outputting the control waveform 731.
  • the waveform data is output from the calculator 725, and the D / A converter provided separately is used.
  • a control waveform 731 may be obtained.
  • the configuration from the A / D converter 724 to the arbitrary waveform generator 726 is not limited to the configuration shown in FIG. 12, and various modifications are possible.
  • a method of generating a time-SOA drive current waveform by the bisection method, which is performed by the calculator 725 will be described in detail.
  • the procedure of the method for generating a time-SOA drive current waveform represents a processing process that can be realized by software arithmetic processing by the calculator 725 or the CPU.
  • FIG. 13 is a diagram for explaining the relationship between the current injected into the SOA and the optical output level in the wavelength swept light source of the present embodiment.
  • the horizontal axis represents the injection current into the SOA, that is, the SOA drive current
  • the vertical axis represents the photodetector output voltage obtained from the photodetector 723 of FIG.
  • the photodetector output voltage corresponds to the light output level from the wavelength swept light source 700.
  • N optical output power data V 0 , V 1 ,..., V N ⁇ 1 that can be acquired by sampling with the A / D converter 724 in one wavelength sweep operation are considered.
  • Each of the N data points corresponds to optical output power data of wavelengths ⁇ 0 , ⁇ 1 ,..., ⁇ N ⁇ 1 included in the entire sweep wavelength band.
  • Each of the optical output power data corresponds one-to-one with the voltage signal 729 of the monitoring signal light in FIG. 12, and corresponds to the voltage of the photodetector output on the vertical axis in FIG.
  • a solid line 805 represents the SOA injected current vs. optical output power characteristics at the wavelength ⁇ n (I-P characteristic).
  • a dotted line 806 shows the SOA injection current versus optical output power characteristic (IP characteristic) at the wavelength ⁇ n ⁇ 1 as a different example of the optical output power data.
  • the optical output power values 804 as a target and V t, the injection current value when V n 803 is obtained to I n.
  • the difference between the current optical output power V n and the target light output power value V t, is defined as an error signal E n as follows.
  • the target optical output power value V t as a numerical value used for controlling the optical output power, it is possible to arbitrarily set a target optical output power value V t.
  • set to a certain value in advance V t using the control flow described below light output performs control for constant the power, further changing the value of the V t Repeat this control, searches for an optimum value. At this time, the search can be performed by measuring the coherence length by the method shown in the first embodiment.
  • the optical output power with respect to the SOA drive current value at a certain time also varies depending on the SOA drive current value injected at the moment before and after the wavelength sweep. Therefore, feedback control is necessary to generate the SOA drive current waveform corresponding to the wavelength sweep. Since the optical output power cannot be obtained when the injection current to the SOA is equal to or less than the oscillation threshold current, the lower limit 801 of the injection current value is the oscillation threshold current Ith , n .
  • the current value I ini, n to be injected initially is set as a value obtained by multiplying the oscillation threshold current by the coefficient a as follows.
  • I ini, n aI th, n formula (8)
  • the value of the coefficient a in the above formula is set as an error signal E n is a negative value when injecting a initial current at each wavelength. That is, at all data points of the N pieces of optical output power data V 0 , V 1 ,..., V N ⁇ 1 , the initial optical output power data is larger than the target optical output power value V t when the initial current is injected. , V n > V t is satisfied.
  • the control step [Delta] I n the time of feedback control of the SOA driving current is defined by the following equation.
  • an allowable error between the optical output power value V n monitored by the photodetector 723 and the target optical output power value V t is defined as D.
  • SOA current value I n at a wavelength lambda n is the initial current value I ini, n defined by Equation (8).
  • M represents that the feedback control is the Mth time.
  • ': ' Indicates that the numerical value on the right side is assigned to the variable on the left side.
  • E n ⁇ 0 (including 0) is shown.
  • E n ⁇ 0 (including 0)
  • the sign of ⁇ I n / 2 M in Expression (10) is positive.
  • the optical gain control means (723 to 726) detects the output light level from the laser resonator, and the output light level is determined based on the detected output light level.
  • the optical gain control signal is generated so as to obtain a predetermined output light level value, and the means 710 for changing the optical gain is controlled by the optical gain control signal.
  • the process for controlling the SOA drive current from the optical output power V n described above is described as a series of controls at a certain wavelength ⁇ n, but in the wavelength swept light source device 720 of the present embodiment, during one wavelength sweep period.
  • the above-mentioned N pieces of optical output power data V 0 , V 1 ,..., V N ⁇ 1 can be acquired.
  • ⁇ until the condition of D is satisfied, will perform the update of the SOA driving current value I n at a wavelength that does not satisfy the condition. That is, each iteration of wavelength sweeping operation, the optimization process of the SOA driving current shown in FIG. 13, the control updates the SOA current value I n to obtain light output power value V t as a target at each wavelength is repeated It is.
  • FIG. 14 is a diagram conceptually illustrating the SOA drive current optimization process including the wavelength sweep time in the wavelength sweep light source device 720 of the present embodiment.
  • a series of progress of the SOA drive current optimization process is shown as time-SOA drive current time waveforms 812, 813, and 814.
  • 14 shows a time-wavelength waveform 820 of the output light from the wavelength swept light source in which the wavelength sweep times are aligned and linked with the time-SOA drive current time waveform.
  • the one-dot chain line 811 at the bottom indicates the time change of the oscillation threshold current of the SOA. Since the time on the horizontal axis in the lower graph of FIG.
  • the oscillation threshold current also corresponds to the wavelength on a one-to-one basis.
  • the top one-dot chain line 812 indicates the time-initial current value I ini, n defined by the equation (8).
  • a plurality of sample points corresponding to the initial current value I ini, n initially set at the start of the optimization process are shown at certain wavelength intervals.
  • the wavelength intervals at which these sample points are acquired may be equal intervals or non-equal intervals.
  • the optimum SOA drive current waveform 815 of the target value corresponding to the optical output level that maximizes the coherence length is shown by a solid line.
  • n is set, when the first feedback control is performed at each wavelength, a new time-SOA drive current waveform 813 (first update) is obtained. That is, the SOA drive current value is updated by Equation (4) for all wavelengths, and the next second wavelength sweep is performed while the SOA is driven by the time waveform of the time-SOA drive current 813.
  • the optical output power data of the second new wavelengths ⁇ 0 , ⁇ 1 ,..., ⁇ N ⁇ 1 can be acquired from the photodetector 723.
  • a time waveform of the updated time-SOA drive current 814 (second update) is obtained.
  • An updated time-SOA drive current 814 is used for the third wavelength sweep.
  • the time-SOA drive current waveform 814 is closer to the time-optimal SOA drive current waveform 815 at each wavelength than the time-SOA drive current waveform 813.
  • the current time waveform 813 and the current time waveform 814 match in the wavelength band near the center. That since the ⁇ condition D is met, updating the second I n were required
  • the SOA drive current-time waveform is obtained by performing the SOA drive current optimization process of FIG. 13 at each wavelength, for example, sample points 817 ⁇ 818 ⁇ 819 at the sweep start wavelength ⁇ 0 . Then, it operates so as to approach the time-optimal SOA drive current waveform 815. Therefore, when paying attention to a certain wavelength, the process of FIG. 14 for controlling the SOA drive current performs discrete control in time in conjunction with or in synchronization with the wavelength sweep.
  • FIG. 15 is a control flowchart of the SOA drive current optimization process including all wavelengths.
  • initial data at each wavelength point is set, and in S906 to S919, the SOA current is updated (S907) and the monitor voltage is acquired (S908) in each wavelength.
  • S917 from S911 N number of wavelengths ⁇ 0, ⁇ 1, ..., with respect to lambda N-1, compared with the tolerance of the error signal E n, such as determining the updated value of the SOA current value I n
  • the above-described operations are sequentially performed.
  • element Cn in sequence C represents whether the need to update the control current value I n at a wavelength lambda n.
  • Cn Fault.
  • AND (C) represents the AND of all the array elements of the array C.
  • FIG. 16 is a diagram showing a time waveform of an actual SOA drive current and an optical output level in the wavelength swept light source of the second embodiment in comparison with the prior art.
  • the upper part of FIG. 16 shows the time-SOA current waveform of 950 when the SOA drive current is constant in the prior art 950 and 951 after executing the optimization process of the SOA drive current according to the present invention.
  • the time-light output level (the output signal of the photodetector) is shown for the corresponding prior art 952 and for the case 953 after executing the SOA drive current optimization process.
  • the SOA current injection is only in the sweep period from the short wavelength to the long wavelength of the wavelength sweep, and as can be seen from FIG. 16, the optical output level is substantially constant with respect to the wavelength change (time change) in the wavelength sweep period. ing.
  • FIG. 17 is a diagram showing an interference fringe waveform of the wavelength swept light of the wavelength swept light source of the second embodiment in comparison with the prior art.
  • the interference fringe waveform obtained by the Michelson interferometer with the optical path length difference of 1 mm with the normalized amplitude is obtained in 954 when the SOA drive current is constant in the prior art, and in 955 after the SOA drive current optimization process is executed according to the present invention. , Respectively.
  • the amplitude of the fringe waveform 955 in FIG. 17 is also constant with respect to time compared to 954 when the SOA drive current is constant.
  • FIG. 18 is a diagram showing the relationship between the value of the optical output power made constant by controlling the SOA current and the coherence length.
  • the value of the optical output power controlled to be constant by the wavelength swept light source of the second embodiment is changed to obtain the coherence length for each optical output power.
  • the coherence length is maximized around 0.35 mW.
  • the coherence length at this time is about 24 mm over the wavelength sweep width of 80 nm, and is about twice the maximum coherence length of 12 mm when the SOA is operated with a constant SOA drive current of the prior art. The value of was obtained.
  • the wavelength with the maximum coherence length has been limited to a narrow range.
  • the wavelength swept light source of the present embodiment can obtain a coherence length that is twice that of the prior art.
  • the bisection method was used as the method for generating the time-SOA drive current waveform in the calculator 725.
  • the optical gain control means (723 to 726) performs temporally discrete control in conjunction with the wavelength sweep, and with respect to the difference between the output light level and the predetermined output light level.
  • the control amount (SOA drive current control step) for changing the optical gain is gradually reduced by half.
  • This bisection method is a simple feedback control method suitable when the control target amount (light output level) is relatively stable with respect to the environmental change of the wavelength swept light source.
  • the computing unit 725 can solve the problem by performing PID control or the like independently for each wavelength ⁇ of the sweep wavelength band instead of the bisection method.
  • the optical gain control means (723 to 726) performs temporally discrete control in conjunction with the wavelength sweep, and with respect to the difference between the output light level and the predetermined output light level.
  • the control for changing the optical gain is performed by the difference, the integration of the difference, and the differential of the difference.
  • N optical output power data V 0 , V 1 ,..., V at all wavelengths during one wavelength sweep period. acquires the n-1, by performing the updating of the SOA driving current I n, in the next wavelength sweep period, it is possible to take the n optical output power data of the new full wavelength.
  • executing an update 10 to 100 times about SOA driving current I n it is possible to reach a target light output level.
  • feedback control may be performed for each predetermined number of sweeps in consideration of the arithmetic processing capability and processing time in the arithmetic unit 725.
  • Repetition frequency of the sweep at the wavelength swept light source because it may also reach over 200kHz, if the actual utilization modes of the wavelength-swept light source, 10 to 100 times about SOA driving current I n a very short time updated in Only. Therefore, it is not necessary to continuously perform feedback control every time the light source wavelength sweep is started.
  • the steps S906 to S919 in the flow of FIG. 15 may be performed at intervals of 100, 1000, or 10,000 sweeps. This implementation interval may be a fixed interval or an irregular interval. For example, it is sufficient if the SOA drive current optimization process of FIGS. 14 and 15 is completed in the warm-up time of the wavelength swept light source device.
  • the series of flows in FIG. 15 may be stopped in an environment where the operation of the wavelength swept light source is stable.
  • the series of flows in FIG. 15 can be performed regularly or irregularly.
  • the SOA drive current optimization process of the present embodiment does not calculate the coherence length as in the first embodiment. Since the optical output level of the light source is monitored, it is directly fed back to the SOA current and controlled to a constant target optical output power value (V t ), so that the operation of the optimization process is much higher than that of the first embodiment. Easy to be.
  • the calculator 725 of the wavelength swept light source device 720 may use an existing CPU or the like.
  • the interferometer as in the first embodiment is used for the wavelength sweep in the normal operation of the wavelength swept light source apparatus. Is unnecessary. With a simpler component and processing, the coherence length of the wavelength swept output light during the wavelength sweep period can always be kept long as in the first embodiment.
  • the wavelength swept light source of the present embodiment has an optical gain value slightly higher than the optical gain in the laser oscillation threshold state, assuming that the relationship between the SOA current and the coherence length as shown in FIG.
  • the means (SOA) for changing the optical gain is controlled.
  • the target optical output power value (V t ) that maximizes the coherence length may be determined based on the actual SOA current-coherence length characteristic.
  • the PSF intensity in the SS-OCT apparatus can be increased, the image detection sensitivity of the SS-OCT apparatus can be increased, and the image quality can be improved.
  • a configuration and a control method for operating the wavelength swept light source by obtaining the SOA current that maintains the coherence length to the maximum, generating the SOA current waveform data, and the like have been described.
  • important elements that influence the coherence length of the output light are the optical deflector 106 shown in FIG. 1 and the optical deflector 703 shown in FIG.
  • an optical deflector suitable for the wavelength swept light source there is an optical deflector using an electro-optic crystal.
  • An optical deflector using an electro-optic crystal has no mechanical moving parts, and the electro-optic effect responds at a high frequency in the range of MHz to GHz, so that high-speed light deflection is possible.
  • various embodiments of an optical deflector that can be suitably used as the wavelength swept light source of the present invention will be shown.
  • KTN (KTa 1-x Nb x O 3 ) crystal or KLTN crystal (K 1-y Li y Ta 1-x Nb x O 3 ) (hereinafter collectively referred to as KTN crystal) is an electro-optic crystal at a low voltage. It is known as a substance having a large electro-optic effect whose refractive index changes greatly. Furthermore, when Ti and Cr materials are used as electrodes, charges can be injected into the crystal, and a high-speed and wide-angle optical deflector can be realized by utilizing an internal electric field generated by the injected charges. Accordingly, when general optical components such as lenses, prisms, and mirrors are used for applications that need to move at high speed, KTN crystals can be used instead of these optical components.
  • FIG. 19 shows a configuration of an optical deflector using a conventional KTN crystal. Electrodes 1102 and 1103 are formed on the two opposing surfaces of the KTN crystal 1101. A voltage is applied from the control voltage source 1104 between the two electrodes 1102 and 1103. Incident light 1105 is incident on a surface (xy plane) perpendicular to the surface on which the electrodes 1102 and 1103 of the KTN crystal 1101 are formed, and propagates in the KTN crystal 1101 in the z-axis direction. By being deflected in the KTN crystal 1101, outgoing light 1106 whose traveling direction is changed in the x-axis direction is obtained. A deflection angle corresponding to the applied voltage from the control voltage source 1104 is obtained.
  • a voltage corresponding to the application of the optical deflector is applied from the control voltage source 1104.
  • the waveform of the applied voltage is a sine wave, a sawtooth wave, or the like, and is applied according to the application of the optical deflector.
  • a voltage of about several hundred volts is applied to the KTN crystal 1101.
  • the charge density N in the KTN crystal approaches spatially uniform when the dielectric constant of the crystal is high, but when the dielectric constant is low, it accumulates only in the vicinity of the cathode when voltage is applied and is not uniform within the crystal. Distribution.
  • the charge density N is not spatially uniform and has a distribution that depends on the position in the crystal such as N (x, y, z), depending on the transmission position of the light deflected in the crystal,
  • the effect of the electro-optic crystal as a lens is different. That is, the lens effect varies depending on the instantaneous value of the drive voltage for light deflection.
  • variable focus lens whose focal length changes in synchronization with the voltage for light deflection is arranged on the emission side of the KTN crystal, and it does not depend on the voltage applied to the KTN crystal and has a constant light collection characteristic.
  • a method for obtaining light has been proposed. This method requires an additional variable focus lens that operates dynamically in addition to the original optical deflector, a power source for operating the variable focus lens, an electric circuit for synchronizing the optical deflector and the variable focus lens, etc. There is a problem that the configuration of the system becomes complicated.
  • the optical deflector suitable for the wavelength swept light source of the present embodiment operates to compensate for a spatial change in the lens effect in the optical deflector having the lens effect.
  • An optical deflector suitable for the wavelength swept light source of this embodiment is an optical deflector in which at least two electrodes are formed on opposite surfaces of an electro-optic crystal, and when a voltage is applied to the at least two electrodes, An optical deflector that causes a non-uniform refractive index distribution in the electro-optic crystal due to an optical effect, thereby bending a light path perpendicular to an electric field formed by the voltage, and an incident side and an exit of the electro-optic crystal An optical element disposed on at least one of the sides, wherein the light condensing characteristics are spatially distributed over a range in which the light is deflected.
  • an optical element in which the condensing characteristics and the optical path length are spatially distributed is arranged at the subsequent stage of the optical deflector, thereby compensating for the change in the deflection direction of the lens effect. Can do.
  • FIG. 20 is a diagram showing a configuration of an optical deflector suitable for a wavelength swept light source according to an embodiment of the present invention.
  • electrodes 1202 and 1203 are formed on two opposing surfaces of the electro-optic crystal 1201, respectively.
  • a voltage for injecting charges into the electro-optic crystal 1201 or deflecting light is applied from the control voltage source 1204 through this electrode pair.
  • KTN potassium tantalate niobate
  • KLTN a crystal doped with lithium
  • KLTN K 1-y Li y Ta 1- x Nb x O 3 (0 ⁇ x ⁇ 1, 0 ⁇ y ⁇ 1)
  • This electric charge distribution forms a non-uniform electric field distribution, and a refractive index distribution corresponding to the change in the electric field distribution is generated.
  • This non-uniform electric field distribution develops a refractive index gradient and bends the path of the light beam orthogonal to the refractive index gradient.
  • the incident light 1205 transmitted through the electro-optic crystal 1201 is deflected and condensed at the same time as the emitted light 1206 or the emitted light. 1208, and after the focus diverges.
  • the trap density N is not always spatially uniform.
  • the lens effect of the electro-optic crystal differs depending on the transmission position of the light deflected in the crystal. That is, the focal position formed by the electro-optic crystal differs depending on the instantaneous value of the voltage applied to the electro-optic crystal, which is an optical deflector.
  • the optical element 1209 having a spatially condensing characteristic that compensates for this variation can provide a deflection operation having a constant focal length regardless of the instantaneous value of the voltage applied to the optical deflector.
  • the beam shape of the emitted light at the time of deflection is the emitted light 1206 or the emitted light 1208 having different focal positions depending on the angle.
  • the emitted light 1208 has a focal position closer to the emission surface of the electro-optic crystal 1201 than the emitted light 1206. This is converted into outgoing light 1210 having the same focal position as the outgoing light 1206 by applying a wavefront change such that the outgoing light beam diverges by the optical element 1209 arranged at the rear stage of the optical deflector 1201.
  • the focal position of the emitted light 1206 is converted into emitted light 1211 similar to the emitted light 1208 by arranging the optical element 1209 so that the focal position of the emitted light 1206 is close to the emission surface of the electro-optic crystal 1201.
  • the optical element 1209 that corrects the divergence angle of the emitted light is spatially optical over the range in which the emitted light is deflected.
  • the characteristics are distributed.
  • the emitted light can have the same focal position at any deflection angle.
  • a laser printer or the like has an advantage that a beam waist can be formed at any position on the surface of the photoreceptor.
  • the outgoing light can be parallel light at any deflection angle.
  • the spectral characteristics of uniform resolution can be obtained without degrading the resolution over the entire spectral band. be able to.
  • FIG. 21 shows a configuration of an optical deflector suitable for a wavelength swept light source according to the fourth embodiment.
  • an optical deflector including an electro-optic crystal 1301 such as KTN or KLTN
  • electrodes 1302 and 1303 are formed on two opposing surfaces of the electro-optic crystal 1301, respectively.
  • Incident light 1305 incident on the electro-optic crystal 1301 passes through the cylindrical concave lens 1307 and then enters the electro-optic crystal 1301.
  • Incident light 1305, which is parallel light, is incident on the electro-optic crystal 1301 while being diffused by the cylindrical concave lens 1307.
  • Incident light 1305 is emitted as emitted light by the lens effect of the electro-optic crystal 1301.
  • the electro-optic crystal 1301 was a KTN crystal having a thickness of 1.2 mm (electrode spacing was also 1.2 mm), and a light beam having a beam diameter of 1.0 mm was incident as linearly polarized light parallel to the electric field.
  • the length of the electro-optic crystal 1301 is 4 mm.
  • Both the reflection film and the antireflection film are formed on the incident / exit end face of the electro-optic crystal 1301, and the distance by which the incident light is deflected by the folded optical path is 12 mm.
  • a sinusoidal AC voltage having an amplitude of ⁇ 360 V and a frequency of 200 kHz is superimposed on the electro-optic crystal 1301 and superimposed on a DC bias voltage of ⁇ 240 V by a control voltage source 1304.
  • the voltage is applied through electrodes 1302 and 1303 made of metal deposited on the crystal surface.
  • a deflection angle of 110 mrad (about 6.5 °) is obtained.
  • the trap density within the crystal is not spatially uniform. Since the trap density in the crystal is not uniform, the effect of the electro-optic crystal as a lens varies depending on the transmission position of light deflected in the crystal.
  • the prism 1309 is not provided, the beam when the outgoing light is deflected has a focal position that differs depending on the deflection angle, and is shaped like the outgoing lights 1306 and 1308.
  • Example 4 a prism 1309 having an angle of 30 ° between the entrance surface and the exit surface is disposed behind the electro-optic crystal 1301 using SF11 glass as a glass material. Since the prism 1309 enlarges the beam diameter of the transmitted light, the spread of the beam after transmission is suppressed by diffraction compared to before entering the prism 1309. The enlargement ratio of the beam diameter by the prism 1309 depends on the incident angle.
  • FIG. 22 shows the incident angle dependence of the beam diameter magnification of the prism in Example 4.
  • Example 4 when the focal position is a deflection angle (voltage) close to the emission surface of the electro-optic crystal 1301 in the deflected state (in the case of the emitted light 1308 in FIG. 21), the incident angle to the prism 1309 is When the deflection angle (voltage) is large and the focal position is far (in the case of the emitted light 1306 in FIG. 21), the incident angle to the prism is arranged to be small. Further, the prism 1309 is arranged such that the apex angle, that is, the direction in which the width of the prism is reduced, is on the side of the deflection angle where the focal position is far.
  • the outgoing light 1308 with a deflection angle (voltage) close to the focal position has an effect that the beam is further expanded and the spread due to diffraction is suppressed compared to the outgoing light 1306 with a deflection angle (voltage) having a large focal length.
  • the outgoing light 1306 becomes outgoing light 1311 by the prism 1309.
  • the prism 1309 whose light collection characteristics are spatially distributed over the range in which the emitted light is deflected suppresses the angle dependency of the focal position, and the AC voltage instantaneously applied from the control voltage source to the optical deflector. Regardless of the value, a deflection operation having a constant focal position is obtained.
  • FIG. 23 shows a configuration of an optical deflector suitable for the wavelength swept light source according to the fifth embodiment.
  • an optical deflector including an electro-optic crystal 1401 such as KTN or KLTN, electrodes 1402 and 1403 are formed on two opposing surfaces of the electro-optic crystal 1401.
  • Incident light 1405 incident on the electro-optic crystal 1401 passes through the cylindrical concave lens 1407 and then enters the electro-optic crystal 1401.
  • the shape and driving conditions of the electro-optic crystal 1401 are the same as those in the fourth embodiment.
  • the light collection effect depending on the deflection angle expressed by the optical deflector is corrected.
  • the free-form surface mirror 1409 is a concave mirror whose curvature radius in the deflection direction changes along the deflection direction, and its curvature decreases from the ground-side electrode 1402 to the control voltage source-side electrode 1403 of the optical deflector. It has become.
  • the curvature of the free-form surface mirror 1409 is small in the radius of curvature of the surface of the deflected outgoing light that is irradiated with the outgoing light when the focal position is at a deflection angle (voltage) close to the outgoing surface of the electro-optic crystal 1401.
  • the radius of curvature of the surface irradiated with the emitted light when the focal position is at a far deflection angle (voltage) is large.
  • FIG. 24 shows the applied voltage dependence of the focal length by only the electro-optic crystal 1401 in the fifth embodiment
  • FIG. 25 shows the relationship between the focal length of the free-form surface mirror and the applied voltage in the fifth embodiment.
  • the curvature of the concave surface irradiated with the emitted light is increased when the focal length is long (the focal position is far from the exit surface), and the focal length of the free-form mirror is shortened to change the focal length with respect to the applied voltage. To compensate. In this way, it can be seen that the curvature of the free-form surface mirror 1409 is set so as to compensate for the variation in the focal length with respect to the applied voltage.
  • the focal length of the free-form surface mirror 1409 is expressed as f (x) as a function of the position, the distance from the focal position of light incident on the free-form surface mirror 1409 at each incident position to the free-form surface mirror 1409 is always f.
  • a free-form surface such as (x) is formed.
  • the light reflected by the free-form surface mirror 1409 becomes parallel light at any deflection angle.
  • parallel light is generated.
  • parallel light is not necessarily important, and it is essential to keep a constant beam diffusion angle within a range in which outgoing light is deflected. It is.
  • Example 5 since a negative DC bias voltage is applied as the applied voltage from the control voltage source 1404, the charge density in the electro-optic crystal 1401 increases near the electrode 403 on the control voltage source side.
  • the relationship between the change in the curvature of the free-form curved mirror 1409 and the electrodes is as described above.
  • a positive DC bias voltage is applied as the applied voltage, the charge density in the vicinity of the ground-side electrode 1402 in the electro-optic crystal 1401 increases. Therefore, the curvature of the free-form surface mirror 1409 may be arranged so as to decrease on the electrode 1402 side.
  • FIG. 26 shows a configuration of an optical deflector suitable for the wavelength swept light source according to the sixth embodiment.
  • an optical deflector including an electro-optic crystal 1501 such as KTN or KLTN
  • electrodes 1502 and 1503 are formed on two opposing surfaces of the electro-optic crystal 1501, respectively.
  • Incident light 1505 incident on the electro-optic crystal 1501 passes through the cylindrical concave lens 1507 and then enters the electro-optic crystal 1501.
  • the shape and driving conditions of the electro-optic crystal 501 are the same as those in the fourth embodiment.
  • Example 6 the condensing effect depending on the deflection angle expressed by the optical deflector is corrected by using a transmissive optical element 1509 in the subsequent stage of the electro-optic crystal 1501.
  • the optical element 1509 is made of SF11 glass as a glass material, a free-form surface is formed on the exit surface, the curvature of the exit surface changes according to the deflection angle of the exit light, and the light is condensed spatially over the range in which the exit light is deflected. The characteristics are distributed.
  • the transmission type optical element 1509 is used in the sixth embodiment, an optical element having a refractive index distribution inside the element may be used so that the same light collecting effect can be obtained.
  • the optical element since the optical element has a flat surface, the optical element can be easily shaped.
  • Such an optical element can provide a condensing characteristic with high spatial resolution as compared with the optical element 1509 that provides the condensing characteristic by the shape of the exit surface by adjusting the internal refractive index distribution.
  • an optical element such as a mirror having a free-form surface is arranged on the emission side of the electro-optic crystal, so that the beam is diffused within the range of deflection of the emitted light of the optical deflector.
  • the corner can be kept constant.
  • FIG. 27 is a diagram illustrating a configuration of an optical deflector suitable for the wavelength swept light source according to the seventh embodiment.
  • an optical deflector including an electro-optic crystal 1601 such as KTN or KLTN
  • electrodes 1602 and 1603 are formed on two opposing surfaces of the electro-optic crystal 1601, respectively.
  • a cylindrical concave lens 1607 is arranged on the emission side of the electro-optic crystal 1601 so that the emitted light 1608 is parallel light. That is, the focal length of the cylindrical concave lens 1607 is that the beam when the focal length is large within the deflection range of the emitted light by the optical deflector (the focal position is a voltage far from the emission surface of the electro-optic crystal 1601) is parallel light. Choose to be.
  • the beam condensing effect expressed by the optical deflector depends on the applied voltage and the deflection angle corresponding to the applied voltage. Therefore, at some deflection angles, the beam condensing characteristic is canceled satisfactorily, but there exists a beam having a deflection angle that cannot completely cancel the condensing characteristic.
  • a spatial phase modulator 1609 is further arranged on the exit side of the cylindrical concave lens 1607.
  • the spatial phase modulator 1609 is an optical element that can independently modulate the phase of each location in the cross section of transmitted light or reflected light.
  • the light deflected by the electro-optic crystal 1601 enters the spatial phase modulator 1609.
  • the amount of phase modulation is continuously controlled so as to greatly delay the phase from the center toward the periphery.
  • the amount of phase modulation is such that, of the deflected outgoing light, the surface irradiated with the outgoing light when the focal position is close to the outgoing surface of the electro-optic crystal 1601 (voltage) is large, and the deflection angle with a far focal position ( The surface irradiated with the emitted light at the time of voltage is small.
  • the output light 1610 having a deflection angle made parallel by the cylindrical concave lens 1607 is not changed in phase by the spatial phase modulator 1609, and the light condensing effect is not canceled by the cylindrical concave lens 1607.
  • a spatially modulated phase modulation is applied to the angular outgoing light 1606 so as to cancel the light collection effect.
  • a holographic optical element or a diffractive optical element having a spatial phase modulation distribution may be used.
  • FIG. 28 shows the result of obtaining the relationship between s and s ′′ when the focal length is changed by this equation.
  • the ratio of the beam waist diameter between the incident light and the outgoing light, that is, the magnification m is expressed as the following formula.
  • FIG. 29 shows the relationship between the beam waist position and the beam diameter magnification when the focal length is changed.
  • the beam waist diameter also changes in accordance with the change in focal length.
  • the condensing characteristic (focal length) of the optical deflector changes according to the deflection angle
  • the beam waist position and beam waist diameter of the emitted light change according to the deflection angle as described above.
  • the characteristics of the compensation optical elements shown in the fourth to seventh embodiments are also changed in accordance with the deflection angle to suppress variations in the beam waist position and the beam diameter.
  • a prism is used as the compensation optical element, which has the same configuration as that of the optical deflector of the fourth embodiment (FIG. 21).
  • the optical path length is longer than in the air, so that the distance between the optical deflector and the beam waist position of the emitted light can be adjusted by the prism. That is, as the beam waist position changes depending on the deflection angle, the optical path length in the prism is set longer when the beam waist position is far from the optical deflector, and when the beam waist position is closer to the optical deflector, the prism is set.
  • the inner optical path length may be set short.
  • the prism apex angle is set so that the prism width is large when the beam waist position is a far deflection angle and the prism width is small when the beam waist position is a near deflection angle. On the deflection angle side where the beam waist position is close. Further, since the magnification of the beam diameter varies depending on the incident angle and the exit angle to the prism, the beam waist diameter can be compensated. Furthermore, the beam waist position and the beam waist diameter can be compensated simultaneously by appropriately determining the apex angle of the prism and the incident angle to the prism.
  • the optical path length and the beam diameter magnification can be set independently to some extent, so that the accuracy of compensation can be improved.
  • the beam waist position and the beam waist diameter of the emitted light can be changed by controlling the optical path length and the lens power until it enters the mirror.
  • the shape of the free-form surface mirror may be determined so as to compensate for a change in either the beam waist position or the beam waist diameter, or the shape of the free-form surface mirror may be determined so as to suppress both changes.
  • the focal length of the free-form surface mirror is shortened with respect to the deflection angle where the beam waist position is far from the optical deflector, and the focal length of the free-form surface mirror is increased with respect to the deflection angle where the beam waist position is close to the optical deflector, good.
  • the beam waist position and beam waist diameter of the emitted light are adjusted by changing the phase modulation amount and refractive index distribution according to the deflection angle. The change depending on the deflection angle can be suppressed.
  • the optical path length and the beam waist diameter can be set independently to some extent, so that the accuracy of compensation can be improved.
  • an optical element with spatially distributed optical characteristics is arranged on the incident side of the optical deflector, depending on the distance between the beam waist position and the optical deflector. Compensation.
  • variation in one or both of the beam waist position and the beam waist diameter of the emitted light can be suppressed within the range of deflection of the emitted light of the optical deflector.
  • the adaptive optics system can be optimized according to the beam waist position and beam waist diameter of incident light.
  • the optical deflector using KTN or KLTN has been described.
  • an optical deflector using a crystal having an electro-optic effect other than these may be used.
  • the optical deflectors of Examples 4 to 8 can be suitably used for the wavelength swept light sources of the first and second embodiments, respectively.
  • the present invention can provide a wavelength swept light source capable of maintaining a long coherence length over a wide swept wavelength band.
  • the image quality of the SS-OCT apparatus can be improved.
  • the present invention can be generally used for an imaging apparatus.
  • it can be used as a wavelength swept light source for an SS-OCT apparatus.

Abstract

SS-OCT装置では、波長掃引光源からの光のコヒーレンス長が長くなるほど干渉信号の強度は上がる。ゲイン媒体に時間的に一定の電流注入を行っていた従来技術の波長掃引光源では、波長掃引をしている全波長帯域に渡って、コヒーレンス長を長い状態に維持できていなかった。SS-OCT装置の出力イメージの品質に直接影響を与えるPSF強度が下がってしまう問題があった。本発明の波長掃引光源は、掃引波長帯域全域に渡り、出力光のコヒーレンス長を長い状態に維持する。波長掃引に対応して、ゲイン媒体の利得を時間と共に変化させ、コヒーレンス長最大の状態に維持される。掃引波長帯域全域に渡り、ゲイン媒体の利得を発振閾値の近傍の状態に維持しゲイン媒体の非飽和利得帯域幅を狭窄化する。コヒーレンス長を長い状態に維持しながら駆動するSOA電流波形データ取得方法、また新規なコヒーレンス長の測定方法、さらに波長掃引光源に好適な光偏向器も開示される。

Description

波長掃引光源、波長掃引光源のための駆動データ作成方法および光偏向器
 本発明は、波長掃引光源に関する。より詳細には、光干渉断層撮影装置で使用できる波長掃引光源に関する。
 波長掃引光源は、光を利用したイメージング技術を用いた光学機器や電子機器に幅広く利用されている。例えばカメラやプリンタ、ファクシミリなどの民生用の電子機器に最もよく利用されており、さらに医療分野にも広がっている。生体内部の断層を非侵襲的にイメージングするために、低干渉光を利用して深さ方向の情報を高分解能で取得することができる光干渉断層撮影(OCT:Optical Coherence Tomography)が注目されてきた。近年では、生体内での分子レベルおよび細胞レベルの生物学的事象を解明するためのモレキュラー(分子)イメージングの分野にも、OCTが利用されている。波長掃引光源は、後述するようにOCT技術を支える重要な基幹要素であり、OCT装置の性能にも大きな影響を与える。
 上述のモレキュラーイメージングでは、光情報を利用してターゲット分子を高感度に検出する方法が主力となっており、OCT装置が利用されている。コヒーレンス干渉を利用して断層イメージを取得するOCT装置には、時間領域OCT(TD-OCT:Time Domain Optical Coherence Tomography)と、フーリエ領域OCT(FD-OCT:Fourier Domain Optical Coherence Tomography)がある。FD-OCTは、さらにスペクトル領域OCT(SD-OCT:Spectral Domain Optical Coherence Tomography)と波長掃引OCT(SS-OCT:Swept-Source Optical Coherence Tomography)に分類される。波長掃引光源を用いたSS-OCTは、高速な応答性の点で特に優れている。
 SS-OCTの応答性は光源の性能によって大きく律速されるため、種々の方式の高速広帯域光源の開発が加速している。また掃引の態様の観点からは、波長掃引光源には、掃引される波数が時間と共に線形的に変化する波数リニアな光源のように、応用分野に応じた動作が要求される。
 特許文献1には、光ファイバループによる波長掃引光源の構成例が示されている。特許文献1の図1を参照すれば、波長掃引光源は、光ファイバループ(11)にゲイン媒体(12)とサーキュレータ(13)を設けている。サーキュレータ(13)から、コリメートレンズ(22)を通して光が空間に取出され、その光がバンドパスフィルタ(24)を通してミラー(25)に到達する。ミラーからの反射光を、バンドパスフィルタ(24)を通してコリメートレンズに再度入射させることで、共振器の往復光路が構成される。バンドパスフィルタ(24)としては、例えば光干渉型の誘電体多層膜フィルタを利用可能であって、フィルタへの光の入射角を変えることによって透過波長を変えることができる。バンドパスフィルタの光軸に対する角度を、ガルバノメータ(26)を用いて高速に変えることによって、高速に波長掃引をすることができる光源を実現できる。
特開2005-347668号公報 特開2011-186218号公報 特開2015-142111号公報
Yoshiaki Yasuno他、"Three-dimensional and high-speed swept-source optical coherence tomography for in vivo investigation of human anterior eye segments,2005年、"OPTICS EXPRESS,Vol.13,No.26,pp.10652-10664
 特許文献1に示されたような波長掃引光源では、通常、周回利得(ゲイン)や発振閾値は波長によって異なる。レーザでは、周回利得スペクトルが発振閾値よりも大きくなる(周回利得スペクトル>発振閾値となる)波長帯域内で発振する。この波長帯域は、発振閾値(=損失)よりも大きくなる周回利得の帯域幅を意味する。以下、この波長帯域幅を非飽和利得帯域幅と呼ぶこととする。
 SS-OCT用の波長掃引光源は、掃引している全波長においてコヒーレンス長が長い状態に維持することが望ましい。これを実現するためには、掃引している全波長に渡って、非飽和利得帯域幅(レーザ発振可能な波長幅)が狭いことが重要である。
 非飽和利得帯域幅は周回利得スペクトルが発振閾値を超える波長幅であるので、全掃引波長に対して、周回利得スペクトルのピークが発振スペクトルに一致するような動作状態にできれば良い。しかしながら、従来技術の波長掃引光源では、ゲイン媒体への注入電流量が時間的に一定であったため、このような動作状態の実現は難しかった。レーザ発振器では、ゲイン媒体への注入電流量が周回利得スペクトルの形状、ピーク値を決める。このため、注入電流量が時間的に一定であった従来技術の波長掃引光源では、掃引する全波長に渡って、周回利得ピーク値が発振閾値と一致するとは限らないためである。この理由は、さらに以下のように説明できる。
 周回利得および発振閾値は、ゲイン媒体の利得、レーザ共振器内のさまざまな光損失、バンドパスフィルタのフィルタリング特性、および、各波長に対して発振に掛けられる時間等から決まる。例えば、特許文献1に開示された波長掃引光源等では、各波長の発振に掛けられる時間、および、バンドパスフィルタのフィルタリング特性は波長によらず一定であるが、ゲイン媒体の利得は掃引波長の中心波長から離れるに従って小さくなっている。したがって、中心波長付近の発振閾値は小さく、中心波長から離れるにしたがって発振閾値は大きくなる。
 従来技術の波長掃引光源では、掃引中のどの波長に対しても、ゲイン媒体への注入電流量を一定としていたので、広帯域に必要な掃引波長幅(例えばSS-OCTでは100nm程度)を確保するためには、大きな発振閾値を持つ波長に合わせてゲイン媒体への電流注入量を決める必要がある。このようにゲイン媒体へ一定電流を注入した場合、発振閾値が小さい波長においては、閾値を超える周回利得の波長幅(非飽和利得帯域幅)が大きい状態となっていた。特許文献1に開示された波長掃引光源では、掃引波長の最短波長および最長波長において発振可能なゲイン媒体への注入電流量としていたため、最短波長および最長波長付近ではコヒーレンス長が長い状態となっている。しかしながら、中心波長に近付くにつれて注入電流量が過剰となって、周回利得のピーク値が発振閾値を超える程度が著しくなる。このため、掃引波長帯域の概ね中央部では、非飽和利得帯域幅が広くなり発振線幅が広がってコヒーレンス長が短い状態となっていた。
 SS-OCT装置では、波長掃引光源からの光のコヒーレンス長が長くなるほど干渉信号の強度は上がることが知られている。干渉信号を周波数変換して得られるOCT画像の1画素すなわち点広がり関数(PSF;point spread function)の強度は、干渉信号のエネルギー(時間積分値)に比例する。したがって、PSFの強度を上げるためには、波長掃引をしている全波長帯域において干渉信号全体を強くする必要があり、掃引する全波長帯域においてコヒーレンス長も長い状態でなければならない。ゲイン媒体に時間的に一定の電流注入を行っていた従来技術の波長掃引光源では、波長掃引をしている全波長帯域に渡って、コヒーレンス長を長い状態に維持できていなかった。その結果、SS-OCT装置の出力イメージの品質に直接影響を与えるPSF強度が下がってしまう問題があった。
 SS-OCT装置においてPSF強度が下がると、検出するイメージのPSFピークとシステムノイズとの差が小さくなるため、対象物からの微小な検出信号がシステムノイズの中に埋もれてしまう。結果として、SS-OCT装置のイメージ品質が低下し、イメージの検出感度が落ちてしまう問題があった。
 本発明は、このような問題に鑑みてなされたものであって、その目的とするところは、広い掃引波長帯域全体に渡って、コヒーレンス長が長い状態を維持することができる光源を提供することにある。また、本発明の波長掃引光源を最適に動作させるための半導体光増幅器(SOA)電流波形データ(光利得制御データ)を作成する方法、コヒーレンス長を求める方法も開示される。また、本発明の波長掃引光源を実施するのに好適な光偏向器についても開示される。
 本発明の1つの実施態様は、発振波長が連続的に変化する光を出力する波長掃引光源において、光を増幅するゲイン媒体、発振波長を変更する手段、および前記ゲイン媒体の光利得を変更する手段を含むレーザ共振器と、波長制御信号によって前記発振波長を変更する手段を制御して、前記発振波長を掃引する波長制御手段と、波長掃引の全波長帯域にわたって、前記レーザ共振器内の光利得が、レーザ発振閾値の状態の光利得よりもわずかに高い光利得値となるように、前記光利得を変更する手段を制御する光利得制御手段とを備えたことを特徴とする波長掃引光源である。
 別の実施態様によれば、前記光利得制御手段は、前記波長制御手段から得られた掃引開始からの複数の時間と、前記複数の時間の各々において波長掃引光のコヒーレンス長が最長となる光利得とを組み合わせた光利得制御データに基づいて、光利得制御信号によって前記光利得を変更する手段を制御するよう構成され得る。
 好ましくは前記発振波長を変更する手段は、回折格子と、前記回折格子からの発振光が入射し、当該入射した前記発振光を垂直に反射するミラーと、前記回折格子への前記発振光の入射角を変更する光偏向器とから構成され得る。
 またさらに別の実施態様によれば、前記光利得制御手段は、前記レーザ共振器からの出力光レベルを検出し、前記検出された出力光レベルに基づいて、前記出力光レベルが、所定の出力光レベル値になるような光利得制御信号を生成し、当該光利得制御信号によって前記光利得を変更する手段を制御するよう構成され得る。
 好ましくは、前記発振波長を変更する手段は、回折格子と、前記回折格子への発振光の入射角を変更する光偏向器とから構成され得る。
 本発明の他の実施態様は、発振波長が連続的に変化する光を出力する波長掃引光源において、光を増幅するゲイン媒体、発振波長を変更する手段、および前記ゲイン媒体の光利得を変更する手段を含むレーザ共振器と、波長制御信号によって前記発振波長を変更する手段を制御して、前記発振波長を掃引する波長制御手段と、前記波長制御手段から得られた掃引開始からの複数の時間と、前記複数の時間の各々において前記波長掃引光のコヒーレンス長が最長となる光利得とを組み合わせた光利得制御データに基づいて、光利得制御信号によって前記光利得を変更する手段とを備えたことを特徴とする波長掃引光源である。
 本発明のさらに異なる実施態様によれば、上述の光偏向器は、電気光学結晶の対向する面に少なくとも2つの電極が形成された光偏向器であって、前記少なくとも2つの電極に電圧を印加すると、電気光学効果により前記電気光学結晶内に一様でない屈折率分布が生じることにより、前記電圧により形成された電界に直交する光の進路を屈曲させる光偏向器と、前記電気光学結晶の入射側と出射側の少なくとも一方に配置された光学素子であって、前記光が偏向する範囲にわたって空間的に集光特性が分布している光学素子とを備えたものとすることができる。
 好ましくは前記光偏向器は、偏向角により集光特性が変化し、前記光学素子は、前記光が偏向する範囲にわたって出射側のビームウェスト位置とビームウェスト径の少なくとも一方の変化を補償するように空間的に光学特性が分布しているものとすることができる。
 以上説明したように、本発明の波長掃引光源およびその駆動データの作成方法、光偏向器によって、波長掃引期間中の出力光のコヒーレンス長を常に長い状態に維持できる。
本発明の第1の実施形態の波長掃引光源の構成を示す図である。 波長掃引光源におけるSOA電流の動作を説明する図である。 本発明においてSOA電流波形データの生成装置を示す図である。 本発明でSOA電流波形データ生成手順のフローを示す図である。 本発明で取得されるコヒーレンス長データ全体を示した表である。 本発明におけるコヒーレンス長の測定装置の構成を示す図である。 本発明におけるコヒーレンス長測定手順のフローチャートである。 本発明の波長掃引光源で、等時間間隔取得したSOA電流の説明図である。 図6に示した干渉計の具体的な構成例を示した図である。 図6に示した干渉計の具体的な別の構成例を示した図である。 図6に示した干渉計の具体的な別の構成例を示した図である。 図6に示した干渉計の具体的な別の構成例を示した図である。 SOA注入電流と光出力パワー、コヒーレンス長の関係を示した図である。 本発明の第2の実施形態の波長掃引光源の構成を示す図である。 制御回路を含む波長掃引光源装置の構成を示す図である。 SOA注入電流と光出力レベルとの関係を説明する図である。 SOA注入電流最適化を波長掃引も含めて説明する図である。 全波長を含むSOA駆動電流の最適化フローチャートである。 第2の実施形態の波長掃引光源の時間SOA駆動電流波形他の図である。 第2の実施形態の波長掃引光源の干渉フリンジ波形の図である。 光出力パワー値とコヒーレンス長との関係を示す図である。 通常のKTN結晶を用いた光偏向器の構成を示した図である。 本発明の第3の実施形態の光偏向器の構成を示した図である。 実施例4の波長掃引光源に好適な光偏向器の構成を示した図である。 実施例4のプリズムのビーム径倍率の入射角依存性を示した図である。 実施例5の波長掃引光源に好適な光偏向器の構成を示した図である。 実施例5の焦点距離の印加電圧依存性を示した図である。 実施例5の自由曲面ミラー焦点距離と印加電圧の関係を示した図である。 実施例6の波長掃引光源に好適な光偏向器の構成を示した図である。 実施例7の波長掃引光源に好適な光偏向器の構成を示した図である。 光学素子の焦点距離と、ガウシアンビーム入射時のビームウェスト位置との関係を示した図である。 焦点距離を変化時のビームウェスト位置とビーム径倍率の関係を示した図である。
 本発明の1つの実施態様は、発振波長が連続的に変化する光を出力する波長掃引光源において、光を増幅するゲイン媒体、発振波長を変更する手段、および前記ゲイン媒体の光利得を変更する手段を含むレーザ共振器と、波長制御信号によって前記発振波長を変更する手段を制御して、前記発振波長を掃引する波長制御手段と、波長掃引の全波長帯域にわたって、前記レーザ共振器内の光利得が、レーザ発振閾値の状態の光利得よりもわずかに高い光利得値となるように、前記光利得を変更する手段を制御する光利得制御手段とを備えた波長掃引光源である。
 光利得制御手段については、異なる方法による構成例(第1の実施形態および第2の実施形態)が開示される。1つは、コヒーレンス長を最大に維持するSOA電流を予め求めて、SOA電流波形データによって光源を動作させるものである。もう1つは、レーザ共振器からの出力光レベルを検出し、検出された出力光レベルに基づいて、出力光レベルが、所定の出力光レベル値になるように光利得(SOA電流)を制御するものである。
 さらに、発明の波長掃引光源を実施するにあたって好適な、光偏向器の様々な構成例についても開示される。
[第1の実施形態]
 本発明の第1の実施形態の波長掃引光源は、掃引波長帯域の全域に渡って、出力光のコヒーレンス長を長い状態に維持できる。波長掃引に対応して、ゲイン媒体の利得を時間と共に変化させることで、コヒーレンス長は最大の状態に維持される。すなわち、掃引波長帯域の全帯域に渡って、各波長においてゲイン媒体の利得を発振閾値の近傍の状態に維持し、ゲイン媒体の非飽和利得帯域幅を狭窄化するように動作する。すなわち、レーザ共振器内の利得(光利得)が、レーザ発振閾値の状態の利得(発振閾値利得または発振利得)よりもわずかに高い利得値となるように、光利得を制御する。従来技術の波長掃引光源において過剰な利得を持っていた掃引波長帯域の概ね中央部で、より周回利得を下げるようにレーザの利得制御手段を制御することで、出力発振光の発振線幅を常に狭い状態に維持する。
 これに限定されないが、例えばレーザ共振器におけるゲイン媒体に注入する半導体光増幅器(SOA:Semiconductor Optical Amplifier)に流すSOA電流によって、ゲイン媒体の周回利得を変化させることができる。本発明の波長掃引光源におけるコヒーレンス長を最大に維持するSOA電流は、SOA電流波形データに基づいて生成される。SOA電流波形データを得るためには、異なるSOA電流に対して、波長掃引光源を掃引した状態で掃引時間情報と対応付けながら、光源出力光の干渉光から各波長に対応したコヒーレンス長を求める。各波長において、異なるSOA電流の中から、最大コヒーレンス長を与えるSOA電流が選択される。波長掃引光源を通常の掃引動作状態に置いたままで、波長掃引光から得られた複数の干渉光に基づいてコヒーレンス長を求め、コヒーレンス長を最大に維持するSOA電流が求められる。コヒーレンス長を最大に維持するSOA電流は、任意の電流波形を出力可能な電流源から供給され、波長掃引光源は、記憶手段(メモリ)に記憶されたSOA電流波形を実現するためのSOA電流波形データに基づいて動作する。レーザ共振器のゲイン媒体の利得の変更は、SOAの駆動電流以外の他の利得可変手段によっても実現できる。
 本発明は、波長掃引光源を、全掃引波長帯域において出力発振光のコヒーレンス長を長い状態に維持しながら駆動するためのSOA電流波形データを求める方法としても実施できる。また、SOA電流波形データを決定するための、本発明の波長掃引光源に適合した本発明に特有の新規なコヒーレンス長の測定方法についても開示される。以下、本発明の波長掃引光源およびその駆動データの作成方法、コヒーレンス長の測定方法について、各実施例と共に説明する。最初に第1の実施形態の波長掃引光源の実施例1として、本発明の波長掃引光源の構成および動作について説明する。
 図1は、本発明の波長掃引光源の構成例を示す図である。本発明の波長掃引光源100は、リットマン型のレーザ共振器101を持つものを例として説明される。後述するように、波長依存の周回利得を持つゲイン媒体を共振器内に有するレーザ発振器であれば、本発明はどのようなタイプのレーザ共振器にも適用可能であることに留意されたい。SOA102から出力された光は、レンズ105によって平行光にコリメートされて、光偏向器106を通して偏向される。光偏向器106は、後述する光偏向器ドライバからのドライブ信号115によって、例えば図1のように図面の上下方向に光路を偏向する。図1では、図の上側に最大に偏向した状態の光路121および下側に最大に偏向した状態の光路122の2つの場合を描いてある。図1の構成では、光路121の状態および光路122の状態が、それぞれ、最長波長の状態および最短波長の状態に対応する。
 光偏向器106によって偏向された光は、回折格子107で回折し、ミラー108に至る。ミラー108に垂直入射した波長の光のみが、上述の光路とは逆方向の光路を辿り、SOA102に戻る。ミラー108からSOA102に戻った光は、SOA102およびカップラ103の間にあるハーフミラー104で反射される。ミラー108およびハーフミラー104の間で往復光路が形成されて、レーザ発振が生じる。発振光の一部は、SOA102からハーフミラー104を透過して、カップラ103と光ファイバ120を介して、出力光としてレーザ共振器101の外部へ出力される。
 上述のように、ミラー108に垂直に入射し反射する波長成分のみで往復する光路が形成され発振状態となる構成となっているため、回折格子107への光の入射角によってレーザ発振波長を変えることができる。SOA102およびミラー108の間の構成は、例えば概ね二次関数形状のバンドパスフィルタを形成し、周回利得の非飽和利得帯域幅を決定する。
 光偏向器ドライバ113は、ドライブ信号115によって、光偏向器106から出射する光の偏向方向(偏向角)を制御する。ドライブ信号115が時間とともに変化するのに応じて、偏向方向も徐々に変化する。波長掃引光源では、発振光の波長が滑らかに一方向に増加または減少して波長掃引されるように、ドライブ信号は周期的に繰り返し変化する。通常、1回の波長掃引と、次の波長掃引との間には所定の間隔がある。すなわち、繰り返し期間の内の一部が実際の波長掃引に使用される。このようにレーザ共振器101では、ドライブ信号115によって時間的に連続に発振波長を変えて、波長掃引動作が実現される。
 光偏向器ドライバ113からは、ドライブ信号と共に、波長掃引に同期した同期信号114が出力される。例えば、短波長から長波長へ波長掃引する場合は、最短波長に相当する偏向角となるタイミングで同期信号114を任意波形電流源111へ出力する。
 任意波形電流源111は、図1には示していない関連付けられたメモリ(記憶手段)内に、波長掃引動作中に使用し、SOA102に印加するSOA電流116を生成するためのSOA電流波形データを備えている。このメモリは、電流波形データを他から読み出して一時的に記憶するものであっても、半永続的に記憶するものであっても良い。電流波形データの内容の詳細については、実施例2および実施例3で説明するが、SOA電流波形データを生成するための情報がメモリに記録される。
 メモリは、任意波形電流源111内に備えられていても良いし、任意波形電流源111外であって波長掃引光源100内のいずれかの場所に備えられていれば良い。任意波形電流源111は、光偏向器ドライバ113からの同期信号114と同期して、メモリに記録された電流波形データまたは情報から、SOA電流波形を生成し、SOA電流波形に従ったSOA電流をSOA102に流す。なお、SOA102にはSOAの温度制御を行う温度制御器が通常接続されているが、ここでは省略してある。これ以降に説明するSOAについても同様である。通常、波長掃引光源では、波長掃引動作の制御を始めとして、光源全体の制御のために、CPUを含む制御部もさらに備えている。制御部は通常メモリを含んでおり、メモリに記憶されまたは読み出されたデータを使用して動作する電子制御機器は、通常、制御部を持っている。したがって、制御部は図1には示されていないことに留意されたい。
 任意波形電流源111に関連付けられたメモリに記録される情報としては、時間-波形そのものを記録しても良いし、時間-波形を生成するための任意の形態のデータを記録しても良い。すなわち、その時間-波形のデータそのもの、時間-波形のデータおよびこれと組み合わせる係数、時間-波形のデータおよびこれと組み合わせるバイアス、または、時間-波形のデータ並びに係数およびバイアスを組み合わせたものを出力しても良い。
 さらに、任意波形電流源111に関連付けられたメモリに記録される情報としては、電流波形をフーリエ変換したデータとすることもできる。この場合も、フーリエ変換したデータを逆フーリエ変換したデータそのもの、逆フーリエ変換したデータおよびこれと組み合わせる係数、逆フーリエ変換したデータおよびこれと組み合わせるバイアス、または、逆フーリエ変換したデータ並びに係数およびバイアスを組み合わせたものを出力しても良い。
 本発明の波長掃引光源を動作させるとき、まず、光偏向器ドライバ113は光偏向器106を偏向動作させる。偏向動作を開始するとともに、任意波形電流源111は、SOA102に印加することになる電流波形データを、任意波形電流源111に関連付けられたメモリへ読み出す。光偏向器ドライバ113は、光偏向器106を偏向動作させる。任意波形電流源111は、メモリに記憶した電流波形データに基づいて、光偏向器ドライバ111からの同期信号114と同期しながら、電流波形データに対応したSOA電流116をSOA102に印加する。本発明の動作原理について、さらに詳細に説明する。
 既に述べたように、SS-OCT装置に使用する波長掃引光源では、全掃引波長帯域においてレーザ共振器からの発振光のコヒーレンス長を長い状態に維持して、発振線幅をできるだけ細くする必要がある。図1の本発明の波長掃引光源の構成では、従来技術と同様に、回折格子107およびミラー108で構成される波長フィルタによって、発振線幅の基本的な狭窄化動作が行われている。これに加えて本発明の波長掃引光源では、周回利得を掃引と同期して動的に調整し、波長掃引中のどの時間においても、非飽和利得帯域幅ができる限り狭くなるようにSOA102の動作を制御する。すなわち、刻々と変化する発振波長において、SOA102を、周回利得が発振閾値よりもわずかに大きい状態とし、非飽和利得帯域幅ができる限り狭くなるようにSOA電流116を設定する。これによって、1回の波長掃引期間中のどの時間においても、発振光のコヒーレンス長が長い状態に維持されるようにする。周回利得の大きさは、波長に依存して異なる。そこで、波長掃引中に波長ととともにSOAの光利得を動的に変えると、1回の波長掃引期間中のどの時間でもコヒーレンス長を長い状態に維持することができる。波長掃引光源では、波長掃引期間中における刻々の波長は、掃引開始からの時間、すなわち掃引時刻と1対1に対応している。したがって、波長に応じて、SOA電流116によって周回利得の設定をすることは、SOA電流116を波長掃引の開始からの時刻と関連付けて、時間―SOA電流波形をどのように決定するかに帰着することになる。本発明の波長掃引光源は、次の実施例2によって得られる時間―SOA電流波形の電流波形データを使用して、SOAを動作させることになる。
 図2は、本発明の波長掃引光源におけるSOA電流の動作を説明する図である。図2の上下のグラフでは、いずれも横軸が掃引開始時刻からの時間を表しており、上下のグラフ間で時間は同期している。図2の上側のグラフは1回の波長掃引期間における発振波長の時間変化を示している。波長掃引開始時点では発振波長は概ね1280nm(波数換算で概ね4.9E6rad/m)で、掃引終了時点で発振波長は概ね1400nm(波数換算で概ね4.5E6rad/m)となっており、掃引波長帯域幅は120nmに達している。図2の下側のグラフは、1回の波長掃引期間において、本発明に従って、SOA電流116を変化させた場合の、時間―SOA電流波形を示している。波長に応じて、SOA電流116を変化させていることがわかる。具体的には、掃引波長帯域の両端(掃引開始および掃引終了時間に対応)の少なくとも一方のSOA電流値よりも、掃引波長帯域の中央部におけるSOA電流を減らしている。これによって、掃引波長帯域の中央部におけるゲイン媒体の周回利得を動的に下げ、波長掃引中のどの時間でも、非飽和利得帯域幅ができる限り狭くなるようSOA102を制御している。図2の上側の図でも明らかなように、掃引時刻と波長との間は1対1に対応する関係となっており、掃引開始時からの時間(時刻)によって、瞬時の発振波長が特定できることがわかる。
 本実施例の図1の波長掃引光源100では、SOA電流を動的に時間変化させながらSOAを駆動するための、時間―SOA電流波形データが既に存在している状態で動作する。このSOA電流波形データを作成する方法および装置については、実施例2および実施例3においてさらに詳細に説明する。
 上述の図1の波長掃引光源100では、レーザ共振器101としてリットマン型外部共振器の構成を例として説明したが、本発明は周回利得を変更する手段を備えることができるものであれば、どのようなタイプのレーザ共振器を持つ波長掃引光源にも適用できる。すなわち、リトロー型外部共振器、ファイバーレーザ、垂直共振器面発光レーザ(VCSEL:Vertical Cavity Surface Emitting Laser)-MEMSレーザ等を備えた波長掃引光源にも適用できる。レーザ共振器が、波長を変更する手段(図1では光偏向器106、回折格子107およびミラー108)と、光利得を変更する手段(図1ではSOA102)とを備えているものであれば、本発明を適用できる。したがって、以下のように一般化が可能である。
 例えば、波長を変更する手段は、図1では光偏向器ドライバ113に相当する波長制御手段から出力される波長制御信号(偏向器ドライブ信号115)を受けて、レーザ発振波長を制御する。波長制御手段が時間とともに連続的に波長を変化させるように制御すれば、波長掃引が可能となる。既に述べたように、図1の本実施例の波長掃引光源の構成例では、光偏向器106の光偏向角度を連続的に変化させることによって発振波長の掃引動作が生じる。
 光偏向器としては、例えば、電気光学材料としてタンタル酸ニオブ酸カリウム(KTN)結晶を用いたKTN偏向器を利用することができる。KTNの他、KLTN(K1-yLiTa1-xNb(0<x<1、0<y<1))、LiNbO(以下、LNという)、LiTaO、LiIO、KNbO、KTiOPO、BaTiO、SrTiO、Ba1-xSrTiO(0<x<1)、Ba1-xSrNb(0<x<1)、Sr0.75Ba0.25Nb、Pb1-yLaTi1-xZr(0<x<1、0<y<1)、Pb(Mg1/3Nb2/3)O-PbTiO、KHPO、KDPO(Dは重水素)、(NH)HPO、BaB、LiB、CsLiB10、GaAs、CdTe、GaP、ZnS、ZnSe、ZnTe、CdS、CdSe、およびZnOの電気光学材料も利用できる。
 また、光利得を変更する手段は、図1では任意波形電流源111に相当する光利得制御手段から出力される光利得制御信号(SOA電流116)を受けて制御され、光利得が制御される。光利得制御手段が時間とともに光利得制御信号を変更すると、それに合わせて、光利得も時間とともに変化する。図1の例では、SOA電流116を時間変動させることによって、SOA102の光利得(ゲイン媒体の周回利得)がSOA電流116に合わせて時間とともに変化する。
 図1で示した波長掃引光源100におけるSOA電流波形データは、1掃引中のSOA電流116の時間推移のデータ(時間およびSOA電流を組み合わせたデータ)となる。波長を変更する手段および光利得を変更する手段からなる波長掃引光源としてより一般化すれば、SOA電流波形データは、1回の波長掃引中の光利得制御信号の時間推移のデータ(時間および光利得制御信号を組み合わせたデータ:光利得制御データ)の一例ということができる。図2において、下側のグラフの縦軸を光利得制御信号値と読み替えれば、光利得制御信号値が波長とともに動的に制御された、光利得制御データを求めることができる。
 したがって本発明の波長掃引光源は、発振波長が連続的に変化する波長掃引光を出力する波長掃引光源において、光を増幅するゲイン媒体、発振波長を変更する手段(106、107、108)、および前記ゲイン媒体の光利得を変更する手段(102)を含むレーザ共振器と、波長制御信号(115)によって前記発振波長を変更する手段を制御して、前記発振波長を掃引する波長制御手段(113)と、前記波長制御手段から得られた掃引開始からの複数の時間と、前記複数の時間の各々において前記波長掃引光のコヒーレンス長が最長となる光利得とを組み合わせた光利得制御データに基づいて、光利得制御信号(116)によって前記光利得を変更する手段を制御する光利得制御手段(111)とを備えたものとして実施できる。
 好ましくは、前記発振波長を変更する手段は、回折格子(107)と、前記回折格子からの発振光が入射し、当該入射した前記発振光を垂直に反射するミラー(108)と、前記回折格子への前記発振光の入射角を変更する光偏向器(106)とから構成することができる。また、前記ゲイン媒体の光利得を変更する手段は、半導体光増幅器(SOA)であって、前記光利得制御信号は前記SOAを駆動するSOA電流であり、光利得制御データは前記複数の時間の各々に対して選択されたSOA電流を表すSOA電流波形データとして実施できる。
 このように、レーザ共振器101が波長を変更する手段および光利得を変更する手段を備えており、それぞれを制御する波長制御手段および光利得制御手段があれば、本発明は、光利得制御データに基づいて動作する波長掃引光源にも適用できる。以後、説明する本発明の実施例についても同様である。上述の実施例では、本発明の波長掃引光源100を動作させるためのSOA電流波形データが既に存在していることを前提に説明したが、次の実施例では、このSOA電流波形データを生成する装置および一連の手順について説明する。
 図3は、本発明の波長掃引光源に関連付けられたメモリに記憶するSOA電流波形データを生成する装置の構成を示す図である。図3のSOA電流波形データを生成する装置は、実施例1の波長掃引光源101の要素を全て含み、これに加えて、SOA電流波形データを取得するための要素を含んでいる。しかしながら、図3における追加の要素については、実施例3で後述するように、干渉計を除く多くの要素がCPU等によるソフトウェア演算処理によって実現可能な処理プロセスを示す機能ブロックを表している。したがって、図3における追加の要素は、干渉計を除けば、プロセッサ以外の特別なハードウェア要素なしにソフトウェア演算処理として実現できる点に留意されたい。具体的には、図3における制御部315は、CPUなどのプロセッサおよび関連付けられたメモリなどのハードウェア要素を含むものである。一方、電流-時間-コヒーレンス長対応データ取得手段313および最大コヒーレンス長対応電流選択手段314等は、図4のフローチャートとともに後述するSOA電流波形データを取得するための一連の手順の要素を便宜的に示したものと理解できる。したがって、電流-時間-コヒーレンス長対応データ取得手段313および最大コヒーレンス長対応電流選択手段314は便宜的に別個のものとして描かれているが、制御部315の制御の下で後述する一連のSOA電流波形データを取得する手順(図4)を実行するものと理解されたい。したがって、制御部315、電流-時間-コヒーレンス長対応データ取得手段313および最大コヒーレンス長対応電流選択手段314間の接続関係、包含関係や矢印等は、制御信号、制御データ、演算結果データなどのやり取りを大まかに概念的に示しているものである点に留意されたい。
 また、電流-時間-コヒーレンス長対応データ取得手段313の中に含まれている時間-コヒーレンス長対応データ取得手段312は、発振出力光のコヒーレンス長を取得する手順を特に1つの明示的な機能ブロックとして表したものである。本発明の波長掃引光源では、動的にSOA電流を変化させるときの各波長におけるSOA電流値を決定するために、発振出力光のコヒーレンス長を取得する必要がある。本発明の波長掃引光源においては、SOA電流波形データ取得の前提として、このコヒーレンス長を求める方法および構成にも特有の新規な側面を含んでおり、これが実施例3として開示される。したがって、時間-コヒーレンス長対応データ取得手段312の構成およびコヒーレンス長を求める動作・手順は、実施例3によって詳細に説明される。
 本実施例では、図3のSOA電流波形データを生成する装置において、SOA電流波形データを求める手順の概要を説明する。図3に示したレーザ共振器101、電流源111および光偏向器ドライバ113は、実施例1の波長掃引光源で使用されたものと同一であって、構成は前述した通りである。ただし、任意波形電流源111については、その使用方法が実施例1の波長掃引光源の場合とは異なる。
 実施例1では本発明の波長掃引光源を実際に光源として利用する際に、既に取得済みのSOA電流波形データに基づいて、波長掃引と同期して時間とともに動的に変化するSOA電流を生成するものとして動作していた。本実施例では、任意波形電流源111は、時間的に変化させずにSOA102に一定値の電流を流す電流源として動作する。簡単のため、以下の説明では電流源111とする。電流源111におけるSOA電流値は、外部からSOA電流設定値データによって設定できる。図3では、電流-時間-コヒーレンス長対応データ取得手段313から設定するものとしているが、これは一例であって、制御部315からSOA電流設定値データを設定しても良い。
 電流-時間-コヒーレンス長対応データ取得手段313は、電流源111にSOA電流値データを入力してSOA102に所定のSOA電流値を設定し、設定したSOA電流値毎に時間-コヒーレンス長対応データを取得する。すなわち、設定したSOA電流値毎に、「時間」および「コヒーレンス長」のデータ組(時間,コヒーレンス長)を求める。ここで言うところの「時間」とは、1回の波長掃引期間における掃引開示時刻からの時間(時刻)を意味している。設定した特定のSOA電流値および特定の時間の条件におけるコヒーレンス長の取得は、時間-コヒーレンス長対応データ取得手段312によって行われる。したがって、時間-コヒーレンス長対応データ取得手段312は、端的には特定条件の下でコヒーレンス長の測定を行う機能ブロックを表している。時間-コヒーレンス長対応データ取得手段312からは、時間-コヒーレンス長対応データ316、すなわち端的には、データ組(時間,コヒーレンス長)316が出力される。電流-時間-コヒーレンス長対応データ取得手段313は、設定したSOA電流値ごとにデータ組(時間,コヒーレンス長)316を取得し、最終的に、「SOA電流値」、「時間」、「コヒーレンス長」の3種類のデータを含むデータ組(SOA電流値,時間,コヒーレンス長)を求める。上述のデータ組の詳細は、後述する図4のフローチャートおよび図5の取得データ例とともに容易に理解できるだろう。
 上述のように、時間-コヒーレンス長対応データ取得手段312は、光偏向器ドライバ113からの同期信号318に同期してレーザ共振器101から波長掃引光を得て、レーザ共振器101からの波長掃引光に対して「時間」毎のコヒーレンス長を測定する。このとき、時間-コヒーレンス長対応データ取得手段312では、SOA102を1つのSOA電流値に設定した状態で、波長掃引光から異なる干渉条件でそれぞれ異なる干渉光を発生させて、1つのSOA設定電流値に対する複数の干渉光データから、「時間」毎にコヒーレンス長が計算される。また、本発明の波長掃引光源のためのSOA電流波形データは、波長掃引動作をしている「波長掃引光」から得られるものである点に注目されたい。すなわち、波長掃引光源を特定の発振波長に設定して、静的に動作させてコヒーレンス長を求めるのではなく、実際の波長掃引動作の条件の下で得られる波長掃引光に基づいて、波長掃引期間の各「時間」に対するコヒーレンス長が計算される。波長掃引期間の各「時間」は、図2の上側のグラフに示されるように、掃引波長帯域内の各波長に1対1に対応している。したがって、時間-コヒーレンス長対応データ取得手段312は、波長掃引光源において、発振波長が動的に刻々と変化する「実際の波長掃引動作状態」の波長掃引光から、波長毎にコヒーレンス長を求めることになる。
 最大コヒーレンス長対応電流選択手段314は、電流-時間-コヒーレンス長対応データ取得手段313によって求められた電流-時間-コヒーレンス長対応データ316を得て、波長掃引期間の各「時間」に対して、コヒーレンス長が最大となるSOA電流値を選択・決定する。言い換えると、「SOA電流値」、「時間」、「コヒーレンス長」の3種類のデータ組(SOA電流値,時間,コヒーレンス長)の内で、同一の「時間」を持つデータの中から、コヒーレンス長が最大であるSOA電流値を選択することになる。最後に、各「時間」において選択したSOA電流値のデータ(時間-最大コヒーレンス長時電流対応データ)、つまり、電流波形データ317を出力する。次に、SOA電流波形データを求める詳細な手順をフローチャートとともに説明する。
 図4は、本発明の波長掃引光源に関連付けられたメモリに記憶するSOA電流波形データを生成する手順のフローチャートを示す図である。以下の説明では、図4のフローチャートの各要素(箱)に付した符号と以下のステップの番号が対応している。
 SOA電流波形データを生成する手順400は、まずステップ401において、各種パラメータ設定する。制御部315から電流-時間-コヒーレンス長対応データ取得手段313に、電流源111に設定するSOA電流値データの初期値、最終値および電流変化ステップが設定される。電流源111に対して設定されるのは、電流ではなくて設定するSOA電流値に対応する設定データであるが、以下の説明では、簡単のためSOA電流値等と呼ぶ。
 手順400は、ステップ402において、SOA電流初期値が設定される。電流-時間-コヒーレンス長対応データ取得手段313は、電流源111に、SOA電流値の初期値を設定する。初期値から開始して、電流変化ステップでSOA電流値を順次変化させ、各SOA電流値に対して、後述するステップ403で、1回の波長掃引における時刻毎のコヒーレンス長を計測する。SOA電流値は、最初に小さい値として徐々に大きくしても良いし、逆でも良い。
 手順400は、ステップ403において、時間-コヒーレンス長対応データを取得する。コヒーレンス長を求める手順の詳細は、実施例3として後述するので、ここでは概要を簡単に述べる。時間-コヒーレンス長対応データ取得手段312は、光偏向器ドライバ113からの同期信号318に同期して、レーザ共振器101から1回の波長掃引分の波長掃引光を取得する。1回の波長掃引分の波長掃引光を取得した後で、実施例3で後述する干渉計の2つのアーム(サンプルアーム、リファレンスアーム)の光路長差を変えて干渉条件を変更し、異なる干渉条件で対応する複数の干渉信号を取得する。取得した複数の干渉光をそれぞれリスケーリングした後に、波長掃引中の時刻毎に光路長差毎のPSFを算出する。その後、各時刻に対して、各時刻に対して、PSFピーク値(強度)が最大値PSFMAXの1/2となるPSFピーク値を持つ2つの光路長差2zPSFHALF1および2zPSFHALF2の間の差分値を求める。この差分値の絶対値2|zPSFHALF1-zPSFHALF2|を、コヒーレンス長Lcとする。ここで、各光路長差2zPSFHALF1、2zPSFHALF2はお互い正負の符号が逆転している。光路長差の正負と測定の仕組みについて以下に述べる。
 リファレンスアームの光路長を基準としたときのサンプルアームの光路長との差を光路長差とすると、サンプルアームの光路長が短い場合には光路長差は負となり、その逆の場合には光路長差は正となる。それぞれの光路長差に対してPSFが計算されるが、光路長差が0に近いものほど干渉が強くなるためPSFピーク値は大きくなり、光路長差が長くなるに従って干渉が弱くなるためPSFピーク値は小さくなる。上述の測定では、光路長差が正負両方とも測定される。このため、PSFピーク値は光路長差が0付近で最大となり、光路長差0の状態から離れるに従って小さくなる。コヒーレンス長は、干渉強度が最大値から半減するまでの可干渉距離として定義されるため、上述のように、PSFピーク値がその最大値から半減する光路長差として求められる。
 もしも、光路長差0を中心に光路長差が正負どちらに対してもPSFピーク値の減衰の仕方が左右対称となる場合は、光路長差が正側または負側だけの測定でも良い。この場合は、光路長差2zが0に近いPSFピーク値(強度)をPSFの最大値PSFMAXとし、そのPSFMAX/2となるPSFピーク値を持つ時の光路長差2zPSFHALFを2倍して、この長さ4zPSFHALFを各時刻におけるコヒーレンス長としても良い。
 ここで光路長差2zが0に近いPSFピーク値をPSFの最大値PSFMAXとする理由は、光路長差2zが0のPSFが計算できないからである。PSFは干渉信号から計算され、実質的には、リスケーリング後の干渉信号の包絡線のフーリエ変換結果がPSFの形状となる。光路長差2zが0の時は、干渉信号の搬送波にあたる周波数は0となり、干渉信号の包絡線を安定に出力することは難しいためである。コヒーレンス長測定においては、光路長差2z=0のPSFを取得する代わりに、光路長差2z=0に近い状態でのPSFを取得することが通常も行われている。最後に、波長掃引中の各時刻およびコヒーレンス長の対応データ(時間-コヒーレンス長対応データ316)を出力する。
 手順400は、ステップ404において、電流-時間-コヒーレンス長対応データを取得する。電流-時間-コヒーレンス長対応データ取得手段313は、ステップ402で出力した時間-コヒーレンス長対応データ316と、電流源111に現在設定中のSOA電流値とを関連付ける。既に、電流-時間-コヒーレンス長対応データが存在していれば、その電流-時間-コヒーレンス長対応データに、このステップで新たにSOA電流と関連付けた時間-コヒーレンス長対応データを追記する。電流-時間-コヒーレンス長対応データが存在していなければ、本ステップで新たにSOA電流と関連付けた時間-コヒーレンス長対応データを、電流-時間-コヒーレンス長対応データとする。既に述べたように、電流-時間-コヒーレンス長対応データは、各々が「SOA電流値」、「時間」、「コヒーレンス長」の3種類のデータを含むデータ組(SOA電流値,時間,コヒーレンス長)である。
 手順400は、ステップ405において、SOA電流が最終値かどうか判断をする。電流-時間-コヒーレンス長対応データ取得手段313は、電流源111に現在設定中のSOA電流値が最終値か否かを判断する。Yesであればステップ407へ、Noであればステップ406へ進む。
 手順400は、ステップ406において、SOA電流値を変更する。電流-時間-コヒーレンス長対応データ取得手段313は、電流源111に設定するSOA電流値を、電流変化ステップに従って現在のSOA電流値から新たなSOA電流値に変更する。そして、ステップ403へ戻り、新たなSOA電流値において、時間-コヒーレンス長対応データ316を取得する。すなわち、データ組(時間,コヒーレンス長)316が得られる。
 手順400は、ステップ407において、時間-最大コヒーレンス長対応電流データ(電流波形データ)が取得される。最大コヒーレンス長対応電流選択手段314は、ステップ404で取得された電流-時間-コヒーレンス長対応データに対して、波長掃引の時刻毎に、コヒーレンス長が最大となるSOA電流値を選択する。そして、時刻毎にそのSOA電流値のデータ(時刻-最大コヒーレンス長時電流対応データ317)、つまり、SOA電流波形データを出力する。
 図5は、図4のフローチャートに示した各ステップを通じて取得されるコヒーレンス長データの全体を簡単化して示した表である。図5に示した表では、表左端の縦方向に、設定したSOA電流値を示している。すなわち、最小SOA電流値は30mAで、最大SOA電流値は100mAであり、電流変化ステップは10mAである。また図5に示した表の最上段に横方向に、コヒーレンス長を求めた波数kを示している。波数は、波長掃引の時刻または波長に1対1に対応しており、実施例3で詳細に説明するように、等時間間隔の中心波数k(0・δt)~k(7・δt)でコヒーレンス長を求めている。
 表5のSOA電流値(mA)および波数kの組み合わせに対して、表5の各セルにおいてコヒーレンス長が示されている。1つのSOA電流値に対して、中心波数がk(0・δt)からk(7・δt)までの各波数に対して、コヒーレンス長が求められる。例えば、SOA電流が100mAである表の最上段のデータでは、中心波数k(0・δt)~k(7・δt)について左から右に9.32~10.65の8個のコヒーレンス長データが求められている。これらの8個のコヒーレンス長データは、図4のフローチャートにおけるステップ404で得られる時間-コヒーレンス長対応データ316に対応する。後述する図8からもわかるように、波数kは波長掃引期間における時間(時刻)と1対1に対応している点に留意されたい。
 上述の8個のコヒーレンス長データを、異なるSOA電流値に対して順次求めることで、表5の全体のセルのコヒーレンス長が求められる。これは、図4のフローチャートにおけるステップ405でSOA電流値が最終値かどうかを判断し、ステップ406でSOA電流値を新たな値に変更して、ステップ404で異なるSOA電流値に対して順次コヒーレンス長の取得が繰り返し実施されることに対応する。
 さらに図4のフローチャートのステップ407に対応して、表5における1つの波数に対応する縦方向に並んだ8つのコヒーレンス長データに対して、コヒーレンス長が最大となるSOA電流値が1つ選択される。例えば、最も左側の列の中心波数がk(0・δt)であるとき、最大のコヒーレンス長11.11は、ハッチングをしたセルに相当しSOA電流が60mAの時に得られる。したがって、波数がk(0・δt)において、コヒーレンス長が最大となるSOA電流値として60mAが選択される。中心波数k(0・δt)~k(7・δt)の各中心波数において、それぞれコヒーレンス長が最大となるSOA電流値が選択される。すなわちハッチングされたセルのSOA電流値が選択される。尚、図5の表において、コヒーレンス長がないセルは、そのSOA電流値および中心波数の条件において、干渉光が安定しないためコヒーレンス長の計算が異常となってしまう場合、発振自体が消失する場合を意味する。
 したがって、本発明の波長掃引光源は、前記光利得制御手段(111)に関連付けられ、前記光利得制御データを記憶したメモリをさらに備え、前記光利得制御データは、前記光利得を変更する手段が、波長掃引期間に渡って、一定値の前記光利得制御信号によって前記光利得を変更する手段を設定光利得値に制御し、複数の異なる値の前記設定光利得値の各々に対して、掃引開始からの前記複数の時間の各々におけるコヒーレンス長をそれぞれ測定し、前記設定光利得値、前記複数の時間の内の1つの時間および対応する前記測定されたコヒーレンス長からなるデータ組を生成し、前記複数の時間の内の前記1つの時間において、前記測定されたコヒーレンス長の中から最大コヒーレンス長を与える設定光利得値を選択し、前記1つの時間の各々と前記選択された設定光利得値から生成されるものとしても実施できる。
 上述の図3~図5の説明ではSOA電流波形データを求める手順を説明したが、これも実施例1と同様により一般化が可能である。図3の光偏向器106を「波長を変更する手段」に、光偏向器ドライバ113を波長制御手段に、偏向器ドライブ信号を波長制御信号に、SOA102を「光利得を変更する手段」に、SOA電流を光利得制御信号または光利得にそれぞれ置き替えても、図3~図5で説明したSOA電流波形データの取得手順の考え方をそのまま適用できる。
 図3の電流源111は、SOA電流値を得て時間的に一定値のSOA電流を流すよう動作するが、SOA102を光利得を変更する手段としたときに、時間的に光利得が一定値となるように光利得変更手段を制御する手段(一定光利得制御手段)を備えれば良い。
 また図3の電流-時間-コヒーレンス長対応データ取得手段313は、電流源111にSOA102のSOA電流を設定する印加電流値を渡すよう動作するが、電流源111の代わりに一定光利得制御手段が光利得を変更する手段を制御する場合は、一定光利得制御手段に光利得値を入力して光利得変更手段の光利得を変え、時間-コヒーレンス長対応データ取得手段312から光利得値毎の時間-コヒーレンス長対応データを得る。そして、光利得毎の時間-コヒーレンス長対応データ、すなわち、光利得-時間-コヒーレンス長対応データを出力する。この場合は、電流-時間-コヒーレンス長対応データ取得手段313を、光利得-時間-コヒーレンス長対応データ取得手段と呼ぶことができる。
 最大コヒーレンス長対応電流選択手段314も同様であり、光利得-時間-コヒーレンス長対応データ取得手段313から得た光利得-時間-コヒーレンス長対応データを得て、1回の波長掃引中の時刻毎に、コヒーレンス長が最大となる光利得を選択する。そして、時刻毎の選択した光利得のデータ(時間-最大コヒーレンス長利得対応データ)、すなわち、光利得制御データを出力する。
 最大コヒーレンス長対応電流選択手段314によって取得される電流-時間-コヒーレンス長対応データのデータ取得の時間間隔δtが、電流源111が制御する電流の時間間隔δt’と異なる場合、これらの時間間隔の間で時間的な整合が取れるようにする。例えば、最大コヒーレンス長対応電流選択手段314内で生成した時間-最大コヒーレンス長電流対応データの時間間隔はδtとなるが、線形的または非線形的な補間によって時間間隔をδt’に作り直すことができる。
 例えば、δt>δt’とすると、電流-コヒーレンス長対応データを時間的に粗く、低頻度で取得することになる。結果として、電流-時間-コヒーレンス長対応データも時間的に粗く取得されることになり、処理時間を短縮できる効果がある。時間-最大コヒーレンス長時電流対応データの波形の最高周波数をfとするとき、δt≦1/(2f)とすれば、全周波数成分を折り返し歪等無く正常に取得することができる。
 次に、図3の時間-コヒーレンス長対応データ取得手段312において、コヒーレンス長Lcを求める手順について、実施例3として詳細に説明する。したがって、実施例1に示した本発明の波長掃引光源は、実施例2に示した手順によって得たSOA電流波形データを利用する。さらに実施例2で得られたSOA電流波形データは、次に説明する実施例3のコヒーレンス長を求める装置および手順を利用して求められるという関係になる。
 最初に述べたように、図3に示した構成において、例えばこれに限定されないがCPUなどを含む制御部315は、図1に示した波長掃引光源がそもそも備えているものである。また、光利得-時間-コヒーレンス長対応データ取得手段313および最大コヒーレンス長対応電流選択手段314のいずれについても、CPUなどを使ったソフトウェアによる演算処理によって実現できる。したがって、本発明の波長掃引光源は、図3の構成要素をすべて備えたものとしても実施できる。実施例2の図4のフローチャートの手順は、波長掃引光を用いて行われるので、波長掃引光源を実際の光源としての用途に使用している最中であっても、SOA電流波形データを取得したり、SOA電流波形データ更新したりすることも可能である。したがって、例えば波長掃引光源を使用する際の温度などの環境条件の変動に応じて、波長掃引光源が動作している環境中にあるときにSOA電流波形データ更新することで、常に波長掃引光源を最適な状態で動作させることもできる。
 図6は、本発明の波長掃引光源において使用するSOA電流波形データを求めるためのコヒーレンス長を測定する装置の構成を示す図である。図6では、図3に示した時間-コヒーレンス長対応データ取得手段312の詳細な構成が点線の領域内に示されている。図6におけるレーザ共振器101、電流源111および光偏向器ドライバ113は、実施例2の図3の場合と同一である。コヒーレンス長を求める手順の概要は、次の通りである。
 時間-コヒーレンス長対応データ取得手段312は、光偏向器ドライバ113からの同期信号318に同期して、レーザ共振器101から1回の波長掃引分の波長掃引光123を取得する。1回の波長掃引分の波長掃引光を取得した後で、干渉計501の2つのアーム(サンプルアーム、リファレンスアーム)光路長差を変えて干渉条件を変更し、異なる干渉条件で対応する複数の干渉信号を取得する。取得した複数の干渉光521をそれぞれリスケーリングした後に、時間-光路長差-PSF対応データ取得手段515が波長掃引中の時刻毎に光路長差毎のPSFを算出する。その後、各時刻に対して、PSFピーク値(強度)が最大値PSFMAXの1/2となるPSFピーク値を持つ2つの光路長差2zPSFHALF1=2qPSFHALF1δzおよび2zPSFHALF2=2qPSFHALF2δzの間の差分値を求め、その差分値の絶対値2|zPSFHALF1-zPSFHALF2|=2δz|qPSFHALF1-qPSFHALF2|をコヒーレンス長Lcとする。ここで、各光路長差2zPSFHALF1、2zPSFHALF2はお互い正負の符号が逆転している。光路長差の正負と測定の仕組みについて以下に述べる。
 リファレンスアームの光路長を基準としたときのサンプルアームの光路長との差を光路長差とすると、サンプルアームの光路長が短い場合には光路長差は負となり、その逆の場合には光路長差は正となる。それぞれの光路長差に対してPSFが計算されるが、光路長差が0に近いものほど干渉が強くなるためPSFピーク値は大きくなり、光路長差が長くなるに従って干渉が弱くなるためPSFピーク値は小さくなる。上述の測定では、光路長差が正負両方とも測定される。このため、PSFピーク値は光路長差が0付近で最大となり、そこから離れるに従い小さくなる。コヒーレンス長は、干渉強度が最大値から半減するまでの可干渉距離として定義されるため、上記の通り、PSFピーク値がその最大値から半減する光路長差として求められる。
 もしも、光路長差2z=0を中心として、光路長差が正負どちらに対してもPSFピーク値の減衰の仕方が左右対称となる場合は、光路長差が正側または負側だけの測定でも良い。この場合は、光路長差2zが0に近いPSFピーク値(強度)をPSFの最大値PSFMAXとし、PSFMAX/2となるPSFピーク値を持つ時の光路長差2zPSFHALFを2倍して、この長さ4zPSFHALFを各時刻におけるコヒーレンス長Lcとしても良い。
 ここで光路長差2zが0に近いPSFピーク値をPSFの最大値PSFMAXとする理由は、光路長差2zが0のときのPSFを計算できないからである。PSFは干渉信号から計算され、実質的には、リスケーリング後の干渉信号の包絡線のフーリエ変換結果がPSFの形状となる。光路長差2zが0の時は、干渉信号の搬送波にあたる周波数は0となり、干渉信号の包絡線を安定に出力することは難しいためである。コヒーレンス長測定においては、光路長差2z=0のPSFを取得する代わりに、0に近い光路長差でのPSFを取得することが通常でも行われている。光路長差2zが0の位置のPSFをPSFの最大値PSFpeak_MAXの代わりに、光路長差の異なる複数のPSFピーク点(2zq、PSFpeak_q)(ただし、qは複数のPSFピーク点を区別する整数)を使って、関数f(2z)フィッティングし、光路長差2z=0となる関数の値f(0)をPSFの最大値PSFpeak_MAXとすることも可能である。この関数f(x)は、たとえば、ガウシアン、ローレンチアン、xを引数とした関数の多項式等が好適である。
 最後に、波長掃引中の各時刻およびコヒーレンス長Lcの対応データ(時間-コヒーレンス長対応データ316)を出力する。
 図6における時間-コヒーレンス長対応データ取得手段312内の要素は、コヒーレンス長Lcを求めるための手順のステップ要素を、機能ブロックとして表したものである。すなわち、干渉計501以外の要素は、CPU等によるソフトウェア演算処理によって実現可能な処理プロセスを示す機能ブロックを表している。したがって、干渉計501を除けば、プロセッサ以外の特別なハードウェア要素なしにソフトウェア演算処理として実現できる点に留意されたい。尚、干渉計501の具体的な構成例は、最後に図9A~図9Dとともに説明する。以下、図6の各処理ブロックの機能および動作について、詳細に説明する。
 レーザ共振器101からの波長掃引光123は、干渉計501に供給される。干渉計501は、レーザ共振器101からの波長掃引光123をファイバコリメータ(FC)502で空間に出射する。出射光は、ビームスプリッタ(BS)503で、参照アーム側およびサンプルアーム側の2つの経路に分けられる。参照アームおよびサンプルアームには、それぞれミラー504およびミラー505が設置されている。それぞれのアームに入射した光は、反射してBS503に戻され、BS503で合波され干渉する。その干渉光はフォトディテクタ(PD)507に入射され、電気信号(干渉信号)521としてPD507から干渉計501の外部に出力される。
 サンプルアーム側のミラー505はミラー位置可動手段506上に設置されており、BS503とサンプルアーム側ミラー505との間の距離を外部から制御できるようになっている。すなわち、参照アーム長およびサンプルアーム長の間の差zを、干渉計501の外部から制御することができる。
 干渉信号取得手段511は、後述する時間-光路長差-PSF対応データ取得手段515からミラー移動完了信号532を受け取った後、光偏向器ドライバ113から得た同期信号に同期して、1回の波長掃引分の干渉信号123を取得し出力する。ミラー移動完了信号532とは、干渉計501内にあるミラー位置可動手段506が、光路長差変更のためにサンプルアーム側ミラーの1回の移動を終えたことを示す信号である。
 リスケ―リング手段512は、時間-波数対応データを使って、干渉信号取得手段から得た干渉信号522をリスケーリングし、そのリスケーリングされた干渉信号523、すなわちi(k)(kは波数)を出力する。リスケーリングについては、非特許文献1に詳説されているが、干渉信号522を再度サンプリング(リサンプリング)して時間軸を伸縮するように動作する。これによって、干渉信号522は、時間に対して周波数が一定の信号に変換される。
 処理時刻発生手段518は、1回の波長掃引に掛かる時間Tと時間ステップδtから、以後のプロセスで処理を実行する時刻情報531を発生する。具体的には、0、δt、2δt、…、(n-1)δtを順番に発生する。ただし、n=T/δtである。
 抽出位置算出手段517は、時間-波数対応データから、処理時刻発生手段518から得た時刻531すなわちpδt(pは0~n-1)に対応する波数、すなわち抽出位置536を求め、これを出力する。具体的には、抽出位置(波数)k(pδt)を出力する。
 処理範囲抽出手段513は、抽出幅Δkおよび抽出位置算出手段518から得た抽出位置k(pδt)を使って、リスケーリング手段512から得たリスケーング済み干渉信号523、すなわちi(k)を、k(pδt)を中心に±Δk/2となる範囲で抽出し、出力する。このとき、たとえば、中心k(pδt)、半値幅±Δk/2となる窓関数(たとえば、方形窓、ガウシアン窓、ハニング窓、ハミング窓、ブラックマン窓等)を使って抽出しても良い。
 窓関数・周波数変換手段514は、処理範囲抽出手段513から得られた干渉信号524に窓関数を掛け、フーリエ変換を使ってパワースペクトル(または、位相を除いた周波数分布、すなわちパワースペクトルの平方根)を計算し、プラス側だけを出力する。この出力された信号525は、点広がり関数(PSF:point spread function)と呼ばれる。
 窓関数については、処理範囲抽出手段513として、ガウシアン窓、ハニング窓、ハミング窓、ブラックマン窓等の中心から離れるに従い強度が小さくなる窓を使っている場合は、使用しなくても良い。また、干渉信号の不連続点に応じて、窓の有無や仕様・種類を決定しても良い。
 時間-光路長差-PSF対応データ取得手段515は、窓関数・周波数変換手段514からPSFデータを得て、処理時刻発生手段518から出力された時刻pδt(pは0~n-1、ただし、n=T/δt)と、干渉計501のアーム長差zの2倍(2z=2qδz;光路長差)と関連づけてPSFデータを保存する。FC502からPD507までの、参照アームおよびサンプルアームの間の光路長差は2zである。
 その後、干渉計501に組み込まれているミラー位置可動手段506にミラー位置信号533を送ってサンプルアーム側ミラー505の位置を変更し、位置移動が完了した後に、ミラー移動完了信号532を干渉信号取得手段511へ出力する。このミラー位置信号533によって変更するサンプルアーム側ミラーの位置は、時間-光路長差-PSF対応データ取得手段515に予め設定したミラー位置範囲±Δz(z=0からの最大距離をΔzとする)および位置ステップδzを使って決定する。例えば、z=qδzとする。ここで、m=Δz/δzであり、q=-m~mの整数である。
 時間-光路長差-PSF対応データ取得手段515がサンプルアーム側ミラー505の位置移動が完了したがどうかを検知する方法としては、ミラー位置可動手段506から移動完了信号532を得て、移動完了したと判断して良い。また、常時ミラー位置可動手段506のサンプルアーム側ミラーの位置を観測して、ミラー位置信号と一致しているか、または、あらかじめ設定した誤差内に収まっている場合に、移動完了したと判断しても良い。さらに、ミラーの移動に十分な時間が経った後に、ミラーが移動完了したと判断しても良い。
 時間-光路長差-PSF対応データ取得手段515は、処理時刻発生手段518から得た時刻pδtにおけるpがn-1に達した時は、時間-光路長差-PSF対応データ534を出力する。出力された時間-光路長差-PSF対応データ534には、時刻が0~(n-1)δtの範囲で、光路長差が-mδz~mδzの範囲のPSFデータが蓄積されている。
 時間-コヒーレンス長対応データ生成手段516は、時間-光路長差-PSF対応データ取得手段515から得た時間-光路長差-PSF対応データ534から、各時刻pδt(pは0~n-1)のコヒーレンス長Lcを算出し、pδtおよびLcを関連付けて保存する。時刻0~(n-1)δtまでの1回の波長掃引分の全てのデータを取得して1つの時間-コヒーレンス長対応データとして保存したら、その時間-コヒーレンス長対応データ316を出力する。
 コヒーレンス長Lcを得る方法の一例としては、時刻pδt毎に各光路長差に対応するPSFのピーク値を求め、まず、PSFピーク値の最大値PSFMAXの1/2となるPSFピーク値を持つ2つの光路長差2zPSFHALF1=2qPSFHALF1δzおよび2zPSFHALF2=2qPSFHALF2δzの差分値を求める。そして、この差分値の絶対値2|zPSFHALF1-zPSFHALF2|=2δz|qPSFHALF1-qPSFHALF2|をコヒーレンス長Lcとする。2z=2qδzでちょうどPSFのピーク値がPSFMAX/2とならない場合には、PSFMAX/2となる2zの位置の両側の位置から数個のPSFのピーク位置およびピーク値を使って線形的または非線形的な補間をし、PSFMAX/2となる2zの位置を算出することもできる。ここで、各光路長差2zPSFHALF1、2zPSFHALF2はお互い正負の符号が逆転している。光路長差の正負とコヒーレンス長Lc測定の仕組みについて以下さらに述べる。
 リファレンスアームの光路長を基準としたときのサンプルアームの光路長との差を光路長差とすると、サンプルアームの光路長が短い場合には光路長差は負となり、その逆の場合には光路長差は正となる。それぞれの光路長差に対してPSFが計算されるが、光路長差が0に近いものほど干渉が強くなるためPSFピーク値は大きくなり、光路長差が長くなるに従って干渉が弱くなるためPSFピーク値は小さくなる。上述の測定では、光路長差が正負両方とも測定される。このため、PSFピーク値は光路長差が0付近で最大となり、光路長差0の状態から離れるに従って小さくなる。コヒーレンス長は、干渉強度が最大値から半減するまでの可干渉距離として定義されるため、上記の通り、PSFピーク値がその最大値から半減する光路長差として求められる。
 もしも、光路長差2z=0を中心として、光路長差が正負どちらに対しても、PSFピーク値の減衰の仕方が左右対称となる場合は、光路長差が正側または負側だけの測定で良い。この場合は、光路長差2zが0に近いPSFのピーク値(強度)をPSFの最大値PSFMAXとし、PSFMAX/2となるPSFピーク値を持つ時の光路長差2zPSFHALFを2倍して、この長さ4zPSFHALFをコヒーレンス長Lcとしても良い。2z=2qδzでちょうどPSFのピーク値がPSFMAX/2とならない場合には、PSFMAX/2となる2zの位置の両側の位置から数個のPSFのピーク位置およびピーク値を使って線形的または非線形的な補間をし、PSFMAX/2となる2zの位置を算出することもできる。このとき、コヒーレンス長Lcはその2zのさらに2倍の4zとなる。
 コヒーレンス長Lcを測定する別の方法としては、光路長差の異なる複数のPSFピーク位置およびピーク値をガウシアンフィッティングし、ガウシアンの半値全幅をコヒーレンス長Lcとすることもできる。ガウシアンの他には、ローレンチアン、関数の多項式等が好適である。
 もし、光路長差2z=0を中心に光路長差が正負どちらに対してもPSFピーク値の減衰の仕方が左右対称となる場合は、光路長差が正側または負側だけの測定で良い。この場合は、光路長差2zが0に近いPSFピーク値(強度)をPSFの最大値PSFMAXとし、PSFMAX/2となるPSFピーク値を持つ時の光路長差2zPSFHALFを2倍して、この長さ4zPSFHALFを各時刻におけるコヒーレンス長としても良い。ここで光路長差2zが0に近いPSFピーク値をPSFの最大値PSFpeak_MAXとする理由は、光路長差2zが0のPSFが計算できないからである。
 PSFは干渉信号から計算され、実質的には、リスケーリング後の干渉信号の包絡線のフーリエ変換結果がPSFの形状となるが、光路長差2zが0の時は、干渉信号の搬送波にあたる周波数は0となり、干渉信号の包絡線を安定に出力することは難しいため、コヒーレンス長測定においては、光路長差2z=0のPSFを取得する代わりに、0に近い光路長差でのPSFを取得することが通常行われている。光路長差2zが0の位置のPSFをPSFの最大値PSFpeak_MAXの代わりに、光路長差の異なる複数のPSFピーク点(2zq、PSFpeak_q)(ただし、qは複数のPSFピーク点を区別する整数)を使って、関数f(2z)フィッティングし、光路長差2z=0となる関数の値f(0)をPSFの最大値PSFpeak_MAXとすることも可能である。この関数f(x)は、たとえば、ガウシアン、ローレンチアン、xを引数とした関数の多項式等が好例である。
 次に、時間-コヒーレンス長対応データ取得手段312によってコヒーレンス長Lcを求める手順を、図7のフローチャートとともに説明する。
 図7は、本発明の波長掃引光源において使用するSOA電流波形データを求めるためのコヒーレンス長を測定する手順を示すフローチャートである。以下の説明では、図7のフローチャートの各要素(箱)に付した符号と以下のステップの番号が対応している。
 コヒーレンス長Lcを求める手順600は、ステップ601において、各種パラメータ設定することで始まる。ステップ601で、下記のパラメータを設定する。
・時間-光路長差-PSF対応データ取得手段515に、ミラー位置の範囲±Δzとミラー位置のステップδzを設定する。
・処理時刻発生手段518に、1掃引時間Tおよび時間ステップδtを設定する。
・リスケ―リング手段512および抽出位置算出手段517に、時間-波数対応データk(t)を設定する。
 手順600は、ステップ602において、初期光路長差を設定する。時間-光路長差-PSF対応データ取得手段515は、ミラー可動手段506に初期のミラー位置信号533を出力する。ミラー位置可動手段506は、そのミラー位置信号533に従った位置(初期位置)にサンプルアーム側ミラー505を移動させる。
 ここで、ミラー位置信号533とは、ミラーの位置zを代表する何らかの信号である。
ミラー位置信号は何らかの電気信号を利用できるが、zを示すものとして、電圧値で示すこともできるし、RS-232CやUSBインタフェースにおいて、zの位置情報をコマンドで送ることもできる。つまり、z=qδzに相当する何らかの電気信号をミラー位置信号と呼ぶ。このように、ミラー位置信号533はz=qδzを表すものである。ここでm=Δz/δzであり、qは-m~mの整数である。昇順の場合はサンプルアームを-mδzの位置に移動させ、降順の場合はmδzの位置に移動させる。時間-光路長差-PSF対応データ取得手段515は、ミラーを移動させた後に、ミラー移動完了信号532を、干渉信号取得手段511へ出力する。
 手順600は、ステップ603において、干渉信号の取得を実行する。干渉信号取得手段511は、時間-光路長差-PSF対応データ取得手段515からミラー移動完了信号532を受け取ると、光偏向器ドライバ113から得た同期信号318に同期して、干渉計501から1回の波長掃引分の干渉信号521を取得し、出力する。
 手順600は、ステップ604において、リスケーリングを実行する。リスケ―リング手段512は、時間-波数対応データを使って、干渉信号取得手段511から得られた1回の波長掃引分の干渉信号521をリスケーリングする。
 手順600は、ステップ605において、初期処理時刻を発生する。処理時刻発生手段518は、1回の波長掃引に掛かる時間Tおよび時間ステップδtから、初期時刻531を発生し出力する。具体的には、例えば時刻情報として、0、δt、2δt、…、(n-1)δt(ただしn=T/δt)のように昇順で時刻情報を発生した場合、初期時刻として0を出力する。
 手順600は、ステップ606において抽出位置を算出する。抽出位置算出手段517は、時間-波数対応データに基づいて、処理時刻発生手段518から得た時刻pδt(pは0~n-1)に対応する波数536、すなわちk(pδt)を求め、これを出力する。ここで抽出位置とは、以後のステップで干渉信号を抽出する時間を意味している点に留意されたい。波長掃引光源の波長掃引期間においては、波長掃引開始からの時間(時刻)は、発振出力光の波長または波数に1対1に対応していることを思い出されたい。
 手順600は、ステップ607において、処理範囲を抽出する。処理範囲抽出手段513は、抽出幅Δk、並びに、抽出位置算出手段517から得た抽出位置536すなわちk(pδt)を使って、リスケーリング手段512から得たリスケーング済み干渉信号523すなわちi(k)を、k(pδt)を中心とした±Δk/2の範囲で抽出し、抽出された干渉信号524を出力する。このとき、例えば中心k(pδt)、半値幅±Δk/2となる窓関数、例えば方形窓、ガウシアン窓、ハニング窓、ハミング窓、ブラックマン窓等を使っても良い。
 手順600は、ステップ608において、抽出された干渉光に窓関数を適用し周波数変換を実行する。窓関数・周波数変換手段514は、処理範囲抽出手段513から得られた干渉信号524に窓関数を掛け、フーリエ変換を使ってパワースペクトル(または、位相を除いた周波数分布、つまり、パワースペクトルの平方根)を計算し、プラス側の周波数の分だけを出力する。この出力された信号525は、点広がり関数(PSF)である。本ステップにおける窓関数の使用については必須ではなく、使用しなくても良い。特に、ステップ607において方形窓以外のものを使用した場合は、本ステップにおける窓関数を使用する必要はない。
 手順600は、ステップ609において、時間-光路長差-PSF対応データを取得し、蓄積する。時間-光路長差-PSF対応データ取得手段515は、窓関数・周波数変換手段514から得たPSFデータ534を、時刻pδt、および、光路長差2z=2qδzと関連付けて、記憶手段に蓄積する。記憶手段は、実施例1で説明したように、波長掃引光源に関連付けられたメモリを使用することができる。
 手順600は、ステップ610において、時間更新を実行する。処理時刻発生手段518は、時刻pδtにおける変数pを更新する。具体的はp+1を更新したpとする。すなわち、(p+1)δtを更新した時刻とする。
 手順600は、ステップ611において、終了時間の判別を実行する。時間-光路長差-PSF対応データ取得手段515は、処理時刻発生手段518が更新した時刻pδtが終了時刻であるか否かを判別する。もしp=n-1でなければ、ステップ606に戻る。もし、p=n-1であれば、ステップ612へ進む。
 手順600は、ステップ612において、光路長差の更新を実施する。時間-光路長差-PSF対応データ取得手段515は、光路長差2zを更新するため、z=qδzの変数qを更新する。昇順の場合はqに1を加えた値を新たなqとなる更新をし、降順であればqから1を引いた値を新たなqとなる更新をする。
 手順600は、ステップ613において、光路長差更新が終了したかの判別を実行する。時間-光路長差-PSF対応データ取得手段515は、z=qδzが最終値となったか否かを判別する。もし、z=qδzが最終値でなければ、z=qδzとするミラー位置信号533をミラー位置可動手段506へ送る。このときミラー位置可動手段506がサンプルアーム側ミラー505をミラー位置信号533に応じた位置に移動させる。ミラーを移動した後、時間-光路長差-PSF対応データ取得手段515はミラー移動完了信号532を出力し、ステップ603へ戻る。もし、z=qδzが最終値であれば、蓄積した時間-光路長差-PSF対応データを出力する。なお、z=qδzの最終値は、zが昇順で更新される場合はmδzであり、降順で更新される場合は-mδzである。ここで、m=Δz/δzである。
 手順600は、ステップ614において、時間-コヒーレンス長対応データを生成する。時間-コヒーレンス長対応データ生成手段516は、時間-光路長差-PSF対応データ取得手段515から得た時間-光路長差-PSF対応データから、各時刻pδt(pは0~n-1)のコヒーレンス長Lcを算出し、pδtおよびLcを関連付けて保存する。時刻0~(n-1)δtまでの1回の波長掃引分の全てのデータを取得して1つの時間-コヒーレンス長対応データ316として保存したら、その時間-コヒーレンス長対応データ316を出力する。ステップ614を完了した段階で、1つのSOA電流に対して、各波数に対するコヒーレンス長が求まることになる。すなわち、図5を再び参照すれば、1つのSOA電流値に対して、中心波数がk(0・δt)からk(7・δt)までの各波数に対して、コヒーレンス長が求められる。
 したがって本発明の波長掃引光源においては、光利得制御データは、前記設定光利得値に設定されたレーザ共振器から出力された、前記波長掃引期間の1回の波長掃引から得られる波長掃引光(123)を、異なる干渉条件で干渉させる干渉計から生成された、前記干渉条件の各々に対応する複数の干渉信号(521)に基づいて、前記複数の干渉信号から求められた前記コヒーレンス長(316)を含むものとしても実施できる。
 図8は、本発明の波長掃引光源において、等時間間隔で求めたSOA電流の動作を説明する図である。図8は、上述の実施例3で求められたコヒーレンス長Lcに基づき、実施例2のSOA電流波形データを求める手順にしたがって、求めたSOA電流波形を示している。既に図2でも同様の説明をしているが、図2では、等波数間隔でコヒーレンス長を計算して、コヒーレンス長が最大となるSOA電流を決定したものである。図2の場合は、時間軸で見るとき、決定されたSOA電流の間隔は不均一である。また、等波長間隔でもコヒーレンス長を計算しても良いが、この場合も、決定されたSOA電流の時間間隔は不均一となってしまう。
 これに対し、図8のSOA電流波形は、実施例3で説明したように等時間間隔となるように、すなわちコヒーレンス長を計算する時刻の間隔(δt)が均一となるように計算時刻を決定している。したがって、0パディングしたデータをフーリエ変換して、SOA電流をより簡便かつ高速に補間することが可能となる。 
 図2および図8では、コヒーレント長を計算する波数kの数は8点だけしか描いてないが、実際には、波長掃引速度に応じて、最大で5000~10000点に対してコヒーレンス長を計算することが好適である。また図5の表では、SOA電流を10mAの電流可変ステップで取っているが、より細かいステップで取ることができるのは言うまでもない。
 図9A~図9Dは、図6に示した時間-コヒーレンス長対応データ取得手段内の干渉計の具体的な構成例を示した図である。図6では、マイケルソン型干渉計501の構成例を示したが、図9A~図9Dのようにマッハツェンダ型の各種構成によっても実現できる。図9Aはマッハツェンダ型(差分出力型、ファイバ干渉計型)を、図9Bはマッハツェンダ型(差分出力型)を、図9Cはマッハツェンダ型(差分出力型、ファイバ干渉計型)を、図9Dはマッハツェンダ型の構成をそれぞれ示す。図9A~図9Dにおいて各略号は、それぞれ、FC:ファイバコリメータ、CPL:カプラ、CCL:サーキュレータ、BPD:バランスドディテクタ、Ms:サンプルアームミラー、M:参照アームミラーを意味する。
 本発明の波長掃引光源におけるコヒーレンス長Lcの改善例を示せば、ガルバノスキャナを使用し、SOA電流が0~600mAの範囲で可変できるレーザ共振器においては、中心波長1310nm±40nmの波長掃引範囲において、20kHzで繰り返し波長掃引し、波長掃引期間が25μsの場合、全掃引波長帯域に渡ってコヒーレンス長Lcを1.7倍にすることができた。
 以上、詳細に述べたように、本発明の波長掃引光源およびその駆動データの作成方法によって、波長掃引期間中における波長掃引出力光のコヒーレンス長を常に長い状態に維持できる。これにより、SS-OCT装置におけるPSF強度を上げて、SS-OCT装置のイメージ検出感度を高くすることができる。
 上述の第1の実施形態の波長掃引光源では、波長掃引の全波長帯域にわたって、レーザ共振器内の光利得が、レーザ発振閾値の状態の光利得よりもわずかに高い光利得値となるように、光利得を変更する手段(SOA)を制御している。しかしながら、レーザ共振器の構成やタイプによっては、レーザ発振閾値の近傍と比べて大幅に高い光利得値において、コヒーレンス長が最大となるような場合もあり得る。このような場合であっても、本実施形態の波長掃引光源のように最大のコヒーレンス長を与えるSOA電流波形を求めることで、掃引波長帯域の全域に渡って、出力光のコヒーレンス長を長い状態に維持できる。
[第2の実施形態]
 上述の第1の実施形態では、波長掃引光源に使用されるSOAのSOA電流波形データを事前に求めて、その取得データに基づいて波長掃引光源を動作させる。したがって、波長掃引光源を通常動作させる前に、このSOA電流波形データを取得する工程が必要である。また、実施例3で説明したSOA電流波形データの取得のためにコヒーレンス長を求める手順は、図6および図9A~図9Dに示した処理ブロック図の干渉計501が必要であり、波長掃引光源の装置としては本来不要な構成要素が必要となる。
 本実施形態では、各掃引波長におけるコヒーレンス長を求めないで、光源の出力レベルに基づいて、より簡素な構成要素と処理によって、第1の実施形態と同様に、波長掃引期間中の波長掃引出力光のコヒーレンス長を常に長い状態に維持する例を示す。
 従来技術について既に述べたように、SS-OCTシステムではコヒーレンス長の1/4よりも深い層を観測しようとした場合に像が不明瞭になる。このため、コヒーレンス長の長い波長掃引光源が求められる。別の言い方をすれば、波長掃引光源では瞬時のレーザ線幅を狭く維持する必要がある。瞬時線幅の逆数はコヒーレンス長の関係にある。
 ここで、SOAの動作について再度検討する。SOAは、広い帯域幅と短いキャリア寿命に特徴があり、光励起を必要とせずに、電流注入により直接電気信号で駆動できる実用性にも優れる。一方で波長掃引光源に限らず一般的な連続発振の半導体レーザにおいても、注入電流の増加によりレーザ線幅が増大することが知られている。
 図10は、波長掃引光源においてSOAへの注入電流と光出力パワーおよびコヒーレンス長の関係を示した図である。横軸にSOAへの注入電流を示し、縦軸にある特定の波長におけるコヒーレンス長および出力光レベルを示している。図10から明らかなように、SOAへの注入電流の増大にともなって出力光レベルが増大するのに対し、コヒーレンス長は逆に低下する。具体的には図10では、特定の波長においてSOA注入電流が25mA以下の発振閾値に近い状態で、コヒーレンス長が急速に伸びている。つまり波長掃引光源においては、出力光レベルが増大する発振閾値付近でSOAを駆動したほうが、より長いコヒーレンス長が得られることを示している。しかしながら、SOAの利得スペクトルに波長依存性があるため、通常の一定電流によるSOA駆動では、このような発振閾値付近でコヒーレンス長が伸張した状態の駆動条件を全波長帯域に渡って実現できないことが問題であった。第1の実施形態の波長掃引光源では、コヒーレンス長を最大化する時間-SOA電流波形データを事前取得することによってこの問題に対処していた。
 本実施形態の波長掃引光源では、波長掃引光源の通常動作中に、出力光レベルをリアルタイムでモニタすることにより、さらに簡単な制御方法によって、全掃引波長帯域に渡って長いコヒーレンス長を維持するように動作する。波長掃引光源内の発振波長を変更する手段(波長選択素子)の動作に連動して、出力光のレベルに基づいてSOAへの注入電流を制御する。これによって、掃引波長帯域内の各波長において発振閾値付近でSOAを駆動することで、光源の瞬時発振線幅を狭め、コヒーレンス長を最大化する。本実施形態の波長掃引光源は、後述する波長掃引と連動した制御ステップを数回以上繰り返せば、少なくとも一定程度はコヒーレンス長が伸張した状態に直ちに近づき、光源の通常動作である波長掃引動作中においてコヒーレンス長を測定する必要はない。後述する光出力パワーの目標値の決定は、波長掃引光源装置の製造調整段階だけで実施可能であって、第1の実施形態の波長掃引光源においてコヒーレンス長を求めるのに必要だった干渉計も不要である。
 図11は、本発明の第2の実施形態の波長掃引光源の構成を示す図である。波長掃引光源700は、構成要素として、光を増幅する光増幅素子(SOA)701、特定の波長帯域の光を選択的にフィルタリングする発振波長を変更する手段(波長選択素子)702、光の一部を分岐し光源に帰還させる例えばハーフミラーなどの光帰還素子705、発振波長を変更する手段702を駆動するドライバ装置709、光増幅素子(SOA)701を駆動するドライバ装置708が含まれる。発振波長を変更する手段(波長選択素子)702に含まれる複数の光学素子は、後述するように光学的に往復光路を形成し、光がこれら素子間を繰返し往復することが可能な光共振器として構成されている。
 光増幅素子(以下、SOA)701から出力された一方の光は、発振波長を変更する手段(以下、波長選択素子)702に入射する。ここでの波長選択素子702は、光偏向器703および回折格子704で構成される。第1の実施形態では、回折格子はリットマン型の配置としていたが、本実施形態の波長掃引光源では、回折格子704は、回折光が入射光と同軸逆方向に回折されるリトロー条件で配置されている。波長選択素子702としては、ファブリペローフィルタなどの別の波長選択原理に基づく素子を用いても良い。 
 光偏向器703は、ガルバノミラー、MEMSミラー、KTN光偏向器などを用いることができる。SOA701のもう一方の出力光は、光帰還素子105において一部が反射され、回折格子704および光帰還素子705の間で構成される共振器に帰還される。光帰還素子705を透過した光は、そのまま出力光としてレンズ706で集光され、光ファイバ707に結合され、光源の出力光712となる。図11の波長掃引光源における共振器は直線型の構成であるが、これは一例であり、それぞれの光学素子が光ファイバで接続されたリング型の共振器で構成されても良い。
 波長掃引光源の共振器において出力光の線幅を狭窄化する動作は、以下のように説明される。SOA701から出力された光は広い光スペクトルを有しており、波長選択素子702においてフィルタリングされ、波長選択素子702が持っているフィルタの形状によって切り出されたスペクトルが、再びSOA701に入射する。この際にSOA710で増幅され、SOA701の反対側の端面から出力さる。光帰還素子705により、光の一部は再びSOA701に戻され、このようにして光共振器内をある一定時間の間に複数回、周回する。波長掃引光源では、波長選択素子702の選択波長が変化するため、この一定時間はある特定の波長をフィルタスペクトルが横切る時間で決まる。この横切る時間内に複数回、波長選択素子102を光が通過することでフィルタリングの累積効果により出力光の線幅が狭窄化される。
 図10でも説明したように、波長掃引光源では、使用するSOA701の利得スペクトルや共振器損失の波長依存性などの影響によって、レーザ発振の閾値に波長依存性が生じる。そのため、コヒーレンス長を長く維持するために最適なSOA駆動電流値は、波長ごとに異なる。本実施形態の波長掃引光源では、波長ごとに異なる最適なSOA駆動電流値にSOAを設定するために、光源からの光出力パワーを監視し、制御する。
 波長掃引光源からの光出力パワーは、SOA701で増幅を受ける回数に依存する。一方で、既に述べたようにコヒーレンス長も波長選択素子702を通過する回数に依存する。このため、光出力パワーおよびコヒーレンス長は、それぞれ共振器における発振光の周回数に関連付けられている。波長選択素子702の透過特性がガウス形状のスペクトルを持っているとすると、その形状関数F(λ)は次式で表される。

Figure JPOXMLDOC01-appb-I000001
ここで、ηは波長選択素子702の最大透過率を、λは中心波長を、Δλはスペクトルの半値全幅を表している。光が光共振器を一周回したときに波長選択素子702以外の部分で生じる正味の利得Gnet(ネットゲイン)と、波長選択素子702で生じる損失とが釣り合う状態が、レーザの発振閾値の状態となる。光共振器の損失は各構成素子の特性で決まっているため、ネットゲインGnetは主にSOA701の利得で変化する。
 波長選択素子702における損失Lは、dBで表すと次式となる。

Figure JPOXMLDOC01-appb-I000002
発振閾値の条件から、ネットゲインGnetと損失Lとの間で次式の関係が求められる。
net+L=0            式(3)
さらに式(3)に式(1)および式(2)を代入すると、以下の関係式が求められる。

Figure JPOXMLDOC01-appb-I000003
ここで、式(4)の左辺(λ-λ)は、ネットゲインGnetが損失Lよりも大きい状態であって発振閾値を上回っている波長帯域の半幅を表している。したがって、発振閾値を上回る波長領域の全幅Δλeffは、式(4)の半幅を2倍して、式(5)を組み合わせ、次式となる。

Figure JPOXMLDOC01-appb-I000004
 式(6)のΔλeffは、レーザ発振に関わる実効的な波長選択幅を表しており、実効波長選択幅とも呼ばれる。ここで注目すべきは、波長選択素子702のスペクトル幅Δλが波長によらず不変であるとみなせる場合(条件1)には、レーザ発振に関わる実効波長選択幅ΔλeffはネットゲインGnetつまりSOA701の駆動電流によって決定されることである。これは、各波長においてネットゲインGnetを一定にすることによって、実効波長選択幅Δλeffを一定に維持できることを意味する。さらに波長選択素子702の選択波長が掃引波長帯域内において一定の速度で変化すると仮定(条件2)すると、実効波長選択幅Δλeffが一定であれば、各波長において共振器内における光の周回数も一定となる。
 上述の2つの条件が成り立てば、共振器内における光の周回数、ネットゲインGnetおよび実効波長選択幅Δλeffはいずれも各波長で一定となる。すなわち、波長掃引光源の光出力パワーおよびコヒーレンス長は、掃引帯域内の波長によらず一定となる。したがって、監視信号光として光出力パワーをモニタしながら、ネットゲインGnetを一定化するようにSOA駆動電流を制御することで、掃引波長帯域内の各波長でコヒーレンス長を一定にすることができる。そのため光偏向器703の制御電圧は、選択波長が掃引波長帯域内において一定の速度で変化する波形とすることが好ましい
 しかし実際には波長選択素子102は、制御信号(波長掃引信号)として正弦波によって駆動されることもある。このような場合は、掃引速度が波長によって異なり、実効波長選択幅Δλeffが一定であったとしても共振器内における光の周回数は波長によって異なる。したがって、光出力パワーを一定化したとしてもコヒーレンス長は必ずしも一定になるとは限らない。しかしながら、掃引波長帯域の中心付近の波長選択素子702の選択波長の変化速度は概ね一定であり、正弦波によって駆動しても、掃引波長帯域の中心付近では一定程度はコヒーレンス長が一定に近い状態が得られる。上記の議論は波長空間で表されているが、これらを周波数空間に置き換えたとしても本発明の効果は同様に得ることができる。
 以上のように、ネットゲインGnetを発振閾値に近づけ光出力パワーを一定化する制御を行うことで、実効波長選択幅Δλeffをさらに狭くし、各波長において、図10に示したように発振閾値状態の近傍でコヒーレンス長を最大化した状態とすることができる。図10のグラフでは発振閾値の状態(点線)に近づくにつれてコヒーレンス長が指数関数的に増大するように見えるが、実際には発振閾値に近づくにつれて発振状態が不安定になる。このため、レーザ縦モード間のモード競合によりコヒーレンス長はかえって劣化する。この不安定性も考慮すれば、ネットゲインGnetやSOA駆動電流値の最適値は、発振閾値よりわずかに高い値(図10の点線の右側)となる。
 光出力パワーの一定化を行うためには、光源の実際の動作である波長掃引に対応して、SOA駆動電流値を刻々と時間的に(対応する波長と同期して)増減させるよう制御する必要がある。本実施形態の波長掃引光源は、この制御回路を含む波長掃引光源装置として実施できる。
 図12は、波長掃引光源およびSOAの注入電流を制御して波長掃引出力光のコヒーレンス長を常に長い状態に維持する制御回路を含む波長掃引光源装置の構成を示す。波長掃引光源装置720は、図11に示した波長掃引光源700を含み、さらに波長掃引光源700からの出力光の一部を分岐して取り出す光分岐素子722、一方の分岐光728の光出力パワーを測定する光検出器723、光検出器723の出力信号をデジタル信号に変換するA/D変換器724、デジタル信号を処理してSOAのドライバ装置708への制御信号となる制御波形データを生成する演算機(例えばCPU)725、演算機で生成した波形信号711を出力する任意波形発生器726を備えている。
 波長掃引光源700からの出力光を光分岐素子722によって一部の光を取出して一方の分岐光は光源出力光727として、他方の分岐光728は監視信号光として光検出器723で検出する。光検出器723では、監視信号光は電圧信号729に変換されA/D変換器724によってサンプリングする。サンプリングされた信号は演算機725に入力され、演算機725では入力された信号と後述する設定値730とを比較し、エラー信号を生成して、エラー信号に応じた制御波形731を生成する。設定値730は、後述の波長毎の閾値電流や、初期電流、許容誤差D、目標光出力パワー値Vなどが含まれ、図示していないメモリや記憶装置から読み出しても良い。制御波形731は任意波形発生器726に入力され、任意波形発生器726は入力波形のデータ列に対し制御信号711としてデータ点間が滑らかに補間された電圧波形の形態で出力する。波長掃引光源700内の電流ドライバ708は、任意波形発生器726から入力される電圧波形をSOA駆動電流波形に変換する機能を有しており、入力電圧波形に対応したSOA駆動電流をSOA701に注入する。発振波長を変更する手段(例えば光偏向器)702を駆動するドライバ装置709は、波長掃引に同期したトリガ信号710を任意波形発生器726に供給する。
 上述の説明では、演算機725を、D/A変換機能を含み制御波形731を出力するものとして説明したが、演算機725から波形データを出力して、別に設けられたD/A変換器によって制御波形731を得ても良い。A/D変換器724から、任意波形発生器726までの構成は、図12に示した構成だけに限られず様々な変形が可能である。次に演算機725によって実施される、2分法による時間-SOA駆動電流波形の生成方法について詳細に説明する。後述の時間-SOA駆動電流波形の生成方法の手順は、演算機725またはCPU等によるソフトウェア演算処理によって実現可能な処理プロセスを示す。図12の構成では、光共振器およびSOAの特性が決まっており、後述する光出力パワーの目標値Vが既知であれば、プロセッサ以外に、第1の実施形態における干渉計などのような特別なハードウェア要素なしに、ソフトウェア演算処理として実現できる。
 図13は、本実施形態の波長掃引光源におけるSOAへの注入電流および光出力レベルの関係を説明する図である。図13のグラフでは、横軸にSOAへの注入電流すなわちSOA駆動電流を示し、縦軸には図12の光検出器723から得られる光検出器出力電圧を示す。光検出器出力電圧は、波長掃引光源700からの光出力レベルに対応する。第2の実施形態の波長掃引光源に特有の、波長掃引動作を実施しながらSOA駆動電流を最適化する制御は、以下のように説明される。
 ここで、1回の波長掃引動作においてA/D変換器724でサンプリングすることで取得できるN個の光出力パワーデータV、V、…、VN-1を考える。このN個の各データ点は、全掃引波長帯域に含まれる波長λ、λ、…、λN-1の光出力パワーデータにそれぞれ対応している。この光出力パワーデータの各々は、図12における監視信号光の電圧信号729に一対一に対応し、図13の縦軸の光検出器出力の電圧に対応する。ここで、N個の波長の内のn番目のデータ点803、つまり波長λにおける光出力パワーデータVについて注目する。
 図13において、実線805は波長λにおけるSOA注入電流対光出力パワー特性(I-P特性)を示している。また点線806は、光出力パワーデータの異なる例として、波長λn-1におけるSOA注入電流対光出力パワー特性(I-P特性)を示している。目標となる光出力パワー値804をVとし、V803が得られたときの注入電流値をIとする。現在の光出力パワーVと目標光出力パワー値Vとの差を、エラー信号Eとして次式のように定義する。
 E=V-V          式(7) 
 目標光出力パワー値Vを、光出力パワーの制御に用いる数値として、目標光出力パワー値Vを任意に設定することが可能である。目標光出力パワー値Vの最適値は、波長掃引光源700のそれぞれの構成要素のゲイン特性、損失特性によって決まるため、光源ごとに存在するはずである。目標光出力パワー値Vの最適値を求めるに際しては、波長掃引光源としての実際の動作をさせる前に、以下に説明する制御フローを使用してあらかじめVをある値に設定して光出力パワーを一定化する制御を行い、さらにこのVの値を変えてこの制御を繰返し実行し、最適値を探索する。このとき、第1の実施形態で示した手法などによってコヒーレンス長を測定することで探索することができる。
 このように決定された最適な目標光出力パワー値Vについては、波長掃引光源700の出力特性が著しく変化しない限りは変化しないと考えられる。このため、波長掃引光源としての実際の動作をさせる前に一旦目標光出力パワー値Vを決定してしまえば、以後は同じ値を用いることができる。したがって、例えば波長掃引光源の製造時に、波長掃引光源装置の製造調整のための外部調整設備として光干渉計などを備えておいて、V値の決定のためにコヒーレンス長を測定できれば十分である。装置の製造・調整工程で目標光出力パワー値Vの決定が行われ、個々の装置にこのV値を記憶させておけば、第1の実施形態のように波長掃引光源装置自体に光干渉計を備える必要はない。
 原理的には各波長において図13に示したI-P特性をそれぞれ取得し、得られたI-P特性および目標光出力パワー値Vから各波長でコヒーレンス長を最大化するための所望の光出力パワーが得られるSOA注入電流値を決定すれば良いはずである。しかしながら実際には、ある時間におけるSOA駆動電流値に対する光出力パワーは、波長掃引において時間的に前後する瞬間に注入されたSOA駆動電流値によっても変化する。したがって波長掃引に対応したSOA駆動電流波形を生成するためには、フィードバック制御が必要である。SOAへの注入電流が発振閾値電流以下では光出力パワーは得られないので、注入電流値の下限801は発振閾値電流Ith、nとなる。SOA駆動電流の最適化プロセスで、初期に注入する電流値Iini、nは、発振閾値電流に係数aを乗算した値として次式の通り設定する。
 Iini,n=aIth,n       式(8)
 上式で係数aの値は、各波長において初期電流を注入したときにエラー信号Eが負の値となるように設定する。つまり、N個の光出力パワーデータV、V、…、VN-1の全てのデータ点において、初期電流を注入時に最初の光出力パワーデータが目標光出力パワー値Vよりも大きい、V>Vの関係が満たされることになる。次に、SOA駆動電流をフィードバック制御するときの制御ステップΔIを、次式のように定義する。


Figure JPOXMLDOC01-appb-I000005

 また光検出器723によってモニタしている光出力パワー値Vと目標光出力パワー値Vとの間の許容誤差をDと定義する。SOA駆動電流の最適化プロセスにおいて最初のフィードバック制御を行う時、波長λにおけるSOA電流値Iは式(8)によって定義した初期電流値Iini,nとなる。SOA駆動電流Iが設定され、波長掃引中のある時点でモニタされた光出力パワーVから求められるエラー信号Eの絶対値が、定められた許容誤差Dより大きいとき、最適化プロセスはIの値を次の波長掃引のために更新するよう制御する必要がある。E<0(0も含む)のときは検出された出力パワーが目標値よりも小さくV>V(0も含む)となるので、注入電流Iは減らすべきであり、次式のようにIの値を更新する。


Figure JPOXMLDOC01-appb-I000006

 上式で、Mはフィードバック制御がM回目のものであることを表している。また、‘:=’は、右辺の数値を左辺の変数に代入することを示している。上の説明ではE<0(0も含む)の場合を示したが、E>0(0も含む)の場合、式(10)の中のΔI/2の符号は正となる。制御ループでは、|E|<Dとなった場合には電流値Iを更新しない。全ての波長においてエラー信号Eの大きさが許容誤差D以下となったとき、SOA駆動電流の最適化プロセスの制御は停止する。
 したがって、本実施形態の波長掃引光源では、光利得制御手段(723~726)は、レーザ共振器からの出力光レベルを検出し、前記検出された出力光レベルに基づいて、前記出力光レベルが、所定の出力光レベル値になるような光利得制御信号を生成し、当該光利得制御信号によって光利得を変更する手段710を制御するよう構成されている。
 上述の光出力パワーVからSOA駆動電流を制御するプロセスは、ある波長λにおける一連の制御として説明しているが、本実施形態の波長掃引光源装置720では、1回の波長掃引期間中に、上述のN個の光出力パワーデータV、V、…、VN-1を取得することができる。そして、|E|<Dの条件が成り立つまでは、条件を満たさない波長においてSOA駆動電流値Iの更新を実行することになる。つまり波長掃引動作を繰り返すたびに、図13に示したSOA駆動電流の最適化プロセスによって、各波長において目標となる光出力パワー値Vを得るためのSOA電流値Iの更新の制御が繰り返される。
 図14は、本実施形態の波長掃引光源装置720におけるSOA駆動電流の最適化プロセスを波長掃引の時間も含めて概念的に説明する図である。図14の上側にはSOA駆動電流の最適化プロセスの一連の進行を、時間-SOA駆動電流の時間波形812、813、814として示している。また図14の下側には、波長掃引時間を揃えて時間-SOA駆動電流の時間波形と連動させて、波長掃引光源からの出力光の時間-波長波形820を示している。時間-SOA駆動電流の時間波形の内、一番下の一点鎖線811は、SOAの発振閾値電流の時間変化を示している。図14の下側のグラフで横軸の時間は光源出力光の波長と一対一に対応しているので、各波長に対して発振閾値電流も一対一に対応していることになる。一番上の一点鎖線812は、式(8)で定義された時間-初期電流値Iini、nを示している。初期電流値Iini、n812上には、ある波長間隔で、最適化プロセスの開始時に最初に設定される初期電流値Iini、nに対応した複数のサンプル点が示されている。これらのサンプル点が取得される波長間隔は、等間隔でも非等間隔でも良い。
 時間-SOA駆動電流の時間波形の内、コヒーレンス長を最大化する光出力レベルに対応する目標値の最適SOA駆動電流波形815が、実線で示されている。初期電流値Iini、nの設定の後で、各波長において1回目のフィードバック制御が行われると新しい時間-SOA駆動電流波形813(1回目の更新)が得られる。すなわち、全波長において式(4)によってSOA駆動電流値が更新され、時間-SOA駆動電流813の時間波形によってSOAを駆動しながら次の2回目の波長掃引が行われる。2回目の波長掃引の間に光検出器723から2回目の新たな波長λ、λ、…、λN-1の光出力パワーデータを取得することができる。新たに全波長の光出力パワーデータが取得され2回目の制御が終わると、更新された時間-SOA駆動電流814(2回目の更新)の時間波形が得られる。3回目の波長掃引のために更新された時間-SOA駆動電流814が使用される。時間-SOA駆動電流波形814は、時間-SOA駆動電流波形813と比べて、各波長においてより時間―最適SOA駆動電流波形815に近づいている。
 また図14の例では、中央付近の波長帯において、電流時間波形813と電流時間波形814が一致している。これは、この中央部の波長帯域では、1回目のSOA駆動電流値Iの更新後で|E|<Dの条件が成立したため、2回目のIの更新が不要であったことを示している。図14に示したように、SOA駆動電流-時間波形は、各波長において、例えば掃引開始波長λにおけるサンプル点817→818→819のように、図13のSOA駆動電流の最適化プロセスが実施されて、時間―最適SOA駆動電流波形815に近づいていくよう動作する。したがってある波長に着目すれば、SOA駆動電流を制御する図14のプロセスは、波長掃引に連動または同期して、時間的に離散的な制御を実施していることになる。
 図15は、全ての波長を含むSOA駆動電流の最適化プロセスの制御フローチャートである。S901からS905において、各波長点における初期データが設定され、S906からS919で、各波長においてSOA電流の更新(S907)およびモニタ電圧取得(S908)が行われる。S911からS917の繰り返しループによって、N個の波長λ、λ、…、λN-1に対して、エラー信号Eの許容誤差との比較、SOA電流値Iの更新値の決定などの上述の動作が順次行われている。
 図15のフローチャートで、配列Cの要素Cnは波長λにおいて制御電流値Iを更新する必要かどうかを表す。制御が必要な場合Cn=Faultとなる。AND(C)は配列Cの全配列要素のANDを表しておりAND(C)=Trueの場合、SOA駆動電流の最適化プロセス制御は停止する。
 図16は、第2の実施形態の波長掃引光源において、実際のSOA駆動電流の時間波形および光出力レベルを従来技術と比較して示した図である。図16の上側には、従来技術のSOA駆動電流一定の場合950、本発明によってSOA駆動電流の最適化プロセスを実行後の場合951の時間―SOA電流波形をそれぞれ示した。また図16の下側には、対応する従来技術の場合952、SOA駆動電流の最適化プロセスを実行後の場合953で、時間-光出力レベル(光検出器の出力信号)をそれぞれ示した。SOA電流注入は、波長掃引の短波長から長波長への掃引区間のみであり、図16からわかるように、波長掃引区間において、光出力レベルは波長変化(時間変化)に対してほぼ一定となっている。
 図17は、第2の実施形態の波長掃引光源の波長掃引光の干渉フリンジ波形を従来技術と比較して示した図である。振幅を規格化した光路長差1mmのマイケルソン干渉計により取得した干渉フリンジ波形を、従来技術のSOA駆動電流一定の場合954、本発明によってSOA駆動電流の最適化プロセスを実行後の場合955で、それぞれ示している。図17のフリンジ波形955の振幅も、SOA駆動電流一定の場合954と比較して、時間に対して一定となっている。
 図18は、SOA電流の制御により一定化された光出力パワーの値とコヒーレンス長の関係を示す図である。第2の実施形態の波長掃引光源で一定に制御する光出力パワーの値を変えて、それぞれの光出力パワーの場合のコヒーレンス長を求めている。図18では、おおよそ0.35mW付近でコヒーレンス長が最大になることがわかる。このときのコヒーレンス長は、80nmの波長掃引幅に渡って概ね約24mmとなっており、従来技術の一定のSOA駆動電流でSOAを動作させた場合のコヒーレンス長最大値である12mmの約2倍の値が得られた。しかも、従来技術の場合ではコヒーレンス長最大となる波長は狭い範囲に限られていた。図18からも明らかなように、本実施形態の波長掃引光源により、従来技術の場合と比べて2倍ものコヒーレンス長を得ることができる。
 第2の実施形態の波長掃引光源では、演算機725における時間-SOA駆動電流波形の生成法として2分法を用いた。2分法では、光利得制御手段(723~726)は、波長掃引に連動して時間的に離散的な制御を実施し、前記出力光レベルと前記所定の出力光レベルとの差分に対して、前記光利得を変更する制御量(SOA駆動電流の制御ステップ)を漸次半減している。この2分法は、波長掃引光源の環境変化に対して制御対象量(光出力レベル)が比較的安定している場合に適した簡便なフィードバック制御方法である。
 一方で温度変化などの環境変化に対して光源特性がより敏感な場合では、光出力レベルを目標値に収束させることができない場合もあり得る。そのような場合には、演算機725によって、2分法に代えて、掃引波長帯域の各波長λについて独立にPID制御など行うことで解決することもできる。PID法によれば、光利得制御手段(723~726)は、波長掃引に連動して時間的に離散的な制御を実施し、前記出力光レベルと前記所定の出力光レベルとの差分に対して、前記光利得を変更する制御(SOA駆動電流の操作)を、当該差分、前記差分の積分、および前記差分の微分の要素によって行う。
 図14でも説明したように、本実施形態のSOA駆動電流の最適化プロセスでは、1回の波長掃引期間の間に、全波長におけるN個の光出力パワーデータV、V、…、VN-1を取得して、SOA駆動電流Iの更新を行って、次の波長掃引期間で、新たな全波長のN個の光出力パワーデータを取ることが可能である。通常、10~100回程度のSOA駆動電流Iの更新を実行すれば、目標の光出力レベルに達することができる。しかしながら、演算機725における演算処理能力・処理時間などを考慮して、所定の掃引回数毎にフィードバック制御を行えば良い。
 波長掃引光源における掃引の繰り返し周波数は、200kHz以上にも達する場合があるので、波長掃引光源の実際の利用態様からすれば、10~100回程度のSOA駆動電流Iの更新は非常に短い時間に過ぎない。したがって、光源の波長掃引開始の度に毎回、連続的にフィードバック制御をする必要はない。例えば、図15のフローのS906~S919のステップを、100、1000または10000回の掃引毎に間隔を置いて実施しても良い。この実施間隔は、一定間隔であっても良いし不定期のものであっても良い。例えば、波長掃引光源装置のウオームアップの時間で、図14および図15のSOA駆動電流の最適化プロセスが完了すれば十分である。
 SOA駆動電流の最適化プロセスが一旦完了すれば、波長掃引光源の動作が安定している環境の場合であれば、図15の一連のフローは停止させたままでも良い。一方、環境変化が大きくて、光源の状態に変化が生じやすい場合には、定期的にまたは不定期に、図15の一連のフローを実施することもできる。本実実施形態のSOA駆動電流の最適化プロセスは、第1の実施形態のように、コヒーレンス長の計算を行わない。光源の光出力レベルをモニタしながら、SOA電流に直接フィードバックして、一定の目標光出力パワー値(V)に制御するだけなので、最適化プロセスの動作は第1の実施形態と比べて非常に簡単である。波長掃引光源装置720の演算機725は、既存のCPUなどを利用すれば良い。装置の製造調整の段階で一定の目標光出力パワー値(V)が事前に決定されていれば、波長掃引光源装置の通常動作の波長掃引のために第1の実施形態のような干渉計は不要である。より簡素な構成要素と処理によって、第1の実施形態と同様に、波長掃引期間中の波長掃引出力光のコヒーレンス長を常に長い状態に維持することができる。
 尚、本実施形態の波長掃引光源では、図10に示したようなSOA電流とコヒーレンス長の関係が成り立つ場合を前提として、レーザ発振閾値の状態の光利得よりもわずかに高い光利得値となるように、光利得を変更する手段(SOA)を制御するものとして説明した。しかしながら、レーザ共振器の構成やタイプによっては、レーザ発振閾値の近傍と比べて大幅に高い光利得値において、コヒーレンス長が最大となるような例外的な場合もあり得る。そのような場合でも、コヒーレンス長が最大となるような目標光出力パワー値(V)を実際のSOA電流-コヒーレンス長特性に基づいて決定すれば良いのであって、図12の波長掃引光源装置の構成や図15の制御フローに変わりはない。本実施形態の波長掃引光源によれば、SS-OCT装置におけるPSF強度を上げて、SS-OCT装置のイメージ検出感度を高くし、イメージ品質を向上させることができる。
[第3の実施形態]
 上述の2つの実施形態では、コヒーレンス長を最大に維持するSOA電流を求め、SOA電流波形データを生成して、波長掃引光源を動作させる構成および制御方法を示してきた。波長掃引光源において、出力光のコヒーレンス長に影響を与える重要な要素は、図1に示した光偏向器106や図11に示した光偏向器703である。波長掃引光源に好適な光偏向器としては、電気光学結晶を用いた光偏向器がある。電気光学結晶を用いた光偏向器は、機械的な可動部分が無く、電気光学効果がMHz~GHzの範囲の高周波数で応答するため、高速な光偏向が可能である。以下では、本発明の波長掃引光源として好適に利用可能な光偏向器の様々な実施例を示す。
 電気光学結晶であるKTN(KTa1-xNb)結晶またはKLTN結晶(K1-yLiTa1-xNb)(以下、まとめてKTN結晶という)は、低い電圧で屈折率が大きく変わる電気光学効果が大きい物質として知られている。さらに、電極としてTi、Cr材料を用いると、結晶内に電荷を注入することができ、注入された電荷によって生じる内部電界を利用して、高速・広角な光偏向器を実現することができる。従って、レンズ、プリズム、ミラーといった一般的な光学部品を、高速で動く必要がある用途に用いる場合に、これら光学部品に代えて、KTN結晶を用いることができる。
 図19は、従来技術のKTN結晶を用いた光偏向器の構成を示す。KTN結晶1101の対向する2つ面のそれぞれに、電極1102、1103が形成されている。2つの電極1102,1103の間には、制御電圧源1104から電圧が印加される。入射光1105は、KTN結晶1101の電極1102、1103が形成された面と直行する面(xy平面)に入射され、KTN結晶1101内をz軸方向に伝搬する。KTN結晶1101内において偏向を受けて、x軸方向に進行方向を変えた出射光1106が得られる。制御電圧源1104からの印加電圧に応じた偏向角が得られる。制御電圧源1104からは、光偏向器の用途に応じた電圧が与えられる。例えば、印加電圧の波形は、正弦波、鋸波などであり、光偏向器の用途に応じて印加される。適切な最大偏向角を得るために、KTN結晶1101へは、概ね数百V程度の電圧を印加する。
 高周波動作時には、印加電圧によって電極から注入された電子の移動距離が電極間の距離より短くなるため、理想的な空間電荷制御状態が実現されず、偏向角が減少する。この問題点に対しては、偏向を生じさせるための印加電圧を印加する前に、KTN結晶1101にバースト状の波形の電圧を印加することによって、結晶中へ電子を注入し、予めトラップ準位に電子を捕捉させる制御法が提案されている(例えば、特許文献2参照)。
 KTN結晶1101中のトラップに、予め電子を充填することによって、偏向を生じさせるための電圧の印加に際しては、KTN結晶1101中に電界の分布または傾斜を生じさせることが可能となり、光偏向を実現することができる。トラップ準位に電子を捕捉させる方法として、DCバイアス電圧に光偏向や走査のための交流電圧を重畳して、KTN結晶1101に印加する方法も提案されている(例えば、特許文献3参照)。
 DCバイアス電圧を印加し続けると、トラップ準位から熱励起等で再放出された電子を補うように電子が負電極から供給されるため、長時間経過しても一定の電子が結晶内に残留する。従って、高周波動作に必要な電子の捕捉量を一定に保持し、長時間経過しても偏向角が減少しない偏向器を実現することができる。
 KTN結晶内の電荷密度Nは、結晶の誘電率が高い場合には空間的に一様に近づくが、誘電率が低い場合には電圧印加時の陰極付近のみに蓄積され、結晶内で不均一な分布となる。電荷密度Nが空間的に一様ではなく、N(x,y,z)のように結晶内の位置に依存するような分布の場合には、結晶内で偏向される光の透過位置によって、電気光学結晶のレンズとしての効果が異なることになる。つまり、光偏向のための駆動電圧の瞬時値によって、レンズ効果が異なることになる。
 これに対し、光偏向のための電圧に同期して焦点距離が変化する可変焦点レンズをKTN結晶の出射側に配置し、KTN結晶への印加電圧に依存せず、集光特性が一定な偏向光を得られる方法が提案されている。この方法では、本来の光偏向器の他に動的に動作する付加的な可変焦点レンズが必要となり、可変焦点レンズを動作させる電源、光偏向器と可変焦点レンズを同期させるための電気回路などの構成が複雑になってしまう問題があった。
 そこで、本実施形態の波長掃引光源に好適な光偏向器は、レンズ効果を有する光偏向器において、レンズ効果の空間的変化を補償するよう動作する。本実施形態の波長掃引光源に好適な光偏向器は、電気光学結晶の対向する面に少なくとも2つの電極が形成された光偏向器であって、前記少なくとも2つの電極に電圧を印加すると、電気光学効果により前記電気光学結晶内に一様でない屈折率分布が生じることにより、前記電圧により形成された電界に直交する光の進路を屈曲させる光偏向器と、前記電気光学結晶の入射側と出射側の少なくとも一方に配置された光学素子であって、前記光が偏向する範囲にわたって空間的に集光特性が分布している光学素子とを備える。本実施形態の光偏向器によれば、光偏向器の後段に集光特性や光路長が空間的に分布している光学素子を配置することにより、レンズ効果の偏向方向における変化を補償することができる。
 以下、図面を参照しながら本発明の第3の実施形態の波長掃引光源に好適な光偏向器の様々な構成について詳細に説明する。
 図20は、本発明の一実施形態にかかる波長掃引光源に好適な光偏向器の構成を示す図である。電気光学効果を利用した光偏向器は、電気光学結晶1201の対向する2つ面のそれぞれに、電極1202、1203が形成されている。制御電圧源1204から、この電極対を介して、電気光学結晶1201内への電荷の注入、または光の偏向のための電圧を印加する。
 電気光学結晶として、例えばタンタル酸ニオブ酸カリウム(KTa1-xNb(0<x<1):以下、KTNという)、またはリチウムをドープした結晶(K1-yLiTa1-xNb(0<x<1、0<y<1):以下、KLTNという)を用いることができる。これらの電気光学結晶は、電圧印加による電界に伴って、結晶に電荷の注入が行なわれる。その結果、結晶内に、注入電荷により形成される空間電荷分布、または、注入電荷がさらに電気光学結晶中に捕捉されて生成されるトラップ電荷分布が生じる。この電荷分布によって一様でない電界分布が形成され、電界分布の変化に応じた屈折率分布が生じる。この一様でない電界分布が屈折率の勾配を発現させ、屈折率の勾配に直交する光線の進路を屈曲させる。このとき、電気光学結晶の中央付近で屈折率が高い非線形な屈折率分布となるので、電気光学結晶1201を透過した入射光1205は、偏向されると同時に集光されて出射光1206または出射光1208となり、焦点以降は発散する。
 光偏向器として用いる電気光学結晶1201の組成分布や、電荷充填のための印加電圧によっては、トラップ密度Nが常に空間的に一様とは限らない。この場合には、結晶内で偏向される光の透過位置によって、電気光学結晶のレンズ効果が異なる。つまり、光偏向器である電気光学結晶への印加電圧の瞬時値によって、電気光学結晶によって形成される焦点位置が異なることになる。この変動を補償するような空間的に集光特性が分布している光学素子1209によって、光偏向器への印加電圧の瞬時値に関わらず、一定の焦点距離を有する偏向動作が得られる。
 例えば、空間的に集光特性が分布している光学素子1209がない場合には、偏向時の出射光のビーム形状は、角度によって焦点位置が異なる出射光1206または出射光1208となる。このとき、出射光1208は、出射光1206に比べて焦点位置が電気光学結晶1201の出射面に近い。これを、光偏向器1201の後段に配置する光学素子1209によって、出射光のビームが発散するような波面変化を与え、出射光1206と同じ焦点位置を有するような出射光1210に変換する。または、出射光1206の焦点位置が電気光学結晶1201の出射面に近くなるような光学素子1209の配置によって、出射光1206の焦点位置を出射光1208と同様の出射光1211に変換する。
 光偏向器として用いる電気光学結晶1201が偏向角により集光特性が変化しているとき、出射光のビームの発散角を補正する光学素子1209は、出射光が偏向する範囲にわたって、空間的に光学特性が分布している。このようなビームの補正によって、出射光はどの偏向角においても等しい焦点位置を有することができる。例えば、レーザプリンタ等では感光体面上のどの位置においてもビームウェストを形成できるという利点が生じる。また、出射光をどの偏向角においても平行光となるようにすることもでき、例えば、分光装置などの応用では、分光帯域全域において分解能が劣化することなく、一様な分解能の分光特性を得ることができる。
 次に、上述の本実施形態の波長掃引光源に好適な光偏向器に関して、より具体的な実施例を述べる。
 図21は、実施例4にかかる波長掃引光源に好適な光偏向器の構成を示す。KTN、KLTNなどの電気光学結晶1301を備える光偏向器は、電気光学結晶1301の対向する2つの面のそれぞれに、電極1302、1303が形成されている。電気光学結晶1301に入射される入射光1305は、シリンドリカル凹レンズ1307を透過した後に、電気光学結晶1301に入射される。平行光である入射光1305は、このシリンドリカル凹レンズ1307によって、電気光学結晶1301へ拡散しながら入射される。入射光1305は、電気光学結晶1301のレンズ効果によって、出射光として出射される。
 電気光学結晶1301は、厚さ1.2mmのKTN結晶を用い(電極間隔も1.2mmとなる)、ビーム直径1.0mmの光線を、電界に平行な直線偏光で入射した。電気光学結晶1301の長さは4mmである。電気光学結晶1301の入出射端面は、反射膜と反射防止膜の両方が形成されており、折り返し光路により入射光が偏向される距離は12mmとなっている。
 電気光学結晶1301に制御電圧源1304により、振幅±360V、周波数200kHzの正弦波状のAC電圧を、-240VのDCバイアス電圧に重畳して印加する。電圧は結晶表面に蒸着された金属からなる電極1302、1303を介して印加されている。このとき、110mrad(約6.5°)の偏向角が得られる。DCバイアス電圧を印加することによって、結晶内のトラップ密度は、空間的に一様ではない。結晶内でのトラップ密度が一様でないために、結晶内で偏向される光の透過位置によって、電気光学結晶のレンズとしての効果が異なる。プリズム1309がない場合には、出射光が偏向した時のビームは、偏向角によって焦点位置が異なり、出射光1306、1308のような形状である。
 実施例4では、電気光学結晶1301の後段に、SF11ガラスを硝材とし、入射面と出射面のなす角が30°のプリズム1309を配置した。プリズム1309は透過する光のビーム径を拡大するので、透過後のビームは、プリズム1309に入射される前に比べて、回折による広がりが抑えられる。このプリズム1309によるビーム径の拡大率は、入射角に依存する。
 図22は、実施例4におけるプリズムのビーム径倍率の入射角依存性示す。実施例4では、偏向されている状態のうち、焦点位置が電気光学結晶1301の出射面に近い偏向角(電圧)のとき(図21において出射光1308の場合)、プリズム1309への入射角が大きくなり、焦点位置が遠い偏向角(電圧)のとき(図21において出射光1306の場合)、プリズムへの入射角が小さくなるように配置している。さらに、プリズム1309の頂角、つまりプリズムの幅が小さくなる方向が、焦点位置が遠い偏向角の側になるように配置する。
 この配置によって、焦点位置が近い偏向角(電圧)の出射光1308は、焦点距離が大きい偏向角(電圧)の出射光1306に比べて、ビームがより拡大され、回折による広がりを抑制する効果がより得られ、出射光1310となる。一方、出射光1306は、プリズム1309によって、出射光1311となる。このように出射光が偏向する範囲にわたって空間的に集光特性が分布しているプリズム1309によって、焦点位置の角度依存性が抑制され、制御電圧源から光偏向器に印加されるAC電圧の瞬時値に関わらず、一定の焦点位置を有する偏向動作が得られる。
 図23は、実施例5にかかる波長掃引光源に好適な光偏向器の構成を示す。KTN、KLTNなどの電気光学結晶1401を備える光偏向器は、電気光学結晶1401の対向する2つの面のそれぞれに、電極1402、1403が形成されている。電気光学結晶1401に入射される入射光1405は、シリンドリカル凹レンズ1407を透過した後に、電気光学結晶1401に入射される。電気光学結晶1401の形状、駆動条件は、実施例4と同じである。
 実施例5では、電気光学結晶1401の後段に、自由曲面ミラー1409を用いることにより、光偏向器によって発現される偏向角に依存した集光効果を補正している。自由曲面ミラー1409は、偏向方向に沿って偏向方向の曲率半径が変化している凹面ミラーであり、その曲率は光偏向器の接地側の電極1402から制御電圧源側の電極1403に向かって小さくなっている。つまり、自由曲面ミラー1409の曲率は、偏向された出射光のうち、焦点位置が電気光学結晶1401の出射面に近い偏向角(電圧)のときの出射光が照射される面の曲率半径は小さく、焦点位置が遠い偏向角(電圧)のときの出射光が照射される面の曲率半径は大きくなっている。
 図24は、実施例5における電気光学結晶1401のみによる焦点距離の印加電圧依存性を示し、図25は、実施例5における自由曲面ミラーの焦点距離と印加電圧の関係を示す。電気光学結晶1401の焦点距離が短い電圧(焦点位置が電気光学結晶1401の出射面に近い電圧)のときの出射光が照射される凹面の曲率は小さくし、自由曲面ミラーの焦点距離を長くして印加電圧に対する焦点距離の変動を補償する。一方、焦点距離が長い電圧(焦点位置が出射面から遠い電圧)のときの出射光が照射される凹面の曲率は大きくし、自由曲面ミラーの焦点距離を短くして印加電圧に対する焦点距離の変動を補償する。このようにして、印加電圧に対する焦点距離の変動を補償するように、自由曲面ミラー1409の曲率が設定されていることがわかる。
 自由曲面ミラー1409の焦点距離は、位置の関数としてf(x)と表すと、各々の入射位置において自由曲面ミラー1409に入射される光の焦点位置から、自由曲面ミラー1409までの距離が常にf(x)となるような自由曲面を形成している。これによって、自由曲面ミラー1409で反射された光は、どの偏向角度においても平行光となる。実施例5では、平行光を生成する例としたが、平行光にすることが必ずしも重要ではなく、出射光が偏向する範囲内で、各々の瞬間のビームの拡散角を一定に保つことが本質である。
 実施例5では、制御電圧源1404からの印加電圧として、マイナスのDCバイアス電圧を印加しているため、電気光学結晶1401内の電荷密度が、制御電圧源側の電極403付近で大きくなるため、自由曲面ミラー1409の曲率の変化と電極の関係を上記のようにしている。印加電圧として、プラスのDCバイアス電圧を印加する場合には、電気光学結晶1401内の接地側の電極1402付近の電荷密度が大きくなる。従って、自由曲面ミラー1409の曲率が電極1402側で小さくなるように配置すれば良い。
 図26は、実施例6にかかる波長掃引光源に好適な光偏向器の構成を示す。KTN、KLTNなどの電気光学結晶1501を備える光偏向器は、電気光学結晶1501の対向する2つの面のそれぞれに、電極1502、1503が形成されている。電気光学結晶1501に入射される入射光1505は、シリンドリカル凹レンズ1507を透過した後に、電気光学結晶1501に入射される。電気光学結晶501の形状、駆動条件は、実施例4と同じである。
 実施例6では、電気光学結晶1501の後段に、透過型の光学素子1509を用いることにより、光偏向器によって発現される偏向角に依存した集光効果を補正している。光学素子1509は、SF11ガラスを硝材とし、出射面において自由曲面が形成され、出射光の偏向角に応じて出射面の曲率が変化しており、出射光が偏向する範囲にわたって空間的に集光特性が分布している。
 実施例6では透過型の光学素子1509を用いたが、同様の集光効果を奏するように、素子の内部に屈折率分布を有する光学素子を用いても良い。この場合、光学素子の形状は、表面が平面となるので、光学素子の形状加工が容易になる。このような光学素子は、内部の屈折率分布を調整することにより、出射面の形状で集光特性を与える光学素子1509に比べ、空間的に分解能の高い集光特性を与えることができる。
 実施例5および実施例6によれば、自由曲面を有するミラーなどの光学素子を、電気光学結晶の出射側に配置することによって、光偏向器の出射光の偏向の範囲内で、ビームの拡散角を一定に保つことができる。この光偏向器を、回折格子などの分散素子と組み合わせて分光装置として用いると、分光帯域全域にわたって分解能が劣化することなく、一様な分解能の分光特性を得ることができる。
 図27は、実施例7にかかる波長掃引光源に好適な光偏向器の構成を示す図である。KTN、KLTNなどの電気光学結晶1601を備える光偏向器は、電気光学結晶1601の対向する2つの面のそれぞれに、電極1602、1603が形成されている。電気光学結晶1601の出射側にシリンドリカル凹レンズ1607を配置して、出射光1608を平行光としている。すなわち、シリンドリカル凹レンズ1607の焦点距離は、光偏向器による出射光の偏向範囲内で焦点距離が大きい偏向角(焦点位置が電気光学結晶1601の出射面から遠い電圧)の時のビームが、平行光になるように選択する。光偏向器によって発現されるビームの集光効果は、印加電圧と、印加電圧に対応する偏向角に依存している。従って、一部の偏向角においては、ビームの集光特性が良好に打ち消されているが、集光特性を打ち消しきれない偏向角のビームが存在している。
 そこで、シリンドリカル凹レンズ1607の出射側に、さらに空間位相変調器1609を配置した。空間位相変調器1609は、透過光または反射光の断面内において、各々場所の位相を独立に変調することができる光学素子である。電気光学結晶1601にて偏向された光は、空間位相変調器1609に入射される。空間位相変調器1609は、中央から周辺部に向かって位相を大きく遅らせるように位相変調量が連続的に制御されている。位相変調量は、偏向された出射光のうち、焦点位置が電気光学結晶1601の出射面に近い偏向角(電圧)のときの出射光が照射される面は大きく、焦点位置が遠い偏向角(電圧)のときの出射光が照射される面は小さくなっている。
 このようにして、シリンドリカル凹レンズ1607によって平行光とされた偏向角の出射光1610に対しては、空間位相変調器1609で位相変化を与えず、シリンドリカル凹レンズ1607によって集光効果が打ち消されなかった偏向角の出射光1606に対して集光効果を打ち消すような空間分布の位相変調を与える。空間位相変調器1609が挿入される前の焦点距離に応じて、各々の出射光の焦点距離を延伸することにより、出射光の偏向する範囲内全域にわたって、光の拡散角を一定に保つことができる。
 空間位相変調器として、空間的な位相変調分布を有するホログラフィック光学素子、回折光学素子を用いても良い。
 レンズ効果を有する光学素子にガウシアンビームを入射する場合、入射側のビームウェストの光学素子からの位置をs、出射側のビームウェストの位置をs''とすると、下記のSelfの式の関係が成り立つ。
Figure JPOXMLDOC01-appb-I000007

はレイリー長、fは光学素子の焦点距離である。この式により焦点距離を変化させた時のs、s''の関係を求めた結果を、図28に示す。
 実施例4~実施例7において、光偏向器(電気光学結晶)への入射光は、ほぼ平行光であると想定しており、入射側のビームウェストの位置は、光偏向器から離れた位置にある。このとき、図28の矢印Aで示すように、光偏向器のレンズパワー(=焦点距離の逆数)により、出射側のビームウェストの位置は、単調に増減する。しかしながら、光偏向器への入射側のビームウェストの位置は様々な場合があるので、図28の矢印Bで示すように、入射側ビームウェスト位置によっては単調に増減しない場合がある。このような場合においても光偏向器のレンズパワーを補正するためには、Selfの式の関係を考慮して、外部に配置する光学素子のレンズ効果を決定することが有効である。
 また、ビームウェスト位置だけでなく、ビームウェスト径も偏向角により変化するため、これを補正することが望まれる。Selfの式によれば、入射光と出射光の間のビームウェスト径の比率、すなわち倍率mは下式のように表される。
Figure JPOXMLDOC01-appb-I000008
 図29は、焦点距離を変化させた時のビームウェスト位置とビーム径倍率の関係を示す。図29の矢印Cで示すように、焦点距離の変化に応じ、ビームウェスト径も変化することがわかる。このように、偏向角に応じて光偏向器の集光特性(焦点距離)が変化し、出射光のビームウェスト位置、ビームウェスト径が上述したように偏向角により変化する。このとき、実施例4から実施例7に示した補償用の光学素子の特性も、偏向角に応じて変化させ、ビームウェスト位置とビーム径の変動を抑制することが望ましい。
 そこで、ビームウェスト位置が光偏向器の出射側にある場合についての補償方法を述べる。ビームウェスト位置が光偏向器の入射側にある場合も同様の方法で補償することができる。
 例えば、実施例4(図21)の光偏向器と同じ構成であって、補償光学素子としてプリズムを使用する場合について説明する。光がプリズム内部を透過するとき、光路長は、空気中よりも長くなるため、光偏向器と出射光のビームウェスト位置との間の距離を、プリズムにより調整することができる。すなわち、偏向角によってビームウェスト位置が変化するのに伴い、ビームウェスト位置が光偏向器から遠い場合にはプリズム内の光路長を長く設定し、ビームウェスト位置が光偏向器から近い場合にはプリズム内の光路長を短く設定すれば良い。
 ビームウェスト位置が偏向角により単調に増減する場合、ビームウェスト位置が遠い偏向角のときプリズムの幅が大きく、ビームウェスト位置が近い偏向角のときプリズムの幅が小さくなるように、プリズムの頂角側をビームウェスト位置が近い偏向角側に配置する。また、プリズムへの入射角度および出射角度によりビーム径の倍率が変わるので、ビームウェスト径を補償することもできる。さらに、プリズムの頂角およびプリズムへの入射角度を適当に決定することにより、ビームウェスト位置およびビームウェスト径を同時に補償することもできる。
 実施例4(図21)の構成に、さらにプリズムを追加する構成とすれば、光路長とビーム径の倍率とをある程度独立に設定することができるため、補償の精度を向上させることができる。
 次に、実施例5(図23)の光偏向器と同じ構成であって、補償光学素子として自由曲面ミラーを使用する場合について説明する。自由曲面ミラーの場合には、ミラーに入射するまでの光路長とレンズパワーとを制御することにより、出射光のビームウェスト位置とビームウェスト径を変化させることができる。ビームウェスト位置、ビームウェスト径のいずれかの変化を補償するように自由曲面ミラーの形状を決定しても良いし、両方の変化を抑制するように自由曲面ミラーの形状を決定しても良い。例えば、ビームウェスト位置が光偏向器から遠い偏向角に対して自由曲面ミラーの焦点距離を短くし、ビームウェスト位置が光偏向器から近い偏向角に対して自由曲面ミラーの焦点距離を長くすれば良い。
 同様にして、空間位相変調器、屈折率分布を有する媒体を用いる場合も、偏向角に応じて位相変調量、屈折率分布を変化させることにより、出射光のビームウェスト位置、ビームウェスト径を調整することができ、偏向角に依存した変化を抑えることができる。
 実施例5(図23)の構成に、さらに自由曲面ミラーを追加する構成とすれば、光路長とビームウェスト径をある程度独立に設定することができるため、補償の精度を向上させることができる。
 ビームウェスト位置が入射側にある場合も、同様にして、光偏向器の入射側に空間的に光学特性が分布している光学素子を配置して、ビームウェスト位置と光偏向器の距離に応じて補償を行えば良い。
 実施例8によれば、光偏向器の出射光の偏向の範囲内で、出射光のビームウェスト位置、ビームウェスト径の一方、または双方のバラツキを抑制することができる。また、入射光のビームウェスト位置、ビームウェスト径に応じて補償光学系を最適化することもできる。この光偏向器を、回折格子などの分散素子と組み合わせて分光装置として用いると、分光帯域全域にわたって分解能が劣化することなく、一様な分解能の分光特性を得ることができる。
 尚、実施例4から実施例8においては、KTNまたはKLTNを用いた光偏向器について説明したが、これら以外の電気光学効果を有する結晶を用いた光偏向器であっても良い。実施例4から実施例8の光偏向器は、第1の実施形態および第2の実施形態の波長掃引光源にそれぞれ好適に利用できる。
 以上、様々な実施形態および実施例とともに詳細に述べたように、本発明によって、広い掃引波長帯域の全体に渡って、コヒーレンス長が長い状態を維持することができる波長掃引光源を提供できる。SS-OCT装置のイメージ品質を向上させることができる。
 本発明は、一般的にイメージング装置に利用することができる。特に、SS-OCT装置の波長掃引光源として利用できる。

Claims (20)

  1.  発振波長が連続的に変化する光を出力する波長掃引光源において、
      光を増幅するゲイン媒体、
      発振波長を変更する手段、および
      前記ゲイン媒体の光利得を変更する手段
     を含むレーザ共振器と、
     波長制御信号によって前記発振波長を変更する手段を制御して、前記発振波長を掃引する波長制御手段と、
     波長掃引の全波長帯域にわたって、前記レーザ共振器内の光利得が、レーザ発振閾値の状態の光利得よりもわずかに高い光利得値となるように、前記光利得を変更する手段を制御する光利得制御手段と
     を備えたことを特徴とする波長掃引光源。
  2.  前記光利得制御手段は、
     前記波長制御手段から得られた掃引開始からの複数の時間と、前記複数の時間の各々において波長掃引光のコヒーレンス長が最長となる光利得とを組み合わせた光利得制御データに基づいて、光利得制御信号によって前記光利得を変更する手段を制御するよう構成されたことを特徴とする請求項1に記載の波長掃引光源。
  3.  前記光利得制御手段に関連付けられ、前記光利得制御データを記憶したメモリをさらに備え、
     前記光利得制御データは、
      前記光利得を変更する手段が、波長掃引期間に渡って、一定値の前記光利得制御信号によって前記光利得を変更する手段を設定光利得値に制御し、
      複数の異なる値の前記設定光利得値の各々に対して、掃引開始からの前記複数の時間の各々におけるコヒーレンス長をそれぞれ測定し、前記設定光利得値、前記複数の時間の内の1つの時間および対応する前記測定されたコヒーレンス長からなるデータ組を生成し、
      前記複数の時間の内の前記1つの時間において、前記測定されたコヒーレンス長の中から最大コヒーレンス長を与える設定光利得値を選択し、
      前記1つの時間の各々と前記選択された設定光利得値から生成されること
     を特徴とする請求項2に記載の波長掃引光源。
  4.  前記光利得制御データは、
     前記設定光利得値に設定されたレーザ共振器から出力された、前記波長掃引期間の1回の波長掃引から得られる波長掃引光に対して、異なる干渉条件で当該波長掃引光から前記干渉条件の各々に対応する複数の干渉信号が生成され、前記干渉条件の各々に対応する前記複数の干渉信号に基づいて求められた前記コヒーレンス長を含むことを特徴とする請求項3に記載の波長掃引光源。
  5.  前記ゲイン媒体の光利得を変更する手段は、半導体光増幅器(SOA)であって、前記光利得制御信号は前記SOAを駆動するSOA電流であり、光利得制御データは前記複数の時間の各々に対して選択されたSOA電流を表すSOA電流波形データであることを特徴とする請求項2乃至4いずれかに記載の波長掃引光源。
  6.  前記発振波長を変更する手段は、
      回折格子と、
      前記回折格子からの発振光が入射し、当該入射した前記発振光を垂直に反射するミラーと、
      前記回折格子への前記発振光の入射角を変更する光偏向器と
    から構成されることを特徴とする請求項2乃至4いずれかに記載の波長掃引光源。
  7.  前記光利得制御手段は、
     前記レーザ共振器からの出力光レベルを検出し、前記検出された出力光レベルに基づいて、前記出力光レベルが、所定の出力光レベル値になるような光利得制御信号を生成し、当該光利得制御信号によって前記光利得を変更する手段を制御するよう構成されたことを特徴とする請求項1に記載の波長掃引光源。
  8.  前記ゲイン媒体の光利得を変更する手段は、半導体光増幅器(SOA)であって、前記光利得制御信号は前記SOAを駆動するSOA電流であることを特徴とする請求項7に記載の波長掃引光源。
  9.  前記光利得制御手段は、波長掃引に連動して時間的に離散的な制御を実施し、前記出力光レベルと前記所定の出力光レベルとの差分に対して、前記光利得の変更量を漸次半減する2分法を用いて制御されることを特徴とする請求項7に記載の波長掃引光源。
  10.  前記光利得制御手段は、波長掃引に連動して時間的に離散的な制御を実施し、前記出力光レベルと前記所定の出力光レベルとの差分に対して、前記光利得の変更量を、当該差分、前記差分の積分、および前記差分の微分の要素によって行うPID法を用いて制御されることを特徴とする請求項7に記載の波長掃引光源。
  11.  前記発振波長を変更する手段は、
      回折格子と、
      前記回折格子への発振光の入射角を変更する光偏向器と
    から構成されることを特徴とする請求項7乃至10いずれかに記載の波長掃引光源。
  12.  前記光偏向器は、
     電気光学結晶の対向する面に少なくとも2つの電極が形成された光偏向器であって、前記少なくとも2つの電極に電圧を印加すると、電気光学効果により前記電気光学結晶内に一様でない屈折率分布が生じることにより、前記電圧により形成された電界に直交する光の進路を屈曲させる光偏向器と、
     前記電気光学結晶の入射側と出射側の少なくとも一方に配置された光学素子であって、前記光が偏向する範囲にわたって空間的に集光特性が分布している光学素子と
     を備えたことを特徴とする請求項6または11に記載の波長掃引光源。
  13.  前記光偏向器は、偏向角により集光特性が変化し、
     前記光学素子は、前記光が偏向する範囲にわたって出射側のビームウェスト位置とビームウェスト径の少なくとも一方の変化を補償するように空間的に光学特性が分布していることを特徴とする請求項12に記載の波長掃引光源。
  14.  前記光学素子は、プリズムであり、前記電気光学結晶からの出射光の焦点位置が前記電気光学結晶の出射面から遠い偏向角の側に、前記プリズムの頂角が配置されていることを特徴とする請求項13に記載の波長掃引光源。
  15.  前記光学素子は、自由曲面ミラーであり、前記電気光学結晶からの出射光の焦点位置が前記電気光学結晶の出射面から遠い偏向角の側に、前記自由曲面ミラーの曲率半径が大きくなるように配置されていることを特徴とする請求項13に記載の波長掃引光源。
  16.  前記光学素子は、透過型の光学素子であり、前記電気光学結晶からの出射光の焦点位置が前記電気光学結晶の出射面から遠い偏向角の側に、前記透過型の光学素子の出射面の曲率半径が大きくなるように配置されていることを特徴とする請求項13に記載の波長掃引光源。
  17.  前記光学素子は、空間位相変調器であり、前記電気光学結晶からの出射光の焦点位置が前記電気光学結晶の出射面から遠い偏向角の側に、前記空間位相変調器の位相変調量が大きくなるように配置されていることを特徴とする請求項13に記載の波長掃引光源。
  18.  前記電気光学結晶の入射側と出射側の少なくとも一方に配置された第2の光学素子であって、前記光が偏向する範囲にわたって出射側のビームウェスト位置およびビームウェスト径の変化を補償するように空間的に光学特性が分布している第2の光学素子をさらに備えたことを特徴とする請求項12乃至17のいずれかに記載の波長掃引光源。
  19.  前記光偏向器は、タンタル酸ニオブ酸カリウム(KTa1-xNbx3(0≦x≦1))結晶、リチウムを添加したK1-yLiyTa1-xNbx3(0≦x≦1、0≦y≦1)結晶のいずれかであることを特徴とする請求項12乃至18のいずれかに記載の波長掃引光源。
  20.  発振波長が連続的に変化する光を出力する波長掃引光源において、
      光を増幅するゲイン媒体、
      発振波長を変更する手段、および
      前記ゲイン媒体の光利得を変更する手段
     を含むレーザ共振器と、
     波長制御信号によって前記発振波長を変更する手段を制御して、前記発振波長を掃引する波長制御手段と、
     前記波長制御手段から得られた掃引開始からの複数の時間と、前記複数の時間の各々において前記波長掃引光のコヒーレンス長が最長となる光利得とを組み合わせた光利得制御データに基づいて、光利得制御信号によって前記光利得を変更する手段と
     を備えたことを特徴とする波長掃引光源。
PCT/JP2017/043446 2016-12-09 2017-12-04 波長掃引光源、波長掃引光源のための駆動データ作成方法および光偏向器 WO2018105549A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
EP17878420.3A EP3553903B1 (en) 2016-12-09 2017-12-04 Wavelength-swept light source, method for generating drive data for wavelength-swept light source, and optical deflector
US16/466,608 US11165219B2 (en) 2016-12-09 2017-12-04 Swept light source and drive data generation method and optical deflector for swept light source
EP23153071.8A EP4191808A1 (en) 2016-12-09 2017-12-04 Swept light source and drive data generation method and optical deflector for swept light source
JP2018554986A JP6774504B2 (ja) 2016-12-09 2017-12-04 波長掃引光源、波長掃引光源のための駆動データ作成方法および光偏向器
US17/380,309 US11721948B2 (en) 2016-12-09 2021-07-20 Swept light source and drive data generation method and optical deflector for swept light source

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2016-239354 2016-12-09
JP2016239354 2016-12-09

Related Child Applications (2)

Application Number Title Priority Date Filing Date
US16/466,608 A-371-Of-International US11165219B2 (en) 2016-12-09 2017-12-04 Swept light source and drive data generation method and optical deflector for swept light source
US17/380,309 Division US11721948B2 (en) 2016-12-09 2021-07-20 Swept light source and drive data generation method and optical deflector for swept light source

Publications (1)

Publication Number Publication Date
WO2018105549A1 true WO2018105549A1 (ja) 2018-06-14

Family

ID=62492324

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2017/043446 WO2018105549A1 (ja) 2016-12-09 2017-12-04 波長掃引光源、波長掃引光源のための駆動データ作成方法および光偏向器

Country Status (4)

Country Link
US (2) US11165219B2 (ja)
EP (2) EP4191808A1 (ja)
JP (1) JP6774504B2 (ja)
WO (1) WO2018105549A1 (ja)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10610096B2 (en) 2016-12-21 2020-04-07 Acucela Inc. Miniaturized mobile, low cost optical coherence tomography system for home based ophthalmic applications
US11357401B2 (en) 2018-06-20 2022-06-14 Acucela Inc. Miniaturized mobile, low cost optical coherence tomography system for home based ophthalmic applications
US11393094B2 (en) 2020-09-11 2022-07-19 Acucela Inc. Artificial intelligence for evaluation of optical coherence tomography images
US11497396B2 (en) 2021-03-24 2022-11-15 Acucela Inc. Axial length measurement monitor
US11684254B2 (en) 2020-08-04 2023-06-27 Acucela Inc. Scan pattern and signal processing for optical coherence tomography
US11730363B2 (en) 2019-12-26 2023-08-22 Acucela Inc. Optical coherence tomography patient alignment system for home based ophthalmic applications
WO2023218782A1 (ja) * 2022-05-12 2023-11-16 浜松ホトニクス株式会社 光源装置及び制御方法
US11911105B2 (en) 2020-09-30 2024-02-27 Acucela Inc. Myopia prediction, diagnosis, planning, and monitoring device
US11974807B2 (en) 2020-08-14 2024-05-07 Acucela Inc. System and method for optical coherence tomography a-scan decurving

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4191808A1 (en) 2016-12-09 2023-06-07 Nippon Telegraph And Telephone Corporation Swept light source and drive data generation method and optical deflector for swept light source
JP2022138828A (ja) * 2021-03-11 2022-09-26 株式会社トプコン スペクトル測定方法及びスペクトル測定装置
CN117119643B (zh) * 2023-10-20 2024-04-02 光越科技(深圳)有限公司 一种扫频光源及扫频光学相干层析成像系统

Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59207756A (ja) * 1983-05-10 1984-11-24 Yokogawa Hewlett Packard Ltd 光信号発生器
JP2003023198A (ja) * 2001-07-09 2003-01-24 Gigaphoton Inc レーザ装置用波長制御装置及び制御方法
JP2005347668A (ja) 2004-06-07 2005-12-15 Sun Tec Kk 波長走査型ファイバレーザ光源
JP2007109923A (ja) * 2005-10-14 2007-04-26 Olympus Corp 光検出装置およびそれを用いた光通信システム
JP2007115900A (ja) * 2005-10-20 2007-05-10 Nec Corp 波長可変光源、波長可変光源モジュール、および波長可変光源の駆動方法
US20080037608A1 (en) * 2005-09-06 2008-02-14 Yan Zhou Light source for swept source optical coherence tomography based on cascaded distributed feedback lasers with engineered band gaps
JP2008529068A (ja) * 2005-01-24 2008-07-31 ソルラブス、 インコーポレイテッド 高速に波長スキャンする小型マルチモードレーザ
JP2009198593A (ja) * 2008-02-19 2009-09-03 Nippon Telegr & Teleph Corp <Ntt> 可変分散補償器
JP2010062426A (ja) * 2008-09-05 2010-03-18 Sun Tec Kk 波長走査型レーザ光源
JP2010175271A (ja) * 2009-01-27 2010-08-12 Sun Tec Kk 光断層画像表示システム
JP2011186218A (ja) 2010-03-09 2011-09-22 Nippon Telegr & Teleph Corp <Ntt> 光偏向器の駆動方法
JP2012151419A (ja) * 2011-01-21 2012-08-09 Nippon Telegr & Teleph Corp <Ntt> 波長可変光源
JP2014042010A (ja) * 2012-07-25 2014-03-06 Canon Inc 波長掃引光源の駆動方法
JP2014126599A (ja) * 2012-12-25 2014-07-07 Nippon Telegr & Teleph Corp <Ntt> ビームスキャナおよび偏向光源
JP2015142111A (ja) 2014-01-30 2015-08-03 日本電信電話株式会社 波長掃引光源
US20150373431A1 (en) * 2014-06-23 2015-12-24 Canon Kabushiki Kaisha Split Spectrum Spectral Domain Optical Coherence Tomography
JP2016032608A (ja) * 2014-07-31 2016-03-10 株式会社ニデック 光コヒーレンストモグラフィー装置、および眼底画像処理プログラム

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004188828A (ja) 2002-12-12 2004-07-08 Sony Corp 可逆性多色記録媒体の記録装置
US20060193354A1 (en) * 2003-02-05 2006-08-31 Yehuda Rosenblatt External Cavity Tunable Laser and Control
JP4770693B2 (ja) 2006-10-17 2011-09-14 セイコーエプソン株式会社 走査型光学装置
EP2662683B1 (en) * 2011-01-05 2016-05-04 Nippon Telegraph And Telephone Corporation Wavelength swept light source
JP5565910B2 (ja) 2011-01-21 2014-08-06 日本電信電話株式会社 光偏向器
US8994954B2 (en) * 2012-12-28 2015-03-31 Axsun Technologies, Inc. System and method for stabilizing mode locked swept laser for OCT medical imaging
JP6040103B2 (ja) 2013-06-06 2016-12-07 浜松ホトニクス株式会社 補償光学システムの対応関係特定方法、補償光学システム、および補償光学システム用プログラム
JP6006693B2 (ja) 2013-08-06 2016-10-12 日本電信電話株式会社 光偏向装置
KR20150054542A (ko) * 2013-11-12 2015-05-20 삼성전자주식회사 광대역 파장 가변광원 시스템 및 이를 적용한 장치
EP4191808A1 (en) 2016-12-09 2023-06-07 Nippon Telegraph And Telephone Corporation Swept light source and drive data generation method and optical deflector for swept light source

Patent Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59207756A (ja) * 1983-05-10 1984-11-24 Yokogawa Hewlett Packard Ltd 光信号発生器
JP2003023198A (ja) * 2001-07-09 2003-01-24 Gigaphoton Inc レーザ装置用波長制御装置及び制御方法
JP2005347668A (ja) 2004-06-07 2005-12-15 Sun Tec Kk 波長走査型ファイバレーザ光源
JP2008529068A (ja) * 2005-01-24 2008-07-31 ソルラブス、 インコーポレイテッド 高速に波長スキャンする小型マルチモードレーザ
US20080037608A1 (en) * 2005-09-06 2008-02-14 Yan Zhou Light source for swept source optical coherence tomography based on cascaded distributed feedback lasers with engineered band gaps
JP2007109923A (ja) * 2005-10-14 2007-04-26 Olympus Corp 光検出装置およびそれを用いた光通信システム
JP2007115900A (ja) * 2005-10-20 2007-05-10 Nec Corp 波長可変光源、波長可変光源モジュール、および波長可変光源の駆動方法
JP2009198593A (ja) * 2008-02-19 2009-09-03 Nippon Telegr & Teleph Corp <Ntt> 可変分散補償器
JP2010062426A (ja) * 2008-09-05 2010-03-18 Sun Tec Kk 波長走査型レーザ光源
JP2010175271A (ja) * 2009-01-27 2010-08-12 Sun Tec Kk 光断層画像表示システム
JP2011186218A (ja) 2010-03-09 2011-09-22 Nippon Telegr & Teleph Corp <Ntt> 光偏向器の駆動方法
JP2012151419A (ja) * 2011-01-21 2012-08-09 Nippon Telegr & Teleph Corp <Ntt> 波長可変光源
JP2014042010A (ja) * 2012-07-25 2014-03-06 Canon Inc 波長掃引光源の駆動方法
JP2014126599A (ja) * 2012-12-25 2014-07-07 Nippon Telegr & Teleph Corp <Ntt> ビームスキャナおよび偏向光源
JP2015142111A (ja) 2014-01-30 2015-08-03 日本電信電話株式会社 波長掃引光源
US20150373431A1 (en) * 2014-06-23 2015-12-24 Canon Kabushiki Kaisha Split Spectrum Spectral Domain Optical Coherence Tomography
JP2016032608A (ja) * 2014-07-31 2016-03-10 株式会社ニデック 光コヒーレンストモグラフィー装置、および眼底画像処理プログラム

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
See also references of EP3553903A4
YOSHIAKI YASUNO ET AL.: "Three-dimensional and high-speed swept-source optical coherence tomography for in vivo investigation of human anterior eye segments", OPTICS EXPRESS, vol. 13, no. 26, 2005, pages 10652 - 10664, XP055506441, DOI: doi:10.1364/OPEX.13.010652

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10952607B2 (en) 2016-12-21 2021-03-23 Acucela Inc. Miniaturized mobile, low cost optical coherence tomography system for home based ophthalmic applications
US10610096B2 (en) 2016-12-21 2020-04-07 Acucela Inc. Miniaturized mobile, low cost optical coherence tomography system for home based ophthalmic applications
US11890053B2 (en) 2016-12-21 2024-02-06 Acucela Inc. Miniaturized mobile, low cost optical coherence tomography system for home based ophthalmic applications
US11627874B2 (en) 2016-12-21 2023-04-18 Acucela Inc. Miniaturized mobile, low cost optical coherence tomography system for home based ophthalmic applications
US11357401B2 (en) 2018-06-20 2022-06-14 Acucela Inc. Miniaturized mobile, low cost optical coherence tomography system for home based ophthalmic applications
US11896308B2 (en) 2018-06-20 2024-02-13 Acucela Inc. Miniaturized mobile, low cost optical coherence tomography system for home based ophthalmic applications
US11576572B2 (en) 2018-06-20 2023-02-14 Acucela Inc. Miniaturized mobile, low cost optical coherence tomography system for home based ophthalmic applications
US11730363B2 (en) 2019-12-26 2023-08-22 Acucela Inc. Optical coherence tomography patient alignment system for home based ophthalmic applications
US11684254B2 (en) 2020-08-04 2023-06-27 Acucela Inc. Scan pattern and signal processing for optical coherence tomography
US11974807B2 (en) 2020-08-14 2024-05-07 Acucela Inc. System and method for optical coherence tomography a-scan decurving
US11393094B2 (en) 2020-09-11 2022-07-19 Acucela Inc. Artificial intelligence for evaluation of optical coherence tomography images
US11798164B2 (en) 2020-09-11 2023-10-24 Acucela Inc. Artificial intelligence for evaluation of optical coherence tomography images
US11620749B2 (en) 2020-09-11 2023-04-04 Acucela Inc. Artificial intelligence for evaluation of optical coherence tomography images
US11911105B2 (en) 2020-09-30 2024-02-27 Acucela Inc. Myopia prediction, diagnosis, planning, and monitoring device
US11497396B2 (en) 2021-03-24 2022-11-15 Acucela Inc. Axial length measurement monitor
US11779206B2 (en) 2021-03-24 2023-10-10 Acucela Inc. Axial length measurement monitor
WO2023218782A1 (ja) * 2022-05-12 2023-11-16 浜松ホトニクス株式会社 光源装置及び制御方法

Also Published As

Publication number Publication date
JP6774504B2 (ja) 2020-10-28
US11721948B2 (en) 2023-08-08
JPWO2018105549A1 (ja) 2019-04-11
US20210408755A1 (en) 2021-12-30
US20200076154A1 (en) 2020-03-05
US11165219B2 (en) 2021-11-02
EP4191808A1 (en) 2023-06-07
EP3553903A4 (en) 2020-08-12
EP3553903B1 (en) 2023-04-12
EP3553903A1 (en) 2019-10-16

Similar Documents

Publication Publication Date Title
WO2018105549A1 (ja) 波長掃引光源、波長掃引光源のための駆動データ作成方法および光偏向器
JP5913135B2 (ja) 波長掃引光源
JP5750138B2 (ja) 波長可変光源を利用した光干渉断層撮影法の方法及び装置
JP6389455B2 (ja) 最適化されたkクロックを有する多速度oct掃引光源
US20120026505A1 (en) Swept Frequency Laser for FD OCT with Intracavity Element and Method of Operation
JP2004037165A (ja) 干渉計装置
JP2013546189A (ja) 大きなコム間隔を有する周波数コム源
KR101642952B1 (ko) 두 파장 레이저 간섭계를 이용한 길이 측정 장치
US10205296B2 (en) Swept light source and method for controlling the same
JP2015142111A (ja) 波長掃引光源
CN104180972A (zh) 一种测量光学腔线宽的装置和方法
JP6180961B2 (ja) 干渉計の較正方法及びこの較正方法を用いた干渉計
Naber et al. Optical techniques for Rydberg physics in lattice geometries: A technical guide
US10161738B2 (en) OCT swept laser with cavity length compensation
JP5296814B2 (ja) 波長掃引光源
JP5357195B2 (ja) 波長可変光源
JP6047059B2 (ja) データ処理装置およびリサンプリング方法
JP2018022041A (ja) 波長掃引光源
JP6621026B2 (ja) 波長掃引光源の評価方法
JP5730076B2 (ja) 光源装置、検査装置及び光干渉断層撮像装置
JP2015090928A (ja) 波長掃引光源
JP6053017B2 (ja) 波長掃引光源
JP2019054223A (ja) 波長掃引光源

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17878420

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2018554986

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2017878420

Country of ref document: EP

Effective date: 20190709