Abstract In our energy-consuming society, understanding magnetization reversal in permanent magnets is crucial for improving energy conversion efficiency between electric energy and mechanical energy. First-order reversal curves (FORCs) have enabled qualitative studies of reversal mechanisms, however, further understanding based on quantitative analysis is still difficult due to the complexity of FORC diagrams. We introduce a machine-learning based approach combining the Gaussian mixture model and the Davies-Bouldin index to separate characteristic features in FORC diagrams of Nd-Fe-B magnets. The clustering method is evaluated using several FORC diagrams obtained at different temperatures, and we also demonstrated hysteresis loop reconstruction using the clustered FORC diagrams.