The present invention relates to a magnetic resonance spectroscopic imaging (MRSI) method, specifically to a magnetic resonance spectroscopic imaging method with up to three spatial dimensions and one spectral dimension. Interleaving dynamically switched magnetic field gradients into the spectroscopic encoding scheme enables multi-region shimming in a single shot to compensate the spatially varying spectral line broadening resulting from local magnetic field gradients. The method also employs sparse spectral sampling with controlled spectral aliasing and nonlinear sampling density to maximize encoding speed, data sampling efficiency and sensitivity.