MALDI mass spectrometry can simultaneously measure hundreds of biomolecules directly from tissue. Using essentially the same technique but different sample preparation strategies, metabolites, lipids, peptides and proteins can be analyzed. Spatially correlated analysis, imaging MS, enables the distributions of these biomolecular ions to be simultaneously measured in tissues. A key advantage of imaging MS is that it can annotate tissues based on their MS profiles and thereby distinguish biomolecularly distinct regions even if they were unexpected or are not distinct using established histological and histochemical methods e.g. neuropeptide and metabolite changes following transient electrophysiological events such as cortical spreading depression (CSD), which are spreading events of massive neuronal and glial depolarisations that occur in one hemisphere of the brain and do not pass to the other hemisphere , enabling the contralateral hemisphere to act as an internal control. A proof-of-principle imaging MS study, including 2D and 3D datasets, revealed substantial metabolite and neuropeptide changes immediately following CSD events which were absent in the protein imaging datasets. The large high dimensionality 3D datasets make even rudimentary contralateral comparisons difficult to visualize. Instead non-negative matrix factorization (NNMF), a multivariate factorization tool that is adept at highlighting latent features, such as MS signatures associated with CSD events, was applied to the 3D datasets. NNMF confirmed that the protein dataset did not contain substantial contralateral differences, while these were present in the neuropeptide dataset.