METHOD AND DEVICE FOR MEASURING THE ALTITUDE OF AN AIRCRAFT IN FLIGHT RELATIVE TO AT LEAST ONE POINT ON THE GROUND

Abstract

A method and a device for measuring altitude of an aircraft relative to a point on the ground, said aircraft carrying a radar system comprising a directional antenna to transmit a radio frequency signal along a aiming axis, including. controlling the transmission of a radiofrequency signal along the axis, calculating received powers as a function of radial distance on a sum channel and an elevation deviation channel, calculating tilt angular deviation values, determining an estimator of the radial distance of the aircraft relative to the point on the ground intercepted by the aiming axis as a function of at least one zero crossing of the angular deviation measurement in a selected area of the angular deviation measurement curve, calculating an aircraft altitude relative to said point on the ground as a function of the estimator of the radial distance and the elevation angle of the aiming axis.

Claims

1. A method for measuring the altitude of an aircraft in flight relative to at least one point on the ground, said aircraft carrying a radar system comprising at least one directional antenna, adapted to transmit at least one radio frequency signal along an aiming axis in a controllable direction, said direction being defined by an elevation angle and a bearing angle, and to receive a reflected radiofrequency signal, the method being implemented by a computing processor and comprising: a) controlling the transmission of a radio frequency signal along a aiming axis having a predetermined elevation angle, b) calculating powers received as a function of a radial distance on two reception channels comprising a first channel, called the sum channel and a second channel, called the elevation deviation channel, c) calculating tilt angle deviation values as a function of the radial distance forming an angular deviation curve, d) determining an estimator of the radial distance of the aircraft relative to a point on the ground intercepted by the aiming axis as a function of at least one zero crossing of the angular deviation in a selected area of the angular deviation curve, and e) calculating an aircraft altitude relative to said ground point intercepted by the aiming axis as a function of the radial distance estimator and the elevation angle of the aiming axis.

2. The method according to claim 1, further comprising computing a barycenter radial distance corresponding to a barycenter associated with a received power curve on the sum channel.

3. The method according to claim 2, wherein said selected angular deviation curve area is defined by a distance interval containing said radial barycenter distance, and wherein the power received on the sum channel is greater than the power received on the elevation deviation channel.

4. A method according to claim 1, wherein the determination of the estimator of the radial distance of the aircraft from the point on the ground intercepted by the aiming axis comprises a determination of a zero-crossing of the angular deviation curve in said selected area, said zero-crossing corresponding to a radial zero-crossing distance, said estimator being equal to said radial zero-crossing distance.

5. The method according to claim 1, further comprising applying filtering to determine said radial distance estimator of the aircraft from the point on the ground intercepted by the aiming axis.

6. The method according to claim 5, wherein said filtering comprises calculating an average value of the radial distances of said selected area of the angular deviation curve, corresponding to angular deviation values less than a predetermined threshold deviation value.

7. The method according to claim 5, wherein said filtering comprises calculating a median value of the radial distances of said selected area of the angular deviation curve, corresponding to angular deviation values less than a predetermined threshold deviation value.

8. The method according to claim 1, comprising a command to transmit radio frequency signals along several aiming axis axes simultaneously, and wherein said steps b) to e) are carried out for each of said aiming axis axes, so as to obtain an altitude profile of a plurality of points on the ground respectively intercepted by each of said aiming axes.

9. A non-transitory memory storing a computer program comprising software instructions which, when implemented by a programmable electronic device, implement an aircraft altitude measurement method according to claim 1.

10. A device for measuring the altitude of an aircraft in flight relative to at least one point on the ground, said aircraft carrying a radar system comprising at least one directional antenna, adapted to transmit at least one radio frequency signal along an aiming axis in a controllable direction, said direction being defined by an elevation angle and a bearing angle, and to receive a reflected radio frequency signal, the device comprising at least one calculation processor configured to implement: a. a module configured to control the transmission of a radio frequency signal along an aiming axis having a predetermined elevation angle, b. a module configured to calculate received powers as a function of a radial distance on two reception channels comprising a first channel called the sum channel and a second channel called the deviation channel, c. a module configured to calculate angular deviation values as a function of radial distances forming an angular deviation curve, d. a module configured to determine an estimator of the radial distance of the aircraft relative to the point on the ground intercepted by the aiming axis as a function of at least one zero crossing of the angular deviation in a selected area of the angular deviation curve, and e. a module configured to calculate an aircraft altitude relative to said point on the ground intercepted by the aiming axis as a function of the estimator of the radial distance and the elevation angle of the aiming axis.

11. The device according to claim 10, further comprising a barycenter calculation module, configured to calculate a barycenter radial distance corresponding to a barycenter associated with a received power curve on the sum channel, said selected area of the angular deviation curve defined by a distance interval comprising said barycenter radial distance, and wherein the received power on the sum channel is greater than the received power on the deviation elevation channel.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0031] Further features and advantages of the invention will be apparent from the description given below, by way of indication and not in any way limiting, with reference to the appended figures, of which:

[0032] FIG. 1 schematically illustrates an application of the invention;

[0033] FIG. 2 is a block diagram of the main functional blocks of a device for measuring the altitude of an aircraft according to an embodiment of the invention;

[0034] FIG. 3 is an example of a graph representing curves of received power on two channels as a function of radial distance;

[0035] FIG. 4 is a flowchart of the main steps of an aircraft altitude measurement process according to an embodiment of the invention, and

[0036] FIG. 5 is an example of a graph representing the angular deviation value as a function of radial distance.

DETAILED DESCRIPTION OF EMBODIMENTS

[0037] FIG. 1 schematically illustrates an embodiment of the invention, and provides a geometric description of the application context of the invention.

[0038] In the example shown in FIG. 1, an aircraft 2, such as an airplane, is shown in flight and moving in a direction (d) represented by a dashed arrow. In the illustrated schematic example, the direction (d) is parallel to the horizontal plane.

[0039] The aircraft 2 is shown in FIG. 1 in a case where it flies over a terrain S, which is irregular and steep.

[0040] The aircraft 2 is equipped with a radar system 4, comprising at least one directional antenna 6, adapted to emit a beam of radio frequency waves or radio frequency signal 8, along an aiming axis A in a controllable direction.

[0041] For example, in one embodiment, the radar system 4 is fixed to the front of the aircraft.

[0042] The direction of the aiming axis A is classically defined by two angles in the orthogonal reference frame (X, Y, Z) shown in FIG. 1, called elevation angle and bearing angle respectively.

[0043] In the example of FIG. 1, the aircraft 2 is moving in a horizontal plane (X, Z), so the aircraft fuselage axis is in the horizontal plane.

[0044] The elevation angle Sant is the angle formed between the aiming axis A and the horizontal plane (X, Z). The bearing angle formed between the aiming axis A and the vertical plane (Y, Z), is not shown in FIG. 1.

[0045] The aiming axis A of the antenna 6 of the onboard radar system 4 intercepts the ground S at a point P, hereafter called the interception point relative to the aiming axis A.

[0046] The method of the invention makes it possible to obtain a precise and robust estimate of the radial distance Dant between the phase center of the transmitting antenna 6 and the point P, on the one hand, and, on the other hand, to deduce from this a measurement of the altitude h at the point P, the distance, vertically from the point P, between the point P and the horizontal plane passing through a reference point of the aircraft, such as the barycenter of the aircraft.

[0047] In a simple way, the following geometrical relation is verified:

h=D.sub.ant×sin (S.sub.ant)   [Math 1]

[0048] where sin( ) is the trigonometric sine function, D.sub.ant is the radial distance value at the estimated intercept point and S.sub.ant is the elevation angle value.

[0049] Precisely estimating the radial distance D.sub.ant in the fixed elevation angle direction S.sub.ant is difficult, since the radio wave beam 8 intercepts an area on the ground around the intercept point P and the radar system thus receives reflected radio waves corresponding to all directions within the beam 8.

[0050] In a known manner, the antenna radiation pattern representing the received power as a function of the estimated radial distance comprises several lobes, and the power is not maximum in the direction of the main lobe. In other words, such an antenna radiation pattern is insufficient for obtaining a robust estimate of the radial distance.

[0051] Furthermore, it should be noted that the power measurements may be noisy, which increases the difficulty of estimating the radial distance D.sub.ant.

[0052] The method of the invention proposes computing an estimator of the radial distance of the aircraft from the ground along the aiming axis.

[0053] FIG. 2 illustrates one embodiment of an onboard system 3 for measuring the altitude of an aircraft in flight according to one embodiment of the invention.

[0054] The onboard system 3 comprises an onboard radar system 4 and a device 10 for measuring the altitude of an aircraft relative to at least one point on the ground, other than a point located at right angles to the vertical of the aircraft.

[0055] The radar system 4 comprises at least one directional antenna 6, whose transmission direction is controllable by a control unit 8.

[0056] For example, the radar system is a mechanically scanned system, and the control unit 8 is a motor adapted to turn the directional antenna in the selected direction.

[0057] In a variant, the radar system 4 is an electronically scanned DFB (Digital Beam Forming) or MIMO radar system.

[0058] The radar system 4 has at least two reception channels, which are respectively a first reception channel, called sum channel, and a second reception channel, called elevation deviation channel.

[0059] These reception channels are known in the field of radars, and particularly in the field of single pulse radars.

[0060] The sum channel, hereafter called the S channel, is a reception channel on which the signals received in each reception unit (such as in each quadrant) of a reception antenna are added.

[0061] For example, a power curve associated with the sum channel represents corresponding power values as a function of measured radial distance values. Each measured distance value corresponds to a sampling time t(i) of the radio frequency signal, converted to radial distance d(i) by the formula:

[00001]$\begin{array}{cc}d\left(i\right)=\frac{c\times t\left(i\right)}{2}& \left[\mathrm{Math}2\right]\end{array}$

[0062] Where c is the speed of light.

[0063] The elevation deviation channel, hereafter called the E channel, is a reception channel on which signals received by some reception units (such as the top two quadrants of a 4-quadrant antenna) of a reception antenna are subtracted from other reception units (such as the bottom two quadrants of a 4-quadrant antenna) of the reception antenna.

[0064] For example, a power curve associated with channel E represents the power values of the radio frequency signal received on channel E for each radial distance d(i) as defined above.

[0065] The powers calculated as a function of radial distance for each of the reception channels are transmitted to an onboard aircraft altitude measurement device 10, according to one embodiment.

[0066] The device 10 is a programmable electronic device, such as an onboard computer.

[0067] The device 10 comprises an electronic computing unit 11, comprising one or more processors, configured to implement computing modules described in more detail below. It also comprises an electronic memory unit 12, configured to store data, such as altitude measurements 14 in association with ground intercept points, forming an altitude profile of a spatial area of the ground, for example, as described in more detail below.

[0068] The calculation modules include in particular a module 16 for calculating the powers as a function of the radial distance for the two reception channels, a module 18 for calculating the angular deviation in elevation, a module 20 for calculating a barycenter associated with the powers received from channel S, a module 22 for calculating a radial distance estimator, and a module 24 for calculating the aircraft altitude relative to a point on the ground corresponding to the estimated radial distance

[0069] As explained in more detail below, the module 22 for calculating a radial distance estimator implements a determination of zero-crossing the angular deviation measurement in an area defined by an interval of radial distances, chosen as a function of the powers on the two reception channels, and the radial distance of barycenter calculated by the barycenter calculation module 20.

[0070] In one embodiment, each of these modules is implemented as software, comprising code instructions executable by the computing unit 11 when the device 10 is powered on.

[0071] According to one embodiment, the programmable electronic device 10 is implemented as an ASIC or FPGA type of programmed board.

[0072] FIG. 3 illustrates an example of received power curves P.sub.S, P.sub.E associated with respective channels S and E as a function of radial distances. The graph in FIG. 3 represents the power in decibels on the ordinate, as a function of the distance in meters on the abscissa.

[0073] In addition, a curve G includes areas G.sub.1 to G.sub.5, in which the power values of the channel S are greater than the power values of channel E, are shown in FIG. 3. It is clear from the Figure that the power curves P.sub.S, P.sub.E obtained in this way are noisy and that it is difficult to obtain a precise assessment of the radial distance corresponding to the aiming axis.

[0074] FIG. 4 is a flowchart of the main steps of a method for measuring the altitude of an aircraft according to one embodiment of the invention.

[0075] The method comprises a first step 30 of controlling the setting of an elevation angle value, S.sub.ant for the aiming axis orientation of a directional antenna of an aircraft-borne radar system (airborne radar system).

[0076] Then, in a power calculation step 32, received power values are calculated and stored, as a function of radial range, for a range interval, with a predetermined sampling step, for each of the first channel (sum channel) and second channel (elevation deviation channel).

[0077] In a variant, in one embodiment, the power is calculated by averaging the powers from received radio frequency signals coming from multiple bearing antenna pointing, within a few degrees.

[0078] The step 32 of calculating the powers of the received radio frequency signals is followed by a step 34 of calculating the angular deviation values as a function of the distance, forming an angular deviation curve. The angular deviation calculation is taken, in a known way, from the amplitudes of the signals received on the sum channel and the elevation deviation channel respectively.

[0079] One example of an angular deviation curve, corresponding to the power curves of FIG. 3, is shown in FIG. 5.

[0080] The elevation angular deviation is expressed in degrees (ordinate of the graph shown in FIG. 5), as a function of the radial distance expressed in meters.

[0081] In a variant, the elevation angular deviation in the field is calculated by averaging the elevation angular deviations from multiple bearing antenna pointing, within a few degrees.

[0082] The method further comprises a step 36 of calculating a barycenter associated with the sum channel, making it possible to obtain a barycenter radial distance value Do.

[0083] The radial barycenter distance is calculated from the powers calculated for the sum channel, for example, by the formula:

[00002]$\begin{array}{cc}{D}_0=\frac{{.Math.}_{i}d\left(i\right)\times Ps\left(i\right)}{{.Math.}_{i}Ps\left(i\right)}& \left[\mathrm{Math}3\right]\end{array}$

[0084] where i is the distance box index, d(i) is the radial distance for the distance box of index i (according to the formula [Math 2] and Ps(i) is the power of the Sum path for the distance box of index i.

[0085] Advantageously, the barycenter associated with the sum channel belongs to the main lobe of the antenna pattern of the onboard radar system antenna.

[0086] The method further comprises a determination 38 of the radial distance(s) for which the angular deviation value is equal to zero in a selected area of the angular deviation curve.

[0087] The selected area is one of the areas G.sub.i in which the power received on the sum channel (S channel) is greater than the power received on the elevation deviation channel (E channel), and more particularly the area defined by an interval of distances [D.sub.min, D.sub.max], including the barycenter radial distance D.sub.0 calculated in the barycenter calculation step 36.

[0088] The distance D.sub.max is the greatest distance of the area in which the power received on the sum channel is greater than the power received on the elevation deviation channel, and including the radial barycenter distance D.sub.0.

[0089] The distance D.sub.max is the estimated maximum radial distance, in one embodiment.

[0090] In FIG. 5, as an example, the selected area Z has been illustrated, as well as the distances D.sub.min, D.sub.max and the radial distance of barycenter D.sub.0 on the x-axis. The chosen area Z corresponds to the area G2 in FIG. 3.

[0091] The radial distance value corresponding to the zero crossing of the angular deviation measurement in the chosen area Z is an estimator of the radial distance D.sub.ant of the aircraft relative to the ground along the aiming axis.

[0092] In one embodiment, the method further comprises, in the event that multiple zero crossings of the angular deviation are detected in the selected area, filtering (step 39) the corresponding radial distance values to determine the radial distance estimator.

[0093] For example, the filtering implemented in the filtering step 39 consists in calculating the average value of the radial distances values corresponding angular deviation values lower in absolute value than a predetermined threshold deviation value, that is, comprised between −ε and ε, with ε being worth a fraction of the antenna aperture in elevation, preferably a few tenths of the antenna aperture in elevation, equal to 0.1°, for example.

[0094] According to one embodiment, the filtering consists in selecting the median value of the radial distances corresponding to angular deviation values lower in absolute value than a predetermined threshold deviation value, comprised between −ε and ε.

[0095] According to another embodiment, in step 38 of determining a zero crossing of the angular deviation in the selected angular deviation area, a polynomial regression (of degree 1 or higher) is applied.

[0096] The method further includes a step 40 of calculating an aircraft altitude as a function of the radial distance estimator D.sub.ant and the elevation angle S.sub.ant of the aiming axis, by applying the formula explained above [Math 1], with the elevation angle S.sub.ant fixed.

[0097] The aircraft altitude value associated with the interception point of the line of the sight P on the ground is stored at step 42. The interception point P is in a predetermined spatial reference frame defined by coordinates, for example, knowing the position of the aircraft in the spatial reference frame.

[0098] In a variant or additionally, knowing the aircraft altitude, it is possible to calculate and store the altitude of the interception point P as a function of the previously calculated altitude.

[0099] Optionally, steps 30 to 42 are repeated for several antenna aiming axis directions, and the altitude and/or altitude values relative to the interception points on the ground are stored, which makes it possible to obtain an altitude profile 44 area on the ground illuminated by the radar system, located in front of the aircraft, for example.

[0100] When the radar system is beam forming by computing, power measurements for several sighting directions are acquired at the same time, which makes it possible to obtain an instantaneous altitude profile 44 of the area on the ground illuminated by the transmitted wave beam, located in front of the aircraft, for example.

[0101] Advantageously, in the landing phase, thanks to the method of the invention, an instantaneous profile of an area on the ground, such as the area located in front of the aircraft, which is a landing strip, for example, is available on board the aircraft. This profile is computed autonomously and the estimated altitudes are precise, with the altitude estimation error being less than one meter.