1. A higher-order Generalized series parallel imaging method for acquiring high spatio-temporal resolution functional magnetic resonance images comprising: performing restricted sampling of a central part of an input image in k-space, so as to acquire data with a rate lower than a Nyquist rate; applying a high-order generalized series (HGS) reconstruction procedure to data acquired as the sampling result to acquire a first reconstructed image; and applying a parallel magnetic resonance reconstruction procedure to the first reconstructed image to acquire a second reconstructed image.