Brevet FR2662265A1 — publié le 19911122

"Dispositif éliminateur d'échos fixes pour échographe ultraso- nore". Description La présente invention concerne un dispositif éli- minateur d'échos fixes pour échographe ultrasonore à émission d'impulsions récurrentes muni d'un circuit d'intercorrélation et d'interpolation, en vue de la détermination de profils de vitesse de milieux en mouvement, comportant au moins un filtre éliminateur d'échos fixes du type passe-bande à lignes à re- tard. Une application particulièrement avantageuse de ce dispositif consiste en la mesure de profils de vitesse d'écou- lements sanguins par échographie ultrasonore. Le problème technique général à résoudre par tout échographe ultrasonore destiné à la mesure et, complémentaire- ment, à l'imagerie des vitesses d'écoulements sanguins est d'éliminer les échos fixes dus à la grande réflectivité des tissus biologiques (ici, les parois des vaisseaux) qui peut dépasser de 40 d B la réflectivité du sang (globules rouges). L'utilisation d'un dispositif d'élimination d'échos fixes est donc indispensable, avant toute estimation du profil de vites- se, pour réduire la dynamique des signaux à traiter et pour mesurer correctement les vitesses d'écoulement sur toute une section des vaisseaux, en particulier les faibles vitesses qui apparaissent au voisinage des parois de vaisseaux et qui sont, de ce fait, masquées par les forts signaux provenant de ces parois. Une solution simple à ce problème technique géné- ral consiste par exemple, dans un dispositif éliminateur d'échos fixes pour échographe ultrasonore, à introduire un filtre, du type à lignes à retard, composé d'une ligne retar- dée d'une période de récurrence en parallèle sur une ligne de retard nul Des coefficients de pondération, respectivement + 1 et -1, sont affectés à ces lignes qui, après pondération, sont additionnées par un additionneur Ce filtre connu réalise donc la différence entre deux lignes échographiques consécutives ce qui, par principe, conduit à une diminution quasi complète des échos produits par les tissus fixes Toutefois, cette techni- que présente l'inconvénient majeur d'atténuer également les signaux correspondant aux faibles vitesses d'écoulement A ti- tre d'exemple, on peut montrer que la réponse du filtre ci- dessus en fonction de la vitesse d'écoulement est telle que, pour une fréquence de récurrence de 1/T = 5 k Hz, et une fré- quence d'émission de Nc = 5 M Hz, un signal correspondant à V = 5 cm/s est atténué de 30 d B. Ceci rend difficile, voire impossible, la mesure des vitesses d'écoulement là o elles sont les plus faibles, c'est-à-dire près des parois des vaisseaux Or, la connaissan- ce de ces vitesses est très importante pour l'étude et le dia- gnostic clinique des artères, par exemple. Il est connu, par ailleurs, que l'on peut augmen- ter la pente du filtre mentionné plus haut, dans la zone de coupure en augmentant son ordre Appliqué à la mesure des vi- tesses par échographie, ceci revient à traiter non plus 2 li- gnes échographiques successives mais davantage, soit M, à l'aide d'un filtre du type à lignes à retard, composé de M li- gnes parallèles i de retard respectif (i-1)T De façon à assu- rer l'élimination des échos fixes, les M lignes sont pondérées par des coefficients dont la somme est nulle Puis toutes les lignes sont additionnées pour ensuite être traitées par un circuit estimateur de vitesse, par corrélation, décalage de phase, etc Cependant, il est difficile avec ce type de fil- tre d'obtenir à la fois une bonne mesure du profil de vitesse de l'écoulement étudié et une cadence d'image élevée Si N est le nombre de lignes à traiter (typiquement N = 16 pour les vaisseaux périphériques et N = 8 pour le coeur), et si M li- gnes sont utilisées pour le fitlre d'élimination des échos fi- xe, il ne reste en réalité que N-M+ 1 lignes utilisables pour l'estimation de vitesse, ou alors, pour avoir le même nombre de lignes à traiter, il faudrait diminuer la cadence d'image. On connaît aussi de la demande de brevet français 2 617 982 au nom de la demanderesse un dispositif éliminateur d'échos fixes pour échographe ultrasonore comportant, associés en parallèle, au moins deux filtres de même structure, conju- gués et déphasés entre eux de (M-1)T, chaque filtre étant com- posé, d'une part, de M lignes parallèles i (i = 1,, M) de retard respectif (i-1)T, T étant la période de récurrence de l'échographe, d'autre part, de moyens de pondération des M li- gnes i, la somme des coefficients affectés auxdites lignes étant nulle, et, enfin, d'un additionneur des lignes i ainsi retardées et pondérées, dispositif gâce auquel on puisse obte- nir une estimation satisfaisante des faibles vitesses d'écou- lements sanguins, notamment à proximtié des parois des vais- seaux, sans diminution importante des lignes à traiter, donc sans réduction de la cadence d'image On obtient ainsi deux filtres d'ordre M-1 ayant en module la même reponse en fré- quence On dispose alors en sortie de ces filtres de: 2 (N-M+ 1) signaux pour estimer les vitesses au lieu de N-M+ 1. Ces signaux ne sont pas tout à fait indépendants, mais avec un estimateur de vitesse utilisant une opération non linéaire disposé en aval (corrélation 1 bit, décalage de phase,), on obtient ainsi une amélioration appréciable Cependant, une corrélation plus fine que la corrélation 1 bit atténuerait les effets de l'amélioration ainsi obtenue Le problème technique précis que se propose de résoudre la présente invention con- siste, tout en conservant un nombre élevé de lignes utilisa- bles pour l'estimation de vitesse, à obtenir un filtrage, pour la réjection des échos fixes, qui soit correct pour les vites- ses faibles à mesurer tout en restant proche de l'optimum pour les vitesses plus élevées. Le but principal de l'invention est de réaliser un dispositif éliminateur d'échos fixes qui effectue un filtrage correct pour toute la plage de vitesses utiles à mesurer typi- quement comprises entre O et 100 cm/s. Ce but est atteint et les inconvénients de l'art antérieur sont atténués ou supprimés grâce au fait que le dis- positif éliminateur d'échos fixes défini en préambule est re- marquable en ce qu'il comporte une pluralité de filtres élimi- nateurs d'échos fixes du type précité en parallèle, dont les bandes passantes sont juxtaposées, des détecteurs d'amplitude reliés chacun à la sortie de chacun desdits filtres et à l'en- trée d'un comparateur, des éléments de retard branchés entre la sortie de chaque filtre et un multiplexeur, ce dernier re- cevant un signal de commande issu dudit comparateur pour sé- lectionner en sortie dudit multiplexeur celui parmi les si- gnaux de sortie desdits filtres qui présente la plus grande amplitude. On obtient ainsi un multifiltrage, chaque filtre étant adapté à une gamme de vitesses particulière Le choix consistant à sélectionner à travers le multiplexeur le filtre qui fournit les données ayant la plus grande amplitude se jus- tifie par le fait qu'on cherche à obtenir une amplitude, pour le signal utile débarrassé des échos fixes, qui se situe net- tement au-dessus du niveau de bruit, pour les vitesses faibles particulièrement De plus, l'utilisation de plusieurs filtres dont les bandes passantes sont juxtaposées, permet de sélec- tionner un filtre dévolu aux vitesses les plus faibles dont la pente est forte, tout en conservant pour ce filtre (et aussi pour les autres filtres en parallèle) un ordre peu élevé En ce qui concerne l'ordre N des filtres on observe expérimenta- lement qu'il n'est pas possible d'améliorer de façon signifi- cative le dispositif suppresseur d'échos fixes selon l'inven- tion lorsqu'on utilise des filtres à 3 coefficients Par ail- leurs, les résultats obtenus avec des filtres à 5 coefficients sont à peu près semblables à ceux obtenus avec des filtres à 4 coefficients quoiqu'au prix d'une plus grande complexité et aussi, selon certains modes de réalisation basés sur le multifiltrage, au risque accru de discontinuités apparaissant dans les profils de vitesses mesurés, raison pour laquelle ces modes de réalisation, dans lesquels les fonctions de suppress- sion d'échos fixes et d'intercorrélation sont totalement dis- sociées, ne sont pas retenus De préférence, les filtres éli- minateurs d'échos fixes selon l'invention sont des filtres à 4 coefficients, à somme des coefficients nulle, cette dernière caractéristique étant la condition première, connue, pour la suppression d'échos fixes, soit: n E ak = ( 1) k=O Selon un mode de réalisation préféré de l'inven- tion, les filtres suppresseurs d'échos fixes sont conçus et sélectionnés pour rendre maximale la fonction A( 0) du rapport signal à bruit entre la sortie et l'entrée desdits filtres: n -jko |k O E ak e A( 0) = k= O ( 2) d E ak 2 k=O avec: 0 = 4 T Nc VT/C o: Nc: fréquence du signal échographique pendant l'impulsion, V:vitesse du flux sanguin correspondant au filtre sélec- tionné, T période de récurrence de l'impulsion de signal échogra- phique, C:vitesse du son dans le milieu analysé, la somme des carrés des coefficients étant normée: n k O ak 2 = K ( 3) k= O K étant une constante. La recherche des filtres s'effectue de la façon suivante: pour chaque vitesse V fixée, on cherche à maximiser le rapport signal à bruit A( 0) en faisant varier les coeffi- cients (ak) en respectant la condition ( 1) On choisit une dynamique sur les coefficients, c'est-à-dire un pas de varia- tion régulier entre une valeur maximale (positive) et une va- leur minimale (négative) de même valeur absolue que la précé- dente, et on considère toutes les combinaisons possibles des coefficients ak (k(t lO,nl) telles que les conditions ( 1) et ( 3) soient vérifiées On retient ensuite pour chaque valeur de vitesse les combinaisons de coefficients qui donnent, pour A( 0), les valeurs les plus élevées Ensuite, pour une plage de vitesses considérée, on retient, parmi les différentes combi- naisons correspondant à autant de points de vitesse de la pla- ge considérée, la combinaison qui est commune à ces différents points Etant donné le petit nombre de coefficients des fil- tres requis (de préférence 4 comme indiqué ci-dessus), cette démarche permet de sélectionner un petit nombre de filtres (de l'ordre de 3 filtres) pour la plage de vitesses utiles. Un mode de réalisation avantageux du dispositif selon l'invention est remarquable en ce qu'il comporte les 3 filtres éliminateurs d'échos fixes à quatre coefficients sui- vants: un filtre pour vitesses faibles ao = 2/110, al = 1/f a 2 = -1/110, a 3 = -2/110, un filtre pour vitesses moyennes ao = 1/2, ai = -1/2, a 2 = -1/2, a 3 = 1/2, un filtre pour vitesses élevées ao = 1/2, ai = -1/2, a 2 = 1/2, a 3 = -1/2. La description qui suit en regard des dessins an- nexés, le tout donné à titre d'exemple non limitatif, fera bien comprendre comment l'invention peut être réalisée. La figure 1 est le schéma synoptique d'un échogra- phe ultrasonore à émission d'impulsions du type profilomètre et comportant un dispositif éliminateur d'échos fixes selon l'invention. La figure 2 est le schéma synoptique du dispositif éliminateur d'échos fixes selon l'invention, associé à un cir- cuit d'intercorrélation et d'interpolation dit: dispositif de mesure des vitesses. La figure 3 est un diagramme de temps permettant d'expliciter certaines fréquences de fonctionnement du dispo- sitif de la figure 2. Les figures 4 a et 4 b représentent des structures de filtres éliminateurs d'échos fixes selon l'invention. La figure 5 est la courbe de réponse d'un filtre selon l'art antérieur et de trois filtres en parallèle selon l'invention. La figure 6 représente des profils de vitesse ob- tenus pour et au voisinage de l'artère poplitée, en a selon l'art antérieur (filtre éliminateur d'écho fixe unique à trois coefficients), en b selon l'invention. La figure 1 montre, sous forme schématique, un dispositif de mesure de la vitese d'organes en mouvement et d'écoulement sanguins Ce dispositif, appelé profilomètre, est fondé sur le principe d'une émission d'impulsions courtes (de l'ordre de 1 ps de largeur) selon une fréquence de récurrence de l'ordre de quelques k Hz (Ni = 1/T, par exemple égale à k Hz) Le profilomètre fait généralement partie d'un appareil d'exploration par échographie ultrasonore plus complexe pou- vant comporter en outre un système de visualisation bidimen- sionnelle des tissus et des flux sanguins analysés Sur la fi- gure 1 on a représenté un transducteur ultrasonore 1 qui émet (flèche 2), pendant chaque impulsion, une onde acoustique cen- trée autour d'une fréquence centrale Nc de 5 M Hz par exem- ple Un signal d'écho à différentes profondeurs des tissus ex- plorés est reçu pendant toute la durée T qui sépare deux im- pulsions consécutives L'émetteur d'impulsions est représenté en 4 La référence 5 désigne un étage séparateur connu, qui permet d'utiliser le même transducteur 1 en émission et en ré- ception En sortie du séparateur 5 on obtient un signal analo- gique Sa sur un conducteur 6 Le signal Sa est traité par une chaîne de réception et de traitement constituée par un ampli- ficateur à gain variable 7, un échantillonneur 8, un disposi- tif éliminateur d'échos fixes 9 et un dispositif de mesure des vitesses 11 La fonction de l'amplificateur à gain variable 7 est la suivante: l'amplitude des échos reçus décroît lorsque la distance des réflecteurs augmente, essentiellement à cause de l'absorption des ultrasons par les tissus traversés Il est os de ce fait nécessaire de compenser cette atténuation pour ré- duire autant que possible la dynamique du signal numérique ap- pliqué à l'éliminateur d'échos fixes 9 Il faut que la loi des variations de gain de l'amplificateur 7 (de l'ordre de 1 d B/ps à 6 M Hz) soit parfaitement reproductible, à mieux que 0,01 d B Après échantillonnage à une fréquence Fsi par exemple égale à 20 M Hz, en 8, ce qui correspond à un suréchantillonna- ge par deux environ, le signal se présente sous la forme d'échantillons Sk(t) en entrée du Suppresseur d'échos fixes 9 L'ensemble constitué par les éléments 9 et 11, qui consti- tue le coeur de l'invention, est représenté plus en détail à la figure 2. L'échantillon Sk(t) est fourni à une pluralité de filtres Fl, F 2, F 3, référencés 13, 14, 15 respectivement, chacun étant un filtre éliminateur d'échos fixes du type passe bande à lignes à retard comme décrit plus loin en référence aux figures 4 a et 4 b Les bandes passantes des filtres Fl, F 2, F 3 sont juxtaposées et l'on peut tracer la courbe de réponse de ces filtres en fonction des vitesses Doppler V à mesurer (voir figure 5), de façon telle que ces différentes courbes de réponse, homothétiques des courbes de réponse en fréquence, se recoupent dans la zone des vitesses faibles, environ comprise entre 10 et 30 cm/s Le signal de sortie des filtres 13, 14, , est un signal échantillonné à la même fréquence Fsi que le signal d'entrée Sk(t), noté Dlj, D 2 j, D 3 j respectivement. Les signaux Dlj, D 2 j, D 3 j sont fournis chacun à un détecteur d'amplitude A 1 ( 16), A 2 ( 17), A 3 ( 18) respectivement et aussi à un élément de retard 19, 21, 22 La sortie de chaque détecteur d'amplitude est reliée à une entrée correspondante d'un compa- rateur 23 et la sortie de chaque élément de retard est reliée à une entrée correspondante d'un multiplexeur 24 Le comparateur 23 sélectionne en séquence l'amplitude la plus élevée parmi les signaux D 1 j, D 2 j, D 3 j et, au moyen d'un si- gnal de commande CO sur un conducteur multiple 25, sélection- ne, dans le multiplexeur 24, le signal Dj de plus grande am- plitude qui est de ce fait fourni à la sortie 26 du multiple- xeur 24, avec la notation DR Outre le signal de commande logique CO qui est élaboré en 23 par des circuits à seuil et des portes logiques dont le choix et l'agencement sont à la portée de l'homme de l'art, la valeur d'amplitude maximale fournie par l'un des détecteurs d'amplitude Al, A 2, A 3 est sé- lectionnée par le comparateur 23 et transmise, sous la forme du signal AR au dispositif de mesure 11 à des fins expli- quées ci-dessous Les éléments de retard 19, 21, 22 sont réa- lisés par exemple sous la forme de registres à décalage, fonc- tionnant à la fréquence d'échantillonnage d'entrée Fsi; leur fonction est de retarder les signaux de sortie des filtres d'un temps au moins égal à l'estimation des amplitudes corres- pondantes, soit le temps nécessaire pour le traitement du si- gnal à travers les détecteurs d'amplitude et le comparateur 23. Pour mieux expliciter le fonctionnement du circuit électronique de la figure 2 on a représenté un générateur d'horloge 27 qui émet quatre signaux d'horloge représentés par leurs valeurs de fréquence: Ni, Fsi, No, Fso respectivement, ces signaux d'horloge étant utilisés de façon classique dans la technique de l'échographie ultrasonore No est la récurren- ce de sortie, choisie comme un sous multiple de Ni, soit No Ni N relation dans laquelle N est le nombre d'impulsions d'émission (ou tirs) pris en compte pour établir un profil de vitesse La corrélation entre paires de signaux successifs est effectuée, et les N-M fonctions de corrélation sont moyennées comme déjà indiqué plus haut (M = 4 dans l'exemple considéré ici) Fso est la fréquence d'échantillonnage de sortie, les différents échantillons obtenus en sortie représentant les points de vi- tesse successifs du profil de vitesse recherché Fso est un sous multiple de la fréquence d'échantillonnage d'entrée; par exemple: Fso = Fsi/4 (voir la figure 3) Dans ces conditions, le signal Fsi assure la synchronisation à la fréquence d'échantillonnage d'entrée dans tous les éléments représentés à la figure 2, mis à part le multiplexeur 24 o la synchroni- sation s'effectue à la fréquence d'échantillonnage de sorite Fso, et mis à part le comparateur 23 qui n'appartient pas à la chaîne directe de traitement du signal d'information reçu en écho La fourniture du signal Fsi aux éléments autres que 23 et 24 n'est pas représentée, figure 2, pour ne pas surcharger le dessin Le signal Ni en entrée des filtres 13, 14, 15 sym- bolise la présence de retards multiples de 1/Ni ( 1/Ni = T) comme expliqué ci-dessous en référence aux figures 4 a et 4 b. Ce signal Ni qui figure aussi en entrée du dispositif de mesu- re de vitesses 11 symbolise dans ce dernier la présence de re- tards 1/Ni nécessaire au moyennage sur N-M corrélations suc- cessives comme déjà décrit De façon connue la corrélation, en 11, s'effectue de préférence comme décrit dans la demande de brevet européen 0 225 667 (incorporée ici par référence) Il s'agit d'une corrélation 1 bit Pour chaque échantillon de vi- tesse à déterminer S'l à partir de deux tirs consécutifs donnés (voir figure 3 du présent fascicule), on définit une fenêtre temporelle de largeur W qui regroupe un nombre prédé- terminé constant d'échantillons d'entrée, par exemple 32, la corrélation étant effectuée bit à bit entre ces deux fenêtres qui sont décalées entre elles exactement de la durée T, c'est-à-dire qui se correspondent entre deux tirs consécu- tifs Pour l'échantillon de vitesse suivant, les deux fenêtres des deux tirs considérés sont décalées d'un nombre prédétermi- né de périodes d'échantillons d'entrée, par exemple 4, ce qui implique un important recouvrement d'une fenêtre à la suivan- te Ceci est symbolisé par la fourniture du signal d'horloge Fso à l'estimateur de vitesse 11 Les détecteurs d'amplitude 1 1 16, 17, 18, le comparateur 23 et le multiplexeur 24 sont quant à eux synchronisés par le signal d'horloge Fso de façon que le séquencement du comparateur et par suite du multiplexeur s'ef- fectue à une fréquence élevée qui est la même que la fréquence des échantillons de vitesse de sortie S'l Cette fréquence élevée permet, au moyen des opérations de corrélation, de moyennage et d'interpolation classiques qui suivent de s'af- franchir de problèmes de discontinuités sur les profils de vi- tesses qui pourraient se manifester lors de la commutation du signal de sortie d'un filtre à un autre. Chaque détecteur d'amplitude calcule l'énergie AM(to) s'effectue pendant la durée d'une fenêtre selon la re- lation: ta +w AM(to) = E D(t) to D(t) étant le signal de sortie du filtre considéré. Le comparateur effectue en permanence la comparai- son entre les signaux AM 1, AM 2, AM 3 obtenus à la cadence Fso en sortie des détecteurs 16, 17, 18, fournit le signal de com- mande CO du multiplexeur et l'amplitude résultante AR(t O) telle que: AR(to) = MaxlAM 1 (to), AM 2 (to), AM 3 (to)l Les deux signaux CO et AR sont échantillonnés à la fréquence de sortie Fso du dispositif; ces deux signaux sont calculés à chaque période de récurrence de sortie 1/No du dispositif Les trois détecteurs d'amplitude et le comparateur sont mis à l'état initial à chaque récurrence de sortie 1/NO. Le signal CO calculé à cette même cadence est fourni au multi- plexeur à chaque période de récurrence d'entrée 1/Ni de façon que le multiplexage soit constant pendant le calcul du profil de vitesse Le signal AR est utilisé de façon connue, dans le dispositif de mesure de vitesses 11, en association avec des circuits à seuils d'amplitude pour confirmer ou infirmer la validité des points de vitesse calculés, c'est-àdire pour effectuer une opération de segmentation du profil de vitesses recherché On notera que les retards apportés par les éléments de retard 19, 21, 22 doivent, chacun, être au moins égaux à la durée W de la fenêtre d'analyse cette durée étant aussi celle qui est nécessaire à l'estimation des amplitudes AM 1, AM 2, AM 3 Ainsi les signaux Dlj, D 2 j, D 3 j retardés sont rendus syn- chrones du signal de commande CO du multiplexeur 24. La figure 4 a représente sous forme théorique la structure de chacun des filtres Fl, F 2, F 3, utilisés selon l'invention La figure 4 b illustre une réalisation pratique de ce filtre Il s'agit d'un filtre éliminateur d'échos fixes à 4 coefficients ao, al, a 2, a 3 (n = 4) On rappelle que les si- gnaux réfléchis par une cible en mouvement vérifient, d'une période de récurrence d'entrée à la suivante: Sk+ 1 (t) = Sk(t-T), T étant le décalage temporel entre signaux successifs du mou- vement de la cible: 2 VT C V: vitesse de la cible. En notation complexe, on obtient: k+ 1 = Sk e-J O 0 = 4 E Nc VT/C Le signal de sortie de chaque filtre s'exprime par: n Dj = E ak Sk, k=O n Dj = So E ak eik O k=O Comme le montre la structure classique des filtres représentés aux figures 4 a et 4 b On démontre que les rapports signal sur bruit d'entrée et de sortie des filtres vérifient alors la re- lation: g N -jk O k=O soit x SNRI ( 4) n SNR Ek ak 2 soit: | N -jk O Ek O ak e SN Rn = A( 0) = O ( 2) SNRI N d E akp Améliorer le rapport signal sur bruit entre l'en- trée et la sortie de chaque filtre Fl, F 2, F 3 s'obtient en rendant maximal le rapport A( 0), tout en respectant la condi- tion d'élimination d'échos fixes n E ak = ( 1) D'autre part, pour obtenir une comparaison d'am- plitude correcte entre les différents filtres recherchés, il convient de normer l'amplitude de sortie par la constante K / n IV E ak 2 = K ( 3) k=O De manière équivalente, il suffit d'imposer la contrainte n E ak 2 = k= O Une autre contrainte que l'on se donne en plus de celles exprimées par les relations ( 1) et ( 3) consiste, pour la recherche des coefficients ak à s'imposer une dynamique sur ces coefficients, comme déjà décrit ci-dessus Alors qu'il y aurait théoriquement un filtre différent à déterminer pour chaque point de vitesse du profil de vitesse à établir, ce qui s'avérerait impraticable pour la conception du circuit, le respect des trois contraintes précitées conduit expérimentale- ment à des groupes de filtres dont les caractéristiques sont proches, pour chaque groupe, et il est possible de choisir pour chaque groupe un filtre particulier dont les coefficients s'expriment de façon simple, ce filtre particulier étant choi- si pour couvrir une gamme de vitesses particulière On est ainsi conduit à sélectionner trois filtres, dont la courbe de réponse combinée, en amplitude en fonction des vitesses est représentée en trait plein à la figure 5, pour une fréquence de récurrence d'entrée Ni de 5 k Hz Cette réponse combinée est constituée par les parties de courbe en trait plein qui se re- joignent deux à deux: RF 1, RF 2, RF 3. * Le filtre F 1, à front raide, est sélectionné par le multiple- xeur 24 pour la mesure des vitesses les plus faibles (jusqu'à 16 cm/s); ses coefficients sont ao = 2/110, al = 1//10, a 2 = -1/110, a 3 = -2/110 Le filtre F 2 est sélectionné pour la mesure des vitesses com- prises entre 16 cm/s et 25 cm/s; ses coefficients sont: ao = 1/2, ai = 1/2, a 2 = -1/2, a 3 = 1/2 Le filtre F 3, sélectionné pour la mesure des vitesses supé- rieures à 25 cm/s, a comme coefficients: ao = 1/2, ai = -1/2, a 2 = 1/2, a 3 = -1/2 Pour les trois filtres, on obtient (voir relation ( 3)): K = 1. On a aussi représenté, à la figure 5, pour un point de vitesse particulier, PV, la façon dont s'opère la sélection du filtre qui fournit en sortie le signal de plus grande amplitude, en l'occurrence AM 2 fournie par le filtre F 2 La réponse combinée des trois filtres RF 1-RF 2- RF 3 pour l'élimination d'échos fixes peut être comparée à la courbe de réponse RF représentée en trait mixte à la figure 5 La courbe RF est la réponse d'un filtre éliminateur d'échos fixes classique, unique, à trois coefficients (ao = 2, al = -1, a 2 = -1) qui n'est satisfaisant ni pour les vitesses faibles (absence d'un front raide), ni pour les vitesses élevées (amplitudes faibles). Sur les figures Ga et Gb, on a représenté des pro- fils de vitesse obtenus pour et au voisinage de l'artère po- plitée, en Ga en utilisant un filtre classique unique à trois coefficients du type indiqué au paragraphe précédent et en 6 b selon l'invention -Les abscisses représentent la profondeur PR en nombre de points, l'abscisse 400 correspondant à la profon- deur 15 mm Les parties de profils situées au voisinage de l'abscisse 280 rendent compte de l'écoulement sanguin dans l'artère poplitée et celle située au voisinage de l'abscisse 100, de l'écoulement sanguin dans uneartériole Le profil de vitesse de la figure 6 b, qui ne présente pas de discontinui- tés, contrairement à celui de la figure Ga, est plus proche de la réalité physique des déplacements dans les milieux analy- sés, par rapport à ce dernier, et montre la supériorité du filtre de type Fl sur le filtre à réponse RF, dans la gamme des vitesses les plus faibles Le dispositif selon l'invention permet aussi, de façon non représentée, de fournir des profils de vitesse convenables pour une gamme de vitesses variant de O à 30 cm/s et plus comme c'est le cas par exemple pour l'écou- lement sanguin dans l'artère carotide. Revendications: 1 Dispositif éliminateur d'échos fixes pour échogra- phe ultrasonore à émission d'impulsions récurrentes muni d'un circuit d'intercorrélation et d'interpolation, en vue de la détermination de profils de vitesse de milieux en mouvement, comportant au moins un filtre éliminateur d'échos fixes du ty- pe passe-bande à lignes à retard, caractérisé en ce qu'il com- porte une pluralité de filtres éliminateurs d'échos fixes du type précité en parallèle, dont les bandes passantes sont jux- taposées, des détecteurs d'amplitude reliés chacun à la sortie de chacun desdits filtres et à l'entrée d'un comparateur, des éléments de retard branchés entre la sortie de chaque filtre et un multiplexeur, ce dernier recevant un signal de commande issu dudit comparateur pour sélectionner en sortie dudit mul- tiplexeur celui parmi les signaux de sortie desdits filtres qui présente la plus grande amplitude. 2: Dispositif éliminateur d'échos fixes selon la re- vendication 1, caractérisé en ce que lesdits filtres élimina- teurs d'échos fixes sont des filtres à 4 coefficients, à somme des coefficients nulle. 3 Dispositif éliminateur d'échos fixes selon la re- vendication 1 ou 2, caractérisé en ce que lesdits filtres éli- minateurs d'échos fixes, dont la somme des N coefficients est nulle, soit: n Ekak = 0, k= O sont conçus et sélectionnés pour rendre maximale la fonction du rapport signal à bruit entre la sortie et l'entrée desdits filtres A( 0): n -jk O I k ak e A( 0) = O V E ak 2 k= O avec: 0 = 41 Nc VT/C ou: Nc: fréquence du signal échographique pendant l'impulsion, V: vitesse du flux sanguin correspondant au filtre sélec- tionné, T: période de récurrence de l'impulsion de signal échogra- phique, C: vitesse du son dans le milieu analysé, la somme des carrés des coefficients étant normée: n E ak 2 = 1 k=O 4 Dispositif éliminateur d'échos fixes selon la re- vendication 2 ou 3, caractérisé en ce qu'il comporte les trois filtres éliminateurs d'échos fixes à quatre coefficients sui- vants: un filtre pour vitesses faibles: ao = 2/10, al = 1/f 10, a 2 = 1/10, a 3 = -2/10, un filtre pour vitesses moyennes: ao = 1/2, al = -1/2, a 2 = -1/2, a 3 = 1/2, un filtre pour vitesses élevées: ao = 1/2, al = 1/2, a 2 = 1/2, a 3 = -1/2. Echographe ultrasonore à émission d'impulsions comportant un dispositif éliminateur d'échos fixes selon l'une des revendications 1 à 4.