Consider the inverse scattering of time-harmonic point sources by an infinite
rough surface which is supposed to be a local perturbation of a plane. A novel
version of reverse time migration is proposed to reconstruct the shape and
location of the rough surface. The method is based on a modified
Helmholtz-Kirchhoff identity associated with a special rough surface, leading
to a modified imaging functional which always reaches a peak on the boundary of
the rough surface for sound-soft case and penetrable case, and hits a nadir on
the boundary of the rough surface for sound-hard case. Numerical experiments
are presented to show the powerful imaging quality, especially for
multi-frequency data.Comment: 29 pages,26 figure