Abstract. The multiscale second order local structure of an image (Hessian) isexamined with the purpose of developing a vessel enhancement filter. A vesselness measure is obtained on the basis of all eigenvalues of the Hessian. This measure is tested on two dimensional DSA and three dimensional aortoiliac and cerebral MRA data. Its clinical utility is shown by the simultaneous noise and background suppression and vessel enhancement in maximum intensity projections and volumetric displays. 1

: This paper describes a method for the enhancement of curvilinear structures such as vessels and bronchi in 3D medical images. A 3D line enhancement filter is developed with the aim of discriminating line structures from other structures and recovering line structures of various widths. The 3D line filter is based on a combination of the eigenvalues of the 3D Hessian matrix. Multi-scale integration is formulated by taking the maximum among single-scale filter responses, and its characteristics are examined to derive criteria for the selection of parameters in the formulation. The resultant multi-scale line-filtered images provide significantly improved segmentation and visualization of curvilinear structures. The usefulness of the method is demonstrated by the segmentation and visualization of brain vessels from MRI (magnetic resonance imaging) and MRA (magnetic resonance angiography), bronchi from a chest CT, and liver vessels (portal veins) from an abdominal CT. Keywords: 3D image ...

Abstract—Interpretation of magnetic resonance angiography (MRA) is problematic due to complexities of vascular shape and to artifacts such as the partial volume effect. We present new methods to assist in the interpretation of MRA. These include methods for detection of vessel paths and for determination of branching patterns of vascular trees. They are based on the ordered region growing (ORG) algorithm that represents the image as an acyclic graph, which can be reduced to a skeleton by specifying vessel endpoints or by a pruning process. Ambiguities in the vessel branching due to vessel overlap are effectively resolved by heuristic methods that incorporate a priori knowledge of bifurcation spacing. Vessel paths are detected at interactive speeds on a 500-MHz processor using vessel endpoints. These methods apply best to smaller vessels where the image intensity peaks at the center of the lumen which, for the abdominal MRA, includes vessels whose diameter is less than 1 cm. Index Terms—Magnetic resonance angiography, skeletonization, visualization. I.

Abstract. Morphological connected set filters for extraction of filamentous details from medical images are developed. The advantages of these filters are that they are shape preserving and do not amplify noise. Two approaches are compared: (i) multi-scale filtering (ii) single-step shape filtering using connected set (or attribute) thinnings. The latter method highlights all filamentous structure in a single filtering stage, regardless of the scale. The second approach is an order of magnitude faster than the first, filtering a 256 3 volume in 41.65 s on a 400 MHz Pentium II. 1

Three-dimensional rotational angiography (3DRA) is a new and promising technique for obtaining high-resolution isotropic 3D images of vascular structures. However, due to the relatively high noise level and the presence of other background structures in clinical 3DRA images, noise reduction is inevitable. In this paper, we evaluate a number of linear and nonlinear diffusion techniques for this purpose. Specifically, we analyze the effects of these techniques on the thresholdbased visualization and quantification of vascular anomalies in 3DRA images. The results of invitro experiments indicate that edge-enhancing anisotropic diffusion filtering is most suitable: the increase in the user-dependency of visualizations and quantifications is considerably less with this technique compared to linear filtering techniques, and it is better at reducing noise near edges than isotropic nonlinear diffusion. However, in view of the memory and computation-time requirements of this technique, the latter scheme may be considered a useful alternative.