FR3099580A1 - Procédé et système de caractérisation ultrasonore non invasive d’un milieu hétérogène - Google Patents
Procédé et système de caractérisation ultrasonore non invasive d’un milieu hétérogène Download PDFInfo
- Publication number
- FR3099580A1 FR3099580A1 FR1908904A FR1908904A FR3099580A1 FR 3099580 A1 FR3099580 A1 FR 3099580A1 FR 1908904 A FR1908904 A FR 1908904A FR 1908904 A FR1908904 A FR 1908904A FR 3099580 A1 FR3099580 A1 FR 3099580A1
- Authority
- FR
- France
- Prior art keywords
- point
- calculated
- output
- input
- medium
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 77
- 238000012512 characterization method Methods 0.000 title claims abstract description 29
- 238000002604 ultrasonography Methods 0.000 title claims description 42
- 230000008569 process Effects 0.000 title claims description 11
- 230000004044 response Effects 0.000 claims abstract description 83
- 239000011159 matrix material Substances 0.000 claims abstract description 56
- 238000004364 calculation method Methods 0.000 claims description 53
- 238000005259 measurement Methods 0.000 claims description 37
- 238000005457 optimization Methods 0.000 claims description 10
- 239000013598 vector Substances 0.000 claims description 10
- 230000010349 pulsation Effects 0.000 claims description 9
- 238000004088 simulation Methods 0.000 claims description 6
- 101710084218 Master replication protein Proteins 0.000 claims description 3
- 101710112083 Para-Rep C1 Proteins 0.000 claims description 3
- 101710112078 Para-Rep C2 Proteins 0.000 claims description 3
- 102100022881 Rab proteins geranylgeranyltransferase component A 1 Human genes 0.000 claims description 3
- 102100022880 Rab proteins geranylgeranyltransferase component A 2 Human genes 0.000 claims description 3
- 101710119887 Trans-acting factor B Proteins 0.000 claims description 3
- 101710119961 Trans-acting factor C Proteins 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims 1
- 238000012360 testing method Methods 0.000 description 11
- 230000005540 biological transmission Effects 0.000 description 10
- 238000009792 diffusion process Methods 0.000 description 10
- 238000003384 imaging method Methods 0.000 description 8
- 230000004075 alteration Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 239000000523 sample Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 5
- 210000004185 liver Anatomy 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000001934 delay Effects 0.000 description 4
- 235000013372 meat Nutrition 0.000 description 4
- 238000012285 ultrasound imaging Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 3
- 238000006731 degradation reaction Methods 0.000 description 3
- 238000001727 in vivo Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 210000000577 adipose tissue Anatomy 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000002091 elastography Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000012623 in vivo measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000003387 muscular Effects 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
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/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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
- A61B8/4488—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer the transducer being a phased array
-
- 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
-
- 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
-
- 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/4463—Signal correction, e.g. distance amplitude correction [DAC], distance gain size [DGS], noise filtering
-
- 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
-
- 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
- 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/52046—Techniques for image enhancement involving transmitter or receiver
- G01S7/52049—Techniques for image enhancement involving transmitter or receiver using correction of medium-induced phase aberration
-
- 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/02—Indexing codes associated with the analysed material
- G01N2291/024—Mixtures
- G01N2291/02483—Other human or animal parts, e.g. bones
-
- 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/02—Indexing codes associated with the analysed material
- G01N2291/024—Mixtures
- G01N2291/02491—Materials with nonlinear acoustic properties
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Gynecology & Obstetrics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
Procédé et système de caractérisation ultrasonore non invasive d’un milieu hétérogène Procédé de caractérisation ultrasonore non invasive d’un milieu hétérogène comprenant une étape de génération d’une série d’ondes ultrasonores incidentes, une étape d’enregistrement d’une matrice de réflexion expérimentale R ui (t) définie entre la base d’émission (i) en entrée et une base de réception (u) en sortie, une étape de détermination d’une réponse REP(r, Δ r) du milieu entre un transducteur virtuel d’entrée (TVin) calculé à partir d’une focalisation en entrée de la matrice de réflexion expérimentale créant une tache focale d’entrée autour d’un premier point (P1) et un transducteur virtuel de sortie (TVout) calculé à partir d’une focalisation en sortie de la matrice de réflexion expérimentale créant une tache focale de sortie autour d’un deuxième point (P2), ladite réponse étant exprimée en fonction d’un point central (PC) de position spatiale (r) dans le milieu situé au milieu des premier et deuxième points (P1, P2). Figure pour l’abrégé : FIGURE 4
Description
La présente description concerne des procédés et systèmes de caractérisation ultrasonore non invasive d’un milieu hétérogène, et s’applique notamment à l’imagerie médicale ou au contrôle non destructif et plus généralement à tous les domaines dans lesquels l’imagerie ultrasonore peut être utilisée.
ETAT DE L’ART
Dans le domaine de l’imagerie acoustique, on cherche à caractériser un milieu totalement ou partiellement inconnu en le sondant de manière active au moyen d’ondes ultrasonores. C’est notamment le principe de l’échographe utilisé en imagerie médicale.
La résolution d'un système d'imagerie acoustique peut être définie comme la capacité à discerner les petits détails d'un objet. En principe, un système d'imagerie acoustique est limité par la diffraction et la résolution théorique est donnée par λ / 2 (où λ est la longueur d'onde du son dans le milieu), ou par l'ouverture angulaire finie du détecteur. En pratique, la résolution est cependant souvent détériorée par les variations de vitesses du son dans le milieu de propagation.
En effet, la plupart du temps en imagerie acoustique, le milieu est considéré comme homogène, avec une vitesse du son c0 constante. Or l’hypothèse d’un milieu homogène n’est pas toujours respectée. Par exemple, dans le cas de l’échographie du foie, la sonde est placée entre les côtes du patient. Les ondes acoustiques traversent une succession de couches de graisse et de muscle avant d’atteindre l’organe visé. Les tissus mous possèdent chacun des propriétés mécaniques différentes. La vitesse du son est donc loin d’être homogène, entre 1450 m/s pour les tissus adipeux et 1600 m/s pour le foie. Les variations de vitesse du son entraînent un déphasage différent des ondes suivant les endroits par lesquels elles se propagent. Il en résulte une aberration du front d’onde acoustique qui conduit à une distorsion de l'image échographique résultante, et donc à une dégradation de sa résolution et de son contraste. Ces aberrations peuvent être telles qu’elles compromettent les résultats de l’examen médical.
Comme illustré sur lesfigures 1A à 1C, les méthodes d’échographie conventionnelle utilisent un réseau 10 de transducteurs piézo-électriques 11 qui peuvent émettre et recevoir des impulsions ultrasonores de manière indépendante. On repère la position de chacun des transducteurs par le vecteur u. Lorsqu’un tel réseau est placé en vis-à-vis d’un milieu 20 que l’on cherche à étudier, celui-ci peut être insonifié et imagé de différentes manières.
La première manière pour générer une image échographique du milieu à étudier est d’émettre une impulsion ultrasonore depuis un des transducteurs du réseau dont la position est repérée par le vecteur (figure 1A,schéma de gauche). Cela donne lieu à une onde incidente cylindrique (ou sphérique) divergente pour un réseau 1D (ou 2D) de transducteurs. Celle-ci est réfléchie par les diffuseurs 21 du milieu 20 et le champ rétrodiffusé est enregistré par chacun des transducteurs 11 en fonction du temps (figure 1A,schéma de droite). En répétant cette opération avec chaque transducteur utilisé successivement comme source, on mesure l’ensemble des réponses impulsionnelles entre chaque transducteur, où le vecteur désigne la position du détecteur. Ces réponses forment la matrice de réflexion exprimée dans la base des transducteurs. L’intérêt d’une telle mesure réside dans le fait que cette matrice contient toute l’information sur le milieu étudié, un ensemble d’opérations matricielles pouvant ensuite lui être appliquées à des fins d’imagerie du milieu, par exemple. En revanche, une telle acquisition suppose que le milieu reste fixe pendant toute la durée des mesures, ce qui peut être très difficile dans le cas d’utilisation in-vivo. De plus, l’énergie émise par un seul élément piézo-électrique est faible, ce qui peut induire un mauvais rapport signal à bruit.
D’autres méthodes sont connues pour générer une image du milieu à étudier dans lesquelles on effectue des émissions focalisées par une technique de formation de voies (ou « beamforming » en langue anglaise). Comme montré par lafigure 1B, schéma de gauche, ces méthodes consistent à appliquer aux transducteurs 11 un jeu de retards appropriés, basés sur un modèle de vitesse homogène, afin de corriger les temps de parcours des ondes pour que toutes les impulsions arrivent ensemble au point focal visé, de position . En raison des limites physiques de la diffraction, on concentre les ultrasons dans une zone délimitée par l’ouverture de la sonde échographique. Afin de construire une image échographique, une étape de focalisation est également effectuée en réception. L’ensemble des échos captés par les éléments 11 du réseau 10 est alors traité pour simuler l’effet d’une lentille en réception, comme décrit sur lafigure 1B,schéma de droite. Les signaux reçus par les transducteurs sont remis en phase en les décalant temporellement. Ces délais sont identiques à ceux appliqués à l’émission. Dans la phase d’émission, tous les signaux interfèrent au point de position . En réception, les signaux provenant de ce même point interfèrent électroniquement par sommation des signaux au temps balistique . Cette sommation donne le résultat final de la focalisation en réception. La méthode illustrée sur lafigure 1B, dite méthode confocale à double focalisation à l’émission et à la réception permet d’imager directement la réflectivité du milieu avec une résolution latérale limitée par la diffraction, une excellente résolution axiale seulement limitée par la durée de l’impulsion initiale et un excellent contraste. Toutefois, cette méthode est chronophage car elle nécessite de focaliser physiquement à l’émission en chacun des points du milieu ou au moins à une profondeur, sur chacune des lignes de l’image.
Une autre technique d’imagerie, développée plus récemment, consiste à générer une image du milieu en insonifiant le milieu avec une série d’ondes planes. Lafigure 1Cillustre le principe de cette échographie dite en ondes planes, décrite par exemple dans l’article de G. Montaldo et al.« Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography”(IEEE Trans. Ultrason., Ferroelect. Freq. Control 56 489-506, 2009). Des retards sont appliqués sur chaque signal à l’émission (figure 1C,schéma de gauche) pour la formation d’un front d’onde incliné d’un angle par rapport au réseau de transducteurs 10. A la réception (figure 1C,schéma de droite), le champ rétrodiffusé par le milieu, est mesuré par tous les capteurs de positions pour une série d’ondes planes incidentes dont on fait varier l’angle d’incidence . L’ensemble de ces réponses forment une matrice de réflexion définie entre la base de Fourier (ou base d’ondes planes) en entrée et la base des transducteurs en sortie. Une fois cette matrice enregistrée, les signaux sont décalés temporellement avant d’être sommés de manière cohérente afin de focaliser numériquement les données à l’émission et à la réception pour chaque point de position . Le nombre d’acquisitions nécessaires pour former une image échographique est réduit par rapport à l’échographie standard (émissions focalisées), et ceci pour un même niveau de contraste et de résolution de l’image échographique.
Lafigure 2illustre l’influence d’aberrations du milieu sur les méthodes d’imagerie échographique conventionnelle (figures 1A à 1C). Les retards initialement déterminés et à appliquer sur chacun des transducteurs du réseau à l’émission et à la réception ne sont pas optimaux pour l’évaluation d’une image du milieu, car ils sont déterminés avec l’hypothèse d’un milieu homogène avec une vitesse du son constante. Une couche aberratrice 22 induit une distorsion du front d’onde incident. A l’émission ou excitation, étape 25, les lois de retard utilisées ne permettent pas de concentrer l’énergie acoustique dans zone délimitée par les limites de la diffraction, zones usuellement appelée de tache focale. A la réception, en étape 26, les lois de retard utilisées ne permettent pas de sélectionner correctement les signaux ultrasonores provenant du point focal du milieu, et mélangent les signaux provenant d’une tache focale également aberrée. Il en résulte une double aberration dans le processus de construction d’image, ce qui dégrade fortement sa résolution. De nouvelles lois de retard peuvent alors être recalculées afin de compenser l’effet de la couche aberratrice en ajoutant une loi de retard supplémentaire aux délais généralement utilisés en formation de voies.
Cependant, ces corrections en aberration ne corrigent pas complètement ni ces aberrations ni la dégradation en résolution. Il existe un besoin de mieux estimer la qualité de la focalisation dans le milieu.
Le document“The van Cittert-Zernike theorem in pulse echo measurements”, (Raoul Mallart and Mathias Fink, J. Acoust. Soc. Am. 90(5), November 1991) a étudié les propriétés statistiques du champ réfléchi par un milieu aléatoire en régime de diffusion simple. Il est notamment montré que, pour une onde incidente focalisée, la covariance spatiale du champ réfléchi est proportionnelle, depuis le champ lointain, à la transformée de Fourier de la fonction d’ouverture en transmission. Autrement dit, ce théorème explique que l’étude des propriétés statistiques du champ réfléchi en champ lointain permet de déterminer la qualité de focalisation de l’onde incidente dans le milieu.
Cependant, cette approche ne fournit qu’une estimation globale et moyenne de la résolution d’une image échographique car elle nécessite de moyenner statistiquement les corrélations du champ réfléchi sur un grand nombre de réalisations du désordre, i.e sur un grand nombre de points de focalisation de l’onde incidente. Elle ne permet pas d’obtenir une évaluation précise et locale de la qualité de focalisation en chaque point de l’image. Par ailleurs, cette approche n’est valable qu’en régime de diffusion simple.
La demande de brevet WO-2010/001027 propose un procédé de sondage par ultrasons capable de séparer la composante de diffusion multiple de la composante de diffusion simple par filtrage d’une matrice de transfert fréquentielle représentative des réponses entre les transducteurs de l’ensemble des transducteurs. Ce procédé permet d’obtenir une information sur la diffusion multiple dans laquelle l’onde réfléchie résulte de plusieurs réflexions successives sur des diffuseurs, et pour laquelle le temps de vol n’est pas directement lié à la distance entre un diffuseur et les transducteurs.
Cependant, cette approche ne permet d’obtenir que l’évolution en profondeur des intensités de diffusion multiple et de diffusion simple. Elle ne permet pas d’accéder aux variations latérales de ces quantités et ne permet donc pas d’accéder pour chaque point de l’image au ratio entre les intensités de diffusion multiple et de diffusion simple. Ainsi, on ne peut obtenir qu’une information moyenne à chaque profondeur sur la fiabilité de l’image échographique, fondée sur le fait que la diffusion simple est prédominante.
Il est donc nécessaire de proposer une méthode palliant chacun des inconvénients précités.
RESUME
La présente description a pour objet, selon un premier aspect, un procédéde caractérisation ultrasonore non invasive d’un milieu hétérogène, le procédé comprenant :
- une étape de génération d’une série d’ondes ultrasonores incidentes dans une zone du milieu hétérogène générée au moyen de transducteurs, ladite série d’ondes ultrasonores incidentes constituant une base d’émission (i) ;
- une étape d’enregistrement d’une matrice de réflexion expérimentaleR ui (t) définie entre la base d’émission (i) en entrée et une base de réception (u) en sortie ;
- une étape de détermination d’une réponse REP(r, Δ r) du milieu entre un transducteur virtuel d’entrée (TVin) de position spatialer in calculé à partir d’une focalisation en entrée de la matrice de réflexion expérimentale correspondant à une tache focale d’entrée autour d’un premier point (P1) et un transducteur virtuel de sortie (TVout) de position spatialer out calculé à partir d’une focalisation en sortie de la matrice de réflexion expérimentale correspondant à une tache focale de sortie autour d’une deuxième point (P2), ladite réponse étant exprimée en fonction d’un point central (PC) de position spatialerdans le milieu, ledit point central (PC) étant situé au milieu des premier et deuxième points (P1, P2) et étant l’origine d’un axe de mesure (AXm) passant par lesdits premier et deuxième point, ledit axe de mesure formant un angle β par rapport à un premier axe (X) du milieu, et le premier point (P1) étant à une coordonnée d’éloignement+ Δ rsur l’axe de mesure, et le deuxième point (P2) étant à une coordonnée d’éloignement- Δ rsur l’axe de mesure.
Grâce à ces dispositions, le procédé permet, de manière très avantageuse, de sonder très localement le milieu en toute direction de l’axe de mesure passant par le premier point et le deuxième point, pour déterminer par les focalisations en entrée et en sortie, une nouvelle matrice de réponses locales REP(r, Δ r) (très riche en informations locales) en tout point du milieu de positionret pour toute direction angulaire β d’analyse.
Cette méthode permet alors avantageusement de définir des paramètres caractéristiques qui donnent des informations locales sur le milieu, paramètres très utiles pour quantifier la qualité de l’image échographique, et qui peuvent être utilisés pour optimiser ces images par calcul sans nécessiter d’itérer de nouvelles émissions et/ou acquisitions, ce qui est un avantage important lors de mesures in vivo.
La matrice de réflexion expérimentaleR ui (t) enregistrée peut être une matrice « réelle », c’est-à-dire composée de coefficients réels dans le domaine temporel, les signaux électriques enregistrés par chacun des transducteurs étant des nombres réels. En variante, cette matrice peut être une matrice « complexe », c’est-à-dire composée de valeurs complexes, par exemple dans le cas d’une démodulation pour une formation de voies en phase et en quadrature (connu en langue anglaise sous la dénomination « beamforming IQ »).
Dans divers modes de réalisation du procédé selon la présente divulgation, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes.
Selon un aspect, le procédé comprend en outre :
- une étape de détermination d’un profil de réponses PR(δr) qui est une pluralité de réponses REP(r, Δ r) calculées pour une pluralité de valeurs de la coordonnée d’éloignement Δ ret pour un même point central (PC), et pour un même axe de mesure (AXm), correspondant à un angle β prédéterminé, δr étant l’éloignement du deuxième point par rapport au point central, c’est-à-dire la valeur telle queΔ r= δr.u β ,u β étant un vecteur unitaire dans la direction de l’axe de mesure AXmdéfini par l’angle β..
Selon un aspect, le procédé comprend en outre :
- une étape de détermination de la résolution w(r) du point central (PC) à partir d’un module du profil de réponses, dans laquelle la résolution w(r) est une largeur du pic dudit module du profil de réponses PR(δr), centré autour de la coordonnée d’éloignement nulle (|Δ r|= 0).
Selon un aspect du procédé, la largeur du pic est estimée à une hauteur qui est une portion de la hauteur maximale dudit pic, ladite portion étant par exemple la moitié de la hauteur maximale du pic.
Selon un aspect, le procédé comprend en outre :
- une étape de détermination d’un critère de focalisation F(r) du point central (PC) à partir de la résolution w(r) et d’une résolution théorique w0(r), ladite résolution théorique w0(r) étant déterminée à partir de la base d’émission (i) en entrée et la base de réception (u) en sortie.
Selon un aspect, le procédé comprend en outre :
- une étape d’optimisation d’image dans laquelle on calcule le critère de focalisation F(r) en une pluralité de point du milieu, et on optimise au moins un paramètre de calcul de focalisation en entrée et/ou de focalisation en sortie en minimisant ou en maximisant une moyenne dudit critère de focalisation F(r) pour ladite pluralité de points.
Selon un aspect du procédé, ledit au moins un paramètre de calcul comprend la vitesse du son dans le milieu.
Selon un aspect du procédé, la résolution théorique est déterminée par une technique comprise dans la liste suivante :
- un premier calcul analytique au point central (PC) pour une pulsation (ω1) et la base d’émission (i) et la base de réception (u), dans lequel la résolution théorique est calculée par l’angle depuis lequel est vu le réseau de transducteur depuis le point central (PC),
- un deuxième calcul analytique au point central (PC) pour une plage de pulsation (Δω), la base d’émission (i) et la base de réception (u), dans lequel la résolution théorique est un calcul intégral sur ladite plage de pulsation de l’angle depuis lequel est vu le réseau de transducteur depuis le point central (PC) pondéré par le spectre fréquentiel de signaux de la matrice de réflexion expérimentaleR ui (t), et
- un troisième calcul de simulation de propagation des ondes, premièrement entre le premier point du milieu correspondant à un transducteur virtuel d’entrée (TVin) et la base d’émission (i), et deuxièmement entre le deuxième point du milieu correspondant à un transducteur virtuel de sortie (TVout) et la base de réception (u), ladite simulation utilisant la réponse REP(r, Δ r) et un modèle de propagation des ondes dans le milieu.
Selon un aspect, le procédé comprend en outre :
- une étape de détermination d’un taux de symétrie α(r) qui est le coefficient de corrélation moyen entre deux réponses réciproques, ladite moyenne étant effectuée pour des valeurs de coordonnées d’éloignement de module supérieures à une résolution prédéterminée wd(r) et/ou étant effectuée pour une plage de valeurs d’angle β ou pour une valeur d’angle prédéterminé βd,
une première réponse REP1(r, Δ r) du milieu entre un transducteur virtuel d’entrée (TVin) de position spatialer in calculé à partir d’une focalisation en entrée de la matrice de réflexion expérimentale correspondant à une tache focale d’entrée autour d’un premier point (P1) et un transducteur virtuel de sortie (TVout) de position spatialer out calculé à partir d’une focalisation en sortie de la matrice de réflexion expérimentale correspondant une tache focale de sortie autour d’un deuxième point (P2), et
une deuxième réponse REP2(r, - Δ r) du milieu entre un transducteur virtuel d’entré (TVin) de position spatialer in calculé à partir d’une focalisation en entrée de la matrice de réflexion expérimentale correspondant à une tache focale d’entrée autour du deuxième point (P2) et un transducteur virtuel de sortie (TVout) de position spatialer out calculé à partir d’une focalisation en sortie de la matrice de réflexion expérimentale correspondant à une tache focale de sortie autour du premier point (P1).
Selon un aspect du procédé, le taux de symétrie est calculé par la formule suivante :
ou par la formule suivante :
dans lesquelles
Re[.] est l’opérateur mathématique de partie réelle,
|.| est l’opérateur mathématique de module,
<.> est l’opérateur mathématique de moyenne, et
* est l’opérateur de conjugaison des nombres complexes.
Selon un aspect, le procédé comprend en outre :
- une étape de détermination d’un premier indicateur de diffusion multipleε (r)qui est calculé par :
α(r) étant le taux de symétrie défini au point central (PC) du milieu.
Selon un aspect, le procédé comprend en outre :
- une étape de détermination d’une intensité afocale Ioff(r) qui est la moyenne d’un module au carré des réponses REP(r, Δ r), la moyenne étant effectuée pour des valeurs d’éloignement de module supérieures à une résolution prédéterminée wd(r) et/ou étant effectuée pour une plage de valeurs d’angle β ou pour une valeur d’angle prédéterminé βd, par exemple calculé par la formule suivante :
- une étape de détermination d’une intensité de diffusion multiple IM(r) qui est le produit entre le taux de symétrie α(r) et l’intensité afocale Ioff(r), c’est-à-dire :
- une étape de détermination d’une intensité de bruit IN(r) étant le produit entre un moins le taux de symétrie α(r) et l’intensité afocale Ioff(r), c’est-à-dire :
de telle sorte que l’on ait la relation suivante :
Selon un aspect du procédé, le premier indicateur de diffusion multiple ε(r) est calculé par un ratio entre l’intensité de diffusion multiple IM(r) et l’intensité de bruit IN(r), c’est-à-dire :
Selon un aspect, le procédé comprend en outre :
- une étape de détermination d’une intensité confocale Ion(r) qui est la valeur d’un module au carré de la réponse REP(r,|Δ r| = 0) pour une valeur de coordonnée d’éloignement nulle (|Δ r|= 0), c’est-à-dire pour un point du milieu pour lequel le premier point (P1), le deuxième point (P2) et le point central (PC) sont confondus,
- une étape de détermination d’une intensité de diffusion simple IS(r) calculé à partir de l’équation suivante :
Selon un aspect, le procédé comprend en outre :
- une étape de détermination d’un deuxième indicateur de diffusion multiple γ(r) qui est calculé par :
Selon un aspect, le procédé comprend en outre :
- une étape de détermination d’une image d’un paramètre de caractérisation locale du milieu, ledit paramètre de caractérisation locale étant déterminé à partir de la réponse REP(r, Δ r).
Selon un aspect du procédé, le paramètre de caractérisation locale est choisi dans une liste comprenant la résolution w(r), le critère de focalisation F(r), le taux de symétrie α(r), le premier indicateur de diffusion multiple ε(r), le deuxième indicateur de diffusion multiple γ(r), l’intensité afocale Ioff(r), l’intensité confocale Ion(r), l’intensité de diffusion multiple IM(r), l’intensité de diffusion simple IS(r), l’intensité de bruit IN(r).
La présente description concerne, selon un deuxième aspect, un systèmede caractérisation ultrasonore non invasive d’un milieu hétérogène configuré pour la mise en œuvre de l’ensemble des exemples de procédés de caractérisation ultrasonore tels que décrits précédemment. Le système de caractérisation ultrasonore selon le deuxième aspect comprend :
-un premier réseau (10) de transducteurs adaptés pour générer une série d’ondes ultrasonores incidentes dans une zone du milieu hétérogène, ladite série d’ondes ultrasonores incidentes constituant une base d’émission (i), et pour enregistrer en fonction du temps les ondes ultrasonores rétrodiffusées par ladite zone ; et
- une unité de calcul (42) couplée au premier réseau de transducteurs et adaptée pour :
- enregistrer d’une matrice de réflexion expérimentaleR ui (t) définie entre la base d’émission (i) en entrée et une base de réception (u) en sortie ;
- déterminer une réponse REP(r, Δ r) du milieu entre un transducteur virtuel d’entrée (TVin) de position spatialer in calculé à partir d’une focalisation en entrée de la matrice de réflexion expérimentale correspondant à une tache focale d’entrée autour d’un premier point (P1) et un transducteur virtuel de sortie (TVout) de position spatialer out calculé à partir d’une focalisation en sortie de la matrice de réflexion expérimentale correspondant à une tache focale de sortie autour d’une deuxième point (P2), ladite réponse étant exprimée en fonction d’un point central (PC) de position spatiale r dans le milieu, ledit point central (PC) étant situé au milieu des premier et deuxième points (P1, P2) et étant l’origine d’un axe de mesure (AXm) passant par lesdits premier et deuxième point, ledit axe de mesure formant un angle β par rapport à un premier axe (X) du milieu, et le premier point (P1) étant à une coordonnée d’éloignement+ Δ rsur l’axe de mesure, et le deuxième point (P2) étant à une coordonnée d’éloignement- Δ rsur l’axe de mesure.
Le système de caractérisation selon la présente description peut comprendre au moins un réseau de transducteurs à la fois émetteur et récepteur, ou plusieurs réseaux de transducteurs, certains étant dédiés à l’émission, d’autres à la réception des ondes ultrasonores.
D’autres avantages et caractéristiques de la technique présentée ci-dessus apparaîtront à la lecture de la description détaillée ci-dessous, présentée de manière non limitative à des fins d’illustration, faite par référence aux figures dans lesquelles :
Lesfigures 1A à 1C(déjà décrites) illustrent des mécanismes d’émission/réception connus pour l’imagerie et la quantification ultrasonore ;
Lafigure 2(déjà décrite) illustre l’impact des aberrations en imagerie ultrasonore, selon l’art antérieur ;
Lafigure 3illustre un exemple de système de caractérisation ultrasonore pour la mise en œuvre des procédés de caractérisation ultrasonore selon la présente description ;
Lafigure 4illustre les définitions utilisées dans le procédé de caractérisation ultrasonore selon la présente description ;
Lafigure 5illustre un exemple d’image représentant le module d’une matrice de réponse REP(r, Δ r) selon le procédé de caractérisation ultrasonore tel qu’en figure 4 ;
Lafigure 6illustre un exemple de profil de réponse PR(δr) correspondant à la matrice de réponse de la figure 5 ;
Lafigure 7illustre à des images échographiques de trois milieux hétérogènes, la partie A de cette figure correspondant à un milieu de test dit « phantom », la partie B de cette figure correspondant à un milieu ayant une couche de viande posée sur le milieu de test, et la partie C de cette figure correspondant à un milieu qui est un foie « in vivo » ;
Lafigure 8illustre une image de résolution w(r) (à gauche) et une image de résolution théoriquew 0 (r)(à droite) établies pour le milieu de test (de la partie A) de la figure 7 ;
Lafigure 9illustre des images du critère de focalisation F(r) établies pour les trois milieux hétérogènes de la figure 7 (respectivement parties A, B, et C) ;
Lafigure 10illustre des images du premier indicateur de diffusion multiple ε(r) de focalisation établies pour les trois milieux hétérogènes de la figure 7 (respectivement parties A, B, et C) ;
Lafigure 11illustre des images du deuxième indicateur de diffusion multiple γ(r) établies pour les trois milieux hétérogènes de la figure 7 (respectivement parties A, B, et C) ;
Lafigure 12illustre trois calculs de vitesses de son optimales, réalisées à partir de trois modèles du milieu B de la figure 7 et à partir du critère de focalisation F(r) défini.
Dans les différents modes de réalisation qui vont être décrits par référence aux figures, des éléments semblables ou identiques portent les mêmes références.
Dans la description détaillée qui suit, seuls certains modes de réalisation sont décrits en détail pour assurer la clarté de l’exposé mais ces exemples ne visent pas à limiter la portée générale des principes ressortant de la présente description.
Les différents modes de réalisation et aspects décrits dans la présente description peuvent être combinés ou simplifiés de multiples manières. En particulier, les étapes des différents procédés peuvent être répétées, interverties, et/ou exécutées en parallèle, sauf précision contraire.
Lafigure 3illustre un exemple d’un système 40 de caractérisation ultrasonore pour la mise en œuvre des procédés de caractérisation ultrasonore d’un milieu hétérogène 20, selon la présente description. Le système 40 comprend au moins un premier réseau 10 de transducteurs 11, par exemple un réseau linéaire ou bidimensionnel ; les transducteurs sont par exemple des transducteurs piézoélectriques ultrasonores pouvant se présenter classiquement sous la forme d'une barrette rigide mise en contact avec le milieu 20. Le réseau de transducteurs fait par exemple partie d’un dispositif de sondage 41 ; le réseau de transducteurs est connecté à une unité de calcul 42, qui peut elle-même être reliée à un dispositif d’affichage 43 ; l’unité de calcul émet et enregistre des signaux électriques vers /à partir de chacun des transducteurs 11. Les transducteurs ultrasonores transforment ensuite ces signaux électriques en ondes ultrasonores et inversement. L’unité de calcul 42 est configurée pour la mise en œuvre d’étapes de calcul ou traitement notamment pour la mise en œuvre d’étapes de procédés selon la présente description. Par convention, on définit un repère spatial du milieu 20, en prenant un premier axe X et un deuxième axe Z perpendiculaire à celui-ci. Par simplification, le premier axe X correspond à la direction dans laquelle les transducteurs 11 sont alignés pour un réseau linéaire, et le deuxième axe Z correspond à la profondeur du milieu 20 par rapport à ce réseau 10 de transducteurs 11. Cette définition peut être étendue à un repère spatial à trois axes dans le cas d’un réseau 10 bidimensionnel.
Dans lafigure 3comme dans la suite de la description, il est fait référence à un réseau de transducteurs pour l’émission et la réception, étant bien entendu que, dans un cas plus général, plusieurs réseaux de transducteurs pourront être utilisés simultanément. Ils pourront être à la fois émetteur et récepteur, ou bien seulement émetteur pour certains et seulement récepteur pour d’autres.
Lorsque dans la présente description, il est fait référence à des étapes de calcul ou traitement pour la mise en œuvre notamment d’étapes de procédés, il est entendu que chaque étape de calcul ou traitement peut être mis en œuvre par logiciel, hardware, firmware, microcode ou toute combinaison appropriée de ces technologies. Lorsqu’un logiciel est utilisé, chaque étape de calcul ou traitement peut être mise en œuvre par des instructions de programme d’ordinateur ou du code logiciel. Ces instructions peuvent être stockées ou transmises vers un support de stockage lisible par un ordinateur (ou unité de calcul) et/ou être exécutées par un ordinateur (ou unité de calcul) afin de mettre en œuvre ces étapes de calcul ou traitement.
Définition de l’analyse d’un point du milieu
La présente description décrit des procédés et systèmes de caractérisation ultrasonore non invasive d’un échantillon hétérogène. Ces procédés et systèmes sont basés sur des définitions représentées enfigure 4: On définit un point central PC de position spatialerdans le repère spatial du milieu, et situé au milieu d’un premier point P1 et d’un deuxième point P2. On définit un axe de mesure AXmpassant par le premier point P1 et le deuxième point P2, et formant un angle β par rapport au premier axe X du réseau de transducteurs 11. Le point central PC est situé sur l’origine de l’axe de mesure AXm(coordonnée d’éloignement nulle sur l’axe de mesure). Le premier point P1 est à une coordonnée d’éloignement- Δ ret le deuxième point P2 est à une coordonnée d’éloignement+ Δ rdu point central PC, origine de l’axe de mesure.
La position spatialeret la coordonnée d’éloignementΔ rsont notée en gras, signifiant que ces éléments sont des vecteurs de position et de décalage par rapport à une position, vecteurs pris dans le repère spatial du milieu (X, Z). Le vecteur de coordonnée d’éloignementΔ rprend ainsi en compte la direction de l’axe de mesure AXm, et son angle β par rapport au premier axe X. D’autres définitions des positions des points relativement aux autres sont possibles et accessibles à tout spécialiste du métier des ultrasons. Notamment, les premier et deuxièmes points peuvent être repérés par une distance |Δ r| et l’angle β, ou par un autre repérage.
Comme représenté enfigure 4, le procédé de caractérisation ultrasonore non invasive comprend :
- une étape de génération d’une série d’ondes ultrasonores incidentes USindans une zone dudit milieu hétérogène, au moyen d’un réseau 10 de transducteurs 11, ladite série d’ondes ultrasonores incidentes étant une base d’émissioni; et
- une étape d’enregistrement d’une matrice de réflexion expérimentaleR ui (t) définie entre la base d’émissionien entrée et une base de réceptionuen sortie ;
- une étape de détermination d’une réponse REP(r, Δ r) du milieu entre un transducteur virtuel d’entrée TVinde position spatialer in au premier point P1 et un transducteur virtuel de sortie TVoutde position spatialer out au deuxième point P2, ladite réponse étant exprimée en fonction du point central PC de position spatialer.
Le point central PC étant au milieu des deux points P1, P2, on a les relations suivantes :
La base d’émissionien entrée étant par exemple une base d’ondes générées chacune par un seul des transducteurs 11 du réseau 10 ou une base d’ondes planes d’inclinaison angulaire θ par rapport à l’axe X, comme décrit précédemment dans la description desfigures 1A à 1C.
La base de réceptionuest usuellement la base des transducteurs 11.
Dans l’étape de détermination d’une réponse REP(r, Δ r), on applique :
- un processus de focalisation en entrée à partir de la matrice de réflexion expérimentaleR ui (t) créant une tache focale d’entrée autour du premier point P1, ladite tache focale d’entrée correspondant au transducteur virtuel d’entrée TVin, et
- un processus de focalisation en sortie à partir de la matrice de réflexion expérimentaleR ui (t) créant une tache focale de sortie autour du deuxième point P2, ladite tache focale de sortie correspondant au transducteur virtuel de sortie TVout.
Le premier point P1 étant relatif au transducteur virtuel d’entrée TVin, il se situe alors à une coordonnée -Δ rsur l’axe de mesure AXmpar rapport au point central PC et le deuxième point P2 étant relatif au transducteur virtuel de sortie TVout, il se situe alors à une coordonnée+ Δ rsur l’axe de mesure AXmpar rapport au point central PC.
Par exemple, un calcul de la réponse REP(r, Δ r) du milieu entre le transducteur virtuel d’entrée TVinet le transducteur virtuel de sortie TVoutpar un processus de focalisation, qui est par exemple un procédé de formation de voie amélioré, qui peut être exprimé par la formule simplifiée suivante :
(Equ. 1)
dans laquelle
Ninest le nombre d’éléments de la base d’émissioni,
Noutest le nombre d’éléments de la base de réceptionuen sortie,
R ui (t) est la matrice de réflexion expérimentale, dont est l’élément de la matrice de réflexion expérimentaleR ui (t) enregistré par le transducteuru out consécutif à l’émissioni in au temps τ.
Le temps τ est la somme du temps de vol à l’aller τinde l’onde ultrasonore entre les transducteurs de la base d’émissioniet le premier point P1 et le temps de vol au retour τoutde l’onde ultrasonore entre le deuxième point P2 et les transducteurs de la base de réceptionu, comme explicité par la formule suivante :
(Equ. 2)
Grâce à ces dispositions, le procédé permet de sonder très localement le milieu en toute direction correspondant à l’axe de mesure AXm, pour en extraire par les focalisations en entrée et en sortie plus d’information locale sur le milieu au point central PC de position spatialer, entre le premier point et deuxième point du milieu hétérogène 20. Cette information locale est entièrement contenue dans les valeurs de la réponse calculée, la réponse REP(r, Δ r) du milieu que l’on peut exploiter pour caractériser chaque point du milieu par exemple en terme de résolution ou en terme de diffusion multiple.
Les réponses REP(r, Δ r) peuvent ainsi être déterminées pour tout ensemble de valeurs de distance d’éloignement |Δ r|réelle, par exemple entre deux bornes telles que - Δrmaxet + Δrmax, ces bornes étant déterminées pour que les transducteurs virtuels en entrée et en sortie restent à l’intérieur du milieu 20. (|Δ r|, β) sont alors les coordonnées polaires du vecteur coordonnées d’éloignementΔ r.
Dans cette précédente convention de description spatiale autour du point central PC de position spatialer, la réponse REP(r, - Δ r) correspond à inverser les positions spatiales des transducteurs virtuels d’entrée et de sortie.
L’ensemble des réponses REP(r, Δ r) peut alors être enregistré dans une matrice de même nom.
Lafigure 5montre une image correspondant à une sous-matrice de la matrice de réponseREP, ladite sous-matrice correspondant à un ensemble de plusieurs points centraux PC de position spatialerdans laquelle la coordonnée d’axe Z est fixée, et l’angle β est nul. Ainsi, dans cette image l’abscisse correspond à une variation selon l’axe X de la position du point central PC et l’ordonnée correspond à la coordonnée d’éloignementΔ rpar rapport à ce point central. Les valeurs des points de cette image (réponse) en dehors de l’axe des abscisses de cette image ont une valeur faible (mais non nulle). Les valeurs des points de cette image (réponse) sur l’axe des abscisses ont une valeur correspondant à l’intensité du point d’image ultrasonore au point central PC.
Une telle image peut être extraite de la réponse REP(r, Δ r) pour des variations de coordonnées d’éloignementΔ rsur un seul axe de mesure AXm ou plusieurs axes de mesures, c’est-à-dire pour une ou plusieurs valeurs d’angle β.
Une image polaire représentant la variation du module de la réponse en fonction de la distance d’éloignement |Δ r|et de l’angle β peut également être construite, ce qui donne une représentation de la variation de la réponse autour d’un point central PC, et donc de de la tâche focale en ce point.
Calcul d’un profil de réponses PR(δr)
Ayant obtenu les réponses du milieu déterminées selon le procédé ci-dessus, on peut effectuer une étape de détermination d’un profil de réponses PR(δr), le profil de réponses étant une pluralité des réponses REP(r, Δ r) calculées pour une pluralité de valeurs de la coordonnée d’éloignementΔ r. Ce profil de réponses PR(δr) est considéré pour un même point central PC de position spatialeret selon un même axe de mesure AXm, correspondant à une même direction d’angle β. La profil de réponse PR(δr) est donc déterminé pour une pluralité d’éloignements δr, l’éloignement δr étant l’abscisse du deuxième point P2 par rapport au point central PC, c’est-à-dire la valeur telle queΔ r= δr.u β ,u β étant un vecteur unitaire dans la direction de l’axe de mesure AXmdéfini par l’angle β. Autrement dit, le profil de réponses PR(δr) est une coupe verticale de l’image de lafigure 5, et ce profil de réponses est la courbe de ladite coupe d’image.
Les réponses REP(r, Δ r) pouvant être des valeurs complexes notamment dans le cas d’utilisation de formulation de focalisation en valeurs complexe, comme cela est connu dans des formations de voies en phase et en quadrature (connu en langue anglaise sous la dénomination « beamforming IQ »). En conséquence, le profil de réponses PR(δr) pourra également être représenté par tout module de ces réponses complexes.
Lafigure 6montre un exemple schématique du carré du module d’un profil de réponses PR(δr) que l’on peut déterminer pour un point central PC du milieu et pour une direction d’angle β : ici pour β = 0 car l’image de la figure 5 est également établie pour cette valeur d’angle.
Cependant, l’angle β peut prendre toute valeur entre zéro (0) et pi (π), et donc une courbe PRcde profil de réponse peut être tracée ou déterminée pour plusieurs valeurs d’angle β.
L’ensemble des profils de réponses PR(δr) ou PR(δr, β) (si l’on prend plusieurs angles, mais nous ne garderons que la position spatiale dans la suite par soucis de simplification de l’exposé) peut être enregistré dans une matrice de même nom.
La sous-matrice représentée enfigure 5est donc l’ensemble des profils de réponses PR(δr) pour un ensemble de point central PC de position spatialerdans laquelle la coordonnée d’axe Z est fixée, et l’angle β est nul (ou constant).
Ce profil de réponses PR(δr) est une représentation de base permettant de déterminer de nouveaux paramètres de caractérisation locale du milieu et/ou de la performance du procédé d’imagerie par ultrason (i.e. de la formation de voie). Nous illustrerons les résultats de ces paramètres de caractérisation par des images du milieu hétérogène, construites avec lesdits paramètres de caractérisation.
Ces paramètres de caractérisation pour trois cas différents de milieux hétérogènes sont exposés plus en détail ci-dessous. Lafigure 7montre des images échographiques de ces trois milieux hétérogènes :
- sur la partie A de cette figure 7 (à gauche), on aperçoit une image d’un milieu de test (milieu A) dit « phantom », qui comprend deux cylindres de rigidité différentes,
- sur la partie B de cette figure 7 (au centre), on aperçoit une image d’un milieu ayant une couche de viande posée sur le milieu de test précédent (milieu B), et
- sur la partie C de cette figure 7 (à droite), on aperçoit une image d’un milieu qui est un foi « in vivo » (milieu C).
Calcul de la résolution du point w(r)
On peut alors effectuer une étape de détermination de la résolution w(r) du point central PC dans la direction de l’axe de mesure AXmd’angle β, à partir d’un module du profil de réponses.
On remarque que le module du profil de réponse comme représenté enfigure 6comprend un pic ou maximum autour de la coordonnée d’éloignementΔ rnulle (| Δ r|= 0). La résolution w(r) peut alors être déterminée par exemple par la largeur de ce pic. La largeur de ce pic est par exemple déterminée à une hauteur qui est une portion de la hauteur maximale de ce pic. Par exemple, cette portion de la hauteur sera la moitié ou un tier (1/3) ou deux tiers (2/3) ou toute autre ratio de la hauteur maximale. La hauteur maximale du pic est en fait l’intensité de l’image échographique au point central PC si on regarde seulement le module au carré du profil de réponse, i.e. |PR(δr=0)|², comme cela est le cas sur l’exemple illustré en figure 6.
On comprend que la résolution dépend du point central PC considéré, mais également de l’angle β.
Par conséquent, la méthode proposée permet d’obtenir en chaque point :
- la résolution axiale pour un angle β de π/2, i.e. par la valeur w(r) prise pour un profil de réponse correspondant à un angle β=π/2, et
- la résolution latérale pour un angle β nul, i.e. par la valeur w(r) prise pour un profil de réponses correspondant à un angle β=0.
La méthode permet de définir en tout point du milieu l’étendu de la tache focale et donc la résolution du procédé d’échographie dans chacune des directions d’angle β.
L’image enpartie gauche de la figure 8montre un exemple de calcul de la résolution w(r,β) en chaque point du milieu de test (milieu A). On remarque que la résolution se dégrade avec la profondeur et lorsque l’on se déplace vers le bord de l’image.
Calcul de la résolution théorique w0(r)
Selon une première variante, la résolution théorique w0(r) est déterminée par un premier calcul analytique au point central (PC) pour une pulsation (ω1), la base d’émission (i) et la base de réception (u) : Elle est calculée par l’angle depuis lequel est vu le réseau de transducteur depuis le point central (PC).
Selon une deuxième variante, la résolution théorique w0(r) est déterminée par un deuxième calcul analytique au point central (PC) pour une plage de pulsation (Δω) et la base d’émission (i) et la base de réception (u). Elle est obtenue par un calcul intégral sur ladite plage de pulsation et sur l’angle depuis lequel est vu le réseau de transducteur depuis le point central (PC) pondéré par le spectre fréquentiel de signaux enregistrés. Ce dernier peut être obtenu en moyennant le module de la transformée de Fourier des éléments de la matrice de réflexion expérimentaleR ui (t).
Selon une troisième variante, la résolution théorique w0(r) est déterminée par un troisième calcul de simulation de propagation des ondes, premièrement entre le premier point du milieu correspondant à un transducteur virtuel d’entrée (TVin) et la base d’émission (i), et deuxièmement entre le deuxième point du milieu correspondant à un transducteur virtuel de sortie (TVout) et la base de réception (u), ladite simulation utilisant la réponse REP(r, Δ r) et un modèle de propagation des ondes dans le milieu.
L’image enpartie droite de la figure 8montre un exemple de calcul de la résolution théorique w0(r) en chaque point du milieu de test (milieu A). On remarque que la résolution se dégrade également avec la profondeur et lorsque l’on se déplace vers le bord de l’image. On remarque que cette image est très similaire à l’image en partie gauche de cette même figure. Le calcul de la résolution effectué par le procédé précédent est donc en accord avec le calcul de la résolution théorique.
Calcul d’un critère de focalisation F(r)
On peut alors effectuer une étape de détermination d’un critère de focalisation F(r) du point central PC à partir de la résolution w(r) et d’une résolution théorique w0(r). La résolution théorique est par exemple déterminée à partie d’une base d’émissionien entrée, une base de réceptionuen sortie et d’une modélisation de la propagation des ondes ultrasonores dans le milieu.
Usuellement, le critère de focalisation F(r) est un ratio desdites résolution et résolution théorique, ou l’inverse (une simple convention). C’est-à-dire, que nous pouvons obtenir :
(Equ. 3)
ou
(Equ. 4)
Lafigure 9illustre des images du critère de focalisation F(r) établies pour les trois milieux hétérogènes de la figure 7 (milieux A, B, et C en correspondance entre les figures 7 et 9).
Une valeur de un (1) pour ce critère de focalisation correspond à une résolution et une résolution théorique identiques (en clair sur la figure). Une valeur proche de zéro (0) pour ce critère de focalisation correspond des valeurs de résolutions divergentes (en sombre sur la figure), c’est-à-dire une focalisation dégradée.
L’image du milieu A illustre une grande homogénéité avec une moyenne de ce critère de focalisation proche de 0,97. Cela signifie que l’image échographique est bien formée et que les hypothèses de focalisation sont correctes. Les images des milieux B et C montrent des dégradations notables correspondant à des hétérogénéités situées en amont de la propagation des ondes ultrasonores : couche de viande en surface pour milieu B et les tissus adipeux ou musculaires pour le foie du milieu C. L’image du milieu C met en relief des zones très dégradées (zones sombres en bas à gauche de cette image) qui signifient que l’image réalisée en figure 7 (partie C) a également les mêmes problèmes aux mêmes endroits. Ici, l’image du critère de focalisation peut indiquer à l’opérateur du système échographique que la partie gauche de son image échographique est de mauvaise qualité, notamment en résolution spatiale. Cela peut l’aider à interpréter son image échographique pour établir un diagnostic préférentiellement sur les zones correctement imagées ou pour modifier sa façon de générer ladite image échographique. Par exemple, cela peut inciter l’opérateur à modifier le positionnement de la sonde (i.e. le réseau de transducteurs) de telle sorte à obtenir un critère de focalisation de bonne qualité dans la/les zone(s) d’intérêt pour le diagnostic médical.
Avantageusement, le système d’imagerie qui met en œuvre la présente technique sera apte à superposer une image échographique à une image du critère de focalisation.
Calcul d’un taux de symétrie α(r)
On peut alors effectuer une étape de détermination d’un taux de symétrie α(r) du point central PC, ce taux de symétrie étant un coefficient de corrélation moyen entre deux réponses réciproques : c’est-à-dire en échangeant les positionsr in etr out du transducteur virtuel en entrée TVinet du transducteur virtuel en sortie TVout. On teste ainsi le critère de réciprocité acoustique du milieu au point central PC de position spatialer.
Pour ce procédé, on détermine :
une première réponse REP1(r, Δ r) du milieu entre un transducteur virtuel d’entrée TVin(position spatialer in ) calculé à partir d’une focalisation en entrée de la matrice de réflexion expérimentale créant une tache focale d’entrée autour du premier point P1 et un transducteur virtuel de sortie TVout(position spatialer out ) calculé à partir d’une focalisation en sortie de la matrice de réflexion expérimentale créant une tache focale de sortie autour du deuxième point P2, et
une deuxième réponse REP2(r, - Δ r) du milieu entre un transducteur virtuel d’entré TVin(position spatialer in ) calculé à partir d’une focalisation en entrée de la matrice de réflexion expérimentale créant une tache focale d’entrée autour du deuxième point P2 et un transducteur virtuel de sortie TVout(position spatialer out ) calculé à partir d’une focalisation en sortie de la matrice de réflexion expérimentale créant une tache focale de sortie autour du premier point P1.
En outre, la moyenne du coefficient de corrélation moyen est effectuée pour des valeurs de coordonnées d’éloignementΔ rde module supérieures à une résolution prédéterminée wd(r) (comme représenté enfigure 6) et/ou étant effectuée pour une plage de valeurs d’angle β ou pour une valeur d’angle prédéterminé βd.
La résolution prédéterminée wd(r) est avantageusement une valeur supérieure à la moitié de la résolution w(r). De préférence, la résolution prédéterminée wd(r) est une valeur supérieure à une fois, deux fois ou trois fois la résolution w(r), pour mieux exclure les valeurs du pic, comme cela est visible en figure 6. Ainsi, la moyenne du coefficient de corrélation selon cette variable, signifie que l’on moyenne le coefficient de corrélation pour des valeurs de coordonnée d’éloignementΔ réloignées du pic, et donc pour des transducteurs virtuels éloignés du point central PC, ce qui permet en effet de tester la réciprocité acoustique de ce point central PC, i.e. la symétrie de la matrice de réponse REP(r, Δ r) par rapport à l’axe d’abscisse (voir exemple d’une telle matrice enfigure 5). Bien que les valeurs au-delà de la résolution prédéterminée soient faibles, la comparaison par une corrélation moyenne de la symétrie permet d’estimer de manière fiable un taux de symétrie α(r) du point central PC.
Avantageusement, la moyenne du coefficient de corrélation est effectuée pour une plage de valeurs de valeurs d’angle β ou pour une valeur d’angle prédéterminé βd. Ainsi, la moyenne du coefficient de corrélation selon cette variable, signifie que l’on moyenne le coefficient de corrélation pour une ou plusieurs valeurs d’angle β, ce qui permet de tester la symétrie angulaire ou selon un secteur angulaire du taux de symétrie autour du point central PC, i.e. la réciprocité de ce point central PC.
Avantageusement, la moyenne du coefficient de corrélation est effectuée pour des valeurs de coordonnées d’éloignementΔ rde module supérieur à une résolution prédéterminée wd(r).et pour une plage de de valeurs d’angle β. Ainsi, la corrélation est moyennée sur des valeurs de coordonnées d’éloignement pour lesquelles la contribution de diffusion simple est nulle, cette dernière n’apparaissant que pour un module de la coordonnée d’éloignement inférieure à la résolution prédéterminée.
Cela permet aussi d'obtenir une estimation plus robuste et plus stable du taux de symétrie α(r).
Par exemple, le taux de symétrie α(r) peut être calculé par la formule de corrélation suivante :
(Equ. 5)
ou par exemple par la formule de corrélation suivante :
(Equ. 6)
dans lesquelles
Re[.] est l’opérateur mathématique de partie réelle,
|.| est l’opérateur mathématique de module,
<.> est l’opérateur mathématique de moyenne, cet opérateur pouvant être effectué selon une ou plusieurs variables (par exemple ci-dessus des valeurs de coordonnées d’éloignement supérieures à une résolution prédéterminée et des valeurs d’angle β), et
* est l’opérateur de conjugaison des nombres complexes.
Plusieurs formules de corrélations peuvent être utilisées et le spécialiste modifiera cette définition selon ses besoins et les caractéristiques du milieu observé.
D’une manière générale, le taux de symétrie α(r) proche de zéro (0) si la propagation des ondes ultrasonores ne se comporte pas de manière réciproque autour du point central PC, et le taux de symétrie α(r) proche de un (1) si la propagation des ondes ultrasonores se comporte de manière symétrique ou réciproque autour du point central PC.
Ainsi, ce taux de symétrie teste la validité de la réciprocité acoustique pour la partie des signaux correspondant à une diffusion multiple. Cela permet de discriminer la diffusion multiple du bruit, le bruit ne respectant pas la propriété de réciprocité acoustique.
Calcul d’un premier indicateur de diffusion multiple
ε
(r)
On peut alors effectuer une étape de détermination d’un premier indicateur de diffusion multiple ε(r) du point central PC, cet indicateur de diffusion multiple étant par exemple calculé par la formule suivante :
(Equ. 7)
dans laquelle
α(r) est le taux de symétrie défini au point central PC du milieu.
Ce premier indicateur de diffusion multiple ε(r) est nul si le taux de symétrie est nul, et il tend vers l’infini si le taux de symétrie est proche de un (1).
Lafigure 10illustre des images du premier indicateur de diffusion multiple ε(r) établies pour les trois milieux hétérogènes de la figure 7 (milieux A, B, et C en correspondance entre les figures 7 et 10). Ces figures sont graduées en décibel (échelle logarithmique).
Une valeur grande pour ce premier indicateur de diffusion multiple correspond à une proportion de diffusion multiple importante (en clair sur la figure). Une valeur faible pour ce critère de focalisation correspond à peu de diffusion multiple ou non significative (en sombre sur la figure). Ainsi, ce premier indicateur de diffusion multiple montre une zone localisée de diffusion des ondes ultrasonores derrière un des cylindres d’hétérogénéité pour le milieu A. Pour les milieux B et C, ce premier indicateur de diffusion multiple met en relief la diffusion dans une majorité du volume de ces milieux, soit en aval de la partie viande du milieu B, soit en aval des tissus adipeux du foie.
Calcul d’une intensité afocale Ioff(r)
On peut alors effectuer une étape de détermination d’une intensité afocale Ioff(r) du point central PC, cette intensité afocale étant une moyenne d’un module au carré des réponses REP(r, Δ r), la moyenne étant effectuée pour des valeurs d’éloignement de module supérieures à une résolution prédéterminée wd(r) (c’est-à-dire |Δ r | >wd(r)) et/ou étant effectuée pour une plage de valeurs d’angle β ou pour une valeur d’angle prédéterminée βd, par exemple calculé par la formule suivante :
(Equ. 8)
avec
(Equ. 9)
dans laquelle
I(r, Δ r) est l’intensité du point, i.e. le module au carré de la réponse en ce point,
|.| est l’opérateur mathématique de module, et
<.> est l’opérateur mathématique de moyenne, cet opérateur pouvant être effectué selon une ou plusieurs variables (par exemple ci-dessus des valeurs de coordonnées d’éloignement supérieures à une résolution prédéterminée et des valeurs d’angle β).
Calcul d’une intensité de diffusion multiple
I
M
(r)
On peut alors effectuer une étape de détermination d’une intensité de diffusion multiple IM(r), cette intensité de diffusion multiple étant le produit entre le taux de symétrie α(r) et l’intensité afocale Ioff(r), c’est-à-dire :
(Equ. 10)
Calcul d’une intensité de bruit IN(r)
On peut alors effectuer une étape de détermination d’une intensité de bruit IN(r), cette intensité de bruit étant le produit entre un moins le taux de symétrie α(r) et l’intensité afocale Ioff(r), c’est-à-dire :
(Equ. 11)
de telle sorte que l’on ait la relation suivante :
(Equ. 12)
Alors, le premier indicateur de diffusion multiple ε(r) peut être également calculé par un ratio entre l’intensité de diffusion multiple IM(r) et l’intensité de bruit IN(r), c’est-à-dire selon la formule suivante :
(Equ. 13)
Calcul d’une intensité confocale Ion(r)
On peut alors effectuer une étape de détermination d’une intensité confocale Ion(r), cette intensité confocale étant alors la valeur d’un module au carré de la réponse REP(r, Δ r=0)) pour une valeur de coordonnée d’éloignementΔ rnulle, i.e. |Δ r | = 0, c’est-à-dire pour un point du milieu pour lequel le premier point P1, le deuxième point P2 et le point central PC sont confondus.
Calcul d’une intensité de diffusion simple IS(r)
On peut alors effectuer une étape de détermination d’une intensité de diffusion simpleI S (r), cette intensité de diffusion simple étant calculée à partir de l’équation suivante :
(Equ. 14)
Calcul d’un deuxième indicateur de diffusion multiple
γ
(r)
On peut alors effectuer une étape de détermination d’un deuxième indicateur de diffusion multiple γ(r), ce deuxième indicateur de diffusion multiple étant calculé par la formule suivante :
(Equ. 15)
Lafigure 11illustre des images du deuxième indicateur de diffusion multiple γ(r) établies pour les trois milieux hétérogènes de la figure 7 (milieux A, B, et C en correspondance entre les figures 7 et 10). Ces figures sont graduées en décibel (échelle logarithmique).
Une valeur grande pour ce premier indicateur de diffusion multiple correspond à une proportion de diffusion multiple importante par rapport à la diffusion simple (en clair sur la figure). Une valeur faible pour ce critère de focalisation correspond à peu de diffusion multiple faible par rapport à la diffusion simple (en sombre sur la figure). Ainsi, ce deuxième indicateur de diffusion multiple montre que le milieu C est très diffusant.
Optimisation du calcul de l’image échographique
Une utilisation possible des calculs précédents est l’optimisation du calcul de l’image échographique.
L’image échographique est par exemple calculée par des focalisations successives en entrée et en sortie de l’ensemble des points du milieu. L’ensemble des points de l’axe d’abscisse de la figure 5 correspond à une ligne de cette image échographique calculée à la profondeur z utilisée pour former l’image de cette figure 5.
Cependant ce calcul d’image échographie dépend fortement d’une hypothèse de milieu homogène dans lequel la vitesse du son (la vitesse de propagation des ondes ultrasonores) est bien connue et constante. Si cette hypothèse est inexacte, les lois de retard de focalisation ne correspondent pas au milieu considéré et la focalisation est imparfaite. La résolution de l’image échographique est alors dégradée. D’autres paramètres de calcul de focalisation influencent éventuellement la focalisation.
On peut alors effectuer une étape d’optimisation d’image dans laquelle on calcule le critère de focalisation F(r) en une pluralité de point du milieu (par exemple une zone prédéterminée de l’image), et on optimise au moins un paramètre de calcul de focalisation en entrée et/ou de focalisation en sortie en minimisant ou en maximisant une moyenne dudit critère de focalisation F(r) pour ladite pluralité de points.
Par exemple, ledit au moins un paramètre de calcul comprend au moins la vitesse du son dans le milieu. Eventuellement, ce paramètre de calcul est la vitesse du son.
Dans le cas de l’optimisation du critère de focalisation F(r) à partir de d’une plage possible de vitesse du son c, on obtient une courbe comme une de celles présentées enfigure 12, et pour laquelle le critère de focalisation F(r) présente un maximum pour une vitesse du son optimale copt. Cette vitesse du son optimale coptcorrespond en fait à une vitesse du son intégrée sur toute l’épaisseur du milieu traversée (entre la sonde et le plan focal de la zone prédéterminée).
Pour améliorer l’estimation de la vitesse du son dans le milieu, il faut par exemple introduire un modèle du milieu multi couches, chaque couche étant supposée homogène avec une vitesse de son cide cette couche. Il faut estimer l’épaisseur de chaque couche, soit avec des connaissancesa prioridu milieu, soit à l’aide d’une première image échographique. On estime alors par optimisation du critère de focalisation la vitesse du son dans la couche la plus externe avec un modèle homogène, puis on estime par l’optimisation la couche suivante la vitesse du son dans la couche suivante avec un modèle à deux couches, et ainsi de suite.
La figure 12 illustre les optimisations effectuées pour une même zone prédéterminée le milieu B en utilisant :
a) un modèle à une couche en courbe B1 avec une vitesse du son optimale c1opt,
b) un modèle à deux couches en courbe B2 avec une vitesse du son optimale c2opt,
c) un modèle à trois couches en courbe B1 avec une vitesse du son optimale c3opt,
La méthode d’optimisation présentée ci-dessus et la modélisation multicouche du milieu permettent donc d’affiner les estimations de vitesse du son, et permettent donc de déterminer l’évolution de la vitesse du son dans la profondeur du milieu.
Claims (18)
- Procédé de caractérisation ultrasonore non invasive d’un milieu hétérogène, le procédé comprenant :
- une étape de génération d’une série d’ondes ultrasonores incidentes dans une zone du milieu hétérogène générée au moyen de transducteurs, ladite série d’ondes ultrasonores incidentes constituant une base d’émission (i) ;
- une étape d’enregistrement d’une matrice de réflexion expérimentaleR u i (t) définie entre la base d’émission (i) en entrée et une base de réception (u) en sortie ;
- une étape de détermination d’une réponse REP(r, Δ r) du milieu entre un transducteur virtuel d’entrée (TVin) de position spatialer in calculé à partir d’une focalisation en entrée de la matrice de réflexion expérimentale correspondant à une tache focale d’entrée autour d’un premier point (P1) et un transducteur virtuel de sortie (TVout) de position spatialer out calculé à partir d’une focalisation en sortie de la matrice de réflexion expérimentale correspondant à une tache focale de sortie autour d’une deuxième point (P2), ladite réponse étant exprimée en fonction d’un point central (PC) de position spatiale r dans le milieu, ledit point central (PC) étant situé au milieu des premier et deuxième points (P1, P2) et étant l’origine d’un axe de mesure (AXm) passant par lesdits premier et deuxième point, ledit axe de mesure formant un angle β par rapport à un premier axe (X) du milieu, et le premier point (P1) étant à une coordonnée d’éloignement+ Δ rsur l’axe de mesure, et le deuxième point (P2) étant à une coordonnée d’éloignement- Δ rsur l’axe de mesure. - Procédé selon la revendication 1, comprenant en outre :
- une étape de détermination d’un profil de réponses PR(δr) qui est une pluralité de réponses REP(r, Δ r) calculées pour une pluralité de valeurs de la coordonnée d’éloignement Δ ret pour un même point central (PC), et pour un même axe de mesure (AXm), correspondant à un angle β prédéterminé, δr étant l’éloignement du deuxième point par rapport au point central, c’est-à-dire la valeur telle queΔ r= δr.u β ,u β étant un vecteur unitaire dans la direction de l’axe de mesure AXmdéfini par l’angle β. - Procédé selon la revendication 2, comprenant en outre :
- une étape de détermination de la résolution w(r) du point central (PC) à partir d’un module du profil de réponses, dans laquelle la résolution w(r) est une largeur du pic dudit module du profil de réponses PR(δr), centré autour de la coordonnée d’éloignement nulle (δr = 0). - Procédé selon la revendication 3, dans lequel la largeur du pic est estimée à une hauteur qui est une portion de la hauteur maximale dudit pic, ladite portion étant par exemple la moitié de la hauteur maximale du pic.
- Procédé selon la revendication 2, comprenant en outre :
- une étape de détermination d’un critère de focalisation F(r) du point central (PC) à partir de la résolution w(r) et d’une résolution théorique w0(r), ladite résolution théorique w0(r) étant déterminée à partir de la base d’émission (i) en entrée et la base de réception (u) en sortie. - Procédé selon la revendication 5, comprenant en outre :
- une étape d’optimisation d’image dans laquelle on calcule le critère de focalisation F(r) en une pluralité de point du milieu, et on optimise au moins un paramètre de calcul de focalisation en entrée et/ou de focalisation en sortie en minimisant ou en maximisant une moyenne dudit critère de focalisation F(r) pour ladite pluralité de points. - Procédé selon la revendication 6, dans laquelle ledit au moins un paramètre de calcul comprend la vitesse du son dans le milieu.
- Procédé selon la revendication 5, dans laquelle la résolution théorique est déterminée par une technique comprise dans la liste comprise dans :
- un premier calcul analytique au point central (PC) pour une pulsation (ω1) et la base d’émission (i) et la base de réception (u), dans lequel la résolution théorique est calculée par l’angle depuis lequel est vu le réseau de transducteur depuis le point central (PC),
- un deuxième calcul analytique au point central (PC) pour une plage de pulsation (Δω) et la base d’émission (i) et la base de réception (u), dans lequel la résolution théorique est un calcul intégral sur ladite plage de pulsation de l’angle depuis lequel est vu le réseau de transducteur depuis le point central (PC) pondéré par le spectre fréquentiel de signaux de la matrice de réflexion expérimentaleR ui (t), et
- un troisième calcul de simulation de propagation des ondes, premièrement entre le premier point du milieu correspondant à un transducteur virtuel d’entrée (TVin) et la base d’émission (i), et deuxièmement entre le deuxième point du milieu correspondant à un transducteur virtuel de sortie (TVout) et la base de réception (u), ladite simulation utilisant la réponse REP(r, Δ r) et un modèle de propagation des ondes dans le milieu. - Procédé selon la revendication 1, comprenant en outre :
- une étape de détermination d’un taux de symétrie α(r) qui est le coefficient de corrélation moyen entre deux réponses réciproques, ladite moyenne étant effectuée pour des valeurs de coordonnées d’éloignement de module supérieures à une résolution prédéterminée wd(r) et/ou étant effectuée pour une plage de valeurs d’angle β ou pour une valeur d’angle prédéterminé βd,
une première réponse REP1(r, Δ r) du milieu entre un transducteur virtuel d’entrée (TVin) de position spatialer in calculé à partir d’une focalisation en entrée de la matrice de réflexion expérimentale correspondant à une tache focale d’entrée autour d’un premier point (P1) et un transducteur virtuel de sortie (TVout) de position spatialer out calculé à partir d’une focalisation en sortie de la matrice de réflexion expérimentale correspondant une tache focale de sortie autour d’un deuxième point (P2), et
une deuxième réponse REP2(r, - Δ r) du milieu entre un transducteur virtuel d’entré (TVin) de position spatialer in calculé à partir d’une focalisation en entrée de la matrice de réflexion expérimentale correspondant à une tache focale d’entrée autour du deuxième point (P2) et un transducteur virtuel de sortie (TVout) de position spatialer out calculé à partir d’une focalisation en sortie de la matrice de réflexion expérimentale correspondant à une tache focale de sortie autour du premier point (P1). - Procédé selon la revendication 9, dans lequel le taux de symétrie est calculé par la formule suivante :
ou par la formule suivante :
dans lesquelles
Re[.] est l’opérateur mathématique de partie réelle,
|.| est l’opérateur mathématique de module,
<.> est l’opérateur mathématique de moyenne, et
* est l’opérateur de conjugaison des nombres complexes. - Procédé selon la revendication 9, comprenant en outre :
- une étape de détermination d’un premier indicateur de diffusion multiple ε(r) qui est calculé par :
α(r) étant le taux de symétrie défini au point central (PC) du milieu. - Procédé selon la revendication 9, comprenant en outre :
- une étape de détermination d’une intensité afocale Ioff(r) qui est la moyenne d’un module au carré des réponses REP(r, Δ r), la moyenne étant effectuée pour des valeurs d’éloignement de module supérieures à une résolution prédéterminée wd(r) et/ou étant effectuée pour une plage de valeurs d’angle β ou pour une valeur d’angle prédéterminé βd, par exemple calculé par la formule suivante :
- une étape de détermination d’une intensité de diffusion multiple IM(r) qui est le produit entre le taux de symétrie α(r) et l’intensité afocale Ioff(r), c’est-à-dire :
- une étape de détermination d’une intensité de bruit IN(r) étant le produit entre un moins le taux de symétrie α(r) et l’intensité afocale Ioff(r), c’est-à-dire :
de telle sorte que l’on ait la relation suivante :
- Procédé selon la revendication 12, dans lequel le premier indicateur de diffusion multiple ε(r) est calculé par un ratio entre l’intensité de diffusion multiple IM(r) et l’intensité de bruit IN(r),c’est-à-dire :
- Procédé selon la revendication 12, comprenant en outre :
- une étape de détermination d’une intensité confocale Ion(r) qui est la valeur d’un module au carré de la réponse REP(r,|Δ r| = 0) pour une valeur de coordonnée d’éloignement nulle (|Δ r |= 0), c’est-à-dire pour un point du milieu pour lequel le premier point (P1), le deuxième point (P2) et le point central (PC) sont confondus,
- une étape de détermination d’une intensité de diffusion simple IS(r) calculé à partir de l’équation suivante :
- Procédé selon les revendications 12 et la revendication 14, comprenant en outre :
- une étape de détermination d’un deuxième indicateur de diffusion multiple γ(r) qui est calculé par :
- Procédé selon l’une quelconque des revendications 1 à 15, comprenant en outre :
- une étape de détermination d’une image d’un paramètre de caractérisation locale du milieu, ledit paramètre de caractérisation locale étant déterminé à partir de la réponse REP(r, Δ r). - Procédé selon la revendication 16, dans lequel le paramètre de caractérisation locale est choisi dans une liste comprenant la résolution w(r), le critère de focalisation F(r), le taux de symétrie α(r), le premier indicateur de diffusion multiple ε(r), le deuxième indicateur de diffusion multiple γ(r), l’intensité afocale Ioff(r), l’intensité confocale Ion(r), l’intensité de diffusion multiple IM(r), l’intensité de diffusion simple IS(r), l’intensité de bruit IN(r).
- Système (40) de caractérisation ultrason ore non invasive d’un milieu hétérogène (20)comprenant :
-un premier réseau (10) de transducteurs adaptés pour générer une série d’ondes ultrasonores incidentes dans une zone du milieu hétérogène, ladite série d’ondes ultrasonores incidentes constituant une base d’émission (i), et pour enregistrer en fonction du temps les ondes ultrasonores rétrodiffusées par ladite zone ; et
- une unité de calcul (42) couplée au premier réseau de transducteurs et adaptée pour :
- enregistrer d’une matrice de réflexion expérimentaleR ui (t) définie entre la base d’émission (i) en entrée et une base de réception (u) en sortie ;
- déterminer une réponse REP(r, Δ r) du milieu entre un transducteur virtuel d’entrée (TVin) de position spatialer in calculé à partir d’une focalisation en entrée de la matrice de réflexion expérimentale correspondant à une tache focale d’entrée autour d’un premier point (P1) et un transducteur virtuel de sortie (TVout) de position spatialer out calculé à partir d’une focalisation en sortie de la matrice de réflexion expérimentale correspondant à une tache focale de sortie autour d’une deuxième point (P2), ladite réponse étant exprimée en fonction d’un point central (PC) de position spatiale r dans le milieu, ledit point central (PC) étant situé au milieu des premier et deuxième points (P1, P2) et étant l’origine d’un axe de mesure (AXm) passant par lesdits premier et deuxième point, ledit axe de mesure formant un angle β par rapport à un premier axe (X) du milieu, et le premier point (P1) étant à une coordonnée d’éloignement+ Δ rsur l’axe de mesure, et le deuxième point (P2) étant à une coordonnée d’éloignement- Δ rsur l’axe de mesure.
Priority Applications (8)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1908904A FR3099580A1 (fr) | 2019-08-02 | 2019-08-02 | Procédé et système de caractérisation ultrasonore non invasive d’un milieu hétérogène |
KR1020227002139A KR20220038678A (ko) | 2019-08-02 | 2020-07-31 | 초음파를 사용하여 이종 매질을 비침습적으로 특성화하는 방법 및 시스템 |
PCT/FR2020/051416 WO2021023933A1 (fr) | 2019-08-02 | 2020-07-31 | Procédé et système de caractérisation ultrasonore non invasive d'un milieu hétérogène |
EP20820244.0A EP4007914A1 (fr) | 2019-08-02 | 2020-07-31 | Procédé et système de caractérisation ultrasonore non invasive d'un milieu hétérogène |
AU2020324570A AU2020324570A1 (en) | 2019-08-02 | 2020-07-31 | Method and system for non-invasively characterising a heterogeneous medium using ultrasound |
US17/631,929 US20240032889A1 (en) | 2019-08-02 | 2020-07-31 | Method and system for non-invasively characterising a heterogeneous medium using ultrasound |
CN202080052379.4A CN114144118B (zh) | 2019-08-02 | 2020-07-31 | 使用超声波非侵入性地表征异质介质的方法和系统 |
JP2022500603A JP7515563B2 (ja) | 2019-08-02 | 2020-07-31 | 超音波を使用して不均一な媒体を非侵襲的に特性化するための方法及びシステム |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1908904 | 2019-08-02 | ||
FR1908904A FR3099580A1 (fr) | 2019-08-02 | 2019-08-02 | Procédé et système de caractérisation ultrasonore non invasive d’un milieu hétérogène |
Publications (1)
Publication Number | Publication Date |
---|---|
FR3099580A1 true FR3099580A1 (fr) | 2021-02-05 |
Family
ID=68807030
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR1908904A Pending FR3099580A1 (fr) | 2019-08-02 | 2019-08-02 | Procédé et système de caractérisation ultrasonore non invasive d’un milieu hétérogène |
Country Status (8)
Country | Link |
---|---|
US (1) | US20240032889A1 (fr) |
EP (1) | EP4007914A1 (fr) |
JP (1) | JP7515563B2 (fr) |
KR (1) | KR20220038678A (fr) |
CN (1) | CN114144118B (fr) |
AU (1) | AU2020324570A1 (fr) |
FR (1) | FR3099580A1 (fr) |
WO (1) | WO2021023933A1 (fr) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4817614A (en) * | 1986-08-20 | 1989-04-04 | Siemens Aktiengesellschaft | Method and apparatus for adaptive focusing in a medical ultrasound imaging apparatus |
US5720289A (en) * | 1994-08-05 | 1998-02-24 | Acuson Corporation | Method and apparatus for a geometric aberration transform in an adaptive focusing ultrasound beamformer system |
WO2010001027A1 (fr) | 2008-06-09 | 2010-01-07 | Centre National De La Recherche Scientifique - Cnrs - | Procede et dispositif de sondage par propagation d'ondes. |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2815717B1 (fr) * | 2000-10-20 | 2003-01-10 | Centre Nat Rech Scient | Procede et dispositif non invasif de focalisation d'ondes acoustiques |
WO2003017843A1 (fr) * | 2001-08-24 | 2003-03-06 | The Brigham And Women's Hospital, Inc. | Therapie par ultrasons |
US9117439B2 (en) * | 2008-03-13 | 2015-08-25 | Supersonic Imagine | Method and apparatus for ultrasound synthetic imagining |
DE102009025463A1 (de) * | 2009-06-12 | 2011-03-10 | Technische Universität Dresden | Anordnung und Verfahren zur Bestimmung von Schichtdicken und Schallgeschwindigkeiten in Medien mit Hilfe von Ultraschall |
ITGE20090070A1 (it) * | 2009-08-31 | 2011-03-01 | Esaote Spa | Metodo e dispositivo per il rilevamento e la visualizzazione di informazioni emodinamiche in particolare del flusso ematico nelle vene, mediante ultrasoni |
JP6000196B2 (ja) | 2012-09-27 | 2016-09-28 | 富士フイルム株式会社 | 超音波診断装置、音速決定方法およびプログラム |
CN103969651A (zh) * | 2014-04-24 | 2014-08-06 | 中国科学院声学研究所 | 自适应声学成像方法 |
JP6724797B2 (ja) | 2017-01-13 | 2020-07-15 | コニカミノルタ株式会社 | 超音波診断装置、および、その制御方法 |
-
2019
- 2019-08-02 FR FR1908904A patent/FR3099580A1/fr active Pending
-
2020
- 2020-07-31 KR KR1020227002139A patent/KR20220038678A/ko active Search and Examination
- 2020-07-31 CN CN202080052379.4A patent/CN114144118B/zh active Active
- 2020-07-31 JP JP2022500603A patent/JP7515563B2/ja active Active
- 2020-07-31 AU AU2020324570A patent/AU2020324570A1/en active Pending
- 2020-07-31 US US17/631,929 patent/US20240032889A1/en active Pending
- 2020-07-31 EP EP20820244.0A patent/EP4007914A1/fr active Pending
- 2020-07-31 WO PCT/FR2020/051416 patent/WO2021023933A1/fr active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4817614A (en) * | 1986-08-20 | 1989-04-04 | Siemens Aktiengesellschaft | Method and apparatus for adaptive focusing in a medical ultrasound imaging apparatus |
US5720289A (en) * | 1994-08-05 | 1998-02-24 | Acuson Corporation | Method and apparatus for a geometric aberration transform in an adaptive focusing ultrasound beamformer system |
WO2010001027A1 (fr) | 2008-06-09 | 2010-01-07 | Centre National De La Recherche Scientifique - Cnrs - | Procede et dispositif de sondage par propagation d'ondes. |
Non-Patent Citations (2)
Title |
---|
G. MONTALDO ET AL.: "Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography", IEEE TRANS. ULTRASON., FERROELECT. FREQ. CONTROL, vol. 56, 2009, pages 489 - 506, XP011255897 |
RAOUL MALLARTMATHIAS FINK: "The van Cittert-Zernike theorem in pulse echo measurements", J. ACOUST. SOC. AM., vol. 90, no. 5, November 1991 (1991-11-01), XP000236651, DOI: 10.1121/1.401867 |
Also Published As
Publication number | Publication date |
---|---|
CN114144118B (zh) | 2024-07-12 |
KR20220038678A (ko) | 2022-03-29 |
CN114144118A (zh) | 2022-03-04 |
US20240032889A1 (en) | 2024-02-01 |
WO2021023933A1 (fr) | 2021-02-11 |
EP4007914A1 (fr) | 2022-06-08 |
AU2020324570A1 (en) | 2022-01-27 |
JP2022542005A (ja) | 2022-09-29 |
JP7515563B2 (ja) | 2024-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3824280B1 (fr) | Procédés et systèmes de caractérisation ultrasonore non invasive d'un milieu hétérogène | |
EP3967239B1 (fr) | Procédé et système de caractérisation ultrasonore d'un milieu | |
EP2440140B1 (fr) | Procede et dispositif ultrasonores pour caracteriser un milieu | |
EP3967241A1 (fr) | Procédé et système de caractérisation ultrasonore d'un milieu | |
EP3967240A1 (fr) | Procédé et système de caractérisation ultrasonore d'un milieu | |
WO2020128344A1 (fr) | Procédé d'imagerie ultrasonore par transformée de fourier bidimensionnelle, programme d'ordinateur et dispositif de sondage à ultrasons correspondants | |
EP2342582B1 (fr) | Procede et dispositif de sondage par propagation d'ondes | |
WO2022058675A1 (fr) | Procédé et système de caractérisation ultrasonore d'un milieu | |
EP3968015A1 (fr) | Procédé et système de caractérisation ultrasonore d'un milieu | |
FR3099580A1 (fr) | Procédé et système de caractérisation ultrasonore non invasive d’un milieu hétérogène | |
CN112823283B (zh) | 通过使用超声非侵入性地表征非均匀介质的方法和系统 | |
WO2023170299A1 (fr) | Procede et dispositif d'estimation locale de la vitesse du son dans une region d'interet d'un milieu | |
WO2018109314A1 (fr) | Procede de traitement de signaux issus d'une acquisition par sondage ultrasonore, programme d'ordinateur et dispositif de sondage a ultrasons correspondants |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 2 |
|
PLSC | Publication of the preliminary search report |
Effective date: 20210205 |
|
PLFP | Fee payment |
Year of fee payment: 3 |
|
PLFP | Fee payment |
Year of fee payment: 4 |
|
PLFP | Fee payment |
Year of fee payment: 5 |