Generative Modelling of 3D in-silico Spongiosa with Controllable Micro-Structural Parameters

Research in vertebral bone micro-structure generally requires costly procedures to obtain physical scans of real bone with a specific pathology under study, since no methods are available yet to generate realistic bone structures in-silico. Here we propose to apply recent advances in generative adversarial networks (GANs) to develop such a method. We adapted style-transfer techniques, which have been largely used in other contexts, in order to transfer style between image pairs while preserving its informational content. In a first step, we trained a volumetric generative model in a progressive manner using a Wasserstein objective and gradient penalty (PWGAN-GP) to create patches of realistic bone structure in-silico. The training set contained 7660 purely spongeous bone samples from twelve human vertebrae (T12 or L1) with isotropic resolution of 164um and scanned with a high resolution peripheral quantitative CT (Scanco XCT). After training, we generated new samples with tailored micro-structure properties by optimizing a vector z in the learned latent space. To solve this optimization problem, we formulated a differentiable goal function that leads to valid samples while compromising the appearance (content) with target 3D properties (style). Properties of the learned latent space effectively matched the data distribution. Furthermore, we were able to simulate the resulting bone structure after deterioration or treatment effects of osteoporosis therapies based only on expected changes of micro-structural parameters. Our method allows to generate a virtually infinite number of patches of realistic bone micro-structure, and thereby likely serves for the development of bone-biomarkers and to simulate bone therapies in advance.

[1]  Carl-Fredrik Westin,et al.  Visualization and Processing of Tensors and Higher Order Descriptors for Multi-Valued Data , 2014 .

[2]  P. Zysset,et al.  Comparative Effects of Teriparatide and Risedronate in Glucocorticoid-Induced Osteoporosis in Men: 18-Month Results of the EuroGIOPs Trial , 2013, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[3]  Örjan Smedby,et al.  Techniques for Computing Fabric Tensors: A Review , 2014, Visualization and Processing of Tensors and Higher Order Descriptors for Multi-Valued Data.

[5]  Simon Osindero,et al.  Conditional Generative Adversarial Nets , 2014, ArXiv.

[6]  Keigo Hirakawa,et al.  Fast Convolutional Distance Transform , 2019, IEEE Signal Processing Letters.

[7]  Guigang Zhang,et al.  Deep Learning , 2016, Int. J. Semantic Comput..

[8]  Jimmy Ba,et al.  Adam: A Method for Stochastic Optimization , 2014, ICLR.

[9]  Alena-Kathrin Schnurr,et al.  Synthesis of CT images from digital body phantoms using CycleGAN , 2019, International Journal of Computer Assisted Radiology and Surgery.

[10]  Eli Konen,et al.  Creating Artificial Images for Radiology Applications Using Generative Adversarial Networks (GANs) - A Systematic Review. , 2020, Academic radiology.

[11]  Leon A. Gatys,et al.  A Neural Algorithm of Artistic Style , 2015, ArXiv.

[12]  Stefan Winkler,et al.  The Unusual Effectiveness of Averaging in GAN Training , 2018, ICLR.

[13]  Örjan Smedby,et al.  Evaluation of the plate-rod model assumption of trabecular bone , 2012, 2012 9th IEEE International Symposium on Biomedical Imaging (ISBI).

[14]  G. Huber,et al.  A new algorithm for estimating the rod volume fraction and the trabecular thickness from in vivo computed tomography. , 2016, Medical physics.

[15]  Léon Bottou,et al.  Wasserstein GAN , 2017, ArXiv.

[16]  Hiroshi Fujita,et al.  Automated Pulmonary Nodule Classification in Computed Tomography Images Using a Deep Convolutional Neural Network Trained by Generative Adversarial Networks , 2019, BioMed research international.

[17]  Aaron C. Courville,et al.  Improved Training of Wasserstein GANs , 2017, NIPS.

[18]  Jaakko Lehtinen,et al.  Progressive Growing of GANs for Improved Quality, Stability, and Variation , 2017, ICLR.

[19]  Hayit Greenspan,et al.  GAN-based Synthetic Medical Image Augmentation for increased CNN Performance in Liver Lesion Classification , 2018, Neurocomputing.