Ultrasound-modulated optical tomography: direct recovery of elasticity distribution from experimentally measured intensity autocorrelation.

Based on an ultrasound-modulated optical tomography experiment, a direct, quantitative recovery of Young's modulus (E) is achieved from the modulation depth (M) in the intensity autocorrelation. The number of detector locations is limited to two in orthogonal directions, reducing the complexity of the data gathering step whilst ensuring against an impoverishment of the measurement, by employing ultrasound frequency as a parameter to vary during data collection. The M and E are related via two partial differential equations. The first one connects M to the amplitude of vibration of the scattering centers in the focal volume and the other, this amplitude to E. A (composite) sensitivity matrix is arrived at mapping the variation of M with that of E and used in a (barely regularized) Gauss-Newton algorithm to iteratively recover E. The reconstruction results showing the variation of E are presented.

A high precision instrument to measure angular and binocular deviation introduced by aircraft windscreens by using a shadow casting technique.

Objects viewed through transparent sheets with residual non-parallelism and irregularity appear shifted and distorted. This distortion is measured in terms of angular and binocular deviation of an object viewed through the transparent sheet. The angular and binocular deviations introduced are particularly important in the context of aircraft windscreens and canopies as they can interfere with decision making of pilots especially while landing, leading to accidents. In this work, we have developed an instrument to measure both the angular and binocular deviations introduced by transparent sheets. This instrument is especially useful in the qualification of aircraft windscreens and canopies. It measures the deviation in the geometrical shadow cast by a periodic dot pattern trans-illuminated by the distorted light beam from the transparent test specimen compared to the reference pattern. Accurate quantification of the shift in the pattern is obtained by cross-correlating the reference shadow pattern with the specimen shadow pattern and measuring the location of the correlation peak. The developed instrument is handy to use and computes both angular and binocular deviation with an accuracy of less than ±0.1 mrad (≈0.036 mrad) and has an excellent repeatability with an error of less than 2%.

Accelerated gradient based diffuse optical tomographic image reconstruction.

PURPOSE: Fast reconstruction of interior optical parameter distribution using a new approach called Broyden-based model iterative image reconstruction (BMOBIIR) and adjoint Broyden-based MOBIIR (ABMOBIIR) of a tissue and a tissue mimicking phantom from boundary measurement data in diffuse optical tomography (DOT). METHODS: DOT is a nonlinear and ill-posed inverse problem. Newton-based MOBIIR algorithm, which is generally used, requires repeated evaluation of the Jacobian which consumes bulk of the computation time for reconstruction. In this study, we propose a Broyden approach-based accelerated scheme for Jacobian computation and it is combined with conjugate gradient scheme (CGS) for fast reconstruction. The method makes explicit use of secant and adjoint information that can be obtained from forward solution of the diffusion equation. This approach reduces the computational time many fold by approximating the system Jacobian successively through low-rank updates. RESULTS: Simulation studies have been carried out with single as well as multiple inhomogeneities. Algorithms are validated using an experimental study carried out on a pork tissue with fat acting as an inhomogeneity. The results obtained through the proposed BMOBIIR and ABMOBIIR approaches are compared with those of Newton-based MOBIIR algorithm. The mean squared error and execution time are used as metrics for comparing the results of reconstruction. CONCLUSIONS: We have shown through experimental and simulation studies that Broyden-based MOBIIR and adjoint Broyden-based methods are capable of reconstructing single as well as multiple inhomogeneities in tissue and a tissue-mimicking phantom. Broyden MOBIIR and adjoint Broyden MOBIIR methods are computationally simple and they result in much faster implementations because they avoid direct evaluation of Jacobian. The image reconstructions have been carried out with different initial values using Newton, Broyden, and adjoint Broyden approaches. These algorithms work well when the initial guess is close to the true solution. However, when initial guess is far away from true solution, Newton-based MOBIIR gives better reconstructed images. The proposed methods are found to be stable with noisy measurement data.

Assessment of ultrasound modulation of near infrared light on the quantification of scattering coefficient.

PURPOSE: To assess the effect of ultrasound modulation of near infrared (NIR) light on the quantification of scattering coefficient in tissue-mimicking biological phantoms. METHODS: A unique method to estimate the phase of the modulated NIR light making use of only time averaged intensity measurements using a charge coupled device camera is used in this investigation. These experimental measurements from tissue-mimicking biological phantoms are used to estimate the differential pathlength, in turn leading to estimation of optical scattering coefficient. A Monte-Carlo model based numerical estimation of phase in lieu of ultrasound modulation is performed to verify the experimental results. RESULTS: The results indicate that the ultrasound modulation of NIR light enhances the effective scattering coefficient. The observed effective scattering coefficient enhancement in tissue-mimicking viscoelastic phantoms increases with increasing ultrasound drive voltage. The same trend is noticed as the ultrasound modulation frequency approaches the natural vibration frequency of the phantom material. The contrast enhancement is less for the stiffer (larger storage modulus) tissue, mimicking tumor necrotic core, compared to the normal tissue. CONCLUSIONS: The ultrasound modulation of the insonified region leads to an increase in the effective number of scattering events experienced by NIR light, increasing the measured phase, causing the enhancement in the effective scattering coefficient. The ultrasound modulation of NIR light could provide better estimation of scattering coefficient. The observed local enhancement of the effective scattering coefficient, in the ultrasound focal region, is validated using both experimental measurements and Monte-Carlo simulations.

Single-resolution and multiresolution extended-Kalman-filter-based reconstruction approaches to optical refraction tomography.

The problem of reconstruction of a refractive-index distribution (RID) in optical refraction tomography (ORT) with optical path-length difference (OPD) data is solved using two adaptive-estimation-based extended-Kalman-filter (EKF) approaches. First, a basic single-resolution EKF (SR-EKF) is applied to a state variable model describing the tomographic process, to estimate the RID of an optically transparent refracting object from noisy OPD data. The initialization of the biases and covariances corresponding to the state and measurement noise is discussed. The state and measurement noise biases and covariances are adaptively estimated. An EKF is then applied to the wavelet-transformed state variable model to yield a wavelet-based multiresolution EKF (MR-EKF) solution approach. To numerically validate the adaptive EKF approaches, we evaluate them with benchmark studies of standard stationary cases, where comparative results with commonly used efficient deterministic approaches can be obtained. Detailed reconstruction studies for the SR-EKF and two versions of the MR-EKF (with Haar and Daubechies-4 wavelets) compare well with those obtained from a typically used variant of the (deterministic) algebraic reconstruction technique, the average correction per projection method, thus establishing the capability of the EKF for ORT. To the best of our knowledge, the present work contains unique reconstruction studies encompassing the use of EKF for ORT in single-resolution and multiresolution formulations, and also in the use of adaptive estimation of the EKF's noise covariances.

Convergence analysis of the Newton algorithm and a pseudo-time marching scheme for diffuse correlation tomography.

We propose a self-regularized pseudo-time marching scheme to solve the ill-posed, nonlinear inverse problem associated with diffuse propagation of coherent light in a tissuelike object. In particular, in the context of diffuse correlation tomography (DCT), we consider the recovery of mechanical property distributions from partial and noisy boundary measurements of light intensity autocorrelation. We prove the existence of a minimizer for the Newton algorithm after establishing the existence of weak solutions for the forward equation of light amplitude autocorrelation and its Fréchet derivative and adjoint. The asymptotic stability of the solution of the ordinary differential equation obtained through the introduction of the pseudo-time is also analyzed. We show that the asymptotic solution obtained through the pseudo-time marching converges to that optimal solution provided the Hessian of the forward equation is positive definite in the neighborhood of optimal solution. The superior noise tolerance and regularization-insensitive nature of pseudo-dynamic strategy are proved through numerical simulations in the context of both DCT and diffuse optical tomography.

Study of turbid media with light: recovery of mechanical and optical properties from boundary measurement of intensity autocorrelation of light.

We discuss the inverse problem associated with the propagation of the field autocorrelation of light through a highly scattering object like tissue. In the first part of the work, we reconstruct the optical absorption coefficient mu(a) and particle diffusion coefficient D(B) from simulated measurements which are integrals of a quantity computed from the measured intensity and intensity autocorrelation g(2)(tau) at the boundary. In the second part we recover the mean square displacement (MSD) distribution of particles in an inhomogeneous object from the sampled g(2)(tau) measured on the boundary. From the MSD, we compute the storage and loss moduli distributions in the object. We have devised computationally easy methods to construct the sensitivity matrices which are used in the iterative reconstruction algorithms for recovering these parameters from the measurements. The results of the reconstruction of mu(a), D(B), MSD and the viscoelastic parameters, which are presented, show reasonably good position and quantitative accuracy.

A pseudo-dynamic sub-optimal filter for elastography under static loading and measurements.

We propose a pseudo-dynamic form of a sub-optimal Kalman filter for elastography of plane-strain models of soft tissues under strictly static deformations and partial measurements. Since the tissue material is nearly incompressible and is thus prone to volumetric locking via standard displacement-based finite element formulations, we use a Cosserat point approach for deriving the static equilibrium equations. A pseudo-dynamical form of the equilibrium equations, with added noise and appropriate augmentation by the discretized shear modulus as additional states, is then adopted as the process equation such that its steady-state solution approaches the static response of the plane-strain model. A fictitious noise of small intensity is also added to the measurement equation and, following linearization of the process equation, a Kalman filter is applied to reconstruct the shear modulus profile. We present several numerical experiments, some of which also bring forth the relative advantages of the proposed approach over a deterministic reconstruction based on a quasi-Newton search.

Implementation of a phase array diffuse optical tomographic imager.

Diffuse optical tomography (DOT) using near-infrared (NIR) light is a promising tool for noninvasive imaging of deep tissue. This technique is capable of quantitative reconstructions of absorption coefficient inhomogeneities of tissue. The motivation for reconstructing the optical property variation is that it, and, in particular, the absorption coefficient variation, can be used to diagnose different metabolic and disease states of tissue. In DOT, like any other medical imaging modality, the aim is to produce a reconstruction with good spatial resolution and accuracy from noisy measurements. We study the performance of a phase array system for detection of optical inhomogeneities in tissue. The light transport through a tissue is diffusive in nature and can be modeled using diffusion equation if the optical parameters of the inhomogeneity are close to the optical properties of the background. The amplitude cancellation method that uses dual out-of-phase sources (phase array) can detect and locate small objects in turbid medium. The inverse problem is solved using model based iterative image reconstruction. Diffusion equation is solved using finite element method for providing the forward model for photon transport. The solution of the forward problem is used for computing the Jacobian and the simultaneous equation is solved using conjugate gradient search. The simulation studies have been carried out and the results show that a phase array system can resolve inhomogeneities with sizes of 5 mm when the absorption coefficient of the inhomogeneity is twice that of the background tissue. To validate this result, a prototype model for performing a dual-source system has been developed. Experiments are carried out by inserting an inhomogeneity of high optical absorption coefficient in an otherwise homogeneous phantom while keeping the scattering coefficient same. The high frequency (100 MHz) modulated dual out-of-phase laser source light is propagated through the phantom. The interference of these sources creates an amplitude null and a phase shift of 180 degrees along a plane between the two sources with a homogeneous object. A solid resin phantom with inhomogeneities simulating the tumor is used in our experiment. The amplitude and phase changes are found to be disturbed by the presence of the inhomogeneity in the object. The experimental data (amplitude and the phase measured at the detector) are used for reconstruction. The results show that the method is able to detect multiple inhomogeneities with sizes of 4 mm. The localization error for a 5 mm inhomogeneity is found to be approximately 1 mm.

Noise-tolerance analysis for detection and reconstruction of absorbing inhomogeneities with diffuse optical tomography using single- and phase-correlated dual-source schemes.

An iterative reconstruction procedure is used to invert intensity data from both single- and phase-correlated dual-source illuminations for absorption inhomogeneities. The Jacobian for the dual source is constructed by an algebraic addition of the Jacobians estimated for the two sources separately. By numerical simulations, it is shown that the dual-source scheme performs superior to the single-source system in regard to (i) noise tolerance in data and (ii) ability to reconstruct smaller and lower contrast objects. The quality of reconstructions from single-source data, as indicated by mean-square error at convergence, is markedly poorer compared to their dual-source counterpart, when noise in data was in excess of 2%. With fixed contrast and decreasing inhomogeneity diameter, our simulations showed that, for diameters below 7 mm, the dual-source scheme has a higher percentage contrast recovery compared to the single-source scheme. Similarly, the dual-source scheme reconstructs to a higher percentage contrast recovery from lower contrast inhomogeneity, in comparison to the single-source scheme.

Application of ultrasound-tagged photons for measurement of amplitude of vibration of tissue caused by ultrasound: theory, simulation, and experiments.

We investigate the modulation of an optical field caused by its interaction with an ultrasound beam in a tissue mimicking phantom. This modulation appears as a modulation in the intensity autocorrelation, which is measured by a photon counting correlator. The factors contributing to the modulation are: 1. amplitude of vibration of the particles of the tissue, 2. refractive index modulation, and 3. absorption coefficient in the region of the tissue intercepted by the ultrasound beam and light. We show in this work that a significant part of the contribution to this modulation comes from displacement of the tissue particles, which in turn is governed by the elastic properties of the tissue. We establish, both through simulations and experiments using an optical elastography phantom, the effects of the elasticity and absorption coefficient variations on the modulation of intensity autocorrelation. In the case where there is no absorption coefficient variation, we suggest that the depth of modulation can be calibrated to measure the displacement of tissue particles that, in turn, can be used to measure the tissue elasticity.

Diffuse optical tomography through solving a system of quadratic equations: theory and simulations.

This paper discusses the iterative solution of the nonlinear problem of optical tomography. In the established forward model-based iterative image reconstruction (MOBIIR) method a linear perturbation equation containing the first derivative of the forward operator is solved to obtain the update vector for the optical properties. In MOBIIR, the perturbation equation is updated by recomputing the first derivative after each update of the optical properties. In the method presented here a nonlinear perturbation equation, containing terms up to the second derivative, is used to iteratively solve for the optical property updates. Through this modification, reconstructions with reasonable contrast recovery and accuracy are obtained without the need for updating the perturbation equation and therefore eliminating the outer iteration of the usual MOBIIR algorithm. To improve the performance of the algorithm the outer iteration is reintroduced in which the perturbation equation is recomputed without re-estimating the derivatives and with only updated computed data. The system of quadratic equations is solved using either a modified conjugate gradient descent scheme or a two-step linearized predictor-corrector scheme. A quick method employing the adjoint of the forward operator is used to estimate the derivatives. By solving the nonlinear perturbation equation it is shown that the iterative scheme is able to recover large contrast variations in absorption coefficient with improved noise tolerance in data. This ability has not been possible so far with linear algorithms. This is demonstrated by presenting results of numerical simulations from objects with inhomogeneous inclusions in absorption coefficient with different contrasts and shapes.

Design, fabrication, and characterization of a tissue-equivalent phantom for optical elastography.

We suitably adapt the design of a tissue-equivalent phantom used for photoacoustic imaging to construct phantoms for optical elastography. The elastography phantom we consider should have optical properties such as scattering coefficient, scattering anisotropy factor, and refractive index; mechanical properties such as storage and loss modulus; and acoustic properties such as ultrasound velocity, attenuation coefficient, and acoustic impedance to match healthy and diseased tissues. The phantom is made of poly (vinyl alcohol) (PVA) and its mechanical, optical, and acoustic properties are tailored by physical cross-linking effected through subjecting a suitable mix of PVA stock and water to a number of freeze-thaw cycles and by varying the degree of hydrolysis in the PVA stock. The optical, mechanical, and acoustic properties of the samples prepared are measured by employing different techniques. The measured variations in the values of optical scattering coefficient, scattering anisotropy factor, and refractive index and storage modulus are found to be comparable to those in normal and diseased breast tissues. The acoustic properties such as sound speed, acoustic attenuation coefficient, and density are found to be close to the average values reported in the literature for normal breast tissue.

Diffuse optical tomography using intensity measurements and the a priori acquired regions of interest: theory and simulations.

Light transmission data collected around an object show large variation with source-detector separation owing to the presence of single or multiple inhomogeneous regions in the object. This variation in the measured intensity is made use of to reconstruct regions of the inhomogeneous inclusions. In addition, it is possible to select a set of data from the above which is most likely least affected by the presence of the inhomogeneity, and estimate reasonably accurately the background optical properties from it. The reconstructed region is found to always contain the inhomogeneity and is of size approximately 140% by area of the inhomogeneity. With the regions to be reconstructed a priori known, a model-based iterative reconstruction procedure for reconstructing the optical properties of the region converged five times faster than without such information. It is also shown that whereas for the full object, a view-based propagation-backpropagation reconstruction procedure failed to converge, owing to large underdeterminacy of the problem, a smaller problem attempting to reconstruct a priori specified regions of interest converged and did so faster than a non-view-based approach for similar objects. Reconstruction results are presented from simulated transmitted intensity data from the following objects with regions of inhomogeneity in both absorption and scattering: (i) single centrally located inhomogeneity, (ii) two off-centred inhomogeneous regions of equal size and contrast (iii) two off-centred inhomogeneous regions of unequal size and equal contrast and (iv) two off-centred inhomogeneous regions of unequal size and contrast. Whereas the model-based iterative image reconstruction procedure gave good convergence in the first, second and third cases, in the fourth case the reconstructions failed to recover the exact numerical value of the optical properties in the higher contrast region.

Reconstruction of optical properties of low-scattering tissue using derivative estimated through perturbation Monte-Carlo method.

An iterative method for the reconstruction of optical properties of a low-scattering object, which uses a Monte-Carlo-based forward model, is developed. A quick way to construct and update the Jacobian needed to reconstruct a discretized object, based on the perturbation Monte-Carlo (PMC) approach, is demonstrated. The projection data is handled either one view at a time, using a propagation-backpropagation (PBP) strategy where the dimension of the inverse problem and consequently the computation time are smaller, or, when this approach failed, using all the views simultaneously with a full dataset. The main observations and results are as follows. 1. Whereas the PMC gives an accurate and quick method for constructing the Jacobian the same, when adapted to update the computed projection data, the data are not accurate enough for use in the iterative reconstruction procedure leading to convergence. 2. The a priori assumption of the location of inhomogeneities in the object reduces the dimension of the problem, leading to faster convergence in all the cases considered, such as an object with multiple inhomogeneities and data handled one view at a time (i.e., the PBP approach). 3. On the other hand, without a priori knowledge of the location of inhomogeneities, the problem was too ill posed for the PBP approach to converge to meaningful reconstructions when both absorption and scattering coefficients are considered as unknowns. Finally, to bring out the effectiveness of this method for reconstructing low-scattering objects, we apply a diffusion equation-based algorithm on a dataset from one of the low-scattering objects and show that it fails to reconstruct object inhomogeneities.

Optical tomographic microscope for quantitative imaging of phase objects.

We describe a tomographic microscope, for imaging phase objects, that makes use of the transport-of-intensity equation to estimate the phase of the transmitted light through the object. The wave-front data from optical fibers are reconstructed with an algorithm that incorporates correction for the ray bending. The reconstructed refractive-index cross sections of the fibers are found to be in agreement with the available values specified in the catalogs.

Enhancement of far-field spatial coherence of a partially coherent light source by source masking.

The role of a sharp autocorrelation phase mask, called the bleached uniformly redundant array, for improving the spatial coherence in the far-field of partially coherent light sources is studied. It is shown both theoretically and experimentally that the input source correlation plays an important role in determining the amount of enhancement introduced by the phase mask.

Coherent optical processing using noncoherent light after source masking.

Coherent optical processing starting with spatially noncoherent illumination is described. Good spatial coherence is introduced in the far field by modulating a noncoherent source when masks with sharp autocorrelation are used. The far-field mutual coherence function of light is measured and it is seen that, for the masks and the source size used here, we get a fairly large area over which the mutual coherence function is high and flat. We demonstrate traditional coherent processing operations such as Fourier transformation and image deblurring when coherent light that is produced in the above fashion is used. A coherence-redundancy merit function is defined for this type of processing system. It is experimentally demonstrated that the processing system introduced here has superior blemish tolerance compared with a traditional processor that uses coherent illumination.

Generation of optimum random pinhole arrays for multiple imaging.

Through an analysis using the transfer function of a pinhole camera, the multiple imaging characteristics of photographic diffusers described by Grover and Tremblay [Appl. Opt. 21, 4500(1982)] is studied. It is found that only one pinhole diameter satisfies the optimum imaging condition for best contrast transfer at any desired spatial frequency. A simple method of generating random pinhole arrays with a controlled pinhole diameter is described. These pinhole arrays are later used to generate high frequency sinusoidal gratings from a coarse grid. The contrast in the final gratings is found to be reasonably high.

Multiple imaging and multichannel optical processing with split lenses.

The use of split lenses for multiple imaging and multichannel optical processing is demonstrated. Conditions are obtained for nonoverlapping of multipled images and avoiding crosstalk in the multichannel processing. Almost uniform intensity across the multipled images is an advantage here, while the low f/No. of the split lens segments puts a limit in the resolution in image processing. Experimental results of multiple imaging and of a few multichannel processing are presented.