In the last decade, activity-dependent strategies for labelling multiple immediate early gene (IEG) ensembles in mice have generated unprecedented insight into the mechanisms of memory encoding, storage, and retrieval. However, few strategies exist for brain-wide mapping of multiple ensembles, including their overlapping population, and none incorporate capabilities for downstream network analysis. Here, we introduce a scalable workflow to analyze traditionally coronally-sectioned datasets produced by activity-dependent tagging systems. Intrinsic to this pipeline is