Développements algorithmiques pour l’amélioration des résultats de l’interférométrie RADAR en milieu urbain

Le suivi des espaces urbanisés et de leurs dynamiques spatio-temporelles représente un enjeu important pour la population urbaine, autant sur le plan environnemental, économique et social. Avec le lancement des satellites portant des radars à synthèse d’ouverture de la nouvelle génération (TerraSAR-X, COSMO-SkyMed, ALOS, RADARSAT-2,Sentinel-1, Constellation RADARSAT), il est possible d’obtenir des séries temporelles d’images avec des résolutions spatiales et temporelles fines. Ces données multitemporelles aident à mieux analyser et décrire les structures urbaines et leurs variations dans l’espace et dans le temps. L’interférométrie par satellite est effectuée en comparant les phases des images RSO prises à différents passages du satellite au-dessus du même territoire. En optant pour des positions du satellite séparées d’une longue ligne de base, l’InSAR mène à la création des modèles numériques d’altitude (MNA). Si cette ligne de base est courte et à la limite nulle, nous avons le cas de l’interférométrie différentielle (DInSAR) qui mène à l’estimation du mouvement possible du terrain entre les deux acquisitions. Pour toutes les deux applications de l’InSAR, deux opérations sont importantes qui garantissent la génération des interférogrammes de qualité. La première est le filtrage du bruit omniprésent dans les phases interférométriques et la deuxième est le déroulement des phases. Ces deux opérations deviennent particulièrement complexes en milieu urbain où au bruit des phases s’ajoutent des fréquents sauts et discontinuités des phases dus à la présence des bâtiments et d’autres structures surélevées. L’objectif de cette recherche est le développement des nouveaux algorithmes de filtrage et de déroulement de phase qui puissent mieux performer que les algorithmes considérés comme référence dans ce domaine. Le but est d’arriver à générer des produits InSAR de qualité en milieu urbain. Concernant le filtrage, nous avons établi un algorithme qui est une nouvelle formulation du filtre Gaussien anisotrope adaptatif. Quant à l’algorithme de déroulement de phase, il est fondé sur la minimisation de l’énergie par un algorithme génétique ayant recours à une modélisation contextuelle du champ de phase. Différents tests ont été effectués avec des images RSO simulées et réelles qui démontrent le potentiel de nos algorithmes qui dépasse à maints égards celui des algorithmes standard. Enfin, pour atteindre le but de notre recherche, nous avons intégré nos algorithmes dans l’environnement du logiciel SNAP et appliqué l’ensemble de la procédure pour générer un MNA avec des images RADARSAT-2 de haute résolution d’un secteur de la Ville de Montréal (Canada) ainsi que des cartes des mouvements du terrain dans la région de la Ville de Mexico (Mexique) avec des images de Sentinel-1 de résolution plutôt moyenne. La comparaison des résultats obtenus avec des données provenant des sources externes de qualité a aussi démontré le fort potentiel de nos algorithmes.The monitoring of urban areas and their spatiotemporal dynamics is an important issue for the urban population, at the environmental, economic, as well as social level. With the launch of satellites carrying next-generation synthetic aperture radars (TerraSAR-X, COSMO-SkyMed, ALOS, RADARSAT-2, Sentinel-1, Constellation RADARSAT), it is possible to obtain time series of images with fine temporal and spatial resolutions. These multitemporal data help to better analyze and describe urban structures, and their variations in space and time. Satellite interferometry is performed by comparing the phases of SAR images taken at different satellite passes over the same territory. By opt-ing for satellite positions separated by a long baseline, InSAR leads to the creation of digital elevation models (DEM). If this baseline is short and, at the limit zero, we have the case of differential interferometry (DInSAR) which leads to the estimation of the possible movement of the land between the two acquisitions. In both InSAR applica-tions, two operations are important that ensure the generation of quality interferograms. The first is the filtering of ubiquitous noise in the interferometric phases and the second is the unwrapping of the phases. These two operations become particularly complex in urban areas where the phase noise is added to the frequent jumps and discontinuities of phases due to the presence of buildings and other raised structures. The objective of this research is the development of new filtering and phase unwrap-ping algorithms that can perform better than algorithms considered as reference in this field. The goal is to generate quality InSAR products in urban areas. Regarding filtering, we have established an algorithm that is a new formulation of the adaptive anisotropic Gaussian filter. As for the phase unwrapping algorithm, it is based on the minimization of energy by a genetic algorithm using contextual modelling of the phase field. Various tests have been carried out with simulated and real SAR images that demonstrated the potential of our algorithms that in many respects exceeds that of standard algorithms. Finally, to achieve the goal of our research, we integrated our algorithms into the SNAP software environment and applied the entire procedure to generate a DEM with high-resolution RADARSAT-2 images from an area of the City of Montreal (Canada) as well as maps of land movement in the Mexico City region (Mexico) with relatively medium-resolution Sentinel-1 images. Comparison of the results with data from external quality sources also demonstrated the strong potential of our algorithms

A Sparsity-Based InSAR Phase Denoising Algorithm Using Nonlocal Wavelet Shrinkage

An interferometric synthetic aperture radar (InSAR) phase denoising algorithm using the local sparsity of wavelet coefficients and nonlocal similarity of grouped blocks was developed. From the Bayesian perspective, the double-l1 norm regularization model that enforces the local and nonlocal sparsity constraints was used. Taking advantages of coefficients of the nonlocal similarity between group blocks for the wavelet shrinkage, the proposed algorithm effectively filtered the phase noise. Applying the method to simulated and acquired InSAR data, we obtained satisfactory results. In comparison, the algorithm outperformed several widely-used InSAR phase denoising approaches in terms of the number of residues, root-mean-square errors and other edge preservation indexes

Intercomparison and Validation of SAR-Based Ice Velocity Measurement Techniques within the Greenland Ice Sheet CCI Project

Ice velocity is one of the products associated with the Ice Sheets Essential Climate Variable. This paper describes the intercomparison and validation of ice-velocity measurements carried out by several international research groups within the European Space Agency Greenland Ice Sheet Climate Change Initiative project, based on space-borne Synthetic Aperture Radar (SAR) data. The goal of this activity was to survey the best SAR-based measurement and error characterization approaches currently in practice. To this end, four experiments were carried out, related to different processing techniques and scenarios, namely differential SAR interferometry, multi aperture SAR interferometry and offset-tracking of incoherent as well as of partially-coherent data. For each task, participants were provided with common datasets covering areas located on the Greenland ice-sheet margin and asked to provide mean velocity maps, quality characterization and a description of processing algorithms and parameters. The results were then intercompared and validated against GPS data, revealing in several cases significant differences in terms of coverage and accuracy. The algorithmic steps and parameters influencing the coverage, accuracy and spatial resolution of the measurements are discussed in detail for each technique, as well as the consistency between quality parameters and validation results. This allows several recommendations to be formulated, in particular concerning procedures which can reduce the impact of analyst decisions, and which are often found to be the cause of sub-optimal algorithm performance

Atmospheric artifacts correction for InSAR using empirical model and numerical weather prediction models

lnSAR has been proved its unprecedented ability and merits of monitoring ground deformation on large scale with centimeter to millimeter scale accuracy. However, several factors affect the reliability and accuracy of its applications. Among them, atmospheric artifacts due to spatial and temporal variations of atmosphere state often pose noise to interferograms. Therefore, atmospheric artifacts m itigalion remains one of the biggest challenges to be addressed in the In SAR community. State-of-the-art research works have revealed atmospheric artifacts can be partially compensated with empirical models, temporal-spatial filtering approach in lnSAR time series, pointwise GPS zenith path delay and numerical weather prediction models. In this thesis, firstly, we further develop a covariance weighted linear empirical model correction method. Secondly, a realistic LOS direction integration approach based on global reanalysis data is employed and comprehensively compared with the conventional method that integrates along zenith direction. Finally, the realistic integration method is applied to local WRF numerical forecast model data. l'vbreover, detailed comparisons between different global reanalysis data and local WRF model are assessed. In terms of empirical models correcting methods, many publications have studied correcting stratified tropospheric phase delay by assuming a linear model between them and topography. However, most of these studies ha\19 not considered the effect of turbulent atmospheric artefacts when adjusting the linear model to data. In this thesis, an improved technique that minimizes the influence of turbulent atmosphere in the model adjustment has been presented. In the proposed algorithm, the model is adjusted to the phase differences of pixels instead of using the unwrapped phase of each pixel. In addition, the different phase differences are weighted as a function of its APS covariance estimated from an empirical variogram to reduce in the model adjustment the impact of pixel pairs with significant turbulent atmosphere. The performance of the proposed method has been validated with both simulated and real Sentinel-1 SAR data in Tenerife island, Spain. Considering methods using meteorological observations to mitigate APS, an accurate realistic com puling strategy utilizing global atmospheric reanalysis data has been implemented. With the approach, the realistic LOS path along satellite and the monitored points is considered, rather than converting from zenith path delay. Com pared with zenith delay based method, the biggest advantage is that it can avoid errors caused by anisotropic atmospheric behaviour. The accurate integration method is validated with Sentinel-1 data in three test sites: Tenerife island, Spain, Almeria, Spain and Crete island, Greece. Compared to conventional zenith method, the realistic integration method shows great improvement. A variety of global reanalysis data are available from different weather forecasting organizations, such as ERA-Interim, ERAS, MERRA2. In this study, the realistic integration mitigation method is assessed on these different reanalysis data. The results show that these data are feasible to mitigate APS to some extent in most cases. The assessment also demonstrates that the ERAS performs the best statistically, compared to other global reanalysis data. l'vbreover, as local numerical weather forecast models have the ability to predict high spatial resolution atmospheric parameters, by using which, it has the potential to achieve APS mitigation. In this thesis, the realistic integration method is also employed on the local WRF model data in Tenerife and Almeria test s ites. However, it turns out that the WRF model performs worse than the original global reanalysis data.Las técnicas lnSAR han demostrado su capacidad sin precedentes y méritos para el monitoreo de la deformaci6n del suelo a gran escala con una precisión centimétrica o incluso milimétrica. Sin embargo, varios factores afectan la fiabilidad y precisión de sus aplicaciones. Entre ellos, los artefactos atmosféricos debidos a variaciones espaciales y temporales del estado de la atm6sfera a menudo añaden ruido a los interferogramas. Por lo tanto, la mitigación de los artefactos atmosféricos sigue siendo uno de los mayores desafíos a abordar en la comunidad lnSAR. Los trabajos de investigaci6n de vanguardia han revelado que los artefactos atmosféricos se pueden compensar parcialmente con modelos empíricos, enfoque de filtrado temporal-espacial en series temporales lnSAR, retardo puntual del camino cenital con GPS y modelos numéricos de predicción meteorológica. En esta tesis, en primer lugar, desarrollamos un método de corrección de modelo empírico lineal ponderado por covarianza. En segundo lugar, se emplea un enfoque realista de integracion de dirección LOS basado en datos de reanálisis global y se compara exhaustivamente con el método convencional que se integra a lo largo de la dirección cenital. Finalmente, el método de integraci6n realista se aplica a los datos del modelo de pronóstico numérico WRF local. Ademas, se evalúan las comparaciones detalladas entre diferentes datos de reanálisis global y el modelo WRF local. En términos de métodos de corrección con modelos empíricos, muchas publicaciones han estudiado la corrección del retraso estratificado de la fase troposférica asumiendo un modelo lineal entre ellos y la topografía. Sin embargo, la mayoría de estos estudios no han considerado el efecto de los artefactos atmosféricos turbulentos al ajustar el modelo lineal a los datos. En esta tesis, se ha presentado una técnica mejorada que minimiza la influencia de la atm6sfera turbulenta en el ajuste del modelo. En el algoritmo propuesto, el modelo se ajusta a las diferencias de fase de los pixeles en lugar de utilizar la fase sin desenrollar de cada pixel. Además, las diferentes diferencias de fase se ponderan en función de su covarianza APS estimada a partir de un variograma empírico para reducir en el ajuste del modelo el impacto de los pares de pixeles con una atm6sfera turbulenta significativa. El rendimiento del método propuesto ha sido validado con datos SAR Sentinel-1 simulados y reales en la isla de Tenerife, España. Teniendo en cuenta los métodos que utilizan observaciones meteorológicas para mitigar APS, se ha implementado una estrategia de computación realista y precisa que utiliza datos de reanálisis atmosférico global. Con el enfoque, se considera el camino realista de LOS a lo largo del satélite y los puntos monitoreados, en lugar de convertirlos desde el retardo de la ruta cenital. En comparación con el método basado en la demora cenital, la mayor ventaja es que puede evitar errores causados por el comportamiento atmosférico anisotrópico. El método de integración preciso se valida con los datos de Sentinel-1 en tres sitios de prueba: la isla de Tenerife, España, Almería, España y la isla de Creta, Grecia. En comparación con el método cenital convencional, el método de integración realista muestra una gran mejora

Atmospheric artifacts correction for InSAR using empirical model and numerical weather prediction models

lnSAR has been proved its unprecedented ability and merits of monitoring ground deformation on large scale with centimeter to millimeter scale accuracy. However, several factors affect the reliability and accuracy of its applications. Among them, atmospheric artifacts due to spatial and temporal variations of atmosphere state often pose noise to interferograms. Therefore, atmospheric artifacts m itigalion remains one of the biggest challenges to be addressed in the In SAR community. State-of-the-art research works have revealed atmospheric artifacts can be partially compensated with empirical models, temporal-spatial filtering approach in lnSAR time series, pointwise GPS zenith path delay and numerical weather prediction models. In this thesis, firstly, we further develop a covariance weighted linear empirical model correction method. Secondly, a realistic LOS direction integration approach based on global reanalysis data is employed and comprehensively compared with the conventional method that integrates along zenith direction. Finally, the realistic integration method is applied to local WRF numerical forecast model data. l'vbreover, detailed comparisons between different global reanalysis data and local WRF model are assessed. In terms of empirical models correcting methods, many publications have studied correcting stratified tropospheric phase delay by assuming a linear model between them and topography. However, most of these studies ha\19 not considered the effect of turbulent atmospheric artefacts when adjusting the linear model to data. In this thesis, an improved technique that minimizes the influence of turbulent atmosphere in the model adjustment has been presented. In the proposed algorithm, the model is adjusted to the phase differences of pixels instead of using the unwrapped phase of each pixel. In addition, the different phase differences are weighted as a function of its APS covariance estimated from an empirical variogram to reduce in the model adjustment the impact of pixel pairs with significant turbulent atmosphere. The performance of the proposed method has been validated with both simulated and real Sentinel-1 SAR data in Tenerife island, Spain. Considering methods using meteorological observations to mitigate APS, an accurate realistic com puling strategy utilizing global atmospheric reanalysis data has been implemented. With the approach, the realistic LOS path along satellite and the monitored points is considered, rather than converting from zenith path delay. Com pared with zenith delay based method, the biggest advantage is that it can avoid errors caused by anisotropic atmospheric behaviour. The accurate integration method is validated with Sentinel-1 data in three test sites: Tenerife island, Spain, Almeria, Spain and Crete island, Greece. Compared to conventional zenith method, the realistic integration method shows great improvement. A variety of global reanalysis data are available from different weather forecasting organizations, such as ERA-Interim, ERAS, MERRA2. In this study, the realistic integration mitigation method is assessed on these different reanalysis data. The results show that these data are feasible to mitigate APS to some extent in most cases. The assessment also demonstrates that the ERAS performs the best statistically, compared to other global reanalysis data. l'vbreover, as local numerical weather forecast models have the ability to predict high spatial resolution atmospheric parameters, by using which, it has the potential to achieve APS mitigation. In this thesis, the realistic integration method is also employed on the local WRF model data in Tenerife and Almeria test s ites. However, it turns out that the WRF model performs worse than the original global reanalysis data.Las técnicas lnSAR han demostrado su capacidad sin precedentes y méritos para el monitoreo de la deformaci6n del suelo a gran escala con una precisión centimétrica o incluso milimétrica. Sin embargo, varios factores afectan la fiabilidad y precisión de sus aplicaciones. Entre ellos, los artefactos atmosféricos debidos a variaciones espaciales y temporales del estado de la atm6sfera a menudo añaden ruido a los interferogramas. Por lo tanto, la mitigación de los artefactos atmosféricos sigue siendo uno de los mayores desafíos a abordar en la comunidad lnSAR. Los trabajos de investigaci6n de vanguardia han revelado que los artefactos atmosféricos se pueden compensar parcialmente con modelos empíricos, enfoque de filtrado temporal-espacial en series temporales lnSAR, retardo puntual del camino cenital con GPS y modelos numéricos de predicción meteorológica. En esta tesis, en primer lugar, desarrollamos un método de corrección de modelo empírico lineal ponderado por covarianza. En segundo lugar, se emplea un enfoque realista de integracion de dirección LOS basado en datos de reanálisis global y se compara exhaustivamente con el método convencional que se integra a lo largo de la dirección cenital. Finalmente, el método de integraci6n realista se aplica a los datos del modelo de pronóstico numérico WRF local. Ademas, se evalúan las comparaciones detalladas entre diferentes datos de reanálisis global y el modelo WRF local. En términos de métodos de corrección con modelos empíricos, muchas publicaciones han estudiado la corrección del retraso estratificado de la fase troposférica asumiendo un modelo lineal entre ellos y la topografía. Sin embargo, la mayoría de estos estudios no han considerado el efecto de los artefactos atmosféricos turbulentos al ajustar el modelo lineal a los datos. En esta tesis, se ha presentado una técnica mejorada que minimiza la influencia de la atm6sfera turbulenta en el ajuste del modelo. En el algoritmo propuesto, el modelo se ajusta a las diferencias de fase de los pixeles en lugar de utilizar la fase sin desenrollar de cada pixel. Además, las diferentes diferencias de fase se ponderan en función de su covarianza APS estimada a partir de un variograma empírico para reducir en el ajuste del modelo el impacto de los pares de pixeles con una atm6sfera turbulenta significativa. El rendimiento del método propuesto ha sido validado con datos SAR Sentinel-1 simulados y reales en la isla de Tenerife, España. Teniendo en cuenta los métodos que utilizan observaciones meteorológicas para mitigar APS, se ha implementado una estrategia de computación realista y precisa que utiliza datos de reanálisis atmosférico global. Con el enfoque, se considera el camino realista de LOS a lo largo del satélite y los puntos monitoreados, en lugar de convertirlos desde el retardo de la ruta cenital. En comparación con el método basado en la demora cenital, la mayor ventaja es que puede evitar errores causados por el comportamiento atmosférico anisotrópico. El método de integración preciso se valida con los datos de Sentinel-1 en tres sitios de prueba: la isla de Tenerife, España, Almería, España y la isla de Creta, Grecia. En comparación con el método cenital convencional, el método de integración realista muestra una gran mejora.Postprint (published version

A Tutorial on Speckle Reduction in Synthetic Aperture Radar Images

Speckle is a granular disturbance, usually modeled as a multiplicative noise, that affects synthetic aperture radar (SAR) images, as well as all coherent images. Over the last three decades, several methods have been proposed for the reduction of speckle, or despeckling, in SAR images. Goal of this paper is making a comprehensive review of despeckling methods since their birth, over thirty years ago, highlighting trends and changing approaches over years. The concept of fully developed speckle is explained. Drawbacks of homomorphic filtering are pointed out. Assets of multiresolution despeckling, as opposite to spatial-domain despeckling, are highlighted. Also advantages of undecimated, or stationary, wavelet transforms over decimated ones are discussed. Bayesian estimators and probability density function (pdf) models in both spatial and multiresolution domains are reviewed. Scale-space varying pdf models, as opposite to scale varying models, are promoted. Promising methods following non-Bayesian approaches, like nonlocal (NL) filtering and total variation (TV) regularization, are reviewed and compared to spatial- and wavelet-domain Bayesian filters. Both established and new trends for assessment of despeckling are presented. A few experiments on simulated data and real COSMO-SkyMed SAR images highlight, on one side the costperformance tradeoff of the different methods, on the other side the effectiveness of solutions purposely designed for SAR heterogeneity and not fully developed speckle. Eventually, upcoming methods based on new concepts of signal processing, like compressive sensing, are foreseen as a new generation of despeckling, after spatial-domain and multiresolution-domain method

Inpainting of Cyclic Data using First and Second Order Differences

Cyclic data arise in various image and signal processing applications such as interferometric synthetic aperture radar, electroencephalogram data analysis, and color image restoration in HSV or LCh spaces. In this paper we introduce a variational inpainting model for cyclic data which utilizes our definition of absolute cyclic second order differences. Based on analytical expressions for the proximal mappings of these differences we propose a cyclic proximal point algorithm (CPPA) for minimizing the corresponding functional. We choose appropriate cycles to implement this algorithm in an efficient way. We further introduce a simple strategy to initialize the unknown inpainting region. Numerical results both for synthetic and real-world data demonstrate the performance of our algorithm.Comment: accepted Converence Paper at EMMCVPR'1

Advanced Geoscience Remote Sensing

Nowadays, advanced remote sensing technology plays tremendous roles to build a quantitative and comprehensive understanding of how the Earth system operates. The advanced remote sensing technology is also used widely to monitor and survey the natural disasters and man-made pollution. Besides, telecommunication is considered as precise advanced remote sensing technology tool. Indeed precise usages of remote sensing and telecommunication without a comprehensive understanding of mathematics and physics. This book has three parts (i) microwave remote sensing applications, (ii) nuclear, geophysics and telecommunication; and (iii) environment remote sensing investigations

Elevation and Deformation Extraction from TomoSAR

3D SAR tomography (TomoSAR) and 4D SAR differential tomography (Diff-TomoSAR) exploit multi-baseline SAR data stacks to provide an essential innovation of SAR Interferometry for many applications, sensing complex scenes with multiple scatterers mapped into the same SAR pixel cell. However, these are still influenced by DEM uncertainty, temporal decorrelation, orbital, tropospheric and ionospheric phase distortion and height blurring. In this thesis, these techniques are explored. As part of this exploration, the systematic procedures for DEM generation, DEM quality assessment, DEM quality improvement and DEM applications are first studied. Besides, this thesis focuses on the whole cycle of systematic methods for 3D & 4D TomoSAR imaging for height and deformation retrieval, from the problem formation phase, through the development of methods to testing on real SAR data. After DEM generation introduction from spaceborne bistatic InSAR (TanDEM-X) and airborne photogrammetry (Bluesky), a new DEM co-registration method with line feature validation (river network line, ridgeline, valley line, crater boundary feature and so on) is developed and demonstrated to assist the study of a wide area DEM data quality. This DEM co-registration method aligns two DEMs irrespective of the linear distortion model, which improves the quality of DEM vertical comparison accuracy significantly and is suitable and helpful for DEM quality assessment. A systematic TomoSAR algorithm and method have been established, tested, analysed and demonstrated for various applications (urban buildings, bridges, dams) to achieve better 3D & 4D tomographic SAR imaging results. These include applying Cosmo-Skymed X band single-polarisation data over the Zipingpu dam, Dujiangyan, Sichuan, China, to map topography; and using ALOS L band data in the San Francisco Bay region to map urban building and bridge. A new ionospheric correction method based on the tile method employing IGS TEC data, a split-spectrum and an ionospheric model via least squares are developed to correct ionospheric distortion to improve the accuracy of 3D & 4D tomographic SAR imaging. Meanwhile, a pixel by pixel orbit baseline estimation method is developed to address the research gaps of baseline estimation for 3D & 4D spaceborne SAR tomography imaging. Moreover, a SAR tomography imaging algorithm and a differential tomography four-dimensional SAR imaging algorithm based on compressive sensing, SAR interferometry phase (InSAR) calibration reference to DEM with DEM error correction, a new phase error calibration and compensation algorithm, based on PS, SVD, PGA, weighted least squares and minimum entropy, are developed to obtain accurate 3D & 4D tomographic SAR imaging results. The new baseline estimation method and consequent TomoSAR processing results showed that an accurate baseline estimation is essential to build up the TomoSAR model. After baseline estimation, phase calibration experiments (via FFT and Capon method) indicate that a phase calibration step is indispensable for TomoSAR imaging, which eventually influences the inversion results. A super-resolution reconstruction CS based study demonstrates X band data with the CS method does not fit for forest reconstruction but works for reconstruction of large civil engineering structures such as dams and urban buildings. Meanwhile, the L band data with FFT, Capon and the CS method are shown to work for the reconstruction of large manmade structures (such as bridges) and urban buildings

Feasibility of Imaging Tissue Electrical Conductivity by Switching Field Gradients with MRI.

Tissue conductivity is a biophysical marker of tissue structure and physiology. Present methods of measuring tissue conductivity are limited. Electrical impedance tomography, and magnetic resonance electrical impedance tomography rely on passing external current through the object being imaged, which prevents its use in most human imaging. Recently, the RF field used for MR excitation has been used to non-invasively measure tissue conductivity. This technique is promising, but conductivity at higher frequencies is less sensitive to tissue structure. Measuring tissue conductivity non-invasively at low frequencies remains elusive. It has been proposed that eddy currents generated during the rise and decay of gradient pulses could act as a current source to map low-frequency conductivity. This work centers on a gradient echo pulse sequence that uses large gradients prior to excitation to create eddy currents. The electric and magnetic fields during a gradient pulse are simulated by a finite-difference time-domain simulation. The sequence is also tested with a phantom and an animal MRI scanner equipped with gradients of high gradient strengths and slew rate. The simulation demonstrates that eddy currents in materials with conductivity similar to biological tissue decay with a half-life on the order of nanoseconds and any eddy currents generated prior to excitation decay completely before influencing the RF signal. Gradient-induced eddy currents can influence phase accumulation after excitation but the effect is too small to image. The animal scanner images show no measurable phase accumulation. Measuring low-frequency conductivity by gradient-induced eddy currents is presently unfeasible