Multivariate analytical routines have become increasingly popular in the study of cerebral function in health and in disease states. Spatial covariance analysis of functional neuroimaging data has been used to identify and validate characteristic topographies associated with specific brain disorders. Voxel-wise correlations can be used to assess similarities and differences that exist between covariance topographies. While the magnitude of the resulting topographical correlations is critical, statistical significance can be difficult to determine in the setting of large data vectors (comprised of over 100,000 voxel weights) and substantial autocorrelation effects. Here, we propose a novel method to determine the p-value of such correlations using pseudo-random network simulations.