Method for interferometric radar measurement

[0001] This application is a continuation in part application of U.S. patent application Ser. No. 09/889,759 filed on Jul. 20, 2001, now abandoned incorporated herein by reference wherein priority is claimed under 35 U.S.C. 120, this application also claims priority under 35 U.S.C. §371 from PCT/DE99/04066 filed on Dec. 22, 1999, and also claims priority under 35 U.S.C. 119 from German Patent Application DE 199020078 filed on Jan. 21, 1999.

[0002] The invention relates to a method for interferometric radar measurements.

[0003] Due to their construction, radar devices are precise range-finding systems, which means that without special measures, a radar device is capable of determining only the distance of a target from the antenna, but not its direction. It is possible to determine only whether a target is present within the lobe of the antenna.

[0004] This problem is eliminated to a large extent in conjunction with the ROSAR or Heli-Radar system known until now by using, for example 16 vertically staggered antennas with an antenna opening angle of, for example 2.5°. It is possible with this type of a system to determine the location of an elevated obstruction etc. within an accuracy of about 2.5° in terms of elevation. However, in this case, targets located at the same distance are also displayed in the same antenna in the same image spot.

[0005] The azimuthal resolution of the known Heli-Radar system amounts to about 0.2° because of a special signal processing. These features are taught in the disclosure in DE 39 22 086 C1. However, the direction of an obstruction and thus the location in space at which this obstruction is located can be determined only with the help of a triangulation, whereby in the simplest case, two locally separated radar installations can be used for this purpose.

[0006] However, it is also possible to use the properties of a coherent radar system and to perform a triangulation with the help of the phase of the emitted signal. Thus, a coherent radar system is used which coherently transmits a signal via a transmitting antenna and receives the echoes scattered back via two locally separated receiving antennas. A coherent evaluation permits a calculation of the phase difference between receiving signals. The direction from which the scattered echoes are received is determined based on the phase difference. Now, once the distance and direction of an “obstruction” have been computed, its location in space can be determined as well. This type of three-dimensional determination of a location with the help of a coherent radar system which comprises one transmitting antenna and two receiving antennas is generally referred to as “radar interferometry” and known for a long time. It is used already for the generation of topographic charts with the help of SAR-systems installed on aircraft, for example by the DOSAR system of the firm Dornier GmbH.

[0007] With these designs, reference is also made to the following published documents pertaining to the further state of the art:

[0008] C. T. Allan, Review Article, Interferometric Synthetic Aperture Radar, in IEEE Geoscience and Remote Sensing Society News Letter, September 1995, p. 6 ff; S. Buckreuss, J. Moreira, H. Rinkel and G. Waller, Advanced SAR Interferometry Study, DLR Bulletin 94, Jun. 10, 1994, Institut für Hochfrequenztechnik, Oberpfaffenhofen.

[0009] The entire prior art known to this date and the state of the art cited above, including the ROSAR principle on which the present invention is based, projects terrain elevations or other elevated obstructions in one plane, so that it is not possible to recognize the elevation of the given obstruction if the reproduced topography of the terrain present is unknown. However, a three-dimensional image is required for controlling flights.

[0010] The present invention is based on the problem of proposing measures on the basis of the ROSAR principle that permit a quasi-three-dimensional image representation of terrain and other obstructions.

[0011] The problem is solved with an arrangement for interferometric radar measurement having a transmitter and two assigned coherent receiving antennas with receiving channels. The transmitter and receiving antennas are arranged on the turnstile of the ROSAR system of a helicopter radar. An additional transmitting/receiving antenna is provided wherein this additional transmitting/receiving antenna is sharply focused in the elevation direction. The transmitter and receiving antennas are arranged at the end of the turnstile. The receiving antennas are also arranged at the end of the turnstile.

[0012] The arrangement includes a process whereby two coherent receiving antennas with receiving channels are assigned to a first transmitter, and the path length difference of the two distances can be calculated to measured receiving point P from the wave length of the transmitted radar signal and of the measured phase difference of the reception echo of both coherent receiving channels. A helicopter operating according to the ROSAR principle is used for the interferometric radar measurement, whereby two coherent receiving antennas are assigned to a transmitter of the ROSAR system arranged on a rotating turnstile on the radar. Additionally, receiving signals of the sharply focused or second transmitting/receiving antenna can be evaluated for determination of the phase difference.

[0013] The sight angle is used for calculating the coordinates of the respective receiving point for representing the image points on the integrated graphic display screen in the ROSAR system. The antennas and the center of the image on the graphic display screen are in a fixed relationship to each other.

[0014] The first transmitter along with the two receiving antennas are used to determine the location of point P. However, this determination of the location of point P may contain inaccuracies or ambiguities because these two receiving antennas have a value of a measured phase difference that is ambiguous wherein this ambiguity can only be determined down to a value ranging between 0 and 2π. Therefore, this ambiguity must be determined by additional measurements. Thus to provide these additional measurements there is the additional or sharply focused transmitting and receiving antenna. This transmitting and receiving antenna is located in the end of the antenna cross and wherein this antenna has a narrow beam. This antenna covers the lower range of the sight angle. The distance to the impact point on the ground can be clearly determined by the receive echo of this additional transmitting/receiving antenna. Thus, this INROSAR-system accepts the distance as a basic value and calculates the further ambiguities based on the rising distance from the continuous phase transitions.

[0015] Thus, this system is a system with a first antenna transmitter with two receiving antennas and then one additional transmitting and receiving antenna that is used to dispel any ambiguities of the signal received by the first antenna transmitter and receiver.

[0023] According to the general idea of the invention, the goal is to obtain in conjunction with a helicopter operating according to the ROSAR-system a quasi-three-dimensional radar image representation for flight guidance by associating with a transmitter located on the rotating rotary cross two coherent receiving antennas with receiving channels.

[0024] The ROSAR-system used heretofore is comprised of, for example 16 transmitters and receivers with their channels for obtaining a three-dimensional image. However, these transmitters and receivers have a directional inaccuracy of about 2.5°. Now, if said ROSAR-system, as mentioned above, is expanded by a highly precise coherent receiving channel, only one transmitter and two coherent receivers instead of the, for example sixteen transmitters and receivers employed until now will be required for obtaining the highly precise three-dimensional radar image. The directional inaccuracy found until now is enhanced by the interferometric principle by about a factor of 100.

[0025] This is explained in the following description of an exemplified embodiment of the invention, which is sketched in FIG. 1.

[0026] A helicopter operating according to the ROSAR principle flies over the surface of the earth at an altitude H. One transmitting antenna and two receiving antennas with associated coherent transmitting and receiving electronics are mounted on the end of the rotating antenna cross. The received echoes are amplified, digitized and processed further.

[0027] The distance between this arrangement as described above, which is referred to in the following as the INROSAR-system, and the impact point P, which is located at a relative altitude h, is referred to as R. The distance from the antenna A₁ of the INROSAR-system to the impact point P amounts to R+ΔR and is therefore by a small amount ΔR greater than the distance R to the antenna A2. The difference ΔR between the two distances can be calculated based on the known wavelength λ of the emitted radar signal and the measured phase difference φ of the receiving echo of the two coherent receiving channels. Impact point Q has a relative altitude lower than impact point P and is on the surface of the Earth E.

[0028] Now, this phase difference φ of the receiving echo is in turn calculated based on the images generated by processing the receiving echo. Each of the two images is present in a complex, digital form, i.e. it comprises a real part and an imaginary part, or equivalent: the amplitude and the phase.

[0029] Now, the phase difference φ follows up to a multiple of p (modulo p) through complex multiplication of the image points of the one image with the conjugated complex image points of the other image, and subsequent formation of the arctangent of the respective real and imaginary parts. The phase difference φ is obtained in this way, and by inserting φ in equation 1, ΔR is then obtained. $\begin{array}{cc}\Delta R=\frac{\lambda }{4\pi }\mathrm{\Delta \Phi }& \left(1\right)\end{array}$

[0030] The phase centers of the two receiving antennas A₁ and A₂ are removed by the length B, the so-called base line. The following results from the cosinus theorem and a few simple angle relations: $\begin{array}{cc}\cos \left(\theta \right)=\frac{{\left(R+\Delta R\right)}^2-R^2-B^2}{2·R·B}& \left(2\right)\end{array}$

[0031] After the sight angle θ has been calculated in equation (2), it is now possible to determine the relative altitude h as follows:
h=H−R·cos(θ)  (3)

[0032] The altitude h is actually not required in connection with the INROSAR system for representing the image points on the graphics display screen DS, but only the sight angle θ is used for calculating the coordinates of an impact point P on the integrated graphics display screen in the Rosar system. Furthermore, whether the angle of inclination of the antenna is known or not is unimportant as well because the representation on the display screen is only a relative representation of the image points with respect to the vertical line in relation to the base line B of the two antennas A₁ and A₂. The representation of the image is in fact dependent upon the position of the helicopter, for example due to the pitching. However, the antennas of the INROSAR-system and the center of the image are always in a fixed relation to each other. The altitude h and the angle of inclination α of the antennas are only required if a topographical chart with an absolute altitude H of the area over which the aircraft is passing is to be generated with the help of said INROSAR-system. The formulas specified above are useful also for a consideration of errors, as will be explained in the following.

[0033] The errors relevant to the INROSAR-system are the phase noise df and the change in the base line B between the phase centers of the antennas A₁ and A₂. The phase noise is composed of the sum of proportions of the different components. The greatest contributions are supplied by the transmitter, the receivers, the system timer and the noise of the A/D-converter. A typical order of magnitude for the entire phase noise δφ of an INROSAR-system amounts to approximately 5°. The change in the base line between the phase centers of the antennas A₁ and A₂ may be caused, for example by heating due to the incidence of sunlight rays. 0.001 m is assumed to be a typical value. The various influences result in a scatter δh of the altitude of the impact point P and thus in a scatter of the sight angle δθ. $\begin{array}{cc}\delta h=\frac{\lambda ·R}{4\pi ·B}\mathrm{\delta \varphi }& \left(4\right)\\ \delta h=-R·\tan \left(\theta \right)\frac{\delta B}{B}& \left(5\right)\end{array}$

[0034] This results in scatter of the sight angle dθ as follows: $\begin{array}{cc}\mathrm{\delta \theta }=\arcsin \left(\frac{\delta h}{R}\right)& \left(6\right)\end{array}$

[0035] In conjunction with an exemplified embodiment according to FIG. 1, the helicopter flies in the normal position, which means that the antennas A₁ and A₂ are positioned vertically one on top of the other. ΔR is determined based on equation (1). The value of the measured phase difference Δφ of the echo from the antennas A₁, and A₂ is ambiguous and can be determined only down to a value ranging between 0 and 2π. This ambiguity of 2π has to be determined by means of additional measurements. Thus, as shown in FIG. 4 an extra transmitter/receiver A₃ complementing the INROSAR conception comprising a transmitting/receiving antenna that has a narrow beam in elevation and covers the lower range of the sight angle can be used. The distance to the impact point on the ground can be clearly determined by the receive echo of this additional transmitting/receiving antenna. The INROSAR-system accepts the distance as a basic value and calculates the further ambiguities based on the rising distance from the continuous phase transitions. The following calculation example supplies the detailed explanations.

[0036] The calculation is based on the situation that the helicopter flies in its normal position. This means that the antennas A₁ and A₂ are vertically arranged one on top of the other.

The following parameters apply:

[0037]

Parameter	Meaning	Value 1/Value 2
H	flight altitude of INROSAR	100 m
R + ΔR	Distance between impact point P and	Example 1: 500.009 m
	antenna A1	Example 2:  500.09 m
R	Distance between impact point P and	500.00 m
	antenna A2
B	Base line between the phase centers	0.15 m
	of the antennas
δB	Error of length of base line B	0.001 m
δφ	Phase noise of the INROSAR system	5°
α	Angle in inclination of antennas A1	90° (vertically)
	and A2
λ	Radar wavelength	0.0090909

[0038] From equation (2) follows: $\begin{array}{cc}\theta =\arccos \left(\frac{{\left(R+\Delta R\right)}^2-R^2-B^2}{2·R·B}\right)& \left(7\right)\end{array}$

[0039] $\begin{array}{c}\mathrm{\theta 1}=\arccos \left(\frac{\left({500.009}^2-{500.000}^2-{0.15}^2\right)}{2·500.000·0.15}\right)\\ =\arccos \left(0.05985\right)\\ =86.57°\end{array}$ $\begin{array}{c}\mathrm{h1}=300-500.00·\cos \left(86.57°\right)\\ =70.09m\end{array}$

[0040] $\begin{array}{c}\mathrm{\theta 1}=\arccos \left(\frac{\left({500.09}^2-{500.00}^2-{0.15}^2\right)}{2·500.000·0.15}\right)\\ =\arccos \left(0.0599904\right)\\ =53.14°\end{array}$ $\begin{array}{c}\mathrm{h1}=300-500.00·\cos \left(53.14°\right)\\ =0.048m\end{array}$

[0041] From the equations (4) and (5) follows for the scatter δh of the altitude h of the impact point P: $\begin{array}{c}\delta h_{\mathrm{\delta \varphi }}=\frac{0.00909·500.00}{4·\pi ·0.15}\left(5^0/{57.3}^0\right)\\ =0.21m\end{array}$ $\mathrm{exactly}\text{:}0.210401168m$ $\begin{array}{c}\delta h_{\delta B}=-500.00·\tan \left({53.14}^0\right)·\frac{0.001}{0.15}\\ =4.45m\end{array}-\mathrm{based}\mathrm{on}\left(2\right)\text{:}=2.035-0.048m$
This results in a scatter of the sight angle δθ as follows:
Due to phase noise, δφ=5°: $\begin{array}{c}\mathrm{\delta \theta }=\arcsin \left(\frac{0.21}{500.00}\right)\\ =0.02°;\end{array}$
and because of errors in the length of base line B by δB=0.001 m: $\begin{array}{c}\mathrm{\delta \Phi }=\arcsin \frac{4.45}{500.00}\\ =0.5°.\end{array}$

[0042] FIG. 2 shows a block diagram of the exemplified embodiment shown in FIG. 1. This block diagram is equipped with the components required for the proposed interferometric radar method and requires no further explanations for the expert in the field.

[0043] However, the process is as follows, transmitting antenna T₁ emits a signal that then reflects back as multiple signals including at least two received signals into Antennas A₁ and A₂ as signals Se1 and Se2 respectively. These signals are then fed into the electronic components as shown in FIG. 2. In addition, data of the height H0 is also fed into these electronic components.

[0044] These signals S_e1 and S_e2 are fed or decoupled into demodulator 1 and 2 respectively wherein they are then converted from analog signals into a digital signals in Analog Digital converters 1 and 2 respectively. Next, the signals undergo a correction of movement errors via input from a kinematic sensor wherein this information is then fed into the respective memories 1 and 2. This information is then processed in a correlator with range curvature correction wherein it is then processed in an interferometry module and then printed to an associated display screen.

[0045] FIG. 3 shows a display of a rotating cross having an antenna lobe which is used to emit and receive signals. The antenna lobe can contain the two antennas and one transmitter as shown in FIG. 1 or it can contain the two transmitters and one antenna in FIG. 4 which also includes the additional sharply focused antenna which includes an additional transmitter and receiver.

[0046] As shown in FIG. 4 this additional sharply focused signal covers a lower region of the associated sight angle wherein as described above, this additional antenna A3 including an additional transmitter and receiver is disposed adjacent to the other antennas A₁ and A₂ and the associated other transmitter T₁. This additional antenna is for correlating and calculating the different ambiguities of the received signals in antennas A₁ and A₂. The distance to the impact point on the ground can be clearly determined by the receive echo of this additional transmitting/receiving antenna. The INROSAR-system accepts the distance as a basic value and calculates the further ambiguities based on the rising distance from the continuous phase transitions.

[0047] Finally, FIG. 5 shows the electronics design similar to that shown in FIG. 2 wherein this design includes the input from the additional antenna and receiver A3 for additional processing.

[0048] For example, a coherent radar signal (See FIG. 4) is transmitted by one antenna A₁ and then received by two coherent receiving antennas A₁ and A₂. The two coherent signals S_e1 and S_e2 are coherently demodulated in Demodulator 1 and Demodulator 2 and then converted from an analog signal to a digital signal in Analog to Digital converter (A/D) converter 1 and A/D converter 2 respectively. This information is then respectively fed into a system for the correction of movement errors. In addition, a kinematic sensor determines the movement of the helicopter (altitude and velocity) and feeds this information in the modules “Correction of Movement Errors 1 and 2”. These digital signals are then corrected and then written in Memory 1 and Memory 2.

[0049] A separate narrow beamed transmit/receive antenna supplies a signal S_e3 which is converted from analog to digital as well. In addition, in a distance module, an unambiguous distance D₀ to the nearest illuminated point on the ground is determined. The distance D₀ and the altitude H₀ are supplied to the Geometry Module. This information is forwarded to the Depth of Focus Module and Distance Intervals Module. Both modules in combination with the Reference Functions Module calculate current valid reference functions. These reference functions are written to Memory 3. The Range curvature comparator module determines the actual range curvature from the current helicopter movement and altitude. The reference functions from memory 3 and the current range curvature are fed in to the respective modules “Correlation with Range Curvature Correction 1 and 2” and corrected to reflect the current helicopter movement and altitude.

[0050] The radar signals from Memory 1 and Memory 2 are also fed into the respective modules “Correlation with Range Curvature Correction 1 and 2”. Here they are correlated with the reference functions to give two 2 dimensional images 1 and 2. The two respective images 1 and 2 are forwarded to the Interferometry Module, where the 3-dimensional image is calculated. The 3-dimensional image is transformed to a quasi-3 dimensional image so that it can be shown on a 2 dimensional display.

[0051] Thus, with this additional signal S_e3 is used to clarify any ambiguities associated with the measurement of the other two signals S_e1 and S_e2 that are read by antennas A₁ and A₂.

[0052] Accordingly, while at least one embodiment of the present invention has been shown and described, it is to be understood that many changes and modifications may be made thereunto without departing from the spirit and scope of the invention as defined in the appended claims.