A 3D shape retrieval method for orthogonal fringe projection based on a combination of variational image decomposition and variational mode decomposition

The orthogonal fringe projection technique has as wide as long practical application nowadays. In this paper, we propose a 3D shape retrieval method for orthogonal composite fringe projection based on a combination of variational image decomposition (VID) and variational mode decomposition (VMD). We propose a new image decomposition model to extract the orthogonal fringe. Then we introduce the VMD method to separate the horizontal and vertical fringe from the orthogonal fringe. Lastly, the 3D shape information is obtained by the differential 3D shape retrieval method (D3D). We test the proposed method on a simulated pattern and two actual objects with edges or abrupt changes in height, and compare with the recent, related and advanced differential 3D shape retrieval method (D3D) in terms of both quantitative evaluation and visual quality. The experimental results have demonstrated the validity of the proposed method.     

An improved windowed Fourier transform for fringe demodulation

A modified algorithm of windowed Fourier transform (WFT) for phase retrieval from electronic speckle-shearing fringe patterns with carriers is proposed. The algorithm is based on the introduction of a fast Fourier transform (FFT) in WFT to reduce computation time for fringe demodulation. Since boundary effects in FFT will influence the accuracy of phase retrieval, the Gerchberg method is employed to extrapolate the fringe pattern at the boundaries to reduce boundary effects. A theoretical analysis of the algorithm is presented. Both simulated and experimental results show that the proposed method has reduced the computation time significantly compared with the convolution method of WFT without sacrificing measurement accuracy.     

Iterative phase retrieval without support

An iterative phase retrieval method for nonperiodic objects has been developed from the charge-flipping algorithm proposed in crystallography. A combination of the hybrid input-output (HIO) algorithm and the flipping algorithm has greatly improved performance. In this combined algorithm the flipping algorithm serves to find the support (object boundary) dynamically, and the HIO part improves convergence and moves the algorithm out of local minima. It starts with a single intensity measurement in the Fourier domain and does not require a priori knowledge of the support in the image domain. This method is suitable for general image recovery from oversampled diffuse elastic x-ray and electron-diffraction intensities. The relationship between this algorithm and the output-output algorithm is elucidated.     

一种基于傅里叶变换的分析载波条纹的新方法

针对传统傅里叶变换法处理光载波干涉条纹图时会有边缘效应产生的问题,提出了一种基于傅里叶变换法的外推延拓方法,并从理论上进行了数学推导。为了验证这种方法的正确性,分别对一维数字信号和二维空间载波条纹图进行了数值模拟,进一步分析了误差产生的原因,并与传统的傅里叶变换法对比。结果表明该法可以有效抑制传统傅里叶变换法处理光载波干涉条纹图时边缘效应所造成的较大误差,在基于空间域相位调制技术的波面干涉测量中,对空间载波条纹图进行处理,可以使相位的计算精度达到3.3 mrad。     

离焦模糊图像中汉字识别技术研究

提出了辐射投影变换(RPT),给出了数学模型和实现方法以及相应的投影函数(RPTF)。为了进一步提高其性能,对原始投影变换进行了优化改进,即累加原始函数的值并作旋转变换,从而得到了改进投影累加函数(IPAF),它可直接用于离焦模糊图像中汉字的识别。将傅里叶变换应用到IPAF上,使描述和识别用的特征量大为减少。采用基于IPAF的傅里叶子(前8个傅里叶系数)作为特征量,对离焦模糊的一级汉字进行识别,总识别率达到了97 4%。     

Asymmetric multiple-image hiding using phase retrieval technique based on amplitude- and phase-truncation in fractional Fourier domain

We propose a multiple-image hiding scheme based on the amplitude- and phase-truncation approach, and phase retrieval iterative algorithm in the fractional Fourier domain. The proposed scheme offers multiple levels of security with asymmetric keys. Multiple input images multiplied with random phase masks are independently fractional Fourier transformed with different orders. The individual keys and common keys are generated by using phase and amplitude truncation of fractional spectrum. After using two fractional Fourier transform, the resultant encrypted image is hided in a host image with phase retrieval iterative algorithm. Using the correct universal keys, individual keys, and fractional orders, one can recover the original image successfully. Computer simulation results with four gray-scale images support the proposed method. To measure the validity of the scheme, we calculated the mean square error between the original and the decrypted images. In this scheme, the encryption process and generation of decryption keys are complicated and should be realized using computer. For decryption, an optoelectronic setup has been suggested.     

Acoustic characteristics analysis of slender cavity system with arbitrary impedance end conditions

A modeling method for the dynamic characteristics analysis of a slender acoustical cavity with impedance end conditions is established. In order to satisfy the continuity requirement at impedance ends for the first order differential of sound pressure, field function is constructed as the standard Fourier series supplemented by boundary smoothed auxiliary polynomials. System characteristic equation is derived by solving the governing differential equation and impedance acoustic boundary of slender acoustical cavity system simultaneously,relevant acoustical modal information is obtained via the state space solution procedure. In numerical simulation, various acoustic variables, such as acoustical modal frequency, sound pressure modal shape, sound pressure response and the particle velocity, are presented for the slender acoustical cavity system with different boundary conditions and compared with those results in the existing literature. The correctness and effectiveness of the proposed method are then fully validated.     

考虑任意阻抗壁面条件管腔结构声场特性分析

针对任意阻抗壁面条件一维管腔声学系统建模,对系统动力学特性进行预报。为了满足阻抗边界条件对声压一阶导数连续性要求,管腔声压函数通过在标准傅里叶级数端点位置引入边界光滑辅助多项式进行构建。结合壁面阻抗声学边界和管腔声学Helmholtz控制微分方程得到强形式标准特征值问题,获得相应的声学模态信息。在数值仿真中,通过算例给出各种边界条件下管腔声学模态频率、声压振型、声压和质点振速频率响应曲线,与现有文献中相关结果进行对比,充分验证了本文求解方法的正确性和有效性,证明该方法可对任意阻抗壁面条件管腔系统声学特性进行准确预报。     

Exact boundary condition for time-dependent wave equation based on boundary integral

An exact non-reflecting boundary conditions based on a boundary integral equation or a modified Kirchhoff-type formula is derived for exterior three-dimensional wave equations. The Kirchhoff-type non-reflecting boundary condition is originally proposed by L. Ting and M.J. Miksis [J. Acoust. Soc. Am. 80 (1986) 1825] and numerically tested by D. Givoli and D. Cohen [J. Comput. Phys. 117 (1995) 102] for a spherically symmetric problem. The computational advantage of Ting–Miksis boundary condition is that its temporal non-locality is limited to a fixed amount of past information. However, a long-time instability is exhibited in testing numerical solutions by using a standard non-dissipative finite-difference scheme. The main purpose of this work is to present a new exact boundary condition and to eliminate the long-time instability. The proposed exact boundary condition can be considered as a limit case of Ting–Miksis boundary condition when the two artificial boundaries used in their method approach each other. Our boundary condition is actually a boundary integral equation on a single artificial boundary for wave equations, which is to be solved in conjunction with the interior wave equation. The new boundary condition needs only one artificial boundary, which can be of any shape, i.e., sphere, cubic surface, etc. It keeps all merits of the original Kirchhoff boundary condition such as restricting the temporal non-locality, free of numerical evaluation of any special functions and so on. Numerical approximation to the artificial boundary condition on cubic surface is derived and three-dimensional numerical tests are carried out on the cubic computational domain.     

A Fourier spectral method for the Navier–Stokes equations with volume penalization for moving solid obstacles

An efficient numerical scheme to compute flows past rigid solid bodies moving through viscous incompressible fluid is presented. Solid obstacles of arbitrary shape are taken into account using the volume penalization method to impose no-slip boundary condition. The 2D Navier–Stokes equations, written in the vorticity-streamfunction formulation, are discretized using a Fourier pseudo-spectral scheme. Four different time discretization schemes of the penalization term are proposed and compared. The originality of the present work lies in the implementation of time-dependent penalization, which makes the above method capable of solving problems where the obstacle follows an arbitrary motion. Fluid–solid coupling for freely falling bodies is also implemented. The numerical method is validated for different test cases: the flow past a cylinder, Couette flow between rotating cylinders, sedimentation of a cylinder and a falling leaf with elliptical shape.     

一种简化的分析光纤布拉格光栅的时域有限差分束传输法

提出了一种便于分析光纤布拉格光栅(FBG)、基于强制边界条件的简化TD-FD-BPM模型.采用交替隐式法(ADIM)和综合道格拉斯(GD)法,离散化时域束传输方程,获得比普通差分格式低二阶的误差。提出并采用强制边界条件,简化了光波导分析中边界处理问题.重点讨论和确定了在此条件下激励源位置、脉冲形状和径向分布函数、连续波激励条件等.在修正的正弦脉冲激励情况下,通过离散傅利叶变换,直接得到较为理想的短FBG反射谱.此法可减少数值计算时间,直观反映光脉冲传输.     

Implementation of nonreflecting boundary conditions for the nonlinear Euler equations

Computationally efficient nonreflecting boundary conditions are derived for the Euler equations with acoustic, entropic and vortical inflow disturbances. The formulation linearizes the Euler equations near the inlet/outlet boundaries and expands the solution in terms of Fourier–Bessel modes. This leads to an ‘exact’ nonreflecting boundary condition, local in space but nonlocal in time, for each Fourier–Bessel mode of the perturbation pressure. The perturbation velocity and density are then calculated using acoustic, entropic and vortical mode splitting. Extension of the boundary conditions to nonuniform swirling flows is presented for the narrow annulus limit which is relevant to many aeroacoustic problems. The boundary conditions are implemented for the nonlinear Euler equations which are solved in space using the finite volume approximation and integrated in time using a MacCormack scheme. Two test problems are carried out: propagation of acoustic waves in an annular duct and the scattering of a vortical wave by a cascade. Comparison between the present exact conditions and commonly used approximate local boundary conditions is made. Results show that, unlike the local boundary conditions whose accuracy depends on the group velocity of the scattered waves, the present conditions give accurate solutions for a range of problems that have a wide array of group velocities. Results also show that this approach leads to a significant savings in computational time and memory by obviating the need to store the pressure field and calculate the nonlocal convolution integral at each point in the inlet and exit boundaries.