Organizations: Juntendo University, Tokyo, Japan · RIKEN Center for Advanced Intelligence Project, Tokyo, Japan · Tokyo University of Agriculture and Technology, Tokyo, Japan
Multi-shell diffusion MRI (dMRI) unlocks more expressive microstructural modeling than single-shell scans, yet its longer acquisition time hinders deployment in large-scale cohorts and time-constrained clinical settings. Synthesizing an unobserved shell from a single-shell input is fundamentally ill-posed and further complicated by protocol mismatch, where source and target gradient direction sets may not align. We propose DTI-SHNet, a single-to-multi-shell synthesis framework that operates in the real symmetric spherical harmonics (SH) coefficient domain and performs spatially aware volumetric regression. Given a source shell, we estimate diffusion tensor imaging (DTI) and use direction-agnostic parametric maps along with a brain mask as conditioning priors to guide a 3D U-Net regressor from source-shell to target-shell SH coefficients. To couple coefficient accuracy with signal fidelity, we introduce a signal consistency regularization that reconstructs signals on randomly sampled canonical directions from predicted coefficients and enforces agreement in the signal domain. Experiments on UK Biobank and Cam-CAN data for b=1000 to b=2000 dMRI synthesis show that DTI-SHNet achieves competitive visual quality compared to advanced methods, while better preserving downstream diffusion measures. Our code is available at https://github.com/xiaovhua/dti-shnet.
Figures & tables
Figure 1: Overview of our synthesis pipeline. Single-shell data at b=b1 are represented as real symmetric SH coefficient volumes Cb1 . DTI-SHNet predicts target coefficient volumes Cb2 guided by DTI priors. Target-shell signals are synthesized for desired directions by SH reconstruction, with an additional signal consistency constraint Lsig applied on canonical directions sampled during training.
Figure 2: Framework of DTI-SHNet. A 3D U-Net regresses target-shell SH coefficients from source-shell coefficients, conditioned on DTI priors.
Figure 3: Qualitative comparison for b=2000 synthesis from b=1000 on Cam-CAN. Representative slices illustrate that DTI-SHNet better preserves anatomical structures and reduces artifacts compared with competing baselines.
Method
UKB
Cam-CAN
PSNR ↑
SSIM ↑
NMSE ↓
PSNR ↑
SSIM ↑
NMSE ↓
Extrapolation [ 2 ]
21.9±0.6
0.825±0.014
0.180±0.028
22.0±0.5
0.805±0.014
0.134±0.023
DDPM [ 5 ]
21.5±0.8
0.824±0.022
0.229±0.030
24.0±0.8
0.883±0.023
0.109±0.021
SH-DNN [ 11 ]
25.6±0.7
0.930±0.007
0.079±0.016
26.7±0.5
0.951±0.005
0.045±0.007
DirUNet [ 4 ]
24.1±0.6
0.895±0.010
0.109±0.022
29.0±0.6
0.968±0.004
0.026±0.004
DTI-SHNet
25.3±0.7
0.922±0.008
0.083±0.018
29.3±0.6
0.970±0.004
0.025±0.004
Table 1: Main results on UKB and Cam-CAN for synthesizing b=2000 from b=1000 .
Condition
UKB
Cam-CAN
Priors
Lsig
PSNR ↑
SSIM ↑
NMSE ↓
Δ FA ↓
Δ MD ↓
PSNR ↑
SSIM ↑
NMSE ↓
Δ FA ↓
Δ MD ↓
×
×
25.33
0.9218
0.0829
0.0235
0.337
27.90
0.9547
0.0341
0.0604
0.903
✓
×
25.42
0.9235
0.0814
0.0224
0.322
28.64
0.9627
0.0288
0.0418
0.656
×
✓
25.29
0.9219
0.0838
0.0233
0.282
28.61
0.9639
0.0290
0.0343
0.467
✓
✓
25.33
0.9216
0.0830
0.0220
0.281
29.32
0.9700
0.0246
0.0319
0.426
Table 3: Ablation study on the DTI-derived priors and the signal consistency loss Lsig for synthesizing b=2000 from b=1000 . ΔMD is scaled by 10−4 .
Figure 4: Qualitative NODDI comparison between ground-truth multi-shell dMRI and synthesized multi-shell dMRI. We show representative maps of the NDI, ODI, and ISO, together with their corresponding histograms.
Diffusion models have achieved remarkable quality in multi-modal MRI synthesis, but their computational cost (hundreds of sampling steps and separate models per modality) limits clinical deployment. We observe that this inefficiency stems from an unnecessary starting point: diffusion begins from pure noise, discarding the structural information already present in available MRI sequences. We propose WFM (Wavelet Flow Matching), which instead learns a direct flow from an informed prior, the mean of conditioning modalities in wavelet space, to the target distribution. Because the source and target share underlying anatomy and differ primarily in contrast, this formulation enables accurate synthesis in just 1-2 integration steps. A single 82M-parameter model with class conditioning synthesizes all four BraTS modalities (T1, T1c, T2, FLAIR), replacing four separate diffusion models totaling 326M parameters. On BraTS 2024, WFM achieves 26.8 dB PSNR and 0.94 SSIM, within 1-2 dB of diffusion baselines, while running 250-1000x faster (0.16-0.64s vs. 160s per volume). This speed-quality trade-off makes real-time MRI synthesis practical for clinical workflows. Code is available at https://github.com/yalcintur/WFM.
Yalcin Tur, Mihajlo Stojkovic, Ulas Bagci
Department of Computer Science, Stanford University, Stanford, CA, USA · Machine and Hybrid Intelligence Lab, Feinberg School of Medicine, Northwestern University, Chicago, IL, USA
Diffusion MRI (dMRI) probes water diffusion within the brain, providing unique information about tissue microstructure and white matter connectivity. High-quality dMRI enables advanced and accurate analyses but demands dense acquisitions, resulting in long scan times and hardware constraints that are often infeasible. Further, dMRI requires dedicated computational methods to support diverse analysis tasks (e.g., parcellation, registration, and synthesis) for both neuroscientific and clinical applications. Here, we present dMRI-F, a foundation model that unifies dMRI data enhancement and downstream tasks across diverse acquisitions and populations, while requiring only minimal acquisition input. dMRI-F is pre-trained with a q-space masked reconstruction task: predicting densely sampled multi-shell data from only six volumes. It learns generalizable latent feature representations of dMRI data, which serve as a core foundation for adapting to diverse downstream tasks, including brain parcellation, brain registration, and cross-modal synthesis. Across nine independently acquired datasets, we demonstrate that dMRI-F outperforms state-of-the-art methods with robust generalization across acquisitions, populations, and tasks. Our results demonstrate dMRI-F's versatility as a general-purpose framework that enables diverse dMRI applications from low-cost, sparse acquisitions.
Xi Zhu, Yijie Li, Xiaofan Wang +4
University of Electronic Science and Technology of China.
Diffusion magnetic resonance imaging is a noninvasive technique that is sensitive to the motion of water molecules in human tissue. The image signal is commonly described in q-space; however, a complete coverage of this space is often infeasible due to hardware and acquisition time limitations. In this work, we propose to represent the entire diffusion signal using an implicit neural representation network with wavelet activation functions. The signal space is encoded using Fourier features and spherical harmonics basis functions for a continuous output distribution of the raw diffusion data, enabling both super-resolution and q-space interpolation. Additionally, we introduce a pretraining strategy using a group template, enabling faster per-subject optimisation when combined with nonlinear registrations of subjects into template space. We tested the method with 10 subjects from the Human Connectome Project and 26 subjects from in-house data (multidimensional diffusion MRI from a Parkinson's disease cohort). The results show that the proposed method yields high-fidelity reconstructions with significantly reduced memory requirements compared to raw diffusion-weighted images. The pretraining scheme showed good generalisability properties on data acquired with very different protocols. We further assess the method on downstream metrics relevant to diffusion MRI data from the learned representations.