Next Article in Journal
Wireless Local Area Networks Threat Detection Using 1D-CNN
Previous Article in Journal
Identification of Aluminothermic Reaction and Molten Aluminum Level through Vision System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Altimetry Method for an Interferometric Radar Altimeter Based on a Phase Quality Evaluation

1
Agency for Defense Development, Yuseong P.O. Box 35, Daejeon 34186, Republic of Korea
2
Department of Radio and Information Communication Engineering, Chungnam National University, 99 Daehak-ro, Yuseong-gu, Daejeon 34134, Republic of Korea
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(12), 5508; https://doi.org/10.3390/s23125508
Submission received: 21 April 2023 / Revised: 26 May 2023 / Accepted: 29 May 2023 / Published: 12 June 2023
(This article belongs to the Section Radar Sensors)

Abstract

:
A radar altimeter (RA) is useful to improve autonomous functions such as landing guidance or navigation control of an aircraft. To ensure more precise and safer flights by aircraft, an interferometric RA (IRA) capable of measuring the angle of a target is required. However, the phase-comparison monopulse (PCM) technique used in IRAs has a problem in that an angular ambiguity arises with respect to a target with multiple reflection points, such as terrain. In this paper, we propose an altimetry method for IRAs that reduces the angular ambiguity by evaluating the quality of the phase. The altimetry method as introduced here is sequentially described based on synthetic aperture radar, a delay/Doppler radar altimeter, and PCM techniques. Finally, a phase quality evaluation method is proposed for use in the azimuth estimation process. Aircraft captive flight test results are presented and analyzed, and the validity of the proposed method is examined.

1. Introduction

Terrain-referenced navigation (TRN) [1,2,3] is a navigation technology that estimates the current location on a map by comparing the measured altitude with the altitude from a digital elevation model (DEM). TRN has the advantage of being able to operate in the face of hostile interference to global positioning systems (GPS) but has the disadvantage of being vulnerable to altitude errors in the employed altimeter and DEM.
An altimeter is a key sensor that determines the performance of TRN. Depending on the method, barometric altimeters, laser altimeters, and radar altimeters can be used, among others. Radar altimeters [4,5,6,7] have the advantage of being able to operate even in poor weather conditions, unlike laser altimeters [8,9], because they have good permeability to clouds and fog. However, due to the wide beam width of the antenna, not only the altitude of the target directly below the beam but also the altitude of the target at the edge of the beam can be considered as that of the target directly below, resulting in altitude errors.
To minimize these errors, researchers have investigated radar altimeters for aircraft to measure not only the range of the target but also the angle of the target [10,11,12,13,14]. Among these altimeters, an interferometric radar altimeter (IRA) [10,11,12] estimates the target angle using the phase comparison monopulse (PCM) technique [15,16] with monopulse radar [17,18]. Unlike conventional radar altimeters that provide only one-dimensional range information, an IRA can provide three-dimensional information (range/elevation/azimuth) about the target and reduce altitude errors.
However, the PCM technique, which estimates the target angle by comparing the phase differences of the target signals between the receiving antennas, assumes a single target with no phase error. Therefore, if the phases of multiple target signals are overlapped or if the phase error increases for some reason, ambiguity arises in the angle estimation and the accuracy decreases [19,20,21]. This angular ambiguity often occurs in IRAs because these devices acquire reflected signals from the terrain with multiple reflectors.
In this paper, we propose an altimetry method that evaluates the phase quality of target signals and estimates the most reliable angle with the lowest level of error. In order to explain the proposed method, the process of deriving three-dimensional information (range/elevation/azimuth) of the target is described in Section 2. In Section 2.1, an altimetry concept for IRAs is introduced. In Section 2.2, the design of the transmitted signal to measure the relative range to the terrain based on radar theory is described. In order to derive the angle (=elevation) of the target, zero Doppler filtering [22] suitable for real-time processing on an aircraft is utilized, as described in Section 2.3.
In Section 2.4, before undertaking range fast Fourier transform (FFT), range cell migration compensation (RCMC) processing in the time domain is presented based on the synthetic aperture radar (SAR) theory [23]. In Section 2.5, an angle (=azimuth) estimation technique based on multiple-baseline interferometry is introduced. In Section 2.6, the limitations of the conventional angle estimation technique are analyzed. In order to overcome these limitations, the method for estimating the most reliable angle by means of a phase quality evaluation is proposed in Section 2.7. In Section 2.8, the validity of the proposed method is verified by analyzing actual aircraft test data. Finally, we conclude the paper in Section 3.

2. Proposed Method for IRAs

2.1. Altimetry Concept of IRAs

The operational concept of TRN with an IRA is shown in Figure 1. In this figure, TRN on a cruising aircraft uses an IRA to measure the vertical range ( D ) to the zero Doppler region within the antenna beam width. It then navigates by deriving the height of the terrain ( H ( = A D ) ) using the altitude from sea level ( A ) and by estimating its current position by comparing H to the height of the DEM [24,25].
However, as shown in Figure 2, the H of the terrain in the zero Doppler region ( P 1 P 3 ) varies as H 2 H 3 depending on the shape of the terrain. Accordingly, if the angle of the reflection point is not known, measurement error of the target range ( R ) occurs. For example, if the shortest range ( R 3 ) at point P 3 is regarded as the vertical range ( D 2 ), the measurement error is R 2 R 3 .
Given that the IRA can estimate the target angle (=azimuth) based on the PCM technique [15], estimating the angle ( θ ) with target range ( R ) can minimize the measurement error, as shown in Figure 2. For example, at point P 3 , the vertical range ( D 3 ) and the horizontal range ( Y 3 ) can be calculated based on R 3 and θ 3 , meaning that TRN can navigate by comparing the height of the DEM with H 3 at point P 3 .

2.2. Design of the Transmitted Signal for Computation of the Target Range

The IRA described in this paper transmits and receives linear frequency modulation (LFM) pulses to measure the altitude of distant terrain, as shown in Figure 1. The transmitted signal ( S t ) can be defined using Equation (1).
S t = A t exp j 2 π f c t + k e 2 t 2
In Equation (1), A t is the magnitude of the transmitted signal, f c is the carrier frequency, and t is the time. k e ( = B / T p ) is the slope of the frequency modulation, which varies depending on the frequency bandwidth ( B ) and the pulse width ( T p ). T p should be designed to be small to reduce the blind zone of the transmission time. However, because T p ( = 1 / B p ) is proportional to the signal-to-noise ratio (SNR) of the received signal, as shown in Equation (2) [26], it is necessary to analyze the SNR before designing the transmitted signal.
S N R = P t · G t · G r · λ 2 · σ · N P R F ( 4 π ) 3 · R 4 · K T 0 · F · B p · L T
In Equation (2), P t is the transmission power, G t is the transmission antenna gain, G r is the receiving antenna gain, λ is the wavelength, σ is the target radar cross section (RCS), N P R F is the pulse repetition frequency (PRF) number in one coherent processing interval (CPI), R is the range to the target, K is the Boltzmann constant, T 0 is the absolute temperature, F is the noise figure (NF), and L T is the total loss including the system loss and the propagation loss.
Figure 3 shows the geometry of the IRA’s reflected signal in the flight path of the cruising aircraft.
In Figure 3, T a ( = t a N / 2 ( t a N / 2 ) ) is the aperture time [23], which represents the time during one CPI, as expressed in Equation (3). The pulse repetition interval (PRI) is equal to 1/PRF, L is the length traveled during T a , V is the velocity of the aircraft, and Δ X is the resolution of the along track.
T a = L V = R λ 2 V · X = P R I · N P R F
If the resolution of the cross track is Δ Y and the inherent reflection coefficient of an object with an arbitrary material is σ 0 , σ can be represented by σ 0 · X · Y . Therefore, the SNR of the terrain signal in Equations (2) and (3) can be expressed as Equation (4). Considering the range resolution ( Δ R ( = c / 2 B ) ) and the angle of the terrain slope ( β ), Δ Y can be simplified to Δ R / sin β or it can be set to an arbitrary representative value.
S N R = P t · G t · G r · λ 2 · σ · Y · P R F ( 4 π R ) 3 · K T 0 · F · B p · L T · 2 V
By calculating a T p that satisfies the minimum SNR for detecting the terrain in Equation (4), the LFM pulse can be designed by combining B and PRF. However, considering the sampling frequency ( f s ) that satisfies the Nyquist sampling theorem, the altitude bin ( Δ H ) that can be processed without aliasing, i.e., without range ambiguity, can be obtained as shown in Equation (5). Thus, it is necessary to consider these parameters in the design stage of the transmitted signal.
Δ H = T p B · f s 2 · c 2

2.3. Employing Conventional Zero Doppler Filtering to Derive the Target Elevation

Regarding the acquisition of the three-dimensional information of the target, the target range computation described earlier is important. In addition to the acquisition of the azimuth, which will be described later, it is necessary to acquire the angle (=elevation; in the three-dimensional (3D) radar, the angle is categorized as azimuth and elevation; since the elevation is used in the DEM to indicate the height, it is expressed as the angle (=elevation) for clarity) of the flight direction. In this paper, instead of acquiring the elevation of the target, the zero Doppler filtering of the delay/Doppler radar altimeter (DDA) technique [22] is employed to derive the elevation region, and the RCMC is processed in the time domain.
Zero Doppler filtering in the DDA technique is a part of the range Doppler algorithm of the SAR technique [23]. It is performed only until the azimuth FFT; then, only the frequency signals in the zero Doppler region shown in Figure 1 are acquired. Compared to the range Doppler algorithm, it is less computationally intensive, which has the advantage of enabling real-time processing in aircraft flying at high speeds. In this regard, the concept of the target’s Doppler frequencies for the flight path is shown in Figure 4.
The Doppler frequency acquired by an aircraft flying at a constant altitude, as shown in Figure 4, is determined by Equation (6) [23]. Because f d is zero in the vertical direction of the aircraft, instead of acquiring the entire terrain signal, zero Doppler filtering acquires only the terrain signal at the zero Doppler point, such as f d 0 ( α = 0 ) , excluding the components f d k ( α 0 ) . In this way, the elevation region to be computed can be derived.
f d = 2 V λ · sin α
In order to perform zero Doppler filtering using the DDA technique, the azimuth FFT of the range Doppler algorithm should be processed beforehand. Therefore, it is necessary to compile a number of data in the along-track direction, which correspond to the amount of data collected during T a in Figure 3. If we plot the detected range ( R k ) for one target at each observation point ( t a k ) parallel to the range (=time t e ) axis while moving along the track, we can plot a range curve approximated by the red line shown in Figure 5. This is the range cell migration (RCM) process of SAR processing, requiring RCMC processing to compensate for it.
The received signal ( S r ) for the transmitted signal in Equation (1) is converted to a baseband signal and is determined by Equation (7). The conventional range Doppler algorithm performed the azimuth FFT on this signal and compensated for the range curve in Figure 5 through RCMC processing [23].
S r t = A r exp j 2 π f c τ t + k e 2 t τ t 2 = A r exp j 2 π f c τ 0 + k e 2 t 2 k e τ 0 t + k e 2 τ 0 2
In Equation (7), A r is the magnitude of the received signal, and τ ( t ) is the delay time of the reflected signal from the target, which can be expressed as Equation (8). In Equation (8), compared to f c , both f d and t are negligibly small. Accordingly, in this paper, we assume that τ ( t ) τ 0 .
τ t = 2 R t c = 2 R k c 2 V · sin α c · t = 2 R k c f d f c · t 2 R k c = τ 0

2.4. RCMC Processing in the Time Domain for Zero Doppler Filtering

The IRA undertakes deramping before converting the received signal ( S r ) to the baseband signal for real-time processing. Thus, S r can be calculated according to Equation (9) instead of Equation (7). In Equation (9), it is difficult to employ conventional RCMC processing as is.
S r t = A r exp j 2 π f c τ 0 + k e τ 0 t + k e 2 τ 0 2
In order to overcome this difficulty, in this paper, we present a method to compensate for the range curve by performing RCMC processing in the time domain before the azimuth FFT, although it is not as accurate as conventional RCMC processing performed in the frequency domain after the azimuth FFT.
The target range ( R k ) detected at each observation point ( t a k ) in Figure 5 can be expressed as Equation (10) [23], where R 0 is the closest range in Figure 5 and R k is the RCM component that forms the range curve indicated by the red line in Figure 5.
R k = R 0 2 + V 2 t a k t a 0 2 R 0 + R k = R 0 + V 2 ( t a k t a 0 ) 2 2 R 0
Decomposing Equation (9) into the azimuth time ( t a ) and the range time ( t e ) in Figure 5 and substituting 2 R k / c of Equation (8) into τ 0 of Equation (9), Equation (9) can be written as Equation (11).
S r t a , t e = A r exp j 4 π R k λ + j 4 π k e c R k t e j 4 π k e c 2 R k 2
Performing the range FFT operation on Equation (11) transforms it into Equation (12).
S r t a , t e = A r sinc f e k e 2 R k c · exp j 4 π R k λ j 4 π k e c 2 R k 2
In Equation (12), f e denotes the range frequency. The range (=time t e ) axis of the received signal was transformed into the frequency domain, and the center of the sinc function is k e 2 R k c .
Given that R k varies at every point ( t a k ) in Equation (10), R k is added to the closest range ( R 0 ) at every point ( t a k ) in Equation (12). Consequently, in this paper, we calculate R k for each observation point ( t a k ) and move the range cell by an amount equal to the calculated value, as shown in Figure 6. However, the proposed method has the advantage of real-time processing, but the accuracy is somewhat decreased because errors arise depending on the accuracy of R 0 and V and due to the error based on the size of R .
In Figure 6, the target signal aligned to the range cell by RCMC processing is the signal on the range axis and the time axis, where the range axis corresponds to the range time ( t e ) and the time axis corresponds to the azimuth time ( t a ) in Figure 5. By performing the azimuth FFT on this time axis, the target signal with target range τ t g t and Doppler frequency f d t g t can be obtained, as shown in Figure 7. The time axis is then converted to the frequency axis, and the target’s elevation region can be derived toward the direct downward direction by means of zero Doppler filtering, which detects only signals with a zero Doppler frequency in order to extract only signals for which f d t g t = 0 .

2.5. Conventional Angle-Estimation Technique Based on Multiple Baseline Interferometry

The computation of the target’s range and elevation to acquire three-dimensional information of the target was described earlier. In this section, we describe the estimation of the target’s azimuth (orthogonal to the direction of flight). The target angle estimation is based on the PCM technique of monopulse radar. However, considering that it is difficult to resolve the angular ambiguity with a typical single-baseline antenna, in this paper, we employ a method based on a multiple-baseline antenna to estimate the azimuth of the target by comparing the phase differences of the target signals received by each antenna, as in several earlier studies [14,27].
In Figure 8, Ant L (Left), Ant C (Center), and Ant R (Right) are placed with interantenna spacings of L 1 and L 2 , respectively, and the reflected signals from the target within the zero Doppler region are received by each antenna. Equation (13) can be derived from the geometry in Figure 8 for range R k and phase k of the received signal.
sin θ k = λ L k k 2 π + m k + L k 2 R k , ( k = 1,2 )
In Equation (13), L k is the interantenna spacing, and θ k is the angle of the target obtained from Ant L and Ant R. m k is a variable representing the periodicity of the measured phase difference k ( = k 0 ) , which has angular ambiguity that occurs at every period (2π) [28].
To resolve this angular ambiguity, the assumption that θ 0 is most reliable when the difference ( θ m i n ) between s i n θ 1 and s i n θ 2 is minimized is used to obtain m k , as shown in Equation (14). The obtained values of s i n θ 1 and s i n θ 2 are arithmetically averaged, as shown in Equation (15), to obtain the estimated angle ( θ ^ ) of θ 0 [18,27]. From the range ( R 0 ) and angle ( θ ^ ) obtained in this way, the vertical range ( D k ( = R 0 cos θ ^ ) ) in the direct downward direction can be estimated.
θ m i n = m i n m 1 , m 2 sin ( θ 1 ) s i n ( θ 2 ) , θ 1,2 θ w i d t h
θ ^ = asin s i n ( θ 1 ) + s i n ( θ 2 ) 2
For the multiple-baseline antenna designed in this paper, the values of θ k ( = asin sin θ k ) for k of Equation (13) is plotted in Figure 9. In Figure 9, the blue line represents the graph of 1 (= the phase difference between Ant R and C) versus θ 1 (when k = 1 ), and the red line represents the graph of 2 (=the phase difference between Ant L and Ant C) versus θ 2 (when k = 2 ).
If k , as indicated by the red point, is −0.4075 (when k = 1 ) and 0.2864 (when k = 2 ), then θ 1 and θ 2 are both 0.745°, resulting in θ m i n = 0 and θ ^ = 0.745°, respectively. If k , as indicated by the blue point, is 0.1107 (when k = 1 ) and 0.6491 (when k = 2 ), then θ 1 = −0.2° and θ 2 = 1.69°, resulting in θ m i n = 1.89° and θ ^ = 0.745°, respectively. Therefore, even with the added error in k , θ ^ corresponding to Equation (14) can be derived according to Equation (15).

2.6. Analysis of the Limitations of the Conventional Angle Estimation Technique

Before estimating the azimuth of the target, the peak point of the range cells within the zero Doppler filtering in Figure 7 is specified as target A with zero elevation, as shown in Figure 10. The azimuth is then estimated according to Equations (13)–(15) from the phase information of the target. In this regard, we analyzed and concluded that the conventional angle estimation technique has limitations in two aspects: the target identification process and the angle (=azimuth) estimation process.
First, regarding the target identification process, unlike monopulse radar, which detects a single target, an IRA receives target signals from multiple reflective points on the terrain, meaning that the peak in the range spectrum may be the signal reflected from a single target with a large RCS but may also be superimposed signals from multiple targets with similar ranges and different angles. Therefore, for a single target signal, the phase error may be small, and the estimated angle may be accurate, but for multiple target signals, the phase error increases, resulting in a large angle error.
Second, regarding the angle estimation process, it appears to be reasonable to derive Equations (14) and (15) based on the assumption that θ ^ is most reliable when the difference between s i n θ 1 and s i n θ 2 is minimized for a multiple-baseline antenna, but it is difficult to determine the extent to which the result of Equation (15) can be trusted. In this regard, Figure 11 is instructive.
In Figure 11, when 1 is 0.2752 and 2 is 0.764, θ m i n is 2.49° in regions 1 and 2, yielding θ ^ = 0.745° in region 1 and θ ^ = −13.295° in region 2. In this case, the angle cannot be estimated using the conventional method because there are two solutions. These two limitations motivated us to consider the phase quality evaluation in this paper.

2.7. Azimuth Estimation Based on a Phase Quality Evaluation

In order to overcome the two limitations of the conventional angle estimation techniques described above, we propose a method that evaluates the phase quality of the target signals in the range spectrum, as shown in Figure 10, before estimating the target angle (=azimuth), after which it selects the target whose azimuth is estimated to be the most reliable. In other words, instead of “calculate signal magnitude → select target → estimate azimuth”, the azimuth of the proposed method is estimated in the order of “calculate signal magnitude → select target group → calculate phase quality and estimate azimuth in advance → evaluate phase quality → select the most reliable target → finally estimate azimuth”.
This method takes into account the contamination of the phase information. Just as Equations (14) and (15) in Section 2.6 are derived from the assumption that θ ^ of the conventional multiple-baseline antenna is most reliable when θ m i n is minimized, the proposed method works on the assumption that the better the phase quality is, the more reliable θ ^ is.
First, we set a threshold larger than the noise level, as shown in Figure 10. We classify the signals in a certain region of the range spectrum (e.g., N/2 to N when the total number of samples in the range spectrum is N) into a group of targets. The azimuths of these targets are then estimated according to Equations (14) and (15), and their phase qualities are evaluated as described below.
In this paper, we precalculate and store the phase difference versus the angle (=azimuth) values in Figure 9 because the phase quality could not be determined by the conventional Equations (14) and (15). For example, if the target azimuth is 0.745°, −0.4075 for φ 1 and 0.2864 for φ 2 are then calculated and stored in the database (DB) ( φ k is referred to as the stored value to distinguish it from the measured value ( k )). Then, from the range spectrum in Figure 10 obtained from the aircraft captive flight test (CFT), the target azimuths of the entire target group are calculated. If the target azimuth ( θ ^ ) is determined to be 0.745°, 1 is equal to 0.1107, and 2 is equal to 0.6491 in Equation (13), the quality of phase q is then calculated as shown in Equation (16).
q = φ 1 1 + φ 2 2
Subsequently, by comparing the phase quality ( q ) of the entire target group, as shown in Equation (17), selecting the mth sample with the minimum value as the target, and estimating the target azimuth ( θ ^ ), the proposed method completes the target azimuth estimation.
q m i n = m i n m q m , m = m 1 , . . , m k , . . , M
The proposed method is adopted for use with targets 1 and 2 in Figure 11. For target 1, θ ^ is 0.745°, 1 is 0.2752, 2 is 0.764, φ 1 is −0.4075, and φ 2 is 0.2864. Thus, according to Equation (16), q is calculated and found to be 1.1603. For target 2, θ ^ is −13.295°, 1 is 0.2752, 2 is 0.764, φ 1 is 0.9424, and φ 2 is 1.227. Then, q is calculated according to Equation (16) and found to be 1.1302. Therefore, for this sample, −13.295° with the smaller q ( = 1.1302 ) is estimated as the target azimuth ( θ ^ ). The subsequent procedure is as described earlier. Figure 12 shows a flow chart of the proposed method.

2.8. Aircraft CFT Results

The IRA was mounted on the underside of an aircraft, as shown in Figure 13, to conduct an aircraft CFT and acquire test data. The test was conducted in an area near Jeollabuk-do, South Korea, and the plane flew along the route shown in Figure 14, with a flight altitude of about 1.5 km above sea level and a flight speed of about 170 knots.
The processing results of the kth scan of the flight test data are shown in Figure 15, Figure 16 and Figure 17 below. The height of the DEM in the orthogonal direction of the flight direction is shown in Figure 15. In this figure, the green line indicates the height of the DEM of the detectable region within the antenna beam width of the IRA, the blue dot is the height of the terrain computed by the conventional method, and the red dot is the height of the terrain computed by the proposed method.
The range spectrum in the zero Doppler region obtained by RCMC processing in the time domain, azimuth FFT, and zero Doppler filtering for the above terrain is shown in Figure 16. This range spectrum is identical to the range spectrum in Figure 10, but in Figure 16, we can see the target (red dot) selected as a result of the phase quality evaluation.
The results of the azimuth estimation of the signals above the threshold in Figure 16 are shown in Figure 17. In Figure 17, the green line shows the target angle (=azimuth) on the DEM, and the black dots show the target angles estimated by comparing the phase quality. Figure 17 shows that there is some similarity between the estimated angles and the angles of the DEM, but there are significant angular errors in that the black dots are far off the green line.
This may be an angular ambiguity inherent in the multiple-baseline antenna-based PCM technique [14], an inherent error in the DEM due to the accuracy of DEM production, or a phase error due to the superposition of the phases of multiple targets in one range cell ( Δ R ).
In this paper, the phase quality containing these errors is compared using Equation (16) to calculate the angle at each point, and the angle ( θ ^ ) is estimated as the target angle by evaluating the phase quality ( q ) according to Equation (17) and selecting the mth sample with the minimum value of q as the target. In Figure 17, the red dot represents the target angle estimated by the proposed method, and the blue dot represents the target angle at the maximum signal point according to the conventional method.

3. Discussion

In this paper, an altitude calculation for one scan was performed for a total of 11,541 scans with the proposed method, and the error was analyzed by comparing this outcome with the DEM results. Figure 18 and Figure 19 show the altimetry results when RCMC processing in the time domain and zero Doppler filtering as described in Section 2.4 are performed, with the azimuth of the maximum signal point assumed to be 0 degrees without the angle estimation. This process is similar to that of the (non-interferometric) radar altimeter but is more accurate because it utilizes RCMC processing and zero Doppler filtering.
Figure 18 shows the altimetry error for each scan plotted against the flight time axis, and Figure 19 shows the result as a histogram. We calculate the distance between the 3D coordinates of the red dots labelled as targets and the 3D coordinates of the light green lines representing the DEM in Figure 17 and take the smallest of these values as the altimetry error. The mean of the error is 13.55 m, and the standard deviation of the error is 27.87 m. Representing this as a Gaussian distribution results in the red graph shown in Figure 19. The Gaussian distribution is obtained using the least-squares approximation of the histogram.
Figure 20 and Figure 21 show the altimetry results of the IRA based on the azimuth estimation of the maximum signal point. Considering the 3 dB beam width of the antenna, it is assumed that the estimated angle is valid if it is within θ w i d t h (see Equation (14)). As a result, the valid rate of the angle is 50.68%, the mean of the error is 12.35 m, and the standard deviation of the error is 25.00 m. Representing this as a Gaussian distribution results in the red graph in Figure 21.
As shown in Figure 21, the histogram is significantly different from the Gaussian distribution because, as previously discussed in Section 2.6, the peak could be the signal reflected from a single target with a large RCS but could also be a superposition of signals from multiple targets with similar ranges and different angles. In other words, the deviation between a signal reflected from a single target and the superposition of signals from multiple targets is large such that the histogram appears to differ considerably from the calculated Gaussian distribution.
The altimetry results of the proposed method used to estimate the azimuth by evaluating the phase quality are shown in Figure 22 and Figure 23. Here, the angle valid rate is 52.95%, the mean of the error is 12.22 m, and the standard deviation of the error is calculated to be 16.41 m. This is represented by a Gaussian distribution, as represented by the red graph in Figure 23.
As shown in Figure 23, the histogram based on the proposed method appears to be similar to a Gaussian distribution, unlike that in Figure 21, and Figure 22 shows fewer extreme outliers and a smaller standard deviation of the error, unlike the outcomes in Figure 20. Table 1 summarizes the results for the three cases above.
In Table 1, the proposed method improves the standard deviation of the IRA. The circular error probability (CEP) of TRN depends on the accuracy of IRA. Thus, the smaller the standard deviation of the IRA, the more accurate the CEP of TRN. On the other hand, the mean values for the proposed and conventional methods are similar. There could be many reasons for this result. One of the reasons is that there may be some offset due to seasonal vegetation when the DEM is constructed and when the IRA measures the terrain.
The proposed method has the advantage of reducing the standard deviation but has the disadvantage of increasing the size of data due to the complexity of the processing procedure compared to the conventional method. However, the added procedures are all after the azimuth FFT, so the volume of data does not increase significantly.

4. Conclusions

In this paper, we proposed an altimetry method for IRA by evaluating the phase quality. Based on the SNR analysis between the terrain and IRA on the cruising aircraft, the LFM pulse is designed to obtain high-resolution range information, the RCMC processing in the time domain is designed based on the zero Doppler filtering of the DDA to obtain the elevation information of the target in the direct downward direction, and the PCM technique based on the conventional multiple baseline interferometry is employed to obtain the azimuth information using the phase difference of the received signal.
Then, the azimuths of a group of targets are computed by comparing the quality of the obtained phases. Finally, the azimuth information is estimated by selecting the most reliable target by evaluating their phase qualities. The test results verified that the proposed method can accurately calculate the altitude of the IRA. In the analysis process of the test results, it was found that the proposed method exhibits smaller standard deviation characteristics compared to the conventional method.

Author Contributions

Conceptualization, J.-S.H.; methodology, J.-S.H. and S.-Y.H.; data curation, J.-S.H.; writing—original draft preparation, J.-S.H.; writing—review and editing, J.-S.H. and S.-Y.H.; visualization, J.-S.H.; supervision, S.-Y.H.; project administration, J.-S.H. and S.-Y.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Carreño, S.; Wilson, P.; Ridao, P.; Petillot, Y. A survey on terrain based navigation for AUVs. In Proceedings of the OCEANS 2010 MTS/IEEE SEATTLE, Seattle, WA, USA, 20–23 September 2010; pp. 1–7. [Google Scholar]
  2. Zhang, Q.Q.; Zhao, L.; Zhou, J.H.; Wang, B.H.; Zhang, Y.P. A real-time airborne terrain aided inertial navigation system and its performance analysis. Adv. Space Res. 2017, 60, 2751–2762. [Google Scholar] [CrossRef]
  3. Sung, C.K.; Lee, S.J. Moment matched Gaussian kernel and region representative likelihood for performance improvement of PMF-based TRN. Int. J. Control Autom. Syst. 2020, 18, 1691–1704. [Google Scholar] [CrossRef]
  4. BiNgöl, H.E.; Akin, B.; KOC, Ö. Radar altimeter as a navigation aid using hierarchical elevation map clustering. In Proceedings of the 2012 IEEE/ION Position, Location and Navigation Symposium, Myrtle Beach, SC, USA, 23–26 April 2012; pp. 377–381. [Google Scholar]
  5. Choi, J.H.; Jang, J.H.; Roh, J.E. Optical-delay-based wide-range FMCW radar altimeter. Electron. Lett. 2014, 50, 1975–1977. [Google Scholar] [CrossRef]
  6. Choi, J.H.; Jang, J.H.; Roh, J.E. Design of an FMCW Radar Altimeter for Wide-Range and Low Measurement Error. IEEE Trans. Instrum. Meas. 2015, 64, 3517–3525. [Google Scholar] [CrossRef]
  7. Videmsek, A.; Haag, M.U.; Bleakley, T. Radar Altimeter Aiding of GNSS for Precision Approach and Landing of RPA. In Proceedings of the 2019 Integrated Communications, Navigation and Surveillance Conference (ICNS), Herndon, VA, USA, 9–10 April 2019; pp. 1–16. [Google Scholar]
  8. Amzajerdian, F.; Pierrottet, D.; Petway, L.; Hines, G.; Roback, V. Lidar system for precision navigation and safe landing on planetary bodies. In Proceedings of the International Symposium on Photoelectronic Detection and Imaging 2011: Laser Sensing and Imaging; and Biological and Medical Applications of Photonics Sensing and Imaging, Beijing, China, 24–26 May 2011; pp. 27–33. [Google Scholar]
  9. Han, K.J.; Sung, C.K.; Yu, M.J.; Park, C.G. Performance improvement of laser TRN using particle filter incorporating Gaussian mixture noise model and novel fault detection algorithm. In Proceedings of the 31st International Technical Meeting of the Satellite Division of the Institute of Navigation, Miami, FL, USA, 24–28 September 2018; pp. 3317–3326. [Google Scholar]
  10. Park, I.C.; Lee, S.I.; Chun, J.H.; Lee, H.J.; Jang, J.H. Simulation study of altitude and angle estimation with an InSAR altimeter. J. Korean Inst. Electromagn. Eng. Sci. 2014, 25, 838–848. [Google Scholar]
  11. Lee, J.; Bang, H. Radial basis function network-based available measurement classification of interferometric radar altimeter for terrain-aided navigation. IET Radar Sonar Navig. 2018, 12, 920–930. [Google Scholar] [CrossRef]
  12. Kim, Y.; Park, J.; Bang, H. Terrain referenced navigation using an interferometric radar altimeter. J. Inst. Navig. 2018, 65, 157–167. [Google Scholar] [CrossRef]
  13. Oh, J.; Sung, C.K.; Lee, J.; Lee, S.W.; Lee, S.J.; Yu, M.J. Accurate measurement calculation method for interferometric radar altimeter-based terrain referenced navigation. Sensors 2019, 19, 1688. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Ha, J.S.; Jang, J.H.; Lee, S.J.; Park, Y.S.; Hong, S.Y. Angle Ambiguity Analysis of Antenna Baseline for Interferometric Radar Altimeter. J. Korean Inst. Electromagn. Sci. 2021, 32, 626–634. [Google Scholar] [CrossRef]
  15. Sherman, S.M.; Barton, D.K. Monopulse Principles and Techniques, 2nd ed.; Artech House: Boston, MA, USA, 2011. [Google Scholar]
  16. Zhang, Y.X.; Liu, Q.F.; Hong, R.J.; Pan, P.P.; Deng, Z.M. A Novel Monopulse Angle Estimation Method for Wideband LFM Radars. Sensors 2016, 16, 817. [Google Scholar] [CrossRef] [PubMed]
  17. Mahafza, B.R. Radar Systems Analysis and Design Using Matlab, 3rd ed.; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  18. Jensen, J. Angle measurement with a phase monopulse radar altimeter. IEEE Trans. Antennas Propag. 1999, 47, 715–724. [Google Scholar] [CrossRef]
  19. Blair, W.D.; Brandt-Pearce, M. Monopulse DOA estimation of two unresolved Rayleigh targets. IEEE Trans. Aerosp. Electron. Syst. 2001, 37, 452–469. [Google Scholar] [CrossRef]
  20. Lee, S.P.; Cho, B.L.; Kim, Y.S. Unambiguous Angle Estimation of Unresolved Targets in Monopulse Radar. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 1170–1177. [Google Scholar] [CrossRef]
  21. Lee, S.H.; Lee, S.J.; Choi, I.O.; Kim, K.T. ICA-based phase-comparison monopulse technique for accurate angle estimation of multiple targets. IET Radar Sonar Navig. 2018, 12, 323–331. [Google Scholar] [CrossRef]
  22. Raney, R. The Delay/Doppler radar altimeter. IEEE Trans. Geosci. Remote Sens. 1998, 36, 1578–1588. [Google Scholar] [CrossRef]
  23. Cumming, I.G.; Wong, F.H. Digital Processing of Synthetic Aperture Radar Data; Artech House: Norwood, MA, USA, 2005. [Google Scholar]
  24. Nam, S.H.; Sung, C.K.; Lee, S.W. 3rd Order RB-PMF based Terrain Aided Navigation Algorithm considering Altitude Error. In Proceedings of the Korean Society for Aeronautical and Space Sciences Fall Conference, Jeju, Republic of Korea, 28 November–1 December 2018; pp. 291–292. [Google Scholar]
  25. Sung, C.K.; Lee, S.J. Performance Comparisons of Terrain Referenced Navigations Based on PMF and RB-PMF Considering Altitude Error. In Proceedings of the Korean Society for Aeronautical and Space Sciences Fall Conference, Jeju, Republic of Korea, 28 November–1 December 2018; pp. 974–975. [Google Scholar]
  26. Skolnik, M.I. Introduction to Radar Systems, 3rd ed.; McGraw-Hill: New York, NY, USA, 2001. [Google Scholar]
  27. Sim, H.K.; Jung, M.A.; Kim, S.C. A performance analysis of phase comparison monopulse algorithm for antenna spacing and antenna array. J. Korean Inst. Commun. Inf. Sci. 2015, 40, 1413–1419. [Google Scholar] [CrossRef] [Green Version]
  28. Lipsky, S.E. Microwave Passive Direction Finding; John Wiley and Sons: New York, NY, USA, 1987. [Google Scholar]
Figure 1. Operational concept of TRN using an IRA.
Figure 1. Operational concept of TRN using an IRA.
Sensors 23 05508 g001
Figure 2. Geometry of altitude error in terrain.
Figure 2. Geometry of altitude error in terrain.
Sensors 23 05508 g002
Figure 3. Geometry of the IRA in the flight path.
Figure 3. Geometry of the IRA in the flight path.
Sensors 23 05508 g003
Figure 4. Doppler frequency of the target for the flight path.
Figure 4. Doppler frequency of the target for the flight path.
Sensors 23 05508 g004
Figure 5. Range curve of the target for the flight path.
Figure 5. Range curve of the target for the flight path.
Sensors 23 05508 g005
Figure 6. RCMC processing in the time domain.
Figure 6. RCMC processing in the time domain.
Sensors 23 05508 g006
Figure 7. Azimuth FFT and zero Doppler filtering.
Figure 7. Azimuth FFT and zero Doppler filtering.
Sensors 23 05508 g007
Figure 8. Geometry of multiple-baseline interferometry.
Figure 8. Geometry of multiple-baseline interferometry.
Sensors 23 05508 g008
Figure 9. Concept of angle estimation based on the phase difference.
Figure 9. Concept of angle estimation based on the phase difference.
Sensors 23 05508 g009
Figure 10. Range spectrum in the zero Doppler region.
Figure 10. Range spectrum in the zero Doppler region.
Sensors 23 05508 g010
Figure 11. Two targets with the same angle difference.
Figure 11. Two targets with the same angle difference.
Sensors 23 05508 g011
Figure 12. Flow chart of the proposed method.
Figure 12. Flow chart of the proposed method.
Sensors 23 05508 g012
Figure 13. IRA and aircraft.
Figure 13. IRA and aircraft.
Sensors 23 05508 g013
Figure 14. Test area and corresponding path.
Figure 14. Test area and corresponding path.
Sensors 23 05508 g014
Figure 15. Test result (terrain height in the cross track).
Figure 15. Test result (terrain height in the cross track).
Sensors 23 05508 g015
Figure 16. Test result (target range in the zero Doppler region).
Figure 16. Test result (target range in the zero Doppler region).
Sensors 23 05508 g016
Figure 17. Test result (target azimuth in the zero Doppler region).
Figure 17. Test result (target azimuth in the zero Doppler region).
Sensors 23 05508 g017
Figure 18. RA error result (max point).
Figure 18. RA error result (max point).
Sensors 23 05508 g018
Figure 19. RA error histogram (max point).
Figure 19. RA error histogram (max point).
Sensors 23 05508 g019
Figure 20. IRA error result (max point).
Figure 20. IRA error result (max point).
Sensors 23 05508 g020
Figure 21. IRA error histogram (max point).
Figure 21. IRA error histogram (max point).
Sensors 23 05508 g021
Figure 22. IRA error result (proposed method).
Figure 22. IRA error result (proposed method).
Sensors 23 05508 g022
Figure 23. IRA error histogram (proposed method).
Figure 23. IRA error histogram (proposed method).
Sensors 23 05508 g023
Table 1. Comparison of test results by the algorithms.
Table 1. Comparison of test results by the algorithms.
MethodsAngle Valid Rate (%)Mean (m)Std (m)
RA (max point)-13.5527.87
IRA (max point)50.6812.3525.00
IRA (proposed)52.9512.2216.41
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ha, J.-S.; Hong, S.-Y. Altimetry Method for an Interferometric Radar Altimeter Based on a Phase Quality Evaluation. Sensors 2023, 23, 5508. https://doi.org/10.3390/s23125508

AMA Style

Ha J-S, Hong S-Y. Altimetry Method for an Interferometric Radar Altimeter Based on a Phase Quality Evaluation. Sensors. 2023; 23(12):5508. https://doi.org/10.3390/s23125508

Chicago/Turabian Style

Ha, Jong-Soo, and Sung-Yong Hong. 2023. "Altimetry Method for an Interferometric Radar Altimeter Based on a Phase Quality Evaluation" Sensors 23, no. 12: 5508. https://doi.org/10.3390/s23125508

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop