A method of magnetic resonance imaging for the acquisition of a three-dimensional dataset of an object of investigation, where spatial encoding by three mutually orthogonal magnetic field gradients is performed such, that signal is readout under a read-gradient in one spatial direction k 1 and spatial encoding in the other two spatial directions k 2, k 3 is performed by applying phase encoding gradients in the other two spatial directions k 2, k 3 prior to signal acquisition and data acquisition is performed in a sequential manner such, that at each acquisition step signal is acquired under said readout gradient, but with different combinations of the two phase encoding gradients, wherein a complete dataset is acquired using a randomized sampling pattern in the k 2 -k 3 plane, which comprises two datasets: a basic dataset and a complementary dataset, wherein the temporal sequence of spatial encoding is applied such, that first a basic dataset containing N 3D(b) phase encoding steps is acquired using a randomized sampling pattern in the two phase encoding directions, where the sampling density of the randomized sampling pattern decreases monotonously from the center of k-space to generate a first dataset suitable for reconstruction of strongly undersampled data by use of the compressed sensing method to generate a first 3D image volume of the object of investigation with a predefined sufficient minimum quality, wherein data acquisition then proceeds by acquiring a complimentary dataset in which N 3D(c) additional phase encoding steps are acquired supplementing the already acquired data such, that the sampling density allows the reconstruction of 3D-image volumes of increasing image quality until data are acquired with a predefined sampling density to reconstruct a final 3D-volume of optimum quality, and wherein monitoring of motion of the object of investigation is applied such that data acquisition is terminated whenever involuntary motion exceeds a predefined threshold, where motion monitoring is performed either using navigators built into the measurement sequence or by extrinsic motion detection devices, characterized in that spatial encoding in the acquisition of the basic as well as the complimentary dataset is performed according to a view ordering scheme in a segmented fashion in L segments with N seg phase encoding steps in each segment, wherein according to the view ordering scheme, the total k-space sampling pattern of the basic or complimentary dataset is subdivided into a predefined number of N seg concentric or rectilinear sectors containing L phase encoding steps each and for the acquisition of each phase encoding step in each segment one phase encoding step is selected from each sector such that each segment contains one phase encoding step from each sector, and that the numbers N seg and L are different between the basic and the complimentary acquisition and chosen such that N seg and L are numbers and that N seg ∗ L≈N 3D(b) and N 3D(c) respectively. wherein a complete dataset is acquired using a randomized sampling pattern in the k 2 -k 3 plane, which comprises two datasets: a basic dataset and a complementary dataset, wherein the temporal sequence of spatial encoding is applied such, that first a basic dataset containing N 3D(b) phase encoding steps is acquired using a randomized sampling pattern in the two phase encoding directions, where the sampling density of the randomized sampling pattern decreases monotonously from the center of k-space to generate a first dataset suitable for reconstruction of strongly undersampled data by use of the compressed sensing method to generate a first 3D image volume of the object of investigation with a predefined sufficient minimum quality, wherein data acquisition then proceeds by acquiring a complimentary dataset in which N 3D(c) additional phase encoding steps are acquired supplementing the already acquired data such, that the sampling density allows the reconstruction of 3D-image volumes of increasing image quality until data are acquired with a predefined sampling density to reconstruct a final 3D-volume of optimum quality, and wherein monitoring of motion of the object of investigation is applied such that data acquisition is terminated whenever involuntary motion exceeds a predefined threshold, where motion monitoring is performed either using navigators built into the measurement sequence or by extrinsic motion detection devices, characterized in that spatial encoding in the acquisition of the basic as well as the complimentary dataset is performed according to a view ordering scheme in a segmented fashion in L segments with N seg phase encoding steps in each segment, wherein according to the view ordering scheme, the total k-space sampling pattern of the basic or complimentary dataset is subdivided into a predefined number of N seg concentric or rectilinear sectors containing L phase encoding steps each and for the acquisition of each phase encoding step in each segment one phase encoding step is selected from each sector such that each segment contains one phase encoding step from each sector, and that the numbers N seg and L are different between the basic and the complimentary acquisition and chosen such that N seg and L are numbers and that N seg ∗ L≈N 3D(b) and N 3D(c) respectively.