A method registers images by first determining spatial gradient magnitudes for each pixel in each image to provide a corresponding energy image. Each energy image is transformed according to motion parameters associated with the energy image. An average sum of the transformed energy images is maximized. The motion parameters are then updated according to the maximized average sum. The above steps are repeated until a termination condition is reached, and then the images can be registered with each other according to the updated motion parameters.