Compound quantitative ultrasonic tomography of long bones using wavelets analysis

Similar documents
L-band compact printed quadrifilar helix antenna with Iso-Flux radiating pattern for stratospheric balloons telemetry

Benefits of fusion of high spatial and spectral resolutions images for urban mapping

SUBJECTIVE QUALITY OF SVC-CODED VIDEOS WITH DIFFERENT ERROR-PATTERNS CONCEALED USING SPATIAL SCALABILITY

Linear MMSE detection technique for MC-CDMA

PANEL MEASUREMENTS AT LOW FREQUENCIES ( 2000 Hz) IN WATER TANK

Optical component modelling and circuit simulation

Power- Supply Network Modeling

A New Approach to Modeling the Impact of EMI on MOSFET DC Behavior

A 100MHz voltage to frequency converter

UML based risk analysis - Application to a medical robot

FeedNetBack-D Tools for underwater fleet communication

Dictionary Learning with Large Step Gradient Descent for Sparse Representations

Nonlinear Ultrasonic Damage Detection for Fatigue Crack Using Subharmonic Component

BANDWIDTH WIDENING TECHNIQUES FOR DIRECTIVE ANTENNAS BASED ON PARTIALLY REFLECTING SURFACES

3D MIMO Scheme for Broadcasting Future Digital TV in Single Frequency Networks

High finesse Fabry-Perot cavity for a pulsed laser

Enhanced spectral compression in nonlinear optical

Wireless Energy Transfer Using Zero Bias Schottky Diodes Rectenna Structures

Indoor MIMO Channel Sounding at 3.5 GHz

Modelling and Hazard Analysis for Contaminated Sediments Using STAMP Model

A sub-pixel resolution enhancement model for multiple-resolution multispectral images

High-frequency acoustic imaging with focused transducer for rapid micro echography of interfaces through buried structures non-perfectly planar

Process Window OPC Verification: Dry versus Immersion Lithography for the 65 nm node

Signal and Noise scaling factors in digital holography

Opening editorial. The Use of Social Sciences in Risk Assessment and Risk Management Organisations

A new radar sensor for cutting height measurements in tree harvesting applications

Enhancement of Directivity of an OAM Antenna by Using Fabry-Perot Cavity

The Galaxian Project : A 3D Interaction-Based Animation Engine

An improved topology for reconfigurable CPSS-based reflectarray cell,

Gis-Based Monitoring Systems.

RFID-BASED Prepaid Power Meter

On the robust guidance of users in road traffic networks

Two Dimensional Linear Phase Multiband Chebyshev FIR Filter

A multi-sine sweep method for the characterization of weak non-linearities ; plant noise and variability estimation.

QPSK-OFDM Carrier Aggregation using a single transmission chain

Concepts for teaching optoelectronic circuits and systems

Dynamic Platform for Virtual Reality Applications

Simulation Analysis of Wireless Channel Effect on IEEE n Physical Layer

Resonance Cones in Magnetized Plasma

Indoor Channel Measurements and Communications System Design at 60 GHz

analysis of noise origin in ultra stable resonators: Preliminary Results on Measurement bench

A generalized white-patch model for fast color cast detection in natural images

Influence of ground reflections and loudspeaker directivity on measurements of in-situ sound absorption

Small Array Design Using Parasitic Superdirective Antennas

A high PSRR Class-D audio amplifier IC based on a self-adjusting voltage reference

Reference Selection for an Active Ultrasound Wild Salmon Monitoring System

Application of the multiresolution wavelet representation to non-cooperative target recognition

Long reach Quantum Dash based Transceivers using Dispersion induced by Passive Optical Filters

Study on a welfare robotic-type exoskeleton system for aged people s transportation.

Sound level meter directional response measurement in a simulated free-field

On the role of the N-N+ junction doping profile of a PIN diode on its turn-off transient behavior

Augmented reality as an aid for the use of machine tools

Analysis of the Frequency Locking Region of Coupled Oscillators Applied to 1-D Antenna Arrays

VR4D: An Immersive and Collaborative Experience to Improve the Interior Design Process

A design methodology for electrically small superdirective antenna arrays

Performance of Frequency Estimators for real time display of high PRF pulsed fibered Lidar wind map

Gate and Substrate Currents in Deep Submicron MOSFETs

Adaptive noise level estimation

Analytic Phase Retrieval of Dynamic Optical Feedback Signals for Laser Vibrometry

Reconfigurable antennas radiations using plasma Faraday cage

New paradigm in design-manufacturing 3Ds chain for training

A Novel Piezoelectric Microtransformer for Autonmous Sensors Applications

Development of an On-Chip Sensor for Substrate Coupling Study in Smart Power Mixed ICs

A STUDY ON THE RELATION BETWEEN LEAKAGE CURRENT AND SPECIFIC CREEPAGE DISTANCE

PMF the front end electronic for the ALFA detector

A modal method adapted to the active control of a xylophone bar

Design of Cascode-Based Transconductance Amplifiers with Low-Gain PVT Variability and Gain Enhancement Using a Body-Biasing Technique

Analogic fiber optic position sensor with nanometric resolution

Towards Decentralized Computer Programming Shops and its place in Entrepreneurship Development

Modelling and Analysis of Static Transmission Error. Effect of Wheel Body Deformation and Interactions between Adjacent Loaded Teeth

Radio direction finding applied to DVB-T network for vehicular mobile reception

Design of an Efficient Rectifier Circuit for RF Energy Harvesting System

INVESTIGATION ON EMI EFFECTS IN BANDGAP VOLTAGE REFERENCES

Floating Body and Hot Carrier Effects in Ultra-Thin Film SOI MOSFETs

Neel Effect Toroidal Current Sensor

Performance comparison of pulse-pair and wavelets methods for the pulse Doppler weather radar spectrum

DUAL-BAND PRINTED DIPOLE ANTENNA ARRAY FOR AN EMERGENCY RESCUE SYSTEM BASED ON CELLULAR-PHONE LOCALIZATION

MODELING OF BUNDLE WITH RADIATED LOSSES FOR BCI TESTING

A technology shift for a fireworks controller

Geometrical Specification Model for Gear - Expression, Metrology and Analysis

Characterization of Few Mode Fibers by OLCI Technique

A Low-Profile Cavity-Backed Dual-Polarized Spiral Antenna Array

NOVEL BICONICAL ANTENNA CONFIGURATION WITH DIRECTIVE RADIATION

An Operational SSL HF System (MILCOM 2007)

Convergence Real-Virtual thanks to Optics Computer Sciences

UV Light Shower Simulator for Fluorescence and Cerenkov Radiation Studies

A notched dielectric resonator antenna unit-cell for 60GHz passive repeater with endfire radiation

Scaled Laboratory Experiments of Shallow Water Acoustic Propagation

Direct optical measurement of the RF electrical field for MRI

High acquisition rate infrared spectrometers for plume measurement

Computational models of an inductive power transfer system for electric vehicle battery charge

Demand Response by Decentralized Device Control Based on Voltage Level

100 Years of Shannon: Chess, Computing and Botvinik

Measures and influence of a BAW filter on Digital Radio-Communications Signals

Stewardship of Cultural Heritage Data. In the shoes of a researcher.

Globalizing Modeling Languages

Probabilistic VOR error due to several scatterers - Application to wind farms

An image segmentation for the measurement of microstructures in ductile cast iron

New Structure for a Six-Port Reflectometer in Monolithic Microwave Integrated-Circuit Technology

A Low-cost Through Via Interconnection for ISM WLP

Transcription:

Compound quantitative ultrasonic tomography of long bones using wavelets analysis Philippe Lasaygues To cite this version: Philippe Lasaygues. Compound quantitative ultrasonic tomography of long bones using wavelets analysis. M. André. Compound quantitative ultrasonic tomography of long bones using wavelets analysis, 2005, United States. Springer, 28, pp.223-229, 2007, <10.1007/1-4020-5721-0_24>. <hal- 00440737> HAL Id: hal-00440737 https://hal.archives-ouvertes.fr/hal-00440737 Submitted on 11 Dec 2009 HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers. L archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

COMPOUND QUANTITATIVE ULTRASONIC TOMOGRAPHY OF LONG BONES USING WAVELETS ANALYSIS Philippe Lasaygues Laboratoire de Mécanique et d'acoustique, UPR CNRS 7051, 31 chemin Joseph AIGUIER 13402 Marseille cedex 20 Abstract: Key words: Compound Quantitative Ultrasonic Tomography (CQUT) is used to long bones imaging. We showed that an iterative tool might be used to provide, from reflection tomography, qualitative images of the shape of the object, and to provide, from transmission tomography, quantitative images of the velocity map. Both tomographies are based on ultrasonic propagation in bones, particularly perturbed by this high-contrasted heterogeneous medium. Reflected and transmitted signal are composed of several packages of waves. We propose a novel algorithm based on the wavelet analysis tool adapted to ultrasonic signals that allow the transmitted signals to be cleaned and filtered and the useful information to be separated from the unwanted noise. Ultrasonic tomography, bone imaging, wavelets analysis 1. INTRODUCTION Quantitative ultrasonic tomography is used to long bones imaging. In previous works, we showed that an iterative tool might be used to provide, from reflection tomography, qualitative images of the shape of the object, and to provide, from transmission tomography, quantitative images of the velocity map. Both tomographies are based on ultrasonic propagation in bones, particularly perturbed by this high-contrasted heterogeneous medium. Reflected and transmitted signal are composed of several packages of waves, which had followed various pathways within the cortical shell of long bones. Signal and image processing have a large important part in the complete process. Transmission tomography is based on the knowledge of the Time-Of- Flight (TOF) measured on the extracted ultrasonic signal and related to the

2 Philippe Lasaygues useful velocities of the ultrasonic waves throughout the shell. To optimise our TOF detection processing, we propose a novel signal processing based on the wavelet analysis adapted to ultrasonic signals that allow the transmitted signals to be cleaned and filtered and the useful information to be separated from the unwanted noise. 2. COMPOUND ULTRASONIC TOMOGRAPHY Ultrasonic Tomography (UT) is based on a linearization of the inverse acoustic scattering problem. It allows perturbations (theoretically small) of a reference medium to be visualized. For long bones, the acoustic impedance of the cortical is highly contrasted compared to the surrounding soft tissues (or water), and so, perturbs ultrasonic propagation (refraction, attenuation and diffraction). Therefore, bone imaging is a non-linear inverse problem with no single solutions. However, for some hypotheses (weak heterogeneous medium, straight lines waves propagation) and a priori knowledge of the analyzed bones, we showed that one solution is possible [1]. Our algorithm combines two tomographies, Ultrasonic Reflection Tomography (URT) and Ultrasonic Transmission Tomography (UTT). URT [2] provides the realistic shape of the bodies and UTT is especially adapted to mapping local velocities inside the area defined. On the final image, the mean geometrical and acoustical parameters may be measured and used to improve the dimensioning and the quantification, ensuring, the complete iterative Compound UT in which the behaviour and the convergence are excellent. Long human bones are effectively irregular hollow tubes (Figure 1) and should support the propagation of more complex waves similar to elastic volume waves. Close to the external water/bone interface, the wave is deviated from its initial position because of the refraction throughout the shell. The corrections of this refraction are based on the knowledge of the dimensions and on a velocity ground level (Snell-Descartes laws).

COMPOUND QUANTITATIVE ULTRASONIC TOMOGRAPHY OF LONG BONES USING WAVELETS ANALYSIS 3 Figure 1. Supposed pathways of the wave propagation throughout cortical long bone for CQUT 3. WAVELET ANALYSIS OF TRANSMITTED SIGNALS For transmission tomography, projections must be constructed from determination of the TOF between both transducers. The experimental ultrasonic acquisition protocol enabled us, according to the compensation procedure, to digitize correctly the transmitted signals [3], composed of several packages of waves which had followed various pathways within the object (cortex and/or medullary and/or cortex) [4]. Various algorithms of existing T.O.F. measurements are based on various definitions of transit times of pulse, such as zero-crossing or thresholding. Our algorithm was based on detecting the arrival time of the signal corresponding to the wave propagating at, or close to, the longitudinal velocity within the shell ("fluid" modeling). The level of the detection accepted for each signal was fixed by a subjective threshold. Then, the number of bad T.O.Fs detected increased because these thresholds are very sensitive to the signal-to-noise- ratio (SNR), and are not correctly adapted to all the values (usually 15 6 or 20 6 ) that we had to automatically treat. Indeed, the choice of the lower limit of the threshold was fixed for each view angle in relation to the noise. If we extracted the useful signal from this noise and if we standardized all the packages with a criterion which did not perturb their temporal position, the errors decrease.

4 Philippe Lasaygues Figure 2. Transmitted signal in water Our idea was to exploit wavelet analysis, which allowed simultaneous joint studies in the time and frequency domains. Two aspects were considered in particular: the significant spectral bandwidth of the emitted signal; that is, the transmitted signal without the object (Figure 2) and the axial resolution of the transducer, brought back to the support of this emitted signal. Then we defined a Region-Of-Interest (R.O.I.) consisting of this time and frequency bandwidth for the wavelet analysis, filtering and synthesis (sum of the wavelet coefficients). Figure 3. Wavelets analysis of a transmitted signal through the medullary hollow zone (Zone 1)

COMPOUND QUANTITATIVE ULTRASONIC TOMOGRAPHY OF LONG BONES USING WAVELETS ANALYSIS 5 Figure 4. Filtering the transmitted signal corresponding to the wave propagating close to the boundary bone/water (zone 3) Figure 5. Filtering transmitted signals in solid cortex area (zone 2) (a - temporal representation; b - wavelet decomposition modulus and visualization of the R.O.I. on the bandwidth [0.8-1.5] MHz (at 25 db) X [115 123.35] µs; c - comparison between the initial signal and its synthesis on the R.O.I.; d zoom, the SNR and the resolution of the firstarriving signal increase.)

6 Philippe Lasaygues In the medullary hollow zone, the processing of the signals is generally easy and the T.O.F. are calculated with a margin for error lower than 2% (Figure 3). In the solid cortex (Figure 4 & 5), signals are much more disturbed. For example, on the Figure 5, three lots of waves were distinguished with different amplitudes and the first-arriving signal was the lowest and the arbitrary first-time could be 115 µs. The upper signal arrived in water after the T.O.F. (123.35 µs) and so, did not correspond to a direct propagation inside the shell. This signal was analyzed and the R.O.I. was defined. Figures 4-d display the curves after filtering. The SNR was higher than before processing, and the signal was filtered for low and high frequencies. The resolution and the quantification of the first two packages had improved. Hence, the T.O.F. detections increased and the false values decreased. The same threshold can be used for all filtering signals for all view angles. The advantage of our algorithm was that the noise reduction, the thresholding and the T.O.F. detection could be made simultaneously on the image of the wavelet coefficients. 4. CQUT OF FEMUR The human sample was a femoral specimen of approximately 32 ± 5 mm for the external diameter and 16 ± 2 mm for the internal diameter. The a priori mean-velocities were 3400 m/s in bone and 1478 m/s in water. In previous studies, the femur was compared to a circular tube, but this approximation (Figure 6-left) was not suitable enough. Figure 6-right displays the final image obtained by CQUT. Figure 6. Human femur CQUT, without (left) and with(right) signal processing

COMPOUND QUANTITATIVE ULTRASONIC TOMOGRAPHY OF LONG BONES USING WAVELETS ANALYSIS 7 The fluid in the internal shape was reconstructed with a correct velocity value ( 1500 m/s) and the dimension of this cavity was 15 17 mm. The greatest external shape (30 34 mm) was exact and the mean-velocity estimated in the cortical shell was 3150 m/s. 5. REFERENCES [1] Ouedraogo E., Lasaygues P., Lefebvre JP., Gindre M., Talmant M. and Laugier P., Contrast and velocity ultrasonic tomography of long bones, Ultrasonic Imaging 24, (2002) 135 146. [2] Lasaygues P., Lefebvre JP., Cancellous and cortical bone imaging by reflected tomography, Ultrasonic Imaging 23, (2001) 55-68 [3] Ouedraogo E, Lasaygues P, Lefebvre J.P., Talmant M, Gindre M., and Laugier P., Multistep compensation technique for ultrasound tomography of bone, Acoustical Imaging, 26, (2001), 153-160 [4] Ouedraogo E. 2002 Caractérisation géométrique et Acoustique des os long par tomographie ultrasonore, Ph.D Thesis University of Cergy-Pontoise