The effect of prolonged spaceflight on cerebrospinal fluid and perivascular spaces of astronauts and cosmonauts.

Long-duration spaceflight induces changes to the brain and cerebrospinal fluid compartments and visual acuity problems known as spaceflight-associated neuro-ocular syndrome (SANS). The clinical relevance of these changes and whether they equally affect crews of different space agencies remain unknown. We used MRI to analyze the alterations occurring in the perivascular spaces (PVS) in NASA and European Space Agency astronauts and Roscosmos cosmonauts after a 6-mo spaceflight on the International Space Station (ISS). We found increased volume of basal ganglia PVS and white matter PVS (WM-PVS) after spaceflight, which was more prominent in the NASA crew than the Roscosmos crew. Moreover, both crews demonstrated a similar degree of lateral ventricle enlargement and decreased subarachnoid space at the vertex, which was correlated with WM-PVS enlargement. As all crews experienced the same environment aboard the ISS, the differences in WM-PVS enlargement may have been due to, among other factors, differences in the use of countermeasures and high-resistive exercise regimes, which can influence brain fluid redistribution. Moreover, NASA astronauts who developed SANS had greater pre- and postflight WM-PVS volumes than those unaffected. These results provide evidence for a potential link between WM-PVS fluid and SANS.

Model-based super-resolution reconstruction for pseudo-continuous Arterial Spin Labeling.

Arterial spin labeling (ASL) is a promising, non-invasive perfusion magnetic resonance imaging technique for quantifying cerebral blood flow (CBF). Unfortunately, ASL suffers from an inherently low signal-to-noise ratio (SNR) and spatial resolution, undermining its potential. Increasing spatial resolution without significantly sacrificing SNR or scan time represents a critical challenge towards routine clinical use. In this work, we propose a model-based super-resolution reconstruction (SRR) method with joint motion estimation that breaks the traditional SNR/resolution/scan-time trade-off. From a set of differently oriented 2D multi-slice pseudo-continuous ASL images with a low through-plane resolution, 3D-isotropic, high resolution, quantitative CBF maps are estimated using a Bayesian approach. Experiments on both synthetic whole brain phantom data, and on in vivo brain data, show that the proposed SRR Bayesian estimation framework outperforms state-of-the-art ASL quantification.

Edge illumination x-ray phase contrast simulations using the CAD-ASTRA toolbox.

Edge illumination x-ray phase contrast imaging (XPCI) provides increased contrast for low absorbing materials compared to attenuation images and sheds light on the material microstructure through dark field contrast. To apply XPCI in areas such as non-destructive testing and inline inspection, where scanned samples are increasingly compared to simulated reference images, accurate and efficient simulation software is required. However, currently available simulators rely on expensive Monte Carlo techniques or wave-optics frameworks, resulting in long simulation times. Furthermore, these simulators are often not optimized to work with computer-aided design (CAD) models, a common and memory-efficient method to represent manufactured objects, hindering their integration in an inspection pipeline. In this work, we address these shortcomings by introducing an edge illumination XPCI simulation framework built upon the recently developed CAD-ASTRA toolbox. CAD-ASTRA allows for the efficient simulation of x-ray projections from CAD models through GPU-accelerated ray tracing and supports ray refraction in a geometric optics framework. The edge illumination implementation is validated and its performance is benchmarked against GATE, a state-of-the-art Monte Carlo simulator, revealing a simulation speed increase of up to three orders of magnitude, while maintaining high accuracy in the resulting images.

Optimal experimental design and estimation for q-space trajectory imaging.

Tensor-valued diffusion encoding facilitates data analysis by q-space trajectory imaging. By modeling the diffusion signal of heterogeneous tissues with a diffusion tensor distribution (DTD) and modulating the encoding tensor shape, this novel approach allows disentangling variations in diffusivity from microscopic anisotropy, orientation dispersion, and mixtures of multiple isotropic diffusivities. To facilitate the estimation of the DTD parameters, a parsimonious acquisition scheme coupled with an accurate and precise estimation of the DTD is needed. In this work, we create two precision-optimized acquisition schemes: one that maximizes the precision of the raw DTD parameters, and another that maximizes the precision of the scalar measures derived from the DTD. The improved precision of these schemes compared to a naïve sampling scheme is demonstrated in both simulations and real data. Furthermore, we show that the weighted linear least squares (WLLS) estimator that uses the squared reciprocal of the noisy signal as weights can be biased, whereas the iteratively WLLS estimator with the squared reciprocal of the predicted signal as weights outperforms the conventional unweighted linear LS and nonlinear LS estimators in terms of accuracy and precision. Finally, we show that the use of appropriate constraints can considerably increase the precision of the estimator with only a limited decrease in accuracy.

Systematic review of reconstruction techniques for accelerated quantitative MRI.

PURPOSE: To systematically review the techniques that address undersampling artifacts in accelerated quantitative magnetic resonance imaging (qMRI). METHODS: A literature search was conducted using the Embase, Medline, Web of Science Core Collection, Coherence Central Register of Controlled Trials, and Google Scholar databases for studies, published before July 2022 proposing reconstruction techniques for accelerated qMRI. Studies are reviewed according to inclusion criteria, and included studies are categorized based on the methodology used. RESULTS: A total of 292 studies included in the review are categorized. A technical overview of each category is provided, and the categories are described in a unified mathematical framework. The distribution of the reviewed studies over time, application domain, and parameters of interest is illustrated. CONCLUSION: An increasing trend in the number of articles that propose new techniques for accelerated qMRI reconstruction indicates the importance of acceleration in qMRI. The techniques are mostly validated for relaxometry parameters and brain scans. The categories of techniques are compared based on theoretical grounds, highlighting existing trends and potential gaps in the field.

ADEPT: Accurate Diffusion Echo-Planar imaging with multi-contrast shoTs.

PURPOSE: To introduce a novel imaging and parameter estimation framework for accurate multi-shot diffusion MRI. THEORY AND METHODS: We propose a new framework called ADEPT (Accurate Diffusion Echo-Planar imaging with multi-contrast shoTs) that enables fast diffusion MRI by allowing diffusion contrast settings to change between shots in a multi-shot EPI acquisition (i.e., intra-scan modulation). The framework estimates diffusion parameter maps directly from the acquired intra-scan modulated k-space data, while simultaneously accounting for shot-to-shot phase inconsistencies. The performance of the estimation framework is evaluated using Monte Carlo simulation studies and in-vivo experiments and compared to that of reference methods that rely on parallel imaging for shot-to-shot phase correction. RESULTS: Simulation and real-data experiments show that ADEPT yields more accurate and more precise estimates of the diffusion metrics in multi-shot EPI data in comparison with the reference methods. CONCLUSION: ADEPT allows fast multi-shot EPI diffusion MRI without significantly degrading the accuracy and precision of the estimated diffusion maps.

Grating designs for cone beam edge illumination X-ray phase contrast imaging: a simulation study.

Edge illumination is an emerging X-ray phase contrast imaging technique providing attenuation, phase and dark field contrast. Despite the successful transition from synchrotron to lab sources, the cone beam geometry of lab systems limits the effectiveness of using conventional planar gratings. The non-parallel incidence of X-rays introduces shadowing effects, worsening with increasing cone angle. To overcome this limitation, several alternative grating designs can be considered. In this paper, the effectiveness of three alternative designs is compared to conventional gratings using numerical simulations. Improvements in flux and contrast are discussed, taking into account practical considerations concerning the implementation of the designs.

Efficient iterative reconstruction with beam shape compensation for THz computed tomography.

Terahertz (THz) computed tomography is an emerging nondestructive and non-ionizing imaging method. Most THz reconstruction methods rely on the Radon transform, originating from x-ray imaging, in which x rays propagate in straight lines. However, a THz beam has a finite width, and ignoring its shape results in blurred reconstructed images. Moreover, accounting for the THz beam model in a straightforward way in an iterative reconstruction method results in extreme demands in memory and in slow convergence. In this paper, we propose an efficient iterative reconstruction that incorporates the THz beam shape, while avoiding the above disadvantages. Both simulation and real experiments show that our approach results in improved resolution recovery in the reconstructed image. Furthermore, we propose a suitable preconditioner to improve the convergence speed of our reconstruction.

Improved diffusion parameter estimation by incorporating T2 relaxation properties into the DKI-FWE model.

The free water elimination (FWE) model and its kurtosis variant (DKI-FWE) can separate tissue and free water signal contributions, thus providing tissue-specific diffusional information. However, a downside of these models is that the associated parameter estimation problem is ill-conditioned, necessitating the use of advanced estimation techniques that can potentially bias the parameter estimates. In this work, we propose the T2-DKI-FWE model that exploits the T2 relaxation properties of both compartments, thereby better conditioning the parameter estimation problem and providing, at the same time, an additional potential biomarker (the T2 of tissue). In our approach, the T2 of tissue is estimated as an unknown parameter, whereas the T2 of free water is assumed known a priori and fixed to a literature value (1573 ms). First, the error propagation of an erroneous assumption on the T2 of free water is studied. Next, the improved conditioning of T2-DKI-FWE compared to DKI-FWE is illustrated using the Cramér-Rao lower bound matrix. Finally, the performance of the T2-DKI-FWE model is compared to that of the DKI-FWE and T2-DKI models on both simulated and real datasets. The error due to a biased approximation of the T2 of free water was found to be relatively small in various diffusion metrics and for a broad range of erroneous assumptions on its underlying ground truth value. Compared to DKI-FWE, using the T2-DKI-FWE model is beneficial for the identifiability of the model parameters. Our results suggest that the T2-DKI-FWE model can achieve precise and accurate diffusion parameter estimates, through effective reduction of free water partial volume effects and by using a standard nonlinear least squares approach. In conclusion, incorporating T2 relaxation properties into the DKI-FWE model improves the conditioning of the model fitting, while only requiring an acquisition scheme with at least two different echo times.

Virtual grating approach for Monte Carlo simulations of edge illumination-based x-ray phase contrast imaging.

The design of new x-ray phase contrast imaging setups often relies on Monte Carlo simulations for prospective parameter studies. Monte Carlo simulations are known to be accurate but time consuming, leading to long simulation times, especially when many parameter variations are required. This is certainly the case for imaging methods relying on absorbing masks or gratings, with various tunable properties, such as pitch, aperture size, and thickness. In this work, we present the virtual grating approach to overcome this limitation. By replacing the gratings in the simulation with virtual gratings, the parameters of the gratings can be changed after the simulation, thereby significantly reducing the overall simulation time. The method is validated by comparison to explicit grating simulations, followed by representative demonstration cases.

Brain ventricular volume changes induced by long-duration spaceflight.

Long-duration spaceflight induces detrimental changes in human physiology. Its residual effects and mechanisms remain unclear. We prospectively investigated the changes in cerebrospinal fluid (CSF) volume of the brain ventricular regions in space crew by means of a region of interest analysis on structural brain scans. Cosmonaut MRI data were investigated preflight (n = 11), postflight (n = 11), and at long-term follow-up 7 mo after landing (n = 7). Post hoc analyses revealed a significant difference between preflight and postflight values for all supratentorial ventricular structures, i.e., lateral ventricle (mean % change ± SE = 13.3 ± 1.9), third ventricle (mean % change ± SE = 10.4 ± 1.1), and the total ventricular volume (mean % change ± SE = 11.6 ± 1.5) (all P < 0.0001), with higher volumes at postflight. At follow-up, these structures did not quite reach baseline levels, with still residual increases in volume for the lateral ventricle (mean % change ± SE = 7.7 ± 1.6; P = 0.0009), the third ventricle (mean % change ± SE = 4.7 ± 1.3; P = 0.0063), and the total ventricular volume (mean % change ± SE = 6.4 ± 1.3; P = 0.0008). This spatiotemporal pattern of CSF compartment enlargement and recovery points to a reduced CSF resorption in microgravity as the underlying cause. Our results warrant more detailed and longer longitudinal follow-up. The clinical impact of our findings on the long-term cosmonauts' health and their relation to ocular changes reported in space travelers requires further prospective studies.

A Bottom-Up Volume Reconstruction Method for Atom Probe Tomography.

This paper describes a reconstruction method for atom probe tomography based on a bottom-up approach accounting for (i) the final tip morphology (which is frequently induced by inhomogeneous evaporation probabilities across the tip surface due to laser absorption, heat diffusion effects, and inhomogeneous material properties), (ii) the limited (and changing) field of view, and (iii) the detector efficiency. The reconstruction starts from the final tip morphology and reverses the evaporation sequence through the pseudo-deposition of defined small reconstruction volumes, which are then stacked together to create the full three-dimensional (3D) tip. The subdivision in small reconstruction volumes allows the scheme to account for the changing tip shape and field of view as evaporation proceeds. Atoms within the same small reconstruction volume are reconstructed at once by placing atoms back onto their possible lattice sites through a trajectory-matching process involving simulated and experimental hit maps. As the ejected ion trajectories are simulated using detailed electrostatic modeling inside the chamber, no simplifications have been imposed on the shape of the trajectories, projection laws, or tip surface. We demonstrate the superior performance of our approach over the conventional reconstruction method (Bas) for an asymmetrical tip shape.

On the generalizability of diffusion MRI signal representations across acquisition parameters, sequences and tissue types: Chronicles of the MEMENTO challenge.

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

Constrained spherical deconvolution of nonspherically sampled diffusion MRI data.

Constrained spherical deconvolution (CSD) of diffusion-weighted MRI (DW-MRI) is a popular analysis method that extracts the full white matter (WM) fiber orientation density function (fODF) in the living human brain, noninvasively. It assumes that the DW-MRI signal on the sphere can be represented as the spherical convolution of a single-fiber response function (RF) and the fODF, and recovers the fODF through the inverse operation. CSD approaches typically require that the DW-MRI data is sampled shell-wise, and estimate the RF in a purely spherical manner using spherical basis functions, such as spherical harmonics (SH), disregarding any radial dependencies. This precludes analysis of data acquired with nonspherical sampling schemes, for example, Cartesian sampling. Additionally, nonspherical sampling can also arise due to technical issues, for example, gradient nonlinearities, resulting in a spatially dependent bias of the apparent tissue densities and connectivity information. Here, we adopt a compact model for the RFs that also describes their radial dependency. We demonstrate that the proposed model can accurately predict the tissue response for a wide range of b-values. On shell-wise data, our approach provides fODFs and tissue densities indistinguishable from those estimated using SH. On Cartesian data, fODF estimates and apparent tissue densities are on par with those obtained from shell-wise data, significantly broadening the range of data sets that can be analyzed using CSD. In addition, gradient nonlinearities can be accounted for using the proposed model, resulting in much more accurate apparent tissue densities and connectivity metrics.

Accelerating in vivo fast spin echo high angular resolution diffusion imaging with an isotropic resolution in mice through compressed sensing.

PURPOSE: Echo planar imaging (EPI) is commonly used to acquire the many volumes needed for high angular resolution diffusion Imaging (HARDI), posing a higher risk for artifacts, such as distortion and deformation. An alternative to EPI is fast spin echo (FSE) imaging, which has fewer artifacts but is inherently slower. The aim is to accelerate FSE such that a HARDI data set can be acquired in a time comparable to EPI using compressed sensing. METHODS: Compressed sensing was applied in either q-space or simultaneously in k-space and q-space, by undersampling the k-space in the phase-encoding direction or retrospectively eliminating diffusion directions for different degrees of undersampling. To test the replicability of the acquisition and reconstruction, brain data were acquired from six mice, and a numerical phantom experiment was performed. All HARDI data were analyzed individually using constrained spherical deconvolution, and the apparent fiber density and complexity metric were evaluated, together with whole-brain tractography. RESULTS: The apparent fiber density and complexity metric showed relatively minor differences when only q-space undersampling was used, but deteriorate when k-space undersampling was applied. Likewise, the tract density weighted image showed good results when only q-space undersampling was applied using 15 directions or more, but information was lost when fewer volumes or k-space undersampling were used. CONCLUSION: It was found that acquiring 15 to 20 diffusion directions with a full k-space and reconstructed using compressed sensing could suffice for a replicable measurement of quantitative measures in mice, where areas near the sinuses and ear cavities are untainted by signal loss.

FleXCT: a flexible X-ray CT scanner with 10 degrees of freedom.

Laboratory based X-ray micro-CT is a non-destructive testing method that enables three dimensional visualization and analysis of the internal and external morphology of samples. Although a wide variety of commercial scanners exist, most of them are limited in the number of degrees of freedom to position the source and detector with respect to the object to be scanned. Hence, they are less suited for industrial X-ray imaging settings that require advanced scanning modes, such as laminography, conveyor belt scanning, or time-resolved imaging (4DCT). We introduce a new X-ray scanner FleXCT that consists of a total of ten motorized axes, which allow a wide range of non-standard XCT scans such as tiled and off-centre scans, laminography, helical tomography, conveyor belt, dynamic zooming, and X-ray phase contrast imaging. Additionally, a new software tool 'FlexRayTools' was created that enables reconstruction of non-standard XCT projection data of the FleXCT instrument using the ASTRA Toolbox, a highly efficient and open source set of tools for tomographic projection and reconstruction.

Small medial femoral condyle morphotype is associated with medial compartment degeneration and distinct morphological characteristics: a comparative pilot study.

PURPOSE: Early-onset degeneration of the knee is linked to genetics, overload, injury, and potentially, knee morphology. The purpose of this study is to explore the characteristics of the small medial femoral condyle, as a distinct knee morphotype, by means of a landmark-based three-dimensional (3D) analysis and statistical parametric mapping. METHODS: Sixteen knees with a small medial femoral condyle (SMC) were selected from a database of patients with distinct knee joint anatomy and 16 gender-matched knees were selected from a control group database. 3D models were generated from the medical imaging. After normalization for size, a set of pre-defined landmark-based parameters was analysed for the femur and tibia. Local shape differences were evaluated by matching all bone surfaces onto each other and comparing the distances to the mean control group bone shape. RESULTS: The small medial condyle group showed a significant association with medial compartment degeneration and had a 4% and 13% smaller medial condyle anteroposteriorly and mediolaterally, whereas the distal femur was 3% wider mediolaterally. The lateral condyle was 2% smaller anteroposteriorly and 8% wider mediolaterally. The complete tibial plateau was 3% smaller mediolaterally and the medial tibial plateau was 6% smaller. CONCLUSION: A new knee morphotype demonstrated an increased risk for medial compartment degeneration and was differentiated from a healthy control group based on the following morphological characteristics: a smaller medial femoral condyle and medial tibial plateau, a wider lateral femoral condyle and a wider distal femur on a smaller tibial plateau. This pilot study suggests a role for the SMC knee morphotype in the multifactorial process of medial compartment degeneration. LEVEL OF EVIDENCE: III.

Fracture patterns in midshaft clavicle fractures.

Current classifications of midshaft clavicle fractures are based on radiography. The aim of the study was to evaluate the fracture pattern of clavicle fractures using 3-dimensional computed tomography (3D CT). A retrospective analysis was performed on CT scans of 65 acute clavicle fractures. Using quantitative 3D CT reconstruction techniques, the fracture of the clavicle was virtually reduced. Based on these reconstructions, a group-based fracture heat map and small fragment heat map, and the location of the most common fracture line were determined. Also, the direction and amount of displacement were evaluated. Three fracture patterns could be distinguished. The primary fracture line in type 1 is going from posteromedial to anterolateral and located between 50% and 68% of the clavicle's length. In type 2, a transverse fracture line is located around 55%, and in type 3, a superolateral to inferomedial line is located between 47% and 56%. Wedged fracture fragments can be seen in types 1 and 2 and are mainly situated inferiorly. The displacement is similar in all types, but the main direction of displacement is specific for the different types (posterior, anterior, inferior). We can conclude that several fracture patterns can be seen in clavicle fractures. Most fractures are located laterally at the midshaft of the clavicle. Wedged segments are mainly located inferiorly, and at the posterior part of the clavicle, no comminution is ever seen. The direction of displacement depends on the fracture pattern.

Supporting measurements or more averages? How to quantify cerebral blood flow most reliably in 5 minutes by arterial spin labeling.

PURPOSE: To determine whether sacrificing part of the scan time of pseudo-continuous arterial spin labeling (PCASL) for measurement of the labeling efficiency and blood T1 is beneficial in terms of CBF quantification reliability. METHODS: In a simulation framework, 5-minute scan protocols with different scan time divisions between PCASL data acquisition and supporting measurements were evaluated in terms of CBF estimation variability across both noise and ground truth parameter realizations taken from the general population distribution. The entire simulation experiment was repeated for a single-post-labeling delay (PLD), multi-PLD, and free-lunch time-encoded (te-FL) PCASL acquisition strategy. Furthermore, a real data study was designed for preliminary validation. RESULTS: For the considered population statistics, measuring the labeling efficiency and the blood T1 proved beneficial in terms of CBF estimation variability for any distribution of the 5-minute scan time compared to only acquiring ASL data. Compared to single-PLD PCASL without support measurements as recommended in the consensus statement, a 26%, 33%, and 42% reduction in relative CBF estimation variability was found for optimal combinations of supporting measurements with single-PLD, free-lunch, and multi-PLD PCASL data acquisition, respectively. The benefit of taking the individual variation of blood T1 into account was also demonstrated in the real data experiment. CONCLUSIONS: Spending time to measure the labeling efficiency and the blood T1 instead of acquiring more averages of the PCASL data proves to be advisable for robust CBF quantification in the general population.

X-ray phase contrast simulation for grating-based interferometry using GATE.

The overall importance of x-ray phase contrast (XPC) imaging has grown substantially in the last decades, in particular with the recent advent of compact lab-based XPC systems. For optimizing the experimental XPC setup, as well as benchmarking and testing new acquisition and reconstruction techniques, Monte Carlo (MC) simulations are a valuable tool. GATE, an open source application layer on top of the Geant4 simulation software, is a versatile MC tool primarily intended for various types of medical imaging simulations. To our knowledge, however, there is no GATE-based academic simulation software available for XPC imaging. In this paper, we extend the GATE framework with new physics-based tools for accurate XPC simulations. Our approach combines Monte Carlo simulations in GATE for modelling the x-ray interactions in the sample with subsequent numerical wave propagation, starting from the GATE output.