A fast and efficient method for reconstructing an image from undersampled, parallel MRI data sets acquired with non-Cartesian trajectories includes the calculation of unsampled k-space data from the acquired k-space data and sets of calculated reconstruction coefficients. To reduce the computation time, only a few reference reconstruction coefficients are calculated using a matrix inversion step and the remaining reconstruction coefficients are produced by interpolating between the reference reconstruction coefficients.