Image reconstruction method for laminar optical tomography with only a single Monte-Carlo simulation

Laminar optical tomography(LOT)is a new mesoscopic functional optical imaging technique.Currently,the forward problem of LOT image reconstruction is generally solved on the basis of Monte-Carlo(MC)methods.However,considering the nonlinear nature of the image reconstruction in LOT with the increasing number of source positions,methods based on MC take too much computation time.This letter develops a fast image reconstruction algorithm based on perturbation MC(pMC)for reconstructing the absorption or scattering image of a slab medium,which is suitable for LOT or other functional optical tomography system with narrow source-detector separation and dense sampling.To calculate the pMC parameters,i.e.,the path length passed by a photon and the collision numbers experienced in each voxel with only one baseline MC simulation,we propose a scheme named as the trajectory translation and target voxel regression(TTTVR)based on the reciprocity principle.To further speed up the image reconstruction procedure,the weighted average of the pMC parameters for all survival photons is adopted and the region of interest(ROI)is extracted from the raw data to save as the prior information of the image reconstruction.The method is applied to the absorption reconstruction of the layered inhomogeneous media.Results demonstrate that the reconstructing time is less than 20 s with the X-Y section of the sample subdivided into 50×50 voxels,and the target size quantitativeness ratio can be obtained in a satisfying accuracy in the source-detector separations of 0.4 and 1.25 mm,respectively.     

融合结构先验信息的稳态扩散光学断层成像重建算法研究

扩散光学断层成像作为一种无辐射损伤、低成本的光学在体成像技术,有着良好的应用前景,但具有空间分辨率低、难以定量的缺陷.为了提高扩散光学断层成像的分辨率,实现光学参数分布的精确重建,基于有限元方法,提出了融合结构先验信息的稳态扩散光学断层成像重建算法.该算法以扩散近似作为成像模型,通过软先验的Laplace 正则化方法引入由MicroCT提供的空间结构信息.采用伴随法计算Jacobian矩阵,Levenberg-Marquardt方法用来进行迭代优化.仿真结果表明该算法不仅能获得精确的光学参数值分布,而且显著地提高了迭代收敛的速度.     

Monte-Carlo simulation of optical trap stiffness measurement

The high precision calibration of optical trap stiffness is the foundation of the weak force measurement in an optical tweezers system. And the accuracy of the trap stiffness measurement is limited by the bandwidth of the acquisition system. In this article, such an influence is analyzed and discussed. The stiffness measuring process using an acquisition system with a finite acquisition time is numerically simulated by using Monte-Carlo method. Then the simulated results are analyzed by thermal motion analysis method to deduce the trap stiffness for different trapping system and for measuring systems with different acquisition time. As a comparison the power spectrum analysis method is used to study the thermal motion of the bead and to compute the trap stiffness for the same acquisition system, from which it is concluded that the bandwidth of the acquisition system is determined by its acquisition time, not the sampling frequency. The influence of the finite acquisition time or the limited bandwidth on the trap stiffness measurement is discussed. The numerical simulation shows that the measured position, which is here the average position within the acquisition time, shifts to the trap center due to the trapping force, which gives an alternative interpretation for the deviation of the measured stiffness from the true trap stiffness.     

Local diffusion regularization method for optical tomography reconstruction by using robust statistics

We formulate a solution to the diffuse optical tomography (DOT) inverse problem as the minimization of an energy functional of the solution and the data. For the solution prior we introduce a local diffusion regularization potential with a threshold based on robust statistics (the Hubert function). We compare results on simulated data for the Hubert function and two other standard regularization functionals, Tikhonov and total variation.     

An efficient and robust reconstruction method for optical tomography with the time-domain radiative transfer equation

An efficient and robust method based on the complex-variable-differentiation method (CVDM) is proposed to reconstruct the distribution of optical parameters in two-dimensional participating media. An upwind-difference discrete-ordinate formulation of the time-domain radiative transfer equation is well established and used as forward model. The regularization term using generalized Gaussian Markov random field model is added in the objective function to overcome the ill-posed nature of the radiative inverse problem. The multi-start conjugate gradient method was utilized to accelerate the convergence speed of the inverse procedure. To obtain an accurate result and avoid the cumbersome formula of adjoint differentiation model, the CVDM was employed to calculate the gradient of objective function with respect to the optical parameters. All the simulation results show that the CVDM is efficient and robust for the reconstruction of optical parameters.     

Image reconstruction in diffuse optical tomography using the coupled radiative transport-diffusion model

The coupled radiative transport-diffusion model can be used as light transport model in situations in which the diffusion equation is not a valid approximation everywhere in the domain. In the coupled model, light propagation is modelled with the radiative transport equation in sub-domains in which the approximations of the diffusion equation are not valid, such as within low-scattering regions, and the diffusion approximation is used elsewhere in the domain. In this paper, an image reconstruction method for diffuse optical tomography based on using the coupled radiative transport-diffusion model is developed. In the approach, absorption and scattering distributions are estimated by minimising a regularised least-squares error between the measured data and solution of the coupled model. The approach is tested with simulations. Reconstructions from different cases including domains with low-scattering regions are shown. The results show that the coupled radiative transport-diffusion model can be utilised in image reconstruction problem of diffuse optical tomography and that it produces as good quality reconstructions as the full radiative transport equation also in the presence of low-scattering regions.     

Monte-Carlo simulation of effective stiffness of time-sharing optical tweezers

<正>The Brownian motion of a polystyrene bead trapped in a time-sharing optical tweezers(TSOT) is numerically simulated by adopting Monte-Carlo technique.By analyzing the Brownian motion signal,the effective stiffness of a TSOT is acquired at different switching frequencies.Simulation results confirm that for a specific laser power and duty ratio,the effective stiffness varies with the frequency at low frequency range,while at high frequency range it keeps constant.Our results reveal that the switching frequency can be used to control the stability of time-sharing optical tweezers in a range.     

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.     

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.     

利用蒙特卡洛法分析红外光学系统的杂散辐射

介绍了杂散辐射的概念及其对光学系统成像质量的影响。以一透视式红外系统为例,通过对系统的光机结构的合理简化,利用杂散辐射分析软件Light tools中已有的蒙特卡洛分析方法,对已建立的红外系统光机结构的三维模型进行光线追迹分析,得到了系统像面的照度分布结果。同时,基于辐射度学相关理论对该系统中各结构表面的热辐射进行了定量计算,分析了结构中不同组成部分对系统像面的杂散辐射影响,提出本文系统中前镜框内表面为主要杂散辐射源。针对分析结果,探讨了红外系统的杂散辐射抑制方法,提出的方法可为后续设计提供参考。     

层析成像图像重建算法综述

介绍了层析成像技术的图像重建算法,并从正向问题数学模型的简化和反向问题数学模型的映射结构的角度比较了各种算法的特点和优劣。研究表明:用本质是线性算法的各种变换方法重建图像存在严重失真,而卷积滤波的引入可以使变换方法的重建效果有所改善;基于导数搜索的迭代算法对初始值依赖性强、收敛速度慢并且容易陷入局部最优解;基于Fourier变换的方法具有本质的局限性;小波变换则可以同时刻画图像时域和频域的细节特征;有限元法通过重建对象像素的智能划分可以简化正问题的复杂性;而具有物理背景的蒙特卡罗法、模拟退火法、遗传算法、粒子滤波法及神经网络法更适合于复杂且非线性的图像重建;智能化、仿生化、并行化以及各种算法的融合是层析成像图像重建算法的发展趋势。     

Image reconstruction algorithm for diffraction enhanced imaging-based computed tomography

A number of reconstruction algorithms for differential phase-contrast computed tomography reconstruct the refractive index decrement using the refraction angles of X-rays. However, these approaches require that the refraction angle at each view be extracted from several raw images captured at the same view, leading to unacceptably long exposure times and huge X-ray doses. We develop an iterative reconstruction algorithm to reconstruct the complex refractive index from raw images gathered only once at each view angle. Using the Tikhonov regularization method as basis, we simultaneously obtain the refractive index decrement and attenuation coefficient by minimizing the cost function. The experimental results show that the proposed algorithm effectively yields reconstructed images of low noise level.     

Image reconstruction in optoacoustic tomography for dispersive acoustic media

Conventional image reconstruction methods for optoacoustic tomography (OAT) assume an idealized, non-dispersive acoustic medium. However, the linear attenuation coefficient and the phase velocity of acoustic waves propagating in soft tissue depend on temporal frequency and satisfy a known dispersion law. These frequency-dependent effects are incorporated into an optoacoustic wave equation, and a corresponding reconstruction method for OAT is developed. The improvement in image fidelity that can be achieved over conventional reconstruction methods is demonstrated by use of computer-simulation studies.     

基于全变分最小化和交替方向法的康普顿散射成像重建算法

康普顿散射成像技术利用射线与物质作用后的散射光子信息对物质的电子密度进行成像.与传统的透射成像方式相比,康普顿散射成像具有系统结构灵活、成像对比度高、辐射剂量低等优势,在无损检测、医疗诊断、安全检查等领域有着广阔的应用前景.但其重建问题是一个非线性的逆问题,通常是不适定的,其解对噪声和测量误差非常敏感.为解决此问题,本文结合全变分最小化正则化方法和交替方向法提出了一种新的康普顿散射成像重建算法.该算法首先将问题对应的TV模型转化为与之等价的带约束的优化问题,然后利用增广拉格朗日乘子法将优化问题分解为两个具有解析解的子问题,并通过交替求解子问题使增广拉格朗日函数达到最小,进而得到重建的图像.在仿真实验中,通过与主流的ASD-POCS方法进行对比,证明了该算法在重建精度和重建效率方面的优势.     

Simplified method for polarization-sensitive optical coherence tomography

We report a method for extracting the birefringence properties of biological samples with micrometer-scale resolution in three dimensions, using a new form of polarization-sensitive optical coherence tomography. The method measures net retardance, net fast axis, and total reflectivity as a function of depth into the sample. Polarization sensing is accomplished by illumination of the sample with at least three separate polarization states during consecutive acquisitions of the same pixel, A scan, or B scan. The method can be implemented by use of non-polarization-maintaining fiber and a single detector. In a calibration test of the system, net retardance was measured with an average error of 7.5 degrees (standard deviation 2.2 degrees ) over the retardance range 0 degrees to 180 degrees , and a fast axis with average error of 4.8 degrees over the range 0 degrees to 180 degrees .     

Sparse Bayesian reconstruction method for multispectral bioluminescence tomography

We present a sparse Bayesian reconstruction method based on multiple types of a priori information for multispectral bioluminescence tomography （BLT）. In the Bayesian approach, five kinds of a priori information are incorporated, reducing the ill-posedness of BLT. Specifically, source sparsity characteristic is considered to promote reconstruction results. Considering the computational burden in the multispectral case, a series of strategies is adopted to improve computational efficiency, such as optimal permissible source region strategy and node model of the finite element method. The performance of the proposed algorithm is validated by a heterogeneous three-dimensional （3D） micron scale computed tomography atlas and a mouse-shaped phantom. Reconstructed results demonstrate the feasibility and effectiveness of the proposed algorithm.     

Two-dimensional scanning probe driven by a solenoid-based single actuator for optical coherence tomography

We propose the two-dimensional scanning probe operating with a single actuator, which is thought useful as a sample probe for optical coherence tomography (OCT). The probe was designed to use a single-body lensed fiber cantilever loaded with an iron-bead and driven by a single-solenoid actuator. Elliptic spiral trace patterns were achieved using off-axis magnetic fields of the solenoid. A three-dimensional OCT image was obtained for a scanning area of 3.8 mm × 3.4 mm at an acquisition speed of 16.7 s/V. Up to 27 Hz B-scan rate, the proposed probe worked well, and 1000 A-scans were made per each B-scan.     

Early-photon guided reconstruction method for time-domain fluorescence lifetime tomography

A reconstruction method guided by early-photon fluorescence yield tomography is proposed for time-domain fluorescence lifetime tomography(FLT) in this study. The method employs the early-arriving photons to reconstruct a fluorescence yield map, which is utilized as a priori information to reconstruct the FLT via all the photons along the temporal-point spread functions. Phantom experiments demonstrate that, compared with the method using all the photons for reconstruction of fluorescence yield and lifetime maps, the proposed method can achieve higher spatial resolution and reduced crosstalk between different targets without sacrificing the quantification accuracy of lifetime and contrast between heterogeneous targets.     

Reconstruction method for samples with refractive index discontinuities in optical diffraction tomography

We present a reconstruction method for samples containing localized refractive index (RI) discontinuities in optical diffraction tomography. Abrupt RI changes induce regional phase perturbations and random spikes, which will be expanded and strengthened by existing tomographic algorithms, resulting in contaminated reconstructions. This method avoids the disturbance by recognition and separation of the discontinuous regions, and recombination of individually reconstructed data. Three-dimensional RI distributions of two fusion spliced optical fibers with different typical discontinuities are demonstrated, showing distinctly detailed structures of the samples as well as the positions and estimated shapes of the discontinuities.     

Micro-spectrometry in the visible range with full-field optical coherence tomography for single absorbing layers

A time-domain full-field OCT adapted to the visible range and with an original configuration using an interferometric objective, that minimizes mechanical vibrations and some settings and that performs imaging without moving the sample, is presented. This setup achieves micrometer scale imaging, 1.5 μm in the axial direction and 1.2 μm in the lateral one. The principle of micro-spectrometry from OCT data by Fourier transform is described and the influence of some key data processing parameters is simulated and discussed. The experimental spectra reconstruction from tomographic data is validated by comparison with transmittance spectra. Imaging and spectra of dyes at a micrometer scale are obtained from the same data volume.