Gunnar Läthén

Blood vessel segmentation using multi-scale quadrature filtering

The segmentation of blood vessels is a common problem in medical imaging and various applications are found in diagnostics, surgical planning, training and more. Among many different techniques, the use of multiple scales and line detectors is a popular approach. However, the typical line filters used are sensitive to intensity variations and do not target the detection of vessel walls explicitly. In this article, we combine both line and edge detection using quadrature filters across multiple scales. The filter result gives well defined vessels as linear structures, while distinct edges facilitate a robust segmentation. We apply the filter output to energy optimization techniques for segmentation and show promising results in 2D and 3D to illustrate the behavior of our method. The conference version of this article received the best paper award in the bioinformatics and biomedical applications track at ICPR 2008.

Modified Gradient Search for Level Set Based Image Segmentation

Level set methods are a popular way to solve the image segmentation problem. The solution contour is found by solving an optimization problem where a cost functional is minimized. Gradient descent methods are often used to solve this optimization problem since they are very easy to implement and applicable to general nonconvex functionals. They are, however, sensitive to local minima and often display slow convergence. Traditionally, cost functionals have been modified to avoid these problems. In this paper, we instead propose using two modified gradient descent methods, one using a momentum term and one based on resilient propagation. These methods are commonly used in the machine learning community. In a series of 2-D/3-D-experiments using real and synthetic data with ground truth, the modifications are shown to reduce the sensitivity for local optima and to increase the convergence rate. The parameter sensitivity is also investigated. The proposed methods are very simple modifications of the basic method, and are directly compatible with any type of level set implementation. Downloadable reference code with examples is available online.

Phase based level set segmentation of blood vessels

The segmentation and analysis of blood vessels has received much attention in the research community. The results aid numerous applications for diagnosis and treatment of vascular diseases. Here we use level set propagation with local phase information to capture the boundaries of vessels. The basic notion is that local phase, extracted using quadrature filters, allows us to distinguish between lines and edges in an image. Noting that vessels appear either as lines or edge pairs, we integrate multiple scales and capture information about vessels of varying width. The outcome is a ¿global¿ phase which can be used to drive a contour robustly towards the vessel edges. We show promising results in 2D and 3D. Comparison with a related method gives similar or even better results and at a computational cost several orders of magnitude less. Even with very sparse initializations, our method captures a large portion of the vessel tree.

Momentum Based Optimization Methods for Level Set Segmentation

Segmentation of images is often posed as a variational problem. As such, it is solved by formulating an energy functional depending on a contour and other image derived terms. The solution of the segmentation problem is the contour which extremizes this functional. The standard way of solving this optimization problem is by gradient descent search in the solution space, which typically suffers from many unwanted local optima and poor convergence. Classically, these problems have been circumvented by modifying the energy functional. In contrast, the focus of this paper is on alternative methods for optimization. Inspired by ideas from the machine learning community, we propose segmentation based on gradient descent with momentum. Our results show that typical models hampered by local optima solutions can be further improved by this approach. We illustrate the performance improvements using the level set framework.

Automatic Tuning of Spatially Varying Transfer Functions for Blood Vessel Visualization

Computed Tomography Angiography (CTA) is commonly used in clinical routine for diagnosing vascular diseases. The procedure involves the injection of a contrast agent into the blood stream to increase the contrast between the blood vessels and the surrounding tissue in the image data. CTA is often visualized with Direct Volume Rendering (DVR) where the enhanced image contrast is important for the construction of Transfer Functions (TFs). For increased efficiency, clinical routine heavily relies on preset TFs to simplify the creation of such visualizations for a physician. In practice, however, TF presets often do not yield optimal images due to variations in mixture concentration of contrast agent in the blood stream. In this paper we propose an automatic, optimization-based method that shifts TF presets to account for general deviations and local variations of the intensity of contrast enhanced blood vessels. Some of the advantages of this method are the following. It computationally automates large parts of a process that is currently performed manually. It performs the TF shift locally and can thus optimize larger portions of the image than is possible with manual interaction. The method is based on a well known vesselness descriptor in the definition of the optimization criterion. The performance of the method is illustrated by clinically relevant CT angiography datasets displaying both improved structural overviews of vessel trees and improved adaption to local variations of contrast concentration.

A Fast Optimization Method for Level Set Segmentation

Level set methods are a popular way to solve the image segmentation problem in computer image analysis. A contour is implicitly represented by the zero level of a signed distance function, and evolved according to a motion equation in order to minimize a cost function. This function defines the objective of the segmentation problem and also includes regularization constraints. Gradient descent search is the de facto method used to solve this optimization problem. Basic gradient descent methods, however, are sensitive for local optima and often display slow convergence. Traditionally, the cost functions have been modified to avoid these problems. In this work, we instead propose using a modified gradient descent search based on resilient propagation (Rprop), a method commonly used in the machine learning community. Our results show faster convergence and less sensitivity to local optima, compared to traditional gradient descent.

Evaluation of transfer function methods in direct volume rendering of the blood vessel lumen

Visualization of contrast enhanced blood vessels in CT angiography data presents a challenge due to varying concentration of the contrast agent. The purpose of this work is to evaluate the correctness (effectiveness) in visualizing the vessel lumen using two different 3D visualization strategies, thereby assessing the feasibility of using such visualizations for diagnostic decisions. We compare a standard visualization approach with a recent method which locally adapts to the contrast agent concentration. Both methods are evaluated in a parallel setting where the participant is instructed to produce a complete visualization of the vessel lumen, including both large and small vessels, in cases of calcified vessels in the legs. The resulting visualizations are thereafter compared in a slice viewer to assess the correctness of the visualized lumen. The results indicate that the participants generally overestimated the size of the vessel lumen using the standard visualization, whereas the locally adaptive method better conveyed the true anatomy. The participants did find the interpretation of the locally adaptive method to be less intuitive, but also noted that this did not introduce any prohibitive complexity in the work flow. The observed trends indicate that the visualized lumen strongly depends on the width and placement of the applied transfer function and that this dependency is inherently local rather than global. We conclude that methods that permit local adjustments, such as the method investigated in this study, can be beneficial to certain types of visualizations of large vascular trees.

Non-ring Filters for Robust Detection of Linear Structures

Many applications in image analysis include the problem of linear structure detection, e.g. segmentation of blood vessels in medical images, roads in satellite images, etc. A simple and efficient solution is to apply linear filters tuned to the structures of interest and extract line and edge positions from the filter output. However, if the filter is not carefully designed, artifacts such as ringing can distort the results and hinder a robust detection. In this paper, we study the ringing effects using a common Gabor filter for linear structure detection, and suggest a method for generating non-ring filters in 2D and 3D. The benefits of the non-ring design are motivated by results on both synthetic and natural images.