The discrete fracture system of a rock mass plays a crucial role in controlling the stability of rock slopes. To fully account for the geometric shape and distribution characteristics of jointed rock masses, terrestrial laser scanning (TLS) was employed to acquire high-resolution point-cloud data, and a developed automatic discontinuity-identification technology was coupled to automatically interpret and characterize geometric information such as orientation, trace length, spacing, and set number of the discontinuities. The discrete element method (DEM) was applied to study the influence of the geometric morphology and distribution characteristics of discontinuities on slope stability by generating a discrete fracture network (DFN) with the same statistical characteristics as the actual discontinuities. Based on slope data from the Yebatan Hydropower Station, a simulation was conducted to verify the applicability of the automatic discontinuity identification technology and the discrete fracture network-discrete element method (DFN-DEM). Various geological parameters, including trace length, persistence, and density, were examined to investigate the morphological evolution and response characteristics of rock slope excavation under different joint combination conditions through simulation. The simulation results indicate that joint parameters affect slope stability, with density having the most significant impact. The impact of joint parameters on stability is relatively small within a reasonable range but becomes significant beyond a certain threshold, further validating that the accuracy of field geological surveys is critical for simulation. This study provides a scientific basis for the construction of complex rock slope models, engineering assessments, and disaster prevention and mitigation, which is of great value in both theory and engineering applications.