MRI is widely used for noninvasive hemodynamic-based functional brain imaging. In traditional spatial encoding, however, gradient switching limits the temporal resolution, which makes it difficult to unambiguously identify possible fast nonhemodynamic changes. In this paper we propose a novel reconstruction approach, called dynamic inverse imaging (InI), that is capable of providing millisecond temporal resolution when highly parallel detection is used. To achieve an order-of-magnitude speedup in generating time-resolved contrast estimates and dynamic statistical parametric maps (dSPMs), the spatial information is derived from an array of detectors rather than by time-consuming gradient-encoding methods. The InI approach was inspired by electroencephalography (EEG) and magnetoencephalography ( Key words: MRI; inverse; parallel MRI; magnetoencephalography; minimum-norm; MNE; dSPM; InI; electroencephalgraphy Single-shot fast MRI with 1-3-s time resolution has been the principal technology for functional MRI (fMRI) (1-4). This time resolution is adequate for observing hemodynamic responses (HDRs), which are secondary to the much faster neural activity (5). Because of its high spatial resolution, noninvasiveness, and flexibility of contrast preparation, MRI has also been considered for its potential to further detect fast, nonhemodynamic functional changes related to human brain activity (6,7). However, efforts to identify suitable MRI contrast mechanisms have been hampered by the lack of required temporal resolution. Ultrafast MRI (100 frames per second or more) is a mandatory tool for testing any novel direct neuronal contrast mechanisms. It may not be possible to discern a contrast change that occurs and then disappears on a 10-ms scale without an imaging technique with similar temporal resolution. Thus, a fast imaging technique with millisecond temporal resolution is indispensable as the prerequisite for allowing MR to probe direct neural contrast mechanisms.In traditional Fourier MRI, the temporal resolution is constrained by the k-space traversing scheme used to spatially encode the data. Echo-planar imaging (EPI) (8) and spiral imaging (9) utilize switching gradients to achieve k-space traversing at a faster imaging rate than gradientecho or spin-echo imaging. With state-of-the-art technology, 2D single-slice T 2 *-weighted images can be obtained in approximately 80 ms by EPI or spiral imaging. This allows whole-head coverage with 3-mm isotropic resolution in 2-3 s. Modest improvements can be made by optimizing the k-space sampling and reconstruction methods. For example, instead of completing the k-space traversal for every measurement, accelerated MR data acquisition can be achieved by various alterations of the k-space trajectories and the associated image reconstruction algorithms, such as partial k-space sampling (10). Alternatively, a priori information-based methods can also improve the temporal resolution of MR dynamic measurements (11) with various implementations, such as keyhole imaging (12), singu...