Accelerating coordinate descent in iterative reconstruction

Scott S. Hsieh, John M. Hoffman, Frederic Noo

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Abstract

Iterative coordinate descent (ICD) is an optimization strategy for iterative reconstruction that is sometimes considered incompatible with parallel compute architectures such as graphics processing units (GPUs). We present a series of modifications that render ICD compatible with GPUs and demonstrate the code on a diagnostic, helical CT dataset. Our reference code is an open-source package, FreeCT ICD, which requires several hours for convergence. Three modifications are used. First, as with our reference code FreeCT ICD, the reconstruction is performed on a rotating coordinate grid, enabling the use of a stored system matrix. Second, every other voxel in the z-is updated direction simultaneously, and the sinogram data is shuffled to coalesce memory access. This increases the parallelism available to the GPU. Third, NS voxels in the xy-plane are updated simultaneously. This introduces possible crosstalk between updated voxels, but because the interaction between non-adjacent voxels is small, small values of NS still converge effectively. We find NS = 16 enables faster reconstruction via greater parallelism, and NS = 256 remains stable but has no additional computational benefit. When tested on a pediatric dataset of size 736x16x14000 reconstructed to a matrix size of 512x512x128 on a single GPU, our implementation of ICD can converge within 10 HU RMS in less than 5 minutes. This suggests that ICD could be competitive with simultaneous update algorithms on modern, parallel compute architectures.

Original languageEnglish (US)
Title of host publicationMedical Imaging 2019
Subtitle of host publicationPhysics of Medical Imaging
EditorsTaly Gilat Schmidt, Guang-Hong Chen, Hilde Bosmans
PublisherSPIE
ISBN (Electronic)9781510625433
DOIs
StatePublished - 2019
EventMedical Imaging 2019: Physics of Medical Imaging - San Diego, United States
Duration: Feb 17 2019Feb 20 2019

Publication series

NameProgress in Biomedical Optics and Imaging - Proceedings of SPIE
Volume10948
ISSN (Print)1605-7422

Conference

ConferenceMedical Imaging 2019: Physics of Medical Imaging
CountryUnited States
CitySan Diego
Period2/17/192/20/19

Keywords

  • GPU programming
  • Iterative reconstruction

ASJC Scopus subject areas

  • Electronic, Optical and Magnetic Materials
  • Biomaterials
  • Atomic and Molecular Physics, and Optics
  • Radiology Nuclear Medicine and imaging

Fingerprint Dive into the research topics of 'Accelerating coordinate descent in iterative reconstruction'. Together they form a unique fingerprint.

Cite this