1. Introduction
Lidar is an important technique for the measurement of range, velocity and vibration, as well as atmospheric parameters [
1J. M. Vaughan, “Coherent laser spectroscopy and Doppler lidar sensing in the atmosphere,” Phys. Scr. T T78(1), 73–81 (1998). [CrossRef]
–
4X. Mao, D. Inoue, S. Kato, and M. Kagami, “Amplitude-modulated laser radar for range and speed measurement in car applications,” IEEE Trans. Intell. Transp. Syst. 13(1), 408–413 (2012). [CrossRef]
]. For velocity detection, one of the commonly used lidar schemes is pulsed Doppler lidar [
3M. Harris, R. I. Yong, F. Köpp, A. Dolfi, and J.-P. Cariou, “Wake vortex detection and monitoring,” Aerosp. Sci. Technol. 6(5), 325–331 (2002). [CrossRef]
,
5A. D. Bouteyre, G. Canat, M. Valla, B. Augère, C. Besson, D. Goular, L. Lombard, J.-P. Cariou, A. Durecu, D. Fleury, L. Bricteux, S. Brousmiche, S. Lugan, and B. Macq, “Pulsed 1.5ìm lidar for axial aircraft wake vortex detection based on high brightness large-core fiber amplifier,” IEEE J. Sel. Top. Quantum Electron. 15, 441–450 (2009).
], where in the transmitter, optical pulses are intensity modulated at a constant frequency
fCW before launched out through a telescope. When the target is moving, the reflected signal acquires a Doppler frequency shift
, where
v is the radial velocity of target and λ is the laser wavelength. This reflected optical signal is detected by an optical receiver, and the frequency
carried on the signal can be obtained through a Fourier process so that the target velocity can be determined. The value of the Doppler frequency shift
fd can be either positive or negative depending on the direction of the velocity (moving towards or away from the observer). Theoretically,
fCW determines the maximum unambiguous velocity that could be measured without spectral aliasing of the received signal. Long optical pulses are usually desired to ensure the accuracy of velocity measurements as the resolution of velocity is determined by
, where
T is the pulse duration. On the other hand, it is well known that short optical pulses are required to achieve high range resolution if the target distance needs to be measured. This is because the target range information is obtained by measuring the roundtrip time of the optical pulse between the lidar and the target [
4X. Mao, D. Inoue, S. Kato, and M. Kagami, “Amplitude-modulated laser radar for range and speed measurement in car applications,” IEEE Trans. Intell. Transp. Syst. 13(1), 408–413 (2012). [CrossRef]
,
5A. D. Bouteyre, G. Canat, M. Valla, B. Augère, C. Besson, D. Goular, L. Lombard, J.-P. Cariou, A. Durecu, D. Fleury, L. Bricteux, S. Brousmiche, S. Lugan, and B. Macq, “Pulsed 1.5ìm lidar for axial aircraft wake vortex detection based on high brightness large-core fiber amplifier,” IEEE J. Sel. Top. Quantum Electron. 15, 441–450 (2009).
]. Clearly, in a pulsed Doppler lidar system, a tradeoff has to be made between the velocity resolution and the range resolution.
In comparison to short pulsed lidar systems, FMCW lidar uses frequency chirped long optical pulses at much lower peak power and is able to provide large dynamic range while maintaining excellent resolution [
6C. J. Karlsson and F. Å. A. Olsson, “Linearization of the frequency sweep of a frequency-modulated continuous-wave semiconductor laser radar and the resulting ranging performance,” Appl. Opt. 38(15), 3376–3386 (1999). [CrossRef] [PubMed]
–
8M.-C. Amann, T. Bosch, M. Lescure, R. Myllylä, and M. Rioux, “Laser ranging: a critical reveiw of usual techniques for distance measurement,” Opt. Eng. 40(1), 10–19 (2001). [CrossRef]
]. For range measurements, saw-tooth or triangular waveforms are often used as the chirping signal in which the modulation frequency linearly shifts between the upper and the lower boundaries
f1 and
f2, as shown in
Fig. 1
[
8M.-C. Amann, T. Bosch, M. Lescure, R. Myllylä, and M. Rioux, “Laser ranging: a critical reveiw of usual techniques for distance measurement,” Opt. Eng. 40(1), 10–19 (2001). [CrossRef]
,
9D. Pierrottet, F. Amzajerdian, L. Petway, B. Barnes, G. Lockard, and M. Rubio, “Linear FMCW laser radar for precision range and vector velocity measurements,” Proc. Meter. Res. Soc. Symp. 1076 (2008).
]. When a saw-tooth waveform is used, on a stationary target, the target range
R can be calculated through the measurement of the frequency difference between the transmitted and the returned modulation signal [
7P. Adany, C. Allen, and R. Hui, “Chirped lidar using simplified homodyne detection,” J. Lightwave Technol. 27(16), 3351–3357 (2009). [CrossRef]
]
often referred to as the de-chirped frequency. In
Eq. (1),
T is the optical pulse duration,
is the chirping bandwidth, and c is the speed of light. The range resolution
only depends on the chirping bandwidth. For a moving target, an additional Doppler frequency shift will be superimposed on
fR . This shift cannot be easily isolated so that a saw-tooth waveform is not suitable for the range and velocity measurement of moving targets.
Fig. 1 Illustration of linear FM chirped signal for FMCW lidar. T’ is the period and T is the pulse width.
The triangular chirping waveform illustrated in
Fig. 1(b) contains both up-ramp and down-ramp segments., A lidar system using this chirp can detect range and velocity independently by calculating the two related frequencies
where,
and
are de-chirped frequencies measured during the up-ramp and the down-ramp sections, respectively. the velocity direction cannot be determined from a simple frequency analysis. Instead, de-chirping for up and down ramp sections must be performed separately within their respective temporal slots [
9D. Pierrottet, F. Amzajerdian, L. Petway, B. Barnes, G. Lockard, and M. Rubio, “Linear FMCW laser radar for precision range and vector velocity measurements,” Proc. Meter. Res. Soc. Symp. 1076 (2008).
].
In general, fine range resolution can be obtained with FMCW lidar by increasing chirping bandwidth. Whereas, a simple continuous wave at a single frequency is preferred in Doppler lidar for velocity measurements, which provides the most accurate results by avoiding possible complications of wideband modulation and the difficulty of maintaining the linearity of chirping. For simultaneous range and velocity measurements, a composite lidar waveform consisting of a wideband chirp and a single frequency tone may help optimize the overall performance.
In the lidar receiver [
3M. Harris, R. I. Yong, F. Köpp, A. Dolfi, and J.-P. Cariou, “Wake vortex detection and monitoring,” Aerosp. Sci. Technol. 6(5), 325–331 (2002). [CrossRef]
,
5A. D. Bouteyre, G. Canat, M. Valla, B. Augère, C. Besson, D. Goular, L. Lombard, J.-P. Cariou, A. Durecu, D. Fleury, L. Bricteux, S. Brousmiche, S. Lugan, and B. Macq, “Pulsed 1.5ìm lidar for axial aircraft wake vortex detection based on high brightness large-core fiber amplifier,” IEEE J. Sel. Top. Quantum Electron. 15, 441–450 (2009).
], coherent heterodyne detection is the most common choice because it has much higher receiver sensitivity compared with direct detection. It is also less susceptible to the phase noise of the optical carrier and the associated signal fading compared to homodyne detection. Heterodyne detection operates at an intermediate frequency
fIF, which is the frequency difference between the transmitted optical signal and the optical local oscillator (LO). A common method of creating this frequency shift in optics is to use an acousto-optic modulator (AOM). Practical AOMs can shift the optical frequency by as much as 1GHz. To avoid spectral aliasing using heterodyne detection, the maximum modulation bandwidth of the lidar system has to be lower than
fIF. This limits the range resolution. This limitation can be overcome by using an electro-optic I/Q modulator which is able to provide carrier-suppressed single-sideband modulation to generate a much higher
fIF for coherent heterodyne lidar system and achieve unprecedented range resolution [
10S. Gao and R. Hui, “Frequency-modulated continuous-wave lidar using I/Q modulator for simplified heterodyne detection,” Opt. Lett. 37(11), 2022–2024 (2012). [CrossRef] [PubMed]
].
Although using an I/Q modulator to generate multi-gigahertz intermediate frequency
fIF for heterodyne lidar helps improving range resolution, it does not make full use of the capability of an I/Q modulator. In fact, an I/Q modulator allows independent RF signals to be loaded onto the upper and the lower optical sidebands [
11Y. Zhang, M. O’Sullivan, and R. Hui, “Digital subcarrier multiplexing for flexible spectral allocation in optical transport network,” Opt. Express 19(22), 21880–21889 (2011). [CrossRef] [PubMed]
]. Generating high frequency
fIF in [
10S. Gao and R. Hui, “Frequency-modulated continuous-wave lidar using I/Q modulator for simplified heterodyne detection,” Opt. Lett. 37(11), 2022–2024 (2012). [CrossRef] [PubMed]
] only addresses one of the two optical sidebands. The other sideband is not used so that only range information was measured in a fiber delay line in reference [
10S. Gao and R. Hui, “Frequency-modulated continuous-wave lidar using I/Q modulator for simplified heterodyne detection,” Opt. Lett. 37(11), 2022–2024 (2012). [CrossRef] [PubMed]
].
In this paper we demonstrate that by applying a wide band chirping on the upper optical sideband and a constant frequency modulation on the lower optical sideband, it is possible to independently measure the target range and the velocity at the same time. A fiber-optic 90° hybrid was used in the coherent receiver to recover the complex optical spectrum composed of the in-phase and quadrature components. The resolution of velocity measurement is no longer dependent on the chirping bandwidth so that the system is more flexible and practical. The direction of velocity can be directly read from the spectrum of the received signal without the need of temporally truncated spectral analysis used in a conventional lidar system based on triangular chirping waveform. Thanks to the wide modulation bandwidth of the electro-optic I/Q modulator this lidar system allows wide chirping bandwidth on the upper sideband to ensure the fine range resolution and high enough modulation frequency on the lower sideband to allow the measurement of large unambiguous vector velocity, while maintaining a relatively simple optical system configuration. Because de-chirping in the receiver requires low-pass filtering [
7P. Adany, C. Allen, and R. Hui, “Chirped lidar using simplified homodyne detection,” J. Lightwave Technol. 27(16), 3351–3357 (2009). [CrossRef]
] after mixing, wide chirping bandwidth does not necessarily degrade lidar detection sensitivity.
2. System description
Figure 2
shows the block diagram of the proposed lidar system. The optical signal and the LO are split from the same CW laser source. The LO branch is directly connected to the receiver and the transmitted signal is modulated through an electro-optic I/Q modulator. The modulating RF signals from a digital arbitrary waveform generator (AWG) are fed to the I and the Q channels of the modulator independently. Optical carrier suppression and sideband rejection ratio are optimized by properly adjusting the DC bias of the modulator. As a result, a constant frequency modulation at
fCW is loaded onto the lower optical sideband and a wideband saw-tooth chirping signal is loaded onto the upper optical sideband. The modulated optical signal is amplified by an erbium doped fiber amplifier (EDFA), and a polarization controller (PC) behind the EDFA adjusts the state of polarization of the optical signal to match the principle axis of the 2x2 fiber-optic polarization beam splitter (PBS) connected to the telescope. A quarter wave plate is placed in front of the telescope to create a 90° polarization rotation for each roundtrip. Thus the optical signal reflected from the target and enters the telescope will pass through the PBS with minimum attenuation before it reaches the receiver through the adjacent fiber port. In the receiver, the optical signal is combined with the LO through a 2x2 fiber-optic 90° hybrid which allows the separation of the inphase and the quadrature components in the optical signal with two photodetectors. After photodetection, electrical signals are sampled by a pair of high speed analog-to-digital converters (ADC) before they are analyzed through digital signal processing.
Fig. 2 Block diagram of the proposed lidar system.
To create a wideband chirp waveform which linearly sweeps from
f1 to
f2 on the upper optical sideband, and a constant frequency modulation at
fCW on the lower optical sideband, the two voltage waveforms used to drive I and Q arms of the I/Q modulator are [
10S. Gao and R. Hui, “Frequency-modulated continuous-wave lidar using I/Q modulator for simplified heterodyne detection,” Opt. Lett. 37(11), 2022–2024 (2012). [CrossRef] [PubMed]
],
where, V
D is the amplitude of the driving voltage signal. When both intensity modulators (I and Q branches) inside the I/Q modulator are biased at null point, and the phase shifter is biased at the quadrature point, the optical field at the output of the modulator is [
11Y. Zhang, M. O’Sullivan, and R. Hui, “Digital subcarrier multiplexing for flexible spectral allocation in optical transport network,” Opt. Express 19(22), 21880–21889 (2011). [CrossRef] [PubMed]
]
where
f0 is the optical frequency of the CW laser, E
s is the constant input optical field and
Vπ is the voltage required for the transfer function to change from the minimum to the maximum for each intensity modulator. For small-signal modulation,
,
Eq. (6) can be linearized to
This is a carrier suppressed optical signal with independent modulating signals carried on the upper and the lower sidebands. Since the optical LO is not modulated, its optical field can be expressed as,
where,
ALO is the field amplitude. The optical field reflected from the target
Esig contains a time delay Δt and a Doppler frequency shift
fd compared to the launched optical field given by
Eq. (7), so that,
Where,
Asig is the amplitude of the reflected optical field. Sending
ELO and
Esig into the two input ports of the 90 degree optical hybrid, the optical fields
E1 and
E2 exiting the two output ports can be described by
Photocurrents of the two photodetectors after the 90° optical hybrid are,
where
is the responsivity of the photodetector, and terms due to direct detection have been neglected for simplicity. Because these two photocurrent signals represent the real and the imaginary parts of the optical signal, the complex optical signal can be recovered in the electrical domain as,
The complex optical spectrum can be obtained through a fast Fourier transform (FFT) on this composite electrical signal
Ifinal, in which the positive and the negative sidebands represented by the 1st and the 2nd terms on the RHS of
Eq. (13) are not redundant. Doppler frequency shift
fd can be obtained directly from the negative sideband, while the positive sideband corresponds to the time-delayed saw-tooth chirping signal. De-chirping can be applied by mixing the positive sideband with the original chirping waveform to find the de-chirped frequency
fR and the target range. In this process, the Doppler frequency shift
fd can be easily subtracted to make the range measurement accurate.
3. Experiment
The transmitter used in the lidar system experiment was built based on a Ciena commercial optical transmitter which was originally designed for 10Gb/s optical transmission with electronic domain pre-compensation [
11Y. Zhang, M. O’Sullivan, and R. Hui, “Digital subcarrier multiplexing for flexible spectral allocation in optical transport network,” Opt. Express 19(22), 21880–21889 (2011). [CrossRef] [PubMed]
]. Two 21.42GS/s digital-to-analog converters (DAC) with 6-bit resolution were equipped in the transmitter, which converted the digital modulation signal to analog format to drive the two RF ports of the I/Q modulator. The amplitude of the RF modulating signal can be varied by adjusting the gain of the two RF power amplifiers after the DACs. This transmitter has a 15-bit long internal digital memory so that the maximum temporal length of the arbitrary waveform generated by this transmitter was 1.5298µs. The bandwidth of the I/Q modulator used in the experiment was 10GHz. An integrated tunable laser assembly (Emcore TTX1994) was used as the CW laser source with 100kHz spectral linewidth, The laser wavelength was set at 1549.65nm in the experiment.
The optical signal from the transmitter was collimated onto a spinning disc through free-space optics. The surface of the disc is covered with retro-reflective tape so that the reflected optical signal can be easily collected. The disc was placed approximately 1.6m away from the telescope. The velocity seen by the laser beam can be adjusted by changing either the beam position along the radius of the disc, or the angle between the laser beam and the disc surface normal. Note that the velocity is zero if the disc surface normal is in the same direction of the laser beam, the sign of the velocity seen by the laser can also be flipped through this angle change. The spinning speed of the disc is approximately 2316 rpm (revolutions per minute) and the beam position on the disc is about 4.5cm off the center.
In the receiver, two wideband photodiodes were used after the 90 degree optical hybrid both with 30GHz bandwidth. A real-time oscilloscope (LeCory 8600A) with 20Gs/s sampling rate and approximately 6GHz analog bandwidth was utilized as the ADCs after the RF pre-amplifiers which digitally record the two detected signals. Then the digital signals were processed using a MATLAB program to extract the velocity and range information.
To simultaneously measuring the speed and the distance, two different modulation signals were loaded independently onto the upper and the lower optical sidebands. The upper sideband was loaded with a saw-tooth linear chirping as shown in
Fig. 1(a) with 4.3GHz chirping bandwidth. A raised cosine window function was applied to the chirping signal to minimize the edge effect on each optical pulse. The pulse repetition period was
, and pulse duration was
. The lower sideband was loaded with a single frequency modulation at
, which corresponds to a maximum measurable non-aliasing velocity of 4.149km/s. For the signal on the lower sideband, the pulse width equals to the repetition period
T’ so that the amplitude of the sinusoid is constant. The selection of modulation frequency at
minimizes the phase discontinuity between data frames. The dashed curve in
Fig. 3
shows the signal optical spectrum measured by an optical spectrum analyzer (OSA) with 0.01 nm spectral resolution. The solid trace in
Fig. 3 was the same optical spectrum but measured with a coherent heterodyne technique by mixing the optical signal with another narrow linewidth laser whose central frequency was approximately 6 GHz away from that of the signal laser. The heterodyne spectrum shown in
Fig. 3 was acquired through an RF spectrum analyzer so that the spectral resolution is much better than that measured with the OSA. The optical spectrum in
Fig. 3 clearly shows a constant frequency component on the left side for velocity measurement, a residual optical carrier in the middle, and a wideband frequency chirp on the right side for range measurement. In addition, the residual modulation sideband on the opposite side of the optical carrier is also visible, but with >20dB power suppression ratio. Note that there are also a few spectral lines in
Fig. 3 which were generated from nonlinear mixing between the residual optical carrier, the single-tone modulation and its unsuppressed sideband. This is attributed to the nonlinear transfer function of the system as well as the saturation of the photodiode since the power of the optical LO in the heterodyne spectral measurement was relatively high.
Fig. 3 Signal optical spectrum measured with an OSA with 0.01nm resolution (dashed line) and with coherent heterodyne detection (solid line).
At the lidar receiver, the optical signal reflected from the spinning disc was combined with LO through the 90° optical hybrid. The two output ports are detected by two photodiodes, digitized and combined to reconstruct the complex optical field as described by
Eq. (13). With the perspective of real-time lidar operation, we truncated the recorded time-domain digital signal into lengths equal to the pulse repetition time
, corresponding to the update rate of approximately 653 MHz. This composite signal is converted into frequency domain through an FFT process.
Figure 4
shows a typical spectrum of the complex electrical signals measured at the lidar receiver. This figure covering the full spectral window demonstrates the digital recovery of the complex optical field spectrum comparable to that shown in
Fig. 3, and the two optical sidebands can be utilized for independent purposes.
Fig. 4 Power spectral density of the detected complex electrical signal described by
Eq. (13).
Velocity detection is accomplished by measuring the Doppler induced frequency shift of the single-tone on the lower optical sideband.
Figure 5
shows the expanded view at that single-tone for the target with a positive and a negative velocity by changing the orientation angle of the spinning disc. In the two measurements, both the radial position of the laser beam on the disc and the distance between the lidar and the disc remain unchanged. The solid line and the dashed line in
Fig. 5 represent the spectra of the received signal with the target moving towards and away from the observer, respectively. Because of the Doppler effect, the measured frequency is different from the original single-tone modulation frequency of
, The precise target velocity information can be obtained by measuring this frequency shift. The bottom and the top horizontal axes of
Fig. 5 indicate the frequency
f and the corresponding velocity calculated from
. The target velocities indicated by the solid and the dashed spectra are + 5.81m/s and −7.62m/s, corresponding to the disc orientation angle of 32° and −44°, respectively. Based on the 3-dB width of the spectra, the resolution of the velocity measurement in this experiment is approximately 0.5m/s, which agrees with the theoretical value predicted by
.
Fig. 5 Zoomed-in view at the single-tone modulation frequency for velocity measurement. Solid and dashed lines represent spectra measured when the target moves toward and away from the observer, respectively.
Figure 6
shows the de-chirped spectra for range measurements. De-chirping was performed by mixing the received signal shown on the upper sideband of
Fig. 4 with the original chirping waveform. However, for a moving target, the Doppler-effect shifts the detected signal by a frequency
fd which is a function of the target velocity. Obviously this Doppler frequency shift has an impact on the de-chirping process. The solid and the dashed lines in
Fig. 6 represent de-chirped spectra obtained for the target with the positive and the negative speed as discussed previously. The distance between the lidar transmitter and the target can be calculated by,
Where
fR is the peak frequency of the de-chirped spectrum,
fd is the Doppler frequency shift which is obtained from the velocity measurement described above, and
B is the chirping bandwidth. Note that in
Fig. 6 there is an extra frequency component at
in each spectrum but at a much lower amplitude. This is originated from the back reflection off the fiber APC (angled physical contact) terminal which leads to the telescope. Since optical signal reflected from the fiber terminal is not affected by the Doppler frequency shift caused by the target velocity, the peak frequency of this component is identical for the two measurements. In fact, the actual distance between the fiber terminal (at the telescope) and the target can simply be determined from the frequency difference between the spectral peaks caused by the target and the fiber terminal as
. This resulted in the target distances (from the fiber terminal) of
and
for the measurements with positive and negative velocities of the target, respectively. Apparently the location of the disc was slightly moved in the process of flipping its orientation angle. Based on the 3dB width of the de-chirped spectrum, the range resolution is determined as approximately 2.8cm. This roughly agrees with the calculated range resolution
, where
is defined as the width between the minima on the two sides of the central peak.
Fig. 6 De-chirped signal spectra for distance measurements. Solid and dashed lines represent spectra measured when the target moves toward and away from the observer, respectively.