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. (C) 2011 American Association of Physicists in Medicine. DOI: 10.1118/1.3531572]

A pseudo-time EnKF incorporating shape based reconstruction for diffuse optical tomography

Purpose: The authors aim at developing a pseudo-time, sub-optimal stochastic filtering approach based on a derivative free variant of the ensemble Kalman filter (EnKF) for solving the inverse problem of diffuse optical tomography (DOT) while making use of a shape based reconstruction strategy that enables representing a cross section of an inhomogeneous tumor boundary by a general closed curve. Methods: The optical parameter fields to be recovered are approximated via an expansion based on the circular harmonics (CH) (Fourier basis functions) and the EnKF is used to recover the coefficients in the expansion with both simulated and experimentally obtained photon fluence data on phantoms with inhomogeneous inclusions. The process and measurement equations in the pseudo-dynamic EnKF (PD-EnKF) presently yield a parsimonious representation of the filter variables, which consist of only the Fourier coefficients and the constant scalar parameter value within the inclusion. Using fictitious, low-intensity Wiener noise processes in suitably constructed ``measurement'' equations, the filter variables are treated as pseudo-stochastic processes so that their recovery within a stochastic filtering framework is made possible. Results: In our numerical simulations, we have considered both elliptical inclusions (two inhomogeneities) and those with more complex shapes (such as an annular ring and a dumbbell) in 2-D objects which are cross-sections of a cylinder with background absorption and (reduced) scattering coefficient chosen as mu(b)(a)=0.01mm(-1) and mu('b)(s)=1.0mm(-1), respectively. We also assume mu(a) = 0.02 mm(-1) within the inhomogeneity (for the single inhomogeneity case) and mu(a) = 0.02 and 0.03 mm(-1) (for the two inhomogeneities case). The reconstruction results by the PD-EnKF are shown to be consistently superior to those through a deterministic and explicitly regularized Gauss-Newton algorithm. We have also estimated the unknown mu(a) from experimentally gathered fluence data and verified the reconstruction by matching the experimental data with the computed one. Conclusions: The PD-EnKF, which exhibits little sensitivity against variations in the fictitiously introduced noise processes, is also proven to be accurate and robust in recovering a spatial map of the absorption coefficient from DOT data. With the help of shape based representation of the inhomogeneities and an appropriate scaling of the CH expansion coefficients representing the boundary, we have been able to recover inhomogeneities representative of the shape of malignancies in medical diagnostic imaging. (C) 2012 American Association of Physicists in Medicine. [DOI: 10.1118/1.3679855]

Linearization and Reconstruction of Non-linear Diffuse Optical Tomographic Image

Diffuse optical tomography (DOT) is one of the ways to probe highly scattering media such as tissue using low-energy near infra-red light (NIR) to reconstruct a map of the optical property distribution. The interaction of the photons in biological tissue is a non-linear process and the phton transport through the tissue is modelled using diffusion theory. The inversion problem is often solved through iterative methods based on nonlinear optimization for the minimization of a data-model misfit function. The solution of the non-linear problem can be improved by modeling and optimizing the cost functional. The cost functional is f(x) = x(T)Ax - b(T)x + c and after minimization, the cost functional reduces to Ax = b. The spatial distribution of optical parameter can be obtained by solving the above equation iteratively for x. As the problem is non-linear, ill-posed and ill-conditioned, there will be an error or correction term for x at each iteration. A linearization strategy is proposed for the solution of the nonlinear ill-posed inverse problem by linear combination of system matrix and error in solution. By propagating the error (e) information (obtained from previous iteration) to the minimization function f(x), we can rewrite the minimization function as f(x; e) = (x + e)(T) A(x + e) - b(T)(x + e) + c. The revised cost functional is f(x; e) = f(x) + e(T)Ae. The self guided spatial weighted prior (e(T)Ae) error (e, error in estimating x) information along the principal nodes facilitates a well resolved dominant solution over the region of interest. The local minimization reduces the spreading of inclusion and removes the side lobes, thereby improving the contrast, localization and resolution of reconstructed image which has not been possible with conventional linear and regularization algorithm.

Structural and Optical Modification in Bi doped As40S60 Thin Films.

The influence of substitution of Bi atom instead of S atoms on the structural and optical properties of thin films of As40S60 are reported. The density is found to be increased with the addition Bi heavy metal into As2S3. The amorphous to polycrystalline structure of the bulk sample is observed for Bi more than 7%. The glass transition temperature is found to be decreased with addition of Bi. The absorption edge shifts to shorter wavelength, thereby decreasing optical band gap of BixAs(40)S(60-x) (x= 0,2 and 4% here) film. The optical parameter change is discussed from the stand point of chemical bonds formed in the films and related to the defect states produced due to incorporation of Bi atoms in place of chalcogenide S atoms.

Region-of-interest diffuse optical tomography system

Diffuse optical tomography (DOT) using near-infrared light is a promising tool for non-invasive imaging of deep tissue. This technique is capable of quantitative reconstruction of absorption (mu(a)) and scattering coefficient (mu(s)) inhomogeneities in the tissue. The rationale for reconstructing the optical property map is that the absorption coefficient variation provides diagnostic information about metabolic and disease states of the tissue. The aim of DOT is to reconstruct the internal tissue cross section with good spatial resolution and contrast from noisy measurements non-invasively. We develop a region-of-interest scanning system based on DOT principles. Modulated light is injected into the phantom/tissue through one of the four light emitting diode sources. The light traversing through the tissue gets partially absorbed and scattered multiple times. The intensity and phase of the exiting light are measured using a set of photodetectors. The light transport through a tissue is diffusive in nature and is modeled using radiative transfer equation. However, a simplified model based on diffusion equation (DE) can be used if the system satisfies following conditions: (a) the optical parameter of the inhomogeneity is close to the optical property of the background, and (b) mu(s) of the medium is much greater than mu(a) (mu(s) >> mu(a)). The light transport through a highly scattering tissue satisfies both of these conditions. A discrete version of DE based on finite element method is used for solving the inverse problem. The depth of probing light inside the tissue depends on the wavelength of light, absorption, and scattering coefficients of the medium and the separation between the source and detector locations. Extensive simulation studies have been carried out and the results are validated using two sets of experimental measurements. The utility of the system can be further improved by using multiple wavelength light sources. In such a scheme, the spectroscopic variation of absorption coefficient in the tissue can be used to arrive at the oxygenation changes in the tissue. (C) 2016 AIP Publishing LLC.

Minimal residual method provides optimal regularization parameter for diffuse optical tomography

The inverse problem in the diffuse optical tomography is known to be nonlinear, ill-posed, and sometimes under-determined, requiring regularization to obtain meaningful results, with Tikhonov-type regularization being the most popular one. The choice of this regularization parameter dictates the reconstructed optical image quality and is typically chosen empirically or based on prior experience. An automated method for optimal selection of regularization parameter that is based on regularized minimal residual method (MRM) is proposed and is compared with the traditional generalized cross-validation method. The results obtained using numerical and gelatin phantom data indicate that the MRM-based method is capable of providing the optimal regularization parameter. (C) 2012 Society of Photo-Optical Instrumentation Engineers (SPIE). DOI: 10.1117/1.JBO.17.10.106015]

A LSQR-type method provides a computationally efficient automated optimal choice of regularization parameter in diffuse optical tomography

Purpose: Developing a computationally efficient automated method for the optimal choice of regularization parameter in diffuse optical tomography. Methods: The least-squares QR (LSQR)-type method that uses Lanczos bidiagonalization is known to be computationally efficient in performing the reconstruction procedure in diffuse optical tomography. The same is effectively deployed via an optimization procedure that uses the simplex method to find the optimal regularization parameter. The proposed LSQR-type method is compared with the traditional methods such as L-curve, generalized cross-validation (GCV), and recently proposed minimal residual method (MRM)-based choice of regularization parameter using numerical and experimental phantom data. Results: The results indicate that the proposed LSQR-type and MRM-based methods performance in terms of reconstructed image quality is similar and superior compared to L-curve and GCV-based methods. The proposed method computational complexity is at least five times lower compared to MRM-based method, making it an optimal technique. Conclusions: The LSQR-type method was able to overcome the inherent limitation of computationally expensive nature of MRM-based automated way finding the optimal regularization parameter in diffuse optical tomographic imaging, making this method more suitable to be deployed in real-time. (C) 2013 American Association of Physicists in Medicine. http://dx.doi.org/10.1118/1.4792459]

Optimum Parameter Selection in Sparse Reconstruction of Frequency-Domain Optical-Coherence Tomography Signals

For a multilayered specimen, the back-scattered signal in frequency-domain optical-coherence tomography (FDOCT) is expressible as a sum of cosines, each corresponding to a change of refractive index in the specimen. Each of the cosines represent a peak in the reconstructed tomogram. We consider a truncated cosine series representation of the signal, with the constraint that the coefficients in the basis expansion be sparse. An l(2) (sum of squared errors) data error is considered with an l(1) (summation of absolute values) constraint on the coefficients. The optimization problem is solved using Weiszfeld's iteratively reweighted least squares (IRLS) algorithm. On real FDOCT data, improved results are obtained over the standard reconstruction technique with lower levels of background measurement noise and artifacts due to a strong l(1) penalty. The previous sparse tomogram reconstruction techniques in the literature proposed collecting sparse samples, necessitating a change in the data capturing process conventionally used in FDOCT. The IRLS-based method proposed in this paper does not suffer from this drawback.

Defect-selective photothermal imaging: a sensitive tool for damage mapping in optical coatings

A defect-selective photothermal imaging system for the diagnostics of optical coatings is demonstrated. The instrument has been optimized for pump and probe parameters, detector performance, and signal processing algorithm. The imager is capable of mapping purely optical or thermal defects efficiently in coatings of low damage threshold and low absorbance. Detailed mapping of minor inhomogeneities at low pump power has been achieved through the simultaneous action of a low-noise fiber optic photothermal beam defection sensor and a common-mode-rejection demodulation (CMRD) technique. The linearity and sensitivity of the sensor have been examined theoretically and experimentally, and the signal to noise ratio improvement factor is found to be about 110 compared to a conventional bicell photodiode. The scanner is so designed that mapping of static or shock sensitive samples is possible. In the case of a sample with absolute absorptance of 3.8 x 10(-4), a change in absorptance of about 0.005 x 10(-4) has been detected without ambiguity, ensuring a contrast parameter of 760. This is about 1085% improvement over the conventional approach containing a bicell photodiode, at the same pump power. The merits of the system have been demonstrated by mapping two intentionally created damage sites in a MgF2 coating on fused silica at different excitation powers. Amplitude and phase maps were recorded for thermally thin and thick cases, and the results are compared to demonstrate a case which, in conventional imaging, would lead to a deceptive conclusion regarding the type and location of the damage. Also, a residual damage profile created by long term irradiation with high pump power density has been depicted.

In situ pump probe optical absorption studies on Sb/As2S3 nanomultilayered film

We report the in situ optical transmission change in the complete visible region of the electromagnetic spectrum to asses the kinetics of photo induced interdiffusion in Sb/As2S3 nanomultilayered film. The interdiffusion of Sb into As2S3 matrix results in the formation of Sb-As2S3 ternary solid solutions which is explained by the change in optical band gap, absorption coefficients and Tauc parameter (B-1/2) with evolution of time. The wavelength dependence of the time constants provides a better description of photo induced effects. The time evolution of the absorption coefficients and optical band gap are significantly faster in this process. (C) 2009 Elsevier B.V. All rights reserved.

Structural, dielectric, impedance and optical properties of CaBi2B2O7 glasses and glass-nanocrystal composites

Optically clear glasses were fabricated by quenching the melt of CaCO3-Bi2O3-B2O3 (in equimolecular ratio). The amorphous and glassy characteristics of the as-quenched samples were confirmed via the X-ray powder diffraction (XRD) and differential scanning calorimetric (DSC) studies These glasses were found to. have high thermal stability parameter (S). The optical transmission studies carried out in the 200-2500 nm wavelength range confirmed both the as-quenched and heat-treated samples to be transparent between 400 nm and 2500 nm. The glass-plates that were heat-treated just above the glass transition temperature (723 K) for 6 h retained approximate to 60% transparency despite having nano-crystallites (approximate to 50-100 nm) of CaBi2B2O7 (CBBO) as confirmed by both the XRD and transmission electron microscopy (TEM) studies. The dielectric properties and impedance characteristics of the as-quenched and heat-treated (723 K/6 h) samples were studied as a function of frequency at different temperatures. Cole-Cole equation was employed to rationalize the impedance data.

Simplified theory of optical nonlinearities in spin-polarized semiconductors

The spin degree of freedom is largely disregarded in existing theories of the density-dependent optical properties of an interacting electron-hole plasma in quasiequilibrium. Here, we extended the pair equation, which is applicable to a bulk semiconductor at elevated temperatures, to calculate optical nonlinearities due to a spin-polarized plasma. We obtained agreement with recent circular dichroism data in laser-excited GaAs by using the plasma density alone as the fitting parameter. The simplicity of our theory, based on the analytical pair-equation formula, makes it ideal for conveniently modelling absorption in a carrier spin-polarized semiconductor.

On self-interacting oligoribonucleotides. I. Absorption and optical rotatory dispersion of 2′-5′- and 3′-5′-oligoguanylic acids

A series of 2′-5′-oligoguanylic acids are prepared by reacting G(cyclic)p with takadiastase T1 ribonuclease and separating the products chromatographically. The 3′-5′-oligoguanylic acids are obtained by separating the products of alkaline degradation of 3′-5′-poly(G). The optical rotatory dispersion and hypochromism of both 2′-5′- and 3′-5′-oligoguanylic acids are studied at two different pH. The optical rotatory dispersion spectrum of 2′-5′-GpG is significantly different from that of 3′-5′-GpG. The magnitude of rotation of the long-wavelength peak of 2′-5′-GpG is larger than that of 3′-5′-GpG. This finding contradicts the explanation that the extra stability and more intense circular dichroism band of other 3′-5′-dinucleoside monophosphates is due to H-bond formation between 2′-OH and either the base or the phosphate oxygen. The end phosphate group has a marked effect on the spectrum of GpG between 230 and 250 mμ. In addition the optical rotatory dispersion spectra of 2′-5′ exhibit strong pH, temperature, and solvent dependence between 230 and 250 mμ. ΔH and AS for order ⇌ disorder transition is estimated to be 9.7 kcal/mole and 35.2 eu, respectively. The optical rotatory dispersion spectra of guanine-rich oligoribonucleotides, GpGpC, GpGpU, GpGpGpC, and GpGpGpU are compared to the calculated optical rotatory dispersion from the semiempirical expression of Cantor and Tinoco, using measured optical rotatory dispersion of dimers. Contrary to previous studies, agreement is found not at all satisfactory. However, optical rotatory dispersion of 3′-5′-GpGpGpC and GpGpGpU can be estimated from the semiempirical expression, if a next-nearest interaction parameter is introduced empirically. Such interaction parameter can be calculated from the measured properties of trinucleotide sequences like GpGpG, GpGpC, and GpGpU, assuming that only the nearest-neighbor interaction is important. The optical rotatory dispersion of single-stranded poly(G) is also predicted. The importance of syn-anti equilibrium and next-nearest-neighbor interaction in oligoguanylic acids is suggested as a probable explanation.

Photo induced optical changes in Sb/As2S3 multilayered film and (As2S3)(0.93)Sb-0.07 film of equal thickness

The increase in optical band gap (photo bleaching) due to light illumination was studied at room temperature as well as at low (4.2 K) temperature for Sb/As2S3 multilayered film of 640 nm thickness by Fourier Transform Infrared Technique. The interdiffusion of Sb into As2S3 matrix results the formation of Sb-As2S3 ternary solid solutions which is explained by the change in optical band gap (E-g), absorption coefficient (alpha), Tauc parameter (B-1/2), Urbach edge (E-e). At the same time, photo darkening phenomena was observed in (As2S3)(0.93)Sb-0.07 film of same thickness both at low and room temperatures. From our X-ray Photoelectron Spectroscopy measurements,we are able to show that some of the As-As, S-S and Sb-Sb bonds are converted into As-S and S-Sb bonds in case of multilayers. We found that the energetically favoured heteropolar bond formation take place by a phonon-assisted mechanism using the lone pair pi electrons of S-2(0). But in case of (As2S3)(0.93)Sb-0.02 film, the homopolar bonds are playing a major role. (C) 2010 Elsevier B.V. All rights reserved.

Temperature dependent optical constants of amorphous Ge2Sb2Te5 thin films

This study focuses on the temperature dependent optical band gap changes in the amorphous Ge2Sb2Te5 (GST) films. The behavior of the amorphous GST thin films at low temperatures has been studied. The band gap increment of around 0.2 eV is observed at low temperature (4.2 K) compared to room temperature (300 K). The band gap changes associated with the temperature are completely reversible. The other optical parameters like Urbach energy and Tauc parameter (B-1/2) are studied for different temperatures and discussed. The observed changes in optical band gap (E-g) are fitting to Fan's one phonon approximation. Phonon energy ((h) over bar omega) corresponding to a frequency of 3.59 THz is derived from Fan's approximation, which is close to the reported value of 3.66 THz. (C) 2010 Elsevier B.V. All rights reserved.