Signal Processing for Synthetic Aperture Sonar Image Enhancement

This thesis contains a description of SAS processing algorithms, offering improvements in Fourier-based reconstruction, motion-compensation, and autofocus. Fourier-based image reconstruction is reviewed and improvements shown as the result of improved system modelling. A number of new algorithms based on the wavenumber algorithm for correcting second order effects are proposed. In addition, a new framework for describing multiple-receiver reconstruction in terms of the bistatic geometry is presented and is a useful aid to understanding. Motion-compensation techniques for allowing Fourier-based reconstruction in widebeam geometries suffering large-motion errors are discussed. A motion-compensation algorithm exploiting multiple receiver geometries is suggested and shown to provide substantial improvement in image quality. New motion compensation techniques for yaw correction using the wavenumber algorithm are discussed. A common framework for describing phase estimation is presented and techniques from a number of fields are reviewed within this framework. In addition a new proof is provided outlining the relationship between eigenvector-based autofocus phase estimation kernels and the phase-closure techniques used astronomical imaging. Micronavigation techniques are reviewed and extensions to the shear average single-receiver micronavigation technique result in a 3 - 4 fold performance improvement when operating on high-contrast images. The stripmap phase gradient autofocus (SPGA) algorithm is developed and extends spotlight SAR PGA to the wide-beam, wide-band stripmap geometries common in SAS imaging. SPGA supersedes traditional PGA-based stripmap autofocus algorithms such as mPGA and PCA - the relationships between SPGA and these algorithms is discussed. SPGA's operation is verified on simulated and field-collected data where it provides significant image improvement. SPGA with phase-curvature based estimation is shown and found to perform poorly compared with phase-gradient techniques. The operation of SPGA on data collected from Sydney Harbour is shown with SPGA able to improve resolution to near the diffraction-limit. Additional analysis of practical stripmap autofocus operation in presence of undersampling and space-invariant blurring is presented with significant comment regarding the difficulties inherent in autofocusing field-collected data. Field-collected data from trials in Sydney Harbour is presented along with associated autofocus results from a number of algorithms

Sonar image interpretation for sub-sea operations

Mine Counter-Measure (MCM) missions are conducted to neutralise underwater explosives. Automatic Target Recognition (ATR) assists operators by increasing the speed and accuracy of data review. ATR embedded on vehicles enables adaptive missions which increase the speed of data acquisition. This thesis addresses three challenges; the speed of data processing, robustness of ATR to environmental conditions and the large quantities of data required to train an algorithm. The main contribution of this thesis is a novel ATR algorithm. The algorithm uses features derived from the projection of 3D boxes to produce a set of 2D templates. The template responses are independent of grazing angle, range and target orientation. Integer skewed integral images, are derived to accelerate the calculation of the template responses. The algorithm is compared to the Haar cascade algorithm. For a single model of sonar and cylindrical targets the algorithm reduces the Probability of False Alarm (PFA) by 80% at a Probability of Detection (PD) of 85%. The algorithm is trained on target data from another model of sonar. The PD is only 6% lower even though no representative target data was used for training. The second major contribution is an adaptive ATR algorithm that uses local sea-floor characteristics to address the problem of ATR robustness with respect to the local environment. A dual-tree wavelet decomposition of the sea-floor and an Markov Random Field (MRF) based graph-cut algorithm is used to segment the terrain. A Neural Network (NN) is then trained to filter ATR results based on the local sea-floor context. It is shown, for the Haar Cascade algorithm, that the PFA can be reduced by 70% at a PD of 85%. Speed of data processing is addressed using novel pre-processing techniques. The standard three class MRF, for sonar image segmentation, is formulated using graph-cuts. Consequently, a 1.2 million pixel image is segmented in 1.2 seconds. Additionally, local estimation of class models is introduced to remove range dependent segmentation quality. Finally, an A* graph search is developed to remove the surface return, a line of saturated pixels often detected as false alarms by ATR. The A* search identifies the surface return in 199 of 220 images tested with a runtime of 2.1 seconds. The algorithm is robust to the presence of ripples and rocks

Multichannel techniques for 3D ISAR

This thesis deals with the challenge of forming 3D target reconstruction by using spatial multi-channel ISAR configurations. The standard output of an ISAR imaging system is a 2D projection of the true three-dimensional target reflectivity onto an image plane. The orientation of the image plane cannot be predicted a priori as it strongly depends on the radar-target geometry and on the target motion, which is typically unknown. This leads to a difficult interpretation of the ISAR images. In this scenario, this thesis aim to give possible solutions to such problems by proposing three 3D processing based on interferometry, beamforming techniques and MIMO InISAR systems. The CLEAN method for scattering centres extraction is extended to multichannel ISAR systems and a multistatic 3D target reconstruction that is based on a incoherent technique is suggested

Multichannel techniques for 3D ISAR

This thesis deals with the challenge of forming 3D target reconstruction by using spatial multi-channel ISAR configurations. The standard output of an ISAR imaging system is a 2D projection of the true three-dimensional target reflectivity onto an image plane. The orientation of the image plane cannot be predicted a priori as it strongly depends on the radar-target geometry and on the target motion, which is typically unknown. This leads to a difficult interpretation of the ISAR images. In this scenario, this thesis aim to give possible solutions to such problems by proposing three 3D processing based on interferometry, beamforming techniques and MIMO InISAR systems. The CLEAN method for scattering centres extraction is extended to multichannel ISAR systems and a multistatic 3D target reconstruction that is based on a incoherent technique is suggested

Three Dimensional Bistatic Tomography Using HDTV

The thesis begins with a review of the principles of diffraction and reflection tomography; starting with the analytic solution to the inhomogeneous Helmholtz equation, after linearization by the Born approximation (the weak scatterer solution), and arriving at the Filtered Back Projection (Propagation) method of reconstruction. This is followed by a heuristic derivation more directly couched in the radar imaging context, without the rigor of the general inverse problem solution and more closely resembling an imaging turntable or inverse synthetic aperture radar. The heuristic derivation leads into the concept of the line integral and projections (the Radon Transform), followed by more general geometries where the plane wave approximation is invalid. We proceed next to study of the dependency of reconstruction on the space-frequency trajectory, combining the spatial aperture and waveform. Two and three dimensional apertures, monostatic and bistatic, fully and sparsely sampled and including partial apertures, with controlled waveforms (CW and pulsed, with and without modulation) define the filling of k-space and concomitant reconstruction performance. Theoretical developments in the first half of the thesis are applied to the specific example of bistatic tomographic imaging using High Definition Television (HDTV); the United States version of DVB-T. Modeling of the HDTV waveform using pseudonoise modulation to represent the hybrid 8VSB HDTV scheme and the move-stop-move approximation established the imaging potential, employing an idealized, isotropic 18 scatterer. As the move-stop-move approximation places a limitation on integration time (in cross correlation/pulse compression) due to transmitter/receiver motion, an exact solution for compensation of Doppler distortion is derived. The concept is tested with the assembly and flight test of a bistatic radar system employing software-defined radios (SDR). A three dimensional, bistatic collection aperture, exploiting an elevated commercial HDTV transmitter, is focused to demonstrate the principle. This work, to the best of our knowledge, represents a first in the formation of three dimensional images using bistatically-exploited television transmitters

Opportunistic radar imaging using a multichannel receiver

Bistatic Synthetic Aperture Radars have a physically separated transmitter and receiver where one or both are moving. Besides the advantages of reduced procurement and maintenance costs, the receiving system can sense passively while remaining covert which offers obvious tactical advantages. In this work, spaceborne monostatic SARs are used as emitters of opportunity with a stationary ground-based receiver. The imaging mode of SAR systems over land is usually a wide-swath mode such as ScanSAR or TOPSAR in which the antenna scans the area of interest in range to image a larger swath at the expense of degraded cross-range resolution compared to the conventional stripmap mode. In the bistatic geometry considered here, the signals from the sidelobes of the scanning beams illuminating the adjacent sub-swath are exploited to produce images with high cross-range resolution from data obtained from a SAR system operating in wide-swath mode. To achieve this, the SAR inverse problem is rigorously formulated and solved using a Maximum A Posteriori estimation method providing enhanced cross-range resolution compared to that obtained by classical burst-mode SAR processing. This dramatically increases the number of useful images that can be produced using emitters of opportunity. Signals from any radar satellite in the receiving band of the receiver can be used, thus further decreasing the revisit time of the area of interest. As a comparison, a compressive sensing-based method is critically analysed and proves more sensitive to off-grid targets and only suited to sparse scene. The novel SAR imaging method is demonstrated using simulated data and real measurements from C-band satellites such as RADARSAT-2 and ESA’s satellites ERS-2, ENVISAT and Sentinel-1A. In addition, this thesis analyses the main technological issues in bistatic SAR such as the azimuth-variant characteristic of bistatic data and the effect of imperfect synchronisation between the non-cooperative transmitter and the receiver