A method for producing an image of an organ or tissue of a subject with a magnetic resonance imaging system comprising the steps of: a) injecting a contrast agent into a blood vessel of the subject which induces a differential between the susceptibilities of the blood vessel and the tissue of interest; b) acquiring MRI image data sets during a first passage of the contrast agent through the tissue of interest using an imaging sequence adapted to collect image data at both short and long echo times; c) computing image maps of relative cerebral blood volume (rCBV) from data in the acquired MRI image data sets; d) computing image maps of absolute cerebral blood volume, cerebral blood flow and mean transit time from data in the acquired image data sets; e) correcting the rCBV, CBV, CBF and MTT image maps for variations to T1, T2 and T2* relaxation rates caused by contrast agent extravasation or recirculation, wherein the correcting of the rCBV image maps includes computing the T2* relaxation rate R2*(t) at each time point from the image data signals collected at the short and long echo times; fitting a model to the R2*(t) data to account for confounding modifications to T2 or T2*; and selecting an arterial input function from the R2*(t) data, or wherein, if step b) comprises the acquisition of T2-weighted spin echo MRI signals, the correcting of the rCBV image maps includes computing the T2 relaxation rate R2(t) at each time point from the image data signals collected at the short and long echo times, fitting a model to the R2(t) data to account for the confounding modifications to T2 and selecting an arterial input function from the R2(t) data; f) creating T1-weighted DCE parameter maps for preselected DCE permeability parameters from the short echo time data, wherein a DCE concentration-time curve is derived from data in the acquired image data sets, wherein a two-step process for eliminating confounding T2* effects from the DCE concentration-time curve is performed, wherein, in a first step, 1/T2*(t) is estimated at each time point from short and long TE echo data and wherein, in a second step, a corrected short TE echo signal is obtained by extrapolating each time point of the short TE echo signal back to TE=0, and wherein a change in a spin lattice relaxation rate ΔR1(t) curve corrected for confounding T2* effects is calculated from the corrected short TE echo signal using an estimate of a pre-contrast spin lattice relaxation rate obtained from a corrected pre-contrast baseline signal and an estimate of a baseline signal determined from a first time point of the collected image data at short and long echo times, the change in the spin lattice relaxation rate ΔR1(t) curve then being used to determine the DCE concentration time curve; g) computing image maps that provide measures of tissue recirculation, tissue leakage volumes, cellular volumes and cellular density. a) injecting a contrast agent into a blood vessel of the subject which induces a differential between the susceptibilities of the blood vessel and the tissue of interest; b) acquiring MRI image data sets during a first passage of the contrast agent through the tissue of interest using an imaging sequence adapted to collect image data at both short and long echo times; c) computing image maps of relative cerebral blood volume (rCBV) from data in the acquired MRI image data sets; d) computing image maps of absolute cerebral blood volume, cerebral blood flow and mean transit time from data in the acquired image data sets; e) correcting the rCBV, CBV, CBF and MTT image maps for variations to T1, T2 and T2* relaxation rates caused by contrast agent extravasation or recirculation, wherein the correcting of the rCBV image maps includes computing the T2* relaxation rate R2*(t) at each time point from the image data signals collected at the short and long echo times; fitting a model to the R2*(t) data to account for confounding modifications to T2 or T2*; and selecting an arterial input function from the R2*(t) data, or wherein, if step b) comprises the acquisition of T2-weighted spin echo MRI signals, the correcting of the rCBV image maps includes computing the T2 relaxation rate R2(t) at each time point from the image data signals collected at the short and long echo times, fitting a model to the R2(t) data to account for the confounding modifications to T2 and selecting an arterial input function from the R2(t) data; f) creating T1-weighted DCE parameter maps for preselected DCE permeability parameters from the short echo time data, wherein a DCE concentration-time curve is derived from data in the acquired image data sets, wherein a two-step process for eliminating confounding T2* effects from the DCE concentration-time curve is performed, wherein, in a first step, 1/T2*(t) is estimated at each time point from short and long TE echo data and wherein, in a second step, a corrected short TE echo signal is obtained by extrapolating each time point of the short TE echo signal back to TE=0, and wherein a change in a spin lattice relaxation rate ΔR1(t) curve corrected for confounding T2* effects is calculated from the corrected short TE echo signal using an estimate of a pre-contrast spin lattice relaxation rate obtained from a corrected pre-contrast baseline signal and an estimate of a baseline signal determined from a first time point of the collected image data at short and long echo times, the change in the spin lattice relaxation rate ΔR1(t) curve then being used to determine the DCE concentration time curve; g) computing image maps that provide measures of tissue recirculation, tissue leakage volumes, cellular volumes and cellular density.