Described here are systems and methods for retrospectively estimating and correcting for rigid-body motion by using a joint optimization technique to jointly solve for motion parameters and the underlying image. This method is implemented for magnetic resonance imaging (“MRI”), but can also be adapted for other imaging modalities.