Image reconstruction in fluorescence molecular tomography with sparsity-initialized maximum-likelihood expectation maximization

Cone-beam X-ray luminescence computed tomography based on MLEM with adaptive FISTA initial image

BACKGROUND AND OBJECTIVE

As an emerging dual-mode optical molecular imaging, cone-beam X-ray luminescence computed tomography (CB-XLCT) has shown potential in early tumor diagnosis and other applications with increased depth and little autofluorescence. However, due to the low transfer efficiency of PNPs to convert X-ray energy to visible or near-infrared (NIR) light and X-ray dose limitation, the signal to noise ratio of projections is quite low, making the quality of CB-XLCT relatively poor.

METHODS

To improve the reconstruction quality of low-counts CB-XLCT imaging, an adaptive reconstruction algorithm (named ADFISTA-MLEM) based on the maximum likelihood expectation estimation (MLEM) framework is proposed for CB-XLCT reconstruction from Poisson distributed projections. In the proposed framework, the image reconstructed by fast iterative shrinkage-thresholding algorithm (FISTA) is used as the initial image for MLEM iterations to improve reconstruction accuracy, in which both the projection noise model and the sparsity constraint of the image could be considered. For relative quantitative imaging, a specific normalization is applied to the projection data and system matrix. To determine the hyperparameter of FISTA, which may be different for different projections, an adaptive strategy (ADFISTA) is then designed for adaptive update of the hyperparameter with reconstructed image in each iteration.

RESULTS AND CONCLUSIONS

Results from numerical simulations and phantom experiments indicate that the proposed framework can obtain superior reconstruction accuracy in terms of contrast to noise ratio and shape similarity. In addition, high intensity-concentration linearity between different probe targets indicates its potential for quantitative CB-XLCT imaging.

Model-based system matrix for iterative reconstruction in sub-diffuse angular-domain fluorescence optical projection tomography

This work concerns a fluorescence optical projection tomography system for low scattering tissue, like lymph nodes, with angular-domain rejection of highly scattered photons. In this regime, filtered backprojection (FBP) image reconstruction has been shown to provide reasonable quality images, yet here a comparison of image quality between images obtained by FBP and iterative image reconstruction with a Monte Carlo generated system matrix, demonstrate measurable improvements with the iterative method. Through simulated and experimental phantoms, iterative algorithms consistently outperformed FBP in terms of contrast and spatial resolution. Moreover, when projection number was reduced, in order to reduce total imaging time, iterative reconstruction suppressed artifacts that hampered the performance of FBP reconstruction (structural similarity of the reconstructed images with "truth" was improved from 0.15 ± 1.2 × 10^-3 to 0.66 ± 0.02); and although the system matrix was generated for homogenous optical properties, when heterogeneity (62.98 cm^-1 variance in µ_s ) was introduced to simulated phantoms, the results were still comparable (structural similarity homo: 0.67 ± 0.02 vs hetero: 0.66 ± 0.02).

In Vivo Visualization of Eye Vasculature Using Super-Resolution Ultrasound Microvessel Imaging

OBJECTIVE

The choroidal vessels, which supply oxygen and nutrient to the retina, may play a pivotal role in eye disease pathogenesis such as diabetic retinopathy and glaucoma. In addition, the retrobulbar circulation that feeds the choroid shows an important pathophysiologic role in myopia and degenerative myopia. Owing to the light-absorbing retinal pigment epithelium (RPE) and optically opaque sclera, choroidal and retrobulbar vasculature were difficult to be observed using clinically accepted optical coherence tomography angiography (OCT-A) technique. Here, we have developed super-resolution ultrasound microvessel imaging technique to visualize the deep ocular vasculature.

METHODS

An 18-MHz linear array transducer with compounding plane wave imaging technique and contrast agent - microbubble was implemented in this study. The centroid intensity of each microbubble was detected using image deconvolution algorithm with spatially variant point spread function, and then accumulated in successive frames in order to reconstruct microvasculature. The image deconvolution technique was first evaluated in a simulation study and experimental flow phantoms. The performance was then validated on normal rabbit eyes in vivo.

RESULTS

The image deconvolution based super-resolution ultrasound microvessel imaging technique shows good performance on either simulation study or flow phantoms. In vivo rabbit eye study indicated that the micron-level choroidal and retrobulbar vessels around the optic nerve head were successfully reconstructed in multiple 2D views and 3D volume imaging.

CONCLUSIONS

Our results demonstrate the capability of using super-resolution ultrasound microvessel imaging technique to image the microvasculature of the posterior pole of the eye. This efficient approach can potentially lead to a routinely performed diagnostic procedure in the field of ophthalmology.

Reconstruction for Fluorescence Molecular Tomography via Adaptive Group Orthogonal Matching Pursuit

OBJECTIVE

Fluorescence molecular tomography (FMT) is a promising medical imaging technology aimed at the non-invasive, specific, and sensitive detection of the distribution of fluorophore. Conventional sparsity prior-based methods of FMT commonly face problems such as over-sparseness, spatial discontinuity, and poor robustness, due to the neglect of the interrelation within the local subspace. To address this, we propose an adaptive group orthogonal matching pursuit (AGOMP) method.

METHODS

AGOMP is based on a novel local spatial-structured sparse regularization, which leverages local spatial interrelations as group sparsity without the hard prior of the tumor region. The adaptive grouped subspace matching pursuit method was adopted to enhance the interrelatedness of elements within a group, which alleviates the over-sparsity problem to some extent and improves the accuracy, robustness, and morphological similarity of FMT reconstruction. A series of numerical simulation experiments, based on digital mouse with both one and several tumors, were conducted, as well as in vivo mouse experiments.

RESULTS

The results demonstrated that the proposed AGOMP method achieved better location accuracy, fluorescent yield reconstruction, relative sparsity, and morphology than state-of-the-art methods under complex conditions for levels of Gaussian noise ranging from 5-25%. Furthermore, the in vivo mouse experiments demonstrated the practical application of FMT with AGOMP.

CONCLUSION

The proposed AGOMP can improve the accuracy and robustness for FMT reconstruction in biomedical application.

Sensitivity Uniformity Ratio as a New Index to Optimize the Scanning Geometry for Fluorescent Molecular Tomography

Background:

Molecular fluorescence imaging is widely used as a noninvasive method to study the cellular and molecular mechanisms. In the optical imaging system, the sensitivity is the change of the intensity received by the detector for the changed optical characteristics (fluorescence) in each sample point. Sensitivity could be considered as a function of imaging geometry. A favor imaging system has a uniform and high-sensitivity coefficient for each point of the sample. In this study, a new parameter was proposed which the optimal angle between the source and detector could be determined based on this parameter.

Methods:

For evaluation of the new method, a two-dimensional mesh with a radius of 20 mm and 5133 nodes was built. In each reconstruction, 0.5-mm fluorescence heterogeneity with a contrast-to-purpose ratio of fluorescence yield of 10 was randomly added at different points of the sample. The source and the detector were simulated in different geometric conditions. The calculations were performed using the NIRFAST and MATLAB software. The relationship between mean squared error (MSE) and sensitivity uniformity ratio (SUR) was evaluated using the correlation coefficient. Finally, based on the new index, an optimal geometrical strategy was introduced.

Results:

There was a negative correlation coefficient (R = −0.78) with the inverse relationship between the SUR and MSE indices. The reconstructed images showed that the better image quality achieved using the optimal geometry for all scanning depths. For the conventional geometry, there is an artifact in the opposite side of the inhomogeneity at the shallow depths, which has been eliminated in the reconstructed images achieved using the optimal geometry.

Conclusion:

The SUR is a powerful computational tool which could be used to determine the optimal angles between the source and detector for each scanning depth.