Abstract:Image reconstruction in soft-field tomography is based on an inverse problem formulation, where a forward model is fitted to the data. In medical applications, where the anatomy presents complex shapes, it is common to use finite element models (FEMs) to represent the volume of interest and solve a partial differential equation that models the physics of the system. Over the last decade, there has been a shifting interest from 2D modeling to 3D modeling, as the underlying physics of most problems are 3D. Although the increased computational power of modern computers allows working with much larger FEM models, the computational time required to reconstruct 3D images on a fine 3D FEM model can be significant, on the order of hours. For example, in electrical impedance tomography (EIT) applications using a dense 3D FEM mesh with half a million elements, a single reconstruction iteration takes approximately 15-20 min with optimized routines running on a modern multi-core PC. It is desirable to accelerate image reconstruction to enable researchers to more easily and rapidly explore data and reconstruction parameters. Furthermore, providing high-speed reconstructions is essential for some promising clinical application of EIT. For 3D problems, 70% of the computing time is spent building the Jacobian matrix, and 25% of the time in forward solving. In this work, we focus on accelerating the Jacobian computation by using single and multiple GPUs. First, we discuss an optimized implementation on a modern multi-core PC architecture and show how computing time is bounded by the CPU-to-memory bandwidth; this factor limits the rate at which data can be fetched by the CPU. Gains associated with the use of multiple CPU cores are minimal, since data operands cannot be fetched fast enough to saturate the processing power of even a single CPU core. GPUs have much faster memory bandwidths compared to CPUs and better parallelism. We are able to obtain acceleration factors of 20 times on a single NVIDIA S1070 GPU, and of 50 times on four GPUs, bringing the Jacobian computing time for a fine 3D mesh from 12 min to 14 s. We regard this as an important step toward gaining interactive reconstruction times in 3D imaging, particularly when coupled in the future with acceleration of the forward problem. While we demonstrate results for EIT, these results apply to any soft-field imaging modality where the Jacobian matrix is computed with the adjoint method.

Gpu Acceleration of Acoustic Radiation Force Impulse Imaging Using Scanning-Mode

Scanning-mode 2D Acoustic Radiation Force Impulse (S2d-Arfi) Imaging Based on GPU Acceleration

Fast and High-Resolution Acoustic Beamforming: A Convolution Accelerated Deconvolution Implementation

Sound Speed Image Reconstruction of Ultrasound Computed Tomography Based on GPU Acceleration

FPGA Acceleration of Image Reconstruction for Real-Time Photoacoustic Tomography

GPU-Based Parallel Computing Methods for Medical Photoacoustic Image Reconstruction

Real-time Photoacoustic and Ultrasound Dual-Modality Imaging System Facilitated with Graphics Processing Unit and Code Parallel Optimization

Parallel Computing For Quantitative Blood Flow Imaging In Photoacoustic Microscopy

Toward Parallel Optimal Computation of Ultrasound Computed Tomography Using GPU

Parallel Accelerated SAR Image Simulation Based on High Frequency Approximation Techniques

Parallel processing of SAFT ultrasonic imaging technique based on GPU

Parallel Implementation of Real-Time ISAR Imaging Based on Graphics Processing Units

Accelerating Image Reconstruction in Three-Dimensional Optoacoustic Tomography on Graphics Processing Units

Comprehensive framework of GPU-accelerated image reconstruction for photoacoustic computed tomography

GPU Accelerated High Intensity Ultrasound Acoustical Power Computation

Motion compensation and autofocus issues for airborne high-resolution SAR based on GPU architecture

Swept-source endoscopic optical coherence tomography real-time imaging system based on GPU acceleration for axial megahertz high-speed scanning

AI Enabled High Frame Rate Portable Ultrasound Imaging Pipeline: Prototype Implementation with GPU Acceleration

Precise Tracking of Impulsive Acoustic Radiation Force Induced Small Displacements for Shear Wave Speed Estimation

A Fast Inverse Synthetic Aperture Radar Imaging Scheme Combining GPU-Accelerated Shooting and Bouncing Ray and Back Projection Algorithm under Wide Bandwidths and Angles

Multi-GPU Jacobian accelerated computing for soft-field tomography