Joint radial trajectory correction for accelerated T2 * mapping on an MR-Linac

BACKGROUND

T2 * mapping can characterize tumor hypoxia, which may be associated with resistance to therapy. Acquiring T2 * maps during MR-guided radiotherapy could inform treatment adaptation by, for example, escalating the dose to resistant sub-volumes.

PURPOSE

The purpose of this work is to demonstrate the feasibility of the accelerated T2 * mapping technique using model-based image reconstruction with integrated trajectory auto-correction (TrACR) for MR-guided radiotherapy on an MR-Linear accelerator (MR-Linac).

MATERIALS AND METHODS

The proposed method was validated in a numerical phantom, where two T2 * mapping approaches (sequential and joint) were compared for different noise levels (0,0.1,0.5,1) and gradient delays ([1, -1] and [1, -2] in units of dwell time for x- and y-axis, respectively). Fully sampled k-space was retrospectively undersampled using two different undersampling patterns. Root mean square errors (RMSEs) were calculated between reconstructed T2 * maps and ground truth. In vivo data was acquired twice weekly in one prostate and one head and neck cancer patient undergoing treatment on a 1.5 T MR-Linac. Data were retrospectively undersampled and T2 * maps reconstructed, with and without trajectory corrections were compared.

RESULTS

Numerical simulations demonstrated that, for all noise levels, T2 * maps reconstructed with a joint approach demonstrated less error compared to an uncorrected and sequential approach. For a noise level of 0.1, uniform undersampling and gradient delay [1, -1] (in units of dwell time for x- and y-axis, respectively), RMSEs for sequential and joint approaches were 13.01 and 9.32 ms, respectively, which reduced to 10.92 and 5.89 ms for a gradient delay of [1, 2]. Similarly, for alternate undersampling and gradient delay [1, -1], RMSEs for sequential and joint approaches were 9.80 and 8.90 ms, respectively, which reduced to 9.10 and 5.40 ms for gradient delay [1, 2]. For in vivo data, T2 * maps reconstructed with our proposed approach resulted in less artifacts and improved visual appearance compared to the uncorrected approach. For both prostate and head and neck cancer patients, T2 * maps reconstructed from different treatment fractions showed changes within the planning target volume (PTV).

CONCLUSION

Using the proposed approach, a retrospective data-driven gradient delay correction can be performed, which is particularly relevant for hybrid devices, where full information on the machine configuration is not available for image reconstruction. T2 * maps were acquired in under 5 min and can be integrated into MR-guided radiotherapy treatment workflows, which minimizes patient burden and leaves time for additional imaging for online adaptive radiotherapy on an MR-Linac.

EEG-based BCI Dataset of Semantic Concepts for Imagination and Perception Tasks

Electroencephalography (EEG) is a widely-used neuroimaging technique in Brain Computer Interfaces (BCIs) due to its non-invasive nature, accessibility and high temporal resolution. A range of input representations has been explored for BCIs. The same semantic meaning can be conveyed in different representations, such as visual (orthographic and pictorial) and auditory (spoken words). These stimuli representations can be either imagined or perceived by the BCI user. In particular, there is a scarcity of existing open source EEG datasets for imagined visual content, and to our knowledge there are no open source EEG datasets for semantics captured through multiple sensory modalities for both perceived and imagined content. Here we present an open source multisensory imagination and perception dataset, with twelve participants, acquired with a 124 EEG channel system. The aim is for the dataset to be open for purposes such as BCI related decoding and for better understanding the neural mechanisms behind perception, imagination and across the sensory modalities when the semantic category is held constant.

Rapid three-dimensional multiparametric MRI with quantitative transient-state imaging

Novel methods for quantitative, transient-state multiparametric imaging are increasingly being demonstrated for assessment of disease and treatment efficacy. Here, we build on these by assessing the most common Non-Cartesian readout trajectories (2D/3D radials and spirals), demonstrating efficient anti-aliasing with a k-space view-sharing technique, and proposing novel methods for parameter inference with neural networks that incorporate the estimation of proton density. Our results show good agreement with gold standard and phantom references for all readout trajectories at 1.5 T and 3 T. Parameters inferred with the neural network were within 6.58% difference from the parameters inferred with a high-resolution dictionary. Concordance correlation coefficients were above 0.92 and the normalized root mean squared error ranged between 4.2 and 12.7% with respect to gold-standard phantom references for T1 and T2. In vivo acquisitions demonstrate sub-millimetric isotropic resolution in under five minutes with reconstruction and inference times < 7 min. Our 3D quantitative transient-state imaging approach could enable high-resolution multiparametric tissue quantification within clinically acceptable acquisition and reconstruction times.

Model-based super-resolution reconstruction of T2 maps

PURPOSE

High-resolution isotropic T₂ mapping of the human brain with multi-echo spin-echo (MESE) acquisitions is challenging. When using a 2D sequence, the resolution is limited by the slice thickness. If used as a 3D acquisition, specific absorption rate limits are easily exceeded due to the high power deposition of nonselective refocusing pulses. A method to reconstruct 1-mm³ isotropic T₂ maps is proposed based on multiple 2D MESE acquisitions. Data were undersampled (10-fold) to compensate for the prolonged scan time stemming from the super-resolution acquisition.

THEORY AND METHODS

The proposed method integrates a classical super-resolution with an iterative model-based approach to reconstruct quantitative maps from a set of undersampled low-resolution data. The method was tested on numerical and multipurpose phantoms, and in vivo data. T₂ values were assessed with a region-of-interest analysis using a single-slice spin-echo and a fully sampled MESE acquisition in a phantom, and a MESE acquisition in healthy volunteers.

RESULTS

Numerical simulations showed that the best trade-off between acceleration and number of low-resolution datasets is 10-fold acceleration with 4 acquisitions (acquisition time = 18 min). The proposed approach showed improved resolution over low-resolution images for both phantom and brain. Region-of-interest analysis of the phantom compartments revealed that at shorter T₂ , the proposed method was comparable with the fully sampled MESE. For the volunteer data, the T₂ values found in the brain structures were consistent across subjects (8.5-13.1 ms standard deviation).

CONCLUSION

The proposed method addresses the inherent limitations associated with high-resolution T₂ mapping and enables the reconstruction of 1 mm³ isotropic relaxation maps with a 10 times faster acquisition.

Compressive source separation: theory and methods for hyperspectral imaging

We propose and analyze a new model for hyperspectral images (HSIs) based on the assumption that the whole signal is composed of a linear combination of few sources, each of which has a specific spectral signature, and that the spatial abundance maps of these sources are themselves piecewise smooth and therefore efficiently encoded via typical sparse models. We derive new sampling schemes exploiting this assumption and give theoretical lower bounds on the number of measurements required to reconstruct HSI data and recover their source model parameters. This allows us to segment HSIs into their source abundance maps directly from compressed measurements. We also propose efficient optimization algorithms and perform extensive experimentation on synthetic and real datasets, which reveals that our approach can be used to encode HSI with far less measurements and computational effort than traditional compressive sensing methods.