Temporal Non-Local Means Filtering Reveals Real-Time Whole-Brain Cortical Interactions in Resting fMRI

Intensity variations over time in resting BOLD fMRI exhibit spatial correlation patterns consistent with a set of large scale cortical networks. However, visualizations of this data on the brain surface, even after extensive preprocessing, are dominated by local intensity fluctuations that obscure larger scale behavior. Our novel adaptation of non-local means (NLM) filtering, which we refer to as temporal NLM or tNLM, reduces these local fluctuations without the spatial blurring that occurs when using standard linear filtering methods. We show examples of tNLM filtering that allow direct visualization of spatio-temporal behavior on the cortical surface. These results reveal patterns of activity consistent with known networks as well as more complex dynamic changes within and between these networks. This ability to directly visualize brain activity may facilitate new insights into spontaneous brain dynamics. Further, temporal NLM can also be used as a preprocessor for resting fMRI for exploration of dynamic brain networks. We demonstrate its utility through application to graph-based functional cortical parcellation. Simulations with known ground truth functional regions demonstrate that tNLM filtering prior to parcellation avoids the formation of false parcels that can arise when using linear filtering. Application to resting fMRI data from the Human Connectome Project shows significant improvement, in comparison to linear filtering, in quantitative agreement with functional regions identified independently using task-based experiments as well as in test-retest reliability.