Brevet FR3038217A1 — publié le 20170106

Procédé d'insonification pour obtenir un champ prédéterminé d'ondes ultrasonores, et procédé de fabrication pour réaliser une lentille ultrasonore à ces fins. La présente invention est relative aux procédés d'insonification pour obtenir un champ prédéterminé d'ondes ultrasonores dans un milieu sensiblement homogène masqué par une barrière osseuse, et aux procédés de fabrication pour réaliser des lentilles ultrasonores à ces fins. Plus particulièrement, l'invention concerne un procédé d'insonification pour obtenir, à partir d'une sonde ultrasonore, un champ objectif prédéterminé d'ondes ultrasonores dans un milieu comprenant une partie interne sensiblement homogène masquée par au moins une barrière aberratrice, le procédé comprenant une étape de simulation au cours de laquelle on simule au moins une propagation d'ondes ultrasonores dans le milieu à partir d'une cartographie de propriétés de propagation des ondes ultrasonores dans le milieu. On entend ici par « partie sensiblement homogène » un milieu ayant des caractéristiques sensiblement homogènes pour la propagation des ondes ultrasonores. Il peut s'agir par exemple du cerveau d'un humain ou animal, auquel cas la barrière aberratrice est constituée par le crâne. Le document W02004019784A1 décrit un exemple d'un tel procédé, dans lequel, à partir de la simulation susmentionnée, on détermine des signaux individuels à émettre par chaque transducteur pour obtenir le champ objectif dans la partie interne du milieu. Ce procédé fonctionne bien mais nécessite en pratique un grand nombre de transducteurs (typiquement 500 à 1000) et une électronique complexe capable de piloter individuellement chaque transducteur. La présente invention a notamment pour but de 35 perfectionner encore ce type de procédé, notamment afin de 3038217 2 permettre sa mise en oeuvre par des moyens simplifiés et peu coûteux. A cet effet, l'invention propose un procédé d'insonification d'un milieu comprenant au moins une 5 barrière aberratrice et une partie interne sensiblement homogène masquée par ladite barrière aberratrice, pour obtenir un champ objectif prédéterminé d'onde ultrasonore en au moins une zone prédéterminée appartenant à ladite partie interne, en faisant émettre lesdites ondes 10 ultrasonores par une sonde ultrasonore placée dans une position prédéterminée par rapport audit milieu à l'opposé de ladite partie interne par rapport à ladite barrière aberratrice, le procédé d'insonification comprenant au moins les étapes suivantes : 15 (c)une étape de calcul de lentille au cours de laquelle on calcule, en utilisant un modèle du milieu comprenant une cartographie de propriétés de propagation des ondes ultrasonores obtenue à partir d'une imagerie réelle du milieu, une lentille ultrasonore qui, lorsqu'elle 20 est interposée entre la sonde ultrasonore et la barrière aberratrice, est adaptée pour générer des aberrations telles que lorsque la sonde ultrasonore émet une onde ultrasonore prédéterminée, ladite onde prédéterminée génère ledit champ objectif prédéterminé d'onde ultrasonore dans 25 ladite zone prédéterminée, (d)une étape de réalisation de la lentille calculée à l'étape (c) de calcul de lentille, (e)une étape d'émission au cours de laquelle on fait émettre ladite onde ultrasonore prédéterminée par 30 ladite sonde ultrasonore au travers de la lentille ultrasonore disposée dans ladite position prédéterminée. Grâce à ces dispositions, on peut obtenir le champ objectif d'ondes ultrasonores de façon très simple, éventuellement avec un nombre réduit de transducteurs 35 (éventuellement un seul). En particulier, le procédé 3038217 3 d'insonification selon l'invention ne requiert pas l'émission d'une forme d'onde particulièrement complexe à partir de la sonde ultrasonore, puisque la sonde peut émettre une onde prédéterminée particulièrement simple, par 5 exemple une onde focalisée à une distance prédéterminée. On peut ainsi, à partir d'une imagerie réelle du milieu telle qu'une imagerie rayons X (CT) ou IRM, calculer rapidement par des moyens informatiques la lentille ultrasonore, puis réaliser cette lentille également 10 rapidement, par exemple par des moyens de fabrication informatisés tels que l'impression tridimensionnelle ou le fraisage numérique. On peut alors insonifier le milieu de façon voulue et personnalisée, avec une sonde ultrasonore standard et des moyens de commande également standard. La 15 lentille est réutilisable un grand nombre de fois pour le même milieu 3 (par exemple pour le même patient). La sonde ultrasonore et ses moyens de commande sont réutilisables pour plusieurs milieux (par exemple pour plusieurs patients), sans changement ni reprogrammation. 20 Le procédé d'insonification est ainsi rendu très facile d'utilisation, et de façon peu coûteuse. Dans des modes de réalisation préférés du procédé d'insonification selon l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des 25 dispositions suivantes : - l'étape (c) de calcul de lentille est précédée par une étape d'imagerie (a) au cours de laquelle est réalisée ladite cartographie ; - l'étape d'imagerie (a) est réalisée par imagerie 30 rayons X ou par IRM ; - ladite cartographie est réalisée à partir de mesures de densité osseuse ; - l'étape (d) de réalisation de la lentille est effectuée par une méthode choisie parmi l'impression 35 tridimensionnelle de la lentille et le fraisage à commande 3038217 4 numérique d'au moins un bloc de matériau pour former la lentille ; - l'étape (d) de réalisation de la lentille comprend les sous-étapes suivantes : 5 (dl) une sous-étape de réalisation d'au moins un moule effectuée par une méthode choisie parmi l'impression tridimensionnelle dudit au moins un moule et le fraisage à commande numérique d'au moins un bloc de matériau pour former ledit au moins un moule, 10 (d2) une sous-étape de moulage dans laquelle on moule ladite lentille ou au moins un composant de ladite lentille dans ledit au moins un moule ; - l'étape (c) de calcul de lentille est précédée par une étape de simulation (b) au cours de laquelle : 15 (p') on calcule une propagation d'une onde ultrasonore prédéterminée dans le milieu depuis la sonde ultrasonore disposée à ladite position prédéterminée, jusqu'à une position choisie de la lentille et on détermine des premiers temps d'arrivée to(x,y) de ladite onde 20 ultrasonore prédéterminée au niveau de ladite position choisie de la lentille, (b2) on simule une rétro-propagation dudit champ objectif prédéterminé d'onde ultrasonore dans le milieu (3), depuis la zone prédéterminée jusqu'à ladite position 25 choisie de la lentille, (b3) on détermine des deuxième temps d'arrivée tl(x,y) de ladite rétro-propagation dudit champ objectif prédéterminé d'onde ultrasonore à ladite position choisie de la lentille, 30 (b4) on détermine une loi de délai At(x,y) à ladite position choisie de la lentille, égale à la somme des premiers temps d'arrivée to(x,y) et des deuxièmes temps d'arrivée tl(x,Y), et au cours de l'étape (c) de calcul de lentille, on 35 calcule la lentille ultrasonore à partir de ladite loi de 3038217 5 délai At(x,y), de façon que lorsque la sonde ultrasonore émet ladite onde ultrasonore prédéterminée au travers de ladite lentille ultrasonore, ladite onde prédéterminée reproduit le champ d'onde voulu à l'endroit cible après 5 traversée de ladite barrière aberratrice ; - les premiers temps d'arrivée to(x,y) et les deuxièmes temps d'arrivée tl(x,y) sont déterminés en différents points de contrôle M(x,y) appartenant à une surface de contrôle prédéterminée située à ladite position 10 choisie de la lentille ; - la surface de contrôle est située entre la sonde ultrasonore et la barrière aberratrice ; - à l'étape (c) de calcul de lentille, on calcule une épaisseur e(x,y) de lentille ultrasonore, en chaque 15 point L(x,y) de la lentille correspondant à un point de contrôle M(x,y), par la formule suivante : e(x,y)= eo + At(x,y)/(1/c-1/c1)(3) où : - c est la célérité de l'onde ultrasonore dans le milieu hors de la barrière aberratrice, 20 - cl est la célérité de l'onde ultrasonore dans la lentille ; - eo est une constante telle que e(x,y) ait une valeur positive et supérieure à une épaisseur minimale prédéterminée en tout point de la lentille ; 25 - l'onde ultrasonore est monochromatique et à l'étape (c) de calcul de lentille (9), on calcule une épaisseur e(x,y) de lentille, en chaque point L(x,y) de la lentille correspondant à un point de contrôle M(x,y), par la formule suivante : 30 e(x,y)= eo + At(x,y)/(1/c-1/cl), modulo T/(1/c-1/c1). où : - c est la célérité de l'onde ultrasonore dans le milieu hors de la barrière aberratrice, - cl est la célérité de l'onde ultrasonore dans la lentille ; - eo est un nombre réel tel que e(x,y) ait une 3038217 6 valeur positive et supérieure à une épaisseur prédéterminée en tout point de la lentille, - T est la période de l'onde ultrasonore ; - l'étape (c) de calcul de lentille tient compte 5 des temps de parcours de ladite onde ultrasonore prédéterminée dans la lentille ultrasonore et des angles d'incidences de ladite onde ultrasonore prédéterminée à la surface de la lentille ultrasonore ; - l'étape (c) de calcul de lentille tient compte 10 de l'ensemble de la propagation de ladite onde ultrasonore prédéterminée, y compris des réverbérations entre la sonde ultrasonore, la lentille ultrasonore et la barrière aberratrice, ce qui permet d'optimiser le rendement de la transmission d'énergie ultrasonore dans la zone 15 prédéterminée ; - l'étape (c) de calcul de lentille se fait itérativement afin d'optimiser le rendement de la transmission d'énergie ultrasonore dans la zone prédéterminée, à l'aide d'un modèle analytique ou 20 numérique ; - le milieu dans lequel on fait propager les ondes ultrasonores est une tête humaine ou animale, la barrière aberratrice étant un crâne et la partie interne étant un cerveau ; 25 - ledit champ prédéterminé d'onde ultrasonore est focalisé dans ladite partie interne en au moins une zone prédéterminée. Par ailleurs, l'invention a également pour objet un procédé de fabrication pour réaliser une lentille 30 ultrasonore aux fins d'insonification d'un milieu comprenant au moins une barrière aberratrice et une partie interne sensiblement homogène masquée par ladite barrière aberratrice, pour obtenir un champ objectif prédéterminé d'onde ultrasonore en au moins une zone prédéterminée 35 appartenant à ladite partie interne, en faisant émettre 3038217 7 lesdites ondes ultrasonores par une sonde ultrasonore placée dans une position prédéterminée par rapport audit milieu à l'opposé de ladite partie interne par rapport à ladite barrière aberratrice, 5 le procédé de fabrication comprenant au moins une étape de calcul de lentille (c) au cours de laquelle on calcule, en utilisant un modèle du milieu comprenant une cartographie de propriétés de propagation des ondes ultrasonores obtenue à partir d'une imagerie réelle du milieu, une lentille 10 ultrasonore qui, lorsqu'elle est interposée entre la sonde ultrasonore et la barrière aberratrice, est adaptée pour générer des aberrations telles que lorsque la sonde ultrasonore émet une onde ultrasonore prédéterminée, ladite onde prédéterminée génère ledit champ objectif prédéterminé 15 d'onde ultrasonore dans ladite zone prédéterminée. Dans des modes de réalisation préférés du procédé de fabrication selon l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes : 20 - le procédé de fabrication comporte en outre une étape (d) de réalisation de la lentille calculée à l'étape (c) de calcul de lentille ; - l'étape (d) de réalisation de la lentille est effectuée par une méthode choisie parmi l'impression 25 tridimensionnelle de la lentille et le fraisage à commande numérique d'au moins un bloc de matériau pour former la lentille ; - l'étape (d) de réalisation de la lentille comprend les sous-étapes suivantes : 30 (dl) une sous-étape de réalisation de moule effectuée par une méthode choisie parmi l'impression tridimensionnelle dudit au moins un moule et le fraisage à commande numérique d'au moins un bloc de matériau pour former ledit au moins un moule, 35 (d2) une sous-étape de moulage dans laquelle on 3038217 8 moule ladite lentille ou au moins un composant de ladite lentille (9) dans ledit au moins un moule. D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description 5 suivante d'une de ses formes de réalisation, donnée à titre d'exemple non limitatif, en regard des dessins joints. Sur les dessins : - la figure 1 est une vue d'ensemble schématique d'un exemple dispositif de génération d'ondes ultrasonores 10 mettant en oeuvre l'invention, comprenant une lentille ultrasonore, - la figure 2 est un schéma fonctionnel du dispositif de la figure 1, - la figure 3 est un schéma illustrant le mode de 15 calcul de la lentille acoustique utilisable dans le dispositif de la figure 1, - et la figure 4 est un schéma illustrant le fonctionnement du dispositif de la figure 1. Sur les différentes figures, les mêmes références 20 désignent des éléments identiques ou similaires. Le dispositif 1 de génération d'ondes ultrasonores représenté sur la figure 1 est adapté pour générer des champs d'ondes ultrasonores prédéterminés dans un milieu 3, par exemple la tête d'un patient P, ou une autre partie du 25 corps. Le milieu 3 comporte au moins une barrière aberratrice 8 (ici le crâne du patient P) et une partie interne sensiblement homogène 7 (ici le cerveau du patient P) masquée par ladite barrière aberratrice 8. Plus généralement, le milieu sensiblement homogène 30 7 pourrait être tout milieu tissulaire d'un patient humain ou de tout autre vertébré, milieu tissulaire ayant des caractéristiques sensiblement homogènes pour la propagation des ondes ultrasonores, tandis que la barrière aberratrice pourrait être toute autre barrière osseuse. A titre 35 d'exemple alternatif, le milieu tissulaire pourrait 3038217 9 éventuellement être constitué par le coeur ou l'ensemble coeur-poumons du patient P, auquel cas la barrière osseuse serait constituée par la cage thoracique. Le dispositif 1 est destiné à générer des ondes 5 ultrasonores dans le cerveau 7 du patient P (ou plus généralement dans la partie interne 7 du milieu 3), à des fréquences par exemple de l'ordre de 0,1 à 10 MHz, notamment de 0,2 à 3 MHz, depuis l'extérieur du crâne 8. Cette génération d'ondes ultrasonores peut être 10 destinée par exemple : - à réaliser une image ou une série d'images échographiques du cerveau, qu'il s'agisse d'imagerie statique ou fonctionnelle, notamment d'imagerie Doppler permettant de visualiser les écoulements sanguins ou 15 d'imagerie thermique permettant de visualiser les échauffements provoqués par un traitement par hyperthermie, - et/ou de réaliser un traitement par hyperthermie, notamment pour détruire des tumeurs, coaguler des hémorragies, activer localement des médicaments thermo20 activables, - et/ou ouvrir localement la barrière hématoencéphalique afin de faire diffuser localement un médicament précédemment injecté par voie intraveineuse - et/ou activer/inhiber le cerveau par 25 neuromodulation ultrasonore. Dans tous les cas, il est nécessaire de pouvoir générer, avec le plus de précision possible, un ou plusieurs champ(s) objectif(s) prédéterminé(s) d'ondes ultrasonores dans le cerveau 7 du patient P, par exemple 30 pour focaliser en un ou plusieurs points du cerveau les ondes ultrasonores émises par une sonde ultrasonore 2, ou pour générer des champs d'ondes plus complexes. Comme représenté sur la figure 4, le champ d'ondes objectif peut par exemple être un champ d'ondes focalisé en une tache 35 focale F dans le cerveau 7, ou éventuellement plusieurs 3038217 10 taches focales. Comme représenté sur la figure 1, la sonde ultrasonore 2 peut être par exemple intégrée dans un casque 4 positionné de façon prédéterminée sur la tête 3 du 5 patient P, ou bien ladite sonde ultrasonore 2 peut être portée par tout autre système de positionnement connu. Dans tous les cas, la sonde ultrasonore 2 est positionnée à l'extérieur du crâne 8 du patient dans une position prédéterminée. 10 Une lentille acoustique 9, dite ici également lentille ultrasonore, également intégrée au casque 4, est interposée entre la sonde ultrasonore 2 et le crâne 8 du patient P. Un gel ou un liquide, éventuellement contenu dans une poche souple (non représentée), peut être 15 interposé entre la lentille ultrasonore 9 et le crâne 8 et/ou entre la sonde ultrasonore 2 et la lentille ultrasonore 9, de façon à assurer une bonne transmission des ondes ultrasonores. Ce gel ou liquide sera ici appelé milieu externe et considéré comme faisant partie du milieu 20 de propagation 3 des ondes ultrasonores. La lentille acoustique 9 peut être réalisée en tout matériau dans lequel la célérité cl des ondes acoustiques de compression est différente de la célérité c desdites ondes dans le cerveau (c est sensiblement la célérité du 25 son dans l'eau, soit environ 1480 m/s). La lentille ultrasonore peut être réalisée par exemple en matière plastique, notamment en polydiméthylsiloxane également appelé PDMS (c1=1030 m/s) ou en polyméthylpentène, connu sous l'appellation TPX® (c1=2090 m/s). 30 Eventuellement, certaines parties de la lentille ultrasonore peuvent être réalisées en un matériau atténuant, réfléchissant ou absorbant les ondes ultrasonores, notamment pour atténuer localement le champ d'ondes ultrasonore produit dans le cerveau 7 ou le crâne 35 8, par exemple dans les zones ne nécessitant pas de 3038217 11 traitement ou pour limiter l'échauffement du crâne. La sonde ultrasonore 2 peut être commandée par un circuit électronique 5 (B) lui-même éventuellement contrôlé par un ordinateur 6 (PC). 5 Comme représenté sur la figure 2, la sonde ultrasonore 2 peut comprendre un ou plusieurs transducteur(s) ultrasonores 11-18. Ces transducteurs sont au nombre de 8 dans l'exemple représenté, mais pourraient être plus nombreux notamment dans des applications 10 d'imagerie, ou moins nombreux. Eventuellement, on peut utiliser un seul transducteur T1, notamment lorsqu'on cherche uniquement à focaliser l'onde ultrasonore notamment aux fins de traitement par effet thermique ou mécanique. La sonde ultrasonore 2 peut éventuellement comporter des 15 transducteurs destinés à l'imagerie ultrasonore et d'autres destinés au traitement par hyperthermie, comme expliqué dans le document W02004019784A1 susmentionné. Comme représenté sur la figure 2, le circuit électronique 5 qui commande la sonde ultrasonore 2 peut 20 comprendre : - une unité centrale électronique 10 (CPU) commandée par l'ordinateur 7, - au moins une mémoire centrale 11 (MEM) reliée à l'unité centrale électronique 10, 25 - un convertisseur analogique-digital A/D1-A/D8 relié à chaque transducteur 11-18 et généralement associé à une mémoire tampon (non représentée), - éventuellement, un circuit spécifique de traitement de signal 12 (DSP). 30 Pour mettre en oeuvre le dispositif décrit précédemment, on conçoit et fabrique la lentille ultrasonore 9 de façon qu'elle permette d'obtenir le champ d'ondes objectif souhaité, par exemple en suivant les étapes suivantes : 35 3038217 12 (a) Détermination d'un modèle tridimensionnel du milieu 3 : On détermine préalablement un modèle tridimensionnel du milieu 3 comprenant une cartographie de 5 propriétés de propagation des ondes ultrasonores dans le milieu. Cette étape comprend généralement une opération d'imagerie, notamment imagerie rayons X ou IRM, notamment afin de déterminer la densité osseuse de la barrière 10 aberratrice 8 en chaque point, ce qui permet de déterminer les paramètres de propagation des ondes ultrasonores en chaque point, comme expliqué dans le document W02004019784A1 susmentionné (notamment densité massique p, célérité c des ondes ultrasonores et/ou absorption 1 15 desdites ondes ultrasonores, en tout point de la barrière aberratrice 8). Le modèle tridimensionnel peut être par exemple chargé dans l'ordinateur 6. Un utilisateur peut positionner sur l'image du milieu 3 le champ d'onde ultrasonore 20 souhaité, par exemple la ou les taches focale(s) F souhaitée(s) (Figure 4). Ce positionnement peut s'effectuer par exemple en visualisant sur l'image du milieu 3 une lésion à traiter ou autre zone d'intérêt, et en repérant cette zone sur l'image à l'aide de l'interface utilisateur 25 de l'ordinateur 6 (par exemple souris, écran tactile ou autre). Dans ce qui suit, on fera référence à « la zone d'intérêt F » par mesure de simplicité, pour désigner une ou plusieurs taches focales ou une ou des zones à champ d'onde objectif plus complexe. 30 (b) simulation : A partir du modèle tridimensionnel susmentionné, on peut simuler la propagation d'une onde ultrasonore à travers la barrière aberratrice pour ensuite déterminer 35 les délais ou déphasages qui devront être induits par la 3038217 13 lentille, puis en déduire la forme à donner à la lentille. Cette étape peut être réalisée par l'ordinateur 6 susmentionné ou un autre ordinateur. Cette étape de simulation peut comporter les sous-5 étapes suivantes : (b1) Détermination des temps d'arrivée sur une surface de contrôle (Si - S'1) : Comme représenté en Si et S'1 sur la figure 3, il est possible de déterminer des temps d'arrivée to(x,y) 10 d'une onde ultrasonore prédéterminée 13 émise par la sonde ultrasonore 2, en différents points de contrôle M(x,y), appartenant préférentiellement à une même surface de contrôle 14 disposée entre la sonde 2 et la zone d'intérêt F susmentionnée. Cette surface de contrôle 14 est située 15 entre la barrière aberratrice 8 et la sonde ultrasonore 2, a la position qui doit être occupée par la lentille acoustique 9 susmentionnée. La surface de contrôle 14 peut par exemple être un plan XY dont les axes X, Y sont perpendiculaires à l'axe Z disposé dans le sens de la 20 profondeur. Cette surface de contrôle pourrait être toute autre surface, plane ou courbe. Ces temps d'arrivée to(x,y) sont aisément déterminés par l'ordinateur 6 ou autre (à une constante de temps Tt° près quelconque commune à toute la surface de 25 contrôle 14), puisque la partie externe du milieu (à l'extérieur de la barrière aberratrice) est considérée comme homogène pour la propagation des ondes ultrasonores. Dans le cas simple et le plus courant où la sonde est sphérique de rayon R, cette détermination peut être 30 faite très simplement sous la forme : to(x,y) = T'a- d(M)/c (2) où : - d(M) est la distance entre le point de contrôle M(x,y) et le centre de courbure de la sonde, - c est la célérité de l'onde ultrasonore dans la 35 partie externe du milieu 3. 3038217 14 Plus généralement, les temps d'arrivée to(x,y) peuvent être déterminés par calcul de propagation de l'onde ultrasonore entre la sonde et les points de contrôle ou par toute méthode connue de propagation d'onde, ce qui est 5 simple en milieu homogène (on peut alors effectuer une propagation dans un modèle simplifié permettant un calcul analytique). Les temps d'arrivée to(x,y) seront dits ci-après premiers temps d'arrivée. 10 (b2) simulation de la rétro-propagation d'onde ultrasonore voulue dans la zone d'intérêt dans le modèle tridimensionnel (52): Comme représenté en S2 sur la figure 3, au cours de cette sous-étape b2, on simule une rétro-propagation du 15 champ objectif prédéterminé d'onde ultrasonore dans le milieu 3, depuis la zone d'intérêt F dans le milieu 3 jusqu'à la surface de contrôle 14 susmentionnée. Cette simulation peut être effectuée par exemple par l'ordinateur 6, par exemple en utilisant une équation 20 d'ondes telle que l'équation (1) ci-après : _ a - ( - (1+2 (r)--) P (i)V 1 1 2P(i'm at - V où i désigne un vecteur p(i,t) P (i) ) c (i)2 at 2 =s (i,t) (1), de position _ du point considéré, p désigne la pression et S désigne les signaux ultrasonores générés par une source ultrasonore éventuellement présente 25 au point considéré. La propagation des ondes ultrasonores dans la barrière aberratrice 8 et la partie interne 7 peut être simulée dans l'ordinateur par différences finies, en discrétisant l'équation (1) ci-dessus. Cette simulation 30 peut également être effectuée par éléments finis, ou par une méthode de diffraction impulsionnelle, ou toute autre méthode connue. (b3) Détermination des deuxièmes temps d'arrivée 3038217 15 sur la surface de contrôle 14 (5'2) : Comme représenté en S'2 sur la figure 3, la simulation susmentionnée permet de déterminer des temps d'arrivée ti (x, y) de l'onde ultrasonore 13c, 13b 5 représentative de la rétro-propagation du champ d'onde objectif prédéterminé auxdits points de contrôle M(x,y) appartenant à la surface de contrôle 14. Les temps d'arrivée tl(x,y) sont comptés par exemple à partir de l'instant d'émission de l'onde ultrasonore 13c 10 dans ladite zone F, x et y étant respectivement les coordonnées de chaque point de contrôle M considéré, sur les axes X, Y. Les temps d'arrivée tl(x,y) sont dits ci-après deuxièmes temps d'arrivée. 15 (b4) Détermination de la loi de délai de la lentille (53) : Comme représenté en S3 sur la figure 3, on détermine une loi de délai At(x,y) ou loi de retard auxdits différents points de contrôle M(x,y), égale à la somme des 20 premiers temps d'arrivée to(x,y) et des deuxièmes temps d'arrivée tl(x,Y)- (c)Calcul de la lentille ultrasonore (54) : 25 A partir de la loi de délai susmentionnés At(x,y), l'ordinateur 6 ou autre calcule ensuite le profil de la lentille ultrasonore 9 qui, lorsqu'elle est interposée entre la sonde ultrasonore 2 et la barrière aberratrice 8, est adaptée pour introduire ledit délai At(x,y), de façon 30 que lorsque la sonde ultrasonore 2 émet ladite onde ultrasonore prédéterminée 13a, ladite onde prédéterminée produit après traversée de la lentille puis de la barrière aberratrice le champ d'onde ultrasonore prédéterminé dans la zone d'intérêt . 35 Plus précisément, on peut calculer une épaisseur e 3038217 16 de la lentille ultrasonore 9, en chaque point de ladite lentille correspondant à un point de contrôle M(x,y) de la surface de contrôle 14, pour que l'épaisseur e de lentille introduise une correction de temps de propagation d'onde 5 ultrasonore correspondant à ladite loi de délai, modulo la période T=1/f de l'onde ultrasonore (f est la fréquence centrale de l'onde ultrasonore) et / ou à une constante de temps To près quelconque commune à toute la surface de contrôle. 10 Par exemple, on peut calculer l'épaisseur e(x,y) de lentille ultrasonore, en chaque point L(x,y) de la lentille correspondant à un point de contrôle M(x,y), par la formule (3) suivante : e(x,y)= At(x,Y)/(1/c-1/c1)+ eo (3) . 15 La constante e0 est un nombre réel quelconque, tel que e(x,y) ait une valeur positive et minimale en tout point de la lentille. Dans le cas de l'émission d'une onde monochromatique, on peut calculer l'épaisseur e(x,y) de 20 lentille ultrasonore, en chaque point L(x,y) de la lentille correspondant à un point de contrôle M(x,y), par la formule (4) suivante : e(x,y)= At(x,y)/(1/c-1/c1)+ eo, modulo T/(1/c-1/c1) (4) . Modulo T/(1/c-1/c1) signifie que e(x,y)= 25 At(x,y)/(1/c-1/c1)+ e0 +n.T/(1/c-1/c1), où n est un entier positif ou négatif tel que e(x ,y) ait une valeur positive. L'entier n peut être commun à toute la lentille 9 (lentille simple) ou éventuellement être différent suivant les régions de la lentille 9 (lentille de Fresnel). 30 On notera que l'étape (c) de calcul de la lentille peut tenir compte au moins des temps de parcours de l'onde ultrasonore prédéterminée dans la lentille ultrasonore 9 et des angles d'incidences de ladite onde ultrasonore prédéterminée à la surface de la lentille ultrasonore 9. 35 L'étape (c) de calcul de lentille peut aussi tenir compte 3038217 17 de l'ensemble de la propagation de ladite onde ultrasonore prédéterminée, y compris des réverbérations entre la sonde ultrasonore 2, la lentille ultrasonore 9 et la barrière aberratrice 8. Eventuellement l'étape (c) de calcul de 5 lentille peut se faire itérativement afin d'optimiser le rendement de la transmission d'énergie ultrasonore dans la zone prédéterminée F, à l'aide d'un modèle analytique ou numérique. 10 (d) Réalisation de la lentille ultrasonore (55) : Une fois calculé le profil de la lentille 9 pour le patient P (ou plus généralement pour une application donnée), celle-ci peut être fabriquée, de préférence par une méthode de fabrication rapide. 15 On peut en particulier utiliser des méthodes de fabrication couramment utilisées en prototypage rapide, par exemple : - par impression tridimensionnelle de la lentille 9, 20 - par fraisage à commande numérique d'un bloc de matériau pour former la lentille 9, - par un procédé en deux sous-étapes comprenant : (dl) une sous-étape de réalisation de moule (non représenté) effectuée par une méthode choisie parmi 25 l'impression tridimensionnelle du moule et le fraisage à commande numérique d'un bloc de matériau pour former le moule, (d2) une sous-étape de moulage dans laquelle on moule ladite lentille 9 dans ledit moule. 30 (e)Utilisation : Une fois la lentille 9 réalisée, on peut la fixer dans le casque 4 susmentionné ou autre support, et ainsi la positionner dans la position prédéterminée qu'elle doit 35 occuper sur le crâne 8 du patient P. Un opérateur peut 3038217 18 alors faire émettre l'onde prédéterminée 13a par la sonde ultrasonore 2 au travers de la lentille ultrasonore 9. Comme représenté sur la figure 4, la lentille 9 introduit des aberrations qui déforment le front d'onde 13b 5 en amont du crâne 8, de façon qu'après passage au travers du crâne 8 ou autre barrière aberratrice, le front d'onde 13c corresponde à l'onde théorique 13 susmentionnée. On obtient ainsi le champ d'ondes objectif, par exemple une focalisation dans la tache focale F. 10 L'onde ainsi émise dans le cerveau ou autre partie interne 7 peut servir à traiter la zone d'intérêt F comme expliqué précédemment. De plus, on peut éventuellement imager le cerveau au ou autre partie interne 7 au moins au voisinage de la zone d'intérêt F, lorsque la sonde 15 ultrasonore 2 comporte des transducteurs ultrasonore d'imagerie, puisque la lentille ultrasonore 9 compense les aberrations dues au crâne ou autre barrière aberratrice 8, aussi bien en émission qu'en réception. Le procédé selon l'invention peut en outre 20 comporter une étape facultative (f) de vérification du positionnement du dispositif, au cours de laquelle on enregistre sur la sonde ultrasonore les échos ultrasonores qui se sont réfléchis sur ladite barrière aberratrice et on les compare aux mêmes signaux réfléchis simulés en 25 utilisant le modèle susmentionné du milieu 3. REVENDICATIONS1. Procédé de fabrication pour réaliser une lentille ultrasonore aux fins d'insonification d'un milieu comprenant au moins un é barrière aberratrice (8) et une partie interne (7) sensiblement homogène masquée par ladite barrière aberratrice (8), pour obtenir un champ objectif prédéterminé d'onde ultrasonore en au moins une zone prédéterminée (F) appartenant à ladite partie interne (7), en faisant émettre lesdites ondes ultrasonores par une sonde ultrasonore (2) placée dans une position prédéterminée par rapport audit milieu (3) à l'opposé de ladite partie interne (7) par rapport à ladite barrière aberratrice (8), le procédé de fabrication comprenant au moins une étape de calcul de lentille (c) au cours de laquelle on calcule, en utilisant un modèle du milieu comprenant une cartographie de propriétés de propagation des ondes ultrasonores obtenue à partir d'une imagerie réelle du milieu (3), une lentille ultrasonore (9) qui, lorsqu'elle est interposée entre la sonde ultrasonore (2) et la barrière aberratrice (8), est adaptée pour générer des aberrations telles que lorsque la sonde ultrasonore (2) émet une onde ultrasonore prédéterminée (13a), ladite onde prédéterminée génère ledit champ objectif prédéterminé d'onde ultrasonore dans ladite zone prédéterminée (F). 2. Procédé de fabrication selon la revendication 1, comportant en outre une étape (d) de réalisation de la lentille (9) calculée à l'étape (c) de calcul de lentille. 3. Procédé de fabrication selon la revendication 2, dans lequel l'étape (d) de réalisation de la lentille (9) est effectuée par une méthode choisie parmi l'impression tridimensionnelle de la lentille (9) et le fraisage à commande numérique d'au moins un bloc de matériau pour former la lentille (9). 3038217 20 4. Procédé de fabrication selon la revendication 2, dans lequel l'étape (d) de réalisation de la lentille (9) comprend les sous-étapes suivantes : (dl) une sous-étape de réalisation de moule 5 effectuée par une méthode choisie parmi l'impression tridimensionnelle dudit au moins un moule et le fraisage à commande numérique d'au moins un bloc de matériau pour former ledit au moins un moule, (d2) une sous-étape de moulage dans laquelle on 10 moule ladite lentille (9) ou au moins un composant de ladite lentille (9) dans ledit au moins un moule. 5. Procédé de fabrication selon l'une quelconque des revendications précédentes, dans lequel l'étape (c) de calcul de lentille est précédée par une étape d'imagerie 15 (a) au cours de laquelle est réalisée ladite cartographie. 6. Procédé de fabrication selon la revendication 5, dans lequel l'étape d'imagerie (a) est réalisée par imagerie rayons X ou par IRM. 7. Procédé de fabrication selon la revendication 6, 20 dans lequel ladite cartographie est réalisée à partir de mesures de densité osseuse. 8. Procédé de fabrication selon l'une quelconque des revendications précédentes, dans lequel l'étape (c) de calcul de lentille est précédée par une étape de simulation 25 (b) au cours de laquelle (b1) on calcule une propagation d'une onde ultrasonore prédéterminée (13a) dans le milieu depuis la sonde ultrasonore (2) disposée à ladite position prédéterminée, jusqu'à une position choisie de la lentille 30 (9) et on détermine des premiers temps d'arrivée to(x,y) de ladite onde ultrasonore prédéterminée au niveau de ladite position choisie de la lentille (9), (b2) on simule une rétro-propagation dudit champ objectif prédéterminé d'onde ultrasonore dans le milieu 35 (3), depuis la zone prédéterminée (F) jusqu'à ladite 3038217 21 position choisie de la lentille (9), (b3) on détermine des deuxième temps d'arrivée ti(x,y) de ladite rétro-propagation dudit champ objectif prédéterminé d'onde ultrasonore à ladite position choisie de la lentille (9), (b4) on détermine une loi de délai ,At(x,y) à ladite position choisie de la lentille (9), égale à la somme des premiers temps d'arrivée to(x,y) et des deuxièmes temps d'arrivée ti(xfY), 10 et au cours de l'étape (c) de calcul de lentille, on calcule la lentille ultrasonore (9) à partir de ladite loi de délai L\t(x,y), de façon que lorsque la sonde ultrasonore (2) émet ladite onde ultrasonore prédéterminée (13a) au travers de ladite lentille ultrasonore, ladite onde 15 prédéterminée reproduise le champ d'ondes objectif dans ladite zone F, après traversé de la barrière aberratrice 8. 9. Procédé de fabrication selon la revendication 8, dans lequel les premiers temps d'arrivée to(x,y) et les deuxièmes temps d'arrivée ti(x,y) sont déterminés en 20 différents points de contrôle M(x,y) appartenant à une surface de contrôle prédéterminée (14) située à ladite position choisie de la lentille (9). 10.Procédé de fabrication selon la revendication 9, dans lequel la surface de contrôle (14) est située entre la 25 sonde ultrasonore (2) et la barrière aberratrice (8). 11. Procédé de fabrication selon l'une quelconque des revendications 9 et 10, à l'étape (c) de calcul de lentille (9), on calcule une épaisseur e(x,y) de lentille ultrasonore, en chaque point L(x,y) de la lentille 30 correspondant à un point de contrôle M(x,y), par la formule suivante : e(x,y)- At(x,y)/(1/c-1/c1)+ eo où : - c est la célérité de l'onde ultrasonore dans le milieu (3) hors de la barrière aberratrice (8), 35 - cl est la célérité de l'onde ultrasonore 3038217 22 dans la lentille ; eo est un nombre réel tel que e(x,y) ait une valeur positive et supérieure à une épaisseur prédéterminée en tout point de la lentille. 12.Procédé de fabrication selon l'une quelconque des revendications 9 et 10, dans lequel l'onde ultrasonore est monochromatique et à l'étape (c) de calcul de lentille (9), on calcule une épaisseur e(x,y) de lentille, en chaque point L(x,y) de la lentille correspondant à un 10 point de contrôle M(x,y), par la formule suivante : e(x,y)= At(x,y)/(1/c-1/c1)+ eo, modulo T/(1/c-1/c1). où : c est la célérité de l'onde ultrasonore dans le milieu hors de la barrière aberratrice, cl est la célérité de l'onde ultrasonore 15 dans la lentille ; eo est un nombre réel tel que e(x,y) ait une valeur positive et supérieure à une épaisseur prédéterminée en tout point de la lentille, - T est la période de l'onde ultrasonore. 20 13.Procédé de fabrication selon l'une quelconque des revendications 8 à 12, dans lequel à la sous-étape (b3), on détermine les premiers temps d'arrivée to(x,y) en simulant une propagation de ladite onde ultrasonore prédéterminée dans un modèle simplifié permettant un calcul 25 analytique. 14.Procédé de fabrication selon l'une quelconque *des revendications précédentes, dans lequel l'étape (c) de calcul de la lentille tient compte des temps de parcours de ladite onde ultrasonore prédéterminée dans la lentille 30 ultrasonore (9) et des angles d'incidences de ladite onde ultrasonore prédéterminée la surface de la lentille ultrasonore (9). 15.Procédé de fabrication selon l'une quelconque des revendications précédentes, dans lequel l'étape (c) de 35 calcul de lentille tient compte de l'ensemble de la 3038217 23 propagation de ladite onde ultrasonore prédéterminée, y compris des réverbérations entre la sonde ultrasonore (2), la lentille ultrasonore (9) et la barrière aberratrice (8). 16.Procédé de fabrication selon l'une quelconque 5 des revendications précédentes, dans lequel l'étape (c) de calcul de lentille se fait itérativement afin d'optimiser le rendement de la transmission d'énergie ultrasonore dans la zone prédéterminée (F), à l'aide d'un modèle analytique ou numérique. 10