Carbonate rocks undergo low-temperature, post-depositional changes, including mineral precipitation, dissolution, or recrystallisation (diagenesis). Unravelling the sequence of these events is time-consuming, expensive, and relies on destructive analytical techniques, yet such characterization is essential to understand their post-depositional history for mineral and energy exploitation and carbon storage. Conversely, hyperspectral imaging offers a rapid, non-destructive method to determine mineralogy, while also providing compositional and textural information. It is commonly employed to differentiate lithology, but it has never been used to discern complex diagenetic phases in a largely monomineralic succession. Using spatial-spectral endmember extraction, we explore the efficacy and limitations of hyperspectral imaging to elucidate multi-phase dolomitization and cementation in the Cathedral Formation (Western Canadian Sedimentary Basin). Spectral endmembers include limestone, two replacement dolomite phases, and three saddle dolomite phases. Endmember distributions were mapped using Spectral Angle Mapper, then sampled and analyzed to investigate the controls on their spectral signatures. The absorption-band position of each phase reveals changes in %Ca (molar Ca/(Ca + Mg)) and trace element substitution, whereas the spectral contrast correlates with texture. The ensuing mineral distribution maps provide meter-scale spatial information on the diagenetic history of the succession that can be used independently and to design a rigorous sampling protocol.