in combination with cognitive tasks entailing sequences of sensory and cognitive processes, event-related acquisition schemes allow using functional MRI to examine not only the topography but also the temporal sequence of cortical activation across brain regions (time-resolved fMRI). In this study, we compared two data-driven methods - fuzzy clustering method (FCM) and independent component analysis (ICA) - in the context of time-resolved fMRI data collected during the performance of a newly devised visual imagery task. We analyzed a multisubject fMRI data set using both methods and compared their results in terms of within and between-subject consistency and spatial and temporal correspondence of obtained maps and time courses. Both FCM and spatial ICA allowed discriminating the contribution of distinct networks of brain regions to the main cognitive stages of the task (auditory perception, mental imagery and behavioural response), with good agreement across methods. Whereas ICA worked optimally on the original time series, averaging with respect to the task onset (and thus introducing some a priori information on the stimulation protocol) was found to be indispensable in the case of FCM. On averaged time series, FCM led to a richer decomposition of the spatio-temporal patterns of activation and allowed a finer separation of the neurocognitive processes subserving the mental imagery task. This study confirms the efficacy of the two examined methods in the data-driven estimation of hemodynamic responses in time-resolved fMRI studies and provides empirical guidelines to their use.