Rock discontinuities significantly impact the mechanical and hydraulic behavior of rock masses. A crucial aspect of rock engineering involves classifying discontinuities with similar orientations into groups. For this purpose, clustering algorithms, such as K-Means and Fuzzy K-Means (FKM), have been employed. However, the outcomes of these algorithms are influenced by the selection of initial cluster centers. This paper proposes an improved FKM algorithm to automatically identify rock discontinuity sets based on the Fisher distribution (FFKM). The algorithm uses the Fisher distribution to generate and select appropriate initial cluster centers. The performance of FFKM was initially validated using a published data set, and its results were compared with other clustering methods commonly used for grouping discontinuities. Results demonstrated the superior performance of FFKM over the FKM algorithm, comparable to other methods. Subsequently, the proposed algorithm was employed to analyze a fracture data set sampled at an open-pit iron mine. The FFKM facilitated identifying the correct number of sets and produced results consistent with the fracture sets observed in the field. Finally, the algorithm was verified using an artificial discontinuity data set, and the results demonstrated that the method correctly identified the number of sets and provided discontinuity sets similar to the original data set. The FFKM algorithm offers significant advantages: it maintains the essential characteristics of the FKM algorithm, effectively addresses the challenge of selecting suitable initial cluster centers, requires only the expected number of discontinuity sets as an input parameter, processes data within an acceptable computation time, serves as a tool for defining the number of discontinuity sets, and mitigates the drawbacks of the stereographic projection method.