Imaging the internal structure of trunks via multiscale phase inversion of ground-penetrating radar data

Ground-penetrating radar (GPR) is a geophysical technique widely used in near-surface non-invasive detecting. It has the ability to obtaining a high-resolution internal structure of living trunks. Full wave inversion (FWI) has been widely used to reconstruct the dielectric constant and conductivity distribution for cross-well application. However, in some cases, the amplitude information is not reliable due to the antenna coupling, radiation pattern and other effects. We present a multiscale phase inversion (MPI) method, which largely matches the phase information by normalizing the magnitude spectrum; in addition, a natural multiscale approach by integrating the input data with different times is implemented to partly mitigate the local minimal problem. Two synthetic GPR datasets generated from a healthy oak tree trunk and from a decayed trunk are tested by MPI and FWI. Field GPR dataset consisting of 30 common shot GPR data are acquired on a standing white oak tree (Quercus alba); the MPI and FWI methods are used to reconstruct the dielectric constant distribution of the tree cross-section. Results indicate that MPI has more tolerance to the starting model, noise level and source wavelet. It can provide a more accurate image of the dielectric constant distribution compared to the conventional FWI.