Brevet FR2655260A1 — publié le 19910607

DESCRIPTION La présente invention concerne un dispositif de me- sure et de visualisation de paramètres physiologiques d'un écoulement sanguin, comprenant une unité de mesure par écho- graphie ultrasonore de la vitesse V(t,z) dudit écoulement en fonction du temps t et de la profondeur z d'exploration. 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 mesure et de la visualisation, aux fins de diagnostic, de paramètres physiologiques caracté- ristiques desdits écoulements. L'étude clinique par échographie ultrasonore des écoulements sanguins est réalisée actuellement avec des appa- reils mettant en oeuvre l'effet Doppler qui rend compte de la différence de fréquence, appelée fréquence Doppler f D, entre l'onde émise par un transducteur piézoélectrique et l'onde re- çue après interaction avec l'écoulement considéré La vitesse V du flux sanguin est liée à f D par la relation f D = 2 (V/C)f E Cos O ( 1) f E étant la fréquence de l'onde ultrasonore émise, 3 l'angle entre le faisceau ultrasonore et la direction de l'écoulement, et C la vitesse de propagation du son. Les dispositifs Doppler connus permettent, notam- ment, d'accéder à des paramètres intéressants particulièrement les praticiens comme le spectogramme des vitesses donnant la distribution de la vitesse d'écoulement V en fonction du temps t Toutefois, le résultat obtenu est entaché d'une certaine -imprécision liée de façon inhérente au principe même de la méthode et qui provient du fait que la fréquence f E émise par le transducteur piézoélectrique présente une certaine distribution qui se répercute sur la vitesse V via la relation ( 1) rappelée ci-dessus. Aussi, le problème technique à résoudre par l'objet de la présente invention est de réaliser un dispositif de me- sure et de visualisation de paramètres physiologiques d'un écoulement sanguin, conforme au préambule, qui permettrait d'obtenir un spectogramme des vitesses plus précis et qui conduirait à la détermination de paramètres physiologiques supplémentaires. La solution au problème technique posé consiste, selon la présente invention, en ce que, la mesure de la vites- se V(t,z) étant indépendante de la fréquence de l'onde ultra- sonore utilisée, ledit dispositif de mesure et de visualisa- tion comprend, en sortie d'une mémoire de stockage des valeurs de la vitesse V(t,z), d'une part, un multiplexeur qui, pour un instant donné t, inscrit dans un registre i d'une mémoire le nombre Ni(t) d'échantillons de mesure présentant la valeur Vi à AV près, AV étant le pas de discrétisation, et, d'autre part, un dispositif d'affichage de la vitesse en fonction du temps, tel que la brillance au point (Vi,t) soit proportion- nelle à Ni(t). Parmi les dispositifs pouvant donner une mesure de vitesses indépendantes de la fréquence ultrasonore, on peut citer ceux qui fonctionnent sur le principe de la corrélation temporelle déjà décrite dans la demande de brevet européen n O O 225 667, et dont l'unité de mesure de la vitesse d'écou- lement comprend un circuit d'intercorrélation délivrant à par- tir de deux échos successifs des valeurs de fonction de corré- lation, et un circuit de multiplexage-interpolation fournis- -sant à partir desdites valeurs des fonctions de corrélation une estimée de la vitesse V(t,z). La précision obtenue sur la mesure de la vitesse d'écoulement permet d'envisager, outre la visualisation du flux sanguin en mode M, par exemple, et la représentation du spectogramme des vitesses, de calculer d'autres paramètres physiologiques jusque là inaccessibles avec les systèmes Doppler conventionnels du fait de leur mauvaise résolution en profondeur Parmi les paramètres dont la présente invention autorise la détermination figurent le diamètre d'écoulement D(t), la vitesse moyenne V(t) et le débit Q(t) du flux san- guin Les moyens techniques de calcul de ces divers paramètres seront précisés en détail dans la description qui va suivre, en regard des dessins annexés, donnés à titre d'exemples non limitatifs. La figure 1 est le schéma d'un dispositif de mesure et de visualisation selon l'invention. Le schéma de la figure 1 montre un dispositif de mesure et de visualisation de paramètres physiologiques d'un écoulement sanguin 10 Ce dispositif comprend un transducteur piézoélectrique 100 qui peut être, par exemple, une barrette multiéléments Un étage 200 d'émission relié au transducteur 100 assure la formation d'un faisceau ultrasonore d'explora- tion tandis qu'une unité 300 de mesure traite les signaux échographiques renvoyés vers le transducteur 100 de façon à fournir une estimée de la vitesse V(tz) de l'écoulement en fonction du temps t et de la profondeur z d'exploration. De façon classique, l'étage d'émission 200 comprend un séquenceur composé d'un oscillateur et d'un diviseur de fréquence qui commande à la fréquence de récurrence 1/T choi- sie un générateur dont les signaux électriques d'excitation sont envoyés vers le transducteur 100 qui les convertit en trains périodiques de signaux impulsionnels ultrasonores Un séparateur 101 des étages d'émission 200 et de mesure 300 est inséré entre le transducteur 100 et lesdits étages 200, 300 et évite l'aveuglement des circuits de mesure par les signaux d'émission. L'unité 300 de mesure de la vitesse V(tz) montrée à la figure 1 comporte un éliminateur 301 d'échos fixes, à deux points par exemple, qui fournit du signal reçu Sn(t Z) un signal dn(tz) débarassé de ses composantes fixes dues aux réflexions spéculaires sur les parois du vaisseau sanguin examiné Le signal dn(t,z) issu de l'éliminateur 301 d'échos fixes est ensuite traité selon la méthode dite d'intercorréla- tion décrite dans la demande de brevet européen N O O 225 667 et qui utilise le fait que les signaux ultrasonores rétrodiffusés par une cible en mouvement sont reliés par l'équation suivante: dn+l(t) = dn(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 le temps supplémentaire nécessaire à l'onde ultrasonore pour parcourir le trajet transducteurcible-transducteur, d'un tir à l'autre Autrement dit: T = 2 VT/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 cherchée. La fonction d'intercorrélation entre dn(t) et dn+l(t) définie par: to+W to+W C (to,u) = I d (t+u)d (t)dt n,n+ 1 to n+ 1 n vérifie: Cn,n+l(to,u)=Cn,n(to,u-T) Le temps to est lié à la profondeur d'exploration z par to= 2 z/C, et W est la largeur de la fenêtre d'intégration. La fonction Cnn(t O,u) est une fonction d'autocor- rélation, et, de ce fait, est maximale pour u=o Ainsi, une mesure du décalage temporel r, 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 A cet effet, la fonction d'inter- corrélation est échantillonné, à un pas d'échantillonnage At, entre Umin=-I At et Umax=I At par pas de 1 de façon à obte- nir 2 I+ 1 valeurs de fonctions de corrélation La valeur maxi- mum de ces 2 I+ 1 valeurs correspondant à u=uo permet de mesurer r en utilisant l'égalité T=uo Le calcul des fonctions de cor- rélation est réalisé par le circuit 302 d'intercorrélation montré à la figure 1. Afin de s'affranchir des erreurs inhérentes à l'échantillonnage dans la détermination du maximum de la fonc- tion de corrélation, il est possible d'utiliser un circuit 303 de multiplexage-interpolation fournissant à partir des valeurs des fonctions de corrélation une estimée plus précise de la vitesse et la valeur du pic de corrélation correspondant La demande de brevet européen n'O 225 667 donne un exemple de ce type de traitement du signal échographique dans lequel la cor- rélation entre signaux est une corrélation dite " 1 bitn en ce sens que les signaux dn+l et dn utilisés précédemment sont réduits au signe du signal ultrasonore On sait que dans ce cas, le pic de la fonction de corrélation est de forme trian- gulaire isocèle La connaissance de cette forme permet à par- tir du point le plus haut et de ses deux voisins de reconsti- tuer complètement, par interpolation linéaire, le pic de cor- rélation et donc de déterminer avec précision l'emplacement de uo. Contrairement aux vélocimètres Doppler tradition- nels, la vitesse d'écoulement V(t,z) ainsi déterminée présente l'avantage d'être insensibles à la dispersion de la fréquence de l'onde ultrasonore utilisée, ce qui permet d'envisager une exploitation beaucoup plus complète des résultats. Les valeurs trouvées pour la vitesse V(t,z) sont stockées dans une mémoire 304 de stockage en vue de traite- ments ultérieurs En particulier, à l'aide de moyens de visua- lisation connus, il est possible d'afficher sur un premier écran 305 une image de l'écoulement étudié selon la présenta- tion, dite en "mode M", de l'évolution au cours du temps du profil de vitesse dans le vaisseau Le codage de la vitesse est celle utilisée dans les sytèmes CFM (Color Flow Mapping) par exemple, le rouge code un sens de déplacement, le bleu le sens inverse et l'intensité de la vitesse est codée par l'intensité de la couleur. -30: Un autre traitement particulièrement important pour le médecin est la réalisation d'un spectogramme de vitesse précis A cet effet, le dispositif de mesure et de visualisa- tion de la figure 1 comprend, en sortie de la mémoire 304 de stockage, d'une part, un multiplexeur 401 qui, pour un instant -t donné, inscrit dans un registre i d'une mémoire 402 le nom- bre Ni(t) d'échantillons de mesure présentant la valeur Vi à AV près, AV étant le pas de discrétisation, et, d'autre part, un dispositif 403 d'affichage de la vitesse en fonction du temps tel que la brillance au point (Vi,t) est proportionnelle à Ni(t). La figure 1 montre également comment peuvent être calculés d'autres paramètres physiologiques caractéristiques de l'écoulement Il s'agit d'abord du diamètre D(t) de l'écoulement dont la mesure met en oeuvre, en sortie de l'éliminateur 301 d'échos fixes, un circuit 404 de calcul de l'énergie locale E(t,z) = E d 2 (t,z) et un circuit 405 de N n calcul du diamètre d'écoulement D(t) = z 2 (t)-zt(t) constitué par un détecteur de seuil de valeur Eo ajustable, zl(t) et z 2 (t) étant définis par E(t,zl(t))=E(t,z 2 (t))= Eo. Connaissant, à chaque instant t, le diamètre D(t) de l'écoulement, il est possible de calculer la vitesse moyenne V(t) du flux sanguin grâce à un circuit 406 constitué par un additionneur donnant EDV(t,z) et par un diviseur par M, M étant le nombre d'échantillons de la vitesse V(t,z) utilisés sur le segment lZl,z 2 l. De même, le débit Q(t) peut être mesuré par un circuit 407 comprenant, d'une part, un circuit 408 de calcul du centre de gravité zo(t)=EDV(t,z) z/V(t), constitué d'un additionneur donnant EDV(t,z)z et d'un diviseur par la vitesse moyenne V(t) calculée par le circuit 406, et, d'autre part, un additionneur 401 donnant EDV(t,z)1 z-zol suivi d'un multiplieur par la constante A=Ucos 8/sin 2 e, 8 étant l'angle de l'écoulement avec le faisceau ultrasonore. REVENDICATIONS 1 Dispositif de mesure et de visualisation de paramètres physiologiques d'un écoulement sanguin ( 10), comprenant une unité de mesure ( 300) par échographie ultrasonore de la vitesse V(t,z) dudit écoulement en fonction du temps t et de la profondeur z d'exploration, caractérisé en ce que, la mesure de la vitesse V(t,z) étant indépendante de la fréquence de l'onde ultrasonore utilisée, ledit dispositif de mesure et de visualisation comprend, en sortie d'une mémoire ( 304) de stockage des valeurs de la vitesse V(t,z), d'une part, un multiplexeur ( 401) qui, pour un instant t donné, inscrit dans un registre i d'une mémoire ( 402) le nombre Ni(t) d'échantillons de mesure présentant la valeur Vi à AV près, AV étant le pas de discrétisation, et, d'autre part, un dispositif ( 403) d'affichage de la vitesse en fonction du temps tel que la brillance au point (Vi,t) est proportionnelle à Ni(t). 2 Dispositif de mesure et de visualisation selon la revendication 1, caractérisé en ce que, comprenant un circuit ( 404) de calcul de l'énergie locale E(t,z), il comporte également un circuit ( 405) de calcul du diamètre d'écoulement D(t)=z 2 (t)-zi(t) constitué par un détecteur de seuil de valeur Eo, z 2 (t) et z 1 (t) étant définis par E(t,zl(t))=E(t,z 2 (t))=Eo. 3 Dispositif de mesure et de visualisation selon la revendication 2, caractérisé en ce qu'il comporte un circuit ( 406) de calcul de la vitesse moyenne V(t) constitué par un additionneur donnant EDV(t,z) et par un diviseur par M, M étant le nombre d'échantillons de mesure sur le segment lzl,z 2 l. 4 Dispositif de mesure et de visualisation selon la revendication 3, caractérisé en ce qu'il comporte un circuit ( 407) de calcul du débit Q(t) comprenant, d'une part, un circuit ( 408) de calcul du centre de gravité zo(t) de l'écoulement, constitué d'un additionneur donnant EDV(t,z)z et d'un diviseur par la vitesse moyenne V(t), et, d'autre part, un additionneur ( 409) donnant EDV(t,z) z-zol suvi d'un multiplieur par la constante A=Ucos 2/sin 2 e, $ étant l'angle de l'écoulement avec le faisceau ultrasonore. Dispositif de mesure et de visualisation selon l'une quelconque des revendications 1 à 4, caractérisé en ce que l'unité de mesure de la vitesse d'écoulement comprend un circuit ( 302) d'intercorrélation délivrant, à partir de deux échos successifs, des valeurs de fonctions de corrélation, et un circuit ( 303) de multiplexageinterpolation fournissant à partir desdites valeurs de corrélation une estimée de la vitesse V(t,z).