Algorithms for suppressing ultrasonic backscattering from material structure

In pulse-echo ultrasonic inspection the backscattering from the material structure appears in the received ultrasonic images as clutter, often referred to as grain noise, which impairs the inspection results. A toolbox including algorithms for suppressing ultrasonic clutter is presented in the paper. Several processing algorithms capable of suppressing grain noise have been proposed, of which the split spectrum processing (SSP) probably is the most renowned. The classical SSP technique applies a filter bank to some frequency band that has to be precisely known in advance, to obtain a set of narrow-band signals that are tested for mutual correlation using some statistical operation. A number of SSP algorithms with different statistical operations are included in the toolbox. A completely different approach is to use explicit statistical models of grain noise and defects and to design an optimal filter based on those models. A simple such algorithm, based on noncoherent detection (NCD) known from communications, is also included in the toolbox. The toolbox, implemented in Matlab, is provided with a user-friendly graphical interface facilitating comparison of the algorithms.     

Quantitative estimation of ultrasonic attenuation in a solid in the immersion case with correction of diffraction effects

This paper presents a method of diffraction correction for the log-spectral difference method to estimate quantitatively attenuation of a solid in the immersion case. The correction method is established based on the angular spectrum approach that is used to calculate the echoes from the front and back surfaces of the immersed solid. An example is given of a copper plate submerged in water and inspected by a linear array with a cylindrically curved surface. The correction method is first applied to a theoretical estimation of attenuation which is linearly dependent on frequency. The results have shown that the evaluated attenuation coefficient is in excellent agreement with the exact value. Then the method is applied to a real situation, in which the results have shown that the method yields reasonable evaluated attenuation values. This work has demonstrated that the method is able to correct effectively the diffraction effect so as to achieve a quantitative estimation of attenuation.     

Nonlinear mapping technique for data visualization and clustering assessment of LIBS data: application to ChemCam data

ChemCam is a remote laser-induced breakdown spectroscopy (LIBS) instrument that will arrive on Mars in 2012, on-board the Mars Science Laboratory Rover. The LIBS technique is crucial to accurately identify samples and quantify elemental abundances at various distances from the rover. In this study, we compare different linear and nonlinear multivariate techniques to visualize and discriminate clusters in two dimensions (2D) from the data obtained with ChemCam. We have used principal components analysis (PCA) and independent components analysis (ICA) for the linear tools and compared them with the nonlinear Sammon’s map projection technique. We demonstrate that the Sammon’s map gives the best 2D representation of the data set, with optimization values from 2.8% to 4.3% (0% is a perfect representation), together with an entropy value of 0.81 for the purity of the clustering analysis. The linear 2D projections result in three (ICA) and five times (PCA) more stress, and their clustering purity is more than twice higher with entropy values about 1.8. We show that the Sammon’s map algorithm is faster and gives a slightly better representation of the data set if the initial conditions are taken from the ICA projection rather than the PCA projection. We conclude that the nonlinear Sammon’s map projection is the best technique for combining data visualization and clustering assessment of the ChemCam LIBS data in 2D. PCA and ICA projections on more dimensions would improve on these numbers at the cost of the intuitive interpretation of the 2D projection by a human operator.     

Maximum a posteriori deconvolution of ultrasonic signals using multiple transducers

A new method for deconvolution of ultrasonic pulse-echo measurements employing multiple-transducer setup is proposed in the paper. An optimal way of estimating the material reflection sequence for a linear signal generation model using maximum a posteriori estimation is proposed. The method combines the measurements from a number of transducers covering different frequency bands yielding an optimal estimate of the reflection sequence. The main idea of this approach is to complement the information unavailable from one transducer in some frequency bands with the information from the other transducers. The method is based on the assumption that the measurements are performed using transducers with identical apertures and apodization, which are located exactly at the same position relative to the test object during the measurement. An error analysis presented in the paper proves that when the above assumptions are fulfilled, the proposed method, by utilizing more data for estimation, consistently yields more accurate reflection sequence estimates than the classical Wiener filter. Experimental evidence is presented using both simulated and real ultrasonic data as a verification of the correctness of the multiple-transducer model and the estimation scheme. An illustration of the advantages of the method is also given using real ultrasonic data.     

A spatial impulse response based method for determining effective geometrical parameters for spherically focused transducers

This paper proposes a novel method for determining effective geometrical parameters (GPs) of spherically focused transducer. Theoretical basis of the method is the spatial impulse response method (SIRM) that is a powerful tool to calculate transient fields from a piston-like transducer. Experimentally, the method is implemented by measuring arrival times of the pulse signals (either transmission or pulse echo signals) and using them to determine the time intervals between the direct and edge wave pulses. Since the time intervals for a given measurement position depend on the transducer's GPs, then the effective GPs can be obtained. The measurements are carried out in the near-field so that the time intervals can be resolved. The proposed method has been applied to a spherically focused transducer operating in pulse echo mode. The results have shown that the effective GPs used in the SIRM give a closer agreement between the theory and the measurements than the nominal GPs. The method is best suited for characterizing broadband transducers with short pulse excitations. With certain modification it can be also applied to narrowband transducers. The method can be also applied to transducers of other shapes, e.g., planar circle, and planar or curved rectangle provided the piston approximation is appropriate.     

Anion concentration-dependent partitioning mechanism in the extraction of uranium into room-temperature ionic liquids

The mode of partitioning of uranyl ion between nitrate-containing aqueous phases and various N,N'-dialkylimidazolium-based room-temperature ionic liquids (RTILs) in the presence of tri-n-butyl phosphate (TBP) is shown to change from an ion-exchange process to one involving extraction of a neutral uranyl-TBP-nitrato complex as the aqueous nitrate concentration is increased. Increasing the hydrophobicity of the RTIL cation eventually leads to nitrato complex extraction as the predominant mode of partitioning, regardless of nitrate concentration.     

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.     

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.     

Identical extraction behavior and coordination of trivalent or hexavalent f-element cations using ionic liquid and molecular solvents

The extraction of both UO2(2+) and trivalent lanthanide and actinide ions (Am3+, Nd3+, Eu3+) by dialkylphosphoric or dialkylphosphinic acids from aqueous solutions into the ionic liquid, 1-decyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide has been studied and compared to extractions into dodecane. Radiotracer partitioning measurements show comparable patterns of distribution ratios for both the ionic liquid/aqueous and dodecane/aqueous systems, and the limiting slopes at low acidity indicate the partitioning of neutral complexes in both solvent systems. The metal ion coordination environment, elucidated from EXAFS and UV-visible spectroscopy measurements, is equivalent in the ionic liquid and dodecane solutions with coordination of the uranyl cation by two hydrogen-bonded extractant dimers, and of the trivalent cations by three extractant dimers. This is the first definitive report of a system where both the biphasic extraction equilibria and metal coordination environment are the same in an ionic liquid and a molecular organic solvent.     

Synthetic aperture imaging using sources with finite aperture: deconvolution of the spatial impulse response

A method for ultrasonic synthetic aperture imaging using finite-sized transducers is introduced that is based on a compact, linear, discrete model of the ultrasonic measurement system developed using matrix formalism. Using this model a time-domain algorithm for deconvolution of the transducer's spatial impulse responses (SIRs) is developed that is based on a minimum mean square error (MMSE) criterion. The algorithm takes the form of a spatiotemporal filter that compensates for the SIRs associated with a finite-sized transducer at every point of the processed image. A major advantage of the proposed method is that it can be used for any transducer, provided that its associated SIRs are known. This is in contrast to the synthetic aperture focusing technique (SAFT), which treats the transducer as a point source. The performance of the method is evaluated with simulations and experiments, performed in water using a linear phased array. The results obtained using the proposed method are compared to those obtained with a classical time-domain SAFT algorithm. For a finite aperture source, it is clearly shown that the resolution obtained using the proposed method is superior to that obtained using the SAFT algorithm.