For resonance image data of an imaged subject, a method that first detects and estimates the dominant motions of k-space data (i.e., the motion vectors) and then constructs a graphical model for each estimated motion vector. The segments of the k-space that are determined to be corrupted by motion are restored by minimizing the energy associated with the corresponding graphical model. Consequently, the MR image of the imaged subject becomes free of motion artifacts.