A Novel, Efficient Algorithm for Subsurface Radar Imaging below a Non-Planar Surface

In classical radar imaging, such as in Earth remote sensing, electromagnetic waves are usually assumed to propagate in free space. However, in numerous applications, such as ground penetrating radar or non-destructive testing, this assumption no longer holds. When there is a multi-material background, the subsurface image reconstruction becomes considerably more complex. Imaging can be performed in the spatial domain or, equivalently, in the wavenumber domain (k-space). In subsurface imaging, to date, objects with a non-planar surface are commonly reconstructed in the spatial domain, by the Backprojection algorithm combined with ray tracing, which is computationally demanding. On the other hand, objects with a planar surface can be reconstructed more efficiently in k-space. However, many non-planar surfaces are partly planar. Therefore, in this paper, a novel concept is introduced that makes use of the efficient k-space-based reconstruction algorithms for partly planar scenarios, too. The proposed algorithm forms an image from superposing sub-images where as many image parts as possible are reconstructed in the wavenumber domain, and only as many as necessary are reconstructed in the spatial domain. For this, a segmentation scheme is developed to determine which parts of the image volume can be reconstructed in the wavenumber domain. The novel concept is verified by measurements, both from monostatic synthetic aperture radar data and multiple-input–multiple-output radar data. It is shown that the computational efficiency for imaging irregularly shaped geometries can be significantly augmented when applying the proposed concept.