Search

The noninvasive reconstruction of optical parameter fields in participating phantom is investigated experimentally based on the measured time-resolved radiative signals, which are detected by the time-correlated single-photon counting system. The discrete ordinate method is employed to solve the forward transient radiative transfer equation to simulate the time-resolved radiative transfer in the physical phantom exposed to ultra-short pulse laser irradiation. On top of that, the sequential quadratic programming algorithm based on the generalized Gaussian Markov random field is applied to solve the inverse problem. To improve the computational efficiency, an adjoint equation technique is adopted to determine the gradient of objective function. Good agreement was obtained between the predicted optical parameter fields and realistic ones. All the reconstruction results demonstrate that the proposed reconstruction methods are proved to be accurate and efficient experimentally.

A reconstruction scheme for the fluorescence tomography is investigated based on the time-domain radiative transfer equation (TD-RTE). Two coupled TD-RTEs, which can provide considerable measurement data, are used as the forward model and solved by the discrete ordinate method. The sequential quadratic programming (SQP) is employed to build the reconstruction scheme for solving the inverse problem. The gradient of objective function is calculated efficiently by the adjoint equation technique. Considering the ill-posed nature of the inverse problem, the regularization term based on the generalized Gaussian Markov random field (GGMRF) model is adopted to enhance the reconstructed image. Influence of the initial guess, contrast, noisy data, and shape of the fluorescent target are analyzed. Simulated results show that the proposed algorithm performs efficiently and accurately on reconstructing the distribution of the fluorescence yield.

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.

The non-destructive testing (NDT) of ceramic materials using mid-infrared ultrashort-pulse laser is investigated in this study. The discrete ordinate method is applied to solve the transient radiative transfer equation in 2D semitransparent medium and the emerging radiative intensity on boundary serves as input for the inverse analysis. The sequential quadratic programming algorithm is employed as the inverse technique to optimize objective function, in which the gradient of objective function with respect to reconstruction parameters is calculated using the adjoint model. Two reticulated porous ceramics including partially stabilized zirconia and oxide-bonded silicon carbide are tested. The retrieval results show that the main characteristics of defects such as optical properties, geometric shapes and positions can be accurately reconstructed by the present model. The proposed technique is effective and robust in NDT of ceramics even with measurement errors.

Simultaneous reconstruction of the radiative parameter distributions in participating media, based on the transient radiative transfer equation (TRTE), was investigated. The discrete ordinate method was employed to solve the direct TRTE for two-dimensional participating media. A multi-start iterative technique, based on the conjugate gradient (CG) method, was proposed to accelerate the convergence speed of the inverse procedure. The maximum a posteriori (MAP) estimation was used to overcome the retrieval difficulties resulting from the ill-posed nature of the inverse problem. The gradient of the objective function was obtained using adjoint differentiation. All the simulation results show that MAP combined with the multi-start iterative CG technique is effective in simultaneously reconstructing the distributions of absorption and scattering coefficients in two-dimensional inhomogeneous media, even with noisy data.

A maximum a posteriori (MAP) estimation based on Bayesian framework is applied to image reconstruction of two-dimensional highly scattering inhomogeneous medium. The finite difference method (FDM) and conjugate gradient (CG) algorithm serve as the forward and inverse solving models, respectively. The generalized Gaussian Markov random field model (GGMRF) is treated as the regularization, and finally the influence of the measurement errors and initial distributions is investigated. Through the test cases, the MAP estimate algorithm is demonstrated to greatly improve the reconstruction results of the optical coefficients.

Reconstructing the distribution of optical parameters in the participating medium based on the frequency-domain radiative transfer equation (FD-RTE) to probe the internal structure of the medium is investigated in the present work. The forward model of FD-RTE is solved via the finite volume method (FVM). The regularization term formatted by the generalized Gaussian Markov random field model is used in the objective function to overcome the ill-posed nature of the inverse problem. The multi-start conjugate gradient (MCG) method is employed to search the minimum of the objective function and increase the efficiency of convergence. A modified adjoint differentiation technique using the collimated radiative intensity is developed to calculate the gradient of the objective function with respect to the optical parameters. All simulation results show that the proposed reconstruction algorithm based on FD-RTE can obtain the accurate distributions of absorption and scattering coefficients. The reconstructed images of the scattering coefficient have less errors than those of the absorption coefficient, which indicates the former are more suitable to probing the inner structure.

Sequential quadratic programming (SQP) is used as an optimization algorithm to reconstruct the optical parameters based on the time-domain radiative transfer equation (TD-RTE). Numerous time-resolved measurement signals are obtained using the TD-RTE as forward model. For a high computational efficiency, the gradient of objective function is calculated using an adjoint equation technique. SQP algorithm is employed to solve the inverse problem and the regularization term based on the generalized Gaussian Markov random field (GGMRF) model is used to overcome the ill-posed problem. Simulated results show that the proposed reconstruction scheme performs efficiently and accurately.