一种改进的光学层析图像重建方法

Alexander D．Klose将联合差分方法用于光学层析图像重建的梯度计算中,但给出的对光学参量的求导算法有局限,他的算法只能实现对边界点光学参量的导数计算,而无法实现对内部点光学参量导数的计算,会导致图像重建失败。在联合差分算法的基础上,研究了针对内部点光学参量的求导方法,给出了一种基于树形结构的对内部点光学参量求导的策略。具体实现时,为了降低计算复杂度,采用近似梯度计算方法。算法的仿真实验结果表明：该方法可以有效地实现对内部点光学参量的导数计算,提出的近似计算方法可降低梯度计算复杂度,提高运算速度,并可得到良好的图像重建质量。     

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.     

Accurate reconstruction of the optical parameter distribution in participating medium based on the frequency-domain radiative transfer equation

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.     

一种光学层析图像的多准则重建方法

光学层析成像是一个病态重建问题,为克服重建过程的病态性,提出将多准则优化理论引入到图像重建中。利用了三个用于光学层析图像重建的准则:平方误差函数、图像熵函数和局部平滑函数。采用向量优化方法将多准则优化问题转化为单准则优化问题求解。为了确定各个目标函数间的权重系数,提出一种动态权重系数求解方法。重建过程目标函数关于光学参量的梯度计算是关键,因此提出一种基于梯度树的计算方法。实验过程中对多准则重建结果和基于平方误差函数的单准则重建结果做了比较,证明该方法能够克服传统的偏重单一目标的单准则重建的不足,有效地重建光学层析图像,提高图像重建质量。     

Numerical developments for short-pulsed Near Infra-Red laser spectroscopy. Part II: inverse treatment

Indirect optical spectroscopy or tomography, that is, mapping of optical properties in scattering and absorption inside a medium given a set of measurements at the boundaries, is highly dependent on the radiative transfer model used to track radiative energy propagation in semi-transparent materials. In the first part of this study, a numerical tool adapted for treating radiative transfer in the frame of short-pulsed laser beam interaction with non-homogeneous matter has been presented. In this paper, it is intended to show how such numerical tools can undergo inversion through adjoint treatment or reverse differentiation.Adjoint models, as well as reverse differentiation, are used in order to allow an efficient computation of the gradient, in the unknown optical parameters space, of an objective or cost function estimating the residual between data obtained at the boundary and predictions by numerical simulations. This gradient is a crucial indication as to update, through line minimization, the set of internal optical properties of the medium.First, the theoretical background of the inverse treatments, both reverse differentiation and adjoint model, for the transient radiative transfer equation model introduced in Part I is developed. Second, different reconstruction configurations are presented. Time-dependent sampling and time filtering effects of the measurements are addressed. Image reconstructions from simulated data are achieved for material phantoms of simple geometry.     

Optical tomography using the time-independent equation of radiative transfer—Part 2: inverse model

Optical tomography is a novel imaging modality that is employed to reconstruct cross-sectional images of the optical properties of highly scattering media given measurements performed on the surface of the medium. Recent advances in this field have mainly been driven by biomedical applications in which near-infrared light is used for transillumination and reflectance measurements of highly scattering biological tissues. Many of the reconstruction algorithms currently utilized for optical tomography make use of model-based iterative image reconstruction (MOBIIR) schemes. The imaging problem is formulated as an optimization problem, in which an objective function is minimized. In the simplest case the objective function is a normalized-squared error between measured and predicted data. The predicted data are obtained by using a forward model that describes light propagation in the scattering medium given a certain distribution of optical properties.In part I of this two-part study, we presented a forward model that is based on the time-independent equation of radiative transfer. Using experimental data we showed that this transport-theory-based forward model can accurately predict light propagation in highly scattering media that contain void-like inclusions. In part II we focus on the details of our image reconstruction scheme (inverse model). A crucial component of this scheme involves the efficient and accurate determination of the gradient of the objective function with respect to all optical properties. This calculation is performed using an adjoint differentiation algorithm that allows for fast calculation of this gradient. Having calculated this gradient, we minimize the objective function with a gradient-based optimization method, which results in the reconstruction of the spatial distribution of scattering and absorption coefficients inside the medium. In addition to presenting the mathematical and numerical background of our code, we present reconstruction results based on experimentally obtained data from highly scattering media that contain void-like regions. These types of media play an important role in optical tomographic imaging of the human brain and joints.     

New developments in frequency domain optical tomography. Part I: Forward model and gradient computation

This two part study introduces new developments in frequency domain optical tomography to take into account the collimated source direction in the computation of both the forward and the adjoint models. The solution method is based on the least square finite element method associated to the discrete ordinates method where no empirical stabilization is needed. In this first part of the study, the solution method of the forward model is highlighted with an easy handling of complex boundary condition through a penalization method. Gradient computation from an adjoint method is developed rigorously in a continuous manner through a lagrangian formalism for the deduction of the adjoint equation and the gradient of the objective function. The proposed formulation can be easily generalized to stationary and time domain optical tomography by keeping the same expressions.     

基于对偶数理论的资料同化新方法

针对变分资料同化中目标泛函梯度计算精度不高且复杂等问题, 提出了一种基于对偶数理论的资料同化新方法, 主要优点是: 能避免复杂的伴随模式开发及其逆向积分, 只需在对偶数空间通过正向积分就能同时计算出目标泛函和梯度向量的值. 首先利用对偶数理论把梯度分析过程转换为对偶数空间中目标泛函计算过程, 简单、高效和高精度地获得梯度向量值; 其次结合典型的最优化方法, 给出了非线性物理系统资料同化问题的新求解算法; 最后对Lorenz 63混沌系统、包含开关的不可微物理模型和抛物型偏微分方程分别进行了资料同化数值实验, 结果表明: 新方法能有效和准确地估计出预报模式的初始条件或物理参数值.     

基于衍射光栅曲率波前传感器的实时波前重构研究

自适应光学系统要求波前传感器能实现动态实时测量,曲率波前传感技术符合这一发展要求。一种新型的基于扭曲衍射光栅的曲率波前传感器在探测装置的实现方法方面具有较大优势,其波前重构已应用于光学度量。根据衍射光学理论,对其探测信号进行数值模拟,并利用Neumann边界条件的Green函数法对其波前重构进行数值模拟。结果表明：Green函数法归结为2矩阵相乘,计算速度快,达到实时重构要求; Green函数法对阶数不高的Zernike多项式重构效果较好;影响重构误差的主要因素是光强梯度的边界噪声。     

Adjoint wall functions: A new concept for use in aerodynamic shape optimization

The continuous adjoint method for the computation of sensitivity derivatives in aerodynamic optimization problems of steady incompressible flows, modeled through the k–ε turbulence model with wall functions, is presented. The proposed formulation leads to the adjoint equations along with their boundary conditions by introducing the adjoint to the friction velocity. Based on the latter, an adjoint law of the wall that bridges the gap between the solid wall and the first grid node off the wall is proposed and used during the solution of the system of adjoint (to both the mean flow and turbulence) equations. Any high Reynolds turbulence model, other than the k–ε one used in this paper, could also profit from the proposed adjoint wall function technique. In the examined duct flow problems, where the total pressure loss due to viscous effects is used as objective function, emphasis is laid on the accuracy of the computed sensitivity derivatives, rather than the optimization itself. The latter might rely on any descent method, once the objective function gradient has accurately been computed.     

Aerodynamic design optimization by using a continuous adjoint method

This paper presents the fundamentals of a continuous adjoint method and the applications of this method to the aerodynamic design optimization of both external and internal flows.General formulation of the continuous adjoint equations and the corresponding boundary conditions are derived.With the adjoint method,the complete gradient information needed in the design optimization can be obtained by solving the governing flow equations and the corresponding adjoint equations only once for each cost function,regardless of the number of design parameters.An inverse design of airfoil is firstly performed to study the accuracy of the adjoint gradient and the effectiveness of the adjoint method as an inverse design method.Then the method is used to perform a series of single and multiple point design optimization problems involving the drag reduction of airfoil,wing,and wing-body configuration,and the aerodynamic performance improvement of turbine and compressor blade rows.The results demonstrate that the continuous adjoint method can efficiently and significantly improve the aerodynamic performance of the design in a shape optimization problem.     

Optical tomography with the discontinuous Galerkin formulation of the radiative transfer equation in frequency domain

Optical tomography is an inverse method of probing semi-transparent media with the help of light sources. The reconstruction of the optical properties usually employs finite volumes or continuous finite elements formulations of light transport as a forward model for the predictions. In a previous study, we have introduced a generalization of the inversion approach with finite elements formulations by using an integral form of the objective function. The novelty is that the surfaces of the detectors are taken into account in the reconstruction and compatibility is obtained for all finite element formulations. This present paper illustrates this new approach by developing a Discontinuous Galerkin formulation as a forward model for an optical tomography application in the frequency domain framework. Numerical tests are performed to gauge the accuracy of the method in recovering optical properties distribution with a gradient-based algorithm where the adjoint method is used to fastly compute the objective function gradient. It is seen that the reconstruction is accurate and can be affected by noise on the measurements as expected. Filtering of the gradient at each iteration of the reconstruction is used to cope with the ill-posed nature of the inverse problem and to improves the quality and accuracy of the reconstruction.     

变分伴随正则化方法从雷达回波反演海洋波导 Ⅰ.理论推导部分

针对传统统计反演算法在雷达回波反演海洋波导(RFC)方面计算量过大的问题,提出一种变分伴随正则化物理反演算法.在变分伴随方法中分别导出切线性模式、伴随方程及伴随边界条件、伴随方程求解表达式、以及泛函梯度的数学表达形式;讨论了伴随方法求泛函梯度对复方程如何协调的问题.考虑到雷达电磁波传播的特征,选择了适当的正则化项来解决反演中的不适定问题.最后给出反演算法实施中的迭代格式.     

动态背景下基于光流场分析的运动目标检测算法

针对现有动态背景下运动目标检测算法的不足,提出一种基于光流场分析的运动目标检测算法.首先根据前背景在光流梯度幅值和光流矢量方向上的差异确定目标的大致边界,然后通过点在多边形内部原理获得边界内部的稀疏像素点,最后以超像素为节点,利用混合高斯模型拟合的表观信息和超像素的时空邻域关系构建马尔可夫随机场模型的能量函数,并通过使目标函数能量最小化得到最终的运动目标检测结果.该算法不需要任何先验假设,能够同时处理动态背景和静态背景两种情况.多组实验结果表明,本文算法在检测的准确性和处理速度上均优于现有算法.     

光强传播方程相位恢复的实验验证

在傍轴近似条件下, 可以利用光强传播方程(ITE)对畸变波前进行相位恢复。对于衍射受限的光学系统, 很难获得相位的边界径向斜率值作为边界条件, 此外, 要获得精确的圆域边界采样值也并非易事。为了克服上述困难, 进一步研究了相位恢复改进方法, 即改变了方程的表示形式、计算区域和边界条件, 并用多重网格法进行求解获得重构相位, 最后再将其修正。为了验证算法的精确性, 搭建了实验系统对算法进行实验验证, 将由CCD探测的光强分布进行计算得到的相位分布与相位恢复算法(PR算法)的结果进行比较, 证明在光强分布非均匀的情况下, 该方法恢复的相位均方根误差能够达到017 λ, 可以适用于波前传感精度要求不是很高的相位恢复。     

Elastic-wave identification of penetrable obstacles using shape-material sensitivity framework

This study deals with elastic-wave identification of discrete heterogeneities (inclusions) in an otherwise homogeneous “reference” solid from limited-aperture waveform measurements taken on its surface. On adopting the boundary integral equation (BIE) framework for elastodynamic scattering, the inverse query is cast as a minimization problem involving experimental observations and their simulations for a trial inclusion that is defined through its boundary, elastic moduli, and mass density. For an optimal performance of the gradient-based search methods suited to solve the problem, explicit expressions for the shape (i.e. boundary) and material sensitivities of the misfit functional are obtained via the adjoint field approach and direct differentiation of the governing BIEs. Making use of the message-passing interface, the proposed sensitivity formulas are implemented in a data-parallel code and integrated into a nonlinear optimization framework based on the direct BIE method and an augmented Lagrangian whose inequality constraints are employed to avoid solving forward scattering problems for physically inadmissible (or overly distorted) trial inclusion configurations. Numerical results for the reconstruction of an ellipsoidal defect in a semi-infinite solid show the effectiveness of the proposed shape-material sensitivity formulation, which constitutes an essential computational component of the defect identification algorithm.     

基于梯度阈值自适应处理的红外图像超分辨率重建

超分辨率图像重建中,Huber马尔可夫随机场模型是一种常用的正则化算子.针对Huber函数中固定梯度阈值引起图像重建效果不佳的问题,本文提出一种梯度阈值自适应处理的红外图像超分辨率重建算法.在最大后验概率理论框架下,构造了基于数据项和正则项的正则化模型;通过迭代的方式,利用中间重建结果不断更新正则化参量,解决了Huber马尔可夫随机场模型中梯度阈值不易选择的难题.实验结果表明,改进算法能够根据局部梯度特征自适应选择相应的正则化参量并找到最优解,较好恢复目标细节的同时有效抑制了图像噪音.     

基于梯度阈值自适应处理的红外图像超分辨率重建

超分辨率图像重建中,Huber马尔可夫随机场模型是一种常用的正则化算子.针对Huber函数中固定梯度阈值引起图像重建效果不佳的问题,本文提出一种梯度阈值自适应处理的红外图像超分辨率重建算法.在最大后验概率理论框架下,构造了基于数据项和正则项的正则化模型;通过迭代的方式,利用中间重建结果不断更新正则化参量,解决了Huber马尔可夫随机场模型中梯度阈值不易选择的难题.实验结果表明,改进算法能够根据局部梯度特征自适应选择相应的正则化参量并找到最优解,较好恢复目标细节的同时有效抑制了图像噪音.     

Geometric significance of the spinor Lie derivative. I

In a previous article, the writer explored the geometric foundation of the generally covariant spinor calculus. This geometric reasoning can be extended quite naturally to include the Lie covariant differentiation of spinors. The formulas for the Lie covariant derivatives of spinors, adjoint spinors, and operators in spin space are deduced, and it is observed that the Lie covariant derivative of an operator in spin space must vanish when taken with respect to a Killing vector. The commutator of two Lie covariant derivatives is calculated; it is noted that the result is consistent with the geometric interpretation of the Jacobi identity for vectors. Lie current conservation is seen to spring from the result that the operator of spinor affine covariant differentiation commutes with the operator of spinor Lie covariant differentiation with respect to a Killing vector. It is shown that differentiations of the spinor field defined geometrically are Lorentz-covariant.