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

Algebraic Representation and Geometric Interpretation of Hilbert Transformed Signals

This thesis covers a fundamental problem of local phase based signal processing: the isotropic generalization of the classical one dimensional analytic signal (D. Gabor) to higher dimensional signal domains. The classical analytic signal extends a real valued one dimensional signal to a complex valued signal by means of the classical 1D Hilbert transform. This signal extension enables the complete analysis of local phase and local amplitude information for each frequency component in the sense of Fourier analysis. In case of two dimensional signal domains, e.g. for images, additional geometric information is required to characterize higher dimensional signals locally. The local geometric information is called orientation, which consists of the main orientation and apex angle for two superimposed one dimensional signals. The problem of two dimensional signal analysis is the fact that in general those signals could consist of an unlimited number of superimposed one dimensional signals with individual orientations. Local phase, amplitude and additional orientation information can be extracted by the monogenic signal (M. Felsberg and G. Sommer) which is always restricted to the subclass of intrinsically one dimensional signals, i.e. the class of signals which only make use of one degree of freedom within the embedding signal domain. In case of 2D images the monogenic signal enables the rotationally invariant analysis of lines and edges. In contrast to the 1D analytic signal the monogenic signal extends all real valued signals of dimension n to a (n+1) - dimensional vector valued monogenic signal by means of the generalized first order Hilbert transform, which is also known as the Riesz transform. The analytic signal and the monogenic signal show that a direct relation between analytical signals and their algebraic representation exists. This fact has motivated the work and the results of this thesis, namely the extension of the 2D monogenic signal to more general 2D analytic signals, their algebraic representation, and their most geometric embedding. In case of more general 2D signals the geometric algebra will be shown to be a natural representation, and the conformal space as the geometric embedding for the signal interpretation. In this thesis we present 2D analytic signals as generalizations of the 2D monogenic signal which now extend the original 2D signal to a multi-vector valued signal in homogeneous conformal space by means of higher order Hilbert transforms, and by means of a so called hybrid matrix geometric algebra representation. The 2D analytic signal and the more general multi-vector signal will be interpreted in conformal space which delivers a descriptive geometric interpretation and algebraic embedding of signals. In case of 2D image signals the 2D analytic signal and the multi-vector signal enable the rotationally invariant analysis of lines, edges, corners and junctions in one unified framework. Furthermore, additional local curvature can be determined by first order generalized Hilbert transforms without the need of derivatives. This so called conformal monogenic signal can be defined for any signal domain