The invention discloses a complex
earth surface ground penetrating radar amplitude-preserving imaging method based on an orthogonal body-fitted grid, and the method comprises the following steps: 1), constructing an
earth surface model according to the ground fluctuation data; constructing an electromagnetic wave underground propagation velocity model according to
geological exploration data; (2) an orthogonal body-fitted grid is generated by means of a Ryskin-Leal method; 3) discretizing the orthogonal body-fitted mesh by using a
finite difference time domain method; 4) performing reverse extrapolation on the received
ground penetrating radar signal, and reversely propagating an electromagnetic
wave field to an initial moment; in the reverse extrapolation process, the spatial position of a reflection interface is calculated through cross-correlation imaging conditions, and an underground imaging profile is obtained; 5) calculating the gradient of a target function, solving an inversion problem by adopting a
conjugate gradient method, and iteratively updating the electromagnetic wave underground propagation velocity model; and 6) generating an imaging profile of the underground structure according to the electromagnetic wave underground propagation velocity model after iterative convergence. According to the invention, the precision, the stability and the anti-
noise performance of underground imaging are obviously improved.