Search

We propose a 3D+t/2D+t registration strategy to relate preoperative CTA to intraoperative X-ray angiography for improved image guidance during coronary interventions. We first derive 4D coronary models from CTA and then align these models both temporally and spatially to the intraoperative images. Temporal alignment is based on aligning ECG signals, whereas the spatial alignment uses a similarity metric based on centerline projection and fuzzy X-ray segmentation. In the spatial alignment step we use information from multiple time points simultaneously and take into account rigid respiratory motion. Evaluation shows improved performance compared to a 3D/2D approach with respect to both registration success and reproducibility.

Despite the growing interest in regression based shape estimation, no study has yet systematically compared different regression methods for shape estimation. We aimed to fill this gap by comparing linear regression methods with a special focus on shapes with landmark position uncertainties. We investigate two scenarios: In the first, the uncertainty of the landmark positions was similar in the training and test dataset, whereas in the second the uncertainty of the training and test data were different. Both scenarios were tested on simulated data and on statistical models of the left ventricle estimating the end-systolic shape from end-diastole with landmark uncertainties derived from the segmentation process, and of the femur estimating the 3D shape from one projection with landmark uncertainties derived from the imaging setup. Results show that in the first scenario linear regression methods tend to perform similar. In the second scenario including estimates of the test shape landmark uncertainty in the regression improved results.

We propose a conditional statistical shape model to predict patient specific cardiac motion from the 3D end-diastolic CTA scan. The model is built from 4D CTA sequences by combining atlas based segmentation and 4D registration. Cardiac motion estimation is, for example, relevant in the dynamic alignment of pre-operative CTA data with intra-operative X-ray imaging. Due to a trend towards prospective electrocardiogram gating techniques, 4D imaging data, from which motion information could be extracted, is not commonly available. The prediction of motion from shape information is thus relevant for this purpose. Evaluation of the accuracy of the predicted motion was performed using CTA scans of 50 patients, showing an average accuracy of 1.1 mm.

Quantitative information about the geometry of the carotid artery bifurcation may help in predicting the development of atherosclerosis. A geodesic active contours based segmentation method combining both gradient and intensity information was developed for semi-automatic, accurate and robust quantification of the carotid bifurcation angle in Black Blood MRA data. The segmentation method was evaluated by comparing its accuracy to inter and intra observer variability on a large dataset that has been acquired as part of a longitudinal population study which investigates the natural progression of carotid atherosclerosis. Furthermore, the method is shown to be robust to initialization differences. The bifurcation angle obtained from the segmented lumen corresponds well with the angle derived from the manual lumen segmentation, which demonstrates that the method has large potential to replace manual segmentations for extracting the carotid bifurcation angle from Black Blood MRA data.

We present a framework for statistical analysis in large cohorts of structural brain connectivity, derived from diffusion weighted MRI. A brain network is defined between subcortical gray matter structures and a cortical parcellation obtained with FreeSurfer. Connectivity is established through minimum cost paths with an anisotropic local cost function and is quantified per connection. The connectivity network potentially encodes important information about brain structure, and can be analyzed using multivariate regression methods. The proposed framework can be used to study the relation between connectivity and e.g. brain function or neurodegenerative disease. As a proof of principle, we perform principal component regression in order to predict age and gender, based on the connectivity networks of 979 middle-aged and elderly subjects, in a 10-fold cross-validation. The results are compared to predictions based on fractional anisotropy and mean diffusivity averaged over the white matter and over the corpus callosum. Additionally, the predictions are performed based on the best predicting connection in the network. Principal component regression outperformed all other prediction models, demonstrating the age and gender information encoded in the connectivity network.

The extraction of coronary artery centerlines from computed tomography angiography data is relevant in clinical practice. A large number of (semi-)automatic methods have therefore been presented for this purpose. However, prior to the work described in this paper no standardized evaluation methodology has been published to reliably evaluate and compare the performance of coronary artery centerline extraction algorithms. This paper describes the deployment of a publicly available database of coronary CTA data with corresponding reference standard derived from manually annotated centerlines, a standardized evaluation framework consisting of well-defined evaluation measures, and an on-line tool for the comparison of coronary CTA centerline extraction techniques.

This paper presents a novel method for segmenting the coronary lumen in CTA data. The method is based on graph cuts, with edge-weights depending on the intensity of the centerline, and robust kernel regression. A quantitative evaluation in 28 coronary arteries from 12 patients is performed by comparing the semi-automatic segmentations to manual annotations. This evaluation showed that the method was able to segment the coronary arteries with high accuracy, compared to manually annotated segmentations, which is reflected in a Dice coefficient of 0.85 and average symmetric surface distance of 0.22 mm.

Generation of a reference standard from multiple manually annotated datasets is a non-trivial problem. This paper discusses the weighted averaging of 3D open curves, which we used to generate a reference standard for vessel tracking data. We show how weighted averaging can be implemented by applying the Mean Shift algorithm to paths, and discuss the details of our implementation. Our approach can handle cases where the observer centerlines take different branches in a natural way. The method has been evaluated on synthetic data, and has been used to generate reference centerlines for evaluation of vessel tracking algorithms.

This paper presents a Bayesian framework for tracking of tubular structures such as vessels. Compared to conventional tracking schemes, its main advantage is its non-deterministic character, which strongly increases the robustness of the method. A key element of our approach is a dedicated observation model for tubular structures in regions with varying intensities. Furthermore, we show how the tracking method can be used to obtain a probabilistic segmentation of the tracked tubular structure. The method has been applied to track the internal carotid artery from CT angiography data of 14 patients (28 carotids) through the skull base. This is a challenging problem, owing to the close proximity of bone, overlap in intensity values of lumen voxels and (partial volume) bone voxels, and the tortuous path of the vessels. The tracking was successful in 25 cases, and the extracted path were found to be close (< 1.0mm) to manually traced paths by two observers.

Tracking of tubular elongated structures is an important goal in a wide range of biomedical imaging applications. A Bayesian tube tracking algorithm is presented that allows to easily incorporate a priori knowledge. Because probabilistic tube tracking algorithms are computationally complex, steps towards a computational efficient implementation are suggested in this paper. The algorithm is evaluated on 2D and 3D synthetic data with different noise levels and clinical CTA data. The approach shows good performance on data with high levels of Gaussian noise.