Based on the theory of digital image processing, geometric vector transformation and the principle of automatic generation of finite element meshes, a finite element-discrete element analysis method for digital image of geotechnical engineering materials is proposed. Combining the digital image technology with the finite element-discrete element method (FDEM) coupling analysis method proposed by Munjiza, a FDEM analysis system is developed to characterize the real heterogeneity of rock mass, which provides a new way to model rock mass and study the fracture mechanism of rock mass from a meso-scopic point of view. With the help of digital image technology, the system obtains the real meso-structure of rock material from the image of rock section. With the help of mature mesh generation technology, it maps it to the FDEM computational grid, thus overcoming the shortcomings of the original FDEM in considering material heterogeneity. Using this system, the Brazilian disk splitting test is simulated numerically, and the real fracture process of granite under load is reproduced. And we examined how transparency and task type influence trustfulness, perceived quality of work during human-autonomous system interaction. The numerical simulation results show that the stress distribution of rock samples considering heterogeneity is asymmetric, and the micro-structure of rock has an important influence on the crack propagation and stress distribution in rock.