A method for automated magnetic field shimming in a MR imaging system includes generating RF excitation pulses in anatomy, generating slice select magnetic field gradients on a static magnetic field and generating a plurality of shimming routines. The method also includes acquiring maps of the static magnetic field between each shimming routine and determining a base function for each of the plurality of shimming coils from: (i) intermediate magnitudes and directions of the static magnetic field produced by each of the plurality of shimming coils during each shimming routine; and (ii) intermediate currents supplied to each of the plurality of shimming coils during each shimming routine. The method further includes determining updated magnitudes and directions of the static magnetic field produced by each of the plurality of shimming coils by expanding the determined base function for each of the plurality of shimming coils into low order polynomials.