Regularization Techniques for Inverse Problem in DOT Applications

Diffuse optical tomography (DOT) is an emerging diagnostic technique which uses near-infra-red light to investigate the optical coefficients distribution in biological tissues. The surface of the tissue is illuminated by light sources, then the outgoing light is measured by detectors placed at various locations on the surface itself. In order to reconstruct the optical coefficients, a mathematical model of light propagation is employed: such model leads to the minimization of the discrepancy between the detected data and the corresponding theoretical field. Due to severe ill-conditioning, regularization techniques are required: common procedures consider mainly \u2113 1-norm (LASSO) and \u2113 2-norm (Tikhonov) regularization. In the present work we investigate two original approaches in this context: The elastic-net regularization, previously used in machine learning problems, and the Bregman procedure. Numerical experiments are performed on synthetic 2D geometries and data, to evaluate the performance of these approaches. The results show that these techniques are indeed suitable choices for practical applications, where DOT is used as a cheap, first-level and almost real-Time screening technique for breast cancer detection

Reconstruction Method for Optical Tomography Based on the Linearized Bregman Iteration with Sparse Regularization

Optical molecular imaging is a promising technique and has been widely used in physiology, and pathology at cellular and molecular levels, which includes different modalities such as bioluminescence tomography, fluorescence molecular tomography and Cerenkov luminescence tomography. The inverse problem is ill-posed for the above modalities, which cause a nonunique solution. In this paper, we propose an effective reconstruction method based on the linearized Bregman iterative algorithm with sparse regularization (LBSR) for reconstruction. Considering the sparsity characteristics of the reconstructed sources, the sparsity can be regarded as a kind of a priori information and sparse regularization is incorporated, which can accurately locate the position of the source. The linearized Bregman iteration method is exploited to minimize the sparse regularization problem so as to further achieve fast and accurate reconstruction results. Experimental results in a numerical simulation and in vivo mouse demonstrate the effectiveness and potential of the proposed method

Improvements in the robustness and accuracy of bioluminescence tomographic reconstructions of distributed sources within small animals

High quality three-dimensional bioluminescence tomographic (BLT) images, if available, would constitute a major advance and provide much more useful information than the two-dimensional bioluminescence images that are frequently used today. To-date, high quality BLT images have not been available, largely because of the poor quality of the data being input into the reconstruction process. Many significant confounds are not routinely corrected for and the noise in this data is unnecessarily large and poorly distributed. Moreover, many of the design choices affecting image quality are not well considered, including choices regarding the number and type of filters used when making multispectral measurements and choices regarding the frequency and uniformity of the sampling of both the range and domain of the BLT inverse problem. Finally, progress in BLT image quality is difficult to gauge owing to a lack of realistic gold-standard references that engage the full complexity and uncertainty within a small animal BLT imaging experiment. Within this dissertation, I address all of these issues. I develop a Cerenkov-based gold-standard wherein a Positron Emission Tomography (PET) image can be used to gauge improvements in the accuracy of BLT reconstruction algorithms. In the process of creating this reference, I discover and describe corrections for several confounds that if left uncorrected would introduce artifacts into the BLT images. This includes corrections for the angle of the animal’s skin surface relative to the camera, for the height of each point on the skin surface relative to the focal plane, and for the variation in bioluminescence intensity as a function of luciferin concentration over time. Once applied, I go on to derive equations and algorithms that when employed are able to minimize the noise in the final images under the constraints of a multispectral BLT data acquisition. These equations and algorithms allow for an optimal choice of filters to be made and for the acquisition time to be optimally distributed among those filtered measurements. These optimizations make use of Barrett’s and Moore-Penrose pseudoinverse matrices which also come into play in a paradigm I describe that can be used to guide choices regarding sampling of the domain and range