Fluorescence molecular tomography using a two-step three-dimensional shape-based reconstruction with graphics processing unit acceleration

Generalized total variation iterative constraint strategy in limited angle optical diffraction tomography

Due to incompleteness of input data inherent to Limited Angle Tomography (LAT), specific additional constraints are usually employed to suppress image artifacts. In this work we demonstrate a new two-stage regularization strategy, named Generalized Total Variation Iterative Constraint (GTVIC), dedicated to semi-piecewise-constant objects. It has been successfully applied as a supplementary module for two different reconstruction algorithms: an X-ray type solver and a diffraction-wise solver. Numerical tests performed on a detailed phantom of a biological cell under conical illumination pattern show significant reduction of axial blurring in the reconstructed refractive index distribution after GTVIC is added. Analogous results were obtained with experimental data.

A reconstruction approach in wavelet domain for fluorescent molecular tomography via rotated sources illumination

Background

Fluorescent molecular tomography (FMT) aims at reconstructing the spatial map of optical and fluorescence parameters from fluence measurements. Basically, solving large-scale matrix equations is computationally expensive for image reconstruction of FMT. Despite the reconstruction quality can be improved with more sources, it may result in higher computational costs for reconstruction. This article presents a novel method in the wavelet domain with rotated sources illumination.

Methods

We use the finite element method for the computation of the forward model. The global inverse problem is solved based on wavelet in conjunction with principal component analysis. The iterative reconstruction is implemented with sources rotated in a certain angle. The original excitation light sources are used to reconstruct the image in the first iteration. Then, upon the sources are rotated by a certain angle, they are employed for the next iteration of reconstruction.

Results

Simulation results demonstrate that our method can considerably reduce the time taken for the computation of inverse problem in FMT. Furthermore, the approach proposed is also shown to largely outperform the traditional method in terms of the precision of inverse solutions.

Conclusions

Our method has the capability to locate the inclusions. The proposed method can significantly speed up the reconstruction process with the high reconstruction quality.

Full domain-decomposition scheme for diffuse optical tomography of large-sized tissues with a combined CPU and GPU parallelization

The common approach to diffuse optical tomography is to solve a nonlinear and ill-posed inverse problem using a linearized iteration process that involves repeated use of the forward and inverse solvers on an appropriately discretized domain of interest. This scheme normally brings severe computation and storage burdens to its applications on large-sized tissues, such as breast tumor diagnosis and brain functional imaging, and prevents from using the matrix-fashioned linear inversions for improved image quality. To cope with the difficulties, we propose in this paper a parallelized full domain-decomposition scheme, which divides the whole domain into several overlapped subdomains and solves the corresponding subinversions independently within the framework of the Schwarz-type iterations, with the support of a combined multicore CPU and multithread graphics processing unit (GPU) parallelization strategy. The numerical and phantom experiments both demonstrate that the proposed method can effectively reduce the computation time and memory occupation for the large-sized problem and improve the quantitative performance of the reconstruction.

High-performance fluorescence molecular tomography through shape-based reconstruction using spherical harmonics parameterization

Fluorescence molecular tomography in the near-infrared region is becoming a powerful modality for mapping the three-dimensional quantitative distributions of fluorochromes in live small animals. However, wider application of fluorescence molecular tomography still requires more accurate and stable reconstruction tools. We propose a shape-based reconstruction method that uses spherical harmonics parameterization, where fluorophores are assumed to be distributed as piecewise constants inside disjointed subdomains and the remaining background. The inverse problem is then formulated as a constrained nonlinear least-squares problem with respect to shape parameters, which decreases ill-posedness because of the significantly reduced number of unknowns. Since different shape parameters contribute differently to the boundary measurements, a two-step and modified block coordinate descent optimization algorithm is introduced to stabilize the reconstruction. We first evaluated our method using numerical simulations under various conditions for the noise level and fluorescent background; it showed significant superiority over conventional voxel-based methods in terms of the spatial resolution, reconstruction accuracy with regard to the morphology and intensity, and robustness against the initial estimated distribution. In our phantom experiment, our method again showed better spatial resolution and more accurate intensity reconstruction. Finally, the results of an in vivo experiment demonstrated its applicability to the imaging of mice.

Three-dimensional photoacoustic tomography based on graphics-processing-unit-accelerated finite element method

Compared with commonly used analytical reconstruction methods, the frequency-domain finite element method (FEM) based approach has proven to be an accurate and flexible algorithm for photoacoustic tomography. However, the FEM-based algorithm is computationally demanding, especially for three-dimensional cases. To enhance the algorithm's efficiency, in this work a parallel computational strategy is implemented in the framework of the FEM-based reconstruction algorithm using a graphic-processing-unit parallel frame named the "compute unified device architecture." A series of simulation experiments is carried out to test the accuracy and accelerating effect of the improved method. The results obtained indicate that the parallel calculation does not change the accuracy of the reconstruction algorithm, while its computational cost is significantly reduced by a factor of 38.9 with a GTX 580 graphics card using the improved method.