
Proceedings Paper
Accelerating ring artifact correction for flat-detector CT using the CUDA frameworkFormat | Member Price | Non-Member Price |
---|---|---|
$17.00 | $21.00 |
Paper Abstract
Ring artifacts often appear in flat-detector CT because of imperfect or defect detector elements or calibration.
In high-spatial resolution CT images reducing such artifacts becomes a necessity. In this paper, we used the
post-processing ring correction in polar coordinates (RCP)1 to eliminate the ring artifacts. The median filter is
applied to the uncorrected images in polar coordinates and ring artifacts are extracted from the original images.
The algorithm has a very high computational cost due to the time-expensive median filtering and coordinate
transformation on CPUs. Graphics processing units (GPUs)ca n be seen as parallel co-processors with high
computational power. All steps of the RCP algorithm were implemented with CUDA2(Compute Unified Device
Architecture, NVIDIA). We introduced a new GPU-based branchless vectorized median (BVM)filter. 3, 4 This
algorithm is based on minmax sorting and keeps track of a sorted array from which values are deleted and to
which new values are inserted. For comparison purpose a modified pivot median filter5 on GPUs was presented,
which compares a pivot element to all other values and recursively finds the median element. We evaluated the
performance of the RCP method using 512 slices, each slice consisted of 512×512 pixels. This post-processing
method efficiently reduces ring artifacts in the reconstructed images and improves image quality. Our CUDAbased
RCP is up to 13.6 times faster than the optimized CPU-based (single core)r outine. Comparing our two
GPU-based median filters showed a performance benefit by roughly 60% when switching from Pivot to BVM
code. The main reason is that the BVM algorithm is branchless and makes use of data-level parallelism. The
BVM method is better suited to the model of modern graphics processing. A multi-GPU solution showed that
the performance scaled nearly linearly.
Paper Details
Date Published: 22 March 2010
PDF: 8 pages
Proc. SPIE 7622, Medical Imaging 2010: Physics of Medical Imaging, 76223A (22 March 2010); doi: 10.1117/12.844254
Published in SPIE Proceedings Vol. 7622:
Medical Imaging 2010: Physics of Medical Imaging
Ehsan Samei; Norbert J. Pelc, Editor(s)
PDF: 8 pages
Proc. SPIE 7622, Medical Imaging 2010: Physics of Medical Imaging, 76223A (22 March 2010); doi: 10.1117/12.844254
Show Author Affiliations
W. Chen, Institute of Medical Physics, Friedrich- Alexander-Univ. Erlangen-Nürnberg (Germany)
D. Prell, Institute of Medical Physics, Friedrich- Alexander-Univ. Erlangen-Nürnberg (Germany)
D. Prell, Institute of Medical Physics, Friedrich- Alexander-Univ. Erlangen-Nürnberg (Germany)
Y, Kyriakou, Institute of Medical Physics, Friedrich- Alexander-Univ. Erlangen-Nürnberg (Germany)
W. A. Kalender, Institute of Medical Physics, Friedrich- Alexander-Univ. Erlangen-Nürnberg (Germany)
W. A. Kalender, Institute of Medical Physics, Friedrich- Alexander-Univ. Erlangen-Nürnberg (Germany)
Published in SPIE Proceedings Vol. 7622:
Medical Imaging 2010: Physics of Medical Imaging
Ehsan Samei; Norbert J. Pelc, Editor(s)
© SPIE. Terms of Use
