Evaluating real-time image reconstruction in diffuse optical tomography using physiologically realistic test data.

We describe a method to increase the speed of convergence for the simultaneous reconstruction of absorption and scattering images in Diffuse Optical Tomography (DOT). We used the diffusion approximation of the radiative transfer equation and the Finite Element Method (FEM) to solve the forward problem. The absorption and reduced scattering images are reconstructed by inverting the distribution of the moments of the time-dependent detected light flux. The inverse problem is solved with an optimization algorithm like ART or Conjugate Gradient. This ill-posed inverse problem can be simplified by using a priori knowledge of the studied objects. In this paper, we consider that DOT is a functional imaging technique that can be complemented by an anatomical imaging technique like Magnetic Resonance Imaging (MRI). We used anatomic information obtained from MRI as prior knowledge to compute optical absorption and scattering images. In a first step, MRI segmented images were only used to mesh our phantoms, with a finer resolution around boundaries. In a second step, we computed optical images with homogeneous properties from the segmented MRI image, in order to initialise our optimisation process. These two initialisations yield better reconstructed images. Reconstruction from simulated and experimental data will be presented.

We describe a method to increase the speed of convergence for the simultaneous reconstruction of absorption and scattering images in Diffuse Optical Tomography (DOT). We used the diffusion approximation of the radiative transfer equation and the Finite Element Method (FEM) to solve the forward problem. The absorption and reduced scattering images are reconstructed by inverting the distribution of the moments of the time-dependent detected light flux. The inverse problem is solved with an optimization algorithm like ART or Conjugate Gradient. This ill-posed inverse problem can be simplified by using <i>a priori</i> knowledge of the studied objects. In this paper, we consider that DOT is a functional imaging technique that can be complemented by an anatomical imaging technique like Magnetic Resonance Imaging (MRI). We used anatomic information obtained from MRI as prior knowledge to compute optical absorption and scattering images. In a first step, MRI segmented images were only used to mesh our phantoms, with a finer resolution around boundaries. In a second step, we computed optical images with homogeneous properties from the segmented MRI image, in order to initialise our optimisation process. These two initialisations yield better reconstructed images. Reconstruction from simulated and experimental data will be presented.

Effective optode configuration for the image reconstruction in diffuse optical tomography

Through application of spatio-temporal regularisation techniques, we demonstrate the real-time three-dimensional dynamic reconstruction of the optical properties of a hemispherical infant head phantom, with moving absorption and scattering targets.

Diffuse optical tomography (DOT) is a non-invasive imaging modality that uses near-infrared light to probe the optical properties of tissue. In conventionally used deterministic methods for DOT inversion, the measurement errors were not taken into account, resulting in unsatisfactory noise robustness and, consequently, affecting the DOT image reconstruction quality. In order to overcome this defect, an extended Kalman filter (EKF)-based DOT reconstruction algorithm was introduced first, which improved the reconstruction results by incorporating a priori information and measurement errors to the model. Further, to mitigate the instability caused by the ill-condition of the observation matrix in the tomographic imaging problem, a new, to the best of our knowledge, estimation algorithm was derived by incorporating Tikhonov regularization to the EKF method. To verify the effectiveness of the EKF algorithm and Tikhonov regularization-based EKF algorithm for DOT imaging, a series of numerical simulations and phantom experiments were conducted, and the experimental results were quantitatively evaluated and compared with two conventionally used deterministic methods involving the algebraic reconstruction technique and Levenberg-Marquardt algorithm. The results show that the two EKF-based algorithms can accurately estimate the location and size of the target, and the imaging accuracy and noise robustness are obviously improved. Furthermore, the Tikhonov regularization-based EKF obtained optimal parameter estimations, especially under the circumstance of low absorption contrast (1.2) and high noise level (10%).

The underdeterminedness of the inverse problems encountered in diffuse optical tomography (DOT) becomes especially severe when detecting breast cancers, because much more variables are needed to be reconstructed due to the big-size. With the addition of ill-condition caused by the diffusive nature of light propagation, the ill-posedness makes it very difficult to improve the image reconstruction. Fortunately, from the anatomy viewpoint, we have known that the cancer is distributed locally and only amounts to a small percentage of the whole breast. This makes it possible to employ the compressive sensing theory to mitigate the ill-posedness, based on the prior knowledge about the sparsity of the signal to be reconstructed. Specifically speaking, sparsity regularizations can be used in DOT to improve the image reconstruction under the premise that un-increase the number of measurements required in the reconstruction. In this paper, we primarily focus on comparing the performances of different kinds of L_p-norm-based regularizations in terms of theory and real effects, respectively. The numerical and phantom experiments have proven that the sparsity regularizations can dramatically improve the image reconstruction. Furthermore, as the p in the L_p-norm decreasing to zero, the solutions become sparser and the corresponding image quality gets higher, with smooth L₀-norm-based regularization providing the highest image quality.

We present high spatial density, multi-modal, parallel-plate Diffuse Optical Tomography (DOT) imaging systems for the purpose of breast tumor detection. One hybrid instrument provides time domain (TD) and continuous wave (CW) DOT at 64 source fiber positions. The TD diffuse optical spectroscopy with PMT- detection produces low-resolution images of absolute tissue scattering and absorption while the spatially dense array of CCD-coupled detector fibers (108 detectors) provides higher-resolution CW images of relative tissue optical properties. Reconstruction of the tissue optical properties, along with total hemoglobin concentration and tissue oxygen saturation, is performed using the TOAST software suite. Comparison of the spatially-dense DOT images and MR images allows for a robust validation of DOT against an accepted clinical modality. Additionally, the structural information from co-registered MR images is used as a spatial prior to improve the quality of the functional optical images and provide more accurate quantification of the optical and hemodynamic properties of tumors. We also present an optical-only imaging system that provides frequency domain (FD) DOT at 209 source positions with full CCD detection and incorporates optical fringe projection profilometry to determine the breast boundary. This profilometry serves as a spatial constraint, improving the quality of the DOT reconstructions while retaining the benefits of an optical-only device. We present initial images from both human subjects and phantoms to display the utility of high spatial density data and multi-modal information in DOT reconstruction with the two systems.

The relationship between conventional reference measurement correction and approximation error modelling in linear reconstruction in diffuse optical tomography is considered.

Malignant breast lesions manifest in a variety of shapes and configurations. Here we propose an algorithm utilizing anisotropic-tumor-prior guided diffuse optical tomography (DOT) reconstructions to recover an unknown tumor’s orientation and shape anisotropy.

Convergence rate calculation of simultaneous iterative reconstruction technique algorithm for diffuse optical tomography image reconstruction: A feasibility study

The image reconstruction for Fourier-domain optical coherence tomography (FD-OCT) could be achieved by iterative methods, which offer a more accurate estimation than the traditional inverse discrete Fourier transform (IDFT) reconstruction. However, the existing iterative methods are mostly A-line-based and are developed on CPU, which causes slow reconstruction. Besides, A-line-based reconstruction makes the iterative methods incompatible with most existing image-level image processing techniques. In this paper, we proposed an iterative method that enables B-scan-based OCT image reconstruction, which has three major advantages: (1) Large-scale parallelism of the OCT dataset is achieved by using GPU acceleration. (2) A novel image-level cross-domain regularizer was developed, such that the image processing could be performed simultaneously during the image reconstruction; an enhanced image could be directly generated from the OCT interferogram. (3) The scalability of the proposed method was demonstrated for 3D OCT image reconstruction. Compared with the state-of-the-art (SOTA) iterative approaches, the proposed method achieves higher image quality with reduced computational time by orders of magnitude. To further show the image enhancement ability, a comparison was conducted between the proposed method and the conventional workflow, in which an IDFT reconstructed OCT image is later processed by a total variation-regularized denoising algorithm. The proposed method can achieve a better performance evaluated by metrics such as signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR), while the speed is improved by more than 30 times. Real-time image reconstruction at more than 20 B-scans per second was realized with a frame size of 4096 (axial) × 1000 (lateral), which showcases the great potential of the proposed method in real-world applications.

Current fluorescence diffuse optical tomography (fDOT) systems can provide large data sets and, in addition, the unknown parameters to be estimated are so numerous that the sensitivity matrix is too large to store. Alternatively, iterative methods can be used, but they can be extremely slow at converging when dealing with large matrices. A few approaches suitable for the reconstruction of images from very large data sets have been developed. However, they either require explicit construction of the sensitivity matrix, suffer from slow computation times, or can only be applied to restricted geometries. We introduce a method for fast reconstruction in fDOT with large data and solution spaces, which preserves the resolution of the forward operator whilst compressing its representation. The method does not require construction of the full matrix, and thus allows storage and direct inversion of the explicitly constructed compressed system matrix. The method is tested using simulated and experimental data. Results show that the fDOT image reconstruction problem can be effectively compressed without significant loss of information and with the added advantage of reducing image noise.

This paper introduces an improved diffusion model which is accurate and fast in computation for the cases of mu(a)/mu'(s) < 0.07 as good as the conventional diffusion model for the cases of mu(a)/mu'(s) < 0.007 for surface measurement, hence more suitable than the conventional model to be the forward model used in the image reconstruction in the diffuse optical tomography. Deviation of the diffusion approximation (DA) on the medium surface is first studied in the Monte Carlo (MC) diffusion hybrid model for reflectance setup. A modification of DA and an improved MC diffusion hybrid model based on this modified DA are introduced. Numerical tests show that for media with relatively strong absorption the present modified diffusion approach can reduce the surface deviation significantly in both the hybrid and pure diffusion model, and consumes nearly no more computation time than the conventional diffusion approach.

We attempt to develop a systematic scheme through adopting high-pass filtering (HPF) to well resolve value-preserved images such as medical images. Our approach is derived from the Poisson maximum a posteriori superresolution algorithm employing the HP filters, where four filters are considered such as two low-pass-filter-combination based filters, wavelet filter, and negative-oriented Laplacian HP filter. The proposed approach is incorporated into the procedure of finite-element-method (FEM)-based image reconstruction for diffuse optical tomography in the direct current domain, posterior to each iteration without altering the original FEM modeling. This approach is justified with various HPF for different cases that breast-like phantoms embedded with two or three inclusions that imitate tumors are employed to examine the resolution performances under certain extreme conditions. The proposed approach to enhancing image resolution is evaluated for all tested cases. A qualitative investigation of reconstruction performance for each case is presented. Following this, we define a set of measures on the quantitative evaluation for a range of resolutions including separation, size, contrast, and location, thereby providing a comparable evaluation to the visual quality. The most satisfactory result is obtained by using the wavelet HP filter, and it successfully justifies our proposed scheme.

Near infrared (NIR) diffused optical tomography (DOT) is emerging as a potential tool of non-invasively diagnosing woman breast cancers, neonatal brain hypoxia, and other human organ diseases. The intensive and worldwide investigations in theory and experiment have revealed the possibility of NIR DOT in providing both anatomical and functional information of biological tissue simultaneously, which is important for distinguishing between healthy and diseased tissues, such as benign and malignant tumors. In this paper, our recent DOT experiments on human lower legs and forearms are presented using our time-resolved measuring system and image reconstruction algorithm based on the modified generalized pulse spectrum technique. It was found that the image quality in the experiments, including both the spatial resolution and the quantitativeness of the targets, was rather poor, and the interior blood vessels undisclosed in the absorption images. To clarify this issue, the influences of target contrast and size on the image reconstruction were investigated with simulated data. We have accordingly obtained the following observations: the quantitativeness of the reconstructed optical properties was prone to be spoiled by the small size ratio and high contrast of the interior targets (such as blood vessels) to the background, and the incompleteness of information embedded in the featured data-types, in addition to the experimental noise, evidently answers for the degradation of the spatial resolution and quantitativeness. It was shown in a further simulative investigation that the image quality could be substantially improved by making full use of the time-resolved data.