Le principe de la radiographie neutronique (nR) repose sur l’atténuation des neutrons à travers la matière qu’ils traversent. Contrairement aux rayons X qui sont diffusés par le nuage d’électrons d’un atome, les neutrons peuvent être absorbés ou dispersés par son noyau. Les neutrons sont sensibles aux éléments légers, tels que l’hydrogène (H), et peuvent donc être utilisés pour radiographier des applications biologiques telles que les tissus animaux 1,2,3,4,5,6,7 ou humains 8,9 et les systèmes souterrains du sol / racines 10,11,12,13,14 ,15. L’imagerie neutronique est une technique complémentaire à l’imagerie par rayons X, qui est capable de détecter des éléments lourds16,17,18. Le nR basé sur l’atténuation est régi par les coefficients d’atténuation linéaires des matériaux dans l’échantillon et par l’épaisseur de l’échantillon, comme décrit par la loi de Beer-Lambert, qui stipule que le faisceau transmis est directement proportionnel à la quantité de matériau et à la longueur du trajet à travers le matériau. Ainsi, la transmittance, T, peut être calculée comme suit:
(1)
où I0 et I sont, respectivement, les intensités incidente et transmise; μ et x sont respectivement le coefficient d’atténuation linéaire et l’épaisseur d’un échantillon homogène. Le coefficient d’atténuation μ est donné par:
(2)
où σ est la section efficace d’atténuation des neutrons de l’échantillon (diffusion et absorption), ρ est sa densité, NA est le nombre d’Avogadro et M est sa masse molaire.
Le contraste en radiographie d’échantillons biologiques utilisant des neutrons de basse énergie (c.-à-d. des énergies inférieures à 0,5 eV) est principalement dû à un changement de densité de H (pour une épaisseur d’échantillon fixe). Cela est dû à la probabilité d’interaction d’un neutron avec le noyau H, qui est plus grande qu’avec les autres noyaux présents dans les échantillons biologiques, et au fait que la densité de l’atome H est primordiale car c’est l’atome le plus abondant dans les échantillons biologiques.
Depuis ses débuts, la nR et la tomodensitométrie neutronique (nCT) ont été largement utilisées pour les applications de matériaux et d’ingénierie 19,20,21,22,23. Les premières expériences de démonstration de la sensibilité des neutrons au H dans des échantillons biologiques ont commencé au milieu des années 195024 avec les mesures de spécimens de plantes. Les travaux se poursuivent tout au long des années 1960 avec, par exemple, la radiographie d’un thorax humain25 ou derats26, dans laquelle l’utilisation d’agents de contraste, tels que l’oxyde de gadolinium (Gd2O3), est explorée. De plus, on a émis l’hypothèse que le contraste entre le tissu tumoral humain et le tissu normal était dû à une augmentation locale de la teneur en H. Au cours de ces essais initiaux, il a été conclu qu’une augmentation du flux de neutrons et de la résolution spatiale améliorerait la qualité du nR et augmenterait probablement sa popularité en tant que technique complémentaire pour des applications industrielles ou biomédicales. Les études les plus récentes comprennent des mesures nR et nCT effectuées sur des échantillons de tissus cancéreux1 et des sections d’organes animaux 2,3,27 pour des applications biomédicales et médico-légales.
Situé au Oak Ridge National Laboratory, Oak Ridge, TN, le High Flux Isotope Reactor (HFIR) est une puissante source de neutrons qui produit des neutrons par réaction de fission. Ces neutrons ont des énergies de l’ordre de 2 MeV et sont « refroidis » dans la piscine du réacteur par des réactions cinétiques avec de l’eau lourde pour atteindre des énergies de l’ordre de 100-300 eV. L’optimisation d’une expérience neutronique, qu’elle soit de diffusion ou d’imagerie, commence par la compréhension de la source de neutrons et des propriétés de la ligne de faisceau telles que l’intensité du faisceau, la distribution de l’énergie et l’effet de fond (neutrons rapides, neutrons retardés, rayons gamma). Dans le hall de guidage froid HFIR où se trouve la ligne de faisceau d’imagerie, les neutrons sont en outre « refroidis » par des interactions cinétiques avec un modérateur H liquide. Ils sont ensuite transportés dans un système de guidage incurvé loin de la ligne de visée de la source, éliminant ainsi les neutrons rapides et la pollution gamma. Comme illustré sur la figure 1, la ligne de faisceau d’imagerie neutroniqueCG-1D 28,29 est placée sur un guide froid, ce qui implique que la gamme d’énergie neutronique varie de quelques meV à quelques dizaines d’eV (dans ce cas, la longueur d’onde de neutrons utilisable correspondante varie de 0,8 à 10 Å) avec un flux compris entre 107 n/(cm2∙s) à la position de l’échantillon. Un système d’ouverture/diffuseur motorisé définit la géométrie du sténopé de l’instrument d’imagerie. Les neutrons parcourent une distance de 6,59 m dans un tube de vol rempli d’hélium (He) avec des fenêtres en aluminium (Al) à chaque extrémité. Les tubes de vol sont utilisés pour transporter les neutrons tout en limitant la diffusion de l’air de sorte que la perte d’intensité du faisceau soit minimale. Pour les mesures décrites dans ce manuscrit, le diffuseur est constitué d’une nano-poudre d’oxyde d’aluminium (Al2O3) de 50 nm d’épaisseur de 1 mm enfermée dans un récipient Al. Le diffuseur réduit les artefacts du faisceau provenant du guide neutronique (qui sont agrandis par la géométrie du sténopé d’une ligne de faisceau d’imagerie), sinon de fortes fluctuations d’intensité horizontales et verticales sont visibles dans la radiographie et la normalisation des données devient difficile. Pour les expériences illustrées ici, les neutrons sont convertis en lumière à l’aide d’un phosphore fluorure/sulfure de zinc de lithium-6 de 25 μm d’épaisseur (6LiF/ZnS:Ag).
L’optimisation de la collimation dépend de la position de l’échantillon au détecteur, de la résolution spatiale requise et du temps d’acquisition. Lorsque l’échantillon se trouve à quelques cm du scintillateur, des collimations élevées (L/D au-dessus de 800, où L est la distance entre l’ouverture du sténopé de diamètre, D et le détecteur) donnent une meilleure résolution spatiale au détriment du flux de neutrons. Une faible collimation (L/D inférieure à 800) est préférable pour les études dynamiques in situ lorsque la résolution temporelle prévaut sur la résolution spatiale. Pour les mesures décrites dans ce manuscrit, la résolution L/D et spatiale était d’environ 355 et 75 μm, respectivement. La résolution temporelle variait en fonction du rapport signal sur bruit (SNR). L’échantillon a été positionné aussi près du scintillateur que possible pour réduire la distorsion géométrique telle que le flou. Des étapes de translation et de rotation sont disponibles pour placer l’échantillon à proximité des détecteurs et effectuer la tomodensitométrie (TDM). CG-1D propose trois types de détecteurs: un dispositif à couplage de charge (CCD) de 2048 pixels x 2048 pixels avec un pas de pixel de 13,5 μm, un détecteur scientifique complémentaire à semi-conducteur à oxyde métallique (sCMOS) de 2560 pixels x 2160 pixels avec un pas de pixel de 6,5 μm et un détecteur à plaque à microcanaux (MCP)30,31 avec 512 pixels x 512 pixels avec une taille de pixel de 55μm. Les neutrons dispersés sont absorbés avec du caoutchouc de bore de ~5 mm d’épaisseur pour protéger la puce du détecteur contre les neutrons visibles. Cette absorption génère des rayons gamma qui peuvent être arrêtés par du plomb (Pb) placé entre le caoutchouc de bore et le détecteur. Chaque détecteur est optimisé pour un champ de vision (FOV) différent ainsi que des résolutions spatiales et temporelles. Pour les mesures du fémur de rat et des poumons de souris, le détecteur CCD a été utilisé pour sa grande capacité de champ de vision (~ 7 cm x 7 cm) et sa résolution spatiale raisonnable d’environ 75 μm. La nCT du système racine/sol de la plante a été réalisée avec la sCMOS, car l’objectif était d’acquérir des nCT le plus rapidement possible au détriment du champ de vision (qui était limité à ~ 5 cm x 4,2 cm); Ainsi, la résolution spatiale en a évidemment souffert. Dans ces détecteurs, les neutrons sont convertis en lumière ou en particule alpha à des fins de détection. La rotation de l’échantillon autour de son axe vertical et l’acquisition de radiographies à des angles de rotation consécutifs permettent l’acquisition de nCT. Le modèle rendu volumétrique en 3 dimensions de l’échantillon étudié est obtenu à l’aide de l’ordinateur portable interne iMARS3D basé sur python Jupyter (FBP), pyMBIR ou d’un logiciel commercial, tous décrits ci-dessous.
Enfin, les neutrons qui n’ont pas interagi avec l’échantillon ou le détecteur sont recueillis en position d’arrêt du faisceau à environ 1 m en aval du système de détection afin de minimiser le bruit de fond. La butée de faisceau CG-1D mesure 0,75 m de large, 0,5 m de haut et 35 mm d’épaisseur et est en B4C en époxy. La butée du faisceau est renforcée avec 10 mm de carbonate de lithium enrichi à 95% (6 Li2CO3) dans un époxy résistant au feu où frappe le faisceau de neutrons, avec une cavité doublée de 6Li, de plomb (Pb) et d’acier conçue pour contenir le taux élevé de rayons gamma secondaires. La butée de poutre est directement fixée au mur de blindage en acier de la ligne de faisceau. Une photographie de la ligne de faisceau CG-1D est donnée à la figure 2.
Trois logiciels de reconstruction ont été utilisés pour reconstruire les trois données expérimentales en 3D, respectivement. La reconstruction de l’échantillon de poumon de souris a été réalisée à l’aide d’Octopus32, un logiciel de reconstruction commercial qui utilise FBP. Le logiciel Octopus se trouve sur un PC serveur et peut être utilisé pour reconstruire les données recueillies à la ligne de faisceau. Un logiciel de reconstruction, nommé iMARS3D, est disponible chez CG-1D. Il est basé sur le code open source TomoPY33 avec des fonctionnalités supplémentaires telles que la correction automatisée de l’inclinaison, les filtres de post-traitement, etc. iMARS3D comprend le prétraitement des données (soustraction de l’arrière-plan et du bruit), le recadrage, le filtrage médian (pour corriger les impacts gamma et les pixels morts), la correction automatisée de la fluctuation de l’intensité du faisceau et la correction de l’inclinaison de l’échantillon. Une fois les sinogrammes créés, d’autres traitements de données tels que la suppression des artefacts annulaires et le lissage sont une option. Les différentes étapes de la reconstruction sont enregistrées sur le serveur d’analyse (puis déplacées dans le dossier partagé de la proposition), tandis que les tranches 2D finales sont immédiatement stockées dans le dossier partagé de la proposition. Le fémur du rat a été reconstruit à l’aide d’iMARS3D. L’échantillon de racine/sol de la plante a été prétraité en filtrant les données médianes à l’aide de TomoPY, suivi d’une correction de l’axe d’inclinaison à l’aide de la bibliothèque SciPy de Python. La reconstruction a été réalisée à l’aide d’un paquet python développé en interne appelé pyMBIR (construit à l’aide de noyaux de la boîte à outils ASTRA34) qui implémente une suite d’algorithmes tomographiques du FBP de base aux techniques avancées de reconstruction itérative basées sur des modèles35 qui peuvent obtenir des reconstructions de haute qualité à partir d’ensembles de données sur les neutrons extrêmement clairsemés et bruyants. Tous les volumes rendus basés sur les outils de reconstruction susmentionnés sont représentés en contraste d’atténuation. Toute la visualisation a été réalisée à l’aide du progiciel commercial de visualisation, de segmentation et d’analyse de données AMIRA36.
Ce manuscrit vise à démontrer la procédure d’utilisation de l’imagerie neutronique (nR et nCT) sur la ligne de faisceau HFIR CG-1D. Cette étude illustre également les capacités actuelles de pointe nR et nCT pour les échantillons biologiques, en particulier un poumon de souris, un os de rat et des systèmes racinaires / sol végétaux. Le poumon de souris a été choisi pour illustrer la complémentarité des neutrons pour mesurer le tissu pulmonaire, alors que les rayons X sont surtout sensibles aux os. L’échantillon osseux, un fémur de rat, avait un implant en titane (Ti), illustrant ainsi le contraste entre l’os et le métal, et la possibilité de voir l’interface os/métal (ce qui est difficile à mesurer avec les rayons X car les métaux les atténuent fortement4). Enfin, le système d’eau plante-racine illustre la capacité tridimensionnelle (3D) de la nCT à mesurer les systèmes racinaires/sol in situ. Il montre en outre les avantages / inconvénients de l’utilisation de nR pour les échantillons biologiques. De toute évidence, cette méthode peut être utilisée en toute sécurité pour mesurer la dynamique de l’eau dans un système racinaire végétal, mais ne peut pas être considérée comme une technique d’imagerie animale ou humaine vivante en raison des risques associés à l’exposition aux rayonnements, limitant ainsi les études à des souris (mortes) ou à des mesures de type pathologie dans lesquelles, par exemple, un échantillon de tissu est réséqué d’un patient (animal ou humain) et préparé par fixation avant d’être mesuré dans un faisceau de neutrons.