The probabilistic topography and inter-individual variability of the pituitary gland (PG) remain undetermined. The absence of a standardized reference atlas hinders research on PG volumetrics. In this study, we aimed at creating maximum probability maps for the anterior and posterior PG in young female adults. We manually delineated the anterior and posterior parts of the pituitary glands in 26 healthy subjects using high-resolution MRI T1 images. A three-step procedure and a cost function-masking approach were employed to optimize spatial normalization for the PG. We generated probabilistic atlases and maximum probability maps, which were subsequently coregistered back to the subjects’ space and compared to manual delineations. Manual measurements led to a total pituitary volume of 705 ± 88 mm³, with the anterior and posterior volumes measuring 614 ± 82 mm³ and 91 ± 20 mm³, respectively. The mean relative volume difference between manual and atlas-based estimations was 1.3%. The global pituitary atlas exhibited an 80% (± 9%) overlap for the DICE index and 67% (± 11%) for the Jaccard index. Similarly, these values were 77% (± 13%) and 64% (± 14%) for the anterior pituitary atlas and 62% (± 21%) and 47% (± 17%) for the posterior PG atlas, respectively. We observed a substantial concordance and a significant correlation between the volume estimations of the manual and atlas-based methods for the global pituitary and anterior volumes. The maximum probability maps of the anterior and posterior PG lay the groundwork for automatic atlas-based segmentation methods and the standardized analysis of large PG datasets.