Blind Deconvolution of Ultrasonic Signals Using High-Order Spectral Analysis and Wavelets

Defect detection by ultrasonic method is limited by the pulse width. Resolution can be improved through a deconvolution process with a priori information of the pulse or by its estimation. In this paper a regularization of the Wiener filter using wavelet shrinkage is presented for the estimation of the reflectivity function. The final result shows an improved signal to noise ratio with better axial resolution.Comment: 8 pages, CIARP 2005, LNCS 377

Improvement of signal analysis for the ultrasonic microscopy

This dissertation describes the improvement of signal analysis in ultrasonic microscopy for nondestructive testing. Specimens with many thin layers, like modern electronic components, pose a particular challenge for identifying and localizing defects. In this thesis, new evaluation algorithms have been developed which enable analysis of highly complex layer-stacks. This is achieved by a specific evaluation of multiple reflections, a newly developed iterative reconstruction and deconvolution algorithm, and the use of classification algorithms with a highly optimized simulation algorithm. Deep delaminations inside a 19-layer component can now not only be detected, but also localized. The new analysis methods also enable precise determination of elastic material parameters, sound velocities, thicknesses, and densities of multiple layers. The highly improved precision of determined reflections parameters with deconvolution also provides better and more conclusive results with common analysis methods.:Kurzfassung......................................................................................................................II Abstract.............................................................................................................................V List ob abbreviations........................................................................................................X 1 Introduction.......................................................................................................................1 1.1 Motivation.....................................................................................................................2 1.2 System theoretical description.....................................................................................3 1.3 Structure of the thesis..................................................................................................6 2 Sound field.........................................................................................................................8 2.1 Sound field measurement............................................................................................8 2.2 Sound field modeling..................................................................................................11 2.2.1 Reflection and transmission coefficients.........................................................11 2.2.2 Sound field modeling with plane waves..........................................................13 2.2.3 Generalized sound field position.....................................................................19 2.3 Receiving transducer signal.......................................................................................20 2.3.1 Calculation of the transducer signal from the sound field...............................20 2.3.2 Received signal amplitude..............................................................................21 2.3.3 Measurement of reference signals..................................................................24 3 Ultrasonic Simulation......................................................................................................27 3.1 State of the art............................................................................................................27 3.2 Simulation approach..................................................................................................28 3.2.1 Sound field measurement based simulation...................................................28 3.2.2 Reference signal based simulation.................................................................30 3.3 Determination of the impulse response.....................................................................31 3.3.1 1D ray-trace algorithm....................................................................................31 3.3.2 2D ray-trace algorithm....................................................................................33 3.3.3 Complexity reduction – optimizations.............................................................35 4 Deconvolution – Determination of reflection parameters............................................38 4.1 State of the art............................................................................................................39 4.1.1 Decomposition techniques..............................................................................39 4.1.2 Deconvolution.................................................................................................41 4.2 Analytic signal investigations for deconvolution.........................................................42 4.3 Single reference pulse deconvolution........................................................................44 4.4 Multi-pulse deconvolution..........................................................................................47 4.4.1 Homogeneous multi-pulse deconvolution.......................................................48 4.4.2 Multi-pulse deconvolution with simulated GSP profile....................................49 5 Reconstruction.................................................................................................................50 5.1 State of the art............................................................................................................50 5.2 Reconstruction approach...........................................................................................51 5.3 Direct material parameter estimation.........................................................................52 5.3.1 Sound velocities and layer thickness..............................................................52 5.3.2 Density, elastic modules and acoustic attenuation.........................................54 5.4 Iterative material parameter determination of a single layer......................................56 5.5 Reconstruction of complex specimens......................................................................60 5.5.1 Material characterization of multiple layers ....................................................60 5.5.2 Iterative simulation parameter optimization with correlation...........................62 5.5.3 Pattern recognition reconstruction of specimens with known base structure. 66 6 Applications and results.................................................................................................71 6.1 Analysis of stacked components................................................................................71 6.2 Time-of-flight and material analysis...........................................................................74 7 Conclusions and perspectives.......................................................................................78 References.......................................................................................................................82 Figures.............................................................................................................................86 Tables...............................................................................................................................88 Appendix..........................................................................................................................89 Acknowledgments.........................................................................................................100 Danksagung...................................................................................................................101Die vorgelegte Dissertation befasst sich mit der Verbesserung der Signalauswertung für die Ultraschallmikroskopie in der zerstörungsfreien Prüfung. Insbesondere bei Proben mit vielen dünnen Schichten, wie bei modernen Halbleiterbauelementen, ist das Auffinden und die Bestimmung der Lage von Fehlstellen eine große Herausforderung. In dieser Arbeit wurden neue Auswertealgorithmen entwickelt, die eine Analyse hochkomplexer Schichtabfolgen ermöglichen. Erreicht wird dies durch die gezielte Auswertung von Mehrfachreflexionen, einen neu entwickelten iterativen Rekonstruktions- und Entfaltungsalgorithmus und die Nutzung von Klassifikationsalgorithmen im Zusammenspiel mit einem hoch optimierten neu entwickelten Simulationsalgorithmus. Dadurch ist es erstmals möglich, tief liegende Delaminationen in einem 19-schichtigem Halbleiterbauelement nicht nur zu detektieren, sondern auch zu lokalisieren. Die neuen Analysemethoden ermöglichen des Weiteren eine genaue Bestimmung von elastischen Materialparametern, Schallgeschwindigkeiten, Dicken und Dichten mehrschichtiger Proben. Durch die stark verbesserte Genauigkeit der Reflexionsparameterbestimmung mittels Signalentfaltung lassen sich auch mit klassischen Analysemethoden deutlich bessere und aussagekräftigere Ergebnisse erzielen. Aus den Erkenntnissen dieser Dissertation wurde ein Ultraschall-Analyseprogramm entwickelt, das diese komplexen Funktionen auf einer gut bedienbaren Oberfläche bereitstellt und bereits praktisch genutzt wird.:Kurzfassung......................................................................................................................II Abstract.............................................................................................................................V List ob abbreviations........................................................................................................X 1 Introduction.......................................................................................................................1 1.1 Motivation.....................................................................................................................2 1.2 System theoretical description.....................................................................................3 1.3 Structure of the thesis..................................................................................................6 2 Sound field.........................................................................................................................8 2.1 Sound field measurement............................................................................................8 2.2 Sound field modeling..................................................................................................11 2.2.1 Reflection and transmission coefficients.........................................................11 2.2.2 Sound field modeling with plane waves..........................................................13 2.2.3 Generalized sound field position.....................................................................19 2.3 Receiving transducer signal.......................................................................................20 2.3.1 Calculation of the transducer signal from the sound field...............................20 2.3.2 Received signal amplitude..............................................................................21 2.3.3 Measurement of reference signals..................................................................24 3 Ultrasonic Simulation......................................................................................................27 3.1 State of the art............................................................................................................27 3.2 Simulation approach..................................................................................................28 3.2.1 Sound field measurement based simulation...................................................28 3.2.2 Reference signal based simulation.................................................................30 3.3 Determination of the impulse response.....................................................................31 3.3.1 1D ray-trace algorithm....................................................................................31 3.3.2 2D ray-trace algorithm....................................................................................33 3.3.3 Complexity reduction – optimizations.............................................................35 4 Deconvolution – Determination of reflection parameters............................................38 4.1 State of the art............................................................................................................39 4.1.1 Decomposition techniques..............................................................................39 4.1.2 Deconvolution.................................................................................................41 4.2 Analytic signal investigations for deconvolution.........................................................42 4.3 Single reference pulse deconvolution........................................................................44 4.4 Multi-pulse deconvolution..........................................................................................47 4.4.1 Homogeneous multi-pulse deconvolution.......................................................48 4.4.2 Multi-pulse deconvolution with simulated GSP profile....................................49 5 Reconstruction.................................................................................................................50 5.1 State of the art............................................................................................................50 5.2 Reconstruction approach...........................................................................................51 5.3 Direct material parameter estimation.........................................................................52 5.3.1 Sound velocities and layer thickness..............................................................52 5.3.2 Density, elastic modules and acoustic attenuation.........................................54 5.4 Iterative material parameter determination of a single layer......................................56 5.5 Reconstruction of complex specimens......................................................................60 5.5.1 Material characterization of multiple layers ....................................................60 5.5.2 Iterative simulation parameter optimization with correlation...........................62 5.5.3 Pattern recognition reconstruction of specimens with known base structure. 66 6 Applications and results.................................................................................................71 6.1 Analysis of stacked components................................................................................71 6.2 Time-of-flight and material analysis...........................................................................74 7 Conclusions and perspectives.......................................................................................78 References.......................................................................................................................82 Figures.............................................................................................................................86 Tables...............................................................................................................................88 Appendix..........................................................................................................................89 Acknowledgments.........................................................................................................100 Danksagung...................................................................................................................10

Signal processing for guided wave structural health monitoring

The importance of Structural Health Monitoring (SHM) in several industrial fields has been continuously growing in the last few years with the increasing need for the development of systems able to monitor continuously the integrity of complex structures. In order to be competitive with conventional non destructive evaluation techniques, SHM must be able to effectively detect the occurrence of damage in the structure, giving information regarding the damage location. Ultrasonic guided waves offer the possibility of inspecting large areas of structures from a small number of sensor positions. However, inspection of complex structures is difficult as the reflections from different features overlap. Therefore damage detection becomes an extremely challenging problem and robust signal processing is required in order to resolve strongly overlapping echoes. In our work we have considered at first the possibility of employing a deconvolution approach for enhancing the resolution of ultrasonic time traces and the potential and the limitations of this approach for reliable SHM applications have been shown. The effects of noise on the bandwidth of the typical signals in SHM and the effects of frequency dependent phase shifts are the main detrimental issues that strongly reduce the performance of deconvolution in SHM applications. The second part of this thesis is concerned with the evaluation of a subtraction approach for SHM when changes of environmental conditions are taken into account. Temperature changes result in imperfect subtraction even for an undamaged structure, since temperature changes modify the mechanical properties of the material and therefore the velocity of propagation of ultrasonic guided waves. Compensation techniques have previously been used effectively to overcome temperature effects, in order to reduce the residual in the subtraction. In this work the performance of temperature compensation techniques has been evaluated also in the presence of other detrimental effects, such as liquid loading and different temperature responses of materials in adhesive joints. Numerical simulations and experiments have been conducted and it has been shown that temperature compensation techniques can cope in principle with non temperature effects. It is concluded that subtraction approach represents a promising method for reliable Structural Health Monitoring. Nonetheless the feasibility of a subtraction approach for SHM depends on environmental conditions

医用超音波における散乱体分布の高解像かつ高感度な画像化に関する研究

Ultrasound imaging as an effective method is widely used in medical diagnosis andNDT (non-destructive testing). In particular, ultrasound imaging plays an important role in medical diagnosis due to its safety, noninvasive, inexpensiveness and real-time compared with other medical imaging techniques. However, in general the ultrasound imaging has more speckles and is low definition than the MRI (magnetic resonance imaging) and X-ray CT (computerized tomography). Therefore, it is important to improve the ultrasound imaging quality. In this study, there are three newproposals. The first is the development of a high sensitivity transducer that utilizes piezoelectric charge directly for FET (field effect transistor) channel control. The second is a proposal of a method for estimating the distribution of small scatterers in living tissue using the empirical Bayes method. The third is a super-resolution imagingmethod of scatterers with strong reflection such as organ boundaries and blood vessel walls. The specific description of each chapter is as follows: Chapter 1: The fundamental characteristics and the main applications of ultrasound are discussed, then the advantages and drawbacks of medical ultrasound are high-lighted. Based on the drawbacks, motivations and objectives of this study are stated. Chapter 2: To overcome disadvantages of medical ultrasound, we advanced our studyin two directions: designing new transducer improves the acquisition modality itself, onthe other hand new signal processing improve the acquired echo data. Therefore, the conventional techniques related to the two directions are reviewed. Chapter 3: For high performance piezoelectric, a structure that enables direct coupling of a PZT (lead zirconate titanate) element to the gate of a MOSFET (metal-oxide semiconductor field-effect transistor) to provide a device called the PZT-FET that acts as an ultrasound receiver was proposed. The experimental analysis of the PZT-FET, in terms of its reception sensitivity, dynamic range and -6 dB reception bandwidth have been investigated. The proposed PZT-FET receiver offers high sensitivity, wide dynamic range performance when compared to the typical ultrasound transducer. Chapter 4: In medical ultrasound imaging, speckle patterns caused by reflection interference from small scatterers in living tissue are often suppressed by various methodologies. However, accurate imaging of small scatterers is important in diagnosis; therefore, we investigated influence of speckle pattern on ultrasound imaging by the empirical Bayesian learning. Since small scatterers are spatially correlated and thereby constitute a microstructure, we assume that scatterers are distributed according to the AR (auto regressive) model with unknown parameters. Under this assumption, the AR parameters are estimated by maximizing the marginal likelihood function, and the scatterers distribution is estimated as a MAP (maximum a posteriori) estimator. The performance of our method is evaluated by simulations and experiments. Through the results, we confirmed that the band limited echo has sufficient information of the AR parameters and the power spectrum of the echoes from the scatterers is properly extrapolated. Chapter 5: The medical ultrasound imaging of strong reflectance scatterers based on the MUSIC algorithm is the main subject of Chapter 5. Previously, we have proposed a super-resolution ultrasound imaging based on multiple TRs (transmissions/receptions) with different carrier frequencies called SCM (super resolution FM-chirp correlation method). In order to reduce the number of required TRs for the SCM, the method has been extended to the SA (synthetic aperture) version called SA-SCM. However, since super-resolution processing is performed for each line data obtained by the RBF (reception beam forming) in the SA-SCM, image discontinuities tend to occur in the lateral direction. Therefore, a new method called SCM-weighted SA is proposed, in this version the SCM is performed on each transducer element, and then the SCM result is used as the weight for RBF. The SCM-weighted SA can generate multiple B-mode images each of which corresponds to each carrier frequency, and the appropriate low frequency images among them have no grating lobes. For a further improvement, instead of simple averaging, the SCM applied to the result of the SCM-weighted SA for all frequencies again, which is called SCM-weighted SA-SCM. We evaluated the effectiveness of all the methods by simulations and experiments. From the results, it can be confirmed that the extension of the SCM framework can help ultrasound imaging reduce grating lobes, perform super-resolution and better SNR(signal-to-noise ratio). Chapter 6: A discussion of the overall content of the thesis as well as suggestions for further development together with the remaining problems are summarized.首都大学東京, 2019-03-25, 博士（工学）首都大学東

Application of High Resolution Inversion of Ultrasonic Data to the Imaging of Multi-Layered Composite Structures

Ultrasonic imaging has evolved from its early application which utilized only amplitude C-scans to more complex techniques which make extensive use of digital signal processing. Techniques, such as one-and two-dimensional deconvolution processing and synthetic aperture focussing techniques (SAFT), are becoming more widely accepted for conventional applications. In general, each of these techniques aims to improve the interpretability of the ultrasonic image by increasing the resolution in one or more dimensions

Sparsity driven ultrasound imaging

An image formation framework for ultrasound imaging from synthetic transducer arrays based on sparsity-driven regularization functionals using single-frequency Fourier domain data is proposed. The framework involves the use of a physics-based forward model of the ultrasound observation process, the formulation of image formation as the solution of an associated optimization problem, and the solution of that problem through efficient numerical algorithms. The sparsity-driven, model-based approach estimates a complex-valued reflectivity field and preserves physical features in the scene while suppressing spurious artifacts. It also provides robust reconstructions in the case of sparse and reduced observation apertures. The effectiveness of the proposed imaging strategy is demonstrated using experimental data

Development, Optimization and Clinical Evaluation Of Algorithms For Ultrasound Data Analysis Used In Selected Medical Applications.

The assessment of soft and hard tissues is critical when selecting appropriate protocols for restorative and regenerative therapy in the field of dental surgery. The chosen treatment methodology will have significant ramifications on healing time, success rate and overall long-time oral health. Currently used diagnostic methods are limited to visual and invasive assessments; they are often user-dependent, inaccurate and result in misinterpretation. As such, the clinical need has been identified for objective tissue characterization, and the proposed novel ultrasound-based approach was designed to address the identified need. The device prototype consists of a miniaturized probe with a specifically designed ultrasonic transducer, electronics responsible for signal generation and acquisition, as well as an optimized signal processing algorithm required for data analysis. An algorithm where signals are being processed and features extracted in real-time has been implemented and studied. An in-depth algorithm performance study has been presented on synthetic signals. Further, in-vitro laboratory experiments were performed using the developed device with the algorithm implemented in software on animal-based samples. Results validated the capabilities of the new system to reproduce gingival assessment rapidly and effectively. The developed device has met clinical usability requirements for effectiveness and performance

New ultrasonic signal processing techniques for NDE applications

New ultrasonic signal processing techniques have been developed for nondestructive evaluation (NDE) applications. This dissertation has two parts. The first part is about the application of the wavelet transform to ultrasonic flaw detection. Wavelet transform is a newly developed signal analysis tool that handles time-localized signals such as an ultrasonic flaw signal quite well. A wavelet transform based signal processing technique has been developed which uses only partial knowledge of the flaw signal waveform that may be obtained from a reference experiment. The detection performance of the proposed technique is found to be comparable to that of the matched filter which requires exact knowledge of the flaw signal waveform and the noise autocorrelation function to obtain good detection performance. The proposed technique based on the wavelet transform can therefore be quite useful in situations where the flaw signal waveform is unknown or partially known. The detection performance of the proposed technique which was evaluated for hard-alpha detection in titanium samples using experimentally obtained grain noise data and simulated flaw data was very close to that of the matched filter;The second part of this dissertation describes a Kalman filter based deconvolution algorithm for ultrasonic signals and its application to material characterization and hard-alpha detection. The Kalman filter based deconvolution algorithm is based on state-space modeling of the ultrasonic measurement system. Since the Kalman filter can handle time-varying systems and non-stationary statistics quite naturally, it is better suited for such situations than the Wiener filter approach. A signal processing technique using Kalman filter based deconvolution algorithm has been developed and applied to characterize materials with different grain sizes and to detect inclusions from host material. The proposed method was tested using experimentally obtained ultrasonic data from pure titanium samples with different grain sizes. The results showed good detection performance for detecting inclusions larger that 4 mm