A method for generating a magnetic resonance (MR) image includes applying a pulse sequence including a quadratic field gradient. A first k-space data set is acquired from each of a plurality of RF coils where each first k-space data set including uniformly undersampled data. A randomly undersampled k-space data set is generated for each RF coil from the first k-space data set. A compressed sensing reconstruction technique is applied to the randomly undersampled k-space data set of each RF coil to generate a second k-space data set for each RF coil where each second k-space data set including uniformly undersampled data. A phase scrambling reconstruction technique is applied to the second k-space data set of each RF coil to generate a low resolution coil image for each RF coil. A MR image is generated by applying a parallel imaging technique to the low resolution coil image and second k-space data set for each RF coil.