A series of velocity encoded MR image frames are acquired. To increase the temporal resolution of the acquired image frames radial projections are acquired and each image frame is highly undersampled. The radial projections for each velocity encoding direction are interleaved throughout the scan and a composite phase image is reconstructed from these and used to reconstruct a velocity image for each image frame in a highly constrained backprojection method.