EP3555611A1 - Method of processing signals arising from an acquisition by ultrasound probing, corresponding computer program and ultrasound-based probing device - Google Patents

Method of processing signals arising from an acquisition by ultrasound probing, corresponding computer program and ultrasound-based probing device

Info

Publication number
EP3555611A1
EP3555611A1 EP17821680.0A EP17821680A EP3555611A1 EP 3555611 A1 EP3555611 A1 EP 3555611A1 EP 17821680 A EP17821680 A EP 17821680A EP 3555611 A1 EP3555611 A1 EP 3555611A1
Authority
EP
European Patent Office
Prior art keywords
matrix
singular
signals
emission
ultrasonic
Prior art date
Legal status (The legal status 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 status listed.)
Withdrawn
Application number
EP17821680.0A
Other languages
German (de)
French (fr)
Inventor
Sébastien ROBERT
Eduardo-Rigoberto LOPEZ VILLAVERDE
Claire PRADA JULIA
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Centre National de la Recherche Scientifique CNRS
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Original Assignee
Centre National de la Recherche Scientifique CNRS
Commissariat a lEnergie Atomique CEA
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
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 Centre National de la Recherche Scientifique CNRS, Commissariat a lEnergie Atomique CEA, Commissariat a lEnergie Atomique et aux Energies Alternatives CEA filed Critical Centre National de la Recherche Scientifique CNRS
Publication of EP3555611A1 publication Critical patent/EP3555611A1/en
Withdrawn legal-status Critical Current

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
    • 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
    • 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/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)
  • Chemical & Material Sciences (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Signal Processing (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Mathematical Physics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

This method (700) of processing ultrasound signals comprises the control (702) of emission transducers for M successive emissions of plane ultrasound waves exhibiting M different angles of emission, the control (702) of N reception transducers to simultaneously receive N temporal measurement signals per emission, and the obtaining (704) of a matrix [MR(t)] of ultrasound temporal signals, each coefficient MRi,j(t) of this matrix representing the measurement signal received by the i-th reception transducer due to the j-th emission. It furthermore comprises a decomposition into singular values (710) of a matrix [FTMR(f)] of frequency signals which is obtained (708) by transforming the matrix [MR(t)], an elimination (712) of some of the singular values and a reconstruction (714) of a denoised matrix [MRU(t)] of temporal signals on the basis of the non-eliminated singular values.

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.  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.
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.  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.
Elle porte plus particulièrement sur un procédé de traitement acquérant les signaux ultrasonores de la façon suivante :  It relates more particularly to a processing method acquiring the ultrasonic signals as follows:
- 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,  control of the emission transducers for M successive emissions of plane ultrasonic waves having M different successive emission angles in M emission zones,
- 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,  control of 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,
- 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.  obtaining a matrix [MR (t)] of ultrasound time signals of size NxM, each MR coefficient (t) of this matrix representing the measurement signal received by the ith receiving transducer due to the ith emission .
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.  Such an acquisition is generally done using a 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.
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. The matrix [MR (t)] of temporal signals obtained by this type of acquisition, generally described as a matrix of plane waves, 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. Given the current computing capabilities of the processors, this processing can be embedded in control instruments for real-time processing.
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.  This type of acquisition, generally referred to as "plane-wave compounding" or "plane-wave imaging", is for example described in the article by 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 image is of high quality and requires only a few tens of ultrasonic shots (usually M is between 10 and 30) for a sensor of N = 128 transmit / receive transducers, against 128 shots or almost for others imaging techniques such as the synthetic focus in all points based on an acquisition of type FMC (English "Full Matrix Capture") or equivalent. In addition, 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.
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. 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. However, whatever the chosen imaging technique, the images obtained in non-destructive testing may have a significant noise level depending on the properties of the medium being probed. For example, 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.
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.  It may thus be desired to design a method for processing ultrasonic signals that makes it possible to overcome at least some of the aforementioned problems and constraints while exploiting the advantageous acquisition of signals by successive transmissions of plane waves.
Il est donc proposé un procédé de traitement de signaux issus d'une acquisition par sondage ultrasonore comportant les étapes suivantes :  It is therefore proposed a signal processing method resulting from an ultrasonic sound acquisition comprising the following steps:
- 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,  control of the emission transducers for M successive emissions of plane ultrasonic waves having M different successive emission angles in M emission zones,
- 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  control of 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
- 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 obtaining a matrix [MR (t)] of ultrasonic time signals of size NxM, each coefficient MRj j (t) of this matrix representing the measurement signal received by the ith receiving transducer due to the jth emission and then
- 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,  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,
- élimination d'une partie des valeurs singulières et vecteurs singuliers associés issus de ladite décomposition en valeurs singulières, et elimination of a portion of the singular values and associated singular vectors resulting from said singular value decomposition, and
- reconstitution d'une matrice débruitée [MRu(t)] de signaux temporels à partir des valeurs singulières et vecteurs singuliers non éliminés.- reconstitution of a denoised matrix [MRu (t)] of time signals from singular values and singular vectors not eliminated.
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. Although 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.
De façon optionnelle :  Optionally:
- 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 the elimination of a part of the singular values and associated singular vectors is done according to a predetermined criterion of distinction between singular values related to defects and singular values related to noise, and
- 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. 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.
De façon optionnelle également, la transformation et la transformation inverse sont des transformations de Fourier discrètes.  Also optionally, the transformation and the inverse transformation are discrete Fourier transformations.
De façon optionnelle également, l'élimination d'une partie des valeurs singulières et vecteurs singuliers associés se fait par :  Optionally also, the elimination of a part of singular values and associated singular vectors is done by:
- 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  comparison of a decrease in amplitudes of the singular values with a theoretical curve of decay of singular values resulting from a theoretical noise matrix whose components are independent Gaussian random variables, then
- suppression des valeurs singulières appartenant, à un coefficient de proportionnalité commun près, à la courbe théorique de décroissance. - removal of the singular values belonging to a coefficient of proportionality common to the theoretical decay curve.
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 : Optionally also, 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:
où Π[0;2] désigne la fonction porte sur l'intervalle [0;2]. where Π [ 0; 2 ] designates the function relates to the interval [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é. Optionally also, an ultrasonic signal processing method according to the invention 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.
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.  Optionally also, 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.
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.  Optionally also, 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.
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.  There is also provided 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.
II est également proposé un dispositif de sondage à ultrasons comportant : It is also proposed an ultrasonic sounding device comprising:
- une sonde comprenant L transducteurs d'émission à ultrasons et N transducteurs de réception à ultrasons, a probe comprising L ultrasonic emission transducers and N ultrasonic reception transducers,
- 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,  control means for the L emission transducers for M successive emissions of plane ultrasonic waves having M different successive emission angles in M emission zones,
- 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  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
- 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,  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,
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, 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,
- élimination d'une partie des valeurs singulières et vecteurs singuliers associés issus de ladite décomposition en valeurs singulières, et elimination of a portion of the singular values and associated singular vectors resulting from said singular value decomposition, and
- reconstitution d'une matrice débruitée [MRu(t)] de signaux temporels à partir des valeurs singulières et vecteurs singuliers non éliminés.- reconstitution of a denoised matrix [MRu (t)] of time signals from singular values and singular vectors not eliminated.
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 : The invention will be better understood with the aid of the description which follows, given solely by way of example and with reference to the appended drawings in which:
- 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, FIG. 1 schematically represents the general structure of an ultrasonic sounding device according to one embodiment of the invention,
- la figure 2 illustre un principe d'émissions successives d'ondes ultrasonores planes mis en œuvre par le dispositif de la figure 1 , FIG. 2 illustrates a principle of successive transmissions of plane ultrasonic waves implemented by the device of FIG. 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,  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,
- 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,  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,
- la figure 6 illustre trois exemples d'images reconstituées sans ou avec application des principes de l'invention, et  FIG. 6 illustrates three examples of images reconstituted without or with application of the principles of the invention, and
- 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.  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.
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. With reference to FIG. 1, a sounding device 100 of an object 102 according to one embodiment of the invention 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.
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. 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. In the mode In the embodiment of FIG. 1, 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. But in another equivalent embodiment, the probe 104 could be in direct contact with the object 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. 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.
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. 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. Thus, in the non-limiting example of FIG. 1, 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. In addition, the number L of transmitters could quite be different from the number N of receivers.
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. 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.
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 à : 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:
- 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. activate the L = N transducers 108i, ..., 108 N as transmitters for M successive emissions of plane ultrasonic waves having M different successive emission angles in M emission regions of the object 102, activating the transducers 108i,..., 108 N as receivers for, following each transmission, simultaneously receiving, by these N receivers and for a predetermined duration of the desired inspection depth, N measurement time signals, measuring in particular echoes due to reflections of each emission considered.
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. 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. As a result, there are as many delay laws as there are successive emissions desired.
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. As illustrated in FIG. 2 in a case where the number M of successive emissions is odd and the emission angles succeed each other with a constant pitch in a symmetrical angular sector with respect to the direction z orthogonal to the transducer network I O8 1 , ..., 108 N , 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. Finally, the last plane wave emission is associated with a delay law T M allowing the emission of a plane wave of emission angle Θ Μ = - with respect to the direction z in a last emission zone ZE M partially situated outside the opening of the probe 104. generally, the jth plane is associated with a T delay law wave transmitting for emitting a wave angle planar emission 0 j = θ + (d - 1). (θ Μ - Θ ^ Μ - 1) with respect to the z direction. In reality, most often M is even: there is no emission at 0 ° and the angular pitch is not 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. Given the acquisition technique used, 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. In particular, 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.
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. As a variant, and as advantageously made possible by the principle of acquiring signals by successive transmissions of plane waves, 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. In particular, when it is a question of detecting defects such as a crack disposed at the bottom of an object to be inspected in non-destructive inspection, 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.
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. In order to improve the quality of the measurement signals used to reconstitute the imaged zone, it is also possible to apply 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.
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. With reference again to FIG. 1, on reception of the signals resulting from each of the M successive transmissions, 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.
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. 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.
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.  Optionally, 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.
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. 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.
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. 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. Although it is known to break down into singular values a 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)]. Indeed, the matrix denoted [PFTMR (f)] of MxM and dimensions defined by the product [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 physical interpretation of singular values and singular vectors is then not the same as in the method of decomposition of the Temporary Turnaround Operator (DORT method) as taught in the article by Prada et al entitled "Eigenmodes of the time reversal operator: a solution to selective multiple-target media ", published in Wave Motion 20, pages 151 -163 (1994). In addition, the matrix [FTMR (f)] is generally neither square nor symmetrical, M being often even significantly less than N in acquisition by emission of plane waves. In this way, a decomposition into singular values does not produce diagonalization as in most works dealing with the matrix of inter-element impulse responses.
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 :  More precisely, the operation of decomposition into singular values makes it possible to estimate the matrices U (f), S (f) and V (f) such that:
[FTMR(f)] = u{f)s{f)v {f) =∑ai {f).ul {f).v {f) , où U(f) = [ui (f) , .... uN(f)] et i=\ [FTMR (f)] = u {f) s {f) v {f) = Σa i {f) u (f) .v {f), where U (f) = [ui (f), .... u N (f)] and 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)] : V = [v ! (f), v M (f)] are orthogonal unitary matrices of respective sizes NxN and MxM, which respectively contain the singular vectors in reception and transmission, where S is a real matrix of size NxM containing K = min (N; M) singular values a, (f) of the matrix [FTMR (f)], ordered decreasingly at a given frequency f σι (ί) ≥ ... ≥ a K (f) ≥ 0, and where min (N) M) is the function that returns the minimum value between N and M (usually M). 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). This elimination is done according to a criterion of distinction between singular values related to defects and singular values related to noise, the first being amplitudes higher than the seconds. Given that a ^ f) ≥ ... ≥ a K (f) ≥ 0, we have to find the integer value P (f) between 1 and K such that a ^ f), a P ( f) (f) can be considered as related to defects to be detected in the object 102 and P (f) +1 (f), a K (f) can be eliminated because they are considered to be related to noise. In the case of small defects ideally spaced from each other, P (f) is equal to the number of defects present in the area of interest inspected. One can then write the matrix [FTMR (f)] in the form of an addition of two matrices [FTMRu (f)] and [FTMR N (f)]:
[FTMR(f)] = [FTMRU (/)] + [FTMRN (/)] , avec [FTMR (f)] = [FTMR U (/)] + [FTMR N (/)], with
Pif)  Pif)
[FTMRU (/)] = ^σ,. (/).«,. (/).v,†(/) , appelée matrice de signal utile, et [FTMR U (/)] = ^ σ ,. (/). ",. (/).v, † (/), called a useful signal matrix, and
i=l [FTMRN ( f)] = ∑ ; ( ( ).ν; ( ) , appelée matrice de bruit. i = 1 [FTMR N (f)] = Σ ; (() .ν; (), called the noise matrix.
;=p(/)+i  ; = P (/) + i
Réduire le rang de la matrice [FTMR(f)] revient ainsi à ne conserver que [FTMRu(f)].  Reducing the rank of the matrix [FTMR (f)] thus amounts to retaining only [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 :  In practice, the determination of P (f) is done by a study of the decay curve of the amplitudes of the singular values. This study can be done on the basis of random matrix theory as taught in:
- 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  Marcenko et al's article, "Distribution of eigenvalues for some sets of random matrices", published in Mathematics of the USSR-Sbornik, vol. 1, No. 4, pp. 457-483 (1967), and
- 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] , the article by Aubry et al entitled "Detection and imaging in a random medium: a matrix method to overcome multiple scattering and aberration", published in Journal of Applied Physics, 106 (4), 044903 (2009). Assuming that all the components of the noise matrix [FTMR N (f)] are independent Gaussian random variables, it can indeed be shown that the singular values related to the noise σ Ρ (ί) + ι (f), a K (f) belong, with a coefficient of proportionality common to a theoretical curve defined by a reciprocal function F "1 (1 -a) itself defined by a function F (a), called distribution of random singular values , with: Fia) = i + 2 arcsin ()] Π [0 . 2] ,
où Π[0;2] désigne la fonction porte sur l'intervalle [0;2]. where Π [ 0; 2 ] designates the function relates to the interval [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. 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.
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.  Thus, the number of singular values calculated by executing instructions 1 30 that do not belong to this theoretical curve defines the rank P (f) of the useful signal matrix [FTMRu (f)]. By way of example, FIG. 3 illustrates in arbitrary unit a.u. ("arbitrary unit") a singular value distribution calculated at a frequency f = 5 MHz, for an acquisition made with a polyethylene sample by emitting M = 64 plane waves between -31.5 ° and +31, 5 °. Two singular values stand out from the theoretical curve, of which one quite clearly, indicating in this example that P (f) = 2 to 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.  To calculate the matrix [FTMRu (f)] of useful signal on several frequencies f, it is sufficient to identify at each desired frequency the number of singular values not belonging to the theoretical curve.
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é. 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.
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.  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.
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. Finally, 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. Alternatively or in addition, the processing performed by instructions 136 could include adaptive and selective focusing.
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é : According to a relatively simple but not always optimal adaptation of the synthetic focus in all points, in each pixel P of the reconstructed digital image, 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:
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. where μ (Ρ) 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.
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 : This method of calculation is often not optimal because, as illustrated in FIG. 4, all the plane waves emitted by the probe 104 do not interact with the D or defects related to the singular values of the useful signal matrix. [FTMRu (f)], that is the singular values a ^ f), a P (f) (f). It may then be advantageous to select a restricted angular sector [9 m1 ; 9 m 2 ] in the set of successive emissions outside which the plane waves do not interact with the above-mentioned defect (s). This reduced angular sector is defined by its minimum (m) and maximum (m 2 ) emission indices where ~ \ ≤m † ≤m 2 ≤M. The module A (P) is then preferably calculated by the following coherent addition:
) μ (Ρ)ρί(Ρ) «[/ [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) . ) μ (Ρ) ρ ί (Ρ) " [/ [t y (P)]. In the case of a defect similar to a point reflector, a clever method for determining the values of ΙΤΗ and m 2 will now be detailed with reference to FIG. 5. It consists in comparing, for at least one of the singular values σι (ί), a P (f) (f), and in particular at least for the first singular value a ^ f) which is the highest, an experimental phase value of the singular vector which is associated with it in emission with a value theoretical phase in the presence of the defect D. 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. In FIG. 5, 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 ).
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 :  Under the experimental conditions given above, the theoretical phase value is a corrected phase calculated as follows to always be negative:
2nfc(XD sin 9j + ZD cos 2nf c (X D sin 9j + Z D cos
avec c la célérité de l'onde plane dans le milieu considéré et 6j son angle d'émission. with c the speed of the plane wave in the medium considered and 6j its emission angle.
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) . 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).
La définition ci-dessus de la phase théorique corrigée donne la courbe en trait continu de la figure 5.  The above definition of the corrected theoretical phase gives the continuous line curve of 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 :  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, <pj1( c) = arg[i 1( c)] . - max [arg[¾ e)] .] . Vy, 1 <y M M, <ρ 1 (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.  This definition of the corrected phase of the singular vector v ^ y gives the M values represented by small circles in 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. It can be seen that the experimental and theoretical values are in good correspondence within a restricted angular sector whose extremities give the values of ΙΤΗ and m 2 . Outside this narrow angular sector, the experimental values are very different from the theoretical values.
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 . In a general way, if the comparison is carried out for several singular values linked to the defect (s) D, 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.
La figure 6 illustre, dans un exemple d'objet sondé comportant un défaut D circulaire central, trois images obtenues :  FIG. 6 illustrates, in an example of a probed object comprising a central circular defect D, three images obtained:
- par une acquisition FMC avec reconstitution de l'image par focalisation synthétique en tous points (image de gauche),  by a FMC acquisition with reconstruction of the image by synthetic focusing in all points (left image),
- par une acquisition et un traitement tels qu'enseignés dans le document WO 2015/092250 A1 (image au centre), et  by acquisition and treatment as taught in WO 2015/092250 A1 (image in the middle), and
- 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).  by acquisition and filtering treatment and selection of a restricted angular sector as taught according to the principles of the present invention (right image).
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.  There is a strong attenuation of the noise (gain of 20 dB in particular between the image in the center and the right one in terms of SNR) and a better visibility of the defect.
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.  With reference to FIG. 7, 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.
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)]. During a step 702, 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)].
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. These sequences are M in number, integer can be much lower than the number N of transducers 108i, ..., 108 N. After each shot, the signals are received on all N transducers, digitized and transmitted to the electronic circuit 1 12.
Au cours d'une étape 704, l'unité de traitement 1 14 exécutant les instructions During a step 704, the processing unit 1 14 executing the 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. 124 records the measurement signals MRj j (t), these signals being digitized to facilitate their subsequent processing. 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.
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. During an optional step 706, 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.
Au cours d'une étape 708, l'unité de traitement 1 14 exécutant les instructions During a step 708, the processing unit 1 14 executing the 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. 128 performs a discrete Fourier transform of the matrix [MR (t)] to obtain the matrix [FTMR (f)] of frequency signals.
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.  During a step 710, the processing unit 1 14 executing the instructions 130 performs a singular value decomposition of the matrix [FTMR (f)], as detailed previously.
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)].  During a step 712, 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)].
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.  During a step 714, 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.
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. During a step 716, 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.
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é.  Finally, during a last step 718, 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.
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é. It clearly appears that an ultrasonic sounding device such as that described above makes it possible to obtain a de-fuzzy and very good quality fault detection (s), better than those obtained by the FMC acquisitions with synthetic focusing image reconstructions. in all respects or by the treatment as taught in document WO 2015/092250 A1. In addition, the teaching of the present patent application can be advantageously combined with that of WO 2015/092250 A1. Note also that the invention is not limited to the embodiment described above. It will be apparent to those skilled in the art that various modifications can be made to the embodiment described above, in the light of the teaching that has just been disclosed.
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.  In particular, the computer program instructions could be replaced by electronic circuits dedicated to the functions performed during the execution of these 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é.  In general terms, in the following claims, the terms used are not to be construed as limiting the claims to the embodiment set forth in this specification, but should be interpreted to include all equivalents that the claims are intended to cover. because of their formulation and whose prediction is within the reach of those skilled in the art by applying his general knowledge to the implementation of the teaching that has just been disclosed.

Claims

REVENDICATIONS
1 . Procédé (700) de traitement de signaux (S) issus d'une acquisition par sondage ultrasonore comportant les étapes suivantes : 1. Method (700) for processing signals (S) resulting from an ultrasonic sound acquisition comprising the following steps:
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), control (702) of the emission transducers (108 1; 108 N ) for M successive emissions of plane ultrasonic waves having M different successive emission angles (θ 1; ..., Θ Μ ) in M zones of emission (ZE 1, ZE M ),
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, control (702) of N receiving transducers (108 1; 108 N ) so as to receive (704) simultaneously and for a predetermined duration, for each transmission, N time measurement signals, in particular measuring echoes due to reflections of the program in question,
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, obtaining (704) a matrix [MR (t)] of ultrasonic time signals of size NxM, each coefficient MRi j (t) of this matrix representing the measurement signal received by the ith receiving transducer due to the -th issue,
caractérisé en ce qu'il comporte en outre les étapes suivantes : characterized in that it further comprises the following steps:
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,  transformation (708) of the matrix [MR (t)] of time signals into a matrix [FTMR (f)] of frequency signals, then decomposition into singular values (710) of the matrix [FTMR (f)] of frequency signals,
é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  elimination (712) of a part of the singular values and associated singular vectors resulting from said singular value decomposition (710), and
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.  reconstitution (712, 714) of a denoised matrix [MRu (t)] of time signals from the singular values and singular vectors not eliminated.
2. Procédé (700) de traitement de signaux ultrasonores (S) selon la revendication 1 , dans lequel :  The method (700) of ultrasonic signal processing (S) according to claim 1, wherein:
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. the elimination (712) of a part of the singular values and associated singular vectors is done according to a predetermined criterion of distinction between singular values related to defects and singular values related to noise, and the reconstitution (712, 714) of the denoised matrix [MRu (t)] of time signals comprises a reconstruction (712) of a denoised matrix [FTMRu (f)] of frequency signals from singular values and singular vectors not eliminated , then an inverse transformation (714) of this denoised matrix [FTMRu (f)] of frequency signals into the denoised matrix [MRu (t)] of time signals.
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.  An ultrasonic signal processing method (700) according to claim 2, wherein the transformation (708) and the inverse transform (714) are discrete Fourier transformations.
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 :  4. Method (700) of ultrasonic signal processing (S) according to claim 2 or 3, wherein the elimination (712) of a portion of singular values and associated singular vectors is by:
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  comparison of a decrease of amplitudes of the singular values with a theoretical curve of decay of singular values resulting from a theoretical noise matrix whose components are independent Gaussian random variables, then
suppression des valeurs singulières appartenant, à un coefficient de proportionnalité commun près, à la courbe théorique de décroissance.  suppression of the singular values belonging to a coefficient of proportionality common to the theoretical decay curve.
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 : 5. Method (700) for processing ultrasonic signals (S) according to claim 4, in which the theoretical decay curve is defined by a reciprocal function F "1 (1-a) itself defined by a function F (a ), called distribution of random singular values, such that:
où Π[0;2] désigne la fonction porte sur l'intervalle [0;2]. where Π [ 0; 2 ] designates the function relates to the interval [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é. The method (700) for processing ultrasonic signals (S) according to any one of claims 1 to 5, further comprising a reconstruction (716, 718) of an image area calculated by calculation, at each point of a plurality of predetermined points of this imaged zone, of a value resulting from a coherent summation of instantaneous values taken respectively by at least a part of the NxM time signals of the matrix [MRu (t)] at flight times respectively corresponding to a passing through the point considered according to a predetermined propagation mode.
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. 7. The method (700) of ultrasonic signal processing (S) according to claim 6, wherein the calculation is done on part of the NxM time signals of the matrix [MRu (t)] in a restricted angular sector (m 1; m 2 ) in all successive issues.
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. The method (700) of ultrasonic signal processing (S) according to claim 7, wherein the restricted angular sector (m, m 2 ) is selected (716) based on a comparison, for at least one singular values not eliminated, an experimental phase value of the singular vector which is associated in emission with a theoretical phase value in the presence of said at least one defect related to this singular vector, so that the plane waves which in are excluded do not interact with at least one defect related to singular values and singular vectors not eliminated.
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).  9. Computer program (1 18) downloadable from a communication network and / or recorded on a computer readable medium and / or executable by a processor (1 14), characterized in that it comprises instructions (120, 124, 126, 128, 130, 132, 134, 136) for carrying out the steps of an ultrasonic signal processing method (700) (S) according to any one of claims 1 to 8, when said program is executed on a computer (1 12).
10. Dispositif (100) de sondage à ultrasons comportant :  An ultrasonic sounding device (100) comprising:
une sonde (104) comprenant L transducteurs d'émission à ultrasons (108i , 108N) et N transducteurs de réception à ultrasons (108i , 108N), a probe (104) comprising L ultrasonic emission transducers (108i, 108 N ) and N ultrasonic receiving transducers (108i, 108 N ),
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), means (1 12) for controlling the L emission transducers for M successive emissions of plane ultrasonic waves having M different successive emission angles (θ 1; ..., Θ Μ ) in M emission zones (ZE 1; ZE M ),
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  means (1 12) for controlling the N receiving transducers so as to receive, simultaneously and for a predetermined duration, for each transmission, N measurement time signals, in particular measuring echoes due to reflections of the emission in question, and
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 : a processor (1 14) for reconstituting a matrix [MR (t)] of ultrasonic time signals of size NxM, each coefficient MRi j (t) of this matrix representing the measurement signal received by the ith receiving transducer due to the j-th emission, characterized in that the processor (1 14) is further configured to perform the following processes:
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.  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 elimination of a part of the values singular and associated singular vectors resulting from said singular value decomposition, and reconstitution of a denoised matrix [MRu (t)] of temporal signals from singular values and singular vectors not eliminated.
EP17821680.0A 2016-12-15 2017-12-05 Method of processing signals arising from an acquisition by ultrasound probing, corresponding computer program and ultrasound-based probing device Withdrawn EP3555611A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1662525A FR3060753B1 (en) 2016-12-15 2016-12-15 METHOD FOR PROCESSING ULTRASONIC SURVEY ACQUISITION SIGNALS, COMPUTER PROGRAM, AND CORRESPONDING ULTRASONIC SURVEY DEVICE
PCT/FR2017/053391 WO2018109314A1 (en) 2016-12-15 2017-12-05 Method of processing signals arising from an acquisition by ultrasound probing, corresponding computer program and ultrasound-based probing device

Publications (1)

Publication Number Publication Date
EP3555611A1 true EP3555611A1 (en) 2019-10-23

Family

ID=58609497

Family Applications (1)

Application Number Title Priority Date Filing Date
EP17821680.0A Withdrawn EP3555611A1 (en) 2016-12-15 2017-12-05 Method of processing signals arising from an acquisition by ultrasound probing, corresponding computer program and ultrasound-based probing device

Country Status (6)

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

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2946753B1 (en) 2009-06-11 2011-07-22 Centre Nat Rech Scient ULTRASONIC METHOD AND DEVICE FOR CHARACTERIZING A MEDIUM
FR2993362B1 (en) 2012-07-12 2016-07-01 Commissariat Energie Atomique METHOD FOR PROCESSING ULTRASONIC SURVEY ACQUISITION SIGNALS, COMPUTER PROGRAM, AND CORRESPONDING ULTRASONIC SURVEY DEVICE
FR3015742B1 (en) 2013-12-20 2016-01-22 Commissariat Energie Atomique METHOD FOR PROCESSING ULTRASONIC SURVEY ACQUISITION SIGNALS, COMPUTER PROGRAM, AND CORRESPONDING ULTRASONIC SURVEY DEVICE

Also Published As

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

Similar Documents

Publication Publication Date Title
EP2872884B1 (en) Method of processing ultrasonic testing data, computer program and corresponding ultrasonic testing apparatus
EP2294452B1 (en) Method and device for measuring using wave propagation
EP3824280B1 (en) Methods and systems for non-invasively characterising a heterogeneous medium using ultrasound
EP0751743A1 (en) Method and device for estimating and characterising bone properties
FR2971342A1 (en) IMAGING DEVICE WITH CADENCE OPTIMIZATION
EP3084416B1 (en) Method for processing signals from an ultrasound probe acquisition, corresponding computer program and ultrasound probe device
FR3114157A1 (en) METHOD AND SYSTEM FOR ULTRASOUND CHARACTERIZATION OF A MEDIUM
EP3899523B1 (en) Method for ultrasound imaging using two-dimensional fourier transform, corresponding ultrasound probe device and computer program
FR3114155A1 (en) METHOD AND SYSTEM FOR ULTRASOUND CHARACTERIZATION OF A MEDIUM
FR3114158A1 (en) METHOD AND SYSTEM FOR ULTRASOUND CHARACTERIZATION OF A MEDIUM
EP3555659B1 (en) Method of acquiring ultrasonic testing signals, and corresponding computer program and ultrasonic testing device
Laroche et al. An inverse approach for ultrasonic imaging by total focusing point for close reflectors separation
FR3060753B1 (en) METHOD FOR PROCESSING ULTRASONIC SURVEY ACQUISITION SIGNALS, COMPUTER PROGRAM, AND CORRESPONDING ULTRASONIC SURVEY DEVICE
FR3114159A1 (en) METHOD AND SYSTEM FOR ULTRASOUND CHARACTERIZATION OF A MEDIUM
FR3114156A1 (en) METHOD AND SYSTEM FOR ULTRASOUND CHARACTERIZATION OF A MEDIUM
WO2021023933A1 (en) Method and system for non-invasively characterising a heterogeneous medium using ultrasound
FR3140439A1 (en) Ultrasound imaging method by multidimensional Fourier transform using two separate multi-element transducers

Legal Events

Date Code Title Description
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: UNKNOWN

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20190603

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

DAV Request for validation of the european patent (deleted)
DAX Request for extension of the european patent (deleted)
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20220701