Calibrated simulations for dynamic focusing of ultrasound through the temporal window

Résumé technique : Simulations calibrées pour la focalisation dynamique des ultrasons à travers la fenêtre temporale

Énoncé du problème
La neuromodulation par ultrasons focalisés (FUS) émerge comme un outil pour les applications cliniques thérapeutiques et l'investigation causale du cerveau. Un défi critique pour le déploiement de réseaux de transducteurs multi-éléments destinés au pilotage dynamique réside dans l'écart entre les simulations acoustiques et les pressions focales mesurées. Les transducteurs à foyer unique, qui reposent souvent sur une valeur de calibration unique, sont simples, mais les réseaux multi-éléments capables de diriger le faisceau dans un espace 3D introduisent des erreurs complexes dues au diaphonie (cross-talk) mécanique ou électrique et à des hypothions inexactes concernant la transmission par le crâne. Les modèles de simulation existants ne parviennent souvent pas à rendre compte de ces facteurs, entraînant des décalages importants (dépassant 20 %) entre les pressions de sortie prédites et réelles. De plus, l'hétérogénéité microstructurale du crâne, particulièrement la diffusion de l'os trabéculaire, est difficile à simuler avec précision avec les résolutions d'imagerie actuelles, ce qui conduit à une sous-estimation des pertes d'atténuation. L'article traite de la nécessité de méthodes de calibration efficaces et précises pour combler l'écart entre la simulation et la réalité afin de garantir la sécurité et la précision des traitements de neuromodulation personnalisés.

Méthodologie
Les auteurs ont développé un cadre de calibration combinant un échantillonnage expérimental parcimonieux et des simulations acoustiques optimisées à travers un sous-espace pilotable en 3D.

• Définition de l'espace cible : Les coordonnées relatives des cibles pour quatre structures cérébrales profondes bilatérales (noyau caudé, amygdale, noyau centromédian du thalamus et noyau subthalamique) ont été extraites de 157 sujets. Ces coordonnées ont défini un espace de calibration en forme de cône tronqué capturant 95 % des cibles potentielles.
• Calibration Simulation-Mesure :
  • Holographie acoustique : Les tentatives initiales d'utilisation de l'holographie pour la cartographie de la pression de la source ont révélé des décalages significatifs (> 20 %) dus à la diaphonie des éléments non modélisée, nécessitant une mesure directe à travers l'espace de pilotage.
  • Échantillonnage parcimonieux : Les auteurs ont déterminé qu'un plan de calibration 2D pouvait être caractérisé avec précision par l'échantillonnage de 50 points avec un espacement équidistant de 6 mm, en ajustant les données de décalage avec une fonction gaussienne anisotrope ($R^2 = 0,83$). Pour l'espace 3D, la sélection de modèles bayésiens et la détection du coude (elbow detection) ont indiqué que 20 points étaient suffisants pour modéliser le décalage, une conclusion qui s'est généralisée à cinq transducteurs ayant des configurations identiques.
  • Axi-symétrie : Bien que les simulations suggèrent des champs axisymétriques, les mesures expérimentales de 8 plans à des pas de 45° ont révélé de légères différences dues à la variabilité de fabrication, suggérant qu'une caractérisation 3D complète est nécessaire pour les réseaux non axisymétriques.
• Optimisation de l'atténuation du crâne :
  • Validation ex vivo : Des fragments de fenêtres temporales ont été extraits de 5 crânes de cadavres (produisant 8 fragments utilisables) et scannés par IRM Zero Echo Time (ZTE). Ces fragments ont été validés par rapport à des extractions virtuelles de 157 sujets vivants via IRM, confirmant la similitude de l'épaisseur et de l'angle d'incidence.
  • Dérivation du coefficient d'atténuation : En utilisant des simulations de renversement temporel (time-reversal) ciblant 40 emplacements dans l'espace de calibration pour chaque fragment, les auteurs ont modifié de manière itérative le coefficient d'atténuation de l'os homogène ($\alpha_0$). La valeur minimisant l'erreur quadratique moyenne entre les pressions de crête simulées et mesurées a été déterminée à 15,3 dB/cm/MHz.
  • Limites d'erreur : L'étude a calculé les erreurs de pression aux 5e et 95e percentiles, trouvant que les pressions réelles pouvaient être 44 % inférieures ou 33 % supérieures aux valeurs simulées en utilisant le coefficient optimisé.

Résultats clés

• Efficacité de l'échantillonnage : L'étude a établi qu'un échantillonnage de champ à haute résolution n'est pas nécessaire pour la calibration ; 50 points pour les plans 2D et 20 points pour les volumes 3D suffisent pour modéliser le décalage entre simulation et réalité à l'aide de fonctions gaussiennes anisotropes.
• Atténuation optimisée : Un coefficient d'atténuation osseuse optimisé de 15,3 dB/cm/MHz a été dérivé pour la fenêtre temporale, s'alignant sur les estimations expérimentales précédentes (13,3–16,6 dB/cm/MHz à 1 MHz).
• Caractérisation du décalage : Le décalage entre la simulation et la mesure contient une composante graduelle à basse fréquence, les plus grands écarts se produisant aux points cibles latéraux en champ proche.
• Généralisabilité : L'approche de calibration et le coefficient d'atténuation optimisé se sont révélés applicables à différents designs de réseaux (spirale de Fermat, cercle concentrique, grille) et à plusieurs unités de transducteurs.

Signification et affirmations
L'article affirme que ce cadre de calibration et les coefficients de la fenêtre temporale optimisés offrent une solution pratique pour améliorer la précision des rapports et des évaluations de sécurité pour la neuromodulation personnalisée. En s'élochant des estimations uniques et conservatrices ou des simulations non calibrées, les auteurs proposent une approche basée sur les données qui :

1. Améliore la sécurité : Capturer la pression de crête réelle plutôt que la position focale prévue garantit des calibrations de sécurité conservatrices pour des paramètres tels que l'indice mécanique (MI) et l'intensité spatiale de crête moyenne pulsée ($I_{SPPA}$).
2. Améliore la précision : Le coefficient d'atténuation optimisé tient compte des pertes de diffusion inhérentes à la microstructure osseuse qui sont typiquement sous-estimées dans les simulations.
3. Facilite l'extensibilité : La méthode d'échantillonnage parcimonieux (20–50 points) rend la calibration réalisable pour des dispositages individuels déployés sur le terrain, répondant aux défis logistiques de la caractérisation de chaque point dans une plage de pilotage.

Les auteurs soulignent que bien que les dérivés thermiques soient également critiques pour la sécurité, les champs de pression calibrés dérivés ici peuvent être directement intégrés dans les simulations thermiques (par exemple, l'équation de biochauffe de Pennes) en ajustant le coefficient d'absorption pour représenter la fraction non diffusante de l'atténuation. L'étude conclut qu'une spécification claire de ces méthodes d'atténuation est essentielle pour guider les efforts de réplication et assurer la sécurité des réseaux d'ultrasons pilotables dans les applications humaines.