acquiring, via control device circuitry using a plurality of reception coils as part of blood oxygenation level dependent (BOLD) functional MR imaging, (i) MR signals using parallel imaging, and (ii) a MR sequence defining a k-space trajectory in which undersampling in at least two k-space directions is performed;

reconstructing, via the control device circuitry using a reconstruction technique for undersampled MR data, the functional MR data set from (i) the MR signals, and (ii) sensitivity information regarding the plurality of reception coils; and

generating, via the control device circuitry, MR images of the imaging region of the brain of the patient based upon the functional MR data set,

wherein the k-space trajectory is selected to enable a controlled aliasing in each one of three spatial dimensions, including a readout direction, and

wherein reconstructing the functional MR data set comprises using a low rank plus sparse reconstruction technique in which dynamic MR data as a space-time matrix is used as a linear superposition of a spatially- and temporally-correlated image background matrix and a sparse dynamic information matrix.