Theo van Walsum

Standardized evaluation methodology and reference database for evaluating coronary artery centerline extraction algorithms.

Efficiently obtaining a reliable coronary artery centerline from computed tomography angiography data is relevant in clinical practice. Whereas numerous methods have been presented for this purpose, up to now no standardized evaluation methodology has been published to reliably evaluate and compare the performance of the existing or newly developed coronary artery centerline extraction algorithms. This paper describes a standardized evaluation methodology and reference database for the quantitative evaluation of coronary artery centerline extraction algorithms. The contribution of this work is fourfold: (1) a method is described to create a consensus centerline with multiple observers, (2) well-defined measures are presented for the evaluation of coronary artery centerline extraction algorithms, (3) a database containing 32 cardiac CTA datasets with corresponding reference standard is described and made available, and (4) 13 coronary artery centerline extraction algorithms, implemented by different research groups, are quantitatively evaluated and compared. The presented evaluation framework is made available to the medical imaging community for benchmarking existing or newly developed coronary centerline extraction algorithms.

Fluid Flow Visualization

This chapter presents an overview of techniques for visualization of fluid flow data. As a starting point, a brief introduction to experimental flow visualization is given. The rest of the chapter concentrates on computer graphics flow visualization. A pipeline model of the flow visualization process is used as a basis for presentation. Conceptually, this process centres around visualization mapping, or the translation of physical flow parameters to visual representations. Starting from a set of standard mappings partly based on equivalents from experimental visualization, a number of data preparation techniques is described, to prepare the flow data for visualization. Next, a number of perceptual effects and rendering techniques are described, and some problems in visual presentation are discussed. The chapter ends with some concluding remarks and suggestions for future development.

Coronary Lumen Segmentation Using Graph Cuts and Robust Kernel Regression

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.

Reproducibility, accuracy, and predictors of accuracy for the detection of coronary atherosclerotic plaque composition by computed tomography: an ex vivo comparison to intravascular ultrasound.

Purpose:To determine the reproducibility, accuracy, and predictors of accuracy of computed tomography (CT) angiography to detect and characterize coronary atherosclerotic plaque as compared with intravascular ultrasound. Methods:Ten ex vivo human coronary arteries were imaged in a moving phantom by dual-source CT (collimation: 0.6 mm, reconstructed slice thickness: 0.4 mm) and intravascular ultrasound (IVUS). Coregistered cross-sections were assessed at 0.4 mm intervals for the presence and composition of atherosclerotic plaque (noncalcified, mixed, and calcified) on CT and IVUS by independent readers to determine reader agreement and diagnostic accuracy. Quantitative measurements of lumen and plaque area, plaque eccentricity, and intimal thickness on IVUS were used to determine predictors for the detection of noncalcified plaque by CT. Results:Within 1002 coregistered cross-sections, the interobserver agreement to detect plaque on CT was K = 0.48, K = 0.42, and K = 1.00 for noncalcified, mixed, and calcified plaque; respectively. The sensitivity and specificity of CT was 57% out of 84% for noncalcified, 32% of 92% for mixed, and 56% of 93% for calcified plaque when compared with IVUS; respectively. Misclassification occurred in 68% of mixed and 43% of noncalcified plaques. The odds of detecting noncalcified plaque in CT independently increased by 56% (95% CI: 47%–77%, P < 0.0001) with every 0.1 mm increase in maximum intimal thickness as measured by IVUS. Detection rate for noncalcified plaques was poor for plaques <1 mm (36%) but excellent for plaques >1 mm maximal intimal thickness (90%). Conclusion:Reader agreement and diagnostic accuracy for the detection of coronary atherosclerotic plaque vary with plaque composition. Intimal thickness independently predicts detection of noncalcified plaque by CT with excellent sensitivity for >1 mm thick plaques.

Bayesian tracking of tubular structures and its application to carotid arteries in CTA

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.

Bayesian tracking of elongated structures in 3D images

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.

Bone Displacement and the Role of Longitudinal Ligaments During Balloon Vertebroplasty in Traumatic Thoracolumbar Fractures

Study Design. In a human cadaveric burst fracture model with and without longitudinal ligament damage, the amount of anterior and posterior bone displacement (ABD, PBD) during balloon vertebroplasty after pedicle-screw instrumentation was investigated quantitatively. Objectives. To investigate, in a burst fracture model with and without longitudinal ligament damage, the amount of ABD, PBD, and cement leakage at various phases during balloon vertebroplasty in combination with pedicle-screw instrumentation. Summary of Background Data. The role of intact longitudinal ligaments in traumatic spine fractures, for prevention of bone retropulsion and subsequent reduction, has been discussed in several studies but is still up for debate. In a recent human cadaveric burst fracture study, inflatable bone tamps and calcium phosphate cement were used for the augmentation of the anterior column after pedicle-screw instrumentation. The additional balloon vertebroplasty procedure was found to be feasible and safe, but no data pertaining to unwarranted bone displacement or cement leakage during the procedure are available for burst fractures with damaged longitudinal ligaments. Methods. Ten thoracic and 10 lumbar burst fractures, with rotation or flexion components, were created, and balloon vertebroplasty with calcium phosphate cement was performed after pedicle-screw instrumentation. Volumetric datasets (using the 3-dimensional (3D) rotational x-ray imaging technique) of the fractures were obtained during the following phases: intact, fractured, after reduction and stabilization with pedicle-screws, after inflation of the balloons, after deflation and removal of the balloons, and after injection of the cement. The amount of ABD and PBD was measured on reconstructed sagittal images and recorded together with the presence of extracorporal cement leakage. The continuity of the longitudinal ligaments was assessed after anatomic dissection. Results. During the balloon vertebroplasty procedure, a significant (P < 0.05) increase of ABD (at both thoracic and lumbar level) and PBD (thoracic level) occurred after inflation of the balloons. After deflation and subsequent injection of the cement, however, the ABD and PBD returned to the preinflation levels. The absolute amount of ABD and PBD (<1 mm) during inflation was considered of little clinical importance. No differences in ABD or PBD were observed for specimens with or without continuity of the corresponding longitudinal ligament, irrespective of the level, at any of the phases during the experiment (P > 0.5 in all cases). A small amount of cement leakage was observed in the psoas compartment of one specimen with intact longitudinal ligaments. Conclusions. It is suggested that balloon vertebroplasty after pedicle-screw instrumentation may safely be used, in terms of bone displacement and cement leakage, in fracture types where damage to longitudinal ligaments is to be expected.

Semiautomatic carotid lumen segmentation for quantification of lumen geometry in multispectral MRI

Quantitative information about the geometry of the carotid artery bifurcation is relevant for investigating the onset and progression of atherosclerotic disease. This paper proposes an automatic approach for quantifying the carotid bifurcation angle, carotid area ratio, carotid bulb size and the vessel tortuosity from multispectral MRI. First, the internal and external carotid centerlines are determined by finding a minimum cost path between user-defined seed points where the local costs are based on medialness and intensity. The minimum cost path algorithm is iteratively applied after curved multi-planar reformatting to refine the centerline. Second, the carotid lumen is segmented using a topology preserving geodesic active contour which is initialized by the extracted centerlines and steered by the MR intensities. Third, the bifurcation angle and vessel tortuosity are automatically extracted from the segmented lumen. The methods for centerline tracking and lumen segmentation are evaluated by comparing their accuracy to the inter- and intra-observer variability on 48 datasets (96 carotid arteries) acquired as part of a longitudinal population study. The evaluation reveals that 94 of 96 carotid arteries are segmented successfully. The distance between the tracked centerlines and the reference standard (0.33 mm) is similar to the inter-observer variation (0.32 mm). The lumen segmentation accuracy (average DSC=0.89, average mean absolute surface distance=0.31 mm) is close to the inter-observer variation (average dice=0.92, average mean surface distance=0.23 mm). The correlation coefficient of manually and automaticly derived bifurcation angle, carotid proximal area ratio, carotid proximal bulb size and vessel totuosity quantifications are close to the correlation of these measures between observers. This demonstrates that the automated method can be used for replacing manual centerline annotation and manual contour drawing for lumen segmentation in MRIs data prior to quantifying the carotid bifurcation geometry.

Accuracy evaluation of direct navigation with an isocentric 3D rotational X-ray system

Minimally invasive interventions are often performed under fluoroscopic guidance. Drawbacks of fluoroscopic guidance are the fact that the presented images are 2D projections and that both the patient and the clinician are exposed to radiation. Image-guided navigation using pre-interventionally acquired 3D MR or CT data is an alternative. However, this often requires invasive anatomical landmark-based, marker-based or surface-based image-to-patient registration. In this paper, a coupling between an image-guided navigation system and an intraoperative C-arm X-ray device with 3D imaging capabilities (3D rotational X-ray (3DRX) system) that enables direct navigation without invasive image-to-patient registration on 3DRX volumes, is described and evaluated. The coupling is established in a one-time preoperative calibration procedure. The individual steps in the registration procedure are explained and evaluated. The acquired navigation accuracy using this coupling is approximately one millimeter.

Feasibility of multimodal deformable registration for head and neck tumor treatment planning

PURPOSE To investigate the feasibility of using deformable registration in clinical practice to fuse MR and CT images of the head and neck for treatment planning. METHOD AND MATERIALS A state-of-the-art deformable registration algorithm was optimized, evaluated, and compared with rigid registration. The evaluation was based on manually annotated anatomic landmarks and regions of interest in both modalities. We also developed a multiparametric registration approach, which simultaneously aligns T1- and T2-weighted MR sequences to CT. This was evaluated and compared with single-parametric approaches. RESULTS Our results show that deformable registration yielded a better accuracy than rigid registration, without introducing unrealistic deformations. For deformable registration, an average landmark alignment of approximatively 1.7 mm was obtained. For all the regions of interest excluding the cerebellum and the parotids, deformable registration provided a median modified Hausdorff distance of approximatively 1 mm. Similar accuracies were obtained for the single-parameter and multiparameter approaches. CONCLUSIONS This study demonstrates that deformable registration of head-and-neck CT and MR images is feasible, with overall a significanlty higher accuracy than for rigid registration.