Vishal Patel

Mesh-based spherical deconvolution: A flexible approach to reconstruction of non-negative fiber orientation distributions

Diffusion-weighted MRI has enabled the imaging of white matter architecture in vivo. Fiber orientations have classically been assumed to lie along the major eigenvector of the diffusion tensor, but this approach has well-characterized shortcomings in voxels containing multiple fiber populations. Recently proposed methods for recovery of fiber orientation via spherical deconvolution utilize a spherical harmonics framework and are susceptible to noise, yielding physically-invalid results even when additional measures are taken to minimize such artifacts. In this work, we reformulate the spherical deconvolution problem onto a discrete spherical mesh. We demonstrate how this formulation enables the estimation of fiber orientation distributions which strictly satisfy the physical constraints of realness, symmetry, and non-negativity. Moreover, we analyze the influence of the flexible regularization parameters included in our formulation for tuning the smoothness of the resultant fiber orientation distribution (FOD). We show that the method is robust and reliable by reconstructing known crossing fiber anatomy in multiple subjects. Finally, we provide a software tool for computing the FOD using our new formulation in hopes of simplifying and encouraging the adoption of spherical deconvolution techniques.

K-SVD for HARDI denoising

Noise is an important concern in high-angular resolution diffusion imaging studies because it can lead to errors in downstream analyses of white matter structure. To address this issue, we investigate a new approach for denoising diffusion-weighted data sets based on the K-SVD algorithm. We analyze its characteristics using both simulated and biological data and compare its performance with existing methods. Our results show that K-SVD provides robust and effective noise reduction and is practical for use in high-volume applications.

LONI MiND: Metadata in NIfTI for DWI

A wide range of computational methods have been developed for reconstructing white matter geometry from a set of diffusion-weighted images (DWIs), and many clinical studies rely on publicly-available implementations of these methods for analyzing DWI datasets. Unfortunately, the poor interoperability between DWI analysis tools often effectively restricts users to the algorithms provided by a single software suite, which may be suboptimal relative to those in other packages, or outdated given recent developments in the field. A major barrier to data portability and the interoperability between DWI analysis tools is the lack of a standard format for representing and communicating essential DWI-related metadata at various stages of post-processing. In this report, we address this issue by developing a framework for storing metadata in NIfTI for DWI (MiND). We utilize the standard NIfTI format extension mechanism to store essential DWI metadata in an extended header for multiple commonly-encountered DWI data structures. We demonstrate the utility of this approach by implementing a full suite of tools for DWI analysis workflows which communicate solely through the MiND mechanism. We also show that the MiND framework allows for simple, direct DWI data visualization, and we illustrate its effectiveness by constructing a group atlas for 330 subjects using solely MiND-centric tools for DWI processing. Our results indicate that the MiND framework provides a practical solution to the problem of interoperability between DWI analysis tools, and it effectively expands the analysis options available to end users.

Scalar connectivity measures from fast-marching tractography reveal heritability of white matter architecture

Recent advances in diffusion-weighted MRI (DWI) have enabled studies of complex white matter tissue architecture in vivo. To date, the underlying influence of genetic and environmental factors in determining central nervous system connectivity has not been widely studied. In this work, we introduce new scalar connectivity measures based on a computationally-efficient fast-marching algorithm for quantitative tractography. We then calculate connectivity maps for a DTI dataset from 92 healthy adult twins and decompose the genetic and environmental contributions to the variance in these metrics using structural equation models. By combining these techniques, we generate the first maps to directly examine genetic and environmental contributions to brain connectivity in humans. Our approach is capable of extracting statistically significant measures of genetic and environmental contributions to neural connectivity.

Mesh-based spherical deconvolution for physically valid fiber orientation reconstruction from diffusion-weighted MRI

High angular resolution diffusion imaging (HARDI) methods have enabled the reconstruction of complex spin diffusion profiles in central nervous system white matter through diffusion-weighted MRI. For recovery of the underlying fiber orientations, conventional spherical deconvolution techniques based on spherical harmonics typically have difficulty producing fiber orientation distributions (FODs) that simultaneously satisfy the physical constraints of being real, symmetric, and non-negative. In this work, we propose a novel approach for HARDI reconstruction that is guaranteed to generate FODs satisfying these constraints. By using a meshed representation of the unit sphere, we formulate the spherical deconvolution as a convex optimization problem and compute the solution using a projected gradient descent algorithm. Flexible regularization is also included in our method to allow for tuning the sharpness of the reconstructed FOD. In our experiments, we present simulated results to examine the effects of varying the regularization parameters, and we illustrate the robustness of our method by applying it to several biological data sets to reconstruct known white matter fiber geometry.

A 3-dimensional DTI MRI-based model of GBM growth and response to radiation therapy

Glioblastoma (GBM) is both the most common and the most aggressive intra-axial brain tumor, with a notoriously poor prognosis. To improve this prognosis, it is necessary to understand the dynamics of GBM growth, response to treatment and recurrence. The present study presents a mathematical diffusion-proliferation model of GBM growth and response to radiation therapy based on diffusion tensor (DTI) MRI imaging. This represents an important advance because it allows 3-dimensional tumor modeling in the anatomical context of the brain. Specifically, tumor infiltration is guided by the direction of the white matter tracts along which glioma cells infiltrate. This provides the potential to model different tumor growth patterns based on location within the brain, and to simulate the tumors response to different radiation therapy regimens. Tumor infiltration across the corpus callosum is simulated in biologically accurate time frames. The response to radiation therapy, including changes in cell density gradients and how these compare across different radiation fractionation protocols, can be rendered. Also, the model can estimate the amount of subthreshold tumor which has extended beyond the visible MR imaging margins. When combined with the ability of being able to estimate the biological parameters of invasiveness and proliferation of a particular GBM from serial MRI scans, it is shown that the model has potential to simulate realistic tumor growth, response and recurrence patterns in individual patients. To the best of our knowledge, this is the first presentation of a DTI-based GBM growth and radiation therapy treatment model.

A framework for secure and decentralized sharing of medical imaging data via blockchain consensus

The electronic sharing of medical imaging data is an important element of modern healthcare systems, but current infrastructure for cross-site image transfer depends on trust in third-party intermediaries. In this work, we examine the blockchain concept, which enables parties to establish consensus without relying on a central authority. We develop a framework for cross-domain image sharing that uses a blockchain as a distributed data store to establish a ledger of radiological studies and patient-defined access permissions. The blockchain framework is shown to eliminate third-party access to protected health information, satisfy many criteria of an interoperable health system, and readily generalize to domains beyond medical imaging. Relative drawbacks of the framework include the complexity of the privacy and security models and an unclear regulatory environment. Ultimately, the large-scale feasibility of such an approach remains to be demonstrated and will depend on a number of factors which we discuss in detail.

Image-driven modeling of the proliferation and necrosis of glioblastoma multiforme

BackgroundThe heterogeneity of response to treatment in patients with glioblastoma multiforme suggests that the optimal therapeutic approach incorporates an individualized assessment of expected lesion progression. In this work, we develop a novel computational model for the proliferation and necrosis of glioblastoma multiforme.MethodsThe model parameters are selected based on the magnetic resonance imaging features of each tumor, and the proposed technique accounts for intrinsic cell division, tumor cell migration along white matter tracts, as well as central tumor necrosis. As a validation of this approach, tumor growth is simulated in the brain of a healthy adult volunteer using parameters derived from the imaging of a patient with glioblastoma multiforme. A mutual information metric is calculated between the simulated tumor profile and observed tumor.ResultsThe tumor progression profile generated by the proposed model is compared with those produced by existing models and with the actual observed tumor progression. Both qualitative and quantitative analyses show that the model introduced in this work replicates the observed progression of glioblastoma more accurately relative to prior techniques.ConclusionsThis image-driven model generates improved tumor progression profiles and may contribute to the development of more reliable prognostic estimates in patients with glioblastoma multiforme.