Skip to main content

Physics-Guided Ultra-Low-Field MRI Enhancement & Simulation

Project description

ULF-Synth

ULF-Synth: Physics-Guided Ultra-Low-Field MRI Enhancement for Pediatric Neuroimaging

arXiv License: MIT Python

Abstract

Ultra-low-field (ULF) MRI enables portable and accessible neuroimaging, but suffers from low signal-to-noise ratio and limited spatial resolution relative to high-field (HF) systems. Acquiring paired ULF–HF data for supervised enhancement is often infeasible, particularly in resource-limited settings. We introduce ULF-Synth, a framework combining: (i) acquisition-based synthesis of realistic ULF images from HF volumes for large-scale paired training, and (ii) a spatial-frequency domain objective that prioritizes recovery of high-frequency anatomical detail. The formulation is architecture-agnostic, consistently improving structural similarity and perceptual fidelity across encoder-decoder, adversarial, and diffusion-based translation models. Trained exclusively on synthetic data, our models generalize to real 64 mT ULF acquisitions, improving multiclass brain segmentation and achieving higher radiologist preference and diagnostic acceptability in a blinded reader study.


Simulation Pipeline

The ULF synthesis module models the key physical phenomena distinguishing ULF from HF acquisitions:

Effect Implementation
1 Signal scaling $(B_{{ULF}}/B_{{HF}})^2$ polarization ratio
2 T2* decay & B0 inhomogeneity Spatially-varying exponential decay from random B0 field maps
3 Thermal noise Gaussian noise scaled to SNR 15–50
4 k-space cropping Reduced resolution (45–55%)
5 k-space undersampling Accelerated acquisition (2×–3×) with center-out sampling
6 B0 off-resonance distortion Phase distortion from random B0 field maps

Qualitative Results

Sample results


Installation

git clone https://github.com/toufiqmusah/ULF-Synth.git
cd ULF-Synth
pip install -e .

This installs ulfsynth and all core dependencies (including PyTorch). The bundled nnUNet translation fork is automatically discovered and installed during setup — no extra steps needed.

Install from PyPI (future)

pip install ulfsynth          # simulation only
pip install ulfsynth[full]    # with enhancement support

CLI

The ulfsynth package provides three commands:

ulfsynth simulate — ULF synthesis from HF volumes

# Single volume
ulfsynth simulate input.nii.gz output.nii.gz

# Folder of NIfTI files
ulfsynth simulate /path/to/hf/scans/ /path/to/ulf/scans/

# Reproducible seed
ulfsynth simulate input.nii.gz output.nii.gz --seed 42

ulfsynth enhance — ULF→HF restoration (requires nnUNet)

# Single volume (CPU)
ulfsynth enhance --device cpu input.nii.gz output.nii.gz

# Folder of NIfTI files (GPU)
ulfsynth enhance /path/to/ulf/scans/ /path/to/enhanced/scans/

# Weights are downloaded from HuggingFace on first use.
# Pre-download: ulfsynth download-weights

ulfsynth download-weights — cache pretrained weights

ulfsynth download-weights

Caches model weights from HuggingFace to ~/.cache/ulfsynth/.


Python API

Simulation

from ulfsynth.simulate import simulate_ulf, simulate_file, simulate_folder, sample_params

# Generate one ULF volume with random parameters
ulf_volume, affine, header, params = simulate_ulf("hf_input.nii.gz")

# With a fixed seed
ulf_volume, affine, header, params = simulate_ulf("hf_input.nii.gz", seed=42)

# Custom parameters
params = sample_params()
params["signal_target"] = 30
ulf_volume, affine, header, params = simulate_ulf("hf_input.nii.gz", params=params)

# Single-file convenience (returns params dict)
params = simulate_file("hf_input.nii.gz", "ulf_output.nii.gz", seed=42)

# Batch folder processing (returns list of params)
results = simulate_folder("hf_scans/", "ulf_scans/", seed=42)

Output preserves the input affine and header metadata.

Enhancement

from ulfsynth.enhance import enhance_file, enhance_folder

# Single file
enhance_file("ulf_input.nii.gz", "enhanced_output.nii.gz", device="cpu")

# Batch folder processing
enhance_folder("ulf_scans/", "enhanced_scans/", device="cuda")

Requires nnUNet (pip install -e src/nn-translation/). Weights are auto-downloaded on first call.


Roadmap

  • Physics-guided ULF synthesis pipeline
  • Pre-trained enhancement weights — ULF→HF restoration models
  • Python package — pip install ulfsynth
  • Docker image — zero-config containerized pipeline

Citation

@misc{musah2026ulfsynth,
  title        = {ULF-Synth: Physics-Guided Ultra-Low-Field MRI Enhancement for Pediatric Neuroimaging},
  author       = {Toufiq Musah and Salvatore Calcagno and Federica Proietto Salanitri and Xiaomeng Li and Maruf Adewole and Marawan Elbatel},
  year         = {2026},
  eprint       = {2605.24625},
  archivePrefix = {arXiv},
  url          = {https://arxiv.org/abs/2605.24625}
}

License

MIT

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

ulfsynth-0.1.10.tar.gz (22.2 MB view details)

Uploaded Source

Built Distribution

If you're not sure about the file name format, learn more about wheel file names.

ulfsynth-0.1.10-py3-none-any.whl (1.4 MB view details)

Uploaded Python 3

File details

Details for the file ulfsynth-0.1.10.tar.gz.

File metadata

  • Download URL: ulfsynth-0.1.10.tar.gz
  • Upload date:
  • Size: 22.2 MB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/6.2.0 CPython/3.12.11

File hashes

Hashes for ulfsynth-0.1.10.tar.gz
Algorithm Hash digest
SHA256 fe826697c0b721416cb93906f8ee070ba55e6589c91b94db6f01c4d696967fcd
MD5 48ee772c4b2f4709b28084d3db076b87
BLAKE2b-256 7b5f19f1f72891762857319c5774f4432a6a173a965848ac9e2d687dd5306964

See more details on using hashes here.

File details

Details for the file ulfsynth-0.1.10-py3-none-any.whl.

File metadata

  • Download URL: ulfsynth-0.1.10-py3-none-any.whl
  • Upload date:
  • Size: 1.4 MB
  • Tags: Python 3
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/6.2.0 CPython/3.12.11

File hashes

Hashes for ulfsynth-0.1.10-py3-none-any.whl
Algorithm Hash digest
SHA256 b9a3e889dd80b855c7dc1ae425b4ef675d39f554db03e4cbc5774011c9773fe3
MD5 2e118711d05299044dbb22aea0bda855
BLAKE2b-256 20333ca9dbd5dc5933d5d6a4ebb2d6fce7fba1c94ca1154456933c3ac4291fb1

See more details on using hashes here.

Supported by

AWS Cloud computing and Security Sponsor Datadog Monitoring Depot Continuous Integration Fastly CDN Google Download Analytics Pingdom Monitoring Sentry Error logging StatusPage Status page