smheart.eu

Code and Data

We’ll make most of our codes (mostly MATLAB/PYTHON) available here. Please be aware that this is academic code, and has mostly not been optimized for efficiency or coding elegance. All codes and data are provided as is, and all warranty is excluded.

You are welcome to provide feedback and report bugs to aurelien.bustin@ihu-liryc.fr.


Primal-Dual Beltrami Regularization

The geometrically derived Beltrami framework, introduced by Sochen, Kimmel and Malladi (1998) offers an ideal compromise between feature preservation and avoidance of staircasing artifacts. Until now, one of the main limiting factors of the Beltrami regularizers was the lack of really efficient optimization schemes. Here, we start from one of the most efficient TV-optimization methods, primal-dual projected gradients, and apply it to the Beltrami functional. Doing so, we achieve better performance than ROF denoising for the basic grey-scale denoising problem, then extend the method to more involved problems such as inpainting, deconvolution, and the color case, all in a straightforward fashion. With the proposed primal-dual projected gradients optimization algorithm, the benefits of the geometric Beltrami regularizer become available at no extra computational cost, compared to state-of-the-art TV/ROF regularizers.

MATLAB Example

The MATLAB scripts/functions are included here to provide a way to illustrate the Beltrami algorithm.

Downloads – MATLAB scripts/functions demonstrating Beltrami denoising

When using the code, please refer to the following publication for documentation and as reference to be cited:

Reference

  1. Zosso D, Bustin A. A Primal-Dual Projected Gradient Algorithm for Efficient Beltrami Regularization, UCLA CAM Report 14-52, 2014. [PDF]

Myocardial T1 Mapping Denoising using Vectorized Beltrami Regularization

MATLAB/MEX Example

Code coming soon …

Reference

  1. Bustin A, Ferry P, Codreanu A, Beaumont M, Liu S, Burschka D, Felblinger J, Brau A, Menini A, Odille F. Impact of Denoising on Precision and Accuracy of Saturation-Recovery-Based Myocardial T1 mapping, Journal of Magnetic Resonance Imaging, 46(5):1377-1388, DOI: 10.1002/jmri.25684 [PDF]

2D/3D single-contrast patch-based denoising

MATLAB/MEX Example

Code coming soon …

References

  1. Bustin A, Ginami G, Cruz G, Correia T, Ismail T, Rashid I, Neji R, Botnar R, Prieto C. Five-minute whole-heart coronary MRA with sub-millimetre isotropic resolution, 100% respiratory scan efficiency, and 3D-PROST reconstruction, Magnetic Resonance in Medicine, 81(1):102-115, 2019, DOI: 10.1002/mrm.27354 [PDF]
  2. Bustin A, Voilliot D, Menini A, Felblinger J, de Chillou C, Burschka D, Bonnemains L, Odille F. Isotropic reconstruction of MR images using 3D patch-based self-similarity learning, IEEE Transactions on Medical Imaging, 37(8):1932-1942, 2018, 10.1109/TMI.2018.2807451 [PDF]

2D/3D multi-contrast patch-based denoising

This technique jointly reconstructs multi-contrast MR images by exploiting the highly redundant information, on a local and non-local scale, and the strong correlation shared between the multiple contrast images. This is achieved by enforcing multi-dimensional low-rank in the under sampled images.

MATLAB/MEX Example

Code coming soon …

Reference

  1. Bustin A, Cruz G, Jaubert O, Lopez K, Botnar R, Prieto C. High-dimensionality under sampled patch-based reconstruction (HD-PROST) for accelerated multi-contrast MRI, Magnetic Resonance in Medicine, 2019, 1-15. DOI: 10.1002/mrm.27694 [PDF]