Objective. Ultrafast ultrasound imaging using planar or diverging waves for transmission is a promising approach for efficient 3D imaging with matrix arrays. This technique has advantages for B-mode imaging and advanced techniques, such as 3D vector flow imaging (VFI). The computation load of the cross-beam technique is associated with the number of transmit anglesmand receive anglesn. The full velocity vector is obtained using the least square fashion. However, the beamforming is repeatedm × ntimes using a conventional time-domain delay-and-sum (DAS) beamformer. In the 3D case, the collection and processing of data from different beams increase the amount of data that must be processed, requiring more storage capacity and processing power. Furthermore, the large computation complexity of DAS is another major concern. These challenges translate into longer computational times, increased complexity in data processing, and difficulty in real-time applications.Approach. In response to this issue, this study proposes a novel Fourier domain beamformer for 3D plane wave imaging, which significantly increases the computational speed. Additionally, a selective compounding strategy is proposed for VFI, which reduces the beamforming process fromm × ntom(wheremandnrepresent the number of transmission and reception, respectively), effectively shortening the processing time. The underlying principle is to decompose the receive wavefront into a series of plane waves with different slant angles. Each slant angle can produce a sub-volume for coherent or selective compounding. This method does not rely on the assumption that the plane wave is perfect and the results show that our proposed beamformer is better than DAS in terms of resolution and image contrast. In the case of velocity estimation, for the Fourier-based method, only Tx angles are assigned in the beamformer and the selective compounding method produces the final image with a specialized Rx angle.Main results. Simulation studies andin vitroexperiments confirm the efficacy of this new method. The proposed beamformer shows improved resolution and contrast performance compared to the DAS beamformer for B-mode imaging, with a suppressed sidelobe level. Furthermore, the proposed technique outperforms the conventional DAS method, as evidenced by lower mean bias and standard deviation in velocity estimation for VFI. Notably, the computation time has been shortened by 40 times, thus promoting the real-time application of this technique. The efficacy of this new method is verified through simulation studies andin vitroexperiments and evaluated by mean bias and standard deviation. Thein vitroresults reveal a better velocity estimation: the mean bias is 2.3%, 3.4%, and 5.0% forvx,vy, andvz, respectively. The mean standard deviation is 1.8%, 1.7%, and 3.4%. With DAS, the evaluated mean bias is 9.8%, 4.6%, and 6.7% and the measured mean standard deviation is 7.5%, 2.5%, and 3.9%.Significance. In this work, we propose a novel Fourier-based method for both B-mode imaging and functional VFI. The new beamformer is shown to produce better image quality and improved velocity estimation. Moreover, the new VFI computation time is reduced by 40 times compared to conventional methods. This new method may pave a new way for real-time 3D VFI applications.