1. A method for magnetic resonance (MR) imaging comprising: acquiring at least two MR images with different motion-weighting originating from a RF and gradient pulse sequence causing signal attenuation from diffusion but not flow (flow-compensated data); acquiring at least two MR images with different motion-weighting originating from a RF and gradient pulse sequence causing signal attenuation from diffusion and flow (non-compensated data); performing a model fit to the flow-compensated and non-compensated data in which at least one of the adjustable parameters are constrained to be the same for both sets of data; and obtaining quantitative information on microscopic flow by extracting at least one parameter of the intravoxel incoherent motion (IVIM) effect from the model fit.