Optical detection of range using lasers, often referenced by a mnemonic, LIDAR, for light detection and ranging, also sometimes called laser RADAR (radio-wave detection and ranging), is used for a variety of applications, from altimetry, to imaging, to collision avoidance. LIDAR provides finer scale range resolution with smaller beam sizes than conventional microwave ranging systems, such as RADAR. Optical detection of range can be accomplished with several different techniques, including direct ranging based on round trip travel time of an optical pulse to an object, and chirped detection based on a frequency difference between a transmitted chirped optical signal and a returned signal scattered from an object, and phase-encoded detection based on a sequence of single frequency phase changes that are distinguishable from natural signals.
To achieve acceptable range accuracy and detection sensitivity, direct long-range LIDAR systems use short pulse lasers with low pulse repetition rate and extremely high pulse peak power. The high pulse power can lead to rapid degradation of optical components. Chirped and phase-encoded LIDAR systems use long optical pulses with relatively low peak optical power. In this configuration, the range accuracy increases with the chirp bandwidth or length and bandwidth of the phase codes rather than the pulse duration, and therefore excellent range accuracy can still be obtained.
Useful optical bandwidths have been achieved using wideband radio frequency (RF) electrical signals to modulate an optical carrier. Recent advances in LIDAR include using the same modulated optical carrier as a reference signal that is combined with the returned signal at an optical detector to produce in the resulting electrical signal a relatively low beat frequency in the RF band that is proportional to the difference in frequencies or phases between the references and returned optical signals. This kind of beat frequency detection of frequency differences at a detector is called heterodyne detection. It has several advantages known in the art, such as the advantage of using RF components of ready and inexpensive availability.
Recent work by current inventors, show a novel arrangement of optical components and coherent processing to detect Doppler shifts in returned signals that provide not only improved range but also relative signed speed on a vector between the LIDAR system and each external object. These systems are called hi-res range-Doppler LIDAR herein. See for example World Intellectual Property Organization (WIPO) publications WO 2018/160240 and WO 2018/144853.
These improvements provide range, with or without target speed, in a pencil thin laser beam of proper frequency or phase content. When such beams are swept over a scene, information about the location and speed of surrounding objects can be obtained. This information is expected to be of value in control systems for autonomous vehicles, such as self driving, or driver assisted, automobiles.
The sampling and processing that provides range accuracy and target speed accuracy involve integration of one or more laser signals of various durations, in a time interval called integration time. To cover a scene in a timely way involves repeating a measurement of sufficient accuracy (involving one or more signals often over one to tens of microseconds) often enough to sample a variety of angles (often on the order of thousands) around the autonomous vehicle to understand the environment around the vehicle before the vehicle advances too far into the space ahead of the vehicle (a distance on the order of one to tens of meters, often covered in a particular time on the order of one to a few seconds). The number of different angles that can be covered in the particular time (often called the cycle or sampling time) depends on the sampling rate. The current inventors have recognized that a tradeoff can be made between integration time for range and speed accuracy, sampling rate, and pattern of sampling different angles, with one or more LIDAR beams, to effectively determine the environment in the vicinity of an autonomous vehicle as the vehicle moves through that environment.
In a first set of embodiments, a method for optimizing a scan pattern of a LIDAR system on an autonomous vehicle includes receiving, on a processor, data that indicates first signal-to-noise ratio (SNR) values of a signal reflected by a target and detected by the LIDAR system based on values of a range of the target, where the first SNR values are for a respective value of a scan rate of the LIDAR system. The method further includes receiving, on the processor, data that indicates second SNR values of the signal based on values of the range of the target, where the second SNR values are for a respective value of an integration time of the LIDAR system. The method further includes receiving, on the processor, data that indicates a first angle and a second angle that defines an angle range of the scan pattern. The method further includes receiving, on the processor, data that indicates a maximum design range of the target at each angle in the angle range. The method further includes determining, for each angle in the angle range, a maximum scan rate of the LIDAR system based on a maximum value among of those scan rates where the first SNR value based on the maximum design range exceeds a minimum SNR threshold. The method further includes determining, for each angle in the angle range, a minimum integration time of the LIDAR system based on a minimum value among of those integration times where the second SNR value based on the maximum design range exceeds a minimum SNR threshold. The method further includes defining, with the processor, a scan pattern of the LIDAR system based on the maximum scan rate and the minimum integration time at each angle in the angle range. The method further includes operating the LIDAR system according to the scan pattern.
In other embodiments, a system or apparatus or computer-readable medium is configured to perform one or more steps of the above methods.
Still other aspects, features, and advantages are readily apparent from the following detailed description, simply by illustrating a number of particular embodiments and implementations, including the best mode contemplated for carrying out the invention. Other embodiments are also capable of other and different features and advantages, and their several details can be modified in various obvious respects, all without departing from the spirit and scope of the invention. Accordingly, the drawings and description are to be regarded as illustrative in nature, and not as restrictive.
Embodiments are illustrated by way of example, and not by way of limitation, in the figures of the accompanying drawings in which like reference numerals refer to similar elements and in which:
A method and apparatus and system and computer-readable medium are described for scanning of LIDAR to support operation of a vehicle. In the following description, for the purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the present invention. It will be apparent, however, to one skilled in the art that the present invention may be practiced without these specific details. In other instances, well-known structures and devices are shown in block diagram form in order to avoid unnecessarily obscuring the present invention.
Notwithstanding that the numerical ranges and parameters setting forth the broad scope are approximations, the numerical values set forth in specific non-limiting examples are reported as precisely as possible. Any numerical value, however, inherently contains certain errors necessarily resulting from the standard deviation found in their respective testing measurements at the time of this writing. Furthermore, unless otherwise clear from the context, a numerical value presented herein has an implied precision given by the least significant digit. Thus a value 1.1 implies a value from 1.05 to 1.15. The term “about” is used to indicate a broader range centered on the given value, and unless otherwise clear from the context implies a broader range around the least significant digit, such as “about 1.1” implies a range from 1.0 to 1.2. If the least significant digit is unclear, then the term “about” implies a factor of two, e.g., “about X” implies a value in the range from 0.5× to 2×, for example, about 100 implies a value in a range from 50 to 200. Moreover, all ranges disclosed herein are to be understood to encompass any and all sub-ranges subsumed therein. For example, a range of “less than 10” for a positive only parameter can include any and all sub-ranges between (and including) the minimum value of zero and the maximum value of 10, that is, any and all sub-ranges having a minimum value of equal to or greater than zero and a maximum value of equal to or less than 10, e.g., 1 to 4.
Some embodiments of the invention are described below in the context of a single front mounted hi-res Doppler LIDAR system on a personal automobile; but, embodiments are not limited to this context. In other embodiments, one or multiple systems of the same type or other high resolution LIDAR, with or without Doppler components, with overlapping or non-overlapping fields of view or one or more such systems mounted on smaller or larger land, sea, air or space vehicles, piloted or autonomous, are employed.
Using an optical phase-encoded signal for measurement of range, the transmitted signal is in phase with a carrier (phase=0) for part of the transmitted signal and then changes by one or more phase changes represented by the symbol Δϕ (so phase=0, Δϕ, 2Δϕ . . . ) for short time intervals, switching back and forth between the two or more phase values repeatedly over the transmitted signal. The shortest interval of constant phase is a parameter of the encoding called pulse duration τ and is typically the duration of several periods of the lowest frequency in the band. The reciprocal, 1/τ, is baud rate, where each baud indicates a symbol. The number N of such constant phase pulses during the time of the transmitted signal is the number N of symbols and represents the length of the encoding. In binary encoding, there are two phase values and the phase of the shortest interval can be considered a 0 for one value and a 1 for the other, thus the symbol is one bit, and the baud rate is also called the bit rate. In multiphase encoding, there are multiple phase values. For example, 4 phase values such as Δϕ*{0, 1, 2 and 3}, which, for Δϕ=π/2 (90 degrees), equals {0, π/2, π and 3π/2}, respectively; and, thus 4 phase values can represent 0, 1, 2, 3, respectively. In this example, each symbol is two bits and the bit rate is twice the baud rate.
Phase-shift keying (PSK) refers to a digital modulation scheme that conveys data by changing (modulating) the phase of a reference signal (the carrier wave). The modulation is impressed by varying the sine and cosine inputs at a precise time. At radio frequencies (RF), PSK is widely used for wireless local area networks (LANs), RF identification (RFID) and Bluetooth communication. Alternatively, instead of operating with respect to a constant reference wave, the transmission can operate with respect to itself changes in phase of a single transmitted waveform can be considered the symbol. In this system, the demodulator determines the changes in the phase of the received signal rather than the phase (relative to a reference wave) itself. Since this scheme depends on the difference between successive phases, it is termed differential phase-shift keying (DPSK). DPSK can be significantly simpler to implement in communications applications than ordinary PSK, since there is no need for the demodulator to have a copy of the reference signal to determine the exact phase of the received signal (thus, it is a non-coherent scheme).
For optical ranging applications, since the transmitter and receiver are in the same device, coherent PSK can be used. The carrier frequency is an optical frequency fc and a RF f0 is modulated onto the optical carrier. The number N and duration τ of symbols are selected to achieve the desired range accuracy and resolution. The pattern of symbols is selected to be distinguishable from other sources of coded signals and noise. Thus, a strong correlation between the transmitted and returned signal is a strong indication of a reflected or backscattered signal. The transmitted signal is made up of one or more blocks of symbols, where each block is sufficiently long to provide strong correlation with a reflected or backscattered return even in the presence of noise. In the following discussion, it is assumed that the transmitted signal is made up of M blocks of N symbols per block, where M and N are non-negative integers.
The observed frequency f′ of the return differs from the correct frequency f=fc+f0 of the return by the Doppler effect given by Equation 1.
Where c is the speed of light in the medium, vo is the velocity of the observer and vs is the velocity of the source along the vector connecting source to receiver. Note that the two frequencies are the same if the observer and source are moving at the same speed in the same direction on the vector between the two. The difference between the two frequencies, Δf=f′−f, is the Doppler shift, ΔfD, which causes problems for the range measurement, and is given by Equation 2.
Note that the magnitude of the error increases with the frequency f of the signal. Note also that for a stationary LIDAR system (vo=0), for an object moving at 10 meters a second (vs=10), and visible light of frequency about 500 THz, then the size of the Doppler shift is on the order of 16 megahertz (MHz, 1 MHz=106 hertz, Hz, 1 Hz=1 cycle per second). In various embodiments described below, the Doppler shift is detected and used to process the data for the calculation of range.
In phase coded ranging, the arrival of the phase coded return is detected in the return signal by cross correlating the transmitted signal or other reference signal with the returned signal, implemented practically by cross correlating the code for a RF signal with a electrical signal from an optical detector using heterodyne detection and thus down-mixing back to the RF band. Cross correlation for any one lag is computed by convolving the two traces, i.e., multiplying corresponding values in the two traces and summing over all points in the trace, and then repeating for each time lag. Alternatively, the cross correlation can be accomplished by a multiplication of the Fourier transforms of each of the two traces followed by an inverse Fourier transform. Efficient hardware and software implementations for a Fast Fourier transform (FFT) are widely available for both forward and inverse Fourier transforms.
Note that the cross correlation computation is typically done with analog or digital electrical signals after the amplitude and phase of the return is detected at an optical detector. To move the signal at the optical detector to a RF frequency range that can be digitized easily, the optical return signal is optically mixed with the reference signal before impinging on the detector. A copy of the phase-encoded transmitted optical signal can be used as the reference signal, but it is also possible, and often preferable, to use the continuous wave carrier frequency optical signal output by the laser as the reference signal and capture both the amplitude and phase of the electrical signal output by the detector.
For an idealized (noiseless) return signal that is reflected from an object that is not moving (and thus the return is not Doppler shifted), a peak occurs at a time Δt after the start of the transmitted signal. This indicates that the returned signal includes a version of the transmitted phase code beginning at the time Δt. The range R to the reflecting (or backscattering) object is computed from the two way travel time delay based on the speed of light c in the medium, as given by Equation 3.
R=c*Δt/2 (3)
For an idealized (noiseless) return signal that is scattered from an object that is moving (and thus the return is Doppler shifted), the return signal does not include the phase encoding in the proper frequency bin, the correlation stays low for all time lags, and a peak is not as readily detected, and is often undetectable in the presence of noise. Thus Δt is not as readily determined; and, range R is not as readily produced.
According to various embodiments of the inventor's previous work, the Doppler shift is determined in the electrical processing of the returned signal; and the Doppler shift is used to correct the cross-correlation calculation. Thus, a peak is more readily found; and, range can be more readily determined.
In some Doppler compensation embodiments, rather than finding ΔfD by taking the spectrum of both transmitted and returned signals and searching for peaks in each, then subtracting the frequencies of corresponding peaks, as illustrated in
As described in more detail in inventor's previous work the Doppler shift(s) detected in the cross spectrum are used to correct the cross correlation so that the peak 135 is apparent in the Doppler compensated Doppler shifted return at lag Δt, and range R can be determined. In some embodiments simultaneous I/Q processing is performed as described in more detail in World Intellectual Property Organization publication WO 2018/144853 entitled “Method and system for Doppler detection and Doppler correction of optical phase-encoded range detection”, the entire contents of which are hereby incorporated by reference as if fully set forth herein. In other embodiments, serial I/Q processing is used to determine the sign of the Doppler return as described in more detail in World Intellectual Property Organization publication WO 2019/014177 entitled “Method and System for Time Separated Quadrature Detection of Doppler Effects in Optical Range Measurements”, the entire contents of which are herby incorporated by reference as if fully set forth herein. In other embodiments, other means are used to determine the Doppler correction; and, in various embodiments, any method known in the art to perform Doppler correction is used. In some embodiments, errors due to Doppler shifting are tolerated or ignored; and, no Doppler correction is applied to the range measurements.
The returned signal is depicted in graph 160 which has a horizontal axis 102 that indicates time and a vertical axis 114 that indicates frequency as in graph 110. The chirp 116 of graph 110 is also plotted as a dotted line on graph 160. A first returned signal is given by trace 166a, which is just the transmitted reference signal diminished in intensity (not shown) and delayed by Δt. When the returned signal is received from an external object after covering a distance of 2R, where R is the range to the target, the returned signal start at the delayed time Δt is given by 2R/c, where c is the speed of light in the medium (approximately 3×108 meters per second, m/s), related according to Equation 3, described above. Over this time, the frequency has changed by an amount that depends on the range, called fR, and given by the frequency rate of change multiplied by the delay time. This is given by Equation 4a.
fR=(f2−f1)/τ*2R/c=2BR/cτ (4a)
The value of fR is measured by the frequency difference between the transmitted signal 116 and returned signal 166a in a time domain mixing operation referred to as de-chirping. So the range R is given by Equation 4b.
R=fRcτ/2B (4b)
Of course, if the returned signal arrives after the pulse is completely transmitted, that is, if 2R/c is greater than τ, then Equations 4a and 4b are not valid. In this case, the reference signal is delayed a known or fixed amount to ensure the returned signal overlaps the reference signal. The fixed or known delay time of the reference signal is multiplied by the speed of light, c, to give an additional range that is added to range computed from Equation 4b. While the absolute range may be off due to uncertainty of the speed of light in the medium, this is a near-constant error and the relative ranges based on the frequency difference are still very precise.
In some circumstances, a spot (pencil beam cross section) illuminated by the transmitted light beam encounters two or more different scatterers at different ranges, such as a front and a back of a semitransparent object, or the closer and farther portions of an object at varying distances from the LIDAR, or two separate objects within the illuminated spot. In such circumstances, a second diminished intensity and differently delayed signal will also be received, indicated on graph 160 by trace 166b. This will have a different measured value off that gives a different range using Equation 4b. In some circumstances, multiple additional returned signals are received.
Graph 170 depicts the difference frequency fR between a first returned signal 166a and the reference chirp 116. The horizontal axis 102 indicates time as in all the other aligned graphs in
A common method for de-chirping is to direct both the reference optical signal and the returned optical signal to the same optical detector. The electrical output of the detector is dominated by a beat frequency that is equal to, or otherwise depends on, the difference in the frequencies of the two signals converging on the detector. A Fourier transform of this electrical output signal will yield a peak at the beat frequency. This beat frequency is in the radio frequency (RF) range of Megahertz (MHz, 1 MHz=106 Hertz=106 cycles per second) rather than in the optical frequency range of Terahertz (THz, 1 THz=1012 Hertz). Such signals are readily processed by common and inexpensive RF components, such as a Fast Fourier Transform (FFT) algorithm running on a microprocessor or a specially built FFT or other digital signal processing (DSP) integrated circuit. In other embodiments, the return signal is mixed with a continuous wave (CW) tone acting as the local oscillator (versus a chirp as the local oscillator). This leads to the detected signal which itself is a chirp (or whatever waveform was transmitted). In this case the detected signal would undergo matched filtering in the digital domain as described in Kachelmyer 1990, the entire contents of which are hereby incorporated by reference as if fully set forth herein, except for terminology inconsistent with that used herein. The disadvantage is that the digitizer bandwidth requirement is generally higher. The positive aspects of coherent detection are otherwise retained.
In some embodiments, the LIDAR system is changed to produce simultaneous up and down chirps. This approach eliminates variability introduced by object speed differences, or LIDAR position changes relative to the object which actually does change the range, or transient scatterers in the beam, among others, or some combination. The approach then guarantees that the Doppler shifts and ranges measured on the up and down chirps are indeed identical and can be most usefully combined. The Doppler scheme guarantees parallel capture of asymmetrically shifted return pairs in frequency space for a high probability of correct compensation.
In some embodiments, two different laser sources are used to produce the two different optical frequencies in each beam at each time. However, in some embodiments, a single optical carrier is modulated by a single RF chirp to produce symmetrical sidebands that serve as the simultaneous up and down chirps. In some of these embodiments, a double sideband Mach-Zehnder intensity modulator is used that, in general, does not leave much energy in the carrier frequency; instead, almost all of the energy goes into the sidebands.
As a result of sideband symmetry, the bandwidth of the two optical chirps will be the same if the same order sideband is used. In other embodiments, other sidebands are used, e.g., two second order sideband are used, or a first order sideband and a non-overlapping second sideband is used, or some other combination.
As described in World Intellectual Property Organization publication WO 2018/160240 entitled “Method and System for Doppler Detection and Doppler Correction of Optical Chirped Range Detection,” the entire contents of which are hereby incorporated by reference as if fully set forth herein, when selecting the transmit (TX) and local oscillator (LO) chirp waveforms, it is advantageous to ensure that the frequency shifted bands of the system take maximum advantage of available digitizer bandwidth. In general, this is accomplished by shifting either the up chirp or the down chirp to have a range frequency beat close to zero.
In order to depict how to use hi-res range-Doppler detection systems, some generic hardware approaches are described.
The transmitted signal is then transmitted to illuminate an area of interest, often through some scanning optics 218. The detector array is a single paired or unpaired detector or a one-dimensional (1D) or two-dimensional (2D) array of paired or unpaired detectors arranged in a plane roughly perpendicular to returned beams 291 from the object. The reference beam 207b and returned beam 291 are combined in zero or more optical mixers 284 to produce an optical signal of characteristics to be properly detected. The frequency, phase or amplitude of the interference pattern, or some combination, is recorded by acquisition system 240 for each detector at multiple times during the signal duration D. The number of temporal samples processed per signal duration or integration time affects the down-range extent. The number or integration time is often a practical consideration chosen based on number of symbols per signal, signal repetition rate and available camera frame rate. The frame rate is the sampling bandwidth, often called “digitizer frequency.” The only fundamental limitations of range extent are the coherence length of the laser and the length of the chirp or unique phase code before it repeats (for unambiguous ranging). This is enabled because any digital record of the returned heterodyne signal or bits could be compared or cross correlated with any portion of transmitted bits from the prior transmission history.
The acquired data is made available to a processing system 250, such as a computer system described below with reference to
Any known apparatus or system may be used to implement the laser source 212, modulators 282a, 282b, beam splitter 216, reference path 220, optical mixers 284, detector array 230, scanning optics 218, or acquisition system 240. Optical coupling to flood or focus on a target or focus past the pupil plane are not depicted. As used herein, an optical coupler is any component that affects the propagation of light within spatial coordinates to direct light from one component to another component, such as a vacuum, air, glass, crystal, mirror, lens, optical circulator, beam splitter, phase plate, polarizer, optical fiber, optical mixer, among others, alone or in some combination.
In an embodiment, monostatic coherent LIDAR performance of the system 200′ is modeled by including system parameters in a so called “link budget” and a LIDAR is operated based on that modeling. A link budget estimates the expected value of the signal to noise ratio (SNR) for various system and target parameters. In one embodiment, on the system side, a link budget includes one or more of output optical power, integration time, detector characteristics, insertion losses in waveguide connections, mode overlap between the imaged spot and the monostatic collection waveguide, and optical transceiver characteristics. In another embodiment, on the target side, a link budget includes one or more of atmospheric characteristics, target reflectivity, and target range.
In the near field 406, a primary driver of the SNR is a diameter of the collimated return beam 291 before it is focused by the collimation optics 229 to the tip 217.
In an embodiment, in the near field 406, as the diameter of the collimated return beam 291 grows at larger target ranges, a diameter of the focused return beam 291 by the collimation optics 229 at the tip 217 shrinks.
While the discussion in relation to
Based on the traces in
Ultimately, beam walk off 419 is a significant inhibitor of SNR in coherent monostatic LIDAR system 200′. In conventional coherent monostatic LIDAR systems, the scan rate of the beam is set at a fixed scan rate over the angle range and the resulting target range, where the fixed scan rate is chosen so that the associated SNR of the fixed scan rate is above an SNR threshold for the entire target range of interest. In conventional coherent LIDAR systems, this results in a relatively low fixed scan rate being used to scan the beam over a scan trajectory 460, which results in large gaps 462 between adjacent scans, as depicted in
In addition to the scan rate of the beam, the SNR of the return beam 291 is affected by the integration time over which the acquisition system 240 and/or processing system 250 samples and processes the return beam 291. In some embodiments, the beam is scanned between discrete angles and is held stationary or almost stationary at discrete angles in the angle range 227 for a respective integration time at each discrete angle. The SNR of the return beam 291 is affected by the value of the integration time and the target range. As previously discussed, the cross-sectional area of the beam increases with target range resulting in increased atmospheric scattering and thus an intensity of the return beam 291 decreases with increasing range at the rate of 1/R2. Accordingly, a longer integration time is needed to achieve the same SNR for a return beam 291 from a longer target range.
Conventional monostatic LIDAR systems typically select a fixed integration time (e.g., 1.6 μs) for the scanning at the range of angles 227 and resulting target ranges, so that the SNR associated with the fixed integration time exceeds a target SNR threshold indicated by horizontal dashed line 452 over the target range. It is here recognized that this results in relatively high integration times being used to scan the beam over the range of angles 227 in conventional monostatic LIDAR systems, integration times that are higher than needed at some angles.
Here is presented a method for operating a LIDAR that addresses this noted drawback, including minimizing the integration time at each angle within the range of angles 227 using the target range at each angle, so to minimize the integration time over at least some angles 227.
In some embodiments a vehicle is controlled at least in part based on data received from a hi-res Doppler LIDAR system mounted on the vehicle when the LIDAR is operated as described herein.
In some embodiments, the vehicle includes ancillary sensors (not shown), such as a GPS sensor, odometer, tachometer, temperature sensor, vacuum sensor, electrical voltage or current sensors, among others well known in the art. In some embodiments, a gyroscope 330 is included to provide rotation information.
In designing the system 301′, a predetermined maximum design range of the beams at each angle is determined and represents a maximum anticipated target range at each angle in the range 326. In other embodiments, the maximum design range at each angle is not predetermined but is regularly measured and updated in the memory of the processing system 250 at incremental time periods. In an embodiment, the first beam 342 is oriented toward the surface 349 (inclined downward) and intersects the surface 349 within some maximum design range from the vehicle 310. Thus, at the first angle, the system 320 does not consider targets positioned below the surface 349. In an example embodiment, the first angle of the first beam 342 is about −15 degrees or in a range from about −25 degrees to about −10 degrees with respect to the arrow 313; and the maximum design range is about 4 meters (m) or within a range from about 1 m to about 10 m, or in a range from about 2 m to about 6 m. In an embodiment, the second beam 346 is inclined upward (e.g., oriented toward the sky) and intersects a ceiling 347 above which the operation of vehicle 310 is not concerned, within some maximum design range from the vehicle 310. Thus, at the second angle the system 320 does not consider targets positioned above the ceiling 347. In an example embodiment, the ceiling 347 is at an altitude of about 12 m or in a range from about 8 m to about 15 m from the surface 349 (e.g., that defines an altitude of 0 m), the second angle of the second beam 346 is about 15 degrees or in a range from about 10 degrees to about 20 degrees with respect to the arrow 313; and the maximum design range is about 7 m or within a range from about 4 m to about 10 m or within a range from about 1 m to about 15 m. In some embodiments, the ceiling 347 altitude depends on an altitude of the LIDAR system 320 (e.g. about 1.5 m or in a range of about 1 m to about 4 m, where the surface 349 defines 0 m). In an embodiment, an intermediate beam 344 between the first beam 342 and second beam 346 is oriented about parallel with the arrow 313 and intersects a target 343 positioned at a maximum design range from the vehicle 310.
Although
In step 601, data is received on a processor that indicates first SNR values of a signal reflected by a target and detected by the LIDAR system based on values of a range of the target, where the first SNR values are for a respective value of a scan rate of the LIDAR system. In an embodiment, in step 601 the data is first SNR values of the focused return beam 291 on the fiber tip 217 in the system 200′. In another embodiment, in step 601 the data is first SNR values of the focused return beam 291 on the detector array 230 in the system 200. In one embodiment, the data includes values of traces 440a and/or trace 440b and/or trace 440c, or the equations or models that produce those traces, and which indicate SNR values of the return beam 291, where each trace 440 is for a respective value of the scan rate of the beam. In some embodiments, the data is not limited to traces 440a, 440b, 440c and includes SNR values of less or more curves than are depicted in
In some embodiments, the first SNR values are predetermined and are received by the processor in step 601. In other embodiments, the first SNR values are measured by the LIDAR system and subsequently received by the processor in step 601. In one embodiment, the data is input in step 601 using an input device 712 and/or uploaded to the memory 704 of the processing system 250 over a network link 778 from a local area network 780, internet 790 or external server 792, either unsolicited or in response to a query.
In step 603, data is received on a processor that indicates second SNR values of a signal reflected by a target and detected by the LIDAR system based on values of a range of the target, where the second SNR values are for a respective value of an integration time of the LIDAR system. In an embodiment, in step 603 the data is second SNR values of the focused return beam 291 in the system 200′ for a respective integration time over which the beam is processed by the acquisition system 240 and/or processing system 250. In one embodiment, the data includes values of trace 450a and/or trace 450b and/or trace 450c and/or trace 450d or the equations or models that produce those traces, and which indicate SNR values of the return beam 291, where each trace 450 is for a respective value of the integration time that the beam is processed by the acquisition system 240 and/or processing system 250. In some embodiments, the data is not limited to traces 450a, 450b, 450c, 450d and includes less or more traces than are depicted in
In some embodiments, the second SNR values are predetermined and are received by the processor in step 603. In other embodiments, the second SNR values are measured by the LIDAR system and subsequently received by the processor in step 603. In one embodiment, the data is input in step 603 using an input device 712 and/or uploaded to the memory 704 of the processing system 250 over a network link 778 from a local area network 780, internet 790 or external server 792, either unsolicited or in response to a query.
In step 605, data is received on a processor that indicates the first angle and the second angle that defines the angle range 326. In one embodiment, in step 605 the first angle and the second angle are input using an input device 712 (e.g. mouse or pointing device 716) and/or are uploaded to the processing system 250 over a network link 778. In an embodiment the first angle is defined as the angle between the first beam 342 and the arrow 313 indicating the direction of travel of the vehicle 310 and the second angle is defined as the angle between the second beam 346 and the arrow 313. In an embodiment, the first angle and second angle are symmetric with respect to the arrow 313, e.g. the first angle and the second angle are equal and opposite to each other. In one embodiment, the first angle is selected so that the first beam 342 is oriented towards the surface 349 and the second angle is selected so that the second beam 346 is oriented away from the surface 349 and towards the ceiling 347, in one or two dimensions. In one embodiment, steps 601, 603 and 605 are simultaneously performed in one step where the data in steps 601, 603 and 605 is received at the processor in one simultaneously step.
In step 607, data is received on a processor that indicates the maximum design range of the target at each angle in the angle range. In an embodiment, the maximum design range is the predetermined maximum range of the target at each angle in the range of angles 326. In one embodiment, in step 607 the maximum design range of the target at each angle is provided by the traces 512a, 512b. The maximum design range data received in step 607 is not limited to the data of the traces 512a, 512b which is based on a specific range of angles 326. In one embodiment, traces 512a, 512b depict that the maximum design range for beams 342, 346 at the first angle and the second angle (e.g. <100 m) is less than the maximum design range of intermediate beam 344 at an angle between the first and second angle (e.g. >100 m). In another embodiment, the data in step 607 is provided over a first range of angles that is greater than the range of angles 326. In one embodiment, the data in step 607 is provided at an incremental angle over the angle range, wherein the incremental angle is selected in a range from about 0.005 degrees to about 0.01 degrees or in a range from about 0.0005 degrees to about 0.01 degrees.
In one example embodiment, the data in step 607 is input using an input device 712 (e.g. mouse or pointing device 716) and/or are uploaded to the processing system 250 over a network link 778. In some embodiments, the maximum design range is predetermined and received during step 607. In other embodiments, the system 200, 200′ is used to measure the maximum design range at each angle in the angle range 326 and the maximum design range at each angle is subsequently received by the processing system 250 in step 607.
In step 609, a maximum scan rate of the LIDAR system is determined at each angle in the angle range 326 so that the SNR of the LIDAR system is greater than a minimum SNR threshold. At each angle, the maximum design range for that angle is first determined based on the received data in step 607. In one embodiment, the maximum design range for the angle is determined using trace 512a or 512b. First SNR values received in step 601 are then determined for the maximum design range at the angle and it is further determined which of these first SNR values exceed the minimum SNR threshold. In one embodiment, values of traces 440a, 440b, 440c are determined for a maximum design range (e.g. about 120 m) and it is further determined that the values of traces 440a, 440b exceeds the minimum SNR threshold 442. Among those first SNR values which exceed the minimum SNR threshold, the first SNR values with the maximum scan rate is selected and the maximum scan rate is determined in step 609 for that angle. In the above embodiment, among the values of the traces 440a, 440b which exceeds the minimum SNR threshold 442 at the maximum design range (e.g. about 120 m), the curve 440b values are selected as the maximum scan rate and the maximum scan rate (e.g. moderate scan rate associated with curve 440b) is determined in step 609 for that angle. In an embodiment,
In step 611, a minimum integration time of the LIDAR system is determined at each angle in the angle range 326 so that the SNR of the LIDAR system is greater than a minimum SNR threshold. At each angle, the maximum design range for that angle is first determined based on the received data in step 607. In one embodiment, the maximum design range for the angle is determined using trace 512a or 512b. Second SNR values received in step 603 are then determined for the maximum design range at the angle and it is further determined which of these second SNR values exceed the minimum SNR threshold. In one embodiment, values of traces 450a, 450b, 450c, 450d are determined for a maximum design range (e.g. about 120 m) and it is further determined that the values of traces 450a, 450b, 450c exceeds the minimum SNR threshold 452. Among those second SNR values which exceed the minimum SNR threshold, the second SNR values with the minimum integration time is selected and the minimum integration time is determined in step 611 for that angle. In the above embodiment, among the values of the curves 450a, 450b, 450c which exceeds the minimum SNR threshold 452 at the maximum design range (e.g. about 120 m), the curve 450c values are selected with the minimum integration time and the minimum integration time (e.g. about 800 ns) is determined in step 611 for that angle.
In step 611, the minimum integration time of the LIDAR system 320 is determined at the angle 345a. The maximum design range (e.g. 30 m) of the beam 343a at the angle 345a is first determined using the data in step 607. Second SNR values from step 603 for the maximum design range are then determined. In an example embodiment, the second SNR values include values of traces 450a, 450b, 450c, 450d. It is then determined which of those second SNR values at the maximum design range exceeds the minimum SNR threshold. In the example embodiment, the values of the traces 450a, 450b, 450c, 450d at the maximum design range (e.g. 30 m) all exceed the minimum SNR threshold 452. It is then determined which of these second SNR values has the minimum integration time and this minimum integration time is determined in step 611 for that angle. In the example embodiment, the curve 450d has the minimum integration time (e.g. about 400 ns) and thus this minimum integration time is used to process the beam 343a at the angle 345a.
In step 613, a determination is made whether additional angles remain in the angle range 326 to perform another iteration of steps 609, 611. In one embodiment, where an initial iteration of steps 609, 611 was at the first angle of the first beam 342 in the range 326, step 613 involves a determination of whether the previous iteration of steps 609, 611 was at or beyond the second angle of the second beam 346 in the angle range 326. In another embodiment, where an initial iteration of the steps 609, 611 was at the second angle of the range 326, step 613 involves a determination of whether the previous iteration of steps 609, 611 was at or beyond the first angle of the angle range 326. If step 613 determines that more angles in the range 326 remain, the method proceeds to block 615. If step 613 determines that no further angles in the range 326 remain, the method proceeds to block 617.
In step 615, a determination is made of a subsequent angle at which to iterate steps 609, 611 after having performed a previous iteration of steps 609, 611 at an initial angle. In an embodiment,
θs=θi+SmIm (5)
In an example embodiment, if the initial angle 345a is −15 degrees, the maximum scan rate is 15 degrees per second and the minimum integration time is 2 μs, then the subsequent angle 345b is about −14.97 degrees using Equation 5. After determining the subsequent angle in step 615, the method proceeds back to block 609 so that steps 609, 611 are iterated at the subsequent angle.
In step 617, after it is determined that no further iterations of steps 609, 611 need to be performed, a scan pattern of the LIDAR system is defined based on the maximum scan rate from step 609 and the minimum integration time from step 611 at each angle in the angle range 326. In one embodiment, the scan pattern includes the maximum scan pattern and minimum integration time for each angle in the angle range 326 between the first angle of the first beam 342 and the second angle of the second beam 346. In an example embodiment, the scan pattern is stored in a memory (e.g. memory 704) of the processing system 250. In another example embodiment, the angle increments between adjacent angles in the scan pattern is determined in step 615, e.g. the angle increment is the spacing between the subsequent angle and the initial angle in step 615. In another example embodiment,
In step 619, the LIDAR system is operated according to the scan pattern determined in step 617. In an embodiment, in step 619 the beam of the LIDAR system is scanned in the angle range 326 over one or more cycles, where the scan rate of the beam at each angle is the maximum scan rate in the scan pattern for that angle and the integration time of the LIDAR system at each angle is the minimum integration time for that angle. In one embodiment, in step 619, at each angle the processing system 250 of the LIDAR system transmits one or more signals to the scanning optics 218 so that the scan rate at each angle is the maximum scan rate of the scan pattern for that angle. Additionally, in one embodiment, in step 619, at each angle the processing system 250 of the LIDAR system adjusts the integration time of the acquisition system 240 and/or processing system 250 for the return beam 291 received at each angle so that the integration time is the minimum integration time of the scan pattern for that angle. This advantageously ensures that the beam is scanned at the maximum scan rate and that the return beam is processed at the shortest integration time at each angle, while ensuring the LIDAR system maintains an adequate SNR at each angle.
In another embodiment, during or after step 619, the processor operates the vehicle 310 based at least in part on the data collected by the LIDAR system during step 619. In one embodiment, the processing system 250 of the LIDAR system and/or the processor 314 of the vehicle 310 transmit one or more signals to the steering and/or braking system of the vehicle based on the data collected by the LIDAR system in step 619. In one example embodiment, the processing system 250 transmits one or more signals to the steering or braking system of the vehicle 310 to control a position of the vehicle 310 in response to the LIDAR data. In other embodiments, the processing system 250 transmits one or more signals to the processor 314 of the vehicle 310 based on the LIDAR data collected in step 619 and the processor 314 in turn transmits one or more signals to the steering and braking system of the vehicle 310.
A sequence of binary digits constitutes digital data that is used to represent a number or code for a character. A bus 710 includes many parallel conductors of information so that information is transferred quickly among devices coupled to the bus 710. One or more processors 702 for processing information are coupled with the bus 710. A processor 702 performs a set of operations on information. The set of operations include bringing information in from the bus 710 and placing information on the bus 710. The set of operations also typically include comparing two or more units of information, shifting positions of units of information, and combining two or more units of information, such as by addition or multiplication. A sequence of operations to be executed by the processor 702 constitutes computer instructions.
Computer system 700 also includes a memory 704 coupled to bus 710. The memory 704, such as a random access memory (RAM) or other dynamic storage device, stores information including computer instructions. Dynamic memory allows information stored therein to be changed by the computer system 700. RAM allows a unit of information stored at a location called a memory address to be stored and retrieved independently of information at neighboring addresses. The memory 704 is also used by the processor 702 to store temporary values during execution of computer instructions. The computer system 700 also includes a read only memory (ROM) 706 or other static storage device coupled to the bus 710 for storing static information, including instructions, that is not changed by the computer system 700. Also coupled to bus 710 is a non-volatile (persistent) storage device 708, such as a magnetic disk or optical disk, for storing information, including instructions, that persists even when the computer system 700 is turned off or otherwise loses power.
Information, including instructions, is provided to the bus 710 for use by the processor from an external input device 712, such as a keyboard containing alphanumeric keys operated by a human user, or a sensor. A sensor detects conditions in its vicinity and transforms those detections into signals compatible with the signals used to represent information in computer system 700. Other external devices coupled to bus 710, used primarily for interacting with humans, include a display device 714, such as a cathode ray tube (CRT) or a liquid crystal display (LCD), for presenting images, and a pointing device 716, such as a mouse or a trackball or cursor direction keys, for controlling a position of a small cursor image presented on the display 714 and issuing commands associated with graphical elements presented on the display 714.
In the illustrated embodiment, special purpose hardware, such as an application specific integrated circuit (IC) 720, is coupled to bus 710. The special purpose hardware is configured to perform operations not performed by processor 702 quickly enough for special purposes. Examples of application specific ICs include graphics accelerator cards for generating images for display 714, cryptographic boards for encrypting and decrypting messages sent over a network, speech recognition, and interfaces to special external devices, such as robotic arms and medical scanning equipment that repeatedly perform some complex sequence of operations that are more efficiently implemented in hardware.
Computer system 700 also includes one or more instances of a communications interface 770 coupled to bus 710. Communication interface 770 provides a two-way communication coupling to a variety of external devices that operate with their own processors, such as printers, scanners and external disks. In general the coupling is with a network link 778 that is connected to a local network 780 to which a variety of external devices with their own processors are connected. For example, communication interface 770 may be a parallel port or a serial port or a universal serial bus (USB) port on a personal computer. In some embodiments, communications interface 770 is an integrated services digital network (ISDN) card or a digital subscriber line (DSL) card or a telephone modem that provides an information communication connection to a corresponding type of telephone line. In some embodiments, a communication interface 770 is a cable modem that converts signals on bus 710 into signals for a communication connection over a coaxial cable or into optical signals for a communication connection over a fiber optic cable. As another example, communications interface 770 may be a local area network (LAN) card to provide a data communication connection to a compatible LAN, such as Ethernet. Wireless links may also be implemented. Carrier waves, such as acoustic waves and electromagnetic waves, including radio, optical and infrared waves travel through space without wires or cables. Signals include man-made variations in amplitude, frequency, phase, polarization or other physical properties of carrier waves. For wireless links, the communications interface 770 sends and receives electrical, acoustic or electromagnetic signals, including infrared and optical signals, that carry information streams, such as digital data.
The term computer-readable medium is used herein to refer to any medium that participates in providing information to processor 702, including instructions for execution. Such a medium may take many forms, including, but not limited to, non-volatile media, volatile media and transmission media. Non-volatile media include, for example, optical or magnetic disks, such as storage device 708. Volatile media include, for example, dynamic memory 704. Transmission media include, for example, coaxial cables, copper wire, fiber optic cables, and waves that travel through space without wires or cables, such as acoustic waves and electromagnetic waves, including radio, optical and infrared waves. The term computer-readable storage medium is used herein to refer to any medium that participates in providing information to processor 702, except for transmission media.
Common forms of computer-readable media include, for example, a floppy disk, a flexible disk, a hard disk, a magnetic tape, or any other magnetic medium, a compact disk ROM (CD-ROM), a digital video disk (DVD) or any other optical medium, punch cards, paper tape, or any other physical medium with patterns of holes, a RAM, a programmable ROM (PROM), an erasable PROM (EPROM), a FLASH-EPROM, or any other memory chip or cartridge, a carrier wave, or any other medium from which a computer can read. The term non-transitory computer-readable storage medium is used herein to refer to any medium that participates in providing information to processor 702, except for carrier waves and other signals.
Logic encoded in one or more tangible media includes one or both of processor instructions on a computer-readable storage media and special purpose hardware, such as ASIC 720.
Network link 778 typically provides information communication through one or more networks to other devices that use or process the information. For example, network link 778 may provide a connection through local network 780 to a host computer 782 or to equipment 784 operated by an Internet Service Provider (ISP). ISP equipment 784 in turn provides data communication services through the public, world-wide packet-switching communication network of networks now commonly referred to as the Internet 790. A computer called a server 792 connected to the Internet provides a service in response to information received over the Internet. For example, server 792 provides information representing video data for presentation at display 714.
The invention is related to the use of computer system 700 for implementing the techniques described herein. According to one embodiment of the invention, those techniques are performed by computer system 700 in response to processor 702 executing one or more sequences of one or more instructions contained in memory 704. Such instructions, also called software and program code, may be read into memory 704 from another computer-readable medium such as storage device 708. Execution of the sequences of instructions contained in memory 704 causes processor 702 to perform the method steps described herein. In alternative embodiments, hardware, such as application specific integrated circuit 720, may be used in place of or in combination with software to implement the invention. Thus, embodiments of the invention are not limited to any specific combination of hardware and software.
The signals transmitted over network link 778 and other networks through communications interface 770, carry information to and from computer system 700. Computer system 700 can send and receive information, including program code, through the networks 780, 790 among others, through network link 778 and communications interface 770. In an example using the Internet 790, a server 792 transmits program code for a particular application, requested by a message sent from computer 700, through Internet 790, ISP equipment 784, local network 780 and communications interface 770. The received code may be executed by processor 702 as it is received, or may be stored in storage device 708 or other non-volatile storage for later execution, or both. In this manner, computer system 700 may obtain application program code in the form of a signal on a carrier wave.
Various forms of computer readable media may be involved in carrying one or more sequence of instructions or data or both to processor 702 for execution. For example, instructions and data may initially be carried on a magnetic disk of a remote computer such as host 782. The remote computer loads the instructions and data into its dynamic memory and sends the instructions and data over a telephone line using a modem. A modem local to the computer system 700 receives the instructions and data on a telephone line and uses an infra-red transmitter to convert the instructions and data to a signal on an infra-red a carrier wave serving as the network link 778. An infrared detector serving as communications interface 770 receives the instructions and data carried in the infrared signal and places information representing the instructions and data onto bus 710. Bus 710 carries the information to memory 704 from which processor 702 retrieves and executes the instructions using some of the data sent with the instructions. The instructions and data received in memory 704 may optionally be stored on storage device 708, either before or after execution by the processor 702.
In one embodiment, the chip set 800 includes a communication mechanism such as a bus 801 for passing information among the components of the chip set 800. A processor 803 has connectivity to the bus 801 to execute instructions and process information stored in, for example, a memory 805. The processor 803 may include one or more processing cores with each core configured to perform independently. A multi-core processor enables multiprocessing within a single physical package. Examples of a multi-core processor include two, four, eight, or greater numbers of processing cores. Alternatively or in addition, the processor 803 may include one or more microprocessors configured in tandem via the bus 801 to enable independent execution of instructions, pipelining, and multithreading. The processor 803 may also be accompanied with one or more specialized components to perform certain processing functions and tasks such as one or more digital signal processors (DSP) 807, or one or more application-specific integrated circuits (ASIC) 809. A DSP 807 typically is configured to process real-world signals (e.g., sound) in real time independently of the processor 803. Similarly, an ASIC 809 can be configured to performed specialized functions not easily performed by a general purposed processor. Other specialized components to aid in performing the inventive functions described herein include one or more field programmable gate arrays (FPGA) (not shown), one or more controllers (not shown), or one or more other special-purpose computer chips.
The processor 803 and accompanying components have connectivity to the memory 805 via the bus 801. The memory 805 includes both dynamic memory (e.g., RAM, magnetic disk, writable optical disk, etc.) and static memory (e.g., ROM, CD-ROM, etc.) for storing executable instructions that when executed perform one or more steps of a method described herein. The memory 805 also stores the data associated with or generated by the execution of one or more steps of the methods described herein.
In the foregoing specification, the invention has been described with reference to specific embodiments thereof. It will, however, be evident that various modifications and changes may be made thereto without departing from the broader spirit and scope of the invention. The specification and drawings are, accordingly, to be regarded in an illustrative rather than a restrictive sense. Throughout this specification and the claims, unless the context requires otherwise, the word “comprise” and its variations, such as “comprises” and “comprising,” will be understood to imply the inclusion of a stated item, element or step or group of items, elements or steps but not the exclusion of any other item, element or step or group of items, elements or steps. Furthermore, the indefinite article “a” or “an” is meant to indicate one or more of the item, element or step modified by the article.
The following references were cited herein, the entire contents of each of which are hereby incorporated by reference as if fully set forth herein, except for terminology inconsistent with that used herein.
This application is a continuation of U.S. patent application Ser. No. 16/725,349, filed Dec. 23, 2019, which is a continuation of International Application No. PCT/US2019/043482 filed Jul. 25, 2019, which claims the benefit of and priority to U.S. Patent Application No. 62/711,893, filed Jul. 30, 2018. The entire disclosures of U.S. patent application Ser. No. 16/725,349, International Application No. PCT/US2019/043482, and U.S. Patent Application No. 62/711,893 are hereby incorporated by reference as if fully set forth herein.
Number | Name | Date | Kind |
---|---|---|---|
1086312 | Blaha | Feb 1914 | A |
5371358 | Chang | Dec 1994 | A |
5371581 | Wangler et al. | Dec 1994 | A |
20100053715 | O'Neill et al. | Mar 2010 | A1 |
20100277712 | Zaugg | Nov 2010 | A1 |
20160327635 | Scheim | Nov 2016 | A1 |
20170150029 | Robertson | May 2017 | A1 |
20170299697 | Swanson | Oct 2017 | A1 |
20180113200 | Steinberg | Apr 2018 | A1 |
Number | Date | Country |
---|---|---|
WO-2018107237 | Jun 2018 | WO |
WO-2018144853 | Aug 2018 | WO |
WO-2018160240 | Sep 2018 | WO |
WO-2019014177 | Jan 2019 | WO |
Entry |
---|
International Preliminary Report on Patentability and Written Opinion on PCT Appl. Ser. No. PCT/US2019/043482 dated Feb. 11, 2021 (9 pages). |
International Search Report and Written Opinion on PCT/US2019/043482, dated Oct. 21, 2019, 10 pages. |
Kachelmyer, “Range-Doppler Imaging with a Laser Radar,” The Lincoln Laboratory Journal, 1990, vol. 3, No. 1, pp. 87-118. |
Number | Date | Country | |
---|---|---|---|
20210255290 A1 | Aug 2021 | US |
Number | Date | Country | |
---|---|---|---|
62711893 | Jul 2018 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 16725349 | Dec 2019 | US |
Child | 17099982 | US | |
Parent | PCT/US2019/043482 | Jul 2019 | US |
Child | 16725349 | US |