Brevet FR2646918A1 — publié le 19901116

DESCRIPTION *DISPOSITIF DE MESURE DE LA VITESSE D'ECOULEMENTS SANGUINS PAR ECHOGRAPHIE ULTRASONORE A CADENCE DE MESURE AMELIOREE". La présente invention concerne un dispositif de me- sure de la vitesse d'écoulements sanguins à partir d'une sé- quence de N signaux échographiques successifs, comprenant un éliminateur d'échos fixes à M points, suivi d'une unité de me- sure de ladite vitesse par corrélation-sommation-interpolation des N-M+1 signaux indépendants issus dudit éliminateur d'échos fixes. L'invention trouve une application particulièrement avantageuse dans le domaine général de l'exploration échogra- phique d'écoulements sanguins dans les vaisseaux, et plus par- ticulièrement dans celui de la visualisation de tels écoule- ments. On connaît de la demande de brevet européen n" O 225 667 (PHF 85/593) un dispositif de mesure conforme au préambule dans lequel ladite unité de mesure est constituée d'une mémoire dans laquelle sont stockés lesdits M-N+1 signaux indépendants, d'un corrélateur fournissant M-N fonctions d'in- tercorrélation des M-N+1 signaux mémorisés, d'un additionneur donnant la fonction d'intercorrélation moyenne, et d'un cir- cuit d'interpolation, qui permet de déterminer le maximum de la fonction d'intercorrélation moyenne, appelée aussi pic de corrélation, dont la position est directement liée à la vites- se d'écoulements recherchée. On peut, en effet, rappeler que le dispositif connu de l'état de la technique utilise le fait que les signaux ultrasonores successifs rétrodiffusés par une cible en mouvement, lorsque l'émission est récurrente à la pé- riode de récurrence T, sont reliés par l'équation suivante: Sn+1(t) = Sn(t-T) Ceci signifie que le signal n+1 est la réplique du signal n précédent à un décalage temporel T près. Ce dernier représente 2646 9 18 le temps supplémentaire nécessaire à l'onde ultrasonore pour parcourir le trajet transducteur-cible-transducteur, d'un tir à l'autre. Autrement dit: T = 2VT/C o V est la vitesse de la cible et C la vitesse du son. On voit qu'une mesure de T permet de mesurer la vitesse V cher- chée. La fonction d'intercorrélation entre Sn(t) et Sn+1(t) définie, sur une fenêtre de largeur W, par la rela- tion: Cn (n+lto,u) = ftO+W Sn (t+u) S (t)dt nn+l to n+1 n vérifie: Cn,n+ 1l(tO,u)=Cnn(to,u-T) Le temps to est lié à la profondeur d'exploration z par to=2z/C. La fonction Cnn(to,u) est une fonction d'autocor- rélation, et, de ce fait, est maximale pour u=o. Ainsi, une mesure du décalage temporel T, et donc de la vitesse V, peut se faire en cherchant pour quel paramè- tre u la fonction Cn,n+1(to,u) est maximale. Si ce maximum est obtenu pour la valeur uo de u, on en déduit T en utilisant l'égalité uo=T et V=uoC/2T. Le dispositif décrit dans la demande de brevet européen mentionnée plus haut a apporté à la mesure de la vi- tesse des écoulements sanguins, et à leur visualisation, une contribution décisive. Il faut toutefois remarquer que la pré- cision du résultat obtenu dépend du nombre de fonctions d'in- tercorrélation dont on fait la moyenne, c'est-à-dire N-M pour le dispositif connu. Dans le but de réduire la variance de la mesure, il y a donc avantage à augmenter le nombre N de si- gnaux à traiter simultanément, ou à diminuer le nombre de points M de l'éliminateur d'échos fixes. 2646 9 18 En pratique, il est guère possible de réduire sensiblement M au-dessous de 3 ou 4. On sait en effet que l'utilisation d'un dispositif d'élimination d'échos fixes est indispensable, avant toute estimation de profil de vitesses, du fait de la grande réflectivité des parois des vaisseaux qui peut dépasser de 40dB la réflectivité du sang (globules blancs). L'élimina- teur d'échos fixes le plus simple est celui dit à 2 points, il est composé d'une ligne retardé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 addition- neur. 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. Cependant, cette technique présente l'incon- vénient majeur d'atténuer également les signaux correspondant aux faibles vitesses d'écoulement. A titre 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 si- gnal correspondant à v=Scm/s est atténué de 30dB. Ceci rend difficile, voire impossible, la mesure des vitesses d'écoule- ment là o elles sont les plus faibles, c'est-à-dire près des parois des vaisseaux. Or, la connaissance de ces vitesses est très importante pour l'étude et le diagnostic clinique des ar- tères, par exemple. C'est pourquoi il est préférable d'utili- ser un éliminateur d'échos à au moins 3 points comme celui dé- crit dans la demande de brevet européen né 0 298 569 (PHF 87/554) qui permet de s'affranchir complètement des échos produits par les tissus fixes sans pour autant réduire de façon excessive les signaux provenant des écoulements de fai- ble vitesse. On en conclut que la seule façon de pouvoir amélio- rer la précision de la mesure serait d'augmenter le nombre N de signaux successifs utilisés pour obtenir une estimée de la vitesse. Toutefois, il est clair que cette solution a pour 2646 9 18 conséquence immédiate de diminuer la cadence de mesure, ou la cadence d'image dans le cas d'une visualisation, celle-ci étant donnée par 1/NT. Aussi, le problème technique à résoudre par l'objet de la présente invention est de réaliser un dispositif de me- sure conforme au préambule qui permettrait d'obtenir une meil- leure précision de mesure sans en altérer la cadence, ou, ce qui revient au même, d'augmenter la cadence de mesure sans di- minuer la précision. O10 La solution au problème technique posé consiste, selon la présente invention, en ce que ladite unité de mesure comporte, d'une part, NF voies parallèles de traitement constituées, chacune, d'un filtre, d'une mémoire de stockage des N-M+1 signaux filtrés, et d'un corrélateur fournissant N-M fonctions d'intercorrélation, et, d'autre part, un addition- neur qui réalise la moyenne des NF(N-M) fonctions d'inter- corrélation obtenues, un circuit d'interpolation appliqué à la fonction d'intercorrélation moyenne fournissant une estimée de la vitesse recherchée. L'effet technique essentiel qui se trouve à la base de l'invention réside donc dans la multipli- cation par un facteur NF du nombre de fonctions d'intercor- rélation indépendantes à partir desquelles est calculée la fonction d'intercorrélation moyenne. Ceci est rendu possible par le fait que les opérations de filtrage conservent les re- lations temporelles entre les signaux et conduisent, par conséquent, aux mêmes résultats en corrélation. L'indépendance des fonctions d'intercorrélation obtenues par des voies paral- lèles de traitement différentes est assurée par l'absence de corrélation entre les NF filtres. Ainsi, pour la même fré- quence de mesure 1JNT, on obtient une estimée de vitesse dont la précision est améliorée d'un facteur fNF. Inversement, si 1/NoT est la cadence de mesure du dispositif connu décrit plus haut donnant la même précision que le dispositif selon 2646 9 18 l'invention, la relation entre N et No est donnée par: NF(N-M)=NoM soit: N=(No-M)/NF+M avec No=15, M=3 et NF=6, on a N=5, ce qui correspond à un gain d'un facteur 3 en cadence de mesure. On remarque, cepen- dant, que lorsqu'on augmente le nombre NF de voies parallè- les de traitement, on obtient une valeur limite pour N égale à M. O10 Le choix des NF filtres est guidé par l'idée de réaliser artificiellement NF transducteurs piézoélectriques différents en appliquant des filtrages différents à un seul transducteur, sans perdre de vue le souci de ne pas dégrader la résolution du signal. Aussi, il est prévu que, dans un mode de réalisation particulier de l'invention, Fe étant la fré- quence d'échantillonnage du dispositif, chaque filtre est dé- fini sur l'intervalle]O,Fe/2] par une réponse en module égale à 1 et une phase aléatoire, et par une réponse nulle pour la fréquence zéro. Chacun des NF filtres est donc caractérisé par sa répartition aléatoire de phase. On peut également re- marquer que ce type de filtres préserve la bande de fréquence des signaux. Enfin, de façon connue, on peut utiliser un mode particulier de mise en oeuvre utilisant la corrélation dite "lbitn, et dans lequel chaque mémoire de stockage ne mémorise que le signe des N-M+l1 signaux issus du filtre correspondant. Dans ce cas, le pic de la fonction d'intercorrélation est de forme triangulaire isocèle. La cQnnaissance de cette forme permet à partir du point le plus haut et de ses deux voisins de reconstituer complètement, par interpolation linéaire, le pic de corrélation et donc de déterminer avec précision l'em- placement uo de son maximum. La description qui va suivre en regard des dessins annexés, donnés à titre d'exemples non limitatifs, fera bien comprendre en quoi consiste l'invention et comment elle peut être réalisée. 2646 9 18 La figure I est le schéma d'un dispositif de mesure selon l'invention. Les figures 2a et 2b donnent la réponse en fréquen- ce et la phase de filtres utilisés dans le dispositif de la figure 1. La figure 3 donne quelques exemples de réponses temporelles de filtres qui ont une réponse en fréquence conforme à celle de la figure 2. La figure 1 montre, sous forme schématique, un dis- positif de mesure de la vitesse d'écoulements sanguins à par- tir d'une séquence de N signaux échographiques successifs Sn(t), n variant par pas de 1 de 1 à N. Ces signaux échogra- phiques proviennent d'une unité 100 d'émission-réception, non représentée, mais comprenant, de façon classique, un transduc- teur piézoélectrique qui convertit en trains périodiques d'im- pulsions ultrasonores de période de récurrence T les signaux électriques d'excitation qu'il reçoit d'un étage d'émission. Un étage de réception fournit à sa sortie les signaux échogra- phiques Sn(t), échantillonnés à la fréquence d'échantillon- nage Fe, renvoyés vers le transducteur piézoélectrique par le milieu soumis à observation. Le dispositif de mesure de vites- se de la figure 1 comporte un éliminateur 200 d'échos fixes à M points qui, de la séquence des N signaux échographiques Sn(t), donne une séquence de N-M+1 signaux di(t) indépen- dants et exempts des composantes de forte amplitude provenant des parois fixes des vaisseaux. Comme indiqué plus haut, un exemple d'éliminateur d'échos fixes à M=3 points qui pourrait convenir au dispositif selon l'invention est donné dans la de- mande de brevet européen n 0 298 569. Ainsi que le montre la figure 1, l'éliminateur 200 d'échos fixes est suivi d'une unité de mesure 300 de la vites- se V de l'écoulement sanguin considéré, par corrélation-somma- tion-interpolation des N-M+1 signaux indépendants di(t) is- sus de l'éliminateur d'échos fixes, conformément à l'enseigne- ment général de la demande de brevet européen né 0 225 667. Comme on peut le voir à la figure 1, l'unité de mesure 300 comporte, en premier lieu, NF voies parallèles de 2646 9 18 traitement notées 310j, avec j=1,2,...,NE. Chacune de ces voies de traitement comprend un filtre 311j dont la réponse H(f) en fonction de la fréquence f est schématisée sur les figures 2a et 2b: le module IH(f)l vaut 0 pour f=O et 1 de la fréquence 0 à la fréquence Fe/2, Fe étant la fréquence d'échantillonnage, tandis que la phase 0(f) est, dans le même intervalle, une fonction aléatoire de la fréquence: 0(f)=Ran(f) Quelques exemples de réponses temporelles h(t) de filtres de ce type sont donnés sur la figure 3. Les quatre filtres montrés sur cette figure ont la propriété d'être complètement décorrélés entre eux, ce qui assure l'entière indépendance des NF voies de traitement les unes par rapport aux autres, cette propriété étant essentielle pour assurer l'efficacité du dispositif selon l'invention. L'opération effectuée par les filtres 311j est une convolution qui, du signal entrant di(t) donne, en sortie, le signal d?(t) égal à: d?(t) = di(t) x h.(t) i 2. J Ensuite, la voie 310j de traitement comporte une mémoire 312j de stockage des valeurs des N-M+1 signaux filtrés d?(t) servant au calcul, par un corrélateur 313j, des N-M fonctions d'intercorrélation définies par: C? (tou-T) = Ci (tou) = di+W(t+u) d?(t)dt. 1. 1,i+1 itoi+l Les signaux di(t) utilisés pour déterminer les fonctions d'intercorrélation peuvent être limités aux seuls signes des signaux euxmêmes. Cette méthode de mise en oeuvre simplifie considérablement les calculs de corrélation et la taille des mémoires 312j dans la mesure o les signaux traités n'occupe qu'un mot d'1 bit. Le filtrage par les filtres 311j n'affecte pas le résultat des corrélations, car la convolution ne modifie pas les relations temporelles entre signaux, en d'autres termes: t +Wdi (t+u) di(t)dt=t+Wd1(t+u) di(t) dt quelque soit j to ' i+1 i' '=to i+l Cette propriété permet d'obtenir effectivement NF(N-M) fonctions d'intercorrélation indépendantes dont la moyenne est effectuée par un additionneur 320. Comme conséquence de l'impédance des fonctions de corrélation, la fonction de corrélation moyenne: J., (t,.T -C(to,u-T)=i.j Cii (tou-T) présente une variance a2 divisée par NF. La fonction C(to,u-T) est finalement traitée de façon classique par un interpolateur 330 qui fournit l'estimée V de la vitesse cherchée en calculant la position uo=T=2VT/C du pic de la corrélation, dont la hauteur CMax est également mesurée en vue d'être utilisée dans un circuit de validation 340 de la valeur estimée V de la vitesse. La validation de la mesure de vitesse est indispensable. Car, en dehors des zones d'écoulement, le signal sortant de l'éliminateur 200 d'échos fixes est essentiellement du bruit. Le résultat alors fourni par le dispositif selon l'invention n'est pas l'indication d'une vitesse nulle, et il devient nécessaire de valider, ou non, ce résultat. Pour cela, on procède à deux comparaisons. D'une part, on calcule l'énergie locale E du signal en sortie de l'éliminateur d'échos fixes: E= N-M+1 E = E d2 (t) i=1 i E est ensuite comparée à un seuil Eo. Si le résultat de cette comparaison est positif, on procède alors à la comparaison du maximum CMax du pic de corrélation à un deuxième seuil El. L'estimée V de la vitesse ne sera validée que si CMax est supérieur à El. La vitesse ainsi validée est traitée pour visualisation, de façon classique, par une unité 400, non représentée, comportant un dispositif de mémorisation, conversion de balayage et codage de couleur. 2646 9 18 REVENDICATIONS 1. Dispositif de mesure de la vitesse d'écoulements sanguins à partir d'une séquence de N signaux échographiques successifs, comprenant un éliminateur (200) d'échos fixes à M points, suivi d'une unité de mesure (300) de ladite vitesse par corrélation-sommation-interpolation des N-M+1 signaux in- dépendants issus dudit éliminateur (200) d'échos fixes, carac- térisé en ce que ladite unité de mesure (300) comporte, d'une part, NF voies parallèles (310j) de traitement constituées, chacune, d'un filtre (311j), d'une mémoire (312j) de stockage des N-M+1 signaux filtrés, et d'un corrélateur (313j) fournis- sant N-M fonctions d'intercorrélation, et, d'autre part, un additionneur (320) qui réalise la moyenne des NF(N-M) fonc- tions d'intercorrélation obtenues, un circuit (330) d'interpo- lation appliqué à la fonction d'intercorrélation moyenne four- nissant une estimée de la vitesse recherchée. 2. Dispositif de mesure selon la revendication 1, ca- ractérisé en ce que, Fe étant la fréquence d'échantillonnage dudit dispositif, chaque filtre (311j) est défini sur l'inter- valle]O,Fe/2] par une réponse en module égale à 1 et une pha- se aléatoire, et par une réponse nulle pour la fréquence zéro. 3. Dispositif de mesure selon l'une des revendications 1 ou 2, caractérisé en ce que chaque mémoire (312j) de stocka- ge ne mémorise que le signe des N-M+1 signaux issus du filtre (311j) correspondant. 4. Dispositif de mesure selon l'une quelconque des re- vendications 1 à 3, caractérisé en ce qu'il dispose d'un cir- cuit (340) de validation de l'estimée obtenue pour ladite vi- tesse.