Torbjorn Vik

Validation and comparison of registration methods for free-breathing 4D lung CT

We have compared and validated image registration methods with respect to the clinically relevant use-case of lung CT max-inhale to max-exhale registration. Four fundamentally different algorithms representing main approaches for image registration were compared using clinical images. Each algorithm was assigned to a different person with extensive working knowledge of its usage. Quantitative and qualitative evaluation is performed. Whereas the methods achieve similar results in target registration error, characteristic differences come to show by closer analysis of the displacement fields.

Robust Pose Estimation and Recognition Using Non-Gaussian Modeling of Appearance Subspaces

We present an original appearance model that generalizes the usual Gaussian visual subspace model to non-Gaussian and nonparametric distributions. It can be useful for the modeling and recognition of images under difficult conditions such as large occlusions and cluttered backgrounds. Inference under the model is efficiently solved using the mean shift algorithm

Filter learning and evaluation of the computer aided visualization and analysis (CAVA) paradigm for pulmonary nodules using the LIDC-IDRI database

We present a simple rendering scheme for thoracic CT datasets which yields a color coding based on local differential geometry features rather than Hounsfield densities. The local curvatures are computed on several resolution scales and mapped onto different colors, thereby enhancing nodular and tubular structures. The rendering can be used as a navigation device to quickly access points of possible chest anomalies, in particular lung nodules and lymph nodes. The underlying principle is to use the nodule enhancing overview as a possible alternative to classical CAD approaches by avoiding explicit graphical markers. For performance evaluation we have used the LIDC-IDRI lung nodule data base. Our results indicate that the nodule-enhancing overview correlates well with the projection images produced from the IDRI expert annotations, and that we can use this measure to optimize the combination of differential geometry filters.

Validation of automatic landmark identification for atlas-based segmentation for radiation treatment planning of the head-and-neck region

Manual contouring of target volumes and organs at risk in radiation therapy is extremely time-consuming, in particular for treating the head-and-neck area, where a single patient treatment plan can take several hours to contour. As radiation treatment delivery moves towards adaptive treatment, the need for more efficient segmentation techniques will increase. We are developing a method for automatic model-based segmentation of the head and neck. This process can be broken down into three main steps: i) automatic landmark identification in the image dataset of interest, ii) automatic landmark-based initialization of deformable surface models to the patient image dataset, and iii) adaptation of the deformable models to the patient-specific anatomical boundaries of interest. In this paper, we focus on the validation of the first step of this method, quantifying the results of our automatic landmark identification method. We use an image atlas formed by applying thin-plate spline (TPS) interpolation to ten atlas datasets, using 27 manually identified landmarks in each atlas/training dataset. The principal variation modes returned by principal component analysis (PCA) of the landmark positions were used by an automatic registration algorithm, which sought the corresponding landmarks in the clinical dataset of interest using a controlled random search algorithm. Applying a run time of 60 seconds to the random search, a root mean square (rms) distance to the ground-truth landmark position of 9.5 ± 0.6 mm was calculated for the identified landmarks. Automatic segmentation of the brain, mandible and brain stem, using the detected landmarks, is demonstrated.

Local motion analysis in 4D lung CT using fast groupwise registration

Automatic tools for the propagation of contours and surfaces in time series of CT images (4DCT) are necessary e.g. for dose calculation in radio-therapy-planning. The state of the art method is the registration of 3-D frames of the sequence on a static reference volume or the registration of consecutive volume pairs. Based on the computed vector fields, the contours and surfaces are then propagated through the entire 4-D series. However, the registrations frequently suffer from common image or gating artefacts in the time series. To overcome the problem, we propose a groupwise 4-D registration approach combined with a motion model to perform the registration of a sequence of images.

Mean shift-based Bayesian image reconstruction into visual subspace

We present a new robust algorithm for reconstructing images into a linear subspace using MAP estimation. The algorithm takes into account the a priori distribution of the subspace variables and the noise is robustly modeled to allow for occlusions. The subspace distribution is estimated using nonparametric density estimation techniques. An efficient optimization scheme based on the mean shift procedure D Comaniciu et al. (2002) and on half-quadratic theory [ D Geman et al. (1992), P Charbonnier et al. (1997)] is developed, making optimization of the MAP function feasible for high-dimensional images. Preliminary results on real images demonstrate the contribution of a priori distribution modeling of sub-space variables, with respect to standard reconstruction methods over linear subspaces.

Repeatability of dose painting by numbers treatment planning in prostate cancer radiotherapy based on multiparametric magnetic resonance imaging

Dose painting by numbers (DPBN) refers to a voxel-wise prescription of radiation dose modelled from functional image characteristics, in contrast to dose painting by contours which requires delineations to define the target for dose escalation. The direct relation between functional imaging characteristics and DPBN implies that random variations in images may propagate into the dose distribution. The stability of MR-only prostate cancer treatment planning based on DPBN with respect to these variations is as yet unknown. We conducted a test-retest study to investigate the stability of DPBN for prostate cancer in a semi-automated MR-only treatment planning workflow. Twelve patients received a multiparametric MRI on two separate days prior to prostatectomy. The tumor probability (TP) within the prostate was derived from image features with a logistic regression model. Dose mapping functions were applied to acquire a DPBN prescription map that served to generate an intensity modulated radiation therapy (IMRT) treatment plan. Dose calculations were done on a pseudo-CT derived from the MRI. The TP and DPBN map and the IMRT dose distribution were compared between both MRI sessions, using the intraclass correlation coefficient (ICC) to quantify repeatability of the planning pipeline. The quality of each treatment plan was measured with a quality factor (QF). Median ICC values for the TP and DPBN map and the IMRT dose distribution were 0.82, 0.82 and 0.88, respectively, for linear dose mapping and 0.82, 0.84 and 0.94 for square root dose mapping. A median QF of 3.4% was found among all treatment plans. We demonstrated the stability of DPBN radiotherapy treatment planning in prostate cancer, with excellent overall repeatability and acceptable treatment plan quality. Using validated tumor probability modelling and simple dose mapping techniques it was shown that despite day-to-day variations in imaging data still consistent treatment plans were obtained.

Annotation-free probabilistic atlas learning for robust anatomy detection in CT images

A fully automatic method generating a whole body atlas from CT images is presented. The atlas serves as a reference space for annotations. It is based on a large collection of partially overlapping medical images and a registration scheme. The atlas itself consists of probabilistic tissue type maps and can represent anatomical variations. The registration scheme is based on an entropy-like measure of these maps and is robust with respect to field-of-view variations. In contrast to other atlas generation methods, which typically rely on a sufficiently large set of annotations on training cases, the presented method requires only the images. An iterative refinement strategy is used to automatically stitch the images to build the atlas. Affine registration of unseen CT images to the probabilistic atlas can be used to transfer reference annotations, e.g. organ models for segmentation initialization or reference bounding boxes for field-of-view selection. The robustness and generality of the method is shown using a three-fold cross-validation of the registration on a set of 316 CT images of unknown content and large anatomical variability. As an example, 17 organs are annotated in the atlas reference space and their localization in the test images is evaluated. The method yields a recall (sensitivity), specificity and precision of at least 96% and thus performs excellent in comparison to competitors.

Precise anatomy localization in CT data by an improved probabilistic tissue type atlas

Automated interpretation of CT scans is an important, clinically relevant area as the number of such scans is increasing rapidly and the interpretation is time consuming. Anatomy localization is an important prerequisite for any such interpretation task. This can be done by image-to-atlas registration, where the atlas serves as a reference space for annotations such as organ probability maps. Tissue type based atlases allow fast and robust processing of arbitrary CT scans. Here we present two methods which significantly improve organ localization based on tissue types. A first problem is the definition of tissue types, which until now is done heuristically based on experience. We present a method to determine suitable tissue types from sample images automatically. A second problem is the restriction of the transformation space: all prior approaches use global affine maps. We present a hierarchical strategy to refine this global affine map. For each organ or region of interest a localized tissue type atlas is computed and used for a subsequent local affine registration step. A three-fold cross validation on 311 CT images with different fields-of-view demonstrates a reduction of the organ localization error by 33%.

Towards a comprehensive CT image segmentation for thoracic organ radiation dose estimation and reporting

Administered dose of ionizing radiation during medical imaging is an issue of increasing concern for the patient, for the clinical community, and for respective regulatory bodies. CT radiation dose is currently estimated based on a set of very simplifying assumptions which do not take the actual body geometry and organ specific doses into account. This makes it very difficult to accurately report imaging related administered dose and to track it for different organs over the life of the patient. In this paper this deficit is addressed in a two-fold way. In a first step, the absorbed radiation dose in each image voxel is estimated based on a Monte-Carlo simulation of X-ray absorption and scattering. In a second step, the image is segmented into tissue types with different radio sensitivity. In combination this allows to calculate the effective dose as a weighted sum of the individual organ doses. The main purpose of this paper is to assess the feasibility of automatic organ specific dose estimation. With respect to a commercially applicable solution and respective robustness and efficiency requirements, we investigated the effect of dose sampling rather than integration over the organ volume. We focused on the thoracic anatomy as the exemplary body region, imaged frequently by CT. For image segmentation we applied a set of available approaches which allowed us to cover the main thoracic radio-sensitive tissue types. We applied the dose estimation approach to 10 thoracic CT datasets and evaluated segmentation accuracy and administered dose and could show that organ specific dose estimation can be achieved.