An iterative multiscaling approach for solving the electromagnetic inverse scattering problem related to the imaging of shallow subsurface targets with the ground-penetrating radar (GPR) is proposed. The approach combines the zooming properties of the multiscaling technique with the reconstruction capabilities of an Inexact-Newton (IN) method developed in <inline-formula> <tex-math notation="LaTeX">L^{p} </tex-math></inline-formula> spaces. It is based on multifrequency processing that allows one to face the ill-posedness of the inverse scattering problem by exploiting the regularization properties of a truncated Landweber (LW) method. Experimental data, extracted from radargrams obtained by the GPR in a real situation, are used for validation. The reconstruction results are also compared with those from competitive alternatives, such as a standard IN method or a state-of-the-art multifrequency Conjugate Gradient (CG)-based approach.