WO2018109314A1 - Procede de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage a ultrasons correspondants - Google Patents

Procede de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage a ultrasons correspondants Download PDF

Info

Publication number
WO2018109314A1
WO2018109314A1 PCT/FR2017/053391 FR2017053391W WO2018109314A1 WO 2018109314 A1 WO2018109314 A1 WO 2018109314A1 FR 2017053391 W FR2017053391 W FR 2017053391W WO 2018109314 A1 WO2018109314 A1 WO 2018109314A1
Authority
WO
WIPO (PCT)
Prior art keywords
matrix
singular
signals
emission
ultrasonic
Prior art date
Application number
PCT/FR2017/053391
Other languages
English (en)
Inventor
Sébastien ROBERT
Eduardo-Rigoberto LOPEZ VILLAVERDE
Claire PRADA JULIA
Original Assignee
Commissariat à l'énergie atomique et aux énergies alternatives
Centre National De La Recherche Scientifique
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Commissariat à l'énergie atomique et aux énergies alternatives, Centre National De La Recherche Scientifique filed Critical Commissariat à l'énergie atomique et aux énergies alternatives
Priority to EP17821680.0A priority Critical patent/EP3555611A1/fr
Priority to JP2019532089A priority patent/JP2020502518A/ja
Priority to CA3046106A priority patent/CA3046106A1/fr
Priority to US16/469,918 priority patent/US11054401B2/en
Publication of WO2018109314A1 publication Critical patent/WO2018109314A1/fr

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/4463—Signal correction, e.g. distance amplitude correction [DAC], distance gain size [DGS], noise filtering
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13—Tomography
    • A61B8/14—Echo-tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/54—Control of the diagnostic device
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/043—Analysing solids in the interior, e.g. by shear waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
    • G01N29/0609—Display arrangements, e.g. colour displays
    • G01N29/0645—Display representation or displayed parameters, e.g. A-, B- or C-Scan
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
    • G01N29/0654—Imaging
    • G01N29/069—Defect imaging, localisation and sizing using, e.g. time of flight diffraction [TOFD], synthetic aperture focusing technique [SAFT], Amplituden-Laufzeit-Ortskurven [ALOK] technique
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/22—Details, e.g. general constructional or apparatus details
    • G01N29/26—Arrangements for orientation or scanning by relative movement of the head and the sensor
    • G01N29/262—Arrangements for orientation or scanning by relative movement of the head and the sensor by electronic orientation or focusing, e.g. with phased arrays
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/46—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by spectral analysis, e.g. Fourier analysis or wavelet analysis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/48—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by amplitude comparison
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8977—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8995—Combining images from different aspect angles, e.g. spatial compounding
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52077—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging with means for elimination of unwanted signals, e.g. noise or interference
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5246—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
    • A61B8/5253—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode combining overlapping images, e.g. spatial compounding
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/10—Number of transducers
    • G01N2291/106—Number of transducers one or more transducer arrays
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8997—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using synthetic aperture techniques
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10132—Ultrasound image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20048—Transform domain processing
    • G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing

Definitions

  • the present invention relates to a method for processing signals derived from ultrasonic sounding acquisition to perform imaging or to a method for processing signals resulting from an acquisition by ultrasonic probing. adaptive and selective focus. It also relates to a computer program and a corresponding ultrasonic sounding device.
  • the invention applies in particular to the field of non-destructive ultrasonic testing, in which the acquisition of ultrasonic signals makes it possible to detect and / or visualize defects in structures, but it can also be applied to any type of detection or ultrasonic ultrasound imaging, particularly in the medical field for the inspection of areas of interest in the human or animal body.
  • N receiving transducers so as to simultaneously receive, for a predetermined duration, for each transmission, N measurement time signals, in particular measuring echoes due to reflections of the emission in question,
  • each MR coefficient (t) of this matrix representing the measurement signal received by the ith receiving transducer due to the ith emission .
  • phased array probing device in which each transducer is both transmitter and receiver, switching between these two modes can be controlled electronically.
  • the sensor may be brought into contact with the object to be probed or remotely, but in the latter case it must be immersed to ensure the transmission of the ultrasonic waves in the object to be probed.
  • This sensor can be linear (1 D) or matrix (2D), with rigid or flexible elements.
  • the matrix [MR (t)] of temporal signals obtained by this type of acquisition can then be the subject of a processing, in particular for providing an image of the area of interest inspected or for the extraction of significant parameters of structural defects in the area of interest inspected.
  • this processing can be embedded in control instruments for real-time processing.
  • plane-wave compounding or “plane-wave imaging”
  • plane-wave imaging This type of acquisition, generally referred to as “plane-wave compounding” or “plane-wave imaging”
  • Montaldo et al entitled “Coherent plane-wave compounding for very high frame rate ultrasound transient elastography, "published in IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control, vol. 56, No. 3, pages 489-506, March 2009. It is mainly considered in the medical field and implemented in some ultrafast ultrasound machines for imaging maps of the elasticity of the human body.
  • the resulting algorithms and the associated applications lend themselves particularly well to a parallelization of the calculations on processors of the GPU type (of the English "Graphie Processing Unit”) implanted in the graphics cards of the computers.
  • the performances of an ultrasound system implementing a "plane-wave compounding" technique can thus reach in practice 10,000 images / s.
  • Another advantage of this acquisition technique lies in the fact that each shot is made by soliciting all the emission transducers so that the energy emitted is high, making this method less sensitive to the phenomena of attenuation, of electronic noise. or structure.
  • This type of acquisition is also used in the patent application WO 2015/092250 A1 by cleverly adapting the principle of synthetic focusing in all points, so as to take advantage of the simplicity of the technique of "plane-wave compounding" in to achieve a high acquisition rate and image quality, in terms of spatial resolution and contrast, related to a focus synthesized in all points of the desired image.
  • This made it possible to consider the "plane-wave compounding" technique for non-destructive testing applications.
  • the images obtained in non-destructive testing may have a significant noise level depending on the properties of the medium being probed.
  • it is an electronic noise when the material is homogeneous and viscoelastic, or a structure noise when the waves are diffused by heterogeneities of the material.
  • Images calculated with focusing techniques in all points are also affected by artifacts, or ghost echoes, related to geometry echoes, for example the echo of a part background near that of a defect.
  • N receiving transducers so as to simultaneously receive, for a predetermined duration, for each transmission, N measurement time signals, in particular measuring echoes due to reflections of the emission in question, and
  • the plane wave matrix [MR (t)] is not of the same nature as the matrixes of inter-element impulse responses usually obtained by conventional acquisition techniques exploiting the synthetic focus in all points, it has been unexpectedly found that a noise filtering method based on a singular value decomposition of a transform of the plane of the plane waves in the frequency domain provides surprising results in terms of noise attenuation.
  • the denoised matrix [MRu (t)] thus obtained makes it possible, in particular, for a picture reconstruction of quality which is markedly improved with respect to what the plane wave matrix [MR (t)] produces without this treatment.
  • the reconstruction of the denoised matrix [MRu (t)] of time signals comprises a reconstruction of a denoised matrix [FTMRu (f)] of frequency signals from the singular values and singular vectors not eliminated, then an inverse transformation of this de-energized matrix [FTMRu (f)] of frequency signals in the denoised matrix [MRu (t)] of time signals.
  • the transformation and the inverse transformation are discrete Fourier transformations.
  • the theoretical decay curve is defined by a reciprocal function F "1 (1 -a) itself defined by a function F (a), called the distribution of random singular values, such that:
  • an ultrasonic signal processing method may further comprise a reconstitution of an imaged area by calculating, at each point of a plurality of predetermined points of this imaged area, a resulting value. a consistent summation of values snapshots taken respectively by at least a part of the NxM time signals of the matrix [MRu (t)] at flight times respectively corresponding to a passage through the point considered according to a predetermined propagation mode.
  • the calculation is done on a part of the NxM time signals of the matrix [MRu (t)] in a restricted angular sector in the set of successive transmissions, this restricted angular sector being selected so that the waves planes that are excluded do not interact with at least one defect related to singular values and singular vectors not eliminated.
  • the restricted angular sector is selected on the basis of a comparison, for at least one of the non-eliminated singular values, of an experimental phase value of the singular vector associated with it in transmission with a value. theoretical phase in the presence of said at least one defect related to this singular vector.
  • a computer program downloadable from a communication network and / or recorded on a computer readable medium and / or executable by a processor, comprising instructions for performing the steps of a signal processing method ultrasound apparatus according to the invention, when said program is executed on a computer.
  • an ultrasonic sounding device comprising:
  • a probe comprising L ultrasonic emission transducers and N ultrasonic reception transducers
  • N receiving transducers means for controlling the N receiving transducers so as to simultaneously receive, for a predetermined duration, for each transmission, N measurement time signals, in particular measuring echoes due to reflections of the emission in question, and
  • a processor for reconstituting a matrix [MR (t)] of ultrasonic time signals of size NxM, each MR coefficient (t) of this matrix representing the measurement signal received by the i-th reception transducer due to the -th issue,
  • the processor being further configured to perform the following processes: transformation of the matrix [MR (t)] of time signals into a matrix [FTMR (f)] of frequency signals, then decomposition into singular values of the matrix [FTMR (f)] of frequency signals,
  • FIG. 1 schematically represents the general structure of an ultrasonic sounding device according to one embodiment of the invention
  • FIG. 2 illustrates a principle of successive transmissions of plane ultrasonic waves implemented by the device of FIG. 1,
  • FIG. 3 graphically illustrates the result of a singular value decomposition of a matrix of frequency signals obtained by transforming a plane wave matrix itself obtained by means of transmissions carried out according to the principle of the figure 2,
  • FIGS. 4 and 5 diagrammatically and graphically illustrate the principle of a restricted angular sector selection for improving the detection of a fault, according to one embodiment of the invention
  • FIG. 6 illustrates three examples of images reconstituted without or with application of the principles of the invention.
  • FIG. 7 illustrates the successive steps of a method for acquiring and processing ultrasonic signals implemented by the device of FIG. 1, according to one embodiment of the invention.
  • a sounding device 100 of an object 102 comprises an ultrasound probe 104 having a housing 106, that is to say an element of indeformable structure which serves as a reference attached to the probe 104, in which are arranged, for example linearly or by matrix, N transducers 108 ; ..., 108 N fixed or mobile arranged in a network.
  • the object 102 is, for example, a mechanical part that one wishes to examine by non-destructive testing or, in a medical context, a part of a human or animal body that one wishes to control in a non-invasive manner.
  • the object 102 is immersed in a liquid, such as water 1 10, and the probe 104 is kept at a distance from the object 102 so that the water 1 10 separates them.
  • the probe 104 could be in direct contact with the object 102.
  • Transducers 108 ; ..., 108 N are adapted to individually transmit ultrasonic waves towards the object 102 in response to control signals identified as C, along principal directions parallel to each other, indicated by arrows in dotted in Figure 1, and in a main plane which is that of the figure.
  • Transducers 108 ; ..., 108 N are further designed to detect echoes of the ultrasonic waves reflected on and in the object 102 and to provide measurement signals identified under the general reference S and corresponding to these echoes.
  • the transducers 108 ; ..., 108 N fulfill both the transmission and reception functions, but different transmitters receivers could also be provided in different and independent boxes while remaining consistent with the principles of the invention.
  • the number L of transmitters could quite be different from the number N of receivers.
  • the sounding device 100 further comprises an electronic circuit 1 12 for controlling the transducers 108 1,..., 108 N of the probe 104 and for processing the measurement signals S.
  • This electronic circuit 1 12 is connected to the probe 104 in order to to transmit to it the control signals C and to receive the measurement signals S.
  • the electronic circuit 1 12 is for example that of a computer. It has a central processing unit 1 14, such as a microprocessor designed to transmit the control signals C to the probe 104 and to receive the measurement signals S from the probe 104, and a memory 116 in which is recorded in particular. a computer program 1 18.
  • the computer program 1 18 first comprises instructions 120 for generating the control signals C of the transducers 108 ; ..., 108 N and receive their echoes. These instructions are more precisely programmed to:
  • the plane ultrasonic waves are obtained on emission by applying to the transducers 108 ; ..., 108 N of the delay laws stored in memory 1 16 in a base 122 of delay laws.
  • Each delay law defines delays to be applied to the transducers 108 ; ..., 108 N in transmission, so as to generate a plane ultrasonic wave at a desired emission angle among the M different successive emission angles.
  • the first plane wave emission is associated with a delay law ⁇ relating to signals emitted by the transducers I O8 1 , ..., 108 N , allowing the emission of a plane wave of emission angle ⁇ 1 with respect to the direction z in a first emission zone ZEi partially located outside the opening of the probe 104.
  • the (M + 1) / 2-th emission of planar wave is associated with a uniform law of delays T (M + i ) / 2 for the emission of a plane wave of emission angle zero with respect to the direction z in a (M + 1) / 2-th emission zone ZE (M + 1) / 2 covering the opening of the probe 104.
  • the area to be imaged must be contained in the union of the M successive emission areas. As a result, this zone can extend beyond the opening of the probe 104, as can be seen in FIG. 2.
  • the imaged zone can take the form of a zone sectoral bounded by the ends of the emission zones of maximum and minimum angles. It is thus possible to obtain an image of S-scan type.
  • the M successive emission angles different to ⁇ ⁇ can be defined around an average direction ⁇ ( ⁇ + ⁇ ) / 2 not perpendicular to the array of transducers 108 ; ..., 108 N.
  • this crack being moreover perpendicular to the transducer array, it is preferable to shift the area to be inspected with respect to the probe 104 and emit around an average of 45 ° for example. The area to be inspected can even be shifted to the point of completely leaving the opening of the probe 104.
  • an apodization of the ultrasonic signals emitted by the transducers 108i,..., 108 N to form a better plane ultrasound wave. quality without distortion due to edge effects.
  • Such apodization is performed on the occasion of each emission spatially on all the transducers using an apodization window such as a trapezoidal amplitude law, Hamming or Blackman-Harris. It has the result of providing a better definition of the successive emission areas.
  • the set S of the NxM measurement time signals received by the N transducers 108 ; ..., 108 N is returned by the probe 104 to the central processing unit 1 14.
  • the computer program 1 18 then further comprises instructions 124 for constructing a matrix [MR (t)] of ultrasonic time signals of size NxM, referred to as the matrix of the plane waves.
  • Each coefficient MR, j (t) of this matrix represents the measurement signal received by the transducer 108, in response to the jth emission.
  • the computer program 1 18 further comprises instructions 126 for temporal filtering of the matrix [MR (t)], this filtering to remove any information at flight times excluded from the zone. of interest in object 102.
  • the computer program 1 18 further comprises instructions 128 for transforming the matrix [MR (t)] into a matrix [FTMR (f)] of frequency signals by Fourier transform, advantageously by discrete Fourier transform after time sampling of the ultrasonic signals forming the coefficients of the matrix [MR (t)], or even more advantageously by calculation of FFT (of the English "Fast Fourier Transform” ) if the number of samples of each coefficient of the matrix [MR (t)] allows it.
  • the computer program 1 18 further comprises instructions 130 for breaking down the matrix [FTMR (f)] of frequency signals into singular values over a frequency band.
  • FTMR (f) matrix of inter-element impulse responses usually obtained by conventional acquisition techniques exploiting synthetic focus in all points
  • this operation is not equivalent when it is applied to a matrix such as the matrix [FTMR (f)].
  • ⁇ [FTMR (f)] where " ⁇ " is the symbol of the conjugate transpose of a matrix , does not represent the time reversal operator in emission as for the matrix of the interelements impulse responses.
  • the computer program 1 18 further includes instructions 132 for reducing the rank of the matrix [FTMR (f)], eliminating a portion of the singular values a, (f).
  • the function F (a) gives values between 0 and 1 on the support interval ⁇ e [0; 2], where ⁇ is the singular value variable. It is also strictly growing. Its inverse function F "1 (a) is thus also strictly increasing, so that the function F ⁇ 1 (1 -a) gives the desired number of singular values decay curve, to a constant of proportionality to adjust it to the experimental curve.
  • M 64 plane waves between -31.5 ° and +31, 5 °.
  • P (f) 2 to 5 MHz.
  • the matrix [FTMRu (f)] thus reconstituted is a denoised matrix of frequency signals, the noise subspace represented by the matrix [FTMR N (f)] having been eliminated.
  • the computer program 1 1 8 further comprises instructions 134 for transforming the matrix [FTMRu (f)] into a denoised matrix [MRu (t)] of time signals by inverse Fourier transform, advantageously by inverse discrete Fourier transform or, even more advantageously, by calculating IFFT (Inverse Fast Fourier Transform) if the number of samples of each coefficient of the matrix [FTMRu (f)] allows it.
  • IFFT Inverse Fast Fourier Transform
  • the computer program 1 18 includes instructions, designated by the general reference 136, for processing the matrix [MRu (t)].
  • the processing carried out by the instructions 1 36 may include a digital image reconstruction of the area of interest in the object 102 by adapting the synthetic focusing principle in all respects, as taught for example in the document WO 2015/092250 A1. . This restores a digital image of the area of interest the quality is better than if the reconstitution had been carried out on the non-denoised matrix [MR (t)]. In particular, the Signal to Noise Report (SNR) is improved.
  • the processing performed by instructions 136 could include adaptive and selective focusing.
  • the module A (P) of a coherent summation involving the NxM temporal signals of the matrix is calculated [ MRu (t)] to NxM flight times calculated according to a predetermined propagation mode, each flight time t being the time taken by the j-th plane wave to be received by the ith receiving transducer through the pixel considered according to the predetermined propagation mode:
  • ⁇ ( ⁇ ) and pj (P) are weighting coefficients respectively in transmission and reception whose expressions depend on the application considered to take into account phenomena or processes such as a filtering of geometry echoes, a attenuation compensation due to spatial spread of waves, etc.
  • This comparison is made at a chosen frequency f c which can be the central operating frequency of the probe 104, a frequency for which the singular value ai (f) takes its highest value, or any other predetermined frequency.
  • f c can be the central operating frequency of the probe 104, a frequency for which the singular value ai (f) takes its highest value, or any other predetermined frequency.
  • the comparison is thus made at 5 MHz for the singular transmission vector v ⁇ f) corresponding to a fault location (s) D with coordinates (X D , Z D ).
  • the coordinates (X D , Z D ) can be determined with the useful signal matrix [FTMRu (f)] by calculating the retro-propagation of the singular vector in reception U i (f) at the frequency f c , as for example taught in the article by Lopez Villaverde et al, entitled “Ultrasonic imaging of defects in coarse-grained steels with the decomposition of the time reversal operator", published in Journal of the Acoustical Society of America, volume 140, No. 1, pages 541 - 550 (2016).
  • the experimental phase value of the singular vector in emission v ⁇ y is also a corrected phase calculated in the following way to be always negative:
  • Vy, 1 ⁇ y M M, ⁇ 1 (c) arg [i 1 (c) ]. - max [arg [ 3 ⁇ 4 e) ].].
  • the values finally retained for ⁇ and m 2 can be respectively the minimum and the maximum of the values found for each of the values singular of the useful signal matrix.
  • the advantageous calculation mode detailed above can also be combined with an adaptation of the synthetic focusing principle in all respects as taught in the document WO 2015/092250 A1.
  • FIG. 6 illustrates, in an example of a probed object comprising a central circular defect D, three images obtained:
  • an exemplary method 700 for acquiring and processing ultrasonic signals that can be implemented by the device 100 of FIG. 1 will now be described according to a preferred embodiment of the invention.
  • the processing unit 1 14 executing the instructions 120 controls the transmission and reception sequences of the transducers 108i,..., 108 N for the acquisition of the measurement signals MRj j (t ) of the matrix [MR (t)].
  • Steps 702 and 704 can be executed simultaneously, i.e., it is not necessary to wait until all shots are fired to begin recording the measurement signals and perform processing such as image reconstruction.
  • the processing unit 114 executing the instructions 126 carries out a temporal filtering of the matrix [MR (t)], this filtering being intended to delete any information at flight times excluded from the area of interest.
  • This step 706 is intended to then facilitate the separation of the two subspaces represented by the [FTMRu (f)] and [FTMR N (f)] matrices, especially when the defects to be imaged are close to a strongly echogenic interface, such as a room floor. It makes it possible to limit the zone to be imaged to a neighborhood close to the defects by excluding in particular the disturbing echogenic interfaces. She finds all her interest in the imagery of cracks forming from the bottom of the object.
  • the processing unit 1 14 executing the instructions 130 performs a singular value decomposition of the matrix [FTMR (f)], as detailed previously.
  • the processing unit 1 14 executing the instructions 132 reduces the rank of the matrix [FTMR (f)] keeping only the useful signal matrix [FTMRu (f)].
  • the processing unit 1 14 executing the instructions 134 performs a discrete inverse Fourier transform of the matrix [FTMRu (f)] to obtain the denoised matrix [MRu (t)] of time signals.
  • the processing unit 1 14 executing the instructions 136 selects, in an optional but advantageous manner, a restricted angular sector in the set of successive transmissions outside which the plane waves do not interact. with the defect (s) to be detected.
  • This narrow angular sector is defined by its minimum emission index (ITH) and maximum (m 2 ) for example according to the method detailed above.
  • the processing unit 1 14 always executing the instructions 136 reconstitutes and displays a digital image of the effective zone of interest by adapting the synthetic focusing principle in all points from the denoised matrix [MRu (t)] in the selected restricted angular sector.
  • the computer program instructions could be replaced by electronic circuits dedicated to the functions performed during the execution of these instructions.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Acoustics & Sound (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Analytical Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Signal Processing (AREA)
  • Radiology & Medical Imaging (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Mathematical Physics (AREA)
  • Quality & Reliability (AREA)
  • Theoretical Computer Science (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

Ce procédé (700) de traitement de signaux ultrasonores comporte la commande (702) de transducteurs d'émission pour M émissions successives d'ondes ultrasonores planes présentant M angles d'émission différents, la commande (702) de N transducteurs de réception pour recevoir simultanément N signaux temporels de mesure par émission, et l'obtention (704) d'une matrice [MR(t)] de signaux temporels ultrasonores, chaque coefficient MRi,j(t) de cette matrice représentant le signal de mesure reçu par le i-ième transducteur de réception dû à la j-ième émission. Il comporte en outre une décomposition en valeurs singulières (710) d'une matrice [FTMR(f)] de signaux fréquentiels obtenue (708) par transformation de la matrice [MR(t)], une élimination (712) d'une partie des valeurs singulières et une reconstitution (714) d'une matrice débruitée [MRU(t)] de signaux temporels à partir des valeurs singulières non éliminées.

Description

PROCEDE DE TRAITEMENT DE SIGNAUX ISSUS D'UNE ACQUISITION PAR SONDAGE ULTRASONORE, PROGRAMME D'ORDINATEUR ET DISPOSITIF DE SONDAGE A ULTRASONS CORRESPONDANTS La présente invention concerne un procédé de traitement de signaux issus d'une acquisition par sondage ultrasonore pour réaliser de l'imagerie ou de la focalisation adaptative et sélective. Elle concerne également un programme d'ordinateur et un dispositif de sondage à ultrasons correspondants.
L'invention s'applique en particulier au domaine du contrôle non destructif par ultrasons, dans lequel l'acquisition de signaux ultrasonores permet de détecter et/ou visualiser des défauts dans des structures, mais elle peut aussi s'appliquer à tout type de détection ou imagerie échographique ultrasonore, notamment au domaine médical pour l'inspection de zones d'intérêt dans le corps humain ou animal.
Elle porte plus particulièrement sur un procédé de traitement acquérant les signaux ultrasonores de la façon suivante :
- commande de L transducteurs d'émission pour M émissions successives d'ondes ultrasonores planes présentant M angles d'émission successifs différents dans M zones d'émission,
- commande de N transducteurs de réception de manière à recevoir simultanément et pendant une durée prédéterminée, pour chaque émission, N signaux temporels de mesure, mesurant en particulier des échos dus à des réflexions de l'émission considérée,
- obtention d'une matrice [MR(t)] de signaux temporels ultrasonores de taille NxM, chaque coefficient MR (t) de cette matrice représentant le signal de mesure reçu par le i-ième transducteur de réception dû à la j-ième émission.
Une telle acquisition se fait généralement à l'aide d'un dispositif de sondage à capteur multiéléments, dans lequel chaque transducteur est à la fois émetteur et récepteur, une commutation entre ces deux modes pouvant être commandée électroniquement. Le capteur peut être mis au contact de l'objet à sonder ou à distance, mais dans ce dernier cas il doit être immergé pour assurer la transmission des ondes ultrasonores dans l'objet à sonder. Ce capteur peut être linéaire (1 D) ou matriciel (2D), à éléments rigides ou flexibles.
La matrice [MR(t)] de signaux temporels obtenue par ce type d'acquisition, généralement qualifiée de matrice des ondes planes, peut alors faire l'objet d'un traitement, notamment pour la fourniture d'une image de la zone d'intérêt inspectée ou pour l'extraction de paramètres significatifs de défauts de structure dans la zone d'intérêt inspectée. Compte tenu des capacités de calcul actuelles des processeurs, ce traitement peut être embarqué dans les instruments de contrôle pour des traitements en temps réel.
Ce type d'acquisition, généralement qualifié de « plane-wave compounding » ou « plane-wave imaging », est par exemple décrit dans l'article de Montaldo et al, intitulé « Cohérent plane-wave compounding for very high frame rate ultrasonography and transient elastography », publié dans IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control, vol. 56, n° 3, pages 489-506, mars 2009. Il est essentiellement envisagé dans le domaine médical et mis en œuvre dans certains échographes ultra-rapides permettant d'imager des cartographies de l'élasticité du corps humain. L'image résultante est de haute qualité et ne nécessite que quelques dizaines de tirs ultrasonores (généralement M est compris entre 10 et 30) pour un capteur de N = 128 transducteurs d'émission/réception, contre 128 tirs ou presque pour d'autres techniques d'imagerie telles que la focalisation synthétique en tous points basée sur une acquisition de type FMC (de l'anglais « Full Matrix Capture ») ou équivalent. De plus, les algorithmes résultants et les applications associées se prêtent particulièrement bien à une parallélisation des calculs sur des processeurs de type GPU (de l'anglais « Graphie Processing Unit ») implantés dans les cartes graphiques des ordinateurs. Les performances d'un échographe implémentant une technique de « plane-wave compounding » peuvent ainsi atteindre en pratique 10 000 images/s. Un autre atout de cette technique d'acquisition réside dans le fait que chaque tir est réalisé en sollicitant tous les transducteurs d'émission de sorte que l'énergie émise est élevée, rendant cette méthode moins sensible aux phénomènes d'atténuation, de bruit électronique ou de structure.
Ce type d'acquisition est également exploité dans la demande de brevet WO 2015/092250 A1 en lui adaptant astucieusement le principe de focalisation synthétique en tous points, de manière à tirer profit de la simplicité de la technique de « plane-wave compounding » en vue d'atteindre une cadence d'acquisition élevée et de la qualité d'image, en termes de résolution spatiale et de contraste, liée à une focalisation synthétisée en tous points de l'image souhaitée. Cela a permis d'envisager la technique de « plane-wave compounding » pour des applications de contrôle non destructif. Mais quelle que soit la technique d'imagerie choisie, les images obtenues en contrôle non destructif peuvent présenter un niveau de bruit significatif en fonction des propriétés du milieu sondé. Par exemple, il s'agit d'un bruit électronique lorsque le matériau est homogène et viscoélastique, ou d'un bruit de structure lorsque les ondes sont diffusées par des hétérogénéités du matériau. Les images calculées avec les techniques de focalisation en tous points sont aussi affectées par des artefacts, ou échos fantômes, liés à des échos de géométrie, par exemple l'écho d'un fond de pièce à proximité de celui d'un défaut.
Il peut ainsi être souhaité de concevoir un procédé de traitement de signaux ultrasonores qui permette de s'affranchir d'au moins une partie des problèmes et contraintes précités tout en exploitant l'acquisition avantageuse des signaux par émissions successives d'ondes planes.
Il est donc proposé un procédé de traitement de signaux issus d'une acquisition par sondage ultrasonore comportant les étapes suivantes :
- commande de L transducteurs d'émission pour M émissions successives d'ondes ultrasonores planes présentant M angles d'émission successifs différents dans M zones d'émission,
- commande de N transducteurs de réception de manière à recevoir simultanément et pendant une durée prédéterminée, pour chaque émission, N signaux temporels de mesure, mesurant en particulier des échos dus à des réflexions de l'émission considérée, et
- obtention d'une matrice [MR(t)] de signaux temporels ultrasonores de taille NxM, chaque coefficient MRjj(t) de cette matrice représentant le signal de mesure reçu par le i-ième transducteur de réception dû à la j-ième émission, puis
- transformation de la matrice [MR(t)] de signaux temporels en une matrice [FTMR(f)] de signaux fréquentiels, puis décomposition en valeurs singulières de la matrice [FTMR(f)] de signaux fréquentiels,
- élimination d'une partie des valeurs singulières et vecteurs singuliers associés issus de ladite décomposition en valeurs singulières, et
- reconstitution d'une matrice débruitée [MRu(t)] de signaux temporels à partir des valeurs singulières et vecteurs singuliers non éliminés.
Bien que la matrice des ondes planes [MR(t)] ne soit pas de même nature que les matrices de réponses impulsionnelles inter-éléments habituellement obtenues par les techniques classiques d'acquisition exploitant la focalisation synthétique en tous points, il a été constaté de façon inattendue qu'une méthode de filtrage de bruit basée sur une décomposition en valeurs singulières d'une transformée de la matrice des ondes planes dans le domaine fréquentiel fournit des résultats surprenants en termes d'atténuation du bruit. La matrice débruitée [MRu(t)] ainsi obtenue permet notamment une reconstitution d'image de qualité nettement améliorée par rapport à ce que produit la matrice des ondes planes [MR(t)] sans ce traitement.
De façon optionnelle :
- l'élimination d'une partie des valeurs singulières et vecteurs singuliers associés se fait selon un critère prédéterminé de distinction entre valeurs singulières liées à des défauts et valeurs singulières liées à du bruit, et
- la reconstitution de la matrice débruitée [MRu(t)] de signaux temporels comporte une reconstruction d'une matrice débruitée [FTMRu(f)] de signaux fréquentiels à partir des valeurs singulières et vecteurs singuliers non éliminés, puis une transformation inverse de cette matrice débruitée [FTMRu(f)] de signaux fréquentiels en la matrice débruitée [MRu(t)] de signaux temporels.
De façon optionnelle également, la transformation et la transformation inverse sont des transformations de Fourier discrètes.
De façon optionnelle également, l'élimination d'une partie des valeurs singulières et vecteurs singuliers associés se fait par :
- comparaison d'une décroissance d'amplitudes des valeurs singulières avec une courbe théorique de décroissance de valeurs singulières issues d'une matrice de bruit théorique dont les composantes sont des variables aléatoires gaussiennes indépendantes, puis
- suppression des valeurs singulières appartenant, à un coefficient de proportionnalité commun près, à la courbe théorique de décroissance.
De façon optionnelle également, la courbe théorique de décroissance est définie par une fonction réciproque F"1 (1 -a) elle-même définie par une fonction F (a) , dite de répartition de valeurs singulières aléatoires, telle que :
Figure imgf000006_0001
où Π[0;2] désigne la fonction porte sur l'intervalle [0;2].
De façon optionnelle également, un procédé de traitement de signaux ultrasonores selon l'invention peut en outre comporter une reconstitution d'une zone imagée par calcul, en chaque point d'une pluralité de points prédéterminés de cette zone imagée, d'une valeur résultant d'une sommation cohérente de valeurs instantanées prises respectivement par au moins une partie des NxM signaux temporels de la matrice [MRu(t)] à des temps de vol correspondant respectivement à un passage par le point considéré selon un mode de propagation prédéterminé.
De façon optionnelle également, le calcul se fait sur une partie des NxM signaux temporels de la matrice [MRu(t)] dans un secteur angulaire restreint dans l'ensemble des émissions successives, ce secteur angulaire restreint étant sélectionné de telle sorte que les ondes planes qui en sont exclues n'interagissent pas avec au moins un défaut lié aux valeurs singulières et vecteurs singuliers non éliminés.
De façon optionnelle également, le secteur angulaire restreint est sélectionné sur la base d'une comparaison, pour au moins l'une des valeurs singulières non éliminées, d'une valeur de phase expérimentale du vecteur singulier qui lui est associé en émission avec une valeur de phase théorique en présence dudit au moins un défaut lié à ce vecteur singulier.
II est également proposé un programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, comprenant des instructions pour l'exécution des étapes d'un procédé de traitement de signaux ultrasonores selon l'invention, lorsque ledit programme est exécuté sur un ordinateur.
II est également proposé un dispositif de sondage à ultrasons comportant :
- une sonde comprenant L transducteurs d'émission à ultrasons et N transducteurs de réception à ultrasons,
- des moyens de commande des L transducteurs d'émission pour M émissions successives d'ondes ultrasonores planes présentant M angles d'émission successifs différents dans M zones d'émission,
- des moyens de commande des N transducteurs de réception de manière à recevoir simultanément et pendant une durée prédéterminée, pour chaque émission, N signaux temporels de mesure, mesurant en particulier des échos dus à des réflexions de l'émission considérée, et
- un processeur de reconstitution d'une matrice [MR(t)] de signaux temporels ultrasonores de taille NxM, chaque coefficient MR (t) de cette matrice représentant le signal de mesure reçu par le i-ième transducteur de réception dû à la j-ième émission,
le processeur étant en outre configuré pour effectuer les traitements suivants : - transformation de la matrice [MR(t)] de signaux temporels en une matrice [FTMR(f)] de signaux fréquentiels, puis décomposition en valeurs singulières de la matrice [FTMR(f)] de signaux fréquentiels,
- élimination d'une partie des valeurs singulières et vecteurs singuliers associés issus de ladite décomposition en valeurs singulières, et
- reconstitution d'une matrice débruitée [MRu(t)] de signaux temporels à partir des valeurs singulières et vecteurs singuliers non éliminés.
L'invention sera mieux comprise à l'aide de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés dans lesquels :
- la figure 1 représente schématiquement la structure générale d'un dispositif de sondage à ultrasons selon un mode de réalisation de l'invention,
- la figure 2 illustre un principe d'émissions successives d'ondes ultrasonores planes mis en œuvre par le dispositif de la figure 1 ,
- la figure 3 illustre graphiquement le résultat d'une décomposition en valeurs singulières d'une matrice de signaux fréquentiels obtenue par transformation d'une matrice d'ondes planes elle-même obtenue à l'aide d'émissions réalisées selon le principe de la figure 2,
- les figures 4 et 5 illustrent schématiquement et graphiquement le principe d'une sélection de secteur angulaire restreint pour améliorer la détection d'un défaut, selon un mode de réalisation de l'invention,
- la figure 6 illustre trois exemples d'images reconstituées sans ou avec application des principes de l'invention, et
- la figure 7 illustre les étapes successives d'un procédé d'acquisition et de traitement de signaux ultrasonores mis en œuvre par le dispositif de la figure 1 , selon un mode de réalisation de l'invention.
En référence à la figure 1 , un dispositif de sondage 100 d'un objet 102 selon un mode de réalisation de l'invention comporte une sonde à ultrasons 104 présentant un boîtier 106, c'est-à-dire un élément de structure indéformable qui sert de référentiel attaché à la sonde 104, dans lequel sont disposés, par exemple linéairement ou matriciellement, N transducteurs 108 ; ... ,108N fixes ou mobiles disposés en réseau.
L'objet 102 est par exemple une pièce mécanique que l'on souhaite examiner par contrôle non destructif ou bien, dans un contexte médical, une partie de corps humain ou animal que l'on souhaite contrôler de manière non invasive. Dans le mode de réalisation de la figure 1 , l'objet 102 est immergé dans un liquide, tel que de l'eau 1 10, et la sonde 104 est maintenue à distance de l'objet 102 afin que l'eau 1 10 les sépare. Mais dans un autre mode de réalisation équivalent, la sonde 104 pourrait être au contact direct de l'objet 102.
Les transducteurs 108 ; ... ,108N sont conçus pour émettre individuellement des ondes ultrasonores en direction de l'objet 102 en réponse à des signaux de commande identifiés sous la référence générale C, selon des directions principales parallèles les unes aux autres, indiquées par des flèches en pointillés sur la figure 1 , et dans un plan principal qui est celui de la figure.
Les transducteurs 108 ; ... ,108N sont en outre conçus pour détecter des échos des ondes ultrasonores se réfléchissant sur et dans l'objet 102 et pour fournir des signaux de mesure identifiés sous la référence générale S et correspondant à ces échos. Ainsi, dans l'exemple non limitatif de la figure 1 , les transducteurs 108 ; ... ,108N remplissent à la fois les fonctions d'émission et de réception, mais des récepteurs différents des émetteurs pourraient également être prévus dans des boîtiers différents et indépendants tout en restant conformes avec les principes de l'invention. De plus, le nombre L d'émetteurs pourrait tout à fait être différent du nombre N de récepteurs.
Le dispositif de sondage 100 comporte en outre un circuit électronique 1 12 de commande des transducteurs 108i , ...,108N de la sonde 104 et de traitement des signaux de mesure S. Ce circuit électronique 1 12 est connecté à la sonde 104 afin de lui transmettre les signaux de commande C et afin de recevoir les signaux de mesure S. Le circuit électronique 1 12 est par exemple celui d'un ordinateur. Il présente une unité centrale de traitement 1 14, telle qu'un microprocesseur conçu pour émettre vers la sonde 104 les signaux de commande C et pour recevoir de la sonde 104 les signaux de mesure S, et une mémoire 1 16 dans laquelle est notamment enregistré un programme d'ordinateur 1 18.
Le programme d'ordinateur 1 18 comporte tout d'abord des instructions 120 pour générer les signaux C de commande des transducteurs 108 ; ... ,108N et recevoir leurs échos. Ces instructions sont plus précisément programmées de manière à :
- activer les L = N transducteurs 108i , ... ,108N en tant qu'émetteurs pour M émissions successives d'ondes ultrasonores planes présentant M angles d'émission successifs différents dans M zones d'émission de l'objet 102, - activer les transducteurs 108i , ... ,108N en tant que récepteurs pour, suite à chaque émission, recevoir simultanément, par ces N récepteurs et pendant une durée prédéterminée de la profondeur d'inspection souhaitée, N signaux temporels de mesure, mesurant en particulier des échos dus à des réflexions de chaque émission considérée.
Les ondes ultrasonores planes sont obtenues à l'émission en appliquant aux transducteurs 108 ; ... ,108N des lois de retards enregistrées en mémoire 1 16 dans une base 122 de lois de retards. Chaque loi de retards définit des retards à appliquer aux transducteurs 108 ; ... ,108N en émission, de manière à engendrer une onde ultrasonore plane à un angle d'émission souhaité parmi les M angles d'émission successifs différents. Il est donc prévu autant de lois de retards que d'émissions successives souhaitées.
Comme illustré sur la figure 2 dans un cas où le nombre M d'émissions successives est impair et où les angles d'émissions se succèdent avec un pas constant dans un secteur angulaire symétrique par rapport à la direction z orthogonale au réseau de transducteurs I O81, ... ,108N, la première émission d'onde plane est associée à une loi de retards ΤΊ portant sur des signaux émis par les transducteurs I O81, ... ,108N, permettant l'émission d'une onde plane d'angle d'émission Θ1 par rapport à la direction z dans une première zone d'émission ZEi partiellement située en dehors de l'ouverture de la sonde 104. La (M+1 )/2-ième émission d'onde plane est associée à une loi de retards T(M+i)/2 uniforme pour l'émission d'une onde plane d'angle d'émission nul par rapport à la direction z dans une (M+1 )/2-ième zone d'émission ZE(M+1)/2 couvrant l'ouverture de la sonde 104. Enfin, la dernière émission d'onde plane est associée à une loi de retards TM permettant l'émission d'une onde plane d'angle d'émission ΘΜ = - par rapport à la direction z dans une dernière zone d'émission ZEM partiellement située en dehors de l'ouverture de la sonde 104. D'une façon générale, la j-ième émission d'onde plane est associée à une loi de retards T permettant l'émission d'une onde plane d'angle d'émission 0j = θ + (j - 1 ).( ΘΜ - Θ^Μ - 1 ) par rapport à la direction z. En réalité, le plus souvent M est pair : il n'y a alors pas d'émission à 0° et le pas angulaire n'est pas constant.
Compte tenu de la technique d'acquisition employée, la zone à imager doit être contenue dans l'union des M zones d'émissions successives. Il en résulte que cette zone peut s'étendre au-delà de l'ouverture de la sonde 104, comme cela est visible sur la figure 2. En particulier, la zone imagée peut prendre la forme d'une zone sectorielle délimitée par les extrémités des zones d'émission d'angles maximal et minimal. On peut ainsi obtenir une image de type S-scan.
En variante, et comme avantageusement rendu possible par le principe d'acquisition des signaux par émissions successives d'ondes planes, les M angles d'émission successifs différents à ΘΜ peuvent être définis autour d'une direction moyenne θ(Μ+ΐ)/2 non perpendiculaire au réseau de transducteurs 108 ; ... ,108N. En particulier, lorsqu'il s'agit de détecter des défauts tels qu'une fissure disposée au fond d'un objet à inspecter en contrôle non destructif, cette fissure étant en outre perpendiculaire au réseau de transducteurs, il est préférable de décaler latéralement la zone à inspecter par rapport à la sonde 104 et d'émettre autour d'une moyenne de 45° par exemple. La zone à inspecter peut même être décalée au point de sortir complètement de l'ouverture de la sonde 104.
De manière à améliorer la qualité des signaux de mesure exploités pour reconstituer la zone imagée, il est en outre possible d'appliquer une apodisation des signaux ultrasonores émis par les transducteurs 108i , ... ,108N pour former une onde ultrasonore plane de meilleure qualité, sans distorsion subie à cause des effets de bords. Une telle apodisation est réalisée à l'occasion de chaque émission spatialement sur l'ensemble des transducteurs à l'aide d'une fenêtre d'apodisation telle qu'une loi d'amplitude trapézoïdale, de Hamming ou de Blackman-Harris. Elle a pour résultat de fournir une meilleure définition des zones d'émission successives.
En référence de nouveau à la figure 1 , à la réception des signaux résultant de chacune des M émissions successives, l'ensemble S des NxM signaux temporels de mesure reçus par les N transducteurs 108 ; ... ,108N est renvoyé par la sonde 104 à l'unité centrale de traitement 1 14.
Le programme d'ordinateur 1 18 comporte alors en outre des instructions 124 pour construire une matrice [MR(t)] de signaux temporels ultrasonores de taille NxM, qualifiée de matrice des ondes planes. Chaque coefficient MR,j(t) de cette matrice représente le signal de mesure reçu par le transducteur 108, en réponse à la j-ième émission.
De façon optionnelle, le programme d'ordinateur 1 18 comporte en outre des instructions 126 pour effectuer un filtrage temporel de la matrice [MR(t)], ce filtrage visant à supprimer toute information se trouvant à des temps de vol exclus de la zone d'intérêt dans l'objet 102.
Le programme d'ordinateur 1 18 comporte en outre des instructions 128 pour transformer la matrice [MR(t)] en une matrice [FTMR(f)] de signaux fréquentiels par transformée de Fourier, avantageusement par transformée de Fourier discrète après échantillonnage temporel des signaux ultrasonores formant les coefficients de la matrice [MR(t)], ou, encore plus avantageusement encore, par calcul de FFT (de l'anglais « Fast Fourier Transform ») si le nombre d'échantillons de chaque coefficient de la matrice [MR(t)] le permet.
Le programme d'ordinateur 1 18 comporte en outre des instructions 130 pour décomposer la matrice [FTMR(f)] de signaux fréquentiels en valeurs singulières sur une bande de fréquences. Bien qu'il soit connu de décomposer en valeurs singulières une matrice de réponses impulsionnelles inter-éléments habituellement obtenues par les techniques classiques d'acquisition exploitant la focalisation synthétique en tous points, cette opération n'est pas équivalente lorsqu'elle est appliquée sur une matrice telle que la matrice [FTMR(f)]. En effet, la matrice notée [PFTMR(f)] de dimensions MxM et définie par le produit [FTMR(f)]†.[FTMR(f)], où «† » est le symbole de la transposée conjuguée d'une matrice, ne représente pas l'opérateur de retournement temporel en émission comme pour la matrice des réponses impulsionnelles interéléments. L'interprétation physique des valeurs singulières et vecteurs singuliers n'est alors pas la même que dans la méthode de Décomposition de l'Opérateur de Retournement Temporel (méthode DORT) telle qu'enseignée dans l'article de Prada et al, intitulé « Eigenmodes of the time reversai operator: a solution to sélective multiple-target média », publié dans Wave Motion 20, pages 151 -163 (1994). De plus, la matrice [FTMR(f)] n'est généralement ni carrée ni symétrique, M étant souvent même nettement inférieur à N en acquisition par émission d'ondes planes. De la sorte, une décomposition en valeurs singulières ne produit pas de diagonalisation comme dans la plupart des travaux portant sur la matrice des réponses impulsionnelles inter-éléments.
Plus précisément, l'opération de décomposition en valeurs singulières permet d'estimer les matrices U(f), S(f) et V(f) telles que :
[FTMR(f)] = u{f)s{f)v {f) =∑ai {f).ul {f).v {f) , où U(f) = [ui (f) , .... uN(f)] et i=\
V=[v! (f), vM(f)] sont des matrices unitaires orthogonales de tailles respectives NxN et MxM, qui contiennent respectivement les vecteurs singuliers en réception et émission, où S est une matrice réelle de taille NxM contenant les K = min(N;M) valeurs singulières a,(f) de la matrice [FTMR(f)], ordonnées de façon décroissante à une fréquence f donnée σι (ί) ≥ ... ≥ aK(f) ≥ 0, et où min(N;M) est la fonction qui retourne la valeur minimale entre N et M (généralement c'est M). Le programme d'ordinateur 1 18 comporte en outre des instructions 132 pour réduire le rang de la matrice [FTMR(f)], en éliminant une partie des valeurs singulières a,(f). Cette élimination se fait selon un critère de distinction entre valeurs singulières liées à des défauts et valeurs singulières liées à du bruit, les premières étant d'amplitudes supérieures aux secondes. Compte tenu du fait que a^f)≥ ...≥ aK(f)≥ 0, il s'agit de trouver la valeur entière P(f) comprise entre 1 et K telle que a^f), aP(f)(f) peuvent être considérées comme liées à des défauts à détecter dans l'objet 102 et aP(f)+1(f) , aK(f) peuvent être éliminées car considérées comme liées à du bruit. Dans le cas de petits défauts idéalement espacés les uns des autres, P(f) est égal au nombre de défauts présents dans la zone d'intérêt inspectée. On peut alors écrire la matrice [FTMR(f)] sous la forme d'une addition de deux matrices [FTMRu(f)] et [FTMRN(f)] :
[FTMR(f)] = [FTMRU (/)] + [FTMRN (/)] , avec
Pif)
[FTMRU (/)] = ^σ,. (/).«,. (/).v,†(/) , appelée matrice de signal utile, et
i=l [FTMRN ( f)] = ∑ ; ( ( ).ν; ( ) , appelée matrice de bruit.
;=p(/)+i
Réduire le rang de la matrice [FTMR(f)] revient ainsi à ne conserver que [FTMRu(f)].
En pratique, la détermination de P(f) se fait par une étude de la courbe de décroissance des amplitudes des valeurs singulières. Cette étude peut se faire sur la base de la théorie des matrices aléatoires telle qu'enseignée dans :
- l'article de Marcenko et al, intitulé « Distribution of eigenvalues for some sets of random matrices », publié dans Mathematics of the USSR-Sbornik, vol. 1 , n°4, pages 457-483 (1967), et
- l'article de Aubry et al, intitulé « Détection and imaging in a random médium: a matrix method to overcome multiple scattering and aberration », publié dans Journal of Applied Physics, 106(4), 044903 (2009). En supposant que toutes les composantes de la matrice de bruit [FTMRN(f)] sont des variables aléatoires gaussiennes indépendantes, il peut en effet être montré que les valeurs singulières liées au bruit σΡ(ί)+ι (f) , aK(f) appartiennent, à un coefficient de proportionnalité commun près, à une courbe théorique définie par une fonction réciproque F"1(1 -a) elle-même définie par une fonction F(a), dite de répartition de valeurs singulières aléatoires, avec : Fia) = i + 2 arcsin ( )] Π[0.2] ,
où Π[0;2] désigne la fonction porte sur l'intervalle [0;2].
La fonction F (a) donne des valeurs comprises entre 0 et 1 sur l'intervalle support σ e [0;2], où σ est la variable de valeur singulière. Elle est par ailleurs strictement croissante. Sa fonction réciproque F"1 (a) est donc également strictement croissante, de sorte que la fonction F~1 (1 -a) donne la courbe de décroissance des valeurs singulières recherchée, à une constante de proportionnalité près pour l'ajuster à la courbe expérimentale.
Ainsi, le nombre de valeurs singulières calculées par exécution des instructions 1 30 n'appartenant pas à cette courbe théorique définit le rang P(f) de la matrice [FTMRu(f)] de signal utile. A titre d'exemple, la figure 3 illustre en unité arbitraire a.u. (de l'anglais « arbitrary unit ») une distribution de valeurs singulières calculées à une fréquence f = 5 MHz, pour une acquisition réalisée avec un échantillon de polyéthylène en émettant M = 64 ondes planes entre -31 ,5° et +31 ,5°. Deux valeurs singulières se détachent de la courbe théorique, dont une assez nettement, indiquant dans cet exemple que P(f) = 2 à 5 MHz.
Pour calculer la matrice [FTMRu(f)] de signal utile sur plusieurs fréquences f, il suffit d'identifier à chaque fréquence souhaitée le nombre de valeurs singulières n'appartenant pas à la courbe théorique.
La matrice [FTMRu(f)] ainsi reconstituée est une matrice débruitée de signaux fréquentiels, le sous-espace de bruit représenté par la matrice [FTMRN(f)] ayant été éliminé.
Le programme d'ordinateur 1 1 8 comporte en outre des instructions 134 pour transformer la matrice [FTMRu(f)] en une matrice débruitée [MRu(t)] de signaux temporels par transformée de Fourier inverse, avantageusement par transformée de Fourier discrète inverse, ou, encore plus avantageusement encore, par calcul de IFFT (de l'anglais « Inverse Fast Fourier Transform ») si le nombre d'échantillons de chaque coefficient de la matrice [FTMRu(f)] le permet.
Enfin, le programme d'ordinateur 1 18 comporte des instructions, désignées par la référence générale 136, de traitement de la matrice [MRu(t)]. Le traitement réalisé par les instructions 1 36 peut inclure une reconstitution d'image numérique de la zone d'intérêt dans l'objet 102 par adaptation du principe de focalisation synthétique en tous points, comme enseigné par exemple dans le document WO 2015/092250 A1 . On reconstitue ainsi une image numérique de la zone d'intérêt dont la qualité est meilleure que si la reconstitution avait été réalisée sur la matrice non débruitée [MR(t)]. Notamment, le Rapport Signal sur Bruit (RSB) est amélioré. En variante ou en complément, le traitement réalisé par les instructions 136 pourrait inclure une focalisation adaptative et sélective.
Selon une adaptation relativement simple mais pas toujours optimale de la focalisation synthétique en tous points, en chaque pixel P de l'image numérique reconstituée, on calcule le module A(P) d'une sommation cohérente impliquant les NxM signaux temporels de la matrice [MRu(t)] à NxM temps de vols calculés selon un mode de propagation prédéterminé, chaque temps de vol t étant le temps mis par la j-ième onde plane pour être reçue par le i-ième transducteur de réception en passant par le pixel considéré selon le mode de propagation prédéterminé :
Figure imgf000015_0001
où μ (Ρ) et pj(P) sont des coefficients de pondération respectivement en émission et en réception dont les expressions dépendent de l'application considérée pour prendre en compte des phénomènes ou traitements tels qu'un filtrage d'échos de géométrie, une compensation d'atténuation due à un étalement spatial des ondes, etc.
Ce mode de calcul n'est souvent pas optimal parce que, comme l'illustre la figure 4, toutes les ondes planes émises par la sonde 104 n'interagissent pas avec le ou les défauts D liés aux valeurs singulières de la matrice de signal utile [FTMRu(f)], c'est-à-dire les valeurs singulières a^f), aP(f)(f). Il peut alors être avantageux de sélectionner un secteur angulaire restreint [9m1 ;9m2] dans l'ensemble des émissions successives à l'extérieur duquel les ondes planes n'interagissent pas avec le ou les défauts précités. Ce secteur angulaire retreint est défini par ses indices d'émissions minimal (m ) et maximal (m2) où ~\≤m†≤m2≤M. Le module A(P) est alors de préférence calculé par l'addition cohérente suivante :
) μ (Ρ)ρί(Ρ) «[/ [ty(P)] . Dans le cas d'un défaut assimilable à un réflecteur ponctuel, une méthode astucieuse pour déterminer les valeurs de ΙΤΗ et m2 va maintenant être détaillée en référence à la figure 5. Elle consiste à comparer, pour au moins l'une des valeurs singulières σι (ί), aP(f)(f), et notamment au moins pour la première valeur singulière a^f) qui est la plus élevée, une valeur de phase expérimentale du vecteur singulier qui lui est associé en émission avec une valeur de phase théorique en présence du ou des défauts D. Cette comparaison est réalisée à une fréquence fc choisie qui peut être la fréquence centrale de fonctionnement de la sonde 104, une fréquence pour laquelle la valeur singulière ai (f) prend sa valeur la plus élevée, ou toute autre fréquence prédéterminée. Sur la figure 5, la comparaison est ainsi réalisée à 5 MHz pour le vecteur singulier en émission v^f) correspondant à un emplacement de défaut(s) D de coordonnées (XD,ZD) .
Dans les conditions expérimentales indiquées ci-dessus, la valeur de phase théorique est une phase corrigée calculée de la façon suivante pour être toujours négative :
2nfc(XD sin 9j + ZD cos
Figure imgf000016_0001
avec c la célérité de l'onde plane dans le milieu considéré et 6j son angle d'émission.
Les coordonnées (XD,ZD) peuvent être déterminées avec la matrice de signal utile [FTMRu(f)] en calculant la rétro-propagation du vecteur singulier en réception Ui (f) à la fréquence fc, comme par exemple enseigné dans l'article de Lopez Villaverde et al, intitulé « Ultrasonic imaging of defects in coarse-grained steels with the décomposition of the time reversai operator », publié dans Journal of the Acoustical Society of America, volume 140, n°1 , pages 541 -550 (2016) .
La définition ci-dessus de la phase théorique corrigée donne la courbe en trait continu de la figure 5.
La valeur de phase expérimentale du vecteur singulier en émission v^y est également une phase corrigée calculée de la façon suivante pour être toujours négative :
Vy, 1 < y≤ M, <pj1( c) = arg[i 1( c)] . - max [arg[¾ e)] .] .
Cette définition de la phase corrigée du vecteur singulier v^y donne les M valeurs représentées par des petits cercles sur la figure 5.
On s'aperçoit que les valeurs expérimentales et théoriques sont en bonne correspondance à l'intérieur d'un secteur angulaire restreint dont les extrémités donnent les valeurs de ΙΤΗ et m2. En dehors de ce secteur angulaire restreint, les valeurs expérimentales s'éloignent très nettement des valeurs théoriques.
D'une façon générale, si la comparaison est effectuée pour plusieurs valeurs singulières liées au(x) défaut(s) D, les valeurs finalement retenues pour ΙΤΗ et m2 peuvent être respectivement le minimum et le maximum des valeurs trouvées pour chacune des valeurs singulières de la matrice de signal utile. Le mode de calcul avantageux détaillé ci-dessus peut également être combiné avec une adaptation du principe de focalisation synthétique en tous points telle qu'enseignée dans le document WO 2015/092250 A1 .
La figure 6 illustre, dans un exemple d'objet sondé comportant un défaut D circulaire central, trois images obtenues :
- par une acquisition FMC avec reconstitution de l'image par focalisation synthétique en tous points (image de gauche),
- par une acquisition et un traitement tels qu'enseignés dans le document WO 2015/092250 A1 (image au centre), et
- par une acquisition et un traitement de filtrage et de sélection d'un secteur angulaire restreint tels qu'enseignés conformément aux principes de la présente invention (image de droite).
On note une forte atténuation du bruit (gain de 20 dB notamment entre l'image au centre et celle de droite en termes de RSB) et une meilleure visibilité du défaut.
En référence à la figure 7, un exemple de procédé 700 d'acquisition et de traitement de signaux ultrasonores que peut mettre en œuvre le dispositif 100 de la figure 1 va à présent être décrit selon un mode de réalisation préféré de l'invention.
Au cours d'une étape 702, l'unité de traitement 1 14 exécutant les instructions 120 commande les séquences d'émissions et de réceptions des transducteurs 108i , ... ,108N pour l'acquisition des signaux de mesure MRjj(t) de la matrice [MR(t)].
Ces séquences sont au nombre de M, nombre entier pouvant être très inférieur au nombre N de transducteurs 108i , ...,108N. Après chaque tir, les signaux sont reçus sur l'ensemble des N transducteurs, numérisés et transmis au circuit électronique 1 12.
Au cours d'une étape 704, l'unité de traitement 1 14 exécutant les instructions
124 enregistre les signaux de mesure MRjj(t), ces signaux étant numérisés pour faciliter leur traitement ultérieur. Les étapes 702 et 704 peuvent être exécutées simultanément, c'est-à-dire qu'il n'est pas nécessaire d'attendre que tous les tirs soient réalisés pour commencer à enregistrer les signaux de mesure et effectuer un traitement tel qu'une reconstitution d'image.
Au cours d'une étape 706 optionnelle, l'unité de traitement 1 14 exécutant les instructions 126 effectue un filtrage temporel de la matrice [MR(t)], ce filtrage visant à supprimer toute information se trouvant à des temps de vol exclus de la zone d'intérêt. Cette étape 706 a pour but de faciliter ensuite la séparation des deux sous- espaces représentés par les matrices [FTMRu(f)] et [FTMRN(f)], en particulier lorsque les défauts à imager sont proches d'une interface fortement échogène, tel un fond de pièce. Elle permet de limiter la zone à imager à un voisinage proche des défauts en excluant notamment les interfaces échogènes perturbatrices. Elle trouve tout son intérêt dans l'imagerie de fissures se formant depuis le fond de l'objet.
Au cours d'une étape 708, l'unité de traitement 1 14 exécutant les instructions
128 effectue une transformée de Fourier discrète de la matrice [MR(t)] pour obtenir la matrice [FTMR(f)] de signaux fréquentiels.
Au cours d'une étape 710, l'unité de traitement 1 14 exécutant les instructions 130 réalise une décomposition en valeurs singulières de la matrice [FTMR(f)], comme détaillé précédemment.
Au cours d'une étape 712, l'unité de traitement 1 14 exécutant les instructions 132 réduit le rang de la matrice [FTMR(f)] en ne conservant que la matrice de signal utile [FTMRu(f)].
Au cours d'une étape 714, l'unité de traitement 1 14 exécutant les instructions 134 effectue une transformée de Fourier inverse discrète de la matrice [FTMRu(f)] pour obtenir la matrice débruitée [MRu(t)] de signaux temporels.
Au cours d'une étape 716, l'unité de traitement 1 14 exécutant les instructions 136 sélectionne, de façon optionnelle mais avantageuse, un secteur angulaire restreint dans l'ensemble des émissions successives à l'extérieur duquel les ondes planes n'interagissent pas avec le ou les défauts à détecter. Ce secteur angulaire retreint est défini par ses indices d'émissions minimal (ITH) et maximal (m2) par exemple selon la méthode détaillée précédemment.
Enfin, au cours d'une dernière étape 718, l'unité de traitement 1 14 exécutant toujours les instructions 136 reconstitue et affiche une image numérique de la zone effective d'intérêt par adaptation du principe de focalisation synthétique en tous points à partir de la matrice débruitée [MRu(t)] dans le secteur angulaire restreint sélectionné.
Il apparaît clairement qu'un dispositif de sondage à ultrasons tel que celui décrit précédemment permet d'obtenir une détection de défaut(s) débruitée et de très bonne qualité, meilleure que celles obtenues par les acquisitions FMC avec reconstitutions d'images par focalisation synthétique en tous points ou par le traitement tel qu'enseigné dans le document le document WO 2015/092250 A1 . De plus, l'enseignement de la présente demande de brevet peut être avantageusement combiné à celui du document WO 2015/092250 A1 . On notera par ailleurs que l'invention n'est pas limitée au mode de réalisation décrit précédemment. Il apparaîtra en effet à l'homme de l'art que diverses modifications peuvent être apportées au mode de réalisation décrit ci-dessus, à la lumière de l'enseignement qui vient de lui être divulgué.
En particulier, les instructions de programme d'ordinateur pourraient être remplacées par des circuits électroniques dédiés aux fonctions réalisées lors de l'exécution de ces instructions.
D'une façon générale, dans les revendications qui suivent, les termes utilisés ne doivent pas être interprétés comme limitant les revendications au mode de réalisation exposé dans la présente description, mais doivent être interprétés pour y inclure tous les équivalents que les revendications visent à couvrir du fait de leur formulation et dont la prévision est à la portée de l'homme de l'art en appliquant ses connaissances générales à la mise en œuvre de l'enseignement qui vient de lui être divulgué.

Claims

REVENDICATIONS
1 . Procédé (700) de traitement de signaux (S) issus d'une acquisition par sondage ultrasonore comportant les étapes suivantes :
commande (702) de L transducteurs d'émission (1081 ; 108N) pour M émissions successives d'ondes ultrasonores planes présentant M angles d'émission successifs différents (θ1 ; ... , ΘΜ) dans M zones d'émission (ZE1 ; ZEM),
commande (702) de N transducteurs de réception (1081 ; 108N) de manière à recevoir (704) simultanément et pendant une durée prédéterminée, pour chaque émission, N signaux temporels de mesure, mesurant en particulier des échos dus à des réflexions de l'émission considérée,
obtention (704) d'une matrice [MR(t)] de signaux temporels ultrasonores de taille NxM, chaque coefficient MRij(t) de cette matrice représentant le signal de mesure reçu par le i-ième transducteur de réception dû à la j-ième émission,
caractérisé en ce qu'il comporte en outre les étapes suivantes :
transformation (708) de la matrice [MR(t)] de signaux temporels en une matrice [FTMR(f)] de signaux fréquentiels, puis décomposition en valeurs singulières (710) de la matrice [FTMR(f)] de signaux fréquentiels,
élimination (712) d'une partie des valeurs singulières et vecteurs singuliers associés issus de ladite décomposition en valeurs singulières (710), et
reconstitution (712, 714) d'une matrice débruitée [MRu(t)] de signaux temporels à partir des valeurs singulières et vecteurs singuliers non éliminés.
2. Procédé (700) de traitement de signaux ultrasonores (S) selon la revendication 1 , dans lequel :
l'élimination (712) d'une partie des valeurs singulières et vecteurs singuliers associés se fait selon un critère prédéterminé de distinction entre valeurs singulières liées à des défauts et valeurs singulières liées à du bruit, et la reconstitution (712, 714) de la matrice débruitée [MRu(t)] de signaux temporels comporte une reconstruction (712) d'une matrice débruitée [FTMRu(f)] de signaux fréquentiels à partir des valeurs singulières et vecteurs singuliers non éliminés, puis une transformation inverse (714) de cette matrice débruitée [FTMRu(f)] de signaux fréquentiels en la matrice débruitée [MRu(t)] de signaux temporels.
3. Procédé (700) de traitement de signaux ultrasonores (S) selon la revendication 2, dans lequel la transformation (708) et la transformation inverse (714) sont des transformations de Fourier discrètes.
4. Procédé (700) de traitement de signaux ultrasonores (S) selon la revendication 2 ou 3, dans lequel l'élimination (712) d'une partie des valeurs singulières et vecteurs singuliers associés se fait par :
comparaison d'une décroissance d'amplitudes des valeurs singulières avec une courbe théorique de décroissance de valeurs singulières issues d'une matrice de bruit théorique dont les composantes sont des variables aléatoires gaussiennes indépendantes, puis
suppression des valeurs singulières appartenant, à un coefficient de proportionnalité commun près, à la courbe théorique de décroissance.
5. Procédé (700) de traitement de signaux ultrasonores (S) selon la revendication 4, dans lequel la courbe théorique de décroissance est définie par une fonction réciproque F"1 (1 -a) elle-même définie par une fonction F(a), dite de répartition de valeurs singulières aléatoires, telle que :
Figure imgf000021_0001
où Π[0;2] désigne la fonction porte sur l'intervalle [0;2].
6. Procédé (700) de traitement de signaux ultrasonores (S) selon l'une quelconque des revendications 1 à 5, comportant en outre une reconstitution (716, 718) d'une zone imagée par calcul, en chaque point d'une pluralité de points prédéterminés de cette zone imagée, d'une valeur résultant d'une sommation cohérente de valeurs instantanées prises respectivement par au moins une partie des NxM signaux temporels de la matrice [MRu(t)] à des temps de vol correspondant respectivement à un passage par le point considéré selon un mode de propagation prédéterminé.
7. Procédé (700) de traitement de signaux ultrasonores (S) selon la revendication 6, dans lequel le calcul se fait sur une partie des NxM signaux temporels de la matrice [MRu(t)] dans un secteur angulaire restreint (m1 ; m2) dans l'ensemble des émissions successives.
8. Procédé (700) de traitement de signaux ultrasonores (S) selon la revendication 7, dans lequel le secteur angulaire restreint (m , m2) est sélectionné (716) sur la base d'une comparaison, pour au moins l'une des valeurs singulières non éliminées, d'une valeur de phase expérimentale du vecteur singulier qui lui est associé en émission avec une valeur de phase théorique en présence dudit au moins un défaut lié à ce vecteur singulier, de telle sorte que les ondes planes qui en sont exclues n'interagissent pas avec au moins un défaut lié aux valeurs singulières et vecteurs singuliers non éliminés.
9. Programme d'ordinateur (1 18) téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur (1 14), caractérisé en ce qu'il comprend des instructions (120, 124, 126, 128, 130, 132, 134, 136) pour l'exécution des étapes d'un procédé (700) de traitement de signaux ultrasonores (S) selon l'une quelconque des revendications 1 à 8, lorsque ledit programme est exécuté sur un ordinateur (1 12).
10. Dispositif (100) de sondage à ultrasons comportant :
une sonde (104) comprenant L transducteurs d'émission à ultrasons (108i , 108N) et N transducteurs de réception à ultrasons (108i , 108N),
des moyens (1 12) de commande des L transducteurs d'émission pour M émissions successives d'ondes ultrasonores planes présentant M angles d'émission successifs différents (θ1 ; ... , ΘΜ) dans M zones d'émission (ZE1 ; ZEM),
des moyens (1 12) de commande des N transducteurs de réception de manière à recevoir simultanément et pendant une durée prédéterminée, pour chaque émission, N signaux temporels de mesure, mesurant en particulier des échos dus à des réflexions de l'émission considérée, et
un processeur (1 14) de reconstitution d'une matrice [MR(t)] de signaux temporels ultrasonores de taille NxM, chaque coefficient MRij(t) de cette matrice représentant le signal de mesure reçu par le i-ième transducteur de réception dû à la j-ième émission, caractérisé en ce que le processeur (1 14) est en outre configuré pour effectuer les traitements suivants :
transformation de la matrice [MR(t)] de signaux temporels en une matrice [FTMR(f)] de signaux fréquentiels, puis décomposition en valeurs singulières de la matrice [FTMR(f)] de signaux fréquentiels élimination d'une partie des valeurs singulières et vecteurs singuliers associés issus de ladite décomposition en valeurs singulières, et reconstitution d'une matrice débruitée [MRu(t)] de signaux temporels à partir des valeurs singulières et vecteurs singuliers non éliminés.
PCT/FR2017/053391 2016-12-15 2017-12-05 Procede de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage a ultrasons correspondants WO2018109314A1 (fr)

Priority Applications (4)

Application Number Priority Date Filing Date Title
EP17821680.0A EP3555611A1 (fr) 2016-12-15 2017-12-05 Procede de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage a ultrasons correspondants
JP2019532089A JP2020502518A (ja) 2016-12-15 2017-12-05 超音波探査による取得に由来する信号を処理する方法、対応するコンピュータプログラムおよび超音波探査装置
CA3046106A CA3046106A1 (fr) 2016-12-15 2017-12-05 Procede de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage a ultrasons correspondants
US16/469,918 US11054401B2 (en) 2016-12-15 2017-12-05 Method for processing signals coming from an acquisition by ultrasound probing, corresponding computer program and ultrasound probing device

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1662525 2016-12-15
FR1662525A FR3060753B1 (fr) 2016-12-15 2016-12-15 Procede de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage a ultrasons correspondants

Publications (1)

Publication Number Publication Date
WO2018109314A1 true WO2018109314A1 (fr) 2018-06-21

Family

ID=58609497

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/FR2017/053391 WO2018109314A1 (fr) 2016-12-15 2017-12-05 Procede de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage a ultrasons correspondants

Country Status (6)

Country Link
US (1) US11054401B2 (fr)
EP (1) EP3555611A1 (fr)
JP (1) JP2020502518A (fr)
CA (1) CA3046106A1 (fr)
FR (1) FR3060753B1 (fr)
WO (1) WO2018109314A1 (fr)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010142927A1 (fr) * 2009-06-11 2010-12-16 Centre National De La Recherche Scientifique Procede et dispositif ultrasonores pour caracteriser un milieu
WO2014009671A1 (fr) * 2012-07-12 2014-01-16 Commissariat à l'énergie atomique et aux énergies alternatives Procédé de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage à ultrasons correspondants
WO2015092250A1 (fr) 2013-12-20 2015-06-25 Commissariat à l'énergie atomique et aux énergies alternatives Procédé de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage à ultrasons correspondants

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010142927A1 (fr) * 2009-06-11 2010-12-16 Centre National De La Recherche Scientifique Procede et dispositif ultrasonores pour caracteriser un milieu
WO2014009671A1 (fr) * 2012-07-12 2014-01-16 Commissariat à l'énergie atomique et aux énergies alternatives Procédé de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage à ultrasons correspondants
WO2015092250A1 (fr) 2013-12-20 2015-06-25 Commissariat à l'énergie atomique et aux énergies alternatives Procédé de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage à ultrasons correspondants

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
AUBRY ET AL.: "Détection and imaging in a random médium: a matrix method to overcome multiple scattering and aberration", JOURNAL OF APPLIED PHYSICS, vol. 106, no. 4, 2009, pages 044903, XP012123911, DOI: doi:10.1063/1.3200962
DEYLAMI ALI MOHADES ET AL: "An improved minimum variance beamforming applied to plane-wave imaging in medical ultrasound", 2016 IEEE INTERNATIONAL ULTRASONICS SYMPOSIUM (IUS), IEEE, 18 September 2016 (2016-09-18), pages 1 - 4, XP032988777, DOI: 10.1109/ULTSYM.2016.7728895 *
ELODIE TIRAN ET AL: "Multiplane wave imaging increases signal-to-noise ratio in ultrafast ultrasound imaging", PHYSICS IN MEDICINE AND BIOLOGY, INSTITUTE OF PHYSICS PUBLISHING, BRISTOL GB, vol. 60, no. 21, 21 October 2015 (2015-10-21), pages 8549 - 8566, XP020290386, ISSN: 0031-9155, [retrieved on 20151021], DOI: 10.1088/0031-9155/60/21/8549 *
LOPEZ VILLAVERDE ET AL.: "Ultrasonic imaging of defects in coarse-grained steels with the décomposition of the time reversai operator", JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA, vol. 140, no. 1, 2016, pages 541 - 550, XP012209490, DOI: doi:10.1121/1.4958683
MARCENKO ET AL.: "Distribution of eigenvalues for some sets of random matrices", MATHEMATICS OF THE USSR-SBORNIK, vol. 1, no. 4, 1967, pages 457 - 483
MONTALDO ET AL.: "Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography", IEEE TRANSACTIONS ON ULTRASONICS FERROELECTRICS AND FREQUENCY CONTROL, vol. 56, no. 3, March 2009 (2009-03-01), pages 489 - 506, XP011255897
PRADA ET AL.: "Eigenmodes of the time reversai operator: a solution to selective multiple-target média", WAVE MOTION, vol. 20, 1994, pages 151 - 163
YING LI ET AL: "Plane-wave imaging using synthetic aperture imaging reconstruction technique with regularized singular-value decomposition (RSVD)", 2016 IEEE INTERNATIONAL ULTRASONICS SYMPOSIUM (IUS), IEEE, 18 September 2016 (2016-09-18), pages 1 - 3, XP032988813, DOI: 10.1109/ULTSYM.2016.7728888 *

Also Published As

Publication number Publication date
FR3060753A1 (fr) 2018-06-22
US11054401B2 (en) 2021-07-06
US20190317054A1 (en) 2019-10-17
CA3046106A1 (fr) 2018-06-21
EP3555611A1 (fr) 2019-10-23
JP2020502518A (ja) 2020-01-23
FR3060753B1 (fr) 2019-07-26

Similar Documents

Publication Publication Date Title
EP2872884B1 (fr) Procédé de traitement de signaux issus d&#39;une acquisition par sondage ultrasonore, programme d&#39;ordinateur et dispositif de sondage à ultrasons correspondants
EP2294452B1 (fr) Procede et dispositif de sondage par propagation d&#39;ondes
EP3824280B1 (fr) Procédés et systèmes de caractérisation ultrasonore non invasive d&#39;un milieu hétérogène
WO1995026160A1 (fr) Procede et dispositif d&#39;evaluation et de caracterisation des proprietes des os
FR2971342A1 (fr) Dispositif d&#39;imagerie avec optimisation de cadence
FR3114157A1 (fr) Procédé et système de caractérisation ultrasonore d’un milieu
FR3114155A1 (fr) Procédé et système de caractérisation ultrasonore d’un milieu
FR3114158A1 (fr) Procédé et système de caractérisation ultrasonore d’un milieu
EP3555659B1 (fr) Procede d&#39;acquisition de signaux par sondage ultrasonore, programme d&#39;ordinateur et dispositif de sondage a ultrasons correspondants
EP3084416B1 (fr) Procédé de traitement de signaux issus d&#39;une acquisition par sondage ultrasonore, programme d&#39;ordinateur et dispositif de sondage à ultrasons correspondants
EP2342582B1 (fr) Procede et dispositif de sondage par propagation d&#39;ondes
EP3899523B1 (fr) Procédé d&#39;imagerie ultrasonore par transformée de fourier bidimensionnelle, programme d&#39;ordinateur et dispositif de sondage à ultrasons correspondants
Laroche et al. An inverse approach for ultrasonic imaging by total focusing point for close reflectors separation
FR3060753B1 (fr) Procede de traitement de signaux issus d&#39;une acquisition par sondage ultrasonore, programme d&#39;ordinateur et dispositif de sondage a ultrasons correspondants
FR3114159A1 (fr) Procédé et système de caractérisation ultrasonore d’un milieu
FR3114156A1 (fr) Procédé et système de caractérisation ultrasonore d’un milieu
EP4007914A1 (fr) Procédé et système de caractérisation ultrasonore non invasive d&#39;un milieu hétérogène
FR3140439A1 (fr) Procédé d’imagerie ultrasonore par transformée de Fourier multidimensionnelle à l’aide de deux transducteurs multiéléments distincts

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17821680

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 3046106

Country of ref document: CA

ENP Entry into the national phase

Ref document number: 2019532089

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2017821680

Country of ref document: EP

Effective date: 20190715