We propose a reliable direct imaging method based on the reverse time migration for finding extended obstacles with phaseless total field data. We prove that the imaging resolution of the method is essentially the same as the imaging results using the scattering data with full phase information when the measurement is far away from the obstacle. The imaginary part of the cross-correlation imaging functional always peaks on the boundary of the obstacle. Numerical experiments are included to illustrate the powerful imaging quality

This work pertains to numerical aspects of a finite element method based discontinuous functions. Our study focuses on the Interior Penalty Discontinuous Galerkin method (IPDGM) because of its high-level of flexibility for solving the full wave equation in heterogeneous media. We assess the performance of IPDGM through a comparison study with a spectral element method (SEM). We show that IPDGM is as accurate as SEM. In addition, we illustrate the efficiency of IPDGM when employed in a seismic imaging process by considering two-dimensional problems involving the Reverse Time Migration.