Abstract:Background: The advancement of x-ray CT into the domains of photon counting spectral imaging and dynamic cardiac and perfusion imaging has created many new challenges and opportunities for clinicians and researchers. To address challenges such as dose constraints and scanning times while capitalizing on opportunities such as multi-contrast imaging and low-dose coronary angiography, these multi-channel imaging applications require a new generation of CT reconstruction tools. These new tools should exploit the relationships between imaging channels during reconstruction to set new image quality standards while serving as a platform for direct translation between the preclinical and clinical domains. Purpose: We outline and demonstrate a new Multi-Channel Reconstruction (MCR) Toolkit for GPU-based analytical and iterative reconstruction of preclinical and clinical multi-energy and dynamic x-ray CT data. To promote open science, open-source distribution of the Toolkit will coincide with the release of this publication (GPL v3; gitlab.oit.duke.edu/dpc18/mcr-toolkit-public). Methods: The MCR Toolkit source code is implemented in C/C++ and NVIDIA's CUDA GPU programming interface, with scripting support from MATLAB and Python. The Toolkit implements matched, separable footprint CT reconstruction operators for projection and backprojection in two geometries: planar, cone-beam CT (CBCT) and 3rd generation, cylindrical multi-detector row CT (MDCT). Analytical reconstruction is performed using filtered backprojection (FBP) for circular CBCT, weighted FBP (WFBP) for helical CBCT, and cone-parallel projection rebinning followed by WFBP for MDCT. Arbitrary combinations of energy and temporal channels are iteratively reconstructed under a generalized multi-channel signal model for joint reconstruction. We solve this generalized model algebraically using the split Bregman optimization method and the BiCGSTAB(l) linear solver interchangeably for both CBCT and MDCT data. Rank-sparse kernel regression (RSKR) and patch-based singular value thresholding (pSVT) are used to regularize the energy and time dimensions, respectively. Under a Gaussian noise model, regularization parameters are estimated automatically from the input data, dramatically reducing algorithm complexity for end users. Multi-GPU parallelization of the reconstruction operators is supported to manage reconstruction times. Results: Denoising with RSKR and pSVT and post-reconstruction material decomposition are illustrated with preclinical and clinical cardiac photon-counting (PC)CT data. A digital MOBY mouse phantom with cardiac motion is used to illustrate single energy (SE), multi-energy (ME), time resolved (TR), and combined multi-energy and time-resolved (METR) helical, CBCT reconstruction. A fixed set of projection data is used across all reconstruction cases to demonstrate the Toolkit's robustness to increasing data dimensionality. Identical reconstruction code is applied to in vivo cardiac PCCT data acquired in a mouse model of atherosclerosis (METR). Clinical cardiac CT reconstruction is illustrated using the XCAT phantom and the DukeSim CT simulator, while dual-source, dual-energy CT reconstruction is illustrated for data acquired with a Siemens Flash scanner. Benchmarking results with NVIDIA RTX 8000 GPU hardware demonstrate 61%-99% efficiency in scaling computation from one to four GPUs for these reconstruction problems. Conclusions: The MCR Toolkit provides a robust solution for temporal and spectral x-ray CT reconstruction problems and was built from the ground up to facilitate translation of CT research and development between preclinical and clinical applications.

ORMIR_XCT: A Python package for high resolution peripheral quantitative computed tomography image processing

TIGRE-VarianCBCT for On-Board Cone-Beam Computed Tomography, an Open-Source Toolkit for Imaging, Dosimetry and Clinical Research

Automated quantitative analysis of peri-articular bone microarchitecture in HR-pQCT knee images

PyTomography: A Python Library for Quantitative Medical Image Reconstruction

Quantitative computed tomography (QCT) as a radiology reporting tool by using optical character recognition (OCR) and macro program

TomoPyUI: a user‐friendly tool for rapid tomography alignment and reconstruction

Optimization of Image Quality and Acquisition Time for Lab-Based X-ray Microtomography Using an Iterative Reconstruction Algorithm

Analysis of Iterative Region-of-interest Image Reconstruction for X-Ray Computed Tomography.

MCR toolkit: A GPU-based toolkit for multi-channel reconstruction of preclinical and clinical x-ray CT data

HR-pQCT cross-calibration using standard vs. Laplace-Hamming binarization approach

The clinical application of high-resolution peripheral computed tomography (HR-pQCT) in adults: state of the art and future directions

A Real-time Image Reconstruction System for Particle Treatment Planning Using Proton Computed Tomography (pCT)

Alrecon: computed tomography reconstruction web application based on Solara

Technical Note: MRQy - An open-source tool for quality control of MR imaging data

MRQy: An Open-Source Tool for Quality Control of MR Imaging Data

Towards the adoption of quantitative computed tomography in the management of interstitial lung disease

Statistical Image Reconstruction for High-Throughput Thermal Neutron Computed Tomography

Technical Report: Quality Assessment Tool for Machine Learning with Clinical CT

XCAT-3.0: A Comprehensive Library of Personalized Digital Twins Derived from CT Scans

IRXCT: Iterative Reconstruction and visualization application for X-ray Computed Tomography

Sparse Matrix-Based HPC Tomography