header advert
Results 1 - 2 of 2
Results per page:
The Bone & Joint Journal
Vol. 99-B, Issue 7 | Pages 927 - 933
1 Jul 2017
Poltaretskyi S Chaoui J Mayya M Hamitouche C Bercik MJ Boileau P Walch G

Aims

Restoring the pre-morbid anatomy of the proximal humerus is a goal of anatomical shoulder arthroplasty, but reliance is placed on the surgeon’s experience and on anatomical estimations. The purpose of this study was to present a novel method, ‘Statistical Shape Modelling’, which accurately predicts the pre-morbid proximal humeral anatomy and calculates the 3D geometric parameters needed to restore normal anatomy in patients with severe degenerative osteoarthritis or a fracture of the proximal humerus.

Materials and Methods

From a database of 57 humeral CT scans 3D humeral reconstructions were manually created. The reconstructions were used to construct a statistical shape model (SSM), which was then tested on a second set of 52 scans. For each humerus in the second set, 3D reconstructions of four diaphyseal segments of varying lengths were created. These reconstructions were chosen to mimic severe osteoarthritis, a fracture of the surgical neck of the humerus and a proximal humeral fracture with diaphyseal extension. The SSM was then applied to the diaphyseal segments to see how well it predicted proximal morphology, using the actual proximal humeral morphology for comparison.


Orthopaedic Proceedings
Vol. 95-B, Issue SUPP_28 | Pages 31 - 31
1 Aug 2013
Mayya M Poltaretskyi S Hamitouche C Chaoui J
Full Access

INTRODUCTION

Automated MRI bone segmentation is one of the most challenging problems in medical imaging. To increase the segmentation robustness, a prior model of the structure could guide the segmentation. Statistical Shape Models (SSMs) are efficient examples for such application. We present an automated SSM construction approach of the scapula bone with an adapted initialisation to address the correspondences problem. Our innovation stems from the derivation of a robust SSM based on Watershed segmentation which steers the elastic registration at some critical zones.

METHODS

The basic idea is to relate only corresponding parts of the shape under investigation. A sample from the samples set is chosen as a common reference (atlas), and the other samples are landmarked and registered to it so that the corresponding points can be identified. The registration has three levels: alignment, rigid and elastic transformations.

To align two scapulae, we define a coordinate system, attach it to each scapula and align both systems. For this, we automatically locate three characteristic points on the scapula's surface. All samples are then scaled to the atlas and the rigid registration is determined by minimising the Euclidian distance between surfaces using Levenberg-Marquadt algorithm.

Afterwards, the samples are locally deformed toward the atlas using directly their landmarks (traditional approach). Unfortunately, landmarks-correspondences could be mismatched at some anatomically complex, “critical,” zones of the scapula. To overcome such a problem, we suggest to 3D-segment these “critical” zones using a 3D Watershed-based method.

Watershed is based on a physical concept of immersion, where it is achieved in a similar way to water filling geographic basins. We believe that this is a natural way to segment the surface of the scapula since it has two large “basins”: the glenoid and the subscapularis fossa. Watershed is followed by geometrical operations to establish eight separated zones on the surface of the scapula.

Once we have the zones, surface-to-surface correspondence is defined and the landmarks' point-to-point correspondences are obtained within each zone pair separately. The elastic registration is then applied on the whole surface via a multi-resolution B-Spline algorithm. The atlas is built through an iterative procedure to eliminate the bias to the initial choice and the correspondences are identified by a reverse registration. Finally, the statistical model can be constructed by performing Principle Component Analysis (PCA).