Frequency-domain optical tomographic image reconstruction algorithm with the simplified spherical harmonics (SP3) light propagation model

Unrolled-DOT: an interpretable deep network for diffuse optical tomography

SIGNIFICANCE

Imaging through scattering media is critical in many biomedical imaging applications, such as breast tumor detection and functional neuroimaging. Time-of-flight diffuse optical tomography (ToF-DOT) is one of the most promising methods for high-resolution imaging through scattering media. ToF-DOT and many traditional DOT methods require an image reconstruction algorithm. Unfortunately, this algorithm often requires long computational runtimes and may produce lower quality reconstructions in the presence of model mismatch or improper hyperparameter tuning.

AIM

We used a data-driven unrolled network as our ToF-DOT inverse solver. The unrolled network is faster than traditional inverse solvers and achieves higher reconstruction quality by accounting for model mismatch.

APPROACH

Our model "Unrolled-DOT" uses the learned iterative shrinkage thresholding algorithm. In addition, we incorporate a refinement U-Net and Visual Geometry Group (VGG) perceptual loss to further increase the reconstruction quality. We trained and tested our model on simulated and real-world data and benchmarked against physics-based and learning-based inverse solvers.

RESULTS

In experiments on real-world data, Unrolled-DOT outperformed learning-based algorithms and achieved over 10× reduction in runtime and mean-squared error, compared to traditional physics-based solvers.

CONCLUSION

We demonstrated a learning-based ToF-DOT inverse solver that achieves state-of-the-art performance in speed and reconstruction quality, which can aid in future applications for noninvasive biomedical imaging.

Three-dimensional iterative reconstruction of pulsed radiation sources using spherical harmonic decomposition

Neutron and x-ray imaging are essential ways to diagnose a pulsed radiation source. The three-dimensional (3D) intensity distribution reconstructed from two-dimensional (2D) radiation images can significantly promote research regarding the generation and variation mechanisms of pulsed radiation sources. Only a few (≤5) projected images at one moment are available due to the difficulty in building imaging systems for high-radiation-intensity and short-pulsed sources. The reconstruction of a 3D source with a minimal number of 2D images is an ill-posed problem that leads to severe structural distortions and artifacts of the image reconstructed by conventional algorithms. In this paper, we present an iterative method to reconstruct a 3D source using spherical harmonic decomposition. Our algorithm improves the representation ability of spherical harmonic decomposition for 3D sources by enlarging the order of the expansion, which is limited in current analytical reconstruction algorithms. Prior knowledge of the source can be included to obtain a reasonable solution. Numerical simulations demonstrate that the reconstructed image quality of the iterative algorithm is better than that of the analytical algorithm. The iterative method can suppress the effect of noise in the integral projection image and has better robustness and adaptability than the analytical method.

Diagnostic Evaluation of Rheumatoid Arthritis (RA) in Finger Joints Based on the Third-Order Simplified Spherical Harmonics (SP3) Light Propagation Model

This work focuses on the evaluation of third-order simplified spherical harmonics (SP3) model-based image reconstruction with respect to its clinical utility to diagnose rheumatoid arthritis (RA). The existing clinical data of 219 fingers was reconstructed for both absorption and scattering maps in fingers by using the reduced-Hessian sequential quadratic programming (rSQP) algorithm that employs the SP3 model of light propagation. The k-fold cross validation method was used for feature extraction and classification of SP3-based tomographic images. The performance of the SP3 model was compared to the DE and ERT models with respect to diagnostic accuracy and computational efficiency. The results presented here show that the SP3 model achieves clinically relevant sensitivity (88%) and specificity (93%) that compare favorably to the ERT while maintaining significant computational advantage over the ERT (i.e., the SP3 model is 100 times faster than the ERT). Furthermore, it is also shown that the SP3 is similar in speed but superior in diagnostic accuracy to the DE. Therefore, it is expected that the method presented here can greatly aid in the early diagnosis of RA with clinically relevant accuracy in near real-time at a clinical setting.

Ultrafast and Ultrahigh-Resolution Diffuse Optical Tomography for Brain Imaging with Sensitivity Equation based Noniterative Sparse Optical Reconstruction (SENSOR)

We introduce a novel image reconstruction method for time-resolved diffuse optical tomography (DOT) that yields submillimeter resolution in less than a second. This opens the door to high-resolution real-time DOT in imaging of the brain activity. We call this approach the sensitivity equation based noniterative sparse optical reconstruction (SENSOR) method. The high spatial resolution is achieved by implementing an asymptotic l ₀-norm operator that guarantees to obtain sparsest representation of reconstructed targets. The high computational speed is achieved by employing the nontruncated sensitivity equation based noniterative inverse formulation combined with reduced sensing matrix and parallel computing. We tested the new method with numerical and experimental data. The results demonstrate that the SENSOR algorithm can achieve 1 mm³ spatial-resolution optical tomographic imaging at depth of ∼60 mean free paths (MFPs) in 20∼30 milliseconds on an Intel Core i9 processor.

High Resolution, Deep Imaging Using Confocal Time-of-Flight Diffuse Optical Tomography

Light scattering by tissue severely limits how deep beneath the surface one can image, and the spatial resolution one can obtain from these images. Diffuse optical tomography (DOT) is one of the most powerful techniques for imaging deep within tissue - well beyond the conventional  ∼ 10-15 mean scattering lengths tolerated by ballistic imaging techniques such as confocal and two-photon microscopy. Unfortunately, existing DOT systems are limited, achieving only centimeter-scale resolution. Furthermore, they suffer from slow acquisition times and slow reconstruction speeds making real-time imaging infeasible. We show that time-of-flight diffuse optical tomography (ToF-DOT) and its confocal variant (CToF-DOT), by exploiting the photon travel time information, allow us to achieve millimeter spatial resolution in the highly scattered diffusion regime ( mean free paths). In addition, we demonstrate two additional innovations: focusing on confocal measurements, and multiplexing the illumination sources allow us to significantly reduce the measurement acquisition time. Finally, we rely on a novel convolutional approximation that allows us to develop a fast reconstruction algorithm, achieving a 100× speedup in reconstruction time compared to traditional DOT reconstruction techniques. Together, we believe that these technical advances serve as the first step towards real-time, millimeter resolution, deep tissue imaging using DOT.