skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: MemXCT: Design, Optimization, Scaling, and Reproducibility of X-Ray Tomography Imaging

Journal Article · · IEEE Transactions on Parallel and Distributed Systems

Here, this work extends our previous research entitled "MemXCT: Memory-centric X-ray CT Reconstruction with Massive Parallelization" that was originally published at SC19 conference (Hidayetoglu et al., 2019) with reproducibility of the computational imaging performance. X-ray computed tomography (XCT) is regularly used at synchrotron light sources to study the internal morphology of materials at high resolution. However, experimental constraints, such as radiation sensitivity, can result in noisy or undersampled measurements. Further, depending on the resolution, sample size and data acquisition rates, the resulting noisy dataset can be in the order of terabytes. Advanced iterative reconstruction techniques can produce high-quality images from noisy measurements, but their computational requirements have made their use an exception rather than the rule. We propose a novel memory-centric approach that avoids redundant computations at the expense of additional memory complexity. We develop a memory-centric iterative reconstruction system, MemXCT, that uses an optimized SpMV implementation with two-level pseudo-Hilbert ordering and multi-stage input buffering. We evaluate MemXCT on various supercomputer architectures involving KNL and GPU. MemXCT can reconstruct a large (11Kx11K) mouse brain tomogram in 10 seconds using 4096 KNL nodes (256K cores). The results presented in our original article at the SC19 were based on large-scale supercomputing resources. The MemXCT application was selected for the Student Cluster Competition (SCC) Reproducibility Challenge and evaluated on a variety of cloud computing resources by universities around the world in the SC20 conference. We summarize the results of the top-ranked SCC Reproducibility Challenge teams and identify the most pertinent measures for ensuring the reproducibility of our experiments in this article.

Research Organization:
Argonne National Laboratory (ANL), Argonne, IL (United States)
Sponsoring Organization:
National Science Foundation (NSF); Semiconductor Research Corporation (SRC); US Department of Defense (DoD); Defense Advanced Research Projects Agency (DARPA); USDOE Office of Science (SC)
Grant/Contract Number:
AC02-06CH11357; OCI-0725070; ACI-1238993; NA0002374
OSTI ID:
1871328
Journal Information:
IEEE Transactions on Parallel and Distributed Systems, Vol. 33, Issue 9; ISSN 1045-9219
Publisher:
IEEECopyright Statement
Country of Publication:
United States
Language:
English

References (44)

Model-based Iterative CT Image Reconstruction on GPUs
  • Sabne, Amit; Wang, Xiao; Kisner, Sherman J.
  • PPoPP '17: 22nd ACM SIGPLAN Symposium on Principles and Practice of Parallel Programming, Proceedings of the 22nd ACM SIGPLAN Symposium on Principles and Practice of Parallel Programming https://doi.org/10.1145/3018743.3018765
conference January 2017
Hybrid MPI-OpenMP Programming for Parallel OSEM PET Reconstruction journal October 2006
Fast tomographic reconstruction on multicore computers journal December 2010
Bestimmung der Absorption des rothen Lichts in farbigen Flüssigkeiten journal January 1852
Pushing the limits for medical image reconstruction on recent standard multicore processors journal June 2012
Automatic CPU-GPU communication management and optimization journal June 2011
Parallel Transposition of Sparse Data Structures conference January 2016
High performance model based image reconstruction journal November 2016
Plug-and-Play priors for model based reconstruction conference December 2013
A comparative study of X-ray tomographic microscopy on shales at different synchrotron facilities: ALS, APS and SLS journal November 2012
Model-based iterative reconstruction for synchrotron X-ray tomography conference May 2014
Special Issue on SC16 Student Cluster Competition Reproducibility Initiative journal December 2017
Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods book January 1994
A fast forward projection using multithreads for multirays on GPUs in medical image reconstruction: Fast forward projection on GPUs journal June 2011
Hyperspectral image reconstruction for x-ray fluorescence tomography journal January 2015
cuMBIR conference June 2018
Analysis of the clustering properties of the Hilbert space-filling curve journal January 2001
A Cache‐Aware Algorithm for PDEs on Hierarchical Data Structures Based on Space‐Filling Curves journal January 2006
Using Peano Curves for Bilevel Display of Continuous-Tone Images journal May 1982
Exploring the Performance Benefit of Hybrid Memory System on HPC Environments conference May 2017
Accelerating advanced MRI reconstructions on GPUs journal October 2008
Algorithm-Enabled Low-Dose Micro-CT Imaging journal March 2011
Fast calculation of the exact radiological path for a three-dimensional CT array: Technical Reports: 3D CT array path calculation journal March 1985
TomoPy: a framework for the analysis of synchrotron tomographic data journal August 2014
A high-throughput x-ray microtomography system at the Advanced Photon Source journal April 2001
Sur une courbe, qui remplit toute une aire plane journal March 1890
Maximum a posteriori estimation of crystallographic phases in X-ray diffraction tomography
  • Gürsoy, Doĝa; Biçer, Tekin; Almer, Jonathan D.
  • Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, Vol. 373, Issue 2043 https://doi.org/10.1098/rsta.2014.0392
journal June 2015
Automatic Generation of Warp-Level Primitives and Atomic Instructions for Fast and Portable Parallel Reduction on GPUs conference February 2019
Cache oblivious matrix multiplication using an element ordering based on a Peano curve journal September 2006
TomoBank: a tomographic data repository for computational x-ray science journal February 2018
Fiji: an open-source platform for biological-image analysis journal June 2012
Time-resolved X-ray Tomography of Gasoline Direct Injection Sprays journal September 2015
Trace: a high-throughput tomographic reconstruction engine for large-scale datasets journal January 2017
A data-parallel algorithm for iterative tomographic image reconstruction conference January 1999
Iterative reconstruction methods in X-ray CT journal April 2012
Quantifying Mesoscale Neuroanatomy Using X-Ray Microtomography journal September 2017
MemXCT: memory-centric X-ray CT reconstruction with massive parallelization
  • Hidayetoğlu, Mert; Biçer, Tekin; de Gonzalo, Simon Garcia
  • SC '19: The International Conference for High Performance Computing, Networking, Storage, and Analysis, Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis https://doi.org/10.1145/3295500.3356220
conference November 2019
On partitioning dynamic adaptive grid hierarchies conference January 1996
CUDA optimization strategies for compute- and memory-bound neuroimaging algorithms journal June 2012
A Study of Energy and Locality Effects Using Space-Filling Curves conference May 2014
Massively parallel 3D image reconstruction
  • Wang, Xiao; Sabne, Amit; Sakdhnagool, Putt
  • SC '17: The International Conference for High Performance Computing, Networking, Storage and Analysis, Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis https://doi.org/10.1145/3126908.3126911
conference November 2017
Querying multi-dimensional data indexed using the Hilbert space-filling curve journal March 2001
The ASTRA Toolbox: A platform for advanced algorithm development in electron tomography journal October 2015
Recent Advances in CT Image Reconstruction journal January 2013

Similar Records

A sparse matrix–vector multiplication based algorithm for accurate density matrix computations on systems of millions of atoms
Journal Article · Thu Feb 15 00:00:00 EST 2018 · Computer Physics Communications · OSTI ID:1871328

Resolution Versus Error for Computational Electron Microscopy
Journal Article · Sat Jul 01 00:00:00 EDT 2017 · Microscopy and Microanalysis · OSTI ID:1871328

Planetary normal mode computation: Parallel algorithms, performance and reproducibility
Journal Article · Mon Jan 11 00:00:00 EST 2021 · IEEE Transactions on Parallel and Distributed Systems · OSTI ID:1871328