Experimental demonstration of an analytic method for image reconstruction in optical diffusion tomography with large data sets

We report the first experimental test of an analytic image reconstruction algorithm for optical tomography with large data sets. Using a continuous-wave optical tomography system with 10(8) source-detector pairs, we demonstrate the reconstruction of an absorption image of a phantom consisting of a highly scattering medium containing absorbing inhomogeneities.     

彩色漫射物体的压缩全息层析成像

漫射物体的压缩全息利用其非相干散射密度函数在统计意义上满足稀疏先验这一假设,可以从多幅散斑图案实现漫射物体的层析重建,避免了散斑和不同平面的散焦图像之间的串扰。将单波长照明条件拓展到红、绿、蓝三色波长,提出一种新的适用于彩色漫射物体的压缩全息层析成像方法,搭建多波长照明条件下漫射物体的层析成像模型,并通过解压缩推理方法有效地分离不同平面的三维非相干密度函数。数值模拟实验表明,所提方法可以在多幅二维彩色散斑图案中成功进行彩色层析漫射物体的压缩重建,有效地抑制了散斑效应以及不同平面的散焦图像之间的相互串扰。     

Three-element phased-array approach to diffuse optical imaging based on postprocessing of continuous-wave data

We present a multielement phased-array approach to diffuse optical imaging based on postprocessing of continuous-wave data for the improvement of spatial resolution. In particular, we present a theoretical and experimental analysis of the performance of a three-element source array in the study of an optically turbid medium with two embedded cylindrical inclusions. We find that the proposed phased-array approach is able to resolve two cylinders with side-to-side separation of 10 mm that are not resolved by the intensity associated with a single light source.     

Acousto-optic-assisted diffuse optical tomography

We introduce and experimentally demonstrate acousto-optic-assisted diffuse optical tomography (DOT) using a holography-based acousto-optic setup. The method is based on probing a scattering medium with a localized acoustical modulation of the phase of the scattered light. The optical properties of the scattering medium are recovered with ultrasound-limited resolution by applying DOT reconstruction methods on a set of the measured intensities of light, modulated at different locations throughout the medium.     

Image-guided optical spectroscopy provides molecular-specific information in vivo: MRI-guided spectroscopy of breast cancer hemoglobin, water, and scatterer size

A multimodality instrument that integrated optical or near-infrared spectroscopy into a magnetic resonance imaging (MRI) breast coil was used to perform a pilot study of image-guided spectroscopy on cancerous breast tissue. These results are believed to be the first multiwavelength spectroscopic images of breast cancer using MRI-guided constraints, and they show the cancer tumor to have high hemoglobin and water values, decreased oxygen saturation, and increased subcellular granularity. The use of frequency-domain diffuse tomography methods at many wavelengths provides the spectroscopy required for recovering maps of absorbers and scattering spectra, but the integration with MRI allows these data to be recovered on an image field that preserves high resolution and fuses the two data sets together. Integration of molecular spectroscopy into standard clinical MRI can be achieved with this approach to spectral tomography.     

In vivo near-infrared spectral detection of pressure-induced changes in breast tissue

A diffuse near-infrared tomography system was used to measure dynamic changes in the absolute optical properties of the human breast that were induced through pressure applied to the tissue surface. Results from five subjects show that absorption and scattering coefficients changed measurably when pressure was increased and that these relative changes correlated with the subjects' body-mass index, indicating that the effect depends on tissue composition. Fitting the absolute absorption and scattering coefficients at six wavelengths to the molar absorption spectra of the three predominant chromophores revealed that both the average total hemoglobin and oxygen saturation increased by 10%, while water concentration decreased by more than 12%. These changes indicate that the pressure-induced variation is likely due to water displacement and vascular volume increase in the region being imaged, for mild application of pressure to the breast. These results suggest that the pressure applied during optical measurements of tissue may alter the tissue physiology, and care should be taken to factor this effect into the design of optical medical instrumentation. In addition, the technique provides a unique approach to measuring tissue elastic changes in vivo in the female breast and may offer a new method for dynamic contrast imaging based on elasto-optical measurements.     

Depth-adaptive regularized reconstruction for reflection diffuse optical tomography

We study reflection diffuse optical tomography using two-dimensional (2D) continuous-wave source-detector arrays on the surface of semi-infinite medium, aiming at imaging the perfusion and the hemoglobin oxygen saturation variation of human cerebral cortex with brain activation. We had previously formulated the inverse problem with Moore-Penrose inversion. When we use simple regularization in this inverse problem, the reconstruction sensitivity decreases markedly with the depth so that the signal in the deep range may be masked by an unwanted signal in the shallow range. In this paper, we propose a depth-adaptive regularized reconstruction, in which we assign a smaller regularization parameter with the depth. We demonstrate improvement of the three-dimensional (3D) reconstruction uniformity using the proposed scheme.     

离体正常乳腺组织350～850 nm波段光谱特性

采用带有积分球附件的紫外/可见/近红外分光光度计测量了离体正常乳腺组织在350～850 nm光谱范围的反射率和透射率,运用反向倍加法得到了离体正常乳腺组织在相应光谱范围的光学参数,分析了正常乳腺组织的光学穿透深度随波长的变化情况。实验结果表明：350～850 nm波段正常乳腺组织的约化散射系数μ′_s大于吸收系数μ_a。μ′_s随着波长的增加而减小,即从350 nm波长值为9.731 mm-1～850 nm波长值为1.476 mm-1。μ_a从350 nm波长值为0.798 mm-1～850 nm波长值为0.102 mm-1,410 nm波长处存在一个吸收峰,其值为0.506 mm-1。光学穿透深度随着波长的增加而增大,从350 nm波长值为0.199 mm-1～850 nm波长值为1.439 mm。基于反向倍加法计算获得乳腺组织的光学参数,采用Monte Carlo模拟得到其相应光谱范围的反射率和透射率,并与实际测量值进行比较,二者的一致性较好。实验结果为乳腺组织的光活检及其光学治疗提供重要参考。     

Analysis of light propagation in highly scattering media by path-length-assigned Monte Carlo simulations

Numerical analysis of optical propagation in highly scattering media is investigated when light is normally incident to the surface and re-emerges backward from the same point. This situation corresponds to practical light scattering setups, such as in optical coherence tomography. The simulation uses the path-length-assigned Monte Carlo method based on an ellipsoidal algorithm. The spatial distribution of the scattered light is determined and the dependence of its width and penetration depth on the path-length is found. The backscattered light is classified into three types, in which ballistic, snake, and diffuse photons are dominant.     

Demonstration of a simple entangling optical gate and its use in bell-state analysis

We demonstrate a new architecture for an optical entangling gate that is significantly simpler than previous realizations, using partially polarizing beam splitters so that only a single optical mode-matching condition is required. We demonstrate operation of a controlled-z gate in both continuous-wave and pulsed regimes of operation, fully characterizing it in each case using quantum process tomography. We also demonstrate a fully resolving, nondeterministic optical Bell-state analyzer based on this controlled-z gate. This new architecture is ideally suited to guided optics implementations of optical gates.     

面向近红外脑功能成像的光学自导引漫射光断层成像方法

为降低近红外脑功能漫射光断层成像(DOT)固有的逆问题病态性,并避免多模态方法的图像配准等问题,提出了基于光学自导引提供先验功能信息的脑功能DOT方法(OT-DOT),并发展了图像重构方法.模拟验证表明:上皮厚度(TLT)已知时,OT-DOT获得的重构量化度(QR)约为传统DOT的4.2倍;当TLT的估计误差小于±10%时,OT-DOT重构的QR值可达92%以上,远远优于传统DOT;噪声鲁棒性测试表明,OT-DOT与传统DOT的噪声鲁棒性相近.利用连续光DOT测量系统的仿体实验重构结果表明,所发展的OT-DOT算法获得的重构结果优于传统DOT算法.     

Optical bistability by nonlinear reflection in a polymeric blend

We have analysed the propagation of two continuous-wave Ar⁺-laser beams travelling in opposite directions in a two-component nonlinear heterogeneous medium (the polymeric blend PMMA-EVA), in a mirror retroreflection geometry. Bistable reflection associated with the laser-induced suppression of scattering is observed. This phenomenon has a thermal origin. Based on light- and heat-transport theory inside the material, a simplified explanation of the experimental evidence is presented, showing the conditions for the onset of optical bistability. Qualitative agreement with our experimental results is found.     

On equivalent diffuse conditions for an internal layer in multislab atmospheric radiation

This paper reports on new discrete ordinates conditions for efficiently solving a set of multislab atmospheric radiation problems characterized by an optically stationary internal layer, i.e. an internal layer whose optical (absorption/scattering) properties and optical thickness do not change from one problem to another in the set. The discrete ordinates conditions reported here are founded in a recently developed spectral nodal method for solving multislab atmospheric radiation problems with anisotropic scattering. We suitably use the optically discretized equations of our recently developed spectral nodal method to derive discrete ordinates diffuse conditions, which model the response—diffuse radiation leaving the layer—of an internal layer to an anisotropic inner source and to diffuse radiation that is incident upon the layer at top and bottom. These conditions can be used to replace an optically stationary internal layer in multislab atmospheric radiation computations, while saving computer resources and without degrading the numerical results.     

高光谱分辨率激光雷达鉴频器的设计与分析

提出了利用Fabry-Perot干涉仪的反射场实现高光谱分辨率激光雷达精细探测大气光学参量的新方法和思路.设计了高光谱分辨率的分光系统,并分析了干涉仪反射场的光谱透过特征曲线.结合高光谱激光雷达探测信号特征,讨论分析了谱分离比和瑞利信号透过率随反射率和腔长的变化曲线,同时结合误差传递公式,建立了仿真分析模型,讨论了回波光束发散角和入射角变化对激光雷达探测结果的影响.结果表明,所提出的Fabry-Perot干涉仪反射场可以实现高光谱分辨率激光雷达探测系统的精细分光,同时探测结果误差随回波光束发散角变化不敏感,控制发散角在10 mrad以内,入射角在1.5 mrad以内时,可以实现气溶胶光学参数廓线的高精度探测.     

Morphological characterization of cells in concentrated suspensions using multispectral diffuse optical tomography

Based on a non-spherical model of particle scattering, we investigate the capabilities and limitations of a T-matrix based inverse algorithm to morphologically characterize cells in concentrated suspensions. Here the cells are modeled as randomly orientated spheroidal particles with homogenous dielectric properties and suspended in turbid media. The inverse algorithm retrieves the geometrical parameters and the concentration of cells simultaneously by inverting the reduced scattering coefficient spectra obtained from multispectral diffuse optical tomography (MS-DOT). Both round and spheroidal cells are tested and the role of multiple and higher order scattering of particles on the performance of the algorithm is evaluated using different concentrations of cells.     

Spectral broadening and stimulated Raman conversion in a continuous-wave optical parametric oscillator

We have confirmed that single-frequency oscillation of a continuous-wave singly resonant optical parametric oscillator is limited to operation below a critical value of the pumping ratio, as predicted by early theoretical treatments. We also report different regimes of spectral broadening as well as stimulated Raman conversion of the signal wave above this critical pump level. We show that spectral broadening may be eliminated by implementing output coupling of the signal wave and demonstrate 8.6 W of total signal and idler output with single-frequency spectra at both wavelengths for 14.5 W of pump power.     

Fast Monte Carlo inversion for extracting the optical properties of tubular tissues

Reconstruction of absorption coefficient μa and scattering coefficient μs is very important for applications of diffuse optical tomography and near infrared spectroscopy. Aiming at the early cancer detection of cervix and stomach, we present a fast inverse Monte-Carlo scheme for extracting μa and μs of a tubular tissue from the measurement on frequency domain. Results show that the computation time for reconstructing one set of μa and μs is less than 1 rain and the relative errors in reconstruction axe less than ±10% for the optical properties of normal cervical tissue and precancerous lesions.     

An efficient Jacobian scaling method for time-domain optical mammography

<正>Time-domain diffuse optical tomography can efficiently reconstruct optical parameters which can be further applied in diagnosing early breast cancer.Nevertheless,the performances of reconstructed imaging are badly influenced by different Jacobian magnitudes of absorption coefficient and reduced scattering coefficient.With the introudction of a relative data type based on generalized pulse spectrum technique, an efficient Jacobian scaling method is proposed.The interrelated simulated validation is also revealed for the enhancing performances.     

Diffuse-reflectance spectroscopy from 500 to 1060 nm by correction for inhomogeneously distributed absorbers

Diffuse-reflectance spectroscopy for measurement of the absorption and scattering coefficients of biological tissue produces reliable results for wavelengths from 650 to 1050 nm. Implicitly, this approach assumes homogeneously distributed absorbers. A correction factor is introduced for inhomogeneous distribution of blood concentrated in discrete cylindrical vessels. This factor extends the applicability of diffusion theory to lower wavelengths. We present measurements of in vivo optical properties in the wavelength range 500-1060 nm.     

Correlation transfer and diffusion of ultrasound-modulated multiply scattered light

We develop a temporal correlation transfer equation (CTE) and a temporal correlation diffusion equation (CDE) for ultrasound-modulated multiply scattered light. These equations can be applied to an optically scattering medium with embedded optically scattering and absorbing objects to calculate the power spectrum of light modulated by a nonuniform ultrasound field. We present an analytical solution based on the CDE and Monte Carlo simulation results for light modulated by a cylinder of ultrasound in an optically scattering slab. We further validate with experimental measurements the numerical calculations for an actual ultrasound field. The CTE and CDE are valid for moderate ultrasound pressures and on a length scale comparable with the optical transport mean-free path. These equations should be applicable to a wide spectrum of conditions for ultrasound-modulated optical tomography of soft biological tissues.