Accelerating simultaneous algebraic reconstruction technique with motion compensation using CUDA-enabled GPU

Computed Tomography as a Characterization Tool for Engineered Scaffolds with Biomedical Applications

The ever-growing field of materials with applications in the biomedical field holds great promise regarding the design and fabrication of devices with specific characteristics, especially scaffolds with personalized geometry and architecture. The continuous technological development pushes the limits of innovation in obtaining adequate scaffolds and establishing their characteristics and performance. To this end, computed tomography (CT) proved to be a reliable, nondestructive, high-performance machine, enabling visualization and structure analysis at submicronic resolutions. CT allows both qualitative and quantitative data of the 3D model, offering an overall image of its specific architectural features and reliable numerical data for rigorous analyses. The precise engineering of scaffolds consists in the fabrication of objects with well-defined morphometric parameters (e.g., shape, porosity, wall thickness) and in their performance validation through thorough control over their behavior (in situ visualization, degradation, new tissue formation, wear, etc.). This review is focused on the use of CT in biomaterial science with the aim of qualitatively and quantitatively assessing the scaffolds’ features and monitoring their behavior following in vivo or in vitro experiments. Furthermore, the paper presents the benefits and limitations regarding the employment of this technique when engineering materials with applications in the biomedical field.

Direct Measure of Electrode Spatial Heterogeneity: Influence of Processing Conditions on Anode Architecture and Performance

In this work, the spatial (in)homogeneity of aqueous processed silicon electrodes using standard poly(acrylic acid)-based binders and slurry preparation conditions is demonstrated. X-ray nanotomography shows segregation of materials into submicron-thick layers depending on the mixing method and starting binder molecular weights. Using a dispersant, or in situ production of dispersant from the cleavage of the binder into smaller molecular weight species, increases the resulting lateral homogeneity while drastically decreasing the vertical homogeneity as a result of sedimentation and separation due to gravitational forces. This data explains some of the variability in the literature with respect to silicon electrode performance and demonstrates two potential ways to improve slurry-based electrode fabrications.

Reconstruction of a high-quality volumetric image and a respiratory motion model from patient CBCT projections

PURPOSE

To develop and evaluate a method of reconstructing a patient- and treatment day- specific volumetric image and motion model from free-breathing cone-beam projections and respiratory surrogate measurements. This Motion-Compensated Simultaneous Algebraic Reconstruction Technique (MC-SART) generates and uses a motion model derived directly from the cone-beam projections, without requiring prior motion measurements from 4DCT, and can compensate for both inter- and intrabin deformations. The motion model can be used to generate images at arbitrary breathing points, which can be used for estimating volumetric images during treatment delivery.

METHODS

The MC-SART was formulated using simultaneous image reconstruction and motion model estimation. For image reconstruction, projections were first binned according to external surrogate measurements. Projections in each bin were used to reconstruct a set of volumetric images using MC-SART. The motion model was estimated based on deformable image registration between the reconstructed bins, and least squares fitting to model parameters. The model was used to compensate for motion in both projection and backprojection operations in the subsequent image reconstruction iterations. These updated images were then used to update the motion model, and the two steps were alternated between. The final output is a volumetric reference image and a motion model that can be used to generate images at any other time point from surrogate measurements.

RESULTS

A retrospective patient dataset consisting of eight lung cancer patients was used to evaluate the method. The absolute intensity differences in the lung regions compared to ground truth were 50.8 ± 43.9 HU in peak exhale phases (reference) and 80.8 ± 74.0 in peak inhale phases (generated). The 50th percentile of voxel registration error of all voxels in the lung regions with >5 mm amplitude was 1.3 mm. The MC-SART was also applied to measured patient cone-beam projections acquired with a linac-mounted CBCT system. Results from this patient data demonstrate the feasibility of MC-SART and showed qualitative image quality improvements compared to other state-of-the-art algorithms.

CONCLUSION

We have developed a simultaneous image reconstruction and motion model estimation method that uses Cone-beam computed tomography (CBCT) projections and respiratory surrogate measurements to reconstruct a high-quality reference image and motion model of a patient in treatment position. The method provided superior performance in both HU accuracy and positional accuracy compared to other existing methods. The resultant reference image and motion model can be combined with respiratory surrogate measurements to generate volumetric images representing patient anatomy at arbitrary time points.

CT brush and CancerZap!: two video games for computed tomography dose minimization

BACKGROUND

X-ray dose from computed tomography (CT) scanners has become a significant public health concern. All CT scanners spray x-ray photons across a patient, including those using compressive sensing algorithms. New technologies make it possible to aim x-ray beams where they are most needed to form a diagnostic or screening image. We have designed a computer game, CT Brush, that takes advantage of this new flexibility. It uses a standard MART algorithm (Multiplicative Algebraic Reconstruction Technique), but with a user defined dynamically selected subset of the rays. The image appears as the player moves the CT brush over an initially blank scene, with dose accumulating with every "mouse down" move. The goal is to find the "tumor" with as few moves (least dose) as possible.

RESULTS

We have successfully implemented CT Brush in Java and made it available publicly, requesting crowdsourced feedback on improving the open source code. With this experience, we also outline a "shoot 'em up game" CancerZap! for photon limited CT.

CONCLUSIONS

We anticipate that human computing games like these, analyzed by methods similar to those used to understand eye tracking, will lead to new object dependent CT algorithms that will require significantly less dose than object independent nonlinear and compressive sensing algorithms that depend on sprayed photons. Preliminary results suggest substantial dose reduction is achievable.

Improved digital breast tomosynthesis images using automated ultrasound

PURPOSE

Digital breast tomosynthesis (DBT) offers poor image quality along the depth direction. This paper presents a new method that improves the image quality of DBT considerably through the a priori information from automated ultrasound (AUS) images.

METHODS

DBT and AUS images of a complex breast-mimicking phantom are acquired by a DBT/AUS dual-modality system. The AUS images are taken in the same geometry as the DBT images and the gradient information of the in-slice AUS images is adopted into the new loss functional during the DBT reconstruction process. The additional data allow for new iterative equations through solving the optimization problem utilizing the gradient descent method. Both visual comparison and quantitative analysis are employed to evaluate the improvement on DBT images. Normalized line profiles of lesions are obtained to compare the edges of the DBT and AUS-corrected DBT images. Additionally, image quality metrics such as signal difference to noise ratio (SDNR) and artifact spread function (ASF) are calculated to quantify the effectiveness of the proposed method.

RESULTS

In traditional DBT image reconstructions, serious artifacts can be found along the depth direction (Z direction), resulting in the blurring of lesion edges in the off-focus planes parallel to the detector. However, by applying the proposed method, the quality of the reconstructed DBT images is greatly improved. Visually, the AUS-corrected DBT images have much clearer borders in both in-focus and off-focus planes, fewer Z direction artifacts and reduced overlapping effect compared to the conventional DBT images. Quantitatively, the corrected DBT images have better ASF, indicating a great reduction in Z direction artifacts as well as better Z resolution. The sharper line profiles along the Y direction show enhancement on the edges. Besides, noise is also reduced, evidenced by the obviously improved SDNR values.

CONCLUSIONS

The proposed method provides great improvement on the quality of DBT images. This improvement makes it easier to locate and to distinguish a lesion, which may help improve the accuracy of the diagnosis using DBT imaging.

Improving filtered backprojection reconstruction by data-dependent filtering

Filtered backprojection, one of the most widely used reconstruction methods in tomography, requires a large number of low-noise projections to yield accurate reconstructions. In many applications of tomography, complete projection data of high quality cannot be obtained, because of practical considerations. Algebraic methods tend to handle such problems better, but are computationally more expensive. In this paper, we introduce a new method that improves the filtered backprojection method by using a custom data-dependent filter that minimizes the projection error of the resulting reconstruction. We show that the computational cost of the new method is significantly lower than that of algebraic methods. Experiments on both simulation and experimental data show that the method is able to produce more accurate reconstructions than filtered backprojection based on popular static filters when presented with data with a limited number of projections or statistical noise present. Furthermore, the results show that the method produces reconstructions with similar accuracy to algebraic methods, but is faster at producing them. Finally, we show that the method can be extended to exploit certain forms of prior knowledge, improving reconstruction accuracy in specific cases.

Medical image processing on the GPU - past, present and future

Graphics processing units (GPUs) are used today in a wide range of applications, mainly because they can dramatically accelerate parallel computing, are affordable and energy efficient. In the field of medical imaging, GPUs are in some cases crucial for enabling practical use of computationally demanding algorithms. This review presents the past and present work on GPU accelerated medical image processing, and is meant to serve as an overview and introduction to existing GPU implementations. The review covers GPU acceleration of basic image processing operations (filtering, interpolation, histogram estimation and distance transforms), the most commonly used algorithms in medical imaging (image registration, image segmentation and image denoising) and algorithms that are specific to individual modalities (CT, PET, SPECT, MRI, fMRI, DTI, ultrasound, optical imaging and microscopy). The review ends by highlighting some future possibilities and challenges.

Evaluation of Algebraic Iterative Image Reconstruction Methods for Tetrahedron Beam Computed Tomography Systems

Tetrahedron beam computed tomography (TBCT) performs volumetric imaging using a stack of fan beams generated by a multiple pixel X-ray source. While the TBCT system was designed to overcome the scatter and detector issues faced by cone beam computed tomography (CBCT), it still suffers the same large cone angle artifacts as CBCT due to the use of approximate reconstruction algorithms. It has been shown that iterative reconstruction algorithms are better able to model irregular system geometries and that algebraic iterative algorithms in particular have been able to reduce cone artifacts appearing at large cone angles. In this paper, the SART algorithm is modified for the use with the different TBCT geometries and is tested using both simulated projection data and data acquired using the TBCT benchtop system. The modified SART reconstruction algorithms were able to mitigate the effects of using data generated at large cone angles and were also able to reconstruct CT images without the introduction of artifacts due to either the longitudinal or transverse truncation in the data sets. Algebraic iterative reconstruction can be especially useful for dual-source dual-detector TBCT, wherein the cone angle is the largest in the center of the field of view.