Gordon L. Kindlmann

Semi-automatic generation of transfer functions for direct volume rendering

Although direct volume rendering is a powerful tool for visualizing complex structures within volume data, the size and complexity of the parameter space controlling the rendering process makes generating an informative rendering challenging. In particular, the specification of the transfer function-the mapping from data values to renderable optical properties-is frequently a time consuming and unintuitive task. Ideally, the data being visualized should itself suggest an appropriate transfer function that brings out the features of interest without obscuring them with elements of little importance. We demonstrate that this is possible for a large class of scalar volume data, namely that where the regions of interest are the boundaries between different materials. A transfer function which makes boundaries readily visible can be generated from the relationship between three quantities: the data value and its first and second directional derivatives along the gradient direction. A data structure we term the histogram volume captures the relationship between these quantities throughout the volume in a position independent, computationally efficient fashion. We describe the theoretical importance of the quantities measured by the histogram volume, the implementation issues in its calculation, and a method for semiautomatic transfer function generation through its analysis. We conclude with results of the method on both idealized synthetic data as well as real world datasets.

Multidimensional transfer functions for interactive volume rendering

Most direct volume renderings produced today employ 1D transfer functions which assign color and opacity to the volume based solely on the single scalar quantity which comprises the data set. Though they have not received widespread attention, multi-dimensional transfer functions are a very effective way to extract materials and their boundaries for both scalar and multivariate data. However, identifying good transfer functions is difficult enough in 1D, let alone 2D or 3D. This paper demonstrates an important class of 3D transfer functions for scalar data, and describes the application of multi-dimensional transfer functions to multivariate data. We present a set of direct manipulation widgets that make specifying such transfer functions intuitive and convenient. We also describe how to use modern graphics hardware to both interactively render with multidimensional transfer functions and to provide interactive shadows for volumes. The transfer functions, widgets and hardware combine to form a powerful system for interactive volume exploration.

DTI measures in crossing-fibre areas: Increased diffusion anisotropy reveals early white matter alteration in MCI and mild Alzheimer's disease

Though mild cognitive impairment is an intermediate clinical state between healthy aging and Alzheimers disease (AD), there are very few whole-brain voxel-wise diffusion MRI studies directly comparing changes in healthy control, mild cognitive impairment (MCI) and AD subjects. Here we report whole-brain findings from a comprehensive study of diffusion tensor indices and probabilistic tractography obtained in a very large population of healthy controls, MCI and probable AD subjects. As expected from the literature, all diffusion indices converged to show that the cingulum bundle, the uncinate fasciculus, the entire corpus callosum and the superior longitudinal fasciculus are the most affected white matter tracts in AD. Significant differences between MCI and AD were essentially confined to the corpus callosum. More importantly, we introduce for the first time in a degenerative disorder an application of a recently developed tensor index, the mode of anisotropy, as well as probabilistic crossing-fibre tractography. The mode of anisotropy specifies the type of anisotropy as a continuous measure reflecting differences in shape of the diffusion tensor ranging from planar (e.g., in regions of crossing fibres from two fibre populations of similar density or regions of kissing fibres) to linear (e.g., in regions where one fibre population orientation predominates), while probabilistic crossing-fibre tractography allows to accurately trace pathways from a crossing-fibre region. Remarkably, when looking for whole-brain diffusion differences between MCI patients and healthy subjects, the only region with significant abnormalities was a region of crossing fibres in the centrum semiovale, showing an increased mode of anisotropy. The only white matter region demonstrating a significant difference in correlations between neuropsychological scores and a diffusion measure (mode of anisotropy) across the three groups was the same region of crossing fibres. Further examination using probabilistic tractography established explicitly and quantitatively that this previously unreported increase of mode and co-localised increase of fractional anisotropy was explained by a relative preservation of motor-related projection fibres (at this early stage of the disease) crossing the association fibres of the superior longitudinal fasciculus. These findings emphasise the benefit of looking at the more complex regions in which spared and affected pathways are crossing to detect very early alterations of the white matter that could not be detected in regions consisting of one fibre population only. Finally, the methods used in this study may have general applicability for other degenerative disorders and, beyond the clinical sphere, they could contribute to a better quantification and understanding of subtle effects generated by normal processes such as visuospatial attention or motor learning.

The transfer function bake-off

A liquid fuel pumping apparatus for supplying fuel to an internal combustion engine includes a member which is movable against the action of a spring system by means of centrifugally operable weights. The member is coupled to a fuel control rod of the pumping apparatus, and the spring system includes first and second springs. The first spring is preloaded so as to be deflected by the weights only when the speed of the engine attains a predetermined value. The other spring is deflected by the weights at lower engine speeds and also acts to transmit the force exerted by the weights to the first spring.

Superquadric tensor glyphs

Tensor field visualization is a challenging task due in part to the multi-variate nature of individual tensor samples. Glyphs convey tensor variables by mapping the tensor eigenvectors and eigenvalues to the orientation and shape of a geometric primitive, such as a cuboid or ellipsoid. Though widespread, cuboids and ellipsoids have problems of asymmetry and visual ambiguity. Cuboids can display misleading orientation for tensors with underlying rotational symmetry. Ellipsoids differing in shape can be confused, from certain viewpoints, because of similarities in profile and shading. This paper addresses the problems of asymmetry and ambiguity with a new tunable continuum of glyphs based on superquadric surfaces. Superquadric tensor glyphs enjoy the necessary symmetry properties of ellipsoids, while also imitating cuboids and cylinders to better convey shape and orientation, where appropriate. The new glyphs are demonstrated on fields of diffusion tensors from the human brain.

Tensorlines: advection-diffusion based propagation through diffusion tensor fields

Tracking linear features through tensor field datasets is an open research problem with widespread utility in medical and engineering disciplines. Existing tracking methods, which consider only the preferred local diffusion direction as they propagate, fail to accurately follow features as they enter regions of local complexity. This shortcoming is a result of partial voluming; that is, voxels in these regions often contain contributions from multiple features. These combined contributions result in ambiguities when deciding local primary feature orientation based solely on the preferred diffusion direction. We introduce a novel feature extraction method which we term tensorline propagation. Our method resolves the above ambiguity by incorporating information about the nearby orientation of the feature, as well as the anisotropic classification of the local tensor. The nearby orientation information is added in the spirit of an advection term in a standard diffusion based propagation technique, and has the effect of stabilizing the tracking. To demonstrate the efficacy of tensorlines, we apply this method to the neuroscience problem of tracking white-matter bundles within the brain.

Orthogonal tensor invariants and the analysis of diffusion tensor magnetic resonance images

This paper outlines the mathematical development and application of two analytically orthogonal tensor invariants sets. Diffusion tensors can be mathematically decomposed into shape and orientation information, determined by the eigenvalues and eigenvectors, respectively. The developments herein orthogonally decompose the tensor shape using a set of three orthogonal invariants that characterize the magnitude of isotropy, the magnitude of anisotropy, and the mode of anisotropy. The mode of anisotropy is useful for resolving whether a region of anisotropy is linear anisotropic, orthotropic, or planar anisotropic. Both tensor trace and fractional anisotropy are members of an orthogonal invariant set, but they do not belong to the same set. It is proven that tensor trace and fractional anisotropy are not mutually orthogonal measures of the diffusive process. The results are applied to the analysis and visualization of diffusion tensor magnetic resonance images of the brain in a healthy volunteer. The theoretical developments provide a method for generating scalar maps of the diffusion tensor data, including novel fractional anisotropy maps that are color encoded for the mode of anisotropy and directionally encoded colormaps of only linearly anisotropic structures, rather than of high fractional anisotropy structures. Magn Reson Med, 2006.

A geometric analysis of diffusion tensor measurements of the human brain

The degree of diffusion tensor anisotropy is often associated with the organization of structural tissues such as white matter. Numerous measures of diffusion anisotropy have been proposed, which could lead to confusion in interpreting and comparing results from different studies. In this study, a new method for representing the diffusion tensor shape, called the three‐phase (3P) plot, is described. This is a graphical technique based upon a barycentric coordinate system, which weights the tensor shape by a combination of linear, cylindrical, and spherical shape factors. This coordinate system can be used to map and potentially segment different tissues based upon the tensor shape. In addition, the 3P plot can be used to examine the shape properties of existing measures of diffusion anisotropy. In this paper, the 3P plot is used to compare four well‐known anisotropy measures: the anisotropy index, the fractional anisotropy, the relative anisotropy, and the volume fraction. Computer simulations and diffusion tensor images of normal brains were obtained to study the properties of this new mapping technique. Magn Reson Med 44:283–291, 2000.

Strategies for direct volume rendering of diffusion tensor fields

Diffusion-weighted magnetic resonance imaging is a relatively new modality capable of elucidating the fibrous structure of certain types of tissue, such as the white matter within the brain. One tool for interpreting this data is volume rendering because it permits the visualization of three dimensional structure without a prior segmentation process. In order to use volume rendering, however, we must develop methods for assigning opacity and color to the data, and create a method to shade the data to improve the legibility of the rendering. Previous work introduced three such methods: barycentric opacity maps, hue-balls (for color), and lit-tensors (for shading). The paper expands on and generalizes these methods, describing and demonstrating further means of generating opacity, color, and shading from the tensor information. We also propose anisotropic reaction-diffusion volume textures as an additional tool for visualizing the structure of diffusion data. The patterns generated by this process can be visualized on their own or they can be used to supplement the volume rendering strategies described in the rest of the paper. Finally, because interpolation between data points is a fundamental issue in volume rendering, we conclude with a discussion and evaluation of three distinct interpolation methods suitable for diffusion tensor MRI data.

Diffusion Tensor Visualization with Glyph Packing

A common goal of multivariate visualization is to enable data inspection at discrete points, while also illustrating larger-scale continuous structures. In diffusion tensor visualization, glyphs are typically used to meet the first goal, and methods such as texture synthesis or fiber tractography can address the second. We adapt particle systems originally developed for surface modeling and anisotropic mesh generation to enhance the utility of glyph-based tensor visualizations. By carefully distributing glyphs throughout the field (either on a slice, or in the volume) into a dense packing, using potential energy profiles shaped by the local tensor value, we remove undue visual emphasis of the regular sampling grid of the data, and the underlying continuous features become more apparent. The method is demonstrated on a DT-MRI scan of a patient with a brain tumor