A method of generating an enhanced perfusion image comprising the use of a blind deconvolution algorithm and the adiabatic approximation to the Johnson and Wilson model (aaJW) and generation of an image, wherein the blind deconvolution algorithm and the aaJW model are used in the generation of values of the following parameters: voxel specific arterial input function cp[t] and voxel specific tissue residue function r[t].