Brevet FR3163830A1 — publié le 20260102

Traitement d’élastographie ultrasonore pour caractériser un milieu d’intérêt La divulgation concerne un procédé de traitement d’élastographie ultrasonore et un dispositif (20) correspondant. Le procédé comprend : obtention de premières données (DT1) représentatives d’une réponse ultrasonore (16), d’un milieu d’intérêt (4), à des impulsions ultrasonores (8) focalisées dans le milieu et à des ondes de sondage (14) émises séquentiellement selon des valeurs (VL) distinctes d’au moins un paramètre (PR) des ondes de sondage ; génération, à partir des premières données respectivement pour chaque dite valeur (VL) de paramètre, de deuxièmes données (DT2) représentatives d’une propagation d’ondes de cisaillement dans le milieu ; estimation, à partir des deuxièmes données associées à chaque dite valeur (VL) de paramètre, d’une célérité individuelle (VI) représentative de la célérité de déplacement des ondes de cisaillement pour respectivement ladite valeur ; et estimation d’au moins une célérité (VC), dite célérité composée, des ondes de cisaillement par une combinaison des célérités individuelles. Figure pour l’abrégé : Figure 6 Traitement d’élastographie ultrasonore pour caractériser un milieu d’intérêt La présente divulgation concerne les examens ou scans, ou observations, utilisant des ondes de type ultrasonore pour caractériser une région d’intérêt et porte plus particulièrement sur le traitement de données caractérisant une réponse d’un tel milieu d’intérêt à un scan ultrasonore en mode d’élastographie. Les dispositifs ultrasonores trouvent aujourd’hui de nombreuses applications, dans les domaines notamment de l'acoustique, de l’étude des matériaux, de l'imagerie médicale et vétérinaire, ainsi que dans le domaine de la biomédecine. De tels dispositifs sont en particulier utilisés couramment pour caractériser et observer des structures internes du corps humain par exemple. Le principe d’un scan ultrasonore repose ainsi sur l'émission, par une sonde, d’ondes ultrasonores qui traversent une région d’intérêt, par exemple des tissus d’un corps examiné, ces ondes se réfléchissant lorsqu'elles rencontrent des éléments particuliers, tels que des interfaces entre différents types de tissus. Les ondes réfléchies sont ensuite captées en retour par la sonde ou un autre dispositif et converties en signaux électriques. Des moyens de traitement analysent alors la réponse ultrasonore du milieu d’intérêt grâce aux signaux pour en générer des données caractérisant le milieu. Ainsi, à partir d’une série d’insonifications du milieu au moyen d’éléments transducteurs ultrasonores, il est possible de recevoir et analyser les échos rétrodiffusés par ce même milieu et d’en déduire des caractéristiques du milieu étudié. En imagerie par ultrasons par exemple, dit aussi échographie, le traitement de ces échos ultrasonores permet de générer diverses informations relatives au milieu étudié, telles qu’une image ou série d’images selon un mode approprié, tel qu’un mode B (c'est-à-dire un mode bidimensionnel d’intensité du signal reçu), un mode Doppler ou un mode d'élastographie par ondes de cisaillement appelé aussi SWE pour « ShearWave® Elastography ». L'élastographie ultrasonore est une technique de caractérisation ultrasonore qui vise à estimer l’élasticité ou la rigidité d’un milieu d’intérêt examiné. Elle fonctionne en focalisant des impulsions ultrasonores dans une zone spécifique d’une région d’intérêt, du corps d’un sujet par exemple, ce qui génère des ondes de cisaillement se propageant latéralement dans ladite région, par exemple à travers les tissus du corps examiné. Pour caractériser ces ondes de cisaillement, des ondes de sondage sont en outre émises, permettant de suivre la propagation de ces ondes de cisaillement en temps réel ou quasi-temps réel. En analysant la vitesse de propagation de ces ondes de cisaillement, il est en effet possible d'estimer la rigidité du milieu d'intérêt. Les tissus ou zones plus rigides permettent aux ondes de cisaillement de se propager plus rapidement, alors que les tissus (ou zones) plus mous (ou molles) ralentissent leur progression ou proposent une progression moins rapide desdites ondes. A partir de cette technique, on peut créer une carte de l'élasticité des tissus du milieu examiné, utile notamment pour détecter certaines conditions médicales, par exemple des anomalies telles que les tumeurs ou les fibroses, qui présentent souvent une rigidité différente de celle des tissus sains. Cette technique permet également d’indiquer un niveau de dureté des tissus, que l’on peut rapprocher en médecine du palper du praticien. Actuellement, les systèmes ultrasonores fonctionnant en mode d'élastographie ultrasonore présentent des performances qui ne sont pas entièrement satisfaisantes en termes notamment de précision et de fiabilité de l’estimation de la vitesse d’ondes de cisaillement se propageant dans un milieu d’intérêt, ce qui limite la capacité de ces systèmes à caractériser le milieu considéré. Ces limitations peuvent entraîner des erreurs ou imprécisions dans l'évaluation de la rigidité tissulaire, voire dégrader la qualité des images dans le cas d’application dans le domaine de l’imagerie, notamment en imagerie médicale. Une mauvaise estimation de la vitesse des ondes de cisaillement peut entraîner une interprétation inexacte de la rigidité des tissus ou de tout ou partie des zones examinées, ce qui limite notamment l'utilisation de l'élastographie ultrasonore dans des applications médicales. Exposé de la divulgation L’un des objets de la présente divulgation est de résoudre au moins l’un des problèmes ou déficiences décrits précédemment. En particulier, un objet de la présente divulgation est de traiter des données d’élastographie ultrasonore de façon plus efficace pour permettre une estimation précise et fiable de la célérité d’ondes de cisaillement se propageant dans un milieu d’intérêt. En particulier, un objet de la présente divulgation est de traiter des données ultrasonores représentatives d’une réponse ultrasonore, d’un milieu d’intérêt, à des impulsions ultrasonores focalisées dans ledit milieu en combinaison avec des ondes de sondage dont on fait varier au moins un paramètre au cours du temps. A cet effet, selon un premier aspect, la présente divulgation vise un procédé de traitement d’élastographie ultrasonore, dit aussi procédé de traitement, le procédé comprenant : a) obtention de premières données représentatives d’une réponse ultrasonore, d’un milieu d’intérêt, à des impulsions ultrasonores focalisées dans ledit milieu en combinaison avec des ondes de sondage émises dans ledit milieu séquentiellement selon des valeurs distinctes d’au moins un paramètre desdites ondes de sondage ; b) génération, à partir des premières données respectivement pour chaque valeur dudit au moins un paramètre, de deuxièmes données représentatives d’une propagation d’ondes de cisaillement dans le milieu d’intérêt sous l’effet des impulsions ultrasonores ; c) estimation, à partir des deuxièmes données associées à chaque valeur dudit au moins un paramètre, d’une célérité individuelle représentative de la célérité de déplacement des ondes de cisaillement pour respectivement ladite valeur ; et d) estimation d’au moins une célérité, dite célérité composée, des ondes de cisaillement par une combinaison des célérités individuelles. Le procédé selon la divulgation peut comporter d’autres caractéristiques qui peuvent être prises séparément ou en combinaison, notamment parmi les modes de réalisation qui suivent qui sont présentés à titre d’illustration uniquement et peuvent être combinés ou associés sauf stipulation contraire. Selon un exemple, ledit au moins un paramètre des ondes de sondage comprend au moins l’un parmi : - un angle d’émission selon lequel les ondes de sondage sont émises en direction du milieu d’intérêt ; et - une fréquence selon laquelle sont émises les ondes de sondage en direction du milieu d’intérêt. Selon un exemple, les ondes de sondage comprennent des ondes planes ultrasonores utilisées pour sonder le milieu d’intérêt. Selon un exemple, le procédé comprend, préalablement à l’obtention a) : e1) formation des ondes de cisaillement dans le milieu d’intérêt par focalisation des impulsions ultrasonores ; et e2) émission, dans le milieu d’intérêt, des ondes de sondage configurées séquentiellement selon les valeurs distinctes dudit au moins un paramètre pour détecter la réponse ultrasonore du milieu d’intérêt pour chaque dite valeur. Selon un exemple, les deuxièmes données générées en b) définissent, pour chaque valeur dudit au moins un paramètre, une séquence d’images représentative de la propagation des ondes de cisaillement causées par les impulsions ultrasonores dans le milieu d’intérêt. Selon un exemple, au cours de l’estimation d), la combinaison des célérités individuelles est déterminée par un calcul d’une moyenne des célérités individuelles pour estimer ladite au moins une célérité composée. Selon un exemple, l’estimation d) comprend : - attribution de poids, en association avec les célérités individuelles, en fonction d’au moins une caractéristique des deuxièmes données générées respectivement pour chaque valeur dudit au moins un paramètre ; et - estimation de ladite au moins une célérité composée par calcul de la moyenne en pondérant les célérités individuelles selon lesdits poids associés. Selon un exemple, l’estimation d) comprend : - détermination de variances spatiales caractérisant respectivement les deuxièmes données associées à chaque valeur dudit au moins un paramètre ; et - détermination, à partir des variances spatiales, de scores de confiance associés respectivement aux célérités individuelles, le poids attribué à chaque célérité individuelle étant fonction du score de confiance associé. Selon un exemple, l’estimation d) comprend : - si un score de confiance associée à des deuxièmes données correspondant à une valeur dudit au moins un paramètre ne satisfait pas à un au moins critère de qualité, le poids attribué à la célérité individuelle associée à ladite valeur est nul de sorte à exclure ladite célérité individuelle de l’estimation d) de ladite au moins une célérité composée. Selon un exemple, le procédé comprend en outre : f) estimation, à partir de ladite au moins une célérité composée, d’une dureté du milieu d’intérêt. Selon un exemple, au cours de l’estimation d), une pluralité de célérités composées sont estimées en différentes positions du milieu d’intérêt, l’estimation d) comprenant en outre : - génération, à partir des célérités composées, de données cartographiques représentatives des célérités composées, et/ou de duretés estimées à partir desdites célérités composées, dans respectivement les différentes positions du milieu d’intérêt. Selon un exemple, le procédé comprend en outre : - affichage des célérités composées, et/ou de duretés estimées à partir desdites célérités composées, sous forme d’au moins une image ou carte. Selon un exemple, le procédé comprend, après une phase préliminaire au cours de laquelle sont exécutées les étapes a) à d) pour chaque dite valeur distinctes dudit au moins un paramètre, au moins une itération périodique d’un cycle comprenant : - estimation d’au moins une célérité individuelle additionnelle par réitération des étapes a) à c) pour au moins l’une parmi les dites valeurs dudit au moins un paramètre ; et - estimation d’au moins une célérité composée des ondes de cisaillement par réitération de l’étape d) en utilisant la célérité individuelle additionnelle et au moins une célérité individuelle estimée au cours d’une itération précédente du cycle et/ou au cours de la phase préliminaire. Selon un deuxième aspect, la présente divulgation peut impliquer un programme d'ordinateur comprenant des instructions qui, lorsque le programme est exécuté par un ordinateur, cause la mise en œuvre du procédé selon le premier aspect. En particulier, les différentes étapes du procédé selon le premier aspect peuvent être définies par des instructions de programmes d’ordinateurs. Un tel programme d’ordinateur peut utiliser n’importe quel langage de programmation ou équivalent, et il peut se trouver sous la forme d’un code source, d’un code objet, ou d’un code intermédiaire entre un code source et un code objet, tel que dans une forme partiellement compilée, ou sous n’importe quelle autre forme souhaitable. Selon un troisième aspect, la présente divulgation concerne un support d’enregistrement (ou support d’informations), lisible par un ordinateur (ou un processeur), sur lequel est enregistré un programme d’ordinateur selon le ce même aspect de la présente divulgation. D’une part, le support d’enregistrement peut être n'importe quel entité ou dispositif capable de stocker le programme, telle qu’au moins une mémoire volatile et/ou non volatile. Par exemple, le support peut comporter un moyen de stockage, tel qu'une mémoire non-volatile réinscriptible, une mémoire ROM, un CD-ROM ou une mémoire ROM de type circuit microélectronique, ou encore un moyen d'enregistrement magnétique ou un disque dur. Cette mémoire peut par exemple comprendre une mémoire de carte graphique (ou carte vidéo), ce type de mémoire étant notamment apte à traiter des données d’image (ou données vidéo). D'autre part, ce support d’enregistrement peut également être un moyen basé sur des signaux transmissible tel qu'un signal électrique ou optique, un tel signal pouvant être acheminé via un câble électrique ou optique, par radio classique ou hertzienne ou par faisceau laser autodirigé ou par d'autres moyens. Le programme d’ordinateur selon la présente divulgation peut être en particulier téléchargé grâce à un réseau filaire ou non, de type local ou non (Bluetooth® par exemple, Wi-Fi, Ethernet, Internet, 4G, 5G ou autres). Alternativement, le support d'enregistrement peut être un circuit intégré dans lequel le programme d’ordinateur est incorporé, le circuit intégré étant adapté pour exécuter ou pour être utilisé dans l'exécution du procédé en question. Selon un quatrième aspect, la présente divulgation concerne un dispositif de traitement d’élastographie ultrasonore, dit aussi dispositif de traitement, configuré pour mettre en œuvre le procédé du premier aspect de la présente divulgation. Selon un exemple, le dispositif de traitement comprend une mémoire associée à un processeur, cette mémoire comprenant un programme d’ordinateur selon la présente divulgation. Selon un exemple, la présente divulgation vise un dispositif de traitement d’élastographie ultrasonore, ledit dispositif comprenant : - un module d’obtention configuré pour obtenir des premières données représentatives d’une réponse ultrasonore, d’un milieu d’intérêt, à des impulsions ultrasonores focalisées dans ledit milieu en combinaison avec des ondes de sondage émises dans ledit milieu séquentiellement selon des valeurs distinctes d’au moins un paramètre desdites ondes de sondage ; - un module de génération configuré pour générer, à partir des premières données respectivement pour chaque valeur dudit au moins un paramètre, de deuxièmes données représentatives d’une propagation d’ondes de cisaillement dans le milieu d’intérêt sous l’effet des impulsions ultrasonores ; - un premier module d’estimation configuré pour estimer, à partir des deuxièmes données associées à chaque valeur dudit au moins un paramètre, une célérité individuelle représentative de la célérité de déplacement des ondes de cisaillement pour respectivement ladite valeur ; et - un deuxième module d’estimation configuré pour estimer au moins une célérité, dite célérité composée, des ondes de cisaillement par une combinaison des célérités individuelles. A noter que les différents modes de réalisation mentionnés ci-avant (ainsi que ceux décrits ci-après) en relation avec le procédé de traitement de l’invention ainsi que les avantages associés s’appliquent de façon analogue au dispositif de traitement de l’invention. Pour chaque étape du procédé de traitement, le dispositif de traitement de l’invention peut comprendre un module correspondant configuré pour réaliser ladite étape. Selon un mode de réalisation, l'invention est mise en œuvre au moyen de composants logiciels et/ou matériels. Dans cette optique, le terme « module » peut correspondre dans ce document aussi bien à un composant logiciel, qu'à un composant matériel ou à un ensemble de composants matériels et/ou logiciels. Selon un cinquième aspect, la présente divulgation concerne un système de traitement comprenant : - le dispositif de traitement selon le quatrième aspect de la divulgation ; et - au moins une sonde ultrasonore configurée pour focaliser des impulsions ultrasonores dans ledit milieu et/ou émettre les ondes de sondage ultrasonores dans ledit milieu. En mettant en œuvre la présente divulgation, il devient avantageusement possible de traiter des données d’élastographie ultrasonore de façon efficace pour permettre une estimation précise et fiable de la célérité d’ondes de cisaillement se propageant dans un milieu d’intérêt. En particulier, la divulgation prévoit le traitement de premières données représentatives d’une réponse ultrasonore, d’un milieu d’intérêt, à des impulsions ultrasonores focalisées dans le milieu en combinaison avec des ondes de sondage dont on fait varier au moins un paramètre au cours du temps. En générant des célérités individuelles associées respectivement à différentes valeurs d’un paramètre des ondes de sondage que l’on fait varier, on peut avantageusement obtenir de multiples estimations de célérités des ondes de cisaillement pour une position donnée du milieu 4. En combinant ces célérités individuelles obtenues pour chaque valeur du paramètre, par exemple par le calcul d’une moyenne ou tout autre combinaison approprié, on peut ainsi obtenir une estimation de qualité de la célérité des ondes de cisaillement se propageant dans un milieu. Les caractéristiques et avantages de la présente divulgation apparaitront plus précisément à la lecture de la description qui va suivre, donnée uniquement à titre d'exemples non limitatifs, et faite en référence aux figures annexées. Notamment, les exemples illustrés dans les figures peuvent être combinés entre eux, sauf incohérence notoire. Brève description des figures D’autres caractéristiques et avantages de la présente divulgation ressortiront de la description des exemples de réalisation non limitatifs de la présente divulgation exposés ci-après, en référence aux figures 1 à 16 annexées, sur lesquelles : FIG. 1 représente schématiquement un système ultrasonore apte à réaliser un scan ultrasonore selon au moins un exemple de réalisation de la présente divulgation ; FIG. 2 représente schématiquement l’émission d’ondes de compression dans une région d’intérêt pour provoquer des ondes de cisaillement, selon au moins un exemple de réalisation de la présente divulgation ; FIG. 3 représente schématiquement l’émission d’ondes de sondage dans une région d’intérêt pour analyser la propagation d’ondes de cisaillement dans ladite région, selon au moins un exemple de réalisation de la présente divulgation ; FIG. 4 représente schématiquement la variation d’un paramètre (angle d’émission) des ondes de sondage émises dans une région d’intérêt, selon au moins un exemple de réalisation de la présente divulgation ; FIG. 5 représente schématiquement la variation d’un paramètre (fréquence) des ondes de sondage émises dans une région d’intérêt, selon au moins un exemple de réalisation de la présente divulgation ; FIG. 6 représente schématiquement un système de traitement, et un dispositif de traitement, selon au moins un exemple de réalisation de la présente divulgation ; FIG. 7 représente schématiquement un système de traitement, et un dispositif de traitement, selon au moins un exemple de réalisation de la présente divulgation ; FIG. 8 représente schématiquement un système de traitement, et un dispositif de traitement, selon au moins un exemple de réalisation de la présente divulgation ; FIG. 9 représente schématiquement, sous forme d’un diagramme, les étapes d’un procédé de traitement mis en œuvre par un dispositif de traitement, selon au moins un exemple de réalisation de la présente divulgation ; FIG. 10 représente schématiquement la mise en œuvre du procédé de traitement de la FIG. 9 , selon au moins un exemple de réalisation de la présente divulgation ; FIG. 11 représente schématiquement le procédé de traitement de la FIG. 9 selon au moins un exemple réalisation de la présente divulgation ; et FIG. 12 représente schématiquement des données générées au cours du procédé de traitement de la FIG. 9 , selon au moins un exemple de réalisation de la présente divulgation ; FIG. 13 , FIG. 14 et FIG. 15 représentent schématiquement la mise en œuvre du procédé de traitement de la FIG. 9 , selon au moins un exemple de réalisation de la présente divulgation ; et FIG. 16 représente schématiquement le procédé de traitement de la FIG. 9 selon au moins un exemple de réalisation de la présente divulgation. La présente divulgation porte sur le traitement de données d’élastographie ultrasonore, c’est-à-dire des données obtenues par une étude (un scan) ultrasonore d’une région d’intérêt en mode d’élastographie ultrasonore, dit aussi mode d’élastographie par ondes de cisaillement ou SWE pour « ShearWave® Elastography ». Un scan au sens de la présente divulgation correspond à un scan ultrasonore (ou scan échographique), une observation ultrasonore (ou échographique) ou encore une phase d’observation par ultrasons (ou échographie), c’est-à-dire une phase d’émission et réception d’ondes ultrasonores en direction du milieu observé, et respectivement en provenance dudit milieu. Un tel scan comprend la propagation d’ondes ultrasonores dans une région d’intérêt et la récupération en retour d’une réponse ultrasonore, c’est-à-dire d’échos échographiques (ou échos ultrasonores) émis par la région d’intérêt. Un scan permet d’obtenir (ou générer) des données ultrasonores représentatives de la région d’intérêt. Ces données peuvent alors être utilisées pour par exemple caractériser le milieu, construire des images ou établir des cartographies représentatives du milieu étudié, bien que d’autres utilisations des données ultrasonores soient possibles. En mode d’élastographie par ondes de cisaillement, un scan ultrasonore permet d’obtenir des données ultrasonores qui, une fois traitée, renseignent sur l’élasticité, la rigidité de la région d’intérêt observée. Le principe de ce mode ultrasonore est rappelé ci-après en référence aux figures. Comme illustré en FIG. 1 selon un exemple, il est possible d’utiliser un système ultrasonore SY1 comprenant un dispositif émetteur 2 et un dispositif de contrôle 5 pour réaliser un scan ultrasonore d’une région (ou milieu) d’intérêt 4, par exemple dans une partie d’un corps d’un sujet (un corps humain ou animal par exemple). Pour ce faire, la sonde ultrasonore 2 est pilotée par le dispositif 5 pour émettre des ondes ultrasonores W1 dans la région d’intérêt 4. Dans l’exemple particulier de l’imagerie ultrasonore, les échos échographiques W2 reçus en retour depuis la région d’intérêt 4 sont utilisés en particulier pour caractériser la région d’intérêt 4. Une image ultrasonore représentative de cette région peut en particulier être générée à partir des échos ultrasonores W2 reçus. Le dispositif émetteur ultrasonore 20 est par exemple une sonde ultrasonore (ou sonde échographique) configurée pour émettre des ondes ultrasonores W1 vers la région d’intérêt 4 en réponse à un ou des signaux électriques SG1 fournis par le dispositif de contrôle 5. Pour ce faire, le dispositif de contrôle 5 peut en particulier comprendre une unité de commande 6 comportant un ou des pulseurs électroniques (dits plus simplement « pulseurs ») ou des amplificateurs linéaires, configurés pour générer le ou les signaux électriques SG1 destinés à piloter la sonde ultrasonore 2. La sonde ultrasonore 2 peut comprendre un ou des transducteurs ultrasonores 2a ( FIG. 2 ), chacun d’eux étant piloté par un signal électrique SG1 délivré par le dispositif de contrôle 5. En réponse aux signaux électriques SG1 reçus, les transducteurs 2a produisent des ondes ultrasonores W1 en direction de la région d’intérêt 4. La transmission des signaux électriques SG1 peut être assurée par une liaison depuis les pulseurs vers les transducteurs 2a de la sonde 2. La sonde ultrasonore 2 peut donc être configurée pour émettre, au cours d’une opération d’émission, des ondes ultrasonores W1 vers la région d’intérêt 4 et pour recevoir en retour, au cours d’une opération de réception, des ondes ultrasonores W2 depuis la région d’intérêt 4, par exemple sous forme d’échos ultrasonores causés par les ondes W1. Ainsi, chaque transducteur 2a peut être configuré pour convertir un signal électrique SG1 fourni par le dispositif de contrôle 5 en ondes ultrasonores W1 et réciproquement. A noter toutefois que des sondes 2 différentes peuvent être utilisées pour les opérations respectivement d’émission et de réception. Les éléments transducteurs 2a peuvent être disposés selon diverses configurations selon le cas, par exemple selon une ligne de transducteurs, une matrice, ou selon un réseau de transducteurs ou toute autre configuration appropriée. Comme illustré dans les figures 2 et 3 selon un exemple, le système SY1 comprenant la sonde 2 et le dispositif de contrôle 5 peut en particulier être configuré pour opérer selon un mode d’élastographie par ondes de cisaillement, dit aussi mode ShearWave ® ou « ShearWave ® Elastography ». De façon générale, le mode d’élastographie ultrasonore comprend une phase d’excitation de la région d’intérêt par émission d’ondes ultrasonores et une phase d’acquisition et traitement au cours de laquelle des données ultrasonores émises en retour par la région d’intérêt sont acquises et traitées. Plus précisément, en élastographie impulsionnelle, la sonde 2 est utilisée pour émettre des impulsions ultrasonores 8 focalisées dans la région d’intérêt 4 ( FIG. 2 ). Ces impulsions ultrasonores 8, typiquement émises sous forme de tirs d’impulsions, forment des ondes de compression causant des contraintes mécaniques dans la région d’intérêt. Dans la région du point focal 10, se produit un phénomène dit de « pression de radiation acoustique » selon lequel un transfert d’énergie causées par les impulsions 8 engendre des déplacements de matière dans le milieu étudié, ce qui génère des ondes de cisaillement 12 se propageant latéralement autour du point focal 10. Grâce à une analyse de ces ondes de cisaillement 12, on peut estimer la vitesse (ou célérité) de ces ondes et en déduire la rigidité et/ou élasticité du milieu. La détection et l’analyse du déplacement des ondes de cisaillement 12 sont réalisées grâce à l’émission d’ondes de sondage 14 dans la région d’intérêt 4 ( FIG. 3 ). Pour ce faire, la sonde ultrasonore 2 peut être configurée pour émettre alternativement des impulsions ultrasonores 8 et des ondes de sondage 14 dans la région d’intérêt 4. Les impulsions ultrasonores 8 excitent la région d’intérêt 4 pour former les ondes de cisaillement 12 tandis que les ondes de sondage 14 sont émises pour obtenir une réponse ultrasonore 16 (sous forme d’échos rétrodiffusés) qui permet de caractériser ces ondes de cisaillement 12. On supposera par la suite que les ondes de sondages 14 utilisées sont des ondes planes, à titre d’exemple, ce qui permet avantageusement d’obtenir divers angles d’émission et donc de couvrir une zone appropriée du milieu. A noter toutefois que le type et les caractéristiques des ondes de sondage 14 peuvent être adaptés selon le besoin et le milieu étudié en particulier. Des exemples de réalisation faisant intervenir, en tant qu’ondes de sondage 14, des ondes ultrasonores divergentes ou des ondes ultrasonores focalisées sont possibles. Les signaux de la réponse ultrasonore 16 peuvent être acquis par la sonde 2 et transmis sous forme de données DT1, dites premières données, à un dispositif de traitement 20, tel que représenté en FIG. 3 selon un exemple de réalisation. Plus le milieu examiné présente une dureté importante et plus la célérité des ondes de cisaillement 12 est élevée. Par l’analyse des premières données DT1, on peut ainsi évaluer la vitesse (ou célérité) de propagation des ondes de cisaillement 12 et en déduire une estimation de la dureté et/ou élasticité de la région d’intérêt 4. En particulier, la célérité, notée c, des ondes de cisaillement 12 est directement liée au module µ de cisaillement par la relation suivante : où ρ est la masse volumique de la région d’intérêt 4 examinée. A titre d’exemple, la masse volumique ρ dans le corps humain est environ égale à 1000 kg/m 3 . En outre, le module µ de cisaillement est directement proportionnel au module de Young E selon l’expression suivante : [Math 1] (equation 1) A partir de la célérité c des ondes de cisaillement 12, on peut donc en déduire le module de Young E, et donc la dureté c, où : [Math 2] (equation 2) Comme décrit par la suite, les célérités des ondes de cisaillement, ainsi que les niveaux de rigidité ou élasticité obtenus, peuvent être exploités de diverses manières, par exemple pour la génération et l’affichage de cartes paramétriques permettant de visualiser des niveaux de dureté dans le milieu d’intérêt 4, ces cartes pouvant le cas échéant être superposées ou associées à une ou plusieurs images ultrasonores relatives au mode B (présentée en nuances de gris). Comme illustré en FIG. 4 , il est possible de faire varier au cours du temps l’angle d’émission θ selon lequel sont émises les ondes de sondage 14, à savoir des ondes planes dans les exemples considérés. En faisant varier cet angle θ, on peut avantageusement améliorer la qualité de l’estimation de la célérité c des ondes de cisaillement 14 dans le milieu considéré. Selon un exemple, on peut ainsi configurer la sonde ultrasonore 2 de sorte à ce qu’elle émette séquentiellement des ondes de sondage 14 selon des valeurs distinctes VL de l’angle d’émission θ, par exemple selon trois valeurs distinctes notées VL1, VL2 et VL3. En moyennant la réponse ultrasonore 16 (figures 3-4) collectée par la sonde 2 selon les différentes valeurs VL1-VL3 d’angle θ, on peut avantageusement obtenir une estimation de qualité exploitable de la célérité des ondes de cisaillement 14, avec un niveau de bruit réduit et une meilleure résolution, par comparaison à une solution basée sur l’analyse de la réponse ultrasonore 16 selon un unique angle d’émission θ (une seule valeur d’angle). Cette technique, dite de moyennage des données ultrasonores, permet de déterminer une estimation de confiance statistique à différentes positions de la région d’intérêt 4. A noter cependant que la technique exposée ci-dessus de moyennage des données ultrasonores ne permet de déterminer qu’une estimation unique de la célérité c des ondes de cisaillement 14 en chaque point étudié, pour un scan ultrasonore donné. Il a été observé que les résultats ne sont pas toujours satisfaisants en termes notamment de précision et de fiabilité, ce qui peut entraver la caractérisation du milieu observé. Lorsque les données ultrasonores moyennées ne permettent pas d’obtenir une estimation de la célérité de suffisamment bonne qualité, un résultat pertinent ne peut être obtenu, ce qui limite les performances de cette technique. Dans le cas par exemple d’applications en imagerie médicale, il a été observé que les données ultrasonores moyennées selon les différentes valeurs d’angle VL ne permettent pas toujours d’obtenir une estimation de célérité de qualité suffisante, conduisant dans certains cas à ignorer le résultat pour certains points dans l’espace de la région d’intérêt 4, ce qui peut se traduire par exemple par une zone dans la région étudiée dans laquelle les informations sont non exploitables ou sans pouvoir obtenir d’information utile dans la ou les images cartographiques produites. La présente divulgation se propose de répondre notamment aux problèmes et contraintes précédemment décrits grâce notamment à un procédé de traitement d’élastographie ultrasonore et à un dispositif de traitement correspondant, visant à caractériser une région d’intérêt à partir de données ultrasonores obtenues selon un mode d’élastographie par ondes de cisaillement. Selon des modes particuliers de réalisation, le procédé de traitement comprend ainsi : a) obtention de premières données représentatives d’une réponse ultrasonore, d’un milieu d’intérêt, à des impulsions ultrasonores focalisées dans ledit milieu en combinaison avec des ondes de sondage émises dans ledit milieu séquentiellement selon des valeurs distinctes d’au moins un paramètre desdites ondes de sondage ; b) génération, à partir des premières données respectivement pour chaque valeur dudit au moins un paramètre, de deuxièmes données représentatives d’une propagation d’ondes de cisaillement dans le milieu d’intérêt sous l’effet des impulsions ultrasonores ; c) estimation, à partir des deuxièmes données associées à chaque valeur dudit au moins un paramètre, d’une célérité individuelle représentatives de la célérité de déplacement des ondes de cisaillement pour respectivement ladite valeur ; et d) estimation d’au moins une célérité, dite célérité composée, des ondes de cisaillement par une combinaison des célérités individuelles. D’autres aspects et avantages de la présente divulgation ressortiront des exemples de réalisation décrits ci-dessous en référence aux figures mentionnées ci-avant. Des exemples de réalisation de la présente divulgation sont décrits ci-après conjointement aux figures 1 à 16, dans lesquels le dispositif de traitement d’élastographie ultrasonore (dit aussi dispositif de traitement) 20 représenté notamment en FIG. 3 est utilisé pour mettre en œuvre le procédé de traitement d’élastographie ultrasonore (dit aussi procédé de traitement) de la présente divulgation. Sauf indication contraire ou incohérence notoire, les éléments communs ou analogues à plusieurs figures portent les mêmes signes de référence et présentent des caractéristiques identiques ou analogues, de sorte que ces éléments communs ne sont généralement pas à nouveau décrits, par souci de concision. Les termes « premier(s) » (ou première(s)), « deuxième(s) », etc.) sont utilisés dans ce document par convention arbitraire pour permettre d’identifier et de distinguer différents éléments (tels que des modules, des données, etc.) mis en œuvre dans les modes de réalisation décrits ci-après. Comme illustré en figures 2 et 3, le dispositif de traitement 20 est configuré pour obtenir des premières données notées DT1, ces données étant représentatives d’une réponse ultrasonore 16 d’une région d’intérêt 4 à des impulsions ultrasonores 8 focalisées dans la région d’intérêt 4 en combinaison avec des ondes des sondage 14. Pour ce faire, les impulsions ultrasonores 8 et les ondes de sondage 14 peuvent être émises alternativement au cours du temps dans la région d’intérêt 4. Les impulsions ultrasonores 8 engendrent les ondes de cisaillement 12 dans la région d’intérêt 4 tandis que l’émission d’ondes de sondage 14 permet, à partir de la réponse du milieu, d’analyser la propagation de ces ondes de cisaillement 14 dans ladite région. On suppose par la suite à titre d’exemple que la sonde ultrasonore 2 (figures 1 à 4) est utilisée pour émettre les impulsions ultrasonores 8 en combinaison avec les ondes de sondage 14. A noter toutefois que les ondes 8 et 14 précitées peuvent être émises par tout moyen approprié, commun ou distinct, que ce soit par un seul dispositif ou des dispositifs en charge respectivement de l’émission des impulsions 8 et des ondes de sondage 14. On suppose par la suite que la sonde 2 forme avec le dispositif de traitement 20 un système, dit aussi système de traitement, noté SY2 (figures 3 et 6). Le dispositif de traitement 20 peut obtenir les premières données DT1 de diverses manières. Selon un exemple, le dispositif de traitement 20 peut ainsi extraire les premières données DT1 depuis une mémoire, interne ou externe, ou encore recevoir ces données DT1 depuis l’extérieur, par exemple depuis le système ultrasonore SY1 comprenant la sonde 2 ou depuis tout autre entité appropriée, tel qu’un serveur par exemple. Le dispositif de traitement 20 peut être configuré pour traiter les premières données ultrasonores ST1 en temps réel ou quasi-temps réel, au fur-et-à-mesure que les données sont reçues depuis la sonde 2 (ou plus généralement depuis le système ultrasonore SY1). Selon un autre exemple, les premières données DT1 sont générées lors d’une phase d’acquisition ou de scan ultrasonore, puis ces données DT1 sont reçues et traitées par le dispositif de traitement 20 lors d’une phase de traitement postérieure à la phase d’acquisition (par exemple ultérieurement à une phase de scan échographique). Les ondes de sondage 14 émises par la sonde ultrasonore 2 sont configurées pour détecter et suivre la propagation des ondes de cisaillement 12 causées par les impulsions ultrasonores 8 au sein de la région d’intérêt 4. On suppose par la suite à titre d’exemple que l’on utilise, en tant qu’ondes de sondage 14 ( FIG. 3 ), des ondes planes ultrasonores ce qui permet avantageusement d’obtenir divers angles d’émission et donc de couvrir rapidement une zone appropriée du milieu, bien que d’autres types d’ondes de sondage soient possibles. Les ondes de sondage 14 émises pas la sonde 2 se caractérisent par au moins un paramètre noté PR, dont le nombre et la nature peuvent être adaptés selon le cas. Les ondes de sondage 14 sont émises séquentiellement dans la région d’intérêt 4 selon des valeurs distinctes VL, dites aussi valeurs de paramètre, dudit au moins un paramètre PR. Autrement dit, la valeur du ou des paramètres PR des ondes de sondage 14 sont modifiées (ou adaptées) au cours du temps de sorte à prendre successivement au moins deux valeurs VL distinctes. Par souci de simplification de l’exposé de la présente divulgation, l’on considère par la suite à titre purement illustratif que la valeur VL d’un seul paramètre PR est modifiée pour générer les premières données DT1, bien que d’autres exemples selon lesquels les valeurs VL d’au moins deux paramètres PR des ondes de sondage 14 sont modifiées sont également possibles. Par souci de simplification de l’exposé de la présente divulgation, l’on considère par la suite à titre purement illustratif que la valeur d’un paramètre PR est modifiée au cours d’un scan ultrasonore de sorte que ce paramètre PR prenne séquentiellement trois valeurs distinctes notées VL1, VL2 et VL3. Ces valeurs VL1, VL2 et VL3 peuvent être appliquées selon un schéma répétitif, par exemple en les appliquant successivement de façon cyclique (selon par exemple une séquence « VL1 - VL2 - VL3 - VL1 - VL2 - etc. »). Selon un exemple illustré en FIG. 4 , l’angle d’émission θ selon lequel sont émises les ondes de sondage 14 dans le (ou en direction du) milieu d’intérêt 4 est utilisé en tant que paramètre PR variant au cours du temps. Autrement dit, la sonde 2 peut être configurée pour faire varier l’angle d’émission θ, en tant que paramètre PR, de sorte à ce que ce dernier prenne séquentiellement les valeurs VL1, VL2 et VL3, conduisant ainsi à la génération d’ondes de cisaillement 12 pour respectivement ces valeurs distinctes VL1, VL2 et VL3. Dans l’exemple représenté en FIG. 4 , la valeur VL2 est fixé à zéro ( FIG. 4 (B)) tandis que les valeurs VL1 et VL3 sont fixées à des valeurs non nulles (figures 4(A) et 4(C)). Les valeurs VL1 et VL3 sont par exemple distinctes l’une de l’autre. L’angle d’émission θ peut par exemple varier de sorte à prendre soit une valeur négative (VL1) soit une valeur positive (VL3), ce qui revient à orienter alternativement les ondes de sondage 14 de part et d’autre d’un axe normal à une surface de la région d’intérêt 4. D’autres configurations d’angles d’émission sont toutefois possibles. A noter que d’autres types de paramètre PR des ondes de sondage 14 sont possibles. Ainsi, selon un exemple illustré en FIG. 5 , la fréquence d’émission f selon laquelle sont émises les ondes de sondage 14 dans le milieu d’intérêt 4 est utilisée en tant que paramètre PR variant au cours du temps. Autrement dit, la sonde 2 peut être configurée pour faire varier la fréquence f des ondes de sondage 14 de sorte à ce que cette fréquence prenne séquentiellement des valeurs VL distinctes, notées respectivement VF1, VF2 et VF3 dans cet exemple, ce qui permet ainsi d’observer les ondes de cisaillement 12 selon différentes fréquences des ondes de sondage 14. Dans l’exemple représenté en FIG. 5 , parties (A), (B) et (C), la valeur VL diminue progressivement depuis VL1 jusqu’à VL3, bien que d’autres exemples soient possibles, tels que notamment la configuration inverse où les valeurs VL croissent progressivement de VL1 à VL3. Selon un exemple, ledit au moins un paramètre PR qui varie séquentiellement au cours de la phase de scan ultrasonore, comprend au moins l’un quelconque parmi, ou les deux parmi : - l’angle d’émission θ selon lequel les ondes de sondage 14 sont émises en direction du milieu d’intérêt 4 ; et - la fréquence f selon laquelle sont émises les ondes de sondage 14 en direction du milieu d’intérêt 4. La FIG. 6 représente le dispositif de traitement 20 coopérant avec la sonde 2 selon un exemple de réalisation. Comme illustré, le dispositif 20 est configuré pour recevoir les premières données DT1 produites par la sonde 2 (ou plus généralement par le système ultrasonore SY1 de la FIG. 1 ), par exemple via une liaison de communication 9 qui peut être une liaison filaire ou sans fil selon le cas. Dans l’exemple considéré, le dispositif de traitement 20 et la sonde ultrasonore 2 forment ensemble un (ou font partie d’un) système de traitement noté SY2 (figures 3 et 6). A titre d’exemple, on considère par la suite que le système ultrasonore SY1 ( FIG. 1 ) et le système de traitement SY2 (figures 3 et 6) forment un seul et même système. Le dispositif de traitement 20 et le dispositif de contrôle 5 ( FIG. 1 ) peuvent par exemple former un seul et même dispositif. Dans ce cas, le dispositif de traitement 20 est configuré d’une part pour piloter la sonde ultrasonore 2 par l’envoi de signaux électrique SG1 de sorte à émettre les ondes 8 et 14, et d’autre part pour traiter les premières données DT1, représentatives de la réponse ultrasonore 16 du milieu 4, captées en retour par la sonde ultrasonore 2. D’autres exemples dans lesquels le dispositif de traitement 20 et le dispositif de contrôle 5 sont distincts sont également possibles. Le dispositif de traitement 20 peut par exemple être un serveur, ou tout autre moyen informatique, apte à obtenir et traiter les premières données DT1. Comme illustré en FIG. 6 selon un exemple, le dispositif de traitement 20 comprend un processeur 22 et une mémoire 24. La mémoire 24 est apte à stoker sous forme d’un programme d’ordinateur PG1 des instructions définissant les étapes des procédés décrits dans la présente divulgation. A ce titre, la mémoire 24 constitue un support d’enregistrement (ou support d’informations) conforme à des modes particuliers de réalisation, lisible par le dispositif de traitement 20, et sur lequel est enregistré un programme d’ordinateur PG1 conforme à des modes particuliers de réalisation. Ce programme d’ordinateur PG1 comporte des instructions pour l’exécution des étapes d’un procédé de traitement (dit aussi procédé de traitement d’élastographie ultrasonore) dont des modes particuliers de réalisation sont décrits dans la présente divulgation. Le processeur 22 est ainsi configuré pour exécuter les instructions du programme d’ordinateur PG1 afin de réaliser des étapes du procédé de traitement de la divulgation. Selon les configurations et types de dispositif informatique considérés, la mémoire 24 peut comprendre de la mémoire volatile (telle que la RAM) et/ou non volatile (telle que de la mémoire ROM, flash, EEPROM, etc., ou tout autre dispositif de stockage et/ou support lisible par ordinateur comme décrit ci-après, ou une combinaison parmi ces exemples de catégories de mémoires. La mémoire 24 peut, par exemple, être gérée en mode DMA (pour « Direct Memory Access » en anglais). Cette mémoire 24 utilisée par le dispositif de traitement 20 peut par exemple comprendre tout ou partie d’une mémoire de carte graphique (ou carte vidéo), ce type de mémoire étant notamment apte à traiter et/ou envoyer des données d’image pouvant être utilisées pour afficher une ou des images sur un écran (ou une unité) d’affichage. La mémoire 24 est apte à stocker diverses données pouvant être utilisées du dispositif de traitement 20 au cours du procédé de traitement de la divulgation, telles que notamment les premières données DT1 susmentionnées, des deuxièmes données ultrasonores DT2 générées à partir des premières données DT1, et des estimations de vitesse, à savoir au moins une célérité individuelle VI (par exemple VI1, VI2 et VI3) représentative de la célérité de déplacement des ondes de cisaillement 12 dans la région d’intérêt 4 et au moins une célérité composée VC des ondes de cisaillement 12. La nature et la fonction de ces données apparaîtront plus précisément ci-après, dans des exemples de réalisation du dispositif de traitement 20 et du procédé de la divulgation. Le dispositif de traitement 20 peut prendre la forme de divers moyens informatiques appropriés (tels que par exemple station de travail, ordinateur, serveur, etc.) comprenant tout ou partie des éléments décrits ci-dessus ainsi qu’éventuellement d'autres éléments non mentionnés. De façon générale, le dispositif de traitement 20 comprend des moyens appropriés pour mettre en œuvre les étapes du procédé de traitement tel que décrit ci-après dans des exemples de réalisation de la présente divulgation. Le système ultrasonore SY2 (figures 3 et 6) peut être un système d’imagerie ultrasonore, par exemple utilisé dans le domaine médical ou vétérinaire. Les images ultrasonores générées par le système SY2 peuvent être soit analysées en temps réel, par exemple par un utilisateur ou un algorithme, et/ou un module d’intelligence artificielle, soit analysées ultérieurement et/ou dans un autre emplacement que celui dans lequel se trouve le système SY2. Le système ultrasonore SY1 peut être un système médical à ultrasons. De même, la sonde 20 peut être une sonde médicale à ultrasons. Aussi, le système de traitement SY2 peut être un système appliqué dans le domaine médical, bien que d’autres applications soient possibles. A titre d’exemple, le système SY2 peut être associé à la sonde à ultrasons 2, afin d’étudier un milieu 4 ( FIG. 2 ), notamment pour collecter et traiter des premières données DT1 caractérisant un tel milieu 4 selon un mode d’élastographie par ondes de cisaillement. Le milieu 4 ainsi observé peut être de diverses natures selon le cas. Il peut s’agir par exemple de métaux, de tissus d’êtres vivants, en particulier des tissus humains ou tissus d’animaux. L’observation d’un milieu 4 comprenant une ou des structures minérales par exemple est également possible (graviers, volcan, cartographie d’un sol, par exemple de fonds marins, etc.). Le système de traitement SY2 peut donc être configuré pour diverses applications, notamment dans les domaines de l'acoustique, de l’étude des matériaux, de l'imagerie médicale, vétérinaire, et/ou de la biomédecine. Le système SY2 peut être mis en œuvre dans des applications autres que l’imagerie ultrasonore, par exemple pour caractériser une région d’intérêt à partir des premières données DT1 sans pour autant générer d’image ultrasonore. La FIG. 7 représente schématiquement un exemple de réalisation non limitatif du dispositif de traitement 20, et plus généralement du système de traitement SY2, tel que précédemment décrit conjointement aux figures 1-6, à savoir en l’espèce un dispositif d'imagerie ultrasonore. Dans cet exemple, le dispositif de traitement 20 est compris dans le corps d’une station de commande, cette dernière étant en outre munie d’une interface de commande et d’un ou plusieurs dispositifs d’affichage. La station de commande peut également comprendre le dispositif de contrôle 5 ( FIG. 1 ) précédemment décrit. La sonde 2 est par exemple reliée de façon communicante avec le dispositif de traitement 20 (et le dispositif de contrôle 5) via une liaison 9 qui prend ici la forme d’un câble de connexion, bien que des variantes soient possibles dans lesquelles par exemple la connexion est assurée via une liaison sans fil. Le système de traitement SY2 peut être configuré pour réaliser une image ultrasonore de diverses natures, par exemple une image dite de type B-mode du milieu M (image en mode B, présentée en niveaux de gris), une image dite Doppler illustrant les mouvements de fluides dans le milieu observé, et/ou une image montrant une caractéristique mécanique du milieu (par exemple des données d'image d'élastographie obtenues à l’aide d’ondes de cisaillement (« ShearWave ® Elastography »). Selon un exemple, le signal électrique SG1 transmis par le dispositif de contrôle 5 cause ainsi l’émission par la sonde ultrasonore 2 d’ondes W, à savoir des impulsions ultrasonores 8 ( FIG. 3 ) en combinaison avec des ondes de sondage 14 ( FIG. 4 ). Comme représenté en FIG. 8 selon un mode de réalisation, le processeur 22 du dispositif de traitement 20, piloté par le programme d’ordinateur PG1 ( FIG. 6 ), peut mettre en œuvre un certain nombre de modules, à savoir : un module d’obtention MD2, un module de génération MD4, un premier module d’estimation MD6 et un deuxième module d’estimation MD8, voire éventuellement aussi un troisième module MD10 d’estimation. Plus spécifiquement, le module d’obtention MD2 est configuré pour obtenir des premières données DT1 représentatives d’une réponse ultrasonore 16, d’un milieu d’intérêt 4, à des impulsions ultrasonores 8 focalisées dans ledit milieu en combinaison avec des ondes de sondage 14 émises dans ledit milieu séquentiellement selon des valeurs distinctes VL (par exemple VL1 à VL3) d’au moins un paramètre PR (par exemple l’angle d’émission θ) desdites ondes de sondage (figures 1-6). Le module de génération MD4 est configuré pour générer, à partir des premières données DT1 respectivement pour chaque valeur VL dudit au moins un paramètre PR, des deuxièmes données DT2 représentatives d’une propagation d’ondes de cisaillement 12 dans le milieu d’intérêt 4 sous l’effet des impulsions ultrasonores 8. Le premier module d’estimation MD6 est configuré pour estimer, à partir des deuxièmes données DT2 associées à chaque valeur VL dudit au moins un paramètre PR, une célérité individuelle VI (par exemple VI1 à VI3, respectivement) représentative de la célérité de déplacement des ondes de cisaillement 12 pour respectivement ladite valeur (par exemple VL1 à VL3, respectivement). Le deuxième module d’estimation MD8 est configuré pour estimer au moins une célérité VC (par exemple VC1 à VC3, respectivement), dite célérité composée, des ondes de cisaillement 12 par une combinaison des célérités individuelles VI. Le troisième module d’estimation MD10 est, le cas échéant, configuré pour estimer, à partir de ladite au moins une célérité composée VC, une dureté (ou élasticité) du milieu d’intérêt 4. La configuration et le fonctionnement des modules MD2, MD3 et MD4 du dispositif de contrôle 10 apparaîtront plus précisément dans les exemples de réalisation décrits ci-après en référence aux figures. Les modules MD2, MD3 et MD4 tels que représentés en FIG. 4 ne représentent qu’un exemple de mise en œuvre non limitatif de la divulgation. De manière générale, pour chaque étape du procédé de traitement de la présente divulgation, le dispositif de traitement 20 de la divulgation peut comprendre un module correspondant configuré pour réaliser ladite étape (et réciproquement). Des modes de réalisation du procédé de traitement de la divulgation sont à présent décrits en référence aux figures 9 à 16. Dans ces exemples, le procédé de traitement est mis en œuvre par le dispositif de traitement 20, et plus généralement par le système de traitement SY2, tels que précédemment décrits conjointement aux figures 1 à 8. Pour ce faire, le dispositif de traitement 20, coopérant avec la sonde ultrasonore 2, peut exécuter le programme d’ordinateur PG1. Au cours d’une étape S2 d’obtention ( FIG. 9 ), le dispositif de traitement 20 obtient des premières données DT1 représentatives d’une réponse ultrasonore 16, d’un milieu (ou région) d’intérêt 4, suite à émission d’impulsions ultrasonores 8 focalisées dans ledit milieu 4 en combinaison avec des ondes de sondage 14 émises dans ledit milieu 4 séquentiellement selon des valeurs distinctes VL d’au moins un paramètre PR des ondes de sondage 14. En d’autres termes, les premières données DT1 représentent des échos ultrasonores 16 émis par le milieu 4 en réponse aux impulsions ultrasonores 8 combinées aux ondes de sondage 14 (figures 2 et 3). Ces ondes de sondage 14 se caractérisent par au moins un paramètre PR qui varie au cours du temps. Ainsi, la valeur VL dudit au moins un paramètre PR varie séquentiellement au cours de la phase de scan ultrasonore de sorte à prendre successivement différentes valeurs. La valeur VL du ou des paramètres PR peut être modifiée sous le contrôle par exemple du dispositif de contrôle 5 ( FIG. 1 ). Plus particulièrement, la réponse ultrasonore 16 définie par les première données DT1 caractérise des ondes de cisaillement 12 provoquées par les impulsions ultrasonores 8 focalisées dans la région d’intérêt 4. Ces impulsions ultrasonores 8, émises par exemple sous forme de tirs d’impulsions, forment des ondes de compression causant des contraintes mécaniques dans la région d’intérêt 4. Dans la région du point focal 10 ( FIG. 3 ), se produit un phénomène de pression de radiation acoustique selon lequel un transfert d’énergie engendre des déplacements de matière dans le milieu examiné, ce qui génère des ondes de cisaillement 12 se propageant en particulier latéralement autour du point focal 10. Les premières données DT1 sont acquises par la sonde ultrasonore 20 au cours d’une phase d’acquisition, correspondant à un ou plusieurs scans ultrasonores en mode d’élastographie par ondes de cisaillement, selon un principe décrit ci-avant en référence aux figures 2 et 3. Ces premières données DT1 sont ensuite transmises, depuis cette sonde 20, au dispositif de traitement 20, par exemple via la liaison de communication 9. Ces premières données DT1 peuvent toutefois être obtenues par le dispositif de traitement 20 par d’autres manières selon le cas. En particulier, le scan ultrasonore peut par exemple être réalisé dans un premier temps, puis le dispositif de traitement 20 peut obtenir ces premières données DT1 ultérieurement, selon une quelconque façon appropriée. Par exemple, le dispositif de traitement 20 peut extraire les premières données DT1 depuis une mémoire, interne ou externe, ou encore recevoir ces données DT1 depuis l’extérieur, par exemple depuis le système ultrasonore SY1 comprenant la sonde 2 ou depuis tout autre entité appropriée, tel qu’un serveur par exemple. On suppose par la suite à titre d’exemple que les ondes de sondage 14 sont des ondes planes ultrasonores, bien que d’autres exemples soient possibles. Selon un autre exemple, les ondes de sondage 14 sont des ondes ultrasonores de type divergentes ou focalisées. Selon un exemple, les ondes de sondage 14 comprennent au moins l’un quelconque parmi les types d’ondes précités. On suppose par la suite à titre purement illustratif que l’angle d’émission θ selon lequel sont émises les ondes de sondage 14 est utilisé en tant que paramètre PR des ondes de sondage 14. A noter toutefois que le procédé de traitement de la divulgation peut être réalisée de façon analogue, selon un exemple de réalisation, en faisant varier la valeur de la fréquence f des ondes de sondage 14 en tant que paramètre PR. Toutefois, dans la suite de la description, à des fins de simplifications, seul l’exemple basé sur l’angle d’émission θ en tant que paramètre PR sera illustré. Ainsi, au cours de la phase d’acquisition, la sonde ultrasonore 2 émet les ondes de sondage 14 en faisant varier au cours du temps la valeur VL de l’angle d’émission θ (figures 4 et 6). Selon l’exemple considéré, la sonde 20 émet séquentiellement les ondes de sondage 14 selon N valeurs distinctes VL de l’angle d’émission θ (où N est un entier au moins égal à 2), ce qui se traduit par une variation de l’orientation des ondes de sondage 14 au cours du temps. A titre d’exemple, on suppose par la suite que N=3 et que l’angle d’émission θ prend successivement les valeurs distinctes VL1, VL2 et VL3. Ces valeurs VL1, VL2 et VL3 peuvent par exemple être appliquées de manière répétée (selon un schéma répétitif), par exemple en les affectant successivement de façon cyclique à l’angle d’émission θ (selon par exemple une séquence « VL1 - VL2 - VL3 - VL1 - VL2 - etc. »). Selon un exemple, préalablement à l’étape S2 d’obtention des premières données DT1 par le dispositif de traitement 20, la sonde 20 (ou plus généralement le système de traitement SY2) réalise les étapes suivantes : - formation des ondes de cisaillement 12 dans le milieu d’intérêt 4 par focalisation des impulsions ultrasonores 16 ; et - émission, dans le milieu d’intérêt 4, des ondes de sondage 14 configurées séquentiellement selon les valeurs distinctes VL1-VL3 du paramètre PR (à savoir l’angle d’émission θ par exemple) pour détecter la réponse ultrasonore 16 du milieu d’intérêt 4 pour chaque valeur VL du paramètre PR. Au cours d’une étape S4 de génération ( FIG. 9 ), le dispositif 20 génère, à partir des premières données DT1, respectivement pour chaque valeur VL du paramètre PR (c’est-à-dire l’angle d’émission θ dans cet exemple), des deuxièmes données DT2 représentatives d’une propagation d’ondes de cisaillement 12 dans le milieu d’intérêt 4 sous l’effet des impulsions ultrasonores 8 (figures 2-3). Autrement dit, le dispositif 20 traite les première données DT1 et, à partir de ce traitement, génère des deuxièmes données DT2 pour respectivement les valeurs VL1, VL2 et VL3 de l’angle d’émission θ des ondes de sondage 14. Ces deuxièmes données DT2 représentent le déplacement des ondes de cisaillement 12 dans le milieu 4. La configuration de ces deuxièmes données peut varier selon le cas. Selon un exemple illustré en FIG. 10 , ces deuxièmes données DT2 définissent, pour chaque valeur VL1, VL2 et VL3 de l’angle d’émission θ, une séquence d’images IM1 représentative de la propagation des ondes de cisaillement 12 causées par les impulsions ultrasonores 8 dans le milieu d’intérêt 4. Le dispositif 20 génère ainsi une pluralité de séquences d’images IM1 en association avec respectivement les différentes valeurs VL de l’angle d’émission θ. Chaque image (ou trame) IM1 peut représenter une position des ondes de cisaillement à un instant donné. Deux images IM1 consécutives représentent alors une estimation de déplacement entre deux acquisitions correspondantes pour une valeur VL donnée du paramètre PR considéré. Ainsi, une série d’image IM1 pour une valeur VL donnée peut représenter un déplacement du milieu (et des ondes de cisaillement 12) au cours du temps. A partir d’un déphasage (ou de changements) entre deux image IM1 successives pour une valeur VL donnée, le dispositif peut déduire un déplacement des ondes de cisaillement 12. A noter toutefois que les images IM1 ainsi produites ne sont pas nécessairement générées, stockées et/ou affichées ou restituées à un utilisateur, bien que cela soit possible. Au cours d’une étape S6 d’estimation ( FIG. 9 ), le dispositif 20 estime, à partir des deuxièmes données DT2 associées à chaque valeur VL du paramètre PR (c’est-à-dire l’angle d’émission θ dans cet exemple), une célérité individuelle VI représentative de la célérité de déplacement des ondes de cisaillement 12 pour respectivement ladite valeur VL. Ainsi, le dispositif 20 détermine, sur la base des deuxièmes données DT2, les célérités individuelles VI1, VI2 et VI3 représentatives des célérités des ondes de cisaillement 12 pour respectivement les valeurs VL1, VL2 et VL3 de l’angle d’émission θ ( FIG. 10 ). Les célérités individuelles VI sont par exemples déterminées (S6, FIG. 9 ) par analyse des images IM1 définies par les deuxièmes données DT2. Par comparaison d’une série d’images IM1 obtenues pour une valeur VL donnée, par exemple en estimant un déphasage (ou des changements) entre deux images IM1 successives pour la valeur VL donnée, le dispositif 20 peut calculer une vitesse VI correspondante en chaque point de la région d’intérêt 4. Au cours d’une étape S8 d’estimation ( FIG. 9 ), le dispositif 20 estime au moins une célérité VC, dite célérité composée, des ondes de cisaillement 12 par une combinaison des célérités individuelles VI précédemment obtenues (en S6). Au moins deux célérités individuelles VI estimées en S6 sont combinées pour déterminer la célérité composée VC. Ainsi, le dispositif 20 combine dans cet exemple les vitesses individuelles VI1, VI2 et VI3, associées respectivement aux valeurs VL1, VL2 et VL3, pour en déduire la célérité composée VC. La FIG. 11 représente un exemple de réalisation des étapes S6 et S8 du procédé de traitement. Le dispositif 20 détermine en S6, par traitement des premières données DT1 reçues, des deuxièmes données DT2 définissant une séquence d’images IM1 correspondant successivement aux valeurs VL1, VL2 et VL3, de façon cyclique, au cours du temps. Le dispositif 20 a besoin d’au moins deux images (ou trames) IM1 pour une valeur VL donnée pour pouvoir les comparer et en déduire une vitesse individuelle VI correspondante pour ladite valeur V. Pour ce faire, le dispositif 20 peut donc stocker en mémoire les images IM1 qu’il détermine au cours du temps et comparer une pluralité d’images IM1 obtenues pour une même valeur VL donnée pour en déduire une vitesse individuelle VL correspondante. Ce processus peut être réalisé de façon analogue pour chaque valeur VL1, VL2 et VL3 afin d’en déduire respectivement des vitesse individuelles VI1, Vl2 et -VI3. Ces vitesses individuelles VI1, Vl2 et -VI3 sont ensuite combinées pour déduire une célérité composée VC au cours de l’étape S8 d’estimation. La manière dont les célérités individuelles VI obtenues en S6 sont combinées (S8) pour déterminer la célérité composée VC peut varier selon le cas. Selon un exemple, la combinaison des célérités individuelles VI1, Vl2 et -VI3 est déterminée par un calcul d’une moyenne des célérités individuelles VI pour estimer la célérité composée VC. On suppose par la suite que le dispositif 20 calcule une célérité composée VC à partir d’une moyenne des célérités individuelles VI1, Vl2 et VI3. Un traitement analogue peut être réalisé par le dispositif 20 pour estimer une pluralité de célérités composées VC. En particulier, des célérités composées VC peuvent être déterminées (S8) en association avec des positions respectives du milieu d’intérêt 4. A noter que des calculs autres que des moyennes, combinant également les célérités individuelles VI pour déterminer la célérité composée VC, sont possibles. La moyenne calculée en S8 peut être une moyenne pondérée selon des poids K attribués respectivement à chaque célérité individuelle VI, ce qui permet avantageusement de réaliser une combinaison intelligente des estimations de vitesse selon différentes valeurs VL du paramètre PR considéré. Selon un exemple, le dispositif 20 attribue (S8a, FIG. 9 ) ainsi des poids (ou coefficients) K1, K2 et K3 en association respectivement avec les célérités individuelles VI1, VI2 et VI3. Ces poids K1-K3 sont déterminés en fonction d’au moins une caractéristique des deuxièmes données DT2 générées respectivement pour chaque valeur VL1, VL2 et VL3 du paramètre PR considéré. Le dispositif 20 estime la célérité composée VC par calcul de la moyenne en pondérant les célérités individuelles VI1, Vl2 et -VI3selon les poids associés K1, K2 et K3. En d’autres termes, la célérité composée VC peut alors être définie selon l’expression suivante : [Math 3] (equation 3) La ou les caractéristiques des deuxièmes données DT2, qui sont utilisées pour déterminer les poids de pondération K1-K3 peuvent être adaptées selon le cas. Ces caractéristiques sont par exemple représentatives d’un niveau de qualité des deuxièmes données DT2, ou encore des images IM1 définies sur la base de ces deuxièmes données DT2 (figures 9-10). Il est ainsi avantageusement possible d’adapter la prépondérance ou poids, dans la moyenne composée VC, des célérités individuelles VI en fonction de la pertinence ou qualité des réponses ultrasonores 16 correspondantes obtenues selon les diverses valeurs VL du paramètre PR considéré (à savoir de l’angle d’émission θ dans cet exemple). On peut par exemple ne prendre en compte que les célérités individuelles VI pertinentes pour chaque pixel des images IM1. Cette pondération est rendue possible par le fait que l’on obtient une pluralité de célérités individuelles VI1, Vl2 et VI3 associées à des valeurs respectives VL1, VL2 et VL3 du paramètre PR, ce qui offre un volume conséquent de données exploitables parmi lesquelles on peut avantageusement sélectionner les célérités individuelles VI les plus pertinentes et adapter le cas échéant l’importance de chaque célérité individuelle VI dans le calcul de la célérité composée VC. Selon un exemple, au moins deux célérités individuelles parmi VI1, Vl2 et VI3 sont sélectionnées en S8 ( FIG. 9 ) pour être prises en compte dans le calcul de la valeur composée VC. Au moins une célérité individuelle VL, par exemple VL1, peut être exclue de la combinaison des célérités individuelles VI réalisée en S8, par exemple s’il est détecté que cette célérité individuelle VL est associée à des deuxièmes données DT2 qui ne respectent pas au moins un critère de qualité prédéterminé. Au moins un poids K peut ainsi être fixé à zéro (K = 0) pour exclure une célérité individuelle VI correspondante du calcul de la célérité composée VC, ce qui permet avantageusement de ne retenir que les célérités individuelles VI les plus pertinentes dans le calcul de la moyenne en S8b ( FIG. 9 ). Selon un exemple, au cours de l’étape S8 d’estimation, le dispositif 20 détermine des variances spatiales caractérisant respectivement les deuxièmes données DT2 associées à chaque valeur VL du paramètre PR considéré. Le dispositif 20 détermine ensuite, à partir des variances spatiales ainsi obtenues, des scores de confiance associés respectivement aux célérités individuelles VI. Un tel score de confiance, associé à une célérité individuelle VI, représente une probabilité que la célérité individuelle VI est fiable (ou acceptable). Le poids K respectif attribué à chaque célérité individuelle VI peut alors être fonction du score de confiance associé. Il est ainsi avantageusement possible de pondérer l’importance donnée à chaque célérité individuelle VI dans la détermination de la célérité composée VC, de façon à obtenir une estimation précise et fiable de la célérité composée des ondes de cisaillement 12. Selon un exemple, au cours de l’étape S8 d’estimation, le dispositif 20 compare chaque score de confiance ainsi obtenu avec une valeur de référence (par exemple une valeur seuil). Si un score de confiance associée à des deuxièmes données DT2 correspondant à une valeur VL ne satisfait pas à un au moins un critère de qualité (fonction du résultat de cette comparaison), le poids K attribué à la célérité individuelle VI associée à cette valeur VL est nul (K est fixé à 0) de sorte à exclure ladite célérité individuelle VI de l’estimation (S8) de la célérité composée VC. Ainsi, si un score de confiance ( FIG. 10 ) correspondant à une valeur VL donnée de l’angle d’émission θ est par exemple inférieure à une valeur seuil, la célérité individuelle VI associée n’est pas prise en compte dans le calcul (S8) de la célérité composée VC, ce qui permet avantageusement de maximiser la fiabilité et la précision du résultat ainsi obtenu. En prenant en compte les scores de confiance dans l’estimation de la célérité composée VI, il est ainsi possible de sélectionner les célérités individuelles VI les plus pertinentes et ainsi avantageusement minimiser la probabilité d’obtenir des « pixels morts » dans les images IM2, c’est-à-dire des pixels pour lesquels les données d’image sont considérés comme non exploitables pour estimer la célérité composée VC des ondes de cisaillement 12. La ou les célérités composées VC déterminées en S8 ( FIG. 9 ) peuvent ensuite être exploitées de diverses manières selon le cas. Selon un exemple, au cours d’une étape S10 d’estimation ( FIG. 9 ), le dispositif 20 réalise une estimation, à partir de la ou les célérités composées VC, d’une dureté (ou rigidité, ou élasticité) du milieu d’intérêt 4. Selon un exemple, le dispositif 20 déclenche une restitution (visuelle, vocale, etc.), par exemple par affichage, de la ou les estimations de célérités composées VC. Les célérités composées VC sont par exemple affichées sur un écran sous forme d’au moins une image IM2. Cet écran peut être compris dans le système de traitement SY2, comme illustré à titre d’exemple en FIG. 7 ; il peut être déporté dans le même environnement ou dissocié et/ou éloigné du système de traitement. Selon un exemple, le dispositif 20 cause l’affichage de la ou les célérités composées, et/ou de la ou les duretés estimées à partir des célérités composées, sous forme d’au moins une image IM2 (par exemple sous la forme d’un élastogramme), par exemple sous la forme d’une séquence d’images IM2 pouvant éventuellement former une vidéo qu’un utilisateur peut visionner, soit immédiatement après obtention de la séquence d’images, sous en quasi simultané soit de manière dissociée. On suppose par la suite que le dispositif 20 génère en sus en S10 des troisièmes données DT3 représentatives de célérités composées VC déterminées en S8 ( FIG. 9 ). La nature de ces troisièmes données DT3 peut varier selon le cas. Ces données DT3 peuvent notamment prendre la forme de données d’image définissant des images IM2 représentatives des célérités composées VC déterminées en S8 à différentes positions de la région d’intérêt 4. Comme illustré en FIG. 12 selon un exemple, les troisièmes données DT3 générées en S10 ( FIG. 9 ) peuvent définir une carte paramétrique (ou image cartographique) IM2 représentant, en différents pixels de ladite carte, la célérité composée VC des ondes de cisaillement 12 et/ou la dureté correspondante du milieu 4. Chaque pixel de cette image cartographique IM2 correspond à une position respective de la région d’intérêt 4. Selon un exemple, au cours de l’estimation S8 ( FIG. 9 ), une pluralité de célérités composées VC est estimée en différentes positions du milieu d’intérêt 4. Le dispositif 20 génère (S10) alors, à partir des célérités composées VC, des données cartographiques DT3 représentatives des célérités composées VC, et/ou de duretés estimées à partir des célérités composées VC, dans respectivement les différentes positions du milieu d’intérêt 4. Ces images IM2 peuvent alors être affichées ou restituées selon l’un des modes par exemple indiqué ci-dessus. Comme déjà décrit en référence à la FIG. 11 , les deuxièmes données DT2 déterminées en S6 ( FIG. 9 ) peuvent définir une séquence d’images IM1 correspondant successivement aux valeurs VL1, VL2 et VL3, ces valeurs VL se répétant de façon cyclique au cours du temps. On peut ainsi répéter les étapes S2 à S10 ( FIG. 9 ) au cours du temps pour déterminer, à partir de chaque nouvelle acquisition ultrasonore, au moins une vitesse composée VC. Les figures 13, 14 et 15 représentent schématiquement le déroulement du procédé de traitement de la FIG. 9 selon un exemple de réalisation. Au cours d’une phase préliminaire PH1 ( FIG. 13 ), les étapes S2 à S8 sont réalisées pour trois valeurs distinctes VL1, VL2 et VL3 de l’angle d’émission θ utilisé comme paramètre variable PR des ondes de sondage 14. Au cours de cette phase préliminaire PH1a, le dispositif 20 détermine (S4) des deuxièmes données DT2 définissant au moins deux images IM1 pour chaque valeur VL1-VL3 du paramètre PR, afin d’en déduire (S6) par comparaison des vitesse individuelles VI1-VI3, respectivement. A des fins de simplification de l’exposé de la divulgation, on suppose par la suite que le dispositif 20 estime chaque vitesse individuelle VI par comparaison de deux images IM1 seulement pour une valeur VL donnée, bien qu’un plus grand nombre d’images IM1 puissent être traitées pour une même valeur VL pour en déduire une vitesse individuelle VI correspondante. Le dispositif 20 peut ainsi estimer une vitesse composée VC, nommée VC0, à partir des vitesse individuelle VI1-VI3 selon l’étape S8 ( FIG. 9 ) précédemment décrite. Après cette phase préliminaire PH1, au moins une itération périodique d’un cycle S15 ( FIG. 14 ) est réalisée par le dispositif 20. Au cours de ce cycle S15, le dispositif 20 détermine de nouvelles images IM1 pour respectivement les valeurs VL1-VL3 et estime de nouvelles célérités individuelles VI1-V3 par réitération des étapes S2 à S6 ( FIG. 9 ) en utilisant les nouvelles image IM1 ainsi que des images IM1 obtenues antérieurement lors de la phase préliminaire PH1. A partir des nouvelles célérités individuelles VI1-VI3 ainsi obtenues, le dispositif 20 peut alors en déduire une nouvelle vitesse composée VC, nommée VC1, en réitérant l’étape S8 d’estimation ( FIG. 9 ) comme déjà décrit. Comme illustré en FIG. 15 , le cycle S15 peut ainsi être répété autant que nécessaire pour mettre à jour la vitesse composée VC au cours du temps, et ce en chaque position cible de la région d’intérêt 4. Par exemple, au cours d’une deuxième itération du cycle S15, le dispositif 20 détermine de nouvelles images IM1 pour respectivement les valeurs VL1, VL2 et VL3 et estime de nouvelles célérités individuelles VI1, Vl2 et V3 par réitération des étapes S2 à S6 ( FIG. 9 ) en utilisant les nouvelles images IM1 ainsi que des anciennes images IM1 obtenues lors de la précédente itération du cycle S15. A partir des nouvelles célérités individuelles VI1, Vl2 etVI3 ainsi obtenues, le dispositif 20 peut alors en déduire une nouvelle vitesse composée VC, nommée VC2, en réitérant l’étape S8 d’estimation ( FIG. 9 ) comme déjà décrit. Ce processus itératif peut être réalisé pour estimer des célérités composées VC en de multiples points (ou positions) de la région d’intérêt 4. De cette manière, on peut avantageusement générer rapidement de nouvelles estimations de célérités composées VC avec un minimum de latence tout en conservant une bonne qualité d’estimation. Cette technique peut être utilisée pour générer en temps réel des estimations de célérité. A chaque nouvelle réponse ultrasonore 16 ( FIG. 3 ) obtenue pour différentes valeurs VL1, VL2 et VL3, on peut déduire des deuxièmes données DT2 et utiliser celles-ci en combinaison avec des deuxièmes données DT2 d’au moins une précédente itération du cycle, de sorte à garantir la robustesse du résultat tout en minimisant les temps de latence. Selon un exemple illustré en figures 13 et 16, on suppose que le dispositif 20 réalise la phase préliminaire PH1 comme précédemment décrit en référence à la FIG. 13 . Une fois cette phase préliminaire PH1 effectuée, au moins une itération périodique d’un cycle S20 est réalisée par le dispositif 20 ( FIG. 6 ). Au cours de ce cycle S20, le dispositif 20 détermine une nouvelle image IM1 pour une valeurs VL donnée, par exemple VL1 comme illustré en FIG. 16 , par réitération des étapes S2 et S4 ( FIG. 9 ) comme précédemment décrit. Le dispositif 20 estime ensuite une nouvelle célérité individuelle VI1 par réitération de l’étape S6 d’estimation ( FIG. 9 ) en utilisant la nouvelle image IM1 ainsi que les anciennes images IM1 précédemment obtenues pour la valeur VL1 lors de la phase préliminaire PH1. Le dispositif 20 peut alors estimer une nouvelle célérité composée VC, nommée VC1, par combinaison de la nouvelle célérité individuelle VI1 avec les anciennes célérités individuelles VI2 et VI3 précédemment obtenues lors de la phase préliminaire PH1 pour respectivement les valeurs VL2 et VL3. Comme illustré en FIG. 16 , le cycle S20 peut ainsi être répété autant que nécessaire pour mettre à jour successivement au moins une vitesse individuelle VI, et ainsi mettre à jour la vitesse composée VC au cours du temps, et ce en chaque position cible de la région d’intérêt 4. Par exemple, au cours d’une deuxième itération du cycle S20, le dispositif 20 détermine une nouvelle image IM1 pour respectivement la valeur VL2 et estime une nouvelle célérité individuelle VI2 par réitération des étapes S2 à S6 ( FIG. 9 ) en utilisant la nouvelle image IM1 obtenu lors de cette deuxième itération ainsi que les anciennes images IM2 obtenues, pour la valeur VL2, au cours respectivement de la précédente itération du cycle S20 et de la phase préliminaire PH1. Cette nouvelle célérité individuelle VI2 peut alors être à nouveau combinée avec les anciennes célérités individuelles VI1 et VI3 utilisées lors de la précédente itération du cycle S20. Ce processus itératif peut être réalisé pour estimer des célérités composées VC en de multiples points de la région d’intérêt 4. De cette manière, on peut avantageusement générer rapidement de nouvelles estimations de célérités composées VC avec un minimum de latence tout en conservant une bonne qualité d’estimation. Cette technique peut être utilisée pour générer en temps réel des estimations de célérité(s). A chaque nouvelle réponse ultrasonore 16 (cf. FIG. 3 ) obtenue pour une valeur VL donnée, on peut combiner au moins une célérité individuelle VI additionnelle avec au moins une célérité individuelle estimée lors d’une itération précédente du cycle, de sorte à garantir la robustesse du résultat tout en minimisant les temps de latence. Les exemples ci-dessus de mise en œuvre du procédé de la présente divulgation sont décrits dans le cas où l’angle d’émission θ varie en tant que paramètre PR. A noter toutefois que la divulgation peut s’appliquer de façon analogue en faisant varier, en tant que paramètre PR, un autre paramètre PR des ondes de sondage 14, par exemple en faisant varier la fréquence f ( FIG. 3 ) des ondes de sondage 14 émises dans le milieu 4 pour suivre la propagation des ondes de cisaillement 12. Selon un exemple, des ondes planes ultrasonores sont utilisées en tant qu’ondes de sondage 14 ( FIG. 3 ) et la fréquence f varie séquentiellement en tant que paramètre PR des ondes de sondage 14. Selon un exemple, d’une part des ondes planes ultrasonores sont utilisées en tant qu’ondes de sondage 14 ( FIG. 3 ), et d’autre part la fréquence f ainsi que l’angle d’émission θ varient séquentiellement en tant que paramètres PR des ondes de sondage 14. Selon un exemple, des ondes divergentes ultrasonores sont utilisées en tant qu’ondes de sondage 14. Cette implémentation est particulièrement adaptée dans le cas où l’on utilise une sonde ultrasonore 2 de profil courbe. Selon un exemple, des ondes focalisées ultrasonores sont utilisées en tant qu’ondes de sondage 14, ce qui permet notamment avantageusement d’observer avec précision une région d’intérêt 4 de relativement petite taille. La présente divulgation permet avantageusement de traiter des données d’élastographie ultrasonore DT1 de façon efficace pour permettre une estimation précise et fiable de la célérité d’ondes de cisaillement 12 se propageant dans un milieu d’intérêt 4. En particulier, la divulgation prévoit le traitement de premières données DT1 représentatives d’une réponse ultrasonore 16, d’un milieu d’intérêt 4, à des impulsions ultrasonores 8 focalisées dans le milieu en combinaison avec des ondes de sondage 16 dont on fait varier au moins un paramètre PR au cours du temps. Comme précédemment décrit, la technique dite de moyennage des données ultrasonores implique de déterminer des données ultrasonores en moyennant des réponses ultrasonores d’un milieu d’intérêt en association avec plusieurs valeurs de l’angle d’émission θ des ondes de sondage 14. Cette technique permet de déterminer une unique estimation de célérité des ondes de cisaillement 12 pour une position donnée du milieu 4, ce qui ne permet pas toujours de caractériser de façon satisfaisante la vitesse des ondes de cisaillement, et donc la rigidité du milieu. Si sa qualité est insuffisante, la célérité individuelle ainsi obtenue ne peut être exploitée de façon satisfaisante, ce qui dégrade le résultat de l’estimation de célérité et entrave l’observation du milieu 4. En cas d’application d’imagerie par exemple, l’image ainsi obtenue risque de comporter des zones (ou pixels) non exploitables, qui devront donc par exemple être représentés sous forme de pixels noirs. Contrairement à la technique précitée de moyennage des données ultrasonores, la présente divulgation repose sur une combinaison de célérités individuelles VI permettant d’en déduire une célérité composée VC. Ainsi, au lieu de moyenner des données ultrasonores représentant des réponses ultrasonores pour différentes valeurs VL d’un paramètre PR des ondes de sondage 14, la divulgation prévoit par exemple de moyenner des célérités individuelles VI estimées individuellement pour respectivement différentes valeurs VL d’un paramètre PR des ondes de sondage 14. Ainsi, en générant des célérités individuelles VI associées respectivement à différentes valeurs VL d’un paramètre PR que l’on fait varier, on peut avantageusement obtenir de multiples estimations de célérités des ondes de cisaillement 12 pour une position donnée du milieu 4. Il a certes été constaté que les besoins en termes de ressources de traitement peuvent être plus élevés que lorsque la technique de moyennage des données ultrasonores est mise en œuvre. En revanche, le principe de la divulgation permet avantageusement d’améliorer considérablement la qualité des estimations ainsi obtenues. La divulgation entraîne une augmentation de la taille des données exploitables, toutes choses étant égales par ailleurs, ce qui peut certes augmenter les besoins en ressources mais permet en contrepartie une meilleure précision et fiabilité de l’estimation. Aussi, en équipant un dispositif de traitement 20 avec des ressources de traitement adaptées (par exemple dans le système ultrasonore SY1), on peut avantageusement atteindre un compromis pertinent entre robustesse de l’estimation de célérité et coût de traitement, notamment en temps et moyens de traitement (mémoire, moyens de calcul…). Un tel compromis est possible et adaptable au regard notamment des développements actuels et futurs des capacités de traitement des systèmes ultrasonores. La robustesse des estimations de célérité peut avantageusement être maximisée en pondérant l’importance des célérités individuelles VI dans le calcul de la célérité composée VC en fonction de la pertinence ou qualité des célérités individuelles VI obtenues respectivement pour chaque valeur VL du paramètre PR considéré. On peut ainsi affecter des poids K différents selon la valeur VL du paramètre PR, voire exclure au moins une célérité individuelle VI (en fixant par exemple son poids K à zéro) pour ne retenir dans le calcul de la célérité composée VC que certaines parmi les célérités individuelles VI. Par une analyse des ondes de cisaillement 12, on peut donc estimer de façon précise et fiable la vitesse de ces ondes et en déduire la rigidité ou élasticité du milieu, par exemple dans des tissus d’un corps humain ou animal (par exemple dans des tissus mous, tels que tissus adipeux, ou encore fibres musculaires). Comme le comprend l’homme du métier, tous les modes de réalisation et variantes décrits ci-avant dont certains ont été simplifiés à dessein pour faciliter les explications, ne constituent que des exemples non limitatifs de mise en œuvre de la présente divulgation. En particulier, l’homme du métier pourra envisager une quelconque adaptation ou combinaison des modes de réalisation et variantes décrits ci-avant, afin de répondre à un besoin particulier. La présente divulgation ne se limite donc pas aux exemples de réalisation décrits ci-avant mais s’étend notamment à un procédé de traitement qui inclurait des étapes secondaires sans pour cela sortir de la portée de la présente divulgation. Il en serait de même d’un système de traitement pour la mise en œuvre d’un tel procédé. Procédé de traitement d’élastographie ultrasonore comprenant : a) obtention (S2) de premières données (DT1) représentatives d’une réponse ultrasonore (16) d’un milieu d’intérêt (4), à des impulsions ultrasonores (8) focalisées dans ledit milieu en combinaison avec des ondes de sondage (14) émises dans ledit milieu séquentiellement selon des valeurs (VL) distinctes d’au moins un paramètre (PR) desdites ondes de sondage ; b) génération (S4), à partir des premières données respectivement pour chaque valeur (VL) dudit au moins un paramètre, de deuxièmes données (DT2) représentatives d’une propagation d’ondes de cisaillement (12) dans le milieu d’intérêt sous l’effet des impulsions ultrasonores ; c) estimation (S6), à partir des deuxièmes données associées à chaque valeur (VL) dudit au moins un paramètre, d’une célérité individuelle (VI) représentative de la célérité de déplacement des ondes de cisaillement pour respectivement ladite valeur ; et d) estimation (S8) d’au moins une célérité (VC), dite célérité composée, des ondes de cisaillement par une combinaison des célérités individuelles. Procédé selon la revendication 1, dans lequel ledit au moins un paramètre des ondes de sondage comprend au moins l’un parmi : - un angle d’émission (θ) selon lequel les ondes de sondage sont émises en direction du milieu d’intérêt ; et - une fréquence (f) selon laquelle sont émises les ondes de sondage en direction du milieu d’intérêt. Procédé selon la revendication 1 ou 2, dans lequel les ondes de sondage (14) comprennent des ondes planes ultrasonores pour sonder le milieu d’intérêt. Procédé selon l’une quelconque des revendications précédentes, dans lequel le procédé comprend, préalablement à l’obtention a) : e1) formation des ondes de cisaillement dans le milieu d’intérêt par focalisation des impulsions ultrasonores ; et e2) émission, dans le milieu d’intérêt, des ondes de sondage configurées séquentiellement selon les valeurs distinctes dudit au moins un paramètre pour détecter la réponse ultrasonore du milieu d’intérêt pour chaque valeur. Procédé selon l’une quelconque des revendications précédentes, dans lequel les deuxièmes données générées en b) définissent, pour chaque valeur dudit au moins un paramètre, une séquence d’images (IM1) représentative de la propagation des ondes de cisaillement causées par les impulsions ultrasonores dans le milieu d’intérêt. Procédé selon l’une quelconque des revendications précédentes, dans lequel au cours de l’estimation d), la combinaison des célérités individuelles est déterminée par un calcul d’une moyenne des célérités individuelles (VI1, Vl2, VI3) pour estimer ladite au moins une célérité composée (VC). Procédé selon la revendication 6, dans lequel l’estimation d) comprend : - attribution de poids (K1, K2, K3), en association avec les célérités individuelles (VI1, VI2, VI3), en fonction d’au moins une caractéristique des deuxièmes données générées respectivement pour chaque valeur (VL1, VI2, VL3) dudit au moins un paramètre ; et - estimation de ladite au moins une célérité composée par calcul de la moyenne en pondérant les célérités individuelles selon lesdits poids associés. Procédé selon la revendication 7, dans lequel l’estimation d) comprend : - détermination de variances spatiales caractérisant respectivement les deuxièmes données associées à chaque valeur dudit au moins un paramètre ; et - détermination, à partir des variances spatiales, de scores de confiance associés respectivement aux célérités individuelles, le poids attribué à chaque célérité individuelle étant fonction du score de confiance associé. Procédé selon la revendication 7 ou 8, dans lequel l’estimation d) comprend : - si un score de confiance associée à des deuxièmes données correspondant à une valeur dudit au moins un paramètre ne satisfait pas à un au moins critère de qualité, le poids attribué à la célérité individuelle associée à ladite valeur est nul de sorte à exclure ladite célérité individuelle de l’estimation d) de ladite au moins une célérité composée. Procédé selon l’une quelconque des revendications précédentes, dans lequel le procédé comprend en outre : f) estimation, à partir de ladite au moins une célérité composée, d’une dureté du milieu d’intérêt. Procédé selon l’une quelconque des revendications précédentes, dans lequel au cours de l’estimation d), une pluralité de célérités composées (VC) sont estimées en différentes positions du milieu d’intérêt, l’estimation d) comprenant en outre : - génération, à partir des célérités composées, de données cartographiques (DT3) représentatives des célérités composées, et/ou de duretés estimées à partir desdites célérités composées, dans respectivement les différentes positions du milieu d’intérêt. Procédé selon l’une quelconque des revendications précédentes, dans lequel le procédé comprend en outre : - affichage des célérités composées, et/ou de duretés estimées à partir desdites célérités composées, sous forme d’au moins une image ou carte (IM2). Procédé selon l’une quelconque des revendications précédentes, dans lequel le procédé comprend, après une phase préliminaire (PH1) au cours de laquelle sont exécutées les étapes a) à d) pour chaque dite valeur distinctes (VL) dudit au moins un paramètre, au moins une itération périodique d’un cycle (S20) comprenant : - estimation d’au moins une célérité individuelle additionnelle par réitération des étapes a) à c) pour au moins l’une parmi les dites valeurs dudit au moins un paramètre ; et - estimation d’au moins une célérité composée des ondes de cisaillement par réitération de l’étape d) en utilisant la célérité individuelle additionnelle et au moins une célérité individuelle estimée au cours d’une itération précédente du cycle et/ou au cours de la phase préliminaire. Programme d’ordinateur (PG1) comportant des instructions pour l’exécution des étapes d’un procédé de configuration selon l’une quelconque des revendications précédentes lorsque ledit programme est exécuté par un ordinateur. Dispositif (20) de traitement d’élastographie ultrasonore comprenant : - un module d’obtention (MD2) configuré pour obtenir des premières données (DT1) représentatives d’une réponse ultrasonore (16), d’un milieu d’intérêt (4), à des impulsions ultrasonores (8) focalisées dans ledit milieu en combinaison avec des ondes de sondage (14) émises dans ledit milieu séquentiellement selon des valeurs distinctes (VL) d’au moins un paramètre (PR) desdites ondes de sondage ; - un module de génération (MD4) configuré pour générer, à partir des premières données respectivement pour chaque valeur dudit au moins un paramètre, de deuxièmes données (DT2) représentatives d’une propagation d’ondes de cisaillement (12) dans le milieu d’intérêt sous l’effet des impulsions ultrasonores ; - un premier module d’estimation (MD6) configuré pour estimer, à partir des deuxièmes données associées à chaque valeur dudit au moins un paramètre, une célérité individuelle (VI) représentative de la célérité de déplacement des ondes de cisaillement pour respectivement ladite valeur ; et - un deuxième module d’estimation (MD8) configuré pour estimer au moins une célérité (VC), dite célérité composée, des ondes de cisaillement par une combinaison des célérités individuelles. Système de traitement (SY2) comprenant : - le dispositif de traitement (20) selon l’une quelconque des revendications précédentes ; et - au moins une sonde ultrasonore (2) configurée pour focaliser des impulsions ultrasonores dans ledit milieu et/ou émettre les ondes de sondage ultrasonores dans ledit milieu.