A method includes receiving a forward spatial encoding polarity magnetic resonance (MR) coil image and a reverse spatial encoding polarity MR coil image generated from data obtained with a magnetic field gradient that is reversed with respect to the magnetic field gradient with which the forward spatial encoding polarity MR coil image is acquired. The method also includes performing an iterative shift map calculation algorithm to determine a pixel shift map corresponding to a minimized difference between the forward and reverse spatial encoding polarity MR coil images, converting the pixel shift map into a magnetic field shift map by determining a magnetic field value corresponding to each pixel in the pixel shift map, and providing the magnetic field shift map as an input to a shim calculation process that includes determining a level of at least one shim current.