Diffuse optical tomography uses near infrared (NIR) light as the probing media to recover the distributions of tissue optical properties. It has a potential to become an adjunct imaging modality for breast and brain imaging, that is capable of providing functional information of the tissue under investigation. As NIR light propagation in the tissue is dominated by scattering, the image reconstruction problem (inverse problem) tends to be non-linear and ill-posed, requiring usage of advanced computational methods to compensate this. Traditional image reconstruction methods in diffuse optical tomography employ l2 –norm based regularization, which is known to remove high frequency noises in the re-constructed images and make them appear smooth. The recovered contrast in the reconstructed image in these type of methods are typically dependent on the iterative nature of the method employed, in which the non-linear iterative technique is known to perform better in comparison to linear techniques. The usage of non-linear iterative techniques in the real-time, especially in dynamical imaging, becomes prohibitive due to the computational complexity associated with them. In the rapid dynamic diffuse optical imaging, assumption of a linear dependency in the solutions between successive frames results in a linear inverse problem. This new frame work along with the l1–norm based regularization can provide better robustness to noise and results in a better contrast recovery compared to conventional l2 –based techniques. Moreover, it is shown that the proposed l1-based technique is computationally efficient compared to its counterpart(l2 –based one). The proposed framework requires a reasonably close estimate of the actual solution for the initial frame and any suboptimal estimate leads to erroneous reconstruction results for the subsequent frames. Modern diffuse optical imaging systems are multi-modal in nature, where diffuse optical imaging is combined with traditional imaging modalities such as MRI, CT, and Ultrasound. A novel approach that can more effectively use the structural information provided by the traditional imaging modalities in these scenarios is introduced, which is based on prior image constrained- l1 minimization scheme. This method has been motivated by the recent progress in the sparse image reconstruction techniques. It is shown that the- l1 based frame work is more effective in terms of localizing the tumor region and recovering the optical property values both in numerical and gelatin phantom cases compared to the traditional methods that use structural information.

Diffuse Optical Tomography is a diagnostic imaging modality where optical parameters such as absorption coefficient, scattering coefficient and refractive index distributions are recovered to form the internal tissue metabolic image. Near-infrared (NIR) light has the potential to be used as a noninvasive means of diagnostic imaging within the human breast. Due to the diffusive nature of light in tissue, computational model-based methods are required for functional imaging. The main goal is to recover the spatial variation of optical properties which shed light on the different metabolic states of tissue and tissue like media. This thesis addresses the issue of quantitative recovery of optical properties of tissue-mimicking phantom and pork tissue using diffuse optical tomography (DOT). The main contribution of the present work is the development of robust, efficient and fast optical property reconstruction algorithms for a direct 3-D DOT imaging system. There are both theoretical and experimental contributions towards the development of an imaging system and procedures to minimize accurate data collection time, overall system automation as well as development of computational algorithms. In nurturing the idea of imaging using NIR light into a fully developed direct 3-D imaging system, challenges from the theoretical and computational aspects have to be met. The recovery of the optical property distribution in the interior of the object from the often noisy boundary measurements on light, is an ill-posed ( and nonlinear) problem. This is particularly true, when one is interested in a direct 3-D image reconstruction instead of the often employed stacking of 2-D cross-sections obtained from solving a set of 2-D DOT problems. In order to render the DOT, a useful diagnostic imaging tool and a robust reconstruction procedure giving accurate and reliable parameter recovery in the scenario, where the number of unknowns far outnumbers the number of independent data sets that can be gathered (for example, the direct 3-D recovery mentioned earlier) is essential. Here, the inversion problem is often solved through iterative methods based on nonlinear optimization for the minimization of a data-model misfit function. An interesting development in this direction has been the development of Broyden’ s and adjoint Broyden’ s methods that avoids direct Jacobian computation in each iteration thereby making the full 3-D a reality. Conventional model based iterative image reconstruction (MoBIIR) algorithm uses Newton’ s and it’s variant methods, where it required repeated evaluation of whole Jacobian, which consumes bulk time in reconstruction process. The explicit secant and adjoint information based fast 2-D/3-D image reconstruction algorithms without repeated evaluation of the Jacobian is proposed in diffuse optical tomography, where the computational time has been decreased many folds by updating the Jacobian successively through low rank update. An alternative route to the iterative solution is attempted by introducing an artificial dynamics in the system and treating the steady-state response of the artificially evolving dynamical system as a solution. The objective is to consider a novel family of pseudo-dynamical 2-D and 3-D systems whose numerical integration in time provides an asymptotic solution to the inverse problem at hand. We convert Gauss-Newton’ s equation for updates into a pseudo-dynamical (PD) form by explicitly adding a time derivative term. As the pseudo-time integration schemes do not need such explicit matrix inversion and depending on the pseudo-time step size, provides for a layer of regularization that in turn helps in superior quality of 2-D and 3-D image reconstruction. A cost effective frequency domain Matlab based 2-D/3-D automated imaging system is designed and built. The complete instrumentation (including PC-based control software) has been developed using a single modulated laser source (wavelength 830nm) and a photo-multiplier tube (PMT). The source and detector fiber change their positions dynamically allowing us to gather data at multiple source and detector locations. The fiber positions are adjusted on the phantom surface automatically for scanning variable size phantoms. A heterodyning scheme was used for reading out the measurement using a lock-in-amplifier. The Matlab program carries out sequence of actions such as instrument control, data acquisition, data organization, data calibration and reconstruction of image. The Gauss-Newton’ s, Broyden’ s, adjoint Broyden’ s and pseudo-time integration algorithms are evaluated using the simulation data as well as data from the experimental DOT system. Validation of the system and the reconstruction algorithms were carried out on a real tissue, a pork tissue with an embedded fat inhomogeneity. The results were found to match the known parameters closely.

Diffuse optical tomography uses near infrared (NIR) light as the probing media to re-cover the distributions of tissue optical properties with an ability to provide functional information of the tissue under investigation. As NIR light propagation in the tissue is dominated by scattering, the image reconstruction problem (inverse problem) is non-linear and ill-posed, requiring usage of advanced computational methods to compensate this. Diffuse optical image reconstruction problem is always rank-deficient, where finding the independent measurements among the available measurements becomes challenging problem. Knowing these independent measurements will help in designing better data acquisition set-ups and lowering the costs associated with it. An optimal measurement selection strategy based on incoherence among rows (corresponding to measurements) of the sensitivity (or weight) matrix for the near infrared diffuse optical tomography is proposed. As incoherence among the measurements can be seen as providing maximum independent information into the estimation of optical properties, this provides high level of optimization required for knowing the independency of a particular measurement on its counterparts. The utility of the proposed scheme is demonstrated using simulated and experimental gelatin phantom data set comparing it with the state-of-the-art methods. The traditional image reconstruction methods employ ℓ2-norm in the regularization functional, resulting in smooth solutions, where the sharp image features are absent. The sparse recovery methods utilize the ℓp-norm with p being between 0 and 1 (0 ≤ p1), along with an approximation to utilize the ℓ0-norm, have been deployed for the reconstruction of diffuse optical images. These methods are shown to have better utility in terms of being more quantitative in reconstructing realistic diffuse optical images compared to traditional methods. Utilization of ℓp-norm based regularization makes the objective (cost) function non-convex and the algorithms that implement ℓp-norm minimization utilizes approximations to the original ℓp-norm function. Three methods for implementing the ℓp-norm were con-sidered, namely Iteratively Reweigthed ℓ1-minimization (IRL1), Iteratively Reweigthed Least-Squares (IRLS), and Iteratively Thresholding Method (ITM). These results in-dicated that IRL1 implementation of ℓp-minimization provides optimal performance in terms of shape recovery and quantitative accuracy of the reconstructed diffuse optical tomographic images. Photoacoustic tomography (PAT) is an emerging hybrid imaging modality combining optics with ultrasound imaging. PAT provides structural and functional imaging in diverse application areas, such as breast cancer and brain imaging. A model-based iterative reconstruction schemes are the most-popular for recovering the initial pressure in limited data case, wherein a large linear system of equations needs to be solved. Often, these iterative methods requires regularization parameter estimation, which tends to be a computationally expensive procedure, making the image reconstruction process to be performed off-line. To overcome this limitation, a computationally efficient approach that computes the optimal regularization parameter is developed for PAT. This approach is based on the least squares-QR (LSQR) decomposition, a well-known dimensionality reduction technique for a large system of equations. It is shown that the proposed framework is effective in terms of quantitative and qualitative reconstructions of initial pressure distribution.