A majorize-minimize (MM) mathematical principle is applied to least squares regularization estimation problems to effect efficient processing of image data sets to provide good quality images. In a ground penetrating radar application, these approaches can reduce processing time and memory use by accounting for a symmetric nature of a given radar pulse, accounting for similar discrete time delays between transmission of a given radar pulse and reception of reflections from the given radar pulse, and accounting for a short duration of the given radar pulse.