Radar Imaging in Challenging Scenarios from Smart and Flexible Platforms

undefine

Ultralight Radar Sensor for Autonomous Operations by Mini- and Micro-UAS

In recent years the boost in operations by mini- and micro-UAS (Unmanned Aircraft Systems, also known as Remotely Piloted Aircraft Systems - RPAS - or simply drones) and the successful miniaturization of electronic components were experienced. Radar sensors demonstrated to have favorable features for these operations. However, despite their ability to provide meaningful information for navigation, sense-and-avoid, and imaging tasks, currently very few radar sensors are exploited onboard or developed for autonomous operations with mini- and micro-UAS. Exploration of indoor complex, dangerous, and not easily accessible environments represents a possible application for mini-UAS based on radar technology. In this scenario, the objective of the thesis is to develop design strategies and processing approaches for a novel ultralight radar sensor able to provide the miniaturized platform with Simultaneous Localization and Mapping (SLAM) capabilities, mainly but not exclusively indoors. Millimeter-wave Interferometric Synthetic Aperture Radar (mmw InSAR) technology has been identified as a key asset. At the same time, testing of commercial lightweight radar is carried out to assess potentialities towards autonomous navigation, sense-and-avoid, and imaging. The two main research lines can be outlined as follows: - Long-term scenario: Development of very compact and ultralight Synthetic Aperture Radar able to provide mini- or micro-UAS with very accurate 3D awareness in indoor or GPS-denied complex and harsh environments. - Short-term scenario: Assessment of true potentialities of current commercial radar sensors in a UAS-oriented scenario. Within the framework of long-term scenario, after a review of state-of-art SAR sensors, Frequency-Modulated Continuous Wave (FMCW) SAR technology has been selected as preferred candidate. Design procedure tailored to this technology and software simulator for operations have been developed in MATLAB environment. Software simulator accounts for the analysis of ambiguous areas in a three-dimensional environment, different SAR focusing algorithms, and a Ray-Tracing algorithm specifically designed for indoor operations. The simulations provided relevant information on actual feasibility of the sensor, as well as mission design characteristics. Additionally, field tests have been carried out at Fraunhofer Institute FHR with a mmw SAR. Processing approaches developed from simulations proved to be effective when dealing with field tests. A very lightweight FMCW radar sensor manifactured by IMST GmbH has been tested for short-term scenario operations. The codes for data acquisition were developed in Python language both for Windows-based and GNU/Linux-based operative systems. The radar provided information on range and angle of targets in the scene, thus being interesting for radar-aided UAS navigation. Multiple-target tracking and radar odometry algorithms have been developed and tested on actual field data. Radar-only odometry provided to be effective under specific circumstances

Analyse temps-frequence et traitement des signaux RSO à haute résolution spatiale pour la surveillance des grands ouvrages d'art

The thesis is composed of two research axis. The first one consists in proposing time-frequency signal processing tools for frequency modulated continuous wave (FMCW) radars used for displacements measurements, while the second one consists in designing a spaceborne synthetic aperture radar (SAR) signal processing methodology for infrastructure monitoring when an external point cloud of the envisaged structure is available. In the first part of the thesis, we propose our solutions to the nonlinearity problem of an X-band FMCW radar designed for millimetric displacement measurements of short-range targets. The nonlinear tuning curve of the voltage controlled oscillator from the transceiver can cause a dramatic resolution degradation for wideband sweeps. To mitigate this shortcoming, we have developed two time warping-based methods adapted to wideband nonlinearities: one estimates the nonlinear terms using the high order ambiguity function, while the other is an autofocus approach which exploits the spectral concentration of the beat signal. Onwards, as the core of the thesis, we propose a novel method for scattering centers detection and tracking in spaceborne SAR images adapted to infrastructure monitoring applications. The method is based on refocusing each SAR image on a provided 3D point cloud of the envisaged infrastructure and identifying the reliable scatterers to be monitored by means of four dimensional (4D) tomography. The refocusing algorithm is compatible with stripmap, spotlight and sliding spotlight SAR images and consists of an azimuth defocusing followed by a modified back-projection algorithm on the given set of points which exploits the time-frequency structure of the defocused azimuth signal. The scattering centers of the refocused image are detected in the 4D tomography framework by testing if the main response is at zero elevation in the local elevation-velocity spectral distribution. The mean displacement velocity is estimated from the peak response on the zero elevation axis, while the displacements time series for detected single scatterers is computed as double phase difference of complex amplitudes.Finally, we present the measurement campaigns carried out on the Puylaurent water-dam and the Chastel landslide using GPS measurements, topographic surveys and laser scans to generate the point clouds of the two structures. The comparison between in-situ data and the results obtained by combining TerraSAR-X data with the generated point clouds validate the developed SAR signal processing chain.Cette thèse s'articule autour de deux axes de recherche. Le premier axe aborde les aspects méthodologiques liés au traitement temps-fréquence des signaux issus d'un radar FMCW (à onde continue modulée en fréquence) dans le contexte de la mesure des déplacements fins. Le second axe est dédié à la conception et à la validation d'une chaîne de traitement des images RSO (radar à synthèse d'ouverture) satellitaire. Lorsqu'un maillage 3D de la structure envisagée est disponible, les traitements proposés sont validés par l'intercomparaison avec les techniques conventionnelles d'auscultation des grands ouvrages d'art.D'une part, nous étudions la correction de la non-linéarité d'un radar FMCW en bande X, à courte portée, conçu pour la mesure des déplacements millimétriques. La caractéristique de commande non linéaire de l'oscillateur à large bande, entraine une perte de résolution à la réception. Afin de pallier cet inconvénient, nous avons développé deux méthodes basées sur le ré-échantillonnage temporel (time warping) dans le cas des signaux à large bande non-stationnaires. La première approche estime la loi de fréquence instantanée non linéaire à l'aide de la fonction d'ambiguïté d'ordre supérieur, tandis que la deuxième approche exploite la mesure de concentration spectrale du signal de battement dans un algorithme d'autofocus radial.D'autre part, nous proposons un cadre méthodologique général pour la détection et le pistage des centres de diffusion dans les images RSO pour la surveillance des grands ouvrages d'art. La méthode est basée sur la ré-focalisation de chaque image radar sur le maillage 3D de l'infrastructure étudiée afin d'identifier les diffuseurs pertinents par tomographie 4D (distance – azimut – élévation – vitesse de déformation). L'algorithme de ré-focalisation est parfaitement compatible avec les images RSO acquises dans les différents modes (« stripmap », « spotlight » et « sliding spotlight ») : dé-focalisation en azimut suivie par rétroprojection modifiée (conditionnée par la structure temps-fréquence du signal) sur l'ensemble donné des points. Dans la pile d'images ré-focalisées, les centres de diffusion sont détectés par tomographie 4D : test de conformité à l'hypothèse d'élévation zéro dans le plan élévation – vitesse de déformation. La vitesse moyenne correspond au maximum à l'élévation zéro, tandis que la série temporelle des déplacements est obtenue par double différence de phase des amplitudes complexes pour chaque diffuseur pertinent.Nous présentons également les campagnes in situ effectuées au barrage de Puylaurent (et glissement de Chastel) : les relevés GPS, topographiques et LIDAR sol employées au calcul des maillages 3D. La comparaison entre les déplacements mesurés in situ et les résultats obtenus par l'exploitation conjointe de la télédétection RSO satellitaires et les maillages 3D valident la chaîne de traitement proposée.Teza cuprinde două axe principale de cercetare. Prima axă abordează aspecte metodologice de prelucraretimp-frecvenţă a semnalelor furnizate de radare cu emisie continuă şi modulaţie de frecvenţă (FMCW)în contextul măsurării deplasărilor milimetrice. În cadrul celei de-a doua axe, este proiectată şi validatăo metodă de prelucrare a imaginilor satelitare SAR (radar cu apertură sintetică) ce este destinatămonitorizării infrastructurii critice şi care se bazează pe existenţa unui model 3D al structurii respective.În prima parte a tezei, sunt investigate soluţii de corecţie a neliniarităţii unui radar FMCW în bandaX destinat măsurării deplasărilor milimetrice. Caracteristica de comandă neliniară a oscilatorului debandă largă determină o degradare a rezoluţiei în distanţă. Pentru a rezolva acest inconvenient, au fostelaborate două metode de corecţie a neliniarităţii, adaptate pentru semnale de bandă largă, ce se bazeazăpe conceptul de reeşantionare neuniformă sau deformare a axei temporare. Prima abordare estimeazăparametrii neliniarităţii utilizând funcţii de ambiguitate de ordin superior, iar cea de-a doua exploateazăo măsură de concentraţie spectrală a semnalului de bătăi într-un algoritm de autofocalizare în distanţă.În a doua parte a lucrării, este propusă o metodologie generală de detecţie şi monitorizare a centrilorde împrăştiere în imagini SAR în scopul monitorizării elementelor de infrastructură critică. Metoda sebazează pe refocalizarea fiecărei imagini radar pe un model 3D al structurii investigate în scopul identificăriicentrilor de împrăştiere pertinenţi (ţinte fiabile ce pot fi monitorizate în timp) cu ajutorul tomografiei SAR4D (distanţă-azimut-elevaţie-viteză de deplasare). Algoritmul de refocalizare este compatibil cu imaginiSAR achiziţionate în moduri diferite (« stripmap », « spotlight » şi « sliding spotlight ») şi constă într-odefocalizare în azimut urmată de o retroproiecţie modificată (condiţionată de structura timp-frecvenţă asemnalului) pe modelul 3D al structurii. Ţintele sunt identificate în stiva de imagini refocalizate cu ajutorultomografiei 4D prin efectuarea unui test de conformitate cu ipoteza că centrii de împrăştiere pertinenţivor avea elevaţie zero în planul local elevaţie-viteză. Viteza medie de deformare corespunde maximuluide pe axa de elevaţie nulă, iar seria temporară a deplasărilor se obţine printr-o dublă diferenţă de fază aamplitudinilor complexe corespunzătoare ţintelor identificate.În final sunt prezentate campaniile de măsurători pe teren efectuate la un baraj şi o alunecare de terendin regiunea Puylaurent (Franţa) destinate obţinerii modelului 3D al celor două elemente de infrastructurăprin măsurători GPS, topografice şi LIDAR. Comparaţia între deformările măsurate pe teren şi rezultateleobţinute prin combinarea imaginilor SAR cu modelele 3D au permis validarea metodologiei propuse

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

Radar Technology

In this book “Radar Technology”, the chapters are divided into four main topic areas: Topic area 1: “Radar Systems” consists of chapters which treat whole radar systems, environment and target functional chain. Topic area 2: “Radar Applications” shows various applications of radar systems, including meteorological radars, ground penetrating radars and glaciology. Topic area 3: “Radar Functional Chain and Signal Processing” describes several aspects of the radar signal processing. From parameter extraction, target detection over tracking and classification technologies. Topic area 4: “Radar Subsystems and Components” consists of design technology of radar subsystem components like antenna design or waveform design

Study of the speckle noise effects over the eigen decomposition of polarimetric SAR data: a review

This paper is focused on considering the effects of speckle noise on the eigen decomposition of the co- herency matrix. Based on a perturbation analysis of the matrix, it is possible to obtain an analytical expression for the mean value of the eigenvalues and the eigenvectors, as well as for the Entropy, the Anisotroopy and the dif- ferent a angles. The analytical expressions are compared against simulated polarimetric SAR data, demonstrating the correctness of the different expressions.Peer ReviewedPostprint (published version

Subsurface radar imaging from space

© Cranfield University, 2018Ground Penetrating Radar (GPR) and Synthetic Aperture Radar (SAR) are two widely used techniques for acquiring radar images. GPR, as its name suggests, produces radar images of the below ground environment. SAR is a remote sensing technique which allows moving radar systems to produce radar images with dramatically improved resolutions over conventional radar systems. Despite their benefits, both GPR and SAR suffer from certain limitations. In the case of GPR, the radar system has to be in close proximity with the subsurface volume being surveyed, which limits the process to relatively small areas that are easily accessible. SAR allows large areas to be surveyed rapidly from large distances, but cannot distinguish buried objects from surface objects. This thesis focuses on a radar technique that offers the opportunity to overcome these limitations and allow subsurface radar imaging of large areas using radar data gathered by remote sensing systems. This novel technique is known as Virtual Bandwidth SAR (VB-SAR). VB-SAR utilises changes in soil moisture over a series of SAR images to differentiate buried objects from objects on the surface. In addition to this differentiation, VB-SAR also allows extremely high (centimetre scale) subsurface range resolutions to be obtained from SAR images with range resolutions measured in metres. This research has experimentally demonstrated the basic feasibility of performing remote subsurface radar imaging with the VB-SAR scheme. Within the laboratory environment a buried target has been successfully imaged using VB-SAR and the fundamentals of VB-SAR have been verified. Dramatic increases in subsurface range resolutions have been demonstrated, as has the ability of the VB-SAR scheme to work correctly over a range of radar frequencies, observation angles and polarisations. This laboratory work has been enabled by use of the Tomographic Profiling (TP) imaging scheme. TP is a synthetic aperture based imaging algorithm, but unlike conventional SAR TP produces images with a constant look angle over the entire imaging scene. This enabled the performance of the VB-SAR imaging scheme to be easily evaluated over a range of look angles using a single radar dataset and simplified the experimental setup. In addition to the experimental work, simulation exercises have been conducted and image processors have been implemented. Simulation, using a simulator created as part of this work, has allowed testing of the VB-SAR scheme in a range of scenarios (sidelooking SAR, different soils, multiple buried targets). The image processor work has implemented a high performance TP processor and a practical VB-SAR imager