This paper describes a procedure used to characterize the three-dimensional ͑3D͒ grain shape of lunar soil and undertake simulations of lunar soil by image-based discrete element method ͑DEM͒. Given that detailed 3D grain-shape information is unavailable for real lunar soil, a simulant material, FJS-1, is used in this study. We use the high-resolution micro X-ray CT system at SPring-8, a synchrotron radiation facility in Japan, to visualize precise 3D images of the granular assembly of FJS-1. A newly developed imageanalysis procedure is then applied to identify individual grains. Using the obtained grain-shape data, a sufficient number of FJS-1 grains are directly modeled for DEM simulation using an efficient modeling scheme. A series of particle flow simulations are then performed with the modeled grains. The resulting slope angles are in good agreement with experimental results. We discuss the effect on the slope angle of grain parameters such as contact stiffness, restitution coefficient, and interparticle friction.