Brevet FR3140259A1 — publié le 20240405

L’invention concerne un procédé d’analyse d’une région d’intérêt d’un milieu viscoélastique, le procédé incluant une étape d’acquisition (200) de signaux représentatifs du déplacement d’au moins une composante propagative de cisaillement et d’une composante de déformation statique, et une étape de traitement (300) pour déterminer au moins une propriété du milieu, l’étape de traitement comprenant les sous-étapes suivantes : détermination d’images de la région d’intérêt, estimation d’images de déplacement par comparaison des images de la région d’intérêt, filtrage des images de déplacement pour atténuer la composante de déformation statique, détermination d’une propriété du milieu à partir des images de déplacement filtrées. Figure à publier avec l’abrégé : Fig. 6 PROCEDE D’ANALYSE D’UN MILIEU PERMETTANT DE REDUIRE LES EFFETS D’ARTEFACTS DUS A DES DEFORMATION STATIQUES DANS LE MILIEU DOMAINE DE L'INVENTION La présente invention concerne le domaine technique général de l’imagerie d’un objet cible, ou d’un milieu diffus tel qu’un tissu biologique humain ou animal. Plus précisément, la présente invention concerne un dispositif et un procédé pour mesurer des propriétés viscoélastiques d’un tissu biologique d’intérêt. Elle s'applique en particulier, mais non exclusivement, à la mesure de paramètres de viscoélasticité du foie d'un humain ou d'un animal, cette mesure étant corrélée à la quantité de fibrose présente dans le foie. ARRIERE PLAN DE L'INVENTION 1. Généra lités Afin de mesurer les propriétés viscoélastiques tissulaires, il est connu d'utiliser l'élastographie par onde de cisaillement. Cette technique consiste à observer le déplacement d'une onde de cisaillement à travers une région d'intérêt (ROI) d’un milieu pour déterminer des paramètres de l’onde de cisaillement dans la ROI, telles que la vitesse de propagation de l’onde de cisaillement, ces paramètres étant directement reliés aux propriétés du milieu analysé (e.g. le module de cisaillement du milieu analysé). Pour ce faire, des ondes de cisaillement sont générées dans une ROI, par exemple en mettant en œuvre une sollicitation mécanique, ou une sollicitation acoustique. Les ondes de cisaillement générées provoquent des déplacements variant dans le temps tout en se déplaçant d’un point de génération à de multiples emplacements dans le milieu analysé. Différentes techniques d’imagerie comme l’imagerie par résonance magnétique (ou « magnetic resonance elastography » en anglais), l’imagerie optique (ou « optical coherence elastography » en anglais) ou l’imagerie ultrasonore (ou « ultrasound elastography » en anglais) peuvent détecter ces déplacements. Dans le cas particulier de l’imagerie ultrasonore, une sonde d’imagerie est utilisée. Une telle sonde d’imagerie comprend un boîtier et un réseau de transducteurs, mécaniquement solidaire du boîtier, pour l’émission d’ondes ultrasonores haute fréquence et la réception des échos acoustiques permettant de suivre le déplacement des ondes de cisaillement. Notamment, le suivi des déplacements induits par les ondes de cisaillement en fonction du temps aux multiples emplacements permet une estimation de la vitesse de cisaillement, qui à son tour permet d’estimer une (ou plusieurs) propriété(s) du milieu analysé. La caractérisation de propriétés viscoélastiques du milieu telles que la rigidité au cisaillement à l'aide de l'estimation de la vitesse de cisaillement a d'importantes applications médicales car ces propriétés sont étroitement liées à l'état du milieu par rapport à une pathologie. Typiquement, au moins une partie d'un tissu peut devenir plus rigide que les tissus environnants indiquant un début ou la présence d'une maladie ou d'un état tel qu'un cancer, une tumeur, une fibrose. Cependant, les techniques d'imagerie de vitesse de cisaillement conventionnelles souffrent d'artefacts dus à la présence de déformations statiques – engendrées par un mouvement de la sonde et/ou par un mouvement d’un organe du patient etc. – et qui dégradent la qualité des signaux acquis par la sonde, ce qui rend difficile la mesure précise des propriétés viscoélastiques du milieu. Un but de la présente invention est de proposer un procédé d’analyse d’une région d'intérêt (ROI) d’un milieu viscoélastique permettant de réduire les effets associés aux artefacts (dus à la présence de déformations statiques) dans la détermination d’une (ou de plusieurs) propriété(s) du milieu viscoélastique. Avant de présenter l’invention, une description des principes fondamentaux associés à l’élastographie par onde de cisaillement va être détaillée afin de mieux comprendre les notions relatives au problème de déformation statique évoqué ci-dessus. 2. Théorie relative à l’élastographie par onde de cisaillement 2.1. Propagation d’ondes en milieu élastique Des principes fondamentaux de la propagation d’ondes en milieu élastique vont maintenant être rappelés. Pour cela, considérons un milieu solide élastique homogène et isotrope. L’équation d’onde en l’absence de force s’écrit de la manière suivante : où : le vecteur désigne le déplacement en tout point, le symbole désigne l’opérateur différentiel vectoriel, le symbole désigne l’opérateur Laplacien, est la densité du milieu, et et sont le premier coefficient de Lamé et le module de cisaillement du solide, respectivement. L’équation d’onde impose la propagation de deux types de composantes propagatives dans le milieu : la composante propagative de compression (ondes P) et la composante propagative de cisaillement (ondes S). Ces composantes propagatives, dans une approximation monochromatique (onde plane) ont des comportements distincts. La composante propagative de compression est à divergence nulle, polarisée longitudinalement et se propage à une vitesse c p telle que : La composante propagative de cisaillement est à rotationnel nul, polarisée transversalement et se propage à une vitesse c s telle que : Ainsi, ces deux composantes propagatives peuvent se propager à des vitesses différentes, selon les valeurs relatives des deux coefficients de Lamé. Dans un tissu biologique, qui est quasi incompressible ( ), les vitesses c p et c s des composantes propagatives de compression et de cisaillement ont des valeurs très différentes, typiquement autour de 1500 m/s pour l’onde compression et autour de 1 à 5 m/s pour l’onde de cisaillement. De plus, les tissus biologiques étant quasi-incompressibles, la composante propagative de compression est d’amplitude négligeable par rapport à la composante propagative de cisaillement et nous pouvons considérer que les déplacements observés dans le cadre des techniques d’élastographie sont uniquement liés à des effets de cisaillement. On comprend donc que la composante propagative de compression est celle exploitée par l’imagerie ultrasonore classique (permettant d’atteindre une résolution spatiale acceptable avec des fréquences de l’ordre du MHz) utilisée pour imager les déplacements dans le cadre de l’invention, et la composante propagative de cisaillement est celle sur laquelle se base l’estimation des paramètres viscoélastiques tissulaires dans les méthodes d’élastographie. 2.2. Elastographies harmonique, transitoire et passive En référence à la , le principe de l’élastographie harmonique se base sur l'utilisation d'une source mécanique pour générer – dans un milieu – une onde mécanique incluant une composante propagative de cisaillement CIS, CIS’ (et une composante propagative de compression CMP, CMP’ dont l’amplitude est négligeable dans les tissus biologiques) à basse fréquence dont la propagation est suivie par imagerie échographique ultrarapide. En particulier la illustre la génération d’une composante propagative de cisaillement CIS et d’une composante propagative de compression CMP lors d’une compression du milieu par une sonde S. Le suivi de la propagation de la composante propagative de cisaillement au cours du temps nous permet de retrouver sa vitesse de propagation et le module de Young du milieu. L’idée sous-jacente est basée sur un modèle rhéologique simple. En effet, si l’on suppose que le milieu est quasi incompressible (approximation valide dans le cas de tissus biologiques), alors le module de cisaillement (et donc la vitesse de l’onde de cisaillement) est proportionnel au module de Young du milieu comme suit : Ainsi, la mesure de la vitesse de la composante propagative de cisaillement permet d’estimer le module de Young du milieu. Le Doppler tissulaire ultrasonore, allié à une échographie ultrarapide pour éviter tout problème d’aliasing, permet de suivre la composante propagative de cisaillement et ainsi d’en mesurer sa vitesse pour retrouver le module de Young. Différentes techniques de génération de la composante propagative de cisaillement ont été explorées dans l’art antérieur. Une première technique utilise une source harmonique externe au réseau de transducteurs pour générer la composante propagative de cisaillement. Le problème principal de cette technique réside dans l’utilisation d’une source externe de vibration qui augmente significativement la complexité de mise en œuvre de la mesure. Une deuxième technique – dite d’élastographie transitoire – popularisée par la méthode “Vibration Controlled Transient Elastography” (VCTE™) implémentée dans le système Fibroscan® (Echosens, Paris, France), résout ce problème en générant la composante propagative de cisaillement à l’aide de la vibration d’un élément piézo-électrique, qui est ensuite utilisé pour suivre la propagation de la composante propagative de cisaillement. La force du dispositif est que, par des effets de symétrie, la composante propagative de cisaillement polarisée longitudinalement, i.e. le terme de champ proche, est isolée des autres et captée par traitement Doppler le long de la ligne enregistrée par l’élément piézoélectrique. Ainsi, par simple traitement Doppler opéré le long de la ligne enregistrée par l’élément piezoéléctrique, la vitesse de la composante propagative de cisaillement est mesurée. Cependant, cette deuxième technique d’élastographie développée dans le Fibroscan ne permet pas de faire une élastographie bi-dimensionnelle et n’est pas adaptée pour capter des changements de dureté au sein du milieu. Une troisième technique, dite ARFI ou Supersonic Shear Wave selon les implémentations, consiste à générer une onde de cisaillement à l'intérieur du tissu lui même en utilisant un faisceau ultrasonore focalisé de forte intensité qui génère par effet non-linéaire une force mécanique dite ‘de radiation’. Le même transducteur ultrasonore est utilisé pour générer l’onde de cisaillement et pour en mesurer le déplacement ultérieurement par acquisition ultrarapide. Cette technique a l’avantage de ne pas générer par elle-même de déformation statique du milieu, mais n'empêche pas les déformations statiques d’élastographie passive (par exemple celles issues des battements cardiaques) décrites ci-après d’interférer avec l’onde générée. Une quatrième technique – dite d’élastographie passive par onde naturelle – consiste à exploiter les ondes de cisaillement induites naturellement par le corps, par exemple par le cœur sur le foie lors d’un cycle cardiaque. La technique originelle se base sur une hypothèse de champ de déplacement diffus : les composantes propagatives de cisaillement se propagent dans toutes les directions. Dans le cadre de cette hypothèse, la dérivée temporelle de l’intercorrélation spatiale du champ de déplacement en un point est liée à des fonctions de Green causale et anti-causale en ce point. Les propriétés spatio-temporelles de ces fonctions de Green peuvent être utilisées pour estimer le module de Young du milieu. Dans ce contexte, la présence de déformations statiques induit un biais dans la mesure de l’intercorrélation et donc dans celle de la dureté. Aucune solution pour la réduction des effets associés aux déformations statiques n’a été développée pour cette troisième technique, bien qu’une telle solution soit nécessaire pour des applications d’élastographie passive du foie où le cœur induit de fortes déformations statiques. 2.3. Problème des déformations statiques Le Doppler tissulaire permet de mesurer le déplacement en chaque point du milieu, en fonction du temps. Dans le cas des élastographies passive et transitoire, ce déplacement est une somme de trois composantes : les deux composantes propagatives de compression et de cisaillement décrites en 2.1; une composante de déformation statique décrite ci-dessous. La composante de déformation statique est induite par une contrainte sur le milieu provoquée : soit par la sonde, soit par un générateur externe de l’onde de cisaillement tel qu’un vibreur externe, soit par un organe. On entend, dans le cadre de la présente invention par « composante de déformation statique » ou « composante de déformation quasi-statique », une déformation du milieu induite par l’application d’une contrainte – par exemple depuis la surface du milieu – et révélatrice de phénomènes de cisaillement non propagatifs. Cette composante de déformation statique étant différente à chaque instant, elle s’ajoute aux composantes propagatives dans le déplacement mesuré. Elle est donc problématique pour toutes les méthodes basées sur la mesure des composantes propagatives, e.g. les techniques d’élastographie passive ou transitoire. Dans les tissus biologiques, cette composante est forte à proximité de la source de la contrainte, et nulle infiniment loin de la source de la contrainte. Une telle composante de déformation statique peut être due aux mouvements de la sonde, et/ou aux déplacements d’un organe – tel que le cœur à chaque battement cardiaque ou les poumons – etc. Dans la suite, le problème de déformation statique sera présenté en référence à une sonde incluant : un boitier, un réseau de transducteurs logé dans le boitier et configuré pour : émettre des ondes ultrasonores d’excitation vers le milieu à analyser (organe, tissu biologique, etc.), et recevoir des échos acoustiques (i.e. ondes ultrasonores réfléchies par les différentes interfaces du milieu à analyser), et un excitateur logé dans le boitier et configuré pour : générer des vibrations, et les transmettre à la sonde afin de produire mécaniquement une onde de cisaillement dont l’observation permet la mesure de propriétés du milieu à analyser. On a illustré à la une sonde S au repos disposée sur le milieu M à analyser. Lorsque l’excitateur est activé pour générer et transmettre des vibrations à la sonde S (pour générer la composante propagative de cisaillement), ces vibrations entraînent des contraintes successives d’appui ( ) et de relâchement ( ) de la sonde S sur le milieu engendrant des mouvements locaux et des déformations statiques. Comme représenté aux figures 2b et 2c, dans le référentiel de la sonde S, le mouvement relatif à la sonde d’un point (tel que le point A) situé à la surface du milieu est nul (par continuité), tandis que l’amplitude du mouvement relatif à la sonde d’un point (tel que le point B, ou le point C) situé sous la surface du milieu augmente avec la profondeur, et tend vers une asymptote correspondant à l’amplitude du mouvement de la sonde pour les points les plus profonds (tel que le point D). Afin de mieux comprendre ce phénomène, il est possible d’utiliser un modèle physique du milieu tel qu’illustré à la et qui consiste en un empilement de ressorts dont les coefficients de raideur CR1, CR2, CR3 sont croissants en fonction de la profondeur, du fait des conditions aux limites du problème. Dans un milieu incompressible, le coefficient de raideur de ces ressorts est lié aux contraintes de cisaillement et dépend du module d’Young du milieu, ainsi qu’à la géométrie du problème. A partir du modèle physique illustré à la figure 3, une étude du déplacement des points B, C, situés à différentes profondeurs dans le milieu peut être effectuée, lorsqu’une force F appliquée par la sonde induit un déplacement d’un point A de la surface du milieu. Par un bilan de force, il est possible de déduire que les compressions et aux points B et C s’expriment par les formules suivantes : On a donc , avec et . Ainsi, le déplacement – dû à l’application de la force proche du point A – mesuré en chaque point A, B, C du milieu varie en fonction de la profondeur. A titre illustratif, on a représenté à la un profil de la composante de déformation statique (i.e. déplacement du milieu dû aux mouvements de la sonde et relatif à la sonde) en fonction de la profondeur. Un tel profil présente globalement une forme logarithmique tendant vers une asymptote. Comme illustré à la , cette composante de déformation statique DS (i.e. déplacement du milieu dû aux mouvements de la sonde) vient se superposer aux composantes propagatives CP (i.e. déplacement du milieu dû en particulier à la composante propagative de cisaillement de la , la composante propagative de compression étant négligeable), ce qui rend difficile l’observation de la composante propagative de cisaillement seule, et donc la mesure de propriétés du milieu à analyser. En effet, l'artefact de la (dû à la composante de déformation statique engendrée par les mouvements de la sonde et/ou par les mouvements d’un organe du patient etc.) est particulièrement gênant du fait que : la fréquence temporelle des déplacements induits par cette composante de déformation statique DS est similaire à celle des déplacements induits par la composante propagative de cisaillement, ce qui rend sa suppression par filtrage purement temporel difficile; l’amplitude des déplacements induits par cette composante de déformation statique DS : est significativement plus élevée que celle des déplacements induits par la composante propagative de cisaillement sous-jacente (un à deux ordres de grandeurs), et augmente avec la profondeur, ce qui est problématique car l’amplitude des déplacements induits par la composante propagative de cisaillement diminue avec la profondeur (par des effets de dispersion et d’absorption). 3. Présentation de l’état de la technique Différentes méthodes ont été utilisées afin de s’abstraire de cet artefact. a) Une première méthode, mise en place dans le système Fibroscan®, consiste à utiliser une excitation large bande, ou de manière équivalente, à retenir la vibration d’une source externe pour la génération de l’onde de cisaillement. Plus précisément, le système Fibroscan® propose d’atténuer l’oscillation de la source externe afin de générer l’onde de cisaillement sur un court instant, comme décrit dans le document WO 2018/078002. Cette atténuation (ou « damping », selon la terminologie anglo-saxonne) permet de faciliter la séparation des déplacements du milieu dus à la composante de déformation statique (qui est ramenée à zéro après un cycle d’oscillation), des déplacements du milieu dus à la propagation de la composante propagative de cisaillement. Ainsi, la composante propagative de cisaillement et la composante de déformation statique sont découplées temporellement. Le cisaillement est observé lorsque la composante de déformation statique est très faible dans le milieu. Le problème majeur de cette méthode est sa faible efficacité énergétique puisque beaucoup d’énergie est gaspillée dans la génération d’une vibration large bande (i.e. pour exciter et amortir la source externe de génération de la composante propagative de cisaillement). b) Une deuxième méthode, décrite dans un document – intitulé “Probe Oscillation Shear Elastography (PROSE): A High Frame-Rate Method for Two-Dimensional Ultrasound Shear Wave Elastography”, au nom de D. C. Mellema et al. – consiste à utiliser la stroboscopie afin d’échantillonner systématiquement la même réalisation de la comosante de déformation statique. Ainsi, la composante de déformation statique ne biaise plus les mesures de la composante propagative de cisaillement. Le problème de cette méthode est qu’elle impose une fréquence d’échantillonnage temporelle (inférieure à deux fois la fréquence de vibration de la source externe) trop faible pour analyser avec efficacité la propagation de la composante propagative de cisaillement en basse fréquence, par exemple à 50Hz. De plus, cette méthode s’est avérée assez peu efficace in-vivo dû à son incapacité à supprimer l’intégralité de la composante de déformation statique, comme décrit dans le document intitulé “Probe Oscillation Shear Wave Elastography: Initial In Vivo Results in Liver” par D. C. Mellema et al. c) Enfin, de nombreuses méthodes basées sur des techniques de filtrage ont été explorées afin de réduire les mouvements parasites (dû aux composantes propagatives de compression mais surtout à celles de cisaillement dans des directions non voulues) de celui dû à la composante propagative de cisaillement d’intérêt, notamment en élastographie par résonance magnétique. Ces méthodes impliquent notamment des filtrages unidimensionnels ou bi-dimensionnels. Un exemple de telles méthodes est décrit dans le document intitulé “Spatio-temporal directional filtering for improved inversion of MR elastography images” par Manduca et al. Un inconvénient de ces méthodes est qu’elles nécessitent l’emploi de transformée de Fourier bi- ou tri-dimensionnelles, ce qui peut s’avérer coûteux computationnellement parlant et poser des problèmes de performance aux limites des films. Par ailleurs, certaines de ces méthodes se basent sur des filtrages directionnels non-linéaires (seuillage dans le domaine fréquentiel) et peuvent engendrer d’importants artefacts, e.g. une propagation d’onde fantôme dans un cas où il n’y a pas de propagation. d) Inspirées de ces méthodes, des techniques de filtrage utilisant un filtre à réponse impulsionnelle finie (FIR) passe bande ont été testées sans succès du fait de leur manque de robustesse dans la séparation de la composante de déformation statique et de la composante propagative de cisaillement. Ces essais infructueux sont décrits dans le document intitulé “Probe Oscillation Shear Wave Elastography: Initial In Vivo Results in Liver” par D.C. Mellema et al. qui enseigne à l’homme du métier que l’utilisation d’un filtre à réponse impulsionnelle finie n’est pas adaptée pour le filtrage de la composante de déformation statique. Des techniques de filtrage utilisant un filtre à réponse impulsionnelle infinie (IIR) ont également été envisagées. Toutefois, les filtre à réponse impulsionnelle infinie (IIR) engendrent des distorsions de phase et ont des effets transitoires, générant d’importants artefacts dans la mesure de vitesse de propagation de la composante propagative de cisaillement. Un but de la présente invention est de proposer un procédé d’analyse d’une région d’intérêt d’un milieu permettant de pallier au moins l’un des inconvénients précités. BREVE DESCRIPTION DE L'INVENTION A cet effet, l’invention propose un procédé d’analyse d’une région d’intérêt d’un milieu viscoélastique, le procédé incluant une phase d’acquisition de signaux représentatifs du déplacement d’au moins une composante propagative de cisaillement et d’une composante de déformation statique, la phase d’acquisition comprenant les étapes suivantes : émettre successivement dans le temps des signaux d’acquisition dans le milieu viscoélastique, détecter et enregistrer des signaux de réception reçus du milieu viscoélastique, chaque signal de réception étant associé à un signal d’acquisition respectif, une phase de traitement pour déterminer au moins une propriété du milieu, ladite phase de traitement comprenant les sous-étapes suivantes : détermination d’images de la région d’intérêt du milieu à partir des signaux de réception, chaque image étant déterminée à partir d’un signal de réception respectif et comprenant des informations représentatives de la région d’intérêt à des instants différents, estimation d’images de déplacement par comparaison des images de la région d’intérêt, chaque image de déplacement comportant des informations représentatives de la composante propagative de cisaillement et d’une composante de déformation statique, filtrage des images de déplacement pour atténuer la composante de déformation statique, détermination de ladite et au moins une propriété du milieu à partir des images de déplacement filtrées. On entend, dans le cadre de la présente invention par « onde mécanique », une onde comprenant : une composante propagative de cisaillement générée soit naturellement (i.e. élastographie passive dans laquelle la composante propagative de cisaillement est générée par un organe tel que le cœur lors d’un cycle cardiaque), soit artificiellement (i.e. élastographie harmonique ou transitoire dans laquelle la composante propagative de cisaillement est générée par la force de radiation ultrasonore, ou par une source externe – tel qu’un excitateur de vibration inertiel du type décrit dans le document WO 2022/084502, et une composante propagative de compression, générée elle aussi soit naturellement soit artificiellement, et dont l’amplitude dans les tissus biologiques est faible relativement à la composante propagative de cisaillement. L’amplitude de cette composante propagative de compression étant très faible relativement à l’amplitude de la composante propagative de cisaillement, celle-ci sera négligée dans la suite de la description. On entend, dans le cadre de la présente invention, par « onde de cisaillement » (« shear wave » selon la terminologie anglosaxonne) ou « composante propagative de cisaillement », une onde « basse fréquence » polarisée transversalement en champ lointain, c’est-à-dire une onde produisant des mouvements du milieu selon une direction perpendiculaire à la direction de propagation de l’onde et dont la fréquence est comprise typiquement entre 10Hz et 1500Hz. On entend, dans le cadre de la présente invention, par « onde de compression » (« compressional wave » selon la terminologie anglosaxonne) ou « composante propagative de compression », une onde « basse fréquence » polarisée longitudinalement en champ lointain, c’est-à-dire une onde produisant des mouvements du milieu selon une direction parallèle à la direction de propagation de l’onde et dont la fréquence est comprise typiquement entre 10Hz et 1500Hz. On entend, dans le cadre de la présente invention, par « composante de déformation statique », les déformations induites par une contrainte sur le milieu. Dans les milieux biologiques, cette composante de déformation statique est liée aux phénomènes de cisaillement. On entend, dans le cadre de la présente invention, par « onde ultrasonore », une onde de compression « haute fréquence » polarisée longitudinalement, c’est-à-dire une onde produisant des mouvements du milieu selon une direction parallèle à la direction de propagation de l’onde et dont la fréquence est comprise typiquement entre 15kHz et 20MHz. Dans le cadre de l’invention, cette onde ultrasonore est utilisée comme moyen d’acquérir les images de déplacement utilisées pour suivre la propagation de la composante propagative de cisaillement. On entend, dans le cadre de la présente invention, par « vitesse de groupe » d’une onde, une vitesse à laquelle une forme globale de l'enveloppe des amplitudes de l'onde, connue sous le nom d’enveloppe de l'onde, se propage dans l'espace. Des aspects préférés mais non limitatifs du procédé d’analyse selon l’invention sont les suivants : l’étape de filtrage peut comprendre l’application d’un filtre à réponse impulsionnelle finie dont les caractéristiques sont définies en fonction : d’une fréquence moyenne de la composante propagative de cisaillement comprise entre 10Hz et 1500Hz, et d’une vitesse maximale de groupe de la composante propagative de cisaillement dans la région d’intérêt ; le filtre à réponse impulsionnelle finie peut comprendre une bande passante et une bande atténuée : la bande passante étant définie au moins par une borne inférieure K-6, ladite borne inférieure K-6 étant égale à un rapport entre la fréquence moyenne de la composante de cisaillement et la vitesse maximale de groupe de la composante de cisaillement, l’atténuation dans la bande passante étant inférieure ou égale à 6 dB, la bande atténuée étant définie au moins par une borne d’arrêt K-30 proportionnelle à la borne inférieure K-6, l’atténuation dans la bande atténuée étant supérieure ou égale à 30dB, et le rapport entre la borne inférieure K-6 et la borne d’arrêt K-30 définissant un coefficient de raideur r du filtre à réponse impulsionnelle finie, ledit coefficient de raideur étant compris entre 1 et 4.5 ; les vitesses minimale et maximale peuvent être déterminées en fonction de duretés minimale et maximale prédéfinies dans la région d’intérêt ; la fréquence moyenne de la composante de cisaillement peut être estimée à partir d’un signal représentatif d’un spectre temporel des images de déplacement ; le procédé peut comprendre en outre une phase d’excitation incluant une étape consistant à générer la composante propagative de cisaillement en appliquant au milieu viscoélastique une excitation ayant la forme d'une impulsion basse fréquence qui présente une fréquence centrale f comprise entre 10Hz et 1500Hz ; le filtre à réponse impulsionnelle finie peut être de type passe-haut dans une gamme de fréquence utile comprise entre 0 et 50 m⁻1 ; ledit filtre à réponse impulsionnelle finie peut être obtenu par combinaison d’un filtre de type passe-tout et d’un filtre de type passe-bas dans ladite gamme de fréquence utile ; le filtre de type passe-bas peut être un filtre passe-bas à fenêtre rectangulaire ou un filtre dont la réponse impulsionnelle est telle que sa valeur maximale est inférieure ou égale à 2 fois sa valeur moyenne. D'autres avantages et caractéristiques du procédé d’analyse selon l’invention ressortiront mieux de la description qui va suivre de plusieurs variantes d’exécution, données à titre d'exemples non limitatifs, à partir des dessins annexés sur lesquels : est une représentation schématique illustrant la génération d’ondes de compression et de cisaillement lors d’une compression d’un milieu, est une représentation schématique d’une sonde au repos positionnée sur un milieu à analyser, est une représentation schématique illustrant la génération de déformations statiques par appui de la sonde sur le milieu à analyser, est une représentation schématique illustrant la génération de déformations statiques par relâchement de la sonde sur le milieu à analyser, est une représentation schématique d’un modèle de ressorts, est une représentation schématique d’un profil de déplacement induit par la composante de déformation statique en fonction d’une direction de propagation de la composante propagative de cisaillement, est une représentation schématique d’un profil de déplacement induit par la superposition de la composante de déformation statique et de la composante propagative de cisaillement, en fonction d’une direction de propagation de la composante propagative de cisaillement, est une représentation schématique d’un profil de déplacement induit par la composante propagative de cisaillement en fonction d’une direction de propagation de la composante propagative de cisaillement, une fois la composante de déformation statique, indiquée DS dans la , atténuée en utilisant la méthode décrite selon l’invention, est une représentation schématique d’un dispositif d’analyse selon l’invention, est une représentation schématique d’une variante de réalisation du procédé d’analyse selon l’invention, est une représentation schématique d’une réponse impulsionnelle d’un exemple de filtre pouvant être utilisé lors d’une sous-étape de filtrage du procédé d’analyse. DESCRIPTION DETAILLEE DE L'INVENTION On va maintenant décrire plus en détails différents modes de réalisation du procédé d’analyse selon l’invention en référence aux figures. Dans ces différentes figures, les éléments équivalents sont désignés par la même référence numérique. Dans la suite, le procédé d’analyse sera décrit en référence au traitement de données acquises en utilisant une sonde d'élastographie ultrasonore permettant : de générer au moins une onde élastique basse fréquence – dite « onde de cisaillement » – dans le milieu, et simultanément à la génération de l’onde basse fréquence : émettre des ondes ultrasonores haute fréquence, et recevoir des échos acoustiques dus aux réflexions des ondes ultrasonores dans le milieu, afin d’observer la propagation de l’onde élastique basse fréquence dans le milieu. Il est bien évident pour l’homme du métier que le procédé d’analyse peut être mis en œuvre avec des techniques d’élastographie passive par onde naturelle – dans lesquelles les ondes de cisaillement sont générées naturellement par le corps. Dans ce cas, la sonde comprend uniquement un réseau de transducteurs pour l’émission d’ondes ultrasonores en fonction d’une profondeur haute fréquence (15kHz - 20MHz) et la réception d’échos acoustiques afin d’observer la propagation d’une (ou plusieurs) onde(s) de cisaillement générée(s) naturellement par le corps. Il est également bien évident pour l’homme du métier que le procédé d’analyse selon l’invention peut être mis en œuvre avec d'autres techniques d'imageries que l'imagerie ultrasonore ("ultrasound elastography" en anglais), notamment l’imagerie par résonance magnétique ("magnetic resonance elastography" en anglais), ou l’imagerie optique ("optical coherence elastography" en anglais). 1. Généralités En référence à la , on a illustré un exemple de dispositif dans lequel le procédé d’analyse d’un milieu décrit dans la suite peut être mis en œuvre. Ce dispositif comprend : une sonde S d’acquisition de signaux, et une unité de commande et de traitement U c pour : le pilotage de la sonde S, et le traitement des signaux acquis par la sonde S. 1. 1 . Sonde d'élastographie impulsionnelle La sonde S comprend : un réseau de transducteurs pour l’émission des ondes ultrasonores et la réception des échos acoustiques, et un excitateur pour la génération de l’onde de cisaillement. L’excitateur peut être un excitateur de vibration inertielle du type décrit dans le document WO 2022/084502. Il permet de générer des vibrations. Une partie fixe de l’excitateur est mécaniquement solidaire du réseau de transducteurs pour transmettre les vibrations à la sonde afin de produire mécaniquement l’onde de cisaillement nécessaire à la mesure des propriétés viscoélastiques du tissu à analyser. L’utilisation d’un excitateur de vibration inertiel permet d’obtenir une sonde dont le poids et l’encombrement sont limités. L’utilisation d’un excitateur de vibration inertiel permet en outre de générer une onde de cisaillement de puissance suffisante pour mesurer les propriétés viscoélastiques tissulaires tout en consommant une énergie inférieure à d’autres solutions (telles que les solutions basées sur la pression de radiation ultrasonore). L’excitateur peut aussi être externe à la sonde, et la mouvoir dans sa totalité, bien que cette solution ne soit pas optimale du point de vue de l’encombrement. Une telle sonde étant connue de l’état de la technique, elle ne sera pas décrite plus en détails dans la suite. 1. 2. Unité de commande et de traitement L’unité de commande et de traitement U c est reliée à la sonde S par des moyens de communication avec ou sans fil. Elle permet de piloter les transducteurs de la sonde S, et de traiter les données acquises par les transducteurs de la sonde S. Elle permet en outre d’activer l’excitateur de vibration inertiel pour la génération d’une (ou plusieurs) onde(s) de cisaillement dans le milieu à analyser qui peut être un organe d’un patient tel que le foie. Plus précisément, l’unité de commande et de traitement Uc permet : de commander à l’excitateur de vibration inertiel une vibration permettant la génération d’une onde de cisaillement dans le milieu à analyser, de commander aux transducteurs l'émission d'ondes ultrasonores vers le milieu à analyser, de commander aux transducteurs la réception des échos réfléchis par le milieu à analyser et leur conversion en signaux de réception, de traiter les signaux de réception. L’unité de commande et de traitement U c peut être composée d’une ou de plusieurs entités physiques distinctes, éventuellement distantes de la sonde S. L’unité de commande et de traitement Uc comprend par exemple : un (ou plusieurs) contrôleur(s) 11, tel qu’un Smartphone, un assistant personnel (ou « PDA », sigle de l’expression anglo-saxonne « Personal Digital Assistant »), ou tout type de terminal mobile connu de l’homme du métier ; et un (ou plusieurs) calculateur(s) 12, tel qu’un/des ordinateur(s), un/des micro-ordinateur(s), une/des stations de travail, et/ou d'autres dispositifs connus de l’homme du métier incluant un/des processeur(s), un/des microcontrôleur(s), un/des automate(s) programmable(s), un/des circuit(s) intégré(s) spécifique(s) d'application, et/ou d'autres circuits programmables, une (ou plusieurs) unité(s) de stockage 13 comportant une (ou plusieurs) mémoire(s) qui peu(ven)t être une mémoire ROM/RAM, une clé USB, une mémoire d’un serveur central. Outre la conservation de données associées à l’analyse d’un milieu, l’unité de stockage 13 permet également de stocker des instructions de code de programmation destinées à exécuter les étapes du procédé d’analyse décrit dans la suite. 1. 3 . Principe de fonctionnement En référence aux figures 5 et 6, le principe de fonctionnement du dispositif d’imagerie est le suivant. Le contrôleur 11 fait générer des vibrations dans le milieu viscoélastique par l’excitateur, au cours d'une étape d'excitation (étape 100). Ceci permet de générer une onde mécanique dans le milieu, incluant une composante propagative de cisaillement et une composante propagative de compression négligeable dans les tissus biologiques. Les déplacements du milieu induits par la propagation l’onde mécanique dans le milieu viscoélastique sont mesurés par échographie, au cours d'une étape d'observation (étape 200) concomitante à l'étape d'excitation (100). A cet effet, le contrôleur 11 pilote l’émission d’ondes ultrasonores (à une fréquence comprise par exemple entre 15kHz et 20MHz) par les transducteurs de la sonde S. Ces ondes ultrasonores pénètrent dans le milieu où elles se réfléchissent sur des particules diffusantes contenues dans ledit milieu – comme par exemple des particules de collagène contenues dans le milieu s’il s’agit d’un tissu humain – ce qui permet de suivre les mouvements du milieu. Les ondes ultrasonores ainsi émises peuvent être « planes » (c'est à dire des ondes dont le front d'onde est rectiligne dans un plan d'analyse (D, Z)) ou spirales, ou divergentes, ou encore focalisées. En particulier, l’étape d'observation (200) comprend les sous-étapes suivantes : le contrôleur 11 commande l’émission par le réseau de transducteurs, d’une succession de tirs d'ondes ultrasonores dans le milieu (à une cadence comprise entre 10 et 10000 tirs par seconde), le contrôleur 11 commande la détection par le réseau de transducteurs, et l’enregistrement (en temps réel) dans l’unité de stockage 13 des signaux acoustiques reçus comprenant les échos générés par les ondes ultrasonores en interagissant avec les particules diffusantes du milieu. Les signaux acoustiques enregistrés dans l’unité de stockage 13 sont ensuite traités (généralement en temps différé) par le calculateur lors d’une étape de traitement 300. Plus précisément, l’étape de traitement 300 comprend une sous-étape de détermination – par un processus classique de formation de voies (ou « beamforming » selon la terminologie anglo-saxonne) – d’images d’une région d’intérêt du milieu contenue dans le champ d’observation (i.e. la zone insonifiée par les ondes ultrasonores). Chaque image est représentative de la région d’intérêt après un tir respectif. Comme les tirs sont réalisés successivement dans le temps, chaque image est représentative de la région d’intérêt à un instant différent. Ces images successives de la région d’intérêt sont ensuite comparées, par corrélation et en particulier par intercorrélation : soit deux à deux soit avec une image de référence (qui peut être une image de déplacement précédemment déterminée ou une image de la région d’intérêt obtenue lors d’une étape préliminaire d’observation), pour déterminer des images de déplacement du milieu. Ces images de déplacement sont représentatives de la propagation de l’onde mécanique dans le milieu. Ces images de déplacement comprennent également une composante de déformation statique, qui peut significativement gêner le suivi de la propagation de l’onde mécanique, utilisée pour la détermination d’une (ou plusieurs) propriété(s) du milieu. C’est pourquoi les inventeurs ont développé une sous-étape de filtrage des signaux acoustiques reçus – et en particulier des images de déplacement – pour supprimer la composante de déformation statique des déplacements mesurés, et ainsi ne conserver que la composante propagative, i.e. l’onde mécanique, dans les images de déplacement. Cette sous-étape de filtrage va maintenant être décrite plus en détails. 2. Sous-étape de filtrage La sous-étape de filtrage comprend l’application d’un filtre à réponse impulsionnelle finie dont les caractéristiques sont définies en fonction : d’une fréquence moyenne de la composante propagative de cisaillement comprise typiquement entre 10Hz et 1500Hz, et d’une vitesse maximale (notamment d’une vitesse maximale de groupe) de la composante propagative de cisaillement dans la région d’intérêt. La fréquence moyenne de la composante propagative de cisaillement dépend de la solution retenue pour la génération de cette composante propagative de cisaillement. Toutefois, dans tous les cas, cette fréquence moyenne est connue ou peut être mesurée. Par exemple, dans le cas de l’utilisation d’un excitateur, celui-ci est configuré pour générer une composante propagative de cisaillement à une fréquence moyenne souhaitée (par exemple à une fréquence moyenne de 50 Hz). Dans ce cas, la fréquence moyenne de la composante propagative cisaillement est donc connue à priori. Si la solution retenue pour la génération de l’onde de cisaillement est naturelle (génération de l’onde de cisaillement par un organe du patient), alors la fréquence moyenne de la composante propagative de cisaillement peut être calculée à partir des images de déplacement par toute technique connue de l’homme du métier. De même, la vitesse de groupe maximale de la composante propagative de cisaillement dans la région d’intérêt dépend du type de région d’intérêt observée. Toutefois, pour chaque type de région d’intérêt, cette vitesse de groupe maximale est connue dans la littérature ou peut être mesurée. Par exemple, si la région d’intérêt consiste en un tissu mou (foie, rate, prostate, sein, muscle etc.), il est connu que la vitesse de groupe maximale de la composante propagative de cisaillement est de l’ordre de 6m/s. Ainsi, les valeurs de la fréquence moyenne et de la vitesse de groupe maximale de la composante propagative de cisaillement sont connues et permettent le dimensionnement du filtre à réponse impulsionnelle finie qui va maintenant être présenté ci-dessous. 2. 1. Caractéristiques du filtre En référence à la , on a illustré la réponse fréquentielle (i.e. variation du gain du filtre en fonction de la fréquence du signal qui lui est appliqué en entrée) d’un exemple de filtre pouvant être utilisé lors de la sous-étape de filtrage. Cette réponse fréquentielle est une courbe représentative d’une atténuation (en dB) en fonction d’une fréquence (en Hz). Dans ce mode de réalisation, le filtre à réponse impulsionnelle finie est un filtre « passe-haut ». Comme illustré à la , cette réponse fréquentielle (ou tracé de gain d’un diagramme de Bode) est composée : d’une bande atténuée BA, d’une bande passante BP, et d’une bande de transition BT entre les bandes atténuée et passante BA, BP. 2 . 1.1 . Bande atténuée La bande atténuée BA correspond à une plage de fréquences dans laquelle un signal appliqué en entrée du filtre apparaît fortement diminué en sortie du filtre : dans cette plage de fréquences, l’amplitude du signal appliqué en entrée du filtre est atténuée (à une valeur proche de zéro) en sortie du filtre. Dans le cas d’un filtre passe-haut, la bande atténuée BA peut être définie par une borne d’arrêt. La borne d’arrêt – ou nombre d’onde d’arrêt – correspond à un rapport entre une fréquence et une vitesse de groupe (représentative de la vitesse de propagation de la composante propagative de cisaillement) exprimé en cycles par mètre. Cette borne d’arrêt permet d’exprimer vers quelle fréquence le filtre cesse de produire son effet « atténué ». Dans le cadre de la présente invention, cet « effet atténué » est considéré produit lorsque l’amplitude d’un signal appliqué en entrée du filtre est divisée par 32 (ou plus) en sortie du filtre (-30dB). Nous utiliserons donc le terme K -30 dans la suite du texte pour caractériser la partie haute de la bande atténuée. Ainsi, la borne d’arrêt K -30 (ou nombre d’onde d’arrêt K -30 ) de la bande atténuée BA est représentative de la fréquence à laquelle l’amplitude d’un signal appliqué en entrée du filtre est divisée par 32 en sortie du filtre (-30dB). Dans cette bande atténuée BA, l’amplitude d’un signal appliqué en entrée du filtre est donc atténuée d’un facteur d’atténuation FA2 supérieur ou égal à 30dB. 2. 1.2. Bande passante La bande passante BP correspond à une plage de fréquences dans laquelle un signal appliqué en entrée du filtre apparaît peu ou pas atténué en sortie du filtre : dans cette plage de fréquences, l’amplitude du signal appliqué en entrée du filtre est sensiblement inchangée (ou faiblement atténuée) en sortie du filtre. Dans le cas d’un filtre passe-haut, la bande passante BP peut être définie par une borne inférieure. Cette borne inférieure – ou nombre d’onde – correspond à un rapport entre une fréquence et une vitesse de groupe (représentative de la vitesse de propagation de la composante propagative de cisaillement) exprimé en cycles par mètre. Cette borne inférieure permet d’exprimer vers quelle fréquence le filtre cesse à produire son effet « passant ». Dans le cadre de la présente invention, cet « effet passant » est considéré produit lorsque l’amplitude d’un signal appliqué en entrée du filtre est divisée par 2 (ou moins) en sortie du filtre (-6dB). Nous utiliserons donc le terme K -6 dans la suite du texte pour caractériser la partie basse de la bande passante. Ainsi, la borne inférieure K -6 de la bande passante BP est représentative de la fréquence à laquelle l’amplitude d’un signal appliqué en entrée du filtre est divisée par 2 en sortie du filtre (-6dB). Dans cette bande passante BP, l’amplitude d’un signal appliqué en entrée du filtre est donc atténuée d’un facteur d’atténuation FA 1 inférieur ou égal à 6dB. La borne inférieure K -6 correspond au rapport entre : la fréquence moyenne de la composante propagative de cisaillement, et la vitesse maximale de groupe de la composante propagative de cisaillement. Comme indiqué précédemment, la fréquence moyenne de la composante propagative de cisaillement est connue. De même, la vitesse de groupe maximale de la composante propagative de cisaillement est connue (de l’ordre de 6 m/s dans le cas d’un tissu mou). Il est donc possible de définir la borne inférieure K -6 de la bande passante BP en fonction du type de tissu analysé et du type d’excitation envisagée pour la génération de la composante propagative de cisaillement. 2.1.3. Bande de transition La bande de transition BT correspond à une plage de fréquences dans laquelle un signal appliqué en entrée n’est ni fortement atténué (bande atténuée BA) ni sensiblement conservé inchangé (bande passante BP). Cette bande de transition BT caractérise la capacité du filtre à conserver ou rejeter plus ou moins abruptement les fréquences. Elle est donc représentative de la « réactivité » du filtre. Dans cette bande de transition, nous caractérisons le filtre par un coefficient de raideur (ou « constante de raideur ») correspondant au rapport entre : la borne inférieure K -6 de la bande passante BP, et la borne d’arrêt K -30 de la bande atténuée BA. Ainsi, le coefficient de raideur se définit comme le rapport suivant : Ce coefficient de raideur illustre la pente de la réponse en amplitude du filtre dans la bande de transition entre un régime où le filtre atténue fortement et un régime où le filtre n’atténue plus. Les inventeurs ont déterminé qu’il était préférable que le coefficient de raideur soit inférieur ou égal à 4.5 pour que le filtre permette une atténuation significative de la composante de déformation statique dans les signaux acoustiques reçus sans impacter de façon inacceptable les fréquences basses de cisaillement correspondant à des vitesses élevées des ondes de cisaillement, et donc des duretés élevées. L’efficacité du filtre dépend de deux critères de la bande de transition : une fréquence de coupure, le coefficient de raideur. La fréquence de coupure du filtre correspond à la fréquence limite de fonctionnement utile du filtre. Plus précisément, la fréquence de coupure permet d’exprimer vers quelle fréquence le filtre commence à produire son effet « passant ». Ainsi, la borne inférieure K -6 de la bande passante BP est représentative de la fréquence de coupure du filtre. 2.2. Variantes du filtre D’autres types de filtres peuvent être utilisés pour atténuer la composante de déformation statique contenue dans les signaux reçus par la sonde. Par exemple, un filtre de type passe-bande peut être utilisé à la place du filtre passe-haut illustré à la . Dans le cas d’un filtre passe-bande, la bande passante BP est en outre définie par une borne supérieure K max égale à un rapport entre la fréquence moyenne de la composante propagative de cisaillement et une vitesse minimale de groupe de la composante propagative de cisaillement. Les valeurs de la fréquence moyenne et de la vitesse de groupe minimale étant connue (mesurable ou connue de la littérature en fonction du tissu à imager), cette borne supérieure K max est aisément déterminable. Dans tous les cas, le filtre choisi présente le comportement d’un filtre de type passe-haut dans une gamme de fréquence utile, par exemple comprise entre entre 0 et 50 m ⁻ 1 . Avantageusement, ce filtre à réponse impulsionnelle finie peut être obtenu par combinaison d’un filtre de type « passe-tout » et d’un filtre de type « passe-bas » dans ladite gamme de fréquence utile. Ceci permet de faciliter l’obtention d’un filtre efficace pour : l'atténuation de la composante de déformation statique, et la conservation de la composante propagative de cisaillement. 2. 3 . Théorie relative à l’invention Une technique de filtrage le long de la direction de la composante de déformation statique est proposée, afin de drastiquement réduire l’impact de la composante de déformation statique. Cette technique s’appuie sur l’idée que la composante de déformation statique et la composante propagative de cisaillement ont des caractéristiques spatiales très différentes : faibles longueurs d’onde (quelques cms au maximum à 50Hz) pour la composante propagative de cisaillement avec des vitesses de groupe variant entre 0.5 m/s et 6 m/s dans les tissus mous comme la rate et le foie. hautes longueurs d’onde (dizaines de m à 50Hz) pour la composante de déformation statique. 2. 3. 1 . Filtre FIR passe-haut Dans un mode de réalisation, la technique de filtrage se base sur un filtre à réponse impulsionnelle finie (FIR) passe-haut (pas nécessairement à réponse de phase linéaire). Ses caractéristiques fréquentielles sont définies en fonction de la fréquence moyenne des ondes de cisaillement étudiées f ainsi que des vitesses de groupe minimales et maximales d’intérêt c min et c max . La fréquence moyenne de l’onde de cisaillement peut être estimée en utilisant la formule ci-dessous : Où : est un signal représentatif du spectre temporel de l’onde de cisaillement (e.g. la transformée de Fourier temporelle des images de déplacement). Des valeurs typiques de fréquence moyenne des ondes de cisaillement sont comprises entre 10Hz et 500Hz. La valeur de fréquence moyenne est le plus souvent supposée connue. Par exemple, dans le cas où l’onde de cisaillement est générée naturellement (élastographie passive), la valeur de la fréquence moyenne peut être mesurée. Dans le cas où l’onde de cisaillement est générée artificiellement (élastographie transitoire) la valeur de la fréquence moyenne est déduite des caractéristiques de la source mécanique utilisée. La vitesse de groupe d’une onde de cisaillement dont les fréquences se situent dans l’intervalle est estimée comme suit : Des valeurs typiques de vitesse de groupe des ondes de cisaillement sont comprises entre 0.5 m/s et 5.5 m/s pour des tissus mous tels que le foie. Nous nous intéressons en particulier à deux valeurs du nombre d’onde : K -6 qui correspond à la valeur du nombre d’onde (en nombre de cycles par mètre) pour laquelle l’atténuation du filtre est égale à 6 dB, K -30 qui correspond à la valeur du nombre d’onde (en nombre de cycles par mètre) pour laquelle l’atténuation du filtre est égale à 30 dB. A titre indicatif, on rappelle que le « nombre d’onde » (également connu sous le nom de « répétence », ou « repetency » en anglais) est une grandeur proportionnelle à l'inverse de la longueur d'onde. Le « nombre d’onde » est l’analogue, dans l’espace, de la fréquence temporelle et caractérise la propagation d’une onde dans l’espace, comme le fait la fréquence temporelle dans le temps. L’atténuation mentionnée ci-dessus est calculée en fonction de la valeur du gain du filtre dans la bande passante comme suit : où : H est l’amplitude de la réponse fréquentielle du filtre d’intérêt et G H est une quantité représentative du gain du filtre dans la bande passante. Cette quantité peut être par exemple le maximum de la valeur absolue de la réponse fréquentielle dans la bande passante ou bien la valeur moyenne de la valeur absolue de la réponse fréquentielle dans la bande passante. Associées à ces deux valeurs, nous définissons le coefficient de raideur de la réponse du filtre passe-haut comme ceci : Ce coefficient de raideur caractérise la pente moyenne de la réponse fréquentielle du filtre entre ces deux valeurs puisque nous avons : A plus haut niveau, il peut être vu comme quantificateur de l’efficacité du filtre passe-haut. Ce filtre est très générique et est réalisable par l’homme du métier en utilisant des outils de synthèse de filtre à réponse impulsionnelle finie, comme les fonctions « firwin » et « firwin2 » de la librairie « scipy signal processing » du langage de programmation Python. Avantageusement, le coefficient de raideur du filtre est choisi afin d’obtenir une efficience suffisante du filtre. Par exemple, la valeur du coefficient de raideur peut être choisi inférieur à 4.5. Un filtre passe-haut de ce type est réalisable par l’homme du métier en construisant un filtre FIR à distorsion de phase. Un exemple de code est donné ci-dessous pour le cas d’une onde de cisaillement dont la fréquence de groupe est supposée à 50Hz. Nous définissons un filtre à réponse impulsionnelle finie de 2.5 cm de long. # Fréquence d'échantillonnage Dz = 0.6523e-3 # Fréquence de l'onde de cisaillement f = 50 # Hz # Définition du filtre passe-haut probe_fir_len = 2.5e-2 probe_half_fir = int(np.round(probe_fir_len / Dz - 1) / 2) probe_fir = -np.ones((2 * probe_half_fir +1, )) / (2 * probe_half_fir + 1) probe_fir[probe_half_fir] = 0 probe_fir[probe_half_fir] = 1 - 1./len(probe_fir) # Réponse fréquentielle du filtre w, h_probe_fir = sig.freqz(probe_fir, fs=1/Dz, worN=8192) # Fréquence de coupure et coefficient de raideur associés cutoff_0 = np.argmin(np.abs(20*np.log10(np.abs(h_probe_fir) / np.abs(h_probe_fir[-1])) + 30)) cutoff_1 = np.argmin(np.abs(20*np.log10(np.abs(h_probe_fir) / np.abs(h_probe_fir[-1])) + 6)) k_6 = w[cutoff_1] k_30 = w[cutoff_0] c_max = f / k_30 c_min = f / k_6 r = k_6 / k_30 Dans ce cas, les grandeurs caractéristiques du filtre sont les suivantes : K -6 = 24.6 m -1 , K -30 = 5.71 m -1 et r = 4.31. Le lecteur appréciera qu’une contrainte de raideur inférieure à 4.5 est difficilement atteignable en utilisant la fonction « firwin2 » puisque la fonction « firwin2 » impose une contrainte de linéarité de phase (ou de manière équivalente d’antisymétrie du filtre) pour assurer que le filtre ne présente pas de distorsion de phase. Un exemple ci-dessous illustre cette difficulté dans un cas similaire à ci-dessus : probe_fir_len = 2.5e-2 probe_half_fir = int(np.round(probe_fir_len / Dz - 1) / 2) numtaps = 2 * probe_half_fir cutoff_1 = k_min cutoff_0 = k_min / 16 hp = sig.firwin2(numtaps, [0.0, cutoff_0, cutoff_1, 1.0 / 2 / Dz], [0, 10**(-30 / 20), 1.0, 1.0], fs = 1 / Dz, antisymmetric=True, window='rect') w, h = sig.freqz(hp, fs=1/Dz, worN=2048*4) cutoff_0_meas = np.argmin(np.abs(20*np.log10(np.abs(h) / np.abs(h[-1])) + 30)) cutoff_1_meas = np.argmin(np.abs(20*np.log10(np.abs(h) / np.abs(h[-1])) + 6)) k_6 = w[cutoff_1_meas] k_30 = w[cutoff_0_meas] c_max = f / k_30 c_min = f / k_6 r = k_6 / k_30 slope = 1. / alpha Dans ce cas, le coefficient de raideur associé est à 16.7. Pour réduire le coefficient de raideur, il est possible de faire varier la longueur du filtre ou d’autres paramètres. Ainsi, l’exemple ci-dessous permet d’atteindre une raideur inférieure à 4.5 en utilisant “firwin” et un filtre de taille 3.2 cm. # Fréquence d'échantillonnage Dz = 0.6523e-3 # Caractéristiques de l'onde de cisaillement f = 50 # Hz c_min = np.sqrt(12 / 3) # m/s c_max = np.sqrt(500 / 3) # m/s k_max = f / c_min # m^-1 k_min = f / c_max # m^-1 # Définition du filtre passe-haut probe_fir_len = 3.2e-2 probe_half_fir = int(np.round(probe_fir_len / Dz - 1) / 2) numtaps = 2 * probe_half_fir hp = sig.firwin(numtaps+1, k_max, width=2*(k_max-k_min), fs = 1 / Dz, pass_zero='highpass') # Réponse fréquentielle du filtre w, h = sig.freqz(hp, fs=1/Dz, worN=2048*4) # Nombres d'onde k_6 et k_30 et coefficient de raideur associé cutoff_0_meas = np.argmin(np.abs(20*np.log10(np.abs(h) / np.abs(h[-1])) + 30)) cutoff_1_meas = np.argmin(np.abs(20*np.log10(np.abs(h) / np.abs(h[-1])) + 6)) k_6 = w[cutoff_1_meas] k_30 = w[cutoff_0_meas] c_max = f / k_30 c_min = f / k_6 r = k_6 / k_30 Dans ce cas, les grandeurs caractéristiques du filtre sont les suivantes : K -6 = 25.1 m -1 , K -30 = 5.6 m -1 et r = 4.48. On peut ainsi générer un filtre selon la spécification de raideur, mais la longueur nécessaire à la réalisation du filtre : 3.2 cm, sera perdue pour la visualisation des vitesses des ondes de cisaillement (raccourcissement de la profondeur de la zone d'intérêt). Bien entendu, d’autres fonctions Matlab ou python (que la fonction « fiwin2 ») peuvent être utilisées. 2.3.2. Filtre FIR passe-bande Dans un mode de réalisation de l’invention, la technique de filtrage se base sur un filtre passe-bande unidimensionnel appliqué le long de la déformation. De manière similaire au filtre passe haut défini au point 2.4.1, nous nous intéressons aux quantités suivantes : K -6 qui correspond à la valeur du nombre d’onde (en nombre de cycles par mètre) pour laquelle l’atténuation du filtre est égale à 6 dB dans la bande de transition inférieure, le coefficient de raideur du filtre. Il est à noter que le coefficient de raideur correspondant sur lequel nous imposons une contrainte similaire à celle du filtre passe-haut (i.e. coefficient de raideur inférieur ou égal à 4.5) est défini dans ce cas sur la bande de transition inférieure du filtre passe-bande. Aucune contrainte n’est imposée sur le coefficient de raideur dans la bande de transition supérieure. La longueur de sa réponse impulsionnelle ne doit pas dépasser quelques cm, une taille typique étant comprise entre 1cm et 3cm, afin de limiter les effets aux limites tout en maintenant une région suffisamment grande dans le milieu d’intérêt. Le lecteur aura compris que de nombreuses modifications peuvent être apportées à l’invention décrite précédemment sans sortir matériellement des nouveaux enseignements et des avantages décrits ici. Par conséquent, toutes les modifications de ce type sont destinées à être incorporées à l’intérieur de la portée des revendications jointes.   Procédé d’analyse d’une région d’intérêt d’un milieu viscoélastique, le procédé incluant une phase d’acquisition (200) de signaux représentatifs du déplacement d’au moins une composante propagative de cisaillement et d’une composante de déformation statique, la phase d’acquisition (200) comprenant les étapes suivantes : émettre successivement dans le temps des signaux d’acquisition dans le milieu viscoélastique, détecter et enregistrer des signaux de réception reçus du milieu viscoélastique, chaque signal de réception étant associé à un signal d’acquisition respectif, caractérisé en ce que le procédé comprend en outre une phase de traitement (300) pour déterminer au moins une propriété du milieu, ladite phase de traitement comprenant les sous-étapes suivantes : détermination d’images de la région d’intérêt du milieu à partir des signaux de réception, chaque image étant déterminée à partir d’un signal de réception respectif et comprenant des informations représentatives de la région d’intérêt à des instants différents, estimation d’images de déplacement par comparaison des images de la région d’intérêt, chaque image de déplacement comportant des informations représentatives de la composante propagative de cisaillement et d’une composante de déformation statique, filtrage des images de déplacement pour atténuer la composante de déformation statique, détermination de ladite et au moins une propriété du milieu à partir des images de déplacement filtrées. Procédé d’analyse selon la revendication 1, dans lequel l’étape de filtrage comprend l’application d’un filtre à réponse impulsionnelle finie dont les caractéristiques sont définies en fonction : d’une fréquence moyenne de la composante propagative de cisaillement comprise entre 10Hz et 1500Hz, et d’une vitesse maximale de groupe de la composante propagative de cisaillement dans la région d’intérêt. Procédé d’analyse selon la revendication 2, dans lequel le filtre à réponse impulsionnelle finie comprend une bande passante (BP) et une bande atténuée (BA) : la bande passante (BP) étant définie au moins par une borne inférieure K -6 , ladite borne inférieure K -6 étant égale à un rapport entre la fréquence moyenne de la composante de cisaillement et la vitesse maximale de groupe de la composante de cisaillement, l’atténuation dans la bande passante étant inférieure ou égale à 6 dB, la bande atténuée (BA) étant définie au moins par une borne d’arrêt K -30 proportionnelle à la borne inférieure K -6 , l’atténuation dans la bande atténuée étant supérieure ou égale à 30dB, et le rapport entre la borne inférieure K -6 et la borne d’arrêt K -30 définissant un coefficient de raideur r du filtre à réponse impulsionnelle finie, ledit coefficient de raideur étant compris entre 1 et 4.5. Procédé d’analyse selon l’une quelconque des revendications 2 ou 3, dans lequel les vitesses minimale et maximale sont déterminées en fonction de duretés minimale et maximale prédéfinies dans la région d’intérêt. Procédé selon l’une quelconque des revendications 2 à 4, dans lequel la fréquence moyenne de la composante de cisaillement est estimée à partir d’un signal représentatif d’un spectre temporel des images de déplacement. Procédé selon l’une des revendications 2 à 4, lequel comprend en outre une phase d’excitation (100) incluant une étape consistant à générer la composante propagative de cisaillement en appliquant au milieu viscoélastique une excitation ayant la forme d'une impulsion basse fréquence qui présente une fréquence centrale f comprise entre 10Hz et 1500Hz. Procédé selon l’une quelconque des revendications 2 à 6, dans lequel le filtre à réponse impulsionnelle finie est de type passe-haut dans une gamme de fréquence utile comprise entre 0 et 50 m⁻1. Procédé selon la revendication 7, dans lequel ledit filtre à réponse impulsionnelle finie est obtenu par combinaison d’un filtre de type passe-tout et d’un filtre de type passe-bas dans ladite gamme de fréquence utile. Procédé selon la revendication 8, dans lequel le filtre de type passe-bas est un filtre passe-bas à fenêtre rectangulaire ou un filtre dont la réponse impulsionnelle est telle que sa valeur maximale est inférieure ou égale à 2 fois sa valeur moyenne.