A new method for the estimation of oceanic mixed-layer depth using shipboard X-band radar images

Shipboard X-band radar images acquired on 24 June 2009 are used to study nonlinear internal wave characteristics in the northeastern South China Sea. The studied images show three nonlinear internal waves in a packet. A method based on the Radon Transform technique is introduced to calculate internal wave parameters such as the direction of propagation and internal wave velocity from backscatter images. Assuming that the ocean is a two-layer finite depth system, we can derive the mixed-layer depth by applying the internal wave velocity to the mixed-layer depth formula. Results show reasonably good agreement with in-situ thermistor chain and conductivity-temperature-depth data sets.

一种快速获取激光拼焊焊缝错配的新算法

本文讨论了基于视觉和结构光的焊缝形貌视觉检测系统的组成原理,针对拼板激光焊接工程中对焊缝检测的实时性要求高的特点,创新性地提出了一种快速获得焊缝错配缺陷的算法,主要步骤为:首先基于有关标准的检测指标阈值的设定,再次,针对在线提取到的图像进行预处理,主要是加窗和中值滤波;最后为Radon变换与错配检测。该算法能减少计算任务,从而快速获得错配检测指标。实验给出了不等厚板拼焊时的线性错配的数值和分布,同时验证了该算法的有效性。

叠前地震数据规则化、重建及噪音压制

I address of reconstruction of spatial irregular sampling seismic data to regular grids. Spatial irregular sampling data impairs results of prestack migration, multiple attenuations, spectra estimation. Prestack 5-D volumes are often divided into sub-sections for further processing. Shot gathers are easy to obtain from irregular sampling volumes. My strategy for reconstruction is as follows: I resort irregular sampling gathers into a form of easy to bin and perform bin regularization, then utilize F-K inversion to reconstruct seismic data. In consideration of poor ability of F-K regularization to fill in large gaps, I sort regular sampling gathers to CMP and proposed high-resolution parabolic Radon transform to interpolate data and extrapolate offsets. To strong interfering noise--multiples, I use hybrid-domain high-resolution parabolic Radon transform to attenuate it. F-K regularization demand ultimately for lower computing costs. I proposed several methods to further improve efficiency of F-K inversion: first I introduce 1D and 2D NFFT algorithm for a rapid calculation of DFT operators; then develop fast 1D and 2D CG method to solve least-square equations, and utilize preconditioner to accelerate convergence of CG iterations; what’s more, I use Delaunay triangulation for weight calculation and use bandlimit frequency and varying bandwidth technique for competitive computation. Numerical 2D and 3D examples are offered to verify reasonable results and more efficiency. F-K regularization has poor ability to fill in large gaps, so I rearrange data as CMP gathers and develop hybrid-domain high-resolution parabolic Radon transforms which be used ether to interpolate null traces and extrapolate near and far offsets or suppress a strong interfere noise: multiples. I use it to attenuate multiples to verify performances of our algorithm and proposed routines for industrial application. Numerical examples and field data examples show a nice performance of our method.

叠前时间偏移：各向异性、多次波衰减及应用

At present, in order to image complex structures more accurately, the seismic migration methods has been developed from isotropic media to the anisotropic media. This dissertation develops a prestack time migration algorithm and application aspects for complex structures systematically. In transversely isotropic media with a vertical symmetry axis (VTI media), the dissertation starts from the theory that the prestack time migration is an approximation of the prestack depth migration, based on the one way wave equation and VTI time migration dispersion relation, by combining the stationary-phase theory gives a wave equation based VTI prestack time migration algorithm. Based on this algorithm, we can analytically obtain the travel time and amplitude expression in VTI media, as while conclude how the anisotropic parameter influence the time migration, and by analyzing the normal moveout of the far offset seismic data and lateral inhomogeneity of velocity, we can update the velocity model and estimate the anisotropic parameter model through the time migration. When anisotropic parameter is zero, this algorithm degenerates to the isotropic time migration algorithm naturally, so we can propose an isotopic processing procedure for imaging. This procedure may keep the main character of time migration such as high computational efficiency and velocity estimation through the migration, and, additionally, partially compensate the geometric divergence by adopting the deconvolution imaging condition of wave equation migration. Application of this algorithm to the complicated synthetic dataset and field data demonstrates the effectiveness of the approach. In the dissertation we also present an approach for estimating the velocity model and anisotropic parameter model. After analyzing the velocity and anisotropic parameter impaction on the time migration, and based on the normal moveout of the far offset seismic data and lateral inhomogeneity of velocity, through migration we can update the velocity model and estimate the anisotropic parameter model by combining the advantages of velocity analysis in isotropic media and anisotropic parameter estimation in VTI media. Testing on the synthetic and field data, demonstrates the method is effective and very steady. Massive synthetic dataset、2D sea dataset and 3D field datasets are used for VTI prestack time migration and compared to the stacked section after NMO and prestack isotropic time migration stacked section to demonstrate that VTI prestack time migration method in this paper can obtain better focusing and less positioning errors of complicated dip reflectors. When subsurface is more complex, primaries and multiples could not be separated in the Radon domain because they can no longer be described with simple functions (parabolic). We propose an attenuating multiple method in the image domain to resolve this problem. For a given velocity model，since time migration takes the complex structures wavefield propagation in to account, primaries and multiples have different offset-domain moveout discrepancies, then can be separated using techniques similar to the prior migration with Radon transform. Since every individual offset-domain common-reflection point gather incorporates complex 3D propagation effects, our method has the advantage of working with 3D data and complicated geology. Testing on synthetic and real data, we demonstrate the power of the method in discriminating between primaries and multiples after prestack time migration, and multiples can be attenuated in the image space considerably.

Fourier Transform Method for Measuring the Micro Amplitu of Water Surface Waves

采用一种非接触的光学方法傅立叶变换莫尔法(Fourier transform method),结合数字图像处理技术,对微幅振荡的水表面波的振幅进行测量.它是对全场中每一个像素点进行测量,比接触测量法具有更高的灵敏度.它为微幅水表面波振幅的测量提供了一种手段.通过将计算机生成的周期性光栅图像经投影机直接投影到被测物体的参考平面,经CCD摄像头、图像板捕捉存储形成数字化的光栅图像,利用傅立叶变换莫尔法处理光栅图像,从而获得包含有水表面波的振幅的相位信息,再经适当的几何变换获得振幅信息.我们在垂直振荡装置上进行了不同激励频率和不同振幅的表面波的振幅测量.

Using wavelet transform to study the Lipschitz local singular exponent in wall turbulence

In this paper, wavelet,transform is introduced to study the Lipschitz local singular exponent for characterising the local singularity behavior of fluctuating velocity in wall turbulence. I, is found that the local singular exponent is negative when the ejections and sweeps of coherent structures occur in a turbulent boundary layer.

A two-step phase-retrieval method in Fourier-transform ghost imaging

A two-step phase-retrieval method, based on Fourier-transform ghost imaging, was demonstrated. For the complex objects, the phase-retrieval process was divided into two steps: first got the complex object's amplitude from the Fourier-transform patterns of the squared object function, then combining with the Fourier-transform patterns of the object function to get the phase. The theoretical basis of this technique is outlined, and the experimental results are presented. (C) 2008 Elsevier B.V. All rights reserved.

Observing the fractional Fourier transform by free-space Fresnel diffraction

The fractional Fourier transform of an object can be observed in the free-space Fresnel diffraction pattern of the object. (C) 1997 Optical Society of America

Scaled fractional Fourier transform and its optical implementation

The scaled fractional Fourier transform is suggested and is implemented optically by one lens for different values of phi and output scale. In addition, physically it relates the FRT with the general lens transform-the optical diffraction between two asymmetrically positioned planes before and after a lens. (C) 1997 Optical Society of America.

Phase-space analysis of wavefront coding imaging systems

We explore the use of the Radon-Wigner transform, which is associated with the fractional Fourier transform of the pupil function, for determining the point spread function (PSF) of an incoherant defocused optical system. Then we introduce these phase-space tools to analyse the wavefront coding imaging system. It is shown that the shape of the PSF for such a system is highly invarient to the defocous-related aberrations except for a lateral shift. The optical transfer function of this system is also investigated briefly from a new understanding of ambiguity function.

Phase unwrapping algorithms, respectively, based on path-following and discrete cosine transform

This paper describes a path-following phase unwrapping algorithm and a phase unwrapping algorithm based on discrete cosine transform (DCT) which accelerates the Computation and suppresses the propagation of noise. Through analysis of fringe pattern with serious noises simulated in mathematic model, we make a contrast between path-following algorithm and DCT algorithm. The advantages and disadvantages or analytical fringe pattern are also given through comparison of two algorithms. Three-dimensional experimental results have been given to prove the validity of these algorithms. Despite DCT phase unwrapping technique robustness and speed in some cases, it cannot be unwrapping inconsistencies phase. The path-following algorithm can be used in automation analysis of fringe patterns with little influence of noise. (c) 2007 Elsevier GmbH. All rights reserved.

Measurement of gain spectrum for Fabry-Perot semiconductor lasers by the Fourier transform method with a deconvolution process

To improve the accuracy of measured gain spectra, which is usually limited by the resolution of the optical spectrum analyzer (OSA), a deconvolution process based on the measured spectrum of a narrow linewidth semiconductor laser is applied in the Fourier transform method. The numerical simulation shows that practical gain spectra can be resumed by the Fourier transform method with the deconvolution process. Taking the OSA resolution to be 0.06, 0.1, and 0.2 nm, the gain-reflectivity product spectra with the difference of about 2% are obtained for a 1550-nm semiconductor laser with the cavity length of 720 pm. The spectra obtained by the Fourier transform method without the deconvolution process and the Hakki-Paoli method are presented and compared. The simulation also shows that the Fourier transform method has less sensitivity to noise than the Hakki-Paoli method.