Adaptive Image Segmentation Reveals Substantial Cortical Bone Remodeling During Early Fracture Repair

The goal of this study was to develop an image analysis algorithm for quantifying the effects of remodeling on cortical bone during early fracture healing. An adaptive thresholding technique with boundary curvature and tortuosity control was developed to automatically identify the endocortical and pericortical boundaries in the presence of high-gradient bone mineral density (BMD) near the healing zone. The algorithm successfully segmented more than 47,000 microCT images from 12 healing ovine osteotomies and intact contralateral tibiae. Resampling techniques were used to achieve data dimensionality reduction on the segmented images, allowing characterization of radial and axial distributions of cortical BMD. Local (transverse slice) and total (whole bone) remodeling scores were produced. These surrogate measures of cortical remodeling derived from BMD revealed that cortical changes were detectable throughout the region covered by callus and that the localized loss of cortical BMD was highest near the osteotomy. Total remodeling score was moderately and significantly correlated with callus volume and mineral composition (r > 0.64, p < 0.05), suggesting that the cortex may be a source of mineral needed to build callus.