To analyze the hazard-causing modes of landslides, this paper proposes a three-dimensional discrete element model reconstruction method that employs an unmanned aerial vehicle survey and multi-electrode resistivity tomography method. To convert the resistivity profile into a material profile, we adopt the peak of the probability density method for material classification and utilize the Haar wavelet transform for image denoising. Subsequently, inverse distance weighting interpolation and the curtain-point method are used to transform two-dimensional profiles into a 3D visualization model. Similarly, the triangular mesh boundary can be extracted from the 3D visualization model using the curtain-point method. A mapping function f including the macroscopic parameters, was defined to populate the particles within the boundaries. Using the iterative method and defining the loss function L for parameter calibration, the targeted 3D discrete element model was constructed after setting the velocity threshold. This method was applied to the Changhe landslide (September 14, 2019) in Gansu Province, China, which had a typical damaged soil layer due to earthquake and rainfall factors. The results indicate that the lower part first exhibits significant displacement, followed by the upper and middle parts, which is consistent with the on-site inspections and UAV findings.