Identifying patterns from the neuroimaging recordings of brain activity related to the unobservable psychological or mental state of an individual can be treated as a unsupervised pattern recognition problem. The main challenges, however, for such an analysis of fMRI data are: a) defining a physiologically meaningful feature-space for representing the spatial patterns across time; b) dealing with the high-dimensionality of the data; and c) robustness to the various artifacts and confounds in the fMRI time-series. In this paper, we present a network-aware feature-space to represent the states of a general network, that enables comparing and clustering such states in a manner that is a) meaningful in terms of the network connectivity structure; b)computationally efficient; c) low-dimensional; and d) relatively robust to structured and random noise artifacts. This feature-space is obtained from a spherical relaxation of the transportation distance metric which measures the cost of transporting ``mass'' over the network to transform one function into another. Through theoretical and empirical assessments, we demonstrate the accuracy and efficiency of the approximation, especially for large problems. While the application presented here is for identifying distinct brain activity patterns from fMRI, this feature-space can be applied to the problem of identifying recurring patterns and detecting outliers in measurements on many different types of networks, including sensor, control and social networks.

Abstract Classically, analysis of the time‐varying data acquired during fMRI experiments is done using static activation maps obtained by testing voxels for the presence of significant activity using statistical methods. The models used in these analysis methods have a number of parameters, which profoundly impact the detection of active brain areas. Also, it is hard to study the temporal dependencies and cascading effects of brain activation from these static maps. In this paper, we propose a methodology to visually analyze the time dimension of brain function with a minimum amount of processing, allowing neurologists to verify the correctness of the analysis results, and develop a better understanding of temporal characteristics of the functional behaviour. The system allows studying time‐series data through specific volumes‐of‐interest in the brain‐cortex, the selection of which is guided by a hierarchical clustering algorithm performed in the wavelet domain. We also demonstrate the utility of this tool by presenting results on a real data‐set.