Elemental and isotopic concentrations are processed using the software HDIP (version 1.6., Teledyne Photon Machines). The background was fitted using a cubic spline function for each mass channel and drift is corrected linearly using NIST610 as a drift monitor. For the elemental concentration maps, each mass sweep is transformed into a quadratic pixel and is calibrated using NIST610 line scan analyses. The assignment of pixels to specific mineral phases (muscovite, biotite, quartz, plagioclase, K-Fsp) is achieved using the labeled segments tool within HDIP. This tool provides flexibility in pixel selection, allowing manual selection, cluster analyses based on a k-means algorithm, or user-defined thresholds using logical operators. In the case of the k-means algorithm, a predefined number of clusters is necessary, which we determined based on the expected mineral phases. For instance, in the GER2B dataset, we used 3 clusters for quartz, plagioclase, and muscovite, while in the GÄD 08 dataset, 4 clusters were utilized to represent quartz, plagioclase, muscovite, and biotite. The k-means algorithm partitions the data points into clusters by assigning each point to the cluster whose mean (centroid) is closest to the point. This iterative process refines cluster assignments by recalculating centroids and reassigning data points until the clusters' positions converge and remain stable. To account for mixed analyses occurring at grain boundaries, an additional cluster per map was added. The clusters were generated by taking into account a specific mass channels ( 28 Si, 85 Rb, 88 Sr, and 87 Sr for GÄD08 and 7 Li, 27Al , 28 Si, 85 Rb, 88 Sr, and 87 Sr for GER2B).