FMCW LIDAR SYSTEM AND METHOD FOR SIMULTANEOUS RANGE AND VELOCITY MEASUREMENT

Abstract

The invention relates to a system and a method for simultaneous range and velocity measurement in an FMCW LiDAR system. A first light source (16) produces first light having a first frequency that varies according to a first chirp rate. A second light source (18) produces second light having a second frequency that is constant or that varies according to a second chirp rate being different from the first chirp rate. Measuring light obtained by combining the first and second light therefore has two different frequency components during a measurement interval. A splitter (22) separates the measuring light into reference light and output light, and a scanning unit (28) directs the output light towards an object (12) and receives input light that is obtained by reflection of the output light at the object (12). A detector (32) detects a superposition of the reference light and the input light. A computing unit (34) computes unambiguously the range and relative velocity by analyzing beat frequencies resulting from the superposition, wherein ambiguities due to Doppler frequency shifts are removed by performing a decision tree analysis.

Claims

1. An FMCW LiDAR system for simultaneous range and velocity measurement, comprising a first light source configured to produce first light having a first frequency that varies according to a first chirp rate, a second light source configured to produce second light having a second frequency that is constant or that varies according to a second chirp rate that is different from the first chirp rate, an optical combiner configured to combine the first light and the second light, thereby obtaining measuring light having at least two different frequency components during a measurement interval, a splitter configured to separate the measuring light into reference light and output light, a scanning unit configured to direct the output light towards an object along different directions and to receive input light that is obtained by reflection of the output light at the object, a detector configured to detect a superposition of the reference light and the input light, and a computing unit configured to compute unambiguously a range to the object and a relative velocity between the system and the object by analyzing beat frequencies resulting from the superposition detected by the detector, wherein the computing unit is further configured to remove ambiguities due to Doppler frequency shifts by performing a decision tree analysis.

2. The system of claim 1, wherein the second frequency is constant, the first frequency varies according to a chirp rate CR1 during a first portion of the measurement interval and according to a chirp rate CR2 during a second portion of the measurement interval, and wherein the computing unit is configured to compute beat frequencies separately for each portion of the measurement interval, and to compute the range and the velocity by analyzing the beat frequencies measured during both portions of the measurement interval.

3. The system of claim 2, wherein the computing unit is configured, when performing the decision tree analysis, to determine, separately for each of the first and second measurement intervals, how many beat frequencies have been measured, and whether there is a beat frequency that occurs in both measurement intervals.

4. The system of claim 2, wherein the chirp rate CR1 is a positive chirp rate so that the frequency increases during the first portion of the measurement interval, and wherein the chirp rate CR2 is a negative chirp rate so that the frequency decreases during the second portion of the measurement interval.

5. The system of claim 1, wherein the system comprises a third light source configured to produce third light having a third frequency that varies according to a third chirp rate that is different from the first chirp rate and the second chirp rate, and wherein the computing unit is configured to compute beat frequencies for the first light, the second light and the third light, and to compute the range and the velocity by analyzing said beat frequencies.

6. The system of claim 5, wherein analyzing the beat frequencies includes the steps of: a) assigning the beat frequencies in different combinations to the first light, the second light and the third light, b) computing, for at least one combination, a preliminary value for the range and preliminary values for the velocity at least for the first light and the second light, c) determining a combination for which at least two preliminary values for the velocity are sufficiently similar, d) determining final values for the range and for the velocity by adopting the values computed for the combination determined in step c).

7. The system of claim 1, comprising an optical circulator connecting the splitter, the scanning unit and the detector so that the output light is directed towards the scanning unit and the input light is directed towards the detector.

8. A method for simultaneous range and velocity measurement in an FMCW LiDAR system, comprising the following steps: a) producing first light having a first frequency that varies according to a first chirp rate; b) producing second light having a second frequency that is constant or that varies according to a second chirp rate that is different from the first chirp rate; c) combining the first light and the second light, thereby obtaining measuring light having at least two different frequency components during a measurement interval; d) separating the measuring light into reference light and output light; e) directing the output light towards an object along different directions and receiving input light that is obtained by reflection of the output light at the object; f) detect a superposition of the reference light and the input light; g) computing unambiguously a range to the object and a relative velocity between the system and the object by analyzing beat frequencies resulting from the superposition detected by the detector, wherein step g) includes the step of removing ambiguities due to Doppler frequency shifts by performing a decision tree analysis.

9. The method of claim 8, wherein the second frequency is constant, the first frequency varies according to a chirp rate CR1 during a first portion of the measurement interval and according to a chirp rate CR2 during a second portion of the measurement interval, beat frequencies are computed separately for each portion of the measurement interval, and the range and the velocity are computed by analyzing the beat frequencies measured during both portions of the measurement interval.

10. The method of claim 9, wherein, when performing the decision tree analysis, it is determined, separately for each of the first and second measurement intervals, how many beat frequencies have been measured, and whether there is a beat frequency that occurs in both measurement intervals.

11. The method of claim 9, wherein the chirp rate CR1 is a positive chirp rate so that the frequency increases during the first portion of the measurement interval, and wherein the chirp rate CR2 is a negative chirp rate so that the frequency decreases during the second portion of the measurement interval.

12. The method of claim 8, comprising the steps of producing third light having a third frequency that varies according to a third chirp rate that is different from the first chirp rate and the second chirp rate, and computing beat frequencies for the first light, the second light and the third light, and computing the range and the velocity by analyzing said beat frequencies.

13. The method of claim 12, wherein analyzing the beat frequencies includes the steps of: assigning the beat frequencies in different combinations to the first light, the second light and the third light, computing, for at least one combination, a preliminary value for the range and preliminary values for the velocity at least for the first light and the second light, determining a combination for which at least two preliminary values for the velocity are sufficiently similar, determining final values for the range and for the velocity by adopting the values computed for the combination determined in step c).

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0058] Various features and advantages of the present invention may be more readily understood with reference to the following detailed description taken in conjunction with the accompanying drawings in which:

[0059] FIG. 1 is a schematic side view of a vehicle approaching an object which is detected by a scanner system according to the invention;

[0060] FIG. 2 is a top view of the scanner system device shown in FIG. 1;

[0061] FIG. 3 schematically shows the layout of a scanning system according to a first embodiment of the invention;

[0062] FIG. 4 is a graph showing the time dependency of the frequency of the two light sources used in the embodiment shown in FIG. 3;

[0063] is FIG. 5a is a graph showing for the prior art the frequencies of the input light and the reference light of an FMCW LiDAR range measurement system during the up-chirp and the down-chirp interval for a zero relative velocity between the scanning system and the object;

[0064] FIG. 5b is a graph showing the beat frequency for the constellation assumed in FIG. 5a;

[0065] FIG. 5c is a graph showing an idealized power spectrum for the constellation assumed in FIGS. 5a and 5b;

[0066] FIGS. 6a to 6c graphs corresponding to the graphs shown in FIGS. 5a to 5c , respectively, but for a non-zero relative velocity;

[0067] FIG. 7 is a graph similar to FIG. 6a, but assuming that the absolute value of the Doppler frequency is larger than the frequency shift due to the range R;

[0068] FIG. 8a is a graph showing the frequencies of the input light and the reference light of the range measurement system according to the first embodiment of the invention during the up-chirp and the down-chirp intervals;

[0069] FIG. 8b is a graph showing the beat frequencies and the Doppler frequency detected by the detector for the frequencies shown in FIG. 8a of the first embodiment during the up-chirp and the down-chirp interval;

[0070] FIGS. 8c and 8d are graphs showing power idealized power spectra after performing a FFT for the frequencies shown in FIG. 8b ;

[0071] FIGS. 9a to 9c are graphs corresponding to FIGS. 8b , 8c and 8c , respectively, for a scenario in which the same beat frequency occurs during the up- and the down-chirp interval;

[0072] FIG. 10 is a graph corresponding to FIG. 8b for a scenario in which slightly different beat frequencies occur during the up- and the down-chirp interval;

[0073] FIGS. 11a to 11c are graphs corresponding to FIGS. 9a to 9c , respectively, for a scenario in which one beat frequency occurs both in the up- and down-chirp intervals and an additional higher beat frequency occurs only in the down-chirp interval;

[0074] FIGS. 12a to 12c are graphs corresponding to FIGS. 11 a to 11c , respectively, for a similar scenario in which the beat frequencies have a different relationship;

[0075] FIG. 12d is a graph corresponding to FIG. 8a for the scenario shown in FIGS. 12a to 12c ;

[0076] FIGS. 13a to 13c are graphs corresponding to FIGS. 9a to 9c , respectively, for a scenario in which two beat frequencies occur both in the up- and down-chirp intervals;

[0077] FIG. 14 is a graph similar to FIG. 4, but for a second embodiment in which both light sources generate light having varying frequencies;

[0078] FIGS. 15a and 15b are graphs corresponding to FIGS. 13a and 13b , respectively, but for the second embodiment and assuming a scenario in which the largest beat frequency is in the down-chirp interval;

[0079] FIG. 16 schematically shows the layout of a scanning system comprising three light sources according to a third embodiment of the invention;

[0080] FIG. 17 shows the effect of Doppler shifts on an ideal power spectrum measured with a scanning system shown in FIG. 16; and

[0081] FIG. 18 is a diagram showing how the signal-to-noise ratio (SNR) improves with the use of two or three light sources compared to the conventional prior art methods.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

[0082] 1. Introduction

[0083] FIG. 1 is a schematic side view showing a vehicle 10 that approaches an object 12 represented by a tree. The vehicle 10 has a scanner system 14 that scans the environment lying ahead of the vehicle 10 with light beams L11, L21, L31, and L41. From range information associated to each light beam, a three-dimensional image of the environment is computationally reconstructed. In addition, the scanner system 14 determines the relative velocity to the object 12. This information is particularly important if the object 12 is not fixed, but moves, too.

[0084] The information computed by the scanner system 14 about the environment lying ahead of the vehicle 10 may be used, for example, to assist the driver of the vehicle 10 in various ways. For example, warning messages may be generated if a collision of the vehicle 10 with the object 12 threatens. If the vehicle 10 drives autonomously, range and velocity information about the environment lying ahead are required by the algorithms that control the vehicle 10.

[0085] As is apparent in FIG. 1, the scanner system 14 emits the light beams L11 to L41 in different directions in a vertical plane (i.e. the plane of the paper in FIG. 1) so that the environment is scanned in the vertical direction. Scanning also takes place simultaneously in the horizontal direction, as FIG. 2 illustrates in a top view of the scanner system 14. Four light beams L11, L12, L13 and L14 are emitted simultaneously or consecutively in a horizontal plane in different directions.

[0086] For the sake of simplicity it is assumed in FIGS. 1 and 2 that only four light beams Ln1 to Ln4 are emitted vertically and horizontally. However, in real applications the scanner system 14 emits many more light beams. For example, k.2n light beams are preferred, wherein n is a natural number which is typically between 7 and 13 and specifies how many beams are emitted in one of k (horizontal or inclined) planes, wherein k is a natural number which is typically between 1 and 256.

[0087] 2. First Embodiment

[0088] FIG. 3 schematically shows the layout of the scanner system 14 according to a first embodiment of the invention. The scanner system 14 is configured as an FMCW LiDAR system and comprises a first light source 16 and a second light source 18. During operation of the scanner system 14, the first light source 16 produces first light having a first frequency f.sub.chirp that periodically varies (“chirps”) between a lower frequency f.sub.1 and a higher frequency f.sub.h. Typically, the first light source 16 comprises a tunable laser light source or a laser light source emitting light with a fixed frequency which is then frequency modulated, as this is well known in the art as such. However, other types of light sources are contemplated as well.

[0089] The graph of FIG. 4 shows how the frequency f.sub.chirp of the first light source 16 varies over time tin this embodiment. Each measurement interval having a chirp duration T is divided into two halves of equal length T/2. During the first interval, the frequency f.sub.chirp increases is linearly with a constant and positive up-chirp rate r.sub.chirp, i.e. df.sub.chirp/dt=r.sub.chirp. The first interval will therefore be referred to in the following as up-chirp interval. During the second interval, the frequency f.sub.chirp decreases linearly with a constant negative down-chirp rate−r.sub.chirp, i.e. df.sub.chirp/dt=−r.sub.chirp. This results in a triangular wave like frequency variation as shown in FIG. 4.

[0090] In this embodiment, the second light source 18 continuously produces second light having a constant frequency f.sub.cw during the up-chirp and the down-chirp intervals. Also the second light source 18 typically comprises a laser light source, but other types of light sources are also contemplated. In FIG. 4, the frequency f.sub.cw is indicated with a broken line.

[0091] During the round-trip time 2R/c, which the first light requires to get to the object 12 and back to the scanner system 14, the frequency change of the first light should be so small that it can always be clearly distinguished from the frequency f.sub.cw of the second light. For a typical distance R=150 m and a chirp rate r.sub.chirp=10.sup.14 s.sup.−2, the wavelength of the first light may be 1550 nm, and the wavelength of the second light may be 1554 nm.

[0092] Both light sources 16, 18 are connected to an optical combiner 20 that combines the first light and the second light to measuring light having two different frequency components, namely the varying frequency f.sub.chirp and the constant frequency f.sub.cw. The optical combiner 20 may be realized as an optical 2:1 multiplexer having two input ports and one output port.

[0093] The measuring light enters a splitter 22 that separates the measuring light into reference light (sometimes also referred to as “local oscillator”) and output light. In this embodiment, the output light passes an optical amplifier 24 and an optical circulator 26 that guides the amplified output light towards a scanning unit 28.

[0094] Instead of using an amplifier 24 that amplifies only the output light, it is also possible, for example, to amplify the measuring light before it enters the splitter 22, to use at least one amplifier arranged in at least one light path between the light sources 16, 18 and the optical combiner 20, or to dispense with the amplifier 24 completely.

[0095] The scanning unit 28 directs the output light towards the object 12—in FIG. 3 represented by a moving car—along different directions, as it has been explained above with reference to FIGS. 1 and 2. While some part of each output light beam is usually absorbed by the object 12, another part is reflected. As long as the surface of the object 12 is not perfectly specular, diffuse reflection occurs which ensures that some output light is reflected at the object 12 such that it propagates back towards the scanning unit 28. The reflected output light, which is referred to in the following as input light, is received by the scanning unit 28 and guided towards the optical circulator 26.

[0096] The optical circulator 26 passes the input light towards a further combiner 30 that combines the reference light, which was separated from the measuring light by the splitter 22, with the input light. A detector 32 arranged behind the further combiner 30 thus detects a superposition of the reference light and the input light. The detector 32 may be configured as a balanced detector, as this is known in the art as such. The electric signals produced by the detector 32 are fed to a computing unit 34 that computes the range R to the object and a relative velocity v between the scanner system 14 and the object 12 by analyzing beat frequencies resulting from the superposition detected by the detector 32.

[0097] Since no light has to be routed from the further combiner 30 towards the splitter 22, using the optical circulator 26 is not mandatory. Often it suffices to use simpler polarization sensitive beam splitting elements instead of the optical circulator 26.

[0098] a) FMCW LiDAR principle

[0099] For a basic explanation of the FMCW LiDAR principle, the second light produced by the second light source 18 will be ignored in the following discussion.

[0100] FIG. 5a illustrates the frequencies f.sub.chirp of the reference light (solid line) and the frequency f.sub.i of the input light (broken line) for the case that there is no relative velocity between the scanner system 14 and the object (v=0). It can be seen that the frequency f.sub.i is shifted to the right on the time axis, which is a consequence of the fact that the input light is delayed by the round-trip time T.sub.p=2R/c it takes for the output light to reach the object 12 and to return to the scanner system 14 as input light.

[0101] Due to the superposition of the reference light and the input light having similar frequencies, the detector 32 detects a beat signal having a beat frequency f.sub.b which is equal to the frequency difference at any given time during the measurement interval of length T, as this is shown in FIG. 5a. Since the delay time T.sub.p is proportional to the beat frequency f.sub.b, the range R to the object can be computed from the beat frequency f.sub.b according to

[00001]$\begin{array}{cc}R=\frac{Tc}{4\Delta f}{f}_b& \left(1\right)\end{array}$

in which f=f.sub.h−f.sub.l is the total change in chirp frequency during the up-chirp or down-chirp interval of length T/2. Sometimes, this difference f is referred to as the bandwidth of the LiDAR system.

[0102] FIG. 5b is a graph that illustrates how the beat frequency f.sub.b changes during a measurement interval of length T. Since the absolute value of the beat frequency f.sub.b is identical during the up-chirp and down-chirp intervals, it suffices to measure the beat frequency once during a measurement interval of length T. FIG. 5b also shows that there is a dead time of length T.sub.p each time the sign of the chirp rate changes. During the dead times no useful measurement can be made.

[0103] The beat frequency f.sub.b can be determined by measuring the intensity I(t) at the detector 32 as a function of time t, followed by a Fast Fourier Transform (FFT) that yields a frequency spectrum P(f) as shown in FIG. 5c for an ideal case, i.e. in the absence of noise and with an infinitely long duration of the measurement interval T. In real systems, it is quite demanding to detect the beat frequency in the noise, because the input light received by the detector 32 has an extremely low intensity.

[0104] b) Doppler effect

[0105] In the following it will be assumed that there is a relative velocity between the scanner system 14 and the object 12 (v≠0). As a result of the relative movement, a Doppler shift occurs that has to be taken into account. The velocity v is defined along the line of light propagation. Sometimes the velocity defined in this manner is referred to as radial velocity.

[0106] As it is illustrated in FIG. 6a, the frequency f.sub.i of the input light (broken line) is not only shifted along the time axis relative to the frequency f.sub.chirp of the reference light, but also vertically as a result of the Doppler shift. Here it is assumed that the velocity v is negative, e.g. because the object approaches the scanner system 14. Consequently, the Doppler frequency f.sub.D is positive, which implies an increase of the frequency f.sub.i of the input light.

[0107] As a result of the Doppler shift, the beat frequency is no longer the same during the up-chirp and the down-chirp interval. Instead, in this example, the beat frequency increases during the up-chirp interval and decreases during the down-chirp interval. The resulting beat frequencies f.sub.b1 and f.sub.b2 are also shown in the graphs of FIG. 6b and FIG. 6c that correspond to FIGS. 5b and 5c illustrating the stationary case without Doppler shift.

[0108] In FIGS. 6a to 6c is assumed that the absolute value of the Doppler frequency |f.sub.D| is smaller than the frequency shift f.sub.R due to the range R. Then f.sub.R and f.sub.D can be computed from the detected beat frequencies f.sub.b1 and f.sub.b2 according to the following equations:

[00002]$\begin{array}{cc}{f}_R=\frac{f_{b1}+f_{b2}}{2}{f}_D=\frac{f_{b2}-f_{b1}}{2}& \left(2\right)\end{array}$

[0109] The range induced frequency shift f.sub.R is thus the arithmetic mean of the detected beat frequencies f.sub.b1 and f.sub.b2, while the Doppler frequency f.sub.D is one half of the difference between the beat frequencies f.sub.b1 and f.sub.b2.

[0110] From the range induced frequency shift f.sub.R the range R can be computed according to equation (1), setting f.sub.b=f.sub.R. The velocity v can be derived from the Doppler frequency f.sub.D according to

[00003]$\begin{array}{cc}v=f_D.Math.\frac{\lambda }{2}& \left(3\right)\end{array}$

with λ being the center wavelength of the output light.

[0111] If the absolute value of the Doppler frequency |f.sub.D| is larger than the frequency shift f.sub.R due to the range R, a situation as exemplarily shown in FIG. 7 occurs. In that case, f.sub.R and f.sub.D can be computed from the detected beat frequencies f.sub.b1 iand f.sub.b2 from equations (4):

[00004]$$\begin{gathered} \begin{array}{cc}{f}_R=.Math.\frac{f_{b2}-f_{b1}}{2}.Math. \\ {f}_D=\{\begin{array}{c}\frac{f_{b2}+f_{b1}}{2}\mathrm{if}{f}_{b2}<f_{b1}\\ -\left(\frac{f_{b2}+f_{b1}}{2}\right)\mathrm{if}{f}_{b2}>f_{b1}\end{array}& \left(4\right)\end{array} \end{gathered}$$

[0112] The problem with this kind of evaluation is that the FFT produces only two frequency peaks as shown in FIG. 6c. However, it is unknown whether the absolute frequency of the Doppler frequency |f.sub.D| is smaller or larger than the frequency shift f.sub.R due to the range R. Therefore, the above distinction between two different cases cannot be reasonably made, and thus it is not clear whether the Doppler frequency f.sub.D is positive or negative.

[0113] This ambiguity is not relevant in RADAR applications because the wavelength and the modulation frequencies are quite different. In RADAR applications the absolute value of the Doppler frequency |f.sub.D| is always significantly smaller than f.sub.R. This is different in LiDAR applications. For example, in a LiDAR scanner system using a wavelength of 1550 nm, the Doppler frequency f.sub.D is about 65 MHz for a speed of 180 km/h, while the beat frequencies (which are quite similar to f.sub.R) are typically between 5 MHz and 100 MHz. In other words, the absolute value of the Doppler frequency |f.sub.D| is just within the range of possible beat frequencies, and there is no way to tell whether the absolute value |f.sub.D| is larger or smaller than f.sub.R. Without this knowledge, it is not possible to use the correct equations for the computation of the range R and the velocity v.

[0114] c) Adding second light with constant wavelength

[0115] As has been explained above with reference to FIGS. 3 and 4, the scanner system 14 of the first embodiment comprises an additional second light source 18 that produces second light having a constant frequency f.sub.cw. FIG. 8a illustrates the frequency f chirp of the first light (bold solid line) and the frequency f.sub.cw of the second light (bold broken line). In addition to what is shown in FIG. 4, FIG. 8a also illustrates the two frequency components (thin lines) of the input light, i.e. after reflection at the object 12, for the case of an approaching object (i.e. v<0). It can be seen that both frequency curves f.sub.chirp and f.sub.cw are shifted “upward”, i.e. towards larger frequencies, as a result of the Doppler frequency f.sub.D which is positive in this example.

[0116] It can be derived from basic equations of electromagnetic theory that interference of the reference light with the input light results in the following Intensity I(t) at the detector 32:

[00005]$\begin{array}{cc}I\left(t\right)\sim A+B\cos \left(\left(r_{chirp}.Math.\frac{2R}{c}+4\pi \frac{v}{\lambda }\right)t\right)+C\cos \left(4\pi .Math.\frac{v}{\left(\lambda -\Delta \lambda \right)}t\right)& \left(5\right)\end{array}$

with A, B and C being constants and λ being the wavelength difference corresponding to the frequency difference ω=(f.sub.chirp−f.sub.cw). In this derivation it has been assumed that cosine terms containing ω are constants, because ω is very large compared to the other terms (e.g. 3,000 GHz, see further above).

[0117] Since λ is small compared to λ, equation (5) can be rewritten as

I(t)˜A+B cos((αR+βv)t)+C cos(βvt)   (6) [0118] with α=r.sub.chirp.Math.2/c and β=4π/λ

[0119] From equation (6) it can be concluded that there are generally two frequency peaks in the Fourier spectrum during the up-chirp interval and two frequency peaks during the down-chirp interval. One of the two beat frequencies depends both on the range R and the velocity v and is the result of the interference between the first light and its reflected portion that is affected by the Doppler effect. It is the same beat frequency that would be observed if the second light was absent.

[0120] The other beat frequency corresponds to the Doppler frequency f.sub.D (see equation (3)) and is thus depends only on the velocity v. This frequency is the result of the interference between the second light with and without Doppler shift.

[0121] FIG. 8b shows the time dependency of the frequency components of the input light in a representation similar to FIGS. 5b and 6b, and FIGS. 8c and 8d show the frequency peaks obtained in a first measurement during the up-chirp interval and in a second measurement during the down-chirp interval, respectively. These peaks occur in the

[0122] Fourier spectrum obtained by a FFT of an electrical signal s(t) produced in the detector 32 and corresponding to I(t) according to equations (5) and (6).

[0123] By comparing FIGS. 8c and 8d, it can be seen that the left peak (smallest frequency) must represent the Doppler frequency f.sub.D, because it is present in both spectra. The other peak in each spectrum must represent the beat frequency f.sub.b1 during up-chirp and f.sub.b2 during down-chirp.

[0124] There are other scenarios in which it is not that easy to identify the frequencies that are represented by the observed peaks. However, by carefully analyzing the different scenarios that may occur, it is possible to exclude unrealistic scenarios so that correct values for the range R and the velocity v can be computed. This will be explained in the next section.

[0125] d) Analysis of Possible Scenarios

[0126] Two Peaks

[0127] There are scenarios in which only two peaks are observed in each Fourier spectrum, a first peak during the up-chirp interval and a second peak during the down-chirp interval. These peaks may be identical or different.

[0128] FIGS. 9a to 9c illustrates a first scenario in which the same frequency component f.sub.b1=f.sub.b2 occurs in the spectrum during the up-chirp and the down-chirp interval. One explanation for observing only one frequency peak is that the velocity v and therefore also the Doppler frequency f.sub.D equals zero. With this assumption, the range R can be computed according to equation (1) with f.sub.b=f.sub.b1=f.sub.b2. In FIG. 9a it is assumed that the values for T, f.sub.b and f result in a range R=50 m.

[0129] Another explanation for observing only one frequency peak would be that both beat frequencies are zero, i.e. f.sub.b1=f.sub.b2=0. This would imply a range R=0 and lead to a velocity v≠0, for example 27 m/s. However, such a result is not plausible, because a high relative speed cannot occur if the object 12 is positioned directly in front of the scanner system 14 (R=0). This solution must therefore be discarded.

[0130] This implies that if two peaks occur in the spectrum, the velocity v must be zero, and the range R is computed according to equation (1).

[0131] FIG. 10 illustrates a two peak scenario in which the two frequency components f.sub.b1 and f.sub.b2 are slightly different. A reasonable explanation for this scenario is that the Doppler frequency f.sub.D is so small such that it cannot be detected. Since |f.sub.D| is much smaller than the beat frequencies, the range R and the velocity v have to be computed according to equations (1) to (3).

[0132] A scenario with two significantly different peaks during up-chirp and down-chirp interval cannot occur. Such peaks could only be a result of a strong Doppler shift, but the Doppler frequency f.sub.D cannot be equal to two different beat frequencies.

[0133] A scenario in which there are two peaks in the up-chirp and no peak in the down-chirp interval (or vice versa) cannot occur. Such invalid measurements could only be explained by artefacts or noise and must be discarded.

[0134] Three Peaks

[0135] If three peaks occur, the situation is more complicated, but there is still an algorithm making it possible to unambiguously determine the range R and the velocity v.

[0136] Type 1

[0137] FIGS. 11a to 11c are graphs similar to FIGS. 9a to 9c that illustrate a scenario with one identical peak both in the up-chirp and in the down-chirp interval and an additional peak in the down-chirp interval. The only reasonable explanation for this scenario is that the Doppler frequency f.sub.D appears both in the up-chirp and the down-chirp interval, but coincides with the beat frequency f.sub.b1 in the up-chirp interval.

[0138] From the lower frequency peak in the down-chirp interval (see FIG. 11c) the absolute value of the Doppler frequency f.sub.D can be derived, but not its sign. However, from equation (6) the condition

f.sub.b1=|α R+βb|=|βb|  (7)

can be derived for the up-chirp interval, because both frequencies have the same magnitude. From this condition it becomes clear that v must be negative, because all other quantities (including a) are positive. Then equation (7) can be rewritten as

f.sub.b1=α R+βv=−βv   (8)

from which the quantities R and v can be computed as

[00006]$\begin{array}{cc}R=\frac{2f_{b1}}{\alpha }\mathrm{and}v=-\frac{f_{b1}}{\beta }& \left(9\right)\end{array}$

[0139] FIG. 11a indicates exemplary values v=−13 m/s and R=50 m.

[0140] From equation (6) follows that also for the second beat frequency f.sub.b2=|α R+βv|. Using equation (9) yields the relationship f.sub.b2=3f.sub.b1.

[0141] As a matter of course, the same considerations apply, mutatis mutandis, if the there are two peaks in the up-chirp interval and only one peak in the down-chirp interval.

[0142] Type 2

[0143] FIGS. 12a to 12c are graphs similar to FIGS. 11a to 11c that illustrate another scenario with one identical peak both in the up-chirp and in the down-chirp interval and an additional peak in the up-chirp interval. The two scenario types differ with regard to their relationship between the beat frequencies f.sub.b1 and f.sub.b2. While in the type 1 scenario f.sub.b2=3f.sub.b1, the relationship f.sub.b2=2f.sub.b1 holds for the type 2 scenario.

[0144] The only reasonable explanation for this scenario is that the frequency shift caused by the range R is completely compensated for by the Doppler shift during the up-chirp interval, as this is illustrated in FIG. 12d. For equation (7) this implies that |α R+βv|=0. As r.sub.chirp and thus also α are positive during the up-chirp interval, the velocity v must be negative.

[0145] The values for R and v are then given by

[00007]$\begin{array}{cc}R=\frac{f_{b2}}{\alpha }\mathrm{and}v=-\frac{f_D}{\beta }& \left(10\right)\end{array}$

[0146] As a matter of course, the same considerations apply, mutatis mutandis, if there are two peaks in the up-chirp interval and only one peak in the down-chirp interval. In that case the velocity v will be positive.

[0147] Four Peaks

[0148] In most measurements there will be four peaks in the spectra, two peaks in the up-chirp interval and two peaks in the down-chirp interval. All scenarios in which there are three or all four peaks in either the up-chirp or down-chirp interval are invalid because there is no reasonable explanation for them. Such measurements have to be discarded.

[0149] If two peaks during the up-chirp interval and two peaks during the down-chirp interval are observed, there are different scenarios that can to be easily distinguished.

[0150] If all four peaks are different, the measurement has to be discarded, because there is no reasonable explanation for such a scenario. The Doppler peak must always be equal in both intervals.

[0151] The same is true for scenarios in which the peaks in the up-chirp interval are the same as in the down-chirp interval. If there is a Doppler shift, the beat frequencies f.sub.b1 and f.sub.b2 must be different.

[0152] The two types of valid four peak measurements have already been discussed above in section 2.b), but will be briefly explained again in the following.

[0153] Type 1

[0154] FIGS. 8a to 8d discussed above illustrate a scenario in which the velocity v and thus the Doppler frequency f.sub.D is small. The Doppler frequency f.sub.D can be clearly identified as the identical peak occurring in both measurements. For small Doppler frequencies, the range R and the velocity v can be computed according to equations (1) to (3). The sign of the v velocity depends on whether f.sub.b1 is larger or less than f.sub.b2. In FIGS. 8a to 8d it is assumed that f.sub.b2>f.sub.b1. According to equation (2), this results in a positive Doppler frequency f.sub.D and thus in a positive velocity v according to equation (3).

[0155] From equation (2) it can be derived that in this case the Doppler frequency f.sub.D must be twice the difference between the beat frequencies f.sub.b2 and f.sub.b1. If this condition is grossly violated, the measurement is invalid and has to be discarded.

[0156] Type 2

[0157] FIGS. 13a to 13c illustrate a scenario in which the velocity v and thus the Doppler frequency f.sub.D is large (|f.sub.D|>f.sub.R). The Doppler frequency f.sub.D can again be clearly identified as the identical peak occurring in both measurements. For large Doppler frequencies, the range R and the velocity v can be computed according to equations (1), (3) and (4). The sign of the v velocity depends on whether f.sub.b1 is larger or less than f.sub.b2. In FIGS. 13a to 13c it is assumed that f.sub.b2>f.sub.b1. According to equation (4), this results in a negative Doppler frequency f.sub.D and thus in a negative velocity v according to equation (3).

[0158] From equation (4) it can be derived that in this case the Doppler frequency f.sub.D must be the arithmetic mean of the beat frequencies f.sub.b2 and f.sub.b1. If this condition is grossly violated, the measurement is invalid and has to be discarded.

[0159] If the Doppler frequency is smaller than both beat frequencies f.sub.b1 and f.sub.b1, this can only be the result of an invalid measurement, because then the Doppler frequency f.sub.D cannot be equal to the arithmetic mean of the beat frequencies f.sub.b2 and f.sub.b1.

[0160] e) Intermediate Result

[0161] It has been demonstrated that by adding second light having a constant wavelength, it is possible to resolve the Doppler ambiguity. In all valid scenarios values for the range R and the velocity v can be is unambiguously computed from the measured spectra during the up-chirp and the down-chirp interval.

[0162] f) Alternative Approaches

[0163] As can be seen in FIG. 8b, for example, the Doppler frequency can be observed continuously at any time during the measurement interval of length T, whereas the beat frequencies are interrupted between the up-chirp and down-chirp intervals. This can be used to clearly distinguish the Doppler frequency f.sub.D which removes any ambiguities.

[0164] However, identifying the Doppler frequency during the short interruptions between up-chirp and down-chirp intervals can be difficult due to noise issues. The identification can be made more reliable by selecting a higher or lower intensity for the second light.

[0165] Another option is to apply a slow amplitude modulation to the second light so that it becomes easier to identify the frequency peaks directly.

[0166] 3. Second Embodiment

[0167] In the first embodiment discussed above the second light source 18 produces second light is having a constant frequency f.sub.cw during the up-chirp and the down-chirp intervals. As it has been demonstrated above, the range R and the velocity v can then be unambiguously computed.

[0168] It is also possible to use a second light source 18′ that produces second light having a frequency f.sub.chirp2 that also varies over time similar to the first light produced by the first light source 16. In that case the frequency change must be different, i.e. for the chirp rates r.sub.chirp1, r.sub.chirp2 of the two light sources 16 and 18 the condition |r.sub.chirp1|≠|r.sub.chirp| must apply, as this is illustrated in FIG. 14. The measurement procedure is then basically the same as described in the context of the first embodiment. Therefore, the intensity I(t) at the detector 32 is measured as a function of time t, followed by a Fast Fourier Transform (FFT) that yields one frequency spectrum P(f) for the up-chirp interval and one spectrum for the down-chirp interval.

[0169] Unlike in the first embodiments, in which also two or three peaks can be observed in the two spectra, there will be always two peaks in each spectrum. Nevertheless, there are two main scenarios with a couple of different sub-scenarios that can then be distinguished.

[0170] a) Largest beat frequency is in down-chirp interval

[0171] If the largest beat frequency occurs in the down-chirp interval, as this is shown in FIGS. 15a and 15b, the velocity v must be negative. The range R and the velocity v can then be computed according to equations (11):

[00008]$\begin{array}{cc}R=c\frac{f_{b4}-f_{b3}}{2.Math.\left(r_{chirp2}-r_{chirp1}\right)}v=-\frac{R.Math.\left(r_{chirp1}+r_{chirp2}\right)-f_{b4}-f_{b{3}_{}}}{f_{b1}+f_{b2}}& \left(11\right)\end{array}$

[0172] b) Largest beat frequency is in up-chirp interval

[0173] If the largest beat frequency occurs in the up-chirp interval, the velocity v must be positive. The range R and the velocity v can then be computed according to equations (12):

[00009]$\begin{array}{cc}R=c\frac{f_{b2}-f_{b1}}{2.Math.\left(r_{chirp2}-r_{chirp1}\right)}v=-\frac{R.Math.\left(r_{chirp1}+r_{chirp2}\right)-f_{b2}-f_{b1}}{f_{b1}+f_{b2}}& \left(12\right)\end{array}$

[0174] For the beat frequencies f.sub.b3 and f.sub.b4 the following definition applies:

[0175] If there is a beat frequency in the down-chirp interval between f.sub.b1 and f.sub.b2, this frequency is f.sub.b4, and the other beat frequency in the down-chirp interval is f.sub.b3. If there is no beat frequency in the down-chirp interval between f.sub.b1 and f.sub.b2, f.sub.b4 is the smaller of the two beat frequencies, and f.sub.b3, and the other beat frequency in the down-chirp interval is f.sub.b3.

[0176] The range R and the velocity v can thus be unambiguously computed.

[0177] As a special case, the two chirp rates may have the same absolute value, but opposite signs, i.e. r.sub.chirp2=−r.sub.chirp1. Equations (12) then simplify correspondingly.

[0178] 4. Third Embodiment

[0179] In the first and second embodiments discussed above, the scanner system 14 comprises two light sources producing light having different chirp rates (in the first embodiment, one of the chirp rates is zero). FIG. 16 shows the architecture of a scanner system 14 that comprises not two, but three light sources 16, 17 and 18 having different chirp rates, and preferably emitting light beams having different center wavelengths as well, although this is not mandatory. With equal center wavelengths, it is difficult to identify the desired beat frequencies among the beat frequencies being a result of the superposition of the three reference light portions. One approach to overcome this problem would be to provide three completely distinct light paths, and to ensure that the light spots on the object do not overlap. However, this significantly adds to the system complexity and costs.

[0180] As will become apparent from the following description, the configuration with three light sources has the benefit that not two, but only one measurement is required per measurement interval of length T. Nevertheless, the range R and the velocity v can be o unambiguously computed as in the other embodiments.

[0181] FIG. 17 shows an ideal power spectrum P(f) (i.e. a spectrum in the absence of any noise and assuming an infinite duration T of the measurement interval) comprising three beat peaks at frequencies f.sub.b1, f.sub.b2 and f.sub.b3 without Doppler effect (solid lines) and at frequencies f.sub.b1D, f.sub.b2D and f.sub.b3D shifted by the Doppler effect (dashed lines). Here it is is assumed that the chirp rate r.sub.chirp of the first light source 16 is 0 Hz so that the beat frequency f.sub.b1 is also 0 Hz. For the other two light sources 17, 18 it is assumed that the beat frequencies f.sub.b2 and f.sub.b1 are 33 MHz and 133 MHz, respectively.

TABLE-US-00001 TABLE 1 Frequency assignment and preliminary range and velocity computation Preliminary Preliminary range velocity computation computation according according Case Frequency Assignment to equation to equations A f.sub.b1D = f.sub.min, f.sub.b2D = f.sub.mid, f.sub.b3D = f.sub.max (13) (16) B f.sub.b1D = f.sub.min, f.sub.b2D = f.sub.mid, f.sub.b3D = f.sub.max (13) (17) C f.sub.b1D = f.sub.min, f.sub.b2D = f.sub.mid, f.sub.b3D = f.sub.max (14) (18) D f.sub.b1D = f.sub.mid, f.sub.b2D = f.sub.min, f.sub.b3D = f.sub.max (14) (18) E f.sub.b1D = f.sub.max, f.sub.b2D = f.sub.min, f.sub.b3D = f.sub.mid (14) (18) F f.sub.b1D = f.sub.max, f.sub.b2D = f.sub.mid, f.sub.b3D = f.sub.min (15) (19) A f.sub.b1D = f.sub.min, f.sub.b2D = f.sub.mid, f.sub.b3D = f.sub.max (13) (16) G f.sub.b1D = f.sub.max, f.sub.b2D = f.sub.min, f.sub.b3D = f.sub.mid (15) (19) H f.sub.b1D = f.sub.mid, f.sub.b2D = f.sub.min, f.sub.b3D = f.sub.max (15) (19) I f.sub.b1D = f.sub.max, f.sub.b2D = f.sub.mid, f.sub.b3D = f.sub.min (15) (20)

[0182] In the presence of a radial velocity v≠0 the beat frequencies shown with solid lines will be differently shifted, which is indicated by thin horizontal arrows in FIG. 17. If the Doppler frequency f.sub.D is −45 MHz, the beat frequency f.sub.b1 of the first light source 16 is shifted to f.sub.b1D=|045| MHz=45 MHz, the beat frequency f.sub.b2 of the second light source 17 is shifted to f.sub.b2D=33-45| MHz=12 MHz, and the beat frequency f.sub.b3 of the third light source 18 is shifted to f.sub.b3D=|133-45| MHz=88 MHz.

[0183] In a real measurement, it is initially not known which Doppler-shifted beat peak belongs to which light source. In particular, it is also not known which beat peaks experience a “frequency reflection” at f=0 Hz. In the example shown, this is the case for the first light source 16 and for the second light source 17, but not for the third light source 18.

[0184] Table 1 shows a first step of a solution scheme according to this embodiment for the problem how the peaks can be correctly assigned to the three light sources. With three light sources 16, 17 and 18, the cases A to I can be distinguished. As long as the beat peaks do not overlap randomly, each power spectrum P(f) obtained by FFT contains three beat peaks having frequencies f.sub.min, f.sub.mid, and f.sub.max in increasing order. The Doppler affected beat peaks associated with the light sources 16, 17 and 18 are designated with f.sub.b1D, f.sub.b2D and f.sub.b3D. The second column in the table contains possible assignments of the measured beat frequencies to the peaks in the power spectrum. The third column refers to the equation for calculating a preliminary range R, and the fourth column refers to the equation for calculating preliminary velocities v.sub.1, v.sub.2 and v.sub.3 for each light source 16, 17, 18.

[0185] In this first part of the solution scheme the following equations for calculating the distance R are used:

[0186] Cases A and B:

[00010]$\begin{array}{cc}R=\frac{c}{2}\frac{\frac{f_{b2D}-\left(f_{b1D}.Math.\frac{{\lambda }_1}{{\lambda }_2}\right)}{r_{chirp1}+r_{chirp2}}}{1-\frac{r_{chirp1}.Math.\left(1+\frac{{\lambda }_1}{{\lambda }_2}\right)}{r_{chirp1}+r_{chirp2}}}& \left(13\right)\end{array}$

[0187] Cases C, D and E:

[00011]$\begin{array}{cc}R=\frac{c}{2}\frac{\frac{f_{b2D}+\left(f_{b1D}.Math.\frac{{\lambda }_1}{{\lambda }_2}\right)}{-r_{chirp1}+r_{chirp2}}}{1+\frac{r_{chirp1}.Math.\left(1-\frac{{\lambda }_1}{{\lambda }_2}\right)}{-r_{chirp1}+r_{chirp2}}}& \left(14\right)\end{array}$

[0188] Cases F, G, H and I:

[00012]$\begin{array}{cc}R=\frac{c}{2}\frac{\frac{f_{b2D}-\left(f_{b1D}.Math.\frac{{\lambda }_1}{{\lambda }_2}\right)}{-r_{chirp1}-r_{chirp2}}}{1+\frac{r_{chirp1}.Math.\left(1+\frac{{\lambda }_1}{{\lambda }_2}\right)}{-r_{chirp1}-r_{chirp2}}}& \left(15\right)\end{array}$

[0189] The following equations for calculating preliminary velocities v.sub.1, v.sub.2 and v.sub.3 are used:

[0190] Case A:

[00013]$\begin{array}{cc}{v}_1={\lambda }_1\left(\frac{1}{2}{f}_{b1D}-\frac{R.Math.r_{chirp1}}{c}\right)v_2={\lambda }_2\left(\frac{1}{2}{f}_{b2D}-\frac{R.Math.r_{chirp2}}{c}\right)v_3={\lambda }_3\left(\frac{1}{2}{f}_{b3D}-\frac{R.Math.r_{chirp3}}{c}\right)& \left(16\right)\end{array}$

[0191] Case B:

[00014]$\begin{array}{cc}{v}_1=-{\lambda }_1.Math.-\frac{1}{2}{f}_{b1D}+\frac{R.Math.r_{chirp1}}{c}.Math.v_2=-{\lambda }_2.Math.-\frac{1}{2}{f}_{b2D}+\frac{R.Math.r_{chirp2}}{c}.Math.v_3=-{\lambda }_3.Math.-\frac{1}{2}{f}_{b3D}+\frac{R.Math.r_{chirp3}}{c}.Math.& \left(17\right)\end{array}$

[0192] Cases C, D and E:

[00015]$\begin{array}{cc}{v}_1=-{\lambda }_1.Math.\frac{1}{2}{f}_{b1D}+\frac{R.Math.r_{chirp1}}{c}.Math.v_2=-{\lambda }_2.Math.-\frac{1}{2}{f}_{b2D}+\frac{R.Math.r_{chirp2}}{c}.Math.v_3=-{\lambda }_3.Math.-\frac{1}{2}{f}_{b3D}+\frac{R.Math.r_{chirp3}}{c}.Math.& \left(18\right)\end{array}$

[0193] Cases F, G and H:

[00016]$\begin{array}{cc}{v}_1=-{\lambda }_1.Math.\frac{1}{2}{f}_{b1D}+\frac{R.Math.r_{chirp1}}{c}.Math.v_2=-{\lambda }_2.Math.\frac{1}{2}{f}_{b2D}+\frac{R.Math.r_{chirp2}}{c}.Math.v_3=-{\lambda }_3.Math.-\frac{1}{2}{f}_{b3D}+\frac{R.Math.r_{chirp3}}{c}.Math.& \left(19\right)\end{array}$

[0194] Case I:

[00017]$\begin{array}{cc}{v}_1=-{\lambda }_1\left|\frac{1}{2}{f}_{b1D}+\frac{R.Math.r_{chirp1}}{c}\right|v_2=-{\lambda }_2\left|\frac{1}{2}{f}_{b2D}+\frac{R.Math.r_{chirp2}}{c}\right|v_3=-{\lambda }_3\left|\frac{1}{2}{f}_{b3D}+\frac{R.Math.r_{chirp3}}{c}\right|& \left(20\right)\end{array}$

[0195] Table 2 shows a second step of the solution scheme according to this embodiment. This table represents a binary decision tree that is used to determine, on the basis of the preliminary values for the velocity v.sub.1, v.sub.2 (see equation 21), which case A to I is present, and how then final values for the range R and the velocity v can to be calculated.

[0196] In this calculation is has been taken into account that the measured beat peaks are not arbitrarily sharp, but can only be measured with the frequency resolution 1/T. Due to this limited frequency resolution, a certain reconstruction error rate must be expected, since the calculations of the preliminary velocities v.sub.1, v.sub.2 and v.sub.3 cannot be compared with each other with arbitrary accuracy.

TABLE-US-00002 TABLE 2 Decision tree for computing final values forR and v Case Frequency Assignment Computation of final values for range R and velocity v A f.sub.b1D = f.sub.min, f.sub.b2D = f.sub.mid, f.sub.b3D = f.sub.max IF (21) = true with X = A, THEN compute R according to (13) and v according to (22), ELSE Case B B f.sub.b1D = f.sub.min, f.sub.b2D = f.sub.mid, f.sub.b3D = f.sub.max IF (21) = true with X = B, THEN compute R according to (13) and v according to (23), ELSE Case C C f.sub.b1D = f.sub.min, f.sub.b2D = f.sub.mid, f.sub.b3D = f.sub.max IF (21) = true with X = C, THEN compute R according to (14) and v according to (24), ELSE Case D D f.sub.b1D = f.sub.mid, f.sub.b2D = f.sub.min, f.sub.b3D = f.sub.max IF (21) = true with X = D, THEN compute R according to (14) and v according to (24), ELSE Case E E f.sub.b1D = f.sub.max, f.sub.b2D = f.sub.min, f.sub.b3D = f.sub.mid IF (21) = true with X = E, THEN compute R according to (14) and v according to (24), ELSE Case F F f.sub.b1D = f.sub.max, f.sub.b2D = f.sub.mid, f.sub.b3D = f.sub.min IF (21) = true with X = F, THEN compute R according to (15) and v according to (24), ELSE Case G G f.sub.b1D = f.sub.max, f.sub.b2D = f.sub.min, f.sub.b3D = f.sub.mid IF (21) = true with X = G, THEN compute R according to (15) and v according to (24), ELSE Case H H f.sub.b1D = f.sub.mid, f.sub.b2D = f.sub.min, f.sub.b3D = f.sub.max IF (21) = true with X = H, THEN compute R according to (15) and v according to (24), ELSE Case I I f.sub.b1D = f.sub.max, f.sub.b2D = f.sub.mid, f.sub.b3D = f.sub.min IF (21) = true with X = I, THEN compute R according to (15) and v according to (24).

[0197] The condition used to determine which case is present is defined as (21):

[00018]$\begin{array}{cc}.Math.v_{1X}-v_{2X}.Math..Math.\frac{2}{3}.Math.\left(\frac{1}{{\lambda }_1}+\frac{1}{{\lambda }_2}+\frac{1}{{\lambda }_3}\right)<\frac{1}{T}\mathrm{AND}.Math.v_{1X}-v_{3X}.Math..Math.\frac{2}{3}.Math.\left(\frac{1}{{\lambda }_1}+\frac{1}{{\lambda }_2}+\frac{1}{{\lambda }_3}\right)<\frac{1}{T}\mathrm{AND}.Math.v_{2X}-v_{3X}.Math..Math.\frac{2}{3}.Math.\left(\frac{1}{{\lambda }_1}+\frac{1}{{\lambda }_2}+\frac{1}{{\lambda }_3}\right)<\frac{1}{T}& \left(21\right)\end{array}$

[0198] This equation mathematically expresses that the preliminary velocity values obtained using the light produces by the first, second and third light source are all very similar.

[0199] For the computation of the final range R, Table 2 refers to equations (13) to (15) recited above, and for the computation of the final velocity v to the following equations (22) to (24):

[00019]$\begin{array}{cc}v={\lambda }_1\left(\frac{1}{2}{f}_{b1D}-\frac{R.Math.r_{chirp1}}{c}\right)& \left(22\right)\end{array}$$\begin{array}{cc}v=-{\lambda }_1.Math.-\frac{1}{2}{f}_{b1D}+\frac{R.Math.r_{chirp1}}{c}.Math.& \left(23\right)\end{array}$$\begin{array}{cc}v=-{\lambda }_1.Math.\frac{1}{2}{f}_{b1D}+\frac{R.Math.r_{chirp1}}{c}.Math.& \left(24\right)\end{array}$

[0200] Cases A to I do not occur with equal frequency. Rather, case A dominates so that the decision tree does not have to be run through for most measurements. This significantly reduces the computational effort required to perform the decision tree analysis.

[0201] For typical LiDAR applications such as autonomously driving vehicles, a range error tolerance of 5 cm and a velocity tolerance of 5 cm/s is acceptable. On the basis of Monte

[0202] Carlo simulations, it can be shown that for such tolerances more than 99.9% of the computations fulfill these tolerances. This percentage can be increased still further by increasing the chirp rate of the light sources 16, 17, 18 and/or the duration T of the measuring interval.

[0203] A part of the few remaining mismeasurements can be identified by the fact that the computed values cannot occur in reality. For example, negative values for the range R or velocities |v|>70 m/s cannot occur in road traffic situations. The rest of the mismeasurements can be identified by plausibility analyses of temporally or spatially neighboring pixels.

[0204] FIG. 18 is a diagram showing how the signal-to-noise ratio (SNR) improves with the use of two or three light sources compared to the conventional prior art method. For the comparison it was assumed that both the total light source power and the total measurement duration are the same for all measurements. According to the invention, measurements are made with a longer chirp duration Tin order to achieve the same beat peak height. In return, a significant time saving is achieved, because only one (in the case of three light sources) or two (in the case of two light sources) FFTs instead of three FFTs have to be performed sequentially.

[0205] In the comparison, a maximum object distance of 300 m has been assumed, which corresponds to a time of flight ToF=2 μs. This unavoidable runtime is the reason why the approach according to the invention results in an improved SNR compared to the prior art approach. This is most easily understood by assuming three sequential FFT measurements each having a duration T=2 μs for the conventional measurement. Even with infinitely high laser power, no signal photons would be obtained from an object at a distance of 300 m, since the measurements are finished before the photons arrive at the scanner system. If, on the other hand, one performs an FFT measurement with 6 μs (=3.Math.T), signal photons reach the measurement system for the duration 3.Math.T-ToF=4 μs. It is therefore advantageous with regard to the SNR to perform fewer FFT measurements with longer chirp durations T. One of the advantages of this embodiment is that the dead time, in which no photons are yet measurable, is significantly reduced, since the dead time occurs with every FFT.

[0206] The diagram in FIG. 18 shows how many times the SNR improves, according to the invention, as a function of the chirp duration T. It can be seen that an FFT measurement with three simultaneously operating light sources, i.e. with three 3 beat peaks in one FFT, always offers the greatest SNR advantage. This must be traded off against the disadvantage of higher system complexity due to the provision of a third light source and the electronics that are required for its control.