1. A functional magnetic resonance imaging method comprising the steps of: (a) detecting a physiological timing profile; (b) measuring cerebral k-space data; (c) ordering the measured cerebral k-space data into a unit cycle based on the detected physiological timing profile; (d) estimating corresponding physiological effects for individual k-space data points by non-linear least-square fitting of cerebral k-space data; (e) removing the estimated physiological effects from the cerebral k-space data by subtraction, providing cerebral k-space data from which the estimated physical effects have been removed; and (f) providing an image from said cerebral k-space data from which the estimated physical effects have been removed.