Related work can be classified in the following two main topics.
2.1. Diffusion Tensor Imaging and Tractography
Conventional structural imaging techniques such as T1-, T2- and proton density-weighted imaging generally create high contrast between major tissue groups in the brain, which are: Gray Matter (GM), White Matter (WM), and Cerebrospinal Fluid (CSF). Such structural imaging techniques are generally well suited for the study of tissue macrostructure yet provide little insight into the orientation of white matter fibers. DWI [
17,
18] is a variant of conventional Magnetic Resonance Imaging (MRI) based on the tissue water diffusion rate, which is better suited for the study of white matter pathways. Although DWI refers to the contrast of the acquired images, DTI is a specific type of modeling (or abstraction) of the DWI datasets, in which diffusivity of water molecules is represented by tensors [
10]. As an in-vivo non-destructive technique that requires no chemical tracers, DTI is presently one of the most promising methods for the study of white matter architecture in living humans. DTI provides quantitative estimates of white matter integrity and orientation by measuring molecular diffusion of water molecules (or Brownian motion). It is based on the phenomenon of diffusion anisotropy in the nerve tissue: water molecules diffuse faster along the neural fiber direction and slower in the fiber-transverse direction.
DTI tensor datasets record, essentially, the spectrum of diffusion strengths per diffusion direction at every point in a 3D scan. Such datasets can be visualized using glyph techniques, which render the local major, medium, and minor eigenvalues of the diffusion tensor, scaled by their corresponding eigenvalues, over a dense sampling of the 3D scanned volume [
29].
Tractography methods improve in the above by reducing occlusion and also explicitly showing the fiber tract paths by essentially integrating streamlines in the major direction of the DTI eigenvector, seeded from a suitably chosen point set [
11,
12]. More elaborate approaches include streamlines defined by tensor deflection [
10], bi-tensor modeling [
19] and streamsurface tracking [
20]. Fiber tracking has several applications, including noninvasive visualization of white matter pathways, segmenting of specific tracts in the brain for image analysis, and relating white matter tract anatomy to brain tumors and lesions in patients who are candidates for neurosurgery [
10].
Tractography essentially delivers large sets of (hundreds of thousands of) 3D fiber trails embedded in the 3D scan volume. Visualizing these next is challenging. Conventional line rendering is a simple and efficient baseline rendering technique to show such datasets that has been employed since fiber tracking was first introduced [
12,
21]. Yet, such naive rendering of the densely seeded and criss-crossing 3D trails creates too much clutter and occlusion for many tasks to be accomplished effectively. Numerous visualization tools have been proposed to leverage interactivity to explore large DTI datasets, e.g., [
22]. However, navigating the large space of parameter settings can be daunting for users. More extensive illumination [
23], ambient occlusion [
24] and alpha-blending [
25] have been applied to line rendering of fiber tracts to further improve their visual exploration. Rendering 3D tubes with Phong shading [
26] is based on more detailed geometric modeling than line rendering, and can produce even higher quality visualizations. The methods by Stoll et al. [
27] and Merhof et al. [
28] produce similar shading using screenspace techniques and incorporate illustrative rendering techniques [
15] as well. Hyperstreamlines [
25,
29] are an extension of cylindrical tubes that provide a richer representation of the DTI field.
Clustering forms separate class of methods that help visualizing large tract sets. Such methods detect trails that are similar in terms of spatial position, diffusion values, and possibly other attributes (e.g., curvature, length) and group them into clusters. These in turn enable one to create a simplified visualization by e.g., rendering each cluster with a different color or reducing the cluster to a simpler representative, which is next visualized. Xu et al. [
30] use the DBSCAN clustering algorithm [
31] to create such clusters and combine them with user specification of regions of interest to create a rich palette of focus-and-context fiber tract visualization. Poco et al. [
32] approach this by reducing every 3D fiber to a high-dimensional feature vector that represents position, geometry, and smoothness. Next, dimensionality reduction [
33] is used to create 2D scatterplots where similar fibers appear as point clusters. This enables one to easily select similar-fiber bundles by simply selecting point clusters in the projection. Comparisons of fiber clustering methods are presented in [
34,
35]. Interestingly, some clustering methods, such as DBSCAN or mean shift [
36], are related to density estimation, which is also used by CUBu [
8], the bundling technique that we adapt for our DTI visualization (see further
Section 2.2). The key difference between our method and clustering is that we
deform the implicitly clustered fibers to simplify the resulting visualization.
Surface reconstruction or extraction from DTI datasets have also been applied as a means for indirect volume visualization [
29]. Merhof et al. [
28] have shown that encompassing fiber tract bundles with isosurfaces yields a preferable representation for use in neurosurgery. Ridge and valley surfaces [
29] are demonstrated to capture the cores of sheet-like tracts. Visualizing such surfaces provides important added value compared to all fiber-based alternatives. Indeed, even if there is a dense surface-like distribution of fiber tracts, creating the appearance of a
surface from this (without gaps) can be challenging when applying line-rendering-based techniques. Besides explicit surface extraction, such surfaces can be also created
implicitly in the image space. The Depth-Dependent Halos (DDH) method [
13,
14,
37] performs this merging of structures quite well, as colinear fibers are visually combined into thicker bundles with illustrative halos. The major downside of DDH is that the black-and-white rendering (a stylistic choice motivated by nonphotorealistic rendering) offers limited scalar visualization capabilities. Our method, discussed next, can reproduce DDH but also add shading and information color-coding to the fiber rendering.
2.2. Trail and Edge-Bundling Techniques
In contrast to the above-mentioned rendering techniques, which modify how a 3D trail set is
depicted, trail and edge-bundling techniques modify the actual
trail set to emphasize structures of interest. Bundling essentially aims to (a) identify trails that are spatially close and similar from the perspective of one or several data attributes; and (b) deform these trails so they get spatially closer, so as to create more visual structure in the ensuing rendering thereof. Applications of edge bundling include graph drawing simplification [
8,
38], trajectory exploration [
7], eye-tracking analysis [
5] and streamline bundling [
39].
Bundling of 2D datasets has been extensively explored. Image-based edge-bundling techniques, including SBEB [
40], ADEB [
5], KDEEB [
38], CUBu [
8], and FFTEB [
41] show that the problem of bundling of 2D trails, as well as the rendering of resulting bundles using a variety of styles, can be efficiently and effectively addressed using image processing approaches that scale well on graphics hardware. Bundling is however far less explored for 3D volumetric trail datasets such as DTI tracts. Bottger et al. [
42,
43] have proposed mean-shift bundling of 3D connectivity graphs obtained through functional Magnetic Resonance Imaging (fMRI). Simplification of anatomical connectivity using edge bundling is a relatively new topic, one example thereof being KDEEB [
38] applied to fiber tracts [
25,
44]. Other examples of fiber tract simplification methods include multi-scale local fiber tract contraction [
16] and two-dimensional neural maps [
45]. To our knowledge, such 3D methods have two main limitations: They cannot, in the same time (1) recover well plausible anatomical structures, such as the mix of surfaces and tube-like structures present in a DTI dataset; and (2) allow interactive visual exploration of DTI datasets of hundreds of thousands of fibers in real time using different rendering styles. For example, the method in [
16], which is technically closest to ours, requires hours to complete on a dataset of 78 K fiber tracts. At a higher level, bundling-and-rendering methods for DTI data have evolved largely separated, with the notable exception of [
16]. This, we believe, is an artificial separation of two method classes which ultimately aim to the same goal—presenting a large, complex, dense, volumetric 3D trail dataset in a
simplified manner to the user.
Summarizing the above, our contribution to the creation of simplified visualization of large DTI trail sets is three-fold:
We present a method for bundling 3D DTI tracts that allows the user to preserve relevant underlying anatomical elements (fiber tubular bundles, sheets, and manifolds) in a controlled manner;
We propose several rendering styles of the 3D bundled structures that emphasize several aspects of interest in the data;
We propose a joint implementation of the above two points that can handle 3D trail datasets of hundreds of thousands of trails in real time on consumer graphics hardware.
As our method builds atop of the backbone of the CUBu technique [
8]. We inherit from CUBu the kernel density-based advection of trails that ultimately creates the bundles, and the full GPU-based bundling computation for computational scalability. As described next, we modify CUBu in several respects: We implement the entire pipeline in 3D (using a volumetric density map,
Section 3.2); allow trail endpoints to bundle under specific constraints (
Section 3.2.1); use the DTI volume anisotropy rather than the CUBu isotropic bundling (
Section 3.2.2); reseed tracts in sparsely populated volume areas rather than bundling a fixed, predefined, trail set (
Section 3.2.3); and propose several rendering modes designed to emphasize structure specific to DTI fiber sheets (
Section 3.3).