3D reconstruction of the scapula from biplanar radiographs

Access to 3D bone models is critical for applications ranging from pre-operative planning to biomechanics studies. This work presents a method for 3D reconstruction of the scapula from biplanar radiographs, which is based on the combination of a parametric model approach in conjunction with a Moving Least Squares (MLS) deformation technique. A parametric scapula model was created by fitting geometric primitives (with their descriptive parameters) to the CT reconstruction of a dry scapula. These geometric primitives were then used to define a set of handles which allow the user to control the as-rigid-as-possible deformation of the template model in real-time, until optimal correspondence between the actual X-ray images and the retro-projection of the deformed model. When applied to 10 dry scapulae, the presented method allowed obtaining reconstructions which were on average within 1mm of the CT-derived model at scapula regions of interest. Morphological parameters such as the glenoid's dimensions and orientation were determined with errors of 1° and less than 1mm, on average. This is of great interest as the current methods used in clinical practice, which are based on 2D-CT, are subject to uncertainties of the order of 5° for glenoid version. This method is of particular interest as it further reduces our dependence to CT for 3D reconstruction of bones and clinical parameter estimation.

[1]  Scott Tashman,et al.  Validation of a new model-based tracking technique for measuring three-dimensional, in vivo glenohumeral joint kinematics. , 2006, Journal of biomechanical engineering.

[2]  Paul J. Besl,et al.  A Method for Registration of 3-D Shapes , 1992, IEEE Trans. Pattern Anal. Mach. Intell..

[3]  Weiliang Xu,et al.  Diagnostic radiograph based 3D bone reconstruction framework: Application to the femur , 2011, Comput. Medical Imaging Graph..

[4]  W. Skalli,et al.  Coupling 2D/3D registration method and statistical model to perform 3D reconstruction from partial x-rays images data , 2009, 2009 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[5]  Kimerly A Powell,et al.  Use of three-dimensional computed tomography for the analysis of the glenoid anatomy. , 2005, Journal of shoulder and elbow surgery.

[6]  G. Hazan,et al.  Variability of measurement of glenoid version on computed tomography scan. , 1999, Journal of shoulder and elbow surgery.

[7]  Gilles Walch,et al.  Glenoid version: how to measure it? Validity of different methods in two-dimensional computed tomography scans. , 2010, Journal of shoulder and elbow surgery.

[8]  Nicola Hagemeister,et al.  3D morphometric analysis of 43 scapulae , 2012, Surgical and Radiologic Anatomy.

[9]  Scott Schaefer,et al.  Image deformation using moving least squares , 2006, ACM Trans. Graph..

[10]  J.A. De Guise,et al.  Surface reconstruction from planar x-ray images using moving least squares , 2008, 2008 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[11]  E. Berthonnaud,et al.  Three-dimensional in vivo displacements of the shoulder complex from biplanar radiography , 2005, Surgical and Radiologic Anatomy.

[12]  K. S. Arun,et al.  Least-Squares Fitting of Two 3-D Point Sets , 1987, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[13]  Guoyan Zheng,et al.  A 2D/3D correspondence building method for reconstruction of a patient-specific 3D bone surface model using point distribution models and calibrated X-ray images , 2009, Medical Image Anal..