Can we improve the prediction of hip fracture by assessing bone structure using shape and appearance modelling?

Copyright 2013 Elsevier B.V., All rights reserved.Peer reviewedPreprin

Image analysis for extracapsular hip fracture surgery

PhD ThesisDuring the implant insertion phase of extracapsular hip fracture surgery, a surgeon visually inspects digital radiographs to infer the best position for the implant. The inference is made by “eye-balling”. This clearly leaves room for trial and error which is not ideal for the patient. This thesis presents an image analysis approach to estimating the ideal positioning for the implant using a variant of the deformable templates model known as the Constrained Local Model (CLM). The Model is a synthesis of shape and local appearance models learned from a set of annotated landmarks and their corresponding local patches extracted from digital femur x-rays. The CLM in this work highlights both Principal Component Analysis (PCA) and Probabilistic PCA as regularisation components; the PPCA variant being a novel adaptation of the CLM framework that accounts for landmark annotation error which the PCA version does not account for. Our CLM implementation is used to articulate 2 clinical metrics namely: the Tip-Apex Distance and Parker’s Ratio (routinely used by clinicians to assess the positioning of the surgical implant during hip fracture surgery) within the image analysis framework. With our model, we were able to automatically localise signi cant landmarks on the femur, which were subsequently used to measure Parker’s Ratio directly from digital radiographs and determine an optimal placement for the surgical implant in 87% of the instances; thereby, achieving fully automatic measurement of Parker’s Ratio as opposed to manual measurements currently performed in the surgical theatre during hip fracture surgery

Computer-aided diagnosis of complications of total hip replacement X-ray images

Hip replacement surgery has experienced a dramatic evolution in recent years supported by the latest developments in many areas of technology and surgical procedures. Unfortunately complications that follow hip replacement surgery remains the most challenging dilemma faced both by the patients and medical experts. The thesis presents a novel approach to segment the prosthesis of a THR surgical process by using an Active Contour Model (ACM) that is initiated via an automatically detected seed point within the enarthrosis region of the prosthesis. The circular area is detected via the use of a Fast, Randomized Circle Detection Algorithm. Experimental results are provided to compare the performance of the proposed ACM based approach to popular thresholding based approaches. Further an approach to automatically detect the Obturator Foramen using an ACM approach is also presented. Based on analysis of how medical experts carry out the detection of loosening and subsidence of a prosthesis and the presence of infections around the prosthesis area, this thesis presents novel computational analysis concepts to identify the key feature points of the prosthesis that are required to detect all of the above three types of complications. Initially key points along the prosthesis boundary are determined by measuring the curvature on the surface of the prosthesis. By traversing the edge pixels, starting from one end of the boundary of a detected prosthesis, the curvature values are determined and effectively used to determine key points of the prosthesis surface and their relative positioning. After the key-points are detected, pixel value gradients across the boundary of the prosthesis are determined along the boundary of the prosthesis to determine the presence of subsidence, loosening and infections. Experimental results and analysis are presented to show that the presence of subsidence is determined by the identification of dark pixels around the convex bend closest to the stem area of the prosthesis and away from it. The presence of loosening is determined by the additional presence of dark regions just outside the two straight line edges of the stem area of the prosthesis. The presence of infections is represented by the determination of dark areas around the tip of the stem of the prosthesis. All three complications are thus determined by a single process where the detailed analysis defer. The experimental results presented show the effectiveness of all proposed approaches which are also compared and validated against the ground truth recorded manually with expert user input

A Hierarchical Method Based on Active Shape Models and Directed Hough Transform for Segmentation of Noisy Biomedical Images; Application in Segmentation of Pelvic X-ray Images

Background Traumatic pelvic injuries are often associated with severe, life-threatening hemorrhage, and immediate medical treatment is therefore vital. However, patient prognosis depends heavily on the type, location and severity of the bone fracture, and the complexity of the pelvic structure presents diagnostic challenges. Automated fracture detection from initial patient X-ray images can assist physicians in rapid diagnosis and treatment, and a first and crucial step of such a method is to segment key bone structures within the pelvis; these structures can then be analyzed for specific fracture characteristics. Active Shape Model has been applied for this task in other bone structures but requires manual initialization by the user. This paper describes a algorithm for automatic initialization and segmentation of key pelvic structures - the iliac crests, pelvic ring, left and right pubis and femurs - using a hierarchical approach that combines directed Hough transform and Active Shape Models. Results Performance of the automated algorithm is compared with results obtained via manual initialization. An error measures is calculated based on the shapes detected with each method and the gold standard shapes. ANOVA results on these error measures show that the automated algorithm performs at least as well as the manual method. Visual inspection by two radiologists and one trauma surgeon also indicates generally accurate performance. Conclusion The hierarchical algorithm described in this paper automatically detects and segments key structures from pelvic X-rays. Unlike various other x-ray segmentation methods, it does not require manual initialization or input. Moreover, it handles the inconsistencies between x-ray images in a clinical environment and performs successfully in the presence of fracture. This method and the segmentation results provide a valuable base for future work in fracture detection

Development of a fully automatic shape model matching (FASMM) system to derive statistical shape models from radiographs: application to the accurate capture and global representation of proximal femur shape

SummaryObjectiveTo evaluate the accuracy and sensitivity of a fully automatic shape model matching (FASMM) system to derive statistical shape models (SSMs) of the proximal femur from non-standardised anteroposterior (AP) pelvic radiographs.DesignAP pelvic radiographs obtained with informed consent and appropriate ethical approval were available for 1105 subjects with unilateral hip osteoarthritis (OA) who had been recruited previously for The arcOGEN Study. The FASMM system was applied to capture the shape of the unaffected (i.e., without signs of radiographic OA) proximal femur from these radiographs. The accuracy and sensitivity of the FASMM system in calculating geometric measurements of the proximal femur and in shape representation were evaluated relative to validated manual methods.ResultsDe novo application of the FASMM system had a mean point-to-curve error of less than 0.9 mm in 99% of images (n = 266). Geometric measurements generated by the FASMM system were as accurate as those obtained manually. The analysis of the SSMs generated by the FASMM system for male and female subject groups identified more significant differences (in five of 17 SSM modes after Bonferroni adjustment) in their global proximal femur shape than those obtained from the analysis of conventional geometric measurements. Multivariate gender-classification accuracy was higher when using SSM mode values (76.3%) than when using conventional hip geometric measurements (71.8%).ConclusionsThe FASMM system rapidly and accurately generates a global SSM of the proximal femur from radiographs of varying quality and resolution. This system will facilitate complex morphometric analysis of global shape variation across large datasets. The FASMM system could be adapted to generate SSMs from the radiographs of other skeletal structures such as the hand, knee or pelvis

Active Shape Modeling of the Hip in the Prediction of Incident Hip Fracture

The objective of this study was to evaluate right proximal femur shape as a risk factor for incident hip fracture using active shape modeling (ASM). A nested case-control study of white women 65 years of age and older enrolled in the Study of Osteoporotic Fractures (SOF) was performed. Subjects (n = 168) were randomly selected from study participants who experienced hip fracture during the follow-up period (mean 8.3 years). Controls (n = 231) had no fracture during follow-up. Subjects with baseline radiographic hip osteoarthritis were excluded. ASM of digitized right hip radiographs generated 10 independent modes of variation in proximal femur shape that together accounted for 95% of the variance in proximal femur shape. The association of ASM modes with incident hip fracture was analyzed by logistic regression. Together, the 10 ASM modes demonstrated good discrimination of incident hip fracture. In models controlling for age and body mass index (BMI), the area under receiver operating characteristic (AUROC) curve for hip shape was 0.813, 95% confidence interval (CI) 0.771–0.854 compared with models containing femoral neck bone mineral density (AUROC = 0.675, 95% CI 0.620–0.730), intertrochanteric bone mineral density (AUROC = 0.645, 95% CI 0.589–0.701), femoral neck length (AUROC = 0.631, 95% CI 0.573–0.690), or femoral neck width (AUROC = 0.633, 95% CI 0.574–0.691). The accuracy of fracture discrimination was improved by combining ASM modes with femoral neck bone mineral density (AUROC = 0.835, 95% CI 0.795–0.875) or with intertrochanteric bone mineral density (AUROC = 0.834, 95% CI 0.794–0.875). Hips with positive standard deviations of ASM mode 4 had the highest risk of incident hip fracture (odds ratio = 2.48, 95% CI 1.68–3.31, p < .001). We conclude that variations in the relative size of the femoral head and neck are important determinants of incident hip fracture. The addition of hip shape to fracture-prediction tools may improve the risk assessment for osteoporotic hip fractures. © 2011 American Society for Bone and Mineral Research