Three‐dimensional multifractal analysis of trabecular bone under clinical computed tomography

摘要

Purpose An adequate understanding of bone structural properties is critical for predicting fragility conditions caused by diseases such as osteoporosis, and in gauging the success of fracture prevention treatments. In this work we aim to develop multiresolution image analysis techniques to extrapolate high‐resolution images predictive power to images taken in clinical conditions. Methods We performed multifractal analysis (MFA) on a set of 17 ex?vivo human vertebrae clinical CT scans. The vertebr? failure loads ( F Failure ) were experimentally measured. We combined bone mineral density (BMD) with different multifractal dimensions, and BMD with multiresolution statistics (e.g., skewness, kurtosis) of MFA curves, to obtain linear models to predict F Failure . Furthermore we obtained short‐ and long‐term precisions from simulated in?vivo scans, using a clinical CT scanner. Ground‐truth data — high‐resolution images — were obtained with a High‐Resolution Peripheral Quantitative Computed Tomography (HRpQCT) scanner. Results At the same level of detail, BMD combined with traditional multifractal descriptors (Lipschitz–H?lder exponents), and BMD with monofractal features showed similar prediction powers in predicting F Failure (87%, adj. R 2 ). However, at different levels of details, the prediction power of BMD with multifractal features raises to 92% (adj. R 2 ) of F Failure . Our main finding is that a simpler but slightly less accurate model, combining BMD and the skewness of the resulting multifractal curves, predicts 90% (adj. R 2 ) of F Failure . Conclusions Compared to monofractal and standard bone measures, multifractal analysis captured key insights in the conditions leading to F Failure . Instead of raw multifractal descriptors, the statistics of multifractal curves can be used in several other contexts, facilitating further research.