This study concentrates on advancing mathematical and computational methodology for radar tomography imaging in which the unknown volumetric velocity distribution of a wave within a bounded domain is to be reconstructed. Our goal is to enable effective simulation and inversion of a large amount of full-wave data within a realistic 2D or 3D geometry. For propagating and inverting the wave, we present a rigorous multigrid-based forward approach which utilizes the finitedifference time-domain method and a nested finite element grid structure. We also introduce and validate a multigrid-based inversion algorithm which allows regularization of the unknown distribution through a coarseto-fine inversion scheme. Using this approach, sparse signals can be effectively inverted, as the coarse fluctuations are reconstructed before the finer ones. Furthermore, the number of nonzero entries in the system matrix can be compressed and, thus, the inversion procedure can be speeded up. As the test scenario, we investigate satellite-based asteroid interior reconstruction. We use both full-wave and projected wave data and estimate the accuracy of the inversion under different error sources: noise and positioning inaccuracies. The results suggest that the present inversion technique allows recovering the interior with a single satellite recording backscattering data. Robust results can be achieved, when the peak-to-peak signal-to-noise ratio is above 10 dB. Furthermore, the robustness for the deep interior part can be enhanced if two satellites can be utilized in the measurements.