1. Puce 1999-2007: Signal Processing courses (30h) in Ecole Centrale de Paris,

  2. Puce 2003, 2005, 2007, 2013, 2015, 2017: SAE courses (3h) on Targets Radar Analysis and Recognition,

  3. Puce 2006, 2009: SAE courses (6h) on Time-Frequency and Wavelet Signals Analysis,

  4. Puce 2004-2008: ENSEA Courses (4h) on Time-Frequency Distributions Applied to Radar Imaging,

  5. Puce 2008, 2009: Student Project supervision in ENSAE,

  6. Puce 1997-2001: Courses (10h) "Wavelet and Time-Frequency Distributions" in Ecole Nationale Supérieure des Télécommunications (ENST). Editing of courses notes. Student Project Supervision,

  7. Puce 1993-1995: ESIEA Signal Processing courses (36h).

  8. Puce 2012, 2013, 2014, 2015 : Courses in Signal Processing, License Pro, IUT d’Orsay (30h)

  9. Puce 2012, 2013: Courses in Numerical Analysis, Licence, IUT Ville Avray (15h)

  10. Puce 2011-2017: Courses in Statistical and Harmonic Analysis, 2nd year, Orsay IUT (30h)


  1. Puce Supervision of Master Internship of Paris 6 University (Mélanie Mahot) from Master Recherche Systèmes de Télécommunications Numériques (STN) with Télécom ParisTech , Paris from March to June 2009 entitled "Robust STAP Techniques in Heterogeneous and non-stationary environnment"

  2. Puce Supervision of Master Internship (Pierre Formont) from Master Recherche Signal Images Parole Télécom of Gipsa-Lab, Grenoble, France from February to July 2009 entitled "Polarimetric Classification in Highly Textured SAR Images"

  3. Puce Supervision of Master Internship of Paris VII University (Frédéric Pascal) from April to July 2003 entitled "Detection in Non-Gaussian Environment",

  4. Puce Supervision of LSS Master Internship (Mohamed Tria) from April to July 2001 entitled "Wavelet SAR Analysis",

  5. Puce Supervision of ENSEA Master Internship (Eloi Thoumie) from April to July 2000 entitled "Wavelet SAR Analysis",

  6. Puce Supervision of ENSEA Master Internship (Emmanuelle Jay) from April to July 1998 entitled "Detection in Non-Gaussian Environment",

  7. Puce Supervision of EPF Master Internship (Véronique Hirsinger) from January to June 1994 entitled: "Analysis and Modelisation of Ground Clutter",

  8. Puce Supervision of ENSEA Master Internship (Christophe Corroyer) from April to July 1992 entitled: "Broad-Band Signals Analysis Using Affine Time-Frequency Distributions".


  1. Puce Supervision of Giovanni Manfredi’s Post doc (2017-2018) in ONERA/DEMR and SONDRA entitled Detection and Micro-Doppler Analysis of Humain Activities under Foliage in L and P-Band.

  2. Puce Supervision of Wajih Ben-Abdallah’s Post doc (2016-2017) in L2S/SONDRA Labs entitled Space-Time ARMA Models for Radar Applications in collaboration with Pascal Bondon from L2S.

  3. Puce Supervision of Pierre Formont’s Post doc (2013-2014) in SONDRA Lab entitled Change Detection in SAR and Hyperspectral Images.

  4. Puce Supervision of Miguel Veganzones’s Post doc (2012-2013) in Gipsa Lab entitled Robust Estimation and Detection Schemes and Classification for Hyperspectral Imaging in collaboration with Jocelyn Chanussot from Gipsa-Lab.

  5. Puce Supervision of Lionel Bombrun’s Post doc (2009-2010) from Gipsa Lab (6 months SONDRA visit) entitled Stochastic Processes for SAR Backscattering  Analysis.

  6. Puce Supervision of Gabriel Vasile’s Post doc (2007-2008) entitled "Statistical Analysis of Polarimetry in Highly textured SAR Images" (CNES grant


  1. Puce Supervision of ONERA/DGA Pierre Formont’s Ph.D. (2009-2012) entitled "Statistical and Geometric Tools for the Classification of Highly Textured Polarimetric SAR Images" in collaboration with Gabriel Vasile from GipsaLab, CNRS, Frédéric Pascal from Supélec/SONDRA and Laurent Ferro-Famil from Institut d'Electronique et de Télécommunications of Rennes University. This thesis has been defended with merit (Mention Très Honorable) in December 2013 (Defense Presentation),

      Synthetic Aperture Radars (SAR) now provide high resolution images of the Earth surface. Traditional classification algorithms are based on a Gaussian assumption for the distribution of the signal, which is no longer valid when the background is heterogeneous, which is particularly the case for polarimetric SAR images, especially in urban areas. A compound Gaussian model, called the SIRV model, allows to take into account these pheno- mena. This thesis is then devoted to studying the impact of this model for the classification of polarimetric SAR images in order to improve the interpretation of classification results in a polarimetric sense, and to propose tools better suited to this model. Indeed, classical techniques using the Gaussian assumption actually use the power information of each pixel much more than the polarimetric information. Furthermore, it is often neces- sary to compute a mean of covariance matrices, usually by taking the standard arithmetical mean. However, the space of covariance matrices has a Riemannian structure, not an Euclidean one, which means this definition of the mean is not correct. We will then present several methods to use the actual polarimetric information thanks to the SIRV model to improve the classification results. The benefit of using a correct, Riemannian definition of the mean will also be demonstrated on simulated and real data. Finally, a preliminary study of an extension of this work to hyperspectral imagery will be presented.

  1. Puce Supervision (50%) of Guilhem Pailloux’s Ph.D. (2007-2009) entitled "Noise Structured Covariance Estimation in Adaptive Detection" in collaboration with Philippe Forster (main advisor) from Ecole Normale Supérieure de Cachan (LESIR/SATIE) and Frédéric Pascal from SONDRA. This thesis has been defended with merit (Mention Très Honorable) in June 2010 (Defense presentation),

        This thesis deals with Radar detection in Gaussian and non-Gaussian noise. In this context, the clutter covariance matrix commonly exhibits a particular persymmetric structure. This structure is exploited into a particular matrix transformation to provide two new covariance matrices estimates for gaussian and non-gaussian noise. We use then this particular linear transformation in order to develop and to study the statistical property of the two new detectors based on these estimates for both Gaussian and non-Gaussian environments. The improvement in terms of detection performances of these new detectors is shown through a lot of simulations and validation on operational data, for both Gaussian and non-Gaussian noise. Moreover this exploitation is extended to space-time adaptive processing and reduced rank technical. All the results confirm then the high interest of taking into account this particular structure in radar detection process compared to classical detection schemes. The case of Toeplitz matrices is also studied. The Toeplitz matrices are a particular class of structured matrices obtained with space-time processing which theoretically allows improving the performance of detectors based on this matrix assumption. In this context, some preliminary results are presented in order to conclude this thesis.

  1. Puce Supervision (50%) of Frédéric Pascal’s Ph.D. (2003-2006) entitled "Detection and Estimation in non-Gaussian Environment" in collaboration with Philippe Forster (main advisor) from GEA (Groupe d’Electromagnétisme Appliqué) and Paris X University and Pascal Larzabal from Ecole Normale Supérieure de Cachan (LESIR/SATIE). This thesis has been defended with Merit (Très Honorable avec les félicitations du jury) in December 2006 in Paris X University (Defense presentation),

This thesis deals with radar detection in impulsive noise contexts. Indeed, under Gaussian assumptions, classical detectors, like Optimum Gaussian Detector, present several failures when the background scatterers are heterogeneous and non-Gaussian, which is the case for ground or sea clutter. Clutter physical models based on compound noise modeling (SIRP, Compound Gaussian Processes) allow to correctly describe the reality (range power variations or clutter transitions areas). However, these models depend on several unknown parameters (covariance matrix, statistical distribution of the texture, disturbance parameters) which have to be estimated. When this estimation scheme is carried out, it is possible to build optimum radar detectors (Generalized Likelihood Ratio Test - Linear Quadratic) associated to this non-Gaussian background. Based on these noise models, this thesis presents a complete analysis of several estimation schemes of the noise covariance matrix, associated to the detection problem. A statistical study of the main covariance matrix estimates used in the literature, is performed. Moreover, an improved estimate is proposed: the Fixed Point estimate, very attractive thanks to its good statistical and detection properties. This thesis also describes detection performance and theoretical properties (texture-CFAR and matrix-CFAR) of the GLRT-LQ detector built with the studied covariance matrix estimates. In particular, the detector invariance to the texture distribution and to the covariance matrix structure is shown. Finally, these detectors are analyzed with simulated data and then, experimented on real ground clutter data.

  1. Puce Supervision (70%) of Mohamed Tria’s Ph.D. (2001-2003) entitled "Analysis of SAR Images using Continuous Multidimensional Wavelet Transform" in collaboration with Messaoud Benidir (main advisor) from Laboratoire des Signaux et Systèmes (LSS) of Supelec. This thesis has been defended (Mention Très Honorable) in November 2005 in Paris XI University (Defense presentation),

The classical SAR (Synthetic Aperture Radar) imaging methods assume that all the reflectors  are isotropic (they behave the same way regardless the aspect under which they are viewed) and white (they have the same properties within the emitted bandwidth). Unfortunately, new  capacities in SAR imaging (large bandwidth, large angular excursions of analysis) make these  assumptions obsolete. The purpose of the thesis work consisted in developing a time-frequency method based on the use of continuous wavelets allowing to highlight the anisotropic and coloration effects of the irradiated points. A great part of the work consisted in exploiting the coloration and the anisotropy of the reflectors to discriminate a target in a SAR image. Finally, the continuous wavelet analysis has been extended to polarimetric in order to improve the classification of random targets and deterministic ones then to interferometry with the aim of obtaining a better estimation of target height estimation.

  1. Puce Supervision (50%) of Emmanuelle Jay’s Ph.D. (1998-2001) entitled "Detection in non-Gaussian Environment" in collaboration with Patrick Duvaut (main advisor) and David Declercq from ETIS Lab in ENSEA. This thesis has been defended with Merit (Très Honorable avec les félicitations du jury) in June 2002 in Cergy-Pontoise University (Defense presentation).

For a long time, radar echoes coming from the various returns of the transmitted signal on many objects of the environment (clutter) have been exclusively modelled by Gaussian vectors. The related optimal detection procedure was then performed by the classical matched filter. Then, the technological improvement of radar systems showed that the true nature of the clutter could not be considered as Gaussian any more. Though the optimality of the matched filter is no more valid in such cases, CFAR techniques (Constant False Alarm Rate) were proposed for this detector in order to adapt the value of the detection threshold to the multiple local variations of the clutter. In spite of their diversity, none of these techniques turned to be either robust or optimal in these situations. With the modelling of the clutter by non-Gaussian complex processes, such as SIRP (Spherically Invariant Random Process), optimal structures of coherent detection have been found. These models describe many non-Gaussian laws, like K-distribution or Weibull laws, and are acknowledged in the literature to model many experimental situations in a relevant way. To identify the law of their characteristic component (namely the texture) without statistical a priori on the model, we propose, in this thesis, to tackle the problem by a Bayesian approach. Two new estimation methods of the texture law emerge from this proposition: the first one is a parametric method, based on a Padé approximation of the moment generating function, and the second one results from a Monte Carlo estimation. These estimations are carried out on reference clutter data and lead to two new optimal detection strategies, respectively named PEOD (Padé Estimated Optimum Detector) and BORD (Bayesian Optimum Detector Radar). The asymptotic expression of the BORD (convergence in law), called the "Asymptotic BORD", is established together with its law. This last result gives access to the optimal theoretical performances of the Asymptotic BORD, and may also applied to the BORD if the data correlation matrix is non-singular. The detection performances of BORD and those of Asymptotic BORD are evaluated on experimental ground clutter data. We obtained results that validate both the relevance of SIRP model for the clutter, the optimality of the BORD and its adaptability to any type of environment.

  1. Puce Supervision (co-direction 20%) of Mélanie Mahot’s Ph.D. (2009-2012) entitled "Robust Covariance Matrix Estimation in Signal Processing" in collaboration with Philippe Forster (main advisor) from Ecole Normale Supérieure de Cachan (LESIR/SATIE) and with Frédéric Pascal from SONDRA. This thesis has been defended with merit (Mention Très Honorable) in December 2012 in Ecole Normale Supérieure of Cachan, (Defense Presentation).

In many signal processing applications, the covariance matrix of the received data must be known. If unknown, it is firstly estimated with some training data. Classically, the background is considered as Gaussian. In such a case, the maximum likelihood estimator is the Sample Covariance Matrix (SCM). However, due to high resolution methods or other new technics, the Gaussian assumption is not valid anymore. Moreover, even when the data are Gaussian, the SCM can be strongly influenced by some disturbances such as missing data and/or outliers. In this thesis, we use a more general model which encompasses a large panel of distributions: the elliptical distributions. Many campagns of measurement have shown that this model leads to a better modelling of the data. In this context, we present more robust and adapted estimators: the M-estimators and Fixed Point Estimator (FPE). Their properties are derived in terms of performance and robustness, and they are compared to the SCM. We show that these estimators can be used instead of the SCM with nearly the same performance when the data are Gaussian, and better performance when the data are non-Gaussian. Theoretical results are emphasized on simulations and on real data in a context of Space Time Adaptive Processing. 

SAR Imaging considers a target as a set of scatterers, isotropic and non-dispersive. This assumption becomes obsolete for new applications radar using a large band and a large angular spreading. The aim is to study the non-stationary behavior by using the basic tool: the time-frequency analysis. The usual time-frequency transforms allow to build hyperimages (images in 4 dimensions) describing the energy versus emitted frequency and observation angle. This information source is related to radar cross-section (RCS). Another information source in radar imaging is polarimetry.The goal is to jointly use these two information sources to characterize scatterers. By applying continuous wavelet and polarimetric coherent decompositions on the Sinclair matrix, some polarimetric hyperimages are built. These hyperimages allow to describe the polarimetric behavior versus the observation angle and the emitted frequency. In some case, they globally describe the backscattering mecanism by its nature, its relative orientation in the horizontal plan, its absolute orientation in the vertical plan. Moreover, they highlight a new phenomena: the polarimetric non-stationary. This would be caused by the fact the radar does not see the same geometry during the processing. From polarimetric hyperimages, a new phenomenological classification is proposed.

  1. Puce Supervision (50%) of Mickaël Duquenoy’s Ph.D. (2004-2007) entitled "Time-Frequency Analysis Applied to Polarimetric SAR Imaging" in collaboration with Eric Pottier (main advisor) and Laurent Ferro-Famil from Institut d'Electronique et de Télécommunications of Rennes University, Rennes. This thesis has been defended with merit (Mention Très Honorable) in October 2009 in IETR, Rennes, France (Defense presentation),

  1. Puce Thesis Director of Joana Frontera-Pons’s Ph.D. (2011-2014) entitled "Robust Detection and Classification for Hyperspectral Imaging" in collaboration with Frédéric Pascal from Supélec/SONDRA. This thesis has been defended with merit (Mention Très Honorable) in December 2014 in University of Paris-Saclay (Supelec), (Defense presentation).

  1. Puce Participation to the supervision (20%) of Julien Totems’s Ph.D. (2007-2010) entitled "Pulsed laser emission and time-frequency processing for coherent detection lidar vibrometry" in collaboration with Nadine Martin (main advisor) from Gipsa Lab, CNRS and Véronique Jolivet (ONERA/DOTA/SLS). This thesis has been defended with merit (Mention Très Honorable) in February 2011 in the University of Grenoble.

The use of pulsed lasers could lead to new functionnalities and increased compacity of lidar systems for remote vibration sensing. However, specific amplitude and phase noises affect the lidar signal, and particularly decrease the performance of a polypulse based emission regime, thought to be promising for very long range measurements.

This work first deals with the experimental characterization of these noise sources in order to properly model them, with a focus on atmospheric turbulence. We then seek to optimize the employed waveform and signal processing in regard of the vibration and noise conditions. An original method is proposed that involves maximum likelihood based estimation of the vibration-induced Doppler shift, and its extraction from a time-frequency representation of the signal. The benefits of this approach are shown in simulation and experimentation, by comparing the performance of various emission modes.

  1. Puce Participation to the supervision of Supelec Chong Chin Yuan’s Ph.D. (2008-2011) entitled "Signal Processing for MIMO Radars : Detection Under Gaussian and non-Gaussian Environments and Application to STAP" in collaboration with Marc Lesturgie from ONERA and Frédéric Pascal from Supélec/SONDRA. This thesis has been defended in University of Paris Sud with merit (Mention Très Honorable) in November 2011 (Defence Presentation),

      A Multiple-Input Multiple Output (MIMO) radar can be broadly defined as a radar system employing multiple transmit waveforms and having the ability to jointly process signals received at multiple receive antennas. In terms of configurations, the antennas can be widely separated or co-located. The first part of the thesis is on detection under Gaussian and non- Gaussian environments using a MIMO radar which contains several widely separated subarrays with one or more elements each. Two different situa- tions are considered. Firstly, we consider that the interference is Gaussian but correlation between subarrays can arise due to insufficient spacing and the imperfect orthogonality of waveforms. Secondly, we consider that the interference is non-Gaussian, a situation which arises under sea and ground clutter and when the resolution is very high. The second part is on the application of MIMO techniques to Space- Time Adaptive Processing (STAP). The coherent MIMO configuration is studied in terms of antenna element distribution and inter-element spac- ing to improve detection and estimation performance. A preliminary study is also done on the use of spatial diversity to improve detection stability w.r.t. target Radar Cross Section (RCS) fluctuations and velocity direction changes.

Hyperspectral imaging (HSI) extends from the fact that for any given material, the amount of emitted radiation varies with wavelength. HSI sensors measure the radiance of the materials within each pixel area at a very large number of contiguous spectral bands and provide image data containing both spatial and spectral information. Classical adaptive detection schemes assume that the background is zero-mean Gaussian or with known mean vector that can be exploited. However, when the mean vector is unknown, as it is the case for hyperspectral imaging, it has to be included in the detection process. We propose in this work an extension of classical detection methods when both covariance matrix and mean vector are unknown.

However, the actual multivariate distribution of the background pixels may differ from the generally used Gaussian hypothesis. The class of elliptical distributions has already been popularized for background characterization in HSI. Although these non-Gaussian models have been exploited for background modeling and detection schemes, the parameters estimation (covariance matrix, mean vector) is usually performed using classical Gaussian-based estimators. We analyze here some robust estimation procedures (M-estimators of location and scale) more suitable when non-Gaussian distributions are assumed. Jointly used with M-estimators, these new detectors allow to enhance the target detection performance in non-Gaussian environment while keeping the same performance than the classical detectors in Gaussian environment. Therefore, they provide a unified framework for target detection and anomaly detection in HSI.

  1. Puce Thesis Director of Uyhour Tan’s Ph.D. (2014-2017) entitled "Colocated MIMO Radar Waveform Optimization" in collaboration with Jean-Paul Guivarch, Claude Adnet from THALES Air Systems and Olivier Rabaste  from ONERA/DEMR/TSI.

  1. Puce Thesis Director of Eugénie Terreaux’s Ph.D. (2015-2018) entitled "Random Matrix Theory for Hyperspectral Imaging" in collaboration with Frédéric Pascal (CentraleSupelec/L2S). This thesis has been defended with merit (Mention Très Honorable) in November 2018 in University of Paris-Saclay (CentraleSupélec). 2018 Best PhD thesis prize in ONERA Physic Domain.

     Hyperspectral imaging generates large data due to the spectral and spatial high resolution, as it is the case for more and more other kinds of applications. For hyperspectral imaging, the data complexity comes from the spectral and spatial heterogeneity, the non-gaussianity of the noise and other physical processes. Nevertheless, this complexity enhances the wealth of collected informations, that need to be processed with adapted methods. Random matrix theory and robust processes are here suggested for hyperspectral imaging application: the random matrix theory is adapted to large data and the robustness enables to better take into account the non-gaussianity of the data. This thesis aims to enhance the model order selection on a hyperspectral image and the unmixing problem. As the model order selection is concerned, three new algorithms are developped, and the last one, more robust, gives better performances. One financial application is also presented. As for the unmixing problem, three methods that take into account the peculierities of hyperspectral imaging are suggested. The random matrix theory is of great interest for hyperspectral image processing, as demonstrated in this thesis. Differents methods developped here can be applied to other field of signal processing requiering the processing of large data.

  1. Puce Thesis Director of Ahmad Bitar’s Ph.D. (2015-2018) entitled "Exploitation of Sparsity for Hyperspectral Target Detection" en collaboration with Cheong Loong Fah from the National University of Singapore. This thesis has been defended with merit (Mention Très Honorable) in June 2018 in University of Paris-Saclay (CentraleSupélec)

       The title of this PhD thesis is formed by three keywords: sparsity, hyperspectral image, and target detection. Sparsity is a word that is used everywhere and in everyday life. It generally means « small in number or amount, often spread over a large area». A hyperspectral image is a three dimensional data cube consisting of a series of images of the same spatial scene in a contiguous and multiple narrow spectral wavelength (color) bands. According to the high spectral dimensionality, target detection is not surprisingly one of the most important applications in hyperspectral imagery. The main objective of this PhD thesis is to answer the question « How and Why can sparsity be exploited for hyperspectral target detection? ». Answering this question has allowed us to develop different target detection methods that mainly take into consideration the heterogeneity of the environment, the fact that the total image area of all the targets is very small relative to the whole image, and the estimation challenge of the covariance matrix (surrounding the test pixel) in large dimensions. The proposed mehods are evaluated on both synthetic and real experiments, the results of which demonstrate their effectiveness for hyperspectral target detection.

  1. Puce Thesis Director of Ammar Mian (2016-2019) entitled "Exploitation of SAR and Hyperspectral Time Series" in collaboration with Guillaume Ginolhac (co-director) and Abdourrahmane Atto from University of Savoie Mont-Blanc (LISTIC).