Factorized Geometrical Autofocus for Synthetic Aperture Radar Processing

This paper describes a factorized geometrical autofocus (FGA) algorithm, specifically suitable for ultrawideband synthetic aperture radar. The strategy is integrated in a fast factorized back-projection chain and relies on varying track parameters step by step to obtain a sharp image; focus measures are provided by an object function (intensity correlation). The FGA algorithm has been successfully applied on synthetic and real (Coherent All RAdio BAnd System II) data sets, i.e., with false track parameters introduced prior to processing, to set up constrained problems involving one geometrical quantity. Resolution (3 dB in azimuth and slant range) and peak-to-sidelobe ratio measurements in FGA images are comparable with reference results (within a few percent and tenths of a decibel), demonstrating the capacity to compensate for residual space variant range cell migration. The FGA algorithm is finally also benchmarked (visually) against the phase gradient algorithm to emphasize the advantage of a geometrical autofocus approach

Autofocus and analysis of geometrical errors within the framework of fast factorized back-projection

This paper describes a Fast Factorized Back-Projection (FFBP) formulation that includes a fully integrated autofocus algorithm, i.e. the Factorized Geometrical Autofocus (FGA) algorithm. The base-two factorization is executed in a horizontal plane, using a Merging (M) and a Range History Preserving (RHP) transform. Six parameters are adopted for each sub-aperture pair, i.e. to establish the geometry stage-by-stage via triangles in 3-dimensional space. If the parameters are derived from navigation data, the algorithm is used as a conventional processing chain. If the parameters on the other hand are varied from a certain factorization step and forward, the algorithm is used as a joint image formation and autofocus strategy. By regulating the geometry at multiple resolution levels, challenging defocusing effects, e.g. residual space-variant Range Cell Migration (RCM), can be corrected. The new formulation also serves another important purpose, i.e. as a parameter characterization scheme. By using the FGA algorithm and its inverse, relations between two arbitrary geometries can be studied, in consequence, this makes it feasible to analyze how errors in navigation data, and topography, affect image focus. The versatility of the factorization procedure is demonstrated successfully on simulated Synthetic Aperture Radar (SAR) data. This is achieved by introducing different GPS/IMU errors and Focus Target Plane (FTP) deviations prior to processing. The characterization scheme is then employed to evaluate the sensitivity, to determine at what step the autofocus function should be activated, and to decide the number of necessary parameters at each step. Resulting FGA images are also compared to a reference image (processed without errors and autofocus) and to a defocused image (processed without autofocus), i.e. to validate the novel approach further

Autofocus and Back-Projection in Synthetic Aperture Radar Imaging.

Spotlight-mode Synthetic Aperture Radar (SAR) imaging has received considerable attention due to its ability to produce high-resolution images of scene reflectivity. One of the main challenges in successful image recovery is the problem of defocusing, which occurs due to inaccuracies in the estimated round-trip delays of the transmitted radar pulses. The problem is most widely studied for far-field imaging scenarios with a small range of look angles since the problem formulation can be significantly simplified under the assumptions of planar wavefronts and one-dimensional defocusing. In practice, however, these assumptions are frequently violated. MultiChannel Autofocus (MCA) is a subspace-based approach to the defocusing problem that was originally proposed for far-field imaging, with a small range of look angles. A key motivation behind MCA is the observation that there exists a low-return region within the recovered image, due to the weak illumination near the edges of the antenna footprint. The strength of the MCA formulation is that it can be easily extended to more realistic scenarios with polar-format data, spherical wavefronts, and arbitrary terrain, due to its flexible linear-algebraic framework. The main aim of this thesis is to devise a more broadly effective autofocus approach by adopting MCA to the aforementioned scenarios. By forming the solution space in a domain where the defocusing effect is truly one-dimensional, we show that drastically improved restorations can be obtained for applications with small to fairly wide ranges of look angles. When the terrain topography is known, we utilize the versatile backprojection-based imaging methods in the model formulations for MCA to accurately account for the underlying geometry. The proposed extended MCA shows reductions in RMSE of up to 50% when the underlying terrain is highly elevated. We also analyze the effects of the filtering step, the amount of wave curvature, the shape of the terrain, and the flight path of the radar, on the reconstructed image via backprojection. Finally, we discuss the selection of low-return constraints and the importance of using terrain elevation within MCA formulation.PHDElectrical Engineering: SystemsUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/135868/1/zzon_1.pd

An Efficient Solution to the Factorized Geometrical Autofocus Problem

This paper describes a new search strategy within the scope of factorized geometrical autofocus (FGA) and synthetic-aperture-radar processing. The FGA algorithm is a fast factorized back-projection formulation with six adjustable geometry parameters. By tuning the flight track step by step and maximizing focus quality by means of an object function, a sharp image is formed. We propose an efficient two-stage approach for the geometrical variation. The first stage is a low-order (few parameters) parallel search procedure involving small image areas. The second stage then combines the local hypotheses into one global autofocus solution, without the use of images. This method has been applied successfully on ultrawideband CARABAS II data. Errors due to a constant acceleration are superposed on the measured track prior to processing, giving a 6-D autofocus problem. Image results, including resolution, peak-to-sidelobe ratio and magnitude values for point-like targets, finally confirm the validity of the strategy. The results also verify the prediction that there are several satisfying autofocus solutions for the same radar data

A Generalized Phase Gradient Autofocus Algorithm

The phase gradient autofocus (PGA) algorithm has seen widespread use and success within the synthetic aperture radar (SAR) imaging community. However, its use and success has largely been limited to collection geometries where either the polar format algorithm (PFA) or range migration algorithm is suitable for SAR image formation. In this work, a generalized phase gradient autofocus (GPGA) algorithm is developed which is applicable with both the PFA and backprojection algorithm (BPA), thereby directly supporting a wide range of collection geometries and SAR imaging modalities. The GPGA algorithm preserves the four crucial signal processing steps comprising the PGA algorithm, while alleviating the constraint of using a single scatterer per range cut for phase error estimation which exists with the PGA algorithm. Moreover, the GPGA algorithm, whether using the PFA or BPA, yields an approximate maxi- mum marginal likelihood estimate (MMLE) of phase errors having marginalized over unknown complex-valued reflectivities of selected scatterers. Also, in this work a new approximate MMLE, termed the max-semidefinite relaxation (Max-SDR) phase estimator, is proposed for use with the GPGA algorithm. The Max-SDR phase estimator provides a phase error estimate with a worst-case approximation bound compared to the solution set of MMLEs (i.e., solution set to the non-deterministic polynomial- time hard (NP-hard) GPGA phase estimation problem). Moreover, in this work a specialized interior-point method is presented for more efficiently performing Max- SDR phase estimation by exploiting low-rank structure typically associated with the GPGA phase estimation problem. Lastly, simulation and experimental results produced by applying the GPGA algorithm with the PFA and BPA are presented

Computational Algorithms for Improved Synthetic Aperture Radar Image Focusing

High-resolution radar imaging is an area undergoing rapid technological and scientific development. Synthetic Aperture Radar (SAR) and Inverse Synthetic Aperture Radar (ISAR) are imaging radars with an ever-increasing number of applications for both civilian and military users. The advancements in phased array radar and digital computing technologies move the trend of this technology towards higher spatial resolution and more advanced imaging modalities. Signal processing algorithm development plays a key role in making full use of these technological developments.In SAR and ISAR imaging, the image reconstruction process is based on using the relative motion between the radar and the scene. An important part of the signal processing chain is the estimation and compensation of this relative motion. The increased spatial resolution and number of receive channels cause the approximations used to derive conventional algorithms for image reconstruction and motion compensation to break down. This leads to limited applicability and performance limitations in non-ideal operating conditions.This thesis presents novel research in the areas of data-driven motion compensation and image reconstruction in non-cooperative ISAR and Multichannel Synthetic Aperture Radar (MSAR) imaging. To overcome the limitations of conventional algorithms, this thesis proposes novel algorithms leading to increased estimation performance and image quality. Because a real-time imaging capability is important in many applications, special emphasis is placed on the computational aspects of the algorithms.For non-cooperative ISAR imaging, the thesis proposes improvements to the range alignment, time window selection, autofocus, time-frequency-based image reconstruction and cross-range scaling procedures. These algorithms are combined into a computationally efficient non-cooperative ISAR imaging algorithm based on mathematical optimization. The improvements are experimentally validated to reduce the computational burden and significantly increase the image quality under complex target motion dynamics.Time domain algorithms offer a non-approximated and general way for image reconstruction in both ISAR and MSAR. Previously, their use has been limited by the available computing power. In this thesis, a contrast optimization approach for time domain ISAR imaging is proposed. The algorithm is demonstrated to produce improved imaging performance under the most challenging motion compensation scenarios. The thesis also presents fast time domain algorithms for MSAR. Numerical simulations confirm that the proposed algorithms offer a reasonable compromise between computational speed and image quality metrics

Advances in Synthetic Aperture Radar from a Wavenumber Perspective

This dissertation examines the wavenumber domain of Synthetic Aperture Radar (SAR) images. This domain is the inverse Fourier transform domain of a SAR image. The dissertation begins with the radar receiver's signal model and develops equations describing the wavenumber domain of a SAR image produced by a generalized bistatic and monostatic SAR system. Then, closed form expressions for bistatic synthetic aperture radar spatial resolution of a generalized system from the wavenumber domain are developed. These spatial resolution equations have not previously appeared in the literature. From these equations, significant resolution is found in both range and cross-range forecasting a forward-scatter bistatic SAR image when the elevation angles of each bistatic platform are significantly different. Next, wavenumber and time domain image formation algorithms are discussed. Developed within this dissertation is a wavenumber preprocessing method that increases the speed of the Back Projection Algorithm (BPA). This preprocessing method takes advantage of deramped SAR radar returns and their polar wavenumber format. This new algorithm is called the Fast Decimated Wavenumber Back Projection Algorithm (FDWBPA). Matlab functions are included to implement this algorithm, simulate bistatic SAR images and process the data from anechoic chamber tests demonstrating forward scatter resolution

Iterative synthetic aperture radar imaging algorithms

Synthetic aperture radar is an important tool in a wide range of civilian and military imaging applications. This is primarily due to its ability to image in all weather conditions, during both the day and the night, unlike optical imaging systems. A synthetic aperture radar system contains a step which is not present in an optical imaging system, this is image formation. This is required because the acquired data from the radar sensor does not directly correspond to the image. Instead, to form an image, the system must solve an inverse problem. In conventional scenarios, this inverse problem is relatively straight forward and a matched lter based algorithm produces an image of suitable image quality. However, there are a number of interesting scenarios where this is not the case. Scenarios where standard image formation algorithms are unsuitable include systems with data undersampling, errors in the system observation model and data that is corrupted by radio frequency interference. Image formation in these scenarios will form the topics of this thesis and a number of iterative algorithms are proposed to achieve image formation. The motivation for these proposed algorithms is primarily from the eld of compressed sensing, which considers the recovery of signals with a low-dimensional structure. The rst contribution of this thesis is the development of fast algorithms for the system observation model and its adjoint. These algorithms are required by large-scale gradient based iterative algorithms for image formation. The proposed algorithms are based on existing fast back-projection algorithms, however, a new decimation strategy is proposed which is more suitable for some applications. The second contribution is the development of a framework for iterative near- eld image formation, which uses the proposed fast algorithms. It is shown that the framework can be used, in some scenarios, to improve the visual quality of images formed from fully sampled data and undersampled data, when compared to images formed using matched lter based algorithms. The third contribution concerns errors in the system observation model. Algorithms that correct these errors are commonly referred to as autofocus algorithms. It is shown that conventional autofocus algorithms, which work as a post-processor on the formed image, are unsuitable for undersampled data. Instead an autofocus algorithm is proposed which corrects errors within the iterative image formation procedure. The proposed algorithm is provably stable and convergent with a faster convergence rate than previous approaches. The nal contribution is an algorithm for ultra-wideband synthetic aperture radar image formation. Due to the large spectrum over which the ultra-wideband signal is transmitted, there is likely to be many other users operating within the same spectrum. These users can produce signi cant radio frequency interference which will corrupt the received data. The proposed algorithm uses knowledge of the RFI spectrum to minimise the e ect of the RFI on the formed image

Radar Imaging in Challenging Scenarios from Smart and Flexible Platforms

undefine

Contribution to ground-based and UAV SAR systems for Earth observation

Mankind's way of life is the main driver of a planetary-scale change that is marked by the growing of human population's demand of energy, food, goods, services and information. As a result, it have emerged new ecological, economical, social and geopolitical concerns. In this scenario, SAR remote sensing is a potential tool that provides unique information about the Earth's properties and processes that can be used to solve societal challenges of local and global dimension. SARs, which are coherent systems that are able to provide high resolution images with weather independence, represent a suitable alternative for EO with diverse applications. Some examples of SAR application areas are topography (DEM generation with interferometry), agriculture (crop classification or soil moisture), or geology (monitoring surface deformation). In this framework, the encompassing objective of the present doctoral work has been part of the implementation and the subsequent evaluation of capabilities of two X-band SAR sensors. On the one hand, the RISKSAR-X radar designed to be operated at ground to monitor small-scale areas of observation and, on the other, the ARBRES-X sensor designed to be integrated into small UAVs. Despite its inherently dissimilar conception, the concurrence of both sensors has been evidenced along this manuscript. By taking advantage of the similarities between them, it has been possible to analogously assess both sensors to obtain conclusions. In this context, the common link has been the development of the polarimetric OtF operation mode of the RISKSAR-X, allowing this sensor to be operated equivalently to the ARBRES-X. Regarding the RISKSAR-X SAR sensor, several hardware contributions have been developed during part of this Ph.D. with the aim of improving the system performance. By endowing the system with the capability to operate in the fully polarimetric OtF acquisition mode, the relative long scanning time has been reduced. It is of great interest since the measured scatterers that present a short term variable reflectivity during the scanning time, such as moving vegetation, may degrade the extracted parameters from the retrieved data and the SAR image reconstruction. During this doctoral activity, it has been studied the image blurring, the decorrelation and the coherence degradation introduced by this effect. Furthermore, a new term in the differential interferometric coherence that takes into account the image blurring has been introduced. Concerning the ARBRES-X SAR system, one of the main objectives pursued during this Ph.D. has been the integration of the sensor into a small UAV MP overcoming restrictions of weight, size and aerodynamics of the platform. The use of this type of platforms is expected to open up new possibilities in airborne SAR remote sensing, since it offers much more versatility than the commonly used fixed wings UAVs. Different innovative flight strategies with this type of platforms have been assessed and some preliminary results have been obtained with the use of the ARBRES-X SAR system. During the course of the present doctoral work, much effort has been devoted to achieve the first experimental repeat-pass interfereometric results obtained with the UAV MP together with the ARBRES-X. Moreover, the sensor has been endowed with fully polarimetric capabilities by applying the improvements developed to the RISKSAR-X radar, which is another example of the duality between both systems. Finally, a vertical and a semicircular aperture have been successfully performed obtaining SLC images of the scenario, which envisages the capability of the UAV MP to perform tomographic images and complete circular apertures in the future. In conclusion, the UAV MP is a promising platform that opens new potentials for several applications, such as repeat-pass interferometry or differential tomography imaging with the realization of almost arbitrary trajectories.El mode de viure de la humanitat és el principal motor d'un canvi a escala planetària que està marcat per la creixent demanda d'energia, d'aliment, de béns, de serveis i d'informació de les poblacions humanes. Com a resultat, han sorgit noves inquietuds ecològiques, econòmiques, socials i geopolítiques. En aquest escenari, la detecció remota SAR és una eina potencial que proporciona informació única sobre les propietats i processos de la Terra que es pot utilitzar per resoldre reptes socials de dimensió local i global. Els SARs, que són sistemes coherents que poden proporcionar imatges d'alta resolució amb independència del temps, representen una alternativa adequada per a l'observació de la Terra. Alguns exemples d'àrees d'aplicació SAR són la topografia (generació de DEM amb interferometria), l'agricultura (classificació de cultius o humitat del sòl) o la geologia (monitoratge de deformació superficial). En aquest context, l'objectiu general del present doctorat ha estat part de la implementació i posterior avaluació de les capacitats de dos sensors SAR de banda X. D'una banda, el radar RISKSAR-X dissenyat per funcionar a terra i monitoritzar àrees d'observació a petita escala i, d'altra, el sensor ARBRES-X dissenyat per ser integrat en petits UAVs. Malgrat la seva concepció inherentment diferent, la concurrència d'ambdós sensors s'ha evidenciat al llarg d'aquest manuscrit. Aprofitant les similituds entre ells, s'han pogut avaluar de forma anàloga els dos sensors per obtenir conclusions. En aquest sentit, el vincle comú ha estat el desenvolupament del mode de funcionament polimètric OtF del RISKSAR-X, permetent que aquest sensor operi de forma equivalent a l'ARBRES-X. Pel que fa al sensor RISKSAR-X, s'han desenvolupat diverses contribucions hardware durant part d'aquest doctorat amb l'objectiu de millorar el rendiment del sistema. En dotar el sistema de la possibilitat d'operar en el mode d'adquisició totalment polarimètric OtF, s'ha reduït el relatiu llarg temps d'escaneig. Això és de gran interès ja que els blancs mesurats que presenten una reflectivitat variable a curt termini, com ara la vegetació en moviment, poden degradar els paràmetres extrets de les dades recuperades i la reconstrucció d'imatges SAR. Durant aquesta activitat doctoral s'ha estudiat el desenfocat de la imatge, la decorrelació i la degradació de la coherència introduïts per aquest efecte. A més, s'ha introduït un nou terme en la coherència interferomètrica diferencial que té en compte el desenfocat de la imatge. Pel que fa al sistema ARBRES-X, un dels principals objectius perseguits durant aquest doctorat ha estat la integració del sensor en un petit UAV MP superant les restriccions de pes, grandària i aerodinàmica de la plataforma. S'espera que l'ús d'aquest tipus de plataformes obri noves possibilitats en la detecció remota SAR aerotransportada, ja que ofereix molta més versatilitat que els UAV d'ales fixes habituals. S'han avaluat diferents estratègies de vol innovadores amb aquest tipus de plataformes i s'han obtingut resultats preliminars amb l'ús del sistema ARBRES-X. Durant el transcurs del present treball, s'ha dedicat molt esforç a assolir els primers resultats experimentals d'interferometria de múltiple passada obtinguts amb l'UAV MP conjuntament amb l'ARBRES-X. A més, el sensor ha estat dotat de capacitats totalment polarimètriques aplicant les millores desenvolupades al radar RISKSAR-X, el qual constitueix un altre exemple de la dualitat entre ambdós sistemes. Finalment, s'han realitzat amb èxit una apertura vertical i semicircular obtenint imatges SLC de l'escenari, el qual permet preveure la capacitat de l'UAV MP per a realitzar imatges tomogràfiques i apertures circulars completes en el futur. En conclusió, l'UAV MP és una plataforma prometedora que obre nous potencials per a diverses aplicacions, com ara la interferometria de múltiple passada o la tomografia diferencial amb la realització de trajectòries gairebé arbitràries.Postprint (published version