Your address will show here +12 34 56 78

Because it is sensitive to the microscopic displacements of water molecules diffusing in the tissues, diffusion-weighted magnetic resonance imaging (DW-MRI) has become the modality of choice to study the anatomy of the white matter at the micrometer level, with applications in the diagnosis and treatment of many neurological disorders and brain injuries. Our new approach at the intersection of physical modeling, numerical simulations and signal processing proposes a way to bridge the gap between the native millimeter resolution of MRI and neuroanatomical indices at the micrometer scale.

We use specific MRI acquisitions known as diffusion-weighted MRI which by design are sensitive to the micrometer-scale diffusing motion of water molecules in the tissues at body temperature. In traditional approaches, forward mathematical models based on analytical formulas are devised to link the geometry of white matter tissue, the diffusion of water molecules in that geometric arrangement and the physics of MRI acquisitions. The free parameters of such a forward model are then adjusted to new MRI measurements by minimizing a cost function, providing indices of the white matter microstructure locally in each voxel. Such indices usually reflect the morphometry of the axons and glial cells making up the white matter.
In our latest work, we presented a forward model in which the MRI signal is not modeled analytically but instead computed by numerical simulations of the random walk of water in a 3D mesh representing the tissue. Such simulations, known as Monte Carlo simulations, provide signals that are more physically accurate than most traditional closed-form analytical models used in the literature. Due to their long computation times, their use has been mainly restricted to validation rather than directly in a model of the signal. In addition, being numerical results rather than mathematical formulas, they do not easily lend themselves to continuous optimization frameworks for parameter estimation. To circumvent these difficulties, we made use of combinatorial and convex optimization on a presimulated collection of MRI signals, or fingerprints, wherein each fingerprint uniquely characterizes a typical microstructural configuration.
On a rat model of spinal cord injury, we showed that the microstructural indices obtained with our model based on in vivo, non-invasive MRI were consistent with histological observations performed after sacrificing the animal. In particular, our model was able to detect axonal loss in the injured part of the spinal cord. The method was then applied on in vivo human brain data. Indices of axon diameter and density in specific white matter pathways were found to agree with previous histological observations in the human brain.
Future work will focus on using 3D geometries which more closely mimic white matter tissue in our presimulated dictionary of fingerprints. The method will also be applied on existing clinical datasets. This work encourages the use of computational models of the diffusion signal and gets us a step closer to quantitative characterization of the white matter microstructure in vivo.
Signal Processing Lab (LTS5), Ecole polytechnique fédérale de Lausanne, Switzerland
Computational Radiology Laboratory, Boston Children’s Hospital, Harvard Medical School, Boston, USA
Institute of NeuroScience, Université catholique de Louvain (Mont-Godinne), Belgium

Twitter (@MaximeTaquet):

Medical Imaging
ImagX.R, the research group on image guided proton therapy, has the pleasure to announce the PhD defence of Sylvain Deffet. He will present his research on proton radiography and discuss how this method could improve treatment outcome in proton therapy.

Proton therapy is an advanced form of radiation therapy which is increasingly used worldwide. Unlike photons, protons deliver a sharp dose at a precise location corresponding to their range in the patient. Thanks to this physical property, proton therapy has the potential to spare healthy tissues better than conventional radiation therapy. However, a consequence of this dosimetric property is that the range of the protons inside the patient must be accurately predicted to deliver the dose as planned. Unfortunately, several uncertainties arising during treatment planning may significantly impact the range of the protons and hence jeopardize dose conformity. To better quantify and potentially reduce the uncertainties, the implementation of imaging techniques that would provide a direct information on the energy reduction of protons in the patient is highly desirable. In the present thesis, such a method is introduced: proton radiography. The system that we propose relies on the use of protons having an energy high enough for them to traverse the patient and stop inside a detector which measures their residual range. By comparing measured and predicted residual ranges, it is possible to quantify range uncertainties in the patient in clinical conditions.

Although the principle of proton radiography is relatively straightforward, a clinical implementation of such a measurement device is a complex issue. This thesis aims at developing a proton radiography system and also at conceiving an acquisition process that could be conveniently implemented in clinics. The methodologies to correctly take advantage of proton radiography measurements are also discussed, in order to further optimize the dose delivered to the patient and improve de facto the treatment outcome.