A method for magnetic resonance, MR, imaging comprising the steps of: 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 as 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 as non-compensated data; the method being characterized in that : a joint analysis of the flow-compensated data and non-compensated data is made comprising performing a model fit to each one of the flow-compensated and non-compensated data in which at least one adjustable model parameter of each model fit is constrained to be the same for both sets of data; and obtaining quantitative information on microscopic flow by extracting at least one parameter of an intravoxel incoherent motion, IVIM, effect from the model fits, wherein a diffusion weighting, b, and a velocity dispersion weighting, α , are adjusted independently of each other between acquiring the flow-compensated and non-compensated data, as means for performing a diffusion-perfusion correlation experiment, wherein α = γ ∫ 0 t E G t ' t 'd t ' = − ∫ 0 t E q t ' d t ' , where yis the nuclear gyromagnetic ratio, G ( t ) is the magnetic field gradient, t E is the echo time and q ( t ) the dephasing factor. 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 as 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 as non-compensated data; the method being characterized in that : a joint analysis of the flow-compensated data and non-compensated data is made comprising performing a model fit to each one of the flow-compensated and non-compensated data in which at least one adjustable model parameter of each model fit is constrained to be the same for both sets of data; and obtaining quantitative information on microscopic flow by extracting at least one parameter of an intravoxel incoherent motion, IVIM, effect from the model fits, wherein a diffusion weighting, b, and a velocity dispersion weighting, α , are adjusted independently of each other between acquiring the flow-compensated and non-compensated data, as means for performing a diffusion-perfusion correlation experiment, wherein α = γ ∫ 0 t E G t ' t 'd t ' = − ∫ 0 t E q t ' d t ' , where yis the nuclear gyromagnetic ratio, G ( t ) is the magnetic field gradient, t E is the echo time and q ( t ) the dephasing factor. a joint analysis of the flow-compensated data and non-compensated data is made comprising performing a model fit to each one of the flow-compensated and non-compensated data in which at least one adjustable model parameter of each model fit is constrained to be the same for both sets of data; and obtaining quantitative information on microscopic flow by extracting at least one parameter of an intravoxel incoherent motion, IVIM, effect from the model fits, wherein a diffusion weighting, b, and a velocity dispersion weighting, α , are adjusted independently of each other between acquiring the flow-compensated and non-compensated data, as means for performing a diffusion-perfusion correlation experiment, wherein α = γ ∫ 0 t E G t ' t 'd t ' = − ∫ 0 t E q t ' d t ' , where yis the nuclear gyromagnetic ratio, G ( t ) is the magnetic field gradient, t E is the echo time and q ( t ) the dephasing factor. <EOS>