Search

• Shape registration was used to map bone shapes onto template bone, allowing the spatial correlation structure of bone mineral density to be studied. • Non-stationary random field of bone mineral density with deterministic age dependence and sex sensitive was explored. • Energy distance between Gaussian and non-Gaussian random fields is convergent, but to non-zero value. • Average bone mineral density can be simulated with only one coefficient in Karhunen-Loève expansion with acceptable accuracy. • Interactive web application BoneGen allows to generate pelvic bone mineral density realisations. Capturing the population variability of bone properties is of paramount importance to biomedical engineering. The aim of the present paper is to describe variability and correlations in bone mineral density with a spatial random field inferred from routine computed tomography data. Random fields were simulated by transforming pairwise uncorrelated Gaussian random variables into correlated variables through the spectral decomposition of an age-detrended correlation matrix. The validity of the random field model was demonstrated in the spatiotemporal analysis of bone mineral density. The similarity between the computed tomography samples and those generated via random fields was analyzed with the energy distance metric. The random field of bone mineral density was found to be approximately Gaussian/slightly left-skewed/strongly right-skewed at various locations. However, average bone density could be simulated well with the proposed Gaussian random field for which the energy distance, i.e., a measure that quantifies discrepancies between two distribution functions, is convergent with respect to the number of correlation eigenpairs. The proposed random field model allows the enhancement of computational biomechanical models with variability in bone mineral density, which could increase the usability of the model and provides a step forward in in-silico medicine. [Display omitted] [ABSTRACT FROM AUTHOR]