A synthetic aperture radar imaging method that combines each radar return pulse with a sinusoid to reduce the radar return pulses to a baseband frequency and deskew each radar return pulse. It includes determining a maximum likelihood estimate (MLE) of residual motion parameters for a dominant scatterer on the ground relative to the airborne radar and correcting for errors in inertial navigation system measurements based on the MLE residual motion parameters. It includes convolving each radar return pulse with its corresponding radar transmission pulse to generate a range compressed image for each radar return pulse and generating a sub-band range profile image for each radar return pulse and its corresponding radar transmission pulse based on the corresponding range compressed image that has been corrected for residual motion. Performing bandwidth extrapolation on each sub-band and subsequently combining the three bands to produce an enhanced resolution image without grating lobes.