Reverse time migration (RTM) is the important intermediate step for focusing the radar diffracted energy of the targets in ground penetrating radar imaging. The conventional RTM algorithm demands a large number of iterative trial experiments and depends on the experts' decision on the estimation of the velocity or relative permittivity of the subsurface medium. Meanwhile, the RTM profile is vulnerable to artifacts, which are composed of noise interference, multiple interferences, arc-shaped clutter, and crosstalk, so it is difficult to inspect visually. Therefore, we propose a RTM method based on accurate velocity estimation and total variation (TV) de-noising to improve the accuracy of the RTM imaging. Firstly, the appropriate migration velocity is obtained automatically by autofocusing metrics to reduce the number of visual inspection times and corrections in the migration processing. Secondly, the TV de-noising strategy based on split Bregman iterative is applied to the RTM profile with the cross-correlation imaging condition, so that the edge of the target can be obvious and the position can be accurate. Then, we apply the proposed method to the simulation data of the pipeline model and the tunnel lining model. All results show that the selected three different autofocusing metrics have unimodality and unbiasedness, which can focus on a single relative permittivity to obtain appropriate migration velocity. Furthermore, the TV de-noising strategy successfully eliminates artifacts, reconstructs contours, enhances the edge sharpness, and improves the quality and accuracy of the radar profile. Finally, we take the field data of the LiuYang River tunnel to verify the applicability of our method, we choose 2–30 as a wider range of relative permittivity based on prior information. Considering the different lateral velocity of the field data, we adopt lateral segmentation processing to improve the quality of the GPR profile. The tests of simulation data and the field data indicate that the proposed method can provide a scientific and effective way for accurate interpretation of GPR data.