Department of Radiology and Diagnostic Imaging, University of Alberta, Edmonton, AB, T6G 2B7, Canada.
Department of Physics, University of Alberta, Edmonton, AB, T6G 2E1, Canada.
Ann Biomed Eng. 2019 Nov;47(11):2178-2187. doi: 10.1007/s10439-019-02310-4. Epub 2019 Jun 19.
In this study, a nonlinear grid-search inversion has been developed to estimate the thickness and elastic velocities of long cortical bones, which are important determinants of bone strength, from axially-transmitted ultrasonic data. The inversion scheme is formulated in the dispersive frequency-phase velocity domain to recover bone properties. The method uses ultrasonic guided waves to retrieve overlying soft tissue thickness, cortical thickness, compressional, and shear-wave velocities of the cortex. The inversion strategy requires systematic examination of a large set of trial dispersion-curve solutions within a pre-defined model space to match the data with minimum cost in a least-squares sense. The theoretical dispersion curves required to solve the inverse problem are computed for bilayered bone models using a semi-analytical finite-element method. The feasibility of the proposed approach was demonstrated by the numerically simulated data for a 1 mm soft tissue-5 mm bone bilayer and ex-vivo data from a bovine femur plate with an overlying 2 mm-thick soft-tissue mimic. The bootstrap method was employed to evaluate the inversion uncertainty and stability. Our results have shown that the cortical thickness and wave speeds could be recovered with fair accuracy.
在这项研究中,开发了一种非线性网格搜索反演方法,用于从轴向透射超声数据估计长皮质骨的厚度和弹性速度,这些是骨骼强度的重要决定因素。反演方案在频散相位速度域中制定,以恢复骨骼特性。该方法使用超声导波来获取覆盖的软组织厚度、皮质厚度、皮质的压缩波和剪切波速度。反演策略需要系统地检查一组预定义模型空间内的大量试验色散曲线解,以便在最小二乘意义上使数据与最小代价相匹配。使用半解析有限元方法为双层骨模型计算求解反问题所需的理论色散曲线。通过对 1 毫米软组织-5 毫米骨双层的数值模拟数据和具有 2 毫米厚软组织模拟的牛股骨板的离体数据进行的数值模拟,验证了所提出方法的可行性。使用自举法评估反演不确定性和稳定性。我们的结果表明,皮质厚度和波速可以以相当高的精度恢复。