A conventional 3D phase contrast acquisition generates images with good spatial resolution, but often gives rise to artifacts due to pulsatile flow. 2D cine phase contrast, on the other hand, can register dynamic flow, but has a poor spatial resolution perpendicular to the imaging plane. A combination of both high spatial and temporal resolution may be advantageous in some cases, both in quantitative flow measurements and in MR angiography. The described 3D cine phase contrast pulse sequence creates a temporally resolved series of 3D data sets with velocity encoded data.