Iterative two-step temporal phase-unwrapping applied to high sensitivity three-dimensional profilometry

Although temporal phase unwrapping method can be applied to solve some problems to measure an object with steep shapes, isolated parts or fringe undersampling in three-dimensional (3D) shape measurement, it needs to acquire and process a sequence of fringe pattern images which would take much time. Servin et al. proposed a 2-step temporal phase unwrapping algorithm, which only needs the 2 extreme phase-maps to achieve exactly the same results as standard temporal unwrapping method. But this method is only validated by computer simulation, shortage of experimental demonstration, its sensitivity coefficient G is limited, and it cannot be used when the G value is larger. We proposed an iterative two-step temporal phase-unwrapping algorithm which is an extension of Servin׳s method. First, add a fringe pattern with an intermediate sensitivity, project the fringe patterns of different sensitivity onto the tested object’s surface, and collect deformed fringe patterns with a CCD camera. Then we obtain the unwrapped phase with larger sensitivity coefficient G by cascading the sensitivity coefficients. And we derive the initial phase conditions of the 2-step temporal phase unwrapping algorithm. Finally, the experimental evaluation is conducted to prove the validity of the proposed method. The results are analyzed and compared with Servin׳s method. The experimental results show that the proposed method can achieve higher sensitivity and more accurate measurement, and it can overcome the main disadvantages encountered by Servin׳s method.     

Weighted robust Basis Function for phase unwrapping

This work presents a robust algorithm for phase unwrapping. The proposed algorithm is based on the expansion of the estimated phase through a linear combination of a set of Basis Functions. We present a novel weighted robust functional which is minimised using a two step strategy. This model allows us to reduce the influence of noise and to remove inconsistent pixels in the estimation of the unwrapped phase. The proposed model assumes that the phase is smooth. Under this assumption, experiments demonstrate that if the phase is corrupted by high levels of noise, our model presents a better performance than state of the art algorithms. For low levels of noise, the results are comparable.     

一种新的层次粒子滤波的目标跟踪方法

提出一种新的层次粒子滤波算法,选择局部区域特征点和颜色信息建立目标模型,引入粒子的二阶采样过程.算法通过粒子的一阶权重更新获得好的初始分布,二阶权重更新保证粒子的高置信度和高的采样效率,当粒子数目小于一定阈值时进行重要性重采样,利用仿射模型对目标区域精确定位及姿态修正.实验表明:改进算法将目标局部特征分布与目标颜色信息相结合,通过二阶采样过程,保证了局部特征跟踪的稳定性,解决了经典理论中误匹配导致的采样点发散问题,在目标部分遮挡情况下也可以完成实时目标跟踪.     

Local zone wise elastic and plastic properties of electron beam welded Ti–6Al–4V alloy using digital image correlation technique: A comparative study between uniform stress and virtual fields method

Joining of materials using welding results in the formation of material zones with varying microstructure across the weld. Extraction of the mechanical properties of those individual heterogeneous zones are important in designing components and structures comprised of welds. In this study, the zone wise local extraction of the elastic and plastic properties of an electron beam welded Ti–6Al–4V titanium alloy has been carried out using both the uniform stress method (USM) and the virtual fields method (VFM) involving digital image correlation (DIC) technique. The surface strain field obtained using DIC technique from a transverse weld specimen tensile testing is used for extracting the zone wise strain evolution. Initially, using uniform stress assumption, zone wise full range stress–strain curves are extracted. In USM methodology, the elastic and plastic material models are fitted to the zone wise stress–strain curves and required parameters are extracted from it. But inherent disadvantage is lot of images need to be processed for the parameter extraction. Recently, VFM is gaining lot of popularity in characterization domain as it is robust, accurate and faster. VFM is based on the principle of virtual work where, the weak form of local equilibrium equations and kinematically admissible virtual displacement fields are utilized for parameter extraction. Hollomon׳s power law is used here as the hardening rule. Young׳s modulus, Poisson׳s ratio, yield stress, strength coefficient and strain hardening exponent are the parameters extracted zone wise using both USM and VFM. A Vicker׳s microhardness measurement is also conducted across the weld zone towards mapping the strength behavior. Fusion zone has reported higher yield strength, strength coefficient and Poisson׳s ratio. Young׳s modulus value is found decreasing from base metal towards the fusion zone. The trend observed in parameter variation across the weld zone obtained by both USM and VFM compares very well. Due to various advantages associated with VFM technique it is generally recommended for parameter extraction.     

Acoustic detection and classification of river boats

We present a robust algorithm to detect the arrival of a boat of a certain type when other background noises are present. It is done via the analysis of its acoustic signature against an existing database of recorded and processed acoustic signals. We characterize the signals by the distribution of their energies among blocks of wavelet packet coefficients. To derive the acoustic signature of the boat of interest, we use the Best Discriminant Basis method. The decision is made by combining the answers from the Linear Discriminant Analysis (LDA) classifier and from the Classification and Regression Trees (CART) that is also accompanied with an additional unit, called Aisles, that reduces false alarms rate. The proposed algorithm is a generic solution for process control that is based on a learning phase (training) followed by an automatic real time detection while minimizing the false alarms rate.     

A linear least squares approach for evaluation of crack tip stress field parameters using DIC

In the present work, an experimental study is carried out to estimate the mixed-mode stress intensity factors (SIF) for different cracked specimen configurations using digital image correlation (DIC) technique. For the estimation of mixed-mode SIF׳s using DIC, a new algorithm is proposed for the extraction of crack tip location and coefficients in the multi-parameter displacement field equations. From those estimated coefficients, SIF could be extracted. The required displacement data surrounding the crack tip has been obtained using 2D-DIC technique. An open source 2D DIC software Ncorr is used for the displacement field extraction. The presented methodology has been used to extract mixed-mode SIF׳s for specimen configurations like single edge notch (SEN) specimen and centre slant crack (CSC) specimens made out of Al 2014-T6 alloy. The experimental results have been compared with the analytical values and they are found to be in good agreement, thereby confirming the accuracy of the algorithm being proposed.     

Parameter estimation of optical fringes with quadratic phase using the fractional Fourier transform

Optical fringes with a quadratic phase are often encountered in optical metrology. Parameter estimation of such fringes plays an important role in interferometric measurements. A novel method is proposed for accurate and direct parameter estimation using the fractional Fourier transform (FRFT), even in the presence of noise and obstacles. We take Newton׳s rings fringe patterns and electronic speckle pattern interferometry (ESPI) interferograms as classic examples of optical fringes that have a quadratic phase and present simulation and experimental results demonstrating the performance of the proposed method.     

Automated image processing algorithm for Young's fringe pattern analysis

An automatic of image processing for the analysis of Young's fringe pattern is proposed. A Young's fringe patttern is first FFT treated to get two patterns of 0 degree and 90 degree phase shifted. Regular two pattern shift algorithm is used to calculate the phase of the Young's pattern. One or two times of phase shift iteration are carried out to improve the phase. Then the phase is least square fitted to get a phase plane. The displacement magnitude and direction are eaasily induced from the phase plane. Because every step is full field treatment, very precise result is achieved. The algorithm is simple and easy to implement. Theory and experiment are presented to prove the proposal. A program is designed to execute the algorithm and the processing is automated by a Personal Computer with an image board.     

Hilbert reconstruction of phase-shifted second-harmonic holographic images

New techniques are presented that make phase-shifting holography viable for second-harmonic generation (SHG) holography with weak object fields. We developed an intrinsic phase shift calibration of SHG holograms, an algorithm that extracts the reference and object intensity directly from a set of phase-shifted holographic data, and a more robust phase-shifting holography reconstruction algorithm based on π-shifted hologram pairs that permits self-calibration of the phase shift and recovery of the complex field through a Hilbert transform.     

Revisiting the thermostatistics of the Grüneisen parameters and applications to quasiharmonic solids

We present an integrated thermostatistical study of the Grüneisen parameter and its application to a quasiharmonic solid. This analysis comprises two parts. The first part considers the connections between different thermodynamic formulations of the Grüneisen parameter, including those arising from the Mie-Grüneisen formula and related equations of state. We also establish the most general consequences of the so-called Grüneisen׳s rule. In the second part, Grüneisen׳s rule is used to establish the thermostatistics of a solid that obeys the quasiharmonic approximation. Both cases are relevant to the application of fundamental microscopic models of the heat capacity of solids.     

Kalman filter based initial guess estimation for digital image correlation

In applications digital image correlation based algorithms often present a basis for analysis of movement/deformation of bodies. The sequence of the obtained images is analyzed for this purpose. Especially, in cases when the body׳s movement/deformation between two successive images is significant, the initial guess can have a major influence on the execution speed of the algorithm. In the worst case it can even cause the divergence of the algorithm. This was the inspiration to develop a new and unique approach for an accurate and reliable determination of an initial guess for each image pixel. Kalman filter has been used for this purpose. It uses past measurements of observed variable(s) for calculations. Beside that it also incorporates state space model of the actual system. This is one of the most important advantages provided by Kalman filter. The determined initial guess by the proposed method is actually close to the true one and it enables fast convergence. Even more important property of this approach is the fact that it is not path-dependant because each image pixel, which is defined in ROI, is tracked through the sequence of images based on its own past measurements and general state space model. Consequently, the proposed method can be used to analyze tasks where discontinuities between image pixels are present. The applied method can be used to predict an initial guess where reference and deformed subsets are related by translational and rotational motion. The advantages mentioned above are verified with numerical and real experiments. The experimental validations are performed by NR (Newton–Raphson) approach which is the most widely used. Beside NR method the presented algorithm is applicable for other registration methods as well. It is used as an addition for calculation of initial guesses in a sequence of deformed images.     

Statistical hadronization and hadronic micro-canonical ensemble II

We present a Monte Carlo calculation of the micro-canonical ensemble of the ideal hadron-resonance gas including all known states up to a mass of about 1.8 GeV and full quantum statistics. The micro-canonical average multiplicities of the various hadron species are found to converge to the canonical ones for moderately low values of the total energy, around 8 GeV, thus bearing out previous analyses of hadronic multiplicities in the canonical ensemble. The main numerical computing method is an importance sampling Monte Carlo algorithm using the product of Poisson distributions to generate multi-hadronic channels. It is shown that the use of this multi-Poisson distribution allows for an efficient and fast computation of averages, which can be further improved in the limit of very large clusters. We have also studied the fitness of a previously proposed computing method, based on the Metropolis Monte Carlo algorithm, for event generation in the statistical hadronization model. We find that the use of the multi-Poisson distribution as proposal matrix dramatically improves the computation performance. However, due to the correlation of subsequent samples, this method proves to be generally less robust and effective than the importance sampling method.Received: 9 July 2004, Revised: 21 July 2004, Published online: 9 November 2004     

An algorithm for calculating the chemical potential in associating and reacting fluids

We present a Monte Carlo algorithm to determine the chemical potential of associating and reacting fluids. The algorithm is based on the fact that the chemical potential of a component is the same in the monomer (unbonded) state as in any bonded state. We demonstrate that the chemical potential of the unbonded specie can be calculated at relatively low cost and with high accuracy. The algorithm is applicable to both homogeneous as well as inhomogeneous systems. We compare the results of the presented algorithm against the findings of Widom's single stage particle insertion method for four commonly encountered inhomogeneous systems of associating fluids in phase equilibria studies. The constancy of the chemical potential throughout an inhomogeneous system under equilibrium is used as an independent test of the algorithm. The uncertainty in chemical potential values over the system for the cases studied was on an average 30 times smaller in the new algorithm, with at least 5 times fewer insertions than in the traditional Widom's method.     

Half-quadratic cost functions for phase unwrapping

We present a generic regularized formulation, based on robust half-quadratic regularization, for unwrapping noisy and discontinuous wrapped phase maps. Two cases are presented: the convex case and the nonconvex case. The unwrapped phase with the convex formulation is unique and robust to noise; however, the convex function solution deteriorates as a result of real discontinuities in phase maps. Therefore we also present a nonconvex formulation that, with a parameter continuation strategy, shows superior performance.     

基于选择采样的高效迭代相位提取算法

在光学精密测量中,相移干涉法应用广泛。常用的相移器件容易出现相移误差,采用等步距相位提取算法会产生测量误差。基于最小二乘的迭代相位提取算法可以有效消除该类相位提取误差,提高测量精度,但是其迭代过程运行时间长,效率低。提出了一种基于选择采样的迭代相位提取算法,先对干涉图像进行等间隔抽样,降低计算量;再根据对比度滤除干涉图像中低质量像素点,防止误差增大,进行最小二乘迭代求解相位。仿真实验对算法进行了分析和验证,在抽样间隔为2时的选择采样方法与所有像素点全部代入计算相比,运行时间从6.687 s降为0.725 s,均方根误差仅为0.032 9。实验结果证明:选择采样的迭代相位提取算法运算时间短、误差小,非常适合高速相移干涉测量应用。     

Improved differential 3D shape retrieval

Phase unwrapping is a complex step in three-dimensional (3D) surface measurement. To simplify the computation process, Martino et al. proposed a differential algorithm. However, it will result in large error when the orthogonal fringes are not in horizontal or vertical direction. To solve this problem, the relationship between projector׳s and camera׳s coordinate systems is introduced. With the data obtained from coordinate transformation, the improved differential algorithm can be used for orthogonal fringes in any direction. Besides that, taking advantage of Fourier differentiation theorem makes operation and calculation simpler. By contrast, the results of experiments show that the proposed method is applicable to the patterns with orthogonal fringes in every direction. In addition, Fourier differentiation theorem effectively increases the speed of differential process.     

Robust Quantum Computing in Decoherence-Free Subspaces with Double-Dot Spin Qubits

We present a scheme for implementing robust quantum gates in decoherence-free subspaces （DFSs） with double-dot spin qubits. Through the resonator-assisted interaction, the controllable interqubit couplings can be achieved only by adjusting the qubit transition frequencies. We construct a set of logic gates on the DFS-encoded qubits to eliminate the collective noise effects, and thus the gate fidelities can be enhanced remarkably. This proposal may offer a potential approach to realize the robust quantum computing with spin qubits.     

A practical computational framework for the multidimensional moment-constrained maximum entropy principle

The maximum entropy principle is a versatile tool for evaluating smooth approximations of probability density functions with a least bias beyond given constraints. In particular, the moment-based constraints are often a common prior information about a statistical state in various areas of science, including that of a forecast ensemble or a climate in atmospheric science. With that in mind, here we present a unified computational framework for an arbitrary number of phase space dimensions and moment constraints for both Shannon and relative entropies, together with a practical usable convex optimization algorithm based on the Newton method with an additional preconditioning and robust numerical integration routine. This optimization algorithm has already been used in three studies of predictability, and so far was found to be capable of producing reliable results in one- and two-dimensional phase spaces with moment constraints of up to order 4. The current work extensively references those earlier studies as practical examples of the applicability of the algorithm developed below.     

Gyrator变换全息图及其在图像加密中的应用

提出了gyrator变换全息图,利用gyrator变换快速算法模拟实现了gyrator变换全息图的产生和再现,并研究了基于相移数字全息的gyrator变换全息图.在此基础上提出了采用正弦相位光栅实现光学图像加密的新方法.该方法利用gyrator变换在相空间的旋转特性,将gyrator变换的角度、光栅的频率及光栅的旋转角度作为加密密钥,并利用两个或两个以上的gyrator变换系统的级联实现图像加密,增加了系统的安全性.依据相移数字全息进行的两个gyrator变换系统级联的仿真实验验证了该方法的可行性、有效性及其良好的安全性能.