Neuropeptides are a chemically diverse class of cell-to-cell signaling molecules that are widely expressed throughout the central nervous system, often in a cell-specific manner. While cell-to-cell differences in neuropeptides is expected, it is often unclear how exactly neuropeptide expression varies among neurons. Here we created a microscopy-guided, high-throughput single cell matrix-assisted laser desorption/ionization mass spectrometry approach to investigate the neuropeptide heterogeneity of individual neurons in the central nervous system of the neurobiological model Aplysia californica, the California sea hare. In all, we analyzed more than 26,000 neurons from 18 animals and assigned 866 peptides from 66 prohormones by mass matching against an in silico peptide library generated from known Aplysia prohormones retrieved from the UniProt database. Louvain-Jaccard (LJ) clustering of mass spectra from individual neurons revealed 40 unique neuronal populations, or LJ clusters, each with a distinct neuropeptide profile. Prohormones and their related peptides were generally found in single cells from ganglia consistent with the prohormones' previously known ganglion localizations. Several LJ clusters also revealed the cellular colocalization of behaviorally related prohormones, such as an LJ cluster exhibiting achatin and neuropeptide Y, which are involved in feeding, and another cluster characterized by urotensin II, small cardiac peptide, sensorin A, and FRFa, which have shown activity in the feeding network or are present in the feeding musculature. This mass spectrometry-based approach enables the robust categorization of large cell populations based on single cell neuropeptide content and is readily adaptable to the study of a range of animals and tissue types.