NiftyPET: High-throughput image reconstruction and analysis¶
NiftyPET is a software platform and a Python namespace package encompassing sub-packages for high-throughput PET image reconstruction, manipulation, processing and analysis with high quantitative accuracy and precision. One of its key applications is brain imaging in dementia with the use of amyloid tracers. See below for the description of the above amyloid PET image reconstructed using NiftyPET, superimposed on the MR T1 weighted image [*].
NiftyPET includes two packages:
nimpa
: https://github.com/NiftyPET/NIMPA (neuro-image manipulation, processing and analysis)nipet
: https://github.com/NiftyPET/NIPET (quantitative PET neuroimaging)
The core routines are written in CUDA C and embedded in Python C extensions to enable user-friendly and high-throughput executions on NVIDIA graphics processing units (GPU). The scientific aspects of this software platform are covered in two open-access publications:
- NiftyPET: a High-throughput Software Platform for High Quantitative Accuracy and Precision PET Imaging and Analysis Neuroinformatics (2018) 16:95. https://doi.org/10.1007/s12021-017-9352-y
- Rapid processing of PET list-mode data for efficient uncertainty estimation and data analysis Physics in Medicine & Biology (2016). https://doi.org/10.1088/0031-9155/61/13/N322
An example application of NiftyPET in the development of novel image reconstruction:
- Faster PET Reconstruction with Non-Smooth Priors by Randomization and Preconditioning Physics in Medicine & Biology (2019). https://doi.org/10.1088/1361-6560/ab3d07
Although, NiftyPET is dedicated to high-throughput image reconstruction and analysis of brain images, it can equally well be used for whole body imaging. Strong emphasis is put on the data, which are acquired using positron emission tomography (PET) and magnetic resonance (MR), especially using the hybrid and simultaneous PET/MR scanners.
This software platform covers the entire processing pipeline, from the raw list-mode (LM) PET data through to the final image statistic of interest (e.g., regional SUV), including LM bootstrapping and multiple independent reconstructions to facilitate voxel-wise estimation of uncertainties.
[*] | The above dynamic transaxial and coronal images show the activity of 18F-florbetapir during the one-hour dynamic acquisition. Note that the signal in the brain white matter dominates over the signal in the grey matter towards the end of the acquisition, which is a typical presentation of a negative amyloid beta (Abeta) scan. |