Development of Time - Frequency Techniques for Sonar Applications

Sonar signal processing comprises of a large number of signal processing algorithms for implementing functions such as Target Detection, Localisation, Classification, Tracking and Parameter estimation. Current implementations of these functions rely on conventional techniques largely based on Fourier Techniques, primarily meant for stationary signals. Interestingly enough, the signals received by the sonar sensors are often non-stationary and hence processing methods capable of handling the non-stationarity will definitely fare better than Fourier transform based methods.Time-frequency methods(TFMs) are known as one of the best DSP tools for nonstationary signal processing, with which one can analyze signals in time and frequency domains simultaneously. But, other than STFT, TFMs have been largely limited to academic research because of the complexity of the algorithms and the limitations of computing power. With the availability of fast processors, many applications of TFMs have been reported in the fields of speech and image processing and biomedical applications, but not many in sonar processing. A structured effort, to fill these lacunae by exploring the potential of TFMs in sonar applications, is the net outcome of this thesis. To this end, four TFMs have been explored in detail viz. Wavelet Transform, Fractional Fourier Transfonn, Wigner Ville Distribution and Ambiguity Function and their potential in implementing five major sonar functions has been demonstrated with very promising results. What has been conclusively brought out in this thesis, is that there is no "one best TFM" for all applications, but there is "one best TFM" for each application. Accordingly, the TFM has to be adapted and tailored in many ways in order to develop specific algorithms for each of the applications.

Probabilistic sonar scan matching for an AUV

This paper proposes MSISpIC, a probabilistic sonar scan matching algorithm for the localization of an autonomous underwater vehicle (AUV). The technique uses range scans gathered with a Mechanical Scanning Imaging Sonar (MSIS), the robot displacement estimated through dead-reckoning using a Doppler velocity log (DVL) and a motion reference unit (MRU). The proposed method is an extension of the pIC algorithm. An extended Kalman filter (EKF) is used to estimate the robot-path during the scan in order to reference all the range and bearing measurements as well as their uncertainty to a scan fixed frame before registering. The major contribution consists of experimentally proving that probabilistic sonar scan matching techniques have the potential to improve the DVL-based navigation. The algorithm has been tested on an AUV guided along a 600 m path within an abandoned marina underwater environment with satisfactory results

Dam wall detection and tracking using a Mechanically Scanned Imaging Sonar

In dam inspection tasks, an underwater robot has to grab images while surveying the wall meanwhile maintaining a certain distance and relative orientation. This paper proposes the use of an MSIS (mechanically scanned imaging sonar) for relative positioning of a robot with respect to the wall. An imaging sonar gathers polar image scans from which depth images (range & bearing) are generated. Depth scans are first processed to extract a line corresponding to the wall (with the Hough transform), which is then tracked by means of an EKF (Extended Kalman Filter) using a static motion model and an implicit measurement equation associating the sensed points to the candidate line. The line estimate is referenced to the robot fixed frame and represented in polar coordinates (rho&thetas) which directly corresponds to the actual distance and relative orientation of the robot with respect to the wall. The proposed system has been tested in simulation as well as in water tank conditions

Pose-based SLAM with probabilistic scan matching algorithm using a mechanical scanned imaging sonar

This paper proposes a pose-based algorithm to solve the full SLAM problem for an autonomous underwater vehicle (AUV), navigating in an unknown and possibly unstructured environment. The technique incorporate probabilistic scan matching with range scans gathered from a mechanical scanning imaging sonar (MSIS) and the robot dead-reckoning displacements estimated from a Doppler velocity log (DVL) and a motion reference unit (MRU). The proposed method utilizes two extended Kalman filters (EKF). The first, estimates the local path travelled by the robot while grabbing the scan as well as its uncertainty and provides position estimates for correcting the distortions that the vehicle motion produces in the acoustic images. The second is an augment state EKF that estimates and keeps the registered scans poses. The raw data from the sensors are processed and fused in-line. No priory structural information or initial pose are considered. The algorithm has been tested on an AUV guided along a 600 m path within a marina environment, showing the viability of the proposed approach

SLAM using an Imaging Sonar for Partially Structured Underwater Environments

In this paper we describe a system for underwater navigation with AUVs in partially structured environments, such as dams, ports or marine platforms. An imaging sonar is used to obtain information about the location of planar structures present in such environments. This information is incorporated into a feature-based SLAM algorithm in a two step process: (I) the full 360deg sonar scan is undistorted (to compensate for vehicle motion), thresholded and segmented to determine which measurements correspond to planar environment features and which should be ignored; and (2) SLAM proceeds once the data association is obtained: both the vehicle motion and the measurements whose correct association has been previously determined are incorporated in the SLAM algorithm. This two step delayed SLAM process allows to robustly determine the feature and vehicle locations in the presence of large amounts of spurious or unrelated measurements that might correspond to boats, rocks, etc. Preliminary experiments show the viability of the proposed approach

Integració d'un sensor sonar d'escombrat lateral en un robot submarí

El grup de Visió per Computador i Robòtica (VICOROB) disposa de varis robots submarins per a la recerca i inspecció subaquàtica. Recentment s’ha adquirit un sensor sonar d’escombrat lateral el qual s’utilitza per realitzar imatges acústiques del fons marí quan aquest es mou principalment a velocitat constant i mantenint el rumb. Els robots del grup VICOROB estan equipats amb diferents tipus de sensors i càmeres per analitzar el fons marí. Aquest sensors són de gran qualitat i permeten conèixer de manera bastant satisfactòria l’entorn a les proximitats del robot. Freqüentment però, aquest sensors estant sotmesos a diferents restriccions depenent de la seva naturalesa de funcionament, de tal manera que es necessària la seva combinació per resoldre determinats problemes en diferents situacions. Amb aquest projecte, es pretén integrar un nou sistema de captura d’imatges sonores del fons marí, en un dels robots. Amb la integració d’aquest nou sensor, s’espera obtenir una opció alternativa els sistemes actuals que pugui aportar informació addicional sobre el fons. Aquest sistema podrà ser utilitzat per realitzar tasques per les quals els altres sensors no estant preparats o bé per complementar informació d’altres sensor

Underwater slam for estructured environments using and imaging sonar

Aquesta tesi tracta sobre el problema de la navegació per a vehicles submarins autònoms que operen en entorns artificials estructurats com ara ports, canals, plataformes marines i altres escenaris similars. A partir d'una estimació precisa de la posició en aquests entorns, les capacitats dels vehicles submarins s'incrementen notablement i s'obre una porta al seu funcionament autònom. El manteniment, inspecció i vigilància d'instal lacions marines són alguns exemples de possibles aplicacions. Les principals contribucions d'aquesta tesi consisteixen per una banda en el desenvolupament de diferents sistemes de localització per a aquelles situacions on es disposa d'un mapa previ de l'entorn i per l'altra en el desenvolupament d'una nova solució al problema de la Localització i Construcció Simultània de Mapes (SLAM en les seves sigles en anglès), la finalitat del qual és fer que un vehicle autònom creï un mapa de l'entorn desconegut que el rodeja i, al mateix temps, utilitzi aquest mapa per a determinar la seva pròpia posició. S'ha escollit un sonar d'imatges d'escaneig mecànic com a sensor principal per a aquest treball tant pel seu relatiu baix cost com per la seva capacitat per produir una representació detallada de l'entorn. Per altra banda, les particularitats de la seva operació i, especialment, la baixa freqúència a la que es produeixen les mesures, constitueixen els principals inconvenients que s'han hagut d'abordar en les estratègies de localització proposades. Les solucions adoptades per aquests problemes constitueixen una altra contribució d'aquesta tesi. El desenvolupament de vehicles autònoms i el seu ús com a plataformes experimentals és un altre aspecte important d'aquest treball. Experiments portats a terme tant en el laboratori com en escenaris reals d'aplicació han proporcionat les dades necessàries per a provar i avaluar els diferents sistemes de localització proposats.

Habitat suitability for marine fishes using presence-only modelling and multibeam sonar

Improved access to multibeam sonar and underwater video technology is enabling scientists to use spatially-explicit, predictive modelling to improve our understanding of marine ecosystems. With the growing number of modelling approaches available, knowledge of the relative performance of different models in the marine environment is required. Habitat suitability of 5 demersal fish taxa in Discovery Bay, south-east Australia, were modelled using 10 presence-only algorithms: BIOCLIM, DOMAIN, ENFA (distance geometric mean [GM], distance harmonic mean [HM], median [M], area-adjusted median [Ma], median + extremum [Me], area-adjusted median + extremum [Mae] and minimum distance [Min]), and MAXENT. Model performance was assessed using kappa and area under curve (AUC) of the receiver operator characteristic. The influence of spatial range (area of occupancy) and environmental niches (marginality and tolerance) on modelling performance were also tested. MAXENT generally performed best, followed by ENFA-GM and -HM, DOMAIN, BIOCLIM, ENFA-M, -Min, -Ma, -Mae and -Me algorithms. Fish with clearly definable niches (i.e. high marginality) were most accurately modelled. Generally, Euclidean distance to nearest reef, HSI-b (backscatter), rugosity and maximum curvature were the most important variables in determining suitable habitat for the 5 demersal fish taxa investigated. This comparative study encourages ongoing use of presence-only approaches, particularly MAXENT, in modelling suitable habitat for demersal marine fishes.

Understanding fish-habitat associations using video observations and sonar imaging

This thesis was motivated by the increasing role of species distribution models in managing marine fishes and their habitats. These models provided new information about fish-habitat associations. However, models are potentially influenced by fish behaviour towards underwater video systems.

Submarine sonar FX

Submarine Sonar FX is a sound effect of a submarine/nautical type of sonar.

Evaluation of four supervised learning methods for benthic habitat mapping using backscatter from multi-beam sonar

An understanding of the distribution and extent of marine habitats is essential for the implementation of ecosystem-based management strategies. Historically this had been difficult in marine environments until the advancement of acoustic sensors. This study demonstrates the applicability of supervised learning techniques for benthic habitat characterization using angular backscatter response data. With the advancement of multibeam echo-sounder (MBES) technology, full coverage datasets of physical structure over vast regions of the seafloor are now achievable. Supervised learning methods typically applied to terrestrial remote sensing provide a cost-effective approach for habitat characterization in marine systems. However the comparison of the relative performance of different classifiers using acoustic data is limited. Characterization of acoustic backscatter data from MBES using four different supervised learning methods to generate benthic habitat maps is presented. Maximum Likelihood Classifier (MLC), Quick, Unbiased, Efficient Statistical Tree (QUEST), Random Forest (RF) and Support Vector Machine (SVM) were evaluated to classify angular backscatter response into habitat classes using training data acquired from underwater video observations. Results for biota classifications indicated that SVM and RF produced the highest accuracies, followed by QUEST and MLC, respectively. The most important backscatter data were from the moderate incidence angles between 30° and 50°. This study presents initial results for understanding how acoustic backscatter from MBES can be optimized for the characterization of marine benthic biological habitats. © 2012 by the authors.

Comparison of a self-processed EM3000 multibeam echosounder dataset with a QTC view habitat mapping and a sidescan sonar imagery, Tamaki Strait, New Zealand

A methodology for automatically processing the data files from an EM3000 multibeam echosounder (Kongsberg Maritime, 300 kHz) is presented. Written in MatLab, it includes data extraction, bathymetry processing, computation of seafloor local slope, and a simple correction of the backscatter along-track banding effect. The success of the latter is dependent on operational restrictions, which are also detailed. This processing is applied to a dataset acquired in 2007 in the Tamaki Strait, New Zealand. The resulting maps are compared with a habitat classification obtained with the acoustic ground-discrimination software QTC View linked to a 200-kHz single-beam echosounder and to the imagery from a 100-kHz sidescan sonar survey, both performed in 2002. The multibeam backscatter map was found to be very similar to the sidescan imagery, quite correlated to the QTC View map on one site but mainly uncorrelated on another site. Hypotheses to explain these results are formulated and discussed. The maps and the comparison to prior surveys are used to draw conclusions on the quality of the code for further research on multibeam benthic habitat mapping.

Granitic coastal geomorphology: Applying integrated terrestrial and bathymetric LiDAR with multibeam sonar to examine coastal landscape evolution

Coasts composed of resistant lithologies such as granite are generally highly resistant to erosion. They tend to evolve over multiple sea level cycles with highstands acting to remove subaerially weathered material. This often results in a landscape dominated by plunging cliffs with shore platforms rarely occurring. The long-term evolution of these landforms means that throughout the Quaternary these coasts have been variably exposed to different sea level elevations which means erosion may have been concentrated at different elevations from today. Investigations of the submarine landscape of granitic coasts have however been hindered by an inability to accurately image the nearshore morphology. Only with the advent of multibeam sonar and aerial laser surveying can topographic data now be seamlessly collected from above and below sea level. This study tests the utility of these techniques and finds that very accurate measurements can be made of the nearshore thereby allowing researchers to study the submarine profile with the same accuracy as the subaerial profile. From a combination of terrestrial and marine LiDAR data with multibeam sonar data, it is found that the morphology of granite domes is virtually unaffected by erosion at sea level. It appears that evolution of these landscapes on the coast is a very slow process with modern sea level acting only to remove subaerially weathered debris. The size and orientation of the joints determines the erosional potential of the granite. Where joints are densely spaced (<2 m apart) or the bedrock is highly weathered can semi-horizontal surfaces form.

Accounting for uncertainty in volumes of seabed change measured with repeat multibeam sonar surveys

Seafloors of unconsolidated sediment are highly dynamic features; eroding or accumulating under the action of tides, waves and currents. Assessing which areas of the seafloor experienced change and measuring the corresponding volumes involved provide insights into these important active sedimentation processes. Computing the difference between Digital Elevation Models (DEMs) obtained from repeat Multibeam Echosounders (MBES) surveys has become a common technique to identify these areas, but the uncertainty in these datasets considerably affects the estimation of the volumes displaced. The two main techniques used to take into account uncertainty in volume estimations are the limitation of calculations to areas experiencing a change in depth beyond a chosen threshold, and the computation of volumetric confidence intervals. However, these techniques are still in their infancy and, as a result, are often crude, seldom used or poorly understood. In this article, we explored a number of possible methodological advances to address this issue, including: (1) using the uncertainty information provided by the MBES data processing algorithm CUBE, (2) adapting fluvial geomorphology techniques for volume calculations using spatially variable thresholds and (3) volumetric histograms. The nearshore seabed off Warrnambool harbour - located in the highly energetic southwest Victorian coast, Australia - was used as a test site. Four consecutive MBES surveys were carried out over a four-months period. The difference between consecutive DEMs revealed an area near the beach experiencing large sediment transfers - mostly erosion - and an area of reef experiencing increasing deposition from the advance of a nearby sediment sheet. The volumes of sediment displaced in these two areas were calculated using the techniques described above, both traditionally and using the suggested improvements. We compared the results and discussed the applicability of the new methodological improvements. We found that the spatially variable uncertainty derived from the CUBE algorithm provided the best results (i.e. smaller confidence intervals), but that similar results can be obtained using as a fixed uncertainty value derived from a reference area under a number of operational conditions.