Brevet FR2617982A1 — publié le 19890113

- 1 _ DISPOSITIF D'ELIMINATION D'ECHOS FIXES POUR ECHOGRAPHE ULTRA- SONORE'. La présente invention concerne un dispositif d'é- limination d'échos fixes pour échographe ultrasonore, compor- tant un filtre du type à lignes à retard. Le dispositif d'élimination d'échos fixes selon l'invention trouve une application particulièrement avantageu- se dans la mesure et l'imagerie des vitesses d'écoulements sanguins par échographie ultrasonore. Le problème technique général à résoudre par tout échographe ultrasonore destiné à la mesure et à 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 dB la réflecti- vité du sang (globules rouges). L'utilisation d'un dispositif d'élimination d'échos fixes est donc indispensable, avant toute estimation du profil de vitesse, pour réduire la dynami- que des signaux à. traiter et pour mesurer correctement les vi- tesses 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 à ce problème technique général est donnée, par exemple, dans la demande de brevet français né 2 580 076 qui décrit un dispositif d'élimination d'échos fixes pour échographe ultrasonore, comportant un filtre, du type à lignes à retard, composé en l'espèce.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, 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 2- é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 5kHz et une fréquence d'é- mission de 5MHz, un signal correspondant à v = 5cm/s est atté- nué de 30dB. 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 augmenter la pente du filtre mentionné plus haut, dans la zone de coupu- re en augmentant son ordre. Appliqué à la mesure des vitesses par échographie, ceci revient à traiter non plus 2 lignes é- chographiques successives mais d'avantage, soit M, à l'aide d'un filtre du type à lignes à retard, composé de M lignes pa- rallèles i de retard respectif (i-1)T. De façon à assurer l'é- limination des échos fixes, les M lignes sont pondérées par des coefficients dont la somme est nulle. Puis toutes les li- gnes sont additionnées pour ensuite être traitées par un cir- cuit estimateur de vitesse, par corrélation, décalage de pha- se, etc... Cependant, il est difficile avec ce type de filtre 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 vais- seaux périphériques et N=8 pour le coeur), et si M lignes sont utilisées pour le filtre d'élimination des échos fixes, il ne reste en réalité que N-M+1 lignes utilisables pour l'estima- tion de vitesse, ou alors, pour avoir le même nombre de lignes à traiter, il faudrait diminuer la cadence d'image. Le problème à résoudre pour l'objet de la présente demande est de proposer un dispositif d'élimination d'échos fixes pour échographe ultrasonore, comportant un filtre - 3 composé, d'une part, de M lignes parallèles i(i=1,...,M) de retard respectif (i-U)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 grâce auquel on puisse ob- tenir une estimation satisfaisante des faibles vitesses d'écoulements sanguins, notamment à proximité des parois des vaisseaux, sans diminution importante des lignes à traiter, donc sans réduction de la cadence d'image. La solution au problème technique posé consiste, selon la présente invention, en ce qu'audit filtre, appelé premier filtre, est associé en parallèle au moins un deuxième filtre de même structure, conjugué du premier filtre et dépha- sé de (M-1)T, la réponse impulsionnelle h2(t) du deuxième filtre étant reliée à la réponse impulsionnelle h1(t) du premier filtre par: h2(t) = hl((M-1)T-t) (1) On obtient ainsi deux filtres d'ordre M-1 ayant, en module, la même réponse en fréquence, condition essentielle pour que ces deux filtres puissent être associés. En effet, de l'équation (1) on tire la relation entre les réponses en fré- quence H1(f) et H2(f): H2(f) = Hl(f) exp(j2Tf(M-1)T) d'o tH2(f)I = IH1(f) l D'autre part, on dispose 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 (corrélation 1 bit, décalage de phase,...), le gain est néanmoins appréciable. La description qui va suivre, en regard des dessins annexés donnés à titre d'exemples non limitatifs, fera bien -4- comprendre en quoi consiste l'invention et comment elle peut être réalisée. La figure 1 est le schéma d'un dispositif d'élimination d'échos fixes selon l'invention. La figure 2 est le diagramme de la réponse en vitesse d'un dispositif d'élimination d'échos fixes à trois points. La figure 3 est un schéma montrant comment le dispositif d'élimination d'échos fixes peut être associé à un estimateur de vitesse. La figure 1 donne le schéma d'un dispositif 10 d'élimination d'échos fixes pour échographe ultrasonore. Ce dispositif se trouve placé, d'une façon générale, entre le circuit de réception du faisceau ultrasonore et l'estimateur de vitesse. Comme le montre la figure 1, le dispositif 10 d'élimination d'échos fixes représenté comporte un filtre 100 composé, d'une part, de M lignes i(i=1,...,M) de retard respectif (i-1)T, T étant la période de récurrence de l'échographe. Ces retards sont obtenus à l'aide des lignes à retard 103. D'autre part, le filtre 100 comprend également des moyens de pondération 1011,...101M des M lignes i. De façon à assurer l'élimination des échos fixes, la somme des coefficients al,... aM affectés auxdites lignes est nulle: M E ai = O i=1 Enfin, un additionneur 102 fait la somme des lignes i ainsi retardées et pondérées. La réponse impulsionnelle hl(t) du filtre 100 est égale à: M hl(t) = E ai6 (t-(i-1)T) i=1 o 6(t) est la fonction de Dirac. - 5 - La réponse du filtre 100 à un signal échographique sI(t) provenant du circuit de réception de l'échographe, I étant le rang du tir correspondant, est donc donnée par: M XI(t) = E aisI-(i-1)(t) i=1 La réponse en vitesse G(v) d'un filtre 100 pour le- quel M vaut 3 et les coefficients a1, a2 et a3 respecti- vement 1, -1/2, -1/2, est montrée à la figure 2. On voit que ce filtre a une sensibilité aux faibles vitesses beaucoup plus grande que le filtre d'ordre 1 correspondant pour lequel a1 = 1 et a2 = -1. En particulier, pour une vitesse de 6cm/s, le gain en réponse est de 6dB. Par contre, on voit que sur N lignes successives à traiter (I=1,...,N), M sont utilisées pour l'élimination des échos fixes et que par conséquent seules N-M+1 serviront effectivement à estimer la vitesse. D'o une précision de mesure moins grande ou l'obligation de diminuer la cadence d'image. Comme on peut le voir à la figure 1, le dispositif d'élimination d'échos fixes comporte également, en parallè- le avec le filtre 100, appelé premier filtre, un deuxième fil- tre 200 qui présente la même structure que ledit premier fil- tre 100. Le deuxième filtre 200 est conjugué du premier filtre et déphasé de (M-1)T, sa réponse impulsionnelle h2(t) étant reliée à hl(t) par: h2(t) = hl((M-1)T-t) (2) Comme, par ailleurs, h2(t) est aussi donné par M h2(t) = E a'i6(t-(i-1)T) i=1 - 6 -- o les coefficients a'i sont ceux fournis par les moyens de pondération 201i, on voit qu'il existe nécessairement une relation entre les coefficients ai et a'i. Par exemple, si M est égal à 3, on a: hl(t) = a16(t) + a26(t-T) + a36(t-2T) h2(t) = a'16(t) + a'26(t-T) + a'36(t-2T) La relation (2) implique: a'1=a3 a'2=a2 a'3=a1 En sortie de l'additionneur 202, le signal fourni par le fil- tre 200 est: M YI(t) = E a'i sI-(i-1)(t) i=1 La réponse en vitesse G(v) du deuxième filtre 200 est la même que celle du premier filtre 100. Avec des coeffi- cients a'l, a'2, a'3 égaux respectivement à -1/2, -1/2, 1, la réponse G(v) du deuxième filtre est donc donnée par la courbe en trait mixte de la figure 2. Lorsque la relation (2) est satisfaite, et comme on l'a montré plus haut, les deux filtres 100 et 200 ont, en mo- dule, la même réponse en fréquence, et donc en énergie. Le dispositif 10 d'élimination d'échos fixes selon l'invention peut donc être considéré comme un filtre unique ayant pour réponse en vitesse et pour module de réponse en énergie respectivement la réponse en vitesse commune G(v) et le module de réponse en énergie commun aux deux filtres 100 et 200. L'avantage du dispositif 10 d'élimination d'échos-fixes montré à la figure 1 est qu'il permet de traiter 2(N-M+1) li- gnes pour estimer la vitesse au lieu de N-M+1 comme c'est le cas pour chacun des filtres 100 et 200 pris séparément. Il est bien clair cependant que les signaux XI(t) et YI(t) provenant des filtres 100 et 200 ne sont pas indé- 7 - pendants dans la mesure o ils sont construits à partir des mêmes signaux d'entrée si(t). Toutefois, ils sont différents du fait queles coefficients de pondération ai et a'i sont différents. Pour tirer parti de cette différence, l'estimateur de vitesses qui traite les signaux XI(t) et YI(t) doit u- tiliser une opération non linéaire. La figure 3 montre un exemple de la façon dont un dispositif 10 d'élimination d'échos fixes conforme à l'inven- tion peut être associé à un estimateur 20 de vitesse. Les si- gnaux XI(t) et YI(t) issus respectivement des filtres 100 et 200 sont retardés d'une période de récurrence T par les li- gnes à retard 203 de sorte que les signaux successifs XI(t), XT_1(t), d'une part, et YI(t), YIl(t), d'autre part, puissent être traités par les circuits, identiques, 204. Ces circuits 204 peuvent être des corrélateurs qui permettent d'extraire les fonctions de corrélation CI et C'I. Un exemple de tels corrélateurs est donné dans la demande de bre- vet français n" 85 17851 qui montre également comment la mesu- re des fonctions de corrélation permet de déterminer les vi- tesses d'écoulement recherchées. Ensuite, ces fonctions de corrélation CI et C'I sont non seulement accumulées sur I (I=1,...,N-M+1) mais également entre elles, dans un accumula- teur 205, ce qui augmente de façon non négligeable le rapport signal à bruit de l'estimateur de vitesse. A la place de corrélateurs, on peut également uti- liser pour les circuits 204 des déphaseurs, le déphasage entre deux lignes échographiques étant fonction de la vitesse d'écoulement. - 8 - REVENDICATIONS: 1. Dispositif d'élimination d'échos fixes pour écho- graphe ultrasonore, comportant un filtre (100) composé, 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 (101i) des M lignes i, la somme des coefficients (ai) affectés auxdites lignes étant nulle, et, enfin, d'un additionneur (102) des lignes i ainsi retardées et pondérées, caractérisé en ce qu'audit fil- tre (100), appelé premier filtre, est associé en parallèle au moins un deuxième filtre (200) de même structure, conjugué du premier filtre (100) et déphasé de (M-1)T, la réponse impul- sionnelle (h2(t)) du deuxième filtre étant reliée à la ré- ponse impulsionnelle (h1(t)) du premier filtre par: h2(t) = hl((M-1)T-t) 2. Dispositif selon la revendication 1, caractérisé en ce que, M étant égal à 3 et lesdits coefficients a1, a2 et a3 étant respectivement 1, -%, -, les réponses impulsion- nelles des premier et deuxième filtres sont respectivement: h1(t) = 6(t) 6(t-T) - 6(t-2T) h2(t) = -%6(t) - %6(t-T) + 6(t-2T) o 6(t) est la fonction de Dirac.