MRI Non-Uniform Fast Fourier Transform (MRI-NUFFT)
2026-06-18 Pierre-Antoine Comby, Guillaume Daval‐Frérot, Asma Tanabene, et al. Zenodo (CERN European Organization for Nuclear Research)

MRI-NUFFT is an open-source Python library that provides state-of-the-art non-Cartesian MRI tools: trajectories, data loading, fast and memory-efficient operators to be used on laptops, clusters, and MRI consoles.

DD-INR: Dynamics-Driven Implicit Neural Representation for Accelerated Whole-Brain Functional MRI Reconstruction
2026-06-09 Qiaoxin Li, Caini Pan, Pierre-Antoine Comby, et al. Lecture Notes in Computer Science

Accelerated acquisition of fMRI enables enhanced detection of neurovascular (BOLD) activity in the brain, but image reconstruction becomes challenging with high k-space undersampling: Task-evoked BOLD signals are small in magnitude, which traditional anatomical MRI reconstruction methods fail to recover, as they favor…

Weakly Convex Ridge Regularization for 3D Non-Cartesian MRI Reconstruction
2026-03-28 German Shâma Wache, Chaithya Giliyar Radhakrishna, Asma Tanabene, Sebastian Neumayer arXiv (Cornell University)

While highly accelerated non-Cartesian acquisition protocols significantly reduce scan time, they often entail long reconstruction delays. Deep learning based reconstruction methods can alleviate this, but often lack stability and robustness to distribution shifts.

On the Impact of Artifacts Induced by Mismatches Between Auto‐Calibration Signal and Accelerated 3D GRE Data at 11.7T
2025-10-18 Joseph Obriot, Franck Mauconduit, Vincent Gras, et al. Magnetic Resonance in Medicine

PURPOSE: field inhomogeneity artifacts arising from remote locations in the FOV and encountered in accelerated 3D gradient-recalled echo (GRE) sequences at ultra-high field, and at providing mitigation strategies. METHODS: Measurements were conducted at 11.7T using a head-shaped phantom and an accelerated 3D GRE…

Bringing GRAPPA to non-Cartesian MRI through SPARKLING: An application to MP-RAGE anatomical MRI
2025-09-16 Chaithya Giliyar Radhakrishna, Alexandre Vignaud, Maxime Bertrait, et al. ISMRM & ISMRT 2025-Annual Meeting & Exhibition

Motivation: Well established clinically used parallel imaging techniques like GRAPPA and CAIPIRINHA are mostly limited to Cartesian sampling, and their extension to non-Cartesian sampling is not direct. Goal(s): We propose to extend the SPARKLING framework to ensure GRAPPA accelerated Cartesian sampling in k-space…

gGRAPPA: A flexible, GPU-accelerated python package for fast and efficient generalized GRAPPA reconstruction
2025-09-16 Maxime Bertrait, Chaithya Giliyar Radhakrishna, Philippe Ciuciu ISMRM & ISMRT 2025-Annual Meeting & Exhibition

Motivation: Existing open-source MR reconstruction tools often fail to efficiently utilize GPU resources and lack support for generalized GRAPPA implementations. Many tools are limited to 2D or 3D reconstruction, and few incorporate advanced techniques such as 2D-CAIPIRINHA, which enhances imaging capabilities.

MRI-NUFFT: Doing non-Cartesian MRI has never been easier
2025-04-15 Pierre-Antoine Comby, Guillaume Daval‐Frérot, Caini Pan, et al. Journal of Open Source Software

MRI-NUFFT is a python package that provides a universal interface to various Non-Uniform Fast Fourier Transform libraries running on CPU or GPU (gpuNUFFT, finufft, cufinufft, pyNFFT), adding compatibility with standard array library (NumPy, CuPy, PyTorch, TensorFlow, etc.) On top of these libraries it extends the…

MRI-NUFFT: An open source Python package to make non-Cartesian MR Imaging easier
2024-11-26 Pierre-Antoine Comby, Guillaume Daval-Frerot, Chaithya Giliyar Radhakrishna, et al. 2025 ISMRM & ISMRT Annual Meeting & Exhibition

Motivation: Non-Cartesian imaging remains complicated to use for MRI due to the high computational cost of the Non-Uniform Fourier Transform for image reconstruction. Goal(s): To provide a uniform interface for reconstructing magnetic resonance images from non-Cartesian k-space data and a collection of non-Cartesian…

Benchmarking common and advanced non-Cartesian trajectories with high acceleration and static off-resonance effects
2024-08-14 Guillaume Daval-Frérot, Chaithya Giliyar Radhakrishna, Fernanda Ponce, et al. ISMRM & ISMRT Annual Meeting & Exhibition

Non-Cartesian sampling patterns allow for highly accelerated MRI exams at the cost of exacerbated and more complex artifacts, each impacting image quality in unique ways. Patient-induced B0 field inhomogeneities can notably cause distortions and blurring, but some trajectories are by design more robust to them.

Non-cartesian 3D-SPARKLING vs cartesian 3D-EPI encoding schemes for functional magnetic resonance Imaging at 7 Tesla
2024-05-13 Zaineb Amor, Philippe Ciuciu, Chaithya Giliyar Radhakrishna, et al. Plos one

The quest for higher spatial and/or temporal resolution in functional MRI (fMRI) while preserving a sufficient temporal signal-to-noise ratio (tSNR) has generated a tremendous amount of methodological contributions in the last decade ranging from Cartesian vs. non-Cartesian readouts, 2D vs.

B0 field distortions monitoring and correction for 3D non-Cartesian fMRI acquisitions using a field camera: Application to 3D-SPARKLING at 7T
2023-08-03 Zaineb Amor, Chaithya Giliyar Radhakrishna, Caroline Le Ster, et al. Joint Annual Meeting ISMRM-ESMRMB & ISMRT 31st Annual Meeting

The benefit of retrospective correction of static and dynamic B0 field imperfections was studied on non-Cartesian compressed sensing-based acquisition patterns by taking as an example 3D-SPARKLING encoding scheme. Long-TR dynamic acquisitions of T2*-weighted MRI (fMRI-like) volumes were acquired at 7T, retrospectively…

Prospects of non-Cartesian 3D-SPARKLING encoding for functional MRI: A preliminary case study for retinotopic mapping
2023-08-03 Zaineb Amor, Chaithya Giliyar Radhakrishna, Guillaume Daval-Frérot, et al. Joint Annual Meeting ISMRM-ESMRMB & ISMRT 31st Annual Meeting

For the first time, non-Cartesian 3D-SPARKLING (Spreading Projection Algorithm for Rapid K-space sampLING) encoding scheme was used for task-based fMRI (retinotopic mapping) at 7T. Additionally, this new acquisition technique was compared with 3D-EPI, considered as one of the reference acquisition schemes in…

Impact of ΔB0 field imperfections correction on BOLD sensitivity in 3D-SPARKLING fMRI data
2023-06-03 Zaineb Amor, Caroline Le Ster, Chaithya Giliyar Radhakrishna, et al. Magnetic Resonance in Medicine

Purpose Static and dynamic field imperfections are detrimental to functional MRI (fMRI) applications, especially at ultra‐high magnetic fields (UHF). In this work, a field camera is used to assess the benefits of retrospectively correcting field perturbations on Blood Oxygen Level Dependent (BOLD) sensitivity in…

Improving spreading projection algorithm for rapid k-space sampling trajectories through minimized off-resonance effects and gridding of low frequencies
2023-05-22 Chaithya Giliyar Radhakrishna, Guillaume Daval-Frérot, Aurélien Massire, et al. Magnetic Resonance in Medicine

Purpose Non‐Cartesian MRI with long arbitrary readout directions are susceptible to off‐resonance artifacts due to patient induced inhomogeneities. This results in degraded image quality with strong signal losses and blurring.

Designing and learning non-Cartesian k-space sampling trajectories for accelerated 3D MRI
2023-04-25 Chaithya Giliyar Radhakrishna theses.fr (ABES)

Conception et apprentissage de trajectoires d'échantillonnage non cartésiennes dans l'espace k pour l'IRM 3D accélérée L'imagerie par résonance magnétique (IRM) est une technique d'imagerie médicale non invasive très répandue qui permet de comprendre les structures des tissus du corps humain. Cependant, sa portée et…

Jointly learning non-cartesian k-space trajectories and reconstruction networks for 2D and 3D MR imaging through projection
2023-01-24 Chaithya Giliyar Radhakrishna, Philippe Ciuciu Bioengineering

Compressed sensing in magnetic resonance imaging essentially involves the optimization of (1) the sampling pattern in k-space under MR hardware constraints and (2) image reconstruction from undersampled k-space data. Recently, deep learning methods have allowed the community to address both problems simultaneously,…

Benchmarking learned non-Cartesian k-space trajectories and reconstruction networks
2022-01-27 Chaithya Giliyar Radhakrishna, Philippe Ciuciu CoRR

We benchmark the current existing methods to jointly learn non-Cartesian k-space trajectory and reconstruction: PILOT, BJORK, and compare them with those obtained from the recently developed generalized hybrid learning (HybLearn) framework. We present the advantages of using projected gradient descent to enforce MR…

NC-PDNet: A density-compensated unrolled network for 2D and 3D non-Cartesian MRI reconstruction
2022-01-18 Zaccharie Ramzi, Chaithya Giliyar Radhakrishna, Jean-Luc Starck, Philippe Ciuciu IEEE Transactions on Medical Imaging

Deep Learning has become a very promising avenue for magnetic resonance image (MRI) reconstruction. In this work, we explore the potential of unrolled networks for non-Cartesian acquisition settings.

Hybrid learning of non-cartesian k-space trajectory and MR image reconstruction networks
2021-10-25 Chaithya Giliyar Radhakrishna, Zaccharie Ramzi, Philippe Ciuciu 2022 IEEE 19th International Symposium on Biomedical Imaging (ISBI)

Compressed sensing (CS) in Magnetic resonance Imaging (MRI) essentially involves the optimization of 1) the sampling pattern in k-space under MR hardware constraints and 2) image reconstruction from the undersampled k-space data. Recently, deep learning methods have allowed the community to address both problems…

Learning the sampling density in 2D SPARKLING MRI acquisition for optimized image reconstruction
2021-08-23 Chaithya Giliyar Radhakrishna, Zaccharie Ramzi, Philippe Ciuciu 2021 29th European Signal Processing Conference (EUSIPCO)

The SPARKLING algorithm was originally developed for accelerated 2D magnetic resonance imaging (MRI) in the compressed sensing (CS) context. It yields non-Cartesian sampling trajectories that jointly fulfill a target sampling density while each individual trajectory complies with MR hardware constraints.

Optimizing full 3D SPARKLING trajectories for high-resolution T2*-weighted Magnetic Resonance Imaging
2021-08-06 Chaithya Giliyar Radhakrishna, Pierre Weiss, Guillaume Daval-Frérot, et al. IEEE Transactions on Medical Imaging

The Spreading Projection Algorithm for Rapid K-space samplING, or SPARKLING, is an optimization-driven method that has been recently introduced for accelerated 2D T2*-w MRI using compressed sensing. It has then been extended to address 3D imaging using either stacks of 2D sampling patterns or a local 3D strategy that…

Calibration-less multi-coil compressed sensing magnetic resonance image reconstruction based on OSCAR regularization
2021-03-19 Loubna El Gueddari, Chaithya Giliyar Radhakrishna, Emilie Chouzenoux, Philippe Ciuciu Journal of Imaging

Over the last decade, the combination of compressed sensing (CS) with acquisition over multiple receiver coils in magnetic resonance imaging (MRI) has allowed the emergence of faster scans while maintaining a good signal-to-noise ratio (SNR). Self-calibrating techniques, such as ESPiRIT, have become the standard…

PySAP: Python Sparse Data Analysis Package for multidisciplinary image processing
2020-06-25 S. Farrens, Antoine Grigis, Loubna El Gueddari, et al. Astronomy and Computing

We present the open-source image processing software package PySAP (Python Sparse data Analysis Package) developed for the COmpressed Sensing for Magnetic resonance Imaging and Cosmology (COSMIC) project. This package provides a set of flexible tools that can be applied to a variety of compressed sensing and image…

A multi-center study on human brain glutathione conformation using magnetic resonance spectroscopy
2018-09-04 Deepika Shukla, Pravat K Mandal, Lars Ersland, et al. Journal of Alzheimer’s Disease

Molecular dynamics simulation and in vitro nuclear magnetic resonance (NMR) studies on glutathione (GSH) indicated existence of closed and extended conformations. The present work in a multi-center research setting reports in-depth analysis of GSH conformers in vivo using a common magnetic resonance spectroscopy (MRS)…