Development of a simple numerical model for trabecular bone structures

Carlos A. Peña-Solórzano, David W. Albrecht, David M. Paganin, Peter C. Harris, Chris J. Hall, Richard B. Bassed, Matthew R. Dimmock

Research output: Contribution to journalArticleResearchpeer-review

Abstract

Purpose: Advances in additive manufacturing processes are enabling the fabrication of surrogate bone structures for applications including use in high-resolution anthropomorphic phantoms. In this research, a simple numerical model is proposed that enables the generation of microarchitecture with similar statistical distribution to trabecular bone.
Methods: A human humerus, radius, ulna, and several vertebrae were scanned on the Imaging and Medical beamline at the Australian Synchrotron and the proposed numerical model was developed through the definition of two complex functions that encode the trabecular thickness and position-dependant spacing to generate volumetric surrogate trabecular structures. The structures reproduced those observed at 19 separate axial locations through the experimental bone volumes. The applicability of the model when incorporating a two-material approximation to absorption- and phase-contrast CT was also investigated through simulation. Results: The synthetic structures, when compared with the real trabecular microarchitecture, yielded an average mean thickness error of 2 μm, and a mean difference in standard deviation of 33 μm for the humerus, 24 μm for the ulna and radius, and 15 μm for the vertebrae. Simulated absorption- and propagation-based phase contrast CT projection data were generated and reconstructed using the derived mathematical simplifications from the two-material approximation, and the phase-contrast effects were successfully demonstrated.
Conclusions: The presented model reproduced trabecular distributions that could be used to generate phantoms for quality assurance and validation processes. The implication of utilizing a two-material approximation results in simplification of the additive manufacturing process and the generation of synthetic data that could be used for training of machine learning applications.

Original languageEnglish
Article number1766
Number of pages11
JournalMedical Physics
DOIs
Publication statusPublished - 2019

Keywords

  • numerical model
  • PB-CT
  • trabecular

Cite this

@article{a0d69ecbb01e4250829e14c43fe3bd75,
title = "Development of a simple numerical model for trabecular bone structures",
abstract = "Purpose: Advances in additive manufacturing processes are enabling the fabrication of surrogate bone structures for applications including use in high-resolution anthropomorphic phantoms. In this research, a simple numerical model is proposed that enables the generation of microarchitecture with similar statistical distribution to trabecular bone. Methods: A human humerus, radius, ulna, and several vertebrae were scanned on the Imaging and Medical beamline at the Australian Synchrotron and the proposed numerical model was developed through the definition of two complex functions that encode the trabecular thickness and position-dependant spacing to generate volumetric surrogate trabecular structures. The structures reproduced those observed at 19 separate axial locations through the experimental bone volumes. The applicability of the model when incorporating a two-material approximation to absorption- and phase-contrast CT was also investigated through simulation. Results: The synthetic structures, when compared with the real trabecular microarchitecture, yielded an average mean thickness error of 2 μm, and a mean difference in standard deviation of 33 μm for the humerus, 24 μm for the ulna and radius, and 15 μm for the vertebrae. Simulated absorption- and propagation-based phase contrast CT projection data were generated and reconstructed using the derived mathematical simplifications from the two-material approximation, and the phase-contrast effects were successfully demonstrated. Conclusions: The presented model reproduced trabecular distributions that could be used to generate phantoms for quality assurance and validation processes. The implication of utilizing a two-material approximation results in simplification of the additive manufacturing process and the generation of synthetic data that could be used for training of machine learning applications.",
keywords = "numerical model, PB-CT, trabecular",
author = "Pe{\~n}a-Sol{\'o}rzano, {Carlos A.} and Albrecht, {David W.} and Paganin, {David M.} and Harris, {Peter C.} and Hall, {Chris J.} and Bassed, {Richard B.} and Dimmock, {Matthew R.}",
year = "2019",
doi = "10.1002/mp.13435",
language = "English",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "American Association of Physicists in Medicine",

}

Development of a simple numerical model for trabecular bone structures. / Peña-Solórzano, Carlos A.; Albrecht, David W.; Paganin, David M.; Harris, Peter C.; Hall, Chris J.; Bassed, Richard B.; Dimmock, Matthew R.

In: Medical Physics, 2019.

Research output: Contribution to journalArticleResearchpeer-review

TY - JOUR

T1 - Development of a simple numerical model for trabecular bone structures

AU - Peña-Solórzano, Carlos A.

AU - Albrecht, David W.

AU - Paganin, David M.

AU - Harris, Peter C.

AU - Hall, Chris J.

AU - Bassed, Richard B.

AU - Dimmock, Matthew R.

PY - 2019

Y1 - 2019

N2 - Purpose: Advances in additive manufacturing processes are enabling the fabrication of surrogate bone structures for applications including use in high-resolution anthropomorphic phantoms. In this research, a simple numerical model is proposed that enables the generation of microarchitecture with similar statistical distribution to trabecular bone. Methods: A human humerus, radius, ulna, and several vertebrae were scanned on the Imaging and Medical beamline at the Australian Synchrotron and the proposed numerical model was developed through the definition of two complex functions that encode the trabecular thickness and position-dependant spacing to generate volumetric surrogate trabecular structures. The structures reproduced those observed at 19 separate axial locations through the experimental bone volumes. The applicability of the model when incorporating a two-material approximation to absorption- and phase-contrast CT was also investigated through simulation. Results: The synthetic structures, when compared with the real trabecular microarchitecture, yielded an average mean thickness error of 2 μm, and a mean difference in standard deviation of 33 μm for the humerus, 24 μm for the ulna and radius, and 15 μm for the vertebrae. Simulated absorption- and propagation-based phase contrast CT projection data were generated and reconstructed using the derived mathematical simplifications from the two-material approximation, and the phase-contrast effects were successfully demonstrated. Conclusions: The presented model reproduced trabecular distributions that could be used to generate phantoms for quality assurance and validation processes. The implication of utilizing a two-material approximation results in simplification of the additive manufacturing process and the generation of synthetic data that could be used for training of machine learning applications.

AB - Purpose: Advances in additive manufacturing processes are enabling the fabrication of surrogate bone structures for applications including use in high-resolution anthropomorphic phantoms. In this research, a simple numerical model is proposed that enables the generation of microarchitecture with similar statistical distribution to trabecular bone. Methods: A human humerus, radius, ulna, and several vertebrae were scanned on the Imaging and Medical beamline at the Australian Synchrotron and the proposed numerical model was developed through the definition of two complex functions that encode the trabecular thickness and position-dependant spacing to generate volumetric surrogate trabecular structures. The structures reproduced those observed at 19 separate axial locations through the experimental bone volumes. The applicability of the model when incorporating a two-material approximation to absorption- and phase-contrast CT was also investigated through simulation. Results: The synthetic structures, when compared with the real trabecular microarchitecture, yielded an average mean thickness error of 2 μm, and a mean difference in standard deviation of 33 μm for the humerus, 24 μm for the ulna and radius, and 15 μm for the vertebrae. Simulated absorption- and propagation-based phase contrast CT projection data were generated and reconstructed using the derived mathematical simplifications from the two-material approximation, and the phase-contrast effects were successfully demonstrated. Conclusions: The presented model reproduced trabecular distributions that could be used to generate phantoms for quality assurance and validation processes. The implication of utilizing a two-material approximation results in simplification of the additive manufacturing process and the generation of synthetic data that could be used for training of machine learning applications.

KW - numerical model

KW - PB-CT

KW - trabecular

UR - http://www.scopus.com/inward/record.url?scp=85062481623&partnerID=8YFLogxK

U2 - 10.1002/mp.13435

DO - 10.1002/mp.13435

M3 - Article

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

M1 - 1766

ER -