L’observation de la Terre par des systèmes d’acquisition d’images permet de suivre l’évolution temporelle de phénomènes naturels tels que les séismes, les volcans ou les mouvements gravitaires. Différentes techniques existent dont l’imagerie satellitaire, la photogrammétrie terrestre et les mesures in-situ. Les séries temporelles d’images issues d’appareils photo automatiques (Time Lapse) sont une source d’informations en plein essor car elles offrent un compromis intéressant en termes de couverture spatiale et de fréquence d’observation pour mesurer les déplacements de surface de zones spécifiques. Cette thèse est consacrée à l’analyse de séries d’images issues de la photographie terrestre et de l’imagerie radar satellitaire pour la mesure du déplacement des glaciers Alpins. Nous nous intéressons en particulier aux problèmes du traitement de Time Lapse stéréo pour le suivi d’objets géophysiques dans des conditions terrain peu favorables à la photogrammétrie. Nous proposons une chaîne de traitement mono-caméra qui comprend les étapes de sélection automatique des images, de recalage et de calcul de champs de déplacement bidimensionnel (2D). L’information apportée par les couples stéréo est ensuite exploitée à l’aide du logiciel MICMAC pour reconstruire le relief et obtenir le déplacement tridimensionnel(3D). Plusieurs couples d’images radar à synthèse d’ouverture (SAR) ont également été traités à l’aide des outils EFIDIR pour obtenir des champs de déplacement 2D dans la géométrie radar sur des orbites ascendantes ou descendantes. La combinaison de mesures obtenues quasi-simultanément sur ces deux types d’orbites permet de reconstruire le déplacement 3D. Ces méthodes ont été mises en oeuvre sur des séries de couples stéréo acquis par deux appareils photo automatiques installés sur la rive droite du glacier d’Argentière et sur des images du satellite TerraSAR-X couvrant le massif du Mont-Blanc. Les résultats sont présentés sur des données acquises lors d’une expérimentation multi-instruments menée en collaboration avec l’IGN à l’automne 2013, incluant le déploiement d’un réseau de Géocubes qui ont fournit des mesures GPS. Elles sont utilisées pour évaluer la précision des résultats obtenus par télédétection proximale et spatiale sur ce type de glacier. / Earth observation by image acquisition systems allows the survey of temporal evolution of natural phenomena such as earthquakes, volcanoes or gravitational movements. Various techniques exist including satellite imagery, terrestrial photogrammetry and in-situ measurements. Image time series from automatic cameras (Time Lapse) are a growing source of information since they offer an interesting compromise in terms of spatial coverage and observation frequency in order to measure surface motion in specific areas. This PhD thesis is devoted to the analysis of image time series from terrestrial photography and satellite radar imagery to measure the displacement of Alpine glaciers. We are particularly interested in Time Lapse stereo processing problems for monitoring geophysical objects in unfavorable conditions for photogrammetry. We propose a single-camera processing chain that includes the steps of automatic photograph selection, coregistration and calculation of two-dimensional (2D) displacement field. The information provided by the stereo pairs is then processed using the MICMAC software to reconstruct the relief and get the three-dimensional (3D) displacement. Several pairs of synthetic aperture radar (SAR) images were also processed with the EFIDIR tools to obtain 2D displacement fields in the radar geometry in ascending or descending orbits. The combination of measurements obtained almost simultaneously on these two types of orbits allows the reconstruction of the 3D displacement. These methods have been implemented on time series of stereo pairs acquired by two automatic cameras installed on the right bank of the Argentière glacier and on TerraSAR-X satellite images covering the Mont-Blanc massif. The results are presented on data acquired during a multi-instrument experiment conducted in collaboration with the French Geographic National Institute (IGN) during the fall of 2013,with a network of Géocubes which provided GPS measurements. They are used to evaluate the accuracy of the results obtained by proximal and remote sensing on this type of glacier.

A quase inacessibilidade de grandes partes do Oceano Austral torna o conhecimento da espessura do gelo marinho limitado. Esta informação é essencial para a determinação do balanço de massa deste componente da criosfera. Na transição do inverno para a primavera de 2006, uma equipe de pesquisadores, coletou uma série de perfis de espessura de gelo marinho no norte e noroeste do mar de Weddell. Eles estavam a bordo do navio de pesquisa alemão Polarstern do Alfred-Wegener-Institut, Helmholtz-Zentrum für Polar- und Meeresforschung (AWI) e utilizaram um sistema composto por um sensor eletromagnético, altímetro laser e um Sistema de Posicionamento Global Diferencial (DGPS) aerotransportado por helicóptero. Simultaneamente, a Agência Espacial Europeia (ESA), adquiriu imagens ENVISAT ASAR WSM da área de pesquisa. O objetivo principal desta tese é determinar a relação existente entre o retroespalhamento medido pelo Radar de Abertura Sintética e a espessura do gelo marinho obtida pelo HEM (Helicopter-borne ElectroMagnetic system) usando os dados citados acima. Utilizamos os programas de computador NEST®, MATLAB®, EXCEL®, ArcGIS®/ArcMAP®, para processar, analisar e selecionar as imagens, para determinar a relação entre o retroespalhamento e as medidas, quase-tempo-coincidentes, de espessura do gelo. Projetamos as trajetórias dos voos sobre as imagens obtidas nas mesmas datas e extraímos os pixels referentes aos locais onde foram medidas as espessuras de gelo marinho. Apropriamos os valores de espessura do gelo para a área de cada pixel sobre o qual se referiam. Uma análise estatística determinou que o parâmetro que melhor representa a espessura do gelo dentro da área do pixel é a média. A regressão linear é a melhor forma de ajuste das relações entre o valor de retroespalhamento do pixel e a espessura do gelo marinho contido na área deste pixel. O coeficiente de correlação linear de Pearson, resultante de análise paramétrica, indica uma forte correlação (0,75) entre retroespalhamento e espessura do gelo marinho. Porém, a análise não paramétrica de Spearman resultou em um coeficiente de correlação baixo (0,06) o que pode indicar que os dados analisados são compostos por duas populações distintas (e.g., gelo de primeiro ano e plurianual). Contudo, a análise não paramétrica de Kolmogorov-Smirnov aventa a possibilidade de que não tenhamos amostrado toda população. Esta seria a razão de não haver valores de espessuras relativas ao intervalo entre -9.21dB e -1.35dB, o que poderia induzir ao baixo valor do coeficiente de correlação na análise de Spearman. Aplicamos a equação linear: y=0,6345x+12,015 às imagens e pudemos separá-las em doze classes: uma para água e onze para gelo marinho, estas com intervalos de um metro. Isto possibilitou a observação de importantes feições oceanográficas como: canais de águas abertas; cristas de compressão, decaimento do gelo marinho, deslocamento de icebergs (como indicadores de correntes, marés e ventos), liberação de gelo por geleiras de maré, desprendimento de iceberg da plataforma de gelo Larsen C, esteiras de ondas e vórtices oceânicos. Com base no exposto, podemos afirmar que, estatisticamente, nossos resultados são robustos e significantes, com nível de confiança entre 95% e 99%. A equação que propomos é um primeiro passo para inferir-se a espessura do gelo marinho a partir de coeficientes de retroespalhamento SAR. / The almost inaccessibility of large parts of the Southern Ocean makes the knowledge on the sea-ice thickness limited. This information is essential for determining the mass balance of this cryosphere component. During the transition from winter to spring 2006, a researcher team performed several sea-ice thickness profiles in the north and northwest of the Weddell Sea. They were aboard the German research vessel Polarstern from the Alfred-Wegener-Institute, Helmholtz-Zentrum für Polar- und Meeresforschung (AWI) and utilizing for those measurements a Helicopter-borne Electromagnetic system (HEM). The HEM is composed of an electromagnetic sensor, laser altimeter, and a Differential Global Positioning System (DGPS). Simultaneously, the European Space Agency (ESA) ENVISAT ASAR WSM acquired images from the research area. The main objective of this thesis is to determine the relationship between the Synthetic Aperture Radar backscattering and the HEM’s sea-ice thickness measurements using data listed above. We used the computer programs: NEST®, MATLAB®, EXCEL®, ArcGIS® / ArcMap®, to process, analyze and select images, in order to determine the relationship between the backscattering and the quasi-time-coincident ice thickness measurements. We projected the flights trajectories on the images obtained on the same dates, extracting the pixels pertaining to the places where the sea-ice thicknesses were determined. We ascribed ice thickness values for each area covered by the pixel. A statistical analysis determined that the best ice thickness parameter within a pixel area is its mean. Linear regression is the best way to adjust the relationship between the pixel backscatter value and the sea-ice thickness contained within the pixel area. The Pearson linear correlation coefficient, resulting from parametric analysis, indicates a strong correlation (0.75) between backscatter and sea-ice thickness. However, the nonparametric Spearman analysis resulted in a low correlation coefficient (0.06), which may indicate that the analyzed data consist of two distinct populations (e.g., first-year and multi-year ice). However, the Kolmogorov-Smirnov nonparametric analysis brought up the possibility that we just have not sampled the entire population. This could explain the no existence of sea-ice thicknesses values on the interval from -9.21dB to -1.35dB, which could lead to the low correlation coefficient in the Spearman analysis. We applied the linear equation: =,+, to the images, separating them into twelve classes: one for water and eleven to sea-ice, the latter in one-meter thickness intervals. The results enabled the observation of important oceanographic features such as open water channels, pressure ridges, sea-ice decay, icebergs motion (as indicators of currents, tides and winds), glaciers discharge, iceberg calving from Larsen C ice shelf, wakes and oceanic eddies. Based on the foregoing, we can say that, statistically, our results are robust and significantly, with a confidence level from 95% to 99%. The proposed equation is a first step to inferring sea-ice thickness from SAR backscatter coefficients.

A quase inacessibilidade de grandes partes do Oceano Austral torna o conhecimento da espessura do gelo marinho limitado. Esta informação é essencial para a determinação do balanço de massa deste componente da criosfera. Na transição do inverno para a primavera de 2006, uma equipe de pesquisadores, coletou uma série de perfis de espessura de gelo marinho no norte e noroeste do mar de Weddell. Eles estavam a bordo do navio de pesquisa alemão Polarstern do Alfred-Wegener-Institut, Helmholtz-Zentrum für Polar- und Meeresforschung (AWI) e utilizaram um sistema composto por um sensor eletromagnético, altímetro laser e um Sistema de Posicionamento Global Diferencial (DGPS) aerotransportado por helicóptero. Simultaneamente, a Agência Espacial Europeia (ESA), adquiriu imagens ENVISAT ASAR WSM da área de pesquisa. O objetivo principal desta tese é determinar a relação existente entre o retroespalhamento medido pelo Radar de Abertura Sintética e a espessura do gelo marinho obtida pelo HEM (Helicopter-borne ElectroMagnetic system) usando os dados citados acima. Utilizamos os programas de computador NEST®, MATLAB®, EXCEL®, ArcGIS®/ArcMAP®, para processar, analisar e selecionar as imagens, para determinar a relação entre o retroespalhamento e as medidas, quase-tempo-coincidentes, de espessura do gelo. Projetamos as trajetórias dos voos sobre as imagens obtidas nas mesmas datas e extraímos os pixels referentes aos locais onde foram medidas as espessuras de gelo marinho. Apropriamos os valores de espessura do gelo para a área de cada pixel sobre o qual se referiam. Uma análise estatística determinou que o parâmetro que melhor representa a espessura do gelo dentro da área do pixel é a média. A regressão linear é a melhor forma de ajuste das relações entre o valor de retroespalhamento do pixel e a espessura do gelo marinho contido na área deste pixel. O coeficiente de correlação linear de Pearson, resultante de análise paramétrica, indica uma forte correlação (0,75) entre retroespalhamento e espessura do gelo marinho. Porém, a análise não paramétrica de Spearman resultou em um coeficiente de correlação baixo (0,06) o que pode indicar que os dados analisados são compostos por duas populações distintas (e.g., gelo de primeiro ano e plurianual). Contudo, a análise não paramétrica de Kolmogorov-Smirnov aventa a possibilidade de que não tenhamos amostrado toda população. Esta seria a razão de não haver valores de espessuras relativas ao intervalo entre -9.21dB e -1.35dB, o que poderia induzir ao baixo valor do coeficiente de correlação na análise de Spearman. Aplicamos a equação linear: y=0,6345x+12,015 às imagens e pudemos separá-las em doze classes: uma para água e onze para gelo marinho, estas com intervalos de um metro. Isto possibilitou a observação de importantes feições oceanográficas como: canais de águas abertas; cristas de compressão, decaimento do gelo marinho, deslocamento de icebergs (como indicadores de correntes, marés e ventos), liberação de gelo por geleiras de maré, desprendimento de iceberg da plataforma de gelo Larsen C, esteiras de ondas e vórtices oceânicos. Com base no exposto, podemos afirmar que, estatisticamente, nossos resultados são robustos e significantes, com nível de confiança entre 95% e 99%. A equação que propomos é um primeiro passo para inferir-se a espessura do gelo marinho a partir de coeficientes de retroespalhamento SAR. / The almost inaccessibility of large parts of the Southern Ocean makes the knowledge on the sea-ice thickness limited. This information is essential for determining the mass balance of this cryosphere component. During the transition from winter to spring 2006, a researcher team performed several sea-ice thickness profiles in the north and northwest of the Weddell Sea. They were aboard the German research vessel Polarstern from the Alfred-Wegener-Institute, Helmholtz-Zentrum für Polar- und Meeresforschung (AWI) and utilizing for those measurements a Helicopter-borne Electromagnetic system (HEM). The HEM is composed of an electromagnetic sensor, laser altimeter, and a Differential Global Positioning System (DGPS). Simultaneously, the European Space Agency (ESA) ENVISAT ASAR WSM acquired images from the research area. The main objective of this thesis is to determine the relationship between the Synthetic Aperture Radar backscattering and the HEM’s sea-ice thickness measurements using data listed above. We used the computer programs: NEST®, MATLAB®, EXCEL®, ArcGIS® / ArcMap®, to process, analyze and select images, in order to determine the relationship between the backscattering and the quasi-time-coincident ice thickness measurements. We projected the flights trajectories on the images obtained on the same dates, extracting the pixels pertaining to the places where the sea-ice thicknesses were determined. We ascribed ice thickness values for each area covered by the pixel. A statistical analysis determined that the best ice thickness parameter within a pixel area is its mean. Linear regression is the best way to adjust the relationship between the pixel backscatter value and the sea-ice thickness contained within the pixel area. The Pearson linear correlation coefficient, resulting from parametric analysis, indicates a strong correlation (0.75) between backscatter and sea-ice thickness. However, the nonparametric Spearman analysis resulted in a low correlation coefficient (0.06), which may indicate that the analyzed data consist of two distinct populations (e.g., first-year and multi-year ice). However, the Kolmogorov-Smirnov nonparametric analysis brought up the possibility that we just have not sampled the entire population. This could explain the no existence of sea-ice thicknesses values on the interval from -9.21dB to -1.35dB, which could lead to the low correlation coefficient in the Spearman analysis. We applied the linear equation: =,+, to the images, separating them into twelve classes: one for water and eleven to sea-ice, the latter in one-meter thickness intervals. The results enabled the observation of important oceanographic features such as open water channels, pressure ridges, sea-ice decay, icebergs motion (as indicators of currents, tides and winds), glaciers discharge, iceberg calving from Larsen C ice shelf, wakes and oceanic eddies. Based on the foregoing, we can say that, statistically, our results are robust and significantly, with a confidence level from 95% to 99%. The proposed equation is a first step to inferring sea-ice thickness from SAR backscatter coefficients.

A quase inacessibilidade de grandes partes do Oceano Austral torna o conhecimento da espessura do gelo marinho limitado. Esta informação é essencial para a determinação do balanço de massa deste componente da criosfera. Na transição do inverno para a primavera de 2006, uma equipe de pesquisadores, coletou uma série de perfis de espessura de gelo marinho no norte e noroeste do mar de Weddell. Eles estavam a bordo do navio de pesquisa alemão Polarstern do Alfred-Wegener-Institut, Helmholtz-Zentrum für Polar- und Meeresforschung (AWI) e utilizaram um sistema composto por um sensor eletromagnético, altímetro laser e um Sistema de Posicionamento Global Diferencial (DGPS) aerotransportado por helicóptero. Simultaneamente, a Agência Espacial Europeia (ESA), adquiriu imagens ENVISAT ASAR WSM da área de pesquisa. O objetivo principal desta tese é determinar a relação existente entre o retroespalhamento medido pelo Radar de Abertura Sintética e a espessura do gelo marinho obtida pelo HEM (Helicopter-borne ElectroMagnetic system) usando os dados citados acima. Utilizamos os programas de computador NEST®, MATLAB®, EXCEL®, ArcGIS®/ArcMAP®, para processar, analisar e selecionar as imagens, para determinar a relação entre o retroespalhamento e as medidas, quase-tempo-coincidentes, de espessura do gelo. Projetamos as trajetórias dos voos sobre as imagens obtidas nas mesmas datas e extraímos os pixels referentes aos locais onde foram medidas as espessuras de gelo marinho. Apropriamos os valores de espessura do gelo para a área de cada pixel sobre o qual se referiam. Uma análise estatística determinou que o parâmetro que melhor representa a espessura do gelo dentro da área do pixel é a média. A regressão linear é a melhor forma de ajuste das relações entre o valor de retroespalhamento do pixel e a espessura do gelo marinho contido na área deste pixel. O coeficiente de correlação linear de Pearson, resultante de análise paramétrica, indica uma forte correlação (0,75) entre retroespalhamento e espessura do gelo marinho. Porém, a análise não paramétrica de Spearman resultou em um coeficiente de correlação baixo (0,06) o que pode indicar que os dados analisados são compostos por duas populações distintas (e.g., gelo de primeiro ano e plurianual). Contudo, a análise não paramétrica de Kolmogorov-Smirnov aventa a possibilidade de que não tenhamos amostrado toda população. Esta seria a razão de não haver valores de espessuras relativas ao intervalo entre -9.21dB e -1.35dB, o que poderia induzir ao baixo valor do coeficiente de correlação na análise de Spearman. Aplicamos a equação linear: y=0,6345x+12,015 às imagens e pudemos separá-las em doze classes: uma para água e onze para gelo marinho, estas com intervalos de um metro. Isto possibilitou a observação de importantes feições oceanográficas como: canais de águas abertas; cristas de compressão, decaimento do gelo marinho, deslocamento de icebergs (como indicadores de correntes, marés e ventos), liberação de gelo por geleiras de maré, desprendimento de iceberg da plataforma de gelo Larsen C, esteiras de ondas e vórtices oceânicos. Com base no exposto, podemos afirmar que, estatisticamente, nossos resultados são robustos e significantes, com nível de confiança entre 95% e 99%. A equação que propomos é um primeiro passo para inferir-se a espessura do gelo marinho a partir de coeficientes de retroespalhamento SAR. / The almost inaccessibility of large parts of the Southern Ocean makes the knowledge on the sea-ice thickness limited. This information is essential for determining the mass balance of this cryosphere component. During the transition from winter to spring 2006, a researcher team performed several sea-ice thickness profiles in the north and northwest of the Weddell Sea. They were aboard the German research vessel Polarstern from the Alfred-Wegener-Institute, Helmholtz-Zentrum für Polar- und Meeresforschung (AWI) and utilizing for those measurements a Helicopter-borne Electromagnetic system (HEM). The HEM is composed of an electromagnetic sensor, laser altimeter, and a Differential Global Positioning System (DGPS). Simultaneously, the European Space Agency (ESA) ENVISAT ASAR WSM acquired images from the research area. The main objective of this thesis is to determine the relationship between the Synthetic Aperture Radar backscattering and the HEM’s sea-ice thickness measurements using data listed above. We used the computer programs: NEST®, MATLAB®, EXCEL®, ArcGIS® / ArcMap®, to process, analyze and select images, in order to determine the relationship between the backscattering and the quasi-time-coincident ice thickness measurements. We projected the flights trajectories on the images obtained on the same dates, extracting the pixels pertaining to the places where the sea-ice thicknesses were determined. We ascribed ice thickness values for each area covered by the pixel. A statistical analysis determined that the best ice thickness parameter within a pixel area is its mean. Linear regression is the best way to adjust the relationship between the pixel backscatter value and the sea-ice thickness contained within the pixel area. The Pearson linear correlation coefficient, resulting from parametric analysis, indicates a strong correlation (0.75) between backscatter and sea-ice thickness. However, the nonparametric Spearman analysis resulted in a low correlation coefficient (0.06), which may indicate that the analyzed data consist of two distinct populations (e.g., first-year and multi-year ice). However, the Kolmogorov-Smirnov nonparametric analysis brought up the possibility that we just have not sampled the entire population. This could explain the no existence of sea-ice thicknesses values on the interval from -9.21dB to -1.35dB, which could lead to the low correlation coefficient in the Spearman analysis. We applied the linear equation: =,+, to the images, separating them into twelve classes: one for water and eleven to sea-ice, the latter in one-meter thickness intervals. The results enabled the observation of important oceanographic features such as open water channels, pressure ridges, sea-ice decay, icebergs motion (as indicators of currents, tides and winds), glaciers discharge, iceberg calving from Larsen C ice shelf, wakes and oceanic eddies. Based on the foregoing, we can say that, statistically, our results are robust and significantly, with a confidence level from 95% to 99%. The proposed equation is a first step to inferring sea-ice thickness from SAR backscatter coefficients.

Le bruit de chatoiement (speckle) lié aux systèmes d'imagerie cohérente a des conséquences sur l'analyse et l'interprétation des images radar à synthèse d'ouverture (RSO). Pour corriger ce défaut, nous profitons de séries temporelles d'images RSO bien recalées. Nous améliorons le filtre adaptatif temporel non-local à l'aide de méthodes performantes de débruitage adaptatif et proposons un filtrage temporel adaptatif basé sur les patchs. Pour réduire le biais du débruitage, nous proposons une méthode originale, rapide et efficace de débruitage multitemporel. L'idée principale de l'approche proposée est d'utiliser l'image dite "de ratio", donnée par le rapport entre l'image et la moyenne temporelle de la pile. Cette image de ratio est plus facile à débruiter qu'une image isolée en raison de sa meilleure stationnarité. Par ailleurs, les structures fines stables dans le temps sont bien préservées grâce au moyennage multitemporel. Disposant d'images débruitées, nous proposons ensuite d'utiliser la méthode du rapport de vraisemblance généralisé simplifié pour détecter les zones de changement ainsi que l'amplitude des changements et les instants de changements intéressants dans de longues séries d'images correctement recalées. En utilisant le partitionnement spectral, on applique le rapport de vraisemblance généralisé simplifié pour caractériser les changements des séries temporelles. Nous visualisons les résultats de détection en utilisant l'échelle de couleur 'jet' et une colorisation HSV. Ces méthodes ont été appliquées avec succès pour étudier des zones cultivées, des zones urbaines, des régions portuaires et des changements dus à des inondations. / The inherent speckle which is attached to any coherent imaging system affects the analysis and interpretation of synthetic aperture radar (SAR) images. To take advantage of well-registered multi-temporal SAR images, we improve the adaptive nonlocal temporal filter with state-of-the-art adaptive denoising methods and propose a patch based adaptive temporal filter. To address the bias problem of the denoising results, we propose a fast and efficient multitemporal despeckling method. The key idea of the proposed approach is the use of the ratio image, provided by the ratio between an image and the temporal mean of the stack. This ratio image is easier to denoise than a single image thanks to its improved stationarity. Besides, temporally stable thin structures are well-preserved thanks to the multi-temporal mean. Without reference image, we propose to use a patch-based auto-covariance residual evaluation method to examine the residual image and look for possible remaining structural contents. With speckle reduction images, we propose to use simplified generalized likelihood ratio method to detect the change area, change magnitude and change times in long series of well-registered images. Based on spectral clustering, we apply the simplified generalized likelihood ratio to detect the time series change types. Then, jet colormap and HSV colorization may be used to vividly visualize the detection results. These methods have been successfully applied to monitor farmland area, urban area, harbor region, and flooding area changes.

Microwave imaging finds numerous applications involving optically obscured targets. One particular area is breast cancer detection, since microwave technology promises fast low-cost image reconstruction without the use of harmful radiation typical of X-ray mammography. However, the success of microwave imaging is hindered by a critical issue, the complex nature of near-field electromagnetic scattering in tissue. To overcome this, specialized image reconstruction algorithms alongside sensitive measurement hardware are required. In this work, real-time near-field microwave imaging algorithms known as quantitative microwave holography and scattered power mapping are explored. They are experimentally demonstrated to identify potential tumor regions in tissue phantoms. Alongside this development, quality control techniques for evaluating microwave hardware are also described. Two new methods for improving the image reconstruction quality are also presented. First, a novel technique, which combines two commonly used mathematical approximations of scattering (the Born and Rytov approximations), is demonstrated yielding improved image reconstructions due to the complimentary nature of the approximations. Second, a range migration algorithm is introduced which enables near-field refocusing of a point-spread function (PSF), which is critical for algorithms that rely on measured PSFs to perform image reconstruction. / Thesis / Doctor of Philosophy (PhD) / Breast cancer remains as one of the highest causes of cancer-related deaths in women in Canada. Though X-ray mammography remains the gold standard for regular breast cancer screening, its use of harmful radiation, painful breast compression, and radiologist dependent evaluation remain as detracting factors for its use. Over the past 40 years, researchers have been exploring the use of microwave technology in place of X-ray mammography. Microwave radiation, used at power levels similar to that of a cellphone, has been demonstrated successfully in simulations of breast scans. However, in experimental evaluations with breast phantoms, the complex scattering path of the radiation through tissue complicates image reconstruction. In this thesis, methods of improving the accuracy of microwave algorithms are explored, alongside new breast phantom structures that replicate well the electrical properties of tissue. The results of this work demonstrate the flexibility of microwave imaging, and the adversities that still need to be overcome for it to begin seeing clinical use.

L'urbanisation et la gestion de l'environnement urbain et sa périphérie deviennent l'un des problèmes les plus cruciaux dans les pays développés et en développement. Dans ces circonstances, les données de télédétection sont une source importante d'information qui reflète les interactions entre les êtres humains et leur environnement. Compte tenu de leur indépendance totale des contraintes logistiques sur le terrain, l'éclairage (lumière du jour) et météorologiques (nuages) conditions, Synthetic Aperture Radar (SAR) les systèmes de satellites peuvent fournir des contributions importantes dans des environnements complexes de reconstruction 3-D. La nouvelle génération de haute résolution SAR capteurs comme COSMOSkyMed, TerraSAR-X RADARSAT-2 a permis d'acquérir des images SAR à haute résolution. Ici, l'attention est mis sur la technique pour l'imagerie 3-D nominée tomographie SAR: à partir d'une pile d'images ont été recueillies en utilisant les données multibaseline effectuées dans la configuration interférométrique, une telle technique permet d'extraire les informations de hauteur formant une ouverture synthétique dans la direction d'élévation afin d'obtenir une résolution sensiblement améliorée. Cette thèse de doctorat se concentre sur les potentialités élevées de techniques tomographiques en 3-D surveillance des changements et la caractérisation des zones complexes et denses bâties en utilisant des estimateurs mono-dimensionnelle de base comme Beamforming, Capon et MUSIC combinée au satellite très haute résolution des images SAR. 2-D et de l'analyse 3-D ont été présentés sur la zone urbaine de Paris en utilisant les données TerraSAR-X à haute résolution et de polarisation unique. Être porté principalement sur les techniques tomographiques 3-D, dans les méthodes de travail 4-D présentés, tels que le Compressive Sensing, ne sont pas pris en compte. Dans un premier temps, l'analyse de la qualité interférométrique de l'ensemble de données transformées a montré de bonnes valeurs de cohérence moyenne et ont permis de détecter des images considérées comme des valeurs aberrantes. L'extraction des tomographies 2-D sur l'azimut différent de profil a montré la capacité de distinguer plus d'un diffuseur à l'intérieur de la même cellule de résolution et de reconstituer les profilés de construction verticaux. Successivement, une caractérisation 3-D globale en terme de bâtiments hauteurs et réflectivité verticale a été réalisée dans le but de développer un outil de suivi des changements des structures simples. En outre, la possibilité de corriger les distorsions géométriques en raison de l'escale (qui affecte fortement ce genre de scénarios) et de déterminer les informations sur le nombre de diffuseurs (jusqu'à trois) et la réflectivité correspondant à l'intérieur d'une cellule de résolution ont été évalués. / The urbanization and the management of urban environment and its periphery become one of the most crucial issues in developed and developing countries. In these circumstances, remote sensing data are an important source of information that reflects interactions between human beings and their environment. Given their complete independence from logistic constraints on the ground, illumination (daylight), and weather (clouds) conditions, Synthetic Aperture Radar (SAR) satellite systems may provide important contributions in complex environments 3-D reconstruction. The new generation of high resolution SAR sensors as COSMO-SkyMed, TerraSAR-X and RADARSAT-2 allowed to acquire high resolution SAR imagery. Here the attention is put on the 3-D imaging technique called SAR Tomography: starting from a stack of images collected using multibaseline data performed in interferometric configuration, such a technique allows to retrieve height information forming a synthetic aperture in the elevation direction in order to achieve a substantially improved resolution. The present PhD thesis is focused on the high potentialities of tomographic techniques in 3-D change monitoring and characterization for complex and dense built-up areas using basic mono-dimensional estimators as Beamforming, Capon and MUSIC combined to very high satellite SAR resolution imagery. 2-D and 3-D analysis have been presented over the urban area of Paris using TerraSAR-X data at high resolution and single polarisation. Being mainly focused on the 3-D tomographic techniques, in the presented work 4-D methods, such as compressive sensing (CS), have not been taken into account. At first, the analysis of the interferometric quality of the processed data set has been performed and results showed good mean coherence values within the entire stack. The extraction of 2-D tomograms over different azimuth-profile has showed the capabilities to distinguish more than one scatterer within the same resolution cell and to reconstruct the vertical building profiles. Successively, a global 3-D characterization both in term of buildings heights and vertical reflectivity has been performed in order to develop a monitoring tool for the changes of single structures. Moreover, the possibility to correct the geometric distortions due to the layover (that strongly affects such kind of scenarios) and to determine the information about the number of scatterers (up to three) and the corresponding reflectivity within one resolution cell have been evaluated. Moreover an innovative time stability analysis of the observed scene have been carried out in order to detect the stable and unstable scatterers. Globally, the investigations showed noisier and sparser point clouds for the Capon method, whereas better capabilities for the Beamforming and MUSIC ones. Indeed, it was possible to detect different scatterers located within the same resolution cell and to resolve pixels affected by the layover. This has lead to perform a good reconstruction of building shape and location and a good estimation of their elevation. The 3-D time stability analysis demonstrated the possibility to monitor the 3-D change depending on the time. Eventually, it is possible to assert that processing high resolution SAR data allows to achieve a strong improvement in 3-D imaging capabilities. It has been demonstrated the potentialities of TomoSAR technique in distortions correction and in 3-D change monitoring using basic mono-dimensional estimators.

A time-domain backprojection processor for airborne synthetic aperture radar (SAR) has been developed at the University of Massachusetts’ Microwave Remote Sensing Lab (MIRSL). The aim of this work is to produce a SAR processor capable of addressing the motion compensation issues faced by frequency-domain processing algorithms, in order to create well focused SAR imagery suitable for interferometry. The time-domain backprojection algorithm inherently compensates for non-linear platform motion, dependent on the availability of accurate measurements of the motion. The implementation must manage the relatively high computational burden of the backprojection algorithm, which is done using modern graphics processing units (GPUs), programmed with NVIDIA’s CUDA language. An implementation of the Non-Equispaced Fast Fourier Transform (NERFFT) is used to enable efficient and accurate range interpolation as a critical step of the processing. The phase of time- domain processed imagery is dif erent than that of frequency-domain imagery, leading to a potentially different approach to interferometry. This general purpose SAR processor is designed to work with a novel, dual-frequency S- and Ka-band radar system developed at MIRSL as well as the UAVSAR instrument developed by NASA’s Jet Propulsion Laboratory. These instruments represent a wide range of SAR system parameters, ensuring the ability of the processor to work with most any airborne SAR. Results are presented from these two systems, showing good performance of the processor itself.

No description available.

Detecting water under vegetation is critical to tracking the status of geological ecosystems like wetlands. Researchers use different methods to estimate water presence, avoiding costly on-site measurements. Optical satellite imagery allows the automatic delineation of water using the concept of the Normalised Difference Water Index (NDWI). Still, optical imagery is subject to visibility conditions and cannot detect water under the vegetation, a typical situation for wetlands. Synthetic Aperture Radar (SAR) imagery works under all visibility conditions. It can detect water under vegetation but requires deep network algorithms to segment water presence, and manual annotation work is required to train the deep models. This project uses DEEPAQUA, a cross-modal knowledge distillation method, to eliminate the manual annotation needed to extract water presence from SAR imagery with deep neural networks. In this method, a deep student model (e.g., UNET) is trained to segment water in SAR imagery. The student model uses the NDWI algorithm as the non-parametric, cross-modal teacher. The key prerequisite is that NDWI works on the optical imagery taken from the exact location and simultaneously as the SAR. Three different deep architectures are tested in this project: UNET, SegNet, and UNET++, and the Otsu method is used as the baseline. Experiments on imagery from Swedish wetlands in 2020-2022 show that cross-modal distillation consistently achieved better segmentation performances across architectures than the baseline. Additionally, the UNET family of algorithms performed better than SegNet with a confidence of 95%. The UNET++ model achieved the highest Intersection Over Union (IOU) performance. However, no statistical evidence emerged that UNET++ performs better than UNET, with a confidence of 95%. In conclusion, this project shows that cross-modal knowledge distillation works well across architectures and removes tedious and expensive manual work hours when detecting water from SAR imagery. Further research could evaluate performances on other datasets and student architectures. / Att upptäcka vatten under vegetation är avgörande för att hålla koll på statusen på geologiska ekosystem som våtmarker. Forskare använder olika metoder för att uppskatta vattennärvaro vilket undviker kostsamma mätningar på plats. Optiska satellitbilder tillåter automatisk avgränsning av vatten med hjälp av konceptet Normalised Difference Water Index (NDWI). Optiska bilder fortfarande beroende av siktförhållanden och kan inte upptäcka vatten under vegetationen, en typisk situation för våtmarker. Synthetic Aperture Radar (SAR)-bilder fungerar under alla siktförhållanden. Den kan detektera vatten under vegetation men kräver djupa nätverksalgoritmer för att segmentera vattennärvaro, och manuellt anteckningsarbete krävs för att träna de djupa modellerna. Detta projekt använder DEEPAQUA, en cross-modal kunskapsdestillationsmetod, för att eliminera det manuella annoteringsarbete som behövs för att extrahera vattennärvaro från SAR-bilder med djupa neurala nätverk. I denna metod tränas en djup studentmodell (t.ex. UNET) att segmentera vatten i SAR-bilder semantiskt. Elevmodellen använder NDWI, som fungerar på de optiska bilderna tagna från den exakta platsen och samtidigt som SAR, som den icke-parametriska, cross-modal lärarmodellen. Tre olika djupa arkitekturer testas i detta examensarbete: UNET, SegNet och UNET++, och Otsu-metoden används som baslinje. Experiment på bilder tagna på svenska våtmarker 2020-2022 visar att cross-modal destillation konsekvent uppnådde bättre segmenteringsprestanda över olika arkitekturer jämfört med baslinjen. Dessutom presterade UNET-familjen av algoritmer bättre än SegNet med en konfidens på 95%. UNET++-modellen uppnådde högsta prestanda för Intersection Over Union (IOU). Det framkom dock inga statistiska bevis för att UNET++ presterar bättre än UNET, med en konfidens på 95%. Sammanfattningsvis visar detta projekt att cross-modal kunskapsdestillation fungerar bra över olika arkitekturer och tar bort tidskrävande och kostsamma manuella arbetstimmar vid detektering av vatten från SAR-bilder. Ytterligare forskning skulle kunna utvärdera prestanda på andra datamängder och studentarkitekturer.