Anticollision radar, especially for an aircraft when taxiing, and anticollision system
10585186 ยท 2020-03-10
Assignee
Inventors
Cpc classification
G01S13/90
PHYSICS
International classification
G01S13/34
PHYSICS
Abstract
A radar attached laterally to airplane fuselage to detect obstacles on a collision course with a portion of the airplane facing the radar. The radar includes an emission antennal channel and reception antennal channels in the same plane. The radar Establishing in a radar coordinate system a first distance/Doppler map allowing echoes to be separated into distance and Doppler resolution cells in reception antennal channels; Establishing a second distance/Doppler map of smaller size by selecting a subset of distance/Doppler resolution cells corresponding to possible positions of targets liable to collide with the airplane; Establishing new distance/Doppler maps by forming beams computationally from the subsets of distance resolution cells; and Temporal integration, in each beam, of successive distance/Doppler maps.
Claims
1. An anticollision radar for a taxiing aircraft, wherein the anticollision radar is attached laterally to a body of an aircraft to detect obstacles on a collision course with a portion of the aircraft, the anticollision radar being oriented along a sighting axis perpendicular to an axis of a fuselage of the aircraft; the anticollision radar including at least one emission antenna channel emitting a wave of a FMCW (frequency modulation continuous-wave) type and a plurality of reception antenna channels placed in a same plane, a form of said wave and a field of angular coverage of said anticollision radar depending on a velocity of the aircraft; said anticollision radar processing at least the following steps: establishing in a coordinate system of the anticollision radar a first distance/Doppler map allowing echoes to be separated into distance and Doppler resolution cells in said plurality of reception antenna channels; establishing a second distance/Doppler map of smaller size than the first distance/Doppler map by selecting from the first distance/Doppler map a subset of distance and Doppler resolution cells corresponding to positions of targets liable to collide with said portion of the aircraft, said selection depending on the velocity of the aircraft; establishing new distance/Doppler maps by forming beams computationally from subsets of the distance/Doppler resolution cells that are retained for each reception antenna channel, so as to separate a beam of angular observation in a plane of the antenna channels into at least two beams utilizing Doppler filtering in unfocused synthetic antenna processing; temporal integration, in each beam formed, of a plurality of successive distance/Doppler maps obtained in the establishing new distance/Doppler maps step in order to obtain an integrated distance/Doppler map for each beam, an integration duration depending on a waveform used; and at least one step of detecting in each of the beams formed on a basis of said integrated distance/Doppler map.
2. The anticollision radar as claimed in claim 1, wherein in the detecting step, an azimuth angle of a target is estimated, for each detected pixel in the integrated distance/Doppler map, by comparing a received level in the distance/Doppler maps of two adjacent formed beams using amplitude-comparison monopulse.
3. The anticollision radar as claimed in claim 1, wherein the detecting step is followed by a step of tracking detected targets in the coordinate system of said anticollision radar.
4. The anticollision radar as claimed in claim 1, wherein the anticollision radar changes coordinate systems, passing from the coordinate system associated with the anticollision radar to a coordinate system associated with said portion of the aircraft, in order to estimate coordinates of targets in the coordinate system associated with said portion, which positions said targets relative to said portion of the aircraft.
5. The anticollision radar as claimed in claim 1, wherein a plane of the reception antenna channels is a horizontal plane.
6. The anticollision radar as claimed in claim 1, wherein a beam forming by the antenna channels comprises an angular field of +/45 in the plane of the antenna channels, relative to a sighting axis, that is perpendicular to an axis of the body of the aircraft.
7. The anticollision radar as claimed in claim 1, wherein antennae of the reception antenna channels in which the beam forming is carried out computationally are spaced apart from one another by more than one half wavelength, and wherein an ambiguity to localize a target being removed by a Doppler frequency of a distance/Doppler cell of said target.
8. The anticollision radar as claimed in claim 1, wherein the first distance/Doppler map is obtained after distance compressing and Doppler compressing said wave of the FMCW type.
9. The anticollision radar as claimed in claim 1, comprising four reception antenna channels.
10. The anticollision radar as claimed in claim 1, wherein said wave of the anticollision radar operates in a millimeter frequency domain.
11. The anticollision radar as claimed in claim 10, wherein said wave of the anticollision radar operates in a frequency band of 76 to 77 GHz.
12. The anticollision radar as claimed in claim 1, wherein said portion is a wingtip of the aircraft.
13. An anticollision system for a taxiing airplane, comprising at least two radars as claimed in claim 1, each installed laterally per side of the body of the airplane forming its fuselage.
14. The anticollision system as claimed in claim 13, wherein each radar is installed in front of a wing.
15. The anticollision system as claimed in claim 14, wherein each radar is attached in a location of a porthole in a vicinity of a cockpit.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) Other features and advantages of the invention will become apparent from the following description, which is given with regard to the appended drawings, which show:
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
DETAILED DESCRIPTION
(15)
(16) The system is composed of two radar sensors 1, 2 according to the invention, which sensors are each installed laterally on one side of the fuselage of the airplane 10. The radars 1, 2 are oriented along a sighting axis 8, 9 substantially perpendicular to the axis of the fuselage so that their beams 3, 4 cover the zones located in front of the wing. Each radar is for example placed in the vicinity of the cockpit so as to detect obstacles on a collision course with the wingtips 5, 6. It may especially be installed in the location of a porthole. The processing according to the invention and which will be described below advantageously makes it possible to avoid having to place the radar on the ends of the wings. The same processing means that it is only necessary to install one radar per side of the fuselage. Each radar is autonomous and for example delivers warning information directly to the pilot, without needing to fuse the information from a plurality of radars.
(17) The field of coverage of the radars is optimized to ensure obstacles are detected as far as possible from and continuously to as close as possible to the wingtip. Thus, a radar 1, 2 possesses a field that is wide in azimuth, for example of about 90. More particularly, on account, on the one hand, of the installation geometry and, on the other hand, of the maximum angular aperture achievable for a static radar antenna, the angular coverage of the radar in the horizontal plane extends for example by +/45 about the sighting axis 8, 9. Under these conditions, the wingtip is protected from a zero distance to a maximum distance Dmax, typically of about 90 meters, corresponding to the radar range Rmax of about 70 meters. On account of the braking distances required to stop the carrier, the anticollision system must for example be effective up to a velocity of about 10 meters/second.
(18)
(19)
(20) The radar comprises at least one emission antennal channel 31 and a plurality of reception antennal channels 32 preferably distributed in the same horizontal plane. In one possible embodiment, the radar includes one emission channel and four antennal channels 320, 321, 322, 323 as illustrated in
(21) The radar such as illustrated in
(22) The radar uses a FMCW (frequency modulation continuous-wave) waveform the characteristics of which are adapted depending on the needs of the airplane application, relative to the required resolution. So as to optimize both the ambiguity function of the radar and the signal-to-noise ratio for various collision hypotheses, the processing adapts the waveform to the current velocity of the carrier.
(23)
(24) According to the invention, the waveform of the radar is set to one of these forms FO.sub.0, FO.sub.1, FO.sub.2, FO.sub.3 depending on the velocity of the carrier. More particularly, the characteristics of the waveform depend on the velocity of the carrier, these characteristics being essentially the duration of the ramp T and the modulation band F. One example of four waveforms FO.sub.0, FO.sub.1, FO.sub.2, FO.sub.3 is given in the appendix.
(25) At constant processing load, it is advantageous at high velocity, i.e. when detection at large distances is desired, to privilege the Doppler resolution of the radar to the detriment of distance resolution. In contrast, at low velocity, it is preferable to use a high distance resolution, separation by velocity not being very discriminate. This aspect will be described below with regard to processing operations carried out on the signal received by the radar.
(26)
(27) In the first phase 51, the radar performs for each antennal channel a step 511 of distance and Doppler compressing, a step 512 of selecting a distance-Doppler zone, a computational beam forming step 37, then a step 514 of integrating the distance-Doppler map.
(28)
(29) The compression is carried out after the signal received has been demodulated with the emission signal and encoded. The compression itself includes two steps. The first step consists in carrying out the distance compression by fast Fourier transform (FFT) 62 to a temporal depth NFF.sub.CD corresponding to the modulation duration of each of the frequency ramps 61. In each distance box 63 a distance value biased by the Doppler effect is obtained. The distance FFT of this first step are therefore followed in a second step by Doppler FFT to correct this bias.
(30) The second step therefore allows the signal received in each of the distance boxes 63 thus formed to be analyzed by carrying out a second fast Fourier transform 64 to a temporal depth corresponding to a number NFFT.sub.Dop of successive ramps. At the end of this second step, for each reception antenna, a distance-Doppler map 65 that is NFF.sub.CDNFFT.sub.Dop in size is obtained. The distance angular processing and the Doppler processing may be inverted, i.e. it is possible to start by carrying out the Doppler FFT and then the distance FFT.
(31) The separation of the targets on the Doppler axis is used to separate the signals according to their direction of arrival, in conformity with the principle of lobe sharpening by Doppler filtering in unfocused synthetic antenna processing. Specifically, each of the columns of this distance-Doppler matrix 65 may be likened to echoes attached to a subbeam in a given azimuthal direction.
(32) As
Vr=Va sin(.sub.az)(1)
(33) where Va is the longitudinal velocity of the taxiing airplane.
(34) Therefore, its Doppler frequency is written:
(35)
(36) where represents wavelength.
(37) Two fixed points separated angularly by an angular deviation az generate two Doppler frequencies separated by:
(38)
(39) Thus, if the observation duration used for the calculation of the Doppler Fourier transform is Te=1/Fd.
(40) Each column of the Doppler distance matrix 65 corresponds to the output of a subbeam typically of width
(41)
centered on the azimuthal angle
(42)
where Fd is the central frequency of the filter in question.
(43) If a fixed point located at an azimuth .sub.az and at a distance R from the radar is considered, the Fourier transform on the Doppler axis makes it possible to achieve a resolution in transverse distance:
(44)
For example, for: =3.9 mm R=40 m Va=10 m/sec Te=50 ms az=0
a transverse resolution of r.sub.t=15.6 cm is obtained.
(45) In the choice of the waveform described above, the invention advantageously uses the fact that the angular resolution, which translates into a resolution in transverse distance as seen from the radar, is inversely proportional to the velocity of the carrier. In the case where the observed target is moving at a radial velocity Vro relative to the radar, the angular separation by Doppler filtering is unchanged, but because of the angle/velocity ambiguity, the direct estimation of the angle from the measurement of the Doppler frequency is biased because:
(46)
where .sub.azestim is the estimation of the azimuthal angle obtained directly from the measurement of the Doppler frequency Fd.
(47) The Doppler compression therefore has the main objective of providing the radar with transverse resolution, radial resolution being obtained by the width of the emitted band.
(48) The angular location must for its part be refined by a complementary means, assuming that it is desired to detect both static and moving targets. This is done by selecting a distance-Doppler zone.
(49)
(50) As indicated above, in order to take advantage of the angular discriminating power of the Doppler axis to concentrate on certain angular zones and to limit computational power, a distance-Doppler zone of interest 81 is selected in this step 512. To select this zone, the processing takes into account the current velocity of the carrier with respect to a distance-Doppler zone of interest, this corresponding, as seen from the radar, to a limited distance swath and a limited angular field. Exemplary selection criteria are described below. Only obstacles located, relative to the fuselage, at a lateral distance in the vicinity of the distance of the wingtip are liable to cause a collision, and must be taken into consideration. For example, the extent of the zone of danger may be considered to extend plus or minus 10 meters on either side of the wingtip. This criteria makes it possible to define a minimum radar swath limiting the zone of processing on the radial distance axis, as seen from the radar. Furthermore, on account of reaction time and braking distance, the minimum longitudinal detection distance relative to the wingtip for a first alarm level depends on the velocity of the carrier. Moreover, the radar processing requires a minimum number of successive detections to reliably detect the presence of an obstacle with a very low false alarm risk. These two constraints allow the upper distance limit of the detection zone to be determined, as seen from the wingtip. This limit is for example about 75 m for a velocity Va of 10 m/sec, and about 20 m for a velocity Va of 2.5 m/s. It makes it possible to determine the first limit of the minimum detection angular field in the horizontal plane, as seen from the radar. Lastly, it is necessary to continue to follow the position of the obstacle until a minimum distance corresponding to the braking distance is reached, in order to verify that said object remains on a collision course. This constraint allows the lower distance limit of the detection zone, as seen from the wingtip, to be determined. This limit is typically 55 m at 10 m/sec and 10 m at 2.5 m/sec. It makes it possible to determine the second limit of the minimum detection angular field in the horizontal plane, as seen from the radar.
(51)
(52) To these detection fields, defined by distance and azimuthal angle limits, correspond directly, for static obstacles, a limited distance and Doppler field in the distance-Doppler map. However, if detection of mobile obstacles is required, it is necessary to widen the Doppler field to take into account the contribution made by the radial velocity of these obstacles, in accordance with relationship (6). Moreover, if it is envisioned to detect mobile objects entering tardily into the field of view of the radar, it is also necessary to increase angularly the limits of the detection field. Thus, according to the invention, the limits of the Doppler field are established depending on the velocity of the carrier 10, for example using the table below.
(53) TABLE-US-00001 Limits of the angular field Velocity of the relative to the Velocity field to Doppler field to carrier sighting axis 8 consider consider 7.5 to 10 m/s 0 to 45 10 m/s to 0 m/s 4 KHz to 0 KHz 5 to 7.5 m/s 35 to 45 10 m/s to 4.3 m/s 3 KHz to 2.1 KHz 2.5 to 5 m/s 45 to 10 5.4 m/s to 3.5 m/s 2.1 KHz to 1.7 KHz 0 to 2.5 m/s 45 to 0 5.4 m/s to 1.8 m/s 2 KHz to 0.9 KHz
(54) The limits of the distance field extend from 30 meters to 70 meters, whatever the velocity of the carrier 10.
(55) The following step 37 is the step of forming azimuthal beams. In this step, the radar therefore establishes new distance/Doppler maps by forming beams computationally from the subsets of those distance resolution cells that are retained for each antennal channel, so as to separate the beam of angular observation in the horizontal plane into a plurality of, at least two, subbeams. The processing of the preceding steps 511, 512 was carried out independently on each of the four reception channels, up to this step of forming beams computationally. In this step, the processing means of the radar focus, by coherent integration on reception, the radar beam in a plurality of different directions in the distance-Doppler maps constructed beforehand, one map corresponding to each reception channel. This operation may be carried out by the operator defined by the following relationship (8) on each of the pixels (i, j) of the input distance-Doppler maps:
(56)
where: N.sub.ant: the number of reception antennae d: the distance between each pair of elementary antennae : the wavelength .sub.az: the direction in which the beam is formed i, j: the indices of the pixels along the Doppler axis (i) and distance axis (j) s(i, j): the signal in amplitude and phase of the distance/Doppler box i, j.
(57) For example, in the exemplary embodiment with 4 reception channels, the beams may typically be formed in the directions: .sub.az=+45, .sub.az=+15, .sub.az=15 and .sub.az=+45, from four distance Doppler maps.
(58) This processing thus delivers four distance/Doppler maps 81 that are (N.sub.doppSel, N.sub.distSel) in size in each of the aforementioned directions.
(59)
(60) By virtue of the processing carried out by a radar according to the invention, it is also possible to improve the resolution of the reception beams including in the case of large deflections, by spacing the reception antennae by a distance larger than /2, this not usually being possible because of the presence of grating lobes that tend to cause large deflections to appear in detections located in the sector opposite their actual position relative to the axis of the antenna. This situation is shown in
(61)
(62) One objective of this step 514 is especially to increase the signal-to-noise ratio by non-coherent temporal integration of the distance-Doppler images issued from the preceding step, i.e. in each of the sighting directions. More particularly,
(63) The integration is carried out on each of the distance and Doppler axes of each map 81, while taking account, during the integration time, of possible distance and Doppler migrations, and while limiting the integration to the temporal movements of pixels possibly corresponding to collision courses: On the Doppler axis: Each (i, j) output pixel is obtained by non-coherent integration of N.sub.intDopp Doppler pixels located in a given distance box of the input map, where N.sub.intDopp is chosen depending on the maximum Doppler migration during the integration time. In other words, an integration time is set and the number of filters the Doppler target possibly moves is observed. It is moreover also possible at this level to limit the integration solely to pixels that have a predefined signal-to-noise ratio that is judged to be enough to correspond to a presumption of presence of a useful signal. On the time axis: All the pixels of the distance/Doppler maps issued from the preceding step are summed non-coherently to a temporal depth defined by a parameter N.sub.NC. This summation may optionally be a rolling summation.
(64) An output pixel x.sub.out (i, j) is thus obtained by non-coherent integration according to the following relationship (9):
x.sub.out(i, j)=.sub.n=0.sup.N.sup.
where: i: index on the Doppler axis j: index on the distance axis n: index on the time axis 110
(65) The duration of the non-coherent integration may be 250 ms to 300 ms depending on the waveform used.
(66) The numerical values applied to the two parameters N.sub.NC and N.sub.intDop are for example the following N.sub.NC=5 (i.e. 250 ms of non-coherent integration); N.sub.intDopp depends on the current waveform and varies in accordance with the following table:
(67) TABLE-US-00002 Number Number of integration Angular resolution of Doppler boxes (N.sub.intDOP) on the (function of V.sub.a) boxes in the Waveform Doppler axis [.sub.cross 70 m] output map FO.sub.0 10 [2, 3] 26 [3 m, 4 m] FO.sub.1 8 [1.8, 2.7] 50 [2.2 m, 3.2 m] FO.sub.2 2 [0.8, 1.8] 102 [1 m, 2 m] FO.sub.3 1 [0.9, 4] 88 [1 m, 5.4 m]
(68) In the very high-resolution case corresponding to the last row of the first column, N.sub.intDopp=1 is obtained because no Doppler migration is possible and the target cannot change Doppler filter. In the first row of the third column, a number of Doppler boxes equal to 26 is obtained, this number being related to the angular resolution (here between 2 and 3) and corresponding to the resolution. One filter corresponds to one given direction. In this case, there are 26 filters of interest for an angular resolution comprised between 2 and 3.
(69) At the end of this processing, a distance-Doppler map 111 is obtained for each sighting direction, corresponding to each beam formed.
(70) The following steps 52, 53, 54 consist in detecting then in characterizing objects present in the processed field, in terms of distance, velocity and angular position.
(71) The first 52 of these steps consists in detecting, in each of the beams formed, the presence of a possible target. The detection is carried out in a conventional way, by comparing each pixel relative to noise, this comparison optionally being completed by validation of a criterion of contrast with neighboring pixels or detection of outlines in the distance Doppler maps 111 issued from the preceding processing steps. This operation is carried out for each beam formed.
(72) The velocity and location of useful targets are estimated in a conventional way in the following steps:
(73) For each pixel of indices i, j (see
(74) The composite distance Doppler frequency F.sub.dist is estimated: The index i of the pixel giving rise to a detection is extracted in order to identify the central frequency of the distance cell 121, and the measurement is refined by an estimation of contrast with adjacent pixels on the distance axis. The Doppler frequency F.sub.dop is estimated: The index j of the pixel giving rise to a detection is extracted in order to identify the central frequency of the distance cell 121, and the measurement is refined by an estimation of contrast with adjacent pixels on the Doppler axis. The radial distance Drad is estimated using the following relationship (10):
(75)
where: T.sub.m is the duration of the ramp of the waveform of index m; and F.sub.m is the modulation band of the ramp of index m.
(76) In the example of the four waveforms in
(77) Thus the radial distance echoes are rectified and bias associated with the Doppler shift corrected, hence the composite frequency term dependent on distance and Doppler shift.
(78) In the example in
(79) Advantageously, the azimuthal angle is estimated for each detected pixel by comparing the received level in the distance/Doppler maps of two adjacent formed beams, for example using a known amplitude monopulse type operator.
(80) The detecting and locating step 52 is for example followed by a tracking step 53. The tracking is achieved by known processing in the coordinate system 131 of the radar. The targets are tracked so as to consolidate the detections, to refine their obtained velocity and position measurements, and to estimate their paths.
(81) In a following step 54 a change of coordinate system is carried out, passing from the coordinate system 131 associated with the radar to the coordinate system 132 associated with the wingtip of the airplane as illustrated in
(82) At this stage, only detections on a collision course are for example retained. The paths positioned in this new coordinate system may then be exploited by the anticollision system in order to be displayed and/or allow a pre-alarm or urgent alarm decision to be taken.
(83) The invention was described for anticollision relative to the wingtips of a taxiing airplane. It also applies to any other portion of airplanes provided that they are located facing the radar.
(84) Lastly, the invention may also apply to a taxiing machine other than an aircraft. The description relating to the fuselage of the airplane then applies to the body of the machine. In particular, the invention may advantageously be applied to taxiing machines with large wingspans, in order to prevent them colliding with various potential obstacles.
APPENDIX
(85) This appendix presents a guide example of choices for the waveforms in
(86) Waveform FO.sub.0: applied for V.sub.a comprised between 7.5 and 10 m/s
(87) F=50 MHz in order to obtain a resolution of 3 m;
(88) F.sub.dopp=8 kHz considering the entirety of the beam (100);
(89) F.sub.r=1.3 F.sub.dopp=10.2 kHz (T.sub.r=98 s) to respect the constraint F.sub.rec>F.sub.dopp;
(90) T=T.sub.r2 s=96 s; loop stabilization rest time
(91) F.sub.ech=625 kHz (T.sub.ech=1.6 s), the maximum sampling frequency of the radar;
(92) In this case, the distance of ambiguity is located at 90 m.
(93) Velocity field at 3 dB (100): 7.7 m/s, 7.7 m/s
(94) Angular resolution: [0.22, 029] where
(95)
and T.sub.int=50 ms
Namely a transverse resolution at 70 m of [0.27 m, 0.35 m]
Waveform FO.sub.1: applied for V.sub.a comprised between 5 and 7.5 m/s
F=65 MHz in order to obtain a resolution of 2.3 m;
F.sub.dopp=6 kHz considering the entirety of the beam (100);
F.sub.r=1.3 F.sub.dopp=7.6 kHz (T.sub.r=128 s) to respect the constraint F.sub.rec>F.sub.dopp;
T=T.sub.r2 s=126 s; loop stabilization rest time;
F.sub.ech=625 kHz (T.sub.ech=1.6 s), the maximum sampling frequency of the radar;
In this case, the distance of ambiguity is located at 90 m.
Velocity field at 3 dB (100): 5.7 m/s, 5.7 m/s
Angular resolution: [0.29, 0.44]
Namely a transverse resolution at 70 m of [0.35 m, 0.53 m]
Waveform FO.sub.2: applied for V.sub.a comprised between 2.5 and 5 m/s
F=100 MHz in order to obtain a resolution of 1.5 m;
F.sub.dopp=4 kHz considering the entirety of the beam (100);
F.sub.r=1.3 F.sub.dopp=5.5 kHz (T.sub.r=180 s) to respect the constraint F.sub.rec>F.sub.dopp;
T=T.sub.r2 s=178 s; loop stabilization rest time;
F.sub.ech=625 kHz (T.sub.ech=1.6 s), the maximum sampling frequency of the radar;
In this case, the distance of ambiguity is located at 83 m.
Velocity field at 3 dB (100): 4.2, +4.2 m/s
Angular resolution: [0.4, 0.9]
Namely a transverse resolution at 70 m of [0.5 m, 1.1 m]
Waveform FO.sub.3: applied for V.sub.a lower than 2.5 m/s
F=200 MHz in order to obtain a resolution of 0.75 m;
F.sub.dopp=2 kHz considering the entirety of the beam (100);
F.sub.r=1.3 F.sub.dopp=2.56 kHz (T.sub.r=195 s) to respect the constraint F.sub.rec>F.sub.dopp;
T=T.sub.r2 s=388 s; loop stabilization rest time;
F.sub.ech=625 kHz (T.sub.ech=1.6 s), the maximum sampling frequency of the radar;
In this case, the distance of ambiguity is located at 90 m.
Velocity field at 3 dB (100): 2, +2 m/s
Angular resolution: [0.90, 5]
Namely a transverse resolution at 70 m of [1.1 m, 6.1 m]