In the context of bone diseases research, recent works have highlighted the crucial role of the osteocyte system. This system, hosted in the lacuno-canalicular network (LCN), plays a key role in the bone remodeling process. However, few data are available on the LCN due to the limitations of current microscopy techniques, and have mainly only been obtained from 2D histology sections. Here we present, for the first time, an automatic method to quantify the LCN in 3D from synchrotron radiation microtomography images. After segmentation of the LCN, two binary images are generated, one of lacunae (hosting the cell body) and one of canaliculi (small channels linking the lacunae). The binary image of lacunae is labeled, and for each object, lacunar descriptors are extracted after calculating the second order moments and the intrinsic volumes. Furthermore, we propose a specific method to quantify the ramification of canaliculi around each lacuna. To this aim, a signature of the numbers of canaliculi at different distances from the lacunar surface is estimated through the calculation of topological parameters. The proposed method was applied to the 3D SR micro-CT image of a human femoral mid-diaphysis bone sample. Statistical results are reported on 399 lacunae and their surrounding canaliculi.