Brevet FR3008802A1 — publié le 20150123

Procédé et dispositif de cartographie de milieux fibreux. DOMAINE DE L'INVENTION La présente invention est relative aux procédés et 5 dispositifs de cartographie de milieux fibreux. ARRIERE-PLAN DE L'INVENTION Un tel procédé a déjà été décrit par exemple par Derode et Fink (Spatial coherence of ultrasonic speckle in composites, Derode A., Fink M., IEEE Trans. Ultrason. 10 Ferroelectr. Freq. Control, 1993;40(6):666-75), qui enseigne de faire émettre successivement des ondes ultrasonores focalisées par une barrette de transducteur disposée à la surface d'un matériau composite, avec plusieurs orientations de la barrette de transducteurs. On 15 calcule pour chaque tir d'onde ultrasonore, une fonction de cohérence spatiale entre les signaux captés par les transducteurs suite à la réverbération de l'onde ultrasonore émise, et on détermine la direction des fibres du matériau composite comme étant la direction de la 20 barrette de transducteurs correspondant au maximum de la fonction de cohérence spatiale. Ce procédé connu convient pour un milieu simple tel qu'un matériau composite dont les fibres sont disposées régulièrement ; il ne convient pas pour étudier un milieu 25 plus complexe tel qu'un tissu biologique. Or, il existe un besoin pour cartographier la structure des tissus biologiques composés de fibres tels que notamment le myocarde, les muscles et le cerveau. Cette structure joue un rôle crucial à la fois dans la fonction 30 mécanique (tissus musculaires) mais aussi électrique (cerveau, muscle, coeur) de ces tissus, et l'orientation des fibres dans l'espace est donc un paramètre très important à déterminer pour le diagnostic et l'exploration fonctionnelle de ces organes. 35 Par exemple, en imagerie cérébrale, il est très important d'identifier les faisceaux de fibres neuronales qui relient entre eux les différentes zones cérébrales. Actuellement, la seule technique capable de fournir une image tridimensionnelle de l'organisation des fibres est 5 l'Imagerie par Résonance Magnétique de tenseur de diffusion (IRM de diffusion). Cette technique très lente à mettre en oeuvre est utilisée pour l'exploration du cerveau adulte, mais reste trop limitée pour imager les organes en mouvements tels que le coeur. Enfin l'IRM n'est pas utilisée 10 pour imager le cerveau de très jeunes enfants, en particulier les bébés prématurés qui peuvent présenter des anomalies du développement cérébral malheureusement impossible à diagnostiquer avec les techniques actuelles. RESUME DE L'INVENTION 15 La présente invention a notamment pour but de pallier ces inconvénients. A cet effet, l'invention propose un procédé de cartographie de milieu fibreux, comprenant : (a) une étape de mesure au cours de laquelle on 20 fait émettre par un ensemble de transducteurs Tif, dans un champ d'observation appartenant à un milieu comportant des fibres, un nombre N d'ondes ultrasonores incidentes non focalisées (c'est-à-dire non focalisées dans le champ d'observation) 1 ayant des fronts d'ondes différents, et on 25 fait capter par les transducteurs Ti7 des signaux RFrawl,ii(t) représentatifs d'ondes ultrasonores réverbérées par le milieu respectivement à partir des ondes incidentes 1, (b) une étape de synthèse de données cohérentes au 30 cours de laquelle on détermine à partir des N ensembles de signaux RFrawl,ii(t) captés, pour un nombre M de points de focalisation fictifs Pk dans le champ d'observation, des signaux cohérents RFcoherentk,ii(t) correspondant aux signaux qui auraient été reçus par les transducteurs Ti7 si 35 une onde focalisée au point Pk avait été émise par lesdits transducteurs, (c)une étape de cartographie des fibres du milieu, au cours de laquelle on détermine une présence et une orientation de fibres en chaque point Pk, en comparant une cohérence spatiale entre les signaux cohérents RFcoherentk,ii(t) selon plusieurs directions. On peut ainsi, très rapidement et aisément, cartographier la structure des tissus biologiques composés de fibres tels que le myocarde et autres muscles, et le cerveau, grâce au fait que les signaux rétrodiffusés contiennent des informations sur la microstructure des tissus qui ne sont pas visible directement sur l'image échographique (B mode). C'est l'analyse de la cohérence spatiale qui révèle l'orientation des fibres car l'anisotropie des tissus se retrouve au niveau de la fonction de cohérence mesurée dans différentes directions. Dans divers modes de réalisation du procédé selon l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes : - au cours de l'étape (a), l'ensemble de transducteurs utilisé est un réseau bidimensionnel de transducteurs ; - au cours de l'étape (c), on détermine une intégrale d'une fonction de cohérence spatiale entre transducteurs dans plusieurs directions, et on détermine la direction des fibres comme étant une direction maximisant ladite intégrale ; - les ondes ultrasonores incidentes sont des ondes planes ayant différentes directions de propagation ; - les ondes ultrasonores incidentes sont des ondes divergentes (émises par la barrette ultrasonore comme si elles provenaient de différents points sources) ; - les ondes ultrasonores incidentes sont émises successivement ; - les ondes ultrasonores incidentes sont codées spatio-temporellement et émises simultanément, puis les ondes réverbérées sont captées simultanément puis séparées par décodage ; - on affiche une image des fibres détectées dans le milieu ; - on détermine une image échographique du champ d'observation et on affiche cette image échographique en superposition avec l'image des fibres ; - on détermine l'image échographique par formation de voies à partir des signaux cohérents déterminés à l'étape (b) ; - le milieu à imager est un tissu humain ou animal (notamment mammifère). Par ailleurs, l'invention a également pour objet un dispositif pour la mise en oeuvre d'un procédé de cartographie tel que défini ci-dessus, comprenant un ensemble de transducteurs Ti7 et des moyens de commande et traitement adaptés pour : (a) faire émettre par l'ensemble de transducteurs Tif, dans un champ d'observation d'un milieu comportant des fibres, un nombre N d'ondes ultrasonores incidentes 1 ayant des fronts d'ondes différents, et faire capter par les transducteurs Ti7 des signaux RFrawl,ii(t) représentatifs d'ondes ultrasonores réverbérées respectivement à partir des ondes incidentes 1, (b) déterminer à partir des N ensembles de signaux RFrawl,ii(t) captés, pour un nombre M de points de focalisation fictifs Pk dans le champ d'observation, des signaux cohérents RFcoherentk,ii(t) correspondant aux signaux qui auraient été reçus par les transducteurs Ti7 si une onde focalisée au point Pk avait été émise par lesdits transducteurs, (c)déterminer une présence et une orientation de fibres en chaque point Pk, en comparant des fonctions de 35 cohérence spatiale entre les signaux RFcoherentk,ii(t) selon plusieurs directions. Avantageusement, l'ensemble de transducteurs est bidimensionnel. BREVE DESCRIPTION DES DESSINS D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description suivante d'un de ses modes de réalisation, donné à titre d'exemple non limitatif, en regard des dessins joints. Sur les dessins : - la figure 1 est une vue schématique d'un dispositif pour la mise en oeuvre d'un procédé selon un mode de réalisation de l'invention, et - la figure 2 est un schéma bloc d'une partie du dispositif de la figure 1. DESCRIPTION DETAILLEE Sur les différentes figures, les mêmes références désignent des éléments identiques ou similaires. La figure 1 montre un exemple de dispositif d'imagerie fonctionnant par émission et réception d'ondes 20 ultrasonores de compression, par exemple de fréquences comprises entre 2 et 40 MHz. Le dispositif d'imagerie représenté sur la figure 1 est adapté pour réaliser une imagerie ultrasonore synthétique d'un champ d'observation la dans un milieu 25 fibreux 1, par exemple des tissus d'un patient, notamment un muscle (myocarde ou autre) ou un cerveau. Le dispositif d'imagerie comporte par exemple : - un réseau 2 de n transducteurs ultrasonores, par exemple un réseau à deux dimensions comprenant par exemple 30 quelques centaines de transducteurs et adapté pour réaliser une image tridimensionnelle (3D) du champ d'observation la ; - une baie électronique 3 ou similaire commandant le réseau de transducteurs 2 et adaptée pour acquérir les 35 signaux captés par ce réseau de transducteurs ; - un ordinateur 4 ou similaire pour commander la baie électronique 3 et visualiser les images ultrasonores obtenues à partir desdits signaux captés. Le réseau 2 de transducteurs peut par exemple être 5 un réseau matriciel plan s'étendant selon deux axes perpendiculaires X, Y, l'axe Z perpendiculaire aux axes X, Y désignant la direction de la profondeur dans le champ d'observation. Dans ce qui suit, les transducteurs seront notés Tif, i et j étant deux indices désignant le rang de 10 chaque transducteur respectivement selon les axes X et Y. Le réseau 2 de transducteurs peut notamment comprendre n1 transducteurs dans la direction X et n2 transducteurs dans la direction Y, avec n= ni*n2. La description qui suit sera faite en prenant pour exemple ce type de réseau 2 de 15 transducteurs, mais d'autres formes de de réseau de transducteurs sont également possibles dans le cadre de la présente invention. Comme représenté sur la figure 2, la baie électronique 3 peut comprendre par exemple : 20 - n convertisseurs analogique / digital 5 (A/Dii) connectés individuellement aux n transducteurs Ti7 du réseau 2 de transducteurs, - n mémoires tampon 6 (Big) respectivement connectées aux n convertisseurs analogique / digital 5, 25 - une unité centrale 8 (CPU) communiquant avec les mémoires tampon 6 et l'ordinateur 4, - une mémoire 9 (MEM) connectée à l'unité centrale 8, - un processeur numérique de signal 10 (DSP) 30 connecté à l'unité centrale 8. Ce dispositif permet de mettre en oeuvre un procédé de cartographie des fibres du milieu 1, qui inclut notamment les trois étapes suivantes, mises en oeuvre par l'unité centrale 8 assistée du processeur numérique de signal 9 : 35 a) Mesure (émission / réception et enregistrement des données brutes), b) Synthèse de données cohérentes, c) Analyse de l'orientation des fibres, d) Eventuellement, détermination d'une image du milieu en B mode et superposition avec la cartographie des fibres. Etape a: Mesure (émission / réception et enregistrement des données brutes) : Le réseau de transducteur est mis en contact avec le milieu 1 et un nombre N d'ondes incidentes ultrasonores est émis dans le milieu 1 par les transducteurs Tif (N peut être compris par exemple entre 2 et 100, notamment entre 5 et 10). Les ondes incidentes en question sont non focalisées (plus précisément, non focalisées dans le champ d'observation) et ont des respectivement fronts d'onde différents, c'est-à-dire des fronts d'onde de formes différentes et ou d'orientation différentes. Avantageusement, les ondes incidentes peuvent être des ondes planes d'inclinaisons toutes différentes, caractérisées par leurs angles d'inclinaison respectifs ax, CXy par rapport à l'axe Z, respectivement dans les plans (X, Z) et (Y, Z), ou encore des ondes divergentes émises comme si elles provenaient de différents points de l'espace. Les ondes incidentes sont généralement des impulsions de moins d'une microseconde, typiquement environ 1 cycle de l'onde ultrasonore à la fréquence centrale. Les tirs d'ondes incidentes peuvent être séparés les uns des autres par exemple d'environ 50 à 200 microsecondes. Chaque onde incidente rencontre dans le milieu 1 des diffuseurs qui réverbèrent l'onde incidente. L'onde ultrasonore réverbérée est captée par les transducteurs Tif du réseau. Le signal ainsi capté par chaque transducteur Tif provient de l'ensemble du milieu 1, puisque l'onde incidente n'est pas focalisée en émission. Les signaux réverbérés captés par les n transducteurs Tii sont alors numérisées par les convertisseurs analogique-digital correspondants A/Dij est mémorisé dans les mémoires tampon correspondantes Big. Les signaux ainsi 5 mémorisés dans les mémoires tampon après chaque tir incidents seront appelées ci-après données brutes RF ("RF" est un terme usuellement utilisé dans le domaine, simplement du fait de la fréquence des ultrasons utilisés). Ces données brutes RF consistent en une matrice de ni*n2 10 signaux temporels RFrawl,ii(t) captés respectivement par les transducteurs Ti7 après le tir 1 d'onde ultrasonore incidente. Après chaque tir 1 d'onde incidente, les signaux mémorisés dans les mémoires tampons Bi7 sont transférés 15 dans la mémoire 9 du processeur de signal 8 aux fins de traitement par ce processeur. A la fin de l'étape (a), la mémoire 9 contient donc N matrices de signaux bruts RF. On notera que les différentes ondes incidentes pourraient être codées spatio-temporellement, de façon à 20 permettre une émission simultanée de tout ou partie des ondes incidentes et une et réception également simultanée des ondes réverbérées, qui sont alors séparées par décodage avant mémorisation. 25 Etape b: Synthèse de données cohérentes RF A partir des N matrices de données brutes RF, un nombre M de matrices de données cohérentes synthétiques RF est calculé par le processeur 8, respectivement en M points Pk(x, y, z) du champ d'observation la (k étant un entier 30 compris entre 1 et M et x, y, z étant les coordonnées du point Pk sur les axes X, Y, Z). Chacune de ces M matrices de données cohérentes synthétiques RF comporte ni*n2 , signaux temporels RFcoherentkii(t) correspondant aux signaux qui seraient captés respectivement par les 35 transducteurs Ti7 si les transducteurs émettaient une onde incidente focalisée au point Pk. Les matrices de données RF cohérentes peuvent être obtenues par exemple en supposant une vitesse de propagation c homogène dans tout le milieu 1 pour les ondes de compression ultrasonores, selon le principe expliqué notamment dans le document EP2101191 ou dans l'article de Montaldo et al. « Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography » (IEEE Trans Ultrason Ferroelectr Freq Control 2009 Mar ; 56(3): 489-506). La direction de propagation de l'onde plane correspondant à chaque tir 1 étant connue, et la vitesse de propagation c étant connue, le processeur 8 peut calculer pour chaque point Pk le temps de propagation rec(1,k) de l'onde incidente 1 jusqu'au point Pk, et le temps de propagation - rec (1,k, i, j) de l'onde réverbérée depuis le point Pk vers le transducteur Tif, donc le temps de trajet total aller-retour r(1,k, 1, j) = 7 + 1 rec(1,P Pk) + - rec 1 f Pkf l f j). Le signal spatialement cohérent pour le transducteur Tij, correspondant au point de focalisation virtuel Pk, est alors calculé selon la formule : RFcoherenCi(t)=IB(1)RFrawhi(r(1, k, i, j)) ( 1 ) où B(1) est une fonction de pondération pour la 25 contribution de chaque tir 1 d'onde incidente (dans les cas courants, les valeurs B(1) peuvent être toutes égales à 1). Les matrices de données cohérentes RFcoherentk peuvent ensuite être éventuellement affinées en corrigeant les effets d'aberrations dans le milieu 1, par exemple comme 30 expliqué dans les documents susmentionnés EP2101191 ou Montaldo et al. Etape c: Analyse de l'orientation des fibres On détermine ensuite, pour chaque matrice RFcoherentk, une cohérence spatiale, indicatrice de la cohérence entre les signaux RFcoherentki.(t) pour un même point Pk. Cette cohérence spatiale peut être mesurée par exemple par une fonction de cohérence spatiale R(m) calculée à partir des corrélations des signaux ck(ij,pq) reçus sur des transducteurs ij et pq, en sommant toutes corrélations entre les paires de transducteurs distants de m éléments dans une direction donnée du plan (X,Y). T2 tu)=I(RFcoheren,j(t)-RFcoheren,j)(RFcoherenG(t)-RFcohereno,) T1 ( 2) Où RFcoherent,,,,, est une moyenne temporelle de RFcoherentk,ii, et Tl, T2 sont deux instants. En considérant uniquement des transducteurs alignés entre eux selon une même direction du plan (X, Y) et en 15 renumérotant ces transducteurs Tq, q allant de 1 à Q, ces corrélations peuvent s'écrire c(p, q) et on obtient : T2 c(p,q)=I(RFcoherene(t)-RFcoheren)(RFcoherentkw(t)-RFcoherenq) T1 (2bis) 20 R(m) = Q m ic(q,q) (3) Q rq_i c(q,q+m) Le théorème de Van Cittert Zernike établit l'allure de cette fonction R(m) dans un milieu de diffuseur aléatoire (donc isotrope) pour un faisceau monochromatique. R(m) est 25 la transformée de Fourier spatiale au carré de la tâche focale. Pour une tâche focale dont l'extension latérale est donnée par une fonction sin(ax)/x, R(m) est un triangle dont le sommet est en m=0 (autocorrélation) et qui s'annule en m=Q. 30 Pour un milieu non isotrope, un surcroit de cohérence spatiale est obtenu lorsque la direction d'alignement des transducteurs est alignée le long des fibres. L'intégrale Sk de cette fonction dans la direction d'alignement considérée du plan (X, Y) donne un paramètre de cohérence spatiale, qui est maximum dans la direction d'alignement des fibres. En calculant ce paramètre de cohérence spatiale dans plusieurs directions d'alignement des transducteurs, on peut connaître la direction produisant le paramètre de cohérence spatiale Sk maximum et donc en déduire la direction des fibres au point Pk. On notera que les fonctions de cohérence spatiales R(m) susmentionnées ou les paramètres de cohérence spatiale Sk pourraient être moyennées sur plusieurs ponts Pk voisins, donc dans un petit volume du champ d'observation autour d'un point d'intérêt. Un autre paramètre de cohérence spatiale possible est le critère de focalisation Ck qui donne le ratio entre l'énergie cohérente et l'énergie incohérente rétrodiffusée. Avec les notations ci-dessus, c'est-à-dire en numérotant de q=1 à Q les transducteurs alignés selon une même direction du plan (X,Y), on a : (IQ q=1RFcoherentke(t-tq 12 ) Qr iRFcoherent4(t-tql2) Ck - (4) Où tq est un retard permettant de remettre en phase tous les signaux RFcoherentkg(t). Comme dans le cas précédent, on calcule ce paramètre 25 de cohérence spatiale selon plusieurs directions pour chaque point Pk, et on détermine la direction des fibres comme la direction maximisant le paramètre Ck. On peut ainsi déterminer une cartographie tridimensionnelle des fibres du milieu dans le champ 30 d'observation la, de façon très rapide. Cette cartographie peut avantageusement être présentée à l'utilisateur du dispositif sous la forme d'images de plans de coupe du milieu 1, affichées par exemple sur l'écran de l'ordinateur 4. Eventuellement, ces images peuvent être calculées en reconstituant une continuité entre les fibres détectées aux différents points Pk. Etape d: formation d'image A partir des matrices RFcoherentk calculées à l'étape (b), on peut éventuellement former une image tridimensionnelle B mode du champ d'observation la, par formation de voie en réception, comme expliqué par exemple dans le document EP2101191 susmentionné. Cette image B mode peut être éventuellement superposée à la cartographie des fibres déterminée à l'étape (c), et on peut afficher sur l'écran de l'ordinateur des images de plans de coupe du champ d'observation montrant à la fois l'image B mode et les fibres superposées à cette image. REVENDICATIONS1. Procédé de cartographie de milieu fibreux, comprenant : (a) une étape de mesure au cours de laquelle on fait émettre par un ensemble (2) de transducteurs Tif, dans un champ d'observation (la) appartenant à un milieu (1) comportant des fibres, un nombre N d'ondes ultrasonores incidentes 1 non focalisées ayant des fronts d'ondes différents, et on fait capter par les transducteurs Ti7 des signaux RFrawl,ii(t) représentatifs d'ondes ultrasonores réverbérées par le milieu respectivement à partir des ondes incidentes 1, (b) une étape de synthèse de données cohérentes au cours de laquelle on détermine à partir des N ensembles de signaux RFrawl,ii(t) captés, pour un nombre M de points de focalisation fictifs Pk dans le champ d'observation, des signaux cohérents RFcoherentk,ii(t) correspondant aux signaux qui auraient été reçus par les transducteurs Ti7 si une onde focalisée au point Pk avait été émise par lesdits transducteurs, (c)une étape de cartographie des fibres du milieu (1), au cours de laquelle on détermine une présence et une orientation de fibres en chaque point Pk, en comparant une cohérence spatiale des signaux cohérents RFcoherentk,ii(t) selon plusieurs directions. 2. Procédé selon la revendication 1, dans lequel au cours de l'étape (a), l'ensemble (2) de transducteurs utilisé est un réseau bidimensionnel de transducteurs. 3. Procédé selon la revendication 2, dans lequel au cours de l'étape (c), on détermine une intégrale de fonctions de cohérence spatiale entre transducteurs dans plusieurs directions, et on détermine la direction des fibres comme étant une direction maximisant ladite intégrale. 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel les ondes ultrasonores incidentes sont des ondes planes ayant différentes directions de propagation. 5. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel les ondes ultrasonores incidentes sont des ondes divergentes. 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel on affiche une image des fibres détectées dans le milieu (1). 7. Procédé selon la revendication 6, dans lequel on détermine une image échographique du champ d'observation (la) et on affiche cette image échographique en superposition avec l'image des fibres. 8. Procédé selon la revendication 7, dans lequel on détermine l'image échographique par formation de voies à partir des signaux cohérents déterminés à l'étape (b). 9. Procédé selon l'une quelconque des revendications précédentes, dans lequel le milieu à imager est un tissu 20 humain ou animal. 10. Dispositif pour la mise en oeuvre d'un procédé de cartographie selon l'une quelconque des revendications précédentes, comprenant un ensemble (2) de transducteurs et des moyens de commande et traitement (8, 4) adaptés 25 pour : (a) faire émettre par l'ensemble (2) de transducteurs Tif, dans un champ d'observation (la) d'un milieu (1) comportant des fibres, un nombre N d'ondes ultrasonores incidentes 1 ayant des fronts d'ondes 30 différents, et faire capter par les transducteurs Ti7 des signaux RFrawl,ii(t) représentatifs d'ondes ultrasonores réverbérées respectivement à partir des ondes incidentes 1, (b) déterminer à partir des N ensembles de signaux RFrawl,ii(t) captés, pour un nombre M de points de 35 focalisation fictifs Pk dans le champ d'observation, dessignaux cohérents RFcoherentk,ii(t) correspondant aux signaux qui auraient été reçus par les transducteurs Ti7 si une onde focalisée au point Pk avait été émise par lesdits transducteurs, (c)déterminer une présence et une orientation de fibres en chaque point Pk, en comparant une cohérence spatiale entre les signaux cohérents RFcoherentk,ii(t) selon plusieurs directions. 11. Dispositif selon la revendication 10, dans 10 lequel l'ensemble (2) de transducteurs est bidimensionnel.