# Data post-processing algorithms for active forward-looking sonar system/Duomenu apdorojimo algoritmo taikymas aktyviosiose priesakiniu hidrolokatoriu sistemose.

IntroductionThis paper focuses on data processing for an active forward-looking sonar (FLS). Digital sonars, in general, are concerned that are transmitting different waveforms on the elements of a uniform linear array. The FLS enables row data processing that extracts target locations used for navigation safety purposes on large ships and vessels. The FLS images represent overlapping areas that can be combined to one image. Numerous different techniques have been developed for the object classification and interpretation of those images. However, many image processing algorithms consider images simply as a collection of points and do not associate them with physical sensors.

The FLS structure previously discussed [1] allows us to run beamforming and image post-processing algorithms automatically. Nevertheless, not only range data should be processed, the angular object information will be important as well. In particular, in the case of limited visibility the real time processing would improve the whole system performance. In practice, detailed information about the location and type of threat in shallow waters, which could not been determined by radar, is safety-critical. FLS applications would provide this information on a long distance. The redundant information in the raw data will be effectively used for data extraction necessary for navigation, image enhancement techniques and for the decrease of an output data stream. For a cost-effective solution wide beam transducers or multichannel reception modules with similar characteristics could be used as an alternative to narrow ones. In this case, an additional beamforming should be implemented. Moreover, the performance of the different beamforming algorithms should be modeled with the real sounding signals. This paper illustrates the efficiency of the proposed beamforming and post-processing algorithms in terms of angular aperture and target detection. Focus will be on the post-processing details and some of the algorithms implemented in FLS software applications.

Angular resolution in the case of spread spectrum signals

The angular resolution of a sonar system is determined by the width of the beam pattern - the narrower it is the better. In the case of electrical scanning (classical phased array) [2], the width of the beam pattern is determined by the spatial length of the sensor array, and the far-field beam pattern of the array when it is complex weighted is given by

D'(f,[f.sub.x]) = [N.summation over n=-N][a.sub.n](f)exp{+j[2[pi][f.sub.x]nd + [[theta].sub.n](f)]}, (1)

where [a.sub.n] (f)is an amplitude window of the array, N amount of sensors and d - distance between the sensors [4]. We know that a phase shift [[theta].sub.n] (f )in radians is equivalent to a time delay [[tau].sub.n] in seconds, that is

[[theta].sub.n](f)=2[pi]f[[tau].sub.n], [[tau].sub.n] = [[theta].sub.n](f) / 2[pi]f. (2)

In the case of long duration signals, the use of timedelay circuits to do beam steering is really not necessary, since beam steering can be accomplished via phase weighting, which can be implemented by using digital signal processing (digital beamforming). For signals with a spread spectrum, the directional diagrams of the sensor array may vary dynamically in time and we must directly compensate all the delays [[tau].sub.n]. It is possible to perform delay compensation and formation of the beam pattern in the frequency domain. In principle we must compute the complex frequency spectrum from the signal in each channel of the array, then multiply all given spectrum values in each channel with the expression exp[- j[[theta].sub.n](f)]. From the results we must take the inverse Fourier' transform and then obtain a delayed signal as a final result.

Several advanced digital algorithms may be used for delay compensation [2]. Instead of computing intensive FDFIB (Frequency Domain Frequency Invariant Beamformer) described before, an advanced Block-Phase (BP) method could be used [3, 4]. This method solves many problems related to beamforming and we can obtain good results in the output of the array. But the question is how to calculate and check the beam pattern in the case of spread spectrum signals. Assume that the output signal of the array in different partial directions [[beta].sub.y] is represented as [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. For a beam pattern we can calculate the output signal for each discrete direction as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)

[[beta].sub.y] = [d.sub.g]*[pi] / pa, (4)

where [p.sub.a] - amount of points in the direction diagram, [d.sub.g] - sampled partial directions dg = [bar.0,1,..., [p.sub.a]] [4].

A three-dimensional image of the two-dimensional signal matrix LX is shown in Fig. 1. Here [[beta].sub.y] = 30[degrees]. Horizontal axis represents the sampled partial directions dg and the vertical axis (0-200) represents the signal samples. We can firmly distinguish the presence of the spread spectrum signal at the (sampled) directions where the beam is steered. Moreover, we notice that the spread spectrum signal at the output of the array is quite stable in time and thus we can use some well-known signal amplitude estimation techniques to derive the beam pattern in polar coordinates.

[FIGURE 1 OMITTED]

[FIGURE 2 OMITTED]

The results of modeling showed that the signals with spread spectrum at large angles would widen the main lobe of the beam pattern approximately 4 times less than the classical phase compensation size. Fig.s 3 to 6 present modeling results of the horizontal far-field beam patterns in case of [[beta].sub.y] = 30[degrees] and [[beta].sub.y] = 10[degrees] (N=81) with conventional and Block-Phase methods.

[FIGURE 3 OMITTED]

[FIGURE 4 OMITTED]

[FIGURE 5 OMITTED]

[FIGURE 6 OMITTED]

Data post-processing algorithms

Synthetic aperture beamforming theory allows us to define location and the beamwidth of the main lobe of each channel [2]. However, the beamwidth depends on the number of physical sensors: the more sensors we use, the narrower the main beam. For example, the width of the synthetic aperture main_lobe [B.sub.0] of the sonar system with

equivalent distance antenna array with the number of array elements of N, distance between array elements d and with the sounding signals of wavelength X could be calculated as

[B.sub.0] [approximately equal to] 2 * [lambda] / d * 1 / N * (5)

Therefore, to reduce the angular aperture and increase the angular resolution the number of array elements should be increased. In practice, this will give N squared increase in the number of calculations for conventional beamforming algorithm. Concerning the FPGA implementation this will result in algorithm complexity that could affect the cost of the equipment.

The sonar applications with wide beam array elements would need additional post-processing. Nevertheless, the proposed additional method does not increase the angular resolution of the system but gives an increase in object location accuracy [5].

The navigation would need real-time information about potentially dangerous underwater objects in shallow water areas. As the basic information acquired forms the prior area maps, so the main thread sources are single floating and unmarked bottom objects. The goal of the discussed sonar application is to detect and locate those objects on a maximum distance. Moreover, object location accuracy could be achieved as the post-processing of the synthetic aperture FLS sonar data. In fact, this allows us to achieve desired antenna array by applying additional beamforming methods.

Furthermore, the sonar system could use the array antenna, which consists of separate narrow beam transducers with known characteristics. In this case, antenna with N elements separated from each other by d would have the beam pattern [F.sub.v]([beta]) and using a rectangular window could be formed as follows

[F.sub.v][beta] = sin([pi]Nd(sin[beta] - sin[[beta].sub.v])/[[lambda].sub.0] / N sin ([pi]d(sin[beta] - sin[[beta].sub.v])/[[lambda].sub.0] (6)

where [[lambda].sub.0] - wavelength of the sounding frequency, [[beta].sub.v] - partial direction [2][6]. However, other standard windows functions (Hamming, Kaiser, Chebyshev etc.) reduce sidelobe levels at the expense of widening the main lobe. Consequently, if the synthetics channel location gives an overlap of the main lobes we could reduce the angular aperture using the signal levels in neighboring channels. Here we consider the case of a single point object observed from the distance [7].

Generally, synthetic aperture is done with channel overlapping (Fig. 7.) and the main lobe of each channel is in [+ or -][pi]. So it could be represented as sin x/x [5].

In this case, the signal levels of an object located at a direction for 0 and 1 channels are:

[x.sub.0] = A * sin[alpha] / [alpha], (7)

[x.sub.1] = A * sin(-[pi] + [alpha]) / (-[pi] + [alpha]), (8)

where A - signal amplitude before the array processing. As a consequence

[alpha] = [pi] [x.sub.1] / [x.sub.1] + [x.sub.0]. (9)

[FIGURE 7 OMITTED]

What is more, this necessary object correction could be easily implemented in post-processing software application. The proposed post-processing algorithm was tested on the models and prototype application. The modeling includes real data from the 6-channel FLS sonar at 490 kHz frequency. Spread spectrum signals are with 32 periods per signal element. Fig.s 8-9 show modeled raw FLS data and post-processed object data. The modeling showed that the object seen in multiple neighboring channels could be successfully localized.

[FIGURE 8 OMITTED]

[FIGURE 9 OMITTED]

Proposed post-processing is implemented in a software defined forward-looking sonar prototype [1]. So the signal processing includes sounding signal generation, data reception from sensors, beamforming, optimal filtering, data post-processing, and image representation (Fig. 10, Fig. 11).

[FIGURE 10 OMITTED]

[FIGURE 11 OMITTED]

Conclusions

Active digital forward-looking sonar could be implemented as a software defined application. The redundant information in the raw data could be effectively used for data extraction necessary for navigation. Additional post-processing in cost-effective solutions could give an increase in object location accuracy. This method has some limitations. Processing of the signal levels in neighboring channels would increase noise sensitivity, especially for distant objects. In practice, array calibration coefficients should be applied at the beamforming stage. Furthermore, additional research of the target classification is required. Future work will extend the processing and data fusion capabilities [7, 8]. The data received should be combined with vessel motion, speed and depth at the measured point.

References

[1.] Berdnikova J., Kozevnikov V., Zamarajev J., Raja A. Advanced Forward-Looking Sonar and Imagery Data Processing // Electronics and Electrical Engineering. Kaunas: Technologija, 2010.--No. 4(100).--P. 55- 58.

[2.] Bhattacharyya A. K. Phased Array Antennas.--WileyInterscience, 2006.--496 p.

[3.] Ruuben T, Derkats J. Some Methods of Signal Processing and Beamforming in Hydrographic Applications // Electronics and Electrical Engineering.--Kaunas: Technologija, 2008.--No. 6(86).--P. 27--32.

[4.] Ruuben T. Wideband Digital Beamforming in Sonar Systems. Doctoral thesis.--Tallinn University of Technology Press, 2009. -140 p.

[5.] Ricker D. W. Echo Signal Processing.--Kluwer Academic Publishers, 2003.--484 p

[6.] Nielsen R. O. Sonar Signal Processing.--Artech House Inc., 1991.--368 p.

[7.] Blondel F., Murton B. J. Handbook of seafloor sonar imagery.--John Wiley & Sons Ltd, 1997.--314 p.

[8.] Bloch I. Information Fusion in Signal and Image Processing.--John Wiley & Sons Inc., 2008.--295 p.

Received 2011 03 13

J. Berdnikova, V. Kozevnikov, T. Ruuben, A. Raja

Department of Radio and Communication Engineering, Tallinn University of Technology, Ehitajate tee 5, 19086 Tallinn, Estonia, phone: +3726202350, e-mail: juliad@lr.ttu.ee

Printer friendly Cite/link Email Feedback | |

Title Annotation: | TELECOMMUNICATIONS ENGINEERING / TELEKOMUNIKACIJU INZINERIJA |
---|---|

Author: | Berdnikova, J.; Kozevnikov, V.; Ruuben, T.; Raja, A. |

Publication: | Elektronika ir Elektrotechnika |

Article Type: | Report |

Geographic Code: | 4EXES |

Date: | Apr 1, 2011 |

Words: | 1896 |

Previous Article: | Optimization of sensor grid aided by modelling of sensor relationships with digital filters/Jutiklio tinklelio taikymo optimizavimas modeliuojant... |

Next Article: | Synthesis of six-conductors symmetrically coupled microstrip line, operating in a normal mode/Simetriskai susietu sesiu laidininku mikrojuostelines... |

Topics: |