Semi-sparse deconvolution robust to uncertainties in the impulse responses

The received signal in ultrasonic pulse-echo inspection can be modeled as a convolution between an impulse response, or prototype, and the reflection sequence that is the impulse characteristic of the inspected object. Deconvolution aims at approximately inverting this process to improve the time resolution so that the overlap between echoes from closely spaced reflectors becomes small. For the relatively high contrast reflection sequences often found in non-destructive testing applications, semi-sparse deconvolution algorithms will typically yield better results than the classical Wiener filter solution. However, this requires a prototype that is a good representative for all echo responses found in the signals. Since, in practice, the prototype is often chosen by the operator directly from the inspection data, the prototype may incidentally be a bad representative for modeling the echoes for the remaining part of the object. Because of the sensitivity to deviations in the prototype this can yield deconvolution results with poor reproducibility. This paper presents a new semi-sparse deconvolution algorithm that is robust to deviations in the prototype. The new robust algorithm is based on a modification of an earlier presented non-robust semi-sparse algorithm. The robustness is obtained by including a stochastic model of the variations in the prototypes to the signal model when deriving the algorithm. Experiments performed using simulated data verify that the robust algorithm is less sensitive to deviations in the prototypes compared to the non-robust version of the algorithm and show that the proposed algorithm yields better estimates than its non-robust version and the Wiener filter in scenarios for which the algorithm was derived. Results using real ultrasonic data further show that the algorithm can be useful in practical scenarios where similar deconvolution results are required from slightly different echoes.     

Determining railgun plasma current distribution using Jansson'smethod to deconvolve B-dot probe signals

The authors present a method of determining the plasma armature current distribution from B-dot probe signals. This method utilizes the knowledge that the B-dot probe signal can be shown to be a convolution of the plasma armature current distribution with the impulse response function of the B-dot probe. Jansson's (1984) relaxation-based nonlinear deconvolution technique is developed for use with recorded B-dot probe data. Results obtained with this technique indicate that significant improvement in estimations can be obtained over previously used linear deconvolution methods     

Regularization for improving the deconvolution in real-time near-field acoustic holography

Near-field acoustic holography is a measuring process for locating and characterizing stationary sound sources from measurements made by a microphone array in the near-field of the acoustic source plane. A technique called real-time near-field acoustic holography (RT-NAH) has been introduced to extend this method in the case of nonstationary sources. This technique is based on a formulation which describes the propagation of time-dependent sound pressure signals on a forward plane using a convolution product with an impulse response in the time-wavenumber domain. Thus the backward propagation of the pressure field is obtained by deconvolution. Taking the evanescent waves into account in RT-NAH improves the spatial resolution of the solution but makes the deconvolution problem "ill-posed" and often yields inappropriate solutions. The purpose of this paper is to focus on solving this deconvolution problem. Two deconvolution methods are compared: one uses a singular value decomposition and a standard Tikhonov regularization and the other one is based on optimum Wiener filtering. A simulation involving monopoles driven by nonstationary signals demonstrates, by means of objective indicators, the accuracy of the time-dependent reconstructed sound field. The results highlight the advantage of using regularization and particularly in the presence of measurement noise.     

Wavelet based noise suppression technique and its application to ultrasonic flaw detection

The noise suppression techniques with wavelet transform (WT) are widely used in non-destructive testing and evaluation (NDT&E), especially in ultrasonics. Complete reconstruction theory with hard or soft thresholds, reconstruction technique based on the singularities of noise and signal, matched filter with an impulse response, and optimal frequency-to-bandwidth ratio of wavelet technique have all been used to analyze ultrasonic signals for noise suppression. But a more simple and effective technique has been pursued for decades. This paper develops a new technique using WT for the right purpose. In this work, WT is treated as a band-pass filter whose central frequency and frequency bandwidth (CF&FB) are determined by the spectra distribution of an ultrasonic signal captured from real testing situation. For the purpose of matching their CF&FB well, a technique for evaluating the optimal scale of a daughter wavelet is carried out too. By acting this daughter wavelet as a band-pass filter, we can obtain excellent de-noising results, even when the signal to noise ratio (SNR) is below -18 dB. The performance of the technique has been done by ultrasonic signals with computer generated white noises. Finally, the experimental verification is performed on a pipeline specimen with man-made small flaws with good results obtained. The results show that the technique is more suitable for processing heavy noised ultrasonic signals, and it can also be used in automatic flaw detection.     

Reflectivity function reconstruction using the frequency deconvolution technique

Deconvolution of ultrasonic echo signals improves resolution and quality of ultrasonic images. A frequency deconvolution algorithm depends on the Fast Fourier transform is proposed for ultrasonic data. The stability of the algorithm and the influence of the truncation effect on the deconvoluted results were investigated with respect to the duration time of reflectivity function reconstruction and the signal to noise ratio. Reliability of the separation of reconstructing the reflectivity of a biological tissue is estimated by frequency deconvolution of the echo ultrasound signals.     

Maximum a posteriori deconvolution of sparse ultrasonic signals using genetic optimization

Deconvolution of sparse spike sequences has received much attention in the field of seismic exploration. In certain situations in ultrasonic non-destructive testing (NDT) of materials, similar conditions as those found in seismic exploration occur. One example is the problem of detecting disbonds in layered aluminum structures. The reflection sequence convolved with the impulse response of the transducer results in masking closely spaced reflections. Deconvolution of these signals may reveal the reflection sequence and thus make the interpretation easier. In this paper we use the Bernoulli-Gaussian (BG) distribution for modeling the signal generation. This relatively simple model allows maximum a posteriori (MAP) estimation of the reflection sequence. A derivation of the MAP criterion is given for clarity. We propose a genetic algorithm for optimizing the MAP criterion. The genetic algorithm approach is motivated by the fact that the criterion is non-convex, implying that the criterion may have more than one local minimum point. The probability of obtaining the global optimal solution is increased by using the proposed genetic algorithm. One of the key features in genetic algorithms, the so-called cross-over operator, has been modified and adapted to the structure of the BG deconvolution problem to improve the efficiency of the search. The algorithm is tested on simulated data using the probability of detection (PD) and probability of false alarm (PFA) as evaluation criteria. The algorithm is also tested on real ultrasonic data from a layered aluminum structure. The results show considerable improvements in the possibility of interpreting the signals.     

超声探头的响应特性及暂态回波模型

建立检测系统的数学模型,可以更好地理解超声检测的物理本质。分析了超声波从产生、介质中传播、缺陷耦合以及接收的全过程,将缺陷回波表示为探头响应函数与缺陷响应的时域卷积。利用空间脉冲响应和基尔霍夫近似建立了超声波与平面型缺陷的耦合模型,用大平面试块底面回波和大平面响应进行反卷积求得了探头的响应函数,并详细分析了探头在不同偏置位置时不同大小缺陷响应的特点,发现缺陷回波由直达回波和边缘回波组成,直达回波和边缘回波极性相反,且直达回波的幅值远远大于边缘回波。      

Efficient array beam forming by spatial filtering for ultrasound B-mode imaging

This paper proposes an efficient array beam-forming method using spatial matched filtering (SMF) for ultrasonic imaging. In the proposed method, ultrasonic waves are transmitted from an array subaperture with fixed transmit focus as in conventional array imaging. At receive, radio frequency echo signals from each receive channel are passed through a spatial matched filter that is constructed based on the system transmit-receive spatial impulse response. The filtered echo signals are then summed without time delays. The filter concentrates and spatially registers the echo energy from each element so that the pulse-echo impulse response of the summed output is focused with acceptably low side lobes. Analytical beam pattern analysis and simulation results using a linear array show that this spatial filtering method can improve lateral resolution and contrast-to-noise ratio as compared with conventional dynamic receive focusing (DRF) methods. Experimental results with a linear array are consistent but point out the need to address additional practical issues. Spatial filtering is equivalent to synthetic aperture methods that dynamically focus on both transmit and receive throughout the field of view. In one common example of phase aberrations, the SMF method was degraded to a degree comparable to conventional DRF methods.     

Method for measuring violin sound radiation based on bowed glissandi and its application to sound synthesis

This work presents a method for measuring and computing violin-body directional frequency responses, which are used for violin sound synthesis. The approach is based on a frame-weighted deconvolution of excitation and response signals. The excitation, consisting of bowed glissandi, is measured with piezoelectric transducers built into the bridge. Radiation responses are recorded in an anechoic chamber with multiple microphones placed at different angles around the violin. The proposed deconvolution algorithm computes impulse responses that, when convolved with any source signal (captured with the same transducer), produce a highly realistic violin sound very similar to that of a microphone recording. The use of motion sensors allows for tracking violin movements. Combining this information with the directional responses and using a dynamic convolution algorithm, helps to improve the listening experience by incorporating the violinist motion effect in stereo.     

Nonlinear dynamic range compression deconvolution

We introduce a dynamic range image compression technique for nonlinear deconvolution; the impulse response of the distortion function and the noisy distorted image are jointly transformed to pump a clean reference beam in a two-beam coupling arrangement. The Fourier transform of the pumped reference beam contains the deconvolved image and its conjugate. In contrast to standard deconvolution approaches, for which noise can be a limiting factor in the performance, this approach allows the retrieval of distorted signals embedded in a very high-noise environment.     

Detecting small flaws near the interface in pulse-echo

Flaw detection near the interface surface is a common problem in many pulse-echo NDT applications due to interference with the interface echo, orders of magnitude above the flaw echoes. Several digital signal processing techniques like deconvolution, Hilbert transform and cepstrum analysis have been proposed to improve axial resolution. However, they require strict linearity, which takes a large portion of the system dynamic range just to fit the interface echo, thus reducing the dynamic range available for flaw detection. This work presents a new alternative based on the time-domain phase analysis of the received signals. Differently from conventional approaches, it works quite well with saturated signals resulting when a high gain is applied to detect small flaws. These can be detected in a range of a fraction of one wavelength from the interface surface, even using narrow-band transducers, as it has been experimentally verified. The method can be easily hardware implemented for real-time processing.     

Minimum entropy deconvolution of pulse-echo signals acquired from attenuative layered media.

In this article deconvolution of ultrasonic pulse-echo data acquired from attenuative layered media is considered. The problem is divided in two subproblems: treating the sparse reflection sequence caused by the layered structure of the media and treating the frequency-dependent attenuation. The first subproblem is solved by means of joint maximum a posteriori estimation of the assumed zero mean, white, nonstationary reflection sequence and its corresponding sequence of unknown standard deviations. This approach leads to an algorithm that seeks minimum entropy solutions for the reflection sequence and therefore the algorithm serves as a novel link between the classical Wiener filter and methods for sparse or minimum entropy deconvolution. The second subproblem is solved by introducing a new signal processing-oriented, linear discrete-time model for frequency-dependent attenuation in isotropic and homogeneous media. The deconvolution algorithm is tested using simulated data and its performance for real normal incidence pulse-echo data from a composite material is also demonstrated. The results show that the algorithm, in combination with the attenuation model, yields estimates that reveal the internal structure of the composite and, thus, simplify the interpretation of the ultrasonic data.     

An alternative approach for measuring the scattered acoustic pressure field of immersed single and multiple cylinders

The form function of an elastic target can be obtained from the scattered signal through a deconvolution process. The deconvolution process uses the signal measured from an acoustically hard target (reference signal) to compensate for the impulse response of the measurement system. In this paper, it is shown that this approach limits the usable frequency range of the signal and leads to inaccuracies in the final results. An alternative approach is proposed in which the reference signal is replaced by the specular echo. A procedure is described for extracting the specular echo from the measured signal even in cases where it is not completely isolated from the resonant components. Modifications are made to the existing deconvolution formulation and it is further extended to be applicable to multiple scattering measurements. Experimental results show that the new approach provides improved accuracy and wider usable frequency range in both single and multiple scattering experiments.     

基于Gabor时频滤波的电磁超声螺栓轴力测量*

高强度螺栓轴向预紧力(轴力)的测量在工程应用中具有重要意义。使用电磁超声波对螺栓轴力检测时对超声回波声时测量精度要求较高,传统互相关估计法对超声回波声时估计易因噪声干扰发生估计错误,无法满足轴力测量精度要求。本文针对互相关法对电磁超声测量信号声时估计存在不准确的问题,提出了Gabor时频滤波法。通过螺栓轴力测量实验采集测量信号,对测量信号进行Gabor变换,在时频域中进行滤波,再对去噪后的信号进行互相关估计测得信号的声时,进而计算出螺栓轴力。实验表明：Gabor时频滤波法能有效地滤除电磁超声信号中的噪声,改善互相关估计的稳定性,提高螺栓轴力测量的准确率。     

激光诱导叶绿素荧光寿命的测量及其特性分析

提出了一种用于评估植物生长状况及环境监测的激光诱导叶绿素荧光寿命测量方法. 采用波长355 nm的激光作为光源激发叶绿素荧光, 由光电倍增管接收其荧光信号, 由于被测叶绿素荧光衰减函数与激光脉冲、仪器响应函数卷积在一起, 根据它们的特性, 运用时间分辨测量法分别测得叶绿素荧光及其背景信号, 并结合一种新型解卷积算法可分离出真实的叶绿素荧光衰减函数, 从而获取叶绿素的荧光寿命. 测试结果表明: 该方法能够实现叶绿素荧光寿命的高精度实时监测, 对不同叶绿素含量的溶液荧光寿命进行了测试, 证明叶绿素含量与其荧光寿命具有相关性, 并且拟合了叶绿素含量与荧光寿命的标定曲线. 关键词： 荧光寿命 激光诱导荧光 时间分辨测量法 叶绿素含量     

Piezoelectric thick film ultrasonic transducers fabricated by a sol-gel spray technique

Thick film broadband ultrasonic transducers (UTs) produced by a sol-gel spray technique and operated below 10 MHz are presented. These UTs are formed by dispersing PZT and LiTaO3 particles, respectively in Al2O3 and PZT sol-gel solution. The 50-100 microm thick films have been deposited on curved steel, flat steel and aluminum substrates and steel rods. Ultrasonic pulse-echo signals with a signal to noise ratio of more than 25 dB are experimentally obtained for the operating temperatures up to 250 degrees C.     

模糊条码图像的正则化复原算法

由条码扫描仪获得条码图像的过程可以用理想条码信号与扫描仪光学系统点扩散函数的卷积模型来描述。反卷积是消除由光学系统点扩散带来的模糊现象的最好办法。为克服反卷积的病态问题,研究了反卷积的正则化方法;针对条码信号的特点,构建了适合于条码信号复原的惩罚项,提出了条码信号的正则化复原算法及其适合于计算机运算的迭代算法。通过实验研究了算法在不同情况下的抗干扰能力。实验结果表明,正则化条码信号复原算法在消除系统点扩散函数的影响的同时能够很好地抑制噪声。     

A hybrid MLS technique for room impulse response estimation

The measurement of room impulse response (RIR) when there are high background noise levels frequently means one must deal with very low signal-to-noise ratios (SNR). If such is the case, the measurement might yield unreliable results, even when synchronous averaging techniques are used. Furthermore, if there are non-linearities in the apparatus or system time variances, the final SNR can be severely degraded. The test signals used in RIR measurement are often disturbed by non-stationary ambient noise components. A novel approach based on the energy analysis of ambient noise - both in the time and in frequency - was considered. A modified maximum length sequence (MLS) measurement technique, referred to herein as the hybrid MLS technique, was developed for use in room acoustics. The technique consists of reducing the noise energy of the captured sequences before applying the averaging technique in order to improve the overall SNRs and frequency response accuracy. Experiments were conducted under real conditions with different types of underlying ambient noises. Results are shown and discussed. Advantages and disadvantages of the hybrid MLS technique over standard MLS technique are evaluated and discussed. Our findings show that the new technique leads to a significant increase in the overall SNR.     

A new technique for the identification of ultrasonic flaw signals using deconvolution

The identification of ultrasonic flaw signals is a difficult task in the angle beam ultrasonic testing of welded joints due to the presence of non-relevant signals from the geometric reflectors such as weld-roots and counter-bores. This paper describes a new approach to identify ultrasonic flaw signals in such a problematic situation. A similarity function is defined as the deconvolution of a target signal by a reference signal. The similarity functions for the same type of flaws/references are symmetric bandlimited impulse-like patterns with larger amplitudes while those for different types of flaws/references are asymmetrical broad patterns with relatively smaller amplitudes. Therefore, ultrasonic signals could be identified by the pattern of the similarity function. In the initial experiments, the proposed technique showed great potential for identifying ultrasonic flaw signals in the inspection of weld joints.