A three dimensional projection reconstruction pulse sequence is employed to acquire velocity encoded NMR data from which an image indicative of spin motion is reconstructed. The velocity encoding is along all three axes and it may include acquisitions at more than one velocity encoding first moment M1. When more than one first moment M1 is acquired, a 1DFT along the velocity encoding axis is performed prior to reconstructing images from the acquired NMR data.