超声层析成像的迭代重建

利用变形Born迭代方法,建立了超声衍射重建算法。在迭代过程中,为了解决超声逆散射问题中的非线性性,需要反复地求解前向散射方程和逆散射方程,以达到全场和未知函数的近似,较好地重建物体内部的断层图象。由于逆散射方程是一个不适定性的方程组,要用正则化方法处理方程的不适定性问题,使迭代方法收敛于问题的真实解,才能成功地应用于较高对比度物体的图象重建问题。用Picard准则对不适定问题进行了分析,给出了通过简单图形．确定模型受噪声污染情况以及正则化方法适用范围的方法。在重建实验中。对建立的图像重建算法进行了实验仿真。达到了较好的效果。     

Solution of convergence problem in ultrasound inverse scattering tomography

When the product of contrast and size of an object, which is to be reconstructed by using the ultrasound inverse scattering tomography algorithm, is large, it is well known that those algorithms fail to converge to a unique global minimum. In order to solve this well known and difficult convergence problem, in this paper we present a new method, which converges to the true solution, for obtaining the scattering potential without using the Born or Rytov approximation. This method converts the nonlinear nature of the problem into a linear one. Through computer simulations we will show the validity of the new approach for high contrast two-dimensional scattering objects which are insonified by an incident ultrasound plane wave. Numerical results show that the reconstruction error is very small for circularly symmetric two-dimensional cylindrical objects whose refractive indices range from small to even sufficiently large values for which the previous inverse scattering algorithms fail to converge.     

Ultrasound imaging using variations of the iterative Born technique [biomedical diagnosis

The iterative Born method is an inverse technique that has been used successfully in ultrasound imaging. However, the calculation cost of the standard iterative Born method is high, and parallel computation is limited to the forward problem. In this work, two methods are introduced to increase the rate of convergence of the iterative Born algorithm. These methods are tested on three different objects. The results are promising, with both algorithms giving accurate results at lower computational cost. The first method, referred to as the coarse resolution initial value (CRIV) method, uses the iterative Born algorithm for a coarse grid to quickly estimate the initial value of the object to be reconstructed. From this initial value, the final image is obtained for a finer grid with additional iterations. The cost of this method is 40% less than that of the iterative Born technique. The second method, the quadriphase source (QS) method, simultaneously uses four single sources, and object reconstruction for each is performed in parallel; the reconstruction results for all four sources then are averaged to obtain the final image. The cost of this method is 20% less than that of the standard iterative Born method. When the object to be reconstructed is of low contrast and/or has a small phase shift, the QS method is very promising because parallel computation can be used to solve both the forward and inverse problems. However, the QS method fails for high contrast objects.     

An iterative solution of the two-dimensional electromagnetic inverse scattering problem

A new method, based on an iterative procedure, for solving the two-dimensional inverse scattering problem is presented. This method employs an equivalent Neumann series solution in each iteration step. The purpose of the algorithm is to provide a general method to solve the two-dimensional imaging problem when the Born and the Rytov approximations break down. Numerical simulations were calculated for several cases where the conditions for the first order Born approximation were not satisfied. The results show that in both high and low frequency cases, good reconstructed profiles and smoothed versions of the original profiles can be obtained for smoothly varying permittivity profiles (lossless) and discontinuous profiles (lossless), respectively. A limited number of measurements around the object at a single frequency with four to eight plane incident waves from different directions are used. The method proposed in this article could easily be applied to the three-dimensional inverse scattering problem, if computational resources are available.     

Density imaging using a multiple-frequency DBIM approach

Current inverse scattering methods for quantitative density imaging have limitations that keep them from practical experimental implementations. In this work, an improved approach, termed the multiple-frequency distorted Born iterative method (MF-DBIM) algorithm, was developed for imaging density variations. The MF-DBIM approach consists of inverting the wave equation by solving for a single function that depends on both sound speed and density variations at multiple frequencies. Density information was isolated by using a linear combination of the reconstructed single-frequency profiles. Reconstructions of targets using MF-DBIM from simulated data were compared with reconstructions using methods currently available in the literature, i.e., the dual-frequency DBIM (DF-DBIM) and T-matrix approaches. Useful density reconstructions, i.e., root mean square errors (RMSEs) less than 30%, were obtained with MF-DBIM even with 2% Gaussian noise in the simulated data and using frequency ranges spanning less than an order of magnitude. Therefore, the MFDBIM approach outperformed both the DF-DBIM method (which has problems converging with noise even an order of magnitude smaller) and the T-matrix method (which requires a ka factor close to unity to achieve convergence). However, the convergence of all the density imaging algorithms was compromised when imaging targets with object functions exhibiting high spatial frequency content.     

Multifrequency near-field acoustic tomography and holography of 3D subbottom inhomogeneities

A method of coherent multifrequency acoustic tomography and holography of spatially localized subbottom inhomogeneities in shallow seas is proposed. This method is based on solving of the near-field inverse scattering problem that makes it possible to realize a subwavelength resolution. It involves the analysis of measurement data obtained by the 2D transversal scanning with the source-receiver system along the sea bottom, over the area of sounded inhomogeneities. The solution begins with the Born approximation, where the original 3D integral equation for the scattered field is reduced to the 1D Fredholm equation of the first kind relative to the depth profile of the lateral spectrum of inhomogeneities. When solving this integral equation for each pair of spectral components, the generalized discrepancy method is in use. Then, corrections to the Born approximation can be obtained in the proposed iterative procedure. For distributed inhomogeneities, the inverse Fourier transform of the retrieved spectrum gives their 3D distribution that can be visualized as tomography images. For solid targets, this spectrum is used to obtain their shape (i.e. to solve the problem of computer holography). Corresponding results of the numerical simulation are presented.     

超声逆散射成像问题中正则化参数研究

研究了二维理想情况下,基于精确场描述的超声逆散射成像问题,先用矩量法将波动方程化为离散形式,分别用BI和DBI算法进行迭代重建。影响整个算法的一个关键因素是散射场方程的正则化求解,具有明显的不适定性。文章基于L曲线法,提出以解的范数和残差变化量的加权形式作为确定正则化参数的依据,在迭代过程根据问题不适定性程度,自适应地调整搜索范围。仿真结果表明,该算法可快速地找到最优正则化参数。     

Comparison of the born iterative method and tarantola's method for an electromagnetic time-domain inverse problem

Two methods of solving the nonlinear two-dimensional electromagnetic inverse scattering problem in the time domain are considered. These are the Born iterative method and the method originally proposed by Tarantola for the seismic reflection inverse problems. The former is based on Born-type iterations on an integral equation, whereby at each iteration the problem is linearized, and its solution is found via a regularized optimization. The latter also uses an iterative method to solve the nonlinear system of equations. Although it linearizes the problem at each stage as well, no optimization is carried out at each iteration; rather the problem as a whole is posed as a (regularized) optimization. Each method is described briefly and its computational complexity is analyzed. Tarantola's method is shown to have a lower numerical complexity compared to the Born iterative method for each iteration, but in the examples considered, required more iterations to converge. Both methods perform well when inverting a smooth profile; however, the Born iterative method gave better results in resolving localized point scatterers.     

基于同伦延拓的冗余机械臂逆运动学优化算法研究

目的 针对偏置冗余机械臂的逆运动学,采用传统数值法存在依赖初始值、奇异位姿收敛性差等问题,提出一种改进数值法。方法 首先将非线性方程组转化为同伦方程组,引入同伦延拓算法能够有效避免依赖初始值的问题,同时能够获取逆运动学解空间。然后考虑奇异位姿,将同伦方程组转化为最小二乘问题,采用Levenberg Marquardt算法对同伦方程组进行路径追踪,以获取逆运动学解空间。最后将关节极限避免问题映射为解空间优化问题,引入二进制改进粒子群优化算法,获得最优逆运动学解。结果 实验结果表明,相较于传统数值法,文中所提数值法针对逆运动学求解具有更高的收敛率、更快的收敛速度,同时二进制改进粒子群算法能够有效避免关节极限问题。结论 采用文中所提数值法求解逆运动学的精度较高,能够满足实时性要求,对于机械臂用于包装作业具有一定的理论意义和工程应用价值。     

二维弹性波散射时域Born近似

本文探讨了非零面积缺陷对一种典型脉冲波（称之为基脉冲）散射的正反问题。利用平面简谐波人射时弹性波散射场的频域Born近似解,通过Fourier变换,给出了散射体（或者说缺陷）对基脉冲入射时的脉冲响应函数（称之为基脉冲响应函数）,并通过Radon变换法建立了散射体特征函数与基脉冲响应函数的关系式,模拟识别结果表明时域Born近似方法对无损检测技术具有潜在的应用价值。     

MUSIC imaging applied to total internal reflection tomography

This paper applies the multiple signal classification (MUSIC) imaging method to determine the locations of a collection of small anisotropic spherical scatterers in the framework of total internal reflection tomography. Multiple scattering between scatterers is considered, and the inverse scattering problem is nonlinear, which, however, is solved by the proposed fast analytical approach where no associated forward problem is iteratively evaluated. The paper also discusses the role of propagating and evanescent waves, the polarization of incidence waves, separation of scatterers from the surface of the substrate, and the level of noise on the resolution of imaging.     

Nonlinearized approach to profile inversion

A method for reconstructing the index of refraction of a bounded inhomogeneous object of known geometric configuration from measured far-field scattering data is presented. This work is an extension of recent results on the direct scattering problem wherein the governing domain integral equation was solved iteratively by a successive relaxation technique. The relaxation parameters were chosen to minimize the residual error at each step. Convergence of this process was established for indices of refraction much larger than required for convergence of the Born approximation. For the inverse problem, the same technique is applied, except is this case both the index of refraction and the field are unknown. Iterative solutions for both unknowns are postulated with two relaxation parameters at each step. They are determined by simultaneously minimizing the residual errors in satisfying the domain integral equation and matching the measured data. This procedure retains the nonlinear relation between the two unknowns. Numerical results are presented for the dielectric slab. The algorithm is shown to be effective in cases where the iterative solution of the direct problem is rapidly convergent and outperforms the Born-based approaches.     

Nonparaxial vector-field modeling of optical coherence tomography and interferometric synthetic aperture microscopy

A large-aperture, electromagnetic model for coherent microscopy is presented and the inverse scattering problem is solved. Approximations to the model are developed for near-focus and far-from-focus operations. These approximations result in an image-reconstruction algorithm consistent with interferometric synthetic aperture microscopy (ISAM): this validates ISAM processing of optical-coherence-tomography and optical-coherence-microscopy data in a vectorial setting. Numerical simulations confirm that diffraction-limited resolution can be achieved outside the focal plane and that depth of focus is limited only by measurement noise and/or detector dynamic range. Furthermore, the model presented is suitable for the quantitative study of polarimetric coherent microscopy systems operating within the first Born approximation.     

Multi-parameter acoustic imaging of uniform objects in inhomogeneous soft tissue

The problem studied in this paper is ultrasound image reconstruction from frequency-domain measurements of the scattered field from an object with contrast in attenuation and sound speed. The case in which the object has uniform but unknown contrast in these properties relative to the background is considered. Background clutter is taken into account in a physically realistic manner by considering an exact scattering model for randomly located small scatterers that vary in sound speed. The resulting statistical characteristics of the interference are incorporated into the imaging solution, which includes application of a total-variation minimization-based approach in which the relative effect of perturbation in sound speed to attenuation is included as a parameter. Convex optimization methods provide the basis for the reconstruction algorithm. Numerical data for inversion examples are generated by solving the discretized Lippman-Schwinger equation for the object and speckle-forming scatterers in the background. A statistical model based on the Born approximation is used for reconstruction of the object profile. Results are presented for a two-dimensional problem in terms of classification performance and compared with minimum-l2-norm reconstruction. Classification using the proposed method is shown to be robust down to a signal-to-clutter ratio of less than 1 dB.     

Application of the sinc basis moment method to the reconstruction of infinite circular cylinders

A solution of the ultrasonic scattering and inverse scattering problem has been obtained by solving the inhomogeneous Helmholtz wave equation by the sinc basis moment method. In this numerical study, the algorithm of S.A. Johnson and M.L. Tracy (1983) has been applied to the reconstruction of an infinite circular cylinder that is subject to an incident cylindrical wave of ultrasound and is surrounded by a homogeneous coupling medium. For weak scattering cylinders, successful reconstructions have been obtained using the known exact solution for the scattered field as the input data for the algorithm. A detailed discussion of sampling requirements for this algorithm is presented, and the threshold derived correlates well with results of a numerical study of variation of the sampling density. Effects of varying object contrast, object size, grid size, sampling density, and method of iteration are investigated. Because the algorithm is slow, optimization of computation is described.     

The diagonalized contrast source inversion approach for elastic wave inversion

This study focuses on the inverse scattering of objects embedded in a homogeneous elastic background. The medium is probed by ultrasonic sources, and the scattered fields are observed along a receiver array. The goal is to retrieve the shape, location, and constitutive parameters of the objects through an inversion procedure. The problem is formulated using a vector integral equation. As is well-known, this inverse scattering problem is nonlinear and ill-posed. In a realistic configuration, this nonlinear inverse scattering problem involves a large number of unknowns, hence the application of full nonlinear inversion approaches such as Gauss-Newton or nonlinear gradient methods might not be feasible, even with present-day computer power. Hence, in this study we use the so-called diagonalized contrast source inversion (DCSI) method in which the nonlinear problem is approximately transformed into a number of linear problems. We will show that, by using a three-step procedure, the nonlinear inverse problem can be handled at the cost of solving three constrained linear inverse problems. The robustness and efficiency of this approach is illustrated using a number of synthetic examples.     

Multistatic-multiview resolution from Born fields for strips in Fresnel zone

We address the inverse scattering problem of estimating the resolution limits achievable in the reconstruction of a dielectric strip object within a two-dimensional and scalar geometry. The scattered field is observed over a bounded rectilinear domain located in the Fresnel zone, and a single-frequency multistatic-multiview configuration is considered. The analysis is performed by casting the problem as the inversion of the linearized scattering operator arising from the Born approximation and by means of its singular-value decomposition. Finally, the role of the geometrical parameters of the measurement configuration is highlighted.     

基于单频数据重构散射障碍的新型线性采样方法

散射及反散射的数学理论与计算一直是应用数学领域中的重要课题，其成果在地质勘探、无损探测、医学成像等领域都具有广泛的应用．线性采样方法（LSM）是近年来反散射理论中一类非常流行的非迭代型重建算法，但是这种方式很难推广到如半空间中障碍反散射等更为复杂的问题中．本文基于单频数据研究Dirichlet障碍反散射问题的数值重建算法．通过构造带有阻尼边界条件的辅助边值问题，提出了一类新型的线性采样方法，并在理论上严格证明了该方法在任意给定的波数下重构障碍形状及位置的有效性．     

A study on the reconstruction of moderate contrast targets using the distorted born iterative method

Previous tomographic methods using ultrasound for reconstructing sound speed and attenuation images suffered from convergence issues for targets with moderate speed of sound contrast. Convergence problems can be overcome by the use of the multiple frequency, distorted Born iterative method (DBIM). The implementation of DBIM for measurement configurations in which receiver positions are fixed was studied, and a novel regularization scheme was developed. The regularization parameter needed to stabilize the inversion process initially was found through the Rayleigh quotient iteration, then relaxed according to the relative residual error between the measured and estimated scattered fields. The DBIM was successfully stabilized for both full and partial receiver angular coverage without a significant loss in spatial resolution. The effects of variable density in the reconstructions were briefly explored through simulations. The ability to reconstruct targets with moderate contrast was validated through experimental measurements. Speed of sound profiles for balloons filled with saline in a background of water were reconstructed using multiple frequency DBIM techniques. The mean squared error for speed of sound reconstructions of the balloon phantoms with 16.4% sound speed contrast was 1.1%.     

Two-dimensional inverse problem of diffusion tomography: the approach applicable for small inclusions

A method for solving the two-dimensional inverse problems of optical diffusion tomography is proposed. The method is especially designed for the imaging of small inclusions embedded in the backgrounds of strongly scattering media. Numerical simulations show that the results are stable with respect to external noise at the boundary of the sample. The location of an inclusion is obtained with an accuracy of the order of several photon transport mean-free paths in the medium in cases both with and without noise in the scattering data used for the solution of the inverse problem.