2,658 research outputs found
Parallel CT image reconstruction based on GPUs
[EN] In X-ray computed tomography (CT) iterative methods are more suitable for the reconstruction of images with high contrast and precision in noisy conditions from a small number of projections. However, in practice, these methods are not widely used due to the high computational cost of their implementation. Nowadays technology provides the possibility to reduce effectively this drawback. It is the goal of this work to develop a fast GPU-based algorithm to reconstruct high quality images from under sampled and noisy projection data.Research supported by ANITRAN Project PROMETEO/2010/039.Flores, LA.; Vidal Gimeno, VE.; Mayo Nogueira, P.; RĂłdenas Escribá, FDA.; VerdĂş MartĂn, GJ. (2014). Parallel CT image reconstruction based on GPUs. Radiation Physics and Chemistry. 95(1):247-250. https://doi.org/10.1016 / j.radphyschem.2013.03.011S24725095
Accelerating iterative CT reconstruction algorithms using Tensor Cores
Tensor Cores are specialized hardware units added to recent NVIDIA GPUs to speed up matrix multiplication-related tasks, such as convolutions and densely connected layers in neural networks. Due to their specific hardware implementation and programming model, Tensor Cores cannot be straightforwardly applied to other applications outside machine learning. In this paper, we demonstrate the feasibility of using NVIDIA Tensor Cores for the acceleration of a non-machine learning application: iterative Computed Tomography (CT) reconstruction. For large CT images and real-time CT scanning, the reconstruction time for many existing iterative reconstruction methods is relatively high, ranging from seconds to minutes, depending on the size of the image. Therefore, CT reconstruction is an application area that could potentially benefit from Tensor Core hardware acceleration. We first studied the reconstruction algorithm's performance as a function of the hardware related parameters and proposed an approach to accelerate reconstruction on Tensor Cores. The results show that the proposed method provides about 5 x increase in speed and energy saving using the NVIDIA RTX 2080 Ti GPU for the parallel projection of 32 images of size 512 x 512. The relative reconstruction error due to the mixed-precision computations was almost equal to the error of single-precision (32-bit) floating- point computations. We then presented an approach for real-time and memory-limited applications by exploiting the symmetry of the system (i.e., the acquisition geometry). As the proposed approach is based on the conjugate gradient method, it can be generalized to extend its application to many research and industrial fields
Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT
Modern CT image reconstruction algorithms rely on projection and back-projection operations to refine an image estimate in iterative image reconstruction. A widely-used state-of-the-art technique is distance-driven projection and back-projection. While the distance-driven technique yields superior image quality in iterative algorithms, it is a computationally demanding process. This has a detrimental effect on the relevance of the algorithms in clinical settings. A few methods have been proposed for enhancing the distance-driven technique in order to take advantage of modern computer hardware. This study explores a two-dimensional extension of the branchless method, which is a technique that does not compromise image quality. The extension of the branchless method is named “pre-projection integration” because it gets a performance boost by integrating the data before the projection and back-projection operations. It was written with Nvidia’s CUDA framework and carefully designed for massively parallel graphics processing units (GPUs). The performance and the image quality of the pre-projection integration method were analyzed. Both projection and back-projection are significantly faster with pre-projection integration. The image quality was analyzed using cone beam CT image reconstruction algorithms within Jeffrey Fessler’s Image Reconstruction Toolbox. Images produced from regularized, iterative image reconstruction algorithms using the pre-projection integration method show no significant artifacts
Four-dimensional tomographic reconstruction by time domain decomposition
Since the beginnings of tomography, the requirement that the sample does not
change during the acquisition of one tomographic rotation is unchanged. We
derived and successfully implemented a tomographic reconstruction method which
relaxes this decades-old requirement of static samples. In the presented
method, dynamic tomographic data sets are decomposed in the temporal domain
using basis functions and deploying an L1 regularization technique where the
penalty factor is taken for spatial and temporal derivatives. We implemented
the iterative algorithm for solving the regularization problem on modern GPU
systems to demonstrate its practical use
Investigation of iterative image reconstruction in three-dimensional optoacoustic tomography
Iterative image reconstruction algorithms for optoacoustic tomography (OAT),
also known as photoacoustic tomography, have the ability to improve image
quality over analytic algorithms due to their ability to incorporate accurate
models of the imaging physics, instrument response, and measurement noise.
However, to date, there have been few reported attempts to employ advanced
iterative image reconstruction algorithms for improving image quality in
three-dimensional (3D) OAT. In this work, we implement and investigate two
iterative image reconstruction methods for use with a 3D OAT small animal
imager: namely, a penalized least-squares (PLS) method employing a quadratic
smoothness penalty and a PLS method employing a total variation norm penalty.
The reconstruction algorithms employ accurate models of the ultrasonic
transducer impulse responses. Experimental data sets are employed to compare
the performances of the iterative reconstruction algorithms to that of a 3D
filtered backprojection (FBP) algorithm. By use of quantitative measures of
image quality, we demonstrate that the iterative reconstruction algorithms can
mitigate image artifacts and preserve spatial resolution more effectively than
FBP algorithms. These features suggest that the use of advanced image
reconstruction algorithms can improve the effectiveness of 3D OAT while
reducing the amount of data required for biomedical applications
- …