Geological maps of Earth typically incorporate field observations of rock lithology, structure, composition, and more. In contrast, conventional planetary geological maps are often made using primarily qualitative morphostratigraphic remote sensing observations of planetary surfaces. However, it is possible to define independent quantitative spectral units (SUs) of planetary surfaces, which potentially contain information about surface composition, grain size, and space weathering exposure. Here, we demonstrate a generic method to combine independently derived geomorphic and SUs, using the Rachmaninoff basin, Mercury, as an example to create a new geostratigraphic map. From this geostratigraphic map, we can infer some compositional differences within geomorphic units, which clarifies and elaborates on the geological evolution of the region.