The mechanical properties and fracture propagation of Longmaxi shale loading under uniaxial compression were measured using eight cylindrical shale specimens (4 mm in diameter and 8 mm in height), with the bedding plane oriented at 0 • and 90 • to the axial loading direction, respectively, by micro computed tomography (micro-CT). Based on the reconstructed three-dimensional (3-D) CT images of cracks, different stages of the crack growth process in the 0 • and 90 • orientation specimen were revealed. The initial crack generally occurred at relatively smaller loading force in the 0 • bedding direction specimen, mainly in the form of tensile splitting along weak bedding planes. Shear sliding fractures were dominant in the specimens oriented at 90 • , with a small number of parallel cracks occurring on the bedding plane. The average thickness and volume of cracks in the 90 • specimen is higher than those for the specimen oriented at 0 • . The geometrical characterization of fractures segmented from CT scan binary images shows that a specific surface area correlates with tortuosity at the different load stages of each specimen. The 3-D box-counting dimension (BCD) calculations can accurately reflect crack evolution law in the shale. The results indicate that the cracks have a more complex pattern and rough surface at an orientation of 90 • , due to crossed secondary cracks and shear failure.Various experimental methods, e.g., electromagnetic radiation (EMR), acoustic emission (AE), micro-seismic (MS), infrared technique, high speed digital imaging, and computerized tomography (CT), have been used to analyze the mechanical behavior of rocks (including shale) during uniaxial or triaxial loading [15][16][17][18][19][20][21]. CT scanning has been widely applied in the investigation of damage propagation since 1997, due to the three-dimensional (3-D) visualization and high-resolution imaging [22][23][24][25][26][27][28]. Kawakata et al. used CT to study the initial mesoscopic damage characteristics of rock materials, and observed spatial fault development in granite undergoing compression [29,30]. In the experiment, only the damaged specimens, rather than the entire process from crack initiation to specimen destruction, were scanned, and then observed by 3-D reconstruction of X-ray CT images. To examine the entire process, Ge et al. investigated rock compression with a real-time CT test [23]. Through the use of the industrial CT images, the rock meso-damage propagation law was revealed. Rock failure started with micro-pore and micro-crack compression, growth, bifurcation, and development, followed by fracture. Sun et al. displayed real-time deformation in backfill body under loading using industrial CT in several slices, and provided insights into the process of faulting [31]. However, the relatively low spatial resolution and image quality of conventional industrial CT is not suitable for fine-grained rocks. Furthermore, the 3-D information was not fully captured or used, which limited the accuracy of further calculation ...