In ground-based astronomical observations or artificial space target detections, images obtained from a ground-based telescope are severely distorted due to atmospheric turbulence. The distortion can be partially compensated by employing adaptive optics (pre-detection compensation), image restoration techniques (post-detection compensation), or a combination of both (hybrid compensation). This paper focuses on the improvement of the most commonly used practical post-processing techniques, Richardson–Lucy (R–L) iteration blind deconvolution, which is studied in detail and improved as follows: First, the total variation (TV) norm is redefined using the Gaussian gradient magnitude and a set scheme for regularization parameter selection is proposed. Second, the Gaussian TV constraint is proposed to impose to the R–L algorithm. Last, the Gaussian TV R–L (GRL) iterative blind deconvolution method is finally presented, in which the restoration precision is visually increased and the convergence property is considerably improved. The performance of the proposed GRL method is tested by both simulation experiments and observed field data.