=Paper= {{Paper |id=Vol-1909/paper13 |storemode=property |title=The Detection of Lengthened Objects by Pulse Altimeter |pdfUrl=https://ceur-ws.org/Vol-1909/paper13.pdf |volume=Vol-1909 |authors=Artem K. Sorokin,Vladimir G. Vazhenin }} ==The Detection of Lengthened Objects by Pulse Altimeter== https://ceur-ws.org/Vol-1909/paper13.pdf
             The Detection of Lengthened Objects by Pulse Altimeter

                                          Artem K. Sorokin, Vladimir G. Vazhenin
                                                 Ural Federal University
                                              Yekaterinburg, Russia, 620004
                                                    sorokinak@urfu.ru



                                                          Abstract
                       The algorithm of lengthened objects detection is based on statistical
                       processing of reflected pulse altimeter signal. The method of maximum
                       posteriori probability is applied to make a decision. Minimum of maximum
                       posteriori probability is accorded with lengthened objects position. This fact
                       is taken into account to detect the objects’ position. This algorithm allows
                       autonomic navigation to be implemented. Two cases of the algorithm
                       implementation are shown in this article. The first method is based on the
                       unsharpened beam regime of radar altimeter and the second is based on the
                       application of the Doppler’s filtration. The results of flight experiments
                       proving the correctness of the created algorithm are introduced in this
                       article.



1    Introduction
The necessity of robustness navigation system design for the unmanned airborne vehicles (UVA) is demonstrated in the
article [1]. It is applicable for the relief navigation system supplement, if the relief of the terrain is mild. While examining
this fact, the onboard pulse radar altimeter [2] with wave length of the carrier frequency in C-band as the sensor was
applied. This algorithm was designed for a nearby horizontal flight where the evolutions are absent.
    While solving the problem of lengthened objects (LO) detection, the properties of different terrain types were
explored and the patterns of the reflected signals amplitude distributions were obtained. The terrains were classified by
the width of the backscattering diagram and the module of the reflection coefficient so that allows to their identification
to be done, as it is made in [1]. On the base of [1, 2, 4] the amplitude distributions of the reflected signals, related to the
known terrains were obtained. The term “the pattern” was introduced. The given distributions allow us to estimate the
compliance rate between patterns and current distributions if they are compared during the flight and model experiments.
The comparison is corresponded to the estimation of the pattern’s crossed square and current distribution (see Fig. 1). As
the result of the square computation, it is possible to obtain the probability of the correct detection and the minimum of
the maximum posteriori probability function corresponding to the border between two changing types of the terrain.
    While modeling, a facet phenomenological model of the reflected signal was applied. It was described in [1].
In this paper the points from [1] are enlarged and applied to LO of the stripe type, the accuracy of the border (between
two typical terrains) detection is also evaluated.
                                                    2          2
                                           U Usin Ucos




                               Figure 1: The algorithm of lengthened objects (LO) detection

2    The Typical Regime of Radar Altimeter
The model experiment is organized according to the Fig. 2. as follows: UVA is moving from the “terrain I” across the
“terrain II” to the “terrain I”, here the “terrain II” is LO. In Fig. 2 it is marked the next: D – the LO width, L – the
exposure spot width, Θ – the LO’s orientation angle. As the result of signal accumulation for each, the UVA’s position is
identified (there were more than 10 000 envelope counts per the UVA’s position), used for the probability estimation
evaluation.




                            Figure 2: The Model experiment of linear reference (LO) detection

    The combinations of typical terrains, which are applicable for LO detection are described in [1].
    The results of the model experiment in the case of “forest/asphalt” are shown in Fig. 3.
    According to the Fig. 3, the real borders are shifted symmetrically to the borders detected by the algorithm. It is
explained by the minimum shift of posteriori probability function to the side of the less reflective terrain. The minima of
the posterior probability detect the border between typical terrains which can be evaluated by the shift error (in Fig. 3 the
shift error is about 15% of the exposure spot width), and the minimum probability of the correct detection is 0,8.
The results of the model experiments for LO detection are introduced in the Table 1 for the case with Θ=[30˚;90˚] and
different signal/noise ratio.
      Figure 3: The results of model experiment for LO, with its width 0,8 of the exposure spot, terrain combination
                                                     “asphalt/forest”

    The Table 1 makes it possible to choose the terrain’s combinations which could be distinguished during the flight by
the pulse altimeter. As the result, the best LO for the algorithm is the “water” surface with the background of the
scattering terrains such as “meadow” or “forest”. This result is similar to the point in [1]. In addition, the real width of the
LO is closely connected to the flight’s height and the narrow objects can be detected at low altitude.

                Table 1: The ability of the lengthened objects detection under the additive noise influence

                      SN=0дБ                   SN=10dB             SN=20dB       SN=30dB
                      D=0,1L      D=0,8L       D=0,1L D=0,8L       D=0,1L D=0,8L D=0,1L D=0,8L
                      FAF         FAF          FAF    FAF          FAF    FAF    FAF    FAF
                      FWF         FWF          FWF    FWF          FWF    FWF    FWF    FWF
                      GWG         GWG                 GWG                 GWG           GWG
                      BAB         BAB                 BAB
                      BWB                      BWB BWB             BWB       BWB       BWB       BWB
                      MAM         MAM                 MAM                    MAM                 MAM
                      MWM         MWM          MWM MWM             MWM       MWM       MWM       MWM
                      SWS         SWS                 SWS                    SWS                 SWS
                                  ABA
                                  WGW                    WGW                 WGW                 WGW
                                  WMW                    WMW                 WMW                 WMW
                                  WSW                    WSW                 WSW                 WSW
                                                         AFA                 AFA                 AFA
                                                         WFW                 WFW                 WFW
                                                                             AGA
                                                         AMA                 AMA                 AMA
                                                         GAG                 GAG                 GAG

    In the Table 1 following abbreviations are used: SN – signal/noise ratio, А – asphalt, B – bushes, C – concrete, F –
forest, G – grass, M – meadow, N – ground, S – snow, W – water. For example, the abbreviation FAF is equivalent to the
combination “Forest-Asphalt-Forest”. Shadowed combinations with high false detection probability of the LO are also
presented in the table. It can happen because of signal and noise nature, being not the matter of importance for the
algorithm it does not matter. These cases have to be excluded from the study.
   Table 1 can be used to choose the LO, for designing recommendations and while choosing the UVA route.

3    The Doppler Filtration
The application of the Doppler filtration (DF) as the method of increasing the algorithm’s spatial resolution for solving
the problem of LO detection is shown in this part of the reserch.
    According to [2], we can observe the curves of equivalent Doppler frequencies on the terrain, which appearance can
be explained by equivalent radial velocities of the UVA in relation to the terrain. We need to find the conditions of DF
application which are necessary for supplying the maximum posteriori probability.
    Three series of experiments were conducted for the test. In the first range of the experiments, the orientation of the
Doppler filter (see fig. 4) was changed and the probability of the correct discrimination of the terrain was estimated [1].




                 Figure 4: Model experiment to find the best conditions of Doppler filtration application

    In Fig. 5 the results of the model experiments are shown. It was obtained that the maximum posteriori probability
could be gained according to the position in Fig. 4, which is marked with Δ1.
    Similarly, the other experiment was conducted. The width of the Doppler’s filter’s stripe was changed to the optimal
angle, which was obtained in the previous experiment. As the result, it was achieved that the minimal d (0,2L) provides
the maximum posteriori probability. Here is d – the Doppler stripe width; L – the exposure spot width.




                 Figure 5: Dependence Pcor.det.(α) for “water” combination and some background terrains
    So, during the test it was obtained that the DF provides better discrimination for terrains if the width of the Doppler
stripe is small (d≤0,2L) and the direction of the exposure is corresponded to nadir direction.
    The application of the Doppler filtration provides the increasing of the spatial resolution for the radar altimeter. It
allows us to detect narrow LO, but simultaneously it can prevent us from the detection if the angle φ in Fig. 6 is small.




                          Figure 6: The model experiment m(φ) for two ratios between d and D

    The third experiment was about LO’s orientation. Fig. 6 shows the relative part of the LO in the Doppler’s filter’s
stripe (m) from the angle of the LO (φ)orientation . The model experiment showed that the maximum influence in the
reflected signal delivered by the LO is oriented co directional to the Doppler stripe. If the angle is φ=30˚, the part of the
LO decreases by half. As the result we suggested that the advantages of the DF are applicable only if φ≥30˚. In this case
maximum posteriori probability gained if φ=90˚.
    The similar to [1] experiment was conducted during the exploration and the results are shown in Table 2. The results
of the experiment similar to the model experiment in [1] are presented in Table II. Here are the number of discriminated
terrains if the signal/noise ratio (sn) and the width of the LO are changed (D).

                       Table 2: The ability of lo detection if the normal noise influences are added

                          SN=0дБ        SN=10dB       SN=20dB       SN=30dB
                          D=0,1L D=0,8L D=0,1L D=0,8L D=0,1L D=0,8L D=0,1L D=0,8L
                          FAF, FAF, FAF,       FAF, FAF, FAF, FAF, FAF,
                          FWF, FWF, FWF, FWF, FWF, FWF, FWF, FWF,
                          GWG, GWG, GWG, GWG, GWG, GWG, GWG, GWG,
                          BWB, BWB, BAB,       BAB, BAB, BAB, BAB, BAB,
                          MWM, MWM, BWB, BWB, BWB, BWB, BWB, BWB,
                          SWS SWS, MAM, MAM, MAM, MAM, MAM, MAM,
                                 WGW, MWM, MWM, MWM, MWM, MWM, MWM,
                                 WMW, SWS, SWS, SWS, SWS, SWS, SWS,
                                 WSW, GAG      ABA, GAG ABA, GAG, ABA,
                                 AFA           WGW,          WGW, SAS      WGW,
                                               WMW,          WMW,          WMW,
                                               WSW,          WSW,          WSW,
                                               WFW,          WFW,          WFW,
                                               AGA,          AGA,          AGA,
                                               AMA,          AMA,          AMA,
                                               GAG,          GAG,          GAG,
                                               WBW           WBW,          WBW,
                                                             ASA           SAS,
                                                                           ASA
    The information in the Table 2 provides the detection of LO with the probability of correct discrimination higher than
0,6 if φ is changed from 30˚ up to 90˚. As the result, it is shown that the change of the signal/noise ratio leads to the
change of typical LO number which can be detected by the algorithm (decreasing leads to the increasing of the LOs), also
the increasing of the LO width leads to the LOs number increasing.

4    The Flight Experiments
The information from the experiments [4] was used during investigations. Two flights were made. The flights had similar
tracks in three heights: 60m, 100m and 300m. The weather conditions of flights were different, in pattern’s flight it was
rainy and in the second flight it was dry. The results of the experiment were processed by the created algorithm from [1].
    To conduct the flight experiment, the information was obtained from pulse radar altimeter [2] synchronized with the
information from satellite navigation system (GPS) and with video record (see Fig. 8). So, we had the ability to fix the
change of the terrain with the accuracy of about a few meters.
    The comparison of the current model experiments’ results with the typical regime for pulse radar altimeter shows that
DF allows us to distinguish approximately twice typical LOs than typical regime.




                                        Figure 7: The scheme of flight experiment

    The program in MATLAB’s script was also designed. The program makes it possible to accomplish the process of
flight experiments according to the created algorithm.
    Furthermore, the information was processed by applying the multistage scheme. At the first stage the information
from the pulse altimeter was prepared (the pulse train was formed, computed quadrature with considering the amplifying
and many flags). At the next stage the histograms of the reflected signals were formed (see Fig. 7). In the pattern
experiment the histograms were identified according to the video record. For the processing the segments of homogenous
typical terrain were chosen. There were no any significant evolutions of the UVA. These segments for different heights
were normalized and after averaging, the histograms for typical terrains were build. As a result, the pattern histograms for
typical terrains with the types: “water”, “forest”, “meadow” and “ground” were received. The patterns for the case of
Doppler filtration were also obtained
    In the second flight experiment the information for all tracks was processed by the designed algorithm, and we have
the time dependent on the correct detection probability (for example, see Fig. 8).
    The plots of the time dependence of the probability were analyzed in accordance to the video record. According to the
Neumann-Pearson criteria, the optimal level, which provided the maximum probability of the correct discrimination, was
evaluated, if the false detection was fixed. The information of Table 3 was obtained as the results of correct and false
detection.
    In the Table 3 next abbreviations were used: level – the optimal level of the correct probability detection without DF;
    level – the optimal level of the correct probability detection with DF;
    < > – the percentage of the objects which were detected correctly for the typical application of the radar;
    < > – the percentage of the objects which were detected correctly for the case of Doppler filtration.
                   Figure 8: The probability of correct detection for different terrains during the flight

               Table 3: The summary table for the results of the exploration during the flight experiments

                                  The terrain type      Water Forest Ground Meadow
                                Level Pcor.det.no.DF   0,4    0,6    0,6    0,6
                                Level Pcor.det.DF      0,4    0,5    0,6    0,5
                                < N%no.DF >            44     27     43     63
                                < N%DF >               60     31     31     54

    According to the Table 3, during the flight experiments it was found out that for typical application of the pulse radar
the detection of the “water” was provided in 44% cases, «forest» – 27% (because «forest» – significantly heterogeneous
terrain), «meadow» – 63%, and «ground» – 43%. In addition, the application of the DF makes it possible to increase the
probability detection for narrow LOs (the width is less than ¼ of the exposure spot’s width) of the “water” type at least
15%, but the probability of correct detection simultaneously and significantly decreases for “ground” and “meadow”,
because of wide backscattering diagram, that leads to decreasing the contribution of these terrains into the reflected
signal.
    So, it is recommended to apply both types of processing: the DF and typical processing. It depends from the type of
terrain which we want to detect.
    In the future, the series of flights above homogenous terrains are planned to be done. The experiments will be
implemented for creating more preciously pattern histograms

5    Conclusion
In this article the significant results were obtained:
1. Two cases of the created algorithm application were examined: the pulse radar altimeter with unsharpened beam and
the narrowed beam by the DF. It was shown that the regime of Doppler filtration allows us to increase the number of
detected terrain’s combinations at least twice. For example, if DF is used, the algorithm can detect next terrain
combinations (with narrow “water”) “Bushes/water/bushes”, “Grass/water/grass” and so on.
2. The flight experiments validated the correctness of the created algorithm. It is obtained, that the combination of
unsharpened beam and DF algorithms allow us to detect more than 60% of the narrow “water” objects, “meadow”,
“ground” and “forest” objects can also be found
3. The results of the current investigation allow us to choose the lengthened objects which are suitable for the UVA
navigation. It is shown that the terrain combinations “water/forest” and “asphalt/forest” provide the maximum of the
correct discrimination probability if LO is oriented to the flight direction with the angle from 30 up to 90 degrees and the
signal/noise ratio is not less than 0dB.
References
1.   Sorokin A. K., Vazhenin V. G., Shalgin V.V. Determination algorithm of underlying surfaces border position by
     pulse altimeter signal 2014 24nd Int. Crimean Conf. “Microwave & Telecommunication Technology”
     (CriMiCo’2014). Sevastopol, 2014, pp. 1163 -1164.
2.   Skolnik M. I.,”Radar Handbook”, Third Edition. The McGraw-Hill Companies, 2008. 1351p.
3.   Ulaby F. T. Handbook of Radar Scattering Statistics for Terrain. ARTECH HOUSE, 1989. 358p.
4.   Melnikov S.A. About the measurements of the terrain by the pulse radar altimeter application [Ob opredelenii
     izmerenij podstilajushhej poverhnosti s ispol'zovaniem vysotomernogo kanala] / S.A. Melnikov, R.V. Gordon, N.N.
     Kalmykov // Radar altimetry – 2013 : the proceedings of the IV All-Russian STC. – Ekaterinburg : FortDialog-Iset,
     2013. – pp. 82-91.