Osteocytes are the most abundant cell in the bone, and have multiple functions including mechanosensing and regulation of bone remodeling activities. Since osteocytes are embedded in the bone matrix, their inaccessibility makes in vivo studies problematic. Therefore, a non-invasive technique with high spatial resolution is desired. The purpose of this study is to investigate the use of third harmonic generation (THG) microscopy as a noninvasive technique for high-resolution imaging of the lacunar-canalicular network (LCN) in live mice. By performing THG imaging in combination with two-and three-photon fluorescence microscopy, we show that THG signal is produced from the bone-interstitial fluid boundary of the lacuna, while the interstitial fluid-osteocyte cell boundary shows a weaker THG signal. Canaliculi are also readily visualized by THG imaging, with canaliculi oriented at small angles relative to the optical axis exhibiting stronger signal intensity compared to those oriented perpendicular to the optical axis (parallel to the image plane). By measuring forward-versus epi-detected THG signals in thinned versus thick bone samples ex vivo, we found that the epi-collected THG from the LCN of intact bone contains a superposition of backward-directed and backscattered forward-THG. As an example of a biological application, THG was used as a label-free imaging technique to study structural variations in the LCN of live mice deficient in both histone deacetylase 4 and 5 (HDAC4, HDAC5). Three-dimensional analyses were performed and revealed statistically significant differences between the HDAC4/5 double knockout and wild type mice in the number of osteocytes per volume and the number of canaliculi per lacunar surface area. These changes in osteocyte density and dendritic projections occurred without differences in lacunar size. This study demonstrates that THG microscopy imaging of the LCN in live mice enables quantitative analysis of osteocytes in animal models without the use of dyes or physical sectioning.