Detection of anomalies in radio tomography of asteroids: Source count and forward errors

Abstract The purpose of this study was to advance numerical methods for radio tomography in which asteroid׳s internal electric permittivity distribution is to be recovered from radio frequency data gathered by an orbiter. The focus was on signal generation via multiple sources (transponders) providing one potential, or even essential, scenario to be implemented in a challenging in situ measurement environment and within tight payload limits. As a novel feature, the effects of forward errors including noise and a priori uncertainty of the forward (data) simulation were examined through a combination of the iterative alternating sequential (IAS) inverse algorithm and finite-difference time-domain (FDTD) simulation of time evolution data. Single and multiple source scenarios were compared in two-dimensional localization of permittivity anomalies. Three different anomaly strengths and four levels of total noise were tested. Results suggest, among other things, that multiple sources can be necessary to obtain appropriate results, for example, to distinguish three separate anomalies with permittivity less or equal than half of the background value, relevant in recovery of internal cavities.

[1]  Harri Hakula,et al.  Conditionally Gaussian Hypermodels for Cerebral Source Localization , 2008, SIAM J. Imaging Sci..

[2]  Mikko Kaasalainen,et al.  Interpretation of lightcurves of atmosphereless bodies. II - Practical aspects of inversion , 1992 .

[3]  L. Kettunen,et al.  Yee‐like schemes on a tetrahedral mesh, with diagonal lumping , 1999 .

[4]  Wlodek Kofman,et al.  Antennas for sounding of a cometary nucleus in the Rosetta mission , 2001 .

[5]  J. Idier.  Bayesian Approach to Inverse Problems: Idier/Bayesian , 2010 .

[6]  Anthony O'Hagan,et al.  Kendall's Advanced Theory of Statistics: Vol. 2B, Bayesian Inference. , 1996 .

[7]  A. Nuttall.  Some windows with very good sidelobe behavior , 1981 .

[8]  Slawomir Wiak,et al.  Optimization and Inverse Problems in Electromagnetism , 2003 .

[9]  D. Braess.  Finite Elements: Theory, Fast Solvers, and Applications in Solid Mechanics , 1995 .

[10]  Sampsa Pursiainen,et al.  Computational methods in electromagnetic biomedical inverse problems , 2008 .

[11]  A. Walther,et al.  Three-dimensional reconstruction of a comet nucleus by optimal control of Maxwell's equations: A contribution to the experiment CONSERT onboard space craft Rosetta , 2010, 2010 IEEE Radar Conference.

[12]  Clifford H. Thurber,et al.  Parameter estimation and inverse problems , 2005 .

[13]  Bert Reubens,et al.  Assessing and analyzing 3D architecture of woody root systems, a review of methods and applications in tree and soil stability, resource acquisition and allocation , 2008, Plant and Soil.

[14]  Kurt M. Strack,et al.  Society of Exploration Geophysicists , 2007 .

[15]  Stephan Ulamec,et al.  Radar Tomography of Asteroids , 2011 .

[16]  Eric L. Mokole,et al.  Ultra-Wideband, Short-Pulse Electromagnetics 7 , 2007 .

[17]  Stephan Ulamec,et al.  ASSERT : a radar Tomography of Asteroids , 2011 .

[18]  V. Utsi,et al.  Advances in long-range GPR systems and their applications to mineral exploration, geotechnical and static correction problems , 2009 .

[19]  A. P. Annan,et al.  Ground-penetrating radar for high-resolution mapping of soil and rock stratigraphy , 1989 .

[20]  Wlodek Kofman,et al.  Consert experiment: description and performances in view of the new target , 2004 .

[21]  Wlodek Kofman,et al.  A two dimensional simulation of the CONSERT experiment (radio tomography of comet Wirtanen) , 1999 .

[22]  Luigi Colangeli,et al.  The new rosetta targets : observations, simulations and instrument performances , 2004 .

[23]  Misac N. Nabighian.  Electromagnetic Methods in Applied Geophysics-Applications : part a and part b , 1991 .

[24]  Misac N. Nabighian,et al.  Electromagnetic Methods in Applied Geophysics , 1988 .

[25]  Mikko Kaasalainen,et al.  Multimodal inverse problems: Maximum compatibility estimateand shape reconstruction , 2011 .

[26]  E. Jakeman,et al.  Non-Gaussian models for the statistics of scattered waves , 1988 .

[27]  Jun S. Liu,et al.  Monte Carlo strategies in scientific computing , 2001 .

[28]  Wlodek Kofman,et al.  Dielectric properties of comet analog refractory materials , 2002 .

[29]  D. Daniels.  Ground Penetrating Radar , 2005 .

[30]  Martin Burger,et al.  Hierarchical Bayesian inference for the EEG inverse problem using realistic FE head models: Depth localization and source separation for focal primary currents , 2012, NeuroImage.

[31]  D. Plettemeier,et al.  The Comet Nucleus Sounding Experiment by Radiowave Transmission (CONSERT): A Short Description of the Instrument and of the Commissioning Stages , 2007 .

[32]  Mikko Kaasalainen,et al.  Inverse problems of generalized projection operators , 2006 .

[33]  T. Hagfors,et al.  A characterization of a comet nucleus interior:: inversion of simulated radio frequency data , 1999 .

[34]  Stephan Ulamec,et al.  Consert for Asteroid - Radar Tomography for Near Earth Asteroid , 2010 .

[35]  Daniela Calvetti,et al.  Hypermodels in the Bayesian imaging framework , 2008 .

[36]  Michael J. S. Belton,et al.  Mitigation of Hazardous Comets and Asteroids , 2011 .

[37]  Vito Pascazio,et al.  3D Markov Random Field in realistic inverse scattering , 2010, 2010 IEEE International Geoscience and Remote Sensing Symposium.

[38]  Daniela Calvetti,et al.  Introduction to Bayesian Scientific Computing: Ten Lectures on Subjective Computing , 2007 .

[39]  Vito Pascazio,et al.  Bayesian Regularization in Nonlinear Imaging: Reconstructions From Experimental Data in Nonlinearized Microwave Tomography , 2011, IEEE Transactions on Geoscience and Remote Sensing.

[40]  F. Harris.  On the use of windows for harmonic analysis with the discrete Fourier transform , 1978, Proceedings of the IEEE.

[41]  James Irving,et al.  Numerical modeling of ground-penetrating radar in 2-D using MATLAB , 2006, Comput. Geosci..

[42]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[43]  Sampsa Pursiainen,et al.  Iterative alternating sequential (IAS) method for radio tomography of asteroids in 3D , 2013 .

[44]  E. Grün,et al.  Comet nucleus sounding experiment by radiowave transmission , 1998 .

[45]  Wlodek Kofman,et al.  A priori information required for a two or three dimensional reconstruction of the internal structure of a comet nucleus (consert experiment) , 2002 .