Karl D. Fritscher

3D image segmentation using combined shape-intensity prior models

ObjectiveStatistical models for medical images have been developed to increase robustness in the segmentation process. In this project, a fully automatic approach to build a statistical shape-intensity model and combine this model with level set segmentation was designed, implemented and tested by applying the algorithm to clinical image data.MethodsBy using a hierarchical registration approach based on mutual information and demons registration, 3D statistical shape-intensity models were created by applying Principal Component Analysis. Using these models in combination with level set segmentation results in a fully automatic modeling and segmentation pipeline.ResultsExamples for shape-intensity models were synthesized and these models were used to automatically segment 3D MRI and CT images. Quantitative evaluation of the framework was performed by comparing automatic segmentation results to segmentation results of medical experts.ConclusionEvaluation tests in which this method was used for the automatic segmentation of femora and cardiac MRI endocardial surfaces are very promising. The implementation of an additional cost function term and the addition of information about the surroundings of an organ in the model are currently under development.

Automatic segmentation of head and neck CT images for radiotherapy treatment planning using multiple atlases, statistical appearance models, and geodesic active contours.

PURPOSE Accurate delineation of organs at risk (OARs) is a precondition for intensity modulated radiation therapy. However, manual delineation of OARs is time consuming and prone to high interobserver variability. Because of image artifacts and low image contrast between different structures, however, the number of available approaches for autosegmentation of structures in the head-neck area is still rather low. In this project, a new approach for automated segmentation of head-neck CT images that combine the robustness of multiatlas-based segmentation with the flexibility of geodesic active contours and the prior knowledge provided by statistical appearance models is presented. METHODS The presented approach is using an atlas-based segmentation approach in combination with label fusion in order to initialize a segmentation pipeline that is based on using statistical appearance models and geodesic active contours. An anatomically correct approximation of the segmentation result provided by atlas-based segmentation acts as a starting point for an iterative refinement of this approximation. The final segmentation result is based on using model to image registration and geodesic active contours, which are mutually influencing each other. RESULTS 18 CT images in combination with manually segmented labels of parotid glands and brainstem were used in a leave-one-out cross validation scheme in order to evaluate the presented approach. For this purpose, 50 different statistical appearance models have been created and used for segmentation. Dice coefficient (DC), mean absolute distance and max. Hausdorff distance between the autosegmentation results and expert segmentations were calculated. An average Dice coefficient of DC = 0.81 (right parotid gland), DC = 0.84 (left parotid gland), and DC = 0.86 (brainstem) could be achieved. CONCLUSIONS The presented framework provides accurate segmentation results for three important structures in the head neck area. Compared to a segmentation approach based on using multiple atlases in combination with label fusion, the proposed hybrid approach provided more accurate results within a clinically acceptable amount of time.

Interactive Medical Image Segmentation Using PDE Control of Active Contours

Segmentation of injured or unusual anatomic structures in medical imagery is a problem that has continued to elude fully automated solutions. In this paper, the goal of easy-to-use and consistent interactive segmentation is transformed into a control synthesis problem. A nominal level set partial differential equation (PDE) is assumed to be given; this open-loop system achieves correct segmentation under ideal conditions, but does not agree with a human experts ideal boundary for real image data. Perturbing the state and dynamics of a level set PDE via the accumulated user input and an observer-like system leads to desirable closed-loop behavior. The input structure is designed such that a user can stabilize the boundary in some desired state without needing to understand any mathematical parameters. Effectiveness of the technique is illustrated with applications to the challenging segmentations of a patellar tendon in magnetic resonance and a shattered femur in computed tomography.

Hip fracture discrimination from dual-energy X-ray absorptiometry by statistical model registration

Although the areal Bone Mineral Density (BMD) measurements from dual-energy X-ray absorptiometry (DXA) are able to discriminate between hip fracture cases and controls, the femoral strength is largely determined by the 3D bone structure. In a previous work a statistical model was presented which parameterizes the 3D shape and BMD distribution of the proximal femur. In this study the parameter values resulting from the registration of the model onto DXA images are evaluated for their hip fracture discrimination ability with respect to regular DXA derived areal BMD measurements. The statistical model was constructed from a large database of QCT scans of females with an average age of 67.8 ± 17.0 years. This model was subsequently registered onto the DXA images of a fracture and control group. The fracture group consisted of 175 female patients with an average age of 66.4 ± 9.9 years who suffered a fracture on the contra lateral femur. The control group consisted of 175 female subjects with an average age of 65.3 ± 10.0 years and no fracture history. The discrimination ability of the resulting model parameter values, as well as the areal BMD measurements extracted from the DXA images were evaluated using a logistic regression analysis. The area under the receiver operating curve (AUC) of the combined model parameters and areal BMD values was 0.840 (95% CI 0.799-0.881), whilst using only the areal BMD values resulted in an AUC of 0.802 (95% CI 0.757-0.848). These results indicate that the discrimination ability of the areal BMD values is improved by supplementing them with the model parameter values, which give a more complete representation of the subject specific shape and internal bone distribution. Thus, the presented method potentially allows for an improved hip fracture risk estimation whilst maintaining DXA as the current standard modality.

Trabecular Bone Analysis in CT and X-Ray Images of the Proximal Femur for the Assessment of Local Bone Quality

Currently, conventional X-ray and CT images as well as invasive methods performed during the surgical intervention are used to judge the local quality of a fractured proximal femur. However, these approaches are either dependent on the surgeons experience or cannot assist diagnostic and planning tasks preoperatively. Therefore, in this work a method for the individual analysis of local bone quality in the proximal femur based on model-based analysis of CT- and X-ray images of femur specimen will be proposed. A combined representation of shape and spatial intensity distribution of an object and different statistical approaches for dimensionality reduction are used to create a statistical appearance model in order to assess the local bone quality in CT and X-ray images. The developed algorithms are tested and evaluated on 28 femur specimen. It will be shown that the tools and algorithms presented herein are highly adequate to automatically and objectively predict bone mineral density values as well as a biomechanical parameter of the bone that can be measured intraoperatively.

Automatic cardiac 4d segmentation using level sets

For the analysis of shape variations of the heart and the cardiac motion in a clinical environment it is necessary to segment a large amount of data in order to be able to build statistically significant models. Therefore it has been the aim of this project to find and develop methods that allow the creation of a fully automatic segmentation pipeline for the segmentation of endocardium and myocardium in ECG-triggered MRI images. For this purpose a combination of a number of image processing techniques, from the fields of segmentation, modeling and image registration have been used and extended to create a segmentation pipeline that reduces the need for supplementary manual correction of the segmented labels to a minimum.

A statistical model of shape and bone mineral density distribution of the proximal femur for fracture risk assessment

This work presents a statistical model of both the shape and Bone Mineral Density (BMD) distribution of the proximal femur for fracture risk assessment. The shape and density model was built from a dataset of Quantitative Computed Tomography scans of fracture patients and a control group. Principal Component Analysis and Horns parallel analysis were used to reduce the dimensionality of the shape and density model to the main modes of variation. The input data was then used to analyze the model parameters for the optimal separation between the fracture and control group. Feature selection using the Fisher criterion determined the parameters with the best class separation, which were used in Fisher Linear Discriminant Analysis to find the direction in the parameter space that best separates the fracture and control group. This resulted in a Fisher criterion value of 6.70, while analyzing the Dual-energy X-ray Absorptiometry derived femur neck areal BMD of the same subjects resulted in a Fisher criterion value of 0.98. This indicates that a fracture risk estimation approach based on the presented model might improve upon the current standard clinical practice.

Assessment of the individual fracture risk of the proximal femur by using statistical appearance models.

PURPOSE Standard diagnostic techniques to quantify bone mineral density (BMD) include dual-energy x-ray absorptiometry (DXA) and quantitative computed tomography. However, BMD alone is not sufficient to predict the fracture risk for an individual patient. Therefore, the development of tools, which can assess the bone quality in order to predict individual biomechanics of a bone, would mean a significant improvement for the prevention of fragility fractures. In this study, a new approach to predict the fracture risk of proximal femora using a statistical appearance model will be presented. METHODS 100 CT data sets of human femur cadaver specimens are used to create statistical appearance models for the prediction of the individual fracture load (FL). Calculating these models offers the possibility to use information about the inner structure of the proximal femur, as well as geometric properties of the femoral bone for FL prediction. By applying principal component analysis, statistical models have been calculated in different regions of interest. For each of these models, the individual model parameters for each single data set were calculated and used as predictor variables in a multilinear regression model. By this means, the best working region of interest for the prediction of FL was identified. The accuracy of the FL prediction was evaluated by using a leave-one-out cross validation scheme. Performance of DXA in predicting FL was used as a standard of comparison. RESULTS The results of the evaluative tests demonstrate that significantly better results for FL prediction can be achieved by using the proposed model-based approach (R=0.91) than using DXA-BMD (R=0.81) for the prediction of fracture load. CONCLUSIONS The results of the evaluation show that the presented model-based approach is very promising and also comparable to studies that partly used higher image resolutions for bone quality assessment and fracture risk prediction.

Evaluation of segmentation methods on head and neck CT: Auto‐segmentation challenge 2015

Purpose Automated delineation of structures and organs is a key step in medical imaging. However, due to the large number and diversity of structures and the large variety of segmentation algorithms, a consensus is lacking as to which automated segmentation method works best for certain applications. Segmentation challenges are a good approach for unbiased evaluation and comparison of segmentation algorithms. Methods In this work, we describe and present the results of the Head and Neck Auto‐Segmentation Challenge 2015, a satellite event at the Medical Image Computing and Computer Assisted Interventions (MICCAI) 2015 conference. Six teams participated in a challenge to segment nine structures in the head and neck region of CT images: brainstem, mandible, chiasm, bilateral optic nerves, bilateral parotid glands, and bilateral submandibular glands. Results This paper presents the quantitative results of this challenge using multiple established error metrics and a well‐defined ranking system. The strengths and weaknesses of the different auto‐segmentation approaches are analyzed and discussed. Conclusions The Head and Neck Auto‐Segmentation Challenge 2015 was a good opportunity to assess the current state‐of‐the‐art in segmentation of organs at risk for radiotherapy treatment. Participating teams had the possibility to compare their approaches to other methods under unbiased and standardized circumstances. The results demonstrate a clear tendency toward more general purpose and fewer structure‐specific segmentation algorithms.

Contour-Driven Atlas-Based Segmentation

We propose new methods for automatic segmentation of images based on an atlas of manually labeled scans and contours in the image. First, we introduce a Bayesian framework for creating initial label maps from manually annotated training images. Within this framework, we model various registration- and patch-based segmentation techniques by changing the deformation field prior. Second, we perform contour-driven regression on the created label maps to refine the segmentation. Image contours and image parcellations give rise to non-stationary kernel functions that model the relationship between image locations. Setting the kernel to the covariance function in a Gaussian process establishes a distribution over label maps supported by image structures. Maximum a posteriori estimation of the distribution over label maps conditioned on the outcome of the atlas-based segmentation yields the refined segmentation. We evaluate the segmentation in two clinical applications: the segmentation of parotid glands in head and neck CT scans and the segmentation of the left atrium in cardiac MR angiography images.